Physical growth and neurodevelopment during the first year of life: a cohort study of the Japan Environment and Children’s Study

Background The association between a slower physical growth and poorer neurodevelopment has been established in infants born preterm or small for gestational age. However, this association is inconsistent in term-born infants, and detailed investigations in infancy, when intervention is most beneficial for improving outcomes, are lacking. We therefore examined this association separately by sex during the first year of life in term-born infants. Methods Using data collected until children reached 12 months old in an ongoing prospective cohort of the Japan Environment and Children’s Study, we analyzed 44,264 boys and 42,541 girls with singleton term-birth. The exposure variables were conditional variables that disentangle linear growth from weight gain relative to linear growth, calculated from the length and weight at birth and 4, 7 and 10 months old. Neurodevelopmental delay was identified using the Japanese-translated version of Ages & Stages Questionnaires, third edition. Results A reduced risk of neurodevelopmental delay at 6 months old was observed in children with a higher birth weight (adjusted relative risks [aRRs]: 0.91 and 0.93, 95 % confidence intervals [95 % CIs]: 0.87–0.96 and 0.88–0.98 in boys and girls, respectively) and increased linear growth between 0 and 4 months old (aRRs: 0.85 and 0.87, 95 % CIs: 0.82–0.88 and 0.83–0.91 in boys and girls, respectively). A reduced risk at 12 months was found in children with an increased linear growth between 0 and 4 months (aRRs: 0.92 and 0.90, 95 % CIs: 0.87–0.98 and 0.84–0.96 in boys and girls, respectively), boys with an increased relative weight gain between 0 and 4 months (aRR: 0.90, 95 % CI: 0.84–0.97), and girls with a higher birth weight (aRR: 0.89, 95 % CI: 0.83–0.96). Conclusions These results suggest that a slow physical growth by four months old may be a predictor of neurodevelopmental delay during infancy. Supplementary Information The online version contains supplementary material available at 10.1186/s12887-021-02815-9.


Background
Prenatal and postnatal growth of the body has been linked to later neurocognitive development. Growth restriction in utero is suggested to occur in infants born small for gestational age while poor growth from birth to term-equivalent age is often observed in infants born preterm. Such growth-faltering infants are at high risk of poor neurodevelopmental outcomes later in life [1,2]. In these infants, a lack of catch-up growth is an important indicator of a poor neurodevelopment [3][4][5]. The association between physical growth and neurodevelopment has also been examined in general populations consisting mostly of term-born infants, but the findings are less consistent in such settings. Indeed, while several studies with limited sample sizes have shown positive results [6][7][8][9][10][11][12][13][14], others have reported no association [15,16]. Furthermore, as outcomes, a majority of these studies investigated the neurocognitive function during preschool through adulthood [6][7][8][9][10][13][14][15][16] with only a few examining those in infancy [8,11,12], during which intervention has the potential to improve the outcomes of both the parents and their children [17,18].
Most of these previous studies have separately assessed anthropometric gains in either length/height or weight between two age points (e.g. birth to six months old). As the length/height and weight are strongly correlated within an individual, assessing these indices per se hampers the determination of which component is more influential. To dissociate the effects of length/height growth from those of weight growth, Adair et al. and Horta et al. used a sophisticated method of 'conditional growth modeling' [19,20]. This modeling determines conditional variables by regressing the current size (i.e. length or weight) against all previous sizes [21,22]. With respect to each period, the variables comprise 'linear growth' and 'weight gain relative to linear growth', which represent length/height change and weight change separated from change in length/height, respectively. Their studies demonstrated a closer relationship of school achievement and intelligence with linear growth than with relative weight again, noting some differences in patterns between sexes.
Using a conditional growth model with a nationwide prospective cohort from the Japan Environment and Children's Study (JECS) that includes over 100,000 children, the present study examined the relationship between physical growth and neurodevelopment separately by sex during the first year of life in term-born infants. To collect neurodevelopmental information from such a large number of children, we used a parent-reported screening tool rather than individualized face-to-face tests. During this period, neurodevelopment is difficult to assess precisely, and important changes can manifest in the individual developmental trajectory. However, a detailed investigation of this relationship during infancy is extremely important for detecting neurodevelopmental delay for early intervention.

Design
The JECS is a nationwide, multicenter, prospective birth cohort study funded by the Ministry of Environment, Japan. The study design details have been described elsewhere [23,24]. In brief, pregnant participants were registered between January 2011 and March 2014 in 15 regional centers covering a wide geographical area in Japan. During pregnancy, data were obtained during the first and second/third trimesters using self-administered questionnaires. Detailed information regarding the mother and child was obtained from medical records transcripts during the first trimester, at the time of delivery, and when the child was 1 month old. After delivery, data were collected at one and six months old and every six months until the child was six years old, then twice a year thereafter via self-reported questionnaires completed by the parents.
The JECS protocol was reviewed and approved by the Ministry of Environment's Institutional Review Board for Epidemiological Studies and by the ethics committees of all participating institutions (#15000141). The ethical approval for this study was an extension of the ethical approval for the JECS protocol. Written informed consent was obtained from all parents.

Participants
In this study, we used the fixed dataset "jecs-an-20180131" released in March 2018. The dataset includes all available data of 104,065 fetuses, linked to their mothers' data, collected until children were 12 months old. As subjects, we selected 89,953 children with live singleton term birth (≥ 37 and < 42 gestational weeks) who had parents of Japanese nationality and complete basic information on mother's parity, child's sex, length and weight at birth (Fig. 1). After excluding 3,148 children who had malformation/severe diseases, we eventually analyzed 86,805 children, including 44,264 boys and 42,541 girls.

Exposure
The exposure variable was body growth from conception to 10 months old. Length and weight at birth, which can be viewed as growth from conception (length and weight = 0) to birth, were obtained from medical records. Subsequent anthropometric data were all derived from the personal health record of the child. Nearly all Japanese children receive a public health checkup at four months old. Their length and weight were thus directly measured by a public health nurse or a hospital nurse and officially documented in the record. They also receive similar checkups at 7 or 10 months old, occasionally at both 7 and 10 months old, depending on residency. Parents were asked to transcribe these anthropometric data at the 4-month checkup in the 6month questionnaire and those at the 7-and 10-month checkups in the 12-month questionnaire. Because the data were not measured at exactly 4, 7 and 10 months, the values were extrapolated from the anthropometric data at birth and those at the 4-, 7-and 10-month checkups [25]. The extrapolated data were referenced to the Japanese growth standards separately by sex and then converted into z-scores (i.e. length-for-age, weightfor-age) [26,27]. In the standards, the mother's parity and child's gestational age are additionally taken into consideration when calculating z-scores at birth.
These z-scored length and weight values at different ages were processed using conditional growth modeling, as in the studies conducted by Adair et al. and Horta et al. [19,20] but additionally including birth length in our study. Conditional variables were obtained by regressing the current size (i.e. length or weight) against all previous sizes [21,22]. At each time point, a conditional variable represents physical growth during a time interval, and a positive value represents a linear growth or relative weight gain faster than predicted in that period. For example, the conditional length at seven months old indicates the linear growth from four to seven months old. To estimate the conditional length, the current length was regressed against all previous lengths and weights. The conditional length at seven months old was derived by regressing the length at seven months old against the length and weight at four months old and birth. To estimate the conditional relative weight, the current weight was regressed against the current length and all previous lengths and weights. For instance, the conditional relative weight at seven months old was calculated by regressing the weight at seven months old against the length at seven months old and the length and weight at four months old and birth. Note that the length and weight at birth per se were treated as conditional variables at birth, as they had no previous anthropometric data to be regressed against. Consequently, our conditional variables included length and weight at birth and conditional length and relative weight at 4, 7 and/or 10 months old.

Outcomes
The outcome measure was neurodevelopmental delay at 6 and 12 months old using the Japanese-translated version of the Ages & Stages Questionnaires: third edition (ASQ-3). This version was prepared through a backtranslation procedure and was approved by the publisher of the original English version [28]. The ASQ-3 can identify infants or young children who require further neurodevelopmental assessments to determine their eligibility for early intervention. The findings of the questionnaire basically agree with those of professionally administered developmental batteries [29,30]. It has been used in clinical and research settings and translated into several languages [31][32][33].
The ASQ-3 assesses five developmental domains: Communication, Gross Motor, Fine Motor, Problem Solving, and Personal-Social domains. For each domain, six skills are described to which parents answer "yes," "sometimes," or "not yet," depending on whether or not their child has demonstrated the described skill. The responses are converted to points, with "yes" receiving 10 points; "sometimes", 5 points; and "not yet", 0 points. The child's score for each developmental domain is the sum of all points received for the items under that domain and ranges from 0 to 60 points. The cut-off score for each domain is defined as two standard deviations below the mean score of large, standardized samples in the United States of America. A child is defined as having a neurodevelopmental delay if a score is at or below the cut-off level in any developmental domain. Although preliminary cut-off scores of the Japanese translation were recently proposed [34], these were not recommended to be used with confidence before 24 months old because of very limited sample sizes. Therefore, the cut-off scores were determined using the same methodologies as the original version, based on available data at 6 (n = 82,994) and 12 months old (n = 78,947) ( Fig. 1), which represent the general Japanese population.

Statistical analyses
To examine which component and periods of physical growth are specifically associated with infant neurodevelopment, we conducted multivariable regression analyses using conditional variables at different ages. Because these variables are not correlated with each other, the variables and covariates were included together in the same model without concerns about collinearity [20]. The covariates were (i) mother's age; (ii) maternal smoking during pregnancy, as recorded in the first trimester; (iii) maternal and (iv) paternal education level (junior high school, high school, and university or graduate school); (v) annual family income (< 4 000 000, 4 000 000-5 999 999, ≥ 6 000 000 JPY); (vi) home speech stimulation at 1 month old (whether a mother did or did not talk to her baby habitually). The 'home speech stimulation' covariate was used instead of the Home Observation for Measurement of the Environment scale [35], which is not employed in the JECS.
All analyses were performed using the R software program, version 4.0.3 (The R Foundation, Vienna, Austria). Quasi-Poisson regression models were used to estimate adjusted relative risks (aRRs) with 95 % confidence intervals (CIs).
In these children, neurodevelopmental delay at 6 months old was observed at similar rates in boys (8.9 %) and girls (8.5 %), whereas that at 12 months old was more frequently identified in boys (16.9 %) than in girls (12.6 %). The baseline characteristics of the children who did not have all data available were comparable to those of children with complete data, but their parents tended to more often be multipara, have a smoking habit and have less education and income than those of complete-data subjects ( Table 1). The comparisons between children with and without anthropometric data using other combinations among 4, 7 and 10 months old showed similar tendencies (Table S1, S2, S3).
Multivariable regression models revealed that a reduced risk of neurodevelopmental delay at 6 months old was found in boys with a higher birth weight (aRR: 0.91, 95 % CI: 0.87-0.96) and an increased conditional length at 4 months old (aRR: 0.85, 95 % CI: 0.82-0.88) ( Table 2), and similarly in girls with a higher birth weight (aRR: 0.93, 95 % CI: 0.88-0.98) and an increased conditional length at 4 months old (aRR: 0.87, 95 % CI: 0.83-0.91) ( Table 3). On analyzing each neurodevelopmental domain of ASQ-3, the conditional length at four months old was associated with all five domains, while the birth weight was associated with some domains, supporting the stronger effect sizes of conditional length at four months old than those of birth weight for both sexes.
With regard to neurodevelopmental delay at 12 months old, we first created a regression model using the data of children who had complete anthropometric data throughout the first year of age and all covariates (Tables 4 and 5, model 1). In boys, a reduced risk was related to an increased conditional length at 4 months old (aRR: 0.92, 95 % CI: 0.87-0.98) and an increased conditional relative weight at 4 months old (aRR: 0.90, 95 % CI: 0.84-0.97) ( Table 4). In girls, a reduced risk was related to a higher birth weight (aRR: 0.89, 95 % CI: 0.83-0.96) and an increased conditional length at 4 months old (aRR: 0.90, 95 % CI: 0.84-0.96) ( Table 5). Analyses of each neurodevelopmental domain revealed that, for both sexes, conditional length at four months old was most consistently associated with several domains compared with other conditional variables, again suggesting the importance of this variable. Conversely,  Data represent the adjusted relative risk (95 % confidence interval). All of the conditional variables and the following covariates were included together in the model: maternal age and smoking, maternal and paternal education, family income, and home speech stimulation. cLength conditional length; crWeight conditional relative weight the risk was largely unrelated to conditional variables at 7 and 10 months old, although an increased risk was paradoxically observed in girls with an increased conditional length at 7 months old (aRR: 1.16, 95 % CI: 1.07-1.26) ( Table 5). We also built 2 other regression models using the data of children with available data from 4 to 7 months old (model 2) and those with available data from 4 to 10 months old (model 3). The results of the two models included some inconsistencies but yielded similar results concerning the importance of conditional variables in early infancy.

Discussion
The present study investigated the relationship between physical growth and neurodevelopmental delay during the first year of life in term-born infants. As indices of physical growth, we discriminated linear growth from relative weight gain, represented as 'conditional length' and 'conditional relative weight', respectively. A reduced risk of neurodevelopmental delay at six months old was observed in boys and girls with a higher birth weight and an increased conditional length at four months old. A reduced risk of neurodevelopmental delay at 12 months old was found in boys and girls with a higher conditional length at 4 months old, boys with a relatively high conditional relative weight at 4 months old and girls with a higher birth weight. These findings indicate that children with a poor growth in early infancy are at an increased risk for neurodevelopmental delay. In contrast, conditional variables at 7 and 10 months old did not increase the risk. This was one of the largest birth cohort studies yet conducted to examine the association between physical growth and the neurodevelopmental outcome. Our large sample size had strong statistical power and revealed subtle but significant associations. Furthermore, body sizes were converted to z-scores using domestically standardized growth curves, and the sizes at birth were further adjusted for the mother's parity and child's gestational age, enabling the precise estimation of the body growth in utero. However, the present study included several limitations. First, the information was mainly collected via self-administered questionnaires. Although the anthropometric data were based on official records, the identified neurodevelopmental delay might be somewhat equivocal, relying on the parent-reported screening test of the Japanese version of ASQ-3. Next, a considerable proportion of children had missing information on physical growth at different time points and neurodevelopmental outcomes and were thus excluded from the analysis, possibly producing selection bias. Finally, the results of our regression analyses were not adjusted for other unmeasured confounders, including the maternal intelligent quotient, which greatly influences the neurocognitive development of offspring.
We found that a higher conditional length at 4 months old was consistently associated with a reduced risk of neurodevelopmental delay at 6 and 12 months of age for both sexes. The results agree with the findings from previous studies using conditional modeling, which reported a greater impact of linear growth than relative weight gain on schooling and intelligence [19,20]. The biological mechanisms underlying this association remains unclear but is likely multifactorial. Prenatal and postnatal nutritional deficiency can induce growth failure of the body and brain, as observed in low-income countries [36]. Although severe malnutrition is less likely to occur in Japan than in other countries, placental dysfunction and maternal excessive dieting may cause a poor nutrient supply to the fetus, while inappropriate feeding habits might produce nutritional insufficiency during infancy. Another possible explanation is hormonal factors, such as the growth hormone/insulin-like growth factor-I pathway and the hypothalamic-pituitary-gonadal axis. These endocrine systems are regulated by nutritional and genetic factors [37][38][39][40][41] and play key roles in both the physical and cerebral growth of the fetus and child [42][43][44][45][46]. Subtle differences in such hormonal levels Data represent the adjusted relative risk (95 % confidence interval). All of the conditional variables and the following covariates were included together in the model: maternal age and smoking, maternal and paternal education, family income, and home speech stimulation. cLength conditional length; crWeight conditional relative weight  Data represent the adjusted relative risk (95 % confidence interval). All of the conditional variables and the following covariates were included together in the model: maternal age and smoking, maternal and paternal education, family income, and home speech stimulation. a,b The values of the variables were different between the models because these variables were estimated from different combinations of the previous data ([model 1] at birth, 4 and 7 months old; [model 3] at birth and 4 months old). cLength conditional length; crWeight conditional relative weight might be responsible for the association between physical growth and neurodevelopment. We also found other associations between physical growth and neurodevelopmental delay during infancy. The association between a higher birth weight and a reduced risk of neurodevelopmental delay at 6 and 12 months old is concordant with the well-known finding that birth weight is a significant determinant of later neurodevelopment [1,[47][48][49]. The boy-specific association between the conditional relative weight at 4 months and neurodevelopmental delay at 12 months was in line with the association between conditional relative weight and school achievement reported by Adair et al. [20]. The aforementioned endocrine factors have differential effects on somatic and cognitive development between sexes [46,50] and might account at least partially for this association. Conversely, higher conditional variables at 7 and 10 months old did not reduce the risk of neurodevelopmental delay at 12 months, suggesting less susceptibility of the developing brain in late infancy than in early infancy. An increased risk at 12 months was also paradoxically observed in girls with an increased conditional length at either 7 or 10 months old in all the regression models. When an infant shows slow growth in early infancy, they may be inclined to be fed more by their parents and thus may grow faster thereafter. Such modification of parents' feeding behavior might yield this paradoxical association, although this information was not collected in the JECS.

Conclusions
Using appropriate conditional modeling, we showed that poor early physical growth was associated with an increased risk of neurodevelopment delay during the first year of life in term-born infants partly in a sexdependent manner. If a child presents with growth failure in early infancy, the infant should be monitored carefully to detect neurodevelopmental delay for early intervention. Ongoing investigations in the JECS cohort may reveal later-term outcomes.