The TT Genotype of the STAT4 rs7574865 Polymorphism Is Associated with High Disease Activity and Disability in Patients with Early Arthritis

Abstract

Background

The number of copies of the HLA-DRB1 shared epitope, and the minor alleles of the STAT4 rs7574865 and the PTPN22 rs2476601 polymorphisms have all been linked with an increased risk of developing rheumatoid arthritis. In the present study, we investigated the effects of these genetic variants on disease activity and disability in patients with early arthritis.

Methodology and Results

We studied 640 patients with early arthritis (76% women; median age, 52 years), recording disease-related variables every 6 months during a 2-year follow-up. HLA-DRB1 alleles were determined by PCR-SSO, while rs7574865 and rs2476601 were genotyped with the Taqman 5′ allelic discrimination assay. Multivariate analysis was performed using generalized estimating equations for repeated measures. After adjusting for confounding variables such as gender, age and ACPA, the TT genotype of rs7574865 in STAT4 was associated with increased disease activity (DAS28) as compared with the GG genotype (β coefficient [95% confidence interval] = 0.42 [0.01–0.83], p = 0.044). Conversely, the presence of the T allele of rs2476601 in PTPN22 was associated with diminished disease activity during follow-up in a dose-dependent manner (CT genotype = −0.27 [−0.56– −0.01], p = 0.042; TT genotype = −0.68 [−1.64– −0.27], p = 0.162). After adjustment for gender, age and disease activity, homozygosity for the T allele of rs7574865 in STAT4 was associated with greater disability as compared with the GG genotype.

Conclusions

Our data suggest that patients with early arthritis who are homozygous for the T allele of rs7574865 in STAT4 may develop a more severe form of the disease with increased disease activity and disability.

About the Authors:

Amalia Lamana

Affiliation: Rheumatology Service, Hospital Universitario de La Princesa, Instituto de Investigación Sanitaria Princesa, Madrid, Spain

Alejandro Balsa

Affiliation: Rheumatology Service, Hospital Universitario La Paz, IdiPaz, Instituto de Investigación Sanitaria La Paz, Madrid, Spain

Blanca Rueda

Affiliation: Departamento de Enfermería, Facultad de Ciencias de la Salud, Universidad de Granada, Granada, Spain

Ana M. Ortiz

Affiliation: Rheumatology Service, Hospital Universitario de La Princesa, Instituto de Investigación Sanitaria Princesa, Madrid, Spain

Laura Nuño

Affiliation: Rheumatology Service, Hospital Universitario La Paz, IdiPaz, Instituto de Investigación Sanitaria La Paz, Madrid, Spain

Maria Eugenia Miranda-Carus

Affiliation: Rheumatology Service, Hospital Universitario La Paz, IdiPaz, Instituto de Investigación Sanitaria La Paz, Madrid, Spain

Maria F. Gonzalez-Escribano

Affiliation: Immunology Service, Hospital Universitario Virgen del Rocio, Sevilla, Spain

Miguel A. Lopez-Nevot

Affiliation: Immunology Service, Hospital Universitario Virgen de las Nieves, Granada, Spain

Dora Pascual-Salcedo

Affiliation: Immunology Service, Hospital Universitario La Paz, IdiPaz, Instituto de Investigación Sanitaria La Paz, Madrid, Spain

Javier Martin

Affiliation: Instituto de Parasitología y Biomedicina Lopez Neyra, CSIC, Granada, Spain

Isidoro González-Álvaro

* E-mail: [email protected]

Affiliation: Rheumatology Service, Hospital Universitario de La Princesa, Instituto de Investigación Sanitaria Princesa, Madrid, Spain

Introduction

Rheumatoid arthritis (RA) is caused by an interaction between genetic and environmental factors that leads to a break in immunological self-tolerance. Twin concordance data suggest that RA has a heritability rate of 50–60% [1]. The HLA region is the main genetic factor accounting for this heritability, particularly the HLA-DRB1 alleles that encode the shared epitope (SE) [2]. In addition, during the past few years, genome-wide association studies and candidate gene association studies have led to the identification of new genes associated with the risk of developing RA [3]. The role of PTPN22 and STAT4 as genetic markers of susceptibility to RA has been extensively replicated in different populations [3]–[5].

In contrast to the clear association between the risk of developing RA and alleles in HLA-DRB1, PTPN22, and STAT4, it remains unclear whether these genetic factors affect disease severity once immunological self-tolerance is broken. Some authors have reported an association between the SE and radiological progression, systemic complications and cardiovascular events in RA [6]–[9]. However, several studies of the role of genetic factors outside the HLA region in the radiological progression of RA have produced contradictory results [10], and little information is available on the impact of genetic factors on disease activity and disability [11].

When our study was designed, only the presence of the SE and the polymorphisms in PTPN22 and STAT4 were clearly associated with an increased risk of developing RA in distinct Spanish populations. Thus, the objective of our study was to determine whether the presence of these genetic markers contributes to a poorer disease course and increased disability in patients with early arthritis (EA).

[Figure omitted. See PDF.]

Table 1. Baseline characteristics of patients with early arthritis.

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

Methods

Ethics Statement

The local ethics committees reviewed and approved the protocols of both EA Clinical Registers, and all the patients signed an informed consent form prior to inclusion in the study.

Objective

The hypothesis of the present study is that the genetic variants that confer a high risk of developing RA are also associated with a worse clinical course of disease. The specific objective was to determine whether the presence of the SE, the minor alleles of STAT4 rs7574865 or PTPN22 rs2476601 polymorphisms affect disease activity and disability in EA patients.

[Figure omitted. See PDF.]

Figure 1. Patients with early arthritis who were homozygous for the minor rs7574865 allele of STAT4 exhibited stronger disease activity during the follow-up.

The data are categorized according to the rs7574865 genotype of STAT4 (upper panel), the rs2476601 genotype of PTPN22 (middle panel), and the number of copies of the shared epitope (SE; lower panel). Data represent the interquartile range (p75, upper edge of the box; p25, lower edge; p50, midline), p95 (line above the box) and p5 (line below the box) of the DAS28 adjusted for gender and age (aDAS28). Dots represent outliers.

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

[Figure omitted. See PDF.]

Table 2. Effect of the presence of the shared epitope and the minor rs7574865 allele in STAT4 or rs2476601 in PTPN22 on disease activity (DAS28) and disability (HAQ) in patients with early arthritis.

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

Participants

The study sample comprised 641 patients with EA who attended the Hospital Universitario (HU) La Paz (355 patients) and HU La Princesa (286 patients). Inclusion criteria included 2 or more swollen joints for at least 4 weeks and symptoms for less than a year. Only data from patients fulfilling the 1987 American College of Rheumatology criteria for RA [12] within the 2-year follow-up (n = 464) or with chronic undifferentiated arthritis (UA: n = 177) were analyzed. The protocol included 5 visits during the two-year follow-up period and at each visit, the following data were collected and entered into an electronic database: clinical and demographic information; disease duration at the beginning of follow-up; Disease Activity Score of 28 tender and swollen joint counts (DAS28) [13]; global disease activity on a 100-mm visual analogue scale assessed by both the patient and the physician; the score from the Health Assessment Questionnaire (HAQ; Spanish version) [14]; and the results of laboratory tests including the erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), rheumatoid factor (RF, assessed by nephelometry; positive >20 IU/ml), and anti-citrullinated peptide antibody (ACPA, measured by enzyme immunoassay: Euro-Diagnostica Immunoscan RA; positive >50 IU/ml).

[Figure omitted. See PDF.]

Figure 2. Effect of the rs7574865 and rs2476601 polymorphisms in STAT4 and PTPN22 on the level of disease activity in patients with early arthritis.

A) STAT4 genotypes: GG, white bars; GT, gray bars; TT, black bars. B) PTPN22 genotypes: CC, white bars; CT, gray bars; TT, black bars. Data represent the proportion of visits in which remission, or low, moderate or high disease activity was recorded in the patients of each genotype. This figure illustrates the findings presented in Table 3, in which the statistical significance of the differences in disease activity associated with the genotype of each polymorphism is shown.

https://doi.org/10.1371/journal.pone.0043661.g002

[Figure omitted. See PDF.]

Table 3. Effect of the presence of the shared epitope and the minor rs7574865 allele of STAT4 or rs2476601 in PTPN22 on the level of disease activity.

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

Description of the Procedures

Genotyping.

DNA was obtained from peripheral blood using standard methods, and the samples were genotyped for the PTPN22 rs2476601 and the STAT4 rs7574865 genetic variants using a predesigned SNP Genotyping Assay (Applied Biosystems, Part numbers: C_16021387_20 and C_29882391_10, respectively: Foster City, California, USA). Polymerase chain reactions (PCRs) were carried out according to the manufacturer’s instructions and after PCR, the genotype of each sample was automatically attributed by measuring allele-specific fluorescence on an ABI Prism 7900 Sequence Detection System using SDS version 2.3 (Applied Biosystems). Duplicate samples and negative controls were included to verify the genotyping accuracy.

[Figure omitted. See PDF.]

Figure 3. Patients with early arthritis who are homozygous for the minor rs7574865 allele of STAT4 exhibit increased disability during follow-up.

The data are categorized according to the rs7574865 genotype of STAT4 (upper panel), the rs2476601 genotype of PTPN22 (middle panel), and the number of copies of the shared epitope (SE; lower panel). The data are presented as the interquartile range (p75, upper edge of the box; p25, lower edge; p50, midline), the p95 (line above box) and p5 (line below the box) of the HAQ adjusted for gender and age (aHAQ). Dots represent outliers.

https://doi.org/10.1371/journal.pone.0043661.g003

HLA-DRB1 was genotyped using a reverse dot-blot kit with sequence-specific oligonucleotide probes (Dynal RELI™ SSO HLA-DRB1 typing kit: Dynal Biotech, Bromborough, UK). When necessary, high-resolution typing of HLA-DRB1_03 samples was performed using Dynal AllSet™ SSP DRB1_03.

Statistical Analysis

Normally distributed quantitative variables were represented as the mean (± standard deviation: SD), while non-normally distributed variables were represented as the median and interquartile range (IQR). Qualitative variables were described using an estimation of the proportions. Variables with a normal distribution were analyzed by the t-test, while the Mann Whitney or Kruskal-Wallis tests were used for variables with a non-normal distribution. A χ² or Fisheŕs exact test was used to compare categorical variables.

Only patients whose DAS28 or HAQ scores were available on at least two visits were considered for analysis. In addition, the genotype of HLA-DRB1, PTPN22 or STAT4 was not available (lack of DNA sample or undetermined) in 5%, 9% and 8.7% of cases, respectively. Therefore, we used data from 1,718 visits (449 patients; 3.8 visits/patient) to analyze disease activity and data from 1,813 visits (475 patients; 3.8 visits/patient) to analyze disability. To identify factors that influenced disease activity (DAS28 as the dependent variable) and disability (HAQ as the dependent variable) during the follow-up, we fitted two population-averaged models by generalized linear models nested by patient and visit using the xtgee command of Stata 10.1 for Windows (StataCorp LP, College Station, Texas, USA). The population-averaged generalized estimating equations were first modeled by adding all variables with a p value <0.15 in the bivariate analysis. The final models were constructed using quasi-likelihood estimation based on the independence model information criterion [15] and Wald tests, removing all variables with p>0.15.

In addition, to confirm the results obtained for disease activity we generated an ordered logistic model using the ologit command of Stata. The dependent variable was the disease activity level, using the cut-off points for DAS28 proposed by Prevoo et al. (<2.6, remission; 2.6 to 3.2, low disease activity; 3.2 to 5.1, moderate disease activity; >5.1 high disease activity [13]). The analysis was modeled as described above for xtgee, with remission considered 0 and low, moderate and high disease activity represented as 1, 2 and 3. The ordered logistic analysis estimates cut-off points that aid the interpretation of the coefficients for each independent variable according to the levels of the dependent variable.

Results

Population

At baseline, patients with definitive RA exhibited increased disease activity and disability, and higher acute phase reactant values than patients with UA (Table 1). Similarly, the percentages of RA patients positive for ACPA, RF and the SE were significantly higher than among UA patients (Table 1). However, UA patients were younger and tended to exhibit shorter disease duration than RA patients (Table 1). No significant differences were observed between the two subgroups in the distribution of the polymorphisms in STAT4 and PTPN22 studied (Table 1).

The STAT4 rs7574865 Polymorphism is Associated with the Evolution of Disease Activity in Early Arthritis Patients

The DAS28 values were significantly higher in elderly and female patients, those fulfilling the RA criteria, those with positive ACPA, and in patients from one of the centers (Table 2). After adjusting for these confounding variables, patients who were homozygous for the minor allele rs7574865 in STAT4 displayed significantly higher DAS28 values (β coefficient 0.42 [95% confidence interval: 0.01–0.83], p = 0.044: Table 2 and Figure 1), whereas the presence of a SE copy had no significant effect on disease activity. By contrast, the presence of the minor allele of rs2476601 in PTPN22 was associated dose-dependently with lower DAS28 values when compared with CC homozygous patients, although this effect was not significant for the TT genotype due to the low number of cases in our population (β coefficient for CT genotype −0.27 [−0.56– −0.01], p = 0.042; TT genotype −0.68 [−1.64– −0.27], p = 0.162: Figure 1 and Table 2).

To confirm our findings and obtain additional information, we analyzed the effects of these variables when disease activity is categorized as remission, or as low, moderate or high disease activity. Patients with the TT genotype of the rs7574865 allele in STAT4 displayed greater disease activity than heterozygous or GG genotype patients (Figure 2A). By contrast, the presence of a T allele of rs2476601 in PTPN22 was associated with more visits at weaker disease activity or remission (Figure 2B). Using this approach, multivariate analysis revealed a stronger association between the ACPA and TT genotype of rs7574865 in STAT4 and increased levels of disease activity, without significant changes in the contribution of the other independent variables when compared with the previous model (gender, age, etc.: Table 3). In addition, the cut-off points obtained in the ordered logistic regression model (rows at the bottom of Table 3) suggest that variables with positive coefficients, such as ACPA, gender or TT genotype of rs7574865 in STAT4, increased disease activity to moderate or high levels, while those with negative coefficients, such as the minor allele of rs2476601 in PTPN22, decreased the disease activity to the level of remission.

The Effect of the STAT4 rs7574865 Polymorphism on Disability Reported by Early Arthritis Patients

As expected, disability was significantly correlated with DAS28 values. Moreover, female gender and older age were associated with significantly higher HAQ values. After adjustment for these confounding variables, we detected a significant dose-dependent effect of the rs7574865 T allele of STAT4 on the HAQ value (Table 2 and Figure 3). Neither the presence of the SE, or of the PTPN22 rs2476601 minor allele, contributed to disability in our population.

Discussion

To the best of our knowledge, this is the first study demonstrating the influence of the STAT4 polymorphism rs7574865 on disease course and disability in EA. STAT4 is a member of the STAT family and it is implicated in the signals mediated by interleukin (IL) 12, IL-23 and type 1 interferons [16]. Although the mechanism by which the minor allele of rs7574865 contributes to the rupture of self-tolerance is not well understood, it is thought to involve dysregulation of helper T cell differentiation towards T_H1 and/or T_H17 phenotypes [16]. Kariuki et al reported that lupus patients with the minor rs7574865 allele of STAT4 exhibit lower levels of serum interferon-α but a comparatively stronger biological responses to this cytokine [17]. Moreover, the minor allele of rs7574865 has also been associated with severe disease in patients with systemic lupus erythematosus [18]. Therefore, the more intense disease activity observed in patients who were homozygous for the rs7574865 T allele may be associated with increased sensitivity to IL-12, IL-23 or interferons.

By contrast, the presence of the minor allele of rs2476601 in PTPN22 was associated with a trend towards diminished disease activity. The functional effect of this polymorphism is open to debate, although most studies suggest that the allele promotes a gain of function and hence, an increased threshold for T cell receptor (TCR) signaling that may lead to defective deletion of autoreactive clones in the thymus [3], [19]. Our results could support this hypothesis, suggesting that once self-tolerance is broken in patients with the rs2476601 T allele, autoreactive lymphocytes are activated less efficiently through TCR.

Surprisingly, the presence of the SE was not associated with increased disease activity or disability, probably due to the inclusion of ACPA in the statistical models. However, excluding this variable revealed no significant association between the presence of the SE and a worse outcome (data not shown). Although unexpected, there may be several explanations for this finding: a) a proportion of SE carriers were non-smokers and thus, they did not develop ACPA; b) some patients may be heterozygous carriers of SE and carry protective DRB1 alleles; or c) SE may be less important in RA patients from southern Europe than northern Europe. Taken together, these findings suggest that ACPA is a better marker of disease severity than SE. As with PTPN22, specific genetic factors may confer a risk of developing the disease, and while some may only modify disease severity, others may be involved in both processes.

Few studies have evaluated the effects of genetic markers that confer a risk of developing RA on disease activity and disability. In one study, the influence of 5 susceptibility genes on disease severity was analyzed, including that of DRB1, STAT4 and PTPN22 [11]. Only a weak association was reported between TRAF1/C5 and disease severity, as assessed by swollen joint counts, with no association of STAT4 or the remaining genes. The conflicting findings of this and our study may be due to the differences in study design. As opposed to the cross-sectional cohort study carried out previously, we performed a longitudinal study with almost 4 observations per patient, providing a more detailed picture of disease activity and disability. Our study involved patients from 2 clinical centers, a confounding variable for which the analysis of disease activity was adjusted. This is particularly important given the poor reproducibility of tender and swollen joint counts [20]. Furthermore, the study of Morgan and coworkers only adjusted the analysis for disease duration, which was longer than in our cohort. However, probably the most significant factor contributing to the divergent findings is the genetic background of the populations analyzed, as TRAF1/C5 and SE had a stronger influence on RA in the British population, and STAT4 in the Spanish population.

In summary, our data suggest that homozygosity for the STAT4 rs7574865 minor allele is associated with a poorer disease course and greater disability than homozygosity for the common allele of this variant. We propose that the rs7574865 polymorphism may be used as a molecular tool to identify patients that may benefit from early aggressive treatment.

Acknowledgments

We thank Dr. Loreto Carmona for her critical reading of the study and Mark Sefton for his editorial assistance. We also thank Ms. Teresa Velasco for her help in the Early Arthritis Clinic and Ms. Belen Diaz-Sanchez for her technical support.

Author Contributions

Conceived and designed the experiments: IG-A AB JM MEM-C. Performed the experiments: AL BR MFG-E MAL-N DP-S. Analyzed the data: IG-A. Wrote the paper: IG-A AL. Collection of clinical data: AMO IG-A AB LN.

Citation: Lamana A, Balsa A, Rueda B, Ortiz AM, Nuño L, Miranda-Carus ME, et al. (2012) The TT Genotype of the STAT4 rs7574865 Polymorphism Is Associated with High Disease Activity and Disability in Patients with Early Arthritis. PLoS ONE 7(8): e43661. https://doi.org/10.1371/journal.pone.0043661

References

1. MacGregor AJ, Snieder H, Rigby AS, Koskenvuo M, Kaprio J, et al. (2000) Characterizing the quantitative genetic contribution to rheumatoid arthritis using data from twins. Arthritis Rheum 43: 30–37.

2. Orozco G, Rueda B, Martin J (2006) Genetic basis of rheumatoid arthritis. Biomed Pharmacother 60: 656–662.

3. Gregersen PK, Olsson LM (2009) Recent advances in the genetics of autoimmune disease. Annu Rev Immunol 27: 363–391.

4. Orozco G, Alizadeh BZ, Delgado-Vega AM, Gonzalez-Gay MA, Balsa A, et al. (2008) Association of STAT4 with rheumatoid arthritis: a replication study in three European populations. Arthritis Rheum 58: 1974–1980.

5. Orozco G, Sanchez E, Gonzalez-Gay MA, Lopez-Nevot MA, Torres B, et al. (2005) Association of a functional single-nucleotide polymorphism of PTPN22, encoding lymphoid protein phosphatase, with rheumatoid arthritis and systemic lupus erythematosus. Arthritis Rheum 52: 219–224.

6. Goronzy JJ, Matteson EL, Fulbright JW, Warrington KJ, Chang-Miller A, et al. (2004) Prognostic markers of radiographic progression in early rheumatoid arthritis. Arthritis Rheum 50: 43–54.

7. Weyand CM, Xie C, Goronzy JJ (1992) Homozygosity for the HLA-DRB1 allele selects for extraarticular manifestations in rheumatoid arthritis. J Clin Invest 89: 2033–2039.

8. de Vries-Bouwstra JK, Goekoop-Ruiterman YP, Verpoort KN, Schreuder GM, Ewals JA, et al. (2008) Progression of joint damage in early rheumatoid arthritis: association with HLA-DRB1, rheumatoid factor, and anti-citrullinated protein antibodies in relation to different treatment strategies. Arthritis Rheum 58: 1293–1298.

9. Gonzalez-Gay MA, Gonzalez-Juanatey C, Lopez-Diaz MJ, Pineiro A, Garcia-Porrua C, et al. (2007) HLA-DRB1 and persistent chronic inflammation contribute to cardiovascular events and cardiovascular mortality in patients with rheumatoid arthritis. Arthritis Rheum 57: 125–132.

10. Marinou I, Maxwell JR, Wilson AG (2010) Genetic influences modulating the radiological severity of rheumatoid arthritis. Ann Rheum Dis 69: 476–482.

11. Morgan AW, Robinson JI, Conaghan PG, Martin SG, Hensor EM, et al. (2010) Evaluation of the rheumatoid arthritis susceptibility loci HLA-DRB1, PTPN22, OLIG3/TNFAIP3, STAT4 and TRAF1/C5 in an inception cohort. Arthritis Res Ther 12: R57.

12. Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, et al. (1988) The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum 31: 315–324.

13. Prevoo ML, van 't Hof MA, Kuper HH, van Leeuwen MA, van de Putte LB, et al. (1995) Modified disease activity scores that include twenty-eight-joint counts. Development and validation in a prospective longitudinal study of patients with rheumatoid arthritis. Arthritis Rheum 38: 44–48.

14. Esteve-Vives J, Batlle-Gualda E, Reig A (1993) Spanish version of the Health Assessment Questionnaire: reliability, validity and transcultural equivalency. Grupo para la Adaptacion del HAQ a la Poblacion Espanola. J Rheumatol 20: 2116–2122.

15. Pan W (2001) Model selection in estimating equations. Biometrics 57: 529–534.

16. Korman BD, Kastner DL, Gregersen PK, Remmers EF (2008) STAT4: genetics, mechanisms, and implications for autoimmunity. Curr Allergy Asthma Rep 8: 398–403.

17. Kariuki SN, Kirou KA, MacDermott EJ, Barillas-Arias L, Crow MK, et al. (2009) Cutting edge: autoimmune disease risk variant of STAT4 confers increased sensitivity to IFN-alpha in lupus patients in vivo. J Immunol 182: 34–38.

18. Taylor KE, Remmers EF, Lee AT, Ortmann WA, Plenge RM, et al. (2008) Specificity of the STAT4 genetic association for severe disease manifestations of systemic lupus erythematosus. PLoS Genet 4: e1000084.

19. Burn GL, Svensson L, Sanchez-Blanco C, Saini M, Cope AP (2011) Why is PTPN22 a good candidate susceptibility gene for autoimmune disease? FEBS Lett 585: 3689–3698.

20. Scott DL, Houssien DA (1996) Joint assessment in rheumatoid arthritis. Br J Rheumatol 35 Suppl 2: 14–18.