Printer Friendly

Differential Expression Profiling of Long Noncoding RNA and mRNA during Osteoblast Differentiation in Mouse.

1. Introduction

Ossification is a tightly regulated process which is performed by specialized cells called osteoblasts differentiated from mesenchymal progenitors. The osteoblast differentiation process is regulated by several key factors and signaling pathways. Runt-related transcription factor 2 (Runx2) is the master switch in the commitment of mesenchymal progenitors to osteoblast lineage [1]. Runx2 is affected by several upstream regulators such as the Wnt/Notch system, Sox9, Msx2, and hedgehog signaling as well as by cofactors such as Osx and Atf4 [1-4]. A few paracrine and endocrine factors, including bone morphogenetic proteins (BMP) and parathyroid hormone, serve as coactivators. Vitamin D and histone deacetylase enzymes coordinate this process more finely [1].

Long noncoding RNAs (lncRNAs) are a class of RNA transcripts longer than 200 nucleotides, lacking open reading frames and protein-coding possibilities [5]. Tens of thousands of lncRNAs have been identified in mammalian genomes in recent decades [6]. Recently, in many studies, lncRNA has emerged as an important regulator in a variety of biological processes, such as epigenetic regulation, chromatin remodeling, genomic imprinting, transcriptional control, and pre-/posttranslational mRNA processing [7-11]. In terms of osteogenesis, several lncRNAs have been found to act as key regulators. One such example is maternal expression gene 3 (MEG3) regulating the expression of Bmp4, Runx2, and Osx in human mesenchymal stem cells [12, 13]. Antidifferentiation ncRNA (ANCR) inhibits Runx2 expression in association with the enhancer of zeste homolog 2 (EZH2); thus, downregulation of ANCR promotes osteoblast differentiation through modulation of EZH2/Runx2 [14].

In order to fully understand the lncRNA biology including its role in osteogenesis, it is necessary to characterize the expression pattern of lncRNA during osteoblast differentiation. In this study, we analyzed RNA-sequencing datasets obtained at nine different time points in the osteoblast differentiation of preosteoblasts isolated from neonatal mouse calvaria, using various bioinformatic approaches. We focused on identifying differentially expressed lncRNAs and mRNAs throughout the process and finding potential markers that exhibited significant changes in expression during the osteoblast differentiation.

2. Materials and Methods

2.1. RNA-Seq Data Processing and Analysis of Differential Gene Expression. We analyzed the RNA-Seq data generated by Kemp et al. [15] (GSE54461, nine time points: 2, 4, 6, 8, 10, 12, 14, 16, and 18 days of osteoblast differentiation). All sequencing reads were aligned to the mouse genome reference (mm10) using the GSNAP alignment tool [16]. Ensembl release 74 annotations were used to measure gene expression. In the case of lncRNAs, NONCODE v4 ( annotations were used. Since NONCODE v4 annotations were based on mm9 mouse genome assembly, the positions of lncRNAs of NONCODE were converted from mm9 to mm10 using the LiftOver utility in UCSC. To assess gene expression, RPKM (reads per kilobase of exon per million mapped reads) values were calculated [17]. Hierarchical clustering of genes expressed in samples at the nine time points was performed using the flashClust R package. To establish the differences in gene expression patterns among nine time points of osteoblast differentiation, we performed differential gene expression analysis using the R package DESeq [18]. The false discovery rate (FDR) was controlled by adjusting p values using the Benjamini-Hochberg algorithm. Differentially expressed genes were defined as those with FDR less than 5% with an absolute value of fold change [greater than or equal to] 2. Similar differential expression analyses were performed for lncRNAs.

2.2. Time Series Analysis of Differential Gene Expression. R package DESeq [18] and edgeR [19] were employed to identify genes that were differentially expressed across the differentiation time period to a significant extent, designated as time series genes. We selected time series genes that displayed significant differential expression with FDR < 5% in both DESeq and edgeR, absolute fold change [greater than or equal to] 2 (between day 2 and at least one other time point), and maximum RPKM [greater than or equal to] 3 across the time series. Details of the methods for each package are described in the supplementary methods. Similarly, time series expression tests were conducted for lncRNAs. Using the RPKM of each time series gene and lncRNA, principal component analysis (PCA) was performed with the aid of the "prcomp" module in R.

2.3. Time Series Gene Clustering. The time series genes identified were clustered using the R package Mfuzz [20] that performs soft clustering based on the fuzzy c-means algorithm. The advantage of soft clustering is that the algorithm clearly reflects the strength of association of an individual gene with a cluster. Average RPKM values (triplicates at each time point) of individual genes were employed as input values for Mfuzz clustering. The number of clusters was set to 12 and the fuzzifier coefficient, M, to 1.5. Heat maps of the clusters were drawn using the R module "heatmap.2" in the "gplots" package [21].

2.4. Weighted Gene Coexpression Network Analysis. Gene coexpression network analysis was performed using the R package "WGCNA" [22]. Details of the methods for constructing gene coexpression network analysis are described in the supplementary methods.

2.5. Gene Ontology (GO) Term and KEGG Pathway Enrichment Analysis. GO terms of each gene were obtained from Ensembl BioMart and KEGG pathways from the KEGG PATHWAY database. Details of the procedure to conduct the enrichment analysis are described in the supplementary methods.

2.6. Analysis of Motif Enrichment. MEME Suite 4.12.20 [23] and the HOmo sapiens COmprehensive MOdel COllection (HOCOMOCO) v11 mouse transcription factor database [24] were used for the identification of motifs in the promoter regions on the lncRNAs identified as potential markers. Details of the procedure to perform the motif enrichment analysis are described in the supplementary methods.

3. Results and Discussion

For convenience, the analysis results of mRNA and lncRNAs were described separately. The overall expression pattern was described first, followed by time series analysis, generation of modules by weighted gene coexpression network analysis, and identification of markers for osteoblast differentiation (Supplementary Figure 1).

3.1. Expression Profile of mRNA during Osteoblast Differentiation. A total of 12 to 33 million reads (at each time point) were processed and mapped against mouse genome reference (mm10) sequence, and uniquely mapped reads (89.01-90.74% of the total reads) were used for further analysis (Table 1). Based on Ensembl release 74 annotations, a total of 28,582 genes had at least one read for whole RNA-Seq datasets. By hierarchical clustering [25], two distinct clusters were generated for the early (2, 4, 6, and 8 days) and late (10, 12, 14, 16, and 18 days) differentiation stages (Figure 1(a)). As expected, the number of differentially expressed genes increased with time (Figure 1(b)).

Comparing the expression of individual genes at different time points using all possible combinations, 46% of the mapped genes (13,130 out of 28,582) were differentially expressed over time (fold-change difference [greater than or equal to] 2, FDR < 0.05). For example, a comparison of gene expression between the two time points selected in the early and late stages (day 2 and day 18) is shown in Figure 1(c). A total of 7238 genes were differentially expressed between days 2 and 18 (4431 upregulated and 2807 downregulated on day 18). GO and KEGG pathway analyses revealed that genes upregulated on day 18 were strongly associated with osteogenesis processes, such as "extracellular matrix binding," "positive regulation of bone mineralization," "collagen type I," "elevation of cytosolic calcium ion concentration," "bone mineralization involved in bone maturation," and "positive regulation of Wnt receptor signaling pathway" (Supplementary Table 1). On the other hand, genes downregulated on day 18 were associated with cell proliferation processes, such as "cell division," "G1/S transition of mitotic cell cycle," "regulation of cell cycle," "M phase of mitotic cell cycle," and "DNA replication and response to DNA damage stimulus" (Supplementary Table 2).

3.2. Dynamic Changes of mRNA Expression over Time of Osteoblast Differentiation. During the differentiation process, several genes are dynamically expressed via complex regulatory mechanisms. To characterize temporal gene expression changes and patterns, we identified genes that were differentially expressed across the time course of osteoblast differentiation. A number of rules were applied to establish significant differential expression of genes at different time points: (1) statistical significance of temporal gene expression changes assessed via edgeR and DESeq methods (FDR < 0.05), (2) absolute fold change between day 2 and at least one other time point [greater than or equal to] 2, and (3) maximum RPKM across the time series [greater than or equal to] 3.

In total, 4058 genes were significantly differentially expressed during osteoblast differentiation. They were associated with osteogenesis processes, such as "extracellular matrix organization," "osteoblast proliferation," "osteoblast differentiation," "positive regulation of osteoblast differentiation," and "actin cytoskeleton and bone mineralization." Principal component analysis (PCA) of these genes revealed that 95% of the variations in gene expression could be explained by the first two principal components (PCs) (Supplementary Figure 2) and that the PCs dominantly separate the datasets according to differentiation stage (days 2-8 versus days 10-18) (Figure 1(d)). To evaluate the osteoblast differentiation stages in the current dataset, we have checked the expression level of the twelve mature osteoblast markers (Bglap, Ibsp, Dmp1, Col13a1, Pthr1, Lifr, Bambi, Dlx3, Hnf1a, Phex, Ptgis, and Cdo1) described in Kalajzic et al. [26]. The nine genes, except three genes (Pthr1, Bambi, and Hnf1a), exhibited very similar expression patterns and showed the highest expression at day 18, which indicate that the cells at day 18 were mature osteoblasts (Supplementary Figure 3).

Osteoblasts undergo several stages before maturation and mineralization of bone matrix. Differentiation of osteoblasts, both in vitro and in vivo, is divided into three stages: (1) cell proliferation, (2) matrix maturation, and (3) matrix mineralization [27]. Here, we examined differences in gene expression patterns at the early and late osteoblast differentiation stages. Based on the current study, gene expression profile analysis demonstrated that osteoblast differentiation could be clearly subdivided into two stages. The early stage corresponded to "cell proliferation" based on the Stein and Jane definition [27]. However, we could not separate the late differentiation stage into "matrix maturation" and "matrix mineralization" stages based on the hierarchical clustering and PCA analysis. In addition, the division of osteoblast differentiation into two stages (early and late stages) is observed in the time series gene clustering analysis and weighted gene coexpression network analysis described below.

3.3. Clustering Analysis of Time Series Genes. Differentially expressed genes identified by time series analysis were clustered according to their expression profiles (RPKM values) using the fuzzy c-means algorithm implemented in R Mfuzz package. Genes with similar time-specific expression patterns were clustered into 12 groups, each containing 191-570 genes (Figure 2(a)). Eight among the 12 clusters (clusters 2, 3, 4, 5, 6, 8, 11, and 12) showed high expression patterns at the early osteoblast differentiation stage. Functional enrichment analysis for each cluster showed that these 8 clusters were highly enriched for genes related to cell proliferation processes. On the other hand, genes included in the remaining 4 clusters (1, 7, 9, and 10) showed high expression patterns at the late osteoblast differentiation stage. These clusters included several genes associated with functional roles in osteogenesis, such as "extracellular matrix organization," "osteoblast differentiation," "metabolic process," "collagen fibril organization," "Wnt-protein binding," and "skeletal system development."

We focused on clusters 1, 7, 9, and 10, each containing 511, 570, 243, and 297 genes with high expression patterns in the late differentiation stage (Figure 2(b)). Assessment of individual genes within the clusters showed that several are involved in osteogenesis. For example, cluster 1 contained Fgfr2 that plays an essential role in skeletal development [28], cluster 7 contained Bmp4 that promotes formation of the bone and cartilage by stimulating differentiation processes of osteoblasts [29], and cluster 9 contained Fgf18 that has been shown to stimulate the proliferation of cultured mouse primary osteoblasts [30]. Expression levels of Fgfr2, Bmp4, and Fgf18 were particularly high between 12 and 18 days of the osteoblast differentiation process.

We determined the common function among the three clusters through GO and KEGG enrichment analysis. As shown in Table 2, the statistically significant GO terms in the clusters 1, 7, and 9 were "collagen metabolism," "subset of extracellular matrix," "regulation of osteoblast differentiation," and "Wnt signaling" which is known to be involved in the regulation of osteoblast lineage cells [31]. The enriched pathways in the KEGG analysis in the clusters 1, 7, and 9, were "cell adhesion," "cytokinecytokine receptor interaction," and "PPAR signaling pathway" which is functionally associated with bone metabolism [32]. In the case of the cluster 10, "extracellular matrix" is the only statistically significant GO/KEGG term. This indicates that the genes belonging to the three clusters (1, 7, and 9) are highly correlated with the late osteoblast stage but the cluster 10 is less correlated with the late osteoblast stage. Hence, together with the coexpression network analysis described below, we used the genes belonging to the three clusters (1, 7, and 9 except 10) to find biomarkers of osteoblast differentiation.

3.4. Construction of a Gene Coexpression Network. Gene coexpression patterns may provide information on gene networks or pathways related to different biological phenomena. We constructed a coexpression network of genes using time series gene expression data. The weighted gene coexpression network analysis (WGCNA) generated a total of 16 network modules (Figure 3(a)). The eigengenes for all 16 modules were calculated and the correlations between eigengenes and differentiation stages evaluated (Figure 3(b)). Genes belonging to 9 modules (midnight blue, purple, red, tan, brown, turquoise, black, blue, and green) were enriched with GO and KEGG pathways associated with the early osteoblast stage (days 2-8), such as "cell proliferation" and "osteoblast differentiation processes." Meanwhile, genes belonging to the remaining 7 modules (cyan, salmon, magenta, green yellow, pink, yellow, and gray) were enriched with GO and KEGG pathways associated with the late osteoblast stage (days 10-18), such as "cell adhesion" and "osteoclast differentiation."

Notably, the genes belonging to the cyan and salmon modules exhibited high trait association correlation (r > 0.6), and they were enriched with the following KEGG pathways and GO terms, besides "cell adhesion" and "osteoclast differentiation," as mentioned above: PPAR signaling, subset of extracellular matrix, calcium metabolism, glucose homeostasis, cell signaling, and regulation of bone mineralization. This indicates that these two modules are strongly associated with late osteoblast differentiation (day 18). Thus, from the two modules (cyanide and salmon), we extracted 36 and 62 hub genes, respectively, with high connectivity (MM > 0.9), and these genes were used to find biomarkers for osteoblast differentiation as described below.

3.5. Identification of Markers Associated with Osteoblast Differentiation. We combined the results from gene clustering analysis and gene coexpression networks to determine novel biomarkers for osteoblast differentiation. Focusing on the late osteoblast differentiation stage, we considered genes belonging to clusters 1, 7, and 9 as well as the hub genes belonging to the cyan and salmon modules, as described above. We set the following cut-off values for biomarker candidates: absolute value of fold changes between early differentiation (average expression from samples of days 2-8) and late differentiation (average expression from samples of days 10-18) must be > 2 and at least one of the samples should have an RPKM value > 5 for the gene.

We identified a total of 72 genes as potentially useful marker candidates for osteoblast differentiation (Supplementary Table 3). Functions associated with these genes include "regulation of bone mineralization (GO:0030500)," "extracellular matrix (GO:0031012)," "collagen (GO:0005581)," "collagen binding (GO:0005518)," and "cell adhesion (GO:0007155)" (Supplementary Table 4). Extracellular matrix, including collagen, is important for bone matrix construction [33], and cell adhesion affects cell-matrix interactions and further cell differentiation [34, 35].

Examples of genes associated with each function are as follows. Potential markers associated with "regulation of bone mineralization" are Ifitm5, Bglap, and Bglap2. Potential markers that are functionally associated with extracellular matrix included Abi3bp, Mfap4, Rarres2, Dcn, Dpt, Cilp, and Sod3 as well as the well-known gene Itgbl1 (integrin, beta-like 1; with EGF-like repeat domains) (Figure 4(a)). Dcn (decorin) is also known as a skeletal-related gene [36]. We identified C1qa, Marco, C1qc, Adipoq, and C1qb as collagen-related genes and identified Abi3bp, Srgn, and Dcn as genes related to "collagen binding." Markers associated with "cell adhesion" included genes such as Cd36, Fbln7, Mfap4, Ibsp, and Dpt. In addition, Marco (macrophage receptor with collagen structure), Ntrk2 (neurotrophic tyrosine kinase receptor type 2), and Mfap4 (microfibrillar-associated protein 4) were identified as novel potential markers (Figure 4(b)). Although we identified potentially useful marker candidates for osteoblast differentiation, our study does have limitations. Since we rigorously filtered out 4058 time series genes based on the clustering analysis and WGCNA analysis (trait correlation and hub genes), the current method did not identify the known biomarkers (Fgfr2, Bmp4, and Fgf18) described in the clustering analysis. In addition, because the candidates were mainly identified based on the bioinformatic observations, their biological relevance would need to be further investigated at the cellular and molecular levels experimentally.

3.6. Expression Profiles of IncRNAs during Osteoblast Differentiation. The NONCODE (v4) database was used for lncRNA annotations, and at least one read was mapped onto 37,431 lncRNAs during osteoblast differentiation. Hierarchical clustering was performed based on the expression profiles in each sample. Upon application of similar approaches for known Ensembl genes to lncRNAs, we identified that 54% of the total lncRNAs (20,167 out of 37,431) were differentially expressed across the stages of osteoblast differentiation. Similar to the expression profile of mRNA described above, the expression patterns of lncRNAs were different between the early differentiation stage (2-8 days) and the late differentiation stage (10-18 days) (Figure 5(a)). The number of differentially expressed lncRNAs tended to increase with time (Figure 5(b)).

3.7. Dynamic Changes in lncRNAs during Osteoblast Differentiation. Examination of the temporal expression changes of lncRNAs revealed that levels of 3948 out of 37,431 lncRNAs were dynamically altered over the course of osteoblast differentiation. In PCA analysis using 3948 time series lncRNAs, 97% of the variations in lncRNA expression patterns could be explained by the first two principal components (PCs) (Figure 5(c)). Some of the 3948 time series lncRNAs highlighted were associated with osteogenic function such as "actin cytoskeleton organization," "extracellular matrix organization," "Wnt signaling pathway," "embryonic skeletal system morphogenesis," "Notch signaling pathway," and "collagen fibril organization."

3.8. Clustering Analysis of Time Series lncRNAs. Using R package Mfuzz [20], differentially expressed 3948 time series lncRNAs were divided into 12 clusters, and each cluster contained 142-595 lncRNAs (Figure 6(a)). The lncRNAs contained in 5 clusters (2, 3, 6, 9, and 11) were highly expressed in early osteoblast differentiation stage. Functional annotation analysis showed that these clusters are associated with "regulation of cell proliferation," "actin cytoskeleton organization," and "cell-cycle-related processes." The remaining seven clusters (1, 4, 5, 7, 8, 10, and 12) included lncRNAs that displayed a higher expression in the late osteoblast differentiation stage and related to "cell adhesion," "cell differentiation," "collagen fibril organization," "Wnt signaling pathway," and "phosphorylation processes."

3.9. Identification of lncRNAs Associated with Osteoblast Differentiation. Expression patterns of lncRNAs belonging to cluster 5 were not only very similar but also significantly increased over time, especially between 10 and 18 days. Thus, we focused on cluster 5 to find the lncRNAs markers associated with late differentiation of osteoblasts. Cut-off values were set for lncRNAs (Figure 6(b), upper panel), as in the procedures in novel biomarker identification for known genes.

In total, 89 lncRNAs were identified as potentially useful markers for osteoblast differentiation (Supplementary Table 5). GO enrichment analysis for these lncRNAs disclosed significant enrichment of "cell adhesion (GO:0007155)" and "Wnt signaling pathway (GO:0016055)," which is important in the differentiation and/or function of osteoblasts [31]. For examples, NONMMUG016555, NONMMUG038646, NONMMUG007704, NONMMUG001799, NONMMUG 005875, NONMMUG002249, NONMMUG012067, NONMMUG035551, NONMMUG015487, and NONMMUG 035553 are related to "cell adhesion." NONMMUG038646, NONMMUG016590, and NONMMUG001799 are related to "Wnt signaling pathway (GO:0016055)." In addition, we identified novel candidate lncRNA markers, and two examples are presented in Figure 6(b) (bottom panel).

Next, to see whether the identified lncRNAs could be regulated by common regulation factors, we searched the sequence motifs in the promoter regions (the upstream 1000 bp nucleotide sequences of lncRNAs) of the 89 lncRNAs with background controls generated from all NONCODE lncRNA promoter regions. We performed the motif enrichment analysis using MEME Suite package [23] and the HOCOMOCO v11 mouse transcription factor database [24] to detect known transcription factor binding motifs that are significantly enriched within the promoter region. A total of five significant motifs were identified (Table 3). The most significant enriched motif (WT1_MOUSE.H11-MO.1.A) was the Wilms tumor protein (Wt1) transcription factor binding motif. Recently, it was shown that Wt1 was expressed during development in the limb tissue of E11.5 to E16.5 mice and the loss of the Wt1 expression affected a reduction of non-haematopoetic MSC cells in the E18.5 hindlimb [37]. These findings suggest that Wt1 may play a functional role in the bone development. The second most significant enriched motif (EGR2_MOUSE.H11MO.0.A) was early growth response 2 (Egr2) transcription factor binding motif. Chandra et al. previously suggested that EGFR-induced Egr2 expression is crucial for osteoprogenitor maintenance and new bone formation [38]. Further research on the multiple layers of regulatory mechanisms, including lncRNA expression and transcription regulations, will be required to fully understand the role of these lncRNAs during bone development.

4. Conclusions

We employed public RNA-Seq data (NCBI GEO accession number GSE54461) to profile mRNA and lncRNA expression and identified a series of biomarkers potentially associated with bone differentiation. Data from the current study showed that a combination of bioinformatic approaches, including (1) time series, (2) clustering, and (3) coexpression network analyses, provides an effective means to identify novel candidate markers associated with osteoblast differentiation. Expression patterns of mRNA and lncRNA during osteoblast differentiation were defined as two distinct stages (early and late osteoblast differentiation). Functional annotation showed that the members in the late stage are involved in osteogenesis processes. Notably, common transcription factor binding motifs were enriched in the identified lncRNA markers. Our findings provide a valuable resource for lncRNA study and understanding of the mechanism of mouse osteoblast differentiation. Further research is needed to determine the functions of lncRNA and mRNA identified as potential markers in this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This work was supported by a grant of the Korea Health Technology R&D project through KHIDI funded by the Ministry of Health and Welfare, Republic of Korea (HI14C0175), the Collaborative Genome Program for Fostering New Post-Genome Industry of the National Research Foundation (NRF) funded by the Ministry of Science and ICT (MSIT) (2017M3C9A6047623), and a grant (2016R1D1A1B03930275) of the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Republic of Korea.

Supplementary Materials

Supplementary Table 1: functional annotation of upregulated genes on day 18 (day 2 versus day 18). Supplementary Table 2: functional annotation of downregulated genes on day 18 (day 2 versus day 18). Supplementary Table 3: known Ensembl genes potentially associated with osteoblast differentiation. Supplementary Table 4: enriched GO terms of known Ensembl genes potentially associated with osteoblast differentiation. Supplementary Table 5: lncRNAs (NONCODE v4) potentially associated with osteoblast differentiation. Supplementary Figure 1: Overview of the methods and procedures to identify biomarkers. Supplementary Figure 2: proportion of variance for principal components. Supplementary Figure 3: expression profiles of known mature osteoblast markers. (Supplementary Materials)


[1] A. Rutkovskiy, K. O. Stenslokken, and I. J. Vaage, "Osteoblast differentiation at a glance," Medical Science Monitor Basic Research, vol. 22, pp. 95-106, 2016.

[2] K. Nakashima, X. Zhou, G. Kunkel et al., "The novel zinc finger-containing transcription factor osterix is required for osteoblast differentiation and bone formation," Cell, vol. 108, no. 1, pp. 17-29, 2002.

[3] F. Otto, A. P. Thornell, T. Crompton et al., "Cbfal, a candidate gene for cleidocranial dysplasia syndrome, is essential for osteoblast differentiation and bone development," Cell, vol. 89, no. 5, pp. 765-771, 1997.

[4] T. Komori, H. Yagi, S. Nomura et al., "Targeted disruption of Cbfal results in a complete lack of bone formation owing to maturational arrest of osteoblasts," Cell, vol. 89, no. 5, pp. 755-764, 1997.

[5] X. Y. Zhao and J. D. Lin, "Long noncoding RNAs: a new regulatory code in metabolic control," Trends in Biochemical Sciences, vol. 40, no. 10, pp. 586-596, 2015.

[6] W. Zhang, R. Dong, S. Diao, J. Du, Z. Fan, and F. Wang, "Differential long noncoding RNA/mRNA expression profiling and functional network analysis during osteogenic differentiation of human bone marrow mesenchymal stem cells," Stem Cell Research & Therapy, vol. 8, no. 1, p. 30, 2017.

[7] P. Carninci, T. Kasukawa, S. Katayama et al., "The transcriptional landscape of the mammalian genome," Science, vol. 309, no. 5740, pp. 1559-1563, 2005.

[8] 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.

[9] J. E. Wilusz, H. Sunwoo, and D. L. Spector, "Long noncoding RNAs: functional surprises from the RNA world," Genes & Development, vol. 23, no. 13, pp. 1494-1504, 2009.

[10] T. Zhang, X. Zhang, K. Han et al., "Analysis of long noncoding RNA and mRNA using RNA sequencing during the differentiation of intramuscular preadipocytes in chicken," PLoS One, vol. 12, no. 2, article e0172389, 2017.

[11] M. Kretz, D. E. Webster, R. J. Flockhart et al., "Suppression of progenitor differentiation requires the long noncoding RNA ANCR," Genes & Development, vol. 26, no. 4, pp. 338-343, 2012.

[12] L. Wang, Y. Wang, Z. Li, Z. Li, and B. Yu, "Differential expression of long noncoding ribonucleic acids during osteogenic differentiation of human bone marrow mesenchymal stem cells," International Orthopaedics, vol. 39, no. 5, pp. 1013-1019, 2015.

[13] W. Zhuang, X. Ge, S. Yang et al., "Upregulation of lncRNA MEG3 promotes osteogenic differentiation of mesenchymal stem cells from multiple myeloma patients by targeting BMP4 transcription," Stem Cells, vol. 33, no. 6, pp. 19851997, 2015.

[14] L. Zhu and P. C. Xu, "Downregulated LncRNA-ANCR promotes osteoblast differentiation by targeting EZH2 and regulating Runx2 expression," Biochemical and Biophysical Research Communications, vol. 432, no. 4, pp. 612-617, 2013.

[15] J. P. Kemp, C. Medina-Gomez, K. Estrada et al., "Phenotypic dissection of bone mineral density reveals skeletal site specificity and facilitates the identification of novel loci in the genetic regulation of bone mass attainment," PLoS Genetics, vol. 10, no. 6, article e1004423, 2014.

[16] T. D. Wu and S. Nacu, "Fast and SNP-tolerant detection of complex variants and splicing in short reads," Bioinformatics, vol. 26, no. 7, pp. 873-881, 2010.

[17] A. Mortazavi, B. A. Williams, K. McCue, L. Schaeffer, and B. Wold, "Mapping and quantifying mammalian transcriptomes by RNA-Seq," Nature Methods, vol. 5, no. 7, pp. 621628, 2008.

[18] S. Anders and W. Huber, "Differential expression analysis for sequence count data," Genome Biology, vol. 11, no. 10, article R106, 2010.

[19] M. D. Robinson, D. J. McCarthy, and G. K. Smyth, "edgeR: a bioconductor package for differential expression analysis of digital gene expression data," Bioinformatics, vol. 26, no. 1, pp. 139-140, 2010.

[20] L. Kumar and M. E. Futschik, "Mfuzz: a software package for soft clustering of microarray data," Bioinformation, vol. 2, no. 1, pp. 5-7, 2007.

[21] G. R. Warnes, B. Bolker, L. Bonebakker et al., "gplots: various R programming tools for plotting data. 2013," R Package Version 2.1, 2014.

[22] P. Langfelder and S. Horvath, "WGCNA: an R package for weighted correlation network analysis," BMC Bioinformatics, vol. 9, no. 1, p. 559, 2008.

[23] R. C. McLeay and T. L. Bailey, "Motif enrichment analysis: a unified framework and an evaluation on ChIP data," BMC Bioinformatics, vol. 11, no. 1, p. 165, 2010.

[24] I. V. Kulakovskiy, I. E. Vorontsov, I. S. Yevshin et al., "HOCOMOCO: towards a complete collection of transcription factor binding models for human and mouse via largescale ChIP-Seq analysis," Nucleic Acids Research, vol. 46, no. D1, pp. D252-D259, 2018.

[25] M. B. Eisen, P. T. Spellman, P. O. Brown, and D. Botstein, "Cluster analysis and display of genome-wide expression patterns," Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 25, pp. 14863-14868, 1998.

[26] I. Kalajzic, A. Staal, W. P. Yang et al., "Expression profile of osteoblast lineage at defined stages of differentiation," The Journal of Biological Chemistry, vol. 280, no. 26, pp. 2461826, 2005.

[27] G. S. Stein and B. L. Jane, Molecular Mechanisms Mediating Developmental and Hormone-Regulated Expression of Genes in Osteoblasts in Cellular and Molecular Biology of Bone, Academic Press, San Diego, CA, USA, 1993.

[28] K. Yu, J. Xu, Z. Liu et al., "Conditional inactivation of FGF receptor 2 reveals an essential role for FGF signaling in the regulation of osteoblast function and bone growth," Development, vol. 130, no. 13, pp. 3063-3074, 2003.

[29] P. Ducy and G. Karsenty, "The family of bone morphogenetic proteins," Kidney International, vol. 57, no. 6, pp. 2207-2214, 2000.

[30] T. Shimoaka, T. Ogasawara, A. Yonamine et al., "Regulation of osteoblast, chondrocyte, and osteoclast functions by fibroblast growth factor (FGF)-18 in comparison with FGF-2 and FGF10," The Journal of Biological Chemistry, vol. 277, no. 9, pp. 7493-7500, 2002.

[31] F. Long, "Building strong bones: molecular regulation of the osteoblast lineage," Nature Reviews Molecular Cell Biology, vol. 13, no. 1, pp. 27-38, 2012.

[32] M. Takano, F. Otsuka, Y. Matsumoto et al., "Peroxisome proliferator-activated receptor activity is involved in the osteoblastic differentiation regulated by bone morphogenetic proteins and tumor necrosis factor-a," Molecular and Cellular Endocrinology, vol. 348, no. 1, pp. 224-232, 2012.

[33] A. Matsugaki, G. Aramoto, T. Ninomiya, H. Sawada, S. Hata, and T. Nakano, "Abnormal arrangement of a collagen/apatite extracellular matrix orthogonal to osteoblast alignment is constructed by a nanoscale periodic surface structure," Biomaterials, vol. 37, pp. 134-143, 2015.

[34] Y. S. Lee and C. M. Chuong, "Adhesion molecules in skeletogenesis: I. Transient expression of neural cell adhesion molecules (NCAM) in osteoblasts during endochondral and intramembranous ossification," Journal of Bone and Mineral Research, vol. 7, no. 12, pp. 1435-1446, 1992.

[35] B. Chen, B. Ji, and H. Gao, "Modeling active mechanosensing in cell-matrix interactions," Annual Review of Biophysics, vol. 44, no. 1, pp. 1-32, 2015.

[36] N. C. Ho, L. Jia, C. C. Driscoll, E. M. Gutter, and C. A. Francomano, "A skeletal gene database," Journal of Bone and Mineral Research, vol. 15, no. 11, pp. 2095-2122, 2000.

[37] S. L. McHaffie, Investigating the Role of Wt1 in Bone and Marrow Biology. Doctoral Dissertation, University of Edinburgh, UK, 2014.

[38] A. Chandra, S. Lan, J. Zhu, V. A. Siclari, and L. Qin, "Epidermal growth factor receptor (EGFR) signaling promotes proliferation and survival in osteoprogenitors by increasing early growth response 2 (EGR2) expression," The Journal of Biological Chemistry, vol. 288, no. 28, pp. 20488-20498, 2013.

Minjung Kim, (1) Youngseok Yu, (1) Ji-Hoi Moon (iD), (1,2) InSong Koh (iD), (3,4) and Jae-Hyung Lee (1,2)

(1) Department of Life and Nanopharmaceutical Sciences, Kyung Hee University, Seoul, Republic of Korea

(2) Department of Maxillofacial Biomedical Engineering, School of Dentistry, and Institute of Oral Biology, Kyung Hee University, Seoul, Republic of Korea

(3) Department of Physiology, College of Medicine, Hanyang University, Seoul, Republic of Korea

(4) Department of Biomedical Informatics, Hanyang University, Seoul, Republic of Korea

Correspondence should be addressed to InSong Koh; and Jae-Hyung Lee;

Received 6 September 2017; Revised 22 December 2017; Accepted 17 January 2018; Published 22 March 2018

Academic Editor: Ignazio Piras

Caption: Figure 2: Clustering analysis of time series genes during osteoblast differentiation. (a) Soft clusters of 4058 time series gene expression data using Mfuzz. The numbers on the x-axis (time, 1~9) correspond to the nine time points of osteoblast differentiation (days 2, 4, 6, 8, 10, 12, 14, 16, and 18). (b) Heat maps and representative gene expression profiles (RNA-Seq) of four major clusters. The color key indicates gene expression values.

Caption: Figure 4: Expression patterns of five genes upregulated during late osteoblast differentiation stages. (a) Expression profiles of two known genes, Dcn (left panel) and Itgbl1 (right panel). (b) Expression profiles of three novel biomarker candidates, Marco (left panel), Ntrk2 (middle panel), Mfap4 (right panel).
Table 1: RNA-Seq read mapping summary.

                                           Uniquely     % of uniquely
SRA ID          Sample     Total reads   mapped reads   mapped reads

SRR1146385    B6_2_rep1    20,605,665     18,353,730        89.07
SRR1146386    B6_2_rep2    21,204,511     18,881,462        89.04
SRR1146387    B6_2_rep3    21,196,674     18,867,579        89.01
SRR1146388    B6_4_rep1    21,042,912     18,976,744        90.18
SRR1146389    B6_4_rep2    21,636,813     19,507,642        90.16
SRR1146390    B6_4_rep3    21,628,526     19,496,233        90.14
SRR1146391    B6_6_rep1    22,012,993     19,750,118        89.72
SRR1146392    B6_6_rep2    22,652,935     20,316,973        89.69
SRR1146393    B6_6_rep3    22,656,580     20,317,497        89.68
SRR1146394    B6_8_rep1    16,068,544     14,395,620        89.59
SRR1146395    B6_8_rep2    16,509,407     14,785,805        89.56
SRR1146396    B6_8_rep3    16,498,155     14,769,672        89.52
SRR1146397    B6_10_rep1   11,685,016     10,532,491        90.14
SRR1146398    B6_10_rep2   12,016,157     10,826,863        90.10
SRR1146399    B6_10_rep3   12,016,476     10,827,771        90.11
SRR1146400    B6_12_rep1   32,109,724     29,035,678        90.43
SRR1146401    B6_12_rep2   33,017,225     29,851,897        90.41
SRR1146402    B6_12_rep3   32,970,446     29,805,326        90.40
SRR1146403    B6_14_rep1   22,996,425     20,867,104        90.74
SRR1146404    B6_14_rep2   23,630,602     21,437,833        90.72
SRR1146405    B6_14_rep3   23,616,942     21,422,044        90.71
SRR1146406    B6_16_rep1   26,956,212     24,189,458        89.74
SRR1146407    B6_16_rep2   27,696,210     24,846,901        89.71
SRR1146408    B6_16_rep3   27,659,416     24,807,073        89.69
SRR1146409    B6_18_rep1   22,811,796     20,549,338        90.08
SRR1146410    B6_18_rep2   23,478,940     21,148,934        90.08
SRR1146411    B6_18_rep3   23,452,037     21,114,152        90.03

Table 2: Enriched GO terms and KEGG pathways of cluster 1, 7, and
9 genes.

Cluster        Gene ontology    p value    KEGG pathway     p value

               Wnt-activated    <0.0001    ECM-receptor     0.0007
                 receptor                  interactions

              Collagen type I   <0.0001   Focal adhesion    0.0022

Cluster 1       Calcium ion     0.0001

                Osteoblast      0.0003

                Wnt-protein     0.0003

                 Collagen       <0.0001      Cytokine-       0.004
                  binding                    cytokine

                 Positive       <0.0001
               regulation of

Cluster 7      Extracellular    <0.0001

               Extracellular    <0.0001

                 Positive       0.0001
               regulation of

                 Collagen       <0.0001   PPAR signaling    <0.0001

                Biomineral      <0.0001     Osteoclast       0.001
                  tissue                  differentiation

Cluster 9      Regulation of    <0.0001    Adipocytokine     0.006
                   bone                      signaling
              mineralization                  pathway

               Extracellular    <0.0001      Cytokine-      0.0089
                  region                     cytokine

               Extracellular    <0.0001

Cluster 10     Extracellular    <0.0001

Table 3: Enriched transcriptional binding motifs in the 89 lncRNA
promoter regions.

                          Transcription               Adjusted
Motif ID                     factor        p value    p value

WT1_MOUSE.H11MO.1.A            Wt1        8.014E-06   0.004246
EGR2_MOUSE.H11MO.0.A          Egr2        5.517E-05   0.02887
ZN281_MOUSE.H11MO.0.A         Zn281       5.737E-05   0.03001
ZIC1_MOUSE.H11MO.0.B          Zic1        7.868E-05   0.04092
SP5_MOUSE.H11MO.1.C            Sp5        9.138E-05   0.04737

Figure 1: Global expression patterns of genes at the nine time
points of osteoblast differentiation. (a) Hierarchical clustering
of the transcriptome over the time period of osteoblast
differentiation. All known genes with RPKM values [greater than or
equal to] 3 were used for analysis. (b) Number of genes showing
up-or downregulation during osteoblast differentiation (fold change
[greater than or equal to] 2 or [less than or equal to] 0.5, FDR <
0.05). (c) Volcano plot showing differentially expressed genes in
red and blue. The x- and y-axes represent the magnitude of fold
changes ([log.sub.2] transformed) and adjusted p value
(-[log.sub.2]) by Benjamini-Hochberg correction, respectively. FC =
fold change (d18/d2). (d) Principal component analysis (PCA) of
4058 time series gene expression profiles for different samples
(three replicates at each time point of osteoblast


                Day 4              Day 6              Day 8

Day 2       899 [up arrow]    1755 [up arrow]    2490 [up arrow]
           374 [down arrow]   647 [down arrow]   980 [down arrow]

Day 4                          691 [up arrow]    1492 [up arrow]
                              286 [down arrow]   646 [down arrow]

Day 6                                             815 [up arrow]
                                                 402 [down arrow]

Day 8

Day 10

Day 12

Day 14

Day 16

                Day 10              Day 12              Day 14

Day 2       3169 [up arrow]     4251 [up arrow]     4139 [up arrow]
           1637 [down arrow]   2402 [down arrow]   2571 [down arrow]

Day 4       2294 [up arrow]     3235 [up arrow]     3256 [up arrow]
           1237 [down arrow]   1950 [down arrow]   2109 [down arrow]

Day 6       1681 [up arrow]     2650 [up arrow]     2742 [up arrow]
           1093 [down arrow]   1805 [down arrow]   1959 [down arrow]

Day 8       896 [up arrow]      1650 [up arrow]     1861 [up arrow]
           619 [down arrow]    1315 [down arrow]   1538 [down arrow]

Day 10                          408 [up arrow]      642 [up arrow]
                               299 [down arrow]    563 [down arrow]

Day 12                                              548 [up arrow]
                                                   540 [down arrow]

Day 14

Day 16

                Day 16              Day 18

Day 2       4385 [up arrow]     4431 [up arrow]
           2551 [down arrow]   2807 [down arrow]

Day 4       3448 [up arrow]     3567 [up arrow]
           2187 [down arrow]   2315 [down arrow]

Day 6       2984 [up arrow]     3098 [up arrow]
           2020 [down arrow]   2214 [down arrow]

Day 8       2061 [up arrow]     2362 [up arrow]
           1525 [down arrow]   1706 [down arrow]

Day 10      885 [up arrow]      1213 [up arrow]
           522 [down arrow]    725 [down arrow]

Day 12      847 [up arrow]      1151 [up arrow]
           669 [down arrow]    701 [down arrow]

Day 14      710 [up arrow]      896 [up arrow]
           508 [down arrow]    501 [down arrow]

Day 16                          758 [up arrow]
                               629 [down arrow]

Note: Table made from bar graph.

Figure 3: WGCNA analysis of 4058 time series gene expression data at
the nine time points of osteoblast differentiation. (a) Hierarchical
cluster trees depict coexpression modules identified using WGCNA.
(b) A heat map plot of module-trait association. Each row
corresponds to the module eigengene, the first principal component
of a module. Each column corresponds to trait, the nine time points
of osteoblast differentiation (days 2, 4, 6, 8, 10, 12, 14, 16, and
18). Each cell contains the corresponding correlation value and p
value (left panel). Each of the sixteen modules contains between 16
and 819 genes (right panel). (c) Functional annotations of cyan and
salmon modules. Gene ontology (GO) and KEGG pathway terms and
corresponding p values are shown.


                                      Module-trait relationships

                                        Day 2            Day 3

MEmidnightblue    [Midnight blue]   -0.24 (0.2)       0.31 (0.1)
MEpurple          [Purple]           0.54 (0.003)     0.23 (0.3)
MEred             [Red]              0.61 (7e-04)     0.33 (0.09)
MEtan             [Tan]              0.65 (2e-04)     0.42 (0.03)
MEbrown           [Brown]            0.63 (4e-04)     0.47 (0.01)
MEturquoise       [Turquoise]        0.56 (0.002)     0.45 (0.02)
MEblack           [Black]            0.29 (0.01)      0.26 (0.2)
MEblue            [Blue]             0.46 (0.01)      0.42 (0.03)
MEgreen           [Green]            0.43 (0.02)      0.04 (0.04)
MEcyan            [Cyan]            -0.32 (0.1)      -0.32 (0.1)
MEsalmon          [Salmon]           -0.3 (0.1)      -0.34 (0.08)
MEmagenta         [Magenta]          -0.9 (1e-10)     -0.3 (0.1)
MEgreenyellow     [Greenyellow]     -0.69 (7e-05)    -0.48 (0.01)
MEpink            [Pink]            -0.81 (3e-07)    -0.41 (0.04)
MEyellow          [Yellow]          -0.77 (3e-06)     -0.4 (0.04)
MEgrey            [Grey]             -0.7 (63-05)     -0.3 (0.1)

                              Module-trait relationships

                      Day 6            Day 8            Day 10

MEmidnightblue      0.62 (6e-04)     0.48 (0.01)    -0.011 (1)
MEpurple            0.34 (0.08)      0.18 (0.4)    -0.0073 (1)
MEred               0.23 (0.2)       0.17 (0.4)      0.016 (0.9)
MEtan               0.19 (0.3)       0.11 (0.6)     -0.083 (0.7)
MEbrown              0.3 (0.1)      0.037 (0.9)      -0.24 (0.2)
MEturquoise         0.36 (0.06)      0.12 (0.5)       -0.2 (0.3)
MEblack             0.46 (0.02)      0.45 (0.02)    0.0087 (1)
MEblue              0.39 (0.04)      0.24 (0.2)       0.11 (0.6)
MEgreen              0.4 (0.04)      0.28 (0.2)      0.087 (0.7)
MEcyan             -0.25 (0.2)      -0.25 (0.2)      -0.15 (0.4)
MEsalmon           -0.32 (0.1)      -0.31 (0.1)      -0.17 (0.4)
MEmagenta          -0.11 (0.6)      -0.28 (0.2)       0.19 (0.3)
MEgreenyellow      -0.21 (0.3)      -0.14 (0.5)       0.37 (0.05)
MEpink            -0.093 (0.6)      -0.12 (0.5)       0.24 (0.2)
MEyellow           -0.15 (0.4)     -0.026 (0.9)       0.17 (0.4)
MEgrey             -0.15 (0.5)      -0.27 (0.2)      -0.19 (0.3)

                    Module-trait relationships

                      Day 12           Day 14

MEmidnightblue     -0.16 (0.4)      -0.26 (0.2)
MEpurple           -0.14 (0.5)      -0.38 (0.05)
MEred              -0.17 (0.4)      -0.29 (0.1)
MEtan              -0.22 (0.3)      -0.41 (0.03)
MEbrown            -0.31 (0.1)       -0.3 (0.1)
MEturquoise         -0.3 (0.1)      -0.32 (0.1)
MEblack            -0.27 (0.2)      -0.36 (0.07)
MEblue             -0.26 (0.2)      -0.36 (0.06)
MEgreen            -0.07 (0.2)      -0.33 (0.09)
MEcyan             0.028 (0.9)       0.13 (0.5)
MEsalmon            0.64 (0.8)        0.3 (0.1)
MEmagenta           0.16 (0.4)       0.17 (0.4)
MEgreenyellow       0.33 (0.09)      0.25 (0.2)
MEpink              0.23 (0.2)       0.26 (0.2)
MEyellow            0.23 (0.02)      0.28 (0.3)
MEgrey              0.12 (0.3)       0.25 (0.2)

                    Module-trait relationships

                      Day 16           Day 18

MEmidnightblue     -0.36 (0.06)     -0.37 (0.05)
MEpurple          -0.099 (0.6)      -0.67 (1e-04)
MEred              -0.33 (0.09)     -0.57 (0.02)
MEtan              -0.21 (0.3)      -0.45 (0.02)
MEbrown            -0.03 (0.1)      -0.28 (0.2)
MEturquoise        -0.31 (0.1)      -0.36 (0.06)
MEblack            -0.36 (0.07)     -0.48 (0.01)
MEblue              0.31 (0.1)      -0.47 (0.01)
MEgreen            -0.39 (0.05)     -0.44 (0.02)
MEcyan              0.33 (0.1)        0.8 (42-07)
MEsalmon            0.44 (0.02)      0.64 (3e-04)
MEmagenta           0.15 (0.5)       0.13 (0.5)
MEgreenyellow        0.2 (0.3)      0.094 (0.6)
MEpink              0.23 (0.3)       0.22 (0.3)
MEyellow             0.3 (0.1)      0.032 (0.1)
MEgrey              0.48 (0.1)       0.13 (0.5)

Module            Number of genes

[Midnight blue]   31
[Purple]          112
[Red]             236
[Tan]             90
[Brown]           660
[Turquoise]       819
[Black]           148
[Blue]            716
[Green]           280
[Cyan]            42
[Salmon]          76
[Magenta]         113
[Greenyellow]     110
[Pink]            129
[Yellow]          480
[Grey]            16
Total             4058


Module     Gene ontology    p value       KEGG pathway       p value

Cyan        Complement      < 0.0001     Complement and      < 0.0001
            activation,                    coagulation
             classical                      cascades

            Integral to     < 0.0001       Osteoclast         0.0019
             membrane                    differentiation

             Collagen        0.0002

          Lipid metabolic    0.0006

           Elevation of      0.0008
            calcium ion

Salmon        Glucose       < 0.0001      Adipocytokine      < 0.0001
            homeostasis                     signaling

            Lipoprotein     < 0.0001     PPAR signaling       0.0001
             particle                        pathway

           Adiponectin-     < 0.0001   Glycerophospholipid    0.0012
             mediated                      metabolism

            Biomineral      < 0.0001      Cell adhesion       0.0115
              tissue                    molecules (CAMs)

             Positive       < 0.0001
           regulation of

           Regulation of    < 0.0001

           Extracellular     0.0002

           cGMP-mediated     0.0005

Figure 5: Global expression patterns of lncRNAs at the nine time
points of osteoblast differentiation. (a) Hierarchical clustering
of lncRNAs across the stages of osteoblast differentiation. (b)
Number of lncRNAs showing up- or downregulation during osteoblast
differentiation (fold change [greater than or equal to] 2 or [less
than or equal to] 0.5, FDR < 0.05). (c) Principal component
analysis (PCA) of 3948 time series lncRNA expression profiles for
different samples (three replicates for each time point of
osteoblast differentiation).


                Day 2              Day 4               Day 6

Day 2      1572 [up arrow]    2512 [up arrow]     3686 [up arrow]
           499 [down arrow]   792 [down arrow]   1317 [down arrow]

Day 4                          935 [up arrow]     1904 [up arrow]
                                    5171         742 [down arrow]

Day 6                                             1315 [up arrow]
                                                 628 [down arrow]

Day 8

Day 10

Day 12

Day 14

Day 16

                 Day 8              Day 10              Day 12

Day 2       3448 [up arrow]     4908 [up arrow]     4549 [up arrow]
           1946 [down arrow]   2758 [down arrow]   2943 [down arrow]

Day 4       2270 [up arrow]     3694 [up arrow]     3477 [up arrow]
           1410 [down arrow]   2351 [down arrow]   2754 [down arrow]

Day 6       1765 [up arrow]     2999 [up arrow]     2904 [up arrow]
           1281 [down arrow]   2409 [down arrow]   2666 [down arrow]

Day 8       833 [up arrow]      1741 [up arrow]     1995 [up arrow]
           1022 [down arrow]   1981 [down arrow]   2572 [down arrow]

Day 10                          538 [up arrow]      696 [up arrow]
                               603 [down arrow]    940 [down arrow]

Day 12                                              1014 [up arrow]
                                                   1410 [down arrow]

Day 14

Day 16

                Day 14              Day 16

Day 2       5109 [up arrow]     4683 [up arrow]
           2970 [down arrow]   3124 [down arrow]

Day 4       3967 [up arrow]     3710 [up arrow]
           2743 [down arrow]   2980 [down arrow]

Day 6       3333 [up arrow]     3183 [up arrow]
           2665 [down arrow]   3094 [down arrow]

Day 8       2205 [up arrow]     2369 [up arrow]
           2371 [down arrow]   3032 [down arrow]

Day 10      946 [up arrow]      1063 [up arrow]
           863 [down arrow]    1386 [down arrow]

Day 12      1510 [up arrow]     1349 [up arrow]
           1484 [down arrow]   1888 [down arrow]

Day 14      1575 [up arrow]     1095 [up arrow]
           1116 [down arrow]   1253 [down arrow]

Day 16                          1103 [up arrow]
                               1770 [down arrow]

Note: Table made from bar graph.

Figure 6: Clustering analysis of time series lncRNAs during
osteoblast differentiation. (a) Soft clusters of 3948 time series
lncRNAs generated by Mfuzz. The numbers on the x-axis (time, 1-9)
correspond to the nine time points of osteoblast differentiation
(days 2, 4, 6, 8, 10, 12, 14, 16, and 18). (b) Heat map of cluster
5. The color key indicates gene expression values (top left panel).
Enriched GO terms and corresponding p values in cluster 5 are shown
(top right panel). Gene expression profiles of two representative
lncRNAs, NONMMUG033994 (bottom left panel) and NONMMUG037688 (bottom
right panel), visualized in the bottom panel. * NONCODE v4 database


Gene ontology                  p value

Cell adhesion                  0.0003
Wnt signaling pathway          0.0018
COPYRIGHT 2018 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2018 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Kim, Minjung; Yu, Youngseok; Moon, Ji-Hoi; Koh, InSong; Lee, Jae-Hyung
Publication:International Journal of Genomics
Date:Jan 1, 2018
Previous Article:Genome-Wide Identification, Phylogeny, and Expression Analysis of ARF Genes Involved in Vegetative Organs Development in Switchgrass.
Next Article:Systems Pharmacology-Based Approach of Connecting Disease Genes in Genome-Wide Association Studies with Traditional Chinese Medicine.

Terms of use | Privacy policy | Copyright © 2022 Farlex, Inc. | Feedback | For webmasters |