A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological datasets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA)1. Here, we performed a genome-wide association study (GWAS) meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ~10 million single nucleotide polymorphisms (SNPs). We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 1012–4. We devised an in-silico pipeline using established bioinformatics methods based on functional annotation5, cis-acting expression quantitative trait loci (cis-eQTL)6, and pathway analyses7–9 – as well as novel methods based on genetic overlap with human primary immunodeficiency (PID), hematological cancer somatic mutations and knock-out mouse phenotypes – to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

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.

Rheumatoid arthritis is a common autoimmune disease with a complex genetic etiology. Here, through a genome-wide association study of rheumatoid arthritis, we identified a polymorphism in CCR6, the gene encoding chemokine (C-C motif) receptor 6 (a surface marker for Th17 cells) at 6q27, that was associated with rheumatoid arthritis susceptibility and was validated in two independent replication cohorts from Japan (rs3093024, a total of 7,069 individuals with rheumatoid arthritis (cases) and 20,727 controls, overall odds ratio = 1.19, P = 7.7 x 10(-19)). We identified a triallelic dinucleotide polymorphism of CCR6 (CCR6DNP) in strong linkage disequilibrium with rs3093024 that showed effects on gene transcription. The CCR6DNP genotype was correlated with the expression level of CCR6 and was associated with the presence of interleukin-17 (IL-17) in the sera of subjects with rheumatoid arthritis. Moreover, CCR6DNP was associated with susceptibility to Graves' and Crohn's diseases. These results suggest that CCR6 is critically involved in IL-17-driven autoimmunity in human diseases.

Rheumatoid arthritis (RA) is a chronic symmetric polyarticular joint disease that primarily affects the small joints of the hands and feet. The inflammatory process is characterized by infiltration of inflammatory cells into the joints, leading to proliferation of synoviocytes and destruction of cartilage and bone. In RA synovial tissue, the infiltrating cells such as macrophages, T cells, B cells and dendritic cells play important role in the pathogenesis of RA. Migration of leukocytes into the synovium is a regulated multi‐step process, involving interactions between leukocytes and endothelial cells, cellular adhesion molecules, as well as chemokines and chemokine receptors. Chemokines are small, chemoattractant cytokines which play key roles in the accumulation of inflammatory cells at the site of inflammation. It is known that synovial tissue and synovial fluid from RA patients contain increased concentrations of several chemokines, such as monocyte chemoattractant protein‐4 (MCP‐4)/CCL13, pulmonary and activation‐regulated chemokine (PARC)/CCL18, monokine induced by interferon‐γ (Mig)/CXCL9, stromal cell‐derived factor 1 (SDF‐1)/CXCL12, monocyte chemotactic protein 1 (MCP‐1)/CCL2, macrophage inflammatory protein 1α (MIP‐1α)/CCL3, and Fractalkine/CXC3CL1. Therefore, chemokines and chemokine‐receptors are considered to be important molecules in RA pathology.

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.