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.

Targeting DNA2 overcomes metabolic reprogramming in multiple myeloma.

DNA damage resistance is a major barrier to effective DNA-damaging therapy in multiple myeloma (MM). To discover mechanisms through which MM cells overcome DNA damage, we investigate how MM cells become resistant to antisense oligonucleotide (ASO) therapy targeting Interleukin enhancer binding factor 2 (ILF2), a DNA damage regulator that is overexpressed in 70% of MM patients whose disease has progressed after standard therapies have failed. Here, we show that MM cells undergo adaptive metabolic rewiring to restore energy balance and promote survival in response to DNA damage activation. Using a CRISPR/Cas9 screening strategy, we identify the mitochondrial DNA repair protein DNA2, whose loss of function suppresses MM cells' ability to overcome ILF2 ASO-induced DNA damage, as being essential to counteracting oxidative DNA damage. Our study reveals a mechanism of vulnerability of MM cells that have an increased demand for mitochondrial metabolism upon DNA damage activation.

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.

Transcriptional Heterogeneity Overcomes Super-Enhancer Disrupting Drug Combinations in Multiple Myeloma.

Multiple myeloma (MM) is a malignancy that is often driven by MYC and that is sustained by IRF4, which are upregulated by super-enhancers. IKZF1 and IKZF3 bind to super-enhancers and can be degraded using immunomodulatory imide drugs (IMiD). Successful IMiD responses downregulate MYC and IRF4; however, this fails in IMiD-resistant cells. MYC and IRF4 downregulation can also be achieved in IMiD-resistant tumors using inhibitors of BET and EP300 transcriptional coactivator proteins; however, in vivo these drugs have a narrow therapeutic window. By combining IMiDs with EP300 inhibition, we demonstrate greater downregulation of MYC and IRF4, synergistic killing of myeloma in vitro and in vivo, and an increased therapeutic window. Interestingly, this potent combination failed where MYC and IRF4 expression was maintained by high levels of the AP-1 factor BATF. Our results identify an effective drug combination and a previously unrecognized mechanism of IMiD resistance. SIGNIFICANCE: These results highlight the dependence of MM on IKZF1-bound super-enhancers, which can be effectively targeted by a potent therapeutic combination pairing IMiD-mediated degradation of IKZF1 and IKZF3 with EP300 inhibition. They also identify AP-1 factors as an unrecognized mechanism of IMiD resistance in MM. See related article by Neri, Barwick, et al., p. 56. See related commentary by Yun and Cleveland, p. 5. This article is featured in Selected Articles from This Issue, p. 4.

ETV4-Dependent Transcriptional Plasticity Maintains MYC Expression and Results in IMiD Resistance in Multiple Myeloma.

Immunomodulatory drugs (IMiD) are a backbone therapy for multiple myeloma (MM). Despite their efficacy, most patients develop resistance, and the mechanisms are not fully defined. Here, we show that IMiD responses are directed by IMiD-dependent degradation of IKZF1 and IKZF3 that bind to enhancers necessary to sustain the expression of MYC and other myeloma oncogenes. IMiD treatment universally depleted chromatin-bound IKZF1, but eviction of P300 and BRD4 coactivators only occurred in IMiD-sensitive cells. IKZF1-bound enhancers overlapped other transcription factor binding motifs, including ETV4. Chromatin immunoprecipitation sequencing showed that ETV4 bound to the same enhancers as IKZF1, and ETV4 CRISPR/Cas9-mediated ablation resulted in sensitization of IMiD-resistant MM. ETV4 expression is associated with IMiD resistance in cell lines, poor prognosis in patients, and is upregulated at relapse. These data indicate that ETV4 alleviates IKZF1 and IKZF3 dependency in MM by maintaining oncogenic enhancer activity and identify transcriptional plasticity as a previously unrecognized mechanism of IMiD resistance. SIGNIFICANCE: We show that IKZF1-bound enhancers are critical for IMiD efficacy and that the factor ETV4 can bind the same enhancers and substitute for IKZF1 and mediate IMiD resistance by maintaining MYC and other oncogenes. These data implicate transcription factor redundancy as a previously unrecognized mode of IMiD resistance in MM. See related article by Welsh, Barwick, et al., p. 34. See related commentary by Yun and Cleveland, p. 5. This article is featured in Selected Articles from This Issue, p. 4.

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.

Tumor burden limits bispecific antibody efficacy through T cell exhaustion averted by concurrent cytotoxic therapy.

BCMA-CD3-targeting bispecific antibodies (BsAb) are a recently developed immunotherapy class which shows potent tumor killing activity in multiple myeloma (MM). Here, we investigated a murine BCMA-CD3-targeting BsAb in the immunocompetent Vk*MYC and its IMiD-sensitive derivative Vk*MYChCRBN models of MM. The BCMA-CD3 BsAb was safe and efficacious in a subset of mice, but failed in those with high-tumor burden, consistent with clinical reports of BsAb in leukemia. The combination of BCMA-CD3 BsAb with pomalidomide expanded lytic T cells and improved activity even in IMiD resistant high-tumor burden cases. Yet, survival was only marginally extended due to acute toxicity and T cell exhaustion, which impaired T cell persistence. In contrast, the combination with cyclophosphamide was safe and allowed for a tempered pro-inflammatory response associated with long-lasting complete remission. Concurrent cytotoxic therapy with BsAb actually improved T cell persistence and function, offering a promising approach to patients with a large tumor burden.

Monosomic loss of MIR15A/MIR16-1 is a driver of multiple myeloma proliferation and disease progression.

The most common genetic abnormality in multiple myeloma (MM) is the deletion of chromosome 13, seen in almost half of newly diagnosed patients. Unlike chronic lymphocytic leukemia, where a recurrent minimally deleted region including MIR15A/MIR16-1 has been mapped, the deletions in MM predominantly involve the entire chromosome and no specific driver gene has been identified. Additional candidate loci include RB1 and DIS3, but while biallelic deletion of RB1 is associated with disease progression, DIS3 is a common essential gene and complete inactivation is not observed. The Vk*MYC transgenic mouse model of MM spontaneously acquires del(14), syntenic to human chromosome 13, and Rb1 complete inactivation, but not Dis3 mutations. Taking advantage of this model, we explored the role in MM initiation and progression of two candidate loci on chromosome 13: RB1 and MIR15A/MIR16-1. Monoallelic deletion of Mir15a/Mir16-1 but not Rb1 was sufficient to accelerate the development of monoclonal gammopathy in wildtype mice, and the progression of MM in Vk*MYC mice, resulting in increased expression of Mir15a/Mir16-1 target genes and plasma cell proliferation, which was similarly observed in patients with MM.

Identification of PIKfyve kinase as a target in multiple myeloma.

The cellular cytotoxicity of APY0201, a PIKfyve inhibitor, against multiple myeloma was initially identified in an unbiased in vitro chemical library screen. The activity of APY0201 was confirmed in all 25 cell lines tested and in 40% of 100 ex vivo patient-derived primary samples, with increased activity in primary samples harboring trisomies and lacking t(11;14). The broad anti-multiple myeloma activity of PIKfyve inhibitors was further demonstrated in confirmatory screens and showed the superior potency of APY0201 when compared to the PIKfyve inhibitors YM201636 and apilimod, with a mid-point half maximal effective concentration (EC50) at nanomolar concentrations in, respectively, 65%, 40%, and 5% of the tested cell lines. Upregulation of genes in the lysosomal pathway and increased cellular vacuolization were observed in vitro following APY0201 treatment, although these cellular effects did not correlate well with responsiveness. We confirm that PIKfyve inhibition is associated with activation of the transcription factor EB, a master regulator of lysosomal biogenesis and autophagy. Furthermore, we established an assay measuring autophagy as a predictive marker of APY0201 sensitivity. Overall, these findings indicate promising activity of PIKfyve inhibitors secondary to disruption of autophagy in multiple myeloma and suggest a strategy to enrich for likely responders.

Clinical Presentation and Gene Expression Profiling of Immunoglobulin M Multiple Myeloma Compared With Other Myeloma Subtypes and Waldenström Macroglobulinemia.

PURPOSE: Multiple myeloma (MM) is a clonal bone marrow disease characterized by the neoplastic transformation of differentiated postgerminal B cells. It is a heterogeneous disease both at the genetic level and in terms of clinical outcome. Immunoglobulin M (IgM) MM is a rare subtype of myeloma. Similar to Waldenström macroglobulinemia (WM), patients with MM experience IgM monoclonal gammopathy; however, both diseases are distinct in terms of treatment and clinical behavior. MATERIALS AND METHODS: To shed light on the presentation of IgM MM, its prognosis, and its gene expression profiling, we identified and characterized 21 patients with IgM MM from our database. RESULTS: One of these patients presented with a rare IgM monoclonal gammopathy of undetermined significance that progressed to smoldering myeloma. The median survival of the 21 patients was 4.9 years, which was comparable to a matched group of patients with non-IgM MM with similar myeloma prognostic factors (age, gender, albumin, creatinine, anemia, lactate dehydrogenase, ß2-microglobulin, cytogenetics abnormalities), but much less than the median survival reported for patients with WM (9 years). We identified a cluster of genes that differ in their expression profile between MM and WM and found that the patients with IgM MM displayed a gene expression profile most similar to patients with non-IgM MM, confirming that IgM MM is a subtype of MM that should be differentiated from WM. CONCLUSION: Because the prognosis of IgM MM and WM differ significantly, an accurate diagnosis is essential. Our gene expression model can assist with the differential diagnosis in controversial cases.

Low expression of hexokinase-2 is associated with false-negative FDG-positron emission tomography in multiple myeloma.

18F-Fluorodeoxyglucose (FDG)-positron emission tomography (PET) and diffusion-weighted magnetic resonance imaging with background signal suppression (DWIBS) are 2 powerful functional imaging modalities in the evaluation of malignant plasma cell (PC) disease multiple myeloma (MM). Preliminary observations have suggested that MM patients with extensive disease according to DWIBS may be reported as being disease-free on FDG-PET ("PET false-negative"). The aim of this study was to describe the proportion of PET false-negativity in a representative set of 227 newly diagnosed MM patients with simultaneous assessment of FDG-PET and DWIBS, and to identify tumor-intrinsic features associated with this pattern. We found the incidence of PET false-negativity to be 11%. Neither tumor load-associated parameters, such as degree of bone marrow PC infiltration, nor the PC proliferation rate were associated with this subset. However, the gene coding for hexokinase-2, which catalyzes the first step of glycolysis, was significantly lower expressed in PET false-negative cases (5.3-fold change, P < .001) which provides a mechanistic explanation for this feature. In conclusion, we demonstrate a relevant number of patients with FDG-PET false-negative MM and a strong association between hexokinase-2 expression and this negativity: a finding which may also be relevant for clinical imaging of other hematological cancers.

Search for rare protein altering variants influencing susceptibility to multiple myeloma.

The genetic basis underlying the inherited risk of developing multiple myeloma (MM) is largely unknown. To examine the impact of rare protein altering variants on the risk of developing MM we analyzed high-coverage exome sequencing data on 513 MM cases and 1,569 healthy controls, performing both single variant and gene burden tests. We did not identify any recurrent coding low-frequency alleles (1-5%) with moderate effect that were statistically associated with MM. In a gene burden analysis we did however identify a promising relationship between variation in the marrow kinetochore microtubule stromal gene KIF18A, which plays a role in control mitotic chromosome positioning dynamics, and risk of MM (P =3.6x10-6). Further analysis showed KIF18A displays a distinct pattern of expression across molecular subgroups of MM as well as being associated with patient survival. Our results inform future study design and provide a resource for contextualizing the impact of candidate MM susceptibility genes.

The varied distribution and impact of RAS codon and other key DNA alterations across the translocation cyclin D subgroups in multiple myeloma.

We examined a set of 805 cases that underwent DNA sequencing using the FoundationOne Heme (F1H) targeted sequencing panel and gene expression profiling. Known and likely variant calls from the mutational data were analyzed for significant associations with gene expression defined translocation cyclin D (TC) molecular subgroups. The spectrum of KRAS, NRAS, and BRAF codon mutations varied across subgroups with NRAS mutations at Q61 codon being common in hyperdiploid (HRD) and t(11;14) myeloma while being rare in MMSET and MAF. In addition, the presence of RAS-RAF mutations was inversely associated with NFκB pathway activation in all subgroups excluding MAF. In the MMSET subgroup, cases with low FGFR3 expression frequently had RAS-RAF mutations. Conditional inference tree analysis determined that mutation and homozygous deletion of TP53, CDKN2C, and RB1 were key prognostic factors associated with adverse outcome in a non-relapse clinical setting. In conclusion, this study highlights the heterogeneity in the distribution and clinical outcomes of RAS codon and other mutations in multiple myeloma dependent upon primary molecular subgroup.

Clonal selection and double-hit events involving tumor suppressor genes underlie relapse in myeloma.

To elucidate the mechanisms underlying relapse from chemotherapy in multiple myeloma, we performed a longitudinal study of 33 patients entered into Total Therapy protocols investigating them using gene expression profiling, high-resolution copy number arrays, and whole-exome sequencing. The study illustrates the mechanistic importance of acquired mutations in known myeloma driver genes and the critical nature of biallelic inactivation events affecting tumor suppressor genes, especially TP53, the end result being resistance to apoptosis and increased proliferation rates, which drive relapse by Darwinian-type clonal evolution. The number of copy number aberration changes and biallelic inactivation of tumor suppressor genes was increased in GEP70 high risk, consistent with genomic instability being a key feature of high risk. In conclusion, the study highlights the impact of acquired genetic events, which enhance the evolutionary fitness level of myeloma-propagating cells to survive multiagent chemotherapy and to result in relapse.

Removing batch effects from purified plasma cell gene expression microarrays with modified ComBat.

BACKGROUND: Gene expression profiling (GEP) via microarray analysis is a widely used tool for assessing risk and other patient diagnostics in clinical settings. However, non-biological factors such as systematic changes in sample preparation, differences in scanners, and other potential batch effects are often unavoidable in long-term studies and meta-analysis. In order to reduce the impact of batch effects on microarray data, Johnson, Rabinovic, and Li developed ComBat for use when combining batches of gene expression microarray data. We propose a modification to ComBat that centers data to the location and scale of a pre-determined, 'gold-standard' batch. This modified ComBat (M-Combat) is designed specifically in the context of meta-analysis and batch effect adjustment for use with predictive models that are validated and fixed on historical data from a 'gold-standard' batch. RESULTS: We combined data from MIRT across two batches ('Old' and 'New' Kit sample preparation) as well as external data sets from the HOVON-65/GMMG-HD4 and MRC-IX trials into a combined set, first without transformation and then with both ComBat and M-ComBat transformations. Fixed and validated gene risk signatures developed at MIRT on the Old Kit standard (GEP5, GEP70, and GEP80 risk scores) were compared across these combined data sets. Both ComBat and M-ComBat eliminated all of the differences among probes caused by systematic batch effects (over 98% of all untransformed probes were significantly different by ANOVA with 0.01 q-value threshold reduced to zero significant probes with ComBat and M-ComBat). The agreement in mean and distribution of risk scores, as well as the proportion of high-risk subjects identified, coincided with the 'gold-standard' batch more with M-ComBat than with ComBat. The performance of risk scores improved overall using either ComBat or M-Combat; however, using M-ComBat and the original, optimal risk cutoffs allowed for greater ability in our study to identify smaller cohorts of high-risk subjects. CONCLUSION: M-ComBat is a practical modification to an accepted method that offers greater power to control the location and scale of batch-effect adjusted data. M-ComBat allows for historical models to function as intended on future samples despite known, often unavoidable systematic changes to gene expression data.

Artesunate overcomes drug resistance in multiple myeloma by inducing mitochondrial stress and non-caspase apoptosis.

Although novel drugs have contributed immensely to improving outcomes of patients with multiple myeloma (MM), many patients develop drug resistance and ultimately succumb to MM. Here, we show that artesunate, an anti-malarial drug, reliably induces cell death in vitro in naïve as well as drug-resistant MM cells at concentrations shown to be safe in humans. Artesunate induced apoptosis predominantly through the non-caspase mediated pathway by primarily targeting mitochondria and causing outer mitochondrial membrane permeabilization that led to cytosolic and subsequent nuclear translocation of mitochondrial proteins apoptosis inducing factor (AIF) and endonuclease G (EndoG). Nuclear translocation of AIF and EndoG was accompanied by low levels of reactive oxygen species (ROS) and increased mitochondrial production of superoxide. These effects were present before apoptosis was evident and were related to intracellular levels of bivalent iron (Fe+2). Artesunate's unique mechanism probably was at least partially responsible for, its ability to act synergistically with multiple anti-myeloma agents. Our findings suggest that artesunate acts through iron to affect the mitochondria and induce low ROS and non-caspase-mediated apoptosis. Its potency, toxicity profile, and synergism with other drugs make it an intriguing new candidate for MM treatment.