Abstract:Background: The regulatory mechanisms of RNA methylation during the processes of osteogenic and adipogenic differentiation of bone marrow mesenchymal stem cells (BMSCs) have yet to be fully understood. The objective of our study was to analyze and validate the contribution of RNA methylation regulators to the mechanisms underlying the osteogenic and adipogenic differentiation of rat BMSCs. Methods: We downloaded the GSE186026 from the Gene Expression Omnibus (GEO). Differentially expressed genes (DEGs) were screened using the DESeq2 package in R software (version 3.6.3). A total of 50 RNA methylation genes obtained from literature review and summary were intersected with the previous DEGs to obtain RNA methylation genes, which have different expressions (RM-DEGs). Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were utilized to reveal the functional enrichment. Quantitative real-time polymerase chain reaction (qRT-PCR) was performed to validate RM-DEGs. Protein-protein interaction network (PPI) analysis and visual analysis were performed using STRING and Cytoscape. RM-DEGs regulatory network was constructed to analyze the top 10 hub genes. The relationship between RM-DEGs, some enriched GO and pathways was also been analyzed. The miRNAs and RM-DEGs regulatory networks were established by using miRWalk and TargetScan. Results: As part of our research, we detected varying levels of expression for m6A regulators Mettl3 and Rbm15, as well as m7G regulators Mettl1 and Wdr4, in relation to osteogenic differentiation, along with m6A regulator Fmr1 in adipogenic differentiation. The protein-protein interaction (PPI) networks were constructed for 49 differentially expressed genes (DEGs) related to RNA methylation during the process of osteogenic differentiation, and 13 DEGs for adipogenic differentiation. Moreover, top10 hub genes were calculated. In osteogenic differentiation, Mettl3 regulated the Wnt pathway and Hippo pathway by regulating Smad3, Rbm15 regulated the Notch pathway by Notch1, Mettl1 regulated the PI3K-Akt pathway by Gnb4. In adipogenic differentiation, Fmr1 regulated the PI3K-Akt pathway by Egfr. M6A methylation sites of Smad3, Notch1 and Gnb4 were predicted, and the results showed that all three genes were possibly methylated by m6A, and more than 9 sites per gene were possibly methylated. Finally, we constructed the regulatory networks of Mettl3, Rbm15, Mettl1, and Wdr4 and 109 miRNAs in osteogenic differentiation, Fmr1 and 118 miRNAs in adipogenic differentiation. Conclusions: Mettl3(m6A), Rbm15(m6A), Wdr4 and Mettl1(m7G) were differentially expressed in osteogenic differentiation, while Fmr1(m6A) was differentially expressed in adipogenic differentiation. These findings offered potential candidates for further research on the involvement of RNA methylation in the osteogenic and adipogenic differentiation of BMSCs.

Exploring regulatory elements in low-methylated regions for gene expression prediction

Regulatory elements in low-methylated regions predict directional change of gene expression

RNA-seq assistant: machine learning based methods to identify more transcriptional regulated genes

Advancing regulatory genomics with machine learning

Using methylation data to improve transcription factor binding prediction

Genome-wide in silico prediction of gene expression

Predicting genome-wide DNA methylation using methylation marks, genomic position, and DNA regulatory elements

Gbdmr: identifying differentially methylated CpG regions in the human genome via generalized beta regressions

A Unified Model for Differential Expression Analysis of RNA-seq Data Via L1-Penalized Linear Regression

EPInformer: A Scalable Deep Learning Framework for Gene Expression Prediction by Integrating Promoter-enhancer Sequences with Multimodal Epigenomic Data

Tissue-specific RNA methylation prediction from gene expression data using sparse regression models

A semisupervised model to predict regulatory effects of genetic variants at single nucleotide resolution using massively parallel reporter assays

A generalized linear model for decomposing cis-regulatory, parent-of-origin, and maternal effects on allele-specific gene expression

EMO: Predicting Non-coding Mutation-induced Up- and Down-regulation of Risk Gene Expression using Deep Learning

Penalised regression improves imputation of cell-type specific expression using RNA-seq data from mixed cell populations compared to domain-specific methods

A comprehensive comparison of residue-level methylation levels with the regression-based gene-level methylation estimations by ReGear

Direction-aware functional class scoring enrichment analysis of Infinium DNA methylation data

Predicting regulatory mutations and their target genes by new computational integrative analysis: A study of follicular lymphoma

Computational Screening of Potential Regulators for Mrna-Protein Expression Level Discrepancy

Integrating distal and proximal information to predict gene expression via a densely connected convolutional neural network

Bioinformatics analysis and validation of RNA methylation-related genes in osteogenic and adipogenic differentiation of rat bone marrow mesenchymal stem cells