A p53 score derived from TP53 CRISPR/Cas9 HMCLs predicts survival and reveals a major role of BAX in the response to BH3 mimetics.

ABSTRACT: To establish a strict p53-dependent gene-expression profile, TP53-/- clones were derived from TP53+/+ and TP53-/mut t(4;14) human myeloma cell lines (HMCLs) using CRISPR/Cas9 technology. From the 17 dysregulated genes shared between the TP53-/- clones from TP53+/+ HMCLs, we established a functional p53 score, involving 13 genes specifically downregulated upon p53 silencing. This functional score segregated clones and myeloma cell lines as well as other cancer cell lines according to their TP53 status. The score efficiently identified samples from patients with myeloma with biallelic TP53 inactivation and was predictive of overall survival in Multiple Myeloma Research Foundation-coMMpass and CASSIOPEA cohorts. At the functional level, we showed that among the 13 genes, p53-regulated BAX expression correlated with and directly affected the MCL1 BH3 mimetic S63845 sensitivity of myeloma cells by decreasing MCL1-BAX complexes. However, resistance to S63845 was overcome by combining MCL1 and BCL2 BH3 mimetics, which displayed synergistic efficacy. The combination of BH3 mimetics was effective in 97% of patient samples with or without del17p. Nevertheless, single-cell RNA sequencing analysis showed that myeloma cells surviving the combination had lower p53 score, showing that myeloma cells with higher p53 score were more sensitive to BH3 mimetics. Taken together, we established a functional p53 score that identifies myeloma cells with biallelic TP53 invalidation, demonstrated that p53-regulated BAX is critical for optimal cell response to BH3 mimetics, and showed that MCL1 and BCL2 BH3 mimetics in combination may be of greater effectiveness for patients with biallelic TP53 invalidation, for whom there is still an unmet medical need.

High-dose melphalan treatment significantly increases mutational burden at relapse in multiple myeloma.

High-dose melphalan (HDM) improves progression-free survival in multiple myeloma (MM), yet melphalan is a DNA-damaging alkylating agent; therefore, we assessed its mutational effect on surviving myeloma cells by analyzing paired MM samples collected at diagnosis and relapse in the IFM 2009 study. We performed deep whole-genome sequencing on samples from 68 patients, 43 of whom were treated with RVD (lenalidomide, bortezomib, and dexamethasone) and 25 with RVD + HDM. Although the number of mutations was similar at diagnosis in both groups (7137 vs 7230; P = .67), the HDM group had significantly more mutations at relapse (9242 vs 13 383, P = .005). No change in the frequency of copy number alterations or structural variants was observed. The newly acquired mutations were typically associated with DNA damage and double-stranded breaks and were predominantly on the transcribed strand. A machine learning model, using this unique pattern, predicted patients who would receive HDM with high sensitivity, specificity, and positive prediction value. Clonal evolution analysis showed that all patients treated with HDM had clonal selection, whereas a static progression was observed with RVD. A significantly higher percentage of mutations were subclonal in the HDM cohort. Intriguingly, patients treated with HDM who achieved complete remission (CR) had significantly more mutations at relapse yet had similar survival rates as those treated with RVD who achieved CR. This similarity could have been due to HDM relapse samples having significantly more neoantigens. Overall, our study identifies increased genomic changes associated with HDM and provides rationale to further understand clonal complexity.

CARD11 gain of function upregulates BCL2A1 expression and promotes resistance to targeted therapies combination in B-cell lymphoma.

A strategy combining targeted therapies is effective in B-cell lymphomas (BCL), such as mantle cell lymphoma (MCL), but acquired resistances remain a recurrent issue. In this study, we performed integrative longitudinal genomic and single-cell RNA-sequencing analyses of patients with MCL who were treated with targeted therapies against CD20, BCL2, and Bruton tyrosine kinase (OAsIs trial). We revealed the emergence of subclones with a selective advantage against OAsIs combination in vivo and showed that resistant cells were characterized by B-cell receptor (BCR)-independent overexpression of NF-κB1 target genes, especially owing to CARD11 mutations. Functional studies demonstrated that CARD11 gain of function not only resulted in BCR independence but also directly increased the transcription of the antiapoptotic BCL2A1, leading to resistance against venetoclax and OAsIs combination. Based on the transcriptional profile of OAsIs-resistant subclones, we designed a 16-gene resistance signature that was also predictive for patients with MCL who were treated with conventional chemotherapy, underlying a common escape mechanism. Among druggable strategies to inhibit CARD11-dependent NF-κB1 transduction, we evaluated the selective inhibition of its essential partner MALT1. We demonstrated that MALT1 protease inhibition led to a reduction in the expression of genes involved in OAsIs resistance, including BCL2A1. Consequently, MALT1 inhibition induced synergistic cell death in combination with BCL2 inhibition, irrespective of CARD11 mutational status, both in vitro and in vivo. Taken together, our study identified mechanisms of resistance to targeted therapies and provided a novel strategy to overcome resistance in aggressive BCL. The OAsIs trial was registered at www.clinicaltrials.gov #NCT02558816.

del(17p) without TP53 mutation confers a poor prognosis in intensively treated newly diagnosed patients with multiple myeloma.

Despite tremendous improvements in the outcome of patients with multiple myeloma in the past decade, high-risk patients have not benefited from the approval of novel drugs. The most important prognostic factor is the loss of parts of the short arm of chromosome 17, known as deletion 17p (del(17p)). A recent publication (on a small number of patients) suggested that these patients are at very high-risk only if del(17p) is associated with TP53 mutations, the so-called "double-hit" population. To validate this finding, we designed a much larger study on 121 patients presenting del(17p) in > 55% of their plasma cells, and homogeneously treated by an intensive approach. For these 121 patients, we performed deep next generation sequencing targeted on TP53. The outcome was then compared with a large control population (2505 patients lacking del(17p)). Our results confirmed that the "double hit" situation is the worst (median survival = 36 months), but that del(17p) alone also confers a poor outcome compared with the control cohort (median survival = 52.8 months vs 152.2 months, respectively). In conclusion, our study clearly confirms the extremely poor outcome of patients displaying "double hit," but also that del(17p) alone is still a very high-risk feature, confirming its value as a prognostic indicator for poor outcome.

Whole genome copy number analysis in search of new prognostic biomarkers in first line treatment of mantle cell lymphoma. A study by the LYSA group.

Mantle cell lymphoma (MCL) is a lymphoproliferative disorder characterized by the t(11;14)(q13;q32) CCND1/IGH translocation. This lymphoma is however extremely heterogeneous in terms of molecular alterations. Moreover, the course of the disease can vary greatly between indolent forms with slow progression and aggressive conditions rapidly pejorative. The identification of early markers allowing to predict individual patients outcome has however been unsuccessful so far. The LyMa trial treated homogeneously a cohort of young MCL patients. This appeared as a good opportunity to search for biomarkers of response to therapy. DNA extracted from diagnostic paraffin-embedded lymph node biopsies from 100 patients with newly diagnosed MCL, homogeneously treated in this prospective clinical trial, were investigated for copy number alterations and copy neutral loss of heterozygosity using the Oncoscan SNP-array scanning the whole genome. An independent confirmatory cohort was used to strengthen the possibly relevant anomalies observed. Here we describe the recurrent anomalies identified with this technique. Deletions of 17p(TP53) and 9p(CDKN2A) were more frequent in refractory or early relapsing patients (10%), but had no significant impact in univariate analysis on progression-free (PFS) or overall survival (OS). Regardless of the presence of TP53 or CDKN2A deletions, gains in 7p22 (8,5%) were associated with better PFS in univariate but not in multivariate analysis including MCL International Prognostic Index and treatment. Gains of 11q(CCDN1), suggesting gains of the CCND1/IGH fusion, were associated with worse OS and PFS in univariate and multivariate analyses. This worse prognosis impact was confirmed by FISH in an independent confirmatory cohort. This work, using a whole genome approach, confirms the broad genomic landscape of MCL and shows that gains of the CCND1/IGH fusion can be considered as a new prognostic structural variant. Genomic abnormalities of prognostic impact could be useful to strengthen or de-escalate treatment schedules or choosing targeted therapies or CART-cells.

Prognostic role of circulating exosomal miRNAs in multiple myeloma.

Exosomes, secreted by several cell types, including cancer cells, can be isolated from the peripheral blood and have been shown to be powerful markers of disease progression in cancer. In this study, we examined the prognostic significance of circulating exosomal microRNAs (miRNAs) in multiple myeloma (MM). A cohort of 156 patients with newly diagnosed MM, uniformly treated and followed, was studied. Circulating exosomal miRNAs were isolated and used to perform a small RNA sequencing analysis on 10 samples and a quantitative reverse transcription polymerase chain reaction (qRT-PCR) array on 156 samples. We studied the relationship between miRNA levels and patient outcomes, including progression-free survival (PFS) and overall survival (OS). We identified miRNAs as the most predominant small RNAs present in exosomes isolated from the serum of patients with MM and healthy controls by small RNA sequencing of circulating exosomes. We then analyzed exosomes isolated from serum samples of 156 patients using a qRT-PCR array for 22 miRNAs. Two of these miRNAs, let-7b and miR-18a, were significantly associated with both PFS and OS in the univariate analysis and were still statistically significant after adjusting for the International Staging System and adverse cytogenetics in the multivariate analysis. Our findings support the use of circulating exosomal miRNAs to improve the identification of patients with newly diagnosed MM with poor outcomes. The results require further validation in other independent prospective MM cohorts.

The MYRACLE protocol study: a multicentric observational prospective cohort study of patients with multiple myeloma.

BACKGROUND: Despite recent advances in the treatment of multiple myeloma, the disease constantly relapses and is still considered as incurable. The current knowledge about the biological mechanisms underlying resistance to the different class of drugs in multiple myeloma remains poor. The primary objective of the MYRACLE (Myeloma Resistance And Clonal Evolution) cohort, a multicenter prospective cohort of patients with multiple myeloma, is to address this limitation. We here describe the study background, design and methods used for this cohort. METHODS/DESIGN: All patients (> 18 year old) diagnosed with de novo or relapsed multiple myeloma and treated in two hematology department from west of France are included in the MYRACLE cohort. Patients provide a signed informed to be included in the study. All subjects are followed until refusal to participate in the study or death. The MYRACLE cohort prospectively collects data on socio-economic status, medical status, imaging, prognosis factors, MM therapies and associated events (resistance, safety issues). Patients also complete standardized quality of life questionnaires. In addition, bone marrow samples will be collected at time of diagnosis and relapses to perform biomarkers analysis and functional assays exploring mechanisms underlying drug resistance. DISCUSSION: The "real-life" MYRACLE cohort offers the opportunity to prospectively collect epidemiological, medical, QoL and biological data from MM patients during the course of the disease (at time of diagnosis and subsequent relapses). At mid-tem, this integrative cohort will be unique at producing a large variety of data that can be used to conceive the most effective personalized therapy for MM patients. Additionally, the MYRACLE cohort will allow integrating the medical care of MM patients in a health and pharmacoeconomic perspective.

Prognostic value of involved/uninvolved free light chain ratio determined by Freelite and N Latex FLC assays for identification of high-risk smoldering myeloma patients.

Background Smoldering multiple myeloma (SMM) is an asymptomatic plasma cell disorder with a high risk of progression to symptomatic multiple myeloma (MM). The serum free light chain (sFLC) ratio is a powerful prognostic factor for SMM: an sFLC ratio ≥8 has been reported to be associated with a high risk of progression to MM, and an sFLC ratio ≥100 has been described as a criterion for ultra-high-risk SMM, and has been integrated into the definition criteria for MM since 2014. However, all recommendations were based on sFLC measured using the first commercialized assay, Freelite™, while other assays are now available. We aimed to evaluate the safety and accuracy of N-Latex sFLC to identify high-risk and ultra-high-risk SMM. Methods The sFLC ratio was measured at diagnosis with both Freelite and N-Latex assays in a cohort of 176 SMM patients on a BN Prospec nephelometer. Demographic, clinical, therapeutic and laboratory data were collected at the time of diagnosis and at follow-up. Results Sixty-two patients (35.2%) progressed to MM within 2 years. Compared to Freelite™ sFLC, N Latex sFLC ratios ≥8 and ≥100 provided similar performances for the identification of high-risk and ultra-high risk SMM patients. Conclusions Our results evidenced that the N-Latex assay could be used for SMM monitoring, like Freelite. However, an N-Latex sFLC ratio ≥70 appears to provide similar performances to a Freelite sFLC ratio ≥100, with a slightly better positive predictive value. Both assays provided accurate identification of high-risk and ultra-high risk SMM patients. These results should be confirmed in an independent study.

Role of additional chromosomal changes in the prognostic value of t(4;14) and del(17p) in multiple myeloma: the IFM experience.

In multiple myeloma, cytogenetic changes are important predictors of patient outcome. In this setting, the most important changes are deletion 17p, del(17p), and translocation of chromosomes 4 and 14, t(4;14), conferring a poor outcome. However, a certain degree of heterogeneity is observed in the survival of these high-risk patients. We hypothesized that other chromosomal changes may impact the outcome. We retrospectively analyzed a large series of 242 patients displaying either t(4;14) (157 patients) or del(17p) (110 patients), 25 patients presenting both abnormalities, using single nucleotide polymorphism array. In patients with t(4;14), del(1p32), del22q, and >30 chromosomal structural changes negatively impacted progression-free survival (PFS). For overall survival (OS), del(13q14), del(1p32), and the number of chromosomal structural changes worsened the prognosis of patients. For patients with del(17p), del6q worsened the prognosis of patients, whereas trisomy 15 and monosomy 14 were found to have a protective effect on PFS. For OS, del(1p32) worsened the prognosis of patients, whereas having >8 numerical changes was found to have a protective effect on survival. This study, which is the largest series of high-risk patients analyzed with the most modern genomic technique, identified 1 main factor negatively impacting survival: del(1p32).

Understanding the role of hyperdiploidy in myeloma prognosis: which trisomies really matter?

The prognosis of multiple myeloma is mainly dependent upon chromosomal changes. The 2 major abnormalities driving poor outcome are del(17p) and t(4;14). However, the outcome of these high-risk patients is not absolutely uniform, with some patients presenting long survival. We hypothesized that these better outcomes might be related to concomitant "good-risk" chromosomal changes exploring hyperdiploidy. We analyzed a large series of 965 myeloma patients, including 168 patients with t(4;14) and 126 patients with del(17p), using high-throughput single-nucleotide polymorphism arrays after plasma cell sorting. As expected, trisomic chromosomes were highly associated. Using the LASSO model, we found that only chromosome 3, when trisomic, was associated with a longer progression-free survival and that 3 trisomies modulated overall survival (OS) in myeloma patients: trisomies 3 and 5 significantly improved OS, whereas trisomy 21 worsened OS. In patients with t(4;14), trisomies 3 and/or 5 seemed to overcome the poor prognosis. For the first time, using a specific modeling approach, we show that not all trisomies display the same prognostic impact. This finding could be important for routine assessment of prognosis in myeloma, and some high-risk patients with a traditional evaluation could in fact be standard-risk patients.

Differential and limited expression of mutant alleles in multiple myeloma.

Recent work has delineated mutational profiles in multiple myeloma and reported a median of 52 mutations per patient, as well as a set of commonly mutated genes across multiple patients. In this study, we have used deep sequencing of RNA from a subset of these patients to evaluate the proportion of expressed mutations. We find that the majority of previously identified mutations occur within genes with very low or no detectable expression. On average, 27% (range, 11% to 47%) of mutated alleles are found to be expressed, and among mutated genes that are expressed, there often is allele-specific expression where either the mutant or wild-type allele is suppressed. Even in the absence of an overall change in gene expression, the presence of differential allelic expression within malignant cells highlights the important contribution of RNA-sequencing in identifying clinically significant mutational changes relevant to our understanding of myeloma biology and also for therapeutic applications.

Mechanisms of resistance to bispecific T-cell engagers in multiple myeloma and their clinical implications.

ABSTRACT: Bispecific T-cell engagers (TCEs) are revolutionizing patient care in multiple myeloma (MM). These monoclonal antibodies, that redirect T cells against cancer cells, are now approved for the treatment of triple-class exposed relapsed/refractory MM (RRMM). They are currently tested in earlier lines of the disease, including in first line. Yet, primary resistance occurs in about one-third of patients with RRMM, and most responders eventually develop acquired resistance. Understanding the mechanisms of resistance to bispecific TCE is thus essential to improve immunotherapies in MM. Here, we review recent studies investigating the clinical and molecular determinants of resistance to bispecific TCE. Resistance can arise from tumor-intrinsic or tumor-extrinsic mechanisms. Tumor-intrinsic resistance involves various alterations leading to the loss of the target antigen, such as chromosome deletions, point mutations, or epigenetic silencing. Loss of major histocompatibility complex (MHC) class I, preventing MHC class I: T-cell receptor (TCR) costimulatory signaling, was also reported. Tumor-extrinsic resistance involves abundant exhausted T-cell clones and several factors generating an immunosuppressive microenvironment. Importantly, some resistance mechanisms impair response to 1 TCE while preserving the efficacy of others. We next discuss the clinical implications of these findings. Monitoring the status of target antigens in tumor cells and their immune environment will be key to select the most appropriate TCE for each patient and to design combination and sequencing strategies for immunotherapy in MM.

The shaping and functional consequences of the dosage effect landscape in multiple myeloma.

BACKGROUND: Multiple myeloma (MM) is a malignant proliferation of plasma B cells. Based on recurrent aneuploidy such as copy number alterations (CNAs), myeloma is divided into two subtypes with different CNA patterns and patient survival outcomes. How aneuploidy events arise, and whether they contribute to cancer cell evolution are actively studied. The large amount of transcriptomic changes resultant of CNAs (dosage effect) pose big challenges for identifying functional consequences of CNAs in myeloma in terms of specific driver genes and pathways. In this study, we hypothesize that gene-wise dosage effect varies as a result from complex regulatory networks that translate the impact of CNAs to gene expression, and studying this variation can provide insights into functional effects of CNAs. RESULTS: We propose gene-wise dosage effect score and genome-wide karyotype plot as tools to measure and visualize concordant copy number and expression changes across cancer samples. We find that dosage effect in myeloma is widespread yet variable, and it is correlated with gene expression level and CNA frequencies in different chromosomes. Our analysis suggests that despite the enrichment of differentially expressed genes between hyperdiploid MM and non-hyperdiploid MM in the trisomy chromosomes, the chromosomal proportion of dosage sensitive genes is higher in the non-trisomy chromosomes. Dosage-sensitive genes are enriched by genes with protein translation and localization functions, and dosage resistant genes are enriched by apoptosis genes. These results point to future studies on differential dosage sensitivity and resistance of pro- and anti-proliferation pathways and their variation across patients as therapeutic targets and prognosis markers. CONCLUSIONS: Our findings support the hypothesis that recurrent CNAs in myeloma are selected by their functional consequences. The novel dosage effect score defined in this work will facilitate integration of copy number and expression data for identifying driver genes in cancer genomics studies. The accompanying R code is available at http://www.canevolve.org/dosageEffect/.

Identification of human leucocyte antigen (HLA)-A*0201-restricted cytotoxic T lymphocyte epitopes derived from HLA-DOß as a novel target for multiple myeloma.

Despite the recent development of effective therapeutic agents against multiple myeloma (MM), new therapeutic approaches, including immunotherapies, remain to be developed. Here we identified novel human leucocyte antigen (HLA)-A*0201 (HLA-A2)-restricted cytotoxic T lymphocyte (CTL) epitopes from a B cell specific molecule HLA-DOß (DOB) as a potential target for MM. By DNA microarray analysis, the HLA-DOB expression in MM cells was significantly higher than that in normal plasma cells. Twenty-five peptides were predicted to bind to HLA-A2 from the amino acid sequence of HLA-DOB. When screened for the immunogenicity in HLA-A2-transgenic mice immunized with HLA-DOB cDNA, 4 peptides were substantially immunogenic. By mass spectrometry analysis of peptides eluted from HLA-A2-immunoprecipitates of MM cell lines, only two epitopes, HLA-DOB232-240 (FLLGLIFLL) and HLA-DOB185-193 (VMLEMTPEL), were confirmed for their physical presence on cell surface. When healthy donor blood was repeatedly stimulated in vitro with these two peptides and assessed by antigen-specific γ-interferon secretion, HLA-DOB232-240 was more immunogenic than HLA-DOB185-193 . Additionally, the HLA-DOB232-240 -specific CTLs, but not the HLA-DOB185-193 -specific CTLs, displayed an major histocompatibility complex class I-restricted reactivity against MM cell lines expressing both HLA-A2 and HLA-DOB. Taken together, based on the physical presence on tumour cell surface and high immunogenicity, HLA-DOB232-240 might be useful for developing a novel immunotherapy against MM.

Chromothripsis identifies a rare and aggressive entity among newly diagnosed multiple myeloma patients.

Multiple myeloma (MM) develops from a premalignant plasma cell proliferative disorder, and with time can progress to a more aggressive disease in extramedullary locations. The gradually clinical evolution is supported by clonal expansion of cells that acquire genetic lesions over years. This model of cancer evolution based on ongoing genomic instability mechanism may apply to development of most MM cases. However, in a small fraction of newly diagnosed MM who relapse quickly and finally die within 2 years, the gradual model appears to be untenable. Analysis of high resolution copy number profiles obtained using single nucleotide polymorphism array data from 764 newly diagnosed MM identified large numbers of genomic rearrangements with the hallmarks of chromothripsis in 1.3% of samples. Moreover, this catastrophic event confers a poor outcome. Because chromothripsis appears to occur in a single crisis, our results suggest that high-risk MM patients use this novel way of cancer evolution.

Translocation t(14;16) and multiple myeloma: is it really an independent prognostic factor?

Many trials in myeloma are stratified on cytogenetic abnormalities. Among them, the most commonly chosen are the t(4;14), the del(17p), and the t(14;16). If data are well established for t(4;14) and del(17p), very few data support the use of t(14;16). To address this issue, we retrospectively analyzed 1003 patients with newly diagnosed myeloma for this abnormality. We identified 32 patients with the t(14;16). Compared with patients lacking the t(14;16), we did not observe any difference in overall survival (P = .28). Moreover, in multivariate analyses, the t(14;16) was not prognostic (P = .39). In conclusion, our data do not support the use of t(14;16)-specific probes in the diagnostic panels of multiple myeloma.

Acquired resistance to a GPRC5D-directed T-cell engager in multiple myeloma is mediated by genetic or epigenetic target inactivation.

Bispecific antibodies targeting GPRC5D demonstrated promising efficacy in multiple myeloma, but acquired resistance usually occurs within a few months. Using a single-nucleus multi-omic strategy in three patients from the MYRACLE cohort (ClinicalTrials.gov registration: NCT03807128 ), we identified two resistance mechanisms, by bi-allelic genetic inactivation of GPRC5D or by long-range epigenetic silencing of its promoter and enhancer regions. Molecular profiling of target genes may help to guide the choice of immunotherapy and early detection of resistance in multiple myeloma.

Exploring the impact of dexamethasone on gene regulation in myeloma cells.

Among glucocorticoids (GCs), dexamethasone (Dex) is widely used in treatment of multiple myelomas. However, despite a definite benefit, all patients relapse. Moreover, the molecular basis of glucocorticoid efficacy remains elusive. To determine genomic response to Dex in myeloma cells, we generated bulk and single-cell multi-omics data and high-resolution contact maps of active enhancers and target genes. We show that a minority of glucocorticoid receptor-binding sites are associated with enhancer activity gains, increased interaction loops, and transcriptional activity. We identified and characterized a predominant enhancer enriched in cohesin (RAD21) and more accessible upon Dex exposure. Analysis of four gene-specific networks revealed the importance of the CTCF-cohesin couple and the synchronization of regulatory sequence openings for efficient transcription in response to Dex. Notably, these epigenomic changes are associated with cell-to-cell transcriptional heterogeneity, in particular, lineage-specific genes. As consequences, BCL2L11-encoding BIM critical for Dex-induced apoptosis and CXCR4 protective from chemotherapy-induced apoptosis are rather up-regulated in different cells. In summary, our work provides new insights into the molecular mechanisms involved in Dex escape.

MicroRNAs 15a/16-1 function as tumor suppressor genes in multiple myeloma.

Multiple myeloma is characterized by frequent chromosomal alterations. Deletion of chr 13, especially band 13q14, is commonly observed in early stages of MM, suggesting the presence of tumor suppressor genes within this region. Here, we functionally validate the role of the microRNAs-15a/16-1 cluster, centered at the deleted region, as TSGs and delineate their downstream target genes in MM. Using "sponge" lentiviral vectors to competitive stably inhibit mature microRNAs in vitro and in vivo, we have documented enhanced proliferative and invasive capacity of cells with stably inhibition of miR-16. Importantly, miR-16 inhibition decreased animal survival in a xenograft model of MM by increasing tumor load and host angiogenesis. Expression profiling analysis of miR-16-deficient cells identified a large number of downstream target genes including FGFR1, PI3KCa, MDM4, VEGFa, as well as secondary affected genes such as JUN and Jag1. We validated designated genes showing binding sites within the conserved 3'-UTR and also within the mRNA coding region as direct miR-16 targets, thus indicating that the miRNAs may have many more targets than anticipated by conventional prediction methods. This loss-of-function system, which mimics the 13q chromosomal deletion, provides a valuable tool to investigate their function in MM pathogenesis and their potential use as therapeutic targets.

In-depth analysis of alternative splicing landscape in multiple myeloma and potential role of dysregulated splicing factors.

Splicing changes are common in cancer and are associated with dysregulated splicing factors. Here, we analyzed RNA-seq data from 323 newly diagnosed multiple myeloma (MM) patients and described the alternative splicing (AS) landscape. We observed a large number of splicing pattern changes in MM cells compared to normal plasma cells (NPC). The most common events were alterations of mutually exclusive exons and exon skipping. Most of these events were observed in the absence of overall changes in gene expression and often impacted the coding potential of the alternatively spliced genes. To understand the molecular mechanisms driving frequent aberrant AS, we investigated 115 splicing factors (SFs) and associated them with the AS events in MM. We observed that ~40% of SFs were dysregulated in MM cells compared to NPC and found a significant enrichment of SRSF1, SRSF9, and PCB1 binding motifs around AS events. Importantly, SRSF1 overexpression was linked with shorter survival in two independent MM datasets and was correlated with the number of AS events, impacting tumor cell proliferation. Together with the observation that MM cells are vulnerable to splicing inhibition, our results may lay the foundation for developing new therapeutic strategies for MM. We have developed a web portal that allows custom alternative splicing event queries by using gene symbols and visualizes AS events in MM and subgroups. Our portals can be accessed at http://rconnect.dfci.harvard.edu/mmsplicing/ and https://rconnect.dfci.harvard.edu/mmleafcutter/ .