Novel Insights Into the Role of N6-Methyladenosine RNA Modification in Bone Pathophysiology

RNA modifications in physiology and pathology: Progressing towards application in clinical settings

The potential to modify individual nucleotides through chemical means in order to impact the electrostatic charge, hydrophobic properties, and base pairing of RNA molecules is harnessed in the medical application of stable synthetic RNAs like mRNA vaccines and synthetic small RNA molecules. These modifications are used to either increase or decrease the production of therapeutic proteins. Additionally, naturally occurring biochemical alterations of nucleotides play a role in regulating RNA metabolism and function, thereby modulating essential cellular processes. Research elucidating the mechanisms through which RNA modifications govern fundamental cellular functions in multicellular organisms has enhanced our comprehension of how irregular RNA modification profiles can lead to human diseases. Collectively, these fundamental scientific findings have unveiled the molecular and cellular functions of RNA modifications, offering new opportunities for therapeutic intervention and paving the way for a variety of innovative clinical strategies.

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.

m6A-related lncRNAs as potential biomarkers and the lncRNA ELFN1-AS1/miR-182-5p/BCL-2 regulatory axis in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoid subtype. However, unsatisfactory survival outcomes remain a major challenge, and the underlying mechanisms are poorly understood. N6-methyladenosine (m6A), the most common internal modification of eukaryotic mRNA, participates in cancer pathogenesis. In this study, m6A-associated long non-coding RNAs (lncRNA) were retrieved from publicly available databases. Univariate, LASSO, and multivariate Cox regression analyses were performed to establish an m6A-associated lncRNA model specific to DLBCL. Kaplan-Meier curves, principal component analysis, functional enrichment analyses and nomographs were used to study the risk model. The underlying clinicopathological characteristics and drug sensitivity predictions against the model were identified. Risk modelling based on the three m6A-associated lncRNAs was an independent prognostic factor. By regrouping patients using our model-based method, we could differentiate patients more accurately for their response to immunotherapy. In addition, prospective compounds that can target DLBCL subtypes have been identified. The m6A-associated lncRNA risk-scoring model developed herein holds implications for DLBCL prognosis and clinical response prediction to immunotherapy. In addition, we used bioinformatic tools to identify and verify the ceRNA of the m6A-associated lncRNA ELFN1-AS1/miR-182-5p/BCL-2 regulatory axis. ELFN1-AS1 was highly expressed in DLBCL and DLBCL cell lines. ELFN1-AS1 inhibition significantly reduced the proliferation of DLBCL cells and promoted apoptosis. ABT-263 inhibits proliferation and promotes apoptosis in DLBCL cells. In vitro and in vivo studies have shown that ABT-263 combined with si-ELFN1-AS1 can inhibit DLBCL progression.

Epitranscriptomic modifications in mesenchymal stem cell differentiation: advances, mechanistic insights, and beyond

RNA modifications, known as the "epitranscriptome", represent a key layer of regulation that influences a wide array of biological processes in mesenchymal stem cells (MSCs). These modifications, catalyzed by specific enzymes, often termed "writers", "readers", and "erasers", can dynamically alter the MSCs' transcriptomic landscape, thereby modulating cell differentiation, proliferation, and responses to environmental cues. These enzymes include members of the classes METTL, IGF2BP, WTAP, YTHD, FTO, NAT, and others. Many of these RNA-modifying agents are active during MSC lineage differentiation. This review provides a comprehensive overview of the current understanding of different RNA modifications in MSCs, their roles in regulating stem cell behavior, and their implications in MSC-based therapies. It delves into how RNA modifications impact MSC biology, the functional significance of individual modifications, and the complex interplay among these modifications. We further discuss how these intricate regulatory mechanisms contribute to the functional diversity of MSCs, and how they might be harnessed for therapeutic applications. The review also highlights current challenges and potential future directions in the study of RNA modifications in MSCs, emphasizing the need for innovative tools to precisely map these modifications and decipher their context-specific effects. Collectively, this work paves the way for a deeper understanding of the role of the epitranscriptome in MSC biology, potentially advancing therapeutic strategies in regenerative medicine and MSC-based therapies.

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.

Epigenetic regulatory mechanism of ADAMTS12 expression in osteoarthritis

BACKGROUND

Osteoarthritis (OA) is a degenerative joint disease with lacking effective prevention targets. A disintegrin and metalloproteinase with thrombospondin motifs 12 (ADAMTS12) is a member of the ADAMTS family and is upregulated in OA pathologic tissues with no fully understood molecular mechanisms.

METHODS

The anterior cruciate ligament transection (ACL-T) method was used to establish rat OA models, and interleukin-1 beta (IL-1β) was administered to induce rat chondrocyte inflammation. Cartilage damage was analyzed via hematoxylin-eosin, Periodic Acid-Schiff, safranin O-fast green, Osteoarthritis Research Society International score, and micro-computed tomography assays. Chondrocyte apoptosis was detected by flow cytometry and TdT dUTP nick-end labeling. Signal transducer and activator of transcription 1 (STAT1), ADAMTS12, and methyltransferase-like 3 (METTL3) levels were detected by immunohistochemistry, quantitative polymerase chain reaction (qPCR), western blot, or immunofluorescence assay. The binding ability was confirmed by chromatin immunoprecipitation-qPCR, electromobility shift assay, dual-luciferase reporter, or RNA immunoprecipitation (RIP) assay. The methylation level of STAT1 was analyzed by MeRIP-qPCR assay. STAT1 stability was investigated by actinomycin D assay.

RESULTS

The STAT1 and ADAMTS12 expressions were significantly increased in the human and rat samples of cartilage injury, as well as in IL-1β-treated rat chondrocytes. STAT1 is bound to the promoter region of ADAMTS12 to activate its transcription. METTL3/ Insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2) mediated N6-methyladenosine modification of STAT1 promoted STAT1 mRNA stability, resulting in increased expression. ADAMTS12 expression was reduced and the IL-1β-induced inflammatory chondrocyte injury was attenuated by silencing METTL3. Additionally, knocking down METTL3 in ACL-T-produced OA rats reduced ADAMTS12 expression in their cartilage tissues, thereby alleviating cartilage damage.

CONCLUSION

METTL3/IGF2BP2 axis increases STAT1 stability and expression to promote OA progression by up-regulating ADAMTS12 expression.

Therapeutic approach of natural products that treat osteoporosis by targeting epigenetic modulation

Osteoporosis (OP) is a metabolic disease that affects bone, resulting in a progressive decrease in bone mass, quality, and micro-architectural degeneration. Natural products have become popular for managing OP in recent years due to their minimal adverse side effects and suitability for prolonged use compared to chemically synthesized products. These natural products are known to modulate multiple OP-related gene expressions, making epigenetics an important tool for optimal therapeutic development. In this study, we investigated the role of epigenetics in OP and reviewed existing research on using natural products for OP management. Our analysis identified around twenty natural products involved in epigenetics-based OP modulation, and we discussed potential mechanisms. These findings highlight the clinical significance of natural products and their potential as novel anti-OP therapeutics.

Discovery of potential biomarkers for osteoporosis diagnosis by individual omics and multi-omics technologies

INTRODUCTION

Global aging has made osteoporosis an increasingly serious public health problem. Osteoporotic fractures seriously affect the quality of life of patients and increase disability and mortality rates. Early diagnosis is important for timely intervention. The continuous development of individual- and multi-omics methods is helpful for the exploration and discovery of biomarkers for the diagnosis of osteoporosis.

AREAS COVERED

In this review, we first introduce the epidemiological status of osteoporosis and then describe the pathogenesis of osteoporosis. Furthermore, the latest progress in individual- and multi-omics technologies for exploring biomarkers for osteoporosis diagnosis is summarized. Moreover, we clarify the advantages and disadvantages of the application of osteoporosis biomarkers obtained using the omics method. Finally, we put forward valuable views on the future research direction of diagnostic biomarkers of osteoporosis.

EXPERT OPINION

Omics methods undoubtedly provide greatly contribute to the exploration of diagnostic biomarkers of osteoporosis; however, in the future, the clinical validity and clinical utility of the obtained potential biomarkers should be thoroughly examined. In addition, the improvement and optimization of the detection methods for different types of biomarkers and standardization of the detection process guarantee the reliability and accuracy of the detection results.

YTHDF1 Enhances Chondrogenic Differentiation by Activating the Wnt/β-Catenin Signaling Pathway

Cartilage is derived from the chondrogenic differentiation of stem cells, for which the regulatory mechanism has not been fully elucidated. N6-methyladenosine (m6A) messenger RNA (mRNA) methylation is the most common posttranscriptional modification in eukaryotic mRNAs and is mediated by m6A regulators. However, whether m6A regulators play roles in chondrogenic differentiation is unknown. Herein, we aim to determine the role of a main m6A reader protein, YTH N6-methyladenosine RNA binding protein 1 (YTHDF1), in chondrogenic differentiation regulation. Western blotting (WB) assays found that the expression of YTHDF1 increased during chondrogenic differentiation of human bone marrow mesenchymal stem cells (hBMSCs). The results of quantitative polymerase chain reaction, WB, immunohistochemistry, and Alcian blue staining revealed that overexpression of YTHDF1 increased cartilage matrix synthesis and the expression of chondrogenic markers when hBMSCs, ATDC5 cells, or C3H10T1/2 cells were induced to undergo chondrogenesis. Conversely, chondrogenesis was clearly inhibited when YTHDF1 was knocked down in hBMSCs, ATDC5 cells, or C3H10T1/2 cells. Further RNA sequencing and molecular biology experiments found that YTHDF1 activated the Wnt/β-catenin signaling pathway during chondrogenic differentiation. Finally, the effects of overexpression and knockdown of YTHDF1 on chondrogenic differentiation were reversed by inhibiting or activating β-catenin activity. Therefore, we demonstrated that YTDHF1 promoted chondrogenic differentiation through activation of the Wnt/β-catenin signaling pathway.

METTL3 Regulates Osteoclast Biological Behaviors via iNOS/NO-Mediated Mitochondrial Dysfunction in Inflammatory Conditions

Excessive differentiation of osteoclasts contributes to the disruption of bone homeostasis in inflammatory bone diseases. Methyltransferase-like 3 (METTL3), the core methyltransferase that installs an N6-methyladenosine (m⁶A) modification on RNA, has been reported to participate in bone pathophysiology. However, whether METTL3-mediated m⁶A affects osteoclast differentiation in inflammatory conditions remains unelucidated. In this study, we observed that the total m⁶A content and METTL3 expression decreased during LPS-induced osteoclastogenesis. After knocking down METTL3, we found reduced levels of the number of osteoclasts, osteoclast-related gene expression and bone resorption area. A METTL3 deficiency increased osteoclast apoptosis and pro-apoptotic protein expression. RNA sequencing analysis showed that differentially expressed genes in METTL3-deficient cells were mainly associated with the mitochondrial function. The expression of the mitochondrial function-related genes, ATP production and mitochondrial membrane potential decreased after METTL3 knockdown. Moreover, the most obviously upregulated gene in RNA-Seq was Nos2, which encoded the iNOS protein to induce nitric oxide (NO) synthesis. METTL3 knockdown increased the levels of Nos2 mRNA, iNOS protein and NO content. NOS inhibitor L-NAME rescued the inhibited mitochondrial function and osteoclast formation while suppressing osteoclast apoptosis in METTL3-silenced cells. Mechanistically, a METTL3 deficiency promoted the stability and expression of Nos2 mRNA, and similar results were observed after m⁶A-binding protein YTHDF1 knockdown. Further in vivo evidence revealed that METTL3 knockdown attenuated the inflammatory osteolysis of the murine calvaria and suppressed osteoclast formation. In conclusion, these data suggested that METTL3 knockdown exacerbated iNOS/NO-mediated mitochondrial dysfunction by promoting a Nos2 mRNA stability in a YTHDF1-dependent manner and further inhibited osteoclast differentiation and increased osteoclast apoptosis in inflammatory conditions.

m6A regulator-mediated methylation modification patterns and immune microenvironment infiltration characterization in osteoarthritis

Osteoarthritis (OA) is a common disease in orthopedics. RNA N6-methyladenosine (m6A) exerts an essential effect in a variety of biological processes in the eukaryotes. In this study, we determined the effect of m6A regulators in the OA along with performing the subtype classification. Differential analysis of OA and normal samples in the database of Gene Expression Omnibus identified 9 significantly differentially expressed m6A regulators. These regulators were monitored by a random forest algorithm so as to evaluate the risk of developing OA disease. On the basis of these 9 moderators, a nomogram was established. The results of decision curve analysis suggested that the patients could benefit from a nomogram model. The OA sample was classified as 2 m6A models through a consensus clustering algorithm in accordance with these 9 regulators. These 2 m6A patterns were then assessed with principal component analysis. We also determined the m6A scores for the 2 m6A patterns and their correlation with immune infiltration. The results indicated that type A had a higher m6A score than type B. Thus, we suggest that the m6A pattern may provide a new approach for diagnose and provide novel ideas for molecular targeted therapy of OA.

Exosomes from dysfunctional chondrocytes affect osteoarthritis in Sprague-Dawley rats through FTO-dependent regulation of PIK3R5 mRNA stability

Aims

Exosomes (exo) are involved in the progression of osteoarthritis (OA). This study aimed to investigate the function of dysfunctional chondrocyte-derived exo (DC-exo) on OA in rats and rat macrophages.

Methods

Rat-derived chondrocytes were isolated, and DCs induced with interleukin (IL)-1β were used for exo isolation. Rats with OA (n = 36) or macrophages were treated with DC-exo or phosphate-buffered saline (PBS). Macrophage polarization and autophagy, and degradation and chondrocyte activity of cartilage tissues, were examined. RNA sequencing was used to detect genes differentially expressed in DC-exo, followed by RNA pull-down and ribonucleoprotein immunoprecipitation (RIP). Long non-coding RNA osteoarthritis non-coding transcript (OANCT) and phosphoinositide-3-kinase regulatory subunit 5 (PIK3R5) were depleted in DC-exo-treated macrophages and OA rats, in order to observe macrophage polarization and cartilage degradation. The PI3K/AKT/mammalian target of rapamycin (mTOR) pathway activity in cells and tissues was measured using western blot.

Results

DC-exo inhibited macrophage autophagy (p = 0.002) and promoted M1 macrophage polarization (p = 0.002). DC-exo at 20 μg/ml induced collagen degradation (p < 0.001) and inflammatory cell infiltration (p = 0.023) in rats. OANCT was elevated in DC (p < 0.001) and in cartilage tissues of OA patients (p < 0.001), and positively correlated with patients’ Kellgren-Lawrence grade (p < 0.001). PIK3R5 was increased in DC-exo-treated cartilage tissues (p < 0.001), and OANCT bound to fat mass and obesity-associated protein (FTO) (p < 0.001). FTO bound to PIK3R5 (p < 0.001) to inhibit the stability of PIK3R5 messenger RNA (mRNA) (p < 0.001) and disrupt the PI3K/AKT/mTOR pathway (p < 0.001).

Conclusion

Exosomal OANCT from DC could bind to FTO protein, thereby maintaining the mRNA stability of PIK3R5, further activating the PI3K/AKT/mTOR pathway to exacerbate OA.

Cite this article: Bone Joint Res 2022;11(9):652–668.

Regulatory Role of N6-Methyladenosine (m6A) Modification in Osteoarthritis

Osteoarthritis (OA) is the most common joint disease, usually occurring in middle-aged and elderly people. However, current treatment for OA in its early stages is ineffective, and drug therapy is often ineffective in slowing the progression of the disease. In fact, a deeper understanding of the underlying molecular mechanisms of OA could help us to better develop effective therapeutic measures. N6-methyladenosine (m6A) is a methylation that occurs at the adenosine N6-position, which is the most common internal modification on eukaryotic mRNAs. The role and mechanisms of m6A in mammalian gene regulation have been extensively studied. The “Writer”, “eraser”, and “reader” proteins are key proteins involved in the dynamic regulation of m6A modifications. Recent studies on post-transcriptional regulation alone have shown that m6a modification has an important role in the development of OA. This paper summarizes the specific regulatory processes of M6A in disease and reviews the role of m6A in OA, describing its pathophysiological role and molecular mechanisms, as well as its future research trends and potential clinical applications in OA.

N6-methyladenosine methyltransferase WTAP-stabilized FOXD2-AS1 promotes the osteosarcoma progression through m6A/FOXM1 axis

Long noncoding RNAs (lncRNAs) play critical roles in tumor progression regulation, including osteosarcoma. Evidence indicates that N⁶-methyladenosine (m⁶A) modification modulates mRNA stability to regulate osteosarcoma tumorigenesis. Here, present research aims to detect the roles of m⁶A-modified lncRNA FOXD2-AS1 in the osteosarcoma pathophysiological process. Clinical data unveiled that osteosarcoma patients with higher FOXD2-AS1 expression had a poorer overall survival rate compared to those with lower FOXD2-AS1 expression. Functional research illuminated that FOXD2-AS1 accelerated the migration, proliferation and tumor growth in vitro and in vivo. Mechanistically, a remarkable m⁶A-modified site was found on the 3ʹ-UTR of FOXD2-AS1, and m⁶A methyltransferase WTAP (Wilms’ tumor 1 associated protein) promoted the methylation modification, thus enhancing the stability of FOXD2-AS1 transcripts. Furthermore, FOXD2-AS1 interacted with downstream target FOXM1 mRNA through m⁶A sites, forming a FOXD2-AS1/m⁶A/FOXM1 complex to heighten FOXM1 mRNA stability. In conclusion, these findings propose a novel regulatory mechanism in which m⁶A-modified FOXD2-AS1 accelerates the osteosarcoma progression through m⁶A manner, which may provide new concepts for osteosarcoma tumorigenesis.

Identification and Molecular Analysis of m6A-circRNAs from Cashmere Goat Reveal Their Integrated Regulatory Network and Putative Functions in Secondary Hair Follicle during Anagen Stage

Simple Summary

Cashmere is a natural, high-end textile material. It is derived from the secondary hair follicle (SHFs) tissue in the skin of cashmere goats. Previous studies have indicated that m⁶A modifications in circRNA molecules play important roles in a variety of cells through multiple mechanisms. However, little information is available on the expression profile and functional regulatory characteristics of m⁶A-modified circRNA (m⁶A-circRNA) in SHFs of cashmere goats. In this study, a total of 15 m⁶A-circRNAs were identified. Six of these m⁶A-circRNAs were revealed to have significantly higher expression in skin at anagen than at telogen. To gain insight into the potential regulatory mechanisms of the anagen up-regulated m⁶A-circRNAs, we constructed the regulatory networks along with related pathways in SHFs of cashmere goats. In addition, we found that the expression trends of four m⁶A-circRNAs in the SHFs during SHF cycles were highly similar to their host genes. However, the expression patterns of two m⁶A-circRNAs were inconsistent with the linear RNAs from their host genes in the SHFs of cashmere goats. These results will provide new insights to elucidate the biological functions and regulatory features of m⁶A-circRNA in SHF development and cashmere growth in goats.

Abstract

N⁶-methyladenosine (m⁶A) is the most abundant modification in linear RNA molecules. Over the last few years, interestingly, many circRNA molecules are also found to have extensive m⁶A modification sites with temporal and spatial specific expression patterns. To date, however, little information is available concerning the expression profiling and functional regulatory characteristics of m⁶A modified circRNAs (m⁶A-circRNAs) in secondary hair follicles (SHFs) of cashmere goats. In this study, a total of fifteen m⁶A-circRNAs were identified and characterized in the skin tissue of cashmere goats. Of these, six m⁶A-circRNAs were revealed to have significantly higher expression in skin at anagen compared with those at telogen. The constructed ceRNA network indicated a complicated regulatory relationship of the six anagen up-regulated m⁶A-circRNAs through miRNA mediated pathways. Several signaling pathways implicated in the physiological processes of hair follicles were enriched based on the potential regulatory genes of the six anagen up-regulated m⁶A-circRNAs, such as TGF-beta, axon guidance, ribosome, and stem cell pluripotency regulatory pathways, suggesting the analyzed m⁶A-circRNAs might be essentially involved in SHF development and cashmere growth in cashmere goats. Further, we showed that four m⁶A-circRNAs had highly similar expression trends to their host genes in SHFs of cashmere goats including m⁶A-circRNA-ZNF638, -TULP4, -DNAJB6, and -CAT. However, the expression patterns of two m⁶A-circRNAs (m⁶A-circRNA-STAM2 and -CAAP1) were inconsistent with the linear RNAs from their host genes in the SHFs of cashmere goats. These results provide novel information for eluci-dating the biological function and regulatory characteristics of the m⁶A-circRNAs in SHF development and cashmere growth in goats.

N6-methyladenosine Modification-Related Long Non-Coding RNAs are Potential Biomarkers for Predicting the Prognosis of Patients With Osteosarcoma

Background: The role of N6-methyladenosine (m6A)-related long non-coding RNAs (lncRNAs) in osteosarcoma (OS) has not been fully studied yet. We aimed to identify m6A-related lncRNAs that could act as prognostic biomarkers for OS. Methods: Pearson correlation was performed to identify m6A-related lncRNAs. Univariate and multivariate Cox regression analyses were performed to construct the risk model and assess whether the risk score was an independent prognostic factor for patients with OS. Gene Set Enrichment Analysis (GSEA) was performed to analyze the functions of genes in high-risk and low-risk groups. StarBase and Cytoscape were used to construct a competing endogenous RNA (ceRNA) network based on m6A-related prognostic lncRNA signature. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to analyze the function of genes involved in the ceRNA network. Results: We extracted 122 common lncRNAs from TCGA and Gene Expression Omnibus (GEO) databases. Pearson correlation results revealed 59 significant m6A-related lncRNAs in The Cancer Genome Atlas (TCGA) database, from which 2 were screened to construct a risk signature in TCGA dataset, which was then validated in the GEO dataset. A corresponding risk score was calculated and shown to be an independent prognostic factor for patients with OS. Enrichment analysis indicated that cell proliferation-related biological processes were more common in the high-risk group, while immune-related biological processes were more common in the low-risk group. Moreover, we established a nomogram that had a good ability to predict the overall survival of patients with OS. Additionally, a ceRNA network based on small nucleolar RNA host gene 7 (SNHG7) and small nucleolar RNA host gene 12 (SNHG12) was constructed, with genes that were enriched in hepatocellular carcinoma, gastric cancer, and non-small-cell lung cancer pathways. Conclusion: Our study revealed the prognostic role of m6A-related lncRNAs in OS and identified SNHG7 and SNHG12 as potential biomarkers for predicting the prognosis of patients with OS. These findings have enriched our understanding of the role of m6A modification in the dysregulation of lncRNAs in OS.

Transcriptome-wide m6A methylome during osteogenic differentiation of human adipose-derived stem cells

Objectives

Adipose-derived stem cells are frequently used for bone regeneration both in vitro and in vivo. N⁶-methyladenosine (m⁶A) is the most abundant post-transcriptional modification on eukaryotic RNAs and plays multifaceted roles in development and diseases. However, the regulatory mechanisms of m⁶A in osteogenic differentiation of human adipose-derived stem cells (hASCs) remain elusive. The present study aimed to build the transcriptome-wide m⁶A methylome during the osteogenic differentiation of hASCs.

Materials and methods

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 m⁶A level was colorimetrically measured, and the expression of m⁶A regulators was confirmed by qRT-PCR and western blot. Moreover, m⁶A MeRIP-seq and RNA-seq were performed to build the transcriptome and m⁶A 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.

Results

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 m⁶A-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.

Conclusions

The total m⁶A level was reduced with osteogenic differentiation of hASCs. The transcriptome-wide m⁶A methylome built in the present study indicated quite a few signaling pathways, and hub genes were influenced by m⁶A modification. Future studies based on these epigenetic clues could promote understanding of the mechanisms of osteogenic differentiation of hASCs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-021-02508-1.