Rheumatoid Arthritis: Plenge RM

Plenge RM. · Division of Rheumatology, Immunology and Allergy, Brigham and Women's, Hospital, Harvard Medical School, Boston, USA. · Curr Opin Rheumatol. · Pubmed #19365266 No free full text.

Abstract: PURPOSE OF REVIEW: Recent human genetic discoveries have increased our understanding of rheumatoid arthritis (RA) susceptibility. These discoveries are reviewed and placed in the context of potential important clinical applications. RECENT FINDINGS: Genome-wide association studies and related methodologies have expanded the number of validated RA risk loci beyond HLA-DRB1 'shared epitope' alleles to include additional major histocompatibility complex risk alleles and more than 10 regions outside the major histocompatibility complex locus. The newly discovered risk alleles are common in the general population, and most have a modest effect on risk of RA (odds ratio approximately 1.15 per copy of each risk allele). Although the actual causal mutation and causal gene for most loci remain to be determined, these studies are beginning to reveal general themes: many risk loci are associated with other autoimmune diseases, many genes fall within discrete biological pathways (e.g., the NF-kappaB signaling pathway), and human genetics can subset disease into clinically meaningful categories (e.g., presence or absence of autoantibodies). SUMMARY: Approximately one-third of the genetic basis of RA can be explained by known risk loci. Future studies need to pinpoint the actual causal mutations, expand the number of risk loci, and translate these discoveries to improve care of patients with RA.

Fernando MM, Stevens CR, Walsh EC, De Jager PL, Goyette P, Plenge RM, Vyse TJ, Rioux JD. · Section of Molecular Genetics and Rheumatology, Faculty of Medicine, Imperial College London, London, United Kingdom. · PLoS Genet. · Pubmed #18437207 links to  free full text

Abstract: The major histocompatibility complex (MHC) is one of the most extensively studied regions in the human genome because of the association of variants at this locus with autoimmune, infectious, and inflammatory diseases. However, identification of causal variants within the MHC for the majority of these diseases has remained difficult due to the great variability and extensive linkage disequilibrium (LD) that exists among alleles throughout this locus, coupled with inadequate study design whereby only a limited subset of about 20 from a total of approximately 250 genes have been studied in small cohorts of predominantly European origin. We have performed a review and pooled analysis of the past 30 years of research on the role of the MHC in six genetically complex disease traits - multiple sclerosis (MS), type 1 diabetes (T1D), systemic lupus erythematosus (SLE), ulcerative colitis (UC), Crohn's disease (CD), and rheumatoid arthritis (RA) - in order to consolidate and evaluate the current literature regarding MHC genetics in these common autoimmune and inflammatory diseases. We corroborate established MHC disease associations and identify predisposing variants that previously have not been appreciated. Furthermore, we find a number of interesting commonalities and differences across diseases that implicate both general and disease-specific pathogenetic mechanisms in autoimmunity.

Plenge RM, Criswell LA. · Division of Rheumatology, Immunology and Allergy, Harvard Medical School, Brigham and Women's Hospital, Boston, MA, USA. · Curr Opin Rheumatol. · Pubmed #18349743 No free full text.

Abstract: PURPOSE OF REVIEW: Tumor necrosis factor alpha (TNF) inhibitors are a mainstay of treatment in rheumatoid arthritis, yet there are no effective clinical or biomarker predictors of which patients will respond. Here we review genetic association studies conducted to search for DNA biomarkers of response to anti-TNF therapy. RECENT FINDINGS: The entirety of genetic association studies to date that focus on response to anti-TNF therapy has been limited to a small number of genetic variants within a few candidate genes (primarily within the major histocompatibility complex region). Moreover, these studies have been conducted in a relatively small number of rheumatoid arthritis patients (approximately 1000 patients across all studies combined). From these studies, no single genetic factor is associated unequivocally with treatment response, although some studies suggest that alleles within the major histocompatibility complex may influence response. SUMMARY: Additional studies are required to investigate the genetic basis of response to anti-TNF therapy. These studies should include an unbiased search of DNA variation across the human genome--now feasible through cost-effective genome-wide association studies--and be conducted in large patient collections powered to detect modest effect sizes.

Gregersen PK, Amos CI, Lee AT, Lu Y, Remmers EF, Kastner DL, Seldin MF, Criswell LA, Plenge RM, Holers VM, Mikuls TR, Sokka T, Moreland LW, Bridges SL, Xie G, Begovich AB, Siminovitch KA. · The Feinstein Institute for Medical Research, North Shore-Long Island Jewish Health System, Manhasset, NY, USA. · Nat Genet. · Pubmed #19503088 No free full text.

Abstract: We conducted a genome-wide association study of rheumatoid arthritis in 2,418 cases and 4,504 controls from North America and identified an association at the REL locus, encoding c-Rel, on chromosome 2p13 (rs13031237, P = 6.01 x 10(-10)). Replication in independent case-control datasets comprising 2,604 cases and 2,882 controls confirmed this association, yielding an allelic OR = 1.25 (P = 3.08 x 10(-14)) for marker rs13031237 and an allelic OR = 1.21 (P = 2.60 x 10(-11)) for marker rs13017599 in the combined dataset. The combined dataset also provides definitive support for associations at both CTLA4 (rs231735; OR = 0.85; P = 6.25 x 10(-9)) and BLK (rs2736340; OR = 1.19; P = 5.69 x 10(-9)). c-Rel is an NF-kappaB family member with distinct functional properties in hematopoietic cells, and its association with rheumatoid arthritis suggests disease pathways that involve other recently identified rheumatoid arthritis susceptibility genes including CD40, TRAF1, TNFAIP3 and PRKCQ.

Lee YC, Raychaudhuri S, Cui J, De Vivo I, Ding B, Alfredsson L, Padyukov L, Costenbader KH, Seielstad M, Graham RR, Klareskog L, Gregersen PK, Plenge RM, Karlson EW. · Brigham and Women's Hospital, Boston, Massachusetts 02115, USA. · Arthritis Rheum. · Pubmed #19404952 No free full text.

Abstract: OBJECTIVE: Previous studies have demonstrated that the PRL -1149 T (minor) allele decreases prolactin expression and may be associated with autoimmune disease. The aim of this study was to determine the role of the PRL -1149 G/T polymorphism (rs1341239) in rheumatoid arthritis (RA) susceptibility. METHODS: We examined the association between PRL -1149 G/T and RA risk in 4 separate study populations, consisting of a total of 3,405 RA cases and 4,111 controls of self-reported white European ancestry. Samples were genotyped using 1 of 3 genotyping platforms, and strict quality control metrics were applied. We tested for association using a 2-tailed Cochran-Mantel-Haenszel additive, fixed-effects model. RESULTS: In the individual populations, odds ratios (ORs) for an association between PRL -1149 T and RA risk ranged from 0.80 to 0.97. In a joint meta-analysis across all 4 populations, the OR for an association between PRL -1149 T and RA risk was 0.90 (95% confidence interval 0.84-0.96, P=0.001). CONCLUSION: Our findings indicate a possible association between the PRL -1149 T allele and decreased RA risk. The effect size is small but similar to ORs for other genetic polymorphisms associated with complex traits, including RA.

Cui J, Taylor KE, Destefano AL, Criswell LA, Izmailova ES, Parker A, Roubenoff R, Plenge RM, Weinblatt ME, Shadick NA, Karlson EW. · Division of Rheumatology, Immunology & Allergy, Brigham and Women's Hospital and Harvard Medical School, Boston, Massachusetts 02115, USA. · Mol Med. · Pubmed #19287509 links to  free full text

Abstract: We carried out a genome-wide association study of genetic predictors of anti-cyclic citrullinated peptide antibody (anti-CCP) level in 531 self-reported non-Hispanic Caucasian Rheumatoid Arthritis (RA) patients enrolled in the Brigham Rheumatoid Arthritis Sequential Study (BRASS). For replication, we then analyzed 289 single nucleotide polymorphisms (SNPs) with P < 0.001 in BRASS in an independent population of 849 RA patients from the North American Rheumatoid Arthritis Consortium (NARAC). BRASS and NARAC samples were genotyped using the Affymetrix 100K and Illumina 550K platforms respectively. Association between SNPs and anti-CCP titer was tested using general linear models. The five most significant SNPs from BRASS all were within the major histocompatibility complex (MHC) region (P < or = 3.5 x 10(-6)). After controlling for the human leukocyte antigen shared epitope (HLA-SE), the top SNPs still yielded P values < 0.0002. In NARAC, a single SNP from the MHC region near BTNL2 and HLA-DRA, rs1980493 (r(2) = 0.85 with the top five SNPs from BRASS), was associated significantly with CCP titer (P = 6.1 x 10(-5)) even after adjustment for the HLA-SE (P = 0.0002). The top SNPs found in BRASS and NARAC had r(2) = 0.46 and 0.64, respectively, to HLA-DRB1 DR3 alleles. These results confirm that the most significant genome region affecting anti-CCP titers in RA is the MHC region. We identified a SNP in moderate linkage disequilibrium (LD) with HLA-DR3, which may influence anti-CCP titer independently of the HLA-SE.

Liao KP, Gunnarsson M, Källberg H, Ding B, Plenge RM, Padyukov L, Karlson EW, Klareskog L, Askling J, Alfredsson L. · Division of Rheumatology, Immunology, and Allergy, Brigham and Women's Hospital, Boston, MA 02114, USA. · Arthritis Rheum. · Pubmed #19248096 No free full text.

Abstract: OBJECTIVE: The co-occurrence of autoimmune diseases such as rheumatoid arthritis (RA) and type 1 diabetes mellitus (DM) has been reported in individuals and families. In this study, the strength and nature of this association were investigated at the population level in a Swedish case-control cohort. METHODS: For this case-control study, 1,419 patients with incident RA diagnosed between 1996 and 2003 were recruited from university, public, and private rheumatology units throughout Sweden; 1,674 matched control subjects were recruited from the Swedish national population registry. Sera from the subjects were tested for the presence of antibodies to cyclic citrullinated peptide (anti-CCP), rheumatoid factor (RF), and the 620W PTPN22 allele. Information on a history of diabetes was obtained by questionnaire, telephone interview, and/or medical record review. The prevalence of type 1 DM and type 2 DM was compared between patients with incident RA and control subjects and further stratified for the presence of anti-CCP, RF, and the PTPN22 risk allele. RESULTS: Type 1 DM was associated with an increased risk of RA (odds ratio [OR] 4.9, 95% confidence interval [95% CI] 1.8-13.1), and this association was specific for anti-CCP-positive RA (OR 7.3, 95% CI 2.7-20.0), but not anti-CCP-negative RA. Further adjustment for the presence of PTPN22 attenuated the risk of anti-CCP-positive RA in patients with type 1 DM to an OR of 5.3 (95% CI 1.5-18.7). No association between RA and type 2 DM was observed. CONCLUSION: The association between type 1 DM and RA is specific for a particular RA subset, anti-CCP-positive RA. The risk of developing RA later in life in patients with type 1 DM may be attributed, in part, to the presence of the 620W PTPN22 allele, suggesting that this risk factor may represent a common pathway for the pathogenesis of these 2 diseases.

Graham RR, Cotsapas C, Davies L, Hackett R, Lessard CJ, Leon JM, Burtt NP, Guiducci C, Parkin M, Gates C, Plenge RM, Behrens TW, Wither JE, Rioux JD, Fortin PR, Graham DC, Wong AK, Vyse TJ, Daly MJ, Altshuler D, Moser KL, Gaffney PM. · Program in Medical and Population Genetics, Broad Institute of Harvard and Massachusetts Institute of Technology, Cambridge, MA 02142, USA. · Nat Genet. · Pubmed #19165918 No free full text.

Abstract: Systemic lupus erythematosus (SLE) is an autoimmune disease influenced by genetic and environmental factors. We carried out a genome-wide association scan and replication study and found an association between SLE and a variant in TNFAIP3 (rs5029939, meta-analysis P = 2.89 x 10(-12), OR = 2.29). We also found evidence of two independent signals near TNFAIP3 associated with SLE, including one previously associated with rheumatoid arthritis (RA). These results establish that variants near TNFAIP3 contribute to differential risk of SLE and RA.

Korman BD, Seldin MF, Taylor KE, Le JM, Lee AT, Plenge RM, Amos CI, Criswell LA, Gregersen PK, Kastner DL, Remmers EF. · National Institute of Arthritis and Musculoskeletal and Skin Diseases, NIH, Bethesda, Maryland 20892-1849, USA. · Arthritis Rheum. · Pubmed #19116934 No free full text.

Abstract: OBJECTIVE: The single-nucleotide polymorphism (SNP) rs11761231 on chromosome 7q has been reported to be sexually dimorphic marker for rheumatoid arthritis (RA) susceptibility in a British population. We sought to replicate this finding and to better characterize susceptibility alleles in the region in a North American population. METHODS: DNA from 2 North American collections of RA patients and controls (1,605 cases and 2,640 controls) was genotyped for rs11761231 and 16 additional chromosome 7q tag SNPs using Sequenom iPlex assays. Association tests were performed for each collection and also separately, contrasting male cases with male controls and female cases with female controls. Principal components analysis (EigenStrat) was used to determine association with RA before and after adjusting for population stratification in the subset of the samples for which there were whole-genome SNP data (772 cases and 1,213 controls). RESULTS: We failed to replicate an association of the 7q region with RA. Initially, rs11761231 showed evidence for association with RA in the North American Rheumatoid Arthritis Consortium (NARAC) collection (P=0.0073), and rs11765576 showed association with RA in both the NARAC (P=0.038) and RA replication (P = 0.0013) collections. These markers also exhibited sex differentiation. However, in the whole-genome subset, neither SNP showed significant association with RA after correction for population stratification. CONCLUSION: While 2 SNPs on chromosome 7q appeared to be associated with RA in a North American cohort, the significance of this finding did not withstand correction for population substructure. Our results emphasize the need to carefully account for population structure to avoid false-positive disease associations.

Ding B, Padyukov L, Lundström E, Seielstad M, Plenge RM, Oksenberg JR, Gregersen PK, Alfredsson L, Klareskog L. · Karolinska Institutet, Stockholm, Sweden. · Arthritis Rheum. · Pubmed #19116921 No free full text.

Abstract: OBJECTIVE: To identify additional variants in the major histocompatibility complex (MHC) region that independently contribute to risk in 2 disease subsets of rheumatoid arthritis (RA) defined according to the presence or absence of antibodies to citrullinated protein antigens (ACPAs). METHODS: In a multistep analytical strategy using unmatched as well as matched analyses to adjust for HLA-DRB1 genotype, we analyzed 2,221 single-nucleotide polymorphisms (SNPs) spanning 10.7 Mb, from 6p22.2 to 6p21.31, across the MHC. For ACPA-positive RA, we analyzed samples from the Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA) and the North American Rheumatoid Arthritis Consortium (NARAC) studies (totaling 1,255 cases and 1,719 controls). For ACPA-negative RA, we used samples from the EIRA study (640 cases and 670 controls). Plink and SAS statistical packages were used to conduct all statistical analyses. RESULTS: A total of 299 SNPs reached locus-wide significance (P<2.3x10(-5)) for ACPA-positive RA, whereas surprisingly, no SNPs reached this significance for ACPA-negative RA. For ACPA-positive RA, we adjusted for known DRB1 risk alleles and identified additional independent associations with SNPs near HLA-DPB1 (rs3117213; odds ratio 1.42 [95% confidence interval 1.17-1.73], Pcombined=0.0003 for the strongest association). CONCLUSION: There are distinct genetic patterns of MHC associations in the 2 disease subsets of RA defined according to ACPA status. HLA-DPB1 is an independent risk locus for ACPA-positive RA. We did not identify any associations with SNPs within the MHC for ACPA-negative RA.

Raychaudhuri S, Remmers EF, Lee AT, Hackett R, Guiducci C, Burtt NP, Gianniny L, Korman BD, Padyukov L, Kurreeman FA, Chang M, Catanese JJ, Ding B, Wong S, van der Helm-van Mil AH, Neale BM, Coblyn J, Cui J, Tak PP, Wolbink GJ, Crusius JB, van der Horst-Bruinsma IE, Criswell LA, Amos CI, Seldin MF, Kastner DL, Ardlie KG, Alfredsson L, Costenbader KH, Altshuler D, Huizinga TW, Shadick NA, Weinblatt ME, de Vries N, Worthington J, Seielstad M, Toes RE, Karlson EW, Begovich AB, Klareskog L, Gregersen PK, Daly MJ, Plenge RM. · Program in Medical and Population Genetics, Broad Institute of Harvard and Massachusetts Institute of Technology, Cambridge, Massachusetts 02142, USA. · Nat Genet. · Pubmed #18794853 No free full text.

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

Liu C, Batliwalla F, Li W, Lee A, Roubenoff R, Beckman E, Khalili H, Damle A, Kern M, Furie R, Dupuis J, Plenge RM, Coenen MJ, Behrens TW, Carulli JP, Gregersen PK. · Biogen Idec Inc., Cambridge, Massachusetts, USA. · Mol Med. · Pubmed #18615156 links to  free full text

Abstract: The prediction of response (or non-response) to anti-TNF treatment for rheumatoid arthritis (RA) is a pressing clinical problem. We conducted a genome-wide association study using the Illumina HapMap300 SNP chip on 89 RA patients prospectively followed after beginning anti-TNF therapy as part of Autoimmune Biomarkers Collaborative Network (ABCoN [Autoimmune Bio-markers Collaborative Network]) patient cohort. Response to therapy was determined by the change in Disease Activity Score (DAS28) observed after 14 wks. We used a two-part analysis that treated the change in DAS28 as a continuous trait and then incorporated it into a dichotomous trait of "good responder" and "nonresponder" by European League Against Rheumatism (EULAR) criteria. We corrected for multiple tests by permutation, and adjusted for potential population stratification using EIGENSTRAT. Multiple single nucleotide polymorphism (SNP) markers showed significant associations near or within loci including: the v-maf musculoaponeurotic fibrosarcoma oncogene homolog B (MAFB) gene on chromosome 20; the type I interferon gene IFNk on chromosome 9; and in a locus on chromosome 7 that includes the paraoxonase I (PON1) gene. An SNP in the IL10 promoter (rs1800896) that was previously reported as associated with anti-TNF response was weakly associated with response in this cohort. Replications of these results in independent and larger data sets clearly are required. We provide a reference list of candidate SNPs (P < 0.01) that can be investigated in future pharmacogenomic studies.

Johnsen AK, Plenge RM, Butty V, Campbell C, Dieguez-Gonzalez R, Gomez-Reino JJ, Shadick N, Weinblatt M, Gonzalez A, Gregersen PK, Benoist C, Mathis D. · Joslin Diabetes Center, Harvard Medical School, Boston, Massachusetts, USA. · Arthritis Rheum. · Pubmed #18576312 links to  free full text

Abstract: OBJECTIVE: It has been suggested that polymorphisms in IL1 are correlated with severe and/or erosive rheumatoid arthritis (RA), but the implicated alleles have differed among studies. The aim of this study was to perform a broad and well-powered search for association between allelic polymorphism in IL1A and IL1B and the susceptibility to or severity of RA. METHODS: Key coding and regulatory regions in IL1A and IL1B were sequenced in 24 patients with RA, revealing 4 novel single-nucleotide polymorphisms (SNPs) in IL1B. These and a comprehensive set of 24 SNPs tagging most of the underlying genetic diversity were genotyped in 3 independent RA case-control sample sets and 1 longitudinal RA cohort, totaling 3,561 patients and 3,062 control subjects. RESULTS: No fully significant associations were observed. Analysis of the discovery case-control sample sets indicated a potential association of IL1B promoter region SNPs with susceptibility to RA (for RA3/A, odds ratio [OR] 1.27, P = 0.0021) or with the incidence of radiographic erosions (for RA4/C, OR 1.56, P = 0.036), but these findings were not replicated in independent case-control samples. No association with rheumatoid factor, anti-cyclic citrullinated peptide, or the Disease Activity Score in 28 joints was found. None of the associations previously observed in other studies were replicated here. CONCLUSION: In spite of a broad and highly powered study, we observed no robust, reproducible association between IL1A/B variants and the susceptibility to or severity of RA in white individuals of European descent. Our results provide evidence that, in the majority of cases, polymorphism in IL1A and IL1B is not a major contributor to genetic susceptibility to RA.

Amos CI, Chen WV, Remmers E, Siminovitch KA, Seldin MF, Criswell LA, Lee AT, John S, Shephard ND, Worthington J, Cornelis F, Plenge RM, Begovich AB, Dyer TD, Kastner DL, Gregersen PK. · Departments of Epidemiology and Biomathematics, University of Texas, M,D, Anderson Cancer Center, 1155 Pressler Street, Unit 1340, Houston, Texas 77030, USA. · BMC Proc. · Pubmed #18466527 links to  free full text

Abstract: ABSTRACT : For Genetic Analysis Workshop 15 Problem 2, we organized data from several ongoing studies designed to identify genetic and environmental risk factors for rheumatoid arthritis. Data were derived from the North American Rheumatoid Arthritis Consortium (NARAC), collaboration among Canadian researchers, the European Consortium on Rheumatoid Arthritis Families (ECRAF), and investigators from Manchester, England. All groups used a common standard for defining rheumatoid arthritis, but NARAC also further selected for a more severe phenotype in the probands. Genotyping and family structures for microsatellite-based linkage analysis were provided from all centers. In addition, all centers but ECRAF have genotyped families for linkage analysis using SNPs and these data were additionally provided. NARAC also had additional data from a dense genotyping analysis of a region of chromosome 18 and results from candidate gene studies, which were provided. Finally, smoking influences risk for rheumatoid arthritis, and data were provided from the NARAC study on this behavior as well as some additional phenotypes measuring severity. Several questions could be evaluated using the data that were provided. These include comparing linkage analysis using single-nucleotide polymorphisms versus microsatellites and identifying credible regions of linkage outside the HLA region on chromosome 6p13, which has been extensively documented; evaluating the joint effects of smoking with genetic factors; and identifying more homogenous subsets of families for whom genetic susceptibility might be stronger, so that linkage and association studies may be more efficiently conducted.

Lee HS, Lee AT, Criswell LA, Seldin MF, Amos CI, Carulli JP, Navarrete C, Remmers EF, Kastner DL, Plenge RM, Li W, Gregersen PK. · The Feinstein Institute for Medical Research, North Shore LIJ Health System, Manhasset, New York 11030, United States of America. · Mol Med. · Pubmed #18309376 links to  free full text

Abstract: Recent evidence suggests that additional risk loci for RA are present in the major histocompatibility complex (MHC), independent of the class II HLA-DRB1 locus. We have now tested a total of 1,769 SNPs across 7.5Mb of the MHC located from 6p22.2 (26.03 Mb) to 6p21.32 (33.59 Mb) derived from the Illumina 550K Beadchip (Illumina, San Diego, CA, USA). For an initial analysis in the whole dataset (869 RA CCP + cases, 1,193 controls), the strongest association signal was observed in markers near the HLA-DRB1 locus, with additional evidence for association extending out into the Class I HLA region. To avoid confounding that may arise due to linkage disequilibrium with DRB1 alleles, we analyzed a subset of the data by matching cases and controls by DRB1 genotype (both alleles matched 1:1), yielding a set of 372 cases with 372 controls. This analysis revealed the presence of at least two regions of association with RA in the Class I region, independent of DRB1 genotype. SNP alleles found on the conserved A1-B8-DR3 (8.1) haplotype show the strongest evidence of positive association (P ~ 0.00005) clustered in the region around the HLA-C locus. In addition, we identified risk alleles that are not present on the 8.1 haplotype, with maximal association signals (P ~ 0.001-0.0027) located near the ZNF311 locus. This latter association is enriched in DRB1*0404 individuals. Finally, several additional association signals were found in the extreme centromeric portion of the MHC, in regions containing the DOB1, TAP2, DPB1, and COL11A2 genes. These data emphasize that further analysis of the MHC is likely to reveal genetic risk factors for rheumatoid arthritis that are independent of the DRB1 shared epitope alleles.

Tian C, Plenge RM, Ransom M, Lee A, Villoslada P, Selmi C, Klareskog L, Pulver AE, Qi L, Gregersen PK, Seldin MF. · Rowe Program in Human Genetics, University of California Davis, Davis, California, United States of America. · PLoS Genet. · Pubmed #18208329 links to  free full text

Abstract: European population genetic substructure was examined in a diverse set of >1,000 individuals of European descent, each genotyped with >300 K SNPs. Both STRUCTURE and principal component analyses (PCA) showed the largest division/principal component (PC) differentiated northern from southern European ancestry. A second PC further separated Italian, Spanish, and Greek individuals from those of Ashkenazi Jewish ancestry as well as distinguishing among northern European populations. In separate analyses of northern European participants other substructure relationships were discerned showing a west to east gradient. Application of this substructure information was critical in examining a real dataset in whole genome association (WGA) analyses for rheumatoid arthritis in European Americans to reduce false positive signals. In addition, two sets of European substructure ancestry informative markers (ESAIMs) were identified that provide substantial substructure information. The results provide further insight into European population genetic substructure and show that this information can be used for improving error rates in association testing of candidate genes and in replication studies of WGA scans.

Plenge RM, Cotsapas C, Davies L, Price AL, de Bakker PI, Maller J, Pe'er I, Burtt NP, Blumenstiel B, DeFelice M, Parkin M, Barry R, Winslow W, Healy C, Graham RR, Neale BM, Izmailova E, Roubenoff R, Parker AN, Glass R, Karlson EW, Maher N, Hafler DA, Lee DM, Seldin MF, Remmers EF, Lee AT, Padyukov L, Alfredsson L, Coblyn J, Weinblatt ME, Gabriel SB, Purcell S, Klareskog L, Gregersen PK, Shadick NA, Daly MJ, Altshuler D. · Program in Medical and Population Genetics, Broad Institute of Harvard and Massachusetts Institute of Technology, Cambridge, Massachusetts 02142, USA. · Nat Genet. · Pubmed #17982456 links to  free full text

Abstract: 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, approximately 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 x 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.

Remmers EF, Plenge RM, Lee AT, Graham RR, Hom G, Behrens TW, de Bakker PI, Le JM, Lee HS, Batliwalla F, Li W, Masters SL, Booty MG, Carulli JP, Padyukov L, Alfredsson L, Klareskog L, Chen WV, Amos CI, Criswell LA, Seldin MF, Kastner DL, Gregersen PK. · National Institute of Arthritis and Musculoskeletal and Skin Diseases, Bethesda, MD, USA. · N Engl J Med. · Pubmed #17804842 links to  free full text

Abstract: BACKGROUND: Rheumatoid arthritis is a chronic inflammatory disease with a substantial genetic component. Susceptibility to disease has been linked with a region on chromosome 2q. METHODS: We tested single-nucleotide polymorphisms (SNPs) in and around 13 candidate genes within the previously linked chromosome 2q region for association with rheumatoid arthritis. We then performed fine mapping of the STAT1-STAT4 region in a total of 1620 case patients with established rheumatoid arthritis and 2635 controls, all from North America. Implicated SNPs were further tested in an independent case-control series of 1529 patients with early rheumatoid arthritis and 881 controls, all from Sweden, and in a total of 1039 case patients and 1248 controls from three series of patients with systemic lupus erythematosus. RESULTS: A SNP haplotype in the third intron of STAT4 was associated with susceptibility to both rheumatoid arthritis and systemic lupus erythematosus. The minor alleles of the haplotype-defining SNPs were present in 27% of chromosomes of patients with established rheumatoid arthritis, as compared with 22% of those of controls (for the SNP rs7574865, P=2.81x10(-7); odds ratio for having the risk allele in chromosomes of patients vs. those of controls, 1.32). The association was replicated in Swedish patients with recent-onset rheumatoid arthritis (P=0.02) and matched controls. The haplotype marked by rs7574865 was strongly associated with lupus, being present on 31% of chromosomes of case patients and 22% of those of controls (P=1.87x10(-9); odds ratio for having the risk allele in chromosomes of patients vs. those of controls, 1.55). Homozygosity of the risk allele, as compared with absence of the allele, was associated with a more than doubled risk for lupus and a 60% increased risk for rheumatoid arthritis. CONCLUSIONS: A haplotype of STAT4 is associated with increased risk for both rheumatoid arthritis and systemic lupus erythematosus, suggesting a shared pathway for these illnesses.

Plenge RM, Seielstad M, Padyukov L, Lee AT, Remmers EF, Ding B, Liew A, Khalili H, Chandrasekaran A, Davies LR, Li W, Tan AK, Bonnard C, Ong RT, Thalamuthu A, Pettersson S, Liu C, Tian C, Chen WV, Carulli JP, Beckman EM, Altshuler D, Alfredsson L, Criswell LA, Amos CI, Seldin MF, Kastner DL, Klareskog L, Gregersen PK. · Broad Institute of Harvard and the Massachusetts Institute of Technology, Cambridge, MA, USA. · N Engl J Med. · Pubmed #17804836 links to  free full text

Abstract: BACKGROUND: Rheumatoid arthritis has a complex mode of inheritance. Although HLA-DRB1 and PTPN22 are well-established susceptibility loci, other genes that confer a modest level of risk have been identified recently. We carried out a genomewide association analysis to identify additional genetic loci associated with an increased risk of rheumatoid arthritis. METHODS: We genotyped 317,503 single-nucleotide polymorphisms (SNPs) in a combined case-control study of 1522 case subjects with rheumatoid arthritis and 1850 matched control subjects. The patients were seropositive for autoantibodies against cyclic citrullinated peptide (CCP). We obtained samples from two data sets, the North American Rheumatoid Arthritis Consortium (NARAC) and the Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA). Results from NARAC and EIRA for 297,086 SNPs that passed quality-control filters were combined with the use of Cochran-Mantel-Haenszel stratified analysis. SNPs showing a significant association with disease (P<1x10(-8)) were genotyped in an independent set of case subjects with anti-CCP-positive rheumatoid arthritis (485 from NARAC and 512 from EIRA) and in control subjects (1282 from NARAC and 495 from EIRA). RESULTS: We observed associations between disease and variants in the major-histocompatibility-complex locus, in PTPN22, and in a SNP (rs3761847) on chromosome 9 for all samples tested, the latter with an odds ratio of 1.32 (95% confidence interval, 1.23 to 1.42; P=4x10(-14)). The SNP is in linkage disequilibrium with two genes relevant to chronic inflammation: TRAF1 (encoding tumor necrosis factor receptor-associated factor 1) and C5 (encoding complement component 5). CONCLUSIONS: A common genetic variant at the TRAF1-C5 locus on chromosome 9 is associated with an increased risk of anti-CCP-positive rheumatoid arthritis.

Karlson EW, Chibnik LB, Cui J, Plenge RM, Glass RJ, Maher NE, Parker A, Roubenoff R, Izmailova E, Coblyn JS, Weinblatt ME, Shadick NA. · Division of Rheumatology, Allergy, and Immunology, Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts, USA. · Ann Rheum Dis. · Pubmed #17666451 links to  free full text

Abstract: BACKGROUND: HLA-DRB1 shared epitope (HLA-SE), PTPN22 and CTLA4 alleles are associated with cyclic citrullinated peptide (CCP) and rheumatoid arthritis (RA). OBJECTIVE: We examined associations between HLA-SE, PTPN22, CTLA4 genotypes and RA phenotypes in a large cohort to (a) replicate prior associations with CCP status, and (b) determine associations with radiographic erosions and age of diagnosis. METHODS: A total of 689 RA patients from the Brigham RA Sequential Study (BRASS) were genotyped for HLA-SE, PTPN22 (rs2476601) and CTLA4 (rs3087243). Association between genotypes and CCP, rheumatoid factor (RF) erosive phenotypes and age at diagnosis were assessed with multivariable models adjusting for age, sex and disease duration. Novel causal pathway analysis was used to test the hypothesis that genetic risk factors and CCP are in the causal pathway for predicting erosions. RESULTS: In multivariable analysis, presence of any HLA-SE was strongly associated with CCP+ (odds ratio (OR) 3.05, 95% CI 2.18-4.25), and RF+ (OR 2.53, 95% CI 1.83-3.5) phenotypes; presence of any PTPN22 T allele was associated with CCP+ (OR 1.81, 95% CI 1.24-2.66) and RF+ phenotypes (OR 1.84, 95% CI 1.27-2.66). CTLA4 was not associated with CCP or RF phenotypes. While HLA-SE was associated with erosive RA phenotype (OR 1.52, 95% CI 1.01-2.17), this was no longer significant after conditioning on CCP. PTPN22 and CTLA4 were not associated with erosive phenotype. Presence of any HLA-SE was associated with an average 3.6 years earlier diagnosis compared with absence of HLA-SE (41.3 vs 44.9 years, p = 0.002) and PTPN22 was associated with a 4.2 years earlier age of diagnosis (39.5 vs 43.6 years, p = 0.002). CTLA4 genotypes were not associated with age at diagnosis of RA. CONCLUSIONS: In this large clinical cohort, we replicated the association between HLA-SE and PTPN22, but not CTLA4 with CCP+ and RF+ phenotypes. We also found evidence for associations between HLA-SE, and PTPN22 and earlier age at diagnosis. Since HLA-SE is associated with erosive phenotype in unconditional analysis, but is not significant after conditioning on CCP, this suggests that CCP is in the causal pathway for predicting erosive phenotype.

Kallberg H, Padyukov L, Plenge RM, Ronnelid J, Gregersen PK, van der Helm-van Mil AH, Toes RE, Huizinga TW, Klareskog L, Alfredsson L, Anonymous00205. · Institute for Environmental Medicine, Karolinska Institutet, Stockholm, Sweden. · Am J Hum Genet. · Pubmed #17436241 links to  free full text

Abstract: Gene-gene and gene-environment interactions are key features in the development of rheumatoid arthritis (RA) and other complex diseases. The aim of this study was to use and compare three different definitions of interaction between the two major genetic risk factors of RA--the HLA-DRB1 shared epitope (SE) alleles and the PTPN22 R620W allele--in three large case-control studies: the Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA) study, the North American RA Consortium (NARAC) study, and the Dutch Leiden Early Arthritis Clinic study (in total, 1,977 cases and 2,405 controls). The EIRA study was also used to analyze interactions between smoking and the two genes. "Interaction" was defined either as a departure from additivity, as interaction in a multiplicative model, or in terms of linkage disequilibrium--for example, deviation from independence of penetrance of two unlinked loci. Consistent interaction, defined as departure from additivity, between HLA-DRB1 SE alleles and the A allele of PTPN22 R620W was seen in all three studies regarding anti-CCP-positive RA. Testing for multiplicative interactions demonstrated an interaction between the two genes only when the three studies were pooled. The linkage disequilibrium approach indicated a gene-gene interaction in EIRA and NARAC, as well as in the pooled analysis. No interaction was seen between smoking and PTPN22 R620W. A new pattern of interactions is described between the two major known genetic risk factors and the major environmental risk factor concerning the risk of developing anti-CCP-positive RA. The data extend the basis for a pathogenetic hypothesis for RA involving genetic and environmental factors. The study also raises and illustrates principal questions concerning ways to define interactions in complex diseases.

Plenge RM, Padyukov L, Remmers EF, Purcell S, Lee AT, Karlson EW, Wolfe F, Kastner DL, Alfredsson L, Altshuler D, Gregersen PK, Klareskog L, Rioux JD. · Broad Institute of MIT and Harvard, Cambridge, MA, USA. · Am J Hum Genet. · Pubmed #16380915 links to  free full text

Abstract: 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.