Skip to main content

Association between genetic risk score and tri-ponderal mass index growth trajectories among different dietary consumption adolescents in a prospective Taiwanese cohort



Single-nucleotide polymorphisms (SNPs) in various genetic loci are associated with childhood obesity; however, their influence on adolescent growth patterns has rarely been explored. This study investigated whether genetic variants could predict tri-ponderal mass index (TMI)-derived growth trajectories and the interaction between genetic and dietary factors.


We conducted Taiwan Puberty Longitudinal Study, a prospective cohort that recruited 1,135 children since 2018. Anthropometric measurements were recorded every three months, while dietary nutrition assessment and biological sampling for genotyping were collected during the first visit. TMI growth trajectory groups were identified using growth mixture modeling. A multinomial logistic regression model for different growth trajectories was used to examine the effect of candidate SNPs, and the most related SNPs were used to establish the genetic risk score. We then explored the effect of the genetic risk score in subgroup analysis according to dietary calories and different dietary consumption patterns.


Three TMI-based growth trajectory groups were identified among adolescents. The “increased weight” trajectory group accounted for approximately 9.7% of the participants. FTO/rs7206790 was associated with the increased weight growth trajectory after adjusting for the baseline TMI and other correlated covariates (OR: 2.13, 95% CI: 1.08–4.21). We generated the genetic risk score using 4 SNPs (FTO/rs7206790, ADCY9/rs2531995, TFAP2B/rs4715210, and TMEM18/rs6548238) and selected the threshold of 10 points to define risk categories. There were 11.66% and 3.24% of participants belonged to the increased weight trajectory in high- and low-risk groups, respectively; and the predictive ability of the genetic risk score was notable among low calories intake participants (OR: 1.90, 95% CI: 1.18–3.05 vs. OR: 1.17, 95% CI: 0.78–1.75 in high calories intake group).


Our results offer a new perspective on the genetic and dietary basis of changes in adolescent obesity status. Individualized interventions for obesity prevention may be considered among high-risk children.


Childhood obesity is a major public health concern with increasing global prevalence. It is associated with many noncommunicable diseases, such as type II diabetes, hypertension, dyslipidemia, and cardiovascular disease [1, 2]. However, some studies have found that individuals who are obese in adulthood but had a below-average weight during childhood have the highest risk of some cardiovascular diseases [3, 4]. Because adolescence is associated with rapid body development, obesity growth trajectories during this period may play a pivotal role in increasing the risk of cardiometabolic dysfunction [5, 6].

Body mass index (BMI) is a simple and common index for defining overweight and obesity in adults. However, the use of this index in adolescents and children is less convenient because the growth chart for the same sex and age is required as a reference. 7 The tri-ponderal mass index (TMI), which is defined as the weight (kg)/height (m3), was proposed as a better indicator of adolescent body fat composition [8]. It is useful for measuring overweight, and obesity in adolescents in clinical practice, and a persistent increase in the TMI growth trajectory during adolescence can predict diabetes in early adulthood more accurately than can BMI [9].

Obesity is believed to result from the interplay of the environment and innate genetic factors [10, 11]. Similar to obesity, several environmental factors are associated with the growth trajectory in childhood, including parental education level [12], a decrease in moderate to vigorous physical activity [13], and early puberty [14]. However, these factors may not completely explain the different growth patterns in adolescents. Researchers have identified more than 1,000 possible obesity-related loci [15] after the publication of the first obesity-related genome-wide association study (GWAS) in 2007, which demonstrated BMI-associated single-nucleotide polymorphisms (SNPs) in the FTO gene [16]. Nevertheless, no GWAS study has been conducted to explore the association between SNPs and TMI growth trajectories during adolescence.

We conducted a prospective longitudinal study to construct TMI growth trajectories among adolescent participants and then investigated the association between 17 obesity-related genetic polymorphisms and the different trajectory groups. Studies have shown that genetic influences on obesity may be modified by dietary fat intake [17, 18]; therefore, we explored whether the effect of genetic variants is altered by the participants’ different proportions of calories and main nutrient intake.

Materials and methods

Study population

The Taiwan Puberty Longitudinal Study (TPLS) is a multidisciplinary, longitudinal project. This project recruited children from pubertal and pediatric endocrine clinics of several hospitals [Taipei Medical University Hospital, Cathay General Hospital (CGH), Taipei Municipal Wanfang Hospital, and National Cheng Kung University Hospital (NCKUH)] beginning in 2018. Girls aged 6–14 years and boys aged 9–17 years were invited to participate and received prospective follow-ups. The exclusion criteria were (1) the presence of hereditary diseases, such as Turner syndrome and Prader–Willi syndrome; (2) the presence of neurological and psychiatric diseases, such as Tourette syndrome or attention-deficit/hyperactivity disorder; and (3) intake of medicines that may affect natural growth during adolescence, such as leuplin (leuprorelin acetate), oral steroids, and ritalin (methylphenidate). Parents or guardians had to complete an informed consent form at the baseline visit. This study was approved by the Institutional Review Board of Taipei Medical University (N201802018), that of CGH (CGH-P108107), and that of NCKUH (B-BR-108-076) and complied with the principles outlined in the Helsinki Declaration.

Weight index measurement and other covariates

The participants’ weight and height were measured using an electronic scale, with shoes taken off, light clothing worn, and participants fasted beforehand. TMI was calculated as weight (kg)/height (m3), and the same measurement was performed at the first visit and then repeated every 3 months. Each participant’s parents or guardians were asked to answer a baseline characteristics questionnaire along with their children, which included questions related to the participants’ birth date, gender, postnatal information (birth weight and breastfeeding condition), family socioeconomic status (family income and parental education level), and environmental factors (home tobacco smoke exposure and incense burning). Moreover, a Chinese version of the International Physical Activity Questionnaire was used to calculate the average exercise time for all children [19].

At the first visit, a pediatric endocrinologist assessed the pubic hair and breast development of each girl and measured testicular volume using the Prader orchidometer for each boy [20]. If the two breasts of a girl were in different developmental stages or the testes of a boy were not identical, the stage of the more advanced side was adopted. Finally, pubertal development was then graded according to the 5-stage scale described by Tanner [21, 22].

Nutrient intake assessment

Twenty-four-hour dietary recall was conducted by trained registered dietitians on two different days after the first recruitment. Energy and nutrient intakes were estimated on the basis of Nutritionist Edition, COFIT Pro, Version 1.0.0, a software package for nutrient analysis that features a Taiwanese food composition table as the nutrient database (Taipei, Taiwan). A validation study of COFIT has been published by our team [23]. The dietary records were coded and linked to this system to calculate nutrient and energy intake.

Candidate SNP genotyping

Participants were requested not to eat or drink for at least 1 h before scraping the brush inside the mouth 10 times. The collection brush was then maintained at − 80 °C and transferred for DNA purification [24]. We applied the DNA extraction protocol using the Gentra Puregene Buccal Cell Kit 140 [24].

Genotyping was performed using the Sequenom MassARRAY iPLEX platform [25] at the National Center for Genome Medicine, Taiwan. In total, 17 candidate obesity-related SNPs identified from our previous unpublished study in a Taiwanese children cohort were selected, and most of them have also been reported to be associated with obesity from literature review (rs574367 [26], rs9356744 [26], rs12597579 [26], rs6567160 [26], rs3817334 [26], rs1555543 [26], rs4715210 [26], rs7206790 [27], rs6548238 [28], rs9939609 [29], rs1421085 [30], rs12463617 [31], rs8053360 [32], rs16858082 [33], rs7498665 [34], and rs4788102 [35]) except rs2531995. After polymerase chain reaction amplification, shrimp alkaline phosphatase incubation, primer extension, and chip dispensing, data acquisition and automated genotype calling were performed through matrix-assisted laser desorption/ionization time-of-flight mass spectrometry.

Statistical analysis

The TMI growth trajectories were constructed with growth mixture modeling (GMM) by using Mplus software version 7. GMM was conducted to explore an individual’s pattern of weight change with latent trajectory classes. It allowed different trajectory groups to vary around different means and let each latent class have its estimates of variances and covariate influences [36]. The participants’ TMI from every trimonthly visit was adopted for growth trajectory identification. Using the full-information maximum likelihood method, missing values were managed in GMM by assuming them to be missing at random [37]. We tested the number of trajectories from two to five and examined the trajectory pattern as linear, quadratic, or cubic to find a suitable model. The Bayesian information criterion (BIC) and the bootstrap likelihood ratio test (BLRT) were used to determine the model’s adequacy. We ensured at least 5% of participants were sorted to each group.

For descriptive analyses, baseline participant characteristics in each growth trajectory group were compared using the chi-square test and analysis of variance for categorical and continuous variables, respectively. Multinomial logistic regression was used and odds ratios (ORs) of genetic variants for specific growth trajectories were reported after adjustment for age, sex, and baseline confounding covariates. Although the parents or guardians of some participants omitted the information concerning family income and education level, we still included these participants with missing data in our analysis using the missing indicators method [38].

To explore the combined effects of obesity-related SNPs on growth trajectory, an additive model was used to determine the genetic risk score (GRS). We selected risk alleles with a p of ≤ 0.15 in regression analysis and weighted them according to their coefficients in the regression model [39]. We then used logistic regression to evaluate the effect of the GRS and the predictive ability among different cut-off points.

For gene-diet interaction investigation, we used all participants’ average daily dietary calories and average percentage of calories from carbohydrate, fat and protein as the cut-off points to define different dietary consumption patterns. The same regression method was applied to explore whether different associations existed between the GRS and increased weight growth trajectory among high- or low-calories/nutrients intake subgroups. Statistical significance was defined by a two-tailed p of < 0.05. All analyses were performed using SAS version 9.3 (SAS Institute).


A total of 1,135 participants in the TPLS, for whom follow-up data for > 6 months were available, were used to establish growth trajectory groups (average follow-up: 4.4 times, 414 days). After excluding the model with < 5% participants in each trajectory and conducting an adequacy assessment by using the BIC and BLRT, we determined three TMI growth trajectory groups (Fig. 1). Among the three groups, most children were categorized as the “stable weight” group (n = 934, 82.3%) and demonstrated almost no change in the TMI during the follow-up period. The TMI decreased in approximately 8.0% (n = 91) of the participants, whereas it increased in the remaining 9.7% (n = 110) children during the observation period.

Fig. 1
figure 1

Three TMI-derived growth trajectory groups identified using growth mixture modeling

Table 1 presents the basic characteristics of each trajectory group. In the stable weight growth pattern group, the baseline TMI was significantly lower and the daily calorie intake per kilogram of weight was higher than that in the other two groups. The participants in the “increased weight” group appeared to have a lower socioeconomic status, which is consistent with the results of a previous study [12]. Concerning pubertal development, the probability of the early Tanner stage was higher in the “decreased weight” group and lower in the increased weight group.

Table 1 Basic characteristics of three TMI-based growth trajectory groups in the Taiwan Puberty Longitudinal Study

We then investigated the effect of the 17 candidate SNPs and found only 1 SNP (FTO rs7206790) to be significantly associated with the increased weight growth trajectory (adjusted OR: 2.13, 95% confidence interval (CI): 1.08–4.21, p-value: 0.03 for additive model; Table 2). Because there were missing information for some participants, for example, around 16.6% participants didn’t provide parental educational level, and the data of family income were not available in 17.4% participants (Table 1), we performed a sensitivity analysis that excluded all participants with missing variables. The adjusted OR was 2.46 for FTO rs7206790 for increase weight growth trajectory, though the p-value was not siginificant due to insufficient power (p-value: 0.54, total participants number: 743). For GRS construction, we selected 4 genes with p ≤ 0.15 in regression analysis to generate the GRS (FTO/ rs7206790, ADCY9/ rs2531995, TFAP2B/ rs4715210, and TMEM18/ rs6548238; Table 2). The number of risk alleles was counted for each SNP, and rs7206790 and rs6548238 were double-weighted according to the coefficients in regression analysis to generate a GRS ranging from 0 to 12. We defined GRS ≥ 10 as the threshold of the high risk group because of the best discrimination ability (Table 3), and 11.66% and 3.24% of participants belonged to the increased weight trajectory in the high- and low-risk groups, respectively (p < 0.001). In addition, we also used 5-fold cross-validation method to test the GRS, and the result is consistent.

Table 2 Association of 17 candidate SNP alleles with the increased weight growth trajectory in the Taiwan Puberty Longitudinal Study
Table 3 The performance of different threshold of GRS to predict increased weight TMI growth trajectory

For the dietary calories and nutrient intake subgroup analysis, the cut-off points for daily calories, carbohydrate, fat and protein intake were 1598.7 Kcal, 47.7%, 33.8% and 15.5%, respectively. Though no significant gene–diet interaction was found, our GRS exhibited relatively different effects between the paired subgroups. For instance, the GRS was associated with an increased risk of the increased weight trajectory in children consuming on a low-calories diet (OR: 1.90, 95% CI: 1.18–3.05, p-value: 0.01), but the influence on children with high-calories intake was not significant (OR: 1.17, 95% CI: 0.78–1.75, p-value: 0.44; Table 4; Fig. 2).

Fig. 2
figure 2

The performance of genetic risk score among different dietary consumption subgroups

Table 4 Odds ratios of GRS for increased weight TMI growth trajectory in different nutrients intake subgroups


By using the GMM statistical method, we categorized the TPLS participants into three TMI-based growth trajectory groups: stable weight, decreased weight, and increased weight. After adjusting for gender, age, baseline TMI and other covariates, rs7206790 showed a significant association with the increased weight growth trajectory. Our study provided evidence for the association between the GRS of the four SNPs (namely rs7206790, rs2531995, rs4715210, and rs6548238) and the increased weight trajectory. Children with GRS ≥ 10 had a significantly increased risk to become obesity during adolescence compared with those with GRS < 10, and the influence seemed more remarkable for those with low calories intake.

Although many obesity-related GWAS have been conducted and have identified many possible genetic loci, relatively few studies have been designed to examine the genetic influence on the change of obese status in children. Through a systematic review of nine twin and five adopted children studies, Silventoinen et al. found that genetic factors had a considerably stronger influence on BMI variation than environmental factors did up to the age of 18 years [11]. Similar to the studies on adults, the strongest association was found in the SNPs of the FTO gene. More than 10 studies with different designs have explored the effects of genes on adolescence obesity since [40]. Warrington et al. have detected that the A allele at FTO/ rs1558902 was associated with the decrease in BMI in infancy, but showed an inverse effect from early childhood [41]. Besides, Graff et al. also found that the association between BMI and obesity-related loci, such as FTO, TFAP2B, and TMEM18, varied by age during adolescence and early adulthood, meaning that genetic consequences may change over time [42]. The genetic finding of this study was compatible with our GRS components, although the dominant SNPs may differ because of the change in ethnic groups.

Compared with normal body weight adolescents, individuals who are obese during pubertal growth have a significantly higher tendency to develop severe obesity in adulthood [43]. Moreover, our previous finding showed that a persistent increase in the TMI during adolescence is an independent risk factor for early-onset diabetes adjusting for childhood obesity status [9]. Obesity development in adolescents is not simply explained by energy balance disruption but may be a complex condition involving genetic, environmental, biological, and behavioral elements. Public obesity prevention strategies are mainly focused on the healthy promotion dimension, such as dietary habit modulation and encouraging physical activity. Because the increase in health awareness and lifestyle modification has been shown to be effective during the adolescent period when school-based interventions are applied [44] these approaches can be readily implemented during adolescence. Our finding that participants with a higher GRS would more easily develop obesity during adolescence, independent of baseline obesity status, is helpful for differentiating young children who are at a higher risk of obesity. Furthermore, the relationship between energy intake and obesity is more complex during adolescence, and our result supports the hypothesis that the genetic influence was more notable among the children with low dietary calories. Our research could help prevent the metabolic disorder problems resulting from obesity despite of low dietary caloriesby way of other lifestyle modification, should the genetic influence ring the alarm bell in the first place.

The main strength of our study is that anthropometric measurements were conducted every 3 months, and thus, we could easily capture different obesity development patterns during early adolescence and construct the growth trajectory. Moreover, we collected several possible confounding factors for the adolescent growth trajectory (including baseline obesity status, pubertal development, and environmental and behavioral factors) from the baseline questionnaire and physical examination, which enhances the accuracy of our results. Furthermore, few researchers have investigated the associations between obesity susceptibility variants and the change in obesity status across childhood. To our best knowledge, ours is the first study to investigate genetic influences on the TMI growth trajectory.

A limitation of our study is that most confounding factors were evaluated at baseline, and these conditions, especially exercise habits, may change and produce information bias. For example, Kwon et al. found that compared with consistently active participants, children who were active at first but decreased moderately to vigorous- intensity physical activity with age were more likely to become obese (OR: 2.77, 97% CI: 1.16–6.58) [13]. Academic performance is highly valued in Taiwanese high school students. In one study that assessed their physical activity, only 5.4% of Taiwanese adolescents (age: 13–15 years) met the World Health Organization recommendation for physical activity each day [45]. Because the overall physical activity performance of Taiwanese adolescents is suboptimal, the influence on the association between exercise habits and growth trajectory patterns may be tolerable. Besides, we used 24-hour dietary record to capture the nutrient intake of children. Even we have collected the information from 2 different days, it might be challenging to evaluate the long-term dietary consumption for adolescents without the follow-up longitudinal data. Moreover, the observation period of growth trajectory was not sufficiently long to distinguish more distinct growth trajectory groups. Nevertheless, the proportion of adolescents in the increased weight TMI growth trajectory group in this study was almost equal to the sum of children in the two different weight gain trajectory groups established using 6-year follow-up records in another independent cohort (9.7% vs. 9.3%) [9]. This result indicates that our TMI trajectory model was relatively stable during adolescence, and our GRS was a suitable tool for predicting children who were at a higher risk of obesity. Finally, we only recruited Taiwanese adolescents in the TPLS; whether the finding of our GRS could be applied to other ethnic populations needs to be confirmed in further studies.

In conclusion, FTO rs7206790 was significantly associated with an increase in the TMI during the adolescent period, and our 4-SNP GRS approach is a practical tool for determining cumulative genetic susceptibility and predicting the increased weight TMI growth trajectory, especially in the children with lower calories intake. However, further research with a larger sample size from a different population is necessary to confirm the gene–diet interaction as it concerns the adolescent growth trajectory.

Availability of data and materials

All of the data are available with reasonable request from the corresponding author.



Bayesian information criterion


bootstrap likelihood ratio test


Body mass index


confidence interval


growth mixture modeling


genetic risk score


genome-wide association study


odds ratio


Single-nucleotide polymorphism


tri-ponderal mass index


Taiwan Puberty Longitudinal Study


  1. Gurnani M, Birken C, Hamilton J. Childhood obesity: causes, Consequences, and management. Pediatr Clin North Am. 2015;62:821–40.

    Article  Google Scholar 

  2. Sahoo K, Sahoo B, Choudhury AK, Sofi NY, Kumar R, Bhadoria AS. Childhood obesity: causes and consequences. J Family Med Prim Care. 2015;4:187–92.

    Article  Google Scholar 

  3. Abraham S, Collins G, Nordsieck M. Relationship of childhood weight status to morbidity in adults. HSMHA Health Rep. 1971;86:273–84.

    Article  CAS  Google Scholar 

  4. Bjerregaard LG, Jensen BW, Angquist L, Osler M, Sorensen TIA, Baker JL. Change in overweight from childhood to early adulthood and risk of type 2 diabetes. N Engl J Med. 2018;378:1302–12.

    Article  Google Scholar 

  5. Attard SM, Herring AH, Howard AG, Gordon-Larsen P. Longitudinal trajectories of BMI and cardiovascular disease risk: the national longitudinal study of adolescent health. Obesity. 2013;21:2180–8.

    Article  Google Scholar 

  6. Norris T, Mansukoski L, Gilthorpe MS, et al. Distinct body mass index trajectories to young-adulthood obesity and their different cardiometabolic consequences. Arterioscler Thromb Vasc Biol. 2021;41:1580–93.

    Article  CAS  Google Scholar 

  7. Barlow SE, Expert C. Expert committee recommendations regarding the prevention, assessment, and treatment of child and adolescent overweight and obesity: summary report. Pediatrics. 2007;120(Suppl 4):164-92.

    Google Scholar 

  8. Peterson CM, Su H, Thomas DM, et al. Tri-ponderal mass index vs body mass index in estimating body fat during adolescence. JAMA Pediatr. 2017;171:629–36.

    Article  Google Scholar 

  9. Wu YF, Fan HY, Chen YC, Kuo KL, Chien KL. Adolescent tri-ponderal mass index growth trajectories and incident diabetes mellitus in early adulthood. J Clin Endocrinol Metab. 2021;106:e2919–27.

    Article  Google Scholar 

  10. Stunkard AJ, Harris JR, Pedersen NL, McClearn GE. The body-mass index of twins who have been reared apart. N Engl J Med. 1990;322:1483–7.

    Article  CAS  Google Scholar 

  11. Silventoinen K, Rokholm B, Kaprio J, Sorensen TI. The genetic and environmental influences on childhood obesity: a systematic review of twin and adoption studies. Int J Obes (Lond). 2010;34:29–40.

    Article  CAS  Google Scholar 

  12. Ng CD, Cunningham SA. In, out, and fluctuating: obesity from adolescence to adulthood. Ann Epidemiol. 2020;41:14–20.

    Article  Google Scholar 

  13. Kwon S, Janz KF, Letuchy EM, Burns TL, Levy SM. Active lifestyle in childhood and adolescence prevents obesity development in young adulthood. Obesity. 2015;23:2462–9.

    Article  CAS  Google Scholar 

  14. Viner RM, Costa S, Johnson W. Patterns of BMI development between 10 and 42 years of age and their determinants in the 1970 british Cohort Study. J Epidemiol Commun Health. 2019;73:79–85.

    Article  Google Scholar 

  15. Buniello A, MacArthur JAL, Cerezo M, et al. The NHGRI-EBI GWAS catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019;47:D1005-D12.

    Article  Google Scholar 

  16. Scuteri A, Sanna S, Chen WM, et al. Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits. PLoS Genet. 2007;3:e115.

    Article  Google Scholar 

  17. Sanchez-Moreno C, Ordovas JM, Smith CE, Baraza JC, Lee YC, Garaulet M. APOA5 gene variation interacts with dietary fat intake to modulate obesity and circulating triglycerides in a Mediterranean population. J Nutr. 2011;141:380–5.

    Article  CAS  Google Scholar 

  18. Alsulami S, Nyakotey DA, Dudek K, et al. Interaction between metabolic genetic risk score and dietary fatty acid intake on central obesity in a Ghanaian Population. Nutrients. 2020;2:12.

    Google Scholar 

  19. Qu NN, Li KJ. [Study on the reliability and validity of international physical activity questionnaire (chinese vision, IPAQ)]. Zhonghua liu xing bing xue za zhi = Zhonghua. liuxingbingxue zazhi. 2004;25:265–8.

    Google Scholar 

  20. Zachmann M, Prader A, Kind HP, Hafliger H, Budliger H. Testicular volume during adolescence. Cross-sectional and longitudinal studies. Helv Paediatr Acta. 1974;29:61–72.

    CAS  Google Scholar 

  21. Tanner JM, Whitehouse RH. Clinical longitudinal standards for height, weight, height velocity, weight velocity, and stages of puberty. Arch Dis Child. 1976;51:170–9.

    Article  CAS  Google Scholar 

  22. Marshall WA, Tanner JM. Variations in pattern of pubertal changes in girls. Arch Dis Child. 1969;44:291–303.

    Article  CAS  Google Scholar 

  23. Wang JS, Hsieh RH, Tung YT, Chen YH, Yang C, Chen YC. Evaluation of a technological image-based dietary assessment tool for children during pubertal growth: a pilot study. Nutrients. 2019;2:11.

    Google Scholar 

  24. Nedel F, Andre Dde A, de Oliveira IO, Tarquinio SB, Demarco FF. Buccal cells submitted to three different storage conditions before DNA extraction. J Appl Oral Sci. 2009;17:113–5.

    Article  CAS  Google Scholar 

  25. Gabriel S, Ziaugra L, Tabbaa D. SNP genotyping using the Sequenom MassARRAY iPLEX platform. Curr Protoc Hum Genet. 2009;2:12.

    Google Scholar 

  26. Lu Y, Loos RJ. Obesity genomics: assessing the transferability of susceptibility loci across diverse populations. Genome Med. 2013;5:55.

    Article  Google Scholar 

  27. Xu Y, Ling J, Yang M, et al. Rs7206790 and rs11644943 in FTO gene are associated with risk of obesity in chinese school-age population. PLoS ONE. 2014;9:e108050.

    Article  Google Scholar 

  28. Dorajoo R, Ong RT, Sim X, et al. The contribution of recently identified adult BMI risk loci to paediatric obesity in a singaporean chinese childhood dataset. Pediatr Obes. 2017;12:e46–50.

    Article  CAS  Google Scholar 

  29. Chang YC, Liu PH, Lee WJ, et al. Common variation in the fat mass and obesity-associated (FTO) gene confers risk of obesity and modulates BMI in the chinese population. Diabetes. 2008;57:2245–52.

    Article  CAS  Google Scholar 

  30. Cha SW, Choi SM, Kim KS, et al. Replication of genetic effects of FTO polymorphisms on BMI in a korean population. Obesity. 2008;16:2187–9.

    Article  CAS  Google Scholar 

  31. El Hajj Chehadeh S, Osman W, Nazar S, et al. Implication of genetic variants in overweight and obesity susceptibility among the young arab population of the United Arab Emirates. Gene. 2020;739:144509.

    Article  Google Scholar 

  32. Liu C, Chu C, Zhang J, et al. IRX3 is a genetic modifier for birth weight, adolescent obesity and transaminase metabolism. Pediatr Obes. 2018;13:141–8.

    Article  CAS  Google Scholar 

  33. Fu J, Li G, Li L, et al. The role of established east asian obesity-related loci on pediatric leptin levels highlights a neuronal influence on body weight regulation in chinese children and adolescents: the BCAMS study. Oncotarget. 2017;8:93593–607.

    Article  Google Scholar 

  34. Renstrom F, Payne F, Nordstrom A, et al. Replication and extension of genome-wide association study results for obesity in 4923 adults from northern Sweden. Hum Mol Genet. 2009;18:1489–96.

    Article  Google Scholar 

  35. Xi B, Shen Y, Reilly KH, et al. Sex-dependent associations of genetic variants identified by GWAS with indices of adiposity and obesity risk in a chinese children population. Clin Endocrinol (Oxf). 2013;79:523–8.

    Article  CAS  Google Scholar 

  36. Jung TWK. An introduction to latent class growth analysis and growth mixture modeling. Soc Personal Psychol Compass. 2008;2:302–17.

    Article  Google Scholar 

  37. Little TD, Jorgensen TD, Lang KM, Moore EW. On the joys of missing data. J Pediatr Psychol. 2014;39:151–62.

    Article  Google Scholar 

  38. Huberman M, Langholz B. Application of the missing-indicator method in matched case-control studies with incomplete data. Am J Epidemiol. 1999;150:1340–5.

    Article  CAS  Google Scholar 

  39. Chowdhury MZI, Turin TC. Variable selection strategies and its importance in clinical prediction modelling. Fam Med Community Health. 2020;8:e000262.

    Article  Google Scholar 

  40. Yilmaz B, Gezmen Karadag M. The current review of adolescent obesity: the role of genetic factors. J Pediatr Endocrinol Metab. 2021;34:151–62.

    Article  CAS  Google Scholar 

  41. Warrington NM, Howe LD, Paternoster L, et al. A genome-wide association study of body mass index across early life and childhood. Int J Epidemiol. 2015;44:700–12.

    Article  Google Scholar 

  42. Graff M, North KE, Richardson AS, et al. BMI loci and longitudinal BMI from adolescence to young adulthood in an ethnically diverse cohort. Int J Obes (Lond). 2017;41:759–68.

    Article  CAS  Google Scholar 

  43. The NS, Suchindran C, North KE, Popkin BM, Gordon-Larsen P. Association of adolescent obesity with risk of severe obesity in adulthood. JAMA. 2010;304:2042–7.

    Article  CAS  Google Scholar 

  44. Group HS, Foster GD, Linder B, et al. A school-based intervention for diabetes risk reduction. N Engl J Med. 2010;363:443–53.

    Article  Google Scholar 

  45. Sheen YJ, Hsu CC, Jiang YD, Huang CN, Liu JS, Sheu WH. Trends in prevalence and incidence of diabetes mellitus from 2005 to 2014 in Taiwan. J Formos Med Assoc. 2019;118(Suppl 2):66–73.

    Article  Google Scholar 

Download references


The authors thank all the participants and parents of the current study.


This study was funded to Chen YC from the Ministry of Science and Technology, Taiwan (Grant Numbers MOST 107-2314-B-038-113-MY3 and MOST 110-2628-B-038-014 and MOST 111-2628-B-038-022). This organization had no role in the study design, data analysis, or writing of this article.

Author information

Authors and Affiliations



CYC carried out cohort recruitment and experiments; WYF and CKL conceived the design of study; WYF performed the data analysis; WYF and CYC interpreted the findings. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yang-Ching Chen.

Ethics declarations

Ethics approval and consent to participate

Taiwan Puberty Longitudinal Study study was approved by the Institutional Review Board of Taipei Medical University (N201802018), that of CGH (CGH-P108107), and that of NCKUH (B-BR-108-076) and complied with the principles outlined in the Helsinki Declaration. All participants offered written consent.

Consent for publication

All authors have agreed to publish this manuscript.

Competing interests

All authors have reported that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, YF., Chien, KL. & Chen, YC. Association between genetic risk score and tri-ponderal mass index growth trajectories among different dietary consumption adolescents in a prospective Taiwanese cohort. Nutr Metab (Lond) 19, 83 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: