Genome-Wide Identification of N6-Methyladenosine (m6A) SNPs Associated With Rheumatoid Arthritis

Genome-wide identification of cell type-specific susceptibility genes for Juvenile dermatomyositis through the analysis of N6-methyladenosine-associated SNPs

Genome-wide association studies (GWASs) have pinpointed genetic loci associated with juvenile dermatomyositis (JDM). Functional genes within the GWAS loci may be cell type-specific, but their identity remains largely unknown. N6-methyladenosine (m6A) plays a pivotal role in regulating various cellular processes and is linked to autoimmune diseases. This study aimed to underscore the potential functional genes within the GWAS loci through the analysis of m6A-associated SNPs (m6A-SNPs), specifically within relevant cell types. JDM-associated m6A-SNPs were identified from the GWAS summary dataset. The correlation between m6A-SNPs and gene expression was assessed through bulk tissue and single-cell eQTL analyses. To further investigate the relationship between gene expression and JDM, Mendelian randomization analysis was employed. Additionally, differential expression analyses were conducted on bulk tissues, as well as single-cell transcriptomic data comprising 6 JDM patients and 11 juvenile controls (99,396 cells). Seven m6A-SNPs associated with JDM were identified. Bulk tissue analysis revealed differential expression of HLA-DPA1, HLA-DPB1, MICB, HLA-A, HLA-F, HLA-DQB2, HLA-DRB5, TAP2, PSMB9, MICA, AIF1, and DDX39B influenced by m6A-SNPs, all showing associations with JDM in both differential expression and Mendelian randomization analyses. In single-cell analysis, the six m6A-SNPs within the HLA locus acted as cell-type-specific eQTLs, correlating with the expression of HLA-A, HLA-B, HLA-C, HLA-DPB1, HLA-DQA1, HLA-DQB1 and HLA-DRB1 in myeloid, T or B cells. Notably, these genes displayed abnormal expression in T, B, and myeloid cells of JDM patients. The present study identified m6A-SNPs within JDM susceptibility genes, shedding light on the intricate interplay between m6A-SNPs, gene expression, and JDM.

Increased m6A RNA methylation and METTL3 expression may contribute to the synovitis progression of rheumatoid arthritis

OBJECTIVE

Rheumatoid arthritis (RA) is an autoimmune disease characterized by synovial hyperplasia and progressive bone destruction. The tumor-like growth of fibroblast-like synoviocytes (FLSs) plays a crucial role in the pathogenesis of RA. The N6 methyladenine (m6A) mRNA methylation modification, regulated by methyltransferases (METTL3) and demethylation enzymes, is a novel epigenetic regulator in the development of RA. However, there is limited research on m6A methylation modifications in RA synovitis and a lack of mechanistic studies on their impact on the function of RA-FLSs.

METHODS

This study utilized clinical synovial tissue specimens and FLSs as research subjects. The m6A methylation level and the expression of methyltransferases and demethylation enzymes were detected. RNA interference and gene overexpression methods were employed to investigate the mechanism of METTL3 in RA-FLSs. The study also examined the proliferation, apoptosis, migration, invasion, and cytokine levels of RA-FLSs, as well as the expression of METTL3 in RA animal models.

RESULTS

In this study, we found that m6A methylation levels were elevated in synovial tissues and FLSs of RA patients. Immunohistochemical staining showed that METTL3 and METTL14 levels were up-regulated in synovial tissues of RA, the mRNA levels of METTL3, METTL14, WTAP, FTO, and ALKBH5 were significantly higher in synovial tissues and FLSs of RA patients. Overexpression of METTL3 could promote the proliferation, migration, and secretion of IL-6, RANKL of RA-FLSs; inhibition of METTL3 expression could inhibit the abnormal proliferation, migration, invasion, and secretion of IL-6, RANKL, at the same time promoted the apoptosis and secretion of OPG, thus inhibited RA-FLSs tumor-like growth. In CIA mice, the use of MTX and STM2457 reduced METTL3 expression, synovial hyperplasia and bone destruction.

CONCLUSION

Abnormal modification of m6A methylation exists in synovial tissues and FLSs of RA patients, and inhibition of METTL3 can reduce synovitis and bone destruction. Our findings suggest that m6A methylation might control FLS-mediated tumor-like phenotype, and be a novel target for RA treatment.

RNA N6-methyladenosine modification in arthritis: New insights into pathogenesis

The commonest type of eukaryotic RNA modification, N6-methyladenosine (m6A), has drawn increased scrutiny in the context of pathological functioning as well as relevance in determination of RNA stability, splicing, transportation, localization, and translation efficiency. The m6A modification plays an important role in several types of arthritis, especially osteoarthritis and rheumatoid arthritis. Recent studies have reported that m6A modification regulates arthritis pathology in cells, such as chondrocytes and synoviocytes via immune responses and inflammatory responses through functional proteins classified as writers, erasers, and readers. The aim of this review was to highlight recent advances relevant to m6A modification in the context of arthritis pathogenesis and detail underlying molecular mechanisms, regulatory functions, clinical applications, and future perspectives of m6A in arthritis with the aim of providing a foundation for future research directions.

Exploring the impact of m6A modification on immune diseases: mechanisms and therapeutic implication

N6-methyladenosine (m6A) is a chemical modification of RNA and has become a widely discussed topic among scientific researchers in recent years. It is distributed in various organisms, including eukaryotes and bacteria. It has been found that m6A is composed of writers, erasers and readers and is involved in biological functions such as splicing, transport and translation of RNA. The balance of the human immune microenvironment is important for human health abnormalities. Increasing studies have found that m6A affects the development of immune diseases such as inflammatory enteritis and systemic lupus erythematosus (SLE) by participating in the homeostatic regulation of the immune microenvironment in vivo. In this manuscript, we introduce the composition, biological function, regulation of m6A in the immune microenvironment and its progression in various immune diseases, providing new targets and directions for the treatment of immune diseases in clinical practice.

IGF2BP3 regulates the expression of RRM2 and promotes the progression of rheumatoid arthritis via RRM2/Akt/MMP-9 pathway

BACKGROUND

Rheumatoid arthritis (RA) is a common inflammatory and autoimmune disease. Ribonucleotide Reductase Regulatory Subunit M2 (RRM2) is a crucial and a rate-limiting enzyme responsible for deoxynucleotide triphosphate(dNTP) production. We have found a high expression level of RRM2 in patients with RA, but the molecular mechanism of its action remains unclear.

METHODS

We analyzed the expression of hub genes in RA using GSE77298 datasets downloaded from Gene Expression Omnibus database. RRM2 and insulin-like growth factor-2 messenger ribonucleic acid (mRNA)-binding protein 3 (IGF2BP3) gene knockdown was achieved by infection with lentiviruses. The expression of RRM2, IGF2BP3, matrix metalloproteinase (MMP)-1, and MMP-9 were detected via western blotting assay. Cell viability was detected via 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. MeRIP-qRT-PCR was performed to test the interaction of IGF2BP3 and RRM2 mRNA via m6A modification. Cell proliferation was determined by clone formation assay. Migration and invasion assays were performed using transwell Boyden chamber.

RESULTS

RRM2 and IGF2BP3 were highly expressed in clinical specimens and tumor necrosis factor alpha (TNF-α) and interleukin (IL)-1β-stimulated synovial cells. RRM2 and IGF2BP3 knockdown inhibited the proliferation, migration, and invasion of MH7A cells. The inhibitory effects of IGF2BP3 knockdown were effectively reversed by simultaneously overexpressing RRM2 in MH7A cells. By analyzing N6-methyladenosine (m6A)2Target database, five m6A regulatory target binding sites for IGF2BP3 were identified in RRM2 mRNA, suggesting a direct relationship between IGF2BP3 and RRM2 mRNA. Additionally, in RRM2 small hairpin (sh)RNA lentivirus-infected cells, the levels of phosphorylated Akt and MMP-9 were significantly decreased compared with control shRNA lentivirus-infected cells.

CONCLUSION

The present study demonstrated that RRM2 promoted the Akt phosphorylation leading to high expression of MMP-9 to promote the migration and invasive capacities of MH7A cells. Overall, IGF2BP promotes the expression of RRM2, and regulates the migration and invasion of MH7A cells via Akt/MMP-9 pathway to promote RA progression.

ALKBH5 gene polymorphisms and risk of neuroblastoma in Chinese children from Jiangsu Province

Background

Neuroblastoma is one of the most common extracranial malignant solid tumors in children. AlkB homolog 5 (ALKBH5) is an RNA N6-methyladenosine (m6A) demethylase that plays a critical role in tumorigenesis and development. We assessed the association between single nucleotide polymorphisms (SNPs) in ALKBH5 and the risk of neuroblastoma in a case-control study including 402 patients and 473 non-cancer controls.

Methods

Genotyping was determined by the TaqMan method. The association between ALKBH5 polymorphisms (rs1378602 and rs8400) and the risk of neuroblastoma was evaluated using the odds ratio (OR) and 95% confidence interval (CI).

Results

We found no strong association of ALKBH5 rs1378602 and rs8400 with neuroblastoma risk. Further stratification analysis by age, sex, primary site, and clinical stage showed that the rs1378602 AG/AA genotype was associated with a lower risk of neuroblastoma in males (adjusted OR = 0.58, 95% CI = 0.35-0.97, p = 0.036) and children with retroperitoneal neuroblastoma (adjusted OR = 0.58, 95% CI = 0.34-0.98, p = 0.040).

Conclusions

ALKBH5 SNPs do not seem to be associated with neuroblastoma risk. More studies are required to confirm this negative result and reveal the relationship between gene polymorphisms of the m6A modifier ALKBH5 and neuroblastoma.

ALKBH5-Mediated RNA m6 A Methylation Regulates the Migration, Invasion, and Proliferation of Rheumatoid Fibroblast-Like Synoviocytes

OBJECTIVE

Fibroblast-like synoviocytes (FLSs) are critical for promoting joint damage in rheumatoid arthritis (RA). N6 -methyladenosine (m6 A) modification plays key roles in various diseases, but its role in the pathogenesis of RA is largely unknown. Here, we investigate increased demethylase ALKBH5 promotion of proliferation, migration, and invasion of RA FLSs via regulating JARID2 expression.

METHODS

ALKBH5 expression in FLSs was evaluated using real-time quantitative polymerase chain reaction (RT-qPCR) and Western blot. 5-ethynyl-2'-deoxyuridine, scratch wound healing, and transwell assays were implemented to determine the role of ALKBH5 on RA FLS proliferation, mobility, and migration. Then, m6 A sequencing combined with RNA sequencing was performed to identify the potential targets of ALKBH5. RNA immunoprecipitation and RNA pulldown were then used to validate the interaction between the protein and messenger RNA (mRNA). Collagen-induced arthritis (CIA) and delayed-type hypersensitivity arthritis (DTHA) models were further established to assess the therapeutic potency of ALKBH5 in vivo.

RESULTS

We demonstrated that ALKBH5 expression was increased in FLSs and synovium from RA. Functionally, ALKBH5 knockdown inhibited the proliferation, migration, and invasion of RA FLSs, whereas overexpression of ALKBH5 displayed the opposite effect. Mechanistically, ALKBH5 mediated m6 A modification in the JARID2 mRNA and enhanced its mRNA stability in cooperation with IGF2BP3. Intriguingly, the severity of arthritis was attenuated in mice with DTHA and ALKBH5 knockout or rats with CIA and intra-articular injection of ALKBH5 short hairpin RNA.

CONCLUSION

Our findings suggest that ALKBH5-mediated m6 A modification is crucial for synovial hyperplasia and invasion in RA. ALKBH5 might be a potential therapeutic target for RA and even for dysregulated fibroblasts in a wide range of diseases.

The role of m6A RNA methylation in autoimmune diseases: Novel therapeutic opportunities

N6-methyladenosine (m6A) modifications, as one of the most common forms of internal RNA chemical modifications in eukaryotic cells, have gained increasing attention in recent years. The m6A RNA modifications exert various crucial roles in various biological processes, such as embryonic development, neurogenesis, circadian rhythms, and tumorigenesis. Recent advances have highlighted that m6A RNA modification plays an important role in immune response, especially in the initiation and progression of autoimmune diseases. In this review, we summarized the regulatory mechanisms of m6A methylation and its biological functions in the immune system and mainly focused on recent progress in research on the potential role of m6A RNA methylation in the pathogenesis of autoimmune diseases, thus providing possible biomarkers and potential targets for the prevention and treatment of autoimmune diseases.

Genome-wide identification of m6A-associated single nucleotide polymorphisms in complex diseases of nervous system

N6-methyladenosine (m6A) is one of the most abundant chemical modifications on RNA and can affect the occurrence and development of diseases. Some studies have shown that the expressions of some m6A-related genes are significantly regulated by single nucleotide variants (SNV). However, the function of m6A-associated single nucleotide polymorphisms (m6A-SNP) remains unclear in multiple sclerosis (MS), Alzheimer's disease (AD) and Parkinson's disease (PD). Here, we identified the disease-associated m6A-SNPs by integrating genome-wide association study (GWAS) and m6A-SNPs from the RMVar database, and confirmed the relationship between these identified m6A-SNPs and their target genes in eQTL analysis and gene differential expression analysis. Finally, 26 genes corresponding to 20 m6A-SNPs with eQTL signals were identified and differentially expressed (P < 0.05) in MS, 15 genes corresponding to 12 m6A-SNPs (P < 1e-04) were differentially expressed in AD, and 27 PD-associated m6A-SNPs that regulated the expression of 31 genes were identified. There were 5 HLA genes with eQTL signals (HLA-DQB1, HLA-DRB1, HLA-DQA1, HLA-DQA2 and HLA-DQB1-AS1) to be detected in the three diseases. In summary, our study provided new insights into understanding the potential roles of these m6A-SNPs in disease pathogenesis as well as therapeutic target.

ALKBH5-mediated m6A demethylation of HS3ST3B1-IT1 prevents osteoarthritis progression

HS3ST3B1-IT1 was identified as a downregulated long noncoding RNA in osteoarthritic cartilage. However, its roles and mechanisms in the pathogenesis of osteoarthritis (OA) are unclear. In this study, we demonstrated that the expressions of HS3ST3B1-IT1 and its maternal gene HS3ST3B1 were downregulated and positively correlated in osteoarthritic cartilage. Overexpression of HS3ST3B1-IT1 significantly increased chondrocyte viability, inhibited chondrocyte apoptosis, and upregulated extracellular matrix (ECM) proteins, whereas HS3ST3B1-IT1 knockdown had the opposite effects. In addition, HS3ST3B1-IT1 significantly ameliorated monosodium-iodoacetate-induced OA in vivo. Mechanistically, HS3ST3B1-IT1 upregulated HS3ST3B1 expression by blocking its ubiquitination-mediated degradation. Knockdown of HS3ST3B1 reversed the effects of HS3ST3B1-IT1 on chondrocyte viability, apoptosis, and ECM metabolism. AlkB homolog 5 (ALKBH5)-mediated N6-methyladenosine (m6A) demethylation stabilized HS3ST3B1-IT1 RNA. Together, our data revealed that ALKBH5-mediated upregulation of HS3ST3B1-IT1 suppressed OA progression by elevating HS3ST3B1 expression, suggesting that HS3ST3B1-IT1/HS3ST3B1 may serve as potential therapeutic targets for OA treatment.

N6-Methyladenosine Directly Regulates CD40L Expression in CD4+ T Lymphocytes

T cell activation is a highly regulated process, modulated via the expression of various immune regulatory proteins including cytokines, surface receptors and co-stimulatory proteins. N⁶-methyladenosine (m⁶A) is an RNA modification that can directly regulate RNA expression levels and it is associated with various biological processes. However, the function of m⁶A in T cell activation remains incompletely understood. We identify m⁶A as a novel regulator of the expression of the CD40 ligand (CD40L) in human CD4⁺ lymphocytes. Manipulation of the m⁶A 'eraser' fat mass and obesity-associated protein (FTO) and m⁶A 'writer' protein methyltransferase-like 3 (METTL3) directly affects the expression of CD40L. The m⁶A 'reader' protein YT521-B homology domain family-2 (YTHDF2) is hypothesized to be able to recognize and bind m⁶A specific sequences on the CD40L mRNA and promotes its degradation. This study demonstrates that CD40L expression in human primary CD4⁺ T lymphocytes is regulated via m⁶A modifications, elucidating a new regulatory mechanism in CD4⁺ T cell activation that could possibly be leveraged in the future to modulate T cell responses in patients with immune-related diseases.

RNA-binding proteins in degenerative joint diseases: A systematic review

RNA-binding proteins (RBPs), which are conserved proteins comprising multiple intermediate sequences, can interact with proteins, messenger RNA (mRNA) of coding genes, and non-coding RNAs to perform different biological functions, such as the regulation of mRNA stability, selective polyadenylation, and the management of non-coding microRNA (miRNA) synthesis to affect downstream targets. This article will highlight the functions of RBPs, in degenerative joint diseases (intervertebral disc degeneration [IVDD] and osteoarthritis [OA]). It will reviews the latest advancements on the regulatory mechanism of RBPs in degenerative joint diseases, in order to understand the pathophysiology, early diagnosis and treatment of OA and IVDD from a new perspective.

Genome-wide identification of RNA modification-related single nucleotide polymorphisms associated with rheumatoid arthritis

BACKGROUND

RNA modification plays important roles in many biological processes, such as gene expression control. The aim of this study was to identify single nucleotide polymorphisms related to RNA modification (RNAm-SNPs) for rheumatoid arthritis (RA) as putative functional variants.

METHODS

We examined the association of RNAm-SNPs with RA in summary data from a genome-wide association study of 19,234 RA cases and 61,565 controls. We performed eQTL and pQTL analyses for the RNAm-SNPs to find associated gene expression and protein levels. Furthermore, we examined the associations of gene expression and circulating protein levels with RA using two-sample Mendelian randomization analysis methods.

RESULTS

A total of 160 RNAm-SNPs related to m⁶A, m¹A, A-to-I, m⁷G, m⁵C, m⁵U and m⁶Am modifications were identified to be significantly associated with RA. These RNAm-SNPs were located in 62 protein-coding genes, which were significantly enriched in immune-related pathways. RNAm-SNPs in important RA susceptibility genes, such as PADI2, SPRED2, PLCL2, HLA-A, HLA-B, HLA-DRB1, HLA-DPB1, TRAF1 and TXNDC11, were identified. Most of these RNAm-SNPs showed eQTL effects, and the expression levels of 26 of the modifiable genes (e.g., PADI2, TRAF1, HLA-A, HLA-DRB1, HLA-DPB1 and HLA-B) in blood cells were associated with RA. Circulating protein levels, such as CFB, GZMA, HLA-DQA2, IL21, LRPAP1 and TFF3, were affected by RNAm-SNPs and were associated with RA.

CONCLUSION

The present study identified RNAm-SNPs in the reported RA susceptibility genes and suggested that RNAm-SNPs may affect RA risk by affecting the expression levels of corresponding genes and proteins.

In silico identification and in vitro expression analysis of breast cancer-related m6A-SNPs

Research on m6A-associated SNPs (m6A-SNPs) has emerged recently due to their possible critical roles in many key biological processes. In this sense, several investigations have identified m6A-SNPs in different diseases. In order to gain a more complete understanding of the role that m6A-SNPs can play in breast cancer, we performed an in silico analysis to identify the m6A-SNPs associated with breast cancer and to evaluate their possible effects. For this purpose, we downloaded SNPs related to breast cancer and a list of m6A-SNPs from public databases in order to identify which ones appear in both. Subsequently, we assessed the identified m6A-SNPs in silico by expression quantitative trait loci (eQTL) analysis and differential gene expression analysis. We genotyped the m6A-SNPs found in the in silico analysis in 35 patients with breast cancer, and we carried out a gene expression analysis experimentally on those that showed differences. Our results identified 981 m6A-SNPs related to breast cancer. Four m6A-SNPs showed an eQTL effect and only three were in genes that presented an altered gene expression. When the three m6A-SNPs were evaluated in the tissue sample of our breast cancer patients, only the m6A-SNP rs76563149 located in ZNF354A gene presented differences in allele frequencies and a low gene expression in breast cancer tissues, especially in luminal B HER2+ subtype. Future investigations of these m6A-SNPs should expand the study in different ethnic groups and increase the sample sizes to test their association with breast cancer and elucidate their molecular function.

Comprehensive Analysis and Functional Characteristics of Differential Expression of N6-Methyladenosine Methylation Modification in the Whole Transcriptome of Rheumatoid Arthritis

N6-methyladenosine (m6A) modification is the most prevalent chemical modification in eukaryotic mRNA and is associated with the development of various immune diseases. However, the role of m6A methylation in rheumatoid arthritis (RA) development is unclear. We preliminarily explored the role of m6A methylation-related mRNAs in RA for its clinical application. The discovery of m6A methylation-modifying genes in this study may provide a fresh perspective on the development of drugs for RA treatment. High-throughput sequencing combined with methylated RNA immunoprecipitation (MeRIP-seq) and RNA sequencing were used to assess whole-transcriptome m6A modifications in the synovium of patients with RA. The relationship between m6A-modified target genes and RA inflammation and macrophages was determined. The expression of the m6A-modified significant transcript-enriched inflammatory signaling pathway was assessed through animal experiments. Differentially expressed m6A genes were correlated with macrophage activation involved in immune response, vascular endothelium, MAPK signaling pathway, PI3K - Akt signaling pathway, and other inflammatory processes. Furthermore, combined analysis with m6A-seq and RNA-seq revealed 120 genes with significant changes in both m6A modification and mRNA expression. We selected the top 3 candidate mRNAs that were upregulated and downregulated simultaneously. The expression of phosphatase and tensin homolog deleted on chromosome ten (PTEN) mRNA and protein in RA patients was lower than that in healthy control (HC). SHC-binding protein 1 (SHCBP1) and neurexophilin-3 (NXPH3) mRNA expressions were increased in RA patients. The expression of M1 macrophages was increased in RA patients. RA markers are such as rheumatoid factor (RF) and peptide containing citrulline (CCP). Further animal experiments showed that the expression of synovial MAPK, PI3K, and Akt1 proteins in the RA model was increased, and the PTEN, p-PTEN protein expression was decreased. PI3K, Akt1, PTEN, and p-PTEN were correlated to RA joint inflammation. This study revealed a unique pattern of differential m6A methylation modifications in RA and concluded that m6A modification is related to the occurrence of RA synovial inflammation.

Alterations of m6A RNA methylation regulators contribute to autophagy and immune infiltration in primary Sjögren's syndrome

N6-methyladenosine (m6A) RNA modification is a new epigenetic regulation mechanism on eukaryotic mRNA. Few autoimmune diseases focused on the role of m6A in their pathogenies, and m6A modulation in the pathological process of primary Sjögren's syndrome (pSS) is still unknown. In this work, three microarray datasets of pSS patients were downloaded from the GEO database: datasets #1 and #2 from the whole peripheral blood (PB) samples, dataset #3 from the labial salivary gland tissue samples, as well as a PB cohort collected from our hospital. Six differentially expressed m6A regulators were identified by comparing the PB dataset #1 of pSS and healthy controls using the Wilcox test and logistic regression analysis. Among them, four (ALKBH5, RBMX, RBM15B, and YTHDF1) were confirmed as down-regulated in PB dataset #2 and in our PB cohort by RT-PCR, and four (ALKBH5, METTL3, RBM15B, and YTHDF1) were confirmed as down-regulated in the dataset #3 of the labial gland tissue. In addition, discrepantly expressed m6A regulators accompanied by diverse immunocytes, including dendritic cells (DCs), T cells, and CD56dim natural killer cells, and among the regulators, ALKBH5 and METTL3 were comprehensively linked with the infiltrated immune cells. Notably, the most enriched autophagy mechanism mediated by m6A was observed in pSS using functional annotation analysis. Ten hub genes were identified using a protein-protein interaction network, and their expression in PB dataset #2 and the expression of three genes (PIK3CA, STAT1, and MAPK3) in the labial gland tissue dataset #3 were confirmed. Our study provides evidence that m6A methylation is widely involved in the immune infiltration and autophagy of pSS, thus contributing to the pathogenesis of this disease and potentially representing a novel therapeutic target.

Predictive value of drug efficacy by M6A modification patterns in rheumatoid arthritis patients

Background

Rheumatoid arthritis is a highly heterogeneous autoimmune disease characterized by unpredictable disease flares and significant differences in therapeutic response to available treatments. One possible reason for poor efficacy is that it cannot be treated accurately due to no optimal stratification for RA patients.

Objective

This study aims to construct an RA classification model by m6A characters and further predict response to medication.

Methods

Twenty m6A regulators were used to construct a random forest diagnosis model, and RNA-seq analysis was employed for external validation. The RNA modification patterns mediated by 20 m6A regulators were systematically evaluated in 1191 RA samples and explored different molecular clusters associated with other immune microenvironment characteristics and biological pathways. Then, we established an m6A score model to quantify the m6A modification patterns. The model was applied to patients at baseline to test the association between m6Ascore and infliximab responsiveness.

Results

The m6A diagnosis model showed good discriminatory ability in distinguishing RA. Patients with RA were classified into three clusters with distinct molecular and cellular signatures. Cluster A displayed strongly activated inflammatory cells and pathways. Specific innate lymphocytes occupied cluster B. Cluster C was mainly enriched in prominent adaptive lymphocytes and NK-mediated cytotoxicity signatures with the highest m6A score. Patients with a low m6Ascore exhibited significantly infliximab therapeutic benefits compared with those with a high m6Ascore (p< 0.05).

Conclusion

Our study is the first to provide a comprehensive analysis of m6A modifications in RA, which provides an innovative patient stratification framework and potentially enables improved therapeutic decisions.

Extracellular vesicles derived from human umbilical cord mesenchymal stem cells alleviate osteoarthritis of the knee in mice model by interacting with METTL3 to reduce m6A of NLRP3 in macrophage

BACKGROUND

Osteoarthritis (OA) is a prevalent degenerative joint disease that not only significantly impairs the quality of life of middle-aged and elderly individuals but also imposes a significant financial burden on patients and society. Due to their significant biological properties, extracellular vesicles (EVs) have steadily received great attention in OA treatment. This study aimed to investigate the influence of EVs on chondrocyte proliferation, migration, and apoptosis and their protective efficacy against OA in mice.

METHODS

The protective impact of EVs derived from human umbilical cord mesenchymal stem cells (hucMSCs-EVs) on OA in mice was investigated by establishing a mouse OA model by surgically destabilizing the medial meniscus (DMM). Human chondrocytes were isolated from the cartilage of patients undergoing total knee arthroplasty (TKA) and cultured with THP-1 cells to mimic the in vivo inflammatory environment. Levels of inflammatory factors were then determined in different groups, and the impacts of EVs on chondrocyte proliferation, migration, apoptosis, and cartilage extracellular matrix (ECM) metabolism were explored. N6-methyladenosine (m6A) level of mRNA and methyltransferase-like 3 (METTL3) protein expression in the cells was also measured in addition to microRNA analysis to elucidate the molecular mechanism of exosomal therapy.

RESULTS

The results indicated that hucMSCs-EVs slowed OA progression, decreased osteophyte production, increased COL2A1 and Aggrecan expression, and inhibited ADAMTS5 and MMP13 overexpression in the knee joint of mice via decreasing pro-inflammatory factor secretion. The in vitro cell line analysis revealed that EVs enhanced chondrocyte proliferation and migration while inhibiting apoptosis. METTL3 is responsible for these protective effects. Further investigations revealed that EVs decreased the m6A level of NLRP3 mRNA following miR-1208 targeted binding to METTL3, resulting in decreased inflammatory factor release and preventing OA progression.

CONCLUSION

This study concluded that hucMSCs-EVs inhibited the secretion of pro-inflammatory factors and the degradation of cartilage ECM after lowering the m6A level of NLRP3 mRNA with miR-1208 targeting combined with METTL3, thereby alleviating OA progression in mice and providing a novel therapy for clinical OA treatment.

The crucial mechanism and therapeutic implication of RNA methylation in bone pathophysiology

Methylation is the most common posttranscriptional modification in cellular RNAs, which has been reported to modulate the alteration of RNA structure for initiating relevant functions such as nuclear translocation and RNA degradation. Recent studies found that RNA methylation especially N6-methyladenosine (m⁶A) regulates the dynamic balance of bone matrix and forms a complicated network in bone metabolism. The modulation disorder of RNA methylation contributes to several pathological bone diseases including osteoporosis (OP), osteoarthritis (OA), rheumatoid arthritis (RA), and so on. In the review, we will discuss advanced technologies for detecting RNA methylation, summarize RNA methylation-related biological impacts on regulating bone homeostasis and pathological bone diseases. In addition, we focus on the promising roles of RNA methylation in early diagnosis and therapeutic implications for bone-related diseases. Then, we aim to establish a theoretical basis for further investigation in this meaningful field.

N(6)-methyladenosine RNA methyltransferase like 3 inhibits extracellular matrix synthesis of endplate chondrocytes by downregulating sex-determining region Y-Box transcription factor 9 expression under tension

OBJECTIVE

Tension stimulation is an important inducer of endplate cartilage degeneration, but the specific regulatory mechanism remains unclear. This study was the first to reveal the mechanism by which methyltransferase-like 3 (METTL3)-mediated N(6)-methyladenosine (m6A) modification affected the extracellular matrix anabolism by tension-induced endplate chondrocytes.

METHOD

We examined the differences in METTL3 expression and m6A methylation levels in human endplate chondrocytes and human cartilage endplate tissues under in vitro tension. The effect on endplate cartilage degeneration was evaluated by manipulating m6A methylation mediated by METTL3 in vivo and in vitro. The effect of METTL3-mediated m6A methylation on the stability of sex-determining region Y-box transcription factor 9 (SOX9) gene expression was determined experimentally.

RESULTS

METTL3 expression and m6A methylation levels were significantly increased in degenerative human endplate cartilage tissue. Similarly, tension stimulation inhibited the ability of human endplate chondrocytes to synthesize extracellular matrix, which was accompanied by an increase in METTL3-mediated m6A methylation. The ability of endplate chondrocytes to resist tension was significantly enhanced by inhibiting METTL3 expression and subsequently downregulating m6A methylation in vitro and in vivo, thereby reducing intervertebral disc degeneration. Furthermore, METTL3 mediated SOX9 RNA methylation and disrupted SOX9 mRNA stability, thereby inhibiting the gene expression of the downstream collagen type II alpha 1 chain.

CONCLUSION

Tension stimulation downregulated SOX9 expression through METTL3-mediated m6A methylation, thereby inhibiting the synthesis of extracellular matrix in endplate chondrocytes.

The roles of N6-methyladenosine methylation in the regulation of bone development, bone remodeling and osteoporosis

N⁶-methyladenosine (m⁶A), a novel epitranscriptomic RNA modification, plays crucial roles in a variety of biological processes and diseases. Recently, there are growing evidence supporting that m⁶A methylation is essential for bone development and homeostasis through the regulation of key genes by regulating RNA stability, localization, turnover and translation efficiency. In this review, we summarized our current understanding of the functional roles of m⁶A methylation and its related regulators in bone development and bone remodeling. These findings will offer new directions and insights on the further investigations of m⁶A methylation in bone biology. Moreover, we also discussed important advances of m⁶A methylation related regulators as potential therapeutic targets, which allows for novel therapeutic strategies on the medications of bone-related diseases including osteoporosis.

The molecular structure and biological functions of RNA methylation, with special emphasis on the roles of RNA methylation in autoimmune diseases

Autoimmune diseases such as rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), and systemic vasculitis are caused by the body's immune response to autoantigens. The pathogenesis of autoimmune diseases is complex. RNA methylation is known to play a key role in disease progression as it regulates almost all aspects of RNA processing, including RNA nuclear export, translation, splicing, and noncoding RNA processing. This review summarizes the mechanisms, molecular structures of RNA methylations and their roles in biological functions. Similar to the roles of RNA methylation in cancers, RNA methylation in RA and SLE involves "writers" that deposit methyl groups to form N6-methyladenosine (m6A) and 5-methylcytosine (m5C), "erasers" that remove these modifications, and "readers" that further affect mRNA splicing, export, translation, and degradation. Recent advances in detection methods have identified N1-methyladenosine (m1A), N6,2-O-dimethyladenosine (m6Am), and 7-methylguanosine (m7G) RNA modifications, and their roles in RA and SLE need to be further studied. The relationship between RNA methylation and other autoimmune diseases has not been reported, and the roles and mechanisms of RNA modifications in these diseases need to be explored in the future.

m6A-mediated upregulation of AC008 promotes osteoarthritis progression through the miR-328-3p‒AQP1/ANKH axis

Long noncoding RNAs (lncRNAs) have emerged as important regulators of osteoarthritis (OA), but the biological roles and clinical significance of most lncRNAs in OA are not fully understood. Microarray analysis was performed to identify differentially expressed lncRNAs, mRNAs, and miRNAs between normal and osteoarthritic cartilage. We found that AC008440.5 (abbreviated AC008), as well as AQP1 and ANKH, were highly expressed in osteoarthritic cartilage, whereas miR-328-3p was expressed at a low level in osteoarthritic cartilage. Functional assays showed that ectopic expression of AC008, AQP1, and ANKH significantly decreased chondrocyte viability and promoted chondrocyte apoptosis and extracellular matrix (ECM) degradation, whereas knockdown of AC008, AQP1, and ANKH resulted in the opposite effects. Moreover, miR-328-3p overexpression increased chondrocyte viability and attenuated chondrocyte apoptosis and ECM degradation, whereas inhibition of miR-328-3p resulted in the opposite effects. Bioinformatics analysis, RNA immunoprecipitation (RIP), and luciferase assays revealed that AC008 functioned as a competing endogenous RNA (ceRNA) to regulate miR-328-3p, which specifically targeted the AQP1 and ANKH genes. In addition, miR-328-3p significantly ameliorated MIA-induced OA, whereas AC008 accelerated OA progression in vivo. Furthermore, fat mass and obesity-associated (FTO)-mediated N6-methyladenosine demethylation downregulated AC008 transcription, while lower FTO expression led to upregulation of AC008 transcription in OA. In conclusion, our data reveal that AC008 plays a critical role in OA pathogenesis via the miR-328-3p‒AQP1/ANKH pathway, suggesting that AC008 may be a potential therapeutic target for OA.

Role of m6A in osteoporosis, arthritis and osteosarcoma (Review)

RNA modification is a type of post-transcriptional modification that regulates important cellular pathways, such as the processing and metabolism of RNA. The most abundant form of methylation modification is RNA N6-methyladenine (m6A), which plays various post-transcriptional regulatory roles in cellular biological functions, including cell differentiation, embryonic development and disease occurrence. Bones play a pivotal role in the skeletal system as they support and protect muscles and other organs, facilitate movement and ensure haematopoiesis. The development and remodelling of bones require a delicate and accurate regulation of gene expression by epigenetic mechanisms that involve modifications of histone, DNA and RNA. The present review discusses the enzymes and proteins involved in mRNA m6A methylation modification and summarises current research progress and the mechanisms of mRNA m6A methylation in common orthopaedic diseases, including osteoporosis, arthritis and osteosarcoma.

METTL3 involves the progression of osteoarthritis probably by affecting ECM degradation and regulating the inflammatory response

We aimed to identify RNA N6-methyladenosine methylation associated genes in osteoarthritis (OA), and to explore possible regulatory mechanisms of these RNA methylation associated genes. Bioinformatics analyses, including differential expression analysis, functional enrichment analysis, verification analysis, and box plot analysis, were conducted based on different datasets from OA and non-OA patients. Gene expression at mRNA and protein levels was determined by quantitative reverse transcription PCR, western blot and immunofluorescence. Interleukin 1β (IL-1β)-treated SW1353 cells was used as cell model. Lentiviral vector was used for over-expression METTL3 in vitro. CCK-8 assay kit was used to determine cell viability and inflammatory cytokines (IL-1α, IL-6, IL-8, IL-10 and TNF-α) was detected using ELISA kits. Bioinformatics analysis showed that METTL3 expression was decreased in OA group, which was confirmed in clinical samples. Expression of METTL3 was also reduced in IL-1β-treated cells. Levels of inflammatory cytokines were obviously reduced in the METTL3 overexpression group, while IL-1β treatment reversed such decrease caused by METTL3 overexpression (p < 0.05). Both METTL3 overexpression and IL-1β treatment promoted expression of p65 protein and p-ERK (p < 0.01). Additionally, increased expression of MMP1 and MMP3, and decreased expression of MMP13, TIMP-1, and TIMP-2 at both mRNA and protein levels were observed in the METTL3 overexpression group when compared with the control group (p < 0.01). Expression of m6A methylation gene METTL3 was reduced in OA. METTL3 is involved in OA probably by regulating the inflammatory response. METTL3 overexpression may affect extracellular matrix degradation in OA by adjusting the balance between TIMPs and MMPs.

Modeling multi-species RNA modification through multi-task curriculum learning

N6-methyladenosine (m6A) is the most pervasive modification in eukaryotic mRNAs. Numerous biological processes are regulated by this critical post-transcriptional mark, such as gene expression, RNA stability, RNA structure and translation. Recently, various experimental techniques and computational methods have been developed to characterize the transcriptome-wide landscapes of m6A modification for understanding its underlying mechanisms and functions in mRNA regulation. However, the experimental techniques are generally costly and time-consuming, while the existing computational models are usually designed only for m6A site prediction in a single-species and have significant limitations in accuracy, interpretability and generalizability. Here, we propose a highly interpretable computational framework, called MASS, based on a multi-task curriculum learning strategy to capture m6A features across multiple species simultaneously. Extensive computational experiments demonstrate the superior performances of MASS when compared to the state-of-the-art prediction methods. Furthermore, the contextual sequence features of m6A captured by MASS can be explained by the known critical binding motifs of the related RNA-binding proteins, which also help elucidate the similarity and difference among m6A features across species. In addition, based on the predicted m6A profiles, we further delineate the relationships between m6A and various properties of gene regulation, including gene expression, RNA stability, translation, RNA structure and histone modification. In summary, MASS may serve as a useful tool for characterizing m6A modification and studying its regulatory code. The source code of MASS can be downloaded from https://github.com/mlcb-thu/MASS.

m6A RNA Methylation in Systemic Autoimmune Diseases-A New Target for Epigenetic-Based Therapy?

The general background of autoimmune diseases is a combination of genetic, epigenetic and environmental factors, that lead to defective immune reactions. This erroneous immune cell activation results in an excessive production of autoantibodies and prolonged inflammation. During recent years epigenetic mechanisms have been extensively studied as potential culprits of autoreactivity. Alike DNA and proteins, also RNA molecules are subjected to an extensive repertoire of chemical modifications. N6-methyladenosine is the most prevalent form of internal mRNA modification in eukaryotic cells and attracts increasing attention due to its contribution to human health and disease. Even though m6A is confirmed as an essential player in immune response, little is known about its role in autoimmunity. Only few data have been published up to date in the field of RNA methylome. Moreover, only selected autoimmune diseases have been studied in respect of m6A role in their pathogenesis. In this review, I attempt to present all available research data regarding m6A alterations in autoimmune disorders and appraise its role as a potential target for epigenetic-based therapies.

Emerging Perspectives of RNA N 6-methyladenosine (m6A) Modification on Immunity and Autoimmune Diseases

N ⁶-methyladenosine (m⁶A) modification, the addition of a methylation decoration at the position of N6 of adenosine, is one of the most prevalent modifications among the over 100 known chemical modifications of RNA. Numerous studies have recently characterized that RNA m⁶A modification functions as a critical post-transcriptional regulator of gene expression through modulating various aspects of RNA metabolism. In this review, we will illustrate the current perspectives on the biological process of m⁶A methylation. Then we will further summarize the vital modulatory effects of m⁶A modification on immunity, viral infection, and autoinflammatory disorders. Recent studies suggest that m⁶A decoration plays an important role in immunity, viral infection, and autoimmune diseases, thereby providing promising biomarkers and therapeutic targets for viral infection and autoimmune disorders.

Regulatory Role of the RNA N6-Methyladenosine Modification in Immunoregulatory Cells and Immune-Related Bone Homeostasis Associated With Rheumatoid Arthritis

Rheumatoid arthritis (RA) is a systemic autoimmune disease for which the etiology has not been fully elucidated. Previous studies have shown that the development of RA has genetic and epigenetic components. As one of the most highly abundant RNA modifications, the N⁶-methyladenosine (m⁶A) modification is necessary for the biogenesis and functioning of RNA, and modification aberrancies are associated with various diseases. However, the specific functions of m⁶A in the cellular processes of RA remain unclear. Recent studies have revealed the relationship between m⁶A modification and immune cells associated with RA. Therefore, in this review, we focused on discussing the functions of m⁶A modification in the regulation of immune cells and immune-related bone homeostasis associated with RA. In addition, to gain a better understanding of the progress in this field of study and provide the proper direction and suggestions for further study, clinical application studies of m⁶A modification were also summarized.

METTL3 promotes IL-1β-induced degeneration of endplate chondrocytes by driving m6A-dependent maturation of miR-126-5p

METTL3 is an important regulatory molecule in the process of RNA biosynthesis. It mainly regulates mRNA translation, alternative splicing and microRNA maturation by mediating m6A‐dependent methylation. Interleukin 1β (IL‐1β) is an important inducer of cartilage degeneration that can induce an inflammatory cascade reaction in chondrocytes and inhibit the normal biological function of cells. However, it is unclear whether IL‐1β is related to METTL3 expression or plays a regulatory role in endplate cartilage degeneration. In this study, we found that the expression level of METTL3 and methylation level of m6A in human endplate cartilage with different degrees of degeneration were significantly different, indicating that the methylation modification of m6A mediated by METTL3 was closely related to the degeneration of human endplate cartilage. Next, through a series of functional experiments, we found that miR‐126‐5p can play a significant role in IL‐1β–induced degeneration of endplate chondrocytes. Moreover, we found that miR‐126‐5p can inhibit the PI3K/Akt signalling pathway by targeting PIK3R2 gene, leading to the disorder of cell vitality and functional metabolism. To further determine whether METTL3 could regulate miR‐126‐5p maturation, we first confirmed that METTL3 can bind the key protein underlying pri‐miRNA processing, DGCR8. Additionally, when METTL3 expression was inhibited, the miR‐126‐5p maturation process was blocked. Therefore, we hypothesized that METTL3 can promote cleavage of pri‐miR‐126‐5p and form mature miR‐126‐5p by combining with DGCR8.

Downregulation of m6 A reader YTHDC2 promotes tumor progression and predicts poor prognosis in non-small cell lung cancer

Background

m⁶A modification affects the pathological progress of many diseases by affecting RNA stability and translocation. YTHDC2, a m⁶A reader, is associated with multiple cancers; however, little is known of its role in non‐small cell lung cancer (NSCLC).

Methods

The GEPIA, Oncomine and GEO databases were analyzed to assess expression of YTHDC2 in NSCLC patients. Quantitative polymerase chain reaction, western blot and immunohistochemistry were used to detect YTHDC2 expression in different NSCLC cell lines (H1299, H460, H292 and A549) and patients. The effects of YTHDC2 on NSCLC cell lines (A549 and H1299) proliferation and migration were employed using CCK8 and transwell assays. The potential target RNAs of YTHDC2 were obtained from the POSTAR database. Functional enrichment analysis of YTHDC2 targeted RNAs was performed using the Metascape database.

Results

GEPIA, Oncomine and GEO databases showed low expression of YTHDC2 in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) patients. YTHDC2 expression was significantly decreased in different NSCLC cell lines and our clinical samples. Moreover, low expression of YTHDC2 was significantly associated with poor differentiation, lymph node metastasis, tumor size and stage. In addition, YTHDC2 could suppress the proliferation and migration ability of A549 and H1299 cell lines. Kaplan‐Meier Plotter database analysis revealed that patients with low level of YTHDC2 had a significantly poor prognosis. Finally, functional enrichment analysis of YTHDC2 targeted RNAs indicated several enriched pathways related to cancer.

Conclusions

These findings elucidate that YTHDC2 suppresses tumorigenesis in NSCLC, indicating that YTHDC2 may be a promising therapeutic target for NSCLC.

Key points

Significant findings of the study

This study demonstrated that the downregulation of YTHDC2 promotes tumor progression and predicts poor prognosis in non‐small cell lung cancer (NSCLC).

What this study adds

YTHDC2 might be a promising therapeutic target for non‐small cell lung cancer (NSCLC).

Multifaceted Functions and Novel Insight Into the Regulatory Role of RNA N6-Methyladenosine Modification in Musculoskeletal Disorders

RNA modifications have emerged as key regulators of transcript expression in diverse physiological and pathological processes. As one of the most prevalent types of RNA modifications, N⁶-methyladenosine (m⁶A) has become the highlight in modulation of various diseases through interfering RNA splicing, translation, nuclear export, and decay. In many cases, the detailed functions of m⁶A in cellular processes and diseases remain unclear. Notably, recent studies have determined the relationship between m⁶A modification and musculoskeletal disorders containing osteosarcoma, osteoarthritis, rheumatoid arthritis, osteoporosis, etc. Herein, this review comprehensively summarizes the recent advances of m⁶A modification in pathogenesis and progression of musculoskeletal diseases. Specifically, the underlying molecular mechanisms, detection technologies, regulatory functions, clinical implications, and future perspectives of m⁶A in musculoskeletal disorders are discussed, with the aim to provide a novel insight into their association.

ALKBH5 gene polymorphisms and Wilms tumor risk in Chinese children: A five-center case-control study

BACKGROUND

Wilms tumor is a frequently diagnosed renal cancer among children with unclear genetic causes. N6-methyladenosine (m⁶ A) modification genes play critical roles in tumorigenesis. However, whether genetic variations of m⁶ A modification genes predispose to Wilms tumor remain unclear. ALKBH5 (AlkB homolog 5), a crucial member of m⁶ A modification genes, encodes a demethylase that functions to reverse m⁶ A RNA methylation.

METHODS

Herein, we evaluated the association of single nucleotide polymorphisms (SNPs) in the m⁶ A modification gene ALKBH5 and Wilms tumor susceptibility in a large multi-center case-control study. A total of 414 Wilms tumor cases and 1199 healthy controls were genotyped for ALKBH5 rs1378602 and rs8400 polymorphisms by TaqMan.

RESULTS

No significant association was detected between these two polymorphisms and Wilms tumor risk. Moreover, 1, 2, and 1-2 protective genotypes (rs1378602 AG/AA or rs8400 GG) did not significantly reduce Wilms tumor risk, compared with risk genotypes only. Stratification analysis revealed a significant relationship between rs1378602 AG/AA genotypes and decreased Wilms tumor risk in children in clinical stage I diseases [adjusted odds ratio (OR) = 0.56, 95% confidence interval (CI) = 0.32-0.98, P = .042]. The presence of 1-2 protective genotypes was correlated with decreased Wilms tumor risk in subgroups of age > 18 months, when compared to the absence of protective genotypes (adjusted OR = 0.74, 95% CI = 0.56-0.98, P = .035).

CONCLUSION

Collectively, our results demonstrate that ALKBH5 SNPs may exert a weak influence on susceptibility to Wilms tumor. This finding increases the understanding of the role of the m⁶ A gene in tumorigenesis of Wilms tumor.

Integration Analysis of m6A-SNPs and eQTLs Associated With Sepsis Reveals Platelet Degranulation and Staphylococcus aureus Infection are Mediated by m6A mRNA Methylation

Sepsis is a major threat with high mortality rate for critically ill patients. Response to pathogen infection by the host immune system is a key biological process involved in the onset and development of sepsis. Heterogeneous host genome variation, especially single nucleotide polymorphisms (SNPs), has long been suggested to contribute to differences in disease progression. However, the function of SNPs located in non-coding regions remains to be elucidated. Recently, m⁶A mRNA modification levels were revealed to differ at SNPs. As m⁶A is a crucial regulator of gene expression, these SNPs might control genes by changing the m⁶A level on mRNA. To investigate the potential role of m⁶A SNPs in sepsis, we integrated m⁶A-SNP and expression quantitative trait loci (eQTLs) data. Analysis revealed 15,720 m⁶A-cis-eQTLs and 381 m⁶A-trans-eQTLs associated with sepsis. We identified 1321 genes as locations of m⁶A-cis-eQTLs. These were enriched in platelet degranulation and Staphylococcus aureus infection pathways, which are vital for the pathophysiological process of sepsis. We conclude that m⁶A modification of mRNA plays a very important role in sepsis, with m⁶A-cis-eQTLs potentially having the most effect on individual variation in sepsis progression.

METTL3 promotes experimental osteoarthritis development by regulating inflammatory response and apoptosis in chondrocyte

OBJECTIVE

This study was to investigate the functional role of RNA methyltransferase METTL3, an enzyme catalyzes the formation of N6-methyladenosine (m6A) on the target mRNA, in the development of osteoarthritis (OA) and the underlying mechanism.

METHODS

Cytokine IL-1β was used to stimulate the chondroprogenitor cell line ATDC5 cells to mimic the inflammatory condition in vitro. The level of METTL3 mRNA and m6A as well as inflammatory cytokines were detected by qRT-PCR. Cell activity was detected by CCK-8. The rate of apoptotic cell was measured by flow cytometry. Western blot was used to detect the levels of NF-κB signaling molecules and collagen in cells. Methylation inhibitor cycloleucine and methyl donor betaine were used to treat collagenase-induced OA mice.

RESULTS

In IL-1β-treated ATDC5 cells, the METTL3 mRNA levels and the percentage of m6A methylated mRNA of total mRNA were increased in a dose-dependent manner. Silencing of METTL3 by shRNA reduced the percentage of IL-1β-induced apoptosis, suppressed IL-1β-induced increased inflammatory cytokines levels and activation of NF-κB signaling in chondrocytes. Moreover, silencing of METTL3 promotes degradation of extracellular matrix (ECM) by reducing the expression of MMP-13 and Coll X, elevating the expression of Aggrecan and Coll II. In a OA mouse model induced by collagenase, injection of methylation inhibitor cycloleucine or methyl donor betaine does not affects METTL3 mRNA expression, but significantly inhibits or promotes the total level of m6A as well as inflammatory condition and ECM degradation, respectively.

CONCLUSION

METTL3 has a functional role in mediates osteoarthritis progression by regulating NF-κB signaling and ECM synthesis in chondrocytes that shed insight on developing preventive and curative strategies for OA by focusing on METTL3 and mRNA methylation.

DRUM: Inference of Disease-Associated m6A RNA Methylation Sites From a Multi-Layer Heterogeneous Network

Recent studies have revealed that the RNA N 6-methyladenosine (m6A) modification plays a critical role in a variety of biological processes and associated with multiple diseases including cancers. Till this day, transcriptome-wide m6A RNA methylation sites have been identified by high-throughput sequencing technique combined with computational methods, and the information is publicly available in a few bioinformatics databases; however, the association between individual m6A sites and various diseases are still largely unknown. There are yet computational approaches developed for investigating potential association between individual m6A sites and diseases, which represents a major challenge in the epitranscriptome analysis. Thus, to infer the disease-related m6A sites, we implemented a novel multi-layer heterogeneous network-based approach, which incorporates the associations among diseases, genes and m6A RNA methylation sites from gene expression, RNA methylation and disease similarities data with the Random Walk with Restart (RWR) algorithm. To evaluate the performance of the proposed approach, a ten-fold cross validation is performed, in which our approach achieved a reasonable good performance (overall AUC: 0.827, average AUC 0.867), higher than a hypergeometric test-based approach (overall AUC: 0.7333 and average AUC: 0.723) and a random predictor (overall AUC: 0.550 and average AUC: 0.486). Additionally, we show that a number of predicted cancer-associated m6A sites are supported by existing literatures, suggesting that the proposed approach can effectively uncover the underlying epitranscriptome circuits of disease mechanisms. An online database DRUM, which stands for disease-associated ribonucleic acid methylation, was built to support the query of disease-associated RNA m6A methylation sites, and is freely available at: www.xjtlu.edu.cn/biologicalsciences/drum.

Comparative study for haplotype block partitioning methods - Evidence from chromosome 6 of the North American Rheumatoid Arthritis Consortium (NARAC) dataset

Haplotype-based methods compete with “one-SNP-at-a-time” approaches on being preferred for association studies. Chromosome 6 contains most of the known genetic biomarkers for rheumatoid arthritis (RA) disease. Therefore, chromosome 6 serves as a benchmark for the haplotype methods testing. The aim of this study is to test the North American Rheumatoid Arthritis Consortium (NARAC) dataset to find out if haplotype block methods or single-locus approaches alone can sufficiently provide the significant single nucleotide polymorphisms (SNPs) associated with RA. In addition, could we be satisfied with only one method of the haplotype block methods for partitioning chromosome 6 of the NARAC dataset? In the NARAC dataset, chromosome 6 comprises 35,574 SNPs for 2,062 individuals (868 cases, 1,194 controls). Individual SNP approach and three haplotype block methods were applied to the NARAC dataset to identify the RA biomarkers. We employed three haplotype partitioning methods which are confidence interval test (CIT), four gamete test (FGT), and solid spine of linkage disequilibrium (SSLD). P-values after stringent Bonferroni correction for multiple testing were measured to assess the strength of association between the genetic variants and RA susceptibility. Moreover, the block size (in base pairs (bp) and number of SNPs included), number of blocks, percentage of uncovered SNPs by the block method, percentage of significant blocks from the total number of blocks, number of significant haplotypes and SNPs were used to compare among the three haplotype block methods. Individual SNP, CIT, FGT, and SSLD methods detected 432, 1,086, 1,099, and 1,322 associated SNPs, respectively. Each method identified significant SNPs that were not detected by any other method (Individual SNP: 12, FGT: 37, CIT: 55, and SSLD: 189 SNPs). 916 SNPs were discovered by all the three haplotype block methods. 367 SNPs were discovered by the haplotype block methods and the individual SNP approach. The P-values of these 367 SNPs were lower than those of the SNPs uniquely detected by only one method. The 367 SNPs detected by all the methods represent promising candidates for RA susceptibility. They should be further investigated for the European population. A hybrid technique including the four methods should be applied to detect the significant SNPs associated with RA for chromosome 6 of the NARAC dataset. Moreover, SSLD method may be preferred for its favored benefits in case of selecting only one method.

Regulatory Role of RNA N6-Methyladenosine Modification in Bone Biology and Osteoporosis

Osteoporosis is a metabolic skeletal disorder in which bone mass is depleted and bone structure is destroyed to the degree that bone becomes fragile and prone to fractures. Emerging evidence suggests that N⁶-methyladenosine (m⁶A) modification, a novel epitranscriptomic marker, has a significant role in bone development and metabolism. M⁶A modification not only participates in bone development, but also plays important roles as writers and erasers in the osteoporosis. M⁶A methyltransferase METTL3 and demethyltransferase FTO involves in the delicate process between adipogenesis differentiation and osteogenic differentiation, which is important for the pathological development of osteoporosis. Conditional knockdown of the METTL3 in bone marrow stem cells (BMSCs) could suppress PI3K-Akt signaling, limit the expression of bone formation-related genes (such as Runx2 and Osterix), restrain the expression of vascular endothelial growth factor (VEGF) and down-regulate the decreased translation efficiency of parathyroid hormone receptor-1 mRNA. Meanwhile, knockdown of the METTL3 significantly promoted the adipogenesis process and janus kinase 1 (JAK1) protein expression via an m⁶A-dependent way. Specifically, there was a negative correlation between METTL3 expression and porcine BMSCs adipogenesis. The evidence above suggested that the relationship between METTL3 expression and adipogenesis was inverse, and osteogenesis was positive, respectively. Similarly, FTO regulated for BMSCs fate determination during osteoporosis through the GDF11-FTO-PPARγ axis, prompting the shift of MSC lineage commitment to adipocyte and inhibiting bone formation during osteoporosis. In this systematic review, we summarize the most up-to-date evidence of m⁶A RNA modification in osteoporosis and highlight the potential role of m⁶A in prevention, treatment, and management of osteoporosis.