Variations in the TNFα gene and their interactions with the IL4R and IL10 genes in relation to hand osteoarthritis
BMC Musculoskeletal Disorders volume 15, Article number: 311 (2014)
The development of osteoarthritis (OA) involves inflammation, but the evidence for participation of genes propagating or inhibiting inflammation in the OA process is inconsistent. We investigated the associations of common variants in the TNFα gene, and their interactions with other cytokine genes, with hand OA among Finnish women.
This cross-sectional study was based on bilateral hand radiographs of 542 female dentists and teachers which were classified according to the presence of OA (radiographic K-L score ≥ 2 in ≥ 3 joints) using reference images. The genotypes were determined by PCR-based methods. The degree of pairwise linkage disequilibrium (LD) and haplotypes were constructed and analyzed by the SNPStats software. The associations between four TNFα SNPs and hand OA were tested using logistic regression adjusting for age, occupation, and BMI, and fitting a log-additive model of inheritance. Gene-gene interactions of TNFα SNPs with IL4R and IL10 SNPs were examined by stratified logistic regression analyses. Possible interactions of the TNFα SNPs with variants in the previously reported IL1β and IL6 genes in influencing hand OA were also explored.
Two TNFα polymorphisms (“-1031” and “-863”) were associated with hand OA (OR = 1.45, 95% CI 1.01-2.07 and 1.55, 1.06-2.25, respectively). These associations retained when adjusting further for IL1β “3954” and IL6 “174”. The TNFα G-A-G haplotype was associated with an increased risk of hand OA (1.61, 1.10-2.37, p = 0.01). Interactions were observed between TNFα “-1031” and IL4R Ser503Pro, TNFα “-1031” and IL10 “-1082”, and TNFα “-863” and IL10 “-1082” SNPs with regard to hand OA (p = 0.012, p = 0.0068, and p = 0.02, respectively). The carriage of the TNFα “-1031” minor allele doubled the risk (2.01, 1.26 - 3.22) only in women with the IL4R Ser/Ser genotype. Similarly, the TNFα “-1031” and “-863” minor alleles were associated with an increased risk of hand OA only in IL10 G/G or A/A homozygotes (2.54, 1.45-4.47 and 2.60, 1.46-4.62, respectively) but not in heterozygotes (G/A).
Our results suggest that the TNF α gene variants play a role in the etiology of hand OA. In addition, the findings are suggestive of a gene-gene interaction of the TNF α with IL4R and IL10 genes.
Osteoarthritis (OA) is the most common joint disorder worldwide and rapidly increasing with ageing populations. OA is a dynamic process involving all the structures within the joint, i.e., cartilage, synovial membrane and subchondral bone. It shows clinical heterogeneity in joint numbers and regions involved. Some patients may have only one site affected (hip, knee, or hand) affected (local OA), while others have clustered joint regions affected in a characteristic distribution (generalized OA) .
The hand is among the most frequently affected site in OA . The prevalence of hand OA is higher in women than in men over the age of 50 . Although the pathogenesis of hand OA is largely unknown, familial aggregation and heritability studies indicate a significant genetic role in addition to the involvement of mechanical (repetitive joint loading) and lifestyle related factors (e.g. obesity) [3, 4].
The development and progression of OA are nowadays believed to involve inflammation [5–8]. Chondrocytes, as well as synovial cells, of OA patients produce increased levels of pro-inflammatory cytokines, which affect metabolism and enhance the catabolism of all joint tissues affected in OA . Among pro-inflammatory cytokines, interleukin-1β (IL1β) and tumor necrosis factor alpha (TNFα) seem prominent and of major importance to cartilage destruction as they are synthesized during the OA process [9–11]. In vivo studies have shown that these cytokines can act independently or in concert with other cytokines (e.g. IL6) in the induction and propagation of inflammation . Synthesis of the IL1β and TNFα is inhibited by anti-inflammatory cytokines such as IL4, IL10 and IL13 . On the other hand, expression of cytokine genes like IL1 and TNFα is up-regulated in OA .
The gene encoding TNFα is located in the class III region of the major histocompatibility complex (MHC) which is the most gene-dense and polymorphic region of the entire genome . TNFα is driving the inflammatory cascade . IL4 in turn is a cytokine produced by T cells, which plays a major role in immunoglobulin E (IgE) production. Its signals are conferred to effector cells through binding to the alpha chain of IL4 receptor (IL4R). IL4 and IL4R are expressed by human articular chondrocytes; data suggest that mechanical stimulation induces the release of IL4 by human chondrocytes after the recognition and transduction of the mechanical signal by integrin [16, 17]. Therefore, the IL4R is an active autocrine or paracrine signaling molecule in a regulatory pathway in the maintenance of human articular cartilage structure and function . Regulation of the structure and function of human articular cartilage occurs by mediating other biochemical responses to mechanical strain, proteoglycan synthesis, or altering the expressions of other extra cellular matrix (ECM) proteins involved in the pathogenesis of OA [16, 17].
So far, the evidence for involvement of genes propagating or inhibiting inflammation in the development or progression of OA is inconsistent, and the observed associations were not replicated in an independent population . Most of the previous studies examined the role of a single gene in OA without taking into consideration the interaction of the genes participating in the regulation of balance between pro- and anti-inflammatory processes. Our group has reported the associations of the IL1 extended haplotype and common IL6 promoter variants with symmetrical DIP OA [19, 20].
The aim of the current study was to investigate the associations of common variants in the TNFα gene and their interactions with variants in the IL4R and IL10 genes in relation to hand OA among middle-aged Finnish women, representing two occupations: dentists and teachers. The possible interactions of the TNFα with variants in the IL1β and IL6 in influencing hand OA were also explored.
Study design and participants’ selection
This was a cross-sectional study, samples of which were taken randomly from two occupational groups. The study participants were identified from the registers of the Finnish Dental Association and the Finnish Teachers Trade Union. Four hundred and thirty-six women aged 45 to 63 were randomly selected from both occupational groups (altogether 872 subjects) by using the place of residence (Helsinki or its neighboring cities) as an inclusion criterion. Of those subjects who received the questionnaires in 2002, 542 (62% of the invited) participated in a clinical examination between October 2002 and March 2003. Of these, 294 (67% of the invited) were dentists and 248 (57% of the invited) teachers. Participation in the study was voluntary and based on informed consent. The study was approved by the Hospital District of Helsinki and Uusimaa Ethics Committee for Research in Occupational Health and Safety.
Hand radiography and image analysis
Both hands of the study participants were radiographed by exposing Kodak X-ray films with Siemens X-ray equipment (48 kV, 10 mA, focus film distance = 115 cm; Siemens, Munich, Germany). The analogue radiographs were evaluated by an experienced radiologist who was blinded to the occupation, age, and all health data of the participants. Each distal interphalangeal (DIP), proximal interphalangeal (PIP), and thumb interphalangeal (IP) joint of both hands was graded separately, and classified for the presence of OA using a modified Kellgren and Lawrence (K-L) system ; the classification criteria were: grade 0 = no OA, grade 1 = doubtful OA, grade 2 = mild OA, grade 3 = moderate OA, and grade 4 = severe OA. The description of reference images used in the classification is given elsewhere . The reliability of the readings was estimated by measuring intra-observer and inter-observer agreements (intraclass correlation) within a limited sample of radiographs and a second participating radiologist. The inter-observer agreement for OA ranged from 0.67 to 0.85 for DIP joints and from 0.39 to 0.61 for PIP joints. The intra-observer agreement for OA ranged from 0.73 to 0.88 for DIP joints and from 0.67 to 0.92 for PIP joints .
Participants who had at least three finger joints with radiographic OA of grade 2 to 4 were classified as having hand OA. Otherwise, the participants were classified as not having hand OA.
Weight was measured without shoes to the accuracy of 0.1 kg. Body mass index [BMI = weight (kg)/height (m)2] was calculated based on weight and self-reported height. BMI data was missing from one participant. Age, occupation, and BMI were considered as possible confounders in the analyses. The variants in the IL1β and IL6 that was previously shown to influence hand OA were also included among the covariates.
Blood samples were taken from each study participant in the clinical examination and stored at +4°C until DNA extraction using a DNA extraction kit (PUREGENE ®DNA Purification Kit; Gentra Systems, Plymouth, MN, USA).
The TNFα “-1031” (rs1799964) and the “-857” (rs1799724) genotypes were determined by the TaqMan® SNP Genotyping Assay (Applied Biosystems, C___7514871_10 and C__11918223_10, respectively).
The TNFα “-863” (rs1800630) genotype was determined and the “-857” genotype re-determined by the Pyrosequencing® PSQ 96MA SNP/SQA system with PyroMark Assay Design self-designed protocol.
The TNFα “-308” (rs1800629) genotype was determined by PCR-RFLP method with Nco I (New England BioLabs (NEB) 10 U/μL) restriction enzyme. The primers were from Ozen et al..
In the TNFα “-1031” locus the T-allele was denoted as the wild type allele and the C-allele as the variant allele, in the “-863” and “-857” loci the C-alleles were denoted as the wild type alleles and the A- and T-allele as the variant alleles, respectively, and in the “-308” locus the G-allele was denoted as the wild type allele and the A-allele as the variant allele.
The IL4R Ser503Pro (1507 T > C, rs1805015) and Ser752Ala (2254 T > G, rs1805016) polymorphisms were genotyped by PCR-based TaqMan® SNP Genotyping Assays (Applied Biosystems, C_234284_1 and C_8903091_10 respectively).
The IL10 “-1082” (rs1800896) SNP was genotyped with primers from Koch et al. .
An additional file has detailed description about the genotyping (see Additional file 1).
For quality control 10% of the genotyped samples were blindly repeated with 100% concordant results. Genotype data was available from all participants.
The potential deviation of the allele frequencies from the Hardy-Weinberg equilibrium (HWE) was tested from controls using the chi-square test. The degree of pairwise linkage disequilibrium (LD) for four TNFα SNPs and two IL4R SNPs were calculated using SNPStats software . Haplotypes were constructed and analyzed by the same software.
Logistic regression analysis was used to test the associations between SNPs and hand OA. For each SNP, a log-additive model of inheritance was fitted. To evaluate whether the observed association between the TNFα and OA was modified by variants in other cytokine genes, gene-gene interactions were tested for all TNFα SNPs by stratified logistic regression analyses.
Both crude and adjusted odds ratios (ORs) and their 95% confidence intervals (CIs) were calculated. The ORs were adjusted for the potential confounding factors, i.e., age (continuous), occupation (dentists vs. teachers), and BMI (continuous). Since the crude and adjusted ORs did not differ significantly, only the adjusted ORs are shown in the results. In addition, the OR’s were further adjusted for genetic variants in the IL1β and IL6.
In addition to exploring whether the effect of the TNFα SNPs on hand OA were independent of the genetic variants in the IL1β and IL6, we estimated the individual and joint effects of the TNFα (“-863”), IL1β and IL6 polymorphisms using the combinations of two dummy (0, 1) variables. First, we calculated the sum of the minor alleles of IL1β and IL6, by summing up the number of minor alleles of two SNPs. This was dichotomized (first dummy variable): 0 = non-carriers of any minor allele of the IL1β and IL6 and 1 = carriers of at least one minor allele. For the TNFα “-863” SNP we used the dominant model, with the homozygous genotype of the major allele as the reference (second dummy variable).
All analyses were hypothesis driven. The statistical significance of the p-value was defined as the 1% level. P-values were adjusted for multiple testing using Sidák’s method . We used SNPStats software  and SPSS 20.0 for the analyses.
The prevalence of hand OA with at least three affected finger joints was 29.5%, being higher among teachers (35.5%) than dentists (24.5%) (Table 1). Participants with hand OA were significantly older and had higher BMI than those without OA.
The genotype frequencies were in HWE in all of the studied polymorphic loci (Table 2). When adjusted for age, occupation and BMI, two TNFα SNPs (“-1031” and “-863”) were associated with hand OA (OR = 1.45, 95% CI 1.01-2.07, p = 0.04 and 1.55, 1.06-2.25, p = 0.02, respectively) (Table 3). Further adjustment for the IL1β and IL6 SNPs had a negligible effect on the observed point estimates, though improving the estimate’s precision (p = 0.03 and p = 0.01, respectively). No statistically significant associations were found between the other two TNFα SNPs and hand OA. Neither were there associations between the SNPs in the IL4R or IL10 and hand OA.
Statistically significant interactions were found between the TNFα “-1031” and IL4R Ser503Pro SNPs, TNFα “-1031” and IL10 “-1082” SNPs, and TNFα “-863” and IL10 “-1082” SNPs and hand OA (p = 0.012, p = 0.0068, and p = 0.02, respectively). The carriage of the TNFα (“-1031”) minor allele was associated with a double risk of hand OA (2.01, 1.26 - 3.22) in women with the IL4R Ser/Ser genotype (Table 4). Similarly, the TNFα “-1031” and “-863” minor alleles were associated with an increased risk of hand OA only in IL10 G/G or A/A homozygotes (2.54, 1.45-4.47 and 2.60, 1.46-4.62, respectively) but not in heterozygotes (G/A).
We also examined the individual and joint effects of the TNFα “-863”, IL1β, and IL6 polymorphisms on hand OA. The risk of hand OA was the highest in the carriers of the minor alleles in all three genes (4.37, 1.84-10.38, p = 0.001). Somewhat lower risks were observed for carriers of the TNFα “-863” minor allele (3.73, 1.28-10.85, p = 0.016) and carriers of minor alleles of IL1β and IL6 SNPs (2.89, 1.29-6.48, p = 0.010).
The degree of pairwise LD between three TNFα SNPs (“-1031”, “-863”, and “-308”) was high; the “-1031” and “-863” SNPs were in complete linkage (D’ = 1, r2 = 0.93, p < 0.0001), and the IL4R Ser503Pro and Ser752Ala polymorphisms were also in a strong LD (D’ = 0.998, r2 = 0.48, p < 0.0001). The three TNFα promoter polymorphisms composed a total of four haplotypes. The most common of these haplotypes was T-C-G (0.69), followed by G-A-G (0.16), T-C-A (0.13), and G-C-G (0.02). The two IL4R polymorphisms, on the other hand, composed three haplotypes, i.e., Ser-Ser (0.87), Ser-Pro (0.10) and Ala-Pro (0.03).
The TNFα G-A-G haplotype was associated with an increased risk of hand OA when adjusted for age, occupation, and BMI (1.61, 1.10-2.37, p = 0.01) (Table 5). There was no difference between participants with and without hand OA in the IL4R haplotype distribution.
We investigated whether the TNF α promoter polymorphisms are associated with hand OA among Finnish women. The minor alleles of the TNFα “-1031” and “-863” loci, as well as their haplotype, were found to be associated with an increased risk of hand OA. The observed associations were independent of the variants in the IL1β and IL6 genes. Furthermore, our findings suggest that the effect of TNFα polymorphisms on hand OA is modified by the variants within the IL4R and IL10 genes.
A traditional paradigm of OA as a “wear and tear” disease leading to the loss of cartilage has been revised. Nowadays, OA is considered a complex disease with inflammatory mediators released by cartilage, bone and synovium . TNFα is one of the most typical pro-inflammatory cytokines that along with IL1β is connected with cartilage destruction. These two cytokines, which are produced by chondrocytes, mononuclear cells, osteoblasts and synovial tissues, induce the production of a number of inflammatory and catabolic factors . Among the four polymorphic loci studied here, only the TNFα “-308” locus has been shown to affect the TNFα protein levels, the minor allele of the SNP was associated with increased TNFα production in response to various stimuli [27, 28]. However, also opposite observations, e.g., no effect on the protein levels or lowered protein levels, have been reported [29, 30]. When studying the above SNPs with F-SNP-program that is freely available on the internet (http://compbio.cs.queensu.ca/F-SNP/) connected to the main databases, and computationally predicting functional SNPs, all four SNPs are predicted to be functional as they seem to be in the transcription factor binding site . However, the protein level alteration by the studied SNPs still remains unsolved and needs to be further studied.
The few studies that examined associations between TNFα polymorphisms and knee or hip OA have given conflicting results . To our knowledge, this is the first study to report an association between the TNFα variant alleles and hand OA. TNFα can act independently or in concert with other cytokines (e.g., IL1β, IL4, IL6, and Il10) to initiate and expand inflammation . Ignorance of the complex interrelationships between pro- and anti-inflammatory cytokines might be the reason for failure to detect an association between the TNFα polymorphisms and OA. It has been suggested that the combined use of information from multiple markers may be more effective to reveal the association between a genomic region and a trait than a single marker analysis .
Previously our group reported the associations of the IL1 and IL6 gene polymorphisms with hand OA [19, 20]. The current findings suggest that TNFα promoter polymorphisms may increase the risk of hand OA independently of the polymorphisms in the IL1 and IL6 genes, and that the effect attributed to combination of variants in all three genes is larger, but less than additive.
When we examined the association of the TNFα polymorphisms with hand OA taking into consideration the variants of other cytokine genes, the IL4R and IL10 polymorphisms appeared to act as effect modifiers.
Vargiolu et al.  reported an association of genetic variants in the coding region of the IL4R gene with hand OA among men and women in the age range of 41 to 84 years. However, we failed to replicate this association among our participants. Differences between the study populations, OA phenotypes, and minor allele frequencies might be the reasons for discrepancies in the findings. Our participants were younger (mean age 53 years) than in the study by Vargiolu and coworkers, which may partly explain the difference of the prevalence of hand OA between their and our study (55.6% and 29.5%, respectively).
Naturally occurring anti-inflammatory cytokines such as IL10 inhibit the synthesis of IL1 and TNFα . The IL10 “-1082” polymorphism has been shown to affect the level of the protein production: the A-allele is connected with a significantly higher protein production than the G-allele . This polymorphism was also associated with rheumatoid arthritis in a meta-analysis . As to its role in DIP OA, no association was found in the Dutch population . Similarly, the IL10 “-1082” did not associate with hand OA in our study.
A major strength of the study was that all study participants were of the ethnically relatively homogenous Finnish origin. Each ethnic group has its own set of environmental and genetic factors that contribute to the disease risk, and differences in allelic frequency often affect our ability to detect a susceptibility allele. The Finnish population is known to be a genetic isolate, which originated from a small founder population some 2000 years ago. Therefore, the Finnish population with the relatively homogenous gene pool  offers an optimal material for association studies.
Another strength of our study is that we analyzed haplotypes in addition to SNPs. Grouping of SNPs in haplotypes generally leads to a stronger association with the phenotype than individual polymorphisms.
A limitation of our study is the relatively small number of participants, leading to reduced power to detect small effects and an increased likelihood of spurious findings. This needs to be considered while interpreting the observed associations. Another obvious limitation is the fact that our study participants were all women and consequently the results cannot be generalized to men.
Our results suggest that variants in the TNF α gene play a role in the etiology of hand OA in Finnish women. In addition, the findings are suggestive of a gene-gene interaction of the TNF α with the IL4R and IL10 genes. However, these findings should be considered with caution until replicated in other study population.
Tumor necrosis factor alpha
Body mass index
Polymerase chain reaction
Single nucleotide polymorphism
Interleukin 1, beta
Kellgren and Lawrence
Restriction fragment length polymorphism
Major histocompatibility complex.
Bijlsma JW, Berenbaum F, Lafeber FP: Osteoarthritis: an update with relevance for clinical practice. Lancet. 2011, 377 (9783): 2115-2126. 10.1016/S0140-6736(11)60243-2.
Loeser RF: Age-related changes in the musculoskeletal system and the development of osteoarthritis. Clin Geriatr Med. 2010, 26 (3): 371-386. 10.1016/j.cger.2010.03.002.
Kalichman L, Hernandez-Molina G: Hand osteoarthritis: an epidemiological perspective. Semin Arthritis Rheum. 2010, 39 (6): 465-476. 10.1016/j.semarthrit.2009.03.001.
Arden N, Nevitt MC: Osteoarthritis: epidemiology. Best Pract Res Clin Rheumatol. 2006, 20 (1): 3-25. 10.1016/j.berh.2005.09.007.
Kapoor M, Martel-Pelletier J, Lajeunesse D, Pelletier JP, Fahmi H: Role of proinflammatory cytokines in the pathophysiology of osteoarthritis. Nat Rev Rheumatol. 2011, 7 (1): 33-42. 10.1038/nrrheum.2010.196.
Farahat MN, Yanni G, Poston R, Panayi GS: Cytokine expression in synovial membranes of patients with rheumatoid arthritis and osteoarthritis. Ann Rheum Dis. 1993, 52 (12): 870-875. 10.1136/ard.52.12.870.
Sellam J, Berenbaum F: The role of synovitis in pathophysiology and clinical symptoms of osteoarthritis. Nat Rev Rheumatol. 2010, 6 (11): 625-635. 10.1038/nrrheum.2010.159.
Berenbaum F: Osteoarthritis as an inflammatory disease (osteoarthritis is not osteoarthrosis!). Osteoarthritis Cartilage. 2013, 21 (1): 16-21. 10.1016/j.joca.2012.11.012.
Krasnokutsky S, Attur M, Palmer G, Samuels J, Abramson SB: Current concepts in the pathogenesis of osteoarthritis. Osteoarthritis Cartilage. 2008, 16 (Suppl 3): S1-S3.
Pelletier JP, Martel-Pelletier J, Abramson SB: Osteoarthritis, an inflammatory disease: potential implication for the selection of new therapeutic targets. Arthritis Rheum. 2001, 44 (6): 1237-1247. 10.1002/1529-0131(200106)44:6<1237::AID-ART214>3.0.CO;2-F.
Martel-Pelletier J, Alaaeddine N, Pelletier JP: Cytokines and their role in the pathophysiology of osteoarthritis. Front Biosci. 1999, 4: D694-D703. 10.2741/Martel.
Page Thomas DP, King B, Stephens T, Dingle JT: In vivo studies of cartilage regeneration after damage induced by catabolin/interleukin-1. Ann Rheum Dis. 1991, 50 (2): 75-80. 10.1136/ard.50.2.75.
Fernandes JC, Martel-Pelletier J, Pelletier JP: The role of cytokines in osteoarthritis pathophysiology. Biorheology. 2002, 39 (1–2): 237-246.
Malemud CJ, Islam N, Haqqi TM: Pathophysiological mechanisms in osteoarthritis lead to novel therapeutic strategies. Cells Tissues Organs. 2003, 174 (1–2): 34-48.
Hajeer AH, Hutchinson IV: TNF-alpha gene polymorphism: clinical and biological implications. Microsc Res Tech. 2000, 50 (3): 216-228. 10.1002/1097-0029(20000801)50:3<216::AID-JEMT5>3.0.CO;2-Q.
Millward-Sadler SJ, Wright MO, Lee H, Nishida K, Caldwell H, Nuki G, Salter DM: Integrin-regulated secretion of interleukin 4: A novel pathway of mechanotransduction in human articular chondrocytes. J Cell Biol. 1999, 145 (1): 183-189. 10.1083/jcb.145.1.183.
Salter DM, Millward-Sadler SJ, Nuki G, Wright MO: Integrin-interleukin-4 mechanotransduction pathways in human chondrocytes. Clin Orthop Relat Res. 2001, Oct (391 Suppl): S49-S60.
Ryder JJ, Garrison K, Song F, Hooper L, Skinner J, Loke Y, Loughlin J, Higgins JP, MacGregor AJ: Genetic associations in peripheral joint osteoarthritis and spinal degenerative disease: a systematic review. Ann Rheum Dis. 2008, 67 (5): 584-591.
Solovieva S, Kamarainen OP, Hirvonen A, Hamalainen S, Laitala M, Vehmas T, Luoma K, Nakki A, Riihimaki H, Ala-Kokko L, Mannikko M, Leino-Arjas P: Association between interleukin 1 gene cluster polymorphisms and bilateral distal interphalangeal osteoarthritis. J Rheumatol. 2009, 36 (9): 1977-1986. 10.3899/jrheum.081238.
Kamarainen OP, Solovieva S, Vehmas T, Luoma K, Riihimaki H, Ala-Kokko L, Mannikko M, Leino-Arjas P: Common interleukin-6 promoter variants associate with the more severe forms of distal interphalangeal osteoarthritis. Arthritis Res Ther. 2008, 10 (1): R21-10.1186/ar2374.
Kellgren JH, Lawrence JS: Radiological assessment of osteo-arthrosis. Ann Rheum Dis. 1957, 16 (4): 494-502. 10.1136/ard.16.4.494.
Solovieva S, Vehmas T, Riihimaki H, Luoma K, Leino-Arjas P: Hand use and patterns of joint involvement in osteoarthritis. A comparison of female dentists and teachers. Rheumatology (Oxford). 2005, 44: 521-528. 10.1093/rheumatology/keh534.
Ozen S, Alikasifoglu M, Bakkaloglu A, Duzova A, Jarosova K, Nemcova D, Besbas N, Vencovsky J, Tuncbilek E: Tumour necrosis factor alpha G → A 238 and G → A 308 polymorphisms in juvenile idiopathic arthritis. Rheumatology (Oxford). 2002, 41 (2): 223-227. 10.1093/rheumatology/41.2.223.
Koch W, Kastrati A, Bottiger C, Mehilli J, von Beckerath N, Schomig A: Interleukin-10 and tumor necrosis factor gene polymorphisms and risk of coronary artery disease and myocardial infarction. Atherosclerosis. 2001, 159 (1): 137-144. 10.1016/S0021-9150(01)00467-1.
Sole X, Guino E, Valls J, Iniesta R, Moreno V: SNPStats: a web tool for the analysis of association studies. Bioinformatics. 2006, 22 (15): 1928-1929. 10.1093/bioinformatics/btl268.
Sidak Z: Rectangular Confidence Regions for the Means of Multivariate Normal Distributions. J Am Stat Assoc. 1967, 62 (318): 626-633.
Higuchi T, Seki N, Kamizono S, Yamada A, Kimura A, Kato H, Itoh K: Polymorphism of the 5′-flanking region of the human tumor necrosis factor (TNF)-alpha gene in Japanese. Tissue Antigens. 1998, 51 (6): 605-612. 10.1111/j.1399-0039.1998.tb03002.x.
Wilson AG, Symons JA, McDowell TL, McDevitt HO, Duff GW: Effects of a polymorphism in the human tumor necrosis factor alpha promoter on transcriptional activation. Proc Natl Acad Sci U S A. 1997, 94 (7): 3195-3199. 10.1073/pnas.94.7.3195.
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 (2): 88-102.
Sharma S, Sharma A, Kumar S, Sharma SK, Ghosh B: Association of TNF haplotypes with asthma, serum IgE levels, and correlation with serum TNF-alpha levels. Am J Respir Cell Mol Biol. 2006, 35 (4): 488-495. 10.1165/rcmb.2006-0084OC.
Lee PH, Shatkay H: F-SNP: computationally predicted functional SNPs for disease association studies. Nucleic Acids Res. 2008, 36 (Database issue): D820-D824.
Han L, Song JH, Yoon JH, Park YG, Lee SW, Choi YJ, Nam SW, Lee JY, Park WS: TNF-alpha and TNF-beta Polymorphisms are Associated with Susceptibility to Osteoarthritis in a Korean Population. Korean J Pathol. 2012, 46 (1): 30-37. 10.4132/KoreanJPathol.2012.46.1.30.
Gauderman WJ, Murcray C, Gilliland F, Conti DV: Testing association between disease and multiple SNPs in a candidate gene. Genet Epidemiol. 2007, 31 (5): 383-395. 10.1002/gepi.20219.
Vargiolu M, Silvestri T, Bonora E, Dolzani P, Pulsatelli L, Addimanda O, Mancarella L, Punzi L, Fioravanti A, Facchini A, Romeo G, Meliconi R: Interleukin-4/interleukin-4 receptor gene polymorphisms in hand osteoarthritis. Osteoarthritis Cartilage. 2010, 18 (6): 810-816. 10.1016/j.joca.2010.02.005.
Turner DM, Williams DM, Sankaran D, Lazarus M, Sinnott PJ, Hutchinson IV: An investigation of polymorphism in the interleukin-10 gene promoter. Eur J Immunogenet. 1997, 24 (1): 1-8. 10.1111/j.1365-2370.1997.tb00001.x.
Zhang J, Zhang Y, Jin J, Li M, Xie K, Wen C, Cheng R, Chen C, Lu J: The -1082A/G polymorphism in the Interleukin-10 gene and the risk of rheumatoid arthritis: a meta-analysis. Cytokine. 2011, 56 (2): 351-355. 10.1016/j.cyto.2011.05.022.
Riyazi N, Kurreeman FA, Huizinga TW, Dekker FW, Stoeken-Rijsbergen G, Kloppenburg M: The role of interleukin 10 promoter polymorphisms in the susceptibility of distal interphalangeal osteoarthritis. J Rheumatol. 2005, 32 (8): 1571-1575.
Peltonen L, Jalanko A, Varilo T: Molecular genetics of the Finnish disease heritage. Hum Mol Genet. 1999, 8 (10): 1913-1923. 10.1093/hmg/8.10.1913.
van Saase JL, van Romunde LK, Cats A, Vandenbroucke JP, Valkenburg HA: Epidemiology of osteoarthritis: Zoetermeer survey. Comparison of radiological osteoarthritis in a Dutch population with that in 10 other populations. Ann Rheum Dis. 1989, 48 (4): 271-280. 10.1136/ard.48.4.271.
Sowers M, Lachance L, Hochberg M, Jamadar D: Radiographically defined osteoarthritis of the hand and knee in young and middle-aged African American and Caucasian women. Osteoarthritis Cartilage. 2000, 8 (2): 69-77. 10.1053/joca.1999.0273.
Haara MM, Manninen P, Kroger H, Arokoski JP, Karkkainen A, Knekt P, Aromaa A, Heliovaara M: Osteoarthritis of finger joints in Finns aged 30 or over: prevalence, determinants, and association with mortality. Ann Rheum Dis. 2003, 62 (2): 151-158. 10.1136/ard.62.2.151.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2474/15/311/prepub
We are grateful to Katariina Luoma, MD (University of Helsinki, Helsinki University Central Hospital, Radiology department, Finland), for performing the second readings of the hand radiographs for reliability analysis. We are also grateful to Sirpa Hyttinen and Mari Kukkonen for referring the genotyping results and the input of the data to the database, and Susanna Lemmelä for valuable advice in the preparation of the manuscript. The study was financially supported by a grant from the Finnish Work Environment Fund (101334) and grant from National Doctoral Programme of Musculoskeletal Disorders and Biomaterials (TBDP). The study sponsors had no role in the study design, collection, analysis and interpretation of data, in the writing of the manuscript, or in the decision to submit the manuscript for publication.
The authors declare that they have no competing interests.
SH carried the main responsibility of the design and performance of the genotyping and data analyses, and preparation of the manuscript; SS participated in the study design, data collection, the design of the data analyses, the interpretation of the results, and preparation of the manuscript; TV participated in the data collection, radiological examinations, data analysis, and preparation of the manuscript; PL-A carried the main responsibility of the overall study design, and participated in the data collection and preparation of the manuscript; AH participated in the conception and design of the study, the interpretation of the data, and preparation of the manuscript. All authors have read and approved the final version of the manuscript.
Electronic supplementary material
About this article
Cite this article
Hämäläinen, S., Solovieva, S., Vehmas, T. et al. Variations in the TNFα gene and their interactions with the IL4R and IL10 genes in relation to hand osteoarthritis. BMC Musculoskelet Disord 15, 311 (2014). https://doi.org/10.1186/1471-2474-15-311