Skip to main content
Download PDF

Whole-genome transcription and DNA methylation analysis of peripheral blood mononuclear cells identified aberrant gene regulation pathways in systemic lupus erythematosus

Arthritis Research & Therapy volume 18, Article number: 162 (2016) Cite this article

Abstract

Background

Recent achievement in genetics and epigenetics has led to the exploration of the pathogenesis of systemic lupus erythematosus (SLE). Identification of differentially expressed genes and their regulatory mechanism(s) at whole-genome level will provide a comprehensive understanding of the development of SLE and its devastating complications, lupus nephritis (LN).

Methods

We performed whole-genome transcription and DNA methylation analysis in PBMC of 30 SLE patients, including 15 with LN (SLE LN+) and 15 without LN (SLE LN−), and 25 normal controls (NC) using HumanHT-12 Beadchips and Illumina Human Methy450 chips. The serum proinflammatory cytokines were quantified using Bio-plex Human Cytokine 27-plex assay. Differentially expressed genes and differentially methylated CpG were analyzed with GenomeStudio, R, and SAM software. The association between DNA methylation and gene expression were tested. Gene interaction pathways of the differentially expressed genes were analyzed by IPA software.

Results

We identified 552 upregulated genes and 550 downregulated genes in PBMC of SLE. Integration of DNA methylation and gene expression profiling showed that 334 upregulated genes were hypomethylated, and 479 downregulated genes were hypermethylated. Pathway analysis on the differential genes in SLE revealed significant enrichment in interferon (IFN) signaling and toll-like receptor (TLR) signaling pathways. Nine IFN- and seven TLR-related genes were identified and displayed step-wise increase in SLE LN− and SLE LN+. Hypomethylated CpG sites were detected on these genes. The gene expressions for MX1, GPR84, and E2F2 were increased in SLE LN+ as compared to SLE LN− patients. The serum levels of inflammatory cytokines, including IL17A, IP-10, bFGF, TNF-α, IL-6, IL-15, GM-CSF, IL-1RA, IL-5, and IL-12p70, were significantly elevated in SLE compared with NC. The levels of IL-15 and IL1RA correlated with their mRNA expression. The upregulation of IL-15 may be regulated by hypomethylated CpG sites in the promotor region of the gene.

Conclusions

Our study has demonstrated that significant number of differential genes in SLE were involved in IFN, TLR signaling pathways, and inflammatory cytokines. The enrichment of differential genes has been associated with aberrant DNA methylation, which may be relevant to the pathogenesis of SLE. Our observations have laid the groundwork for further diagnostic and mechanistic studies of SLE and LN.

Background

Systemic lupus erythematosus (SLE) is a potentially severe autoimmune disease characterized by intermittent episodes of increased disease activity that require treatment with immunosuppressive agents [1]. Renal involvement affects more than 60 % of patients with SLE, and is a major contributor to morbidity and mortality [2]. Nephritis remains one of the most devastating complications of lupus, and is an important cause of chronic renal failure despite the availability of new therapeutic regimens [3]. Optimizing the management of lupus nephritis (LN) is therefore important, both to reduce the healthcare burden to society and to improve the outcome of patients [4]. SLE should be closely monitored for renal manifestations, early diagnosis and treatment are essential for renal preservation [5]. The most commonly used marker of renal disease is urinary protein. Spot urinary protein-to-creatinine ratio is currently recommended as a marker. However, it is not sufficiently sensitive to detect early nephritis. Conventional markers of active renal disease are serum levels of anti-dsDNA, anti-C1q antibody, and complement levels. However, the sensitivity of these markers varied greatly depending on the assay system used. The specificity ranges between 50 and 75 % for active renal disease [6, 7]. A better understanding of the pathogenesis of LN is an important step in identifying more biomarkers and targeted therapeutic approaches [8].

Substantial research has helped define the pathogenic mechanisms of SLE renal manifestations [9]. A number of genetic polymorphisms have been associated with lupus, it is now recognized that multiple genes, each making only a small contribution, account for the risk of developing SLE. Several genome-wide scans have allowed major advances in the identification of genetic regions predisposing to SLE, bringing the total number of validated loci to more than 50 [10]. However, many nongenetic factors (such as environment) have now been linked to lupus, environmentally induced epigenetics may also play a critical role in autoimmune immunopathology [11]. Global epigenetic modification, such as DNA methylation and chromatin modification, are directly influenced by the environment, and play an important role in the pathogenesis of SLE [12]. Epigenetic alterations have an important effect on gene expression regulation. Many studies have investigated the association between genetics, DNA methylation, and gene expression in SLE [13–16]. In previous work, hypomethylation of global genomic DNA has been demonstrated in several autoimmune-related gene promoters in isolated T cells from SLE. Genes with aberrant DNA demethylation identified in lupus T cells include ITGAL, CD70, PRF1, and the X chromosome gene CD40LG [17–21]. Similar data are not yet available for LN, and most of the association studies to date have yielded inconsistent results across populations [22].

In the present study, we analyzed the DNA methylation distribution and the relevant mRNA expression variation in peripheral blood mononuclear cells (PBMC) from SLE patients with or without LN using whole-genome DNA methylation and gene expression arrays. We also profiled the levels of inflammatory cytokines in sera of SLE patients by bead-based multiplex assays. The integration of genome-wide gene transcription and DNA methylation expended our insight into the mechanisms affecting epigenetic alteration, gene expression, and LN susceptibility. Combined with the cytokine expression, our data provide a unique resource for identifying biomarkers that may be useful in the diagnosis of lupus nephritis.

Methods

SLE patients and controls

We studied 30 SLE patients (15 SLE with lupus nephritis (LN+) and 15 SLE without lupus nephritis (LN−)) and 25 normal controls (NC). Patients and controls were matched for age, sex, and ethnicity. All SLE patients meet the American College of Rheumatology (ACR) classification criteria for SLE [23]. The clinical features of SLE patients included in this study are shown in Table 1. SLE LN+ is defined as clinical and laboratory manifestations that meet ACR criteria (persistent proteinuria >0.5 g/day or greater than 3+ by dipstick, and/or urine cellular casts including red blood cells (RBC), hemoglobin, granular, tubular, or mixed) [2, 23]. A review of the ACR criteria has recommended that “active urinary sediment” (>5 RBCs/high-power field (HPF), >5 white blood cells (WBC)/HPF in the absence of infection) can be substituted for cellular casts [24]. The patients enrolled in the study did not have lymphopenia or cytopenia in previous history. This study was approved by the institutional review board at Xiangya Hospital, Central South University of China (Changsha, Hunan, China). All the participants in the study signed a written informed consent prior to participation [2].

Table 1 Clinical manifestation and laboratory data of SLE patients and normal controls
Full size table

DNA and RNA isolation

Peripheral blood samples were obtained from SLE patients and normal controls. The PBMC were isolated from heparinized blood by density gradient centrifugation over Ficoll-Paque Plus (GE Healthcare, Piscataway, NJ, USA). Total RNA was isolated from PBMC by standard phenol–chloroform extraction using Trizol reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) according to the manufacturer’s instructions and the concentration was measured on Nanodrop ND-1000 Spectrophotometers (Thermo Fisher Scientific, Waltham, MA, USA). RNA quality was checked on Bioanalyzer Nanochip (Agilent Technologies, Santa Clara, CA, USA) and the samples with RNA integrity number (RIN) >7 were required for microarray analysis. Genomic DNA was isolated from whole blood cells using genomic DNA extraction kits (Life Technologies, Gaithersburg, MD, USA) and DNA integrity was analyzed by agarose gel electrophoresis.

Genome-wide DNA methylation analysis

DNA methylation status of 485,000 CpG sites across the whole genome was analyzed using the Illumina HumanMethylation 450 BeadChip array. The array covers 99 % of RefSeq genes, with an average of 17 CpG sites per gene region distributed across the promoter, 5′UTR, first exon, gene body, and 3′UTR. It covers 96 % of CpG islands, with additional coverage in island shores and the regions flanking them (Illumina, Inc., San Diego, CA, USA). Genomic DNA (1 μg) extracted from PBMC was bisulfite converted using EZ DNA Methylation kit (Zymo Research Corp, Orange, CA, USA. cat. no. D5004), and utilizing a cyclic denaturation step during the conversion reaction. Four microliters of bisulfite-converted DNA were used for hybridization on Infinium HumanMethylation 450 BeadChip kit. According to the Illumina Infinium HD Methylation protocol, the following steps were performed, whole-genome amplification, end-point fragmentation, precipitation, and resuspension. The resuspended samples were hybridized onto HumanMethylation 450 BeadChips at 48 °C for 18 hours. Then unhybridized and nonspecifically hybridized DNA were washed away, followed by a single nucleotide extension and repeated rounds of staining. Finally, the BeadChip was washed, coated, and scanned. After scanning, the intensities of images were extracted using GenomeStudio Methylation module software (Illumina, Inc.). The methylation score for each CpG was represented as a β value according to the fluorescent intensity ratio. A β value may take any value between 0 (nonmethylated) and 1 (completely methylated). Expression values were extracted, filtered, and normalized using an R package (RnBeads) [25]. The normalized methylation data were analyzed by significance analysis of microarrays (SAM, Stanford University) [26]. The significant differentially methylated sites with FDR <0.05, fold change >1.2 were selected.

Gene expression studies

Illumina HumanHT-12 v4.0 Expression Beadchips (Illumina, Inc.) were used for mRNA transcription profiling. The platform contains 47,323 transcripts. Each RNA sample was amplified using the Illumina TotalPrep RNA Amplification kit (Thermal Fisher Scientific) with biotin-UTP (Enzo Life Sciences, Inc., Farmingdale, NY, USA) labeling. The TotalPrep Illumina RNA Amplification kit uses T7-oligo(dT) primer to generate single-stranded cDNA followed by a second-strand synthesis to generate double-stranded cDNA, purifies through spin column, uses T7 RNA polymerase to synthesize biotin-labeled cRNA, and purifies again. The cRNA was then quantified using the ND-1000 Spectrophotometer (Thermo Fisher Scientific). A total of 1.5 μg cRNA from each sample was hybridized on each array using standard Illumina protocols with streptavidin-Cy3 (Amersham Biosciences Corp., Piscataway, NJ, USA) being used for detection. Slides were scanned on an Illumina HiScan scanner and analyzed using GenomeStudio (Illumina, Inc.). The raw files of Illumina HumanHT-12 v4.0 Expression array were extracted and normalized (Quantile) in GenomeStudio Module. Normalized raw data were filtered by detection p value and median centered by using the Genespring software (Silicon Genetics, Redwood City, CA, USA). Expression data was performed by t test without any error correction on the samples. Genes were considered statistically significant at p < 0.05.

Multiplex cytokine assays

Diluted plasma specimens were prepared for analysis in a 96-well plate utilizing the Bio-plex Human Cytokine 27-plex assay (Bio-Rad Laboratories, Veenendaal, The Netherlands) according to the manufacturer’s instructions. The 27-plex assay kit contains beads conjugated with monoclonal antibodies specific for interleukin (IL)-17A, regulated on activation, normal T cell expressed and secreted (RANTES), platelet-derived growth factor-BB (PDGF-BB), IL-2, IL-4, IL-13, basic fibroblast growth factor (bFGF), macrophage inflammatory protein 1 alpha (MIP-1α), IFN-gamma-inducible protein 10 (IP-10), monocyte chemoattractant protein 1 (MCP-1), granulocyte colony-stimulating factor (G-CSF), Eotaxin, tumor necrosis factor alpha (TNF-α), IL-6, IL-9, IL-15, granulocyte-macrophage colony-stimulating factor (GM-CSF), IL-1RA, IL-5, IL-12p70, interferon (IFN)-ϒ, IL-1β, IL-7, IL-8, MIP-1β and vascular endothelial growth factor (VEGF). Analyses were quantified using a Magpix analytical test instrument, which utilizes xMAP technology (Luminex Corp., Austin, TX, USA) and xPONENT 4.2 software (Luminex Corp.). Concentrations of cytokines (pg/ml) were determined on the basis of the fit of a standard curve for mean fluorescence intensity versus pg/ml.

Real-time PCR validation of the differential gene expression

Differentially expressed genes were validated using Taqman assays on a 7900HT Fast Real-Time PCR system (Applied Biosystems, Foster City, CA, USA). The assay was performed by using a TaqMan RNA-to-CT 1-Step kit (Applied Biosystems) in a total volume of 20 μl, which contained a final concentration of 900 nM sense and antisense primers, 250 nM Taqman gene probe, 1 × TaqMan RT Enzyme Mix, and 1 × TaqMan RT-PCR Mix. The cDNA amplification were monitored using 7900HT Fast Real-Time PCR system under the conditions of 48 °C for 15 min, 95 °C for 10 min, and 40 cycles of 95 °C for 15 s and 60 °C for 1 min. This assay was carried out in triplicate for each sample, including a no-template control. The relative quantity (RQ) of the gene expression in each sample was calculated by normalizing to housekeeping gene GAPDH.

Pathway analysis

The Ingenuity Pathway Analysis (IPA) tool was used to identify gene networks, functions, and canonical pathways. The IPA software was used for its ability to analyze mRNA data in the context of known biological response and regulatory networks. The score for each network was derived from a p value that indicates the expected likelihood of the focus genes being present in a network compared with that expected by chance. The software determines the probability that each biological function assigned to that data set was due to chance alone (http://www.ingenuity.com/).

Statistical analysis

Data are shown as mean values and standard errors of the means. The univariate comparisons were done using a one-way ANOVA or two-sample t test, and multivariable adjustments were done using ANCOVAR to simultaneously adjust for age, SLE Disease Activity Index (SLEDAI) score, and the health status satisfaction score. All p values were from two-sided tests. Correlations between pairs of continuous variables were done using Pearson’s R. Dichotomized variables were compared using Fisher’s exact chi-square test. Bonferroni correction was applied in some of the analyses. p values of < 0.05 were considered significant. All analyses were performed using SAS version 9.3 software (SAS Institute Inc., Cary, NC, USA). Graphics were carried out using Prism version 6.0 (GraphPad, San Diego, CA, USA).

Results

Aberrant mRNA transcription in SLE

First, we investigated the gene expression profiles in PBMC of 30 SLE patients (SLE), including 15 SLE with LN (SLE LN+) and 15 SLE without LN (SLE LN−), and 25 matched normal controls (NC) using Illumina Human HT-12 v4.0 Beadchips. By comparing the gene transcription levels between SLE (including SLE LN+ and SLE LN−) and NC, we identified 879 upregulated and 855 downregulated genes which are significant (p < 0.05) in SLE as compared with NC (Fig. 1a). Further analysis was performed to compare the gene expression differences between SLE LN+ vs. NC and SLE LN− vs. NC. There are 1000 upregulated and 1028 downregulated in SLE LN+ vs. NC, and 834 upregulated and 806 downregulated genes in SLE LN− vs. NC (Fig. 1a). By Venn diagram analysis, we were able to identify 552 upregulated genes and 550 downregulated genes which are common in all three comparisons (Fig. 1a, Additional file 1: Table S1, and Additional file 2: Table S2), indicating that these common genes are consistently differentially regulated in SLE as compared to the NC group. The profile of the differentially expressed genes in all samples was shown as a heat map in Fig. 1b. The gene expression microarray data has been submitted to the GEO database with accession number GSE81622. The list of the differentially expression genes in each comparison was shown in Additional file 1: Table S1 and Additional file 2: Table S2.

Fig. 1
figure 1

Differentially expressed genes in the PBMC of SLE patients compared with normal controls. a Venn diagram shows the number of upregulated genes (left panel) and downregulated genes (right panel) in three comparisons: SLE vs. NC, SLE LN− vs. NC, and SLE LN+ vs. NC. There are 552 upregulated and 550 downregulated genes in common in all three comparisons. b Heat map shows hierarchical clustering of the commonly up- and downregulated genes in the three comparisons. NC normal controls, SLE systemic lupus erythematosus, SLE LN + SLE with lupus nephritis, SLE LN − SLE without lupus nephritis

Full size image

Aberrant DNA methylation in SLE

Next, we performed a global DNA methylation analysis on the same sets of SLE patients and normal controls using a gene chip containing 485,000 CpG sites across the whole genome. Figure 2 showed the numbers of differentially methylated CpG sites in three comparisons: SLE vs. NC, SLE LN+ vs. NC, and SLE LN− vs. NC (Fig. 2a). There are 1813 hypermethylated sites and 3785 hypomethylated sites which are common in all three comparisons (Fig. 2a, Additional file 3: Table S3 and Additional file 4: Table S4). The distribution and location of the hypo- and hypermethylated sites are shown is Fig. 2b. Nearly half of the differentially methylated sites (42.95 % of the hypomethylated and 42.99 % for hypermethylated sites) are located on or in the vicinity of the promoter regions, including TSS200 and TSS1500, which is highly suggestive of methylation as a regulatory mechanism on gene transcription. A high percentage of altered methylation sites are also clustered within the gene body, 5′UTR, and first exon (Fig. 2b). The list of the differentially methylated CpG sites in each comparison was shown in Additional file 3: Table S3 and Additional file 4: Table S4. The DNA methylation for all samples has been submitted to the GEO database with accession number GSE81622.

Fig. 2
figure 2

Differentially methylated CpG sites in PBMC of SLE patients compared with normal controls by DNA methylation analysis. a Venn diagram showing the number of hypomethylated (left panel) and hypermethylated CpG sites (right panel) in three comparisons, SLE vs. NC, SLE LN+ vs. NC, and SLE LN− vs. NC. A total of 1813 hypomethylated CpG sites and 3785 hypermethylated CpG sites are in common in all three comparisons. b The distribution of differentially methylated CpG sites on the genome. The bar charts showed the genomic distribution of hypomethylated (left panel) and hypermethylated (right panel) probes across different genomic regions: TSS1500 and TSS200 (probes located within 1500 and 200 bp from transcription start site, respectively); 5′UTR region; first exon; gene body and 3′UTR region. c Heat maps depicting cluster of the 1065 hypomethylation sites in 333 of 552 upregulated genes, which have one or more hypomethylated CpG sites, and 3024 hypermethylation sites in 487 of 550 downregulated genes, which have one or more hypermethylated sites. NC normal controls, SLE systemic lupus erythematosus, SLE LN + SLE with lupus nephritis, SLE LN − SLE without lupus nephritis

Full size image

Correlation of gene expression and DNA methylation in SLE

In order to determine the potential effects of the aberrant DNA methylation on the abnormal gene expression in SLE, we investigated the methylation status of the 552 commonly upregulated genes and 550 commonly downregulated genes in SLE. Among the 552 upregulated genes, 333 genes (60.5 %) have one or more hypomethylated CpG sites (total 1065 CpG site in 333 genes as shown in Fig. 2c, left panel; and Additional file 5: Table S5); while among the 550 downregulated genes, hypermethylated CpG sites were found in 479 genes (87.1 %) (Fig. 2c, right panel; and Additional file 5: Table S5). Our data suggested that DNA methylation may play an important role in the regulation of gene expression in the development of SLE.

Functional analysis on the genes differentially expressed in SLE

To decipher the molecular pathways associated with dysregulated genes in SLE, we performed Ingenuity Pathway Analysis (IPA) on the differentially expressed genes which are common in all three comparisons. The molecular pathways most significantly correlated with the 552 upregulated genes are showed in Fig. 3a. The top five pathways enriched in the upregulated genes including interferon signaling (1.48E-07), role of pattern recognition receptors in recognition of bacteria and viruses (1.09E-06), altered T cell and B cell signaling in rheumatology arthritis (9.43E-05), cell cycle: G2/M DNA damage checkpoint regulation (2.56E-04) and toll-like receptor signaling (6.58E-04) (Fig. 3a). Figure 3b exhibited the molecular pathways related with the 550 downregulated genes, including natural killer cell signaling (2.84E-19), T cell receptor signaling (3.89E-10), iCOS-iCOSL signaling in T helper cells (2.19E-09), cytotoxic T lymphocyte-mediated apoptosis of target cells (7.38E-08) and role of NFAT in regulation of the immune response (8.89E-08) (Fig. 3b).

Fig. 3
figure 3

Gene interaction networks among the common up- and downregulated genes in the PBMC of SLE using IPA software. a Top ten pathways showing inter-relationships among 552 upregulated genes in SLE. b Top ten pathways showing inter-relationships among the 550 downregulated genes in SLE. The top ten pathways are represented by the red bar. The lines connecting the bars depict possible gene interactions in the network. The expected likelihood of the gene involved in the pathway was indicated by p value in the network

Full size image

DNA hypomethylation related with the overexpression of IFN-related genes in SLE

Among the molecular pathways modulated by aberrantly expressed genes in SLE, type I IFN (IFN-I) signaling pathway showed highest association with the upregulated genes in SLE. The IFN-related genes which are highly expressed in SLE were significantly enriched in this pathway (Fig. 4a). Nine of the IFN-I-related genes, including IFI35, IFIT1, IFIT3, IFITM2, IFITM3, IFNGR2, myxovirus (influenza) resistance 1 (MX1), OAS1 and STAT1, showed a stepwise increase from NC to SLE LN− to SLE LN+, although the difference between SLE LN− and SLE LN+ on these genes (except MX1) was not statistically significant (Fig. 4b), suggesting that the dysregulation of IFN pathway genes may play an important role in the pathogenesis of SLE. The gene expression data on some IFN genes were further confirmed by real-time PCR and data was shown in Additional file 6: Figure S1. Based on the average expression of these nine IFN genes, we calculated an IFN score for each sample. The mean IFN score for SLE patients (15.23 ± 5.95) was significantly greater than NC (9.00 ± 2.94) (p < 0.0001). The IFN score for SLE LN+ (22.57 ± 12.84) was slightly higher than SLE LN− (16.13 ± 8.36) but not statistically significant (p > 0.1) (Fig. 4c).

Fig. 4
figure 4

Upregulated genes in PBMC of SLE were significantly enriched in the interferon signaling pathway. a A diagram showing the interferon signaling pathway and upregulated genes involved in the network. b Array-based mRNA expression represented by normalized intensity in SLE LN+, SLE LN−, and normal controls.*Indicates p < 0.05 as compared to normal controls. c Interferon scores for SLE patients and controls. The IFN score for each sample was derived from the combined fold change of relative intensity for each gene. The mean intensity in normal controls for each gene was calculated and used as baseline intensity for the gene. Then the intensity of mRNA expression for each gene in patients or controls was normalized with the baseline intensity and resulted in the fold change of relative intensity. The interferon score was a sum of relative intensity for the nine genes (shown in Fig. 5b) that were significantly enriched in the interferon pathway. d Hierarchical clustering of the differentially methylated CpG sites on the upregulated genes associated with the interferon signaling pathway. Hypomethylated CpG sites were identified in all the upregulated IFN genes. e The CpG sites in MX1, IFIT1 and IFIT3 genes were hypomethylated in SLE patients comparing with NC. *Indicates p < 0.05. f The mRNA expression of IFN genes (IFT1, IFT3 and MX1) were reversely correlated with the methylation status of the CpG sites in PBMC. NC normal controls, SLE systemic lupus erythematosus, SLE LN + SLE with lupus nephritis, SLE LN − SLE without lupus nephritis

Full size image

In order to understand if the increased expression of IFN signaling genes in SLE were regulated by altered DNA methylation, we examined the methylation status of the CpG sites on the upregulated IFN genes in all NC and SLE samples. Among the 47 CpG sites in the nine IFN-related genes, 41 CpG sites exhibited decreased methylation (hypomethylation) in SLE patients compared with NC on mean β value (p < 0.05), including two sites on four genes (IFI35, IFITM2, IFIT1 and OAS1), three on IFIT3, five on IFITM3, six sites on STAT1, eight sites on IFNGR2 and 11 sites on MX1 (Fig. 4d). Shown in Fig. 4e are the six CpG sites which exhibited the greatest methylation change on the IFN genes between SLE and NC, four on MX1 (cg22862203, cg26312951, cg16785077, cg21549285), one on IFIT3 (cg06188083) and one on IFIT1 (cg05552874) (Fig. 4e). The methylation of these sites reversely correlated with the transcription levels of the corresponding genes in PBMC (Fig. 4e, f).

DNA hypomethylation associated with higher transcription of the genes in the toll-like receptor signaling pathway in SLE

In addition to interferon signaling pathway, we noted that the expression of several genes related to the toll-like receptor signaling (TLR) pathway were also elevated in SLE. TLRs play an essential role in the innate immune system and regulate a pro-/anti-inflammatory balance (Fig. 5a). We found that mRNAs for cluster of differentiation 14 (CD14), eukaryotic translation initiation factor 2-alpha kinase 2 (EIF2AK2), IL1RN, interleukin-1 receptor-associated kinase 3 (IRAK3), TLR1, TLR4, TLR7 and TLR8 were significantly upregulated in SLE compared with NC (Fig. 5b). The expression of the TLR7 and EIF2AK2 genes were higher in SLE LN+ as compared to SLE LN−, but the change was not statistically significant (Fig. 5b).

Fig. 5
figure 5

Upregulated genes in PBMC of SLE were significantly enriched in the toll-like receptor signaling pathway. a A diagram showing the toll-like receptor signaling pathway and upregulated genes involved in the network. b Array-based mRNA expression represented by normalized intensity in SLE LN+, SLE LN−, and normal controls.*Indicates p < 0.05 as compared to normal controls. c The CpG sites in EIF2AK2 and IRAK3 genes were hypomethylated in SLE patients compared with NC. *Indicates p < 0.05. NC normal controls, SLE systemic lupus erythematosus, SLE LN + SLE with lupus nephritis, SLE LN − SLE without lupus nephritis

Full size image

We then examined the DNA methylation status on the CpG sites of these differentially expressed genes. There are multiple CpG sites (eight on IRAK3, two on EIF2AK2, four on TLR4, three on IL-1RA, two on TLR8 and one on TLR7), which showed different methylation between SLE and NC. The most significantly hypomethylated sites were found on EIF2AK2 (cg14126601 and cg17326313) and IRAK3 (cg19995654, cg12866960, and cg07914866) (Fig. 5c). The analysis showed that the increase in gene expression was inversely correlated with the reduction in DNA methylation level (decreased ∆β) on the same genes, suggesting that DNA hypomethylation may be involved in regulation of gene expression for the TLR pathway in SLE.

Aberrant mRNA expression in SLE with LN

In order to identify the genes which may play a role in the development of LN, we compared the gene expression profiles in PBMC between SLE LN+ and SLE LN− patients. From the 552 common upregulated and 550 common downregulated genes in SLE, we only identified three genes (MX1, G protein-coupled receptor 84 (GPR84), E2F transcription factor 2 (E2F2)), which were significantly increased in SLE LN+ compared with both NC and SLE LN− (p < 0.05, Fig. 6). We then checked the methylation status of the three upregulated genes and no significant methylation change was identified on these genes between SLE LN+ and SLE LN− patients.

Fig. 6
figure 6

Upregulated genes in PBMC of SLE LN+ patients. Three genes were identified to be significantly upregulated in SLE LN+ compared with SLE LN− and NC (p < 0.05). Expression value represented by normalized intensity in SLE LN+, SLE LN−, and NC.*Indicates p < 0.05. NC normal controls, SLE systemic lupus erythematosus, SLE LN + SLE with lupus nephritis, SLE LN − SLE without lupus nephritis

Full size image

Aberrant level of serum cytokine/chemokine in SLE and their transcriptional regulation

Cytokines are immunoregulatory mediators which play important roles in the pathogenesis of SLE. In an effort to understand the transcriptional regulation on serum cytokines, we performed quantitative analysis on 27 cytokines/chemokines in the sera of SLE and NC using a bead-based multiplex assay. Fifteen of the 27 serum factors, including IL-2, IL-4, IL-5, IL-6, IL-10, IL-12p70, IL-15, IL-17A, bFGF, MCP-1, IP-10, G-CSF, GM-CSF, IL-1RA and TNF-α were significantly higher in SLE LN− and SLE LN+, compared with NC (Fig. 7a, Additional file 7: Table S6). However, two cytokines, PDGF-BB and RANTES, were significantly lower in the two SLE groups compared with NC (Fig. 7a). There was no significant difference on the serum levels of the cytokines between SLE LN+ patients and their SLE LN− counterparts (Fig. 7a, Additional file 7: Table S6). At the transcription level, we found that IL-15 mRNA was elevated in SLE, and RANTES (chemokine (C-C motif) ligand 5 (CCL5)) mRNA was decreased in SLE, compared with NC (Fig. 7b), indicating that the serum level of these two cytokines corresponded with their mRNA expression in PBMC. We further checked the DNA methylation status on IL-15 and CCL5, and identified six of 11 methylation sites (cg01856970, cg02934500, cg25546588, cg16655388, cg26269613, and cg20060523) on IL-15, which are located at 5′UTR, TSS1500, and TSS 200 regions, were hypomethylated in both SLE LN+ and SLE LN− patients in relation to NC (Fig. 7c). We also found three CpG sites (cg00447324, cg08656816, and cg19411729) at the 5′UTR region of CCL5 (RANTES) gene were hypermethylated in SLE compared with controls (Fig. 7c).

Fig. 7
figure 7

Quantification of 27 cytokines and chemokines in sera of SLE patients and controls using bead-based multiplex assays. a Seventeen cytokines showed significant difference in SLE patients compared with NC, 15 were increased and two were decreased in both SLE LN− and SLE LN+ compared with NC (* p < 0.05). b mRNA expression level in PBMC for the cytokines which showed significant difference in sera of SLE patients. IL-15 mRNA was significantly upregulated and CCL5 mRNA were downregulated in PBMC of SLE patients compared with NC (* p < 0.05). c Six CpG sites in the IL-15 gene were hypomethylated and three CpG sites in CCL5 (RANTES) gene were hypomethylated in PBMC of SLE compared with NC. (* p < 0.05). NC normal controls, SLE systemic lupus erythematosus, SLE LN + SLE with lupus nephritis, SLE LN − SLE without lupus nephritis

Full size image

Discussion

Multiple studies on a variety of cell types have demonstrated altered gene expression and DNA methylation patterns in autoimmune diseases, providing important new directions toward understanding the pathogenesis of immune-mediated diseases [27, 28]. A number of studies have compared gene expression in PBMC from patients with SLE versus healthy individuals or other autoimmune conditions, which provide substantial advantages in the search for diagnostic and prognostic biomarkers in autoimmune disease [29, 30]. But LN remains a major challenge and one of the most severe manifestations of SLE. Early diagnosis and intervention are essential for favorable patient outcomes [31]. It would be of great help to identify biomarkers so that we could assess LN activity or kidney damage, and predict an imminent exacerbation of renal disease accurately [2].

In this study, we performed an unbiased genome-wide gene transcription and DNA methylation study in PBMC from SLE patients with LN and without LN and compared with matched normal controls. By cross comparison of different SLE groups with NC, we identified a group of genes (552 up- and 550 downregulated) which are consistently differentially expressed in SLE as compared with controls. DNA methylation analysis on the whole-genome level distinguished a large number of CpG sites which were differentially methylated (1813 hypermethylated and 3785 hypomethylated sites) in SLE compared to NC. With integrated analysis of the gene transcription and DNA methylation data, we found that 60.5 % of the upregulated genes contain hypomethylated CpG sites and 87.1 % downregulated genes harbor hypermethylated CpG sites in their regulatory regions of the genome, suggesting that DNA methylation may have regulatory effect on most of the differentially expressed genes in SLE. Profiling mRNA expression and DNA methylation in the same set of SLE and NC samples enable us to investigate the complexity of gene dysregulation in the context of this disease.

Further analysis was performed to elucidate the potential molecular pathways modulated by the aberrantly expressed and methylated genes in SLE. The current analysis was mainly focused on the common upregulated genes with DNA hypomethylation. To our expectation, the upregulated genes in SLE were significantly enriched in immune regulatory pathways such as IFN signaling pathway, TLR signaling pathway, and proinflammatory signaling pathways, which have been reported to be related with pathogenesis in SLE and other autoimmune disorders [32–34]. It is well known that the overexpression of type I IFN-related genes in peripheral blood immune cells is a major SLE signature [35]. However, the pathogenic role of IFNs and their exact regulation in the development of SLE is still under debate. By exploring the globe gene transcription and DNA methylation in PBMC of SLE with different clinical manifestations (with or without LN), our study provided further evidence that the overexpression of a cluster of IFN genes was involved in the activation of the IFN signaling pathway. In addition, we noticed a stepwise increased expression of nine IFN genes from NC to SLE LN− to SLE LN+, suggesting a pathogenic role of the IFNs in the development of SLE and LN, even though the difference between SLE LN− and SLE LN+ was not statistically significant for all IFNs except for MX1. Interestingly, almost all the overexpressed IFN genes in SLE were associated with multiple hypomethylated CpG sites in the regulatory region of the genes, indicating epigenetic modification may play a major role in the regulation of IFN gene expression. Previous studies have reported the hypomethylation of IFN genes in the CD4+ T cells in SLE [13]. Our study detected the hypomethylated CpG sites in association with upregulation of IFN genes in PBMC of SLE, suggesting the hypomethylation may be a general regulating mechanism in all immune cells in PBMC, including T, B, macrophages, and dendritic cells. Further studies need to be performed on different cell types and different stages of the diseases to elucidate the exact molecular mechanisms of IFN gene regulation.

Toll-like receptor pathway is another molecular pathway identified in our study to be enriched with dysregulated genes in SLE. Seven upregulated genes in both SLE LN− and SLE LN+, including TLR1, TLR4, TLR7, TLR8, CD14, EIF2AK2 and IRAK3, were involved in the regulation of TLR pathways. TLR regulation plays an essential role in the pathogenesis of SLE [36]. It is believed that abnormal activation of TLR7/TLR8 or TLR9 by self RNA or DNA in antigen-presenting cells may cause breach of immune tolerance and results in cytokine, IFN-I, and autoantibody production in T, B, and myeloid cells [37–39]. TLR4 have also been shown to induce TNF-α production in monocytes of SLE patients [40, 41]. Substantial studies have been focused on the TLR signaling pathway for exploration of therapeutic target using adjuvants or inhibitors [42, 43]. Our data revealed the aberrant DNA methylation of TLR and its signaling pathway in SLE, implying a novel therapeutic strategy by intervention of methylation in SLE.

Serum cytokines/chemokines are important biomarkers reflecting immune cell activation and inflammatory reaction during the development of SLE [44]. In this study, we measured the serum levels of 27 pre-inflammatory cytokines/chemokines in SLE patients with or without lupus nephritis and compared them with healthy controls. Our analysis identified 15 cytokines/chemokines which significantly increased in the sera of both SLE LN− and SLE LN+ patients, as compared with NC. The increased levels of Th1 cytokines (IL-2, IL-12, IL-15, TNF-a), and Th2 cytokines (IL-4, IL-5, IL-6, IL-10,) reflecting the activation of both Th1 and Th2 T helper cell subsets, which is in agreement with previous studies showing that IL-12, IL-10, IL-6 and TNF-α were correlated with active SLE [45–48]. We also observed the increased level of IL-17, which is secreted by Th17 cells and involved in T cell activation and autoantibodies in SLE [49–51]. It is not surprising that we could not find significant difference on any of the serum cytokines between SLE LN+ and SLE LN− in our sample set. Although several cytokines, including MCP-1, IP-10 and IL-6, have been reported to have increased expression in the kidney and urine of SLE patients with active LN compared with no renal SLE, the serum level were not different between SLE patients with or without LN [52], implying that the serum levels of inflammatory cytokines/chemokines are not reflecting the LN activity in SLE patients.

Our data also showed that RANTES (CCL5) and platelet-derived growth factor-BB (PDGF-BB) were significantly lower in both SLE LN− and SLE LN+ compared with NC. RANTES, as a proinflammatory chemokine, has been extensively studied and reported to be increased in SLE patients [53, 54]. The discrepancy of our data with previous studies may be due to sample variation and different method used. PDGF-BB is a platelet-derived factor involved in cell proliferation. Its role in SLE has not been reported before.

The regulation of the serum cytokines/chemokines in SLE on the genome level was not fully understood. By analyzing the mRNA expression and the DNA methylation of the cytokine/chemokine genes, we were able to find out if the increased serum level of cytokines in SLE is associated with their mRNA transcription and DNA methylation in the PBMC of the patients. To our surprise, we only found two out of the 17 differentially expressed cytokines, the IL-15 and RANTES, in which the mRNA expression in PBMC correspond with their serum levels. The discrepancy between serum protein level and mRNA for most of the cytokines/chemokines may be due to the cytokine gene expression regulation is cell-type specific and the differential expression may not be detected in PBMC, which is a mixture of different cell types. On the other hand, the increased mRNA transcription on IL-15 in SLE PBMC suggesting the transcription of this cytokine may be activated in multiple cell types. Indeed, IL-15 can be produced by macrophage/monocytes, activated NK cells, CD4+ T cells, and B cells, and involved in the pathogenesis of SLE and other autoimmune diseases [55–57]. Furthermore, the hypomethylation on several CpG sites in the regulatory region of this gene suggested a strong epigenetic regulation of IL-15 in SLE. In addition, our study showed that the lower level of RANTES in serum is correlated with lower mRNA expression in PBMC in SLE patients. Interestingly, the hypermethylated CpG sites at the 5′UTR region of CCL5 (RANTES) gene in SLE patients also indicated a suppressive effect of epigenetic regulation on this gene. This data is currently under validation in a larger sample cohort.

LN is one of the most devastating complications of SLE, but the molecular basis for this disease is not yet clear. Gene expression profiling on PBMC and different subsets of blood cells have identified candidate genes that might contribute to the pathogenesis of SLE and its complications [7, 58, 59]. However, the role of the differentially expressed genes identified from peripheral blood of SLE patients in the pathogenesis of LN is still under debate. By comparing the gene expression profiles in PBMC from SLE LN− and SLE LN+ patients, we identified three genes which are significantly elevated in SLE LN+ patients, namely, MX1, GPR84 and E2F2. MX1 is one of the IFN-I-inducible genes which have been extensively studied in SLE. Increased MX1 gene expression has been detected in both PBMC and renal intrinsic cells of lupus nephritis patients suggesting its role in the pathogenesis of LN [60, 61]. Our data consolidated previous finding by detecting the upregulation of MX1 on transcription in PBMC and in serum in SLE LN+ patients. Furthermore, by analyzing DNA methylation of this gene, we postulated that the higher expression of MX1 in SLE may be regulated by hypomethylation of the CpGs in the regulatory region.

E2F2 and GPR84 are novel biomarkers which have not been previously reported to be associated with SLE or LN. E2F2 is a member of E2F transcription factor family. As a cell cycle activator, E2F2 can regulate the gene expression that govern the cell cycle and DNA synthesis, and thus cell proliferation [62, 63]. There was a study that showed that E2F2 can promote the expression of microRNA let-7a by binding to its promoter and let-7a has been shown to contribute to hyperplasia and the proinflammatory response in SLE [64]. Therefore, the increased expression of E2F2 may associate with the inflammation in SLE and LN. GPR84 is a G protein-coupled receptor for sensing free fatty acid and active inflammatory reaction. Activation of GPR84 by its ligand could elicit chemotaxis of polymorphonuclear leukocytes (PMNs) and macrophages and amplified lipopolysaccharide (LPS)-stimulated production of the proinflammatory cytokine IL-8 from PMNs and TNF-α from macrophages [65]. Further analysis on the expression of these two genes in different cell types and the molecular pathways is warranted to elucidate their role in the pathogenesis of SLE and LN.

Although our current analysis was primarily focused on the genes which were upregulated and hypomethylated in SLE, the identification of a large amount of downregulated genes which were hypermethylated in the SLE genome by globe gene expression and DNA methylation analysis implying that downregulation or silencing of gene expression in PBMC may play an equally important role in SLE and LN [66–70]. Previous studies have demonstrated that the downregulation of IL-2 [35, 67, 71], IL-17 [66], and Notch-1 [69] in SLE T lymphocytes was due to the binding cAMP-responsive element modulator (CREM)-α to the promoter of the genes followed by enhanced CpG DNA methylation. We therefore examined methylation status on the CpG sites of these genes in our study. By comparing the mean β value between NC and SLE (both with LN and without LN involvement), we found that four out of 70 sites on Notch-1 gene, two out of ten sites on IL-17A gene and one out of two sites on IL-2 gene were hypermethylated in SLE (p < 0.05) (data not shown). Further analysis on the SLE downregulation gene data set is going on with the expectation to find out more dysregulated pathways modulated by DNA hypermethylation.

There are some limitations in the current study. First, the gene expression and DNA methylation profiling were only performed in PBMC. Given the highly complex nature of this disease, different cell subsets in peripheral blood may have different transcription and methylation profiles during the development of the diseases. However, previous studies using different cell types have identified specific gene signatures from different cell types. These cell type-specific genes allowed a correct classification of PBMC independent from their heterogenic cellular composition. Therefore, based on the expression of cell type-specific signature genes, we should be able to normalize the expression profile of PBMC from all individuals [72]. On the other hand, the advantage of using PBMC includes easy accessibility and less blood volume requirement. The result from PBMC has reflected a general picture of profiling alteration in whole blood. It is therefore ideal for biomarker screening [73]. The second limitation is that the sample size in this study is relatively small. The statistical power may not be strong enough to draw conclusions on some of the differential genes and pathways. The current data will need to be further validated in a larger sample cohort. Third, functional analysis will need to be conducted to confirm the genetic and epigenetic interactions on the dysregulation of genes and pathways in SLE.

Conclusions

In conclusion, we integrated the gene expression, DNA methylation, and cytokine profiling data in PBMC of SLE patients with and without lupus nephritis. Our study has demonstrated that significant number of differential genes in SLE were involved in IFN, TLR signaling pathways, and inflammatory cytokines. The enrichment of differential genes has been associated with aberrant DNA methylation, which may be relevant to the pathogenesis of SLE. Our observations have laid the groundwork for further diagnostic and mechanistic studies of SLE and LN.

Abbreviations

ACR, American College of Rheumatology; bFGF, basic fibroblast growth factor; CCL5, chemokine (C-C motif) ligand 5; CD14, cluster of differentiation 14; E2F2, E2F transcription factor 2; EIF2AK2, eukaryotic translation initiation factor 2-alpha kinase 2; G-CSF, granulocyte colony-stimulating factor; GM-CSF, granulocyte-macrophage colony-stimulating factor; GPR84, G protein-coupled receptor 84; HPF, high-power field; IFN, interferon; IL, interleukin; IP-10, IFN-gamma-inducible protein 10; IPA, Ingenuity Pathway Analysis; IRAK, interleukin-1 receptor-associated kinase; LN, lupus nephritis; MCP-1, monocyte chemoattractant protein; MIP-1, macrophage inflammatory protein 1; MX1, myxovirus (influenza) resistance 1; NC, normal controls; PBMC, peripheral blood mononuclear cell; PDGF, platelet-derived growth factor-BB; RANTES, regulated on activation, normal T cell expressed and secreted; SLE, systemic lupus erythematosus; SLEDAI, SLE Disease Activity Index; SLE LN+, SLE with lupus nephritis; SLE LN−, SLE without lupus nephritis; TLR, toll-like receptor; TNF-α, tumor necrosis factor alpha

References

  1. Frieri M. Mechanisms of disease for the clinician: systemic lupus erythematosus. Ann Allergy Asthma Immunol. 2013;110(4):228–32.

    Article  PubMed  Google Scholar 

  2. Hahn BH, McMahon MA, Wilkinson A, Wallace WD, Daikh DI, Fitzgerald JD, Karpouzas GA, Merrill JT, Wallace DJ, Yazdany J, et al. American College of Rheumatology guidelines for screening, treatment, and management of lupus nephritis. Arthritis Care Res. 2012;64(6):797–808.

  3. Koutsokeras T, Healy T. Systemic lupus erythematosus and lupus nephritis. Nat Rev Drug Discov. 2014;13(3):173–4.

    Article  CAS  PubMed  Google Scholar 

  4. Lorenz G, Desai J, Anders HJ. Lupus nephritis: update on mechanisms of systemic autoimmunity and kidney immunopathology. Curr Opin Nephrol Hypertens. 2014;23(3):211–7.

    Article  CAS  PubMed  Google Scholar 

  5. Castrejon I, Tani C, Jolly M, Huang A, Mosca M. Indices to assess patients with systemic lupus erythematosus in clinical trials, long-term observational studies, and clinical care. Clin Exp Rheumatol. 2014;32(5 Suppl 85):S-85–95.

    PubMed  Google Scholar 

  6. Weening JJ, D'Agati VD, Schwartz MM, Seshan SV, Alpers CE, Appel GB, Balow JE, Bruijn JA, Cook T, Ferrario F, et al. The classification of glomerulonephritis in systemic lupus erythematosus revisited. Kidney Int. 2004;65(2):521–30.

  7. Li Y, Fang X, Li QZ. Biomarker profiling for lupus nephritis. Genomics Proteomics Bioinformatics. 2013;11(3):158–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Borchers AT, Leibushor N, Naguwa SM, Cheema GS, Shoenfeld Y, Gershwin ME. Lupus nephritis: a critical review. Autoimmun Rev. 2012;12(2):174–94.

    Article  PubMed  Google Scholar 

  9. Lisnevskaia L, Murphy G, Isenberg D. Systemic lupus erythematosus. Lancet. 2014;384(9957):1878–88.

    Article  PubMed  Google Scholar 

  10. Raj P, Rai E, Song R, Khan S, Wakeland BE, Viswanathan K, Arana C, Liang C, Zhang B, Dozmorov I, et al. Regulatory polymorphisms modulate the expression of HLA class II molecules and promote autoimmunity. Elife. 2016;5. doi: 10.7554/eLife.12089

  11. Picascia A, Grimaldi V, Pignalosa O, De Pascale MR, Schiano C, Napoli C. Epigenetic control of autoimmune diseases: from bench to bedside. Clin Immunol. 2015;157(1):1–15.

    Article  CAS  PubMed  Google Scholar 

  12. Miao CG, Yang JT, Yang YY, Du CL, Huang C, Huang Y, Zhang L, Lv XW, Jin Y, Li J. Critical role of DNA methylation in the pathogenesis of systemic lupus erythematosus: new advances and future challenges. Lupus. 2014;23(8):730–42.

  13. Absher DM, Li X, Waite LL, Gibson A, Roberts K, Edberg J, Chatham WW, Kimberly RP. Genome-wide DNA methylation analysis of systemic lupus erythematosus reveals persistent hypomethylation of interferon genes and compositional changes to CD4+ T-cell populations. PLoS Genet. 2013;9(8), e1003678.

  14. Coit P, Jeffries M, Altorok N, Dozmorov MG, Koelsch KA, Wren JD, Merrill JT, McCune WJ, Sawalha AH. Genome-wide DNA methylation study suggests epigenetic accessibility and transcriptional poising of interferon-regulated genes in naive CD4+ T cells from lupus patients. J Autoimmun. 2013;43:78–84.

  15. Jeffries MA, Dozmorov M, Tang Y, Merrill JT, Wren JD, Sawalha AH. Genome-wide DNA methylation patterns in CD4+ T cells from patients with systemic lupus erythematosus. Epigenetics. 2011;6(5):593–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Zhao M, Liu S, Luo S, Wu H, Tang M, Cheng W, Zhang Q, Zhang P, Yu X, Xia Y, et al. DNA methylation and mRNA and microRNA expression of SLE CD4+ T cells correlate with disease phenotype. J Autoimmun. 2014;54:127–36.

  17. Zhang Y, Zhao M, Sawalha AH, Richardson B, Lu Q. Impaired DNA methylation and its mechanisms in CD4(+)T cells of systemic lupus erythematosus. J Autoimmun. 2013;41:92–9.

    Article  PubMed  Google Scholar 

  18. Lu Q, Kaplan M, Ray D, Zacharek S, Gutsch D, Richardson B. Demethylation of ITGAL (CD11a) regulatory sequences in systemic lupus erythematosus. Arthritis Rheum. 2002;46(5):1282–91.

    Article  CAS  PubMed  Google Scholar 

  19. Lu Q, Wu A, Richardson BC. Demethylation of the same promoter sequence increases CD70 expression in lupus T cells and T cells treated with lupus-inducing drugs. J Immunol. 2005;174(10):6212–9.

    Article  CAS  PubMed  Google Scholar 

  20. Kaplan MJ, Lu Q, Wu A, Attwood J, Richardson B. Demethylation of promoter regulatory elements contributes to perforin overexpression in CD4+ lupus T cells. J Immunol. 2004;172(6):3652–61.

    Article  CAS  PubMed  Google Scholar 

  21. Lu Q, Wu A, Tesmer L, Ray D, Yousif N, Richardson B. Demethylation of CD40LG on the inactive X in T cells from women with lupus. J Immunol. 2007;179(9):6352–8.

    Article  CAS  PubMed  Google Scholar 

  22. Ramos PS, Brown EE, Kimberly RP, Langefeld CD. Genetic factors predisposing to systemic lupus erythematosus and lupus nephritis. Semin Nephrol. 2010;30(2):164–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Tan EM, Cohen AS, Fries JF, Masi AT, McShane DJ, Rothfield NF, Schaller JG, Talal N, Winchester RJ. The 1982 revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum. 1982;25(11):1271–7.

  24. Dooley MA, Aranow C, Ginzler EM. Review of ACR renal criteria in systemic lupus erythematosus. Lupus. 2004;13(11):857–60.

    Article  CAS  PubMed  Google Scholar 

  25. Assenov Y, Muller F, Lutsik P, Walter J, Lengauer T, Bock C. Comprehensive analysis of DNA methylation data with RnBeads. Nat Methods. 2014;11(11):1138–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98(9):5116–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Somers EC, Richardson BC. Environmental exposures, epigenetic changes and the risk of lupus. Lupus. 2014;23(6):568–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Lee HS, Bae SC. What can we learn from genetic studies of systemic lupus erythematosus? Implications of genetic heterogeneity among populations in SLE. Lupus. 2010;19(12):1452–9.

    Article  CAS  PubMed  Google Scholar 

  29. Nikpour M, Dempsey AA, Urowitz MB, Gladman DD, Barnes DA. Association of a gene expression profile from whole blood with disease activity in systemic lupus erythaematosus. Ann Rheum Dis. 2008;67(8):1069–75.

    Article  CAS  PubMed  Google Scholar 

  30. Baechler EC, Batliwalla FM, Karypis G, Gaffney PM, Ortmann WA, Espe KJ, Shark KB, Grande WJ, Hughes KM, Kapur V, et al. Interferon-inducible gene expression signature in peripheral blood cells of patients with severe lupus. Proc Natl Acad Sci U S A. 2003;100(5):2610–5.

  31. Anders HJ, Appel GB. Lupus nephritis: Implications of the new ACR lupus nephritis guidelines. Nat Rev Nephrol. 2012;8(9):500–1.

    Article  PubMed  Google Scholar 

  32. Rönnblom L, Eloranta ML. The interferon signature in autoimmune diseases. Curr Opin Rheumatol. 2013;25(2):248–53.

    Article  PubMed  Google Scholar 

  33. Ma K, Li J, Fang Y, Lu L. Roles of B cell-intrinsic TLR signals in systemic lupus erythematosus. Int J Mol Sci. 2015;16(6):13084–105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Oates JC, Gilkeson GS. The biology of nitric oxide and other reactive intermediates in systemic lupus erythematosus. Clin Immunol. 2006;121(3):243–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Hagberg N, Rönnblom L. Systemic lupus erythematosus--a disease with a dysregulated type I interferon system. Scand J Immunol. 2015;82(3):199–207.

    Article  CAS  PubMed  Google Scholar 

  36. Guggino G, Giardina AR, Ciccia F, Triolo G, Dieli F, Sireci G. Are Toll-like receptors and decoy receptors involved in the immunopathogenesis of systemic lupus erythematosus and lupus-like syndromes? Clin Dev Immunol. 2012;2012:135932.

    Article  PubMed  Google Scholar 

  37. Guo Y, Chai Q, Zhao Y, Li P, Qiao J, Huang J. Increased activation of toll-like receptors-7 and -8 of peripheral blood mononuclear cells and upregulated serum cytokines in patients with pediatric systemic lupus erythematosus. Int J Clin Exp Med. 2015;8(11):20472–80.

    PubMed  PubMed Central  Google Scholar 

  38. Chen JQ, Szodoray P, Zeher M. Toll-like receptor pathways in autoimmune diseases. Clin Rev Allergy Immunol. 2016;50(1):1–17.

    Article  CAS  PubMed  Google Scholar 

  39. Mohan C, Putterman C. Genetics and pathogenesis of systemic lupus erythematosus and lupus nephritis. Nat Rev Nephrol. 2015;11(6):329–41.

    Article  CAS  PubMed  Google Scholar 

  40. Marques CP, Maor Y, de Andrade MS, Rodrigues VP, Benatti BB. Possible evidence of systemic lupus erythematosus and periodontal disease association mediated by Toll-like receptors 2 and 4. Clin Exp Immunol. 2016;183(2):187–92.

    Article  CAS  PubMed  Google Scholar 

  41. Carvalheiro T, Gomes D, Pinto LA, Ines L, Lopes A, Henriques A, Pedreiro S, Martinho A, Trindade H, Young HA, et al. Sera from patients with active systemic lupus erythematosus patients enhance the toll-like receptor 4 response in monocyte subsets. J Inflamm. 2015;12:38.

  42. Hennessy EJ, Parker AE, O'Neill LA. Targeting Toll-like receptors: emerging therapeutics? Nat Rev Drug Discov. 2010;9(4):293–307.

    Article  CAS  PubMed  Google Scholar 

  43. Li J, Wang X, Zhang F, Yin H. Toll-like receptors as therapeutic targets for autoimmune connective tissue diseases. Pharmacol Ther. 2013;138(3):441–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Yu SL, Kuan WP, Wong CK, Li EK, Tam LS. Immunopathological roles of cytokines, chemokines, signaling molecules, and pattern-recognition receptors in systemic lupus erythematosus. Clin Dev Immunol. 2012;2012:715190.

    Article  PubMed  PubMed Central  Google Scholar 

  45. McCarthy EM, Smith S, Lee RZ, Cunnane G, Doran MF, Donnelly S, Howard D, O'Connell P, Kearns G, Ni Gabhann J, et al. The association of cytokines with disease activity and damage scores in systemic lupus erythematosus patients. Rheumatology. 2014;53(9):1586–94.

  46. Qiu F, Song L, Yang N, Li X. Glucocorticoid downregulates expression of IL-12 family cytokines in systemic lupus erythematosus patients. Lupus. 2013;22(10):1011–6.

    Article  CAS  PubMed  Google Scholar 

  47. Yang X, Sun B, Wang H, Yin C, Wang X, Ji X. Increased serum IL-10 in lupus patients promotes apoptosis of T cell subsets via the caspase 8 pathway initiated by Fas signaling. J Biomed Res. 2015;29(3):232–40.

    PubMed  Google Scholar 

  48. Umare V, Pradhan V, Nadkar M, Rajadhyaksha A, Patwardhan M, Ghosh KK, Nadkarni AH. Effect of proinflammatory cytokines (IL-6, TNF-alpha, and IL-1beta) on clinical manifestations in Indian SLE patients. Mediators Inflamm. 2014;2014:385297.

  49. Dong G, Ye R, Shi W, Liu S, Wang T, Yang X, Yang N, Yu X. IL-17 induces autoantibody overproduction and peripheral blood mononuclear cell overexpression of IL-6 in lupus nephritis patients. Chin Med J (Engl). 2003;116(4):543–8.

  50. Chang SH, Dong C. Signaling of interleukin-17 family cytokines in immunity and inflammation. Cell Signal. 2011;23(7):1069–75.

    Article  CAS  PubMed  Google Scholar 

  51. Abdel Galil SM, Ezzeldin N, El-Boshy ME. The role of serum IL-17 and IL-6 as biomarkers of disease activity and predictors of remission in patients with lupus nephritis. Cytokine. 2015;76(2):280–7.

    Article  CAS  PubMed  Google Scholar 

  52. Misra R, Gupta R. Biomarkers in lupus nephritis. Int J Rheum Dis. 2015;18(2):219–32.

    Article  PubMed  Google Scholar 

  53. Kulkarni O, Anders HJ. Chemokines in lupus nephritis. Front Biosci. 2008;13:3312–20.

    Article  CAS  PubMed  Google Scholar 

  54. Lu MM, Wang J, Pan HF, Chen GM, Li J, Cen H, Feng CC, Ye DQ. Increased serum RANTES in patients with systemic lupus erythematosus. Rheumatol Int. 2012;32(5):1231–3.

  55. Tagaya Y, Bamford RN, DeFilippis AP, Waldmann TA. IL-15: a pleiotropic cytokine with diverse receptor/signaling pathways whose expression is controlled at multiple levels. Immunity. 1996;4(4):329–36.

    Article  CAS  PubMed  Google Scholar 

  56. Aringer M, Stummvoll GH, Steiner G, Koller M, Steiner CW, Hofler E, Hiesberger H, Smolen JS, Graninger WB. Serum interleukin-15 is elevated in systemic lupus erythematosus. Rheumatology. 2001;40(8):876–81.

  57. Baranda L, de la Fuente H, Layseca-Espinosa E, Portales-Perez D, Nino-Moreno P, Valencia-Pacheco G, Abud-Mendoza C, Alcocer-Varela J, Gonzalez-Amaro R. IL-15 and IL-15R in leucocytes from patients with systemic lupus erythematosus. Rheumatology. 2005;44(12):1507–13.

  58. Frangou EA, Bertsias GK, Boumpas DT. Gene expression and regulation in systemic lupus erythematosus. Eur J Clin Invest. 2013;43(10):1084–96.

    Article  CAS  PubMed  Google Scholar 

  59. Becker AM, Dao KH, Han BK, Kornu R, Lakhanpal S, Mobley AB, Li QZ, et al. SLE peripheral blood B cell, T cell and myeloid cell transcriptomes display unique profiles and each subset contributes to the interferon signature. PLoS One. 2013;8(6), e67003.

  60. Feng X, Wu H, Grossman JM, Hanvivadhanakul P, FitzGerald JD, Park GS, et al. Association of increased interferon-inducible gene expression with disease activity and lupus nephritis in patients with systemic lupus erythematosus. Arthritis Rheum. 2006;54(9):2951–62.

    Article  CAS  PubMed  Google Scholar 

  61. Watanabe S, Imaizumi T, Tsuruga K, Aizawa T, Ito T. Glomerular expression of myxovirus resistance protein 1 in human mesangial cells: possible activation of innate immunity in the pathogenesis of lupus nephritis. Nephrology. 2013;18(12):833–7.

    Article  CAS  PubMed  Google Scholar 

  62. Chen L, Yu JH, Lu ZH, Zhang W. E2F2 induction in related to cell proliferation and poor prognosis in non-small cell lung carcinoma. Int J Clin Exp Pathol. 2015;8(9):10545–54.

    PubMed  PubMed Central  Google Scholar 

  63. Zhou S, Sun H, Zheng B, Li R, Zhang W. Cell cycle transcription factor E2F2 mediates non-stress temperature response of AtHSP70-4 in Arabidopsis. Biochem Biophys Res Commun. 2014;455(3-4):139–64.

    Article  CAS  PubMed  Google Scholar 

  64. Chafin CB, Regna NL, Caudell DL, Reilly CM. MicroRNA-let-7a promotes E2F-mediated cell proliferation and NFkB activation in vitro. Cell Mol Immunol. 2014;11(1):79–83.

    Article  CAS  PubMed  Google Scholar 

  65. Suzuki M, Takaishi S, Nagasaki M, Onozawa Y, Iino I, Maeda H. Medium-chain fatty acid-sensing receptor, GPR84, is a proinflammatory receptor. J Biol Chem. 2013;288(15):10684–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Hedrich CM, Rauen T, Kis-Toth K, Kyttaris VC, Tsokos GC. cAMP-responsive element modulator alpha (CREMalpha) suppresses IL-17F protein expression in T lymphocytes from patients with systemic lupus erythematosus (SLE). J Biol Chem. 2012;287(7):4715–25.

    Article  CAS  PubMed  Google Scholar 

  67. Hedrich CM, Crispin JC, Rauen T, Ioannidis C, Apostolidis SA, Lo MS, Kyttaris VC, Tsokos GC. cAMP response element modulator alpha controls IL2 and IL17A expression during CD4 lineage commitment and subset distribution in lupus. Proc Natl Acad Sci U S A. 2012;109(41):16606–11.

  68. Hedrich CM, Crispin JC, Tsokos GC. Epigenetic regulation of cytokine expression in systemic lupus erythematosus with special focus on T cells. Autoimmunity. 2014;47(4):234–41.

    Article  CAS  PubMed  Google Scholar 

  69. Rauen T, Grammatikos AP, Hedrich CM, Floege J, Tenbrock K, Ohl K, Kyttaris VC, Tsokos GC. cAMP-responsive element modulator alpha (CREMalpha) contributes to decreased Notch-1 expression in T cells from patients with active systemic lupus erythematosus (SLE). J Biol Chem. 2012;287(51):42525–32.

  70. Sunahori K, Nagpal K, Hedrich CM, Mizui M, Fitzgerald LM, Tsokos GC. The catalytic subunit of protein phosphatase 2A (PP2Ac) promotes DNA hypomethylation by suppressing the phosphorylated mitogen-activated protein kinase/extracellular signal-regulated kinase (ERK) kinase (MEK)/phosphorylated ERK/DNMT1 protein pathway in T-cells from controls and systemic lupus erythematosus patients. J Biol Chem. 2013;288(30):21936–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Hedrich CM, Rauen T, Tsokos GC. cAMP-responsive element modulator (CREM)alpha protein signaling mediates epigenetic remodeling of the human interleukin-2 gene: implications in systemic lupus erythematosus. J Biol Chem. 2011;286(50):43429–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Kyogoku C, Smiljanovic B, Grün JR, Biesen R, Schulte-Wrede U, Häupl T, et al. Cell-specific type I IFN signatures in autoimmunity and viral infection: what makes the difference? PLoS One. 2013;8(12), e83776.

    Article  PubMed  PubMed Central  Google Scholar 

  73. Gilbert M, Punaro M. Blood gene expression profiling in pediatric systemic lupus erythematosus and systemic juvenile idiopathic arthritis: from bench to bedside. Pediatr Rheumatol Online J. 2014;12:16.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank the participants and the staff at the Division of Rheumatology, Xiangya Hospital for making this study possible. We thank Yun Lian at Microarray Core, UT Southwestern Medical Center for the microarray data submission to the GEO database.

This work was supported by grants from the National Natural Science Foundation of China (81270852, 81373206, 81401357) and the Scientific Research Fund of Xiangya Hospital (2013Q10).

Funding

This work was supported by grants from National Natural Science Foundation of China (81270852, 81373206, 81401357) and Scientific Research Fund of Xiangya Hospital (2013Q10).

Availability of supporting data

The microarray data for gene expression and methylation has been submitted to the GEO database with accession number GSE81622.

Additional supplementary tables and figures were also submitted with the manuscript.

Authors’ contributions

HZ participated in study design, performed data analysis, and drafted the manuscript. WM participated in data analysis, data interpretation, and revised the manuscript. HL participated in sample collection, cytokine data acquisition, and drafted the manuscript. TC participated in methylation data acquisition, data analysis, and drafted the manuscript. SL participated in the gene expression data acquisition, data analysis, and drafted the manuscript. IR participated in real-time PCR and cytokine data acquisition, data analysis, and drafted the manuscript. XZ coordinated the sample collection, conceived the study, and revised the manuscript. QZL conceived the study, participated in the study design, and revised and finalized the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no conflict of interest relating to the conduct of this study or the publication of this manuscript.

Consent for publication

All authors read and approved the final manuscript and consent for publication in Arthritis Research and Therapy.

Ethics approval and consent to participate

This study was approved by the Institutional Review Board at Xiangya Hospital, Central South University of China (Changsha, Hunan, China). Informed consent was obtained from all participants included in studies.

Author information

Author notes

    Authors and Affiliations

    1. Department of Rheumatology, Xiangya Hospital, Central South University, 87 Xiangya Road, Changsha, Hunan, 410008, People’s Republic of China

      Honglin Zhu, Hui Luo, Xiaoxia Zuo & Quan-Zhen Li

    2. Department of Immunology, University of Texas Southwestern Medical Center, 5323 Harry Hines Blvd., Dallas, TX, 75390, USA

      Honglin Zhu, Wentao Mi, Hui Luo, Tao Chen, Shengxi Liu, Indu Raman & Quan-Zhen Li

    Authors
    1. Honglin Zhu
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Wentao Mi
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Hui Luo
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Tao Chen
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Shengxi Liu
      View author publications

      You can also search for this author in PubMed Google Scholar

    6. Indu Raman
      View author publications

      You can also search for this author in PubMed Google Scholar

    7. Xiaoxia Zuo
      View author publications

      You can also search for this author in PubMed Google Scholar

    8. Quan-Zhen Li
      View author publications

      You can also search for this author in PubMed Google Scholar

    Corresponding authors

    Correspondence to Xiaoxia Zuo or Quan-Zhen Li.

    Additional files

    Additional file 1: Table S1.

    Upregulated genes in PBMC of SLE patients compared with normal controls by microarray analysis. S1-1 showed the complete list of 552 upregulated genes in SLE. Three group comparisons in gene expression, SLE (LN− and LN+) vs. NC (S1-2), SLE LN− vs. NC (S1-3), and SLE LN+ vs. NC (S1-4), were performed to identify upregulated genes by fold change (FC) > 1.2 and p < 0.05 within each group. The 552 common upregulated genes in SLE were then selected from each of the three group comparisons. S1-5 showed gene expression of 552 common upregulated genes represented by normalized intensity in SLE. (XLSX 694 kb)

    Additional file 2: Table S2.

    Downregulated genes in PBMC of SLE patients compared with normal controls by microarray analysis. S2-1 showed the complete list of 550 downregulated genes in SLE. Three group comparisons, SLE (LN− and LN+) vs. NC (S2-2), SLE LN− vs. NC (S2-3), and SLE LN+ vs. NC (S2-4), were performed to identify downregulated genes by fold change (FC) > 1.2 and p < 0.05 within each group. The 550 common downregulated genes in SLE were then selected from the genes identified in each of the three group comparisons. S2-5 showed gene expression of 550 common downregulated genes represented by normalized intensity in SLE. (XLSX 710 kb)

    Additional file 3: Table S3.

    Hypermethylated CpG sites in PBMC of SLE patients comparing normal controls by genome-wide DNA methylation analysis. Three group comparisons, SLE LN+ vs. NC (S3-1), SLE LN− vs. NC (S3-2), and SLE (LN− and LN+) vs. NC (S3-3), were performed to identify hypermethylated DNA CpG sites by fold change (FC) > 1.2 and q value < 5% within each group. The 1813 common hypermethylated sites identified in each of the three group comparisons were then selected and listed in S3-4. (XLSX 814 kb)

    Additional file 4: Table S4.

    Hypomethylated CpG sites in PBMC of SLE patients comparing normal controls by genome-wide DNA methylation analysis. Three group comparisons, SLE LN+ vs. NC (S4-1), SLE LN− vs. NC (S4-2), and SLE (LN− and LN+) vs. NC (S4-3), were performed to identify hypomethylated CpG sites by fold change (FC) > 1.2 and q value < 5% within each comparison. The 3785 common hypomethylated sites identified in each of the three group comparisons were then selected and listed in S4-4. (XLSX 1843 kb)

    Additional file 5: Table S5.

    Hypomethylated or hypermethylated CpG sites associated with the commonly up- or downregulated genes in SLE. S5-1 showed a complete list of 1065 hypomethylated CpG sites that were associated with 333 upregulated genes in SLE, and the normalized β values of each site for NC and SLE samples. S5-2 showed 3024 hypermethylated CpG sites that were associated with 478 downregulated genes, and the normalized β values of each site for NC and SLE samples. (XLSX 2264 kb)

    Additional file 6: Figure S1.

    Validation on the expression of five IFN-related genes in PBMC of a different cohort of SLE patients and controls by real-time PCR. The Ct value for IFI35, IFI44, MX1, STAT1, and IFIT1 was normalized with housekeeping gene GAPDH for each sample. Relative quantity (RQ) represents the relative level of gene expression in SLE and NC. (*Indicates significant difference between NC and SLE, p < 0.05). (PDF 496 kb)

    Additional file 7: Table S6.

    Profile 27 cytokine/chemokine from sera of SLE patients and normal controls. Serum cytokines were quantitatively measured using Bio-Plex Pro™ Human Cytokine 27-plex assay kits. The data represented the mean concentration for each cytokine in the sera of NC, SLE LN−, or SLE LN+ patients. (PDF 479 kb)

    Rights and permissions

    Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

    Reprints and permissions

    About this article

    Check for updates. Verify currency and authenticity via CrossMark

    Cite this article

    Zhu, H., Mi, W., Luo, H. et al. Whole-genome transcription and DNA methylation analysis of peripheral blood mononuclear cells identified aberrant gene regulation pathways in systemic lupus erythematosus. Arthritis Res Ther 18, 162 (2016). https://doi.org/10.1186/s13075-016-1050-x

    Download citation

    • Received:

    • Accepted:

    • Published:

    • DOI: https://doi.org/10.1186/s13075-016-1050-x

    Keywords

    Arthritis Research & Therapy

    ISSN: 1478-6362

    Contact us