Comparison of Long Noncoding RNA and mRNA Expression Profiles in Mesenchymal Stem Cells Derived from Human Periodontal Ligament and Bone Marrow

Dong, Rui; Du, Juan; Wang, Liping; Wang, Jinsong; Ding, Gang; Wang, Songlin; Fan, Zhipeng

doi:https://doi.org/10.1155/2014/317853

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2014 | Article ID 317853 | https://doi.org/10.1155/2014/317853

Comparison of Long Noncoding RNA and mRNA Expression Profiles in Mesenchymal Stem Cells Derived from Human Periodontal Ligament and Bone Marrow

Rui Dong,¹Juan Du,¹Liping Wang,¹Jinsong Wang,^2,3Gang Ding,⁴Songlin Wang,^2,3and Zhipeng Fan¹

Academic Editor: Jinhua Yu

Received07 Jan 2014

Accepted16 Feb 2014

Published27 Mar 2014

Abstract

Mesenchymal stem cells (MSCs) in different anatomic locations possess diverse biological activities. Maintaining the pluripotent state and differentiation depend on the expression and regulation of thousands of genes, but it remains unclear which molecular mechanisms underlie MSC diversity. Thus, potential MSC applications are restricted. Long noncoding RNAs (lncRNAs) are implicated in the complex molecular circuitry of cellular processes. We investigated differences in lncRNA and mRNA expression profiles between bone marrow stem cells (BMSCs) and periodontal ligament stem cells (PDLSCs) with lncRNA microarray assays and bioinformatics analysis. In PDLSCs, numerous lncRNAs were significantly upregulated () or downregulated () compared to BMSCs. Furthermore, 1,578 mRNAs were differentially expressed. These genes implicated cellular pathways that may be associated with MSC characteristics, including apoptosis, MAPK, cell cycle, and Wnt signaling pathway. Signal-net analysis indicated that phospholipase C beta 4, filamin B beta, calcium/calmodulin-dependent protein kinase II gamma, and the ionotropic glutamate receptor, AMPA 1, had the highest betweenness centrality among significant genes in the differential gene profile network. A comparison between the coding-noncoding gene coexpression networks of PDLSCs and BMSCs identified chemokine (C-X-C motif) ligand 12 as a core regulatory factor in MSC biology. These results provided insight into the mechanisms underlying MSC biology.

1. Introduction

Stem cells are undifferentiated cells that can either self-renew or differentiate to produce mature progeny cells [1, 2]. The two major categories are embryonic and adult stem cells. Adult stem cells are undifferentiated cells found in specialized tissues and organs of adults. Compared to embryonic stem cells, adult stem cells that exist in various organs of the body are easily accessible, and their use is less controversial in terms of ethics [3, 4]. Mesenchymal stem cells (MSCs) have been identified as mesoderm-derived stromal cells that can differentiate into various mesoderm-type cell lineages. MSCs hold significant promise for tissue regeneration, due to their potential for self-renewal and multilineage differentiation [5–7]. Humans have abundant adult MSCs available for use in cell-based tissue engineering. MSCs from various tissues, including bone marrow, periosteum, skeletal muscle, and adipose tissue, have similar epitope profiles, but significant differences have been observed in MSC properties; that is, MSCs vary in their differentiation, proliferation, and migration potentials according to the tissue source [8–12]. Traditionally, bone-marrow-derived MSCs (BMSCs) have been studied for bone regeneration applications. BMSCs are a population of multipotent, nonhematopoietic marrow-derived cells that are easily expanded in culture and differentiate into cells with an osteogenic phenotype [13, 14]. BMSC transplantations have enhanced periodontal tissue regeneration and bone formation [15, 16]. Interestingly, Hu and colleagues investigated whether BMSCs might give rise to different types of epithelial cells, and they tested their potential for serving as a source of ameloblasts. Those results showed, for the first time, that BMSCs could be reprogrammed to become ameloblast-like cells [17]. Thus, BMSCs offered a novel approach for tooth-tissue engineering; they could be induced to become both mesenchymal and epithelial cells in tooth applications [17]. However, scientists disagree on whether BMSCs are ideal seeding cells for tooth engineering. Jing pointed out that the differentiation ability of BMSCs decreases significantly with increasing age of the donor [18]. In the past few decades, several new populations of MSCs have been isolated from dental and craniofacial tissues on the basis of their stem cell properties. These new populations included stem cells derived from the periodontal ligament (PDLSCs), from dental pulp, and from apical papilla, among others [19–24]. When transplanted into animals, these dental tissue-derived stem cells could generate bone/dentin-like mineralized tissue, and they were capable of repairing tooth defects and regenerating periodontal tissue [21, 25, 26]. In contrast to BMSCs, these cells were easily accessible, and they were more intimately associated with dental tissues [3]. Although dental tissue-derived MSCs and BMSCs are regulated by similar factors and share a common protein expression profile, these populations differ significantly in their proliferative ability and developmental potentials in vitro. Furthermore, importantly, they differ in their ability to develop into distinct tissues representative of the microenvironments from which they were derived in vivo. For example, BMSCs formed only bone tissue in the mouse model when treated in the same manner as the dental tissue-derived stem cells [19, 27]. However, the chondrogenic and adipogenic potentials of dental tissue-derived MSCs appeared to be weaker than those of BMSCs [22, 28]. Conversely, the neurogenicity of dental tissue-derived stem cells may be more potent than that of BMSCs, probably due to their neural crest origin [22, 28].

From the time that dental stem cells were first identified, they have been spotlighted in the dental tissue engineering field. Recently, numerous investigators have attempted to use these cells for dental tissue regeneration and assess their potential in preclinical applications [26, 29]. However, little is known about the characteristics of dental stem cells and the molecular mechanism underlying their diverse biological activities; thus, their potential application is restricted. Clues on the molecules that control MSC biology can be obtained by comparing molecular expression in MSCs with different biological activities. The development of microarray methods for large-scale analyses of mRNA gene expression has made it possible to search systematically for key molecules [30, 31]. With the introduction of these genome-wide research techniques, various groups have attempted to describe and compare the gene expression patterns of specialized adult stem cells [32–34]. Long, noncoding RNAs (lncRNAs) are transcribed RNA molecules longer than 200 nucleotides. LncRNAs have been shown to have comprehensive functions in both normal development and disease states [35]. Many studies have revealed that lncRNAs exert important roles in biological processes, including roles in cell differentiation, transcription, imprinting, chromatin modification, and others [36, 37]. Specifically, previous studies have demonstrated that lncRNAs are extremely important for controlling cell or tissue differentiation [38–40].

In this study, we investigated differences in lncRNA and mRNA expression profiles between PDLSCs and BMSCs with microarray assays and bioinformatics analyses. Our results provided useful information for elucidating the different mechanisms that govern MSCs derived from different tissues.

2. Materials and Methods

2.1. Cell Culture

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.” We collected impacted, third molars with immature roots from 3 healthy male patients (18–20 years old) under approved guidelines set by the Beijing Stomatological Hospital, Capital Medical University, after obtaining informed patient consent. Molars were removed, disinfected with 75% ethanol, and then washed with PBS. PDLSCs were isolated from each sample, cultured, and identified as previously described [21]. Briefly, PDLSCs were separated from the periodontal ligament in the middle one-third of the root. Then, the tissue was digested in a solution of 3 mg/mL collagenase type I (Worthington-Biochem, USA) and 4 mg/mL dispase (Roche, Germany) for 1 h at 37°C. Single-cell suspensions were obtained by passing the cells through a 70 μm strainer (Falcon, BD Labware, USA). Three separate PDLSC cultures were grown in a humidified, 5% CO₂ incubator at 37°C in alpha-modified Eagle’s medium (α-MEM; Invitrogen, California, USA) supplemented with 15% fetal bovine serum (FBS; Invitrogen), 2 mmol/L glutamine, 100 U/mL penicillin, and 100 μg/mL streptomycin (Invitrogen).

BMSCs derived from 18–20-year-old males () were obtained from Cyagen Biosciences (Guangzhou, China). Three separate BMSC cultures were grown in a humidified, 5% CO₂ incubator at 37°C, in Dulbecco’s MEM (Invitrogen), supplemented with 15% FBS (Invitrogen), 2 mmol/L glutamine, 100 U/mL penicillin, and 100 μg/mL streptomycin (Invitrogen). The culture medium was changed every 3 days. All MSCs were used in subsequent experiments after 3–5 passages.

2.2. Microarray Detection

MSCs were grown to 90% confluence; then, the BMSCs () and PDLSCs () were briefly rinsed with PBS, lysed, and total RNA was isolated with Trizol reagents (Invitrogen). rRNA was removed from total RNA and purified RNA was amplified and transcribed to produce fluorescent cRNA. Reverse transcription was performed along the entire length of the transcripts, without the 3′ bias, with a random priming method. cDNA was labeled and hybridized to the GeneChip Human Gene 2.0 ST Array (Affymetrix), according to the manufacturer’s protocol. After hybridization, washing, and staining, the chip was scanned according to the manufacturer’s instructions. Microarray experiments were performed at Genminix Informatic Ltd. (Shanghai, China), a microarray service certified by Affymetrix.

2.3. Real-Time RT-PCR Analysis

Real-time, reverse transcription-PCR (RT-PCR) was used to verify the differential expression of genes that were detected on the microarray. Total RNA was isolated from MSCs with Trizol reagents (Invitrogen). For real-time RT-PCR, 2 μg aliquots of RNA as template were combined with random hexamers and reverse transcriptase, according to the manufacturer’s protocol (Invitrogen). Real-time PCR reactions were performed with the QuantiTect SYBR Green PCR kit (Qiagen, Germany) and an iCycler iQ Multicolor Real-Time PCR Detection System. The relative level of gene expression was calculated with the method, as previously described [41]. Primers used for amplifying specific genes are shown in Table 1.

2.4. Bioinformatics Analysis

Differentially expressed genes were selected with the TwoClassDif method [9, 42, 43]. Gene ontology (GO) analysis was applied to analyze the main functions of differentially expressed genes. Gene ontology is the key functional classification method used at NCBI. GO can organize genes into hierarchical categories and uncover gene regulatory networks on the basis of biological processes and molecular functions [17, 44]. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, significantly changed pathways were identified and connected in a pathway network (Path-net), where connections were based on the relationship between these pathways. This approach was previously used to summarize the pathway interactions among genes that were differentially expressed under the influence of disease, and it revealed why certain pathways were activated [45].

Based on the GO and KEGG pathway analyses, we established an interactions repository (Signal-net) derived from KEGG to show the core genes that played an important role in this MSC gene network [46, 47]. To determine the interactions among genes, we constructed a coding-noncoding gene coexpression network (CNC network), which has also been called a gene coexpression network. This CNC network was based on a correlation analysis that evaluated associations between differentially expressed lncRNAs and mRNAs [45]. We calculated the Pearson correlation for each pair of genes and used the most significantly correlated pairs to construct the network [48]. The purpose of network structure analysis was to locate core regulatory factors (genes). In the network, the core regulatory factors were those connected to large numbers of adjacent genes, and, thus, they exhibited the greatest degrees of connectivity. In considering different networks, we evaluated the core regulatory factors by the degree of difference they showed in their roles in the PDLSC and BMSC networks [49], which was measured with the variable Diffk (difference in normalized connectivities).

2.5. Statistics

All statistical calculations were performed with SPSS10 statistical software. Statistical analyses included comparisons with the t-test, Fisher’s exact test, χ²test, and the Pearson correlation, as appropriate; values less than 0.05 were considered statistically significant.

3. Results

3.1. Comparison of lncRNA and mRNA Expression Profiles between PDLSCs and BMSCs

To reveal the molecular mechanisms underlying MSCs derived from different tissues, we screened the gene expression patterns in PDLSCs and BMSCs with the human GeneChip microarray method. Because we included only three samples in each group, we applied the RVM t-test, which can effectively raise the degrees of freedom in analyses of small sample sizes to filter the genes that were differentially expressed in PDLSCs and BMSCs. After determining significant differences and the false discovery rate (FDR) in the analysis, the differentially expressed genes were selected according to the value threshold. Hierarchical clustering showed systematic variations in the expression of lncRNAs and mRNAs between PDLSCs and BMSCs. From the microarray data, a comparison of lncRNA expression levels between PDLSCs and BMSCs identified an average of 970 lncRNAs that were significantly differentially expressed (see Supplementary Table 1 in Supplementary Material available online at http://dx.doi.org/10.1155/2014/317853); of those, 457 were upregulated and 513 were downregulated in the PDLSCs compared to the BMSCs. In addition, a total of 1,578 mRNAs were differentially expressed in the PDLSCs and BMSCs (Supplementary Table 2); of those, 862 were upregulated and 716 were downregulated in the PDLSCs compared to the BMSCs.

To confirm the reliability of the microarray data, we randomly selected six lncRNAs among the 970 differentially expressed lncRNAs and analyzed their expression with real-time RT-PCR. These data confirmed that, compared to BMSCs, PDLSCs showed increased expression of the lncRNAs coded as NR_045555, NR_027621, and NR_033651, and decreased expression of the lncRNAs coded as NR_037182, NR_037595, and XR_111050 (Figure 1). Similarly, we randomly selected six mRNAs among the 1,578 differentially expressed mRNAs and analyzed their expression with real-time RT-PCR. These data confirmed that the mRNAs BARX1, S100A4, WNT2B, and IGFBP5 were increased and that the mRNAs HOXA9 and HOXC8 were decreased in PDLSCs compared to BMSCs (Figure 2). The expression levels of these 12 genes were consistent with the microarray results; thus, these results confirmed the reliability of the microarray data.

Figure 1

Real-time RT-PCR results show differential lncRNA expression levels in stem cells derived from bone marrow (BMSCs) or periodontal ligament tissue (PDLSCs). The lncRNAs coded as NR_045555, NR_027621, and NR_033651 showed increased expression in PDLSCs, and the lncRNAs coded as NR_037182, XR_111050, and NR_037595 showed decreased expression in PDLSCs compared to BMSCs. GAPDH was used as an internal control. Student’s t-test was performed to determine statistical significance; all error bars represent s.d. ( tissue samples); .

3.2. Bioinformatics Analysis of BMSC and PDLSC Microarray Data

Next, we performed a bioinformatics analysis to discover the key factors that controlled MSC functions. First, a GO analysis was applied to analyze the main functions of the differentially expressed genes according to gene ontology, which is the key functional classification used by NCBI. According to the threshold, the analysis determined which GOs were significantly differently regulated between PDLSCs and BMSCs with a value and FDR < 0.05. The negative logarithm of the value (-LgP) was used to represent the correlation between gene expression and the relevant biological process. The GO analysis identified 166 genes that were significantly upregulated and 104 that were downregulated among all differentially expressed genes in PDLSCs (data not shown). The results clearly showed which important functions were involved with the differentially expressed genes. The top five upregulated GO functions (upGOs) were related to the response to the mitotic cell cycle, the phase of the mitotic cell cycle, mitotic prometaphase, the cell cycle checkpoint, and mitotic sister chromatid segregation (Supplementary Figure 1). The top five downregulated GO functions (downGOs) were related to the anterior/posterior pattern, embryonic skeletal system morphogenesis, signal transduction, cochlea morphogenesis, and blood vessel remodeling (Supplementary Figure 2).

Based on the KEGG database, we identified the pathways that mediated the functions of the differentially expressed genes. We identified a total of 67 pathways that showed significant differences due to differential gene expression; changes in 31 pathways involved upregulated genes and changes in 36 pathways involved downregulated genes (Supplementary Figures 3 and 4). We performed Path-net analysis to generate an interaction network that included these significantly changed pathways (Figure 3). The top 3 upregulated pathways were apoptosis, MAPK, and cell cycle signaling. The top 3 downregulated pathways were focal adhesion, Wnt, and adherens junction signaling. In addition, cytokine-cytokine receptor interactions and pathways related to cancer were up-/downregulated. These data suggested that these pathways may play key roles in the different core epigenetic mechanisms of PDLSCs and BMSCs.

Figure 3

The interaction network of significant pathways (Path-net) in stem cells. Pathways that were significantly different between PDLSCs and BMSCs were connected in a Path-net diagram to show the relationships between these pathways. The role of each pathway in the network was measured by counting its connections to upstream and downstream pathways, known as in-degree (upstream connections), out-degree (downstream connections), or degree (all connections). A high degree pathway indicated that it regulated or was regulated by many other pathways, which implied an important role in the signaling network. The circles represent the pathways; blue represents downregulated pathways, red represents upregulated pathways, and yellow represents up- and downregulated pathways. The lines indicate interactions between pathways.

We performed a Signal-net analysis to further investigate the global network, based on the significantly regulated GOs and pathways. With Signal-net, we screened important candidate genes involved in the differences between PDLSCs and BMSCs (Figure 4). In the Signal-net analysis, the genes are characterized by measuring their “betweenness centrality,” the number of times a node is located in the shortest path between 2 other nodes. This measure reflects the importance of a node in a graphic network relative to other nodes. The four most important differentially expressed genes were identified in the network (Supplementary Table 3); these were phospholipase C beta 4 (PLCβ4), filamin B beta (FLNB), calcium/calmodulin-dependent protein kinase II gamma (CAMK2G), and the ionotropic glutamate receptor, AMPA 1 (GRIA1).

Finally, we used a coding-noncoding gene coexpression (CNC) network to evaluate the interactions among genes and identify the core regulatory genes in the network. Based on our previous results, we built CNC networks to identify the interactions among the differentially expressed lncRNAs and mRNAs in PDLSCs and BMSCs [45]. We used 65 lncRNAs and 208 mRNAs to build the CNC network for PDLSCs and 75 lncRNAs and 187 mRNAs to build the network for BMSCs. In the CNC networks, each mRNA could correlate with one to tens of lncRNAs and vice versa. We used the CNC networks to implicate the interregulation of lncRNAs and mRNAs in the different molecular mechanisms of PDLSCs and BMSCs (Supplementary Figures 5 and 6). In the CNC network of PDLSCs, 17 genes showed a degree ≥ 59 and a clustering coefficient ≥ 0.6. This indicated that these genes, including 4 lncRNAs and 13 mRNAs (Table 2), played important roles in the network. In the CNC network of BMSCs, 20 mRNAs showed a degree ≥ 29 and a clustering coefficient ≥ 0.7. This indicated that (Table 3) these genes played important roles in the network. According to the Diffk values () for these networks, 16 genes (Table 4), including 2 lncRNAs and 14 mRNAs, showed different connectivities between PDLSCs and BMSCs, indicating that their roles were different in core pathways that governed MSC functions. The top three mRNAs were chemokine (C-X-C motif) ligand 12 (CXCL12), integrin alpha 2 (ITGA2, CD49B), and cell division cycle 20 homolog (CDC20), which were upregulated. The two lncRNAs identified () were FR020479 and FR191603; the former was downregulated and the latter was upregulated.

4. Discussion

The presence of different MSCs in dental and craniofacial tissues has encouraged clinical studies to investigate tissue regeneration in orofacial and periodontal regions [50, 51]. In the past few decades, MSC-mediated tissue regeneration has made surprising progress [25, 26, 52]. However, bone marrow has remained the principal source of MSCs for most preclinical and clinical applications. Interestingly, the MSCs from different anatomic locations possess diverse biological activities [8–12]. The challenge lies in identifying the specific genes that are associated with distinct MSC functions. To that end, in the present study, we identified lncRNAs and mRNAs that were differentially expressed in PDLSCs and BMSCs.

We identified 970 differentially expressed lncRNAs and 1,578 differentially expressed mRNAs in BMSCs and dental tissue-derived MSCs. This information may be useful for further studies on gene functions and regulation mechanisms in MSCs. Furthermore, we found that several of the upregulated genes in PDLSCs may be associated with PDLSC characteristics. For instance, BARX1, a transcription factor expressed in the mesenchyme of molar primordia, is involved in the regulation of tooth morphogenesis, in the development of tooth and craniofacial mesenchyme that originates from the neural crest [53–55], and possibly, in the regulation of MSC differentiation.

To identify the key factors that regulated MSC functions, we applied bioinformatics analyses to classify the microarray data. The GO analysis revealed specific functional pathways that were enriched in the genes responsible for the divergent features of PDLSCs and BMSCs. These differentially expressed genes were subsequently organized into hierarchical categories based on pertinent biological processes. A high degree pathway interacted with a high number of other pathways, which implied an important role in cell biological features. Further pathway analyses indicated that apoptosis, MAPK, cytokine-cytokine receptor interaction, focal adhesion, pathways in cancer, Wnt, cell cycle, and adherens junctions signaling pathways were involved in the diverse biological activities of PDLSCs and BMSCs. It is well known that these pathways play an important role in regulating cellular apoptosis, survival, and differentiation.

To identify important genes involved in the different epigenetic mechanisms of PDLSCs and BMSCs, we performed Signal-net analysis on the significantly regulated GOs and pathways. This analysis revealed that PLCβ4, FLNB, CAMK2G, and GRIA1 exhibited the most betweenness centrality. PLCβ4 and CAMK2G were upregulated in PDLSCs. It was reported that PLCβ4 was highly expressed in the retina and the cerebellum, where calcium plays an important role in the transduction of extracellular signals [56–58]. Moreover, CAMK2G is activated by intracellular calcium/calmodulin [59]. Thus, the Signal-net analysis results suggested that these genes were important in calcium-sensitive signaling cascades that regulate cell function. In addition, FLNB regulates intracellular communication and signaling by linking the protein actin to the cell membrane. This activity allows direct communication between the cell membrane and the cytoskeletal network, which provides a means to control and guide proper skeletal development [60, 61].

The CNC network comparisons indicated that CXCL12 was a core regulatory factor, which may be involved in the diverse biological activities of PDLSCs and BMSCs. CXCL12, also known as stromal cell-derived factor-1, stimulates migration by rearranging the actin cytoskeleton, increasing focal adhesion, and stimulating matrix metalloproteinase production in MSCs [62, 63]. Thus, CXCL12 can recruit MSC to participate in the regeneration of injured tissues [64]. Presumably, MSC migration is mediated through an intracellular pathway, for example, the MAPK/ERK signaling pathways [62]. Our results were consistent with previous reports and may also be applicable to the differentiation mechanisms previously described in MSCs.

Additionally, we identified some lncRNAs that were differentially expressed in PDLSCs and BMSCs, for example, FR020479 and FR191603. Previous studies demonstrated that lncRNAs may function by controlling the transcriptional regulation of neighboring coding genes [65, 66]. Identifying differentially expressed nearby coding mRNAs may enhance our understanding of the function of lncRNAs in MSCs. However, further studies must be performed to investigate that hypothesis.

5. Conclusion

This study provided comprehensive profiles of mRNA and lncRNA expression in PDLSCs and BMSCs, two tissue-derived MSCs. In addition, potential regulatory mechanisms were identified with bioinformatics analyses. Although more studies are required to demonstrate the precise role and mechanisms of these lncRNAs and mRNAs, the genomic data we identified with microarray analyses may increase our understanding of MSC biology.

Conflict of Interests

The authors declare that there is no potential conflict of interests.

Acknowledgments

This work was supported by Grants from the National Natural Science Foundation of China (81271101 to Zhipeng Fan, 81271100 to Rui Dong, and 81170931 to Juan Du), the Program for New Century Excellent Talents in the University (NCET-12-0611 to Zhipeng Fan), High-Level Talents of the Beijing Health System (2013-3-035 to Zhipeng Fan and 2013-3-034 to Rui Dong), the Beijing Natural Science Foundation (7112057 to Rui Dong), and the Beijing Nova Program (2011083 to Rui Dong).

Supplementary Materials

Expression profiles of lncRNA and mRNA between bone marrow stem cells (BMSCs) and periodontal ligament stem cells (PDLSCs) were investigated with lncRNA microarray assays and bioinformatics analysis. In PDLSCs, 970 lncRNAs that were significantly differentially expressed compared to BMSCs (Supplementary Table 1). Furthermore, 1,578 mRNAs were differentially expressed in the PDLSCs and BMSCs (Supplementary Table 2). For bioinformatics analysis, the results of the GO analysis showed which important functions were involved with the differentially expressed genes, including the top upregulated and downregulated GO functions (upGOs and downGOs) (Supplementary Figures 1 and 2). Then based on the KEGG database, there were identified 67 pathways that showed significant differences due to differential gene expression (Supplementary Figures 3 and 4), which play key roles in the different core epigenetic mechanisms of PDLSCs and BMSCs. To further investigate the global network, differentially expressed genes were identified by a Signal-net analysis (Supplementary Table 3). Finally, coding-noncoding gene coexpression (CNC) networks were used to implicate the interregulation of lncRNAs and mRNAs in the different molecular mechanisms of PDLSCs and BMSCs (Supplementary Figures 5 and 6).

Supplementary Materials

References

J. Cai, M. L. Weiss, and M. S. Rao, “In search of “stemness’,” Experimental Hematology, vol. 32, no. 7, pp. 585–598, 2004.
View at: Publisher Site | Google Scholar
A. J. Wagers and I. L. Weissman, “Plasticity of adult stem cells,” Cell, vol. 116, no. 5, pp. 639–648, 2004.
View at: Publisher Site | Google Scholar
C. Morsczeck, G. Schmalz, T. E. Reichert, F. Völlner, K. Galler, and O. Driemel, “Somatic stem cells for regenerative dentistry,” Clinical Oral Investigations, vol. 12, no. 2, pp. 113–118, 2008.
View at: Publisher Site | Google Scholar
A. Mauron and M. E. Jaconi, “Stem cell science: current ethical and policy issues,” Clinical Pharmacology and Therapeutics, vol. 82, no. 3, pp. 330–333, 2007.
View at: Publisher Site | Google Scholar
B. N. Jahagirdar and C. M. Verfaillie, “Multipotent adult progenitor cell and stem cell plasticity,” Stem Cell Reviews, vol. 1, no. 1, pp. 53–59, 2005.
View at: Google Scholar
D. G. Phinney and D. J. Prockop, “Concise review: mesenchymal stem/multipotent stromal cells: the state of transdifferentiation and modes of tissue repair—current views,” Stem Cells, vol. 25, no. 11, pp. 2896–2902, 2007.
View at: Publisher Site | Google Scholar
J. Feng, A. Mantesso, C. de Bari, A. Nishiyama, and P. T. Sharp, “Dual origin of mesenchymal stem cells contributing to organ growth and repair,” Proceedings of the National Academy of Sciences of the United States of America, vol. 108, no. 16, pp. 6503–6508, 2011.
View at: Publisher Site | Google Scholar
H. Yang, L. N. Gao, Y. An et al., “Comparison of mesenchymal stem cells derived from gingival tissue and periodontal ligament in different incubation conditions,” Biomaterials, vol. 34, no. 29, pp. 7033–7047, 2013.
View at: Publisher Site | Google Scholar
S. Yu, J. Long, J. Yu et al., “Analysis of differentiation potentials and gene expression profiles of mesenchymal stem cells derived from periodontal ligament and Wharton's jelly of the umbilical cord,” Cells Tissues Organs, vol. 197, no. 3, pp. 209–223, 2013.
View at: Publisher Site | Google Scholar
O. Hayashi, Y. Katsube, M. Hirose, H. Ohgushi, and H. Ito, “Comparison of osteogenic ability of rat mesenchymal stem cells from bone marrow, periosteum, and adipose tissue,” Calcified Tissue International, vol. 82, no. 3, pp. 238–247, 2008.
View at: Publisher Site | Google Scholar
Y. Sakaguchi, I. Sekiya, K. Yagishita, and T. Muneta, “Comparison of human stem cells derived from various mesenchymal tissues: superiority of synovium as a cell source,” Arthritis and Rheumatism, vol. 52, no. 8, pp. 2521–2529, 2005.
View at: Publisher Site | Google Scholar
Y. Zhang, C.-D. Li, X.-X. Jiang, H.-L. Li, P.-H. Tang, and N. Mao, “Comparison of mesenchymal stem cells from human placenta and bone marrow,” Chinese Medical Journal, vol. 117, no. 6, pp. 882–887, 2004.
View at: Google Scholar
A. I. Caplan, “Mesenchymal stem cells: cell-based reconstructive therapy in orthopedics,” Tissue Engineering, vol. 11, no. 7-8, pp. 1198–1211, 2005.
View at: Publisher Site | Google Scholar
S. P. Bruder, K. H. Kraus, V. M. Goldberg, and S. Kadiyala, “The effect of implants loaded with autologous mesenchymal stem cells on the healing of canine segmental bone defects,” Journal of Bone and Joint Surgery A, vol. 80, no. 7, pp. 985–996, 1998.
View at: Google Scholar
Y. M. Lee, Y. J. Seol, Y. T. Lim et al., “Tissue-engineered growth of bone by marrow cell transplantation using porous calcium metaphosphate matrices,” Journal of Biomedical Materials Research, vol. 54, pp. 216–223, 2001.
View at: Google Scholar
H. Kawaguchi, A. Hirachi, N. Hasegawa et al., “Enhancement of periodontal tissue regeneration by transplantation of bone marrow mesenchymal stem cells,” Journal of Periodontology, vol. 75, no. 9, pp. 1281–1287, 2004.
View at: Publisher Site | Google Scholar
B. Hu, F. Unda, S. Bopp-Kuchler et al., “Bone marrow cells can give rise to ameloblast-like cells,” Journal of Dental Research, vol. 85, no. 5, pp. 416–421, 2006.
View at: Publisher Site | Google Scholar
W. Jing, L. Wu, Y. Lin, L. Liu, W. Tang, and W. Tian, “Odontogenic differentiation of adipose-derived stem cells for tooth regeneration: necessity, possibility, and strategy,” Medical Hypotheses, vol. 70, no. 3, pp. 540–542, 2008.
View at: Publisher Site | Google Scholar
S. Gronthos, M. Mankani, J. Brahim, P. G. Robey, and S. Shi, “Postnatal human dental pulp stem cells (DPSCs) in vitro and in vivo,” Proceedings of the National Academy of Sciences of the United States of America, vol. 97, no. 25, pp. 13625–13630, 2000.
View at: Publisher Site | Google Scholar
M. Miura, S. Gronthos, M. Zhao et al., “SHED: stem cells from human exfoliated deciduous teeth,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 10, pp. 5807–5812, 2003.
View at: Publisher Site | Google Scholar
B.-M. Seo, M. Miura, S. Gronthos et al., “Investigation of multipotent postnatal stem cells from human periodontal ligament,” The Lancet, vol. 364, no. 9429, pp. 149–155, 2004.
View at: Publisher Site | Google Scholar
W. Sonoyama, Y. Liu, T. Yamaza et al., “Characterization of the apical papilla and its residing stem cells from human immature permanent teeth: a pilot study,” Journal of Endodontics, vol. 34, no. 2, pp. 166–171, 2008.
View at: Publisher Site | Google Scholar
R. d'Aquino, A. Graziano, M. Sampaolesi et al., “Human postnatal dental pulp cells co-differentiate into osteoblasts and endotheliocytes: a pivotal synergy leading to adult bone tissue formation,” Cell Death and Differentiation, vol. 14, no. 6, pp. 1162–1171, 2007.
View at: Publisher Site | Google Scholar
G. Laino, R. D'Aquino, A. Graziano et al., “A new population of human adult dental pulp stem cells: a useful source of living autologous fibrous bone tissue (LAB),” Journal of Bone and Mineral Research, vol. 20, no. 8, pp. 1394–1402, 2005.
View at: Publisher Site | Google Scholar
G. Ding, Y. Liu, W. Wang et al., “Allogeneic periodontal ligament stem cell therapy for periodontitis in swine,” Stem Cells, vol. 28, no. 10, pp. 1829–1838, 2010.
View at: Publisher Site | Google Scholar
Y. Liu, Y. Zheng, G. Ding et al., “Periodontal ligament stem cell-mediated treatment for periodontitis in miniature swine,” Stem Cells, vol. 26, no. 4, pp. 1065–1073, 2008.
View at: Publisher Site | Google Scholar
S. Batouli, M. Miura, J. Brahim et al., “Comparison of stem-cell-mediated osteogenesis and dentinogenesis,” Journal of Dental Research, vol. 82, no. 12, pp. 976–981, 2003.
View at: Google Scholar
W. Zhang, X. F. Walboomers, S. Shi, M. Fan, and J. A. Jansen, “Multilineage differentiation potential of stem cells derived from human dental pulp after cryopreservation,” Tissue Engineering, vol. 12, no. 10, pp. 2813–2823, 2006.
View at: Publisher Site | Google Scholar
S.-H. Kim, K.-H. Kim, B.-M. Seo et al., “Alveolar bone regeneration by transplantation of periodontal ligament stem cells and bone marrow stem cells in a canine peri-implant defect model: a pilot study,” Journal of Periodontology, vol. 80, no. 11, pp. 1815–1823, 2009.
View at: Publisher Site | Google Scholar
D. J. Duggan, M. Bittner, Y. Chen, P. Meltzer, and J. M. Trent, “Expression profiling using cDNA microarrays,” Nature Genetics, vol. 21, no. 1, pp. 10–14, 1999.
View at: Publisher Site | Google Scholar
P. O. Brown and D. Botstein, “Exploring the new world of the genome with DNA microarrays,” Nature Genetics, vol. 21, no. 1, pp. 33–37, 1999.
View at: Publisher Site | Google Scholar
S. Shi, P. G. Robey, and S. Gronthos, “Comparison of human dental pulp and bone marrow stromal stem cells by cDNA microarray analysis,” Bone, vol. 29, no. 6, pp. 532–539, 2001.
View at: Publisher Site | Google Scholar
N. Inomata, H. Tomita, Z. Ikezawa, and H. Saito, “Differential gene expression profile between cord blood progenitor-derived and adult progenitor-derived human mast cells,” Immunology Letters, vol. 98, no. 2, pp. 265–271, 2005.
View at: Publisher Site | Google Scholar
T. Schilling, R. Küffner, L. Klein-Hitpass, R. Zimmer, F. Jakob, and N. Schütze, “Microarray analyses of transdifferentiated mesenchymal stem cells,” Journal of Cellular Biochemistry, vol. 103, no. 2, pp. 413–433, 2008.
View at: Publisher Site | Google Scholar
C. P. Ponting, P. L. Oliver, and W. Reik, “Evolution and functions of long noncoding RNAs,” Cell, vol. 136, no. 4, pp. 629–641, 2009.
View at: Publisher Site | Google Scholar
J. T. Kung, D. Colognori, and J. T. Lee, “Long noncoding RNAs: past, present, and future,” Genetics, vol. 193, no. 3, pp. 651–669, 2013.
View at: Publisher Site | Google Scholar
S. Ghosal, S. Das, and J. Chakrabarti, “Long noncoding RNAs: new players in the molecular mechanism for maintenance and differentiation of pluripotent stem cells,” Stem Cells and Development, vol. 22, no. 16, pp. 2240–2253, 2013.
View at: Publisher Site | Google Scholar
M. Cesana, D. Cacchiarelli, I. Legnini et al., “A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA,” Cell, vol. 147, no. 2, pp. 358–369, 2011.
View at: Publisher Site | Google Scholar
W. Hu, J. R. Alvarez-Dominguez, and H. F. Lodish, “Regulation of mammalian cell differentiation by long non-coding RNAs,” EMBO Reports, vol. 13, no. 11, pp. 971–983, 2012.
View at: Publisher Site | Google Scholar
M. Kretz, D. E. Webster, R. J. Flockhart et al., “Suppression of progenitor differentiation requires the long noncoding RNA ANCR,” Genes and Development, vol. 26, no. 4, pp. 338–343, 2012.
View at: Publisher Site | Google Scholar
S. A. Bustin, “Absolute quantification of mrna using real-time reverse transcription polymerase chain reaction assays,” Journal of Molecular Endocrinology, vol. 25, no. 2, pp. 169–193, 2000.
View at: Google Scholar
G. W. Wright and R. M. Simon, “A random variance model for detection of differential gene expression in small microarray experiments,” Bioinformatics, vol. 19, no. 18, pp. 2448–2455, 2003.
View at: Publisher Site | Google Scholar
R. Clarke, H. W. Ressom, A. Wang et al., “The properties of high-dimensional data spaces: implications for exploring gene and protein expression data,” Nature Reviews Cancer, vol. 8, no. 1, pp. 37–49, 2008.
View at: Publisher Site | Google Scholar
M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene ontology: tool for the unification of biology,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000.
View at: Publisher Site | Google Scholar
M. A. Pujana, J.-D. J. Han, L. M. Starita et al., “Network modeling links breast cancer susceptibility and centrosome dysfunction,” Nature Genetics, vol. 39, no. 11, pp. 1338–1349, 2007.
View at: Publisher Site | Google Scholar
C. Li and H. Li, “Network-constrained regularization and variable selection for analysis of genomic data,” Bioinformatics, vol. 24, no. 9, pp. 1175–1182, 2008.
View at: Publisher Site | Google Scholar
J. D. Zhang and S. Wiemann, “KEGGgraph: a graph approach to KEGG PATHWAY in R and bioconductor,” Bioinformatics, vol. 25, no. 11, pp. 1470–1471, 2009.
View at: Publisher Site | Google Scholar
C. Prieto, A. Risueño, C. Fontanillo, and J. de Las Rivas, “Human gene coexpression landscape: confident network derived from tissue transcriptomic profiles,” PLoS ONE, vol. 3, no. 12, Article ID e3911, 2008.
View at: Publisher Site | Google Scholar
M. R. J. Carlson, B. Zhang, Z. Fang, P. S. Mischel, S. Horvath, and S. F. Nelson, “Gene connectivity, function, and sequence conservation: predictions from modular yeast co-expression networks,” BMC Genomics, vol. 7, article 40, 2006.
View at: Publisher Site | Google Scholar
K. G. Silvério, B. B. Benatti, M. Z. Casati, E. A. Sallum, and F. H. Nociti Jr., “Stem cells: potential therapeutics for periodontal regeneration,” Stem Cell Reviews, vol. 4, no. 1, pp. 13–19, 2008.
View at: Publisher Site | Google Scholar
L. Rimondini and S. Mele, “Stem cell technologies for tissue regeneration in dentistry,” Minerva Stomatologica, vol. 58, no. 10, pp. 483–500, 2009.
View at: Google Scholar
A. Giuliani, A. Manescu, M. Langer et al., “Three years after transplants in human mandibles, histological and in-line holotomography revealed that stem cells regenerated a compact rather than a spongy bone: biological and clinical implications,” Stem Cells Translational Medicine, vol. 2, no. 4, pp. 316–324, 2013.
View at: Publisher Site | Google Scholar
S. M. Sperber and I. B. Dawid, “barx1 is necessary for ectomesenchyme proliferation and osteochondroprogenitor condensation in the zebrafish pharyngeal arches,” Developmental Biology, vol. 321, no. 1, pp. 101–110, 2008.
View at: Publisher Site | Google Scholar
I. Miletich, G. Buchner, and P. T. Sharpe, “Barx1 and evolutionary changes in feeding,” Journal of Anatomy, vol. 207, no. 5, pp. 619–622, 2005.
View at: Google Scholar
T. A. Mitsiadis, M.-L. Mucchielli, S. Raffo, J.-P. Proust, P. Koopman, and C. Goridis, “Expression of the transcription factors Otlx2, Barx1 and Sox9 during mouse odontogenesis,” European Journal of Oral Sciences, vol. 106, supplement 1, pp. 112–116, 1998.
View at: Google Scholar
C. W. Taylor, “Controlling calcium entry,” Cell, vol. 111, no. 6, pp. 767–769, 2002.
View at: Publisher Site | Google Scholar
Y. Kido, C. T. Gordon, S. Sakazume et al., “Further characterization of atypical features in auriculocondylar syndrome caused by recessive PLCB4 mutations,” American Journal of Medical Genetics A, vol. 161, no. 9, pp. 2339–2346, 2013.
View at: Publisher Site | Google Scholar
A. M. Lyon and J. J. Tesmer, “Structural insights into phospholipase C-beta function,” Molecular Pharmacology, vol. 84, no. 4, pp. 488–500, 2013.
View at: Publisher Site | Google Scholar
M. Mattingsdal, A. A. Brown, S. Djurovic et al., “Pathway analysis of genetic markers associated with a functional MRI faces paradigm implicates polymorphisms in calcium responsive pathways,” Neuroimage, vol. 70, pp. 143–149, 2013.
View at: Publisher Site | Google Scholar
J. Lu, G. Lian, R. Lenkinski et al., “Filamin B mutations cause chondrocyte defects in skeletal development,” Human Molecular Genetics, vol. 16, no. 14, pp. 1661–1675, 2007.
View at: Publisher Site | Google Scholar
T. P. Stossel, J. Condeelis, L. Cooley et al., “Filamins as integrators of cell mechanics and signalling,” Nature Reviews Molecular Cell Biology, vol. 2, no. 2, pp. 138–145, 2001.
View at: Publisher Site | Google Scholar
H. Gao, W. Priebe, J. Glod, and D. Banerjee, “Activation of signal transducers and activators of transcription 3 and focal adhesion kinase by stromal cell-derived factor 1 is required for migration of human mesenchymal stem cells in response to tumor cell-conditioned medium,” Stem Cells, vol. 27, no. 4, pp. 857–865, 2009.
View at: Publisher Site | Google Scholar
J. Y. Hsieh, T. S. Huang, S. M. Cheng et al., “miR-146a-5p circuitry uncouples cell proliferation and migration, but not differentiation, in human mesenchymal stem cells,” Nucleic Acids Research, vol. 41, no. 21, pp. 9753–9763, 2013.
View at: Publisher Site | Google Scholar
T. Ponomaryov, A. Peled, I. Petit et al., “Induction of the chemokine stromal-derived factor-1 following DNA damage improves human stem cell function,” Journal of Clinical Investigation, vol. 106, no. 11, pp. 1331–1339, 2000.
View at: Google Scholar
T. Hung and H. Y. Chang, “Long noncoding RNA in genome regulation: prospects and mechanisms,” RNA Biology, vol. 7, no. 5, pp. 582–585, 2010.
View at: Google Scholar
H. Ma, Y. Hao, X. Dong et al., “Molecular mechanisms and function prediction of long noncoding RNA,” TheScientificWorldJOURNAL, vol. 2012, Article ID 541786, 11 pages, 2012.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2014 Rui Dong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

3214

Downloads

2082

Citations