NSD2 drives t(4;14) myeloma cell dependence on adenylate kinase 2 by diverting one-carbon metabolism to the epigenome.

Chromosomal translocation (4;14), an adverse prognostic factor in multiple myeloma (MM), drives overexpression of the histone methyltransferase NSD2. A genome-wide CRISPR screen in MM cells identified adenylate kinase 2 (AK2), an enzyme critical for high energy phosphate transfer from the mitochondria, as an NSD2-driven vulnerability. AK2 suppression in t(4;14) MM cells decreased NADP(H) critical for conversion of ribonucleotides to deoxyribonucleosides, leading to replication stress, DNA damage and apoptosis. Driving a large genome-wide increase in chromatin methylation, NSD2 overexpression depletes S-adenosylmethionine (SAM), compromising synthesis of creatine from its precursor guanidinoacetate. Creatine supplementation restored NADP(H) levels, reduced DNA damage and rescued AK2-deficient t(4;14) MM cells. As the creatine phosphate shuttle constitutes an alternative means for mitochondrial high energy phosphate transport, these results indicate that NSD2-driven creatine depletion underlies the hypersensitivity of t(4;14) MM cells to AK2 loss. Furthermore, AK2 depletion in t(4;14) cells impaired protein folding in the endoplasmic reticulum consistent with impaired utilization of mitochondrial ATP. Accordingly, AK2 suppression increased sensitivity of MM cells to proteasome inhibition. These findings delineate a novel mechanism in which aberrant transfer of carbon to the epigenome creates a metabolic vulnerability, with direct therapeutic implications for t(4;14) MM.

Targeting high-risk multiple myeloma genotypes with optimized anti-CD70 CAR-T cells.

Despite the success of BCMA-targeting CAR-Ts in multiple myeloma, patients with high-risk cytogenetic features still relapse most quickly and are in urgent need of additional therapeutic options. Here, we identify CD70, widely recognized as a favorable immunotherapy target in other cancers, as a specifically upregulated cell surface antigen in high risk myeloma tumors. We use a structure-guided design to define a CD27-based anti-CD70 CAR-T design that outperforms all tested scFv-based CARs, leading to >80-fold improved CAR-T expansion in vivo. Epigenetic analysis via machine learning predicts key transcription factors and transcriptional networks driving CD70 upregulation in high risk myeloma. Dual-targeting CAR-Ts against either CD70 or BCMA demonstrate a potential strategy to avoid antigen escape-mediated resistance. Together, these findings support the promise of targeting CD70 with optimized CAR-Ts in myeloma as well as future clinical translation of this approach.

KDM6A Regulates Immune Response Genes in Multiple Myeloma.

The histone H3K27 demethylase KDM6A is a tumor suppressor in multiple cancers, including multiple myeloma (MM). We created isogenic MM cells disrupted for KDM6A and tagged the endogenous protein to facilitate genome wide studies. KDM6A binds genes associated with immune recognition and cytokine signaling. Most importantly, KDM6A binds and activates NLRC5 and CIITA encoding regulators of Major Histocompatibility Complex (MHC) genes. Patient data indicate that NLRC5 and CIITA, are downregulated in MM with low KDM6A expression. Chromatin analysis shows that KDM6A binds poised and active enhancers and KDM6A loss led to decreased H3K27ac at enhancers, increased H3K27me3 levels in body of genes bound by KDM6A and decreased gene expression. Reestablishing histone acetylation with an HDAC3 inhibitor leads to upregulation of MHC expression, offering a strategy to restore immunogenicity of KDM6A deficient tumors. Loss of Kdm6a in murine RAS-transformed fibroblasts led to increased growth in vivo associated with decreased T cell infiltration. Statement of significance: We show that KDM6A participates in immune recognition of myeloma tumor cells by directly regulating the expression of the master regulators of MHC-I and II, NLRC5 and CIITA. The expression of these regulators can by rescued by the HDAC3 inhibitors in KDM6A-null cell lines.

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.

Integrated epigenetic and transcriptional single-cell analysis of t(11;14) multiple myeloma and its BCL2 dependency.

ABSTRACT: The translocation t(11;14) occurs in 20% of patients with multiple myeloma (MM) and results in the upregulation of CCND1. Nearly two-thirds of t(11;14) MM cells are BCL2 primed and highly responsive to the oral BCL2 inhibitor venetoclax. Although it is evident that this unique sensitivity to venetoclax depends on the Bcl-2 homology domain 3- proapoptotic protein priming of BCL2, the biology underlying t(11;14) MM dependency on BCL2 is poorly defined. Importantly, the epigenetic regulation of t(11;14) transcriptomes and its impact on gene regulation and clinical response to venetoclax remain elusive. In this study, by integrating assay for transposase-accessible chromatin by sequencing (ATAC-seq) and RNA-seq at the single-cell level in primary MM samples, we have defined the epigenetic regulome and transcriptome associated with t(11;14) MM. A B-cell-like epigenetic signature was enriched in t(11;14) MM, confirming its phylogeny link to B-cell rather than plasma cell biology. Of note, a loss of a B-cell-like epigenetic signature with a gain of canonical plasma cell transcription factors was observed at the time of resistance to venetoclax. In addition, MCL1 and BCL2L1 copy number gains and structural rearrangements were linked to venetoclax resistance in patients with t(11;14) MM. To date, this is the first study in which both single-cell (sc) ATAC-seq and scRNA-seq analysis are integrated into primary MM cells to obtain a deeper resolution of the epigenetic regulome and transcriptome associated with t(11;14) MM biology and venetoclax resistance.

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.

Integrated Genomic Analysis of Primary Prostate Tumor Foci and Corresponding Lymph Node Metastases Identifies Mutations and Pathways Associated with Metastasis.

Prostate cancer is a highly heterogeneous disease and mortality is mainly due to metastases but the initial steps of metastasis have not been well characterized. We have performed integrative whole exome sequencing and transcriptome analysis of primary prostate tumor foci and corresponding lymph node metastases (LNM) from 43 patients enrolled in clinical trial. We present evidence that, while there are some cases of clonally independent primary tumor foci, 87% of primary tumor foci and metastases are descended from a common ancestor. We demonstrate that genes related to oxidative phosphorylation are upregulated in LNM and in African-American patients relative to White patients. We further show that mutations in TP53, FLT4, EYA1, NCOR2, CSMD3, and PCDH15 are enriched in prostate cancer metastases. These findings were validated in a meta-analysis of 3929 primary tumors and 2721 metastases and reveal a pattern of molecular alterations underlying the pathology of metastatic prostate cancer. We show that LNM contain multiple subclones that are already present in primary tumor foci. We observed enrichment of mutations in several genes including understudied genes such as EYA1, CSMD3, FLT4, NCOR2, and PCDH15 and found that mutations in EYA1 and CSMD3 are associated with a poor outcome in prostate cancer.

Genome-scale functional genomics identify genes preferentially essential for multiple myeloma cells compared to other neoplasias.

Clinical progress in multiple myeloma (MM), an incurable plasma cell (PC) neoplasia, has been driven by therapies that have limited applications beyond MM/PC neoplasias and do not target specific oncogenic mutations in MM. Instead, these agents target pathways critical for PC biology yet largely dispensable for malignant or normal cells of most other lineages. Here we systematically characterized the lineage-preferential molecular dependencies of MM through genome-scale clustered regularly interspaced short palindromic repeats (CRISPR) studies in 19 MM versus hundreds of non-MM lines and identified 116 genes whose disruption more significantly affects MM cell fitness compared with other malignancies. These genes, some known, others not previously linked to MM, encode transcription factors, chromatin modifiers, endoplasmic reticulum components, metabolic regulators or signaling molecules. Most of these genes are not among the top amplified, overexpressed or mutated in MM. Functional genomics approaches thus define new therapeutic targets in MM not readily identifiable by standard genomic, transcriptional or epigenetic profiling analyses.

CD8+ T cell activation in cancer comprises an initial activation phase in lymph nodes followed by effector differentiation within the tumor.

Improvements in tumor immunotherapies depend on better understanding of the anti-tumor T cell response. By studying human tumor-draining lymph nodes (TDLNs), we found that activated CD8+ T cells in TDLNs shared functional, transcriptional, and epigenetic traits with TCF1+ stem-like cells in the tumor. The phenotype and TCR overlap suggested that these TDLN cells were precursors to tumor-resident stem-like CD8+ T cells. Murine tumor models revealed that tumor-specific CD8+ T cells were activated in TDLNs but lacked an effector phenotype. These stem-like cells migrated into the tumor, where additional co-stimulation from antigen-presenting cells drove effector differentiation. This model of CD8+ T cell activation in response to cancer is different from that of canonical CD8+ T cell activation to acute viruses, and it proposes two stages of tumor-specific CD8+ T cell activation: initial activation in TDLNs and subsequent effector program acquisition within the tumor after additional co-stimulation.

Inhibition of NADPH Oxidase-ROS Signal using Hyaluronic Acid Nanoparticles for Overcoming Radioresistance in Cancer Therapy.

Upregulation of NADPH oxidases (NOXs) in cancer cells leads to chronic increase in intracellular reactive oxygen species (ROS) and adaptation to a high ROS level for cell survival and, thereby, low sensitivity to radiotherapy. To overcome resistance to radiotherapy, we have developed a bioactive and CD44 targeted hyaluronic acid nanoparticle encapsulated with an NOX inhibitor, GKT831 (HANP/GKT831). We found that HANP/GKT831 had stronger inhibitory effects on ROS generation and cell proliferation than that of GKT831 alone in cancer cells. Systemic delivery of HANP/GKT831 led to the targeted accumulation in breast cancer patient derived xenograft (PDX) tumors in nude mice. Importantly, the combination of systemic delivery of HANP/GKT831 with a low dose of local radiotherapy significantly enhanced tumor growth inhibition in breast cancer PDX models. Our results showed that HANP/GKT831 primed tumor cells to radiation-induced DNA damage and cell death by downregulation of DNA repair function and oncogenic signal pathways.

Therapeutic implications of mitochondrial stress-induced proteasome inhibitor resistance in multiple myeloma.

The connections between metabolic state and therapy resistance in multiple myeloma (MM) are poorly understood. We previously reported that electron transport chain (ETC) suppression promotes sensitivity to the BCL-2 antagonist venetoclax. Here, we show that ETC suppression promotes resistance to proteasome inhibitors (PIs). Interrogation of ETC-suppressed MM reveals integrated stress response-dependent suppression of protein translation and ubiquitination, leading to PI resistance. ETC and protein translation gene expression signatures from the CoMMpass trial are down-regulated in patients with poor outcome and relapse, corroborating our in vitro findings. ETC-suppressed MM exhibits up-regulation of the cystine-glutamate antiporter SLC7A11, and analysis of patient single-cell RNA-seq shows that clusters with low ETC gene expression correlate with higher SLC7A11 expression. Furthermore, erastin or venetoclax treatment diminishes mitochondrial stress-induced PI resistance. In sum, our work demonstrates that mitochondrial stress promotes PI resistance and underscores the need for implementing combinatorial regimens in MM cognizant of mitochondrial metabolic state.

Multivariant Transcriptome Analysis Identifies Modules and Hub Genes Associated with Poor Outcomes in Newly Diagnosed Multiple Myeloma Patients.

The molecular mechanisms underlying chemoresistance in some newly diagnosed multiple myeloma (MM) patients receiving standard therapies (lenalidomide, bortezomib, and dexamethasone) are poorly understood. Identifying clinically relevant gene networks associated with death due to MM may uncover novel mechanisms, drug targets, and prognostic biomarkers to improve the treatment of the disease. This study used data from the MMRF CoMMpass RNA-seq dataset (N = 270) for weighted gene co-expression network analysis (WGCNA), which identified 21 modules of co-expressed genes. Genes differentially expressed in patients with poor outcomes were assessed using two independent sample t-tests (dead and alive MM patients). The clinical performance of biomarker candidates was evaluated using overall survival via a log-rank Kaplan-Meier and ROC test. Four distinct modules (M10, M13, M15, and M20) were significantly correlated with MM vital status and differentially expressed between the dead (poor outcomes) and the alive MM patients within two years. The biological functions of modules positively correlated with death (M10, M13, and M20) were G-protein coupled receptor protein, cell-cell adhesion, cell cycle regulation genes, and cellular membrane fusion genes. In contrast, a negatively correlated module to MM mortality (M15) was the regulation of B-cell activation and lymphocyte differentiation. MM biomarkers CTAG2, MAGEA6, CCND2, NEK2, and E2F2 were co-expressed in positively correlated modules to MM vital status, which was associated with MM's lower overall survival.

PDZ Proteins SCRIB and DLG1 Regulate Myeloma Cell Surface CD86 Expression, Growth, and Survival.

Despite advances in the treatment of multiple myeloma in the past decades, the disease remains incurable, and understanding signals and molecules that can control myeloma growth and survival are important for the development of novel therapeutic strategies. One such molecule, CD86, regulates multiple myeloma cell survival via its interaction with CD28 and signaling through its cytoplasmic tail. Although the CD86 cytoplasmic tail has been shown to be involved in drug resistance and can induce molecular changes in multiple myeloma cells, its function has been largely unexplored. Here, we show that CD86 cytoplasmic tail has a role in trafficking CD86 to the cell surface. This is due in part to a PDZ-binding motif at its C-terminus which is important for proper trafficking from the Golgi apparatus. BioID analysis revealed 10 PDZ domain-containing proteins proximal to CD86 cytoplasmic tail in myeloma cells. Among them, we found the planar cell polarity proteins, SCRIB and DLG1, are important for proper CD86 surface expression and the growth and survival of myeloma cells. These findings indicate a mechanism by which myeloma cells confer cellular survival and drug resistance and indicate a possible motif to target for therapeutic gain. IMPLICATIONS: These findings demonstrate the importance of proper trafficking of CD86 to the cell surface in myeloma cell survival and may provide a new therapeutic target in this disease.

Aberrant Extrafollicular B Cells, Immune Dysfunction, Myeloid Inflammation, and MyD88-Mutant Progenitors Precede Waldenstrom Macroglobulinemia.

Waldenstrom macroglobulinemia (WM) and its precursor IgM gammopathy are distinct disorders characterized by clonal mature IgM-expressing B-cell outgrowth in the bone marrow. Here, we show by high-dimensional single-cell immunogenomic profiling of patient samples that these disorders originate in the setting of global B-cell compartment alterations, characterized by expansion of genomically aberrant extrafollicular B cells of the nonmalignant clonotype. Alterations in the immune microenvironment preceding malignant clonal expansion include myeloid inflammation and naïve B- and T-cell depletion. Host response to these early lesions involves clone-specific T-cell immunity that may include MYD88 mutation-specific responses. Hematopoietic progenitors carry the oncogenic MYD88 mutations characteristic of the malignant WM clone. These data support a model for WM pathogenesis wherein oncogenic alterations and signaling in progenitors, myeloid inflammation, and global alterations in extrafollicular B cells create the milieu promoting extranodal pattern of growth in differentiated malignant cells. SIGNIFICANCE: These data provide evidence that growth of the malignant clone in WM is preceded by expansion of extrafollicular B cells, myeloid inflammation, and immune dysfunction in the preneoplastic phase. These changes may be related in part to MYD88 oncogenic signaling in pre-B progenitor cells and suggest a novel model for WM pathogenesis. This article is highlighted in the In This Issue feature, p. 549.

Ten-eleven translocation protein 1 modulates medulloblastoma progression.

BACKGROUND: Medulloblastoma (MB) is the most common malignant pediatric brain tumor that originates in the cerebellum and brainstem. Frequent somatic mutations and deregulated expression of epigenetic regulators in MB highlight the substantial role of epigenetic alterations. 5-hydroxymethylcytosine (5hmC) is a highly abundant cytosine modification in the developing cerebellum and is regulated by ten-eleven translocation (TET) enzymes. RESULTS: We investigate the alterations of 5hmC and TET enzymes in MB and their significance to cerebellar cancer formation. We show total abundance of 5hmC is reduced in MB, but identify significant enrichment of MB-specific 5hmC marks at regulatory regions of genes implicated in stem-like properties and Nanog-binding motifs. While TET1 and TET2 levels are high in MBs, only knockout of Tet1 in the smoothened (SmoA1) mouse model attenuates uncontrolled proliferation, leading to a favorable prognosis. The pharmacological Tet1 inhibition reduces cell viability and platelet-derived growth factor signaling pathway-associated genes. CONCLUSIONS: These results together suggest a potential key role of 5hmC and indicate an oncogenic nature for TET1 in MB tumorigenesis, suggesting it as a potential therapeutic target for MBs.

Chromatin Accessibility Identifies Regulatory Elements Predictive of Gene Expression and Disease Outcome in Multiple Myeloma.

PURPOSE: Multiple myeloma is a malignancy of plasma cells. Extensive genetic and transcriptional characterization of myeloma has identified subtypes with prognostic and therapeutic implications. In contrast, relatively little is known about the myeloma epigenome. EXPERIMENTAL DESIGN: CD138+CD38+ myeloma cells were isolated from fresh bone marrow aspirate or the same aspirate after freezing for 1-6 months. Gene expression and chromatin accessibility were compared between fresh and frozen samples by RNA sequencing (RNA-seq) and assay for transpose accessible chromatin sequencing (ATAC-seq). Chromatin accessible regions were used to identify regulatory RNA expression in more than 700 samples from newly diagnosed patients in the Multiple Myeloma Research Foundation CoMMpass trial (NCT01454297). RESULTS: Gene expression and chromatin accessibility of cryopreserved myeloma recapitulated that of freshly isolated samples. ATAC-seq performed on a series of biobanked specimens identified thousands of chromatin accessible regions with hundreds being highly coordinated with gene expression. More than 4,700 of these chromatin accessible regions were transcribed in newly diagnosed myelomas from the CoMMpass trial. Regulatory element activity alone recapitulated myeloma gene expression subtypes, and in particular myeloma subtypes with immunoglobulin heavy chain translocations were defined by transcription of distal regulatory elements. Moreover, enhancer activity predicted oncogene expression implicating gene regulatory mechanisms in aggressive myeloma. CONCLUSIONS: These data demonstrate the feasibility of using biobanked specimens for retrospective studies of the myeloma epigenome and illustrate the unique enhancer landscapes of myeloma subtypes that are coupled to gene expression and disease progression.

Oncolytic herpes simplex virus infects myeloma cells in vitro and in vivo.

Because most patients with multiple myeloma (MM) develop resistance to current regimens, novel approaches are needed. Genetically modified, replication-competent oncolytic viruses exhibit high tropism for tumor cells regardless of cancer stage and prior treatment. Receptors of oncolytic herpes simplex virus 1 (oHSV-1), NECTIN-1, and HVEM are expressed on MM cells, prompting us to investigate the use of oHSV-1 against MM. Using oHSV-1-expressing GFP, we found a dose-dependent increase in the GFP+ signal in MM cell lines and primary MM cells. Whereas NECTIN-1 expression is variable among MM cells, we discovered that HVEM is ubiquitously and highly expressed on all samples tested. Expression of HVEM was consistently higher on CD138+/CD38+ plasma cells than in non-plasma cells. HVEM blocking demonstrated the requirement of this receptor for infection. However, we observed that, although oHSV-1 could efficiently infect and kill all MM cell lines tested, no viral replication occurred. Instead, we identified that oHSV-1 induced MM cell apoptosis via caspase-3 cleavage. We further noted that oHSV-1 yielded a significant decrease in tumor volume in two mouse xenograft models. Therefore, oHSV-1 warrants exploration as a novel potentially effective treatment option in MM, and HVEM should be investigated as a possible therapeutic target.

Venetoclax sensitivity in multiple myeloma is associated with B-cell gene expression.

Venetoclax is a highly potent, selective BCL2 inhibitor capable of inducing apoptosis in cells dependent on BCL2 for survival. Most myeloma is MCL1-dependent; however, a subset of myeloma enriched for translocation t(11;14) is codependent on BCL2 and thus sensitive to venetoclax. The biology underlying this heterogeneity remains poorly understood. We show that knockdown of cyclin D1 does not induce resistance to venetoclax, arguing against a direct role for cyclin D1 in venetoclax sensitivity. To identify other factors contributing to venetoclax response, we studied a panel of 31 myeloma cell lines and 25 patient samples tested for venetoclax sensitivity. In cell lines, we corroborated our previous observation that BIM binding to BCL2 correlates with venetoclax response and further showed that knockout of BIM results in decreased venetoclax sensitivity. RNA-sequencing analysis identified expression of B-cell genes as enriched in venetoclax-sensitive myeloma, although no single gene consistently delineated sensitive and resistant cells. However, a panel of cell surface makers correlated well with ex vivo prediction of venetoclax response in 21 patient samples and may serve as a biomarker independent of t(11;14). Assay for transposase-accessible chromatin sequencing of myeloma cell lines also identified an epigenetic program in venetoclax-sensitive cells that was more similar to B cells than that of venetoclax-resistant cells, as well as enrichment for basic leucine zipper domain-binding motifs such as BATF. Together, these data indicate that remnants of B-cell biology are associated with BCL2 dependency and point to novel biomarkers of venetoclax-sensitive myeloma independent of t(11;14).