Novel patterns of complex structural variation revealed across thousands of cancer genome graphs

Abstract

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g. deletion, translocation) or complex (e.g. chromothripsis, chromoplexy) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,833 tumor whole genome sequences (WGS), we introduce three complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early replicating regions and superenhancers, and are enriched in breast and ovarian cancers. Rigma comprise "chasms" of low-JCN deletions at late-replicating fragile sites in esophageal and other gastrointestinal (GI) adenocarcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold back inversions that are enriched in acral but not cutaneous melanoma and associated with a previously uncharacterized mutational process of non-APOBEC kataegis. Clustering of tumors according to genome graph-derived features identifies subgroups associated with DNA repair defects and poor prognosis. Cancer genomics | Structural variation | DNA rearrangements | Mutational Processes | Genome graphs | Copy number alterationsCorrespondence: mski@mskilab.org K. Hadi, X.Yao, J. Behr et al. | bioRχiv | November 12, 2019 | 1-8

Nanopore sequencing of DNA concatemers reveals higher-order features of chromatin structure

“…As predicted from our Pore-C data, these three probes were proximal and amplified in the HCC1954 cell line, and distinct and diploid in ANA51 (Fig 3B, top). These results confirm the presence of a com- //github.com/mskilab/JaBbA) 45 . Analysis of binned read-depth (bottom track) and junction calls yielded the genome graph (second track from bottom).…”

“…High-order chromatin contacts resolve the allelic structure of a complex cancer rearrangement. Structural DNA variants are common in cancer genomes, yielding complex loci harboring many copies of genomic intervals and rearrangement junctions in cis or trans phase 45,46 .…”

“…Junction and read depth data were input into a junction-balanced genome graph inference algorithm JaBbA (https://github. com/mskilab/JaBbA) 45 . JaBbA integrates read depth and structural variation information for a given genome and produces a junction-balanced genome graph by solving a constrained optimization problem 45 .…”

“…com/mskilab/JaBbA) 45 . JaBbA integrates read depth and structural variation information for a given genome and produces a junction-balanced genome graph by solving a constrained optimization problem 45 . The junction-balanced graph allows querying of non-contiguous regions of genome that become connected through rearrangements through the R package gGnome (https://github.com/mskilab/ gGnome).…”

Enhancer hijacking determines intra- and extrachromosomal circularMYCNamplicon architecture in neuroblastoma

“…As we were only aiming at the rearrangements common to all amplicons, we only considered breakpoints with more than 50 variant-support reads ('allele depth'). gGnome 40 was used to represent these data as a genome graph with nodes being breakpoint-free genomic intervals and edges being rearrangements ('alternate edge') or connections in the reference genomes ('reference edge'). We considered only nodes with high copy number, i.e.…”