J. Bart A. Crusius

Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci

To identify new genetic risk factors for rheumatoid arthritis, we conducted a genome-wide association study meta-analysis of 5,539 autoantibody-positive individuals with rheumatoid arthritis (cases) and 20,169 controls of European descent, followed by replication in an independent set of 6,768 rheumatoid arthritis cases and 8,806 controls. Of 34 SNPs selected for replication, 7 new rheumatoid arthritis risk alleles were identified at genome-wide significance (P < 5 × 10−8) in an analysis of all 41,282 samples. The associated SNPs are near genes of known immune function, including IL6ST, SPRED2, RBPJ, CCR6, IRF5 and PXK. We also refined associations at two established rheumatoid arthritis risk loci (IL2RA and CCL21) and confirmed the association at AFF3. These new associations bring the total number of confirmed rheumatoid arthritis risk loci to 31 among individuals of European ancestry. An additional 11 SNPs replicated at P < 0.05, many of which are validated autoimmune risk alleles, suggesting that most represent genuine rheumatoid arthritis risk alleles.

Common variants at CD40 and other loci confer risk of rheumatoid arthritis

To identify rheumatoid arthritis risk loci in European populations, we conducted a meta-analysis of two published genome-wide association (GWA) studies totaling 3,393 cases and 12,462 controls. We genotyped 31 top-ranked SNPs not previously associated with rheumatoid arthritis in an independent replication of 3,929 autoantibody-positive rheumatoid arthritis cases and 5,807 matched controls from eight separate collections. We identified a common variant at the CD40 gene locus (rs4810485, P = 0.0032 replication, P = 8.2 × 10−9 overall, OR = 0.87). Along with other associations near TRAF1 (refs. 2,3) and TNFAIP3 (refs. 4,5), this implies a central role for the CD40 signaling pathway in rheumatoid arthritis pathogenesis. We also identified association at the CCL21 gene locus (rs2812378, P = 0.00097 replication, P = 2.8 × 10−7 overall), a gene involved in lymphocyte trafficking. Finally, we identified evidence of association at four additional gene loci: MMEL1-TNFRSF14 (rs3890745, P = 0.0035 replication, P = 1.1 × 10−7 overall), CDK6 (rs42041, P = 0.010 replication, P = 4.0 × 10−6 overall), PRKCQ (rs4750316, P = 0.0078 replication, P = 4.4 × 10−6 overall), and KIF5A-PIP4K2C (rs1678542, P = 0.0026 replication, P = 8.8 × 10−8 overall).

Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci

Ankylosing spondylitis is a common, highly heritable inflammatory arthritis affecting primarily the spine and pelvis. In addition to HLA-B*27 alleles, 12 loci have previously been identified that are associated with ankylosing spondylitis in populations of European ancestry, and 2 associated loci have been identified in Asians. In this study, we used the Illumina Immunochip microarray to perform a case-control association study involving 10,619 individuals with ankylosing spondylitis (cases) and 15,145 controls. We identified 13 new risk loci and 12 additional ankylosing spondylitis–associated haplotypes at 11 loci. Two ankylosing spondylitis–associated regions have now been identified encoding four aminopeptidases that are involved in peptide processing before major histocompatibility complex (MHC) class I presentation. Protective variants at two of these loci are associated both with reduced aminopeptidase function and with MHC class I cell surface expression.

Genetic variants at CD28, PRDM1, and CD2/CD58 are associated with rheumatoid arthritis risk

To discover new rheumatoid arthritis (RA) risk loci, we systematically examined 370 SNPs from 179 independent loci with P < 0.001 in a published meta-analysis of RA genome-wide association studies (GWAS) of 3,393 cases and 12,462 controls. We used Gene Relationships Across Implicated Loci (GRAIL), a computational method that applies statistical text mining to PubMed abstracts, to score these 179 loci for functional relationships to genes in 16 established RA disease loci. We identified 22 loci with a significant degree of functional connectivity. We genotyped 22 representative SNPs in an independent set of 7,957 cases and 11,958 matched controls. Three were convincingly validated: CD2-CD58 (rs11586238, P = 1 × 10−6 replication, P = 1 × 10−9 overall), CD28 (rs1980422, P = 5 × 10−6 replication, P = 1 × 10−9 overall) and PRDM1 (rs548234, P = 1 × 10−5 replication, P = 2 × 10−8 overall). An additional four were replicated (P < 0.0023): TAGAP (rs394581, P = 0.0002 replication, P = 4 × 10−7 overall), PTPRC (rs10919563, P = 0.0003 replication, P = 7 × 10−7 overall), TRAF6-RAG1 (rs540386, P = 0.0008 replication, P = 4 × 10−6 overall) and FCGR2A (rs12746613, P = 0.0022 replication, P = 2 × 10−5 overall). Many of these loci are also associated to other immunologic diseases.

TNF-α promoter polymorphisms, production and susceptibility to multiple sclerosis in different groups of patients

TNF-alpha production in whole blood cultures upon stimulation with LPS was determined in 179 individuals from 61 families in order to characterise the magnitude of inherited differences in TNF-alpha production. The three families characterised by highest TNF production showed 7.1 +/- 0.3 ng TNF/ml upon culture with 10 ng LPS and 10.2 +/- 0.2 ng TNF/ml upon culture with 1000 ng LPS. in contrast to the three families characterised by the lowest TNF production that showed a production of 1.6 +/- 0.1 ng TNF upon culture with 10 ng and 2.5 +/- 0.2 ng/ml upon culture with 1000 ng LPS/ml. This difference could not be attributed to the promoter polymorphisms -308 G to A. -238 G to A or -376 G to A, although the -238 GA donors produced 2.1 +/- 0.9 ng TNF upon culture with 10 ng endotoxin compared to 3.2 +/- 2.2 ng TNF for the -238 GG donors. In line with these results the frequency of the -238 GG genotype was increased in hospitalized MS patients in a nursing home (100% 238GG, n = 57) compared to MS patients in an outpatients clinic (94% 238GG, n = 98) or Dutch controls (90% 238GG, n = 180). These results suggest that the -238 GG genotype is differently distributed in hospitalized MS patients in a nursing home.

Genetic analysis of innate immunity in Crohn's disease and ulcerative colitis identifies two susceptibility loci harboring CARD9 and IL18RAP.

The two main phenotypes of inflammatory bowel disease (IBD)--Crohns disease (CD) and ulcerative colitis (UC)--are chronic intestinal inflammatory disorders with a complex genetic background. Using a three-stage design, we performed a functional candidate-gene analysis of innate immune pathway in IBD. In phase I, we typed 354 SNPs from 85 innate immunity genes in 520 Dutch IBD patients (284 CD, 236 UC) and 808 controls. In phase II, ten autosomal SNPs showing association at p < 0.006 in phase I were replicated in a second cohort of 545 IBD patients (326 CD, 219 UC) and 360 controls. In phase III, four SNPs with p < 0.01 in the combined phase I and phase II analysis were genotyped in an additional 786 IBD samples (452 CD, 334 UC) and 768 independent controls. Joint analysis of 1851 IBD patients (1062 CD, 789 UC) and 1936 controls demonstrated strong association to the IL18RAP rs917997 SNP for both CD and UC (p(IBD) 1.9 x 10(-8); OR 1.35). Association in CD is independently supported by the Crohns disease dataset of the Wellcome Trust Case Control Consortium (imputed SNP rs917997, p = 9.19 x 10(-4)). In addition, an association of the CARD9 rs10870077 SNP to CD and UC was observed (p(IBD) = 3.25 x 10(-5); OR 1.21). Both genes are located in extended haplotype blocks on 2q11-2q12 and 9q34.3, respectively. Our results indicate two IBD loci and further support the importance of the innate immune system in the predisposition to both CD and UC.

Effective Detection of Human Leukocyte Antigen Risk Alleles in Celiac Disease Using Tag Single Nucleotide Polymorphisms

Background The HLA genes, located in the MHC region on chromosome 6p21.3, play an important role in many autoimmune disorders, such as celiac disease (CD), type 1 diabetes (T1D), rheumatoid arthritis, multiple sclerosis, psoriasis and others. Known HLA variants that confer risk to CD, for example, include DQA1*05/DQB1*02 (DQ2.5) and DQA1*03/DQB1*0302 (DQ8). To diagnose the majority of CD patients and to study disease susceptibility and progression, typing these strongly associated HLA risk factors is of utmost importance. However, current genotyping methods for HLA risk factors involve many reactions, and are complicated and expensive. We sought a simple experimental approach using tagging SNPs that predict the CD-associated HLA risk factors. Methodology Our tagging approach exploits linkage disequilibrium between single nucleotide polymorphism (SNPs) and the CD-associated HLA risk factors DQ2.5 and DQ8 that indicate direct risk, and DQA1*0201/DQB1*0202 (DQ2.2) and DQA1*0505/DQB1*0301 (DQ7) that attribute to the risk of DQ2.5 to CD. To evaluate the predictive power of this approach, we performed an empirical comparison of the predicted DQ types, based on these six tag SNPs, with those executed with current validated laboratory typing methods of the HLA-DQA1 and -DQB1 genes in three large cohorts. The results were validated in three European celiac populations. Conclusion Using this method, only six SNPs were needed to predict the risk types carried by >95% of CD patients. We determined that for this tagging approach the sensitivity was >0.991, specificity >0.996 and the predictive value >0.948. Our results show that this tag SNP method is very accurate and provides an excellent basis for population screening for CD. This method is broadly applicable in European populations.

Rheumatoid Arthritis Risk Allele PTPRC Is Also Associated With Response to Anti-Tumor Necrosis Factor alpha Therapy

OBJECTIVE Anti-tumor necrosis factor alpha (anti-TNF) therapy is a mainstay of treatment in rheumatoid arthritis (RA). The aim of the present study was to test established RA genetic risk factors to determine whether the same alleles also influence the response to anti-TNF therapy. METHODS A total of 1,283 RA patients receiving etanercept, infliximab, or adalimumab therapy were studied from among an international collaborative consortium of 9 different RA cohorts. The primary end point compared RA patients with a good treatment response according to the European League Against Rheumatism (EULAR) response criteria (n = 505) with RA patients considered to be nonresponders (n = 316). The secondary end point was the change from baseline in the level of disease activity according to the Disease Activity Score in 28 joints (triangle upDAS28). Clinical factors such as age, sex, and concomitant medications were tested as possible correlates of treatment response. Thirty-one single-nucleotide polymorphisms (SNPs) associated with the risk of RA were genotyped and tested for any association with treatment response, using univariate and multivariate logistic regression models. RESULTS Of the 31 RA-associated risk alleles, a SNP at the PTPRC (also known as CD45) gene locus (rs10919563) was associated with the primary end point, a EULAR good response versus no response (odds ratio [OR] 0.55, P = 0.0001 in the multivariate model). Similar results were obtained using the secondary end point, the triangle upDAS28 (P = 0.0002). There was suggestive evidence of a stronger association in autoantibody-positive patients with RA (OR 0.55, 95% confidence interval [95% CI] 0.39-0.76) as compared with autoantibody-negative patients (OR 0.90, 95% CI 0.41-1.99). CONCLUSION Statistically significant associations were observed between the response to anti-TNF therapy and an RA risk allele at the PTPRC gene locus. Additional studies will be required to replicate this finding in additional patient collections.

Genome-wide association study and gene expression analysis identifies CD84 as a predictor of response to etanercept therapy in rheumatoid arthritis.

Anti-tumor necrosis factor alpha (anti-TNF) biologic therapy is a widely used treatment for rheumatoid arthritis (RA). It is unknown why some RA patients fail to respond adequately to anti-TNF therapy, which limits the development of clinical biomarkers to predict response or new drugs to target refractory cases. To understand the biological basis of response to anti-TNF therapy, we conducted a genome-wide association study (GWAS) meta-analysis of more than 2 million common variants in 2,706 RA patients from 13 different collections. Patients were treated with one of three anti-TNF medications: etanercept (n = 733), infliximab (n = 894), or adalimumab (n = 1,071). We identified a SNP (rs6427528) at the 1q23 locus that was associated with change in disease activity score (ΔDAS) in the etanercept subset of patients (P = 8×10−8), but not in the infliximab or adalimumab subsets (P>0.05). The SNP is predicted to disrupt transcription factor binding site motifs in the 3′ UTR of an immune-related gene, CD84, and the allele associated with better response to etanercept was associated with higher CD84 gene expression in peripheral blood mononuclear cells (P = 1×10−11 in 228 non-RA patients and P = 0.004 in 132 RA patients). Consistent with the genetic findings, higher CD84 gene expression correlated with lower cross-sectional DAS (P = 0.02, n = 210) and showed a non-significant trend for better ΔDAS in a subset of RA patients with gene expression data (n = 31, etanercept-treated). A small, multi-ethnic replication showed a non-significant trend towards an association among etanercept-treated RA patients of Portuguese ancestry (n = 139, P = 0.4), but no association among patients of Japanese ancestry (n = 151, P = 0.8). Our study demonstrates that an allele associated with response to etanercept therapy is also associated with CD84 gene expression, and further that CD84 expression correlates with disease activity. These findings support a model in which CD84 genotypes and/or expression may serve as a useful biomarker for response to etanercept treatment in RA patients of European ancestry.

Analysis of SNPs with an effect on gene expression identifies UBE2L3 and BCL3 as potential new risk genes for Crohn's disease

Genome-wide association studies (GWAS) for Crohns disease (CD) have identified loci explaining approximately 20% of the total genetic risk of CD. Part of the other genetic risk loci is probably partly hidden among signals discarded by the multiple testing correction needed in the analysis of GWAS data. Strategies for finding these hidden loci require large replication cohorts and are costly to perform. We adopted a strategy of selecting SNPs for follow-up that showed a correlation to gene expression [cis-expression quantitative trait loci (eQTLs)] since these have been shown more likely to be trait-associated. First we show that there is an overrepresentation of cis-eQTLs in the known CD-associated loci. Then SNPs were selected for follow-up by screening the top 500 SNP hits from a CD GWAS data set. We identified 10 cis-eQTL SNPs. These 10 SNPs were tested for association with CD in two independent cohorts of Dutch CD patients (1539) and healthy controls (2648). In a combined analysis, we identified two cis-eQTL SNPs that were associated with CD rs2298428 in UBE2L3 (P=5.22x10(-5)) and rs2927488 in BCL3 (P=2.94x10(-4)). After adding additional publicly available data from a previously reported meta-analysis, the association with rs2298428 almost reached genome-wide significance (P=2.40x10(-7)) and the association with rs2927488 was corroborated (P=6.46x10(-4)). We have identified UBE2L3 and BCL3 as likely novel risk genes for CD. UBE2L3 is also associated with other immune-mediated diseases. These results show that eQTL-based pre-selection for follow-up is a useful approach for identifying risk loci from a moderately sized GWAS.