- Regular article
- Open Access
FTO variant is not associated with osteoarthritis in the Chinese Han population: replication study for a genome-wide association study identified risk loci
Journal of Orthopaedic Surgery and Researchvolume 13, Article number: 65 (2018)
Osteoarthritis is the most prevalent form of arthritis worldwide and is the major cause of pain and loss of function in elderly people. A signal of the fat mass and obesity-associated (FTO) gene had been reported in a genome-wide association study of osteoarthritis. The FTO polymorphism (rs8044769) might exert its effect on osteoarthritis through obesity, because it was reported as a body mass index-associated single-nucleotide polymorphism. And replication studies showed inconsistent results for this association. Our present study is to check the association of rs8044769 with osteoarthritis and body mass index in Chinese Han population.
A case-control association study was conducted by using 890 osteoarthritis cases and 844 controls in Chinese Han population. rs8044769 was genotyped in all subjects. Allelic and genotypic frequencies were compared between osteoarthritis cases and control subjects. Associations between rs8044769 and body mass index, and body mass index and osteoarthritis were also assessed.
No significant difference was detected in genotype or allele distribution between osteoarthritis cases and controls (P > 0.05). Stratification by gender and body mass index revealed negative association between rs8044769 and osteoarthritis. We did not find any solid association between rs8044769 and higher body mass index. Meanwhile, we demonstrated that higher body mass index (body mass index ≥ 25) was associated with osteoarthritis.
Our present study suggested that rs8044769 was not associated with osteoarthritis susceptibility or higher body mass index, and higher body mass index was a risk factor for osteoarthritis in the Chinese Han population. We also proposed that stratification by clinical parameters was crucial to reduce false-positive result in OA association studies.
Osteoarthritis (OA; OMIM 165720) is the most common degenerative disease of the joints among the elderly. Hip and knee OA are the main causes of pain, functional loss, and depressed quality of life in both developed and developing countries. OA, like other complex chronic disease, develops through a combination of environmental and genetic risk factors .
Obesity has been consistently associated with an augmented risk of OA, but its role in the progression of OA has not been well established. Both biomechanical and metabolic factors are likely to link obesity and OA. Biomechanical overload on joints can accelerate cartilage degeneration, while metabolic dysfunction caused by obesity may have a more severe influence on OA than body weight . Recent studies showed adipose tissue could increase cartilage degeneration through a secretion of pro-inflammatory mediators and synergistic regulation with other inflammatory cytokines .
The fat mass and obesity-associated (FTO) gene encodes a 2-oxoglutarate-dependent nucleic acid demethylase. Its structure has been reported, but the exact function of this gene remains unknown [4, 5]. Genetic studies have showed variants in FTO are associated with an increased body mass index (BMI) which is an obesity-related measurement in several studies [6,7,8,9,10,11,12]. A genome-wide association study (GWAS) named arcOGEN study had found that a FTO single-nucleotide polymorphism (SNP) rs8044769 was in the strongest association with female OA close to genome-wide significance (P = 6.85 × 10− 8). This significance was degenerative after BMI adjustment, hinting that this FTO polymorphism may exert its influence on OA through obesity [13, 14]. The association was replicated in a study on North American Caucasians . A study on Chinese Han population denied this association, but the subject number of the Chinese study was a little bit small . Until now, the association between FTO and OA remains debatable.
Therefore, we evaluated the potential association of the FTO polymorphism (rs8044769) with OA in the Chinese Han population with more cases and controls. We also evaluated the association of this polymorphism with BMI because of the known effect of FTO on BMI.
All patients were recruited from the department of Sports Medicine and Adult Reconstructive Surgery, Drum Tower Hospital, School of Medicine, Nanjing University and were diagnosed as having knee OA according to the previous study . All the radio-graphics were estimated by the Kellgren/Lawrence grading system . Other knee diseases like rheumatoid arthritis, post-traumatic or post-septic arthritis, and developmental dysplasia of the knee were excluded. We set the youngest age for the controls which is the same with the cases. The control subjects were consecutively enrolled at the Physical Examination Center, Drum Tower Hospital Affiliated to the Medical school of Nanjing University. They had no sign or symptom of arthritis or joint diseases. All the study populations were non-institutionalized Han Chinese living in or around Nanjing. This study was approved by the ethics review committee of the Nanjing University and conducted according to the Declaration of Helsinki principles. All participants were informed with letter of consent.
Height and weight were measured (without shoes and in light clothes) by medical technicians. BMI was calculated by taking a person’s weight and dividing it by their height squared. Subjects were classified as underweight (BMI < 18.5 kg/m2), normal (18.5 kg/m2 ≤ BMI < 25 kg/m2), overweight (25 kg/m2 ≤ BMI < 30 kg/m2), and obese (BMI ≥ 30 kg/m2) set by the criteria of World Health Organization .
DNA was extracted from peripheral blood leukocytes using the NucleoSpin Blood QuickPure Kit (MACHEREY-NAGEL GmbH & Co. KG, Düren, German). The SNP was genotyped using the Taqman assay in an ABI 7500 system (Applied Biosystems, Foster City, CA, USA). Genotyping was done by laboratory personnel blinded to case status. Three percent of samples were duplicated, and the concordance rate was 100%.
Characteristics of the OA and the control groups were evaluated by using chi-square and Mann-Whitney test. The Hardy–Weinberg equilibrium (HWE) of genotypes in OA and control groups was tested by using chi-square test. The differences of genotype and allele distributions were evaluated by using Chi-square test. The combined effect for the risk allele was evaluated by Mantel-Haenszel test. BMI among different genotypes of FTO polymorphism were compared using Kruskal-Wallis test. The Potential associations between BMI and OA were evaluated by using unconditional logistic regression modeling. All P value reported as statistical significance was set at < 0.05. Odds ratio (OR) was calculated with 95% confidence interval (95% CI). Calculations were performed using SPSS 12.5S (SPSS, Chicago, USA) or Stata 8.2 (Stata, College Station, TX, USA).
The basic characteristics of the study population were showed in Table 1. On average, OA cases were older and had greater BMI than controls. Females were prevalent in the OA group. The mean age was 62.51 ± 11.43 years and mean BMI was 25.76 ± 3.69 kg/m2 in cases compared to 54.07 ± 11.60 years and 24.91 ± 3.04 kg/m2 in controls, respectively. Genotype and allele frequencies for FTO polymorphism were not significantly different between OA cases and controls (Table 2). We combined our data with the data from the previous Chinese study and found no association between the risk allele with OA (Table 2) . We next examined potential association of rs8044769 with OA when stratified by gender and BMI. No significant difference was observed in any comparisons between OA cases and controls (Table 3).
Associations between rs8044769 and BMI were also evaluated. Genotype and allele frequencies for FTO polymorphism were not significantly different between overweight and obese cases and controls (Table 4). No significant associations were observed between overweight and obese cases and controls in females, males, OA cases and OA controls, respectively. We also compared BMI among different genotypes of rs8044769 in Han Chinese population. We found no significant difference of BMI among three genotypes of rs8044769. We then examined associations between rs8044769 and BMI in females, males, OA cases, and controls, respectively, and found no significant differences of BMI among three genotypes of rs8044769 (Table 5).
We last examined potential associations between BMI and OA (Table 6). Through the logistic regression modeling, we demonstrated a higher BMI (overweight and obese, BMI ≥ 25) was a risk factor for OA in the Chinese Han population (Additional files 1 and 2).
In our study population, the mean age of OA patients is significantly higher than that of controls, and the females are significantly more prevalent in OA cases. The mean BMI in OA patients was significantly higher than that in controls, and the logistic regression analysis also showed that BMI was associated with OA. Our data conformed to the previous reports that age, female gender, and overload were risk factors for OA [20, 21].
The association of this FTO polymorphism (rs8044769) with OA was not replicated in the Chinese Han population. The previous GWAS found rs8044769 in FTO was one of the signals close to genome-wide significance for OA in a European population. Our results did not show any positive associations between rs8044769 and OA, even after stratification by gender and BMI. We combined our data with the data from another Chinese study and found no associations between the risk allele of rs8044769 and OA either.
The association of the FTO polymorphism (rs8044769) with BMI was not replicated in the Chinese Han population either. The previous studies showed that FTO was associated with BMI [6,7,8,9,10,11,12], and two of them referred to rs8044769 as BMI-associated SNP [11, 12]. In the association analysis of rs8044769 with BMI, we found no association between rs8044769 and BMI. And the BMIs corresponding different genotypes of rs8044769 showed no significant difference.
There are several possible reasons for the negative association between rs8044769 and OA. The sample size of our study might be insufficient. But regarding the odds ratio of the risk allele in our data and the previous Chinese study, it is difficult to get a significant difference by enlarging sample size. The gender distribution was not matched in OA patients and controls. In our study the genotype distributions did not vary a lot between females and males, and no significant association was found after stratification by gender.
We think the ethnic difference might be an important reason for the negative result in both the association study between rs8044769 and OA and the association study between rs8044769 and BMI, as the allele frequency and genotype distribution rs8044769 were different between Chinese Han population and Caucasian population [13,14,15,16]. The negative association in our population cannot deny the association of rs8044769 with OA and BMI in other populations.
An extended study of the GWAS study revealed that the association between rs8044769 and OA was degenerated after stratification by BMI, and the association between rs8044769 and BMI was identified in this study . There have been several reports for the association between FTO and BMI [6,7,8,9,10,11,12], and two of them referred to rs8044769 as BMI associated SNP [11, 12]. So, we suspected that the association between rs8044769 and BMI may contribute to the association between this SNP and OA in the GWAS study. However, in our study, rs8044769 was not associated with BMI, and it might be another reason for the negative association between rs8044769 and OA.
From our data we cannot get a conclusion that FTO is not associated with BMI in Chinese population as the allele and genotype differences of rs8044769 do not equal to changes of FTO function or FTO expression. Some other genetic variants in FTO might be associated with BMI or OA in Chinese population, and further studies are needed to evaluate the association of FTO with BMI and OA.
The associations among rs8044769, BMI, and OA drew attentions to some concern on current GWAS results. Shortcomings like the imprecise phenotype definitions, chosen criterion, and limited study populations in the GWAS may lead to some false-positive association [22,23,24,25]. Stratification by clinical parameters such as BMI can reduce the risk of false-positive result, but the cost-effectiveness is uneconomical in the primary study. Future studies in OA genetics are required using more accurate phenotype definition, improved case/control selecting standardization, and increased sample sizes as well as extending studies to different populations. In addition, functional characterization and studies of the risk loci can help us understand better of the OA genetics.
Our results were consistent with a previous study in Chinese population and had enlarged studied subjects. We also conducted a meta-analysis together with the previous study in Chinese population to offer a solid evidence that rs8044769 of FTO was not associated with OA risk or BMI in Chinese population. By comparing with the result of the extended study of the GWAS study, we proposed that stratification by clinical parameters was crucial to reduce false-positive result in OA association studies, especially in GWAS studies with large-scale subjects.
Some limitations should be noted in our study. First, the number of study subjects of our study was still limited, although it was much larger than that in the previous Chinese study . The sample size had a statistical power of 0.27 to detect an effect size of 1.10 with an α-level of 0.05 for the association of the SNP rs8044769 with OA risk. Secondly, we did not find detailed allele or genotype data for rs8044769 in the Caucasian reports, so we made a meta-analysis of only Chinese data. Considering the ethnic difference, our data cannot be generalized as a global result. Further studies in Chinese and other populations are needed to evaluate the association of this SNP with OA and BMI.
The present study showed that the FTO polymorphism (rs8044769) was not associated with OA susceptibility in the Chinese Han population. This polymorphism was not associated with BMI in the Chinese Han population either. Higher BMI including overweight and obese is a risk factor for OA in the Chinese Han population. Further studies are necessary to identify the association of FTO with OA in different populations and to proclaim better of the OA pathogenesis.
- 95% CI:
95% confidence interval
Body mass index
- FTO :
Fat mass and obesity-associated
Genome-wide association study
Valdes AM, Spector TD. Genetic epidemiology of hip and knee osteoarthritis. Nat Rev Rheumatol. 2011;7:23–32.
Killock D. Osteoarthritis: the influence of obesity on OA-does size matter or is metabolic dysfunction more important? Nat Rev Rheumatol. 2012;8:61.
Issa RI, Griffin TM. Pathobiology of obesity and osteoarthritis: integrating biomechanics and inflammation. Pathobiol Aging Age Relat Dis. 2012;2:17470.
Gerken T, Girard CA, Tung YC, Webby CJ, Saudek V, Hewitson KS, et al. The obesity-associated FTO gene encodes a 2-oxoglutarate-dependent nucleic acid demethylase. Science. 2007;318:1469–72.
Han Z, Niu T, Chang J, Lei X, Zhao M, Wang Q, et al. Crystal structure of the FTO protein reveals basis for its substrate specificity. Nature. 2010;464:1205–9.
Frayling TM, Timpson NJ, Weedon MN, Zeggini E, Freathy RM, Lindgren CM, et al. A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science. 2007;316:889–94.
Scott LJ, Mohlke KL, Bonnycastle LL, Willer CJ, Li Y, Duren WL, et al. A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science. 2007;316:1341–5.
Song Y, You NC, Hsu YH, Howard BV, Langer RD, Manson JE, et al. FTO polymorphisms are associated with obesity but not diabetes risk in postmenopausal women. Obesity (Silver Spring). 2008;16:2472–80.
Zeggini E, Weedon MN, Lindgren CM, Frayling TM, Elliott KS, Lango H, et al. Replication of genome-wide association signals in UK samples reveals risk loci for type 2 diabetes. Science. 2007;316:1336–41.
Yang J, Loos RJ, Powell JE, Medland SE, Speliotes EK, Chasman DI, et al. FTO genotype is associated with phenotypic variability of body mass index. Nature. 2012;490:267–72.
Zhao J, Bradfield JP, Li M, Wang K, Zhang H, Kim CE, et al. The role of obesity-associated loci identified in genome-wide association studies in the determination of pediatric BMI. Obesity (Silver Spring). 2009;17:2254–7.
Wing MR, Ziegler JM, Langefeld CD, Roh BH, Palmer ND, Mayer-Davis EJ, et al. Analysis of FTO gene variants with obesity and glucose homeostasis measures in the multiethnic insulin resistance atherosclerosis study cohort. Int J Obes. 2011;35:1173–82.
Zeggini E, Panoutsopoulou K, Southam L, Day-Williams AG, Lopes MC, Boraska V, et al. Identification of new susceptibility loci for osteoarthritis (arcOGEN): a genome-wide association study. Lancet. 2012;380:815–23.
Panoutsopoulou K, Metrustry S, Doherty SA, Laslett LL, Maciewicz RA, Hart DJ, et al. The effect of FTO variation on increased osteoarthritis risk is mediated through body mass index: a mendelian randomisation study. Ann Rheum Dis. 2014;73:2082–6.
Yau MS, Yerges-Armstrong LM, Liu Y, Lewis CE, Duggan DJ, Renner JB, et al. Genome-wide association study of radiographic knee osteoarthritis in North American Caucasians. Arthritis Rheumatol. 2017;69:343–51.
Wang Y, Chu M, Rong J, Xing B, Zhu L, Zhao Y, et al. No association of the single nucleotide polymorphism rs8044769 in the fat mass and obesity-associated gene with knee osteoarthritis risk and body mass index: a population-based study in China. Bone Joint Res. 2016;5:169–74.
Jiang Q, Shi D, Yi L, Ikegawa S, Wang Y, Nakamura T, et al. Replication of the association of the aspartic acid repeat polymorphism in the asporin gene with knee-osteoarthritis susceptibility in Han Chinese. J Hum Genet. 2006;51:1068–72.
Kellgren JH, Lawrence JS. Radiological assessment of osteoarthritis. Ann Rheum Dis. 1963;22:237–55.
The World Health Organization. Physical status: the use and interpretation of anthropometry. Report of a WHO expert committee. World Health Organ Tech Rep Ser. 1995;854:1–452.
Lespasio MJ, Piuzzi NS, Husni ME, Muschler GF, Guarino A, Mont MA. Knee osteoarthritis: a primer. Perm J. 2017;21:16–183.
Du H, Chen SL, Bao CD, Wang XD, Lu Y, Gu YY, et al. Prevalence and risk factors of knee osteoarthritis in Huang-Pu District, Shanghai. China Rheumatol Int. 2005;25:585–90.
Kerkhof HJ, Meulenbelt I, Akune T, Arden NK, Aromaa A, Bierma-Zeinstra SM, et al. Recommendations for standardization and phenotype definitions in genetic studies of osteoarthritis: the TREAT-OA consortium. Osteoarthr Cartil. 2011;19:254–64.
Loughlin J. Osteoarthritis: all types of trouble—defining OA in the genomicera. Nat Rev Rheumatol. 2011;7:200–1.
McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ioannidis JP, et al. Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008;9:356–69.
Todd JA. Statistical false-positive or true disease pathway? Nat Genet. 2006;38:731–3.
We thank Qiting Ge (Laboratory for Bone and Joint Diseases, Australia-China Joint Centre, Model Animal Research Center, Nanjing University, Nanjing, Jiangsu, 210009, China.) for the sample preparation and all the subjects who were enrolled in this study for the donation of their sample and clinical information.
This work was supported by the National Natural Science Foundation of China (81420108021).
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This study was approved by the ethics review committee of the Nanjing University and conducted according to the Declaration of Helsinki principles. All participants were informed with letter of consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Genotype and allele frequencies of FTO polymorphism (rs8044769) for association analysis on OA in Han Chinese population when stratified by gender and BMI. (DOC 47 kb)
Table S2. Genotype and allele frequencies of FTO polymorphism (rs8044769) for association analysis on BMI in Han Chinese population when stratified by gender and OA status. (DOC 44 kb)