Boise, Lawrence

Abstract Although the human genome encodes a vast repertoire of noncoding RNAs that regulate gene expression, the noncoding genome remains underexplored due to technical challenges. Specifically, during transcriptomic sequencing data alignment, the overlap between noncoding and coding loci can create ambiguous read alignments that are subsequently discarded from downstream analysis. For this reason, most of the noncoding genome is excluded from standard genomic annotations used for sequencing alignment. To address this challenge and enable concurrent profiling of the coding and noncoding transcriptome, we systematically integrated standard coding (GENCODE) and noncoding (LncBook) genome annotations, preserving coding gene annotations and removing overlapping noncoding regions. The resulting integrated genome annotation expanded the number of annotated noncoding genes from 40,785 to 138,296 while preserving all coding genes and reducing ambiguous read assignment. To evaluate the utility of our integrated genome annotation for uncovering novel, biologically relevant noncoding RNAs (ncRNAs), we realigned CD138-positive bulk RNA-seq (N = 942) and CD138-negative single-cell RNA-seq (N = 478) data from the MMRF CoMMpass study, generating a comprehensive coding-noncoding atlas of the myeloma bone marrow microenvironment with noncoding genes representing 51% of highly variable genes and displaying significant cell type specificity. Tumor expression profiling based on this integrated profiling identified 15 clusters, including two enriched for amp(1q21) or t(4;14) and associated with shorter progression-free survival (PFS). Differential expression and systematic filtering yielded 19 candidate high-risk ncRNAs, including previously uncharacterized ENSG00000310209 , which was associated with poor PFS (HR = 1.141, P = 0.0025), increased IRF4 activity, Wnt pathway activation, CCL5 signaling, and the accumulation of anergic-like CD8+ T cells. These findings establish integrated coding-noncoding analysis as a strategic approach for discovering functional ncRNAs from transcriptomic sequencing data.

Abstract Multiple myeloma (MM) with t(11;14)(CCND1;IGH) remains the only subset sensitive to the BCL2 inhibitor venetoclax. However, not all patients with t(11;14)(CCND1;IGH) respond to treatment, and some progress early after initial response. To investigate this, we examined 44 whole-genome and whole-exome sequencing data samples from 34 patients with t(11;14) MM treated with venetoclax. The presence of mutations in the RAS pathway was strongly associated with shortened progression-free survival (PFS) and was validated in an independent cohort of 21 patients with MM. The presence of 1q gain was also associated with shorter PFS in patients without RAS mutations. In 10 patients with paired prevenetoclax and postvenetoclax treatment samples, postvenetoclax progression was recurrently driven by the selection of genomic events in the BCL2/MCL1 and RAS pathways and of high-risk features (eg, loss of TP53 and CDKN2C). Overall, our study shows that comprehensive genomic profiling can identify most mechanisms underlying resistance to BCL2 inhibition in t(11;14)(CCND1;IGH) MM.

Abstract Transcription factors and their cofactors are major and selective nononcogene dependencies in multiple myeloma cells. By performing a gain-of-function perturbation screen in human multiple myeloma cell lines, we identified the inhibitor of DNA binding (ID) genes as putative suppressors of multiple myeloma cell fitness. Among them, ID2 was found to be downregulated in multiple myeloma patient cells and acted as a tumor suppressor by directly binding and repressing the basic helix-loop-helix factor TCF3, also known as E2A. Lower ID2 expression in multiple myeloma cells conferred a proliferative advantage by increasing TCF3 activity, leading to a dependency on this transcription factor. In contrast, ID2 overexpression reduced TCF3 binding to DNA, which resulted in cell-cycle arrest and a halt in multiple myeloma cell proliferation. The myeloma bone marrow milieu supported this process by further decreasing the expression of ID2 and enhancing TCF3 activity, partly via IL6, revealing a mechanism by which the tumor microenvironment affects multiple myeloma cell behavior. Significance: Multiple myeloma cells exploit the oncogenic and proliferative potential of TCF3 by downregulating the transcriptional regulator ID2, a process facilitated by the bone marrow microenvironment.

Introduction: Multiple myeloma (MM) with t(11;14)(q13;q32), resulting in CCND1;IGH rearrangement, remains the only subset amenable to targeted therapy with the BCL2 inhibitor, Venetoclax. Nevertheless, not all patients with this translocation respond to treatment and some progress early after initial response. The mechanisms behind this clinical heterogeneity are largely unknown. Methods: To investigate mechanisms of heterogeneity in response to Venetoclax, we interrogated 44 whole genome sequencing (WGS) and 24 whole exome sequencing data from a training cohort including 33 patients with t(11;14) MM treated with Venetoclax at the Mayo Clinic. In a validation cohort, we extracted single nucleotide variants (SNV) and indels from 17 MM patients with t(11;14) with available RNAseq data treated with Venetoclax at Emory University. Results: In the training cohort the median number of lines of therapy was 4 and Venetoclax drug regimens included combinations with dexamethasone, proteasome inhibitors, anti-CD38 agents, and immunomodulatory agents. After a median follow up of 30 months, the median progression free survival (PFS) was 9 months, and 13 patients (39%) achieved a complete remission (CR) to therapy, while 9 (27%) were refractory (i.e. no response or progressed within the first 60 days). For our driver discovery analysis, we integrated SNV, indels, copy number variants (CNV) and structural variants (SV). Deletion 17p (del17p) or TP53 mutation, MCL1 gain/amp, and hyper-APOBEC did not correlate with PFS in our cohort (p=0.59, p=0.71, and p=0.05, respectively). In contrast, mutations involving the MAPK pathway (NRAS, KRAS, BRAF, FGFR3) displayed the strongest association with short PFS, which was retained after adjusting for del17p/TP53 mutations, and 1q gain/amp (HR: 4.76, CI: 1.8-12.9, p=0.002). The association between suboptimal response and resistance to Venetoclax in MM patients with MAPK pathway mutations (n=7) was also observed in our validation cohort from Emory University (p=0.046). This is in line with what has been reported in functional studies in acute myeloid leukemia, where RAS pathway promotes MCL1 under the selective pressure of Venetoclax (Zhang et al. Signal Transduct Target Ther 2022). While MAPK mutations explained most cases of refractory patients (86%), the mechanisms through which initially responding patients subsequently progressed remain unclear. To investigate this, we interrogated 10 patients with pre- and post-treatment samples collected at disease progression (PD). We identified genomic events acquired/selected at PD of genes involved in the anti-apoptotic MCL1/BCL2 and MAPK pathways in 6/10 cases. One patient was refractory to Venetoclax and had double minutes involving BCL2, which caused multiple focal amplifications (>30 copies), likely reflecting extrachromosomal DNA. Another patient progressed with an acquired, focal deletion on BCL2L11. In two patients, the progressing clones harbored mutations in MAPK pathways (one BRAF and the other WNK2). Two patients who relapsed after sustained CR had major branching evolution where a previously undetectable focal deletion of PMAIP1 (NOXA) was seen at PD. Importantly, one of these deletions was a biallelic event. Venetoclax sensitivity assays conducted on PMAIP1 knockout in the MM t(11;14) cell line, KMS12PE, demonstrated significantly reduced cell death across all drug concentrations (p<0.05), validating the loss of NOXA as a mechanism of resistance to Venetoclax in MM. Finally, similarly with what was previously reported in chronic lymphocyte leukemia treated with Venetoclax (e.g., Herling et al Nat Comm 2018), CDKN2A/B monoallelic and biallelic focal deletions were selected post-Venetoclax, at PD in 4 patients while TP53 SNV/deletions were selected from subclonal to clonal in 2 patients at PD. Conclusion: In this study, we demonstrate that mutations involving the MAPK pathway identify a population of MM patients with t(11;14)(CCND1;IGH) who are unlikely to respond to Venetoclax. We also show that genomic events involving the MAPK and BCL2/MCL1 pathways are primarily responsible for disease progression after initial response. These findings highlight the critical role of comprehensive genomic profiling through WGS in identifying mechanisms of resistance to targeted therapies in MM, aiding patient selection for Venetoclax, and guiding the development of novel drug combinations.

Background: Multiple myeloma (MM) is an incurable plasma cell malignancy with several recurrent primary cytogenetic abnormalities including t(11,14)(q13;q32), resulting in increased expression of CCND1. Although 40% of patients with t(11;14) MM respond favorably to the BCL-2 inhibitor venetoclax (ven), 60% fail to achieve an objective response. Previous studies have suggested upregulation of PMAIP1, encoding the pro-apoptotic NOXA protein, in association with ven sensitivity (Leblay, Blood, 2024). Here we sought to evaluate the role for NOXA in ven response in t(11;14) MM. Methods: We analyzed whole exome, whole genome and RNA sequencing data of CD138+ tumor cells from 1,359 patients with MM in the CoMMpass and Mayo Clinic cohorts. The histone mark H3K27ac was assessed in t(11;14)+ KMS12PE and U266B1 human myeloma cell lines using ChIP-seq. Gene expression and protein levels of PMAIP1 (NOXA) were also detected by qRT-PCR and Western blot, respectively. Live cell detection was evaluated by Annexin V and PI using flow cytometry and by YoYo3 dye using IncuCyte following treatment with ven. Results: The t(11;14) MM subtype was associated with a 1.6-5.5-fold increased expression of PMAIP1 compared to non-t(11,14) MM in the CoMMpass and Mayo Clinic cohorts (p<0.0001). These findings are consistent with a significant correlation between CCND1 and PMAIP1 expression in both cohorts (p<0.001). To evaluate whether cyclin D1 is directly associated with increasing expression of PMAIP1 in t(11;14) MM, we generated a cyclin D1 knock out (KO) version of the U266B1 line using CRISPR/cas9. Two independent U266B1 KO cell lines showed reduced expression of PMAIP1 and NOXA protein compared to the WT U266B1 consistent with a role for cyclin D1 in regulating PMAIP1 expression. In addition, the KO lines had a 4.8-fold reduced H3K27ac signal at the PMAIP1 locus in comparison to the WT control suggesting a role for cyclin D1 in promoting H3K27ac. We previously observed increased expression of CCND1 in association with the G allele compared to the A allele of rs9344 in both the CoMMpass and Mayo Clinic cohorts (Leukemia, in press). We next asked whether the G risk allele of rs9344 within CCND1 impacted gene expression and H3k27ac levels at the PMAIP1 locus. To study the allele specific effect of CCND1, we used the CRISPR/cas9 edited KMS12PE WT (rs9344 AA) and Knock In (KI) (rs9344 GG) cell lines. In the KI cells, the G allele was accompanied by an ~2-fold increased expression of PMAIP1 and NOXA protein detected by western blot (<0.05). Increased PMAIP1 expression in KI lines was also accompanied by a 2.2-fold increased H3K27ac further suggesting a role for the G allele of rs9344 in promoting H3K27ac and gene expression at the PMAIP1 locus. Finally, we evaluated the consequence of PMAIP1 expression levels in relation to ven response. Increased expression of PMAIP1 and elevated levels of NOXA protein in the KMS12PE KI lines was associated with increased sensitivity to ven resulting in ~50% increase in cell death of KI vs. WT lines (p<0.02). Conclusion: We demonstrate a role for cyclin D1 in regulating PMAIP1 expression through H3K27ac resulting in increased NOXA protein. Levels of NOXA are associated with the G allele of rs9344 and increased venetoclax response.

Abstract 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

Abstract 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

Abstract Dendritic cells (DC) are essential for the long term survival of plasma cells (PC) in the bone marrow and for the maintenance of persistent antibody titers. But the molecular basis for these interactions is not very clear. Both normal and malignant PC (like myeloma) express the T-cell co-stimulatory molecule CD28. Our data show that CD28 expressed on PC play a major role in CD80/CD86 (B7) mediated DC induction of growth factors such as IL-6 and immunosuppressive enzyme indoleamine 2,3 dioxygenase (IDO) as well as in antibody or DC mediated protection of PC in vitro. While IL-6 protected PC in vitro, IDO induced the formation of regulatory T cells in culture. Use of myeloma drug lenalidomide to inhibit DC-B7 expression or anti-CD28(Fab)2 fragments to block CD28-B7 interactions reduced DC ability to protect MM cells or produce IL-6/IDO. Moreover, just like in T-cell models, CD28 activation in malignant PC triggers similar molecular events such as tyrosine phosphorylation and association of p85. We believe that downstream of CD28, while PI3K/Akt governs survival, the Vav1/SLP76 pathway regulate PC proliferation and these form the molecular basis for a bi-directional molecular bridge between the PC and bone marrow microenvironment. Targeting this bridge could form the basis for future therapeutic strategies. Our data suggests a role for CD28-B7 in plasma cell survival by modulating the bone marrow microenvironment.

Abstract We have previously shown (AAI2009) that dendritic cells (DC) contribute to the survival of long lived plasma cells in the bone marrow, which are important to maintain persistent antibody titres. What is not clear are the molecular interactions that are involved. Both normal and malignant plasma cells(like myeloma cells) express the T-cell co-stimulatory molecule CD28, and both interact with DCs in the bone marrow microenvironment. CD28 activation on MM triggers similar molecular partners downstream as observed in T-cells. Gene expression analysis of normal and malignant plasma cells shows that CD28 expression correlates with poor prognostic subgroups that are associated with specific genetic lesions. In our data, while activation of CD28 with antibodies or DC-CD80/CD86 (B7) protects myeloma cells against cell death, the ligation of B7 on DCs induces DC production of plasma cell survival factor IL-6, and immunosuppressive factor indoleamine 2,3-dioxygenase (IDO). Use of Lenalidomide (a common myeloma drug) to inhibit DC-B7 expression or anti-CD28(Fab)2 fragments to block CD28-B7 interactions reduced DC ability to protect MM cells. Experiments with αCD28(Fab)2 or inhibitors of notch pathway suggest a cross talk that regulates IL-6 production. IDO activity in DC-MM cocultures, suppressed T-cell proliferation in MLRs while an IDO inhibitor alleviated this. Our data suggests a role for CD28-B7 in plasma cell survival by modulating the bone marrow microenvironment.