BioMed Research International
Volume 2015 (2015), Article ID 490295, 9 pages
http://dx.doi.org/10.1155/2015/490295

Research Article

TRAF1/C5 but Not PTPRC Variants Are Potential Predictors of Rheumatoid Arthritis Response to Anti-Tumor Necrosis Factor Therapy

Helena Canhão,1,2,3 Ana Maria Rodrigues,1,2 Maria José Santos,1,4 Diana Carmona-Fernandes,¹ Bruno F. Bettencourt,5,6 Jing Cui,³ Fabiana L. Rocha,5,6 José Canas Silva,⁴ Joaquim Polido-Pereira,1,2 José Alberto Pereira Silva,² José António Costa,⁷ Domingos Araujo,⁷ Cândida Silva,⁸ Helena Santos,⁸ Cátia Duarte,⁹ Rafael Cáliz,¹⁰ Ileana Filipescu,¹¹ Fernando Pimentel-Santos,12,13 Jaime Branco,12,13 Juan Sainz,¹⁴ Robert M. Plenge,3,15 Daniel H. Solomon,3,16 Jácome Bruges-Armas,5,6 José António P. Da Silva,⁹ João Eurico Fonseca,1,2 and Elizabeth W. Karlson³

¹Rheumatology Research Unit, Instituto de Medicina Molecular, Edifício Egas Moniz, Faculdade de Medicina da Universidade de Lisboa, Avenida Egas Moniz, 1649-028 Lisbon, Portugal
²Rheumatology Department, Hospital de Santa Maria (CHLN), Lisbon Academic Medical Centre, 1649-035 Lisbon, Portugal
³Division of Rheumatology, Allergy and Immunology, Section of Clinical Sciences, Brigham and Women’s Hospital, Boston, MA 02115, USA
⁴Rheumatology Department, Hospital Garcia de Orta, 2805-267 Almada, Portugal
⁵SEEBMO, Hospital de Santo Espirito da Ilha Terceira, Azores, 9700-049 Angra do Heroísmo, Portugal
⁶Genetics & Arthritis Research Group (GARG), Institute for Molecular and Cell Biology (IBMC), 4150-180 Oporto, Portugal
⁷Rheumatology Department, Conde de Bertiandos Hospital (ULSAM), 4990-041 Ponte de Lima, Portugal
⁸Instituto Português de Reumatologia, Lisbon, Portugal
⁹Rheumatology Department, Centro Hospitalar da Universidade de Coimbra, 3000-076 Coimbra, Portugal
¹⁰Rheumatology Department, Virgen de las Nieves University Hospital, Granada, Spain
¹¹Rheumatology Department, University of Medicine and Pharmacy “Iuliu Hatieganu”, Cluj-Napoca, Romania
¹²CEDOC, Faculdade de Ciências Médicas, Universidade Nova de Lisboa, Lisbon, Portugal
¹³Rheumatology Department, Hospital de Egas Moniz (CHLO), 1349-019 Lisbon, Portugal
¹⁴Genomic Oncology Area, GENYO Centre for Genomics and Oncological Research, Pfizer/University of Granada/Andalusian Regional Government, PTS Granada, Granada, Spain
¹⁵The Broad Institute of Harvard and MIT, Cambridge, MA, USA
¹⁶Division of Pharmacoepidemiology, Brigham and Women’s Hospital, Boston, MA, USA

Received 25 July 2014; Accepted 26 November 2014

Academic Editor: María E. Esteban

Copyright © 2015 Helena Canhão et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. The aim of our work was to replicate, in a Southern European population, the association reported in Northern populations between PTPRC locus and response to anti-tumor necrosis factor (anti-TNF) treatment in rheumatoid arthritis (RA). We also looked at associations between five RA risk alleles and treatment response. Methods. We evaluated associations between anti-TNF treatment responses assessed by DAS28 change and by EULAR response at six months in 383 Portuguese patients. Univariate and multivariate linear and logistic regression analyses were performed. In a second step to confirm our findings, we pooled our population with 265 Spanish patients. Results. No association was found between PTPRC rs10919563 allele and anti-TNF treatment response, neither in Portuguese modeling for several clinical variables nor in the overall population combining Portuguese and Spanish patients. The minor allele for RA susceptibility, rs3761847 SNP in TRAF1/C5 region, was associated with a poor response in linear and logistic univariate and multivariate regression analyses. No association was observed with the other allellic variants. Results were confirmed in the pooled analysis. Conclusion. This study did not replicate the association between PTPRC and the response to anti-TNF treatment in our Southern European population. We found that TRAF1/C5 risk RA variants potentially influence anti-TNF treatment response.

1. Introduction

Rheumatoid arthritis (RA) is an inflammatory, chronic, and disabling disease. Methotrexate is the most widely used disease modifying antirheumatic drug (DMARD) in RA treatment. However, for refractory and severe cases, anti-tumor necrosis factor (anti-TNF) therapy has become a cornerstone of RA treatment strategy [1]. These drugs have revolutionized RA treatment and prognosis in the last 10–15 years. Nevertheless, only approximately one-third of patients achieve remission and the other third will eventually fail to respond [2]. In a multifactorial and polygenic disease like RA, it is expected that response to treatment may be influenced by genetic, clinical, and biological factors [3]. The identification of predictors of response is of crucial importance to optimize the cost-effective use of expensive medications, such as anti-TNF therapy. Large registries collecting information on sociodemographic characteristics, disease activity, functional status, and treatments have allowed the study of clinical predictors of response [4–7].

Genetic variants associated with RA susceptibility include the HLA-DRB1 region containing shared epitope alleles (SE), which is also associated with severity [8, 9]. Outside the major histocompatibility complex (MHC), PTPN22, TRAF1/C5, and TNFAIP3 loci were the most consistently associated with susceptibility and TRAF1/C5 region also with RA severity [10, 11] and noncardiovascular mortality in some populations [12].

In the nineties, studies performed to look at associations between treatment response and the presence of SE indicated that the response to disease modifying antirheumatic drugs (DMARDs) such as methotrexate (MTX), in combination or monotherapy [13], and more recently with leflunomide [14], etanercept [15], and infliximab [16], was better in the presence of SE.

In recent years, several studies of potential associations between anti-TNF treatment response and polymorphisms in the promoter region of the TNF gene (positions −308 and −238), and other related genes such as lymphotoxin-α and TNF receptors, showed contradictory results [17–20]. With the increasing knowledge on RA pathophysiology and genomewide studies demonstrating that loci related with TNF signaling pathways such as the NF-kB signaling pathway (TRAF1/C5, TNFAIP3, and REL) and other pathological processes such as enhanced citrullination (PADI4) may increase RA risk, it is compelling to explore how those loci could also influence the anti-TNF response [21–28].

Cui et al., analyzing thirty-one risk allele variants, found that the major allele (G) of the rs10919563 PTPRC locus, which is a known predictor of RA risk, was associated with an increased response to anti-TNF therapy, with stronger association in seropositive patients (either anticitrullinated peptides antibodies, ACPA, and/or rheumatoid factors (RF)) [29]. The authors did not find any association with treatment response among the other thirty RA-associated risk alleles studied. In that multicohort study, potential associations between response and HLA-DRB1 were not assessed. One study from BRAGGSS, UK, showed no association between HLA-DRB1 and PTPN22 variants and response to anti-TNF treatment [30]. However, another study from the UK confirmed the association between PTPRC variants and response in the entire cohort, reporting no significance in the ACPA positive group alone [31].

The challenge over the next years will be to identify the RA stages in which genetic variants exert their maximum influence and also to unveil their clinical significance and usefulness as potential therapeutic targets or biomarkers [10].

In this study we aimed to replicate in a Southern European population the association between rs10919563 PTPRC variants and the response to anti-TNF treatment found in previous studies. We also aimed to test whether HLA-DRB1 and other five selected RA susceptibility genes may influence the response to anti-TNF treatment—that is, potential associations between anti-TNF treatment response and risk RA loci related with NF-kB signaling pathway (TRAF1/C5, TNFAIP3, and REL), citrullination (PADI4), and the genetic variants inside the MHC (HLA-DRB1 *04 high-resolution (4-digit) genotyping) and outside the MHC (PTPN22 locus) with the strongest association with RA risk. The analyses were modeled adjusting for clinical variables that influenced treatment response.

2. Material and Methods

2.1. Patients

Primary analyses were performed upon Reuma.pt, the National Register for Rheumatic Diseases from the Portuguese Society of Rheumatology (SPR) established in 2008, which captures more than 90% of patients treated with biological therapies managed in rheumatology departments across Portugal [32]. The register is linked to the Biobanco-IMM [33]. Blood samples were collected from November 2010 up to May 2011 at six major centers. Information on disease activity and treatments has been collected by rheumatologists at every infusion for intravenous drugs and every 3 months for subcutaneous biologic therapies. The decision to initiate and maintain the treatment was guided by the SPR’s recommendations [34]. RA patients fulfilling the American College of Rheumatology (ACR) 1987 revised criteria [35] were eligible for this study whether they were treated with an anti-TNF agent as the first biologic therapy, had a follow-up of at least six months, or had a blood sample collected for DNA assessment. Patients with self-reported non-Caucasian ancestry and those with missing values for DAS28 at baseline or at six months were excluded. Reuma.pt was approved by National Board of Data Protection and Health National Directorate.

In a second step for confirming our findings, we pooled our Portuguese sample with Spanish RA patients from the Rheumatology Department of Virgen de las Nieves (Granada, Spain) and Reina Sofia (Córdoba, Spain) Hospitals, selected with the criteria described above.

The study was conducted in accordance with the regulations governing clinical trials such as the Declaration of Helsinki and was approved by the Hospitals’ Ethics Committees. Patients signed an informed consent for research use of their clinical data and blood samples.

2.2. DNA Extraction and Genotyping

DNA was extracted from whole blood using the QIAamp DNA Blood Mini Kit combined with the automated extraction device QIAcube (QIAGEN, Hilden, Germany), according to the manufacturer’s instructions.

[PTPRC] rs10919563 (A/G) was previously reported in association with anti-TNF treatment response [29, 31]. The other five SNP markers were selected based on (1) relevance for RA biologic NF-kB pathway, [6q23-TNFAIP3] rs10499194 (C/T), [REL] rs13031237 (G/T), and [TRAF1/C5] rs3761847 (A/G); (2) role on citrullination process, [PADI4] rs2240340 (C/T); and (3) strong association with RA risk, [PTPN22] rs2476601 (A/G). HLA-DRB1 *04 high-resolution genotyping was performed by PCR sequence-specific primers (SSP) using Olerup SSP DRB1*04 typing Kit (Olerup SSP AB) as described in the manufacturers’ protocols.

Samples were run with Luminex xMAP system (Tepnel Lifecodes). Allele call was obtained with Quicktype for Lifematch 2.6 software. The samples were genotyped using Taqman SNP genotyping assays (Applied Biosystems, Foster City, USA) as described in the manufacturers’ protocols.

For purposes of quality control, 95% of sample threshold and 95% genotyping success threshold were used. Exclusion criteria also included a minor allele frequency <0.1 and deviation from Hardy-Weinberg equilibrium. In the end, HLA-DRB1 *04 and six allele variants were analyzed.

2.3. Statistical Analysis

2.3.1. Outcome Measures and Covariates

Our primary outcome was the change in disease activity score in 28 joints including erythrocyte sedimentation rate measures (DAS28ESR) between the drug start date and six months of treatment (in our study, a positive variation means a decrease in disease activity at six months in comparison with baseline visit) [36]. Secondary outcome was the proportion of nonresponders versus good responders (excluding moderate responders) defined by the EULAR response criteria at six months [37]. Nonresponse was defined by an absolute change in DAS28 ≤ 0.6 or a change in DAS28 between 0.6 and 1.2 with a DAS28 at six months >5.1. Good response was defined as change in DAS28 of >1.2 and DAS28 at six months ≤3.2. In this study, nonresponse was the reference category, with logistic regression analysis modeling the probability of achieving good response. For both outcomes, covariates’ coefficient > 0 or odds ratio (OR) > 1 predicted favorable response.

The predictors of interest were the minor allele variants of PTPRC and the five SNPs described above as well as the presence of SE (DRB1 *0401/04/05/08 and DRB1 *1001), as previously described in the RA Portuguese population [8].

Covariates collected at drug start date (baseline visit) were gender, age, age at diagnosis, disease duration, years of education, smoking (ever/never), RF (positive/negative), ACPA (CCP2, positive/negative), extra-articular manifestations (yes/no), concomitant therapy with corticosteroids (yes/no), any disease modifying antirheumatic drugs (DMARDs) including MTX (yes/no), DAS28, Health Assessment Questionnaire (HAQ), and physician’s global assessment of disease activity (PhGA). We tested the association between these variables and change in DAS28 at six months by univariate analyses. Then, we built a multivariate linear model with the significant baseline clinical covariates (at a ). The variables that remained significant after adjustment entered the multivariate model for each SNP.

In a second step, in order to confirm our results, we pooled the Spanish population and looked at associations between the SNPs and treatment response.

2.3.2. Primary Analysis

Each SNP was tested for an association with the anti-TNF response, taking anti-TNF drugs as a group. Univariate linear regression analyses for the primary outcome and logistic regression for the secondary outcome were performed, using additive models. Homozygotes for major alleles were classified as 0, 1 for heterozygotes, and 2 for minor allele homozygotes. The response was also modeled with multivariate models including significant baseline clinical predictors for treatment response.

2.3.3. Secondary Analyses

For the loci that presented a significant association with treatment response in the primary analysis, we tested the same relationship by stratifying the patients in ACPA positive and ACPA negative groups.

The study had >80% power to detect a change of >0.6 in the DAS28 score (considered a clinically meaningful change [37]) for allele frequencies > 0.1 and a type I error of <0.05.

There were no assumptions about the direction of effect on treatment response. Results were considered significant for a two sided value <0.05. Analyses were performed using SAS version 9.3 (SAS Institute, Cary, NC, USA).

3. Results

PTPRC, additional five SNP markers, and HLA-DRB1 genotyping were assessed in 416 Portuguese patients. Twenty-seven reported a non-Caucasian ancestry and were excluded. Six were excluded for absence of DAS28 at baseline or at six months, leaving 383 patients for analysis. Table 1 describes the baseline demographic and clinical characteristics for the 383 RA Portuguese patients included in the study.

Table 1: Baseline demographic and clinical characteristics of the 383 rheumatoid arthritis patients treated with anti-TNF drugs.

TNF inhibitor

    Adalimumab
79 (20.6)

    Etanercept
139 (36.3)

    Golimumab
10 (2.6)

    Infliximab
155 (40.5)

Age (years)
52.5 (12.2)

Disease duration (years)
10.7 (8.9)

Female
343 (89.5)

Rheumatoid factor
290 (75.7)

ACPA
278 (72.6)

Extra-articular manifestations
91 (23.7)

Smoking-ever ()
71 (19.3)

Education (years)
7.1 (4.6)

DMARDs ()
346 (91.8)

MTX ()
310 (82.2)

Corticosteroids ()
277 (73.5)

DAS28 ESR
5.77 (1.1)

Physician global assessment (mm)
56.1 (16.5)

Health Assessment Questionnaire
1.45 (0.58)

Values shown are means (SD) or (%).

TNF: tumor necrosis factor; ACPA: anticitrullinated peptides antibodies; DMARDs: disease modifying antirheumatic drugs; MTX: methotrexate; DAS: disease activity score; ESR: erythrocyte sedimentation rate.

At six months, 119 (31.1%) patients were classified as good responders, 175 (45.7%) as moderate responders, and 89 (23.2%) as nonresponders according to the EULAR response criteria [37].

Number of years completed at school, HAQ, and DAS28 at baseline were found to have a significant association with treatment response at six months and were included in the multivariate models (Table 2).

Table 2: Multivariate model of baseline demographic and clinical variables as predictors of response to anti-TNF treatment at 6 months.

Baseline variables
Coefficient ( value)

Age (years)
−0.01 (0.25)

Disease duration (years)
0.01 (0.24)

Female gender
−0.01 (0.95)

ACPA (positive)
−0.07 (0.64)

Smoking (ever)
−0.06 (0.72)

Higher education (years)
0.03 (0.03)*

Corticosteroids (yes)
−0.16 (0.26)

DMARDs (yes)
0.05 (0.83)

Extra-artic. manif. (yes)
0.10 (0.49)

HAQ
−0.40 (0.002)*

PhGA (mm)
−0.001 (0.68)

DAS
0.60 (<0.001)*

value <0.05.

The probability of response to anti-tumor necrosis factor therapy at 6 months was modeled in a multivariate linear regression analysis.

The change in disease activity score assessing 28 joints between the baseline visit and the visit after 6 months of therapy was the continuous outcome.

ACPA: anticitrullinated peptides antibodies; DMARDs: disease modifying antirheumatic drugs; Extra-artic. manif.: extra-articular manifestations; HAQ: Health Assessment Questionnaire; PhGA: physician global assessment; DAS: disease activity score.

We did not find association between PTPRC rs10919563 and anti-TNF treatment response (Table 3).

Table 3: Minor allele genetic variants as predictors of response, analyzed by univariate and multivariate linear and logistic regression models.

SNP
Gene
Chr Position (bp)
Minor allele
MAF
Absolute change in DAS
EULAR good response versus nonresponse

(good responders = 119, nonresponders = 89)

Univariate additive linear model
Multivariate additive linear model
Univariate additive logistic model
Multivariate additive logistic model

Coef (95% CI)
Coef (95% CI)
OR (95% CI)
OR (95% CI)

value
value
value
value

rs2476601
PTPN22
1
A
0.11
0.28 (−0.01, 0.58)
0.12 (−0.15, 0.40)
1.68 (0.88, 3.21)
1.71 (0.82, 3.54)

114,179,091
0.058
0.38
0.12
0.15

rs2240340
PADI4
1
T
0.46
−0.11 (−0.29, 0.07)
−0.09 (−0.25, 0.08)
0.79 (0.54, 1.16)
0.81 (0.52, 1.24)

17,507,279
0.22
0.30
0.22
0.33

rs13031237
REL
2
T
0.45
−0.09 (−0.26, 0.09)
−0.14 (−0.31, 0.02)
0.91 (0.62, 1.33)
0.77 (0.50, 1.19)

60,989,633
0.34
0.09
0.62
0.24

rs10499194
TNFAIP3
6
T
0.29
−0.05 (−0.26, 0.15)
−0.06 (−0.26, 0.12)
0.79 (0.50, 1.25)
0.86 (0.52, 1.40)

138,002,637
0.60
0.50
0.31
0.55

rs10919563
PTPRC
1 
A
0.11
−0.05 (−0.34, 0.23)
0.05 (−0.21, 0.32) 
1.21 (0.65, 2.22) 
1.61 (0.79, 3.26) 

196,967,065
0.72
0.69
0.55
0.19

Values were significant for .

The analyses of response were modeled for the minor alleles. Additive models were used taking the homozygote for the major allele as the reference variable.

383 patients were included in the primary analysis (the outcome was the absolute change in the disease activity score (DAS28) between the baseline and the 6-month visit).

208 patients were assessed in the secondary analyses; 119 were good responders and 89 nonresponders according to the EULAR response criteria.

Values were presented as regression coefficient (coef.), 95% confidence intervals (CI), and values () for linear regression analyses and as odds ratio (OR), 95% confidence intervals (CI), and values () for logistic regression analyses.

In multivariate models, covariates included were number of years completed at school, health assessment questionnaire, and disease activity score at baseline.

SNP: single nucleotide polymorphism; Chr: chromosome; MAF: minor allele frequency.

Univariate analysis looking for association between the other five risk alleles and anti-TNF treatment response demonstrated a worse response for minor (G) allele, rs3761847 SNP, in TRAF1/C5 region either measured by a change in DAS28 at six months (coefficient (coef.) −0.24; 95% confidence interval (CI) −0.43, −0.06; value of 0.009) or by the proportion of good responders versus nonresponders at six months (OR 0.61; CI 0.41, 0.92; value of 0.018). After adjusting for years of school completed, HAQ, and DAS28 at baseline, in two multivariate models, this association remained significant for a value < 0.05 (Table 4). In the univariate linear model, rs3761847 accounted for 1.75% of the change in DAS28 at six months. The multivariate model with the clinical covariates described above explained 23.9% of this change. When rs3761847 was added to the clinical model, the increased to 25.4%.

Table 4: Association of the rs3761847 single nucleotide polymorphism of TRAF1/C5 locus with the response to anti-TNF treatment.

SNP
Ch
Position Bp
Genotype
Count
MAF
Change in DAS
Absolute change in DAS
EULAR good response versus nonresponse

(good = 119, non = 89)

Linear regression models
Logistic regression models

Univariate
Multivariate
Univariate
Multivariate

rs3761847
9
122,730,060
11
167

1.95 (1.26)
Coef. −0.24
Coef. −0.23
OR 0.61
OR 0.58

12
165
0.35
1.82 (1.31)
CI −0.43, −0.06 
CI −0.40, −0.06 
CI 0.41, 0.92 
CI 0.37, 0.91 

22
51

1.38 (1.19)
0.009
0.009
0.018
0.019

Values were significant for .

The analyses of response were modeled for the minor (G) allele, which is the risk allele for RA. Additive models were used with the homozygote for the two major alleles as the reference variable.

383 patients were included in the primary analysis (the outcome was the absolute change in the disease activity score (DAS28) between the baseline and the 6 month visit).

208 patients were assessed in the secondary analyses, 119 were good responders and 89 non-responders according to the EULAR response criteria.

Values were presented as regression coefficient (coef.), 95% confidence intervals (CI) and -values () for linear regression analyses and as odds ratio (OR), 95% confidence intervals (CI) and -values () for logistic regression analyses.

In multivariate models covariates included were number of years completed at school, health assessment questionnaire (HAQ) and disease activity score at baseline.

SNP: single nucleotide polymorphism; Ch: chromosome; Genotype 1 = major allele, 2 = minor allele; MAF: minor allele frequency; DAS: disease activity score.

There were no significant relationships between the other SNP markers tested and response to therapy (Table 3). Similarly, SE was not associated with response to anti-TNF therapy in our population (coef. 0.21, value 0.21; OR 1.24, CI 0.56, 2.71).

In a secondary stratified analysis, the relationship between rs3761847 and treatment response in 278 ACPA positive and 105 ACPA negative patients was tested, but no association was detected (ACPA positive group, coef. −0.16, value 0.09 and ACPA negative group, coef. −0.41, value 0.04).

We replicated these findings pooling 265 Spanish RA patients and testing the association between the risk alleles and anti-TNF treatment response (results not shown).

4. Discussion

Our study did not replicate the association previously published between PTPRC rs10919563 variant and the response to anti-TNF therapy in patients with RA. The analyses suggest that in our Southern European population, the minor (G) allele rs3761847 in the TRAF1/C5 locus might have an association with poor response to anti-TNF treatment at six months. These results were consistent using either the absolute change in DAS28 or the proportion of good/non-responders as outcomes in univariate and multivariate models adjusted for clinical predictors of response. The other four RA susceptibility loci tested and the HLADRB1 were not associated with anti-TNF response.

In 2007, TRAF1/C5 was identified as a risk locus for RA by Plenge and colleagues in a GWAS [21] and by Kurreeman et al. in a candidate gene approach [38]. More recently, it was shown to be a marker of disease severity [10, 11] and in one study of noncardiovascular mortality [12]. No previous studies have reported the association with response to anti-TNF therapy. Due to a high level of linkage disequilibrium between the genes encoding TNF receptor associated factor 1 and complement component 5, it is currently not possible to assure which of these two genes at 9q33.2 encloses the causal variant. Both are possible candidates. The protein, TNF receptor-associated factor 1 (TRAF1), is a member of the TNF receptor (TNFR) associated factor (TRAF) protein family and is encoded by the TRAF gene. TRAF proteins associate with and mediate the signal transduction from various receptors of the TNFR superfamily. TRAF1 and TRAF2 form a heterodimeric complex, which is required for TNF mediated activation of MAPK8/JNK and NF-kB. The protein complex interacts with inhibitor-of-apoptosis proteins and mediates the antiapoptotic signals from TNF receptors. TRAF 1 is a negative regulator of TNF receptor and Toll-like receptor signaling and may contribute to the proliferation of T cells. rs3761847 is located at the upstream of TRAF1 and the downstream of the complement fraction C5 [39]. The clinical and biologic data for C5 are equally relevant. The complement pathway has been implicated in the pathogenesis of rheumatoid arthritis for more than 30 years. C5 cleavage generates the proinflammatory anaphylatoxin C5a, as well as C5b, which initiates the generation of the membrane-attack complex. C5-deficient mice are resistant to inflammatory arthritis in models with a dominant humoral component [21, 40]. It is compelling to hypothesize that this variant or other causative variants at this locus may influence the function or expression levels of TRAF1 and/or C5, affecting RA susceptibility, severity, and anti-TNF treatment response. Functional studies are warranted to confirm these hypotheses. Nevertheless, our result could be that false positive and studies with other populations are required to confirm the replication of these findings.

We did not find any significant associations between response to treatment and the presence of SE, neither with the other four RA risk allele variants related with NF-kB signaling pathway nor with citrullination, chosen for its high association with RA susceptibility. We were also not able to replicate the previously reported PTPRC association with treatment response in our Southern European population. Plant et al. reported a value of 0.04 for PTPRC association with anti-TNF treatment response and no significance for the ACPA positive patients in the stratified analysis. This result was strengthened with a meta-analysis combining their data with the Cui et al. study [29]. Although our sample size was large enough to detect allele variants association with response with a power > 80% for minor allele frequency (MAF) > 0.1, PTPRC showed a MAF of 0.11 which might have made difficult the detection of association. The different genetic background of our population could also account for the lack of replication. In a recent study, we and two Japanese groups also failed to replicate the results found in Northern European populations [41].

We are far from understanding the genetic mechanisms that underlie treatment response in patients with RA, as demonstrated by the large proportion of the variance explained by clinical factors compared to that explained by a single SNP. The ultimate goal for genetic, laboratory, and clinical predictors of treatment response studies is personalizing treatment and medicine practice by identifying biomarkers and specific phenotypes clinically useful for improving the therapeutic strategy. The identification of individual predictors may contribute to building complex algorithms aimed at improving the prediction of a better/worse response to anti-TNF drugs and other classes of biologic therapies.

5. Conclusions

We were not able to replicate the previously reported PTPRC rs10919563 association with treatment response in our Southern European RA population.

The minor (G) allele of rs3761847 in the TRAF1/C5 locus, which is a susceptibility factor for RA related to TNF signaling, was associated with a poor response to anti-TNF treatment at six months, using either the absolute change in DAS28 or the proportion of good responders and nonresponders as outcomes. This association was also observed after adjustment in a multivariate model with baseline clinical predictors of response.

We did not find any significant associations between response to treatment and the presence of SE, neither with the other four RA risk allele variants related with NF-kB signaling pathway nor with citrullination, chosen for its high association with susceptibility.

Additional studies in other populations are necessary to confirm the relevance of these findings.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by a grant from Harvard-Portugal Program HMSP-ICS/SAU-ICT/0002/2010. Daniel H. Solomon received support for this work from the NIH (K24-AR-055989). Elizabeth W. Karlson received support for this work from NIH (K24-AR-AR0524). Reuma.pt received unrestricted grants from Abbott, Bristol-Myers Squibb, Merck Sharp and Dohme, Pfizer, Roche, and UCB Pharma.

References

P. C. Taylor and M. Feldmann, “Anti-TNF biologic agents: still the therapy of choice for rheumatoid arthritis,” Nature Reviews Rheumatology, vol. 5, no. 10, pp. 578–582, 2009.

R. Maini, E. W. St Clair, F. Breedveld et al., “Infliximab (chimeric anti-tumour necrosis factor α monoclonal antibody) versus placebo in rheumatoid arthritis patients receiving concomitant methotrexate: a randomised phase III trial,” The Lancet, vol. 354, no. 9194, pp. 1932–1939, 1999.

I. B. McInnes and G. Schett, “The pathogenesis of rheumatoid arthritis,” The New England Journal of Medicine, vol. 365, no. 23, pp. 2205–2219, 2011.

M. L. Hetland, I. J. Christensen, U. Tarp et al., “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 and Rheumatism, vol. 62, no. 1, pp. 22–32, 2010.

K. L. Hyrich, K. D. Watson, A. J. Silman, and D. P. M. Symmons, “Predictors of response to anti-TNF-α therapy among patients with rheumatoid arthritis: results from the British Society for Rheumatology Biologics Register,” Rheumatology, vol. 45, no. 12, pp. 1558–1565, 2006.

H. Canhão, A. M. Rodrigues, A. F. Mourão et al., “Comparative effectiveness and predictors of response to tumour necrosis factor inhibitor therapies in rheumatoid arthritis,” Rheumatology, vol. 51, no. 11, pp. 2020–2026, 2012.

L. Mancarella, F. Bobbio-Pallavicini, F. Ceccarelli et al., “Good clinical response, remission, and predictors of remission in rheumatoid arthritis patients treated with tumor necrosis factor-α blockers: the GISEA study,” The Journal of Rheumatology, vol. 34, no. 8, pp. 1670–1673, 2007.

D. Ligeiro, J. E. Fonseca, O. Abade et al., “Influence of human leucocyte antigen-DRB1 on the susceptibility to rheumatoid arthritis and on the production of anti-cyclic citrullinated peptide antibodies in a Portuguese population,” Annals of the Rheumatic Diseases, vol. 66, no. 2, pp. 246–248, 2007.

T. Barnetche, A. Constantin, A. Cantagrel, A. Cambon-Thomsen, and P.-A. Gourraud, “New classification of HLA-DRB1 alleles in rheumatoid arthritis susceptibility: a combined analysis of worldwide samples,” Arthritis Research and Therapy, vol. 10, article R26, 2008.

A. W. Morgan, J. I. Robinson, P. G. Conaghan et al., “Evaluation of the rheumatoid arthritis susceptibility loci HLA-DRB1, PTPN22, OLIG3/TNFAIP3, STAT4 and TRAF1/C5 in an inception cohort,” Arthritis Research and Therapy, vol. 12, no. 2, article R57, 2010.

J. A. B. Van Nies, R. B. Marques, S. Trompet et al., “TRAF1/C5 polymorphism is not associated with increased mortality in rheumatoid arthritis: two large longitudinal studies,” Arthritis Research and Therapy, vol. 12, no. 2, article R38, 2010.

V. F. Panoulas, J. P. Smith, P. Nightingale, and G. D. Kitas, “Association of the TRAF1/C5 locus with increased mortality, particularly from malignancy or sepsis, in patients with rheumatoid arthritis,” Arthritis and Rheumatism, vol. 60, no. 1, pp. 39–46, 2009.

J. R. O'Dell, B. S. Nepom, C. Haire et al., “HLA-DRB1 typing in rheumatoid arthritis: Predicting response to specific treatments,” Annals of the Rheumatic Diseases, vol. 57, no. 4, pp. 209–213, 1998.

G. Saruhan-Direskeneli, M. Inanc, I. Fresko et al., “The role of HLA-DRB1 shared epitope alleles in predicting short-term response to leflunomide in rheumatoid arthritis,” Rheumatology, vol. 46, no. 12, pp. 1842–1844, 2007.

L. A. Criswell, R. F. Lum, K. N. Turner et al., “The influence of genetic variation in the HLA-DRB1 and LTA-TNF regions on the response to treatment of early rheumatoid arthritis with methotrexate or etanercept,” Arthritis & Rheumatism, vol. 50, no. 9, pp. 2750–2756, 2004.

A. Martinez, M. Salido, G. Bonilla et al., “Association of the major histocompatibility complex with response to infliximab therapy in rheumatoid arthritis patients,” Arthritis and Rheumatism, vol. 50, no. 4, pp. 1077–1082, 2004.

J. E. Fonseca, T. Carvalho, M. Cruz et al., “Polymorphism at position −308 of the tumour factor α gene and rheumatoid arthritis pharmacogenetics,” Annals of the Rheumatic Diseases, vol. 64, no. 5, pp. 793–794, 2005.

S. Pavy, E. J. M. Toonen, C. Miceli-Richard et al., “Tumour necrosis factor α-308GA polymorphism is not associated with response to TNFα blockers in Caucasian patients with rheumatoid arthritis: systematic review and meta-analysis,” Annals of the Rheumatic Diseases, vol. 69, no. 6, pp. 1022–1028, 2010.

M. Fabris, E. Di Poi, S. Sacco, G. Damante, L. Sinigaglia, and G. Ferraccioli, “TNF-α gene polymorphisms in rheumatoid arthritis patients treated with anti-TNF-α agents: preliminary results,” Reumatismo, vol. 54, no. 1, pp. 19–26, 2002.

L. Padyukov, J. Lampa, M. Heimbürger et al., “Genetic markers for the efficacy of tumour necrosis factor blocking therapy in rheumatoid arthritis,” Annals of the Rheumatic Diseases, vol. 62, no. 6, pp. 526–529, 2003.

R. M. Plenge, M. Seielstad, L. Padyukov et al., “TRAF1-C5 as a risk locus for rheumatoid arthritis—a genomewide study,” The New England Journal of Medicine, vol. 357, no. 12, pp. 1199–1209, 2007.

E. A. Stahl, S. Raychaudhuri, E. F. Remmers et al., “Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci,” Nature Genetics, vol. 42, pp. 509–514, 2010.

P. R. Burton, D. G. Clayton, L. R. Cardon et al., “Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls,” Nature, vol. 447, no. 7145, pp. 661–678, 2007.

R. M. Plenge, L. Padyukov, E. F. Remmers et al., “Replication of putative candidate-gene associations with rheumatoid arthritis in >4,000 samples from North America and Sweden: association of susceptibility with PTPN22, CTLA4, and PADI4,” The American Journal of Human Genetics, vol. 77, no. 6, pp. 1044–1060, 2005.

A. Barton, W. Thomson, X. Ke et al., “Re-evaluation of putative rheumatoid arthritis susceptibility genes in the post-genome wide association study era and hypothesis of a key pathway underlying susceptibility,” Human Molecular Genetics, vol. 17, no. 15, pp. 2274–2279, 2008.

W. Thomson, A. Barton, X. Ke et al., “Rheumatoid arthritis association at 6q23,” Nature Genetics, vol. 39, no. 12, pp. 1431–1433, 2007.

O. Snir, D. Gomez-Cabrero, A. Montes et al., “Non-HLA genes PTPN22, CDK6 and PADI4 are associated with specific autoantibodies in HLA-defined subgroups of rheumatoid arthritis,” Arthritis Research & Therapy, vol. 16, no. 5, article 414, 2014.

J. Varadé, R. Palomino-Morales, N. Ortego-Centeno et al., “Analysis of the rel polymorphism rs13031237 in autoimmune diseases,” Annals of the Rheumatic Diseases, vol. 70, no. 4, pp. 711–712, 2011.

J. Cui, S. Saevarsdottir, B. Thomson et al., “Rheumatoid arthritis risk allele PTPRC is also associated with response to anti-tumor necrosis factor α therapy,” Arthritis & Rheumatism, vol. 62, no. 7, pp. 1849–1861, 2010.

C. Potter, K. L. Hyrich, A. Tracey et al., “Association of rheumatoid factor and anti-cyclic citrullinated peptide positivity, but not carriage of shared epitope or PTPN22 susceptibility variants, with anti-tumour necrosis factor response in rheumatoid arthritis,” Annals of the Rheumatic Diseases, vol. 68, no. 1, pp. 69–74, 2009.

D. Plant, R. Prajapati, K. L. Hyrich et al., “Replication of association of the PTPRC gene with response to anti-tumor necrosis factor therapy in a large UK cohort,” Arthritis and Rheumatism, vol. 64, no. 3, pp. 665–670, 2012.

H. Canhão, A. Faustino, F. Martins et al., “Reuma.pt: the rheumatic diseases portuguese register,” Acta Reumatologica Portuguesa, vol. 36, pp. 45–56, 2011.

Biobanco-IMM, http://www.imm.fm.ul.pt/web/imm/biobank

J. E. Fonseca, M. Bernardes, H. Canhão et al., “Portuguese guide lines for the use of biological agents in rheumatoid arthritis—october 2011 update,” Acta Reumatologica Portuguesa, vol. 36, no. 4, pp. 385–388, 2011.

F. C. Arnett, S. M. Edworthy, D. A. Bloch et al., “The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis,” Arthritis and Rheumatism, vol. 31, no. 3, pp. 315–324, 1988.

DAS 28 calculator, /http://www.das-score.nl/das28/en

A. M. van Gestel, M. L. L. Prevoo, M. A. van’t Hof, M. H. van Rijswijk, L. B. A. van de Putte, and P. L. C. M. van Riel, “Development and validation of the European League against rheumatism response criteria for rheumatoid arthritis,” Arthritis & Rheumatism, vol. 39, no. 1, pp. 34–40, 1996.

F. A. Kurreeman, L. Padyukov, R. B. Marques et al., “A candidate gene approach identifies the TRAF1/C5 region as a risk factor for rheumatoid arthritis,” PLoS Medicine, vol. 4, no. 9, article e278, 2007.

M. Rothe, S. C. Wong, W. J. Henzel, and D. V. Goeddel, “A novel family of putative signal transducers associated with the cytoplasmic domain of the 75 kDa tumor necrosis factor receptor,” Cell, vol. 78, no. 4, pp. 681–692, 1994.

Y. Wang, J. Kristan, L. Hao, C. S. Lenkoski, Y. Shen, and L. A. Matis, “A role for complement in antibody-mediated inflammation: C5-deficient DBA/1 mice are resistant to collagen-induced arthritis,” Journal of Immunology, vol. 164, no. 8, pp. 4340–4347, 2000.

J. Cui, E. A. Stahl, S. Saevarsdottir et al., “Genome-wide association study and gene expression analysis identifies CD84 as a predictor of response to etanercept therapy in rheumatoid arthritis,” PLoS Genetics, vol. 9, no. 3, Article ID e1003394, 2013.