Disease Markers
Volume 2020 (2020), Article ID 8193895, 11 pages
https://doi.org/10.1155/2020/8193895

Research Article

Decreased Peripheral Blood ALKBH5 Correlates with Markers of Autoimmune Response in Systemic Lupus Erythematosus

Qing Luo,¹ Biqi Fu,² Lu Zhang,¹ Yang Guo,¹ Zikun Huang,¹ and Junming Li¹

¹Department of Clinical Laboratory, The First Affiliated Hospital of Nanchang University, Nanchang, Jiangxi 330006, China
²Department of Rheumatology, The First Affiliated Hospital of Nanchang University, Nanchang, Jiangxi, China

Correspondence should be addressed to Zikun Huang; 491353062@qq.com and Junming Li; lisir361@163.com

Received 22 November 2019; Revised 2 June 2020; Accepted 10 June 2020; Published 25 June 2020

Academic Editor: Rudy Foddis

Copyright © 2020 Qing Luo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Although it has been proved that the epigenetic modification of DNA and histones is involved in the pathogenesis of systemic lupus erythematosus (SLE), there is no study to explore whether the modification of N6-methyladenosine (m6A) in RNA is involved. In this study, the mRNA levels of m6A “writers” (METTL3, MTEEL14, and WTAP), “erasers” (FTO and ALKBH5), and “readers” (YTHDF2) in peripheral blood were determined by quantitative reverse transcription-polymerase chain reaction (qRT-PCR). The results demonstrated that the mRNA levels of METTL3, WTAP, FTO, ALKBH5, and YTHDF2 in peripheral blood from SLE patients were significantly decreased. The levels of ALKBH5 mRNA in SLE patients were associated with anti-dsDNA, antinucleosome, rash, and ulceration. Multivariate logistic regression analysis showed that the level of ALKBH5 mRNA in peripheral blood is a risk factor of SLE (). Moreover, our results suggested that there was a positive correlation between m6A“writers” (METTL3 and WTAP), “erasers” (FTO and ALKBH5), and “readers” (YTHDF2) in SLE patients. This study suggests that the mRNA level of ALKBH5 in peripheral blood may be involved in the pathogenesis of SLE.

1. Introduction

The systemic lupus erythematosus (SLE) is a chronic and incurable autoimmune disease characterized by intermittent episodes of increased disease activity that require treatment with immunosuppressive agents [1, 2]. Although there have been many studies trying to elucidate the pathogenesis of SLE, so far it has not been fully elucidated. Existing studies have demonstrated that the dysfunction of immune cells such as T cells, B cells, monocytes, neutrophils, and dendritic cells plays important roles in the pathogenesis of SLE [3–6]. Further elucidation of the aetiology of SLE is of great significance for the development of possible targeted and individualized therapy for SLE [7].

In recent years, epigenetic modifications have been demonstrated to play an important role in the genesis and development of SLE [8, 9]. N6-methyladenosine (m6A) modification is the most prevalent and evolutionarily conserved modification which occurs in nearly all types of RNAs and in most organisms [10]. This modification can be installed by adenosine methyltransferases, reversed by demethylases, and recognized by some RNA-binding proteins [11]. M6A methyltransferase complex, known as the m6A “writer” complex, contains methyltransferase-like 3 (METTL3), methyltransferase-like 14 (METTL14), and Wilms tumor 1-associating protein (WTAP), which functions by depositing the m6As in mammalian mRNA [12]. Fat mass and obesity-associated protein (FTO) and a-ketoglutarate-dependent dioxygenase alkB homolog 5 (ALKBH5) are selective demethylases capable of regulating gene expression and cell fate through oxidative removal of the methyl group in m6A-containing substrates, usually considered as m6A “erasers” [13]. Some RNA-binding proteins such as YT521-B homology domains 2 (YTHDF2) which can recognize m6A modification, decode the methylation code, and finally transform them into diverse functional signals are called m6A “readers” [14].

Recent studies have demonstrated that m6A modification is associated with various human diseases [15, 16]. However, there is no study to characterize m6A modification in patients with SLE. To investigate whether m6A modification plays a role in the genesis and development of SLE, the mRNA levels of METTL3, MTEEL14, WTAP, FTO, ALKBH5, and YTHDF2 in peripheral blood were detected in SLE patients and analyzed for their correlation with clinical variables.

2. Methods

2.1. Patient Variables and Controls

A total of 51 patients that fulfilled the revised American College of Rheumatology criteria for SLE [17] were recruited from the First Affiliated Hospital of Nanchang University from 2018.10 to 2019.3. Among them, 40 patients were new-onset SLE that first-time diagnosis of SLE and no history of immunosuppressive drug or corticosteroid use before recruitment. Among all new-onset SLE patients, 7 patients were reexamined after 15 days of regular treatment by using glucocorticoids and immunosuppressive agents. The other 11 patients were revisiting SLE patients receiving treatment. Disease activity was assessed by the SLE disease activity index (SLEDAI) [18]. 38 healthy controls (CON) without a clinical diagnosis of any inflammatory or autoimmune diseases and without relation to patients of autoimmune disease were enrolled from the First Affiliated Hospital of Nanchang University. In addition, 51 patients fulfilled the revised ACR 2010 criteria for rheumatoid arthritis (RA) [19], 30 patients were infected with hepatitis B virus (HBV) and 27 patients with tuberculosis (TB) were recruited from the First Affiliated Hospital of Nanchang University. The demographic characteristics of the study population are shown in Table 1. The study had approval from the Ethics Committee of the First Affiliated Hospital of Nanchang University (052) and complied with the Helsinki Declaration. All participants provided signed informed consent before they entered this study.

Table 1: Clinical characteristics of SLE patients, RA patients, HBV-infected patients, TB patients, and CON.

Categories
SLE patients
CON
RA patients
HBV
TB patients

51
38
51
30
27

Females (%)
92.16
78.95
82.35
93.33
100

Age (years)

SLEDAI score

Anti-dsDNA (IU/mL)

Anti-ENA

Anti-Sm (%)
29.73

Anti-Ro52 (%)
70.27

Anti-nRNP/Sm (%)
59.46

Anti-RIB-P (%)
24.32

Antinucleosome (%)
32.43

Anti-SSA (%)
64.86

Anti-SSB (%)
18.92

C3 (g/L)

C4 (g/L)

IgG (g/L)

ESR (mm/h)

CRP (mg/L)

WBC (10⁹/L)

RBC (10¹²/L)

HGB (g/L)

HCT (L/L)

PLT (10⁹/L)

Lymphocytes (10⁹/L)

Lymphocytes (%)

Monocytes (10⁹/L)

Monocytes (%)

Neutrophils (10⁹/L)

Neutrophils (%)

Clinical features

Fever (%)
44.44

Rash (%)
36.11

Alopecia (%)
17.14

Arthritis (%)
47.22

NPLE (%)
11.11

Ulceration (%)
11.43

Pleuritis (%)
45.71

Pericarditis (%)
34.29

LN (%)
54.90

SLE compared to CON. Anti-dsDNA: anti-double-stranded DNA; Anti-ENA: antiextractable nuclear antigen; Anti-nRNP/Sm: antinuclear ribonucleoprotein/Smith antibody; Anti-RIB-P: anti-ribosomal P-protein antibody; Anti-Sm: anti-Smith antibody; Anti-SSA: anti-Sjögren syndrome A antigen antibody; Anti-SS-B: anti-Sjögren syndrome B antigen antibody; HBV: hepatitis B virus (HBV); C3: complement 3; C4: complement 4; HC: healthy controls; CRP: C-reactive protein; ESR: erythrocyte sedimentation rate; HCT: hematocrit; HGB: hemoglobin; LN: lupus nephritis; IgG: immunoglobulin G; L: lymphocyte count; L%: lymphocyte percentage; M: monocyte count; M%: monocyte percentage; N: neutrophil count; N%: neutrophil percentage; NPLE: neuropathic lupus erythematosus; PLT: platelet count; RBC: red blood cell count; RA: rheumatoid arthritis; SLE: systemic lupus erythematosus; SLEDAI: SLE disease activity index; TB: tuberculosis; WBC: white blood cell count.

2.2. Blood Sample Collection and Total RNA Isolation

The peripheral blood collection and total RNA isolation were performed as described previously [20].

2.3. QRT-PCR Analysis

Complementary DNA (cDNA) was obtained from purified total RNA samples by reverse transcription using a PrimeScript™ RT reagent kit (Takara Bio Inc., Japan). METTL3, METTL14, WTAP, ALKBH5, FTO, and YTHDF2 transcripts were quantified with an ABI 7500 Real-time PCR System (Applied Biosystems; Thermo Fisher Scientific, Inc.) by using SYBR® Premix Ex Taq™ II (Takara Bio Inc., Japan). The sequences of amplification primers for METTL3, METTL14, WTAP, ALKBH5, FTO, YTHDF2, and GAPDH are listed in Table 2. The relative levels of transcripts were derived by the 2-ΔCt method [20, 21].

Table 2: The amplification primers sequences.

Gene name
Sequence (5-3)

METTL3
F: AAGCTGCACTTCAGACGAAT

R: GGAATCACCTCCGACACTC

METTL14
F: AGAAACTTGCAGGGCTTCCT

R: TCTTCTTCATATGGCAAATTTTCTT

WTAP
F: GGCGAAGTGTCGAATGCT

R: CCAACTGCTGGCGTGTCT

ALKBH5
F: CCCGAGGGCTTCGTCAACA

R: CGACACCCGAATAGGCTTGA

FTO
F: TGGGTTCATCCTACAACGG

R: CCTCTTCAGGGCCTTCAC

YTHDF2
F: GGCAGCACTGAAGTTGGG

R: CTATTGGAAGCCACGATGTTA

GAPDH
F: TGCACCACCAACTGCTTAGC

R: GGCATGGACTGTGGTCATGAG

METTL3: methyltransferase-like 3; METTL14: methyltransferase-like 14; WTAP: Wilms tumor 1-associating protein; ALKBH5: a-ketoglutarate-dependent dioxygenase alkB homolog 5; FTO: fat mass and obesity-associated protein; YTHDF2: YT521-B homology domains 2.

2.4. Blood Routine, Urinary Routine, Serum Inflammatory Indicators, and Autoantibody Determination

The items in blood routine and urinary routine were detected using Sysmex-x2100 (SYSMEX, Kobe, Japan) and SysmexUF-1000i (SYSMEX, Japan). Erythrocyte sedimentation rate (ESR) was determined by the BlessLBY-XC40B automatic ESR analyzer (Bless, Beijing, China) according to the manufacturer’s instructions. The levels of C-reactive protein (CRP), immunoglobulin G (IgG), complement 3 (C3), and complement 4 (C4) in serum were determined by nephelometry methods using IMMAGE800 system (Beckman, CA, USA). Anti-dsDNA of IgG in serum was measured by commercially available enzyme-linked immunosorbent assay (ELISA) kits (Kexin, Shanghai, China). Antiextractable nuclear antigen (Anti-ENA) was determined using line immunoassay kits (Euroimmun, Luebeck, Germany) according to the manufacturer’s instructions.

2.5. Statistical Analysis

Student’s -test or the Mann-Whitney -test was used to compare the data according to the normality. The Spearman method was used for correlation analysis. Logistic regression analysis was used for evaluating the risk factor. All the data were analyzed by GraphPad Prism version 5.0 (GraphPad Software, San Diego, CA) and SPSS version 17.0 (SPSS Inc., Chicago, IL). values < 0.05 were considered statistically significant.

3. Results

3.1. Differential Expression Screening of METTL3, METTL14, WTAP, ALKBH5, FTO, and YTHDF2 in Peripheral Blood of SLE Patients

To identify the expression levels of METTL3, METTL14, WTAP, ALKBH5, FTO, and YTHDF2 in patients with SLE, the mRNA levels of these genes in peripheral blood were determined by qRT-PCR in SLE patients and HCs. Data showed that the mRNA levels of METTL3, WTAP, ALKBH5, and FTO in peripheral blood of SLE patients significantly decreased compared to HCs (all ) (Figure 1). The mRNA level of YTHDF2 trends to decrease in SLE patients, but a significant difference was not reached () (Figure 1). The mRNA level of METTL14 was unchanged () (Figure 1).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 1: Differential expression screening of METTL3, METTL14, WTAP, ALKBH5, FTO, and YTHDF2 in peripheral blood of SLE patients. (a) The level of METTL3 in peripheral blood was significantly decreased in SLE patients () compared to healthy controls (CON, ). (b) The level of METTL14 showed no significant difference between SLE patients () and CON (). (c) The level of WTAP in peripheral blood was significantly decreased in SLE patients () compared to CON (). (d) The level of ALKBH5 in peripheral blood was significantly decreased in SLE patients () compared to CON (). (e) The level of FTO in peripheral blood was significantly decreased in SLE patients () compared to CON (). (f) The level of YTHDF2 trends to decline in SLE patients, but a significant difference was not reached.

3.2. Differential Expression Validation of METTL3, WTAP, ALKBH5, FTO, and YTHDF2

To further verify the differential expression of METTL3, WTAP, ALKBH5, FTO, and YTHDF2 in peripheral blood in the screening stage, an independent cohort including 23 SLE patients and 12 HCs were recruited and detected for the mRNA levels of these genes in peripheral blood. The data from all subjects, including all 51 SLE patients and 38 HCs in screening stage and validation stage, demonstrated that the mRNA levels of METTL3, WTAP, ALKBH5, and FTO in peripheral blood were significantly lower in SLE patients than those in HCs (all ) (Figure 2). Inconsistent with the screening data, the level of YTHDF2 was also significantly decreased in SLE patients than that in HCs () (Figure 2).

(a)

(b)

(c)

(d)

(e)

Figure 2: Differential expression validation of peripheral blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2. (a) The level of METTL3 in peripheral blood was significantly decreased in SLE patients () compared to healthy controls (CON, ). (b) The level of WTAP in peripheral blood was significantly decreased in SLE patients () compared to CON (). (c) The level of ALKBH5 in peripheral blood was significantly decreased in SLE patients () compared to CON (). (d) The level of FTO in peripheral blood was significantly decreased in SLE patients () compared to CON (). (e) The level of YTHDF2 in peripheral blood was significantly decreased in SLE patients () compared to CON ().

3.3. Correlations between the mRNA Levels of Peripheral Blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2 and Autoantibodies in SLE Patients

Production of multiple autoantibodies such as anti-double-stranded DNA (anti-dsDNA) and antiextractable nuclear antigens (Anti-ENAs) is one of the important characteristics of SLE. In 51 SLE patients, 44 SLE patients were tested for anti-dsDNA and 37 SLE patients were tested for anti-ENAs. The correlations between the mRNA levels of peripheral blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2 and autoantibodies were investigated in these SLE patients. As shown in Figure 3(a), the level of peripheral blood ALKBH5 was significantly decreased in patients with positive anti-dsDNA compared to patients with negative anti-dsDNA. Furthermore, the level of peripheral blood ALKBH5 in SLE patients negatively correlated with the level of anti-dsDNA (, ) (Figure 3(b)). And, the level of peripheral blood ALKBH5 was significantly decreased in patients with positive antinucleosome compared to patients with negative antinucleosome () (Figure 3(c)). The level of peripheral blood ALKBH5 did not correlate with other anti-ENAs (data not shown). No obvious correlation was observed between the levels of peripheral blood METTL3, WTAP, FTO, and YTHDF2 and autoantibodies in SLE patients (data not shown).

(a)

(b)

(c)

Figure 3: Correlations between the mRNA level of peripheral blood ALKBH5 and autoantibodies in SLE patients. (a) The mRNA level of peripheral blood ALKBH5 was significantly decreased in SLE patients with positive anti-double-stranded DNA (anti-dsDNA) compared to SLE patients with negative anti-dsDNA. (b) The mRNA level of peripheral blood ALKBH5 in SLE patients negatively correlated with the level of anti-dsDNA. (c) The level of peripheral blood ALKBH5 was significantly decreased in SLE patients with positive antinucleosome compared to SLE patients with negative antinucleosome.

3.4. Decreased mRNA Level of ALKBH5 in Peripheral Blood Was a Risk Factor for SLE

The aforementioned results demonstrated that the levels of peripheral blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2 in SLE patients were decreased and the level of ALKBH5 in peripheral blood was associated with autoantibody production. Considering that the production of autoantibodies plays an important role in the pathogenesis of SLE, the “enter method” of multivariate logistic regression was used to explore whether the expression of peripheral blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2 was risk factors for SLE. As shown in Table 3, the equations about the levels of peripheral blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2 were obtained, . The results demonstrated that the decreased expression of ALKBH5 in peripheral blood was a risk factor for SLE (), while other molecules are not (). Furthermore, the mRNA level of peripheral blood ALKBH5 in other inflammatory disorders including RA, TB, and HBV-infected patients was detected and compared with that in SLE patients. Results showed that the mRNA level of peripheral blood ALKBH5 in SLE patients was significantly lower compared to RA patients () (Figure 4(a)), HBV-infected patients () (Figure 4(b)), and TB patients () (Figure 4(c)).

Table 3: The expression of METTL3, WTAP, ALKBH5, FTO, and YTHDF2 in equation.

B
S.E
Wald
df

Exp (B)

METTL3
-326.128
272.636
1.431
1
0.232
0.000

WTAP
-74.188
66.854
1.231
1
0.267
0.000

ALKBH5
-414.645
108.988
14.474
1
0.000
0.000

FTO
-15.042
58.027
0.067
1
0.795
0.000

YTHDF2
-0.106
14.726
0.000
1
0.994
0.899

Constant
6.197
1.408
19.373
1
0.000
491.119

METTL3: methyltransferase-like 3; WTAP: Wilms tumor 1-associating protein; ALKBH5: a-ketoglutarate-dependent dioxygenase alkB homolog 5; FTO: fat mass and obesity-associated protein; YTHDF2: YT521-B homology domains 2.

(a)

(b)

(c)

Figure 4: Validation of the abnormal expression of peripheral blood ALKBH5 in SLE patients, rheumatoid arthritis (RA) patients, hepatitis B virus- (HBV-) infected patients, and tuberculosis (TB) patients. (a) The peripheral blood level of ALKBH5 in SLE patients () was significantly lower compared to RA patients (). (b) The peripheral blood level of ALKBH5 in SLE patients () was significantly lower compared to HBV-infected patients (). (c) The peripheral blood level of ALKBH5 in SLE patients () was significantly lower compared to TB patients ().

3.5. Correlations between the Levels of Peripheral Blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2 and Clinical Variables in SLE Patients

The aforementioned results demonstrate that the level of peripheral blood ALKBH5 in SLE patients was negatively correlated with the level of anti-dsDNA. Anti-dsDNA was proved to be involved in the genesis of SLE and can be served as markers for evaluating ongoing disease activity in SLE and early relapse of SLE. And, the level of anti-dsDNA in SLE patients was found to negatively correlate with the level of C3 (, ) (Supplement Fig 1A) and C4 (, ) (Supplement Fig 1B) and positively correlate with SLEDAI (, ) (Supplement Fig 1C). Thus, the correlation test was performed to evaluate the correlations between the clinical variables of SLE and the levels of peripheral blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2. Data showed that the mRNA levels of peripheral blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2 in all SLE patients (including all new-onset and revisiting SLE patients) did not correlate with SLEDAI, CRP, ESR, IgG, C3, C4, WBC, RBC, HGB, HCT, PLT, L, M, or N (data not shown). However, further analysis in the new-onset SLE patients showed that the mRNA level of peripheral blood METTL3 positively correlated with PLT (, ) (Figure 5(a)), the mRNA level of peripheral blood ALKBH5 positively correlated with WBC (, ) (Figure 5(b)), and the mRNA level of peripheral blood WTAP positively correlated with M (, ) (Figure 5(c)).

(a)

(b)

(c)

Figure 5: Correlations between peripheral blood METTL3, ALKBH5 and WTAP and clinical variables in new-onset SLE patients. (a) The level of peripheral blood METTL3 in new-onset SLE patients positively correlated with PLT (). (b) The level of peripheral blood ALKBH5 in new-onset SLE patients positively correlated with WBC (). (c) The level of peripheral blood WTAP in new-onset SLE patients positively correlated with M ().

Next, the relationship between the levels of these m6A modification-related molecules and drug therapy of SLE was analyzed. However, we found no significant difference between new-onset SLE patients and revisiting SLE patients (data not shown). Subsequently, a 15-day follow-up evaluation in 7 new-onset SLE patients who received regular treatment was performed, but there was still no difference (data not shown).

Moreover, the relationship between the clinical symptoms of SLE including LN, NPLE, arthritis, fever, rash, alopecia, ulceration, pleuritis, pericarditis, and the mRNA levels of these m6A modification-related molecules was analyzed. As shown in Figure 6, the mRNA levels of peripheral blood METTL3 and WTAP in SLE patients with alopecia was significantly decreased than that in SLE patients without alopecia (all ); the mRNA level of peripheral blood ALKBH5 in SLE patients with rash and ulceration was significantly increased than that in SLE patients without rash and ulceration, respectively (all ). No relationship was found between the levels of peripheral blood FTO, YTHDF2, and SLE symptoms (Figures 6(d) and 6(e)).

(a)

(b)

(c)

(d)

(e)

Figure 6: Correlations between peripheral blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2 and clinical symptoms of SLE. (a) The level of peripheral blood METTL3 in SLE patients with alopecia was significantly decreased than that in SLE patients without alopecia. (b) The level of peripheral blood WTAP in SLE patients with alopecia was significantly decreased than that in SLE patients without alopecia. (c) The levels of peripheral blood ALKBH5 in SLE patients with rash and ulceration were significantly increased than those in SLE patients without rash and ulceration, respectively. (d) No correlation was found between peripheral blood FTO and clinical symptoms of SLE. (e) No correlation was found between peripheral blood YTHDF2 and clinical symptoms of SLE.

3.6. Correlations between the Levels of These m6A Modification-Related Molecules in SLE Patients

Considering that METTL3, WTAP, ALKBH5, FTO, and YTHDF2 are all key molecules in the process of m6A modification, we therefore want to know whether they are interrelated. As shown in Figure 5, the mRNA level of peripheral blood METTL3 positively correlated with the mRNA levels of WTAP (, ) (Figure 7(a)), ALKBH5 (, ) (Figure 7(b)), FTO (, ) (Figure 7(c)), and YTHDF2 (, ) (Figure 7(d)); the mRNA level of peripheral blood WTAP positively correlated with the mRNA level of FTO (, ) (Figure 7(e)); and the mRNA level of peripheral blood YTHDF2 positively correlated with the mRNA levels of ALKBH5 (, ) (Figure 7(f)) and FTO (, ) (Figure 7(g)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 7: Interrelations between peripheral blood METTL3, WTAP, ALKBH5, FTO, and YTHDF2 in SLE patients. (a) The level of peripheral blood METTL3 positively correlated with WTAP. (b) The level of peripheral blood METTL3 positively correlated with ALKBH5. (c) The level of peripheral blood METTL3 positively correlated with FTO. (d) The level of peripheral blood METTL3 positively correlated with YTHDF2. (e) The level of peripheral blood WTAP positively correlated with FTO. (f) The level of peripheral blood YTHDF2 positively correlated with ALKBH5. (g) The level of peripheral blood YTHDF2 positively correlated with FTO.

4. Discussion

Epigenetic alterations, including DNA methylation and histone modifications, have been reported to contribute to SLE progression [8, 22, 23]. Zhao et al. showed that RFX1 regulates CD70 and CD11a expression in lupus T cells by recruiting the histone methyltransferase SUV39H1, triggering autoimmune responses [24]. In recent years, the fields of m6A modification and epitranscriptomics have attracted much attention [25]. As an important form of posttranscriptional gene regulation, m6A modification was found to be involved in almost all aspects of RNA metabolism, which is widely involved in the T cell response to HIV infection [26], type I interferon production [27], and T cell differentiation and homeostasis [28]. Thus, it invites the speculation that m6A modification may be involved in the genesis and progression of SLE [13].

As the regulators of m6A methylation, METTL3, MTEEL14, WTAP, FTO, ALKBH5, and YTHDF2 play pivotal roles in the dynamic regulation of m6A modification. Nettersheim et al. showed that the mRNA levels of METTL3, MTEEL14, WTAP, and ALKBH5 were significantly decreased in testicular germ cell tumors and the m6A modification in RNA of testicular germ cell tumors was controlled by METTL3, ALKBH5, and YTHDF2 [29]. Tian et al. indicated that the mRNA expressions of METTL3 and FTO during the osteogenic induction of BMSCs were increased [30]. In addition, Yang et al. demonstrated that the mRNA levels of FTO, METTL3, METTL14, and WTAP were increased in patients with type-2 diabetes [16]. These studies indicated m6A methylation regulators including m6A methylation writers, erasers, and readers which were changed synchronously. In this study, we also found that the mRNA levels of m6A methylation writers (METTL3 and WTAP), erasers (ALKBH5 and FTO), and reader (YTHDF2) were all decreased in peripheral blood from SLE patients. In accordance with the results of Yang et al. [16], our results manifested that the mRNA levels of peripheral blood m6A methylation writers positively correlated with those of erasers and reader, and the mRNA levels of peripheral blood m6A methylation erasers positively correlated with that of m6A methylation reader in SLE patients. Considering the fact that the collaboration among writers, erasers, and readers sets up the m6A threshold and perturbs that m6A threshold leads to uncontrolled expression/activity of virulence gene and results in the occurrence and development of diseases [31], our study suggested that the decreased METTL3, WTAP, ALKBH5, FTO, and YTHDF2 may play an important role in the pathogenesis of SLE.

SLE is a systemic autoimmune disease characterized by elevated autoimmune antibodies, such as anti-dsDNA and antinucleosome. It is well known that anti-dsDNA and antinucleosome are pathogenic autoantibodies involved in organ damage of SLE [32–34]. In this study, we found that the level of peripheral blood ALKBH5 was significantly decreased in SLE patients compared to HCs, RA patients, TB patients, and HBV-infected patients. Furthermore, the level of peripheral blood ALKBH5 significantly decreased in SLE patients with positive anti-dsDNA or antinucleosome compared to SLE patients with negative anti-dsDNA or antinucleosome. And, the level of peripheral blood ALKBH5 in SLE patients was found to be negatively correlated with the level of anti-dsDNA. In addition, the level of peripheral blood ALKBH5 in SLE patients was associated with clinical symptoms of SLE such as rash, ulceration, and leukopenia. Logistic regression analysis revealed that decreased mRNA level of ALKBH5 in peripheral blood was a risk factor for SLE.

Considering that anti-dsDNA and antinucleosome are involved in the pathogenesis of SLE and are specific diagnostic biomarkers of SLE, our results suggested that the expression of peripheral blood ALKBH5 may be involved in the pathogenesis of SLE and may also serve as a potential biomarker for SLE diagnosis and severity evaluation.

In conclusion, it is the first study to our knowledge that explored the mRNA levels of peripheral blood METTL3, METTL14, WTAP, ALKBH5, FTO, and YTHDF2 in SLE patients. In this study, we found that the mRNA levels of METTL3, WTAP, ALKBH5, FTO, and YTHDF2 in peripheral blood from SLE patients were significantly decreased. Additionally, this study established the correlations between the mRNA level of ALKBH5 and the production of autoantibodies, clinical features, which might improve our understanding of the role of ALKBH5 in SLE.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

QL and JML designed the study. BQF recruited the participants. JL and ZKH did the experimental work. LZ and YG collected the data and participated in the analysis of the data, supervised by QL which performed statistical analyses. QL, JML, JL, and ZKH drafted the manuscript, and all authors revised it critically for intellectual content. All authors have approved the final version. Qing Luo and Biqi Fu contributed equally to this study.

Acknowledgments

This work was supported by the Key Research and Development Plan Project of Jiangxi Province (20181BBG70013), the Science and Technology Plan Project of the Education Department of Jiangxi Province (GJJ170008), the National Natural Science Foundation of China (81360459, 81560344, and 81660277), the Jiangxi Provincial Natural Science Foundation of China (20171BAB205113 and 20171ACB20032), and the Foundation for Distinguished Young Scientists of Jiangxi Province of China (20171BCB23087).

Supplementary Materials

Supplement Figure 1: correlations between anti-double-stranded DNA (anti-dsDNA) and complement 3 (C3), complement 4 (C4), and systemic lupus erythematosus disease activity index (SLEDAI) in systemic lupus erythematosus (SLE) patients. (a) The level of anti-dsDNA negatively correlated with that of C3. (b) The level of anti-dsDNA negatively correlated with that of C4. (c) The level of anti-dsDNA positively correlated with SLEDAI. (Supplementary Materials)

References

M. Frieri, “Mechanisms of disease for the clinician: systemic lupus erythematosus,” Annals of Allergy, Asthma & Immunology, vol. 110, no. 4, pp. 228–232, 2013.

D. P. D'Cruz, M. A. Khamashta, and G. R. Hughes, “Systemic lupus erythematosus,” Lancet, vol. 369, no. 9561, pp. 587–596, 2007.

S. A. Jenks and I. Sanz, “Altered B cell receptor signaling in human systemic lupus erythematosus,” Autoimmunity Reviews, vol. 8, no. 3, pp. 209–213, 2009.

J. C. Crispin, V. C. Kyttaris, C. Terhorst, and G. C. Tsokos, “T cells as therapeutic targets in SLE,” Nature Reviews Rheumatology, vol. 6, no. 6, pp. 317–325, 2010.

Q. Luo, Z. Huang, J. Ye et al., “PD-L1-expressing neutrophils as a novel indicator to assess disease activity and severity of systemic lupus erythematosus,” Arthritis Research & Therapy, vol. 18, no. 1, p. 47, 2016.

G. C. Tsokos, “Systemic lupus erythematosus,” The New England Journal of Medicine, vol. 365, no. 22, pp. 2110–2121, 2011.

G. M. C. La Paglia, M. C. Leone, G. Lepri et al., “One year in review 2017: systemic lupus erythematosus,” Clinical and Experimental Rheumatology, vol. 35, no. 4, pp. 551–561, 2017.

H. Zhu, W. Mi, H. Luo et al., “Whole-genome transcription and DNA methylation analysis of peripheral blood mononuclear cells identified aberrant gene regulation pathways in systemic lupus erythematosus,” Arthritis Research & Therapy, vol. 18, no. 1, p. 162, 2016.

X. Wang, C. Zhang, Z. Wu, Y. Chen, and W. Shi, “CircIBTK inhibits DNA demethylation and activation of AKT signaling pathway via miR-29b in peripheral blood mononuclear cells in systemic lupus erythematosus,” Arthritis Research & Therapy, vol. 20, no. 1, p. 118, 2018.

X. Wang, Z. Lu, A. Gomez et al., “N⁶-methyladenosine-dependent regulation of messenger RNA stability,” Nature, vol. 505, no. 7481, pp. 117–120, 2014.

N. Pinello, S. Sun, and J. J. Wong, “Aberrant expression of enzymes regulating m6A mRNA methylation: implication in cancer,” Cancer Biology & Medicine, vol. 15, no. 4, pp. 323–334, 2018.

G. Tuncel and R. Kalkan, “Importance of m N⁶-methyladenosine (m⁶A) RNA modification in cancer,” Medical Oncology, vol. 36, no. 4, p. 36, 2019.

L. J. Li, Y. G. Fan, R. X. Leng, H. F. Pan, and D. Q. Ye, “Potential link between m⁶A modification and systemic lupus erythematosus,” Molecular Immunology, vol. 93, pp. 55–63, 2018.

A. Maity and B. Das, “N6-methyladenosine modification in mRNA: machinery, function and implications for health and diseases,” The FEBS Journal, vol. 283, no. 9, pp. 1607–1630, 2016.

S. Yang, J. Wei, Y. H. Cui et al., “m⁶A mRNA demethylase FTO regulates melanoma tumorigenicity and response to anti-PD-1 blockade,” Nature Communications, vol. 10, no. 1, p. 2782, 2019.

Y. Yang, F. Shen, W. Huang et al., “Glucose is involved in the dynamic regulation of m6A in patients with type 2 diabetes,” The Journal of Clinical Endocrinology and Metabolism, vol. 104, no. 3, pp. 665–673, 2019.

E. M. Tan, A. S. Cohen, J. F. Fries et al., “The 1982 revised criteria for the classification of systemic lupus erythematosus,” Arthritis and Rheumatism, vol. 25, no. 11, pp. 1271–1277, 1982.

C. Bombardier, D. D. Gladman, M. B. Urowitz et al., “Derivation of the SLEDAI. A disease activity index for lupus patients,” Arthritis and Rheumatism, vol. 35, no. 6, pp. 630–640, 1992.

D. Aletaha, T. Neogi, A. J. Silman et al., “2010 rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative,” Arthritis and Rheumatism, vol. 62, no. 9, pp. 2569–2581, 2010.

Q. Luo, L. Zhang, X. Li et al., “Identification of circular RNAs hsa_circ_0044235 in peripheral blood as novel biomarkers for rheumatoid arthritis,” Clinical and Experimental Immunology, vol. 194, no. 1, pp. 118–124, 2018.

F. Zhang, L. Wu, J. Qian et al., “Identification of the long noncoding RNA NEAT1 as a novel inflammatory regulator acting through MAPK pathway in human lupus,” Journal of Autoimmunity, vol. 75, pp. 96–104, 2016.

H. Wu, M. Zhao, L. Tan, and Q. Lu, “The key culprit in the pathogenesis of systemic lupus erythematosus: aberrant DNA methylation,” Autoimmunity Reviews, vol. 15, no. 7, pp. 684–689, 2016.

Y. Renaudineau and P. Youinou, “Epigenetics and autoimmunity, with special emphasis on methylation,” The Keio Journal of Medicine, vol. 60, no. 1, pp. 10–16, 2011.

M. Zhao, X. Wu, Q. Zhang et al., “RFX1 regulates CD70 and CD11a expression in lupus T cells by recruiting the histone methyltransferase SUV39H1,” Arthritis Research & Therapy, vol. 12, no. 6, p. R227, 2010.

Y. Yang, W. Huang, J. T. Huang et al., “Increased N⁶-methyladenosine in Human Sperm RNA as a Risk Factor for Asthenozoospermia,” Scientific Reports, vol. 6, no. 1, article 24345, 2016.

G. Lichinchi, S. Gao, Y. Saletore et al., “Dynamics of the human and viral m⁶A RNA methylomes during HIV-1 infection of T cells,” Nature Microbiology, vol. 1, no. 4, p. 16011, 2016.

Q. Zheng, J. Hou, Y. Zhou, Z. Li, and X. Cao, “The RNA helicase DDX46 inhibits innate immunity by entrapping m⁶A-demethylated antiviral transcripts in the nucleus,” Nature Immunology, vol. 18, no. 10, pp. 1094–1103, 2017.

H. B. Li, J. Tong, S. Zhu et al., “m⁶A mRNA methylation controls T cell homeostasis by targeting the IL-7/STAT5/SOCS pathways,” Nature, vol. 548, no. 7667, pp. 338–342, 2017.

D. Nettersheim, D. Berger, S. Jostes, G. Kristiansen, G. Lochnit, and H. Schorle, “N6-Methyladenosine detected inRNAof testicular germ cell tumors is controlled byMETTL3,ALKBH5,YTHDC1/F1/F2, andHNRNPCas writers, erasers, and readers,” Andrology, vol. 7, no. 4, pp. 498–506, 2019.

C. Tian, Y. Huang, Q. Li, Z. Feng, and Q. Xu, “Mettl3 regulates osteogenic differentiation and alternative splicing of Vegfa in bone marrow mesenchymal stem cells,” International Journal of Molecular Sciences, vol. 20, no. 3, p. 551, 2019.

S. Panneerdoss, V. K. Eedunuri, P. Yadav et al., “Cross-talk among writers, readers, and erasers of m6A regulates cancer growth and progression,” Science Advances, vol. 4, no. 10, article eaar8263, 2018.

S. Yung and T. M. Chan, “Mechanisms of kidney injury in lupus nephritis—the role of anti-dsDNA antibodies,” Frontiers in Immunology, vol. 6, p. 475, 2015.

J. A. Gómez-Puerta, R. W. Burlingame, and R. Cervera, “Anti-chromatin (anti-nucleosome) antibodies: diagnostic and clinical value,” Autoimmunity Reviews, vol. 7, no. 8, pp. 606–611, 2008.

J. J. Manson, A. Ma, P. Rogers et al., “Relationship between anti-dsDNA, anti-nucleosome and anti-alpha-actinin antibodies and markers of renal disease in patients with lupus nephritis: a prospective longitudinal study,” Arthritis Research & Therapy, vol. 11, no. 5, p. R154, 2009.