Abstract:Introduction: Structural variants (SVs) play a role in proliferation of plasma cells that characterizes multiple myeloma (MM). Several hallmark SVs have been identified in MM, including translocations involving the immunoglobulin heavy chain locus and several partner genes, as well as deletions and duplications affecting key genes such as TP53, RB1, and MYC. Complex SVs such as templated insertions, chromothripsis, and chromoplexy are associated with disease progression but have not yet systematically been studied. Methods: Paired short-read (Illumina, average depth 75x), long-read HiFi (PacBio, 16x), optical genome mapping (OGM; Bionano Genomics, 341x) and Micro-C (Cantata Bio, 45x) data were sequenced from 8 relapsed MM patient-derived xenograft samples and 2 cell lines, covering t(11;14), t(4;14), and hyperdiploid subgroups. Reads were aligned to the hg38 reference genome while complex SVs (chromothripsis, chromoplexy and templated insertion) were systematically identified and accurately reconstructed by combining short-read (Manta), HiFi (pbsv), OGM, and Micro-C data (EagleC). Chromosome organization including loops, topologically associated domains (TADs) and A/B compartments were called from Micro-C data (Mustache and Juicer). Haplotype phasing was conducted using small variants detected from HiFi data. De novo assembly of the genome was conducted using HiFi and Micro-C data and compared to the reference genome. Reconstructed derivative chromosome interactions due to complex SVs were further identified by aligning Micro-C reads to the newly assembled genomes. Previously identified super-enhancers from H3K27Ac ChIP-seq from MM cell lines were mapped onto these SVs and expression of target genes assessed by RNA-seq. Results: 1,340 SVs were found in 10 samples using OGM of which 828 were identified in short-read WGS data. 16 complex SVs including chromothripsis, chromoplexy and templated insertions were detected in 10 samples, and we describe two (Fig.1) in detail here. A complex templated insertion consisting of a t(2;6;8) was found in a t(11;14) PDX sample with copy number gains at breakpoints which could not be resolved from short-read data alone. OGM data revealed a single contig consisting of DNA from chromosomes 2, 6, and 8. Micro-C data revealed interactions between these chromosomes, specifically the gained regions on chromosomes 2 and 8 showed more interactions with chromosome 6 than with other chromosomes, indicating insertions into chromosome 6. As a result, both MYC at chromosome 8 and GALM at chromosome 2 were inserted into the TXNDC5 locus, containing a known MM super-enhancer, leading to the over-expression of both MYC and GALM (log 2FC=1.7 and 1.0) in this sample. A de novo genome assembly indicated neo-TAD formation encompassing the super-enhancer and both genes. In coMMpass dataset, high expression of GALM and MYC are significantly associated with inferior overall (p=0.002 and 0.05, logrank test) and progression-free survival (p=0.01 and 0.09). Another complex templated insertion consisting of a t(3;4;16) was identified in a hyperdiploid PDX sample, again with copy number gains at breakpoints. Although the short-read data indicated an assembly similar to the previous templated insertion SV, the OGM and Micro-C data indicated a 3-way chromosome swap where the short arm of chromosome 3 was translocated to the short arm of chromosome 4, the short arm of chromosome 4 was translocated to the long arm of chromosome 16, and the long arm of chromosome 16 was translocated to the short arm of chromosome 3. Separate contigs for the t(4;16), t(3;16) and t(3;4) were confirmed by OGM data and only 2-way chromosomal interactions were seen in the Micro-C data. This complex event led to numerous oncogenes juxtaposed to super-enhancers, which caused over-expression of target genes such as ACD and RANBP10 (log 2FC=1.0 and 0.4), both of which were on chromosome 16 and are significantly associated with inferior overall (p=0.02 and p=0.004) and progression-free survival (p=0.003 and p=0.005). Conclusion: By utilizing multiomic techniques we were able to reconstruct complex SV events, showing different mechanisms of generating templated insertion events and identifying their biological consequences. We demonstrated how these complex events could potentially facilitate tumor proliferation via engaging the super-enhancers and regulating the expression of key oncogenes.

Resolving the chromatin impact of mosaic variants with targeted Fiber-seq

A personalized multi-platform assessment of somatic mosaicism in the human frontal cortex

Cell-type-informed genotyping of mosaic focal epilepsies reveals cell-autonomous and non-cell-autonomous disease-associated transcriptional programs

Comprehensive identification of somatic nucleotide variants in human brain tissue

Comparative analysis of CRISPR/Cas9-targeted nanopore sequencing approaches in repeat expansion disorders

Deaminase-assisted single-molecule and single-cell chromatin fiber sequencing

Accurate detection of mosaic variants in sequencing data without matched controls

Efficient Cnv Breakpoint Analysis Reveals Unexpected Structural Complexity and Correlation of Dosage-Sensitive Genes with Clinical Severity in Genomic Disorders

Massively parallel single-cell mitochondrial DNA genotyping and chromatin profiling

Amplification-free, CRISPR-Cas9 Targeted Enrichment and SMRT Sequencing of Repeat-Expansion Disease Causative Genomic Regions

Abstract 7026: Multi-omic genomic mapping with long read sequencing

Abstract LB081: Concurrent characterization of somatic nucleotide mutations, 5mC, and 5hmC methylation in clinical research samples using a single nanopore sequencing assay

Unraveling Diverse Mechanisms of Complex Structural Variant Interactions through Multiomic Data in Multiple Myeloma

Genotype-to-phenotype mapping of somatic clonal mosaicism via single-cell co-capture of DNA mutations and mRNA transcripts

Characterization of FMR1 Repeat Expansion and Intragenic Variants by Indirect Sequence Capture

iFlpMosaics enable the multispectral barcoding and high-throughput comparative analysis of mutant and wild-type cells

Sequence-to-expression approach to identify etiological non-coding DNA variations in P53 and cMYC-driven diseases

High Quality Phasing Using Linked-Read Whole Genome Sequencing of Patient Cohorts Informs Genetic Understanding of Complex Traits

Enhanced Detection and Genotyping of Disease-Associated Tandem Repeats Using HMMSTR and Targeted Long-Read Sequencing

An update on the neurological short tandem repeat expansion disorders and the emergence of long-read sequencing diagnostics

Abstract 5400: Accurate identification of single nucleotide variants in whole genome amplified single cells