Abstract
Introduction Research showing Place of Birth (POB) predicts excess weight gain and obesity risk among Latino adults has not prompted similar research in Latino children, although childhood is a critical period for preventing obesity.
Methods Longitudinal cohort observational study on public school children self-identified by parent/guardian as Latino in grades K-12 for school years 2006-07 through 2016-17 with measured weight and height (n= 570,172students; 3,10,3642observations). POB reported by parent/guardian was categorized as continental US (not NYC) (n=295,693), NYC (n=166,361), South America (n=19,452), Central America (n=10,241), Dominican Republic (n=57,0880), Puerto Rico (n=9,687) and Mexico (n=9,647). Age- and sex-specific BMI percentiles were estimated based on established growth charts. Data were analyzed in 2020.
Results Prevalence of obesity was highest among US (non-NYC)-born girls (21%) and boys (27%), followed by NYC-born girls (19%) and boys (25%). Among girls, South Americans (9%) had the lowest prevalence of all levels of obesity, while Puerto Ricans (19%) and Dominicans (15%) had the highest prevalence. Among boys, South Americans also had the lowest prevalence of all levels of obesity (15%), while Puerto Ricans (22%) and Mexicans (21%) had the highest. In adjusted models, obesity risk was highest in US (non-NYC)-born children, followed by children born in NYC (p<0.001). Immigrant Latino children exhibited an advantage, particularly Dominicans, South Americans and Puerto Ricans, and even after controlling for individual and neighborhood sociodemographic features including linguistic isolation and poverty.
Conclusions The heterogeneity of obesity risk among Latino children highlights the importance of POB.
Introduction
Childhood obesity in the U.S. has been increasing since the 1980s, and Latino children continue to be disproportionately affected.1 The latest national estimates from 2013-2015 show that Latino youth aged 2 to 19 years have a 23.6% prevalence of obesity and 8.4% prevalence of extreme obesity, compared to 14.7% and 4.4%, respectively, among their non-Latino white counterparts.2-4 Widening disparities nationwide have been documented among Latino children relative to white children.5 At the same time Latino children account for 25.6% of the US child population, and childhood obesity is predictive of a myriad of chronic health conditions in adulthood including asthma, arthritis, and poorer cardiometabolic and psychological risk profiles.6-8
Research in the U.S. has typically characterized obesity prevalence among Latino children by comparing them in aggregate to other racial/ethnic groups.1, 2, 5 This is despite research documenting a Latino “immigrant advantage,” or better-than-expected obesity profiles of immigrants versus those born in the U.S., despite lower SES profiles.9-11 Obesity research among adults has also highlighted the importance of Place of Birth (POB),12 including genetic, environmental, and socio-cultural factors13, 14 that may further drive subgroup differences among Latino children. Mexican American children have demonstrated lower obesity prevalence among children born in Mexico versus the U.S., particularly among girls.15-17 Other research on Latino children corroborates sex differences in obesity risk irrespective of POB.18, 19 To date, limited work has studied obesity prevalence among large and diverse samples of Latino children by POB and sex.
In New York City (NYC), one of the most urban and diverse areas in the U.S., statistically significant decreases are reported in the prevalence of obesity among Latino children, from 21.9% in 2006-07 to 20.2% in 2016-17.20 Yet, significant disparities remain among Latinos relative to their non-Latino white counterparts despite a comprehensive approach to curb obesity by the NYC Departments of Health and Mental Hygiene (DOHMH) and Education (DOE) among vulnerable communities of color.21-23 In this paper, we examined heterogeneity in obesity risk among Latino children attending NYC public schools by POB, while accounting for important individual and neighborhood level SES and sociocultural factors associated with obesity risk among Latino adults.24, 25 Identifying differences in obesity risk among Latino children by POB is important for forecasting trends in prevalence and related outcomes among this growing population.
Methods
Participants
Data for this study were drawn from the NYC Fitnessgram dataset jointly managed by NYC Department of Education (DOE) and DOHMH, and have been described elsewhere.20, 26-29 Briefly, teachers collected child-level student height and weight annually in NYC public schools. Students’ height and weight measurements were taken annually during physical education classes among K-12 students as part of the NYC Fitnessgram curriculum using a standard protocol. The study population included all NYC Latino students enrolled in a general education public school during the 2006-07 through 2016-17 school years with at least one year of height and weight measurements. Age was defined as student age on December 31st of the given school year based on NYC Fitnessgram measurements that are collected annually within 3 months of this date based on testing schedules.30
All child-level student demographic data were drawn from NYC DOE student enrollment records linked to Fitnessgram data by a unique identifier. The DOE uses demographic information, including Latino ethnicity, only for programmatic, instructional, and administrative planning and decision-making.31 Parents and guardians were asked to identify whether the child was “Hispanic” (defined as “a person of Mexican, Puerto Rican, Cuban, Central or South American, or other Spanish origin”) regardless of race. Hereafter, we use the term Latino since our sample included children from, or descended from, Spanish-speaking people from Latin America. This study was approved by the City University of New York (IRB File #2015–0582) and DOHMH (Protocol # 14–019) Institutional Review Boards.
Measures
Independent Variable
Parents/guardians reported all demographics, including child POB, which was used to derive seven categories: continental U.S. (not NYC or U.S. territories, including Puerto Rico) (n=295,693 unique students), NYC (n=166,361 unique students), South America (n=19,452 unique students), Central America (n=10,241 unique students), Dominican Republic (n=57,0880 unique students), and Mexico (n=9,647 unique students). These subgroups were selected based on prior literature documenting important POB differences in obesity risk in national samples of adults and NYC.12, 33 These documented differences also drove the decision to separate Puerto Rican children from their continental U.S. counterparts (n=9,687), despite the fact that Puerto Rican children are all U.S.-born by definition. Cuban children (n=282) were excluded due to insufficient sample size to generate reliable estimates.
Dependent Variable
The primary outcome of interest was age- and sex-specific BMI percentiles in accordance with the Centers for Disease Control and Prevention (CDC) 34 growth charts.35 Baseline age in months was calculated from the measurement date and students’ date of birth was drawn from school enrollment records. Extreme or biologically implausible values (BIV) were identified for height, weight, weight-for-height, and BMI using CDC’s age- and sex-specific criteria.34 An observation identified as BIV for a student in a single school year was excluded only for that school year (n=1,756 observations). The final analytic sample included 570,172 students or 3,10,3642 observations from 2006-07 through 2016-17 school years, with individual children having 1 to 11 repeated annual observations and 56% of the sample having at least 5 repeated annual observations. Weight status was defined according to CDC’s growth chart-derived norms for sex and age in months and used to compute the BMI percentile for each child as follows; underweight (BMI<5th percentile), normal (5th percentile≤ BMI<85th percentile), overweight (85th percentile≤ BMI< 95th), obese (BMI≥ 95th percentile). Obesity was stratified in accordance with previous reports34, 36 as follows; class I (BMI ≥95th percentile), class II (BMI >120% of the 95th percentile for age and sex or a BMI of ≥35), class III (BMI ≥140% of the 95th percentile for age and sex or a BMI of ≥40 or greater).
Covariates
To categorize students in terms of SES, individual student household poverty (high vs. low) was based on student eligibility/non-eligibility for free/reduced-price school meals through the National School Lunch Program, which provides meal assistance according to household income at or below 185% of the federal poverty level.37 Area-based SES was defined according to the American Community Survey (ACS) 2012-2016 data as the percentage of households in the student’s home zip code receiving food aid from the Supplemental Nutrition Assistance Program.38 Neighborhood linguistic isolation was included given its significant association with BMI and obesity-related outcomes among Latinos in NYC;25, 39 it was measured based on the percentage of Spanish-speaking only households in the home zip code drawing from the ACS 2012–2016 and then categorized into quartiles.
Statistical Analyses
Data were analyzed in 2020. Descriptive statistics were computed to summarize sample characteristics. Means and SDs were generated for all continuous data (BMI percentile, area poverty, and linguistic isolation), while categorical data (sex, POB, weight category, household poverty status, categorized area poverty and linguistic isolation) were reported as frequencies and percentages for all students. Continuous variables were also generated across POB, and weight categories were computed by sex and POB. Chi-square tests of homogeneity were performed to compare the distribution of weight categories across sex and POB subgroups.
Next, repeated measures mixed models were used, where repeated observations were nested within individual children who in turn were nested within home zip codes. These models estimated the association between individual-level POB and BMI percentile over time, taking into account individual- and area-level factors for all students, and also stratified by sex. Additional models were run to test whether slopes in BMI percentile over time (calendar year) were significantly different across the POB subgroups. For all models, random intercepts for area effects (home zip code) and random slopes for time point of observation (school year, to examine period cohort effects) were included with participants and zip codes as the subjects for levels 2 and 3, respectively, to account for between-child and between-zip code variations. A first-order autoregressive covariance structure was used to represent the correlated repeated measurements over time within participants in all models.
Adjusted models included sex (unstratified models), baseline age at time of BMI measurement (continuous variable), household poverty status (binary variable), home area poverty (categorical variable), linguistic isolation (categorical variable), and time (an integer value increasing from 0 to 10 corresponding to the number of repeated observations or years that each child was observed in the dataset) as covariates. Models testing whether slopes in BMI percentile over time (calendar year) were significantly different across the POB subgroups also included calendar year as a categorical variable in the models. A p-value of <0.05 was considered statistically significant, and all p-values were 2-sided. Statistical analyses were performed using SAS v.9.4.40
Results
Descriptive characteristics for Latino children in grades K-12 (nchildren=809,418; nobservations=3,103,642) appear in Table 1. The majority, 52%, of the analytic sample was born in the continental US or NYC (29%). Of the children born outside the continental US, most were from Dominican Republic (10%), followed by South America (3%), Puerto Rico (2%), Mexico (2%), and Central America (2%). There was an equal proportion of females and males, with a high household (81%) and neighborhood (44%) poverty rate, and half of the sample resided in areas with high (quartile 3) or very high (quartile 4) linguistic isolation. More than half of the children (57%) were in the healthy weight category, 21% experienced overweight, 19% were in one of the three obesity categories, and 2% experienced underweight.
Supplemental Figure 1-2 and Table 2 illustrate the six weight categories for all children by POB and sex, taking into account Fitnessgram measurements collected annually from each child. Among Latina girls born in the continental US but not NYC, 21% experienced overweight, 14% were affected by Class I obesity, 5% were affected by Class II, and 2% were affected by Class III. NYC-born girls showed similar prevalence rates for all obesity categories (13%, 4%, and 2%, respectively). Girls born outside the US had a similar or lower proportion of overweight (not including youth with obesity) compared with those in the US, except for Mexican (24%) and Central American girls (22%). South and Central American, Dominican, Puerto Rican, and Mexican girls had a lower proportion of obesity in all categories compared to their US and NYC-born counterparts except for Classes I and II obesity in Puerto Rican girls (13% and 4%, respectively), which was consistent with NYC-born girls.
Latino boys (Supplemental figure 2 Table 2) had a higher proportion of overweight and obesity compared to girls (chi-square test comparing weight categories across POB p<.001 for all). Specifically, among Latino boys born in the continental US but not NYC, 19% experienced overweight, and 18%, 7%, and 2% experienced Classes I - III obesity, respectively. Boys born in NYC showed comparable prevalence for all obesity categories (17%, 6%, and 2%, respectively). Among children born outside the US, boys generally had a higher proportion of overweight but a lower proportion of Class I obesity (ranging from 12-16%) and Class II (ranging from 3-4%) than girls. The sole exception was Puerto Ricans, where boys had lower overweight (17.8%) than girls (19.1%) but higher Class II obesity (5.6%) than girls (4.2%). Puerto Rican boys also had a comparable prevalence of Class III compared to continental US- or NYC-born children (∼2%). All other subgroups born outside the continental US or NYC had ∼1% of boys in the Class III category.
Crude and adjusted repeated measures mixed models for the association between POB and BMI percentile for all children and stratified by sex are presented in Tables 3 and 4, respectively (Type III Fixed Effects Estimates p<.001 for all). Models adjusted for sex (unstratified model), baseline age at time of BMI measurement, household poverty status, home area poverty, linguistic isolation, calendar year, and time as covariates showed the strongest association between POB and BMI percentile in Dominican (b=-3.67, 95%CI: −3.79, −3.56), South American (b=-3.21, 95%CI: −3.40, −3.02) and Puerto Rican (b=-3.15, 95%CI: −3.44, −2.87) children, relative to the reference group (US, non-NYC-born; Table 3). After stratifying by sex, Dominican (b=-3.42, 95%CI: −3.58, −3.26), and South American (b=-3.24, 95%CI: −3.51, −2.98) girls demonstrated significantly lower BMI percentile than US non-NYC-born girls. Among boys, Puerto Ricans (b=-4.18, 95%CI: −4.59, −3.78), Dominicans (b=-3.94, 95%CI: −4.11, −3.78), Central Americans (b=-3.39, 95%CI: −3.73, −3.04) and South Americans (b=-3.19, 95%CI: −3.46, −2.92) all showed lower BMI percentile than US non-NYC born children. Also, the largest disparities by sex within the same POB subgroup were observed in Central American (b=-1.94, 95%CI: −2.30, −1.58 vs. b=-3.39, 95%CI: −3.73, −3.04) and Puerto Rican (b=-2.16, 95%CI: −2.56, −1.78 vs. b=-4.18, 95%CI: −4.59, −3.78) girls and boys, respectively.
Additional models examining whether slopes over time (calendar year) were significantly different across the POB subgroups showed that Puerto Rican children had the greatest decline in BMI percentile over time (b=-2.38, 95%CI: −3.31, −1.44), followed by Dominican (b=-1.99, 95%CI: −2.33, −1.65) and Central American b=-1.70, 95%CI: 2.42, −0.98) children; (Type III Fixed Effects Estimates p<.001 for all; Table S1 and Supplemental Figure 3). Models stratified by sex showed similar findings among boys (b=-2.79, 95%CI: −4.13, −1.44, b=-2.40, 95%CI: - 2.89, −1.91 and b=-2.31, 95%CI: −3.30, −1.32 for Puerto Rican, Dominican and Central American children, respectively). However, among girls, Puerto Rican (b=2.53, 95%CI: −3.83, 1.24), Dominican (b=-2.17, 95%CI: −2.65, −1.70) and South American (b=-2.07, 95%CI: −2.95, −1.19) children showed the greatest declines in BMI percentile.
Discussion
These findings show that Latino children attending NYC public schools in grades K-12 have a high burden of excess weight and obesity risk; however, this burden is not distributed evenly across POB. Specifically, children born in the continental US had the highest proportion of excess weight and obesity risk, followed by children born in NYC. Children born outside the US did exhibit an immigrant advantage, but this was most pronounced among Dominicans, South Americans, and Puerto Ricans and, to a less degree, among Central Americans and Mexicans. This early immigrant advantage among Latino children could potentially translate into protective effects documented among immigrant adults for metabolic diseases.10, 41, 42 However, prevalence of obesity by POB among Latino adults revealed the highest prevalence among Puerto Ricans and the lowest among South Americans.4 This suggests that the early immigrant advantage observed in children may not confer the same benefits into adulthood for all Latinos, and a robust line of sociological research suggests the importance of understanding the social and economic factors that create segmented patterns of experiences by POB.44-46 The potential for obesity risk to follow similar segmented patterns by POB should be explored in future studies using longitudinal designs.
The literature has failed to address the heterogeneity in obesity risk among Latino children of different national or cultural origin. For example, previous surveillance data from the NYC DOHMH showed that the obesity prevalence among all Latino NYC public high school students was 17% among male and 14% among female students,47 masking the risk and advantage among some subgroups. This might explain why a recent study using data with grades K-8 showed that the decrease in obesity in terms of both absolute and relative prevalence remained smaller among Latino children than white children.20 Combining Latino children into one large category obscures relevant differences for understanding risk in excess weight and obesity in this growing segment of the population. Although research and practitioners have a growing recognition of differing health patterns and profiles among the various Latino groups, to date, there are almost no data on obesity of subgroups such as Dominicans, who constitute the fifth largest Latino group in the US.48 An exception is the Hispanic Community Children’s Health Study/Study of Latino Youth, though prevalence of obesity was stratified by those of Mexican background vs. Non-Mexican background.18
Further understanding of obesity risk among Latino children will likely require an exploration of other broad social determinants of health that might explain some of the patterns observed in this dataset. This study accounted for neighborhood-level sociodemographic and cultural factors as covariates, including area poverty derived from the proportion of residents receiving food assistance, as well as the proportion of residents speaking only Spanish in the home as a proxy for linguistic isolation. Others have pointed to a wider array of neighborhood and social factors that may be critical in obesity prevention throughout the life course,14 and factors that may be particularly important for Latino communities such as strategies that focus more on behavior change.49 Contextualizing these factors within a theoretical framework will be important to move beyond documenting differences; rather, we need to move into the next phase of research to explore the risk and resilience mechanisms that are shaped by complex sociocultural and SES factors.50, 51
This study has some limitations. First, the results do not include private, charter, and special education schools, which constitute around 18%, 10%, and 2% of elementary and middle school children, respectively, in NYC; since these children do not participate in NYC Fitnessgram. Second, not all NYC public school children participate in Fitnessgram every year, although since 2010 schools have been incentivized to collect data from a minimum of 85% of students who do not have a testing waiver. The majority of observations in this study came from younger children in elementary and middle school (60%), indicating that high school children disproportionately missed the NYC Fitnessgram assessment. However, the size and heterogeneity of the complete sample, and the ability to stratify analyses based on key demographic factors while maintaining statistical power, are strengths. Third, US-born and NYC-born children may live in households that are still deeply connected to countries throughout Latin America and this dataset does not allow for this exploration. Fourth, this dataset does not allow for a disaggregation of Central and South Americans.
Conclusions
The NYC public school system is the largest in the country, serving 1.1 million children and representing all countries in Latin America.52 Findings presented in this study highlight important disparities by POB and sex after controlling for individual- and neighborhood-level factors. This work combats misperceptions of weight in this population and calls for interventions that promote positive attitudes and behaviors toward achieving a healthy weight.53 Future research should also examine factors in the U.S. that influence immigration waves, neighborhood composition, intergenerational patterns, and that may be working in tandem to exacerbate obesity risk and/or resiliency among Latino children in different settings.
Data Availability
This surveillance data is not available without permission from the NYC Department of Health