Genomic landscape and chronological reconstruction of driver events in multiple myeloma

The genomic landscape of Vk*MYC myeloma highlights shared pathways of transformation between mice and humans

Multiple myeloma (MM) is a heterogeneous disease characterized by frequent MYC translocations. Sporadic MYC activation in the germinal center of genetically engineered Vk*MYC mice is sufficient to induce plasma cell tumors in which a variety of secondary mutations are spontaneously acquired and selected over time. Analysis of 119 Vk*MYC myeloma reveals recurrent copy number alterations, structural variations, chromothripsis, driver mutations, apolipoprotein B mRNA-editing enzyme, catalytic polypeptide (APOBEC) mutational activity, and a progressive decrease in immunoglobulin transcription that inversely correlates with proliferation. Moreover, we identify frequent insertional mutagenesis by endogenous retro-elements as a murine specific mechanism to activate NF-kB and IL6 signaling pathways shared with human MM. Despite the increased genomic complexity associated with progression, advanced tumors remain dependent on MYC. In summary, here we credential the Vk*MYC mouse as a unique resource to explore MM genomic evolution and describe a fully annotated collection of diverse and immortalized murine MM tumors.

The genomic profiling of high-risk smoldering myeloma patients treated with an intensive strategy unveils potential markers of resistance and progression

Smoldering multiple myeloma (SMM) precedes multiple myeloma (MM). The risk of progression of SMM patients is not uniform, thus different progression-risk models have been developed, although they are mainly based on clinical parameters. Recently, genomic predictors of progression have been defined for untreated SMM. However, the usefulness of such markers in the context of clinical trials evaluating upfront treatment in high-risk SMM (HR SMM) has not been explored yet, precluding the identification of baseline genomic alterations leading to drug resistance. For this reason, we carried out next-generation sequencing and fluorescent in-situ hybridization studies on 57 HR and ultra-high risk (UHR) SMM patients treated in the phase II GEM-CESAR clinical trial (NCT02415413). DIS3, FAM46C, and FGFR3 mutations, as well as t(4;14) and 1q alterations, were enriched in HR SMM. TRAF3 mutations were specifically associated with UHR SMM but identified cases with improved outcomes. Importantly, novel potential predictors of treatment resistance were identified: NRAS mutations and the co-occurrence of t(4;14) plus FGFR3 mutations were associated with an increased risk of biological progression. In conclusion, we have carried out for the first time a molecular characterization of HR SMM patients treated with an intensive regimen, identifying genomic predictors of poor outcomes in this setting.

Fluorescence in situ hybridization reveals the evolutionary biology of minor clone of gain/amp(1q) in multiple myeloma

Growing evidence suggests that gain or amplification [gain/amp(1q)] accumulates during disease progression of multiple myeloma (MM). Previous investigations have indicated that small gain/amp(1q) subclones present at the time of diagnosis may evolve into dominant clones upon MM relapse. However, the influence of a minor clone of gain/amp(1q) on MM survival, as well as the correlation between different clonal sizes of gain/amp(1q) and the chromosomal instability (CIN) of MM, remains poorly understood. In this study, we analyzed fluorescence in situ hybridization (FISH) results of 998 newly diagnosed MM (NDMM) patients. 513 patients were detected with gain/amp(1q) at diagnosis. Among these 513 patients, 55 had a minor clone (≤20%) of gain/amp(1q). Patients with a minor clone of gain/amp(1q) displayed similar survival outcomes compared to those without gain/amp(1q). Further analysis demonstrated patients with a minor clone of gain/amp(1q) exhibited a clonal architecture similar to those without gain/amp(1q). Lastly, our results showed a significant increase in the clonal size of the minor clone of gain/amp(1q), frequently observed in MM. These findings suggested that a minor clone of gain/amp(1q) might represent an earlier stage in the pathogenesis of gain/amp(1q) and propose a "two-step" process in the clonal size changes of gain/amp(1q) in MM.

Genomic Classification and Individualized Prognosis in Multiple Myeloma

PURPOSE

Outcomes for patients with newly diagnosed multiple myeloma (NDMM) are heterogenous, with overall survival (OS) ranging from months to over 10 years.

METHODS

To decipher and predict the molecular and clinical heterogeneity of NDMM, we assembled a series of 1,933 patients with available clinical, genomic, and therapeutic data.

RESULTS

Leveraging a comprehensive catalog of genomic drivers, we identified 12 groups, expanding on previous gene expression-based molecular classifications. To build a model predicting individualized risk in NDMM (IRMMa), we integrated clinical, genomic, and treatment variables. To correct for time-dependent variables, including high-dose melphalan followed by autologous stem-cell transplantation (HDM-ASCT), and maintenance therapy, a multi-state model was designed. The IRMMa model accuracy was significantly higher than all comparator prognostic models, with a c-index for OS of 0.726, compared with International Staging System (ISS; 0.61), revised-ISS (0.572), and R2-ISS (0.625). Integral to model accuracy was 20 genomic features, including 1q21 gain/amp, del 1p, TP53 loss, NSD2 translocations, APOBEC mutational signatures, and copy-number signatures (reflecting the complex structural variant chromothripsis). IRMMa accuracy and superiority compared with other prognostic models were validated on 256 patients enrolled in the GMMG-HD6 (ClinicalTrials.gov identifier: NCT02495922) clinical trial. Individualized patient risks were significantly affected across the 12 genomic groups by different treatment strategies (ie, treatment variance), which was used to identify patients for whom HDM-ASCT is particularly effective versus patients for whom the impact is limited.

CONCLUSION

Integrating clinical, demographic, genomic, and therapeutic data, to our knowledge, we have developed the first individualized risk-prediction model enabling personally tailored therapeutic decisions for patients with NDMM.

Impact of Clonal Heterogeneity in Multiple Myeloma

Multiple myeloma is characterized by a highly heterogeneous disease distribution within the bone marrow-containing skeletal system. In this review, we introduce the molecular mechanisms underlying clonal heterogeneity and the spatio-temporal evolution of myeloma. We discuss the clinical impact of clonal heterogeneity, which is thought to be one of the biggest obstacles to overcome therapy resistance and to achieve cure.

Chromothripsis in hematologic malignancies

Chromotrypsis, a phenomenon resulting from catastrophic mitotic errors and genomic instability, is defined by the occurrence of multiple DNA double-strand breaks in one or more chromosomes, subsequently subject to error-prone repair mechanisms. This unique process results in extensive rearrangements in the affected chromosomes, leading to loss of tumor suppressor function, the creation of fusion genes, and/or activation of oncogenes. The importance of chromothripsis in cancer, especially in the field of hematologic disorders, underscores the intricate interplay between genomic instability and the genesis of alterations that contribute to cancer. This accentuates the critical need to unravel these complex processes for the targeted development of specific therapeutic interventions. This review delves into the analysis of chromothripsis cases in various hematologic diseases, such as leukemia, lymphoma, and myeloma, with the aim of unveiling its profound impact on patient prognosis. Furthermore, the study explores the intricate molecular mechanisms underlying chromothripsis and investigates its consequences.

Approach to High-Risk Multiple Myeloma

Improving the outcome of high-risk myeloma (HRMM) is a key therapeutic aim for the next decade. To achieve this aim, it is necessary to understand in detail the genetic drivers underlying this clinical behavior and to target its biology therapeutically. Advances have already been made, with a focus on consensus guidance and the application of novel immunotherapeutic approaches. Cases of HRMM are likely to have impaired prognosis even with novel strategies. However, if disease eradication and minimal disease states are achieved, then cure may be possible.

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

The complexity of Multiple Myeloma (MM) is driven by several genomic aberrations, interacting with disease-related and/or -unrelated factors and conditioning patients' clinical outcome. Patient's prognosis is hardly predictable, as commonly employed MM risk models do not precisely partition high- from low-risk patients, preventing the reliable recognition of early relapsing/refractory patients. By a dimensionality reduction approach, here we dissect the genomic landscape of a large cohort of newly diagnosed MM patients, modelling all the possible interactions between any MM chromosomal alterations. We highlight the presence of a distinguished cluster of patients in the low-dimensionality space, with unfavorable clinical behavior, whose biology was driven by the co-occurrence of chromosomes 1q CN gain and 13 CN loss. Presence or absence of these alterations define MM patients overexpressing either CCND2 or CCND1, fostering the implementation of biology-based patients' classification models to describe the different MM clinical behaviors.

Impact of Rare Structural Variant Events in Newly Diagnosed Multiple Myeloma

PURPOSE

Whole-genome sequencing (WGS) of patients with newly diagnosed multiple myeloma (NDMM) has shown recurrent structural variant (SV) involvement in distinct regions of the genome (i.e., hotspots) and causing recurrent copy-number alterations. Together with canonical immunoglobulin translocations, these SVs are recognized as "recurrent SVs." More than half of SVs were not involved in recurrent events. The significance of these "rare SVs" has not been previously examined.

EXPERIMENTAL DESIGN

In this study, we utilize 752 WGS and 591 RNA sequencing data from patients with NDMM to determine the role of rare SVs in myeloma pathogenesis.

RESULTS

Ninety-four percent of patients harbored at least one rare SV event. Rare SVs showed an SV class-specific enrichment within genes and superenhancers associated with outlier gene expression. Furthermore, known myeloma driver genes recurrently impacted by point mutations were dysregulated by rare SVs.

CONCLUSIONS

Overall, we demonstrate the association of rare SVs with aberrant gene expression supporting a potential driver role in myeloma pathogenesis.

Mechanisms of the role of proto-oncogene activation in promoting malignant transformation of mature B cells

Malignant tumors continue to pose a significant threat to human life and safety and their development is primarily due to the activation of proto-oncogenes and the inactivation of suppressor genes. Among these, the activation of proto-oncogenes possesses greater potential to drive the malignant transformation of cells. Targeting oncogenes involved in the malignant transformation of tumor cells has provided a novel approach for the development of current antitumor drugs. Several preclinical and clinical studies have revealed that the development pathway of B cells, and the malignant transformation of mature B cells into tumors have been regulated by oncogenes and their metabolites. Therefore, summarizing the key oncogenes involved in the process of malignant transformation of mature B cells and elucidating the mechanisms of action in tumor development hold significant importance for the clinical treatment of malignant tumors.

Novel Insights into the Initiation, Evolution, and Progression of Multiple Myeloma by Multi-Omics Investigation

The evolutionary history of multiple myeloma (MM) includes malignant transformation, followed by progression to pre-malignant stages and overt malignancy, ultimately leading to more aggressive and resistant forms. Over the past decade, large effort has been made to identify the potential therapeutic targets in MM. However, MM remains largely incurable. Most patients experience multiple relapses and inevitably become refractory to treatment. Tumor-initiating cell populations are the postulated population, leading to the recurrent relapses in many hematological malignancies. Clonal evolution of tumor cells in MM has been identified along with the disease progression. As a consequence of different responses to the treatment of heterogeneous MM cell clones, the more aggressive populations survive and evolve. In addition, the tumor microenvironment is a complex ecosystem which plays multifaceted roles in supporting tumor cell evolution. Emerging multi-omics research at single-cell resolution permits an integrative and comprehensive profiling of the tumor cells and microenvironment, deepening the understanding of biological features of MM. In this review, we intend to discuss the novel insights into tumor cell initiation, clonal evolution, drug resistance, and tumor microenvironment in MM, as revealed by emerging multi-omics investigations. These data suggest a promising strategy to unravel the pivotal mechanisms of MM progression and enable the improvement in treatment, both holistically and precisely.

IGH::NSD2 Fusion Gene Transcript as Measurable Residual Disease Marker in Multiple Myeloma

Multiple myeloma (MM) is the second most common hematological malignancy. Approximately 15% of MM patients are affected by the t(4;14) translocation resulting in the IGH::NSD2 fusion transcript. Breakage occurs in three major breakpoint regions within the NSD2 gene (MB4-1, MB4-2, and MB4-3), where MB4-1 leads to the production of full-length protein, while truncated proteins are expressed in the other two cases. Measurable residual disease (MRD) has been conclusively established as a crucial prognostic factor in MM. The IGH::NSD2 fusion transcript can serve as a sensitive MRD marker. Using bone marrow (BM) and peripheral blood (PB) samples from 111 patients, we developed a highly sensitive quantitative real-time PCR (qPCR) and digital PCR (dPCR) system capable of detecting fusion mRNAs with a sensitivity of up to 1:100,000. PB samples exhibited sensitivity three orders of magnitude lower compared to BM samples. Patients with an MB4-2 breakpoint demonstrated significantly reduced overall survival (p = 0.003). Our novel method offers a simple and sensitive means for detecting MRD in a substantial proportion of MM patients. Monitoring may be carried out even from PB samples. The literature lacks consensus regarding survival outcomes among patients with different NSD2 breakpoints. Our data align with previous findings indicating that patients with the MB4-2 breakpoint type tend to exhibit unfavorable overall survival.

Clearing soluble MIC reverses the impaired function of natural killer cells from patients with multiple myeloma

BACKGROUND

Major histocompatibility complex (MHC) class I chain-related protein (MIC) is a stress-induced ligand released from multiple myeloma (MM) cells during progression, and soluble MIC impairs natural killer group 2D (NKG2D) activating receptor-mediated recognition and function of natural killer (NK) cells. However, whether clearing soluble MIC with a monoclonal antibody (mAb) can restore NK cell activity of MM patients remains undetermined.

METHODS

We analyzed The Cancer Genome Atlas (TCGA) Multiple Myeloma Research Foundation (MMRF) CoMMpass data set to examine the prognostic significance of MIC expression in MM. We examined the level of soluble MIC in paired peripheral blood (PB) and bone marrow (BM) plasma of patients with MM at diagnosis by ELISA. We evaluated the correlation between the level of soluble MIC and immunophenotype of NK cells from MM patients by multicolor flow cytometry. We also generated MIC-overexpressing MM cell line and characterized the cytotoxic function of patient NK cells in the presence of soluble MIC, and examined the impact of clearing soluble MIC with a humanized mAb (huB10G5).

RESULTS

We characterize the importance of MICA in MM by revealing the significantly better overall survival of patients with high MICA expression from TCGA MMRF CoMMpass data set. The level of soluble MICA is more highly elevated in MM than in precursor stages, and the concentration of soluble MICA is higher in BM plasma than in PB. The concentration of soluble MICA in BM was correlated with myeloma burden, while it was negatively correlated with the frequency of NKG2D+ NK cells in diagnostic BM aspirates of MM patients. Soluble MICA downregulated NKG2D expression and decreased cytotoxicity of MM patient NK cells ex vivo, which were reversed by a humanized soluble MIC-clearing mAb (huB10G5) with enhanced degranulation of NK cells.

CONCLUSIONS

Our findings indicate targeting soluble MIC with huB10G5 might be a viable therapeutic approach to promote NKG2D-dependent cellular immunotherapy outcome in MM.

Clinicopathological and genetic landscape of plasmablastic lymphoma in Taiwan

Plasmablastic lymphoma (PBL) is an aggressive large B-cell lymphoma with a terminal B-cell differentiation phenotype and is frequently associated with immunodeficiency. We aimed to investigate the clinicopathological and immunophenotypic features, genetic alterations, and mutational landscape of PBL in Taiwan. We retrospectively recruited 26 cases. Five (5/18; 28%) patients were HIV-positive and 21 (81%) presented extranodally. There were two morphological groups: one with purely monomorphic large cells (85%) and the other comprising large cells admixed with plasmacytic cells (15%). Phenotypically, the tumors expressed MYC (8/10; 80%), CD138 (20/26; 77%), and MUM1 (20/20; 100%), but not CD20 (n = 26; 0%). Fourteen (54%) cases were positive for EBV by in situ hybridization; the EBV-positive cases were more frequently HIV infected (p = 0.036), with extranodal presentation (p = 0.012) and CD79a expression (p = 0.012), but less frequent light chain restriction (p = 0.029). Using fluorescence in situ hybridization, we identified 13q14 deletion, MYC rearrangement, and CCND1 rearrangement in 74%, 30%, and 5% cases, respectively, without any cases having rearranged BCL6 or IGH::FGFR3 fusion. In the 15 cases with adequate tissue for whole exome sequencing, the most frequent recurrent mutations were STAT3 (40%), NRAS (27%), and KRAS (20%). In conclusion, most PBL cases in Taiwan were HIV-unrelated. Around half of the cases were positive for EBV, with distinct clinicopathological features. Deletion of chromosome 13q14 was frequent. The PBL cases in Taiwan showed recurrent mutations involving JAK-STAT, RAS-MAPK, epigenetic regulation, and NOTCH signaling pathways, findings similar to that from the West.

Discovery and Evaluation of C6-Substituted Pyrazolopyrimidine-Based Bisphosphonate Inhibitors of the Human Geranylgeranyl Pyrophosphate Synthase and Evaluation of Their Antitumor Efficacy in Multiple Myeloma, Pancreatic Ductal Adenocarcinoma, and Colorectal Cancer

Novel C6-substituted pyrazolo[3,4-d]pyrimidine- and C2-substituted purine-based bisphosphonate (C6-PyraP-BP and C2-Pur-BP, respectively) inhibitors of the human geranylgeranyl pyrophosphate synthase (hGGPPS) were designed and evaluated for their ability to block the proliferation of multiple myeloma (MM), pancreatic ductal adenocarcinoma (PDAC), and colorectal cancer (CRC) cells. Pyrazolo[3,4-d]pyrimidine analogs were identified that induce selective intracellular target engagement leading to apoptosis and downregulate the prenylation of Rap-1A in MM, PDAC, and CRC cells. The C6-PyraP-BP inhibitor RB-07-16 was found to exhibit antitumor efficacy in xenograft mouse models of MM and PDAC, significantly reducing tumor growth without substantially increasing liver enzymes or causing significant histopathologic damage, usually associated with hepatotoxicity. RB-07-16 is a metabolically stable compound in cross-species liver microsomes, does not inhibit key CYP 450 enzymes, and exhibits good systemic circulation in rat. Collectively, the current studies provide encouraging support for further optimization of the pyrazolo[3,4-d]pyrimidine-based GGPPS inhibitors as potential human therapeutics for various cancers.

Genomic and immune determinants of resistance to anti-CD38 monoclonal antibody-based therapy in relapsed refractory multiple myeloma

Anti-CD38 antibody therapies have transformed multiple myeloma (MM) treatment. However, a large fraction of patients inevitably relapses. To understand this, we investigated 32 relapsed MM patients treated with daratumumab, lenalidomide, and dexamethasone (Dara-Rd; NCT03848676 ). Whole genome sequencing (WGS) before and after treatment pinpointed genomic drivers associated with early progression, including RPL5 loss and APOBEC mutagenesis. Flow cytometry on 202 blood samples, collected every three months until progression for 31 patients, revealed distinct immune changes significantly impacting clinical outcomes. Progressing patients exhibited significant depletion of CD38+ NK cells, persistence of T cell exhaustion, and reduced depletion of T-reg cells over time. These findings underscore the influence of immune composition and daratumumab-induced immune changes in promoting MM resistance. Integrating genomics and flow cytometry unveiled associations between adverse genomic features and immune patterns. Overall, this study sheds light on the intricate interplay between genomic complexity and the immune microenvironment driving resistance to Dara-Rd.

The successive emergence of ERVL-MaLRs in primates

Although the ERVL-mammalian-apparent LTR retrotransposons (MaLRs) are the fourth largest family of transposable elements in the human genome, their evolutionary history and relationship have not been thoroughly studied. In this study, through RepeatMasker annotations of some representative species and construction of phylogenetic tree by sequence similarity, all primate-specific MaLR members are found to descend from MLT1A1 retrotransposon. Comparative genomic analysis, transposition-in-transposition inference, and sequence feature comparisons consistently show that each MaLR member evolved from its predecessor successively and had a limited activity period during primate evolution. Accordingly, a novel MaLR member was discovered as successor of MSTB1 in Tarsiiformes. At last, the identification of candidate precursor and intermediate THE1A elements provides further evidence for the previously proposed arms race model between ZNF430/ZNF100 and THE1B/THE1A. Taken together, this study sheds light on the evolutionary history of MaLRs and can serve as a foundation for future research on their interactions with zinc finger genes, gene regulation, and human health implications.

High-throughput electron tomography identifies centriole over-elongation as an early event in plasma cell disorders

Plasma cell disorders are clonal outgrowths of pre-malignant or malignant plasma cells, characterized by extensive chromosomal aberrations. Centrosome abnormalities are a major driver of chromosomal instability in cancer but their origin, incidence, and composition in primary tumor cells is poorly understood. Using cutting-edge, semi-automated high-throughput electron tomography, we characterized at nanoscale 1386 centrioles in CD138pos plasma cells from eight healthy donors and 21 patients with plasma cell disorders, and 722 centrioles from different control populations. In plasma cells from healthy individuals, over-elongated centrioles accumulated with age. In plasma cell disorders, centriole over-elongation was notably frequent in early, pre-malignant disease stages, became less pronounced in overt multiple myeloma, and almost entirely disappeared in aggressive plasma cell leukemia. Centrioles in other types of patient-derived B cell neoplasms showed no over-elongation. In contrast to current belief, centriole length appears to be highly variable in long-lived, healthy plasma cells, and over-elongation and structural aberrations are common in this cell type. Our data suggest that structural centrosome aberrations accumulate with age in healthy CD138pos plasma cells and may thus play an important role in early aneuploidization as an oncogenic driver in plasma cell disorders.

Genomic and immune signatures predict clinical outcome in newly diagnosed multiple myeloma treated with immunotherapy regimens

Despite improving outcomes, 40% of patients with newly diagnosed multiple myeloma treated with regimens containing daratumumab, a CD38-targeted monoclonal antibody, progress prematurely. By integrating tumor whole-genome and microenvironment single-cell RNA sequencing from upfront phase 2 trials using carfilzomib, lenalidomide and dexamethasone with daratumumab ( NCT03290950 ), we show how distinct genomic drivers including high APOBEC mutational activity, IKZF3 and RPL5 deletions and 8q gain affect clinical outcomes. Furthermore, evaluation of paired bone marrow profiles, taken before and after eight cycles of carfilzomib, lenalidomide and dexamethasone with daratumumab, shows that numbers of natural killer cells before treatment, high T cell receptor diversity before treatment, the disappearance of sustained immune activation (that is, B cells and T cells) and monocyte expansion over time are all predictive of sustained minimal residual disease negativity. Overall, this study provides strong evidence of a complex interplay between tumor cells and the immune microenvironment that is predictive of clinical outcome and depth of treatment response in patients with newly diagnosed multiple myeloma treated with highly effective combinations containing anti-CD38 antibodies.

Gene interaction network analysis in multiple myeloma detects complex immune dysregulation associated with shorter survival

The plasma cell cancer multiple myeloma (MM) varies significantly in genomic characteristics, response to therapy, and long-term prognosis. To investigate global interactions in MM, we combined a known protein interaction network with a large clinically annotated MM dataset. We hypothesized that an unbiased network analysis method based on large-scale similarities in gene expression, copy number aberration, and protein interactions may provide novel biological insights. Applying a novel measure of network robustness, Ollivier-Ricci Curvature, we examined patterns in the RNA-Seq gene expression and CNA data and how they relate to clinical outcomes. Hierarchical clustering using ORC differentiated high-risk subtypes with low progression free survival. Differential gene expression analysis defined 118 genes with significantly aberrant expression. These genes, while not previously associated with MM, were associated with DNA repair, apoptosis, and the immune system. Univariate analysis identified 8/118 to be prognostic genes; all associated with the immune system. A network topology analysis identified both hub and bridge genes which connect known genes of biological significance of MM. Taken together, gene interaction network analysis in MM uses a novel method of global assessment to demonstrate complex immune dysregulation associated with shorter survival.

Novel Immunotherapies and Combinations: The Future Landscape of Multiple Myeloma Treatment

In multiple myeloma impressive outcomes have improved with the introduction of new therapeutic approaches, mainly those including naked monoclonal antibodies such as daratumumab and isatuximab. However, moving to earlier lines of therapy with effective anti-myeloma drugs led to an increase in the number of patients who developed multi-refractoriness to them early on. Currently, triple- or multi-refractory MM represents an unmet medical need, and their management remains a complicated challenge. The recent approval of new immunotherapeutic approaches such as conjugated monoclonal antibodies, bispecific antibodies, and CAR T cells could be a turning point for these heavily pretreated patients. Nevertheless, several issues regarding their use are unsolved, such as how to select patients for each strategy or how to sequence these therapies within the MM therapeutic landscape. Here we provide an overview of the most recent data about approved conjugated monoclonal antibody belantamab, mafodotin, bispecific antibody teclistamab, and other promising compounds under development, mainly focusing on the ongoing clinical trials with monoclonal antibody combination approaches in advanced and earlier phases of MM treatment.

Resolving therapy resistance mechanisms in multiple myeloma by multiomics subclone analysis

Intratumor heterogeneity as a clinical challenge becomes most evident after several treatment lines, when multidrug-resistant subclones accumulate. To address this challenge, the characterization of resistance mechanisms at the subclonal level is key to identify common vulnerabilities. In this study, we integrate whole-genome sequencing, single-cell (sc) transcriptomics (scRNA sequencing), and chromatin accessibility (scATAC sequencing) together with mitochondrial DNA mutations to define subclonal architecture and evolution for longitudinal samples from 15 patients with relapsed or refractory multiple myeloma. We assess transcriptomic and epigenomic changes to resolve the multifactorial nature of therapy resistance and relate it to the parallel occurrence of different mechanisms: (1) preexisting epigenetic profiles of subclones associated with survival advantages, (2) converging phenotypic adaptation of genetically distinct subclones, and (3) subclone-specific interactions of myeloma and bone marrow microenvironment cells. Our study showcases how an integrative multiomics analysis can be applied to track and characterize distinct multidrug-resistant subclones over time for the identification of molecular targets against them.

Clonal tracking in cancer and metastasis

The eradication of many cancers has proven challenging due to the presence of functionally and genetically heterogeneous clones maintained by rare cancer stem cells (CSCs), which contribute to disease progression, treatment refractoriness, and late relapse. The characterization of functional CSC activity has necessitated the development of modern clonal tracking strategies. This review describes viral-based and CRISPR-Cas9-based cellular barcoding, lineage tracing, and imaging-based approaches. DNA-based cellular barcoding technology is emerging as a powerful and robust strategy that has been widely applied to in vitro and in vivo model systems, including patient-derived xenograft models. This review also highlights the potential of these methods for use in the clinical and drug discovery contexts and discusses the important insights gained from such approaches.

Aiming for the cure in myeloma: Putting our best foot forward

Frontline therapy for multiple myeloma (MM) is evolving to include novel combinations that can achieve unprecedented deep response rates. Several treatment strategies exist, varying in induction regimen composition, use of transplant and or consolidation and maintenance. In this sea of different treatment permutations, the overarching theme is the powerful prognostic factors of disease risk and achievement of minimal residual disease (MRD) negativity. MM has significant inter-patient variability that requires treatment to be individualized. Risk-adapted and response-adapted strategies which are increasingly being explored to define the extent and duration of therapy, and eventually aim for functional curability. In addition, with T-cell redirection therapies rapidly revolutionizing myeloma treatments, the current standard of care for myeloma will change. This review analyzes the current relevant literature in upfront therapy for fit myeloma patients and provides suggestions for treatment approach while novel clinical trials are maturing.

Scientific techniques in adolescent and young adult classic Hodgkin lymphoma

Understanding the tumor microenvironment and genomic landscape is crucial for better prediction of treatment outcomes and developing novel therapies in Hodgkin lymphoma (HL). Recent advancements in genomics have enabled researchers to gain deeper insights into the genomic characteristics of HL at both single-cell resolution and the whole genome level. The use of noninvasive methods such as liquid biopsies and formalin-fixed paraffin-embedded-based imaging techniques has expanded the possibilities of applying cutting-edge analyses to routine clinically available samples. Collaborative efforts between adult and pediatric group are imperative to translate novel findings into routine patient care.

Immunophenotypic correlates of sustained MRD negativity in patients with multiple myeloma

The role of the immune microenvironment in maintaining disease remission in patients with multiple myeloma (MM) is not well understood. In this study, we comprehensively profile the immune system in patients with newly diagnosed MM receiving continuous lenalidomide maintenance therapy with the aim of discovering correlates of long-term treatment response. Leveraging single-cell RNA sequencing and T cell receptor β sequencing of the peripheral blood and CyTOF mass cytometry of the bone marrow, we longitudinally characterize the immune landscape in 23 patients before and one year after lenalidomide exposure. We compare patients achieving sustained minimal residual disease (MRD) negativity to patients who never achieved or were unable to maintain MRD negativity. We observe that the composition of the immune microenvironment in both the blood and the marrow varied substantially according to both MRD negative status and history of autologous stem cell transplant, supporting the hypothesis that the immune microenvironment influences the depth and duration of treatment response.

Intratumor microbiota: Occult participants in the microenvironment of multiple myeloma

More recently, microbiota was detected in several tumorous tissues including multiple myeloma (MM), but the roles of which is still under-studied as paucity of research on tumor biology. Moreover, we also detected the presence of microbiota in the bone marrow of patients with MM by 2bRAD-M sequencing technology, which is an incurable hematological malignancy characterized by accumulation of abnormal plasma cells in the bone marrow. However, the roles of intratumor microbiota in tumor disease remains poorly understood. In this review, we critically reviewed recent literature about microbiota in the tumorigenesis and progression of MM. Importantly, we proposed that the emergence of microbiota in the microenvironment of multiple myeloma may be attributed to microbial dysbiosis and impaired intestinal barrier, due to the increased prevalence of MM in patients with obesity and diabetes, of which the characteristic phenotype is gut microbial dysbiosis and impaired intestinal barrier. When the intestinal barrier is damaged, dysbiotic microbiota and their metabolites, as well as dysregulated immune cells, may participate in the reshaping of the local immune microenvironment, and play pivotal roles in the tumorigenesis and development of multiple myeloma, probably by migrating to the bone marrow microenvironment from intestine. We also discuss the emerging microbiological manipulation strategies to improve long-term outcomes of MM, as well as the prospective of the state-of-the-art techniques to advance our knowledge about the biological implication in the microbiome in MM.

Mechanisms of antigen escape from BCMA- or GPRC5D-targeted immunotherapies in multiple myeloma

B cell maturation antigen (BCMA) target loss is considered to be a rare event that mediates multiple myeloma (MM) resistance to anti-BCMA chimeric antigen receptor T cell (CAR T) or bispecific T cell engager (TCE) therapies. Emerging data report that downregulation of G-protein-coupled receptor family C group 5 member D (GPRC5D) protein often occurs at relapse after anti-GPRC5D CAR T therapy. To examine the tumor-intrinsic factors that promote MM antigen escape, we performed combined bulk and single-cell whole-genome sequencing and copy number variation analysis of 30 patients treated with anti-BCMA and/or anti-GPRC5D CAR T/TCE therapy. In two cases, MM relapse post-TCE/CAR T therapy was driven by BCMA-negative clones harboring focal biallelic deletions at the TNFRSF17 locus at relapse or by selective expansion of pre-existing subclones with biallelic TNFRSF17 loss. In another five cases of relapse, newly detected, nontruncating, missense mutations or in-frame deletions in the extracellular domain of BCMA negated the efficacies of anti-BCMA TCE therapies, despite detectable surface BCMA protein expression. In the present study, we also report four cases of MM relapse with biallelic mutations of GPRC5D after anti-GPRC5D TCE therapy, including two cases with convergent evolution where multiple subclones lost GPRC5D through somatic events. Immunoselection of BCMA- or GPRC5D-negative or mutant clones is an important tumor-intrinsic driver of relapse post-targeted therapies. Mutational events on BCMA confer distinct sensitivities toward different anti-BCMA therapies, underscoring the importance of considering the tumor antigen landscape for optimal design and selection of targeted immunotherapies in MM.

Systematic characterization of therapeutic vulnerabilities in Multiple Myeloma with Amp1q reveals increased sensitivity to the combination of MCL1 and PI3K inhibitors

The development of targeted therapy for patients with Multiple Myeloma (MM) is hampered by the low frequency of actionable genetic abnormalities. Gain or amplification of chr1q (Amp1q) is the most frequent arm-level copy number gain in patients with MM, and it is associated with higher risk of progression and death despite recent advances in therapeutics. Thus, developing targeted therapy for patients with MM and Amp1q stands to benefit a large portion of patients in need of more effective management. Here, we employed large-scale dependency screens and drug screens to systematically characterize the therapeutic vulnerabilities of MM with Amp1q and showed increased sensitivity to the combination of MCL1 and PI3K inhibitors. Using single-cell RNA sequencing, we compared subclones with and without Amp1q within the same patient tumors and showed that Amp1q is associated with higher levels of MCL1 and the PI3K pathway. Furthermore, by isolating isogenic clones with different copy number for part of the chr1q arm, we showed increased sensitivity to MCL1 and PI3K inhibitors with arm-level gain. Lastly, we demonstrated synergy between MCL1 and PI3K inhibitors and dissected their mechanism of action in MM with Amp1q.

Molecular genetic aberrations in the pathogenesis of multiple myeloma

Multiple myeloma (MM) is the second most common form of blood cancer characterized by clonal expansion of malignant plasma cells within the bone marrow. MM is a complex, progressive, and highly heterogeneous malignancy, which occurs via a multistep transformation process involving primary and secondary oncogenic events. Recent advances in molecular techniques have further expanded our understanding of the mutational landscape, clonal composition, and dynamic evolution patterns of MM. The first part of this review describes the key oncogenic events involved in the initiation and progression of MM, together with their prognostic impact. The latter part highlights the most prominent findings concerning genomic aberrations promoted by gene expression profiling (GEP) and next-generation sequencing (NGS) in MM. This review provides a concise understanding of the molecular pathogenesis of the MM genome and the importance of adopting emerging molecular technology in future clinical management of MM.

The Vk*MYC Mouse Model recapitulates human multiple myeloma evolution and genomic diversity

Despite advancements in profiling multiple myeloma (MM) and its precursor conditions, there is limited information on mechanisms underlying disease progression. Clincal efforts designed to deconvolute such mechanisms are challenged by the long lead time between monoclonal gammopathy and its transformation to MM. MM mouse models represent an opportunity to overcome this temporal limitation. Here, we profile the genomic landscape of 118 genetically engineered Vk*MYC MM and reveal that it recapitulates the genomic heterogenenity and life history of human MM. We observed recurrent copy number alterations, structural variations, chromothripsis, driver mutations, APOBEC mutational activity, and a progressive decrease in immunoglobulin transcription that inversely correlates with proliferation. Moreover, we identified frequent insertional mutagenesis by endogenous retro-elements as a murine specific mechanism to activate NF-kB and IL6 signaling pathways shared with human MM. Despite the increased genomic complexity associated with progression, advanced tumors remain dependent on MYC expression, that drives the progression of monoclonal gammopathy to MM.

Single-cell profiling in multiple myeloma: insights, problems, and promises

In a short time, single-cell platforms have become the norm in many fields of research, including multiple myeloma (MM). In fact, the large amount of cellular heterogeneity in MM makes single-cell platforms particularly attractive because bulk assessments can miss valuable information about cellular subpopulations and cell-to-cell interactions. The decreasing cost and increasing accessibility of single-cell platform, combined with breakthroughs in obtaining multiomics data for the same cell and innovative computational programs for analyzing data, have allowed single-cell studies to make important insights into MM pathogenesis; yet, there is still much to be done. In this review, we will first focus on the types of single-cell profiling and the considerations for designing a single-cell profiling experiment. Then, we will discuss what have learned from single-cell profiling about myeloma clonal evolution, transcriptional reprogramming, and drug resistance, and about the MM microenvironment during precursor and advanced disease.

Minimal Residual Disease in Multiple Myeloma: Past, Present, and Future

Responses to treatment have improved over the last decades for patients with multiple myeloma. This is a consequence of the introduction of new drugs that have been successfully combined in different clinical contexts: newly diagnosed, transplant-eligible or ineligible patients, as well as in the relapsed/refractory setting. However, a great proportion of patients continue to relapse, even those achieving complete response, which underlines the need for updated response criteria. In 2014, the international myeloma working group established new levels of response, prompting the evaluation of minimal residual disease (MRD) for those patients already in complete or stringent complete response as defined by conventional serological assessments: the absence of tumor plasma cells in 100,000 total cells or more define molecular and immunophenotypic responses by next-generation sequencing and flow cytometry, respectively. In this review, we describe all the potential methods that may be used for MRD detection based on the evidence found in the literature, paying special attention to their advantages and pitfalls from a critical perspective.

Pathogenesis, clinical characteristics and personalized managements of multiple myeloma with chromosome 1 abnormalities

Multiple myeloma (MM) is a biologically heterogeneous malignancy defined by the proliferation of monoclonal plasma cells. Despite the tremendous advancement in MM treatment over the past decades, relapse remains a major problem which is inevitable for most patients. In particular, a partial of patients with early relapse and poor outcomes are classified as a high-risk group. Apart from the clinical stage, genetic aberrations are now recognized as important prognostic factors for identifying high-risk patients. Chromosome 1 abnormalities (C1As), particularly 1q21 gain or amplification, have been identified as common genetic aberrations in patients with MM and are often considered unfavorable prognostic markers for progression-free survival and overall survival. However, more effective therapeutic approaches are still needed to overcome the negative impact of C1As. Therefore, we summarize the prevalence, pathogenesis, clinical significance and present therapeutic condition of C1As in MM, and attempt to conclude the precise and personalized management for patients with C1As.

Tracking the evolution of therapy-related myeloid neoplasms using chemotherapy signatures

Patients treated with cytotoxic therapies, including autologous stem cell transplantation, are at risk for developing therapy-related myeloid neoplasms (tMN). Preleukemic clones (ie, clonal hematopoiesis [CH]) are detectable years before the development of these aggressive malignancies, although the genomic events leading to transformation and expansion are not well defined. Here, by leveraging distinctive chemotherapy-associated mutational signatures from whole-genome sequencing data and targeted sequencing of prechemotherapy samples, we reconstructed the evolutionary life-history of 39 therapy-related myeloid malignancies. A dichotomy was revealed, in which neoplasms with evidence of chemotherapy-induced mutagenesis from platinum and melphalan were hypermutated and enriched for complex structural variants (ie, chromothripsis), whereas neoplasms with nonmutagenic chemotherapy exposures were genomically similar to de novo acute myeloid leukemia. Using chemotherapy-associated mutational signatures as temporal barcodes linked to discrete clinical exposure in each patient's life, we estimated that several complex events and genomic drivers were acquired after chemotherapy was administered. For patients with prior multiple myeloma who were treated with high-dose melphalan and autologous stem cell transplantation, we demonstrate that tMN can develop from either a reinfused CH clone that escapes melphalan exposure and is selected after reinfusion, or from TP53-mutant CH that survives direct myeloablative conditioning and acquires melphalan-induced DNA damage. Overall, we revealed a novel mode of tMN progression that is not reliant on direct mutagenesis or even exposure to chemotherapy. Conversely, for tMN that evolve under the influence of chemotherapy-induced mutagenesis, distinct chemotherapies not only select preexisting CH but also promote the acquisition of recurrent genomic drivers.

Molecular Evolution of Classic Hodgkin Lymphoma Revealed Through Whole-Genome Sequencing of Hodgkin and Reed Sternberg Cells

The rarity of malignant Hodgkin and Reed Sternberg (HRS) cells in classic Hodgkin lymphoma (cHL) limits the ability to study the genomics of cHL. To circumvent this, our group has previously optimized fluorescence-activated cell sorting to purify HRS cells. Using this approach, we now report the whole-genome sequencing landscape of HRS cells and reconstruct the chronology and likely etiology of pathogenic events leading to cHL. We identified alterations in driver genes not previously described in cHL, APOBEC mutational activity, and the presence of complex structural variants including chromothripsis. We found that high ploidy in cHL is often acquired through multiple, independent chromosomal gains events including whole-genome duplication. Evolutionary timing analyses revealed that structural variants enriched for RAG motifs, driver mutations in B2M, BCL7A, GNA13, and PTPN1, and the onset of AID-driven mutagenesis usually preceded large chromosomal gains. This study provides a temporal reconstruction of cHL pathogenesis.

SIGNIFICANCE

Previous studies in cHL were limited to coding sequences and therefore not able to comprehensively decipher the tumor complexity. Here, leveraging cHL whole-genome characterization, we identify driver events and reconstruct the tumor evolution, finding that structural variants, driver mutations, and AID mutagenesis precede chromosomal gains. This article is highlighted in the In This Issue feature, p. 171.

[Latest Findings on the Mechanism of the Interaction Between Multiple Myeloma Cells and Bone Marrow Microenvironment]

Multiple myeloma (MM) is a hematologic malignancy of terminally differentiated plasma cells. The mechanisms of the pathogenesis and progression of MM include genetic abnormalities of the MM cells and the interaction between MM cells and bone marrow microenvironment (BMME). MM cells start malignant proliferation in BMME and contribute to the pathogenesis and progression of MM through direct or indirect interactions between cells and the extracellular matrix. Exploring the mechanism of interaction between MM cells and the microenvironment is crucial to improving our understanding of the pathogenesis and progression of MM and early diagnosis and treatment. In addition, the metabolic reprogramming of tumors is one of the key issues of oncology research. Herein, we summarized published findings on the the altered metabolic reprogramming of MM and the characteristics of MM metabolic-microbial interactions in order to gain an in-depth understanding of MM pathogenesis and progression and drug resistance mechanisms, and ultimately to explore for new strategies for MM treatment.

Chromosomal Instability in Genome Evolution: From Cancer to Macroevolution

The integrity of the genome is crucial for the survival of all living organisms. However, genomes need to adapt to survive certain pressures, and for this purpose use several mechanisms to diversify. Chromosomal instability (CIN) is one of the main mechanisms leading to the creation of genomic heterogeneity by altering the number of chromosomes and changing their structures. In this review, we will discuss the different chromosomal patterns and changes observed in speciation, in evolutional biology as well as during tumor progression. By nature, the human genome shows an induction of diversity during gametogenesis but as well during tumorigenesis that can conclude in drastic changes such as the whole genome doubling to more discrete changes as the complex chromosomal rearrangement chromothripsis. More importantly, changes observed during speciation are strikingly similar to the genomic evolution observed during tumor progression and resistance to therapy. The different origins of CIN will be treated as the importance of double-strand breaks (DSBs) or the consequences of micronuclei. We will also explain the mechanisms behind the controlled DSBs, and recombination of homologous chromosomes observed during meiosis, to explain how errors lead to similar patterns observed during tumorigenesis. Then, we will also list several diseases associated with CIN, resulting in fertility issues, miscarriage, rare genetic diseases, and cancer. Understanding better chromosomal instability as a whole is primordial for the understanding of mechanisms leading to tumor progression.

Antibody-Drug Conjugates for Multiple Myeloma: Just the Beginning, or the Beginning of the End?

Multiple myeloma is a malignancy of immunoglobulin-secreting plasma cells that is now often treated in the newly diagnosed and relapsed and/or refractory settings with monoclonal antibodies targeting lineage-specific markers used either alone or in rationally designed combination regimens. Among these are the anti-CD38 antibodies daratumumab and isatuximab, and the anti-Signaling lymphocytic activation molecule family member 7 antibody elotuzumab, all of which are used in their unconjugated formats. Single-chain variable fragments from antibodies also form a key element of the chimeric antigen receptors (CARs) in the B-cell maturation antigen (BCMA)-targeted CAR T-cell products idecabtagene vicleucel and ciltacabtagene autoleucel, which are approved in the advanced setting. Most recently, the bispecific anti-BCMA and T-cell-engaging antibody teclistamab has become available, again for patients with relapsed/refractory disease. Another format into which antibodies can be converted to exert anti-tumor efficacy is as antibody-drug conjugates (ADCs), and belantamab mafodotin, which also targets BCMA, represented the first such agent that gained a foothold in myeloma. Negative results from a recent Phase III study have prompted the initiation of a process for withdrawal of its marketing authorization. However, belantamab remains a drug with some promise, and many other ADCs targeting either BCMA or other plasma cell surface markers are in development and showing potential. This contribution will provide an overview of some of the current data supporting the possibility that ADCs will remain a part of our chemotherapeutic armamentarium against myeloma moving forward, and also highlight areas for future development.

High-dose melphalan treatment significantly increases mutational burden at relapse in multiple myeloma

High-dose melphalan (HDM) improves progression-free survival in multiple myeloma (MM), yet melphalan is a DNA-damaging alkylating agent; therefore, we assessed its mutational effect on surviving myeloma cells by analyzing paired MM samples collected at diagnosis and relapse in the IFM 2009 study. We performed deep whole-genome sequencing on samples from 68 patients, 43 of whom were treated with RVD (lenalidomide, bortezomib, and dexamethasone) and 25 with RVD + HDM. Although the number of mutations was similar at diagnosis in both groups (7137 vs 7230; P = .67), the HDM group had significantly more mutations at relapse (9242 vs 13 383, P = .005). No change in the frequency of copy number alterations or structural variants was observed. The newly acquired mutations were typically associated with DNA damage and double-stranded breaks and were predominantly on the transcribed strand. A machine learning model, using this unique pattern, predicted patients who would receive HDM with high sensitivity, specificity, and positive prediction value. Clonal evolution analysis showed that all patients treated with HDM had clonal selection, whereas a static progression was observed with RVD. A significantly higher percentage of mutations were subclonal in the HDM cohort. Intriguingly, patients treated with HDM who achieved complete remission (CR) had significantly more mutations at relapse yet had similar survival rates as those treated with RVD who achieved CR. This similarity could have been due to HDM relapse samples having significantly more neoantigens. Overall, our study identifies increased genomic changes associated with HDM and provides rationale to further understand clonal complexity.

Past, Present, and a Glance into the Future of Multiple Myeloma Treatment

Multiple myeloma (MM) is a challenging hematological cancer which typically grows in bone marrow. MM accounts for 10% of hematological malignancies and 1.8% of cancers. The recent treatment strategies have significantly improved progression-free survival for MM patients in the last decade; however, a relapse for most MM patients is inevitable. In this review we discuss current treatment, important pathways for proliferation, survival, immune suppression, and resistance that could be targeted for future treatments.

Amp(1q) and tetraploidy are commonly acquired chromosomal abnormalities in relapsed multiple myeloma

Long-term disease control in multiple myeloma (MM) is typically an unmet medical need, and most patients experience multiple relapses. Fluorescence in situ hybridization (FISH) is the standard technique to detect chromosomal abnormalities (CAs), which are important to estimate the prognosis of MM and the allocation of risk adapted therapies. In advanced stages, the importance of CAs needs further investigation. From 148 MM patients, two or more paired samples, at least one of which was collected at relapse, were analyzed by FISH. Using targeted next-generation sequencing, we molecularly investigated samples harboring relapse-associated CAs. Sixty-one percent of the patients showed a change in the cytogenetic profile during the disease course, including 10% who acquired high-risk cytogenetics. Amp(1q) (≥4 copies of 1q21), driven by an additional increase in copy number in patients who already had 3 copies of 1q21, was the most common acquired CA with 16% affected patients. Tetraploidy, found in 10% of the samples collected at the last time-point, was unstable over the course of the disease and was associated with TP53 lesions. Our results indicate that cytogenetic progression is common in relapsed patients. The relatively high frequency of amp(1q) suggests an active role for this CA in disease progression.

Therapy-selected clonal hematopoiesis and its role in myeloid neoplasms

Therapy-related myeloid neoplasms (t-MN) account for approximately 10-15% of all myeloid neoplasms and are associated with poor prognosis. Genomic characterization of t-MN to date has been limited in comparison to the considerable sequencing efforts performed for de novo myeloid neoplasms. Until recently, targeted deep sequencing (TDS) or whole exome sequencing (WES) have been the primary technologies utilized and thus limited the ability to explore the landscape of structural variants and mutational signatures. In the past decade, population-level studies have identified clonal hematopoiesis as a risk factor for the development of myeloid neoplasms. However, emerging research on clonal hematopoiesis as a risk factor for developing t-MN is evolving, and much is unknown about the progression of CH to t-MN. In this work, we will review the current knowledge of the genomic landscape of t-MN, discuss background knowledge of clonal hematopoiesis gained from studies of de novo myeloid neoplasms, and examine the recent literature studying the role of therapeutic selection of CH and its evolution under the effects of antineoplastic therapy. Finally, we will discuss the potential implications on current clinical practice and the areas of focus needed for future research into therapy-selected clonal hematopoiesis in myeloid neoplasms.

Molecular Diagnostics of Plasma Cell Neoplasms

Genetic characterization of myeloma at diagnosis by interphase fluorescence in situ hybridization and next-generation sequencing (NGS) can assist with risk stratification and treatment planning. Measurable residual disease (MRD) status after treatment, as evaluated by next-generation flow cytometry or NGS on bone marrow aspirate material, is one of the most important predictors of prognosis. Less-invasive tools for MRD assessment such as liquid biopsy approaches have also recently emerged as potential alternatives.

Whole-genome analysis identifies novel drivers and high-risk double-hit events in relapsed/refractory myeloma

Large-scale analyses of genomic data from patients with newly diagnosed multiple myeloma (ndMM) have been undertaken, however, large-scale analysis of relapsed/refractory MM (rrMM) has not been performed. We hypothesize that somatic variants chronicle the therapeutic exposures and clonal structure of myeloma from ndMM to rrMM stages. We generated whole-genome sequencing (WGS) data from 418 tumors (386 patients) derived from 6 rrMM clinical trials and compared them with WGS from 198 unrelated patients with ndMM in a population-based case-control fashion. We identified significantly enriched events at the rrMM stage, including drivers (DUOX2, EZH2, TP53), biallelic inactivation (TP53), noncoding mutations in bona fide drivers (TP53BP1, BLM), copy number aberrations (CNAs; 1qGain, 17pLOH), and double-hit events (Amp1q-ISS3, 1qGain-17p loss-of-heterozygosity). Mutational signature analysis identified a subclonal defective mismatch repair signature enriched in rrMM and highly active in high mutation burden tumors, a likely feature of therapy-associated expanding subclones. Further analysis focused on the association of genomic aberrations enriched at different stages of resistance to immunomodulatory agent (IMiD)-based therapy. This analysis revealed that TP53, DUOX2, 1qGain, and 17p loss-of-heterozygosity increased in prevalence from ndMM to lenalidomide resistant (LENR) to pomalidomide resistant (POMR) stages, whereas enrichment of MAML3 along with immunoglobulin lambda (IGL) and MYC translocations distinguished POM from the LEN subgroup. Genomic drivers associated with rrMM are those that confer clonal selective advantage under therapeutic pressure. Their role in therapy evasion should be further evaluated in longitudinal patient samples, to confirm these associations with the evolution of clinical resistance and to identify molecular subsets of rrMM for the development of targeted therapies.

Impact of rare structural variant events in newly diagnosed multiple myeloma

Whole genome sequencing (WGS) of newly diagnosed multiple myeloma patients (NDMM) has shown recurrent structural variant (SV) involvement in distinct regions of the genome (i.e. hotspots) and causing recurrent copy number alterations. Together with canonical immunoglobulin translocations, these SVs are recognized as "recurrent SVs". More than half SVs were not involved in recurrent events. The significance of these "rare SVs" has not been previously examined. In this study, we utilize 752 WGS and 591 RNA-seq data from NDMM patients to determine the role of rare SVs in myeloma pathogenesis. 94% of patients harbored at least one rare SV event. Rare SVs showed an SV-class specific enrichment within genes and superenhancers associated with outlier gene expression. Furthermore, known myeloma driver genes recurrently impacted by point mutations were dysregulated by rare SVs. Overall, we demonstrate the association of rare SVs with aberrant gene expression supporting a driver role in myeloma pathogenesis.

SIGNIFICANCE

Characterization of multiple myeloma genome revealed that more than half structural variants are not involved in recurrent events. Here, we demonstrate that these rare SVs hold potential for myeloma pathogenesis through their gene expression impact. Rare SVs contribute to MM heterogeneity and have implications for development of individualized treatment.

Integrative Analysis of Bulk RNA-Seq and Single-Cell RNA-Seq Unveils Novel Prognostic Biomarkers in Multiple Myeloma

Molecular heterogeneity has great significance in the disease biology of multiple myeloma (MM). Thus, the analysis combined single-cell RNA-seq (scRNA-seq) and bulk RNA-seq data were performed to investigate the clonal evolution characteristics and to find novel prognostic targets in MM. The scRNA-seq data were analyzed by the Seurat pipeline and Monocle 2 to identify MM cell branches with different differentiation states. Marker genes in each branch were uploaded to the STRING database to construct the Protein-Protein Interaction (PPI) network, followed by the detection of hub genes by Cytoscape software. Using bulk RNA-seq data, Kaplan-Meier (K-M) survival analysis was then carried out to determine prognostic biomarkers in MM. A total of 342 marker genes in two branches with different differentiation states were identified, and the top 20 marker genes with the highest scores in the network calculated by the MCC algorithm were selected as hub genes in MM. Furthermore, K-M survival analysis revealed that higher NDUFB8, COX6C, NDUFA6, USMG5, and COX5B expression correlated closely with a worse prognosis in MM patients. Moreover, ssGSEA and Pearson analyses showed that their expression had a significant negative correlation with the proportion of Tcm (central memory cell) immune cells. Our findings identified NDUFB8, COX6C, NDUFA6, USMG5, and COX5B as novel prognostic biomarkers in MM, and also revealed the significance of genetic heterogeneity during cell differentiation in MM prognosis.

Prevalence, mutational spectrum and clinical implications of clonal hematopoiesis of indeterminate potential in plasma cell dyscrasias

Clonal hematopoiesis of indeterminate potential (CHIP) is common both in healthy individuals and patients with hematological cancers. Recent studies have showed worse prognosis for patients with multiple myeloma (MM) and non-Hodgkin lymphoma undergoing stem cell transplant, that have concomitant presence of CHIP. Data regarding the clinical and biological role of CHIP in plasma cell dyscrasias (PCDs) is rapidly increasing. However, the prevalence and prognostic implication of CHIP in patients with MM outside of the transplant setting, and in those with other more indolent PCDs remains elusive. Here we explored the prevalence and clinical implications of CHIP detected through next-generation sequencing in 209 patients with PCDs including MM, light chain (AL) amyloidosis (ALA), monoclonal gammopathy of undetermined significance (MGUS), and smoldering multiple myeloma (SMM). To avoid attributing the mutations to the plasma cell clone, CHIP was defined as the presence of DNMT3A, TET2, or ASXL1 mutations in the peripheral blood or bone marrow (DTA-CH). The prevalence of DTA-CH was 19% in patients with PCDs, with no difference between each PCD. TET2 (23%) and DNMT3A (22%), were the most frequently mutated genes. DTA-CH correlated with older age in MM (P = .001) and MGUS/SMM (P = 0.0007), as well as with coronary artery disease or congestive heart failure in MM (P = .03). DTA-CH did not predict worse OS or PFS in either MM or ALA, nor it predict higher risk of progression to MM in patients with MGUS/SMM. Our results overall further elucidate the prevalence and mutational spectrum of CHIP in PCDs, providing more information regarding the clinical relevance of CHIP in this patient population.

A high-throughput electron tomography workflow reveals over-elongated centrioles in relapsed/refractory multiple myeloma

Electron microscopy is the gold standard to characterize centrosomal ultrastructure. However, production of significant morphometrical data is highly limited by acquisition time. We therefore developed a generalizable, semi-automated high-throughput electron tomography strategy to study centrosome aberrations in sparse patient-derived cancer cells at nanoscale. As proof of principle, we present electron tomography data on 455 centrioles of CD138pos plasma cells from one patient with relapsed/refractory multiple myeloma and CD138neg bone marrow mononuclear cells from three healthy donors as a control. Plasma cells from the myeloma patient displayed 122 over-elongated centrioles (48.8%). Particularly mother centrioles also harbored gross structural abnormalities, including fragmentation and disturbed microtubule cylinder formation, while control centrioles were phenotypically unremarkable. These data demonstrate the feasibility of our scalable high-throughput electron tomography strategy to study structural centrosome aberrations in primary tumor cells. Moreover, our electron tomography workflow and data provide a resource for the characterization of cell organelles beyond centrosomes.