PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(12): e14326.
Published online 2010 December 15. doi:  10.1371/journal.pone.0014326
PMCID: PMC3002281

Genetic Variants in Toll-Like Receptors Are Not Associated with Rheumatoid Arthritis Susceptibility or Anti-Tumour Necrosis Factor Treatment Outcome

Janine Santos, Editor

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.

Conclusion

We conclude that genetic variants in TLRs do not play a major role in susceptibility for developing RA nor in anti-TNF treatment outcome in a Caucasian population.

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).

Table 1
Results of the case control association analysis in the Dutch discovery set.

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).

Figure 1
Meta-analysis of case control studies for TLR5 rs2072493 and TLR7 rs3853839.

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).

Table 2
Patient characteristics of RA patients treated with anti-TNF medication.
Table 3
Results of the association analysis for anti-TNF treatment outcome after three months of treatment initiation.

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.

(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.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

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.

References

1. Coenen MJ, Gregersen PK. Rheumatoid arthritis: a view of the current genetic landscape. Genes Immun. 2009;10:101–111. [PMC free article] [PubMed]
2. Raychaudhuri S, Thomson BP, Remmers EF, Eyre S, Hinks A, et al. Genetic variants at CD28, PRDM1 and CD2/CD58 are associated with rheumatoid arthritis risk. Nat Genet. 2009;41:1313–1318. [PMC free article] [PubMed]
3. Ospelt C, Brentano F, Rengel Y, Stanczyk J, Kolling C, et al. 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. 2008;58:3684–3692. [PubMed]
4. Roelofs MF, Joosten LA, bdollahi-Roodsaz S, van Lieshout AW, Sprong T, et al. 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. 2005;52:2313–2322. [PubMed]
5. van der Heijden I, Wilbrink B, Tchetverikov I, Schrijver IA, Schouls LM, et al. Presence of bacterial DNA and bacterial peptidoglycans in joints of patients with rheumatoid arthritis and other arthritides. Arthritis Rheum. 2000;43:593–598. [PubMed]
6. Akira S, Takeda K, Kaisho T. Toll-like receptors: critical proteins linking innate and acquired immunity. Nat Immunol. 2001;2:675–680. [PubMed]
7. Lee EY, Yim JJ, Lee HS, Lee YJ, Lee EB, et al. Dinucleotide repeat polymorphism in intron II of human Toll-like receptor 2 gene and susceptibility to rheumatoid arthritis. Int J Immunogenet. 2006;33:211–215. [PubMed]
8. Sanchez E, Orozco G, Lopez-Nevot MA, Jimenez-Alonso J, Martin J. Polymorphisms of toll-like receptor 2 and 4 genes in rheumatoid arthritis and systemic lupus erythematosus. Tissue Antigens. 2004;63:54–57. [PubMed]
9. Radstake TR, Franke B, Hanssen S, Netea MG, Welsing P, et al. 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. 2004;50:999–1001. [PubMed]
10. Sheedy FJ, Marinou I, O'Neill LA, Wilson AG. The Mal/TIRAP S180L and TLR4 G299D polymorphisms are not associated with susceptibility to, or severity of, rheumatoid arthritis. Ann Rheum Dis. 2008;67:1328–1331. [PubMed]
11. Kang ES, Lee J. Genotypic analysis of Asp299Gly and Thr399Ile polymorphism of Toll-like receptor 4 in systemic autoimmune diseases of Korean population. Rheumatol Int. 2007;27:887–889. [PubMed]
12. Kilding R, Akil M, Till S, Amos R, Winfield J, et al. A biologically important single nucleotide polymorphism within the toll-like receptor-4 gene is not associated with rheumatoid arthritis. Clin Exp Rheumatol. 2003;21:340–342. [PubMed]
13. Jaen O, Petit-Teixeira E, Kirsten H, Ahnert P, Semerano L, et al. No evidence of major effects in several Toll-like receptor gene polymorphisms in rheumatoid arthritis. Arthritis Res Ther. 2009;11:R5. [PMC free article] [PubMed]
14. Etem EO, Elyas H, Ozgocmen S, Yildirim A, Godekmerdan A. The investigation of toll-like receptor 3, 9 and 10 gene polymorphisms in Turkish rheumatoid arthritis patients. Rheumatol Int. 2010 DOI: 10.1007/s00296-010-1472-8. [PubMed]
15. Kuuliala K, Orpana A, Leirisalo-Repo M, Kautiainen H, Hurme M, et al. Polymorphism at position +896 of the toll-like receptor 4 gene interferes with rapid response to treatment in rheumatoid arthritis. Ann Rheum Dis. 2006;65:1241–1243. [PMC free article] [PubMed]
16. Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, et al. Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci U S A. 2009;106:9362–9367. [PubMed]
17. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007;447:661–678. [PMC free article] [PubMed]
18. Plenge RM, Seielstad M, Padyukov L, Lee AT, Remmers EF, et al. TRAF1-C5 as a risk locus for rheumatoid arthritis–a genomewide study. N Engl J Med. 2007;357:1199–1209. [PMC free article] [PubMed]
19. De Rycke L, Vandooren B, Kruithof E, De Keyser F, Veys EM, et al. 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. 2005;52:2146–2158. [PubMed]
20. Welsing PM, van Riel PL. The Nijmegen inception cohort of early rheumatoid arthritis. J Rheumatol. 2004;Suppl 69:14–21. [PubMed]
21. Toonen EJ, Coenen MJ, Kievit W, Fransen J, Eijsbouts AM, et al. 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. 2008;67:1174–1177. [PubMed]
22. Hoogendoorn EH, Hermus AR, de VF, Ross HA, Verbeek AL, et al. Thyroid function and prevalence of anti-thyroperoxidase antibodies in a population with borderline sufficient iodine intake: influences of age and sex. Clin Chem. 2006;52:104–111. [PubMed]
23. Kastbom A, Bratt J, Ernestam S, Lampa J, Padyukov L, et al. Fcgamma receptor type IIIA genotype and response to tumor necrosis factor alpha-blocking agents in patients with rheumatoid arthritis. Arthritis Rheum. 2007;56:448–452. [PubMed]
24. Enevold C, Radstake TR, Coenen MJ, Fransen J, Toonen EJ, et al. Multiplex screening of 22 single-nucleotide polymorphisms in 7 Toll-like receptors: an association study in rheumatoid arthritis. J Rheumatol. 2010;37:905–910. [PubMed]
25. Taylor JD, Briley D, Nguyen Q, Long K, Iannone MA, et al. Flow cytometric platform for high-throughput single nucleotide polymorphism analysis. Biotechniques. 2001;30:661–669. [PubMed]
26. Ye F, Li MS, Taylor JD, Nguyen Q, Colton HM, et al. Fluorescent microsphere-based readout technology for multiplexed human single nucleotide polymorphism analysis and bacterial identification. Hum Mutat. 2001;17:305–316. [PubMed]
27. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–575. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science