Skip Navigation
Skip to contents

JPMPH : Journal of Preventive Medicine and Public Health

OPEN ACCESS
SEARCH
Search

Articles

Page Path
HOME > J Prev Med Public Health > Volume 59(1); 2026 > Article
Original Article
Exposure to Benzene, Toluene, Ethylbenzene, and Xylenes and Risk of Depression: A Cross-sectional Study of a National Sample of Korean Adults
Hyun-Wook Parkorcid, Byung-Sun Choiorcid, Bomi Parkorcid, Wanhyung Leeorcid, Weon-Young Leeorcid
Journal of Preventive Medicine and Public Health 2026;59(1):95-104.
DOI: https://doi.org/10.3961/jpmph.25.522
Published online: November 6, 2025
  • 993 Views
  • 90 Download

Department of Preventive Medicine, Chung-Ang University College of Medicine, Seoul, Korea

Corresponding author: Weon-Young Lee, Department of Preventive Medicine, Chung-Ang University College of Medicine, 84 Heukseouk-ro, Dongjak-gu, Seoul 06974, Korea E-mail: wylee@cau.ac.kr
• Received: June 30, 2025   • Revised: October 6, 2025   • Accepted: October 16, 2025

Copyright © 2026 The Korean Society for Preventive Medicine

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

prev next
  • Objectives:
    Benzene, toluene, ethylbenzene, and xylenes (BTEX) are co-occurring neurotoxicants that are structurally similar aromatic hydrocarbons sharing common metabolic pathways and mechanisms of toxicity. This study investigated the effects of BTEX exposure on depression and aimed to identify the primary contributors to depression risk.
  • Methods:
    We conducted a cross-sectional analysis of 1733 adults from the 2020–2021 Korea National Health and Nutrition Examination Survey. BTEX exposure was quantified based on urinary metabolite concentrations: S-phenylmercapturic acid (SPMA) for benzene, benzylmercapturic acid for toluene, the sum of phenylglyoxylic acid and mandelic acid for ethylbenzene, and methylhippuric acid for xylenes. Depression was defined according to self-reported physician diagnosis. Logistic regression was used to estimate the risk associated with individual chemicals, while weighted quantile sum (WQS) regression was employed to assess mixture effects and identify the primary toxicant. Sensitivity analyses were performed to address potential confounding by tobacco smoke.
  • Results:
    Urinary SPMA concentrations were significantly higher in individuals with depression. Logistic regression revealed a significant association between SPMA levels and depression (odds ratio, 2.62; 95% confidence interval, 1.34 to 5.13). Although the overall BTEX mixture was not significantly associated with depression after covariate adjustment in WQS models, SPMA consistently emerged as the major contributor. The association between SPMA and depression remained robust in sensitivity analyses excluding participants exposed to tobacco smoke.
  • Conclusions:
    Benzene exposure was associated with an increased risk of depression in the general Korean adult population. Therefore, strengthened environmental regulations on benzene could help reduce the public health burden of depression.
Depression is a major global public health concern, clinically defined as a mood disorder characterized by persistent sadness, emptiness, or irritability, often accompanied by somatic and cognitive changes that substantially impair daily functioning [1]. It ranks among the leading causes of disability and disease burden worldwide [2], with an estimated global prevalence of approximately 5% among adults [3]. Thus, identifying the determinants of depression remains a critical research priority.
The etiology of depression involves a complex interplay of multiple factors. Genetic predisposition, epigenetic modifications, and psychological and social stressors have long been recognized as contributing elements [1,4,5]. In recent years, increasing attention has been directed toward environmental determinants [5]. Growing evidence implicates air pollution as a significant risk factor that may influence both the prevalence and severity of depression [5-7].
Among air pollutants, benzene and its organic derivatives—toluene, ethylbenzene, and xylenes—collectively referred to as BTEX [8], are ubiquitous chemicals emitted from industrial processes, fuel combustion, tobacco smoke, consumer products, and building materials [9,10]. Tobacco smoke, in particular, is the primary source of benzene exposure in the general population and contributes to exposure to other BTEX compounds. Their presence is especially notable indoors, where concentrations frequently exceed outdoor levels and where people spend most of their time [9].
BTEX compounds have demonstrated neurotoxic effects even at environmental exposure levels [8,11] and have been shown to induce oxidative stress and disrupt neurotransmitter systems [11-14]. These mechanisms correspond to biological pathways linking air pollutant exposure to depression risk [7]. Specifically, BTEX exposure has been associated with mood disorders, cognitive impairments, and depressive symptoms in both occupational and environmental settings [8,14-20]. Given their widespread occurrence and potential health impacts, numerous countries have implemented specific regulations to limit BTEX exposure [21,22].
However, most previous studies have focused on high-exposure scenarios, such as occupational environments or substance abuse [15-17], while more recent studies using National Health and Nutrition Examination Survey (NHANES) data have explored associations at environmental levels in the United States population [18-20]. Evidence from other populations, however, remains limited, leaving a critical knowledge gap regarding the effects of low-level, chronic exposure. To address this gap, we analyzed data from the nationally representative Korea National Health and Nutrition Examination Survey (KNHANES). Given that BTEX compounds commonly co-exist in daily environments and possess plausible additive neurotoxic potential even at environmental concentrations [8], this study aimed to assess the effects of BTEX exposure on depression risk and to identify the principal contributors within the mixture.
Study Design and Participants
KNHANES is a national surveillance system that has assessed the health and nutritional status of Koreans since 1998 [23]. To reflect increased time spent indoors during the coronavirus disease 2019 (COVID-19) pandemic, KNHANES incorporated indoor air quality (IAQ) monitoring and biomarker assessments of environmentally hazardous substances [24,25]. From July 2020 to August 2021, approximately 1200 households were sampled for IAQ measurements, stratified by season, region, and housing type. Biomarker analyses were performed on about 2000 individuals aged 19 years or older from these households who consented to biological sampling.
Of the 1966 adults who participated in the survey, we excluded those who did not complete the depression questionnaire (n=94), those with urinary creatinine levels <0.3 g/L or >3.0 g/L (n=106) [26], and those with missing covariate data (n=33). The final analytic sample included 1733 participants.
Exposure Assessment
To assess individual BTEX exposure, we used well-established urinary biomarkers of exposure: S-phenylmercapturic acid (SPMA) for benzene [27], benzylmercapturic acid (BMA) for toluene [28], phenylglyoxylic acid (PGA) plus mandelic acid (MA) for ethylbenzene [29], and methylhippuric acid (MHA) for xylenes [30]. Urinary BTEX metabolites were quantified using high-performance liquid chromatography (Nexera XR LC-20AD System; Shimadzu, Kyoto, Japan) coupled with mass spectrometry (Triple Quad API 5500; Sciex, Framingham, MA, USA). The lower limits of detection (LLODs) for SPMA, BMA, PGA, MA, 2-MHA, and 3-MHA and 4-MHA were 0.251, 0.010, 1.078, 0.252, 0.985, and 5.337 μg/L, respectively. Values below the LLOD were imputed as LLOD/√2, consistent with KNHANES analytical guidelines. Metabolite concentrations were adjusted for urinary creatinine to account for variations in dilution due to urine output [26]. Because of their right-skewed distributions, all creatinine-adjusted metabolite concentrations were log-transformed prior to statistical analyses.
Outcome Definition
Depression was defined as a self-reported physician diagnosis. During structured interviews conducted by trained personnel, participants were asked whether they had ever been diagnosed with depression by a physician.
Covariates
Data on socio-demographic, lifestyle, and health-related factors were collected through standardized interviews. Covariates were selected based on established associations with depression reported in previous studies and included age, sex, education level, occupation type, marital status, household income, residential area, smoking status, drinking status, physical activity, obesity, and subjective health status.
Age was categorized as 19–39 years, 40–59 years, or ≥60 years, and sex as male or female. Education level was classified as <high school, high school, or >high school. Occupation was grouped as white-collar (managers, professionals and related workers, clerks, service workers, sales workers), blue-collar (skilled agricultural/forestry/fishery workers, craft and related trades workers, equipment/machine operators and assemblers, elementary workers), or unemployed. The unemployed category included individuals without economic activity, such as homemakers and students, though further differentiation was not possible with available data. Marital status was categorized as never married, widowed/divorced/separated, or married/living with a partner. Household income was divided into quartiles. Residential area was classified as urban or rural. Smoking status was defined as having smoked at least 100 cigarettes in one’s lifetime. Drinking behavior was based on binge drinking within the past year, defined as consuming ≥7 shots of liquor or ≥5 cans of beer for male, and ≥5 shots of liquor or ≥3 cans of beer for female. Physical activity was dichotomized according to global health recommendations for adults (≥150 minutes of moderate-intensity or 75 minutes of vigorous-intensity aerobic activity per week, or an equivalent combination). Obesity was defined as a body mass index ≥25 kg/m². Subjective health status was categorized as very good/good, fair, or poor/very poor.
Statistical Analysis
To account for the KNHANES complex sampling design, sample weights, strata, and cluster variables were applied to all analyses, except for the weighted quantile sum (WQS) regression. Descriptive statistics were stratified by depression status. Categorical variables were expressed as frequencies (weighted proportions) and compared using the Rao–Scott χ2 test. BTEX metabolite concentrations were summarized as weighted geometric means (GMs) and geometric standard deviations (GSDs) and compared using 2-sample t-tests on log-transformed values. Spearman correlation coefficients were calculated to evaluate co-exposure patterns among BTEX metabolites.
Logistic regression was used to examine associations between each BTEX metabolite and depression, modeling exposures both as continuous variables and as categorical variables based on tertiles of their distributions. A series of progressively adjusted models were constructed: Model 1 was unadjusted; Model 2 was adjusted for age and sex; Model 3 additionally adjusted for education level, occupation, marital status, household income, residential area, smoking, drinking, physical activity, obesity, and subjective health status; and Model 4 further adjusted for the other BTEX metabolites. To assess potential non-linearity, restricted cubic spline regression with 3 knots was applied to the continuous metabolite variables [31].
Given that BTEX compounds frequently co-occur due to shared sources and may exert additive neurotoxic effects [8-10], we also evaluated the joint effect of BTEX exposure using WQS regression. This approach constructs a single weighted index representing the overall mixture effect, thereby mitigating issues of multicollinearity and multiple comparisons, while identifying the chemicals contributing most strongly to the association [32,33]. Analyses were conducted using the same covariate sets as the logistic regression models: unadjusted; adjusted for age and sex; and fully adjusted for all covariates. The dataset was randomly divided into a 40% training set (used to estimate metabolite weights via bootstrap, n=1000) and a 60% validation set (used to test the significance of the WQS index regression coefficient) [32]. To address potential instability from a single split, repeated holdout validation was performed by randomly partitioning the dataset 100 times and averaging results to obtain final estimates [33].
Because tobacco smoke is a major source of environmental BTEX exposure, we conducted sensitivity analyses restricting the sample to participants who reported neither current smoking nor environmental tobacco smoke (ETS) exposure, defined as indoor exposure at home, work, or public places within the past 7 days.
All analyses were conducted using R version 4.3.2 (R Foundation for Statistical Computing, Vienna, Austria). The ‘gWQS’ package was used for WQS regression. Statistical significance was defined as p-value <0.05 (two-sided).
Ethics Statement
We conducted a cross-sectional analysis using publicly available data from the KNHANES, which was approved by the Institutional Review Board of the Korea Disease Control and Prevention Agency (2018-01-03-2C-A; 2018-01-03-5C-A).
Participant Characteristics
Participant characteristics are summarized in Table 1. Among the 1733 participants, 85 (5.5%) reported a physician-confirmed diagnosis of depression. The depression group included a lower proportion of participants aged 40–59 years (19.7%) than the non-depression group (41.8%) and a higher proportion of females (73.2 vs. 49.5%). Individuals with depression were less likely to have education beyond high school (27.0%), to be married or living with a partner (50.8%), or to belong to the highest household income quartile (16.6%). Conversely, participants with depression were more likely to be unemployed (51.4%) and to reside in urban areas (94.3%). No significant differences were observed between the groups in smoking status, drinking behavior, physical activity, or obesity. However, a significantly higher proportion of individuals with depression reported “poor” or “very poor” subjective health status (39.2%) compared with those without depression (13.6%).
Distribution of Benzene, Toluene, Ethylbenzene, and Xylenes Metabolites
Table 2 presents the weighted distribution of urinary BTEX metabolites. The GMs (GSDs) of each metabolite (μg/g creatinine) were as follows: SPMA 0.77 (1.71), BMA 5.01 (2.37), PGA+MA 404.44 (1.70), and MHA 127.62 (2.20). These levels were substantially lower than the Biological Exposure Indices (BEIs) recommended by the American Conference of Governmental Industrial Hygienists (ACGIH), which are intended to represent occupational exposure concentrations unlikely to cause adverse effects in nearly all workers: 25 μg/g creatinine for SPMA, 150 000 μg/g creatinine for PGA+MA, and 300 000 μg/g creatinine for MHA [26]. BMA was not compared because it is not included among ACGIH-designated BEI biomarkers [26].
Only SPMA concentrations were significantly higher in the depression group (GM, 0.97; GSD, 1.70) compared with the non-depression group (GM, 0.76; GSD, 1.70) (p<0.01). No significant between-group differences were found for other BTEX metabolites.
Correlations Among Benzene, Toluene, Ethylbenzene, and Xylenes Metabolites
Figure 1 displays the correlations among BTEX metabolites. Significant positive correlations were observed across all pairs. Specifically, SPMA was weakly correlated with BMA (r=0.38) and PGA+MA (r=0.35); BMA correlated with PGA+MA (r=0.31); and PGA+MA correlated with MHA (r=0.35).
Association Between Individual Benzene, Toluene, Ethylbenzene, and Xylenes Metabolites and Depression
Table 3 summarizes the associations between BTEX metabolites and depression. SPMA exhibited a significant positive association with depression risk. In continuous models, each unit increase in ln-transformed SPMA was associated with higher odds of depression: unadjusted (odds ratio [OR], 2.27; 95% CI, 1.38 to 3.72; p<0.01), age-adjusted and sex-adjusted: (OR, 2.70; 95% CI, 1.54 to 4.74, p<0.01), and fully adjusted (OR 2.35; 95% CI, 1.21 to 4.54; p=0.01). The association remained statistically significant after additional adjustment for the other BTEX metabolites (OR, 2.62; 95% CI, 1.34 to 5.13; p<0.01). No statistically significant nonlinearity was observed for any BTEX metabolite (p for non-linearity >0.05). Categorical analyses based on tertiles of metabolite distributions also showed a consistent positive association between SPMA and depression. None of the other BTEX metabolites were significantly associated with depression risk in any of the models.
Association Between Benzene, Toluene, Ethylbenzene, and Xylenes Mixture and Depression
Figure 2 presents the results of the WQS regression analysis assessing BTEX mixture exposure in relation to depression. The BTEX mixture was significantly associated with an increased risk of depression in the unadjusted model (OR, 1.51; 95% CI, 1.12 to 2.03). However, the association was attenuated and became non-significant after adjusting for age and sex (OR, 1.34; 95% CI, 0.91 to 1.99) and in the fully adjusted model (OR, 1.12; 95% CI, 0.76 to 1.67).
Despite the loss of statistical significance in adjusted models, SPMA consistently emerged as the dominant contributor to the overall mixture effect, accounting for 69.8% of the total weight in the unadjusted model, 62.1% in the age-adjusted and sex-adjusted model, and 76.1% in the fully adjusted model.
Sensitivity Analysis
Given that both logistic regression and WQS regression identified SPMA as the principal compound associated with depression, we performed sensitivity analyses to test the robustness of this association while controlling for the potential confounding effects of tobacco smoke. Logistic regression analyses for SPMA were repeated in 2 subgroups: participants who were not current smokers (n=1472) and participants who were neither current smokers nor exposed to ETS (n=1286) (Supplemental Material 1). The significant positive association between urinary SPMA and depression persisted in both subgroups across all models, confirming the robustness of the findings.
In this nationally representative study of Korean adults, we found that higher exposure to benzene, assessed through its urinary metabolite SPMA, was significantly and consistently associated with increased odds of depression. This association remained robust after adjusting for a comprehensive set of potential confounders, including co-exposure to other BTEX compounds. Although analysis of the BTEX mixture using WQS regression did not reach statistical significance after covariate adjustment, SPMA consistently emerged as the predominant contributor to depression risk. Importantly, the association between benzene exposure and depression persisted even among individuals unexposed to tobacco smoke.
Our findings are consistent with recent research examining the effects of volatile organic compound (VOC) exposure on depression. Cross-sectional studies using data from the United States NHANES have reported significant associations between exposure to VOC mixtures and depression [18-20]. Although the specific VOCs identified as key contributors varied across studies, BTEX compounds were consistently among them [18-20]. These population-based findings extend earlier occupational health research, which has long documented that high-level exposure to organic solvents containing BTEX is associated with elevated rates of mood disorders [15-17].
Notably, NHANES studies that measured parent compounds in blood consistently identified benzene as a primary contributor [18,19], aligning with our finding that benzene is a major component of the BTEX effect on depression. In contrast, another NHANES study using urinary biomarkers identified PGA and 3,4-MHA as the principal contributors, likely explained by their higher urinary concentrations (GM: PGA, 2.56 ng/0.01 mg creatinine; 3,4-MHA, 4.67 ng/0.01 mg creatinine) compared with those observed in our study (GM: PGA, 218.01 μg/g creatinine; 3,4-MHA, 105.87 μg/g creatinine) [20]. However, that study’s biomarker panel did not include a specific urinary indicator of benzene exposure, limiting direct comparison [20]. Our findings are further supported by evidence from the UK Biobank, a large-scale cohort study that reported a significant association between long-term ambient benzene exposure and increased risk of major depression [34]. Similarly, research conducted among Gulf state residents found positive associations between blood BTEX levels (particularly benzene) and central nervous system symptoms [10]. Animal studies corroborate these findings, demonstrating that benzene exposure decreases hippocampal acetylcholine and midbrain dopamine and norepinephrine, neurotransmitters implicated in the pathophysiology of depression [35].
A key aspect of our analytical strategy was the evaluation of BTEX compounds as a chemical mixture. This approach is critical because BTEX compounds often originate from shared sources such as fuel combustion and industrial emissions [9,10], leading to correlated exposures in the general population. Moreover, due to their structural similarities, BTEX compounds are believed to exert additive neurotoxic effects even at environmental levels, as noted by the Agency for Toxic Substances and Disease Registry [8]. To address these interdependencies, we employed WQS regression, a statistical method specifically developed for mixture analysis under the assumption of additive effects [32,33]. WQS regression is particularly useful because it mitigates the multicollinearity inherent in correlated exposure data, yielding a single, interpretable estimate of the mixture’s overall effect [32,33]. Another strength of this approach is its capacity to identify the individual chemicals contributing most strongly to the association by assigning empirically derived weights to each component [32]. This enabled us not only to estimate the combined risk of BTEX exposure but also to pinpoint benzene as the principal driver.
The absence of a statistically significant BTEX mixture effect after covariate adjustment may reflect the relatively low exposure levels in our study population [20,26], which could result in subclinical neuropsychiatric effects below the threshold of detection. Additionally, using physician-diagnosed depression, as opposed to symptom-based screening tools employed in NHANES studies [18-20], may have limited sensitivity for identifying mild or subthreshold depressive states. The reduced sample size after adjustment for multiple covariates may also have constrained statistical power to detect modest associations. Nevertheless, the consistent pattern implicating benzene exposure across multiple analytical frameworks underscores its potential role in depression etiology. Larger and longitudinal studies are warranted to confirm these findings and to clarify the biological pathways through which low-level benzene exposure may contribute to depressive disorders.
This study has several notable strengths. To our knowledge, it is the first to demonstrate a potential association between BTEX exposure and depression in a large, nationally representative sample of Korean adults, enhancing the generalizability of the findings. Second, we employed specific urinary metabolites as biomarkers of exposure, providing an integrated measure of internal dose from all exposure routes and sources [27-30]. The use of SPMA—the most specific and widely validated biomarker for benzene exposure [27,36]—adds credibility to our results. However, interpretation of some biomarkers warrants caution. BMA is a metabolite of benzyl alcohol, which is both an oxygenated product of toluene and a solvent commonly used in paints, dyes, and cosmetics [28]. PGA and MA are metabolites of both ethylbenzene and styrene [29]. Thus, caution is needed when interpreting their specificity as exposure indicators. Third, the application of multiple statistical models produced consistent results, reinforcing the conclusion that benzene is the principal toxicant associated with depression.
Nonetheless, several limitations should be acknowledged. First, the cross-sectional design limits the ability to infer causality. Second, exposure assessment was based on a single spot urine sample. Because the urinary biomarkers analyzed have short biological half-lives, they primarily reflect recent BTEX exposure [37]. This limitation may introduce exposure misclassification; however, the consistency of our findings after excluding individuals exposed to tobacco smoke—a major and variable source of acute exposure [38,39]—supports the robustness of our results. Third, depression status was determined using self-reported physician diagnoses. Although this measure has been shown to correspond closely to clinical diagnoses [40], it may underestimate prevalence by missing undiagnosed cases. Future studies would benefit from incorporating standardized diagnostic interviews or validated symptom scales, such as the Patient Health Questionnaire-9, which was included in KNHANES only for limited survey years and therefore had restricted applicability. Fourth, WQS regression assumes directional homogeneity, meaning that all exposures are expected to have effects in the same direction or none at all. In our data, PGA+MA, although not statistically significant, exhibited an opposite direction of association compared with the other compounds, potentially violating this assumption. Therefore, results from mixture analyses should be interpreted with caution. Finally, despite adjusting for a comprehensive set of covariates, residual confounding due to unmeasured factors, such as genetic vulnerability or psychosocial stress [1,4,5], cannot be entirely excluded.
In conclusion, this study provides novel evidence from a nationally representative Korean population that environmental exposure to benzene is associated with an increased risk of depression, independent of tobacco smoke exposure. The robustness of this association across multiple models and sensitivity analyses underscores benzene’s potential role as a neurotoxicant contributing to depressive disorders. These findings highlight the need for strengthened environmental policies and public health interventions aimed at reducing benzene exposure as a means to help mitigate the growing mental health burden of depression. Future research should prioritize longitudinal designs with repeated biomarker assessments to establish causal relationships and better characterize long-term exposure–response dynamics over time.
Supplemental materials are available at https://doi.org/10.3961/jpmph.25.522.

Conflict of Interest

The authors have no conflicts of interest associated with the material presented in this paper.

Funding

None.

Acknowledgements

This research was supported by the Chung-Ang University Graduate Research Scholarship in 2023.

Author Contributions

Conceptualization: Park HW, Choi BS, Park B. Data curation: Park HW, Lee WY. Formal analysis: Park HW, Choi BS, Lee W, Lee WY. Funding acquisition: None. Methodology: Park B, Lee W, Lee WY. Project administration: Park B, Lee WY. Visualization: Park HW, Choi BS, Lee W. Writing – original draft: Park HW. Writing – review & editing: Park HW, Choi BS, Park B, Lee W, Lee WY.

Figure. 1.
Correlations among urinary benzene, toluene, ethylbenzene, and xylenes metabolite levels. SPMA, S-phenyl-mercapturic acid; BMA, benzyl-mercapturic acid; PGA, phenyl-glyoxylic acid; MA, mandelic acid; MHA, methylhippuric acid.
jpmph-25-522f1.jpg
Figure. 2.
ORs for the association between benzene, toluene, ethylbenzene, and xylenes (BTEX) mixture exposure and depression risk, and the weighted contributions of individual BTEX metabolites estimated by weighted quantile sum regression. (A) Model 1: unadjusted. (B) Model 2: adjusted for age and sex. (C) Model 3: adjusted for age, sex, education, occupation, marital status, household income, residential area, smoking status, drinking status, physical activity, obesity, and subjective health status. SPMA, S-phenyl-mercapturic acid; BMA, benzyl-mercapturic acid; PGA, phenyl-glyoxylic acid; MA, mandelic acid; MHA, methylhippuric acid; OR, odds ratio; CI, confidence interval.
jpmph-25-522f2.jpg
Table 1.
Baseline characteristics of the study population
Characteristics Total (n=1733) Non-depression (n=1648) Depression (n=85) p-value1
Age (y) <0.01
 19–39 411 (31.5) 385 (30.5) 26 (48.8)
 40–59 597 (40.6) 583 (41.8) 14 (19.7)
 ≥60 725 (28.0) 680 (27.8) 45 (31.4)
Sex <0.01
 Male 786 (49.2) 767 (50.5) 19 (26.8)
 Female 947 (50.8) 881 (49.5) 66 (73.2)
Education 0.03
 <High school 480 (19.4) 447 (18.9) 33 (27.4)
 High school 547 (36.5) 521 (35.9) 26 (45.7)
 >High school 706 (44.2) 680 (45.1) 26 (27.0)
Occupation <0.01
 White collar 633 (41.9) 609 (42.1) 24 (38.1)
 Blue collar 421 (23.0) 408 (23.7) 13 (10.5)
 Unemployed 679 (35.2) 631 (34.2) 48 (51.4)
Marital status 0.02
 Never married 273 (21.8) 256 (21.1) 17 (33.9)
 Widowed/Divorced/Separated 242 (10.7) 225 (10.4) 17 (15.4)
 Married/Living with partner 1218 (67.5) 1167 (68.5) 51 (50.8)
Household income (quartile) 0.01
 Lowest 299 (12.0) 278 (11.7) 21 (17.1)
 Second 427 (24.7) 404 (24.3) 23 (30.8)
 Third 453 (27.0) 427 (26.5) 26 (35.6)
 Highest 554 (36.3) 539 (37.5) 15 (16.6)
Residence area <0.01
 Urban 1411 (87.0) 1338 (86.5) 73 (94.3)
 Rural 322 (13.0) 310 (13.5) 12 (5.7)
Smoking status 0.36
 No 1086 (59.9) 1026 (59.6) 60 (65.8)
 Yes 647 (40.1) 622 (40.4) 25 (34.2)
Drinking status 0.09
 No 984 (48.8) 923 (48.1) 61 (60.0)
 Yes 749 (51.2) 725 (51.9) 24 (40.0)
Physical activity 0.54
 Insufficient 986 (53.6) 934 (53.9) 52 (49.3)
 Sufficient 747 (46.4) 714 (46.1) 33 (50.7)
Obesity 0.42
 No 1103 (64.0) 1048 (63.8) 55 (68.6)
 Yes 630 (36.0) 600 (36.2) 30 (31.4)
Subjective health status <0.01
 Very good/Good 543 (31.3) 525 (32.0) 18 (19.8)
 Fair 899 (53.7) 865 (54.4) 34 (41.1)
 Bad/Very bad 291 (15.0) 258 (13.6) 33 (39.2)

Values are presented as frequency (weighted proportion, %).

1 Using the Rao–Scott chi-square test.

Table 2.
Weighted distribution of urinary creatinine-adjusted concentrations of BTEX metabolites (μg/g creatinine)
Variables ACGIH BEI Total (n=1733)
Non-depression (n=1648)
Depression (n=85)
p-value1
GM (GSD) T1 T2 GM (GSD) T1 T2 GM (GSD) T1 T2
SPMA 25 0.77 (1.71) 0.59 0.96 0.76 (1.70) 0.58 0.96 0.97 (1.70) 0.73 1.35 <0.01
BMA N/A 5.01 (2.37) 3.21 6.71 4.97 (2.36) 3.21 6.69 5.69 (2.51) 3.46 6.89 0.29
PGA+MA 150 000 404.44 (1.70) 325.67 477.14 405.51 (1.69) 323.07 479.54 386.40 (1.84) 357.83 468.95 0.53
MHA 300 000 127.62 (2.20) 85.68 143.70 127.60 (2.21) 85.49 141.48 127.98 (2.01) 92.57 156.17 0.98

BTEX, benzene, toluene, ethylbenzene, and xylenes; ACGIH, American Conference of Governmental Industrial Hygienists; BEI, biological exposure indices; SPMA, S-phenyl-mercapturic acid; BMA, benzyl-mercapturic acid; PGA, phenyl-glyoxylic acid; MA, mandelic acid; MHA, methylhippuric acid; GM, geometric mean; GSD, geometric standard deviation; T1, first tertile; T2, second tertile; N/A, not available

1 Calculated for log-transformed metabolite concentrations between the non-depression and depression groups using the 2-sample t-test.

Table 3.
Logistic regression analysis of the association between BTEX exposure and depression1
Variables Continuous p-value p-value for non-linearity2 1st tertile 2nd tertile p-value 3rd tertile p-value
SPMA
 Model 1 2.27 (1.38, 3.72) <0.01 0.53 1.00 (reference) 2.71 (1.19, 6.18) 0.02 3.25 (1.51, 7.02) <0.01
 Model 2 2.70 (1.54, 4.74) <0.01 0.27 1.00 (reference) 3.44 (1.42, 8.33) <0.01 4.25 (1.76, 10.27) <0.01
 Model 3 2.35 (1.21, 4.54) 0.01 0.30 1.00 (reference) 3.65 (1.34, 9.96) 0.01 3.50 (1.26, 9.77) 0.02
 Model 4 2.62 (1.34, 5.13) <0.01 0.21 1.00 (reference) 3.54 (1.38, 9.05) <0.01 3.99 (1.51, 10.56) <0.01
BMA
 Model 1 1.19 (0.88, 1.60) 0.26 0.59 1.00 (reference) 0.98 (0.51, 1.88) 0.95 1.08 (0.53, 2.17) 0.84
 Model 2 1.12 (0.80, 1.56) 0.51 0.63 1.00 (reference) 0.88 (0.46, 1.68) 0.69 0.86 (0.39, 1.93) 0.72
 Model 3 1.05 (0.76, 1.44) 0.77 0.58 1.00 (reference) 0.83 (0.39, 1.77) 0.63 0.73 (0.32, 1.64) 0.44
 Model 4 1.00 (0.72, 1.40) 1.00 0.59 1.00 (reference) 0.70 (0.34, 1.45) 0.34 0.62 (0.29, 1.37) 0.21
PGA+MA
 Model 1 0.84 (0.50, 1.43) 0.52 0.86 1.00 (reference) 2.04 (1.05, 3.97) 0.04 1.05 (0.50, 2.20) 0.89
 Model 2 0.79 (0.41, 1.52) 0.49 0.99 1.00 (reference) 1.90 (0.96, 3.76) 0.06 1.07 (0.48, 2.38) 0.86
 Model 3 0.70 (0.35, 1.40) 0.31 0.84 1.00 (reference) 1.58 (0.78, 3.17) 0.20 0.85 (0.38, 1.90) 0.69
 Model 4 0.55 (0.26, 1.18) 0.12 0.90 1.00 (reference) 1.32 (0.67, 2.62) 0.42 0.63 (0.29, 1.37) 0.24
MHA
 Model 1 1.00 (0.73, 1.39) 0.98 0.29 1.00 (reference) 1.18 (0.60, 2.33) 0.63 1.59 (0.79, 3.20) 0.19
 Model 2 1.12 (0.79, 1.59) 0.54 0.53 1.00 (reference) 1.05 (0.54, 2.06) 0.87 1.82 (0.91, 3.64) 0.09
 Model 3 1.01 (0.70, 1.45) 0.98 0.30 1.00 (reference) 1.27 (0.63, 2.55) 0.50 1.78 (0.87, 3.68) 0.12
 Model 4 1.10 (0.72, 1.66) 0.66 0.28 1.00 (reference) 1.19 (0.62, 2.28) 0.59 1.76 (0.85, 3.67) 0.13

Values are presented as odds ratio (95% confidence interval).

BTEX, benzene, toluene, ethylbenzene, and xylenes; SPMA, S-phenyl-mercapturic acid; BMA, benzyl-mercapturic acid; PGA, phenyl-glyoxylic acid; MA, mandelic acid; MHA, methylhippuric acid.

1 Model 1 is unadjusted; Model 2 is adjusted for age and sex; Model 3 is adjusted for age, sex, education, occupation, marital status, household income, residence area, smoking status, drinking status, physical activity, obesity, and subjective health status; Model 4 is further adjusted for the other BTEX metabolites in addition to all covariates in Model 3.

2 Using restricted cubic spline regression with 3 knots applied to continuous metabolite concentrations.

Figure & Data

References

    Citations

    Citations to this article as recorded by  

      • PubReader PubReader
      • Cite
        CITE
        export Copy
        Close
      • XML DownloadXML Download
      Figure
      • 0
      • 1
      Exposure to Benzene, Toluene, Ethylbenzene, and Xylenes and Risk of Depression: A Cross-sectional Study of a National Sample of Korean Adults
      Image Image
      Figure. 1. Correlations among urinary benzene, toluene, ethylbenzene, and xylenes metabolite levels. SPMA, S-phenyl-mercapturic acid; BMA, benzyl-mercapturic acid; PGA, phenyl-glyoxylic acid; MA, mandelic acid; MHA, methylhippuric acid.
      Figure. 2. ORs for the association between benzene, toluene, ethylbenzene, and xylenes (BTEX) mixture exposure and depression risk, and the weighted contributions of individual BTEX metabolites estimated by weighted quantile sum regression. (A) Model 1: unadjusted. (B) Model 2: adjusted for age and sex. (C) Model 3: adjusted for age, sex, education, occupation, marital status, household income, residential area, smoking status, drinking status, physical activity, obesity, and subjective health status. SPMA, S-phenyl-mercapturic acid; BMA, benzyl-mercapturic acid; PGA, phenyl-glyoxylic acid; MA, mandelic acid; MHA, methylhippuric acid; OR, odds ratio; CI, confidence interval.
      Exposure to Benzene, Toluene, Ethylbenzene, and Xylenes and Risk of Depression: A Cross-sectional Study of a National Sample of Korean Adults
      Characteristics Total (n=1733) Non-depression (n=1648) Depression (n=85) p-value1
      Age (y) <0.01
       19–39 411 (31.5) 385 (30.5) 26 (48.8)
       40–59 597 (40.6) 583 (41.8) 14 (19.7)
       ≥60 725 (28.0) 680 (27.8) 45 (31.4)
      Sex <0.01
       Male 786 (49.2) 767 (50.5) 19 (26.8)
       Female 947 (50.8) 881 (49.5) 66 (73.2)
      Education 0.03
       <High school 480 (19.4) 447 (18.9) 33 (27.4)
       High school 547 (36.5) 521 (35.9) 26 (45.7)
       >High school 706 (44.2) 680 (45.1) 26 (27.0)
      Occupation <0.01
       White collar 633 (41.9) 609 (42.1) 24 (38.1)
       Blue collar 421 (23.0) 408 (23.7) 13 (10.5)
       Unemployed 679 (35.2) 631 (34.2) 48 (51.4)
      Marital status 0.02
       Never married 273 (21.8) 256 (21.1) 17 (33.9)
       Widowed/Divorced/Separated 242 (10.7) 225 (10.4) 17 (15.4)
       Married/Living with partner 1218 (67.5) 1167 (68.5) 51 (50.8)
      Household income (quartile) 0.01
       Lowest 299 (12.0) 278 (11.7) 21 (17.1)
       Second 427 (24.7) 404 (24.3) 23 (30.8)
       Third 453 (27.0) 427 (26.5) 26 (35.6)
       Highest 554 (36.3) 539 (37.5) 15 (16.6)
      Residence area <0.01
       Urban 1411 (87.0) 1338 (86.5) 73 (94.3)
       Rural 322 (13.0) 310 (13.5) 12 (5.7)
      Smoking status 0.36
       No 1086 (59.9) 1026 (59.6) 60 (65.8)
       Yes 647 (40.1) 622 (40.4) 25 (34.2)
      Drinking status 0.09
       No 984 (48.8) 923 (48.1) 61 (60.0)
       Yes 749 (51.2) 725 (51.9) 24 (40.0)
      Physical activity 0.54
       Insufficient 986 (53.6) 934 (53.9) 52 (49.3)
       Sufficient 747 (46.4) 714 (46.1) 33 (50.7)
      Obesity 0.42
       No 1103 (64.0) 1048 (63.8) 55 (68.6)
       Yes 630 (36.0) 600 (36.2) 30 (31.4)
      Subjective health status <0.01
       Very good/Good 543 (31.3) 525 (32.0) 18 (19.8)
       Fair 899 (53.7) 865 (54.4) 34 (41.1)
       Bad/Very bad 291 (15.0) 258 (13.6) 33 (39.2)
      Variables ACGIH BEI Total (n=1733)
      Non-depression (n=1648)
      Depression (n=85)
      p-value1
      GM (GSD) T1 T2 GM (GSD) T1 T2 GM (GSD) T1 T2
      SPMA 25 0.77 (1.71) 0.59 0.96 0.76 (1.70) 0.58 0.96 0.97 (1.70) 0.73 1.35 <0.01
      BMA N/A 5.01 (2.37) 3.21 6.71 4.97 (2.36) 3.21 6.69 5.69 (2.51) 3.46 6.89 0.29
      PGA+MA 150 000 404.44 (1.70) 325.67 477.14 405.51 (1.69) 323.07 479.54 386.40 (1.84) 357.83 468.95 0.53
      MHA 300 000 127.62 (2.20) 85.68 143.70 127.60 (2.21) 85.49 141.48 127.98 (2.01) 92.57 156.17 0.98
      Variables Continuous p-value p-value for non-linearity2 1st tertile 2nd tertile p-value 3rd tertile p-value
      SPMA
       Model 1 2.27 (1.38, 3.72) <0.01 0.53 1.00 (reference) 2.71 (1.19, 6.18) 0.02 3.25 (1.51, 7.02) <0.01
       Model 2 2.70 (1.54, 4.74) <0.01 0.27 1.00 (reference) 3.44 (1.42, 8.33) <0.01 4.25 (1.76, 10.27) <0.01
       Model 3 2.35 (1.21, 4.54) 0.01 0.30 1.00 (reference) 3.65 (1.34, 9.96) 0.01 3.50 (1.26, 9.77) 0.02
       Model 4 2.62 (1.34, 5.13) <0.01 0.21 1.00 (reference) 3.54 (1.38, 9.05) <0.01 3.99 (1.51, 10.56) <0.01
      BMA
       Model 1 1.19 (0.88, 1.60) 0.26 0.59 1.00 (reference) 0.98 (0.51, 1.88) 0.95 1.08 (0.53, 2.17) 0.84
       Model 2 1.12 (0.80, 1.56) 0.51 0.63 1.00 (reference) 0.88 (0.46, 1.68) 0.69 0.86 (0.39, 1.93) 0.72
       Model 3 1.05 (0.76, 1.44) 0.77 0.58 1.00 (reference) 0.83 (0.39, 1.77) 0.63 0.73 (0.32, 1.64) 0.44
       Model 4 1.00 (0.72, 1.40) 1.00 0.59 1.00 (reference) 0.70 (0.34, 1.45) 0.34 0.62 (0.29, 1.37) 0.21
      PGA+MA
       Model 1 0.84 (0.50, 1.43) 0.52 0.86 1.00 (reference) 2.04 (1.05, 3.97) 0.04 1.05 (0.50, 2.20) 0.89
       Model 2 0.79 (0.41, 1.52) 0.49 0.99 1.00 (reference) 1.90 (0.96, 3.76) 0.06 1.07 (0.48, 2.38) 0.86
       Model 3 0.70 (0.35, 1.40) 0.31 0.84 1.00 (reference) 1.58 (0.78, 3.17) 0.20 0.85 (0.38, 1.90) 0.69
       Model 4 0.55 (0.26, 1.18) 0.12 0.90 1.00 (reference) 1.32 (0.67, 2.62) 0.42 0.63 (0.29, 1.37) 0.24
      MHA
       Model 1 1.00 (0.73, 1.39) 0.98 0.29 1.00 (reference) 1.18 (0.60, 2.33) 0.63 1.59 (0.79, 3.20) 0.19
       Model 2 1.12 (0.79, 1.59) 0.54 0.53 1.00 (reference) 1.05 (0.54, 2.06) 0.87 1.82 (0.91, 3.64) 0.09
       Model 3 1.01 (0.70, 1.45) 0.98 0.30 1.00 (reference) 1.27 (0.63, 2.55) 0.50 1.78 (0.87, 3.68) 0.12
       Model 4 1.10 (0.72, 1.66) 0.66 0.28 1.00 (reference) 1.19 (0.62, 2.28) 0.59 1.76 (0.85, 3.67) 0.13
      Table 1. Baseline characteristics of the study population

      Values are presented as frequency (weighted proportion, %).

      Using the Rao–Scott chi-square test.

      Table 2. Weighted distribution of urinary creatinine-adjusted concentrations of BTEX metabolites (μg/g creatinine)

      BTEX, benzene, toluene, ethylbenzene, and xylenes; ACGIH, American Conference of Governmental Industrial Hygienists; BEI, biological exposure indices; SPMA, S-phenyl-mercapturic acid; BMA, benzyl-mercapturic acid; PGA, phenyl-glyoxylic acid; MA, mandelic acid; MHA, methylhippuric acid; GM, geometric mean; GSD, geometric standard deviation; T1, first tertile; T2, second tertile; N/A, not available

      Calculated for log-transformed metabolite concentrations between the non-depression and depression groups using the 2-sample t-test.

      Table 3. Logistic regression analysis of the association between BTEX exposure and depression1

      Values are presented as odds ratio (95% confidence interval).

      BTEX, benzene, toluene, ethylbenzene, and xylenes; SPMA, S-phenyl-mercapturic acid; BMA, benzyl-mercapturic acid; PGA, phenyl-glyoxylic acid; MA, mandelic acid; MHA, methylhippuric acid.

      Model 1 is unadjusted; Model 2 is adjusted for age and sex; Model 3 is adjusted for age, sex, education, occupation, marital status, household income, residence area, smoking status, drinking status, physical activity, obesity, and subjective health status; Model 4 is further adjusted for the other BTEX metabolites in addition to all covariates in Model 3.

      Using restricted cubic spline regression with 3 knots applied to continuous metabolite concentrations.


      JPMPH : Journal of Preventive Medicine and Public Health
      TOP