Faculty Bibliography

Although several recent studies have characterized structural variants (SVs) in germline and cancer genomes independently, the genomic contexts of these SVs have not been comprehensively compared. We examined similarities and differences between 2 million germline and 115 thousand tumor SVs from a cohort of 963 patients from The Cancer Genome Atlas. We found significant differences in features related to their genomic sequences and localization that suggest differences between SV-generating processes and selective pressures. For example, our results show that features linked to transposon-mediated processes are associated with germline SVs, while somatic SVs more frequently show features characteristic of chromoanagenesis. These genomic differences enabled us to develop a classifier-the Germline and Tumor Structural Variant or "the great GaTSV" -that accurately distinguishes between germline and cancer SVs in tumor samples that lack a matched normal sample.

There is currently no method to distinguish between germline and somatic structural variants (SVs) in tumor samples that lack a matched normal sample. In this study, we analyzed several features of germline and somatic SVs from a cohort of 974 patients from The Cancer Genome Atlas (TCGA). We identified a total of 21 features that differed significantly between germline and somatic SVs. Several of the germline SV features were associated with each other, as were several of the somatic SV features. We also found that these associations differed between the germline and somatic classes, for example, we found that somatic inversions were more likely to be longer events than their germline counterparts. Using these features we trained a support vector machine (SVM) classifier on 555,849 TCGA SVs to computationally distinguish germline from somatic SVs in the absence of a matched normal. This classifier had an ROC curve AUC of 0.984 when tested on an independent test set of 277,925 TCGA SVs. In this dataset, we achieved a positive predictive value (PPV) of 0.81 for an SV called somatic by the classifier being truly somatic. We further tested the classifier on a separate set of 7,623 SVs from pediatric high-grade gliomas (pHGG). In this non-TCGA cohort, our classifier achieved a PPV of 0.828, showing robust performance across datasets.

In solid tumor oncology, circulating tumor DNA (ctDNA) is poised to transform care through accurate assessment of minimal residual disease (MRD) and therapeutic response monitoring. To overcome the sparsity of ctDNA fragments in low tumor fraction (TF) settings and increase MRD sensitivity, we previously leveraged genome-wide mutational integration through plasma whole-genome sequencing (WGS). Here we now introduce MRD-EDGE, a machine-learning-guided WGS ctDNA single-nucleotide variant (SNV) and copy-number variant (CNV) detection platform designed to increase signal enrichment. MRD-EDGE^SNV uses deep learning and a ctDNA-specific feature space to increase SNV signal-to-noise enrichment in WGS by ~300× compared to previous WGS error suppression. MRD-EDGE^CNV also reduces the degree of aneuploidy needed for ultrasensitive CNV detection through WGS from 1 Gb to 200 Mb, vastly expanding its applicability within solid tumors. We harness the improved performance to identify MRD following surgery in multiple cancer types, track changes in TF in response to neoadjuvant immunotherapy in lung cancer and demonstrate ctDNA shedding in precancerous colorectal adenomas. Finally, the radical signal-to-noise enrichment in MRD-EDGE^SNV enables plasma-only (non-tumor-informed) disease monitoring in advanced melanoma and lung cancer, yielding clinically informative TF monitoring for patients on immune-checkpoint inhibition.

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.

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-year-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-year-old boy with a thalamic tumor that underwent multiple resections, with no actionable driver detected using targeted next-generation sequencing (NGS). Whole-genome and RNAseq 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.

Homologous recombination (HR) deficiency is associated with DNA rearrangements and cytogenetic aberrations¹. Paradoxically, the types of DNA rearrangements that are specifically associated with HR-deficient cancers only minimally affect chromosomal structure². 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.

UNLABELLED: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.

Fanconi anaemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink repair resulting in chromosome breakage^1-3. The FA repair pathway protects against endogenous and exogenous carcinogenic aldehydes^4-7. Individuals with FA are hundreds to thousands fold more likely to develop head and neck (HNSCC), oesophageal and anogenital squamous cell carcinomas⁸ (SCCs). Molecular studies of SCCs from individuals with FA (FA SCCs) are limited, and it is unclear how FA SCCs relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or infection with human papillomavirus⁹ (HPV). Here, by sequencing genomes and exomes of FA SCCs, we demonstrate that the primary genomic signature of FA repair deficiency is the presence of high numbers of structural variants. Structural variants are enriched for small deletions, unbalanced translocations and fold-back inversions, and are often connected, thereby forming complex rearrangements. They arise in the context of TP53 loss, but not in the context of HPV infection, and lead to somatic copy-number alterations of HNSCC driver genes. We further show that FA pathway deficiency may lead to epithelial-to-mesenchymal transition and enhanced keratinocyte-intrinsic inflammatory signalling, which would contribute to the aggressive nature of FA SCCs. We propose that the genomic instability in sporadic HPV-negative HNSCC may arise as a result of the FA repair pathway being overwhelmed by DNA interstrand crosslink damage caused by alcohol and tobacco-derived aldehydes, making FA SCC a powerful model to study tumorigenesis resulting from DNA-crosslinking damage.

Myoepithelial carcinomas (MECs) of soft tissue are rare and aggressive tumors affecting young adults and children, but their molecular landscape has not been comprehensively explored through genome sequencing. Here, we present the whole-exome sequencing (WES), whole-genome sequencing (WGS), and RNA sequencing findings of two MECs. Patients 1 and 2 (P1, P2), both male, were diagnosed at 27 and 37 yr of age, respectively, with shoulder (P1) and inguinal (P2) soft tissue tumors. Both patients developed metastatic disease, and P2 died of disease. P1 tumor showed a rhabdoid cytomorphology and a complete loss of INI1 (SMARCB1) expression, associated with a homozygous SMARCB1 deletion. The tumor from P2 showed a clear cell/small cell morphology, retained INI1 expression and strong S100 positivity. By WES and WGS, tumors from both patients displayed low tumor mutation burdens, and no targetable alterations in cancer genes were detected. P2's tumor harbored an EWSR1::KLF15 rearrangement, whereas the tumor from P1 showed a novel ASCC2::GGNBP2 fusion. WGS evidenced a complex genomic event involving mainly Chromosomes 17 and 22 in the tumor from P1, which was consistent with chromoplexy. These findings are consistent with previous reports of EWSR1 rearrangements (50% of cases) in MECs and provide a genetic basis for the loss of SMARCB1 protein expression observed through immunohistochemistry in 10% of 40% of MEC cases. The lack of additional driver mutations in these tumors supports the hypothesis that these alterations are the key molecular events in MEC evolution. Furthermore, the presence of complex structural variant patterns, invisible to WES, highlights the novel biological insights that can be gained through the application of WGS to rare cancers.