BACKGROUND: Juvenile idiopathic arthritis (JIA) is the most common chronic rheumatic condition of childhood. Genetic association studies have revealed several JIA susceptibility loci with the strongest effect size observed in the human leukocyte antigen (HLA) region. Genome-wide association studies have augmented the number of JIA-associated loci, particularly for non-HLA genes. The aim of this study was to identify new associations at non-HLA loci predisposing to the risk of JIA development in Estonian patients. METHODS: We performed genome-wide association analyses in an entire JIA case-control sample (All-JIA) and in a case-control sample for oligoarticular JIA, the most prevalent JIA subtype. The entire cohort was genotyped using the Illumina HumanOmniExpress BeadChip arrays. After imputation, 16,583,468 variants were analyzed in 263 cases and 6956 controls. RESULTS: We demonstrated nominal evidence of association for 12 novel non-HLA loci not previously implicated in JIA predisposition. We replicated known JIA associations in CLEC16A and VCTN1 regions in the oligoarticular JIA sample. The strongest associations in the All-JIA analysis were identified at PRKG1 (P = 2,54 × 10-6), LTBP1 (P = 9,45 × 10-6), and ELMO1 (P = 1,05 × 10-5). In the oligoarticular JIA analysis, the strongest associations were identified at NFIA (P = 5,05 × 10-6), LTBP1 (P = 9,95 × 10-6), MX1 (P = 1,65 × 10-5), and CD200R1 (P = 2,59 × 10-5). CONCLUSION: This study increases the number of known JIA risk loci and provides additional evidence for the existence of overlapping genetic risk loci between JIA and other autoimmune diseases, particularly rheumatoid arthritis. The reported loci are involved in molecular pathways of immunological relevance and likely represent genomic regions that confer susceptibility to JIA in Estonian patients. Key Points • Juvenile idiopathic arthritis (JIA) is the most common childhood rheumatic disease with heterogeneous presentation and genetic predisposition. • Present genome-wide association study for Estonian JIA patients is first of its kind in Northern and Northeastern Europe. • The results of the present study increase the knowledge about JIA risk loci replicating some previously described associations, so adding weight to their relevance and describing novel loci. • The study provides additional evidence for the existence of overlapping genetic risk loci between JIA and other autoimmune diseases, particularly rheumatoid arthritis.

OBJECTIVES: Nearly 110 susceptibility loci for rheumatoid arthritis (RA) with modest effect sizes have been identified by population-based genetic association studies, suggesting a large number of undiscovered variants behind a highly polygenic genetic architecture of RA. Here, we performed the largest-ever trans-ancestral meta-analysis with the aim to identify new RA loci and to better understand RA biology underlying genetic associations. METHODS: Genome-wide RA association summary statistics in three large case-control collections consisting of 311 292 individuals of Korean, Japanese and European populations were used in an inverse-variance-weighted fixed-effects meta-analysis. Several computational analyses using public omics resources were conducted to prioritise causal variants and genes, RA variant-implicating features (tissues, pathways and transcription factors) and potentially repurposable drugs for RA treatment. RESULTS: We identified 11 new RA susceptibility loci that explained 6.9% and 1.8% of the single-nucleotide polymorphism-based heritability in East Asians and Europeans, respectively, and confirmed 71 known non-human leukocyte antigens (HLA) susceptibility loci, identifying 90 independent association signals. The RA variants were preferentially located in binding sites of various transcription factors and in cell type-specific transcription-activation histone marks that simultaneously highlighted the importance of CD4+ T-cell activation and the potential role of non-immune organs in RA pathogenesis. A total of 615 plausible effector genes, based on gene-based associations, expression-associated variants and chromatin interaction, included targets of drugs approved for RA treatments and potentially repurposable drugs approved for other indications. CONCLUSION: Our findings provide useful insights regarding RA genetic aetiology and variant-driven RA pathogenesis.

OBJECTIVE: Genome-wide association studies (GWAS) in rheumatoid arthritis (RA) have discovered over 100 RA loci, explaining patient-relevant RA pathogenesis but showing a large fraction of missing heritability. As a continuous effort, we conducted GWAS in a large Korean RA case-control population. METHODS: We newly generated genome-wide variant data in two independent Korean cohorts comprising 4068 RA cases and 36 487 controls, followed by a whole-genome imputation and a meta-analysis of the disease association results in the two cohorts. By integrating publicly available omics data with the GWAS results, a series of bioinformatic analyses were conducted to prioritise the RA-risk genes in RA loci and to dissect biological mechanisms underlying disease associations. RESULTS: We identified six new RA-risk loci (SLAMF6, CXCL13, SWAP70, NFKBIA, ZFP36L1 and LINC00158) with pmeta<5×10-8 and consistent disease effect sizes in the two cohorts. A total of 122 genes were prioritised from the 6 novel and 13 replicated RA loci based on physical distance, regulatory variants and chromatin interaction. Bioinformatics analyses highlighted potentially RA-relevant tissues (including immune tissues, lung and small intestine) with tissue-specific expression of RA-associated genes and suggested the immune-related gene sets (such as CD40 pathway, IL-21-mediated pathway and citrullination) and the risk-allele sharing with other diseases. CONCLUSION: This study identified six new RA-associated loci that contributed to better understanding of the genetic aetiology and biology in RA.

Aim: To identify novel genetic variants predisposing to elevation of Alanine aminotransferase (ALT) in rheumatoid arthritis (RA) patients after initiation of methotrexate (MTX) treatment. Patients & methods: We performed genome-wide association studies in 198 RA patients starting MTX. Outcomes were maximum level of ALT and ALT >1.5-times the upper level of normal within the first 6 months of treatment. Results: RAVER2 (rs72675408) was significantly associated with maximum level of ALT (p = 4.36 × 10-8). This variant is in linkage disequilibrium with rs72675451, which is associated with differential expression of JAK1 and RAVER2. Conclusion: We found an association between ALT elevation and genetic variants that may regulate the expression of JAK1 and RAVER2. JAK1 encodes a janus kinase involved in the pathogenesis of RA.

The overwhelming majority of participants in current genetic studies are of European ancestry. To elucidate disease biology in the East Asian population, we conducted a genome-wide association study (GWAS) with 212,453 Japanese individuals across 42 diseases. We detected 320 independent signals in 276 loci for 27 diseases, with 25 novel loci (P < 9.58 × 10-9). East Asian-specific missense variants were identified as candidate causal variants for three novel loci, and we successfully replicated two of them by analyzing independent Japanese cohorts; p.R220W of ATG16L2 (associated with coronary artery disease) and p.V326A of POT1 (associated with lung cancer). We further investigated enrichment of heritability within 2,868 annotations of genome-wide transcription factor occupancy, and identified 378 significant enrichments across nine diseases (false discovery rate < 0.05) (for example, NKX3-1 for prostate cancer). This large-scale GWAS in a Japanese population provides insights into the etiology of complex diseases and highlights the importance of performing GWAS in non-European populations.

BACKGROUND: Biopharmaceutical products (BPs) are widely used to treat autoimmune diseases, but immunogenicity limits their efficacy for an important proportion of patients. Our knowledge of patient-related factors influencing the occurrence of antidrug antibodies (ADAs) is still limited. METHODS AND FINDINGS: The European consortium ABIRISK (Anti-Biopharmaceutical Immunization: prediction and analysis of clinical relevance to minimize the RISK) conducted a clinical and genomic multicohort prospective study of 560 patients with multiple sclerosis (MS, n = 147), rheumatoid arthritis (RA, n = 229), Crohn's disease (n = 148), or ulcerative colitis (n = 36) treated with 8 different biopharmaceuticals (etanercept, n = 84; infliximab, n = 101; adalimumab, n = 153; interferon [IFN]-beta-1a intramuscularly [IM], n = 38; IFN-beta-1a subcutaneously [SC], n = 68; IFN-beta-1b SC, n = 41; rituximab, n = 31; tocilizumab, n = 44) and followed during the first 12 months of therapy for time to ADA development. From the bioclinical data collected, we explored the relationships between patient-related factors and the occurrence of ADAs. Both baseline and time-dependent factors such as concomitant medications were analyzed using Cox proportional hazard regression models. Mean age and disease duration were 35.1 and 0.85 years, respectively, for MS; 54.2 and 3.17 years for RA; and 36.9 and 3.69 years for inflammatory bowel diseases (IBDs). In a multivariate Cox regression model including each of the clinical and genetic factors mentioned hereafter, among the clinical factors, immunosuppressants (adjusted hazard ratio [aHR] = 0.408 [95% confidence interval (CI) 0.253-0.657], p < 0.001) and antibiotics (aHR = 0.121 [0.0437-0.333], p < 0.0001) were independently negatively associated with time to ADA development, whereas infections during the study (aHR = 2.757 [1.616-4.704], p < 0.001) and tobacco smoking (aHR = 2.150 [1.319-3.503], p < 0.01) were positively associated. 351,824 Single-Nucleotide Polymorphisms (SNPs) and 38 imputed Human Leukocyte Antigen (HLA) alleles were analyzed through a genome-wide association study. We found that the HLA-DQA1*05 allele significantly increased the rate of immunogenicity (aHR = 3.9 [1.923-5.976], p < 0.0001 for the homozygotes). Among the 6 genetic variants selected at a 20% false discovery rate (FDR) threshold, the minor allele of rs10508884, which is situated in an intron of the CXCL12 gene, increased the rate of immunogenicity (aHR = 3.804 [2.139-6.764], p < 1 × 10-5 for patients homozygous for the minor allele) and was chosen for validation through a CXCL12 protein enzyme-linked immunosorbent assay (ELISA) on patient serum at baseline before therapy start. CXCL12 protein levels were higher for patients homozygous for the minor allele carrying higher ADA risk (mean: 2,693 pg/ml) than for the other genotypes (mean: 2,317 pg/ml; p = 0.014), and patients with CXCL12 levels above the median in serum were more prone to develop ADAs (aHR = 2.329 [1.106-4.90], p = 0.026). A limitation of the study is the lack of replication; therefore, other studies are required to confirm our findings. CONCLUSION: In our study, we found that immunosuppressants and antibiotics were associated with decreased risk of ADA development, whereas tobacco smoking and infections during the study were associated with increased risk. We found that the HLA-DQA1*05 allele was associated with an increased rate of immunogenicity. Moreover, our results suggest a relationship between CXCL12 production and ADA development independent of the disease, which is consistent with its known function in affinity maturation of antibodies and plasma cell survival. Our findings may help physicians in the management of patients receiving biotherapies.

Certolizumab pegol (CZP) is a PEGylated Fc-free tumor necrosis factor (TNF) inhibitor antibody approved for use in the treatment of rheumatoid arthritis (RA), Crohn’s disease, psoriatic arthritis, axial spondyloarthritis and psoriasis. In a clinical trial of patients with severe RA, CZP improved disease symptoms in approximately half of patients. However, variability in CZP efficacy remains a problem for clinicians, thus, the aim of this study was to identify genetic variants predictive of CZP response. We performed a genome-wide association study (GWAS) of 302 RA patients treated with CZP in the REALISTIC trial to identify common single nucleotide polymorphisms (SNPs) associated with treatment response. Whole-exome sequencing was also performed for 74 CZP extreme responders and non-responders within the same population, as well as 1546 population controls. No common SNPs or rare functional variants were significantly associated with CZP response, though a non-significant enrichment in the RA-implicated KCNK5 gene was observed. Two SNPs near spondin-1 and semaphorin-4G approached genome-wide significance. The results of the current study did not provide an unambiguous predictor of CZP response.

BACKGROUND: In recent years, research has consistently proven the occurrence of genetic overlap across autoimmune diseases, which supports the existence of common pathogenic mechanisms in autoimmunity. The objective of this study was to further investigate this shared genetic component. METHODS: For this purpose, we performed a cross-disease meta-analysis of Immunochip data from 37,159 patients diagnosed with a seropositive autoimmune disease (11,489 celiac disease (CeD), 15,523 rheumatoid arthritis (RA), 3477 systemic sclerosis (SSc), and 6670 type 1 diabetes (T1D)) and 22,308 healthy controls of European origin using the R package ASSET. RESULTS: We identified 38 risk variants shared by at least two of the conditions analyzed, five of which represent new pleiotropic loci in autoimmunity. We also identified six novel genome-wide associations for the diseases studied. Cell-specific functional annotations and biological pathway enrichment analyses suggested that pleiotropic variants may act by deregulating gene expression in different subsets of T cells, especially Th17 and regulatory T cells. Finally, drug repositioning analysis evidenced several drugs that could represent promising candidates for CeD, RA, SSc, and T1D treatment. CONCLUSIONS: In this study, we have been able to advance in the knowledge of the genetic overlap existing in autoimmunity, thus shedding light on common molecular mechanisms of disease and suggesting novel drug targets that could be explored for the treatment of the autoimmune diseases studied.

Hundreds of genomic loci have been associated with a significant number of immune-mediated diseases, and a large proportion of these associated loci are shared among traits. Both the molecular mechanisms by which these loci confer disease susceptibility and the extent to which shared loci are implicated in a common pathogenesis are unknown. We therefore sought to dissect the functional components at loci shared between two autoimmune diseases: coeliac disease (CeD) and rheumatoid arthritis (RA). We used a cohort of 12 381 CeD cases and 7827 controls, and another cohort of 13 819 RA cases and 12 897 controls, all genotyped with the Immunochip platform. In the joint analysis, we replicated 19 previously identified loci shared by CeD and RA and discovered five new non-HLA loci shared by CeD and RA. Our fine-mapping results indicate that in nine of 24 shared loci the associated variants are distinct in the two diseases. Using cell-type-specific histone markers, we observed that loci which pointed to the same variants in both diseases were enriched for marks of promoters active in CD14+ and CD34+ immune cells (P < 0.001), while loci pointing to distinct variants in one of the two diseases showed enrichment for marks of more specialized cell types, like CD4+ regulatory T cells in CeD (P < 0.0001) compared with Th17 and CD15+ in RA (P = 0.0029).

Rheumatoid arthritis (RA) is characterised by chronic synovial joint inflammation. Treatment has been revolutionised by tumour necrosis factor alpha inhibitors (TNFi) but each available drug shows a significant non-response rate. We conducted a genome-wide association study of 1752 UK RA TNFi-treated patients to identify predictors of change in the Disease Activity Score 28 (DAS28) and subcomponents over 3-6 months. The rs7195994 variant at the FTO gene locus was associated with infliximab response when looking at a change in the swollen joint count (SJC28) subcomponent (p = 9.74 × 10-9). Capture Hi-C data show chromatin interactions in GM12878 cells between rs2540767, in high linkage disequilibrium with rs7195994 (R2 = 0.9) and IRX3, a neighbouring gene of FTO. IRX3 encodes a transcription factor involved in adipocyte remodelling and is regarded as the obesity gene at the FTO locus. Importantly, the rs7195994 association remained significantly associated following adjustment for BMI. In addition, using capture Hi-C data we showed interactions between TNFi-response associated variants and 16 RA susceptibility variants.

OBJECTIVE: To investigate the gene-environment interaction between smoking and single nucleotide polymorphisms (SNPs), using Immunochip material, on the risk of developing either of two serologically defined subsets of RA. METHODS: Interaction between smoking and 133,648 genetic markers from the Immunochip was examined for two RA subsets, defined by the presence or absence of ACPA. A total of 1590 ACPA-positive and 891 ACPA-negative cases were compared with 1856 controls in the Swedish Epidemiological Investigation of RA (EIRA) case-control study. Logistic regression models were used to determine the presence of interaction. The proportion attributable to interaction was calculated for each smoking-SNP pair. Replication was carried out in an independent dataset from northern Sweden. To further validate and extend the results, interaction analysis was also performed using genome-wide association studies data on EIRA individuals. RESULTS: In ACPA-positive RA, 102 SNPs interacted significantly with smoking, after Bonferroni correction. All 102 SNPs were located in the HLA region, mainly within the HLA class II region, 51 of which were replicated. No additional loci outside chromosome 6 were identified in the genome-wide association studies validation. After adjusting for HLA-DRB1 shared epitope, 15 smoking-SNP pairs remained significant for ACPA-positive RA, with 8 of these replicated (loci: BTNL2, HLA-DRA, HLA-DRB5, HLA-DQA1, HLA-DOB and TAP2). For ACPA-negative RA, no smoking-SNP pairs passed the threshold for significance. CONCLUSION: Our study presents extended gene variation patterns involved in gene-smoking interaction in ACPA-positive, but not ACPA-negative, RA. Notably, variants in HLA-DRB1 and those in additional genes within the MHC class II region, but not in any other gene regions, showed interaction with smoking.

The human genome, which includes thousands of genes, represents a big data challenge. Rheumatoid arthritis (RA) is a complex autoimmune disease with a genetic basis. Many single-nucleotide polymorphism (SNP) association methods partition a genome into haplotype blocks. The aim of this genome wide association study (GWAS) was to select the most appropriate haplotype block partitioning method for the North American Rheumatoid Arthritis Consortium (NARAC) dataset. The methods used for the NARAC dataset were the individual SNP approach and the following haplotype block methods: the four-gamete test (FGT), confidence interval test (CIT), and solid spine of linkage disequilibrium (SSLD). The measured parameters that reflect the strength of the association between the biomarker and RA were the P-value after Bonferroni correction and other parameters used to compare the output of each haplotype block method. This work presents a comparison among the individual SNP approach and the three haplotype block methods to select the method that can detect all the significant SNPs when applied alone. The GWAS results from the NARAC dataset obtained with the different methods are presented. The individual SNP, CIT, FGT, and SSLD methods detected 541, 1516, 1551, and 1831 RA-associated SNPs respectively, and the individual SNP, FGT, CIT, and SSLD methods detected 65, 156, 159, and 450 significant SNPs respectively, that were not detected by the other methods. Three hundred eighty-three SNPs were discovered by the haplotype block methods and the individual SNP approach, while 1021 SNPs were discovered by all three haplotype block methods. The 383 SNPs detected by all the methods are promising candidates for studying RA susceptibility. A hybrid technique involving all four methods should be applied to detect the significant SNPs associated with RA in the NARAC dataset, but the SSLD method may be preferred because of its advantages when only one method was used.

OBJECTIVE: To investigate the genetic background influencing the development of cardiovascular (CV) disease in patients with rheumatoid arthritis (RA). METHODS: We performed a genome-wide association study (GWAS) in which, after quality control and imputation, a total of 6,308,944 polymorphisms across the whole genome were analyzed in 2,989 RA patients of European origin. Data on subclinical atherosclerosis, obtained through assessment of carotid intima-media thickness (CIMT) and presence/absence of carotid plaques by carotid ultrasonography, were available for 1,355 individuals. RESULTS: A genetic variant of the RARB gene (rs116199914) was associated with CIMT values at the genome-wide level of significance (minor allele [G] β coefficient 0.142, P = 1.86 × 10-8 ). Interestingly, rs116199914 overlapped with regulatory elements in tissues related to CV pathophysiology and immune cells. In addition, biologic pathway enrichment and predictive protein-protein relationship analyses, including suggestive GWAS signals of potential relevance, revealed a functional enrichment of the collagen biosynthesis network related to the presence/absence of carotid plaques (Gene Ontology no. 0032964; false discovery rate-adjusted P = 4.01 × 10-3 ). Furthermore, our data suggest potential influences of the previously described candidate CV risk loci NFKB1, MSRA, and ZC3HC1 (P = 8.12 × 10-4 , P = 5.94 × 10-4 , and P = 2.46 × 10-4 , respectively). CONCLUSION: The present findings strongly suggest that genetic variation within RARB contributes to the development of subclinical atherosclerosis in patients with RA.

Using the Immunochip custom SNP array, which was designed for dense genotyping of 186 loci identified through genome-wide association studies (GWAS), we analyzed 11,475 individuals with rheumatoid arthritis (cases) of European ancestry and 15,870 controls for 129,464 markers. We combined these data in a meta-analysis with GWAS data from additional independent cases (n = 2,363) and controls (n = 17,872). We identified 14 new susceptibility loci, 9 of which were associated with rheumatoid arthritis overall and five of which were specifically associated with disease that was positive for anticitrullinated peptide antibodies, bringing the number of confirmed rheumatoid arthritis risk loci in individuals of European ancestry to 46. We refined the peak of association to a single gene for 19 loci, identified secondary independent effects at 6 loci and identified association to low-frequency variants at 4 loci. Bioinformatic analyses generated strong hypotheses for the causal SNP at seven loci. This study illustrates the advantages of dense SNP mapping analysis to inform subsequent functional investigations.

IMPORTANCE: Posttraumatic stress disorder (PTSD) is a prevalent, serious public health concern, particularly in the military. The identification of genetic risk factors for PTSD may provide important insights into the biological foundation of vulnerability and comorbidity. OBJECTIVE: To discover genetic loci associated with the lifetime risk for PTSD in 2 cohorts from the Army Study to Assess Risk and Resilience in Servicemembers (Army STARRS). DESIGN, SETTING, AND PARTICIPANTS: Two coordinated genome-wide association studies of mental health in the US military contributed participants. The New Soldier Study (NSS) included 3167 unique participants with PTSD and 4607 trauma-exposed control individuals; the Pre/Post Deployment Study (PPDS) included 947 unique participants with PTSD and 4969 trauma-exposed controls. The NSS data were collected from February 1, 2011, to November 30, 2012; the PDDS data, from January 9 to April 30, 2012. The primary analysis compared lifetime DSM-IV PTSD cases with trauma-exposed controls without lifetime PTSD. Data were analyzed from March 18 to December 27, 2015. MAIN OUTCOMES AND MEASURES: Association analyses for PTSD used logistic regression models within each of 3 ancestral groups (European, African, and Latino American) by study, followed by meta-analysis. Heritability and genetic correlation and pleiotropy with other psychiatric and immune-related disorders were estimated. RESULTS: The NSS population was 80.7% male (6277 of 7774 participants; mean [SD] age, 20.9 [3.3] years); the PPDS population, 94.4% male (5583 of 5916 participants; mean [SD] age, 26.5 [6.0] years). A genome-wide significant locus was found in ANKRD55 on chromosome 5 (rs159572; odds ratio [OR], 1.62; 95% CI, 1.37-1.92; P = 2.34 × 10-8) and persisted after adjustment for cumulative trauma exposure (adjusted OR, 1.64; 95% CI, 1.39-1.95; P = 1.18 × 10-8) in the African American samples from the NSS. A genome-wide significant locus was also found in or near ZNF626 on chromosome 19 (rs11085374; OR, 0.77; 95% CI, 0.70-0.85; P = 4.59 × 10-8) in the European American samples from the NSS. Similar results were not found for either single-nucleotide polymorphism in the corresponding ancestry group from the PPDS sample, in other ancestral groups, or in transancestral meta-analyses. Single-nucleotide polymorphism-based heritability was nonsignificant, and no significant genetic correlations were observed between PTSD and 6 mental disorders or 9 immune-related disorders. Significant evidence of pleiotropy was observed between PTSD and rheumatoid arthritis and, to a lesser extent, psoriasis. CONCLUSIONS AND RELEVANCE: In the largest genome-wide association study of PTSD to date, involving a US military sample, limited evidence of association for specific loci was found. Further efforts are needed to replicate the genome-wide significant association with ANKRD55-associated in prior research with several autoimmune and inflammatory disorders-and to clarify the nature of the genetic overlap observed between PTSD and rheumatoid arthritis and psoriasis.

AIM: Bone erosion is a major problem worsening quality of rheumatoid arthritis (RA) patients' lives. However, causal factors responsible for bone erosion in RA have remained unclear. We aimed to examine genetic variants conferring bone erosion in RA using a Korean genome-wide association study (GWAS) and to search for possible biological mechanisms underlying the development of bone erosion. METHOD: We obtained genome-wide single nucleotide polymorphism (SNP) data for 711 Korean RA patients using Illumina HapMap 550v3/660W arrays. Associations between SNPs and bone erosion status based on the Steinbrocker staging system were examined using multivariate logistic regression. Cell-type-specific enrichment of the epigenomic chromatin annotation H3K4me3 at the bone erosion associated variants was further investigated using National Institute of Health Roadmap Epigenomics data. RESULTS: As we tested the associations between 439 289 SNPs and bone erosion in 385 patients with erosive RA and 326 with non-erosive RA, none of the tested SNPs reached the genome-wide significance threshold, although many loci showed modest genetic effect on bone erosion status with suggestive association (e.g., rs2741200 [P = 3.75 × 10-6 ] in the SLA-TG locus and rs12422918 [P = 4.13 × 10-6 ] in SRGAP1). However, the top-ranked SNPs and their linked proxies, which were mostly located in non-coding variants, were significantly co-localized with the highly tissue-specific regulatory marker H3K4me3 in CD8+ memory T-cells (P = 0.014). CONCLUSION: Although, there was no large-effect variants associated with bone erosion in our GWAS, we have shown that CD8+ memory T-cells may have relevance with bone erosion in patients with RA through the analysis of ChiP-seq data.

OBJECTIVE: Genetic factors underlying susceptibility to rheumatoid arthritis (RA) in Arab populations are largely unknown. This genome-wide association study (GWAS) was undertaken to explore the generalizability of previously reported RA loci to Arab subjects and to discover new Arab-specific genetic loci. METHODS: The Genetics of Rheumatoid Arthritis in Some Arab States Study was designed to examine the genetics and clinical features of RA patients from Jordan, the Kingdom of Saudi Arabia, Lebanon, Qatar, and the United Arab Emirates. In total, >7 million single-nucleotide polymorphisms (SNPs) were tested for association with RA overall and with seropositive or seronegative RA in 511 RA cases and 352 healthy controls. In addition, replication of 15 signals was attempted in 283 RA cases and 221 healthy controls. A genetic risk score of 68 known RA SNPs was also examined in this study population. RESULTS: Three loci (HLA region, intergenic 5q13, and 17p13 at SMTNL2/GGT6) reached genome-wide significance in the analyses of association with RA and with seropositive RA, and for all 3 loci, evidence of independent replication was demonstrated. Consistent with the findings in European and East Asian populations, the association of RA with HLA-DRB1 amino acid position 11 conferred the strongest effect (P = 4.8 × 10-16 ), and a weighted genetic risk score of previously associated RA loci was found to be associated with RA (P = 3.41 × 10-5 ) and with seropositive RA (P = 1.48 × 10-6 ) in this population. In addition, 2 novel associations specific to Arab populations were found at the 5q13 and 17p13 loci. CONCLUSION: This first RA GWAS in Arab populations confirms that established HLA-region and known RA risk alleles contribute strongly to the risk and severity of disease in some Arab groups, suggesting that the genetic architecture of RA is similar across ethnic groups. Moreover, this study identified 2 novel RA risk loci in Arabs, offering further population-specific insights into the pathophysiology of RA.

OBJECTIVE: Rheumatoid factor (RF) is a well-established diagnostic and prognostic biomarker in rheumatoid arthritis (RA). However, ∼20% of RA patients are negative for this anti-IgG antibody. To date, only variation at the HLA-DRB1 gene has been associated with the presence of RF. This study was undertaken to identify additional genetic variants associated with RF positivity. METHODS: A genome-wide association study (GWAS) for RF positivity was performed using an Illumina Quad610 genotyping platform. A total of 937 RF-positive and 323 RF-negative RA patients were genotyped for >550,000 single-nucleotide polymorphisms (SNPs). Association testing was performed using an allelic chi-square test implemented in Plink software. An independent cohort of 472 RF-positive and 190 RF-negative RA patients was used to validate the most significant findings. RESULTS: In the discovery stage, a SNP in the IRX1 locus on chromosome 5p15.3 (SNP rs1502644) showed a genome-wide significant association with RF positivity (P = 4.13 × 10(-8) , odds ratio [OR] 0.37 [95% confidence interval (95% CI) 0.26-0.53]). In the validation stage, the association of IRX1 with RF was replicated in an independent group of RA patients (P = 0.034, OR 0.58 [95% CI 0.35-0.97] and combined P = 1.14 × 10(-8) , OR 0.43 [95% CI 0.32-0.58]). CONCLUSION: To our knowledge, this is the first GWAS of RF positivity in RA. Variation at the IRX1 locus on chromosome 5p15.3 is associated with the presence of RF. Our findings indicate that IRX1 and HLA-DRB1 are the strongest genetic factors for RF production in RA.

The aim of this study was to identify genetic variants associated with rheumatoid arthritis (RA) risk in black South Africans. Black South African RA patients (n = 263) were compared with healthy controls (n = 374). Genotyping was performed using the Immunochip, and four-digit high-resolution human leukocyte antigen (HLA) typing was performed by DNA sequencing of exon 2. Standard quality control measures were implemented on the data. The strongest associations were in the intergenic region between the HLA-DRB1 and HLA-DQA1 loci. After conditioning on HLA-DRB1 alleles, the effect in the rest of the extended major histocompatibility (MHC) diminished. Non-HLA single nucleotide polymorphisms (SNPs) in the intergenic regions LOC389203|RBPJ, LOC100131131|IL1R1, KIAA1919|REV3L, LOC643749|TRAF3IP2, and SNPs in the intron and untranslated regions (UTR) of IRF1 and the intronic region of ICOS and KIAA1542 showed association with RA (p < 5 × 10(-5)). Of the SNPs previously associated with RA in Caucasians, one SNP, rs874040, locating to the intergenic region LOC389203|RBPJ was replicated in this study. None of the variants in the PTPN22 gene was significantly associated. The seropositive subgroups showed similar results to the overall cohort. The effects observed across the HLA region are most likely due to HLA-DRB1, and secondary effects in the extended MHC cannot be detected. Seven non-HLA loci are associated with RA in black South Africans. Similar to Caucasians, the intergenic region between LOC38920 and RBPJ is associated with RA in this population. The strong association of the R620W variant of the PTPN22 gene with RA in Caucasians was not replicated since this variant was monomorphic in our study, but other SNP variants of the PTPN22 gene were also not associated with RA in black South Africans, suggesting that this locus does not play a major role in RA in this population.

OBJECTIVE: Methotrexate (MTX) is the drug of first choice for the treatment of rheumatoid arthritis (RA), but is effective only in around 60% of the patients. Identification of genetic markers to predict response is essential for effective treatment within a critical window period of 6 months after diagnosis, but have been hitherto elusive. In this study, we used genome-wide genotype data to identify the potential risk variants associated with MTX (poor)response in a north Indian RA cohort. MATERIALS AND METHODS: Genome-wide genotyping data for a total of 457 RA patients [297 good (DAS28-3≤3.2) and 160 poor (DAS28-3≥5.1) responders] on MTX monotherapy were tested for association using an additive model. Support vector machine and genome-wide pathway analysis were used to identify additional risk variants and pathways. All risk loci were imputed to fine-map the association signals and identify causal variant(s) of therapeutic/diagnostic relevance. RESULTS: Seven novel suggestive loci from genome-wide (P≤5×10(-5)) and three from support vector machine analysis were associated with MTX (poor)response. The associations of published candidate genes namely DHFR (P=0.014), FPGS (P=0.035), and TYMS (P=0.005) and purine and nucleotide metabolism pathways were reconfirmed. Imputation, followed by bioinformatic analysis indicated possible interaction between two reversely oriented overlapping genes namely ENOSF1 and TYMS at the post-transcriptional level. CONCLUSION: In this first ever genome-wide analysis on MTX treatment response in RA patients, 10 new risk loci were identified. These preliminary findings warrant replication in independent studies. Further, TYMS expression at the post-transcriptional level seems to be probably regulated through an antisense-RNA involving the 6-bp ins/del marker in the overlapping segment at 3'UTR of TYMS-ENOSF1, a finding with impending pharmacogenetic applications.

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.

OBJECTIVE: To investigate differences in genetic risk factors for rheumatoid arthritis (RA) in Han Chinese as compared with Europeans. METHODS: A genome-wide association study was conducted in China with 952 patients and 943 controls, and 32 variants were followed up in 2,132 patients and 2,553 controls. A transpopulation meta-analysis with results from a large European RA study was also performed to compare the genetic architecture across the 2 ethnic remote populations. RESULTS: Three non-major histocompatibility complex (non-MHC) loci were identified at the genome-wide significance level, the effect sizes of which were larger in anti-citrullinated protein antibody (ACPA)-positive patients than in ACPA-negative patients. These included 2 novel variants, rs12617656, located in an intron of DPP4 (odds ratio [OR] 1.56, P = 1.6 × 10(-21) ), and rs12379034, located in the coding region of CDK5RAP2 (OR 1.49, P = 1.1 × 10(-16) ), as well as a variant at the known CCR6 locus, rs1854853 (OR 0.71, P = 6.5 × 10(-15) ). The analysis of ACPA-positive patients versus ACPA-negative patients revealed that rs12617656 at the DPP4 locus showed a strong interaction effect with ACPAs (P = 5.3 × 10(-18) ), and such an interaction was also observed for rs7748270 at the MHC locus (P = 5.9 × 10(-8) ). The transpopulation meta-analysis showed genome-wide overlap and enrichment in association signals across the 2 populations, as confirmed by prediction analysis. CONCLUSION: This study has expanded the list of alleles that confer risk of RA, provided new insight into the pathogenesis of RA, and added empirical evidence to the emerging polygenic nature of complex trait variation driven by common genetic variants.

Celiac disease is a common autoimmune disorder characterized by an intestinal inflammation triggered by gluten, a storage protein found in wheat, rye and barley. Similar to other autoimmune diseases such as type 1 diabetes, psoriasis and rheumatoid arthritis, celiac disease is the result of an immune response to self-antigens leading to tissue destruction and production of autoantibodies. Common diseases like celiac disease have a complex pattern of inheritance with inputs from both environmental as well as additive and non-additive genetic factors. In the past few years, Genome Wide Association Studies (GWAS) have been successful in finding genetic risk variants behind many common diseases and traits. To complement and add to the previous findings, we performed a GWAS including 206 trios from 97 nuclear Swedish and Norwegian families affected with celiac disease. By stratifying for HLA-DQ, we identified a new genome-wide significant risk locus covering the DUSP10 gene. To further investigate the associations from the GWAS we performed pathway analyses and two-locus interaction analyses. These analyses showed an over-representation of genes involved in type 2 diabetes and identified a set of candidate mechanisms and genes of which some were selected for mRNA expression analysis using small intestinal biopsies from 98 patients. Several genes were expressed differently in the small intestinal mucosa from patients with celiac autoimmunity compared to intestinal mucosa from control patients. From top-scoring regions we identified susceptibility genes in several categories: 1) polarity and epithelial cell functionality; 2) intestinal smooth muscle; 3) growth and energy homeostasis, including proline and glutamine metabolism; and finally 4) innate and adaptive immune system. These genes and pathways, including specific functions of DUSP10, together reveal a new potential biological mechanism that could influence the genesis of celiac disease, and possibly also other chronic disorders with an inflammatory component.

OBJECTIVE: The number of confirmed rheumatoid arthritis (RA) loci currently stands at 32, but many lines of evidence indicate that expansion of existing genome-wide association studies (GWAS) enhances the power to detect additional loci. This study was undertaken to extend our previous RA GWAS in a UK cohort, adding more independent RA cases and healthy controls, with the aim of detecting novel association signals for susceptibility to RA in a homogeneous UK cohort. METHODS: A total of 3,223 UK RA cases and 5,272 UK controls were available for association analyses, with the extension adding 1,361 cases and 2,334 controls to the original GWAS data set. The genotype data for all RA cases were imputed using the Impute program version 2. After stringent quality control thresholds were applied, 3,034 cases and 5,271 controls (1,831,729 single-nucleotide polymorphisms [SNPs]) were available for analysis. Association testing was performed using Plink software. RESULTS: The analyses indicated a suggestive association with susceptibility to RA (P < 0.0001) for 6 novel RA loci that have been previously found to be associated with other autoimmune diseases; these 6 SNPs were genotyped in independent samples. Two of the associated loci were validated, one of which was associated with RA at genome-wide levels of significance in the combined analysis, identifying a novel RA locus at 22q12 (P = 6.9 × 10(-9) ). In addition, most of the previously known RA susceptibility loci were confirmed to be associated with RA, and for 16 of the loci, the strength of the association was increased. CONCLUSION: This study identified a new RA locus mapping to 22q12. These results support the notion that increasing the power of GWAS enhances novel gene discovery.

Rheumatoid arthritis (RA) is a chronic inflammatory disorder with a polygenic mode of inheritance. This study examined the hypothesis that runs of homozygosity (ROHs) play a recessive-acting role in the underlying RA genetic mechanism and identified RA-associated ROHs. Ours is the first genome-wide homozygosity association study for RA and characterized the ROH patterns associated with RA in the genomes of 2,000 RA patients and 3,000 normal controls of the Wellcome Trust Case Control Consortium. Genome scans consistently pinpointed two regions within the human major histocompatibility complex region containing RA-associated ROHs. The first region is from 32,451,664 bp to 32,846,093 bp (-log10(p)>22.6591). RA-susceptibility genes, such as HLA-DRB1, are contained in this region. The second region ranges from 32,933,485 bp to 33,585,118 bp (-log10(p)>8.3644) and contains other HLA-DPA1 and HLA-DPB1 genes. These two regions are physically close but are located in different blocks of linkage disequilibrium, and ∼40% of the RA patients' genomes carry these ROHs in the two regions. By analyzing homozygote intensities, an ROH that is anchored by the single nucleotide polymorphism rs2027852 and flanked by HLA-DRB6 and HLA-DRB1 was found associated with increased risk for RA. The presence of this risky ROH provides a 62% accuracy to predict RA disease status. An independent genomic dataset from 868 RA patients and 1,194 control subjects of the North American Rheumatoid Arthritis Consortium successfully validated the results obtained using the Wellcome Trust Case Control Consortium data. In conclusion, this genome-wide homozygosity association study provides an alternative to allelic association mapping for the identification of recessive variants responsible for RA. The identified RA-associated ROHs uncover recessive components and missing heritability associated with RA and other autoimmune diseases.