In multiple myeloma malignant plasma cells expand within the bone marrow. Since this site is well-perfused, a rapid dissemination of “fitter” clones may be anticipated. However, an imbalanced distribution of multiple myeloma is frequently observed in medical imaging. Here, we perform multi-region sequencing, including iliac crest and radiology-guided focal lesion specimens from 51 patients to gain insight into the spatial clonal architecture. We demonstrate spatial genomic heterogeneity in more than 75% of patients, including inactivation of CDKN2C and TP53, and mutations affecting mitogen-activated protein kinase genes. We show that the extent of spatial heterogeneity is positively associated with the size of biopsied focal lesions consistent with regional outgrowth of advanced clones. The results support a model for multiple myeloma progression with clonal sweeps in the early phase and regional evolution in advanced disease. We suggest that multi-region investigations are critical to understanding intra-patient heterogeneity and the evolutionary processes in multiple myeloma.

Key Points: Hits in driver genes and bi-allelic events affecting tumor suppressors increase apoptosis resistance and proliferation rate–driving relapse. Excessive biallelic inactivation of tumor suppressors in high-risk cases highlights the need for TP53-independent therapeutic approaches.

The impact of MYC locus aberrations on the outcome of multiple myeloma (MM) patients is still a matter of debate. The aim of this study was to further investigate their influence on the survival of MM patients treated with high-dose chemotherapy. Our data suggest that the favorable prognosis factor

It is understood that cancer is a clonal disease initiated by a single cell, and that metastasis, which is the spread of cancer from the primary site, is also initiated by a single cell. The seemingly natural capability of cancer to adapt dynamically in a Darwinian manner is a primary reason for therapeutic failures. Survival advantages may be induced by cancer therapies and also occur as a result of inherent cell and microenvironmental factors. The selected "more fit" clones outmatch their competition and then become dominant in the tumor via propagation of progeny. This clonal expansion leads to relapse, therapeutic resistance and eventually death. The goal of this study is to develop and demonstrate a more detailed clonality approach by utilizing integrative genomics.

Gene expression profiling (GEP) has significantly contributed to the elucidation of the molecular heterogeneity of multiple myeloma plasma cells (MMPC) and only recently it has been recommended for risk stratification. Prior to GEP MMPC need to be enriched resulting in an inability to immediately freeze bone marrow aspirates or use RNA stabilization reagents. As a result in multi-center MM trials sample processing delay due to shipping may be an important confounder of molecular analyses and risk stratification based on GEP data. We compared GEP data of 145 in-house and 246 shipped samples and detected 3301 down-regulated and 3501 up-regulated genes in shipped samples. For 3994 genes we confirmed differential expression in an independent set of 85 in-house and 97 shipped samples. Differentially expressed genes were enriched in processes like ribosome biogenesis, cell cycle, and apoptosis. Among GEP based risk predictors the IFM-15 seemed to underestimate high risk in shipped samples, whereas the GEP70 and the EMC-92 gene signatures were more robust. In order to provide a tool to assess the “shipping effect” in public repositories, we generated a 17-gene predictor for shipped samples with a 10-fold cross validation error rate of 0.06 for the training set and an error rate of 0.15 for the validation set. Sample processing delay significantly influences GEP of MMPC, implying it should be avoided if samples were used for risk stratification.