Skip to main content

Differential long noncoding RNA/mRNA expression profiling and functional network analysis during osteogenic differentiation of human bone marrow mesenchymal stem cells

Abstract

Background

Mesenchymal stem cells (MSCs) are the most promising cell types for bone regeneration and repair due to their osteogenic potential. MSC differentiation is precisely regulated and orchestrated by the mechanical and molecular signals from the extracellular environment, involving complex pathways regulated at both the transcriptional and post-transcriptional levels. However, the potential role of long noncoding RNA (lncRNA) in the osteogenic differentiation of human MSCs remains largely unclear.

Methods

Here, we undertook the survey of differential coding and noncoding transcript expression profiling and functional network analysis during osteogenic differentiation of human bone marrow mesenchymal stem cells (BMSCs) using human whole transcriptome microarray. The key pathways, mRNAs, and lncRNAs controlling osteogenic differentiation of BMSCs were identified by further bioinformatic analysis. The role of lncRNA in the osteogenic differentiation of MSCs was verified by lncRNA overexpression or knockdown methods.

Results

A total of 1269 coding transcripts with 648 genes significantly upregulated and 621 genes downregulated, and 1408 lncRNAs with 785 lncRNAs significantly upregulated and 623 lncRNAs downregulated were detected along with osteogenic differentiation. Bioinformatic analysis identified that several pathways may be associated with osteogenic differentiation potentials of BMSCs, such as the MAPK signaling pathway, the Jak-STAT signaling pathway, the Toll-like receptor signaling pathway, and the TGF-beta signaling pathway, etc. Bioinformatic analysis also revealed 13 core regulatory genes including seven mRNAs (GPX3, TLR2, BDKRB1, FBXO5, BRCA1, MAP3K8, and SCARB1), and six lncRNAs (XR_111050, NR_024031, FR374455, FR401275, FR406817, and FR148647). Based on the analysis, we identified one lncRNA, XR_111050, that could enhance the osteogenic differentiation potentials of MSCs.

Conclusions

The potential regulatory mechanisms were identified using bioinformatic analyses. We further predicted the interactions of differentially expressed coding and noncoding genes, and identified core regulatory factors by co-expression networks during osteogenic differentiation of BMSCs. Our results could lead to a better understanding of the molecular mechanisms of genes and lncRNAs, and their cooperation underlying MSC osteogenic differentiation and bone formation. We identified that one lncRNA, XR_111050, could be a potential target for bone tissue engineering.

Background

Mesenchymal stem cells (MSCs) are the most promising cell types for bone regeneration and repair due to their osteogenic potential [1, 2]. Bone marrow mesenchymal stem cells (BMSCs) are considered the gold standard for use in bone tissue regeneration among MSCs, and they are an important source of multipotent progenitor cells with self-renewal capacity that can differentiate into osteoblasts, chondrocytes, and adipocytes with great potential for clinic applications [1, 3]. BMSC differentiation is precisely regulated and orchestrated by the mechanical and molecular signals from the extracellular environment, involving complex pathways regulated at both the transcriptional and post-transcriptional levels [46].

Long noncoding RNAs (lncRNAs) are a class of nonprotein coding genes. Recent advances in high-throughput technology and computational methods allow for an unprecedented analysis of such transcripts. Tens of thousands of lncRNAs transcribed from mammalian genomes have been identified [712]. Recent studies have shown the regulatory role and functional diversity of lncRNAs. Many functional lncRNAs show tissue-specific expression patterns and distinct subcellular localizations [9, 11, 13, 14]. Growing evidence suggests that lncRNAs are important modulators of gene expression via their interaction with DNA, RNA, or protein [1525], while determining the function of individual lncRNAs remains a challenge. Some lncRNAs have been shown to regulate osteogenic differentiation of MSCs [2631]. However, the particular functions of lncRNAs during osteogenic differentiation of human BMSCs still remain largely unclear.

A previous study found that osteogenic differentiation markers were enhanced and that the stemness marker NANOG was downregulated at day 7 of osteogenic differentiation induction in BMSCs, indicated that day 7 of the induction was the transit point of the onset of osteogenic differentiation in BMSCs [32]. In the present study, we aimed to investigate the key mRNAs and lncRNAs for the transit point of osteogenic differentiation in MSCs. Here, we performed microarray analysis to identify the differential coding and long noncoding transcript expression profiling between noninduced and day 7 osteogenic-differentiated human BMSCs. A total of 1269 coding transcripts and 1408 lncRNAs were identified with differential expression during osteogenic differentiation in BMSCs. We further predicted interactions between coding and noncoding genes, and identified core regulatory factors by bioinformatic analysis. Our results also highlight the significant involvement of one lncRNA, XR_111050, in the positive regulation of osteogenic differentiation of BMSCs.

Results

The profiling of differentially expressed mRNAs and lncRNAs during osteogenic differentiation of BMSCs

To identify the differentially expressed lncRNAs and mRNAs during osteogenic differentiation of MSCs, we first applied the RVM t test to filter the genes that were differentially expressed, and then the differentially expressed genes with 1.5-fold changes were selected according to the p value threshold false discovery rate (FDR) for subsequent analysis. A total of 1269 coding transcripts with differential expression were identified during osteogenic differentiation (p < 0.05, FDR < 0.05); of these, 648 were upregulated and 621 were downregulated in BMSCs at 7 days after osteogenic induction compared to the uninduced BMSCs (Additional file 1: Table S1). In addition, a total of 1408 lncRNAs were differentially expressed after osteogenic differentiation (p < 0.05, FDR < 0.05); of these, 785 were upregulated and 623 were downregulated in BMSCs at 7 days after osteogenic induction compared to the uninduced BMSCs (Additional file 2: Table S2).

To confirm the reliability of the microarray data, we randomly selected six differentially expressed mRNAs (FGF7, FOXO1, STAT2, TMEM156, FOXM1, and FOXC2), and nine differentially expressed lncRNAs (XR_111050, NR_003255, NR_027621, NR_002196, NR_045555, NR_024593, NR_034115, NR_034181, and NR_037595) to analyze their expression with real-time polymerase chain reaction (RT-PCR). The RT-PCR results confirmed that the expression levels of these randomly selected 15 genes were consistent with the microarray results (Table 1).

Table 1 The real-time polymerase chain reaction (RT-PCR) results showed differential mRNA and lncRNA expression levels during osteogenic differentiation of bone marrow mesenchymal stem cells

Bioinformatic analysis of microarray data during osteogenic differentiation of BMSCs

To identify the key factors that controlled osteogenic differentiation of BMSCs, we performed gene ontology (GO) analysis and pathway analysis. First, GO analysis was performed to enrich the significant functions from the differentially expressed genes obtained during osteogenic differentiation of BMSCs. We obtained 171 GO functions according to p value and FDR (p < 0.05, FDR < 0.05) including 99 upregulated GO functions and 72 downregulated GO functions (Additional file 3: Table S3). The negative logarithm of the p value (–LgP) was used to represent the correlation between gene expression and the relevant biological process. Some important upregulated GO functions may be related to osteogenic differentiation, including response to stimulus, DNA-dependent transcription, ion transport, cell adhesion, and skeletal system development, and some important downregulated GO functions that were related to osteogenic differentiation were cell cycle, cell division, mitosis, DNA replication, and DNA-dependent transcription (Fig. 1a; Additional file 3: Table S3). We obtained 331 significantly upregulated genes and 297 significantly downregulated genes from enriched GO functions (p <0.01, FDR < 0.05).

Fig. 1
figure 1

Significant gene ontology (GO) analyses and pathways of differentially expressed genes during osteogenic differentiation. a The significant GO of differentially expressed genes during osteogenic differentiation. b The significant pathways of differentially expressed genes during osteogenic differentiation. The y axis shows the GO or pathway category and the x axis shows the negative logarithm of the p value (–LgP). A larger –LgP indicated a smaller p value for the difference. ECM extracellular matrix, TGF transforming growth factor

We then enriched the significantly changed pathways that mediated the functions of the differentially expressed genes based on the KEGG database, and identified a total of 67 significant pathways related to differential gene expression which may play key roles during osteogenic differentiation of BMSCs (p < 0.05); of these, 37 pathways involved upregulated genes and 30 pathways involved downregulated genes (Fig. 1b; Additional file 4: Table S4). Some important pathways that may related to osteogenic differentiation were enriched from upregulated genes, including cytokine-cytokine receptor interaction, extracellular matrix (ECM)-receptor interaction, the Toll-like receptor signaling pathway, and the chemokine signaling pathway. Some important pathways may be related to osteogenic differentiation enriched from downregulated genes were cell cycle, DNA replication, pyrimidine metabolism, and purine metabolism (Fig. 1b; Additional file 4: Table S4).

Among the significant pathways related to differential genes, a total of 170 upregulated genes and 152 downregulated genes was identified according to the p value (p < 0.05; Additional file 5: Table S5). The Path-Net was the interaction net of the significant pathways of the differential expression genes, and was built according to the interaction among pathways of the KEGG database to find the interaction between the significant pathways directly and systemically. After performing Path-net analysis to generate an interaction network of these significantly changed pathways, we could identify the relationship between the pathways enriched from differential expression genes and identify why a certain pathway was activated (Fig. 2). Our results identified the crosstalk of key pathways which may have an important role in the differentiation regulation of BMSCs. The crosstalk among the MAPK signaling pathway, the Jak-STAT signaling pathway, the Toll-like receptor signaling pathway, and the transforming growth factor (TGF)-beta signaling pathway were identified to a high degree (Fig. 2; Additional file 6: Table S6). We observed a high degree of increased activity of apoptosis, the cytokine-cytokine receptor interaction, the TGF-beta signaling pathway, the Jak-STAT signaling pathway, and the Toll-like receptor signaling pathway in a cluster of upregulating genes, and a decreased activity of the MAPK signaling pathway, the cell cycle, and the p53 signaling pathway in a cluster of downregulating genes (Fig. 2; Additional file 6: Table S6).

Fig. 2
figure 2

The interaction network of significant pathways (Path-net). The role of each pathway in the network was measured by counting its connections to upstream and downstream pathways. A pathway with a high degree implied that it played an important role in the signaling network. Blue represents downregulated pathways, red represents upregulated pathways, and yellow represents up- and downregulated pathways. The lines indicate interactions between pathways. ECM extracellular matrix, TGF transforming growth factor

We constructed an interactions repository by a Signal-net analysis (Signal-net) based on the significantly regulated GOs and pathways to identify 92 core genes during osteogenic differentiation of BMSCs according to the degree of gene interaction (betweenness centrality) (Fig. 3; Additional file 7: Table S7).

Fig. 3
figure 3

The interaction network of differentially expressed genes (Signal-net). In the Signal-net, the genes are characterized by measuring their “betweenness centrality”, the number of times a node is located in the shortest path between two other nodes. This measure reflects the importance of a node in a network relative to another. The circles represent important functional genes (red: upregulated genes; blue: downregulated genes); the circle size represents the degree of interaction (betweenness centrality), and lines indicate the interactions

To obtain core mRNAs and lncRNA genes during osteogenic differentiation of MSCs, we calculated the Pearson correlation for each pair of genes and used the most significantly correlated pairs to construct a coding-noncoding gene co-expression (CNC) network to identify the interactions and significance degree amongst the differentially expressed lncRNAs and mRNAs according to clustering coefficient and degree (Additional file 8: Figure S1 and Additional file 9: Figure S2). The closely correlated mRNAs and lncRNAs were identified. One mRNA was correlated with one or more lncRNAs and vice versa. We identified 13 core regulatory genes including seven mRNAs (GPX3, TLR2, BDKRB1, FBXO5, BRCA1, MAP3K8, and SCARB1) and six lncRNAs (XR_111050, NR_024031, FR374455, FR401275, FR406817, and FR148647) according to clustering coefficient and degree value (Additional file 10: Table S8; clustering coefficient ≥0.9 and degree >7).

The lncRNA XR_111050 enhanced the osteogenic differentiation potentials of MSCs

In the CNC networks, among the six lncRNAs, XR_111050 showed the most significant clustering coefficient. The RT-PCR results showed that the expression of XR_111050 increased along with osteogenic differentiation of BMSCs (Fig. 4a). To evaluate the effect of XR_111050 on osteogenesis, we inserted the XR_111050 sequence into a lentiviral vector, which was then transduced into BMSCs via lentiviral infection. Enhanced expression of XR_111050 was verified with RT-PCR in BMSCs (Fig. 4b). Two weeks after osteogenic induction, our results showed that XR_111050 overexpression enhanced osteogenic differentiation of BMSCs as shown by Alizarin red staining and quantitative calcium measurements compared with cells infected with the empty vector (Fig. 4c and d). RT-PCR analysis of XR_111050-overexpressing BMSCs showed stronger expression of osteogenic markers (including COL1A2 at 0, 3, and 10 days, OCN and OPN at 0, 3, and 7 days, and BSP at 7 and 14 days after osteogenic induction) compared with the control group (Fig. 4e–h). We also examined expression of the key transcription factors RUNX2 and OSX for regulating osteogenic differentiation. RT-PCR results showed that XR_111050-overexpressing BMSCs had a significantly higher RUNX2 mRNA level compared with cells infected with the empty vector (Fig. 4i). On the other hand, OSX expression levels did not significantly differ between groups (data not shown). We designed and introduced a short-hairpin RNA (shRNA) to target XR_111050. After selection, the knockdown efficiency (80%) was verified by RT-PCR (Fig. 5a). Two weeks after osteogenic induction, Alizarin red staining and quantitative calcium measurements revealed that mineralization was markedly decreased after depletion of XR_111050 in BMSCs (Fig. 5b and c).

Fig. 4
figure 4

XR_111050 enhanced the osteogenic differentiation of BMSCs. The BMSCs were transduced with XR_111050 lentiviral vector or with empty vector. a XR_111050 expression increased during osteogenic differentiation as evaluated by qRT-PCR. b The enhanced XR_111050 expression in BMSCs was verified with qRT-PCR. The increased XR_111050 expression enhanced mineralization of BMSCs as shown by c Alizarin red staining and d calcium quantitative analysis. qRT-PCR results showed that enhanced XR_111050 expression resulted in upregulated expressions of COL1A2 (e), OCN (f), BSP (g) and OPN (h) after osteogenic induction. i The enhanced XR_111050 expression increased RUNX2 expression in BMSCs. GAPDH was used as an internal control. Student’s t test was performed to determine statistical significance. Error bars represent SD (n = 3). *p < 0.05, **p < 0.01. BSP bone sialoprotein, COL1A2 collagen a-2(I) chain, d days, GADPH glyceraldehyde-3-phosphate dehydrogenase, OPN osteopontin, RUNX2 runt-related transcription factor 2, W weeks

Fig. 5
figure 5

XR_111050 knockdown inhibited osteogenic differentiation of BMSCs. The BMSCs was transduced with short-hairpin RNAs (shRNA) to silence XR_111050 (XR_111050sh) or with control shRNA (Contsh). a The reduced XR_111050 expression by shRNA was evaluated by qRT- PCR. XR_111050 knockdown inhibited mineralization, as shown by results of Alizarin red staining (b) and calcium quantitative analysis (c). Error bars represent SD (n = 3). **p < 0.01. GADPH glyceraldehyde-3-phosphate dehydrogenase

We then detected the expression of BDKRB2, MAP3K8, SCARB1, TLR2, FR000997, FR374455, NR_002744, NR_027293, and NR_034115 in BMSCs after overexpression of XR_111050 to investigate possible correlations. The RT-PCR results showed that overexpression of XR_111050 upregulated the expression of BDKRB2, MAP3K8, TLR2, and FR000997 (Fig. 6ad), and downregulated the expression of NR_002744 and NR_034115 compared with control group (Fig. 6e and f), but the expression of SCARB1, FR374455, and NR_027293 were not changed (data not shown).

Fig. 6
figure 6

XR_111050 regulated the expression of BDKRB2, MAP3K8, TLR2, FR000997, NR_002744, and NR_034115 in BMSCs. qRT-PCR results showed that enhanced XR_111050 expression resulted in upregulated expression of BDKRB2 (a), MAP3K8 (b), TLR2 (c), and FR000997 (d), and downregulated expression of NR_034115 (e) and NR_002744 (f). GAPDH was used as an internal control. Student’s t test was performed to determine statistical significance. Error bars represent SD (n = 3). *p < 0.05, **p < 0.01. GADPH glyceraldehyde-3-phosphate dehydrogenase

To confirm whether XR_111050 had similar functions in other MSCs, we detected the XR_111050 expression in periodontal ligament stem cells (PDLSCs) and the RT-PCR results showed that the expression of XR_111050 increased along with osteogenic differentiation of PDLSCs (Additional file 11: Figure S3A). Then we overexpressed XR_111050 in PDLSCs via lentiviral infection. Enhanced expression of XR_111050 was verified with RT-PCR in PDLSCs (Additional file 11: Figure S3B). Two weeks after osteogenic induction, Alizarin red staining and quantitative calcium measurements revealed that mineralization was markedly increased after overexpression of XR_111050 in PDLSCs (Additional file 11: Figure S3C,D).

Discussion

MSC differentiation is precisely regulated and orchestrated by mechanical and molecular signals that have increasingly attracted great attention in recent years [4, 6, 3338]. Therefore, understanding the molecular mechanisms underlying MSC differentiation and modulating lineage commitment of MSCs is of importance for therapeutic purposes [39].

Intensive studies have demonstrated that a number of critical signaling pathways are involved in regulating the lineage commitment of MSCs, including TGF-β/BMP signaling, Wnt signaling, Hedgehogs, Notch, and FGFs [4, 4043]. A large range of biological processes have been reported to be associated with lncRNAs, which are correlated with various human diseases and life processes. Emerging evidence is revealing that lncRNAs regulate differentiation of stem cells [44]. Some lncRNAs have been reported to be key regulators of osteogenic differentiation of MSCs in recent studies [2731, 45, 46]. However, the role of lncRNAs on the lineage commitment of MSCs remains largely unknown.

Here, we organized these studies to outline the profiling of mRNAs and lncRNAs during osteogenic differentiation of BMSCs. By comparing the lncRNA and mRNA expressions of differentiated and undifferentiated BMSCs, we identified 1269 differentially expressed mRNAs and 1408 differentially expressed lncRNAs. Our results provided a dynamic variation profiling in the expressions of lncRNAs and mRNAs during osteogenic differentiation of BMSCs, and implied that upregulated lncRNAs and mRNAs were predominant during osteogenic differentiation of BMSCs. The differentially expressed mRNAs and lncRNAs obtained were used as candidates to screen the key genes controlling osteogenic differentiation of BMSCs by further bioinformatic analysis.

We screened for the important pathways involved in the differentiation mechanism of BMSCs. Among the important pathways, the cytokine-cytokine receptor interaction pathway was demonstrated to be the dominant pathway and 28 genes were enriched. Of these, TNFSF10 was identified as being significantly upregulated, which functions as a negative regulator of bone resorption by activating TNFRSF11B. The interaction or correlation of these pathways were built with Path-Net (the interaction net of the significant pathways of the differential expression genes) according to the interaction amongst pathways of the KEGG database to find the interaction among the significant pathways directly and systemically. A pathway with a high degree implied that it played an important role in the signaling network. The Path-net analysis results indicated that the MAPK signaling pathway, the Jak-STAT signaling pathway, the Toll-like receptor signaling pathway, and the TGF-beta signaling pathway may have important roles in the differentiation regulation of BMSCs to a high degree. The results are consistent with previous reports [4, 41, 4750]. The Path-net identified in our results also provided insights into the role of crosstalk among key signaling pathways regulating osteogenic differentiation of BMSCs. Integration of signaling into signaling networks through cross-talk with other signal transduction pathways contributes to induction of specific activation and differentiation pathways [51]. We identified that the Toll-like receptor signaling pathway was involved in osteogenic differentiation of BMSCs by targeting apoptosis, cytokine-cytokine receptor interaction, the Jak-STAT signaling pathway, and the MAPK signaling pathway. Recent studies implied that signal transduction crosstalk is regulated in a dynamic manner, and that dysregulation of signal transduction crosstalk may contribute to pathogenesis [52]. For example, our results showed that significantly upregulated IL6 was associated with multiple signaling pathways, including the cytokine-cytokine receptor interaction, the Toll-like receptor signaling pathway, and the Jak-STAT signaling pathway, which can enhance osteogenic differentiation [4, 53].

We then predicted some core regulating factors using bioinformatic analysis of differentially expressed mRNAs and lncRNAs. We constructed a potential interaction between mRNAs and lncRNAs by CNC network, and identified key genes, such as GPX3, TLR2, BDKRB1, MAP3K8, and SCARB1, and lncRNAs, such as XR_111050 and NR_024031. TLR2 and MAP3K8, as component of the Toll-like receptor and MAPK signaling pathways which is upregulated during osteogenic differentiation of BMSCs in this study, may participate in osteogenic differentiation via positive or negative regulation of lncRNAs. Studies have suggested that the osteogenic differentiation seems to be enhanced by TLR2 activation in human BMSCs [50, 54, 55], which may work through NF-κB activation [55, 56], while MAP3K8 is an important mediator of intracellular mechanotransduction in human MSCs which is related to shear stress that can induce osteogenic differentiation of human MSCs [57, 58]. These data may provide clues to the potential role of interaction between lncRNAs and TLR signaling or MAPK signaling on osteogenic differentiation of BMSCs. GPX3, BDKRB1, and SCARB1 are involved in bone formation and absorption [5961]. GPX3 belongs to the glutathione peroxidase family, which catalyzes the reduction of organic hydroperoxides and hydrogen peroxide by glutathione, and thereby protect cells against oxidative damage [62]. A previous investigation discovered that elderly osteoporotic women with hip fractures had increased expression of GPX3, suggesting increased GPX3 antioxidative activity in bone samples [59]. Bdkrb1 plays an important role in inflammation and healing; Bdkrb1 knockout mice displayed increased bone loss, and bone marrow cells from Bdkrb1 knockout mice exhibited enhanced differentiation into osteoclasts [60]. The Scarb1 gene product is a high-density lipoprotein receptor which was shown to influence bone metabolism; Scarb1 deficiency promoted osteoblastogenesis but stunted terminal osteocyte differentiation [61]. When studying NR_024031, also named DANCR, researchers found that DANCR promoted the expression of IL-6 and TNF-α, and DANCR-induced IL-6 and TNF-α had bone-resorbing activity indicating that DANCR is involved in the pathology of osteoporosis and may be a biomarker for postmenopausal osteoporosis [63]. DANCR can inhibit the Wnt/β-catenin signal pathway and suppress the odontogenic differentiation in human dental pulp cells [64]. Moreover, DANCR can be activated by Sox4 and promotes cell proliferation and chondrogenesis in synovium-derived MSCs [65]. These reports suggest that GPX3, BDKRB1, SCARB1, and DANCR may play important roles in the osteogenic differentiation of BMSCs.

In the present study, we selected one candidate lncRNA, XR_111050, to investigate its function on the osteogenic differentiation of BMSCs. Our experiments showed that XR_111050 resulted in increased mineralization in BMSCs. Upregulated osteogenic marker (COL1A2, OCN, BSP, and OPN) mRNA levels verified that XR_111050 enhanced in vitro osteogenic differentiation in BMSCs. The activation of multiple transcription factors, such as RUNX2 and OSX, is associated with MSC osteogenic differentiation [66]. Our study found that XR_111050 stimulated the expression of one key transcription factor, RUNX2, but not the expression of another, OSX, which is involved in osteogenic differentiation. In addition, our results also confirmed that XR_111050 could enhance mineralization in PDLSCs. Taken together, our findings identified that XR_111050 was an enhancer of the osteogenic differentiation regulation of MSCs. XR_111050 is predicted to be a human hypothetical LOC100509635 miscRNA, a general term for a series of miscellaneous small RNA, which serves a variety of functions, including some enzyme-like catalysis and RNA processing [67]. Previous studies demonstrated that lncRNAs may function by controlling the transcriptional regulation of neighboring coding genes in a cis or trans manner [16]. Identifying differentially expressed nearby coding mRNA may enhance our understanding of the function of lncRNAs on the BMSCs. In this study, we found that XR_111050 is closely correlated with BDKRB2, MAP3K8, SCARB1, TLR2, FR000997, FR374455, NR_002744, NR_027293, and NR_034115 by lncRNA–mRNA CNC network analysis. We have detected their expression in BMSCs after overexpression of XR_111050 to investigate their correlations and possible mechanism. The results showed that overexpression of XR_111050 enhanced the expression of BDKRB2, MAP3K8, TLR2, and FR000997, and inhibited the expression of NR_002744 and NR_034115. Our microarray data showed that the expression of BDKRB2, MAP3K8, TLR2, and FR000997 was increased, and that of NR_002744 and NR_034115 was decreased after osteogenic differentiation in BMSCs, consistent with the results of overexpression of XR_111050, indicating that BDKRB2, MAP3K8, TLR2, and FR000997 might be positive regulators and NR_002744 and NR_034115 might be negative regulators for osteogenic differentiation of BMSCs. In conclusion, these results implied that XR_111050 may play a role in the osteogenic differentiation regulation of MSCs via a large number of cooperators or downstream targets, such as DKRB2, MAP3K8, TLR2, FR000997, NR_002744, and NR_034115. However, further studies must be performed to investigate this hypothesis.

Conclusions

We identified the differential expression profiling of coding and long noncoding transcripts during osteogenic differentiation of BMSCs with human whole transcriptome microarray. The potential regulatory mechanisms were identified with bioinformatic analyses. We further predicted the interactions of differentially expressed coding and noncoding genes, and identified core regulatory factors by co-expression networks during osteogenic differentiation of BMSCs. Our results also suggested that one lncRNA, XR_111050, is required for osteogenic differentiation of MSCs. Our studies provide a base for further understanding the precise role and mechanisms of these lncRNAs and mRNAs controlling osteogenic differentiation of MSCs and provide an important resource for investigation of MSC biology. In addition, our obtained lncRNAs with differential expression during MSC osteogenic differentiation could be potential targets for bone tissue engineering.

Methods

Cell culture and osteogenic differentiation

BMSCs derived from 18- to 20-year-old males (n = 3) were obtained from Cyagen Biosciences (Guangzhou, China). Human impacted third molars were collected from healthy male patients (16–20 years old). Teeth were first disinfected with 75% ethanol and then washed with phosphate-buffered saline. PDLSCs were isolated, cultured, and identified as previously described [68]. Briefly, PDLSCs were separated from periodontal ligament in the middle one-third of the root. Subsequently, the periodontal ligament tissue was digested in a solution of 3 mg/mL collagenase type I (Worthington Biochemical Corp., Lakewood, NJ, USA) and 4 mg/mL dispase (Roche Diagnostics Corp., Indianapolis, IN, USA) for 1 h at 37 °C. Single-cell suspensions were obtained by cell passage through a 70-μm strainer (Falcon, BD Labware, Franklin Lakes, NJ, USA). MSC cultures were grown in a humidified, 5% CO2 incubator at 37 °C in DMEM alpha modified Eagle’s medium (Invitrogen, Carlsbad, CA, USA), supplemented with 15% fetal bovine serum (FBS; Invitrogen, Carlsbad, CA, USA), 2 mmol/L glutamine, 100 U/mL penicillin, and 100 μg/mL streptomycin (Invitrogen, Carlsbad, CA, USA). The culture medium was changed every 3 days. All MSCs were used in subsequent experiments after 3–5 passages. For osteogenesis differentiation, MSCs were seeded at a density of 2.0 × 105 cells/well into six-well plates with routine medium. When cells reached 80–90% confluence, the medium was changed, and cells were grown in osteogenic differentiation medium with the STEMPRO osteogenesis differentiation Kit (Invitrogen). The osteogenic differentiation medium was replaced every 3 days.

Total RNA isolation and microarray hybridization

BMSCs were grown in osteogenic differentiation medium for 7 days using the StemPro osteogenesis differentiation kit (Invitrogen); the osteogenic differentiation medium was changed every 3 days. Three total RNA samples from BMSCs of three different individuals were extracted separately using TRIzol and the RNeasy mini kit (Qiagen, Germany). RNA quality and quantity were confirmed by multiImager and spectrophotometer (Meriton, China). Each total RNA sample was further purified using the RNeasy mini kit and RNase-Free DNase Set (QIAGEN, GmBH, Germany).

Microarray analysis was performed using Human Transcriptome Array (Affymetrix), covering more than 285,000 coding and noncoding transcripts (>245,000 coding transcripts, >40,000 noncoding transcripts). For each experimental group, three biological replicates (three BMSCs RNA samples from three different individuals) were hybridized. The RNA labeling and microarray hybridization were carried out according to the Affymetrix expression analysis technical manual (by Genminix Informatics Ltd., Shanghai, China). The arrays were scanned using the GeneChip® scanner 3000 (Affymetrix, CA, USA). The microarray raw data were normalized with Affymetrix Expression Console software using the MAS5 statistical algorithm.

Quantitative RT-PCR

qRT-PCR was used to verify the differential expression of genes that were detected on the microarray and to detect the expression of osteogenic differentiation markers. At days 0, 3, 7, 10, and 14 after osteogenic induction, BMSCs were harvested for RNA extraction. Total RNA was isolated from BMSCs with Trizol reagents (Invitrogen). For qRT-PCR, 2 μg aliquots of RNA as a template were combined with random hexamers and reverse transcriptase, according to the manufacturer’s protocol (Invitrogen). qRT-PCR was performed using QuantiTect SYBR Green PCR kit (Qiagen, Germany) and an iCycler iQ Multicolor Real-Time PCR Detection System. Primer sequences used for qRT-PCR are listed in Additional file 12 (Table S9). Relative mRNA levels were calculated using the 2−ΔΔCt method. The ΔCts were obtained from Ct normalized with GAPDH. Pearson’s correlation coefficient was further calculated for each gene using the normalized data to quantify the consistency between microarray experiments and qRT-PCR (p < 0.05 and R > 0.9).

Bioinformatic analysis

The bioinformatic analysis was performed (by Genminix Informatics Ltd., Shanghai, China) as described previously [68]. Briefly, after normalized data were compared and filtered, the differentially expressed genes for the control and experiment group were firstly filtered using RVM t test (Two Class Dif) (p < 0.05) for following analysis. Hierarchical clusters were performed by EPCLUST.

GO analysis was applied to analyze the main function of the differential expression genes according to the GO using the two-sided Fisher’s exact test and χ 2 test for classifying the GO category (p < 0.05, FDR < 0.01). Enrichment of the GO category was calculated for significance of the function. Go-map analysis was made to identify the interaction net of the significant GOs of the differential expression genes.

Pathway analysis was used to discover the significant pathway of the differential genes according to KEGG, Biocarta and Reatome, using Fisher’s exact test and χ 2 test (p < 0.05, FDR < 0.01). The enrichment Re of pathways was calculated. The interaction net of the significant pathways of the differential expression genes (Path-Net) was built to find the interaction among the significant pathways directly and systemically according to the interaction among pathways of the KEGG database.

The gene-gene interaction network (Signal-Net) from differentially expressed genes was constructed to identify the molecular networks between the two genes based on the interaction database from KEGG.

We built a coding-noncoding gene co-expression network (CNC) to further identify the interactions and locate core regulatory factors (genes) among genes according to the normalized signal intensity of specific expression genes.

Plasmid construction and viral infection

The plasmids were constructed with standard methods; all structures were verified by appropriate restriction digestion and/or sequencing. Human full-length XR-111050 cDNA from BMSCs was produced with a standard PCR protocol. This sequence was subcloned into the LV5 lentiviral vector (Genepharma Company, Suzhou, China). Short-hairpin RNAs (shRNA) with the complementary sequences of the target genes were subcloned into the LV3 lentiviral vector (Genepharma Company, Suzhou, China). For viral infections, MSCs were plated overnight and then infected with lentiviruses in the presence of polybrene (6 μg/mL; Sigma-Aldrich, St. Louis, MO, USA) for 6 h. After 48 h, infected cells were selected with 1 μg/mL puromycin for 7 days. The target sequences for the shRNA were: LV3 shRNA (Consh), 5’-TTCTCCGAACGTGTCACGTTTC-3’; XR-111050 shRNA (XR-111050sh), 5’-GGACGTGTCTTTCAGGGAAAG-3’.

Alizarin red detection and quantitative calcium analysis

MSCs were grown in osteogenic differentiation medium using the StemPro osteogenesis differentiation kit (Invitrogen). For detecting mineralization, cells were induced for 2 weeks, fixed with 70% ethanol, and stained with 2% Alizarin red (Sigma-Aldrich, St. Louis, MO, USA). To quantitatively determine calcium, Alizarin red was destained with 10% cetylpyridinium chloride in 10 mM sodium phosphate for 30 min at room temperature. The concentration was determined by measuring the absorbance at 562 nm on a multiplate reader and comparing to a standard calcium curve with calcium dilutions in the same solution. The final calcium level in each group was normalized to the total protein concentration detected in a duplicate plate.

Statistics

All statistical calculations were performed with SPSS11 statistical software. Statistical analyses included comparisons with the t test, Fisher’s exact test, \( \chi \) 2 test, and the Pearson correlation, as appropriate; p values less than 0.05 were considered statistically significant.

Abbreviations

BMSC:

Bone marrow mesenchymal stem cell

CNC:

Coding-noncoding gene co-expression

ECM:

Extracellular matrix

FDR:

False discovery rate

GO:

Gene ontology

lncRNA:

Long noncoding RNA

MSC:

Mesenchymal stem cell

Path-net:

Interaction network of significant pathways

PDLSC:

Periodontal ligament stem cell

(q)RT- PCR:

(Quantitative) real-time polymerase chain reaction

Signal-net:

Gene-gene interaction network

TGF:

Transforming growth factor

References

  1. Bianco P, Riminucci M, Gronthos S, Robey PG. Bone marrow stromal stem cells: nature, biology, and potential applications. Stem Cells. 2001;19(3):180–92.

    Article  CAS  PubMed  Google Scholar 

  2. Yousefi AM, James PF, Akbarzadeh R, Subramanian A, Flavin C, Oudadesse H. Prospect of stem cells in bone tissue engineering: a review. Stem Cells Int. 2016;2016:6180487.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Bianco P, Robey PG, Simmons PJ. Mesenchymal stem cells: revisiting history, concepts, and assays. Cell Stem Cell. 2008;2(4):313–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Chen Q, Shou P, Zheng C, Jiang M, Cao G, Yang Q, et al. Fate decision of mesenchymal stem cells: adipocytes or osteoblasts? Cell Death Differ. 2016;23(7):1128–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Deng P, Zhou C, Alvarez R, Hong C, Wang CY. Inhibition of IKK/NF-kappaB signaling enhances differentiation of mesenchymal stromal cells from human embryonic stem cells. Stem Cell Rep. 2016;6(4):456–65.

    Article  CAS  Google Scholar 

  6. Yi S, Yu M, Yang S, Miron RJ, Zhang Y. Tcf12, a member of basic helix-loop-helix transcription factors, mediates bone marrow mesenchymal stem cell osteogenic differentiation in vitro and in vivo. Stem Cells. 2017;35(2):386–97.

  7. Iyer MK, Niknafs YS, Malik R, Singhal U, Sahu A, Hosono Y, et al. The landscape of long noncoding RNAs in the human transcriptome. Nat Genet. 2015;47(3):199–208.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Gupta R, Ahn R, Lai K, Mullins E, Debbaneh M, Dimon M, et al. Landscape of long noncoding RNAs in psoriatic and healthy skin. J Invest Dermatol. 2016;136(3):603–9.

    Article  CAS  PubMed  Google Scholar 

  9. Necsulea A, Soumillon M, Warnefors M, Liechti A, Daish T, Zeller U, et al. The evolution of lncRNA repertoires and expression patterns in tetrapods. Nature. 2014;505(7485):635–40.

    Article  CAS  PubMed  Google Scholar 

  10. Guttman M, Amit I, Garber M, French C, Lin MF, Feldser D, et al. Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature. 2009;458(7235):223–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Yan X, Hu Z, Feng Y, Hu X, Yuan J, Zhao SD, et al. Comprehensive genomic characterization of long non-coding RNAs across human cancers. Cancer Cell. 2015;28(4):529–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Casero D, Sandoval S, Seet CS, Scholes J, Zhu Y, Ha VL, et al. Long non-coding RNA profiling of human lymphoid progenitor cells reveals transcriptional divergence of B cell and T cell lineages. Nat Immunol. 2015;16(12):1282–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Lennox KA, Behlke MA. Cellular localization of long non-coding RNAs affects silencing by RNAi more than by antisense oligonucleotides. Nucleic Acids Res. 2016;44(2):863–77.

    Article  CAS  PubMed  Google Scholar 

  14. Gloss BS, Dinger ME. The specificity of long noncoding RNA expression. Biochim Biophys Acta. 2016;1859(1):16–22.

    Article  CAS  PubMed  Google Scholar 

  15. Qi W, Song X, Li L. Long non-coding RNA-guided regulation in organisms. Sci China Life Sci. 2013;56(10):891–6.

    Article  CAS  PubMed  Google Scholar 

  16. Luo S, Lu JY, Liu L, Yin Y, Chen C, Han X, et al. Divergent lncRNAs regulate gene expression and lineage differentiation in pluripotent cells. Cell Stem Cell. 2016;18(5):637–52.

    Article  CAS  PubMed  Google Scholar 

  17. Lee S, Kopp F, Chang TC, Sataluri A, Chen B, Sivakumar S, et al. Noncoding RNA NORAD regulates genomic stability by sequestering PUMILIO proteins. Cell. 2016;164(1–2):69–80.

    Article  CAS  PubMed  Google Scholar 

  18. Zhou J, Yang L, Zhong T, Mueller M, Men Y, Zhang N, et al. H19 lncRNA alters DNA methylation genome wide by regulating S-adenosylhomocysteine hydrolase. Nat Commun. 2015;6:10221.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Huang W, Thomas B, Flynn RA, Gavzy SJ, Wu L, Kim SV, et al. DDX5 and its associated lncRNA Rmrp modulate TH17 cell effector functions. Nature. 2015;528(7583):517–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Pefanis E, Wang J, Rothschild G, Lim J, Chao J, Rabadan R, et al. Noncoding RNA transcription targets AID to divergently transcribed loci in B cells. Nature. 2014;514(7522):389–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Chen L, Kostadima M, Martens JH, Canu G, Garcia SP, Turro E, et al. Transcriptional diversity during lineage commitment of human blood progenitors. Science. 2014;345(6204):1251033.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Ard R, Tong P, Allshire RC. Long non-coding RNA-mediated transcriptional interference of a permease gene confers drug tolerance in fission yeast. Nat Commun. 2014;5:5576.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Yoon JH, Abdelmohsen K, Gorospe M. Posttranscriptional gene regulation by long noncoding RNA. J Mol Biol. 2013;425(19):3723–30.

    Article  CAS  PubMed  Google Scholar 

  24. Yang L, Lin C, Jin C, Yang JC, Tanasa B, Li W, et al. lncRNA-dependent mechanisms of androgen-receptor-regulated gene activation programs. Nature. 2013;500(7464):598–602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Sauvageau M, Goff LA, Lodato S, Bonev B, Groff AF, Gerhardinger C, et al. Multiple knockout mouse models reveal lincRNAs are required for life and brain development. eLife. 2013;2:e01749.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Wang L, Wang Y, Li Z, Li Z, Yu B. Differential expression of long noncoding ribonucleic acids during osteogenic differentiation of human bone marrow mesenchymal stem cells. Int Orthop. 2015;39(5):1013–9.

    Article  PubMed  Google Scholar 

  27. Jia Q, Jiang W, Ni L. Down-regulated non-coding RNA (lncRNA-ANCR) promotes osteogenic differentiation of periodontal ligament stem cells. Arch Oral Biol. 2015;60(2):234–41.

    Article  CAS  PubMed  Google Scholar 

  28. Zuo C, Wang Z, Lu H, Dai Z, Liu X, Cui L. Expression profiling of lncRNAs in C3H10T1/2 mesenchymal stem cells undergoing early osteoblast differentiation. Mol Med Rep. 2013;8(2):463–7.

    PubMed  Google Scholar 

  29. Huang Y, Zheng Y, Jia L, Li W. Long noncoding RNA H19 promotes osteoblast differentiation via TGF-β1/Smad3/HDAC signaling pathway by deriving miR-675. Stem Cells. 2015;33(12):3481–92.

    Article  CAS  PubMed  Google Scholar 

  30. Li H, Zhang Z, Chen Z, Zhang D. Osteogenic growth peptide promotes osteogenic differentiation of mesenchymal stem cells mediated by LncRNA AK141205-induced upregulation of CXCL13. Biochem Biophys Res Commun. 2015;466(1):82–8.

    Article  CAS  PubMed  Google Scholar 

  31. Zhuang W, Ge X, Yang S, Huang M, Zhuang W, Chen P, et al. Upregulation of lncRNA MEG3 promotes osteogenic differentiation of mesenchymal stem cells from multiple myeloma patients by targeting BMP4 transcription. Stem Cells. 2015;33(6):1985–97.

    Article  CAS  PubMed  Google Scholar 

  32. Tan J, Lu J, Huang W, Dong Z, Kong C, Li L, et al. Genome-wide analysis of histone H3 lysine9 modifications in human mesenchymal stem cell osteogenic differentiation. PLoS One. 2009;4(8):e6792.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Fan C, Jia L, Zheng Y, Jin C, Liu Y, Liu H, et al. MiR-34a promotes osteogenic differentiation of human adipose-derived stem cells via the RBP2/NOTCH1/CYCLIN D1 coregulatory network. Stem Cell Rep. 2016;7(2):236–48.

    Article  CAS  Google Scholar 

  34. Chen Q, Shou P, Zhang L, Xu C, Zheng C, Han Y, et al. An osteopontin-integrin interaction plays a critical role in directing adipogenesis and osteogenesis by mesenchymal stem cells. Stem Cells. 2014;32(2):327–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Zhang P, Liu Y, Jin C, Zhang M, Lv L, Zhang X, et al. Histone H3K9 acetyltransferase PCAF is essential for osteogenic differentiation through bone morphogenetic protein signaling and may be involved in osteoporosis. Stem Cells. 2016;34(9):2332–41.

    Article  CAS  PubMed  Google Scholar 

  36. Yue R, Zhou BO, Shimada IS, Zhao Z, Morrison SJ. Leptin receptor promotes adipogenesis and reduces osteogenesis by regulating mesenchymal stromal cells in adult bone marrow. Cell Stem Cell. 2016;18(6):782–96.

    Article  CAS  PubMed  Google Scholar 

  37. Luo J, Yang Z, Ma Y, Yue Z, Lin H, Qu G, et al. LGR4 is a receptor for RANKL and negatively regulates osteoclast differentiation and bone resorption. Nat Med. 2016;22(5):539–46.

    Article  CAS  PubMed  Google Scholar 

  38. Liu Z, Yao X, Yan G, Xu Y, Yan J, Zou W, et al. Mediator MED23 cooperates with RUNX2 to drive osteoblast differentiation and bone development. Nat Commun. 2016;7:11149.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Forbes SJ, Rosenthal N. Preparing the ground for tissue regeneration: from mechanism to therapy. Nat Med. 2014;20(8):857–69.

    Article  CAS  PubMed  Google Scholar 

  40. Shi Y, Chen J, Karner CM, Long F. Hedgehog signaling activates a positive feedback mechanism involving insulin-like growth factors to induce osteoblast differentiation. Proc Natl Acad Sci U S A. 2015;112(15):4678–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Zhen G, Wen C, Jia X, Li Y, Crane JL, Mears SC, et al. Inhibition of TGF-β signaling in mesenchymal stem cells of subchondral bone attenuates osteoarthritis. Nat Med. 2013;19(6):704–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Jiang Q, Du J, Yin X, Shan Z, Ma Y, Ma P, et al. Shh signaling, negatively regulated by BMP signaling, inhibits the osteo/dentinogenic differentiation potentials of mesenchymal stem cells from apical papilla. Mol Cell Biochem. 2013;383(1–2):85–93.

    Article  CAS  PubMed  Google Scholar 

  43. Liu N, Shi S, Deng M, Tang L, Zhang G, Liu N, et al. High levels of β-catenin signaling reduce osteogenic differentiation of stem cells in inflammatory microenvironments through inhibition of the noncanonical Wnt pathway. J Bone Miner Res. 2011;26(9):2082–95.

    Article  CAS  PubMed  Google Scholar 

  44. Hu S, Shan G. LncRNAs in stem cells. Stem Cells Int. 2016;2016:2681925.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Song WQ, Gu WQ, Qian YB, Ma X, Mao YJ, Liu WJ. Identification of long non-coding RNA involved in osteogenic differentiation from mesenchymal stem cells using RNA-Seq data. Genet Mol Res. 2015;14(4):18268–79.

    Article  CAS  PubMed  Google Scholar 

  46. Qu Q, Fang F, Wu B, Hu Y, Chen M, Deng Z, et al. Potential role of long non-coding RNA in osteogenic differentiation of human periodontal ligament stem cells. J Periodontol. 2016;87(7):e127–37.

    Article  PubMed  Google Scholar 

  47. Greenblatt MB, Shim JH, Glimcher LH. Mitogen-activated protein kinase pathways in osteoblasts. Annu Rev Cell Dev Biol. 2013;29:63–79.

    Article  CAS  PubMed  Google Scholar 

  48. James AW. Review of signaling pathways governing MSC osteogenic and adipogenic differentiation. Scientifica (Cairo). 2013;2013:684736.

    Google Scholar 

  49. Huang E, Zhu G, Jiang W, Yang K, Gao Y, Luo Q, et al. Growth hormone synergizes with BMP9 in osteogenic differentiation by activating the JAK/STAT/IGF1 pathway in murine multilineage cells. J Bone Miner Res. 2012;27(7):1566–75.

    Article  CAS  PubMed  Google Scholar 

  50. Qi C, Xiaofeng X, Xiaoguang W. Effects of toll-like receptors 3 and 4 in the osteogenesis of stem cells. Stem Cells Int. 2014;2014:917168.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Ivashkiv LB. Cross-regulation of signaling by ITAM-associated receptors. Nat Immunol. 2009;10(4):340–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Hu X, Chen J, Wang L, Ivashkiv LB. Crosstalk among Jak-STAT, Toll-like receptor, and ITAM-dependent pathways in macrophage activation. J Leukoc Biol. 2007;82(2):237–43.

    Article  CAS  PubMed  Google Scholar 

  53. Bastidas-Coral AP, Bakker AD, Zandieh-Doulabi B, Kleverlaan CJ, Bravenboer N, Forouzanfar T, et al. Cytokines TNF-alpha, IL-6, IL-17 F, and IL-4 differentially affect osteogenic differentiation of human adipose stem cells. Stem Cells Int. 2016;2016:1318256.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Kim HS, Shin TH, Yang SR, Seo MS, Kim DJ, Kang SK, et al. Implication of NOD1 and NOD2 for the differentiation of multipotent mesenchymal stem cells derived from human umbilical cord blood. PLoS One. 2010;5(10):e15369.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Pevsner-Fischer M, Morad V, Cohen-Sfady M, Rousso-Noori L, Zanin-Zhorov A, Cohen S, et al. Toll-like receptors and their ligands control mesenchymal stem cell functions. Blood. 2007;109(4):1422–32.

    Article  CAS  PubMed  Google Scholar 

  56. Cho HH, Shin KK, Kim YJ, Song JS, Kim JM, Bae YC, et al. NF-kappaB activation stimulates osteogenic differentiation of mesenchymal stem cells derived from human adipose tissue by increasing TAZ expression. J Cell Physiol. 2010;223(1):168–77.

    CAS  PubMed  Google Scholar 

  57. Yourek G, McCormick SM, Mao JJ, Reilly GC. Shear stress induces osteogenic differentiation of human mesenchymal stem cells. Regen Med. 2010;5(5):713–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Glossop JR, Cartmell SH. Effect of fluid flow-induced shear stress on human mesenchymal stem cells: differential gene expression of IL1B and MAP3K8 in MAPK signaling. Gene Expr Patterns. 2009;9(5):381–8.

    Article  CAS  PubMed  Google Scholar 

  59. Foger-Samwald U, Vekszler G, Horz-Schuch E, Salem S, Wipperich M, Ritschl P, et al. Molecular mechanisms of osteoporotic hip fractures in elderly women. Exp Gerontol. 2016;73:49–58.

    Article  PubMed  Google Scholar 

  60. Goncalves-Zillo TO, Pugliese LS, Sales VM, Mori MA, Squaiella-Baptistao CC, Longo-Maugeri IM, et al. Increased bone loss and amount of osteoclasts in kinin B1 receptor knockout mice. J Clin Periodontol. 2013;40(7):653–60.

    Article  CAS  PubMed  Google Scholar 

  61. Martineau C, Kevorkova O, Brissette L, Moreau R. Scavenger receptor class B, type I (Scarb1) deficiency promotes osteoblastogenesis but stunts terminal osteocyte differentiation. Physiol Rep. 2014;2(10):e12117.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Brigelius-Flohe R, Maiorino M. Glutathione peroxidases. Biochim Biophys Acta. 2013;1830(5):3289–303.

    Article  CAS  PubMed  Google Scholar 

  63. Tong X, Gu PC, Xu SZ, Lin XJ. Long non-coding RNA-DANCR in human circulating monocytes: a potential biomarker associated with postmenopausal osteoporosis. Biosci Biotechnol Biochem. 2015;79(5):732–7.

    Article  CAS  PubMed  Google Scholar 

  64. Chen L, Song Z, Huang S, Wang R, Qin W, Guo J, et al. lncRNA DANCR suppresses odontoblast-like differentiation of human dental pulp cells by inhibiting wnt/β-catenin pathway. Cell Tissue Res. 2016;364(2):309–18.

    Article  CAS  PubMed  Google Scholar 

  65. Zhang L, Chen S, Bao N, Yang C, Ti Y, Zhou L, et al. Sox4 enhances chondrogenic differentiation and proliferation of human synovium-derived stem cell via activation of long noncoding RNA DANCR. J Mol Histol. 2015;46(6):467–73.

    Article  CAS  PubMed  Google Scholar 

  66. Lian JB, Stein GS, Javed A, van Wijnen AJ, Stein JL, Montecino M, et al. Networks and hubs for the transcriptional control of osteoblastogenesis. Rev Endocr Metab Disord. 2006;7(1-2):1–16.

    Article  CAS  PubMed  Google Scholar 

  67. St Laurent G, Wahlestedt C, Kapranov P. The landscape of long noncoding RNA classification. Trends Genet. 2015;31(5):239–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Dong R, Du J, Wang L, Wang J, Ding G, Wang S, et al. Comparison of long noncoding RNA and mRNA expression profiles in mesenchymal stem cells derived from human periodontal ligament and bone marrow. BioMed Res Int. 2014;2014:317853.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by grants from the National Natural Science Foundation of China (81625005 to ZF, 81570936 to RD), the Dalian Science and Technology Project (2015E12SF160 to FW), the Beijing Key Laboratory Foundation of Tooth Regeneration and Function Reconstruction the Program for “Hundred-Thousand-Ten thousand” Talents in Beijing (2014006 to ZF), and the Chinese Medicine Science and Technology Foundation of Beijing (JJ2015-10 to ZF).

Availability of data and materials

The microarray datasets supporting the conclusions of this article are available. All other datasets supporting this article have been uploaded as part of the supplementary material.

Authors’ contributions

WZ was responsible for the collection and assembly of data, data analysis and interpretation, and final approval of the manuscript. RD was responsible for conception and design, the collection and assembly of data, data analysis and interpretation, financial support, and final approval of the manuscript. SD was responsible for the collection and assembly of data, drafting the manuscript, and final approval of the manuscript. JD was responsible for the provision of study material, collection of data, technical expertise, and revising the manuscript. ZF was responsible for conception and design, manuscript writing, financial support, and final approval of the manuscript. FW was responsible for the collection and assembly of data, data analysis and interpretation, manuscript writing, financial support, and final approval of the manuscript. All authors have read and approved the final version of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethical approval and consent to participate

All research involving human stem cells complied with the International Society for Stem Cell Research “Guidelines for the Conduct of Human Embryonic Stem Cell Research”. Human impacted third molars were collected under approved guidelines set by the Beijing Stomatological Hospital, Capital Medical University, with informed consent.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Zhipeng Fan or Fu Wang.

Additional files

Additional file 1: Table S1.

The differentially expressed mRNAs during osteogenic differentiation of BMSCs. (XLS 139 kb)

Additional file 2: Table S2.

The differentially expressed lncRNAs during osteogenic differentiation of BMSCs. (XLS 185 kb)

Additional file 3: Table S3.

The upregulated and downregulated GO functions during osteogenic differentiation of BMSCs. (XLS 54 kb)

Additional file 4: Table S4.

The upregulated and downregulated pathways during osteogenic differentiation of BMSCs. (XLS 33 kb)

Additional file 5: Table S5.

The upregulated and downregulated genes in significant pathways during osteogenic differentiation of BMSCs. (XLS 111 kb)

Additional file 6: Table S6.

The significant pathways enriched by Path-net during osteogenic differentiation of BMSCs. (XLS 25 kb)

Additional file 7: Table S7.

The 92 core genes enriched by Signal-net during osteogenic differentiation of BMSCs. (XLS 37 kb)

Additional file 8: Figure S1.

The co-expression network of mRNAs and lncRNAs in BMSCs before osteogenic differentiation (CNC). Circles represent upregulated (red) genes, and downregulated (blue) genes in BMSCs. The lncRNA genes are encircled in yellow. The lines represent the regulatory relationships between genes (solid lines represent positive correlations, dotted lines represent negative correlations). The circle size represents the degree of centrality. (PNG 3798 kb)

Additional file 9: Figure S2.

The co-expression network of mRNAs and lncRNAs in BMSCs after osteogenic differentiation (CNC). Circles represent upregulated (red) genes, and downregulated (blue) genes in BMSCs. The lncRNA genes are encircled in yellow. The lines represent the regulatory relationships between genes (solid lines represent positive correlations, dotted lines represent negative correlations). The circle size represents the degree of centrality. (PNG 4510 kb)

Additional file 10: Table S8.

The 13 core regulatory genes by CNC during osteogenic differentiation of BMSCs. (XLS 27 kb)

Additional file 11: Figure S3.

XR_111050 enhanced the osteogenic differentiation of PDLSCs. The PDLSCs were transduced with XR_111050 lentiviral vector or with empty vector. (A) XR_111050 expression increased during osteogenic differentiation evaluated by qRT-PCR. (B) The enhanced XR_111050 expression in PDLSCs was verified with qRT-PCR. (C,D) The enhanced XR_111050 expression enhanced mineralization of PDLSCs shown by Alizarin red staining (C) and calcium quantitative analysis (D). Error bars represent SD (n = 3). **p < 0.01. (TIFF 1312 kb)

Additional file 12: Table S9.

The primer sequences used for qRT-PCR. (XLS 30 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, W., Dong, R., Diao, S. et al. Differential long noncoding RNA/mRNA expression profiling and functional network analysis during osteogenic differentiation of human bone marrow mesenchymal stem cells. Stem Cell Res Ther 8, 30 (2017). https://doi.org/10.1186/s13287-017-0485-6

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13287-017-0485-6

Keywords