Skip to main content
Download PDF

YTHDF2-Mediated m6A methylation inhibition by miR27a as a protective mechanism against hormonal osteonecrosis in BMSCs

BMC Musculoskeletal Disorders volume 25, Article number: 359 (2024) Cite this article

Abstract

Background

With the increasing incidence of steroid-induced necrosis of the femoral head (SNFH), numerous scholars have investigated its pathogenesis. Current evidence suggests that the imbalance between lipogenesis and osteoblast differentiation in bone marrow mesenchymal stem cells (BMSCs) is a key pathological feature of SNFH. MicroRNAs (miRNAs) have strong gene regulatory effects and can influence the direction of cell differentiation. N6-methyladenosine (m6A) is a prevalent epigenetic modification involved in diverse pathophysiological processes. However, knowledge of how miRNAs regulate m6A-related factors that affect BMSC differentiation is limited.

Objective

We aimed to investigate the role of miR27a in regulating the expression of YTHDF2 in BMSCs.

Methods

We compared miR27a, YTHDF2, and total m6A mRNA levels in SNFH-affected and control BMSCs. CCK-8 and TUNEL assays were used to assess BMSC proliferation and apoptosis. Western blotting and qRT‒PCR were used to measure the expression of osteogenic (ALP, RUNX2, and OCN) and lipogenic (PPARγ and C/EBPα) markers. Alizarin Red and Oil Red O staining were used to quantify osteogenic and lipogenic differentiation, respectively. miR27a was knocked down or overexpressed to evaluate its impact on BMSC differentiation and its relationship with YTHDF2. Bioinformatics analyses identified YTHDF2 as a differentially expressed gene in SNFH (ROC analysis) and revealed potential signaling pathways through GSEA. The effects of YTHDF2 silencing on the lipogenic and osteogenic functions of BMSCs were assessed.

Results

miR27a downregulation and YTHDF2 upregulation were observed in the SNFH BMSCs. miR27a knockdown/overexpression modulated YTHDF2 expression, impacting BMSC differentiation. miR27a silencing decreased m6A methylation and promoted osteogenic differentiation, while YTHDF2 silencing exerted similar effects. GSEA suggested potential signaling pathways associated with YTHDF2 in SNFH.

Conclusion

miR27a regulates BMSC differentiation through YTHDF2, affecting m6A methylation and promoting osteogenesis. This finding suggests a potential therapeutic target for SNFH.

Peer Review reports

Introduction

Necrosis of the femoral head is a pathological condition with multiple causes that substantially impacts quality of life [1]. In the USA and Europe, 2–10% of hip replacements are attributed to this condition, while in Asia, this figure is 60% [2]. Among nontraumatic cases, steroid-induced necrosis of the femoral head (SNFH) has the highest incidence and is steadily increasing [3]. According to the latest statistics, without effective intervention, 80% of patients with hormonal femoral head necrosis experience collapse of the femoral head, leading to loss of femoral head function and seriously impacting the patient’s quality of life [4]. SNFH often occurs after heavy, long-term use of hormonal drugs, but its exact pathogenesis remains unclear [5]. To date, numerous studies have been conducted on SNFH from both epigenetic and molecular biology perspectives, demonstrating that SNFH-associated miRNAs are related to the osteoblast differentiation of bone marrow mesenchymal stem cells (BMSCs) involved in lipogenesis, suggesting the potential of miRNAs as biomarkers and therapeutic targets for SNFH [6, 7].

N6-methyladenosine (m6A), a prevalent epigenetic modification, involves the methylation of the nitrogen atom at the 6th position of RNA adenosine [8]. This modification strongly impacts biological processes such as alternative splicing, RNA transport, translation, and degradation. Dysregulation of m6A methylation has been linked to various diseases [9, 10]. The enzymes and proteins associated with RNA m6A methylation can be classified into three groups: m6A transferase, m6A desmethyltransferase, and m6A-binding proteins [11]. The regulatory effects of m6A extend across diverse bodily functions, including germ cell development, embryogenesis, cellular differentiation, metabolism, immune response, and tumorigenesis, as evidenced by numerous studies [12]. Despite growing research on the role of m6A in various diseases, its involvement in SNFH remains underexplored.

Our previous work demonstrated the protective effect of miR27a in SNFH and its potential regulatory role in BMSCs [13]. In bone-related diseases, miR27a was shown to promote osteogenesis by regulating the TGF-β/Smad7 signaling pathway in osteoblasts [14]. miR27a may promote the osteogenic differentiation of bone marrow mesenchymal stem cells in hormonal osteonecrosis of the femoral head by targeting and regulating the PI3K gene [15]. Therefore, miR27a has high potential in osteogenic differentiation. This study aimed to elucidate how miR27a regulates SNFH-derived BMSCs through the m6A-related factor YTHDF2.

Materials and methods

Clinical specimens

Bone marrow and serum samples were collected from patients with SNFH and femoral neck fracture (10 patients each) at the Fifth Affiliated Hospital of Xinjiang Medical University from 2019 to 2023. The inclusion criteria for patients with SNFH were as follows: Steinberg stage IV disease according to X-ray or MRI, a history of glucocorticoid use, no history of hip trauma, no history of renal bone disease, and no history of other abnormal bone metabolic disorders, such as ankylosing spondylitis. The inclusion criteria for patients with femoral neck fracture were patients who needed to undergo one or both hip arthroplasty procedures for Garden stage II, III or IV femoral neck fracture and patients who did not have a history of alcohol consumption, cigarette smoking, glucocorticoid use, or other bone metabolic disorders. Information on the clinical samples is shown in Table 1.

Table 1 Clinical characteristics of patients
Full size table

Extraction of RNA from clinical BMSC samples

BMSCs from clinical samples (experimental and control groups) were harvested and lysed in 200 µL of TRIzol reagent (Invitrogen) to quantify RNA expression. After incubation, precooled chloroform was added, and RNA was subsequently extracted and purified. The expression levels were then determined in 96-well plates via an enzyme labeling system.

Modeling of BMSCs for SNFH

Human BMSC progeny were procured from Procell (Wuhan, China), cultivated, and propagated to the third generation. Subsequently, human BMSCs were subjected to methylprednisolone treatment (Pfizer, USA) to establish a model of SNFH BMSCs.

m6A total RNA quantification

For the quantitative analysis of total m6A RNA in the serum samples, we used the EpiQuik m6A RNA Methylation Kit. Initially, the binding solution was prepared. Subsequently, 80 µL of antibody was added to each well of the 96-well plates, and the plates were incubated in an incubator at 37 °C for 90 min. Next, 50 µL of antibody was added to each well, followed by the addition of a fluorescent solution. The plates were then incubated at room temperature in the dark for 4 min. The absorbance was subsequently measured at 450 nm using an enzyme marker.

Cell cultures and transfection

Initially, the BMSCs were resuspended in DMEM containing 100 IU/ml penicillin and 100 µg/ml streptomycin. The cells were then incubated in a CO2 incubator at 37 °C. YTHDF2, negative control plasmid (NC), siRNA-YTHDF2, and siRNA-NC were obtained from Gemma Pharmaceuticals (Shanghai, China). Transfection of plasmids and siRNAs into BMSCs was carried out using Lipofectamine 2000 (Invitrogen), and the transfection efficiency was assessed using Western blotting (WB) and quantitative real-time polymerase chain reaction (qRT‒PCR).

Quantitative real-time PCR

The primers were designed using Primer Premier 5.0 software, with the U6 gene serving as the internal reference. The primer sequences were sent to Shanghai Biological Company for synthesis, and the synthetic sequences are presented in Table 1. The concentration and purity of the extracted RNA were assessed through spectrophotometry. Reverse transcription of cDNA was carried out following the instructions of the reverse transcription kit, and the reverse transcription process utilized Power SYBR Green (TaKaRa, China). The results were analyzed using a real-time fluorescence quantitative PCR instrument. Thresholds and baselines were adjusted based on the internal reference data to determine the Ct value for each group of samples. The validity of the Ct value was confirmed through melting curve analysis. Finally, the data were analyzed using the 2−ΔΔCT method. The primer sequences are shown in Table 2.

Table 2 PCR Primer sequence et al. in this Study
Full size table

Western blot (WB)

Proteins were extracted from BMSCs and bone marrow tissues using RIPA lysis buffer on ice. The primary antibody was diluted with primary antibody diluent (1:1000), and the desired bands were excised based on marker positions, followed by incubation with the corresponding primary antibody. The secondary antibody (1:2000) was diluted with TBST buffer, and the membrane was incubated with the target bands for 2 h at room temperature on a shaker. β-Actin was used as the internal reference protein. Subsequently, the bands were photographed and subjected to quantitative analysis of gray values using ImageJ software.

Download of microarray data

The SNFH dataset GSE123568 was downloaded from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) [16, 17]. Data cleaning, correction, and normalization were performed using the log2 function in R software (version 4.2.1).

Expression analysis of differentially expressed genes (DEGs)

The differential analysis of gene expression levels between the disease and normal groups was conducted using the R software “limma” package [18]. Differentially expressed genes (DEGs) were screened based on the following criteria: |log2FC| > 1 and p value < 0.05. For visualization of the differential expression patterns of m6A-regulated genes in the GSE123568 dataset, heatmaps were generated using the “ggplot2” [19] and “pheatmap” packages in R software. ROC curves were constructed for the identified YTHDF2-associated DEGs using the “pROC” package in R software (version 1.16.2) [20].

GSEA

The median values of the biomarkers were stratified into groups characterized by high and low expression levels. Subsequently, the samples within the high- and low-expression groups were analyzed for discrepancies utilizing the “limma” package. The logFC values obtained were sorted in descending order, followed by the implementation of gene set enrichment analysis (GSEA). The enrichment analysis of GSEA was conducted utilizing the “clusterProfiler” package in R software [21], with a significance threshold set at P < 0.05. This analysis aimed to identify common functions and pathways associated with many genes within the set of differentially expressed genes.

CCK-8 assay

The BMSCs were grouped and individually prepared as single-cell suspensions in 96-well plates. Except for the blank group, each well received 100 µL of cell suspension containing 2 × 103 cells/well. The plates were then placed in a CO2 incubator, allowing the cells to adhere to the well walls. At the specified time points, the 96-well plate was removed, and 10 µL of CCK-8 solution was added to each well, followed by incubation for 2 h. Finally, the absorbance values were measured at 490 nm using an enzyme-linked immunoassay.

TUNEL assay

The TUNEL reaction solution was applied to cover the cells once the grouped BMSCs had adhered to the slide. The DAPI dye storage solution was then diluted 100-fold with PBS to create a working solution of 10 g/ml. A drop of the DAPI working solution was applied to the cell slice, ensuring coverage of all cells. Drops of antifluorescence quenching sealer were added to the cell slides, and coverslips were used to seal them. Finally, the slides were cleaned with deionized water and observed under a microscope.

Osteogenic induction (OI) and Alizarin Red staining

Different groups of BMSCs were induced using osteogenic induction medium supplemented with 50 mmol/L ascorbic acid, 100 nmol/L dexamethasone, 10 mmol/L β-glycerophosphate, and 15% FBS for 15 days. Following the induction period, the BMSCs were fixed in 4% paraformaldehyde for 20 min and subjected to three washes with PBS. Subsequently, the cells were stained with 4 mg/ml Alizarin Red (pH 8.3) for 30 min, followed by three additional PBS washes. The stained cells were photographed, and optical density (OD) values were measured.

Lipogenic induction and oil red O staining

Distinct subgroups of BMSCs were cultivated in lipid-forming induction medium supplemented with 10 µg/ml insulin, 1 µmol/L dexamethasone, 0.5 mmol/L MIBMX, and 0.1 mmol/L indomethacin for 15 days. Subsequently, the BMSCs were fixed with paraformaldehyde for 1 h and stained for 30 min in the absence of light using an Oil Red O solution. The cell morphology was then observed and imaged under a microscope.

Statistical analysis

The data were analyzed utilizing SPSS 27.0 software. Differences between two groups were assessed using t tests, and one-way ANOVA was used to compare differences among multiple groups. For the bioinformatic analysis experiments, statistical methods provided by the R package in the R language were utilized. The experiment was conducted in triplicate, and a P value < 0.05 was considered to indicate statistical significance.

Results

Downregulation of miR27a and dysregulated m6A modification via YTHDF2 in SNFH BMSCs

To assess the expression levels of miR27a and YTHDF2 in patients with SNFH, we obtained bone marrow specimens from patients with SNFH in a clinical setting. The analysis revealed lower expression of miR27a in the SNFH group than in the control group (Fig. 1A). Subsequently, the expression of total m6A mRNA was examined in both the experimental and control groups, revealing a significant increase in m6A methylation levels in the SNFH group (Fig. 1B). In the selective assay of the m6A regulator YTHDF2, an increase in its expression was observed in the SNFH group (Fig. 1C and D). Furthermore, the expression of osteogenic factors (ALP, RUNX2, and OCN) and lipogenic factors (PPARγ and C/EBPα) in both the experimental and control groups was assessed using WB and PCR analyses. The results demonstrated a decrease in osteogenic capacity and an increase in lipogenic capacity in the SNFH group (Fig. 1E-J). Taken together, these findings suggest that the occurrence of SNFH is linked to the reduced expression of miR27a and the abnormally high expression of m6A involving YTHDF2.

Fig. 1
figure 1

m6A methylation involving miR27a and YTHDF2 is responsible for SNFH. A: miR-27a expression in the SNFH group (n = 10) and control group (n = 10). B: Quantification of total m6A-methylated mRNAs in the SNFH and control groups. C, D: YTHDF2 expression was determined in the SNFH and control groups using WB and qRT‒PCR. E-J: The expression of osteogenesis-related factors (ALP, RUNX2, and OCN) and lipogenesis-related factors (PPARγ and C/EBPα) was determined in the SNFH group and control group using WB and qRT‒PCR. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001

Full size image

miR27a protects BMSCs and negatively regulates YTHDF2

To explore the role of miR27a in regulating BMSCs, we performed targeted knockdown and overexpression experiments using a model system. The effective manipulation of miR27a levels was confirmed by qRT‒PCR (Fig. 2A). Notably, miR27a knockdown/overexpression modulated YTHDF2 expression via a negative regulatory relationship, as revealed by WB and qRT‒PCR (Fig. 2B, C). miR27a has a strong gene regulatory effect, and many studies have demonstrated that miR27a regulates the expression of target genes mainly through direct action on these genes. Therefore, we predict that miR27a directly targets YTHDF2 to negatively regulate its expression. Using CCK-8 and TUNEL assays, we demonstrated that miR27a has a protective effect on BMSCs by promoting cell proliferation and division while inhibiting apoptosis (Fig. 3A, B). Further analysis of lipogenic and osteogenic markers in each group confirmed that miR27a promoted BMSC osteogenic differentiation and suppressed BMSC lipogenic potential (Fig. 3C-H).

Fig. 2
figure 2

miR27a Overexpression/Inhibition and Its Impact on YTHDF2. A: Efficacy of miR27a knockdown and overexpression determined via qRT‒PCR. B, C: The relationship between miR27a and YTHDF2 expression determined via qRT‒PCR and WB. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001. ns, not significant

Full size image
Fig. 3
figure 3

miR27a Overexpression/Inhibition and Its Impact on BMSCs. A: The proliferative effect of miR27a on BMSCs was evaluated using a CCK-8 assay. B: The apoptotic effect of miR27a on BMSCs was assessed using TUNEL assays, scale bar = 50 μm. C-H: The expression of osteogenesis-related factors (ALP, RUNX2, and OCN) and lipogenesis-related factors (PPARγ and C/EBPα) in each group was evaluated via WB and qRT‒PCR. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001. ns, not significant

Full size image

The importance of YTHDF2 in SNFH

To explore the role of m6A modification in SNFH, we conducted bioinformatics analysis on the GEO dataset GSE123567. A heatmap was generated to visualize m6A-associated regulators, which indicated their increased expression in patients with SNFH (Fig. 4A). A comparison of m6A-related regulators between the disease and control groups revealed significant differences, with YTHDF2 exhibiting statistical significance at P < 0.05 (Fig. 4B). To validate the clinical utility of YTHDF2, we conducted a receiver operating characteristic (ROC) curve analysis, yielding an area under the curve (AUC) value of 0.75. Thus, YTHDF2 can be considered a potential biomarker for SNFH (Fig. 4C). Gene set enrichment analysis (GSEA) of YTHDF2 indicated its involvement in enriched pathways such as “CELL_CYCLE”, “OLFACTORY_TRANSDUCTION”, “RIBOSOME”, and “RNA_DEGRADATION”, while showing the inhibition of the “OLFACTORY_TRANSDUCTION” pathway in SNFH (Fig. 4D). Therefore, the upregulation of YTHDF2 plays a crucial role in the pathogenesis of SNFH.

Fig. 4
figure 4

The role of YTHDF2 in SNFH. A: Heatmap representing differentially expressed m6A-related genes in the SNFH gene database GSE123568. The X-axis shows the gene set of the experimental and control samples, and the Y-axis is the set of m6A-related genes, where the color shading represents the magnitude of differential expression. B: Statistical analysis of m6A-related genes using the Wilcoxon signed-rank test (P < 0.05, significant). C: Operating characteristic curve for YTHDF2 in GSE123568. D: GSEA enrichment analysis for YTHDF2 as a single gene in SNFH. The Y-axis shows the ES score, where a positive ES indicates that a functional gene set is enriched in the front of the sorted sequence, and a negative ES indicates that a functional gene set is enriched in the back of the sorted sequence. In general, four pathways were promoted, and one pathway was inhibited

Full size image

Knockdown of YTHDF2 promotes the osteoblast differentiation of BMSCs

To investigate the impact of YTHDF2 silencing, we downregulated its expression (Fig. 5A, B). This intervention was observed to increased the proliferation and inhibit the apoptosis of BMSCs (Fig. 5C, D), indicating that decreased YTHDF2 expression promotes the survival of BMSCs in patients with SNFH. Notably, the silencing of YTHDF2 promoted osteogenic differentiation and diminished lipogenic differentiation in BMSCs (Fig. 5E-J), consistent with our earlier identification of the miR27a-mediated negative regulation of YTHDF2. Furthermore, the knockdown of YTHDF2 promoted the osteogenic differentiation of BMSCs, as further validated through Oil Red O staining and Alizarin Red staining (Fig. 5K). In summary, these findings suggest that the reduced expression of YTHDF2 plays a pivotal role in miR27a-mediated BMSC differentiation.

Fig. 5
figure 5

Increased osteogenesis of BMSCs by silencing YTHDF2 in the SNFH group. A, B: The efficiency of YTHDF2 knockdown was determined by qRT‒PCR and WB. C: The proliferative effects of YTHDF2 silencing on BMSCs were evaluated using a CCK-8 assay. D: The apoptotic effects of YTHDF2 silencing on BMSCs were assessed by TUNEL staining (scale bar = 50 μm). E-J: The expression levels of osteogenesis-related factors (ALP, RUNX2, and OCN) and lipogenesis-related factors (PPARγ and C/EBPα) were determined by Western blotting and qRT‒PCR. K: Alizarin red staining after osteogenic induction. Oil red O staining after lipogenic induction (scale bar = 20 μm). *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001

Full size image

Discussion

With the widespread use of hormonal drugs, the incidence of SNFH has been increasing over the past few years. Although SNFH is a common condition in orthopedics, there is currently no effective treatment for this condition [22]. Therefore, it has become particularly urgent and important to clarify the pathogenesis of SNFH [23]. With increasing research into the mechanism of SNFH, the fate of BMSCs has been found to play an important role in the development of disease [24]. Under orthopedic conditions, the osteogenic differentiation of BMSCs toward lipid differentiation directly impacts the course of the disease. However, the proliferation of BMSCs and their differentiation direction are influenced by various factors. Accordingly, we sought to identify the key regulators of BMSC differentiation.

Current evidence suggests that miR27a plays an important regulatory role in gene expression and is involved in gene modification at the molecular level in various disease cells. This molecule can reportedly regulate the expression of multiple genes and plays a crucial role in the development of diseases [25]. Numerous studies have shown that miR27a plays a critical role in various diseases, including hyperlipidemia [26], breast cancer [27], and rectal cancer [28]. miR27a may utilized as a potential biomarker for these diseases. In the context of orthopedic diseases, studies have demonstrated that miR27a expression is closely associated with rheumatoid arthritis [29]. Furthermore, high miR27a expression can delay bone loss and counteract bone loss caused by radiation and aging [30]. Interestingly, Hsu et al. [31] reported that miR27a can inhibit osteoclasts by affecting the Squstm1/p62 signaling pathway, which may have a preventive effect on osteoporosis. The present study revealed that reduced osteogenic levels in clinical SNFH samples were associated with low miR27a expression. Decreases in BMSC survival and osteogenic capacity were observed after miR27a knockdown. Therefore, our findings are consistent with previous studies demonstrating the protective effect of miR27a on orthopedic diseases, particularly in relation to BMSC survival.

m6A methylation is the most prevalent and abundant RNA modification, and it modifies a wide range of RNAs, including mRNAs, rRNAs, tRNAs, and miRNAs [32]. In degenerative disc disease, osteoporosis, osteosarcoma, and osteoarthritis, noncoding RNAs play an important role in these diseases by interacting with m6A methylation-associated factors to regulate various biological processes [33]. Many studies have shown that m6A can act on target RNAs to trigger relevant pathological changes, where m6A serves as an upstream factor and the relevant RNA acts as a downstream signaling molecule. For example, in cervical cancer, m6A methylation was shown to alter circRNA modifications, thereby increasing the risk of cervical cancer metastasis9. Liu et al. [34] discovered that methylation of noncoding RNAs due to m6A methylation in breast cancer can be utilized to predict relevant therapeutic drug targets and target sites for RNA interactions. Although the regulatory influence of N6-methyladenosine (m6A) on RNA has been extensively documented, microRNAs (miRNAs) offer an additional layer of post-transcriptional control. These noncoding RNA molecules bind to specific target sequences in messenger RNAs (mRNAs), resulting in either decreased translation or mRNA degradation and ultimately regulating gene expression [35].

In the present study, we conducted experiments to demonstrate the role of miR27a in regulating the expression of YTHDF2, a key regulator of m6A methylation, and its impact on the lipogenic/osteogenic differentiation of BMSCs from individuals with SNFH. By specifically knocking down and overexpressing miR27a, we observed that the overexpression of miR27a led to a decrease in YTHDF2 expression, while its knockdown resulted in an increase in YTHDF2 expression. These results indicate that miR27a regulates the expression of YTHDF2, which plays a crucial role in the regulation of lipogenic/osteogenic differentiation in BMSCs from patients with SNFH. Many studies have confirmed that miR27a mainly targets downstream genes through direct targeting; therefore, we predicted that miR27a regulates the expression of YTHDF2 through direct targeting. In conclusion, we predicted that miR27a targeting of YTHDF2-mediated m6A methylation has an important effect on the regulation of lipogenic/osteogenic differentiation in BMSCs from patients with SNFH.

Proteins with the YTF domain act as “readers” in m6A methylation, influencing target mRNA splicing, nuclear translocation, stability, translation, and RNA decay by recognizing and translating m6A sites on various transcripts [36]. The primary function of YTHDF1 is to regulate the translation of modified genes, while YTHDF2 mainly contributes to the degradation of target genes. YTHDF3 collaborates with the former two proteins to either accelerate translation or degradation of the target transcriptome [37]. YTHDF2 stabilizes the transcription of maintenance protein 2 (MCM2) and maintenance protein 5 (MCM5) in microchromosomes, promoting hepatocellular carcinoma progression in hepatitis B virus-infected individuals [38]. In NK cells, YTHDF2 promotes the cellular functions of NK cells and creates the necessary conditions for NK cell survival and proliferation through the formation of a STAT5-YTHDF2 positive feedback loop [39]. Numerous studies have confirmed the involvement of YTHDF2 in the development and progression of various diseases, including glioblastoma [40], liver cancer [41], and malignant glioblastoma [42], highlighting that YTHDF2 plays a crucial role in various diseases. Using bioinformatics research methods, we identified YTHDF2 as a differentially expressed gene in SNFH with clinical diagnostic value. GSEA revealed that YTHDF2 primarily regulates BMSCs in SNFH by promoting cell cycle pathways. Herein, we found that BMSCs exhibited increased proliferative and osteogenic capacities following YTHDF2 silencing, indicating that high YTHDF2 expression promotes the occurrence and development of SNFH.

Currently, the alteration of microcirculation in the femoral head is an important cause of pathological changes and disease progression in patients with SNFH and has been widely recognized by the academic community. Research has revealed that there are many reasons for the deterioration of microcirculation in the femoral head, such as infiltration of immune cells [43] and abnormal deposition of lipids [44]. m6A methylation is strongly associated with the infiltration of immune cells in patients with SNFH and can reduce lipid deposition through the expression of m6A-related genes to improve microcirculation in the femoral head. In the present study, by investigating the effect of low expression of YTHDF2 on reducing the lipogenic differentiation of BMSCs, we expected that a reduction in femoral head lipid deposition would improve the microcirculation. However, further experiments are needed to provide more sufficient evidence for this claim.

At the cellular level, osteoblast aging and disruption of metabolic homeostasis are characteristic pathological changes in SNFH. Many studies have shown that m6A is involved in the regulation of aging in a wide range of cells; for example, the demethylase FTO can delay ovarian aging by reducing the level of methylation in granulosa cells [45]. The regulation of m6A methylation also plays an important role in orthopedic-related diseases, such as osteoporosis and lumbar disc degeneration [46]. Unfortunately, the effect of m6A methylation on osteoblast aging was not further investigated in this study, but we hope to alleviate or treat SNFH by promoting the osteogenic differentiation of MSCs and bone homeostasis in the femoral head.

In summary, we provided compelling evidence that miR27a exerts a protective effect on SNFH by negatively regulating YTHDF2 expression, leading to a reduction in m6A methylation in BMSCs. Currently, the clinical diagnosis of hormonal osteonecrosis of the femoral head is mainly based on imaging to assess the blood supply of the femoral head to predict the risk of hormonal osteonecrosis of the femoral head. However, as early imaging changes in hormonal osteonecrosis of the femoral head are not obvious, the clinical diagnosis often misses the optimal treatment period. Therefore, the expression of miR27a and other related genes can be determined to increase the clinical evidence for the early diagnosis of hormonal osteonecrosis of the femoral head. These findings offer novel insights into the early diagnosis and treatment of SNFH.

Data availability

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

References

  1. BATTISTON B, FULCHIGNONI C, FACCENDA C et al. Microsurgery in femoral head osteonecrosis [J]. Minerva Orthop, 2023, 74.

  2. KONARSKI W, POBOŻY T, KONARSKA K, et al. Osteonecrosis related to Steroid and Alcohol Use—An update on pathogenesis [J]. Healthcare. 2023;11:1846.

    Article  PubMed  PubMed Central  Google Scholar 

  3. FANG L, ZHANG G, WU Y et al. SIRT6 Prevents Glucocorticoid-Induced Osteonecrosis of the Femoral Head in Rats [J]. Oxidative Medicine and Cellular Longevity, 2022, 2022.

  4. ZHUANG J, WANG J, ZHANG B et al. The prophylactic effects of naringin on steroid-induced early-stage osteonecrosis in rats: a preliminary study [J]. Cellular and molecular biology (Noisy-le-Grand, France), 2023, 69: 94–104.

  5. CHEN C-Y, RAO S-S, YUE T, et al. Glucocorticoid-induced loss of beneficial gut bacterial extracellular vesicles is associated with the pathogenesis of osteonecrosis [J]. Sci Adv. 2022;8:eabg8335.

    Article  PubMed  Google Scholar 

  6. NAN K, ZHANG Y, ZHANG X, et al. Exosomes from miRNA-378-modified adipose-derived stem cells prevent glucocorticoid-induced osteonecrosis of the femoral head by enhancing angiogenesis and osteogenesis via targeting miR-378 negatively regulated suppressor of fused (Sufu) [J]. Stem Cell Research & Therapy; 2021. p. 12.

  7. LIU Y, ZONG Y, SHAN H, et al. MicroRNA-23b-3p participates in steroid-induced osteonecrosis of the femoral head by suppressing ZNF667 expression [J]. Steroids. 2020;163:108709.

    Article  CAS  PubMed  Google Scholar 

  8. ZHANG H, SHI X, HUANG T, et al. Dynamic landscape and evolution of m6A methylation in human [J]. Nucleic acids research; 2020. p. 48.

  9. LIANG L, ZHU Y, LI J et al. ALKBH5-mediated m6A modification of circCCDC134 facilitates cervical cancer metastasis by enhancing HIF1A transcription [J]. J Experimental Clin Cancer Res, 2022, 41.

  10. KUMARI R, RANJAN P, SULEIMAN Z et al. mRNA modifications in cardiovascular biology and disease: with a focus on m6A modification [J]. Cardiovascular Res, 2021, 118.

  11. ZHU Z M, HUO F C, ZHANG J et al. Crosstalk between m6A modification and alternative splicing during cancer progression [J]. Clin Translational Med, 2023, 13.

  12. JAYASREE P, DUTTA S, KAREMORE P, et al. Crosstalk between m6A RNA methylation and miRNA Biogenesis in Cancer: an unholy Nexus [J]. Molecular Biotechnology; 2023.

  13. CUI Y, HUANG T, ZHANG Z et al. The potential effect of BMSCs with miR-27a in improving steroid-induced osteonecrosis of the femoral head [J]. Sci Rep, 2022, 12.

  14. BAI Y, LIU Y, JIN S et al. Expression of microRNA–27a in a rat model of osteonecrosis of the femoral head and its association with TGF–β/Smad7 signalling in osteoblasts [J]. Int J Mol Med, 2018, 43.

  15. TANG J, YU H, WANG Y et al. miR-27a promotes osteogenic differentiation in glucocorticoid-treated human bone marrow mesenchymal stem cells by targeting PI3K [J]. J Mol Histol, 2021, 52.

  16. BARRETT T, TROUP D. Nucleic Acids Res. 2011;39:D1005–10. WILHITE S, et al. NCBI GEO: Archive for functional genomics data sets - Update [J].

  17. JIA Y, ZHANG Y, LI S et al. Identification and assessment of novel dynamic biomarkers for monitoring non-traumatic osteonecrosis of the femoral head staging [J]. Clin Translational Med, 2023, 13.

  18. RITCHIE M E, PHIPSON B. WU D, Limma powers differential expression analyses for RNA-sequencing and microarray studies [J]. (1362–4962 (Electronic)).

  19. WICKHAM H. ggplot2 [M]. 2016.

  20. ROBIN X, TURCK N FAU - HAINARD A, HAINARD A FAU - TIBERTI N et al. pROC: an open-source package for R and S + to analyze and compare ROC curves [J]. (1471–2105 (Electronic)).

  21. YU G, WANG L-G, HAN Y, et al. clusterProfiler: an R Package for comparing Biological themes among Gene clusters [J]. OMICS. 2012;16:284–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. LI W, LI W, ZHANG W, et al. Exogenous melatonin ameliorates steroid-induced osteonecrosis of the femoral head by modulating ferroptosis through GDF15-mediated signaling [J]. Stem Cell Research & Therapy; 2023. p. 14.

  23. XU H, ZENG Q, ZOU K et al. Glucocorticoid-induced activation of NOX/ROS/NF-κB signaling in MSCs contributes to the development of GONFH [J]. Apoptosis, 2023, 28.

  24. JIA B, JIANG Y, YAO Y et al. Baicalin attenuates dexamethasone-induced apoptosis of bone marrow mesenchymal stem cells by activating the hedgehog signaling pathway [J]. Chin Med J, 2023, Publish Ahead of Print.

  25. ÇOKER GÜRKAN A, KOYUNCU K, OBAKAN YERLIKAYA P et al. miR27a, a fine-tuning molecule, interacts with growth hormone (GH) signaling and ornithine decarboxylase (ODC) via targeting STAT5 [J]. Amino Acids, 2022, 54.

  26. XIANG Y, MAO L, ZUO M-L et al. The Role of MicroRNAs in Hyperlipidemia: From Pathogenesis to Therapeutical Application [J]. Mediators of Inflammation, 2022, 2022: 1–10.

  27. SANGUANSIN S, SAELEE P, KRITSIRIWUTHINAN K, et al. The Association of pre-miR27a Gene Polymorphism and Clinicopathological Data in Thai breast Cancer patients [J]. Asian Pac J cancer Prevention: APJCP. 2023;24:2055–9.

    Article  CAS  Google Scholar 

  28. WU CS, GUO Z J, ZHAO Y F, et al. Relationship between single nucleotide polymorphisms in the genes of microRNA biosynthesis proteins and the prognosis in colorectal cancer [J]. Chin J Cancer Prev Treat. 2020;27:1070–4.

    Google Scholar 

  29. LIU W, SHENG L, NIE L et al. Functional interaction between long non-coding RNA and microRNA in rheumatoid arthritis [J]. J Clin Lab Anal, 2020, 34.

  30. CHANDRA A, LAGNADO A, FARR J et al. Bone marrow adiposity in models of Radiation- and aging-related bone loss is dependent on Cellular senescence [J]. J Bone Miner Res, 2022, 37.

  31. WANG S, MARUYAMA E, MARTINEZ J et al. MicroRNA-27a is essential for bone remodeling by modulating p62-mediated osteoclast signaling [M]. 2022.

  32. QIN Y, LI L. Role of m6A RNA methylation in cardiovascular disease (review) [J]. Int J Mol Med. 2020;46:1958–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. HAN J, KONG H, WANG, X-Q et al. Novel insights into the interaction between N6-methyladenosine methylation and noncoding RNAs in musculoskeletal disorders [J]. Cell Prolif, 2022, 55.

  34. LIU Y, LENG P, LIU Y, et al. Crosstalk between methylation and ncRNAs in breast Cancer: therapeutic and diagnostic implications [J]. Int J Mol Sci. 2022;23:15759.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. HUSSEN B, HIDAYAT H. SALIHI A, et al. MicroRNA: a signature for cancer progression [J]. Volume 138. Biomedicine & Pharmacotherapy; 2021. p. 111528.

  36. MA S, SUN B. YTHDF2 orchestrates tumor-associated macrophage reprogramming and controls antitumor immunity through CD8 T cells [J]. Nat Immunol. 2023;24:1–12.

    Article  Google Scholar 

  37. SHI R, YING S, LI Y, et al. Linking the YTH domain to cancer: the importance of YTH family proteins in epigenetics [J]. Volume 12. Cell Death & Disease; 2021. p. 346.

  38. YANG Y, YAN Y, YIN J, et al. O-GlcNAcylation of YTHDF2 promotes HBV-related hepatocellular carcinoma progression in an N6-methyladenosine-dependent manner [J]. Signal Transduct Target Therapy. 2023;8:63.

    Article  CAS  Google Scholar 

  39. MA S, YAN J, BARR T et al. The RNA m6A reader YTHDF2 controls NK cell antitumor and antiviral immunity [J]. J Exp Med, 2021, 218.

  40. FANG R, CHEN X, ZHANG S et al. EGFR/SRC/ERK-stabilized YTHDF2 promotes cholesterol dysregulation and invasive growth of glioblastoma [J]. Nat Commun, 2021, 12.

  41. HOU J, ZHANG H. YTHDF2 reduction fuels inflammation and vascular abnormalization in hepatocellular carcinoma [J]. Mol Cancer. 2019;18:163.

    Article  PubMed  PubMed Central  Google Scholar 

  42. DIXIT D, PRAGER B, GIMPLE R et al. The RNA m6A reader YTHDF2 maintains Oncogene expression and is a targetable dependency in Glioblastoma Stem cells [J]. Cancer Discov, 2020, 11.

  43. FANG W, PENG P, LIN K et al. m6A methylation modification and immune infiltration analysis in osteonecrosis of the femoral head [J]. J Orthop Surg Res, 2024, 19.

  44. CHENG C, ZHANG H, ZHENG J et al. METTL14 benefits the mesenchymal stem cells in patients with steroid-associated osteonecrosis of the femoral head by regulating the m6A level of PTPN6 [J]. Aging, 2021, 13.

  45. JIANG Z-X, WANG Y-N, LI Z-Y, et al. The m6A mRNA demethylase FTO in granulosa cells retards FOS-dependent ovarian aging [J]. Cell Death & Disease; 2021. p. 12.

  46. LI G, MA L, HE S et al. WTAP-mediated mA modification of lncRNA NORAD promotes intervertebral disc degeneration [J]. Nat Commun, 2022, 13.

Download references

Funding

National Natural Science Foundation of China, No. 81960396 (to CY).

Author information

Authors and Affiliations

  1. The Fifth Clinical Medical College of Xinjiang Medical University, Xinjiang Uygur Autonomous Region, Urumqi, 830011, China

    Tianyi Yuan, Hongjiang Liu, Maimaitiyibubaji Abudoukadier, Zengqiang Yang & Zhiheng Zhou

  2. Department of Orthopedic Center, The Fifth Affiliated Hospital of Xinjiang Medical University, Xinjiang Uygur Autonomous Region, Urumqi, 830011, China

    Yong Cui

Authors
  1. Tianyi Yuan
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Hongjiang Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Maimaitiyibubaji Abudoukadier
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Zengqiang Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Zhiheng Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yong Cui
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Tianyi Yuan was responsible for the experimental manipulation, experimental design, and paper writing; Hongjiang Liu was responsible for the experimental manipulation and data organization; Maimaitiyibubaji was responsible for the experimental manipulation and data analysis; Zengqiang Yang was responsible for the collection of clinical samples; Zhiheng Zhou was responsible for the statistical analysis. Yong Cui was responsible for experimental design, supervision, and proofreading.All authors reviewed the manuscript.

Corresponding author

Correspondence to Yong Cui.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of the Fifth Affiliated Hospital of Xinjiang Medical University.All procedures were conducted after obtaining the patients’ informed consent and were approved by the hospital’s ethical committee (Approval ID: XYDWFYLS-2019-06), adhering to the principles of medical ethics.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yuan, T., Liu, H., Abudoukadier, M. et al. YTHDF2-Mediated m6A methylation inhibition by miR27a as a protective mechanism against hormonal osteonecrosis in BMSCs. BMC Musculoskelet Disord 25, 359 (2024). https://doi.org/10.1186/s12891-024-07481-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-024-07481-3

Keywords

BMC Musculoskeletal Disorders

ISSN: 1471-2474

Contact us