, Byung-Sun Choi
, Bomi Park
, Wanhyung Lee
, Weon-Young Lee
Department of Preventive Medicine, Chung-Ang University College of Medicine, Seoul, Korea
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.
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.
| 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 |
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.
| 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.
| Characteristics | Total (n=1733) | Non-depression (n=1648) | Depression (n=85) | p-value |
|---|---|---|---|---|
| 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-value |
||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| 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-linearity |
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 frequency (weighted proportion, %). Using the Rao–Scott chi-square test.
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
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.