Structural variants shape driver combinations and outcomes in pediatric high-grade glioma

Effective targeting of PDGFRA-altered high-grade glioma with avapritinib

PDGFRA is crucial to tumorigenesis and frequently genomically altered in high-grade glioma (HGG). In a comprehensive dataset of pediatric HGG (n = 261), we detect PDGFRA mutations and/or amplifications in 15% of cases, suggesting PDGFRA as a therapeutic target. We reveal that the PDGFRA/KIT inhibitor avapritinib shows (1) selectivity for PDGFRA inhibition, (2) distinct patterns of subcellular effects, (3) in vitro and in vivo activity in patient-derived HGG models, and (4) effective blood-brain barrier penetration in mice and humans. Furthermore, we report preliminary clinical real-world experience using avapritinib in pediatric and young adult patients with predominantly recurrent/refractory PDGFRA-altered HGG (n = 8). Our early data demonstrate that avapritinib is well tolerated and results in radiographic response in 3/7 cases, suggesting a potential role for avapritinib in the treatment of HGG with specific PDGFRA alterations. Overall, these translational results underscore the therapeutic potential of PDGFRA inhibition with avapritinib in HGG.

Amplified dosage of the NKX2-1 lineage transcription factor controls its oncogenic role in lung adenocarcinoma

Amplification-mediated oncogene overexpression is a critical and widespread driver event in cancer, yet our understanding of how amplification and dosage mediate oncogene regulation is limited. Here, we find that the most significant focal amplification event in lung adenocarcinoma (LUAD) targets a lineage "super-enhancer" near the NKX2-1 lineage transcription factor. The NKX2-1 super-enhancer is targeted by focal and co-amplification with NKX2-1 and controls NKX2-1 expression and regulation. We find that NKX2-1 directly controls enhancer accessibility to drive a lineage-addicted state in LUAD. We precisely map the effects of NKX2-1 dosage modulation upon both overexpression and knockdown and identify both linear and non-linear regulation by NKX2-1 dosage. We find that NKX2-1 is a widespread dependency in LUAD cell lines and that NKX2-1 confers persistence to EGFR inhibitors. Our data suggest a defining role for dosage in the oncogenic regulation of amplified NKX2-1 and that amplified NKX2-1 lineage addiction defines LUAD tumors.

Machine learning-predicted chromatin organization landscape across pediatric tumors

Structural variants (SVs) are increasingly recognized as important contributors to oncogenesis through their effects on 3D genome folding. Recent advances in whole-genome sequencing have enabled large-scale profiling of SVs across diverse tumors, yet experimental characterization of their individual impact on genome folding remains infeasible. Here, we leveraged a convolutional neural network, Akita, to predict disruptions in genome folding caused by somatic SVs identified in 61 tumor types from the Children's Brain Tumor Network dataset. Our analysis reveals significant variability in SV-induced disruptions across tumor types, with the most disruptive SVs coming from lymphomas and sarcomas, metastatic tumors, and germline cell tumors. Dimensionality reduction of disruption scores identified five recurrently disrupted regions enriched for high-impact SVs across multiple tumors. Some of these regions are highly disrupted despite not being highly mutated, and harbor tumor-associated genes and transcriptional regulators. To further interpret the functional relevance of high-scoring SVs, we integrated epigenetic data and developed a modified Activity-by-Contact scoring approach to prioritize SVs with disrupted genome contacts at active enhancers. This method highlighted highly disruptive SVs near key oncogenes, as well as novel candidate loci potentially implicated in tumorigenesis. These findings highlight the utility of machine learning for identifying novel SVs, loci, and genetic mechanisms contributing to pediatric cancers. This framework provides a foundation for future studies linking SV-driven regulatory changes to cancer pathogenesis.

A sequence context-based approach for classifying tumor structural variants without paired normal samples

Although several recent studies have characterized structural variants (SVs) in germline and cancer genomes independently, the genomic contexts of these SVs have not been comprehensively compared. We examined similarities and differences between 2 million germline and 115 thousand tumor SVs from a cohort of 963 patients from The Cancer Genome Atlas. We found significant differences in features related to their genomic sequences and localization that suggest differences between SV-generating processes and selective pressures. For example, our results show that features linked to transposon-mediated processes are associated with germline SVs, while somatic SVs more frequently show features characteristic of chromoanagenesis. These genomic differences enabled us to develop a classifier-the Germline and Tumor Structural Variant or "the great GaTSV" -that accurately distinguishes between germline and cancer SVs in tumor samples that lack a matched normal sample.

(B)On(e)-cohistones and the epigenetic alterations at the root of bone cancer

Identification of mutations in histones in a number of human neoplasms and developmental syndromes represents the most compelling evidence to date for a causal role of epigenetic perturbations in human disease. In most cases, these mutations have gain of function properties that cause deviation from normal developmental processes leading to embryo defects and/or neoplastic transformation. These exciting discoveries represent a step-change in our understanding of the role of chromatin (dys)regulation in development and disease. However, the mechanisms of action of oncogenic histone mutations (oncohistones) remain only partially understood. Here, we critically assess existing literature on oncohistones focussing mainly on bone neoplasms. We show how it is possible to draw parallels with some of the cell-autonomous mechanisms of action described in paediatric brain cancer, although the functions of oncohistones in bone tumours remain under-investigated. In this respect, it is becoming clear that histone mutations targeting the same residues display, at least in part, tissue-specific oncogenic mechanisms. Furthermore, it is emerging that cancer cells carrying oncohistones can modify the surrounding microenvironment to support growth and/or alter differentiation trajectories. A better understanding of oncohistone function in different neoplasms provide potential for identification of signalling that could be targeted therapeutically. Finally, we discuss some of the main concepts and future directions in this research area, while also drawing possible connections and parallels with other cancer epigenetic mechanisms.

Dual targeting of histone deacetylases and MYC as potential treatment strategy for H3-K27M pediatric gliomas

Diffuse midline gliomas (DMGs) are aggressive and fatal pediatric tumors of the central nervous system that are highly resistant to treatments. Lysine to methionine substitution of residue 27 on histone H3 (H3-K27M) is a driver mutation in DMGs, reshaping the epigenetic landscape of these cells to promote tumorigenesis. H3-K27M gliomas are characterized by deregulation of histone acetylation and methylation pathways, as well as the oncogenic MYC pathway. In search of effective treatment, we examined the therapeutic potential of dual targeting of histone deacetylases (HDACs) and MYC in these tumors. Treatment of H3-K27M patient-derived cells with Sulfopin, an inhibitor shown to block MYC-driven tumors in vivo, in combination with the HDAC inhibitor Vorinostat, resulted in substantial decrease in cell viability. Moreover, transcriptome and epigenome profiling revealed synergistic effect of this drug combination in downregulation of prominent oncogenic pathways such as mTOR. Finally, in vivo studies of patient-derived orthotopic xenograft models showed significant tumor growth reduction in mice treated with the drug combination. These results highlight the combined treatment with PIN1 and HDAC inhibitors as a promising therapeutic approach for these aggressive tumors.

New facets in the chromatin-based regulation of genome maintenance

The maintenance of genome integrity by DNA damage response machineries is key to protect cells against pathological development. In cell nuclei, these genome maintenance machineries operate in the context of chromatin, where the DNA wraps around histone proteins. Here, we review recent findings illustrating how the chromatin substrate modulates genome maintenance mechanisms, focusing on the regulatory role of histone variants and post-translational modifications. In particular, we discuss how the pre-existing chromatin landscape impacts DNA damage formation and guides DNA repair pathway choice, and how DNA damage-induced chromatin alterations control DNA damage signaling and repair, and DNA damage segregation through cell divisions. We also highlight that pathological alterations of histone proteins may trigger genome instability by impairing chromosome segregation and DNA repair, thus defining new oncogenic mechanisms and opening up therapeutic options.

3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression

Cancer genomes are composed of many complex structural alterations on chromosomes and extrachromosomal DNA (ecDNA), making it difficult to identify non-coding enhancer regions that are hijacked to activate oncogene expression. Here, we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking. HAPI analysis of HiChIP data from 34 cancer cell lines identified enhancer hijacking events that activate both known and potentially novel oncogenes such as MYC, CCND1, ETV1, CRKL, and ID4. Furthermore, we found enhancer hijacking among multiple oncogenes from different chromosomes, often including MYC, on the same complex amplicons such as ecDNA. We characterized a MYC-ERBB2 chimeric ecDNA, in which ERBB2 heavily hijacks MYC's enhancers. Notably, CRISPRi of the MYC promoter led to increased interaction of ERBB2 with MYC enhancers and elevated ERBB2 expression. Our HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation.

3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression

Cancer genomes are composed of many complex structural alterations on chromosomes and extrachromosomal DNA (ecDNA), making it difficult to identify non-coding enhancer regions that are hijacked to activate oncogene expression. Here, we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking. HAPI analysis of HiChIP data from 34 cancer cell lines identified enhancer hijacking events that activate both known and potentially novel oncogenes such as MYC, CCND1 , ETV1 , CRKL , and ID4 . Furthermore, we found enhancer hijacking among multiple oncogenes from different chromosomes, often including MYC , on the same complex amplicons such as ecDNA. We characterized a MYC - ERBB2 chimeric ecDNA, in which ERBB2 heavily hijacks MYC 's enhancers. Notably, CRISPRi of the MYC promoter led to increased interaction of ERBB2 with MYC enhancers and elevated ERBB2 expression. Our HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation.

Comparison of germline and somatic structural variants in cancers reveal systematic differences in variant generating and selection processes

Although several recent studies have characterized structural variants (SVs) in germline and cancer genomes, the features of SVs in these different contexts have not been directly compared. We examined similarities and differences between 2 million germline and 115 thousand tumor SVs from a cohort of 963 patients from The Cancer Genome Atlas (TCGA). We found significant differences in features related to their genomic sequences and localization that suggest differences between SV-generating processes and selective pressures. For example, we found that transposon-mediated processes shape germline much more than somatic SVs, while somatic SVs more frequently show features characteristic of chromoanagenesis. These differences were extensive enough to enable us to develop a classifier - "the great GaTSV" - that accurately distinguishes between germline and cancer SVs in tumor samples that lack a matched normal sample.

Aberrant DNA repair reveals a vulnerability in histone H3.3-mutant brain tumors

Pediatric high-grade gliomas (pHGG) are devastating and incurable brain tumors with recurrent mutations in histone H3.3. These mutations promote oncogenesis by dysregulating gene expression through alterations of histone modifications. We identify aberrant DNA repair as an independent mechanism, which fosters genome instability in H3.3 mutant pHGG, and opens new therapeutic options. The two most frequent H3.3 mutations in pHGG, K27M and G34R, drive aberrant repair of replication-associated damage by non-homologous end joining (NHEJ). Aberrant NHEJ is mediated by the DNA repair enzyme polynucleotide kinase 3'-phosphatase (PNKP), which shows increased association with mutant H3.3 at damaged replication forks. PNKP sustains the proliferation of cells bearing H3.3 mutations, thus conferring a molecular vulnerability, specific to mutant cells, with potential for therapeutic targeting.

Extrachromosomal circular MiR-17-92 amplicon promotes HCC

BACKGROUND AND AIMS

Extrachromosomal circular DNAs (eccDNAs) are prevalent in cancer genomes and emerge as a class of crucial yet less characterized oncogenic drivers. However, the structure, composition, genome-wide frequency, and contribution of eccDNAs in HCC, one of the most fatal and prevalent cancers, remain unexplored. In this study, we provide a comprehensive characterization of eccDNAs in human HCC and demonstrate an oncogenic role of microRNA (miRNA)-17-92-containing eccDNAs in tumor progression.

APPROACH AND RESULTS

Using the circle-sequencing method, we identify and characterize more than 230,000 eccDNAs from 4 paired samples of HCC tumor and adjacent nontumor liver tissues. EccDNAs are highly enriched in HCC tumors, preferentially originate from certain chromosomal hotspots, and are correlated with differential gene expression. Particularly, a series of eccDNAs carrying the miRNA-17-92 cluster are validated by outward PCR and Sanger sequencing. Quantitative PCR analyses reveal that miRNA-17-92-containing eccDNAs, along with the expression of their corresponding miRNAs, are elevated in HCC tumors and associated with poor outcomes and the age of HCC patients. More intriguingly, exogenous expression of artificial DNA circles harboring the miR-17-92 cluster, which is synthesized by the ligase-assisted minicircle accumulation method, can significantly accelerate HCC cell proliferation and migration.

CONCLUSIONS

These findings delineate the genome-wide eccDNAs profiling of HCC and highlight the functional significance of miRNA-containing eccDNAs in tumorigenesis, providing insight into HCC pathogenesis and cancer therapy, as well as eccDNA and miRNA biology.

SMARCB1 loss activates patient-specific distal oncogenic enhancers in malignant rhabdoid tumors

Malignant rhabdoid tumor (MRT) is a highly malignant and often lethal childhood cancer. MRTs are genetically defined by bi-allelic inactivating mutations in SMARCB1, a member of the BRG1/BRM-associated factors (BAF) chromatin remodeling complex. Mutations in BAF complex members are common in human cancer, yet their contribution to tumorigenesis remains in many cases poorly understood. Here, we study derailed regulatory landscapes as a consequence of SMARCB1 loss in the context of MRT. Our multi-omics approach on patient-derived MRT organoids reveals a dramatic reshaping of the regulatory landscape upon SMARCB1 reconstitution. Chromosome conformation capture experiments subsequently reveal patient-specific looping of distal enhancer regions with the promoter of the MYC oncogene. This intertumoral heterogeneity in MYC enhancer utilization is also present in patient MRT tissues as shown by combined single-cell RNA-seq and ATAC-seq. We show that loss of SMARCB1 activates patient-specific epigenetic reprogramming underlying MRT tumorigenesis.

Ppm1d truncating mutations promote the development of genotoxic stress-induced AML

Hematopoietic stem cells (HSCs) ensure blood cell production during the life-time of an organism, and to do so they need to balance self-renewal, proliferation, differentiation, and migration in a steady state as well as in response to stress or injury. Importantly, aberrant proliferation of HSCs leads to hematological malignancies, and thus, tight regulation by various tumor suppressor pathways, including p53, is essential. Protein phosphatase magnesium-dependent 1 delta (PPM1D) is a negative regulator of p53 and promotes cell survival upon induction of genotoxic stress. Truncating mutations in the last exon of PPM1D lead to the production of a stable, enzymatically active protein and are commonly associated with clonal hematopoiesis. Using a transgenic mouse model, we demonstrate that truncated PPM1D reduces self-renewal of HSCs in basal conditions but promotes the development of aggressive AML after exposure to ionizing radiation. Inhibition of PPM1D suppressed the colony growth of leukemic stem and progenitor cells carrying the truncated PPM1D, and remarkably, it provided protection against irradiation-induced cell growth. Altogether, we demonstrate that truncated PPM1D affects HSC maintenance, disrupts normal hematopoiesis, and that its inhibition could be beneficial in the context of therapy-induced AML.

Somatic structural variation signatures in pediatric brain tumors

Brain cancer is the leading cause of cancer-related death in children. Somatic structural variations (SVs), large-scale alterations in DNA, remain poorly understood in pediatric brain tumors. Here, we detect a total of 13,199 high-confidence somatic SVs in 744 whole-genome sequences of pediatric brain tumors from the Pediatric Brain Tumor Atlas. The somatic SV occurrences have tremendous diversity among the cohort and across different tumor types. We decompose mutational signatures of clustered complex SVs, non-clustered complex SVs, and simple SVs separately to infer their mutational mechanisms. Our finding of many tumor types carrying unique sets of SV signatures suggests that distinct molecular mechanisms shape genome instability in different tumor types. The patterns of somatic SV signatures in pediatric brain tumors are substantially different from those in adult cancers. The convergence of multiple SV signatures on several major cancer driver genes implies vital roles of somatic SVs in disease progression.

Dosage amplification dictates oncogenic regulation by the NKX2-1 lineage factor in lung adenocarcinoma

Amplified oncogene expression is a critical and widespread driver event in cancer, yet our understanding of how amplification-mediated elevated dosage mediates oncogenic regulation is limited. Here, we find that the most significant focal amplification event in lung adenocarcinoma (LUAD) targets a lineage super-enhancer near the NKX2-1 lineage transcription factor. The NKX2-1 super-enhancer is targeted by focal and co-amplification with NKX2-1, and activation or repression controls NKX2-1 expression. We find that NKX2-1 is a widespread dependency in LUAD cell lines, where NKX2-1 pioneers enhancer accessibility to drive a lineage addicted state in LUAD, and NKX2-1 confers persistence to EGFR inhibitors. Notably, we find that oncogenic NKX2-1 regulation requires expression above a minimum dosage threshold-NKX2-1 dosage below this threshold is insufficient for cell viability, enhancer remodeling, and TKI persistence. Our data suggest that copy-number amplification can be a gain-of-function alteration, wherein amplification elevates oncogene expression above a critical dosage required for oncogenic regulation and cancer cell survival.

The 3D chromatin landscape of rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a pediatric soft tissue cancer with a lack of precision therapy options for patients. We hypothesized that with a general paucity of known mutations in RMS, chromatin structural driving mechanisms are essential for tumor proliferation. Thus, we carried out high-depth in situ Hi-C in representative cell lines and patient-derived xenografts (PDXs) to define chromatin architecture in each major RMS subtype. We report a comprehensive 3D chromatin structural analysis and characterization of fusion-positive (FP-RMS) and fusion-negative RMS (FN-RMS). We have generated spike-in in situ Hi-C chromatin interaction maps for the most common FP-RMS and FN-RMS cell lines and compared our data with PDX models. In our studies, we uncover common and distinct structural elements in large Mb-scale chromatin compartments, tumor-essential genes within variable topologically associating domains and unique patterns of structural variation. Our high-depth chromatin interactivity maps and comprehensive analyses provide context for gene regulatory events and reveal functional chromatin domains in RMS.

K27M in canonical and noncanonical H3 variants occurs in distinct oligodendroglial cell lineages in brain midline gliomas

Canonical (H3.1/H3.2) and noncanonical (H3.3) histone 3 K27M-mutant gliomas have unique spatiotemporal distributions, partner alterations and molecular profiles. The contribution of the cell of origin to these differences has been challenging to uncouple from the oncogenic reprogramming induced by the mutation. Here, we perform an integrated analysis of 116 tumors, including single-cell transcriptome and chromatin accessibility, 3D chromatin architecture and epigenomic profiles, and show that K27M-mutant gliomas faithfully maintain chromatin configuration at developmental genes consistent with anatomically distinct oligodendrocyte precursor cells (OPCs). H3.3K27M thalamic gliomas map to prosomere 2-derived lineages. In turn, H3.1K27M ACVR1-mutant pontine gliomas uniformly mirror early ventral NKX6-1⁺/SHH-dependent brainstem OPCs, whereas H3.3K27M gliomas frequently resemble dorsal PAX3⁺/BMP-dependent progenitors. Our data suggest a context-specific vulnerability in H3.1K27M-mutant SHH-dependent ventral OPCs, which rely on acquisition of ACVR1 mutations to drive aberrant BMP signaling required for oncogenesis. The unifying action of K27M mutations is to restrict H3K27me3 at PRC2 landing sites, whereas other epigenetic changes are mainly contingent on the cell of origin chromatin state and cycling rate.

The landscape of tumor cell states and spatial organization in H3-K27M mutant diffuse midline glioma across age and location

Histone 3 lysine27-to-methionine (H3-K27M) mutations most frequently occur in diffuse midline gliomas (DMGs) of the childhood pons but are also increasingly recognized in adults. Their potential heterogeneity at different ages and midline locations is vastly understudied. Here, through dissecting the single-cell transcriptomic, epigenomic and spatial architectures of a comprehensive cohort of patient H3-K27M DMGs, we delineate how age and anatomical location shape glioma cell-intrinsic and -extrinsic features in light of the shared driver mutation. We show that stem-like oligodendroglial precursor-like cells, present across all clinico-anatomical groups, display varying levels of maturation dependent on location. We reveal a previously underappreciated relationship between mesenchymal cancer cell states and age, linked to age-dependent differences in the immune microenvironment. Further, we resolve the spatial organization of H3-K27M DMG cell populations and identify a mitotic oligodendroglial-lineage niche. Collectively, our study provides a powerful framework for rational modeling and therapeutic interventions.

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.