Skip to main content
  • Research article
  • Open access
  • Published:

Reduced genetic influence on childhood obesity in small for gestational age children

Abstract

Background

Children born small-for-gestational-age (SGA) are at increased risk of developing obesity and metabolic diseases later in life, a risk which is magnified if followed by accelerated postnatal growth. We investigated whether common gene variants associated with adult obesity were associated with increased postnatal growth, as measured by BMI z-score, in children born SGA and appropriate for gestational age (AGA) in the Auckland Birthweight Collaborative.

Methods

A total of 37 candidate SNPs were genotyped on 547 European children (228 SGA and 319 AGA). Repeated measures of BMI (z-score) were used for assessing obesity status, and results were corrected for multiple testing using the false discovery rate.

Results

SGA children had a lower BMI z-score than non-SGA children at assessment age 3.5, 7 and 11 years. We confirmed 27 variants within 14 obesity risk genes to be individually associated with increasing early childhood BMI, predominantly in those born AGA.

Conclusions

Genetic risk variants are less important in influencing early childhood BMI in those born SGA than in those born AGA, suggesting that non-genetic or environmental factors may be more important in influencing childhood BMI in those born SGA.

Peer Review reports

Background

Worldwide more than 30 million infants per year are born small for gestational age (SGA) defined by birth weight below the 10th percentile, according to gestational age [1]. These infants are at increased risk of mortality and morbidity in early life as well as greater later risk of developing obesity, dyslipidemia, hypertension and diabetes.

It has been suggested that growth restriction in the womb resulting from poor maternal nutrition, can program the individual’s physiology and metabolism predisposing them to long-term susceptibility to metabolic diseases [2, 3]. Alternative hypotheses explaining the link between low birthweight and higher risk of adult metabolic diseases include the foetal insulin hypothesis [4] that suggests common genetic variants which influence insulin secretion or action may both reduce foetal growth (as insulin is a key foetal growth factor) and increase the risk of type 2 diabetes in later life.

Candidate gene and genome wide association studies have led to the discovery of over 50 gene variants known to be associated with adult obesity [5, 6]. Two recent studies involving the Avon Longitudinal Study of Parents and Children [7], and the 1946 British Birth Cohort study [8], have found that a genetic predisposition score of 8 and 11 adult obesity susceptibility alleles were associated with greater infancy and childhood weight, height and BMI gain. The obesity risk allele score showed a weaker association with BMI than with weight until 3.5 years, due to concurrent gains in both length and weight during infancy. While combining established genetic variants for BMI into a risk-allele score maximised statistical power, this approach masked possible heterogeneity in effects of individual variants.

The detrimental effect of poor fetal growth on long-term metabolic health is observed to be magnified if followed by accelerated postnatal growth [9]. The contribution of genetic factors to such accelerated postnatal growth in those born SGA is unknown. We hypothesized that certain obesogenic gene variants would have a disproportionately higher impact in increasing early childhood BMI in those born SGA than those born appropriate for gestational age (AGA), and therefore explain part of the link between low birthweight and subsequent risk of obesity. This is based on the assumption that underlying metabolic programming in those born SGA would enhance the biological metabolic effect of any Mendelian inherited gene variant. We explored the childhood growth of those born SGA and AGA and investigated the association of candidate obesity genes with the BMI z-score in the Auckland Birthweight Collaborative (ABC) study children, in relation to SGA status.

Methods

The Auckland birthweight collaborative study

The ABC Study has been described in more detail elsewhere [10]. The study principally investigates factors that relate to childhood growth and development. All children included in this study were born at 37 weeks of gestation or more. Infants were defined as SGA if their birthweight was equal to or below the sex-specific 10th percentile for gestational age in the New Zealand population. Infants were defined as AGA if their birthweight was greater than the sex-specific 10th percentile for gestational age in the New Zealand population.

Subjects

Eligible for inclusion in the ABC study were all SGA infants and a selection of one in sixteen AGA infants, born in the Waitemata Health or Auckland Healthcare regions in New Zealand between 16 October 1995 and 12 August 1996. All SGA children and a selection of one in eight AGA children, who were born between 12 August 1996 and 30 November 1997 in the Auckland Healthcare region, were also eligible for inclusion. Data have been collected at birth (phase 1), 12 months (phase 2), 3.5 years (phase 3), 7 years (phase 4), and 11 years (phase 5). This study utilises phases 3, 4, and 5 data. A total of 547 children (228 SGA and 319 AGA, phase 5) have completed BMI and 542 children have genotype data. The sample was restricted to children who were born to European mothers because of the poor response rate from non-European ethnic groups at earlier phases. DNA was extracted from the blood/buccal samples using Qiagen’s DNA extraction kit and following the manufacturer’s instructions. The study received ethical approval from the Northern Regional Ethics Committee. The 3.5 year old children’s study had the ethic approval number of 99/097; the study of 7 year old children is AKX/02/00/319; and the study of 11 year old children is NYT/06/00/112. Signed consent for the study was given by the parents of the children at all phases and also given by the child at age 11 years.

Outcome measure

Body Mass Index (weight/height2) is a common method of assessing obesity status that was used in this study. BMI scores were converted to z-scores using the LMS method to ensure normal distribution of the data [11].

SNP selection

The SNPs were selected from a systematic literature search to identify genetic variants demonstrating association with obesity using the criteria previously described [12]. We included 37 obesity SNPs identified from published candidate gene or genome-wide association studies, located in 21 genes. Genotypic data were not available for nine SNPs: rs1387153 (MTNR1B), rs6020339 (CTNNBL1), rs2388399 (PFKP), rs477181 (MC4R), rs3820152 (ADIPOR1), rs7561317 (TMEM18), rs7647305 (ETV5 – DGKG), rs2844479 (NCR3 - AIF1), and rs10487818 (NAMPT). These SNPs either could not be multiplexed into our sequenom assays, failed or did not pass our quality control measure for inclusion in the analysis. Table 1 provides the list of genes/SNPs selected for investigation including the tested allele and its frequency.

Table 1 Obesity SNPs for association with BMI in the Auckland Birthweight Collaborative study samples (Morgan et al.[12])

Genotyping

Genotyping was performed with the MassARRAY and iPlex systems of the Sequenom genotyping platform (Sequenom, San Diego, CA), which uses the MALDI-TOF primer extension assay [13, 14], according to manufacturer’s recommendations. Assays were optimized in 24 samples consisting of 20 reference Centre d’Etude du Polymorphisme Humain (CEPH) samples and 4 blanks. All sample plates contained cases, controls, blanks, CEPH and duplicate samples. Quality control measures included independent double genotyping, blind to sample identity and blind to the other caller, and where available, comparison of our CEPH genotypes to those in the HapMap (http://www.hapmap.org).

Statistical analysis

Generalised linear model was used to test the linearity of genotype-phenotype relationship for quantitative traits. Since the outcomes of interest, BMI z-score, was a repeated measure, the GENMOD procedure with the REPEATED statement in SAS was used to fit a generalised estimating equations (GEE) model. For linearity of the genotype-phenotype relationship for quantitative traits, additive model was used, each SNP was coded 0, 1, and 2 for each tested allele [15]. The association with SNPs was carried out in four ways due to the disproportionate sampling of AGA and SGA at birth (phase 1): using stratified by (1) AGA and (2) SGA children, (3) using a main effect model with adjustment for the SGA status, and (4) interaction with SGA status. False discovery rate (FDR) q-values were calculated for the total number of tests (148 tests; 37 SNPs times four analyses) using the ‘qvalue’ package in R [16, 17] with the null hypothesis fraction estimated by the smoothing method of Storey and Tibshirani [18]. A prespecified FDR threshold was set at 0.05. Table 2 provides the results of the analysis of association between SNPs and BMI z-score with both q-values and unadjusted p-values; a SNP is called significant if the q-value is less than the prespecified 0.05. SAS (V9.2 SAS Institute., Cary, NC, USA) and R (R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, http://www.R-project.org) were used for statistical analyses. This resulted in a p-value of 0.05 as being statistically significant.

Table 2 Analysis of association between SNPs and repeated measures of BMI z-score

Results

Growth patterns of ABC participants according to SGA status

SGA children had significantly lower BMI z-score than AGA children at all three assessments conducted at 3.5, 7 and 11 years of age (Table 3 and Figure 1). BMI z-score showed no significant interaction between SGA status and age at assessment (p = 0.254) (Table 4 and Figure 2).

Table 3 Differences in birth weight and BMI z-score
Figure 1
figure 1

Birth weight by SGA status.

Table 4 BMI z-score Growth by SGA status
Figure 2
figure 2

BMI z-score by SGA status.

Genotype data

A total of 37 candidate SNPs for obesity were successfully genotyped (Table 1), and their relationship with BMI z-score was tested. There was no indication of a deviation from Hardy Weinberg equilibrium (HWE), except for one SNP (rs7903146) which was excluded from further analysis. Using a prespecified FDR threshold at α = 0.05 [17], we determined that 27 SNPs were significantly associated with the BMI z-score in at least one of the four analyses: stratified by AGA and SGA children, after adjustment for SGA status, and interaction with SGA status (Table 2).

Variants associated with BMI z-score in AGA children

A total of 16 variants: rs2286385 (ADIPOR2), rs6013029 and rs6020712 (CTNNBL1), rs7138803 (FAIM2), rs8050136 and rs9939609 (FTO), rs1042725 (HMGA2), rs12970134, rs17700633, rs17782313, rs4450508, and rs502933 (MC4R), rs2293855 (MTMR9), rs10508503 (PTER), rs10913469 (SEC16B), and rs7498665 (SH2B1) in 10 genes were found to be significantly associated with BMI z-score in AGA children after correction for multiple testing (Table 2). The impact on BMI z-score increase varied with the tested variant allele from 0.100 (rs7498665 SH2B1) to 0.309 (rs6020712 CTNNBL1).

Variants associated with BMI z-score in SGA children

Ten variants in seven genes were significantly associated with BMI z-score in SGA children after correction for multiple testing (Table 2). The impact on BMI z-score increase varied from 0.133 (rs2815752 NEGR1) to 0.458 (rs6020395 CTNNBL1). In contrast to those born AGA, none of the 5 variants in MC4R were associated with increasing BMI z-score in those born SGA.

Variants associated with BMI z-score after adjustment for SGA status

A total of 20 variants: rs2286385 (ADIPOR2), rs6265 (BDNF), rs16986921, rs6013029, rs6020395, rs6020712, rs6020846, rs6096781, and rs6125962 (CTNNBL1), rs6499640, rs8050136, and rs9939609 (FTO), rs1042725 (HMGA2), rs11084753 (KCTD15), rs12970134, rs17782313, rs4450508, and rs502933 (MC4R), rs10508503 (PTER), and rs10913469 (SEC16B) showed significant association after adjustment SGA status and also correction for multiple testing. BMI z-score increased with all tested obesity risk variant alleles ranging from 0.087 (rs1042725 HMGA2) to 0.303 (rs6013029 CTNNBL1).

Variants associated with interaction between BMI z-score and SGA status

A total of eight variants: rs1424233 (MAF), rs12970134 and rs17782313, (MC4R), rs2293855 (MTMR9), rs2568958 and rs2815752 (NEGR1), rs10508503 (PTER), and rs7498665 (SH2B1) showed significant interaction with SGA status. Two of these variant (MC4R and PTER) also showed a significant association with BMI z-score among the whole cohort with adjustment for SGA status, however, their effect was only observed among those born AGA, with no effect in SGA. Of the remaining obesity risk variants showing a significant interaction with SGA status, two variants in NEGR1 were associated with increased BMI in those born SGA but no effect in those born AGA, while SH2B1 obesity risk variant showed increased BMI in those born SGA and reduced BMI in those born AGA. There was a significant interaction between the MAF variant and SGA status, and reduced BMI in the SGA group (Table 2). However, there were no significant effects on BMI when the association was tested within those born SGA and AGA.

Discussion

We confirmed 27 variants within 14 obesity genes to be individually associated with increasing early childhood BMI in those born AGA, most of which were also associated with increasing early childhood BMI in the whole cohort after adjusting for SGA status. Variants associated with the interaction of BMI z-score and SGA status were largely driven by their predominant association with increasing BMI in either the AGA group (MC4R, PTER, MTMR9), or the SGA group (NEGR1, SH2B1), rather than a consistently larger effect seen in the SGA group as hypothesised. None of the 5 MC4R obesity risk variants associated with increasing BMI in AGA were significant in those born SGA. However, all three of the FTO obesity risk variants showed an effect of increasing BMI in those born SGA, one of which was also significant in those born AGA.

Although the biological functions of the majority of the obesity gene variants are still being investigated, most of the obesity risk variants are thought to act centrally on increasing hyperphagic drive [19]. The other physiological mechanisms by which these gene variants could be acting to influence postnatal growth would be through a change in control of body’s energy partitioning from lean to fat tissue, or shift in metabolic efficiency through thermogenesis.

Our New Zealand SGA birth cohort did not show evidence of accelerated postnatal growth, in terms of their BMI z-score, which remained lower than those born AGA at assessment age 3.5, 7 and 11 years, however, they may still have increased postnatal fat deposition relative to lean body or skeletal mass. There is evidence from studies conducted in Switzerland [20] and Spain [21] that pre-pubertal children born SGA have more body fat and less lean tissue than weight or age matched controls. Nonetheless, the higher cardiovascular risks observed with accelerated postnatal growth in those born SGA have been demonstrated to occur with a rapid tempo of childhood gain in BMI rather than body fat [22].

Our results do not suggest that the known obesity risk gene variants have a disproportionately higher impact on increasing childhood BMI in those born SGA compared to AGA. However, until the age of 11 years, those born SGA remained smaller than AGA children, and it remains to be seen whether they develop a greater prevalence of obesity in adulthood. This work leads us to believe that non-genetic or environmental factors are more likely to be important in influencing childhood BMI in those born SGA.

Conclusions

In summary, our longitudinal data with measurement of BMI at 3.5, 7 and 11 years of age, confirm that several of the adult obesity risk gene variants (ADIPOR2, BCDIN3D_FAIM2, BDNF, CTNNBL1, FTO, HMGA2 KCTD15, MTMR9 MC4R, PTER, SEC16B, and SH2B1) were also associated with childhood growth in New Zealand. Most of the association between obesity variants with childhood BMI was observed within children born AGA rather than the SGA, suggesting that non-genetic or environmental factors may be more important in influencing childhood BMI in those born SGA.

Abbreviations

ABC:

Auckland Birthweight Collaborative

AGA:

Appropriate-for-gestational-age

SGA:

Small-for-gestational-age

SNPs:

Single nucleotide polymorphisms

CEPH:

Centre d’Etude du Polymorphisme Humain

HWE:

Hardy Weinberg Equilibrium

MALDI-TOF:

Matrix Assisted Laser Desorption/Ionization – Time of Flight.

References

  1. Lee PA, Chernausek SD, Hokken-Koelega AC, Czernichow P, International Small for Gestational Age Advisory B: International Small for Gestational Age Advisory Board consensus development conference statement: management of short children born small for gestational age, April 24-October 1, 2001. Pediatrics. 2003, 111 (6 Pt 1): 1253-1261.

    Article  PubMed  Google Scholar 

  2. Barker DJ: Maternal nutrition, fetal nutrition, and disease in later life. Nutrition. 1997, 13 (9): 807-813.

    Article  CAS  PubMed  Google Scholar 

  3. Bruce KD, Hanson MA: The developmental origins, mechanisms, and implications of metabolic syndrome. J Nutr. 2010, 140 (3): 648-652.

    Article  CAS  PubMed  Google Scholar 

  4. Hattersley AT, Tooke JE: The fetal insulin hypothesis: an alternative explanation of the association of low birthweight with diabetes and vascular disease. Lancet. 1999, 353 (9166): 1789-1792.

    Article  CAS  PubMed  Google Scholar 

  5. Frayling TM, Timpson NJ, Weedon MN, Zeggini E, Freathy RM, Lindgren CM, Perry JR, Elliott KS, Lango H, Rayner NW, et al: A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science. 2007, 316 (5826): 889-894.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Silventoinen K, Pietilainen KH, Tynelius P, Sorensen TI, Kaprio J, Rasmussen F: Genetic and environmental factors in relative weight from birth to age 18: the Swedish young male twins study. Int J Obes. 2007, 31 (4): 615-621.

    Article  CAS  Google Scholar 

  7. Elks CE, Loos RJ, Sharp SJ, Langenberg C, Ring SM, Timpson NJ, Ness AR, Davey Smith G, Dunger DB, Wareham NJ, et al: Genetic markers of adult obesity risk are associated with greater early infancy weight gain and growth. PLoS Med. 2010, 7 (5): e1000284-

    Article  PubMed  PubMed Central  Google Scholar 

  8. Elks CE, Loos RJ, Hardy R, Wills AK, Wong A, Wareham NJ, Kuh D, Ong KK: Adult obesity susceptibility variants are associated with greater childhood weight gain and a faster tempo of growth: the 1946 British Birth Cohort Study. Am J Clin Nutr. 2012, 95 (5): 1150-1156.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Ong KK, Loos RJ: Rapid infancy weight gain and subsequent obesity: systematic reviews and hopeful suggestions. Acta Paediatr. 2006, 95 (8): 904-908.

    Article  PubMed  Google Scholar 

  10. Thompson JM, Clark PM, Robinson E, Becroft DM, Pattison NS, Glavish N, Pryor JE, Wild CJ, Rees K, Mitchell EA: Risk factors for small-for-gestational-age babies: The Auckland Birthweight Collaborative Study. J Paediatr Child Health. 2001, 37 (4): 369-375.

    Article  CAS  PubMed  Google Scholar 

  11. Cole TJ: The LMS method for constructing normalized growth standards. Eur J Clin Nutr. 1990, 44 (1): 45-60.

    CAS  PubMed  Google Scholar 

  12. Morgan AR, Thompson JM, Murphy R, Black PN, Lam WJ, Ferguson LR, Mitchell EA: Obesity and diabetes genes are associated with being born small for gestational age: results from the Auckland Birthweight Collaborative study. BMC Med Genet. 2010, 11: 125-

    Article  PubMed  PubMed Central  Google Scholar 

  13. Jurinke C, van den Boom D, Cantor CR, Koster H: The use of MassARRAY technology for high throughput genotyping. Adv Biochem Eng Biotechnol. 2002, 77: 57-74.

    CAS  PubMed  Google Scholar 

  14. Storm N, Darnhofer-Patel B, van den Boom D, Rodi CP: MALDI-TOF mass spectrometry-based SNP genotyping. Methods Mol Biol. 2003, 212: 241-262.

    CAS  PubMed  Google Scholar 

  15. Balding DJ: A tutorial on statistical methods for population association studies. Nat Rev Genet. 2006, 7 (10): 781-791.

    Article  CAS  PubMed  Google Scholar 

  16. Dabney A, Storey JD, Warnes GR: qvalue: Q-value estimation for false discovery rate control. R package version 1.26.0. http://www.bioconductor.org/packages/release/bioc/html/qvalue.html,

  17. Storey JD: A direct approch to false discovery rates. J R Stat Soc Ser B (Stat Methodol). 2002, 64 (3): 479-498.

    Article  Google Scholar 

  18. Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100 (16): 9440-9445.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Choquet H, Meyre D: Genetics of Obesity: What have we Learned?. Current genomics. 2011, 12 (3): 169-179.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Jornayvaz FR, Selz R, Tappy L, Theintz GE: Metabolism of oral glucose in children born small for gestational age: evidence for an impaired whole body glucose oxidation. Metabolism. 2004, 53 (7): 847-851.

    Article  CAS  PubMed  Google Scholar 

  21. Ibanez L, Ong K, Dunger DB, de Zegher F: Early development of adiposity and insulin resistance after catch-up weight gain in small-for-gestational-age children. J Clin Endocrinol Metab. 2006, 91 (6): 2153-2158.

    Article  CAS  PubMed  Google Scholar 

  22. Krishnaveni GV, Veena SR, Wills AK, Hill JC, Karat SC, Fall CH: Adiposity, insulin resistance and cardiovascular risk factors in 9-10-year-old Indian children: relationships with birth size and postnatal growth. J Dev Orig health Dis. 2010, 1 (6): 403-411.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references

Acknowledgements

The initial study was funded by the Health Research Council of New Zealand. The 12 month postal questionnaire was funded by Hawkes Bay Medical Research Foundation. The 3.5 year follow-up study was funded by Child Health Research Foundation, Becroft Foundation and Auckland Medical Research Foundation. The 7 year follow-up study was funded by Child Health Research Foundation. The 11 year follow-up was funded by Child Health Research Foundation and the National Heart Foundation. The genetic component of this study was funded by Child Health Research Foundation. We thank Professors T. Lumley, C.J. Wild and C.M. Triggs from the Department of Statistics, University of Auckland, for statistical consultation. EA Mitchell and JMD Thompson are supported by CureKids. The 7 year follow-up study was conducted in the Children’s Research Centre which is supported in part by the Starship Foundation and the Auckland District Health Board. We acknowledge the assistance of Gail Gillies, Barbara Rotherham, and Helen Nagels for contacting or assessing the participants. We sincerely thank the parents and children for participating in this study. DY Han, AR Morgan, WJ Lam, and LR Ferguson are supported by Nutrigenomics New Zealand which is a collaboration between AgResearch Ltd., Plant & Food Research and The University of Auckland with funding through the Foundation for Research Science and Technology.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Dug Yeo Han, Edwin A Mitchell or Lynnette R Ferguson.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

DYH performed the data analyses. DYH and RM wrote the manuscript. ARM and WJL performed genotyping and contributed to the writing of the manuscript. LRF, CRW, and KEW contributed to the writing of the manuscript. JMDT had responsibility for outcome assessment and contributed to the writing of the manuscript. EAM supervised the design and execution of the study and contributed to the writing of the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Han, D.Y., Murphy, R., Morgan, A.R. et al. Reduced genetic influence on childhood obesity in small for gestational age children. BMC Med Genet 14, 10 (2013). https://doi.org/10.1186/1471-2350-14-10

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2350-14-10

Keywords