Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective.

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

Engineered MED12 mutations drive leiomyoma-like transcriptional and metabolic programs by altering the 3D genome compartmentalization.

Nearly 70% of Uterine fibroid (UF) tumors are driven by recurrent MED12 hotspot mutations. Unfortunately, no cellular models could be generated because the mutant cells have lower fitness in 2D culture conditions. To address this, we employ CRISPR to precisely engineer MED12 Gly44 mutations in UF-relevant myometrial smooth muscle cells. The engineered mutant cells recapitulate several UF-like cellular, transcriptional, and metabolic alterations, including altered Tryptophan/kynurenine metabolism. The aberrant gene expression program in the mutant cells is, in part, driven by a substantial 3D genome compartmentalization switch. At the cellular level, the mutant cells gain enhanced proliferation rates in 3D spheres and form larger lesions in vivo with elevated production of collagen and extracellular matrix deposition. These findings indicate that the engineered cellular model faithfully models key features of UF tumors and provides a platform for the broader scientific community to characterize genomics of recurrent MED12 mutations.

Enhancer Coamplification and Hijacking Promote Oncogene Expression in Liposarcoma.

SIGNIFICANCE: Comprehensive profiling of the enhancer landscape and 3D genome structure in liposarcoma identifies extensive enhancer-oncogene coamplification and enhancer hijacking events, deepening the understanding of how oncogenes are regulated in cancer.

Engineered MED12 mutations drive uterine fibroid-like transcriptional and metabolic programs by altering the 3D genome compartmentalization.

Uterine fibroid (UF) tumors originate from a mutated smooth muscle cell (SMC). Nearly 70% of these tumors are driven by hotspot recurrent somatic mutations in the MED12 gene; however, there are no tractable genetic models to study the biology of UF tumors because, under culture conditions, the non-mutant fibroblasts outgrow the mutant SMC cells, resulting in the conversion of the population to WT phenotype. The lack of faithful cellular models hampered our ability to delineate the molecular pathways downstream of MED12 mutations and identify therapeutics that may selectively target the mutant cells. To overcome this challenge, we employed CRISPR knock-in with a sensitive PCR-based screening strategy to precisely engineer cells with mutant MED12 Gly44, which constitutes 50% of MED12 exon two mutations. Critically, the engineered myometrial SMC cells recapitulate several UF-like cellular, transcriptional and metabolic alterations, including enhanced proliferation rates in 3D spheres and altered Tryptophan/kynurenine metabolism. Our transcriptomic analysis supported by DNA synthesis tracking reveals that MED12 mutant cells, like UF tumors, have heightened expression of DNA repair genes but reduced DNA synthesis rates. Consequently, these cells accumulate significantly higher rates of DNA damage and are selectively more sensitive to common DNA-damaging chemotherapy, indicating mutation-specific and therapeutically relevant vulnerabilities. Our high-resolution 3D chromatin interaction analysis demonstrates that the engineered MED12 mutations drive aberrant genomic activity due to a genome-wide chromatin compartmentalization switch. These findings indicate that the engineered cellular model faithfully models key features of UF tumors and provides a novel platform for the broader scientific community to characterize genomics of recurrent MED12 mutations and discover potential therapeutic targets.

Breakage fusion bridge cycles drive high oncogene copy number, but not intratumoral genetic heterogeneity or rapid cancer genome change.

Oncogene amplification is a major driver of cancer pathogenesis. Breakage fusion bridge (BFB) cycles, like extrachromosomal DNA (ecDNA), can lead to high copy numbers of oncogenes, but their impact on intratumoral heterogeneity, treatment response, and patient survival are not well understood due to difficulty in detecting them by DNA sequencing. We describe a novel algorithm that detects and reconstructs BFB amplifications using optical genome maps (OGMs), called OM2BFB. OM2BFB showed high precision (>93%) and recall (92%) in detecting BFB amplifications in cancer cell lines, PDX models and primary tumors. OM-based comparisons demonstrated that short-read BFB detection using our AmpliconSuite (AS) toolkit also achieved high precision, albeit with reduced sensitivity. We detected 371 BFB events using whole genome sequences from 2,557 primary tumors and cancer lines. BFB amplifications were preferentially found in cervical, head and neck, lung, and esophageal cancers, but rarely in brain cancers. BFB amplified genes show lower variance of gene expression, with fewer options for regulatory rewiring relative to ecDNA amplified genes. BFB positive (BFB (+)) tumors showed reduced heterogeneity of amplicon structures, and delayed onset of resistance, relative to ecDNA(+) tumors. EcDNA and BFB amplifications represent contrasting mechanisms to increase the copy numbers of oncogene with markedly different characteristics that suggest different routes for intervention.

Subtype-specific 3D genome alteration in acute myeloid leukaemia.

Acute myeloid leukaemia (AML) represents a set of heterogeneous myeloid malignancies, and hallmarks include mutations in epigenetic modifiers, transcription factors and kinases1-5. The extent to which mutations in AML drive alterations in chromatin 3D structure and contribute to myeloid transformation is unclear. Here we use Hi-C and whole-genome sequencing to analyse 25 samples from patients with AML and 7 samples from healthy donors. Recurrent and subtype-specific alterations in A/B compartments, topologically associating domains and chromatin loops were identified. RNA sequencing, ATAC with sequencing and CUT&Tag for CTCF, H3K27ac and H3K27me3 in the same AML samples also revealed extensive and recurrent AML-specific promoter-enhancer and promoter-silencer loops. We validated the role of repressive loops on their target genes by CRISPR deletion and interference. Structural variation-induced enhancer-hijacking and silencer-hijacking events were further identified in AML samples. Hijacked enhancers play a part in AML cell growth, as demonstrated by CRISPR screening, whereas hijacked silencers have a downregulating role, as evidenced by CRISPR-interference-mediated de-repression. Finally, whole-genome bisulfite sequencing of 20 AML and normal samples revealed the delicate relationship between DNA methylation, CTCF binding and 3D genome structure. Treatment of AML cells with a DNA hypomethylating agent and triple knockdown of DNMT1, DNMT3A and DNMT3B enabled the manipulation of DNA methylation to revert 3D genome organization and gene expression. Overall, this study provides a resource for leukaemia studies and highlights the role of repressive loops and hijacked cis elements in human diseases.

Subtype-associated epigenomic landscape and 3D genome structure in bladder cancer.

Muscle-invasive bladder cancers are characterized by their distinct expression of luminal and basal genes, which could be used to predict key clinical features such as disease progression and overall survival. Transcriptionally, FOXA1, GATA3, and PPARG are shown to be essential for luminal subtype-specific gene regulation and subtype switching, while TP63, STAT3, and TFAP2 family members are critical for regulation of basal subtype-specific genes. Despite these advances, the underlying epigenetic mechanisms and 3D chromatin architecture responsible for subtype-specific regulation in bladder cancer remain unknown. RESULT: We determine the genome-wide transcriptome, enhancer landscape, and transcription factor binding profiles of FOXA1 and GATA3 in luminal and basal subtypes of bladder cancer. Furthermore, we report the first-ever mapping of genome-wide chromatin interactions by Hi-C in both bladder cancer cell lines and primary patient tumors. We show that subtype-specific transcription is accompanied by specific open chromatin and epigenomic marks, at least partially driven by distinct transcription factor binding at distal enhancers of luminal and basal bladder cancers. Finally, we identify a novel clinically relevant transcription factor, Neuronal PAS Domain Protein 2 (NPAS2), in luminal bladder cancers that regulates other subtype-specific genes and influences cancer cell proliferation and migration. CONCLUSION: In summary, our work identifies unique epigenomic signatures and 3D genome structures in luminal and basal urinary bladder cancers and suggests a novel link between the circadian transcription factor NPAS2 and a clinical bladder cancer subtype.

Genome-wide colocalization of RNA-DNA interactions and fusion RNA pairs.

Fusion transcripts are used as biomarkers in companion diagnoses. Although more than 15,000 fusion RNAs have been identified from diverse cancer types, few common features have been reported. Here, we compared 16,410 fusion transcripts detected in cancer (from a published cohort of 9,966 tumor samples of 33 cancer types) with genome-wide RNA-DNA interactions mapped in two normal, noncancerous cell types [using iMARGI, an enhanced version of the mapping of RNA-genome interactions (MARGI) assay]. Among the top 10 most significant RNA-DNA interactions in normal cells, 5 colocalized with the gene pairs that formed fusion RNAs in cancer. Furthermore, throughout the genome, the frequency of a gene pair to exhibit RNA-DNA interactions is positively correlated with the probability of this gene pair to present documented fusion transcripts in cancer. To test whether RNA-DNA interactions in normal cells are predictive of fusion RNAs, we analyzed these in a validation cohort of 96 lung cancer samples using RNA sequencing (RNA-seq). Thirty-seven of 42 fusion transcripts in the validation cohort were found to exhibit RNA-DNA interactions in normal cells. Finally, by combining RNA-seq, single-molecule RNA FISH, and DNA FISH, we detected a cancer sample with EML4-ALK fusion RNA without forming the EML4-ALK fusion gene. Collectively, these data suggest an RNA-poise model, where spatial proximity of RNA and DNA could poise for the creation of fusion transcripts.