The prevalence of chronic-degenerative diseases is related to the presence of obesity, excess energy, fat and salt intake and reduced physical activity, which are elements of the so-called nutritional transition(Reference Yajnik1). However, it has been demonstrated that diseases that are typical of adulthood may originate during intra-uterine life through inadequate fetal development(Reference De Blasio, Gatford and McMillen2–Reference Mi, Cheng and Zhao9). According to the gestational phase during which intra-uterine growth restriction occurs, different systems will be affected, predisposing to diabetes mellitus type 2, arterial hypertension, obesity and dyslipidaemia(Reference Barker10–Reference Yajnik13).
The theory known as the ‘Fetal origins of disease’ or the ‘Thrifty phenotype hypothesis’ suggests that individuals who are exposed to nutritional deficiency in the uterus and during the first years of life adapt themselves to this deficiency. In relation to the pancreas, the adaptation mechanism consists of decreases in the number and activity of β-cells and harm to insulin-dependent tissues, especially muscles. Therefore, to survive during periods of nutritional deprivation, fetuses change their metabolic-endocrine routes and become insulin resistant. If, over the course of life, individuals with lower weights at birth are subjected to situations of increased food intake and decreased physical activity, this adaptation may become evident, thereby predisposing towards reduced insulin resistance, obesity and diabetes mellitus type 2(Reference Eriksson, Forsén and Tuomilehto6, Reference Yudkin11, Reference Yajnik and Deshmukh14, Reference Barker, Winter and Osmmond15).
In Brazil, as in other developing countries like India, the prevalence of chronic-degenerative diseases including metabolic diseases is increasing, representing an important cause of death even in young adults(Reference Ramachandran, Snehalatha and Latha16, 17).
The present study investigated the relationship between birth weight and insulin resistance in a cohort of 506 Brazilian children aged 5–8 years, involving an age group that has been assessed in few similar studies(Reference Yajnik, Fall and Vaydia18–Reference Joglekar, Fall and Deshpande24).
Materials and methods
This is an epidemiological cohort study initially involving 865 children who were followed up from birth until 5–8 years of age(Reference Rondó, Lemos and Pereira25). The children were born to women who took part in a cohort study in Jundiaí city, Brazil, between 1997 and 2000. The women attended antenatal care in twelve health units and five hospitals in Jundiaí city, and were followed before the 16th week of pregnancy to the birth of their babies. All participants were insured by the National Health Service (SUS) that assists low-income families. Women with chronic infectious diseases, metabolic diseases, cardiopathy, mental diseases, hypertension/pre-eclampsia/eclampsia, vaginal bleeding and multiple deliveries were not included in the study. The children were classified after birth as being small for gestational age (SGA), adequate for gestational age and large for gestational age(Reference Williams, Creasy and Cunningham26). Details of the study are given elsewhere(Reference Rondó, Ferreira and Nogueira27).
Out of the 865 mother–infant pairs from the previous cohort(Reference Rondó, Ferreira and Nogueira27), 745 were located and invited to participate in the present study, resulting in a final sample of 649 children, whose parents signed a free and informed consent statement and answered a general questionnaire.
The study was carried out between November 2004 and September 2006 and consisted of two phases. (1) Taking the information from the questionnaire of the first cohort study into account, the mothers who at that time were living in the city of Jundiaí and nearby municipalities were located. The mothers were invited to participate in the present study through telephone contact or by visiting their homes if they did not have a telephone. Next, a home visit was made, during which the objectives of the study were explained, and the ethical consent statement was signed by the child's parents or guardians. At the time of this home visit, a general questionnaire was applied in order to evaluate socio-economic and demographic data and any morbidity among the children. (2) The participants were contacted again by telephone to arrange to collect anthropometric data (weight, height and waist circumference) and blood samples in order to determine glucose, total cholesterol, HDL cholesterol (HDL-c), LDL cholesterol (LDL-c) and TAG concentrations.
The anthropometric data were obtained in accordance with the recommendations of Cameron(Reference Cameron28) and Jelliffe & Jelliffe(Reference Jelliffe and Jelliffe29). The children were weighed, after been fasting for 10–12 h, by a portable electronic scale (Sohnle®, model 7500, Murrhardt, Germany), with accuracy of 100 g. Their height was measured using a SECA® stadiometer (Leicester Portable Height measure model, Hamburg, Germany), with accuracy of 0·1 cm. The waist circumference measurements were obtained using a tape measure (Stanley®, model 34 103, New Britain, CT, USA) of accuracy 0·1 cm. The BMI was classified according to the Centers for Disease Control and Prevention(30) and International Obesity Task Force classifications(Reference Cole, Bellizzi and Flegal31). Mean values of waist circumference were compared with a British standard(Reference McCarthy, Jarrett and Crawley32).
Immediately after collection of the blood samples for determination of the glucose and insulin concentrations and the lipid profile, the material was centrifuged at 1500 rpm for 10 min (Tomy centrifuge, model IC-15AN, Tominaga Works Ltd, Tokyo, Japan) and then stored in a refrigerator at 0°C for a period of 24–48 h, for the biochemical tests' determinations.
Glucose was obtained by an enzymatic method using the Bayer® ADVIA 1200 clinical chemistry system (Pittsburgh, PA, USA) for glucose hexokinase. Fasting glycaemia concentration was based on the criteria currently adopted for diagnosing diabetes mellitus type 2 and was the same for adults and children(33). Insulin was assayed in serum by a chemiluminescence method using the Immulite 2000 apparatus (Immulite®, DPC, Los Angeles, CA, USA). To determine insulin resistance, the homeostasis model assessment (HOMA) method was used(Reference Matthews, Hosker and Rudenski34) by the formula: HOMA = fasting glucose (mmol/l) × fasting insulin (μU/ml)/22·5. Insulin resistance was considered to be present when HOMA was ≥ 2·5(Reference Valerio, Licenziati and Iannuzzi8, Reference D'Annunzio, Vanelli and Meschi35).
Assays for total cholesterol and HDL-c were performed by a colorimetric enzymatic method, using the Bayer® ADVIA 1200 clinical chemistry system. The concentrations of LDL-c were determined by the Friedwald(Reference Friedwald, Levy and Fredrickson36) formula. TAG concentrations were determined photometrically following an enzymatic reaction using the Bayer® ADVIA 1200 clinical chemistry system. Total cholesterol and fractions and TAG concentrations were classified in accordance with the recommendations of the First Guidelines for the Prevention of Atherosclerosis during Childhood and Adolescence(Reference Giuliano, Caramelli and Pellanda37).
The categorical variables were analysed by the χ2 test with Yates correction, applying the Fisher's exact test when the frequency was less than five. The Student's t test or the Mann–Whitney test was used to compare mean values of pairs of independent variables. The relationship between birth weight and the HOMA index, and the variables of interest were determined using Pearson's or Spearman's correlation. The impact of birth weight and other independent variables (selected according to their importance, correlation coefficient and biological collinearity) on the HOMA index was determined by multiple linear regression, using the backward stepwise selection method. The outcome – HOMA index – was analysed separately for each independent variable. The variables with descriptive level of P ≤ 0·20 (birth weight, age, sex, BMI, waist circumference, HDL-c, LDL-c, TAG and per capita income) were selected for entry into the multiple linear regression model.
In addition to the use of birth weight as a continuous variable, a second linear multiple regression model tested the association between birth weight by gestational age (in three categories: SGA, adequate for gestational age and large for gestational age) and HOMA, controlling for the variables cited above. Moreover, biological models were investigated in order to explain the causal pathway of the relationship between birth weight and insulin resistance.
In the final models, the variables with a P ≤ 0·05 were considered as being statistically significant. The software Epi Info 6.04 (Atlanta, GA, USA) and Stata version 9 (College Station, TX, USA) were used for data storage and for statistical analysis.
The present study was carried out in accordance with the Declaration of Helsinki of the World Medical Association, and it was approved by the Research Ethics Committee of the School of Public Health, University of São Paulo.
Results
The biochemical parameters were determined in the second phase of the study involving a total of 506 children. There was a loss of 32 % children, considering all the mothers who were located (n 745 mothers) in the first phase of the study, and a loss of 41·5 % children, considering the mothers who participated in the original cohort (n 865 mothers).
There was a significant difference between the mean age of the children included in the cohort (6·5 (sd 0·7) years) and the children who did not conclude the study (5·9 (sd 0·8) years; P = 0·01). However, there was no difference in relation to per capita income, birth weight or mother's schooling level. The characteristics of the children at birth are presented in Table 1. With regard to gestational age, 95·5 % were born at term ( ≥ 37 and < 42 weeks of gestation) and 53·6 % were female.
GA, gestational age.
Table 2 shows that most of the children (53·4 %) were between 6 and 7 years old. Approximately, 67·0 % of the participants had a per capita income of less than one Brazilian minimum monthly salary (one Brazilian minimum monthly salary = approximately US$ 77). According to the Centers for Disease Control and Prevention BMI classification(30), 45 (8·9 %) of the children were underweight, 346 (68·6 %) were normal, 72 (14·3 %) were at risk of becoming overweight and 41 (8·2 %) were overweight. Sixty-six children (twenty-one males and forty-five females) (12·9 %) were overweight, and twenty-seven children (twelve males and fifteen females; 5·3 %) were obese by the International Obesity Task Force classification(Reference Cole, Bellizzi and Flegal31). The percentages were lower to studies involving Mexican and British children with comparable ages(Reference Moraes, Beltran and Mondini38, Reference Jebb, Rennie and Cole39). Mean waist circumferences of the children were higher than the 75th percentile of a British standard(Reference McCarthy, Jarrett and Crawley32). High fasting glycaemia (5·6–6·9 mmol/l (700–1000 mg/l)) was observed in sixty-one children (12·0 %) and nineteen (3·7 %) presented an abnormal HOMA index ( ≥ 2·5 mmol/l μU/ml). Total cholesterol was high ( ≥ 4·4 mmol/l; 1700 mg/l)(Reference Giuliano, Caramelli and Pellanda37) in 168 children (33·2 %). One hundred and two children (20·2 %) presented HDL-c that was lower than recommended (1·2 mmol/l; 450 mg/l) and twenty-two (4·3 %) were in the group of children with abnormal LDL-c ( ≥ 3·4 mmol/l; 1300 mg/l).
BMMS, Brazilian minimum monthly salary (R$350·00 = approximately US$ 77 per capita income); HOMA, homeostasis model assessment; HDL-c, HDL cholesterol; LDL-c, LDL cholesterol.
Table 3 shows the correlations between birth weight, gestational age, demographic, socio-economic, anthropometric and biochemical factors. Birth weight was significantly correlated with all the anthropometric parameters investigated.
GA, gestational age; WC, waist circumference; TC, total cholesterol; HOMA, homeostasis model assessment; Income, per capita income.
All the variables were assessed in childhood, with the exception of BW and GA.
Table 4 presents multiple linear regression models taking the HOMA index as the dependent variable and birth weight (considered as a continuous variable and categorised by gestational age as SGA, adequate for gestational age and large for gestational age) as the independent variable. In the first and second models, a backward stepwise selection method was utilised, controlling for the following variables: age; sex; BMI; waist circumference; HDL-c; LDL-c; TAG; per capita income. In the first model, the HOMA index was negatively associated with birth weight (P < 0·001) and positively associated with waist circumference (P < 0·001). In the second model, the HOMA index was positively associated with SGA (P = 0·027) and waist circumference (P < 0·001).
SGA, small for gestational age; LGA, large for gestational age.
* First model: R 2 adjusted = 0·16; second model: R 2 adjusted = 0·12; third model: R 2 adjusted = 0·11; fourth model: R 2 adjusted = 0·12; fifth model: R 2 adjusted = 0·13; sixth model: R 2 adjusted = 0·14.
There was a positive association between birth weight and waist circumference (P < 0·001), but no interaction between the two variables (P>0·05). In a univariate linear regression analysis, birth weight did not present a statistically significant association with the HOMA index (P = 0·06; Table 3).
Tu et al. (Reference Tu, West and Ellison40) reported that an inference regarding the existence and direction of an association between birth weight and an outcome occurring later in life depends on whether one adjusts or does not adjust for current weight in setting up the regression model. Therefore, based on Tu et al. (Reference Tu, West and Ellison40) study, and trying to investigate biological models, we carried out a regression analysis assessing the association between birth weight and the HOMA index adjusting the model for body weight and also for BMI. There were statistically significant associations between birth weight, weight in childhood and the HOMA index, and birth weight, BMI in childhood and the HOMA index (third and fourth models in Table 4). The first model described in Table 4 was also adjusted for weight and BMI in childhood (fifth and sixth models in Table 4); the associations between the HOMA index, birth weight and waist circumference remained significant.
Discussion
The children assessed in this cohort presented a considerable prevalence of low ( < 2500 g), deficient or insufficient (2500–2999 g) weight (Table 1); a factor of concern with regard to their potential growth and development and risk of morbidity–mortality(Reference Wallace, Levy and Matthews41).
According to the multiple linear regression analysis utilised in the present study, birth weight (as a continuous variable and categorised as SGA) and insulin resistance, assessed by the HOMA index, were associated. A reduction of 100 g in birth weight resulted in an increase of 0·026 mmol/l μU/ml in the HOMA index. The mean value of HOMA for these children was 1·0 (sd = 0·8; Table 2). Possibly, the increase in the HOMA index can have a clinical significance for the children. According to a large review on validation of the HOMA index, the parameter has an appropriate use in cohort and epidemiological studies for assessing insulin resistance(Reference Fowden, Albrecht and Pepe42).
Yajnig et al. (Reference Yajnik, Fall and Vaydia18) reported an association between insulin resistance and birth weight in 4-year-old Indian children. Four years later, his group of researchers studied the same children (n 477) and found out that the lower birth weight was associated with increased HOMA, and that the highest levels of insulin resistance syndrome variables and total and LDL cholesterol were in children of low birth weight (LBW) but high fat mass at 8 years(Reference Bavdekar, Yajnik and Fall19). Wilkin et al. (Reference Wilkin, Voss and Metcalf20) and Jeffery et al. (Reference Jeffery, Metcalf and Hosking21) reported, in the EarlyBird study, correlations between insulin resistance and weight in approximately 300 children at 5 and 8 years of age, but none with birth weight. The Pune Maternal Nutrition Study(Reference Joglekar, Fall and Deshpande24) assessed the associations of size and body proportions at birth, and growth during infancy and childhood, to body composition and CVD risk factors (insulin resistance, blood pressure, glucose tolerance and plasma lipids) in 698 children at the age of 6 years. The authors reported that smaller mid-upper arm circumference at 6 months predicted higher insulin resistance (P < 0·001), but larger mid-upper arm circumference at 1 year predicted higher systolic blood pressure (P < 0·001). Gupta et al. (Reference Gupta, Gupta and Pareek23) reported an association between insulin resistance and birth weight in 134 Indian children and adolescents from 5 to 16 years of age, controlling for waist circumference, waist:hip ratio, age, sex and BMI. Giapros et al. (Reference Giapros, Evagelidou and Challa22) observed that large for gestational age-born children aged 6·5–8 years had higher HOMA and adiponectin concentration than adequate for gestational age children.
The possible explanations for the association between birth weight and insulin resistance are the following: (a) most of the development of the pancreatic islets occurs in the uterus and, in human subjects, the mass of β-cells increases more than 130-fold between the 12th week of gestation and the 5th month after birth(Reference Yajnik44). Studies on animals have shown that nutritional deprivation within the uterus or just after birth may reduce the number of β-cells and the secretion of insulin, while noting that protein concentrations are directly related to the production of β-cells by the pancreas(Reference Hovi, Anderson and Eriksson45). In addition to lower insulin production, it can also be considered that LBW is related to glucose intolerance(Reference Krentz46); (b) under unfavourable conditions of intra-uterine nutrition, there is a reduction in the proportion of lean mass relative to body fat. This pattern favours insulin resistance as a mean of adaptation(Reference Molist-Brunet, Jimeno-Mollet and Franch-Nadal47). These abnormalities would cause structural changes that could persist into adulthood. The main strategy for fetal adaptation to deficient nutrition is selective redistribution of nutrients for the formation of fundamental organs like the brain (brain sparing). This mechanism would result in underdevelopment of the skeleton and liver, a lower proportion of muscle mass, diminished mass of β-cells from the pancreas and abnormalities in vascular tissues(Reference Yajnik, Godbole and Lubree48, Reference Levitt, Lambert and Woods49). A study carried out in South Africa(Reference McKeigue, Shah and Marmot50) involving 132 young adults indicates that the ‘fetal origins’ expression of the chronic disease phenotype is not dependent on birth weight alone, but on its interaction with subsequent fat accumulation, including measurements of waist circumference.
The present study also shows a positive association between waist circumference and the HOMA index. Mean waist circumferences of the children involved in the study were higher than mean waist circumference measurements of British children at the same age(Reference McCarthy, Jarrett and Crawley32). However, there were not many overweight or obese children, according to the Centers for Disease Control and Prevention or International Obesity Task Force classifications. The same phenomenon has been observed in Asia: central obesity even in the absence of obesity(Reference Yajnik51). Molist-Brunet et al. (Reference Molist-Brunet, Jimeno-Mollet and Franch-Nadal47) observed that as waist circumference increased the risk that insulin resistance would occur was higher. Yajnik(Reference Yajnik51), involving 8-year-old children, also found a relationship between HOMA, insulin, waist circumference and birth weight. A study carried out on the Caribbean island of Barbados showed an association between LBW girls with a family history of diabetes and higher HOMA. LBW girls without a family history of diabetes had a smaller waist circumference than normal birth weight girls. In contrast, LBW girls with a family history of diabetes had similar mean waist circumference than normal birth weight girls with a family history of diabetes(Reference Chambers, Tull and Fraser52).
As stated by Barker(Reference Barker, Winter and Osmmond15), the explanation for the accumulation of centrally located fat is that, during intra-uterine life, energy saving may occur in order to maintain brain development and functioning. These savings would generate a positive energy balance that, in adulthood, would cause an accumulation of fat, especially in the central region of the body. This centralisation of adiposity is related to insulin resistance and lower development of muscle mass, and it consists of a substrate enabling metabolic speed in providing glucose through stress-related metabolic pathways, with the aim of maintaining adequate glucose levels, particularly in the brain(Reference Bjorntorp53).
The fetal origins hypothesis, as related to insulin resistance, affirms the existence of a causal path from prenatal factors to insulin resistance. From the public health viewpoint, we want to know whether interventions undertaken in pregnancy might be of benefit to the health of the offspring later in life.
According to the present results, it is difficult to conclude that waist circumference is a potential pathway between birth weight and insulin resistance. Waist circumference in childhood may be influenced not only by birth weight(Reference Barker, Winter and Osmmond15), but probably by prenatal factors acting through effects on fetal growth (as measured by birth weight)(Reference Li, Manor and Power54) and factors occurring in childhood. Poor prenatal nutrition, for example, may act as a common cause of lower birth weight and higher waist circumference, determining insulin resistance.
Despite the non-significant result of the association between birth weight and the HOMA index (in the univariate linear regression analysis), a P-value of 0·06 has probably a biological value, considering the number of factors that could intervene in the association throughout prenatal life to childhood. However, if we consider that this association is not statistically significant, and that the associations between birth weight and waist circumference and waist circumference and insulin resistance are significant, we can probably suspect, according to Hernán et al. (Reference Hernán, Hernández-Díaz and Werler55), that lower birth weight is a common cause of higher waist circumference and insulin resistance.
Weinberg(Reference Weinberg56) asserted that aetiological inferences cannot be trusted to a regression package and that decisions about confounding and when to adjust can have a qualitative impact on our conclusions. According to the author, while we often cannot construct a complete causal pathway, we can gain important guidance from the attempt, clarify the assumptions we are making and then focus on those epidemiologic questions that are answerable.
Epidemiological studies have been conducted over the last decades in an attempt to determine the relationship between birth weight and outcomes in later life, most of them retrospective studies involving adult populations. Few prospective studies have been carried out on children(Reference Yajnik, Fall and Vaydia18–Reference Joglekar, Fall and Deshpande24). To our knowledge, there is no Brazilian cohort study investigating the association between birth weight and metabolic abnormalities in this age group.
In conclusion, the results of the present study showed increased insulin resistance in apparently healthy, young children, who had lower weight at birth and higher measurements of waist circumference. There is a need to carry out large cohort epidemiological studies assessing several factors throughout pregnancy and childhood to a better understanding of the causal pathway between prenatal factors and insulin resistance. Some questions that now are quite difficult to answer in terms of causal pathways will probably be better understood, and help to develop public health policies that adopt preventive measures to promote adequate maternal-fetal and child development and enable early diagnosis of metabolic abnormalities.
Acknowledgements
The authors gratefully acknowledge Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP; Grant No. 04/04 109-8), for financing the present study, Prof. José Maria P. Souza, and the field workers Alessandra L. Santos, Rosemary F. Ferreira and Josimara F. Moura. There are no conflicts of interest to declare.
J. O. L. collected data, participated in the statistical analysis and interpretation of data, and writing of the paper. P. H. C. R. designed the study protocol, secured funding, facilitated data collection, participated in the statistical analysis and interpretation of data, and writing of the paper. J. A. P. collected data, participated in the statistical analysis and interpretation of data, and writing of the paper. R. G. O. participated in data collection, statistical analysis and writing of the paper. M. B. S. F. participated in the statistical analysis, interpretation of data and writing of the paper. P. B. S. participated in the statistical analysis, interpretation of data and writing of the paper.