Identification of Key Genes and Pathways in Myeloma side population cells by Bioinformatics Analysis

Background: Multiple myeloma (MM) is the second most common hematological malignancy, which is still incurable and relapses inevitably, highlighting further understanding of the possible mechanisms. Side population (SP) cells are a group of enriched progenitor cells showing stem-like phenotypes with a distinct low-staining pattern with Hoechst 33342. Compared to main population (MP) cells, the underlying molecular characteristics of SP cells remain largely unclear. This bioinformatics analysis aimed to identify key genes and pathways in myeloma SP cells to provide novel biomarkers, predict MM prognosis and advance potential therapeutic targets. Methods: The gene expression profile GSE109651 was obtained from Gene Expression Omnibus database, and then differentially expressed genes (DEGs) with P-value <0.05 and |log2 fold-change (FC)| > 2 were selected by the comparison of myeloma light-chain (LC) restricted SP (LC/SP) cells and MP CD138+ cells. Subsequently, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis, protein-protein interaction (PPI) network analysis were performed to identify the functional enrichment analysis of the DEGs and screen hub genes. Cox proportional hazards regression was used to select the potential prognostic DEGs in training dataset (GSE2658). The prognostic value of the potential prognostic genes was evaluated by Kaplan-Meier curve and validated in another external dataset (MMRF-CoMMpass cohort from TCGA). Results: Altogether, 403 up-regulated and 393 down-regulated DEGs were identified. GO analysis showed that the up-regulated DEGs were significantly enriched in innate immune response, inflammatory response, plasma membrane and integral component of membrane, while the down-regulated DEGs were mainly involved in protoporphyrinogen IX and heme biosynthetic process, hemoglobin complex and erythrocyte differentiation. KEGG pathway analysis suggested that the DEGs were significantly enriched in osteoclast differentiation, porphyrin and chlorophyll metabolism and cytokine-cytokine receptor interaction. The top 10 hub genes, identified by the plug-in cytoHubba of the Cytoscape software using maximal clique centrality (MCC) algorithm, were ITGAM, MMP9, ITGB2, FPR2, C3AR1, CXCL1, CYBB, LILRB2, HP and FCER1G. Modules and corresponding GO enrichment analysis indicated that myeloma LC/SP cells were significantly associated with immune system, immune response and cell cycle. The predictive value of the prognostic model including TFF3, EPDR1, MACROD1, ARHGEF12, AMMECR1, NFATC2, HES6, PLEK2 and SNCA was identified, and validated in another external dataset (MMRF-CoMMpass cohort from TCGA). Conclusions: In conclusion, this study provides reliable molecular biomarkers for screening, prognosis, as well as novel therapeutic targets for myeloma LC/SP cells.


Introduction
Multiple myeloma (MM) is a B-cell malignancy characterized by the aberrant expansion of clonal plasma cells within bone marrow, which is the second most common hematological malignancy [1]. Despite Ivyspring International Publisher remarkable progress of biology and recent development of novel therapy [2,3], MM continues to remain incurable due to the emergence of drug resistance and frequent relapses, highlighting the further understanding of the possible mechanisms.
Cancer stem cells (CSCs) are thought to have the distinctive properties of constituting a small fraction of tumor cells with self-renewal capacity and be able to propagate the disease [4,5]. Besides, CSCs are considered to be more resistant to chemo-and radiotherapy and have better DNA repair mechanisms and increased anti-apoptotic activity, just like hematopoietic stem cells [6]. Previously, CSCs have been identified in MM [7] by the evidence that the CD138 -/CD19 + fraction of MM has a greater clonogenic potential and the phenotype of a memory B-cell (CD19 + , CD27 + ), resulting in the development of refractory clones and disease relapse [8]. Then, it has been defined that possible stem cell populations include light-chain restricted cells with a CD138 -/CD19 + /CD27 + phenotype [7,9,10], CD138 + /CD34 + / B7 -H1 + subpopulations [11] and CD38 ++ /CD45plasma cells [12,13]. Despite these phenotypes have been described, the distinct CSCs marker in MM is still controversial.
Side population (SP) cells, first described by Goodell et al. [14], are a group of enriched progenitor cells showing stem-like phenotypes and a distinct low-staining pattern with Hoechst 33342, and have been widely used as a unique source for studying CSCs in the absence of specific markers [4,[15][16][17][18][19][20][21]. Although lots of previous studies have explored the stem-like properties and tumorigenicity of myeloma SP cells, a better understanding of SP cells still remains largely unclear [10,[22][23][24][25]. Thus, it is vital to elucidate the key molecular characteristics expressed within myeloma SP cells.
It is generally known that gene expression profiling analysis based on microarray technology enables the possibilities for identifying certain disease-related biomarkers. Recently, many studies have been carried out on the base of microarray data profiles to identify the pathogenesis of MM [26][27][28]. Nevertheless, the key molecular characteristics of myeloma SP cells in comparison to MP cells have not yet been explored. This bioinformatics analysis was performed to elucidate key candidate genes and pathways in myeloma SP cells, provide novel biomarkers, predict MM prognosis and advance potential therapeutic targets.
In this study, we downloaded microarray dataset GSE109651 (Zhan et al., 2018) from Gene Expression Omnibus database (https://www.ncbi. nlm.nih.gov/geo/), which is a public functional genomics data repository with array-and sequence-based data. By comparing myeloma light-chain (LC) restricted SP (LC/SP) cells with myeloma MP cells based on R software and Bioconductor, differentially expressed genes (DEGs) were identified. Gene Ontology (GO) analysis, Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis and protein-protein interaction (PPI) network analysis were performed to identify the functional enrichment analysis of the DEGs and screen hub genes. Subsequently, we constructed a prognostic model to predict survivals of MM patients. This study provides reliable molecular biomarkers for screening, prognosis, as well as novel therapeutic targets for LC/SP cells of MM.

Microarray data profile
The gene expression dataset GSE109651 was obtained from GEO database. The microarray data of GSE109651, based on the GPL570 platform ([HG-U133_Plus_2] Affymetrix Human Genome U133A Plus 2.0 Array) and normalized using the MAS5 algorithm of the Affymetrix expression console version1.1 software (Affymetrix), includes 7-paired LC/SP cells and MP CD138 + cells of myeloma bone marrow from 7 diagnosed MM patients isolated by fluorescence-activated cell sorting (FACS) using Hoechst 33342 and CD138 antibody. To perform survival analysis, GSE2658 dataset of 559 MM patients and TCGA MM RNA sequencing dataset (MMRF-CoMMpass) of 787 cases with MM including clinicopathological information were downloaded from GEO database and TCGA (https://tcga-data. nci.nih.gov/) databases, respectively. For the retrospective cohort, the patients' characteristics were estimated by Pearson test χ 2 or Fisher's exact test, indicating no significant statistical difference.

Data processing and identification of DEGs of GSE109651
Firstly, we detected the quality of raw data by R statistical software (version 3.6.3, https://www.rproject.org/), including a quality control overview diagram based on the "simpleaffy" package, weights and residuals plot, relative log expression (RLE) boxplot and normalized unscaled standard errors (NUSE) box plot based on the "affyPLM" and "RColorBrewer" packages, RNA degradation curve based on the "affy" package and clustering analysis diagram based on the "gcrma", "graph" and "affycoretools" packages.
Then, DEGs between LC/SP cells and MP CD138 + cells of MM were identified by an empirical Bayes method based on the "limma" package in R.
The process included six main steps: construction of a gene expression matrix, construction of an experimental design matrix, construction of a contrast matrix, fitting of a linear model, Bayes test, and generation of results. In this study, genes with P-value < 0.05 and |log2 fold-change (FC)| > 2 were defined as DEGs.

GO and KEGG pathway enrichment analysis of DEGs
To explore the functional roles of the above DEGs, DAVID database (https://david.ncifcrf.gov/) was used to perform GO term enrichment analysis of molecular function (MF), biological process (BP), and cellular component (CC) and KEGG pathway enrichment analysis. P-value < 0.05 was considered as the cut-off criterion.

PPI network construction and modular analysis
Search Tool for the Retrieval of Interacting Genes (STRING) database (http://www.string-db.org/) was used to construct the PPI network. The visualization and analysis of the PPI network were based on Cytoscape software version 3.7.2. Then, the plug-ins Molecular Complex Detection (MCODE) and Biological Network Gene Ontology tool (BiNGO) in Cytoscape software were used to screen significant modules of the PPI network (the parameters were set to default) and perform GO analysis that the module genes were significantly enriched in.

Identification of hub genes
The plug-in cytoHubba in Cytoscape was used to identify key (hub) genes among the above DEGs by maximal clique centrality (MCC) computing method. The hub genes were selected to discuss their function and effect on myeloma LC/SP cells.

Survival analysis
DEGs significantly associated with myelomaspecific survival in the training dataset (GSE2658) were identified using univariate Cox proportional hazards analysis with P-value < 0.01 by "survival" package [29]. Then, the final genes significantly correlated with survival at a P-value of less than 0.05 were identified by multivariate Cox proportional hazards analysis. Subsequently, the risk score on the base of the aforementioned candidate genes and survival information was calculated as follows: Risk score = ∑β i × ExpGene i (β i was the coefficient value and ExpGene i was the gene expression level). According to the median risk score, the cohort was dichotomized into low-risk and high-risk group, then survival time was compared by the Kaplan-Meier analysis and the log-rank test with a P-value of less than 0.01. Another external dataset (MMRF-CoMMpass cohort) was used to assess the prognostic value through a process similar to the training dataset.

Identification of DEGs
The gene expression dataset GSE109651 included 7-paired LC/SP cells samples and MP CD138 + cells samples of myeloma bone marrow. On the basis of cut-off criterion of DEGs described previously, there were 796 DEGs in LC/SP cells compared with MP CD138 + cells of myeloma bone marrow, among which 393 DEGs were significantly down-regulated and 403 DGEs were significantly up-regulated. The volcano plot of DEGs was shown in Figure 1. The expression heat map of the top 100 DEGs (including 52 significantly down-regulated genes and 48 significantly up-regulated genes) was depicted in Figure 2, which could effectively distinguish LC/SP cells from MP CD138 + cells and might function as biomarker and target of MM. The detailed information of the top 10 DEGs was shown in Table 1.

GO term enrichment analysis of DEGs
To explore the functional roles of the DEGs, we performed GO enrichment analysis of up-regulated and down-regulated DEGs by using DAVID gene annotation tool. It turned out that obvious differences were enriched in BPs, MFs and CCs among the 796 DEGs. For BPs, the up-regulated DEGs were primarily enriched in immune response, including innate, adaptive immune response and T cell differentiation involved in immune response, suggesting that these DEGs could significantly associate with the immune system of myeloma LC/SP cells. Besides, these genes were also significantly enriched in inflammatory response, leukotriene metabolic process and neutrophil chemotaxis. The down-regulated DEGs were significantly enriched in protoporphyrinogen IX biosynthetic process, heme biosynthetic process and erythrocyte differentiation, indicating that the down-regulated DEGs may be relevant to the development and differentiation of erythrocytes. In the CCs group, the up-regulated DEGs were significantly involved in plasma membrane, integral component of membrane and extracellular space. In addition, down-regulated genes were largely enriched in the extracellular exosome and hemoglobin complex. Regarding MFs category, the up-regulated genes were mainly enriched in the binding of carbohydrate, calcium ion and arachidonic acid. Moreover, the most significantly enriched GO terms for down-regulated genes were immunoglobulin receptor binding, oxygen transporter activity and NAD activity. GO enrichment analysis results were displayed in Figure  3 and Table 2.

KEGG pathway enrichment analysis of DEGs
According to the KEGG pathway enrichment analysis of up-and down-regulated DEGs, the upregulated DEGs were mainly enriched in osteoclast differentiation, cytokine-cytokine receptor interaction, Staphylococcus aureus infection, leukocyte transendothelial migration and cell adhesion molecules. Furthermore, enrichment of downregulated DEGs was mostly in the porphyrin and chlorophyll metabolism, hematopoietic cell lineage and metabolic pathways. KEGG analysis results were displayed in Figure 4 and Figure 5, and the detailed analysis results of the top 5 pathways were shown in Table 3.

PPI network construction and modular analysis
On the base of STRING online database and Cytoscape software, we established a PPI network of these DEGs in myeloma LC/SP cells, with 610 nodes and 2922 edges identified, including 288 up-regulated and 322 down-regulated genes. Then, PPI module analysis was implemented by plug-ins MCODE in Cytoscape, and three significant modules were identified from the whole network. The top 3 modules with high scores were selected for display: module 1 contained 53 nodes and 684 edges ( Figure 6A), module 2 contained 21 nodes and 197 edges ( Figure  6B) and module 3 contained 13 nodes and 60 edges ( Figure 6C). Subsequently, corresponding GO term enrichment analysis was performed by plug-ins BiNGO in Cytoscape. Genes in Module 1 were significantly enriched in the defense response, immune system process and immune response. Moreover, genes in Module 2 were mainly enriched in cell cycle phase, cell cycle and M phase of mitotic cell cycle. Additionally, genes in module 3 were primarily enriched in G-protein coupled receptor protein, signaling pathway and chemotaxis. The detailed information of the top 3 modules was shown in Table  4.

Selection of hub genes from the PPI network
Among the previously described DEGs, significant hub genes were identified by plug-in cytoHubba of Cytoscape using MCC algorithm. The top 10 hub genes were ITGAM, MMP9, ITGB2, FPR2, C3AR1, CXCL1, CYBB, LILRB2, HP and FCER1G. The top 10 hub genes and their most relevant functions were displayed in Table 5.

Survival analysis of DEGs
The result of univariate Cox analysis showed 76 survival related genes (Table S1; 19 up-regulated and 57 down-regulated) (P-value < 0.01).
Afterwards, the 76 genes were fitted into the multivariate Cox proportional hazards analysis, and 9 genes including TFF3, EPDR1, MACROD1, ARHGEF12, AMMECR1, NFATC2, HES6, PLEK2 and SNCA were identified with P-value < 0.05 (Table 6). The prognostic models of training and validation dataset containing 9 genes were constructed by discriminating the low-risk group from the high-risk group based on the respective median risk score. Kaplan-Meier curve showed that high-risk group had worse survival compared to the low-risk group in both training and validation dataset (Figure 7).

Discussion
In the present study, we performed a bioinformatics analysis to identify DEGs between myeloma LC/SP cells and MP CD138 + cells to explore the molecular characteristics of LC/SP cells. Based on the gene expression profile, we screened a total of 796 DEGs, including 403 up-regulated and 393 down-regulated genes. Subsequently, deeper exploration of these DEGs were performed by bioinformatics methods, including GO and KEGG pathway enrichment analysis, PPI network construction and modules analysis, selection of hub genes and survival analysis. GO enrichment analysis demonstrated that the up-regulated DGEs were significantly enriched in innate immune response, inflammatory response, plasma membrane and integral component of membrane. Firstly, according to our enrichment analysis, up-regulated DEGs were most enriched in innate immune response. As reported by Grivennikov SI et al., components of innate immunity such as macrophages, and DCs can either induce anti-tumor immune responses or promote tumor growth and progression depending on their morphological and phenotypic subtypes [30]. In multiple solid tumor models, the presence of tumor infiltrating macrophages (TAM) in tumor lesions can promote "stemness" property of cancer cells [31]. However, for MM, the association between SP cells and innate immune response has not been explained clearly yet. With regard to up-regulated DEGs enriched in immune response, it is now well-established that FGR, CXCL1, NLRC4 and S100A9 influence the pathogenesis of cancer by modulating immune responses and promoting progression, aggressiveness and cell survival [32][33][34][35]. Besides, up-regulated DGEs were also significantly enriched in inflammatory response. Then, the up-regulated DEGs were significantly involved in plasma membrane and integral component of membrane. It has been reported that the unique and specific makeup and arrangement of cell membranes of cancer cells are critical for cells to survive, grow and proliferate [36]. The enrichment analysis indicated that myeloma LC/SP cells may have unique plasma membrane and integral component of membrane compared to MP cells, and targeting the uniqueness may lead to the reduction of SP cells. Additionally, the downregulated DEGs were significantly enriched in protoporphyrinogen IX and heme biosynthetic process, hemoglobin complex and erythrocyte differentiation, indicating that the down-regulated DEGs may be relevant to the development and differentiation of erythrocytes. Moreover, extracellular exosome was also enriched significantly. Exosomes are membranous structures that carry signaling molecules and regarded as important mediators of inter-cellular communication in health and disease [37]. Studies have revealed a strong cross-talk between the MM cells and their microenvironment in the bone marrow, which leads to the final phenotype of a typical MM patient [38,39]. This result demonstrated that extracellular exosome may function significantly in this small fraction of MM cells.
KEGG pathway enrichment analysis showed some DEGs were significantly enriched in osteoclast differentiation, porphyrin and chlorophyll metabolism and cytokine-cytokine receptor interaction. Osteolytic bone disease is the hallmark of MM, which deteriorates the quality of life of myeloma patients. It has been demonstrated that increased osteoclast activity is one of the important mechanisms [40]. Among the DEGs, some studies found that PIK3CG, LILRB2 and CYBB could regulate the differentiation of osteoclast, which highlighted the possible biological significance of LC/SP cells in osteoclast differentiation [41][42][43]. ALAS2, significantly enriched in porphyrin and chlorophyll metabolism pathway, plays a key role in erythropoiesis by regulation of erythroid heme synthesis [44]. As already described in GO analysis, parts of the down-regulated DEGs may associate with the erythropoiesis, indicating that myeloma LC/SP cells may impair the erythropoiesis. CXCL1, significantly enriched in cytokine-cytokine receptor interaction pathway, could result in the enhancement of MM cell viability and migration [45]. Additionally, Staphylococcus aureus infection was enriched in KEGG pathway analysis. It's well known that infectious complications are a frequent cause of morbidity and mortality of MM [46]. A prospective study observed the rate of infections varied in different phases of MM, and the most infections were clinically diagnosed as pneumonia and bronchopneumonia caused by Haemophilus influenzae or Streptococcus pneumonia in early-stage MM [47]. As to Staphylococcus aureus infection, a recent study showed that Staphylococcus aureus bacteremia (SAB) may be an early prognostic indicator of cancer because of the phenomenon that patients with SAB were more likely to die from cancer than the general population [48]. Furthermore, an association between SAB and risk of multiple myeloma was described [49]. According to our results, it is presumable that myeloma SP cells may be relevant to infectious complications, especially the SAB infection, which opens a fundamental direction to understand infections for patients suffering from MM. However, the results need to be confirmed in further basic and clinical research. In brief, the enriched GO and KEGG pathways partly clarified the specific molecular characteristics of myeloma SP cells.
Then a PPI network of these DEGs in myeloma LC/SP cells was established, containing 610 nodes and 2922 edges. In the network, we selected three significant modules through the degree of importance and corresponding GO term enrichment analysis was performed, which indicated that myeloma LC/SP cells was significantly associated with immune system process, immune response and cell cycle, basically consistent with what we mentioned above. Subsequently, ten significant hub genes have been identified, containing ITGAM, MMP9, ITGB2, FPR2, C3AR1, CXCL1, CYBB, LILRB2, HP and FCER1G.
ITGAM was identified as the top 1 hub gene and had the highest degree of connectivity.
ITGAM encodes CD11b, a component of the macrophage-1 antigen complex (Mac1, also known as complement receptor 3 [CR3]), which together with CD18, form Mac-1 or CR3, a protein that mediates leukocyte adhesion, migration, and phagocytosis in different cells [50][51][52][53]. CD11b contributes to the phagocytosis of opsonized particles, including apoptotic cells and immune complex [53]. What's more, CD11b is defined as a marker for myeloid-derived suppressor cells, which is reported to be harnessed by malignant cells to restrain antitumor immunity and promote malignant expansion or refractoriness to treatment [54][55][56]. It has been considered as a poor prognostic factor in MM [57] and AML patients [58][59][60][61]. But association between SP cells and CD11b remains unclear, and it is presumable that CD11b may participate in the regulation of biology of LC/SP cells and its up-regulation may promote expansion of MM. ITGB2 produces a protein, known as CD18, which is a cell surface marker expressed on lymphocytes [62] and is involved in cell adhesion and cell-surface mediated signaling [63]. It has been demonstrated that mutation in the ITGB2 gene could lead to leukocyte adhesion deficiency [64]. And its expression in CLL cells predicts disease progression [65]. In MM cell line, ITGB2 is overexpressed in vincristine resistant cell line [66]. Nonetheless, the correlation of drug resistance and ITGB2 requires further analysis.
As a key adhesion receptor, integrin CD11b/ CD18 meditates leukocyte migration and immune functions [67]. Recently, several studies have investigated that the adhesion and angiogenesis system is vital to propagate MM progression with a vicious cycle by the endothelial-MM interaction. β integrin has been described to participate in the homing and adhesion of endothelial progenitor cells to sites of vascular remodeling [68,69]. It has been uncovered that some integrins were detected in high levels in MM, while in non-detectable levels in nonactive MM and MGUS patients, suggesting the adhesion molecules support the interactions between MM and the microvasculature and facilitate disease progression [70]. Furthermore, junctional adhesion molecule A has been identified as a key mediator of MM progression by promoting MM-associated angiogenesis and an independent prognostic factor for both newly diagnosed MM and relapsed/ refractory MM [71,72]. Similarly, our enrichment analysis of up-regulated DEGs had identified positive regulation of angiogenesis and cell adhesion, demonstrating that myeloma SP cells may be relevant to angiogenesis and cell adhesion to propagate MM progression.
The significantly up-regulated MMP9 gene (matrix metallopeptidase 9), one of the most widely investigated matrix metalloproteinases, is a significant protease which plays vital roles in many biological processes and cancer cell invasion, metastasis and angiogenesis [73]. Recently, MMP9 has been identified as a potential biomarker for several cancers [74][75][76][77][78]. As far as MM concerned, previous reports indicated that the expression of MMP9 in MM cells promote MM invasion [79][80][81][82][83][84], which may highlight the role of increased neovascularization in MM progression. In fact, it has been testified that angiogenesis, which is linked to aberrant expression of pro-angiogenic and down-regulation of antiangiogenic genes [85], is a feature of MM progression through the transition from MGUS to MM, and plays a role in medullary and extramedullary dissemination [86,87]. Recently, several angiogenic factors in active MM have been discovered, like VEGF (Vascular endothelial growth factor), FGF-2 (Fibroblast growth factor-2), HGF (Hepatocyte growth factor), MMP-2/9 and so on [88]. Additionally, data have shown that mTORC2 is involved in MM angiogenesis [89], and activation of the PI3K/AKT/mTOR pathway regulates pro-angiogenic factors of MMP-9 [90]. Besides, Notch signaling has been investigated in the cross talk between endothelial cells and MM cells to enable angiogenesis [91]. Consistently, in our GO enrichment analysis, the up-regulated DEGs were enriched in positive regulation of angiogenesis, suggesting the role of angiogenesis of myeloma SP cells. In a conclusion, neovascularization and positive regulation of angiogenesis may be regarded as potential factors in modulating MM progression and deserving prognostic roles. Nevertheless, its biological mechanisms have not clearly revealed yet. In the future, additional studies are needed to further confirm the mechanisms of angiogenesis in myeloma SP cells.
Extramedullary disease of MM remains a key area of therapeutic challenge, and the expression of adhesion molecules and changes in angiogenesis concerning mostly VEGF, MMP-9 and others are involved in extramedullary spread of MM cells [92,93]. Plasma cells from extramedullary plasmacytomas showed angiogenesis related expression [94]. What's more, neovascularization can promote the growth of plasmacytomas [95]. These discoveries supported the idea that increased angiogenesis could facilitate malignant plasma cells growth outside the BM microenvironment. However, there are still unsolved questions on extramedullary myeloma involvement, especially on the relevant association with myeloma SP cells, which require further study.
In our survival analysis, to evaluate the association between the DEGs and clinical survival of MM patients and predict the prognosis of MM patients, we revealed 9 DEGs including TFF3, EPDR1, MACROD1, ARHGEF12, AMMECR1, NFATC2, HES6, PLEK2 and SNCA to be significantly associated with survival and established a survival prediction model based on the 9 genes. Stratified by risk score, a significantly different clinical outcome of MM patients were showed by the Kaplan-Meier curve in both training and validation datasets. However, further investigation of these genes in clinical research is warranted.
There are some limitations in our study. Firstly, the identification of DEGs profile was performed without external validation of other databases because of the absence of available data about SP cells compared to MP cells in MM. Second, we didn't evaluate the correlation of the prognostic model with clinicopathological characteristics. Thirdly, our study was only analyzed based on bioinformatics analysis. Hence, further investigations are warranted to validate the results and enhance our understanding of the biological role of these genes in MM.

Conclusions
To sum up, we performed a comprehensive bioinformatics analysis on microarray data of myeloma LC/SP cells. DEGs were identified to be significantly enriched in various pathways, especially positive regulation of angiogenesis and cell adhesion. The results of this study increase our understanding of novel biomarkers of myeloma LC/SP cells, prediction of MM prognosis and potential therapeutic targets. Nevertheless, further relevant studies are needed to confirm the identified DEGs and pathways in LC/SP cells of MM.