Research on developmental origins of health and disease has pointed to the importance of body size and composition at birth, which reflect the adaptation of the fetus to its intra-uterine environment, and is linked to the risk of subsequent metabolic diseases( Reference Forsen, Eriksson and Tuomilehto 1 ). Neonatal body composition has been a focus for many research studies as it may be a modifiable factor for health and diseases later in life( Reference Ibanez, Ong and Dunger 2 ). Body composition can be measured by a variety of methods, which differ in accuracy, feasibility, cost and complexity. However, some methods are not suitable for neonates, including under-water weighing, stable isotope dilution and dual-energy X-ray absorptiometry (DXA), whereas others methods (such as MRI) are limited by their cost, especially for large epidemiological studies.
The measurement of body composition in neonates or infants is challenging, as the distribution of fluid in the body changes rapidly during the first few days of life, with 5–10 % reduction in body weight attributable to the change in total body water (TBW) and body composition( Reference Podratz, Broughton and Gustafson 3 , Reference Lorenz 4 ). In addition, highly accurate techniques (e.g. MRI, stable isotope dilution) in neonates and infants generally require high compliance from the subject. This can be difficult to achieve in neonates and infants. Air-displacement plethysmography (ADP) for healthy infants has been validated in several studies using the 2H dilution method for body water and a four-compartment body composition model as the reference methods( Reference Urlando, Dempster and Aitkens 5 – Reference Ellis, Yao and Shypailo 7 ). ADP does not require the infant to be restrained during the 2-min measurement. However, the infant ADP device is expensive, large and not readily portable and generally not available outside of major hospitals and research facilities. Therefore, an inexpensive, portable and reliable method to measure body composition is clearly required at this age.
Bioelectrical impedance analysis (BIA) is a simple, non-invasive and useful method for estimating body composition in epidemiological and clinical research( Reference Heitmann 8 ). BIA measures impedance of the body to an imperceptible, harmless electric current transmitted through electrodes placed on the hands and the feet. Impedance (Z) of a conductor is proportional to L2/Z. Thus, impedance of the whole body can be used to estimate the volume of TBW in the aqueous tissue compartments. Fat-free mass (FFM) can be calculated from TBW by using a hydration fraction for FFM, usually assumed to be 0·732( Reference Sheng and Huggins 9 ). Fat mass (FM) can be estimated by subtracting FFM from body weight.
As BIA was first introduced as a method for estimating body composition, many studies have demonstrated its validity to measure body water in healthy populations of older children and adults with normal fluid distribution( Reference Heitmann 8 , Reference Kyle, Bosaeus and De Lorenzo 10 – Reference Kriemler, Puder and Zahner 17 ). However, research on its usefulness for predicting body composition in neonates is limited, especially in Asians. Although a few small studies have reported the use of BIA in infants, most of them have not compared the performance of equations developed from BIA neither with that of simple anthropometry nor did they use an independent validation group( Reference Dung, Fusch and Armbrust 18 – Reference Tang, Ridout and Modi 21 ). In addition, the studies were on infants requiring intensive care or on low birth weight infants in whom the body composition profile may be different from that of healthy infants. Lingwood et al.( Reference Lingwood, Storm van Leeuwen and Carberry 22 ) developed prediction equations for FFM based on BIA in a cohort of seventy-seven healthy infants in the first few months of life, and reported that the contribution of the BIA parameter in their model was not statistically significant and did not improve the prediction over that using weight alone until their study infants reached 3 months of age. In a recent study on predicting FFM using BIA reference to that of ADP also suggested the potential use of BIA in infants, and the subjects used to develop the prediction equation were of mixed age range – that is, from 0 to 6 months( Reference Wibaek, Kaestel and Skov 23 ). These studies on predicting FFM in infants using impedance from BIA were conducted in the Western population, and thus may not necessarily be applicable to the Asian population. A prediction equation for FFM using BIA in Asian infants would be a valuable tool, especially in large cohorts or longitudinal cohort studies.
Our aims were to develop a prediction equation for FFM during the early neonatal period based on BIA and to validate that prediction using the ADP device, PEA POD infant body composition system, in Asian neonates. The applicability of this equation will also be tested in the independent validation group in this cohort as well as in an independent Western cohort of children from a study conducted at the University of Queensland.
Methods
Study design and setting
The neonates were participants of Growing Up in Singapore Towards Healthy Outcomes (GUSTO), a prospective birth cohort study in Singapore( Reference Soh, Tint and Gluckman 24 ). Pregnant women attending maternity units of two major public hospitals in Singapore, the KK Women’s and Children’s Hospital (KKH) and National University Hospital (NUH), were recruited in their first trimester between June 2009 and September 2010. All the participants were Singapore citizens or permanent residents, intended to deliver either in KKH or NUH and to reside in Singapore for the next 5 years. Both participant and her partner were of Chinese, Malay or Indian ethnicity with homogeneous parental ethnic background. Participants receiving chemotherapy or psychotropic drugs and those having type I diabetes mellitus were excluded. Recruited women were followed-up throughout pregnancy and delivery, and their children were followed-up postnatally.
Subjects
Participating neonates were born between November 2009 and May 2011 to GUSTO mothers. A total of 173 healthy singleton neonates who completed the BIA and ADP measurements at birth (0–3 d), and 140 who completed them at week 2 (5–17 d), were included in this analysis. In total, fifty-two neonates were tested at both ages. Most neonates were born at term, although four neonates in each group were born at 35–37 completed weeks of gestation.
Ethics
Study ethics approvals were obtained from the Centralized Institutional Review Board, SingHealth and Institutional Review Board of the National Health Care Group of Singapore.
Anthropometry
All anthropometric measurements at birth were carried out within 24 h after delivery (median: 1 d) and at week 2 (median: 10 d), respectively. The weight (W) of the neonate was measured using an integrated scale in the PEA POD infant body composition system. Recumbent length (L) was measured to the nearest 0·1 cm using a SECA infant mat (SECA 210 Mobile Measuring Mat; SECA Corp.). Recumbent length was measured in duplicate for reliability, and a third measurement was performed if the difference between the first two measurements was >1 cm. The average of the 2 or 2 closest of the 3 measurements was used for the analysis.
Body composition measurement by PEA POD Infant body composition system
The PEA POD Infant Body Composition System version 3.1.0 (Cosmed) was used to measure body composition – that is, FM and FFM. The system was calibrated each day before measuring the neonates. A cylinder with known volume was used to calibrate the chamber, and a 2-kg weight was used to calibrate the scale. Clothing was removed and a tight-fitting cap was placed on the head to minimise the amount of air trapped in the hair. Items on the neonate that could not be removed – for example, the umbilical clamp or the hospital identification bracelet – were used to tare the scale for body weight and volume measurements. The neonate was then placed on the scale to measure body mass, and then placed inside the chamber for body volume measurement, which required approximately 2 min. Weight and volume measurements by the PEA POD system are very precise with CV for repeated measurements of 0 and 0·02–0·09 %, respectively( Reference Urlando, Dempster and Aitkens 5 ). Infant percent body fat was computed by software integral to the PEA POD system. Body mass and body volume of the neonate were used to calculate body density (D B). The raw D B was adjusted for thoracic volume and surface area artefact. The values for the densities of FM and FFM were used to calculate % fat and % FFM( Reference Fomon, Haschke and Ziegler 25 , Reference Butte, Hopkinson and Wong 26 ). The density of fat, 0·9007 kg/l, remains constant throughout life. This value is used by the PEA POD for fat density estimation( Reference Butte, Hopkinson and Wong 26 ). Contrary to fat, the density of FFM changes from birth onwards. Thus, age- and sex-specific densities of FFM based on multi-compartment studies are used by the PEA POD( Reference Fomon, Haschke and Ziegler 25 , Reference Butte, Hopkinson and Wong 26 ).
ADP is suggested to be more accurate for body fat measurement than DXA( Reference Fields and Goran 27 ). There is good agreement for within-subject measurements using the PEA POD with a mean differences of <1 %( Reference Ellis, Yao and Shypailo 7 , Reference Roggero, Gianni and Amato 28 ). Likewise, there is good agreement between measurements of % fat in term and preterm neonates using ADP and measurements using labelled water and a four-compartment model. A study on forty-nine healthy infants demonstrated no significant difference in the mean % fat measured by ADP and that by the four-compartment model (r 2 0·73, standard error of the estimate (see)=3·7 % fat)( Reference Ellis, Yao and Shypailo 7 ). A similar study compared % fat by ADP and H2 18O dilution in seventy preterm infants and nine term infants and also reported good agreement (r 2 0·63, see=1·65 % fat)( Reference Roggero, Gianni and Amato 28 ). The regression between the two different % fat measurements did not deviate significantly from the line of identity for both studies.
Bioelectrical impedance analysis
Bioimpedance was measured using the ImpediMed SFB7 (ImpediMed), a single-channel, tetra-polar bioimpedance spectroscopy device and its electrode 292-STE single-tab electrodes. The device measures resistance and reactance at 256 frequencies between 3 and 1000 kHz. We chose to use resistance at 50 kHz as this frequency is used in single-frequency analysers. Therefore, the prediction equation developed will be applicable to predict FFM in infants using the more common single-frequency analysers.
Neonates were placed in the supine position. Insulating materials (e.g. thin cotton blanket) were used to prevent skin-to-skin contact among the body parts. Sense electrodes were placed on the dorsum of the wrist and dorsum of the ankle at the level of the styloid process and the medial malleolus. Source electrodes were placed on the palm and sole over the metacarpals and metatarsals. Good connection between electrodes and the lead as well as skin contact of the electrodes were ensured. Ten consecutive BIA measurements were completed within 2 min.
Statistical analysis
A random number generator was used to divide the subjects at each study visit (birth and week 2) into two groups: a model development group and a model validation group. The model development group included approximately two-thirds of the subjects and was used to develop prediction equations for FFM (n 116 and 96 for birth and week 2, respectively). The model validation group comprised approximately one-third of the subjects (n 57 and 46 at birth and week 2, respectively). FFM prediction equations based on BIA or simple anthropometry developed from ‘model development groups’ were applied to model validation groups. Subsequently, these predicted FFM were compared with FFM measured by ADP (FFMADP) at respective time points.
For prediction of FFM, either simple anthropometric measurements alone or in combination with the impedance quotient, L2/R50 (cm2/Ω), were used as predictors; R50 is the resistance at 50 kHz and was used for all analyses. In the model development group, the prediction equation for FFM was developed by multiple linear regression analysis using FFMADP as the dependent variable; predictors were W (kg), sex (S); (male=1, female=2) together with L (cm) alone or L2/R50. Weight, length and sex were used in the prediction as these are significant predictors at least at one of the 2 time points (for length and sex). In addition, weight and length reflect the body size of the neonate. Sex was included as it is well accepted that there is sex difference in body composition; R50 is the resistance at 50 kHz and was used for all analyses, and 50 kHz is the frequency most commonly used in single-frequency impedance analysers.
Standardised regression coefficients were used to assess the independent contribution of each variable to the dependent variable (FFM). The developed prediction equation was then used to derive FFM in the validation group within the GUSTO cohort as an independent assessment of the performance of the prediction equation.
Agreement between FFMADP and the predicted FFM by GUSTO equation (FFMGUSTO) either using the impedance quotient (W+S+L2/R50) or simple anthropometry (W+S+L) was assessed using the intra-class correlation coefficient (ICC). Bland–Altman analysis was then used to compare FFMADP and FFM derived from prediction equations. The bias or mean difference of FFM between the two methods was used to determine underestimation (negative bias) or overestimation (positive bias) of the prediction equation. The possible extent of the variation was determined by limits of agreement (LOA) (mean of differences (sd 1·96)) in which the agreement between two methods would lie in 95 % of the repeated samples. The slope of the regression line defined whether there is a systematic variation of the difference with the mean of FFM values by two methods. Analysis was carried out using IBM SPSS 22 and MedCalc version 15 (MedCalc Software). Scattered plots and Bland–Altman plots were produced using MedCalc version 15.
The GUSTO prediction equations developed at birth and week 2 (FFMGUSTO) were also applied to an independent cohort from the University of Queensland to cross-validate the performance of the developed prediction equation in an independent Western population. This Queensland cohort included seventy-seven neonates at birth, fifty-four at 6 weeks, fifty-five at 3 months and fifty-three at 4·5 months of age( Reference Lingwood, Storm van Leeuwen and Carberry 22 ). Intra-class correlations between FFMADP and FFMGUSTO in the Queensland cohort, as well as bias and LOA, were compared with correlations between FFMADP and predicted FFM using the Queensland cohort’s own equation developed at birth.
Results
Table 1 presents the characteristics of the neonates at birth and at week 2 after delivery. The physical characteristics and body composition of the neonates were not different in the model development and validation groups at either time point.
* Age group at birth: I: day 0, II: day 1, III: day 2–3. Age group at week 2: I: <1 week, II: 1–2 weeks, III: >2 weeks.
Significant difference (P<0·05). P values were based on between-group comparisons of two sample t tests for continuous variables and χ 2 test for categorical variables.
Resistance was negatively associated with FFMADP both at birth and at week 2. The strength of the association between FFMADP and R50 was low at birth (r −0·204, P=0·007) but stronger at week 2 (r −0·438, P<0·001). The correlations between FFMADP and weight were r 0·947 and 0·946 (both P<0·001) and between FFMADP and length were r 0·750 and 0·753 (both P<0·001) at birth and week 2, respectively.
Birth
Model development group
The linear regression analysis revealed that weight and sex are both significant predictors of FFM at birth (Table 2). Weight was the strongest contributor to the prediction of FFM. The standardised coefficient for L2/R50 was greater than that for L in the models, but contribution of L2/R50 was not significant in the model including weight and sex. Weight explained 90·6 % of the variance in FFM at birth. Adding sex to the model resulted in minimal improvement in predicting FFM, whereas adding L or L2/R50 led to no improvement.
L2/R50, impedance quotients.
* P<0·05, ** P<0·01 statistically significant standardised regression coefficient from multiple linear regression.
Validation group
The scatterplot for the associations between FFMADP and predicted or FFM(W+S+L) are shown in Fig. 1(a) and (b), respectively. The ICC between FFMADP and was 0·970 (P<0·001) and for FFM (W+S+L) it was 0·971 (P<0·001). Fig. 2 shows the Bland–Altman plots for the validation group. The plots show the mean of FFMADP and predicted FFM either by W+S+L2/R50 (Fig. 2(a)) or W+S+L (Fig. 2(b)) on the x-axis and the difference between the two methods on the y-axis. The mean bias of the predicted was 0·03 kg – that is, 0·7 % of FFM at birth (2·77 kg). The LOA were 0·22 to −0·16 kg. The root mean square error (RMSE) of the prediction equation using W+S+L2/R50 at birth was 0·10 kg. The mean bias of FFM prediction by W+S+L was 0·11 kg, larger than that of the model using W+S+L2/R50. The LOA were +0·30 to −0·08 kg. RMSE for this prediction by W+S+L was 0·15 kg. No significant relationship was observed between the mean and the difference of measured v. predicted FFM.
Week 2
Model development group
Unlike the prediction at birth, the prediction models at week 2 revealed that weight and either length or L2/R50 were significant predictors of neonatal FFM at week 2 but sex was not (Table 2). Adding L2/R50 to the model explained 89·9 % of the variance in FFM at week 2 in neonates compared with 87·4 % with weight alone and 88·8 % with W+S+L.
Validation group
The scatterplots between FFMADP and that predicted FFM by W+S+L2/R50 or W+S+L in the validation group are shown in Fig. 1(c) and (d), respectively. The ICC between the FFMADP and and FFM(W+S+L) were 0·963 and 0·966 (P<0·001), respectively.
The Bland–Altman analysis shows that the mean bias of FFM (kg) at week 2 by ADP and W+S+L2/R50 was −0·01 kg – that is, 0·3 % of FFM at week 2 (2·80 kg). The LOA for the GUSTO equation at week 2 were +0·18 to −0·20 kg (Fig. 3(a)). Bias for prediction of FFM(W+S+L) was slightly larger (−0·03 kg, 0·9 % of FFM). The LOA of this prediction were 0·15 to −0·21 kg (Fig. 3(b)). RMSE for prediction was similar to prediction using L2/R50 at 0·096 kg. No significant relationship was observed between the mean and the difference of measured v. predicted FFM.
Cross-validation with the independent cohort
The GUSTO prediction equations using simple anthropometry or L2/R50 were applied to an independent cohort from the University of Queensland. When GUSTO prediction equations at birth were applied to the Queensland cohort at birth (Table 3), the ICC, bias and LOA of the predictions were very close to those values obtained in own GUSTO cohort, and again inclusion of the impedance quotient did not improve the prediction. However, when the Queensland equations were applied to the GUSTO cohort, the bias increased slightly, although the LOA were similar.
LOA, limit of agreement; ICC, intra-class correlation coefficient; GUSTO, Growing Up in Singapore Towards Healthy Outcomes; UQ, Queensland cohort.
* Bias and LOA were based on Bland–Altman analysis.
When GUSTO predictions equation developed at 2 weeks were applied at 6 weeks in the Queensland cohort, the bias and LOA were slightly smaller compared with results using the Queensland 6-week equation. However, when the GUSTO 2-week equation was applied to the Queensland cohort at 3 months and 4·5 months of age, the bias was larger and LOA were wider compared with results using the Queensland’s age-appropriate equations (Table 4). Agreement worsened with increasing age.
LOA, limit of agreement; GUSTO, Growing Up in Singapore Towards Healthy Outcomes; UQ, Queensland cohort.
* Bias and LOA were based on Bland–Altman analysis.
Discussion
A simple and reliable method is required to measure body composition in newborn infants participating in large epidemiological longitudinal studies( Reference Yajnik, Fall and Coyaji 29 ). Our findings add to the limited available evidence on the validity of BIA estimates of body composition in neonates. We developed prediction equations for FFM using BIA (W+S+L2/R50), as well as simple anthropometric measures (W+S+L), and validated these against the criterion method ADP (PEA POD). At birth, we found that L2/R50 did not contribute significantly to the prediction of FFM. Unlike the findings at birth, the prediction then improved substantially over the subsequent week or two. At 1–2 weeks of age, both L and L2/R50 became significant predictors of FFM. Their contributions to FFM prediction were similar and much greater than the non-significant contribution of sex at this age. The rapid changes in body water status and body composition in the first few days of life may impact on hydration status of FFM, thus confounding the conversion of TBW to FFM( Reference Collins, Reid and Makrides 30 ). In the Queensland cohort, L2/R50 did not become a significant predictor of FFM until 4·5 months, and its contribution was similar to the significant contribution of sex( Reference Lingwood, Storm van Leeuwen and Carberry 22 ). In addition, this finding of significant contribution of L2/R50 in GUSTO may imply that the hydration status in Asian infants is different from the Western populations studied previously.
Although the L2/R50 contributed significantly to the prediction at 2 weeks of age, it did not appear to contribute much more than that provided by simple anthropometry. There is a debate in the literature as to whether L2/R50 performs better than simple anthropometry in predicting FFM. BIA is widely used in adults, adolescents and older children but much less frequently in infants. A few small studies have suggested that BIA is useful for estimating TBW or FFM in infants, but they were performed in low birth weight infants or infants requiring intensive care( Reference Collins, Reid and Makrides 30 ), whose fluid distributions were likely to be different from healthy term infants. Moreover, those studies did not validate their equations in an independent group, nor did they compare prediction using BIA with the use of simple anthropometric measures alone( Reference Mayfield, Uauy and Waidelich 19 , Reference Tang, Ridout and Modi 21 ). Raghavan et al.( Reference Raghavan, Super and Chatburn 20 ) reported high correlation between predicted TBW using BIA and that based on a dilution technique, but the prediction equation did not perform better than did weight alone. Tang et al.( Reference Tang, Ridout and Modi 21 ) reported that weight and impedance index together improved an association (r 2) slightly v. weight alone but they did not test that finding in a validation group. Lingwood et al.( Reference Lingwood, Storm van Leeuwen and Carberry 22 ) concluded that it did not improve the prediction compared with simple anthropometry in infants younger than 3 months of age, whereas a slight improvement was present with the prediction using L2/R50 at 3 and 4·5 month of age. The study of Ethiopean infants by Wibaek et al.( Reference Wibaek, Kaestel and Skov 23 ) , although with a wide age range of infants from 0–6 months, also found that the L2/R50 improved prediction for infants older than 3 months of age. Therefore, our findings consistent with the previous findings suggest that BIA has limited value in very early life in predicting FFM, above and beyond that of anthropometry (weight and length) alone.
The failure of impedance measurements to improve the prediction of FFM at early infancy means that without DXA or ADP, prediction is limited to anthropometric measures: weight, length and skinfold thickness (SFT). However, weight or length has inherent limitations for predicting body composition, because these measures do not differentiate between fat and lean mass. The use of weight alone assumes that FM and FFM have similar proportions over time and in different populations, which is not true. Even individuals with similar weight and length can have very different body composition( Reference Yajnik and Yudkin 31 ). Measurements of SFT have been shown to correlate with body fat measured in older children and infants( Reference Goran, Driscoll and Johnson 32 – Reference Schmelzle and Fusch 35 ). A previous study in GUSTO showed that prediction of FM using subscapular SFT together with W, S and gestational age (GA), referenced to FMADP, only slightly improved prediction compared with prediction using W, S and GA only (r 2 0·811 v. 0·774)( Reference Aris, Soh and Tint 36 ). However, the validity of SFT in infants is being questioned due to rapid changes in hydration status and the variability in skinfold compressibility among neonates( Reference Davies and Lucas 13 , Reference Lohman 37 ). In addition, training of study team members, standardisation and quality control are a challenge, especially for large longitudinal studies.
Yajnik et al.( Reference Yajnik, Fall and Coyaji 29 , Reference Yajnik, Lubree and Rege 38 ) have described the thin-fat phenotype of South Asian neonates: lower birth weight but greater adiposity, reflecting different intra-uterine trajectories in fat v. lean tissues. Studies on total or regional adiposity in Indian infants compared with European infants of UK suggested that Indian-born babies have greater adiposity compared with their Western counterparts( Reference Yajnik, Fall and Coyaji 29 , Reference Nightingale, Rudnicka and Owen 39 , Reference Modi, Thomas and Uthaya 40 ). Other studies have shown that this phenotype at birth persists into childhood and later in life( Reference Krishnaveni, Hill and Veena 41 – Reference Dudeja, Misra and Pandey 43 ). In addition, the first few weeks of life may be a critical period for the development of obesity later in life( Reference Stettler and Iotova 44 – Reference Botton, Heude and Maccario 48 ). These findings may also reflect a different trajectory of FM and FFM growth in Asian populations. It would be of great value if a simple method could be developed to predict FFM beyond that of simple anthropometry during infancy.
When GUSTO prediction equations developed at birth and week 2 were applied to the respective age groups in the Queensland cohort, the performance of the prediction equations seemed to be similar to its own GUSTO cohort. When the Queensland prediction equations developed at birth were applied to the same age group of the GUSTO cohort, the results were similar with a slightly larger bias and wider LOA. Therefore, although the performance of the GUSTO equation in the Queensland Cohort may suggest that prediction equations are transferrable, results from application of the Queensland equations to the GUSTO cohort suggest that FFM prediction equations are best used in the population in which they were originally developed.
When GUSTO FFM prediction equations developed at 2 weeks were applied to 6-week-old infants of the Queensland cohort, the bias and LOA were slightly smaller compared with results using the Queensland 6-week equation. However, when the GUSTO 2-week equations were applied to 3- and 4·5-month-old infants of the Queensland cohort, the bias was larger and LOA were wider than when using Queensland’s equations at 3 months and 4·5 months of age, and agreement decreased with age. This finding implies that prediction equations are best applied in the age group for which they were developed.
A limitation of this study is that FFM density data produced by Fomon et al. ( Reference Fomon, Haschke and Ziegler 25 ) and Butte et al. ( Reference Butte, Hopkinson and Wong 26 ) were based on Western populations. It is possible that FFM density of Asian infants will be different, which may impact on body composition data derived from the PEA POD.
In conclusion, BIA appears to have limited use in predicting FFM in the first few weeks of life compared with simple anthropometry. There is still a need for population- and age-appropriate FFM prediction equations. Future studies should assess the age at which electrical impedance provides a significant improvement over simple anthropometry in predicting body composition, especially in Asian populations.
Acknowledgements
This study acknowledges the contribution of the rest of the GUSTO study group, which includes the following: Pratibha Agarwal, Arijit Biswas, Choon Looi Bong, Birit F. P. Broekman, Shirong Cai, Jerry Kok Yen Chan, Yiong Huak Chan, Cornelia Yin Ing Chee, Helen Y. H. Chen, Yin Bun Cheung, Audrey Chia Chai Kiat Chng, Mary Foong-Fong Chong, Shang Chee Chong, Mei Chien Chua, Chun Ming Ding, Eric Andrew Finkelstein, Doris Fok, Marielle V. Fortier, Anne Eng Neo Goh, Yam Thiam Daniel Goh, Joshua J. Gooley, Wee Meng Han, Mark Hanson, Christiani Jeyakumar Henry, Joanna D. Holbrook, Chin-Ying Hsu, Hazel Inskip, Jeevesh Kapur, Kenneth Kwek, Ivy Yee-Man Lau, Bee Wah Lee, Ngee Lek, Sok Bee Lim, Yen-Ling Low, Iliana Magiati, Lourdes Mary Daniel, Michael J. Meaney, Cheryl Ngo, Krishnamoorthy Naiduvaje, Wei Wei Pang, Anqi Qiu, Boon Long Quah, Victor Samuel Rajadurai, Mary Rauff, Salome A. Rebello, Jenny L. Richmond, Anne Rifkin-Graboi, Lynette Pei-Chi Shek, Allan Sheppard, Borys Shuter, Leher Singh, Walter Stunkel, Lin Lin Su, Kok Hian Tan, Oon Hoe Teoh, Hugo P. S. van Bever, Rob M. van Dam, Inez Bik Yun Wong, P. C. Wong, George Seow Heong Yeo. All the authors thank the GUSTO study team, study participants, parents and their neonates for their valuable contribution in this study.
This research is supported by the Singapore National Research Foundation under its Translational and Clinical Research (TCR) Flagship Programme and administered by the Singapore Ministry of Health’s National Medical Research Council (NMRC), Singapore – NMRC/TCR/004-NUS/2008; NMRC/TCR/012-NUHS/2014. Additional funding was provided by the Singapore Institute for Clinical Sciences, Agency for Science Technology and Research (A*STAR), Singapore.
All the authors were involved in all parts of the study and approved the final version of the manuscript. S. M. S., K. M. G., P. D. G., Y.-S. C., F. Y., Y. S. L. designed and led the GUSTO study. M.-T. T. undertook the experimental studies, undertook the sample analyses and wrote the manuscript. L. C. W. and B. L. contributed to data analysis, statistical analysis, interpretation and manuscript preparation. L. C. W., B. L., M. S. K. and Y. S. L. provided important intellectual contribution to the revision of the manuscript. A. C., S. E. S. and I. M. A. contributed to the data analysis and revision of the manuscript.
P. D. G., K. M. G., Y. S. L. and Y.-S. C. have received re-imbursement for speaking at conferences sponsored by companies selling nutritional products. P. D. G., K. M. G. and Y. S. C. are part of an academic consortium that has received research funding from Abbott Nutrition, Nestle and Danone. No other potential conflicts of interest relevant to this article were reported. KMG is supported by the National Institute for Health Research through the NIHR Southampton Biomedical Research Centre and by the European Union’s Seventh Framework Programme (FP7/2007-2013), project EarlyNutrition under grant agreement no 289346. L. C. W. provides consultancy services to ImpediMed Ltd. ImpediMed Ltd had no involvement in the conduct of this research or in the preparation of this manuscript.