Emerging role of m6A modification in osteogenesis of stem cells

Unraveling Key m6A Modification Regulators Signatures in Postmenopausal Osteoporosis through Bioinformatics and Experimental Verification

OBJECTIVE

Bone marrow mesenchymal stem cells (BMSCs) show significant potential for osteogenic differentiation. However, the underlying mechanisms of osteogenic capability in osteoporosis-derived BMSCs (OP-BMSCs) remain unclear. This study aims to explore the impact of YTHDF3 (YTH N6-methyladenosine RNA binding protein 3) on the osteogenic traits of OP-BMSCs and identify potential therapeutic targets to boost their bone formation ability.

METHODS

We examined microarray datasets (GSE35956 and GSE35958) from the Gene Expression Omnibus (GEO) to identify potential m6A regulators in osteoporosis (OP). Employing differential, protein interaction, and machine learning analyses, we pinpointed critical hub genes linked to OP. We further probed the relationship between these genes and OP using single-cell analysis, immune infiltration assessment, and Mendelian randomization. Our in vivo and in vitro experiments validated the expression and functionality of the key hub gene.

RESULTS

Differential analysis revealed seven key hub genes related to OP, with YTHDF3 as a central player, supported by protein interaction analysis and machine learning methodologies. Subsequent single-cell, immune infiltration, and Mendelian randomization studies consistently validated YTHDF3's significant link to osteoporosis. YTHDF3 levels are significantly reduced in femoral head tissue from postmenopausal osteoporosis (PMOP) patients and femoral bone tissue from PMOP mice. Additionally, silencing YTHDF3 in OP-BMSCs substantially impedes their proliferation and differentiation.

CONCLUSION

YTHDF3 may be implicated in the pathogenesis of OP by regulating the proliferation and osteogenic differentiation of OP-BMSCs.

METTL3-mediated m6A modification of SOX4 regulates osteoblast proliferation and differentiation via YTHDF3 recognition

N6-methyladenosine (m6A), the most prevalent internal modification in mRNA, is related to the pathogenesis of osteoporosis (OP). Although methyltransferase Like-3 (METTL3), an m6A transferase, has been shown to mitigate OP progression, the mechanisms of METTL3-mediated m6A modification in osteoblast function remain unclear. Here, fluid shear stress (FSS) induced osteoblast proliferation and differentiation, resulting in elevated levels of METTL3 expression and m6A modification. Through Methylated RNA Immunoprecipitation Sequencing (MeRIP-seq) and Transcriptomic RNA Sequencing (RNA-seq), SRY (Sex Determining Region Y)-box 4 (SOX4) was screened as a target of METTL3, whose m6A-modified coding sequence (CDS) regions exhibited binding affinity towards METTL3. Further functional experiments demonstrated that knockdown of METTL3 and SOX4 hampered osteogenesis, and METTL3 knockdown compromised SOX4 mRNA stability. Via RNA immunoprecipitation (RIP) assays, we further confirmed the direct interaction between METTL3 and SOX4. YTH N6-Methyladenosine RNA Binding Protein 3 (YTHDF3) was identified as the m6A reader responsible for modulating SOX4 mRNA and protein levels by affecting its degradation. Furthermore, in vivo experiments demonstrated that bone loss in an ovariectomized (OVX) mouse model was reversed through the overexpression of SOX4 mediated by adeno-associated virus serotype 2 (AAV2). In conclusion, our research demonstrates that METTL3-mediated m6A modification of SOX4 plays a crucial role in regulating osteoblast proliferation and differentiation through its recognition by YTHDF3. Our research confirms METTL3-m6A-SOX4-YTHDF3 as an essential axis and potential mechanism in OP.

Advanced glycation end products impair bone marrow mesenchymal stem cells osteogenesis in periodontitis with diabetes via FTO-mediated N6-methyladenosine modification of sclerostin

BACKGROUND

Diabetes mellitus (DM) and periodontitis are two prevalent diseases with mutual influence. Accumulation of advanced glycation end products (AGEs) in hyperglycemia may impair cell function and worsen periodontal conditions. N6-methyladenosine (m6A) is an important post-transcriptional modification in RNAs that regulates cell fate determinant and progression of diseases. However, whether m6A methylation participates in the process of periodontitis with diabetes is unclear. Thus, we aimed to investigate the effects of AGEs on bone marrow mesenchymal stem cells (BMSCs), elucidate the m6A modification mechanism in diabetes-associated periodontitis.

METHODS

Periodontitis with diabetes were established by high-fat diet/streptozotocin injection and silk ligation. M6A modifications in alveolar bone were demonstrated by RNA immunoprecipitation sequence. BMSCs treated with AGEs, fat mass and obesity associated (FTO) protein knockdown and sclerostin (SOST) interference were evaluated by quantitative polymerase chain reaction, western blot, immunofluorescence, alkaline phosphatase and Alizarin red S staining.

RESULTS

Diabetes damaged alveolar bone regeneration was validated in vivo. In vitro experiments showed AGEs inhibited BMSCs osteogenesis and influenced the FTO expression and m6A level in total RNA. FTO knockdown increased the m6A levels and reversed the AGE-induced inhibition of BMSCs differentiation. Mechanically, FTO regulated m6A modification on SOST transcripts, and AGEs affected the binding of FTO to SOST transcripts. FTO knockdown accelerated the degradation of SOST mRNA in presence of AGEs. Interference with SOST expression in AGE-treated BMSCs partially rescued the osteogenesis by activating Wnt Signaling.

CONCLUSIONS

AGEs impaired BMSCs osteogenesis by regulating SOST in an m6A-dependent manner, presenting a promising method for bone regeneration treatment of periodontitis with diabetes.

WTAP-mediated m6A modification modulates bone marrow mesenchymal stem cells differentiation potential and osteoporosis

An imbalance in the differentiation potential of bone marrow mesenchymal stem cells (BMSCs) is an important pathogenic mechanism underlying osteoporosis (OP). N6-methyladenosine (m⁶A) is the most common post-transcriptional modification in eukaryotic cells. The role of the Wilms' tumor 1-associated protein (WTAP), a member of the m⁶A functional protein family, in regulating BMSCs differentiation remains unknown. We used patient-derived and mouse model-derived samples, qRT-PCR, western blot assays, ALP activity assay, ALP, and Alizarin Red staining to determine the changes in mRNA and protein levels of genes and proteins associated with BMSCs differentiation. Histological analysis and micro-CT were used to evaluate developmental changes in the bone. The results determined that WTAP promoted osteogenic differentiation and inhibited adipogenic differentiation of BMSCs. We used co-immunoprecipitation (co-IP), RNA immunoprecipitation (RIP), methylated RNA immunoprecipitation (MeRIP), RNA pulldown, and dual-luciferase assay to explore the direct mechanism. Mechanistically, the expression of WTAP increased during osteogenic differentiation and significantly promoted pri-miR-181a and pri-miR-181c methylation, which was recognized by YTHDC1, and increased the maturation to miR-181a and miR-181c. MiR-181a and miR-181c inhibited the mRNA expression of SFRP1, promoting the osteogenic differentiation of BMSCs. Our results demonstrated that the WTAP/YTHDC1/miR-181a and miR-181c/SFRP1 axis regulated the differentiation fate of BMSCs, suggesting that it might be a potential therapeutic target for osteoporosis.

Yth m6A RNA-Binding Protein 1 Regulates Osteogenesis of MC3T3-E1 Cells under Hypoxia via Translational Control of Thrombospondin-1

Peri-implantitis is a major factor affecting implant prognosis, and the specific anatomy of the peri-implant area makes it more vulnerable to the local hypoxic environment caused by inflammation. N6-methyladenosine (m⁶A) plays a vital role in a multitude of biological processes, and its main "reader" Yth m⁶A RNA-binding protein 1 (YTHDF1) is suggested to affect osteogenic differentiation. However, the mechanism underlying the effect of YTHDF1 on osteogenic differentiation under hypoxic conditions remains unclear. To address this question, we examined the expression of YTHDF1 under hypoxia and observed that hypoxia suppressed osteogenic differentiation but promoted the expression of YTHDF1. Then we knocked down YTHDF1 and found decreased levels of osteogenic-related markers, alkaline phosphatase (ALP) activity, and alizarin red staining (ARS) under normoxia or hypoxia treatment. Bioinformatics analysis identified Thrombospondin-1 (THBS1) might be a downstream factor of YTHDF1. The results revealed that YTHDF1 enhanced the stability of THBS1 mRNA, and immunofluorescence assays found co-localization with YTHDF1 and THBS1 under hypoxia. Loss of function studies showed knocking down YTHDF1 or THBS1 exacerbated the osteogenic inhibition caused by hypoxia. All data imply that hypoxia suppresses osteogenic differentiation and promotes the expression of YTHDF1, which translationally regulates THBS1 in an m⁶A-dependent manner, potentially counteracting hypoxia-induced osteogenic inhibition through the YTHDF1/THBS1 pathway. The results of this study reveal for the first time the molecular mechanism of the regulation of osteogenic differentiation by YTHDF1 under hypoxia and suggest that YTHDF1, together with its downstream factor THBS1, may be critical targets to counteract osteogenic inhibition under hypoxic conditions, providing promising therapeutic strategy for the hypoxia-induced bone loss in peri-implantitis.

Downregulation of METTL14 improves postmenopausal osteoporosis via IGF2BP1 dependent posttranscriptional silencing of SMAD1

Osteoporosis (OP) tends to occur in postmenopausal women, making them prone to fractures. N6-methyladenosine (m6A) methylation plays a crucial role in OP. Herein, we aimed to explore the effects of METTL14 on osteogenesis and the underlying mechanism. Osteogenic differentiation was assessed through osteoblast markers expression, cell proliferation, ALP activity, and mineralization, which were detected by qRT-PCR, CCK-8, EdU assay, ALP staining assay, and ARS staining assay, respectively. Osteoporosis was evaluated in OVX mice using qRT-PCR, microcomputed tomography, and H&E staining assay. The levels of METTL14 and SMAD1 were measured using qRT-PCR and western blot, and their interaction was assessed using RIP and luciferase reporter assay. M6A methylation was analyzed using the Me-RIP assay. The results indicated that m6A, METTL14, and SMAD1 levels were downregulated in patients with OP and OVX mice, and upregulated in osteogenic BMSCs. Knockdown of METTL14 suppressed osteogenesis of BMSCs and reduced bone mass of OVX mice. Moreover, silencing of METTL14 positively related to SMAD1 and inhibited m6A modification of SMAD1 by suppressing its stability. IGF2BP1 was identified as the methylation reader, and which knockdown reversed the upregulation induced by SMAD1. Overexpression of SMAD1 reversed the suppression of osteogenic differentiation induced by METTL14 knockdown. In conclusion, interference with METTL14 inhibited osteogenic differentiation of BSMCs by m6A modification of SMAD1 in an IGFBP1 manner, suggesting that METTL14 might be a novel approach for improving osteoporosis.

Comprehensive analysis of the m6A-related molecular patterns and diagnostic biomarkers in osteoporosis

BACKGROUND

N6-methyladenosine (m6A) modification is a critical epigenetic modification in eukaryotes and involves several biological processes and occurrences of diseases. However, the roles and regulatory mechanisms of m6A regulators in osteoporosis (OP) remain unclear. Thus, the purpose of this study is to explore the roles and mechanisms of m6A regulators in OP.

METHODS

The mRNA and microRNA (miRNA) expression profiles were respectively obtained from GSE56815, GSE7158, and GSE93883 datasets in Gene Expression Omnibus (GEO). The differential expression of 21 m6A regulators between high-bone mineral density (BMD) and low-BMD women was identified. Then, a consensus clustering of low-BMD women was performed based on differentially expressed (DE)-m6A regulators. The m6A-related differentially expressed genes (DEGs), the differentially expressed miRNAs (DE-miRNAs), and biological functions were investigated. Moreover, a weighted gene co-expression network analysis (WGCNA) was constructed to identify the OP-related hub modules, hub genes, and the functional pathways. Then, an m6A regulator-target-pathway network and the competing endogenous RNA (ceRNA) network in key modules were constructed. A least absolute shrinkage and selection operation (LASSO) Cox regression model and a Support Vector Machine-Recursive Feature Elimination (SVM-RFE) model were constructed to identify the candidate genes for OP prediction. The receiver operator characteristic (ROC) curves were used to validate the performances of predictive models and candidate genes.

RESULTS

A total of 10,520 DEGs, 13 DE-m6A regulators, and 506 DE-miRNAs between high-BMD and low-BMD women were identified. Two m6A-related subclusters with 13 DE-m6A regulators were classified for OP. There were 5,260 m6A-related DEGs identified between two m6A-related subclusters, the PI3K-Akt, MAPK, and immune-related pathways, and bone metabolism was mainly enriched in cluster 2. Cell cycle-related pathways, RNA methylation, and cell death-related pathways were significantly involved in cluster 1. Five modules were identified as key modules based on WGCNA, and an m6A regulator-target gene-pathway network and the ceRNA network were constructed in module brown. Moreover, three m6A regulators (FTO, YTHDF2, and CBLL1) were selected as the candidate genes for OP.

CONCLUSION

M6A regulators play an important role in the occurrences and diagnosis of OP.