Search

Extrachromosomal DNA (ecDNA) is a common mode of oncogene amplification but is challenging to analyze. Here, we adapt CRISPR-CATCH, in vitro CRISPR-Cas9 treatment and pulsed field gel electrophoresis of agarose-entrapped genomic DNA, previously developed for bacterial chromosome segments, to isolate megabase-sized human ecDNAs. We demonstrate strong enrichment of ecDNA molecules containing EGFR, FGFR2 and MYC from human cancer cells and NRAS ecDNA from human metastatic melanoma with acquired therapeutic resistance. Targeted enrichment of ecDNA versus chromosomal DNA enabled phasing of genetic variants, identified the presence of an EGFRvIII mutation exclusively on ecDNAs and supported an excision model of ecDNA genesis in a glioblastoma model. CRISPR-CATCH followed by nanopore sequencing enabled single-molecule ecDNA methylation profiling and revealed hypomethylation of the EGFR promoter on ecDNAs. We distinguished heterogeneous ecDNA species within the same sample by size and sequence with base-pair resolution and discovered functionally specialized ecDNAs that amplify select enhancers or oncogene-coding sequences.

Focal amplifications (FA) can mediate targeted therapy resistance in cancer. Understanding the structure and dynamics of FAs is critical for designing treatments that overcome plasticity-mediated resistance. We developed a melanoma model of dual MAPK inhibitor (MAPKi) resistance that bears BRAFV600 amplifications through either extrachromosomal DNA (ecDNA)/double minutes (DM) or intrachromosomal homogenously staining regions (HSR). Cells harboring BRAFV600E FAs displayed mode switching between DMs and HSRs, from both de novo genetic changes and selection of preexisting subpopulations. Plasticity is not exclusive to ecDNAs, as cells harboring HSRs exhibit drug addiction–driven structural loss of BRAF amplicons upon dose reduction. FA mechanisms can couple with kinase domain duplications and alternative splicing to enhance resistance. Drug-responsive amplicon plasticity is observed in the clinic and can involve other MAPK pathway genes, such as RAF1 and NRAS. BRAF FA-mediated dual MAPKi–resistant cells are more sensitive to proferroptotic drugs, extending the spectrum of ferroptosis sensitivity in MAPKi resistance beyond cases of dedifferentiation.Significance:Understanding the structure and dynamics of oncogene amplifications is critical for overcoming tumor relapse. BRAF amplifications are highly plastic under MAPKi dosage challenges in melanoma, through involvement of de novo genomic alterations, even in the HSR mode. Moreover, BRAF FA-driven, dual MAPKi–resistant cells extend the spectrum of resistance-linked ferroptosis sensitivity.This article is highlighted in the In This Issue feature, p. 873

Supplementary Table from Plasticity of Extrachromosomal and Intrachromosomal BRAF Amplifications in Overcoming Targeted Therapy Dosage Challenges

Supplementary Figure from Plasticity of Extrachromosomal and Intrachromosomal BRAF Amplifications in Overcoming Targeted Therapy Dosage Challenges

Focal amplifications (FA) can mediate targeted therapy resistance in cancer. Understanding the structure and dynamics of FAs is critical for designing treatments that overcome plasticity-mediated resistance. We developed a melanoma model of dual MAPK inhibitor (MAPKi) resistance that bears BRAFV600 amplifications through either extrachromosomal DNA (ecDNA)/double minutes (DM) or intrachromosomal homogenously staining regions (HSR). Cells harboring BRAFV600E FAs displayed mode switching between DMs and HSRs, from both de novo genetic changes and selection of preexisting subpopulations. Plasticity is not exclusive to ecDNAs, as cells harboring HSRs exhibit drug addiction–driven structural loss of BRAF amplicons upon dose reduction. FA mechanisms can couple with kinase domain duplications and alternative splicing to enhance resistance. Drug-responsive amplicon plasticity is observed in the clinic and can involve other MAPK pathway genes, such as RAF1 and NRAS. BRAF FA-mediated dual MAPKi–resistant cells are more sensitive to proferroptotic drugs, extending the spectrum of ferroptosis sensitivity in MAPKi resistance beyond cases of dedifferentiation.Significance:Understanding the structure and dynamics of oncogene amplifications is critical for overcoming tumor relapse. BRAF amplifications are highly plastic under MAPKi dosage challenges in melanoma, through involvement of de novo genomic alterations, even in the HSR mode. Moreover, BRAF FA-driven, dual MAPKi–resistant cells extend the spectrum of resistance-linked ferroptosis sensitivity.This article is highlighted in the In This Issue feature, p. 873

SUMMARYExtrachromosomal circular DNA (ecDNA) is an important driver of aggressive tumor growth, promoting high oncogene copy number, intratumoral heterogeneity, accelerated evolution of drug resistance, enhancer rewiring, and poor outcome. ecDNA has been reported in medulloblastoma (MB), the most common malignant pediatric brain tumor, but the ecDNA landscape and its association with specific MB subgroups, its impact on enhancer rewiring, and its potential clinical implications, are not known. We assembled a retrospective cohort of 468 MB patient samples with available whole genome sequencing (WGS) data covering the four major MB subgroups WNT, SHH, Group 3 and Group 4. Using computational methods for the detection and reconstruction of ecDNA1, we find ecDNA in 82 patients (18%) and observe that ecDNA+ MB patients are more than twice as likely to relapse and three times as likely to die of disease. In addition, we find that individual medulloblastoma tumors often harbor multiple ecDNAs, each containing different amplified oncogenes along with co-amplified non-coding regulatory enhancers. ecDNA was substantially more prevalent among 31 analyzed patient-derived xenograft (PDX) models and cell lines than in our patient cohort. By mapping the accessible chromatin and 3D conformation landscapes of MB tumors that harbor ecDNA, we observe frequent candidate “enhancer rewiring” events that spatially link oncogenes with co-amplified enhancers. Our study reveals the frequency and diversity of ecDNA in a subset of highly aggressive tumors and suggests enhancer rewiring as a frequent oncogenic mechanism of ecDNAs in MB. Further, these results demonstrate that ecDNA is a frequent and potent driver of poor outcome in MB patients.

Summary Optical mapping (OM) provides single-molecule readouts of fluorescently labeled sequence motifs on long fragments of DNA, resolved to nucleotide-level coordinates. With the advent of microfluidic technologies for analysis of DNA molecules, it is possible to inexpensively generate long OM data (>150 kbp) at high coverage. In addition to scaffolding for de novo assembly, OM data can be aligned to a reference genome for identification of genomic structural variants. We introduce FaNDOM (Fast Nested Distance Seeding of Optical Maps)—an optical map alignment tool that greatly reduces the search space of the alignment process. On four benchmark human datasets, FaNDOM was significantly (4–14×) faster than competing tools while maintaining comparable sensitivity and specificity. We used FaNDOM to map variants in three cancer cell lines and identified many biologically interesting structural variants, including deletions, duplications, gene fusions and gene-disrupting rearrangements. FaNDOM is publicly available at https://github.com/jluebeck/FaNDOM., Highlights • FaNDOM is a fast open-source aligner for OM data • It utilizes a novel filtering strategy to reduce the search space of alignment • The method enables discovery of large, complex genomic structural variants • Structural variants suggested by FaNDOM include gene fusions and gene disruptions, The bigger picture Optical mapping (OM) is a rapidly maturing strategy for detecting large-scale rearrangements in genomes, leveraging ultra-long fragments of DNA imaged at very high depth of coverage (>100×). OM data reflect an orthogonal strategy to DNA sequencing, instead utilizing image-based detection of fluorescent tags associated with specific DNA motifs. The resulting data can be aligned back to the reference genome for discovery of genomic rearrangements and karyotypic abnormalities. Existing methods, however, are computationally demanding, making discovery harder. We present a novel method, FaNDOM, for alignment of OM data to the reference genome, and the additional discovery of structural variants. FaNDOM utilizes fast filtering algorithms based on constructing graph-based chains of seed matches, achieving orders of magnitude speedup, while maintaining high sensitivity, enabling a more comprehensive search of complex structural variations involving hundreds of kbp., Optical mapping data is an orthogonal technique to DNA sequencing for the identification of genomic structural variants (SVs). We present a method, FaNDOM, which performs fast alignment of optical mapping data to the reference genome for identification of SVs. FaNDOM utilizes a novel filtering strategy, vastly reducing the search space of the alignment process, enabling rapid discovery of biologically interesting events.

BackgroundThe varied clinical manifestations and outcomes in patients with SARS-CoV-2 infections implicate a role of host-genetics in the predisposition to disease severity. This is supported by evidence that is now emerging, where initial reports identify common risk factors and rare genetic variants associated with high risk for severe/ life-threatening COVID-19. Impressive global efforts have focused on either identifying common genetic factors utilizing short-read sequencing data in Genome-Wide Association Studies (GWAS) or whole-exome and genome studies to interrogate the human genome at the level of detecting single nucleotide variants (SNVs) and short indels. However, these studies lack the sensitivity to accurately detect several classes of variants, especially large structural variants (SVs) including copy number variants (CNVs), which account for a substantial proportion of variation among individuals. Thus, we investigated the host genomes of individuals with severe/life-threatening COVID-19 at the level of large SVs (500bp-Mb level) to identify events that might provide insight into the inter-individual clinical variability in clinical course and outcomes of COVID-19 patients.MethodsOptical genome mapping using Bionano’s Saphyr® system was performed on thirty-seven severely ill COVID-19 patients admitted to intensive care units (ICU). To extract candidate SVs, three distinct analyses were undertaken. First, an unbiased whole-genome analysis of SVs was performed to identify rare/unique genic SVs in these patients that did not appear in population datasets to determine candidate loci as decisive predisposing factors associated with severe COVID-19. Second, common SVs with a population frequency filter was interrogated for possible association with severe COVID-19 based on literature surveys. Third, genome-wide SV enrichment in severely ill patients versus the general population was investigated by calculating odds ratios to identify top-ranked genes/loci. Candidate SVs were confirmed using qPCR and an independent bioinformatics tool (FaNDOM).ResultsOur patient-centric investigation identified 11 SVs involving 38 genes implicated in three key host-viral interaction pathways: (1) innate immunity and inflammatory response, (2) airway resistance to pathogens, and (3) viral replication, spread, and RNA editing. These included seven rare/unique SVs (not present in the control dataset), identified in 24.3% (9/37) of patients, impacting up to 31 genes, of whichSTK26andDPP4are the most promising candidates. A duplication partially overlappingSTK26was corroborated with data showing upregulation of this gene in severely ill patients. Further, using a population frequency filter of less than 20% in the Bionano control dataset, four SVs involving seven genes were identified in 56.7% (21/37) of patients.ConclusionThis study is the first to systematically assess and highlight SVs’ potential role in the pathogenesis of COVID-19 severity. The genes implicated here identify novel SVs, especiallySTK26, and extend previous reports involving innate immunity and type I interferon response in the pathogenesis of COVID-19. Our study also shows that optical genome mapping can be a powerful tool to identify large SVs impacting disease outcomes with split survival and add valuable genomic information to the existing sequencing-based technology databases to understand the inter-individual variability associated with SARS-CoV-2 infections and COVID-19 mortality.

Extrachromosomal DNA (ecDNA) is prevalent in human cancers and mediates high expression of oncogenes through gene amplification and altered gene regulation1. Gene induction typically involves cis-regulatory elements that contact and activate genes on the same chromosome2,3. Here we show that ecDNA hubs—clusters of around 10–100 ecDNAs within the nucleus—enable intermolecular enhancer–gene interactions to promote oncogene overexpression. ecDNAs that encode multiple distinct oncogenes form hubs in diverse cancer cell types and primary tumours. Each ecDNA is more likely to transcribe the oncogene when spatially clustered with additional ecDNAs. ecDNA hubs are tethered by the bromodomain and extraterminal domain (BET) protein BRD4 in a MYC-amplified colorectal cancer cell line. The BET inhibitor JQ1 disperses ecDNA hubs and preferentially inhibits ecDNA-derived-oncogene transcription. The BRD4-bound PVT1 promoter is ectopically fused to MYC and duplicated in ecDNA, receiving promiscuous enhancer input to drive potent expression of MYC. Furthermore, the PVT1 promoter on an exogenous episome suffices to mediate gene activation in trans by ecDNA hubs in a JQ1-sensitive manner. Systematic silencing of ecDNA enhancers by CRISPR interference reveals intermolecular enhancer–gene activation among multiple oncogene loci that are amplified on distinct ecDNAs. Thus, protein-tethered ecDNA hubs enable intermolecular transcriptional regulation and may serve as units of oncogene function and cooperative evolution and as potential targets for cancer therapy. Extrachromosomal DNA (ecDNA) congregates in clusters called ecDNA hubs that promote intermolecular interactions between gene-regulatory regions and thereby amplify the expression of oncogenes such as MYC in cancer cell lines.

Oncogene amplification, a major driver of cancer pathogenicity, is often mediated through focal amplification of genomic segments. Recent results implicate extrachromosomal DNA (ecDNA) as the primary driver of focal copy number amplification (fCNA) - enabling gene amplification, rapid tumor evolution, and the rewiring of regulatory circuitry. Resolving an fCNA’s structure is a first step in deciphering the mechanisms of its genesis and the fCNA’s subsequent biological consequences. We introduce a computational method, AmpliconReconstructor (AR), for integrating optical mapping (OM) of long DNA fragments (>150 kb) with next-generation sequencing (NGS) to resolve fCNAs at single-nucleotide resolution. AR uses an NGS-derived breakpoint graph alongside OM scaffolds to produce high-fidelity reconstructions. After validating its performance through multiple simulation strategies, AR reconstructed fCNAs in seven cancer cell lines to reveal the complex architecture of ecDNA, a breakage-fusion-bridge and other complex rearrangements. By reconstructing the rearrangement signatures associated with an fCNA’s generative mechanism, AR enables a more thorough understanding of the origins of fCNAs., Focal copy number amplifications (fCNAs), which drive cancer pathogenicity, arise by a number of mechanisms and can be challenging to call. Here the authors present AmpliconReconstructor for precise and scalable fCNA reconstruction using optical mapping and next-generation sequencing data.

Oncogene amplification, a major driver of cancer pathogenicity, is often mediated through focal amplification of genomic segments. Recent results implicate extrachromosomal DNA (ecDNA) as the primary mechanism driving focal copy number amplification (fCNA) - enabling gene amplification, rapid tumor evolution, and the rewiring of regulatory circuitry. Resolving an fCNA’s structure is a first step in deciphering the mechanisms of its genesis and the subsequent biological consequences. Here, we introduce a powerful new computational method, AmpliconReconstructor (AR), for integrating optical mapping (OM) of long DNA fragments (>150kb) with next-generation sequencing (NGS) to resolve fCNAs at single-nucleotide resolution. AR uses an NGS-derived breakpoint graph alongside OM scaffolds to produce high-fidelity reconstructions. After validating performance by extensive simulations, we used AR to reconstruct fCNAs in seven cancer cell lines to reveal the complex architecture of ecDNA, breakage-fusion-bridge cycles, and other complex rearrangements. By distinguishing between chromosomal and extrachromosomal origins, and by reconstructing the rearrangement signatures associated with a given fCNA’s generative mechanism, AR enables a more thorough understanding of the origins of fCNAs, and their functional consequences.