Genomic Complexity Correlates with the Degree of Marrow Independence of Malignant Plasma Cells in the Context of Extramedullary Disease
Authors: Maclachlan Kylee H , Garces Juan-Jose , Shekarkhand Tala , Rajeeve Sridevi , Hashmi Hamza , Hassoun Hani , Hultcrantz Malin , Korde Neha , Tan Carlyn Rose , Mailankody Sham , Simhal Anish , Chung David J. , Shah Gunjan L. , Scordo Michael , Giralt Sergio , Poos Alexandra Maria , John Lukas , Zhang Yanming , Cimera Robert , Arcila Maria E. , Arango Ossa Juan E , Bhanot Umesh , Jain Mala , Merghoub Tala , Landau Heather J. , Lesokhin Alexander M. , Usmani Saad Z. , Weinhold Niels , Maura Francesco , Shah Urvi A. 
Published: December 2024

Background: Growth of multiple myeloma (MM) beyond the bone marrow (BM) represents clinically aggressive disease and is associated with inferior therapy response. We hypothesize that genomic changes underlie the increasing BM independence seen from focal intra-medullary lesions (IMD) to para-medullary (PMD) and extra-medullary disease (EMD). To address this hypothesis, we have analyzed the largest sequencing cohort of IMD/PMD/EMD to date, comparing with paired BM samples where available.
Methods: Genomic assessment of all available IMD/PMD/EMD from Memorial Sloan Kettering Cancer Center (MSK) and the University of Heidelberg (DKFZ) were included. IMD was defined by MRI appearance, PMD by local extension through bone and EMD as non-contiguous with bone. Comprehensive genomic characterization combined SNP-array, extended targeted sequencing (MSK-IMPACT-Heme), whole exome (WES) and whole genome sequencing (WGS). New data was combined with our published pre-therapy BM WGS (Maura Nat Can 2023), PMD/EMD WGS from autopsy (Landau Nat Comm 2020) and spatially distinct WES (Rasche Nat Comm 2022).
In total, 84 IMD, 105 PMD and 55 EMD samples were assessed; 128 pre-therapy and 126 post-therapy. 134 patients had paired BM samples (113 coincident and 26 distant in time), with 12 patients having &gt; 1 PMD/EMD biopsy. Including our previously published data, the cohort comprised 444 samples from 271 patients; 125 WGS, 66 WES, 145 MSK-IMPACT-Heme and 25 SNP-array. Analyses included GISTIC for significant copy number aberrations (CNA), dndscv for positively selected mutational drivers and mmsig for mutational signatures.
Results: Compared to paired BM and the independent control set of BM samples, PMD and EMD had considerably more CNA, while IMD was more similar to BM. Significant CNAs in PMD/EMD included gains affecting 1q, 8q (MYC), 20p and 20q, with focal gains at 1q21, 6p25 and 18q21. Deletions frequently affected 4p, 8p and 17p (TP53), and numerous focal deletions were detailed. While there was only a trend to difference in 1p/1q between BM and IMD, gain/amp1q was observed in 38% BM, 56% PMD, 78% of EMD overall, and in 83% of post-therapy EMD samples.
Prevalent mutations in PMD/EMD involved KRAS, NRAS and TP53, with KRAS being enriched in PMD (38%), while EMD was enriched for NRAS (40%) and TP53 (33%). PMD/EMD were also enriched for loss-of-function mutations (p=0.007). In contrast, there was no significant difference in any of KRAS/NRAS/BRAF or TP53 between BM and IMD. Dndscv analysis revealed multiple genes positively selected as driver events in PMD/EMD, including epigenetic modifiers (KMT2A, KMT2C, TET2, SMARCA4), immune response (PRMD1, LTB), proliferation (STAT3, NF1, BTG1), adhesion (FAT1), and DNA-damage response (TP53, ATM, BRCA2).
Paired BM-PMD/EMD showed 2 patterns; some had increased mutational burden in PMD/EMD, while others had retained mutations but increased CNA. Within those with &gt; 1 PMD/EMD sample, diverging clones were defined by a variety of mutations (including KRAS, TP53, STAT3, IGF1R, DNMT3A), CNA affecting multiple chromosomes (including 8q; MYC, 17p; TP53), noncanonical translocations, and CNA consistent with chromothripsis.
In the WGS data, we note spatial divergence in mutational signature contribution and emergence of complex structural variants. Several patients with prior melphalan exposure have SBS99 evident in several biopsies, consistent with single cell expansion from MM cells surviving transplant and subsequently seeding in multiple sites.
Conclusions: PMD and EMD demonstrate multiple features of genomic complexity when compared with BM-based myeloma, which include emerging copy number aberrations, mutational burden and complex structural variants. EMD are more complex in general than PMD, while IMD shows only trends to increased genomic aberration compared with BM. Ongoing analyses include an expansion in WGS samples, and correlation of genomic features with clinical response to therapy.