You are here

  • Home
  • Archive
  • Volume 75, Issue 11
  • Genome-wide DNA methylation analysis in multiple tissues in primary Sjögren's syndrome reveals regulatory effects at interferon-induced genes

Article Text

Article menu
Basic and translational research
Extended report
Genome-wide DNA methylation analysis in multiple tissues in primary Sjögren's syndrome reveals regulatory effects at interferon-induced genes
  1. Juliana Imgenberg-Kreuz1,
  2. Johanna K Sandling1,2,
  3. Jonas Carlsson Almlöf1,
  4. Jessica Nordlund1,
  5. Linnea Signér2,
  6. Katrine Braekke Norheim3,
  7. Roald Omdal3,
  8. Lars Rönnblom2,
  9. Maija-Leena Eloranta2,
  10. Ann-Christine Syvänen1,
  11. Gunnel Nordmark2
  1. 1Molecular Medicine and Science for Life Laboratory, Department of Medical Sciences, Uppsala University, Uppsala, Sweden
  2. 2Rheumatology and Science for Life Laboratory, Department of Medical Sciences, Uppsala University, Uppsala, Sweden
  3. 3Clinical Immunology Unit, Department of Internal Medicine, Stavanger University Hospital, Stavanger, Norway
  1. Correspondence to Dr Gunnel Nordmark, Rheumatology and Science for Life Laboratory, Department of Medical Sciences, Uppsala University, Uppsala SE-751 85, Sweden; Gunnel.Nordmark{at}medsci.uu.se

Abstract

Objectives Increasing evidence suggests an epigenetic contribution to the pathogenesis of autoimmune diseases, including primary Sjögren's Syndrome (pSS). The aim of this study was to investigate the role of DNA methylation in pSS by analysing multiple tissues from patients and controls.

Methods Genome-wide DNA methylation profiles were generated using HumanMethylation450K BeadChips for whole blood, CD19+ B cells and minor salivary gland biopsies. Gene expression was analysed in CD19+ B cells by RNA-sequencing. Analysis of genetic regulatory effects on DNA methylation at known pSS risk loci was performed.

Results We identified prominent hypomethylation of interferon (IFN)-regulated genes in whole blood and CD19+ B cells, including at the genes MX1, IFI44L and PARP9, replicating previous reports in pSS, as well as identifying a large number of novel associations. Enrichment for genomic overlap with histone marks for enhancer and promoter regions was observed. We showed for the first time that hypomethylation of IFN-regulated genes in pSS B cells was associated with their increased expression. In minor salivary gland biopsies we observed hypomethylation of the IFN-induced gene OAS2. Pathway and disease analysis resulted in enrichment of antigen presentation, IFN signalling and lymphoproliferative disorders. Evidence for genetic control of methylation levels at known pSS risk loci was observed.

Conclusions Our study highlights the role of epigenetic regulation of IFN-induced genes in pSS where replication is needed for novel findings. The association with altered gene expression suggests a functional mechanism for differentially methylated CpG sites in pSS aetiology.

  • Sjøgren's Syndrome
  • B cells
  • Autoimmunity
  • Gene Polymorphism

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/

https://doi.org/10.1136/annrheumdis-2015-208659

Statistics from Altmetric.com

    Request Permissions

    If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

    Introduction

    Primary Sjögren's syndrome (pSS) is a chronic autoimmune disease characterised by inflammation of salivary and lacrimal glands. Systemic manifestations such as arthritis, pulmonary or renal involvement may also occur.1 The aetiology of pSS is multifactorial, where both genetic and environmental factors are thought to contribute to disease development. However, the molecular mechanisms underlying pSS remain largely elusive.2 Activation of the type I interferon (IFN) system with elevated expression of type I IFN-regulated genes in blood cells and minor salivary glands, a so-called IFN signature, has been demonstrated in pSS.3–7 B cell activation in pSS is reflected by autoantibody synthesis and an increased risk of non-Hodgkin's lymphoma, most commonly of the B cell type.8

    For pSS, to date, about 10 genetic risk loci have been identified through genome-wide association studies (GWAS).9 ,10 These variants only explain a limited proportion of the susceptibility to pSS, and the functional consequence of associated single nucleotide polymorphisms (SNPs) remains unclear for most loci.11 Increasing evidence suggests an epigenetic contribution to the pathogenesis of autoimmune diseases, including pSS.12 Epigenetic modifications constitute an additional layer of genomic regulation, and may serve as a dynamic link between genotype, environment and phenotype, for example by modulating gene expression. Methylation of the DNA base cytosine (5mC) can be studied in large sample sets with array-based methods such as the Illumina HumanMethylation450K array (HM450K), which allows for quantification of DNA methylation of 485 577 CpG sites across the human genome.13

    In pSS, so far, two studies with relatively small sample sizes have applied the HM450K array.14 ,15 Altorok et al14 report a number of differentially methylated CpG sites (DMCs) in naïve CD4+ T cells from 11 patients and 11 controls, with implications in immune responses and lymphocyte activation, including hypomethylation of IFN-regulated genes. Miceli-Richard et al15 found more prominent methylation changes in CD19+ B cells than in CD4+ T cells in 26 patients compared with 22 controls.

    In order to evaluate the role of DNA methylation in pSS in a comprehensive manner we performed an epigenome-wide association study (EWAS) in whole blood, CD19+ B-cells and minor salivary gland biopsy samples from patients and controls using the HM450K array. To further explore the functional role of the DMCs we intersected the disease-associated CpG sites with publicly available chromatin state data16 and performed gene expression analysis in CD19+ B cells. Finally, genetic regulation of methylation at pSS risk loci was investigated.

    Patients and methods

    For full details of methods see online supplementary text.

    Patients and controls

    A total of 108 Caucasian patients with pSS from the rheumatology clinic at the Uppsala University Hospital, Sweden, were included in the study, all fulfilling the American European Consensus Group (AECG) criteria.17 Whole blood was collected from 100 patients, CD19+ B cells from 24 patients and minor salivary gland biopsies from 15 patients. As control samples, DNA from whole blood from 400 healthy blood donors from the Uppsala Bioresource and CD19+ B cells from 47 donors were analysed.18 Only controls falling within the main European cluster in our previous study were included.19 Minor salivary gland biopsies obtained from 13 individuals that were examined for a possible pSS diagnosis, where the biopsies showed no inflammation and serology was negative for autoantibodies, served as control biopsies (table 1).

    View this table:
    Table 1

    Characteristics of patients with primary Sjögren's syndrome (pSS) and control individuals

    Genome-wide methylation analysis

    Genomic DNA from whole blood, CD19+ B cells and minor salivary gland biopsies was isolated using standard procedures. DNA methylation levels of 485 577 CpG sites were determined on the HM450K BeadChip (Illumina, San Diego, California, USA). Signal intensities were parsed into the Minfi R package for Subset-quantile Within Array Normalization.20–22 The post-quality control (QC) data set comprised 388 971 CpG sites. To determine differential methylation between patients with pSS and controls a linear regression model containing cell count estimates,23 ,24 age and sex as covariates was fitted. DMCs with a Bonferroni-adjusted threshold of p<1.3×10−7 were considered significant.

    Gene expression profiling of CD19+ B cells

    Expression analysis on CD19+ B cells (n=16 patients, n=23 controls) was conducted using the TruSeq stranded mRNA sample preparation kit followed by sequencing on a HiSeq2500 instrument (Illumina). QC was conducted using RNA-seQC.25 Reads were mapped with Tophat2 and analysis of differential gene expression was performed using the Cufflinks pipeline.26 ,27

    Methylation quantitative trait loci analysis

    Methylation levels were tested in PLINK for genotype association at loci that have previously shown an association with pSS with genome-wide significance.9 ,28 Quality controlled genotype data for 135 503 probes generated on the Infinium ImmunoChip (Illumina) were available for 382 of the healthy control individuals in our study. All CpG sites within a gene locus plus 100 kb flanking regions were tested against all genotypes within the same region. A Bonferroni corrected p<1.24×10−7 was considered statistically significant.

    Results

    Differential methylation in whole blood

    First, we investigated the difference in methylation levels in whole blood between patients with pSS and controls. We used reference DNA methylation signatures of flow sorted blood cells types to estimate cell counts and found reduced CD4+ and CD8+ T cells in pSS, while CD19+ B cell proportions were estimated to be similar between patients and controls (see online supplementary figure S1). We identified 11 785 (6171 hypomethylated and 5614 hypermethylated) DMCs annotated to 5623 unique genes (see figure 1 and online supplementary table S1).

    Figure 1
    Figure 1

    Manhattan plot of the genome-wide DNA methylation analysis in whole blood from patients with primary Sjögren's syndrome (pSS) and control individuals. The figure shows the −log10 p values of the association of CpG sites with pSS against the chromosomal position of the investigated sites. The hypomethylated CpG sites with a difference in mean methylation-β >0.1 are highlighted by their gene symbol annotation.

    An average difference in β-values of >0.1 between cases and controls was identified at 12 of the 11 785 DMCs, of which 11 DMCs, annotated to seven different genes, were hypomethylated in pSS (table 2). The most pronounced difference in methylation was detected for a CpG site annotated to MX dynamin-like GTPase 1 (MX1) (also referred to as MxA). MX1 is a key mediator of human antiviral immune responses and is induced by type I and type II IFNs.29 Two additional CpG sites in MX1 were found to be distinctly hypomethylated in patients. Of interest, among the top hypomethylated DMCs in pSS we further note CpG sites in the IFN-induced genes IFI44L, PARP9, PLSCR1, IFIT1, IFITM1 and HLA-A, meaning that all of the top hypomethylated sites in pSS are IFN-regulated (table 2). In addition, we detected a large number of DMCs in IFN-induced genes with a difference in methylation <0.1, for example, STAT4, NFAT5, ELF1, OAS1-3 and TREX1. Multiple DMCs with a difference in methylation <0.1 were also observed in the human leucocyte antigen (HLA) region, both major histocompatibility complex (MHC) class I and class II, the majority being hypomethylated in the patients (see online supplementary table S1). We identified one hypermethylated CpG site with an average difference in methylation-β of >0.1 annotated to EBF4, which is a transcription factor belonging to the Olf-1/EBF family with central implications in neural development and B cell maturation (table 2).30

    View this table:
    Table 2

    Differentially methylated CpG sites (DMCs) with the largest difference in mean methylation values between patients with primary Sjögren's syndrome (pSS) and control individuals in whole blood (Δβ-value >0.1)

    We then analysed the 12 top differentially methylated sites in patients with pSS stratified on the presence of anti-Sjögren's Syndrome antigen A (SSA) and/or anti-Sjögren's Syndrome antigen B (SSB)-antibodies. When analysing antibody-positive patients versus controls, a more prominent difference in mean methylation was seen, indicating that the difference in mean methylation for the IFN-induced genes is mainly driven by the antibody-positive patients (see online supplementary table S2). Analysing only patients with pSS (n=57) with early disease defined as ≤3 years from diagnosis to blood sampling, almost identical results of difference in mean methylation were seen as in the analysis of all 100 cases versus controls (data not shown).

    Association analysis of sex-chromosomal CpG sites was conducted separately for female and male individuals. We identified 85 X chromosomal CpG sites (out of 11 232 X chromosomal sites included on the HM450K array), annotated to 56 unique genes to be differentially methylated in female patients compared with female controls, with DMCs in notable genes such as VSIG4, TLR8, CD40L as well as in several microRNAs (miRNAs) (see online supplementary table S3). There were no DMCs in male individuals on the X chromosome and the Y chromosome.

    Pathway analysis of the 500 most significantly associated DMCs in whole blood identified antigen presentation, IFN signalling and graft-versus-host disease signalling as the top canonical pathways (see online supplementary table S4). The strongest gene-set enrichment in disease or function annotation of DMCs was observed for lymphohaematopoietic cancer (p=6×10−13) (see online supplementary table S5). Given the large number of DMCs between patients and controls, we analysed global DNA methylation levels in whole blood, CD19+ B cells and minor salivary gland biopsies and found no difference between patients and controls (see online supplementary figure S2).

    Functional genomic distribution and overlap with chromatin marks

    In general, DMCs were enriched in CpG island shelves and open sea regions and depleted in CpG islands and shores (figure 2A). Investigating the distribution of hypomethylated and hypermethylated DMCs separately, hypomethylated DMCs were over-represented in 5′-untranslated region (UTR), whereas hypermethylated DMCs were more than twofold enriched in 3′-UTR and moderately enriched in gene bodies (figure 2B). Analysing the intersection of pSS associated CpG sites with chromatin marks revealed that DMCs with hypomethylation in patients were enriched in enhancers (H3K4me1 and H3K27ac) and accessible chromatin (DNase I hypersensitive sites, DHS) compared with the distribution of all probes on the array. In contrast, hypermethylated DMCs were depleted for these modifications and also largely under-represented in the active promoter mark H3K4me3. On the other hand, DMCs with hypermethylation in patients were enriched for H3K36me3, which marks an actively transcribed gene body (figure 2C).

    Figure 2
    Figure 2

    Functional genomic distribution of differentially methylated CpG sites (DMCs) in whole blood from patients with primary Sjögren's syndrome (pSS) and controls. This figure illustrates the functional genomic distribution of all associated autosomal DMCs (first panel), of the hypomethylated DMCs (second panel) and of the hypermethylated DMCs (third panel) annotated (A) in relation to CpG island context, (B) in relation to gene region and (C) in relation to chromatin states in CD19+ B and CD3+ T reference cells. The colour scale indicates fold-enrichment (orange) or fold-depletion (blue) of the DMCs in each functionally annotated region. The bold numbers represent annotations to which the DMCs significantly differ compared with the distribution of probes on the HM450K array (post-QC probe set) (Bonferroni corrected X2-test p value <0.005). TSS1500, 1500 bp upstream of transcription start site (TSS); TSS200, 200 bp upstream of TSS; ncRNA, non-coding RNA; NA, probe not annotated to a defined gene property.

    Differential methylation and mRNA expression in CD19+ B cells

    Next, we analysed primary CD19+ B cells from 24 patients and 47 healthy controls and found 453 DMCs, (98 hypomethylated and 355 hypermethylated, annotated to 303 unique genes) (see online supplementary table S6). The top associated DMCs are shown in table 3. Similar to whole blood, several IFN-induced genes showed prominent hypomethylation at multiple CpG sites in CD19+ B cells from patients with pSS. In order to investigate whether differential methylation was associated with gene expression, gene expression analysis was performed in CD19+ B cells from a subset of patients and controls. Significantly upregulated expression was observed for all of the eight IFN-induced genes exhibiting DMCs with hypomethylation >0.2 (table 3). In contrast, for the two genes with hypermethylated DMCs >0.2, no significant association with gene expression was observed.

    View this table:
    Table 3

    Top DMCs and the corresponding differential gene expression in CD19+ B cells in patients with primary Sjögren's syndrome (pSS) and healthy controls

    Differential methylation in minor salivary gland biopsies

    Finally, we studied DNA methylation directly in the primary target organ of the disease, the salivary gland. Given the described IFN signature in the minor salivary glands, we hypothesised there might be epigenetic changes in IFN-induced genes.4 ,5 We found 45 DMCs annotated to 19 unique genes, where the most significant DMC showed hypomethylation in OAS2, which encodes a member of the IFN-inducible 2′-5′A synthetase family, involved in the innate immune response to viral infections (see online supplementary table S7).31

    Genetic regulation of DNA methylation at pSS risk loci

    To investigate whether the effects of established pSS risk alleles are mediated through changes in DNA methylation, we analysed genetic variants in 382 of our control individuals at the following pSS GWAS loci: DDX6-CXCR5, FAM167A-BLK, IL12A, IRF5-TNPO3, STAT4, TNIP1, and within the HLA region, for association with methylation levels in whole blood. Evidence for genetic regulation of DNA methylation, that is, significant methylation quantitative trait loci (metQTL), was identified for all pSS risk loci (in total 36 679 metQTL). Table 4 shows the metQTL for the pSS associated SNPs.9 Apart from the HLA region, the most significant metQTL was observed for the IRF5-TNPO3 SNP rs4728142 and CpG site cg04864179, both located in the IRF5 promoter. Interestingly, the IRF5-TNPO3 locus is also the non-HLA gene locus most significantly associated with pSS.9 Methylation levels in blood at cg04864179 were also directly associated with pSS in our EWAS (p=4.8×10−13, mean methylation-β difference 0.04).

    View this table:
    Table 4

    Methylation quantitative trait loci (metQTL) of genetic variants associated with primary Sjögren's syndrome (pSS) at genome-wide significance and the DNA methylation association of the most significant metQTL CpG site

    Discussion

    Here we report a comprehensive analysis of DNA methylation in pSS; in whole blood, primary CD19+ B cells and minor salivary gland biopsies. The most prominent finding was hypomethylation of IFN-regulated genes, including multiple associated CpG sites in MX1 and IFI44L, replicating previous reports in pSS.14 ,15 Notably, all top DMCs were annotated to IFN-induced genes. We also identified numerous significant, but smaller methylation differences between patients and controls most of which are novel (methylation-β value differences <0.1, online supplementary table S1). We show for the first time that hypomethylation of IFN-regulated genes in pSS CD19+ B cells is correlated with increased gene expression. PSS belongs to the systemic autoimmune diseases, together with systemic lupus erythematosus (SLE) which display an IFN signature, with several possible mechanisms underlying the IFN activation.32–34 Recently, MX1 has been suggested as a potential biomarker for disease activity and type I IFN bioactivity in pSS,35 and IFI44L is described as an indicator gene of the type I IFN signature.7 Hypomethylation in IFI44L, which was initially reported by Altorok et al14 in pSS CD4+ T cells and subsequently by Miceli-Richard et al15 in CD19+ B cells, has also been detected in multiple cell types from patients with SLE.36–39 Absher et al37 found hypomethylation of IFN-induced genes in naïve, memory and regulatory CD4+ T cells, CD19+ B cells and CD14+ monocytes from patients with SLE, with either active or quiescent disease. They conclude that epigenetic changes occur in progenitor cells independent of IFN activity. Coit et al39 reported hypomethylation of IFN-induced genes in neutrophils from patients with SLE and speculate that the exposure to IFN during the disease course may induce methylation differences that will increase the responsiveness to IFN. Taken together, epigenetic changes in different cell types in pSS and SLE makes them poised for type I IFN expression, although the exact mechanisms remain to be elucidated. In concordance with others, we found that the hypomethylation of IFN-induced genes is largely driven by the antibody-positive patients.15 Our pathway analysis of DMCs also revealed IFN signalling among the top associated pathways. Monoclonal antibodies interfering with the IFN signalling pathway currently under development in SLE might be of interest in future clinical trials in selected patients with pSS.

    Polymorphisms in genes of the type I IFN system have shown associations with pSS.9 ,10 ,40 We therefore investigated whether genetic variants in known pSS associated loci mediate disease risk by influencing methylation levels at target CpG sites. We identified significant associations of genetic variants in DDX6-CXCR5, FAM167A-BLK, IL12A, IRF5-TNPO3, STAT4, TNIP1, and within the HLA region with DNA methylation, indicating that pSS GWAS risk alleles have the potential to affect DNA methylation levels.

    For insights into functional mechanisms of DMCs in pSS we studied the regional genomic distribution of associated sites. We found distinct enrichment patterns for overlap with chromatin marks, where hypomethylated DMCs were enriched for enhancer marks and regions of open chromatin, while hypermethylated DMCs were highly under-represented in these functional regions. Hypomethylated CpG sites were also more prominent in promoter regions, indicating a putative transcriptional activation of genes with hypomethylation in these regions.

    Interestingly, distinct hypomethylation and increased gene expression of the PARP9 (BAL1) gene were found in CD19+ B cells. Whereas BAL1 has been shown to be overexpressed in diffuse large B cell lymphomas,41 none of our patients sampled for CD19+ B cells had a previous lymphoma. Disease annotation analysis of the top DMCs in whole blood also pointed to an extensive enrichment of genes associated with lymphoproliferative disorders, including B cell non-Hodgkin's lymphoma. This is intriguing since pSS is the autoimmune disease which displays the highest risk for lymphoma development.42 Whether aberrant methylation in target genes contributes to lymphoma development is an important topic for future studies.

    The strongest genetic association for pSS is found in the HLA region.9 Although none of the most significantly associated DMCs were annotated to HLA genes, we nonetheless identified a considerable number of significant signals within this region with the majority being hypomethylated in patients with pSS, possibly pointing to an increased expression of alternatively spliced antigenic transcripts in pSS.43 ,44 Analysing X chromosomal DNA methylation we observed several DMCs in genes with central roles in connecting innate and adaptive immunity, such as TLR8 and CD40L, as well as in miRNAs.45 Altered miRNA expression constitutes another epigenetic mechanism implicated in pSS pathogenesis,46 and the potential role of miRNA-223 hypomethylation in pSS pathogenesis warrants further investigation.

    In our study we used a well known reference based method for cell type estimation of whole blood samples.23 ,24 In the purified CD19+ B cells compared with the whole blood analysis, we noted a larger mean difference in methylation levels between cases and controls for many associated CpG sites, including at MX1 and IFI44L, perhaps indicating the advantage of a single cell type in the analysis. However, the smaller number of individuals from which purified cells could be obtained, compared with DNA extracted from whole blood, still meant that fewer DMCs were detected studying B cells. Due to the low absolute number of B cells in serum from patients with pSS we were not able to assess methylation patterns in different B cell subtypes. Minor salivary gland biopsies consist of different cell types including epithelial and acinar cells in both patients with pSS and controls, but also inflammatory cells in the patients’ biopsies. Therefore we cannot deduce the cell types that are responsible for the difference in methylation between diseased and normal glands, and the results must be interpreted with caution. Nevertheless, it is noteworthy that the strongest association was found in OAS2, an IFN-induced gene involved in the innate immune response.31 The minor salivary glands are known targets for IFN where both increased IFN-α levels and an IFN signature have been demonstrated.3–5

    In conclusion, our study of epigenetic profiles in multiple tissues in pSS using a large collection of patients and controls has replicated the previously reported hypomethylation of IFN-regulated genes in pSS and identified numerous new associations. We report hypomethylation in regulatory enhancer and promotor regions and show for the first time that hypomethylation of IFN-regulated genes in B cells corresponds to an increase in gene expression. Evidence for genetic control of methylation levels at known pSS risk loci is presented. Independent replication in cells from patients with pSS and controls will be required to confirm these novel findings. Studying the epigenetic basis of pSS will hopefully increase our understanding of the disease mechanisms and guide the search for novel and more specific therapeutic targets.

    Acknowledgments

    The authors thank Rezvan Kiani and Karolina Tandre for collecting samples from patients and controls. Genotyping, epigenotyping and RNA-sequencing were performed at the SNP&SEQ Technology Platform at the National Genomics Infrastructure (NGI) hosted by Science for Life Laboratory in Uppsala, Sweden (http://www.genotyping.se; http://www.sequencing.se). The authors thank Andrei Alexsson, Tomas Axelsson, Anna Haukkala, Maria Hägglund, Charlotta Jakobsson and Anders Lundmark for excellent technical assistance. The authors especially thank all patients and blood donors who contributed samples to this study.

    References

    1. ↵
      1. Jonsson R,
      2. Bolstad AI,
      3. Brokstad KA, et al
      . Sjogren's syndrome—a plethora of clinical and immunological phenotypes with a complex genetic background. Ann N Y Acad Sci 2007;1108:433–47. doi:10.1196/annals.1422.046
      OpenUrlCrossRefPubMedWeb of Science
    2. ↵
      1. Nocturne G,
      2. Mariette X
      . Advances in understanding the pathogenesis of primary Sjögren's syndrome. Nat Rev Rheumatol 2013;9:544–56. doi:10.1038/nrrheum.2013.110
      OpenUrlCrossRefPubMed
    3. ↵
      1. Båve U,
      2. Nordmark G,
      3. Lövgren T, et al
      . Activation of the type I interferon system in primary Sjögren's syndrome: a possible etiopathogenic mechanism. Arthritis Rheum 2005;52:1185–95. doi:10.1002/art.20998
      OpenUrlCrossRefPubMedWeb of Science
    4. ↵
      1. Gottenberg JE,
      2. Cagnard N,
      3. Lucchesi C, et al
      . Activation of IFN pathways and plasmacytoid dendritic cell recruitment in target organs of primary Sjögren's syndrome. Proc Natl Acad Sci USA 2006;103:2770–5. doi:10.1073/pnas.0510837103
      OpenUrlAbstract/FREE Full Text
    5. ↵
      1. Hjelmervik TO,
      2. Petersen K,
      3. Jonassen I, et al
      . Gene expression profiling of minor salivary glands clearly distinguishes primary Sjögren's syndrome patients from healthy control subjects. Arthritis Rheum 2005;52:1534–44. doi:10.1002/art.21006
      OpenUrlCrossRefPubMedWeb of Science
    6. ↵
      1. Emamian ES,
      2. Leon JM,
      3. Lessard CJ, et al
      . Peripheral blood gene expression profiling in Sjögren's syndrome. Genes Immun 2009;10:285–96. doi:10.1038/gene.2009.20
      OpenUrlCrossRefPubMedWeb of Science
    7. ↵
      1. Brkic Z,
      2. Maria NI,
      3. van Helden-Meeuwsen CG, et al
      . Prevalence of interferon type I signature in CD14 monocytes of patients with Sjögren's syndrome and association with disease activity and BAFF gene expression. Ann Rheum Dis 2013;72:728–35. doi:10.1136/annrheumdis-2012-201381
      OpenUrlAbstract/FREE Full Text
    8. ↵
      1. Theander E,
      2. Henriksson G,
      3. Ljungberg O, et al
      . Lymphoma and other malignancies in primary Sjögren's syndrome: a cohort study on cancer incidence and lymphoma predictors. Ann Rheum Dis 2006;65:796–803. doi:10.1136/ard.2005.041186
      OpenUrlAbstract/FREE Full Text
    9. ↵
      1. Lessard CJ,
      2. Li H,
      3. Adrianto I, et al
      . Variants at multiple loci implicated in both innate and adaptive immune responses are associated with Sjögren's syndrome. Nat Genet 2013;45:1284–92. doi:10.1038/ng.2792
      OpenUrlCrossRefPubMed
    10. ↵
      1. Li Y,
      2. Zhang K,
      3. Chen H, et al
      . A genome-wide association study in Han Chinese identifies a susceptibility locus for primary Sjögren's syndrome at 7q11.23. Nat Genet 2013;45:1361–5. doi:10.1038/ng.2779
      OpenUrlCrossRefPubMed
    11. ↵
      1. Ice JA,
      2. Li H,
      3. Adrianto I, et al
      . Genetics of Sjögren's syndrome in the genome-wide association era. J Autoimmun 2012;39:57–63. doi:10.1016/j.jaut.2012.01.008
      OpenUrlCrossRefPubMedWeb of Science
    12. ↵
      1. Konsta OD,
      2. Thabet Y,
      3. Le Dantec C, et al
      . The contribution of epigenetics in Sjögren's Syndrome. Front Genet 2014;5:71. doi:10.3389/fgene.2014.00071
      OpenUrlPubMed
    13. ↵
      1. Bibikova M,
      2. Barnes B,
      3. Tsan C, et al
      . High density DNA methylation array with single CpG site resolution. Genomics 2011;98:288–95. doi:10.1016/j.ygeno.2011.07.007
      OpenUrlCrossRefPubMedWeb of Science
    14. ↵
      1. Altorok N,
      2. Coit P,
      3. Hughes T, et al
      . Genome-wide DNA methylation patterns in naive CD4+ T cells from patients with primary Sjögren's syndrome. Genomics 2014;66:731–9. doi:10.1002/art.38264
      OpenUrl
    15. ↵
      1. Miceli-Richard C,
      2. Wang-Renault SF,
      3. Boudaoud S, et al
      . Overlap between differentially methylated DNA regions in blood B lymphocytes and genetic at-risk loci in primary Sjogren's syndrome. Ann Rheum Dis 2016;75:93340.
      OpenUrlAbstract/FREE Full Text
    16. ↵
      1. Bernstein BE,
      2. Stamatoyannopoulos JA,
      3. Costello JF, et al
      . The NIH Roadmap Epigenomics Mapping Consortium. Nat Biotechnol 2010;28:1045–8. doi:10.1038/nbt1010-1045
      OpenUrlCrossRefPubMedWeb of Science
    17. ↵
      1. Vitali C,
      2. Bombardieri S,
      3. Jonsson R, et al
      . Classification criteria for Sjogren's syndrome: a revised version of the European criteria proposed by the American-European Consensus Group. Ann Rheum Dis 2002;61:554–8. doi:10.1136/ard.61.6.554
      OpenUrlAbstract/FREE Full Text
    18. ↵
      1. Nordmark G,
      2. Wang C,
      3. Vasaitis L, et al
      . Association of genes in the NF-kB pathway with antibody-positive primary Sjogren's syndrome. Scand J Immunol 2013;78:447–54. doi:10.1111/sji.12101
      OpenUrlCrossRefPubMed
    19. ↵
      1. Berggren O,
      2. Alexsson A,
      3. Morris DL, et al
      . IFN-alpha production by plasmacytoid dendritic cell associations with polymorphisms in gene loci related to autoimmune and inflammatory diseases. Hum Mol Genet 2015;24:3571–81. doi:10.1093/hmg/ddv095
      OpenUrlAbstract/FREE Full Text
    20. ↵
      R Core Development Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing, 2011.
    21. ↵
      1. Aryee MJ,
      2. Jaffe AE,
      3. Corrada-Bravo H, et al
      . Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 2014;30:1363–9. doi:10.1093/bioinformatics/btu049
      OpenUrlAbstract/FREE Full Text
    22. ↵
      1. Maksimovic J,
      2. Gordon L,
      3. Oshlack A
      . SWAN: Subset-quantile within array normalization for illumina infinium HumanMethylation450 BeadChips. Genome Biol 2012;13:R44. doi:10.1186/gb-2012-13-6-r44
      OpenUrlCrossRefPubMed
    23. ↵
      1. Houseman EA,
      2. Accomando WP,
      3. Koestler DC, et al
      . DNA methylation arrays as surrogate measures of cell mixture distribution. BMC bioinformatics 2012;13:86. doi:10.1186/1471-2105-13-86
      OpenUrlCrossRefPubMed
    24. ↵
      1. Reinius LE,
      2. Acevedo N,
      3. Joerink M, et al
      . Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PLoS ONE 2012;7:e41361. doi:10.1371/journal.pone.0041361
      OpenUrlCrossRefPubMed
    25. ↵
      1. DeLuca DS,
      2. Levin JZ,
      3. Sivachenko A, et al
      . RNA-SeQC: RNA-seq metrics for quality control and process optimization. Bioinformatics 2012;28:1530–2. doi:10.1093/bioinformatics/bts196
      OpenUrlAbstract/FREE Full Text
    26. ↵
      1. Trapnell C,
      2. Roberts A,
      3. Goff L, et al
      . Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc 2012;7:562–78. doi:10.1038/nprot.2012.016
      OpenUrlCrossRefPubMed
    27. ↵
      1. Kim D,
      2. Pertea G,
      3. Trapnell C, et al
      . TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 2013;14:R36. doi:10.1186/gb-2013-14-4-r36
      OpenUrlCrossRefPubMed
    28. ↵
      1. Purcell S,
      2. Neale B,
      3. Todd-Brown K, et al
      . PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559–75. doi:10.1086/519795
      OpenUrlCrossRefPubMed
    29. ↵
      1. Haller O,
      2. Kochs G
      . Human MxA protein: an interferon-induced dynamin-like GTPase with broad antiviral activity. J Interferon Cytokine Res 2011;31:79–87. doi:10.1089/jir.2010.0076
      OpenUrlCrossRefPubMedWeb of Science
    30. ↵
      1. Wang SS,
      2. Betz AG,
      3. Reed RR
      . Cloning of a novel Olf-1/EBF-like gene, O/E-4, by degenerate oligo-based direct selection. Mol Cell Neurosci 2002;20:404–14. doi:10.1006/mcne.2002.1138
      OpenUrlCrossRefPubMedWeb of Science
    31. ↵
      1. Hovanessian AG,
      2. Justesen J
      . The human 2′-5′oligoadenylate synthetase family: unique interferon-inducible enzymes catalyzing 2′-5′ instead of 3′-5′ phosphodiester bond formation. Biochimie 2007;89:779–88. doi:10.1016/j.biochi.2007.02.003
      OpenUrlCrossRefPubMedWeb of Science
    32. ↵
      1. Nordmark G,
      2. Eloranta ML,
      3. Ronnblom L
      . Primary Sjogren's syndrome and the type I interferon system. Curr Pharm Biotechnol 2012;13:2054–62. doi:10.2174/138920112802273290
      OpenUrlCrossRefPubMed
    33. ↵
      1. Nezos A,
      2. Gravani F,
      3. Tassidou A, et al
      . Type I and II interferon signatures in Sjogren's syndrome pathogenesis: contributions in distinct clinical phenotypes and Sjogren's related lymphomagenesis. J Autoimmun 2015;63:47–58. doi:10.1016/j.jaut.2015.07.002
      OpenUrlCrossRefPubMed
    34. ↵
      1. Rönnblom L,
      2. Alm GV,
      3. Eloranta ML
      . The type I interferon system in the development of lupus. Semin Immunol 2011;23:113–21. doi:10.1016/j.smim.2011.01.009
      OpenUrlCrossRefPubMed
    35. ↵
      1. Maria NI,
      2. Brkic Z,
      3. Waris M, et al
      . MxA as a clinically applicable biomarker for identifying systemic interferon type I in primary Sjogren's syndrome. Ann Rheum Dis 2014;73:1052–9. doi:10.1136/annrheumdis-2012-202552
      OpenUrlAbstract/FREE Full Text
    36. ↵
      1. Coit P,
      2. Jeffries M,
      3. Altorok N, et al
      . 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. doi:10.1016/j.jaut.2013.04.003
      OpenUrlCrossRefPubMedWeb of Science
    37. ↵
      1. Absher DM,
      2. Li X,
      3. Waite LL, et al
      . 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:e1003678. doi:10.1371/journal.pgen.1003678
      OpenUrlCrossRefPubMed
    38. ↵
      1. Zhao M,
      2. Liu S,
      3. Luo S, et al
      . DNA methylation and mRNA and microRNA expression of SLE CD4+ T cells correlate with disease phenotype. J Autoimmun 2014;54:127–36. doi:10.1016/j.jaut.2014.07.002
      OpenUrlCrossRefPubMed
    39. ↵
      1. Coit P,
      2. Yalavarthi S,
      3. Ognenovski M, et al
      . Epigenome profiling reveals significant DNA demethylation of interferon signature genes in lupus neutrophils. J Autoimmun 2015;58:59–66. doi:10.1016/j.jaut.2015.01.004
      OpenUrlCrossRefPubMed
    40. ↵
      1. Nordmark G,
      2. Kristjansdottir G,
      3. Theander E, et al
      . Association of EBF1, FAM167A(C8orf13)-BLK and TNFSF4 gene variants with primary Sjögren's syndrome. Genes Immun 2011;12:100–9. doi:10.1038/gene.2010.44
      OpenUrlCrossRefPubMedWeb of Science
    41. ↵
      1. Camicia R,
      2. Bachmann SB,
      3. Winkler HC, et al
      . BAL1/ARTD9 represses the anti-proliferative and pro-apoptotic IFNg-STAT1-IRF1-p53 axis in diffuse large B-cell lymphoma. J Cell Sci 2013;126(Pt 9):1969–80. doi:10.1242/jcs.118174
      OpenUrlAbstract/FREE Full Text
    42. ↵
      1. Baecklund E,
      2. Smedby KE,
      3. Sutton LA, et al
      . Lymphoma development in patients with autoimmune and inflammatory disorders—what are the driving forces? Semin Cancer Biol 2014;24:61–70. doi:10.1016/j.semcancer.2013.12.001
      OpenUrlCrossRefPubMedWeb of Science
    43. ↵
      1. Lev Maor G,
      2. Yearim A,
      3. Ast G
      . The alternative role of DNA methylation in splicing regulation. Trends Genet 2015;31:274–80. doi:10.1016/j.tig.2015.03.002
      OpenUrlCrossRefPubMed
    44. ↵
      1. Graves M,
      2. Benton M,
      3. Lea R, et al
      . Methylation differences at the HLA-DRB1 locus in CD4+ T-Cells are associated with multiple sclerosis. Mult Scler 2013;20:1033–41. doi:10.1177/1352458513516529
      OpenUrlPubMed
    45. ↵
      1. Haneklaus M,
      2. Gerlic M,
      3. O'Neill LA, et al
      . miR-223: infection, inflammation and cancer. J Intern Med 2013;274:215–26. doi:10.1111/joim.12099
      OpenUrlCrossRefPubMed
    46. ↵
      1. Alevizos I,
      2. Alexander S,
      3. Turner RJ, et al
      . MicroRNA expression profiles as biomarkers of minor salivary gland inflammation and dysfunction in Sjögren's syndrome. Arthritis Rheum 2011;63:535–44. doi:10.1002/art.30131
      OpenUrlCrossRefPubMedWeb of Science

    Supplementary materials

    • Supplementary Data

      This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

    Footnotes

    • Handling editor Tore K Kvien

    • JI-K and JKS contributed equally.

    • Contributors JI-K, JKS, KBN, RO, A-CS and GN designed the study; LS, LR, M-LE and GN collected patient and control material and clinical data; JI-K and JKS performed the experiments; JI-K, JKS, JCA and JN analysed the data; JI-K, JKS and GN drafted the manuscript and all authors read and accepted the final version of the manuscript. JI-K and JKS contributed equally to the study.

    • Funding This study was supported by grants from the Knut and Alice Wallenberg Foundation, the Swedish Research Council for Medicine and Health (Dnr 521-2014-2263 ACS and Dnr 521-2013-2830 LR), the Gustav V: 80-year Foundation, Combine, the Swedish Society of Medicine and the Swedish Rheumatism Association. JKS was supported by a Swedish Research Council postdoc grant (Dnr 350-2012-256). The SNP&SEQ Technology Platform is supported by the Swedish Research Council (VR-RFI), Science for Life Laboratory and the Knut and Alice Wallenberg Foundation.

    • Competing interests None declared.

    • Patient consent Obtained.

    • Ethics approval The study was approved by the Regional Ethics board in Uppsala No. 97358, 217/2006 and 013/2009.

    • Provenance and peer review Not commissioned; externally peer reviewed.

    • Data sharing statement Normalised or raw intensity data of the HM450K BeadChips are available upon request from the authors on a collaborative basis.