Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets.

Circulating proteins have important functions in inflammation and a broad range of diseases. To identify genetic influences on inflammation-related proteins, we conducted a genome-wide protein quantitative trait locus (pQTL) study of 91 plasma proteins measured using the Olink Target platform in 14,824 participants. We identified 180 pQTLs (59 cis, 121 trans). Integration of pQTL data with eQTL and disease genome-wide association studies provided insight into pathogenesis, implicating lymphotoxin-α in multiple sclerosis. Using Mendelian randomization (MR) to assess causality in disease etiology, we identified both shared and distinct effects of specific proteins across immune-mediated diseases, including directionally discordant effects of CD40 on risk of rheumatoid arthritis versus multiple sclerosis and inflammatory bowel disease. MR implicated CXCL5 in the etiology of ulcerative colitis (UC) and we show elevated gut CXCL5 transcript expression in patients with UC. These results identify targets of existing drugs and provide a powerful resource to facilitate future drug target prioritization.

FLT3 stop mutation increases FLT3 ligand level and risk of autoimmune thyroid disease.

Autoimmune thyroid disease is the most common autoimmune disease and is highly heritable1. Here, by using a genome-wide association study of 30,234 cases and 725,172 controls from Iceland and the UK Biobank, we find 99 sequence variants at 93 loci, of which 84 variants are previously unreported2-7. A low-frequency (1.36%) intronic variant in FLT3 (rs76428106-C) has the largest effect on risk of autoimmune thyroid disease (odds ratio (OR) = 1.46, P = 2.37 × 10-24). rs76428106-C is also associated with systemic lupus erythematosus (OR = 1.90, P = 6.46 × 10-4), rheumatoid factor and/or anti-CCP-positive rheumatoid arthritis (OR = 1.41, P = 4.31 × 10-4) and coeliac disease (OR = 1.62, P = 1.20 × 10-4). FLT3 encodes fms-related tyrosine kinase 3, a receptor that regulates haematopoietic progenitor and dendritic cells. RNA sequencing revealed that rs76428106-C generates a cryptic splice site, which introduces a stop codon in 30% of transcripts that are predicted to encode a truncated protein, which lacks its tyrosine kinase domains. Each copy of rs76428106-C doubles the plasma levels of the FTL3 ligand. Activating somatic mutations in FLT3 are associated with acute myeloid leukaemia8 with a poor prognosis and rs76428106-C also predisposes individuals to acute myeloid leukaemia (OR = 1.90, P = 5.40 × 10-3). Thus, a predicted loss-of-function germline mutation in FLT3 causes a reduction in full-length FLT3, with a compensatory increase in the levels of its ligand and an increased disease risk, similar to that of a gain-of-function mutation.

Association of rheumatoid factor, anti-citrullinated protein antibodies and shared epitope with clinical response to initial treatment in patients with early rheumatoid arthritis: data from a randomised controlled trial.

OBJECTIVES: To investigate whether rheumatoid factor (RF), anti-citrullinated protein antibodies (ACPAs) and shared epitope (SE) allele-related genetic markers associate with treatment response to abatacept, certolizumab pegol or tocilizumab versus active conventional treatment (ACT). METHODS: Patients with treatment-naïve early rheumatoid arthritis were randomised in the NORD-STAR trial to ACT, certolizumab pegol, abatacept or tocilizumab, all with methotrexate. Centralised laboratory analyses for ACPA, RF and SE were performed. Clinical Disease Activity Index remission was analysed longitudinally with logistic generalised estimating equations. Differences in treatment effect across RF, ACPA and SE subgroups were assessed with interaction terms at 24 and 48 weeks, adjusted for sex, country, age, body mass index, Disease Activity Score of 28 joints based on C-reactive protein and smoking. RESULTS: In total, 778 patients were included. At 24 weeks, abatacept treatment showed a better response than ACT in the RF and/or ACPA-positive subgroups, but this effect was not significantly different from the negative subgroups. By 48 weeks, abatacept treatment showed better response regardless of RF/ACPA status. No differences were found across RF, ACPA, SE allele, valine at amino acid position 11 or valine-arginine-alanine haplotype subgroups for any biological treatment at 48 weeks. CONCLUSIONS: Based on this randomised controlled trial, abatacept treatment was associated with a better response than ACT in the RF and/or ACPA-positive subgroup at 24 weeks, but this was no longer seen at 48 weeks; adding SE allele-related genetic markers did not strengthen the association. Moreover, ACPA, RF and SE allele-related genotypes were not, alone or in combination, associated with clinical responses of importance sufficiently strongly to warrant implementation in clinical practice. TRIAL REGISTRATION NUMBER: EudraCT 2011-004720-35; ClinicalTrials.gov NCT01491815.

Identification of early risk factors for anti-citrullinated-protein-antibody positive rheumatoid arthritis-a prospective cohort study.

OBJECTIVE: Individuals positive for anti-cyclic-peptide-antibodies (anti-CCP) and musculoskeletal complaints (MSK-C) are at risk for developing rheumatoid arthritis (RA). In this study we aimed to investigate factors involved in arthritis progression. METHODS: Anti-CCP2-positive individuals with MSK-C referred to a rheumatologist were recruited. Individuals lacked arthritis at clinical and ultrasound examination and were followed for ≥three years or until clinical arthritis diagnosis. Blood samples from inclusion were analyzed for; nine anti-citrullinated-protein-antibody (ACPA) reactivities (citrullinated α-1-enolase, fibrinogen, filaggrin, histone, vimentin and tenascin peptides); 92 inflammation-associated proteins; and HLA-shared epitope alleles. Cox regression was applied to the data to identify independent predictors in a model. RESULTS: 267 individuals were included with median follow up of 49 months (IQR: 22-60). 101 (38%) developed arthritis after median 14 months (IQR: 6-27). The analysis identified that presence of at least one ACPA reactivity (HR 8.0, 95% CI 2.9-22), ultrasound detected tenosynovitis (HR 3.4, 95% CI 2.0-6.0), IL6 levels (HR 1.5, 95% CI 1.2-1.8) and IL15-Rα levels (HR 0.6, 95% CI 0.4-0.9) are significant independent predictors for arthritis progression in a prediction model (Harrell's C 0.76 [SE 0.02], AUC 0.82 [95% CI 0.76-0.89], cross-validated AUC 0.70 [95% CI 0.56-0.85]). CONCLUSION: We propose a high-Risk-RA phase characterized by presence of ACPA reactivity, tenosynovitis, IL6, and IL15-Rα and suggest that these factors need to be further investigated for their biological effects and clinical values, to identify individuals at particular low risk and high risk for arthritis progression.

Occupational inhalable agents constitute major risk factors for rheumatoid arthritis, particularly in the context of genetic predisposition and smoking.

OBJECTIVES: To assess the effects of occupational inhalable exposures on rheumatoid arthritis (RA) development and their interactions with smoking and RA-risk genes, stratifying by presence of anticitrullinated protein antibodies (ACPA). METHODS: Data came from the Swedish Epidemiological Investigation of RA, consisting of 4033 incident RA cases and 6485 matched controls. Occupational histories were retrieved, combining with a Swedish national job-exposure matrix, to estimate exposure to 32 inhalable agents. Genetic data were used to define Genetic Risk Score (GRS) or carrying any copy of human leucocyte antigen class II shared epitope (HLA-SE) alleles. Associations were identified with unconditional logistical regression models. Attributable proportion due to interaction was estimated to evaluate presence of interaction. RESULTS: Exposure to any occupational inhalable agents was associated with increased risk for ACPA-positive RA (OR 1.25, 95% CI 1.12 to 1.38). The risk increased as number of exposed agents increased (Ptrend<0.001) or duration of exposure elongated (Ptrend<0.001). When jointly considering exposure to any occupational inhalable agents, smoking and high GRS, a markedly elevated risk for ACPA-positive RA was observed among the triple-exposed group compared with those not exposed to any (OR 18.22, 95% CI 11.77 to 28.19). Significant interactions were found between occupational inhalable agents and smoking/genetic factors (high GRS or HLA-SE) in ACPA-positive RA. CONCLUSIONS: Occupational inhalable agents could act as important environmental triggers in RA development and interact with smoking and RA-risk genes leading to excessive risk for ACPA-positive RA. Future studies are warranted to assess preventive strategies aimed at reducing occupational hazards and smoking, especially among those who are genetically vulnerable.

Low copy numbers of complement C4 and C4A deficiency are risk factors for myositis, its subgroups and autoantibodies.

BACKGROUND: Idiopathic inflammatory myopathies (IIM) are a group of autoimmune diseases characterised by myositis-related autoantibodies plus infiltration of leucocytes into muscles and/or the skin, leading to the destruction of blood vessels and muscle fibres, chronic weakness and fatigue. While complement-mediated destruction of capillary endothelia is implicated in paediatric and adult dermatomyositis, the complex diversity of complement C4 in IIM pathology was unknown. METHODS: We elucidated the gene copy number (GCN) variations of total C4, C4A and C4B, long and short genes in 1644 Caucasian patients with IIM, plus 3526 matched healthy controls using real-time PCR or Southern blot analyses. Plasma complement levels were determined by single radial immunodiffusion. RESULTS: The large study populations helped establish the distribution patterns of various C4 GCN groups. Low GCNs of C4T (C4T=2+3) and C4A deficiency (C4A=0+1) were strongly correlated with increased risk of IIM with OR equalled to 2.58 (2.28-2.91), p=5.0×10-53 for C4T, and 2.82 (2.48-3.21), p=7.0×10-57 for C4A deficiency. Contingency and regression analyses showed that among patients with C4A deficiency, the presence of HLA-DR3 became insignificant as a risk factor in IIM except for inclusion body myositis (IBM), by which 98.2% had HLA-DR3 with an OR of 11.02 (1.44-84.4). Intragroup analyses of patients with IIM for C4 protein levels and IIM-related autoantibodies showed that those with anti-Jo-1 or with anti-PM/Scl had significantly lower C4 plasma concentrations than those without these autoantibodies. CONCLUSIONS: C4A deficiency is relevant in dermatomyositis, HLA-DRB1*03 is important in IBM and both C4A deficiency and HLA-DRB1*03 contribute interactively to risk of polymyositis.

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.

Repurposing antidiabetic drugs for rheumatoid arthritis: results from a two-sample Mendelian randomization study.

Despite increasing therapeutic options to treat rheumatoid arthritis (RA), many patients fail to reach treatment targets. The use of antidiabetic drugs like thiazolidinediones has been associated with lower RA risk. We aimed to explore the repurposing potential of antidiabetic drugs in RA prevention by assessing associations between genetic variation in antidiabetic drug target genes and RA using Mendelian randomization (MR). A two-sample MR design was used to estimate the association between the antidiabetic drug and RA risk using summary statistics from genome-wide association studies (GWAS). We selected independent genetic variants from the gene(s) that encode the target protein(s) of the investigated antidiabetic drug as instruments. We extracted the associations of instruments with blood glucose concentration and RA from the UK Biobank and a GWAS meta-analysis of clinically diagnosed RA, respectively. The effect of genetic variation in the drug target(s) on RA risk was estimated by the Wald ratio test or inverse-variance weighted method. Insulin and its analogues, thiazolidinediones, and sulfonylureas had valid genetic instruments (n = 1, 1, and 2, respectively). Genetic variation in thiazolidinedione target (gene: PPARG) was inversely associated with RA risk (odds ratio [OR] 0.38 per 0.1mmol/L glucose lowering, 95% confidence interval [CI] 0.20-0.73). Corresponding ORs (95%CIs) were 0.83 (0.44-1.55) for genetic variation in the targets of insulin and its analogues (gene: INSR), and 1.12 (0.83, 1.49) 1.25 (0.78-2.00) for genetic variation in the sulfonylurea targets (gene: ABCC8 and KCNJ11). In conclusion, genetic variation in the thiazolidinedione target is associated with a lower RA risk. The underlying mechanisms warrant further exploration.

Distinct HLA Associations with Rheumatoid Arthritis Subsets Defined by Serological Subphenotype.

Rheumatoid arthritis (RA) is the most common immune-mediated arthritis. Anti-citrullinated peptide antibodies (ACPA) are highly specific to RA and assayed with the commercial CCP2 assay. Genetic drivers of RA within the MHC are different for CCP2-positive and -negative subsets of RA, particularly at HLA-DRB1. However, aspartic acid at amino acid position 9 in HLA-B (Bpos-9) increases risk to both RA subsets. Here we explore how individual serologies associated with RA drive associations within the MHC. To define MHC differences for specific ACPA serologies, we quantified a total of 19 separate ACPAs in RA-affected case subjects from four cohorts (n = 6,805). We found a cluster of tightly co-occurring antibodies (canonical serologies, containing CCP2), along with several independently expressed antibodies (non-canonical serologies). After imputing HLA variants into 6,805 case subjects and 13,467 control subjects, we tested associations between the HLA region and RA subgroups based on the presence of canonical and/or non-canonical serologies. We examined CCP2(+) and CCP2(-) RA-affected case subjects separately. In CCP2(-) RA, we observed that the association between CCP2(-) RA and Bpos-9 was derived from individuals who were positive for non-canonical serologies (omnibus_p = 9.2 × 10-17). Similarly, we observed in CCP2(+) RA that associations between subsets of CCP2(+) RA and Bpos-9 were negatively correlated with the number of positive canonical serologies (p = 0.0096). These findings suggest unique genetic characteristics underlying fine-specific ACPAs, suggesting that RA may be further subdivided beyond simply seropositive and seronegative.

Early prediction of clinical response to anti-TNF treatment using multi-omics and machine learning in rheumatoid arthritis.

OBJECTIVES: Advances in immunotherapy by blocking TNF have remarkably improved treatment outcomes for Rheumatoid arthritis (RA) patients. Although treatment specifically targets TNF, the downstream mechanisms of immune suppression are not completely understood. The aim of this study was to detect biomarkers and expression signatures of treatment response to TNF inhibition. METHODS: Peripheral blood mononuclear cells (PBMCs) from 39 female patients were collected before anti-TNF treatment initiation (day 0) and after 3 months. The study cohort included patients previously treated with MTX who failed to respond adequately. Response to treatment was defined based on the EULAR criteria and classified 23 patients as responders and 16 as non-responders. We investigated differences in gene expression in PBMCs, the proportion of cell types and cell phenotypes in peripheral blood using flow cytometry and the level of proteins in plasma. Finally, we used machine learning models to predict non-response to anti-TNF treatment. RESULTS: The gene expression analysis in baseline samples revealed notably higher expression of the gene EPPK1 in future responders. We detected the suppression of genes and proteins following treatment, including suppressed expression of the T cell inhibitor gene CHI3L1 and its protein YKL-40. The gene expression results were replicated in an independent cohort. Finally, machine learning models mainly based on transcriptomic data showed high predictive utility in classifying non-response to anti-TNF treatment in RA. CONCLUSIONS: Our integrative multi-omics analyses identified new biomarkers for the prediction of response, found pathways influenced by treatment and suggested new predictive models of anti-TNF treatment in RA patients.

Molecular mimicry between Anoctamin 2 and Epstein-Barr virus nuclear antigen 1 associates with multiple sclerosis risk.

Multiple sclerosis (MS) is a chronic inflammatory, likely autoimmune disease of the central nervous system with a combination of genetic and environmental risk factors, among which Epstein-Barr virus (EBV) infection is a strong suspect. We have previously identified increased autoantibody levels toward the chloride-channel protein Anoctamin 2 (ANO2) in MS. Here, IgG antibody reactivity toward ANO2 and EBV nuclear antigen 1 (EBNA1) was measured using bead-based multiplex serology in plasma samples from 8,746 MS cases and 7,228 controls. We detected increased anti-ANO2 antibody levels in MS (P = 3.5 × 10-36) with 14.6% of cases and 7.8% of controls being ANO2 seropositive (odds ratio [OR] = 1.6; 95% confidence intervals [95%CI]: 1.5 to 1.8). The MS risk increase in ANO2-seropositive individuals was dramatic when also exposed to 3 known risk factors for MS: HLA-DRB1*15:01 carriage, absence of HLA-A*02:01, and high anti-EBNA1 antibody levels (OR = 24.9; 95%CI: 17.9 to 34.8). Reciprocal blocking experiments with ANO2 and EBNA1 peptides demonstrated antibody cross-reactivity, mapping to ANO2 [aa 140 to 149] and EBNA1 [aa 431 to 440]. HLA gene region was associated with anti-ANO2 antibody levels and HLA-DRB1*04:01 haplotype was negatively associated with ANO2 seropositivity (OR = 0.6; 95%CI: 0.5 to 0.7). Anti-ANO2 antibody levels were not increased in patients from 3 other inflammatory disease cohorts. The HLA influence and the fact that specific IgG production usually needs T cell help provides indirect evidence for a T cell ANO2 autoreactivity in MS. We propose a hypothesis where immune reactivity toward EBNA1 through molecular mimicry with ANO2 contributes to the etiopathogenesis of MS.

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.

HLA-B*27 is significantly enriched in Nordic patients with psoriatic arthritis mutilans.

OBJECTIVES: The genetic contribution to psoriatic disease is substantial with a dominating influence of the HLA region. The profile of HLA class I genotypes likely contributes to shaping clinical phenotypes. Herein we aimed to explore such genotypes in cohorts of closely characterised subsets of psoriatic disease with special focus on psoriatic arthritis mutilans (PAM), a severe and rare form of psoriatic arthritis (PsA). METHODS: Cohorts of patients with the diagnosis of psoriasis vulgaris with or without arthritis (n=1217), psoriasis without arthritis (n=534), psoriatic arthritis without mutilating disease (n=337) and psoriatic arthritis mutilans (n=63) were collected and genotyped for HLA class I and II genes, with standardised methodologies. Cases were compared with a healthy control population (n=2468). Case-only and case-control association tests were performed to address the hypothesis of genetic contribution to clinical phenotypes. RESULTS: The presence of HLA-B*27 was strikingly increased in PAM (45%) compared with PsA without mutilating disease (13%) and with healthy controls (13%). However, within the PAM population, HLA-B*27 did not correlate with clinical markers such as number of mutilating joints, radiographic scoring, disease duration and age of disease onset. CONCLUSIONS: HLA-B*27 emerges as an important genotype marker for PAM.

Genetic overlap between autoimmune diseases and non-Hodgkin lymphoma subtypes.

Epidemiologic studies show an increased risk of non-Hodgkin lymphoma (NHL) in patients with autoimmune disease (AD), due to a combination of shared environmental factors and/or genetic factors, or a causative cascade: chronic inflammation/antigen-stimulation in one disease leads to another. Here we assess shared genetic risk in genome-wide-association-studies (GWAS). Secondary analysis of GWAS of NHL subtypes (chronic lymphocytic leukemia, diffuse large B-cell lymphoma, follicular lymphoma, and marginal zone lymphoma) and ADs (rheumatoid arthritis, systemic lupus erythematosus, and multiple sclerosis). Shared genetic risk was assessed by (a) description of regional genetic of overlap, (b) polygenic risk score (PRS), (c)"diseasome", (d)meta-analysis. Descriptive analysis revealed few shared genetic factors between each AD and each NHL subtype. The PRS of ADs were not increased in NHL patients (nor vice versa). In the diseasome, NHLs shared more genetic etiology with ADs than solid cancers (p = .0041). A meta-analysis (combing AD with NHL) implicated genes of apoptosis and telomere length. This GWAS-based analysis four NHL subtypes and three ADs revealed few weakly-associated shared loci, explaining little total risk. This suggests common genetic variation, as assessed by GWAS in these sample sizes, may not be the primary explanation for the link between these ADs and NHLs.

Soluble urokinase plasminogen activator receptor (suPAR) levels predict damage accrual in patients with recent-onset systemic lupus erythematosus.

OBJECTIVE: The soluble urokinase plasminogen activator receptor (suPAR) has potential as a prognosis and severity biomarker in several inflammatory and infectious diseases. In a previous cross-sectional study, suPAR levels were shown to reflect damage accrual in cases of systemic lupus erythematosus (SLE). Herein, we evaluated suPAR as a predictor of future organ damage in recent-onset SLE. METHODS: Included were 344 patients from the Systemic Lupus International Collaborating Clinics (SLICC) Inception Cohort who met the 1997 American College of Rheumatology classification criteria with 5-years of follow-up data available. Baseline sera from patients and age- and sex-matched controls were assayed for suPAR. Organ damage was assessed annually using the SLICC/ACR damage index (SDI). RESULTS: The levels of suPAR were higher in patients who accrued damage, particularly those with SDI≥2 at 5 years (N = 32, 46.8% increase, p = 0.004), as compared to patients without damage. Logistic regression analysis revealed a significant impact of suPAR on SDI outcome (SDI≥2; OR = 1.14; 95% CI 1.03-1.26), also after adjustment for confounding factors. In an optimized logistic regression to predict damage, suPAR persisted as a predictor, together with baseline disease activity (SLEDAI-2K), age, and non-Caucasian ethnicity (model AUC = 0.77). Dissecting SDI into organ systems revealed higher suPAR levels in patients who developed musculoskeletal damage (SDI≥1; p = 0.007). CONCLUSION: Prognostic biomarkers identify patients who are at risk of acquiring early damage and therefore need careful observation and targeted treatment strategies. Overall, suPAR constitutes an interesting biomarker for patient stratification and for identifying SLE patients who are at risk of acquiring organ damage during the first 5 years of disease.

Genome-wide meta-analysis reveals shared new loci in systemic seropositive rheumatic diseases.

OBJECTIVE: Immune-mediated inflammatory diseases (IMIDs) are heterogeneous and complex conditions with overlapping clinical symptoms and elevated familial aggregation, which suggests the existence of a shared genetic component. In order to identify this genetic background in a systematic fashion, we performed the first cross-disease genome-wide meta-analysis in systemic seropositive rheumatic diseases, namely, systemic sclerosis, systemic lupus erythematosus, rheumatoid arthritis and idiopathic inflammatory myopathies. METHODS: We meta-analysed ~6.5 million single nucleotide polymorphisms in 11 678 cases and 19 704 non-affected controls of European descent populations. The functional roles of the associated variants were interrogated using publicly available databases. RESULTS: Our analysis revealed five shared genome-wide significant independent loci that had not been previously associated with these diseases: NAB1, KPNA4-ARL14, DGQK, LIMK1 and PRR12. All of these loci are related with immune processes such as interferon and epidermal growth factor signalling, response to methotrexate, cytoskeleton dynamics and coagulation cascade. Remarkably, several of the associated loci are known key players in autoimmunity, which supports the validity of our results. All the associated variants showed significant functional enrichment in DNase hypersensitivity sites, chromatin states and histone marks in relevant immune cells, including shared expression quantitative trait loci. Additionally, our results were significantly enriched in drugs that are being tested for the treatment of the diseases under study. CONCLUSIONS: We have identified shared new risk loci with functional value across diseases and pinpoint new potential candidate loci that could be further investigated. Our results highlight the potential of drug repositioning among related systemic seropositive rheumatic IMIDs.

Sequencing of the MHC region defines HLA-DQA1 as the major genetic risk for seropositive rheumatoid arthritis in Han Chinese population.

OBJECTIVE: The strong genetic contribution of the major histocompatibility complex (MHC) region to rheumatoid arthritis (RA) has been generally attributed to human leukocyte antigen (HLA)-DRB1. However, due to the high polymorphisms and linkage disequilibrium within MHC, it is difficult to define novel and/or independent genetic risks using conventional HLA genotyping or chip-based microarray technology. This study aimed to identify novel RA risk variants by performing deep sequencing for MHC. METHODS: We first conducted target sequencing for the entire MHC region in 357 anticitrullinated protein antibodies (ACPA)-positive patients with RA and 1001 healthy controls, and then performed HLA typing in an independent case-control cohort consisting of 1415 samples for validation. All study subjects were Han Chinese. Genetic associations for RA susceptibility and severity were analysed. Comparative modelling was constructed to predict potential functions for the newly discovered RA association variants. RESULTS: HLA-DQα1:160D conferred the strongest and independent susceptibility to ACPA-positive RA (p=6.16×10-36, OR=2.29). DRß1:37N had an independent protective effect (p=5.81×10-16, OR=0.49). As predicted by comparative modelling, the negatively charged DQα1:160D stabilises the dimer of dimers, thus may lead to an increased T cell activation. The negatively charged DRß1:37N encoding alleles preferentially bind with epitope P9 arginine, thus may result in a decreased RA susceptibility. CONCLUSIONS: We provide the first evidence that HLA-DQα1:160D, instead of HLA-DRB1*0405, is the strongest and independent genetic risk for ACPA-positive RA in Han Chinese. Our study also illustrates the value of deep sequencing for fine-mapping disease risk variants in the MHC region.

Association of response to TNF inhibitors in rheumatoid arthritis with quantitative trait loci for CD40 and CD39.

OBJECTIVES: We sought to investigate whether genetic effects on response to TNF inhibitors (TNFi) in rheumatoid arthritis (RA) could be localised by considering known genetic susceptibility loci for relevant traits and to evaluate the usefulness of these genetic loci for stratifying drug response. METHODS: We studied the relation of TNFi response, quantified by change in swollen joint counts ( Δ SJC) and erythrocyte sedimentation rate ( Δ ESR) with locus-specific scores constructed from genome-wide assocation study summary statistics in 2938 genotyped individuals: 37 scores for RA; scores for 19 immune cell traits; scores for expression or methylation of 93 genes with previously reported associations between transcript level and drug response. Multivariate associations were evaluated in penalised regression models by cross-validation. RESULTS: We detected a statistically significant association between Δ SJC and the RA score at the CD40 locus (p=0.0004) and an inverse association between Δ SJC and the score for expression of CD39 on CD4 T cells (p=0.00005). A previously reported association between CD39 expression on regulatory T cells and response to methotrexate was in the opposite direction. In stratified analysis by concomitant methotrexate treatment, the inverse association was stronger in the combination therapy group and dissipated in the TNFi monotherapy group. Overall, ability to predict TNFi response from genotypic scores was limited, with models explaining less than 1% of phenotypic variance. CONCLUSIONS: The association with the CD39 trait is difficult to interpret because patients with RA are often prescribed TNFi after failing to respond to methotrexate. The CD39 and CD40 pathways could be relevant for targeting drug therapy.

Focused HLA analysis in Caucasians with myositis identifies significant associations with autoantibody subgroups.

OBJECTIVES: Idiopathic inflammatory myopathies (IIM) are a spectrum of rare autoimmune diseases characterised clinically by muscle weakness and heterogeneous systemic organ involvement. The strongest genetic risk is within the major histocompatibility complex (MHC). Since autoantibody presence defines specific clinical subgroups of IIM, we aimed to correlate serotype and genotype, to identify novel risk variants in the MHC region that co-occur with IIM autoantibodies. METHODS: We collected available autoantibody data in our cohort of 2582 Caucasian patients with IIM. High resolution human leucocyte antigen (HLA) alleles and corresponding amino acid sequences were imputed using SNP2HLA from existing genotyping data and tested for association with 12 autoantibody subgroups. RESULTS: We report associations with eight autoantibodies reaching our study-wide significance level of p<2.9×10-5. Associations with the 8.1 ancestral haplotype were found with anti-Jo-1 (HLA-B*08:01, p=2.28×10-53 and HLA-DRB1*03:01, p=3.25×10-9), anti-PM/Scl (HLA-DQB1*02:01, p=1.47×10-26) and anti-cN1A autoantibodies (HLA-DRB1*03:01, p=1.40×10-11). Associations independent of this haplotype were found with anti-Mi-2 (HLA-DRB1*07:01, p=4.92×10-13) and anti-HMGCR autoantibodies (HLA-DRB1*11, p=5.09×10-6). Amino acid positions may be more strongly associated than classical HLA associations; for example with anti-Jo-1 autoantibodies and position 74 of HLA-DRB1 (p=3.47×10-64) and position 9 of HLA-B (p=7.03×10-11). We report novel genetic associations with HLA-DQB1 anti-TIF1 autoantibodies and identify haplotypes that may differ between adult-onset and juvenile-onset patients with these autoantibodies. CONCLUSIONS: These findings provide new insights regarding the functional consequences of genetic polymorphisms within the MHC. As autoantibodies in IIM correlate with specific clinical features of disease, understanding genetic risk underlying development of autoantibody profiles has implications for future research.