Comprehensive patient-level classification and quantification of driver events in TCGA PanCanAtlas cohorts

There is a growing need to develop novel therapeutics for targeted treatment of cancer. The prerequisite to success is the knowledge about which types of molecular alterations are predominantly driving tumorigenesis. To shed light onto this subject, we have utilized the largest database of human cancer mutations–TCGA PanCanAtlas, multiple established algorithms for cancer driver prediction (2020plus, CHASMplus, CompositeDriver, dNdScv, DriverNet, HotMAPS, OncodriveCLUSTL, OncodriveFML) and developed four novel computational pipelines: SNADRIF (Single Nucleotide Alteration DRIver Finder), GECNAV (Gene Expression-based Copy Number Alteration Validator), ANDRIF (ANeuploidy DRIver Finder) and PALDRIC (PAtient-Level DRIver Classifier). A unified workflow integrating all these pipelines, algorithms and datasets at cohort and patient levels was created. We have found that there are on average 12 driver events per tumour, of which 0.6 are single nucleotide alterations (SNAs) in oncogenes, 1.5 are amplifications of oncogenes, 1.2 are SNAs in tumour suppressors, 2.1 are deletions of tumour suppressors, 1.5 are driver chromosome losses, 1 is a driver chromosome gain, 2 are driver chromosome arm losses, and 1.5 are driver chromosome arm gains. The average number of driver events per tumour increases with age (from 7 to 15) and cancer stage (from 10 to 15) and varies strongly between cancer types (from 1 to 24). Patients with 1 and 7 driver events per tumour are the most frequent, and there are very few patients with more than 40 events. In tumours having only one driver event, this event is most often an SNA in an oncogene. However, with increasing number of driver events per tumour, the contribution of SNAs decreases, whereas the contribution of copy-number alterations and aneuploidy events increases.

Shp2 in uterine stromal cells critically regulates on time embryo implantation and stromal decidualization by multiple pathways during early pregnancy

Approximately 75% of failed pregnancies are considered to be due to embryo implantation failure or defects. Nevertheless, the explicit signaling mechanisms governing this process have not yet been elucidated. Here, we found that conditional deletion of the Shp2 gene in mouse uterine stromal cells deferred embryo implantation and inhibited the decidualization of stromal cells, which led to embryonic developmental delay and to the death of numerous embryos mid-gestation, ultimately reducing female fertility. The absence of Shp2 in stromal cells increased the proliferation of endometrial epithelial cells, thereby disturbing endometrial epithelial remodeling. However, Shp2 deletion impaired the proliferation and polyploidization of stromal cells, which are distinct characteristics of decidualization. In human endometrial stromal cells (hESCs), Shp2 expression gradually increased during the decidualization process. Knockout of Shp2 blocked the decidual differentiation of hESCs, while Shp2 overexpression had the opposite effect. Shp2 knockout inhibited the proliferation of hESCs during decidualization. Whole gene expression profiling analysis of hESCs during the decidualization process showed that Shp2 deficiency disrupted many signaling transduction pathways and gene expression. Analyses of hESCs and mouse uterine tissues confirmed that the signaling pathways extracellular regulated protein kinases (ERK), protein kinase B (AKT), signal transducer and activator of transcription 3 (STAT3) and their downstream transcription factors CCAAT/enhancer binding protein β (C/EBPβ) and Forkhead box transcription factor O1 (FOXO-1) were involved in the Shp2 regulation of decidualization. In summary, these results demonstrate that Shp2 plays a crucial role in stromal decidualization by mediating and coordinating multiple signaling pathways in uterine stromal cells. Our discovery possibly provides a novel key regulator of embryo implantation and novel therapeutic target for pregnancy failure.

Activation of the ubiquitin-proteasome system contributes to oculopharyngeal muscular dystrophy through muscle atrophy

Oculopharyngeal muscular dystrophy (OPMD) is a late-onset disorder characterized by progressive weakness and degeneration of specific muscles. OPMD is due to extension of a polyalanine tract in poly(A) binding protein nuclear 1 (PABPN1). Aggregation of the mutant protein in muscle nuclei is a hallmark of the disease. Previous transcriptomic analyses revealed the consistent deregulation of the ubiquitin-proteasome system (UPS) in OPMD animal models and patients, suggesting a role of this deregulation in OPMD pathogenesis. Subsequent studies proposed that UPS contribution to OPMD involved PABPN1 aggregation. Here, we use a Drosophila model of OPMD to address the functional importance of UPS deregulation in OPMD. Through genome-wide and targeted genetic screens we identify a large number of UPS components that are involved in OPMD. Half dosage of UPS genes reduces OPMD muscle defects suggesting a pathological increase of UPS activity in the disease. Quantification of proteasome activity confirms stronger activity in OPMD muscles, associated with degradation of myofibrillar proteins. Importantly, improvement of muscle structure and function in the presence of UPS mutants does not correlate with the levels of PABPN1 aggregation, but is linked to decreased degradation of muscle proteins. Oral treatment with the proteasome inhibitor MG132 is beneficial to the OPMD Drosophila model, improving muscle function although PABPN1 aggregation is enhanced. This functional study reveals the importance of increased UPS activity that underlies muscle atrophy in OPMD. It also provides a proof-of-concept that inhibitors of proteasome activity might be an attractive pharmacological approach for OPMD.

Paternally expressed gene 3 (Pw1/Peg3) promotes sexual dimorphism in metabolism and behavior

The paternally expressed gene 3 (Pw1/Peg3) is a mammalian-specific parentally imprinted gene expressed in stem/progenitor cells of the brain and endocrine tissues. Here, we compared phenotypic characteristics in Pw1/Peg3 deficient male and female mice. Our findings indicate that Pw1/Peg3 is a key player for the determination of sexual dimorphism in metabolism and behavior. Mice carrying a paternally inherited Pw1/Peg3 mutant allele manifested postnatal deficits in GH/IGF dependent growth before weaning, sex steroid dependent masculinization during puberty, and insulin dependent fat accumulation in adulthood. As a result, Pw1/Peg3 deficient mice develop a sex-dependent global shift of body metabolism towards accelerated adiposity, diabetic-like insulin resistance, and fatty liver. Furthermore, Pw1/Peg3 deficient males displayed reduced social dominance and competitiveness concomitant with alterations in the vasopressinergic architecture in the brain. This study demonstrates that Pw1/Peg3 provides an epigenetic context that promotes male-specific characteristics through sex steroid pathways during postnatal development.

Repurposing of the enhancer-promoter communication underlies the compensation of Mesp2 by Mesp1

Organisms are inherently equipped with buffering systems against genetic perturbations. Genetic compensation, the compensatory response by upregulating another gene or genes, is one such buffering mechanism. Recently, a well-conserved compensatory mechanism was proposed: transcriptional adaptation of homologs under the nonsense-mediated mRNA decay pathways. However, this model cannot explain the onset of all compensatory events. We report a novel genetic compensation mechanism operating over the Mesp gene locus. Mesp1 and Mesp2 are paralogs located adjacently in the genome. Mesp2 loss is partially rescued by Mesp1 upregulation in the presomitic mesoderm (PSM). Using a cultured PSM induction system, we reproduced the compensatory response in vitro and found that the Mesp2-enhancer is required to promote Mesp1. We revealed that the Mesp2-enhancer directly interacts with the Mesp1 promoter, thereby upregulating Mesp1 expression upon the loss of Mesp2. Of note, this interaction is established by genomic arrangement upon PSM development independently of Mesp2 disruption. We propose that the repurposing of this established enhancer-promoter communication is the mechanism underlying this compensatory response for the upregulation of the adjacent gene.

WAPO-A1 is the causal gene of the 7AL QTL for spikelet number per spike in wheat

Improving our understanding of the genes regulating grain yield can contribute to the development of more productive wheat varieties. Previously, a highly significant QTL affecting spikelet number per spike (SNS), grain number per spike (GNS) and grain yield was detected on chromosome arm 7AL in multiple genome-wide association studies. Using a high-resolution genetic map, we established that the A-genome homeolog of WHEAT ORTHOLOG OF APO1 (WAPO-A1) was a leading candidate gene for this QTL. Using mutants and transgenic plants, we demonstrate in this study that WAPO-A1 is the causal gene underpinning this QTL. Loss-of-function mutants wapo-A1 and wapo-B1 showed reduced SNS in tetraploid wheat, and the effect was exacerbated in wapo1 combining both mutations. By contrast, spikes of transgenic wheat plants carrying extra copies of WAPO-A1 driven by its native promoter had higher SNS, a more compact spike apical region and a smaller terminal spikelet than the wild type. Taken together, these results indicate that WAPO1 affects SNS by regulating the timing of terminal spikelet formation. Both transgenic and wapo1 mutant plants showed a wide range of floral abnormalities, indicating additional roles of WAPO1 on wheat floral development. Previously, we found three widespread haplotypes in the QTL region (H1, H2 and H3), each associated with particular WAPO-A1 alleles. Results from this and our previous study, show that the WAPO-A1 allele in the H1 haplotype (115-bp deletion in the promoter) is expressed at significantly lower levels in the developing spikes than the alleles in the H2 and H3 haplotypes, resulting in reduced SNS. Field experiments also showed that the H2 haplotype is associated with the strongest effects in increasing SNS and GNS (H2>H3>H1). The H2 haplotype is already present in most modern common wheat varieties but is rare in durum wheat, where it might be particularly useful to improve grain yield.

Analysis of HubP-dependent cell pole protein targeting in Vibrio cholerae uncovers novel motility regulators

In rod-shaped bacteria, the emergence and maintenance of long-axis cell polarity is involved in key cellular processes such as cell cycle, division, environmental sensing and flagellar motility among others. Many bacteria achieve cell pole differentiation through the use of polar landmark proteins acting as scaffolds for the recruitment of functional macromolecular assemblies. In Vibrio cholerae a large membrane-tethered protein, HubP, specifically interacts with proteins involved in chromosome segregation, chemotaxis and flagellar biosynthesis. Here we used comparative proteomics, genetic and imaging approaches to identify additional HubP partners and demonstrate that at least six more proteins are subject to HubP-dependent polar localization. These include a cell-wall remodeling enzyme (DacB), a likely chemotaxis sensory protein (HlyB), two presumably cytosolic proteins of unknown function (VC1210 and VC1380) and two membrane-bound proteins, named here MotV and MotW, that exhibit distinct effects on chemotactic motility. We show that while both ΔmotW and ΔmotV mutants retain monotrichous flagellation, they present significant to severe motility defects when grown in soft agar. Video-tracking experiments further reveal that ΔmotV cells can swim in liquid environments but are unable to tumble or penetrate a semisolid matrix, whereas a motW deletion affects both tumbling frequency and swimming speed. Motility suppressors and gene co-occurrence analyses reveal co-evolutionary linkages between MotV, a subset of non-canonical CheV proteins and flagellar C-ring components FliG and FliM, whereas MotW regulatory inputs appear to intersect with specific c-di-GMP signaling pathways. Together, these results reveal an ever more versatile role for the landmark cell pole organizer HubP and identify novel mechanisms of motility regulation.

Bacterial recognition by PGRP-SA and downstream signalling by Toll/DIF sustain commensal gut bacteria in Drosophila

The gut sets the immune and metabolic parameters for the survival of commensal bacteria. We report that in Drosophila, deficiency in bacterial recognition upstream of Toll/NF-κB signalling resulted in reduced density and diversity of gut bacteria. Translational regulation factor 4E-BP, a transcriptional target of Toll/NF-κB, mediated this host-bacteriome interaction. In healthy flies, Toll activated 4E-BP, which enabled fat catabolism, which resulted in sustaining of the bacteriome. The presence of gut bacteria kept Toll signalling activity thus ensuring the feedback loop of their own preservation. When Toll activity was absent, TOR-mediated suppression of 4E-BP made fat resources inaccessible and this correlated with loss of intestinal bacterial density. This could be overcome by genetic or pharmacological inhibition of TOR, which restored bacterial density. Our results give insights into how an animal integrates immune sensing and metabolism to maintain indigenous bacteria in a healthy gut.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

Chromatin profiling reveals heterogeneity in clinical isolates of the human pathogen Aspergillus fumigatus

Invasive Pulmonary Aspergillosis, which is caused by the filamentous fungus Aspergillus fumigatus, is a life-threatening infection for immunosuppressed patients. Chromatin structure regulation is important for genome stability maintenance and has the potential to drive genome rearrangements and affect virulence and pathogenesis of pathogens. Here, we performed the first A. fumigatus global chromatin profiling of two histone modifications, H3K4me3 and H3K9me3, focusing on the two most investigated A. fumigatus clinical isolates, Af293 and CEA17. In eukaryotes, H3K4me3 is associated with active transcription, while H3K9me3 often marks silent genes, DNA repeats, and transposons. We found that H3K4me3 deposition is similar between the two isolates, while H3K9me3 is more variable and does not always represent transcriptional silencing. Our work uncovered striking differences in the number, locations, and expression of transposable elements between Af293 and CEA17, and the differences are correlated with H3K9me3 modifications and higher genomic variations among strains of Af293 background. Moreover, we further showed that the Af293 strains from different laboratories actually differ in their genome contents and found a frequently lost region in chromosome VIII. For one such Af293 variant, we identified the chromosomal changes and demonstrated their impacts on its secondary metabolites production, growth and virulence. Overall, our findings not only emphasize the influence of genome heterogeneity on A. fumigatus fitness, but also caution about unnoticed chromosomal variations among common laboratory strains.