Introduction
Obesity is a serious global public health challenge representing one of the leading causes of morbidity, mortality, and health care expenditure worldwide.Reference Lehnert, Sonntag, Konnopka, Riedel-Heller and Konig1,2 The global age-standardized prevalence of obesity in children increased more than eight-fold between 1975 and 2016; in 2016, there were 124 million children aged 5–19 worldwide who were obese.3 The consequences of childhood obesity are extensive and include multiple medical, social, and psychological comorbidities.Reference Rankin, Matthews and Cobley4,Reference Reinehr5 Moreover, there is considerable evidence to demonstrate that obesity in children persists into adulthood, with obese children being more than five times more likely to be obese as adults than nonobese children [pooled relative risk 5.21, 95% CI 4.50–6.02].Reference Simmonds, Burch and Llewellyn6
Although the rise in obesity rates has historically been thought to be largely driven by increased caloric intake and insufficient physical activity, there is mounting evidence that exposure to synthetic chemicals in our environment, called obesogens, plays an important role in the global obesity epidemic.Reference Heindel and Blumberg7 Moreover, studies suggest that exposures that occur during early life may have profound life-long impacts on body weight homeostasis.Reference La Merrill and Birnbaum8 To date, epidemiological studies have shown an association between exposure to environmental chemicals including pesticides, phthalates, flame retardants, bisphenol A, and polycyclic aromatic hydrocarbons (PAHs) and an increased risk of obesity.Reference La Merrill and Birnbaum8–Reference De Long and Holloway12
PAHs have long been pollutants of concern as known human carcinogens, and more recently because of mounting evidence of endocrine disrupting action. PAHs are mainly generated by incomplete combustion of organic materials and are ubiquitous in the environment; routes of exposure include inhalation, ingestion, and skin absorption. Metabolism of PAHs produces hydroxylated PAHs, which have estrogen-like structure and properties.Reference Schultz and Sinks13 Experimental work has demonstrated that exposure to benzo[a]pyrene, the most well-studied PAH, increases fat mass and leads to weight gain.Reference Irigaray, Ogier and Jacquenet14 PAHs may also alter serotonin signaling and impact feeding behaviours,Reference Stephanou, Konstandi, Pappas and Marselos15 as well as limit lipolysis directly through incorporation into adipocytes.Reference Irigaray, Lacomme, Mejean and Belpomme16 PAHs have been shown to activate early growth response genes and bind to peroxisome proliferator-activated receptors in cell culture systems,Reference Kim, Yamaguchi and Lee17 suggesting potential pathways beyond aryl hydrocarbon receptor binding by which PAHs could mediate impacts on a variety of biological systems, including lipid homeostasis and adipogenesis.Reference Grun and Blumberg18
There is some, albeit limited, evidence to suggest a link between PAH exposure and obesity in children aged 6–19 years.Reference Poursafa, Dadvand and Amin19–Reference Kim, Kam and Lee21 Interestingly, this association was stronger in the younger (i.e., 6–11 years of age) children versus adolescents (i.e., 12–19 years of age).Reference Scinicariello and Buser20 However, despite longitudinal data suggesting that body weight from 2 to 6 years old is the most critical for prediction and realization of adult overweight,Reference de Kroon, Renders, van Wouwe, van and Hirasing22 the evaluation of the association between PAH levels and obesity in children younger than 6 years old has not been determined. In the present study, we used data from the Canadian Health Measures Survey (CHMS, 2009–2015) to examine the association between urinary PAH metabolites and obesity in children aged 3–18 years.
Methods
Data source
The data are from the second (2009–2011), third (2012–2013), and fourth (2014–2015) cycles of the CHMS. The CHMS is an ongoing survey designed to provide comprehensive direct health measures at the national level and collects information from community-dwelling individuals. Full-time members of the Canadian Armed Forces and people living on reserves or in other Aboriginal settlements, in institutions, and in some remote regions are excluded (collectively representing approximately 3% of the Canadian population). The CHMS involves an in-person household interview and a subsequent visit to a mobile examination center (MEC). The household interview gathers general demographic and socioeconomic data and detailed health, nutrition, and lifestyle information. At the MEC, blood and urine are collected and direct physical measurements are taken, including height, weight, and waist circumference (WC). Detailed information about the CHMS is available online.23–25
Cycles 2–4 collected information for the substances of interest from people aged 3–79; the overall response rates were 55.5%, 51.7%, and 53.7%, respectively, resulting in 6395, 5785, and 5794 respondents. Participants in one cycle were not eligible to participate in other cycles. The three cycles combined had a total of 4445 respondents aged 3–18. For the present study, 551 respondents were excluded because of missing values for at least one of the PAH metabolites of interest, 21 were then excluded because of missing values for one of the outcomes, and a further 206 were excluded because of missing values for at least one of the covariates, leaving an analytical sample size of 3667. Comparison of the characteristics of in-scope respondents who were not retained (NR) in the analysis to those who were retained (R) in the analysis revealed only minor differences. In the NR group, there was a slightly higher proportion of people aged 12–18 (NR = 51%, R = 45%), they were slightly more likely to be male (NR = 54%, R = 51%) or to eat fruits or vegetables fewer than five times per day (NR = 73%, R = 68%), and had a higher average WC (NR = 69 cm, R = 67 cm).
Measures
PAH metabolites
For the analysis of the PAH metabolites, urine samples were hydrolyzed using β-glucuronidase enzymatic solution and extracted with an organic solvent at neutral pH. The extracts were evaporated, derivatized with N-methyl-N-(trimethylsilyl)-trifluoroacetamide, and analyzed using an Agilent 7890 gas chromatograph coupled with an Agilent 7000B triple-quad tandem mass spectrometer operating in electron impact ionization mode. Analytes were quantified using multiple reaction monitoring. Further information about the laboratory analysis of the PAH metabolites and the quality control measures used in cycles 2–4 of the CHMS is available online.26–28 Limits of detection (LOD) for all metabolites included in this study are presented in Supplementary Table S1. Those below the LOD were assigned a value equal to the highest detection limit of all three cycles, divided by the square root of 2. For this study, the following groups of PAH metabolites were derived: fluorene (FLUO; sum of 2-hydroxyfluorene, 3-hydroxyfluorene, and 9-hydroxyfluorene); naphthalene (NAP; sum of 1-hydroxynaphthalene and 2-hydroxynaphthalene); phenanthrene (PHEN; sum of 1-hydroxyphenanthrene, 2-hydroxyphenanthrene, 3-hydroxyphenanthrene, 4-hydroxyphenanthrene, and 9-hydroxyphenanthrene); and pyrene (PYR; 1-hydroxypyrene). Total PAH was the sum of all groups by their molarity. Quartile ranges were as follows: NAP (µg/l) Q1: ≤2.88; Q2: 2.88–5.14; Q3: 5.14–9.30; Q4: >9.30. FLUO (µg/l) Q1: ≤0.26; Q2: 0.26–0.44; Q3: 0.44–0.70; Q4: >0.70. PHEN (µg/l) Q1: <0.18; Q2: 0.18–0.31; Q3: 0.31–0.52; Q4: >0.52. PYR (µg/l) Q1: ≤0.06; Q2: 0.06–0.11; Q3: 0.11–0.19; Q4: >0.19. Total PAH (nmol/l) Q1: ≤24.52; Q2: 24.52–42.08; Q3: 42.08–73.35; Q4: >73.35. Note that 3-hydroxybenzo(a)pyrene, 2-hydroxychrysene, 3-hydroxychrysene, 4-hydroxychrysene, 6-hydroxychrysene, and 3-hydroxyfluoranthene were also assessed, but as the percent below LOD exceeded 98% for each, they were excluded from the study.
Creatinine
Creatinine was measured in urine using the colorimetric end-point Jaffe method. An alkaline solution of sodium picrate reacts with creatinine in urine to form a red Janovski complex using Microgenics DRI® Creatinine-Detect® reagents (#917). The absorbance was read at 505 nm on a Hitachi 917 chemistry autoanalyzer. The LOD was 4.0 mg/dl in cycle 2 and 5.0 mg/dl in cycles 3 and 4.26–28
Cotinine
Free cotinine was measured in urine and recovered by solid-phase extraction in a 96-well plate format on an automated robotic workstation.29 The LOD was 1 µg/l in cycle 2 and 1.1 µg/l in cycles 3 and 4.26–28
Height
Height was measured to the nearest 0.1 cm using a ProScale M150 digital stadiometer (Accurate Technology Inc., Fletcher, NC, USA).
Weight
Weight was measured to the nearest 0.1 kg with a Mettler Toledo VLC with Panther Plus terminal scale (Mettler Toledo Canada, Mississauga, Canada).
Waist circumference
WC was measured to the nearest 0.1 cm, directly on the landmarked skin with a flexible, inelastic measuring tape with a tension meter attached. In keeping with the NIH protocol (NIH, 2000), the measure was taken at the highest point of the iliac crest.
Waist-to-height ratio
Waist-to-height ratio (WHtR) is a measure of fat distribution and primarily identifies those with abdominal obesity.Reference Browning, Hsieh and Ashwell30,Reference Ashwell, Lejeune and McPherson31 It was calculated as measured WC in centimeters divided by measured height in centimeters.
Body mass index percentile
Body mass index (BMI) was calculated as measured weight in kilograms divided by measured height in meters squared (kg/m2). BMI percentile was derived from the respondent’s BMI-for-age-and-sex Z-score, based on a set of cutoffs specified by the World Health Organization.Reference de Onis, Onyango and Borghi32
Central obesity
Central obesity categories corresponded to: not at risk (ratio of waist to height below 0.50), at risk (ratio of waist to height from 0.50 to less than 0.55), and high probability (ratio of waist to height at or above 0.55).Reference Mehta33
Covariates
We controlled for the following potential confounders of the association between PAHs and obesity: ethnicity, education, income, diet, and exposure to environmental tobacco smoke (ETS). White or nonwhite ethnicity was based on respondents’ answer to which racial or cultural group they belong. Highest level of household education was defined as less than secondary school graduation, secondary school graduation, or higher. Household income (n = 1176 cases were imputed23–25) was adjusted for household size and was categorized according to lowest household income quintile versus above the lowest quintile. Diet quality was captured through reported fruit and vegetable consumptionReference Garriguet34 and was categorized as less than five times per day versus five or more times per day. This was derived from the sum of the frequency of daily consumption of 100% fruit juices; fruit; tomatoes or tomato sauce (excluding tomato paste, ketchup, or pizza sauce); lettuce or green leafy salad; potatoes (including baked, boiled, mashed, or in potato salad, but excluding sweet potatoes); and spinach, mustard greens, or collards (excluding kale); and “other” types of vegetables, excluding those already mentioned. ETS exposure categories were derived from measured urinary cotinine. The not exposed category included those whose cotinine level fell below the LOD. The low, medium, and high ETS exposure categories were based on weighted tertiles among those whose cotinine level was above the LOD.
Statistical analysis
Descriptive statistics were used to examine the characteristics of the study population. Each group of PAH metabolites and total PAH was categorized into quartiles based on the weighted distribution in the study population. Multivariable linear models were run separately to estimate the association between each group of PAHs (in quartiles) and total PAH metabolites (in quartiles), with BMI percentile, WC, and WHtR, respectively, adjusting for age, sex, ethnicity, education, income quintile, diet, creatinine, ETS exposure, and cycle of CHMS. These models were run for the entire study population and stratified by age group (3–5, 6–11, and 12–18 years). Beta coefficients and predicted marginals (referred to as model-adjusted means) were estimated, and the linear trend in these means across each group of PAH metabolites was tested. A multinomial logistic regression model estimated adjusted odds ratios (ORs) for risk of central obesity and high probability of central obesity (each compared with no risk of central obesity) in association with each group of PAH metabolites and total PAH metabolites. Like other studies that have examined PAHs and anthropometric measures,Reference Scinicariello and Buser20,Reference Kim, Kam and Lee21 urinary creatinine was included as an independent variable in all models to account for variation in dilution in spot urinary samples. Because three cycles of data were pooled for this analysis, all models included a CHMS cycle indicator. The CHMS sample design of collection sites selected within regions results in a limited number of degrees of freedom available for variance estimation.23–25 For this study, 35 degrees of freedom were available and specified in all procedure statements. The multinomial regression could not be run separately by age group due to this limited number of degrees of freedom. To further account for the survey’s complex sampling design, all analyses were weighted using the combined CHMS survey weight generated for cycles 2–4, and variance estimation and significance testing based on Satterthwaite statistics were done using the replicate weights generated by Statistics Canada.23–25,35 Confidence intervals (95%) and P-values were presented for all model estimates. The data were analyzed with SAS 9.3 and SAS-Callable SUDAAN 11.0. The code used to conduct analyses is provided in Supplementary Table S2.
Results
Characteristics of the full study population aged 3–18 and according to age group are presented in Table 1. The study population was 51% male, and 85% came from households with higher than secondary school graduation. The proportion of individuals eating fruits and vegetables fewer than five times per day was highest in the 12–18 age group. On the basis of cotinine levels, 83% of the study population was not exposed to ETS, with a higher proportion of those aged 3–5 and 6–11 not exposed compared with those aged 12–18 (86% and 87% vs. 78%). On average, the population was just above the 60th BMI percentile for age and sex, had an average WHtR of 0.47, and 23% were at risk of or had a high probability of central obesity (i.e., WHtR > = 0.50). NAP metabolites accounted for the bulk of total PAH metabolites, and mean levels of all PAH metabolites increased with age.
ETS, environmental tobacco smoke; FLUO, fluorine; LOD, limits of detection; NAP, naphthalene; PAH, polycyclic aromatic hydrocarbon; PHEN, phenanthrene; PYR, pyrene; WCH, waist to height.
Source: Combined 2009–2011, 2012–2013, and 2014–2015 Canadian Health Measures Survey.
PAH means are geomeans.
BMI percentile values are based on WHO growth reference values specific to sex and age in months.
Central obesity categories correspond to: not at risk (ratio of waist to height below 0.50), at risk (ratio of waist to height from 0.50 to less than 0.55), and high probability (ratio of waist to height at or above 0.55).
ETS exposure category based on measured the urinary cotinine values. Not exposed: below the LOD (<1.1 ng/ml), low exposure (1.1–5 ng/ml), medium exposure (>5–24 ng/ml), and high exposure (>24 ng/ml).
We assessed three measures of obesity: two previously assessed for ages 6–19 in relation to urinary PAHs, BMI, and WC,Reference Scinicariello and Buser20,Reference Kim, Kam and Lee21 and one not previously assessed WHtR. After covariate adjustment, BMI percentile for age and sex increased monotonically across quartiles of NAP and total PAH metabolites in the full population of 3–18 year olds (Supplementary Table S3). On average, those in the highest quartile for the NAP metabolites had a predicted BMI percentile of 66th compared with 56th for those in the lowest quartile (Fig. 1). This positive gradient was observed for the 6–11 and 12–18 age groups, separately, but not for those aged 3–5. There was very little difference in BMI percentile across quartiles of FLUO, PHEN, and PYR metabolites.
WC increased monotonically across the NAP metabolites and total PAH metabolites for the total population and for each age group after covariate adjustment (Supplementary Table S4). WC varied little across the other metabolites. On average, there was a positive linear trend (P < 0.01) in predicted WC across the NAP and total PAH metabolites for all age groups (Fig. 2).
Like WC, after covariate adjustment WHtR increased monotonically across quartiles of NAP metabolites and total PAH metabolites for the total population and each age group (Supplementary Table S5) with a corresponding positive linear trend in average predicted WHtR (Fig. 3). Unlike the other outcomes, WHtR also increased across quartiles of PHEN metabolites for those aged 3–5 and those aged 6–11. Among the total population, those in the highest quartile for NAP and total PAH metabolites had three times greater odds of having central obesity versus no central obesity compared with those in the lowest quartile (Table 2). The results also suggest a potential monotonic increase in the odds of having central obesity versus having no risk across the quartiles of PHEN metabolites.
CHMS, Canadian Health Measures Survey; ETS, environmental tobacco smoke; FLUO, fluorine; LOD, limits of detection; NAP, naphthalene; PAH, polycyclic aromatic hydrocarbon; PHEN, phenanthrene; PYR, pyrene.
Source: Combined 2009–2011, 2012–2013, and 2014–2015 Canadian Health Measures Survey.
Not at risk of central obesity = ratio of waist circumference to height <0.50, at risk of central obesity = 0.50≤ ratio of waist circumference to height <0.55, high probability of central obesity = ratio of waist circumference to height ≥0.55.
Models adjusted for age (continuous), sex, education, diet, creatinine, ETS exposure category (two-level), and cycle of CHMS.
ETS exposure category based on the measured urinary cotinine values. Not exposed: below the LOD (<1.1 ng/ml) and exposed (1.1 ng/ml or higher).
Discussion
In utero and early-life exposures may be critical determinants of later onset disease, including metabolic syndrome and obesity.Reference La Merrill and Birnbaum8,Reference De Long and Holloway12,Reference Legler, Hamers and van Eck van der Sluijs-van de Bor36,Reference de Kroon and van de Bor37 As these appear to represent windows of vulnerability to environmental exposures with long-term consequences, there is considerable interest in assessing risks associated with exposure to contaminants to which the population is widely exposed. We show for the first time that urinary PAH metabolites are associated with WHtR, a proxy for central adipose tissue and predictive indicator of health risks associated with obesity, in children as young as 3–5 years old. These associations were independent of exposure to ETS and robust to control for a variety of covariates. The findings build upon previous work reporting associations between urinary PAH metabolites and obesity in children aged 6–19,Reference Poursafa, Dadvand and Amin19–Reference Kim, Kam and Lee21 and contribute to the literature highlighting young children as a population of concern.
Childhood obesity is associated with increased risk of morbidity and mortality in adulthood from cardiovascular disease, hypertension, and type 2 diabetes and carries substantial societal and economic implications.Reference Reinehr5,Reference Grant-Guimaraes, Feinstein, Laber and Kosoy38 Importantly, many of the cardiometabolic risk factors for obesity-related comorbidities (i.e., dyslipidemia, elevated blood pressure, and impaired glucose tolerance) are evident in overweight and obese children.Reference Anderson, Lebovic and Hamilton39,Reference Ali, Cerjak and Kent40 While BMI is a widely accepted measure of obesity, there has been increasing interest in WHtR, a marker of central obesity, including evidence that it may be a better indicator of whole body fat and visceral adipose tissue,Reference Swainson, Batterham, Tsakirides, Rutherford and Hind41 as well as a better predictor of cardiometabolic risk,Reference Browning, Hsieh and Ashwell30,Reference Ashwell, Gunn and Gibson42 when compared to BMI or WC. Indeed, a recent meta-analysis covering 300,000 adults showed that WHtR was a significantly better predictor of adverse cardiometabolic outcomes than WC and BMI in both sexes.Reference Ashwell, Gunn and Gibson42 Results from a number of studies support use of WHtR in children, with some indicating that it outperforms BMI as an indicator of cardiometabolic risk.Reference Nambiar, Truby, Davies and Baxter43–Reference Graves, Garnett and Cowell46 To our knowledge, WHtR has not been previously assessed in relation to urinary PAH metabolites. Our data, which show stronger associations between urinary PAH metabolites and WHtR in children aged 3–5 than were seen for BMI or WC alone, suggest that WHtR may be a more sensitive indicator of potential obesogenic effects in a pediatric population.
Children may be more vulnerable to environmental contaminants for a number of reasons: potential greater exposure because of behavioral patterns and proximity to pollutant sources, greater body burden per unit body mass, and the potential for higher sensitivity due to developmental stage. In a study using data from the National Health and Nutrition Examination Survey (NHANES), significantly higher levels of urinary PAH metabolites were observed in children aged 6–11 compared to 12–19 and adults,Reference Li, Sandau and Romanoff47 consistent with a greater body burden in younger children. Intake of NAP, normalized by subject mass, has been calculated to be over five times higher for a 10-kg child compared to a 70-kg adult.Reference Preuss, Angerer and Drexler48 Effects of exposure to PAHs may start prenatally, as suggested by associations between maternal PAH exposure and childhood obesity.Reference Rundle, Hoepner and Hassoun49 Urinary PAH metabolites are associated with higher levels of markers of inflammation and oxidative stress during pregnancy,Reference Ferguson, McElrath and Pace50 indicating potential systemic impacts of PAH exposure that could impact pregnancy and development. Emerging evidence that PAHs may alter levels of reproductive hormones in umbilical cord seraReference Yin, Tang, Chen, Li and Liu51 provides a further link between maternal exposures and potential prenatal impacts. Thus associations observed in early childhood could possibly reflect effects of both prenatal and postnatal exposures.
In line with previous studies, measures of obesity were most strongly associated with NAP and total PAH metabolites.Reference Scinicariello and Buser20,Reference Kim, Kam and Lee21 NAP metabolites comprise the bulk of total urinary PAH metabolites. Given the abundance of NAP, and the fact that it exists almost entirely in the gas phase, NAP metabolites have been proposed to serve as surrogates for exposure to airborne PAHs.Reference Rappaport, Waidyanatha and Serdar52 Urinary PAH metabolites correspond well with even relatively low levels in ambient air pollution, indicating that they can serve as biomarkers of exposure.Reference Castano-Vinyals, D’Errico, Malats and Kogevinas53 While the majority of significant positive associations were with NAP metabolites, consistent with inhalation exposures, in 3–5 year olds, WHtR was also significantly associated with PHEN and to a lesser extent FLUO, where dietary exposure cannot be discounted.Reference Li, Sandau and Romanoff47,Reference Ramesh, Walker and Hood54 Importantly, the association between urinary PAH metabolites and obesity measures remained after adjusting for exposure to ETS, suggesting a more generalized effect of environmental and dietary exposures. Air pollution has itself been associated with increased prevalence of cardiometabolic risk factors.Reference Yang, Qian and Li55,Reference Poursafa, Mansourian, Motlagh, Ardalan and Kelishadi56
Strengths of the present population-based study include the inclusion of 3–5 year olds, extending the evaluation of associations with PAH metabolites to a younger and possibly more sensitive population with significant exposure to PAHs. We assessed several measures of obesity, including WHtR, which had not previously been assessed in relation to PAH metabolites. The CHMS included a rich covariate dataset enabling adjustment for potential confounders such as diet and ETS exposure. A significant limitation is the cross-sectional nature of the data, which precludes determination of the directionality of effects. Urinary PAH metabolites were only assessed a single time, and thus, it is possible that proximal exposures mask long-term trends from chronic exposure that might be more relevant to long-term biological effects. However, a recent study in New York City showed that repeated measure of urinary PAH metabolites exhibited overall variation of less than 5%.Reference Jung, Liu and Lovinsky-Desir57 As in any population study, other unmeasured factors may influence or drive the association between PAH metabolites and obesity; these include other air pollutants that may co-vary with PAH exposure. Moreover, PAHs that may co-vary with those assessed hereReference Li, Sandau and Romanoff47 but whose metabolites were not present in measurable quantities in urine (e.g., 3-hydroxybenzo(a)pyrene, 2-hydroxychrysene, 3-hydroxychrysene, 4-hydroxychrysene, 6-hydroxychrysene, 3-hydroxyfluoranthene) could nevertheless have biologically relevant effects. Limited sample sizes prevented us from assessing whether ETS exposure level was an effect modifier of the association between PAH metabolites and our measures of obesity. However, analyses conducted in the subpopulation not exposed to ETS (as assessed by cotinine) showed clearly that effects of NAP and total PAH metabolites were not dependent upon ETS exposure (Supplementary Table S6).
Conclusions
In summary, we found that urinary PAH metabolites were associated with measures of obesity in children as young as 3–5 years. While there is considerable interest in evaluating the potential role of early-life exposures in disease, PAHs and air pollution more broadly remain understudied risk factors for later-life obesity and its associated comorbidities. Given the potential importance of such exposures, further research is warranted to examine potential obesogenic and endocrine disrupting effects of PAHs in infants and young children.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S2040174419000825.
Acknowledgements
Thank you to Alain Filiatreault for assistance with preparation of figures.
Financial support
This work was supported by Statistics Canada and Health Canada. Funding sources had no involvement in study design; collection, analysis, and interpretation of data, nor the writing of the publication.
Conflicts of interest
None.