Genetics of rheumatoid arthritis contributes to biology and drug discovery.

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study 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. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2 - 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses--as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout 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.

Splicing variant of WDFY4 augments MDA5 signalling and the risk of clinically amyopathic dermatomyositis.

OBJECTIVES: Idiopathic inflammatory myopathies (IIMs) are a heterogeneous group of rare autoimmune diseases in which both genetic and environmental factors play important roles. To identify genetic factors of IIM including polymyositis, dermatomyositis (DM) and clinically amyopathic DM (CADM), we performed the first genome-wide association study for IIM in an Asian population. METHODS: We genotyped and tested 496 819 single nucleotide polymorphism for association using 576 patients with IIM and 6270 control subjects. We also examined the causal mechanism of disease-associated variants by in silico analyses using publicly available data sets as well as by in in vitro analyses using reporter assays and apoptosis assays. RESULTS: We identified a variant in WDFY4 that was significantly associated with CADM (rs7919656; OR=3.87; P=1.5×10-8). This variant had a cis-splicing quantitative trait locus (QTL) effect for a truncated WDFY4isoform (tr-WDFY4), with higher expression in the risk allele. Transexpression QTL analysis of this variant showed a positive correlation with the expression of NF-κB associated genes. Furthermore, we demonstrated that both WDFY4 and tr-WDFY4 interacted with pattern recognition receptors such as TLR3, TLR4, TLR9 and MDA5 and augmented the NF-κB activation by these receptors. WDFY4 isoforms also enhanced MDA5-induced apoptosis to a greater extent in the tr-WDFY4-transfected cells. CONCLUSIONS: As CADM is characterised by the appearance of anti-MDA5 autoantibodies and severe lung inflammation, the WDFY4 variant may play a critical role in the pathogenesis of CADM.

Functional variants in NFKBIE and RTKN2 involved in activation of the NF-κB pathway are associated with rheumatoid arthritis in Japanese.

Rheumatoid arthritis is an autoimmune disease with a complex etiology, leading to inflammation of synovial tissue and joint destruction. Through a genome-wide association study (GWAS) and two replication studies in the Japanese population (7,907 cases and 35,362 controls), we identified two gene loci associated with rheumatoid arthritis susceptibility (NFKBIE at 6p21.1, rs2233434, odds ratio (OR) = 1.20, P = 1.3 × 10(-15); RTKN2 at 10q21.2, rs3125734, OR = 1.20, P = 4.6 × 10(-9)). In addition to two functional non-synonymous SNPs in NFKBIE, we identified candidate causal SNPs with regulatory potential in NFKBIE and RTKN2 gene regions by integrating in silico analysis using public genome databases and subsequent in vitro analysis. Both of these genes are known to regulate the NF-κB pathway, and the risk alleles of the genes were implicated in the enhancement of NF-κB activity in our analyses. These results suggest that the NF-κB pathway plays a role in pathogenesis and would be a rational target for treatment of rheumatoid arthritis.

A genome-wide association study identified AFF1 as a susceptibility locus for systemic lupus eyrthematosus in Japanese.

Systemic lupus erythematosus (SLE) is an autoimmune disease that causes multiple organ damage. Although recent genome-wide association studies (GWAS) have contributed to discovery of SLE susceptibility genes, few studies has been performed in Asian populations. Here, we report a GWAS for SLE examining 891 SLE cases and 3,384 controls and multi-stage replication studies examining 1,387 SLE cases and 28,564 controls in Japanese subjects. Considering that expression quantitative trait loci (eQTLs) have been implicated in genetic risks for autoimmune diseases, we integrated an eQTL study into the results of the GWAS. We observed enrichments of cis-eQTL positive loci among the known SLE susceptibility loci (30.8%) compared to the genome-wide SNPs (6.9%). In addition, we identified a novel association of a variant in the AF4/FMR2 family, member 1 (AFF1) gene at 4q21 with SLE susceptibility (rs340630; P = 8.3×10(-9), odds ratio = 1.21). The risk A allele of rs340630 demonstrated a cis-eQTL effect on the AFF1 transcript with enhanced expression levels (P<0.05). As AFF1 transcripts were prominently expressed in CD4(+) and CD19(+) peripheral blood lymphocytes, up-regulation of AFF1 may cause the abnormality in these lymphocytes, leading to disease onset.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases.

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

An association analysis of HLA-DRB1 with systemic lupus erythematosus and rheumatoid arthritis in a Japanese population: effects of *09:01 allele on disease phenotypes.

Objective. To re-evaluate the roles of HLA-DRB1 alleles in susceptibility to SLE and RA and their effects on autoantibody status in large-scale Japanese cohorts. Methods. A total of 656 SLE, 2410 RA and 911 control subjects, who were all Japanese, were genotyped for HLA-DRB1 alleles using sequence-specific oligonucleotide probes. The association of alleles with disease susceptibility was tested by logistic regression analysis and by the relative predispositional effect method. The association with autoantibody status was examined by the standard χ(2) test. Results. HLA-DRB1*15:01, *09:01, *08:02 and *04:01 were significantly associated with SLE susceptibility, while shared epitope (SE) alleles and DRB1*09:01 were associated with RA susceptibility. The compound heterozygote of DRB1*09:01/*15:01 conferred an increased risk for SLE compared with the homozygotes for DRB1*09:01 and *15:01 and was associated with earlier onset of disease, whereas the compound effect of DRB1-SE/*09:01 was not clear in RA. DRB1*09:01 was significantly associated with the appearance of anti-Sm antibody in SLE as well as ACPA in RA, while protectively associated with anti-dsDNA antibody in SLE. No significant interaction was observed between DRB1*09:01 and smoking status for the appearance of ACPA, unlike that observed in SE alleles in RA. Conclusion. We identified HLA-DRB1 alleles associated with SLE and RA in a Japanese population and demonstrated a shared susceptibility of DRB1*09:01 between the diseases as well as its effect on autoantibody production.

Upregulated Fcrl5 disrupts B cell anergy and causes autoimmune disease.

B cell anergy plays a critical role in maintaining self-tolerance by inhibiting autoreactive B cell activation to prevent autoimmune diseases. Here, we demonstrated that Fc receptor-like 5 (Fcrl5) upregulation contributes to autoimmune disease pathogenesis by disrupting B cell anergy. Fcrl5-a gene whose homologs are associated with human autoimmune diseases-is highly expressed in age/autoimmunity-associated B cells (ABCs), an autoreactive B cell subset. By generating B cell-specific Fcrl5 transgenic mice, we demonstrated that Fcrl5 overexpression in B cells caused systemic autoimmunity with age. Additionally, Fcrl5 upregulation in B cells exacerbated the systemic lupus erythematosus-like disease model. Furthermore, an increase in Fcrl5 expression broke B cell anergy and facilitated toll-like receptor signaling. Thus, Fcrl5 is a potential regulator of B cell-mediated autoimmunity by regulating B cell anergy. This study provides important insights into the role of Fcrl5 in breaking B cell anergy and its effect on the pathogenesis of autoimmune diseases.

Regulatory polymorphisms in EGR2 are associated with susceptibility to systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is an autoimmune disease induced by the combinations of environmental and genetic factors. Recently, mice in which the early growth response 2 (EGR2) gene, a zinc-finger transcription factor, is conditionally knocked out in CD2(+) T cells have been shown to develop a lupus-like autoimmune disease. Here, we evaluated if polymorphisms in the EGR2 gene influence SLE susceptibility in humans. We first analyzed the effect of SNPs in the EGR2 region on EGR2 expression, and a significant positive correlation with expression was identified in an SNP located at the 5' flanking region of EGR2 (rs10761670, R=0.23, P=0.00072). We then performed a case-control association study using three sets of SLE cohorts by genotyping 14 tag SNPs in the EGR2 gene region. A peak of association with SLE susceptibility was observed for rs10761670 [Pooled: OR = 1.23 (95% CI 1.10-1.37), P=0.00023). This SNP was also associated with susceptibility to rheumatoid arthritis (RA) [OR = 1.15 (95% CI 1.05-1.26), P = 0.0019), suggesting that EGR2 is a common risk factor for SLE and RA. Among the SNPs in complete linkage disequilibrium with rs10761670 (r(2) = 1.0), two SNPs (rs1412554 and rs1509957) affected the binding of transcription factors and transcriptional activity in vitro, suggesting that they may be candidates of causal regulatory variants in this region. Therefore, EGR2 is a genetic risk factor for SLE, in which increased gene expression may contribute to SLE pathogenesis.

PADI4 polymorphism predisposes male smokers to rheumatoid arthritis.

OBJECTIVE: To elucidate the differential role of peptidyl arginine deiminase 4 (PADI4) polymorphism in rheumatoid arthritis (RA) between Asian and European populations, possible gene-environmental interactions among the PADI4 polymorphism, sex and smoking status were analysed. METHODS: Three independent sets of case-control samples were genotyped for single-nucleotide polymorphisms in PADI4; Japanese samples (first set, 1019 RA patients, 907 controls; second set, 999 RA patients, 1128 controls) using TaqMan assays and Dutch samples (635 RA patients, 391 controls) using Sequenom MassARRAY platform. The association of PADI4 with RA susceptibility was evaluated by smoking status and sex in contingency tables and logistic regression models. RESULTS: In the first set of Japanese samples, PADI4 polymorphism (rs1748033) showed a greater risk in men (OR(allele) 1.39; 95% CI 1.10 to 1.76; p(trend)=0.0054) than in women and in ever-smokers (OR(allele) 1.25; 95% CI 1.02 to 1.53; p(trend)=0.032) than in never-smokers. Moreover, the highest risk was seen in male ever-smokers (OR(allele) 1.46; 95% CI 1.12 to 1.90; p(trend)=0.0047). Similar trends were observed in the second set of Japanese samples as well as in Dutch samples. CONCLUSION: PADI4 polymorphism highly predisposes male smokers to RA, and the genetic heterogeneity observed between Asian and European populations may be partly explained by differences in smoking prevalence among men.

FCRL3, an autoimmune susceptibility gene, has inhibitory potential on B-cell receptor-mediated signaling.

A polymorphism that up-regulates the expression of Fc receptor-like 3 (FCRL3) gene has recently been described as predisposing for several human autoimmune diseases. FCRL3 is preferentially expressed on B cells and is unique in displaying both an ITAM and an ITIM in the cytosolic domain, suggesting signaling functions. Herein, we show that FCRL3 potentially inhibits BCR-mediated signaling, using murine FcgammaRIIB/human FCRL3 chimeric protein. Coligation of the chimeric protein with BCR leads to phosphorylation of tyrosine residues in the cytosolic domain. This coligation inhibits cell tyrosine phosphorylation and calcium mobilization in addition to activation-induced cell death mediated by BCR signaling. Mutational analysis showed the tyrosine residues in two potential ITIMs at 662 and 692 offer the main contributions to this inhibition, which is further supported by strong associations of SH-2 domain-containing phosphatases with the following phosphotyrosine motifs: SHIP with the ITIM-like motif at 662; and SHP-1 and -2 with the canonical ITIM at 692. These results, together with previous genetic data, suggest that augmented inhibition of BCR-mediated signaling by FCRL3 with the disease-risk genotype alter the activation threshold and promote tolerance breakdown in B cells.

Anti-citrullinated peptide/protein antibody (ACPA)-negative RA shares a large proportion of susceptibility loci with ACPA-positive RA: a meta-analysis of genome-wide association study in a Japanese population.

INTRODUCTION: Although susceptibility genes for anti-citrullinated peptide/protein antibodies (ACPA)-positive rheumatoid arthritis (RA) have been successfully discovered by genome-wide association studies (GWAS), little is known about the genetic background of ACPA-negative RA. We intended to elucidate genetic background of ACPA-negative RA. METHOD: We performed a meta-analysis of GWAS comprising 670 ACPA-negative RA and 16,891 controls for 1,948,138 markers, followed by a replication study of the top 35 single nucleotide polymorphisms (SNPs) using 916 cases and 3,764 controls. Inverse-variance method was applied to assess overall effects. To assess overlap of susceptibility loci between ACPA-positive and -negative RA, odds ratios (ORs) of the 21 susceptibility markers to RA in Japanese were compared between the two subsets. In addition, SNPs were stratified by the p-values in GWAS meta-analysis for either ACPA-positive RA or ACPA-negative RA to address the question whether weakly-associated genes were also shared. The correlations between ACPA-positive RA and the subpopulations of ACPA-negative RA (rheumatoid factor (RF)-positive and RF-negative subsets) were also addressed. RESULTS: Rs6904716 in LEMD2 of the human leukocyte antigen (HLA) locus showed a borderline association with ACPA-negative RA (overall p = 5.7 × 10(-8)), followed by rs6986423 in CSMD1 (p = 2.4 × 10(-6)) and rs17727339 in FCRL3 (p = 1.4 × 10(-5)). ACPA-negative RA showed significant correlations of ORs with ACPA-positive RA for the 21 susceptibility SNPs and non-HLA SNPs with p-values far from significance. These significant correlations with ACPA-positive RA were true for ACPA-negative RF-positive and ACPA-negative RF-negative RA. On the contrary, positive correlations were not observed between the ACPA-negative two subpopulations. CONCLUSION: Many of the susceptibility loci were shared between ACPA-positive and -negative RA.

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.

The association of a nonsynonymous single-nucleotide polymorphism in TNFAIP3 with systemic lupus erythematosus and rheumatoid arthritis in the Japanese population.

OBJECTIVE: Genome-wide association (GWA) studies in systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) in Caucasian populations have independently identified risk variants in and near the tumor necrosis factor alpha (TNFalpha)-induced protein 3 gene (TNFAIP3), which is crucial for the regulation of TNF-mediated signaling and Toll-like receptor signaling. The aim of this study was to assess the role of TNFAIP3 in the development of SLE and RA in Japanese subjects. METHODS: We selected 2 single-nucleotide polymorphisms (SNPs) from previous GWA studies. Rs2230926 is a nonsynonymous SNP in TNFAIP3 and is associated with SLE, while rs10499194 is an intergenic SNP associated with RA. We then performed 2 independent sets of SLE case-control comparisons (717 patients and 1,362 control subjects) and 3 sets of RA case-control comparisons (3,446 patients and 2,344 control subjects) using Japanese subjects. We genotyped SNPs using TaqMan assays. RESULTS: We observed a significant association between rs2230926 and an increased risk of SLE and RA in the Japanese population (for SLE, odds ratio [OR] 1.92, 95% confidence interval [95% CI] 1.53-2.41, P = 1.9 x 10(-8); for RA, OR 1.35, 95% CI 1.18-1.56, P = 2.6 x 10(-5)). The intergenic SNP rs10499194 was also associated with SLE and RA, while the risk allele for RA in Caucasians was protective against the diseases in our population. CONCLUSION: We demonstrated a significant association between the nonsynonymous variant in TNFAIP3 and the risk for SLE and RA in the Japanese population. TNFAIP3, similar to STAT4 and IRF5, may be a common genetic risk factor for SLE and RA that is shared between the Caucasian and Japanese populations.

A regulatory variant in CCR6 is associated with rheumatoid arthritis susceptibility.

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.

Contribution of a haplotype in the HLA region to anti-cyclic citrullinated peptide antibody positivity in rheumatoid arthritis, independently of HLA-DRB1.

OBJECTIVE: To examine the risk of anti-cyclic citrullinated peptide (anti-CCP) antibody positivity in rheumatoid arthritis (RA) patients carrying certain haplotypes in the HLA region. METHODS: A total of 1,389 Japanese patients with RA were genotyped for 30 single-nucleotide polymorphisms (SNPs) in the HLA region using commercial oligonucleotide arrays (from Perlegen or Affymetrix) as well as for HLA-DRB1 alleles using a sequence-specific polymerase chain reaction method. Stepwise logistic regression was used to select from among the 30 SNPs the ones that represented a risk of anti-CCP antibody positivity. Haplotypes of the selected SNPs were inferred using an expectation-maximization algorithm. Associations of individual SNPs were evaluated with the Cochran-Armitage test for trend. DRB1 alleles and haplotypes were evaluated with the chi-square test. Heterogeneities of risks among the shared epitope (SE) and non-SE HLA-DRB1 alleles were examined using the exact test. Haplotype associations that were independent of individual HLA-DRB1 alleles were evaluated using the likelihood ratio test. RESULTS: Significant associations were found for 9 SNPs (smallest P value being 2.4x10(-8)) and in 4 HLA-DRB1 alleles (smallest P value being 2.0x10(-10) in DRB1*0405). Stepwise logistic regression selected 4 SNPs (rs9262638, rs7775228, rs4713580, and rs9277359). Among the 16 inferred haplotypes of these 4 SNPs, 6 indicated significant associations (smallest P value being 1.9x10(-11)). Risks among SE and non-SE alleles were significantly heterogeneous (P=0.0095 and P=9.8x10(-9), respectively), indicating the importance of stratification with individual DRB1 alleles rather than SE alleles. Conditional analysis of the risk associated with individual DRB1 alleles identified a risk haplotype that was independent of DRB1 (odds ratio 2.00 [95% confidence interval 1.44-2.79], P=2.6x10(-5)). CONCLUSION: Heterogeneous risks of anti-CCP antibody positivity were confirmed among SE and non-SE alleles in our patient population. A risk haplotype in the HLA region that is independent of HLA-DRB1 was confirmed.