- Open Access
Waist-hip ratio related genetic loci are associated with risk of impaired fasting glucose in Chinese children: a case control study
Nutrition & Metabolism volume 15, Article number: 34 (2018)
The meta-analyses of genome-wide association studies identified several waist-hip ratio (WHR) related loci in individuals of European ancestry. Since the pattern of fat distribution and the relationship between fat distribution and glucose metabolism disturbance in Chinese are different from those in Europeans, the present study aimed to explore the individual and cumulative effects of WHR-related loci on glycemic phenotypes in Chinese children.
A total of 2030 children were recruited from two independent studies. Eleven single nucleotide polymorphisms (SNPs) were selected and genotyped using matrix-assisted laser desorption ionization time of flight mass spectrometry (MALDI-TOF MS). Logistic regression and linear regression model were used to examine the association of 11 SNPs and genetic risk score (GRS) with impaired fasting glucose (IFG) and fasting plasma glucose (FPG), respectively.
Three SNPs (rs6795735, rs984222 and rs1011731) were nominally associated with IFG (all P < 0.05). Each WHR-increasing (C) allele of rs6795735 (ADAMTS9) was associated with a 40.1% increased risk of IFG (OR = 1.401, 95% CI = 1.131–1.735, P = 0.002), which remained significant after Bonferroni correction. We observed no association of both weighted and unweighted GRS with FPG and IFG (all P > 0.05).
We identified individual effects of rs6795735 (ADAMTS9), rs984222 (TBX15-WARS2), and rs1011731 (DNM3-PIGC) on glycemic phenotypes in Chinese children for the first time. The study suggests that genetic predisposition to central obesity is associated with impaired fasting glucose, providing more evidence for the pathogenesis of diabetes.
Obesity is associated with type 2 diabetes mellitus (T2D) [1, 2]. The prevalence of obesity and T2D has increased during the past few decades, particularly in China [3, 4]. As the metabolic consequences of increasing fat mass are disproportionately attributable to the extent of central adiposity, measures of central adiposity, such as waist-hip ratio (WHR) and waist circumference, predict the risk of serious adverse metabolic outcomes better than body mass index (BMI) . Among the indices involved in central obesity, WHR is of particular interest, as a measure of body fat distribution, because it integrates the adverse metabolic risk associated with increasing waist circumference with the more protective role of gluteal fat deposition .
A large meta-analysis of genome-wide association studies (GWAS) reported 14 single nucleotide polymorphisms (SNPs) associated with WHR independent of BMI . By analyzing the associations between these loci and metabolic traits, they identified one locus (at ADAMTS9) reaching nominal significance for fasting plasma glucose (FPG) and three loci (at ADAMTS9, NISCH-STAB1 and ITPR2-SSPN) reaching nominal significance for T2D (all P < 0.05). Afterwards, some studies also examined the associations of 14 WHR-related loci with glucose-related metabolic traits in other populations. For example, Burgdorf et al.  found that WHR-increasing allele (G) of rs4846567 (LYPLAL1/SLC30A10) was nominally associated with decreased concentration of FPG (P = 0.036). Similarly, Van Vliet-Ostaptchouk et al.  observed that the WHR-increasing allele at CPEB4 (rs6861681) was significantly associated with lower FPG after adjusting for WHR (P = 2.3 × 10− 4), but the genetic risk score (GRS) for WHR was not associated with glucose-related traits. Another study by Huang et al.  showed that the genetic predisposition to central obesity, estimated using GRS based on the 14 WHR loci, was linearly related to higher T2D risk (OR = 1.04, 95% CI = 1.01–1.07, P = 0.01). All of these studies were replicated among European adult populations. However, the effects of the 14 WHR-related loci on glycemic phenotypes have not been analyzed in non-European ethnic groups. What’s more, the effects in children, which might be different from adults, have not been evaluated.
Therefore, we conducted the present study in a case-control sample of Chinese children to explore the individual and cumulative effects of WHR-related loci on glycemic phenotypes. Since the pattern of fat distribution and the relationship between fat distribution and metabolic disturbance in Chinese might be different from those in Europeans [10,11,12], our study will thus provide more evidence for the effects of these loci on glycemic phenotypes in a population with different ethnicity and at a younger age.
We conducted a case-control study in a total of 2030 children from two independent study groups, including 705 obese cases (BMI ≥ the 95th BMI percentile for Chinese children) and 1325 non-obese controls (BMI ≤ the 95th BMI percentile) recruited from the urban regions of Beijing, China. The first group came from the study on adolescent lipids, insulin resistance and candidate genes (ALIR) [13,14,15]. The second study group was from the Comprehensive Prevention project for Overweight and Obese Adolescents (CPOOA) with physical exercise and healthy nutrition as instruments [14,15,16]. The study design, recruitment of participants, and data collection of the two study groups was very similar, except for the age differences (ALIR study: children aged 14–17 years old; CPOOA study: children aged 7–18 years old), which had been described in detail previously [14, 15].
Both the ALIR study and CPOOA study were approved by the Ethics Committee of Peking University Health Science Center. Written informed consent was provided by all participants and, in the case of minors, their parents. Studies were performed according to the Declaration of Helsinki.
Anthropometric and glycemic measurements
Anthropometric measurements, including height, weight, waist and hip circumferences were conducted according to standard protocols [14, 15]. BMI was calculated as weight (kg)/height (m2). For children aged 7–18 years, we used uniform BMI percentile criteria, which were determined in a representative Chinese population , to screen the obese and non-obese participants. Children with an age- and sex-specific BMI ≥ 95th percentile were defined as obese; those with an age- and sex-specific BMI less than the 95th percentile were defined as non-obese. Individuals who were underweight, according to the Chinese national screening criteria for malnutrition of school-age children and adolescents , and those with any vital organs diseases, such as heart disease, liver diseases, pulmonary diseases and kidney disease, were excluded.
Blood samples were drawn after a 12 h overnight fast. FPG concentrations were analyzed by a biochemical auto-analyzer (Hitachi 7060, Tokyo, Japan). Individuals with a FPG ≥ 5.6 mmol/L were considered as impaired fasting glucose (IFG) .
SNPs selection and genotyping
The GWAS meta-analysis reported 14 SNPs associated with WHR . Power calculation was performed using Quanto software (University of Southern California, Los Angeles, CA). With the assumed effect size (OR = 0.14) and effect allele frequency ≥ 0.12, statistical power to detect a positive association would be ≥75%, given our sample size. Hence, 12 SNPs with the minor allele frequency ≥ 0.12 in the Chinese population (Hapmap database: http://hapmap.ncbi.nlm.nih.gov) were initially selected and genotyped for the present study. Since Hardy-Weinberg equilibrium (HWE) was violated (P < 0.001) for rs718314, it was excluded from subsequent data analyses, leaving 11 SNPs for our present study. Genomic DNAs of all children were extracted from blood leukocytes by the phenol-chloroform extraction. Genotyping of the 11 SNPs was conducted with a MassARRAY System (Sequenom, San Diego, CA). The genotyping process has been described in detail in our previous studies [20, 21]. The call rates of the 11 SNPs ranged from 98.9% to 99.9%. We also performed genotyping on randomly selected duplicated samples, and the concordance rates were 100%.
The genotype data of the non-obese and obese groups were tested for deviation from HWE separately. In order to estimate the effect of population subdivision, F-statistics (FST) was calculated [22, 23]. A F ST value between 0 and 0.05 suggests little genetic differentiation; a value between 0.05 and 0.15, moderate differentiation; a value between 0.15 and 0.25, large differentiation; and values above 0.25, very large differentiation .
A logistic regression model was used to examine the independent effect of each SNP on risk of IFG under an additive model adjusted for study group, sex, age and age squared, without or with adjustment for age- and sex-specific WHR-Z scores. WHR-Z scores were calculated using the means and standard deviations of our study population. A general linear model was used to examine the independent effect of each SNP on the FPG under an additive model with adjustment for study group, sex, age and age squared, without or with adjustment for age- and sex-specific WHR-Z scores.
To explore cumulative effects of the 11 SNPs on risk of IFG and FPG, weighted and unweighted GRS was implemented. None of the 11 SNPs was found to be in linkage disequilibrium (LD) with each other at an r2 < 0.05, which was calculated in Han Chinese population with LDLink (https://analysistools.nci.nih.gov/LDlink/). Unweighted GRS was calculated by adding the number of effect alleles of 11 SNPs carried by each individual. Weighted GRS was calculated using reported effect sizes (β-coefficients) in the original study, with the following equation: weighted GRS = (β1× SNP1 + β2× SNP2 +. .. + βn × SNPn) × (total number of SNPs / sum of the β-coefficients). The genotypes of each SNP were coded as 0, 1 and 2 according to the number of effect alleles. Effect alleles in the current study were consistent with those in the original study. Logistic regression and linear regression were used to study the associations between GRS and risk of IFG and FPG, respectively.
Due to the reported sexual dimorphism of WHR , we also performed sex-specific analysis for the single SNP analyses and GRS association. A two-sided P value < 0.05 was considered as nominal significance. Adjustment was made for multiple testing using Bonferroni correction for 11 SNPs and 2 models (with and without adjusting for the WHR z-scores), i.e. resulting in 0.00227 (0.05 divided by 22) as statistical significance (two-sided). Statistical analyses were performed using SPSS 18.0 (SPSS Inc., Chicago, IL).
The general characteristics of the study populations are shown in Table 1. The difference in age and FPG between obese and non-obese groups was not significant (P > 0.05), while the differences in gender, BMI, and waist circumference between the two groups were significant (P < 0.05). The comparisons of individuals with and without IFG were shown in the Additional file 1: Table S1.
Table 2 shows the genotype information of the 11 SNPs in 2030 Chinese children. HWE was fulfilled for all SNPs both among non-obese children and obese children (P > 0.05). FST values between the population of original study and our sample are also shown in Table 2. We identified 5 SNPs (rs1011731, rs10195252, rs6795735, rs1294421, rs1055144 in or near DNM3-PIGC, GRB14, ADAMTS9, LY86, NFE2L3, respectively) with moderate genetic differences between the original and our study groups (0.05 ≤ FST ≤ 0.15), while the remaining 6 SNPs had similar allele frequencies compared to the previous data  (FST ≤ 0.05).
The associations between the SNPs and risk of IFG in 2030 Chinese children are shown in Table 3. Three SNPs were found to be associated with the risk of IFG independent of WHR. Each WHR-increasing (G) allele of rs984222 (TBX15-WARS2) was associated with a 24.9% reduced risk of IFG (OR = 0.751, 95% CI = 0.620–0.910, P = 0.004); each WHR-increasing (G) allele of rs1011731 (DNM3-PIGC) was associated with a 38.3% increased risk of IFG (OR = 1.383, 95% CI = 1.063–1.799, P = 0.016); each WHR-increasing (C) allele of rs6795735 (ADAMTS9) was associated with a 40.1% increased risk of IFG (OR = 1.401, 95% CI = 1.131–1.735, P = 0.002), with additional adjustment for age- and sex-specific WHR-Z scores. The associations of rs6795735 (ADAMTS9) with the risk of IFG remained significant after Bonferroni correction. However, we observed no association of both weighted and unweighted GRS with IFG (all P > 0.05).
The associations between the SNPs and FPG in 2030 children are shown in Table 4. A nominally significant association was observed between rs984222 (TBX15-WARS2) and FPG. After additional adjustment for age- and sex-specific WHR-Z scores, each WHR-increasing (G) allele of rs984222 was associated with 0.032 mmol/L decrease in FPG (β = − 0.032, S.E. = 0.016, nominal P = 0.045). Both weighted and unweighted GRS was not associated with FPG concentrations, adjusted for study group, sex, age, age squared, with and without age- and sex-specific WHR-Z scores (all P > 0.05).
Except a nominally significant association between rs10195252 (GRB14) and WHR (P = 0.04), we did not observe associations for the remaining 10 SNPs and WHR in 2030 children (all P > 0.05; Additional file 1: Table S2). In addition, significant association of these SNPs with BMI was also not detected (all P > 0.10; Additional file 1: Table S2).
Further stratified analysis for sex showed that rs984222 (TBX15-WARS2) was nominally associated with IFG only in boys, while rs6795735 (ADAMTS9) and rs4846567 (LYPLAL1) were nominally associated with IFG only in girls (all P < 0.05, Additional file 1: Table S3). However, only rs6795735 (ADAMTS9) showed marked differences in sex-specific effect size (P = 0.001, Additional file 1: Table S3). The effect of weighted and unweighted GRS showed no sexual dimorphism in Chinese children (all P > 0.05, Additional file 1: Table S3).
To our knowledge, this is the first study investigating association of the 11 GWAS derived WHR loci with FPG and risk of IFG in Chinese children. We identified three SNPs (rs6795735, rs984222 and rs1011731) showed individual effect on glycemic phenotypes. After Bonferroni correction, the individual effect of rs6795735 (ADAMTS9) on the risk of IFG remained significant.
ADAMTS9 encodes a member of the ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs) protein family, implicating in the cleavage of proteoglycans, the control of organ shape during development, and the inhibition of angiogenesis . ADAMTS9 was found to be involved in insulin signaling, which might have direct impact on glycemic phenotypes . The SNP rs6795735 located in intron ADAMTS9 was found to be associated with a 40.1% increased risk of IFG after additional adjustment for age- and sex-specific WHR-Z scores in our study (P = 0.002). This SNP was originally found to be associated with WHR in the GWAS . The original study also identified the SNP to be associated with high density lipoprotein cholesterol (Z score = − 2.486, P = 0.013), fasting glucose (Z score = 2.031, P = 0.042) and the risk of T2D (OR = 1.124, P = 0.002). Later in 2012, another GWAS conducted by the DIAbetes Genetics Replication and Meta-analysis (DIAGRAM) Consortium reported that rs6795735 was significantly associated with T2D . Afterwards, three studies explored the associations between rs6795735 and fat distribution or metabolic phenotypes, but did not report significant associations [7, 8, 27]. Our study is the first to validate the effect of rs6795735 on glycemic phenotypes in Chinese children. Our results suggest possible ethnic differences for effect of rs6795735 on glycemic phenotypes, demonstrating the value of conducting genetic epidemiology studies in populations with different ethnicity. However, more large-scaled association studies and functional studies are needed.
SNP rs984222 located near TBX15 and WARS2 was associated with WHR in the original study . Each G allele increased the WHR by 0.034 on average (P = 8.69 × 10− 25). However, no association between rs984222 and metabolic traits (including glycaemic traits and blood lipids) was identified . Subsequently, Burgdorf et al.  and Van Vliet-Ostaptchouk et al.  studied the effect of rs984222 on glycaemic traits in European populations. Hotta et al.  studied its impact on body fat distribution in Japanese. However, none of them found evidence supporting its association with these metabolic phenotypes. Therefore, the present study is the first to identify associations between rs984222 and glycemic phenotypes, independent of WHR. The inconsistent results may be explained by population differences. We conducted our study in children and adolescents, whose developmental and metabolic processes might have different characteristics compared with that of adults. Additionally, different ethnic backgrounds between Chinese and European population may partly account for the inconsistent effect of rs984222 on glycaemic traits.
The association of rs984222 and glycemic phenotypes could be explained by the evidence from functional studies of TBX15 gene. TBX15 gene encodes a transcription factor of a phylogenetically conserved family that regulates a variety of developmental processes. Tbx15 involved in craniofacial and limb development in the mouse [28, 29], and plays an important role in adipocyte differentiation and function . Gesta et al.  found that Tbx15 is highly differentially expressed between visceral and subcutaneous fat in both humans and rodents. Differential expression of Tbx15 between fat depots plays an important role in the interdepot differences in adipocyte differentiation, triglyceride accumulation, and mitochondrial function, which may contribute to the risk of diabetes and metabolic diseases .
Moreover, our study found that the WHR-decreasing allele of rs984222 was significantly associated with increased risk of IFG. This finding is somehow paradoxical. Similarly, van Vliet-Ostaptchouk et al.  found that the WHR-increasing allele at CPEB4 (rs6861681) was associated with lower fasting glucose after adjusting for WHR. The paradoxical associations have also been reported in other studies [7, 32, 33], which may be explained by the genetic pleiotropism, that is, some proteins play different roles in independent biological pathways . Studies had shown that some obese individuals remain metabolically healthy despite having excessive accumulation of body fat, while a number of normal-weight individuals are metabolically unhealthy [34, 35]. The mechanism causing the inconsistent results needs to be further clarified by functional studies and comprehensively assessments of genetic associations with BMI, WHR, glycemic phenotypes and T2D risk.
Dynamin 3 (DNM3) is a member of the dynamin family of enzymes. Dominant-negative mutant in transfected dynamin enzymes could promote GLUT6 and GLUT8 to be accumulated on the surface of rat adipose cells . PIGC encodes a subunit of the enzyme involving in lipid biosynthesis . In the present study, rs1011731 (DNM3-PIGC) was nominally associated with IFG, and the relationship was not significant after Bonferroni correction. However, in the original GWAS, rs1011731 was reported to be associated with only WHR, while not associated with FPG and T2D . The following replication studies didn’t found its association with glycaemic traits, either [7, 8]. More studies are required to investigate the association between rs1011731 and glycaemic traits, especially in East Asian population.
When analyses were conducted in a sex-stratified manner, we observed that the effect of rs6795735 (ADAMTS9) on IFG showed sexual dimorphism in Chinese children. In the original GWAS meta-analysis, seven of the 14 WHR-associated loci exhibited significant sexual dimorphism. And effect sizes were numerically greater in women than in men. Consistent with the original study, we also found that rs6795735 was associated with IFG only in girls, with effect sizes numerically greater in girls than in boys. Sex-specific genetic effects have long been acknowledged. However, the underlying molecular mechanism remains unclear. Some researches believed that the heterogeneity might due to the differences between boys and girls in body composition, patterns of weight gain, hormone biology, and the susceptibility to certain genetic and environmental factors . More studies are required before a clear understanding of the underlying molecular mechanisms.
It was known that WHR is a measure of central obesity and associated with diabetes and prediabetes (including impaired glucose tolerance (IGT) and IFG) [6, 38]. Subjects with IFG have a prominent deficiency in insulin secretion, suggesting an association with insufficient pancreatic B-cell function . We studied the correlations between WHR-related loci found by previous GWAS and WHR, BMI, or glycemic phenotypes (IFG and FPG). Although we failed to confirm the associations between these loci and WHR or BMI, we identified individual effects of the 3 SNPs (rs984222, rs6795735 and rs1011731) on glycemic phenotypes. The different findings between our study and previous studies may be partly explained by the fact that the pattern of fat distribution in Chinese differs from those in Europeans, i.e. Chinese individuals have excessive body fat at normal BMI or WHR [10,11,12]. In addition, we supposed that genetic predisposition to central obesity might be associated with risk of impaired pancreatic islet function, and hypothesized the effect of these loci on glycemic phenotypes might be relevant at earlier stage than that on WHR in our children sample. It had been shown that both insulin resistance and visceral fat deposition (especially liver and epicardial fat) could be the trigger of the metabolic syndrome, which may cause glycemic alteration ahead of measurable WHR changes .
One highlight of the present study was that our study sample was children. Many common complex diseases observed in adults have their developmental origin in childhood [41, 42]. Therefore, it is of great significance to understand how the WHR variants derived from adults operate in children and if they confer risk for childhood obesity and other obesity-related diseases.
The present study had limitations. Firstly, we selected several WHR related genetic loci which were recently identified by the GWAS by Heid and colleagues . There were other WHR related loci reported since the publication of Heid et al. [43,44,45,46], and the addition of more WHR variants to the GRS models would make our results more accurate. However, we think this study could help to provide evidences for the effects of these variants identified by Heid and colleagues. Secondly, the sample size was limited. However, the power of rs6795735, which showed significant association with IFG after adjustment for Bonferroni correction, was 92%, suggesting that the associations might not be false-positive.
In conclusion, we found individual effects of rs6795735 (ADAMTS9), rs984222 (TBX15-WARS2), and rs1011731 (DNM3-PIGC) on glycemic phenotypes. These associations were independent of WHR. Our findings indicated that genetic susceptibility to central obesity might be associated with risk of impaired pancreatic islet function, providing more evidence for the pathogenesis of diabetes.
- ADAMTS9 :
A disintegrin-like and metalloproteinase with thrombospondin type 1 motif
Adolescent lipids, insulin resistance and candidate genes
Body mass index
Comprehensive Prevention project for overweight and Obese Adolescents
- DNM3 :
Fasting plasma glucose
- F ST :
Genetic risk score
Genome-wide association studies
Hardy-Wei nberg equilibrium
Impaired fasting glucose
Impaired glucose tolerance
Polymerase chain reaction
- PIGC :
Phosphatidylinositol glycan anchor biosynthesis, class C
School-based Comprehensive Intervention on Childhood Obesity
Single nucleotide polymorphism
Type 2 diabetes.
- TBX15 :
- WARS2 :
Tryptophanyl-tRNA synthetase 2
Zhao X, Xi B, Shen Y, Wu L, Hou D, Cheng H, et al. An obesity genetic risk score is associated with metabolic syndrome in Chinese children. Gene. 2014;535:299–302.
Ng MC, Hester JM, Wing MR, Li J, Xu J, Hicks PJ, et al. Genome-wide association of BMI in African Americans. Obesity (Silver Spring). 2012;20:622–7.
Song Y, Ma J, Wang HJ, Wang Z, Hu P, Zhang B, et al. Secular trends of obesity prevalence in Chinese children from 1985 to 2010: urban-rural disparity. Obesity (Silver Spring). 2015;23:448–53.
Kong X, Zhang X, Zhao Q, He J, Chen L, Zhao Z, et al. Obesity-related genomic loci are associated with type 2 diabetes in a Han Chinese population. PLoS One. 2014;9:e104486.
Lindgren CM, Heid IM, Randall JC, Lamina C, Steinthorsdottir V, Qi L, et al. Genome-wide association scan meta-analysis identifies three loci influencing adiposity and fat distribution. PLoS Genet. 2009;5:e1000508.
Heid IM, Jackson AU, Randall JC, Winkler TW, Qi L, Steinthorsdottir V, et al. Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution. Nat Genet. 2010;42:949–60.
Burgdorf KS, Gjesing AP, Grarup N, Justesen JM, Sandholt CH, Witte DR, et al. Association studies of novel obesity-related gene variants with quantitative metabolic phenotypes in a population-based sample of 6,039 Danish individuals. Diabetologia. 2012;55:105–13.
van Vliet-Ostaptchouk JV, den Hoed M, Luan J, Zhao JH, Ong KK, van der Most PJ, et al. Pleiotropic effects of obesity-susceptibility loci on metabolic traits: a meta-analysis of up to 37,874 individuals. Diabetologia. 2013;56:2134–46.
Huang T, Qi QB, Zheng Y, Ley SH, Manson JE, Hu FB, et al. Genetic predisposition to central obesity and risk of type 2 diabetes: two independent cohort studies. Diabetes Care. 2015;38:1306–11.
Lear SA, James PT, Ko GT, Kumanyika S. Appropriateness of waist circumference and waist-to-hip ratio cutoffs for different ethnic groups. Eur J Clin Nutr. 2010;64:42–61.
Wang J, Thornton JC, Russell M, Burastero S, Heymsfield S, Pierson RJ. Asians have lower body mass index (BMI) but higher percent body fat than do whites: comparisons of anthropometric measurements. Am J Clin Nutr. 1994;60:23–8.
Deurenberg P, Deurenberg YM, Wang J, Lin FP, Schmidt G. The impact of body build on the relationship between body mass index and percent body fat. Int J Obes Relat Metab Disord. 1999;23:537–42.
Wang HJ, Ma J, Yin ZD, Xu YQ, Pan YP. Study on the relationship between obesity and metabolic syndrome in adolescents. Chinese Journal of School Health. 2007;28:878–82.
Wang HJ, Zhang H, Zhang SW, Pan YP, Ma J. Association of the common genetic variant upstream of INSIG2 gene with obesity related phenotypes in Chinese children and adolescents. Biomed Environ Sci. 2008;21:528–36.
Wang D, Ma J, Zhang S, Hinney A, Hebebrand J, Wang Y, et al. Association of the MC4R V103I polymorphism with obesity: a Chinese case-control study and meta-analysis in 55,195 individuals. Obesity (Silver Spring). 2010;18:573–9.
Ma J, Zhang SW, Wu SX, Zhang H, Li SS. Prevalence of metabolic syndrome among children with different nutritional status in Beijing. Chinese Journal of School Health. 2008;29:103–5.
Ji CY. Report on childhood obesity in China (1)--body mass index reference for screening overweight and obesity in Chinese school-age children. Biomed Environ Sci. 2005;18:390–400.
National Health and Family Planning Commission of PRC. WS/T 456–2014, Screening standard for malnutrition of school-age children and adolescents. China: NHFPC, 2014. In.
Genuth S, Alberti KG, Bennett P, Buse J, Defronzo R, Kahn R, et al. Follow-up report on the diagnosis of diabetes mellitus. Diabetes Care. 2003;26:3160–7.
Wang HJ, Hinney A, Song JY, Scherag A, Meng XR, Grallert H, et al. Association of common variants identified by recent genome-wide association studies with obesity in Chinese children: a case-control study. BMC Med Genet. 2016;7
Meng XR, Song JY, Ma J, Liu FH, Shang XR, Guo XJ, et al. Association study of childhood obesity with eight genetic variants recently identified by genome-wide association studies. Pediatr Res. 2014;76:310–5.
Duan S, Zhang W, Cox NJ, Dolan ME. FstSNP-HapMap3: a database of SNPs with high population differentiation for HapMap3. Bioinformation. 2008;3:139–41.
Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.
Balloux F, Lugon-Moulin N. The estimation of population differentiation with microsatellite markers. Mol Ecol. 2002;11:155–65.
Clark ME, Kelner GS, Turbeville LA, Boyer A, Arden KC, Maki RA. ADAMTS9, a novel member of the ADAM-TS/ metallospondin gene family. Genomics. 2000;67:343–50.
Morris AP, Voight BF, Teslovich TM, Ferreira T, Segre AV, Steinthorsdottir V, et al. Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Genet. 2012;44:981–90.
Hotta K, Kitamoto A, Kitamoto T, Mizusawa S, Teranishi H, So R, et al. Replication study of 15 recently published loci for body fat distribution in the Japanese population. J Atheroscler Thromb. 2013;20:336–50.
Singh MK, Petry M, Haenig B, Lescher B, Leitges M, Kispert A. The T-box transcription factor Tbx15 is required for skeletal development. Mech Dev. 2005;122:131–44.
Agulnik SI, Papaioannou VE, Silver LM. Cloning, mapping, and expression analysis of TBX15, a new member of the T-box gene family. Genomics. 1998;51:68–75.
Gesta S, Bezy O, Mori MA, Macotela Y, Lee KY, Kahn CR. Mesodermal developmental gene Tbx15 impairs adipocyte differentiation and mitochondrial respiration. Proc Natl Acad Sci U S A. 2011;108:2771–6.
Gesta S, Bluher M, Yamamoto Y, Norris AW, Berndt J, Kralisch S, et al. Evidence for a role of developmental genes in the origin of obesity and body fat distribution. Proc Natl Acad Sci U S A. 2006;103:6676–81.
Okada Y, Kubo M, Ohmiya H, Takahashi A, Kumasaka N, Hosono N, et al. Common variants at CDKAL1 and KLF9 are associated with body mass index in east Asian populations. Nat Genet. 2012;44:302–6.
Zhu J, Zong G, Lu L, Gan W, Ji L, Hu R, et al. Association of genetic predisposition to obesity with type 2 diabetes risk in Han Chinese individuals. Diabetologia. 2014;57:1830–3.
Pajunen P, Kotronen A, Korpi-Hyovalti E, Keinanen-Kiukaanniemi S, Oksa H, Niskanen L, et al. Metabolically healthy and unhealthy obesity phenotypes in the general population: the FIN-D2D survey. BMC Public Health. 2011;11:754.
Stefan N, Kantartzis K, Machann J, Schick F, Thamer C, Rittig K, et al. Identification and characterization of metabolically benign obesity in humans. Arch Intern Med. 2008;168:1609–16.
Lisinski I, Schurmann A, Joost HG, Cushman SW, Al-Hasani H. Targeting of GLUT6 (formerly GLUT9) and GLUT8 in rat adipose cells. Biochem J. 2001;358:517–22.
Wisniewski AB, Chernausek SD. Gender in childhood obesity: family environment, hormones, and genes. Gend Med. 2009;6(Suppl 1):76–85.
Siddiquee T, Bhowmik B, Karmaker RK, Chowdhury A, Mahtab H, Azad KA, et al. Association of general and central obesity with diabetes and prediabetes in rural Bangladeshi population. Diabetes Metab Syndr. 2015;9:247–51.
Hong J, Gui MH, Gu WQ, Zhang YF, Xu M, Chi ZN, et al. Differences in insulin resistance and pancreatic B-cell function in obese subjects with isolated impaired glucose tolerance and isolated impaired fasting glucose. Diabet Med. 2008;25:73–9.
Avogaro A. Insulin resistance: trigger or concomitant factor in the metabolic syndrome. Panminerva Med. 2006;48:3–12.
Barker DJ. Sir Richard Doll Lecture Developmental origins of chronic disease. Public Health. 2012;126:185–9.
Gluckman PD, Hanson MA, Cooper C, Thornburg KL. Effect of in utero and early-life conditions on adult health and disease. N Engl J Med. 2008;359:61–73.
Wen W, Kato N, Hwang JY, Guo X, Tabara Y, Li H, et al. Genome-wide association studies in east Asians identify new loci for waist-hip ratio and waist circumference. Sci Rep. 2016;6:17958.
Shungin D, Winkler TW, Croteau-Chonka DC, Ferreira T, Locke AE, Magi R, et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature. 2015;518:187–96.
Yoneyama S, Guo Y, Lanktree MB, Barnes MR, Elbers CC, Karczewski KJ, et al. Gene-centric meta-analyses for central adiposity traits in up to 57 412 individuals of European descent confirm known loci and reveal several novel associations. Hum Mol Genet. 2014;23:2498–510.
Liu CT, Monda KL, Taylor KC, Lange L, Demerath EW, Palmas W, et al. Genome-wide association of body fat distribution in African ancestry populations suggests new loci. PLoS Genet. 2013;9:e1003681.
The study was supported by grants from Beijing Natural Science Foundation [grant number: 7162106] and National Natural Science Foundation of China [grant number: 81573170]. We thank all the children and adolescents and their parents for their participation.
The study was funded by grants from Beijing Natural Science Foundation [grant number: 7162106] and National Natural Science Foundation of China [grant number: 81573170].
Availability of data and materials
The datasets generated and/or analysed during the current study are not publicly available, since ethics approval and participants’ consent does not allow public sharing of data, but are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Both the ALIR study and CPOOA study were approved by the Ethics Committee of Peking University Health Science Center. Written informed consent was provided by all participants and, in the case of minors, their parents. Studies were performed according to the Declaration of Helsinki.
Consent for publication
Written informed consent was provided by all participants and, in the case of minors, their parents.
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. General characteristics of study subjects with and without IFG, Table S2. Associations of 11 SNPs with WHR and BMI in 2030 Chinese children, Table S3 Sex-dependent association of 11 SNPs and GRS with IFG in 2030 Chinese children. (DOCX 25 kb)