RNA Methylation in Systemic Lupus Erythematosus

Differential lncRNA profiles of blood plasma-derived exosomes from systemic lupus erythematosus

INTRODUCTION

Long non-coding RNAs (lncRNAs) dysregulation is key in the pathogenesis of systemic lupus erythematosus (SLE), but the role of exosomal lncRNAs in SLE has not been well studied. We elucidated the profiles of plasma exosomal lncRNAs expression in patients with SLE and predictd their potential clinical significance in SLE.

METHODS

In the screening stage, six newly diagnosed and untreated patients with SLE and six healthy controls were examined by high-throughput sequencing technology, and differential exosomal lncRNA profiles were constructed. In the validation phase, two differentially selected exosomal lncRNAs from 20 patients each with active and stable SLE and 20 healthy controls were verified with RT-qPCR. The correlation between the selected exosomal lncRNAs and SLE clinical indicators was examined. The diagnostic value of the selected exosomal lncRNAs in SLE was analyzed by the receiver operator characteristic (ROC) curve.

RESULTS

Exosomes were successfully extracted from the patients and controls. Sequencing-phase sequencing demonstrated 528 upregulated lncRNAs and 7491 downregulated lncRNAs. In the validation stage, exosomal LINC00667 and DANCR were significantly upregulated in the patients, and positively correlated with Systemic Lupus Erythematosus Disease Activity Index 2000 (SLEDAI-2 K). Exosomal DANCR expression between the active and stable SLE patients was different. The area under the curve(AUC) of exosomal LINC00667 and DANCR for SLE diagnosis was 0.815 and 0.759, respectively.

CONCLUSIONS

Exosomal LINC00667 and DANCR were upregulated in SLE, and might be new biomarkers thereof. Exosomal DANCR was associated with SLE activity.

NSUN2 mediates distinct pathways to regulate enterovirus 71 replication

Increasing evidences suggest that the methyltransferase NSUN2 catalyzes 5-methylcytosine (m5C) modifications on viral RNAs, which are essential for the replication of various viruses. Despite the function of m5C deposition is well characterized, other potential roles of NSUN2 in regulating viral replication remain largely unknown. In this study, the m5C modified residues catalyzed by NSUN2 on enterovirus 71 (EV71) RNAs were mapped. NSUN2, along with m5C modifications, played multiple roles during the EV71 life cycle. Functional m5C modified nucleotides increased the translational efficiency and stability of EV71 RNAs. Additionally, NSUN2 was found to target the viral protein VP1 for binding and promote its stability by inhibiting the ubiquitination. Furthermore, both viral replication and pathogenicity in mice were largely attenuated when functional m5C residues were mutated. Taken together, this study characterizes distinct pathways mediated by NSUN2 in regulating EV71 replication, and highlights the importance of its catalyzed m5C modifications on EV71 RNAs for the viral replication and pathogenicity.

RNA methylation patterns, immune characteristics, and autophagy-related mechanisms mediated by N6-methyladenosine (m6A) regulatory factors in venous thromboembolism

Recent studies have found a link between deep vein thrombosis and inflammatory reactions. N6-methyladenosine (m6A), a crucial element in immunological regulation, is believed to contribute to the pathophysiology of venous thromboembolism (VTE). However, how the m6A-modified immune microenvironment is involved in VTE remains unclear. In the present study, we identified a relationship between VTE and the expression of several m6A regulatory elements by analyzing peripheral blood samples from 177 patients with VTE and 88 healthy controls from public GEO databases GSE19151 and GSE48000. We used machine learning to identify essential genes and constructed a diagnostic model for VTE using multivariate logistic regression. Unsupervised cluster analysis revealed a marked difference between m6A modification patterns in terms of immune cell infiltration, inflammatory reactivity, and autophagy. We identified two m6A-related autophagy genes (i.e., CHMP2B and SIRT1) and the crucial m6A regulator YTHDF3 using bioinformatics. We also examined two potential mechanisms through which YTHDF3 may affect VTE. m6A modification, immunity, and autophagy are closely linked in VTE, offering novel mechanistic and therapeutic insights.

Reduced mitochondrial-encoded NADH dehydrogenase 6 gene expression drives inflammatory CD4+T cells in patients with systemic lupus erythematosus

Abnormal mitochondrial function has been implicated in the progression of systemic lupus erythematosus (SLE), the prototypical autoimmune disease, yet the underlying cause remains unclear. In this study, mitochondrial-encoded NADH dehydrogenase 6 gene (MT-ND6) was identified as having increased m6A methylation and decreased expression in peripheral blood mononuclear cells of SLE patients by MeRIP-seq analysis. MT-ND6 expression was negatively correlated with SLE disease activity index score and 24-h urine protein level, and lower in patients with positive anti-Sm or anti-dsDNA antibodies. With the reduction of MT-ND6 levels, CD4+ T cells in SLE patients exhibited mitochondrial dysfunction, as evidenced by increased levels of reactive oxygen species (ROS) and mitochondrial ROS and insufficient ATP production. Accordingly, in vitro MT-ND6 silencing induced abnormalities in the above mitochondrial indicators in CD4+ T cells, and promoted the development of both transcription and inflammatory factors in these cells. In contrast, treatment with targeted mitochondrial antioxidants largely counteracted the silencing effect of MT-MD6. Thus, reduced MT-ND6 in SLE patients may lead to mitochondrial dysfunction through ROS overproduction, thereby promoting inflammatory CD4+ T cells.

Characteristics of n6-methyladenosine (m6A) regulators and role of FTO/TNC in scleroderma

BACKGROUND

m6A regulators have important roles in a variety of autoimmune diseases, but their potential function in scleroderma, a refractory connective tissue disease, remains unclear. Tenascin C (TNC) is known to be a factor promoting collagen deposition in the development of scleroderma, but the regulatory relationship between TNC and m6A regulators is unknown.

METHODS

We extracted GSE33463 data consisting of forty-one healthy controls and sixty-one patients with scleroderma, and we analyzed the expression levels of twenty-one m6A regulators as well as the associations between them. In addition, we obtained random forest (RF) and nomogram models to predict the likehood of scleroderma. Next, we categorized the m6Aclusters and geneclusters by consensus clustering, and we performed an immune cell infiltration analysis for each cluster. Finally, we injected adenoviruses into a bleomycin (BLM)-induced mouse model of scleroderma, which was used to overexpress FTO and TNC. We assess the extent of skin fibrosis in the mice samples using pathology stains and measuring their hydroxyproline content and collagen mRNA.

RESULTS

We initially identified fourteen differentially expressed m6A regulators (WTAP, RBM15, CBLL1, FTO, ALKBH5, YTHDC1, YTHDC2, YTHDF1, YTHDF2, YTHDF3, RBMX, HNRNPC, IGFBP1 and IGFBP2). We found ALKBH5 to be positively associated with CBLL1 and RBM15, and FTO to be negatively associated with WTAP. In addition, we identified four m6A regulators (CBLL1, IGFBP1, YTHDF2 and IGFBP2) using a RF model, and we designed a nomogram model with those variables that proved reliable according to the calibration curve and clinical impact curve. We found that the m6Acluster A was correlated with Type 1 T helper cell infiltration and the genecluster A was correlated with regulatory T cell infiltration. Finally, we showed that FTO overexpression downregulated the m6A and mRNA levels of TNC, and alleviated skin fibrosis in the mouse model of scleroderma. Thus, our overexpression experiments provide preliminary evidence suggesting that TNC is an adverse factor in scleroderma.

CONCLUSION

Our approach might be useful as a new and accurate scleroderma diagnosis method. Moreover, our results suggested that FTO/TNC might be a novel scleroderma therapeutic target.

Identification of m6A-Related Biomarkers in Systemic Lupus Erythematosus: A Bioinformation-Based Analysis

Background

Systemic Lupus Erythematosus (SLE), a prototypical autoimmune disorder, presents a challenge due to the absence of reliable biomarkers for discerning organ-specific damage within SLE. A growing body of evidence underscores the pivotal involvement of N6-methyladenosine (m6A) in the etiology of autoimmune conditions.

Methods

The datasets, which primarily encompassed the expression profiles of m6A regulatory genes, were retrieved from the Gene Expression Omnibus (GEO) repository. The optimal model, selected from either Random Forest (RF) or Support Vector Machine (SVM), was employed for the development of a predictive nomogram model. To identify pivotal genes associated with SLE, a comprehensive screening process was conducted utilizing LASSO, SVM-RFE, and RF techniques. Within the realm of SLE susceptibility, Weighted Gene Co-expression Network Analysis (WGCNA) was harnessed to delineate relevant modules and hub genes. Additionally, MeRIP-qPCR assays were performed to elucidate key genes correlated with m6A targets. Furthermore, a Mendelian randomization study was conducted based on genome-wide association studies to assess the causative influence of MMP9 on ischemic stroke (IS), which is not only a severe cerebrovascular event but also a common complication of SLE.

Results

Twelve m6A regulatory genes was identified, demonstrating statistical significance (p < 0.05) and utilized for constructing a nomogram model using the RF algorithm. EPSTI1, USP18, HP, and MMP9, as the hub genes, were identified. MMP9 uniquely correlates with m6A modification and was causally linked to an increased risk of IS, as indicated by our inverse variance weighting analysis showing an odds ratio of 1.0134 (95% CI=1.0004-1.0266, p = 0.0440).

Conclusion

Our study identified twelve m6A regulators, shedding light on the molecular mechanisms underlying SLE risk genes. Importantly, our analysis established a causal relationship between MMP9, a key m6A-related gene, and ischemic stroke, a common complication of SLE, thereby providing critical insights for presymptomatic diagnostic approaches.

Diagnostic Values of METTL1-Related Genes and Immune Characteristics in Systemic Lupus Erythematosus

Purpose

Methyltransferase like 1 (METTL1) regulates epitranscriptomes via the m7G modification in mammalian mRNA and microRNA. Systemic lupus erythematosus (SLE) is caused by abnormal immune reactivity and has diverse clinical manifestations. RNA methylation as a mechanism to regulate gene expression is widely implicated in immune regulation. However, the role of m7G in immune response of SLE has not been extensively studied.

Patients and Methods

Expression of METTL1 was identified in the public dataset GSE122459 and validated in an independent cohort of SLE patients. We investigated the association between METTL1-expression and clinical manifestations of SLE. Subsequently, differentially expressed genes (DEG) that were correlated with METTL1-expression in GSE122459 were used for functional enrichment analysis. The correlation between infiltrating immune cells and METTL1, as well as candidate biomarkers identified to be correlated with either METTL1 or immune cell infiltration were assessed by single-sample GSEA. Potential mechanisms were explored with Gene ontology and KEGG pathway enrichment. Diagnostic performances of candidate biomarkers in SLE were analyzed.

Results

The mRNA and protein expression of METTL1 in SLE patients were significantly decreased in both datasets. METTL1-coexpressed DEGs were enriched in several key immune-related pathways. Activated CD8 T cells, activated CD4 T cells, memory B cells and type 2 helper T cells were different between patients with high and low METTL1 expression. Further, activated CD8 T-cells, activated CD4 T-cells, memory B-cells were correlated with METTL1. The genes of LAMP3, CD83, PDCD1LG2, IGKVD3D-20, IGKV5-2, IGKV2D-30, IGLV3-19 and IGLV4-60 were identified as candidate targets that were correlated with immune cell proportion. Moreover, LAMP3, CD83, and PDCD1LG2 expression were of diagnostic value in SLE as indicated by ROC analysis.

Conclusion

Our findings suggested that METTL1 and its candidate targets LAMP3, CD83, PDCD1LG2 may be used for diagnosing SLE and could be explored for developing targeted molecular therapy for SLE.

Drug Repositioning Identifies Six Drug Candidates for Systemic Autoimmune Diseases by Integrative Analyses of Transcriptomes from Scleroderma, Systemic Lupus Erythematosus, and Sjogren's Syndrome

The mechanisms of systemic autoimmune diseases (ADs) are still not clearly understood. Understanding the etiology of systemic ADs and identifying new therapeutic targets require a systems science approach. Using publicly available transcriptome data and bioinformatic analysis, we investigated the differential gene expression profiles of patients with scleroderma, systemic lupus erythematosus, and Sjogren's syndrome. Of these common differentially expressed gene signatures, 208 were regulated in the same direction (either upregulated or downregulated in all datasets) and used for drug repositioning. Six small molecule drug candidates (KU-0063794, YM-155 [sepantronium bromide], MST-312 [telomerase inhibitor IX], PLX-4720, ZM 336372, and 528116.cdx [PIK-75]) were discovered by drug repositioning as potential therapeutics for systemic ADs. The Search Tool for Chemical Interactions was used to find the anticipated target genes of the repositioned molecules. The PI3K/AKT pathway topped the list of common enriched pathways with the most anticipated target genes of the six repositioned small molecules. We also report here the molecular docking findings on the binding affinity between the repositioned drug candidates and genes from the protein-protein interaction network modules of anticipated target genes. Taken together, this study provides new insights and opens up new possibilities on both pathogenesis and treatment of systemic ADs through drug repositioning.

Reduction of Methyltransferase-like 3-Mediated RNA N6-Methyladenosine Exacerbates the Development of Psoriasis Vulgaris in Imiquimod-Induced Psoriasis-like Mouse Model

N6-methyladenosine (m⁶A) methylation is the most pervasive and intensively studied mRNA modification, which regulates gene expression in different physiological processes, such as cell proliferation, differentiation, and inflammation. Studies of aberrant m⁶A in human diseases such as cancer, obesity, infertility, neuronal disorders, immune diseases, and inflammation are rapidly evolving. However, the regulatory mechanism and physiological significance of m⁶A methylation in psoriasis vulgaris are still poorly understood. In this study, we found that m⁶A methylation and Methyltransferase-like 3 (METTL3) were both downregulated in psoriatic skin lesions and were negatively correlated with Psoriasis Area and Severity Index (PASI) scores. Inhibiting m⁶A methylation by knocking down Mettl3 promoted the development of psoriasis and increased its severity in imiquimod-induced psoriasis-like model mice. Our results indicate a critical role of METTL3- mediated m⁶A methylation in the pathogenesis of psoriasis vulgaris.

N6-methyladenine RNA Methylation Epigenetic Modification and Kidney Diseases

RNA methylation modification is a rapidly developing field in epigenetics. N6-methyladensine (m⁶A) is the most common internal modification in eukaryotic mRNA. m⁶A group regulates RNA splicing, stability, translocation, and translation. Enzymes catalyzing this process were termed as writers, erasers, and readers. Recent studies have focused on exploring the role of RNA methylation in human diseases. RNA methylation modifications, particularly m⁶A, play important roles in the pathogenesis of kidney diseases. In this review, we provide a brief description of m⁶A and summarize the impact of m⁶A on acute and chronic kidney disease (CKD) and possible future study directions for this research.

m5C Regulator-mediated methylation modification clusters contribute to the immune microenvironment regulation of multiple myeloma

Background: Multiple myeloma (MM) is a hematological malignancy in which plasma cells proliferate abnormally. 5-methylcytosine (m5C) methylation modification is the primary epigenetic modification and is involved in regulating the occurrence, development, invasion, and metastasis of various tumors; however, its immunological functions have not been systematically described in MM. Thus, this study aimed to clarify the significance of m5C modifications and how the immune microenvironment is linked to m5C methylation in MM.

Method: A total of 483 samples (60 healthy samples, 423 MM samples) from the Gene Expression Omnibus dataset were acquired to assess the expression of m5C regulators. A nomogram model was established to predict the occurrence of MM. We investigated the impact of m5C modification on immune microenvironment characteristics, such as the infiltration of immunocytes and immune response reactions. We then systematically evaluated three different m5C expression patterns to assess immune characteristics and metabolic functional pathways and established m5C-related differentially expressed genes (DEGs). In addition, biological process analysis was performed and an m5C score was constructed to identify potentially significant immunological functions in MM.

Result: Differential expressions of m5C regulators were identified between healthy and MM samples. The nomogram revealed that m5C regulators could predict higher disease occurrence of MM. We identified three distinct m5C clusters with unique immunological and metabolic characteristics. Among the three different m5C clusters, cluster C had more immune characteristics and more metabolism-related pathways than clusters A and B. We analyzed 256 m5C-related DEGs and classified the samples into three different m5C gene clusters. Based on the m5C and m5C gene clusters, we calculated m5C scores and classified each patient into high- and low-m5C score groups.

Conclusion: Our study demonstrated that m5C modification is involved in and contributes to the diversity and complexity of the immune microenvironment, which offers promise for the development of accurate therapeutic strategies.

METTL3 inhibits inflammation of retinal pigment epithelium cells by regulating NR2F1 in an m6A-dependent manner

N6-metyladenosine (m6A) RNA methylation has been proven to be involved in diverse biological processes, but its potential roles in the development of lipopolysaccharide (LPS) induced retinal pigment epithelium (RPE) inflammation have not been revealed. In this study, we explored the effects and underlying mechanisms of methyltransferase-like 3 (METTL3) in LPS stimulated RPE cells. Proliferation of METTL3-silenced RPE cells was examined by Cell counting kit-8 (CCK8) and 5-Ethynyl-2´-Deoxyuridine (Edu). Expression of tight junction proteins ZO-1 and Occludin, and secretion of inflammatory factors interleukins (IL)-1, 6 and 8 were detected by Western blotting or Enzyme-linked immunosorbent assay (ELISA). RNA sequencing and methylated RNA immunoprecipitation (MeRIP) sequencing were used to analyze the target gene nuclear receptor subfamily 2 group F member 1 (NR2F1) of METTL3. Our results showed that both human RPE (hRPE) cells and ARPE19 cells exhibited inhibited proliferation, tight junction protein expression, and increased inflammatory factor secretion after METTL3 silencing. Mechanistically, we found that NR2F1, as a METTL3-methylated target gene, inhibits Occludin level and promotes IL-6 secretion of RPE cells in an m6A-dependent manner. Interestingly, NR2F1 deficiency reversed the decreased Occludin expression and increased IL-6 secretion in METTL3-defective RPE cells. In conclusion, our study revealed that METTL3 attenuates RPE cell inflammation by methylating NR2F1, suggesting the critical role of METTL3 in RPE cells.

N6-methyladenosine-dependent modification of circGARS acts as a new player that promotes SLE progression through the NF-κB/A20 axis

Background

Certain circRNAs could be used as biomarkers to determine the risk of development and/or severity of systemic lupus erythematosus, and their new function in the regulation of gene expression has motivated us to investigate their role in SLE

Methods

Experimental methods including qRT-PCR, RNA immunoprecipitation (RIP), pulldown, dual luciferase reporter assay, RNA interference and cell transfection, RNA fluorescence in situ hybridization, western blotting, and mass spectrometry were used to assessed circGARS (hsa_circRNA_0009000) for immune functions and defined mechanisms by which circGARS promotes the progression in SLE.

Results

Our results demonstrated that the levels of circGARS was remarkably upregulated in SLE and correlated with clinicopathological features. CircGARS directly combined with microRNA-19a (miR-19a). Functionally, circGARS downregulated the expression of TNFAIP3 (A20, tumor necrosis factor alpha-induced protein 3) to mediate the activation of immune responses that were regulated by the nuclear factor-κB (NF-κB) pathway as a negative feedback mechanism. In addition, miR-19a regulated A20 (TNFAIP3) degradation by downregulating the expression of YTH N6-methyladenosine RNA-binding protein 2 (YTHDF2).

Conclusions

The circGARS sponges miR-19a to regulate YTHDF2 expression to promote SLE progression through the A20/NF-κB axis and may act as an independent biomarker to help the treatment of SLE patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13075-022-02732-x.

Involvement of N6-methyladenosine modifications of long noncoding RNAs in systemic lupus erythematosus

BACKGROUND

LncRNAs are potential biomarkers for SLE, but the epigenetic regulatory mechanisms of N6-methyladenosine (m6A) modification in SLE remain largely unclear.

METHODS

In this study, we established m6A modification profile and investigated the potential roles of m6A-related lncRNAs in SLE. The m6A modification profile of SLE was established using MeRIP-seq. Four potential m6A related-lncRNAs (linc02446, linc01410, Xist, and PSMB8-AS1) were selected for validation using qRT-PCR, and their expression and association with clinical characteristics with SLE were evaluated.

RESULTS

Overall, m6A level was lower in patients with SLE than in controls. Compared with controls, the expression of the two m6A related-lncRNAs (Xist and PSMB8-AS1) was downregulated in patients with SLE (all P < 0.05); the linc02446 was up-regulated in PBMCs of patients with SLE (Z=-2.738, P = 0.006), while it was not differentially expressed in T cells (Z=-0.387, P = 0.699). No significant alteration in linc01410 expression was observed in patients (Z=-0.940, P = 0.347). The lower expression levels of Xist and PSMB8-AS1 were associated with many clinical manifestations in patients with SLE (all P < 0.05). Additionally, mRNAs co-expressed with m6A related-lncRNAs (Xist, linc02446, and PSMB8-AS1) also participated in SLE.

CONCLUSION

These results suggest that m6A methylation and m6A related-lncRNAs might be involved in the pathogenesis of SLE. Thus, our findings provide some clues on the potential function of lncRNAs that m6A modification may target in novel therapeutic or diagnostic strategies for SLE.