A survey of ficolin-3 activity in Systemic Lupus Erythematosus reveals a link to hematological disease manifestations and autoantibody profile.

The complement system plays a central role in the pathogenesis of Systemic Lupus Erythematosus (SLE), but most studies have focused on the classical pathway. Ficolin-3 is the main initiator of the lectin pathway of complement in humans, but its role in systemic autoimmune disease has not been conclusively determined. Here, we combined biochemical and genetic approaches to assess the contribution of ficolin-3 to SLE risk and disease manifestations. Ficolin-3 activity was measured by a functional assay in serum or plasma samples from Swedish SLE patients (n = 786) and controls matched for age and sex (n = 566). Genetic variants in an extended 300 kb genomic region spanning the FCN3 locus were analyzed for their association with ficolin-3 activity and SLE manifestations in a Swedish multicenter cohort (n = 985). Patients with ficolin-3 activity in the highest tertile showed a strong enrichment in an SLE cluster defined by anti-Sm/DNA/nucleosome antibodies (OR 3.0, p < 0.001) and had increased rates of hematological disease (OR 1.4, p = 0.078) and lymphopenia (OR = 1.6, p = 0.039). Genetic variants associated with low ficolin-3 activity mapped to an extended haplotype in high linkage disequilibrium upstream of the FCN3 gene. Patients carrying the lead genetic variant associated with low ficolin-3 activity had a lower frequency of hematological disease (OR 0.67, p = 0.018) and lymphopenia (OR 0.63, p = 0.031) and fewer autoantibodies (p = 0.0019). Loss-of-function variants in the FCN3 gene were not associated with SLE, but four (0.5 %) SLE patients developed acquired ficolin-3 deficiency where ficolin-3 activity in serum was depleted following diagnosis of SLE. Taken together, our results provide genetic and biochemical evidence that implicate the lectin pathway in hematological SLE manifestations. We also identify lectin pathway activation through ficolin-3 as a factor that contributes to the autoantibody response in SLE.

B cell polygenic risk scores associate with anti-dsDNA antibodies and nephritis in systemic lupus erythematosus.

OBJECTIVE: B cell function and autoantibodies are important in SLE pathogenesis. In this work, we aimed to investigate the impact of cumulative SLE B cell genetics on SLE subphenotype and autoantibody profile. METHODS: Female patients with SLE (n=1248) and healthy controls (n=400) were genotyped using Illumina's Global Screening Array. Two polygenic risk scores (PRSs), one representing B cell genes and the other B cell activation genes, were calculated for each individual using risk loci for SLE in genes assigned to B cell-related pathways according to the Kyoto Encyclopedia of Genes and Genomes, Gene Ontology and Reactome Databases. RESULTS: Double-stranded DNA (dsDNA) antibodies were more prevalent among patients with a high compared with a low SLE B cell PRS (OR 1.47 (1.07 to 2.01), p=0.018), and effect sizes were augmented in patients with human leucocyte antigen (HLA) risk haplotypes HLA-DRB1*03:01 and HLA-DRB1*15:01 (DRB1*03/15 -/- (OR 0.99 (0.56 to 1.77), p=0.98; DRB1*03/15 +/- or -/+ (OR 1.64 (1.06 to 2.54), p=0.028; and DRB1*03/15 +/+ (OR 4.47 (1.21 to 16.47), p=0.024). Further, a high compared with a low B cell PRS was associated with low complement levels in DRB1*03/15 +/+ patients (OR 3.92 (1.22 to 12.64), p=0.022). The prevalence of lupus nephritis (LN) was higher in patients with a B cell activation PRS above the third quartile compared with patients below (OR 1.32 (1.00 to 1.74), p=0.048). CONCLUSIONS: High genetic burden related to B cell function is associated with dsDNA antibody development and LN. Assessing B cell PRSs may be important in order to determine immunological pathways influencing SLE and to predict clinical phenotype.

Distinct HLA associations with autoantibody-defined subgroups in idiopathic inflammatory myopathies.

BACKGROUND: In patients with idiopathic inflammatory myopathies (IIM), autoantibodies are associated with specific clinical phenotypes suggesting a pathogenic role of adaptive immunity. We explored if autoantibody profiles are associated with specific HLA genetic variants and clinical manifestations in IIM. METHODS: We included 1348 IIM patients and determined the occurrence of 14 myositis-specific or -associated autoantibodies. We used unsupervised cluster analysis to identify autoantibody-defined subgroups and logistic regression to estimate associations with clinical manifestations, HLA-DRB1, HLA-DQA1, HLA-DQB1 alleles, and amino acids imputed from genetic information of HLA class II and I molecules. FINDINGS: We identified eight subgroups with the following dominant autoantibodies: anti-Ro52, -U1RNP, -PM/Scl, -Mi2, -Jo1, -Jo1/Ro52, -TIF1γ or negative for all analysed autoantibodies. Associations with HLA-DRB1∗11, HLA-DRB1∗15, HLA-DQA1∗03, and HLA-DQB1∗03 were present in the anti-U1RNP-dominated subgroup. HLA-DRB1∗03, HLA-DQA1∗05, and HLA-DQB1∗02 alleles were overrepresented in the anti-PM/Scl and anti-Jo1/Ro52-dominated subgroups. HLA-DRB1∗16, HLA-DRB1∗07 alleles were most frequent in anti-Mi2 and HLA-DRB1∗01 and HLA-DRB1∗07 alleles in the anti-TIF1γ subgroup. The HLA-DRB1∗13, HLA-DQA1∗01 and HLA-DQB1∗06 alleles were overrepresented in the negative subgroup. Significant signals from variations in class I molecules were detected in the subgroups dominated by anti-Mi2, anti-Jo1/Ro52, anti-TIF1γ, and the negative subgroup. INTERPRETATION: Distinct HLA class II and I associations were observed for almost all autoantibody-defined subgroups. The associations support autoantibody profiles use for classifying IIM which would likely reflect underlying pathogenic mechanisms better than classifications based on clinical symptoms and/or histopathological features. FUNDING: See a detailed list of funding bodies in the Acknowledgements section at the end of the manuscript.

Mer-tyrosine kinase: a novel susceptibility gene for SLE related end-stage renal disease.

OBJECTIVE: Lupus nephritis (LN) is a common and severe manifestation of SLE. The genetic risk for nephritis and progression to end-stage renal disease (ESRD) in patients with LN remains unclear. Herein, we aimed to identify novel genetic associations with LN, focusing on subphenotypes and ESRD. METHODS: We analysed genomic data on 958 patients with SLE (discovery cohort: LN=338) with targeted sequencing data from 1832 immunological pathway genes. We used an independent multiethnic cohort comprising 1226 patients with SLE (LN=603) as a replication dataset. Detailed functional annotation and functional epigenomic enrichment analyses were applied to predict functional effects of the candidate variants. RESULTS: A genetic variant (rs56097910) within the MERTK gene was associated with ESRD in both cohorts, meta-analysis OR=5.4 (2.8 to 10.6); p=1.0×10-6. We observed decreased methylation levels in peripheral blood cells from SLE patients with ESRD, compared with patients without renal SLE (p=2.7×10-4), at one CpG site (cg16333401) in close vicinity to the transcription start site of MERTK and located in a DNAse hypersensitivity region in T and B cells. Rs56097910 is linked to altered MERTK expression in kidney tissue in public eQTL databases. Two loci were replicated for association with proliferative LN: PRDM1 (rs6924535, pmeta=1.6×10-5, OR=0.58) and APOA1BP (NAXE) (rs942960, pmeta=1.2×10-5, OR=2.64). CONCLUSION: We identified a novel genetic risk locus, MERTK, associated with SLE-ESRD using the data from two large SLE cohorts. Through DNA methylation analysis and functional annotation, we showed that the risk could be mediated through regulation of gene expression. Our results suggest that variants in the MERTK gene are important for the risk of developing SLE-ESRD and suggest a role for PRDM1 and APOA1BP in proliferative LN.

Strong Association of Combined Genetic Deficiencies in the Classical Complement Pathway With Risk of Systemic Lupus Erythematosus and Primary Sjögren's Syndrome.

OBJECTIVE: Complete genetic deficiency of the complement component C2 is a strong risk factor for monogenic systemic lupus erythematosus (SLE), but whether heterozygous C2 deficiency adds to the risk of SLE or primary Sjögren's syndrome (SS) has not been studied systematically. This study was undertaken to investigate potential associations of heterozygous C2 deficiency and C4 copy number variation with clinical manifestations in patients with SLE and patients with primary SS. METHODS: The presence of the common 28-bp C2 deletion rs9332736 and C4 copy number variation was examined in Scandinavian patients who had received a diagnosis of SLE (n = 958) or primary SS (n = 911) and in 2,262 healthy controls through the use of DNA sequencing. The concentration of complement proteins in plasma and classical complement function were analyzed in a subgroup of SLE patients. RESULTS: Heterozygous C2 deficiency-when present in combination with a low C4A copy number-substantially increased the risk of SLE (odds ratio [OR] 10.2 [95% confidence interval (95% CI) 3.5-37.0]) and the risk of primary SS (OR 13.0 [95% CI 4.5-48.4]) when compared to individuals with 2 C4A copies and normal C2. For patients heterozygous for rs9332736 with 1 C4A copy, the median age at diagnosis was 7 years earlier in patients with SLE and 12 years earlier in patients with primary SS when compared to patients with normal C2. Reduced C2 levels in plasma (P = 2 × 10-9 ) and impaired function of the classical complement pathway (P = 0.03) were detected in SLE patients with heterozygous C2 deficiency. Finally, in a primary SS patient homozygous for C2 deficiency, we observed low levels of anti-Scl-70, which suggests a risk of developing systemic sclerosis or potential overlap between primary SS and other systemic autoimmune diseases. CONCLUSION: We demonstrate that a genetic pattern involving partial deficiencies of C2 and C4A in the classical complement pathway is a strong risk factor for SLE and for primary SS. Our results emphasize the central role of the complement system in the pathogenesis of both SLE and primary SS.

Complement C4 Copy Number Variation is Linked to SSA/Ro and SSB/La Autoantibodies in Systemic Inflammatory Autoimmune Diseases.

OBJECTIVE: Copy number variation of the C4 complement components, C4A and C4B, has been associated with systemic inflammatory autoimmune diseases. This study was undertaken to investigate whether C4 copy number variation is connected to the autoimmune repertoire in systemic lupus erythematosus (SLE), primary Sjögren's syndrome (SS), or myositis. METHODS: Using targeted DNA sequencing, we determined the copy number and genetic variants of C4 in 2,290 well-characterized Scandinavian patients with SLE, primary SS, or myositis and 1,251 healthy controls. RESULTS: A prominent relationship was observed between C4A copy number and the presence of SSA/SSB autoantibodies, which was shared between the 3 diseases. The strongest association was detected in patients with autoantibodies against both SSA and SSB and 0 C4A copies when compared to healthy controls (odds ratio [OR] 18.0 [95% confidence interval (95% CI) 10.2-33.3]), whereas a weaker association was seen in patients without SSA/SSB autoantibodies (OR 3.1 [95% CI 1.7-5.5]). The copy number of C4 correlated positively with C4 plasma levels. Further, a common loss-of-function variant in C4A leading to reduced plasma C4 was more prevalent in SLE patients with a low copy number of C4A. Functionally, we showed that absence of C4A reduced the individuals' capacity to deposit C4b on immune complexes. CONCLUSION: We show that a low C4A copy number is more strongly associated with the autoantibody repertoire than with the clinically defined disease entities. These findings may have implications for understanding the etiopathogenetic mechanisms of systemic inflammatory autoimmune diseases and for patient stratification when taking the genetic profile into account.

Contribution of Rare Genetic Variation to Disease Susceptibility in a Large Scandinavian Myositis Cohort.

OBJECTIVE: Idiopathic inflammatory myopathies (IIMs) are a heterogeneous group of complex autoimmune conditions characterized by inflammation in skeletal muscle and extramuscular compartments, and interferon (IFN) system activation. We undertook this study to examine the contribution of genetic variation to disease susceptibility and to identify novel avenues for research in IIMs. METHODS: Targeted DNA sequencing was used to mine coding and potentially regulatory single nucleotide variants from ~1,900 immune-related genes in a Scandinavian case-control cohort of 454 IIM patients and 1,024 healthy controls. Gene-based aggregate testing, together with rare variant- and gene-level enrichment analyses, was implemented to explore genotype-phenotype relations. RESULTS: Gene-based aggregate tests of all variants, including rare variants, identified IFI35 as a potential genetic risk locus for IIMs, suggesting a genetic signature of type I IFN pathway activation. Functional annotation of the IFI35 locus highlighted a regulatory network linked to the skeletal muscle-specific gene PTGES3L, as a potential candidate for IIM pathogenesis. Aggregate genetic associations with AGER and PSMB8 in the major histocompatibility complex locus were detected in the antisynthetase syndrome subgroup, which also showed a less marked genetic signature of the type I IFN pathway. Enrichment analyses indicated a burden of synonymous and noncoding rare variants in IIM patients, suggesting increased disease predisposition associated with these classes of rare variants. CONCLUSION: Our study suggests the contribution of rare genetic variation to disease susceptibility in IIM and specific patient subgroups, and pinpoints genetic associations consistent with previous findings by gene expression profiling. These features highlight genetic profiles that are potentially relevant to disease pathogenesis.

DNA Methylation-Based Interferon Scores Associate With Sub-Phenotypes in Primary Sjögren's Syndrome.

Primary Sjögren's syndrome (pSS) is an autoimmune inflammatory disease with profound clinical heterogeneity, where excessive activation of the type I interferon (IFN) system is considered one of the key mechanisms in disease pathogenesis. Here we present a DNA methylation-based IFN system activation score (DNAm IFN score) and investigate its potential associations with sub-phenotypes of pSS. The study comprised 100 Swedish patients with pSS and 587 Swedish controls. For replication, 48 patients with pSS from Stavanger, Norway, were included. IFN scores were calculated from DNA methylation levels at the IFN-induced genes RSAD2, IFIT1 and IFI44L. A high DNAm IFN score, defined as > meancontrols +2SDcontrols (IFN score >4.4), was observed in 59% of pSS patients and in 4% of controls (p=1.3x10-35). Patients with a high DNAm IFN score were on average seven years younger at symptom onset (p=0.017) and at diagnosis (p=3x10-3). The DNAm IFN score levels were significantly higher in pSS positive for both SSA and SSB antibodies compared to SSA/SSB negative patients (pdiscovery=1.9x10-8, preplication=7.8x10-4). In patients positive for both SSA subtypes Ro52 and Ro60, an increased score was identified compared to single positive patients (p=0.022). Analyzing the discovery and replication cohorts together, elevated DNAm IFN scores were observed in pSS with hypergammaglobulinemia (p=2x10-8) and low C4 (p=1.5x10-3) compared to patients without these manifestations. Patients < 70 years with ongoing lymphoma at DNA sampling or lymphoma at follow-up (n=7), presented an increased DNAm IFN score compared to pSS without lymphoma (p=0.025). In conclusion, the DNAm-based IFN score is a promising alternative to mRNA-based scores for identification of patients with activation of the IFN system and may be applied for patient stratification guiding treatment decisions, monitoring and inclusion in clinical trials.

Variants in BANK1 are associated with lupus nephritis of European ancestry.

The genetic background of lupus nephritis (LN) has not been completely elucidated. We performed a case-only study of 2886 SLE patients, including 947 (33%) with LN. Renal biopsies were available from 396 patients. The discovery cohort (Sweden, n = 1091) and replication cohort 1 (US, n = 962) were genotyped on the Immunochip and replication cohort 2 (Denmark/Norway, n = 833) on a custom array. Patients with LN, proliferative nephritis, or LN with end-stage renal disease were compared with SLE without nephritis. Six loci were associated with LN (p < 1 × 10-4, NFKBIA, CACNA1S, ITGA1, BANK1, OR2Y, and ACER3) in the discovery cohort. Variants in BANK1 showed the strongest association with LN in replication cohort 1 (p = 9.5 × 10-4) and proliferative nephritis in a meta-analysis of discovery and replication cohort 1. There was a weak association between BANK1 and LN in replication cohort 2 (p = 0.052), and in the meta-analysis of all three cohorts the association was strengthened (p = 2.2 × 10-7). DNA methylation data in 180 LN patients demonstrated methylation quantitative trait loci (meQTL) effects between a CpG site and BANK1 variants. To conclude, we describe genetic variations in BANK1 associated with LN and evidence for genetic regulation of DNA methylation within the BANK1 locus. This indicates a role for BANK1 in LN pathogenesis.

Interaction between the STAT4 rs11889341(T) risk allele and smoking confers increased risk of myocardial infarction and nephritis in patients with systemic lupus erythematosus.

OBJECTIVE: To investigate how genetics influence the risk of smoking-related systemic lupus erythematosus (SLE) manifestations. METHODS: Patients with SLE (ndiscovery cohort=776, nreplication cohort=836) were genotyped using the 200K Immunochip single nucleotide polymorphisms (SNP) Array (Illumina) and a custom array. Sixty SNPs with SLE association (p<5.0×10-8) were analysed. Signal transducer and activator of transcription 4 (STAT4) activation was assessed in in vitro stimulated peripheral blood mononuclear cells from healthy controls (n=45). RESULTS: In the discovery cohort, smoking was associated with myocardial infarction (MI) (OR 1.96 (95% CI 1.09 to 3.55)), with a greater effect in patients carrying any rs11889341 STAT4 risk allele (OR 2.72 (95% CI 1.24 to 6.00)) or two risk alleles (OR 8.27 (95% CI 1.48 to 46.27)).Smokers carrying the risk allele also displayed an increased risk of nephritis (OR 1.47 (95% CI 1.06 to 2.03)). In the replication cohort, the high risk of MI in smokers carrying the risk allele and the association between the STAT4 risk allele and nephritis in smokers were confirmed (OR 6.19 (95% CI 1.29 to 29.79) and 1.84 (95% CI 1.05 to 3.29), respectively).The interaction between smoking and the STAT4 risk allele resulted in further increase in the risk of MI (OR 2.14 (95% CI 1.01 to 4.62)) and nephritis (OR 1.53 (95% CI 1.08 to 2.17)), with 54% (MI) and 34% (nephritis) of the risk attributable to the interaction. Levels of interleukin-12-induced phosphorylation of STAT4 in CD8+ T cells were higher in smokers than in non-smokers (mean geometric fluorescence intensity 1063 vs 565, p=0.0063).Lastly, the IL12A rs564799 risk allele displayed association with MI in both cohorts (OR 1.53 (95% CI 1.01 to 2.31) and 2.15 (95% CI 1.08 to 4.26), respectively). CONCLUSIONS: Smoking in the presence of the STAT4 risk gene variant appears to increase the risk of MI and nephritis in SLE. Our results also highlight the role of the IL12-STAT4 pathway in SLE-cardiovascular morbidity.

Molecular pathways in patients with systemic lupus erythematosus revealed by gene-centred DNA sequencing.

OBJECTIVES: Systemic lupus erythematosus (SLE) is an autoimmune disease with extensive heterogeneity in disease presentation between patients, which is likely due to an underlying molecular diversity. Here, we aimed at elucidating the genetic aetiology of SLE from the immunity pathway level to the single variant level, and stratify patients with SLE into distinguishable molecular subgroups, which could inform treatment choices in SLE. METHODS: We undertook a pathway-centred approach, using sequencing of immunological pathway genes. Altogether 1832 candidate genes were analysed in 958 Swedish patients with SLE and 1026 healthy individuals. Aggregate and single variant association testing was performed, and we generated pathway polygenic risk scores (PRS). RESULTS: We identified two main independent pathways involved in SLE susceptibility: T lymphocyte differentiation and innate immunity, characterised by HLA and interferon, respectively. Pathway PRS defined pathways in individual patients, who on average were positive for seven pathways. We found that SLE organ damage was more pronounced in patients positive for the T or B cell receptor signalling pathways. Further, pathway PRS-based clustering allowed stratification of patients into four groups with different risk score profiles. Studying sets of genes with priors for involvement in SLE, we observed an aggregate common variant contribution to SLE at genes previously reported for monogenic SLE as well as at interferonopathy genes. CONCLUSIONS: Our results show that pathway risk scores have the potential to stratify patients with SLE beyond clinical manifestations into molecular subsets, which may have implications for clinical follow-up and therapy selection.

Contributions of de novo variants to systemic lupus erythematosus.

By performing whole-genome sequencing in a Swedish cohort of 71 parent-offspring trios, in which the child in each family is affected by systemic lupus erythematosus (SLE, OMIM 152700), we investigated the contribution of de novo variants to risk of SLE. We found de novo single nucleotide variants (SNVs) to be significantly enriched in gene promoters in SLE patients compared with healthy controls at a level corresponding to 26 de novo promoter SNVs more in each patient than expected. We identified 12 de novo SNVs in promoter regions of genes that have been previously implicated in SLE, or that have functions that could be of relevance to SLE. Furthermore, we detected three missense de novo SNVs, five de novo insertion-deletions, and three de novo structural variants with potential to affect the expression of genes that are relevant for SLE. Based on enrichment analysis, disease-affecting de novo SNVs are expected to occur in one-third of SLE patients. This study shows that de novo variants in promoters commonly contribute to the genetic risk of SLE. The fact that de novo SNVs in SLE were enriched to promoter regions highlights the importance of using whole-genome sequencing for identification of de novo variants.

Genetic and clinical basis for two distinct subtypes of primary Sjögren's syndrome.

OBJECTIVES: Clinical presentation of primary Sjögren's syndrome (pSS) varies considerably. A shortage of evidence-based objective markers hinders efficient drug development and most clinical trials have failed to reach primary endpoints. METHODS: We performed a multicentre study to identify patient subgroups based on clinical, immunological and genetic features. Targeted DNA sequencing of 1853 autoimmune-related loci was performed. After quality control, 918 patients with pSS, 1264 controls and 107 045 single nucleotide variants remained for analysis. Replication was performed in 177 patients with pSS and 7672 controls. RESULTS: We found strong signals of association with pSS in the HLA region. Principal component analysis of clinical data distinguished two patient subgroups defined by the presence of SSA/SSB antibodies. We observed an unprecedented high risk of pSS for an association in the HLA-DQA1 locus of odds ratio 6.10 (95% CI: 4.93, 7.54, P=2.2×10-62) in the SSA/SSB-positive subgroup, while absent in the antibody negative group. Three independent signals within the MHC were observed. The two most significant variants in MHC class I and II respectively, identified patients with a higher risk of hypergammaglobulinaemia, leukopenia, anaemia, purpura, major salivary gland swelling and lymphadenopathy. Replication confirmed the association with both MHC class I and II signals confined to SSA/SSB antibody positive pSS. CONCLUSION: Two subgroups of patients with pSS with distinct clinical manifestations can be defined by the presence or absence of SSA/SSB antibodies and genetic markers in the HLA locus. These subgroups should be considered in clinical follow-up, drug development and trial outcomes, for the benefit of both subgroups.

Function of multiple sclerosis-protective HLA class I alleles revealed by genome-wide protein-quantitative trait loci mapping of interferon signalling.

Interferons (IFNs) are cytokines that are central to the host defence against viruses and other microorganisms. If not properly regulated, IFNs may contribute to the pathogenesis of inflammatory autoimmune, or infectious diseases. To identify genetic polymorphisms regulating the IFN system we performed an unbiased genome-wide protein-quantitative trait loci (pQTL) mapping of cell-type specific type I and type II IFN receptor levels and their responses in immune cells from 303 healthy individuals. Seven genome-wide significant (p < 5.0E-8) pQTLs were identified. Two independent SNPs that tagged the multiple sclerosis (MS)-protective HLA class I alleles A*02/A*68 and B*44, respectively, were associated with increased levels of IFNAR2 in B and T cells, with the most prominent effect in IgD-CD27+ memory B cells. The increased IFNAR2 levels in B cells were replicated in cells from an independent set of healthy individuals and in MS patients. Despite increased IFNAR2 levels, B and T cells carrying the MS-protective alleles displayed a reduced response to type I IFN stimulation. Expression and methylation-QTL analysis demonstrated increased mRNA expression of the pseudogene HLA-J in B cells carrying the MS-protective class I alleles, possibly driven via methylation-dependent transcriptional regulation. Together these data suggest that the MS-protective effects of HLA class I alleles are unrelated to their antigen-presenting function, and propose a previously unappreciated function of type I IFN signalling in B and T cells in MS immune-pathogenesis.

High genetic risk score is associated with early disease onset, damage accrual and decreased survival in systemic lupus erythematosus.

OBJECTIVES: To investigate associations between a high genetic disease risk and disease severity in patients with systemic lupus erythematosus (SLE). METHODS: Patients 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. RESULTS: SLE 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. CONCLUSIONS: A 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.

Shared and Unique Patterns of DNA Methylation in Systemic Lupus Erythematosus and Primary Sjögren's Syndrome.

Objectives: To perform a cross-comparative analysis of DNA methylation in patients with systemic lupus erythematosus (SLE), patients with primary Sjögren's syndrome (pSS), and healthy controls addressing the question of epigenetic sharing and aiming to detect disease-specific alterations. Methods: DNA extracted from peripheral blood from 347 cases with SLE, 100 cases with pSS, and 400 healthy controls were analyzed on the Human Methylation 450k array, targeting 485,000 CpG sites across the genome. A linear regression model including age, sex, and blood cell type distribution as covariates was fitted, and association p-values were Bonferroni corrected. A random forest machine learning classifier was designed for prediction of disease status based on DNA methylation data. Results: We established a combined set of 4,945 shared differentially methylated CpG sites (DMCs) in SLE and pSS compared to controls. In pSS, hypomethylation at type I interferon induced genes was mainly driven by patients who were positive for Ro/SSA and/or La/SSB autoantibodies. Analysis of differential methylation between SLE and pSS identified 2,244 DMCs with a majority of sites showing decreased methylation in SLE compared to pSS. The random forest classifier demonstrated good performance in discerning between disease status with an area under the curve (AUC) between 0.83 and 0.96. Conclusions: The majority of differential DNA methylation is shared between SLE and pSS, however, important quantitative differences exist. Our data highlight neutrophil dysregulation as a shared mechanism, emphasizing the role of neutrophils in the pathogenesis of systemic autoimmune diseases. The current study provides evidence for genes and molecular pathways driving common and disease-specific pathogenic mechanisms.

Genetic variations in A20 DUB domain provide a genetic link to citrullination and neutrophil extracellular traps in systemic lupus erythematosus.

OBJECTIVES: Genetic 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. METHODS: CRISPR/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. RESULTS: Genetic 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. CONCLUSIONS: We 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.

Circulating Levels of Interferon Regulatory Factor-5 Associates With Subgroups of Systemic Lupus Erythematosus Patients.

Systemic Lupus Erythematosus (SLE) is a heterogeneous autoimmune disease, which currently lacks specific diagnostic biomarkers. The diversity within the patients obstructs clinical trials but may also reflect differences in underlying pathogenesis. Our objective was to obtain protein profiles to identify potential general biomarkers of SLE and to determine molecular subgroups within SLE for patient stratification. Plasma samples from a cross-sectional study of well-characterized SLE patients (n = 379) and matched population controls (n = 316) were analyzed by antibody suspension bead array targeting 281 proteins. To investigate the differences between SLE and controls, Mann-Whitney U-test with Bonferroni correction, generalized linear modeling and receiver operating characteristics (ROC) analysis were performed. K-means clustering was used to identify molecular SLE subgroups. We identified Interferon regulating factor 5 (IRF5), solute carrier family 22 member 2 (SLC22A2) and S100 calcium binding protein A12 (S100A12) as the three proteins with the largest fold change between SLE patients and controls (SLE/Control = 1.4, 1.4, and 1.2 respectively). The lowest p-values comparing SLE patients and controls were obtained for S100A12, Matrix metalloproteinase-1 (MMP1) and SLC22A2 (padjusted = 3 × 10-9, 3 × 10-6, and 5 × 10-6 respectively). In a set of 15 potential biomarkers differentiating SLE patients and controls, two of the proteins were transcription factors, i.e., IRF5 and SAM pointed domain containing ETS transcription factor (SPDEF). IRF5 was up-regulated while SPDEF was found to be down-regulated in SLE patients. Unsupervised clustering of all investigated proteins identified three molecular subgroups among SLE patients, characterized by (1) high levels of rheumatoid factor-IgM, (2) low IRF5, and (3) high IRF5. IRF5 expressing microparticles were analyzed by flow cytometry in a subset of patients to confirm the presence of IRF5 in plasma and detection of extracellular IRF5 was further confirmed by immunoprecipitation-mass spectrometry (IP-MS). Interestingly IRF5, a known genetic risk factor for SLE, was detected extracellularly and suggested by unsupervised clustering analysis to differentiate between SLE subgroups. Our results imply a set of circulating molecules as markers of possible pathogenic importance in SLE. We believe that these findings could be of relevance for understanding the pathogenesis and diversity of SLE, as well as for selection of patients in clinical trials.

Whole-genome sequencing identifies complex contributions to genetic risk by variants in genes causing monogenic systemic lupus erythematosus.

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.