TRAF1/C5polymorphism is not associated with increased mortality in rheumatoid arthritis: two large longitudinal studies
© van Nies et al.; licensee BioMed Central Ltd. 2010
Received: 28 September 2009
Accepted: 5 March 2010
Published: 5 March 2010
Recently an association between a genetic variation in TRAF1/C5 and mortality from sepsis or cancer was found in rheumatoid arthritis (RA). The most prevalent cause of death, cardiovascular disease, may have been missed in that study, since patients were enrolled at an advanced disease stage. Therefore, we used an inception cohort of RA patients to investigate the association between TRAF1/C5 and cardiovascular mortality, and replicate the findings on all-cause mortality. As TRAF1/C5 associated mortality may not be restricted to RA, we also studied a large cohort of non-RA patients.
615 RA patients from the Leiden Early Arthritis Clinic (EAC) (mean follow-up 7.6 years) were genotyped for rs10818488. In addition 5634 persons enrolled in the PROspective Study of Pravastatin in the Elderly at Risk (mean follow-up 3.2 years) were genotyped for rs2416808 (R2 >0.99 with rs10818488). The life/death status was determined and for the deceased persons the cause of death was ascertained. Cox proportional hazards and regression models were used to assess hazard ratios (HR) and 95% confidence intervals (CI).
Seventy-seven RA patients died. The main death causes in RA patients were cardiovascular diseases (37.7%), cancer (28.6%) and death due to infections (9.1%). No association was observed between the rs10818488 susceptible genotype AA and cardiovascular mortality (HR 1.08 95%CI 0.54 to 2.15) and all-cause mortality (HR 0.81 95%CI 0.27 to 2.43). Similar findings were observed for rs2416808 susceptible genotype GG in the non-RA cohort (HR 0.99; 95%CI 0.79 to 1.25 and HR 0.89; 95%CI 0.64 to 1.25, respectively).
The TRAF1/C5 region is not associated with an increased mortality risk.
Patients with rheumatoid arthritis (RA) have an increased mortality risk. A recent review, studying data from 84 unique cohorts, showed that mortality rates in RA patients were 1.5 to 1.6 fold higher than in the general population . The attributed causes of death in RA patients are identical to those in the general population , cardiovascular disease being the primary cause of death followed by cancer and infection. Age, sex and most clinical markers that are related to a more severe destructive disease course (among others number of inflamed joints, C reactive protein (CRP) and presence of erosions) are also associated with higher mortality risks .
Genetic risk factors for mortality in RA are scarcely investigated. Presence of the human leukocyte antigen (HLA)-DRB1 shared epitope alleles are reported to be associated with an increased mortality and, in particular, mortality related to cardiovascular disease (CVD) [2–4]. Although a lot of progress has been made in the field of genetics of RA-susceptibility, the HLA-shared epitope alleles still constitute the most powerful genetic risk factor to developing RA. Well-replicated non-HLA risk factors are PTPN22, TNFAIP3, and TRAF1/C5 [5–7]. TRAF1/C5 associated with several autoimmune diseases other than RA such as juvenile idiopathic arthritis (JIA) and systemic lupus erythomatodes [8, 9].
A recent study analyzed the TRAF1/C5 variant, rs3761847, in relation to mortality in RA and observed an increased mortality risk for the susceptible genotype GG compared with the non-susceptible AA genotype (hazard ratio (HR) = 3.96, 95% confidence interval (CI) = 1.24 to 12.6) . Such an observation is relevant because it indicates that genetic factors, such as TRAF1/C5, contribute to the increased mortality in RA. The causes of death were also investigated. Intriguingly, no increased death rate for CVD was found among GG homozygotes, whereas increased death rates were found for cancer and infections. As this study enrolled patients with a median disease duration of 10 (range 4 to 18) years, which were subsequently followed up for two to four years, CVD-related deaths occurring earlier in the disease course may have been missed. To further explore the association between the TRAF1/C5 locus and CVD-related mortality in RA, we studied a cohort of 615 early RA patients that were followed from disease onset to 14 years onwards. In addition, we investigated whether the association with all-cause, infectious and cancer-related mortality Panoulas and colleagues observed could be replicated . Finally, as an association between TRAF1/C5 and mortality may not be restricted to RA patients, we also analysed a large cohort of non-RA patients.
Materials and methods
Early arthritis clinic cohort
The Leiden early arthritis cohort (EAC) is an inception cohort consisting of patients with recent-onset arthritis referred to the Department of Rheumatology of the Leiden University Medical Center from 1993 onwards . Patients were included when arthritis was observed by a rheumatologist. For the present study, patients were selected who fulfilled the American College of Rheumatology 1987 revised criteria for RA within the first year of follow up and had DNA samples available (n = 615) . Written informed consent was obtained from all participants. The study was approved by the appropriate local institutional review board. At inclusion, a physical examination was performed and blood samples were taken to determine CRP, immunoglobulin (Ig) M rheumatoid factor (RF; by ELISA) and anti-cyclic citrullinated peptide (CCP) 2 antibodies (Immunoscan RA Mark 2; Euro-Diagnostica, Arnhem, the Netherlands). Anti-CCP2 positivity had a cut-off level of 25 arbitrary units, according to manufacturer's instructions. Patients also filled in a Health Assessment Questionnaire (HAQ)  and radiographs of hands and feet were taken and scored by one experienced reader (the intraclass-observer correlation coefficients was 0.91), using the Sharp-van der Heijde method .
Cohort of non-RA patients
As well as the RA patients, 5634 participants with available DNA from the PROspective Study of Pravastatin in the Elderly at Risk (PROSPER) were studied. In short, PROSPER is a randomized double-blind placebo-controlled trial that assessed whether pravastatin treatment in elderly men and women diminished the risk of major vascular events . Participants were screened and enrolled in Scotland (Glasgow), Ireland (Cork), and the Netherlands (Leiden) between December 1997 and May 1999. Written informed consent was obtained from all participants. The study was approved by the appropriate local institutional review boards of all centers. The primary endpoint in the study was a combination of death from coronary heart disease (CHD), non-fatal myocardial infarction (MI), and fatal or non-fatal stroke. After three years of follow up, CVD and transient ischemic attacks (TIAs) were less prevalent in the group treated with pravastatin .
rs10818488 was genotyped in the RA cohort and rs2416808, which is in complete linkage disequilibrium with rs10818488 (r2>0.99), was genotyped in the non-RA cohort. Genotypings were performed using the MassArray matrix-assisted laser desorption ionisation time-of-flight mass spectrometry, according to the protocols recommended by the manufacturer (Sequenom, San Diego, CA, USA). Each 384-well plate contained at least 4 positive (CEPH DNA) and 4 negative controls, to check for assay performance and contaminations, respectively. SpectroCaller software (Sequenom, San Diego, CA, USA) supplied by the manufacturer was used to automatically identify the genotypes. Clusters were checked and all doubtful calls were manually evaluated. Ten percent of the genotypes were performed in duplicate and the error rate was below 1%. Allele frequencies and Hardy-Weinberg equilibrium consistency were determined with Haploview . Both SNPs were in agreement with Hardy-Weinberg equilibrium. Both rs10818488 and rs2416808 are in complete linkage-disequilibrium with rs3761847 genotyped by Panoulas and colleagues (r2>0.99 data from Hapmap, and Kurreeman and colleagues) [5, 6, 10].
Notification of death
In the EAC, patients were followed longitudinally from the moment of their inclusion until 1 April, 2008, or death. All RA patients were tracked nationally using the civic registries (Gemeentelijke Basis Administratie) to ascertain life or death status. Causes of death for RA patients were obtained from Statistics Netherlands  and coded according to the International Classification of Diseases 10th revision of the World Health Organization . In the non-RA cohort, patients were followed for 3.2 years and the causes of death in this period were obtained from post-mortem reports and/or certification of death. All endpoints were adjudicated by a study endpoint committee.
Analysis of data
Data are expressed as mean (± standard deviation (SD)) with a 95% CI for continuous variables and as proportions for categorical variables. Differences in baseline patient characteristics between the single nucleotide protocol (SNP) genotypes were compared using a one-way analysis of variance test or Kruskal Wallis test for continuous variables and the chi-squared test for nominal variables. Associations between genotype and mortality were tested with univariate cox regression analyses and log rank tests.
In a study in RA patients (unpublished) we observed that age, CRP level and Sharp-van der Heijde score were independently associated with mortality. In order to investigate whether TRAF1/C5 is associated with mortality after adjusting for gender, HAQ-score and other previously found risk factors, a multivariate cox regression analysis was performed.
Subjects from the cohort of non-RA patients who withdrew consent or died during the study were censored at the date of death or at the last date of follow up. Cox regression analyses in the cohort of non-RA patients were adjusted for gender, age, pravastatin or placebo use, and country.
Panoulas and colleagues observed a HR of 3.96. The RA cohort has a power of 99.6% to identify such an association with C5/TRAF1 based on the observed genotype frequencies and an alpha of 0.05. In the non-RA cohort this power was 100%. Assuming that the finding by Panoulas and colleagues was affected by the winners curse and the true HR would be lower, for example a HR of 1.5, then the power of the RA cohort and non-RA cohort to observe an association was 28% and 97.5%, respectively.
All statistical analyses were performed using Statistical Package for Social Sciences version 16.0 (SPSS, Chicago, IL, USA). In all tests, P values below 0.05 were considered to be significant.
Baseline characteristics of the RA patients per genotype of rs10818488
n = 180 (29%)
n = 324 (53%)
n = 111 (18%)
56.2 ± 15.2
56.3 ± 15.5
57.3 ± 16.5
Females, n (%)
Past or current smoker, n (%)
Swollen joint count
9.0 ± 6.6
9.4 ± 7.3
8.0 ± 6.0
25.5 ± 3.7
25.7 ± 3.8
25.7 ± 4.0
RF positive, n (%)
Anti-CCP positive, n (%)
28.1 ± 29.0
31.6 ± 36.1
25.8 ± 30.1
1.05 ± 0.70
1.10 ± 0.72
1.03 ± 0.77
Total Sharp- van der Heijde score, median (IQR 25-75)
Baseline characteristics of the non-RA patients per genotype of rs2416808
n = 1783 (32%)
n = 2725 (48%)
n = 1123 (20%)
Females, n (%)
Current smoker, n (%)
Systolic blood pressure (mmHg)
Diastolic blood pressure (mmHg)
Total cholesterol (mmol/L)
History of diabetes, n (%)
History of hypertension, n (%)
History of myocardial infarction, n (%)
History of stroke or TIA, n (%)
History of vascular disease, n (%)
History of angina, n (%)
History of claudication, n (%)
Mortality in RA cohort
Hazard ratios for genotypes rs10818488 and rs2416808 in the RA-cohort and Non-RA cohort in univariate cox regression analysis
Hazard ratios RA cohort
AA vs GG
HR (95% CI)
Hazards ratios non-RA cohort
GG vs AA
HR (95% CI)
HR (95% CI)
Multivariate cox regression analyses were performed to assess whether the TRAF1/C5 susceptible genotype was associated with mortality after adjustments for other known risk factors for mortality. Also after adjustments, no significant association between rs10818488 genotype and all-cause (HR = 0.91, 95% CI = 0.61 to 1.37), cardiovascular (HR = 0.97, 95% CI = 0.54 to 1.76), cancer (HR = 0.66, 95% CI = 0.28 to 1.56) or infectious-related mortality (HR = 0.54, 95% CI = 0.09 to 3.13) was observed.
Mortality in non-RA cohort
In the non-RA cohort, 586 participants died, 32% were AA homozygote, 48% AG heterozygote and 20% GG homozygote for rs2416808. Also, no significant association was found between the TRAF1/C5 variants and all-cause mortality (HR = 0.99, 95% CI = 0.89 to 1.11), CVD-related (HR = 0.95, 95% CI = 0.80 to 1.12) and cancer-related mortality (HR = 1.00, 95% CI = 0.82 to 1.21; Table 3). Analyzing only the placebo-treated group did not change the results (data not shown). Similarly, also in the non-RA cohort after adjustments, no significant association was found between rs2416808 and all-cause (HR = 0.99, 95% CI = 0.88 to 1.11), cardiovascular (HR = 0.94, 95% CI = 0.80 to 1.11) or cancer-related mortalities (HR = 1.00, 95% CI = 0.82 to 1.22) in a multivariate cox regression analysis.
The present study was performed to investigate the relation between a well-replicated genetic RA susceptibility factor, TRAF1/C5, and the mortality risk in a cohort of RA patients. No evidence for an association of this risk factor with mortality was observed.
Panoulas and colleagues recently observed an increased mortality risk for carriers of the TRAF1/C5 susceptibility risk genotype in RA patients . Both polymorphisms analyzed in the present study, rs10818488 and rs2416808, are in complete linkage-disequilibrium with the rs3761847 SNP genotyped by Panoulas and colleagues (r2>0.99 data from Hapmap, and Kurreeman and colleagues) [5, 6]. However, we could not replicate their finding, despite having a larger cohort (615 vs 400 RA patients), a longer follow-up duration (mean 7.6 vs 2.6 years) and a higher number of events (77 vs 23 deaths). In addition, Panoulas and colleagues reported an increased risk of death due to cancer and sepsis, but not CVD. It is possible that the cohort design from Panoulas and colleagues was less suitable to study association with CVD-related mortality, because patients were enrolled at an advanced disease stage. In this setting, a possible relation between TRAF1/C5 and cardiovascular mortality may be missed, because patients dying of CVD earlier in the disease course are not part of their cohort. This hypothesis was corroborated by our results in the RA cohort, which showed that most of the CVD-related deaths were concentrated in the first 10 years of disease. However, also in our EAC, no association was found between TRAF1/C5 locus and cardiovascular mortality in RA.
In order to further unravel the eventual association between TRAF1/C5 locus and mortality, we hypothesised that this risk would not be restricted to RA. Therefore, we also genotyped TRAF1/C5 in a large cohort of elderly people that were prospectively followed in PROSPER. This large cohort was well-powered to detect also small HR (e.g. 1.5). Nonetheless, TRAF1/C5 did not confer an increased mortality risk.
This strengthened our findings of an absent association between TRAF1/C5 and risk of death.
In the end, one might argue that any cohort that enrolls older participants is inappropriate to study mortality, because the risk genotype may cause an early death. In this sense, prospective studies that follow the participants over a period of decades may be more suitable to answer this question. Nevertheless, the rs10818488 genotype frequencies in our JIA cohort (GG = 30%, AG = 53%, AA = 17%), which has a mean age of 6.4 years at inclusion, did not differ from the frequencies in the adult EAC cohort . Also the minor allele frequency in the non-RA cohort is similar to those in the general population . In addition, both cohorts were in agreement with Hardy-Weinberg equilibrium. This indicates that there was no selection of the protective genotype in the RA or non-RA cohort. This makes an early death in the risk genotype carriers unlikely and further supports the lack of association of these TRAF1/C5 variants with mortality.
In the RA-cohort, mortality due to infections were relatively infrequent, the number of deaths attributed to infections is insufficient to make definite conclusions on the association between TRAF1/C5 and this specific cause of mortality.
The important causes of death observed in RA patients in the present study are similar to the main causes of death in the general Dutch population . The frequency of CVD-related mortality itself is reported to be higher in RA patients than in healthy individuals . Results from the Nurses Health Study revealed that women with RA had a relative risk of 1.8 for fatal MI . Unfortunately, we were unable to test whether the frequency of cardiovascular death was also increased in our EAC RA cohort, because we do not have mortality information on an age- and gender-matched Dutch control population.
In conclusion, TRAF1/C5 polymorphisms predisposing to RA susceptibility are not associated with all-cause mortality or cardiovascular- or cancer-related mortality in RA and in an elderly cohort of persons without RA.
cyclic citrullinated peptide
coronary heart disease
C reactive protein
early arthritis clinic
enzyme-linked immunosorbent assay
health assessment questionnaire
human leukocyte antigen
juvenile idiopathic arthritis
PROspective Study of Pravastatin in the Elderly at Risk
single nucleotide protocol
transient ischaemic attack.
The authors thank the Dutch Central Bureau for Statistics (CBS) for providing the data on the survival of the general Dutch population and the death causes of RA patients. This work is supported by The Netherlands Organisation for Health Research and Development, Centre for Medical Systems Biology and the Dutch Arthritis Association and the European Union-funded FP7-integrated project Masterswitch no. 223404 and AutoCure.
- Sokka T, Abelson B, Pincus T: Mortality in rheumatoid arthritis: 2008 update. Clin Exp Rheumatol. 2008, 26: S35-S61.PubMedGoogle Scholar
- Mattey DL, Thomson W, Ollier WE, Batley M, Davies PG, Gough AK, Devlin J, Prouse P, James DW, Williams PL, Dixey J, Winfield J, Cox NL, Koduri G, Young A: Association of DRB1 shared epitope genotypes with early mortality in rheumatoid arthritis: results of eighteen years of followup from the early rheumatoid arthritis study. Arthritis Rheum. 2007, 56: 1408-1416. 10.1002/art.22527.View ArticlePubMedGoogle Scholar
- Farragher TM, Goodson NJ, Naseem H, Silman AJ, Thomson W, Symmons D, Barton A: Association of the HLA-DRB1 gene with premature death, particularly from cardiovascular disease, in patients with rheumatoid arthritis and inflammatory polyarthritis. Arthritis Rheum. 2008, 58: 359-369. 10.1002/art.23149.PubMed CentralView ArticlePubMedGoogle Scholar
- Gonzalez-Gay MA, Gonzalez-Juanatey C, Lopez-Diaz MJ, Pineiro A, Garcia-Porrua C, Miranda-Filloy JA, Ollier WE, Martin J, Llorca J: HLA-DRB1 and persistent chronic inflammation contribute to cardiovascular events and cardiovascular mortality in patients with rheumatoid arthritis. Arthritis Rheum. 2007, 57: 125-132. 10.1002/art.22482.View ArticlePubMedGoogle Scholar
- Kurreeman FA, Padyukov L, Marques RB, Schrodi SJ, Seddighzadeh M, Stoeken-Rijsbergen G, Helm-van Mil van der AH, Allaart CF, Verduyn W, Houwing-Duistermaat J, Alfredsson L, Begovich AB, Klareskog L, Huizinga TW, Toes RE: A candidate gene approach identifies the TRAF1/C5 region as a risk factor for rheumatoid arthritis. PLoS Med. 2007, 4: e278-10.1371/journal.pmed.0040278.PubMed CentralView ArticlePubMedGoogle Scholar
- Plenge RM, Seielstad M, Padyukov L, Lee AT, Remmers EF, Ding B, Liew A, Khalili H, Chandrasekaran A, Davies LR, Li W, Tan AK, Bonnard C, Ong RT, Thalamuthu A, Pettersson S, Liu C, Tian C, Chen WV, Carulli JP, Beckman EM, Altshuler D, Alfredsson L, Criswell LA, Amos CI, Seldin MF, Kastner DL, Klareskog L, Gregersen PK: TRAF1-C5 as a risk locus for rheumatoid arthritis--a genomewide study. N Engl J Med. 2007, 357: 1199-1209. 10.1056/NEJMoa073491.PubMed CentralView ArticlePubMedGoogle Scholar
- Helm-van Mil van der AH, Huizinga TW: Advances in the genetics of rheumatoid arthritis point to subclassification into distinct disease subsets. Arthritis Res Ther. 2008, 10: 205-10.1186/ar2384.PubMed CentralView ArticlePubMedGoogle Scholar
- Albers HM, Kurreeman FA, Houwing-Duistermaat JJ, Brinkman DM, Kamphuis SS, Girschick HJ, Wouters C, Van Rossum MA, Verduijn W, Toes RE, Huizinga TW, Schilham MW, ten Cate R: The TRAF1/C5 region is a risk factor for polyarthritis in juvenile idiopathic arthritis. Ann Rheum Dis. 2008, 67: 1578-1580. 10.1136/ard.2008.089060.View ArticlePubMedGoogle Scholar
- Kurreeman FA, Goulielmos GN, Alizadeh BZ, Rueda B, Houwing-Duistermaat J, Sanchez E, Bevova M, Radstake TR, Vonk MC, Galanakis E, Ortego N, Verduyn W, Zervou MI, Consortium S, Roep BO, Dema B, Espino L, Urcelay E, Boumpas DT, Berg van den LH, Wijmenga C, Koeleman BP, Huizinga TW, Toes RE, Martin J: The TRAF1-C5 region on chromosome 9q33 is associated with multiple autoimmune diseases. Ann Rheum Dis. 2009.Google Scholar
- Panoulas VF, Smith JP, Nightingale P, Kitas GD: Association of the TRAF1/C5 locus with increased mortality, particularly from malignancy or sepsis, in patients with rheumatoid arthritis. Arthritis Rheum. 2009, 60: 39-46. 10.1002/art.24176.View ArticlePubMedGoogle Scholar
- van Aken J, van Bilsen JH, Allaart CF, Huizinga TW, Breedveld FC: The Leiden Early Arthritis Clinic. Clin Exp Rheumatol. 2003, 21: S100-S105.PubMedGoogle Scholar
- Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, Medsger TA, Mitchell DM, Neustadt DH, Pinals RS, Schaller JG, Sharp JT, Wilder RL, Hunder GG: The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988, 31: 315-324. 10.1002/art.1780310302.View ArticlePubMedGoogle Scholar
- Heijde van der DM: How to read radiographs according to the Sharp/van der Heijde method. J Rheumatol. 1999, 26: 743-745.PubMedGoogle Scholar
- Heijde van der DM, van Riel PL, Putte van De LB: Sensitivity of a Dutch Health Assessment Questionnaire in a trial comparing hydroxychloroquine vs. sulphasalazine. Scand J Rheumatol. 1990, 19: 407-412. 10.3109/03009749009097629.View ArticlePubMedGoogle Scholar
- Shepherd J, Blauw GJ, Murphy MB, Cobbe SM, Bollen EL, Buckley BM, Ford I, Jukema JW, Hyland M, Gaw A, Lagaay AM, Perry IJ, Macfarlane PW, Meinders AE, Sweeney BJ, Packard CJ, Westendorp RG, Twomey C, Stott DJ: The design of a prospective study of Pravastatin in the Elderly at Risk (PROSPER). PROSPER Study Group. PROspective Study of Pravastatin in the Elderly at Risk. Am J Cardiol. 1999, 84: 1192-1197. 10.1016/S0002-9149(99)00533-0.View ArticlePubMedGoogle Scholar
- Shepherd J, Blauw GJ, Murphy MB, Bollen EL, Buckley BM, Cobbe SM, Ford I, Gaw A, Hyland M, Jukema JW, Kamper AM, Macfarlane PW, Meinders AE, Norrie J, Packard CJ, Perry IJ, Stott DJ, Sweeney BJ, Twomey C, Westendorp RG, PROSPER study group, PROspective Study of Pravastatin in the Elderly at Risk: Pravastatin in elderly individuals at risk of vascular disease (PROSPER): a randomised controlled trial. Lancet. 2002, 360: 1623-1630. 10.1016/S0140-6736(02)11600-X.View ArticlePubMedGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.View ArticlePubMedGoogle Scholar
- Statistics Netherlands. [http://statline.cbs.nl/StatWeb/]
- World Health Organization. [http://apps.who.int/classifications/apps/icd/icd10online/]
- NCBI. [http://www.ncbi.nlm.nih.gov/]
- Solomon DH, Karlson EW, Rimm EB, Cannuscio CC, Mandl LA, Manson JE, Stampfer MJ, Curhan GC: Cardiovascular morbidity and mortality in women diagnosed with rheumatoid arthritis. Circulation. 2003, 107: 1303-1307. 10.1161/01.CIR.0000054612.26458.B2.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.