Search

The bone marrow (BM) microenvironment actively promotes multiple myeloma pathogenesis, and therapies targeting both cancer cells and the niche are highly effective. We were interested in identifying novel signaling pathways supporting multiple myeloma–BM cross-talk. Mutations in the transmembrane receptor Roundabout 1 (ROBO1) were recently identified in patients with multiple myeloma; however, their functional consequences are uncertain. Through protein structure–function studies, we discovered that ROBO1 is necessary for multiple myeloma adhesion to BM stromal and endothelial cells and that ROBO1 knockout (KO) compromises BM homing and engraftment in a disseminated mouse model. ROBO1 KO significantly decreases multiple myeloma proliferation in vitro and intra- and extramedullary tumor growth in vivo. Mechanistically, the ROBO1 C-terminus is cleaved in a ligand-independent fashion and is sufficient to promote multiple myeloma proliferation. Vice versa, mutants lacking the cytoplasmic domain, including the human-derived G674* truncation, act dominantly negative. Interactomic and RNA-sequencing studies suggest that ROBO1 may be involved in RNA processing, supporting further studies.Significance:ROBO1 is highly expressed in t(4;14) multiple myeloma and supports homing and dissemination to the BM niche. ROBO1 knockout causes reduced tumor growth in intramedullary and extramedullary myeloma animal models, while the ROBO1 C-terminus is cleaved in multiple fragments and it is necessary and sufficient to sustain myeloma proliferation.

Supplementary Figure 1 showing ROBO1 expression in primary MM cells versus normal donor bone marrow plasma cells via gene expression profiling and immunohistochemistry. Supplementary Figure 2 showing distribution of ROBO1 expression in the CoMMpass database; ROBO1 expression in CoMMpass database according to cytogenetics; ROBO1 expression in IFM170 database according to cytogenetics; ROBO1 expression in t(4;14) cell lines WT or KO for MMSET. Supplementary Figure 3 showing ROBO1 protein down regulation after ROBO1 knock down via shRNA. Supplementary Figure 4 showing effects of shRNA4- and shRNA-5 mediated ROBO1 KD on cell viability across a panel of MM cell lines with western confirming proper protein down regulation. Supplementary Figure 5 showing successful ROBO1 KO across several KMS11 monoclones, but not H929 monoclones. Supplementary Figure 6 showing expression of ROBO1 FL construct in polyclonally and monoclonally-transduced OPM2 cells. Supplementary Figure 7 showing expression of FL ROBO1 addback in distinct KMS11 ROBO1 KO monoclones. Supplementary Figure 8 showing radiographic images of WT versus KO ROBO1 mu-SCID mice. Supplementary Figure 9 showing adhesion defect in ROBO1 KO KMS11 cells that is fully rescued by FL ROBO1 addback. Supplementary Figure 10 showing Kaplan-Meier survival curve of disseminated mouse model injected with ROBO1 KO versus FL OPM2 cells. Supplementary Figure 11 showing cleaved cytosolic ROBO1 fragments in 293T cells transduced with FL ROBO1 tagged with C-LAP tag and in WT MM.1S cells. Supplementary Figure 12 showing interacting partners of ROBO1 cytosolic domain as deconvoluted via STRING database. Supplementary Figure 13 showing differentially expressed genes in ROBO1 WT versus KO plasmacytoma represented via Volcano plot. Supplementary Figure 14 showing gene set enrichment analysis of differentially expressed genes in ROBO1 WT versus KO plasmacytoma. Supplementary Figure 15 showing schema of patient-derived ROBO1 mutation and expression of G674 truncation in WT and KO ROBO1 OPM2 cells. Supplementary Table 1 showing clinical annotation of newly diagnosed MM parents used in Supplementary figure 1.

Supplementary Methods

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.

Immunoglobulin M (IgM) multiple myeloma (MM) is a rare disease subgroup. Its differentiation from other IgM-producing gammopathies such as Waldenström macroglobulinemia (WM) has not been well characterized but is essential for proper risk assessment and treatment. In this study, we investigated genomic and transcriptomic characteristics of IgM-MM samples using whole-genome and transcriptome sequencing to identify differentiating characteristics from non–IgM-MM and WM. Our results suggest that IgM-MM shares most of its defining structural variants and gene-expression profiling with MM, but has some key characteristics, including t(11;14) translocation, chromosome 6 and 13 deletion as well as distinct molecular and transcription-factor signatures. Furthermore, IgM-MM translocations were predominantly characterized by VHDHJH recombination-induced breakpoints, as opposed to the usual class-switching region breakpoints; coupled with its lack of class switching, these data favor a pre–germinal center origin. Finally, we found elevated expression of clinically relevant targets, including CD20 and Bruton tyrosine kinase, as well as high BCL2/BCL2L1 ratio in IgM-MM, providing potential for targeted therapeutics.

Long noncoding RNAs (lncRNAs) can drive tumorigenesis and are susceptible to therapeutic intervention. Here, we used a large-scale CRISPR interference viability screen to interrogate cell-growth dependency to lncRNA genes in multiple myeloma (MM) and identified a prominent role for the miR-17-92 cluster host gene (MIR17HG). We show that an MIR17HG-derived lncRNA, named lnc-17-92, is the main mediator of cell-growth dependency acting in a microRNA- and DROSHA-independent manner. Lnc-17-92 provides a chromatin scaffold for the functional interaction between c-MYC and WDR82, thus promoting the expression of ACACA, which encodes the rate-limiting enzyme of de novo lipogenesis acetyl-coA carboxylase 1. Targeting MIR17HG pre-RNA with clinically applicable antisense molecules disrupts the transcriptional and functional activities of lnc-17-92, causing potent antitumor effects both in vitro and in vivo in 3 preclinical animal models, including a clinically relevant patient-derived xenograft NSG mouse model. This study establishes a novel oncogenic function of MIR17HG and provides potent inhibitors for translation to clinical trials.

SUMMARYLong noncoding RNAs (lncRNA) can drive the tumorigenesis and be susceptible to therapeutic intervention. To define the landscape of therapeutically actionable lncRNA dependencies in multiple myeloma (MM), we coupled our extensive lncRNA transcriptomic profile with lncRNA targeted CRISPR interference viability screen and identified RNA Regulator of Lipogenesis (RROL) as a leading lncRNA dependency in MM. RROL shares its origin with the microRNA locus MIR17HG, however supports the proliferation and survival of MM cells in a microRNA- and DROSHA- independent manner. We found that RROL provides a chromatin scaffold for the functional interaction between c-MYC and WDR82 to promote the regulation of the lipogenic pathways via the transcriptional control of the rate-limiting enzyme ACC1 in MM cells. Inhibition of RROL with clinically applicable antisense molecules disrupts its transcriptional and functional activities causing potent anti-tumor effects both in vitro and in vivo in two pre-clinical animal models. This study establishes lncRNA RROL as a therapeutically actionable dependency with a unique mechanism of action in support of myeloma cell growth.

The bone marrow (BM) microenvironment actively promotes multiple myeloma pathogenesis, and therapies targeting both cancer cells and the niche are highly effective. We were interested in identifying novel signaling pathways supporting multiple myeloma–BM cross-talk. Mutations in the transmembrane receptor Roundabout 1 (ROBO1) were recently identified in patients with multiple myeloma; however, their functional consequences are uncertain. Through protein structure–function studies, we discovered that ROBO1 is necessary for multiple myeloma adhesion to BM stromal and endothelial cells and that ROBO1 knockout (KO) compromises BM homing and engraftment in a disseminated mouse model. ROBO1 KO significantly decreases multiple myeloma proliferation in vitro and intra- and extramedullary tumor growth in vivo. Mechanistically, the ROBO1 C-terminus is cleaved in a ligand-independent fashion and is sufficient to promote multiple myeloma proliferation. Vice versa, mutants lacking the cytoplasmic domain, including the human-derived G674* truncation, act dominantly negative. Interactomic and RNA-sequencing studies suggest that ROBO1 may be involved in RNA processing, supporting further studies. Significance: ROBO1 is highly expressed in t(4;14) multiple myeloma and supports homing and dissemination to the BM niche. ROBO1 knockout causes reduced tumor growth in intramedullary and extramedullary myeloma animal models, while the ROBO1 C-terminus is cleaved in multiple fragments and it is necessary and sufficient to sustain myeloma proliferation.

Introduction Long noncoding RNAs (lncRNAs) are profoundly dysregulated in multiple myeloma (MM), yet their tumor-promoting features remain to be defined. Here, we have coupled extensive transcriptomic profiling of patient samples (n=360) with a genome-wide CRISPR interference (CRISPRi) viability screen (4 cell lines) to generate a map of lncRNA "dependencies" in MM; moreover, we describe the molecular and functional role – as well as its therapeutic potential – of the screen top hit RNA Regulator of Lipogenesis (RROL), a nuclear-retained lncRNA generated by alternative splicing of MIR17HG. Methods Loss-of-function (LOF) studies using antisense oligonucleotides (ASO) confirmed the RROL dependency in a large panel of MM cell lines and CD138+ patient cells, in vitro and in vivo in animal models. This effect was not antagonized in DROSHA KO MM cells or by ectopic expression of MIR17HG-derived miR-17-92, indicating a microRNA-independent role of RROL. Results Transcriptomic analysis after RROL depletion indicated a significant impact on the de novo lipogenesis (DNL) gene networks. RROL occupancy at DNL gene loci was detected by chromatin isolation by RNA precipitation followed by qRT-PCR (ChIRP-qPCR) and confirmed by DUAL RNA FISH analysis. At functional level, we demonstrated that RROL depletion reduces DNL in MM cells using an unbiased lipidomic profiling as well as by measuring the incorporation of C14-radiolabeled glucose into lipids. Importantly, we also proved that exogenous palmitate, the downstream product of DNL pathway, can significantly rescue the growth inhibitory effect of RROL depletion in MM cells; thus implicating a role for DNL in the growth promoting effect of this lncRNA. RNA-Protein Pull Down (RPPD) and in vivo RNA yeast three-hybrid (Y3H) assays led to identify c-MYC (MYC) as a relevant protein interactor of RROL; a finding validated by RIP-qPCR. Coupling RNA FISH with immunofluorescence, we co-detected RROL and MYC at the ACC1 promoter. Importantly, neither RROL or MYC could occupy ACC1 promoter and exert regulatory control in the absence of the other factor. Moreover, using in vitro (Co-IP/MS) and in vivo (BioID) assays, we identified WDR82 as an RROL-dependent MYC partner implicated in the transcriptional control of ACC1 expression. Finally, to therapeutically antagonize RROL, we have screened >100 ASOs and provided the basis to develop a first-in-class RROL inhibitor. This inhibitor has shown a very strong anti-MM activity in vitro (IC50 Conclusion In conclusion, we here report a unique regulatory function of a novel lncRNA supporting MM cell growth via its control of the lipogenic metabolic axis. The ongoing development of RROL inhibitors may allow clinical application of this unique targeted therapy in MM.

Background Multiple Myeloma (MM) is a disease driven by numerous genetic and epigenetic alterations, however true drivers of the disease have yet to be identified. Methods To identify new dependencies and actionable therapeutic targets in MM we integrated gene expression and genetic dependency (CRISPR KO). Results We found that many of the specific and potent dependencies in MM are transcription factors, especially those establishing plasma cell identity. Among others, the POU2AF1 gene, which encodes the OCA-B/BOB-1, a B cell transcriptional coactivator protein, represented the most striking dependency in MM. Although BOB-1 is expressed throughout B-cell development, we found it to be highly expressed in CD138+ plasma cells from patients with precursor conditions (MGUS and SMM) as well and established MM compared to normal plasma cells (NPC). Loss-of-function studies using shRNA, siRNAs as well as antisense GapMers specific for BOB-1 confirmed significant impact on MM cell viability. Transcriptomic analysis by RNA-sequencing revealed a small set of genes commonly modulated in MM cell lines upon BOB-1 depletion, including the XBP1- BHLHA15 axis involved in lipid synthesis and unfolded protein response (UPR). Interestingly, among the genes most significantly upregulated by BOB-1 depletion was heme oxygenase 1 (HMOX1), that was affected via the NRF2/Keap1 pathway. We observed that HMOX1 expression is significantly lower in MM cells from patients compared to normal plasma cells and correlates with poor clinical outcome, suggesting important role in MM. Moreover, we found that siRNA depletion of HMOX1 reverted the inhibition of MM cell growth caused by BOB-1 KD, confirming significant role for HMOX1 in the BOB-1 addiction observed in MM cells. Next, we performed gene set enrichment analysis (GSEA) and observed ribosome biogenesis pathways and mRNA translation and elongation processes, along with WNT and senescence pathways, to be significantly enriched among genes modulated by BOB-1 depletion in MM cells. Since high protein load is a feature of MM, we evaluated the role of BOB-1 in the translational efficiency of MM cells. In MM cell lines, BOB-1 knockdown decreased de novo protein synthesis, while its overexpression significantly enhances protein synthesis compared to control cells. As MM is characterized by excess production of monoclonal immunoglobulins, we evaluated impact of BOB-1 perturbation on intracellular kappa and lambda light chains production. We observed changes in the intracellular abundance of the light chains with BOB-1 modulation in all MM cell lines tested. As a result, BOB-1 depletion was associated with induction of resistance to proteasome inhibition. Conclusion In conclusion, we report BOB1 as a specific dependency in MM cells with potential role on modulating the protein load/capacity balance in MM cells and therefore the sensitivity to proteasome inhibition.

BCMA targeting chimeric antigen receptor (CAR) T cell therapy has shown deep and durable responses in multiple myeloma. However, relapse following therapy is frequently observed, and mechanisms of resistance remain ill-defined. Here, we perform single cell genomic characterization of longitudinal samples from a patient who relapsed after initial CAR T cell treatment with lack of response to retreatment. We report selection, following initial CAR T cell infusion, of a clone with biallelic loss of BCMA acquired by deletion of one allele and a mutation that creates an early stop codon on the second allele. This loss leads to lack of CAR T cell proliferation following the second infusion and is reflected by lack of soluble BCMA in patient serum. Our analysis suggests the need for careful detection of BCMA gene alterations in multiple myeloma cells from relapse following CAR T cell therapy., Relapse following BCMA targeted CAR T-cell therapy is frequently observed in patients with multiple myeloma (MM). Here, by single cell transcriptome profiling on serially collected bone marrow samples, the authors report biallelic loss of BCMA as the mechanism of resistance underlying both relapse and lack of response to a second CAR T infusion in a patient with MM.