METTL3-Mediated lncRNA m(6)A Modification in the Osteogenic Differentiation of Human Adipose-Derived Stem Cells Induced by NEL-Like 1 Protein.

This study aimed to explore the regulatory mechanism of methyltransferase3 (METTL3) -mediated long non-coding RNA (lncRNA) N6-methyladenosine (m6A) modification in the osteogenic differentiation of human adipose-derived stem cells (hASCs) induced by NEL-like 1 protein (NELL-1). Methylated RNA immunoprecipitation sequencing (MeRIP-seq) and high- throughput sequencing for RNA (RNA-seq) were performed on hASCs. Osteogenic ability was detected by alkaline phosphatase (ALP) staining, Alizarin Red S(ARS) staining, ALP quantification and Quantitative real-time polymerase chain reaction analysis (qRT-PCR). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis predicted the osteogenesis-related pathways enriched for the lncRNAs and identified the target lncRNAs. After overexpression and knockdown of METTL3, methylated RNA immunoprecipitation-qPCR (MeRIP-qPCR) and qRT-PCR were used to detect the levels of m6A modification and the expression of the target lncRNA, and the binding of both was confirmed by RNA binding protein immunoprecipitation (RIP) assay. The effects of lncRNA and METTL3 on phosphorylation of the key proteins of the pathway were detected by western blot analysis. In vitro experiments showed that METTL3 can promote osteogenic differentiation and that its expression level is upregulated. KEGG pathway analysis predicted that lncRNAs with differentially upregulated methylated peaks were enriched mostly in the mitogen-activated protein kinase (MAPK) signaling pathway, in which Serine/threonine protein kinase 3 (STK3) was the predicted target gene of the lncRNA RP11-44 N12.5. The m6A modification and expression of RP11-44 N12.5 were both regulated by METTL3. Subsequently, lncRNA RP11-44 N12.5 and METTL3 were found to regulate the phosphorylation levels of three key proteins in the MAPK signaling pathway, ERK, JNK and p38. This study shows, for the first time, that METTL3 can activate the MAPK signaling pathway by regulating the m6A modification and expression of a lncRNA, thereby enhancing the osteogenic differentiation of hASCs.

Song Y ，Pan Y ，Wu M ，Sun W ，Luo L ，Zhao Z ，Liu J ... -  《-》

Transcriptome-wide m(6)A methylome during osteogenic differentiation of human adipose-derived stem cells.

Adipose-derived stem cells are frequently used for bone regeneration both in vitro and in vivo. N6-methyladenosine (m6A) is the most abundant post-transcriptional modification on eukaryotic RNAs and plays multifaceted roles in development and diseases. However, the regulatory mechanisms of m6A in osteogenic differentiation of human adipose-derived stem cells (hASCs) remain elusive. The present study aimed to build the transcriptome-wide m6A methylome during the osteogenic differentiation of hASCs. hASCs were harvested after being cultured in a basic or osteogenic medium for 7 days, and the osteogenic differentiation was validated by alkaline phosphatase (ALP) and Alizarin Red S staining, ALP activity assay, and qRT-PCR analysis of ALP, RUNX2, BGLAP, SPP1, SP7, and COL1A1 genes. The m6A level was colorimetrically measured, and the expression of m6A regulators was confirmed by qRT-PCR and western blot. Moreover, m6A MeRIP-seq and RNA-seq were performed to build the transcriptome and m6A methylome. Furthermore, bioinformatic analyses including volcano plots, Venn plots, clustering analysis, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, gene sets enrichment analysis, and protein-protein interaction analysis were conducted. In total, 1145 differentially methylated peaks, 2261 differentially expressed genes, and 671 differentially methylated and expressed genes (DMEGs) were identified. GO and KEGG pathway analyses conducted for these DMEGs revealed extensive and osteogenic biological functions. The "PI3K-Akt signaling pathway"; "MAPK signaling pathway"; "parathyroid hormone synthesis, secretion, and action"; and "p53 signaling pathway" were significantly enriched, and the DMEGs in these pathways were identified as m6A-specific key genes. A protein-protein interaction network based on DMEGs was built, and VEGFA, CD44, MMP2, HGF, and SPARC were speculated as the hub DMEGs. The total m6A level was reduced with osteogenic differentiation of hASCs. The transcriptome-wide m6A methylome built in the present study indicated quite a few signaling pathways, and hub genes were influenced by m6A modification. Future studies based on these epigenetic clues could promote understanding of the mechanisms of osteogenic differentiation of hASCs.

Sun W ，Song Y ，Xia K ，Yu L ，Huang X ，Zhao Z ，Liu J ... -  《Stem Cell Research & Therapy》

microRNA expression profiles and the potential competing endogenous RNA networks in NELL-1-induced human adipose-derived stem cell osteogenic differentiation.

Studies have indicated that Nel-like molecule-1 (NELL-1) was an osteoblast-specific cytokine and some specific microRNAs (miRNAs) could serve as competing endogenous RNA (ceRNA) to partake in osteogenic differentiation of human adipose-derived stem cells (hASCs). The aim of this study was to explore the potential functional mechanisms of recombinant human NELL-1 protein (rhNELL-1) during hASCs osteogenic differentiation. rhNELL-1 was added to osteogenic medium to activate osteogenic differentiation of hASCs. High-throughput RNA sequencing (RNA-Seq) was performed and validated by real-time quantitative polymerase chain reaction. Gene ontology functional annotation and Kyoto Encyclopedia of Genes and Genomes pathway analysis were performed to detect the functions of differentially expressed miRNAs and genes. Coding-noncoding gene co-expression network and ceRNA networks were constructed to predict the potential regulatory role of miRNAs. A total of 1010 differentially expressed miRNAs and 1762 differentially expressed messenger RNAs (mRNAs) were detected. miRNA-370-3p, bone morphogenetic protein 2 (BMP2), and parathyroid hormone like hormone (PTHLH) were differentially expressed during NELL-1-induced osteogenesis. Bioinformatic analyses demonstrated that these differentially expressed miRNAs and mRNAs enriched in Rap1 signaling pathway, PI3K-Akt signaling pathway, p53 signaling pathway, Glucagon signaling pathway, and hypoxia-inducible factor-1 signaling pathway, which were important pathways related to osteogenic differentiation. In addition, miRNA-370-3p and has-miR-485-5p were predicted to interact with circ0001543, circ0002405, and ENST00000570267 in ceRNA networks. Based on the gain or loss of functional experiments by transfection, the results showed that miR-370-3p was a key regulator in osteogenic differentiation by targeting BMP2 and disturbing the expression of PTHLH, and participated in NELL-1-stimulated osteogenesis. The present study provided the primary data and evidence for further exploration on the roles of miRNAs and ceRNAs during NELL-1-induced ossification of hASCs.

Yu L ，Cen X ，Xia K ，Huang X ，Sun W ，Zhao Z ，Liu J ... -  《-》

Comprehensive analysis of lncRNA-miRNA-mRNA networks during osteogenic differentiation of bone marrow mesenchymal stem cells.

Long non-coding RNA (lncRNA) plays crucial role in osteogenic differentiation of bone marrow mesenchymal stem cells (BMMSCs), involving in regulation of competing endogenous RNA (ceRNA) mechanisms and conduction of signaling pathways. However, its mechanisms are poorly understood. This study aimed to investigate lncRNAs, miRNAs and mRNAs expression profiles in rat BMMSCs (rBMMSCs) osteogenic differentiation, screen the potential key lncRNA-miRNA-mRNA networks, explore the putative functions and identify the key molecules, as the basis of studying potential mechanism of rBMMSCs osteogenic differentiation driven by lncRNA, providing molecular targets for the management of bone defect. High-throughput RNA sequencing (RNA-seq) was used to determine lncRNAs, miRNAs, and mRNAs expression profiles at 14-day rBMMSCs osteogenesis. The pivotal lncRNA-miRNA and miRNA-mRNA networks were predicted from sequencing data and bioinformatic analysis, and the results were exported by Cytoscape 3.9.0 software. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used for functional exploration. Real-time quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was performed to validate lncRNAs, miRNAs and mRNAs. rBMMSCs were identified, and the osteogenic and adipogenic differentiation ability were detected. A total of 8634 lncRNAs were detected by RNA-seq, and 1524 differential expressed lncRNAs, of which 812 up-regulated and 712 down-regulated in osteo-inductive groups compared with control groups. 30 up-regulated and 61 down-regulated miRNAs, 91 miRNAs were differentially expressed in total. 2453 differentially expressed mRNAs including 1272 up-expressed and 1181 down-expressed were detected. 10 up-regulated lncRNAs were chosen to predict 21 down-regulated miRNAs and 650 up-regulated mRNAs. 49 lncRNA-miRNA and 1515 miRNA-mRNA interactive networks were constructed. GO analysis showed the most important enrichment in cell component and molecular function were "cytoplasm" and "protein binding", respectively. Biological process related to osteogenic differentiation such as "cell proliferation", "wound healing", "cell migration", "osteoblast differentiation", "extracellular matrix organization" and "response to hypoxia" were enriched. KEGG analysis showed differentially expressed genes were mainly enriched in "PI3K-Akt signaling pathway", "Signaling pathway regulating pluripotency of stem cells", "cGMP-PKG signaling pathway", "Axon guidance" and "Calcium signaling pathway". qRT-PCR verified that lncRNA Tug1, lncRNA AABR07011996.1, rno-miR-93-5p, rno-miR-322-5p, Sgk1 and Fzd4 were consistent with the sequencing results, and 4 lncRNA-miRNA-mRNA networks based on validations were constructed, and enrichment pathways were closely related to "PI3K-Akt signaling pathway", "Signaling pathway regulating pluripotency of stem cells" and "Wnt signaling pathway". lncRNAs, miRNAs and mRNAs expression profiles provide clues for future studies on their roles for BMMSCs osteogenic differentiation. Furthermore, lncRNA-miRNA-mRNA networks give more information on potential new mechanisms and targets for management on bone defect.

Liu J ，Yao Y ，Huang J ，Sun H ，Pu Y ，Tian M ，Zheng M ，He H ，Li Z ... -  《BMC GENOMICS》

N6-methyladenosine promotes osteogenic differentiation of PDLSCs from periodontitis patients.

This study aimed to investigate the mechanism of N6-methyladenosine (m6A) in the osteogenic differentiation of periodontal ligament stem cells (PDLSCs) from periodontitis patients. Differentially m6A-methylated lncRNA/mRNA profiles were detected by a m6A epitranscriptomic microarray. Bioinformatics analysis was performed by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis. The transfection efficiency of the lentivirus was detected. The osteogenic activity of PDLSCs from periodontitis patients (PPDLSCs) was assessed. The microarray results showed that 275 lncRNAs and 1292 mRNAs were significantly differentially methylated between PPDLSCs and PDLSCs from healthy people. Among those lncRNAs, lncRNA4114 (transcript_ID: ENST00000444114) showed both reduced m6A methylation levels and expression levels in PPDLSCs. Further bioinformatics analysis predicted that the differentially methylated mRNAs were mainly involved in cell interaction, stem cell pluripotency, and osteogenic differentiation signals. Then, overexpression of methyltransferase like 3 (METTL3) promoted the osteogenic differentiation of PPDLSCs, while knocking down METTL3 showed an inhibitory effect. Furthermore, METTL3 overexpression promotes the stability of lncRNA4114 to upregulate the expression level. Moreover, lncRNA4114 overexpression promoted the osteogenic differentiation of PPDLSCs. METTL3 promotes the osteogenic differentiation of PPDLSCs by regulating the stability of lncRNA4114.

Zhang X ，Liu J ，Gao J ，Sun W ，Chen X ，Wang X ，Qin W ，Jin Z ... -  《-》