Leonid Padyukov

Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci

To identify new genetic risk factors for rheumatoid arthritis, we conducted a genome-wide association study meta-analysis of 5,539 autoantibody-positive individuals with rheumatoid arthritis (cases) and 20,169 controls of European descent, followed by replication in an independent set of 6,768 rheumatoid arthritis cases and 8,806 controls. Of 34 SNPs selected for replication, 7 new rheumatoid arthritis risk alleles were identified at genome-wide significance (P < 5 × 10−8) in an analysis of all 41,282 samples. The associated SNPs are near genes of known immune function, including IL6ST, SPRED2, RBPJ, CCR6, IRF5 and PXK. We also refined associations at two established rheumatoid arthritis risk loci (IL2RA and CCL21) and confirmed the association at AFF3. These new associations bring the total number of confirmed rheumatoid arthritis risk loci to 31 among individuals of European ancestry. An additional 11 SNPs replicated at P < 0.05, many of which are validated autoimmune risk alleles, suggesting that most represent genuine rheumatoid arthritis risk alleles.

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, 3, 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.

A large-scale replication study identifies TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 as risk loci for systemic lupus erythematosus

Genome-wide association studies have recently identified at least 15 susceptibility loci for systemic lupus erythematosus (SLE). To confirm additional risk loci, we selected SNPs from 2,466 regions that showed nominal evidence of association to SLE (P < 0.05) in a genome-wide study and genotyped them in an independent sample of 1,963 cases and 4,329 controls. This replication effort identified five new SLE susceptibility loci (P < 5 × 10−8): TNIP1 (odds ratio (OR) = 1.27), PRDM1 (OR = 1.20), JAZF1 (OR = 1.20), UHRF1BP1 (OR = 1.17) and IL10 (OR = 1.19). We identified 21 additional candidate loci with P≤ 1 × 10−5. A candidate screen of alleles previously associated with other autoimmune diseases suggested five loci (P < 1 × 10−3) that may contribute to SLE: IFIH1, CFB, CLEC16A, IL12B and SH2B3. These results expand the number of confirmed and candidate SLE susceptibility loci and implicate several key immunologic pathways in SLE pathogenesis.

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, ∼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 × 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.

Replication of Putative Candidate-Gene Associations with Rheumatoid Arthritis in >4,000 Samples from North America and Sweden: Association of Susceptibility with PTPN22, CTLA4, and PADI4

Candidate-gene association studies in rheumatoid arthritis (RA) have lead to encouraging yet apparently inconsistent results. One explanation for the inconsistency is insufficient power to detect modest effects in the context of a low prior probability of a true effect. To overcome this limitation, we selected alleles with an increased probability of a disease association, on the basis of a review of the literature on RA and other autoimmune diseases, and tested them for association with RA susceptibility in a sample collection powered to detect modest genetic effects. We tested 17 alleles from 14 genes in 2,370 RA cases and 1,757 controls from the North American Rheumatoid Arthritis Consortium (NARAC) and the Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA) collections. We found strong evidence of an association of PTPN22 with the development of anti-citrulline antibody-positive RA (odds ratio [OR] 1.49; P=.00002), using previously untested EIRA samples. We provide support for an association of CTLA4 (CT60 allele, OR 1.23; P=.001) and PADI4 (PADI4_94, OR 1.24; P=.001) with the development of RA, but only in the NARAC cohort. The CTLA4 association is stronger in patients with RA from both cohorts who are seropositive for anti-citrulline antibodies (P=.0006). Exploration of our data set with clinically relevant subsets of RA reveals that PTPN22 is associated with an earlier age at disease onset (P=.004) and that PTPN22 has a stronger effect in males than in females (P=.03). A meta-analysis failed to demonstrate an association of the remaining alleles with RA susceptibility, suggesting that the previously published associations may represent false-positive results. Given the strong statistical power to replicate a true-positive association in this study, our results provide support for PTPN22, CTLA4, and PADI4 as RA susceptibility genes and demonstrate novel associations with clinically relevant subsets of RA.

Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis

Epigenetic mechanisms integrate genetic and environmental causes of disease, but comprehensive genome-wide analyses of epigenetic modifications have not yet demonstrated robust association with common diseases. Using Illumina HumanMethylation450 arrays on 354 anti-citrullinated protein antibody–associated rheumatoid arthritis cases and 337 controls, we identified two clusters within the major histocompatibility complex (MHC) region whose differential methylation potentially mediates genetic risk for rheumatoid arthritis. To reduce confounding factors that have hampered previous epigenome-wide studies, we corrected for cellular heterogeneity by estimating and adjusting for cell-type proportions in our blood-derived DNA samples and used mediation analysis to filter out associations likely to be a consequence of disease. Four CpGs also showed an association between genotype and variance of methylation. The associations for both clusters replicated at least one CpG (P < 0.01), with the rest showing suggestive association, in monocyte cell fractions in an independent cohort of 12 cases and 12 controls. Thus, DNA methylation is a potential mediator of genetic risk.

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

Genome-wide association identifies multiple ulcerative colitis susceptibility loci

Ulcerative colitis is a chronic, relapsing inflammatory condition of the gastrointestinal tract with a complex genetic and environmental etiology. In an effort to identify genetic variation underlying ulcerative colitis risk, we present two distinct genome-wide association studies of ulcerative colitis and their joint analysis with a previously published scan, comprising, in aggregate, 2,693 individuals with ulcerative colitis and 6,791 control subjects. Fifty-nine SNPs from 14 independent loci attained an association significance of P < 10−5. Seven of these loci exceeded genome-wide significance (P < 5 × 10−8). After testing an independent cohort of 2,009 cases of ulcerative colitis and 1,580 controls, we identified 13 loci that were significantly associated with ulcerative colitis (P < 5 × 10−8), including the immunoglobulin receptor gene FCGR2A, 5p15, 2p16 and ORMDL3 (orosomucoid1-like 3). We confirmed association with 14 previously identified ulcerative colitis susceptibility loci, and an analysis of acknowledged Crohns disease loci showed that roughly half of the known Crohns disease associations are shared with ulcerative colitis. These data implicate approximately 30 loci in ulcerative colitis, thereby providing insight into disease pathogenesis.

Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis

The genetic association of the major histocompatibility complex (MHC) to rheumatoid arthritis risk has commonly been attributed to alleles in HLA-DRB1. However, debate persists about the identity of the causal variants in HLA-DRB1 and the presence of independent effects elsewhere in the MHC. Using existing genome-wide SNP data in 5,018 individuals with seropositive rheumatoid arthritis (cases) and 14,974 unaffected controls, we imputed and tested classical alleles and amino acid polymorphisms in HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1 and HLA-DRB1, as well as 3,117 SNPs across the MHC. Conditional and haplotype analyses identified that three amino acid positions (11, 71 and 74) in HLA-DRβ1 and single–amino-acid polymorphisms in HLA-B (at position 9) and HLA-DPβ1 (at position 9), which are all located in peptide-binding grooves, almost completely explain the MHC association to rheumatoid arthritis risk. This study shows how imputation of functional variation from large reference panels can help fine map association signals in the MHC.

Immunity to citrullinated proteins in rheumatoid arthritis.

Antibodies to citrullinated proteins (ACPA), i.e., to peptides posttranslationally modified by the conversion of arginine to citrulline, are specific serological markers for rheumatoid arthritis (RA). Studies on anticitrulline immunity, summarized in this review, demonstrate that the criterion-based syndrome RA should be subdivided into at least two distinct subsets (ACPA-positive and ACPA-negative disease). A new etiological model is proposed for ACPA-positive RA, built on MHC class II-dependent activation of adaptive immunity. Fundamentals of this model include the following: (a) ACPA antedate onset of arthritis; (b) ACPA may aggravate arthritis in rodents; (c) ACPA are triggered in the context of genes that confer susceptibility to RA (HLA-DRB1 SE) and by environmental agents triggering RA (smoking or bacterial stimuli); (d) ACPA may complex with citrullinated proteins present in target tissue as part of a multistep process for arthritis development. The model provides a new basis for molecular studies on the pathogenesis of ACPA-positive arthritis.