Skip to main content
Download PDF
Download ePub

Urinary polycyclic aromatic hydrocarbon excretion and regional body fat distribution: evidence from the U.S. National Health and Nutrition Examination Survey 2001–2016

Environmental Health volume 21, Article number: 75 (2022) Cite this article

Abstract

Background

Polycyclic aromatic hydrocarbons (PAHs) are environmental pollutants that may contribute to the etiology of obesity. However, it is unclear whether PAHs from environmental sources are associated with regional body fat distribution, and whether the association varies across racial/ethnic groups who may have differential PAH exposure patterns.

Objectives

To examine correlations between PAHs and body fat distribution, and potential racial/ethnic differences among U.S. adults.

Methods

Ten PAHs were measured in spot urine samples from 2691 non-smoking adults (age ≥ 20 years) in the NHANES 2001–2016. Dual-energy X-ray absorptiometry was used to measure fat mass percent (FM%). Partial Pearson correlation coefficients (r) with multivariable adjustment were used to assess PAH-FM% associations.

Results

In the total population, 1-naphthalene, 3-fluorene, and 1-pyrene were inversely correlated with total FM% or trunk FM% (adjusted r ranged: − 0.06 to − 0.08), while 2-naphthalene, 9-fluorene, and 4-phenanthrene were positively correlated with the FM% measurements (r: 0.07–0.11). PAH levels are highest among non-Hispanic Blacks, followed by Hispanics and Whites and some of the correlations were different by these races/ethnicities. Among non-Hispanic Whites, no PAH was correlated with FM%. In contrast, 9-fluorene was positively correlated with total FM% (r = 0.20) and trunk FM% (r = 0.22) among Blacks, and 4-phenanthrene was positively correlated with total FM% (r = 0.23) and trunk FM% (r = 0.24) among Hispanics (P-interaction: 0.010–0.025).

Discussion

In this US adult population, certain PAHs are significantly associated with higher body fat contents among non-Hispanic Blacks and Hispanics but not non-Hispanic Whites, suggesting that minority groups might be particularly susceptible to PAH’s obesogenic effects or the effects of other factors that determine the PAH exposure levels. Alternatively, differences in body composition may contribute to differential PAH metabolism in minority groups. Future studies are warranted to explore the racial/ethnic disparity in PAH exposures, drivers of these exposure differences, and mechanisms through which PAHs may influence body composition by races/ethnicities.

Peer Review reports

Introduction

Polycyclic aromatic hydrocarbons (PAHs) are environmental pollutants that may exert detrimental effects on human health [1]. PAHs are generated through incomplete combustion of organic materials (e.g., coal, oil, gas, wood, garbage, and tobacco) [2]. Common nonoccupational routes of PAH exposures include inhalation (e.g., air contaminated with motor vehicle exhaust, cigarette smoke, wood smoke), diet (e.g., grilled/charred meats or foods, contaminated flour and bread products, processed and pickled foods, and contaminated water and cow’s milk), and skin absorption (e.g., soil and dust) [2, 3]. After absorption, PAHs are transported by chylomicrons and stored mostly in adipose tissue, kidney, and liver [3, 4].

PAHs are potential carcinogens [2], and emerging evidence has suggested that exposure to PAHs may also contribute to the development of metabolic disorders, including diabetes [5, 6], dyslipidemia [6], hypertension [6, 7], and cardiovascular disease [7]. Hydroxylated PAHs, the main metabolic products of PAHs, are endocrine disruptors that mimic estrogens and inhibit thyroid and androgen functions [8,9,10,11]. In an animal study, PAH exposure impaired adipose tissue lipolysis and resulted in increased body weight and fat mass in mice [12]. Epidemiological studies among children and adults also showed that urinary excretion of PAHs was positively associated with higher levels of body mass index (BMI), waist circumference (WC), or waist-to-height ratio [6, 13,14,15,16,17]. However, evidence is limited regarding the associations between PAH exposures and body fat distribution in adults.

Body fat distribution has been recognized to play a significant role in the etiology of metabolic diseases [18]. There is mounting evidence that the biological functions of adipose tissue may differ by locations [18,19,20]. To the best of our knowledge, no study has assessed the associations between individual PAHs and various fat depots (e.g., trunk fat, leg fat, lean mass). In addition, it is known that fat deposition and body fat distribution may differ by race/ethnicity [21]. A recent study conducted among U.S. non-smokers observed the highest levels of urinary PAHs among non-Hispanic Blacks [22], highlighting an environmental exposure disparity in PAHs that could contribute to possible differences in related health outcomes. However, given these exposure differences, no study has examined whether the associations between PAHs and body fatness vary by racial/ethnic groups.

To fill these knowledge gaps, we examined the direction and strength of correlations of urinary PAHs with different body anthropometric measurements, including fat mass (total fat, trunk fat, leg fat), lean mass (bone mineral density, total lean mass), BMI, and WC among a large sample of U.S. adults from the National Health and Nutrition Examination Survey (NHANES) 2001–2016. We explicitly compared correlations of PAHs among major racial/ethnic groups (non-Hispanic Whites, non-Hispanic Blacks, and Hispanics) and hypothesized that PAHs are differentially associated with body fat distribution by race/ethnicity.

Methods

Study population

The NHANES is a nation-wide representative survey assessing the health and nutritional status of residents in the United States. The detailed study design has been described previously [23]. Briefly, for each biennial NHANES survey a complex sampling process was used to randomly select U.S. residents who are representative of the civilian non-institutionalized population. Participants completed interviews at their homes with trained health professionals and underwent physical examinations in a mobile examination center. During the physical examination, a single spot urine sample was collected from NHANES participants aged six years and older. Written informed consents were obtained from all participants. The study protocol was approved by the institutional review board at the Centers for Disease Control and Prevention (Atlanta, Georgia) [23].

The current study combined data from six NHANES cycles (2001–2002, 2003–2004, 2005–2006, 2011–2012, 2013–2014, 2015–2016) among participants ≥20 years old. We excluded two cycles (2007–2008, 2009–2010) because the measurements on body fat distribution were unavailable. Since cigarette smoking is a predominant exposure source of PAHs among smokers, and smoking is known to modulate body weight [24], the association between PAHs and body fat distribution can be largely driven by smoking status. To avoid the strong confounding of cigarette smoking on the associations between PAHs and body fat contents, we restricted the current analysis to non-smokers. In addition, participants below 20 years old (n = 27,842), without medical examination data (n = 2284) or complete DXA scan (n = 10,214), with medical conditions that may influence body fat distribution (kidney disease, physical limitations, diabetes, cardiovascular disease, asthma, pulmonary disease, and cancer) (n = 2504), and with missing values in smoking status (n = 3) were excluded. Of the 61,049 surveyed residents, 7531 individuals had measurements of both PAH levels and body fat distribution. Of them, 2691 non-smokers were included in the main analysis after exclusions. The flowchart of the study design is shown in Supplemental Fig. S1.

Measurements of urinary PAH metabolites

In NHANES study, urinary monohydroxylated metabolites of PAHs were measured in approximately one-third subsample of all participants 6 years and older in the six NHANES cycles. Ten urinary PAHs were measured, including two naphthalene metabolites (1-napthalene, 2-napthalene), three fluorene metabolites (2-fluorene, 3-fluorene, 9-fluorene), four phenanthrene metabolites (1-phenanthrene, 2-phenanthrene, 3-phenanthrene, 4-phenanthrene), and one pyrene metabolite (1-pyrene). In addition to the ten individual PAHs, we created four other variables: the sum of the molar mass of all PAH metabolites, total naphthalene metabolites, total phenanthrene metabolites, and total phenanthrene metabolites, respectively. Of note, the number of PAHs measured in each NHANES cycle varied: it was nine in 2001–2002, ten in 2003–2004/2005–2006/2011–2012, and six in 2013–2014/2015–2016 cycles. We included PAHs measured in at least two survey cycles to ensure a reasonable sample size.

The detailed procedure of laboratory measurement has been described previously [25]. In brief, PAHs were measured in the single spot urine sample by capillary gas chromatography and high-resolution mass spectrometry. To minimize the variability of urinary PAH levels due to differential dilutions, creatinine-adjusted PAHs (dividing PAH concentrations [nanograms per liter urine] by creatinine concentrations [grams per liter urine]) were used in the current analysis as recommended previously [26]. Most PAHs had a detection rate of > 99%, and the detection rate for 4-phenanthrene and 9-fluorene was 89 and 93%, respectively. The coefficient of variation (CV) of PAH measurements ranged from 5 to 13%.

Dual-energy X-ray absorptiometry (DXA) measurements

DXA scans were administered to eligible survey participants 8 years of age and older in the NHANES mobile examination centers. Females with a positive pregnancy test or those who reported being pregnant at the time of the exam were excluded from the DXA examination. Individuals who reported taking tests with radiographic contrast material in the past 72 hours, participants in nuclear medicine studies in the past 3 days, or those who had a self-reported weight (> 300 pounds) or height (> 6′5″) over the DXA table limit were also excluded from the DXA examination [27]. The whole body DXA was taken with a Hologic QDR-4500A fan-beam densitometer (Hologic, Inc., Bedford, Massachusetts) and provided fat mass and lean mass measurements for the total body, both arms, both legs, the trunk, and the head. Hologic Discovery software was used to derive fat mass and lean mass including total fat, trunk fat, leg fat, bone mineral density, and total lean mass. Multiple imputations were conducted to impute missing data for attenuating any potential biases due to missing DXA data. Specifically, missing readings for DXA data were imputed five times by the sequential regression multivariate imputation [27]. Fat mass percentage (FM%) for total body, trunk and legs were calculated as: FM% = fat mass (kg)/total mass (kg) × 100%.

Measurements of body mass index and waist circumference

Body measurements such as weight (kg), height (cm), and WC (cm) were measured by trained health professionals following standardized protocols [28]. BMI was calculated as weight divided by height squared (kg/m2).

Covariates

Information on demographics, lifestyle, medical history, and dietary intake was collected using survey questionnaires during the face-to-face interview [29]. Demographic information included age (years), gender (men, women), race/ethnicity (non-Hispanic White, non-Hispanic Black, Hispanic, or others), education (high school or below, any college, and college graduate or above), and poverty income ratio (< and ≥ 1). The poverty income ratio was used as an indicator of socioeconomic status and was calculated by dividing family income by the poverty threshold adjusted for family size and inflation. Lifestyle factors consisted of alcohol use (nondrinkers, 1–3 drinks/day, and ≥ 4 drinks/day) and moderate-to-vigorous physical activity (yes, no). Dietary intake included total calorie intake (kcal/day) and protein intake (gram/day). Serum cotinine levels (ng/mL), an indicator for environmental tobacco smoke among non-smokers [30], were measured by the isotope-dilution high-performance liquid chromatography/atmospheric pressure chemical ionization tandem mass spectrometric method with a lower detection limit of 0.015 ng/ml. C-reactive protein (CRP) levels were quantified by latex-enhanced nephelometry, with a lower detection limit of 0.1 mg/L.

Statistical analysis

Since the number of missing values varies among PAHs, sample size may vary for individual PAHs. Log-transformation (e base) was applied for creatinine-adjusted PAHs and FM% (total fat, trunk fat, leg fat, trunk/leg ratio), bone mineral density, total lean mass, BMI, and WC to improve the normality. In this cross-sectional analysis, we calculated partial Pearson correlation coefficients (r) weighted by the NHANES sampling weight with adjustment of potential confounding factors to examine correlations between log-transformed PAHs and FM%, bone mineral density, total lean mass, BMI, and WC among non-smokers. Covariates adjusted included age, gender, race/ethnicity, education, poverty income ratio, moderate-to-vigorous physical activity, alcohol use, total calorie intake, protein intake, and serum cotinine and CRP levels. Missing values of continuous covariates were replaced with median values. For categorical variables, a missing indicator variable was created for missing values. To examine the potential differences by race/ethnicity, we performed stratified analyses by main racial/ethnic groups (non-Hispanic White, non-Hispanic Black, and Hispanic). Correlation coefficients of creatinine-adjusted PAHs with FM% were transformed to Fisher z scores and compared using the Student’s t-test between two groups. We also compared correlations of PAHs with trunk FM% versus leg FM% using the Wolfe’s method [31, 32]. Due to the multiple-imputation procedure, each participant has five sets of measured and imputed values of body fat; correlations were calculated within each DXA dataset, and combined into a single composite estimate using the method of Rubin and Schenker [33].

To account for multiple comparisons, False Discovery Rate (FDR) Benjamini-Hochberg Procedure was used for the correction of P values. Statistical significance was determined by a two-sided FDR P value smaller than 0.05. All data were analyzed in RStudio version 1.4.1106 (Rstudio, PBC) and SAS version 9.4 (SAS Institute, Cary, NC).

Results

Characteristics of non-smoking adults stratified by race/ethnicity are shown in Table 1. In the total population, the median age was 37.9 years (interquartile range: 28.0–48.2) and 47.4% were men. In addition, 35.6% were non-Hispanic Whites, 21.7% were non-Hispanic Blacks, and 28.5% were Hispanics. Compared with non-Hispanic Whites, non-Hispanic Blacks and Hispanics were younger and had lower education and lower income. In addition, non-Hispanic Blacks had the highest urinary levels for all ten PAHs compared to non-Hispanic whites and Hispanics. The levels of three out of ten PAHs (i.e., 2-naphthalene, 2-fluorene, and 1-pyrene) were higher among Hispanics compared to non-Hispanic Whites; the remaining seven PAHs had similar levels between non-Hispanic Whites and Hispanics. Non-Hispanic Blacks and Hispanics had higher BMI and leg FM% than non-Hispanic Whites. In addition, Hispanics had higher total FM% and trunk FM% and lower total lean mass compared to non-Hispanic Whites and non-Hispanic Blacks (Table 1).

Table 1 Characteristics of study population among non-smokers aged 20 years stratified by race/ethnicity
Full size table

The adjusted correlation coefficients between the log-transformed PAHs ranged between 0.40 and 0.92 (Supplemental Table S1). The adjusted correlations of PAHs in the same class were comparable with those for PAHs in different classes. Fluorene metabolites had higher correlations with phenanthrene metabolites and pyrene metabolites compared with naphthalene metabolites.

Four out of ten PAHs were significantly correlated with total FM% after the FDR correction: 1-naphthalene, 3-fluorene, and 1-pyrene were inversely correlated with higher total FM% (adjusted r ranged from − 0.06 to − 0.07) and 4-phenanthrene was positively correlated with total FM% (adjusted r = 0.09; all FDR Ps < 0.05) (Table 2). In addition, six out of ten PAHs were significantly correlated with trunk FM% (Table 2). While 1-naphthalene, 3-fluorene, and 1-pyrene were negatively correlated with trunk FM% (adjusted r ranged from − 0.06 to − 0.08), 2-naphthalene, 9-fluorene, and 4-phenanthrene were positively correlated with trunk FM% (adjusted r ranged from 0.07 to 0.11; all FDR Ps < 0.05). These PAHs were generally correlated with the trunk/leg ratio, BMI, and WC, but were not correlated with leg FM%, bone mineral density and total lean mass. Five out of six correlation coefficients of PAHs that were significantly correlated with trunk FM% were statistically stronger than correlation coefficients for leg FM%. Total naphthalene metabolites, total phenanthrene metabolites, total phenanthrene metabolites, and total PAHs were not correlated with body fat distribution.

Table 2 Partial Pearsona correlation between creatinine-adjusted PAHs and body fat percentages among non-smokers aged 20 years and older
Full size table

The strength of PAH-FM% associations was statistically significantly different among three racial/ethnic groups (Tables 3, 4 and 5 & Fig. 1). Specifically, among non-Hispanic Whites, no PAH was correlated with total FM%, trunk FM%, or leg FM% (Table 3). Among non-Hispanic Blacks, 9-fluorene was significantly correlated with higher total FM% (adjusted r = 0.20) and trunk FM% (adjusted r = 0.22; both FDR Ps < 0.05) (Table 4). Among Hispanics, 4-phenanthrene was significantly correlated with higher total FM% (adjusted r = 0.23) and trunk FM% (adjusted r = 0.24; both FDR Ps < 0.05) (Table 5). The concentrations (geometric mean [standard error]) of 9-fluorene and 4-phenanthrene were higher among non-Hispanic Blacks (341.4 [19.3] ng/L and 28.1 [1.56] ng/L) compared with Hispanics (223.9 [14.8] ng/L and 22.6 [1.44] ng/L) and non-Hispanic Whites (232.4 [13.2] ng/L and 21.0 [1.11] ng/L) (Table 1). When compared across racial/ethnic groups, the correlations of 9-fluorene and 4-phenanthrene with total FM% and trunk FM% were statistically significantly stronger among non-Hispanic Blacks and Hispanics than those of non-Hispanic Whites (P for interaction ranged from 0.010 to 0.025) (Tables 3, 4 and 5 & Fig. 1).

Table 3 Partial Pearsona correlation between creatinine-adjusted PAHs and body fat percentages among non-smokers with non-Hispanic White ethnicity aged 20 years and older
Full size table
Table 4 Partial Pearsona correlation between creatinine-adjusted PAHs and body fat percentages among non-smokers with non-Hispanic Black ethnicity aged 20 years and older
Full size table
Table 5 Partial Pearsona correlation between creatinine-adjusted PAHs and body fat percentages among non-smokers with Hispanic ethnicity aged 20 years and older
Full size table
Fig. 1
figure 1

Correlations between selected PAHs that were significantly associated with body fat distribution stratified by race/ethnicity among non-smokers

Full size image

Discussion

In this large sample of US non-smoking adults, we found various correlation patterns between urinary PAH levels and regional body fat distribution. While 1-naphthalene, 3-fluorene, and 1-pyrene were inversely correlated with higher body fat contents, we observed positive correlations for 2-naphthalene, 9-fluorene, and 4-phenanthrene with the same body composition variables. Relevant to exposure disparities, higher PAH concentrations were observed for non-Hispanic Blacks and Hispanics compared to non-Hispanic Whites. Perhaps linked to these higher PAH levels, we found that positive correlations of 9-fluorene and 4-phenanthrene with total FM% and trunk FM% were statistically significantly stronger among non-Hispanic Blacks and Hispanics compared to non-Hispanic Whites. Of PAHs that were associated with body fat, five out of six PAHs had significantly stronger correlations with trunk FM% compared to leg FM%.

To the best of our knowledge, the current analysis is among the first to assess correlations between individual PAHs and body fat distribution among non-smokers, adjusting for a variety of potential confounders. Prior studies in adults evaluating individual PAHs and obesity were scarce [6, 16, 17]. In addition, these studies included both smokers and non-smokers, which can impact the association between PAHs and body fat distribution [24]. Also, previous studies used BMI as the body fatness measure, which could not differentiate between body fat and muscle weight, nor capture visceral fat [34]. In an analysis based on the NHANES 2013–2014 data, total PAHs were not correlated with BMI among adults aged 20–80 years old [16]. A study based on the Korean National Environmental Health Survey found that 2-naphthalene and total PAH were positively correlated with BMI among 3787 adults [17]. In another study among 4765 adults from the 2001–2008 NHANES study, 2-phenanthrene had a monotonic dose-response relationship with higher BMI [6], while 1-naphthalene was inversely correlated with BMI [6]. In the current analysis, we observed the positive correlation of 2-naphthalene and body fat as well as the negative correlation of 1-naphthalene with body fat, but not the positive correlation between 2-phenanthrene and body fat. The differences in our results compared to others might be due to the fact that our study only included non-smokers, where the PAH concentrations tend to be lower and thus the association may differ from the prior study that included smokers [6], who have higher levels of PAHs.

In addition to overall body fat contents, we also examined regional fat distribution. Interestingly, we observed consistently stronger correlations of PAHs with trunk fat compared with leg fat, which may suggest that PAHs impact trunk fat and leg fat differently through differed pathophysiological pathways. Compared with the lower-body fat, the upper-body fat is more active on fat metabolism [35] and has higher exchange of blood flow [36], and thus may favor the storage in and the release of PAHs from fat in this region. Alternatively, the observed differences in the correlations between PAHs and body fat distribution may suggest stronger influences of PAHs on fat metabolism in the trunk area. Similarly, other endocrine-disrupting chemicals have also shown stronger correlations with trunk fat compared with leg fat [32]. In addition to obesity, PAHs are also associated with other cardiometabolic diseases such as diabetes [5, 6], dyslipidemia [6], hypertension [6, 7], and cardiovascular disease [7]. Furthermore, the accumulation of trunk fat, but not leg fat, is associated with increased metabolic diseases [18,19,20]. Therefore, our results are in line with the notion that PAH exposures may potentially link to cardiometabolic diseases through the accumulation of visceral fat. Nevertheless, more epidemiological and mechanistic studies are needed to further elucidate these findings.

PAHs may influence adiposity through a few pathways. First, PAHs are structurally similar to estrogen and have estrogenic properties [37, 38], and may disrupt estrogen-mediated pathways [10, 11]. Phenanthrene and fluoranthene metabolites have been shown to have antiandrogenic effects [39]. Naphthalene metabolites may act as thyroid hormone receptor antagonists [9]. Furthermore, PAHs could alter dopamine and serotonin signaling, which may subsequently impact eating behaviors [40]. PAHs can also accumulate in tissues with high fat contents and inhibit lipolysis, thus resulting in increased fat mass accumulation and weight gain [12, 41]. It is worth noting that we observed different correlation patterns between individual PAHs and body fat distribution. A previous study suggested that PAHs may have different actions on the estrogen-responsive genes [10], which might explain the observed different impact on body fat distribution in the current study. However, this study did not include individual PAHs that were assessed in the current study, and further studies are needed to elucidate the mechanisms under varying correlations of PAHs on body fat distribution.

Interestingly, we observed racial/ethnic differences in the correlations between PAHs and body fat contents: PAHs, particularly 9-fluorene and 4-phenanthrene, had statistically significantly stronger correlations with total FM% and trunk FM% among non-Hispanic Blacks and Hispanics compared to non-Hispanic Whites. The heterogeneity may be partially explained by the racial/ethnic differences in environmental exposures noted by higher concentrations of PAHs among non-Hispanic blacks compared to non-Hispanic whites. While previous studies have found that racial/ethnic differences in total FM%, trunk FM%, or visceral fat, as well as lower metabolic rate and fat oxidation [42,43,44] in non-Hispanic Blacks compared to non-Hispanic whites, these differences have not been evaluated in the context of disparate environmental exposures that may be associated with these health outcomes. Indeed, differences in lifestyle factors may drive these racial/ethnic exposure differences. Prior studies suggested that the main exposure to 9-fluorene and 4-phenanthrene was through diet such as grilled/charred meats or foods, contaminated flour, and bread products, processed and pickled foods, and contaminated water and cow’s milk [2, 3, 45]. Moreover, second-hand smoking at home has also shown to be associated with higher levels of fluorene and phenanthrene metabolites [46]. Compared with non-Hispanic Whites, non-Hispanic Blacks tended to consume more meat prepared at higher temperatures (e.g., grilled, charred, fried) [47,48,49]. In addition, non-Hispanic Blacks and Hispanics are more likely to reside within socioeconomically disadvantaged neighborhoods and thus are disproportionally exposed to higher levels of air pollution (from traffic, factory, construction, cooking, etc.) [50] and secondhand smoke [51] compared to non-Hispanic Whites. These disproportionately higher exposure to sources of PAHs could be further modified by factors such as psychosocial stress [52] that are prevalent among non-Hispanic Black and Hispanic populations with implications on health outcomes linked to adiposity. As such, our findings of racial/ethnic differences might indicate that non-Hispanic Blacks and Hispanics are exposed to higher PAH exposures and may be particularly susceptible to PAH’s obesogenic effects due to higher exposure to PAH sources. It is also likely that there are other race/ethnic-related factors that modulate the associations between these environmental toxicants and obesity-related outcomes, although more research is needed to further elucidate potential mechanisms underlying these racial/ethnic disparities.

Our study was based in a representative population of the US adults. A major strength of our study is the evaluation of associations between PAHs and comprehensive measurements of fat depots, including fat mass and lean mass. We also controlled for a wide range of confounding factors, such as age, gender, diet intake, and serum levels of cotinine and CRP. In addition, we assessed the potential effect modification by race/ethnicity on the associations between PAHs and body fat distribution.

Nevertheless, our study had several limitations. First, the current study used the cross-sectional data, thus the temporal relations cannot be determined. In addition, we could not rule out the possibility that the increased PAH levels may reflect the increased storage and release of PAHs from adipose tissues. Second, the current study did not include measurements of some high-molecular-weight PAHs, such as the well-studied benzo[a]pyrene. High-molecular-weight PAHs were mainly excreted via bile and thus were below the detection limit in the urinary samples [45]. Third, the current study used a single urine sample to measure levels of PAHs; therefore, the long-term exposure of PAHs may not be sufficiently reflected. A prior study showed that measurements of other lipophilic chemicals in a single urine spot test is moderately predictive of long-term exposure over 3 months [53]. Future studies are warranted to replicate the findings using repeated assessments of PAHs in prospective cohort studies to offer potential insights into the temporal relationships between PAHs and body fat distribution.

Conclusion

In conclusion, various correlation patterns between urinary PAH levels and regional body fat distribution were observed in this large sample of US adult non-smokers. In comparison with non-Hispanic Whites, the positive correlations of 9-fluorene and 4-phenanthrene with total FM% and trunk FM% were statistically significantly stronger among non-Hispanic Blacks and Hispanics who had higher PAH exposures. These results suggest that these racial/ethnic groups might be particularly susceptible to the detrimental effects of PAH exposures and future studies are warranted to further explore the mechanisms underlying the racial/ethnic disparity observed in the current study. Such studies could aid in informing strategies for interventions to reduce exposure sources of PAHs in racial/ethnic disparities that contribute to obesity-related health outcomes.

Availability of data and materials

The data used in the current study is openly accessible in the NHANES website: https://www.cdc.gov/nchs/nhanes/index.htm.

References

  1. Marzooghi S, Di Toro DM. A critical review of polycyclic aromatic hydrocarbon phototoxicity models. Environ Toxicol Chem. 2017;36(5):1138–48.

    Article  CAS  Google Scholar 

  2. Polycyclic Aromatic Hydrocarbons (PAHs) Factsheet. [https://www.cdc.gov/biomonitoring/PAHs_FactSheet.html]. Accessed 20 Mar 2022.

  3. Toxicological Profile for Polycyclic Aromatic Hydrocarbons (PAHs). [https://www.atsdr.cdc.gov/toxprofiles/tp69.pdf]. Accessed 20 Mar 2022.

  4. Harris KL, Banks LD, Mantey JA, Huderson AC, Ramesh A. Bioaccessibility of polycyclic aromatic hydrocarbons: relevance to toxicity and carcinogenesis. Expert Opin Drug Metab Toxicol. 2013;9(11):1465–80.

    Article  CAS  Google Scholar 

  5. Stallings-Smith S, Mease A, Johnson TM, Arikawa AY. Exploring the association between polycyclic aromatic hydrocarbons and diabetes among adults in the United States. Environ Res. 2018;166:588–94.

    Article  CAS  Google Scholar 

  6. Ranjbar M, Rotondi MA, Ardern CI, Kuk JL. Urinary Biomarkers of Polycyclic Aromatic Hydrocarbons Are Associated with Cardiometabolic Health Risk. PLoS One. 2015;10(9):e0137536.

    Article  CAS  Google Scholar 

  7. Poursafa P, Moosazadeh M, Abedini E, Hajizadeh Y, Mansourian M, Pourzamani H, et al. A Systematic Review on the Effects of Polycyclic Aromatic Hydrocarbons on Cardiometabolic Impairment. Int J Prev Med. 2017;8:19.

    Article  Google Scholar 

  8. Vinggaard AM, Hnida C, Larsen JC. Environmental polycyclic aromatic hydrocarbons affect androgen receptor activation in vitro. Toxicology. 2000;145(2–3):173–83.

    Article  CAS  Google Scholar 

  9. Sun H, Shen OX, Xu XL, Song L, Wang XR. Carbaryl, 1-naphthol and 2-naphthol inhibit the beta-1 thyroid hormone receptor-mediated transcription in vitro. Toxicology. 2008;249(2–3):238–42.

    Article  CAS  Google Scholar 

  10. Gozgit JM, Nestor KM, Fasco MJ, Pentecost BT, Arcaro KF. Differential action of polycyclic aromatic hydrocarbons on endogenous estrogen-responsive genes and on a transfected estrogen-responsive reporter in MCF-7 cells. Toxicol Appl Pharmacol. 2004;196(1):58–67.

    Article  CAS  Google Scholar 

  11. Sievers CK, Shanle EK, Bradfield CA, Xu W. Differential action of monohydroxylated polycyclic aromatic hydrocarbons with estrogen receptors α and β. Toxicol Sci. 2013;132(2):359–67.

    Article  CAS  Google Scholar 

  12. Irigaray P, Ogier V, Jacquenet S, Notet V, Sibille P, Méjean L, 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(7):1362–72.

    Article  CAS  Google Scholar 

  13. Bushnik T, Wong SL, Holloway AC, Thomson EM. Association of urinary polycyclic aromatic hydrocarbons and obesity in children aged 3-18: Canadian Health Measures Survey 2009-2015. J Dev Orig Health Dis. 2020;11(6):623–31.

    Article  Google Scholar 

  14. 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:354–60.

    Article  CAS  Google Scholar 

  15. Scinicariello F, Buser MC. Urinary polycyclic aromatic hydrocarbons and childhood obesity: NHANES (2001-2006). Environ Health Perspect. 2014;122(3):299–303.

    Article  CAS  Google Scholar 

  16. Hendryx M, Luo J. Latent class analysis of the association between polycyclic aromatic hydrocarbon exposures and body mass index. Environ Int. 2018;121(Pt 1):227–31.

    Article  CAS  Google Scholar 

  17. Lee I, Park H, Kim MJ, Kim S, Choi S, Park J, et al. Exposure to polycyclic aromatic hydrocarbons and volatile organic compounds is associated with a risk of obesity and diabetes mellitus among Korean adults: Korean National Environmental Health Survey (KoNEHS) 2015-2017. Int J Hyg Environ Health. 2021;240:113886.

    Article  CAS  Google Scholar 

  18. Tchkonia T, Thomou T, Zhu Y, Karagiannides I, Pothoulakis C, Jensen MD, et al. Mechanisms and metabolic implications of regional differences among fat depots. Cell Metab. 2013;17(5):644–56.

    Article  CAS  Google Scholar 

  19. Vasan SK, Osmond C, Canoy D, Christodoulides C, Neville MJ, Di Gravio C, et al. Comparison of regional fat measurements by dual-energy X-ray absorptiometry and conventional anthropometry and their association with markers of diabetes and cardiovascular disease risk. Int J Obes. 2018;42(4):850–7.

    Article  CAS  Google Scholar 

  20. Chen GC, Arthur R, Iyengar NM, Kamensky V, Xue X, Wassertheil-Smoller S, et al. Association between regional body fat and cardiovascular disease risk among postmenopausal women with normal body mass index. Eur Heart J. 2019;40(34):2849–55.

    Article  CAS  Google Scholar 

  21. Stults-Kolehmainen MA, Stanforth PR, Bartholomew JB, Lu T, Abolt CJ, Sinha R. DXA estimates of fat in abdominal, trunk and hip regions varies by ethnicity in men. Nutr Diabetes. 2013;3(3):e64.

    Article  CAS  Google Scholar 

  22. Hudson-Hanley B, Smit E, Branscum A, Hystad P, Kile ML. Trends in urinary metabolites of polycyclic aromatic hydrocarbons (PAHs) in the non-smoking U.S. population, NHANES 2001–2014. Chemosphere. 2021;276:130211.

    Article  CAS  Google Scholar 

  23. National Health and Nutrition Examination Survey: Plan and Operations, 1999–2010. [https://www.cdc.gov/nchs/data/series/sr_01/sr01_056.pdf]. Accessed 20 Mar 2022.

  24. Chiolero A, Faeh D, Paccaud F, Cornuz J. Consequences of smoking for body weight, body fat distribution, and insulin resistance. Am J Clin Nutr. 2008;87(4):801–9.

    Article  CAS  Google Scholar 

  25. Romanoff LC, Li Z, Young KJ, Blakely NC 3rd, Patterson DG Jr, Sandau CD. Automated solid-phase extraction method for measuring urinary polycyclic aromatic hydrocarbon metabolites in human biomonitoring using isotope-dilution gas chromatography high-resolution mass spectrometry. J Chromatogr B Anal Technol Biomed Life Sci. 2006;835(1–2):47–54.

    Article  CAS  Google Scholar 

  26. Li Z, Romanoff LC, Lewin MD, Porter EN, Trinidad DA, Needham LL, et al. Variability of urinary concentrations of polycyclic aromatic hydrocarbon metabolite in general population and comparison of spot, first-morning, and 24-h void sampling. J Expo Sci Environ Epidemiol. 2010;20(6):526–35.

    Article  CAS  Google Scholar 

  27. National Health and Nutrition Examination Survey: Technical Documentation for the 1999-2004. Dual Energy X-Ray Absorptiometry (DXA) Multiple Imputation Data files. [https://wwwn.cdc.gov/Nchs/Data/Nhanes/Dxa/dxa_techdoc.pdf]. Accessed 20 Mar 2022.

  28. Lohman TG, Roche AF, Martorell R. Anthropometric standardization reference manual. Champaign: Human Kinetics Books; 1988.

    Google Scholar 

  29. National Health and Nutrition Examination Survey Questionnaires, Datasets, and Related Documentation. [https://wwwn.cdc.gov/nchs/nhanes/]. Accessed 20 Mar 2022.

  30. Cotinine Factsheet. [https://www.cdc.gov/biomonitoring/Cotinine_FactSheet.html]. Accessed 20 Mar 2022.

  31. Rosner B. Fundamentals of biostatistics. Boston: Brooks/Cole, Cengage Learning; 2011.

    Google Scholar 

  32. Zong G, Grandjean P, Wu H, Sun Q. Circulating persistent organic pollutants and body fat distribution: Evidence from NHANES 1999-2004. Obesity (Silver Spring). 2015;23(9):1903–10.

    Article  CAS  Google Scholar 

  33. Rubin DB, Schenker N. Multiple Imputation for Interval Estimation From Simple Random Samples With Ignorable Nonresponse. J Am Stat Assoc. 1986;81(394):366–74.

    Article  Google Scholar 

  34. Romero-Corral A, Lopez-Jimenez F, Sierra-Johnson J, Somers VK. Differentiating between body fat and lean mass-how should we measure obesity? Nat Clin Pract Endocrinol Metab. 2008;4(6):322–3.

    Article  Google Scholar 

  35. Santosa S, Jensen MD. Why are we shaped differently, and why does it matter? Am J Phys Endocrinol Metab. 2008;295(3):E531–5.

    Article  CAS  Google Scholar 

  36. Wronska A, Kmiec Z. Structural and biochemical characteristics of various white adipose tissue depots. Acta physiologica (Oxford, England). 2012;205(2):194–208.

    Article  CAS  Google Scholar 

  37. Schultz TW, Sinks GD. Xenoestrogenic gene expression: structural features of active polycyclic aromatic hydrocarbons. Environ Toxicol Chem. 2002;21(4):783–6.

    Article  CAS  Google Scholar 

  38. Wenger D, Gerecke AC, Heeb NV, Schmid P, Hueglin C, Naegeli H, et al. In vitro estrogenicity of ambient particulate matter: contribution of hydroxylated polycyclic aromatic hydrocarbons. J Appl Toxicol. 2009;29(3):223–32.

    Article  CAS  Google Scholar 

  39. Chang CS, Liao SS. Topographic recognition of cyclic hydrocarbons and related compounds by receptors for androgens, estrogens, and glucocorticoids. J Steroid Biochem. 1987;27(1–3):123–31.

    Article  CAS  Google Scholar 

  40. Stephanou P, Konstandi M, Pappas P, Marselos M. Alterations in central monoaminergic neurotrasmission induced by polycyclic aromatic hydrocarbons in rats. Eur J Drug Metab Pharmacokinet. 1998;23(4):475–81.

    Article  CAS  Google Scholar 

  41. 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(1):35–9.

    Article  CAS  Google Scholar 

  42. Carroll JF, Chiapa AL, Rodriquez M, Phelps DR, Cardarelli KM, Vishwanatha JK, et al. Visceral Fat, Waist Circumference, and BMI: Impact of Race/ethnicity. Obesity. 2008;16(3):600–7.

    Article  Google Scholar 

  43. Rahman M, Temple JR, Breitkopf CR, Berenson AB. Racial differences in body fat distribution among reproductive-aged women. Metab Clin Exp. 2009;58(9):1329–37.

    Article  CAS  Google Scholar 

  44. Weyer C, Snitker S, Bogardus C, Ravussin E. Energy metabolism in African Americans: potential risk factors for obesity. Am J Clin Nutr. 1999;70(1):13–20.

    Article  CAS  Google Scholar 

  45. Li Z, Sandau CD, Romanoff LC, Caudill SP, Sjodin A, Needham LL, et al. Concentration and profile of 22 urinary polycyclic aromatic hydrocarbon metabolites in the US population. Environ Res. 2008;107(3):320–31.

    Article  CAS  Google Scholar 

  46. Jain RB. Trends and concentrations of selected polycyclic aromatic hydrocarbons in general US population: Data from NHANES 2003–2008. Cogent Environ Sci. 2015;1(1):1031508.

    Article  Google Scholar 

  47. Wang Y, Beydoun MA, Caballero B, Gary TL, Lawrence R. Trends and correlates in meat consumption patterns in the US adult population. Public Health Nutr. 2010;13(9):1333–45.

    Article  Google Scholar 

  48. García-Jiménez CI, Mishra AK. Role of ethnicity in consumption of meat products. Appl Econ Lett. 2011;18(7):665–9.

    Article  Google Scholar 

  49. Bogen KT, Keating GA. U.S. dietary exposures to heterocyclic amines*. J Expo Sci Environ Epidemiol. 2001;11(3):155–68.

    Article  CAS  Google Scholar 

  50. Tessum CW, Apte JS, Goodkind AL, Muller NZ, Mullins KA, Paolella DA, et al. Inequity in consumption of goods and services adds to racial–ethnic disparities in air pollution exposure. Proc Natl Acad Sci. 2019;116(13):6001–6.

    Article  CAS  Google Scholar 

  51. Shastri SS, Talluri R, Shete S. Disparities in Secondhand Smoke Exposure in the United States: National Health and Nutrition Examination Survey 2011-2018. JAMA Intern Med. 2021;181(1):134–7.

    Article  Google Scholar 

  52. Williams DR. Race, Stress, and Mental Health: Findings from the Commonwealth Minority Health Survey. In: Hogue C, Hargraves MA, Collins KS, editors. Minority Health in America: Findings and Policy Implication From the Commonwealth Fund Minority Health Survey. Baltimore: Johns-Hopkins University Press; 2000. p. 209–43.

    Google Scholar 

  53. Hauser R, Meeker JD, Park S, Silva MJ, Calafat AM. Temporal variability of urinary phthalate metabolite levels in men of reproductive age. Environ Health Perspect. 2004;112(17):1734–40.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We acknowledge and thank members of the National Center for Health Statistics (NCHS) of the Centers for Disease Control (CDC) and Prevention and the participants who contributed time and data to the National Health and Nutrition Examination Survey.

Funding

This research was supported by NIH grant ES022981. TJT was supported by NIEHS grant P30ES000002. YW is supported by the NIH T32 HL 098048 postdoctoral fellowship.

Author information

Authors and Affiliations

  1. Department of Nutrition, Harvard T.H. Chan School of Public Health, 677 Huntington Avenue, Boston, MA, 02115, USA

    Yeli Wang, Lu Zhu & Qi Sun

  2. Department of Environmental Health, Harvard T.H. Chan School of Public Health, 655 Huntington Avenue, Boston, MA, 02115, USA

    Tamarra James-Todd

  3. Department of Epidemiology, Harvard T.H. Chan School of Public Health, 677 Huntington Avenue, Boston, MA, 02115, USA

    Tamarra James-Todd & Qi Sun

  4. Channing Division of Network Medicine, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, 181 Longwood Avenue, Boston, MA, 02115, USA

    Qi Sun

Authors
  1. Yeli Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Lu Zhu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Tamarra James-Todd
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Qi Sun
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

YW wrote the manuscript draft. YW and LZ conducted the statistical analysis. All authors reviewed the manuscript critically and approved the content of the manuscript.

Corresponding author

Correspondence to Qi Sun.

Ethics declarations

Ethics approval and consent to participate

Written informed consents were obtained from all participants. The study protocol was approved by the institutional review board at the Centers for Disease Control and Prevention (Atlanta, Georgia).

Consent for publication

All authors gave consent for the paper to be published in Environmental Health.

Competing interests

None.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplemental Table S1.

Pearson correlation among PAHs among non-smokers. Supplemental Figure S1. Flowchart of the study design of the current analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, Y., Zhu, L., James-Todd, T. et al. Urinary polycyclic aromatic hydrocarbon excretion and regional body fat distribution: evidence from the U.S. National Health and Nutrition Examination Survey 2001–2016. Environ Health 21, 75 (2022). https://doi.org/10.1186/s12940-022-00890-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12940-022-00890-8

Keywords

  • Polycyclic aromatic hydrocarbons
  • non-smoker
  • racial/ethnic difference
  • obesity
  • body fat distribution
  • NHANES

Environmental Health

ISSN: 1476-069X