Analysis of mutational signatures is becoming routine in cancer genomics, with implications for pathogenesis, classification, prognosis, and even treatment decisions. However, the field lacks a consensus on analysis and result interpretation. Using whole-genome sequencing of multiple myeloma (MM), chronic lymphocytic leukemia (CLL) and acute myeloid leukemia, we compare the performance of public signature analysis tools. We describe caveats and pitfalls of de novo signature extraction and fitting approaches, reporting on common inaccuracies: erroneous signature assignment, identification of localized hyper-mutational processes, overcalling of signatures. We provide reproducible solutions to solve these issues and use orthogonal approaches to validate our results. We show how a comprehensive mutational signature analysis may provide relevant biological insights, reporting evidence of c-AID activity among unmutated CLL cases or the absence of BRCA1/BRCA2-mediated homologous recombination deficiency in a MM cohort. Finally, we propose a general analysis framework to ensure production of accurate and reproducible mutational signature data.

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 2 nd -3 rd 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 premalignant 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.

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.

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 2 nd-3 rd 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 premalignant 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.

Chromothripsis is detectable in 20–30% of newly diagnosed multiple myeloma (NDMM) patients and is emerging as a new independent adverse prognostic factor. In this study we interrogate 752 NDMM patients using whole genome sequencing (WGS) to investigate the relationship of copy number (CN) signatures to chromothripsis and show they are highly associated. CN signatures are highly predictive of the presence of chromothripsis (AUC = 0.90) and can be used identify its adverse prognostic impact. The ability of CN signatures to predict the presence of chromothripsis is confirmed in a validation series of WGS comprised of 235 hematological cancers (AUC = 0.97) and an independent series of 34 NDMM (AUC = 0.87). We show that CN signatures can also be derived from whole exome data (WES) and using 677 cases from the same series of NDMM, we are able to predict both the presence of chromothripsis (AUC = 0.82) and its adverse prognostic impact. CN signatures constitute a flexible tool to identify the presence of chromothripsis and is applicable to WES and WGS data.

Writing-review and editing: all the authors All the authors approved the final version of the manuscript and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Supplemental Digital Content is available in the text

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.