Importance: Irrespective of their genetic makeup, children living in an ideal home environment that supports healthy growth have similar growth potential. However, whether this potential is true for children residing at higher altitudes remains unknown. Objective: To investigate whether altitude is associated with increased risk of linear growth faltering and evaluate the implications associated with the use of the 2006 World Health Organization growth standards, which have not been validated for populations residing 1500 m above sea level. Design, Settings, and Participants: Analysis of 133 nationally representative demographic and health cross-sectional surveys administered in 59 low- A nd middle-income countries using local polynomial and multivariate regression was conducted. A total of 964299 height records from 96552 clusters at altitudes ranging from-372 to 5951 m above sea level were included. Demographic and Health Surveys were conducted between 1992 and 2018. Exposures: Residence at higher altitudes, above and below 1500 m above sea level, and in ideal home environments (eg, access to safe water, sanitation, and health care). Main Outcomes and Measures: The primary outcome was child linear growth deficits expressed in length-for-age/height-for-age z scores (HAZ). Associations between altitude and height among all children and those residing in ideal home environments were assessed. Child growth trajectories above and below 1500 m above sea level were compared and the altitude-mediated height deficits were estimated using multivariable linear regression. Results: In 2010, a total of 842 million people in the global population (approximately 12%) lived 1500 m above sea level or higher, with 67% in Asia and Africa. Eleven percent of the sample was children who resided 1500 m above sea level or higher. These children were born at shorter length and remained on a lower growth trajectory than children residing in areas less than 1500 m above sea level. The negative association between altitude and HAZ was approximately linear through most part of the altitude distribution, indicating no clear threshold for an abrupt decrease in HAZ. A 1000-m above sea level increase in altitude was associated with a 0.163-unit (95% CI,-0.205 to-0.120 units) decrease in HAZ after adjusting for common risk factors using multivariable linear regressions. The HAZ distribution of children residing in ideal home environments was similar to the 2006 World Health Organization HAZ distribution, but only up to 500 m above sea level. Conclusions and Relevance: The findings of this study suggest that residing at a higher altitude may be associated with child growth slowing even for children living in ideal home environments. Interventions addressing altitude-mediated growth restrictions during pregnancy and early childhood should be identified and implemented..
While researchers have used different definitions for high altitude,16 we used the 1500 m above sea level threshold because sites located above this altitude were not considered eligible in the Multicentre Growth Reference Study. The analysis took place in stages. Surveys were conducted between 1992 and 2018. We first estimated the number of people residing 1500 m above sea level in each country. We then compared child growth trajectories at lower than 1500 m above sea level and 1500 m or more above sea level altitudes using local polynomial regression methods. Previous literature has attributed high-altitude growth deficits to poorer nutrition, health, and socioeconomic conditions at high-altitude localities.13,17 Considering these issues, we accounted for the role of confounding factors in 3 ways. First, we used linear regression methods to quantify the child height deficit associated with altitude after adjusting for confounding factors. Second, we used the same multivariable regression models to assess how altitude is associated with immediate causes of malnutrition: diets and disease.4 Third, we restricted the sample to children who resided in ideal home environments8 and used regression methods to assess the association between altitude and child height within this sample. The protocols and questionnaires of DHS surveys have been reviewed and approved by ICF institutional review board and the institutional review boards of the host countries. This study followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guideline for cross-sectional studies. The country-level population data disaggregated by altitude were based on data from the Center for International Earth Science Information Network (CIESIN), Columbia University.18 We used nationally representative and publicly available cross-sectional DHS survey data from 59 low- and middle-income countries. We used all available DHS surveys that collected anthropometric measures for children younger than 5 years and restricted the analysis to surveys that recorded the altitude of the survey cluster or the global positioning system coordinates of the cluster that permitted us to obtain the altitude from external sources. A cluster was equivalent to a village in rural areas or a neighborhood in urban areas. Linear growth (faltering) was expressed in HAZ that measures the distance in height to the median child in the 2006 WHO growth standard.7 For convenience, we used the term HAZ to refer to both length for age (<24 months) and height for age (≥24 months). Children with biologically implausible height measures (HAZ 6) were excluded from the analysis. A total of 964 299 height records (51% boys, 49% girls) from 96 552 clusters from 133 surveys administered in 59 countries were used (eTable 1 in the Supplement). The clusters’ altitudes ranged from −372 to 5951 m above sea level. To examine HAZ age trajectories in low- and middle-income countries,19 we used local polynomial regressions that regressed the child’s HAZ on their age in months. We ran separate regressions for children residing in clusters lower than 1500 m above sea level (n = 857 858) and children residing in clusters 1500 m or more above sea level (n = 106 441). We compared the regression lines and their 95% CIs between the 2 groups to assess whether the differences in growth trajectories were statistically different from zero. The weights in these regressions were based on the Epanechnikov kernel-density function, and the bandwidth was selected using the rule-of-thumb method.20 The HAZ of children was then regressed on the altitude of the DHS cluster in which the child resided using a local polynomial regression, where the weights and the bandwidth were selected using the methods described above. Unadjusted and adjusted linear regression models were used to quantify the HAZ deficit per increments of 1000 m above sea level in the altitude of the DHS cluster in which the child resided (eFigure 1 in the Supplement). Adjusted models controlled for differences in biological and underlying causes of linear growth faltering4,19,21: child age (set of binary variables for each month) and sex, maternal age and level of education (years), household wealth (access to electricity and ownership of radio, television, refrigerator, bicycle, motorcycle, car, and improved floor material), and binary variables capturing access to improved water and sanitation (eTable 2 in the Supplement). We also controlled for residence in a rural area and used subnational region (highest administrative unit in each country) fixed effects to control for economic, political, climatic, and other factors shared by the residents in the same administrative area. We had 1348 survey-specific subnational regions in the full sample. Both unadjusted and adjusted regressions were based on the full sample of children aged 0 to 59 months as well as subsamples consisting of different age groups: 0 to 5, 6 to 11, 12 to 23, and 24 to 59 months. In adjusted models, children with missing control variable values were omitted from the sample. We clustered our SEs at the subnational region level.22 We used multivariable regression models to assess how altitude was associated with immediate causes of undernutrition: inadequate dietary intake and disease.4 As measures of inadequate dietary intake, we used prevalence of exclusive breastfeeding (children aged 0-6 months), dietary diversity (age 6-23 months), and minimum acceptable diet (age 6-23 months) as dependent variables. For disease risk, we used incidence of diarrhea, fever, or cough in the 2 weeks preceding the survey (age 0-59 months). More details are provided in eAppendix 1 and eTable 3 in the Supplement. The adjusted logistic regressions were based on the same set of control variables as described above and the SEs were clustered at the subnational region level. The estimated HAZ deficits per 1000-m above sea level increments in altitude are expressed as odds ratios (ORs). In the final part of the analysis, we restricted the sample to children residing in ideal home environments. Using the Karra et al guidelines,8 children were classified as having lived in an ideal home environment based on the following criteria: (1) singleton birth; (2) access to safe water and sanitation; (3) living in a house with finished floors, parents owned a television and a car; (4) born to highly educated mothers (>13 years of schooling); (5) born in a hospital; and (6) received BCG and first diphtheria, pertussis, and tetanus vaccinations (eAppendix 2 in the Supplement). Unadjusted and adjusted linear models were used to regress the HAZ of children living in ideal environments on the altitude of the cluster in which the child resided. Adjusted models controlled for child sex and age. We used heteroskedasticity robust SEs with this ideal home environment sample. We ran multiple robustness checks to assess the sensitivity of our results. We assessed whether the multivariable regression results were robust to replacing the continuous altitude variable with a binary variable capturing 1500-m above sea level clusters. Climatic patterns differ in low- and high-altitude locations and may also directly affect linear growth. While our subnational fixed effects controlled for climatic differences between regions, we augmented the model with data on long-term average rainfall and temperature in the cluster. We explored sensitivity by adding maternal height as a control in the multivariable regression model. We restricted the sample to children whose mother had lived in the same cluster at least since conception. There were 18 countries that had no clusters in areas 1500 m or more above sea level. We assessed sensitivity by omitting these 18 countries from the sample. We explored sensitivity of the regression results by omitting each country from the sample. We assessed the association between altitude and stunting (HAZ <−2) using the same regression methods as described above. Stata, version 16.1 (StataCorp LLC) was used for all analyses.