Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

The genetic evolution of acral melanoma

Acral melanoma is an aggressive type of melanoma with unknown origins. It is the most common type of melanoma in individuals with dark skin and is notoriously challenging to treat. We examine exome sequencing data of 139 tissue samples, spanning different progression stages, from 37 patients. We find that 78.4% of the melanomas display clustered copy number transitions with focal amplifications, recurring predominantly on chromosomes 5, 11, 12, and 22. These complex genomic aberrations are typically shared across all progression stages of individual patients. TERT activating alterations also arise early, whereas MAP-kinase pathway mutations appear later, an inverted order compared to the canonical evolution. The punctuated formation of complex aberrations and early TERT activation suggest a unique mutational mechanism that initiates acral melanoma. The marked intratumoral heterogeneity, especially concerning MAP-kinase pathway mutations, may partly explain the limited success of therapies for this melanoma subtype.

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.

Integrated analyses highlight interactions between the three-dimensional genome and DNA, RNA and epigenomic alterations in metastatic prostate cancer

The impact of variations in the three-dimensional structure of the genome has been recognized, but solid cancer tissue studies are limited. Here, we performed integrated deep Hi-C sequencing with matched whole-genome sequencing, whole-genome bisulfite sequencing, 5-hydroxymethylcytosine (5hmC) sequencing and RNA sequencing across a cohort of 80 biopsy samples from patients with metastatic castration-resistant prostate cancer. Dramatic differences were present in gene expression, 5-methylcytosine/5hmC methylation and in structural variation versus mutation rate between A and B (open and closed) chromatin compartments. A subset of tumors exhibited depleted regional chromatin contacts at the AR locus, linked to extrachromosomal circular DNA (ecDNA) and worse response to AR signaling inhibitors. We also identified topological subtypes associated with stark differences in methylation structure, gene expression and prognosis. Our data suggested that DNA interactions may predispose to structural variant formation, exemplified by the recurrent TMPRSS2-ERG fusion. This comprehensive integrated sequencing effort represents a unique clinical tumor resource.

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.

Transcriptional immune suppression and up-regulation of double-stranded DNA damage and repair repertoires in ecDNA-containing tumors

Extrachromosomal DNA is a common cause of oncogene amplification in cancer. The non-chromosomal inheritance of ecDNA enables tumors to rapidly evolve, contributing to treatment resistance and poor outcome for patients. The transcriptional context in which ecDNAs arise and progress, including chromosomally-driven transcription, is incompletely understood. We examined gene expression patterns of 870 tumors of varied histological types, to identify transcriptional correlates of ecDNA. Here, we show that ecDNA-containing tumors impact four major biological processes. Specifically, ecDNA-containing tumors up-regulate DNA damage and repair, cell cycle control, and mitotic processes, but down-regulate global immune regulation pathways. Taken together, these results suggest profound alterations in gene regulation in ecDNA-containing tumors, shedding light on molecular processes that give rise to their development and progression.

Quantifying spatial CXCL9 distribution with image analysis predicts improved prognosis of triple-negative breast cancer

Background: The C-X-C motif chemokine ligand 9 (CXCL9) plays a pivotal role in tumor immunity by recruiting and activating immune cells. However, the relationship between CXCL9 expression and prognosis in triple-negative breast cancer (TNBC) is unclear. Methods: We investigated CXCL9 mRNA expression, clinicopathological features, and prognosis in TNBC patients. We also used computational image analysis to quantify and assess the distribution of CXCL9 protein in the tumor core (TC) and invasive margin (IM). Results: CXCL9 mRNA expression was significantly higher in TNBC tumors compared to normal tissue (p < 0.001) and was associated with smaller tumors (p = 0.022) and earlier stages (p = 0.033). High CXCL9 mRNA expression was correlated with improved overall survival (OS) in three independent cohorts (all p < 0.05). In a separate analysis, low CXCL9 protein expression was associated with increased lymph node metastasis (p = 0.018 and p = 0.036). High CXCL9 protein expression in the TC, IM, or both was associated with prolonged OS (all p < 0.001). Conclusion: High CXCL9 expression, at both the mRNA and protein levels, is associated with improved prognosis in TNBC patients. CXCL9 expression in the TC and/or IM may be an independent prognostic factor.

Receptor tyrosine kinase inhibition leads to regression of acral melanoma by targeting the tumor microenvironment

Acral melanoma (AM) is an aggressive melanoma variant that arises from palmar, plantar, and nail unit melanocytes. Compared to non-acral cutaneous melanoma (CM), AM is biologically distinct, has an equal incidence across genetic ancestries, typically presents in advanced stage disease, is less responsive to therapy, and has an overall worse prognosis. Independent analysis of published genomic and transcriptomic sequencing identified that receptor tyrosine kinase (RTK) ligands and adapter proteins are frequently amplified, translocated, and/or overexpressed in AM. To target these unique genetic changes, a zebrafish acral melanoma model was exposed to a panel of narrow and broad spectrum multi-RTK inhibitors, revealing that dual FGFR/VEGFR inhibitors decrease acral-analogous melanocyte proliferation and migration. The potent pan-FGFR/VEGFR inhibitor, Lenvatinib, uniformly induces tumor regression in AM patient-derived xenograft (PDX) tumors but only slows tumor growth in CM models. Unlike other multi-RTK inhibitors, Lenvatinib is not directly cytotoxic to dissociated AM PDX tumor cells and instead disrupts tumor architecture and vascular networks. Considering the great difficulty in establishing AM cell culture lines, these findings suggest that AM may be more sensitive to microenvironment perturbations than CM. In conclusion, dual FGFR/VEGFR inhibition may be a viable therapeutic strategy that targets the unique biology of AM.

CoRAL accurately resolves extrachromosomal DNA genome structures with long-read sequencing

Extrachromosomal DNA (ecDNA) is a central mechanism for focal oncogene amplification in cancer, occurring in approximately 15% of early stage cancers and 30% of late-stage cancers. EcDNAs drive tumor formation, evolution, and drug resistance by dynamically modulating oncogene copy-number and rewiring gene-regulatory networks. Elucidating the genomic architecture of ecDNA amplifications is critical for understanding tumor pathology and developing more effective therapies. Paired-end short-read (Illumina) sequencing and mapping have been utilized to represent ecDNA amplifications using a breakpoint graph, where the inferred architecture of ecDNA is encoded as a cycle in the graph. Traversals of breakpoint graph have been used to successfully predict ecDNA presence in cancer samples. However, short-read technologies are intrinsically limited in the identification of breakpoints, phasing together of complex rearrangements and internal duplications, and deconvolution of cell-to-cell heterogeneity of ecDNA structures. Long-read technologies, such as from Oxford Nanopore Technologies, have the potential to improve inference as the longer reads are better at mapping structural variants and are more likely to span rearranged or duplicated regions. Here, we propose CoRAL (Complete Reconstruction of Amplifications with Long reads), for reconstructing ecDNA architectures using long-read data. CoRAL reconstructs likely cyclic architectures using quadratic programming that simultaneously optimizes parsimony of reconstruction, explained copy number, and consistency of long-read mapping. CoRAL substantially improves reconstructions in extensive simulations and 9 datasets from previously-characterized cell-lines as compared to previous short-read-based tools. As long-read usage becomes wide-spread, we anticipate that CoRAL will be a valuable tool for profiling the landscape and evolution of focal amplifications in tumors.

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.

Chromoplexy Is a Frequent Early Clonal Event in EWSR1-Rearranged Round Cell Sarcomas That Can Be Detected Using Clinically Validated Targeted Sequencing Panels

Chromoplexy is a phenomenon defined by large-scale chromosomal chained rearrangements. A previous study observed chromoplectic events in a subset of Ewing sarcomas (ES), which was linked to an increased relapse rate. Chromoplexy analysis could potentially facilitate patient risk stratification, particularly if it could be detected with clinically applied targeted next-generation sequencing (NGS) panels. Using DELLY, a structural variant (SV) calling algorithm that is part of the MSK-IMPACT pipeline, we characterized the spectrum of SVs in EWSR1-fused round cell sarcomas, including 173 ES and 104 desmoplastic small round cell tumors (DSRCT), to detect chromoplexy and evaluate its association with clinical and genomic features. Chromoplectic events were detected in 31% of the ES cases and 19% of the DSRCT cases. EWSR1 involvement accounted for 76% to 93% of these events, being rearranged with diverse noncanonical gene partners across the genome, involving mainly translocations but also intrachromosomal deletions and inversions. A major breakpoint cluster was located on EWSR1 exons 8-13. In a subset of cases, the SVs disrupted adjacent loci, forming deletion bridges. Longitudinal sequencing and breakpoint allele fraction analysis showed that chromoplexy is an early event that remains detectable throughout disease progression and likely develops simultaneously with the driver fusion. The presence of chromoplexy was validated in an external ES patient cohort with whole exome sequencing. Chromoplexy was significantly more likely to be present in cases that were metastatic at presentation. Together, this study identifies chromoplexy as a frequent genomic alteration in diverse EWSR1-rearranged tumors that can be captured by targeted NGS panels.

SIGNIFICANCE

Chromoplexy is detectable using targeted NGS in a substantial portion of EWSR1-rearranged round cell sarcomas as an early and persistent clonal event, expanding the genomic complexity of fusion-associated sarcomas.

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.

Severus: accurate detection and characterization of somatic structural variation in tumor genomes using long reads

Most current studies rely on short-read sequencing to detect somatic structural variation (SV) in cancer genomes. Long-read sequencing offers the advantage of better mappability and long-range phasing, which results in substantial improvements in germline SV detection. However, current long-read SV detection methods do not generalize well to the analysis of somatic SVs in tumor genomes with complex rearrangements, heterogeneity, and aneuploidy. Here, we present Severus: a method for the accurate detection of different types of somatic SVs using a phased breakpoint graph approach. To benchmark various short- and long-read SV detection methods, we sequenced five tumor/normal cell line pairs with Illumina, Nanopore, and PacBio sequencing platforms; on this benchmark Severus showed the highest F1 scores (harmonic mean of the precision and recall) as compared to long-read and short-read methods. We then applied Severus to three clinical cases of pediatric cancer, demonstrating concordance with known genetic findings as well as revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

Childhood cancer mutagenesis caused by transposase-derived PGBD5

Genomic rearrangements are a hallmark of most childhood tumors, including medulloblastoma, one of the most common brain tumors in children, but their causes remain largely unknown. Here, we show that PiggyBac transposable element derived 5 (Pgbd5) promotes tumor development in multiple developmentally accurate mouse models of Sonic Hedgehog (SHH) medulloblastoma. Most Pgbd5-deficient mice do not develop tumors, while maintaining normal cerebellar development. Ectopic activation of SHH signaling is sufficient to enforce cerebellar granule cell progenitor-like cell states, which exhibit Pgbd5-dependent expression of distinct DNA repair and neurodevelopmental factors. Mouse medulloblastomas expressing Pgbd5 have increased numbers of somatic structural DNA rearrangements, some of which carry PGBD5-specific sequences at their breakpoints. Similar sequence breakpoints recurrently affect somatic DNA rearrangements of known tumor suppressors and oncogenes in medulloblastomas in 329 children. This identifies PGBD5 as a medulloblastoma mutator and provides a genetic mechanism for the generation of oncogenic DNA rearrangements in childhood cancer.

Whole-genome doubling is a double-edged sword: the heterogeneous role of whole-genome doubling in various cancer types

Whole-genome doubling (WGD), characterized by the duplication of an entire set of chromosomes, is commonly observed in various tumors, occurring in approximately 30-40% of patients with different cancer types. The effect of WGD on tumorigenesis varies depending on the context, either promoting or suppressing tumor progression. Recent advances in genomic technologies and large-scale clinical investigations have led to the identification of the complex patterns of genomic alterations underlying WGD and their functional consequences on tumorigenesis progression and prognosis. Our comprehensive review aims to summarize the causes and effects of WGD on tumorigenesis, highlighting its dualistic influence on cancer cells. We then introduce recent findings on WGD-associated molecular signatures and genetic aberrations and a novel subtype related to WGD. Finally, we discuss the clinical implications of WGD in cancer subtype classification and future therapeutic interventions. Overall, a comprehensive understanding of WGD in cancer biology is crucial to unraveling its complex role in tumorigenesis and identifying novel therapeutic strategies. [BMB Reports 2024; 57(3): 125-134].

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.

ATR blocks telomerase from converting DNA breaks into telomeres

Telomerase, the enzyme that maintains telomeres at natural chromosome ends, should be repressed at double-strand breaks (DSBs), where neotelomere formation can cause terminal truncations. We developed an assay to detect neotelomere formation at Cas9- or I-SceI-induced DSBs in human cells. Telomerase added telomeric repeats to DSBs, leading to interstitial telomeric repeat insertions or the formation of functional neotelomeres accompanied by terminal deletions. The threat that telomerase poses to genome integrity was minimized by ataxia telangiectasia and Rad3-related (ATR) kinase signaling, which inhibited telomerase at resected DSBs. In addition to acting at resected DSBs, telomerase used the extruded strand in the Cas9 enzyme-product complex as a primer for neotelomere formation. We propose that although neotelomere formation is detrimental in normal human cells, it may allow cancer cells to escape from breakage-fusion-bridge cycles.

Aneuploidy and complex genomic rearrangements in cancer evolution

Mutational processes that alter large genomic regions occur frequently in developing tumors. They range from simple copy number gains and losses to the shattering and reassembly of entire chromosomes. These catastrophic events, such as chromothripsis, chromoplexy and the formation of extrachromosomal DNA, affect the expression of many genes and therefore have a substantial effect on the fitness of the cells in which they arise. In this review, we cover large genomic alterations, the mechanisms that cause them and their effect on tumor development and evolution.

Urine-derived bladder cancer organoids (urinoids) as a tool for cancer longitudinal response monitoring and therapy adaptation

BACKGROUND

Bladder cancer is one of the most common cancer types worldwide. Generally, research relies on invasive sampling strategies.

METHODS

Here, we generate bladder cancer organoids directly from urine (urinoids). In this project, we establish 12 urinoid lines from 22 patients with non-muscle and muscle-invasive bladder tumours, with an efficiency of 55%.

RESULTS

The histopathological features of the urinoids accurately resemble those of the original bladder tumours. Genetically, there is a high concordance of single nucleotide polymorphisms (92.56%) and insertions & deletions (91.54%) between urinoids and original tumours from patient 4. Furthermore, these urinoids show sensitivity to bladder cancer drugs, similar to their tissue-derived organoid counterparts. Genetic analysis of longitudinally generated tumoroids and urinoids from one patient receiving systemic immunotherapy, identify alterations that may guide the choice for second-line therapy. Successful treatment adaptation was subsequently demonstrated in the urinoid setting.

CONCLUSION

Therefore, urinoids can advance precision medicine in bladder cancer as a non-invasive platform for tumour pathogenesis, longitudinal drug-response monitoring, and therapy adaptation.

LRTK: a platform agnostic toolkit for linked-read analysis of both human genome and metagenome

BACKGROUND

Linked-read sequencing technologies generate high-base quality short reads that contain extrapolative information on long-range DNA connectedness. These advantages of linked-read technologies are well known and have been demonstrated in many human genomic and metagenomic studies. However, existing linked-read analysis pipelines (e.g., Long Ranger) were primarily developed to process sequencing data from the human genome and are not suited for analyzing metagenomic sequencing data. Moreover, linked-read analysis pipelines are typically limited to 1 specific sequencing platform.

FINDINGS

To address these limitations, we present the Linked-Read ToolKit (LRTK), a unified and versatile toolkit for platform agnostic processing of linked-read sequencing data from both human genome and metagenome. LRTK provides functions to perform linked-read simulation, barcode sequencing error correction, barcode-aware read alignment and metagenome assembly, reconstruction of long DNA fragments, taxonomic classification and quantification, and barcode-assisted genomic variant calling and phasing. LRTK has the ability to process multiple samples automatically and provides users with the option to generate reproducible reports during processing of raw sequencing data and at multiple checkpoints throughout downstream analysis. We applied LRTK on linked reads from simulation, mock community, and real datasets for both human genome and metagenome. We showcased LRTK's ability to generate comparative performance results from preceding benchmark studies and to report these results in publication-ready HTML document plots.

CONCLUSIONS

LRTK provides comprehensive and flexible modules along with an easy-to-use Python-based workflow for processing linked-read sequencing datasets, thereby filling the current gap in the field caused by platform-centric genome-specific linked-read data analysis tools.

The importance of escalating molecular diagnostics in patients with low-grade pediatric brain cancer

Pilocytic astrocytomas are the most common pediatric brain tumors, typically presenting as low-grade neoplasms. We report two cases of pilocytic astrocytoma with atypical tumor progression. Case 1 involves a 12-yr-old boy with an unresectable suprasellar tumor, negative for BRAF rearrangement but harboring a BRAF p.V600E mutation. He experienced tumor size reduction and stable disease following dabrafenib treatment. Case 2 describes a 6-yr-old boy with a thalamic tumor that underwent multiple resections, with no actionable driver detected using targeted next-generation sequencing. Whole-genome and RNA-seq analysis identified an internal tandem duplication in FGFR1 and RAS pathway activation. Future management options include FGFR1 inhibitors. These cases demonstrate the importance of escalating molecular diagnostics for pediatric brain cancer, advocating for early reflexing to integrative whole-genome sequencing and transcriptomic profiling when targeted panels are uninformative. Identifying molecular drivers can significantly impact treatment decisions and improve patient outcomes.

Most large structural variants in cancer genomes can be detected without long reads

Short-read sequencing is the workhorse of cancer genomics yet is thought to miss many structural variants (SVs), particularly large chromosomal alterations. To characterize missing SVs in short-read whole genomes, we analyzed 'loose ends'-local violations of mass balance between adjacent DNA segments. In the landscape of loose ends across 1,330 high-purity cancer whole genomes, most large (>10-kb) clonal SVs were fully resolved by short reads in the 87% of the human genome where copy number could be reliably measured. Some loose ends represent neotelomeres, which we propose as a hallmark of the alternative lengthening of telomeres phenotype. These pan-cancer findings were confirmed by long-molecule profiles of 38 breast cancer and melanoma cases. Our results indicate that aberrant homologous recombination is unlikely to drive the majority of large cancer SVs. Furthermore, analysis of mass balance in short-read whole genome data provides a surprisingly complete picture of cancer chromosomal structure.

Acral melanoma: new insights into the immune and genomic landscape

Acral melanoma is a rare subtype of melanoma that arises on the non-hair bearing skin of the nail bed, palms of the hand and soles of the feet. It is unique among melanomas in not being linked to ultraviolet radiation (UVR) exposure from the sun, and, as such, its incidence is similar across populations who are of Asian, Hispanic, African and European origin. Although research into acral melanoma has lagged behind that of sun-exposed cutaneous melanoma, recent studies have begun to address the unique genetics and immune features of acral melanoma. In this review we will discuss the latest progress in understanding the biology of acral melanoma across different ethnic populations and will outline how these new discoveries can help to guide the therapeutic management of this rare tumor.

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.

Two High-Quality Cygnus Genome Assemblies Reveal Genomic Variations Associated with Plumage Color

As an exemplary model for examining molecular mechanisms responsible for extreme phenotypic variations, plumage color has garnered significant interest. The Cygnus genus features two species, Cygnus olor and Cygnus atratus, that exhibit striking disparities in plumage color. However, the molecular foundation for this differentiation has remained elusive. Herein, we present two high-quality genomes for C. olor and C. atratus, procured using the Illumina and Nanopore technologies. The assembled genome of C. olor was 1.12 Gb in size with a contig N50 of 26.82 Mb, while its counterpart was 1.13 Gb in size with a contig N50 of 21.91 Mb. A comparative analysis unveiled three genes (TYR, SLC45A2, and SLC7A11) with structural variants in the melanogenic pathway. Notably, we also identified a novel gene, PWWP domain containing 2A (PWWP2A), that is related to plumage color, for the first time. Using targeted gene modification analysis, we demonstrated the potential genetic effect of the PWWP2A variant on pigment gene expression and melanin production. Finally, our findings offer insight into the intricate pattern of pigmentation and the role of polygenes in birds. Furthermore, these two high-quality genome references provide a comprehensive resource and perspective for comparative functional and genetic studies of evolution within the Cygnus genus.

Transcriptional immune suppression and upregulation of double stranded DNA damage and repair repertoires in ecDNA-containing tumors

Extrachromosomal DNA is a common cause of oncogene amplification in cancer. The non-chromosomal inheritance of ecDNA enables tumors to rapidly evolve, contributing to treatment resistance and poor outcome for patients. The transcriptional context in which ecDNAs arise and progress, including chromosomally-driven transcription, is incompletely understood. We examined gene expression patterns of 870 tumors of varied histological types, to identify transcriptional correlates of ecDNA. Here we show that ecDNA containing tumors impact four major biological processes. Specifically, ecDNA containing tumors upregulate DNA damage and repair, cell cycle control, and mitotic processes, but downregulate global immune regulation pathways. Taken together, these results suggest profound alterations in gene regulation in ecDNA containing tumors, shedding light on molecular processes that give rise to their development and progression.

Subclonal Somatic Copy-Number Alterations Emerge and Dominate in Recurrent Osteosarcoma

Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease.

SIGNIFICANCE

The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.

Improving Hi-C contact matrices using genome graphs

Three-dimensional chromosome structure plays an important role in fundamental genomic functions. Hi-C, a high-throughput, sequencing-based technique, has drastically expanded our comprehension of 3D chromosome structures. The first step of Hi-C analysis pipeline involves mapping sequencing reads from Hi-C to linear reference genomes. However, the linear reference genome does not incorporate genetic variation information, which can lead to incorrect read alignments, especially when analyzing samples with substantial genomic differences from the reference such as cancer samples. Using genome graphs as the reference facilitates more accurate mapping of reads, however, new algorithms are required for inferring linear genomes from Hi-C reads mapped on genome graphs and constructing corresponding Hi-C contact matrices, which is a prerequisite for the subsequent steps of the Hi-C analysis such as identifying topologically associated domains and calling chromatin loops. We introduce the problem of genome sequence inference from Hi-C data mediated by genome graphs. We formalize this problem, show the hardness of solving this problem, and introduce a novel heuristic algorithm specifically tailored to this problem. We provide a theoretical analysis to evaluate the efficacy of our algorithm. Finally, our empirical experiments indicate that the linear genomes inferred from our method lead to the creation of improved Hi-C contact matrices. These enhanced matrices show a reduction in erroneous patterns caused by structural variations and are more effective in accurately capturing the structures of topologically associated domains.

Tracing cancer evolution and heterogeneity using Hi-C

Chromosomal rearrangements can initiate and drive cancer progression, yet it has been challenging to evaluate their impact, especially in genetically heterogeneous solid cancers. To address this problem we developed HiDENSEC, a new computational framework for analyzing chromatin conformation capture in heterogeneous samples that can infer somatic copy number alterations, characterize large-scale chromosomal rearrangements, and estimate cancer cell fractions. After validating HiDENSEC with in silico and in vitro controls, we used it to characterize chromosome-scale evolution during melanoma progression in formalin-fixed tumor samples from three patients. The resulting comprehensive annotation of the genomic events includes copy number neutral translocations that disrupt tumor suppressor genes such as NF1, whole chromosome arm exchanges that result in loss of CDKN2A, and whole-arm copy-number neutral loss of homozygosity involving PTEN. These findings show that large-scale chromosomal rearrangements occur throughout cancer evolution and that characterizing these events yields insights into drivers of melanoma progression.

ChArmTelo Enables Large-Scale Chromosome Arm-Level Telomere Analysis across Human Populations and in Cancer Patients

Telomeres are structures protecting chromosome ends. However, a scalable and cost-effective method to investigate chromosome arm-level (ChArm) telomeres (Telos) in large-scale projects is still lacking, hindering intensive investigation of high-resolution telomeres across cancers and other diseases. Here, ChArmTelo, the first computational toolbox to analyze telomeres at chromosome arm level in human and other animal species, using 10X linked-read and similar technologies, is presented. ChArmTelo currently consists of two algorithms, TeloEM and TeloKnow, for arm-level telomere length (TL) analysis. The algorithms are demonstrated by comprehensive analysis of chromosome arm-level telomere lengths (chArmTLs) in nearly 400 whole genome sequencing samples (WGS) from human populations and animals, including healthy and cancer samples. Notably, considerable performance improvement contributed by using the latest complete telomere-to-telomere reference genome (CHM13v2), compared to hg38, is shown. ChArmTelo reveals population-specific chArmTL differences and liver cancer signatures of chArmTLs and that DNA replication origin disruption may contribute to cancer by affecting TLs. Importantly, ChArmTelo can be readily applied to tens of thousands of cancer and healthy samples with published WGS data.

A pan-cancer landscape of pathogenic somatic copy number variations

OBJECTIVE

Copy number variations (CNVs) play crucial roles in physiological and pathological processes, including cancer. However, the functional implications of somatic CNVs in tumor progression and evolution remain unclear. This study focuses on identifying CNV alterations with high pathogenic potential that drive and sustain tumorigenesis, distinguishing them from passenger alterations that accumulate during tumor growth. Our goal is to explore the variability of CNVs across different tumor types and infer their impact on tumor cell functions.

METHODS

Starting from 7352 copy number profiles across 33 different cancer types, we infer the pathogenicity of each CNV and perform both intra- and inter-tumor analyses to predict the functional impact of different genomic patterns. We evaluate the actionability of genes belonging to altered regions and we correlate the presence of pathogenic regions with genome instability patterns and patients' survival.

RESULTS

Our analysis uncovered large heterogeneity among different tumors suggesting in many cases distinct genetic drivers of tumorigenesis. Recurrent genomic alterations frequently coincide with dysfunctional homologous recombination pathways and negative regulation of the immune system. In certain tumors, the number of pathogenic CNVs emerged as a prognostic biomarker, highlighting their significance in cancer progression.

CONCLUSION

This study contributes to elucidate the functional impact of pathogenic CNVs in tumor progression and sheds light on their potential as prognostic markers in specific cancer types.

Scientific techniques in adolescent and young adult classic Hodgkin lymphoma

Understanding the tumor microenvironment and genomic landscape is crucial for better prediction of treatment outcomes and developing novel therapies in Hodgkin lymphoma (HL). Recent advancements in genomics have enabled researchers to gain deeper insights into the genomic characteristics of HL at both single-cell resolution and the whole genome level. The use of noninvasive methods such as liquid biopsies and formalin-fixed paraffin-embedded-based imaging techniques has expanded the possibilities of applying cutting-edge analyses to routine clinically available samples. Collaborative efforts between adult and pediatric group are imperative to translate novel findings into routine patient care.

The genetic evolution of acral melanoma

Acral melanoma is an aggressive type of melanoma with unknown origins, arising on the sole, palm, or nail apparatus. It is the most common type of melanoma in individuals with dark skin and is notoriously challenging to treat. Our study examined exome sequencing data from 139 tissue samples, spanning different progression stages, collected from 37 patients. We found that 78.4% of the melanomas displayed one or more clustered copy number transitions with focal amplifications, recurring predominantly on chromosomes 5, 11, 12, and 22. These genomic "hailstorms" were typically shared across all progression stages within individual patients. Genetic alterations known to activate TERT also arose early. By contrast, mutations in the MAP-kinase pathway appeared later during progression, often leading to different tumor areas harboring non-overlapping driver mutations. We conclude that the evolutionary trajectories of acral melanomas substantially diverge from those of melanomas on sun-exposed skin, where MAP-kinase pathway activation initiates the neoplastic cascade followed by immortalization later. The punctuated formation of hailstorms, paired with early TERT activation, suggests a unique mutational mechanism underlying the origins of acral melanoma. Our findings highlight an essential role for telomerase, likely in re-stabilizing tumor genomes after hailstorms have initiated the tumors. The marked genetic heterogeneity, in particular of MAP-kinase pathway drivers, may partly explain the limited success of targeted and other therapies in treating this melanoma subtype.

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.

Decoding Cancer Evolution: Integrating Genetic and Non-Genetic Insights

The development of cancer begins with cells transitioning from their multicellular nature to a state akin to unicellular organisms. This shift leads to a breakdown in the crucial regulators inherent to multicellularity, resulting in the emergence of diverse cancer cell subpopulations that have enhanced adaptability. The presence of different cell subpopulations within a tumour, known as intratumoural heterogeneity (ITH), poses challenges for cancer treatment. In this review, we delve into the dynamics of the shift from multicellularity to unicellularity during cancer onset and progression. We highlight the role of genetic and non-genetic factors, as well as tumour microenvironment, in promoting ITH and cancer evolution. Additionally, we shed light on the latest advancements in omics technologies that allow for in-depth analysis of tumours at the single-cell level and their spatial organization within the tissue. Obtaining such detailed information is crucial for deepening our understanding of the diverse evolutionary paths of cancer, allowing for the development of effective therapies targeting the key drivers of cancer evolution.

Evolution of structural rearrangements in prostate cancer intracranial metastases

Intracranial metastases in prostate cancer are uncommon but clinically aggressive. A detailed molecular characterization of prostate cancer intracranial metastases would improve our understanding of their pathogenesis and the search for new treatment strategies. We evaluated the clinical and molecular characteristics of 36 patients with metastatic prostate cancer to either the dura or brain parenchyma. We performed whole genome sequencing (WGS) of 10 intracranial prostate cancer metastases, as well as WGS of primary prostate tumors from men who later developed metastatic disease (n = 6) and nonbrain prostate cancer metastases (n = 36). This first whole genome sequencing study of prostate intracranial metastases led to several new insights. First, there was a higher diversity of complex structural alterations in prostate cancer intracranial metastases compared to primary tumor tissues. Chromothripsis and chromoplexy events seemed to dominate, yet there were few enrichments of specific categories of structural variants compared with non-brain metastases. Second, aberrations involving the AR gene, including AR enhancer gain were observed in 7/10 (70%) of intracranial metastases, as well as recurrent loss of function aberrations involving TP53 in 8/10 (80%), RB1 in 2/10 (20%), BRCA2 in 2/10 (20%), and activation of the PI3K/AKT/PTEN pathway in 8/10 (80%). These alterations were frequently present in tumor tissues from other sites of disease obtained concurrently or sequentially from the same individuals. Third, clonality analysis points to genomic factors and evolutionary bottlenecks that contribute to metastatic spread in patients with prostate cancer. These results describe the aggressive molecular features underlying intracranial metastasis that may inform future diagnostic and treatment approaches.

Long-molecule scars of backup DNA repair in BRCA1- and BRCA2-deficient cancers

Homologous recombination (HR) deficiency is associated with DNA rearrangements and cytogenetic aberrations1. Paradoxically, the types of DNA rearrangements that are specifically associated with HR-deficient cancers only minimally affect chromosomal structure2. Here, to address this apparent contradiction, we combined genome-graph analysis of short-read whole-genome sequencing (WGS) profiles across thousands of tumours with deep linked-read WGS of 46 BRCA1- or BRCA2-mutant breast cancers. These data revealed a distinct class of HR-deficiency-enriched rearrangements called reciprocal pairs. Linked-read WGS showed that reciprocal pairs with identical rearrangement orientations gave rise to one of two distinct chromosomal outcomes, distinguishable only with long-molecule data. Whereas one (cis) outcome corresponded to the copying and pasting of a small segment to a distant site, a second (trans) outcome was a quasi-balanced translocation or multi-megabase inversion with substantial (10 kb) duplications at each junction. We propose an HR-independent replication-restart repair mechanism to explain the full spectrum of reciprocal pair outcomes. Linked-read WGS also identified single-strand annealing as a repair pathway that is specific to BRCA2 deficiency in human cancers. Integrating these features in a classifier improved discrimination between BRCA1- and BRCA2-deficient genomes. In conclusion, our data reveal classes of rearrangements that are specific to BRCA1 or BRCA2 deficiency as a source of cytogenetic aberrations in HR-deficient cells.

The screening, identification, design and clinical application of tumor-specific neoantigens for TCR-T cells

Recent advances in neoantigen research have accelerated the development of tumor immunotherapies, including adoptive cell therapies (ACTs), cancer vaccines and antibody-based therapies, particularly for solid tumors. With the development of next-generation sequencing and bioinformatics technology, the rapid identification and prediction of tumor-specific antigens (TSAs) has become possible. Compared with tumor-associated antigens (TAAs), highly immunogenic TSAs provide new targets for personalized tumor immunotherapy and can be used as prospective indicators for predicting tumor patient survival, prognosis, and immune checkpoint blockade response. Here, the identification and characterization of neoantigens and the clinical application of neoantigen-based TCR-T immunotherapy strategies are summarized, and the current status, inherent challenges, and clinical translational potential of these strategies are discussed.

A Clonal Evolution Simulator for Planning Somatic Evolution Studies

Somatic evolution plays a key role in development, cell differentiation, and normal aging, but also in diseases such as cancer. Understanding mechanisms of somatic mutability and how they can vary between cell lineages will likely play a crucial role in biological discovery and medical applications. This need has led to a proliferation of new technologies for profiling single-cell variation, each with distinctive capabilities and limitations that can be leveraged alone or in combination with other technologies. The enormous space of options for assaying somatic variation, however, presents unsolved informatics problems with regard to selecting optimal combinations of technologies for designing appropriate studies for any particular scientific questions. Versatile simulation tools are needed to explore and optimize potential study designs if researchers are to deploy multiomic technologies most effectively. In this study, we present a simulator allowing for the generation of synthetic data from a wide range of clonal lineages, variant classes, and sequencing technology choices, intended to provide a platform for effective study design in somatic lineage analysis. Users can input various properties of the somatic evolutionary system, mutation classes, and biotechnology options, and then generate samples of synthetic sequence reads and their corresponding ground truth parameters for a given study design. We demonstrate the utility of the simulator for testing and optimizing study designs for various experimental queries.

Characteristics of DNA macro-alterations in breast cancer with liver metastasis before treatment

BACKGROUND

Whole-genome doubling (WGD) has been observed in 30% of cancers, followed by a highly complex rearranged karyotype unfavourable to breast cancer's outcome. However, the macro-alterations that characterise liver metastasis in breast cancer(BC) are poorly understood. Here, we conducted a whole-genome sequencing analysis of liver metastases to explore the status and the time frame model of these macro-alterations in pre-treatment patients with metastatic breast cancer.

RESULTS

Whole-genome sequencing was conducted in 11 paired primary tumours, lymph node metastasis, and liver metastasis fresh samples from four patients with late-stage breast cancer. We also chose five postoperative frozen specimens from patients with early-stage breast cancer before any treatment as control. Surprisingly, all four liver metastasis samples were classified as WGD + . However, the previous study reported that WGD happened in 30% of cancers and 2/5 in our early-stage samples. WGD was not observed in the two separate primary tumours and one lymph node metastasis of one patient with metastatic BC, but her liver metastasis showed an early burst of bi-allelic copy number gain. The phylogenetic tree proves her 4 tumour samples were the polyclonal origin and only one WGD + clone metastasis to the liver. Another 3 metastatic BC patients' primary tumour and lymph node metastasis experienced WGD as well as liver metastasis, and they all showed similar molecular time-frame of copy number(CN) gain across locations within the same patient. These patients' tumours were of monoclonal origin, and WGD happened in a founding clone before metastasis, explaining that all samples share the CN-gain time frame. After WGD, the genomes usually face instability to evolve other macro-alterations. For example, a greater quantity and variety of complex structural variations (SVs) were detected in WGD + samples. The breakpoints were enriched in the chr17: 39 Mb-40 Mb tile, which contained the HER2 gene, resulting in the formation of tyfonas, breakage-fusion-bridge cycles, and double minutes. These complex SVs may be involved in the evolutionary mechanisms of the dramatic increase of HER2 copy number.

CONCLUSION

Our work revealed that the WGD + clone might be a critical evolution step for liver metastasis and favoured following complex SV of breast cancer.

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.

ERα-associated translocations underlie oncogene amplifications in breast cancer

Focal copy-number amplification is an oncogenic event. Although recent studies have revealed the complex structure1-3 and the evolutionary trajectories4 of oncogene amplicons, their origin remains poorly understood. Here we show that focal amplifications in breast cancer frequently derive from a mechanism-which we term translocation-bridge amplification-involving inter-chromosomal translocations that lead to dicentric chromosome bridge formation and breakage. In 780 breast cancer genomes, we observe that focal amplifications are frequently connected to each other by inter-chromosomal translocations at their boundaries. Subsequent analysis indicates the following model: the oncogene neighbourhood is translocated in G1 creating a dicentric chromosome, the dicentric chromosome is replicated, and as dicentric sister chromosomes segregate during mitosis, a chromosome bridge is formed and then broken, with fragments often being circularized in extrachromosomal DNAs. This model explains the amplifications of key oncogenes, including ERBB2 and CCND1. Recurrent amplification boundaries and rearrangement hotspots correlate with oestrogen receptor binding in breast cancer cells. Experimentally, oestrogen treatment induces DNA double-strand breaks in the oestrogen receptor target regions that are repaired by translocations, suggesting a role of oestrogen in generating the initial translocations. A pan-cancer analysis reveals tissue-specific biases in mechanisms initiating focal amplifications, with the breakage-fusion-bridge cycle prevalent in some and the translocation-bridge amplification in others, probably owing to the different timing of DNA break repair. Our results identify a common mode of oncogene amplification and propose oestrogen as its mechanistic origin in breast cancer.

Deterministic evolution and stringent selection during preneoplasia

The earliest events during human tumour initiation, although poorly characterized, may hold clues to malignancy detection and prevention1. Here we model occult preneoplasia by biallelic inactivation of TP53, a common early event in gastric cancer, in human gastric organoids. Causal relationships between this initiating genetic lesion and resulting phenotypes were established using experimental evolution in multiple clonally derived cultures over 2 years. TP53 loss elicited progressive aneuploidy, including copy number alterations and structural variants prevalent in gastric cancers, with evident preferred orders. Longitudinal single-cell sequencing of TP53-deficient gastric organoids similarly indicates progression towards malignant transcriptional programmes. Moreover, high-throughput lineage tracing with expressed cellular barcodes demonstrates reproducible dynamics whereby initially rare subclones with shared transcriptional programmes repeatedly attain clonal dominance. This powerful platform for experimental evolution exposes stringent selection, clonal interference and a marked degree of phenotypic convergence in premalignant epithelial organoids. These data imply predictability in the earliest stages of tumorigenesis and show evolutionary constraints and barriers to malignant transformation, with implications for earlier detection and interception of aggressive, genome-instable tumours.

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.

Tracking the evolution of therapy-related myeloid neoplasms using chemotherapy signatures

Patients treated with cytotoxic therapies, including autologous stem cell transplantation, are at risk for developing therapy-related myeloid neoplasms (tMN). Preleukemic clones (ie, clonal hematopoiesis [CH]) are detectable years before the development of these aggressive malignancies, although the genomic events leading to transformation and expansion are not well defined. Here, by leveraging distinctive chemotherapy-associated mutational signatures from whole-genome sequencing data and targeted sequencing of prechemotherapy samples, we reconstructed the evolutionary life-history of 39 therapy-related myeloid malignancies. A dichotomy was revealed, in which neoplasms with evidence of chemotherapy-induced mutagenesis from platinum and melphalan were hypermutated and enriched for complex structural variants (ie, chromothripsis), whereas neoplasms with nonmutagenic chemotherapy exposures were genomically similar to de novo acute myeloid leukemia. Using chemotherapy-associated mutational signatures as temporal barcodes linked to discrete clinical exposure in each patient's life, we estimated that several complex events and genomic drivers were acquired after chemotherapy was administered. For patients with prior multiple myeloma who were treated with high-dose melphalan and autologous stem cell transplantation, we demonstrate that tMN can develop from either a reinfused CH clone that escapes melphalan exposure and is selected after reinfusion, or from TP53-mutant CH that survives direct myeloablative conditioning and acquires melphalan-induced DNA damage. Overall, we revealed a novel mode of tMN progression that is not reliant on direct mutagenesis or even exposure to chemotherapy. Conversely, for tMN that evolve under the influence of chemotherapy-induced mutagenesis, distinct chemotherapies not only select preexisting CH but also promote the acquisition of recurrent genomic drivers.

Dynamics of synthetic yeast chromosome evolution shaped by hierarchical chromatin organization

Synthetic genome evolution provides a dynamic approach for systematically and straightforwardly exploring evolutionary processes. Synthetic Chromosome Rearrangement and Modification by LoxP-mediated Evolution (SCRaMbLE) is an evolutionary system intrinsic to the synthetic yeast genome that can rapidly drive structural variations. Here, we detect over 260 000 rearrangement events after the SCRaMbLEing of a yeast strain harboring 5.5 synthetic yeast chromosomes (synII, synIII, synV, circular synVI, synIXR and synX). Remarkably, we find that the rearrangement events exhibit a specific landscape of frequency. We further reveal that the landscape is shaped by the combined effects of chromatin accessibility and spatial contact probability. The rearrangements tend to occur in 3D spatially proximal and chromatin-accessible regions. The enormous numbers of rearrangements mediated by SCRaMbLE provide a driving force to potentiate directed genome evolution, and the investigation of the rearrangement landscape offers mechanistic insights into the dynamics of genome evolution.

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.

Molecular Evolution of Classic Hodgkin Lymphoma Revealed Through Whole-Genome Sequencing of Hodgkin and Reed Sternberg Cells

The rarity of malignant Hodgkin and Reed Sternberg (HRS) cells in classic Hodgkin lymphoma (cHL) limits the ability to study the genomics of cHL. To circumvent this, our group has previously optimized fluorescence-activated cell sorting to purify HRS cells. Using this approach, we now report the whole-genome sequencing landscape of HRS cells and reconstruct the chronology and likely etiology of pathogenic events leading to cHL. We identified alterations in driver genes not previously described in cHL, APOBEC mutational activity, and the presence of complex structural variants including chromothripsis. We found that high ploidy in cHL is often acquired through multiple, independent chromosomal gains events including whole-genome duplication. Evolutionary timing analyses revealed that structural variants enriched for RAG motifs, driver mutations in B2M, BCL7A, GNA13, and PTPN1, and the onset of AID-driven mutagenesis usually preceded large chromosomal gains. This study provides a temporal reconstruction of cHL pathogenesis.

SIGNIFICANCE

Previous studies in cHL were limited to coding sequences and therefore not able to comprehensively decipher the tumor complexity. Here, leveraging cHL whole-genome characterization, we identify driver events and reconstruct the tumor evolution, finding that structural variants, driver mutations, and AID mutagenesis precede chromosomal gains. This article is highlighted in the In This Issue feature, p. 171.