Abstract
Objective. In rheumatoid arthritis (RA), many genetic markers, such as the shared-epitope (SE) alleles, are described in association with radiographic progression, but limited data are available on undifferentiated arthritis (UA). We investigated whether single-nucleotide polymorphisms (SNP) and haplotypes in immune response genes and HLA class II alleles are associated with radiographic progression in patients with early UA.
Methods. Progression of radiographic damage was determined in white Dutch patients with early UA after 2 years of followup. Severe progression was defined as an increase in Sharp/van der Heijde Score ≥ 5 points after 2 years of followup. The remainder was classified as mild. These SNP were genotyped by Taqman technology: tumor necrosis factor (TNF) −1031, −863, −857, −308, −238; lymphotoxin-α (LTA) +368, +252; interleukin 10 (IL10) −2849, −1082, −819; IL1A −889, IL1B −31, +3953; and IL1RN +2018. Carriage of SE alleles and HLA-DQA1*05-DQB1*02 haplotype was established. These markers were analyzed in relation to radiographic progression.
Results. Forty-eight out of 151 patients with early UA had severe radiographic progression. Severe radiographic progression was associated with an increased carrier frequency of SE alleles (OR 5.12, 95% CI 2.0–13.1, p < 0.001) and IL10 GGC haplotype (OR 2.8, 95% CI 1.4–5.8, p = 0.003). Mild radiographic progression was associated with the HLA-DQA1*05-DQB1*02 haplotype (OR 0.3, 95% CI, 0.1–0.8, p = 0.013) and with allele TNF −308A (OR 0.4, 95% CI, 0.2–0.9, p = 0.02).
Conclusion. The SE and the IL10 GGC haplotype are associated with severe progression of radiographic damage, in contrast to the DQA1*05-DQB1*02 haplotype and the TNF −308A allele, which are associated with mild radiographic progression in early UA.
The progress of rheumatoid arthritis (RA) is characterized by a large variation in severity of inflammation and joint destruction. Joint damage often begins in the early stages of the disease1. Among patients with early arthritis there is a large group that cannot be classified into any established disease category2,3. The term “undifferentiated arthritis” (UA) was introduced to distinguish these unclassifiable arthritides from RA and other inflammatory rheumatic diseases. The frequency of UA in early arthritis cohorts ranges from 20% to 45%3 and UA can result in severe radiographic damage4,5. Of the patients with UA studied in Jansen, et al, 42% showed radiographic progression with an increase in Sharp-van der Heijde score (SHS) of at least 10 points within 1 year4. In contrast to RA, in which associations with single-nucleotide polymorphisms (SNP) have been extensively studied, fewer data are available on the genetic characteristics of patients with UA6–9.
The major histocompatibility complex (MHC) class II alleles, referred to as the shared-epitope (SE) alleles, contribute to 30%–50% of the genetic predisposition to RA10,11: HLA-DRB1*0401, *0404, *0405, *0408, *0101, *1001, and *1402. Moreover, these SE alleles contribute to the development of antibodies against cyclic citrullinated peptide (ACPA)7 and are associated with a worse radiographic outcome in RA12,13. In contrast with SE, HLA-DRB1*03 is more frequently present in the ACPA-negative patients with RA compared to ACPA-positive patients (OR 1.6, 95% CI 1.2–2.1) and is associated with a less severe disease14.
Apart from MHC class II alleles, SNP in genes coding for cytokines such as tumor necrosis factor-α (TNF), lymphotoxin-α (LTA), interleukin 1α (IL1A), IL1B, and the natural IL-1 receptor antagonist IL1RN play a pivotal role in the regulation of inflammation and seem to be important in the onset and progression of RA as well as UA.
TNF-α is a proinflammatory cytokine with a key function in the pathogenesis of RA15. Excessive production of TNF-α can incite synovial inflammation and proliferation as well as degradation of articular cartilage and bone, which causes joint damage. Patients with RA receiving anti-TNF treatment have a large decrease in disease activity and have significantly less radiographic progression compared to a placebo group16. The carrier frequency of the TNF −308A allele is associated with an increase in TNF-α concentrations in peripheral blood mononuclear cells cultured with T cell activators17.
LTA (formerly TNF-ß), a homolog of TNF-α, has similar biological activity as TNF-α and shares one of its receptors18. Therefore LTA also might be a susceptibility gene for RA and for UA. The TNF and LTA genes are arranged in tandem and map within the MHC centromeric to HLA-B and telomeric to the MHC class III genes.
Other potential candidates in the onset of UA are polymorphisms at the IL1 family loci IL1A, IL1B, and IL1RN. The activity of IL-1, which initiates and perpetuates inflammatory and destructive processes in the rheumatoid joint19, is normally balanced by endogenous inhibitors, e.g., IL-1ra20,21. It is known that an imbalance between IL-1 and IL-1ra may be a predisposing factor for RA20, most likely reflected by homozygosity for the allele 2 of the IL1RN 86-bp variable number of tandem repeats, which is in strong linkage disequilibrium with IL1RN +2018C allele22. A relation between radiographic progression and SNP at the IL1 cluster in RA was also demonstrated23,24.
The antiinflammatory cytokine IL-10 is negatively correlated with the progression of joint destruction in RA25. The IL10 −2849 promoter polymorphism is associated with autoantibody production and subsequent joint damage in RA26. Also, IL10 −592CC and IL10 −1082GG have been associated with more severe radiographic damage only in patients who are ACPA-negative27, while another study showed opposite results28.
We investigated whether SNP in cytokine genes associated with arthritis as well as the presence of SE alleles and the DQA1*05-DQB1*02 haplotype were associated with the progression of radiographic damage in patients with UA.
MATERIALS AND METHODS
The study population consisted of patients with early UA who were included in the Early Arthritis Cohort (EAC) of the Jan van Breemen Institute, a large rheumatology clinic in Amsterdam, The Netherlands. Since 1995, newly referred patients are included and followed prospectively in the EAC if they are aged ≥ 18 years, have peripheral arthritis in ≥ 2 joints, and a symptom duration of ≤ 3 years. Patients were excluded if they had been treated with a disease-modifying antirheumatic drug (DMARD) or if they had spondyloarthropathy, reactive arthritis, crystal-induced arthropathy, systemic lupus erythematosus, Sjögren’s syndrome, or osteoarthritis. The choice of treatment after inclusion was made by the rheumatologist. All included patients were Dutch and white, defined as patients with at least 3 Dutch white grandparents. The local Ethics Committee approved the study protocol, and all patients gave written informed consent.
Measurements
Data were collected on demographics, symptom duration, and disease activity in 28 joints (DAS28)29. The laboratory assessments included erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), IgM-rheumatoid factor (IgM-RF), and ACPA. Radiographs of the hands and feet were also made. Data on demographics and symptom duration were collected at baseline and data on the other variables were collected annually (Table 1).
Dependent variable
Radiographic damage was assessed with the SHS30 by 2 experienced rheumatologists who were blinded to the other variables and had a good intraclass correlation coefficient of 0.95. Severe progression was defined as an increase of the SHS of at least 5 points between baseline and 2 years of followup31. The remainder was classified as mild progression. After 2 years of followup, patients who fulfilled the 1987 American College of Rheumatology (ACR) criteria for RA32 were classified as having RA, and the remainder were classified as UA.
Genotyping
Genomic DNA was extracted from EDTA blood samples using a standard method. Genotyping of promoter SNP was performed at positions −1031 (db SNP ID: rs1799964), −863 (rs1800630), −857 (rs1799724), −308 (rs1800629), and −238 (rs361525) in the TNF gene and of intronic SNP at positions +252 (rs909253) and +368 (rs746868) in the LTA gene. Further, we studied IL1 family genes at position −889 (rs1800587) within the promoter of the IL1A gene, at position −31 in the promoter region (rs1143627), and another at position +3953 in exon 5 (rs1143634) of the IL1B gene and intron 2 SNP IL1RN +2018 (rs419598). Additionally, we studied IL10 promoter SNP IL10 −2849 (rs6703630), IL10 −1082 (rs1800896), and IL10 −819 (rs1800871).
SNP were genotyped with Taqman 5’ allelic discrimination technology (Applied Biosystems Inc., Foster City, CA, USA)33. Sequences of primers, fluorescent reporter dye-labeled probes, and reaction conditions were recorded.
Carriage of 1 or 2 HLA-DRB1 SE copies (i.e., SE+) was inferred from HLA-DQA1-DQB1 haplotypes in very strong linkage disequilibrium with HLA-DRB1 alleles in whites34,35. Polymerase chain reaction-amplified exon 2 DQA1 and DQB1 amplicons were generated for genotyping in a combined single-stranded conformation polymorphism/heteroduplex assay by a semiautomated electrophoresis and gel staining method on the Phastsystem (Amersham Pharmacia Biotech, Uppsala, Sweden). The method has been validated using a panel of reference DNA against the Dynal Allset+ high-resolution typing kit (Dynal A.S., Oslo, Norway)33.
These haplotypes could be reliably characterized in homozygous and heterozygous states33: HLA-DQA1*03-DQB1*0301 and DQA1*-03-DQB1*0302 (both DR4-related SE), DQA1*0101-DQB1*0501 (DR1 and DR10-related SE), and DQA1*0501-DQB1*02 (strictly DRB1*0301-related).
The SE+ status as established by the DQA1 and DQB1 PCR-SSCP/HD assay was validated in 87 established patients with RA by high-resolution sequence-based typing of HLA-DRB1 exon 2. To contain the SE, the HLA-DRB1*0101, *0102, *0401, *0404, *0405, *0408, *0410, and *1001 alleles were taken. The technique was used to correctly classify SE carriage in 99% of the patients. Independent confirmation in a second cohort of DRB1-typed Dutch patients with RA showed that with DQA1 and DQB1 typing, only 2 out of 167 patients (1.2%) would have been incorrectly classified36.
Statistical analysis
Haplotypes for SNP in the TNF-LTA region and the IL10 gene were inferred by use of PHASE 2.0, which implements a Bayesian algorithm to estimate the haplotype frequencies (http://www.stat.washington.edu/stephens/software.html).
Differences in baseline characteristics between patients with severe and mild radiographic progression were analyzed with Student’s t test, Mann-Whitney U test, and chi-squared test, as appropriate. Logistic regression analyses were used to determine associations between radiographic severity and the carrier frequency of alleles and haplotypes. Results are presented as OR with 95% CI, and the p value. If a significant effect was found for a single SNP, logistic regression was used to correct for linkage disequilibrium between distinct alleles. Logistic regression was also used for the determination of confounding of baseline variables. P values < 0.05 were considered statistically significant. Statistical analyses were performed using SPSS 15.0 Software (SPSS, Chicago, IL, USA).
RESULTS
A total of 151 patients with UA were eligible for our study. A minority of the patients (32%) were classified as having severe radiographic progression and the remainder (68%) as mild radiographic progression. At baseline, the mean age and symptom duration were comparable for both groups (Table 1). The mean DAS28, median CRP, mean ESR, percentage IgM-RF positivity, and ACPA positivity were higher in the severe disease group. The severe group had a higher SHS at baseline. After 2 years of followup, the number of DMARD used was significantly lower in the mild group compared to the severe group, while 50% of the patients in the severe group and 18% in the mild group fulfilled the 1987 ACR criteria for RA at that time.
In a univariate analysis, the carrier frequency of TNF −308A allele was significantly lower in the severe group (OR 0.43, 95% CI 0.20–0.90, p = 0.02). Carrier frequencies of other TNF or LTA minor alleles did not differ between patients with severe and mild radiographic progression (Table 2). The SE frequency was significantly lower in carriers of the TNF −308A allele (p < 0.001). Focusing on the TNF-LTA haplotypes, only the carrier frequency of haplotype E, tagged by TNF −308A allele, was significantly lower in the severe group (Table 3). A strong linkage disequilibrium was found between the TNF −308A allele and HLA-DQA1*05-DQB1*02 haplotype: 72.1% of the TNF −308A carriers were carriers of the HLA-DQA1*05-DQB1*02 haplotype. After stratification for radiographic progression, similar percentages were found (mild 53.8% vs progressive 77.1%).
Table 4 shows the carrier frequencies of (combinations of) HLA-DQA1*-DQB1* haplotypes. The numbers of HLA-DR4-related SE copies (OR 4.24, 95% CI 2.05–8.79, p = 0.001) and DR1-, DR4-, or DR10-related SE copies (OR 3.78, 95% CI 1.66–8.60, p = 0.001) were positively associated with an increased risk of severe progression of radiographic damage. Carrier frequency of HLA-DQA1*05-DQB1*02 was significantly lower in the severe disease group (OR 0.30, 95% CI 0.13–0.71, p = 0.013).
SNP in IL1 family genes and IL10 were similarly distributed among the groups (Table 5). The IL10 −1082A allele showed a trend toward a lower carrier frequency in the severe disease group compared to the mild group. Focusing on the IL10 haplotypes, only the carrier frequency of the GGC haplotype of the IL10 −2849G/A, IL10 −1082G/A, and IL10 −819C/T SNP was increased significantly in the severe group compared to the mild group (OR 2.84, 95% CI 1.40–5.75, p = 0.003; Table 6).
In order to detect additive or synergistic effects or interaction between HLA-DQA1*05-DQB1*02 and the IL10 GGC haplotype, logistic regression was performed. The DQA1*05-DQB1*02 haplotype and the IL10 GGC haplotype were both independently associated with progression of radiographic damage. The IL10 GGC haplotype was positively associated, while the DQA1*05-DQB1*02 haplotype was negatively associated.
To explore whether variables measured at baseline act as confounders of the SNP association, multivariate logistic regression was performed. The number of SE, SHS at baseline, CRP at baseline, and the number of DMARD used were not confounders. But ACPA positivity, IgM-RF positivity, DAS28 at baseline, and sex were marked as confounders. In a multivariate logistic model with IL10 GGC haplotype and HLA-DQA1*05-DQB1*02 haplotype, adjusted for these 4 confounders, HLA-DQA1*05-DQB1*02 haplotype was no longer a significant contribution in the model (IL10 GGC haplotype: OR 7.79, 95% CI 2.19–29.01, p = 0.002; and HLA-DQA1*05-DQB1*02: OR 0.33, 95% CI 0.09–1.19, p = 0.091).
DISCUSSION
In patients with early UA, we have shown the IL10 GGC haplotype to be related to severe radiographic progression, in contrast with the DQA1*05-DQB1*02 haplotype [strictly DRB1*0301 (DR3) related]34,35, which is associated with mild progression of radiographic damage after 2 years. The TNF −308A allele was also associated with mild progression, but was not independent from the HLA-DQA1*05-DQB1*02 haplotype. The higher carrier frequency of the TNF −308A allele in the mild disease group was rather unexpected because carriers of the TNF −308A were described as having an increased TNF-α production17.
Contradictory results are reported of association between the TNF −308 SNP and erosive outcome in RA37,38. One study in only IgM-RF-seropositive patients with early RA described an association between the TNF −308 AA+GA genotypes and a higher radiographic progression rate over 5 years, adjusted for the presence of SE39, while our result is in accord with a reported association between the TNF −308GG genotype and erosions in RA40. Other studies found no association between TNF −308 polymorphisms and radiographic progression, probably because of use of the Larsen score, which is a less sensitive method in an early stage of the disease38,41. Similarly, in a recent prospective study of patients with early-onset RA, the erosive damage at 1 year was not influenced by this polymorphism42.
In our study, baseline characteristics such as tender and swollen joints and the number of DMARD taken during the 2 years of followup were similar between carriers and non-carriers of the TNF −308A allele. It therefore seems unlikely that a difference in treatment may account for the reported association between TNF −308A carriage and mild radiographic progression. Therefore, it is hard to explain why patients with this TNF −308A allele, which is known for its association with a high production of proinflammatory TNF-α and a higher radiographic progression in RA, would have milder progression in our UA group. The fact that the TNF −308A carrier frequency was increased in patients with mild UA might be explained by the linkage disequilibrium with the DQA1*05-DQB1*02 haplotype, which is associated with a mild disease and might overrule the effect of the TNF −308A allele.
The carrier frequency of the HLA-DR4-related SE copies and DR1-, DR4-, or DR10-related SE copies was positively associated with severe progression, which is well known, and was independent from ACPA status at baseline. The lower carrier frequency of the DQA1*0501-DQB1*02 haplotype (in very strong linkage disequilibrium with the HLA-DRB1*03 allele) in the severe disease group compared to the mild group might be explained by the absence of HLA alleles carrying the SE, which are known to be associated with a more severe disease, or by the presence of protective HLA-DQ haplotypes9,43. In contrast, it was reported that the frequency of HLA-DRB1*03 in 88 Finnish patients with RA was associated with features of severe disease, such as extraarticular disease (p < 0.01) and prostheses in large joints (p < 0.05), but this was only a small group with a different genetic background44.
The TNF −308A allele and DQA1*0501-DQB1*02 were not independently associated with radiographic progression; this linkage disequilibrium was described before in RA39. Further, HLA-DRB1*0301 was significantly associated with ACPA-negative RA, while the absence of ACPA is associated with less radiographic progression in early RA14,45. The results in our patients with early arthritis are different from an early RA population, but the association between carriage of the TNF −308A allele, supported by the linkage disequilibrium with the DQA1*0501-DQB1*02 haplotype (DRB1*0301-related), are in accord with other reports in RA.
The carrier frequencies of other TNF alleles at −238, −857, −863, and −1031 and the LTA alleles at +368 and +252 were comparable for both radiographic groups. Some studies reported higher frequencies of the TNF −238GG genotype in patients with RA with severe disease46,47, others reported no effect of the TNF −238 polymorphism on radiographic progression48. No associations were found for the SNP TNF −238, −857, −863, and RA, although in a Taiwanese population, ethnic differences may have accounted for these discrepancies49.
The SNP in other cytokine genes (IL1A −889, IL1B −31, IL1B +3953, IL1RN +2018, IL10 −2849, IL10 −1082, and IL10 −819) were not associated with radiographic progression in UA, which is in agreement with earlier results on IL1B −31, IL1B +3953, IL1RN +2018, and IL10 −1082 observed in patients with early RA50. Of note, Cantagrel, et al50 reported that patients who were SE-positive and also carried the minor allele IL1B +3953T showed an increased risk of erosive disease, while patients carrying this allele and who were SE-negative did not. We did not find that result in our study. Although these observations are very interesting, they have to be considered preliminary and will need to be confirmed.
Our results showed a significant association of the promoter IL10 GGC haplotype with severe radiographic damage, but showed only a nonsignificant association with the IL10 −1082 SNP. This association was independent of the other associations found with progression of radiographic damage and therefore the influence of the IL10 haplotype at an early stage of the disease might be high. A recent report described that the IL10 −1082A/G polymorphism was associated with the production of rheumatoid factor28 but not with radiographic progression in early RA. Others described an association between the IL10 −2849 A/G promoter polymorphism and autoantibody production and joint damage in RA26. The fact that in our study on radiographic progression in early UA the association with the IL10 promoter GGC haplotype reached significance but these polymorphisms individually did not might point to the presence of another allele with influence on the GGC haplotype in the haplotype block encompassing the IL10 gene51.
Our study has some limitations. First, the population size was limited, which with stratification in the analyses may cause a lack of power to detect differences between the groups, resulting in a type II error. Larger studies might detect the effects of the described SNP and haplotypes better. Second, we are aware that the associations found in this study have to be viewed with caution, because bias due to multiple testing might exist.
Radiographic progression in early UA seems to be enhanced in patients carrying the GGC haplotype of the IL10 −2849G/A, IL10 −1082G/A, or IL10 −819C/T SNP, as well as the SE, and less severe in patients carrying the DQA1*05-DQB1*02 haplotype or the TNF −308A allele.
Acknowledgments
The authors thank E. de Wit-Taen and V. van der Lugt, clinical research assistants, for collecting patient data, R. Heijmans for genotyping the SNP, and W. Bos for valuable comments on the manuscript.
Footnotes
- Accepted for publication February 9, 2010.