Hostname: page-component-76fb5796d-45l2p Total loading time: 0 Render date: 2024-04-25T14:20:13.972Z Has data issue: false hasContentIssue false

Association of urinary polycyclic aromatic hydrocarbons and obesity in children aged 3–18: Canadian Health Measures Survey 2009–2015

Published online by Cambridge University Press:  06 December 2019

Tracey Bushnik
Affiliation:
Health Analysis Division, Statistics Canada, Ottawa, Canada
Suzy L. Wong
Affiliation:
Health Analysis Division, Statistics Canada, Ottawa, Canada
Alison C. Holloway
Affiliation:
Department of Obstetrics and Gynecology, McMaster University, Hamilton, Canada
Errol M. Thomson*
Affiliation:
Environmental Health Science and Research Bureau, Health Canada, Ottawa, Canada
*
Address for correspondence: Errol Thomson, Environmental Health Science and Research, Bureau, Health Canada, 0802B Tunney’s Pasture, Ottawa, Ontario, K1A 0K9, Canada. Email: errol.thomson@canada.ca
Rights & Permissions [Opens in a new window]

Abstract

Polycyclic aromatic hydrocarbons (PAHs) may contribute to obesity. Childhood obesity is a strong predictor of adult obesity and morbidity; however, the relationship between PAHs and obesity in young children (e.g., aged 3–5) has not been studied. We examined the association between urinary PAH metabolites and measures of obesity in children. We analyzed data from 3667 children aged 3–18 years who participated in the Canadian Health Measures Survey (CHMS, 2009–2015). We ran separate multivariable linear models to estimate the association between quartiles of PAH metabolites and each of body mass index (BMI) percentile, waist circumference (WC), and waist-to-height ratio (WHtR) in the total population, as well as in the age subgroups 3–5, 6–11, and 12–18, adjusting for age, sex, ethnicity, education, income quintile, diet, creatinine, and exposure to environmental tobacco smoke. A multinomial logistic regression model estimated adjusted odds ratios for risk of central obesity. BMI, WC, and WHtR were positively associated with total PAH and naphthalene metabolites in the total population aged 3–18 and in age groups 6–11 and 12–18. In 3–5 year olds, WHtR, but not BMI, was significantly associated with total PAH, naphthalene, and phenanthrene metabolites. Overall, those in the highest quartile for naphthalene or total PAH metabolites had three times greater odds of having central obesity compared with those in the lowest quartile. Urinary PAH metabolites are associated with WHtR, an indicator of central obesity and predictor of health risks associated with obesity, in children as young as 3–5.

Type
Original Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproductionin any medium, provided the original work is properly cited.
Copyright
© Crown Copyright, Statistics Canada and Health Canada 2019. Published with Permission

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 Birnbaum8Reference 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 Amin19Reference 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.2325

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.2628 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.2628

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.2628

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 imputed2325) 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.2325 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.2325,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.

Table 1. Characteristics, by age group, household population aged 3–18, 2007–2015

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.

Fig. 1. Multivariable linear regression association between BMI percentile and quartiles of urinary levels of PAH metabolites, by age group, household population aged 3–18, 2007–2015. Model-adjusted means and 95% confidence intervals are presented for each PAH quartile. The P-value is based on the t-statistic from the test of linear trend in means with a positive gradient across a given set of PAH quartiles.

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).

Fig. 2. Multivariable linear regression association between waist circumference (cm) and quartiles of urinary levels of PAH metabolites, by age group, household population aged 3–18, 2007–2015. Model-adjusted means and 95% confidence intervals are presented for each PAH quartile. The P-value is based on the t-statistic from the test of linear trend in means with a positive gradient across a given set of PAH quartiles.

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.

Fig. 3. Multivariable linear regression association between ratio of waist circumference to height and quartiles of urinary levels of PAH metabolites, by age group, household population aged 3–18, 2007–2015. Model-adjusted means and 95% confidence intervals are presented for each PAH quartile. The P-value is based on the t-statistic from the test of linear trend in means with a positive gradient across a given set of PAH quartiles.

Table 2. Odds ratios (95% CI) from multinomial logistic regression models of association between quartiles of urinary levels of PAH metabolites and risk of central obesity and high probability of central obesity versus not at risk of central obesity, household population aged 3–18, 2007–2015

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 Amin19Reference 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 Baxter43Reference 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.

References

Lehnert, T, Sonntag, D, Konnopka, A, Riedel-Heller, S, Konig, HH. Economic costs of overweight and obesity. Best Pract Res Clin Endocrinol Metab. 2013; 27, 105115.CrossRefGoogle ScholarPubMed
World Health Organisation (WHO). Obesity and Overweight. Fact Sheet, 2018. World Health Organisation, Geneva, Switzerland. http://www.who.int/mediacentre/factsheets/fs311/en/.Google Scholar
NCD Risk Factor Collaboration (NCD-RisC). Worldwide trends in body-mass index, underweight, and obesity from 1975 to 2016: a pooled analysis of 2416 population-based measurement studies in 128-9 million children, adolescents, and adults.NCD Risk Factor Collaboration (NCD-RisC). Lancet. 2017; 390(10113), 26272642. doi: 10.1016/S0140-6736(17)32129-3. Epub. 2017.CrossRefGoogle Scholar
Rankin, J, Matthews, L, Cobley, S, et al. Psychological consequences of childhood obesity: psychiatric comorbidity and prevention. Adolesc Health Med Ther. 2016; 7, 125146.CrossRefGoogle ScholarPubMed
Reinehr, T. Long-term effects of adolescent obesity: time to act. Nat Rev Endocrinol. 2018; 14, 183188.CrossRefGoogle Scholar
Simmonds, M, Burch, J, Llewellyn, A, et al. The use of measures of obesity in childhood for predicting obesity and the development of obesity-related diseases in adulthood: a systematic review and meta-analysis. Health Technol Assess. 2015; 19, 1336.CrossRefGoogle ScholarPubMed
Heindel, JJ, Blumberg, B. Environmental obesogens: mechanisms and controversies. Annu Rev Pharmacol Toxicol. 2019; 59, 89106.CrossRefGoogle ScholarPubMed
La Merrill, M, Birnbaum, LS. Childhood obesity and environmental chemicals. Mt Sinai J Med. 2011; 78, 2248.CrossRefGoogle ScholarPubMed
Hu, H, Kan, H, Kearney, GD, Xu, X. Associations between exposure to polycyclic aromatic hydrocarbons and glucose homeostasis as well as metabolic syndrome in nondiabetic adults. Sci Total Environ. 2015; 505, 5664.CrossRefGoogle ScholarPubMed
Ranciere, F, Lyons, JG, Loh, VH, et al. Bisphenol A and the risk of cardiometabolic disorders: a systematic review with meta-analysis of the epidemiological evidence. Environ Health. 2015; 14, 46.CrossRefGoogle ScholarPubMed
James-Todd, TM, Huang, T, Seely, EW, Saxena, AR. The association between phthalates and metabolic syndrome: the National health and nutrition examination Survey 2001-2010. Environ Health. 2016; 15, 52.Google ScholarPubMed
De Long, NE, Holloway, AC. Early-life chemical exposures and risk of metabolic syndrome. Diabetes Metab Syndr Obes. 2017; 10, 101109.CrossRefGoogle ScholarPubMed
Schultz, TW, Sinks, GD. Xenoestrogenic gene expression: structural features of active polycyclic aromatic hydrocarbons. Environ Toxicol Chem. 2002; 21, 783786.CrossRefGoogle ScholarPubMed
Irigaray, P, Ogier, V, Jacquenet, S, et al. Benzo[a]pyrene impairs beta-adrenergic stimulation of adipose tissue lipolysis and causes weight gain in mice. A novel molecular mechanism of toxicity for a common food pollutant. FEBS J. 2006; 273, 13621372.CrossRefGoogle ScholarPubMed
Stephanou, P, Konstandi, M, Pappas, P, Marselos, M. Alterations in central monoaminergic neurotransmission induced by polycyclic aromatic hydrocarbons in rats. Eur J Drug Metab Pharmacokinet. 1998; 23, 475481.Google ScholarPubMed
Irigaray, P, Lacomme, S, Mejean, L, Belpomme, D. Ex vivo study of incorporation into adipocytes and lipolysis-inhibition effect of polycyclic aromatic hydrocarbons. Toxicol Lett. 2009; 187, 3539.CrossRefGoogle ScholarPubMed
Kim, JH, Yamaguchi, K, Lee, SH, et al. Evaluation of polycyclic aromatic hydrocarbons in the activation of early growth response-1 and peroxisome proliferator activated receptors. Toxicol Sci. 2005; 85, 585593.CrossRefGoogle ScholarPubMed
Grun, F, Blumberg, B. Environmental obesogens: organotins and endocrine disruption via nuclear receptor signaling. Endocrinology. 2006; 147, S50S55.CrossRefGoogle ScholarPubMed
Poursafa, P, Dadvand, P, Amin, MM, et al. Association of polycyclic aromatic hydrocarbons with cardiometabolic risk factors and obesity in children. Environ Int. 2018; 118, 203210.CrossRefGoogle ScholarPubMed
Scinicariello, F, Buser, MC. Urinary polycyclic aromatic hydrocarbons and childhood obesity: NHANES (2001-2006). Environ Health Perspect. 2014; 122, 299303.CrossRefGoogle Scholar
Kim, HW, Kam, S, Lee, DH. Synergistic interaction between polycyclic aromatic hydrocarbons and environmental tobacco smoke on the risk of obesity in children and adolescents: the U.S. National health and nutrition examination Survey 2003–2008. Environ Res. 2014; 135, 354360.CrossRefGoogle ScholarPubMed
de Kroon, ML, Renders, CM, van Wouwe, JP, van, BS, Hirasing, RA. The terneuzen birth cohort: BMI changes between 2 and 6 years correlate strongest with adult overweight. PLoS One. 2010; 5, e9155.CrossRefGoogle ScholarPubMed
Statistics Canada. Canadian Health Measures Survey (CHMS) Data User Guide: Cycle 2, November 2012. Available at http://www23.statcan.gc.ca/imdb-bmdi/documents/5071_D4_T9_V1-eng.htm, 2012.Google Scholar
Statistics Canada. Canadian Health Measures Survey (CHMS) Data User Guide: Cycle 3, November 2014. Available at http://www.statcan.gc.ca, 2014.Google Scholar
Statistics Canada. Canadian Health Measures Survey (CHMS) Data User Guide: Cycle 4, October 2017. Available at http://www.statcan.gc.ca, 2017.Google Scholar
Health Canada. Second Report on Human Biomonitoring of Environmental Chemicals in Canada:Results of the Canadian Health Measures Survey Cycle 2 (2009-2011) (Health Canada, Catalogue H128-1/10/601-1E-PDF) Ottawa: Minister of Health, 2013.Google Scholar
Health Canada. Third Report on Human Biomonitoring of Environmental Chemicals in Canada:Results of the Canadian Health Measures Survey Cycle 3 (2012-2013) (Health Canada, Catalogue H128-1/10/601-2E-PDF) Ottawa: Minister of Health, 2015.Google Scholar
Health Canada. Fourth Report on Human Biomonitoring of Environmental Chemicals in Canada:Results of the Canadian Health Measures Survey Cycle 4 (2014-2015) (Health Canada, Catalogue H128-1/10/601-4E-PDF) Ottawa: Minister of Health, 2017.Google Scholar
Institut national de santé publique du Québec. Analytical Method for the Determination of Cotinine in Urine by HPLC-MS-MS Robotic Workstation Method (C-550), Condensed version for CHMS, 2009. Québec:Institut national de santé publique du Québec.Google Scholar
Browning, LM, Hsieh, SD, Ashwell, M. A systematic review of waist-to-height ratio as a screening tool for the prediction of cardiovascular disease and diabetes: 0.5 could be a suitable global boundary value. Nutr Res Rev. 2010; 23, 247269.CrossRefGoogle ScholarPubMed
Ashwell, M, Lejeune, S, McPherson, K. Ratio of waist circumference to height may be better indicator of need for weight management. BMJ. 1996; 312, 377.CrossRefGoogle ScholarPubMed
de Onis, M, Onyango, AW, Borghi, E, et al. Development of a WHO growth reference for school-aged children and adolescents. Bull World Health Organ. 2007; 85, 660667.CrossRefGoogle ScholarPubMed
Mehta, SK. Waist circumference to height ratio in children and adolescents. Clin Pediatr (Phila). 2015; 54, 652658.CrossRefGoogle ScholarPubMed
Garriguet, D. Diet quality in Canada. Health Rep. 20(3); 2009, 4152.Google ScholarPubMed
Statistics Canada. Instructions for combining multiple cycles of Canadian Health Measures Survey (CHMS) data, 2017. Available at http://www.statcan.gc.ca.Google Scholar
Legler, J, Hamers, T, van Eck van der Sluijs-van de Bor, M. et al. The OBELIX project: early life exposure to endocrine disruptors and obesity. Am J Clin Nutr. 2011; 94, 1933S1938S.CrossRefGoogle ScholarPubMed
de Kroon, ML, van de Bor, M. Obesogenic effects of endocrine disruptors, what do we know from animal and human studies? Environ Int. 2014; 70, 1524.Google Scholar
Grant-Guimaraes, J, Feinstein, R, Laber, E, Kosoy, J. Childhood overweight and obesity. Gastroenterol Clin North Am. 2016; 45, 715728.CrossRefGoogle ScholarPubMed
Anderson, LN, Lebovic, G, Hamilton, J, et al. Body mass index, waist circumference, and the clustering of cardiometabolic risk factors in early childhood. Paediatr Perinat Epidemiol. 2016; 30, 160170.CrossRefGoogle ScholarPubMed
Ali, O, Cerjak, D, Kent, JW, et al. Obesity, central adiposity and cardiometabolic risk factors in children and adolescents: a family-based study. Pediatr Obes. 2014; 9, e58e62.CrossRefGoogle ScholarPubMed
Swainson, MG, Batterham, AM, Tsakirides, C, Rutherford, ZH, Hind, K. Prediction of whole-body fat percentage and visceral adipose tissue mass from five anthropometric variables. PLoS One. 2017; 12, e0177175.CrossRefGoogle ScholarPubMed
Ashwell, M, Gunn, P, Gibson, S. Waist-to-height ratio is a better screening tool than waist circumference and BMI for adult cardiometabolic risk factors: systematic review and meta-analysis. Obes Rev. 2012; 13, 275286.CrossRefGoogle ScholarPubMed
Nambiar, S, Truby, H, Davies, PS, Baxter, K. Use of the waist-height ratio to predict metabolic syndrome in obese children and adolescents. J Paediatr Child Health. 2013; 49, E281E287.CrossRefGoogle ScholarPubMed
Khoury, M, Manlhiot, C, McCrindle, BW. Role of the waist/height ratio in the cardiometabolic risk assessment of children classified by body mass index. J Am Coll Cardiol. 2013; 62, 742751.CrossRefGoogle ScholarPubMed
Mokha, JS, Srinivasan, SR, Dasmahapatra, P, et al. Utility of waist-to-height ratio in assessing the status of central obesity and related cardiometabolic risk profile among normal weight and overweight/obese children: the Bogalusa Heart Study. BMC Pediatr. 2010; 10, 73.CrossRefGoogle ScholarPubMed
Graves, L, Garnett, SP, Cowell, CT, et al. Waist-to-height ratio and cardiometabolic risk factors in adolescence: findings from a prospective birth cohort. Pediatr Obes. 2014; 9, 327338.CrossRefGoogle ScholarPubMed
Li, Z, Sandau, CD, Romanoff, LC, et al. Concentration and profile of 22 urinary polycyclic aromatic hydrocarbon metabolites in the US population. Environ Res. 2008; 107, 320331.CrossRefGoogle ScholarPubMed
Preuss, R, Angerer, J, Drexler, H. Naphthalene–an environmental and occupational toxicant. Int Arch Occup Environ Health. 2003; 76, 556576.CrossRefGoogle ScholarPubMed
Rundle, A, Hoepner, L, Hassoun, A, et al. Association of childhood obesity with maternal exposure to ambient air polycyclic aromatic hydrocarbons during pregnancy. Am J Epidemiol. 2012; 175, 11631172.CrossRefGoogle ScholarPubMed
Ferguson, KK, McElrath, TF, Pace, GG, et al. Urinary polycyclic aromatic hydrocarbon metabolite associations with biomarkers of inflammation, angiogenesis, and oxidative stress in pregnant women. Environ Sci Technol. 2017; 51, 46524660.CrossRefGoogle Scholar
Yin, S, Tang, M, Chen, F, Li, T, Liu, W. Environmental exposure to polycyclic aromatic hydrocarbons (PAHs): the correlation with and impact on reproductive hormones in umbilical cord serum. Environ Pollut. 2017; 220, 14291437.CrossRefGoogle ScholarPubMed
Rappaport, SM, Waidyanatha, S, Serdar, B. Naphthalene and its biomarkers as measures of occupational exposure to polycyclic aromatic hydrocarbons. J Environ Monit. 2004; 6, 413416.CrossRefGoogle ScholarPubMed
Castano-Vinyals, G, D’Errico, A, Malats, N, Kogevinas, M. Biomarkers of exposure to polycyclic aromatic hydrocarbons from environmental air pollution. Occup Environ Med. 2004; 61, e12.CrossRefGoogle ScholarPubMed
Ramesh, A, Walker, SA, Hood, DB, et al. Bioavailability and risk assessment of orally ingested polycyclic aromatic hydrocarbons. Int J Toxicol. 2004; 23, 301333.CrossRefGoogle ScholarPubMed
Yang, BY, Qian, ZM, Li, S, et al. Long-term exposure to ambient air pollution (including PM1) and metabolic syndrome: the 33 Communities Chinese Health Study (33CCHS). Environ Res. 2018; 164, 204211.CrossRefGoogle Scholar
Poursafa, P, Mansourian, M, Motlagh, ME, Ardalan, G, Kelishadi, R. Is air quality index associated with cardiometabolic risk factors in adolescents? The CASPIAN-III Study. Environ Res. 2014; 134, 105109.CrossRefGoogle ScholarPubMed
Jung, KH, Liu, B, Lovinsky-Desir, S, et al. Time trends of polycyclic aromatic hydrocarbon exposure in New York City from 2001 to 2012: assessed by repeat air and urine samples. Environ Res. 2014; 131, 95103.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Characteristics, by age group, household population aged 3–18, 2007–2015

Figure 1

Fig. 1. Multivariable linear regression association between BMI percentile and quartiles of urinary levels of PAH metabolites, by age group, household population aged 3–18, 2007–2015. Model-adjusted means and 95% confidence intervals are presented for each PAH quartile. The P-value is based on the t-statistic from the test of linear trend in means with a positive gradient across a given set of PAH quartiles.

Figure 2

Fig. 2. Multivariable linear regression association between waist circumference (cm) and quartiles of urinary levels of PAH metabolites, by age group, household population aged 3–18, 2007–2015. Model-adjusted means and 95% confidence intervals are presented for each PAH quartile. The P-value is based on the t-statistic from the test of linear trend in means with a positive gradient across a given set of PAH quartiles.

Figure 3

Fig. 3. Multivariable linear regression association between ratio of waist circumference to height and quartiles of urinary levels of PAH metabolites, by age group, household population aged 3–18, 2007–2015. Model-adjusted means and 95% confidence intervals are presented for each PAH quartile. The P-value is based on the t-statistic from the test of linear trend in means with a positive gradient across a given set of PAH quartiles.

Figure 4

Table 2. Odds ratios (95% CI) from multinomial logistic regression models of association between quartiles of urinary levels of PAH metabolites and risk of central obesity and high probability of central obesity versus not at risk of central obesity, household population aged 3–18, 2007–2015

Supplementary material: File

Bushnik et al. supplementary material

Bushnik et al. supplementary material 1

Download Bushnik et al. supplementary material(File)
File 74.8 KB