ORIGINAL RESEARCH article

Front. Immunol., 29 June 2018

Sec. Autoimmune and Autoinflammatory Disorders

Volume 9 - 2018 | https://doi.org/10.3389/fimmu.2018.01473

B Cell-Related Circulating MicroRNAs With the Potential Value of Biomarkers in the Differential Diagnosis, and Distinguishment Between the Disease Activity and Lupus Nephritis for Systemic Lupus Erythematosus

  • 1. Department of Nephrology, Ruijin Hospital, Shanghai Jiaotong University School of Medicine, Shanghai, China

  • 2. School of the 2nd Clinical Medical Sciences, Wenzhou Medical University, Wenzhou, Zhejiang, China

  • 3. Department of Microbiology and Immunology, Institute of Molecular Virology and Immunology, Institute of Tropical Medicine, Wenzhou Medical University, Wenzhou, Zhejiang, China

  • 4. Department of Nephrology of the First Affiliated Hospital, Wenzhou Medical University, Wenzhou, Zhejiang, China

  • 5. Department of Rheumatology of the First Affiliated Hospital, Wenzhou Medical University, Wenzhou, Zhejiang, China

  • 6. Department of Laboratory Medicine of the Second Affiliated Hospital, Wenzhou Medical University, Wenzhou, Zhejiang, China

Abstract

Our understanding of circulating microRNAs (miRNAs) related to systemic lupus erythematosus (SLE) remains very limited. In this study, we screened SLE-specific miRNAs in plasma from 42 B cell-related miRNAs by using miRNA PCR Array. The selected miRNAs were first confirmed in plasma samples from 50 SLE patients, 16 rheumatoid arthritis (RA) patients, and 20 healthy donors using qRT-PCR. We then investigated the relationship between expressions of the selected miRNAs and SLE clinical indicators. As a result, 14 miRNAs (miR-103, miR-150, miR-20a, miR-223, miR-27a, miR-15b, miR-16, miR-181a, miR-19b, miR-22, miR-23a, miR-25, miR-92a, and miR-93) were significantly decreased in the plasma of SLE patients compared with healthy controls (P < 0.05) and could act as the diagnostic signature to distinguish SLE patients from healthy donors. Six miRNAs (miR-92a, miR-27a, miR-19b, miR-23a, miR-223, and miR-16) expressed in plasma were significantly lower in SLE patients than in RA patients (P < 0.05), revealing the potentially diagnostic signature to distinguish SLE patients from RA patients. Furthermore, the downregulated expression of miR-19b, miR-25, miR-93, and miR-15b was associated with SLE disease activity (P < 0.05) while miR-15b and miR-22 expressions were significantly lower in SLE patients with low estimate glomerular filtration rate (eGFR < 60 ml/min/1.73 m2) (P < 0.05). The diagnostic potential of miR-15b for SLE disease activity and lupus nephritis (LN) with low eGFR was validated on an independent validation set with 69 SLE patients and a cross-validation set with 80 SLE patients. In summary, the signature of circulating miRNAs will provide novel biomarkers for the diagnosis of SLE and evaluation of disease activity and LN.

Introduction

Systemic lupus erythematosus (SLE) is a clinically heterogeneous autoimmune disease which affects multiple organs and systems and causes significant morbidity and mortality (1). Recently, several criteria have been developed to diagnose SLE or determine SLE activity at an early phase (2). However, owing to its pathogenesis and the exact etiology that has not been completely elucidated, to discover novel biomarkers for early diagnosis of this disease and predicting the therapeutic outcome is very urgent, which enables clinicians to treat SLE patients with the most optimally biologic therapy as early as possible.

MicroRNAs (miRNAs) are a class of non-coding small RNAs approximately having 19–25 nucleotides. miRNAs play key roles in regulating post-transcriptional gene expression by complementary pairing to their target messenger RNAs (3, 4). Over the last decade, published studies have provided strong evidence for a connection between expression of dysregulated miRNAs and development of several systemic autoimmune diseases, including SLE, which gives us new insights into the pathogenesis of SLE and a new opportunity to find novel diagnostic or therapeutic targets (5–13). In addition, miRNAs were present in serum or plasma in a remarkably stable form, and prevented degradation from endogenous RNase activity (14–17). Thus, cell-free circulating miRNAs will be important biomarkers for SLE diagnosis.

It has been demonstrated that uncontrolled over-activated B cells through their maturation into antibody-producing plasma cells are central to the pathogenesis and development of SLE (9, 18). B cells also regulate T-cell activity and immune response by acting as antigen-presenting cells or via the production of co-stimulatory molecules and proinflammatory cytokines. Therefore, B cell-related miRNAs may be an attractive target for SLE diagnosis and sequent evaluation of therapeutic outcome.

We have previously screened 72-miRNAs profiling differently expressed in B cells of SLE patients using Affymetrix miRNA 2.0 array (19). In this study, expression of 42 selected miRNAs in a large set of plasma specimens from patients with SLE and rheumatoid arthritis (RA) as well as from healthy controls (HCs) were analyzed using qRT-PCR aiming to identify the miRNAs that could potentially serve as novel serum-based biomarkers for the diagnosis of SLE and further distinguishment between disease activity and lupus nephritis (LN).

Materials and Methods

Participants and Study Design

Fifty SLE patients and 16 RA cases receiving clinical care at Department of Rheumatology, the First Affiliated Hospital of Wenzhou Medical University, were enrolled in this study between March 2013 and September of 2015. The plasma of 30 age- and sex-matched HCs were recruited from the Wenzhou local blood bank diagnosed without arthralgia, heart failure, renal failure, or autoimmune disease and free from other inflammatory conditions. Sixty-nine patients in the independent validation cohort were enrolled from March 2013 to September 2015. Additional 80 SLE patients in cross-validation set were enrolled between March 2013 and September 2015 from the Department of Nephrology, Ruijin Hospital of Shanghai Jiao Tong University School of Medicine. The research protocol had been approved by the Medical Ethics Committee of First Affiliated Hospital of Wenzhou Medical University and Ruijin Hospital of Shanghai Jiao Tong University School of Medicine. All participants were informed and gave their written consent to participate in the study. All SLE/RA patients fulfilled the American College of Rheumatology criteria for SLE/RA (20, 21). Disease activity was assessed by the Systematic Lupus Activity Measure and Systemic Lupus Erythematosus Disease Activity Index (SLEDAI) at the time of plasma collection. The SLE patients with SLEDAI >15 were polled together and set as highly active SLE group, those SLE patients with SLEDAI <4 were polled and set as stable SLE group. Modification of diet in renal disease formula was used to estimate glomerular filtration rate (eGFR). All the patient’s clinical features and treatment drugs were collected: drugs for SLE include steroids, mycophenolate mofetil, cyclophosphamide, azathioprine, MTX, tacrolimus, cyclosporine A, and Leflunomide.

The study was designed as a three-phase epidemiological approach to identify the novel circulating miRNAs associated with the disease status of SLE (Figure 1). (1) 10 plasma samples randomly selected from the enrolled SLE patients and 10 plasmas from the age- and sex-matched HCs were used to screen SLE-specific circulating miRNAs from 42 SLE B cell-associated miRNAs using customized miRNAs qPCR array. (2) Different expression of candidate circulating miRNAs was confirmed in the plasma samples comprised of 50 SLE patients, 20 HCs, and 16 patients with RA using qRT-PCR. (3) The identified miRNAs related to the disease activity and LN of SLE were re-validated on two independent SLE cohorts. The Wenzhou validation cohort included 69 SLE patients, and the Shanghai cross-validation cohort consisted of 80 SLE patients. Basic information of all participants was shown in Tables 1 and 2.

Figure 1

Table 1

Clinical characteristicSLERAHC
Number of patients501620
Sex, male/female10/403/136/14
Age37.18 ± 15.0647.62 ± 10.43
Systemic Lupus Erythematosus Disease Activity Index (SLEDAI), mean ± SD (range)6.98 ± 5.38 (0–20)
Erythrocyte sedimentation rate15 (39)14 (16)
The stage of chronic kidney disease
Early stage (CKD 1~2 stage)29 (34)
Late stage (CKD 3~5 stage)5 (34)
Autoantibodies and complement
ANA40 (40)6 (11)
Anti-double-stranded DNA17 (40)0 (11)
Low C3 level32 (46)1 (12)
Low C4 level26 (41)2 (12)
Low P level4 (24)2 (8)
C-reactive protein (mg/l)3.34 ± 6.3241.01 ± 63.17
Treatment drugs
Steroids50
Immunosuppressive drugs28

Clinical details of the patients with systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), and healthy controls (HCs).

In the phase of candidate microRNA selection, the 10 samples were from 50 samples of SLE patients.

Table 2

Clinical characteristicIndependent validation setIndependent validation set
Number of patients6980
Sex, male/female9/6013/67
Age32.97 ± 11.7138.96 ± 14.66
Incipient14 (69)30 (80)
Systemic Lupus Erythematosus Disease Activity Index (SLEDAI), mean ± SD (range)13.30 ± 6.82 (0–30)10.66 ± 7.97 (0–33)
Erythrocyte sedimentation rate35 (61)
The stage of chronic kidney disease
Early stage (CKD 1~2 stage)64 (69)60 (80)
Late stage (CKD 3~5 stage)5 (69)20 (80)
Autoantibodies and complement
ANA69 (69)71
Anti-double-stranded DNA42 (58)54
Low C3 or C4 level52 (67)41 (79)
Low C4 level50 (63)53
Low P level7 (37)2
C-reactive protein (mg/l)12.69 ± 18.8723.98 ± 54.74
Treatment drugs
Steroids6777
Immunosuppressive drugs4661

Clinical features of systemic lupus erythematosus cases on independent validation set from Wenzhou and cross-validation set from Shanghai.

Disease manifestations were defined and scored according to the SLEDAI.

Preparation of Plasma Samples and RNA Isolation

Blood samples were collected into EDTA-2K containing tubes. Blood cells were removed by a two-step centrifugation protocol performed at 4°C (3,000 rpm for 10 min, then 12,000 rpm for 10 min) to obtain platelet-poor, cell-free plasma. All samples were stored at −80°C until analyzed. Plasma RNA was isolated using Trizol® LS Reagent (Ambion, US) and Direct-zol™ RNA kit (Cat#2060&2062, Zymo Research, US) according to the manufacturer’s protocol. To normalize possible sample-to-sample variation caused by RNA isolation, 1 nmol (total volume of 1 µl) of synthetic C. elegans miRNAs cel-miR-39 (MIMAT0000010, Seq: 5′-UCACCGGGUGUAAAUCAGCUUG-3′, GenePharma, China) was added to each denatured sample. The isolated RNAs were digested by Dnase I cocktail to remove the residual DNA, and then were collected in 25 μl DNase/RNase-free water using Direct-zol™ RNA kit. Concentration of RNA was measured by NanoDrop 2000 (Thermo Fisher scientific, US). Isolated RNA was kept at −80°C or immediately used for reverse transcription.

Screening of B Cell-Related miRNAs in SLE Plasma Using the Customized qRT-PCR miRNA Array

Two independent experiments were performed to detect the different expressions of B cell-related miRNAs in SLE plasma using the customized miRNAs array (CT Bioscience, Jiangsu, China). To eliminate the difference among different samples as well as possible, five plasma samples from SLE patients or HC were mixed as one pool. The panel of the customized miRNA array contained 42 B cell-related miRNAs selected from our previous study (19). First, aliquot 100 ng of plasma RNA was added the poly A tail and then transcribed into cDNA according to the manufacturer’s protocol. qRT-PCR detection using the customer PCR chip contained 42 selected miRNA panels was performed on an Applied BioSystems 7500 Real-Time PCR system (Life Technologies). Data were automatically analyzed with SDS Relative Quantification Software version 2.2.2 (Life Technologies). U6, let-7d and let-7g were set as internal controls for plasma miRNA (22), and cel-miR-39 was set as the external control for plasma RNA isolation. The miRNAs that the Ct values were greater than 35 or with unqualified dissolution curve were removed from all datasets. Average Ct values of the duplicate analysis of each miRNA were then subtracted from the average Ct value of the cel-miR-39 or the internal control of U6, let-7d and let-7g for that particular sample, yielding the ΔCt values. The relative expression level of each miRNA was measured through the equation 2−ΔCt. The experiment was repeated three times.

Detection of miRNAs Expression

Plasma miRNA expressions of 119 samples from the First Affiliated Hospital of Wenzhou Medical University were detected. First, 100 ng plasma RNA was reversely transcribed using miScript II RT kit (Qiagen, Valencia, German) according to the manufacturer’s protocol. Stem-loop Real-time qPCR was performed using the miScript SYBR Green PCR Kit (Cat#218073, Qiagen, German) according to the manufacturer’s protocol. The expressions of miRNAs in the 80 plama samples from Ruijin Hospital were detected using mirVana qRT-PCR miRNA Detection Kit (Cat#AM1558, Ambion, US). All RT-qPCR reactions were carried out on an Applied BioSystems 7500 Real-Time PCR system (Life Technologies). All reactions ran in triple. The data were analyzed through the comparative threshold cycle (Ct) method. For expression analysis, the experiment was designed to use the cel-miR-39 as the external control; therefore, the relative quantification of plasma miRNAs was calculated using the equation: amount of target miRNA expression = 2−ΔCt, ΔCt = Ctplasma miRNAs − Ctcel-miR-39. The results were magnified 10,000 times.

Statistical Analysis

Data were presented as the median and quartile. SPSS 22.0 software was applied for statistical analysis. Statistical significance between groups was determined by Mann–Whitney U-test, Student’s t-test, chi-square test, Fisher’s exact test, Pearson product-moment correlation coefficient, and logistic regression analysis with Dunnett’s multiple comparison test as appropriate. P value less than 0.05 was considered statistically significant. Receiver operating characteristic (ROC) curve analysis, plotting the true positive rate (sensitivity) versus the false positive rate (1 − specificity) at various threshold settings were performed for plasma miRNAs, and the areas under curve (AUC) was calculated with SPSS22.0. The maximum of the sum of true positive rate and false positive rate were calculated, and cutoff value with higher specificity was selected. Expression graphs and ΔCt values were analyzed using GraphPad Prism version 5.04 software.

Results

Screening and Verifying SLE-Specific Circulating miRNAs

According to our previous results of miRNAs profiling differentially expressed in B cells of SLE patients and combined with reported literature, 42 B cell-associated miRNAs (Tables S1–S4 in Supplementary Material) were selected to customize an miRNAs panel array for detecting plasma miRNAs associated with SLE. Results from the two set of experiments revealed that expression of 15 miRNAs had twofold change in the plasma between SLE patients and HC. One miRNA, miR-155, was upregulated in the plasma of SLE patients, and the other 14 miRNAs (miR-103, miR-150, miR-20a, miR-223, miR-27a, miR-15b, miR-16, miR-181a, miR-19b, miR-22, miR-23a, miR-25, miR-92a, and miR-93) were downregulated (Table 3). Based on the screening results of SLE-specific circulating miRNAs, we further analyzed the expressions of 15 selected miRNAs in plasma of 50 SLE patients, 16 RA patients, and 20 HCs. As shown in Figure 2, expressions of 14 miRNAs (miR-103, miR-150, miR-20a, miR-223, miR-27a, miR-15b, miR-16, miR-181a, miR-19b, miR-22, miR-23a, miR-25, miR-92a, and miR-93) were significantly downregulated in the plasma of SLE patients compared with those in HCs (P < 0.05), consistent with the result of our customized miRNAs panel array screening. The expression of miR-155 showed no significant difference in the plasma between SLE patients and HCs (P > 0.05). Furthermore, except for miR-181a, 13 miRNAs were also significantly downregulated in the plasma of RA patients compared with HCs (P < 0.05). In addition, six miRNAs (miR-23a, miR-92a, miR-223, miR-27a, miR-16, and miR-19b) in the plasma of SLE patients were significantly lower than those in RA patients (P < 0.05).

Table 3

miRNAFold change (SLE/HC, 2−Δct) A
miRNAFold change (SLE/HC, 2−Δct) B
U6let-7dlet-7gcel-miR-39U6let-7dlet-7gcel-miR-39
miR-1555.591.984.433.03miR-15512.9914.9212.1214.92
miR-1030.610.210.480.33miR-1030.080.10.080.1
miR-1500.510.180.410.28miR-1500.70.810.650.81
miR-20a0.630.220.50.34miR-20a0.610.70.570.7
miR-2230.650.230.510.35miR-2230.280.320.260.32
miR-27a0.570.20.450.31miR-27a0.40.460.370.46
miR-15b0.030.010.030.02miR-15b0.0030.0040.0030.004
miR-160.010.010.010.01miR-160.010.010.010.01
miR-181a0.090.030.070.04miR-181a0.280.320.260.32
miR-19b0.190.070.150.1miR-19b0.160.180.150.18
miR-220.0100.010.01miR-220.10.110.090.11
miR-23a0.220.070.170.12miR-23a0.10.110.090.11
miR-250.160.050.120.08miR-250.460.530.430.53
miR-92a0.040.010.030.02miR-92a0.040.050.040.05
miR-930.030.010.020.01miR-930.160.180.150.18

Differential expression of B cells-associated microRNAs (miRNAs) in the plasma of systemic lupus erythematosus (SLE) patients and healthy control (HC).

miRNA relative expression in the two cohorts’ plasma of SLE and HC. Cel-miR-39 is an external reference, U6 is the most common internal reference, and let-7d and let-7g are common internal reference. The fold change in the columns is the ratio of expression in patients with SLE versus controls when using different normalization. All P values were less than 0.05.

Figure 2

ROC Curve Analysis of 14 Candidate Plasma miRNAs as SLE Diagnosis

Based on the results of SLE-specific circulating miRNA, we investigated whether miRNAs could be used as a new diagnostic marker of SLE using the ROC curve. As shown in Figure 3 and Table 4, the AUC for miR-150, miR-92a, miR-27a, miR-19b, miR-25, miR-23a, miR-93, miR-181a, miR-22, miR-223, miR-16, miR-20a, miR-15b, and miR-103 in predicting SLE from HCs were 0.832, 0.954, 0.948, 0.837, 0.951, 0.881, 0.823, 0.708 0.891, 0.93, 0.951, 0.887, 0.922, and 0.941, respectively. The diagnostic sensitivity for SLE distinguish from HCs was 0.61~0.97 and the specificity was 0.61~1. In addition, the AUCs for miR-92a, miR-27a, miR-19b, miR-223, miR-23a, and miR-16 in predicting SLE from RA were 0.665, 0.873, 0.679, 0.699, 0.761, and 0.733, respectively (Figure 4; Table 5). The diagnostic sensitivity for SLE distinguish from RA patients was 0.59~0.87 and the specificity was 0.71~0.94.

Figure 3

Table 4

miRNAROC area (AUC)95% upper bound95% lower boundSpecificitySensitivity
miR-1500.8320.9460.7180.950.614
miR-92a0.9541.0030.9050.90.898
miR-27a0.94810.89610.818
miR-19b0.8370.9410.7330.8420.791
miR-250.9511.0060.89610.861
miR-23a0.8810.9610.8020.8890.745
miR-930.8230.9570.6880.850.786
miR-181a0.7080.8690.5480.6110.833
miR-220.8910.9680.8150.90.827
miR-2230.930.9920.8690.850.884
miR-160.9511.020.8830.9440.974
miR-20a0.8871.0290.7460.8330.974
miR-15b0.9220.9840.8590.950.824
miR-1030.9411.0020.88110.821

Receiver operating characteristic (ROC) curve (AUC) for 14 microRNAs (miRNAs) in discriminating systemic lupus erythematosus patients from healthy control.

Figure 4

Table 5

miRNAROC area (AUC)95% upper bound95% lower boundSpecificitySensitivity
miR-92a0.6650.7900.5390.9380.592
miR-27a0.8730.9660.7790.7730.867
miR-19b0.6790.8180.5400.7210.625
miR-2230.6990.8500.5490.7210.688
miR-23a0.7610.9000.6220.7060.812
miR-160.7330.8820.5830.8750.625

Receiver operating characteristic (ROC) curve (AUC) for six microRNAs (miRNAs) in discriminating systemic lupus erythematosus patients from rheumatoid arthritis patients.

Correlation Between Plasma miRNAs Expression and SLE Clinical Variables

We next investigated the relationship between miRNAs expressions in plasma and established clinical variables of SLE patients (Table 6). There were only several miRNAs associated with some anti-nuclear antibodies although 14 selected miRNAs were SLE B cell-associated miRNAs. The expression levels of miR-19b, miR-22, miR-23a, and miR-25 in SLE patients with positive anti-β2GP1 were significantly higher than negative anti-β2GP1 patients (P < 0.05). The expression level of miR-15b was significantly higher in SLE patients with positive anti-SSB than those with negative anti-SSB (P = 0.014). There was no correlation between the expression of other miRNA and anti-U1RNP, anti-histone, anti-SSA, anti-Sm, and anti-Ribp in the plasma of SLE. Analysis of the relationship between expression of miRNAs in plasma and disease activity of SLE patients revealed that four miRNAs (miR-15b, miR-93, miR-25, and miR-19b) in the active group were significantly different from the stable group. Expressions of these four miRNAs in active state of SLE were all significantly lower than those in the stable SLE group (Figure 5A). According to the analysis of ROC curve (Figure 5B), AUCs for miR-15b, miR-93, miR-25, and miR-19b in determining active SLE patients were 0.619 (95% CI: 0.521–0.860), 0.702 (95% CI: 0.544–0.860), 0.688 (95% CI: 0.518–0.859), and 0.673 (95% CI: 0.521–0.825), respectively. In addition, data indicated that expressions of two plasma miRNAs (miR-15b and miR-22) were associated with renal injury of SLE. miR-15b and miR-22 were significantly lower in SLE patients with low eGFR (eGFR < 60 ml/min/1.73 m2) than those with better renal function (eGFR > 60 ml/min/1.73 m2) (P = 0.01 and 0.012, respectively). Analysis of ROC curve showed that the AUCs for miR-15b and miR-22 in determining SLE patients with renal damage (eGFR < 60 ml/min/1.73 m2) were 0.872 (95% CI: 0.714–1.000) and 0.859 (95% CI: 0.677–1.000), with sensitivity of 0.800 or 1.000 and specificity of 0.800 or 0.600, respectively (Figure 6). It did not have any significant association between expressions of plasma miRNAs and other SLE clinical variables (Table 7), except for miR-23a, miR-25, miR-16, miR-15b, miR-150, and miR-223. miR-23a was associated with gender, miR-25 was associated with CRP, and miR-15b, miR-150, and miR-223 were associated with ALT. Expression of plasma miR-16 was significantly lower in recurrent SLE patients than in SLE patients with first onset (P = 0.045). Analysis of treatment drugs showed that all SLE patients used steroids and 28 patients used immunosuppressive drugs among 50 SLE patients.

Table 6

miRNAAnti-β2GP1
Negative (n = 8)Positive (n = 16)P value
miR-19b58.85 (3.99, 163.99)376.67 (108.53, 644.82)0.049
miR-22165.04 (19.78, 472.04)345.58 (170.32, 520.86)0.043
miR-23a70.92 (17.95, 199.81)141.87 (108.64, 175.11)0.031
miR-25136.57 (22.04, 436.89)419.55 (308.56, 530.54)0.016

miRNAAnti-SSB
Negative (n = 20)Positive (n = 19)P value

miR-15b90.78 (29.75, 241.64)369.38 (108.46, 653.18)0.014

Correlation of 14 microRNAs (miRNAs) with anti-nuclear antibodies in systemic lupus erythematosus patients.

The results are magnified 10,000 times.

Figure 5

Figure 6

Table 7

miRNASex
NMaleNFemaleP value
miR-23a849.45 (13.65, 146.85)3812.22 (2.72, 43.87)0.046

miRNACRP (mg/l)
N<8N>8P value

miR-2518101.34 (9.50, 559.01)615.19 (0.23, 164.39)0.046

miRNAALT (U/l)
N5–55N>55P value

miR-15b35187.16 (34.70, 485.44)3479.41 (411.09, 5,287.22)0.088
miR-1503220.54 (8.41, 64.96)2166.42 (94.92, 237.91)0.048
miR-223277.36 (1.73, 33.85)21,258.80 (59.06, 2,458.75)0.048

miRNAIncipient
NYesNNoP value

miR-16328.14 (17.71, 265.16)281.12 (0.28, 12.80)0.045

Correlation of 14 microRNA (miRNA) expression with other clinical features of systemic lupus erythematosus.

Double Validation of miR-15b as a Biomarker for SLE Patients With Disease Activity and LN

Based on the above analysis, miR-15b could act as a potential disease biomarker for SLE patients with disease activity and LN with low eGFR. Thus, to validate the results, we collected 69 SLE samples to further determine the expression of miR-15b. Results showed that miR-15b in the expression of SLE activity group was significant different from the stable group (P = 0.039, Figure 7A). According to the analysis of ROC curve (Figure 7B), miR-15b that separate active and stable state of the SLE patients with an AUC of the ROC curve of 0.665 (95% CI: 0.519–0.797), which determined the sensitivity (0.810) and specificity (0.604). Comparison of plasma miRNA expression level in patients with renal function revealed that the miR-15b in patients with low eGFR (eGFR < 60 ml/min/1.73 m2) was significantly decreased (P = 0.033); the analysis of ROC curve (Figure 7C) showed miR-15b to determine the AUC was 0.820 (95% CI: 0.662–0.978), which determined the sensitivity (0.609) and specificity (1.000). In addition, results obtained from another 80 SLE samples also showed that miR-15b in the expression of SLE activity group was significantly different from it in stable group (P = 0.01, Figure 7D). ROC curve analysis showed that AUC of miR-15b in predicting activity of the SLE patients was 0.696 (95% CI: 0.560–0.832), which determined the sensitivity at 0.767 and specificity at 0.586 (Figure 7E). Comparison plasma miRNA expression in patients with different eGFR level, miR-15b in patients with low eGFR (eGFR < 60 ml/min/1.73 m2) was significantly decreased (P = 0.028); based on the analysis of ROC curve, AUC of miR-15b in predicting LN with low GFR was 0.675 (95% CI: 0.508–0.843), which determined the sensitivity at 0.878 and specificity at 0.438 (Figure 7F). The results of treatment drugs showed that in the cohort of 69/80 SLE patients from Wenzhou/Shanghai, 67/77 patients used steroids and 46/61 patients used immunosuppressive drugs.

Figure 7

Discussion

MicroRNA is very stable in blood circulation, and can be preserved for a long time without degradation in vitro, which makes miRNA very suitable for using as a new serum marker for a variety of diseases (11). The stability can be partially explained by the discovery of lipoprotein complexes, including small membrane vesicles of endocytic origin (called exosomes or microvesicles), high-density lipoprotein, and circulating Argonaute 2 complexes (23–25). Here, we used an ROC curve to evaluate the results and found that 14 miRNAs with clinical diagnostic value in the diagnosis of SLE are different from the HCs, six miRNAs with clinical diagnostic value are different between SLE and RA patients. The differential expression of miRNAs in plasma of SLE patients might be used as new biomarkers. Moreover, we analyzed the correlation between SLE-specific circulating miRNA and clinical features of patients with SLE, including autoantibodies and other related indicators. Results indicated that expressions of miR-19b, miR-25, miR-93, and miR-15b in SLE activity group were significantly lower than those in the SLE stable group. In addition, further analysis indicated that expressions of miR-15b and miR-22 correlated significantly with renal damage. Further diagnostic potential of miR-15b for SLE disease activity and LN was selected and validated on an independent validation set with 69 SLE patients and a cross-validation set with 80 SLE patients.

Circulating miRNA has been reported to be differentially expressed in SLE patients, but data are inconsistent in different studies, we reviewed them in PubMed and collected together in Table 8. Carlsen et al. screened 45 plasma miRNAs in two independent cohorts using quantitative RT-PCR assays (26). They found that seven miRNAs were statistically significant in plasma from SLE patients. Expression of miR-142-3p and miR-181a was increased while expression of miR-106a, miR-17, miR-20a, miR-203, and miR-92a was decreased. A 4-miRNA signature was diagnostic of SLE, and a specific subset of miRNA profiles was associated with nephritis. Wang et al. screened the plasma miRNA profiles of SLE patients by miRNA microarrays and identified eight miRNAs by RT-qPCR (27). They found that circulating levels of miR-16, miR-223, miR-23a, miR-15b, miR-150, and miR-25 were upregulated in patients with SLE in contrast to that expressions of circulating miR-155 and miR-92a were downregulated. Surprisingly, SLE-associated miR-155 was not significantly dysregulated in the plasma of both SLE and RA patients, compared with HCs. Using miRNA PCR array, Kim et al. found that nine miRNAs were differentially expressed between the SLE and HC (28). Three miRNAs (hsa-miR-30e-5p, hsa-miR-92a-3p, and hsa-miR-223-3p) were significantly upregulated in plasma of SLE patients (P = 0.048, P = 0.039, and P = 0.046, respectively). In our study, we found the plasma levels of 14 miRNAs were downregulated. The change of miR-92a in patients with SLE was consistent with previous studies (26, 27). We, together with Carlsen and colleagues, showed that expression of miR-20a in the SLE patients was downregulated. But expression of miR-155 in plasma was not significantly different between SLE patients and HCs.

Table 8

miRNAThis studyKim et al. (28)Wang et al. (27)Carlsen et al. (26)Steen et al. (29)
Experiment subjectsPlasmaPlasmaPlasmaPlasmaPlasma
Number (SLE/HC)199/2070/4030/20131/14329/40
miR-155NS↓NS
miR-150↓NSNSNS↓
miR-20a↓↓
miR-223↓↑↑NS↑
miR-15b↓NS
miR-16↓↑NS
miR-181a↓↑
miR-23a↓NS
miR-25↓NSNS
miR-92a↓↑NS↓

The comparison of circulating miRNAs in our study with other SLE studies.

SLE/HC means SLE versus healthy control; number showed in the table means that we counted all the numbers of experiment subjects together.

NS, no significance; miRNAs, microRNAs; SLE, systemic lupus erythematosus; HC, healthy control.

Using miRNAs sequencing, Navarro-Quiroz et al. comprehensively examine the plasma abundance of miRNAs in patients with class II, III, and IV lupus nephritis (LN II, LN III, and LN IV) compared with the plasma miRNA levels in lupus patients with no nephritis (LNN) or control healthy individuals (CTL) (30, 31). They found that expression of 89 miRNAs were significantly different in patients with LN, compared with CTL individuals. 17 miRNAs were differently abundant between patients with LN and LNN group. Compared with CTL individuals, the plasma level of miR-150 was upregulated in patients with LN III while miR-19b and miR-16 were downregulated in the patients with LN II. Comparison of LN II with LNN or LN IV with CTL, expression levels of miR-92a were all upregulated. In Wang and colleagues’ study of serum samples from HCs and early/late stage LN patients were used to analyze the expression of miRNAs by microarray. MiR-223 was upregulated and miR-22 was downregulated in the early LN patients versus HCs (32). On the other hand, miR-223, miR-19b, miR-22, miR-23a, miR-25, miR-92a, and miR-93 were all significantly downregulated in LN patients at the late stage compared to HCs (32).

Recently, several published studies have shown that miR-15b is involved in the pathologies of autoimmune disease. Liu et al. identified miR-15b as an important factor in Th17-associated effects and determined that expression of miR-15b was significantly downregulated in multiple sclerosis patients and in mice with experimental autoimmune encephalomyelitis (33). They also found that O-linked N-acetylglucosamine transferase is a potential target of miR-15b, enabling it to affect the transcriptional regulation of retinoic acid-related orphan receptor γT through O-linked N-acetylglucosamine glycosylation of NF-κB. Singh et al. reported that miR-15b/16 enhances the induction of regulatory T cells by regulating the expression of Rictor and mammalian target of rapamycin (34). Ren et al. used a B-cell profiling chip analysis and found that CyclinD3 was related to SLE and significantly elevated in SLE B cells (35). Furthermore, they demonstrated that activation of TLR7 dramatically increased CyclinD3 expression but significantly decreased miR-15b in B cells in vitro. They identified further that CyclinD3 is a direct target of miR-15b. Our results also showed that plasma miR-15b was decreased in patients with SLE, and was positively correlated with disease activity and LN with low eGFR. Thus, a larger sample size including different pathological types of LN, in a well-designed study will be promising for understanding of the functional role of miR-15b in future.

In conclusion, we have demonstrated here that differential expression of circulating miRNAs is a typical character between SLE patients and HCs. The plasma miR-15b correlated with SLE activity and LN may serve as a biomarker for SLE.

Data Availability Statements

All relevant data are contained within the manuscript, all datasets for this study are included in the manuscript and the supplementary files. All the raw data of this manuscript are available by the authors, without undue reservation, to any qualified researcher.

Statements

Ethics statement

This study was carried out in accordance with the recommendations of the ethical standards of the responsible committee on human experimentation (the hospital of the First Affiliated Hospital of Wenzhou Medical University, and Ruijin Hospital of Shanghai Jiao Tong University School of Medicine, China). The protocol was approved by the hospital of the First Affiliated Hospital of Wenzhou Medical University, and Ruijin Hospital of Shanghai Jiao Tong University School of Medicine, China. All subjects gave written informed consent in accordance with the Declaration of Helsinki.

Author contributions

HZ, XH, and LY: performed the experiments, analyzed and interpreted the data, and drafted the manuscript. GG: performed the experiments and statistical analysis. XL, CC, LS, and BL: acquired the data and material support. XX and NC: made contribution to the conception and design, analyzed and interpreted the data, supervised the study, provided the project funding, revised the manuscript, and finally approved the version of the manuscript for publication.

Funding

This work is supported by grants from the Zhejiang Provincial Natural Science Foundation of China (LY12H05003 and LY14H050006), the Zhejiang Science and Technology Bureau (2012C33126,2016C33214), and the Wenzhou Science and Technology Bureau (Y20170744). These sponsors provide the funding for the experiments and the collection of specimens.

Acknowledgments

The authors thank all the patients and their families, physicians, and technicians who were involved in the study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at https://www.frontiersin.org/articles/10.3389/fimmu.2018.01473/full#supplementary-material.

References

  • 1

    FrieriM. Mechanisms of disease for the clinician: systemic lupus erythematosus. Ann Allergy Asthma Immunol (2013) 110(4):228–32.10.1016/j.anai.2012.12.010

  • 2

    AletahaDNeogiTSilmanAJFunovitsJFelsonDTBinghamCOIIIet al2010 Rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum (2010) 62(9):2569–81.10.1002/art.27584

  • 3

    CarthewRWSontheimerEJ. Origins and mechanisms of miRNAs and siRNAs. Cell (2009) 136(4):642–55.10.1016/j.cell.2009.01.035

  • 4

    KimVNHanJSiomiMC. Biogenesis of small RNAs in animals. Nat Rev Mol Cell Biol (2009) 10(2):126–39.10.1038/nrm2632

  • 5

    PauleyKMChaSChanEK. MicroRNA in autoimmunity and autoimmune diseases. J Autoimmun (2009) 32(3–4):189–94.10.1016/j.jaut.2009.02.012

  • 6

    O’ConnellRMRaoDSChaudhuriAABoldinMPTaganovKDNicollJet alSustained expression of microRNA-155 in hematopoietic stem cells causes a myeloproliferative disorder. J Exp Med (2008) 205(3):585–94.10.1084/jem.20072108

  • 7

    TaganovKDBoldinMPChangKJBaltimoreD. NF-kappaB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. Proc Natl Acad Sci U S A (2006) 103(33):12481–6.10.1073/pnas.0605298103

  • 8

    O’ConnellRMKahnDGibsonWSRoundJLScholzRLChaudhuriAAet alMicroRNA-155 promotes autoimmune inflammation by enhancing inflammatory T cell development. Immunity (2010) 33(4):607–19.10.1016/j.immuni.2010.09.009

  • 9

    AnolikJH. B cell biology: implications for treatment of systemic lupus erythematosus. Lupus (2013) 22(4):342–9.10.1177/0961203312471576

  • 10

    MiaoCGYangYYHeXHuangCHuangYZhangLet alThe emerging role of microRNAs in the pathogenesis of systemic lupus erythematosus. Cell Signal (2013) 25(9):1828–36.10.1016/j.cellsig.2013.05.006

  • 11

    ShenNLiangDTangYde VriesNTakPP. MicroRNAs – novel regulators of systemic lupus erythematosus pathogenesis. Nat Rev Rheumatol (2012) 8(12):701–9.10.1038/nrrheum.2012.142

  • 12

    YanSYimLYLuLLauCSChanVS. MicroRNA regulation in systemic lupus erythematosus pathogenesis. Immune Netw (2014) 14(3):138–48.10.4110/in.2014.14.3.138

  • 13

    QuBShenN. miRNAs in the pathogenesis of systemic lupus erythematosus. Int J Mol Sci (2015) 16(5):9557–72.10.3390/ijms16059557

  • 14

    WilleitPZampetakiADudekKKaudewitzDKingAKirkbyNSet alCirculating microRNAs as novel biomarkers for platelet activation. Circ Res (2013) 112(4):595–600.10.1161/CIRCRESAHA.111.300539

  • 15

    MitchellPSParkinRKKrohEMFritzBRWymanSKPogosova-AgadjanyanELet alCirculating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci U S A (2008) 105(30):10513–8.10.1073/pnas.0804549105

  • 16

    WeberJABaxterDHZhangSHuangDYHuangKHLeeMJet alThe microRNA spectrum in 12 body fluids. Clin Chem (2010) 56(11):1733–41.10.1373/clinchem.2010.147405

  • 17

    MurataKYoshitomiHTanidaSIshikawaMNishitaniKItoHet alPlasma and synovial fluid microRNAs as potential biomarkers of rheumatoid arthritis and osteoarthritis. Arthritis Res Ther (2010) 12(3):R86.10.1186/ar3013

  • 18

    GregersenJWJayneDR. B-cell depletion in the treatment of lupus nephritis. Nat Rev Nephrol (2012) 8(9):505–14.10.1038/nrneph.2012.141

  • 19

    FangCChaoshengCHuidiZXiangyangXXiaokaiDLiSet alScreening of the miRNAs related to disease activity in peripheral blood B lymphocytes of systemic lupus erythematosus patients. J Wenzhou Med Coll (2013) 2:106–11.10.13771/j.cnki.33-1386/r.2013.02.011

  • 20

    HochbergMC. Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum (1997) 40(9):1725.10.1002/art.1780400928

  • 21

    ArnettFCEdworthySMBlochDAMcShaneDJFriesJFCooperNSet alThe American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum (1988) 31(3):315–24.10.1002/art.1780310302

  • 22

    ChenXLiangHGuanDWangCHuXCuiLet alA combination of Let-7d, Let-7g and Let-7i serves as a stable reference for normalization of serum microRNAs. PLoS One (2013) 8(11):e79652.10.1371/journal.pone.0079652

  • 23

    TabetFVickersKCCuesta TorresLFWieseCBShoucriBMLambertGet alHDL-transferred microRNA-223 regulates ICAM-1 expression in endothelial cells. Nat Commun (2014) 5:3292.10.1038/ncomms4292

  • 24

    VickersKCPalmisanoBTShoucriBMShamburekRDRemaleyAT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol (2011) 13(4):423–33.10.1038/ncb2210

  • 25

    CortezMABueso-RamosCFerdinJLopez-BeresteinGSoodAKCalinGA. MicroRNAs in body fluids – the mix of hormones and biomarkers. Nat Rev Clin Oncol (2011) 8(8):467–77.10.1038/nrclinonc.2011.76

  • 26

    CarlsenALSchetterAJNielsenCTLoodCKnudsenSVossAet alCirculating microRNA expression profiles associated with systemic lupus erythematosus. Arthritis Rheum (2013) 65(5):1324–34.10.1002/art.37890

  • 27

    WangHPengWOuyangXLiWDaiY. Circulating microRNAs as candidate biomarkers in patients with systemic lupus erythematosus. Transl Res (2012) 160(3):198–206.10.1016/j.trsl.2012.04.002

  • 28

    KimBSJungJYJeonJYKimHASuhCH. Circulating hsa-miR-30e-5p, hsa-miR-92a-3p, and hsa-miR-223-3p may be novel biomarkers in systemic lupus erythematosus. HLA (2016) 88(4):187–93.10.1111/tan.12874

  • 29

    SteenSOIversenLVCarlsenALBurtonMNielsenCTJacobsenSet alThe circulating cell-free microRNA profile in systemic sclerosis is distinct from both healthy controls and systemic lupus erythematosus. J Rheumatol (2015) 42(2):214–21.10.3899/jrheum.140502

  • 30

    Navarro-QuirozEPacheco-LugoLNavarro-QuirozRLorenziHEspana-PucciniPDiaz-OlmosYet alProfiling analysis of circulating microRNA in peripheral blood of patients with class IV lupus nephritis. PLoS One (2017) 12(11):e0187973.10.1371/journal.pone.0187973

  • 31

    Navarro-QuirozEPacheco-LugoLLorenziHDiaz-OlmosYAlmendralesLRicoEet alHigh-throughput sequencing reveals circulating miRNAs as potential biomarkers of kidney damage in patients with systemic lupus erythematosus. PLoS One (2016) 11(11):e0166202.10.1371/journal.pone.0166202

  • 32

    WangWMouSWangLZhangMShaoXFangWet alUp-regulation of serum MiR-130b-3p level is associated with renal damage in early lupus nephritis. Sci Rep (2015) 5:12644.10.1038/srep12644

  • 33

    LiuRMaXChenLYangYZengYGaoJet alMicroRNA-15b suppresses Th17 differentiation and is associated with pathogenesis of multiple sclerosis by targeting O-GlcNAc transferase. J Immunol (2017) 198(7):2626–39.10.4049/jimmunol.1601727

  • 34

    SinghYGardenOALangFCobbBS. MicroRNA-15b/16 enhances the induction of regulatory T cells by regulating the expression of rictor and mTOR. J Immunol (2015) 195(12):5667–77.10.4049/jimmunol.1401875

  • 35

    RenDLiuFDongGYouMJiJHuangYet alActivation of TLR7 increases CCND3 expression via the downregulation of miR-15b in B cells of systemic lupus erythematosus. Cell Mol Immunol (2016) 13(6):764–75.10.1038/cmi.2015.48

Summary

Keywords

systemic lupus erythematosus, circulating microRNA, diagnosis, biomarker, miR-15b

Citation

Zhang H, Huang X, Ye L, Guo G, Li X, Chen C, Sun L, Li B, Chen N and Xue X (2018) B Cell-Related Circulating MicroRNAs With the Potential Value of Biomarkers in the Differential Diagnosis, and Distinguishment Between the Disease Activity and Lupus Nephritis for Systemic Lupus Erythematosus. Front. Immunol. 9:1473. doi: 10.3389/fimmu.2018.01473

Received

17 March 2018

Accepted

13 June 2018

Published

29 June 2018

Volume

9 - 2018

Edited by

J. Michelle Kahlenberg, University of Michigan, United States

Reviewed by

Ziaur S. M. Rahman, Penn State Milton S. Hershey Medical Center, United States; Scott A. Tenenbaum, University at Albany, United States

Updates

Copyright

*Correspondence: Nan Chen, ; Xiangyang Xue,

†These authors have contributed equally to this work.

Specialty section: This article was submitted to Autoimmune and Autoinflammatory Disorders, a section of the journal Frontiers in Immunology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics