Biological and prognostic impact of apobec-induced mutations in the spectrum of plasma cell dyscrasias

In multiple myeloma (MM), whole exome sequencing (WES) studies have revealed four mutational signatures: two associated with aberrant activities of APOBEC cytidine deaminases (Signatures #2 and #13) and two clock-like signatures associated with "cancer age" (Signatures #1 and #5). Mutational signatures have not been investigated systematically in larger series, nor in other primary plasma cell dyscrasias such as monoclonal gammopathy of unknown significance (MGUS) or primary plasma cell leukemia (pPCL). Finally, while APOBEC activity has been correlated to increased mutational burden and poor-prognosis MAF/MAFB translocations in MM at diagnosis, this has never been confirmed in multivariate analysis in an independent series. To answer these questions, we mined 1151 MM samples from public WES datasets, including samples from the IA9 public release of the CoMMpass trial. The CoMMpass data were generated as part of the Multiple Myeloma Research Foundation Personalized Medicine Initiatives. We also analyzed 6 MGUS/Smoldering MM as well as 5 previously published pPCLs. Extraction of mutational signatures was performed using the NNMF algorithm as previously described (Alexandrov et al. Nature 2013). NNMF in the whole cohort extracted the known 4 signatures pertaining to distinct mutational processes: the two clock-like processes (signatures #1 and #5) and aberrant APOBEC deaminase activity (signatures #2 and #13). While the clock-like processes were more prominent in the cohort as a whole (median 70%, range 0-100%), the APOBEC showed a heterogeneous contribution, more visible in samples with the highest mutation burden. In fact, the absolute and relative contribution of APOBEC activity to the mutational repertoire correlated with the overall number of mutations (r=0.71, p= < 0.0001). As previously described, APOBEC contribution was significantly enriched among MM patients with t(14;16) and with t(14;20) (p<0.001), but the association between relative APOBEC contribution and mutational load remained significant across all cytogenetic subgroups with the exception of t(11;14). In the MGUS/SMM series, APOBEC contribution was generally low. Conversely, APOBEC activity was preponderant in three out of five pPCL samples, all of them characterized by the t(14;16)( IGH / MAF); in the remaining two pPCL the absolute number of APOBEC mutations was similar to MM. Overall, the APOBEC contribution was characterized by a progressive increment from MGUS/SMM to MM and pPCL. We next went on to investigate the prognostic impact of APOBEC signatures at diagnosis. Patients with APOBEC contribution in the 4th quartile had shorter PFS (2-y PFS 47% vs 66%, p<0.0001) and OS (2-y OS 70% vs 85%, p=0.0033) than patients in quartiles 1-3 (Figure 1a-b). This was independent from the association of APOBEC activity with MAF translocations and higher mutational burden, as shown by multivariate analysis with Cox regression (Figure 1c-d). ISS stage III was the only other variable that retained its independent prognostic value for both PFS and OS. We therefore combined both variables and found that co-occurrence of ISS III and APOBEC 4th quartile identifies a fraction of high-risk patients with 2-y OS of 53.8% (95% CI 36.6%-79%), while their simultaneous absence identifies long term survivors with 2-y OS of 93.3% (95% CI 89.6-97.2%). In this study, we provided a global overview on the contribution of mutational processes in the largest whole exome series of plasma cell dyscrasias investigated to date by NNMF. We propose that cases with high APOBEC activity may represent a novel prognostic subgroup that is transversal to conventional cytogenetic subgroups, advocating for closer integration of next-generation sequencing studies and clinical annotation to confirm this finding in independent series

A novel 3D mesenchymal stem cell model of the multiple myeloma bone marrow niche : biologic and clinical applications

Specific niches within the tumor bone marrow (BM) microenvironment afford a sanctuary for multiple myeloma (MM) clones due to stromal cell-tumor cell interactions, which confer survival advantage and drug resistance. Defining the sequelae of tumor cell interactions within the MM niches on an individualized basis may provide the rationale for personalized therapies. To mimic the MM niche, we here describe a new 3D co-culture ex-vivo model in which primary MM patient BM cells are co-cultured with mesenchymal stem cells (MSC) in a hydrogel 3D system. In the 3D model, MSC with conserved phenotype (CD73+CD90+CD105+) formed compact clusters with active fibrous connections, and retained lineage differentiation capacity. Extracellular matrix molecules, integrins, and niche related molecules including N-cadherin and CXCL12 are expressed in 3D MSC model. Furthermore, activation of osteogenesis (MMP13, SPP1, ADAMTS4, and MGP genes) and osteoblastogenic differentiation was confirmed in 3D MSC model. Co-culture of patient-derived BM mononuclear cells with either autologous or allogeneic MSC in 3D model increased proliferation of MM cells, CXCR4 expression, and SP cells. We carried out immune profiling to show that distribution of immune cell subsets was similar in 3D and 2D MSC model systems. Importantly, resistance to novel agents (IMiDs, bortezomib, carfilzomib) and conventional agents (doxorubicin, dexamethasone, melphalan) was observed in 3D MSC system, reflective of clinical resistance. This 3D MSC model may therefore allow for studies of MM pathogenesis and drug resistance within the BM niche. Importantly, ongoing prospective trials are evaluating its utility to inform personalized targeted and immune therapy in MM

A high-risk, Double-Hit, group of newly diagnosed myeloma identified by genomic analysis

Patients with newly diagnosed multiple myeloma (NDMM) with high-risk disease are in need of new treatment strategies to improve the outcomes. Multiple clinical, cytogenetic, or gene expression features have been used to identify high-risk patients, each of which has significant weaknesses. Inclusion of molecular features into risk stratification could resolve the current challenges. In a genome-wide analysis of the largest set of molecular and clinical data established to date from NDMM, as part of the Myeloma Genome Project, we have defined DNA drivers of aggressive clinical behavior. Whole-genome and exome data from 1273 NDMM patients identified genetic factors that contribute significantly to progression free survival (PFS) and overall survival (OS) (cumulative R2 = 18.4% and 25.2%, respectively). Integrating DNA drivers and clinical data into a Cox model using 784 patients with ISS, age, PFS, OS, and genomic data, the model has a cumlative R2 of 34.3% for PFS and 46.5% for OS. A high-risk subgroup was defined by recursive partitioning using either a) bi-allelic TP53 inactivation or b) amplification (≥4 copies) of CKS1B (1q21) on the background of International Staging System III, comprising 6.1% of the population (median PFS = 15.4 months; OS = 20.7 months) that was validated in an independent dataset. Double-Hit patients have a dire prognosis despite modern therapies and should be considered for novel therapeutic approaches

A novel 3D mesenchymal stem cell model of the multiple myeloma bone marrow niche: biologic and clinical applications

Specific niches within the tumor bone marrow (BM) microenvironment afford a sanctuary for multiple myeloma (MM) clones due to stromal cell-tumor cell interactions, which confer survival advantage and drug resistance. Defining the sequelae of tumor cell interactions within the MM niches on an individualized basis may provide the rationale for personalized therapies. To mimic the MM niche, we here describe a new 3D co-culture ex-vivo model in which primary MM patient BM cells are co-cultured with mesenchymal stem cells (MSC) in a hydrogel 3D system. In the 3D model, MSC with conserved phenotype (CD73+CD90+CD105+) formed compact clusters with active fibrous connections, and retained lineage differentiation capacity. Extracellular matrix molecules, integrins, and niche related molecules including N-cadherin and CXCL12 are expressed in 3D MSC model. Furthermore, activation of osteogenesis (MMP13, SPP1, ADAMTS4, and MGP genes) and osteoblastogenic differentiation was confirmed in 3D MSC model. Co-culture of patient-derived BM mononuclear cells with either autologous or allogeneic MSC in 3D model increased proliferation of MM cells, CXCR4 expression, and SP cells. We carried out immune profiling to show that distribution of immune cell subsets was similar in 3D and 2D MSC model systems. Importantly, resistance to novel agents (IMiDs, bortezomib, carfilzomib) and conventional agents (doxorubicin, dexamethasone, melphalan) was observed in 3D MSC system, reflective of clinical resistance. This 3D MSC model may therefore allow for studies of MM pathogenesis and drug resistance within the BM niche. Importantly, ongoing prospective trials are evaluating its utility to inform personalized targeted and immune therapy in MM

Non-overlapping Control of Transcriptome by Promoter- and Super-Enhancer-Associated Dependencies in Multiple Myeloma

Summary: The relationship between promoter proximal transcription factor-associated gene expression and super-enhancer-driven transcriptional programs are not well defined. However, their distinct genomic occupancy suggests a mechanism for specific and separable gene control. We explored the transcriptional and functional interrelationship between E2F transcription factors and BET transcriptional co-activators in multiple myeloma. We found that the transcription factor E2F1 and its heterodimerization partner DP1 represent a dependency in multiple myeloma cells. Global chromatin analysis reveals distinct regulatory axes for E2F and BETs, with E2F predominantly localized to active gene promoters of growth and/or proliferation genes and BETs disproportionately at enhancer-regulated tissue-specific genes. These two separate gene regulatory axes can be simultaneously targeted to impair the myeloma proliferative program, providing an important molecular mechanism for combination therapy. This study therefore suggests a sequestered cellular functional control that may be perturbed in cancer with potential for development of a promising therapeutic strategy. : Uncontrolled proliferation is a hallmark of tumorigenesis and is associated with perturbed transcriptomic profile. Fulciniti et al. explored the interrelationship between E2F transcription factors and BET transcriptional co-activators in multiple myeloma, reporting the existence of two distinct regulatory axes that can be synergistically targeted to impact myeloma growth and survival. Keywords: multiple myeloma, transcription factor, transcriptional regulation, super-enhancer

Genome-Wide Somatic Alterations in Multiple Myeloma Reveal a Superior Outcome Group

International audiencePURPOSE Multiple myeloma (MM) is accompanied by heterogeneous somatic alterations. The overall goal of this study was to describe the genomic landscape of myeloma using deep whole-genome sequencing (WGS) and develop a model that identifies patients with long survival. METHODS We analyzed deep WGS data from 183 newly diagnosed patients with MM treated with lenalidomide, bortezomib, and dexamethasone (RVD) alone or RVD 1 autologous stem cell transplant (ASCT) in the IFM/DFCI 2009 study (ClinicalTrials.gov identifier: NCT01191060). We integrated genomic markers with clinical data. RESULTS We report significant variability in mutational load and processes within MM subgroups. The timeline of observed activation of mutational processes provides the basis for 2 distinct models of acquisition of mutational changes detected at the time of diagnosis of myeloma. Virtually all MM subgroups have activated DNA repair-associated signature as a prominent late mutational process, whereas APOBEC signature targeting C.G is activated in the intermediate phase of disease progression in high-risk MM. Importantly, we identify a genomically defined MM subgroup (17% of newly diagnosed patients) with low DNA damage (low genomic scar score with chromosome 9 gain) and a superior outcome (100% overall survival at 69 months), which was validated in a large independent cohort. This subgroup allowed us to distinguish patients with low-and high-risk hyperdiploid MM and identify patients with prolongation of progression-free survival. Genomic characteristics of this subgroup included lower mutational load with significant contribution from age-related mutations as well as frequent NRAS mutation. Surprisingly, their overall survival was independent of International Staging System and minimal residual disease status. CONCLUSION This is a comprehensive study identifying genomic markers of a good-risk group with prolonged survival. Identification of this patient subgroup will affect future therapeutic algorithms and research planning

The spectrum of neutrophilic dermatoses associated with monoclonal gammopathy: Association with IgA isotype and inflammatory profile

International audienc

Genomic patterns of progression in smoldering multiple myeloma

International audienceWe analyzed whole genomes of unique paired samples from smoldering multiple mye-loma (SMM) patients progressing to multiple myeloma (MM). We report that the genomic landscape, including mutational profile and structural rearrangements at the smoldering stage is very similar to MM. Paired sample analysis shows two different patterns of progression: a "static progression model", where the subclonal architecture is retained as the disease progressed to MM suggesting that progression solely reflects the time needed to accumulate a sufficient disease burden; and a "spontaneous evolution model", where a change in the subclonal composition is observed. We also observe that activation-induced cytidine dea-minase plays a major role in shaping the mutational landscape of early subclinical phases, while progression is driven by APOBEC cytidine deaminases. These results provide a unique insight into myelomagenesis with potential implications for the definition of smoldering disease and timing of treatment initiation