A prognostic neural epigenetic signature in high-grade glioma.

Neural-tumor interactions drive glioma growth as evidenced in preclinical models, but clinical validation is limited. We present an epigenetically defined neural signature of glioblastoma that independently predicts patients' survival. We use reference signatures of neural cells to deconvolve tumor DNA and classify samples into low- or high-neural tumors. High-neural glioblastomas exhibit hypomethylated CpG sites and upregulation of genes associated with synaptic integration. Single-cell transcriptomic analysis reveals a high abundance of malignant stemcell-like cells in high-neural glioblastoma, primarily of the neural lineage. These cells are further classified as neural-progenitor-cell-like, astrocyte-like and oligodendrocyte-progenitor-like, alongside oligodendrocytes and excitatory neurons. In line with these findings, high-neural glioblastoma cells engender neuron-to-glioma synapse formation in vitro and in vivo and show an unfavorable survival after xenografting. In patients, a high-neural signature is associated with decreased overall and progression-free survival. High-neural tumors also exhibit increased functional connectivity in magnetencephalography and resting-state magnet resonance imaging and can be detected via DNA analytes and brain-derived neurotrophic factor in patients' plasma. The prognostic importance of the neural signature was further validated in patients diagnosed with diffuse midline glioma. Our study presents an epigenetically defined malignant neural signature in high-grade gliomas that is prognostically relevant. High-neural gliomas likely require a maximized surgical resection approach for improved outcomes.

Conumee 2.0: enhanced copy-number variation analysis from DNA methylation arrays for humans and mice.

MOTIVATION: Copy-number variations (CNVs) are common genetic alterations in cancer and their detection may impact tumor classification and therapeutic decisions. However, detection of clinically relevant large and focal CNVs remains challenging when sample material or resources are limited. This has motivated us to create a software tool to infer CNVs from DNA methylation arrays which are often generated as part of clinical routines and in research settings. RESULTS: We present our R package, conumee 2.0, that combines tangent normalization, an adjustable genomic binning heuristic, and weighted circular binary segmentation to utilize DNA methylation arrays for CNV analysis and mitigate technical biases and batch effects. Segmentation results were validated in a lung squamous cell carcinoma dataset from TCGA (n = 367 samples) by comparison to segmentations derived from genotyping arrays (Pearson's correlation coefficient of 0.91). We further introduce a segmented block bootstrapping approach to detect focal alternations that achieved 60.9% sensitivity and 98.6% specificity for deletions affecting CDKN2A/B (60.0% and 96.9% for RB1, respectively) in a low-grade glioma cohort from TCGA (n = 239 samples). Finally, our tool provides functionality to detect and summarize CNVs across large sample cohorts. AVAILABILITY AND IMPLEMENTATION: Conumee 2.0 is available under open-source license at: https://github.com/hovestadtlab/conumee2.

Temporal change of DNA methylation subclasses between matched newly diagnosed and recurrent glioblastoma.

The longitudinal transition of phenotypes is pivotal in glioblastoma treatment resistance and DNA methylation emerged as an important tool for classifying glioblastoma phenotypes. We aimed to characterize DNA methylation subclass heterogeneity during progression and assess its clinical impact. Matched tissues from 47 glioblastoma patients were subjected to DNA methylation profiling, including CpG-site alterations, tissue and serum deconvolution, mass spectrometry, and immunoassay. Effects of clinical characteristics on temporal changes and outcomes were studied. Among 47 patients, 8 (17.0%) had non-matching classifications at recurrence. In the remaining 39 cases, 28.2% showed dominant DNA methylation subclass transitions, with 72.7% being a mesenchymal subclass. In general, glioblastomas with a subclass transition showed upregulated metabolic processes. Newly diagnosed glioblastomas with mesenchymal transition displayed increased stem cell-like states and decreased immune components at diagnosis and exhibited elevated immune signatures and cytokine levels in serum. In contrast, tissue of recurrent glioblastomas with mesenchymal transition showed increased immune components but decreased stem cell-like states. Survival analyses revealed comparable outcomes for patients with and without subclass transitions. This study demonstrates a temporal heterogeneity of DNA methylation subclasses in 28.2% of glioblastomas, not impacting patient survival. Changes in cell state composition associated with subclass transition may be crucial for recurrent glioblastoma targeted therapies.

Epigenetic neural glioblastoma enhances synaptic integration and predicts therapeutic vulnerability.

Neural-tumor interactions drive glioma growth as evidenced in preclinical models, but clinical validation is nascent. We present an epigenetically defined neural signature of glioblastoma that independently affects patients' survival. We use reference signatures of neural cells to deconvolve tumor DNA and classify samples into low- or high-neural tumors. High-neural glioblastomas exhibit hypomethylated CpG sites and upregulation of genes associated with synaptic integration. Single-cell transcriptomic analysis reveals high abundance of stem cell-like malignant cells classified as oligodendrocyte precursor and neural precursor cell-like in high-neural glioblastoma. High-neural glioblastoma cells engender neuron-to-glioma synapse formation in vitro and in vivo and show an unfavorable survival after xenografting. In patients, a high-neural signature associates with decreased survival as well as increased functional connectivity and can be detected via DNA analytes and brain-derived neurotrophic factor in plasma. Our study presents an epigenetically defined malignant neural signature in high-grade gliomas that is prognostically relevant.

Mapping pediatric brain tumors to their origins in the developing cerebellum.

BACKGROUND: Distinguishing the cellular origins of childhood brain tumors is key for understanding tumor initiation and identifying lineage-restricted, tumor-specific therapeutic targets. Previous strategies to map the cell-of-origin typically involved comparing human tumors to murine embryonal tissues, which is potentially limited due to species-specific differences. The aim of this study was to unravel the cellular origins of the 3 most common pediatric brain tumors, ependymoma, pilocytic astrocytoma, and medulloblastoma, using a developing human cerebellar atlas. METHODS: We used a single-nucleus atlas of the normal developing human cerebellum consisting of 176 645 cells as a reference for an in-depth comparison to 4416 bulk and single-cell transcriptome tumor datasets, using gene set variation analysis, correlation, and single-cell matching techniques. RESULTS: We find that the astroglial cerebellar lineage is potentially the origin for posterior fossa ependymomas. We propose that infratentorial pilocytic astrocytomas originate from the oligodendrocyte lineage and MHC II genes are specifically enriched in these tumors. We confirm that SHH and Group 3/4 medulloblastomas originate from the granule cell and unipolar brush cell lineages. Radiation-induced gliomas stem from cerebellar glial lineages and demonstrate distinct origins from the primary medulloblastoma. We identify tumor genes that are expressed in the cerebellar lineage of origin, and genes that are tumor specific; both gene sets represent promising therapeutic targets for future study. CONCLUSION: Based on our results, individual cells within a tumor may resemble different cell types along a restricted developmental lineage. Therefore, we suggest that tumors can arise from multiple cellular states along the cerebellar "lineage of origin."

Modeling epigenetic lesions that cause gliomas.

Epigenetic lesions that disrupt regulatory elements represent potential cancer drivers. However, we lack experimental models for validating their tumorigenic impact. Here, we model aberrations arising in isocitrate dehydrogenase-mutant gliomas, which exhibit DNA hypermethylation. We focus on a CTCF insulator near the PDGFRA oncogene that is recurrently disrupted by methylation in these tumors. We demonstrate that disruption of the syntenic insulator in mouse oligodendrocyte progenitor cells (OPCs) allows an OPC-specific enhancer to contact and induce Pdgfra, thereby increasing proliferation. We show that a second lesion, methylation-dependent silencing of the Cdkn2a tumor suppressor, cooperates with insulator loss in OPCs. Coordinate inactivation of the Pdgfra insulator and Cdkn2a drives gliomagenesis in vivo. Despite locus synteny, the insulator is CpG-rich only in humans, a feature that may confer human glioma risk but complicates mouse modeling. Our study demonstrates the capacity of recurrent epigenetic lesions to drive OPC proliferation in vitro and gliomagenesis in vivo.

Ultraviolet radiation shapes dendritic cell leukaemia transformation in the skin.

Tumours most often arise from progression of precursor clones within a single anatomical niche. In the bone marrow, clonal progenitors can undergo malignant transformation to acute leukaemia, or differentiate into immune cells that contribute to disease pathology in peripheral tissues1-4. Outside the marrow, these clones are potentially exposed to a variety of tissue-specific mutational processes, although the consequences of this are unclear. Here we investigate the development of blastic plasmacytoid dendritic cell neoplasm (BPDCN)-an unusual form of acute leukaemia that often presents with malignant cells isolated to the skin5. Using tumour phylogenomics and single-cell transcriptomics with genotyping, we find that BPDCN arises from clonal (premalignant) haematopoietic precursors in the bone marrow. We observe that BPDCN skin tumours first develop at sun-exposed anatomical sites and are distinguished by clonally expanded mutations induced by ultraviolet (UV) radiation. A reconstruction of tumour phylogenies reveals that UV damage can precede the acquisition of alterations associated with malignant transformation, implicating sun exposure of plasmacytoid dendritic cells or committed precursors during BPDCN pathogenesis. Functionally, we find that loss-of-function mutations in Tet2, the most common premalignant alteration in BPDCN, confer resistance to UV-induced cell death in plasmacytoid, but not conventional, dendritic cells, suggesting a context-dependent tumour-suppressive role for TET2. These findings demonstrate how tissue-specific environmental exposures at distant anatomical sites can shape the evolution of premalignant clones to disseminated cancer.

Multiomic neuropathology improves diagnostic accuracy in pediatric neuro-oncology.

The large diversity of central nervous system (CNS) tumor types in children and adolescents results in disparate patient outcomes and renders accurate diagnosis challenging. In this study, we prospectively integrated DNA methylation profiling and targeted gene panel sequencing with blinded neuropathological reference diagnostics for a population-based cohort of more than 1,200 newly diagnosed pediatric patients with CNS tumors, to assess their utility in routine neuropathology. We show that the multi-omic integration increased diagnostic accuracy in a substantial proportion of patients through annotation to a refining DNA methylation class (50%), detection of diagnostic or therapeutically relevant genetic alterations (47%) or identification of cancer predisposition syndromes (10%). Discrepant results by neuropathological WHO-based and DNA methylation-based classification (30%) were enriched in histological high-grade gliomas, implicating relevance for current clinical patient management in 5% of all patients. Follow-up (median 2.5 years) suggests improved survival for patients with histological high-grade gliomas displaying lower-grade molecular profiles. These results provide preliminary evidence of the utility of integrating multi-omics in neuropathology for pediatric neuro-oncology.

Unified rhombic lip origins of group 3 and group 4 medulloblastoma.

Medulloblastoma, a malignant childhood cerebellar tumour, segregates molecularly into biologically distinct subgroups, suggesting that a personalized approach to therapy would be beneficial1. Mouse modelling and cross-species genomics have provided increasing evidence of discrete, subgroup-specific developmental origins2. However, the anatomical and cellular complexity of developing human tissues3-particularly within the rhombic lip germinal zone, which produces all glutamatergic neuronal lineages before internalization into the cerebellar nodulus-makes it difficult to validate previous inferences that were derived from studies in mice. Here we use multi-omics to resolve the origins of medulloblastoma subgroups in the developing human cerebellum. Molecular signatures encoded within a human rhombic-lip-derived lineage trajectory aligned with photoreceptor and unipolar brush cell expression profiles that are maintained in group 3 and group 4 medulloblastoma, suggesting a convergent basis. A systematic diagnostic-imaging review of a prospective institutional cohort localized the putative anatomical origins of group 3 and group 4 tumours to the nodulus. Our results connect the molecular and phenotypic features of clinically challenging medulloblastoma subgroups to their unified beginnings in the rhombic lip in the early stages of human development.

FOXR2 Is an Epigenetically Regulated Pan-Cancer Oncogene That Activates ETS Transcriptional Circuits.

Forkhead box R2 (FOXR2) is a forkhead transcription factor located on the X chromosome whose expression is normally restricted to the testis. In this study, we performed a pan-cancer analysis of FOXR2 activation across more than 10,000 adult and pediatric cancer samples and found FOXR2 to be aberrantly upregulated in 70% of all cancer types and 8% of all individual tumors. The majority of tumors (78%) aberrantly expressed FOXR2 through a previously undescribed epigenetic mechanism that involves hypomethylation of a novel promoter, which was functionally validated as necessary for FOXR2 expression and proliferation in FOXR2-expressing cancer cells. FOXR2 promoted tumor growth across multiple cancer lineages and co-opted ETS family transcription circuits across cancers. Taken together, this study identifies FOXR2 as a potent and ubiquitous oncogene that is epigenetically activated across the majority of human cancers. The identification of hijacking of ETS transcription circuits by FOXR2 extends the mechanisms known to active ETS transcription factors and highlights how transcription factor families cooperate to enhance tumorigenesis. SIGNIFICANCE: This work identifies a novel promoter that drives aberrant FOXR2 expression and delineates FOXR2 as a pan-cancer oncogene that specifically activates ETS transcriptional circuits across human cancers. See related commentary by Liu and Northcott, p. 2977.

Systematic detection of m6A-modified transcripts at single-molecule and single-cell resolution.

Epigenetic modifications control the stability and translation of mRNA molecules. Here, we present a microscopy-based platform for quantifying modified RNA molecules and for relating the modification patterns to single-cell phenotypes. We directly capture mRNAs from cell lysates on oligo-dT-coated coverslips, then visually detect and sequence individual m6A-immunolabled transcripts without amplification. Integration of a nanoscale device enabled us to isolate single cells on the platform, and thereby relate single-cell m6A modification states to gene expression signatures and cell surface markers. Application of the platform to MUTZ3 leukemia cells revealed a marked reduction in cellular m6A levels as CD34+ leukemic progenitors differentiate to CD14+ myeloid cells. We then coupled single-molecule m6A detection with fluorescence in situ hybridization (FISH) to relate mRNA and m6A levels of individual genes to single-cell phenotypes. This single-cell multi-modal assay suite can empower investigations of RNA modifications in rare populations and single cells.

Mutational mechanisms shaping the coding and noncoding genome of germinal center derived B-cell lymphomas.

B cells have the unique property to somatically alter their immunoglobulin (IG) genes by V(D)J recombination, somatic hypermutation (SHM) and class-switch recombination (CSR). Aberrant targeting of these mechanisms is implicated in lymphomagenesis, but the mutational processes are poorly understood. By performing whole genome and transcriptome sequencing of 181 germinal center derived B-cell lymphomas (gcBCL) we identified distinct mutational signatures linked to SHM and CSR. We show that not only SHM, but presumably also CSR causes off-target mutations in non-IG genes. Kataegis clusters with high mutational density mainly affected early replicating regions and were enriched for SHM- and CSR-mediated off-target mutations. Moreover, they often co-occurred in loci physically interacting in the nucleus, suggesting that mutation hotspots promote increased mutation targeting of spatially co-localized loci (termed hypermutation by proxy). Only around 1% of somatic small variants were in protein coding sequences, but in about half of the driver genes, a contribution of B-cell specific mutational processes to their mutations was found. The B-cell-specific mutational processes contribute to both lymphoma initiation and intratumoral heterogeneity. Overall, we demonstrate that mutational processes involved in the development of gcBCL are more complex than previously appreciated, and that B cell-specific mutational processes contribute via diverse mechanisms to lymphomagenesis.

Extended-representation bisulfite sequencing of gene regulatory elements in multiplexed samples and single cells.

The biological roles of DNA methylation have been elucidated by profiling methods based on whole-genome or reduced-representation bisulfite sequencing, but these approaches do not efficiently survey the vast numbers of non-coding regulatory elements in mammalian genomes. Here we present an extended-representation bisulfite sequencing (XRBS) method for targeted profiling of DNA methylation. Our design strikes a balance between expanding coverage of regulatory elements and reproducibly enriching informative CpG dinucleotides in promoters, enhancers and CTCF binding sites. Barcoded DNA fragments are pooled before bisulfite conversion, allowing multiplex processing and technical consistency in low-input samples. Application of XRBS to single leukemia cells enabled us to evaluate genetic copy number variations and methylation variability across individual cells. Our analysis highlights heterochromatic H3K9me3 regions as having the highest cell-to-cell variability in their methylation, likely reflecting inherent epigenetic instability of these late-replicating regions, compounded by differences in cell cycle stages among sampled cells.

Genome-wide programmable transcriptional memory by CRISPR-based epigenome editing.

A general approach for heritably altering gene expression has the potential to enable many discovery and therapeutic efforts. Here, we present CRISPRoff-a programmable epigenetic memory writer consisting of a single dead Cas9 fusion protein that establishes DNA methylation and repressive histone modifications. Transient CRISPRoff expression initiates highly specific DNA methylation and gene repression that is maintained through cell division and differentiation of stem cells to neurons. Pairing CRISPRoff with genome-wide screens and analysis of chromatin marks establishes rules for heritable gene silencing. We identify single guide RNAs (sgRNAs) capable of silencing the large majority of genes including those lacking canonical CpG islands (CGIs) and reveal a wide targeting window extending beyond annotated CGIs. The broad ability of CRISPRoff to initiate heritable gene silencing even outside of CGIs expands the canonical model of methylation-based silencing and enables diverse applications including genome-wide screens, multiplexed cell engineering, enhancer silencing, and mechanistic exploration of epigenetic inheritance.

Selective Modulation of a Pan-Essential Protein as a Therapeutic Strategy in Cancer.

Cancer dependency maps, which use CRISPR/Cas9 depletion screens to profile the landscape of genetic dependencies in hundreds of cancer cell lines, have identified context-specific dependencies that could be therapeutically exploited. An ideal therapy is both lethal and precise, but these depletion screens cannot readily distinguish between gene effects that are cytostatic or cytotoxic. Here, we use a diverse panel of functional genomic screening assays to identify NXT1 as a selective and rapidly lethal in vivo relevant genetic dependency in MYCN-amplified neuroblastoma. NXT1 heterodimerizes with NXF1, and together they form the principal mRNA nuclear export machinery. We describe a previously unrecognized mechanism of synthetic lethality between NXT1 and its paralog NXT2: their common essential binding partner NXF1 is lost only in the absence of both. We propose a potential therapeutic strategy for tumor-selective elimination of a protein that, if targeted directly, is expected to cause widespread toxicity. SIGNIFICANCE: We provide a framework for identifying new therapeutic targets from functional genomic screens. We nominate NXT1 as a selective lethal target in neuroblastoma and propose a therapeutic approach where the essential protein NXF1 can be selectively eliminated in tumor cells by exploiting the NXT1-NXT2 paralog relationship.See related commentary by Wang and Abdel-Wahab, p. 2129.This article is highlighted in the In This Issue feature, p. 2113.

Sarcoma classification by DNA methylation profiling.

Sarcomas are malignant soft tissue and bone tumours affecting adults, adolescents and children. They represent a morphologically heterogeneous class of tumours and some entities lack defining histopathological features. Therefore, the diagnosis of sarcomas is burdened with a high inter-observer variability and misclassification rate. Here, we demonstrate classification of soft tissue and bone tumours using a machine learning classifier algorithm based on array-generated DNA methylation data. This sarcoma classifier is trained using a dataset of 1077 methylation profiles from comprehensively pre-characterized cases comprising 62 tumour methylation classes constituting a broad range of soft tissue and bone sarcoma subtypes across the entire age spectrum. The performance is validated in a cohort of 428 sarcomatous tumours, of which 322 cases were classified by the sarcoma classifier. Our results demonstrate the potential of the DNA methylation-based sarcoma classification for research and future diagnostic applications.

GLI3 Is Associated With Neuronal Differentiation in SHH-Activated and WNT-Activated Medulloblastoma.

Glioma-associated oncogene homolog 3 (GLI3), whose main function is to inhibit GLI1, has been associated with neuronal differentiation in medulloblastoma. However, it is not clear what molecular subtype(s) show increased GLI3 expression. GLI3 levels were assessed by immunohistochemistry in 2 independent cohorts, including a total of 88 cases, and found to be high in both WNT- and SHH-activated medulloblastoma. Analysis of bulk mRNA expression data and single cell RNA sequencing studies confirmed that GLI1 and GLI3 are highly expressed in SHH-activated medulloblastoma, whereas GLI3 but not GLI1 is highly expressed in WNT-activated medulloblastoma. Immunohistochemical analysis has shown that GLI3 is expressed inside the neuronal differentiated nodules of SHH-activated medulloblastoma, whereas GLI1/2 are expressed in desmoplastic areas. In contrast, GLI3 is diffusely expressed in WNT-activated medulloblastoma, whereas GLI1 is suppressed. Our data suggest that GLI3 may be a master regulator of neuronal differentiation and morphology in these subgroups.

Single-Cell RNA-Seq Reveals Cellular Hierarchies and Impaired Developmental Trajectories in Pediatric Ependymoma.

Ependymoma is a heterogeneous entity of central nervous system tumors with well-established molecular groups. Here, we apply single-cell RNA sequencing to analyze ependymomas across molecular groups and anatomic locations to investigate their intratumoral heterogeneity and developmental origins. Ependymomas are composed of a cellular hierarchy initiating from undifferentiated populations, which undergo impaired differentiation toward three lineages of neuronal-glial fate specification. While prognostically favorable groups of ependymoma predominantly harbor differentiated cells, aggressive groups are enriched for undifferentiated cell populations. The delineated transcriptomic signatures correlate with patient survival and define molecular dependencies for targeted treatment approaches. Taken together, our analyses reveal a developmental hierarchy underlying ependymomas relevant to biological and clinical behavior.

Single-cell lineage analysis reveals genetic and epigenetic interplay in glioblastoma drug resistance.

BACKGROUND: Tumors can evolve and adapt to therapeutic pressure by acquiring genetic and epigenetic alterations that may be transient or stable. A precise understanding of how such events contribute to intratumoral heterogeneity, dynamic subpopulations, and overall tumor fitness will require experimental approaches to prospectively label, track, and characterize resistant or otherwise adaptive populations at the single-cell level. In glioblastoma, poor efficacy of receptor tyrosine kinase (RTK) therapies has been alternatively ascribed to genetic heterogeneity or to epigenetic transitions that circumvent signaling blockade. RESULTS: We combine cell lineage barcoding and single-cell transcriptomics to trace the emergence of drug resistance in stem-like glioblastoma cells treated with RTK inhibitors. Whereas a broad variety of barcoded lineages adopt a Notch-dependent persister phenotype that sustains them through early drug exposure, rare subclones acquire genetic changes that enable their rapid outgrowth over time. Single-cell analyses reveal that these genetic subclones gain copy number amplifications of the insulin receptor substrate-1 and substrate-2 (IRS1 or IRS2) loci, which activate insulin and AKT signaling programs. Persister-like cells and genomic amplifications of IRS2 and other loci are evident in primary glioblastomas and may underlie the inefficacy of targeted therapies in this disease. CONCLUSIONS: A method for combined lineage tracing and scRNA-seq reveals the interplay between complementary genetic and epigenetic mechanisms of resistance in a heterogeneous glioblastoma tumor model.