Exploring the landscape of focal amplifications in cancer using AmpliconArchitect

Extrachromosomal circular DNA landscape of breast cancer with lymph node metastasis

Breast cancer (BC) is a complex disease with diverse manifestations, often resulting in lymph node metastasis (LNM) and impacting patient prognosis. Extrachromosomal circular DNA (eccDNA) has emerged as a key player in tumorigenesis, yet its contribution to BC LNM remains elusive. Here, we examined primary tumors and matched LNM tissues from 19 BC patients using the Circle-Seq method. We identified a median count of 44,682 eccDNA in primary tumor tissues and 38,057 in their paired LNM tissues. Furthermore, a ladder-like size distribution is observed in both primary tumor and LNM tissues. Meanwhile, similar repeat sequence distribution and GC content are identified from both primary tissue and LNM tissues. Finally, we found that eccDNA from both groups are flanked with palindromic trinucleotide motifs. These observations indicate that eccDNA of primary tumor and LNM tissues are from similar chromosomal origins. However, a subset of miRNA-associated eccDNA displayed selective enrichment in metastatic lesions, such as miR-6730 and miR-548AA1 genes. This observation implicates the function of miRNA-related eccDNA in the metastatic cascade. Our study uncovers the potential significance of these unique eccDNA molecules, shedding light on their role in cancer metastasis.

Novel molecular subtypes of intracranial germ cell tumors expand therapeutic opportunities

BACKGROUND

Intracranial germ cell tumors (IGCTs) are a rare group of malignancies that are clinically classified as germinomas and nongerminomatous germ cell tumors (NGGCTs). Previous studies have found that somatic mutations involving the mitogen-activated protein kinase/mTOR signaling pathway are common early events. However, a comprehensive genomic understanding of IGCTs is still lacking.

METHODS

We established a cohort including over 100 IGCTs and conducted genomic and transcriptomic sequencing.

RESULTS

We identified novel recurrent driver genomic aberrations, including USP28 truncation mutations and high-level copy number amplification of KRAS and CRKL caused by replication of extrachromosomal DNA. Three distinct subtypes associated with unique genomic and clinical profiles were identified with transcriptome analysis: Immune-hot, MYC/E2F, and SHH. Both immune-hot and MYC/E2F were predominantly identified in germinomas and shared similar mutations involving the RAS/MAPK signaling pathway. However, the immune-hot group showed an older disease onset age and a significant immune response. MYC/E2F was characterized by a younger disease onset age and increased genomic instability, with a higher proportion of tumors showing whole-genome doubling. Additionally, the SHH subtype was mostly identified in NGGCTs.

CONCLUSIONS

Novel genomic aberrations and molecular subtypes were identified in IGCTs. These findings provide molecular basis for the potential introduction of new treatment strategies in this setting.

Imaging extrachromosomal DNA (ecDNA) in cancer

Extrachromosomal DNA (ecDNA) are circular regions of DNA that are found in many cancers. They are an important means of oncogene amplification, and correlate with treatment resistance and poor prognosis. Consequently, there is great interest in exploring and targeting ecDNA vulnerabilities as potential new therapeutic targets for cancer treatment. However, the biological significance of ecDNA and their associated regulatory control remains unclear. Light microscopy has been a central tool in the identification and characterisation of ecDNA. In this review we describe the different cellular models available to study ecDNA, and the imaging tools used to characterise ecDNA and their regulation. The insights gained from quantitative imaging are discussed in comparison with genome sequencing and computational approaches. We suggest that there is a crucial need for ongoing innovation using imaging if we are to achieve a full understanding of the dynamic regulation and organisation of ecDNA and their role in tumourigenesis.

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.

Structure and transcription of integrated HPV DNA in vulvar carcinomas

HPV infections are associated with a fraction of vulvar cancers. Through hybridization capture and DNA sequencing, HPV DNA was detected in five of thirteen vulvar cancers. HPV16 DNA was integrated into human DNA in three of the five. The insertions were in introns of human NCKAP1, C5orf67, and LRP1B. Integrations in NCKAP1 and C5orf67 were flanked by short direct repeats in the human DNA, consistent with HPV DNA insertions at sites of abortive, staggered, endonucleolytic incisions. The insertion in C5orf67 was present as a 36 kbp, human-HPV-hetero-catemeric DNA as either an extrachromosomal circle or a tandem repeat within the human genome. The human circularization/repeat junction was defined at single nucleotide resolution. The integrated viral DNA segments all retained an intact upstream regulatory region and the adjacent viral E6 and E7 oncogenes. RNA sequencing revealed that the only HPV genes consistently transcribed from the integrated viral DNAs were E7 and E6*I. The other two HPV DNA+ tumors had coinfections, but no evidence for integration. HPV-positive and HPV-negative vulvar cancers exhibited contrasting human, global gene expression patterns partially overlapping with previously observed differences between HPV-positive and HPV-negative cervical and oropharyngeal cancers. A substantial fraction of the differentially expressed genes involved immune system function. Thus, transcription and HPV DNA integration in vulvar cancers resemble those in other HPV-positive cancers. This study emphasizes the power of hybridization capture coupled with DNA and RNA sequencing to identify a broad spectrum of HPV types, determine human genome integration status of viral DNAs, and elucidate their structures.

Whole genome sequencing of HER2-positive metastatic extramammary Paget's disease: a case report

BACKGROUND

Extramammary Paget's disease (EMPD) is a rare cancer that occurs within the epithelium of the skin, arising predominantly in areas with high apocrine gland concentration such as the vulva, scrotum, penis and perianal regions. Here, we aim to integrate clinicopathological data with genomic analysis of aggressive, rapidly-progressing de novo metastatic EMPD responding to HER2-directed treatment in combination with other agents, to attain a more comprehensive understanding of the disease landscape.

METHODS

Immunohistochemical staining on the scrotal wall tumor and bone marrow metastasis demonstrated HER2 overexpression. Whole genome sequencing of the tumor and matched blood was performed.

RESULTS

Notable copy number gains (log2FC > 0.9) on chromosomes 7 and 8 were detected (n = 81), with 92.6% of these unique genes specifically located on chromosome 8. Prominent cancer-associated genes include ZNF703, HOOK3, DDHD2, LSM1, NSD3, ADAM9, BRF2, KAT6A and FGFR1. Interestingly, ERBB2 gene did not exhibit high copy number gain (log2FC = 0.4) although 90% of tumor cells stained HER2-positive. Enrichment in pathways associated with transforming growth factor-beta (TGFβ) (FDR = 0.0376, Enrichment Ratio = 8.12) and fibroblast growth factor receptor (FGFR1) signaling (FDR = 0.0082, Enrichment Ratio = 2.3) was detected. Amplicon structure analysis revealed that this was a simple-linear amplification event.

CONCLUSION

Whole genome sequencing revealed the underlying copy number variation landscape in HER2-positive metastatic EMPD. The presence of alternative signalling pathways and genetic variants suggests potential interactions with HER2 signalling, which possibly contributed to the HER2 overexpression and observed response to HER2-directed therapy combined with other agents in a comprehensive treatment regimen.

Disentangling oncogenic amplicons in esophageal adenocarcinoma

Esophageal adenocarcinoma is a prominent example of cancer characterized by frequent amplifications in oncogenes. However, the mechanisms leading to amplicons that involve breakage-fusion-bridge cycles and extrachromosomal DNA are poorly understood. Here, we use 710 esophageal adenocarcinoma cases with matched samples and patient-derived organoids to disentangle complex amplicons and their associated mechanisms. Short-read sequencing identifies ERBB2, MYC, MDM2, and HMGA2 as the most frequent oncogenes amplified in extrachromosomal DNAs. We resolve complex extrachromosomal DNA and breakage-fusion-bridge cycles amplicons by integrating of de-novo assemblies and DNA methylation in nine long-read sequenced cases. Complex amplicons shared between precancerous biopsy and late-stage tumor, an enrichment of putative enhancer elements and mobile element insertions are potential drivers of complex amplicons' origin. We find that patient-derived organoids recapitulate extrachromosomal DNA observed in the primary tumors and single-cell DNA sequencing capture extrachromosomal DNA-driven clonal dynamics across passages. Prospectively, long-read and single-cell DNA sequencing technologies can lead to better prediction of clonal evolution in esophageal adenocarcinoma.

Somatic structural variants drive distinct modes of oncogenesis in melanoma

The diversity of structural variants (SVs) in melanoma and how they impact oncogenesis are incompletely known. We performed harmonized analysis of SVs across melanoma histologic and genomic subtypes, and we identified distinct global properties between subtypes. These included the frequency and size of SVs and SV classes, their relation to chromothripsis events, and the impact on cancer-related genes of SVs that alter topologically associated domain (TAD) boundaries. Following our prior identification of double-stranded break repair deficiency in a subset of triple-wild-type cutaneous melanoma, we identified MRE11 and NBN loss-of-function SVs in melanomas with this mutational signature. Experimental knockouts of MRE11 and NBN, followed by olaparib cell viability assays in melanoma cells, indicated that dysregulation of each of these genes may cause sensitivity to PARP inhibitors in cutaneous melanomas. Broadly, harmonized analysis of melanoma SVs revealed distinct global genomic properties and molecular drivers, which may have biological and therapeutic impact.

Acquired Cross-Resistance in Small Cell Lung Cancer due to Extrachromosomal DNA Amplification of MYC Paralogs

Small cell lung cancer (SCLC) presents as a highly chemosensitive malignancy but acquires cross-resistance after relapse. This transformation is nearly inevitable in patients but has been difficult to capture in laboratory models. Here, we present a preclinical system that recapitulates acquired cross-resistance, developed from 51 patient-derived xenograft (PDX) models. Each model was tested in vivo against three clinical regimens: cisplatin plus etoposide, olaparib plus temozolomide, and topotecan. These drug-response profiles captured hallmark clinical features of SCLC, such as the emergence of treatment-refractory disease after early relapse. For one patient, serial PDX models revealed that cross-resistance was acquired through MYC amplification on extrachromosomal DNA (ecDNA). Genomic and transcriptional profiles of the full PDX panel revealed that MYC paralog amplifications on ecDNAs were recurrent in relapsed cross-resistant SCLC, and this was corroborated in tumor biopsies from relapsed patients. We conclude that ecDNAs with MYC paralogs are recurrent drivers of cross-resistance in SCLC.

SIGNIFICANCE

SCLC is initially chemosensitive, but acquired cross-resistance renders this disease refractory to further treatment and ultimately fatal. The genomic drivers of this transformation are unknown. We use a population of PDX models to discover that amplifications of MYC paralogs on ecDNA are recurrent drivers of acquired cross-resistance in SCLC. This article is featured in Selected Articles from This Issue, p. 695.

Characterization, biogenesis model, and current bioinformatics of human extrachromosomal circular DNA

Human extrachromosomal circular DNA, or eccDNA, has been the topic of extensive investigation in the last decade due to its prominent regulatory role in the development of disorders including cancer. With the rapid advancement of experimental, sequencing and computational technology, millions of eccDNA records are now accessible. Unfortunately, the literature and databases only provide snippets of this information, preventing us from fully understanding eccDNAs. Researchers frequently struggle with the process of selecting algorithms and tools to examine eccDNAs of interest. To explain the underlying formation mechanisms of the five basic classes of eccDNAs, we categorized their characteristics and functions and summarized eight biogenesis theories. Most significantly, we created a clear procedure to help in the selection of suitable techniques and tools and thoroughly examined the most recent experimental and bioinformatics methodologies and data resources for identifying, measuring and analyzing eccDNA sequences. In conclusion, we highlighted the current obstacles and prospective paths for eccDNA research, specifically discussing their probable uses in molecular diagnostics and clinical prediction, with an emphasis on the potential contribution of novel computational strategies.

Full-spectral genome analysis of natural killer/T cell lymphoma highlights impacts of genome instability in driving its progression

BACKGROUND

Natural killer/T cell lymphoma (NKTCL) is a clinically and genetically heterogeneous disease with poor prognosis. Genome sequencing and mutation characterization provides a powerful approach for patient stratification, treatment target discovery, and etiology identification. However, previous studies mostly concentrated on base-level mutations in primary NKTCL, whereas the large-scale genomic alterations in NKTCL and the mutational landscapes in relapsed/refractory NKTCL remain largely unexplored.

METHODS

Here, we assembled whole-genome sequencing and whole-exome sequencing data from 163 patients with primary or relapsed/refractory NKTCL and compared their somatic mutational landscapes at both nucleotide and structure levels.

RESULTS

Our study not only confirmed previously reported common NKTCL mutational targets like STAT3, TP53, and DDX3X but also unveiled several novel high-frequency mutational targets such as PRDM9, DST, and RBMX. In terms of the overall mutational landscape, we observed striking differences between primary and relapsed/refractory NKTCL patient groups, with the latter exhibits higher levels of tumor mutation burden, copy number variants (CNVs), and structural variants (SVs), indicating a strong signal of genomic instability. Complex structural rearrangements such as chromothripsis and focal amplification are also significantly enriched in relapsed/refractory NKTCL patients, exerting a substantial impact on prognosis. Accordingly, we devised a novel molecular subtyping system (i.e., C0-C4) with distinct prognosis by integrating potential driver mutations at both nucleotide and structural levels, which further provides an informative guidance for novel treatments that target these specific driver mutations and genome instability as a whole.

CONCLUSIONS

The striking differences underlying the mutational landscapes between the primary and relapsed/refractory NKTCL patients highlight the importance of genomic instability in driving the progression of NKTCL. Our newly proposed molecular subtyping system is valuable in assisting patient stratification and novel treatment design towards a better prognosis in the age of precision medicine.

Extrachromosomal DNA in cancer

Extrachromosomal DNA (ecDNA) has recently been recognized as a major contributor to cancer pathogenesis that is identified in most cancer types and is associated with poor outcomes. When it was discovered over 60 years ago, ecDNA was considered to be rare, and its impact on tumour biology was not well understood. The application of modern imaging and computational techniques has yielded powerful new insights into the importance of ecDNA in cancer. The non-chromosomal inheritance of ecDNA during cell division results in high oncogene copy number, intra-tumoural genetic heterogeneity and rapid tumour evolution that contributes to treatment resistance and shorter patient survival. In addition, the circular architecture of ecDNA results in altered patterns of gene regulation that drive elevated oncogene expression, potentially enabling the remodelling of tumour genomes. The generation of clusters of ecDNAs, termed ecDNA hubs, results in interactions between enhancers and promoters in trans, yielding a new paradigm in oncogenic transcription. In this Review, we highlight the rapid advancements in ecDNA research, providing new insights into ecDNA biogenesis, maintenance and transcription and its role in promoting tumour heterogeneity. To conclude, we delve into a set of unanswered questions whose answers will pave the way for the development of ecDNA targeted therapeutic approaches.

Deep whole-genome analysis of 494 hepatocellular carcinomas

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

scCircle-seq unveils the diversity and complexity of extrachromosomal circular DNAs in single cells

Extrachromosomal circular DNAs (eccDNAs) have emerged as important intra-cellular mobile genetic elements that affect gene copy number and exert in trans regulatory roles within the cell nucleus. Here, we describe scCircle-seq, a method for profiling eccDNAs and unraveling their diversity and complexity in single cells. We implement and validate scCircle-seq in normal and cancer cell lines, demonstrating that most eccDNAs vary largely between cells and are stochastically inherited during cell division, although their genomic landscape is cell type-specific and can be used to accurately cluster cells of the same origin. eccDNAs are preferentially produced from chromatin regions enriched in H3K9me3 and H3K27me3 histone marks and are induced during replication stress conditions. Concomitant sequencing of eccDNAs and RNA from the same cell uncovers the absence of correlation between eccDNA copy number and gene expression levels, except for a few oncogenes, including MYC, contained within a large eccDNA in colorectal cancer cells. Lastly, we apply scCircle-seq to one prostate cancer and two breast cancer specimens, revealing cancer-specific eccDNA landscapes and a higher propensity of eccDNAs to form in amplified genomic regions. scCircle-seq is a scalable tool that can be used to dissect the complexity of eccDNAs across different cell and tissue types, and further expands the potential of eccDNAs for cancer diagnostics.

Machine learning-based extrachromosomal DNA identification in large-scale cohorts reveals its clinical implications in cancer

The clinical implications of extrachromosomal DNA (ecDNA) in cancer therapy remain largely elusive. Here, we present a comprehensive analysis of ecDNA amplification spectra and their association with clinical and molecular features in multiple cohorts comprising over 13,000 pan-cancer patients. Using our developed computational framework, GCAP, and validating it with multifaceted approaches, we reveal a consistent pan-cancer pattern of mutual exclusivity between ecDNA amplification and microsatellite instability (MSI). In addition, we establish the role of ecDNA amplification as a risk factor and refine genomic subtypes in a cohort from 1015 colorectal cancer patients. Importantly, our investigation incorporates data from four clinical trials focused on anti-PD-1 immunotherapy, demonstrating the pivotal role of ecDNA amplification as a biomarker for guiding checkpoint blockade immunotherapy in gastrointestinal cancer. This finding represents clinical evidence linking ecDNA amplification to the effectiveness of immunotherapeutic interventions. Overall, our study provides a proof-of-concept of identifying ecDNA amplification from cancer whole-exome sequencing (WES) data, highlighting the potential of ecDNA amplification as a valuable biomarker for facilitating personalized cancer treatment.

Aberrant MET Receptor Tyrosine Kinase Signaling in Glioblastoma: Targeted Therapy and Future Directions

Brain tumors represent a heterogeneous group of neoplasms characterized by a high degree of aggressiveness and a poor prognosis. Despite recent therapeutic advances, the treatment of brain tumors, including glioblastoma (GBM), an aggressive primary brain tumor associated with poor prognosis and resistance to therapy, remains a significant challenge. Receptor tyrosine kinases (RTKs) are critical during development and in adulthood. Dysregulation of RTKs through activating mutations and gene amplification contributes to many human cancers and provides attractive therapeutic targets for treatment. Under physiological conditions, the Met RTK, the hepatocyte growth factor/scatter factor (HGF/SF) receptor, promotes fundamental signaling cascades that modulate epithelial-to-mesenchymal transition (EMT) involved in tissue repair and embryogenesis. In cancer, increased Met activity promotes tumor growth and metastasis by providing signals for proliferation, survival, and migration/invasion. Recent clinical genomic studies have unveiled multiple mechanisms by which MET is genetically altered in GBM, including focal amplification, chromosomal rearrangements generating gene fusions, and a splicing variant mutation (exon 14 skipping, METex14del). Notably, MET overexpression contributes to chemotherapy resistance in GBM by promoting the survival of cancer stem-like cells. This is linked to distinctive Met-induced pathways, such as the upregulation of DNA repair mechanisms, which can protect tumor cells from the cytotoxic effects of chemotherapy. The development of MET-targeted therapies represents a major step forward in the treatment of brain tumours. Preclinical studies have shown that MET-targeted therapies (monoclonal antibodies or small molecule inhibitors) can suppress growth and invasion, enhancing the efficacy of conventional therapies. Early-phase clinical trials have demonstrated promising results with MET-targeted therapies in improving overall survival for patients with recurrent GBM. However, challenges remain, including the need for patient stratification, the optimization of treatment regimens, and the identification of mechanisms of resistance. This review aims to highlight the current understanding of mechanisms underlying MET dysregulation in GBM. In addition, it will focus on the ongoing preclinical and clinical assessment of therapies targeting MET dysregulation in GBM.

eccDNA-pipe: an integrated pipeline for identification, analysis and visualization of extrachromosomal circular DNA from high-throughput sequencing data

Extrachromosomal circular DNA (eccDNA) is currently attracting considerable attention from researchers due to its significant impact on tumor biogenesis. High-throughput sequencing (HTS) methods for eccDNA identification are continually evolving. However, an efficient pipeline for the integrative and comprehensive analysis of eccDNA obtained from HTS data is still lacking. Here, we introduce eccDNA-pipe, an accessible software package that offers a user-friendly pipeline for conducting eccDNA analysis starting from raw sequencing data. This dataset includes data from various sequencing techniques such as whole-genome sequencing (WGS), Circle-seq and Circulome-seq, obtained through short-read sequencing or long-read sequencing. eccDNA-pipe presents a comprehensive solution for both upstream and downstream analysis, encompassing quality control and eccDNA identification in upstream analysis and downstream tasks such as eccDNA length distribution analysis, differential analysis of genes enriched with eccDNA and visualization of eccDNA structures. Notably, eccDNA-pipe automatically generates high-quality publication-ready plots. In summary, eccDNA-pipe provides a comprehensive and user-friendly pipeline for customized analysis of eccDNA research.

Innovative insights into extrachromosomal circular DNAs in gynecologic tumors and reproduction

Originating but free from chromosomal DNA, extrachromosomal circular DNAs (eccDNAs) are organized in circular form and have long been found in unicellular and multicellular eukaryotes. Their biogenesis and function are poorly understood as they are characterized by sequence homology with linear DNA, for which few detection methods are available. Recent advances in high-throughput sequencing technologies have revealed that eccDNAs play crucial roles in tumor formation, evolution, and drug resistance as well as aging, genomic diversity, and other biological processes, bringing it back to the research hotspot. Several mechanisms of eccDNA formation have been proposed, including the breakage-fusion-bridge (BFB) and translocation-deletion-amplification models. Gynecologic tumors and disorders of embryonic and fetal development are major threats to human reproductive health. The roles of eccDNAs in these pathological processes have been partially elucidated since the first discovery of eccDNA in pig sperm and the double minutes in ovarian cancer ascites. The present review summarized the research history, biogenesis, and currently available detection and analytical methods for eccDNAs and clarified their functions in gynecologic tumors and reproduction. We also proposed the application of eccDNAs as drug targets and liquid biopsy markers for prenatal diagnosis and the early detection, prognosis, and treatment of gynecologic tumors. This review lays theoretical foundations for future investigations into the complex regulatory networks of eccDNAs in vital physiological and pathological processes.

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.

Single cell multi-omics reveal intra-cell-line heterogeneity across human cancer cell lines

Human cancer cell lines have long served as tools for cancer research and drug discovery, but the presence and the source of intra-cell-line heterogeneity remain elusive. Here, we perform single-cell RNA-sequencing and ATAC-sequencing on 42 and 39 human cell lines, respectively, to illustrate both transcriptomic and epigenetic heterogeneity within individual cell lines. Our data reveal that transcriptomic heterogeneity is frequently observed in cancer cell lines of different tissue origins, often driven by multiple common transcriptional programs. Copy number variation, as well as epigenetic variation and extrachromosomal DNA distribution all contribute to the detected intra-cell-line heterogeneity. Using hypoxia treatment as an example, we demonstrate that transcriptomic heterogeneity could be reshaped by environmental stress. Overall, our study performs single-cell multi-omics of commonly used human cancer cell lines and offers mechanistic insights into the intra-cell-line heterogeneity and its dynamics, which would serve as an important resource for future cancer cell line-based studies.

Methods, bioinformatics tools and databases in ecDNA research: An overview

Extrachromosomal DNA (ecDNA), derived from chromosomes, is a cancer-specific circular DNA molecule. EcDNA drives tumor initiation and progression, which is associated with poor clinical outcomes and drug resistance in a wide range of cancers. Although ecDNA was first discovered in 1965, tremendous technological revolutions in recent years have provided crucial new insights into its key biological functions and regulatory mechanisms. Here, we provide a thorough overview of the methods, bioinformatics tools, and database resources used in ecDNA research, mainly focusing on their performance, strengths, and limitations. This study can provide important reference for selecting the most appropriate method in ecDNA research. Furthermore, we offer suggestions for the current bioinformatics analysis of ecDNA and provide an outlook to the future research.

Circular extrachromosomal DNA promotes tumor heterogeneity in high-risk medulloblastoma

Circular extrachromosomal DNA (ecDNA) in patient tumors is an important driver of oncogenic gene expression, evolution of drug resistance and poor patient outcomes. Applying computational methods for the detection and reconstruction of ecDNA across a retrospective cohort of 481 medulloblastoma tumors from 465 patients, we identify circular ecDNA in 82 patients (18%). Patients with ecDNA-positive medulloblastoma were more than twice as likely to relapse and three times as likely to die within 5 years of diagnosis. A subset of tumors harbored multiple ecDNA lineages, each containing distinct amplified oncogenes. Multimodal sequencing, imaging and CRISPR inhibition experiments in medulloblastoma models reveal intratumoral heterogeneity of ecDNA copy number per cell and frequent putative 'enhancer rewiring' events on ecDNA. This study reveals the frequency and diversity of ecDNA in medulloblastoma, stratified into molecular subgroups, and suggests copy number heterogeneity and enhancer rewiring as oncogenic features of ecDNA.

ReConPlot: an R package for the visualization and interpretation of genomic rearrangements

MOTIVATION

Whole-genome sequencing studies of human tumours have revealed that complex forms of structural variation, collectively known as complex genome rearrangements (CGRs), are pervasive across diverse cancer types. Detection, classification, and mechanistic interpretation of CGRs requires the visualization of complex patterns of somatic copy number aberrations (SCNAs) and structural variants (SVs). However, there is a lack of tools specifically designed to facilitate the visualization and study of CGRs.

RESULTS

We present ReConPlot (REarrangement and COpy Number PLOT), an R package that provides functionalities for the joint visualization of SCNAs and SVs across one or multiple chromosomes. ReConPlot is based on the popular ggplot2 package, thus allowing customization of plots and the generation of publication-quality figures with minimal effort. Overall, ReConPlot facilitates the exploration, interpretation, and reporting of CGR patterns.

AVAILABILITY AND IMPLEMENTATION

The R package ReConPlot is available at https://github.com/cortes-ciriano-lab/ReConPlot. Detailed documentation and a tutorial with examples are provided with the package.

ATACAmp: a tool for detecting ecDNA/HSRs from bulk and single-cell ATAC-seq data

BACKGROUND

High oncogene expression in cancer cells is a major cause of rapid tumor progression and drug resistance. Recent cancer genome research has shown that oncogenes as well as regulatory elements can be amplified in the form of extrachromosomal DNA (ecDNA) or subsequently integrated into chromosomes as homogeneously staining regions (HSRs). These genome-level variants lead to the overexpression of the corresponding oncogenes, resulting in poor prognosis. Most existing detection methods identify ecDNA using whole genome sequencing (WGS) data. However, these techniques usually detect many false positive regions owing to chromosomal DNA interference.

RESULTS

In the present study, an algorithm called "ATACAmp" that can identify ecDNA/HSRs in tumor genomes using ATAC-seq data has been described. High chromatin accessibility, one of the characteristics of ecDNA, makes ATAC-seq naturally enriched in ecDNA and reduces chromosomal DNA interference. The algorithm was validated using ATAC-seq data from cell lines that have been experimentally determined to contain ecDNA regions. ATACAmp accurately identified the majority of validated ecDNA regions. AmpliconArchitect, the widely used ecDNA detecting tool, was used to detect ecDNA regions based on the WGS data of the same cell lines. Additionally, the Circle-finder software, another tool that utilizes ATAC-seq data, was assessed. The results showed that ATACAmp exhibited higher accuracy than AmpliconArchitect and Circle-finder. Moreover, ATACAmp supported the analysis of single-cell ATAC-seq data, which linked ecDNA to specific cells.

CONCLUSIONS

ATACAmp, written in Python, is freely available on GitHub under the MIT license: https://github.com/chsmiss/ATAC-amp . Using ATAC-seq data, ATACAmp offers a novel analytical approach that is distinct from the conventional use of WGS data. Thus, this method has the potential to reduce the cost and technical complexity associated ecDNA analysis.

Extrachromosomal DNA (ecDNA): Unveiling its role in cancer progression and implications for early detection

Extrachromosomal DNA (ecDNA) is a special class of circular DNA in eukaryotes, which is independent of conventional chromosomes. These circular molecules play important roles in biology, especially in cancer biology. The emergence of sequencing technologies such as CCDA-seq and Amplicon Architect has led to a progressive unraveling of the mystery of ecDNA. Consequently, insights into its function and potential applications have begun to surface. Among these studies, the most noteworthy research pertains to cancer-related investigations into ecDNA. Numerous studies have underscored the significance of ecDNA in the pathogenesis of cancer and its role in accelerating cancer evolution. This review provides an overview of the source, structure, and function of ecDNA, while compiling recent advancements in ecDNA in the field of cancer. Nonetheless, further research is imperative to determine its effectiveness and specificity as a biomarker for early cancer detection.

MGA-seq: robust identification of extrachromosomal DNA and genetic variants using multiple genetic abnormality sequencing

Genomic abnormalities are strongly associated with cancer and infertility. In this study, we develop a simple and efficient method - multiple genetic abnormality sequencing (MGA-Seq) - to simultaneously detect structural variation, copy number variation, single-nucleotide polymorphism, homogeneously staining regions, and extrachromosomal DNA (ecDNA) from a single tube. MGA-Seq directly sequences proximity-ligated genomic fragments, yielding a dataset with concurrent genome three-dimensional and whole-genome sequencing information, enabling approximate localization of genomic structural variations and facilitating breakpoint identification. Additionally, by utilizing MGA-Seq, we map focal amplification and oncogene coamplification, thus facilitating the exploration of ecDNA's transcriptional regulatory function.

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.

Sequence Characterization of Extra-Chromosomal Circular DNA Content in Multiple Blackgrass (Alopecurus myosuroides) Populations

Alopecurus myosuroides (blackgrass) is a problematic weed of Western European winter wheat, and its success is largely due to widespread multiple-herbicide resistance. Previous analysis of F2 seed families derived from two distinct blackgrass populations exhibiting equivalent non-target site resistance (NTSR) phenotypes shows resistance is polygenic and evolves from standing genetic variation. Using a CIDER-seq pipeline, we show that herbicide-resistant (HR) and herbicide-sensitive (HS) F3 plants from these F2 seed families as well as the parent populations they were derived from carry extra-chromosomal circular DNA (eccDNA). We identify the similarities and differences in the coding structures within and between resistant and sensitive populations. Although the numbers and size of detected eccDNAs varied between the populations, comparisons between the HR and HS blackgrass populations identified shared and unique coding content, predicted genes, and functional protein domains. These include genes related to herbicide detoxification such as Cytochrome P450s, ATP-binding cassette transporters, and glutathione transferases including AmGSTF1. eccDNA content was mapped to the A. myosuroides reference genome, revealing genomic regions at the distal end of chromosome 5 and the near center of chromosomes 1 and 7 as regions with a high number of mapped eccDNA gene density. Mapping to 15 known herbicide-resistant QTL regions showed that the eccDNA coding sequences matched twelve, with four QTL matching HS coding sequences; only one region contained HR coding sequences. These findings establish that, like other pernicious weeds, blackgrass has eccDNAs that contain homologs of chromosomal genes, and these may contribute genetic heterogeneity and evolutionary innovation to rapidly adapt to abiotic stresses, including herbicide treatment.

Identifying predictors of glioma evolution from longitudinal sequencing

Clonal evolution drives cancer progression and therapeutic resistance. Recent studies have revealed divergent longitudinal trajectories in gliomas, but early molecular features steering posttreatment cancer evolution remain unclear. Here, we collected sequencing and clinical data of initial-recurrent tumor pairs from 544 adult diffuse gliomas and performed multivariate analysis to identify early molecular predictors of tumor evolution in three diffuse glioma subtypes. We found that CDKN2A deletion at initial diagnosis preceded tumor necrosis and microvascular proliferation that occur at later stages of IDH-mutant glioma. Ki67 expression at diagnosis was positively correlated with acquiring hypermutation at recurrence in the IDH-wild-type glioma. In all glioma subtypes, MYC gain or MYC-target activation at diagnosis was associated with treatment-induced hypermutation at recurrence. To predict glioma evolution, we constructed CELLO2 (Cancer EvoLution for LOngitudinal data version 2), a machine learning model integrating features at diagnosis to forecast hypermutation and progression after treatment. CELLO2 successfully stratified patients into subgroups with distinct prognoses and identified a high-risk patient group featured by MYC gain with worse post-progression survival, from the low-grade IDH-mutant-noncodel subtype. We then performed chronic temozolomide-induction experiments in glioma cell lines and isogenic patient-derived gliomaspheres and demonstrated that MYC drives temozolomide resistance by promoting hypermutation. Mechanistically, we demonstrated that, by binding to open chromatin and transcriptionally active genomic regions, c-MYC increases the vulnerability of key mismatch repair genes to treatment-induced mutagenesis, thus triggering hypermutation. This study reveals early predictors of cancer evolution under therapy and provides a resource for precision oncology targeting cancer dynamics in diffuse gliomas.

FLED: a full-length eccDNA detector for long-reads sequencing data

Reconstructing the full-length sequence of extrachromosomal circular DNA (eccDNA) from short sequencing reads has proved challenging given the similarity of eccDNAs and their corresponding linear DNAs. Previous sequencing methods were unable to achieve high-throughput detection of full-length eccDNAs. Herein, a novel algorithm was developed, called Full-Length eccDNA Detection (FLED), to reconstruct the sequence of eccDNAs based on the strategy that combined rolling circle amplification and nanopore long-reads sequencing technology. Seven human epithelial and cancer cell line samples were analyzed by FLED and over 5000 full-length eccDNAs were identified per sample. The structures of identified eccDNAs were validated by both Polymerase Chain Reaction (PCR) and Sanger sequencing. Compared to other published nanopore-based eccDNA detectors, FLED exhibited higher sensitivity. In cancer cell lines, the genes overlapped with eccDNA regions were enriched in cancer-related pathways and cis-regulatory elements can be predicted in the upstream or downstream of intact genes on eccDNA molecules, and the expressions of these cancer-related genes were dysregulated in tumor cell lines, indicating the regulatory potency of eccDNAs in biological processes. The proposed method takes advantage of nanopore long reads and enables unbiased reconstruction of full-length eccDNA sequences. FLED is implemented using Python3 which is freely available on GitHub (https://github.com/FuyuLi/FLED).

Deciphering complex breakage-fusion-bridge genome rearrangements with Ambigram

Breakage-fusion-bridge (BFB) is a complex rearrangement that leads to tumor malignancy. Existing models for detecting BFBs rely on the ideal BFB hypothesis, ruling out the possibility of BFBs entangled with other structural variations, that is, complex BFBs. We propose an algorithm Ambigram to identify complex BFB and reconstruct the rearranged structure of the local genome during the cancer subclone evolution process. Ambigram handles data from short, linked, long, and single-cell sequences, and optical mapping technologies. Ambigram successfully deciphers the gold- or silver-standard complex BFBs against the state-of-the-art in multiple cancers. Ambigram dissects the intratumor heterogeneity of complex BFB events with single-cell reads from melanoma and gastric cancer. Furthermore, applying Ambigram to liver and cervical cancer data suggests that the BFB mechanism may mediate oncovirus integrations. BFB also exists in noncancer genomics. Investigating the complete human genome reference with Ambigram suggests that the BFB mechanism may be involved in two genome reorganizations of Homo Sapiens during evolution. Moreover, Ambigram discovers the signals of recurrent foldback inversions and complex BFBs in whole genome data from the 1000 genome project, and congenital heart diseases, respectively.

Small extrachromosomal circular DNAs as biomarkers for multi-cancer diagnosis and monitoring

BACKGROUND

Small extrachromosomal circular DNAs (eccDNAs) have the potential to be cancer biomarkers. However, the formation mechanisms and functions of small eccDNAs selected in carcinogenesis are not clear, and whether the small eccDNA profile in the plasma of cancer patients represents that in cancer tissues remains to be elucidated.

METHODS

A novel sequencing workflow based on the nanopore sequencing platform was used to sequence naturally existing full-length small eccDNAs in tissues and plasma collected from 25 cancer patients (including prostate cancer, hepatocellular carcinoma and colorectal cancer), and from an independent validation cohort (including 7 cancer plasma and 14 healthy plasma).

RESULTS

Compared with those in non-cancer tissues, small eccDNAs detected in cancer tissues had a significantly larger number and size (P = 0.040 and 2.2e-16, respectively), along with more even distribution and different formation mechanisms. Although small eccDNAs had different general characteristics and genomic annotation between cancer tissues and the paired plasma, they had similar formation mechanisms and cancer-related functions. Small eccDNAs originated from some specific genes had great multi-cancer diagnostic value in tissues (AUC ≥ 0.8) and plasma (AUC > 0.9), especially increasing the accuracy of multi-cancer prediction of CEA/CA19-9 levels. The high multi-cancer diagnostic value of small eccDNAs originated from ALK&ETV6 could be extrapolated from tissues (AUC = 0.804) to plasma and showed high positive predictive value (100%) and negative predictive value (82.35%) in a validation cohort.

CONCLUSIONS

As independent and stable circular DNA molecules, small eccDNAs in both tissues and plasma can be used as ideal biomarkers for cost-effective multi-cancer diagnosis and monitoring.

Variation of extrachromosomal circular DNA in cancer cell lines

The presence of oncogene carrying eccDNAs is strongly associated with carcinogenesis and poor patient survival. Tumour biopsies and in vitro cancer cell lines are frequently utilized as models to investigate the role of eccDNA in cancer. However, eccDNAs are often lost during the in vitro growth of cancer cell lines, questioning the reproducibility of studies utilizing cancer cell line models. Here, we conducted a comprehensive analysis of eccDNA variability in seven cancer cell lines (MCA3D, PDV, HaCa4, CarC, MIA-PaCa-2, AsPC-1, and PC-3). We compared the content of unique eccDNAs between triplicates of each cell line and found that the number of unique eccDNA is specific to each cell line, while the eccDNA sequence content varied greatly among triplicates (∼ 0-1% eccDNA coordinate commonality). In the PC-3 cell line, we found that the large eccDNA (ecDNA) with MYC is present in high-copy number in an NCI cell line isolate but not present in ATCC isolates. Together, these results reveal that the sequence content of eccDNA is highly variable in cancer cell lines. This highlights the importance of testing cancer cell lines before use, and to enrich for subclones in cell lines with the desired eccDNA to get relatively pure population for studying the role of eccDNA in cancer.

Single-cell third-generation sequencing-based multi-omics uncovers gene expression changes governed by ecDNA and structural variants in cancer cells

BACKGROUND

Cancer cells often exhibit large-scale genomic variations, such as circular extrachromosomal DNA (ecDNA) and structural variants (SVs), which have been highly correlated with the initiation and progression of cancer. Currently, no adequate method exists to unveil how these variations regulate gene expression in heterogeneous cancer cell populations at a single-cell resolution.

METHODS

Here, we developed a single-cell multi-omics sequencing method, scGTP-seq, to analyse ecDNA and SVs using long-read sequencing technologies.

RESULTS AND CONCLUSIONS

We demonstrated that our method can efficiently detect ecDNA and SVs and illustrated how these variations affect transcriptomic changes in various cell lines. Finally, we applied and validated this method in a clinical sample of hepatocellular carcinoma (HCC), demonstrating a feasible way to monitor the evolution of ecDNA and SVs during cancer progression.

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.

Integrative profiling of extrachromosomal circular DNA in placenta and maternal plasma provides insights into the biology of fetal growth restriction and reveals potential biomarkers

Introduction: Fetal growth restriction (FGR) is a placenta-mediated pregnancy complication that predisposes fetuses to perinatal complications. Maternal plasma cell-free DNA harbors DNA originating from placental trophoblasts, which is promising for the prenatal diagnosis and prediction of pregnancy complications. Extrachromosomal circular DNA (eccDNA) is emerging as an ideal biomarker and target for several diseases. Methods: We utilized eccDNA sequencing and bioinformatic pipeline to investigate the characteristics and associations of eccDNA in placenta and maternal plasma, the role of placental eccDNA in the pathogenesis of FGR, and potential plasma eccDNA biomarkers of FGR. Results: Using our bioinformatics pipelines, we identified multi-chromosomal-fragment and single-fragment eccDNA in placenta, but almost exclusively single-fragment eccDNA in maternal plasma. Relative to that in plasma, eccDNA in placenta was larger and substantially more abundant in exons, untranslated regions, promoters, repetitive elements [short interspersed nuclear elements (SINEs)/Alu, SINEs/mammalian-wide interspersed repeats, long terminal repeats/endogenous retrovirus-like elements, and single recognition particle RNA], and transcription factor binding motifs. Placental multi-chromosomal-fragment eccDNA was enriched in confident enhancer regions predicted to pertain to genes in apoptosis, energy, cell growth, and autophagy pathways. Placental eccDNA-associated genes whose abundance differed between the FGR and control groups were associated with immunity-related gene ontology (GO) terms. The combined analysis of plasma and placental eccDNA-associated genes in the FGR and control groups led to the identification of potential biomarkers that were assigned to the GO terms of the epigenetic regulation of gene expression and nutrient-related processes, respectively. Conclusion: Together, our results highlight links between placenta functions and multi-chromosomal-fragment and single-fragment eccDNA. The integrative analysis of placental and plasma eccDNA confirmed the potential of these molecules as disease-specific biomarkers of FGR.

Extrachromosomal DNA (ecDNA) in cancer: mechanisms, functions, and clinical implications

Extrachromosomal DNA (ecDNA) is circular DNA that plays an important role in the development and heterogeneity of cancer. The rapid evolution of methods to detect ecDNA, including microscopic and sequencing approaches, has greatly enhanced our knowledge of the role of ecDNA in cancer development and evolution. Here, we review the molecular characteristics, functions, mechanisms of formation, and detection methods of ecDNA, with a focus on the potential clinical implications of ecDNA in cancer. Specifically, we consider the role of ecDNA in acquired drug resistance, as a diagnostic and prognostic biomarker, and as a therapeutic target in the context of cancer. As the pathological and clinical significance of ecDNA continues to be explored, it is anticipated that ecDNA will have broad applications in the diagnosis, prognosis, and treatment of patients with cancer.

Acquired Cross-resistance in Small Cell Lung Cancer due to Extrachromosomal DNA Amplification of MYC paralogs

Small cell lung cancer (SCLC) presents as a highly chemosensitive malignancy but acquires cross-resistance after relapse. This transformation is nearly inevitable in patients but has been difficult to capture in laboratory models. Here we present a pre-clinical system that recapitulates acquired cross-resistance in SCLC, developed from 51 patient-derived xenografts (PDXs). Each model was tested for in vivo sensitivity to three clinical regimens: cisplatin plus etoposide, olaparib plus temozolomide, and topotecan. These functional profiles captured hallmark clinical features, such as the emergence of treatment-refractory disease after early relapse. Serially derived PDX models from the same patient revealed that cross-resistance was acquired through a MYC amplification on extrachromosomal DNA (ecDNA). Genomic and transcriptional profiles of the full PDX panel revealed that this was not unique to one patient, as MYC paralog amplifications on ecDNAs were recurrent among cross-resistant models derived from patients after relapse. We conclude that ecDNAs with MYC paralogs are recurrent drivers of cross-resistance in SCLC.

SIGNIFICANCE

SCLC is initially chemosensitive, but acquired cross-resistance renders this disease refractory to further treatment and ultimately fatal. The genomic drivers of this transformation are unknown. We use a population of PDX models to discover that amplifications of MYC paralogs on ecDNA are recurrent drivers of acquired cross-resistance in SCLC.

Genome-wide characterization of extrachromosomal circular DNA in gastric cancer and its potential role in carcinogenesis and cancer progression

Extrachromosomal circular DNAs (eccDNAs) carrying random genomic segments are broadly found across different cancer types, but their molecular functions and impact in gastric cancer (GC) are rarely known. In this study, we aimed to investigate the potential role of eccDNA in GC. Using the Circle-seq strategy, we observed the eccDNA abundance in gastric cancer tissues (GCT) was aberrantly higher than that of normal adjacent tissues (NAT). The high abundance of eccDNAs carrying oncogene-segments in GCT may represent the DNA damage products of amplified oncogenes. Analysis of GCT over-represented eccDNA carrying enhancer (eccEnhancer) based on data from FANTOM5 project combined with TCGA database suggested the GC over-represented eccEnhancers may contribute to development of GC. GC over-represented eccDNAs carrying pre-miRNA (eccMIR) were enriched to multiple cancer-relevant signal pathways by KEGG analysis. We then synthesized the top six GC over-represented eccMIRs and found four of them enabled high expression of miRNAs and down-regulation of miRNA-target genes in MGC803 cells. Furthermore, we observed the inheritance of GC over-represented eccMIRs benefited host cell proliferation and promoted the aggressive features of host cells. Altogether, this study revealed the GC over-represented eccDNAs carrying functional genomic segments were related to the carcinogenesis of GC and presented the capability to facilitate cancer progression, suggesting the cancerous eccDNAs may serve as a dynamic reservoir for genome plasticity and rapid adaptive evolution of cancer. Therefore, blocking the pathways for eccDNAs generation may provide a novel therapeutic strategy for the treatment of gastric cancer.

Immortalization and transformation of primary cells mediated by engineered ecDNAs

Focal gene amplifications are among the most common cancer-associated mutations, but their evolution and contribution to tumorigenesis have proven challenging to recapitulate in primary cells and model organisms. Here we describe a general approach to engineer large (>1 Mbp) focal amplifications mediated by extrachromosomal circular DNAs (ecDNAs, also known as "double minutes") in a spatiotemporally controlled manner in cancer cell lines and in primary cells derived from genetically engineered mice. With this strategy, ecDNA formation can be coupled with expression of fluorescent reporters or other selectable markers to enable the identification and tracking of ecDNA-containing cells. We demonstrate the feasibility of this approach by engineering MDM2-containing ecDNAs in near-diploid human cells, showing that GFP expression can be used to track ecDNA dynamics under physiological conditions or in the presence of specific selective pressures. We also apply this approach to generate mice harboring inducible Myc - and Mdm2 -containing ecDNAs analogous to those spontaneously occurring in human cancers. We show that the engineered ecDNAs rapidly accumulate in primary cells derived from these animals, promoting proliferation, immortalization, and transformation.

Oncogenic structural aberration landscape in gastric cancer genomes

Structural variants (SVs) are responsible for driver events in gastric cancer (GC); however, their patterns and processes remain poorly understood. Here, we examine 170 GC whole genomes to unravel the oncogenic structural aberration landscape in GC genomes and identify six rearrangement signatures (RSs). Non-random combinations of RSs elucidate distinctive GC subtypes comprising one or a few dominant RS that are associated with specific driver events (BRCA1/2 defects, mismatch repair deficiency, and TP53 mutation) and epidemiological backgrounds. Twenty-seven SV hotspots are identified as GC driver candidates. SV hotspots frequently constitute complexly clustered SVs involved in driver gene amplification, such as ERBB2, CCNE1, and FGFR2. Further deconstruction of the locally clustered SVs uncovers amplicon-generating profiles characterized by super-large SVs and intensive segmental amplifications, contributing to the extensive amplification of GC oncogenes. Comprehensive analyses using adjusted SV allele frequencies indicate the significant involvement of extra-chromosomal DNA in processes linked to specific RSs.

Advancements in Focal Amplification Detection in Tumor/Liquid Biopsies and Emerging Clinical Applications

Focal amplifications (FAs) are crucial in cancer research due to their significant diagnostic, prognostic, and therapeutic implications. FAs manifest in various forms, such as episomes, double minute chromosomes, and homogeneously staining regions, arising through different mechanisms and mainly contributing to cancer cell heterogeneity, the leading cause of drug resistance in therapy. Numerous wet-lab, mainly FISH, PCR-based assays, next-generation sequencing, and bioinformatics approaches have been set up to detect FAs, unravel the internal structure of amplicons, assess their chromatin compaction status, and investigate the transcriptional landscape associated with their occurrence in cancer cells. Most of them are tailored for tumor samples, even at the single-cell level. Conversely, very limited approaches have been set up to detect FAs in liquid biopsies. This evidence suggests the need to improve these non-invasive investigations for early tumor detection, monitoring disease progression, and evaluating treatment response. Despite the potential therapeutic implications of FAs, such as, for example, the use of HER2-specific compounds for patients with ERBB2 amplification, challenges remain, including developing selective and effective FA-targeting agents and understanding the molecular mechanisms underlying FA maintenance and replication. This review details a state-of-the-art of FA investigation, with a particular focus on liquid biopsies and single-cell approaches in tumor samples, emphasizing their potential to revolutionize the future diagnosis, prognosis, and treatment of cancer patients.

Combined direct/indirect detection allows identification of DNA termini in diverse sequencing datasets and supports a multiple-initiation-site model for HIV plus-strand synthesis

Replication of genetic material involves the creation of characteristic termini. Determining these termini is important to refine our understanding of the mechanisms involved in maintaining the genomes of cellular organisms and viruses. Here we describe a computational approach combining direct and indirect readouts to detect termini from next-generation short-read sequencing. While a direct inference of termini can come from mapping the most prominent start positions of captured DNA fragments, this approach is insufficient in cases where the DNA termini are not captured, whether for biological or technical reasons. Thus, a complementary (indirect) approach to terminus detection can be applied, taking advantage of the imbalance in coverage between forward and reverse sequence reads near termini. A resulting metric ("strand bias") can be used to detect termini even where termini are naturally blocked from capture or ends are not captured during library preparation (e.g., in tagmentation-based protocols). Applying this analysis to datasets where known DNA termini are present, such as from linear double-stranded viral genomes, yielded distinct strand bias signals corresponding to these termini. To evaluate the potential to analyze a more complex situation, we applied the analysis to examine DNA termini present early after HIV infection in a cell culture model. We observed both the known termini expected based on standard models of HIV reverse transcription (the U5-right-end and U3-left-end termini) as well as a signal corresponding to a previously described additional initiation site for plus-strand synthesis (cPPT [central polypurine tract]). Interestingly, we also detected putative terminus signals at additional sites. The strongest of these are a set that share several characteristics with the previously characterized plus-strand initiation sites (the cPPT and 3' PPT [polypurine tract] sites): (i) an observed spike in directly captured cDNA ends, an indirect terminus signal evident in localized strand bias, (iii) a preference for location on the plus-strand, (iv) an upstream purine-rich motif, and (v) a decrease in terminus signal at late time points after infection. These characteristics are consistent in duplicate samples in two different genotypes (wild type and integrase-lacking HIV). The observation of distinct internal termini associated with multiple purine-rich regions raises a possibility that multiple internal initiations of plus-strand synthesis might contribute to HIV replication.

Circlehunter: a tool to identify extrachromosomal circular DNA from ATAC-Seq data

In cancer, extrachromosomal circular DNA (ecDNA), or megabase-pair amplified circular DNA, plays an essential role in intercellular heterogeneity and tumor cell revolution because of its non-Mendelian inheritance. We developed circlehunter ( https://github.com/suda-huanglab/circlehunter ), a tool for identifying ecDNA from ATAC-Seq data using the enhanced chromatin accessibility of ecDNA. Using simulated data, we showed that circlehunter has an F1 score of 0.93 at 30× local depth and read lengths as short as 35 bp. Based on 1312 ecDNAs predicted from 94 publicly available datasets of ATAC-Seq assays, we found 37 oncogenes contained in these ecDNAs with amplification characteristics. In small cell lung cancer cell lines, ecDNA containing MYC leads to amplification of MYC and cis-regulates the expression of NEUROD1, resulting in an expression pattern consistent with the NEUROD1 high expression subtype and sensitive to Aurora kinase inhibitors. This showcases that circlehunter could serve as a valuable pipeline for the investigation of tumorigenesis.

Advances in sequencing-based studies of microDNA and ecDNA: Databases, identification methods, and integration with single-cell analysis

Extrachromosomal circular DNA (eccDNA) is a class of circular DNA molecules that originate from genomic DNA but are separate from chromosomes. They are common in various organisms, with sizes ranging from a few hundred to millions of base pairs. A special type of large extrachromosomal DNA (ecDNA) is prevalent in cancer cells. Research on ecDNA has significantly contributed to our comprehension of cancer development, progression, evolution, and drug resistance. The use of next-generation (NGS) and third-generation sequencing (TGS) techniques to identify eccDNAs throughout the genome has become a trend in current research. Here, we briefly review current advances in the biological mechanisms and applications of two distinct types of eccDNAs: microDNA and ecDNA. In addition to presenting available identification tools based on sequencing data, we summarize the most recent efforts to integrate ecDNA with single-cell analysis and put forth suggestions to promote the process.

Integrated analysis of single-cell chromatin state and transcriptome identified common vulnerability despite glioblastoma heterogeneity

In 2021, the World Health Organization reclassified glioblastoma, the most common form of adult brain cancer, into isocitrate dehydrogenase (IDH)-wild-type glioblastomas and grade IV IDH mutant (G4 IDHm) astrocytomas. For both tumor types, intratumoral heterogeneity is a key contributor to therapeutic failure. To better define this heterogeneity, genome-wide chromatin accessibility and transcription profiles of clinical samples of glioblastomas and G4 IDHm astrocytomas were analyzed at single-cell resolution. These profiles afforded resolution of intratumoral genetic heterogeneity, including delineation of cell-to-cell variations in distinct cell states, focal gene amplifications, as well as extrachromosomal circular DNAs. Despite differences in IDH mutation status and significant intratumoral heterogeneity, the profiled tumor cells shared a common chromatin structure defined by open regions enriched for nuclear factor 1 transcription factors (NFIA and NFIB). Silencing of NFIA or NFIB suppressed in vitro and in vivo growths of patient-derived glioblastomas and G4 IDHm astrocytoma models. These findings suggest that despite distinct genotypes and cell states, glioblastoma/G4 astrocytoma cells share dependency on core transcriptional programs, yielding an attractive platform for addressing therapeutic challenges associated with intratumoral heterogeneity.

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.

Multi-organ landscape of therapy-resistant melanoma

Metastasis and failure of present-day therapies represent the most common causes of mortality in patients with cutaneous melanoma. To identify the underlying genetic and transcriptomic landscapes, in this study we analyzed multi-organ metastases and tumor-adjacent tissues from 11 rapid autopsies after treatment with MAPK inhibitor (MAPKi) and/or immune checkpoint blockade (ICB) and death due to acquired resistance. Either treatment elicits shared genetic alterations that suggest immune-evasive, cross-therapy resistance mechanisms. Large, non-clustered deletions, inversions and inter-chromosomal translocations dominate rearrangements. Analyzing data from separate melanoma cohorts including 345 therapy-naive patients and 35 patients with patient-matched pre-treatment and post-acquired resistance tumor samples, we performed cross-cohort analyses to identify MAPKi and ICB as respective contributors to gene amplifications and deletions enriched in autopsy versus therapy-naive tumors. In the autopsy cohort, private/late mutations and structural variants display shifted mutational and rearrangement signatures, with MAPKi specifically selecting for signatures of defective homologous-recombination, mismatch and base-excision repair. Transcriptomic signatures and crosstalks with tumor-adjacent macroenvironments nominated organ-specific adaptive pathways. An immune-desert, CD8+-macrophage-biased archetype, T-cell exhaustion and type-2 immunity characterized the immune contexture. This multi-organ analysis of therapy-resistant melanoma presents preliminary insights with potential to improve therapeutic strategies.

Blocking Genomic Instability Prevents Acquired Resistance to MAPK Inhibitor Therapy in Melanoma

Blocking cancer genomic instability may prevent tumor diversification and escape from therapies. We show that, after MAPK inhibitor (MAPKi) therapy in patients and mice bearing patient-derived xenografts (PDX), acquired resistant genomes of metastatic cutaneous melanoma specifically amplify resistance-driver, nonhomologous end-joining (NHEJ), and homologous recombination repair (HRR) genes via complex genomic rearrangements (CGR) and extrachromosomal DNAs (ecDNA). Almost all sensitive and acquired-resistant genomes harbor pervasive chromothriptic regions with disproportionately high mutational burdens and significant overlaps with ecDNA and CGR spans. Recurrently, somatic mutations within ecDNA and CGR amplicons enrich for HRR signatures, particularly within acquired resistant tumors. Regardless of sensitivity or resistance, breakpoint-junctional sequence analysis suggests NHEJ as critical to double-stranded DNA break repair underlying CGR and ecDNA formation. In human melanoma cell lines and PDXs, NHEJ targeting by a DNA-PKCS inhibitor prevents/delays acquired MAPKi resistance by reducing the size of ecDNAs and CGRs early on combination treatment. Thus, targeting the causes of genomic instability prevents acquired resistance.

SIGNIFICANCE

Acquired resistance often results in heterogeneous, redundant survival mechanisms, which challenge strategies aimed at reversing resistance. Acquired-resistant melanomas recurrently evolve resistance-driving and resistance-specific amplicons via ecDNAs and CGRs, thereby nominating chromothripsis-ecDNA-CGR biogenesis as a resistance-preventive target. Specifically, targeting DNA-PKCS/NHEJ prevents resistance by suppressing ecDNA/CGR rearrangements in MAPKi-treated melanomas. This article is highlighted in the In This Issue feature, p. 799.