Skip to main content

Identification of candidate targets and mechanisms involved in miRNA regulation in multiple myeloma

Abstract

Background

Multiple myeloma (MM) is a complex disease affected by many factors. The recognition of miRNA networks is helpful for specific detection and personalised treatment.

Methods

mRNA expression profiles were obtained from GSE39754 and GSE87830, and differentially expressed mRNAs (DEmRs) between MM and controls were identified. The intersection of the two sets of DEmRs in GSE39754 and GSE87830 was identified as common mRNAs, and enrichment analysis was subsequently performed. Moreover, we analysed differentially expressed miRNAs (DEmiRs) between MM and controls in GSE87830. A regulatory network of target mRNAs related to the overall survival of MM patients was then constructed.

Results

In this study, a total of 356 common mRNAs were identified that were significantly enriched in neutrophil-mediated immunity, Th17 cell differentiation and PI3K-Akt signalling pathways. Moreover, we identified 103 DEmiRs and predicted 91 differentially expressed mRNAs as target mRNAs. Cox regression analysis was used to screen 14 target mRNAs that significantly affected the survival of MM patients. In the constructed integrated regulatory network, HIF1A and THBS1 were found to participate in Th17 cell differentiation and PI3K-Akt signalling pathways.

Conclusion

These findings improve the understanding of the regulatory mechanisms of MM. Genes that are part of integrated regulatory networks may represent candidate targets for MM treatment.

Introduction

Multiple myeloma (MM) is a common malignant tumour derived from plasma cells in the bone marrow [1]. Since plasma cells develop to the final functional stage of B lymphocytes, MM can also be classified as B lymphocyte lymphoma [2]. Currently, the World Health Organisation classifies MM as a type of B cell lymphoma [3]. MM is characterised by abnormal proliferation of bone marrow plasma cells accompanied by excessive production of monoclonal immunoglobulin or light chain (M protein). A small number of patients may have MM that does not produce M protein [4]. The occurrence of MM is often accompanied by extensive osteolytic damage, hypercalcaemia, anaemia and kidney damage [5]. The production of normal immunoglobulin is inhibited, so various bacterial infections are prone to occur. According to statistics, the mortality rate of MM is 48.5% [6]. The 2019 Cancer Survey showed that amongst 32,110 patients with myeloma in the USA, 12,960 patients died, and the incidence and mortality of the disease were much higher in men than in women [7]. Many effective treatments have been developed in clinical studies over the past few decades [8]. MM easily relapses, causing MM to still be defined as an incurable disease [9]. Therefore, it is particularly important to explore new treatment methods for MM and identify specific targets for MM treatment.

MicroRNAs (miRNAs) are a class of noncoding single-stranded RNA molecules with a length of approximately 22 nucleotides. They are widely expressed in animals and plants and play important roles in both physiological and pathological processes [10]. To date, studies have revealed that changes in miRNA expression are present in a variety of tumour types [11]. Moreover, studies have shown that miRNAs have a regulatory effect on tumour progression [12]. Meanwhile, miRNAs can also represent early warning indicators of tumour occurrence, treatment and prognosis [13]. In MM research, investigators have revealed the application of miRNAs in the diagnosis of tumour development and potential tumour treatment, as well as part of the mechanism for evaluating the prognostic status of tumours. In addition, some researchers have studied the effects of miRNAs on tumours in other animal models [14]. These studies also indicated that the application and transformation of miRNAs in MM has important potential [15]. However, although researchers have revealed a large number of miRNAs associated with MM, there are still no specific miRNAs for the diagnosis, treatment or prognosis of MM. For that reason, we collected RNA sequencing data from normal and MM tissues from a public database to determine changes in differential genes. This study provides a theoretical basis for future research on miRNAs in the treatment of MM.

Materials and methods

Data collection

All data were obtained from the Gene Expression Omnibus (GEO) database. GSE39754 included the mRNA expression profile of CD138 purified myeloma plasma cells from 170 newly diagnosed MM patients and CD138 purified plasma cells from 6 healthy donors. GSE87830 included mRNA and miRNA expression profiles of CD138 purified myeloma plasma cells from 95 newly diagnosed MM patients and 4 noncancer samples.

Differential analysis

The limma R package [16] was utilised to identify differentially expressed mRNAs (DEmRs) and differentially expressed miRNAs (DEmiRs) between MM and controls. For these analyses, |log2-fold change (FC) | >1 and P < 0.05 were considered screening conditions.

Functional and pathway enrichment

For enrichment of mRNAs, the clusterProfiler R package [17] was used. These included Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Biological process (BP), cellular component (CC) and molecular function (MF) were included in GO terms. Enrichment results with P < 0.05 were used as the threshold of statistical significance.

Target prediction and Cox regression

The miRTarget database was used to predict the target mRNAs regulated by the identified miRNAs. Cox regression analysis for genes in the MMRF-CoMMpass database was performed to identify mRNAs that significantly influenced MM patient prognosis.

Results

To identify dysregulated mRNAs in MM patients, we performed differential analysis of gene expression profiles of MM and control samples. We obtained a total of 2808 differentially expressed mRNAs in GSE39754 (Fig. 1A). In GSE87830, we obtained a total of 1045 differentially expressed mRNAs (Fig. 1B). Amongst the two sets of DEmRs, we identified 356 common mRNAs with potentially stronger associations with MM (Fig. 1C). The expression of common mRNAs exhibited distinct differences between MM and control samples (Fig. 1D).

Fig. 1
figure 1

Identification of differentially expressed mRNAs between multiple myeloma and controls. A Volcano plot of differentially expressed mRNAs between multiple myeloma and controls in GSE39754. Red indicates upregulation in MM, and blue indicates downregulation. B Volcano plot of differentially expressed mRNAs for multiple myeloma and controls in GSE87830. Red indicates upregulation in MM, and blue indicates downregulation. C Venn diagram of differentially expressed mRNAs between two DEmRs groups. The intersection includes the common mRNAs. D Expression heatmap of common mRNAs in multiple myeloma and control samples of GSE39754

Enrichment of common mRNAs

To identify the biological functions and signalling pathways associated with MM, we performed enrichment analysis on the common mRNAs. GO results indicated that common mRNAs were primarily enriched in biological processes of neutrophil activation involved in the immune response, neutrophil degranulation and neutrophil-mediated immunity (Fig. 2A). In addition, ficolin-1-rich granules, tertiary granules and ficolin-1-rich granule lumens were primarily enriched in cellular components. In molecular function, low-density lipoprotein receptor activity, lipoprotein particle receptor activity and low-density lipoprotein particle binding were enriched. Moreover, KEGG signalling pathway analysis revealed that cell adhesion molecules (CAMs), Th17 cell differentiation and PI3K-Akt signalling pathways were significantly regulated by common genes (Fig. 2B).

Fig. 2
figure 2

GO functions and KEGG signalling pathways involved in common mRNAs. A The primary biological processes, cellular components and molecular functions enriched by common mRNAs. B KEGG pathway in which common mRNAs are involved

Construction of a miRNA regulatory network

In the differential analysis results for GSE87830, we also identified 103 miRNAs that were differentially expressed between MM and controls (Fig. 3A). These included 27 upregulated DEmiRs and 76 downregulated DEmiRs (Fig. 3B). Using the miRTarget database, we predicted 4255 target mRNAs for 103 DEmiRs. By comparison to common mRNAs, we identified 91 differentially expressed target mRNAs (Fig. 3C). Furthermore, we constructed a regulatory network of DEmiRNA target mRNAs (Fig. 3D).

Fig. 3
figure 3

MiRNA-regulated mRNA related to multiple myeloma. A Differentially expressed miRNAs between multiple myeloma and controls in GSE87830. B Volcano plot of differentially expressed miRNAs for multiple myeloma and controls. Red shows upregulation, and blue shows downregulation. C Intersection of target mRNAs predicted using the miRTarget database and common mRNAs. D DEmiRs-target common mRNA regulated network. Triangles are miRNAs and ellipses are mRNAs. Red indicates upregulated expression in MM, and blue indicates downregulated expression

Identification of key regulatory networks

Furthermore, using Cox regression analysis, we identified 14 mRNAs amongst the target mRNAs that significantly affected survival in MM patients (Fig. 4A). Amongst them, TCF7L2, PLD1 and PLA2G7 were protective factors, and other mRNAs were risk factors for prognosis of MM patients (Table 1). Finally, we performed a comprehensive analysis of the regulatory network of 14 mRNAs, revealing that hsa−miR−199a−5p regulated HIF1A, and hsa−miR−491−5p regulated THBS1, affecting most KEGG pathways (Fig. 4B). These included significantly enriched Th17 cell differentiation and PI3K−Akt signalling pathways.

Fig. 4
figure 4

Construction of key miRNA regulatory network. A Intersection genes of target mRNAs and mRNAs significantly influencing MM patient prognosis. B Comprehensive network of target mRNAs regulated by miRNAs involved in KEGG signalling pathways

Table 1 Cox regression analysis of target 14 mRNAs with significant impact on MM patient prognosis

Discussion

In this study, 356 common mRNAs and 103 DEmiRs were identified, and further examination revealed that these DEmiRs were enriched in the miR-199a-5p/HIF1A and miR-491-5p/THBS1 signalling axes. In addition, some DEmiRs were also enriched in Th17 cell differentiation and PI3K-Akt signalling pathways.

The occurrence of MM changes the expression of many genes, and studying these genes may reveal the mechanism of MM. On the other hand, it also provides new detection standards for early tumour warning, treatment and prognostic evaluation. Screening exogenous agonists and inhibitors for target genes is a major trend in tumour treatment research. Therefore, it is particularly important to screen for key genes that regulate tumorigenesis. In this study, we identified 356 common mRNAs and 103 DEmiRs that may be related to MM. Research on these genes may play a role in the treatment of MM in the future.

The role of miRNAs in cancer progression and their potential application value has received researchers’ attention. Calin et al. reported the association between miRNAs and cancer and identified a dynamic relationship between miRNAs and chromosomes in different types of cancer [18,19,20]. Due to the relatively stable structure of miRNAs, they have been leveraged as potential tools for cancer diagnosis [21]. Previous studies have also revealed changes in the expression of miRNAs in MM patients and may serve an early warning and play a diagnostic role in the process of MM [15]. For example, miR-34a effectively distinguished between healthy and MM tissues, suggesting that miR-34a can be used as an indicator to diagnose MM [22]. In this study, DEmiRs were enriched in the miR-199a-5p/HIF1A and miR-491-5p/THBS1 signalling axes. In past studies, miR-199a-5p and miR-491-5p were identified as tumour suppressors [23, 24]. However, the effects of miR-199a-5p and miR-491-5p on MM are still unknown. Our results strongly suggested that the miR-199a-5p/HIF1A and miR-491-5p/THBS1 signalling axes may represent key signalling axes in the diagnosis and treatment of MM.

Th17 cells are a subset of helper T cells involved in the process of inflammation and autoimmune disease [25]. In addition, Th17 cells play a regulatory role in tumour progression. Studies have shown that in patients with gastric cancer, Th17 cells are significantly increased in the peripheral blood [26]. Moreover, expression of Th17 cells in peripheral blood changes with tumour progression, which also suggests that Th17 cells can also be used as indicators of tumour prognosis [27]. Moreover, studies have shown that Th17 cells also mediate the occurrence of MM, and the effect of Th17 cells on MM can be eliminated through miR-21 [28]. In this study, DEmiRs were enriched in Th17 cell differentiation, further illustrating the role of Th17 cells in MM.

The PI3K-Akt signalling pathway is an axial pathway that regulates cell proliferation and survival and plays an important role in maintaining vital signs and pathological processes [29, 30]. Similarly, the PI3K-Akt signalling pathway is also involved in tumour progression [31]. Previous studies have revealed that many miRNAs regulate MM through the PI3K-Akt signalling pathway [32,33,34,35]. In this study, DEmiRs were significantly enriched in the PI3K-Akt signalling pathway, further supporting the role of miRNAs in MM.

Our study also has some limitations. First, our key results were missing validation in clinical samples. We will utilise a large number of tissue samples for experimental validation in future studies. Thus, the molecular function of key candidate markers remains to be further explored.

Conclusions

This study preliminarily explored miRNAs as potential diagnostic markers in MM. Further research revealed that DEmiRs were primarily enriched in the miR-199a-5p/HIF1A and miR-491-5p/THBS1 signalling axes. In addition, some DEmiRs were also enriched in Th17 cell differentiation and PI3K-Akt signalling pathways. This study confirmed the role of miRNAs in MM by analysing differences in miRNAs and provides a theoretical basis for future research and application of MM.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

References

  1. Ludwig H, Boccadoro M, Moreau P, San-Miguel J, Cavo M, Pawlyn C, et al. Recommendations for vaccination in multiple myeloma: a consensus of the European Myeloma Network. Leukemia. 2021;35:31–44.

    Article  Google Scholar 

  2. Bayly E, Nguyen V, Binek A, Piggin A, Baldwin K, Westerman D, et al. Validation of a modified pre-lysis sample preparation technique for flow cytometric minimal residual disease assessment in multiple myeloma, chronic lymphocytic leukemia, and B-non Hodgkin lymphoma. Cytometry B Clin Cytom. 2020;98:385–98.

    Article  CAS  Google Scholar 

  3. Sabattini E, Bacci F, Sagramoso C, Pileri SA. WHO classification of tumours of haematopoietic and lymphoid tissues in 2008: an overview. Pathologica. 2010;102:83–7.

    CAS  PubMed  Google Scholar 

  4. Kazandjian D. Multiple myeloma epidemiology and survival: a unique malignancy. Semin Oncol. 2016;43:676–81.

    Article  Google Scholar 

  5. Mimura N, Hideshima T, Anderson KC. Novel therapeutic strategies for multiple myeloma. Exp Hematol. 2015;43:732–41.

    Article  CAS  Google Scholar 

  6. Neuse CJ, Lomas OC, Schliemann C, Shen YJ, Manier S, Bustoros M, et al. Genome instability in multiple myeloma. Leukemia. 2020;34:2887–97.

    Article  Google Scholar 

  7. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA Cancer J Clin. 2019;69:7–34.

    Article  Google Scholar 

  8. Manier S, Salem KZ, Park J, Landau DA, Getz G, Ghobrial IM. Genomic complexity of multiple myeloma and its clinical implications. Nat Rev Clin Oncol. 2017;14:100–13.

    Article  CAS  Google Scholar 

  9. Smith EL, Mailankody S, Staehr M, Wang X, Senechal B, Purdon TJ, et al. BCMA-targeted CAR T-cell therapy plus radiotherapy for the treatment of refractory myeloma reveals potential synergy. Cancer Immunol Res. 2019;7:1047–53.

    Article  CAS  Google Scholar 

  10. Lu TX, Rothenberg ME. MicroRNA. J Allergy Clin Immunol. 2018;141:1202–7.

    Article  CAS  Google Scholar 

  11. Lee YS, Dutta A. MicroRNAs in cancer. Annu Rev Pathol. 2009;4:199–227.

    Article  CAS  Google Scholar 

  12. Mishra S, Yadav T, Rani V. Exploring miRNA based approaches in cancer diagnostics and therapeutics. Crit Rev Oncol Hematol. 2016;98:12–23.

    Article  Google Scholar 

  13. Qadir MI, Faheem A. miRNA: a diagnostic and therapeutic tool for pancreatic cancer. Crit Rev Eukaryot Gene Expr. 2017;27:197–204.

    Article  Google Scholar 

  14. Zhao JJ, Carrasco RD. Crosstalk between microRNA30a/b/c/d/e-5p and the canonical Wnt pathway: implications for multiple myeloma therapy. Cancer Res. 2014;74:5351–8.

    Article  CAS  Google Scholar 

  15. Gao SS, Wang YJ, Zhang GX, Zhang WT. Potential diagnostic value of circulating miRNA for multiple myeloma: a meta-analysis. J Bone Oncol. 2020;25:100327.

    Article  Google Scholar 

  16. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.

    Article  Google Scholar 

  17. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16:284–7.

    Article  CAS  Google Scholar 

  18. Bhardwaj A, Singh S, Singh AP. MicroRNA-based cancer therapeutics: big hope from small RNAs. Mol Cell Pharmacol. 2010;2:213–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, et al. Frequent deletions and down-regulation of micro-RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A. 2002;99:15524–9.

    Article  CAS  Google Scholar 

  20. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, et al. MicroRNA expression profiles classify human cancers. Nature. 2005;435:834–8.

    Article  CAS  Google Scholar 

  21. Zhu B, Ju S, Chu H, Shen X, Zhang Y, Luo X, et al. The potential function of microRNAs as biomarkers and therapeutic targets in multiple myeloma. Oncol Lett. 2018;15:6094–106.

    PubMed  PubMed Central  Google Scholar 

  22. Kubiczkova L, Kryukov F, Slaby O, Dementyeva E, Jarkovsky J, Nekvindova J, et al. Circulating serum microRNAs as novel diagnostic and prognostic biomarkers for multiple myeloma and monoclonal gammopathy of undetermined significance. Haematologica. 2014;99:511–8.

    Article  CAS  Google Scholar 

  23. Li DJ, Wang X, Yin WH, Niu K, Zhu W, Fang N. MiR-199a-5p suppresses proliferation and invasion of human laryngeal cancer cells. Eur Rev Med Pharmacol Sci. 2020;24:12200–7.

    PubMed  Google Scholar 

  24. Lu L, Cai M, Peng M, Wang F, Zhai X. miR-491-5p functions as a tumor suppressor by targeting IGF2 in colorectal cancer. Cancer Manag Res. 2019;11:1805–16.

    Article  CAS  Google Scholar 

  25. Hou F, Li Z, Ma D, Zhang W, Zhang Y, Zhang T, et al. Distribution of Th17 cells and Foxp3-expressing T cells in tumor-infiltrating lymphocytes in patients with uterine cervical cancer. Clin Chim Acta. 2012;413:1848–54.

    Article  CAS  Google Scholar 

  26. Yu Q, Lou XM, He Y. Preferential recruitment of Th17 cells to cervical cancer via CCR6-CCL20 pathway. PLoS One. 2015;10:e0120855.

    Article  Google Scholar 

  27. Zhang W, Tian X, Mumtahana F, Jiao J, Zhang T, Croce KD, et al. The existence of Th22, pure Th17 and Th1 cells in CIN and cervical cancer along with their frequency variation in different stages of cervical cancer. BMC Cancer. 2015;15:717.

    Article  CAS  Google Scholar 

  28. Rossi M, Altomare E, Botta C, Gallo Cantafio ME, Sarvide S, Caracciolo D, et al. miR-21 antagonism abrogates Th17 tumor promoting functions in multiple myeloma. Leukemia. 2021;35:823–34.

    Article  CAS  Google Scholar 

  29. Xia P, Xu XY. PI3K/Akt/mTOR signaling pathway in cancer stem cells: from basic research to clinical application. Am J Cancer Res. 2015;5:1602–9.

    PubMed  PubMed Central  Google Scholar 

  30. Ersahin T, Tuncbag N, Cetin-Atalay R. The PI3K/AKT/mTOR interactive pathway. Mol Biosyst. 2015;11:1946–54.

    Article  CAS  Google Scholar 

  31. Alzahrani AS. PI3K/Akt/mTOR inhibitors in cancer: at the bench and bedside. Semin Cancer Biol. 2019;59:125–32.

    Article  CAS  Google Scholar 

  32. Zhu B, Chen H, Zhang X, Pan Y, Jing R, Shen L, et al. Serum miR-30d as a novel biomarker for multiple myeloma and its antitumor role in U266 cells through the targeting of the MTDH/PI3K/Akt signaling pathway. Int J Oncol. 2018;53:2131–44.

    CAS  PubMed  Google Scholar 

  33. Zi Y, Zhang Y, Wu Y, Zhang L, Yang R, Huang Y. Downregulation of microRNA253p inhibits the proliferation and promotes the apoptosis of multiple myeloma cells via targeting the PTEN/PI3K/AKT signaling pathway. Int J Mol Med. 2021;47(3):8.

  34. He W, Fu Y, Zheng Y, Wang X, Liu B, Zeng J. Diallyl thiosulfinate enhanced the anti-cancer activity of dexamethasone in the side population cells of multiple myeloma by promoting miR-127-3p and deactivating the PI3K/AKT signaling pathway. BMC Cancer. 2021;21:125.

    Article  CAS  Google Scholar 

  35. Zhao Y, Xie Z, Lin J, Liu P. MiR-144-3p inhibits cell proliferation and induces apoptosis in multiple myeloma by targeting c-Met. Am J Transl Res. 2017;9:2437–46.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

No

Author information

Authors and Affiliations

Authors

Contributions

YY designed the study and drafted the manuscript. RD was responsible for the collection and analysis of the experimental data. RW revised the manuscript critically for important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Rui Wang.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of Lianshui People’s Hospital Affiliated to Kangda College of Nanjing Medical University. Signed written informed consents were obtained from the patients and/or guardians.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, Y., Ding, R. & Wang, R. Identification of candidate targets and mechanisms involved in miRNA regulation in multiple myeloma. World J Surg Onc 20, 23 (2022). https://doi.org/10.1186/s12957-021-02482-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12957-021-02482-1

Keywords