Quantifying intra-tumoral genetic heterogeneity of glioblastoma toward precision medicine using MRI and a data-inclusive machine learning algorithm.

BACKGROUND AND OBJECTIVE: Glioblastoma (GBM) is one of the most aggressive and lethal human cancers. Intra-tumoral genetic heterogeneity poses a significant challenge for treatment. Biopsy is invasive, which motivates the development of non-invasive, MRI-based machine learning (ML) models to quantify intra-tumoral genetic heterogeneity for each patient. This capability holds great promise for enabling better therapeutic selection to improve patient outcome. METHODS: We proposed a novel Weakly Supervised Ordinal Support Vector Machine (WSO-SVM) to predict regional genetic alteration status within each GBM tumor using MRI. WSO-SVM was applied to a unique dataset of 318 image-localized biopsies with spatially matched multiparametric MRI from 74 GBM patients. The model was trained to predict the regional genetic alteration of three GBM driver genes (EGFR, PDGFRA and PTEN) based on features extracted from the corresponding region of five MRI contrast images. For comparison, a variety of existing ML algorithms were also applied. Classification accuracy of each gene were compared between the different algorithms. The SHapley Additive exPlanations (SHAP) method was further applied to compute contribution scores of different contrast images. Finally, the trained WSO-SVM was used to generate prediction maps within the tumoral area of each patient to help visualize the intra-tumoral genetic heterogeneity. RESULTS: WSO-SVM achieved 0.80 accuracy, 0.79 sensitivity, and 0.81 specificity for classifying EGFR; 0.71 accuracy, 0.70 sensitivity, and 0.72 specificity for classifying PDGFRA; 0.80 accuracy, 0.78 sensitivity, and 0.83 specificity for classifying PTEN; these results significantly outperformed the existing ML algorithms. Using SHAP, we found that the relative contributions of the five contrast images differ between genes, which are consistent with findings in the literature. The prediction maps revealed extensive intra-tumoral region-to-region heterogeneity within each individual tumor in terms of the alteration status of the three genes. CONCLUSIONS: This study demonstrated the feasibility of using MRI and WSO-SVM to enable non-invasive prediction of intra-tumoral regional genetic alteration for each GBM patient, which can inform future adaptive therapies for individualized oncology.

Integrated proteogenomic characterization of glioblastoma evolution.

The evolutionary trajectory of glioblastoma (GBM) is a multifaceted biological process that extends beyond genetic alterations alone. Here, we perform an integrative proteogenomic analysis of 123 longitudinal glioblastoma pairs and identify a highly proliferative cellular state at diagnosis and replacement by activation of neuronal transition and synaptogenic pathways in recurrent tumors. Proteomic and phosphoproteomic analyses reveal that the molecular transition to neuronal state at recurrence is marked by post-translational activation of the wingless-related integration site (WNT)/ planar cell polarity (PCP) signaling pathway and BRAF protein kinase. Consistently, multi-omic analysis of patient-derived xenograft (PDX) models mirror similar patterns of evolutionary trajectory. Inhibition of B-raf proto-oncogene (BRAF) kinase impairs both neuronal transition and migration capability of recurrent tumor cells, phenotypic hallmarks of post-therapy progression. Combinatorial treatment of temozolomide (TMZ) with BRAF inhibitor, vemurafenib, significantly extends the survival of PDX models. This study provides comprehensive insights into the biological mechanisms of glioblastoma evolution and treatment resistance, highlighting promising therapeutic strategies for clinical intervention.

The Epigenetic Evolution of Glioma Is Determined by the IDH1 Mutation Status and Treatment Regimen.

Tumor adaptation or selection is thought to underlie therapy resistance in glioma. To investigate longitudinal epigenetic evolution of gliomas in response to therapeutic pressure, we performed an epigenomic analysis of 132 matched initial and recurrent tumors from patients with IDH-wildtype (IDHwt) and IDH-mutant (IDHmut) glioma. IDHwt gliomas showed a stable epigenome over time with relatively low levels of global methylation. The epigenome of IDHmut gliomas showed initial high levels of genome-wide DNA methylation that was progressively reduced to levels similar to those of IDHwt tumors. Integration of epigenomics, gene expression, and functional genomics identified HOXD13 as a master regulator of IDHmut astrocytoma evolution. Furthermore, relapse of IDHmut tumors was accompanied by histologic progression that was associated with survival, as validated in an independent cohort. Finally, the initial cell composition of the tumor microenvironment varied between IDHwt and IDHmut tumors and changed differentially following treatment, suggesting increased neoangiogenesis and T-cell infiltration upon treatment of IDHmut gliomas. This study provides one of the largest cohorts of paired longitudinal glioma samples with epigenomic, transcriptomic, and genomic profiling and suggests that treatment of IDHmut glioma is associated with epigenomic evolution toward an IDHwt-like phenotype. SIGNIFICANCE: Standard treatments are related to loss of DNA methylation in IDHmut glioma, resulting in epigenetic activation of genes associated with tumor progression and alterations in the microenvironment that resemble treatment-naïve IDHwt glioma.

Integrated molecular and multiparametric MRI mapping of high-grade glioma identifies regional biologic signatures.

Sampling restrictions have hindered the comprehensive study of invasive non-enhancing (NE) high-grade glioma (HGG) cell populations driving tumor progression. Here, we present an integrated multi-omic analysis of spatially matched molecular and multi-parametric magnetic resonance imaging (MRI) profiling across 313 multi-regional tumor biopsies, including 111 from the NE, across 68 HGG patients. Whole exome and RNA sequencing uncover unique genomic alterations to unresectable invasive NE tumor, including subclonal events, which inform genomic models predictive of geographic evolution. Infiltrative NE tumor is alternatively enriched with tumor cells exhibiting neuronal or glycolytic/plurimetabolic cellular states, two principal transcriptomic pathway-based glioma subtypes, which respectively demonstrate abundant private mutations or enrichment in immune cell signatures. These NE phenotypes are non-invasively identified through normalized K2 imaging signatures, which discern cell size heterogeneity on dynamic susceptibility contrast (DSC)-MRI. NE tumor populations predicted to display increased cellular proliferation by mean diffusivity (MD) MRI metrics are uniquely associated with EGFR amplification and CDKN2A homozygous deletion. The biophysical mapping of infiltrative HGG potentially enables the clinical recognition of tumor subpopulations with aggressive molecular signatures driving tumor progression, thereby informing precision medicine targeting.

Integrative multi-omics networks identify PKCδ and DNA-PK as master kinases of glioblastoma subtypes and guide targeted cancer therapy.

Despite producing a panoply of potential cancer-specific targets, the proteogenomic characterization of human tumors has yet to demonstrate value for precision cancer medicine. Integrative multi-omics using a machine-learning network identified master kinases responsible for effecting phenotypic hallmarks of functional glioblastoma subtypes. In subtype-matched patient-derived models, we validated PKCδ and DNA-PK as master kinases of glycolytic/plurimetabolic and proliferative/progenitor subtypes, respectively, and qualified the kinases as potent and actionable glioblastoma subtype-specific therapeutic targets. Glioblastoma subtypes were associated with clinical and radiomics features, orthogonally validated by proteomics, phospho-proteomics, metabolomics, lipidomics and acetylomics analyses, and recapitulated in pediatric glioma, breast and lung squamous cell carcinoma, including subtype specificity of PKCδ and DNA-PK activity. We developed a probabilistic classification tool that performs optimally with RNA from frozen and paraffin-embedded tissues, which can be used to evaluate the association of therapeutic response with glioblastoma subtypes and to inform patient selection in prospective clinical trials.

LZTR1 Mutation Mediates Oncogenesis through Stabilization of EGFR and AXL.

LZTR1 is the substrate-specific adaptor of a CUL3-dependent ubiquitin ligase frequently mutated in sporadic and syndromic cancer. We combined biochemical and genetic studies to identify LZTR1 substrates and interrogated their tumor-driving function in the context of LZTR1 loss-of-function mutations. Unbiased screens converged on EGFR and AXL receptor tyrosine kinases as LZTR1 interactors targeted for ubiquitin-dependent degradation in the lysosome. Pathogenic cancer-associated mutations of LZTR1 failed to promote EGFR and AXL degradation, resulting in dysregulated growth factor signaling. Conditional inactivation of Lztr1 and Cdkn2a in the mouse nervous system caused tumors in the peripheral nervous system including schwannoma-like tumors, thus recapitulating aspects of schwannomatosis, the prototype tumor predisposition syndrome sustained by LZTR1 germline mutations. Lztr1- and Cdkn2a-deleted tumors aberrantly accumulated EGFR and AXL and exhibited specific vulnerability to EGFR and AXL coinhibition. These findings explain tumorigenesis by LZTR1 inactivation and offer therapeutic opportunities to patients with LZTR1-mutant cancer. SIGNIFICANCE: EGFR and AXL are substrates of LZTR1-CUL3 ubiquitin ligase. The frequent somatic and germline mutations of LZTR1 in human cancer cause EGFR and AXL accumulation and deregulated signaling. LZTR1-mutant tumors show vulnerability to concurrent inhibition of EGFR and AXL, thus providing precision targeting to patients affected by LZTR1-mutant cancer. This article is highlighted in the In This Issue feature, p. 517.

Glioma progression is shaped by genetic evolution and microenvironment interactions.

The factors driving therapy resistance in diffuse glioma remain poorly understood. To identify treatment-associated cellular and genetic changes, we analyzed RNA and/or DNA sequencing data from the temporally separated tumor pairs of 304 adult patients with isocitrate dehydrogenase (IDH)-wild-type and IDH-mutant glioma. Tumors recurred in distinct manners that were dependent on IDH mutation status and attributable to changes in histological feature composition, somatic alterations, and microenvironment interactions. Hypermutation and acquired CDKN2A deletions were associated with an increase in proliferating neoplastic cells at recurrence in both glioma subtypes, reflecting active tumor growth. IDH-wild-type tumors were more invasive at recurrence, and their neoplastic cells exhibited increased expression of neuronal signaling programs that reflected a possible role for neuronal interactions in promoting glioma progression. Mesenchymal transition was associated with the presence of a myeloid cell state defined by specific ligand-receptor interactions with neoplastic cells. Collectively, these recurrence-associated phenotypes represent potential targets to alter disease progression.

Regulated interaction of ID2 with the anaphase-promoting complex links progression through mitosis with reactivation of cell-type-specific transcription.

Tissue-specific transcriptional activity is silenced in mitotic cells but it remains unclear whether the mitotic regulatory machinery interacts with tissue-specific transcriptional programs. We show that such cross-talk involves the controlled interaction between core subunits of the anaphase-promoting complex (APC) and the ID2 substrate. The N-terminus of ID2 is independently and structurally compatible with a pocket composed of core APC/C subunits that may optimally orient ID2 onto the APCCDH1 complex. Phosphorylation of serine-5 by CDK1 prevented the association of ID2 with core APC, impaired ubiquitylation and stabilized ID2 protein at the mitosis-G1 transition leading to inhibition of basic Helix-Loop-Helix (bHLH)-mediated transcription. The serine-5 phospho-mimetic mutant of ID2 that inefficiently bound core APC remained stable during mitosis, delayed exit from mitosis and reloading of bHLH transcription factors on chromatin. It also locked cells into a "mitotic stem cell" transcriptional state resembling the pluripotent program of embryonic stem cells. The substrates of APCCDH1 SKP2 and Cyclin B1 share with ID2 the phosphorylation-dependent, D-box-independent interaction with core APC. These results reveal a new layer of control of the mechanism by which substrates are recognized by APC.

Inhibition of ERK/MAPK signaling as potential therapy to prevent optic pathway glioma in infants with neurofibromatosis type 1.

Pediatric low-grade gliomas (pLGGs) arise primarily at early stages of development. The molecular mechanisms of pLGG gliomagenesis are unclear, as is the progenitor cell of origin. In this issue of Developmental Cell, Jecrois et al. show that NF1-associated optic pathway gliomas originate from migrating glial progenitors that have distinct MEK/ERK dependency.

Structural Model for Recruitment of RIT1 to the LZTR1 E3 Ligase: Evidences from an Integrated Computational Approach.

Leucine-zipper transcription regulator 1 (LZTR1) is a highly mutated tumor suppressor gene, involved in the pathogenesis of several cancer types and developmental disorders. In proteasomal degradation, it acts as an adaptor protein responsible for the recognition and recruitment of substrates to be ubiquitinated in Cullin3-RING ligase E3 (CRL3) machinery. LZTR1 belongs to the BTB-Kelch family, a multi-domain protein where the Kelch propeller plays as the substrate recognition region and for which no experimental structure has been solved. Recently, large effort mutational analyses pointed to the role of disease-associated LZTR1 mutations in the RAS/MAPK signaling pathway and RIT1, a small Ras-related GTPase protein, has been identified by mass spectroscopy to interact with LZTR1. Hence, a better understanding of native structure, molecular mechanism, and substrate specificity would help clarifying the role of LZTR1 in pathological diseases, thus promoting advancement in the development of novel therapeutic strategies. Here, we address the interaction model between adaptor LZTR1 and substrate RIT1 by applying an integrated computational approach, including molecular modeling and docking techniques. We observe that the interaction model LZTR1-RIT1 is stabilized by an electrostatic bond network established between the two protein surfaces, which is reminiscent of homologous ubiquitin ligases complexes. Then, running MD simulations, we characterize differential conformational dynamics of the multi-domain LZTR1, offering interesting implications on the mechanistic role of specific point mutations. We identify G248R and R283Q as damaging mutations involved in the recognition process of the substrate RIT1 and R412C as a possible allosteric mutation from the Kelch to the C-term BTB-domain. Our findings provide important structural insights on targeting CRL3s for drug discovery.

Pathway-based classification of glioblastoma uncovers a mitochondrial subtype with therapeutic vulnerabilities.

The transcriptomic classification of glioblastoma (GBM) has failed to predict survival and therapeutic vulnerabilities. A computational approach for unbiased identification of core biological traits of single cells and bulk tumors uncovered four tumor cell states and GBM subtypes distributed along neurodevelopmental and metabolic axes, classified as proliferative/progenitor, neuronal, mitochondrial and glycolytic/plurimetabolic. Each subtype was enriched with biologically coherent multiomic features. Mitochondrial GBM was associated with the most favorable clinical outcome. It relied exclusively on oxidative phosphorylation for energy production, whereas the glycolytic/plurimetabolic subtype was sustained by aerobic glycolysis and amino acid and lipid metabolism. Deletion of the glucose-proton symporter SLC45A1 was the truncal alteration most significantly associated with mitochondrial GBM, and the reintroduction of SLC45A1 in mitochondrial glioma cells induced acidification and loss of fitness. Mitochondrial, but not glycolytic/plurimetabolic, GBM exhibited marked vulnerability to inhibitors of oxidative phosphorylation. The pathway-based classification of GBM informs survival and enables precision targeting of cancer metabolism.

Dysregulation of Principal Cell miRNAs Facilitates Epigenetic Regulation of AQP2 and Results in Nephrogenic Diabetes Insipidus.

BACKGROUND: MicroRNAs (miRNAs), formed by cleavage of pre-microRNA by the endoribonuclease Dicer, are critical modulators of cell function by post-transcriptionally regulating gene expression. METHODS: Selective ablation of Dicer in AQP2-expressing cells (DicerAQP2Cre+ mice) was used to investigate the role of miRNAs in the kidney collecting duct of mice. RESULTS: The mice had severe polyuria and nephrogenic diabetes insipidus, potentially due to greatly reduced AQP2 and AQP4 levels. Although epithelial sodium channel levels were decreased in cortex and increased in inner medulla, amiloride-sensitive sodium reabsorption was equivalent in DicerAQP2Cre+ mice and controls. Small-RNA sequencing and proteomic analysis revealed 31 and 178 significantly regulated miRNAs and proteins, respectively. Integrated bioinformatic analysis of the miRNAome and proteome suggested alterations in the epigenetic machinery and various transcription factors regulating AQP2 expression in DicerAQP2Cre+ mice. The expression profile and function of three miRNAs (miR-7688-5p, miR-8114, and miR-409-3p) whose predicted targets were involved in epigenetic control (Phf2, Kdm5c, and Kdm4a) or transcriptional regulation (GATA3, GATA2, and ELF3) of AQP2 were validated. Luciferase assays could not demonstrate direct interaction of AQP2 or the three potential transcription factors with miR-7688-5p, miR-8114, and miR-409-3p. However, transfection of respective miRNA mimics reduced AQP2 expression. Chromatin immunoprecipitation assays demonstrated decreased Phf2 and significantly increased Kdm5c interactions at the Aqp2 gene promoter in DicerAQP2Cre+ mice, resulting in decreased RNA Pol II association. CONCLUSIONS: Novel evidence indicates miRNA-mediated epigenetic regulation of AQP2 expression.

Modulation of Nogo receptor 1 expression orchestrates myelin-associated infiltration of glioblastoma.

As the clinical failure of glioblastoma treatment is attributed by multiple components, including myelin-associated infiltration, assessment of the molecular mechanisms underlying such process and identification of the infiltrating cells have been the primary objectives in glioblastoma research. Here, we adopted radiogenomic analysis to screen for functionally relevant genes that orchestrate the process of glioma cell infiltration through myelin and promote glioblastoma aggressiveness. The receptor of the Nogo ligand (NgR1) was selected as the top candidate through Differentially Expressed Genes (DEG) and Gene Ontology (GO) enrichment analysis. Gain and loss of function studies on NgR1 elucidated its underlying molecular importance in suppressing myelin-associated infiltration in vitro and in vivo. The migratory ability of glioblastoma cells on myelin is reversibly modulated by NgR1 during differentiation and dedifferentiation process through deubiquitinating activity of USP1, which inhibits the degradation of ID1 to downregulate NgR1 expression. Furthermore, pimozide, a well-known antipsychotic drug, upregulates NgR1 by post-translational targeting of USP1, which sensitizes glioma stem cells to myelin inhibition and suppresses myelin-associated infiltration in vivo. In primary human glioblastoma, downregulation of NgR1 expression is associated with highly infiltrative characteristics and poor survival. Together, our findings reveal that loss of NgR1 drives myelin-associated infiltration of glioblastoma and suggest that novel therapeutic strategies aimed at reactivating expression of NgR1 will improve the clinical outcome of glioblastoma patients.

A map of tumor-host interactions in glioma at single-cell resolution.

BACKGROUND: Single-cell RNA sequencing is the reference technique for characterizing the heterogeneity of the tumor microenvironment. The composition of the various cell types making up the microenvironment can significantly affect the way in which the immune system activates cancer rejection mechanisms. Understanding the cross-talk signals between immune cells and cancer cells is of fundamental importance for the identification of immuno-oncology therapeutic targets. RESULTS: We present a novel method, single-cell Tumor-Host Interaction tool (scTHI), to identify significantly activated ligand-receptor interactions across clusters of cells from single-cell RNA sequencing data. We apply our approach to uncover the ligand-receptor interactions in glioma using 6 publicly available human glioma datasets encompassing 57,060 gene expression profiles from 71 patients. By leveraging this large-scale collection we show that unexpected cross-talk partners are highly conserved across different datasets in the majority of the tumor samples. This suggests that shared cross-talk mechanisms exist in glioma. CONCLUSIONS: Our results provide a complete map of the active tumor-host interaction pairs in glioma that can be therapeutically exploited to reduce the immunosuppressive action of the microenvironment in brain tumor.

Temporospatial genomic profiling in glioblastoma identifies commonly altered core pathways underlying tumor progression.

BACKGROUND: Tumor heterogeneity underlies resistance and disease progression in glioblastoma (GBM), and tumors most commonly recur adjacent to the surgical resection margins in contrast non-enhancing (NE) regions. To date, no targeted therapies have meaningfully altered overall patient survival in the up-front setting. The aim of this study was to characterize intratumoral heterogeneity in recurrent GBM using bulk samples from primary resection and recurrent samples taken from contrast-enhancing (EN) and contrast NE regions. METHODS: Whole exome and RNA sequencing were performed on matched bulk primary and multiple recurrent EN and NE tumor samples from 16 GBM patients who received standard of care treatment alone or in combination with investigational clinical trial regimens. RESULTS: Private mutations emerge across multi-region sampling in recurrent tumors. Genomic clonal analysis revealed increased enrichment in gene alterations regulating the G2M checkpoint, Kras signaling, Wnt signaling, and DNA repair in recurrent disease. Subsequent functional studies identified augmented PI3K/AKT transcriptional and protein activity throughout progression, validated by phospho-protein levels. Moreover, a mesenchymal transcriptional signature was observed in recurrent EN regions, which differed from the proneural signature in recurrent NE regions. CONCLUSIONS: Subclonal populations observed within bulk resected primary GBMs transcriptionally evolve across tumor recurrence (EN and NE regions) and exhibit aberrant gene expression of common signaling pathways that persist despite standard or targeted therapy. Our findings provide evidence that there are both adaptive and clonally mediated dependencies of GBM on key pathways, such as the PI3K/AKT axis, for survival across recurrences.

Proline Hydroxylation Primes Protein Kinases for Autophosphorylation and Activation.

Activation of dual-specificity tyrosine-phosphorylation-regulated kinases 1A and 1B (DYRK1A and DYRK1B) requires prolyl hydroxylation by PHD1 prolyl hydroxylase. Prolyl hydroxylation of DYRK1 initiates a cascade of events leading to the release of molecular constraints on von Hippel-Lindau (VHL) ubiquitin ligase tumor suppressor function. However, the proline residue of DYRK1 targeted by hydroxylation and the role of prolyl hydroxylation in tyrosine autophosphorylation of DYRK1 are unknown. We found that a highly conserved proline in the CMGC insert of the DYRK1 kinase domain is hydroxylated by PHD1, and this event precedes tyrosine autophosphorylation. Mutation of the hydroxylation acceptor proline precludes tyrosine autophosphorylation and folding of DYRK1, resulting in a kinase unable to preserve VHL function and lacking glioma suppression activity. The consensus proline sequence is shared by most CMGC kinases, and prolyl hydroxylation is essential for catalytic activation. Thus, formation of prolyl-hydroxylated intermediates is a novel mechanism of kinase maturation and likely a general mechanism of regulation of CMGC kinases in eukaryotes.

Exploring the Molecular Crosstalk between Pancreatic Bud and Mesenchyme in Embryogenesis: Novel Signals Involved.

Pancreatic organogenesis is a multistep process that requires the cooperation of several signaling pathways. In this context, the role of pancreatic mesenchyme is important to define the epithelium development; nevertheless, the precise space-temporal signaling activation still needs to be clarified. This study reports a dissection of the pancreatic embryogenesis, highlighting the molecular network surrounding the epithelium-mesenchyme interaction. To investigate this crosstalk, pancreatic epithelium and surrounding mesenchyme, at embryonic day 10.5, were collected through laser capture microdissection (LCM) and characterized based on their global gene expression. We performed a bioinformatic analysis to hypothesize crosstalk interactions, validating the most promising genes and verifying the precise localization of their expression in the compartments, by RNA in situ hybridization (ISH). Our analyses pointed out also the c-Met gene, a very well-known factor involved in stimulating motility, morphogenesis, and organ regeneration. We also highlighted the potential crosstalk between Versican (Vcan) and Syndecan4 (Sdc4) since these genes are involved in pancreatic tissue repair, strengthening the concept that the same signaling pathways required during pancreatic embryogenesis are also involved in tissue repair. This finding leads to novel strategies for obtaining functional pancreatic stem cells for cell replacement therapies.

Circulating microRNAs expression profile in newly diagnosed and imatinib treated chronic phase - chronic myeloid leukemia.

Chronic myeloid leukemia (CML) is a stem cell derived malignant disorder result of translocation t(9;22)(q34;q11) called Philadelphia chromosome (Ph+). microRNAS (miRNAs) are involved in several biological processes, altering the progression of various pathologies, including CML. This study evaluated whether circulating miRNAs display differential expression profiles in peripheral blood of CML-Chronic Phase (CML-CP) patients newly diagnosed in comparison with CML-CP treated with imatinib. We obtained peripheral blood samples from CML-CP Ph+ patients divided among group 1 (untreated newly diagnosed) and group 2 (treated with imatinib). A pool of total leukocytes from healthy donors was considered as control group. Expression analyses were performed for 768 miRNAs by RT-qPCR array. Bioinformatic tools were used to identify significant pathways and interaction networks. We found 80 deregulated miRNAs between the groups and, according to bioinformatic analysis, they are involved in different pathways, including molecular mechanisms of cancer. The study allows better understanding of disease molecular behavior, and it is useful for possible monitoring CML treatment and prognostic biomarkers identification.

The molecular landscape of glioma in patients with Neurofibromatosis 1.

Neurofibromatosis type 1 (NF1) is a common tumor predisposition syndrome in which glioma is one of the prevalent tumors. Gliomagenesis in NF1 results in a heterogeneous spectrum of low- to high-grade neoplasms occurring during the entire lifespan of patients. The pattern of genetic and epigenetic alterations of glioma that develops in NF1 patients and the similarities with sporadic glioma remain unknown. Here, we present the molecular landscape of low- and high-grade gliomas in patients affected by NF1 (NF1-glioma). We found that the predisposing germline mutation of the NF1 gene was frequently converted to homozygosity and the somatic mutational load of NF1-glioma was influenced by age and grade. High-grade tumors harbored genetic alterations of TP53 and CDKN2A, frequent mutations of ATRX associated with Alternative Lengthening of Telomere, and were enriched in genetic alterations of transcription/chromatin regulation and PI3 kinase pathways. Low-grade tumors exhibited fewer mutations that were over-represented in genes of the MAP kinase pathway. Approximately 50% of low-grade NF1-gliomas displayed an immune signature, T lymphocyte infiltrates, and increased neo-antigen load. DNA methylation assigned NF1-glioma to LGm6, a poorly defined Isocitrate Dehydrogenase 1 wild-type subgroup enriched with ATRX mutations. Thus, the profiling of NF1-glioma defined a distinct landscape that recapitulates a subset of sporadic tumors.

Transcriptional landscape of mouse-aged ovaries reveals a unique set of non-coding RNAs associated with physiological and environmental ovarian dysfunctions.

The progressive and physiological decline in ovarian function depends on the rate of follicular loss by atresia, contributing to the reduction in ovarian reserve. Genetics and environmental factors play important roles in ovarian senescence and in the onset of ovarian dysfunctions such as diminished ovarian reserve. A better understanding of the mechanisms underlying ovarian aging and their regulation by genetic and environmental factors is needed to evaluate ovarian reserve and to predict fertility potential by identification of more accurate and less invasive markers. We report transcriptomic data (i) implicating novel (e.g. EIF2 signalling) and well-known pathways (e.g. TGFß signalling), and (ii) defining a unique set of non-coding RNA (ncRNA), both associated with ovarian function. The latter includes miRNAs (e.g. Mir143 and Mir145), snoRNAs (e.g. Snord16a and Snora34), and one lncRNA (Gas5), which are differentially expressed in middle-aged ovaries (12 months) vs young-aged (3 months) from CD1 mice. Experimental analysis confirms that ovary lifespan varies across genetic backgrounds in mice and, genetics influences the response to environmental perturbations such as diet. Moreover, the identified ncRNAs were verified in a model of reproductive dysfunction promoted by the environmental toxicant ethylenthiourea. We also report the increase of miRNA143 and miRNA145 in follicular fluid of women with diminished ovarian reserve. Their levels inversely correlate with the hormonal profile and with the number of the oocytes recruited upon hormonal stimulation. Overall, we report a transcriptomic signature for ovarian dysfunction in vivo that provides a valuable resource for translational research in human reproductive aging.