Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Polymorphisms in STAT4, PTPN2, PSORS1C1 and TRAF3IP2 Genes Are Associated with the Response to TNF Inhibitors in Patients with Rheumatoid Arthritis

  • Paola Conigliaro ,

    Contributed equally to this work with: Paola Conigliaro, Cinzia Ciccacci

    Affiliation Clinic of Rheumatology, Allergology and Clinical Immunology, Department of “Medicina dei Sistemi”, University of Rome Tor Vergata, Rome, Italy

  • Cinzia Ciccacci ,

    Contributed equally to this work with: Paola Conigliaro, Cinzia Ciccacci

    Affiliation Department of Biomedicine and Prevention, Genetics Section, University of Rome Tor Vergata, Rome, Italy

  • Cristina Politi,

    Affiliation Department of Biomedicine and Prevention, Genetics Section, University of Rome Tor Vergata, Rome, Italy

  • Paola Triggianese ,

    triggianese@med.uniroma2.it

    Affiliation Clinic of Rheumatology, Allergology and Clinical Immunology, Department of “Medicina dei Sistemi”, University of Rome Tor Vergata, Rome, Italy

  • Sara Rufini,

    Affiliation Department of Biomedicine and Prevention, Genetics Section, University of Rome Tor Vergata, Rome, Italy

  • Barbara Kroegler,

    Affiliation Clinic of Rheumatology, Allergology and Clinical Immunology, Department of “Medicina dei Sistemi”, University of Rome Tor Vergata, Rome, Italy

  • Carlo Perricone,

    Affiliation Reumatologia, Dipartimento di Medicina Interna e Specialità Mediche, Sapienza Università di Roma, Rome, Italy

  • Andrea Latini,

    Affiliation Department of Biomedicine and Prevention, Genetics Section, University of Rome Tor Vergata, Rome, Italy

  • Giuseppe Novelli,

    Affiliation Department of Biomedicine and Prevention, Genetics Section, University of Rome Tor Vergata, Rome, Italy

  • Paola Borgiani,

    Affiliation Department of Biomedicine and Prevention, Genetics Section, University of Rome Tor Vergata, Rome, Italy

  • Roberto Perricone

    Affiliation Clinic of Rheumatology, Allergology and Clinical Immunology, Department of “Medicina dei Sistemi”, University of Rome Tor Vergata, Rome, Italy

Abstract

Objective

Rheumatoid Arthritis (RA) is a progressive autoimmune disease characterized by chronic joint inflammation and structural damage. Remission or at least low disease activity (LDA) represent potentially desirable goals of RA treatment. Single nucleotide polymorphisms (SNPs) in several genes might be useful for prediction of response to therapy. We aimed at exploring 4 SNPs in candidate genes (STAT4, PTPN2, PSORS1C1 and TRAF3IP2) in order to investigate their potential role in the response to therapy with tumor necrosis factor inhibitors (TNF-i) in RA patients.

Methods

In 171 RA patients we investigated the following SNPs: rs7574865 (STAT4), rs2233945 (PSORS1C1), rs7234029 (PTPN2) and rs33980500 (TRAF3IP2). Remission, LDA, and EULAR response were registered at 6 months and 2 years after initiation of first line TNF-i [Adalimumab (ADA) and Etanercept (ETN)].

Results

STAT4 variant allele was associated with the absence of a good/moderate EULAR response at 2 years of treatment in the whole RA group and in ETN treated patients. The PTPN2 SNP was associated with no good/moderate EULAR response at 6 months in ADA treated patients. Patients carrying PSORS1C1 variant allele did not reach LDA at 6 months in both the whole RA group and ETN treated patients. TRAF3IP2 variant allele was associated with the lack of LDA and remission achievement at 6 months in all RA cohort while an association with no EULAR response at 2 years of treatment occurred only in ETN treated patients.

Conclusions

For the first time, we reported that SNPs in STAT4, PTPN2, PSORS1C1, and TRAF3IP2 are associated with response to TNF-i treatment in RA patients; however, these findings should be validated in a larger population.

Introduction

Rheumatoid Arthritis (RA) is a progressive autoimmune disease characterized by chronic joint inflammation and structural damage [1]. The management of RA has undergone significant changes with the current “treat to target” strategy [2]. The introduction of biological disease modifying anti-rheumatic drugs (bDMARDs) has changed the face of RA with remission or at least low disease activity (LDA) as achievable goals [3, 4]. Predictive biomarkers of response to therapy with bDMARDs could enable selection of the optimal treatment for the individual patients. Evidence assessed the value of age, gender, concomitant drugs, body mass index, or smoking status for predicting response to treatment [57]. Moreover, RA disease duration, disease activity, functional status, presence of autoantibodies [rheumatoid factor (RF) and anti-citrullinated peptide antibodies (ACPA)], and previous therapies can influence drug response [811]. Genetic inter-individual variability can also contribute to the differences in the response to treatment: some single nucleotide polymorphisms (SNPs) showed an association with bDMARDs response and might be useful for prediction, although few associations have been replicated [1215].

Some genes already known to be involved in RA susceptibility [16] could also be involved in the variability of the response to tumor necrosis factor (TNF)-inhibitors (TNF-i) drugs [15]. Among the known loci associated with RA, the signal transducer and activator of transcription 4 (STAT4) could be one of most interesting candidate genes to study in relation to drug response [17]. A recent meta-analysis demonstrated that SNPs in STAT4 confer susceptibility to RA in total subjects and in major ethnic groups. Moreover, this association was not dependent on RF and ACPA positivity [18]. The protein tyrosine phosphatase non-receptor 2 (PTPN2) is one of the newly investigated genes being recently reported as linked to the pathogenesis of RA [1921]. Additional genetic associations with the development of RA have been suggested for PSORS1C1/CDSN and TRAF3 Interacting Protein 2 (TRAF3IP2) genes that are well-known susceptibility genes for psoriasis and psoriatic arthritis [2123].

Thus, the aim of our study was to investigate the potential role of SNPs in STAT4, PTPN2, PSORS1C1, and TRAF3IP2 as predictors of remission and LDA in a cohort of RA patients treated with first line TNF-i.

Material and Methods

Patients

Medical records of RA patients referred to the Rheumatology Outpatient Clinic at the Department of “Medicina dei Sistemi” (“Policlinico Tor Vergata”, Rome, Italy), were retrospectively analyzed (time frame of the enrollment January 2008-December 2013). Patients were included in the study if they fulfilled the following inclusion criteria: the 2010 American College of Rheumatology (ACR)/European League Against Rheumatism (EULAR) classification criteria for RA [24], ≥ 18 years of age, inadequate response to at least one conventional synthetic (cs)DMARD, including Methotrexate, naïve for biologic treatment. Patients were excluded from the study if they showed impairment of hepatic/renal function, alcohol abuse, recent infection (with the last infection >3 month ago), ongoing history of malignancy (with interval malignancy-free >5 years) or ongoing pregnancy, and if they had missing or incomplete data in the follow-up visits. Therefore, the study included 171 RA patients of Caucasian origin. Patients received recommended doses of TNF-i: subcutaneous injection of Adalimumab (ADA) at 40 mg bi-weekly or Etanercept (ETN) at 50 mg every week. Disease activity and clinical response to therapy were assessed using Simplified Disease Activity Index (SDAI; LDA: ≤ 11, remission: ≤ 3.3) [25], disease activity score on 28 joints [DAS28 based on C-reactive protein (CRP)], and EULAR response criteria [25, 26]. The clinical and laboratory findings were evaluated at baseline and every 3 months from the start of TNF-i therapy; data of LDA, remission and EULAR response were registered at 6 months and 2 years after the beginning of the TNF-i treatment. Laboratory assessment included CRP, RF and ACPA. CRP and RF levels were assessed by nephelometry (normal range, 0–3 mg/L and 0–10 IU, respectively). ACPA were detected with a commercial third generation automated chemiluminescent kit: values >20 IU were considered positive. Peripheral blood samples were obtained at the time of the first medical evaluation from all included RA patients in order to perform the genetic analyses. All patients were naïve for biologic treatments at the time of blood sampling. Samples were stored at -80°C until they were analyzed. Written informed consent was obtained from patients. The study protocol was approved by the local ethics committee of the “Policlinico Tor Vergata” in Rome (Italy).

DNA extraction and genotyping

Genomic DNA was isolated from peripheral blood mononuclear cells using a Qiagen blood DNA mini kit. We have investigated the following SNPs, localized in the genes reported in parenthesis: rs7574865 (STAT4), rs7234029 (PTPN2), rs2233945 (PSORS1C1), and rs33980500 (TRAF3IP2). Genotyping was performed by allelic discrimination assay by TaqMan technology (Applied Biosystems, Foster City, CA, USA) and ABI PRISM 7000. Each assay was run including samples with known genotypes.

Statistical analysis

The Hardy–Weinberg equilibrium was verified for all SNPs by the Pearson χ2 test. We evaluated a possible correlation between the genetic variants and the SDAI LDA, SDAI remission and EULAR response, at 6 months and 2 years from the beginning of the TNF-i treatment. Differences in genotypes frequencies between groups of patients were evaluated by Pearson χ2 test or by the Fisher’s Exact test, where appropriate. Odds ratios (ORs) with 95% CI were calculated. A multivariate logistic regression analysis was used to correct the p-value for sex, csDMARDS and ACPA/RF positivity. All statistical analyses were performed by the SPSS program ver. 19 (IBM Corp, Armonk, NY USA). Two-tailed P values less than 0.05 were considered statistically significant.

Results

A total of 171 RA patients were included in the study, of whom 62.6% (n = 107) were treated with ETN and 37.4% (n = 64) were treated with ADA. Clinical and demographic data of the population are described in Table 1. Patients had longstanding disease in 72.5% of the cases. RF and ACPA were positive in 69.6% and 74.3% of patients, respectively. Mean SDAI at the beginning of the treatment was 27.6 ± 14. Patients with RA receiving concomitant csDMARDs comprised 77.2%. After 6 months of TNF-i treatment, SDAI remission was achieved in 26.7% of the whole RA population, SDAI-LDA was reached in 54% and a good-moderate EULAR response in 73.3% of patients. After 2 years of treatment, SDAI remission was achieved in 29.8%, SDAI-LDA in 63.5%, and a good-moderate EULAR response was reached in 77.9% of patients. At 2 years we observed 65 dropouts patients (38% of the whole study population) because of adverse events (n = 7, 10.7%), secondary failure (n = 33, 50.7%) and concomitant conditions (n = 25, 38.4%). No differences in demographic, clinical data and response to treatment were detected between subgroups of patients treated with ETN or ADA.

thumbnail
Table 1. Demographic and clinical data of 171 patients with Rheumatoid Arthritis included in the study.

https://doi.org/10.1371/journal.pone.0169956.t001

Associations of genetic variants with response to TNF-i treatment

We analyzed four SNPs in four candidate genes to investigate their possible role on TNF-i treatment response. In particular, we compared the genotypes distribution in relation to SDAI LDA (achieved vs not achieved), SDAI remission (achieved vs not achieved), and EULAR response (good/moderate vs no response). All the analyses were performed considering the clinical evaluations during follow up at 6 months and at 2 years after the treatment starting. Firstly, we performed a primary analysis considering the whole cohort of RA patients, independently from the specific administered TNF-i drug (Table 2). In accordance with this preliminary analysis the TRAF3IP2 SNP was associated with no achievement of LDA and remission at 6 months (P = 0.035 and OR = 0.36, P = 0.013 and OR = 0.11, respectively). These associations were confirmed after multiple correction for sex, csDMARDS and ACPA/RF positivity (Padj = 0.03 and Padj = 0.02 respectively). STAT4 SNP was associated with no EULAR response at 2 years of treatment (P = 0.05, OR = 0.38). Patients carrying PSORS1C1 variant allele did not reach LDA at 6 months (P = 0.002, OR = 0.35). The association with lack of achievement of LDA at 6 months was confirmed by the multiple correction (Padj = 0.003, ORadj = 0.36).

thumbnail
Table 2. Association analysis between TRAF3IP2, STAT4, PSORS1C1 and PTPN2 polymorphisms and response to TNF-inhibitors treatment in RA patients.

https://doi.org/10.1371/journal.pone.0169956.t002

In a second step, we repeated the analysis with each drug considered separately and we observed that the associations were drug-specific, except for the TRAF3IP2 SNP (Table 3 and Table 4). Indeed, TRAF3IP2 SNP showed a trend for association with lack of achievement of remission in ETN and ADA treated patients at 6 months of treatment after adjusting for sex, csDMARDS and ACPA/RF positivity (ETN: Padj = 0.08, ORadj = 0.14; ADA: Padj = 0.09, ORadj = 0.14).

thumbnail
Table 3. Association between analyzed polymorphisms and response to Etanercept treatment in RA patients.

https://doi.org/10.1371/journal.pone.0169956.t003

thumbnail
Table 4. Association between analysed polymorphisms and response to Adalimumab treatment in RA patients.

https://doi.org/10.1371/journal.pone.0169956.t004

With regard to ETN treatment, we confirmed an association between EULAR response and STAT4 SNP: patients carrying the variant allele showed a worse EULAR response at 2 years with P = 0.013 and OR = 0.16 (confirmed also after correction for sex, DMARDS and ACPA/RF positivity, Padj = 0.02 and ORadj = 0.14). The lack of achievement of LDA was also observed for patients carrying the variant allele of PSORS1C1 SNP (P = 0.012, OR = 0.35; after correction Padj = 0.023, ORadj = 0.37). Regarding the ADA treatment, we observed an association between EULAR response and TRAF3IP2 SNP at 2 years of treatment (P = 0.027, OR = 0.14): patients carrying the variant allele had a worse response to treatment. The PSORS1C1 SNP was associated with SDAI remission (P = 0.024) at 2 years of treatment, but the association was not confirmed after correction. Lastly, the PTPN2 SNP resulted associated with a worse EULAR response at 6 months of ADA treatment (P = 0.038, OR = 0.26), data partially confirmed also after multiple correction (Padj = 0.06, ORadj = 0.27).

The association of PSORS1C1 with a good-moderate EULAR response, detected in the preliminary analysis in the whole cohort, was not confirmed considering each drug separately. Nonetheless, it is worth to highlight that genotypes carrying the variant allele presented an higher frequency in patients with a good/moderate response at 2 years with both drugs (ORadj = 2.35 in the case of ETN; ORadj = 4.24 in the case of ADA), even if this difference did not reach the statistical significance (Padj = 0.21 and Padj = 0.23 respectively).

Fig 1 summarizes the odds ratio of the relevant associations observed between genetic variants and treatment outcome both in the whole RA population and in the two subgroups stratified by the drug (Fig 1).

thumbnail
Fig 1. Odds ratio for associations between genetic variants and response to TNF-i treatment.

(A) Associations in the whole cohort of Rheumatoid Arthritis patients at 6 months of treatment. [I] SDAI LDA (p 0.003); [II] SDAI remission (p 0.02); [III] SDAI LDA (p 0.03). (B) Associations in patients treated with Adalimumab (grey line) and Etanercept (black lines) considered separately. [I] EULAR response at 6 months of follow up (p 0.06); [II] SDAI LDA at 6 months of follow up (p 0.023,); [III] EULAR response at 2 years of follow up (p 0.02). Multivariate logistic regression analysis was used to correct the p-value for sex, csDMARDS and ACPA/RF positivity. Odds ratios (OR) with 95% CI were reported.

https://doi.org/10.1371/journal.pone.0169956.g001

Discussion

TNF-i have shown a good efficacy in the treatment of chronic inflammatory arthropathies, including RA [27]. However, a consistent part of patients does not reach the therapeutic targets of remission and LDA during TNF-i treatment. Evidence suggested the presence of a confined period of time, defined “window of opportunity”, in which RA patients are more susceptible to treatment [28, 29]. Moreover, a prolonged disease activity and a long disease duration at treatment initiation have been associated with unfavorable outcomes in RA [29, 30]. Furthermore, RA can cause progressive disability that leads to high direct and indirect costs for the health system [31]. Therefore, the identification of specific predicting factors of response to a specific treatment would be enormously useful in the clinical practice to select the patient that would benefit or not from the treatment [27, 32]. Candidate genes encoding proteins involved in the immune response have been fairly investigated to search for a possible association with TNF-i response in several autoimmune diseases, including RA [3336]. However, validated genomic biomarkers currently do not significantly allow the identification of non-responders before treatment in RA [15].

In this study, we evaluated the potential role of SNPs in STAT4, PTPN2, PSORS1C1, and TRAF3IP2 genes on the response to ETN and ADA treatment in RA patients. We selected these genes on the base of our previous studies that aimed to verify the association of common variants with different autoimmune diseases susceptibility [21, 3739]. In our previous study, STAT4, one of the most associated gene with RA susceptibility, related with a higher susceptibility to develop RA and with ACPA positivity, while SNPs in PSORS1C1 and PTPN2 genes were differently associated with joint damage in RA, even if we did not observe an association with RA susceptibility [21]. We considered as a candidate gene also TRAF3IP2 that encodes for Act1, an IL-17R adaptor protein sharing intracellular signal transduction molecules with the TNF-α signaling pathway, serving both as negative regulator of adaptive immunity and as a positive signaling adaptor in IL17-mediated immune responses [40, 41].

The results of the present study showed the presence of some associations of the investigated SNPs with remission, LDA, and EULAR response. Several lines of evidence in the literature demonstrated that ADA and ETN exert a different effect on innate and adaptive immune cell population in RA treated patients despite their similar effectiveness [4244]. Our preliminary data seems to suggest that the associations may be drug specific. However, given the small number of analyzed subjects, these results could be attributable to a low power of the study and need to be further investigated in larger samples. Regarding the response to ETN treatment, we observed a worse EULAR response in patients carrying the variant allele of STAT4 rs7574865 after 2 years of treatment and a lack of LDA achievement at 6 months of treatment in patients carrying the variant allele of PSORS1C1 rs2233945 SNP. These two associations were not present in the group of patients treated with ADA in which, on the contrary, we observed a worse EULAR response in patients carrying the variant allele of TRAF3IP2 after 2 years of treatment and the variant allele of PTPN2 after 6 months of treatment.

A recent study reported that the T allele of rs7574865 is significantly associated with higher levels of STAT4 mRNA and protein expression in a population of patients with early arthritis [45]. Since STAT4 is involved in the signaling of IL-12, IL-23, and IFN-γ, it has been suggested that patients carrying the rs7574865 minor allele might show stronger T helper (Th)1 and Th17 cytokine responses that are represented in RA [4648]. ETN seems able to downregulate both the Th1 and Th17 [49], therefore we could hypothesize that patients carrying the variant allele are less sensitive to the ETN effect. On the contrary, we could not demonstrate an association with ADA treatment that might be explained by previous findings suggesting that ADA seems to increase STAT4 activation in CD4+ T cells from RA patients [50].

We found that the variant allele of PTPN2 (rs7234029) was associated with a worse EULAR response at 6 months of ADA treatment. The same polymorphism, located in an intron region, was associated with a poor prognosis in a Portuguese population of RA patients treated with DMARDs and biologics [51]. PTPN2 was significantly overexpressed in synovial tissue samples from RA patients [52]. Interestingly, we observed a high frequency of bone erosions in patients carrying the variant allele [21]; this could be congruent with the worse response to treatment in patients carrying such variant. Although there are no functional studies on rs7234029 SNP, an in silico analysis revealed that it modulates potentially the binding sites of several transcription factors involved in inflammation [53]. Therefore, further analysis should firstly replicate our preliminary data on large cohort of patients and then define the potential effect of this SNP on different drug’s mechanism of action.

PSORS1C1 can affect IL-17 secretion that plays important roles in synovial inflammation and bone destruction in RA [54]. A significant increase in expression of PSORS1C1 in RA synovial tissues has been also described [55]. The functional role of rs2233945 is not known; the SNP is located in an intronic region and therefore we can speculate that it could be involved in expression regulation. Patients carrying the variant allele of this SNP treated with ETN are less likely to achieve LDA at 6 months suggesting that this genetic variation may interfere with the drug mechanism of action possibly increasing the IL-17 response. However, an opposite effect was suggested at 2 years of treatment in the whole RA population and in particular in the ADA group of patients. This contradictory result may be attributed to either a different mechanism of action of ADA in those patients carrying the variant allele, or to a bias related to the small number of analyzed patients and the dropouts at 2 years. Functional studies on this SNP will give important information on its putative role in TNF-i treated patients.

Regarding TRAF3IP2, we observed that remission and LDA at 6 months of treatment were not achieved in patients carrying the variant genotype, considering patients treated with ETN and ADA together, but this result was not replicated in the two subgroups stratified by the drug. TRAF3IP2 product interacts with TRAF proteins; in the TRAF gene family, TRAF1 is a negative regulator of TNF receptor that was identified as a risk locus for RA in a GWAS [56]. The rs33980500 SNP decreases the binding with TRAF6 and this could also alter TRAF2 and TRAF5 protein interactions within the IL-17R signaling pathway, leading to increased neutrophil chemotaxis and an enhanced immune response [57]. Recently, TRAF1 has been investigated in relation to TNF-i response by Canhao et al., who reported in Portuguese RA patients an association of the minor (G) allele of rs3761847 in the TRAF1/C5 with a poor response to TNF-i treatment at 6 months [58]. However, another study in a Greek population did not confirm the association of TRAF1 with TNF-i treatment [59].

This study presents some limitations. First, we have analyzed only one SNP for each selected gene. Indeed, it is possible that other SNPs (or a combination of SNPs) in these genes (and in other genes) could play a role in response to TNF-i treatment. Moreover, the small number of investigated patients may represent a bias that can confound the interpretation of the results, especially in the analysis concerning the two subgroups stratified by the drug. Given the limited number of patients, we have not performed the multiple comparisons correction, and therefore our results should be considered as preliminary data. Further studies are necessary to replicate our findings in a larger cohort of RA patients together with functional studies in order to confirm and better explore the contribution of these SNPs in the treatment response.

Author Contributions

  1. Conceptualization: RP PC CC C. Perricone.
  2. Data curation: PC BK CC.
  3. Formal analysis: CC SR C. Perricone.
  4. Investigation: PC PT CC C. Perricone.
  5. Methodology: PC PT CC SR AL C. Politi.
  6. Project administration: RP PB.
  7. Resources: PC PT BK RP GN PB.
  8. Software: CC C. Perricone.
  9. Supervision: RP PB.
  10. Validation: PC CC.
  11. Visualization: PC PT CC.
  12. Writing – original draft: PC PT CC C. Politi.
  13. Writing – review & editing: PC PT RP CC PB GN.

References

  1. 1. Combe B. (2009) Progression in early rheumatoid arthritis. Best Pract Res Clin Rheumatol 23: 59–69. pmid:19233046
  2. 2. Smolen JS, Breedveld FC, Burmester GR, Bykerk V, Dougados M, Emery P et al. (2016) Treating rheumatoid arthritis to target: 2014 update of the recommendations of an international task force. Ann Rheum Dis 75: 3–15. pmid:25969430
  3. 3. Smolen JS, Aletaha D, McInnes IB. (2016) Rheumatoid arthritis. Lancet. 388(10055):2023–2038. pmid:27156434
  4. 4. Conigliaro P, Chimenti MS, Triggianese P, Ballanti E, Sunzini F, Duca I, Perricone R. (2016) Remission and low disease activity in a cohort of real-life patients with rheumatoid arthritis treated with first-line anti-tumour necrosis factor. J Int Med Res 44(1 suppl): 90–94. pmid:27683148
  5. 5. Daïen CI, Morel J. (2014) Predictive factors of response to biological disease modifying antirheumatic drugs: towards personalized medicine. Mediators Inflamm 2014:386148. pmid:24523570
  6. 6. Liu Y, Hazlewood GS, Kaplan GG, Eksteen B, Barnabe C. (2016) The Impact of Obesity on Remission and Disease Activity in Rheumatoid Arthritis: A Systematic Review and Meta-Analysis. Arthritis Care Res (Hoboken)
  7. 7. Quintana-Duque MA, Rondon-Herrera F, Mantilla RD, Calvo-Paramo E, Yunis JJ, Varela-Nariño A, et al. (2016) Predictors of remission, erosive disease and radiographic progression in a Colombian cohort of early onset rheumatoid arthritis: a 3-year follow-up study. Clin Rheumatol 35: 1463–1473. pmid:27041382
  8. 8. Gottenberg JE, Ravaud P, Cantagrel A, Combe B, Flipo RM, Schaeverbeke T et al. (2012) Positivity for anti-cyclic citrullinated peptide is associated with a better response to abatacept: data from the 'Orencia and Rheumatoid Arthritis' registry. Ann Rheum Dis 71: 1815–1819. pmid:22615458
  9. 9. Scrivo R, Conigliaro P, Riccieri V, Di Franco M, Alessandri C, Spadaro A, et al. (2015) Distribution of IL-10 family cytokines in serum and synovial fluid of patients with inflammatory arthritis reveals different contribution to systemic and joint inflammation. Clin Exp Immunol 179: 300–308. pmid:25178435
  10. 10. Katchamart W, Johnson S, Lin HJ, Phumethum V, Salliot C, Bombardier C. (2010) Predictors for remission in rheumatoid arthritis patients: A systematic review. Arthritis Care Res (Hoboken) 62: 1128–1143.
  11. 11. Alessandri C, Conti F, Conigliaro P, Mancini R, Massaro L, Valesini G. (2009) Seronegative Autoimmune Diseases. Ann N Y Acad Sci 1173: 52–59. pmid:19758132
  12. 12. Plant D, Prajapati R, Hyrich KL, Morgan AW, Wilson AG, Isaacs JD, et al. (2012) Replication of association of the PTPRC gene with response to anti-tumor necrosis factor therapy in a large UK cohort. Arthritis Rheum 64: 665–670. pmid:21952740
  13. 13. Acosta-Colman I, Palau NF, Tornero JF, Fernandez-Nebro AF, Blanco FF, Gonzalez-Alvaro IF, et al. (2013) GWAS replication study confirms the association of PDE3A-SLCO1C1 with anti-TNF therapy response in rheumatoid arthritis. Pharmacogenomics 14: 727–734. pmid:23651021
  14. 14. Sode J, Vogel U, Bank S, Andersen PS, Hetland ML, Locht H, et al. (2015) Genetic Variations in Pattern Recognition Receptor Loci Are Associated with Anti-TNF Response in Patients with Rheumatoid Arthritis. PLoS One 10: e0139781. pmid:26440629
  15. 15. Sieberts SK, Zhu F, García-García J, Stahl E, Pratap A, Pandey G, et al. (2016) Crowdsourced assessment of common genetic contribution to predicting anti-TNF treatment response in rheumatoid arthritis. Nat Commun 7: 12460. pmid:27549343
  16. 16. Nabi G, Akhter N, Wahid M, Bhatia K, Mandal RK, Dar SA, et al. (2016) Meta-analysis reveals PTPN22 1858C/T polymorphism confers susceptibility to rheumatoid arthritis in Caucasian but not in Asian population. Autoimmunity 49: 197–210. pmid:26763276
  17. 17. Scott IC, Rijsdijk F, Walker J, Quist J, Spain SL, Tan R, et al. (2015) Do Genetic Susceptibility Variants Associate with Disease Severity in Early Active Rheumatoid Arthritis? J Rheumatol 42: 1131–1140. pmid:25979711
  18. 18. Elshazli R, Settin A. (2015) Association of PTPN22 rs2476601 and STAT4 rs7574865 polymorphisms with rheumatoid arthritis: A meta-analysis update. Immunobiology 220: 1012–1024. pmid:25963842
  19. 19. Aradi B, Kato M, Filkova M, Karouzakis E, Klein K, Scharl M, et al. (2015) Protein tyrosine phosphatase nonreceptor type 2: an important regulator of lnterleukin-6 production in rheumatoid arthritis synovial fibroblasts. Arthritis Rheumatol 67: 2624–2633. pmid:26139109
  20. 20. Solus JF, Chung CP, Oeser A, Li C, Rho YH, Bradley KM, et al. (2015) Genetics of serum concentration of IL-6 and TNFα in systemic lupus erythematosus and rheumatoid arthritis: a candidate gene analysis. Clin Rheumatol 34: 1375–1382. pmid:25652333
  21. 21. Ciccacci C, Conigliaro P, Perricone C, Rufini S, Triggianese P, Politi C, et al. (2016) Polymorphisms in STAT4, IL10, PSORS1C1, PTPN2 and MIR146A genes are differently associated with prognostic factors in Italian patients affected by Rheumatoid Arthritis. Clin Exp Immunol 186: 157–163. pmid:27342690
  22. 22. Sun H, Xia Y, Wang L, Wang Y, Chang X. (2013) PSORS1C1 may be involved in rheumatoid arthritis. Immunol Lett 153: 9–14. pmid:23769905
  23. 23. Stuart PE, Nair RP, Tsoi LC, Tejasvi T, Das S, Kang HM, et al. (2015) Genome-wide Association Analysis of Psoriatic Arthritis and Cutaneous Psoriasis Reveals Differences in Their Genetic Architecture. Am J Hum Genet 97: 816–836. pmid:26626624
  24. 24. Aletaha D, Neogi T, Silman AJ, Funovits J, Felson DT, Bingham CO III, et al. (2010) 2010 rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum 62: 2569–2581. pmid:20872595
  25. 25. Smolen JS, Breedveld FC, Schiff MH, Kalden JR, Emery P, Eberl G, et al. (2003) A simplified disease activity index for rheumatoid arthritis for use in clinical practice. Rheumatology 42: 244–257. pmid:12595618
  26. 26. van Gestel AM, Prevoo ML, van 't Hof MA, van Rijswijk MH, van de Putte LB, van Riel PL. (1996) 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 39: 34–40. pmid:8546736
  27. 27. Cuppen BV, Welsing PM, Sprengers JJ, Bijlsma JW, Marijnissen AC, van Laar JM, et al. (2016) Personalized biological treatment for rheumatoid arthritis: a systematic review with a focus on clinical applicability. Rheumatology 55: 826–839. pmid:26715775
  28. 28. Cush JJ. (2007) Early rheumatoid arthritis—is there a window of opportunity? J Rheumatol Suppl. 80:1–7. pmid:17985417
  29. 29. van Nies JA, Tsonaka R, Gaujoux-Viala C, Fautrel B, van der Helm-van Mil AH. (2015) Evaluating relationships between symptom duration and persistence of rheumatoid arthritis: does a window of opportunity exist? Results on the Leiden early arthritis clinic and ESPOIR cohorts. Ann Rheum Dis 74: 806–812. pmid:25561360
  30. 30. Chimenti MS, Triggianese P, Conigliaro P, Candi E, Melino G, Perricone R. (2015) The interplay between inflammation and metabolism in rheumatoid arthritis. Cell Death Dis 6:e1887. pmid:26379192
  31. 31. Wallman JK, Eriksson JK, Nilsson JÅ, Olofsson T, Kristensen LE, Neovius M, et al. (2016) Costs in Relation to Disability, Disease Activity, and Health-related Quality of Life in Rheumatoid Arthritis: Observational Data from Southern Sweden. J Rheumatol 43: 1292–1299. pmid:27252420
  32. 32. Conigliaro P, Chimenti MS, Triggianese P, Sunzini F, Novelli L, Perricone C, et al. (2016) Autoantibodies in inflammatory arthritis. Autoimmun Rev 15: 673–683. pmid:26970491
  33. 33. Tong Q, Zhao L, Qian XD, Zhang LL, Xu X, Dai SM, et al. (2013) Association of TNF-α polymorphism with prediction of response to TNF blockers in spondyloarthritis and inflammatory bowel disease: a meta-analysis. Pharmacogenomics 14: 1691–700. pmid:24192118
  34. 34. Prieto-Pérez R, Almoguera B, Cabaleiro T, Hakonarson H, Abad-Santos F. (2016) Association between Genetic Polymorphisms and Response to Anti-TNFs in Patients with Inflammatory Bowel Disease. Int J Mol Sci 17: 225. pmid:26861312
  35. 35. Oliver J, Plant D, Webster AP, Barton A. (2015) Genetic and genomic markers of anti-TNF treatment response in rheumatoid arthritis. Biomark Med 9: 499–512. pmid:26079957
  36. 36. Morales-Lara MJ, Cañete JD, Torres-Moreno D, Hernández MV, Pedrero F, Celis R, et al. (2012) Effects of polymorphisms in TRAILR1 and TNFR1A on the response to anti-TNF therapies in patients with rheumatoid and psoriatic arthritis. Joint Bone Spine 79: 591–596. pmid:22480748
  37. 37. Ciccacci C, Biancone L, Di Fusco D, Ranieri M, Condino G, Giardina E, et al. (2013) TRAF3IP2 gene is associated with cutaneous extraintestinal manifestations in inflammatory bowel disease. J Crohns Colitis 7: 44–52. pmid:22445837
  38. 38. Perricone C, Ciccacci C, Ceccarelli F, Di Fusco D, Spinelli FR, Cipriano E, et al. (2013) TRAF3IP2 gene and systemic lupus erythematosus: association with disease susceptibility and pericarditis development. Immunogenetics 65: 703–709. pmid:23836313
  39. 39. Ciccacci C, Perricone C, Ceccarelli F, Rufini S, Di Fusco D, Alessandri C, et al. (2014) A multilocus genetic study in a cohort of Italian SLE patients confirms the association with STAT4 gene and describes a new association with HCP5 gene. PLoS One 9:e111991. pmid:25369137
  40. 40. Iwakura Y, Ishigame H, Saijo S, Nakae S. (2011) Functional specialization of interleukin-17 family members. Immunity 34: 149–162. pmid:21349428
  41. 41. Urabe S, Isomoto H, Ishida T, Maeda K, Inamine T, Kondo S, et al. (2015) Genetic Polymorphisms of IL-17F and TRAF3IP2 Could Be Predictive Factors of the Long-Term Effect of Infliximab against Crohn's Disease. Biomed Res Int 2015: 416838. pmid:26558270
  42. 42. Conigliaro P, Triggianese P, Perricone C, Chimenti MS, Di Muzio G, Ballanti E, et al. (2014) Restoration of peripheral blood natural killer and B cell levels in patients affected by rheumatoid and psoriatic arthritis during etanercept treatment. Clin Exp Immunol 177: 234–243. pmid:24666401
  43. 43. Furst DE, Wallis R, Broder M, Beenhouwer DO. (2006) Tumor necrosis factor antagonists: different kinetics and/or mechanisms of action may explain differences in the risk for developing granulomatous infection. Semin Arthritis Rheum 36: 159–167. pmid:16884970
  44. 44. Prajapati R, Plant D, Barton A. (2011) Genetic and genomic predictors of anti-TNF response. Pharmacogenomics 12: 1571–1585. pmid:22044414
  45. 45. Lamana A, López-Santalla M, Castillo-González R, Ortiz AM, Martín J, García-Vicuña R, González-Álvaro I. (2015) The Minor Allele of rs7574865 in the STAT4 Gene Is Associated with Increased mRNA and Protein Expression. PLoS One 10:e0142683. pmid:26569609
  46. 46. Korman BD, Kastner DL, Gregersen PK, Remmers EF. (2008) STAT4: genetics, mechanisms, and implications for autoimmunity. Curr allergy asthma rep 8: 398–403. pmid:18682104
  47. 47. Frucht DM, Aringer M, Galon J, Danning C, Brown M, Fan S, et al. (2000) STAT4 is expressed in activated peripheral blood monocytes, dendritic cells, and macrophages at sites of Th1-mediated inflammation. J Immunol 164: 4659–4664. pmid:10779770
  48. 48. Patakas A, Benson RA, Withers DR, Conigliaro P, McInnes IB, Brewer JM, Garside P. (2012) Th17 effector cells support B cell responses outside of germinal centres. PLos One. 7: e49715. pmid:23166752
  49. 49. Antiga E, Volpi W, Cardilicchia E, Maggi L, Filì L, Manuelli C, et al. (2012) Etanercept downregulates the Th17 pathway and decreases the IL-17+/IL-10+ cell ratio in patients with psoriasis vulgaris. J Clin Immunol 32: 1221–1232. pmid:22699761
  50. 50. Aerts NE, Ebo DG, Bridts CH, Stevens WJ, De Clerck LS. (2010) T cell signal transducer and activator of transcription (STAT) 4 and 6 are affected by adalimumab therapy in rheumatoid arthritis. Clin Exp Rheumatol 28: 208–214. pmid:20483042
  51. 51. Mourão AF, Santos MJ, Mendonça S, Oliveira-Ramos F, Salgado M, Estanqueiro P, et al. (2015) Genetic Predictors of Poor Prognosis in Portuguese Patients with Juvenile Idiopathic Arthritis: Data from Reuma.pt. J Immunol Res 2015: 706515. pmid:26504858
  52. 52. Aradi B, Kato M, Filkova M, Karouzakis E, Klein K, Scharl M, et al. (2015) Protein tyrosine phosphatase nonreceptor type 2: an important regulator of lnterleukin-6 production in rheumatoid arthritis synovial fibroblasts. Arthritis Rheumatol 67: 2624–2633. pmid:26139109
  53. 53. Glas J, Wagner J, Seiderer J, Olszak T, Wetzke M, Beigel F, et al. (2012) PTPN2 gene variants are associated with susceptibility to both Crohn's disease and ulcerative colitis supporting a common genetic disease background. PLoS One 7: e33682. pmid:22457781
  54. 54. Moon YM, Yoon BY, Her YM, Oh HJ, Lee JS, Kim KW, et al. (2012) IL-32 and IL-17 interactand have the potential to aggravate osteoclastogenesis in rheumatoid arthritis. Arthritis Res Ther 14: R246. pmid:23148681
  55. 55. Sun H, Xia Y, Wang L, Wang Y, Chang X. (2013) PSORS1C1 may be involved in rheumatoid arthritis. Immunol Lett 153: 9–14. pmid:23769905
  56. 56. Plenge RM, Seielstad M, Padyukov L, Lee AT, Remmers EF, Ding B, et al. (2007) TRAF1-C5 as a risk locus for rheumatoid arthritis—a genomewide study. N Engl J Med 357: 1199–1209. pmid:17804836
  57. 57. Doyle MS, Collins ES, FitzGerald OM, Pennington SR. (2012) New insight into the functions of the interleukin-17 receptor adaptor protein Act1 in psoriatic arthritis. Arthritis Res Ther 14: 226. pmid:23116200
  58. 58. Canhão H, Rodrigues AM, Santos MJ, Carmona-Fernandes D, Bettencourt BF, Cui J, et al. (2015) TRAF1/C5 but not PTPRC variants are potential predictors of rheumatoid arthritis response to anti-tumor necrosis factor therapy. Biomed Res Int 2015: 490295. pmid:25834819
  59. 59. Zervou MI, Myrthianou E, Flouri I, Plant D, Chlouverakis G, Castro-Giner F, et al. (2013) Lack of association of variants previously associated with anti-TNF medication response in rheumatoid arthritis patients: results from a homogeneous Greek population. PLoS One 8: e74375. pmid:24040234