Introduction
Notable progress has been made in the reduction of child undernutrition in recent decades. However, child stunting remains an important global health issue [1]. In 2014, 157 million out of 667 million children under the age of 5 years were estimated to be stunted, over half of whom were in Asia [1]. Understood to be a direct result of poor nutrition and chronic infection, stunting is associated with a number of later life adverse outcomes, including cognitive impairment, lower educational achievement, adverse pregnancy outcomes among women and increased cardiometabolic disease risk [Reference Adair, Fall, Osmond, Stein, Martorell, Ramirez-Zea, Sachdev, Dahly, Bas, Norris, Micklesfield, Hallal and Victora2–Reference Black, Victora, Walker, Bhutta, Christian, de Onis, Ezzati, Grantham-McGregor, Katz, Martorell and Uauy6]. Given its high burden and far-reaching consequences, comprehensive and well-informed strategies addressing stunting must remain a priority [Reference Black, Victora, Walker, Bhutta, Christian, de Onis, Ezzati, Grantham-McGregor, Katz, Martorell and Uauy6–Reference Hoddinott, Alderman, Behrman, Haddad and Horton8].
Until recently, stunting was generally understood to occur during the first 1000 days of life, with little potential for recovery in later years [Reference Onis and Branca4, Reference Black, Victora, Walker, Bhutta, Christian, de Onis, Ezzati, Grantham-McGregor, Katz, Martorell and Uauy6]. Most research and programmes to date addressing stunting and associated risk factors have thus focused on children up to 5 years of age [Reference Aguayo, Badgaiyan and Paintal9–Reference Varela-Silva, Azcorra, Dickinson, Bogin and Frisancho20]. However, a growing body of evidence from diverse populations indicates the potential to transition between stunted and non-stunted status up till at least 15 years of age [Reference Adair21–Reference Prentice, Ward, Goldberg, Jarjou, Moore, Fulford and Prentice29]. Recent research also suggests cognitive gains among children and adolescents who show catch-up growth to normal height at older ages, as well as deficits among older children and adolescents who transition to being stunted [Reference Mendez and Adair24–Reference Lundeen, Behrman, Crookston, Dearden, Engle, Georgiadis, Penny and Stein28]. More data are needed to better characterise linear growth trajectories and stunting throughout childhood, assess influences on the same, and understand subsequent later life outcomes. This includes a clearer assessment of potential differences in stunting classification by growth reference, given the implications this may have in terms of comparing estimates of stunting burden across populations and informing resource planning for strategies addressing stunting [Reference Partap, Young, Allotey, Sandhu and Reidpath30], and an exploration of potential determinants. Such evidence would facilitate a more comprehensive understanding of the risk of stunting and potential for recovery across childhood, and contribute to a clearer picture of implications for population health and productivity. We sought to address the current gap in understanding of stunting among older children and adolescents by assessing the prevalence and sociodemographic correlates of stunting in 6759 children and adolescents aged 6–19 years in Segamat, Malaysia.
Subjects and methods
Study population and measures of interest
We used data from two rounds of enumeration (2012–2013 and 2013–2014) and a cross-sectional health survey (2013–2014) conducted by the South East Asia Community Observatory health and demographic surveillance system (SEACO HDSS) operating in Segamat, Malaysia [Reference Allotey, Reidpath, Devarajan, Rajagobal, Yasin, Arunachalam, Imelda, Soyiri, Davey and Jahan31, Reference Partap, Young, Allotey, Soyiri, Jahan, Komahan, Devarajan, Sandhu and Reidpath32]. Enumerations included all consenting individuals living within the area covered by SEACO, whilst the health survey collected measures from individuals aged 6 years and above using standardised tools; information was linked across databases [33–37]. Ethical approval for all data collections was obtained from the Monash University Human Research Ethics Committee, and informed consent was obtained from all participants prior to data collection.
Analyses included all children and adolescents with information on age, sex, ethnicity and height (N = 6759). Children were defined as individuals aged <10 years, whilst adolescents were defined as individuals aged 10–19 years [38–Reference de Onis, Onyango, Borghi, Siyam, Nishida and Siekmann41]. Information from the 2012–2013 enumeration was used to match children and adolescents to their parents and determine birth order. Data from both enumerations were used to identify the child's household indices including number of household members, bedrooms, bathrooms and living areas, type of toilet, whether the toilet was shared with other households, main source of drinking water and main method of garbage disposal. All other information was obtained from the health survey. This included the child's age, sex and ethnicity, the child's and parents’ current height and weight (measured using a Transtek digital weighing scale with height gauge, model GBS-721), parents’ age, and parents’ and household head's education.
Variable transformation
We used the World Health Organization (WHO) 2007 and the Centers for Disease Control and Prevention (CDC) 2000 child growth references to express height-for-age z-score and classify stunting among children and adolescents [Reference de Onis, Onyango, Borghi, Siyam, Nishida and Siekmann41, Reference Kuczmarski, Ogden, Guo, Grummer-Strawn, Flegal, Mei, Wei, Curtin, Roche and Johnson42]. These references are distinct in terms of the reference populations chosen and methods used to derive height-for-age reference curves across childhood and adolescence, and in terms of their recommended cut-offs to define stunting [Reference de Onis, Onyango, Borghi, Siyam, Nishida and Siekmann41–Reference Borghi, De Onis, Garza, Van den Broeck, Frongillo, Grummer-Strawn, Van Buuren, Pan, Molinari and Martorell43]. Conversion of absolute height to reference-specific height was done using the Zanthro package in Stata [Reference Vidmar, Cole and Pan44]. The WHO 2007 reference classifies stunting as height <−2 standard deviations (s.d.) from the median [Reference de Onis, Onyango, Borghi, Siyam, Nishida and Siekmann41], whilst the CDC 2000 reference classifies stunting as height <5th percentile [Reference Kuczmarski, Ogden, Guo, Grummer-Strawn, Flegal, Mei, Wei, Curtin, Roche and Johnson42, 45]. Body mass index (BMI) was calculated as weight (kg) divided by the square of the height (m), and BMI-for-age status was also expressed relative to the WHO and CDC references [Reference de Onis, Onyango, Borghi, Siyam, Nishida and Siekmann41, Reference Kuczmarski, Ogden, Guo, Grummer-Strawn, Flegal, Mei, Wei, Curtin, Roche and Johnson42, 45] using the Zanthro package [Reference Vidmar, Cole and Pan44]. The WHO 2007 reference classifies underweight, normal BMI, overweight and obesity as BMI-for-age <−2 s.d., −2 to 2 s.d., >2 to 3 s.d. and >3 s.d., respectively [Reference de Onis, Onyango, Borghi, Siyam, Nishida and Siekmann41]. The CDC 2000 reference defines underweight, normal BMI, overweight and obesity as <5th percentile, 5th to 84th percentile, 85th to 94th percentile, and ⩾95th percentile, respectively [Reference Kuczmarski, Ogden, Guo, Grummer-Strawn, Flegal, Mei, Wei, Curtin, Roche and Johnson42, 45]. Maternal BMI was calculated as described above, and maternal underweight was defined as BMI <18.5 kg/m2. Parental height was assessed as a continuous variable and also categorised into five groups, as has been done previously (160 and above, 155–159, 150–154, 145–149 and <145 cm) [Reference Danaei, Andrews, Sudfeld, Fink, McCoy, Peet, Sania, Smith Fawzi, Ezzati and Fawzi12]. The number of bedrooms, bathrooms and living areas in the household was divided by the number of members to obtain the ratio of rooms to people.
Statistical analysis
We first examined crude stunting prevalence obtained with stunting classified using the WHO v. the CDC reference in the full study population (N = 6759). Prevalence was calculated overall and across demographic subgroups (sex, age and ethnicity). We used Cohen's κ statistic to assess the agreement in classification of stunting between the WHO and CDC references, overall and across demographic subgroups. κ < 0.6 was defined as poor agreement, 0.6 to <0.8 as moderate agreement, 0.8 to <0.9 as good agreement and ⩾0.9 as excellent agreement [Reference McHugh46]. We examined within-household and within-sub-district clustering of stunting in the full population using both references. We then compared adjusted estimates for stunting prevalence using the WHO and CDC references, overall and across demographic subgroups, generated using mixed-effects Poisson regression adjusted for sex, age and ethnicity, and for clustering at the household level.
We then examined associations between stunting and a range of sociodemographic indices in a subset of the population with complete information on all measures of interest, including individual (age, sex, ethnicity, birth order and BMI-for-age status), maternal (height and current underweight) and household measures (bedrooms, bathrooms and living areas per household member, type of toilet, sharing of toilet with other households, main source of drinking water and main method of garbage disposal) (N = 3791). Initially, we assessed differences in measures of interest between children and adolescents included in regressions v. those who were excluded. Following this, we examined sociodemographic measures of interest across categories of stunting among included children and adolescents. Differences in distributions of measures by inclusion status or by stunting status were assessed using Pearson's χ2 tests for categorical variables (Fisher's exact tests for cell counts less than five) and Student's t-test for continuous variables. We then examined univariable associations between each measure and stunting (using Poisson regression). The linearity of continuous measures was assessed by comparing models with continuous v. categorised measures using likelihood ratio tests.
Following this, we used mixed-effects Poisson regression models with robust standard errors to assess the relation between sociodemographic indices and stunting expressed using the CDC reference. Regression models were built with progressive adjustment for groups of variables, starting with individual measures, and then additionally including maternal and then household measures. Fully-adjusted models included all individual, maternal and household variables described above, and were additionally adjusted for clustering at the household level. To confirm trends and obtain more concise estimates of risk, we additionally examined models (i) including maternal height as a continuous variable, and (ii) including variables with collapsed categories for maternal height, type of toilet and main source of drinking water.
Effect modification of stunting risk by sex was explored by comparing minimally-adjusted mixed-effects Poisson regression models (age, ethnicity, sex, measure of interest) with those including an additional term for interaction (sex × measure of interest) using likelihood ratio tests. Furthermore, we performed a number of sensitivity analyses to assess the robustness of associations observed. Crude and fully adjusted associations between maternal age and stunting were assessed, and the effect of its inclusion on other associations was examined. We similarly examined education of the mother, father and head of household, and paternal height. We assessed the specificity of associations by checking the effect of inclusion of the following orthogonal variables in regression models: maternal heart rate, maternal diastolic blood pressure, availability of Internet from the provider Streamyx in the household, and the ratio of motorcycles in the household to the number of household members. Analysis of complete records was undertaken for these sensitivity analyses.
Additionally, we examined the relation between stunting risk and the co-occurrence of associated sociodemographic indices identified from the preceding regression models. These indices were binarised to facilitate clearer interpretation, and included Malay ethnicity, child underweight, maternal height <145 cm, living in a household with an unprotected water source, living in a household with a bucket, hanging or no latrine, and living in a household with a shared toilet. Mixed-effects Poisson regression with robust standard errors was undertaken with stunting as the outcome and the number of associated indices as the primary exposure. The final model was adjusted for all other sociodemographic indices of interest, and for clustering at the household level.
Finally, we checked associations observed using mixed-effects logistic regression models with robust standard errors, and with mixed-effects linear regression models; final models were built identically to Poisson regression models above. We also ran all analyses using the WHO reference to express height-for-age and classify stunting, to assess comparability in the direction and effect size of estimates of association.
All analyses were performed using Stata 14 (Statacorp, Texas, USA).
Results
Analyses covered a total of 6759 children and adolescents (49.3% boys). Children aged 6–9 years comprised 26.2% of the population, and approximately 37% each were aged 10–14 and 15–19 years. The majority of children and adolescents were of Malay ethnicity, followed by Chinese and then Indian ethnicity (67.3, 19.9 and 9.8%, respectively). There was no difference in age or ethnicity distributions by sex (Table 1).
Table 1. Demographic characteristics of study population
Differences in distributions across categories between girls and boys were compared using Pearson's χ2 test.
We first assessed the agreement between the CDC and WHO references in stunting classification and subsequent estimation of prevalence. Stunting prevalence was notably higher when using the CDC reference compared with the WHO reference, with only moderate agreement between the two (crude prevalence using WHO reference: 19.1%, using CDC reference: 29.2%, κ for agreement: 0.73) (Table 2, online Supplementary Table S1). Differences in prevalence estimates, and demographic patterns in prevalence, were generally consistent between references (Table 2, online Supplementary Table S1). Stunting was strongly clustered within households, regardless of the reference used (online Supplementary Table S2).
Table 2. Prevalence of stunting amongst children in Segamat, Malaysia according to the World Health Organization 2007 v. Centers for Disease Control 2000 reference
WHO, World Health Organization; CDC, Centers for Disease Control.
Estimates are for adjusted prevalence, based on mixed-effects Poisson regression models including sex, age and ethnicity and adjusted for clustering at the household level.
a Agreement in classification of stunting between references was calculated using Cohen's κ coefficient.
We further explored the relation between stunting and a range of sociodemographic indices in a subset of the population with complete information on all measures (N = 3791). Children and adolescents included in this subset were generally comparable to those who were excluded in terms of measures of interest (online Supplementary Table S3). Similar to the full study population, 28.5% of children and adolescents were stunted (CDC reference) (Table 3). Compared with children and adolescents of normal height, stunted children and adolescents were marginally older, and were less commonly of Chinese or Indian ethnicity, overweight or obese, of birth order four or above, and with mothers of height 150 cm or above (p < 0.001 for all). Stunted children and adolescents were also less likely to be living in a household with a flush toilet or a toilet not shared with other households (p = 0.009 and p = 0.025, respectively), or a household where the main source of drinking water was piped into the yard or house (p < 0.001) (Table 3).
Table 3. Sociodemographic characteristics of study population across categories of stunting
BMI, body mass index.
Descriptive analyses presented for a subset of children with complete information on all sociodemographic variables of interest above (N = 3791), whose data were used for regression analyses (Table 4, Figs 1–3).
Differences in distributions across categories between normal height and stunted were compared using Pearson's χ2 test, or Fisher's exact test for variables with cell frequencies <5.
a Stunting and BMI-for-age were classified using the Centers for Disease Control 2000 reference.
b Living areas include dining rooms but not kitchens.
c Bore hole toilet both with or without cover.
d Unprotected sources: unprotected dug well, or water taken directly from pond or stream.
e Protected sources: protected dug well or spring, or water from bottles or tanker truck.
We examined associations between sociodemographic indices and child stunting (CDC reference) using mixed-effects Poisson regression. In fully adjusted analyses, stunting risk increased by 8% per year increase in age [fully adjusted risk ratio (RR): 1.08, 95% confidence interval (95% CI) 1.06–1.09, p < 0.001). Children and adolescents of Indian or Chinese ethnicity had an approximately 30% decreased risk of stunting compared with those of Malay ethnicity (RR for Indian ethnicity: 0.66, 95% CI 0.52–0.83, p < 0.001; for Chinese ethnicity: 0.68, 95% CI 0.55–0.83, p < 0.001) (Table 4, Fig. 1, online Supplementary Table S4). Overweight children and adolescents also had a reduced risk of stunting compared with those of normal BMI, whilst underweight was positively associated with stunting risk (RR for overweight: 0.79, 95% CI 0.69–0.90, p = 0.001; for underweight: 1.19, 95% CI 1.02–1.39, p = 0.030). Children and adolescents with mothers of height <145 cm had a 53% increased risk of stunting compared with those whose mothers were 160 cm or taller (RR 1.53, 95% CI 1.21–1.92, p < 0.001). We found no evidence of associations between sex, birth order or current maternal underweight and child stunting risk (Table 4).
Table 4. Relative risk of stuntinga associated with sociodemographic indices
BMI, body mass index.
Estimates based on mixed-effects Poisson regression models adjusted for all other variables above, and for clustering at the household level.
a Stunting and BMI-for-age were classified using the Centers for Disease Control 2000 reference.
b Living areas include dining rooms but not kitchens.
c Bore hole toilet both with or without cover.
d Unprotected sources: unprotected dug well, or water taken directly from pond or stream.
e Protected sources: protected dug well or spring, or water from bottles or tanker truck.
Fig. 1. Association of stunting with selected individual, parental and household water and sanitation indices – stunting expressed using the (a) Centers for Disease Control and Prevention 2000 reference and (b) World Health Organization 2007 reference. Estimates presented from Poisson regression models described in online Supplementary Tables S5 and S10. RR, risk ratio; 95% CI, 95% confidence interval (error bars). Flush toilet includes that connected with septic tank or sewerage system; unprotected drinking water sources: unprotected dug well, or water taken directly from pond or stream; protected sources: public standpipe, protected dug well or spring, or water from bottles or tanker truck.
Additionally, children and adolescents whose main source of drinking water was pumped into the house or yard had an approximately 70% reduced risk of stunting v. those with access to unprotected sources, with smaller inverse associations for other protected water sources (RR for water piped into yard: 0.25, 95% CI 0.15–0.39, p < 0.001; for other protected sources: 0.52, 95% CI 0.38–0.73, p < 0.001). Furthermore, children and adolescents with a flush toilet in the household had an approximately 60% decreased risk of stunting v. those with a bucket, hanging or no latrine (RR for flush toilet with septic tank: 0.38, 95% CI 0.18–0.84, p = 0.017) (Table 4, Fig. 1, online Supplementary Table S4). Consistent trends were observed for related sanitation indices (RR for shared household toilets v. not shared: 1.72, 95% CI 1.07–2.74, p = 0.025; RR per unit increase in bathrooms per household member: 0.67, 95% CI 0.46–0.98, p = 0.036, respectively) (Table 4). We found no evidence of effect modification of associations between measures of interest and stunting risk by sex (online Supplementary Table S6).
We confirmed the robustness and specificity of observed associations through multiple additional checks and sensitivity analyses. Similar magnitudes and directions of association were observed when analyses were repeated with height-for-age and stunting expressed according to the WHO reference (Fig. 1, online Supplementary Table S7). The additional inclusion of maternal age, paternal height, or parental or household head's education to fully adjusted models did not materially change observed associations (Fig. 2, online Supplementary Tables S8–S11), with no evidence of association between parental or head of household's education and stunting risk (online Supplementary Table S9). Furthermore, we observed no association between stunting and theoretically uncorrelated maternal and household variables in fully adjusted regression models, with no effect of their inclusion on other associations (Fig. 2, online Supplementary Table S12). Trends in associations when using linear or logistic regression were consistent with those observed for Poisson regression models (online Supplementary Table S5, data not shown for logistic regression).
Fig. 2. Summary of effect of sensitivity analyses on associations between key indices and stunting risk. (a) Original regression model (online Supplementary Table S5; N = 3791); (b) model with maternal age included (N = 3337); (c) model with maternal education included (N = 3625); (d) model with paternal height included (N = 2723); (e) model with orthogonal variables included (N = 3434). RR, risk ratio; 95% CI, 95% confidence interval (error bars); BMI, body mass index. Stunting and BMI-for-age were defined using the Centers for Disease Control and Prevention 2000 reference. Flush toilet includes that connected with septic tank or sewerage system; unprotected drinking water sources: unprotected dug well, or water taken directly from pond or stream; protected sources: public standpipe, protected dug well or spring, or water from bottles or tanker truck. Orthogonal variables: maternal heart rate, maternal diastolic blood pressure, household Streamyx Internet and motorcycles per household member.
We further examined the relation between stunting and the co-occurrence of associated sociodemographic indices identified in the previous regressions. Indices included Malay ethnicity, child underweight, maternal height <145 cm, unprotected drinking water source, none or hanging or bucket latrine in the household, and shared household toilet. Thirty-four (0.9%) children and adolescents had a maximum of three of any of these measures (online Supplementary Table S13). The risk of stunting increased consistently with the number of co-occurring indices, with a greater than twofold increased risk among children and adolescents with three indices v. none (RR 2.48, 95% CI 1.51–4.09, p < 0.001) (Fig. 3, online Supplementary Table S14).
Fig. 3. Relative risk of stunting associated with the number of co-occurring risk factors. Stunting was classified using the Centers for Disease Control and Prevention 2000 reference. Dashed line represents null association (relative risk of 1); error bars denote 95% confidence interval. Risk factors included: Malay ethnicity, child underweight, maternal height <145 cm, living in a household with unprotected drinking water source (unprotected dug well, or water taken directly from pond or stream), living in household with no bucket or hanging latrine and living in a household with a shared toilet. Estimates are based on Poisson regression models adjusted for the child's age, sex and body mass index-for-age status, birth order, maternal current underweight, rooms per household member, household's main method of garbage disposal and for clustering at the household level.
Discussion
In this large population of Malaysian children and adolescents, we observed a high burden of stunting, with notable discrepancies in prevalence estimated using two international height references. Stunting was independently associated with indices related to nutritional deficits and poor household environment; associations were robust to a series of sensitivity analyses and to the height reference used. The risk of stunting increased with the number of co-occurring indices. Our results provide strong support for a framework for strategies addressing stunting across childhood, and highlight the need for consensus on a single definition of stunting among older children and adolescents to streamline monitoring efforts.
Evidence on the burden and determinants of stunting among children and adolescents aged over 5 years is limited [Reference Svefors, Rahman, Ekstrom, Khan, Lindstrom, Persson and Selling K22, Reference Fetuga, Ogunlesi, Adekanmbi and Alabi47–Reference Senbanjo, Oshikoya, Odusanya and Njokanma52]. This includes research on the comparability of height references in the classification of stunting and subsequent estimation of prevalence. We identified only three other studies comparing the WHO and CDC references among older children and adolescents, based on smaller populations of narrower age ranges than in this study [Reference Fetuga, Ogunlesi, Adekanmbi and Alabi47–Reference Mushtaq, Gull, Mushtaq, Abdullah, Khurshid, Shahid, Shad and Akram49]. All reported discrepancies in stunting classification by reference, although smaller than those observed here [Reference Fetuga, Ogunlesi, Adekanmbi and Alabi47–Reference Mushtaq, Gull, Mushtaq, Abdullah, Khurshid, Shahid, Shad and Akram49]. The large differences in prevalence estimates by reference in our study population have important implications for decisions relating to targeting of interventions and health resource planning. Furthermore, our results indicate a fundamental gap in the current understanding of normal and compromised linear growth at older ages, including whether this is distinct to population. Further longitudinal research is required in this area, along with a more immediate consensus on the universal use of one reference to harmonise monitoring of stunting in this older age group across populations.
The associations between stunting and environmental and nutritional indices observed in our study are consistent with previous work, which has mainly been undertaken in younger children [Reference Aguayo, Badgaiyan and Paintal9–Reference Varela-Silva, Azcorra, Dickinson, Bogin and Frisancho20]. Evidence from older children and adolescents is scarce, generally based on pre-adolescents and on smaller populations than examined here, with limited exploration of environmental or nutritional measures [Reference Svefors, Rahman, Ekstrom, Khan, Lindstrom, Persson and Selling K22, Reference Dearden, Schott, Crookston, Humphries, Penny and Behrman50–Reference Senbanjo, Oshikoya, Odusanya and Njokanma52]. Among younger children, poor sanitation and drinking water facilities are understood to contribute to stunting through increased risk of acute and chronic enteric infections, leading to malabsorption of nutrients essential to growth [Reference Black, Victora, Walker, Bhutta, Christian, de Onis, Ezzati, Grantham-McGregor, Katz, Martorell and Uauy6, Reference Cumming and Cairncross53]. Measures such as underweight or wasting are established indicators of undernutrition resulting from insufficient nutrient intake or absorption, and are thought to at least partly share biological pathways leading to stunting [Reference Black, Victora, Walker, Bhutta, Christian, de Onis, Ezzati, Grantham-McGregor, Katz, Martorell and Uauy6, Reference Khara and Dolan54]. Our results suggest that poor environment and nutrition are also important determinants of linear growth and stunting in later childhood. Importantly, our study indicates that these indices contribute to stunting risk independently from the potential genetic, intrauterine and shared environmental influences of both parents’ height [Reference Ozaltin, Hill and Subramanian55–Reference Garza, Borghi, Onyango and de Onis62]. Targeting these modifiable factors throughout childhood may thus be a useful approach to improving linear growth in this population, with potential for further subsequent intergenerational benefits. On the other hand, the implications of the increased risk of stunting observed among Malay children and adolescents are unclear. There has been debate regarding whether ethnicity-specific differences observed in child growth beyond 5 years of age may be a result of modifiable environmental exposures or genetic differences in growth potential [Reference de Onis, Onyango, Borghi, Siyam, Nishida and Siekmann41]. Whilst we accounted for potential differences in household environment and BMI in this study, we were unable to account for other ethnicity-specific factors such as diet that may additionally explain such differences in risk. Our observations point to an important area of inquiry that may be best addressed through additional studies, including qualitative or genetic research.
Certain observations in our study stand somewhat in contrast to the totality of previous evidence on determinants of stunting risk. Measures of overcrowding and family size [Reference Rakotomanana, Gates, Hildebrand and Stoecker17, Reference Mushtaq, Gull, Khurshid, Shahid, Shad and Siddiqui51, Reference Senbanjo, Oshikoya, Odusanya and Njokanma52], maternal age and underweight [Reference Danaei, Andrews, Sudfeld, Fink, McCoy, Peet, Sania, Smith Fawzi, Ezzati and Fawzi12, Reference Rachmi, Agho, Li and Baur15, Reference Saxton, Rath, Nair, Gope, Mahapatra, Tripathy and Prost18, Reference Black, Allen, Bhutta, Caulfield, de Onis, Ezzati, Mathers and Rivera56], and parental education [Reference Black, Victora, Walker, Bhutta, Christian, de Onis, Ezzati, Grantham-McGregor, Katz, Martorell and Uauy6, Reference Aguayo, Nair, Badgaiyan and Krishna10, Reference Akombi, Agho, Hall, Merom, Astell-Burt and Renzaho11, Reference Fregonese, Siekmans, Kouanda, Druetz, Ly, Diabate and Haddad14–Reference Rah, Cronin, Badgaiyan, Aguayo, Coates and Ahmed16, Reference Torlesse, Cronin, Sebayang and Nandy19] or related measures of socioeconomic status [Reference Zhang, Becares and Chandola63, Reference Kien, Lee, Nam, Oh, Giang and Van Minh64] have previously been implicated in poor linear growth, but were not associated with stunting here. As has been observed with risk factors in other studies, it may be that these measures are not important determinants of stunting in our study population [Reference Aguayo, Badgaiyan and Paintal9, Reference Aguayo, Nair, Badgaiyan and Krishna10, Reference Rachmi, Agho, Li and Baur15, Reference Dearden, Schott, Crookston, Humphries, Penny and Behrman50]. Importantly, we observed increased stunting risk with a low bathroom to household member ratio, and with sharing of toilets between households, suggesting that overcrowding may be relevant in more extreme contexts in this population. As such, these results reiterate the importance of population-specific research to inform suitable strategies addressing stunting.
Valuable insights for the design and targeting of interventions may be gained by examining the co-occurrence of risk factors and their collective effect on stunting risk. Yet, to our knowledge, this has not been explored in earlier studies. In our study population, out of six measures associated with stunting, a small group of children and adolescents had the maximum of three co-occurring indices, and stunting risk increased generally consistently with the number of co-occurring indices. Our population size limited our ability to assess in greater detail the risk of stunting associated with specific combinations of indices, which could be more informative for the purposes of intervention targeting. Nonetheless, we highlight here a potentially useful approach to identify sub-populations at greater risk of stunting through examining the extent and patterns of distribution of co-occurring risk factors.
Child undernutrition is being increasingly recognised as a public health issue of concern in Malaysia, with recent estimates suggesting a notable and rising burden of stunting among children under 5 years of age [65]. Both governmental and non-governmental organisations have acknowledged the growing need to implement targeted strategies to address undernutrition among children in Malaysia [65–67]. Importantly, these would have to be undertaken in coordination with interventions addressing the simultaneously increasing burden of overweight and obesity among children [66, Reference Cheah, Wan Muda, Mohd Hussin and Thon68–72]. Given the high overweight and obesity prevalence observed among children and adolescents in our study population, our results provide further support for the need for a continued research and programmatic focus on the double burden of malnutrition throughout childhood in Malaysia.
To our knowledge, ours is among the largest of the few studies to date that have examined the prevalence and determinants of stunting in children and adolescents older than 5 years of age [Reference Svefors, Rahman, Ekstrom, Khan, Lindstrom, Persson and Selling K22, Reference Fetuga, Ogunlesi, Adekanmbi and Alabi47–Reference Senbanjo, Oshikoya, Odusanya and Njokanma52]. Our study provides strong evidence for associations between environmental and nutritional indices and stunting risk: observed associations were specific, robust to adjustment for multiple covariates and independent of the height reference used to express stunting. For household sanitation, we assessed three measures capturing distinct aspects of the availability and quality of facilities, all of which were independently and consistently related to stunting risk. As this is a cross-sectional analysis, we are unable to make definitive statements regarding causality, and to more clearly distinguish between the relative contribution to stunting risk of current v. historical nutritional and environmental exposures, which may be correlated [Reference Dearden, Schott, Crookston, Humphries, Penny and Behrman50]. Furthermore, we were unable to consider more detailed measures capturing dimensions relating to nutrition, such as diet or food security. However, in light of the previous research suggesting the potential for transition between stunted and non-stunted status at older ages [Reference Adair21–Reference Prentice, Ward, Goldberg, Jarjou, Moore, Fulford and Prentice29], we believe our study provides key information regarding potential factors that may inform such transitions. Importantly, our study indicates that these are no different from factors influencing linear growth at younger ages, suggesting that broadening the current conceptual view of stunting to one that extends across childhood may ultimately allow for greater public health gains.
To conclude, we report independent associations between stunting and nutritional and household environmental indices in this population of older children and adolescents from Malaysia, with notable increases in risk associated with the co-occurrence of indices. Importantly, we observed a high burden of stunting in this population, with large differences in stunting classification by height reference. In the context of the increasing prevalence of overweight and obesity in this and similar populations [Reference Ng, Fleming, Robinson, Thomson, Graetz, Margono, Mullany, Biryukov, Abbafati, Abera, Abraham, Abu-Rmeileh, Achoki, AlBuhairan, Alemu, Alfonso, Ali, Ali, Guzman, Ammar, Anwari, Banerjee, Barquera, Basu, Bennett, Bhutta, Blore, Cabral, Nonato, Chang, Chowdhury, Courville, Criqui, Cundiff, Dabhadkar, Dandona, Davis, Dayama, Dharmaratne, Ding, Durrani, Esteghamati, Farzadfar, Fay, Feigin, Flaxman, Forouzanfar, Goto, Green, Gupta, Hafezi-Nejad, Hankey, Harewood, Havmoeller, Hay, Hernandez, Husseini, Idrisov, Ikeda, Islami, Jahangir, Jassal, Jee, Jeffreys, Jonas, Kabagambe, Khalifa, Kengne, Khader, Khang, Kim, Kimokoti, Kinge, Kokubo, Kosen, Kwan, Lai, Leinsalu, Li, Liang, Liu, Logroscino, Lotufo, Lu, Ma, Mainoo, Mensah, Merriman, Mokdad, Moschandreas, Naghavi, Naheed, Nand, Narayan, Nelson, Neuhouser, Nisar, Ohkubo, Oti, Pedroza, Prabhakaran, Roy, Sampson, Seo, Sepanlou, Shibuya, Shiri, Shiue, Singh, Singh, Skirbekk, Stapelberg, Sturua, Sykes, Tobias, Tran, Trasande, Toyoshima, van de Vijver, Vasankari, Veerman, Velasquez-Melendez, Vlassov, Vollset, Vos, Wang, Wang, Weiderpass, Werdecker, Wright, Yang, Yatsuya, Yoon, Yoon, Zhao, Zhou, Zhu, Lopez, Murray and Gakidou71, Reference Pell, Allotey, Evans, Hardon, Imelda, Soyiri and Reidpath73], our data on stunting point to the complex impact of malnutrition. Our study supports a framework for stunting across childhood, and highlights the need for consensus on a single height reference to monitor stunting at older ages. It provides a basis for further longitudinal research on the exact definition, determinants and consequences of linear growth and stunting throughout childhood in diverse populations. Such work will be fundamental to the design of strategies that may more comprehensively and effectively address childhood stunting worldwide.
Author ORCIDs
Uttara Partap, 0000-0002-2531-1804; Pascale Allotey, 0000-0002-6942-5774; Daniel D. Reidpath, 0000-0002-8796-0420
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/gheg.2019.1.
Financial support
SEACO is funded by the office of the Vice Provost Research, Monash University Australia; the office of the Deputy Dean Research, Faculty of Medicine, Nursing and Health Sciences, Monash University Australia; the Monash University Malaysia Campus and the Jeffrey Cheah School of Medicine and Health Sciences. SEACO is an associate member of the INDEPTH Network. This work was supported by the Wellcome Trust (grant number 098051). MS is supported by the National Institute for Health Research Cambridge Biomedical Research Centre (UK). The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care. UP is supported by the Dr Herchel Smith Fellowship.
Conflict of interest
None.
Ethical standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.
Open access
