Meta-analysis identified the TNFA -308G > A promoter polymorphism as a risk factor for disease severity in patients with rheumatoid arthritis
Arthritis Research & Therapy volume 14, Article number: R264 (2012)
The goal of this study is to investigate whether the -308G > A promoter polymorphism in the tumor necrosis factor alpha (TNFA) gene is associated with disease severity and radiologic joint damage in a large cohort of patients with rheumatoid arthritis (RA).
A long-term observational early RA inception cohort (n = 208) with detailed information about disease activity and radiologic damage after 3, 6 and 9 years of disease was genotyped for the TNFA -308G > A promoter polymorphism (rs1800629). A longitudinal regression analysis was performed to assess the effect of genotype on RA disease severity and joint damage. Subsequently, a meta-analysis, including all publically available data, was performed to further test the association between joint erosions and the TNFA polymorphism. To learn more about the mechanism behind the effect of the polymorphism, RNA isolated from peripheral blood from RA patients (n = 66) was used for TNFA gene expression analysis by quantitative PCR.
Longitudinal regression analysis with correction for gender and disease activity showed a significant difference in total joint damage between GG and GA+AA genotype groups (P = 0.002), which was stable over time. The meta-analysis, which included 2,053 patients, confirmed an association of the genetic variant with the development of erosions (odds ratio 0.78, 95% CI 0.62, 0.98). No significant differences in TNFA gene expression were observed for the different genotypes, confirming earlier findings in healthy individuals.
Our data confirm that the TNFA -308G > A promoter polymorphism is associated with joint damage in patients with RA. This is not mediated by differences in TNFA gene expression between genotypes.
TNF plays a key role in the pathogenesis of rheumatoid arthritis (RA) and other auto-immune diseases. Several genetic variants in the TNFA gene, which codes for TNFα, have been investigated in relation to disease severity, but none of them appears to be highly specific and/or sensitive. The most thoroughly investigated TNFA polymorphism is the -308G > A (rs1800629) located in the promoter of the gene, though studies focusing on this polymorphism have yielded conflicting results, [1–14].
Three studies [3, 6, 7] showed an association of the G allele with more pronounced radiologic joint damage. Other studies reported a worse radiologic outcome in patients carrying the A allele [1, 9, 13], or no association between the genetic variant and radiologic progression [2, 8–12, 14]. Given the contradictory results in these relatively small, individual studies, there is a need for a more comprehensive approach to evaluate the presence or absence of an association between the TNFA -308G > A polymorphism and the severity of RA and radiological joint damage.
As the -308G > A variant is located within the TNFA promoter, several studies have also investigated the polymorphism in relation to TNFA transcription, which could influence TNFα production and in turn, disease activity and severity [2, 15–23]. An important study by Knight et al., using the haploChIP method for high-throughput screening of common DNA polymorphisms that (might) effect gene regulation in vivo, reported that the -308G > A polymorphism was not associated with TNFA expression in human B-cells . A recent meta-analysis in healthy individuals (n = 1,825) shows that TNFα mRNA and protein levels are not influenced by the polymorphism .
Here, we first investigated the possible associations of the -308G > A promoter polymorphism with disease activity, and radiologic joint damage after 3, 6 and 9 years follow-up in a large (n = 208), well-characterized cohort of patients with early RA . Next, we performed a meta-analysis to draw more definitive conclusions on the role of the TNFA -308G > A polymorphism in the severity of RA and radiological joint damage. This meta-analysis included a total of 2,053 patients and combined our data with previously published studies. To evaluate possible mechanisms underlying an association, we investigated gene expression in the different genotype groups in a cohort of 66 individuals with active RA.
Material and methods
The Commissie Mensgebonden Onderzoek (CMO) Regio Arnhem Nijmegen of the Radboud University Nijmegen Medical Centre approved the study (CMO number 2004/014). All patients had provided written informed consent prior to participation in the study. All clinical investigations were conducted according to the principles expressed in the Declaration of Helsinki.
All patients in this study fulfilled the 1987 American College of Rheumatology (ACR) criteria for RA and attended the Department of Rheumatology of the Radboud University Nijmegen Medical Centre or the St. Maartenshospital in Nijmegen, The Netherlands.
The impact of the TNFA -308G > A polymorphism on disease severity was analyzed using DNA isolated from blood of 208 patients from a long-term observational inception cohort with early RA . Patients in this cohort have a disease duration < 1 year and no prior use of DMARDs at enrollment. Disease activity was assessed using the 28-joint disease activity score (DAS28) at each 3-monthly visit. Radiologic joint damage at baseline and after 1, 2, 3, 6 and 9 years of disease duration was measured using radiographs of the hands and feet, read in chronological order by one of four raters according to the Ratingen score. The latter is a validated modification of the Larsen score and evaluates joint surface destruction, graded from 0 to 5, in 38 separate hand and foot joints (range 0 - 190) . The intra-class correlation coefficient (ICC) for inter-rater reliability was 0.85; this was tested previously with the four raters in 10 patients over 9 years.
TNFA expression analysis was assessed in RNA from blood of 66 additional patients with active RA as defined by DAS28 > 3.2. All patients attended the same hospital, and blood for RNA isolation was always collected between 9 and 10 o'clock in the morning to avoid potential circadian fluctuations.
All DNA and RNA analyses were performed in a Coördinatie Commissie ter bevordering van de Kwaliteitsbeheersing van het Laboratoriumonderzoek (CCKL)-accredited laboratory at the department of Human Genetics at the Radboud University Nijmegen Medical Centre in Nijmegen.
Genotyping of the -308G > A TNFA polymorphism
Genomic DNA was extracted from peripheral venous blood (10 ml/sample) according to standard protocols . Both patient cohorts, thus, a total of 274 patients, were genotyped for the -308G > A (rs1800629) promoter polymorphism using a TaqMan genotyping assay (ID: C7514879_10) on the 7500 Fast Real-Time PCR system (Applied Biosystems, Foster City, CA, USA). Results were analyzed using Allelic Discrimination software version 1.4 (Applied Biosystems).
RNA isolation, synthesis of cDNA and quantitative PCR
RNA was isolated from whole blood within 30 minutes after collection, using the RNeasy midi kit (Qiagen Benelux BV Venlo, The Netherlands). The quality and quantity of the purified RNA was checked on a NanoDrop spectrophotometer (Nanodrop technologies, Montchanin, DE, USA). Degradation of RNA was controlled for by agarose gel electrophoresis. Per sample, 1 μg of total RNA was reverse-transcribed using 200 ng random hexanucleotides (Invitrogen, Breda, The Netherlands) and 200 units of moloney murine leukemia virus (M-MLV) Reverse Transcriptase (Invitrogen). For quantitative expression analysis a predesigned and validated gene-specific probe-based TaqMan gene expression assay (ID: Hs00174128_m1) was used according to the manufacturer's protocol (Applied Biosystems). Samples were run on the 7500 Fast Real-Time PCR System (Applied Biosystems) using standard protocols. The general housekeeping gene B2M (Beta-2-microglobulin, TaqMan gene expression assay ID: Hs99999907_m1) was used as the endogenous control. Threshold cycle numbers (referred to as Ct) were obtained using the 7500 System SDS software version 1.4 (Applied Biosystems). All samples were measured twice, and duplicate samples with a SD larger than 0.5 were excluded from the analysis. The relative quantity (RQ) of the gene-specific mRNA was calculated from the average value of the ΔCt ((TNFA Ct) - (endogenous control gene Ct)). Differences in expression between two groups of samples were calculated by the 2ΔΔCt method .
Hardy-Weinberg equilibrium (HWE) was tested in both patient samples using the chi square test. Baseline differences between genotype groups were analyzed using the chi square test, Student's t-test or Wilcoxon test, as appropriate. Longitudinal analysis of radiologic joint damage in patients with RA was performed using generalized linear regression models with random coefficients (mixed models), allowing for analysis of repeated measures within a patient. The joint damage score was the dependent variable, genotype the independent variable, under addition of time, time squared and confounders. Regression assumptions were checked by evaluating the distribution of residuals and by a plot of observed and predicted joint damage scores. Mean differences in TNFA expression (ΔCt; (TNFA Ct) - (endogenous control gene Ct)) were analyzed using the independent Student's t-test.
A P-value < 0.05 was considered statistically significant in each situation. Analyses were performed using SPSS for Windows, version 14.0 (SPSS, Chicago, IL, USA) and SAS version 8.2 (SAS Institute Inc., Cary, NC, USA).
Pubmed (up to March 31, 2012) was searched to identify publications investigating the effect of the -308G > A polymorphism on joint damage. Using the search terms: 1) RA AND TNF AND (polymorphism or -308); 2) RA AND joint damage AND (polymorphism or TNF or -308 and TNF; 3) RA and erosions in combination with polymorphism and TNF or -308. The literature search yielded 292 papers, 37 of which reported on relevant primary analysis of the TNFA -308G > A polymorphism. In order to be included in this meta-analysis, studies had to 1) include RA patients according to the ACR criteria, and 2) report on disease severity (presence or absence of erosions). Review manager 5 was used to perform a meta-analysis including 11 studies (10 published ones and this study) [1–10].
Homogeneity of odds ratios (ORs) among cohorts was calculated using the Breslow-Day method, and pooled ORs were calculated under a fixed and random effects model (Mantel-Haenszel meta-analysis).
Relation of TNFA -308G > A with disease activity and radiologic joint damage
The frequency of the A allele in the inception cohort was 18.3%, genotype frequencies were 66.3% for the GG group, 30.8% for the GA group and 2.9% for the AA group. Due to the small AA genotype group, the GA and AA genotype groups were merged and analyzed as one group. Demographic and disease characteristics were similar in the GG and GA+AA genotype groups at baseline (Table 1). No significant differences in disease activity were observed between genotype groups at baseline (P = 0.38) and after 3 (P = 0.85), 6 (P = 0.1) and 9 (P = 0.12) years of follow-up.
In the longitudinal regression analysis (mixed model) with radiologic joint damage as the dependent variable, a significant difference between the GG and GA+AA groups in the progression of joint damage was observed (Table 2 and Figure 1). As shown in Table 2, the progression of joint damage in the patients was 1.6 points per year, with progression rate slightly declining (-0.13 points) over the years of study. At baseline, the joint damage score of the genotype GA+AA group was on average 5 points lower (P = 0.002) than the score of the GG genotype group, corrected for age and DAS28 at baseline. There was no significant interaction between genotype and time of study, that is, the difference in joint damage score between genotype groups persisted during follow-up (P = 0.61), although in Figure 1 it appears that the joint damage scores of both groups converge at year-9.
Since TNF might uncouple the relation between disease activity and joint damage , we assessed whether the use of TNF blocking agents was frequent in this patient population. The number of patients with past or present treatment with these agents was low, and was similar in both genotype groups (n = 20 (14%) and n = 12 (17%)) in the GG and GA+AA groups, respectively.
A total of 2,053 RA patients from 10 published studies, as well as the current one, were included in a meta-analysis to further confirm the observed association of TNF -308G > A with joint damage, in this case defined as the presence or absence of erosions. Table 3 presents an overview of characteristics and the definition of erosions of the studies included in the meta-analysis. The analysis indeed confirmed that the -308G > A polymorphism is statistically significantly associated with the presence of erosions (P = 0.04, OR 0.78, 95% CI 0.62, 0.98) (Figure 2), with patients carrying the GA/AA genotypes having less erosions. As we observed moderate heterogeneity, we also applied a random effects meta-analysis, producing a similar OR (0.86, 95% CI 0.58, 1.29).
TNFA expression and the role of the -308G > A promoter polymorphism
To assess whether the difference in radiologic joint damage between G-homozygotes and A-allele carriers could be explained at the level of expression of this cytokine, we measured the RNA expression levels of TNFA dependent on -308G > A genotype in a sample of 66 patients with active RA. Among the sixty-six patients, forty-nine carried the GG genotype, fifteen were GA-heterozygotes, and two were A-homozygotes. Again, the GA and AA genotypes were merged and analyzed as one group. No difference in TNFA expression was observed between the GG group when compared to the GA+AA group (P = 0.420) (Figure 3).
Based on earlier literature, we tested the hypothesis that in RA, disease severity defined by the average DAS28 and radiologic joint damage, may be influenced by the TNFA -308G > A genotype. Our own study and the meta-analysis indeed confirm that patients carrying the common GG genotype have worse radiologic outcomes. This could not be explained by differences in either disease activity, genotype-dependent TNFA expression, or the use of TNF blocking agents.
Our results show that the GG genotype is associated with a worse radiologic outcome, in spite of the fact that in our early inception cohort the radiologic progression after 9 years was similar in both genotypes. It would be of great interest to see if these results could be replicated in other early RA cohorts though the number of published studies hitherto is too small for a meta-analysis.
Many studies investigating the relationship between the -308G > A promoter polymorphism and severity of RA and radiological joint damage report contradictive results [1–14]. Therefore, a meta-analysis offers a good alternative for evaluating the presence or absence of an association in situations where individual studies are inconclusive. Our meta-analysis indeed confirmed that the -308G > A polymorphism is associated with the presence of erosions, with patients carrying the GA/AA genotypes having fewer erosions. Although the result of the meta-analysis is consistent with the results of our primary study, it is important to mention that there was significant heterogeneity in the studies included in the meta-analysis (I2 = 53%). This could be caused by differences in study design, mean duration of RA, sample size and ethnicity. An additional reason is the assessment of the presence of erosions, which was assessed in radiographs of different joints (see Table 3), though most studies included erosions in the hands. Also the definition of erosive disease differed between the studies, but the outcome used in the current study was the presence of any erosion in the joints examined. As we observed heterogeneity we also performed a random-effects meta-analysis. This yielded a similar OR but, as expected, a larger CI. Even so, both analyses point to a modest effect of the polymorphism on the development of joint damage.
In this study we also assessed whether the difference in radiologic joint damage between G- homozygotes and A-allele carriers could be explained by TNFA allele-specific differences in mRNA expression. Many studies have focused on the potential allele-specific expression of the TNFA gene and TNFα production. Helmig and co-workers reported that the -308G allele is associated with TNFA expression in a cohort of 178 healthy subjects . In contrast, a recent meta-analysis in which 1,825 healthy subjects were investigated, reported no association between the -308G > A polymorphism and TNF protein and mRNA levels . However, this might still be the case for persons suffering from a disease. For instance, the study of Muoz-Valle, et al. (50 patients with osteoarthritis), and the study of Oregon-Romero, et al. (50 patients with RA) both report an association between the -308G allele and higher TNFA expression levels compared to the -308A allele [21, 22]. Nonetheless, our study does not support this hypothesis and shows that the effect of the -308G > A variant on joint damage is not explained by allele-dependent expression in RA. Linkage disequilibrium (LD) could be an overall factor to explain this incongruent result if we assume that the -308G > A variant is not functional. In fact, LD is strong in this area of the genome, and it is difficult to study the role of a single nucleotide polymorphism (SNP) in this region. Nevertheless, a published haplotype analysis including three SNPs (-857C > T, -308G > A, and -238G > A) in the promoter of the TNFA gene also showed no association between the SNPs and circulating TNF levels or TNFA mRNA expression . It is known that the LD in the TNFA region even extends much further and is in strong LD with the shared epitope (SE) . This could explain some of the inconsistency in the findings, although the results of a study of Khanna and co-workers indicates an effect of the -308G > A polymorphism on joint damage progression independently of the SE . Our patients are not genotyped for the SE and we could not stratify for SE.
Our data support an association between the TNFA -308GG genotype and joint damage, but not with disease activity or TNFA expression. Our meta-analysis shows that the effects of this polymorphism are modest, which makes it unsuitable as a single predictive marker for disease severity.
American College of Rheumatology
Coördinatie Commissie ter bevordering van de Kwaliteitsbeheersing van het Laboratoriumonderzoek
28-joint disease activity score
disease-modifying antirheumatic drug
intra-class correlation coefficient
moloney murine leukemia virus
polymerase chain reaction
single nucleotide polymorphism
tumor necrosis factor.
Mosaad YM, Abdelsalam A, El-Bassiony SR: Association of tumour necrosis factor-alpha -308 G/A promoter polymorphism with susceptibility and disease profile of rheumatoid arthritis. Int J Immunogenet. 2011, 38: 427-433. 10.1111/j.1744-313X.2011.01028.x.
Brinkman BM, Huizinga TW, Kurban SS, van der Velde , Schreuder GM, Hazes JM, Breedveld FC, Verweij CL: Tumour necrosis factor alpha gene polymorphisms in rheumatoid arthritis: association with susceptibility to, or severity of, disease?. Br J Rheumatol. 1997, 36: 516-521. 10.1093/rheumatology/36.5.516.
Hussein YM, Mohamed RH, Pasha HF, El-Shahawy EE, Alzahrani SS: Association of tumor necrosis factor alpha and its receptor polymorphisms with rheumatoid arthritis in female patients. Cell Immunol. 2011, 271: 192-196. 10.1016/j.cellimm.2011.06.023.
Reneses S, Gonzalez-Escribano MF, Fernandez-Suarez A, Pestana L, Davila B, Wichmann I, Garcia A: The value of HLA-DRB1 shared epitope, -308 tumor necrosis factor-alpha gene promoter polymorphism, rheumatoid factor, anti-citrullinated peptide antibodies, and early erosions for predicting radiological outcome in recent-onset rheumatoid arthritis. J Rheumatol. 2009, 36: 1143-1149. 10.3899/jrheum.081075.
Ates O, Hatemi G, Hamuryudan V, Topal-Sarikaya A: Tumor necrosis factor-alpha and interleukin-10 gene promoter polymorphisms in Turkish rheumatoid arthritis patients. Clin Rheumatol. 2008, 27: 1243-1248. 10.1007/s10067-008-0893-1.
Nemec P, Pavkova-Goldbergova M, Stouracova M, Vasku A, Soucek M, Gatterova J: Polymorphism in the tumor necrosis factor-alpha gene promoter is associated with severity of rheumatoid arthritis in the Czech population. Clin Rheumatol. 2008, 27: 59-65. 10.1007/s10067-008-0924-y.
Rezaieyazdi Z, Afshari JT, Sandooghi M, Mohajer F: Tumour necrosis factor a -308 promoter polymorphism in patients with rheumatoid arthritis. Rheumatol Int. 2007, 28: 189-191. 10.1007/s00296-007-0444-0.
Kazkaz L, Marotte H, Hamwi M, Angelique CM, Roy P, Mougin B, Miossec P: Rheumatoid arthritis and genetic markers in Syrian and French populations: different effect of the shared epitope. Ann Rheum Dis. 2007, 66: 195-201.
Khanna D, Wu H, Park G, Gersuk V, Gold RH, Nepom GT, Wong WK, Sharp JT, Reed EF, Paulus HE, Tsao BP: Association of tumor necrosis factor alpha polymorphism, but not the shared epitope, with increased radiographic progression in a seropositive rheumatoid arthritis inception cohort. Arthritis Rheum. 2006, 54: 1105-1116. 10.1002/art.21750.
Pawlik A, Florczak M, Ostanek L, Brzosko M, Brzosko I, Szklarz BG: TNF-alpha -308 promoter polymorphism in patients with rheumatoid arthritis. Scand J Rheumatol. 2005, 34: 22-26.
Rojas-Villarraga A, Diaz FJ, Calvo-Paramo E, Salazar JC, Iglesias-Gamarra A, Mantilla RD, Anaya JM: Familial disease, the HLA-DRB1 shared epitope and anti-CCP antibodies influence time at appearance of substantial joint damage in rheumatoid arthritis. J Autoimmun. 2009, 32: 64-69. 10.1016/j.jaut.2008.11.004.
Lacki JK, Moser R, Korczowska I, Mackiewicz S, Muller W: TNF-alpha gene polymorphism does not affect the clinical and radiological outcome of rheumatoid arthritis. Rheumatol Int. 2000, 19: 137-140. 10.1007/s002960050117.
Fonseca JE, Cavaleiro J, Teles J, Sousa E, Andreozzi VL, Antunes M, maral-Turkman MA, Canhao H, Mourao AF, Lopes J, Caetano-Lopes J, Weinmann P, Sobral M, Nero P, Saavedra MJ, Malcata A, Cruz M, Melo R, Braña A, Miranda L, Patto JV, Barcelos A, da Silva JC, Santos LM, Figueiredo G, Rodrigues M, Jesus H, Quintal A, Carvalho T, da Silva JA, et al: Contribution for new genetic markers of rheumatoid arthritis activity and severity: sequencing of the tumor necrosis factor-alpha gene promoter. Arthritis Res Ther. 2007, 9: R37-10.1186/ar2173.
Wilson AG, de VN, van de Putte LB, Duff GW: A tumour necrosis factor alpha polymorphism is not associated with rheumatoid arthritis. Ann Rheum Dis. 1995, 54: 601-603. 10.1136/ard.54.7.601.
Bouma G, Crusius JB, Oudkerk PM, Kolkman JJ, von Blomberg BM, Kostense PJ, Giphart MJ, Schreuder GM, Meuwissen SG, Pena AS: Secretion of tumour necrosis factor alpha and lymphotoxin alpha in relation to polymorphisms in the TNF genes and HLA-DR alleles. Relevance for inflammatory bowel disease. Scand J Immunol. 1996, 43: 456-463. 10.1046/j.1365-3083.1996.d01-65.x.
Hajeer AH, Hutchinson IV: Influence of TNFalpha gene polymorphisms on TNFalpha production and disease. Hum Immunol. 2001, 62: 1191-1199. 10.1016/S0198-8859(01)00322-6.
Abraham LJ, Kroeger KM: Impact of the -308 TNF promoter polymorphism on the transcriptional regulation of the TNF gene: relevance to disease. J Leukoc Biol. 1999, 66: 562-566.
Kroeger KM, Abraham LJ: Identification of an AP-2 element in the -323 to -285 region of the TNF-alpha gene. Biochem Mol Biol Int. 1996, 40: 43-51.
Kroeger KM, Carville KS, Abraham LJ: The -308 tumor necrosis factor-alpha promoter polymorphism effects transcription. Mol Immunol. 1997, 34: 391-399. 10.1016/S0161-5890(97)00052-7.
Braun N, Michel U, Ernst BP, Metzner R, Bitsch A, Weber F, Rieckmann P: Gene polymorphism at position -308 of the tumor-necrosis-factor-alpha (TNF-alpha) in multiple sclerosis and it's influence on the regulation of TNF-alpha production. Neurosci Lett. 1996, 215: 75-78.
Oregon-Romero E, Vazquez-Del MM, Ruiz-Quezada SL, Navarro-Hernandez RE, Rangel-Villalobos H, Martinez-Bonilla G, Bernard-Medina AG, rmendariz-Borunda J, Garcia-Banuelos J, Munoz-Valle JF: Tumor necrosis factor alpha-308 and -238 polymorphisms in rheumatoid arthritis. Association with messenger RNA expression and sTNF-alpha. J Investig Med. 2008, 56: 937-943.
Munoz-Valle JF, Oregon-Romero E, Rangel-Villalobos H, Martinez-Bonilla GE, Castaneda-Saucedo E, Salgado-Goytia L, Leyva-Vazquez MA, Illades-Aguiar B, arcon-Romero LD, Espinoza-Rojo M, Parra-Rojas I: High expression of TNF alpha is associated with -308 and -238 TNF alpha polymorphisms in knee osteoarthritis. Clin Exp Med. 2012
Helmig S, Aliahmadi N, Stephan P, Dohrel J, Schneider J: TNF-alpha -308 genotypes are associated with TNF-alpha and TGF-beta(1) mRNA expression in blood leucocytes of humans. Cytokine. 2011, 53: 306-310. 10.1016/j.cyto.2010.11.018.
Knight JC, Keating BJ, Rockett KA, Kwiatkowski DP: In vivo characterization of regulatory polymorphisms by allele-specific quantification of RNA polymerase loading. Nat Genet. 2003, 33: 469-475. 10.1038/ng1124.
Mekinian A, Tamouza R, Pavy S, Gestermann N, Ittah M, Mariette X, Miceli-Richard C: Functional study of TNF-alpha promoter polymorphisms: literature review and meta-analysis. Eur Cytokine Netw. 2011, 22: 88-102.
Welsing PM, van Riel PL: The Nijmegen inception cohort of early rheumatoid arthritis. J Rheumatol Suppl. 2004, 69: 14-21.
Rau R, Wassenberg S, Herborn G, Stucki G, Gebler A: A new method of scoring radiographic change in rheumatoid arthritis. J Rheumatol. 1998, 25: 2094-2107.
Miller SA, Dykes DD, Polesky HF: A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988, 16: 1215-10.1093/nar/16.3.1215.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
van den Berg WB, van Riel PL: Uncoupling of inflammation and destruction in rheumatoid arthritis: myth or reality?. Arthritis Rheum. 2005, 52: 995-999. 10.1002/art.20981.
Jawaheer D, Li W, Graham RR, Chen W, Damle A, Xiao X, Monteiro J, Khalili H, Lee A, Lundsten R, Begovich A, Bugawan T, Erlich H, Elder JT, Criswell LA, Seldin MF, Amos CI, Behrens TW, Gregersen PK: Dissecting the genetic complexity of the association between human leukocyte antigens and rheumatoid arthritis. Am J Hum Genet. 2002, 71: 585-594. 10.1086/342407.
We thank all patients for their participation in this study and Caesar Roseboom for his work concerning the genotyping. This work was supported by a personal grant to Marieke Coenen from the Netherlands Organization for Scientific Research (grant 916.76.020).
The authors declare that they have no competing interests.
ET collected DNA and RNA samples from RA patients, carried out and interpreted the results of the molecular genetic studies, participated in the study design and wrote the manuscript. PB participated in the study design and evaluated the manuscript. JF carried out the longitudinal regression analysis and evaluated the manuscript. AdB helped perform the molecular genetic studies, interpreted the results and evaluated the manuscript. AE was involved in collecting DNA and RNA samples from RA patients and evaluated the manuscript. PM helped in collecting data for the meta-analysis and evaluated the manuscript. HM helped in collecting data for the meta-analysis and evaluated the manuscript. HS interpreted the results and evaluated the manuscript. PvR participated in the study design, interpreted the results and evaluated the manuscript. BF participated in the study design, was involved in interpreting the results and evaluated the manuscript. MC was involved in collecting DNA and RNA samples from RA patients, participated in the study design, collected data for the meta-analysis, interpreted the results and wrote the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Toonen, E.J., Barrera, P., Fransen, J. et al. Meta-analysis identified the TNFA -308G > A promoter polymorphism as a risk factor for disease severity in patients with rheumatoid arthritis. Arthritis Res Ther 14, R264 (2012). https://doi.org/10.1186/ar4110