Figures
Abstract
Background
Several studies point to a role of Toll-like receptors (TLRs) in the development of rheumatoid arthritis (RA). We investigated if genetic variants in TLR genes are associated with RA and response to tumour necrosis factor blocking (anti-TNF) medication.
Methodology and Principal Findings
22 single nucleotide polymorphisms (SNPs) in seven TLR genes were genotyped in a Dutch cohort consisting of 378 RA patients and 294 controls. Significantly associated variants were investigated in replication cohorts from The Netherlands, United Kingdom and Sweden (2877 RA patients and 2025 controls). 182 of the Dutch patients were treated with anti-TNF medication. Using these patients and a replication cohort (269 Swedish patients) we analysed if genetic variants in TLR genes were associated with anti-TNF outcome. In the discovery phase of the study we found a significant association of SNPs rs2072493 in TLR5 and rs3853839 in TLR7 with RA disease susceptibility. Meta-analysis of discovery and replication cohorts did not confirm these findings. SNP rs2072493 in TLR5 was associated with anti-TNF outcome in the Dutch but not in the Swedish population.
Citation: Coenen MJH, Enevold C, Barrera P, Schijvenaars MMVAP, Toonen EJM, Scheffer H, et al. (2010) Genetic Variants in Toll-Like Receptors Are Not Associated with Rheumatoid Arthritis Susceptibility or Anti-Tumour Necrosis Factor Treatment Outcome. PLoS ONE 5(12): e14326. https://doi.org/10.1371/journal.pone.0014326
Editor: Janine Santos, University of Medicine and Dentistry of New Jersey, United States of America
Received: July 12, 2010; Accepted: November 22, 2010; Published: December 15, 2010
Copyright: © 2010 Coenen et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by The Danish Biotechnology program (C. Enevold) and The Lundbeck Foundation (C. Enevold) and a personal grant to M. Coenen from the Netherlands Organisation for Scientific Research (grant 916.76.020). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Rheumatoid arthritis (RA) is a severe chronic inflammatory disorder leading to joint damage. The causes of RA are largely unknown, however, the role of genetic factors is evident, with the MHC region as the major contributor [1]. To date, more than 20 non-MHC regions explaining approximately one third of the genetic contribution to RA have been identified [2].
The role of Toll-like receptors (TLRs) in the development of RA is supported by several studies. Synovial fibroblasts of RA patients constitutively express TLR 1-6, and it has been demonstrated that TLR 2, 3, 4 and 7 are up-regulated in RA synovial tissue compared to that of osteoarthritis patients and healthy controls [3]. Furthermore, the ligands for TLRs have been detected in synovium of RA patients [4],[5]. Since TLRs are potent activators of pro-inflammatory cytokines, including tumour necrosis factor (TNF) alpha, this also makes them interesting candidates for treatment outcome prediction, especially for TNF-neutralizing therapy [6].
Several groups explored the contribution of genetic variants in TLR genes to the development of RA. A Korean study showed an association of a dinucleotide repeat in TLR2 with RA [7]. One other study focusing on two SNPs in TLR2 resulting in amino acid substitutions (Arg677Trp and Arg753Gln) could not detect an association between these rare variants and RA [8]. Although conflicting results have been reported on the functional variant Asp299Gly (rs4986790) in TLR4 most studies showed no association of this variant with RA susceptibility [8]–[13]. A French group investigated 10 SNPs in TLR1, 2, 4, 6 and 9 in 100 families but did not find evidence for an association of variants in these genes with RA, autoantibody production or erosions [13]. A Turkish study, investigating variants in TLR3, 9 and 10 in 100 patients, showed an association between a variant in TLR9 and RA susceptibility [14]. So far, only one small pharmacogenetic study on TLRs has been published [15]. The authors did not demonstrate an association between the Asp299Gly variant in TLR4 and the response to disease modifying anti-rheumatic drugs (DMARDs).
Although most studies performed to date point to a lack of association between genetic variants in genes coding for TLRs and RA disease susceptibility or treatment response, these have either been very small or did not comprehensively test potentially functional SNPs in the genes. Therefore we performed an association study of TLRs testing involvement in RA pathogenesis and anti-TNF treatment response including seven TLR genes, three of which (TLR5, 7 and 8) have not been investigated in connection with RA before. The twenty-two SNPs investigated were selected for potential effects on gene function or regulation.
Results
Case control analysis
In the discovery phase of our study we included 378 RA patients and 294 controls and genotyped 22 SNPs in seven TLR genes. Genotyping failed for 51 samples (22 cases and 26 controls). Six SNPs were excluded from the analysis, based on a low MAF (Table 1).
Two of the investigated polymorphisms, rs2072493 in TLR5 and rs3853839 in TLR7, showed nominally significant allelic association with RA (p≤0.01; Table 1). The minor allele of both SNPs conferred a protective effect for the development of RA. A stratified analysis for either the shared epitope or rheumatoid arthritis. We did not identify SNPs that were specifically associated with one of the RA subgroups. The odds ratios for SNPs rs2072493 and rs3853839 were comparable to the odds ratios observed in the complete patient population (data not shown). The frequencies of rs2072493 and rs3853839 were investigated in three additional cohorts from The Netherlands, Sweden and the UK, consisting of 2877 patients and 2025 controls. This analysis showed no association of the SNPs in TLR5 and TLR7 with RA aetiology either in the separate cohorts or in a meta-analysis (Figure 1).
Anti-TNF treatment response
A subsample of the Dutch patients were treated with anti-TNF medication and had EULAR response information available at 3 months (n = 182, Table 2). Results of the association analysis are shown in Table 3. TLR5 rs2072493 showed a significant association with anti-TNF treatment outcome (p = 0.003). In total, 182 patients were included in the analysis. The EULAR response criteria were used to asses association of a genotype with anti-TNF response. The test statistic and p-value were from Fisher's Exact Test for most SNPs, those indicated with * were based on Pearson Chi-square and p-value. A medication specific analysis showed that none of the investigated SNPs was associated with a specific anti-TNF (Table S1). Validation of association found with TLR5 in 269 RA patients treated with anti-TNF from Sweden did not confirm our initial findings (Fisher's exact test p = 0.63).
The differences observed between the Dutch and Swedish cohort were not related to the anti-TNF used (Fisher's exact test Swedish cohort for infliximab (p = 0.75) and etanercept (p = 0.49)).
Discussion
We found no evidence for association between 22 SNPs in TLR genes and RA disease susceptibility or anti-TNF treatment response, in keeping with the majority of the other candidate gene studies investigating TLRs in RA patients [7]–[13]. Furthermore, genome-wide association studies (GWAS) do not report associations between TLR genes and RA disease susceptibility (www.genome.gov/gwastudies). Accessed December 31th 2009), suggesting that genetic variants in these genes are not a major contributor to disease development [16]. However, some TLR genes were poorly covered in the initial GWAS (e.g. no SNPs in TLR5 and 9 were included in these studies) [17],[18]. The TLR genes remain strong candidates for susceptibility genes in RA given their crucial role in mediating inflammatory signalling. One possibility is that these genes might harbour rare variants with large effect sizes, which would not be well tagged with the common SNPs selected for investigation in the current study. With the introduction of next generation sequencing technologies it should become possible to analyse the genes in a thorough manner but much larger sample sizes will be required to ensure power to detect association with rare variants, even if effect sizes are larger.
We could not identify an association of the investigated SNPs with anti-TNF response. We were unable to replicate the initial finding in the Dutch cohort. Differences in ethnicity, gender, RF and co-medication might influence the observed discrepancies. Besides the study has low power to detect an association due to the small patient population. Despite this the function of TLRs makes them likely candidate genes, since they directly modulate TNF expression [6]. In addition, De Rijcke and colleagues showed that anti-TNF treatment in spondylarthropathy results in a reduced expression of TLR2 and TLR4 [19]. The TLR pathway should therefore be investigated in more detail before drawing definite conclusions about the contribution of these genes to anti-TNF treatment outcome.
In conclusion, our study does not support an association between genetic variants in TLR genes and RA disease susceptibility or anti-TNF treatment response.
Materials and Methods
Ethics Statement
The study was approved by the ethics committees of the participating institutions. For the Dutch samples the “Commissie Mensgebonden Onderzoek (CMO) Regio Arnhem Nijmegen” of the Radboud University Nijmegen Medical Centre approved the study (CMO number 2004/014). All Dutch patients and controls gave written informed consent prior to participation in the study. Verbal consent was received from all patients participating in the EIRA study and it was registered at clinical journals according to Swedish law. Written consent was received from healthy control individuals. The ethics committee of Karolinska Institutet approved the study. All UK patients and controls were recruited with ethical approval and provided written informed consent (North-West Multi-Centre Research Ethics Committee (MREC 99/8/84) and the University of Manchester Committee on the Ethics of Research on Human Beings).
Patients and controls
All patients in this study fulfilled the American College of Rheumatology (ACR) criteria. In the discovery phase 378 Dutch RA patient participating in the Nijmegen Inception cohort and 294 anonymous blood donors were included [20]. Cohorts from The Netherlands, Sweden and England were included in the replication phase. The Dutch replication cohort consisted of 313 RA patients participating in the Dutch Rheumatoid Arthritis Monitoring (DREAM) registry [21]. The Dutch controls (n = 527) were part of the Nijmegen Biomedical Study [22].
The Swedish patients (n = 2158) and controls (sex, age and residence area matched, n = 1068) were part of the Epidemiological Investigation of Rheumatoid arthritis (EIRA) study. The UK patients (n = 406) were recruited as part of the Arthritis Research Campaign National Repository of patients and families with one patient from each family selected for genotyping. UK controls without a history of inflammatory arthritis (n = 430) were recruited from General Practitioner records or from blood donors.
Anti-TNF treatment outcome
Patients participating in the DREAM registry were included in the study of anti-TNF treatment response. The registry includes patients with RA who start with anti-TNF treatment, according to the Dutch recommendations (Disease Activity Score (DAS28) >3.2 and previous failure on at least two disease-modifying antirheumatic drugs (DMARDs), one of which has to be methotrexate (MTX)) [21]. The DAS28 at treatment start and after 3 months of treatment were used to calculate the clinical response according to the EULAR criteria [23]. Association analysis was performed for the total patient group and a subgroup analysis for the patients treated with infliximab and adalimumab. Nominally significant findings (p≤0.01) were validated in 269 anti-TNF treated RA patients from Sweden [23].
SNP genotyping
Twenty-two SNPs located in TLR2, 3, 4, 5, 7, 8 and 9 were genotyped in all Dutch patients and Dutch anonymous blood donors (Table 1). SNP selection criteria are described previously [24]. SNPs were genotyped using two multiplexed bead-based assays on a Luminex 100IS flow cytometer (Luminex Corporation, Austin, TX, USA). The tests were based on procedures previously described, with some modifications [24]–[26].
Replication studies were performed using TaqMan® SNP genotyping assays for SNPs rs2072493 in TLR5 (assay ID C__22273027_10) and rs3853839 in TLR7 (assay ID C___2259573_10).
Statistical analysis
For quality control reasons we excluded SNPs that had a minor allele frequency (MAF) <0.05. Hardy-Weinberg Equilibrium (HWE) was tested using the controls (only females were used for X-chromosomal SNPs), no deviations were observed. The whole genome association analysis toolset (PLINK) was used for association analysis at the allelic level [27]. A stratified analysis was performed for RF positive and negative as well as shared epitope positive and negative patients. SNPs with a nominal p-value ≤0.01 in the discovery sample were included in the replication phase. Review Manager 5 (Review Manager (RevMan) Version 5.0. Copenhagen: The Nordic Cochrane Centre, The Cochrane Collaboration, 2008) was used to perform a meta-analysis. Homogeneity of odds ratios (OR) among cohorts was calculated using the Breslow-Day method, and pooled ORs were calculated under a fixed effects model (Mantel-Haenszel meta-analysis).
To assess the effect of SNPs on anti-TNF response, a two-tailed Fisher's exact or Chi-square test was used (where applicable) for the comparison of the EULAR response scores between the three genotype groups in a 3×3 table in SPSS (SPSS for Windows, version 16.0). X-chromosomal SNPs were analyzed for males and females separately.
Supporting Information
Table S1.
Results of the association analysis for anti-TNF treatment outcome after three months of treatment initiation.
https://doi.org/10.1371/journal.pone.0014326.s001
(0.07 MB DOC)
Acknowledgments
We thank all the patients and controls for participation and members of EIRA group for support. The Nijmegen Biomedical Study is a population-based survey conducted at the Department of Epidemiology and Biostatistics and the Department of Clinical Chemistry of the Radboud University Nijmegen Medical Centre. Principal investigators of the Nijmegen Biomedical Study are L.A.L.M. Kiemeney, M. den Heijer, A.L.M. Verbeek, D.W. Swinkels and B. Franke.
Author Contributions
Conceived and designed the experiments: MJC CE PB EJT BF KB TR. Performed the experiments: CE MMS EJT. Analyzed the data: MJC CE EJT LP AK AB WK. Contributed reagents/materials/analysis tools: MJC HS LP AK LK AB WK MJR TLJ DWS PvR BF KB. Wrote the paper: MJC CE TR.
References
- 1. Coenen MJ, Gregersen PK (2009) Rheumatoid arthritis: a view of the current genetic landscape. Genes Immun 10: 101–111.
- 2. Raychaudhuri S, Thomson BP, Remmers EF, Eyre S, Hinks A, et al. (2009) Genetic variants at CD28, PRDM1 and CD2/CD58 are associated with rheumatoid arthritis risk. Nat Genet 41: 1313–1318.
- 3. Ospelt C, Brentano F, Rengel Y, Stanczyk J, Kolling C, et al. (2008) Overexpression of toll-like receptors 3 and 4 in synovial tissue from patients with early rheumatoid arthritis: toll-like receptor expression in early and longstanding arthritis. Arthritis Rheum 58: 3684–3692.
- 4. Roelofs MF, Joosten LA, bdollahi-Roodsaz S, van Lieshout AW, Sprong T, et al. (2005) The expression of toll-like receptors 3 and 7 in rheumatoid arthritis synovium is increased and costimulation of toll-like receptors 3, 4, and 7/8 results in synergistic cytokine production by dendritic cells. Arthritis Rheum 52: 2313–2322.
- 5. van der Heijden I, Wilbrink B, Tchetverikov I, Schrijver IA, Schouls LM, et al. (2000) Presence of bacterial DNA and bacterial peptidoglycans in joints of patients with rheumatoid arthritis and other arthritides. Arthritis Rheum 43: 593–598.
- 6. Akira S, Takeda K, Kaisho T (2001) Toll-like receptors: critical proteins linking innate and acquired immunity. Nat Immunol 2: 675–680.
- 7. Lee EY, Yim JJ, Lee HS, Lee YJ, Lee EB, et al. (2006) Dinucleotide repeat polymorphism in intron II of human Toll-like receptor 2 gene and susceptibility to rheumatoid arthritis. Int J Immunogenet 33: 211–215.
- 8. Sanchez E, Orozco G, Lopez-Nevot MA, Jimenez-Alonso J, Martin J (2004) Polymorphisms of toll-like receptor 2 and 4 genes in rheumatoid arthritis and systemic lupus erythematosus. Tissue Antigens 63: 54–57.
- 9. Radstake TR, Franke B, Hanssen S, Netea MG, Welsing P, et al. (2004) The Toll-like receptor 4 Asp299Gly functional variant is associated with decreased rheumatoid arthritis disease susceptibility but does not influence disease severity and/or outcome. Arthritis Rheum 50: 999–1001.
- 10. Sheedy FJ, Marinou I, O'Neill LA, Wilson AG (2008) The Mal/TIRAP S180L and TLR4 G299D polymorphisms are not associated with susceptibility to, or severity of, rheumatoid arthritis. Ann Rheum Dis 67: 1328–1331.
- 11. Kang ES, Lee J (2007) Genotypic analysis of Asp299Gly and Thr399Ile polymorphism of Toll-like receptor 4 in systemic autoimmune diseases of Korean population. Rheumatol Int 27: 887–889.
- 12. Kilding R, Akil M, Till S, Amos R, Winfield J, et al. (2003) A biologically important single nucleotide polymorphism within the toll-like receptor-4 gene is not associated with rheumatoid arthritis. Clin Exp Rheumatol 21: 340–342.
- 13. Jaen O, Petit-Teixeira E, Kirsten H, Ahnert P, Semerano L, et al. (2009) No evidence of major effects in several Toll-like receptor gene polymorphisms in rheumatoid arthritis. Arthritis Res Ther 11: R5.
- 14. Etem EO, Elyas H, Ozgocmen S, Yildirim A, Godekmerdan A (2010) The investigation of toll-like receptor 3, 9 and 10 gene polymorphisms in Turkish rheumatoid arthritis patients. Rheumatol Int. DOI:https://doi.org/10.1007/s00296-010-1472-8.
- 15. Kuuliala K, Orpana A, Leirisalo-Repo M, Kautiainen H, Hurme M, et al. (2006) Polymorphism at position +896 of the toll-like receptor 4 gene interferes with rapid response to treatment in rheumatoid arthritis. Ann Rheum Dis 65: 1241–1243.
- 16. Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, et al. (2009) Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci U S A 106: 9362–9367.
- 17. (2007) Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447: 661–678.
- 18. Plenge RM, Seielstad M, Padyukov L, Lee AT, Remmers EF, et al. (2007) TRAF1-C5 as a risk locus for rheumatoid arthritis–a genomewide study. N Engl J Med 357: 1199–1209.
- 19. De Rycke L, Vandooren B, Kruithof E, De Keyser F, Veys EM, et al. (2005) Tumor necrosis factor alpha blockade treatment down-modulates the increased systemic and local expression of Toll-like receptor 2 and Toll-like receptor 4 in spondylarthropathy. Arthritis Rheum 52: 2146–2158.
- 20. Welsing PM, van Riel PL (2004) The Nijmegen inception cohort of early rheumatoid arthritis. J Rheumatol Suppl 69: 14–21.
- 21. Toonen EJ, Coenen MJ, Kievit W, Fransen J, Eijsbouts AM, et al. (2008) The tumour necrosis factor receptor superfamily member 1b 676T>G polymorphism in relation to response to infliximab and adalimumab treatment and disease severity in rheumatoid arthritis. Ann Rheum Dis 67: 1174–1177.
- 22. Hoogendoorn EH, Hermus AR, de VF, Ross HA, Verbeek AL, et al. (2006) Thyroid function and prevalence of anti-thyroperoxidase antibodies in a population with borderline sufficient iodine intake: influences of age and sex. Clin Chem 52: 104–111.
- 23. Kastbom A, Bratt J, Ernestam S, Lampa J, Padyukov L, et al. (2007) Fcgamma receptor type IIIA genotype and response to tumor necrosis factor alpha-blocking agents in patients with rheumatoid arthritis. Arthritis Rheum 56: 448–452.
- 24. Enevold C, Radstake TR, Coenen MJ, Fransen J, Toonen EJ, et al. (2010) Multiplex screening of 22 single-nucleotide polymorphisms in 7 Toll-like receptors: an association study in rheumatoid arthritis. J Rheumatol 37: 905–910.
- 25. Taylor JD, Briley D, Nguyen Q, Long K, Iannone MA, et al. (2001) Flow cytometric platform for high-throughput single nucleotide polymorphism analysis. Biotechniques 30: 661–669.
- 26. Ye F, Li MS, Taylor JD, Nguyen Q, Colton HM, et al. (2001) Fluorescent microsphere-based readout technology for multiplexed human single nucleotide polymorphism analysis and bacterial identification. Hum Mutat 17: 305–316.
- 27. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.