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.

3D Genome of Multiple Myeloma Reveals Spatial Genome Disorganization Associated with Copy Number Variations.

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

Single Cell Whole Genome Measurement of Genomic Heterogeneity and Clonal Evolution in Multiple Myeloma

Insights into the Basis of Chromosomal Imbalances During the Clonal Evolution of Multiple Myeloma Using Snp Array Analysis.

Abstract 5676: Spatial 3D genome organization reveals intratumor heterogeneity in primary glioblastoma samples

Integration of Global SNP-based Mapping and Expression Arrays Reveals Key Regions, Mechanisms, and Genes Important in the Pathogenesis of Multiple Myeloma.

[Detection of the Cytogenetic Aberrations in Multiple Myeloma by Using Microrray Comparative Genomic Hybridization].

Single-molecule analysis reveals widespread structural variation in multiple myeloma

Differential Chromatin Organization between t(4;14) and Non-t(4;14) Multiple Myeloma Driven By the Histone Methyltransferase NSD2

Multi-dimensional scaling techniques unveiled gain1q&loss13q co-occurrence in Multiple Myeloma patients with specific genomic, transcriptional and adverse clinical features

Single-cell sequencing analysis of multiple myeloma heterogeneity and identification of new theranostic targets

A single-cell atlas characterizes dysregulation of the bone marrow immune microenvironment associated with outcomes in multiple myeloma

Application of Hi-C technology in three-dimensional genomics research and disease pathogenesis analysis.

Integrative Analysis of Genome, 3D Genome, and Transcriptome Alterations of Clinical Lung Cancer Samples.

Structural Variations of the 3D Genome Architecture in Cervical Cancer Development

Chromosomal Deletions and Amplifications in Multiple Myeloma Detected by 500k Single Nucleotide Polymorphism Array Analysis.

Single-Cell Transcriptomes Combining with Consecutive Genomics Reveal Clonal Evolution and Gene Regulatory Networks in Relapsed and Refractory Multiple Myeloma

3D Genomic Organization in Cancers

A Retrospective Analysis of Cytogenetic and Clinical Characteristics in Patients with Multiple Myeloma

Clonal Heterogeneity and Evolutionary Phylogeny of Critical Cytogenetic Aberrations in Multiple Myeloma

Three-dimensional chromatin landscapes in MLLr AML