Genetic and epigenetic features of bilateral Wilms tumor predisposition in patients from the Children's Oncology Group AREN18B5-Q.

Developing synchronous bilateral Wilms tumor suggests an underlying (epi)genetic predisposition. Here, we evaluate this predisposition in 68 patients using whole exome or genome sequencing (n = 85 tumors from 61 patients with matched germline blood DNA), RNA-seq (n = 99 tumors), and DNA methylation analysis (n = 61 peripheral blood, n = 29 non-diseased kidney, n = 99 tumors). We determine the predominant events for bilateral Wilms tumor predisposition: 1)pre-zygotic germline genetic variants readily detectable in blood DNA [WT1 (14.8%), NYNRIN (6.6%), TRIM28 (5%), and BRCA-related genes (5%)] or 2)post-zygotic epigenetic hypermethylation at 11p15.5 H19/ICR1 that may require analysis of multiple tissue types for diagnosis. Of 99 total tumor specimens, 16 (16.1%) have 11p15.5 normal retention of imprinting, 25 (25.2%) have 11p15.5 copy neutral loss of heterozygosity, and 58 (58.6%) have 11p15.5 H19/ICR1 epigenetic hypermethylation (loss of imprinting). Here, we ascertain the epigenetic and genetic modes of bilateral Wilms tumor predisposition.

Genome-wide mapping of cancer dependency genes and genetic modifiers of chemotherapy in high-risk hepatoblastoma.

A lack of relevant genetic models and cell lines hampers our understanding of hepatoblastoma pathogenesis and the development of new therapies for this neoplasm. Here, we report an improved MYC-driven hepatoblastoma-like murine model that recapitulates the pathological features of embryonal type of hepatoblastoma, with transcriptomics resembling the high-risk gene signatures of the human disease. Single-cell RNA-sequencing and spatial transcriptomics identify distinct subpopulations of hepatoblastoma cells. After deriving cell lines from the mouse model, we map cancer dependency genes using CRISPR-Cas9 screening and identify druggable targets shared with human hepatoblastoma (e.g., CDK7, CDK9, PRMT1, PRMT5). Our screen also reveals oncogenes and tumor suppressor genes in hepatoblastoma that engage multiple, druggable cancer signaling pathways. Chemotherapy is critical for human hepatoblastoma treatment. A genetic mapping of doxorubicin response by CRISPR-Cas9 screening identifies modifiers whose loss-of-function synergizes with (e.g., PRKDC) or antagonizes (e.g., apoptosis genes) the effect of chemotherapy. The combination of PRKDC inhibition and doxorubicin-based chemotherapy greatly enhances therapeutic efficacy. These studies provide a set of resources including disease models suitable for identifying and validating potential therapeutic targets in human high-risk hepatoblastoma.

The Genetic and Epigenetic Features of Bilateral Wilms Tumor Predisposition: A Report from the Children's Oncology Group AREN18B5-Q Study.

This study comprehensively evaluated the landscape of genetic and epigenetic events that predispose to synchronous bilateral Wilms tumor (BWT). We performed whole exome or whole genome sequencing, total-strand RNA-seq, and DNA methylation analysis using germline and/or tumor samples from 68 patients with BWT from St. Jude Children's Research Hospital and the Children's Oncology Group. We found that 25/61 (41%) of patients evaluated harbored pathogenic or likely pathogenic germline variants, with WT1 (14.8%), NYNRIN (6.6%), TRIM28 (5%) and the BRCA-related genes (5%) BRCA1, BRCA2, and PALB2 being most common. Germline WT1 variants were strongly associated with somatic paternal uniparental disomy encompassing the 11p15.5 and 11p13/WT1 loci and subsequent acquired pathogenic CTNNB1 variants. Somatic coding variants or genome-wide copy number alterations were almost never shared between paired synchronous BWT, suggesting that the acquisition of independent somatic variants leads to tumor formation in the context of germline or early embryonic, post-zygotic initiating events. In contrast, 11p15.5 status (loss of heterozygosity, loss or retention of imprinting) was shared among paired synchronous BWT in all but one case. The predominant molecular events for BWT predisposition include pathogenic germline variants or post-zygotic epigenetic hypermethylation at the 11p15.5 H19/ICR1 locus (loss of imprinting). This study demonstrates that post-zygotic somatic mosaicism for 11p15.5 hypermethylation/loss of imprinting is the single most common initiating molecular event predisposing to BWT. Evidence of somatic mosaicism for 11p15.5 loss of imprinting was detected in leukocytes of a cohort of BWT patients and long-term survivors, but not in unilateral Wilms tumor patients and long-term survivors or controls, further supporting the hypothesis that post-zygotic 11p15.5 alterations occurred in the mesoderm of patients who go on to develop BWT. Due to the preponderance of BWT patients with demonstrable germline or early embryonic tumor predisposition, BWT exhibits a unique biology when compared to unilateral Wilms tumor and therefore warrants continued refinement of its own treatment-relevant biomarkers which in turn may inform directed treatment strategies in the future.

Identification and targeting of a HES1-YAP1-CDKN1C functional interaction in fusion-negative rhabdomyosarcoma.

Rhabdomyosarcoma (RMS), a cancer characterized by features of skeletal muscle, is the most common soft-tissue sarcoma of childhood. With 5-year survival rates among high-risk groups at < 30%, new therapeutics are desperately needed. Previously, using a myoblast-based model of fusion-negative RMS (FN-RMS), we found that expression of the Hippo pathway effector transcriptional coactivator YAP1 (YAP1) permitted senescence bypass and subsequent transformation to malignant cells, mimicking FN-RMS. We also found that YAP1 engages in a positive feedback loop with Notch signaling to promote FN-RMS tumorigenesis. However, we could not identify an immediate downstream impact of this Hippo-Notch relationship. Here, we identify a HES1-YAP1-CDKN1C functional interaction, and show that knockdown of the Notch effector HES1 (Hes family BHLH transcription factor 1) impairs growth of multiple FN-RMS cell lines, with knockdown resulting in decreased YAP1 and increased CDKN1C expression. In silico mining of published proteomic and transcriptomic profiles of human RMS patient-derived xenografts revealed the same pattern of HES1-YAP1-CDKN1C expression. Treatment of FN-RMS cells in vitro with the recently described HES1 small-molecule inhibitor, JI130, limited FN-RMS cell growth. Inhibition of HES1 in vivo via conditional expression of a HES1-directed shRNA or JI130 dosing impaired FN-RMS tumor xenograft growth. Lastly, targeted transcriptomic profiling of FN-RMS xenografts in the context of HES1 suppression identified associations between HES1 and RAS-MAPK signaling. In summary, these in vitro and in vivo preclinical studies support the further investigation of HES1 as a therapeutic target in FN-RMS.

Molecular evolution and the decline of purifying selection with age.

Life history theory predicts that the intensity of selection declines with age, and this trend should impact how genes expressed at different ages evolve. Here we find consistent relationships between a gene's age of expression and patterns of molecular evolution in two mammals (the human Homo sapiens and the mouse Mus musculus) and two insects (the malaria mosquito Anopheles gambiae and the fruit fly Drosophila melanogaster). When expressed later in life, genes fix nonsynonymous mutations more frequently, are more polymorphic for nonsynonymous mutations, and have shorter evolutionary lifespans, relative to those expressed early. The latter pattern is explained by a simple evolutionary model. Further, early-expressed genes tend to be enriched in similar gene ontology terms across species, while late-expressed genes show no such consistency. In humans, late-expressed genes are more likely to be linked to cancer and to segregate for dominant disease-causing mutations. Last, the effective strength of selection (Ne s) decreases and the fraction of beneficial mutations increases with a gene's age of expression. These results are consistent with the diminishing efficacy of purifying selection with age, as proposed by Medawar's classic hypothesis for the evolution of senescence, and provide links between life history theory and molecular evolution.

Latent cellular analysis robustly reveals subtle diversity in large-scale single-cell RNA-seq data.

Single-cell RNA sequencing (scRNA-seq) is a powerful tool for characterizing the cell-to-cell variation and cellular dynamics in populations which appear homogeneous otherwise in basic and translational biological research. However, significant challenges arise in the analysis of scRNA-seq data, including the low signal-to-noise ratio with high data sparsity, potential batch effects, scalability problems when hundreds of thousands of cells are to be analyzed among others. The inherent complexities of scRNA-seq data and dynamic nature of cellular processes lead to suboptimal performance of many currently available algorithms, even for basic tasks such as identifying biologically meaningful heterogeneous subpopulations. In this study, we developed the Latent Cellular Analysis (LCA), a machine learning-based analytical pipeline that combines cosine-similarity measurement by latent cellular states with a graph-based clustering algorithm. LCA provides heuristic solutions for population number inference, dimension reduction, feature selection, and control of technical variations without explicit gene filtering. We show that LCA is robust, accurate, and powerful by comparison with multiple state-of-the-art computational methods when applied to large-scale real and simulated scRNA-seq data. Importantly, the ability of LCA to learn from representative subsets of the data provides scalability, thereby addressing a significant challenge posed by growing sample sizes in scRNA-seq data analysis.

Metabolic signaling directs the reciprocal lineage decisions of αß and γδ T cells.

The interaction between extrinsic factors and intrinsic signal strength governs thymocyte development, but the mechanisms linking them remain elusive. We report that mechanistic target of rapamycin complex 1 (mTORC1) couples microenvironmental cues with metabolic programs to orchestrate the reciprocal development of two fundamentally distinct T cell lineages, the αß and γδ T cells. Developing thymocytes dynamically engage metabolic programs including glycolysis and oxidative phosphorylation, as well as mTORC1 signaling. Loss of RAPTOR-mediated mTORC1 activity impairs the development of αß T cells but promotes γδ T cell generation, associated with disrupted metabolic remodeling of oxidative and glycolytic metabolism. Mechanistically, we identify mTORC1-dependent control of reactive oxygen species production as a key metabolic signal in mediating αß and γδ T cell development, and perturbation of redox homeostasis impinges upon thymocyte fate decisions and mTORC1-associated phenotypes. Furthermore, single-cell RNA sequencing and genetic dissection reveal that mTORC1 links developmental signals from T cell receptors and NOTCH to coordinate metabolic activity and signal strength. Our results establish mTORC1-driven metabolic signaling as a decisive factor for reciprocal αß and γδ T cell development and provide insight into metabolic control of cell signaling and fate decisions.

Radical remodeling of the Y chromosome in a recent radiation of malaria mosquitoes.

Y chromosomes control essential male functions in many species, including sex determination and fertility. However, because of obstacles posed by repeat-rich heterochromatin, knowledge of Y chromosome sequences is limited to a handful of model organisms, constraining our understanding of Y biology across the tree of life. Here, we leverage long single-molecule sequencing to determine the content and structure of the nonrecombining Y chromosome of the primary African malaria mosquito, Anopheles gambiae We find that the An. gambiae Y consists almost entirely of a few massively amplified, tandemly arrayed repeats, some of which can recombine with similar repeats on the X chromosome. Sex-specific genome resequencing in a recent species radiation, the An. gambiae complex, revealed rapid sequence turnover within An. gambiae and among species. Exploiting 52 sex-specific An. gambiae RNA-Seq datasets representing all developmental stages, we identified a small repertoire of Y-linked genes that lack X gametologs and are not Y-linked in any other species except An. gambiae, with the notable exception of YG2, a candidate male-determining gene. YG2 is the only gene conserved and exclusive to the Y in all species examined, yet sequence similarity to YG2 is not detectable in the genome of a more distant mosquito relative, suggesting rapid evolution of Y chromosome genes in this highly dynamic genus of malaria vectors. The extensive characterization of the An. gambiae Y provides a long-awaited foundation for studying male mosquito biology, and will inform novel mosquito control strategies based on the manipulation of Y chromosomes.