Abstract

Background and Aim. Overweight prevalence in children and adolescents shows great variability which is related to individual-level and environmental-level factors. The present study aimed to determine the prevalence of and factors associated with overweight in Peruvian children and adolescents living at different altitudes. Methods. 8568 subjects, aged 6–16 y, from the sea level, Amazon, and high-altitude regions were sampled. Overweight was identified using BMI; biological maturation and physical fitness were measured; school characteristics were assessed via an objective audit. Results. Overweight prevalence decreased with age (28.3% at 6 y to 13.9% at 16 y); it was higher in girls (21.7%) than boys (19.8%) and was higher at the sea level (41.3%), compared with Amazon (18.8%) and high-altitude (6.3%) regions. Approximately 79% of the variance in overweight was explained by child-level characteristics. In Model 1, all child-level predictors were significant (); in Model 2, six out of nine added school-level predictors (number of students, existence of policies and practices for physical activity, multisports-roofed, duration of Physical Education classes, and extracurricular activities) were significant (); in Model 3, subjects living at high altitudes were less likely to be overweight than those living at the sea level. Conclusions. Child- and school-level variables played important roles in explaining overweight variation. This information should be taken into account when designing more efficient strategies to combat the overweight and obesity epidemic.

1. Introduction

The prevalence of childhood overweight and obesity has systematically increased in developed and developing countries over the last few decades [1]. It is considered a global epidemic and a public health crisis [2], given the adverse health consequences throughout the life course [3] as well as the related economic costs placing an undue strain on healthcare systems [4].

A systematic analysis [1] of available data from 1980 to 2013 showed that the combined prevalence of overweight and obesity has substantially increased in children and adolescents in developed countries from 16.9% to 23.8% in boys and from 16.2% to 22.6% in girls; in developing countries, this increase was from 8.1% to 12.9% in boys and from 8.4% to 13.4% in girls. This prevalence has also increased rapidly in low- and middle-income countries [5]. Recently, Yang et al. [6] reported overall prevalences of 13.4% for underweight and 21% for overweight (including obesity) among young adolescents in 58 low- and middle-income countries.

It has been suggested that increases in childhood overweight have been steeper in Latin American countries [7]. Currently, more than 20% (approximately 42.5 million) of Latin American children aged 0 to 19 yrs are overweight or obese [8], and Peru is no exception to the growing problem of childhood undernutrition and overnutrition [9]. For example, Pajuelo-Ramírez et al. [10] investigated the prevalence of overweight (including obesity) and chronic malnutrition in 6- to 9-year-old Peruvian children and showed values of 21.5% and 17.8%, respectively.

The development of overweight and obesity is multifaceted, involving biological and environmental factors [11]. Furthermore, girls tend to be more overweight and obese than boys [12]; overweight tends to be higher in older age-groups [13], and those ahead in their maturity status are more likely to be overweight and/or obese [14]. High levels of physical fitness have been shown to be an important contributing factor in maintaining healthier body weights during childhood [15].

It is likely that context-specific information is important to explain why children vary in their weight status [16], and this is apparently evident in Peru. Its territory spreads across three main areas (sea level, Amazon region, and high altitude) with their multifaceted geographical and socioeconomic features. Furthermore, Peru is facing different stages of nutritional [17] and epidemiological [18] transitions. Available information [19] shows differences in overweight and obesity prevalences in children living in different regions, with a very high prevalence in Lima (the capital city). Tarqui-Mamani et al. [20], using data from children aged 5–13 yrs, concluded that living in the urban areas and in Metropolitan Lima was related to a greater likelihood of being overweight than living in the Amazon region or at high altitudes. These results were consistent to those reported by Álvarez-Dongo et al. [21] in children and adolescents aged from 5 to 19 yrs of age.

Children and adolescents spend most of their awake time at school, and identifying school-related factors that may be associated to overweight/obesity is of importance. It has recently been suggested that the school context has links with child/adolescent overweight and obesity [22]. For example, O’Malley et al. [23] showed that USA schools with a high concentration of students from low socioeconomic status households have a greater chance of having higher proportions of obesity students. On the other hand, Pallan et al. [24] reported that, in United Kingdom schools, the only school-level variable associated with BMI z-score was time spent in Physical Education classes (minutes/week).

Very few studies have investigated the combined links of individual- and environmental-level correlates of childhood and adolescent overweight, and this information is not available in Peruvian youth. Furthermore, using a multilevel approach to tackle this complex issue could provide new insights for the development of intervention programs to prevent childhood overweight and obesity. Therefore, the present study has the following aims: (1) to determine the prevalence of overweight in Peruvian children and adolescents by age, sex, and geographical area of residence; (2) to examine the importance of biological characteristics (age, sex, maturation, and physical fitness total score), school-level contexts, and geographical area of residence in explaining variation in BMI categories.

2. Materials and Methods

2.1. Communities

Geographical heterogeneity in Peru is largely expressed in three areas with different altitudes: (1) sea level: in the coast, Barranco (58 m) is one of the 43 districts of Lima Province and is located on the shore of the Pacific Ocean; (2) Amazon region: in the jungle area, La Merced and San Ramon (751 m) are districts with geographical continuity and integrate the Chanchamayo province; (3) high altitude: Junín is localized on a plateau at 4107 m in the central part of the Andes mountain. Geographic, demographic, socioeconomic, and educational characteristics of these areas [25] are reported in Table 1.

2.2. Study Participants

The present cross-sectional sample comes from the Peruvian Health and Optimist Growth Study which, in summary, scrutinizes the main sources of variation in physical growth, motor development, and health in Peruvian children and adolescents living at different altitudes [26]. For this paper, we used a sample of 8568 children and adolescents (3914 boys and 4654 girls, aged 6–16 yrs) from 31 public schools located at the sea level, in the Amazon region, and at high altitude. The analytical sample was selected based on the availability of complete data of anthropometry, biological maturation, and three physical fitness tests (standing long jump, shuttle run, and 12 minute run test); furthermore, we did not consider subjects classified as underweight (ntotal = 384; 4.5% ot the overall sample) because their frequency within each age category (2.2% at age 6, 3.6% at age 7, 6.2% at age 8, 3.7% at age 9, 3.8% at age 10, 4.5% at age 11, 4.1% at age 12, 5.7% at age 13, 4.3% at age 14, 2.6% at age 15, and 4.4% at age 16) and sex (from 3.5% in boys to 4.6% in girls) was very low.

Based on information obtained from the Ministry of Education [27] on the distribution of students enrolled in the 2009-10 academic year (date of data collection), it was decided to intentionally select 31 public education institutions located in urban areas, of which 13 belonged to kindergarten, 10 at the primary school, and 8 at the secondary school. A total of 38% of the total number of students enrolled in the study, that is, 4 out of every 10 students who regularly attended school. The sample only included native children and adolescents from their respective regions (i.e., no migrants were included), regardless of whether their parents or grandparents were migrants or not. We gathered information about their birth place and current place of residence or address residence and cross-checked this information with the participants’ identity card. All assessments were carried out in Barranco between November and December (2009) and April and July (2010), in La Merced and San Ramon between May and August, and in Junín between September and October of 2009.

After initial political and educational contacts with local authorities in each city, formal permission was obtained from schools’ governmental bodies to participate in the study. Informed consent was provided by parents/legal guardians. The Ethics Committee of the National University of Education Enrique Guzmán y Valle (UNE EGyV) as well as all school directors approved the project. Table 2 shows the number of children and adolescents by age, sex, and geographical area of residence. The information about sample size and frequencies for BMI categories (normal weight and overweight) according to IOTF [29] and WHO [30] cutoff points, by age, sex, and geographical area of residence can be found in Supplementary Tables 1 and 2, respectively (please see these tables in the Supplementary Material for comprehensive data analysis).

2.3. Child-Level Correlates
2.3.1. Anthropometry

All measurements were made according to standardized techniques [31]. Height and sitting height were measured with the head positioned to the Frankfurt plane to the nearest 0.1 cm with a portable stadiometer (Sanny, Model ES-2060); body weight was measured at the nearest 0.1 kg using a digital scale (Pesacon, Model IP68). Body mass index (BMI) was computed using the standard formula (weight (kg)/height (m)2), and subjects were classified as either normal weight or overweight (including obesity) according to the International Obesity Task Force (IOTF) [29] and World Health Organization (WHO) [30] cutoff points. A dummy coding was created, with normal weight coded as 0 (reference category) and overweight (including obesity) coded as 1.

2.3.2. Biological Maturation

Biological maturation was assessed with the maturity offset [32] which estimates the distance, in decimal years, each subject is from the age at peak height velocity (PHV). Formulas are sex-specific and use chronological age, height, sitting height, and leg length. A positive (+) maturity offset indicates the number of estimated years a child is beyond PHV; a negative (−) maturity offset is the number of estimated years a child is before the PHV, whereas a zero value indicates that a child is experiencing his/her PHV.

2.3.3. Physical Fitness

Three measures of physical fitness were used. Muscular strength (standing long jump), agility (shuttle run), and cardiorespiratory fitness (12 minute run) were obtained from two well-known test batteries, AAHPERD [33] and EUROFIT [34]. All tests were performed according to their standardized protocols: (1) Standing long jump test: each subject stood with feet apart behind the takeoff line and, without a preparatory run, were instructed to jump as far as possible. The maximum jumping distance was recorded. Two trials were given, and the best score was recorded as the maximum jumping distance in centimetres. (2) Shuttle run test: each subject performed five cycles (round-trip) at maximum speed between two lines separated by five meters. (3) 12 minute run test: in a previously delimited field, schoolchildren ran/walked the maximum possible distance in 12 minutes. All physical fitness test results were converted to z-scores using a grand mean centering, and a total physical fitness score (PFz) was obtained by summing all individual z-scores as suggested [35].

2.4. Sociogeographic Context

Geographical areas of residence were located at different altitudes: sea level (58 m), Amazon region (751 m), and high altitude (4107 m). A coding scheme was generated with two dummies: sea level, the reference category, was coded as 0 0, high altitude was 1 0, and Amazon region was 0 1.

2.5. School-Level Correlates

The description and inventory of characteristics associated with the school contexts were obtained via an objective school audit, using a modified and locally adapted version of the Healthy Eating and Physical Activity modules of the Healthy School Planner designed by the Joint Consortium for School Health [36]. This audit was led by the experienced research team members and maps various domains: (1) school characterization (school size: number of students was divided by 10 so as to have a more interpretable odds ratio; school setting: mixed (the reference category) and urban); (2) policies and practices for physical activity (specifically, the existence of policies and practices issued by the state, school board, or any other government agency to promote physical activity, health, and well-being of students and are organized by the school) (“no policies and practices” was the reference category); (3) physical structure of the school (playground area: “with obstacles” was the reference category and multisports-roofed: “yes” was the reference category); (4) Physical Education classes (frequency: “one” was the reference category and duration: “90 min” was the reference category); and (5) extracurricular activities (“yes” was the reference category).

2.6. Data Quality Control

To certify data quality, the following procedures were used: (1) training of all team members by experienced researchers in the correct use of the technical procedures of body measurement and physical fitness test protocols; (2) conducting a pilot study to assess the quality of data collection; (3) retesting of a random sample of 29 children and adolescents; and (4) estimating the interobserver technical error of the measurement for anthropometric data (height TEM = 0.2 cm; sitting height TEM = 0.1 cm; body mass TEM = 0.1 kg) and ANOVA-based intraclass correlation coefficient (ICC) for physical fitness tests (standing long jump: (ICC = 0.85), shuttle run (ICC = 0.81), and 12 minute run (ICC = 0.79)).

2.7. Statistical Analysis

Exploratory analyses, outlier identification, and normality checks were completed prior to further analyses. Differences in frequencies of overweight by age, sex, and geographical area of residence were computed using a chi-squared (χ2) test. Differences in means and frequencies of child-level variables between geographical areas of residence were computed using analysis of variance (ANOVA), as well as χ2 test, respectively. Additionally, the Tukey HSD test was used for multiple comparisons. All analyses were performed in SPSS software version 24.0.

Because our data had a hierarchical structure, i.e., children nested within schools, we used a two-level random intercept (multilevel) logistic regression model for BMI categories (normal weight versus overweight). A three-level structure was considered—subjects nested within schools which are nested within regions. However, with only three regions, it is not recommended to treat region as a level in a multilevel analysis [37]. Instead, dummy variables for region were included as predictors of BMI categories and were included in the fixed part of the model.

A series of nested models were fit to explain variation in BMI categories using the deviance statistic as a measurement of relative fit [38]; differences in deviances follow a χ2 distribution with degrees of freedom equal to the difference in the number of estimated parameters from both models. It is expected that better fitting models have lower deviance values. Modelling was done in a sequential manner as generally advocated [37, 38]. Firstly, a null model (M0) was fit to the data to compute the intraclass correlation coefficient to estimate the variance accounted for by the random school effects in BMI categories. Secondly, Model 1 (M1) using child-level BMI predictors (age, sex, age-by-sex interaction, maturity offset, and PFz) was fit. Thirdly, in the Model 2 (M2), school-level covariates were added (number of students, school setting, policies and practices for physical activity, playground area, multisports-roofed, frequency and duration of Physical Education classes, and extracurricular activities). In the last model, Model 3 (M3), the predictor of geographical area of residence was added. Covariates were centered at their grand mean when necessary as is common in multilevel analysis [38]. Results are presented as odds ratios and their 95% confidence intervals (95% CI). The multilevel analyses was performed in SuperMix software [39], allowing a simultaneous estimation of all model parameters using maximum likelihood procedures. Statistical significance was set at 5%.

3. Results

Table 3 provides information about the prevalence of overweight (including obesity) among Peruvian children and adolescents, by age, sex, and geographical area of residence, based on IOTF and WHO cutoff points. Using a χ2 test, significant differences () in the prevalence of overweight were found by age (χ2 = 125.53, ), sex (χ2 = 4.68, ), and geographical area of residence (χ2 = 782.51, ). In the total sample, the prevalence of overweight was 20.8% (95% CI = 20.0–21.7). Further, the prevalence decreased with age (from 28.3%, 95% CI = 24.8–32.2 at 6 yrs to 13.9%, 95% CI = 10.4–17.4 at 16 yrs); was higher in girls (21.7%, 95% CI = 20.6–22.9) than in boys (19.8%, 95% CI = 18.5–21.0); and was higher among sea-level children and adolescents (41.3%, 95% CI = 39.1–43.7), followed by the Amazon region (18.8%, 95% CI = 17.6–19.9) and high altitude (6.3%, 95% CI = 5.3–7.3).

Table 4 provides descriptive statistics for the child-level variables. Statistically significant differences () were found for all variables among subjects living at different altitudes. High-altitude subjects are relatively older than their peers from the other areas (F = 62.45, ; 11.6 versus 11.2 and 10.8 yrs). On average, Amazon region subjects lag behind their peers in biological maturation (F = 75.19, ), and sea-level children are more physically fit compared with their peers from high altitude (F = 1371.35, ). The frequency of girls is higher than boys in the three geographical regions (sea level: girls 61.1%, boys 38.9%; Amazon region: girls 52.7%, boys 47.3%; high altitude: girls 51.4%, boys 48.6%) (χ2 = 48.69, ).

Descriptive statistics for the school-level variables are presented in Table 5. The number of students in the schools ranged from 96 to 1200; 83.3% of the schools were located in an urban region; 16.7% had policies, and 38.9% had practices for physical activity. Furthermore, 83.3% of schools had a recreation space without obstacles and 66.7% did not have multisports-roofed complexes. In 94.4% of schools, Physical Education classes were conducted once per week and 55.6% of them lasted >90 min. Finally, 77.8% of schools provided extracurricular activities for schoolchildren.

The results of the multilevel models, based on IOTF cutoff points, are presented in Table 6. The results are very similar (in magnitude and direction) for both criteria (for WHO results, please see Supplementary Table 3 in the Supplementary Material for comprehensive data analysis). The M0 indicated that school-level effects expressed by the intraclass correlation coefficient was 0.209, meaning that ∼21% of the total variance in BMI categories among subjects was at the school level and that ∼79% of the variance was explained by child-level characteristics, respectively.

Results from M1 (child-level predictors) showed that younger subjects (OR = 0.16, 95% CI = 0.14–0.18) as well as those ahead in their maturity status were more likely to be overweight (OR = 14.19, 95% CI = 12.00–16.78). The age-by-sex interaction was statistically significant (OR = 0.74; 95% CI = 0.70–0.78) revealing that boys tend to be less likely to be overweight than girls as they get older. Furthermore, being more physically fit (OR = 0.80; 95% CI = 0.76–0.84) protects subjects from being overweight. There was a significant reduction in deviance from M0 to M1 (from 7927.97 to 6524.56; χ2(5) = 1403.41, ), showing the better fit of M1.

In M2, school-level predictors were added. Six out of nine predictors were statistically significant (). Children from schools with a higher number of students (OR = 1.02; 95% CI = 1.01–1.03) and without multisports-roofed (OR = 2.44; 95% CI = 1.53–3.89) were more likely to be overweight. Surprisingly, children and adolescents from schools with physical activity policies (OR = 5.09; 95% CI = 2.61–9.94) and practices (OR = 2.83; 95% CI = 1.83–4.38), with Physical Education classes with more than 90 minutes of duration (OR = 3.02; 95% CI = 1.84–4.96), as well as with extracurricular activities (OR = 0.59; 95% CI = 0.37–0.92), had a greater likelihood of being overweight. There was a reduction in deviance from M1 to M2 (from 6524.56 to 6480.93; χ2(9) = 43.63, ), showing the better fit of M2.

In the last model (M3), geographical area of residence was added and results showed that subjects living at high altitude (OR = 0.08; 95% CI = 0.04–0.18) were less likely to be overweight than those living at the sea level. In this model, subject-level covariates remained similar as in M1 and M2, but the importance of some school-level covariates changed. Children and adolescents from schools with a lower number of students (OR = 0.97; 95% CI = 0.96–0.98), with playground areas with obstacles (OR = 0.09; 95% CI = 0.04–0.24), with multisports-roofed (OR = 0.40; 95% CI = 0.24–0.66), and with only one Physical Education class per week (OR = 0.03; 95% CI = 0.01–0.13), were more likely to be overweight. There was a reduction in deviance from M2 to M3 (from 6480.93 to 6453.66; χ2(2) = 27.27, ), showing the better fit of M3.

4. Discussion

Since the last century, Peru has been experiencing epidemiological and nutritional transitions with unforeseen changes in the prevalences of overweight and obesity. Here, we reported a combined prevalence of 20.8% for overweight and obesity in children aged 6 to 16 yrs, as well as significantly lower prevalences at higher ages. Moreover, we found that children and adolescents living at the sea level had the highest overweight prevalence, followed by those from the Amazon region and high altitudes. Álvarez-Dongo et al.[21], using WHO cutoff points, reported a combined prevalence of 24.4% in children aged 5 to 9 yrs, decreasing to 14.2% in youth aged 10 to 19 yrs, which are smaller than what we reported using the same cutoff points within the age range of 6 to 9 yrs (31.6%) and 10 to 16 yrs (20.7%). In addition, they also reported a similar trend, with the highest prevalence of overweight and obesity found among children and adolescents living in the coastal regions. These geographical differences may be due to the differences in urbanization levels and lower poverty observed in coastal districts, which may lead to changes in lifestyles, namely, food habits and physical activity levels [40].

We also observed sex differences in the prevalence of childhood/adolescence overweight, favouring boys. Available systematic review [41] revealed similar trends within and between countries, emphasizing that in developing countries, these sex disparities are more exacerbated and are commonly observed before and during puberty. Several factors may contribute to this sex difference: firstly, boys and girls differ in body composition as well as in weight gain patterns [12]; secondly, boys tend to carry out heavy work and physical activities with greater energy expenditure, whereas girls are more responsible for household tasks, involving less energy expenditure [42].

A major finding of our study was that the main fraction of the total variance in BMI categories (∼79%) was explained by child-level characteristics. These results are congruent with those obtained in the International Study of Childhood Obesity, Lifestyle and the Environment (ISCOLE), where 90% of the variance in obesity was explained at the childhood level among 9- to 11-year-old children from 12 countries [43]. Further, ISCOLE demonstrated that 92% to 94% of the variance in BMI and waist circumference, respectively, was explained at the childhood level, with the remainder being explained at the school and site levels [44].

We also showed that younger children were more likely to be overweight, corroborating Cremm et al.’s [45] study with Brazilian children linking their results to the fact that younger children were more prone to consume healthier foods less frequently and engage in less moderate-to-vigorous physical activities or sports as compared to older children. Boys tended to be less likely to be overweight than girls as they get older (the age-by-sex interaction was statistically significant). There is a tendency for girls to surpass boys in their overweight trends, which are also linked to the fact that girls exhibit, on average, higher amounts of body fat and boys higher amounts of muscle mass [12]. These sex differences in body composition are evident from fetal life but emerge primarily during prepubertal and pubertal phases of life. Girls enter puberty earlier—they experience their midgrowth spurt and their adolescent spurt at earlier ages and undergo a more rapid pubertal transition, whereas boys have a substantially longer growth period [46]. Furthermore, the action of sex steroid hormones (estrogen) begins earlier in girls, preparing them for menarche, which results in a greater accumulation of body fat [47], predisposing them to a greater risk of being overweight.

We also showed that children and adolescents ahead in their maturity status were more likely to be overweight, which agrees with [48] data in Portuguese children. This is expected since previous research has shown that the timing of biological maturation further seems to contribute to greater adiposity development and early maturing youths usually are taller and heavier and have higher BMI than their later maturing peers [49]. As such, our results were expected since those closer to their PHV had higher BMI values and, consequently, had a greater likelihood of being overweight. However, to the best of our knowledge, we were not able to identify information about biological maturation and its associations with body weight in Peruvian children and adolescents to compare with our results. We were only able to find reports on girls [50] based on their age at menarche, showing that those living at the sea level had a higher BMI than those living at high altitude; their age at menarche occurred earlier which may also explain the higher prevalence of overweight previously reported among sea-level girls compared with their peers from the other regions.

Physical fitness is an important health marker [51], and there is emerging evidence showing that adequate physical fitness levels are associated with the prevention of overweight during childhood and adolescence, playing an important protective role in its rate of development [52]. Our results confirmed this link, showing that muscular strength, agility, and cardiorespiratory fitness have negative associations with overweight status. This trend was also found by McGavock et al. [53] in that low cardiorespiratory fitness and reductions in physical fitness over time were positively related to weight gain and increased risk of overweight in children 6–15 yrs old. Similarly, Stodden et al. [54] showed that the development of adequate levels of motor competence and health-related fitness may be a key component to prevent “unhealthy” weight gain in children and adolescents. Therefore, we emphasize the important role of physical fitness as a protective factor to reduce overweight/obesity when developing efficient intervention programs in children and adolescents.

Contextual information appears to be important in explaining why children vary in their weight status. Living at the sea level, particularly in the region of Barranco, was related to a greater likelihood of being overweight compared with high altitude, which may indicate that overweight is not yet a nationally distributed problem in this country [19]. One possible explanation for this discrepancy may be due to distinct lifestyle patterns in the three geographical areas. Barranco has a high density population of 10.132 inhabitants per km2, increasing public insecurity and the lack of public recreational infrastructures available to schoolchildren, contributing to the adoption of sedentary lifestyles. In addition, the existence of several fast food restaurants, convenience stores, and marketing of unhealthy food choices increases the likelihood of children acquiring unhealthy eating habits. Conversely, Chanchamayo and Junín have a lower population density (10.08 inhabitants per km2 and 12.22 inhabitnat per km2, respectively), and the majority of the schoolchildren participate in recreational activities, housework, and supporting agricultural activities in the field, showing higher levels of physical activity compared with their peers from the sea level. Additionally, the market penetration is less common and the main livelihoods are some tourism, mainly in the Amazon region, but predominantly stockbreeding and nondiversified agriculture, which limit the nutritional intake and, consequently, affect children and adolescent’s health and weight status [42].

On the other hand, it has been suggested that lower obesity rates among highlanders as compared to lowlanders can also be explained by direct effects of living at high-altitude on human physiology (e.g., loss of appetite due to increased leptin concentrations [55] and increased thermogenesis due to cold temperatures that increase thyroid hormone [56] and catecholamine levels [57]). Thus, it is possible that living at high altitude has a weight-lowering effect per se, independent of various confounders factors, including multiple expressions of families and socioeconomic status as proposed by Woolcott et al. [58].

Additionally, lower obesity rates among high-altitude children and adolescents may also be considered in terms of the “reverse causality” theory. It suggests that overweight/obese individuals might tend to migrate to lower altitudes where the natural environment is friendlier in terms of performing their daily tasks/chores and activities and/or because they have easier access to healthcare services [58, 59]. However, our sample only included native children and adolescents from their respective regions (i.e., no migrants were identified).

It is widely accepted that schools are an important factor in understanding variability in childhood overweight [22]. Our results confirmed this, given that ∼21% of the total variance in BMI categories was attributed to the schools children attended. In previous studies [23, 24], lower values of variance (ranging from 2% to 5%) were found to explain school-level differences in the variability in the odds of a student being overweight. This discrepancy can be due to the higher dissimilarity observed across school environments in developing countries, such as in Peru, when compared to developed countries.

We found that children attending schools with a lower number of students were more likely to be overweight, whereas those attending schools with playground areas without obstacles and without multisports-roofed were less likely to be overweight. Because school size was defined according to the number of enrolled students, schools with more students apparently have more space and facilities for children to play during their recess time and have greater potential to increase children’s physical activity levels [60], which could also promote positive effects on their weight status. Furthermore, it is possible that a playground without obstacles represents a greater free space for children to play, acting as a protective agent from overweight during childhood.

In an effort to fight childhood overweight and obesity, school-based healthy eating and physical activity programs have been implemented across the globe, providing opportunities to enhance children’s future health and well-being [61]. Surprisingly, we found an inverse effect in our results probably because the effectiveness of these programs is not yet very well established. For example, in a systematic review of intervention studies, Campbell et al. [62] tracked only seven studies on prevention of childhood obesity: four revealed effective programs, whereas three did not. In fact, the lack of evidence exploring how policies and practices within the school system are associated with an individual student’s risk of being overweight continues to represent an important gap in the available literature [61]. Finally, Peruvian students from schools with two Physical Education classes per week, with a maximum of 90 minutes’ duration, were less likely to be overweight. In Peru, physical education classes are mandatory by government authorities, and all students are engaged in at least 1 class per week. As such, it is of no surprise that students who have 2 classes per week were less likely to be overweight.

This study is not without limitations. Firstly, the cross-sectional nature of the design does not allow a dynamic analysis of intraindividual changes and interindividual differences in overweight/obesity unfolding during childhood and adolescence. Yet, no study is apparently available in South America sampling hundreds/thousands of children and adolescents from different sociogeographic areas. Secondly, notwithstanding the size of our sample, it is not representative of the Peruvian population and care must be taken when trying to generalize our results. Thirdly, we were not able to obtain information about dietary habits and objective measures of physical activity. Fourthly, our results were not adjusted for the putative influence of family income and/or insurance status. Despite these shortcomings, the study has several merits: firstly, we sampled children aged 6 to 16 yrs which represent an important growth and developmental time window; secondly, we used standard methods and highly reliable data at the individual and school levels; thirdly, we relied on multilevel modelling to capture the complexity of the nested information—children within their schools; fourthly, the inclusion of three different geographical areas of residence permitted a more comprehensive interpretation of differences in childhood.

5. Conclusions

In summary, the present study showed that the prevalence of overweight in Peruvian children and adolescents decreased with age, was higher in girls, and was also higher at the sea level compared to the Amazon region and high altitude. Child-level variables (age, sex, maturation, PFz, and geographical area of residence), as well as school-level variables (number of students, policies and practices for physical activity, playground area, multisports-roofed, and frequency and duration of Physical Education classes), played important roles in explaining variations in childhood overweight. This information should be taken into account when designing more efficient strategies to combat the overweight and obesity epidemic.

Data Availability

The datasets generated and analyzed during the current study are not publicly available due to privacy laws associated with children’s data but are available with a data-sharing agreement as approved by the Ethics Committee of the National University of Education Enrique Guzmán y Valle, Lima, Peru.

Disclosure

The authors alone are responsible for the content and writing of the paper.

Conflicts of Interest

The authors report no conflicts of interest.

Authors’ Contributions

CS, DH, and JM analyzed the data. CS, JM, DH, PTK, and AB drafted the manuscript. AB and JM designed the study and directed its implementation. AB and JM provided necessary logistical support. AB, DH, OV, RG, and PTK provided extensive critical comments and did the final editing of the manuscript.

Acknowledgments

The authors thank the schoolchildren of Barranco, Junín, San Ramon, and La Merced who participated in this study, as well as their parents. We also thank UNE’s Physical Education students and teachers for their high professionalism in all steps of the data collection. Finally, we would like to recognize the late Emeritus Professor Gaston Beunen who was one of the prime movers of the Healthy and Optimistic Growth Study. This work was supported by the Foundation for Science and Technology, Portugal (Grant/Award Number SFRH/BD/141112/2018), and POCH of European Social Found (FSE).

Supplementary Materials

Supplementary Table 1: sample size and frequencies (n (%)) for BMI categories (normal weight and overweight) according to IOTF cutoff points, by age, sex, and geographical area of residence. Supplementary Table 2: sample size and frequencies (n (%)) for BMI categories (normal weight and overweight) according to WHO cutoff points, by age, sex, and geographical area of residence. Supplementary Table 3: multilevel modelling results: odds ratios (OR) and 95% confidence intervals (95% CI) for child- and school-level characteristics. All these supplementary materials were referenced at appropriate sections in the manuscript. (Supplementary Materials)