Next-Generation Sequencing for Clinical Management of Multiple Myeloma : Ready for Prime Time?

Personalized treatment is an attractive strategy that promises increased efficacy with reduced side effects in cancer. The feasibility of such an approach has been greatly boosted by next-generation sequencing (NGS) techniques, which can return detailed information on the genome and on the transcriptome of each patient's tumor, thus highlighting biomarkers of response or druggable targets that may differ from case to case. However, while the number of cancers sequenced is growing exponentially, much fewer cases are amenable to a molecularly-guided treatment outside of clinical trials to date. In multiple myeloma, genomic analysis shows a variety of gene mutations, aneuploidies, segmental copy-number changes, translocations that are extremely heterogeneous, and more numerous than other hematological malignancies. Currently, in routine clinical practice we employ reduced FISH panels that only capture three high-risk features as part of the R-ISS. On the contrary, recent advances have suggested that extending genomic analysis to the full spectrum of recurrent mutations and structural abnormalities in multiple myeloma may have biological and clinical implications. Furthermore, increased efficacy of novel treatments can now produce deeper responses, and standard methods do not have enough sensitivity to stratify patients in complete biochemical remission. Consequently, NGS techniques have been developed to monitor the size of the clone to a sensitivity of up to a cell in a million after treatment. However, even these techniques are not within reach of standard laboratories. In this review we will recapitulate recent advances in multiple myeloma genomics, with special focus on the ones that may have immediate translational impact. We will analyze the benefits and pitfalls of NGS-based diagnostics, highlighting crucial aspects that will need to be taken into account before this can be implemented in most laboratories. We will make the point that a new era in myeloma diagnostics and minimal residual disease monitoring is close and conventional genetic testing will not be able to return the required information. This will mandate that even in routine practice NGS should soon be adopted owing to a higher informative potential with increasing clinical benefits

Application of next-generation sequencing for the genomic characterization of patients with smoldering myeloma

Genomic analysis could contribute to a better understanding of the biological determinants of the evolution of multiple myeloma (MM) precursor disease and an improved definition of high-risk patients. To assess the feasibility and value of next-generation sequencing approaches in an asymptomatic setting, we performed a targeted gene mutation analysis and a genome-wide assessment of copy number alterations (CNAs) by ultra-low-pass whole genome sequencing (ULP-WGS) in six patients with monoclonal gammopathy of undetermined significance and 25 patients with smoldering MM (SMM). Our comprehensive genomic characterization highlighted heterogeneous but substantial values of the tumor fraction, especially in SMM; a rather high degree of genomic complexity, in terms of both mutations and CNAs, and inter-patient variability; a higher incidence of gene mutations and CNAs in SMM, confirming ongoing evolution; intraclonal heterogeneity; and instances of convergent evolution. ULP-WGS of these patients proved effective in revealing the marked genome-wide level of their CNAs, most of which are not routinely investigated. Finally, the analysis of our small SMM cohort suggested that chr(8p) deletions, the DNA tumor fraction, and the number of alterations may have clinical relevance in the progression to overt MM. Although validation in larger series is mandatory, these findings highlight the promising impact of genomic approaches in the clinical management of SMM

Integrative analysis of the genomic and transcriptomic landscape of double-refractory multiple myeloma

In multiple myeloma, novel treatments with proteasome inhibitors (PIs) and immunomodulatory agents (IMiDs) have prolonged survival but the disease remains incurable. At relapse, next-generation sequencing has shown occasional mutations of drug targets but has failed to identify unifying features that underlie chemotherapy resistance. We studied 42 patients refractory to both PIs and IMiDs. Whole-exome sequencing was performed in 40 patients, and RNA sequencing (RNA-seq) was performed in 27. We found more mutations than were reported at diagnosis and more subclonal mutations, which implies ongoing evolution of the genome of myeloma cells during treatment. The mutational landscape was different from that described in published studies on samples taken at diagnosis. The TP53 pathway was the most frequently inactivated (in 45% of patients). Conversely, point mutations of genes associated with resistance to IMiDs were rare and were always subclonal. Refractory patients were uniquely characterized by having a mutational signature linked to exposure to alkylating agents, whose role in chemotherapy resistance and disease progression remains to be elucidated. RNA-seq analysis showed that treatment or mutations had no influence on clustering, which was instead influenced by karyotypic events. We describe a cluster with both amp(1q) and del(13) characterized by CCND2 upregulation and also overexpression of MCL1, which represents a novel target for experimental treatments. Overall, high-risk features were found in 65% of patients. However, only amp(1q) predicted survival. Gene mutations of IMiD and PI targets are not a preferred mode of drug resistance in myeloma. Chemotherapy resistance of the bulk tumor population is likely attained through differential, yet converging evolution of subclones that are overall variable from patient to patient and within the same patient

Timing the initiation of multiple myeloma

The evolution and progression of multiple myeloma and its precursors over time is poorly understood. Here, we investigate the landscape and timing of mutational processes shaping multiple myeloma evolution in a large cohort of 89 whole genomes and 973 exomes. We identify eight processes, including a mutational signature caused by exposure to melphalan. Reconstructing the chronological activity of each mutational signature, we estimate that the initial transformation of a germinal center B-cell usually occurred during the first 2nd-3rd decades of life. We define four main patterns of activation-induced deaminase (AID) and apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) mutagenesis over time, including a subset of patients with evidence of prolonged AID activity during the pre-malignant phase, indicating antigen-responsiveness and germinal center reentry. Our findings provide a framework to study the etiology of multiple myeloma and explore strategies for prevention and early detection

mmsig: a fitting approach to accurately identify somatic mutational signatures in hematological malignancies

Mutational signatures have emerged as powerful biomarkers in cancer patients, with prognostic and therapeutic implications. Wider clinical utility requires access to reproducible algorithms, which allow characterization of mutational signatures in a given tumor sample. Here, we show how mutational signature fitting can be applied to hematological cancer genomes to identify biologically and clinically important mutational processes, highlighting the importance of careful interpretation in light of biological knowledge. Our newly released R package mmsig comes with a dynamic error-suppression procedure that improves specificity in important clinical and biological settings. In particular, mmsig allows accurate detection of mutational signatures with low abundance, such as those introduced by APOBEC cytidine deaminases. This is particularly important in the most recent mutational signature reference (COSMIC v3.1) where each signature is more clearly defined. Our mutational signature fitting algorithm mmsig is a robust tool that can be implemented immediately in the clinic

Waldenström macroglobulinemia whole genome reveals prolonged germinal center activity and late copy number aberrations

The genomic landscape of Waldenström macroglobulinemia (WM) is characterized by somatic mutations in MYD88, present from the precursor stages. Using the comprehensive resolution of whole genome sequencing (WGS) in 14 CD19-selected primary WM samples; comparing clonal and subclonal mutations revealed that germinal center (GC) mutational signatures SBS9 (poly-eta) and SBS84 (AID) have sustained activity, suggesting that the interaction betweenWMand the GC continues over time. Expanding our cohort size with 33 targeted sequencing samples, we interrogated the WM copy number aberration (CNA) landscape and chronology. Of interest, CNA prevalence progressively increased in symptomatic WM and relapsed disease when compared with stable precursor stages, with stable precursors lacking genomic complexity. Two MYD88 wild-type WGS contained a clonal gain affecting chromosome 12, which is typically an early event in chronic lymphocytic leukemia. Molecular time analysis demonstrated that both chromosomal 12 gain events occurred early in cancer development whereas other CNA changes tend to occur later in the disease course and are often subclonal. In summary, WGS analysis in WM allows the demonstration of sustained GC activity over time and allows the reconstruction of the temporal evolution of specific genomic features. In addition, our data suggest that, although MYD88-mutations are central to WM clone establishment and can be observed in precursor disease, CNA may contribute to later phases, and may be used as a biomarker for progression risk from precursor conditions to symptomatic disease. © 2023 by The American Society of Hematology

Early Relapse Risk in Newly Diagnosed Multiple Myeloma Patients Characterized by Next-Generation Sequencing

Purpose: Duration of first remission is important for the survival of patients with multiple myeloma. Experimental design: From the CoMMpass study (NCT01454297), 926 patients with newly diagnosed multiple myeloma, characterized by next-generation sequencing, were analyzed to evaluate those who experienced early progressive disease (PD; time to progression, TTP 6418 months). Results: After a median follow-up of 39 months, early PD was detected in 191/926 (20.6%) patients, 228/926 (24.6%) patients had late PD (TTP >18 months), while 507/926 (54.8%) did not have PD at the current follow-up. Compared with patients with late PD, patients with early PD had a lower at least very good partial response rate (47% vs. 82%, P < 0.001) and more frequently acquired double refractoriness to immunomodulatory drugs (IMiD) and proteasome inhibitors (PI; 21% vs. 8%, P < 0.001). Patients with early PD were at higher risk of death compared with patients with late PD and no PD (HR, 3.65; 95% CI, 2.7-4.93; P < 0.001), showing a dismal median overall survival (32.8 months). In a multivariate logistic regression model, independent factors increasing the early PD risk were TP53 mutation (OR, 3.78, P < 0.001), high lactate dehydrogenase levels (OR, 3.15, P = 0.006), \u3bb-chain translocation (OR, 2.25, P = 0.033), and IGLL5 mutation (OR, 2.15, P = 0.007). Carfilzomib-based induction (OR, 0.15, P = 0.014), autologous stem-cell transplantation (OR, 0.27, P < 0.001), and continuous therapy with PIs and IMiDs (OR, 0.34, P = 0.024) mitigated the risk of early PD. Conclusions: Early PD identifies a high-risk multiple myeloma population. Further research is needed to better identify baseline features predicting early PD and the optimal treatment approaches for patients at risk