Identification and integrated analysis of differentially expressed lncRNAs and circRNAs reveal the potential ceRNA networks during PDLSC osteogenic differentiation

Background: Researchers have been exploring the molecular mechanisms underlying the control of periodontal ligament stem cell (PDLSC) osteogenic differentiation. Recently, long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) were shown to function as competitive endogenous RNAs (ceRNAs) to regulate the effect of microRNAs (miRNAs) on their target genes during cell differentiation. However, comprehensive identification and integrated analysis of lncRNAs and circRNAs acting as ceRNAs during PDLSC osteogenic differentiation have not been performed. Results: PDLSCs were derived from healthy human periodontal ligament and cultured separately with osteogenic induction and normal media for 7 days. Cultured PDLSCs were positive for STRO-1 and CD146 and negative for CD31 and CD45. Osteo-induced PDLSCs showed increased ALP (alkaline phosphatase) activity and up-regulated expression levels of the osteogenesis-related markers ALP, Runt-related transcription factor 2 and osteocalcin. Then, a total of 960 lncRNAs and 1456 circRNAs were found to be differentially expressed by RNA sequencing. The expression profiles of eight lncRNAs and eight circRNAs were measured with quantitative real-time polymerase chain reaction and were shown to agree with the RNA-seq results. Furthermore, the potential functions of lncRNAs and circRNAs as ceRNAs were predicted based on miRanda and were investigated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. In total, 147 lncRNAs and 1382 circRNAs were predicted to combine with 148 common miRNAs and compete for miRNA binding sites with 744 messenger RNAs. These mRNAs were predicted to significantly participate in osteoblast differentiation, the MAPK pathway, the Wnt pathway and the signaling pathways regulating pluripotency of stem cells. Among them, lncRNAs coded as TCONS_00212979 and TCONS_00212984, as well as circRNA BANP and circRNA ITCH, might interact with miRNA34a and miRNA146a to regulate PDLSC osteogenic differentiation via the MAPK pathway. Conclusions: This study comprehensively identified lncRNAs/circRNAs and first integrated their potential ceRNA function during PDLSC osteogenic differentiation. These findings suggest that specific lncRNAs and circRNAs might function as ceRNAs to promote PDLSC osteogenic differentiation and periodontal regeneration.