Corrigendum: Protein arginine methyltransferase 1 is a therapeutic vulnerability in multiple myeloma

[This corrects the article DOI: 10.3389/fimmu.2023.1239614.].

Multiomic Mapping of Acquired Chromosome 1 Copy-Number and Structural Variants to Identify Therapeutic Vulnerabilities in Multiple Myeloma

PURPOSE

Chromosome 1 (chr1) copy-number abnormalities (CNA) and structural variants (SV) are frequent in newly diagnosed multiple myeloma (NDMM) and are associated with a heterogeneous impact on outcomes, the drivers of which are largely unknown.

EXPERIMENTAL DESIGN

A multiomic approach comprising CRISPR, gene mapping of CNAs and SVs, methylation, expression, and mutational analysis was used to document the extent of chr1 molecular variants and their impact on pathway utilization.

RESULTS

We identified two distinct groups of gain(1q): focal gains associated with limited gene-expression changes and a neutral prognosis, and whole-arm gains, which are associated with substantial gene-expression changes, complex genetics, and an adverse prognosis. CRISPR identified a number of dependencies on chr1 but only limited variants associated with acquired CNAs. We identified seven regions of deletion, nine of gain, three of chromothripsis (CT), and two of templated insertion (TI), which contain a number of potential drivers. An additional mechanism involving hypomethylation of genes at 1q may contribute to the aberrant gene expression of a number of genes. Expression changes associated with whole-arm gains were substantial and gene set enrichment analysis identified metabolic processes, apoptotic resistance, signaling via the MAPK pathway, and upregulation of transcription factors as being key drivers of the adverse prognosis associated with these variants.

CONCLUSIONS

Multiple layers of genetic complexity impact the phenotype associated with CNAs on chr1 to generate its associated clinical phenotype. Whole-arm gains of 1q are the critically important prognostic group that deregulate multiple pathways, which may offer therapeutic vulnerabilities.

A Sweet Voice: Acute Febrile Neutrophilic Dermatosis of the Larynx

SIGNIFICANCE STATEMENT

Acute febrile neutrophilic dermatosis (Sweet syndrome) is a rare idiopathic condition characterized by fever and whole-body rash of tender erythematous plaques of unknown etiology. Otorhinolaryngologic manifestations of the disease can be severe, yet they are sparsely reported in the literature. We present the first documented case of laryngeal involvement of Sweet syndrome.

Identifying novel mechanisms of biallelic TP53 loss refines poor outcome for patients with multiple myeloma

Biallelic TP53 inactivation is the most important high-risk factor associated with poor survival in multiple myeloma. Classical biallelic TP53 inactivation has been defined as simultaneous mutation and copy number loss in most studies; however, numerous studies have demonstrated that other factors could lead to the inactivation of TP53. Here, we hypothesized that novel biallelic TP53 inactivated samples existed in the multiple myeloma population. A random forest regression model that exploited an expression signature of 16 differentially expressed genes between classical biallelic TP53 and TP53 wild-type samples was subsequently established and used to identify novel biallelic TP53 samples from monoallelic TP53 groups. The model reflected high accuracy and robust performance in newly diagnosed relapsed and refractory populations. Patient survival of classical and novel biallelic TP53 samples was consistently much worse than those with mono-allelic or wild-type TP53 status. We also demonstrated that some predicted biallelic TP53 samples simultaneously had copy number loss and aberrant splicing, resulting in overexpression of high-risk transcript variants, leading to biallelic inactivation. We discovered that splice site mutation and overexpression of the splicing factor MED18 were reasons for aberrant splicing. Taken together, our study unveiled the complex transcriptome of TP53, some of which might benefit future studies targeting abnormal TP53.

Identifying 1q amplification and PHF19 expressing high-risk cells associated with relapsed/refractory multiple myeloma

Multiple Myeloma is an incurable plasma cell malignancy with a poor survival rate that is usually treated with immunomodulatory drugs (iMiDs) and proteosome inhibitors (PIs). The malignant plasma cells quickly become resistant to these agents causing relapse and uncontrolled growth of resistant clones. From whole genome sequencing (WGS) and RNA sequencing (RNA-seq) studies, different high-risk translocation, copy number, mutational, and transcriptional markers have been identified. One of these markers, PHF19, epigenetically regulates cell cycle and other processes and has already been studied using RNA-seq. In this study a massive (325,025 cells and 49 patients) single cell multiomic dataset was generated with jointly quantified ATAC- and RNA-seq for each cell and matched genomic profiles for each patient. We identified an association between one plasma cell subtype with myeloma progression that we have called relapsed/refractory plasma cells (RRPCs). These cells are associated with 1q alterations, TP53 mutations, and higher expression of PHF19. We also identified downstream regulation of cell cycle inhibitors in these cells, possible regulation of the transcription factor (TF) PBX1 on 1q, and determined that PHF19 may be acting primarily through this subset of cells.

Protein arginine methyltransferase 1 is a therapeutic vulnerability in multiple myeloma

Multiple myeloma (MM) is a devastating plasma cell malignancy characterized by the expansion of aberrant monoclonal plasma cells in the bone marrow, leading to severe clinical manifestations and poor prognosis, particularly in relapsed/refractory cases. Identifying novel therapeutic targets is crucial to improve treatment outcomes in these patients. In this study, we investigated the role of the protein arginine methyltransferase 1 (PRMT1) in MM pathogenesis and explored its potential as a therapeutic target. We observed that PRMT1, responsible for most asymmetric di-methylation in cells, exhibited the highest expression among PRMT family members in MM cell lines and primary MM cells. Importantly, PRMT1 expression was significantly elevated in relapsed/refractory patients compared to newly diagnosed patients. High expression of PRMT1 expression was strongly associated with poor prognosis. We found that genetic or enzymatic inhibition of PRMT1 impaired MM cell growth, induced cell cycle arrest, and triggered cell death. Treatment with MS023, a potent PRMT type I inhibitor, demonstrated a robust inhibitory effect on the viability of primary cells isolated from newly diagnosed and proteasome inhibitor-relapsed/refractory patients in a dose-dependent manner. Suppression of PRMT1 downregulated genes related to cell division and upregulated genes associated with apoptosis pathway. We also found that genes related to immune response and lymphocyte activation were significantly upregulated in PRMT1-suppressed cells. Notably, the activation status of T cells was strikingly enhanced upon co-culturing with PRMT1-KO MM cells. In vivo studies using a xenograft model revealed that targeting PRMT1 by either CRISPR/Cas9-mediated knockout or MS023 treatment significantly attenuated MM tumor growth and prolonged the survival of tumor-bearing mice. Histological analysis further confirmed increased apoptotic cell death in MS023-treated tumors. Collectively, our findings establish PRMT1 as an indispensable and novel therapeutic vulnerability in MM. The elevated expression of PRMT1 in relapsed/refractory patients underscores its potential as a target for overcoming treatment resistance. Moreover, our results highlight the efficacy of MS023 as a promising therapeutic agent against MM, offering new avenues for therapeutic approaches in relapsed/refractory MM.

Molecular characterization and clinical outcome of B-cell precursor acute lymphoblastic leukemia with IG-MYC rearrangement

Rarely, immunophenotypically immature B-cell precursor acute lymphoblastic leukemia (BCP-ALL) carries an immunoglobulin- MYC rearrangement (IG-MYC-r). This can result in diagnostic confusion with Burkitt lymphoma/leukemia and use of individualized treatment schedules of unproven efficacy. Here we compare the molecular characteristics of these conditions and investigate historic clinical outcome data. We identified 90 cases registered in a national BCP-ALL clinical trial/registry. When present, diagnostic material underwent cytogenetic, exome, methylome and transcriptome analyses. The outcomes analyzed were 3-year event-free survival and overall survival. IG-MYC-r was identified in diverse cytogenetic backgrounds, co-existing with either established BCP-ALL-specific abnormalities (high hyperdiploidy, n=3; KMT2A-rearrangement, n=6; iAMP21, n=1; BCR-ABL1, n=1); BCL2/BCL6-rearrangements (n=15); or, most commonly, as the only defining feature (n=64). Within this final group, precursor-like V(D)J breakpoints predominated (8/9) and KRAS mutations were common (5/11). DNA methylation identified a cluster of V(D)J-rearranged cases, clearly distinct from Burkitt leukemia/lymphoma. Children with IG-MYC-r within that subgroup had a 3-year event-free survival of 47% and overall survival of 60%, representing a high-risk BCP-ALL. To develop effective management strategies this group of patients must be allowed access to contemporary, minimal residual disease-adapted, prospective clinical trial protocols.

Alternative splicing in multiple myeloma is associated with the non-homologous end joining pathway

Alternative splicing plays a pivotal role in tumorigenesis and proliferation. However, its pattern and pathogenic role has not been systematically analyzed in multiple myeloma or its subtypes. Alternative splicing profiles for 598 newly diagnosed myeloma patients with comprehensive genomic annotation identified primary translocations, 1q amplification, and DIS3 events to have more differentially spliced events than those without. Splicing levels were correlated with expression of splicing factors. Moreover, the non-homologous end joining pathway was an independent factor that was highly associated with splicing frequency as well as an increased number of structural variants. We therefore identify an axis of high-risk disease encompassing expression of the non-homologous end joining pathway, increase structural variants, and increased alternative splicing that are linked together. This indicates a joint pathogenic role for DNA damage response and alternative RNA processing in myeloma.

The spatio-temporal evolution of multiple myeloma from baseline to relapse-refractory states

Deciphering Multiple Myeloma evolution in the whole bone marrow is key to inform curative strategies. Here, we perform spatial-longitudinal whole-exome sequencing, including 140 samples collected from 24 Multiple Myeloma patients during up to 14 years. Applying imaging-guided sampling we observe three evolutionary patterns, including relapse driven by a single-cell expansion, competing/co-existing sub-clones, and unique sub-clones at distinct locations. While we do not find the unique relapse sub-clone in the baseline focal lesion(s), we show a close phylogenetic relationship between baseline focal lesions and relapse disease, highlighting focal lesions as hotspots of tumor evolution. In patients with ≥3 focal lesions on positron-emission-tomography at diagnosis, relapse is driven by multiple distinct sub-clones, whereas in other patients, a single-cell expansion is typically seen (p < 0.01). Notably, we observe resistant sub-clones that can be hidden over years, suggesting that a prerequisite for curative therapies would be to overcome not only tumor heterogeneity but also dormancy.

High-resolution simulations of chromatin folding at genomic rearrangements in malignant B cells provide mechanistic insights into proto-oncogene deregulation

Genomic rearrangements are known to result in proto-oncogene deregulation in many cancers, but the link to 3D genome structure remains poorly understood. Here, we used the highly predictive heteromorphic polymer (HiP-HoP) model to predict chromatin conformations at the proto-oncogene CCND1 in healthy and malignant B cells. After confirming that the model gives good predictions of Hi-C data for the nonmalignant human B cell-derived cell line GM12878, we generated predictions for two cancer cell lines, U266 and Z-138. These possess genome rearrangements involving CCND1 and the immunoglobulin heavy locus (IGH), which we mapped using targeted genome sequencing. Our simulations showed that a rearrangement in U266 cells where a single IGH super-enhancer is inserted next to CCND1 leaves the local topologically associated domain (TAD) structure intact. We also observed extensive changes in enhancer-promoter interactions within the TAD, suggesting that it is the downstream chromatin remodeling which gives rise to the oncogene activation, rather than the presence of the inserted super-enhancer DNA sequence per se. Simulations of the IGH-CCND1 reciprocal translocation in Z-138 cells revealed that an oncogenic fusion TAD is created, encompassing CCND1 and the IGH super-enhancers. We predicted how the structure and expression of CCND1 changes in these different cell lines, validating this using qPCR and fluorescence in situ hybridization microscopy. Our work demonstrates the power of polymer simulations to predict differences in chromatin interactions and gene expression for different translocation breakpoints.

Perspectives on the Risk-Stratified Treatment of Multiple Myeloma

The multiple myeloma treatment landscape has changed dramatically. This change, paralleled by an increase in scientific knowledge, has resulted in significant improvement in survival. However, heterogeneity remains in clinical outcomes, with a proportion of patients not benefiting from current approaches and continuing to have a poor prognosis. A significant proportion of the variability in outcome can be predicted on the basis of clinical and biochemical parameters and tumor-acquired genetic variants, allowing for risk stratification and a more personalized approach to therapy. This article discusses the principles that can enable the rational and effective development of therapeutic approaches for high-risk multiple myeloma.

Myeloma Genome Project Panel is a Comprehensive Targeted Genomics Panel for Molecular Profiling of Patients with Multiple Myeloma

PURPOSE

We designed a comprehensive multiple myeloma targeted sequencing panel to identify common genomic abnormalities in a single assay and validated it against known standards.

EXPERIMENTAL DESIGN

The panel comprised 228 genes/exons for mutations, 6 regions for translocations, and 56 regions for copy number abnormalities (CNA). Toward panel validation, targeted sequencing was conducted on 233 patient samples and further validated using clinical FISH (translocations), multiplex ligation probe analysis (MLPA; CNAs), whole-genome sequencing (WGS; CNAs, mutations, translocations), or droplet digital PCR (ddPCR) of known standards (mutations).

RESULTS

Canonical immunoglobulin heavy chain translocations were detected in 43.2% of patients by sequencing, and aligned with FISH except for 1 patient. CNAs determined by sequencing and MLPA for 22 regions were comparable in 103 samples and concordance between platforms was R2 = 0.969. Variant allele frequency (VAF) for 74 mutations were compared between sequencing and ddPCR with concordance of R2 = 0.9849.

CONCLUSIONS

In summary, we have developed a targeted sequencing panel that is as robust or superior to FISH and WGS. This molecular panel is cost-effective, comprehensive, clinically actionable, and can be routinely deployed to assist risk stratification at diagnosis or posttreatment to guide sequencing of therapies.

Genetic subtypes of smoldering multiple myeloma are associated with distinct pathogenic phenotypes and clinical outcomes

Smoldering multiple myeloma (SMM) is a precursor condition of multiple myeloma (MM) with significant heterogeneity in disease progression. Existing clinical models of progression risk do not fully capture this heterogeneity. Here we integrate 42 genetic alterations from 214 SMM patients using unsupervised binary matrix factorization (BMF) clustering and identify six distinct genetic subtypes. These subtypes are differentially associated with established MM-related RNA signatures, oncogenic and immune transcriptional profiles, and evolving clinical biomarkers. Three genetic subtypes are associated with increased risk of progression to active MM in both the primary and validation cohorts, indicating they can be used to better predict high and low-risk patients within the currently used clinical risk stratification models.

Incidence of Concomitant Semicircular Canal Dehiscence With Otosclerosis

Objective

The concurrence of otosclerosis and superior semicircular canal dehiscence (SSCD) presents a diagnostic challenge and failure to differentiate between these 2 diagnoses results in mischaracterization and unsuccessful surgery. The objective of this study is to identify the incidence of SSCD in patients who have computed tomography (CT) evidence of otosclerosis.

Study Design

Retrospective chart review.

Setting

Tertiary referral hospital.

Patients

Adults with CT scan of the temporal bone diagnosed with radiological unilateral or bilateral fenestral otosclerosis from January 1995 to April 2018.

Methods

Retrospective review of patient imaging from a multi-center tertiary-referral health system from January 1995 to April 2018. Imaging was reviewed to quantify the incidence of SSCD among patients with CT-diagnosed bilateral fenestral otosclerosis. Poor quality imaging was excluded from review.

Results

One-thousand two-hundred eight patients (1214 CT scans) were identified with otosclerosis, of which 373 were diagnosed with fenestral otosclerosis (663 ears) with imaging of sufficient quality for review. This population was predominantly female (57.2%) with bilateral fenestral otosclerosis (78%). Of these, 23 ears (3.5%) had definitive evidence of SSCD, with an additional 15 ears (2.3%) with possible radiographic evidence of SSCD. There was no significant difference in laterality between the SSCD and otosclerosis.

Conclusions

Among 373 patients with fenestral otosclerosis per CT temporal bone imaging at a tertiary referral hospital, as many as 8.3% of patients had radiographic evidence of SSCD. Given this incidence, it continues to be important to consider SSCD when diagnosing and treating otosclerosis.

Structural variants shape the genomic landscape and clinical outcome of multiple myeloma

Deciphering genomic architecture is key to identifying novel disease drivers and understanding the mechanisms underlying myeloma initiation and progression. In this work, using the CoMMpass dataset, we show that structural variants (SV) occur in a nonrandom fashion throughout the genome with an increased frequency in the t(4;14), RB1, or TP53 mutated cases and reduced frequency in t(11;14) cases. By mapping sites of chromosomal rearrangements to topologically associated domains and identifying significantly upregulated genes by RNAseq we identify both predicted and novel putative driver genes. These data highlight the heterogeneity of transcriptional dysregulation occurring as a consequence of both the canonical and novel structural variants. Further, it shows that the complex rearrangements chromoplexy, chromothripsis and templated insertions are common in MM with each variant having its own distinct frequency and impact on clinical outcome. Chromothripsis is associated with a significant independent negative impact on clinical outcome in newly diagnosed cases consistent with its use alongside other clinical and genetic risk factors to identify prognosis.

Impact of Etiological Cytogenetic Abnormalities on the Depth of Immunoparesis and Survival in Newly Diagnosed Multiple Myeloma

INTRODUCTION/BACKGROUND

Immunoparesis, or low polyclonal immunoglobulin levels, is commonly seen in multiple myeloma (MM), and is associated with poor clinical outcomes. MM can be divided into subgroups with distinct biology and outcomes based on etiologic cytogenetic abnormalities. These include hyperdiploidy and translocations of t(11;14), t(4;14), t(14;16), and t(14;20), with the latter 3 associated with high-risk disease. We hypothesized that the different etiologic cytogenetic abnormalities drive bone marrow microenvironmental changes, resulting in different degrees of immunoparesis, and subgroup-dependent effects on clinical outcomes.

MATERIALS AND METHODS

We performed a retrospective review of 985 newly diagnosed patients enrolled in the Myeloma IX and XI trials. Immunoglobulin levels, survival outcomes, and infection rates were evaluated for each cytogenetic subgroup.

RESULTS

A significant proportion of patients with high-risk t(4;14), t(14;16), or t(14;20) had suppressed polyclonal immunoglobulins compared to standard-risk patients with hyperdiploidy or t(11;14). The clinical impact of immunoparesis depended on the cytogenetic subgroup, with the degree of IgM suppression effecting progression-free and overall survival only in the hyperdiploid subgroup. There was no significant difference in infection rates amongst the etiologic subgroups.

CONCLUSION

These findings demonstrate that the etiologic cytogenetic subgroup influences the degree and clinical impact of immunoparesis. This suggests that the underlying cytogenetic abnormality affects remodeling of the bone marrow plasma cell niche, resulting in suppressed normal plasma cell function, and low immunoglobulin levels.

Diagnostic Evidence GAuge of Single cells (DEGAS): a flexible deep transfer learning framework for prioritizing cells in relation to disease

We propose DEGAS (Diagnostic Evidence GAuge of Single cells), a novel deep transfer learning framework, to transfer disease information from patients to cells. We call such transferrable information "impressions," which allow individual cells to be associated with disease attributes like diagnosis, prognosis, and response to therapy. Using simulated data and ten diverse single-cell and patient bulk tissue transcriptomic datasets from glioblastoma multiforme (GBM), Alzheimer's disease (AD), and multiple myeloma (MM), we demonstrate the feasibility, flexibility, and broad applications of the DEGAS framework. DEGAS analysis on myeloma single-cell transcriptomics identified PHF19high myeloma cells associated with progression. Availability: https://github.com/tsteelejohnson91/DEGAS .

Epigenomic translocation of H3K4me3 broad domains over oncogenes following hijacking of super-enhancers

Chromosomal translocations are important drivers of haematological malignancies whereby proto-oncogenes are activated by juxtaposition with enhancers, often called enhancer hijacking. We analyzed the epigenomic consequences of rearrangements between the super-enhancers of the immunoglobulin heavy locus (IGH) and proto-oncogene CCND1 that are common in B cell malignancies. By integrating BLUEPRINT epigenomic data with DNA breakpoint detection, we characterized the normal chromatin landscape of the human IGH locus and its dynamics after pathological genomic rearrangement. We detected an H3K4me3 broad domain (BD) within the IGH locus of healthy B cells that was absent in samples with IGH-CCND1 translocations. The appearance of H3K4me3-BD over CCND1 in the latter was associated with overexpression and extensive chromatin accessibility of its gene body. We observed similar cancer-specific H3K4me3-BDs associated with hijacking of super-enhancers of other common oncogenes in B cell (MAF, MYC, and FGFR3/NSD2) and T cell malignancies (LMO2, TLX3, and TAL1). Our analysis suggests that H3K4me3-BDs can be created by super-enhancers and supports the new concept of epigenomic translocation, in which the relocation of H3K4me3-BDs from cell identity genes to oncogenes accompanies the translocation of super-enhancers.

Insights into high-risk multiple myeloma from an analysis of the role of PHF19 in cancer

Despite  improvements in outcome, 15-25% of newly diagnosed multiple myeloma (MM) patients have treatment resistant high-risk (HR) disease with a poor survival. The lack of a genetic basis for HR has focused attention on the role played by epigenetic changes. Aberrant expression and somatic mutations affecting genes involved in the regulation of tri-methylation of the lysine (K) 27 on histone 3 H3 (H3K27me3) are common in cancer. H3K27me3 is catalyzed by EZH2, the catalytic subunit of the Polycomb Repressive Complex 2 (PRC2). The deregulation of H3K27me3 has been shown to be involved in oncogenic transformation and tumor progression in a variety of hematological malignancies including MM. Recently we have shown that aberrant overexpression of the PRC2 subunit PHD Finger Protein 19 (PHF19) is the most significant overall contributor to HR status further focusing attention on the role played by epigenetic change in MM. By modulating both the PRC2/EZH2 catalytic activity and recruitment, PHF19 regulates the expression of key genes involved in cell growth and differentiation. Here we review the expression, regulation and function of PHF19 both in normal and the pathological contexts of solid cancers and MM. We present evidence that strongly implicates PHF19 in the regulation of genes important in cell cycle and the genetic stability of MM cells making it highly relevant to HR MM behavior. A detailed understanding of the normal and pathological functions of PHF19 will allow us to design therapeutic strategies able to target aggressive subsets of MM.

Chromothripsis as a pathogenic driver of multiple myeloma

Analysis of the genetic basis for multiple myeloma (MM) has informed many of our current concepts of the biology that underlies disease initiation and progression. Studying these events in further detail is predicted to deliver important insights into its pathogenesis, prognosis and treatment. Information from whole genome sequencing of structural variation is revealing the role of these events as drivers of MM. In particular, we discuss how the insights we have gained from studying chromothripsis suggest that it can be used to provide information on disease initiation and that, as a consequence, it can be used for the clinical classification of myeloma precursor diseases allowing for early intervention and prognostic determination. For newly diagnosed MM, the integration of information on the presence of chromothripsis has the potential to significantly enhance current risk prediction strategies and to better characterize patients with high-risk disease biology. In this article we summarize the genetic basis for MM and the role played by chromothripsis as a critical pathogenic factor active at early disease phases.

PHF19 inhibition as a therapeutic target in multiple myeloma

Epigenetic deregulation is increasingly recognized as a contributing pathological factor in multiple myeloma (MM). In particular tri-methylation of H3 lysine 27 (H3K27me3), which is catalyzed by PHD finger protein 19 (PHF19), a subunit of the Polycomb Repressive Complex 2 (PRC2), has recently shown to be a crucial mediator of MM tumorigenicity. Overexpression of PHF19 in MM has been associated with worse clinical outcome. Yet, while there is mounting evidence that PHF19 overexpression plays a crucial role in MM carcinogenesis downstream mechanisms remain to be elucidated. In the current study we use a functional knock down (KD) of PHF19 to investigate the biological role of PHF19 and show that PHF19KD leads to decreased tumor growth in vitro and in vivo. Expression of major cancer players such as bcl2, myc and EGR1 were decreased upon PHF19KD further underscoring the role of PHF19 in MM biology. Additionally, our results highlighted the prognostic impact of PHF19 overexpression, which was significantly associated with worse survival. Overall, our study underscores the premise that targeting the PHF19-PRC2 complex would open up avenues for novel MM therapies.

Cryopreservation Preserves Cell-Type Composition and Gene Expression Profiles in Bone Marrow Aspirates From Multiple Myeloma Patients

Single-cell RNA sequencing reveals gene expression differences between individual cells and also identifies different cell populations that are present in the bulk starting material. To obtain an accurate assessment of patient samples, single-cell suspensions need to be generated as soon as possible once the tissue or sample has been collected. However, this requirement poses logistical challenges for experimental designs involving multiple samples from the same subject since these samples would ideally be processed at the same time to minimize technical variation in data analysis. Although cryopreservation has been shown to largely preserve the transcriptome, it is unclear whether the freeze-thaw process might alter gene expression profiles in a cell-type specific manner or whether changes in cell-type proportions might also occur. To address these questions in the context of multiple myeloma clinical studies, we performed single-cell RNA sequencing (scRNA-seq) to compare fresh and frozen cells isolated from bone marrow aspirates of six multiple myeloma patients, analyzing both myeloma cells (CD138⁺) and cells constituting the microenvironment (CD138−). We found that cryopreservation using 90% fetal calf serum and 10% dimethyl sulfoxide resulted in highly consistent gene expression profiles when comparing fresh and frozen samples from the same patient for both CD138⁺ myeloma cells (R ≥ 0.96) and for CD138^– cells (R ≥ 0.9). We also demonstrate that CD138^– cell-type proportions showed minimal alterations, which were mainly related to small differences in immune cell subtype sensitivity to the freeze-thaw procedures. Therefore, when processing fresh multiple myeloma samples is not feasible, cryopreservation is a useful option in single-cell profiling studies.

Differential RNA splicing as a potentially important driver mechanism in multiple myeloma

Disruption of the normal splicing patterns of RNA is a major factor in the pathogenesis of a number of diseases. Increasingly research has shown the strong influence that splicing patterns can have on cancer progression. Multiple Myeloma is a molecularly heterogeneous disease classified by the presence of key translocations, gene expression profiles and mutations but the splicing patterns in MM remains largely unexplored. We take a multifaceted approach to define the extent and impact of alternative splicing in MM. We look at the spliceosome component, SF3B1, with hotspot mutations (K700E and K666T/Q) shown to result in an increase in alternative splicing in other cancers. We discovered a number of differentially spliced genes in comparison of the SF3B1 mutant and wild type samples that included, MZB1, DYNLL1, TMEM14C and splicing related genes DHX9, CLASRP, and SNRPE. We identified a broader role for abnormal splicing showing clear differences in the extent of novel splice variants in the different translocation groups. We show that a high number of novel splice loci is associated with adverse survival and an ultra-high risk group. The enumeration of patterns of alternative splicing has the potential to refine MM classification and to aid in the risk stratification of patients.

The molecular make up of smoldering myeloma highlights the evolutionary pathways leading to multiple myeloma

Smoldering myeloma (SMM) is associated with a high-risk of progression to myeloma (MM). We report the results of a study of 82 patients with both targeted sequencing that included a capture of the immunoglobulin and MYC regions. By comparing these results to newly diagnosed myeloma (MM) we show fewer NRAS and FAM46C mutations together with fewer adverse translocations, del(1p), del(14q), del(16q), and del(17p) in SMM consistent with their role as drivers of the transition to MM. KRAS mutations are associated with a shorter time to progression (HR 3.5 (1.5–8.1), p = 0.001). In an analysis of change in clonal structure over time we studied 53 samples from nine patients at multiple time points. Branching evolutionary patterns, novel mutations, biallelic hits in crucial tumour suppressor genes, and segmental copy number changes are key mechanisms underlying the transition to MM, which can precede progression and be used to guide early intervention strategies.

Progression from asymptomatic smoldering multiple myeloma (SMM) to symptomatic Multiple Myeloma occurs at different rates in different patients. Here, the authors report fewer NRAS and FAM46C mutations and adverse translocations in SMM compared to MM, while KRAS mutations are associated with a shorter time to progression.

Monitoring treatment response and disease progression in myeloma with circulating cell-free DNA

Circulating cell-free DNA (cfDNA) has the potential to capture spatial genetic heterogeneity in myeloma (MM) patients. We assessed whether cfDNA levels vary according to risk status defined by the 70 gene expression profile (GEP70). cfDNA levels in 77 patients were significantly higher in the GEP70 high-risk (HR) group compared to the low-risk (LR) group and correlated weakly with clinical markers including lactate dehydrogenase, β₂ -microglobulin, and ISS. Patients with high cfDNA levels were associated with a worse PFS (hazard ratio 6.4; 95% CI of ratio 1.9-22) and OS (hazard ratio 4.4; 95% CI of ratio 1.2-15.7). Circulating tumor DNA (ctDNA) was elevated in the HR group and ctDNA correlated strongly with GEP70 risk score (Spearman r = .69, P = .0027). cfDNA concentrations were significantly elevated between days 3-5 after chemotherapy before falling back to baseline levels. ctDNA in two patients showed a similar spike in levels between days 3 and 5 after chemotherapy with a concomitant increase in allele fraction of KRAS mutations. We assessed cfDNA levels in 25 patients with smoldering myeloma with serial samples and showed increased allele fraction of mutated KRAS at progression in cfDNA. Our study shows that cfDNA is a dynamic tool to capture genetic events in myeloma.

Bone marrow microenvironments that contribute to patient outcomes in newly diagnosed multiple myeloma: A cohort study of patients in the Total Therapy clinical trials

BACKGROUND

The tumor microenvironment (TME) is increasingly appreciated as an important determinant of cancer outcome, including in multiple myeloma (MM). However, most myeloma microenvironment studies have been based on bone marrow (BM) aspirates, which often do not fully reflect the cellular content of BM tissue itself. To address this limitation in myeloma research, we systematically characterized the whole bone marrow (WBM) microenvironment during premalignant, baseline, on treatment, and post-treatment phases.

METHODS AND FINDINGS

Between 2004 and 2019, 998 BM samples were taken from 436 patients with newly diagnosed MM (NDMM) at the University of Arkansas for Medical Sciences in Little Rock, Arkansas, United States of America. These patients were 61% male and 39% female, 89% White, 8% Black, and 3% other/refused, with a mean age of 58 years. Using WBM and matched cluster of differentiation (CD)138-selected tumor gene expression to control for tumor burden, we identified a subgroup of patients with an adverse TME associated with 17 fewer months of progression-free survival (PFS) (95% confidence interval [CI] 5-29, 49-69 versus 70-82 months, χ2 p = 0.001) and 15 fewer months of overall survival (OS; 95% CI -1 to 31, 92-120 versus 113-129 months, χ2 p = 0.036). Using immunohistochemistry-validated computational tools that identify distinct cell types from bulk gene expression, we showed that the adverse outcome was correlated with elevated CD8+ T cell and reduced granulocytic cell proportions. This microenvironment develops during the progression of premalignant to malignant disease and becomes less prevalent after therapy, in which it is associated with improved outcomes. In patients with quantified International Staging System (ISS) stage and 70-gene Prognostic Risk Score (GEP-70) scores, taking the microenvironment into consideration would have identified an additional 40 out of 290 patients (14%, premutation p = 0.001) with significantly worse outcomes (PFS, 95% CI 6-36, 49-73 versus 74-90 months) who were not identified by existing clinical (ISS stage III) and tumor (GEP-70) criteria as high risk. The main limitations of this study are that it relies on computationally identified cell types and that patients were treated with thalidomide rather than current therapies.

CONCLUSIONS

In this study, we observe that granulocyte signatures in the MM TME contribute to a more accurate prognosis. This implies that future researchers and clinicians treating patients should quantify TME components, in particular monocytes and granulocytes, which are often ignored in microenvironment studies.

The functional epigenetic landscape of aberrant gene expression in molecular subgroups of newly diagnosed multiple myeloma

Background

Multiple Myeloma (MM) is a hematological malignancy with genomic heterogeneity and poor survival outcome. Apart from the central role of genetic lesions, epigenetic anomalies have been identified as drivers in the development of the disease.

Methods

Alterations in the DNA methylome were mapped in 52 newly diagnosed MM (NDMM) patients of six molecular subgroups and matched with loci-specific chromatin marks to define their impact on gene expression. Differential DNA methylation analysis was performed using DMAP with a ≥10% increase (hypermethylation) or decrease (hypomethylation) in NDMM subgroups, compared to control samples, considered significant for all the subsequent analyses with p<0.05 after adjusting for a false discovery rate.

Results

We identified differentially methylated regions (DMRs) within the etiological cytogenetic subgroups of myeloma, compared to control plasma cells. Using gene expression data we identified genes that are dysregulated and correlate with DNA methylation levels, indicating a role for DNA methylation in their transcriptional control. We demonstrated that 70% of DMRs in the MM epigenome were hypomethylated and overlapped with repressive H3K27me3. In contrast, differentially expressed genes containing hypermethylated DMRs within the gene body or hypomethylated DMRs at the promoters overlapped with H3K4me1, H3K4me3, or H3K36me3 marks. Additionally, enrichment of BRD4 or MED1 at the H3K27ac enriched DMRs functioned as super-enhancers (SE), controlling the overexpression of genes or gene-cassettes.

Conclusions

Therefore, this study presents the underlying epigenetic regulatory networks of gene expression dysregulation in NDMM patients and identifies potential targets for future therapies.

KTP Laser Treatment of Early Glottic Cancer: A Multi-Institutional Retrospective Study

OBJECTIVES

The primary objectives were to report oncologic outcomes of transoral laser microsurgery with potassium-titanyl-phosphate (KTP) laser (TLM-KTP) ablation of early glottic cancer (EGC). The secondary objectives were to report vocal outcomes and to analyze factors that might influence outcomes.

METHODS

A multi-institutional, retrospective analysis of consecutive patients treated for T1 or T2 glottic squamous cell carcinoma undergoing TLM-KTP ablation with at least 2 years of follow-up was performed. Patients with prior radiation or surgery for laryngeal disease were excluded.

PRIMARY OUTCOME MEASURES INCLUDED

surgical failures requiring radiation or laryngectomy, disease-specific survival (DSS), and overall survival (OS). Secondary outcome measures included: pre- and postoperative Voice Handicap Index-10 (VHI-10) scores. The effects of smoking status, stage, and anterior commissure involvement on outcomes were analyzed.

RESULTS

Overall 88 patients met inclusion criteria (83% male, 79.5% current or former smokers). Mean age was 68 (standard deviation (SD): 12). Mean follow-up was 39.5 months (SD: 15.3). Staging included 50 T1a, 21 T1b, and 20 T2 tumors, including three metachronous second primaries. Radiation and/or laryngectomy avoidance was achieved in 87/88 (98.9%) of patients, inclusive of 24 patients requiring KTP re-treatments. Two patients had biopsy-proven recurrence (2.3%), but only 21 of 24 re-treated patients received a formal biopsy. No patients died from laryngeal cancer. DSS and OS were 100% and 92.3%, respectively. The mean VHI-10 scores were 19.3 preoperatively, 3.8 at 6-months postop, and 3.8 at 2-years postop. Smokers had a longer interval to re-treatment (P = .03), patients with T2 lesions had a shorter interval to re-treatment (0.02), and patients with T2 lesions presented with worse initial VHI-10 scores (0.002).

CONCLUSIONS

A multi-institutional, retrospective case series of TLM-KTP ablation of EGC demonstrated excellent oncologic outcomes when close surveillance and proactive re-treatments were utilized. Disease-specific survival, overall survival, and vocal function were excellent. Additional studies are necessary to further analyze the merits and risks of this treatment approach.

The Chromosome 13 Conundrum in Multiple Myeloma

In this issue of Blood Cancer Discovery, Chesi and colleagues have performed a series of mouse experiments, combined with patient sample analysis, to delineate the role of del(13) in multiple myeloma. They identify loss of the miRNA cluster MIR15A/16-1 as critical for myelomagenesis and progression of disease. See related article by Chesi et al., p. 68.

TarPan: an easily adaptable targeted sequencing panel viewer for research and clinical use

BACKGROUND

The study of cancer genomics continually matures as the number of patient samples sequenced increases. As more data is generated, oncogenic drivers for specific cancer types are discovered along with their associated risks. This in turn leads to potential treatment strategies that pave the way to precision medicine. However, significant financial and analytical barriers make it infeasible to sequence the entire genome of every patient. In contrast, targeted sequencing panels give reliable information on relevant portions of the genome at a fiscally responsible cost. Therefore, we have created the Targeted Panel (TarPan) Viewer, a software tool, to investigate this type of data.

RESULTS

TarPan Viewer helps investigators understand data from targeted sequencing data by displaying the information through a web browser interface. Through this interface, investigators can easily observe copy number changes, mutations, and structural events in cancer samples. The viewer runs in R Shiny with a robust SQLite backend and its input is generated from bioinformatic algorithms reliably described in the literature. Here we show the results from using TarPan Viewer on publicly available follicular lymphoma, breast cancer, and multiple myeloma data. In addition, we have tested and utilized the viewer internally, and this data has been used in high-impact peer-reviewed publications.

CONCLUSIONS

We have designed a flexible, simple to setup viewer that is easily adaptable to any type of cancer targeted sequencing, and has already proven its use in a research laboratory environment. Further, we believe with deeper sequencing and/or more targeted application it could be of use in the clinic in conjunction with an appropriate targeted sequencing panel as a cost-effective diagnostic test, especially in cancers such as acute leukemia or diffuse large B-cell lymphoma that require rapid interventions.

A tailored molecular profiling programme for children with cancer to identify clinically actionable genetic alterations

BACKGROUND

For children with cancer, the clinical integration of precision medicine to enable predictive biomarker-based therapeutic stratification is urgently needed.

METHODS

We have developed a hybrid-capture next-generation sequencing (NGS) panel, specifically designed to detect genetic alterations in paediatric solid tumours, which gives reliable results from as little as 50 ng of DNA extracted from formalin-fixed paraffin-embedded (FFPE) tissue. In this study, we offered an NGS panel, with clinical reporting via a molecular tumour board for children with solid tumours. Furthermore, for a cohort of 12 patients, we used a circulating tumour DNA (ctDNA)-specific panel to sequence ctDNA from matched plasma samples and compared plasma and tumour findings.

RESULTS

A total of 255 samples were submitted from 223 patients for the NGS panel. Using FFPE tissue, 82% of all submitted samples passed quality control for clinical reporting. At least one genetic alteration was detected in 70% of sequenced samples. The overall detection rate of clinically actionable alterations, defined by modified OncoKB criteria, for all sequenced samples was 51%. A total of 8 patients were sequenced at different stages of treatment. In 6 of these, there were differences in the genetic alterations detected between time points. Sequencing of matched ctDNA in a cohort of extracranial paediatric solid tumours also identified a high detection rate of somatic alterations in plasma.

CONCLUSION

We demonstrate that tailored clinical molecular profiling of both tumour DNA and plasma-derived ctDNA is feasible for children with solid tumours. Furthermore, we show that a targeted NGS panel-based approach can identify actionable genetic alterations in a high proportion of patients.

Transcriptome-wide association study of multiple myeloma identifies candidate susceptibility genes

BACKGROUND

While genome-wide association studies (GWAS) of multiple myeloma (MM) have identified variants at 23 regions influencing risk, the genes underlying these associations are largely unknown. To identify candidate causal genes at these regions and search for novel risk regions, we performed a multi-tissue transcriptome-wide association study (TWAS).

RESULTS

GWAS data on 7319 MM cases and 234,385 controls was integrated with Genotype-Tissue Expression Project (GTEx) data assayed in 48 tissues (sample sizes, N = 80-491), including lymphocyte cell lines and whole blood, to predict gene expression. We identified 108 genes at 13 independent regions associated with MM risk, all of which were in 1 Mb of known MM GWAS risk variants. Of these, 94 genes, located in eight regions, had not previously been considered as a candidate gene for that locus.

CONCLUSIONS

Our findings highlight the value of leveraging expression data from multiple tissues to identify candidate genes responsible for GWAS associations which provide insight into MM tumorigenesis. Among the genes identified, a number have plausible roles in MM biology, notably APOBEC3C, APOBEC3H, APOBEC3D, APOBEC3F, APOBEC3G, or have been previously implicated in other malignancies. The genes identified in this TWAS can be explored for follow-up and validation to further understand their role in MM biology.

An acquired high-risk chromosome instability phenotype in multiple myeloma: Jumping 1q Syndrome

Patients with multiple myeloma (MM) accumulate adverse copy number aberrations (CNAs), gains of 1q21, and 17p deletions during disease progression. A subset of these patients develops heightened 1q12 pericentromeric instability and jumping translocations of 1q12 (JT1q12), evidenced by increased copy CNAs of 1q21 and losses in receptor chromosomes (RC). To understand the progression of these aberrations we analyzed metaphase cells of 50 patients with ≥4 CNAs of 1q21 by G-banding, locus specific FISH, and spectral karyotyping. In eight patients with ≥5 CNAs of 1q21 we identified a chromosome instability phenotype similar to that found in ICF syndrome (immunodeficiency, centromeric instability, and facial anomalies). Strikingly, the acquired instability phenotype identified in these patients demonstrates the same transient structural aberrations of 1q12 as those found in ICF syndrome, suggesting similar underlying pathological mechanisms. Four types of clonal aberrations characterize this phenotype including JT1q12s, RC deletions, 1q12-21 breakage-fusion-bridge cycle amplifications, and RC insertions. In addition, recurring transient aberrations include 1q12 decondensation and breakage, triradials, and 1q micronuclei. The acquired self-propagating mobile property of 1q12 satellite DNA drives the continuous regeneration of 1q12 duplication/deletion events. For patients demonstrating this instability phenotype, we propose the term "Jumping 1q Syndrome."

Clonal evolution in myeloma: the impact of maintenance lenalidomide and depth of response on the genetics and sub-clonal structure of relapsed disease in uniformly treated newly diagnosed patients

The emergence of treatment resistant sub-clones is a key feature of relapse in multiple myeloma. Therapeutic attempts to extend remission and prevent relapse include maximizing response and the use of maintenance therapy. We used whole exome sequencing to study the genetics of paired samples taken at presentation and at relapse from 56 newly diagnosed patients, following induction therapy, randomized to receive either lenalidomide maintenance or observation as part of the Myeloma XI trial. Patients included were considered high risk, relapsing within 30 months of maintenance randomization. Patients achieving a complete response had predominantly branching evolutionary patterns leading to relapse, characterized by a greater mutational burden, an altered mutational profile, bi-allelic inactivation of tumor suppressor genes, and acquired structural aberrations. Conversely, in patients achieving a partial response, the evolutionary features were predominantly stable with a similar mutational and structural profile seen at both time points. There were no significant differences between patients relapsing after lenalidomide maintenance versus observation. This study shows that the depth of response is a key determinant of the evolutionary patterns seen at relapse. This trial is registered at clinicaltrials.gov identifier: 01554852.

Microhomology-mediated end joining drives complex rearrangements and overexpression of MYC and PVT1 in multiple myeloma

MYC is a widely acting transcription factor and its deregulation is a crucial event in many human cancers. MYC is important biologically and clinically in multiple myeloma, but the mechanisms underlying its dysregulation are poorly understood. We show that MYC rearrangements are present in 36.0% of newly diagnosed myeloma patients, as detected in the largest set of next generation sequencing data to date (n=1,267). Rearrangements were complex and associated with increased expression of MYC and PVT1, but not other genes at 8q24. The highest effect on gene expression was detected in cases where the MYC locus is juxtaposed next to super-enhancers associated with genes such as IGH, IGK, IGL, TXNDC5/BMP6, FAM46C and FOXO3. We identified three hotspots of recombination at 8q24, one of which is enriched for IGH-MYC translocations. Breakpoint analysis indicates primary myeloma rearrangements involving the IGH locus occur through non-homologous end joining, whereas secondary MYC rearrangements occur through microhomology-mediated end joining. This mechanism is different to lymphomas, where non-homologous end joining generates MYC rearrangements. Rearrangements resulted in overexpression of key genes and chromatin immunoprecipitation-sequencing identified that HK2, a member of the glucose metabolism pathway, is directly over-expressed through binding of MYC at its promoter.

Daratumumab and dexamethasone is safe and effective for triple refractory myeloma patients: final results of the IFM 2014-04 (Etoile du Nord) trial

Single agent daratumumab has shown clinical activity in relapsed, refractory multiple myeloma (RRMM). The Intergroupe Francophone du Myélome 2014-04 trial was designed to further investigate daratumumab in combination with dexamethasone in triple RRMM patients. Patients received daratumumab infusions in combination with weekly dexamethasone until disease progression or unacceptable toxicity. Fifty-seven patients were included in the trial and evaluable for response. The overall response rate and the clinical benefit rate were 33% (n = 19) and 48% (n = 27), respectively. Five (8·8%) patients achieved a very good partial response or better. The median time to response was 4 weeks. For responding patients, the median progression-free survival was 6·6 months, compared to 3·7 months (3·0-5·5) for those with a minimal or stable disease. The median overall survival (OS) for all patients was 16·7 months (11·2-24·0). For responding patients, the median OS was 23·23 months, whereas that of patients with progressive disease was 2·97 months. The incidence of infusion-related reactions was 37%; all cases were manageable and did not lead to dose reduction or permanent treatment discontinuation. These data demonstrate that treatment with daratumumab and dexamethasone results in a meaningful long-term benefit with an acceptable safety profile for patients with triple RRMM.

Amplitude modulation detection and temporal modulation cutoff frequency in normal hearing infants

The goal of this study was to determine if temporal modulation cutoff frequency was mature in three-month-old infants. Normal-hearing infants and young adults were tested in a single-interval forced-choice observer-based psychoacoustic procedure. Two parameters of the temporal modulation transfer function (TMTF) were estimated to separate temporal resolution from amplitude modulation sensitivity. The modulation detection threshold (MDT) of a broadband noise amplitude modulated at 10 Hz estimated the y-intercept of the TMTF. The cutoff frequency of the TMTF, measured at a modulation depth 4 dB greater than the MDT, provided an estimate of temporal resolution. MDT was obtained in 27 of 33 infants while both MDT and cutoff frequency was obtained in 15 infants and in 16 of 16 adults. Mean MDT was approximately 10 dB poorer in infants compared to adults. In contrast, mean temporal modulation cutoff frequency did not differ significantly between age groups. These results suggest that temporal resolution is mature, on average, by three months of age in normal hearing children despite immature sensitivity to amplitude modulation. The temporal modulation cutoff frequency approach used here may be a feasible way to examine development of temporal resolution in young listeners with markedly immature sensitivity to amplitude modulation.

Poor overall survival in hyperhaploid multiple myeloma is defined by double-hit bi-allelic inactivation of TP53

Hyperhaploid multiple myeloma is a rare numerical aberration group defined by a range of 24-34 chromosomes, which is associated with a poor prognosis with a 5-year survival rate of 23%. Hyperhaploid patient samples (n=8) were sequenced and copy number and mutations identified. Samples had a median of 13 monosomies (range 12-14), which in general were those not associated with trisomies in hyperdiploid samples. The chromosomes traditionally trisomic in hyperdiploid myeloma were disomic in hyperhaploid myeloma with retention of heterodisomy. We examined the hyperhaploid samples for frequently mutated genes and found that 8/8 (100%) hyperhaploid samples had a mutation in TP53, exceeding the overall rate of mutation in newly diagnosed patients (5.5%), indicating an oncogenic dependency in this group. All samples with TP53 mutation also had monosomy of chromosome 17, indicating bi-allelic inactivation of TP53. As such, this high risk group is part of double-hit myeloma.

Author Correction: Identification of multiple risk loci and regulatory mechanisms influencing susceptibility to multiple myeloma

The original version of this Article contained an error in the spelling of a member of the PRACTICAL Consortium, Manuela Gago-Dominguez, which was incorrectly given as Manuela Gago Dominguez. This has now been corrected in both the PDF and HTML versions of the Article. Furthermore, in the original HTML version of this Article, the order of authors within the author list was incorrect. The PRACTICAL consortium was incorrectly listed after Richard S. Houlston and should have been listed after Nora Pashayan. This error has been corrected in the HTML version of the Article; the PDF version was correct at the time of publication.

Genetic correlation between multiple myeloma and chronic lymphocytic leukaemia provides evidence for shared aetiology

The clustering of different types of B-cell malignancies in families raises the possibility of shared aetiology. To examine this, we performed cross-trait linkage disequilibrium (LD)-score regression of multiple myeloma (MM) and chronic lymphocytic leukaemia (CLL) genome-wide association study (GWAS) data sets, totalling 11,734 cases and 29,468 controls. A significant genetic correlation between these two B-cell malignancies was shown (Rg = 0.4, P = 0.0046). Furthermore, four of the 45 known CLL risk loci were shown to associate with MM risk and five of the 23 known MM risk loci associate with CLL risk. By integrating eQTL, Hi-C and ChIP-seq data, we show that these pleiotropic risk loci are enriched for B-cell regulatory elements and implicate B-cell developmental genes. These data identify shared biological pathways influencing the development of CLL and, MM and further our understanding of the aetiological basis of these B-cell malignancies.

Combination of flow cytometry and functional imaging for monitoring of residual disease in myeloma

The iliac crest is the sampling site for minimal residual disease (MRD) monitoring in multiple myeloma (MM). However, the disease distribution is often heterogeneous, and imaging can be used to complement MRD detection at a single site. We have investigated patients in complete remission (CR) during first-line or salvage therapy for whom MRD flow cytometry and the two imaging modalities positron emission tomography (PET) and diffusion-weighted magnetic resonance imaging (DW-MRI) were performed at the onset of CR. Residual focal lesions (FLs), detectable in 24% of first-line patients, were associated with short progression-free survival (PFS), with DW-MRI detecting disease in more patients. In some patients, FLs were only PET positive, indicating that the two approaches are complementary. Combining MRD and imaging improved prediction of outcome, with double-negative and double-positive features defining groups with excellent and dismal PFS, respectively. FLs were a rare event (12%) in first-line MRD-negative CR patients. In contrast, patients achieving an MRD-negative CR during salvage therapy frequently had FLs (50%). Multi-region sequencing and imaging in an MRD-negative patient showed persistence of spatially separated clones. In conclusion, we show that DW-MRI is a promising tool for monitoring residual disease that complements PET and should be combined with MRD.

Identification of multiple risk loci and regulatory mechanisms influencing susceptibility to multiple myeloma

Genome-wide association studies (GWAS) have transformed our understanding of susceptibility to multiple myeloma (MM), but much of the heritability remains unexplained. We report a new GWAS, a meta-analysis with previous GWAS and a replication series, totalling 9974 MM cases and 247,556 controls of European ancestry. Collectively, these data provide evidence for six new MM risk loci, bringing the total number to 23. Integration of information from gene expression, epigenetic profiling and in situ Hi-C data for the 23 risk loci implicate disruption of developmental transcriptional regulators as a basis of MM susceptibility, compatible with altered B-cell differentiation as a key mechanism. Dysregulation of autophagy/apoptosis and cell cycle signalling feature as recurrently perturbed pathways. Our findings provide further insight into the biological basis of MM.