Introduction

Among children under five years in low- and middle-income countries (LMICs), an estimated 250 million (44%) are at risk of not reaching their full developmental potential [, ]. A child’s developmental potential is highly influenced by individual factors including nutritional status and perinatal risk factors such as prematurity, in addition to the conditions of the home environment such as exposure to stimulation and chronic poverty []. The early years of life are a critical period of development as a brain’s neural pathways develop most rapidly and proliferate from conception to age two []. Stressors during this period including prematurity, undernutrition, trauma or stress, contribute to long-term effects on the brain’s structure []. In countries in which the prevalence of risk factors for poor childhood development outcomes such as malnutrition and nutrient deficiency are high, higher rates of disability and developmental delay are often observed [].

Children born with risk factors for developmental delay such as prematurity are a high-risk group [, ] whose outcomes in LMICs are not well understood []. Fifteen million babies worldwide are born preterm each year, with more than 80% of them in Africa and south Asia []. In total, nearly a quarter of infants born in sub-Saharan Africa in 2015 had low birth weight (LBW) []. The combined issues of prematurity, LBW and small for gestational age (SGA) are difficult to differentiate in LMICs where precise information on gestational age is often unavailable []. Even though healthcare in resource-limited countries is improving, infants born preterm and LBW are still at high risk of poor outcomes, such as death [, , ]; and for those who survive, nutritional deficiencies and developmental challenges [, ]. Kerstjens et al. found that both moderate and early preterm infants had more frequent problems with fine motor, communication, and personal-social functioning compared with full-term infants []. Prematurity and LBW contribute to intellectual disability and neurodevelopmental deficits leading to infant and childhood morbidity, particularly intellectual and learning disabilities [, ].

Rwanda has worked to reduce neonatal mortality through newborn survival initiatives, including a National Neonatal Care Protocol [], and the establishment of neonatal care units (NCUs) in every public hospital to care for sick and small newborns []. Recent data from Rwanda estimate that 6% of newborns are LBW (birth weight less than 2.5 kg) and 10% are preterm [, ]. Through the efforts to improve care for sick and small newborns, more preterm and/or LBW babies are surviving into childhood, but little is known about their long-term outcomes. The Rwanda Ministry of Health (MOH) has invested in Pediatric Developmental Clinics (PDCs) for the developmental and nutritional follow-up of high-risk babies, including those born preterm and LBW [].

Kirk et al. found poor health, nutrition and development outcomes among children born preterm and LBW when they were ages 1 to 3 years in rural Rwanda []. To our knowledge, there is no study in Rwanda comparing nutritional and developmental outcomes in early childhood of infants born with prematurity and LBW (preterm and/or LBW) to those infants born at normal gestational age and normal birthweight (term/NBW). In this study, we aim to extend knowledge by describing and assessing the differences in development at 24 to 36 months of age between these groups and to assess factors associated with poor development in preterm and/or LBW and term/NBW toddler-age children in rural Rwanda.

Methods

Study setting, design and data collection

We conducted a cross-sectional analysis, using secondary data from two studies conducted in rural communities across 11 of 30 districts in Rwanda. Data from the original two datasets were combined and are described below.

The first study collected data in 2014 on a random sample of children aged 0–11 and 24–36 months located in rural communities in 10 districts that were pre-selected based on high rates of stunting and poverty to receive the implementation of the Early Childhood Development and Family (ECD&F) intervention by Imbuto Foundation and UNICEF Rwanda (“UNICEF study”) []. The data were collected through the collaboration of the Harvard T.H. Chan School of Public Health, University of Rwanda, and Partners In Health/Inshuti Mu Buzima (PIH/IMB) from August to November 2014 to establish a baseline of children’s development prior to ECD&F program implementation. Potential participants were identified through the enumeration of households using community health workers (CHWs) registers, as CHWs are responsible for tracking all children under age five in Rwandan villages []. Following enumeration of age-eligible children, a random sample of children was identified using a random number generator in Microsoft Excel; these sampled households were then visited for household data collection resulting in a total sample of 884 children. More detailed methods are described elsewhere [].

The second study collected data from hospital registers of infants born preterm and/or LBW who were discharged alive from the Neonatal Care Unit at Rwinkwavu District Hospital (RDH) between October 2011 and October 2013 and who received routine community-based care (“Preterm/LBW study”) []. RDH is a Ministry of Health hospital located in Kayonza District that is supported by PIH/IMB. Children were excluded if they had congenital malformations, hypoxic ischemic encephalopathy, or were not preterm and/or LBW. Using geographic and family information from the hospital records, household data were collected on children alive in November to December 2014, when the children were between 1–3 years (n = 86). More detailed methods are described elsewhere [].

In the combined dataset for this study, we included all children from both databases aged 24–36 months at the time of assessment. We excluded children missing developmental status data which is the primary outcome of the study.

The same data collector teams and nutrition and development assessment tools were utilized to collect household data in both studies. Data were merged into one file for analysis; any child with status of preterm and/or LBW was included in the preterm and/or LBW group, whether the child was from the UNICEF study or the preterm/LBW study.

Definitions

Our main outcome, “developmental delay”, was defined using Ages and Stages Questionnaire version 3 (ASQ-3), a parent-reported developmental screening instrument that has been successfully used in Rwanda for providing estimates of the proportion of children whose development may be of concern [, , ]. The ASQ-3 forms are specific to the child’s age, with 21 total forms ranging from 2 months to 60 months (with corrections up to 24 months for prematurity) which ask the caregiver to report on 30 questions based on the child’s age. We defined “developmental delay” as a child whose ASQ-3 score fell below cut-points for their age in any one of five ASQ-3 domains, including gross motor, fine motor, communication, problem solving and personal-social [].

We defined “preterm” as a dichotomous variable; babies born before 37 weeks of gestational age were defined as preterm []. Babies born with a weight under 2.5kg were defined as “LBW” []. In our study, children with either preterm or LBW were defined as “preterm and/or LBW”. We defined children born at term and with normal birth weight as “term/NBW” []. In the UNICEF study, children whose caregivers reported that they were born LBW or premature were categorized as preterm and/or LBW for analysis.

Nutritional status was collected using standard anthropometric procedures and cut-offs from World Health Organization Growth Standards []. “Stunting” was defined as height/length for age z-score < –2 standard deviations (SDs) below the mean, “wasting” defined as weight-for-height/length age z-score < –2 SDs, and “underweight” defined as low weight-for-age z-score < –2 SDs. Stunting and wasting were considered “severe” if z-scores were less than 3 SDs below the mean. Mid-Upper Arm Circumference (MUAC) measurements, which is a proxy measure of nutrient reserve in muscles and fat, were reported as a continuous variable [].

“Ubudehe” is a community-based ranking of relative wealth used in Rwandan national statistics. Ubudehe, at the time of study, was categorized into 6 rankings, with 1 being lowest wealth status and 6 being highest wealth status [].

We observed an interaction between premature and/or LBW and stunting with respect to developmental delay, so we reported these factors as a combined 4-item categorical variable: “preterm and/or LBW and stunting” equals 0 if neither premature and/or LBW nor stunted, 1 if “not premature and/or LBW but stunted”, 2 if “preterm and/or LBW but not stunted”, and 3 if “both premature and/or LBW and stunted”.

Data analysis

We described patient demographics, clinical characteristics, nutritional status and developmental status using frequencies and percentages for categorical data and median and interquartile ranges for continuous data. We reported the proportion of poor developmental outcomes with 95% confidence intervals with a precision of 0.2. We assessed the relationship between patient demographics, clinical characteristics, nutritional status outcomes (stunted, underweight, wasted) and developmental (developmental delay, no delay) using Chi-squared test or Fisher’s Exact test for categorical variables and Wilcoxon rank sum or Kruskal Wallis tests for continuous variables.

We used logistic regression to identify factors associated with poor developmental outcomes built using backward stepwise procedures for factors associated at α = 0.20 in bivariate analyses. Because we observed small numbers in some cells, we used penalized likelihood logistic regression for the final model to provide robust standard errors, confidence intervals and p-values []. All factors significant at the α = 0.05 significance level were retained in the final model. The data were analyzed using Stata v.15.1 (Stata Corp, College Station, TX, USA).

Ethics

The preterm/LBW study received approval from the National Health Research Committee, Rwanda National Ethics Committee (RNEC), the Rwanda Ministry of Health, and the Boston Children’s Hospital Institutional Review Board (IRB). The UNICEF study received approval from RNEC and the Harvard T.H. Chan School of Public Health IRB. Because this is a secondary analysis of deidentified data, additional ethics review was not required; technical review was provided by the PIH/IMB research committee.

Results

Among 970 eligible children in the combined dataset, 445 (46%) met inclusion criteria of ages 24–36 months, of which 40 (9.0%) were born preterm and/or LBW and 405 (91.0%) were born term/NBW (Figure 1). We observed more females than males (55.5%), and 23 (57.5%) preterm and/or LBW and 224 (55.3%) term/NBW were females (Table 1). Of 441 children with nutritional status recorded, 217 (49.2%) had stunting, 10 (2.3%) had wasting, and 64 (14.5%) were underweight. Among 40 children born preterm and/or LBW, 30 (79.0%) had stunting, 4 (10.8%) had wasting and 15 (38.5%) were underweight compared to 187 (46.4%) with stunting, 6 (1.5%) with wasting and 49 (12.2%) underweight children born term/NBW. Based on ASQ-3 scores, 234 (52.6%) children had developmental delay. Twenty-seven (67.5%) children born preterm and/or LBW and 207 (51.1%) children born term/NBW had developmental delay according to the ASQ-3.

Figure 1 

Flow diagram of children eligible for this study, and their nutritional and developmental characteristics, from two datasets of (1) a random sample of children born term/normal birth weight (“UNICEF study”) who were aged 24 to 36 months when data was collected in 2014; and (2) children born preterm and/or low birthweight (“Preterm/LBW study”) who were discharged from the neonatal care unit (NCU) at Rwinkwavu District Hospital (RDH) between October 2011 and October 2013 and were aged 24 to 36 months when data was collected in 2014.

Table 1

Sociodemographic and clinical characteristics of infants born preterm and/or low birth weight (“Preterm/LBW study”) and infants born at term or normal birth weight (NBW) (“UNICEF study”) in Rwanda, at 24–36 months in 2014; N = 445 unless otherwise stated.

Born preterm and/or LBWBorn term/NBW

N = 40N = 405

n%n%

Child age (months)
    24–261025.08220.3
    27–29820.012330.4
    30–321435.010225.2
    33–36820.09824.2
Child sex
    Male1742.518144.7
    Female2357.522455.3
Province
    East3895.08220.3
    Kigali00.07117.5
    North00.08320.5
    South12.58320.5
    West12.58621.2
Caregiver’s relationship, N = 437
    Mother3895.037394.0
    Father00.0133.3
    Other25.0112.8
Caregiver’s education, N = 424
    None completed1025.06516.9
    Primary completed2665.029476.6
    Secondary or higher completed410.0256.5
Caregiver’s marital status, N = 394
    Not married00.06518.4
    Married1640.022663.8
    Living with partner (not married)1742.56317.8
    Divorced/widowed717.500.0
Household wealth (Ubudehe categories)
    Extremely poor00.0164.0
    Moderately poor922.58621.2
    Poor2357.520650.9
    Not poor00.010.3
    Unknown or N/A820.09623.7
Other children in household
    No other children in household820.08420.7
    1–3 other children in household2152.525763.5
    4 or more other children in household1127.56415.8
MUAC (mm) [mean, SD], N = 442143.511.615011.7
Stunted, N = 441
    Normal (z score ≥ –2 SD)821.121653.6
    Moderate (z score < –2 SD)3079.018746.4
Wasted, N = 436
    Normal (z score ≥ –2 SD)3389.239398.5
    Moderate (z score < –2 SD)410.861.5
Underweight, N = 440
    Normal (z score ≥ –2 SD)2461.535287.8
    Moderate (z score < –2 SD)1538.54912.2
ASQ-3 development status
    On-track1332.519848.9
    Developmental delay2767.520751.1

In bivariate analysis (Table 2) factors associated with developmental delay at an alpha of 0.20 included age of child at assessment (p = 0.087), child sex (p = 0.121), premature birth (p = 0.048), stunting (p = 0.002), wasting (p = 0.018) and underweight (p = 0.003). These factors were included in our full model.

Table 2

Factors associated with developmental delay on ASQ-3 in infants born preterm and/or low birth weight (“Preterm/LBW study”) and infants born at term/normal birth weight (NBW) (“UNICEF study”) in Rwanda, at 24–36 months in 2014; N = 445.

ASQ Delay, N = 234ASQ No delay, N = 211p-value

n%n%

Child age (months)0.087
    24–265863.03437.0
    27–297053.46146.6
    30–325345.76354.3
    33–365350.05350.0
Child sex0.121
    Male9648.510251.5
    Female13855.910944.1
Province0.210
    East7058.35041.7
    Kigali3954.93245.1
    North4351.84048.2
    South3541.74958.3
    West4754.04046.0
Caregiver’s relationship, N = 4370.288
    Mother21552.319647.7
    Father538.5861.5
    Other969.2430.8
Caregiver’s education, N = 4240.691
    None completed4256.03344.0
    Primary completed16350.915749.1
    Secondary or higher completed1655.21344.8
Caregiver’s marital status, N = 3940.464
    Not married3655.42944.6
    Married11848.812451.2
    Living with partner (not married)4455.03645.0
    Divorced/widowed571.4228.6
Household wealth (Ubudehe categories)0.368
    Extremely poor1168.8531.3
    Moderately poor4547.45052.6
    Poor12052.410947.6
    Not poor00.01100.0
    Ubudehe unknown or N/A5855.84644.2
Other children in household0.827
    No other children in household5155.44144.6
    1–3 other children in household14451.813448.2
    4 and more children in household3952.03648.0
Prematurity and/or low birth weight0.048
    Yes2767.51332.5
    No20751.119848.9
MUAC (mean, SD)148.011.6151.011.80.007
Stunted, N = 4410.002
    Normal (z score ≥ –2 SD)10245.512254.5
    Moderate (z score < –2 SD)13160.48639.6
Wasted, N = 4360.018
    Normal (z score ≥ –2 SD)22252.120447.9
    Moderate (z score < –2 SD)990.0110.0
Underweight, N= 4400.003
    Normal (z score ≥ –2 SD)18850.018850.0
    Moderate (z score < –2 SD)4570.31929.7

Factors associated with developmental delay in the final model included child age and stunting and history of prematurity and/or LBW, and wasting (Table 3). Older age at assessment was associated with a reduced odds of developmental delay; this was only statistically significant for the children aged 30 to 32 months compared to children in the youngest age group (24–26 months) (OR 0.55, 95% CI 0.31–0.99). Having stunting or a history of preterm and/or LBW was associated with higher odds of developmental delay in all categories compared to children with neither concern. Children with stunting but no history of prematurity and/or LBW had double the odds of developmental delay (OR 2.05, 95% CI 1.37–3.07). Children who had both histories of preterm and/or LBW and stunting had increased odds of delay compared to normal birthweight non-stunted peers, but this was not statistically significant (OR 1.68, 95% CI 0.75–3.78). Children with wasting had borderline statistically significant higher odds of developmental delay (OR 5.79, 95% CI 0.98–34.39).

Table 3

Multivariable logistic regression model with odds ratios (OR), p-value and confidence intervals (CI) for developmental delay as measured by ASQ.

ASQ Screen for Developmental Delay

FullAdjusted

OR95% CIp-valueOR95% CIp-value

Child age (months)
    24–26RefRefRefRefRefRef
    27–290.71[0.41–1.26]0.2440.71[0.40–1.24]0.229
    30–320.57[0.32–1.01]0.0550.55[0.31–0.99]0.045
    33–360.63[0.35–1.15]0.1340.61[0.34–1.10]0.098
Child sex
    MaleRefRefRef
    Female1.36[0.91–2.02]0.132
Interaction: prematurity and/or low birth weight (LBW) and stunting
    Not preterm and/or LBW and not stuntedRefRefRef
    Not preterm and/or LBW but stunted1.89[1.23–2.91]0.0042.05[1.37–3.07]<0.001
    Preterm and/or LBW but not stunted15.93[0.87–290.44]0.06216.13[0.89–291.97]0.060
    Preterm and/or LBW and stunted1.38[0.60–3.21]0.4511.68[0.75–3.78]0.206
MUAC (mean, SD)0.99[0.98–1.01]0.522
Wasted
    Normal (z score ≥ –2 SD)RefRefRefRefRefRef
    Moderate (z score < –2 SD)4.10[0.65–25.91]0.1345.79[0.98–34.39]0.053
Underweight
    Normal (z score ≥ –2 SD)RefRefRef
    Moderate (z score < –2 SD)1.33[0.67–2.66]0.415

Discussion

We found high rates of developmental delay in our study population (52.6%), and significantly more developmental delay among children born preterm and/or LBW (67.5%) compared to term/NBW born children (51.1%) at age 2–3 years. However, we also found the high prevalence of stunting in our population (49.2%) as the major driver of developmental delay in children ages 24–36 months.

While we expected preterm and/or LBW children to have a higher prevalence of developmental delay [, ] compared to term/NBW children, we also found high rates of developmental delay among term/NBW children. We found that compared to non-stunted term babies, stunted children without history of preterm and/or LBW had a significantly higher prevalence of development delay, independent of other factors. Small size at birth is one of the leading predictors of stunting globally [], and thus the high rate of growth faltering and subsequent developmental delay even among children without this perinatal risk factor in our sample is notable. This finding is similar to findings from Ethiopia that children born with intrauterine growth restriction and children born normal birth weight showed similar anthropometric outcomes (underweight and stunting) by age two []. The subsequent strong association between stunting and developmental delay is well known. A meta-analysis conducted in 15 LMICs showed that severe stunting was negatively associated with overall development across all countries []. Similarly, a study in sub-Saharan African countries showed that among the determinants of developmental delay studied, stunting was the strongest predictor of developmental delay [].

We found significantly higher rates of developmental delay at 2–3 years of age among infants born preterm and/or LBW. This increase in developmental delay was also seen in neighboring Burundi, where 27% of LBW children had some form of developmental impairment based on a screening questionnaire at age two after discharge from a neonatal care unit []. In Malawi, infants born preterm were more likely to screen positive for disability on a parent-reported assessment tool and had higher rates of developmental delay at 18 months compared to term infants (22.8% versus 10.9%) []. Similar high prevalence of developmental delay has been shown in other sub-Saharan African countries, including Ghana and Central African Republic, where 24% and 51%, respectively, of children in the general population aged 2–4 years, screened positive for disability []. These rates are also dramatically higher than the rates of delay seen in high-income countries [] and higher than we expected given that most children born premature are expected to have caught up developmentally to their term peers by age two []. However, the disparity across settings is not surprising given that a baby born in a high-income country before 32 weeks’ gestational age has twice the chance of surviving the neonatal period without later neurodevelopmental impairment than if that baby is born in South Asia or sub-Saharan Africa [].

We expected rates of stunting to be high in our study, given the high prevalence of stunting in rural Rwanda (41%) [] and the eastern African region (26%–57.7%)[], however, the 42.9% prevalence of stunting observed in our study exceeds the WHO “very high” threshold for stunting [], and deserves particular attention. Rates of stunting were especially high in children with a history of preterm and/or LBW (79%) and higher than the national average among term/NBW children as well. While the high prevalence of stunting among preterm and/or LBW children is not unexpected, this finding is alarming and consistent with other findings from this region. Van den Boogaard et al. found an 80% prevalence of stunting among infants born LBW in Burundi, compared to a 50% prevalence nationally []. Stunting and prematurity/LBW are closely linked: in a meta-analysis by Christian et al. it was found that increased odds of childhood stunting was associated with prematurity, SGA, and LBW []. Globally, SGA at birth is one of the biggest risk factors for stunting in early childhood []. However, accurately capturing SGA is difficult in LMIC contexts where accurate information on gestational age is often unknown [], which is why in our study we were not able to assess SGA given limitations in precision of gestational age measurements. Because all of our children with preterm and/or LBW who were not stunted had developmental delay, we were also unable to look at these variables separately.

The prevalence of wasting among the children with a history of term/NBW (1.5%) was similar to that in Rwanda (2%) [], but more than five times the national average among the preterm and/or LBW children (10.8%) in our study. Among all children in the study who were wasted, there was a borderline significant association with developmental delay. Exclusive breastfeeding, timely introduction and adequacy of complementary feeding, and access to effective treatment of wasting are needed for all children, but with particular emphasis on those infants born preterm and/or LBW to reduce developmental delay which could reduce the prevalence of severe wasting by 61.4%, and could also decrease malnutrition-related deaths [].

While being preterm and/or LBW was significantly associated with delay in bivariate analysis, it was not a statistically significant independent predictor of poor developmental outcomes in the final model. We suspect that the lack of significance was due to our small sample size of children categorized as preterm and/or LBW, and very high prevalence of stunting, rather than a finding actually contrary to most published literature. All of the children with preterm and/or LBW who were not stunted had developmental delay per the ASQ-3. A 2017 study assessing factors associated with developmental delay among the preterm and/or LBW study population found that children who were SGA and very LBW (birth weight less than 1.5kg) had worse developmental scores than children who were born at a normal size for gestational age or at higher weight []. However, in our study, term/NBW children were also found to have surprisingly high rates of stunting and developmental delay. Together, these findings suggest that interventions to reduce stunting and malnutrition, as well as those to prevent prematurity and LBW will be of benefit in addressing developmental delay. Targeted health and nutrition interventions are needed during the preconception period, pregnancy, lactation and the first 2 years of life, also known as the first 1,000 days [, ], to prevent preterm and/or LBW, stunting and to improve outcomes for children in the early months and years of life. There is strong global evidence for complementary interventions that are also needed to support children across Rwanda to promote positive parenting and early stimulation to support children to achieve their optimal development []. Some of these interventions include: preconception education to promote parenting skills for child development promotion, provision of complementary foods support, ensuring neighborhoods are clean and safe, and access to health and educational services []. In Rwanda, some promising models are emerging such as Sugira Muryango (Strengthen the Family)—a home-visiting intervention with a successfully completed pilot []. Projects that can be brought to scale will be needed to improve outcomes for children across Rwanda, and special attention should be paid to those born preterm and/or LBW throughout their early years.

Our study has several limitations. First, there are some differences between the specific indicators collected by the UNICEF dataset of term/NBW infants and the preterm and/or LBW dataset, and so only those indicators available in both datasets are included in this study. However, the data available for the two databases is sufficiently similar to allow comparison of key socio-demographic characteristics; our primary outcomes of child development and nutritional status were measured using the same tools and methods, and the data were collected by the same study team. Second, there is no tool for developmental evaluation that has been developed specifically for Rwandan children, and so estimates of developmental delay in our study are based on standard Western cut-points on the ASQ-3. However, the ASQ-3 has been adapted and translated into Kinyarwanda, the primary language spoken in all regions of Rwanda, and is used widely in sub-Saharan Africa. Third, there were few children in the study with prematurity and LBW; a small sample size likely affected our ability to detect risk differences between our groups. Variables such as prematurity and caregiver education that are relevant in larger studies might be less so in this study due to small sample size. Larger studies in this higher-risk population would be of great benefit to understanding risks specific to Rwanda that may be applicable to other LMICs. In the UNICEF study, mothers may have been unable to report LBW or prematurity accurately due to lack of available information. However, previous studies in similar settings have found no significant difference in parental reporting of birth size with hospital records 2–3 years after birth []. Further, mothers were better able to report size at birth when she had a facility-based delivery [], which is 91% of all births in Rwanda []. Therefore, we are confident of maternal recall of size at birth and any errors in recall of birth weight most likely would have attenuated our findings. Lastly, precise estimation of gestational age in Rwanda is a major challenge, as is the case in other LMICs, therefore we are not able to determine the difference between preterm and appropriate weight for gestational age from newborns who had intrauterine growth restriction.

Conclusions

The alarmingly high rates of developmental delay and stunting among toddlers in the study of both children with a history of preterm and/or LBW and those born term/NBW warrants attention to the nutritional and developmental needs for these children, and country-level strategies that prioritize prevention of prematurity and intrauterine growth restriction as part of stunting reduction efforts. Since stunting and wasting were both independently associated with developmental delay, a special focus of resources on developmental interventions to support children with malnutrition, and children born preterm and/or LBW who are more likely to be malnourished, would be warranted. The interaction between prematurity and/or LBW and stunting with developmental delay demonstrates the close linkages between these factors and strategies, and interventions should focus on all children at risk for these conditions especially for children in poor, rural communities.