Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci - PubMed (original) (raw)

Meta-Analysis

doi: 10.1038/ng.582. Epub 2010 May 9.

Soumya Raychaudhuri, Elaine F Remmers, Gang Xie, Stephen Eyre, Brian P Thomson, Yonghong Li, Fina A S Kurreeman, Alexandra Zhernakova, Anne Hinks, Candace Guiducci, Robert Chen, Lars Alfredsson, Christopher I Amos, Kristin G Ardlie; BIRAC Consortium; Anne Barton, John Bowes, Elisabeth Brouwer, Noel P Burtt, Joseph J Catanese, Jonathan Coblyn, Marieke J H Coenen, Karen H Costenbader, Lindsey A Criswell, J Bart A Crusius, Jing Cui, Paul I W de Bakker, Philip L De Jager, Bo Ding, Paul Emery, Edward Flynn, Pille Harrison, Lynne J Hocking, Tom W J Huizinga, Daniel L Kastner, Xiayi Ke, Annette T Lee, Xiangdong Liu, Paul Martin, Ann W Morgan, Leonid Padyukov, Marcel D Posthumus, Timothy R D J Radstake, David M Reid, Mark Seielstad, Michael F Seldin, Nancy A Shadick, Sophia Steer, Paul P Tak, Wendy Thomson, Annette H M van der Helm-van Mil, Irene E van der Horst-Bruinsma, C Ellen van der Schoot, Piet L C M van Riel, Michael E Weinblatt, Anthony G Wilson, Gert Jan Wolbink, B Paul Wordsworth; YEAR Consortium; Cisca Wijmenga, Elizabeth W Karlson, Rene E M Toes, Niek de Vries, Ann B Begovich, Jane Worthington, Katherine A Siminovitch, Peter K Gregersen, Lars Klareskog, Robert M Plenge

Affiliations

Meta-Analysis

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

Eli A Stahl et al. Nat Genet. 2010 Jun.

Abstract

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 x 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.

PubMed Disclaimer

Conflict of interest statement

The authors declare that we have no competing financial interest.

Figures

Figure 1

Figure 1. Association with RA risk across 4 loci

Regional association plots show strength of association (-Log(_P_-value)) versus chromosomal position (kilobases, kb) for all SNPs across 1 Megabase (Mb) regions centered on the newly validated SNPs (labeled). PGWAS values are plotted with red/white diamonds for all SNPs, shaded by the degree of LD (r2, see legend) with the validated SNP (larger red diamond). Poverall in combined analysis of GWAS and replication collections is plotted with the blue diamond. Local recombination rates estimated from HapMap CEU (centi-Morgans per Mb, blue line) are plotted against the secondary y-axis, showing recombination hotspots across the region. Labeled green arrows below the plots indicate genes and their orientations. (a) 2p14, SPRED2 locus. (b) 5q11, _IL6ST_-ANKRD55 locus. (c) 5q21, C5orf13 locus. (d) 10p15, IL2RA locus.

Figure 1

Figure 1. Association with RA risk across 4 loci

Regional association plots show strength of association (-Log(_P_-value)) versus chromosomal position (kilobases, kb) for all SNPs across 1 Megabase (Mb) regions centered on the newly validated SNPs (labeled). PGWAS values are plotted with red/white diamonds for all SNPs, shaded by the degree of LD (r2, see legend) with the validated SNP (larger red diamond). Poverall in combined analysis of GWAS and replication collections is plotted with the blue diamond. Local recombination rates estimated from HapMap CEU (centi-Morgans per Mb, blue line) are plotted against the secondary y-axis, showing recombination hotspots across the region. Labeled green arrows below the plots indicate genes and their orientations. (a) 2p14, SPRED2 locus. (b) 5q11, _IL6ST_-ANKRD55 locus. (c) 5q21, C5orf13 locus. (d) 10p15, IL2RA locus.

Figure 1

Figure 1. Association with RA risk across 4 loci

Regional association plots show strength of association (-Log(_P_-value)) versus chromosomal position (kilobases, kb) for all SNPs across 1 Megabase (Mb) regions centered on the newly validated SNPs (labeled). PGWAS values are plotted with red/white diamonds for all SNPs, shaded by the degree of LD (r2, see legend) with the validated SNP (larger red diamond). Poverall in combined analysis of GWAS and replication collections is plotted with the blue diamond. Local recombination rates estimated from HapMap CEU (centi-Morgans per Mb, blue line) are plotted against the secondary y-axis, showing recombination hotspots across the region. Labeled green arrows below the plots indicate genes and their orientations. (a) 2p14, SPRED2 locus. (b) 5q11, _IL6ST_-ANKRD55 locus. (c) 5q21, C5orf13 locus. (d) 10p15, IL2RA locus.

Figure 1

Figure 1. Association with RA risk across 4 loci

Regional association plots show strength of association (-Log(_P_-value)) versus chromosomal position (kilobases, kb) for all SNPs across 1 Megabase (Mb) regions centered on the newly validated SNPs (labeled). PGWAS values are plotted with red/white diamonds for all SNPs, shaded by the degree of LD (r2, see legend) with the validated SNP (larger red diamond). Poverall in combined analysis of GWAS and replication collections is plotted with the blue diamond. Local recombination rates estimated from HapMap CEU (centi-Morgans per Mb, blue line) are plotted against the secondary y-axis, showing recombination hotspots across the region. Labeled green arrows below the plots indicate genes and their orientations. (a) 2p14, SPRED2 locus. (b) 5q11, _IL6ST_-ANKRD55 locus. (c) 5q21, C5orf13 locus. (d) 10p15, IL2RA locus.

Figure 2

Figure 2. Previously validated autoimmune SNPs tested in our replication study

Eighteen SNPs tested in our replication samples were in LD (defined as r2 > 0.3) with a validated autoimmune risk allele. Of these, 5 were validated as RA risk alleles in our study (Poverall < 5×10−8, inner most circle), 6 were suggestive associations (_Preplication_ <0.05 but _Poverall_ > 5×10−8), and 6 demonstrated no evidence of association in our replication samples (Preplication ≥0.05). For the 12 SNPs with suggestive or no evidence of association, each SNP is plotted by the strength of association with RA risk in the replication samples; those closer to the inner circle have more significant Preplication. All of the RA risk alleles confer risk in the same direction as the validated autoimmune risk alleles (when the same allele or a near perfect proxy was tested). We include the following as “autoimmune” diseases in our study, listed on the outside of the circle, although these reflect diseases along the autoimmuneinflammatory spectrum: systemic lupus erythematosus (SLE), celiac disease, Crohn’s disease, multiple sclerosis (MS), psoriasis, and type 1 diabetes (T1D); other autoimmune diseases are not included (e.g., autoimmune thyroiditis). The haplotype tagged by the IL2RA SNP, rs706778, is associated with T1D and MS,; the SH2B3 SNP is associated with both T1D and Celiac disease,. The CCR6 SNP rs3093023 is in partial LD with a SNP associated with Crohn’s disease (rs2301436, r2 = 0.48). The AFF3 SNP has an equivocal association with T1D (where the associated SNP is rs9653442). The IL2-IL21 SNP tested in our study, rs13119723, is in LD (r2=0.67) with a SNP previously implicated in both Celiac disease and RA (rs6822844),, but only in partial LD (r2=0.09) with a T1D SNP (rs4505848). The CD19-NFATC2IP SNP tested in our study, rs8045689, was selected because of GRAIL; it is in partial LD (r2=0.38) with a SNP associated with T1D (rs4788084). The TNIP1 SNP in our study, rs6889239, is in strong LD with an SLE SNP (rs7708392, r2=0.91), but not in LD with another TNIP1 SNP associated with psoriasis (rs17728338, r2<0.01). The ZEB1 SNP (rs2793108) was from the May 2009 release of T1D base, although this SNP did not appear in a subsequent publication. The PXK SLE SNP was tested in this study but not shown, as it is in weak LD with the RA risk SNP (r2 = 0.15 between rs6445975 and rs13315591); the RA SNP was selected because of _P_GWAS<10−6, not because of its association with another autoimmune disease. Note that there are SNPs associated with RA and other autoimmune diseases not shown; we only include those SNPs tested as part of the current study.

Similar articles

Cited by

References

    1. Silman AJ, Pearson JE. Epidemiology and genetics of rheumatoid arthritis. Arthritis Res. 2002;4(Suppl 3):S265–72. - PMC - PubMed
    1. Klareskog L, Catrina AI, Paget S. Rheumatoid arthritis. Lancet. 2009;373:659–72. - PubMed
    1. Suzuki A, et al. Functional haplotypes of PADI4, encoding citrullinating enzyme peptidylarginine deiminase 4, are associated with rheumatoid arthritis. Nat Genet. 2003;34:395–402. - PubMed
    1. Begovich AB, et al. A missense single-nucleotide polymorphism in a gene encoding a protein tyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis. Am J Hum Genet. 2004;75:330–7. - PMC - PubMed
    1. Kurreeman FA, et al. A candidate gene approach identifies the TRAF1/C5 region as a risk factor for rheumatoid arthritis. PLoS Med. 2007;4:e278. - PMC - PubMed

Publication types

MeSH terms

Substances

Grants and funding

LinkOut - more resources