Dear Editor,
Multiple myeloma (MM) is a rarely curable plasma cell malignancy of the bone marrow (BM) that provides a permissive tumor microenvironment (TME), supporting tumor cell growth and dissemination and conferring therapy resistance. The TME contains cellular components, specifically stromal cells, osteoclasts, osteoblasts, osteocytes, endothelial and immune cells, and a non-cellular component, the extracellular matrix (ECM). In cancer, the ECM is an important determiner of cell fate and composition of the TME. Recent research has coined the term “matrisome” for the ensemble of genes encoding ECM proteins and ECM-associated proteins and defined gene sets for core matrisome (approximately 274 genes) and matrisome-associated genes (approximately 1027 genes), including secreted modifiers [1]. Specified cancer matrisomes regulate proliferation, migration, and survival [1]. Hence, changes in ECM composition, integrity, abundance, biomechanical properties, and related signal transduction contribute to tumor progression and outcome in patients [1]. In MM, the ECM bidirectionally interacts with MM cells and co-inhabitants of tumor cell/metastatic niches. Notably, expression levels of genes in MM cells involved in the interaction with the TME have been linked to better (BMP6 [2]) or worse (ANXA2, LGALS1 [3, 4]) survival of patients. We and others have shown that NOTCH signaling alters the TME through juxtacrine signaling between signal-sending cells such as MM or stromal cells expressing the ligands, and signal-receiving cells expressing the receptors [5, 6]. Whether deregulated NOTCH signaling in MM cells controls the expression of genes that dysregulate ECM composition in the BM niche and have prognostic significance, is unknown. Here, we correlated transcriptome profiles of NOTCH-depleted MM cells with recently published matrisome libraries, to identify NOTCH-regulated genes that belong to the matrisome and are related to patient survival.
We transduced human RPMI 8226 and MM.1S cells with shRNAs for knockdown of NOTCH(N)1 and N2 receptors. Efficiency and specificity of depletion were validated by qPCR, flow cytometry analysis, and immunoblotting (Supplementary Fig. 1). RPMI 8226 and MM.1S cells showed different levels of the intracellular cleaved domain of N2 (N2IC), indicating a variable strength of N2 activation (Supplementary Fig. 2). In addition, N1 and N2 depleted MM cells were less viable and more sensitive to bortezomib, melphalan, and lenalidomide (Supplementary Fig. 1), confirming that these receptors control growth and drug resistance. High-throughput transcriptome profiling revealed decreased NOTCH target gene expression of HES4, HES7 in RPMI 8226, and HEY2, HEYL in MM.1S cells (Supplementary Table 1). Many of the 19,720 analyzed genes were significantly up- or downregulated by at least one of the shRNAs in RPMI 8226 cells (shN1: 2761 up and 2758 down; shN2: 3028 up and 3355 down; cut-off: padj < 0.01), whereas in MM.1S cells less genes were regulated (shN1: 503 up and 1032 down; shN2: 823 up and 1641 down; cut-off: padj < 0.01, Supplementary Table 1, Supplementary Fig. 3). Among the top 20 genes commonly regulated after N1 and N2 depletion, we identified nine matrisome-associated genes in RPMI 8226 cells: (i) down - CXCL9, CXCL10, CCL8, MMP13, TNFSF13B, TNFSF10, and (ii) up - LEFTY2, SERPINE1, ZP1 (Fig. 1). In MM.1S cells, two out of three commonly upregulated genes (CLEC7A, TGFA) encode matrisome-associated proteins (Supplementary Fig. 3).
These findings prompted the systematic search for NOTCH-driven matrisome genes within the entire gene expression data set. To this end, we used the MatrisomeDB database that provides live cross-referencing to gene and protein databases for every ECM and ECM-associated gene [1]. GO analysis revealed that N1 and N2 regulate both core matrisome genes and a series of matrisome-associated genes in RPMI 8226 (Table 1) and MM.1S cells (Supplementary Table 2). Overall comparison between up- or downregulated genes showed that expression of 14 and 34 matrisome genes is commonly regulated by N1 and N2 in both cell lines (Supplementary Table 1). QPCR analysis or immunoblotting demonstrated that lower levels of N1 or N2 correlate with lower levels of HPSE in NOTCH-depleted RPMI 8226 and MM.1S cells (Supplementary Fig. 2). The same trend of low HPSE expression can be found in the MM cell line AMO-1, which is characterized by low NOTCH levels. Similarly, lower levels of MMP13, S100A6, IGF1 correspond to decreased NOTCH levels in RPMI 8226 cells, and low levels of MMP13 to low NOTCH levels in AMO-1 cells. In contrast, lower levels of N1 and N2 are associated with higher protein levels of EMID1, TGFBI, and C1QA in RPMI 8226 cells (Supplementary Fig. 2). These data confirm the reliability of the transcriptome analysis and the regulation of matrisome genes through N1 and N2 in MM cells.
In parallel, we performed gene set enrichment analysis (GSEA) to determine expression changes in gene sets after N1 and N2 depletion. Both receptors activate immune system-associated and cytokine activity signatures such as cytokine receptor binding and inflammatory response in RPMI 8226 cells (Fig. 1) or leukocyte cell-cell adhesion, toll-like receptor signaling pathway, or chemoattractant activity in MM.1S cells (Supplementary Table 3), confirming that NOTCH controls a cytokine network, defining a supportive TME in MM [6]. Interestingly, N1 regulates genes associated with bone remodeling and resorption in RPMI 8226 (GO:0046850; GO:0045124), and N1 and N2 control genes associated with osteoclast differentiation in MM.1S (GO:0045670) cells including RUNX2 or SPP1 (Fig. 1; Supplementary Table 3). RUNX2 and SPP1 (osteopontin) control bone homeostasis in skeletal precursors. In MM, RUNX2 may similarly control ECM-modifying genes such as MMP13 and SPP1, and RUNX2 expression correlated with severity of the disease [7]. High levels of MM cell-derived MMP13 enhance the osteolytic activity of osteoclasts and correlate with bone lesions in MM patients [8]. Of note, SPP1 is upregulated as part of a prognostic cancer core matrisome signature identified by transcriptomics and proteomics in breast and colon cancer [1].
To correlate expression levels of ECM genes and patient survival, we first analyzed gene expression in samples of BM plasma cells (BMPC) from healthy donors, in patient samples from monoclonal gammopathy of undetermined significance (MGUS), smoldering (sMM), untreated MM, and in samples of human MM cell lines (HMCL) [9]. Next, we determined their association with survival in MM patients (Fig. 1). We focused the analysis on the 64 matrisome genes regulated in RPMI 8226 cells (Table 1), since we found the same classes of ECM glycoproteins, regulators, or secreted factors regulated in MM.1 S cells (Supplementary Table 2). Seven out of the 64 matrisome genes, TGFBI, C1QA, S100A6, IGF1, HPSE, CXCL12, and CXCL8, showed an association with progression-free and overall survival (Fig. 1, Supplementary Fig. 4) in a previously published cohort of MM patients (n = 387) [9]. TGFBI is an N2-driven target gene with low expression being associated with adverse overall survival (Fig. 1). Accordingly, a global DNA hypermethylation analysis linked the methylation status of TGFBI to an unfavorable prognosis [10]. We further identified C1QA as a novel N2-regulated ECM gene. High levels of C1QA were associated with better prognosis of MM patients (Fig. 1). C1QA encodes the A-chain polypeptide of serum complement subcomponent C1q binding to immunoglobulins complexed to antigen and initiating the complement cascade [11]. In skin cutaneous melanoma, C1QA is a novel prognostic biomarker that has a function as a core TME-related gene [11]. Similarly, high levels of NOTCH-driven EMI domain containing 1 (EMID1) correlate with tumor-infiltrating immune cells and are associated with a favorable prognosis in lung adenocarcinoma [12]. S100A6 is a Ca2+-binding protein that belongs to the S100 family controlling cell growth, differentiation, and survival in cancer and cancer stem cells [13]. S100A6 binds to ECM-associated proteins such as LUM, PRELP, IGFBP-1, and high serum levels are positively correlated with cancer progression of gastric, non-small cell lung, ovarian, and urinary bladder cancer [13]. We showed that S100A6 is downregulated after N2 knockdown in MM cells, and high levels were associated with adverse prognosis of MM patients (Fig. 1). Moreover, S100 proteins are classical binding partners of ANXA2, and in pancreatic cancer the interaction between S1006 and ANXA2 promotes motility and invasiveness of cancer cells [14]. In addition, N2 controls the expression of the ECM regulator HPSE that cleaves heparan sulfate glycosaminoglycans from proteoglycan core proteins to small oligosaccharides [15]. HPSE promotes shedding of syndecan-1 from the MM cell surface, modulates the expression of proteases, alters histone acetylation and gene expression, and promotes tumor growth, angiogenesis, and metastasis of MM cells [15].
Taken together, our data show that both NOTCH receptors participate in the transcriptional control of ECM glycoproteins (TGFBI), ECM-affiliated proteins (C1QA), ECM regulators (HPSE) and secreted factors (S100A6, IGF1) in MM cells in vitro, proofing to be of prognostic significance in clinical settings. Our data confirm that the TME and ECM represented by a tumor-associated matrisome contain potential biomarkers and support findings in omental metastasis of ovarian cancer, in which a 22-matrisome gene and protein signature has been identified, predicting overall survival in solid cancers such as breast, head, and neck squamous cell carcinoma, non-small-cell lung adenoma, kidney clear cell carcinoma, hepatocellular carcinoma, colon cancer or pancreatic ductal adenocarcinoma [1]. In MM, similar signatures with prognostic significance should be refined and may confer impact in diagnostic/prognostic classification and the characterization of therapeutic targets as in colorectal cancer [1]. However, further studies are required to comprehensively answer the question how expression changes in NOTCH-driven matrisome-associated proteins in the BM niche promote MM growth and dissemination.
References
Socovich AM, Naba A. The cancer matrisome: from comprehensive characterization to biomarker discovery. Semin Cell Dev Biol. 2019;89:157–66.
Seckinger A, Meissner T, Moreaux J, Goldschmidt H, Fuhler GM, Benner A, et al. Bone morphogenic protein 6: a member of a novel class of prognostic factors expressed by normal and malignant plasma cells inhibiting proliferation and angiogenesis. Oncogene. 2009;28:3866–79.
Glavey SV, Naba A, Manier S, Clauser K, Tahri S, Park J, et al. Proteomic characterization of human multiple myeloma bone marrow extracellular matrix. Leukemia. 2017;31:2426–34.
Seckinger A, Meissner T, Moreaux J, Depeweg D, Hillengass J, Hose K, et al. Clinical and prognostic role of annexin A2 in multiple myeloma. Blood. 2012;120:1087–94.
Schwarzer R, Kaiser M, Acikgoez O, Heider U, Mathas S, Preissner R, et al. Notch inhibition blocks multiple myeloma cell-induced osteoclast activation. Leukemia. 2008;22:2273–7.
Colombo M, Mirandola L, Chiriva-Internati M, Basile A, Locati M, Lesma E, et al. Cancer cells exploit notch signaling to redefine a supportive cytokine milieu. Front Immunol. 2018;9:1823.
Trotter TN, Li M, Pan Q, Peker D, Rowan PD, Li J, et al. Myeloma cell-derived Runx2 promotes myeloma progression in bone. Blood. 2015;125:3598–608.
Fu J, Li S, Feng R, Ma H, Sabeh F, Roodman GD, et al. Multiple myeloma-derived MMP-13 mediates osteoclast fusogenesis and osteolytic disease. J Clin Invest. 2016;126:1759–72.
Seckinger A, Hillengass J, Emde M, Beck S, Kimmich C, Dittrich T, et al. CD38 as immunotherapeutic target in light chain amyloidosis and multiple myeloma-association with molecular entities, risk, survival, and mechanisms of upfront resistance. Front Immunol. 2018;9:1676.
Kaiser MF, Johnson DC, Wu P, Walker BA, Brioli A, Mirabella F, et al. Global methylation analysis identifies prognostically important epigenetically inactivated tumor suppressor genes in multiple myeloma. Blood. 2013;122:219–26.
Liang Z, Pan L, Shi J, Zhang L. C1QA, C1QB, and GZMB are novel prognostic biomarkers of skin cutaneous melanoma relating tumor microenvironment. Sci Rep. 2022;12:20460.
Shao Y, Zheng Z, Li S, Yang G, Qi F, Fei F. Upregulation of EMID1 accelerates to a favorable prognosis and immune infiltration in lung adenocarcinoma. J Oncol. 2022;2022:5185202.
Lesniak W, Wilanowski T, Filipek A. S100A6 - focus on recent developments. Biol Chem. 2017;398:1087–94.
Nedjadi T, Kitteringham N, Campbell F, Jenkins RE, Park BK, Navarro P, et al. S100A6 binds to annexin 2 in pancreatic cancer cells and promotes pancreatic cancer cell motility. Br J Cancer. 2009;101:1145–54.
Purushothaman A, Sanderson RD. Heparanase: a dynamic promoter of myeloma progression. Adv Exp Med Biol. 2020;1221:331–49.
Acknowledgements
This work was supported by the Interdisciplinary Center for Clinical Research (IZKF), Medical Faculty of Würzburg, grant B-409 (F.J., R.E.), Z-6 (T.B.), and the German Research Foundation (DFG) grants JU426/6-1 (F.J., 370022528), JU426/10-1 and EB447/10-1 (F.J., R.E, 491715122). This publication was supported by the Open Access Publication Fund of the University of Wuerzburg.
Author information
Authors and Affiliations
Contributions
Contribution: FJ developed the concept; DSM and FJ designed the experiments; DSM, JAK, WHC., and MK performed the experiments; MK, CPA, FJ, DH, AS, RE, and FJ provided technical support, advice, and supervision; DSM, JAK, SB, WHC, MK, TB, DH, AS, RE, and F.J. analyzed the data; DSM and FJ wrote the manuscript. All authors approved the final version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Maichl, D.S., Kirner, J.A., Beck, S. et al. Identification of NOTCH-driven matrisome-associated genes as prognostic indicators of multiple myeloma patient survival. Blood Cancer J. 13, 134 (2023). https://doi.org/10.1038/s41408-023-00907-6
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41408-023-00907-6