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.

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.

Objective. STAT4 encodes a transcriptional factor that transmits signals induced by several key cytokines, and it might be a key molecule in the development of autoimmune diseases. Recently, a STAT4 haplotype was reported to be associated with rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) in Caucasian populations. This was replicated in a Korean RA population. Interestingly, the degree of risk of RA susceptibility with the STAT4 haplotype was similar in the Caucasian and Korean populations. The present study was undertaken to investigate the effect of STAT4 on susceptibility to RA and SLE in the Japanese.Methods. We performed an association study using 3 independent Japanese RA case-control populations (total 3,567 cases and 2,199 controls) and 3 independent Japanese SLE populations (total 591 cases). All samples were genotyped using the TaqMan fluorogenic 5 nuclease assay for single-nucleotide polymorphism (SNP) rs7574865, which tags the susceptibility haplotype. The association of the SNP with disease susceptibility in each case-control study was calculated using Fisher's exact test, and the results were combined, using the Mantel-Haenszel method, to obtain combined odds ratios (ORs).Results. We observed a significant association of the STAT4 polymorphism with susceptibility to both RA and SLE. The combined ORs for RA and SLE, respectively, were 1.27 (P ‫؍‬ 8.4 ؋ 10 ؊9 ) and 1.61 (P ‫؍‬ 2.1 ؋ 10 ؊11 ) for allele frequency distribution; these ORs were quite similar to those previously observed in the Caucasian population.Conclusion. We conclude that STAT4 is associated with RA and SLE in the Japanese. Our results indicate that STAT4 is a common genetic risk factor for autoimmune diseases, with similar strength across major racial groups.

Recent evidence suggests that a substantial portion of complex disease risk alleles modify gene expression in a cell-specific manner. To identify candidate causal genes and biological pathways of immune-related complex diseases, we conducted expression quantitative trait loci (eQTL) analysis on five subsets of immune cells (CD4 T cells, CD8 T cells, B cells, natural killer (NK) cells and monocytes) and unfractionated peripheral blood from 105 healthy Japanese volunteers. We developed a three-step analytical pipeline comprising (i) prediction of individual gene expression using our eQTL database and public epigenomic data, (ii) gene-level association analysis and (iii) prediction of cell-specific pathway activity by integrating the direction of eQTL effects. By applying this pipeline to rheumatoid arthritis data sets, we identified candidate causal genes and a cytokine pathway (upregulation of tumor necrosis factor (TNF) in CD4 T cells). Our approach is an efficient way to characterize the polygenic contributions and potential biological mechanisms of complex diseases.

Ossification of the posterior longitudinal ligament (OPLL) of the spine is a subset of "bone-forming" diseases, characterized by ectopic ossification in the spinal ligaments. OPLL is a common disorder among elderly populations in eastern Asia and is the leading cause of spinal myelopathy in Japan. We performed a genomewide linkage study with 142 affected sib pairs, to identify genetic loci related to OPLL. In multipoint linkage analysis using GENEHUNTER-PLUS, evidence of linkage to OPLL was detected on chromosomes 1p, 6p, 11q, 14q, 16q, and 21q. The best evidence of linkage was detected near D21S1903 on chromosome 21q22.3 (maximum Zlr=3.97); therefore, the linkage region was extensively investigated for linkage disequilibrium with single-nucleotide polymorphisms (SNPs) covering 20 Mb. One hundred fifty positional candidate genes lie in the region, and 600 gene-based SNPs were genotyped. There were positive allelic associations with seven genes (P<.01) in 280 patients and 210 controls, and four of the seven genes were clustered within a region of 750 kb, approximately 1.2 Mb telomeric to D21S1903. Extensive linkage disequilibrium and association studies of the four genes indicated that SNPs in the collagen 6A1 gene (COL6A1) were strongly associated with OPLL (P=.000003 for the SNP in intron 32 [-29]). Haplotype analysis with three SNPs in COL6A1 gave a single-point P value of.0000007. Identification of the locus of susceptibility to OPLL by genomewide linkage and linkage disequilibrium studies permits us to investigate the pathogenesis of the disease, which may lead to the development of novel therapeutic tools.

Objective. The peptidylarginine deiminase type 4 gene (PADI4) was recently reported to be associated with rheumatoid arthritis (RA) in a Japanese population. The presence of a single-nucleotide polymorphism (SNP) located in intron 3 of PADI4 provided the strongest evidence of this association. Moreover, functional haplotypes that affect stability of transcripts were identified. However, subsequent research failed to confirm the observed association in a UK population. The present study was undertaken to further investigate the association of PADI4 with RA, using a series of population-based samples from subjects with the same ethnic background as the subjects in the original study.Methods. DNA samples were obtained from 1,230 Japanese RA patients and 948 ethnically matched controls. Genotyping was performed using 5 allele discrimination assays. All samples were genotyped for 3 SNPs on PADI4 (padi4_94, padi4_104, and padi4_102), which comprised the reported haplotypes. Chi-square testing was performed for a case-control study and the PENHAPLO program was used for haplotype estimation.Results. All tested SNPs were found to show significant differences in frequency between cases and controls (P ‫؍‬ 0.010-0.0008), which confirmed the association observed in the original study. Odds ratios calculated for allele frequencies were 1.23, 1.21, and 1.36 in padi4_94, padi4_104, and padi4_102 respectively.Conclusion. Replication of association in individual samples strongly suggests that PADI4 is a true susceptibility gene for RA.