Reasearch Awards nomination

Email updates

Keep up to date with the latest news and content from Arthritis Research & Therapy and BioMed Central.

Open Access Research article

Gene-environment interaction between body mass index and transforming growth factor beta 1 (TGFβ1) gene in knee and hip osteoarthritis

Stella G Muthuri1, Sally Doherty2, Weiya Zhang2, Rose A Maciewicz3, Kenneth R Muir4 and Michael Doherty2*

Author affiliations

1 Division of Epidemiology and Public Health, University of Nottingham, Clinical Sciences Building, Nottingham City Hospital, Nottingham, NG5 1PB, UK

2 Academic Rheumatology, University of Nottingham, Clinical Sciences Building, Nottingham City Hospital, Nottingham, NG5 1PB, UK

3 Respiratory & Inflammation, AstraZeneca, Charnwood R&D, Loughborough, Leicestershire, LE11 5RH, UK

4 Heath Services Research Institute, Warwick Medical School and Warwick University, Coventry CV4 7AL, UK

For all author emails, please log on.

Citation and License

Arthritis Research & Therapy 2013, 15:R52  doi:10.1186/ar4214


The electronic version of this article is the complete one and can be found online at: http://arthritis-research.com/content/15/2/R52


Received:22 May 2012
Revisions received:1 March 2013
Accepted:5 April 2013
Published:18 April 2013

© 2013 Muthuri et al.; licensee 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.

Abstract

Introduction

The objective was to investigate potential gene-environment interaction between body mass index (BMI) and each of eight TGFβ1 polymorphisms in knee and hip osteoarthritis (OA).

Methods

We conducted a case-control study of Caucasian men and women aged 45 to 86 years from Nottingham, United Kingdom (Genetics of OA and Lifestyle (GOAL) study). Cases had clinically severe symptoms and radiographic knee or hip OA; controls had no symptoms and no radiographic knee/hip OA. We used logistic regression to investigate the association of TGFβ1 polymorphisms and OA when stratifying by BMI. Knee and hip OA were analyzed separately with adjustment for potential confounders. Additive and multiplicative interactions were examined.

Results

2,048 cases (1,042 knee OA, 1,006 hip OA) and 967 controls were studied. For hip OA, the highest risk was in overweight (BMI ≥25 kg/m2) individuals with the variant allele of single-nucleotide polymorphism (SNP) rs1800468 (odds ratio (OR) 2.21, 95% confidence interval (CI) 1.55, 3.15). Evaluation of gene-environment interaction indicated significant synergetic interaction (relative excess risk due to interaction (RERI) = 0.93, synergy index (SI) = 4.33) with an attributable proportion due to interaction (AP) of 42% (AP = 0.42; 95% CI 0.16, 0.68). Multiplicative interaction was also significant (OR for interaction (ORINT) = 2.27, P = 0.015). For knee OA, the highest risk was in overweight individuals with homozygous genotype 11 of SNP rs2278422 (OR = 6.95, P <0.001). In contrast, the variant allele indicated slightly lower risks (OR = 4.72, P <0.001), a significant antagonistic interaction (RERI = -2.66, SI = 0.59), AP = -0.56 (95%CI -0.94, -0.17) and a significant multiplicative interaction (ORINT = 0.47, P = 0.013).

Conclusion

TGFβ1 gene polymorphisms interact with being overweight to influence the risk of large joint OA.

Introduction

Osteoarthritis (OA) is a complex disorder with genetic and environmental risk factors both contributing to its development and progression [1]. Obesity is one of the strongest environmental risk factors for knee OA [2] and is considered to be moderately associated with hip OA [3]. Several chromosomal loci and gene variations have also been reported to influence OA disease processes [4,5]. Transforming growth factor beta 1 (TGFβ1) gene is one potential candidate gene for OA. TGFβ1 is considered essential for cartilage integrity and is found at high levels in normal cartilage but has reduced expression in cartilage in OA [6]. Furthermore, increased endogenous TGFβ1 expression is reported during osteophyte formation [6,7], which reportedly stabilise joints in OA and reduce structural progression [8].

Gene association studies have reported independent associations with peripheral OA and spinal osteophytosis. Within the TGFβ1 gene, single nucleotide polymorphisms (SNPs), rs2278422 and rs8179181, have been found to have a possible role in susceptibility to knee and hip OA in a British Caucasian population [9], and a variation on position 29 (Leu10Pro) or SNP rs1982073 has been implicated with spinal osteophytosis in Japanese postmenopausal women [10]. Furthermore, increased TGFβ1 expression has been associated with body mass index (BMI) and abdominal adipose tissue in morbid obesity [11]. The mechanism by which body weight may exert its effect on TGFβ1 is unclear. Evidence from in vitro experiments with human and animal culture cells as well as in vivo animal studies have shown an increased expression of TGFβ1 in response to mechanical stimuli in a number of cell and tissue types [12,13], implying that excessive biomechanical loading of joints may stimulate the up-regulation of endogenous TGFβ1 and thus, predispose to the development or progression of OA.

Environmental and genetic factors may both play a role in large joint OA but no epidemiological studies have described possible gene environmental interactions. Evidence from twin studies estimate that genetic factors influence 39% to 65% of hand and knee OA in women and about 60% of hip OA [5]. This implies that OA is not only under the influence of genetic factors but also environmental factors and that these may modify the penetrance of OA susceptibility genes, even those showing Mendelian gene inheritance. Cartilage loss and osteophytes are central features in OA and weight-bearing may impact this, therefore, we asked whether there were any interactions in polymorphisms in genes involved in joint tissue integrity (for example, TGFβ1) and environmental factors influencing weight-bearing (for example, weight). We utilised data from the Genetics of OA and Lifestyle (GOAL) cohort, where both body weight [14] and polymorphisms in TGFβ1 [9] have been shown to be independent risk factors, to investigate whether being overweight interacts with TGFβ1 gene polymorphisms to influence the risk of hip or knee OA.

Materials and methods

Participants

The GOAL study was designed as a case-control study to characterise risk factors in large joint OA. The study was approved by the Nottingham Research Ethics Committee and fully informed consent was obtained. Participants were unrelated Caucasian men and women, aged between 45 and 86 years and resident in Nottinghamshire. The recruitment of cases and controls took place between 2002 and 2006 and details of the study are published elsewhere [9,15,16]. In brief, cases were selected if they had clinically significant symptomatic OA of the hip or knee, sufficient to warrant hospital referral. They were recruited if they had either undergone joint replacement, were on the orthopaedic waiting list or had been referred with symptomatic knee OA to the Nottingham knee OA clinic. Controls were recruited from hospital intravenous urography (IVU) waiting lists and frequency-matched to cases by age (± 2 years) and gender. Controls with no evidence of hip OA on review of their IVU radiographs were invited to take part in the study. Only individuals that met the study inclusion criteria were sampled [9,15,16]. The participation rate of eligible cases and controls was 62% and 56%, respectively. This included 1,042 knee OA cases, 1,006 hip OA cases and 1,123 non OA controls.

Data collection

The GOAL study collected data using an interview-administered questionnaire and clinical examination. The questionnaire collected information on socio-demographic factors, employment history, occupational activity, and significant injury, and also contained detailed questions on other risk factors for OA. Weight (kg) and height (cm) were measured by a trained research nurse during the clinical examination. New knee, hand and pelvis radiographs were taken at the clinical examination unless the participant had undergone radiography not more than two years prior to the study or had undergone total joint replacement (TJR). Radiographic assessment and grading for features of OA have been described in detail elsewhere [9,15,16].

Genotyping

Genomic DNA was extracted from whole blood using Gentra PureGene kit (Qiagen). TaqMan® allelic discrimination genotyping method was applied to detect genotype polymorphisms at the TGFβ1 locus as described previously [9].

Exposure variables

Body mass index (BMI) was calculated in kg/m2 [14]. Possible confounding factors included age, sex, bone mineral density (BMD) , nodal OA, significant joint injury and occupational risk factors. Calcaneal bone density (single heel determined by hand dominance) in g/cm2 was measured using DXA (Norland Apollo 501A00Z) and the age-adjusted z-score was used for BMD in the analysis. The presence of interphalangeal nodes was determined during the clinical examination and nodal OA was defined as Heberden's and/or Bouchard's nodes present in at least two x-rays of each hand [16]. Participants self-reported previous significant joint injury to the knee and hip joint. Also included in the interviewer-administered questionnaire were detailed questions about jobs held since leaving school. For each job reported, information was sought on tasks performed on an average working day that involved 12 specified occupational activities such as kneeling and squatting as well as the weekly frequency of lifting different levels of weights. For this analysis, the longest held occupation was used. We truncated occupational exposure for patients who had undergone TJR so that the longest held job prior to TJR was selected. This was also truncated in all controls to the longest job held three years earlier before the recruitment date [14]. Established occupational risk factors for knee and hip OA including kneeling, squatting, carrying out heavy work whilst standing (≥1 hour), lifting 25 kg (10+ times per week), lifting 50 kg or 100 kg (1+ times per week), were then scored for each individual. Subsequently, each score was summed (maximum score 6) to estimate occupational risks.

Statistical analysis

Analyses were conducted separately on knee and hip OA using SPSS version 18. For the present analysis, all knee OA (n = 1,042) and hip OA (n = 1,006) patients, and controls (n = 967) who had no symptoms and no radiographic evidence of hip or knee OA (defined as Kellgren-Lawrence score (KL) ≤2 at the knee and KL ≤3 at the hip joint) were included.

Participants were stratified into four strata according to genotype (G) and BMI (E): G-/E-, G+/E-, G-/E+ and G+/E+. Genotype 11 was defined as negative (coded as 0) and genotype 12 or 22 was defined as positive (coded as 1). BMI <25 kg/m2 was defined as negative (coded as 0) and BMI ≥25 kg/m2 was defined as positive (coded as 1). The ddds ratio (OR) and 95% CI was calculated for G+/E-, G-/E+ and G+/E+ compared with G-/E-. The study was powered to detect an OR for interaction (ORINT) of 5 based on the multiplicative model, that is, OR11 > OR01 × OR10. Unconditional logistic regression was used to adjust for confounding. Two models were developed: one to only adjust for age (continuous variable) and gender (model 1); the other to adjust for all putative risk factors including age, gender, BMD (continuous variable), significant joint injury (yes/no), nodal OA (yes/no) and occupational risk factors (no risk factors, and 1, 2 and 3 or more knee occupational risk factors). All statistical tests were two-sided (Wald's c² statistic), with 5% significance level.

To examine interaction, between gene polymorphisms and BMI, multiplicative and additive models of interaction were used. To test for multiplicative interactions, ORs were first examined across different strata and then an interaction term (genotype*BMI) was included in a logistic regression model. Multivariate models were used to adjust for covariates listed in association analyses and a P-value below 0.05 (likelihood ratio test, LRT) was used to indicate a significant interaction. In addition, we assessed departure from additivity of absolute effects as recommended by Rothman [17]. Andersson [18] has described a detailed calculation method of additive interaction, an SAS programme, including instructions for how to adjust the programme for use in SPSS, and an Excel calculator available from EpiNET [19] which we used to estimate the three indicators of interaction and 95% CI. The method follows the approach proposed by Lundberg et al. [20] to calculate three measures of additive interactions: relative excess risk due to interaction (RERI), attributable proportion (AP) and synergy index (SI), and a method suggested by Hosmer and Lemeshow [21] to estimate 95% CI. We considered RERI and AP to be equal to 0 and SI equal to 1 to indicate absence of additive interactions [18]. Conversely, additive interaction is considered present if RERI and AP do not equal 0 and SI exceeds unity. Furthermore, if RERI is greater than 0, this denotes a synergetic interaction, which implies that the combined action between two exposures in an additive model is greater than the sum of the individual effects. However, if RERI is less than 0, it implies an antagonistic interaction - meaning that in the presence of two exposures in an additive model, the action of one exposure variable diminishes the effect of the other [17].

Results

Table 1 shows the characteristics of the OA cases and controls. Overall, OA cases compared to controls were significantly older, had a higher BMD, and were more likely to be obese, to have nodal OA, to report previous injury and to have exposure to occupational risk factors.

Table 1. Characteristics of the study population

Stratified analysis

Table 2 shows the results of stratified analysis of TGFβ1 polymorphisms, being overweight (BMI ≥25 kg/m2) and the risk of knee OA. For example, compared with the reference group, (that is, normal weight (BMI<25 kg/m2) individuals carrying the wild-type genotype (11) of SNP rs2278422, the risk for knee OA was approximately 7.0-fold higher (OR = 6.95, 95% CI 4.22, 11.47) in overweight individuals only, 1.4-fold (OR = 1.43, 95% CI 0.83, 2.43) in individuals carrying the variant allele (12/22) only and 4.7-fold higher (OR = 4.72. 95% CI 2.92, 7.64) when both risk factors were present, suggesting an antagonistic interaction between the variant allele of this polymorphism and being overweight.

Table 2. Risk estimates for the association of TGFβ1 SNPs with knee OA by being overweight (BMI >25k g/m2)

Similar stratified analysis for hip OA is presented in Table 3. For example, for SNP rs1800468, the risk for hip OA was significantly elevated in overweight individuals only (OR = 1.71, 95% CI 1.33, 2.18) but there was no significant association with genotype 12/22 only (P = 0.071). However, there was a 2-fold (OR = 2.21, 95% CI 1.55, 3.15) risk for hip OA when both risk factors were present, suggesting a synergistic interaction between the variant allele of this polymorphism and being overweight.

Table 3. Risk estimates for the association of TGFβ1 SNPs with hip OA by being overweight (BMI >25 kg/m2)

Interaction term in the models

Table 4 shows a summary of significant results from the analysis of interaction between TGFβ1 polymorphisms and BMI in OA (see Additional file 1 for results on all SNPs). For knee OA, a significant antagonistic interaction on both the multiplicative and additive scale was observed for SNP rs2278422. Assuming a multiplicative scale, the OR for interaction was 0.47 (95% CI 0.26, 0.86) whilst under an additive scale, the AP was -0.61 (95% CI -0.95, -0.27), indicating that both factors acted antagonistically (RERI = -2.66, SI = 0.59) in relation to risk of knee OA, so that the AP to knee OA was 61% lower than expected from the addition of separate effects of genotype 12/22 and being overweight. Additive interaction was also observed between SNP rs11466321 and being overweight but this effect became null after adjusting for potential confounders. Furthermore, no statistically significant multiplicative interaction (ORINT = 1.08, P = 0.851) was found for this SNP (Table 4).

Table 4. Significant results from gene-environment interactions of TGFβ1 SNPs and BMI in OA

Additional file 1. Additive and multiplicative gene-environment interactions in knee and hip osteoarthritis. This file contains two tables with results on additive and multiplicative gene-environment interactions among transforming growth factor (TGF)β1 single nucleotide polymorphisms (SNPs), body mass index (BMI) and osteoarthritis (OA). Table S1: Gene-environment interactions among TGFβ1 SNPs, BMI and knee OA. Table S2: Gene-environment interactions among TGFβ1 SNPs, BMI and hip OA.

Format: DOC Size: 151KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

For hip OA, an additive interaction was observed between SNP rs1800468 and being overweight. Both factors acted synergistically (RERI = 0.93, SI = 4.33) to increase the risk of hip OA, so that among overweight persons carrying the variant allele, 42% (AP = 0.42; 95% CI 0.16, 0.68) of hip OA risk was attributable to the action of both exposures as compared with the contribution of each of the two risk factors added to each other. A significant multiplicative interaction (ORINT = 2.27, P = 0.015) for this SNP was also found (Table 4).

Discussion

This is the first study to examine possible interactions between increased BMI and TGFβ1 gene polymorphisms for knee and hip OA. We found significant additive and multiplicative interaction between being overweight and the variant allele of TGFβ1 SNP rs2278422 in knee OA but in hip OA these interactions were indicated by TGFβ1 SNP rs1800468.

Obesity is a known risk for OA and it has been postulated that in addition to the impact on joint mechanics [22,23] there may be a link between dysfunctional metabolism and joint damage via adipose-secreted cytokines termed adipokines, which are known regulators of metabolic homeostasis [24,25]. In vitro experiments with human and animal culture cells as well as in vivo animal studies have shown an increased expression of TGFβ1 in response to mechanical stimuli in a number of cell and tissue types [6,7,12,13]. Furthermore, experimental models have found excess endogenous TGFβ in association with both chondrocyte synthesis and osteophyte formation [6,7]. Leptin may also play a key role in the development of OA mainly through mechanisms that modulate TGFβ function in maintaining cartilage and bone integrity [24,25]. There is evidence of increased leptin expressions in synovial fluid and in cartilage and osteophytes of patients with OA [24]. Therefore, it is plausible that excessive biomechanical loading of joints or stimulation by adipokines such as leptin may stimulate the up-regulation of endogenous TGFβ1, leading to formation of new fibrocartilage that then undergoes endochondral ossification to become osteophytes. This hypothesis may partly explain the antagonistic interaction that was observed in this study between being overweight and the variant allele of the TGFβ1 polymorphism rs2278422 in knee OA (Table 4). A previously reported association using the GOAL population found the heterozygous genotype of the TGFβ1 polymorphism rs2278422 to be associated with reduced risk of knee OA [9].

On the other hand, we found increased BMI to interact with the variant allele of SNP rs1800468 to increase the risk of hip OA (Table 4). This finding suggests that a greater body weight increased the risk of hip OA in this population, but overweight individuals with the variant allele of TGFβ1 SNP rs1800468 appeared to have an even greater risk (Table 3). Our previous association analysis did not find a statistically significant association between SNP rs1800468 and hip OA, further demonstrating that on its own, it did not influence risk for hip OA in this population [9]. These different findings at the knee and hip further support consideration of OA at these two sites to be discrete subsets in terms of risk factor profile.

The functional relevance of these TGFβ1 variants remains unknown. The location of SNP rs1800468 in the promoter genomic region (5'UTR) suggests that it may play a role in gene regulation and signal sequence. SNP rs2278422 resides in intron 8, a region of unknown biological importance. Nevertheless, we cannot exclude the possibility that the observed interactions may be due to linkage disequilibrium (LD). For example, SNP rs1800468 has been found to be in LD with SNP rs1800469 [26] and SNP rs1800469 is also reported to be in strong LD with rs1982073 [27]. Therefore, it is conceivable that these interactions may result from this SNP being tightly linked to another susceptibility genetic variant of biological importance.

Some of the strengths of this study include: a well-characterised cohort; a sufficient sample size to detect modest interactions in a case-control design; and radiographic data on the control population. In addition, confounding resulting from population stratification was minimised, as all participants were Caucasian and were recruited from a homogenous geographical area.

However, this study also has several caveats. First, we recruited hospital cases and controls which may have led to selection bias. Cases were patients with clinically severe symptoms of OA sufficient to warrant hospital referral, so these findings may not be extrapolated to the full spectrum of OA patients. Second, controls were recruited from IVU lists, and hence are more likely to suffer from kidney disease, including malignancy. Although none of these conditions have any known positive or negative association with OA, controls with these diagnoses may be positively or negatively associated with the TGFβ1 gene and result in exaggerated risks. Nevertheless, this was unlikely to have affected our findings since the TGFβ1 polymorphisms we examined did not deviate from Hardy-Weinberg Equilibrium (HWE) in the control population. Third, the study measured cross-sectional weight and height at the end-stage of the disease. Therefore, it is not possible to determine whether being overweight preceded OA or was the consequence of it. Despite this limitation, BMI provided a more reliable proxy indicator for mechanical stresses compared to physical activity and occupational risk factors, which were more prone to recall bias. Fourth, we carried out a number of multiple statistical tests but no adjustments were made, which may have resulted in false positive significant interactions. A common approach used in many studies to attempt to reduce the number of polymorphisms tested is to restrict gene-environment interaction analysis only to gene variants that show significant association with the outcome. However, this approach gives little consideration to a large environmental risk component such as BMI in OA. Instead, we have described and provided the tests of significance used as recommended by Pernerger [28] as a way of dealing with multiple testing. Finally, our results were not replicated in a separate independent dataset, therefore, future studies are needed to validate these findings in different populations and study settings.

Conclusions

In conclusion, we have identified for the first time potential interaction between TGFβ1 and being overweight in large joint OA. The results were supported by both multiplicative and additive models and appeared to differ between knee and hip OA. The results require replication in other populations and if confirmed, further molecular studies are warranted to better understand the underlying mechanisms responsible for such interactions in OA.

Abbreviations

AP: attributable proportion; BMD: bone mineral density; BMI: body mass index; GOAL: study Genetics of OA and Lifestyle study; HWE; Hardy-Weinberg equilibrium; IVU: intravenous urography; KL: Kellgren-Lawrence; LD: linkage disequilibrium; LRT: likelihood ratio test; OA: osteoarthritis; OR: odds ratio; ORINT: odds ratio for interaction; RERI: relative excess risk due to interaction; SI: synergy index; SNP: single nucleotide polymorphism; TGFβ1: transforming growth factor beta 1; TJR: total joint replacement; UTR: untranslated region.

Competing interests

Rose A Maciewicz owns stock or stock options in AstraZeneca and has submitted patent finding on the OA gene. All other authors declare no competing interests.

Authors' contributions

SGM participated in the design of the study, collected, analyzed and interpreted the data, and drafted the manuscript. SD participated in the design of the study, acquisition of data, and critically edited and revised the manuscript. WZ participated in the design of the study and analysis and interpretation of the data, and critically edited and revised the manuscript. RAM conceived and designed the study, interpreted the data, and critically edited and revised the manuscript. KM and MD conceived and designed the study, participated in the analysis and interpretation of the data, and critically edited and revised the manuscript. All authors read and approved the final manuscript for publication.

Acknowledgements

We are grateful to AstraZeneca, Macclesfield, UK for financial support for the GOAL study, for assisting with genotyping and maintaining the DNA collection, and to Arthritis Research, UK for infrastructure support during sample collection (ICAC grant 14851).

Infrastructure support for Academic Rheumatology, University of Nottingham, was provided by the Arthritis Research, UK (ICAC grant 14851). Support for the initial assembly of the Genetics of Osteoarthritis and Lifestyle (GOAL) Study was provided by AstraZeneca, Macclesfield, UK.

References

  1. Felson DT, Lawrence RC, Dieppe PA, Hirsch R, Helmick CG, Jordan JM, Kington RS, Lane NE, Nevitt MC, Zhang Y, Sowers M, McAlindon T, Spector TD, Poole AR, Yanovski SZ, Ateshian G, Sharma L, Buckwalter JA, Brandt KD, Fries JF.: Osteoarthritis: new insights. Part 1: the disease and its risk factors.

    Ann Intern Med 2000, 133:635-646. OpenURL

  2. Blagojevic M, Jinks C, Jeffery A, Jordan KP: Risk factors for onset of osteoarthritis of the knee in older adults: a systematic review and meta-analysis.

    Osteoarthritis Cartilage 2009, 18:24-33. OpenURL

  3. Lievense AM, Bierma-Zeinstra SMA, Verhagen AP, van Baar ME, Verhaar JAN, Koes BW: Influence of obesity on the development of osteoarthritis of the hip: a systematic review.

    Rheumatology 2002, 41:1155-1162. OpenURL

  4. Ryder JJ, Garrison K, Song F, Hooper L, Skinner J, Loke Y, Loughlin J, Higgins JPT, MacGregor AJ: Genetic associations in peripheral joint osteoarthritis and spinal degenerative disease: a systematic review.

    Ann Rheum Dis 2008, 67:584-591. OpenURL

  5. Spector TD, MacGregor AJ: Risk factors for osteoarthritis: genetics.

    Osteoarthritis Cartilage 2004, 12(Suppl A):S39-44. OpenURL

  6. Blaney Davidson EN, Vitters EL, van der Kraan PM, van den Berg WB: Expression of transforming growth factor-beta (TGFbeta) and the TGFbeta signalling molecule SMAD-2P in spontaneous and instability-induced osteoarthritis: role in cartilage degradation, chondrogenesis and osteophyte formation.

    Ann Rheum Dis 2006, 65:1414-1421. OpenURL

  7. van Beuningen HM, Glansbeek HL, van der Kraan PM, van den Berg WB: Osteoarthritis-like changes in the murine knee joint resulting from intra-articular transforming growth factor-beta injections.

    Osteoarthritis Cartilage 2000, 8:25-33. OpenURL

  8. Pottenger LA, Phillips FM, Draganich LF: The effect of marginal osteophytes on reduction of varus-valgus instability in osteoarthritic knees.

    Arthritis Rheum 1990, 33:853-858. OpenURL

  9. Limer KL, Tosh K, Bujac SR, McConnell R, Doherty S, Nyberg F, Zhang W, Doherty M, Muir KR, Maciewicz RA: Attempt to replicate published genetic associations in a large, well-defined osteoarthritis case-control population (the GOAL study).

    Osteoarthritis Cartilage 2009, 17:782-789. OpenURL

  10. Yamada Y, Okuizumi H, Miyauchi A, Takagi Y, Ikeda K, Harada A: Association of transforming growth factor beta1 genotype with spinal osteophytosis in Japanese women.

    Arthritis Rheum 2000, 43:452-460. OpenURL

  11. Alessi MC, Bastelica D, Morange P, Berthet B, Leduc I, Verdier M, Geel O, Juhan-Vague I: Plasminogen activator inhibitor 1, transforming growth factor-beta1, and BMI are closely associated in human adipose tissue during morbid obesity.

    Diabetes 2000, 49:1374-1380. OpenURL

  12. Klein-Nulend J, Roelofsen J, Sterck JG, Semeins CM, Burger EH: Mechanical loading stimulates the release of transforming growth factor-beta activity by cultured mouse calvariae and periosteal cells.

    J Cell Physiol 1995, 163:115-119. OpenURL

  13. Robbins JR, Evanko SP, Vogel KG: Mechanical loading and TGF-beta regulate proteoglycan synthesis in tendon.

    Arch Biochem Biophys 1997, 342:203-211. OpenURL

  14. Holliday KL, McWilliams DF, Maciewicz RA, Muir KR, Zhang W, Doherty M: Lifetime body mass index, other anthropometric measures of obesity and risk of knee or hip osteoarthritis in the GOAL case-control study.

    Osteoarthritis Cartilage 2011, 19:37-43. OpenURL

  15. Doherty M, Courtney P, Doherty S, Jenkins W, Maciewicz RA, Muir K, Zhang W: Nonspherical femoral head shape (pistol grip deformity), neck shaft angle, and risk of hip osteoarthritis: a case-control study.

    Arthritis Rheum 2008, 58:3172-3182. OpenURL

  16. Zhang W, Robertson J, Doherty S, Liu JJ, Maciewicz RA, Muir KR, Doherty M: Index to ring finger length ratio and the risk of osteoarthritis.

    Arthritis Rheum 2008, 58:137-144. OpenURL

  17. Rothman KJ: Epidemiology: An Introduction. New York: Oxford University Press; 2002. OpenURL

  18. Andersson T, Alfredsson L, Kallberg H, Zdravkovic S, Ahlbom A: Calculating measures of biological interaction.

    Eur J Epidemiol 2005, 20:575-579. OpenURL

  19. Excel sheet to calculate measures of biological interaction [http://www.epinet.se/Epidemiologicaltools.htm] webcite

    [Accessed: 16 April 2013]

  20. Lundberg M, Fredlund P, Hallqvist J, Diderichsen F: A SAS program calculating three measures of interaction with confidence intervals.

    Epidemiology 1996, 7:655-656. OpenURL

  21. Hosmer DW, Lemeshow S: Confidence interval estimation of interaction.

    Epidemiology 1992, 3:452-456. OpenURL

  22. Felson DT: Does excess weight cause osteoarthritis and, if so, why?

    Ann Rheum Dis 1996, 55:668-670. OpenURL

  23. Sharma L, Lou C, Cahue S, Dunlop DD: The mechanism of the effect of obesity in knee osteoarthritis: the mediating role of malalignment.

    Arthritis Rheum 2000, 43:568-575. OpenURL

  24. Dumond H, Presle N, Terlain B, Mainard D, Loeuille D, Netter P, Pottie P: Evidence for a key role of leptin in osteoarthritis.

    Arthritis Rheum 2003, 48:3118-3129. OpenURL

  25. Loeser RF: Systemic and local regulation of articular cartilage metabolism: where does leptin fit in the puzzle?

    Arthritis Rheum 2003, 48:3009-3012. OpenURL

  26. Grainger DJ, Heathcote K, Chiano M, Snieder H, Kemp PR, Metcalfe JC, Carter ND, Spector TD: Genetic control of the circulating concentration of transforming growth factor type beta1.

    Hum Mol Genet 1999, 8:93-97. OpenURL

  27. Langdahl BL, Carstens M, Stenkjaer L, Eriksen EF, Langdahl BL, Carstens M, Stenkjaer L, Eriksen EF: Polymorphisms in the transforming growth factor beta 1 gene and osteoporosis.

    Bone 2003, 32:297-310. OpenURL

  28. Perneger TV: What's wrong with Bonferroni adjustments.

    BMJ 1998, 316:1236-1238. OpenURL