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.

Pathogenic antibody response to glucose-6-phosphate isomerase targets a modified epitope uniquely exposed on joint cartilage.

OBJECTIVES: To identify the arthritogenic B cell epitopes of glucose-6-phosphate isomerase (GPI) and their association with rheumatoid arthritis (RA). METHODS: IgG response towards a library of GPI peptides in patients with early RA, pre-symptomatic individuals and population controls, as well as in mice, were tested by bead-based multiplex immunoassays and ELISA. Monoclonal IgG were generated, and the binding specificity and affinity were determined by ELISA, gel size exclusion chromatography, surface plasma resonance and X-ray crystallography. Arthritogenicity was investigated by passive transfer experiments. Antigen-specific B cells were identified by peptide tetramer staining. RESULTS: Peptide GPI293-307 was the dominant B cell epitope in K/BxN and GPI-immunised mice. We could detect B cells and low levels of IgM antibodies binding the GPI293-307 epitopes, and high affinity anti-GPI293-307 IgG antibodies already 7 days after GPI immunisation, immediately before arthritis onset. Transfer of anti-GPI293-307 IgG antibodies induced arthritis in mice. Moreover, anti-GPI293-307 IgG antibodies were more frequent in individuals prior to RA onset (19%) than in controls (7.5%). GPI293-307-specific antibodies were associated with radiographic joint damage. Crystal structures of the Fab-peptide complex revealed that this epitope is not exposed in native GPI but requires conformational change of the protein in inflamed joint for effective recognition by anti-GPI293-307 antibodies. CONCLUSIONS: We have identified the major pathogenic B cell epitope of the RA-associated autoantigen GPI, at position 293-307, exposed only on structurally modified GPI on the cartilage surface. B cells to this neo-epitope escape tolerance and could potentially play a role in the pathogenesis of RA.

Immune complex-mediated neutrophil activation in patients with polymyalgia rheumatica.

OBJECTIVE: Neutrophils are important in host defence. However, neutrophils are also linked to inflammation and organ damage. The purpose of this study was to assess whether markers of neutrophil activation are increased in PMR. METHODS: Levels of immune complexes (IC), calprotectin and neutrophil extracellular traps (NETs) were measured in plasma of healthy individuals (n = 30) and patients with PMR (n = 60), at flare and upon treatment with glucocorticoids using ELISA. Plasma-mediated neutrophil activation was assessed in presence of an FcγRIIA inhibitory antibody (IV.3). RESULTS: Plasma levels of calprotectin and NETs were elevated in PMR (P < 0.001). Mechanistically, neutrophil activation was driven by ICs, present in plasma, able to up-regulate neutrophil activation markers CD66b and CD11b (P < 0.0001) in an FcγRIIA-dependent manner (P < 0.01). Of note, circulating levels of IC correlated with plasma induced CD66b and CD11b (r = 0.51, P = 0.004, and r = 0.46, P = 0.01, respectively) and decreased after glucocorticoid therapy. In contrast to NETs, calprotectin significantly decreased after glucocorticoid therapy (P < 0.001) and was higher in PMR without overlapping GCA compared with patients with overlapping disease (P = 0.014). Interestingly, musculoskeletal involvement was associated with elevated levels of calprotectin before initiation of glucocorticoid therapy (P = 0.036). CONCLUSIONS: Neutrophil activation, including NET formation, is increased in PMR, through IC-mediated engagement of FcγRIIA. Clinically, neutrophil activation is associated with musculoskeletal involvement, with calprotectin, but not NETs, being a biomarker of treatment response in PMR patients. In all, IC-mediated neutrophil activation is a central process in PMR pathogenesis identifying potential novel therapeutic targets (FcγRIIA), as well as soluble markers for disease monitoring (calprotectin).

Genome-wide investigation of persistence to treatment with methotrexate in early rheumatoid arthritis.

OBJECTIVES: To investigate the influence of genetic factors on persistence to treatment of early rheumatoid arthritis (RA) with methotrexate (MTX) monotherapy. METHODS: We conducted a genome-wide association study (GWAS) in a sample of 3902 Swedish early RA patients initiating MTX in DMARD-monotherapy as their first ever DMARD. The outcome, short- and long-term persistence to this treatment, was defined as remaining on MTX at one and at three years, respectively, with no additional DMARDs added. As genetic predictors, we investigated individual SNPs, and a polygenic risk score (PRS) based on SNPs associated with RA risk. The SNP-based heritability of persistence was estimated overall and by RA serostatus. RESULTS: No individual SNP reached genome-wide significance (p < 5e-8), neither for persistence at one nor at three years. The RA PRS was not significantly associated with persistence at one (RR = 0.98 (0.96-1.01)) nor three years (RR = 0.96 (0.93-1.00)). The heritability for persistence was estimated to be 0.45 (0.15-0.75) at one year and 0.14 (0-0.40) at three years. Results in seropositive RA were comparable to those in the analysis of RA overall, while heritability estimates and PRS RRs were attenuated towards the null in seronegative RA. CONCLUSIONS: Despite being the largest GWAS on an MTX treatment outcome to date, no genome-wide significant associations were detected. The modest heritability observed, coupled with the broad spread of suggestively associated loci, indicate that genetic influence is of polygenic nature. Nevertheless, persistence to MTX monotherapy was lower in patients with a greater genetic disposition, per the PRS, towards RA.

Stratified genetic analysis reveals sex differences in MPO-ANCA-associated vasculitis.

OBJECTIVE: To identify and genetically characterize subgroups of patients with ANCA-associated vasculitides (AAV) based on sex and ANCA subtype. METHODS: A previously established SNP dataset derived from DNA sequencing of 1853 genes and genotyping of 1088 Scandinavian cases with AAV and 1589 controls was stratified for sex and ANCA subtype and analysed for association with five top AAV SNPs. rs9274619, a lead variant at the HLA-DQB1/HLA-DQA2 locus previously associated with AAV positive for myeloperoxidase (MPO)-ANCA, was analysed for association with the cumulative disease involvement of ten different organ systems. RESULTS: rs9274619 showed a significantly stronger association to MPO-ANCA-positive females than males [P = 2.0 × 10-4, OR = 2.3 (95% CI 1.5, 3.5)], whereas proteinase 3 (PR3)-ANCA-associated variants rs1042335, rs9277341 (HLA-DPB1/A1) and rs28929474 (SERPINA1) were equally associated with females and males with PR3-ANCA. In MPO-ANCA-positive cases, carriers of the rs9274619 risk allele were more prone to disease engagement of eyes [P = 0.021, OR = 11 (95% CI 2.2, 205)] but less prone to pulmonary involvement [P = 0.026, OR = 0.52 (95% CI 0.30, 0.92)]. Moreover, AAV with both MPO-ANCA and PR3-ANCA was associated with the PR3-ANCA lead SNP rs1042335 [P = 0.0015, OR = 0.091 (95% CI 0.0022, 0.55)] but not with rs9274619. CONCLUSIONS: Females and males with MPO-ANCA-positive AAV differ in genetic predisposition to disease, suggesting at least partially distinct disease mechanisms between the sexes. Double ANCA-positive AAV cases are genetically similar to PR3-ANCA-positive cases, providing clues to the clinical follow-up and treatment of these patients.

Pulmonary fibrosis in relation to genetic loci in an inception cohort of patients with early rheumatoid arthritis from northern Sweden.

OBJECTIVES: Pulmonary manifestations in RA are common comorbidities. Interstitial lung disease (ILD), both idiopathic and in RA, has been associated with several genetic variants. We assessed pulmonary fibrosis (PF) in an inception cohort of RA patients in relation to genetic variants and disease-related factors. METHODS: A total of 1466 early RA patients were consecutively included and followed prospectively from the index date until death or 31 December 2016. Clinical and laboratory data and treatment were continuously registered according to the Swedish Rheumatology Quality Register. DNA was available from 1184 patients and 571 151 genome-wide single-nucleotide polymorphisms (SNPs) were analysed. Thirteen identified genetic variants were extracted. At follow-up, the patients answered a questionnaire regarding disease progression and lung involvement that was validated by reviewing medical records and analysing radiological examinations. RESULTS: The prevalence of PF was 5.6% and the annualized incidence rate was 5.0/1000 (95% CI 3.80, 6.54). Four SNPs were associated with PF in RA: rs35705950 [MUC5B; OR 2.5 (95% CI 1.5, 4.0), adjusted P-value = 0.00016, q-value = 0.0021]; rs111521887 [TOLLIP; OR 1.9 (95% CI 1.3, 2.8), adjusted P-value = 0.0014, q-value = 0.0092]; rs2609255 [FAM13A; OR 1.7 (95% CI 1.1, 2.5), adjusted P-value = 0.013, q-value = 0.055] and rs2736100 [TERT; OR 1.5 (95% CI 1.0, 2.2), adjusted P-value = 0.046, q-value = 0.15]. Older age and RF positivity were associated with increased risk, while MTX treatment was associated with a lower risk of PF. CONCLUSIONS: Development of PF in an inception cohort of RA patients was associated with 4 of 12 ILD risk genes. RA-related factors except for age at diagnosis and RF positivity were of limited importance in PF development.

Anti-citrullinated protein antibody specificities and pulmonary fibrosis in relation to genetic loci in early rheumatoid arthritis.

OBJECTIVES: Pulmonary manifestations in RA are common comorbidities, but the underlying mechanisms are largely unknown. The added value of a multiplex of ACPA and genetic risk markers was evaluated for the development of pulmonary fibrosis (PF) in an inception cohort. METHODS: A total of 1184 patients with early RA were consecutively included and followed prospectively from the index date until death or 31 December 2016. The presence of 21 ACPA fine specificities was analysed using a custom-made microarray chip (Thermo Fisher Scientific, Uppsala, Sweden). Three SNPs, previously found related to PF were evaluated, rs2609255 (FAM13A), rs111521887 (TOLLIP) and rs35705950 (MUC5B). ACPA and genetic data were available for 841 RA patients, of whom 50 developed radiologically defined PF. RESULTS: In unadjusted analyses, 11 ACPA specificities were associated with PF development. In multiple variable analyses, six ACPA specificities were associated with increased risk of PF: vimentin (Vim)60-75, fibrinogen (Fib)ß62-78 (72), Fibα621-635, Bla26, collagen (C)II359-369 and F4-CIT-R (P < 0.01 to P < 0.05). The number of ACPA specificities was also related to PF development (P < 0.05 crude and adjusted models). In multiple variable models respectively adjusted for each of the SNPs, the number of ACPA specificities (P < 0.05 in all models), anti-Vim60-75 (P < 0.05, in all models), anti-Fibß62-78 (72) (P < 0.001 to P < 0.05), anti-CII359-369 (P < 0.05 in all models) and anti-F4-CIT-R AQ4 (P < 0.01 to P < 0.05), anti-Fibα621-635 (P < 0.05 in one) and anti-Bla26 (P < 0.05 in two) were significantly associated with PF development. CONCLUSION: The development of PF in an inception cohort of RA patients was associated with both presence of certain ACPA and the number of ACPA specificities and risk genes.

Identification and functional characterization of a novel susceptibility locus for small vessel vasculitis with MPO-ANCA.

OBJECTIVE: To identify and characterize genetic loci associated with the risk of developing ANCA-associated vasculitides (AAV). METHODS: Genetic association analyses were performed after Illumina sequencing of 1853 genes and subsequent replication with genotyping of selected single nucleotide polymorphisms in a total cohort of 1110 Scandinavian cases with granulomatosis with polyangiitis or microscopic polyangiitis, and 1589 controls. A novel AAV-associated single nucleotide polymorphism was analysed for allele-specific effects on gene expression using luciferase reporter assay. RESULTS: PR3-ANCA+ AAV was significantly associated with two independent loci in the HLA-DPB1/HLA-DPA1 region [rs1042335, P = 6.3 × 10-61, odds ratio (OR) 0.10; rs9277341, P = 1.5 × 10-44, OR 0.22] and with rs28929474 in the SERPINA1 gene (P = 2.7 × 10-10, OR 2.9). MPO-ANCA+ AAV was significantly associated with the HLA-DQB1/HLA-DQA2 locus (rs9274619, P = 5.4 × 10-25, OR 3.7) and with a rare variant in the BACH2 gene (rs78275221, P = 7.9 × 10-7, OR 3.0), the latter a novel susceptibility locus for MPO-ANCA+ granulomatosis with polyangiitis/microscopic polyangiitis. The rs78275221-A risk allele reduced luciferase gene expression in endothelial cells, specifically, as compared with the non-risk allele. CONCLUSION: We identified a novel susceptibility locus for MPO-ANCA+ AAV and propose that the associated variant is of mechanistic importance, exerting a regulatory function on gene expression in specific cell types.

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.

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.

Anti-neutrophil cytoplasmic antibodies predate symptom onset of ANCA-associated vasculitis. A case-control study.

OBJECTIVES: Anti-neutrophil cytoplasmic autoantibodies [ANCA) are important for diagnosis of ANCA-associated vasculitides (AAV). The timing of antibody development is not well established. To investigate the development of proteinase 3 (PR3)- and myeloperoxidase (MPO)-ANCA, blood samples collected before onset of symptoms of AAV were analysed. METHODS: To identify AAV patients with blood samples predating symptoms, the National Patient Register and Cause of Death register were scrutinized for ICD codes for AAV and linked to the registers of five biobanks. Diagnoses of AAV and time point for symptom onset were confirmed by reviewing 504 case-record. Eighty-five AAV cases (34 males, 51 females) with samples >1 month < 10 years from AAV symptom onset and two controls matched for sex, age, and sampling time for each case were included. Samples were screened using ELISAs for ANCA and further analysed for PR3-or MPO- specificities. RESULTS: In ANCA-screen 35.7% of the pre-symptomatic cases and 3.5% of controls tested positive (p < 0.01). 26.2% of the cases were PR3-ANCA+ and 10.7% MPO-ANCA+. Median (Q1-Q3) predating time for PR3-ANCA+ was 2.7 (0.3-7.7) years and MPO-ANCA+ 2.0 (0.9-3.5) years. PR3-ANCA was demonstrated in samples up to nine years before symptom onset. At symptom onset predating PR3-ANCA+ cases were younger than PR3-ANCA- (P < 0.01), and MPO-ANCA+ were older than MPO-ANCA- (p < 0.05). Predating MPO-ANCA+ cases vs. MPO-ANCA- and vs. PR3-ANCA+ cases had more often at symptoms onset manifestations from lungs, kidneys or peripheral nervous system (p < 0.01 and p < 0.05, respectively). CONCLUSION: The PR3-and MPO-ANCAs are present years before AAV symptom onset and represent distinct diseases.

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.

Osteoprotegerin and osteocalcin are associated with atherosclerosis in patients with rheumatoid arthritis: a prospective cohort study.

OBJECTIVES: Patients with rheumatoid arthritis (RA) have an accelerated progression of atherosclerosis. The aim of this study was to examine the associations between subclinical atherosclerosis, assessed by intima-media thickness (IMT), and regulators of bone formation, markers of bone turnover and bone mineral density (BMD) in patients with RA. METHODS: Patients with new-onset RA (n=79), aged ≤60 years at diagnosis, were consecutively included in a study of development of atherosclerosis. Ultrasound measurement of IMT of the common carotid artery was undertaken at inclusion (T0) and after 11 years (T11) (n=54). Bone turnover biomarkers were examined in samples collected at T0 and T11. BMD was assessed at T11. RESULTS: In patients with RA, osteocalcin (OCN) and osteoprotegerin (OPG) measured at T11 were significantly associated with IMT at T11, adjusted for systolic blood pressure (SBP) and age. BMD at T11 and the bone turnover markers procollagen type 1 N-terminal propeptide (P1NP) and carboxy-terminal crosslinked C-terminal telopeptide (CTX) were not associated with IMT. OPG, OCN and sclerostin at T0 were significantly associated with IMT at T11, and OPG and OCN at T0 were associated with change in IMT from T0 to T11. The associations between IMT and bone biomarkers were stronger in patients with joint erosions at onset of RA, than in patients with non-erosive disease. CONCLUSIONS: Atherosclerosis in patients with RA is associated with OPG and OCN, but not with BMD or markers reflecting ongoing bone turnover, indicating that atherosclerosis is not associated with bone turnover per se.

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.

Smoking cessation is associated with lower disease activity and predicts cardiovascular risk reduction in rheumatoid arthritis patients.

OBJECTIVES: Smoking is a major risk factor for the development of both cardiovascular disease (CVD) and RA and may cause attenuated responses to anti-rheumatic treatments. Our aim was to compare disease activity, CVD risk factors and CVD event rates across smoking status in RA patients. METHODS: Disease characteristics, CVD risk factors and relevant medications were recorded in RA patients without prior CVD from 10 countries (Norway, UK, Netherlands, USA, Sweden, Greece, South Africa, Spain, Canada and Mexico). Information on CVD events was collected. Adjusted analysis of variance, logistic regression and Cox models were applied to compare RA disease activity (DAS28), CVD risk factors and event rates across categories of smoking status. RESULTS: Of the 3311 RA patients (1012 former, 887 current and 1412 never smokers), 235 experienced CVD events during a median follow-up of 3.5 years (interquartile range 2.5-6.1). At enrolment, current smokers were more likely to have moderate or high disease activity compared with former and never smokers (P < 0.001 for both). There was a gradient of worsening CVD risk factor profiles (lipoproteins and blood pressure) from never to former to current smokers. Furthermore, former and never smokers had significantly lower CVD event rates compared with current smokers [hazard ratio 0.70 (95% CI 0.51, 0.95), P = 0.02 and 0.48 (0.34, 0.69), P < 0.001, respectively]. The CVD event rates for former and never smokers were comparable. CONCLUSION: Smoking cessation in patients with RA was associated with lower disease activity and improved lipid profiles and was a predictor of reduced rates of CVD events.

On the presence of HLA-SE alleles and ACPA-IgG variable domain glycosylation in the phase preceding the development of rheumatoid arthritis.

OBJECTIVE: Anti-citrullinated protein antibodies (ACPA) in rheumatoid arthritis (RA) patients display a unique feature defined by the abundant presence of N-linked glycans within the variable domains (V-domains). Recently, we showed that N-glycosylation sites, which are required for the incorporation of V-domain glycans, are introduced following somatic hypermutation. However, it is currently unclear when V-domain glycosylation occurs. Further, it is unknown which factors might trigger the generation of V-domain glycans and whether such glycans are relevant for the transition towards RA. Here, we determined the presence of ACPA-IgG V-domain glycans in paired samples of pre-symptomatic individuals and RA patients. METHODS: ACPA-IgG V-domain glycosylation was analysed using ultra-high performance liquid chromatography (UHPLC) in paired samples of pre-symptomatic individuals (median interquartile range (IQR) pre-dating time: 5.8 (5.9) years; n=201; 139 ACPA-positive and 62 ACPA-negative) and RA patients (n=99; 94 ACPA-positive and 5 ACPA-negative). RESULTS: V-domain glycans on ACPA-IgG were already present up to 15 years before disease in pre-symptomatic individuals and their abundance increased closer to symptom onset. Noteworthy, human leucocyte antigen class II shared epitope (HLA-SE) alleles associated with the presence of V-domain glycans on ACPA-IgG. CONCLUSION: Our observations indicate that somatic hypermutation of ACPA, which results in the incorporation of N-linked glycosylation sites and consequently V-domain glycans, occurs already years before symptom onset in individuals that will develop RA later in life. Moreover, our findings provide first evidence that HLA-SE alleles associate with ACPA-IgG V-domain glycosylation in the pre-disease phase and thereby further refine the connection between HLA-SE and the development of ACPA-positive RA.

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.

DNA methylation mapping identifies gene regulatory effects in patients with systemic lupus erythematosus.

OBJECTIVES: Systemic 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. METHODS: DNA 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. RESULTS: We 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. CONCLUSIONS: Our 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.

Novel gene variants associated with cardiovascular disease in systemic lupus erythematosus and rheumatoid arthritis.

OBJECTIVES: Patients with systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) have increased risk of cardiovascular disease (CVD). We investigated whether single nucleotide polymorphisms (SNPs) at autoimmunity risk loci were associated with CVD in SLE and RA. METHODS: Patients with SLE (n=1045) were genotyped using the 200K Immunochip SNP array (Illumina). The allele frequency was compared between patients with and without different manifestations of CVD. Results were replicated in a second SLE cohort (n=1043) and in an RA cohort (n=824). We analysed publicly available genetic data from general population, performed electrophoretic mobility shift assays and measured cytokine levels and occurrence of antiphospholipid antibodies (aPLs). RESULTS: We identified two new putative risk loci associated with increased risk for CVD in two SLE populations, which remained after adjustment for traditional CVD risk factors. An IL19 risk allele, rs17581834(T) was associated with stroke/myocardial infarction (MI) in SLE (OR 2.3 (1.5 to 3.4), P=8.5×10-5) and RA (OR 2.8 (1.4 to 5.6), P=3.8×10-3), meta-analysis (OR 2.5 (2.0 to 2.9), P=3.5×10-7), but not in population controls. The IL19 risk allele affected protein binding, and SLE patients with the risk allele had increased levels of plasma-IL10 (P=0.004) and aPL (P=0.01). An SRP54-AS1 risk allele, rs799454(G) was associated with stroke/transient ischaemic attack in SLE (OR 1.7 (1.3 to 2.2), P=2.5×10-5) but not in RA. The SRP54-AS1 risk allele is an expression quantitative trait locus for four genes. CONCLUSIONS: The IL19 risk allele was associated with stroke/MI in SLE and RA, but not in the general population, indicating that shared immune pathways may be involved in the CVD pathogenesis in inflammatory rheumatic diseases.

The anti-carbamylated protein antibody response is of overall low avidity despite extensive isotype switching.

Objective: To better understand the contribution of autoantibodies in RA and the biology of their responses, we evaluated the avidity of the anti-carbamylated protein (anti-CarP) antibody response. Methods: The avidity of anti-CarP antibody, ACPA and anti-tetanus toxoid IgG were determined using elution assays. Anti-CarP IgG avidity was measured in sera of 107 RA patients, 15 paired SF and serum samples and 8 serially sampled sera before and after disease onset. Results: The avidity of anti-CarP IgG is low compared with the avidity of anti-tetanus toxoid IgG present in the same sera. Likewise, although less pronounced, anti-CarP also displayed a lower avidity as compared with the avidity of ACPA IgG. No difference in anti-CarP IgG avidity is observed between ACPA positive or ACPA negative patients. Anti-CarP IgG avidity is higher in anti-CarP IgM-negative compared with IgM-positive individuals. Furthermore, the anti-CarP avidity in serum is higher than in SF. Using samples of individuals that over time developed RA we observed no anti-CarP avidity maturation in the years before disease onset. In contrast to ACPA avidity, the anti-CarP avidity is not associated with severity of joint destruction. Conclusion: The anti-CarP response is of overall low avidity, even lower than the ACPA IgG avidity, and does not show apparent avidity maturation before or around disease onset. Overall, isotype switch and avidity maturation seem to be uncoupled as isotype switch occurs without avidity maturation, pointing towards a commonality in the regulation of both autoantibody responses as opposed to the pathways governing recall responses.