- Primary research
- Open Access
Identification of key biomarkers associated with cell adhesion in multiple myeloma by integrated bioinformatics analysis
Cancer Cell International volume 20, Article number: 262 (2020)
Multiple Myeloma (MM) is a hematologic malignant disease whose underlying molecular mechanism has not yet fully understood. Generally, cell adhesion plays an important role in MM progression. In our work, we intended to identify key genes involved in cell adhesion in MM.
First, we identified differentially expressed genes (DEGs) from the mRNA expression profiles of GSE6477 dataset using GEO2R with cut-off criterion of p < 0.05 and [logFC] ≥ 1. Then, GO and KEGG analysis were performed to explore the main function of DEGs. Moreover, we screened hub genes from the protein–protein interaction (PPI) network analysis and evaluated their prognostic and diagnostic values by the PrognoScan database and ROC curves. Additionally, a comprehensive analysis including clinical correlation analysis, GSEA and transcription factor (TF) prediction, pan-cancer analysis of candidate genes was performed using both clinical data and mRNA expression data.
First of all, 1383 DEGs were identified. Functional and pathway enrichment analysis suggested that many DEGs were enriched in cell adhesion. 180 overlapped genes were screened out between the DEGs and genes in GO terms of cell adhesion. Furthermore, 12 genes were identified as hub genes based on a PPI network analysis. ROC curve analysis demonstrated that ITGAM, ITGB2, ITGA5, ITGB5, CDH1, IL4, ITGA9, and LAMB1 were valuable biomarkers for the diagnosis of MM. Further study demonstrated that ITGA9 and LAMB1 revealed prognostic values and clinical correlation in MM patients. GSEA and transcription factor (TF) prediction suggested that MYC may bind to ITGA9 and repress its expression and HIF-1 may bind to LAMB1 to promote its expression in MM. Additionally, pan-cancer analysis showed abnormal expression and clinical outcome associations of LAMB1 and ITGA9 in multiple cancers.
In conclusion, ITGA9 and LAMB1 were identified as potent biomarkers associated with cell adhesion in MM.
MM is known as a plasma cell malignancy with an unlimited proliferation of abnormal plasma cells in the BM and high levels of monoclonal protein in the blood and urine . The advance of novel chemotherapeutic agents and cell therapy brings great improvements for MM patients. However, it remains largely incurable because of relapse and resistance [2, 3]. Therefore, searching for novel prognostic biomarkers and therapeutic targets is crucial for MM patients.
Previous studies have demonstrated the close interaction between myeloma cells and the BM microenvironment. The BM microenvironment included ECM components that composed of fibronectin, laminin, collagen and a variety of cellular residents that consisted of bone marrow stromal cells, hematopoietic stem cells, immune cells, etc. . MM cells and the BM microenvironment interacted through a complex network mediated by cytokines and adhesion molecules, including the integrins, the immunoglobulin superfamily, the selectins, the cadherins, and the proteoglycans. The cross-linking activated several known regulatory pathways, which involved in the survival, proliferation, migration, homing, as well as drug resistance of MM cells, playing a crucial role in MM development .
Integrins are heterodimeric membrane glycoproteins expressing on the surface of many types of cells, serving as CAMs and the major receptors for the ECM [5, 6]. Their expression can vary considerably between normal and cancer tissues. Studies showed that integrin expression levels were correlated with pathological outcomes, including patient survival and metastasis in many cancer types. Furthermore, Integrins regulated a series of cellular functions such as cell proliferation, migration, invasion, and survival which play a crucial role in the cancer progression [7, 8]. ITGA9 is one of the least studied integrins among the 24-member integrin family in human. Furthermore, studies have suggested that ITGA9 had the abnormal expression in numerous cancers and has been found to be crucial for a number of biological processes in many types of cancers, such as breast carcinoma, melanoma and lung cancer [9,10,11,12].
Laminins are extracellular heterotrimeric basement membrane glycoproteins . They are composed of three polypeptide chains, named as α, β, and γ . Laminins are involved in diverse physiological and pathological processes, including involvement in basement membrane assembly, neurite outgrowth, and promotion of cell adhesion, migration, protease activity, proliferation and angiogenesis in cancers [15, 16]. The link between cancer cells and laminins is vital in tumor invasion and metastasis. Invading tumor cells interact with laminins to acquire more metastatic potential . LAMB1 is expressed in most tissues and is one of the 3 chains constituting laminin 1. Studies have shown that LAMB1 was shown to be a potential biomarker in some cancers [18, 19].
However, there is rare research on the role which ITGA9 and LAMB1 play in myeloma development and metastasis. With the development of gene expression profiling array and second-generation high-throughput sequencing technology, people analyze the data on gene expression profiles to screen DEGs to find candidate biomarkers and potential therapeutic targets in MM [20, 21].In this study, we aim to find potential biomarkers in MM. First, we analyzed a microarray profile to obtain DEGs between active MM and normal plasma cells. Then, 12 genes were identified as hub genes based on a PPI network analysis. Finally, through integrative bioinformatic analysis, we identified that ITGA9 and LAMB1 were correlated with cell adhesion in MM.
Materials and methods
Microarray data information
Three gene expression profiles (GSE6477, GSE2658 and GSE136324) were from the GEO database (http://www.ncbi.nlm.nih.gov/geo). The GSE6477 dataset, which included samples from 147 patients with different stages of plasma cell neoplasm and 15 NDs, was performed on the Affymetrix Human Genome U133A Array platform. The GSE2658 dataset was composed of 559 samples of MM patients. 426 patients with MM was obtained from the GSE136324 dataset. Both of them were performed on the Affymetrix Human Genome U133 Plus 2.0 Array. All the above-mentioned profiling datasets were accessible online with no conflict.
Identification of DEGs in MM
The online analysis tool GEO2R was used to screen DEGs of the GSE6477 dataset. In this study, we set “p < 0.05 and [logFC] ≥ 1” as the cut-off criterion. A heatmap of the top 100 significantly changed DEGs were drawn using the Multiple Experiment Viewer software.
GO and KEGG enrichment analysis
GO enrichment analysis (http://www.geneontology.org/) and KEGG pathway enrichment analysis (http://www.genome.jp/kegg/pathway.html) were used to illuminate the biological functions and pathways associated with DEGs , which were both integrated in the DAVID (http://david.abcc.ncifcrf.gov/) program. Enriched GO and KEGG terms with p < 0.05 were considered as statistical significance. The ClueGO and CluePedia plug‐in Cytoscape software version 3.7.1 (http://www.cytoscape.org/) were used to analyze the pathways interrelation.
PPI network construction
STRING (http://string-db.org) online database is used to predict the PPI information. The cut-off value for the filtration criteria was set at the confidence score > 0.7. Then, the PPI network was visualized by Cytoscape. The MCODE plugin in Cytoscape was applied to select modules of the PPI network.
The ROC curves were used to explore the sensitivity and specificity of DEGs for MM diagnosis using GraphPad Prism Software (Version V, La Jolla, CA, USA).
The PrognoScan database (http://dna00.bio.kyutech.ac.jp/PrognoScan/) was used to investigate the prognostic impact of selected genes in different cancers. According to the median expression of a particular gene, the patients were split into high and low expression groups. The OS and DSS of GC patients was evaluated using a KM plot. A Cox p value<0.05 was considered as statistically significant.
Analysis of expression in various cancers
The mRNA expression levels of LAMB1 and ITGA9 in various cancers and their normal tissue counterparts were analyzed using the Oncomine database (https://www.oncomine.org/resource/login.html). The relationship between two specific gene mRNA expression levels in various cancers were analyzed using the GEPIA (Gene Expression Profiling Interactive Analysis) database (http://gepia.cancer-pku.cn/).
GSEA version 2.2.1 software was used to analyze genes function from the GSEA website MSIGDB database (http://software.broadinstitute.org/gsea/msigdb). The default weighted enrichment method was applied for enrichment analysis. The random combination was set for 1000 times. The analysis was performed with following settings: FDR<0.25, NOM p-value<0.05 and |NES|>1. The enrichment map plugin in Cytoscape was applied to visualize geneset enrichment results.
An independent sample t-test and one-way ANOVA were used to compare the statistical significance between two or more samples, respectively. Two-tailed p-values < 0.05 were considered as statistically significance. The correlations analysis was assessed by the Spearman correlation coefficient and Chi square test. Statistical analysis was performed by the SPSS 22.0 and GraphPad Prism Software.
Identification of DEGs
We acquired mRNA expression profiles of plasma cell samples from 101 patients of new or relapse multiple myeloma and 15 NDs from the GSE6477. Gene expression distribution of profiles was matched via boxplot analysis (Fig. 1a). Statistical analysis was required to validate the comparison between NDs and MM patients. We applied the GEO2R to identify DEGs associated with MM. 1383 DEGs (538 upregulated genes and 845 downregulated genes) were identified in the MM samples compared to NDs’ samples. We set the cut‐off criterion as p < 0.05 and [logFC] ≥ 1 (Fig. 1b). The top 100 genes were clustered in the heatmap between patients with MM patients and NDs (Fig. 1c).
Functional and pathway enrichment analysis
We performed GO and KEGG enrichment analysis to investigate the functions of DEGs using DAVID. The top GO (Fig. 2a–c, Table 1) and KEGG (Fig. 2d, Table 2) terms for DEGs were shown. For biological process (BP), DEGs were mainly enriched in immune response, inflammatory response, ECM organization, leukocyte migration, cell adhesion (Fig. 2a). DEGs in molecular function (MF) were significantly associated with protein binding, ECM structural constituent, serine type endopeptidase activity, protease binding and receptor activity (Fig. 2b). The cellular components (CC) analysis indicated that proteins encoded by DEGs were mostly located in the extracellular exosome, extracellular space, extracellular region and ECM (Fig. 2c). KEGG enrichment analysis showed that CAMs, proteoglycans in cancer, ECM-receptor interaction and PI3K/Akt signaling pathway were significantly enriched in DEGs. The top KEGG pathways for DEGs were shown (Fig. 2d). Subsequently, we analyzed a pathway interrelation and related genes by examining KEGG enrichment analysis results in ClueGO and CluePedia (Fig. 2e). These results suggested that cell adhesion might play a crucial part in the malignant progression of MM.
PPI network and modular analysis of selected genes
To identify key genes in the cell adhesion processes, 180 overlapped genes in DEGs and genes involved with the cell adhesion were screened, including 123 mRNAs and 57 non-coding RNAs (Fig. 3a). Then, we constructed a PPI network containing 123 protein, with 123 nodes and 429 edges based on the STRING database and Cytoscape (Fig. 3c). 12 hub genes whose degree values ≥ 15 for further analysis were chosen, including ITGB1, FN1, ITGB3, ITGAM, PTPRC, ITGB2, ITGA5, ITGB5, CDH1, IL4, ITGA9 and LAMB1 (Fig. 3b). Additionally, Cytoscape displayed a total of eight modules in the default MCODE settings for modular analysis (Fig. 3d–k, Table 3).
Diagnostic and prognostic values of hub genes in MM
First of all, we performed the ROC curve analysis among 12 hub genes based on the GSE2658. The results showed that ITGAM, ITGB2, ITGA5, ITGB5, CDH1, IL4, ITGA9, and LAMB1 achieved an AUC value of > 0.7, demonstrating that these eight genes have high sensitivity and specificity for MM, suggesting they can be served as biomarkers for the diagnosis of MM (Fig. 4a–h).
Furthermore, we evaluated the influence of all aforementioned 8 hub genes on clinical prognosis using PrognoScan database based on the GSE2658 through K-M curve and log-rank test to identify whether these genes were concerned with the survival of MM patients. ITGA9 (Fig. 4i) and LAMB1 (Fig. 4j) were significantly associated with DSS in MM by evaluating the correlation between these gene expressions and survival rates.
The relationship between ITGA9/LAMB1 mRNA expression and clinical characters of patients with MM
First of all, the relationship between ITGA9 or LAMB1 mRNA expression and the clinical parameters was analyzed based on both clinical data and mRNA expression data of GSE136324 (Table 4, n = 426). The results showed that LAMB1 expression was negatively correlated with β2-MG (Fig. 5a, p = 0.024, spearman correlation coefficient r = -0.110). Patients with low LAMB1 expression had a tendency of being in high R-ISS stage (Fig. 5c, p = 0.001). No significant relationship was found between LAMB1 mRNA expression and malignant plasma cells (PCs) in BM (Fig. 5b, p = 0.9427) and GEP groups (Fig. 5d, p = 0.1785).
Furthermore, there was a negative correlation between ITGA9 expression and malignant PCs in BM (Fig. 5f, p < 0.001, spearman correlation coefficient r = −0.165). Compared with in the standard group, MM patients in the GEP high group had a lower ITGA9 expression (Fig. 5h, p = 0.0039). ITGA9 mRNA expression had no correlation with β-2MG (Fig. 5e, p = 0.323) and R-ISS stage (Fig. 5g, p = 0.281).
Function and signaling pathways analysis of ITGA9 and LAMB1
To investigate the function of ITGA9 and LAMB1 on MM progression, we performed a comprehensive analysis including GSEA, TF prediction and Spearman correlation analysis.
Firstly, an enrichment map was constructed using genesets which related to cell adhesion with p‐value < 0.05 (Fig. 6a). Genesets including GO calcium dependent cell cell adhesion via plasma membrane cell adhesion molecules, GO cell adhesion via plasma membrane adhesion molecules, KEGG ECM receptor interaction, NABA basement membranes, NABA core matrisome were enriched in patients with low LAMB1 expression in the GSE2658, suggesting that LAMB1 was involved in the cell adhesion in myeloma cells (Fig. 6b–f). Furthermore, cancer microenvironment- dn, MANALO hypoxia dn, ELVIDGE hypoxia dn genesets were significantly enriched in patients with low LAMB1 expression in the GSE2658 (Fig. 6g–i). For futher study, we found the positive correlation between LAMB1 and HIF-1 mRNA expression in a varity of cancers using GEPIA database (Fig. 6j, p = 0.00018, spearman correlation coefficient r = 0.038). The result of TF prediction by the online site PROMO (http://alggen.lsi.upc.es/cgi-bin/promo_v3/promo/promoinit.cgi?dirDB=TF_8.3/) strengthened the correlation that HIF-1 may bind to the promoter of LAMB1 (Fig. 6k). Therefore, LAMB1 which may be regulated by HIF-1 played a vital role in myeloma cell adhesion.
ITGA9 also was a cricial cell adhesion molecule since GSEA results showed that genesets including KEGG ECM receptor interaction, GO cell adhesion via plasma membrane adhesion molecules, NABA ECM glycoproteins, GO calcium dependent cell cell adhesion via plasma membrane cell adhesion molecules, GO protein complex involved in cell adhesion, NABA core matrisome, cell–cell adhesion were enriched in patients with low ITGA9 expression in the GSE2658 (Fig. 7b–h). An enrichment map was constructed using genesets which related to cell adhesion with p‐value < 0.05 (Fig. 7a). Results demonstrated that SCHLOSSER myc targets repressed by serum and hallmark myc targets v1 genesets were enriched in patients with low ITGA9 expression in the GSE2658 (Fig. 7i–j). There was a negative correlation between the mRNA expression of ITGA9 and MYC in multiple cancers based on the GEPIA database (Fig. 7k, p = 7.7e-19, spearman correlation coefficient r = −0.09). We predicted that MYC may bind to the promoter of ITGA9 and repressed its expression based on TF prediction by PROMO (Fig. 7l). In conclusion, ITGA9 was a vital cell adhesion molecule in myeloma which may be negtive regulated by MYC.
ITGA9 and LAMB1 expression and clinical outcome association in different cancers
We explored the expression and prognosis of ITGA9 and LAMB1 mRNA in different cancers based on the Oncomine database and the PrognoScan database respectively.
LAMB1 is highly expressed in most cancers while low in some cancers including breast cancer, leukemia, ovarian cancer and prostate cancer (Fig. 8a). Furthermore, LAMB1 was significantly associated with OS and DSS respectively in AML, DLCBL, glioma, NSCLC, colorectal cancer (Fig. 8b–g) and colorectal cancer, breast cancer (Fig. 8h–i).
Compared with ND, the transcript levels of ITGA9 indicated significant low expression in bladder cancer, brain and central nervous system cancer, breast cancer, leukemia, liver cancer, lung cancer and etc., suggesting that the down-regulation of ITGA9 was common in various types of cancer (Fig. 9a). ITGA9 was significantly associated with OS and DSS respectively in AML, B cell lymphoma, lung cancer, colorectal cancer and esophagus cancer (Fig. 9b–f) and colorectal cancer, breast cancer (Fig. 9g–h). The details including hazard ratio (HR) with 95% confidence intervals(CI) and p-values were shown in the Table 5.
MM is a common hematological malignancy so it is vital to investigate the molecular mechanisms. Microarray has been widely used to analyze the expression changes of genes in MM and predict the potential biomarkers.
In this study, we analyzed gene mRNA expression data from BM plasma cells in GSE6477 dataset using GEO2R and revealed that there were 1383 DEGs between MM patients and NDs, consisting of 538 upregulated genes and 845 downregulated genes with p < 0.05 and [logFC] ≥ 1. Furthermore, we performed GO and KEGG enrichment analyses to explore main function of DEGs. GO enrichment analysis in BPs recognized the functional enrichment of DEGs in the immune response, inflammatory response, ECM organization, leukocyte migration, and cell adhesion. KEGG analysis showed enrichment of hematopoietic cell lineage, complement and coagulation cascades, CAMs, proteoglycans in cancer, ECM-receptor interaction and PI3K/Akt signaling pathway. It was quite clear that the loss of cell–cell adhesion was an important event for acquiring the invasive and metastatic phenotype. These results indicated the importance of cell adhesion in MM progression.
We screened overlapped genes between DEGs and genes in GO terms of cell adhesion and constructed a PPI network. The 12 hub genes were listed: ITGB1, FN1, ITGB3, ITGAM, PTPRC, ITGB2, ITGA5, ITGB5, CDH1, IL4, ITGA9, and LAMB1. Among the 12 hub genes, ITGAM, ITGB2, ITGA5, ITGB5, CDH1, IL4, ITGA9 and LAMB1 had diagnostic value, which can distinguish MM from normal people. The mRNA expression of ITGA9 and LAMB1 were significantly associated with the OS and DSS of MM patients. These results suggest that ITGA9 and LAMB1 exhibit both diagnostic and prognostic values in MM. Further investigations regarding to the role of LAMB1 and ITGA9 were required.
Here, we found that LAMB1 was abnormally expressed and was associated with OS and DSS in many cancers, which were consistent with previous results. LAMB1 has a high protein level in high-grade gliomas, suggesting a possible correlation with tumor progression . What’s more, LAMB1 was identified to take part in cell attachment and have the capacity to inhibit metastasis. In prostate cancer, LAMB1 was shown to be involved in cell motility and invasion into the surrounding ECM . In our study, low LAMB1 expression was significantly associated with high β2-MG concentration and high R-ISS stage in MM. GSEA results showed that LAMB1 was involoved in cell adhesion and may be induced by hypoxia. TF prediction and the correlation analysis suggested that HIF-1 may bind to LAMB1 promoter to increase its transcription.
ITGA9 abnormal expression was found in many cancers and was likely to correlate with higher grade cancers . For example, sequencing analysis also found that ITGA9 was significantly down-regulated in cervical squamous cell carcinoma . Genetic variation and epigenetic modification of ITGA9 are related with the tumorgenicity and progression of colorectal cancer . In our study, we found down-regulated ITGA9 was associated with poor outcome in MM. Myeloma patients with low ITGA9 expression had more tendency of having the higher number of malignant PCs in BM and becoming the higher GEP group. GSEA results showed that ITGA9 was a vital cell adhesion molecule and may be repressed by MYC in myeloma. Based on the TF prediction and the correlation analysis, we predicted MYC may bind to the promoter of ITGA9 for transcriptional repression.
We identified eight hub genes, including ITGAM, ITGB2, ITGA5, ITGB5, CDH1, IL4, ITGA9, and LAMB1 to be potential diagnostic markers in MM. Further study demonstrated ITGA9 and LAMB1 which correlated with clinical characters and prognosis may play important roles in the cell adhesion and can be regulated by different TFs in MM.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.
Area under the curve
Cell adhesion molecules
Database for annotation, visualization and integrated discovery
Differentially expressed gene
False discovery rate
Gene expression omnibus
Gene expression profile
Gene expression profiling interactive analysis
Gene enrichment analysis
Integrin subunit alpha
Integrin subunit beta
Kyoto encyclopedia of genes and genomes
Laminin subunit beta 1
Molecular complex detection
Monoclonal gammopathy of undetermined significance
Normalized enrichment score
Protein tyrosine phosphatase receptor type c
Revised-international staging system
Receiver operatoring characteristic
Smoldering multiple myeloma
Search tool for the retrieval of interacting gene
Chng WJ, et al. Targeted therapy in multiple myeloma. Cancer Control. 2005;12(2):91–104.
Robak P, et al. Drug resistance in multiple myeloma. Cancer Treat Rev. 2018;70:199–208.
Leleu X. New hope for relapsed and refractory multiple myeloma. Lancet Oncol. 2013;14(11):1028–9.
Neri P, Bahlis NJ. Targeting of adhesion molecules as a therapeutic strategy in multiple myeloma. Curr Cancer Drug Targets. 2012;12(7):776–96.
Katz BZ. Adhesion molecules–The lifelines of multiple myeloma cells. Semin Cancer Biol. 2010;20(3):186–95.
Hynes RO. Integrins: versatility, modulation, and signaling in cell adhesion. Cell. 1992;69(1):11–25.
Seguin L, et al. Integrins and cancer: regulators of cancer stemness, metastasis, and drug resistance. Trends Cell Biol. 2015;25(4):234–40.
Desgrosellier JS, Cheresh DA. Integrins in cancer: biological implications and therapeutic opportunities. Nat Rev Cancer. 2010;10(1):9–22.
Hoye AM, et al. The newcomer in the integrin family: integrin alpha9 in biology and cancer. Adv Biol Regul. 2012;52(2):326–39.
Allen MD, et al. Clinical and functional significance of alpha9beta1 integrin expression in breast cancer: a novel cell-surface marker of the basal phenotype that promotes tumour cell invasion. J Pathol. 2011;223(5):646–58.
Zhang J, et al. MicroRNA-125b suppresses the epithelial-mesenchymal transition and cell invasion by targeting ITGA9 in melanoma. Tumour Biol. 2016;37(5):5941–9.
Hibi K, et al. Aberrant upregulation of a novel integrin alpha subunit gene at 3p21.3 in small cell lung cancer. Oncogene. 1994;9(2):611–9.
Timpl R, Brown JC. The laminins. Matrix Biol. 1994;14(4):275–81.
Burgeson RE, et al. A new nomenclature for the laminins. Matrix Biol. 1994;14(3):209–11.
Ryan MC, et al. The functions of laminins: lessons from in vivo studies. Matrix Biol. 1996;15(6):369–81.
Givant-Horwitz V, Davidson B, Reich R. Laminin-induced signaling in tumor cells. Cancer Lett. 2005;223(1):1–10.
Engbring JA, Kleinman HK. The basement membrane matrix in malignancy. J Pathol. 2003;200(4):465–70.
Lin Q, et al. Analysis of colorectal cancer glyco-secretome identifies laminin beta-1 (LAMB1) as a potential serological biomarker for colorectal cancer. Proteomics. 2015;15(22):3905–20.
Alinezhad S, et al. Validation of novel biomarkers for prostate cancer progression by the combination of bioinformatics, clinical and functional studies. PLoS ONE. 2016;11(5):e0155901.
Angarica VE, Del Sol A. Bioinformatics tools for genome-wide epigenetic research. Adv Exp Med Biol. 2017;978:489–512.
Kulasingam V, Diamandis EP. Strategies for discovering novel cancer biomarkers through utilization of emerging technologies. Nat Clin Pract Oncol. 2008;5(10):588–99.
Senchenko VN, et al. Novel tumor suppressor candidates on chromosome 3 revealed by NotI-microarrays in cervical cancer. Epigenetics. 2013;8(4):409–20.
Gerashchenko GV, et al. Screening of epigenetic and genetic disturbances of human chromosome 3 genes in colorectal cancer. Ukr Biokhim Zh (1999). 2009;81(4):81–7.
Chen Q, et al. MiR-124-5p inhibits the growth of high-grade gliomas through posttranscriptional regulation of LAMB1. Neuro Oncol. 2014;16(5):637–51.
Ethics approval and consent to participate
Consent for publication
This is not applicable for this review.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Peng, Y., Wu, D., Li, F. et al. Identification of key biomarkers associated with cell adhesion in multiple myeloma by integrated bioinformatics analysis. Cancer Cell Int 20, 262 (2020). https://doi.org/10.1186/s12935-020-01355-z
- Bioinformatics analysis
- Multiple myeloma