Whole-genome transcription and DNA methylation analysis of peripheral blood mononuclear cells identified aberrant gene regulation pathways in systemic lupus erythematosus
- Honglin Zhu†1, 2,
- Wentao Mi†2,
- Hui Luo1, 2,
- Tao Chen2,
- Shengxi Liu2,
- Indu Raman2,
- Xiaoxia Zuo1Email author and
- Quan-Zhen Li1, 2Email author
© The Author(s). 2016
Received: 26 February 2016
Accepted: 20 June 2016
Published: 13 July 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.
KeywordsSystemic lupus erythematous Whole-genome transcription DNA methylation Multiplex cytokine assay Peripheral blood mononuclear cells Lupus nephritis
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
Clinical manifestation and laboratory data of SLE patients and normal controls
NC (n = 25)
SLE LN− (n = 15)
SLE LN+ (n = 15)
SLE LN− vs. NC
SLE LN+ vs. NC
SLE LN+ vs. SLE LN−
Age (mean + SD)
32.2 ± 9.04
28.2 ± 6.7
30.13 ± 7.68
Disease duration (months)
4.90 ± 4.18
3.33 ± 4.21
9.06 ± 4.06
14.93 ± 4.65
1997 ACR classification criteria for SLE
Pleuritis or pericarditis
Positive antinuclear antibody (ANA)
6.64 ± 0.70
6.45 ± 0.90
5.00 ± 0.69
4.74 ± 0.76
1.54 ± 0.21
1.22 ± 0.18
534.5 ± 280.1
293.2 ± 196.2
133.2 ± 147.4
88.02 ± 86.48
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
Aberrant DNA methylation in SLE
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
DNA hypomethylation related with the overexpression of IFN-related genes in SLE
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
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
Aberrant level of serum cytokine/chemokine in SLE and their transcriptional regulation
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
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.
Open AccessThis 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.
- Frieri M. Mechanisms of disease for the clinician: systemic lupus erythematosus. Ann Allergy Asthma Immunol. 2013;110(4):228–32.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- Koutsokeras T, Healy T. Systemic lupus erythematosus and lupus nephritis. Nat Rev Drug Discov. 2014;13(3):173–4.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.Google Scholar
- Li Y, Fang X, Li QZ. Biomarker profiling for lupus nephritis. Genomics Proteomics Bioinformatics. 2013;11(3):158–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Borchers AT, Leibushor N, Naguwa SM, Cheema GS, Shoenfeld Y, Gershwin ME. Lupus nephritis: a critical review. Autoimmun Rev. 2012;12(2):174–94.View ArticlePubMedGoogle Scholar
- Lisnevskaia L, Murphy G, Isenberg D. Systemic lupus erythematosus. Lancet. 2014;384(9957):1878–88.View ArticlePubMedGoogle Scholar
- 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: https://doi.org/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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- Dooley MA, Aranow C, Ginzler EM. Review of ACR renal criteria in systemic lupus erythematosus. Lupus. 2004;13(11):857–60.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Somers EC, Richardson BC. Environmental exposures, epigenetic changes and the risk of lupus. Lupus. 2014;23(6):568–76.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- Anders HJ, Appel GB. Lupus nephritis: Implications of the new ACR lupus nephritis guidelines. Nat Rev Nephrol. 2012;8(9):500–1.View ArticlePubMedGoogle Scholar
- Rönnblom L, Eloranta ML. The interferon signature in autoimmune diseases. Curr Opin Rheumatol. 2013;25(2):248–53.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Oates JC, Gilkeson GS. The biology of nitric oxide and other reactive intermediates in systemic lupus erythematosus. Clin Immunol. 2006;121(3):243–50.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- Chen JQ, Szodoray P, Zeher M. Toll-like receptor pathways in autoimmune diseases. Clin Rev Allergy Immunol. 2016;50(1):1–17.View ArticlePubMedGoogle Scholar
- Mohan C, Putterman C. Genetics and pathogenesis of systemic lupus erythematosus and lupus nephritis. Nat Rev Nephrol. 2015;11(6):329–41.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- Hennessy EJ, Parker AE, O'Neill LA. Targeting Toll-like receptors: emerging therapeutics? Nat Rev Drug Discov. 2010;9(4):293–307.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- Chang SH, Dong C. Signaling of interleukin-17 family cytokines in immunity and inflammation. Cell Signal. 2011;23(7):1069–75.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Misra R, Gupta R. Biomarkers in lupus nephritis. Int J Rheum Dis. 2015;18(2):219–32.View ArticlePubMedGoogle Scholar
- Kulkarni O, Anders HJ. Chemokines in lupus nephritis. Front Biosci. 2008;13:3312–20.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- Frangou EA, Bertsias GK, Boumpas DT. Gene expression and regulation in systemic lupus erythematosus. Eur J Clin Invest. 2013;43(10):1084–96.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar