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.

Prediction of the 3D cancer genome from whole-genome sequencing using InfoHiC

The 3D genome prediction in cancer is crucial for uncovering the impact of structural variations (SVs) on tumorigenesis, especially when they are present in noncoding regions. We present InfoHiC, a systemic framework for predicting the 3D cancer genome directly from whole-genome sequencing (WGS). InfoHiC utilizes contig-specific copy number encoding on the SV contig assembly, and performs a contig-to-total Hi-C conversion for the cancer Hi-C prediction from multiple SV contigs. We showed that InfoHiC can predict 3D genome folding from all types of SVs using breast cancer cell line data. We applied it to WGS data of patients with breast cancer and pediatric patients with medulloblastoma, and identified neo topologically associating domains. For breast cancer, we discovered super-enhancer hijacking events associated with oncogenic overexpression and poor survival outcomes. For medulloblastoma, we found SVs in noncoding regions that caused super-enhancer hijacking events of medulloblastoma driver genes (GFI1, GFI1B, and PRDM6). In addition, we provide trained models for cancer Hi-C prediction from WGS at https://github.com/dmcb-gist/InfoHiC , uncovering the impacts of SVs in cancer patients and revealing novel therapeutic targets.

Chromatin interaction maps identify oncogenic targets of enhancer duplications in cancer

As a major type of structural variants, tandem duplication plays a critical role in tumorigenesis by increasing oncogene dosage. Recent work has revealed that noncoding enhancers are also affected by duplications leading to the activation of oncogenes that are inside or outside of the duplicated regions. However, the prevalence of enhancer duplication and the identity of their target genes remains largely unknown in the cancer genome. Here, by analyzing whole-genome sequencing data in a non-gene-centric manner, we identify 881 duplication hotspots in 13 major cancer types, most of which do not contain protein-coding genes. We show that the hotspots are enriched with distal enhancer elements and are highly lineage-specific. We develop a HiChIP-based methodology that navigates enhancer-promoter contact maps to prioritize the target genes for the duplication hotspots harboring enhancer elements. The methodology identifies many novel enhancer duplication events activating oncogenes such as ESR1, FOXA1, GATA3, GATA6, TP63, and VEGFA, as well as potentially novel oncogenes such as GRHL2, IRF2BP2, and CREB3L1 In particular, we identify a duplication hotspot on Chromosome 10p15 harboring a cluster of enhancers, which skips over two genes, through a long-range chromatin interaction, to activate an oncogenic isoform of the NET1 gene to promote migration of gastric cancer cells. Focusing on tandem duplications, our study substantially extends the catalog of noncoding driver alterations in multiple cancer types, revealing attractive targets for functional characterization and therapeutic intervention.

Comprehensive evaluation and guidance of structural variation detection tools in chicken whole genome sequence data

BACKGROUND

Structural variations (SVs) are widespread across genome and have a great impact on evolution, disease, and phenotypic diversity. Despite the development of numerous bioinformatic tools, commonly referred to as SV callers, tailored for detecting SVs using whole genome sequence (WGS) data and employing diverse algorithms, their performance necessitates rigorous evaluation with real data and validated SVs. Moreover, a considerable proportion of these tools have been primarily designed and optimized using human genome data. Consequently, their applicability and performance in Avian species, characterized by smaller genomes and distinct genomic architectures, remain inadequately assessed.

RESULTS

We performed a comprehensive assessment of the performance of ten widely used SV callers using population-level real genomic data with the validated five common types of SVs. The performance of SV callers varies with the types and sizes of SVs. As compared with other tools, GRIDSS, Lumpy, Wham, and Manta present better detection accuracy. Pindel can detect more small SVs than others. CNVnator and CNVkit can detect more medium and large copy number variations. Given the poor consistency among different SV callers, the combination calling strategy is not recommended. All tools show poor ability in the detection of insertions (especially with size > 150 bp). At least 50× read depth is required to detect more than 80% of the SVs for most tools.

CONCLUSIONS

This study highlights the importance and necessity of using real sequencing data, rather than simulated data only, with validated SVs for SV caller evaluation. Some practical guidance and suggestions are provided for SV detection in future researches.

Modeling lung adenocarcinoma metastases using patient-derived organoids

Approximately 50% of patients with surgically resected early-stage lung cancer develop distant metastasis. At present, there is no in vivo metastasis model to investigate the biology of human lung cancer metastases. Using well-characterized lung adenocarcinoma (LUAD) patient-derived organoids (PDOs), we establish an in vivo metastasis model that preserves the biologic features of human metastases. Results of whole-genome and RNA sequencing establish that our in vivo PDO metastasis model can be used to study clonality and tumor evolution and to identify biomarkers related to organotropism. Investigation of the response of KRASG12C PDOs to sotorasib demonstrates that the model can examine the efficacy of treatments to suppress metastasis and identify mechanisms of drug resistance. Finally, our PDO model cocultured with autologous peripheral blood mononuclear cells can potentially be used to determine the optimal immune-priming strategy for individual patients with LUAD.

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 ∼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, in which the inferred architecture of ecDNA is encoded as a cycle in the graph. Traversals of breakpoint graphs 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 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 Complete Reconstruction of Amplifications with Long reads (CoRAL) 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 10 data sets from previously characterized cell lines compared with previous short- and long-read-based tools. As long-read usage becomes widespread, we anticipate that CoRAL will be a valuable tool for profiling the landscape and evolution of focal amplifications in tumors.

The interplay of mutagenesis and ecDNA shapes urothelial cancer evolution

Advanced urothelial cancer is a frequently lethal disease characterized by marked genetic heterogeneity1. In this study, we investigated the evolution of genomic signatures caused by endogenous and external mutagenic processes and their interplay with complex structural variants (SVs). We superimposed mutational signatures and phylogenetic analyses of matched serial tumours from patients with urothelial cancer to define the evolutionary dynamics of these processes. We show that APOBEC3-induced mutations are clonal and early, whereas chemotherapy induces mutational bursts of hundreds of late subclonal mutations. Using a genome graph computational tool2, we observed frequent high copy-number circular amplicons characteristic of extrachromosomal DNA (ecDNA)-forming SVs. We characterized the distinct temporal patterns of APOBEC3-induced and chemotherapy-induced mutations within ecDNA-forming SVs, gaining new insights into the timing of these mutagenic processes relative to ecDNA biogenesis. We discovered that most CCND1 amplifications in urothelial cancer arise within circular ecDNA-forming SVs. ecDNA-forming SVs persisted and increased in complexity, incorporating additional DNA segments and contributing to the evolution of treatment resistance. Oxford Nanopore Technologies long-read whole-genome sequencing followed by de novo assembly mapped out CCND1 ecDNA structure. Experimental modelling of CCND1 ecDNA confirmed its role as a driver of treatment resistance. Our findings define fundamental mechanisms that drive urothelial cancer evolution and have important therapeutic implications.

53BP1 deficiency leads to hyperrecombination using break-induced replication (BIR)

Break-induced replication (BIR) is mutagenic, and thus its use requires tight regulation, yet the underlying mechanisms remain elusive. Here we uncover an important role of 53BP1 in suppressing BIR after end resection at double strand breaks (DSBs), distinct from its end protection activity, providing insight into the mechanisms governing BIR regulation and DSB repair pathway selection. We demonstrate that loss of 53BP1 induces BIR-like hyperrecombination, in a manner dependent on Polα-primase-mediated end fill-in DNA synthesis on single-stranded DNA (ssDNA) overhangs at DSBs, leading to PCNA ubiquitination and PIF1 recruitment to activate BIR. On broken replication forks, where BIR is required for repairing single-ended DSBs (seDSBs), SMARCAD1 displaces 53BP1 to facilitate the localization of ubiquitinated PCNA and PIF1 to DSBs for BIR activation. Hyper BIR associated with 53BP1 deficiency manifests template switching and large deletions, underscoring another aspect of 53BP1 in suppressing genome instability. The synthetic lethal interaction between the 53BP1 and BIR pathways provides opportunities for targeted cancer treatment.

FindCSV: a long-read based method for detecting complex structural variations

BACKGROUND

Structural variations play a significant role in genetic diseases and evolutionary mechanisms. Extensive research has been conducted over the past decade to detect simple structural variations, leading to the development of well-established detection methods. However, recent studies have highlighted the potentially greater impact of complex structural variations on individuals compared to simple structural variations. Despite this, the field still lacks precise detection methods specifically designed for complex structural variations. Therefore, the development of a highly efficient and accurate detection method is of utmost importance.

RESULT

In response to this need, we propose a novel method called FindCSV, which leverages deep learning techniques and consensus sequences to enhance the detection of SVs using long-read sequencing data. Compared to current methods, FindCSV performs better in detecting complex and simple structural variations.

CONCLUSIONS

FindCSV is a new method to detect complex and simple structural variations with reasonable accuracy in real and simulated data. The source code for the program is available at https://github.com/nwpuzhengyan/FindCSV .

53BP1 deficiency leads to hyperrecombination using break-induced replication (BIR)

Break-induced replication (BIR) is mutagenic, and thus its use requires tight regulation, yet the underlying mechanisms remain elusive. Here we uncover an important role of 53BP1 in suppressing BIR after end resection at double strand breaks (DSBs), distinct from its end protection activity, providing insight into the mechanisms governing BIR regulation and DSB repair pathway selection. We demonstrate that loss of 53BP1 induces BIR-like hyperrecombination, in a manner dependent on Polα-primase-mediated end fill-in DNA synthesis on single-stranded DNA (ssDNA) overhangs at DSBs, leading to PCNA ubiquitination and PIF1 recruitment to activate BIR. On broken replication forks, where BIR is required for repairing single-ended DSBs (seDSBs), SMARCAD1 displaces 53BP1 to facilitate the localization of ubiquitinated PCNA and PIF1 to DSBs for BIR activation. Hyper BIR associated with 53BP1 deficiency manifests template switching and large deletions, underscoring another aspect of 53BP1 in suppressing genome instability. The synthetic lethal interaction between the 53BP1 and BIR pathways provides opportunities for targeted cancer treatment.

Tracking clonal evolution of drug resistance in ovarian cancer patients by exploiting structural variants in cfDNA

Drug resistance is the major cause of therapeutic failure in high-grade serous ovarian cancer (HGSOC). Yet, the mechanisms by which tumors evolve to drug resistant states remains largely unknown. To address this, we aimed to exploit clone-specific genomic structural variations by combining scaled single-cell whole genome sequencing with longitudinally collected cell-free DNA (cfDNA), enabling clonal tracking before, during and after treatment. We developed a cfDNA hybrid capture, deep sequencing approach based on leveraging clone-specific structural variants as endogenous barcodes, with orders of magnitude lower error rates than single nucleotide variants in ctDNA (circulating tumor DNA) detection, demonstrated on 19 patients at baseline. We then applied this to monitor and model clonal evolution over several years in ten HGSOC patients treated with systemic therapy from diagnosis through recurrence. We found drug resistance to be polyclonal in most cases, but frequently dominated by a single high-fitness and expanding clone, reducing clonal diversity in the relapsed disease state in most patients. Drug-resistant clones frequently displayed notable genomic features, including high-level amplifications of oncogenes such as CCNE1, RAB25, NOTCH3, and ERBB2. Using a population genetics Wright-Fisher model, we found evolutionary trajectories of these features were consistent with drug-induced positive selection. In select cases, these alterations impacted selection of secondary lines of therapy with positive patient outcomes. For cases with matched single-cell RNA sequencing data, pre-existing and genomically encoded phenotypic states such as upregulation of EMT and VEGF were linked to drug resistance. Together, our findings indicate that drug resistant states in HGSOC pre-exist at diagnosis and lead to dramatic clonal expansions that alter clonal composition at the time of relapse. We suggest that combining tumor single cell sequencing with cfDNA enables clonal tracking in patients and harbors potential for evolution-informed adaptive treatment decisions.

Extrachromosomal circular DNA orchestrates genome heterogeneity in urothelial bladder carcinoma

Rationale: Extrachromosomal circular DNA is a hallmark of cancer, but its role in shaping the genome heterogeneity of urothelial bladder carcinoma (UBC) remains poorly understood. Here, we comprehensively analyzed the features of extrachromosomal circular DNA in 80 UBC patients. Methods: We performed whole-genome/exome sequencing (WGS/WES), Circle-Seq, single-molecule real-time (SMRT) long-read sequencing of circular DNA, and RNA sequencing (RNA-Seq) on 80 pairs of tumor and AT samples. We used our newly developed circular DNA analysis software, Circle-Map++ to detect small extrachromosomal circular DNA from Circle-Seq data. Results: We observed a high load and significant heterogeneity of extrachromosomal circular DNAs in UBC, including numerous single-locus and complex chimeric circular DNAs originating from different chromosomes. This includes highly chimeric circular DNAs carrying seven oncogenes and circles from nine chromosomes. We also found that large tumor-specific extrachromosomal circular DNAs could influence genome-wide gene expression, and are detectable in time-matched urinary sediments. Additionally, we found that the extrachromosomal circular DNA correlates with hypermutation, copy number variation, oncogene amplification, and clinical outcome. Conclusions: Overall, our study provides a comprehensive extrachromosomal circular DNA map of UBC, along with valuable data resources and bioinformatics tools for future cancer and extrachromosomal circular DNA research.

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.

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.

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.

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.

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.

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.

The Pathology according to p53 Pathway

BACKGROUND

Observations play a pivotal role in the progress of science, including in pathology. The cause of a disease such as cancer is analyzed by breaking it down into smaller organs, tissues, cells, and molecules. The current standard cancer diagnostic procedure, microscopic observation, relies on preserved morphological characteristics. In contrast, molecular analyses explore oncogenic pathway activation that leads to genetic mutations and aberrant protein expression. Such molecular analyses could potentially identify therapeutic targets and has gained considerable attention in clinical oncology.

SUMMARY

This review summarizes the cardinal biomarkers of the p53 pathway, p53, p16, and mouse double minute 2 (MDM2), in the context of traditional surgical pathology and emerging genomic oncology. The p53 pathway, which is dysregulated in more than a half of all cancers, can be applied in several diagnostic settings. A four-classification model of immunophenotype for p53 pathway gene status, tumor types with a high frequency of abnormalities for each p53 pathway gene, and a minimal p53 pathway immunohistochemical panel is also described.

KEY MESSAGES

Immunohistochemistry of oncogenic signals should be interpreted according to molecular findings based on genomic oncology, in addition to the microscopic findings of diagnostic pathology.

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.