Coordinated inheritance of extrachromosomal DNAs in cancer cells.

The chromosomal theory of inheritance dictates that genes on the same chromosome segregate together while genes on different chromosomes assort independently1. Extrachromosomal DNAs (ecDNAs) are common in cancer and drive oncogene amplification, dysregulated gene expression and intratumoural heterogeneity through random segregation during cell division2,3. Distinct ecDNA sequences, termed ecDNA species, can co-exist to facilitate intermolecular cooperation in cancer cells4. How multiple ecDNA species within a tumour cell are assorted and maintained across somatic cell generations is unclear. Here we show that cooperative ecDNA species are coordinately inherited through mitotic co-segregation. Imaging and single-cell analyses show that multiple ecDNAs encoding distinct oncogenes co-occur and are correlated in copy number in human cancer cells. ecDNA species are coordinately segregated asymmetrically during mitosis, resulting in daughter cells with simultaneous copy-number gains in multiple ecDNA species before any selection. Intermolecular proximity and active transcription at the start of mitosis facilitate the coordinated segregation of ecDNA species, and transcription inhibition reduces co-segregation. Computational modelling reveals the quantitative principles of ecDNA co-segregation and co-selection, predicting their observed distributions in cancer cells. Coordinated inheritance of ecDNAs enables co-amplification of specialized ecDNAs containing only enhancer elements and guides therapeutic strategies to jointly deplete cooperating ecDNA oncogenes. Coordinated inheritance of ecDNAs confers stability to oncogene cooperation and novel gene regulatory circuits, allowing winning combinations of epigenetic states to be transmitted across cell generations.

Extrachromosomal DNA amplifications in cancer.

Extrachromosomal DNA (ecDNA) amplification is an important driver alteration in cancer. It has been observed in most cancer types and is associated with worse patient outcome. The functional impact of ecDNA has been linked to its unique properties, such as its circular structure that is associated with altered chromatinization and epigenetic regulatory landscape, as well as its ability to randomly segregate during cell division, which fuels intercellular copy number heterogeneity. Recent investigations suggest that ecDNA is structurally more complex than previously anticipated and that it localizes to specialized nuclear bodies (hubs) and can act in trans as an enhancer for genes on other ecDNAs or chromosomes. In this Review, we synthesize what is currently known about how ecDNA is generated and how its genetic and epigenetic architecture affects proto-oncogene deregulation in cancer. We discuss how recently identified ecDNA functions may impact oncogenesis but also serve as new therapeutic vulnerabilities in cancer.

Reconstructing extrachromosomal DNA structural heterogeneity from long-read sequencing data using Decoil.

Circular extrachromosomal DNA (ecDNA) is a form of oncogene amplification found across cancer types and associated with poor outcome in patients. ecDNA can be structurally complex and can contain rearranged DNA sequences derived from multiple chromosome locations. As the structure of ecDNA can impact oncogene regulation and may indicate mechanisms of its formation, disentangling it at high resolution from sequencing data is essential. Even though methods have been developed to identify and reconstruct ecDNA in cancer genome sequencing, it remains challenging to resolve complex ecDNA structures, in particular amplicons with shared genomic footprints. We here introduce Decoil, a computational method that combines a breakpoint-graph approach with LASSO regression to reconstruct complex ecDNA and deconvolve co-occurring ecDNA elements with overlapping genomic footprints from long-read nanopore sequencing. Decoil outperforms de novo assembly and alignment-based methods in simulated long-read sequencing data for both simple and complex ecDNAs. Applying Decoil on whole-genome sequencing data uncovered different ecDNA topologies and explored ecDNA structure heterogeneity in neuroblastoma tumors and cell lines, indicating that this method may improve ecDNA structural analyses in cancer.

ecDNA hubs drive cooperative intermolecular oncogene expression.

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.

Reflection of neuroblastoma intratumor heterogeneity in the new OHC-NB1 disease model.

Accurate modeling of intratumor heterogeneity presents a bottleneck against drug testing. Flexibility in a preclinical platform is also desirable to support assessment of different endpoints. We established the model system, OHC-NB1, from a bone marrow metastasis from a patient diagnosed with MYCN-amplified neuroblastoma and performed whole-exome sequencing on the source metastasis and the different models and passages during model development (monolayer cell line, 3D spheroid culture and subcutaneous xenograft tumors propagated in mice). OHC-NB1 harbors a MYCN amplification in double minutes, 1p deletion, 17q gain and diploid karyotype, which persisted in all models. A total of 80-540 single-nucleotide variants (SNVs) was detected in each sample, and comparisons between the source metastasis and models identified 34 of 80 somatic SNVs to be propagated in the models. Clonal reconstruction using the combined copy number and SNV data revealed marked clonal heterogeneity in the originating metastasis, with four clones being reflected in the model systems. The set of OHC-NB1 models represents 43% of somatic SNVs and 23% of the cellularity in the originating metastasis with varying clonal compositions, indicating that heterogeneity is partially preserved in our model system.

Central memory phenotype drives success of checkpoint inhibition in combination with CAR T cells.

The immunosuppressive microenvironment in solid tumors is thought to form a barrier to the entry and efficacy of cell-based therapies such as chimeric antigen receptor (CAR) T cells. Combining CAR T cell therapy with checkpoint inhibitors has been demonstrated to oppose immune escape mechanisms in solid tumors and augment antitumor efficacy. We evaluated PD-1/PD-L1 signaling capacity and the impact of an inhibitor of this checkpoint axis in an in vitro system for cancer cell challenge, the coculture of L1CAM-specific CAR T cells with neuroblastoma cell lines. Fluorescence-activated cell sorting-based analyses and luciferase reporter assays were used to assess PD-1/PD-L1 expression on CAR T and tumor cells as well as CAR T cell ability to kill neuroblastoma cells. Coculturing neuroblastoma cell lines with L1CAM-CAR T cells upregulated PD-L1 expression on neuroblastoma cells, confirming adaptive immune resistance. Exposure to neuroblastoma cells also upregulated the expression of the PD-1/PD-L1 axis in CAR T cells. The checkpoint inhibitor, nivolumab, enhanced L1CAM-CAR T cell-directed killing. However, nivolumab-enhanced L1CAM-CAR T cell killing did not strictly correlate with PD-L1 expression on neuroblastoma cells. In fact, checkpoint inhibitor success relied on strong PD-1/PD-L1 axis expression in the CAR T cells, which in turn depended on costimulatory domains within the CAR construct, and more importantly, on the subset of T cells selected for CAR T cell generation. Thus, T cell subset selection for CAR T cell generation and CAR T cell prescreening for PD-1/PD-L1 expression could help determine when combination therapy with checkpoint inhibitors could improve treatment efficacy.

Forward genetic screen of human transposase genomic rearrangements.

BACKGROUND: Numerous human genes encode potentially active DNA transposases or recombinases, but our understanding of their functions remains limited due to shortage of methods to profile their activities on endogenous genomic substrates. RESULTS: To enable functional analysis of human transposase-derived genes, we combined forward chemical genetic hypoxanthine-guanine phosphoribosyltransferase 1 (HPRT1) screening with massively parallel paired-end DNA sequencing and structural variant genome assembly and analysis. Here, we report the HPRT1 mutational spectrum induced by the human transposase PGBD5, including PGBD5-specific signal sequences (PSS) that serve as potential genomic rearrangement substrates. CONCLUSIONS: The discovered PSS motifs and high-throughput forward chemical genomic screening approach should prove useful for the elucidation of endogenous genome remodeling activities of PGBD5 and other domesticated human DNA transposases and recombinases.

Intercellular extrachromosomal DNA copy-number heterogeneity drives neuroblastoma cell state diversity.

Neuroblastoma exhibits significant inter- and intra-tumor genetic heterogeneity and varying clinical outcomes. Extrachromosomal DNAs (ecDNAs) may drive this heterogeneity by independently segregating during cell division, leading to rapid oncogene amplification. While ecDNA-mediated oncogene amplification is linked to poor prognosis in various cancers, the effects of ecDNA copy-number heterogeneity on intermediate phenotypes are poorly understood. Here, we leverage DNA and RNA sequencing from the same single cells in cell lines and neuroblastoma patients to investigate these effects. By analyzing ecDNA amplicon structures, we reveal extensive intercellular ecDNA copy-number heterogeneity. We also provide direct evidence of how this heterogeneity influences the expression of cargo genes, including MYCN and its downstream targets, and the overall transcriptional state of neuroblastoma cells. Our findings highlight the role of ecDNA copy number in promoting rapid adaptability of cellular states within tumors, underscoring the need for ecDNA-specific treatment strategies to address tumor formation and adaptation.

Elimusertib has Antitumor Activity in Preclinical Patient-Derived Pediatric Solid Tumor Models.

The small-molecule inhibitor of ataxia telangiectasia and Rad3-related protein (ATR), elimusertib, is currently being tested clinically in various cancer entities in adults and children. Its preclinical antitumor activity in pediatric malignancies, however, is largely unknown. We here assessed the preclinical activity of elimusertib in 38 cell lines and 32 patient-derived xenograft (PDX) models derived from common pediatric solid tumor entities. Detailed in vitro and in vivo molecular characterization of the treated models enabled the evaluation of response biomarkers. Pronounced objective response rates were observed for elimusertib monotherapy in PDX, when treated with a regimen currently used in clinical trials. Strikingly, elimusertib showed stronger antitumor effects than some standard-of-care chemotherapies, particularly in alveolar rhabdomysarcoma PDX. Thus, elimusertib has strong preclinical antitumor activity in pediatric solid tumor models, which may translate to clinically meaningful responses in patients.

Passenger Gene Coamplifications Create Collateral Therapeutic Vulnerabilities in Cancer.

DNA amplifications in cancer do not only harbor oncogenes. We sought to determine whether passenger coamplifications could create collateral therapeutic vulnerabilities. Through an analysis of >3,000 cancer genomes followed by the interrogation of CRISPR-Cas9 loss-of-function screens across >700 cancer cell lines, we determined that passenger coamplifications are accompanied by distinct dependency profiles. In a proof-of-principle study, we demonstrate that the coamplification of the bona fide passenger gene DEAD-Box Helicase 1 (DDX1) creates an increased dependency on the mTOR pathway. Interaction proteomics identified tricarboxylic acid (TCA) cycle components as previously unrecognized DDX1 interaction partners. Live-cell metabolomics highlighted that this interaction could impair TCA activity, which in turn resulted in enhanced mTORC1 activity. Consequently, genetic and pharmacologic disruption of mTORC1 resulted in pronounced cell death in vitro and in vivo. Thus, structurally linked coamplification of a passenger gene and an oncogene can result in collateral vulnerabilities. SIGNIFICANCE: We demonstrate that coamplification of passenger genes, which were largely neglected in cancer biology in the past, can create distinct cancer dependencies. Because passenger coamplifications are frequent in cancer, this principle has the potential to expand target discovery in oncology. This article is featured in Selected Articles from This Issue, p. 384.

Doxorubicin Changes the Spatial Organization of the Genome around Active Promoters.

In this study, we delve into the impact of genotoxic anticancer drug treatment on the chromatin structure of human cells, with a particular focus on the effects of doxorubicin. Using Hi-C, ChIP-seq, and RNA-seq, we explore the changes in chromatin architecture brought about by doxorubicin and ICRF193. Our results indicate that physiologically relevant doses of doxorubicin lead to a local reduction in Hi-C interactions in certain genomic regions that contain active promoters, with changes in chromatin architecture occurring independently of Top2 inhibition, cell cycle arrest, and differential gene expression. Inside the regions with decreased interactions, we detected redistribution of RAD21 around the peaks of H3K27 acetylation. Our study also revealed a common structural pattern in the regions with altered architecture, characterized by two large domains separated from each other. Additionally, doxorubicin was found to increase CTCF binding in H3K27 acetylated regions. Furthermore, we discovered that Top2-dependent chemotherapy causes changes in the distance decay of Hi-C contacts, which are driven by direct and indirect inhibitors. Our proposed model suggests that doxorubicin-induced DSBs cause cohesin redistribution, which leads to increased insulation on actively transcribed TAD boundaries. Our findings underscore the significant impact of genotoxic anticancer treatment on the chromatin structure of the human genome.

Long-read sequencing of diagnosis and post-therapy medulloblastoma reveals complex rearrangement patterns and epigenetic signatures.

Cancer genomes harbor a broad spectrum of structural variants (SVs) driving tumorigenesis, a relevant subset of which escape discovery using short-read sequencing. We employed Oxford Nanopore Technologies (ONT) long-read sequencing in a paired diagnostic and post-therapy medulloblastoma to unravel the haplotype-resolved somatic genetic and epigenetic landscape. We assembled complex rearrangements, including a 1.55-Mbp chromothripsis event, and we uncover a complex SV pattern termed templated insertion (TI) thread, characterized by short (mostly <1 kb) insertions showing prevalent self-concatenation into highly amplified structures of up to 50 kbp in size. TI threads occur in 3% of cancers, with a prevalence up to 74% in liposarcoma, and frequent colocalization with chromothripsis. We also perform long-read-based methylome profiling and discover allele-specific methylation (ASM) effects, complex rearrangements exhibiting differential methylation, and differential promoter methylation in cancer-driver genes. Our study shows the advantage of long-read sequencing in the discovery and characterization of complex somatic rearrangements.

Mutational topography reflects clinical neuroblastoma heterogeneity.

Neuroblastoma is a pediatric solid tumor characterized by strong clinical heterogeneity. Although clinical risk-defining genomic alterations exist in neuroblastomas, the mutational processes involved in their generation remain largely unclear. By examining the topography and mutational signatures derived from all variant classes, we identified co-occurring mutational footprints, which we termed mutational scenarios. We demonstrate that clinical neuroblastoma heterogeneity is associated with differences in the mutational processes driving these scenarios, linking risk-defining pathognomonic variants to distinct molecular processes. Whereas high-risk MYCN-amplified neuroblastomas were characterized by signs of replication slippage and stress, homologous recombination-associated signatures defined high-risk non-MYCN-amplified patients. Non-high-risk neuroblastomas were marked by footprints of chromosome mis-segregation and TOP1 mutational activity. Furthermore, analysis of subclonal mutations uncovered differential activity of these processes through neuroblastoma evolution. Thus, clinical heterogeneity of neuroblastoma patients can be linked to differences in the mutational processes that are active in their tumors.

Single-cell profiling of alveolar rhabdomyosarcoma reveals RAS pathway inhibitors as cell-fate hijackers with therapeutic relevance.

Rhabdomyosarcoma (RMS) is a group of pediatric cancers with features of developing skeletal muscle. The cellular hierarchy and mechanisms leading to developmental arrest remain elusive. Here, we combined single-cell RNA sequencing, mass cytometry, and high-content imaging to resolve intratumoral heterogeneity of patient-derived primary RMS cultures. We show that the aggressive alveolar RMS (aRMS) subtype contains plastic muscle stem-like cells and cycling progenitors that drive tumor growth, and a subpopulation of differentiated cells that lost its proliferative potential and correlates with better outcomes. While chemotherapy eliminates cycling progenitors, it enriches aRMS for muscle stem-like cells. We screened for drugs hijacking aRMS toward clinically favorable subpopulations and identified a combination of RAF and MEK inhibitors that potently induces myogenic differentiation and inhibits tumor growth. Overall, our work provides insights into the developmental states underlying aRMS aggressiveness, chemoresistance, and progression and identifies the RAS pathway as a promising therapeutic target.

Coordinated inheritance of extrachromosomal DNA species in human cancer cells.

The chromosomal theory of inheritance has dominated human genetics, including cancer genetics. Genes on the same chromosome segregate together while genes on different chromosomes assort independently, providing a fundamental tenet of Mendelian inheritance. Extrachromosomal DNA (ecDNA) is a frequent event in cancer that drives oncogene amplification, dysregulated gene expression and intratumoral heterogeneity, including through random segregation during cell division. Distinct ecDNA sequences, herein termed ecDNA species, can co-exist to facilitate intermolecular cooperation in cancer cells. However, how multiple ecDNA species within a tumor cell are assorted and maintained across somatic cell generations to drive cancer cell evolution is not known. Here we show that cooperative ecDNA species can be coordinately inherited through mitotic co-segregation. Imaging and single-cell analyses show that multiple ecDNAs encoding distinct oncogenes co-occur and are correlated in copy number in human cancer cells. EcDNA species are coordinately segregated asymmetrically during mitosis, resulting in daughter cells with simultaneous copy number gains in multiple ecDNA species prior to any selection. Computational modeling reveals the quantitative principles of ecDNA co-segregation and co-selection, predicting their observed distributions in cancer cells. Finally, we show that coordinated inheritance of ecDNAs enables co-amplification of specialized ecDNAs containing only enhancer elements and guides therapeutic strategies to jointly deplete cooperating ecDNA oncogenes. Coordinated inheritance of ecDNAs confers stability to oncogene cooperation and novel gene regulatory circuits, allowing winning combinations of epigenetic states to be transmitted across cell generations.

Parallel sequencing of extrachromosomal circular DNAs and transcriptomes in single cancer cells.

Extrachromosomal DNAs (ecDNAs) are common in cancer, but many questions about their origin, structural dynamics and impact on intratumor heterogeneity are still unresolved. Here we describe single-cell extrachromosomal circular DNA and transcriptome sequencing (scEC&T-seq), a method for parallel sequencing of circular DNAs and full-length mRNA from single cells. By applying scEC&T-seq to cancer cells, we describe intercellular differences in ecDNA content while investigating their structural heterogeneity and transcriptional impact. Oncogene-containing ecDNAs were clonally present in cancer cells and drove intercellular oncogene expression differences. In contrast, other small circular DNAs were exclusive to individual cells, indicating differences in their selection and propagation. Intercellular differences in ecDNA structure pointed to circular recombination as a mechanism of ecDNA evolution. These results demonstrate scEC&T-seq as an approach to systematically characterize both small and large circular DNA in cancer cells, which will facilitate the analysis of these DNA elements in cancer and beyond.

Defining the landscape of circular RNAs in neuroblastoma unveils a global suppressive function of MYCN.

Circular RNAs (circRNAs) are a regulatory RNA class. While cancer-driving functions have been identified for single circRNAs, how they modulate gene expression in cancer is not well understood. We investigate circRNA expression in the pediatric malignancy, neuroblastoma, through deep whole-transcriptome sequencing in 104 primary neuroblastomas covering all risk groups. We demonstrate that MYCN amplification, which defines a subset of high-risk cases, causes globally suppressed circRNA biogenesis directly dependent on the DHX9 RNA helicase. We detect similar mechanisms in shaping circRNA expression in the pediatric cancer medulloblastoma implying a general MYCN effect. Comparisons to other cancers identify 25 circRNAs that are specifically upregulated in neuroblastoma, including circARID1A. Transcribed from the ARID1A tumor suppressor gene, circARID1A promotes cell growth and survival, mediated by direct interaction with the KHSRP RNA-binding protein. Our study highlights the importance of MYCN regulating circRNAs in cancer and identifies molecular mechanisms, which explain their contribution to neuroblastoma pathogenesis.

Shedding Light on ecDNA Behavior Using CRISPR-Based Live-Cell Imaging.

Extrachromosomal oncogene amplification on extrachromosomal DNA (ecDNA) has emerged as a hallmark of many cancers. In this issue, Yi and colleagues developed a CRISPR-based method for imaging ecDNA in live cells, termed ecTag. Using ecTag, the authors reveal important features of ecDNA in cancer cells such as their random mitotic segregation and clustering into transcriptionally active hubs after mitosis.See related article by Yi et al., p. 468.

The genomic and spatial mobility of extrachromosomal DNA and its implications for cancer therapy.

Extrachromosomal DNA (ecDNA) amplification has been observed in at least 30 different cancer types and is associated with worse patient outcomes. This has been linked to increased oncogene dosage because both oncogenes and associated enhancers can occupy ecDNA. New data challenge the view that only oncogene dosage is affected by ecDNA, and raises the possibility that ecDNA could disrupt genome-wide gene expression. Recent investigations suggest that ecDNA localizes to specialized nuclear bodies (hubs) in which they can act in trans as ectopic enhancers for genes on other ecDNA or chromosomes. Moreover, ecDNA can reintegrate into the genome, possibly further disrupting the gene regulatory landscape in tumor cells. In this Perspective, we discuss the emerging properties of ecDNA and highlight promising avenues to exploit this new knowledge for the development of ecDNA-directed therapies for cancer.