Front. Oncol.Frontiers in OncologyFront. Oncol.2234-943XFrontiers Media S.A.10.3389/fonc.2023.1105196OncologyOriginal ResearchDeveloping mRNA signatures as a novel prognostic biomarker predicting high risk multiple myelomaWangJing123*†GuoLili4†LvChenglan2†ZhouMin2†WanYuan3*1Department of Oncology and Hematology, Yizheng Hospital of Nanjing Drum Tower Hospital Group, Yizheng, China2Department of Hematology, The Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing, China3The Pq Laboratory of BiomeDx/Rx, Department of Biomedical Engineering, Binghamton University, State University of New York (SUNY), Binghamton, NY, United States4The Second Hospital of Nanjing, Nanjing University of Chinese Medicine, Nanjing, China
Edited by: Zhiqiang Liu, Tianjin Medical University, China
Reviewed by: Chandraditya Chakraborty, Dana-Farber Cancer Institute and Harvard Medical School, United States; Xin Li, Houston Methodist Research Institute, United States; Ying Xie, Tianjin Medical University, China; Jinsong Hu, Xi’an Jiaotong University, China
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Background
Multiple myeloma (MM) remains an essentially incurable disease. This study aimed to establish a predictive model for estimating prognosis in newly diagnosed MM based on gene expression profiles.
Methods
RNA-seq data were downloaded from the Multiple Myeloma Research Foundation (MMRF) CoMMpass Study and the Genotype-Tissue Expression (GTEx) databases. Weighted gene coexpression network analysis (WGCNA) and protein-protein interaction network analysis were performed to identify hub genes. Enrichment analysis was also conducted. Patients were randomly split into training (70%) and validation (30%) datasets to build a prognostic scoring model based on the least absolute shrinkage and selection operator (LASSO). CIBERSORT was applied to estimate the proportion of 22 immune cells in the microenvironment. Drug sensitivity was analyzed using the OncoPredict algorithm.
Results
A total of 860 newly diagnosed MM samples and 444 normal counterparts were screened as the datasets. WGCNA was applied to analyze the RNA-seq data of 1589 intersecting genes between differentially expressed genes and prognostic genes. The blue module in the PPI networks was analyzed with Cytoscape, and 10 hub genes were identified using the MCODE plug-in. A three-gene (TTK, GINS1, and NCAPG) prognostic model was constructed. This risk model showed remarkable prognostic value. CIBERSORT assessment revealed the risk model to be correlated with activated memory CD4 T cells, M0 macrophages, M1 macrophages, eosinophils, activated dendritic cells, and activated mast cells. Furthermore, based on OncoPredict, high-risk MM patients were sensitive to eight drugs.
Conclusions
We identified and constructed a three-gene-based prognostic model, which may provide new and in-depth insights into the treatment of MM patients.
Multiple myeloma (MM) is the second most common malignant hematological disease, accounting for approximately 10% of all hematologic malignancies (1). Our understanding of the biological mechanisms underlying the development of MM has advanced greatly. Survival in multiple myeloma has improved significantly in the last decade (2), but MM remains an essentially incurable disease. It is therefore of great clinical importance to find novel molecular markers for molecular targeted treatment of MM. MM is an invariably fatal disease with a highly heterogeneous outcome because of heterogeneous genomes. Transcriptome and proteome maps will accelerate the discovery of new therapeutic targets based on disease biology and the identification of biomarkers to guide therapeutic decisions in MM. Gene expression profiling (GEP) is a useful tool to estimate the aggressiveness of MM and will help to make individualized therapeutic decisions (3). Many different gene expression-based prognostic signatures have been reported for MM in the last decade (4–8). In 2011, the Multiple Myeloma Research Foundation (MMRF) CoMMpass Study was initiated, which gathered information on close to 1200 patients aged 27 to 93 years and followed up on a biannual basis for at least 8 years (9). The main objectives of our study were to establish universal prognostic gene signatures to enable the stratification of newly diagnosed MM patients at higher risk based on the MMRF CoMMpass study, which may provide new biomarkers serving as druggable targets for the treatment of high-risk MM patients.
Materials and methodsData acquisition and preprocessing
We obtained the RNA-seq transcriptome data of 860 newly diagnosed MM samples and 444 normal samples from the MMRF CoMMpass and Genotype-Tissue Expression (GTEx) database. Normalized read counts to these assemblies were calculated using transcripts per kilobase million (TPM) values. Clinical information (age, sex, ISS stage, survival time and status) was also collected from MMRF. Then, the differentially expressed genes (DEGs) were selected by using the “limma” package (10) with | log FC | ≥ 1 and p value < 0.05.
Identification of prognostic genes
Through univariate Cox analysis, the association between expression levels of genes and MM patients’ overall survival (OS) was explored. Genes with P < 0.05 based on the univariate analysis were identified as PGs. A Venn diagram was drawn to select intersecting genes between PGs and DEGs.
Functional enrichment analysis
Gene Ontology (GO) functional enrichment analysis was performed using the R package “clusterProfiler” (11).
Weighted gene co-expression network analysis
The construction of the gene coexpression network was completed using the “WGCNA” package (12) to explore the correlation of intersecting genes and search for important interacting gene modules. The correlation of gene expression profiles and module eigengenes was represented by module membership.
Generation and analysis of a protein-protein interaction network
A PPI network of module eigengenes was constructed using the String database (https://www.string-db.org/). Molecular Complex Detection (MCODE) in Cytoscape (Version 3.8.0,RRID : SCR_003032) was applied to screen hub genes with MCODE score >3 and number of nodes >4.
Construction and verification of the risk model
A total of 756 patients with MM were randomly split into a training cohort (n = 530) and a validation cohort (n = 226). LASSO Cox performs collinearity processing on the filtered Hub genes. The coefficient of more important survival-associated genes is determined when the adding appropriate penalty (lambda) is the minimum. An optimal cutoff was identified via the method of maximally selected rank statistics to develop a prognosis classifier for MM patients in training cohort. The GEO database, including GSE4581 and GSE57317, validated the prognostic model. A receiver operating characteristic (ROC) curve was applied to evaluate the predictive performance of our risk model compared with other models (4–6).
Construction and prediction of the nomogram
A clinical nomogram was developed to predict OS using the “rms” package. The calibration curve was applied to evaluate the consistency between the nomograms. Decision curve analysis (DCA) was used to evaluate the nomogram.
Gene set enrichment analysis (RRID : SCR_003199)
Reactome GSEA pathway (13) interaction analysis was run to compare the gene expression profiles of different risk groups.
Tumor-infiltrating immune cell analysis
The relative abundance of 22 TIIC subpopulations was estimated by applying CIBERSORT (14) in the MM high- and low-risk groups.
Predictions for drug sensitivity analysis
The R package “OncoPredict” (15) of 198 drugs was used to predict in vivo drug responses in high-risk MM patients.
Statistical analysis
All statistical analyses were performed using packages developed in R 4.1.1 (R Project for Statistical Computing, RRID : SCR_001905). A p value of less than 0.05 was used for statistical significance.
ResultsIdentification of DEGs
A heatmap of DEGs according to gene expression in 860 MM samples and 444 normal counterparts is shown in Figure 1A. In total, 4672 overexpressed genes and 937 underexpressed genes were identified (Figure 1B). GO functional enrichment analysis of DEGs is shown in Figure 1C.
(A) Heatmap was used to visualize DEGs. (B) A volcano plot was generated to visualize DEGs. (C) GO functional enrichment analysis of DEGs. (D) A Venn diagram was used to identify intersecting genes between PGs and DEGs.
Intersecting genes between DEGs and PGs
We analyzed the correlation between the expression of each gene and the overall survival of MM patients to find 9564 PGs. A total of 1589 genes were classified as intersecting genes between DEGs and PGs (Figure 1D).
WGCNA
A total of 1589 intersecting genes were included in WGCNA. The soft threshold of β=5 was identified to construct a scale-free network (Figure 2A). Under the clustering criteria of mergecutheight and minModuleSize of 0.25 and 30, respectively, a total of 2 modules (blue and turquoise) were obtained (Figure 2B). Module blue was the most relevant module for trait (Figures 2C, D).
WGCNA. (A) The soft threshold power of WGCNA was determined based on the scale-free fitting index R2. (B) Clustering dendrograms of genes based on coexpression network analysis. (C) A heatmap showing the correlation between gene modules and clinical features. (D) A scatter plot of gene significance versus the module membership in the blue module.
PPI network construction and hub gene selection
A total of 119 genes found in the blue module were imported into the STRING database to obtain the interaction relationships. In total, there were 10 hub genes from the blue modules using the MCODE plug-in of Cytoscape for subsequent analysis. Visualization results are shown in Figure 3A. GO analysis showed that the hub genes were mainly involved in mitotic nuclear division, as illustrated in Figure 3B.
(A) Ten hub genes identified by the MCODE plug-in. (B) GO functional enrichment analysis of hub genes.
Risk prediction model construction and validation
In the training set, the dimensionality reduction of hub genes was implemented by Lasso regression, as presented in Figures 4A, B, and a prognostic model was constructed based on three genes (TTK, GINS1, and NCAPG). The identified risk scoring equation is as follows: Risk score = TTK×0.0770053593671592+ GINS1×1.1810301963703+ NCAPG×1.09399274693476. There were significant differences between the survival of patients with different risks in the training set, validation set and whole MMRF cohort (p < 0.001, Figures 4C-E). The same results are shown in GSE4581 cohort and GSE57317 cohort (Figures 4F, G). ROC curves were applied for the prediction accuracy of different models (Figure 4H).
(A) LASSO for risk factor screening. (B) LASSO variable trajectory diagram. (C) KM curves of the training set. (D) KM curves of the validation set. (E) KM curves of the whole MMRF CoMMpass. (F) KM curves of the GSE4581 cohort. (G) KM curves of the GSE57317 cohort. (H) ROC curve for validation.
Nomogram construction and validation
A nomogram was constructed according to the contributions of age, sex, ISS staging, and risk status, as shown in Figure 5A. The risk status was evaluated by the risk scoring equation. The calibration curve showed that the predicted probability of the nomogram was consistent with the actual observed probability (Figure 5B). The DCA curve demonstrated that the nomogram was clinically useful (Figure 5C).
(A) Nomogram-based prognosis prediction model. (B) Calibration curve of the nomogram. (C) Decision curve analysis of the nomogram.
TIIC analysis
The differential abundance of TIICs is shown in the heatmap (Figure 6A). Six TIICs (activated memory CD4 T cells, M0 macrophages, M1 macrophages, eosinophils, activated dendritic cells, and activated mast cells) showed significant differences between different risk groups (Figure 6B). The results also showed that the 3 genes (TTK, GINS and NCAPG) were correlated with TIICs (Figure 6C).
(A) Heatmap showing the proportion of TIICs in MM samples. (B) Violin plot showing the ratio differentiation of TIICs between MM samples with low or high risk. (C) Correlation of the expression of three genes with the expression of TIICs.
GSEA
Reactome GSEA pathway interaction analysis revealed that 22 significant pathways were dramatically changed in the high-risk group (Figure 7), involving cancer, the immune system, metabolism, and signal transduction.
GSEA results showed that the reactome pathway was enriched in MM samples from the high-risk group.
Suitable drugs for high-risk patients
The OncoPredict algorithm showed that high-risk MM patients were more sensitive to 8 drugs, including Nutlin-3a (MDM2 inhibitor), SB216763 (GSK3 inhibitor), oxaliplatin (platinum anticancer drug), olaparib (PARP inhibitor), irinotecan (TopoI inhibitor), BMS-754807 (IGF-1R/IR inhibitor), AZD8055 (mTOR inhibitor), and camptothecin (TopoI inhibitor) (Figure 8).
Predicted sensitivity scores of drugs that are candidate therapeutic agents for high-risk MM patients (*p<0.05, **p<0.01, ***p<0.001, ****p<0.0001).
Discussion
An in-depth understanding of the biological mechanisms underlying carcinogenesis in MM is vital to find strategies for MM treatment. To the best of our knowledge, there are few studies on the role of three genes (TTK, GINS1, and NCAPG) in MM. TTK is an integral part of the spindle assembly checkpoint, which is responsible for maintaining the integrity of the genome (16). TTK is hardly detectable in normal tissues except the testis and placenta. TTK is an indicator of poor prognosis in breast cancer (17), hepatocellular carcinoma (18), lung cancer (19), and glioma (20). Inhibition of TTK could induce MM cell (AMO-1) apoptosis and deregulate the proliferation in vitro (21). GINS1 is a member of the GINS complex that plays a vital role in DNA replication (22). GINS1 was found to be related to poor prognosis in breast and liver cancer (23, 24). GINS1 physically interacts with TOP2A (25), which might be a highly significant predictor of response to proteasome inhibitors (26). NCAPG, a subunit of condensin I, is responsible for the condensation or structure of mitotic chromosomes (27). The expression level of NCAPG has been shown to be closely related to the prognosis of tumors (28–31). NCAPG involves in the regulation of different signaling pathways, such as PI3K/AKT (32), NF-κB (33) and SRC/STAT3 (34) signaling pathway. However, the effects of GINS1 and NCAPG in MM remain unknown.
The CIBERSORT algorithm was used to estimate immune cell infiltration in the bone marrow (BM) microenvironment. The proliferation, progression, and survival of malignant plasma cells (PCs) in MM are highly affected by the BM microenvironment (35). BM microenvironment is highly enriched for suppressive immune cells such as MDSCs, Tregs, pDCs, Bregs, N2 neutrophils, M2 macrophages, which leads to effector cell dysfunction and lack of persistence (36). We found that M2 macrophages and regulatory T cells (Tregs) had apparent abnormal infiltration in MM, contributing to the immune evasion, proliferation, and drug resistance of MM cells (37). We also found that these immunosuppressive cells are in bidirectional and multidirectional crosstalk inhibited memory effector T-cell populations. Eosinophils (Eos) can contribute to the proliferation of malignant PCs in MM (38). Suzuki et al. study showed that the median time to next treatment (TTNT) in the elevated eosinophil group was significantly longer than that in the nonelevated group (40.3 months vs 8.4 months; P = .017) in relapsed or refractory myeloma patients treated with lenalidomide (39). Mast cell (MC) accumulation correlates with increased neovascularization in MM (40). Raised IL-6 levels can be caused by the presence of MC (41), and it has been identified as the key growth and survival factor for myeloma cells (42). Tumor lesions in MM are highly infiltrated by dendritic cells (43). TTK could affect the tumor microenvironment (TME) by affecting the number of immune cells (44). No literature is available on the role of GINS1 and NCAPG in the TME.
Reactome GSEA pathway interaction analysis revealed that several pathways were dramatically changed in the high-risk myeloma patients. The IL-21 signaling pathway in myeloma cells involves phosphorylation of Erk1/2, Jak1, and Stat3 (45). The upregulation of chaperone-mediated autophagy (CMA) is a potential mechanism of resistance to bortezomib (46). The pentose phosphate pathway protects against oxidative stress-mediated late apoptosis/necrosis of multiple myeloma cells (47).
Our study showed that 8 drugs were effective in high-risk MM patients. Nutlin-3 can disrupt the p53-MDM2 interaction and activate p53. Nutlin-3 with bortezomib may increase clinical responsiveness to bortezomib-based therapy (48). The chemical inhibitor SB216763 leads to a reduction in MM cell growth and augments the response of MM cells to the cytotoxic effects of bortezomib (49). Oxaliplatin triggers bona fide “immunogenic cell death” (ICD), as it provokes a premortem endoplasmic reticulum stress response (50), so it can be used as an inducer of ICD and as a modulator of the TME (51). Drugs that induce excessive amounts of ER stress, such as proteasome inhibitors and novel ER stressors, are predicted to be very effective in targeting MM cells (52). The PARP1 inhibitor olaparib impaired MM cell viability in vitro and was effective against MM in vivo xenografts (53). Efficacy needs to be further evaluated in clinical trials.
Certainly, our research still has some limitations. First, only transcriptome data were included for the MM prognosis study. Second, some clinical data are missing, such as cytogenetic profile and type of measurable disease. Third, the carcinogenic mechanism of the three genes in MM warrants intensive study. Fourth, we used algorithm analyses to predict drug sensitivity but did not verify it.
In conclusion, our findings may improve the understanding of the factors that influence development and prognosis in MM, which may present a new strategy.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.
Ethics statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
JW and YW conceived and designed the study. JW, LG, CL, and MZ performed the electronic search, data collection, abstraction, and data entry. JW, LG, MZ, and CL were statistical advisers. JW, and YW were responsible for the overall direction of the text and discussion. JW and YW had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. All authors contributed to the article and approved the submitted version.
Funding
The work was partially supported by National Cancer Institute (1R01CA230339 and 1R37CA255948). The sponsors had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
ReferencesRajkumarSVKumarS. Multiple myeloma current treatment algorithms. Blood Cancer J (2020) 10(9):94. doi: 10.1038/s41408-020-00359-2RajkumarSV. Multiple myeloma: 2022 update on diagnosis, risk stratification, and management. Am J Hematol (2022) 97(8):1086–107. doi: 10.1002/ajh.26590NCCN. Clinical practice guidelines in oncology, multiple myeloma. version 4.2021 (NCCN).WangWXuSWZhuXYGuoQYZhuMMaoXL. Identification and validation of a novel RNA-binding protein-related gene-based prognostic model for multiple myeloma. Front Genet (2021) 12:665173. doi: 10.3389/fgene.2021.665173ChenXLiuLChenMXiangJWanYLiX. A five-gene risk score model for predicting the prognosis of multiple myeloma patients based on gene expression profiles. Front Genet (2021) 12:785330. doi: 10.3389/fgene.2021.785330LiuYYangNPengXLiuGZhongHLiuL. One-lincRNA and five-mRNA based signature for prognosis of multiple myeloma patients undergoing proteasome inhibitors therapy. BioMed Pharmacother (2019) 118:109254. doi: 10.1016/j.biopha.2019.109254WangJZuoYLvCZhouMWanY. N6-methyladenosine regulators are potential prognostic biomarkers for multiple myeloma. IUBMB Life (2022) 75(2):137–48. doi: 10.1002/iub.2678ChenCLiYMiaoPXuYXieYChenZ. Tumor immune cell infiltration score based model predicts prognosis in multiple myeloma. Sci Rep (2022) 12(1):17082. doi: 10.1038/s41598-022-21763-7ShawML. The MMRF: understanding and sharing the genomic landscape of multiple myeloma. Am J Manag Care (2021) 27(9 Spec No.):SP392. doi: 10.37765/ajmc.2021.88821RitchieMEPhipsonBWuDHuYLawCWShiW. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43(7):e47. doi: 10.1093/nar/gkv007WuTHuEXuSChenMGuoPDaiZ. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (Camb) (2021) 2(3):100141. doi: 10.1016/j.xinn.2021.100141LangfelderPHorvathS. WGCNA: an r package for weighted correlation network analysis. BMC Bioinf (2008) 9:559. doi: 10.1186/1471-2105-9-559JassalBMatthewsLViteriGGongCLorentePFabregatA. The reactome pathway knowledgebase. Nucleic Acids Res (2020) 48(D1):D498–503. doi: 10.1093/nar/gkz1031ChenBKhodadoustMSLiuCLNewmanAMAlizadehAA. Profiling tumor infiltrating immune cells with CIBERSORT. Methods Mol Biol (2018) 1711:243–59. doi: 10.1007/978-1-4939-7493-1_12MaeserDGruenerRFHuangRS. oncoPredict: an r package for predicting in vivo or cancer patient drug response and biomarkers from cell line screening data. Brief Bioinform (2021) 22(6):bbab260. doi: 10.1093/bib/bbab260MusacchioA. The molecular biology of spindle assembly checkpoint signaling dynamics. Curr Biol (2015) 25(20):R1002–18. doi: 10.1016/j.cub.2015.08.051GaoYHQuSSCaoLQYaoM. TTK predicts triple positive breast cancer prognosis and regulates tumor proliferation and invasion. Neoplasma (2022) 69(2):274–82. doi: 10.4149/neo_2021_210421N541LiuXLiaoWYuanQOuYHuangJ. TTK activates akt and promotes proliferation and migration of hepatocellular carcinoma cells. Oncotarget (2015) 6(33):34309–20. doi: 10.18632/oncotarget.5295TsaiYMWuKLChangYYHungJYChangWAChangCY. Upregulation of Thr/Tyr kinase increases the cancer progression by neurotensin and dihydropyrimidinase-like 3 in lung cancer. Int J Mol Sci (2020) 21(5):1640. doi: 10.3390/ijms21051640ChenSWangYNiCMengGShengX. HLF/miR-132/TTK axis regulates cell proliferation, metastasis and radiosensitivity of glioma cells. BioMed Pharmacother (2016) 83:898–904. doi: 10.1016/j.biopha.2016.08.004de BoussacHBruyerAJourdanMMaesARobertNGourzonesC. Kinome expression profiling to target new therapeutic avenues in multiple myeloma. Haematologica (2020) 105(3):784–95. doi: 10.3324/haematol.2018.208306ChangYPWangGBermudezVHurwitzJChenXS. Crystal structure of the GINS complex and functional insights into its role in DNA replication. Proc Natl Acad Sci U S A. (2007) 104(31):12685–90. doi: 10.1073/pnas.0705558104Nieto-JiménezCAlcaraz-SanabriaAPáezRPérez-PeñaJCorrales-SánchezVPandiellaA. DNA-Damage related genes and clinical outcome in hormone receptor positive breast cancer. Oncotarget (2016) 8(38):62834–41. doi: 10.18632/oncotarget.10886LianYFLiSSHuangYLWeiHChenDMWangJL. Up-regulated and interrelated expressions of GINS subunits predict poor prognosis in hepatocellular carcinoma. Biosci Rep (2018) 38(6):BSR20181178. doi: 10.1042/BSR20181178YangHLiuXZhuXZhangMWangYMaM. GINS1 promotes the proliferation and migration of glioma cells through USP15-mediated deubiquitination of TOP2A. iScience (2022) 25(9):104952. doi: 10.1016/j.isci.2022.104952RealeAKhongTMithraprabhuSSavvidouIHockingJBerginK. TOP2A expression predicts responsiveness to carfilzomib in myeloma and informs novel combinatorial strategies for enhanced proteasome inhibitor cell killing. Leuk Lymphoma. (2021) 62(2):337–47. doi: 10.1080/10428194.2020.1832659MurphyLASargeKD. Phosphorylation of CAP-G is required for its chromosomal DNA localization during mitosis. Biochem Biophys Res Commun (2008) 377(3):1007–11. doi: 10.1016/j.bbrc.2008.10.114GuoZYZhuZT. NCAPG is a prognostic biomarker associated with vascular invasion in hepatocellular carcinoma. Eur Rev Med Pharmacol Sci (2021) 25(23):7238–51. doi: 10.26355/eurrev_202112_27417AraiTOkatoAYamadaYSugawaraSKurozumiAKojimaS. Regulation of NCAPG by miR-99a-3p (passenger strand) inhibits cancer cell aggressiveness and is involved in CRPC. Cancer Med (2018) 7(5):1988–2002. doi: 10.1002/cam4.1455DongMXuTCuiXLiHLiXXiaW. NCAPG upregulation mediated by four microRNAs combined with activation of the p53 signaling pathway is a predictor of poor prognosis in patients with breast cancer. Oncol Lett (2021) 21(4):323. doi: 10.3892/ol.2021.12585WuYLinYPanJTuXXuYLiH. NCAPG promotes the progression of lung adenocarcinoma via the TGF-β signaling pathway. Cancer Cell Int (2021) 21(1):443. doi: 10.1186/s12935-021-02138-wGongCAiJFanYGaoJLiuWFengQ. NCAPG promotes the proliferation of hepatocellular carcinoma through PI3K/AKT signaling. Onco Targets Ther (2019) 12:8537–52. doi: 10.2147/OTT.S217916TangFYuHWangXShiJChenZWangH. NCAPG promotes tumorigenesis of bladder cancer through NF-κB signaling pathway. Biochem Biophys Res Commun (2022) 622:101–7. doi: 10.1016/j.bbrc.2022.07.007JiangLRenLChenHPanJZhangZKuangX. NCAPG confers trastuzumab resistance via activating SRC/STAT3 signaling pathway in HER2-positive breast cancer. Cell Death Dis (2020) 11(7):547. doi: 10.1038/s41419-020-02753-xGiannakoulasNNtanasis-StathopoulosITerposE. The role of marrow microenvironment in the growth and development of malignant plasma cells in multiple myeloma. Int J Mol Sci (2021) 22(9):4462. doi: 10.3390/ijms22094462HoMXiaoAYiDZanwarSBianchiG. Treating multiple myeloma in the context of the bone marrow microenvironment. Curr Oncol (2022) 29(11):8975–9005. doi: 10.3390/curroncol29110705UckunFM. Overcoming the immunosuppressive tumor microenvironment in multiple myeloma. Cancers (Basel). (2021) 13(9):2018. doi: 10.3390/cancers13092018WongTWKitaHHansonCAWaltersDKArendtBKJelinekDF. Induction of malignant plasma cell proliferation by eosinophils. PloS One (2013) 8(7):e70554. doi: 10.1371/journal.pone.0070554SuzukiKNishiwakiKGunjiTKatoriMMasuokaHYanoS. Elevated eosinophil level predicted long time to next treatment in relapsed or refractory myeloma patients treated with lenalidomide. Cancer Med (2020) 9(5):1694–702. doi: 10.1002/cam4.2828RibattiDTammaRVaccaA. Mast cells and angiogenesis in human plasma cell malignancies. Int J Mol Sci (2019) 20(3):481. doi: 10.3390/ijms20030481DevetzoglouMVyzoukakiRKokonozakiMXekalouAPappaCAPapadopoulouA. High density of tryptase-positive mast cells in patients with multiple myeloma: correlation with parameters of disease activity. Tumour Biol (2015) 36(11):8491–7. doi: 10.1007/s13277-015-3586-9JasrotiaSGuptaRSharmaAHalderAKumarL. Cytokine profile in multiple myeloma. Cytokine (2020) 136:155271. doi: 10.1016/j.cyto.2020.155271ChauhanDSinghAVBrahmandamMCarrascoRBandiMHideshimaT. Functional interaction of plasmacytoid dendritic cells with multiple myeloma cells: a therapeutic target. Cancer Cell (2009) 16(4):309–23. doi: 10.1016/j.ccr.2009.08.019ChenKXingJYuWXiaYZhangYChengF. Identification and validation of hub genes associated with bladder cancer by integrated bioinformatics and experimental assays. Front Oncol (2021) 11:782981. doi: 10.3389/fonc.2021.782981BrenneATRoTBWaageASundanABorsetMHjorth-HansenH. Interleukin-21 is a growth and survival factor for human myeloma cells. Blood (2002) 99(10):3756–62. doi: 10.1182/blood.V99.10.3756NikesitchNRebeiroPHoLLPothulaSWangXMKhongT. The role of chaperone-mediated autophagy in bortezomib resistant multiple myeloma. Cells (2021) 10(12):3464. doi: 10.3390/cells10123464YinLKosugiMKufeD. Inhibition of the MUC1-c oncoprotein induces multiple myeloma cell death by down-regulating TIGAR expression and depleting NADPH. Blood (2012) 119(3):810–6. doi: 10.1182/blood-2011-07-369686OoiMGHaydenPJKotoulaVMcMillinDWCharalambousEDaskalakiE. Interactions of the Hdm2/p53 and proteasome pathways may enhance the antitumor activity of bortezomib. Clin Cancer Res (2009) 15(23):7153–60. doi: 10.1158/1078-0432.CCR-09-1071PiazzaFManniSTubiLQMontiniBPavanLColpoA. Glycogen synthase kinase-3 regulates multiple myeloma cell growth and bortezomib-induced cell death. BMC Canc (2010) 10:526. doi: 10.1186/1471-2407-10-526KeppOMengerLVacchelliELocherCAdjemianSYamazakiT. Crosstalk between ER stress and immunogenic cell death. Cytokine Growth Factor Rev (2013) 24(4):311–8. doi: 10.1016/j.cytogfr.2013.05.001ZhuHShanYGeKLuJKongWJiaC. Oxaliplatin induces immunogenic cell death in hepatocellular carcinoma cells and synergizes with immune checkpoint blockade therapy. Cell Oncol (Dordr). (2020) 43(6):1203–14. doi: 10.1007/s13402-020-00552-2Serrano-Del ValleAAnelANavalJMarzoI. Immunogenic cell death and immunotherapy of multiple myeloma. Front Cell Dev Biol (2019) 7:50. doi: 10.3389/fcell.2019.00050CaraccioloDSciontiFJuliGAltomareEGolinoGTodoertiK. Exploiting MYC-induced PARPness to target genomic instability in multiple myeloma. Haematologica (2021) 106(1):185–95. doi: 10.3324/haematol.2019.240713