FTO gene polymorphisms and obesity risk: a meta-analysis

1Obesity is becoming a worldwide epidemic in modern society. It is prevalent in individuals of both genders and of all ages, socio-economic strata, and ethnic groups. It is estimated that the total number of overweight adults has reached more than 1.1 billion worldwide, including 312 million obese individuals and that about 10% of children are classified as overweight or obese. The worldwide obesity epidemic is mainly caused by environmental factors, including excessive food intake and a lack of physical activity. Obesity may lead to cardiovascular disease, type 2 diabetes, and several cancers. Overweight and inactivity account for an estimated one-quarter to a one-third of cancers of the breast, colon, endometrium, kidney, and esophagus.

2In 2007, FTO (fat mass and obesity-associated gene) was first discovered in a genome-wide association study (GWAS) for type 2 diabetes, and, almost simultaneously, two other teams independently reported that the FTO gene was associated with obesity (or obesity-related traits) in a GWAS and a genetic association study. The FTO gene, which is located on chromosome 16q12.2 and has nine exons, is strongly conserved across various vertebrate species (for example, fish and chicken) and emerged 450 million years ago. FTO is mainly expressed in the hypothalamus and encodes a 2-oxoglutarate-dependent nucleic acid demethylase. It may play important roles in the management of energy homeostasis, nucleic acid demethylation, and the regulation of body fat masses by lipolysis.

3The FTO gene has recently attracted much attention in obesity research. Previous genetic association-based studies have shown that SNPs in the FTO gene are associated with increased body mass index (BMI), and/or other metabolic-related traits, such as higher fasting insulin, glucose, triglycerides, lower HDL cholesterol, waist circumference, and weight. For example, Scuteri et al. showed that SNP rs9930506 within FTO was associated with BMI, total body weight, and hip circumference; many variants have been reported to be associated with obesity, including rs9939609, rs1421085, rs8050136, rs17817449, and rs1121980. These single nucleotide polymorphisms (SNPs) lie within a 47-kilobase linkage disequilibrium (LD) block encompassing parts of the first two introns and exon 2 of FTO. The transcriptional start site of the retinitis pigmentosa GTPase regulator-interacting protein 1-like (RPGRIP1L) gene is also near the 5' end of FTO. Based on this observation, Fawcett et al. argued that the association signal could be due to a correlation between FTO intronic SNPs and RPGRIP1L.

4Of the SNPs that were reported to associate with obesity, rs9939609 has been of particular interest since it was discovered by Frayling et al. The associations of other SNPs in the FTO gene with obesity (or overweight) have been replicated in large European populations. In some different ethnic groups, however, such as the Chinese and Oceanic populations, no association was observed between rs9939609 and obesity by either genetic association studies or GWAS. However, as more studies are reported, particularly with respect to various ethnicities, there is a need to determine more precisely the effect sizes in each major racial group and to investigate the minor allele frequency (MAF) and the LD patterns of the SNPs across different ethnicities. In addition, some reports have claimed ethnic-specific associations with alternative SNPs, and to that end there has been a degree of confusion. Therefore, we conducted meta-analyses of all available data.

5We searched PubMed, MEDLINE, Web of Science, EMBASE, and BIOSIS Preview up to October 2010 for studies concerning the association between FTO polymorphisms and obesity (or obesity-related traits including body weight, leptin levels, subcutaneous fat, fat mass, hip and waist circumference, lean mass, body height, and BMI). There are no limits on language. Two search themes were combined using the Boolean operator "and". The first theme was ("obesity" (Mesh) OR "overweight" (Mesh) OR "body mass index" (Mesh) OR "obesity, morbid" (Mesh) OR "morbid obesity" OR "morbid obesities" OR "BMI" OR "body weight" OR "leptin levels" OR "subcutaneous fatness" OR "fat mass" OR "hip circumference" OR "waist circumference" OR "lean mass" OR "body height"). The second theme was ("rs9939609" OR "rs8050136" OR "rs1421085" OR "rs17817449" OR "rs1121980" OR "FTO" OR "FTO protein, human" (Substance Name) OR "fat mass and obesity associated" (Substance Name)). Meta-analysis articles were not excluded because several original studies were often combined in meta-analysis articles on genetic association studies.

6Genetic association studies and GWASs with case-control subjects in which the case subjects were obese and the control subjects were healthy were included. At least two studies had to be available for each SNP. Studies were excluded if the control subjects violated the Hardy-Weinberg Equilibrium (HWE).

7Two reviewers (SHP and YMZ) independently searched the databases. The search results were then evaluated by five reviewers (SHP, YMZ, FYX, XBR, and XBL). Disagreements were resolved by discussions among the reviewers. Information on gender, author name, country, ethnicity, year of publication, mean age of examination, mean BMI (calculated as weight in kilograms divided by height in meters squared), and genotypes (for example, TT/TA/AA) were retrieved. MAF and P-values of the HWE were calculated from the above genotype data. Adult individuals with BMI ≥30 kg/m2 were defined as case subjects of obesity, and individuals with BMI <30 kg/m2 were considered control subjects. The criteria for childhood obesity described by Cole et al. were adopted.

8To maximize the number of studies included, analyses were conducted by calculating the ORs under per-allele comparison. We found that an additive genetic model was always employed in most of the genetic association studies concerning the association between FTO and obesity risk. Therefore, the overall effects were also calculated in an additive genetic model in this meta-analysis. For comparison, the ORs in dominant and recessive genetic models were also calculated. Recently, Zintzaras reported a novel method to calculate the generalized odds ratio (ORG), which may overcome the shortcomings of multiple model testing or erroneous model specification. Thus, the ORG calculations were also performed.

9All statistical tests were performed with Stata 10.0 software (Stata, College Station, TX, USA). The effect sizes were calculated by the DerSimonian-Laird method, which is a random effect model because a modest or higher heterogeneous status was found in both allelic comparison and additive genetic models.

10The Cochran's χ2 test (Q test) was used to evaluate heterogeneity between studies, and a threshold P-value of 0.1 was considered statistically significant. The I2 test was also conducted to evaluate heterogeneity between studies. Begg's and Egger's tests were both used to test the significance of the publication bias, with a P-value < 0.1 considered as significant. We used the software developed by Guo et al. to test (using exact test) the significance of the HWE in control subjects, with a threshold P-value of 0.05 and 0.001 considered statistically significant for candidate gene association studies and GWAS, respectively. Heterogeneity is considered higher when I2 >50% and much higher when I2 >75%, and higher heterogeneity is a common phenomenon in genetic association studies. To estimate heterogeneity, meta-regression was performed with the mean age of control subjects (Age_Control) and the mean BMI of control subjects (BMI_Control) as the covariates.

11To identify the source of the heterogeneity between studies, we performed sensitivity analyses by including and excluding some studies. Sensitivity analyses were done sequentially for all of the SNPs and all of the studies (or by some subgroups of the studies).

12We sub-grouped the studies into six groups (Caucasian, Asian, Hispanic, African, South American, and Mixed) on rs9939609, four groups (Caucasian, Asian, African, and Mixed) on rs8050136, and two groups (Caucasian, Asian) on rs1121980.

13The LD patterns of the SNPs were investigated using the HapMap Database (http://hapmap.ncbi.nlm.nih.gov/) and HaploView software (Whitehead Institute for Biomedical Research, Cambridge, USA).

14According to the search strategy, 170, 142, 160, 167, and 171 articles (a total of 810 articles) were retrieved from PubMed, MEDLINE, Web of Science, EMBASE, and BIOSIS Preview, respectively. After the first screening, in which the abstracts or titles were read, 753 articles were excluded, and 59 articles underwent further review. After reviewing these articles, 30 additional articles were excluded, which left a total of 27 articles for inclusion in this meta-analysis. According to the PRISMA guidelines, the flow diagram is shown in Figure 1. A total of 21 articles included 30 studies on rs9939609; five articles included eleven studies on rs1421085; seven articles included nine studies on; rs8050136, two articles included six studies on rs17817449; and three articles included three studies on rs1121980. Not counting the overlapping literature, 27 articles were obtained (Figure 1). In the 30 studies concerning rs9939609, subgroup analyses were performed. The 30 studies were divided into six subgroups according to ethnicity, as follows: studies with Caucasian populations, Asian populations, Hispanic groups, South American, African, and studies using mixed ethnic populations. In all of the included studies, the genotype distributions in the control subjects are consistent with the HWE. Three studies by Song et al. (concerning rs9939609 and rs8050136) and one study by Price et al. were excluded due to the deviation from the HWE (with a P-value of less than 0.05).

15Unfortunately, some articles were excluded because the associations were between one or more SNPs and obesity-related traits rather than obesity itself (for example, ref.). There are five included articles with OR values and other detailed information but without genotype data. We attempted to obtain the genotyped data from these studies by e-mail, but we received no response from the authors.

16The clinical characteristics of the included studies are shown in Table 1. The genotype data, as well as the ORs under both per-allele comparisons and the additive genetic model, are shown in Table 2. The ORs under dominant and recessive genetic models and the ORG are shown in Table S1 (see Additional file 1). These studies were published between 2007 and 2010. A total of 59 studies relating to all five SNPs were included, involving a total of 41,734 obesity cases and 69,837 healthy controls.

17We estimated the MAF in the five polymorphisms from the control subjects of all studies identified for inclusion in the present meta-analysis. Across all studies, the MAFs ranged between 11% and 45%, 40% and 46%, 11% and 44%, 36% and 60%, and 21% and 44% for rs9939609, rs1421085, rs8050136, rs17817449, and rs1121980, respectively.

18Under per-allele comparison, the OR is not available from the study by Li et al., leaving 29 studies for further consideration. A total of 21 out of 29 studies reported a significant, positive association between obesity and the rs9939609 genotype (Table 2, Figure 2). A significant association between obesity and rs9939609 was detected, with an overall OR of 1.31 (95% CI: 1.26 to 1.36) under per-allele comparison, and there is evidence of heterogeneity (I2 = 44.0%). The Begg's test (P = 0.39), and Egger's test (P = 0.17) provided no evidence of publication bias.

19When stratifying the data by ethnicity, no evidence of a significant association was observed in the three studies with Hispanic ethnic groups (summary OR: 1.33, 95% CI: 0.84 to 2.10) under per-allele comparison, with the largest heterogeneity in all six ethnic groups (I2 = 85.2%, P = 0.001) (Figure 2). We observed overall ORs of 1.30 (95% CI: 1.24 to 1.37) and 1.35 (95% CI: 1.27 to 1.44) in the Caucasian and Asian ethnic groups for which more studies had been performed under per-allele comparison.

20By using meta-regression, we detected a significant correlation between the mean control BMI and effect size in an allelic comparison (P = 0.03) (see Additional file 1: Figure S1).

21Very similar results were obtained using an additive genetic model, with an overall OR of 1.31 (95% CI: 1.25 to 1.37) and evidence of heterogeneity (I2 = 52.8%) (Table 2, Additional file 1: Figure S2).

22A total of 11 out of 12 studies reported a significant, positive association between obesity and the rs1421085 genotype (Table 2, Figure 3A). Also under per-allele comparison, a significant association between obesity risk and rs1421085 was found (overall OR = 1.43, 95% CI: 1.33 to 1.53), and there is evidence of heterogeneity (I2 = 38.1%). Begg's test (P = 0.76) and Egger's test (P = 0.84) provided no evidence of publication bias.

23By using meta-regression, no significant correlation was found between the mean control BMI and effect size in an allelic comparison (P = 0.50).

24Overall, 11 studies have investigated the association of rs8050136 with obesity. Of the 11 studies, the OR (under per-allele comparison) of the study by Li et al. is not available. Of the remaining 10 studies, 5 reported a significant, positive association between obesity and the rs8050136 genotype (Table 2, Figure 3B). A significant association between obesity risk and rs8050136 was found (overall OR = 1.25, 95% CI: 1.13 to 1.38) under per-allele comparison, and there is evidence of higher heterogeneity (I2 = 72.0%). Begg's test (P = 1.0) and Egger's test (P = 0.87) provided no evidence of publication bias.

25By using meta-regression, no significant correlation was found between the mean control BMI and effect size in an allelic comparison (P = 0.60).

26All six of the studies involving rs17817449 reported a significant, positive association between obesity and rs17817449 (Table 2, Figure 4A). A significant association between obesity risk and rs17817449 was found (overall OR = 1.54, 95% CI: 1.41 to 1.68) for obesity and rs17817449 under per-allele comparison, and there is evidence of higher heterogeneity (I2 = 45.2%). Begg's test (P = 0.85) and Egger's test (P = 0.90) provided no evidence of publication bias.

27By using meta-regression, no significant correlation was found between the mean control BMI and effect size in an allelic comparison (P = 0.57).

28All three studies reported a significant, positive association between obesity and the rs1121980 SNP (Table 2, Figure 4B). A significant association between obesity risk and rs1121980 was found (overall OR = 1.34, 95% CI: 1.10 to 1.62) under per-allele comparison, with evidence of greater heterogeneity (I2 = 84.5%). Begg's test (P = 0.12) provided no evidence of publication bias, whereas the funnel plot (not shown) and Egger's test (P = 0.09) showed evidence of publication bias.

29By using meta-regression, no significant correlation was found between the mean control BMI and effect size in an allelic comparison (P = 0.28).

30To our knowledge, this meta-analysis is the first one to investigate the associations between FTO polymorphisms and obesity risk across different ethnic groups. This meta-analysis investigated the associations between five FTO polymorphisms (rs9939609, rs1421085, rs8050136, rs17817449, and rs1121980) and obesity risk in 41,734 cases and 69,837 controls from 59 studies, counting the cases and control subjects from each study only once. We found significant evidence for a modest increase in the risk of obesity associated with the five polymorphisms in various ethnic populations (Figures 2, 3, 4). However, subgroup analyses showed that in some ethnic populations, for example, rs9939609 in Hispanic and African, rs1121980 in Caucasian, and rs8050136 in Asian and African, significant associations were not found between the SNPs and obesity risk (see Additional file 1: Figure S3). These results may be due to the differences in MAF or LD patterns across different ethnic populations.

31Interestingly, we found that MAF values are very similar across all five SNPs. For example, the MAFs in the study by Price et al. are 0.36 for rs9939609 and 0.36 for rs17817449. There are also ethnic variations of the MAF in the five polymorphisms. Caucasian, US Hispanics, and US Africans have higher MAFs (>31%), whereas Asian populations (for example, Chinese and Japanese) and South Americans have lower MAFs (between 11% and 20%). When stratifying the studies on rs9939609 into subgroups, the MAFs are different in various populations. For example, the MAF in control subjects is 0.34 to 0.44 in Caucasians, 0.11 to 0.20 in Asians, 0.31 to 0.37 in Hispanics, and 0.17 in South Americans, indicating that the MAF in Asian and South American populations is less than half of that in Caucasian and Hispanic populations. The population differences in MAF and LD structure probably have arisen through evolutionary divergence.

32For rs9939609, we conducted sensitivity analyses. After the exclusion of the three children's studies, there is a small increase in overall heterogeneity (I2 = 44.6%, P = 0.008 when exclusion of the three studies vs. I2 = 44.0%, P = 0.006 when inclusion of all available studies) (see Additional file 1: Figure S4). Similarly, a slight increase in overall heterogeneity was observed after the exclusion of the two GWASs (I2 = 51.6%, P = 0.003 with the exclusion of the two studies vs. I2 = 44.0%, P = 0.006 with the inclusion of all available studies). These results indicate that the children's studies or GWASs are not the sources of the heterogeneity. On the other hand, with the exclusion of the study by Price et al., the heterogeneity decreased rapidly (I2 = 0.0%, P = 0.626 with the exclusion of the study vs. I2 = 44.0%, P = 0.006 with the inclusion of all available studies) in the Caucasian subgroup, indicating that the study by Price et al. is the main source of heterogeneity in the Caucasian subgroup. By using an additive genetic model, coupled with either the inclusion or exclusion of the study by Li et al., the I2 value in the heterogeneity test declined sharply (61.7% vs. 15.6%). Similarly, removing the study by Chang et al. caused the I2 value to increase a little (61.7% vs. 67.0%), suggesting that the study by Li et al. is the main source of heterogeneity in the Asian ethnic group.

33All of the studies on rs1421085 are from Caucasian populations, but when the five studies were excluded (two from Switzerland, two from Germany, and one from the USA) and the remaining seven studies were analyzed (six from France and one from Belgium), the I2 decreased sharply from 38.1% to 0.0%. Similarly, the same tendency was found with rs17817449 (data not shown), suggesting that geographic stratification in an ethnic subgroup is one of the sources of the between-study heterogeneity.

34We investigated the LD structure of the five SNPs in all available ethnic populations (Figure 5). We found that there are substantial differences in the LD structures in various ethnic populations. Similar ethnic populations showed similar patterns (Figure 5), including CEU with TSI; YRI, MKK, and LWK with ASW; and CHB and CHD with JPT. All five SNPs are in strong LD with each other in Caucasians (Figure 5). For example, in one Caucasian population, all five SNPs showed strong LD (r2 ≥ 0.83), and in this case, rs9939609 may be selected as a TagSNP for the other four SNPs. Similarly, in three different populations, including Caucasian, East Asian and Gujarati Indians, four of the five SNPs (rs9939609, rs8050136, rs1421085, and rs17817449) are also in strong LD, indicating that rs9939609 can be used as a TagSNP for the other three SNPs.

35Although the five SNPs are in strong LD with each other in Caucasians, there are obvious differences in the effect size across the five SNPs: 1.3 (1.24 to 1.37) for rs9939609, 1.43 (1.33 to 1.53) for rs1421085, 1.41 (1.13 to 1.76) for rs8050136, 1.54 (1.41 to 1.68) for rs17819449, and 1.36 (0.96 to 1.92) for rs1121980 (Figures 2, 3, 4). The first four SNPs reported significant associations, and the last one reported non-significant associations. On the other hand, the five SNPs are not in strong LD across all the ethnic groups. Therefore, no SNP can be treated as tagSNP across all the ethnic groups. Considering these facts, we showed all the results concerning all the five SNPs without selection of tagSNPs. In fact, some studies with very small sample size were included in this meta-analysis, for example, only 240 subjects were included in the study by Song et al. (Song_AP, Table 2). This fact may be one of the reasons for the inconsistent results across the five SNPs in Caucasians. To address this issue, more studies with larger sample sizes are needed to provide more precise evidence.

36The heterogeneities of all five SNPs are higher in the meta-analysis. There are many reasons for heterogeneity in genetic association studies, and these reasons can be divided into two categories: (1) differential biases due to population stratification, misclassification of clinical outcome, differences in BMI, genotyping errors and over-estimation of genetic effects in the first study and (2) differences in the pattern of LD structure across populations. Using sensitivity analyses, we showed that population stratification may be the reason for higher heterogeneity in studies focused on rs1421085, indicating that population stratification can also exist within a subgroup of the same ethnic background (Caucasians). In the meta-regression analysis, a significant association was detected between the mean control BMI and effect size on rs9939609, suggesting that the difference in BMI value between the subjects may also be a source of heterogeneity.

37Analyses of statistical power constitute a crucial step in the design of genetic association studies. Whereas conventional statistical power calculations for case-control studies disregard many basic elements of analytic complexity and can seriously underestimate true sample size requirements. Many methods have been proposed to address the sample size issue in genetic association study. However, owing to many factors, such as financial pressure, difficulties in sample collection, or neglect of the importance of sample size, too few samples were employed in many genetic association studies, for example, Andreasen_c with total subjects of 377 and Song_AP with total subjects of 240.

38Inherent limitations of this meta-analysis should be pointed out, including the following: (1) Based on the statistical tests, the heterogeneities are higher in the results of this meta-analysis. There may be many potential reasons for the higher heterogeneities, including differences in genetic susceptibility across ethnic groups and measurement error of BMI. (2) Gene-gene and gene-environment interactions could not be addressed in this meta-analysis. Actually, obesity is a complex trait, and many genes are related to obesity, such as MC4R, MAF, and NPC1. On the other hand, lifestyle factors, including diet and physical inactivity, are important contributors to weight gain and obesity (3) It may be unreliable to use Begg's and Egger's test as the criteria for publication bias. Although publication bias was not found, it may exist due to the fact that some studies with non-significant associations were not published, and some articles published in other languages (not English) were not obtained and included in this meta-analysis. (4) The Haplotype of the polymorphisms was not performed due to the lack of the related information.

39This meta-analysis suggests that FTO may represent a low-penetrance susceptible gene for obesity risk. Individual studies with large sample sizes are needed to further evaluate the associations between the polymorphisms and obesity risk in various ethnic populations.

40BMI: body mass index; CI: confidence interval; GWAS: genome-wide association study; HWE: Hardy-Weinberg Equilibrium; LD: linkage disequilibrium; MAF: minor allele frequency; OR: odds ratio; ORG: generalized odds ratio; SNP: single nucleotide polymorphism.

41The authors declare that they have no competing interests.

42SHP, YMZ and MDL contributed to the design of the study. SHP, XBR, and XBL carried out the screening procedure. SHP, YMZ and FYX performed the statistical analysis. SHP drafted the manuscript. All authors read and approved the final manuscript.

43The pre-publication history for this paper can be accessed here:

44http://www.biomedcentral.com/1741-7015/9/71/prepub