Methylation of the chromatin modifier KMT2D by SMYD2 contributes to therapeutic response in hormone-dependent breast cancer.

Activating mutations in PIK3CA are frequently found in estrogen-receptor-positive (ER+) breast cancer, and the combination of the phosphatidylinositol 3-kinase (PI3K) inhibitor alpelisib with anti-ER inhibitors is approved for therapy. We have previously demonstrated that the PI3K pathway regulates ER activity through phosphorylation of the chromatin modifier KMT2D. Here, we discovered a methylation site on KMT2D, at K1330 directly adjacent to S1331, catalyzed by the lysine methyltransferase SMYD2. SMYD2 loss attenuates alpelisib-induced KMT2D chromatin binding and alpelisib-mediated changes in gene expression, including ER-dependent transcription. Knockdown or pharmacological inhibition of SMYD2 sensitizes breast cancer cells, patient-derived organoids, and tumors to PI3K/AKT inhibition and endocrine therapy in part through KMT2D K1330 methylation. Together, our findings uncover a regulatory crosstalk between post-translational modifications that fine-tunes KMT2D function at the chromatin. This provides a rationale for the use of SMYD2 inhibitors in combination with PI3Kα/AKT inhibitors in the treatment of ER+/PIK3CA mutant breast cancer.

Metastatic site influences driver gene function in pancreatic cancer.

Driver gene mutations can increase the metastatic potential of the primary tumor1-3, but their role in sustaining tumor growth at metastatic sites is poorly understood. A paradigm of such mutations is inactivation of SMAD4 - a transcriptional effector of TGFß signaling - which is a hallmark of multiple gastrointestinal malignancies4,5. SMAD4 inactivation mediates TGFß's remarkable anti- to pro-tumorigenic switch during cancer progression and can thus influence both tumor initiation and metastasis6-14. To determine whether metastatic tumors remain dependent on SMAD4 inactivation, we developed a mouse model of pancreatic ductal adenocarcinoma (PDAC) that enables Smad4 depletion in the pre-malignant pancreas and subsequent Smad4 reactivation in established metastases. As expected, Smad4 inactivation facilitated the formation of primary tumors that eventually colonized the liver and lungs. By contrast, Smad4 reactivation in metastatic disease had strikingly opposite effects depending on the tumor's organ of residence: suppression of liver metastases and promotion of lung metastases. Integrative multiomic analysis revealed organ-specific differences in the tumor cells' epigenomic state, whereby the liver and lungs harbored chromatin programs respectively dominated by the KLF and RUNX developmental transcription factors, with Klf4 depletion being sufficient to reverse Smad4's tumor-suppressive activity in liver metastases. Our results show how epigenetic states favored by the organ of residence can influence the function of driver genes in metastatic tumors. This organ-specific gene-chromatin interplay invites consideration of anatomical site in the interpretation of tumor genetics, with implications for the therapeutic targeting of metastatic disease.

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.

Jak2V617F Reversible Activation Shows Its Essential Requirement in Myeloproliferative Neoplasms.

Gain-of-function mutations activating JAK/STAT signaling are seen in the majority of patients with myeloproliferative neoplasms (MPN), most commonly JAK2V617F. Although clinically approved JAK inhibitors improve symptoms and outcomes in MPNs, remissions are rare, and mutant allele burden does not substantively change with chronic therapy. We hypothesized this is due to limitations of current JAK inhibitors to potently and specifically abrogate mutant JAK2 signaling. We therefore developed a conditionally inducible mouse model allowing for sequential activation, and then inactivation, of Jak2V617F from its endogenous locus using a combined Dre-rox/Cre-lox dual-recombinase system. Jak2V617F deletion abrogates MPN features, induces depletion of mutant-specific hematopoietic stem/progenitor cells, and extends overall survival to an extent not observed with pharmacologic JAK inhibition, including when cooccurring with somatic Tet2 loss. Our data suggest JAK2V617F represents the best therapeutic target in MPNs and demonstrate the therapeutic relevance of a dual-recombinase system to assess mutant-specific oncogenic dependencies in vivo. SIGNIFICANCE: Current JAK inhibitors to treat myeloproliferative neoplasms are ineffective at eradicating mutant cells. We developed an endogenously expressed Jak2V617F dual-recombinase knock-in/knock-out model to investigate Jak2V617F oncogenic reversion in vivo. Jak2V617F deletion abrogates MPN features and depletes disease-sustaining MPN stem cells, suggesting improved Jak2V617F targeting offers the potential for greater therapeutic efficacy. See related commentary by Celik and Challen, p. 701. This article is featured in Selected Articles from This Issue, p. 695.

An epigenetic barrier sets the timing of human neuronal maturation.

The pace of human brain development is highly protracted compared with most other species1-7. The maturation of cortical neurons is particularly slow, taking months to years to develop adult functions3-5. Remarkably, such protracted timing is retained in cortical neurons derived from human pluripotent stem cells (hPSCs) during in vitro differentiation or upon transplantation into the mouse brain4,8,9. Those findings suggest the presence of a cell-intrinsic clock setting the pace of neuronal maturation, although the molecular nature of this clock remains unknown. Here we identify an epigenetic developmental programme that sets the timing of human neuronal maturation. First, we developed a hPSC-based approach to synchronize the birth of cortical neurons in vitro which enabled us to define an atlas of morphological, functional and molecular maturation. We observed a slow unfolding of maturation programmes, limited by the retention of specific epigenetic factors. Loss of function of several of those factors in cortical neurons enables precocious maturation. Transient inhibition of EZH2, EHMT1 and EHMT2 or DOT1L, at progenitor stage primes newly born neurons to rapidly acquire mature properties upon differentiation. Thus our findings reveal that the rate at which human neurons mature is set well before neurogenesis through the establishment of an epigenetic barrier in progenitor cells. Mechanistically, this barrier holds transcriptional maturation programmes in a poised state that is gradually released to ensure the prolonged timeline of human cortical neuron maturation.

Disruption of SUV39H1-Mediated H3K9 Methylation Sustains CAR T-cell Function.

Suboptimal functional persistence limits the efficacy of adoptive T-cell therapies. CD28-based chimeric antigen receptors (CAR) impart potent effector function to T cells but with a limited lifespan. We show here that the genetic disruption of SUV39H1, which encodes a histone-3, lysine-9 methyl-transferase, enhances the early expansion, long-term persistence, and overall antitumor efficacy of human CAR T cells in leukemia and prostate cancer models. Persisting SUV39H1-edited CAR T cells demonstrate improved expansion and tumor rejection upon multiple rechallenges. Transcriptional and genome accessibility profiling of repeatedly challenged CAR T cells shows improved expression and accessibility of memory transcription factors in SUV39H1-edited CAR T cells. SUV39H1 editing also reduces expression of inhibitory receptors and limits exhaustion in CAR T cells that have undergone multiple rechallenges. Our findings thus demonstrate the potential of epigenetic programming of CAR T cells to balance their function and persistence for improved adoptive cell therapies. SIGNIFICANCE: T cells engineered with CD28-based CARs possess robust effector function and antigen sensitivity but are hampered by limited persistence, which may result in tumor relapse. We report an epigenetic strategy involving disruption of the SUV39H1-mediated histone-silencing program that promotes the functional persistence of CD28-based CAR T cells. See related article by López-Cobo et al., p. 120. This article is featured in Selected Articles from This Issue, p. 5.

Development of Resistance to Type II JAK2 Inhibitors in MPN Depends on AXL Kinase and Is Targetable.

PURPOSE: Myeloproliferative neoplasms (MPN) dysregulate JAK2 signaling. Because clinical JAK2 inhibitors have limited disease-modifying effects, type II JAK2 inhibitors such as CHZ868 stabilizing inactive JAK2 and reducing MPN clones, gain interest. We studied whether MPN cells escape from type ll inhibition. EXPERIMENTAL DESIGN: MPN cells were continuously exposed to CHZ868. We used phosphoproteomic analyses and ATAC/RNA sequencing to characterize acquired resistance to type II JAK2 inhibition, and targeted candidate mediators in MPN cells and mice. RESULTS: MPN cells showed increased IC50 and reduced apoptosis upon CHZ868 reflecting acquired resistance to JAK2 inhibition. Among >2,500 differential phospho-sites, MAPK pathway activation was most prominent, while JAK2-STAT3/5 remained suppressed. Altered histone occupancy promoting AP-1/GATA binding motif exposure associated with upregulated AXL kinase and enriched RAS target gene profiles. AXL knockdown resensitized MPN cells and combined JAK2/AXL inhibition using bemcentinib or gilteritinib reduced IC50 to levels of sensitive cells. While resistant cells induced tumor growth in NOD/SCID gamma mice despite JAK2 inhibition, JAK2/AXL inhibition largely prevented tumor progression. Because inhibitors of MAPK pathway kinases such as MEK are clinically used in other malignancies, we evaluated JAK2/MAPK inhibition with trametinib to interfere with AXL/MAPK-induced resistance. Tumor growth was halted similarly to JAK2/AXL inhibition and in a systemic cell line-derived mouse model, marrow infiltration was decreased supporting dependency on AXL/MAPK. CONCLUSIONS: We report on a novel mechanism of AXL/MAPK-driven escape from type II JAK2 inhibition, which is targetable at different nodes. This highlights AXL as mediator of acquired resistance warranting inhibition to enhance sustainability of JAK2 inhibition in MPN.

A reversible epigenetic memory of inflammatory injury controls lineage plasticity and tumor initiation in the mouse pancreas.

Inflammation is essential to the disruption of tissue homeostasis and can destabilize the identity of lineage-committed epithelial cells. Here, we employ lineage-traced mouse models, single-cell transcriptomic and chromatin analyses, and CUT&TAG to identify an epigenetic memory of inflammatory injury in the pancreatic acinar cell compartment. Despite resolution of pancreatitis, our data show that acinar cells fail to return to their molecular baseline, with retention of elevated chromatin accessibility and H3K4me1 at metaplasia genes, such that memory represents an incomplete cell fate decision. In vivo, we find this epigenetic memory controls lineage plasticity, with diminished metaplasia in response to a second insult but increased tumorigenesis with an oncogenic Kras mutation. The lowered threshold for oncogenic transformation, in turn, can be restored by blockade of MAPK signaling. Together, we define the chromatin dynamics, molecular encoding, and recall of a prolonged epigenetic memory of inflammatory injury that impacts future responses but remains reversible.

An alternative NURF complex sustains acute myeloid leukemia by regulating the accessibility of insulator regions.

Efficient treatment of acute myeloid leukemia (AML) patients remains a challenge despite recent therapeutic advances. Here, using a CRISPRi screen targeting chromatin factors, we identified the nucleosome-remodeling factor (NURF) subunit BPTF as an essential regulator of AML cell survival. We demonstrate that BPTF forms an alternative NURF chromatin remodeling complex with SMARCA5 and BAP18, which regulates the accessibility of a large set of insulator regions in leukemic cells. This ensures efficient CTCF binding and boundary formation between topologically associated domains that is essential for maintaining the leukemic transcriptional programs. We also demonstrate that the well-studied PHD2-BROMO chromatin reader domains of BPTF, while contributing to complex recruitment to chromatin, are dispensable for leukemic cell growth. Taken together, our results uncover how the alternative NURF complex contributes to leukemia and provide a rationale for its targeting in AML.

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.

KMT2D suppresses Sonic hedgehog-driven medulloblastoma progression and metastasis.

The major cause of treatment failure and mortality among medulloblastoma patients is metastasis intracranially or along the spinal cord. The molecular mechanisms driving tumor metastasis in Sonic hedgehog-driven medulloblastoma (SHH-MB) patients, however, remain largely unknown. In this study we define a tumor suppressive role of KMT2D (MLL2), a gene frequently mutated in the most metastatic ß-subtype. Strikingly, genetic mouse models of SHH-MB demonstrate that heterozygous loss of Kmt2d in conjunction with activation of the SHH pathway causes highly penetrant disease with decreased survival, increased hindbrain invasion and spinal cord metastasis. Loss of Kmt2d attenuates neural differentiation and shifts the transcriptional/chromatin landscape of primary and metastatic tumors toward a decrease in differentiation genes and tumor suppressors and an increase in genes/pathways implicated in advanced stage cancer and metastasis (TGFß, Notch, Atoh1, Sox2, and Myc). Thus, secondary heterozygous KMT2D mutations likely have prognostic value for identifying SHH-MB patients prone to develop metastasis.

Interferon response and epigenetic modulation by SMARCA4 mutations drive ovarian tumor immunogenicity.

Cell-intrinsic mechanisms of immunogenicity in ovarian cancer (OC) are not well understood. The presence of damaging mutations in the SWI/SNF chromatin remodeling complex, such as the SMARCA4 (BRG1) catalytic subunit, has been associated with improved response to ICB, however the mechanism by which this occurs is unclear. The aim of this current study was to examine the alterations in tumor cell-intrinsic and extrinsic immune signaling caused by SMARCA4 loss. Using OC models with loss-of-function mutations in SMARCA4 , we found that SMARCA4 loss resulted in increased cancer cell-intrinsic immunogenicity, characterized by upregulation of long-terminal RNA repeats such as endogenous retroviruses, increased expression of interferon-stimulated genes, and upregulation of antigen presentation machinery. Notably, this response was dependent on IRF3 signaling, but was independent of the type I interferon receptor. Mice inoculated with cancer cells bearing SMARCA4 loss demonstrated increased activation of cytotoxic T cells and NK cells in the tumor microenvironment as well as increased infiltration with activated dendritic cells. These results were recapitulated when animals bearing SMARCA4- proficient tumors were treated with a BRG1 inhibitor, suggesting that modulation of chromatin remodeling through targeting SMARCA4 may serve as a strategy to reverse immune evasion in OC.

GABA Regulates Electrical Activity and Tumor Initiation in Melanoma.

Oncogenes can initiate tumors only in certain cellular contexts, which is referred to as oncogenic competence. In melanoma, whether cells in the microenvironment can endow such competence remains unclear. Using a combination of zebrafish transgenesis coupled with human tissues, we demonstrate that GABAergic signaling between keratinocytes and melanocytes promotes melanoma initiation by BRAFV600E. GABA is synthesized in melanoma cells, which then acts on GABA-A receptors in keratinocytes. Electron microscopy demonstrates specialized cell-cell junctions between keratinocytes and melanoma cells, and multielectrode array analysis shows that GABA acts to inhibit electrical activity in melanoma/keratinocyte cocultures. Genetic and pharmacologic perturbation of GABA synthesis abrogates melanoma initiation in vivo. These data suggest that GABAergic signaling across the skin microenvironment regulates the ability of oncogenes to initiate melanoma. SIGNIFICANCE: This study shows evidence of GABA-mediated regulation of electrical activity between melanoma cells and keratinocytes, providing a new mechanism by which the microenvironment promotes tumor initiation. This provides insights into the role of the skin microenvironment in early melanomas while identifying GABA as a potential therapeutic target in melanoma. See related commentary by Ceol, p. 2128. This article is featured in Selected Articles from This Issue, p. 2109.

Exportin 1 inhibition prevents neuroendocrine transformation through SOX2 down-regulation in lung and prostate cancers.

In lung and prostate adenocarcinomas, neuroendocrine (NE) transformation to an aggressive derivative resembling small cell lung cancer (SCLC) is associated with poor prognosis. We previously described dependency of SCLC on the nuclear transporter exportin 1. Here, we explored the role of exportin 1 in NE transformation. We observed up-regulated exportin 1 in lung and prostate pretransformation adenocarcinomas. Exportin 1 was up-regulated after genetic inactivation of TP53 and RB1 in lung and prostate adenocarcinoma cell lines, accompanied by increased sensitivity to the exportin 1 inhibitor selinexor in vitro. Exportin 1 inhibition prevented NE transformation in different TP53/RB1-inactivated prostate adenocarcinoma xenograft models that acquire NE features upon treatment with the aromatase inhibitor enzalutamide and extended response to the EGFR inhibitor osimertinib in a lung cancer transformation patient-derived xenograft (PDX) model exhibiting combined adenocarcinoma/SCLC histology. Ectopic SOX2 expression restored the enzalutamide-promoted NE phenotype on adenocarcinoma-to-NE transformation xenograft models despite selinexor treatment. Selinexor sensitized NE-transformed lung and prostate small cell carcinoma PDXs to standard cytotoxics. Together, these data nominate exportin 1 inhibition as a potential therapeutic target to constrain lineage plasticity and prevent or treat NE transformation in lung and prostate adenocarcinoma.

Chromatin regulation of transcriptional enhancers and cell fate by the Sotos syndrome gene NSD1.

Nuclear receptor-binding SET-domain protein 1 (NSD1), a methyltransferase that catalyzes H3K36me2, is essential for mammalian development and is frequently dysregulated in diseases, including Sotos syndrome. Despite the impacts of H3K36me2 on H3K27me3 and DNA methylation, the direct role of NSD1 in transcriptional regulation remains largely unknown. Here, we show that NSD1 and H3K36me2 are enriched at cis-regulatory elements, particularly enhancers. NSD1 enhancer association is conferred by a tandem quadruple PHD (qPHD)-PWWP module, which recognizes p300-catalyzed H3K18ac. By combining acute NSD1 depletion with time-resolved epigenomic and nascent transcriptomic analyses, we demonstrate that NSD1 promotes enhancer-dependent gene transcription by facilitating RNA polymerase II (RNA Pol II) pause release. Notably, NSD1 can act as a transcriptional coactivator independent of its catalytic activity. Moreover, NSD1 enables the activation of developmental transcriptional programs associated with Sotos syndrome pathophysiology and controls embryonic stem cell (ESC) multilineage differentiation. Collectively, we have identified NSD1 as an enhancer-acting transcriptional coactivator that contributes to cell fate transition and Sotos syndrome development.

Epigenetic dysregulation from chromosomal transit in micronuclei.

Chromosomal instability (CIN) and epigenetic alterations are characteristics of advanced and metastatic cancers1-4, but whether they are mechanistically linked is unknown. Here we show that missegregation of mitotic chromosomes, their sequestration in micronuclei5,6 and subsequent rupture of the micronuclear envelope7 profoundly disrupt normal histone post-translational modifications (PTMs), a phenomenon conserved across humans and mice, as well as in cancer and non-transformed cells. Some of the changes in histone PTMs occur because of the rupture of the micronuclear envelope, whereas others are inherited from mitotic abnormalities before the micronucleus is formed. Using orthogonal approaches, we demonstrate that micronuclei exhibit extensive differences in chromatin accessibility, with a strong positional bias between promoters and distal or intergenic regions, in line with observed redistributions of histone PTMs. Inducing CIN causes widespread epigenetic dysregulation, and chromosomes that transit in micronuclei experience heritable abnormalities in their accessibility long after they have been reincorporated into the primary nucleus. Thus, as well as altering genomic copy number, CIN promotes epigenetic reprogramming and heterogeneity in cancer.

MLL3 regulates the CDKN2A tumor suppressor locus in liver cancer.

Mutations in genes encoding components of chromatin modifying and remodeling complexes are among the most frequently observed somatic events in human cancers. For example, missense and nonsense mutations targeting the mixed lineage leukemia family member 3 (MLL3, encoded by KMT2C) histone methyltransferase occur in a range of solid tumors, and heterozygous deletions encompassing KMT2C occur in a subset of aggressive leukemias. Although MLL3 loss can promote tumorigenesis in mice, the molecular targets and biological processes by which MLL3 suppresses tumorigenesis remain poorly characterized. Here, we combined genetic, epigenomic, and animal modeling approaches to demonstrate that one of the mechanisms by which MLL3 links chromatin remodeling to tumor suppression is by co-activating the Cdkn2a tumor suppressor locus. Disruption of Kmt2c cooperates with Myc overexpression in the development of murine hepatocellular carcinoma (HCC), in which MLL3 binding to the Cdkn2a locus is blunted, resulting in reduced H3K4 methylation and low expression levels of the locus-encoded tumor suppressors p16/Ink4a and p19/Arf. Conversely, elevated KMT2C expression increases its binding to the CDKN2A locus and co-activates gene transcription. Endogenous Kmt2c restoration reverses these chromatin and transcriptional effects and triggers Ink4a/Arf-dependent apoptosis. Underscoring the human relevance of this epistasis, we found that genomic alterations in KMT2C and CDKN2A were associated with similar transcriptional profiles in human HCC samples. These results collectively point to a new mechanism for disrupting CDKN2A activity during cancer development and, in doing so, link MLL3 to an established tumor suppressor network.

Oncogenic IDH mutations increase heterochromatin-related replication stress without impacting homologous recombination.

Oncogenic mutations in isocitrate dehydrogenases 1 and 2 (IDH1/2) produce 2-hydroxyglutarate (2HG), which inhibits dioxygenases that modulate chromatin dynamics. The effects of 2HG have been reported to sensitize IDH tumors to poly-(ADP-ribose) polymerase (PARP) inhibitors. However, unlike PARP-inhibitor-sensitive BRCA1/2 tumors, which exhibit impaired homologous recombination, IDH-mutant tumors have a silent mutational profile and lack signatures associated with impaired homologous recombination. Instead, 2HG-producing IDH mutations lead to a heterochromatin-dependent slowing of DNA replication accompanied by increased replication stress and DNA double-strand breaks. This replicative stress manifests as replication fork slowing, but the breaks are repaired without a significant increase in mutation burden. Faithful resolution of replicative stress in IDH-mutant cells is dependent on poly-(ADP-ribosylation). Treatment with PARP inhibitors increases DNA replication but results in incomplete DNA repair. These findings demonstrate a role for PARP in the replication of heterochromatin and further validate PARP as a therapeutic target in IDH-mutant tumors.

Epigenetic plasticity cooperates with cell-cell interactions to direct pancreatic tumorigenesis.

The response to tumor-initiating inflammatory and genetic insults can vary among morphologically indistinguishable cells, suggesting as yet uncharacterized roles for epigenetic plasticity during early neoplasia. To investigate the origins and impact of such plasticity, we performed single-cell analyses on normal, inflamed, premalignant, and malignant tissues in autochthonous models of pancreatic cancer. We reproducibly identified heterogeneous cell states that are primed for diverse, late-emerging neoplastic fates and linked these to chromatin remodeling at cell-cell communication loci. Using an inference approach, we revealed signaling gene modules and tissue-level cross-talk, including a neoplasia-driving feedback loop between discrete epithelial and immune cell populations that was functionally validated in mice. Our results uncover a neoplasia-specific tissue-remodeling program that may be exploited for pancreatic cancer interception.

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.