Urinary Volatile Organic Compound Metabolites are Associated with High Blood Pressure Among Non-smoking Participants in the National Health and Nutrition Examination Survey (2011-2018) ======================================================================================================================================================================================== * Katlyn E. McGraw * Arce Domingo-Relloso * Daniel W. Riggs * Danielle N. Medgyesi * Raghavee Neupane * Jeanette A. Stingone * Tiffany R. Sanchez ## ABSTRACT **Background** Volatile organic compounds (VOCs) are ubiquitous environmental pollutants. Exposure to VOCs is associated with cardiovascular disease (CVD) risk factors, including elevated blood pressure (BP) in susceptible populations. However, research in the general population, particularly among non-smoking adults, is limited. We hypothesized that higher VOC exposure is associated with higher BP and hypertension, among non-smokers. **Methods** We included four cycles of data (2011-2018) of non-smoking adults (n=4,430) from the National Health and Nutrition Examination Survey (NHANES). Urinary VOC metabolites were measured by ultra-performance liquid chromatography–mass spectrometry, adjusted for urine dilution, and log-transformed. We estimated mean differences in BP using linear models and prevalence ratio of stage 2 hypertension using modified Poisson models with robust standard errors. Models were adjusted for age, sex, race and ethnicity, education, body mass index, estimated glomerular filtration rate and NHANES cycle. **Results** Participants were 54% female, with a median age of 48 years, 32.3% had hypertension, and 7.9% had diabetes. The mean differences (95% CI) in systolic BP were 1.61 (0.07, 3.15) and 2.46 (1.01, 3.92) mmHg when comparing the highest to lowest quartile of urinary acrolein (CEMA) and 1,3-butadiene (DHBMA) metabolites. The prevalence ratios (PR) for hypertension were 1.06 (1.02, 1.09) and 1.05 (1.01, 1.09) when comparing the highest to lowest quartiles of urinary acrolein (CEMA) and 1,3-butadiene (DHBMA), respectively. **Conclusions** Exposure to VOCs may be relevant yet understudied environmental contributors to CVD risk in the non-smoking, US population. Keywords * Volatile organic compounds * blood pressure * mixtures * acrolein * butadiene ## INTRODUCTION Volatile organic compounds (VOCs) are an under regulated class of chemicals characterized by low molecular weight and high vapor pressure, properties which lead to rapid vaporization in ambient air.1–4 Ambient VOCs are critical in the formation of harmful secondary air pollutants such as ozone and PM2.5.5, 6 However, VOCs emitted directly from sources are just as harmful and include tobacco smoke, vehicle exhaust, industrial solvents, manufacture of consumer goods, and oil and gas refining.7, 8 Indoor, VOC levels can be up to ten times higher than ambient VOC levels.8 Additionally, some VOCs are produced endogenously through processes like lipid peroxidation, glycation, and amino acid oxidation.9–11 Therefore, VOCs are ubiquitous environmental exposures. Urinary biomarkers of VOC exposure provide proximal, personal-level assessments of all routes of exposure to VOCs, including ambient, indoor, and endogenous sources. VOC metabolites from human biomonitoring have been associated with physician-diagnosed cardiovascular disease (CVD),12 metabolic syndrome,13 CVD risk factors,14–17 and ischemic heart disease mortality,18 particularly for acrolein, benzene, 1,3-butadiene, crotonaldehyde, and styrene. However, these studies identified associations in small, geographically homogenous cohorts, and in CVD-susceptible populations, including individuals who smoke. The relation between VOC biomarkers and blood pressure among non-smoking adults has not yet been elucidated in a large, representative sample of the U.S. population, such as the National Health and Nutrition Examination Survey (NHANES). In this study, we evaluated the association between participants’ VOC metabolite levels in urine with measured blood pressure (BP) levels and hypertension defined as systolic BP ≥140 mmHg, diastolic BP ≥90 mmHg, or self-reported use of BP medications19 using data from a representative sample of the non-smoking U.S. adult population participating in NHANES from 2011-2018. We hypothesized that exposure to VOCs, specifically urinary metabolites of acrolein, benzene, 1,3-butadiene, crotonaldehyde, and styrene, would be associated with higher BP levels and prevalence of hypertension. Additionally, we conducted an exploratory assessment of all urinary VOC metabolites measured in NHANES and mixtures of urinary VOC metabolites. ## METHODS ### Study Population Led by the National Center for Health Statistics (NCHS) at the Centers for Disease Control and Prevention, NHANES is a series of continuous surveys designed to assess the health and nutritional status of adults and children in the U.S. NHANES is designed as a multiyear, stratified, clustered four-stage sample of noninstitutionalized civilians with fixed sample-size targets for sampling domains defined by age, sex, race and ethnicity, and socioeconomic status, with data released in 2-year cycles. Participants gave informed consent of the survey process and their rights as a participant, and the survey was approved by the NCHS Review Board.20 Questionnaires were administered in-home followed by standardized health examinations in specially equipped mobile examination centers. Deidentified and detailed health datasets are publicly available on the NHANES website.21 We combined NHANES data from four cycles (2011-2018) to create a larger and more geographically diverse sample, including all cycles with available data on BP and urinary VOC metabolites. ### Exclusion Criteria Of the 39,156 NHANES participants from 2011-2018 combined, 23,825 were 18 years of age and older. We focused on adults because blood pressure varies in childhood and adolescence. There were 7,626 participants with measured urinary VOC metabolites and BP examination data because 1/3 subsample at each cycle are analyzed for environmental chemicals. Because tobacco smoke is a primary source of VOC exposure, self-reported cigarette smoking, and serum cotinine data were used to determine smoking status. We excluded 1,405 individuals who reported currently smoking cigarettes, an additional 451 individuals with serum cotinine levels >10 µg/L, reflecting recent nicotine exposure, and 480 missing self-reported cigarette smoking status or cotinine data, for a final sample of 5,290 non-smoking participants. We further excluded an additional 860 participants with missing data, including 492 missing priority urinary VOCs, 172 missing systolic blood pressure (SBP), 11 missing diastolic blood pressure (DBP), 81 missing urine creatinine, 42 missing BMI, 35 missing diabetes status, 21 missing estimated glomerular filtration rate (eGFR), 3 missing triglycerides (TG), 2 missing high density lipoproteins (HDL), and 1 missing education, for a total of 4,430 participants included in our analysis (Figure S1). ### Urinary VOC Metabolites Urine specimens were collected at mobile examination centers, frozen at -20°C, and shipped to the Division of Environmental Health Laboratory Sciences, National Center for Environmental Health for analysis.22 Ultra-performance liquid chromatography coupled with electrospray tandem mass spectrometry was used to measure 35 VOC metabolites of 20 parent compounds in urine in nanograms per milliliter (ng/mL).23 We report results for 5 hypothesized, priority parent VOCs and their 9 associated urinary metabolites: acrolein (CEMA – N-Acetyl-S-(2-carboxyethyl)-L-cysteine and 3HPMA – N-Acetyl-S-(3-hydroxypropyl)-L-cysteine), benzene (MU – t,t-Muconic Acid and PMA – N-Acetyl-S-(phenyl)-L-cysteine), 1,3-butadiene (DHBMA – N-Acetyl-S-(3,4-dihydroxybutyl)-L-cysteine and MHBMA3 – N-Acetyl-S-(4-hydroxy-2-butenyl)-L-cysteine), crotonaldehyde (HPMMA – N-Acetyl-S-(3-hydroxypropyl-1-methyl)-L-cysteine), and styrene (PGA – N-Acetyl-S-(phenyl)-L-cysteine and MA – Mandelic acid). PGA and MA are metabolites of ethylbenzene and styrene. See Table S1 for all parent compounds, long and short metabolite names, NHANES cycle availability, and limits of detection (LOD). Values below the LOD were imputed with the LOD/√2. Although the benzene metabolite, MU, was only available for one cycle (2017-2018, n=1,046), we included MU in primary analyses for comparison to the other benzene metabolite available at the other three cycles, PMA (2011-2016, n=3,207). To correct for urine dilution, we divided all urinary VOC metabolites by urine creatinine and reported levels as nanogram of VOC metabolite per milligrams of creatinine (ng/mg). ### Blood Pressure and Hypertension BP measurements were collected at mobile examination centers using a sphygmomanometer and appropriately sized cuff for the participant’s arm circumference. Three BP readings were collected after five minutes of seated rest. A fourth measurement was collected if a measurement was interrupted or incomplete. All recorded measurements of SBP and DBP were averaged and used as continuous outcome variables. For hypertension, our primary endpoint was defined by having either a mean SBP ≥140 mmHg, a mean DBP ≥90 mmHg, or self-reported use of BP medications , also known as stage 2 hypertension.19 We also conducted sensitivity analyses using two alternate definitions: 1) stage 1 hypertension, where participants were hypertensive if mean SBP ≥130 mmHg, mean DBP ≥80 mmHg, or were taking BP medications, and 2) if individuals self-reported physician-diagnosed hypertension. ### Covariates Age, sex, race and ethnicity, education, and household income were acquired from self-reported questionnaires. Race and ethnicity were classified as Mexican American, Other Hispanic, Non-Hispanic White, Non-Hispanic Black, Non-Hispanic Asian, and Other Races, including multiple races. We reclassified education as less than a high school education, a high school education or GED equivalent, some college, or bachelor’s degree or more. Household income categories were reduced to $0-$24,999, $25,000-$54,999, $55,000-$74,999, and ≥$75,000. Body mass index (BMI, kg/m2) was collected during examination at mobile examination center. Estimated glomerular filtration rate (eGFR) was calculated using the Chronic Kidney Disease Epidemiology Collaboration equation, 24 which estimates GFR for age, sex, and serum creatinine in µmol/L. The eGFR can influence metabolite excretion in urine and was therefore used for adjustment in our models. Total cholesterol and high-density lipoprotein (HDL) cholesterol were measured in blood collected at mobile examination centers. Diabetes status was determined based on self-reported responses that the participant was ever told by a doctor they had diabetes or were taking medications to lower blood sugar. Prevalent CVD status was determined based on self-reported responses that the participant was ever told by a doctor they had congestive heart failure, coronary heart disease, angina or angina pectoris, heart attack, or stroke. ### Statistics Combined NHANES survey cycles and weights produce estimates representative of the U.S. civilian noninstitutionalized population at the midpoint of the combined survey period. We constructed new sample weights for the combined cycles (4 cycles, 2011-2018) by multiplying 2-year subsample weights for environmental chemicals by 1/4 (for 4 NHANES cycles) as described.25 To characterize our study population, we stratified by hypertension status and continuous and categorical variables of interest using t-tests and χ² tests. We used linear regression models for normally distributed BP and modified Poisson regression with robust standard errors for dichotomous hypertension outcomes to estimate prevalence ratios of hypertension.26 Urine VOC metabolite levels were modeled as: (1) quartiles (to compare each of the highest three quartiles to the lowest quartile), (2) per interquartile range (IQR) on log-transformed levels (to compare the 75th to the 25th percentile), and (3) log-transformed levels with restricted quadratic splines (to evaluate the flexible dose-response relationship). Because NHANES represents a younger, more general population, we chose to adjust for sociodemographic factors age, sex, race and ethnicity, education, eGFR, and NHANES cycle year and BMI in our main model, Model 1, as to not adjust away potential mediating CVD risk factors based on literature reviews for BP and hypertension.27, 28 In Model 2, we further adjust for traditional CVD risk factors including total cholesterol, HDL, triglycerides, diabetes status, and BP medications. We present 9 priority urinary VOC metabolites of parent VOCs acrolein, benzene, 1,3-butadiene, crotonaldehyde, and styrene, from Model 1, in the main text. Model estimates are population weighted mean differences in BP or population weighted prevalence ratio of hypertension with 95% confidence intervals (CI). All central tendency estimates, proportions and effect estimates are population weighted. ### Secondary analyses In secondary analyses, we examined the other 12 VOC metabolites measured in NHANES. We used Wald tests and conducted subgroup analysis to assess effect modification by subgroups of age, sex, and race and ethnicity. To address potential co-pollutant exposure and collinearity, we used hierarchical Bayesian Kernel Machine Regression (BKMR). Because BKMR cannot account for complex survey weighting, we consider these exploratory analyses. BKMR is a kernel-regression based machine learning method that characterizes the exposure response function of multiple predictors on a health outcome while other predictors are fixed to a specific percentile.29, 30 The method also allows us to examine statistical interactions between VOC metabolites within the mixture and joint associations between the whole mixture and health outcomes. The hierarchical BKMR uses hierarchical variable selection to create posterior inclusion probabilities (PIPs) that quantify the relative importance by selecting 1) at the group level and 2) by selecting exposures of each group.29 All metabolites with the same parent compound were grouped in hierarchical analysis except for MHBMA3 and HPMMA which were more highly correlated than DHBMA and MHBMA3 (Figure S2). We excluded the benzene metabolites, MU and PMA, due to missing cycles. BKMRs were run with 10,000 iterations and the tuning parameter r.jump value of 0.1 for continuous SBP and DBP and dichotomous hypertension to achieve higher acceptance rates. ### Sensitivity Analyses We conducted several sensitivity analyses. Because a number of participants reported use of BP lowering medications, we assessed differences between measured and underlying BP.31 For individuals self-reporting use of BP lowering medications we imputed underlying BP values by adding a 10 mm Hg or 5 mm Hg constant to measured SBP or DBP, respectively.31 Additionally, we assessed differences between Stage 1, Stage 2, and physician-diagnosed hypertension. Data analysis was performed in R (version 3.1.3)32 using the nhanesA33, tidyverse34 and survey package35 to account for the complex survey design and sampling weights. ## RESULTS The characteristics of NHANES participants by hypertension status are shown in Table 1. Participants were 54% female, with a median age of 48 years, 32.3% had hypertension, and 17.9% had diabetes mellitus. Participants with hypertension were older and more likely to be non-Hispanic White or non-Hispanic Black, had higher BMI, lower eGFR, lower HDL cholesterol, higher triglycerides, higher total cholesterol, and were more likely to have self-reported physician diagnosed CVD and diabetes mellitus. Additionally, participants with hypertension had higher median levels of all priority urinary VOC metabolites compared to participants without hypertension (Table 2 and Table S2), except for the PMA benzene metabolite. View this table: [Table 1.](http://medrxiv.org/content/early/2024/07/19/2024.07.18.24310671/T1) Table 1. Participant characteristics across categories of hypertension, restricted to non-smoking participants with available blood pressure examination data and urinary VOCs, n=4,430 NHANES 2011-2018 participants. Continuous variables are reported as median [interquartile range, IQR] and categorical variables as number (%). Percentages are population weighted. View this table: [Table 2.](http://medrxiv.org/content/early/2024/07/19/2024.07.18.24310671/T2) Table 2. Median and interquartile ranges of VOC metabolite levels (ng/mg creatinine) across categories of hypertension, restricted to participants with available blood pressure examination data and urinary VOCs, among 4,430 NHANES 2011-2018 participants. MU was measured in n=1,046 NHANES 2017-2018 participants and PMA was measured in n=3,207 NHANES 2011-2016 participants. Results for priority VOC metabolites and BP outcomes are shown in Table 3. For acrolein, the urinary metabolite CEMA was associated with higher SBP and a higher prevalence of hypertension. When comparing the highest to lowest quartile, the mean difference (MD) (95%CI) was 1.61 (0.07, 3.15) mmHg for SBP, and the prevalence ratio (PR) (95%CI) of hypertension was 1.06 (1.02, 1.09). Estimates for SBP and hypertension were also statistically significant per log IQR-change and, although attenuated, generally remained significant after further adjustment for CVD risk factors. In the restricted quadratic spline models, there were clear positive dose-response relationships with SBP and PR of hypertension observed for CEMA above 34.20 ng/mg creatinine (Figures 1 and 2). For the other acrolein metabolite, 3HPMA, results were positive but not significant for SBP (MD [95%CI]: 1.28 [-0.33, 2.90] mmHg); however, the association with hypertension was statistically significant (PR [95%CI]: 1.05 [1.01, 1.09]). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/19/2024.07.18.24310671/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2024/07/19/2024.07.18.24310671/F1) Figure 1. Mean differences (MD) (95% confidence interval) of systolic blood pressure (SBP) by urinary VOC levels (ng/mg creatinine) modeled as restricted quadratic splines for priority VOC metabolites among nonsmoking NHANES 2011-2018 participants (n=4,430). Lines (shaded areas) represent the MD (95%CI) of SBP and RR (95% CI) of hypertension by VOC metabolites modeled as restricted quadratic splines for log transformed VOC metabolite distributions with knots at 10th, 50th, and 90th percentiles. The reference value was set at the 10th percentile. Models were adjusted for age, sex, race and ethnicity, education, BMI, eGFR, and NHANES cycle year. The histograms in the background represent the distribution of each VOC metabolite (ng/mg creatinine). See Figure S2 for DBP. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/19/2024.07.18.24310671/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2024/07/19/2024.07.18.24310671/F2) Figure 2. Prevalence ratios (PR) (95% confidence interval, CI) of hypertension by urinary VOC levels (ng/mg creatinine) modeled as restricted quadratic splines among nonsmoking NHANES 2011-2018 participants (n=4,430). Lines (shaded areas) represent the PR (95%CI) of hypertension by VOC metabolites modeled as restricted quadratic splines for log transformed VOC metabolite distributions with knots at 10th, 50th, and 90th percentiles. The reference value was set at the 10th percentile. Models were adjusted for age, sex, race and ethnicity, education, BMI, eGFR, and NHANES cycle year. The histograms in the background represent the distribution of each VOC metabolite (ng/mg creatinine). See Figure S2 for DBP. View this table: [Table 3.](http://medrxiv.org/content/early/2024/07/19/2024.07.18.24310671/T3) Table 3. Mean differences (95% confidence interval) of systolic blood pressure (SBP), diastolic blood pressure (DBP), and prevalence ratios – PR (95% confidence interval) of hypertension by levels (ng/mg creatinine) of urinary VOC metabolites, restricted to participants with available blood pressure examination data and urinary VOCs, n=4,430 NHANES 2011-2018 participants. MU was measured in n=1,046 NHANES 2017-2018 participants and PMA was measured in n=3,207 NHANES 2011-2016 participants. Model 1 was adjusted for age, sex, race and ethnicity, education, BMI, eGFR, and NHANES cycle year. Model 2 was additionally adjusted for HDL-cholesterol, triglycerides, total cholesterol, diabetes status, and antihypertensive medication. Likewise, for 1,3-butadiene, the urinary metabolite DHBMA was associated with higher SBP and a higher prevalence of hypertension. When comparing the highest to lowest quartile of DHBMA, the MD in SBP was 2.46 (1.01, 3.92) mmHg and the PR of hypertension was 1.05 (1.01, 1.09). Estimates for SBP and hypertension were also statistically significant per log IQR-change and, although attenuated, generally remained significant after further adjustment for CVD risk factors. The linear association apparent in the quartile models was also observed with the restricted quadratic spline models, with clear positive dose-response relationships with SBP and PR of hypertension observed for DHBMA above 252 ng/mg creatinine (Figures 1 and 2). For the other 1,3-butadiene metabolite, MHBMA3, results were less consistent, but there was an association with hypertension per log IQR of MHBMA3 (PR [95%CI]: 1.01 [1.00, 1.02]). For the other priority metabolites, results were less consistent. The benzene metabolite, MU, was associated with higher SBP per log IQR of MU (MD [95%CI]: 0.76 [0.05, 1.48] mmHg), but not with higher prevalence of hypertension. When comparing the highest quartile to the lowest quartile, the crotonaldehyde metabolite, HPMMA, and the styrene metabolite, MA, were associated with higher prevalence of hypertension (PR [95%CI]: 1.06 [1.02, 1.10], and 1.05 [1.01, 1.08]), respectively. Estimates for HPMMA and MA and SBP trended positive but were not statistically significant. Additionally, we found no significant results for any of the priority VOC metabolites and DBP (Table 3 and Figure S3). Results for the other urinary VOC metabolites available in NHANES (i.e., parent compounds: acrylamide, acrylonitrile, 1-bromopropane, carbon disulfide, cyanide, n,n-dimethylformamide, isoprene, propylene oxide, toluene, and xylene) are reported in Table S3. In Wald tests for interaction, we found significant interactions between age and all priority VOC metabolites except for the benzene metabolite PMA (Table S4). In models stratified by age groups for priority VOCs, the associations with SBP were stronger for the age subgroup 40-64 for metabolites of benzene (MU), 1,3-butadiene (DHBMA), and styrene (PGA). For hypertension, associations were stronger for the age group 40-64 and ≥65 for metabolites of acrolein (CEMA) and crotonaldehyde (HPMMA), the age group 40-64 for the metabolites of benzene (MU) and styrene (PGA), and the age group ≥65 for metabolites of 1,3-butadiene (MHBMA3) and styrene (MA). We found significant interactions between sex and metabolites of acrolein, 1,3-butadiene, and crotonaldehyde. By sex, the association with SBP was stronger in females for metabolites of acrolein (CEMA, 3HPMA), 1,3-butadiene (DHBMA) and crotonaldehyde (HPMMA). We found significant interactions between race and ethnicity for all VOC metabolites except for the benzene metabolite PMA and the styrene metabolite PGA. The association with SBP was stronger for Mexican Americans for 1,3-butadiene (DHBMA), for Non-Hispanic Asians for acrolein (CEMA), and for Non-Hispanic Blacks for styrene (MA). The association with hypertension was stronger for Non-Hispanic Whites for metabolites of acrolein (3HPMA), benzene (MU), crotonaldehyde (HPMMA) and styrene (MA). Patterns for other VOCs were inconsistent. In sensitivity analyses, we found that effect sizes of associations with underlying BP were consistent with the main analyses of original, measured BP (data not shown). Similarly, the use of Stage 2 or physician diagnosed definitions of hypertension did not significantly change effect sizes; however, use of stage 1 hypertension outcome attenuated associations, except for the styrene metabolite (MA) (Figure S4). Urinary VOC metabolites originating from the same parent compound were highly correlated except for 1,3-butadiene metabolites, where MHBMA3 was more highly correlated with crotonaldehyde metabolite, HPMMA (Figure S2). We found no significant association between the VOC metabolite mixture and SBP (Figure S5). However, the VOC mixture and overall risk of hypertension was significant and positive. Conditional posterior inclusion probabilities were higher for acrolein metabolite, CEMA, and 1,3-butadiene metabolite, DHBMA, for SBP and CEMA for hypertension (Table S5). ## DISCUSSION In this nationally representative study of non-smoking adults, we found that the acrolein metabolite, CEMA, and the 1,3-butadiene metabolite, DHBMA, were consistently and significantly associated with higher SBP and higher prevalence of hypertension in quartile, log IQR, and restricted spline models. The other priority VOC metabolites were less consistent across BP outcomes. Estimates were attenuated after adjustment for cardiovascular risk factors. Exploratory hierarchical BKMR analyses supported acrolein and 1,3-butadiene as important contributors to SBP and hypertension risk and the overall mixture of VOC metabolites (excluding benzene) was associated with hypertension. Elevated BP is an important risk factor for the development of CVD and is associated with the highest population attributable fraction for CVD deaths in the U.S. (40.6%).36 Systolic hypertension is usually the result of arterial thickening and stiffening caused either by atherosclerosis, medial degradation, or impaired endothelial-mediated vasodilation,37 processes that have been shown to be sensitive to VOCs in experimental studies.14, 38 Taken together, these observations suggest that exposure to VOCs, particularly acrolein and 1,3-butadiene, may be significant contributors to the development of CVD, in the general non-smoking population. Among NHANES 2011-2018 participants the acrolein metabolite, CEMA, was associated with higher SBP and prevalence of hypertension and the 3HPMA metabolite was associated with higher prevalence of hypertension. Acrolein is a pesticide, combustion byproduct, precursor in the production of other chemicals, and byproduct of lipid peroxidation, glycation, and amino acid oxidation; 3HPMA and CEMA account for about 30% of the dose from exposure.39 Epidemiological evidence shows acrolein metabolites are associated with SBP among nonsmoking adults, including in a cohort of 308 adults in Louisville, KY of mean age 51.7 years (MD [95%CI] CEMA: 1.2 [-0.3, 3.3] mmHg and 3HPMA: 0.98 [0.04, 1.9] mmHg)14 and a cohort of 778 Black participants in Jackson, MS of mean age 51.3 years (CEMA: 1.6 [0.4, 2.7] mmHg and 3HPMA: 0.8 [0.01, 1.6] mmHg).40 However, associations between acrolein metabolites and hypertension were null among participants in Jackson, MS. A 2022 analysis of ambient VOC exposure in the Sister Study (n=47,467) found that ambient acrolein exposure was associated with a greater risk of hypertension, particularly among never smoking women 35-74 years of age.41 Likewise, in a nested case-cohort study of non-smoking individuals in northeastern Iran (n=1,198), there was an increased hazard ratio of ischemic heart disease mortality associated with metabolites of acrolein (CEMA, 2.11 (1.31, 3.40) and 3HPMA, 2.99 (1.60, 5.59)) when comparing the highest tertile to the lowest.18 In this study, we also found the 1,3-butadiene metabolite, DHBMA, was positively associated with higher SBP and greater risk of hypertension, and MHBMA3 was associated with higher prevalence of hypertension. Most human exposure from 1,3-butadiene occurs from incomplete combustion of fossil fuels and burning of biomass.42 DHBMA is the most abundant metabolite in urine due to its natural background, which may be caused by endogenous sources.43, 44 It is unknown whether high background levels of DHBMA are due to 1,3-butadiene exposures or endogenous sources like lipid peroxidation. Additionally, crotonaldehyde is produced endogenously from the metabolism of 1,3-butadiene and may contribute to higher levels of DHBMA.45–47 Although MHBMA3 is a sensitive marker of 1,3-butadiene exposure, DHBMA is not a sensitive marker,48 and does not decrease with smoking cessation, compared to MHBMA3.49 Therefore, DHBMA may be a relevant biomarker for CVD risk,50 but may not be indicative of 1,3-butadiene exposure. In the Sister Study (n=47,467), ambient 1,3-butadiene at the census tract level was associated with 1.05 (95% CI 1.00, 1.10) greater risk of hypertension.41 Additionally, there was a higher ischemic heart disease mortality hazard ratio associated with 1,3-butadiene [DHBMA, 2.49 (1.50, 4.15)] in a recent study of non-smoking northeastern Iranians (n=1,198).18 However, associations were null among participants in the Louisville, KY and Jackson, MS studies.14, 40 More studies of low-level 1,3-butadiene exposures are required to better assess the sensitivity of DHBMA as an exposure biomarker. The other priority metabolites were associated with either higher SBP or higher prevalence of hypertension, but not both. The benzene metabolite, MU (but not PMA), was associated with higher SBP. This may be due to smaller sample size compared with the other VOC metabolites. Limited, consistent evidence exists of urinary benzene metabolites and ambient benzene and cardiovascular outcomes due to ongoing analytical issues.51, 52 In the Sister Study, ambient benzene was significantly associated with a 1.06 (95% CI 1.01, 1.10) greater risk of hypertension among never smokers.41 The crotonaldehyde metabolite, HPMMA, was associated with higher prevalence of hypertension. Among nonsmoking Black participants in Jackson, MS, HPMMA was associated with higher SBP (MD: 1.2 [95%CI: 0.08, 1.8] mm Hg),40 but associations with hypertension were null. Similarly, results for HPMMA and BP were null in the Louisville, KY cohort.14 The styrene metabolite, MA (but not PGA), was associated with higher prevalence of hypertension. Consistent with our results, the Consortium of Safe Labor showed that ambient styrene was associated with greater odds of hypertension and higher SBP and DBP among normotensive women.53, 54 Additionally, two longitudinal studies of occupational co-exposure to 1,3-butadiene and styrene have shown increased standardized mortality ratios for arteriosclerotic heart disease in non-White populations;55–57 there was no data on only styrene exposure. However, in the recent case-cohort study in Iran, there was an increased risk of ischemic heart disease mortality associated with both metabolites of styrene [MA, 1.94 (1.18, 3.20) and PGA, 1.54 (1.00, 2.37)].18 ### Limitations & Future Direction This large study of VOC metabolites and blood pressure in a representative population of U.S. adults included 4,430 nonsmoking participants and evaluated 35 different VOC metabolites in urine. To minimize error from multiple testing, we prioritized nine VOC metabolites from five parent VOCs, acrolein, benzene, 1,3-butadiene, crotonaldehyde, and styrene, based on prior evidence. All VOC metabolites with >50% observations above LOD available in NHANES, however, were explored in secondary analyses. Previous studies have measured VOCs in air and grouped them by structure, functional group, or as a measure of total VOCs; here, we used stable urinary metabolites of VOC exposure, a robust measure of short-term exposure. The existing literature has been limited by the sample size of non-smokers with measured VOC metabolites. By combining 4 cycles of NHANES we increased the non-smoking sample size resulting in the largest study of VOC metabolites and BP among the general, non-smoking US population. Because metabolites from the same parent VOCs are highly correlated, we used hierarchical BKMR to confirm single-pollutant models after taking correlated VOC co-exposures into account, as well as determined the effects of VOC mixtures on BP and hypertension. These mixture analyses were only exploratory as a weighting scheme for complex survey designs has not yet been developed for BKMR. Other limitations include the cross-sectional design and potential for exposure misclassification. Data from other cohorts designed to study cardiovascular health, like the Multi-Ethnic Study of Atherosclerosis, may help elucidate differences in cardiovascular risk among VOC exposure groups, particularly with a longitudinal design to evaluate changes in BP levels and hypertension risk over time. ## PERSPECTIVES This study in NHANES found that metabolites of acrolein and 1,3-butadiene were associated with higher SBP and risk of hypertension, in the non-smoking, general US population. Additional research is needed to characterize the contribution of non-smoking related environmental VOC exposures to cardiovascular risk and to develop policy and interventions to curb preventable exposure to VOCs and further reduce the risk of CVD. ## NOVELTY AND RELEVANCE ### What is New? * This is the largest biomonitoring study of VOC exposure and elevated blood pressure in the general, non-smoking U.S. population. * VOCs present in the environment and/or endogenous VOCs are associated with higher systolic BP and greater prevalence of hypertension. ### What is Relevant? * VOC exposure, independent of cigarette smoking, may contribute to the growing prevalence of cardiovascular diseases ### Clinical/Pathophysiological Implications * Reduction of VOC exposure may improve cardiovascular health. ## Data Availability All data is available for download from the National Center for Health Statistics (NCHS), Centers for Disease Control and Prevention (CDC) website. [https://wwwn.cdc.gov/nchs/nhanes/](https://wwwn.cdc.gov/nchs/nhanes/) ## SOURCES OF FUNDING Work in the authors’ laboratories is supported by NIEHS grants P42ES023716, P42ES033719, P30ES009089, R01ES029967, and T32ES00732221. ## DISCLOSURES The authors have no conflicts of interest to disclose. ## SUPPLEMENTAL MATERIAL Tables S1-S5 Figures S1-S5 ## ACKNOWLEDGEMENTS Data and materials produced by Federal agencies are in the public domain and may be reproduced without permission. However, the authors would like to thank participants in the National Health and Nutrition Examination Surveys (NHANES). ## NON-STANDARD ABBREVIATIONS AND ACRONYMS 3HPMA : N-Acetyl-S-(3-hydroxypropyl)-L-cysteine BMI : Body mass index BKMR : Bayesian Kernel Machine Regression BP : Blood pressure CEMA : N-Acetyl-S-(2-carboxyethyl)-L-cysteine CI : Confidence intervals CVD : Cardiovascular disease DBP : Diastolic blood pressure DHBMA : N-Acetyl-S-(3,4-dihydroxybutyl)-L-cysteine eGFR : Estimated glomerular filtration HDL : High density lipoprotein cholesterol HPMMA : N-Acetyl-S-(3-hydroxypropyl-1-methyl)-L-cysteine MA : Mandelic acid MHBMA3 : N-Acetyl-S-(4-hydroxy-2-butenyl)-L-cysteine MU : t,t-Muconic Acid NHANES : National Health and Nutrition Examination Survey PGA : Phenylglyoxylic acid PIPs : Posterior inclusion probabilities PMA – N-Acetyl-S-(phenyl)-L-cysteine SBP : Systolic blood pressure VOCs : Volatile organic compounds * Received July 18, 2024. * Revision received July 18, 2024. * Accepted July 19, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## REFERENCES 1. 1.McDonald BC, de Gouw JA, Gilman JB, Jathar SH, Akherati A, Cappa CD, et al. Volatile chemical products emerging as largest petrochemical source of urban organic emissions. Science (New York, N.Y.). 2018;359:760–764 2. 2.Landrigan PJ, Fuller R, Acosta NJR, Adeyi O, Arnold R, Basu NN, et al. The lancet commission on pollution and health. Lancet (London, England). 2018;391:462–512 3. 3.Armor JN. Volatile organic compounds: An overview. Environmental catalysis. American Chemical Society; 1994:298–300. 4. 4.United States Environmental Protection Agency U. Volatile organic compounds emissions. Report on the Environment. 5. 5.Jerrett M, Burnett RT, Pope CA, Ito K, Thurston G, Krewski D, et al. Long-term ozone exposure and mortality. New England Journal of Medicine. 2009;360:1085–1095 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa0803894&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19279340&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264051000006&link_type=ISI) 6. 6.Jimenez JL, Canagaratna MR, Donahue NM, Prevot ASH, Zhang Q, Kroll JH, et al. Evolution of organic aerosols in the atmosphere. Science (New York, N.Y.). 2009;326:1525 7. 7.Danforth CG, Portier C, Ensor KB, Hopkins L, Evans B, Quist AJL, et al. Development and demonstration of a data visualization platform of short-term guidelines for ambient air levels of benzene during disaster response in houston, tx. Integrated Environmental Assessment and Management.n/a 8. 8.U.S. Environmental Protection Agency. Technical overview of volatile organic compounds. Indoor Air Quality. 2017;2019 9. 9.Kehrer JP, Biswal SS. The molecular effects of acrolein. Toxicol Sci. 2000;57:6–15 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/toxsci/57.1.6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10966506&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000165492700003&link_type=ISI) 10. 10.Esterbauer H, Schaur RJ, Zollner H. Chemistry and biochemistry of 4-hydroxynonenal, malonaldehyde and related aldehydes. Free Radic Biol Med. 1991;11:81–128 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0891-5849(91)90192-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1937131&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991GB75300012&link_type=ISI) 11. 11.Uchida K. Role of reactive aldehyde in cardiovascular diseases. Free Radic Biol Med. 2000;28:1685–1696 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0891-5849(00)00226-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10946210&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000088924800003&link_type=ISI) 12. 12.Xu X, Freeman NC, Dailey AB, Ilacqua VA, Kearney GD, Talbott EO. Association between exposure to alkylbenzenes and cardiovascular disease among national health and nutrition examination survey (nhanes) participants. International journal of occupational and environmental health. 2009;15:385–391 13. 13.Shim YH, Ock JW, Kim Y-J, Kim Y, Kim SY, Kang D. Association between heavy metals, bisphenol a, volatile organic compounds and phthalates and metabolic syndrome. Int J Environ Res Public Health. 2019;16:671 14. 14.McGraw KE, Riggs DW, Rai S, Navas-Acien A, Xie Z, Lorkiewicz P, et al. Exposure to volatile organic compounds – acrolein, 1,3-butadiene, and crotonaldehyde – is associated with vascular dysfunction. Environmental Research. 2021:110903 15. 15.McGraw KEaK, Stacey L. and Riggs, Daniel W. and Rai, Shesh N. and DeJarnett, Natasha K. and Xie, Zhengzhi and Keith, Rachel J. and Oshunbade, Adebamike A. and Hall, Michael E. and Shimbo, Daichi and Bhatnagar, Aruni. . Exposure to volatile organic compounds is associated with hypertension in black adults: The jackson heart study. . Available at SSRN. 2022 16. 16.DeJarnett N, Conklin DJ, Riggs DW, Myers JA, O’Toole TE, Hamzeh I, et al. Acrolein exposure is associated with increased cardiovascular disease risk. Journal of the American Heart Association. 2014;3 17. 17.Abplanalp W, DeJarnett N, Riggs DW, Conklin DJ, McCracken JP, Srivastava S, et al. Benzene exposure is associated with cardiovascular disease risk. PloS one. 2017;12:e0183602 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0183602&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28886060&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) 18. 18.Nalini M, Poustchi H, Bhandari D, Chang CM, Blount BC, Wang L, et al. Volatile organic compounds and mortality from ischemic heart disease: A case-cohort study. American Journal of Preventive Cardiology. 2024:100700 19. 19.Whelton PK, Carey RM, Aronow WS, Casey DE, Collins KJ, Dennison Himmelfarb C, et al. 2017 acc/aha/aapa/abc/acpm/ags/apha/ash/aspc/nma/pcna guideline for the prevention, detection, evaluation, and management of high blood pressure in adults: A report of the american college of cardiology/american heart association task force on clinical practice guidelines. Hypertension. 2018;71:e13–e115 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/HYP.0000000000000065&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29146535&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) 20. 20.U.S. Department of Health and Human Services. National health and nutrition examination survey: Plan and operations, 1999-2010. 2013 21. 21.Centers for Disease Control and Prevention (CDC). National health and nutrition examination survey data. 2011-2018 22. 22.centers for Disease Control and Prevention (CDC). Laboratory procedure manual. Antimony, Arsenic, Barium, Beryllium, Cadmium, Cesium, Cobalt, Lead, Manganese, Molybdenum, Platinum, Strontium, Thallium, Tin, Tungsten, and Uranium. 2012:63 23. 23.Centers for Disease Control and Prevention C. Laboratory procedure manual: Volatile organic compound (vocs) metabolites. 2011 24. 24.Inker LA, Eneanya ND, Coresh J, Tighiouart H, Wang D, Sang Y, et al. New creatinine- and cystatin c–based equations to estimate gfr without race. New England Journal of Medicine. 2021;385:1737–1749 [CrossRef](http://medrxiv.org/lookup/external-ref?access\_num=10.1056/NEJMOA2102953/SUPPL_FILE/NEJMOA2102953_DISCLOSURES.PDF&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) 25. 25.Centers for Disease Control and Prevention (CDC). Module 3: Weighting. National Health and Nutrition Examination Survey. 2022 26. 26.Zou G. A modified poisson regression approach to prospective studies with binary data. American Journal of Epidemiology. 2004;159:702–706 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwh090&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15033648&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000220484900009&link_type=ISI) 27. 27.Fox ER, Samdarshi TE, Musani SK, Pencina MJ, Sung JH, Bertoni AG, et al. Development and validation of risk prediction models for cardiovascular events in black adults: The jackson heart study cohort. JAMA Cardiology. 2016;1:15–25 28. 28.Forde AT, Sims M, Muntner P, Lewis T, Onwuka A, Moore K, et al. Discrimination and hypertension risk among african americans in the jackson heart study. Hypertension. 2020;76:715–723 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/HYPERTENSIONAHA.119.14492&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) 29. 29.Bobb JF, Valeri L, Claus Henn B, Christiani DC, Wright RO, Mazumdar M, et al. Bayesian kernel machine regression for estimating the health effects of multi-pollutant mixtures. Biostatistics. 2015;16:493–508 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biostatistics/kxu058&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25532525&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) 30. 30.Bobb JF. Bayesian kernel machine regression [r package bkmr version 0.2.0]. 2017 31. 31.Balakrishnan P, Beaty T, Young JH, Colantuoni E, Matsushita K. Methods to estimate underlying blood pressure: The atherosclerosis risk in communities (aric) study. PloS one. 2017;12:e0179234–e0179234 [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) 32. 32.R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2019 33. 33. Endres CJ. Nhanesa: Nhanes data retrieval. 2019 34. 34.Wickham, et al. Welcome to the tidyverse. Journal of Open Source Software. 2019;4 35. 35.Lumley T. Complex surveys: A guide to analysis using r: A guide to analysis using r. 2010 36. 36.Yang Q, Cogswell ME, Flanders WD, Hong Y, Zhang Z, Loustalot F, et al. Trends in cardiovascular health metrics and associations with all-cause and cvd mortality among us adults. JAMA. 2012;307:1273–1283 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2012.339&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22427615&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) 37. 37.Mitchell GF, Conlin PR, Dunlap ME, Lacourcière Y, Arnold JM, Ogilvie RI, et al. Aortic diameter, wall stiffness, and wave reflection in systolic hypertension. Hypertension. 2008;51:105–111 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/HYPERTENSIONAHA.107.099721&link_type=DOI) 38. 38.Conklin DJ. Acute cardiopulmonary toxicity of inhaled aldehydes: Role of trpa1. Annals of the New York Academy of Sciences. 2016;1374:59–67 39. 39.Stevens JF, Maier CS. Acrolein: Sources, metabolism, and biomolecular interactions relevant to human health and disease. Molecular nutrition & food research. 2008;52:7–25 40. 40.McGraw KE, Konkle SL, Riggs DW, Rai SN, DeJarnett N, Xie Z, et al. Exposure to volatile organic compounds is associated with hypertension in black adults: The jackson heart study. Environmental Research. 2023;223:115384 41. 41.Xu J, Niehoff NM, White AJ, Werder EJ, Sandler DP. Fossil-fuel and combustion-related air pollution and hypertension in the sister study. Environ Pollut. 2022;315:120401 42. 42.Walker VE, Degner A, Carter EW, Nicklas JA, Walker DM, Tretyakova N, et al. 1,3-butadiene metabolite 1,2,3,4 diepoxybutane induces DNA adducts and micronuclei but not t(9;22) translocations in human cells. Chemico-Biological Interactions. 2019;312:108797 43. 43.Boogaard PJ, van Sittert NJ, Megens HJ. Urinary metabolites and haemoglobin adducts as biomarkers of exposure to 1,3-butadiene: A basis for 1,3-butadiene cancer risk assessment. Chem Biol Interact. 2001;135-136:695–701 44. 44.van Sittert NJ, Megens HJ, Watson WP, Boogaard PJ. Biomarkers of exposure to 1,3-butadiene as a basis for cancer risk assessment. Toxicol Sci. 2000;56:189–202 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/toxsci/56.1.189&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10869468&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000087942900024&link_type=ISI) 45. 45.Conklin DJ, Ogunwale MA, Chen Y, Theis WS, Nantz MH, Fu X-A, et al. Electronic cigarette-generated aldehydes: The contribution of e-liquid components to their formation and the use of urinary aldehyde metabolites as biomarkers of exposure. Aerosol Sci Technol. 2018;52:1219–1232 [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) 46. 46.Grant RL, Jenkins AF. Use of in vivo and in vitro data to derive a chronic reference value for crotonaldehyde based on relative potency to acrolein. Journal of toxicology and environmental health. Part B, Critical reviews. 2015;18:327–343 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/10937404.2015.1081574&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26580244&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) 47. 47.National Center for Biotechnology Information. Crotonaldehyde.CID=447466 48. 48.Carrieri M, Bartolucci GB, Paci E, Sacco P, Pigini D, Zaratin L, et al. Validation of a radial diffusive sampler for measuring occupational exposure to 1,3-butadiene. Journal of Chromatography A. 2014;1353:114–120 49. 49.Carmella SG, Chen M, Han S, Briggs A, Jensen J, Hatsukami DK, et al. Effects of smoking cessation on eight urinary tobacco carcinogen and toxicant biomarkers. Chem Res Toxicol. 2009;22:734–741 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/tx800479s&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19317515&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265412300018&link_type=ISI) 50. 50.Lin C-Y, Lee H-L, Jung W-T, Sung F-C, Su T-C. The association between urinary levels of 1,3-butadiene metabolites, cardiovascular risk factors, microparticles, and oxidative stress products in adolescents and young adults. Journal of Hazardous Materials. 2020;396:122745 51. 51.Tevis DS, Willmore A, Bhandari D, Bowman B, Biren C, Kenwood BM, et al. Large differences in urinary benzene metabolite s-phenylmercapturic acid quantitation: A comparison of five lc–ms-ms methods. Journal of Analytical Toxicology. 2020;45:657–665 52. 52.National Health and Nutrition Examination Survey. Volatile organic compound (voc) metabolites - urine (uvoc_g). 2011-2012 Data Documentation, Codebook, and Frequencies. 2011-2012 53. 53.Mannisto T, Mendola P, Liu D, Leishear K, Sherman S, Laughon SK. Acute air pollution exposure and blood pressure at delivery among women with and without hypertension. Am J Hypertens. 2015;28:58–72 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ajh/hpu077&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24795401&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) 54. 54.Zhu Y, Zhang C, Liu D, Ha S, Kim SS, Pollack A, et al. Ambient air pollution and risk of gestational hypertension. Am J Epidemiol. 2017;186:334–343 55. 55.Matanoski GM, Santos-Burgoa C, Schwartz L. Mortality of a cohort of workers in the styrene-butadiene polymer manufacturing industry (1943-1982). Environ Health Perspect. 1990;86:107–117 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1289/ehp.9086107&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2401250&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DR35500019&link_type=ISI) 56. 56.Divine BJ, Hartman CM. Mortality update of butadiene production workers. Toxicology. 1996;113:169–181 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0300-483X(96)03442-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8901896&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VQ39100025&link_type=ISI) 57. 57.Matanoski GM, Tao X. Styrene exposure and ischemic heart disease: A case-cohort study. American Journal of Epidemiology. 2003;158:988–995 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwg247&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14607807&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F19%2F2024.07.18.24310671.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000186544200010&link_type=ISI)