Abstract:INTRODUCTION: Multiple myeloma (MM) evolution is complex and heterogeneous. Hyperdiploidy (HRD) and translocations affecting the immunoglobulin heavy chain (IGH) locus are historically considered initiating genomic events of MM. The potential impact of genomic events acquired before known MM initiating events has never been addressed. METHODS: To investigate whether genomic deletions are acquired before and after known MM initiating events we interrogated whole-genome sequencing (WGS) data from patients with newly diagnosed MM (NDMM, n=319) and relapsed MM (RRMM, n=59). Our newly developed analytical and chronological workflow integrating single nucleotide variants (SNV), structural variants (SV), and copy number variants (CNV) comprises two main steps. First, we estimated the molecular time (i.e. corrected ratio between duplicated and non-duplicated clonal SNV) of large clonal chromosomal duplications and identified the earliest set of chromosomal gains in each patient. Second, within each early large chromosomal gain we identified clonal SV mediating CNV loss. A deletion on a gain can generate three possible scenarios: 1) one of the duplicated alleles is lost after the gain (i.e. post-gain) causing a CNV jump from 3:1 to 2:1 (total alleles : minor alleles); 2) there is a deletion before the duplication (i.e., pre-gain) causing a CNV jump from 3:1 to 1:0. 3) the deletion occurs on the minor, non-duplicated allele, causing a CNV jump from 3:1 to 2:0. Timing the deletion in relation with the chromosomal duplication is impossible in this scenario. RESULTS: Molecular time data were successfully generated for 249/319 (78%) NDMM and 51/59 (86%) RRMM patients. Restricting our analysis to NDMM with HRD without canonical IGH translocations, 16/170 (9%) of patients acquired deletions before the earliest multi-chromosomal gains, suggesting MM precursors can acquire deletions before HRD. Investigating the entire series, post-gain deletions were observed in 126/417 (30%) samples considering both early and late time windows. Leveraging a background model, we demonstrated that pre-gain events involved more tumor-suppressor genes (TSG) than expected by chance, including TCF3, ATM and TRAF3. In contrast, oncogenes were involved less than expected. In post-gain deletions, both oncogenes and TSG were involved more than expected by chance . To validate and assess the impact of deletions on driver gene expression we investigated WGS and RNAsequencing data from the MMRF CoMMpass study (n=752 NDMM). Because of the low coverage not allowing for molecular time estimation, we limited our analysis to HRD without IGH translocations. Pre-gain and post-gain deletions were observed in 47/431 (11%) and 225/431 (52%) of patients, respectively. We defined loss and gain of function events based on whether the expression level of a gene affected by the deletion was in the first or fourth quartile, respectively. Pre-gain deletions were mostly associated with downregulation of TSG expression [n=31 events in 13/431 (3%) patients]. In contrast,post-gain deletions had a more heterogenous impact with loss- and gain-of-function events. Gain-of-function events were driven by two main mechanisms: 1) the deletion joined an oncogene with a distal regulatory region inducing its overexpression [n=44 events in 6/431 (1.3%) patients]; 2) the deletion caused a new and expressed fusion [n=231 events in 109/431 (25%) patients], resulting in either loss- or gain-of-function. Surprisingly, post-gain deletions had also a major impact on TSG expression. In 60 patients (14%), we observed 146 post-gain deletions where the affected tumor suppressor gene (TSG) expression was downregulated to the level of cases with monoallelic and biallelic deletions, even though two alleles were retained. Finally, to validate these findings, we investigated 16 RRMM patients with available WGS, scATAC-seq and scRNA-seq and observed additional evidence of TSG downregulation after both pre- and post-gain deletions. CONCLUSION: Leveraging a large cohort of NDMM we show that somatic deletions can be acquired before HRD trisomies that are assumed to be initiating events. Furthermore, post-gain deletions emerged as a new mechanism inducing TSG down-regulation, despite an apparently diploid gene status.

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

High-Resolution Genomic Profiles Identify Novel Genes And/Or Chromosomal Regions with Prognostic and Oncogenic Significance in Myeloma Patients

Prognostic Significance Of Copy-Number Alterations In Multiple Myeloma

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

Use Of High-Density Snp-Array Analysis To Identify Novel Chromosomal Abnormalities That Predict Survival In Multiple Myeloma

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].

Identification of Collaborating Oncogeneic Events Leading to Disease Progression in Myeloma Cases with A T(4;14) and T(11;14) Using Snp and Gene Expression Arrays.

Retrospective Analysis of Cytogenetic and Clinical Characteristics on 65 Patients with Multiple Myeloma

Prognostic Power of Abnormal Cytogenetics for Multiple Myeloma: a Multicenter Study in China.

Fluorescence in Situ Hybridization Analysis of Chromosome Aberrations in 60 Chinese Patients with Multiple Myeloma

Status of Chromosome 13 in Multiple Myeloma: Integrated Approach Using SNP Mapping Array and Gene Expression Array.

Snp Array Analysis In Hematologic Malignancies: Avoiding False Discoveries

Revealing Novel Mechanisms Underlying Inactivation of Tumor Suppressor Genes on Duplicated Chromosomes in Multiple Myeloma

Single Cell Multi-Omic Correlation of Single Nucleotide Variants, Copy Number Variation and Surface Epitopes for Clonal Profiling of Myeloma

Chromosomal defects in multiple myeloma

Bacterial Artificial Chromosome Array‐based Comparative Genomic Hybridization Using Paired Formalin‐fixed, Paraffin‐embedded and Fresh Frozen Tissue Specimens in Multiple Myeloma

Peripheral Blood Single-Cell Genotyping and Phenotyping in Multiple Myeloma Reveals Shared Mutations across Multiple Hematopoietic Cell Lineages

Allele-Specific Amplification in Cancer Revealed by SNP Array Analysis.

Acquisition of New in Vivo Genomic Change over Time in Multiple Myeloma.

Detection of Recurrent Cytogenetic Aberrations in Multiple Myeloma: A Comparison Between MLPA and Ifish