Meta-analysis identifies nine new loci associated with rheumatoid arthritis in the Japanese population

Rheumatoid arthritis is a common autoimmune disease characterized by chronic inflammation. We report a meta-analysis of genome-wide association studies (GWAS) in a Japanese population including 4,074 individuals with rheumatoid arthritis (cases) and 16,891 controls, followed by a replication in 5,277 rheumatoid arthritis cases and 21,684 controls. Our study identified nine loci newly associated with rheumatoid arthritis at a threshold of P < 5.0 × 10(-8), including B3GNT2, ANXA3, CSF2, CD83, NFKBIE, ARID5B, PDE2A-ARAP1, PLD4 and PTPN2. ANXA3 was also associated with susceptibility to systemic lupus erythematosus (P = 0.0040), and B3GNT2 and ARID5B were associated with Graves' disease (P = 3.5 × 10(-4) and 2.9 × 10(-4), respectively). We conducted a multi-ancestry comparative analysis with a previous meta-analysis in individuals of European descent (5,539 rheumatoid arthritis cases and 20,169 controls). This provided evidence of shared genetic risks of rheumatoid arthritis between the populations.     

A functional variant in FCRL3, encoding Fc receptor-like 3, is associated with rheumatoid arthritis and several autoimmunities

Rheumatoid arthritis is a common autoimmune disease with a complex genetic etiology. Here we identify a SNP in the promoter region of FCRL3, a member of the Fc receptor-like family, that is associated with susceptibility to rheumatoid arthritis (odds ratio = 2.15, P = 0.00000085). This polymorphism alters the binding affinity of nuclear factor-kappaB and regulates FCRL3 expression. We observed high FCRL3 expression on B cells and augmented autoantibody production in individuals with the disease-susceptible genotype. We also found associations between the SNP and susceptibility to autoimmune thyroid disease and systemic lupus erythematosus. FCRL3 may therefore have a pivotal role in autoimmunity.     

A nonsynonymous functional variant in integrin-alpha(M) (encoded by ITGAM) is associated with systemic lupus erythematosus

We identified and replicated an association between ITGAM (CD11b) at 16p11.2 and risk of systemic lupus erythematosus (SLE) in 3,818 individuals of European descent. The strongest association was at a nonsynonymous SNP, rs1143679 (P = 1.7 x 10(-17), odds ratio = 1.78). We further replicated this association in two independent samples of individuals of African descent (P = 0.0002 and 0.003; overall meta-analysis P = 6.9 x 10(-22)). The genetic association between ITGAM and SLE implicates the alpha(M)beta2-integrin adhesion pathway in disease development.     

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 x 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 x 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 x 10(-7) overall), CDK6 (rs42041, P = 0.010 replication, P = 4.0 x 10(-6) overall), PRKCQ (rs4750316, P = 0.0078 replication, P = 4.4 x 10(-6) overall), and KIF5A-PIP4K2C (rs1678542, P = 0.0026 replication, P = 8.8 x 10(-8) overall).     

A common haplotype of interferon regulatory factor 5 (IRF5) regulates splicing and expression and is associated with increased risk of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a complex autoimmune disease characterized by activation of the type I interferon (IFN) pathway. Here we convincingly replicate association of the IFN regulatory factor 5 (IRF5) rs2004640 T allele with SLE in four independent case-control cohorts (P = 4.4 x 10(-16)) and by family-based transmission disequilibrium test analysis (P = 0.0006). The rs2004640 T allele creates a 5' donor splice site in an alternate exon 1 of IRF5, allowing expression of several unique IRF5 isoforms. We also identify an independent cis-acting variant associated with elevated expression of IRF5 and linked to the exon 1B splice site. Haplotypes carrying the variant associated with elevated expression and lacking the exon 1B donor site do not confer risk of SLE. Thus, a common IRF5 haplotype driving elevated expression of multiple unique isoforms of IRF5 is an important genetic risk factor for SLE, establishing a causal role for type I IFN pathway genes in human autoimmunity.     

Two independent alleles at 6q23 associated with risk of rheumatoid arthritis

To identify susceptibility alleles associated with rheumatoid arthritis, we genotyped 397 individuals with rheumatoid arthritis for 116,204 SNPs and carried out an association analysis in comparison to publicly available genotype data for 1,211 related individuals from the Framingham Heart Study. After evaluating and adjusting for technical and population biases, we identified a SNP at 6q23 (rs10499194, approximately 150 kb from TNFAIP3 and OLIG3) that was reproducibly associated with rheumatoid arthritis both in the genome-wide association (GWA) scan and in 5,541 additional case-control samples (P = 10(-3), GWA scan; P < 10(-6), replication; P = 10(-9), combined). In a concurrent study, the Wellcome Trust Case Control Consortium (WTCCC) has reported strong association of rheumatoid arthritis susceptibility to a different SNP located 3.8 kb from rs10499194 (rs6920220; P = 5 x 10(-6) in WTCCC). We show that these two SNP associations are statistically independent, are each reproducible in the comparison of our data and WTCCC data, and define risk and protective haplotypes for rheumatoid arthritis at 6q23.     

Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM, PXK, KIAA1542 and other loci

Systemic lupus erythematosus (SLE) is a common systemic autoimmune disease with complex etiology but strong clustering in families (lambda(S) = approximately 30). We performed a genome-wide association scan using 317,501 SNPs in 720 women of European ancestry with SLE and in 2,337 controls, and we genotyped consistently associated SNPs in two additional independent sample sets totaling 1,846 affected women and 1,825 controls. Aside from the expected strong association between SLE and the HLA region on chromosome 6p21 and the previously confirmed non-HLA locus IRF5 on chromosome 7q32, we found evidence of association with replication (1.1 x 10(-7) < P(overall) < 1.6 x 10(-23); odds ratio = 0.82-1.62) in four regions: 16p11.2 (ITGAM), 11p15.5 (KIAA1542), 3p14.3 (PXK) and 1q25.1 (rs10798269). We also found evidence for association (P < 1 x 10(-5)) at FCGR2A, PTPN22 and STAT4, regions previously associated with SLE and other autoimmune diseases, as well as at > or =9 other loci (P < 2 x 10(-7)). Our results show that numerous genes, some with known immune-related functions, predispose to SLE.     

Rheumatoid arthritis susceptibility loci at chromosomes 10p15, 12q13 and 22q13

The WTCCC study identified 49 SNPs putatively associated with rheumatoid arthritis at P = 1 x 10(-4) - 1 x 10(-5) (tier 3). Here we show that three of these SNPs, mapping to chromosome 10p15 (rs4750316), 12q13 (rs1678542) and 22q13 (rs3218253), are also associated (trend P = 4 x 10(-5), P = 4 x 10(-4) and P = 4 x 10(-4), respectively) in a validation study of 4,106 individuals with rheumatoid arthritis and an expanded reference group of 11,238 subjects, confirming them as true susceptibility loci in individuals of European ancestry.     

FCGR3B copy number variation is associated with susceptibility to systemic, but not organ-specific, autoimmunity

Naturally occurring variation in gene copy number is increasingly recognized as a heritable source of susceptibility to genetically complex diseases. Here we report strong association between FCGR3B copy number and risk of systemic lupus erythematosus (P = 2.7 x 10(-8)), microscopic polyangiitis (P = 2.9 x 10(-4)) and Wegener's granulomatosis in two independent cohorts from the UK (P = 3 x 10(-3)) and France (P = 1.1 x 10(-4)). We did not observe this association in the organ-specific Graves' disease or Addison's disease. Our findings suggest that low FCGR3B copy number, and in particular complete FCGR3B deficiency, has a key role in the development of systemic autoimmunity.     

Common variants at TRAF3IP2 are associated with susceptibility to psoriatic arthritis and psoriasis

Psoriatic arthritis (PsA) is an inflammatory joint disease that is distinct from other chronic arthritides and which is frequently accompanied by psoriasis vulgaris (PsV) and seronegativity for rheumatoid factor. We conducted a genome-wide association study in 609 German individuals with PsA (cases) and 990 controls with replication in 6 European cohorts including a total of 5,488 individuals. We replicated PsA associations at HLA-C and IL12B and identified a new association at TRAF3IP2 (rs13190932, P = 8.56 × 10?1?). TRAF3IP2 was also associated with PsV in a German cohort including 2,040 individuals (rs13190932, P = 1.95 × 10?3). Sequencing of the exons of TRAF3IP2 identified a coding variant (p.Asp10Asn, rs33980500) as the most significantly associated SNP (P = 1.13 × 10?2?, odds ratio = 1.95). Functional assays showed reduced binding of this TRAF3IP2 variant to TRAF6, suggesting altered modulation of immunoregulatory signals through altered TRAF interactions as a new and shared pathway for PsA and PsV.     

A genome-wide association study for celiac disease identifies risk variants in the region harboring IL2 and IL21

We tested 310,605 SNPs for association in 778 individuals with celiac disease and 1,422 controls. Outside the HLA region, the most significant finding (rs13119723; P = 2.0 x 10(-7)) was in the KIAA1109-TENR-IL2-IL21 linkage disequilibrium block. We independently confirmed association in two further collections (strongest association at rs6822844, 24 kb 5' of IL21; meta-analysis P = 1.3 x 10(-14), odds ratio = 0.63), suggesting that genetic variation in this region predisposes to celiac disease.     

Male-pattern baldness susceptibility locus at 20p11

We conducted a genome-wide association study for androgenic alopecia in 1,125 men and identified a newly associated locus at chromosome 20p11.22, confirmed in three independent cohorts (n = 1,650; OR = 1.60, P = 1.1 x 10(-14) for rs1160312). The one man in seven who harbors risk alleles at both 20p11.22 and AR (encoding the androgen receptor) has a sevenfold-increased odds of androgenic alopecia (OR = 7.12, P = 3.7 x 10(-15)).     

Loci on 20q13 and 21q22 are associated with pediatric-onset inflammatory bowel disease

Inflammatory bowel disease (IBD) is a common inflammatory disorder with complex etiology that involves both genetic and environmental triggers, including but not limited to defects in bacterial clearance, defective mucosal barrier and persistent dysregulation of the immune response to commensal intestinal bacteria. IBD is characterized by two distinct phenotypes: Crohn's disease (CD) and ulcerative colitis (UC). Previously reported GWA studies have identified genetic variation accounting for a small portion of the overall genetic susceptibility to CD and an even smaller contribution to UC pathogenesis. We hypothesized that stratification of IBD by age of onset might identify additional genes associated with IBD. To that end, we carried out a GWA analysis in a cohort of 1,011 individuals with pediatric-onset IBD and 4,250 matched controls. We identified and replicated significantly associated, previously unreported loci on chromosomes 20q13 (rs2315008[T] and rs4809330[A]; P = 6.30 x 10(-8) and 6.95 x 10(-8), respectively; odds ratio (OR) = 0.74 for both) and 21q22 (rs2836878[A]; P = 6.01 x 10(-8); OR = 0.73), located close to the TNFRSF6B and PSMG1 genes, respectively.     

Common variants in DVWA on chromosome 3p24.3 are associated with susceptibility to knee osteoarthritis

Susceptibility to osteoarthritis, the most common human arthritis, is known to be influenced by genetic factors. Through a genome-wide association study using approximately 100,000 SNPs, we have identified a previously unknown gene on chromosome 3p24.3, DVWA, which is associated with susceptibility to knee osteoarthritis. Expressed specifically in cartilage, DVWA encodes a 276-amino-acid protein with two regions corresponding to the von Willebrand factor type A domain (VWA domain). Several DVWA SNPs are significantly associated with knee osteoarthritis in two independent Japanese case-control cohorts. This association was replicated in a Japanese population cohort and a Han Chinese case-control cohort (combined P = 7.3 x 10(-11)). DVWA protein binds to beta-tubulin, and the binding is influenced by two highly associated missense SNPs (rs11718863 and rs7639618) located in the VWA domain. The Tyr169-Cys260 isoform of DVWA, which is overrepresented in knee osteoarthritis, showed weaker interaction. Our findings reveal a new paradigm for study of osteoarthritis etiology and pathogenesis.     

Autoimmune-associated lymphoid tyrosine phosphatase is a gain-of-function variant

A SNP in the gene PTPN22 is associated with type 1 diabetes, rheumatoid arthritis, lupus, Graves thyroiditis, Addison disease and other autoimmune disorders. T cells from carriers of the predisposing allele produce less interleukin-2 upon TCR stimulation, and the encoded phosphatase has higher catalytic activity and is a more potent negative regulator of T lymphocyte activation. We conclude that the autoimmune-predisposing allele is a gain-of-function mutant.     

Loss-of-function variant in DNASE1L3 causes a familial form of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a complex autoimmune disease that causes substantial morbidity. As is typical for many other multifactorial disorders, much of the heritability of SLE remains unknown. We identified a rare autosomal recessive form of SLE, in which autozygome analysis revealed a null mutation in the DNASE1L3 gene. The DNASE1L3-related SLE we describe was always pediatric in onset and correlated with a high frequency of lupus nephritis. Our findings confirm the critical role of impaired clearance of degraded DNA in SLE pathogenesis.     

Mutations in the gene encoding the 3'-5' DNA exonuclease TREX1 are associated with systemic lupus erythematosus

TREX1 acts in concert with the SET complex in granzyme A-mediated apoptosis, and mutations in TREX1 cause Aicardi-Goutières syndrome and familial chilblain lupus. Here, we report monoallelic frameshift or missense mutations and one 3' UTR variant of TREX1 present in 9/417 individuals with systemic lupus erythematosus but absent in 1,712 controls (P = 4.1 x 10(-7)). We demonstrate that two mutant TREX1 alleles alter subcellular targeting. Our findings implicate TREX1 in the pathogenesis of SLE.