Lack of validation of genetic variants associated with anti–tumor necrosis factor therapy response in rheumatoid arthritis: a genome-wide association study replication and meta-analysis
- Ana Márquez†1Email author,
- Aida Ferreiro-Iglesias†2,
- Cristina L Dávila-Fajardo1, 3,
- Ariana Montes2,
- Dora Pascual-Salcedo4,
- Eva Perez-Pampin2,
- Manuel J Moreno-Ramos5,
- Rosa García-Portales6,
- Federico Navarro7,
- Virginia Moreira7,
- César Magro8,
- Rafael Caliz9,
- Miguel Angel Ferrer9,
- Juan José Alegre-Sancho10,
- Beatriz Joven11,
- Patricia Carreira11,
- Alejandro Balsa12,
- Yiannis Vasilopoulos13,
- Theologia Sarafidou13,
- José Cabeza-Barrera3,
- Javier Narvaez14,
- Enrique Raya8,
- Juan D Cañete15,
- Antonio Fernández-Nebro16,
- María del Carmen Ordóñez16,
- Arturo R de la Serna17,
- Berta Magallares17,
- Juan J Gomez-Reino2, 18,
- Antonio González2 and
- Javier Martín1
© Márquez et al.; licensee BioMed Central Ltd. 2014
Received: 9 October 2013
Accepted: 25 February 2014
Published: 11 March 2014
In this study, our aim was to elucidate the role of four polymorphisms identified in a prior large genome-wide association study (GWAS) in which the investigators analyzed the responses of patients with rheumatoid arthritis (RA) to treatment with tumor necrosis factor inhibitors (TNFi). The authors of that study reported that the four genetic variants were significantly associated. However, none of the associations reached GWAS significance, and two subsequent studies failed to replicate these associations.
The four polymorphisms (rs12081765, rs1532269, rs17301249 and rs7305646) were genotyped in a total of 634 TNFi-treated RA patients of Spanish Caucasian origin. Four outcomes were evaluated: changes in the Disease Activity Score in 28 joints (DAS28) after 6 and 12 months of treatment and classification according to the European League Against Rheumatism (EULAR) response criteria at the same time points. Association with DAS28 changes was assessed by linear regression using an additive genetic model. Contingency tables of genotype and allele frequencies between EULAR responder and nonresponder patients were compared. In addition, we combined our data with those of previously reported studies in a meta-analysis including 2,998 RA patients.
None of the four genetic variants showed an association with response to TNFi in any of the four outcomes analyzed in our Spanish patients. In addition, only rs1532269 yielded a suggestive association (P = 0.0033) with the response to TNFi when available data from previous studies were combined in the meta-analysis.
Our data suggest that the rs12081765, rs1532269, rs17301249 and rs7305646 genetic variants do not have a role as genetic predictors of TNFi treatment outcomes.
Rheumatoid arthritis (RA) is a systemic autoimmune disease characterized by chronic inflammation of the synovial joints resulting in joint destruction, polyarthritis and functional disability. This inflammatory condition affects approximately 1% of the Caucasian population, making it a significant cause of comorbidity and mortality .
In recent years, the use of tumor necrosis factor inhibitors (TNFi) has resulted in an improvement in the treatment of RA patients by reducing both inflammation and joint damage [2–4], and their clinical use has become widespread. However, a percentage of patients do not respond adequately to this therapy; therefore, the current use of these agents is based on a trial-and-error approach [5, 6]. Given the adverse effects and the high cost of this type of therapy, the establishment of pharmacogenetic markers to predict the response to TNFi treatment is a highly desirable goal.
Recently, researchers in pharmacogenetic studies have reported several genetic variants associated with clinical response to treatment with TNFi [7–11]. However, to date, only the PTPRC and PDE3A-SLCO1C1 loci have been associated in more than a single study [12–14].
In 2011, Plant et al.  performed a genome-wide association study (GWAS) in a large number of RA patients from the United Kingdom treated with TNFi. These investigators used a three-stage study design. The meta-analysis combining stages 1, 2 and 3 cohorts yielded four single-nucleotide polymorphisms (SNPs) putatively associated with the TNFi response at 6 months, although these associations did not reach the GWAS significance level. Two of the genetic variants mapped within genes, PDZ domain-containing 2 (PDZD2) and eyes absent homolog 4 (EYA4), and two polymorphisms mapped to intergenic regions on chromosomes 1 and 12. However, two subsequent GWASs conducted in European RA patients, whose treatment response was evaluated at 14 weeks, failed to replicate association with any of the four loci[9, 10].
The aim of our study was to assess the role of the four genetic variants identified by Plant et al. with regard to their association with response to TNFi using a large number of RA Spanish patients, as well as to conduct a meta-analysis including previous data.
Baseline characteristics of the rheumatoid arthritis patients a
Collection 1 (N= 438)
Collection 2 (N= 196)
Age, mean ± SD years
61.0 ± 12.04
56.3 ± 14.77
Age at diagnosis, mean ± SD years
44.95 ± 13.03
42.99 ± 13.69
Females, n (%)
Disease duration, mean ± SD years
16.45 ± 8.34
12.93 ± 8.47
Rheumatoid factor–positive, n (%)
Anti-CCP-positive, n (%)
Smoking status, n (%)
Health status, mean ± SD
DAS28 at baseline
5.86 ± 1.12
5.36 ± 1.13
Treatment, n (%)
Previous biologic agents
Anti-TNF drugs, n (%)
EULAR-defined response at 6 months, n (%)
EULAR-defined response at 12 months, n (%)
The Disease Activity Score in 28 joints (DAS28) was measured at baseline and at 6 and 12 months after the first TNFi infusion. Two scales were considered to assess the efficacy of the TNFi therapy. First, the absolute change in DAS28 (∆DAS28 = DAS28end – DAS28baseline) at 6 and 12 months of follow-up. Second, patients were classified as responders (good and moderate) or nonresponders at 6 and 12 months according to the European League Against Rheumatism (EULAR) response criteria .
Genomic DNA was extracted from peripheral white blood cells or saliva using standard procedures. Four SNPs—rs1532269 and rs17301249, intronic polymorphisms mapped within PDZD2 and EYA4, respectively, and rs12081765 and rs7305646 located at intergenic regions on chromosomes 1 and 12, respectively—were genotyped using a single-base extension technology (SNaPshot Multiplex Kit; Applied Biosystems, Foster City, CA, USA) in a multiplex PCR experiment (KAPA2G Fast HotStart; Kapa Biosystems, Wilmington, MA, USA) in collection 1 and using TaqMan allelic discrimination assays on a 7900HT Fast Real-Time PCR System, both purchased from Applied Biosystems, in collection 2. A deviation from Hardy–Weinberg equilibrium (HWE) was detected for the rs1532269 polymorphism in collection 2, so that SNP was genotyped in this sample set using the same methodology used for collection 1.
Power calculations were performed using Quanto version 1.2.4 software . All SNPs were tested for deviations from HWE. The association between SNPs and responses to TNFi was evaluated in two ways. In the first method, linear regression analysis using ∆DAS28 as the continuous dependent variable was performed under an additive genetic model using Plink version 1.07 statistical software . A t-test was used to identify polymorphisms associated with the response. In the second method, genotype and allele frequencies between EULAR-defined responder and nonresponder patients were compared. Plink was used to create 2 × 2 or 2 × 3 contingency tables and to perform a χ2 test and/or Fisher’s exact test. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated according to Woolf’s method. Because our present study is a replication study, no correction was applied to the obtained P-values when TNFi response was evaluated at 6 months. In the analyses involving the TNFi efficacy at 12 months, however, P-values were corrected by using the Benjamini–Hochberg step-up procedure to control for false discovery rates (FDRs) in multiple testing . The results were considered statistically significant when P-values were lower than 0.05.
Clinical variables previously identified as being independent predictors of efficacy of TNFi, including age, gender, smoking status, rheumatoid factor status, anti–cyclic citrullinated peptide antibody (anti-CCP) status, DAS28 at baseline, concurrent and previous treatment and TNFi, were assessed for association with treatment response in a multivariate regression analysis using STATISTICA version 7.0 software (StatSoft, Tulsa, OK, USA) and Plink software in collections 1 and 2, respectively. Only baseline DAS28, gender and TNFi were associated with the efficacy of the therapy. Accordingly, analyses were adjusted for these three variables.
The analysis of the combined data from our study and the previous reports [8–10] was performed using Plink. Heterogeneity between studies was assessed using Cochran’s Q and I2 statistics . Pooled analyses were performed by using the Mantel–Haenszel test under the fixed-effects model or the DerSimonian–Laird test under the random-effects model when heterogeneity was detected. The results were considered to be statistically significant when P-values were lower than 5 × 105.00E−08.
All of the four studied polymorphisms conformed to HWE expectations (P > 0.01). The genotyping success rate was higher than 95%.
Association of the four single-nucleotide polymorphisms with changes in Disease Activity Score in 28 joints at 6 and 12 months in Spanish rheumatoid arthritis patients a
Collection 1 (N= 419)
Collection 2 (N= 193)
Collection 1 (N= 314)
Collection 2 (N= 134)
In the subsequent analysis in collection 2, none of the tested polymorphisms showed an association with ∆DAS28 at 6 months (Table 2) (P = 0.995, P = 0.830, P = 0.458 and P = 0.661 for rs12081765, rs1532269, rs17301249 and rs7305646, respectively) or in the stratified analysis according to the EULAR-defined response (see Additional file 1: Table S1). When TNFi efficacy was evaluated at 12 months, the rs1532269 polymorphism showed an association with ∆DAS28 at that time point (Table 2) (P = 0.022, β = 0.335); however, statistical significance was lost after correction using the Benjamini–Hochberg step-up procedure for FDR (PFDR = 0.087). No association of this SNP with the EULAR-defined response at the 12-month time point was observed (see Additional file 1: Table S2).
Pooled analysis of the tested polymorphisms in the two Spanish cohorts
Meta-analysis of all available studies
Meta-analysis of association of four single-nucleotide polymorphisms with changes in Disease Activity Score in 28 joints in all available studies a
Plant et al. 
Cohort 1 (n = 566)
Cohort 2 (n = 379)
Cohort 3 (n = 341)
Krintel et al. 
Copenhagen cohort (n = 196)
Umicevic Mirkov et al. 
Dutch stage 1 cohort (n = 882)
Collection 1 (n = 438)
Collection 2 (n = 196)
N = 2,998
Our results make it unlikely that any of the four polymorphisms identified by Plant et al. could be used as genetic predictors of TNFi treatment outcomes, because they were not associated in our large Spanish RA study. This lack of association represents a very relevant finding because, to the best of our knowledge, our present study is the first in which the association of these SNPs was analyzed with the same treatment outcomes used for their identification. In addition, the combined analysis with the three previous studies included in our meta-analysis [8–10] showed only a suggestive association of one of the four polymorphisms (even weaker than that reported in the study by Plant et al. ). These findings seem to exclude effects of sufficient magnitude to be useful in predicting response to treatment.
The lack of replication in our study could be ascribed to multiple differences between studies. It is commonly impossible to identify one of them as being more relevant than the others. Genetic differences between populations are an unlikely explanation of the results, given that the allele frequencies of the four tested polymorphisms were very similar between studies. Clinical differences between the patients with RA included in the different reports are possible and difficult to exclude. In this regard, it has already been mentioned that Plant et al.  evaluated the response to TNFi at 6 months, whereas the two subsequent studies used the response at 14 weeks. However, this difference does not apply to our study in which evaluation at 6 months evidenced negative results. A common cause of discrepant results is the overestimation of the true effect in the discovery study. This phenomenon has been characterized as the “winner’s curse,” and it has been very prevalent in genetic association studies . It should be noted that the four SNPs studied by Plant et al.  showed the highest effects in the discovery cohort (which was the only one with a clear association between these four polymorphisms and the clinical response), whereas the three replication studies showed lower effect sizes (β-values less different from zero), thus supporting this possibility. Indeed, significant heterogeneity between studies was observed in the meta-analysis of three of the four analyzed genetic variants. Interestingly, this heterogeneity disappeared when the discovery cohort of Plant et al. was removed . Therefore, variables other than the presence of the four SNPs considered herein could have influenced the efficacy of TNFi in this cohort, accounting for its singularity.
Other GWASs of responses to TNFi treatment in RA have been published [7, 9–11]. This approach represents an important step forward in the understanding of the influence of genetic variability on the efficacy of this therapy. Only one of the observed associations has been found to reach the GWAS statistical significance level, however, and only after combination with data derived from replication studies . This highlights the important role of validation studies in determining the status of the remaining GWAS findings. It is to be expected that these combined efforts will produce useful insights.
The association of four polymorphisms (rs12081765, rs1532269, rs17301249 and rs7305646) previously identified as being associated with TNFi treatment response was not confirmed in the present study. Our results indicate that these four genetic variants are not useful predictors of response to TNFi in patients with RA.
Absolute change in Disease Activity Score in 28 joints
American College of Rheumatology
Cyclic citrullinated peptide
Disease Activity Score in 28 joints
European League Against Rheumatism
Genome-wide association study
Tumor necrosis factor.
The authors thank Sofía Vargas and Sonia García (Instituto de Parasitología y Biomedicina López-Neyra, Consejo Superior de Investigaciones Científicas (CSIC), Granada, Spain) and Carmen Pena (Instituto de Investigación Sanitaria, Santiago de Compostela, Spain) for their excellent technical assistance, Manolo Calaza (Instituto de Investigación Sanitaria, Santiago de Compostela, Spain) for help with data management and all the patients for kindly consenting to participate in the study.
- Scott DL, Wolfe F, Huizinga TW: Rheumatoid arthritis. Lancet. 2010, 376: 1094-1108. 10.1016/S0140-6736(10)60826-4.View ArticlePubMedGoogle Scholar
- Breedveld FC, Weisman MH, Kavanaugh AF, Cohen SB, Pavelka K, van Vollenhoven R, Sharp J, Perez JL, Spencer-Green GT: The PREMIER study: a multicenter, randomized, double-blind clinical trial of combination therapy with adalimumab plus methotrexate versus methotrexate alone or adalimumab alone in patients with early, aggressive rheumatoid arthritis who had not had previous methotrexate treatment. Arthritis Rheum. 2006, 54: 26-37. 10.1002/art.21519.View ArticlePubMedGoogle Scholar
- Emery P, Breedveld FC, Hall S, Durez P, Chang DJ, Robertson D, Singh A, Pedersen RD, Koenig AS, Freundlich B: Comparison of methotrexate monotherapy with a combination of methotrexate and etanercept in active, early, moderate to severe rheumatoid arthritis (COMET): a randomised, double-blind, parallel treatment trial. Lancet. 2008, 372: 375-382. 10.1016/S0140-6736(08)61000-4.View ArticlePubMedGoogle Scholar
- Smolen JS, Van Der Heijde DM, St Clair EW, Emery P, Bathon JM, Keystone E, Maini RN, Kalden JR, Schiff M, Baker D, Han C, Han J, Bala M, Active-Controlled Study of Patients Receiving Infliximab for the Treatment of Rheumatoid Arthritis of Early Onset (ASPIRE) Study Group: Predictors of joint damage in patients with early rheumatoid arthritis treated with high-dose methotrexate with or without concomitant infliximab: results from the ASPIRE trial. Arthritis Rheum. 2006, 54: 702-710. 10.1002/art.21678.View ArticlePubMedGoogle Scholar
- Gibbons LJ, Hyrich KL: Biologic therapy for rheumatoid arthritis: clinical efficacy and predictors of response. BioDrugs. 2009, 23: 111-124. 10.2165/00063030-200923020-00004.View ArticlePubMedGoogle Scholar
- Hetland ML, Christensen IJ, Tarp U, Dreyer L, Hansen A, Hansen IT, Kollerup G, Linde L, Lindegaard HM, Poulsen UE, Schlemmer A, Jensen DV, Jensen S, Hostenkamp G, Østergaard M, All Departments of Rheumatology in Denmark: Direct comparison of treatment responses, remission rates, and drug adherence in patients with rheumatoid arthritis treated with adalimumab, etanercept, or infliximab: results from eight years of surveillance of clinical practice in the nationwide Danish DANBIO registry. Arthritis Rheum. 2010, 62: 22-32. 10.1002/art.27227.View ArticlePubMedGoogle Scholar
- 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: Genome-wide association scan identifies candidate polymorphisms associated with differential response to anti-TNF treatment in rheumatoid arthritis. Mol Med. 2008, 14: 575-581.PubMed CentralView ArticlePubMedGoogle Scholar
- Plant D, Bowes J, Potter C, Hyrich KL, Morgan AW, Wilson AG, Isaacs JD, Barton A: Genome-wide association study of genetic predictors of anti-tumor necrosis factor treatment efficacy in rheumatoid arthritis identifies associations with polymorphisms at seven loci. Arthritis Rheum. 2011, 63: 645-653. 10.1002/art.30130.PubMed CentralView ArticlePubMedGoogle Scholar
- Krintel SB, Palermo G, Johansen JS, Germer S, Essioux L, Benayed R, Badi L, Østergaard M, Hetland ML: Investigation of single nucleotide polymorphisms and biological pathways associated with response to TNFα inhibitors in patients with rheumatoid arthritis. Pharmacogenet Genomics. 2012, 22: 577-589. 10.1097/FPC.0b013e3283544043.View ArticlePubMedGoogle Scholar
- Umiċeviċ Mirkov M, Cui J, Vermeulen SH, Stahl EA, Toonen EJ, Makkinje RR, Lee AT, Huizinga TW, Allaart R, Barton A, Mariette X, Miceli CR, Criswell LA, Tak PP, de Vries N, Saevarsdottir S, Padyukov L, Bridges SL, van Schaardenburg DJ, Jansen TL, Dutmer EA, van de Laar MA, Barrera P, Radstake TR, van Riel PL, Scheffer H, Franke B, Brunner HG, Plenge RM, Gregersen PK: Genome-wide association analysis of anti-TNF drug response in patients with rheumatoid arthritis. Ann Rheum Dis. 2013, 72: 1375-1381. 10.1136/annrheumdis-2012-202405.View ArticlePubMedGoogle Scholar
- Cui J, Stahl EA, Saevarsdottir S, Miceli C, Diogo D, Trynka G, Raj T, Mirkov MU, Canhao H, Ikari K, Terao C, Okada Y, Wedrén S, Askling J, Yamanaka H, Momohara S, Taniguchi A, Ohmura K, Matsuda F, Mimori T, Gupta N, Kuchroo M, Morgan AW, Isaacs JD, Wilson AG, Hyrich KL, Herenius M, Doorenspleet ME, Tak PP, Crusius JB: Genome-wide association study and gene expression analysis identifies CD84 as a predictor of response to etanercept therapy in rheumatoid arthritis. PLoS Genet. 2013, 9: e1003394-10.1371/journal.pgen.1003394.PubMed CentralView ArticlePubMedGoogle Scholar
- Acosta-Colman I, Palau N, Tornero J, Fernández-Nebro A, Blanco F, González-Alvaro I, Cañete JD, Maymó J, Ballina J, Fernández-Gutiérrez B, Olivé A, Corominas H, Erra A, Canela-Xandri O, Alonso A, López Lasanta M, Tortosa R, Julià A, Marsal S: GWAS replication study confirms the association of PDE3A-SLCO1C1 with anti-TNF therapy response in rheumatoid arthritis. Pharmacogenomics. 2013, 14: 727-734. 10.2217/pgs.13.60.View ArticlePubMedGoogle Scholar
- Cui J, Saevarsdottir S, Thomson B, Padyukov L, van der Helm-van Mil AH, Nititham J, Hughes LB, de Vries N, Raychaudhuri S, Alfredsson L, Askling J, Wedrén S, Ding B, Guiducci C, Wolbink GJ, Crusius JB, van der Horst-Bruinsma IE, Herenius M, Weinblatt ME, Shadick NA, Worthington J, Batliwalla F, Kern M, Morgan AW, Wilson AG, Isaacs JD, Hyrich K, Seldin MF, Moreland LW, Behrens TW: Rheumatoid arthritis risk allele PTPRC is also associated with response to anti-tumor necrosis factor α therapy. Arthritis Rheum. 2010, 62: 1849-1861.PubMed CentralPubMedGoogle Scholar
- Plant D, Prajapati R, Hyrich KL, Morgan AW, Wilson AG, Isaacs JD, Barton A, Biologics in Rheumatoid Arthritis Genetics and Genomics Study Syndicate: Replication of association of the PTPRC gene with response to anti-tumor necrosis factor therapy in a large UK cohort. Arthritis Rheum. 2012, 64: 665-670. 10.1002/art.33381.PubMed CentralView ArticlePubMedGoogle Scholar
- Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, Medsger TA, Mitchell DM, Neustadt DH, Pinals RS, Schaller JG, Sharp JT, Wilder RL, Hunder GG: The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988, 31: 315-324. 10.1002/art.1780310302.View ArticlePubMedGoogle Scholar
- van Gestel AM, Prevoo ML, Van’ t Hof MA, van Rijswijk MH, van de Putte LB, van Riel PL: Development and validation of the European League against Rheumatism response criteria for rheumatoid arthritis, comparison with the preliminary American College of Rheumatology and the World Health Organization/International League against Rheumatism criteria. Arthritis Rheum. 1996, 39: 34-40. 10.1002/art.1780390105.View ArticlePubMedGoogle Scholar
- Gauderman WJ, Morrison JM: QUANTO 1.1: A computer program for power and sample size calculations for genetic-epidemiology studies (2006). [http://hydra.usc.edu/gxe/]
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.PubMed CentralView ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995, 57: 289-300.Google Scholar
- Higgins JP, Thompson SG, Deeks JJ, Altman DG: Measuring inconsistency in meta-analyses. BMJ. 2003, 327: 557-560. 10.1136/bmj.327.7414.557.PubMed CentralView ArticlePubMedGoogle Scholar
- Lohmueller KE, Pearce CL, Pike M, Lander ES, Hirschhorn JN: Meta-analysis of genetic association studies supports a contribution of common variants to susceptibility to common disease. Nat Genet. 2003, 33: 177-182. 10.1038/ng1071.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/ applies to the data made available in this article, unless otherwise stated.