We analyzed a comprehensive set of single-nucleotide polymorphisms (SNPs) and length polymorphisms in the interferon regulatory factor 5 (IRF5) gene for their association with the autoimmune disease systemic lupus erythematosus (SLE) in 485 Swedish patients and 563 controls. We found 16 SNPs and two length polymorphisms that display association with SLE (P < 0.0005, OR > 1.4). Using a Bayesian model selection and averaging approach we identified parsimonious models with exactly two variants of IRF5 that are independently associated with SLE. The variants of IRF5 with the highest posterior probabilities (1.00 and 0.71, respectively) of being causal in SLE are a SNP (rs10488631) located 3' of IRF5, and a novel CGGGG insertion-deletion (indel) polymorphism located 64 bp upstream of the first untranslated exon (exon 1A) of IRF5. The CGGGG indel explains the association signal from multiple SNPs in the IRF5 gene, including rs2004640, rs10954213 and rs729302 previously considered to be causal variants in SLE. The CGGGG indel contains three or four repeats of the sequence CGGGG with the longer allele containing an additional SP1 binding site as the risk allele for SLE. Using electrophoretic mobility shift assays we show increased binding of protein to the risk allele of the CGGGG indel and using a minigene reporter assay we show increased expression of IRF5 mRNA from a promoter containing this allele. Increased expression of IRF5 protein was observed in peripheral blood mononuclear cells from SLE patients carrying the risk allele of the CGGGG indel. We have found that the same IRF5 allele also confers risk for inflammatory bowel diseases and multiple sclerosis, suggesting a general role for IRF5 in autoimmune diseases.

Systemic lupus erythematosus (SLE) is the prototype autoimmune disease where genes regulated by type I interferon (IFN) are over-expressed and contribute to the disease pathogenesis. Because signal transducer and activator of transcription 4 (STAT4) plays a key role in the type I IFN receptor signaling, we performed a candidate gene study of a comprehensive set of single nucleotide polymorphism (SNPs) in STAT4 in Swedish patients with SLE. We found that 10 out of 53 analyzed SNPs in STAT4 were associated with SLE, with the strongest signal of association (P = 7.1 × 10−8) for two perfectly linked SNPs rs10181656 and rs7582694. The risk alleles of these 10 SNPs form a common risk haplotype for SLE (P = 1.7 × 10−5). According to conditional logistic regression analysis the SNP rs10181656 or rs7582694 accounts for all of the observed association signal. By quantitative analysis of the allelic expression of STAT4 we found that the risk allele of STAT4 was over-expressed in primary human cells of mesenchymal origin, but not in B-cells, and that the risk allele of STAT4 was over-expressed (P = 8.4 × 10−5) in cells carrying the risk haplotype for SLE compared with cells with a non-risk haplotype. The risk allele of the SNP rs7582694 in STAT4 correlated to production of anti-dsDNA (double-stranded DNA) antibodies and displayed a multiplicatively increased, 1.82-fold risk of SLE with two independent risk alleles of the IRF5 (interferon regulatory factor 5) gene.

Systemic lupus erythematosus (SLE) patients have increased levels of interferon-alfa (IFN-alpha) in the circulation but a reduced number of functionally intact natural IFN-alpha producing cells (IPC) in peripheral blood. In search for tissue localisation of activated IPC, we investigated skin biopsies from SLE patients for the occurrence of such cells. Eleven SLE patients with inflammatory skin lesions and six healthy controls were biopsied. An immunohistochemical technique (IH) and in situ hybridisation (ISH) were used to detect intracellular IFN-alpha protein and IFN-alpha mRNA, respectively. In all 11 biopsies from SLE lesions, a high number of IPC were detected by IH. In the nonlesional SLE biopsies we could also demonstrate IPC in 10/11 patients. In 6/11 SLE patients, IFN-alpha mRNA containing cells could be detected in the specimens. A low number of IPC were detected in 1/6 healthy controls by IH, but no ISH positive cells were seen. Our results demonstrate that SLE patients have active IPC in both dermal lesions and in noninflammatory skin. A recruitment of IPC from blood to peripheral tissues may explain the low number of circulating natural IPC in SLE patients. Because the type I IFN system is involved in the SLE disease process, these results are of interest for the understanding of the pathogenesis in SLE.

ObjectivesSystemic lupus erythematosus (SLE) is a chronic autoimmune condition with heterogeneous presentation and complex aetiology where DNA methylation changes are emerging as a contributing factor. In order to discover novel epigenetic associations and investigate their relationship to genetic risk for SLE, we analysed DNA methylation profiles in a large collection of patients with SLE and healthy individuals.MethodsDNA extracted from blood from 548 patients with SLE and 587 healthy controls were analysed on the Illumina HumanMethylation 450 k BeadChip, which targets 485 000 CpG sites across the genome. Single nucleotide polymorphism (SNP) genotype data for 196 524 SNPs on the Illumina ImmunoChip from the same individuals were utilised for methylation quantitative trait loci (cis-meQTLs) analyses.ResultsWe identified and replicated differentially methylated CpGs (DMCs) in SLE at 7245 CpG sites in the genome. The largest methylation differences were observed at type I interferon-regulated genes which exhibited decreased methylation in SLE. We mapped cis-meQTLs and identified genetic regulation of methylation levels at 466 of the DMCs in SLE. The meQTLs for DMCs in SLE were enriched for genetic association to SLE, and included seven SLE genome-wide association study (GWAS) loci: PTPRC (CD45), MHC-class III, UHRF1BP1, IRF5, IRF7, IKZF3 and UBE2L3. In addition, we observed association between genotype and variance of methylation at 20 DMCs in SLE, including at the HLA-DQB2 locus.ConclusionsOur results suggest that several of the genetic risk variants for SLE may exert their influence on the phenotype through alteration of DNA methylation levels at regulatory regions of target genes.

ObjectivesTo investigate associations between a high genetic disease risk and disease severity in patients with systemic lupus erythematosus (SLE).MethodsPatients with SLE (n=1001, discovery cohort and n=5524, replication cohort) and healthy controls (n=2802 and n=9859) were genotyped using a 200K Immunochip single nucleotide polymorphism array. A genetic risk score (GRS) was assigned to each individual based on 57 SLE risk loci.ResultsSLE was more prevalent in the high, compared with the low, GRS-quartile (OR 12.32 (9.53 to 15.71), p=7.9×10–86 and OR 7.48 (6.73 to 8.32), p=2.2×10–304 for the discovery and the replication cohorts, respectively). In the discovery cohort, patients in the high GRS-quartile had a 6-year earlier mean disease onset (HR 1.47 (1.22 to 1.75), p=4.3×10–5), displayed higher prevalence of damage accrual (OR 1.47 (1.06 to 2.04), p=2.0×10–2), renal disorder (OR 2.22 (1.50 to 3.27), p=5.9×10–5), anti-dsDNA (OR 1.83 (1.19 to 2.81), p=6.1×10–3), end-stage renal disease (ESRD) (OR 5.58 (1.50 to 20.79), p=1.0×10–2), proliferative nephritis (OR 2.42 (1.30 to 4.49), p=5.1×10–3), anti-cardiolipin-IgG (OR 1.89 (1.13 to 3.18), p=1.6×10–2), anti-β2-glycoprotein-I-IgG (OR 2.29 (1.29 to 4.06), p=4.8×10–3) and positive lupus anticoagulant test (OR 2.12 (1.16 to 3.89), p=1.5×10–2) compared with patients in the low GRS-quartile. Survival analysis showed earlier onset of the first organ damage (HR 1.51 (1.04 to 2.25), p=3.7×10–2), first cardiovascular event (HR 1.65 (1.03 to 2.64), p=2.6×10–2), nephritis (HR 2.53 (1.72 to 3.71), p=9.6×10–7), ESRD (HR 6.78 (1.78 to 26.86), p=6.5×10–3) and decreased overall survival (HR 1.83 (1.02 to 3.30), p=4.3×10–2) in high to low quartile comparison.ConclusionsA high GRS is associated with increased risk of organ damage, renal dysfunction and all-cause mortality. Our results indicate that genetic profiling may be useful for predicting outcomes in patients with SLE.

Sjögren’s syndrome (SS) is a common, autoimmune exocrinopathy distinguished by keratoconjunctivitis sicca and xerostomia. Patients frequently develop serious complications including lymphoma, pulmonary dysfunction, neuropathy, vasculitis, and debilitating fatigue. Dysregulation of type I interferon (IFN) pathway is a prominent feature of SS and is correlated with increased autoantibody titers and disease severity. To identify genetic determinants of IFN pathway dysregulation in SS, we performed cis-expression quantitative trait locus (eQTL) analyses focusing on differentially expressed type I IFN-inducible transcripts identified through a transcriptome profiling study. Multiple cis-eQTLs were associated with transcript levels of 2'-5'-oligoadenylate synthetase 1 (OAS1) peaking at rs10774671 (PeQTL = 6.05 × 10−14). Association of rs10774671 with SS susceptibility was identified and confirmed through meta-analysis of two independent cohorts (Pmeta = 2.59 × 10−9; odds ratio = 0.75; 95% confidence interval = 0.66–0.86). The risk allele of rs10774671 shifts splicing of OAS1 from production of the p46 isoform to multiple alternative transcripts, including p42, p48, and p44. We found that the isoforms were differentially expressed within each genotype in controls and patients with and without autoantibodies. Furthermore, our results showed that the three alternatively spliced isoforms lacked translational response to type I IFN stimulation. The p48 and p44 isoforms also had impaired protein expression governed by the 3' end of the transcripts. The SS risk allele of rs10774671 has been shown by others to be associated with reduced OAS1 enzymatic activity and ability to clear viral infections, as well as reduced responsiveness to IFN treatment. Our results establish OAS1 as a risk locus for SS and support a potential role for defective viral clearance due to altered IFN response as a genetic pathophysiological basis of this complex autoimmune disease.

ObjectivesGenetic variations in TNFAIP3 (A20) de-ubiquitinase (DUB) domain increase the risk of systemic lupus erythematosus (SLE) and rheumatoid arthritis. A20 is a negative regulator of NF-κB but the role of its DUB domain and related genetic variants remain unclear. We aimed to study the functional effects of A20 DUB-domain alterations in immune cells and understand its link to SLE pathogenesis.MethodsCRISPR/Cas9 was used to generate human U937 monocytes with A20 DUB-inactivating C103A knock-in (KI) mutation. Whole genome RNA-sequencing was used to identify differentially expressed genes between WT and C103A KI cells. Functional studies were performed in A20 C103A U937 cells and in immune cells from A20 C103A mice and genotyped healthy individuals with A20 DUB polymorphism rs2230926. Neutrophil extracellular trap (NET) formation was addressed ex vivo in neutrophils from A20 C103A mice and SLE-patients with rs2230926.ResultsGenetic disruption of A20 DUB domain in human and murine myeloid cells did not give rise to enhanced NF-κB signalling. Instead, cells with C103A mutation or rs2230926 polymorphism presented an upregulated expression of PADI4, an enzyme regulating protein citrullination and NET formation, two key mechanisms in autoimmune pathology. A20 C103A cells exhibited enhanced protein citrullination and extracellular trap formation, which could be suppressed by selective PAD4 inhibition. Moreover, SLE-patients with rs2230926 showed increased NETs and increased frequency of autoantibodies to citrullinated epitopes.ConclusionsWe propose that genetic alterations disrupting the A20 DUB domain mediate increased susceptibility to SLE through the upregulation of PADI4 with resultant protein citrullination and extracellular trap formation.

Systemic lupus erythematosus (SLE, OMIM 152700) is a systemic autoimmune disease with a complex etiology. The mode of inheritance of the genetic risk beyond familial SLE cases is currently unknown. Additionally, the contribution of heterozygous variants in genes known to cause monogenic SLE is not fully understood. Whole-genome sequencing of DNA samples from 71 Swedish patients with SLE and their healthy biological parents was performed to investigate the general genetic risk of SLE using known SLE GWAS risk loci identified using the ImmunoChip, variants in genes associated to monogenic SLE, and the mode of inheritance of SLE risk alleles in these families. A random forest model for predicting genetic risk for SLE showed that the SLE risk variants were mainly inherited from one of the parents. In the 71 patients, we detected a significant enrichment of ultra-rare ( ≤ 0.1%) missense and nonsense mutations in 22 genes known to cause monogenic forms of SLE. We identified one previously reported homozygous nonsense mutation in the C1QC (Complement C1q C Chain) gene, which explains the immunodeficiency and severe SLE phenotype of that patient. We also identified seven ultra-rare, coding heterozygous variants in five genes ( C1S, DNASE1L3, DNASE1, IFIH1 , and RNASEH2A ) involved in monogenic SLE. Our findings indicate a complex contribution to the overall genetic risk of SLE by rare variants in genes associated with monogenic forms of SLE. The rare variants were inherited from the other parent than the one who passed on the more common risk variants leading to an increased genetic burden for SLE in the child. Higher frequency SLE risk variants are mostly passed from one of the parents to the offspring affected with SLE. In contrast, the other parent, in seven cases, contributed heterozygous rare variants in genes associated with monogenic forms of SLE, suggesting a larger impact of rare variants in SLE than hitherto reported. Electronic supplementary material The online version of this article (10.1007/s00439-018-01966-7) contains supplementary material, which is available to authorized users.