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)
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).
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.
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.
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.
Systemic lupus erythematosus (SLE) is a potentially severe autoimmune disease characterized by intermittent episodes of increased disease activity that require treatment with immunosuppressive agents . Renal involvement affects more than 60 % of patients with SLE, and is a major contributor to morbidity and mortality . 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 . 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 . SLE should be closely monitored for renal manifestations, early diagnosis and treatment are essential for renal preservation . 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 .
Substantial research has helped define the pathogenic mechanisms of SLE renal manifestations . 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 . 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 . 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 . 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 .
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.
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 . 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 . 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 .
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) . The normalized methylation data were analyzed by significance analysis of microarrays (SAM, Stanford University) . 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.
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/).
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).
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.
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.
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).
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).
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).
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.
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).
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 . 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 .
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 . 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 . 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 . 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 . 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 , 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 . 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 . 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 , and Notch-1  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 . 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 . 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.
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.
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
Frieri M. Mechanisms of disease for the clinician: systemic lupus erythematosus. Ann Allergy Asthma Immunol. 2013;110(4):228–32.
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.
Koutsokeras T, Healy T. Systemic lupus erythematosus and lupus nephritis. Nat Rev Drug Discov. 2014;13(3):173–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.
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.
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.
Li Y, Fang X, Li QZ. Biomarker profiling for lupus nephritis. Genomics Proteomics Bioinformatics. 2013;11(3):158–65.
Borchers AT, Leibushor N, Naguwa SM, Cheema GS, Shoenfeld Y, Gershwin ME. Lupus nephritis: a critical review. Autoimmun Rev. 2012;12(2):174–94.
Lisnevskaia L, Murphy G, Isenberg D. Systemic lupus erythematosus. Lancet. 2014;384(9957):1878–88.
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
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Dooley MA, Aranow C, Ginzler EM. Review of ACR renal criteria in systemic lupus erythematosus. Lupus. 2004;13(11):857–60.
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.
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.
Somers EC, Richardson BC. Environmental exposures, epigenetic changes and the risk of lupus. Lupus. 2014;23(6):568–76.
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.
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.
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.
Anders HJ, Appel GB. Lupus nephritis: Implications of the new ACR lupus nephritis guidelines. Nat Rev Nephrol. 2012;8(9):500–1.
Rönnblom L, Eloranta ML. The interferon signature in autoimmune diseases. Curr Opin Rheumatol. 2013;25(2):248–53.
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.
Oates JC, Gilkeson GS. The biology of nitric oxide and other reactive intermediates in systemic lupus erythematosus. Clin Immunol. 2006;121(3):243–50.
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.
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.
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.
Chen JQ, Szodoray P, Zeher M. Toll-like receptor pathways in autoimmune diseases. Clin Rev Allergy Immunol. 2016;50(1):1–17.
Mohan C, Putterman C. Genetics and pathogenesis of systemic lupus erythematosus and lupus nephritis. Nat Rev Nephrol. 2015;11(6):329–41.
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.
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.
Hennessy EJ, Parker AE, O'Neill LA. Targeting Toll-like receptors: emerging therapeutics? Nat Rev Drug Discov. 2010;9(4):293–307.
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.
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.
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.
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.
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.
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.
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.
Chang SH, Dong C. Signaling of interleukin-17 family cytokines in immunity and inflammation. Cell Signal. 2011;23(7):1069–75.
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.
Misra R, Gupta R. Biomarkers in lupus nephritis. Int J Rheum Dis. 2015;18(2):219–32.
Kulkarni O, Anders HJ. Chemokines in lupus nephritis. Front Biosci. 2008;13:3312–20.
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.
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.
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.
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.
Frangou EA, Bertsias GK, Boumpas DT. Gene expression and regulation in systemic lupus erythematosus. Eur J Clin Invest. 2013;43(10):1084–96.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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).
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.
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.
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.
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)
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)
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)
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)
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)
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)
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)
About this article
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