Developmental trends in early childhood and their predictors from an Indian birth cohort.

BACKGROUND: Early childhood developmental pattern analyses not only project future cognition potential, but also identify potential risks for possible intervention. The current study evaluates developmental trends in the first 3 years of life and their predictors in a low and middle income country setting.
METHODS: Trends of early childhood development at 6, 15, 24 and 36 months of age and their predictors were explored in a longitudinal community-based birth cohort study in an urban slum in Vellore, South India. Development was assessed using the Bayley Scales of Infant and Toddler Development-III (BSID-III).
RESULTS: The birth cohort enrolled 251 children with 94, 91, 91 and 87% follow-up at 6, 15, 24 and 36 months respectively. Child development domains of cognition, language, motor and social skills showed a significant decline in scores between 6 and 36 months of age. Higher socioeconomic position (SEP) and nurturing home environment contributed to increase in cognition scores by 1.9 and 0.9 units respectively. However, stunting caused a decline in cognition scores by 1.7 units. Higher maternal cognition, higher SEP, and caregiver responsivity positively contributed to language change over time, while higher maternal depression contributed negatively. An enriching home environment, growth parameters and blood iron status had positive association with change in motor skills.
CONCLUSIONS: A triple intervention plan to enhance home environment and nurturance, early childhood nutrient supplementation, and maternal education and well-being might prevent child developmental decline in high risk settings.

Background

It is estimated that around 250 million young children living in the low and middle income countries (LMIC) risk poor development and sub-optimal livelihood potential [1, 2] because of exposure to twin risks of extreme poverty and stunting [3]. Describing developmental trajectory patterns in the LMIC not only helps in understanding developmental variations different from conventional standards, but also aids in identifying sensitive periods and specific risk factors including adverse childhood experiences. Region-specific analysis of developmental trajectory patterns can in addition help to custom-make specific developmental intervention, which can be later incorporated into public health policy.

Sensitive periods for optimal childhood development include the time period from conception to 2 years of age considered the first-1000 days of life [4] and from two to 5 years of age, the second-1000 days of life [5, 6]. Large multi-national birth cohort studies have shown complex interactions among early childhood infections, nutritional intake, family socio-economic position and home environment, affecting cognition at 2 years [7] as well as 5 years of age [8]. Risk factors including foetal growth restriction manifesting as low birth weight [9], absolute poverty [2], early childhood stunting [2, 5], micronutrient deficiency [10], poor sanitation [11], diarrhoea [12, 13], maternal depression [14] and sub-optimal home environment [15, 16] can adversely impact early child development with persisting effects on later school cognition, learning and behaviour [1, 2, 5].

Early childhood developmental trajectory analyses have been predominantly restricted to cognition [9, 17, 18], with a few evaluating other developmental domains such as social skills [19]. Group based trajectory modelling of a large Chinese birth cohort showed four different patterns of cognition paths in the first 24 months of life [9], while longitudinal social profiles of Japanese children between two and 5 years of age showed three different trajectories [19]. Development of neuropsychological functions in Indian children was described as “non-linear, heterogeneous and protracted” in a cross-sectional study [20]. Indian birth cohort follow-up studies, as part of large multinational projects, have predominantly evaluated early childhood infection/growth in relation to later age cognition [12, 21, 22] barring a few evaluating early childhood growth-development linkage [23]. There have been few studies exploring evolution of all domains of development over time in the early years, especially in a LMIC setting. The current study was planned to evaluate the trends of child developmental domains of cognition, language, social and motor skills in the first 3 years of life and analyse their predictors. It is hypothesized that socio-economic position, home environment, maternal cognition and micronutrient status will influence developmental domains.

Methods

Settings and subjects

The present study was done as an independent sub-analysis of a large prospective, multinational longitudinal birth cohort study conducted across eight different nations across the world -‘The Etiology, Risk Factors and Interactions of Enteric Infections and Malnutrition and the Consequences for Child Health and Development (MAL-ED) Network’ [24]. The Indian study site was a densely populated urban slum in Vellore, South India and covered a population of 12,000 [25]. Children were enrolled at birth by consecutive sampling and details of enrollment, exclusion criteria, follow up and details are available in previously published articles from the same birth cohort [26-28]. The original birth cohort enrollment and consequent follow-ups were approved by the Institutional Review Board of Christian Medical College, Vellore and children were recruited at each stage after informed parental consent.

Measures

Bayley scales of infant and toddler development-III (BSID-III)

The Bayley Scales of Infant and Toddler Development-III (BSID-III) assesses development in the domains of motor, language, cognition and social skills between 1 and 42 months of age [29], and details of administration are available with proposed methodology [30], as well as community level conduct already published [27, 28]. BSID-III was administered at 6, 15, 24 and 36 months of age. Individual domain quotient was calculated from the raw score using child’s chronological age.

The WAMI measure for socio-economic position

The WAMI measure with components of access to improved water and sanitation, assets, maternal education and total household income is a simplified measure of socio-economic position (SEP) developed during the MAL-ED study, [31] and details are provided in other published articles from the same cohort [26-28].

The HOME scale

The Home Observation for the Measurement of the Environment (HOME) scale (Infant/Toddler version) analyses the home environment of the child including stimulation and support and is considered the gold standard with good psychometrics [32, 33]. Both physical and interactive home environment are observed and scored in this measure, specifically mother-child interactions including responsiveness.

The modified version has six subscales consisting of total 48 items with subscales of appropriate play materials, avoidance of restriction and punishment, organization of the environment, parental involvement, responsiveness to parent, and variety in daily stimulation [34]. The HOME scale was administered by a trained social worker, who observed the home environment for 45–60 min at 6, 24 and 36 months of child’s age. Supplementary information was obtained using a mother/caregiver interview as per the MAL-ED study protocol [35]. Further details about this measure is provided in another published article from the same birth cohort [26].

Maternal cognition was assessed by Raven’s progressive matrices, a scale of non-verbal reasoning at 6–8 months of child’s age [36]. The Self Reporting Questionnaire-20 (SRQ-20), developed by the World Health Organisation to assess depressive symptoms in low-resource settings, was used to assess maternal psychological disturbances at 1,6, 15, 24 and 36 months of child age [37]. The total SRQ score was calculated using 16 items, which yielded one factor structure [38].

Blood collection

Blood samples were collected at 7, 15, 24, 36 and 60 months; and details of collection, testing mechanisms and methods are published [28]. Samples were tested for hemoglobin (7, 15, 24, and 36 months); ferritin assays (at 7, 15 and 24 months) and blood lead levels (at 15, 24 and 36 months). As inflammation can influence serum ferritin levels, both transferrin (R) and ferritin (F) levels were utilized to find total body iron levels using the formula, where positive values indicated iron reserves [39]:

Data entry

Data entry was made in the double entry database system managed by Data Coordinating Centre (DCC) of the MAL-ED study. Data collection forms filled by the field workers were authenticated by the field supervisor before entering into the database [24].

Statistical analysis

The outcome variable, BSID scores, assessed at 6, 15, 24 and 36 months of age was summarized using mean and standard deviation (SD) under specific domains namely - cognition, language, motor and socio-emotional domain. Independent predictors such as maternal depression score, maternal cognition and various domains of HOME inventory scale were summarized as mean and SD. Height-for-age and weight-for age scores below − 2 SD were categorized as stunted and underweight respectively. Values of blood iron measured at 7, 15 and 24 months were considered for analysis and the missing values were replaced with the average of measurements available at the other time points. Blood lead levels measured at 15, 24 and 36 months were averaged to obtain mean lead level. Domain-wise score of cognitive development was compared across the time points using repeated measures ANOVA test. Factors predicting cognitive development scores across the four time points were measured using the generalized estimating equations (GEE) – population averaged model based on the following equation,

Cognition development = αi + β1 Gender + β2 Stunted + β3 Underweight + β4 Mother’s cognition scores+ β5 SEP + β6 Body iron + β7 Mean blood lead + β8 Mother’s depression scores + β9–14 Domains of HOME Inventory scale + ε.

where ‘i’ and ‘j’ refers to children and time points of measurements respectively, and, εij representing the random error. Multicollinearity of the covariates included in the model was tested using variance inflation factor and none of the variables were observed to be collinear. GEE model was used to adjust for clustering at the subject level, because of repeated measurements on children. Subject ID of study children was used as the clustering variable in the model. We specified exchangeable correlation structure in the model, based on the quasi Information Criteria (QIC) and robust standard errors were estimated [40]. Beta co-efficients along with 95% confidence interval have been reported for the independent predictors. Model fit was assessed using Wald statistics and p less than 0.05 was considered as statistical significance. Data analysis was done using Stata version 13 (StataCorp. 2013. Stata Statistical Software. Release 13. College station, TX: StataCorp LP).

Results

In the original birth-cohort, 251 newborns were enrolled after registering 301 pregnant mothers. 50 children were excluded as per the exclusion criteria: another sibling registered in the MAL-ED study (n = 8), medical comorbidities in children (n = 7), multiple pregnancy (n = 1), pre-existing plan to migrate (n = 5), mother not available for consent (n = 9), combination of two or more of the above mentioned reasons (n = 10) and mothers/parents refused participation (n = 10). Subsequent follow ups at 6, 15, 24 and 36 months had 235, 229, 228 and 218 children respectively. Migration was the main cause for the loss to follow-up, as illustrated in other published articles from the same cohort [28].

The Vellore birth-cohort had a mean birth weight (SD) of 2.89 (0.44) kg and 17% of newborns weighed less than 2.5 kg. More than 80% of families had monthly income less than 5000 Indian rupees (70 USD). The cohort had girl predominance at birth (56%) as well as subsequent follow-ups. Cohort characteristics at 6 and 36 months are summarized in Table 1, with details of 6, 15, 24 and 36 months in Supplementary Table 1.

Table 1

Baseline characteristics of the birth cohort established in Vellore in 2010

Variables	6 months (n = 235)	36 months (n = 218)
Sex (%)	Female: 129 (54.89%)	116 (53.21%)
HOME total score, mean (SD)	40.77 (3.24)	40.27 (3.18)
Domain-wise score of HOME inventory score, mean (SD)
Emotional and verbal responsivity of caregiver,	10.60 (1.10)	11 (0.07)
Avoidance of restriction and punishment	4.68 (0.52)	4.12 (0.46)
Caregiver promotes child development	4.67 (0.69)	4.96 (0.20)
Organization of physical and temporal environment	10.32 (1.37)	9.49 (1.49)
Provision of appropriate play materials	2.24 (0.69)	2.66 (0.48)
Opportunities for variety in daily stimulation	5.38 (1.57)	5.75 (1.56)
Cleanliness of child	2.87 (0.43)	2.28 (0.85)
WAMI score, mean (SD)	0.453 (0.15)	0.54 (0.17)
Maternal depression (SRQ) score, mean (SD)	4.53 (3.74)	3.82 (3.53)
Maternal cognition raw score, mean (SD)	43.90 (10.49)	–

Developmental domain quotients of cognition, language, motor and social skills significantly differed between 6 and 36 months of age (Table 2) with both cognition and social domains showing decline of more than 15 points (Fig. 1).

Table 2

Mean (SD) domain quotient scores of Bayley Scales of Infant Development -III represented domain-wise across 6, 15, 24 and 36 months in MAL-ED children (N = 216)

Domains	6 months	15 months	24 months	36 months	P-value
Cognition	102.43 (9.93)	101.48 (9.57)	92.52 (6.91)	86.53 (4.55)	< 0.0001
Language	98.0 (8.38)	98.10 (8.90)	98.05 (8.98)	93.13 (4.44)	< 0.0001
Motor	106.04 (14.08)	100.11 (8.51)	103.26 (8.69)	99.19 (7.19)	< 0.0001
Social	115.69 (13.4)	116.74 (11.64)	113.66 (9.98)	98.66 (4.13)	< 0.0001

Developmental quotient scores scores under each domain were compared across time points using repeated measures ANOVA test

Fig. 1

Domain wise scores obtained in Bayley scales of infant development assessed between 6 and 36 months of age in children of MAL-ED cohort, Vellore

Analysing factors responsible for the change in cognition scores between 6 and 36 months, higher SEP, body iron levels and HOME factors of parental responsivity and provision of appropriate play materials positively predicted cognition (Table 3). The maximum positive association was with SEP followed by HOME factors, while stunting had the maximum negative impact. Higher maternal cognition, higher SEP, and caregiver responsivity positively contributed to language change over time, while higher maternal depression contributed negatively. Maximum effect size on language change was seen with SEP.

Table 3

Factors associated with child development scores measured at 6, 15, 24 and 36 months in Vellore cohort of MAL-ED study (N = 216)

Factors	Adjusted beta co-efficient with 95% CI
	Cognitive domaine	Language domaine	Motor domaine
Female	0.51 (−0.76–1.77)	1.14 (− 0.07–2.35)	1.54 (− 0.09–3.16)
Mother’s cognition score	0.03 (− 0.03–0.09)	0.09 (0.03–0.15)	0.09 (0.01–0.17)
Stunted (HAZ < − 2SD)a	− 1.66 (− 3.18 - -0.14)	0.13 (− 1.41–1.66)	− 1.95 (− 3.63 - -0.26)
Underweight (WAZ < − 2SD)b	−0.54 (− 1.91–0.82)	−1.34 (− 2.80–0.13)	−2.41 (− 4.12 - -0.70)
SEP (WAMI index)c
Middle (33rd – 66th percentile)	1.08 (− 0.55–2.71)	1.39 (−0.14–2.79)	0.83 (− 0.96–2.63)
High (>66th percentile)	1.86 (0.28–3.43)	1.73 (0.22–3.24)	1.59 (− 0.35–3.52)
Body iron levels	0.18 (0.06–0.29)	0.05 (− 0.05–0.14)	0.30 (0.17–0.43)
Mean blood lead levelsd	−0.03 (− 0.13–0.07)	−0.07 (− 0.18–0.04)	−0.04 (− 0.18–0.11)
Mother’s depression scores	0 (− 0.18–0.17)	− 0.23 (− 0.41 - − 0.04)	-0.04 (−0.25–0.17)
Domain scores of HOME Inventory scale
Emotional and verbal responsivity of caregiver	0.92 (0.13–1.97)	1.0 (0.23–1.87)	1.30 (− 0.03–2.63)
Avoidance of restriction and punishment	0.68 (− 0.56–1.92)	0.21 (− 1.03–1.45)	1.91 (0.46–3.36)
Caregiver promotes child development	1.03 (−0.89–2.94)	1.09 (− 0.47–2.66)	0.68 (− 1.63–2.98)
Organization of physical and temporal environment	0.11 (−0.32–0.55)	0.39 (− 0.49–1.26)	− 0.12 (− 0.73–0.49)
Provision of appropriate play materials	0.98 (0.06–1.94)	0.39 (− 0.49–1.26)	1.38 (0.09–2.66)
Opportunities for variety in daily stimulation	−0.02 (− 0.46–0.41)	0.14 (− 0.12–0.83)	−0.40 (− 0.94–0.15)

aChildren with height-for-age Z scores < − 2 SD were classified as stunted

bChildren with weight-for-age Z scores < − 2 SD were classified as underweight

cWAMI index was based on water and sanitation facility, household assets, mother’s education, and total income

d Mean blood lead levels were obtained by taking average of lead values measured at 15, 24 and 36 months of age

eDependent variables considered for analysis were cognition, language, and motor domains of BSID-III scores. Regression co-efficients were derived using generalized estimating equations and the model fitness was assessed for each domain separately using Wald statistics and P value was noted as < 0.0001

Factors positively contributing to change in motor scores included higher maternal cognition, higher weight and height for age scores, higher body iron levels and positive home factors of avoidance of restriction and punishment, and availability of appropriate play materials. Both HOME factors and growth parameters had high association with change in motor skills. Body iron levels positively influenced social domain change over time with beta-coefficient (95% CI) of 0.20 (0.05–0.34), while blood lead levels had a negative influence with beta-coefficient (95% CI) of − 0.22 (− 0.34- -0.09) in univariate analyses. There were no significant associations in multivariate analysis, thus not reported in Table 3.

Discussion

This prospective longitudinal birth cohort follow-up analysis from urban Vellore evaluated developmental trends in the domains of cognition, language, motor and social skills; and predictors for individual domain change over time. All developmental domain quotients showed a decreasing trend over time with cognition and language domains dropping more than 15 points between 6 and 36 months of age. Children in the highest tertile of SEP had better cognition and language scores over time. Stunted children had poorer cognition and motor scores over time. Home environmental factors of parental responsivity affected changes in cognition and language scores, while punishment avoidance and toy availability influenced the motor score change. Higher maternal cognition affected changes in both language and motor domains, while maternal depression adversely affected language scores over time. Mean body iron levels was associated with changes in both cognition and motor skills. There was a good level of follow-up in the current study with 94% at 6 months, 91% at 24 months and 87% at 36 months.

Child development evolves over time and developmental measures such as BSID-III have shown weak predictability of school age cognition [41], motor skills [42], and behaviour [43]. Despite this, BSID is the most common developmental tool used across the globe, as it is a sensitive measure of child development [29]. Though child developmental process and sequence are uniform across the world, there can be differences in the rate of achievement between populations. Trends of all developmental domains in an urban low-income setting as shown in this study can be useful for clinical practice in such settings, academic understanding, and analysis of risks contributing to any setback. Though the current analysis has not done a complete trajectory exploration, the birth cohort showed a decreasing trend in all developmental domains’ quotients. This is in discordance with another large rural birth cohort study reported from China, where more than 90% showed an increasing trend in cognitive development between 6 and 24 months of age [9]. The Vellore birth cohort is predominantly a low-income urban slum setting with additional environmental and nutritional challenges, which have been explored in the current analysis. Similar results were reported from another study conducted in rural India, where recruited infants showed a decline in scores in fine motor, receptive and expressive language skills and visual reception over a 6-month follow-up period [23]. A previous birth cohort follow up done in other urban slums of Ramnaickanpalayam, Chinnallapuram and Kaspa in Vellore had also revealed decline in cognition scores between three and seven years of age [44].

The SEP including household wealth affects linear growth [45, 46], development [2] and cognition [15, 47] in children. Childhood SEP is shown to influence neural development particularly of the language and executive function areas in the brain [48]. Better SEP can result in improved nutrition, better sanitation, less infection and enhanced interactive experiences, all of which can improve child development especially language as shown in the current study.

Stunting had independent associations with cognition and motor scores despite correction with SEP, HOME status and micronutrient levels in the current study. Stunting has been shown to affect development and cognition; thus integrated interventions have been suggested to improve both linear growth and child development [2, 5]. A rural study conducted in Telangana, South India found similar results where height z-scores at enrollment had a positive association with all child development domains, the effect of which was attenuated by a nurturing home environment [23].

In the current study, positive home environmental factors influenced developmental gains over time, despite correction with SEP. Home environment can mediate the relationship between SEP and early childhood development [15, 16] as well as between linear growth and development [23]. The two components of home environment – physical and relational modulate each other, with the sensitive and interactive relational factor overriding the negative effects of a sub-optimal home milieu. The effect of physical overcrowding on childhood development and cognition can be mediated by maternal responsiveness [49].

Maternal factors can influence child development, cognition and behaviour directly and indirectly through better home environmental factors. The effect of maternal depression can even start in the ante-natal period not only through altered placental functions, but also by causing foetal epigenetic changes and stress reactivity [14]. Maternal depression as well as lack of breastfeeding practices can lead to adverse paediatric development [50]. Maternal cognition also affects childhood development and cognition through genetic factors as well as enriching interacting experiences [8, 15, 49].

Iron deficiency in early childhood, when there is an increased iron demand to optimize neuronal maturation, neurotransmitter synthesis, mitochondrial function and other iron dependent enzymes [51], can be detrimental as shown in the current analysis. This early childhood iron deficiency can have persisting effects on cognition in later life as evidenced by another analysis on the same cohort showing cumulative iron deficiency negatively impacting verbal, performance and processing speed domains of cognition at 5 years of age [28]. That assessment had shown that more than 40% of children had iron deficiency at 15 and 24 months of age. Later childhood iron supplementation improves body iron stores, but not the persisting effect of early onset iron deficiency on cognition [51, 52].

There are limitations for the current study including a comparatively small sample size. Though BSID-III is a sensitive and adapted measure, there can be variations (within normal limits) in developmental achievements. However the low-income urban slum setting of the cohort, its high early childhood iron deficiency and a significant decline in developmental quotients over time warrant an assessment as in the current study. Strengths of the study include good follow-up of a longitudinal birth cohort, strong data granularity during early childhood period, standardized developmental, home environment and SEP assessments, blood results from a national reference laboratory and rigorous quality control measures internally and externally.

Conclusions

This longitudinal birth cohort follow-up study has shown a significant decline in developmental scores in cognition, language, motor and social domains between 6 and 36 months of age in a LMIC urban slum setting. Positive home environmental factors in addition to SEP enhanced developmental outcomes. Both maternal components of cognition and depression showed effects on progress of child development. Linear growth as well as iron deficiency independently influenced development.

Findings of the current study can be relevant for other high risk settings worldwide. The study confirms three large domains associated with childhood development - nurturing environment and SEP under environmental factors, maternal factors of education and depression, and nutritional factors of iron deficiency and appropriate linear growth. Information about appropriate early childhood nutrition and nurturing activities can be provided in both antenatal and immunization clinics to empower mothers. Community level support groups for new mothers can help and support mothers in the post-partum period, while iron supplementation and fortification can optimize early childhood iron deficiency. A triple intervention plan to enhance home environment and nurturance, maternal education and well-being, and early childhood nutrient supplementation might prevent child developmental decline in high risk settings.

Additional file 1: Supplementary Table 1. Comparison of baseline characteristics of the birth cohort established in Vellore in 2010.