Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Public Health, 12 November 2025

Sec. Environmental Health and Exposome

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1648619

Association between occupational heat exposure and early renal dysfunction among Chinese petrochemical workers: a combined machine learning and WQS modeling study

  • 1Department of Preventive Medicine, School of Public Health, Fujian Medical University, Fuzhou, China
  • 2Key Laboratory of Environment and Health of Fujian Higher Education Institutes, School of Public Health, Fujian Medical University, Fuzhou, China
  • 3Institute of Population Medicine, Fujian Medical University, Fuzhou, China
  • 4Quangang Hospital, Quanzhou, China
  • 5Department of Epidemiology and Health Statistics, School of Public Health, Fujian Medical University, Fuzhou, China
  • 6The First Affiliated Hospital of Fujian Medical University, Fuzhou, China

Objective: To investigate the association between occupational heat exposure and hyperuricemia among petrochemical workers.

Methods: We retrospectively analyzed the association between workplace heat exposure and hyperuricemia by using 10 years of occupational health examination records from 2,312 petrochemical workers in Fujian Province, China. Generalized linear models (GLMs) were employed to estimate the effects of individual exposures. Weighted quantile sum (WQS) regression model was used to evaluate the combined effects of multiple occupational exposures and to identify the relative contribution of each exposure factor. A hyperuricemia risk prediction model was developed using the LightGBM machine-learning algorithm, with feature importance assessed using SHAP (SHapley Additive exPlanations) values.

Results: Occupational heat exposure was significantly associated with an increased risk of hyperuricemia (OR = 1.68, 95% CI: 1.28–2.20). In the GLM analysis, co-exposure to heat with benzene (OR = 1.93, 95% CI 1.05–3.55), H2S (OR = 3.38, 95% CI 1.94–5.88), gasoline (OR = 2.58, 95% CI 1.49–4.48), acid anhydride (OR = 2.21, 95% CI 1.09–4.48) and CO (OR = 2.14, 95% CI 1.16–3.97) further increased the risk (all p < 0.05), suggesting synergistic effects. The WQS analysis indicated that in the mixed occupational hazards exposure, heat exposure (49.2%) contributing nearly half the effect to the overall effect. The LightGBM machine learning model identified length of service, age, BMI, gender, and heat exposure as the main predictors of hyperuricemia. The SHAP analysis confirmed heat exposure as a key independent contributor alongside length of service.

Conclusion: Occupational heat exposure in petrochemical settings is significantly associated with hyperuricemia, suggesting potential early renal dysfunction risk. Integrating machine learning–based predictive models into workplace health surveillance may facilitate the early identification and management of high-risk workers. However, causal inference remains limited by the retrospective design and potential residual confounding, underscoring the need for prospective studies to validate and extend these findings.

Graphical abstract
Flowchart illustrating the relationship between occupational heat exposure, hyperuricemia, and early renal dysfunction among petrochemical workers. It includes factors like inflammation and oxidative stress leading to renal issues. Analysis involved 2312 workers using GLM and WQS models. Machine learning methods, including LASSO and SHAP values, were used for feature selection and performance evaluation.

Graphical Abstract.

1 Introduction

The petrochemical industry involves the processing and transformation of petroleum-derived chemical raw materials, where high temperatures are essential for operations like catalytic cracking and hydrocracking (1, 2). Thus, occupational heat exposure represents a significant workplace hazard. In recent years, climate change-induced extreme heat events have further exacerbated the burden of workplace heat exposure (3), which may significantly increase health risks among petrochemical workers (4, 5). The kidneys, which play a vital role in maintaining fluid and electrolyte homeostasis and in excreting metabolic waste products, are particularly vulnerable to heat stress (6). Under heat stress, peripheral vasodilation and sweating lead to substantial fluid loss, resulting in dehydration and reduced urinary excretion, which in turn elevates serum uric acid levels. Moreover, heat stress induces systemic oxidative stress and inflammatory responses that further impair renal function (7), thereby contributing to the development of both acute kidney injury (AKI) and chronic kidney disease (CKD) (8).

An increasing body of evidence has linked occupational heat exposure to renal dysfunction (911). Archived case reports from the US Occupational Safety and Health Administration (OSHA) indicate that heat-related acute kidney injury (HR-AKI) occurs across diverse occupational settings, both indoors and outdoors (12). However, research focusing on the renal health effects of heat exposure among petrochemical workers remains limited. Our preliminary analysis identified hyperuricemia as the most common abnormality related to impaired kidney function among petrochemical workers. As the end product of purine metabolism, uric acid contributes to intracellular oxidation, endothelial dysfunction, renal fibrosis, and glomerulosclerosis (13, 14). Furthermore, hyperuricemia is an established independent risk factor for the development of incident CKD and rapid renal function decline, as well as a biomarker for early-stage renal dysfunction (15, 16), being closely associated with declining eGFR, albuminuria, and kidney failure (17). Evidence also suggests that treating hyperuricemia in its early stage may delay or even prevent the onset of CKD (18). Therefore, hyperuricemia was selected as the primary indicator of early renal dysfunction in this study. Existing epidemiological studies have reported a higher prevalence of chronic kidney disease (17.9%) among residents near refineries compared to the general population (12.3%) (19), along with elevated uric acid levels in oilfield workers and local wildlife (20, 21). These findings highlight the need to investigate the association between occupational heat exposure and hyperuricemia among petrochemical workers to better inform targeted heat-related health protection strategies.

Occupational exposures in the petrochemical industry often involve co-exposure to multiple hazards, such as benzene, H2S and other chemicals, which may interact synergistically and amplify health risks. Existing evidence suggests that the combination of thermal and chemical exposures in petrochemical operations contributes to increased vulnerability to a range of diseases (22, 23). Traditional single-exposure statistical models are often inadequate for capturing the complex interactions of multiple exposures. By contrast, WQS regression can generate exposure indices that identify the most influential risk factors while accounting for interaction effects among exposures (24). Furthermore, the integration of machine learning with SHAP values offers robust methodologies for investigating nonlinear and high-dimensional relationships inherent in chemical and environmental mixtures (25, 26). Therefore, by combining multiple statistical models and machine learning algorithms, this study may yield valuable insights for developing prevention strategies and promoting the occupational health of petrochemical workers. In this study, we analyzed the association between occupational heat exposure and hyperuricemia among petrochemical workers using multiple statistical models and machine learning methodologies, with aims to provide evidence-based guidance for occupational health management and mitigate the burden of heat-related renal impairment in high-risk industrial settings.

2 Materials and methods

2.1 Study design and participants

This retrospective study was conducted using occupational health examination records collected between January 2013 and December 2022 from the Quangang Petrochemical Industrial Park (QPIP), Quanzhou, Fujian Province, China. Established in 2005, QPIP covers an area of about 30 square kilometers and has a population of around 360,000. Its industrial chain mainly involves ethylene, propylene, carbon tetrachloride, benzene and paraxylene (27).

Under the Law on the Prevention and Control of Occupational Diseases, employers are required to provide regular health examinations for workers exposed to occupational hazards, typically conducted from March to August each year. This uniform schedule reduces seasonal variability between individuals. According to the examination protocol, all workers were instructed to fast overnight before the examination and to remain fasting on the morning of the test day. In this study, historical occupational health records were obtained from the Minnan Branch of the First Affiliated Hospital of Fujian Medical University. It is the sole government-designated hospital authorized to conduct occupational health examinations in the region. All data were de-identified before delivery. All raw data underwent standardized processing, including variable recoding, value assignment, labeling, and logical consistency check. To ensure accuracy and completeness, the processed dataset was further validated through consultations with occupational health experts and enterprise management personnel.

The extracted data includes the following information: (1) Demographic characteristics, such as gender, age, and lifestyle factors (including smoking/drinking frequency); (2) General physical examination, such as height, weight, and body mass index (BMI); (3) Pre-existing chronic conditions (e.g., hyperglycemia, hypertension); and (4) Occupational exposure profile, such as length of service, types of occupational hazards (e.g., heat, benzene, methanol, gasoline, acid anhydrides, carbon monoxide (CO), hydrogen sulfide (H2S), ammonia (NH3), and noise). Routine monitoring of these occupational hazard factors was conducted by Fujian Center for Prevention and Control of Occupational Disease and Chemical Poisoning, which classified and reported exposures according to the National Occupational Hazard Detection Criteria. These records were incorporated into workers’ occupational health files. The study was approved by the Medical Ethics Committee of Fujian Medical University (Approval No.: Fuyi Medical Ethics Review No. 111).

Workers were included in the analysis if they met the following criteria: (1) Aged 18–65 years with ≥1 year of continuous frontline work in petrochemical production; and (2) With complete occupational exposure documentation. The following workers were excluded from data analysis: (1) Pre-existing diagnosis of CKD; (2) Missing key occupational exposure and outcome data in health records; (3) Temporary or short-term rotational work of < 1 year; (4) Comorbid thyroid disorders, hepatic dysfunction, or pre-existing renal impairment; and (5) History of severe cardiovascular diseases or malignant tumors. Among the occupational health examination records obtained from 2013 to 2022, not all workers participated in every examination, and to avoid within-person correlation, we structured the dataset such that each worker contributed only one observation: the first record with hyperuricemia or, if none, the most recent record without hyperuricemia. Thus, the dataset did not include repeated observations from the same individual. Finally, 2,312 petrochemical workers were included in the analysis.

2.2 Definitions of exposure variables and health outcomes

In this study, heat exposure is a categorical hazard (variable) documented in the occupational health examination records. This classification was based on occupational monitoring data collected annually by the Fujian Center for Prevention and Control of Occupational Disease and Chemical Poisoning. The classification thresholds were defined in accordance with GBZ 2.2–2007: Measurement of physical factors in the workplace Part 7. High-temperature work is defined as an operation in which the average WBGT index at the worksite is ≥25 °C during production activities (28, 29).

Body mass index (BMI): According to cut-off points for Chinese adults, overweight and obesity were defined as 24 ≤ BMI < 28 and BMI ≥ 28 kg/m2, respectively (30).

Frequency of drinking was divided into three conditions: (1) Never drinking was defined as not consuming any alcohol or alcoholic beverages for at least 6 months prior to the survey; (2) Often drinking was defined as consuming alcohol at least twice a week, with an intake greater than 50 g per occasion, for more than half a year; and (3) Occasionally drinking was between often drinking and never drinking.

Smoking frequency was divided into three conditions: (1) Never smoking was defined as not smoking for at least 6 months; (2) Often smoking was defined as smoking ≥1 cigarette per day or ≥7 cigarettes per week for at least half a year within the past year; and (3) Occasional smoking is between often smoking and never smoking.

Hyperuricemia was defined as serum uric acid ≥ 420 μmol/L (7.0 mg/dL) in males and serum uric acid ≥ 360 μmol/L (6.0 mg/dL) in females (31).

Hypertension: According to The Chinese Guidelines for Prevention and Treatment of Hypertension (2024 Revision) (32), it is defined as systolic blood pressure (SBP) ≥ 140 mmHg and/or diastolic blood pressure (DBP) ≥ 90 mmHg without any antihypertensive drugs, or those who have a previous history of hypertension and are currently taking antihypertensive drugs.

Hyperglycemia: According to The Chinese Guidelines for the Prevention and Treatment of Type 2 Diabetes criteria (33), fasting blood glucose>7.0 mmol/L; glycosylated hemoglobin level>6.5%, or self-reported formal institutional diagnosis of diabetes, or currently receiving hypoglycemic drugs.

2.3 Statistical analysis

Descriptive analysis was conducted to examine the characteristics of occupational hazards exposure and the main occupational health issues among petrochemical workers. Continuous variables were presented as mean ± standard deviation (X̅±SD), and comparisons between groups were conducted using independent sample t-tests. Categorical data are presented as frequencies and percentages (N, %), with group comparisons conducted using the chi-square test.

Three-stage stepwise adjusted GLM models were established to investigate the associations between each occupational hazard and hyperuricemia. Model 1 (unadjusted) included the single occupational hazard factor without adjusting for any covariates. Model 2 (basic adjustment) was adjusted for gender, age, and BMI. Model 3 (fully adjusted) was further adjusted for additional confounders, including length of service, hyperglycemia, hypertension, smoking, and alcohol. Covariate selection was guided by a DAG developed from prior literature and study variables. In addition, we explored potential unmeasured confounders between occupational hazards and hyperuricemia by calculating E-value. Additional stratified analyses were conducted by gender, age, BMI, and length of service to explore potential effect modification.

To evaluate the joint effects of combined exposure to heat and other occupational exposures, GLMs were performed to estimate the multiplicative interactions. Additive interactions were assessed using the relative excess risk of interaction (RERI), the attributable proportion of interaction (AP), and the synergy index (SI). The estimated interaction effects and their 95% confidence intervals (95% CI) were visualized using forest plots.

To investigate the combined effects of mixed occupational hazards on hyperuricemia and assess the effect contributions of individual factors, a WQS regression model was constructed incorporating 9 occupational hazards (34). Given the primary variables were categorical, WQS regression is well-suited for the analysis of categorical or ordinal exposure indicators. The number of quantiles was set to null. The dataset was randomly split into training (60%) and validation (40%) subsets (random seed = 1800). Each exposure variable was assigned a weight using 200 bootstrap samplings to screen factors that contribute significantly to the outcome, constraining the overall effect direction to be positive (β > 0). Weight distributions were visualized using kernel density plots. To further test the robustness of the results, we conducted negative-direction analysis of WQS model and varied the train/validation split and random seed in the WQS model.

In the initial model comparison step, nine machine learning algorithms were evaluated to identify the optimal classifier. The dataset was randomly partitioned into training (70%) and validation (30%) subsets to assess the predictive performance of each model. To prevent data leakage, all preprocessing steps, including missing data imputation, feature scaling, and feature selection (collinearity diagnosis and LASSO regression) were performed strictly within each training fold during 10-fold cross-validation, random seeds:1000. Key predictors identified by LASSO were then used to train nine machine learning models: Logistic Regression (LR), XGBoost Classifier (XGB), LightGBM Classifier (LGBM), Random Forest Classifier (RF), AdaBoost Classifier (AdaBoost), Gaussian Naive Bayes (GaussianNB), Complement Naive Bayes (ComplementNB), Multilayer Perceptron Classifier (MLP), Support Vector Classifier (SVC). Model performance was evaluated using receiver operating characteristic (ROC), area under the curve (AUC), negative predictive value, precision, recall, sensitivity, F1 score, and decision curve analysis (DCA). In further developing a hyperuricemia prediction model using the LightGBM classifier, the dataset was randomly divided into training (70%) and test (30%) subsets, and the validation set was generated from the training subset through a 10-fold cross-validation procedure. Model interpretability was enhanced using SHAP values to quantify the contribution of each predictor.

All analyses were performed using R (version 4.2.3), python (version 3.11.4), and SAS (version 9.4). Statistical tests were conducted using two-sided tests, and a p < 0.05 was considered statistically significant.

3 Results

3.1 Comparison of characteristics between hyperuricemia and control groups

In this study, among the 2,312 workers, 1,390 (60.1%) had hyperuricemia. Compared to those without hyperuricemia, affected workers had significantly higher BMI (24.23 vs. 22.98 kg/m2) and a higher proportion of males (83.96% vs. 51.84%). As shown in Table 1, workers with hyperuricemia were more likely exposed to occupational hazards, including heat, NH3, benzene, methanol, acid anhydrides, CO, and noise (all p < 0.05).

Table 1
www.frontiersin.org

Table 1. Characteristics of petrochemical workers between groups with or without hyperuricemia, n (%).

3.2 Association between occupational hazards and hyperuricemia

A three-stage stepwise GLM was used to analyze each hazard’s association with hyperuricemia. As shown in Figure 1, in unadjusted analyses, exposure to heat (OR = 2.73, 95% CI: 2.17–3.45), NH3 (OR = 2.91, 95% CI: 2.16–3.91), benzene (OR = 1.88, 95% CI: 1.50–2.37), methanol (OR = 1.99, 95% CI: 1.55–2.56), acid anhydrides (OR = 2.06, 95% CI: 1.59–2.66), CO (OR = 1.87, 95% CI: 1.49–2.34) and noise (OR = 1.27, 95% CI: 1.07–1.51) was associated with a higher risk of hyperuricemia (p < 0.05). After adjustment for age, gender, and BMI, exposure to heat (OR = 1.58, 95% CI: 1.21–2.06) remained statistically significant. In the fully adjusted model, occupational heat exposure remained a significant predictor of hyperuricemia (OR = 1.68, 95% CI: 1.28–2.20), while exposure to gasoline (OR = 0.60, 95% CI: 0.48–0.75) and H2S (OR = 0.72, 95% CI: 0.58–0.89) were associated with lower odds of hyperuricemia.

Figure 1
Forest plot showing odds ratios (OR) with 95% confidence intervals (CI) for various variables across three models. Variables include Heat, NH₃, Benzene, Methanol, H₂S, Gasoline, Acid anhydride, CO, and Noise. Model 1 shows higher ORs, especially for Heat and NH₃. Models 2 and 3 show decreasing ORs for most variables. Confidence intervals are depicted as error bars on a low to high-risk scale.

Figure 1. Odds ratios for hyperuricemia associated with occupational hazards. Model 1 (unadjusted); Model 2 (adjusted for gender, age, and BMI); and Model 3 (adjusted for gender, age, and BMI, length of service, hyperglycemia, hypertension, smoking, and alcohol consumption). *p < 0.05.

3.3 The contribution proportions of various occupational exposures to hyperuricemia

Figure 2 presents the weights of each hazard in the WQS model. Occupational heat exposure had the highest weight (49.2%), indicating it was the dominant contributor to the combined risk. Methanol and NH3 were the next largest contributors (22.1 and 11.8%, respectively).

Figure 2
Bar chart showing the weights of different exposure factors. Heat has the highest weight at 49.2%, followed by Methanol at 22.1%, NH3 at 11.8%, Benzene at 9.9%, Noise at 4.4%, Acid anhydride at 2.4%, and CO at 0.2%. Gasoline and H2S both have a weight of 0%.

Figure 2. WQS regression weights for each occupational exposure. Covariates were adjusted for gender, age, length of service, smoking, alcohol consumption, and BMI. The bar heights indicate each hazard’s contribution to the hyperuricemia risk index.

3.4 Joint effects of combined exposure to heat and other occupational hazards on hyperuricemia

In this study, we further examined the joint effects of heat and other occupational exposures on hyperuricemia. In the fully adjusted Model 3, results of multiplicative interaction analyses showed that heat combined with benzene (OR = 1.93, 95% CI 1.05–3.55), H2S (OR = 3.38, 95% CI 1.94–5.88), gasoline (OR = 2.58, 95% CI 1.49–4.48), acid anhydride (OR = 2.21, 95% CI 1.09–4.48) and CO (OR = 2.14, 95% CI 1.16–3.97) significantly increased the risk of hyperuricemia (all p < 0.05), as shown in Figure 3, while additive interaction analyses showed that these joint exposures did not exhibit significant positive additive effects (Supplementary Figure S3; Supplementary Table S3).

Figure 3
Forest plot comparing the odds ratios (OR) and 95% confidence intervals (CI) for various variables across three models. Variables include Heat#NH₃, Benzene, Methanol, H₂S, Gasoline, Acid anhydride, CO, and Noise, with corresponding ORs and CIs for each model. The plot illustrates risk levels, with <0.01 p-values significant for many comparisons, mostly showing elevated risk except for Heat#CO in Model 1. Horizontal bars indicate CI range, with markers for each model’s OR. Arrow marks low to high risk.

Figure 3. Joint effects of combined exposure to heat and other occupational hazards on hyperuricemia. Model 1 (unadjusted); Model 2 (adjusted for gender, age, and BMI); and Model 3 (adjusted for gender, age, and BMI, length of service, hyperglycemia, hypertension, smoking, and alcohol consumption). *p < 0.05.

3.5 Development of hyperuricemia risk prediction model using machine learning algorithms

3.5.1 Dataset split

To develop the hyperuricemia risk prediction model using machine learning algorithms, the dataset was split into a training set (n = 1,618) and a validation set (n = 694) in a 7:3 ratio. With the exception of noise exposure, no significant differences were observed in key demographic or baseline characteristics between the two datasets (p > 0.05), indicating good comparability. Detailed information is provided in Supplementary Table S4.

3.5.2 Feature selection using LASSO regression

Variables showing significant differences in the training set were entered into a LASSO regression to identify the most influential predictors of hyperuricemia. The optimal λ value corresponding to the minimum standard error distance was 0. 009, yielding 17 non-zero coefficient variables (Supplementary Table S6). The selected variables by LASSO regression include BMI, age, gender, smoking, alcohol drinking, length of service, benzene, noise, gasoline, and heat, used for model prediction and model construction. The results of the LASSO feature selection are shown in Figure 4.

Figure 4
Two plots labeled A and B. Plot A shows coefficient paths against log Lambda with various lines representing different coefficients. Plot B displays binomial deviance against log Lambda with red dots and error bars. Both plots have vertical dotted lines indicating reference points.

Figure 4. Feature selection of LASSO regression. (A) Shows the coefficient changes of LASSO regression; and (B) shows the fitting effect of LASSO regression.

3.5.3 Model performance comparison

The results revealed that the Random Forest model achieved the highest AUC in the training set, while the LightGBM model performed best in the validation set. However, discrepancies between training and validation results suggested potential overfitting in the Random Forest model, whereas LightGBM exhibited superior generalizability and stability. Based on these findings, the LightGBM model was selected for subsequent SHAP-based interpretability analysis. Detailed results are provided in Supplementary Tables S7–S9 and Supplementary Figures S19,S20.

3.5.4 LightGBM model training and evaluation

The ROC curve illustrates that LightGBM model achieved AUC value of 0.99 ± 0.01 for the training set, 0.86 ± 0.026 for the validation set, and 0.86 for the test set. The AUC difference between the validation and test sets was less than 10%, indicating good model generalization. The calibration curve showed good agreement between predicted probabilities and observed outcomes, with minor deviations at lower and midrange probability levels. Decision curve analysis (DCA) revealed that, across a wide range of threshold probabilities (0.1–0.7), the LightGBM model yielded a greater net benefit than both the “treat-all” and “treat-none” strategies. The confusion matrix indicated minimal risk of overfitting in the training set but a slightly elevated false-positive rate in the validation set. In addition, the KS statistics confirmed excellent discriminatory capacity and generalizability (KS = 0.60), with the optimal classification performance observed at a threshold of 0.65. Detailed results are shown in Figure 5.

Figure 5
Panel A shows the ROC curve for training data with high sensitivity and specificity. Panel B and Panel C display similar ROC curves for validation and test data. Panel D presents the learning curve for LGBMClassifier with consistent performance on training and validation sets. Panel E shows a calibration plot indicating model reliability. Panel F illustrates a test decision curve comparing different treatment strategies. Panels G and H provide confusion matrices showing predicted versus true labels. Panel I presents a KS statistic plot for test data, differentiating class distributions.

Figure 5. Consolidated performance evaluation of nine machine learning models. (A) Is ROC curve (Training set); (B) is ROC curve (Validation set); (C) is ROC curve (Test set); (D) is learning curve; (E) is calibration plots; (F) is test decision curve; (G) is confusion matrix of training set; (H) is confusion matrix of test set; and (I) is KS statistical diagram.

3.5.5 Quantifying the contribution of predictive features using SHAP values

The results showed that the SHAP values for BMI and heat were positively associated with hyperuricemia, suggesting an increased risk of hyperuricemia with higher levels of these variables, whereas length of service and age showed a negative association. The bar chart represents the relative importance of each feature in the LightGBM model, with the top five most important features were length of service, age, BMI, gender, and heat exposure. Heat exposure was strongly occupational hazards associated with an increased risk of hyperuricemia. The detailed SHAP analysis results are shown in Figure 6.

Figure 6
Two charts display SHAP values related to a model’s output impact. Chart A shows the impact of features like length of service, age, and BMI with a color gradient indicating feature values from low (blue) to high (red). Chart B displays average SHAP values, ranking features by impact magnitude, with length of service and age as most influential.

Figure 6. SHAP values quantify the contribution of predicted features. (A) The summary of SHAP values for variables in the LightGBM model. The horizontal axis represents the contribution of each feature to the model, while the vertical axis ranks the features based on the total sum of their SHAP values. The color indicates the feature value, with red corresponding to high values and blue corresponding to low values; and (B) is the histogram of feature importance for the LightGBM Model.

4 Discussion

This study examined the association between occupational heat exposure and hyperuricemia among petrochemical workers. Our findings suggested that workplace heat exposure was associated with an increased risk of hyperuricemia. Using WQS regressions, we identified heat exposure as the most influential occupational hazard contributing to hyperuricemia risk. Additionally, a machine learning-based prediction model demonstrated high accuracy in forecasting hyperuricemia risk, with SHAP analysis confirmed the strong contribution of heat exposure.

Exposure to high temperatures has emerged as a significant public health and occupational health concern, particularly due to its association with renal impairment. Prolonged exposure to process-generated heat exacerbates thermal strain and predisposes workers to dehydration and kidney injury (35). Heat stress is a potential contributing factor to chronic kidney disease of unknown etiology (CKDu) in agricultural communities, which is increasingly recognized as a major cause of chronic kidney disease (36). This CKDu, also observed in Central America, often referred to as Mesoamerican nephropathy, has been epidemiologically and experimentally linked to recurrent dehydration and chronic heat exposure (37). Mounting evidence suggests that this condition may represent a distinct form of heat stress nephropathy (HSN), potentially exacerbated by global warming (38, 39). Furthermore, a myriad of chronic kidney disease cases has surfaced among agricultural workers and other individuals engaged in manual labor across various regions of the world (40). The disease is strongly correlated with the conditions of working and living in hot environments (41). Therefore, focusing on early indicators of renal impairment, such as hyperuricemia, may provide a valuable approach for early detection and prevention of heat-related kidney injury.

Petrochemical workers were routinely exposed to multiple occupational hazards such as heat, noise, and chemical toxins during the production and processing stages. Those occupational hazards may impact workers’ health (42). Salem et al. has reported elevated levels of urea, creatinine, and uric acid among Libyan petrochemical workers, suggesting subclinical or overt kidney dysfunction (43). Consistent with these findings, our results also revealed a high prevalence of hyperuricemia among petrochemical workers, which were significantly associated with heat exposure. Several factors may contribute to this high hyperuricemia prevalence. First, the study population consists of frontline occupational workers (male-dominated), who may inherently represent a higher-risk group and generally have higher uric acid levels. Second, the petrochemical industrial park is located in a coastal region, where dietary patterns often involve high consumption of seafoods (e.g., fish, shrimp, crab, and shellfish), leading to increased purine intake and elevated serum uric acid levels. Uric acid, the final product of purine metabolism in the body, is primarily excreted by the kidneys. Under normal physiological conditions, serum uric acid levels are maintained within a certain range. However, heat exposure may disrupt this balance through dehydration, reduced renal perfusion, and impaired excretory function (44), thereby increasing the risk of hyperuricemia. In workplaces that strictly follow occupational production norms and supervision of safety departments, the likelihood of petrochemical workers experiencing rapid and severe renal impairment is relatively low, while hyperuricemia is one of the early manifestations of renal dysfunction (45). Accordingly, attention to changes in uric acid levels and the onset of hyperuricemia may hold greater practical significance for safeguarding kidney health within this occupational population.

In this study, results of multistage GLM regression suggested a significant association between heat exposure and hyperuricemia, consistent with findings from salt pan workers, sugarcane field workers, and construction workers (46, 47). Despite differences in climate and work processes, heat exposure consistently elevates uric acid levels through biological mechanisms involving inflammation, oxidative stress, altered metabolism, and cytokine activation (4850). Experimental evidence also supports this link, as animal studies have demonstrated that chronic heat exposure induces elevated serum uric acid (51). Interestingly, workers with hyperuricemia in our study tended to be younger and had shorter lengths of service. This unexpected finding may be attributed to a healthy worker survivor bias, whereby individuals unable to adapt to hazard exposures often leave high-exposure jobs earlier, and hiring restrictions might create a workforce with inherently higher health baseline levels (52, 53). Additionally, younger workers are frequently assigned to frontline positions, leading to higher exposure intensity despite shorter tenures (54), and longer-tenured workers may possess greater self-protective awareness due to extended training and job rotations, potentially reducing their risk of hyperuricemia.

In assessing joint effects of combined exposure to heat and other occupational hazards on hyperuricemia, significant multiplicative interaction effects were observed. Interestingly, the combined exposure to heat and gasoline or H2S, which both associated with higher hyperuricemia risk, despite their opposing individual effects. However, positive additive interactions were not observed, this finding consistent with Chen Y et al. (29), their research show that a multiplicative interaction rather than additive interaction between heat exposure and dust exposure in the development of hyperuricemia in steelworkers. A possible explanation of our study may lie in interaction mechanisms that modify physiological stress responses. While gasoline exposure alone may appear inversely associated with hyperuricemia, potentially due to residual confounding by work location, better ventilation, or health-based job selection, its co-occurrence with heat exposure may intensify oxidative and metabolic stress. Heat exposure promotes dehydration, reduced renal perfusion, and accumulation of reactive oxygen species (ROS) (55), which can impair uric acid excretion. Combined exposure to gasoline, containing volatile organic compounds (VOCs), may further exacerbate oxidative damage and inflammatory pathways, thereby amplifying renal tubular dysfunction (56, 57). This synergistic effect could result in a multiplicative increase in hyperuricemia risk, even when gasoline alone appears protective. Additionally, workers in high-temperature units often handle multiple chemical processes simultaneously, leading to greater overall exposure intensity and cumulative metabolic burden, which may also contribute to this combined effect. Krishnamurthy et al. found the potential synergistic effect between oil exposure and heat (58). The similar joint effect was observed for the combined exposure to heat and H₂S. This apparent paradox may reflect complex biological interactions and workplace exposure patterns. At low concentrations, H₂S can exert transient antioxidant effects by reduce the production of ROS (59). However, under conditions of heat stress, dehydration, and hypoxia, H₂S metabolism may shift toward pro-oxidative and cytotoxic pathways, enhancing oxidative stress and renal tubular injury (60). Furthermore, workers exposed to both heat and H₂S are typically engaged in frontline production units with higher overall exposure intensity and workload, which could amplify metabolic strain and renal burden (61). Therefore, combined exposure to heat and additional occupational hazards may further increase hyperuricemia risks among workers.

Traditional single-variable regression approaches often fall short in capturing the combined effects of multifactorial exposures. In contrast, the WQS model allows for estimation of cumulative effects while weighing each exposure’s relative contribution (62). Previous studies have used the WQS model to analyze the effects of mixed factors on hyperuricemia and to identify the most influential contributing factors (63). In this study, our WQS analysis identified heat and methanol exposure as the dominant contributors to hyperuricemia risk. The WQS model provides a better understanding of the relative impact of individual hazards on health outcomes and offers more precise strategies for pollutant control and management. Given the complexity of occupational exposures in petrochemical environments, predictive modeling plays an increasingly important role in occupational health surveillance. Machine learning models, particularly those that capture non-linear relationships and interactions, offer substantial advantages over traditional statistical approaches (6466). Anttila et al. has already explored the risk prediction modeling in petrochemical workers, with findings indicating that exposure to hydrocarbons in crude oil may increase the risk of kidney cancer by up to threefold (67). In this study, we comparison nine supervised learning models to predict hyperuricemia risk. After cross-validation and performance evaluation, the LightGBM model was selected for its superior balance of accuracy and generalizability.

The LightGBM model demonstrated strong predictive performance, but the AUC results suggest potential overfitting. In practical terms, this predictive model may be incorporated into occupational health surveillance systems to support dynamic risk stratification and provide early indications of workers at higher risk of hyperuricemia. Such integration could help inform timely preventive measures. Nevertheless, further external validation and real-world implementation studies are warranted before real-world application. While machine learning algorithms offer strong predictive performance, they are often criticized for their lack of transparency (68). To address this, we incorporated SHAP analysis, which has been widely adopted for interpreting machine learning models. By calculating the marginal contribution of each feature to the prediction outcome, SHAP effectively quantified the impact of individual variables (69). SHAP results confirmed that BMI and heat exposure are dominant risk factors, offering actionable insights for workplace surveillance and intervention prioritization.

Although this study identified an association between occupational heat exposure and hyperuricemia among petrochemical workers, several limitations should be acknowledged. First, this research was derived from retrospective health examination records, and we adopted a complete-case analysis approach for all main exposure and outcome variables, which may reduce sample size. The inability to infer causality is also inherently limited. Second, our study did not capture information on dietary habits, hydration patterns, genetic predisposition, medication use, metabolic factors, use of personal protective equipment, medication use, and acclimatization, which may influence the assessment of heat-hyperuricemia association. Third, heat exposure in this study primarily refers to process-generated heat from petrochemical production. Workers’ occupational health examinations were conducted during the warmer months (March to August) each year, which may potentially overestimate the association between heat exposure and hyperuricemia. Nevertheless, the uniform health examination period across all workers helps to minimize between-individual seasonal variability and thus strengthens internal comparability. However, it should be noted that the lack of year-round health examination data may limit our ability to fully characterize seasonal variations in the heat–hyperuricemia relationship. Fourth, during the additive interaction analysis in this study, some Synergy Index (SI) values appeared as NaN, likely because the odds ratio (OR) for a single exposure equaled 1, resulting in a zero denominator in the SI formula. To obtain valid estimates, we followed the software recommendation and set recode = TRUE to reverse-code the exposures. This adjustment ensured computational validity but may have altered the directionality of the additive interaction estimates (negative). Therefore, the opposite directions observed between additive and multiplicative interactions in this study should be interpreted with caution, as they may partly reflect a computational artifact rather than a genuine antagonistic interaction on the additive scale. Finally, the absence of creatinine/eGFR, albuminuria, temperature monitoring during the production process as well as relevant information regarding the workers’ production departments, has precluded a more refined assessment of the association between heat exposure (in terms of intensity and duration) and additional biomarkers of early renal dysfunction.

To address these limitations, our ongoing prospective cohort study has been designed to incorporate structured occupational and dietary questionnaires to systematically collect data on exposure histories and lifestyle factors, as well as surveying nearby residents as a control group to better evaluate the health effects of occupational exposures. Furthermore, we are also actively collaborating with government monitoring agencies to obtain accurate measurements of exposure levels and duration for various occupational hazards, and model dose–response relationships using methods such as restricted cubic splines. These efforts aim to enhance the robustness of future analyses and to inform evidence-based strategies for protecting the health of petrochemical workers.

5 Conclusion

Occupational heat exposure is a prevalent hazard in the petrochemical industry and is significantly associated with hyperuricemia. By integrating epidemiological and machine learning methods, this study quantified both the independent and combined effects of heat and other hazards on kidney health and identified key predictors of hyperuricemia risk. Findings of this study highlight the need for targeted preventive measures, including heat stress reduction and proactive monitoring high-risk workers. Future research should explore the biological mechanisms underlying these associations and externally validate predictive models in diverse occupational settings to better prevent heat-related early renal dysfunction.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by The Medical Ethics Committee of Fujian Medical University (Approval No.: Fuyi Medical Ethics Review No. 111). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

QL: Visualization, Formal analysis, Writing – original draft, Project administration, Data curation, Conceptualization. CW: Writing – review & editing. ML: Writing – review & editing. YZ: Writing – review & editing. YC: Data curation, Formal analysis, Methodology, Writing – review & editing. SD: Writing – review & editing. RX: Writing – review & editing. ZL: Writing – review & editing. WY: Project administration, Resources, Writing – review & editing. WZ: Project administration, Writing – review & editing. JX: Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Acknowledgments

The authors express their gratitude to the data custodian for the provision of workers’ occupational health examination data and the Minnan Branch of the First Affiliated Hospital of Fujian Medical University for their generous support.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2025.1648619/full#supplementary-material

References

1. Mahalkar, V, Kumar, S, and Singhal, S. Long term and short-term occupational health risks associated with petroleum industry in India. World J Advan Eng Technol Sci. (2022) 5:054–61. doi: 10.30574/wjaets.2022.5.2.0044

Crossref Full Text | Google Scholar

2. Benson, C, Dimopoulos, C, Argyropoulos, CD, Mikellidou, CV, and Boustras, GJS. Assessing the common occupational health hazards and their health risks among oil and gas workers. Saf Sci. (2021) 140:105284. doi: 10.1016/j.ssci.2021.105284

Crossref Full Text | Google Scholar

3. Di Blasi, C, Marinaccio, A, Gariazzo, C, Taiano, L, Bonafede, M, Leva, A, et al. Effects of temperatures and heatwaves on occupational injuries in the agricultural sector in Italy. Int J Environ Res Public Health. (2023) 20:2781. doi: 10.3390/ijerph20042781

PubMed Abstract | Crossref Full Text | Google Scholar

4. El Khayat, M, Halwani, DA, Hneiny, L, Alameddine, I, Haidar, MA, and Habib, RR. Impacts of climate change and heat stress on farmworkers’ health: a scoping review. Front Public Health. (2022) 10:782811. doi: 10.3389/fpubh.2022.782811

PubMed Abstract | Crossref Full Text | Google Scholar

5. Ding, R, Wang, Y, Yu, Y, Meng, X, Gong, Q, Tang, Y, et al. Effectiveness of a multifaceted intervention (TEMP) for prevention of occupational heat-related illness among outdoor workers in the power grid ondustry: a cluster randomized controlled trial. Environ Health Perspect. (2025) 133:17001. doi: 10.1289/ehp14172

Crossref Full Text | Google Scholar

6. Chu, L, Chen, K, Crowley, S, and Dubrow, R. Associations between short-term temperature exposure and kidney-related conditions in New York state: the influence of temperature metrics across four dimensions. Environ Int. (2023) 173:107783. doi: 10.1016/j.envint.2023.107783

PubMed Abstract | Crossref Full Text | Google Scholar

7. Luna-Cerón, E, Pherez-Farah, A, Krishnan-Sivadoss, I, and Guerrero-Beltrán, CE. Molecular challenges and opportunities in climate change-induced kidney diseases. Biomolecules. (2024) 14:251. doi: 10.3390/biom14030251

PubMed Abstract | Crossref Full Text | Google Scholar

8. Sanders, AP, and Jayasundara, N. Environmental pollution, heat stress and kidney health: a need for integrated assessment. Nat Rev Nephrol. (2023) 19:145–6. doi: 10.1038/s41581-023-00680-6

PubMed Abstract | Crossref Full Text | Google Scholar

9. Wang, W, Wang, F, Yang, C, Wang, J, Liang, Z, Zhang, F, et al. Associations between heat waves and chronic kidney disease in China: the modifying role of land cover. Environ Int. (2024) 186:108657. doi: 10.1016/j.envint.2024.108657

PubMed Abstract | Crossref Full Text | Google Scholar

10. Zhang, Z, Heerspink, HJL, Chertow, GM, Correa-Rotter, R, Gasparrini, A, Jongs, N, et al. Ambient heat exposure and kidney function in patients with chronic kidney disease: a post-hoc analysis of the DAPA-CKD trial. Lancet Planetary Health. (2024) 8:e225–33. doi: 10.1016/s2542-5196(24)00026-3

PubMed Abstract | Crossref Full Text | Google Scholar

11. Kim, E, Kim, H, Kim, YC, and Lee, JP. Association between extreme temperature and kidney disease in South Korea, 2003-2013: stratified by sex and age groups. Sci Total Environ. (2018) 642:800–8. doi: 10.1016/j.scitotenv.2018.06.055

PubMed Abstract | Crossref Full Text | Google Scholar

12. Shi, DS, Weaver, VM, Hodgson, MJ, and Tustin, AW. Hospitalised heat-related acute kidney injury in indoor and outdoor workers in the USA. Occup Environ Med. (2022) 79:184–91. doi: 10.1136/oemed-2021-107933

PubMed Abstract | Crossref Full Text | Google Scholar

13. Copur, S, Demiray, A, and Kanbay, M. Uric acid in metabolic syndrome: does uric acid have a definitive role? Eur J Intern Med. (2022) 103:4–12. doi: 10.1016/j.ejim.2022.04.022

PubMed Abstract | Crossref Full Text | Google Scholar

14. Cheng, Y, Zhang, H, Zheng, H, Yin, H, Wang, Y, Wang, H, et al. Association between serum uric acid/HDL-cholesterol ratio and chronic kidney disease: a cross-sectional study based on a health check-up population. BMJ Open. (2022) 12:e066243. doi: 10.1136/bmjopen-2022-066243

PubMed Abstract | Crossref Full Text | Google Scholar

15. Gaubert, M, Bardin, T, Cohen-Solal, A, Diévart, F, Fauvel, JP, Guieu, R, et al. Hyperuricemia and hypertension, coronary artery disease, kidney disease: from concept to practice. Int J Mol Sci. (2020) 21:4066. doi: 10.3390/ijms21114066

PubMed Abstract | Crossref Full Text | Google Scholar

16. Cao, X, Wu, L, and Chen, Z. The association between elevated serum uric acid level and an increased risk of renal function decline in a health checkup cohort in China. Int Urol Nephrol. (2018) 50:517–25. doi: 10.1007/s11255-017-1732-6

PubMed Abstract | Crossref Full Text | Google Scholar

17. Sharma, G, Dubey, A, Nolkha, N, and Singh, JA. Hyperuricemia, urate-lowering therapy, and kidney outcomes: a systematic review and meta-analysis. Therapeutic Advan Musculoskeletal Dis. (2021) 13:661. doi: 10.1177/1759720x211016661

PubMed Abstract | Crossref Full Text | Google Scholar

18. Li, L, Yang, C, Zhao, Y, Zeng, X, Liu, F, and Fu, P. Is hyperuricemia an independent risk factor for new-onset chronic kidney disease?: a systematic review and meta-analysis based on observational cohort studies. BMC Nephrol. (2014) 15:122. doi: 10.1186/1471-2369-15-122

PubMed Abstract | Crossref Full Text | Google Scholar

19. Okoye, O, Carnegie, E, and Mora, L. Ambient air pollution and chronic kidney disease risk in Deltan communities: a policy brief, 2023. F1000Res. (2024) 13:265. doi: 10.12688/f1000research.145904.2

PubMed Abstract | Crossref Full Text | Google Scholar

20. Husseiny, HSJZJP Sciences A. Effect of heavy metals on kidney function among petroleum industry workers in the Kurdistan region. Zanco J Pure Appl Sci. (2023) 35:230–7. doi: 10.21271/ZJPAS.35.4.22

Crossref Full Text | Google Scholar

21. Ajarem, JS, Hegazy, AK, Allam, GA, Allam, AA, Maodaa, SN, and Mahmoud, AM. Impact of petroleum industry on goats in Saudi Arabia: heavy metal accumulation, oxidative stress, and tissue injury. Environ Sci Pollut Res Int. (2023) 30:2836–49. doi: 10.1007/s11356-022-22309-0

PubMed Abstract | Crossref Full Text | Google Scholar

22. Chen, WH, Chen, ZB, Yuan, CS, Hung, CH, and Ning, SK. Investigating the differences between receptor and dispersion modeling for concentration prediction and health risk assessment of volatile organic compounds from petrochemical industrial complexes. J Environ Manag. (2016) 166:440–9. doi: 10.1016/j.jenvman.2015.10.050

PubMed Abstract | Crossref Full Text | Google Scholar

23. Fedan, JS, Thompson, JA, Russ, KA, Dey, RD, Reynolds, JS, Kashon, ML, et al. Biological effects of inhaled crude oil vapor. II. Pulmonary effects. Toxicol Appl Pharmacol. (2022) 450:116154. doi: 10.1016/j.taap.2022.116154

PubMed Abstract | Crossref Full Text | Google Scholar

24. Luo, KH, Tu, HP, Yang, CH, Yang, CC, Chen, TH, and Chuang, HY. Use of generalized weighted quantile sum regressions of tumor necrosis factor alpha and kidney function to explore joint effects of multiple metals in blood. Int J Environ Res Public Health. (2022) 19:7399. doi: 10.3390/ijerph19127399

PubMed Abstract | Crossref Full Text | Google Scholar

25. Chen, H, Wang, M, and Li, J. Exploring the association between two groups of metals with potentially opposing renal effects and renal function in middle-aged and older adults: evidence from an explainable machine learning method. Ecotoxicol Environ Saf. (2024) 269:115812. doi: 10.1016/j.ecoenv.2023.115812

PubMed Abstract | Crossref Full Text | Google Scholar

26. Wang, C, Feng, L, and Qi, Y. Explainable deep learning predictions for illness risk of mental disorders in Nanjing, China. Environ Res. (2021) 202:111740. doi: 10.1016/j.envres.2021.111740

PubMed Abstract | Crossref Full Text | Google Scholar

27. Zhang, Y, Chen, Y, Su, Q, Huang, X, Li, Q, Yang, Y, et al. The use of machine and deep learning to model the relationship between discomfort temperature and labor productivity loss among petrochemical workers. BMC Public Health. (2024) 24:3269. doi: 10.1186/s12889-024-20713-4

PubMed Abstract | Crossref Full Text | Google Scholar

28. China. MoHotPsRo, GBZT189.7–2007-measurement of physical factors in the workplace part 7: High temperature. Beijing: China Standards Press (2007).

Google Scholar

29. Chen, Y, Yang, Y, Zheng, Z, Wang, H, Wang, X, Si, Z, et al. Influence of occupational exposure on hyperuricemia in steelworkers: a nested case-control study. BMC Public Health. (2022) 22:1508. doi: 10.1186/s12889-022-13935-x

PubMed Abstract | Crossref Full Text | Google Scholar

30. Chen, K, Shen, Z, Gu, W, Lyu, Z, Qi, X, Mu, Y, et al. Prevalence of obesity and associated complications in China: a cross-sectional, real-world study in 15.8 million adults. Diabetes Obes Metab. (2023) 25:3390–9. doi: 10.1111/dom.15238

PubMed Abstract | Crossref Full Text | Google Scholar

31. Chen, D, Lu, C, Chen, K, Liu, T, Li, Y, Shan, Z, et al. Association between anthropometric indices and hyperuricemia: a nationwide study in China. Clin Rheumatol. (2024) 43:907–20. doi: 10.1007/s10067-024-06884-w

PubMed Abstract | Crossref Full Text | Google Scholar

32. Wang, JG. Chinese guidelines for the prevention and treatment of hypertension (2024 revision). J Geriatric Cardiol: JGC. (2025) 22:1–149. doi: 10.26599/1671-5411.2025.01.008

PubMed Abstract | Crossref Full Text | Google Scholar

33. Jia, W, Weng, J, Zhu, D, Ji, L, Lu, J, Zhou, Z, et al. Standards of medical care for type 2 diabetes in China 2019. Diabetes Metab Res Rev. (2019) 35:e3158. doi: 10.1002/dmrr.3158

PubMed Abstract | Crossref Full Text | Google Scholar

34. Li, S, Guo, B, Jiang, Y, Wang, X, Chen, L, Wang, X, et al. Long-term exposure to ambient PM2.5 and its components associated with diabetes: evidence from a large population-based cohort from China. Diabetes Care. (2023) 46:111–9. doi: 10.2337/dc22-1585

PubMed Abstract | Crossref Full Text | Google Scholar

35. Ebi, KL, Capon, A, Berry, P, Broderick, C, de Dear, R, Havenith, G, et al. Hot weather and heat extremes: health risks. Lancet (London, England). (2021) 398:698–708. doi: 10.1016/s0140-6736(21)01208-3

PubMed Abstract | Crossref Full Text | Google Scholar

36. Francis, A, Harhay, MN, Ong, ACM, Tummalapalli, SL, Ortiz, A, Fogo, AB, et al. Chronic kidney disease and the global public health agenda: an international consensus. Nat Rev Nephrol. (2024) 20:473–85. doi: 10.1038/s41581-024-00820-6

PubMed Abstract | Crossref Full Text | Google Scholar

37. Fischer, RSB, Vangala, C, Mandayam, S, Chavarria, D, García-Trabanino, R, Garcia, F, et al. Clinical markers to predict progression from acute to chronic kidney disease in Mesoamerican nephropathy. Kidney Int. (2018) 94:1205–16. doi: 10.1016/j.kint.2018.08.020

PubMed Abstract | Crossref Full Text | Google Scholar

38. Roncal-Jimenez, C, García-Trabanino, R, Barregard, L, Lanaspa, MA, Wesseling, C, Harra, T, et al. Heat stress nephropathy from exercise-induced uric acid crystalluria: a perspective on Mesoamerican nephropathy. Am J Kidney Dis. (2016) 67:20–30. doi: 10.1053/j.ajkd.2015.08.021

PubMed Abstract | Crossref Full Text | Google Scholar

39. Glaser, J, Lemery, J, Rajagopalan, B, Diaz, HF, García-Trabanino, R, Taduri, G, et al. Climate change and the emergent epidemic of CKD from heat stress in rural communities: the case for heat stress nephropathy. Clin J American Society Nephrol: CJASN. (2016) 11:1472–83. doi: 10.2215/cjn.13841215

PubMed Abstract | Crossref Full Text | Google Scholar

40. Redmon, JH, Levine, KE, Lebov, J, Harrington, J, and Kondash, AJ. A comparative review: chronic kidney disease of unknown etiology (CKDu) research conducted in Latin America versus Asia. Environ Res. (2021) 192:110270. doi: 10.1016/j.envres.2020.110270

PubMed Abstract | Crossref Full Text | Google Scholar

41. Johnson, RJ, Wesseling, C, and Newman, LS. Chronic kidney disease of unknown cause in agricultural communities. N Engl J Med. (2019) 380:1843–52. doi: 10.1056/NEJMra1813869

PubMed Abstract | Crossref Full Text | Google Scholar

42. O’Callaghan-Gordo, C, Orta-Martínez, M, and Kogevinas, MJEH. Health effects of non-occupational exposure to oil extraction. Environ Health. (2016) 15:56–4. doi: 10.1186/s12940-016-0140-1

PubMed Abstract | Crossref Full Text | Google Scholar

43. Salem, GM, Shaboun, S, Algamodei, YM, Almalyan, MF, Althwadi, EM, Zaid, AA, et al. Effect of occupational exposure on hematological and biochemical parameters in workers at oil and gas companies. Mediterr J Pharm Pharm Sci. (2024) 2:95–102. doi: 10.5281/zenodo.6399948

Crossref Full Text | Google Scholar

44. Wang, X, Wang, C, Wang, Z, Li, W, Sun, W, Zhang, F, et al. Antioxidant effect of taurine on chronic heat-stressed broilers. Adv Exp Med Biol. (2022) 1370:161–9. doi: 10.1007/978-3-030-93337-1_16

PubMed Abstract | Crossref Full Text | Google Scholar

45. Sampa, MB, Hossain, MN, Hoque, MR, Islam, R, Yokota, F, Nishikitani, M, et al. Blood uric acid prediction with machine learning: model development and performance comparison. JMIR Med Inform. (2020) 8:e18331. doi: 10.2196/18331

PubMed Abstract | Crossref Full Text | Google Scholar

46. Wesseling, C, Aragón, A, González, M, Weiss, I, Glaser, J, Rivard, CJ, et al. Heat stress, hydration and uric acid: a cross-sectional study in workers of three occupations in a hotspot of Mesoamerican nephropathy in Nicaragua. BMJ Open. (2016) 6:e011034. doi: 10.1136/bmjopen-2016-011034

PubMed Abstract | Crossref Full Text | Google Scholar

47. Venugopal, V, Lennqvist, R, Latha, PK, Shanmugam, R, Krishnamoorthy, M, Selvaraj, N, et al. Occupational heat stress and kidney health in salt pan workers. Kidney Int Reports. (2023) 8:1363–72. doi: 10.1016/j.ekir.2023.04.011

PubMed Abstract | Crossref Full Text | Google Scholar

48. Hansson, E, Glaser, J, Jakobsson, K, Weiss, I, Wesseling, C, Lucas, RAI, et al. Pathophysiological pathophysiological mechanisms by which heat stress potentially induces kidney inflammation and chronic kidney disease in sugarcane workers. Nutrients. (2020) 12:1639. doi: 10.3390/nu12061639

Crossref Full Text | Google Scholar

49. Butler-Dawson, J, Krisher, L, Yoder, H, Dally, M, Sorensen, C, Johnson, RJ, et al. Evaluation of heat stress and cumulative incidence of acute kidney injury in sugarcane workers in Guatemala. Int Arch Occup Environ Health. (2019) 92:977–90. doi: 10.1007/s00420-019-01426-3

PubMed Abstract | Crossref Full Text | Google Scholar

50. Schlader, ZJ, Hostler, D, Parker, MD, Pryor, RR, Lohr, JW, Johnson, BD, et al. The potential for renal injury elicited by physical work in the heat. Nutrients. (2019) 11:2087. doi: 10.3390/nu11092087

PubMed Abstract | Crossref Full Text | Google Scholar

51. Ma, B, He, X, Lu, Z, Zhang, L, Li, J, Jiang, Y, et al. Chronic heat stress affects muscle hypertrophy, muscle protein synthesis and uptake of amino acid in broilers via insulin like growth factor-mammalian target of rapamycin signal pathway. Poult Sci. (2018) 97:4150–8. doi: 10.3382/ps/pey291

PubMed Abstract | Crossref Full Text | Google Scholar

52. Neophytou, AM, Picciotto, S, Hart, JE, Garshick, E, Eisen, EA, and Laden, F. A structural approach to address the healthy-worker survivor effect in occupational cohorts: an application in the trucking industry cohort. Occup Environ Med. (2014) 71:442–7. doi: 10.1136/oemed-2013-102017

PubMed Abstract | Crossref Full Text | Google Scholar

53. Lee, WJ, Jeong, J, Bang, YJ, and Kim, YM. Healthy worker hire and survivor effects in a cohort of medical radiation workers. Int J Epidemiol. (2024) 53:dyae130. doi: 10.1093/ije/dyae130

PubMed Abstract | Crossref Full Text | Google Scholar

54. Le Moual, N, Kauffmann, F, Eisen, EA, and Kennedy, SMThe healthy worker effect in asthma: work may cause asthma, but asthma may also influence work. Am J Respir Crit Care Med. (2008) 177:4–10. doi: 10.1164/rccm.200703-415PP

Crossref Full Text | Google Scholar

55. Lacy, B, and Rahman, MS. Interactive effects of high temperature and pesticide exposure on oxidative status, apoptosis, and renin expression in kidney of goldfish: molecular and cellular mechanisms of widespread kidney damage and renin attenuation. J Applied Toxicol: JAT. (2022) 42:1787–806. doi: 10.1002/jat.4357

PubMed Abstract | Crossref Full Text | Google Scholar

56. Ekpenyong, CE, and Asuquo, AE. Recent advances in occupational and environmental health hazards of workers exposed to gasoline compounds. Int J Occup Med Environ Health. (2017) 30:1–26. doi: 10.13075/ijomeh.1896.00800

PubMed Abstract | Crossref Full Text | Google Scholar

57. Rizk, AA, Abd El-Wahab, EW, El-Marakby, FA, and El-Gazzar, RM. Assessment of oxidative stress among refueling workers in an Egyptian setting. Environ Sci Pollut Res Int. (2020) 27:18099–108. doi: 10.1007/s11356-020-08359-2

PubMed Abstract | Crossref Full Text | Google Scholar

58. Krishnamurthy, J, Engel, LS, Wang, L, Schwartz, EG, Christenbury, K, Kondrup, B, et al. Neurological symptoms associated with oil spill response exposures: results from the Deepwater horizon oil spill coast guard cohort study. Environ Int. (2019) 131:104963. doi: 10.1016/j.envint.2019.104963

PubMed Abstract | Crossref Full Text | Google Scholar

59. Li, X, Jiang, K, Ruan, Y, Zhao, S, Zhao, Y, He, Y, et al. Hydrogen sulfide and its donors: keys to unlock the chains of nonalcoholic fatty liver disease. Int J Mol Sci. (2022) 23:12202. doi: 10.3390/ijms232012202

PubMed Abstract | Crossref Full Text | Google Scholar

60. Mallardi, D, Chimienti, G, Maqoud, F, Orlando, A, Drago, S, Malerba, E, et al. The dual role of exogenous hydrogen sulfide (H(2)S) in intestinal barrier mitochondrial function: insights into cytoprotection and cytotoxicity under non-stressed conditions. Antioxidants. (2025) 14:384. doi: 10.3390/antiox14040384

PubMed Abstract | Crossref Full Text | Google Scholar

61. Chicas, RC, Wang, Y, Jennifer Weil, E, Elon, L, Xiuhtecutli, N, M, CH, et al. The impact of heat exposures on biomarkers of AKI and plasma metabolome among agricultural and non-agricultural workers. Environ Int. (2023) 180:108206. doi: 10.1016/j.envint.2023.108206

Crossref Full Text | Google Scholar

62. Yu, L, Liu, W, Wang, X, Ye, Z, Tan, Q, Qiu, W, et al. A review of practical statistical methods used in epidemiological studies to estimate the health effects of multi-pollutant mixture. Environ Pollut. (2022) 306:119356. doi: 10.1016/j.envpol.2022.119356

PubMed Abstract | Crossref Full Text | Google Scholar

63. Niu, Z, Chen, T, Duan, Z, Han, S, Shi, Y, Yu, W, et al. Associations of exposure to phthalate with serum uric acid and hyperuricemia risk, and the mediating role of systemic immune inflammation. Ecotoxicol Environ Saf. (2024) 287:117269. doi: 10.1016/j.ecoenv.2024.117269

PubMed Abstract | Crossref Full Text | Google Scholar

64. Deng, J, Ji, W, Liu, H, Li, L, Wang, Z, Hu, Y, et al. Development and validation of a machine learning-based framework for assessing metabolic-associated fatty liver disease risk. BMC Public Health. (2024) 24:2545. doi: 10.1186/s12889-024-19882-z

PubMed Abstract | Crossref Full Text | Google Scholar

65. Kristinsson, Æ, Gu, Y, Rasmussen, SM, Mølgaard, J, Haahr-Raunkjær, C, Meyhoff, CS, et al. Prediction of serious outcomes based on continuous vital sign monitoring of high-risk patients. Comput Biol Med. (2022) 147:105559. doi: 10.1016/j.compbiomed.2022.105559

PubMed Abstract | Crossref Full Text | Google Scholar

66. Deng, Y, Ma, Y, Fu, J, Wang, X, Yu, C, Lv, J, et al. A dynamic machine learning model for prediction of NAFLD in a health checkup population: a longitudinal study. Heliyon. (2023) 9:e18758. doi: 10.1016/j.heliyon.2023.e18758

PubMed Abstract | Crossref Full Text | Google Scholar

67. Anttila, A, Pokhrel, A, Heikkilä, P, Viinanen, R, and Pukkala, E. Kidney cancer risk in oil refining in Finland: a nested case-referent study. J Occup Environ Med. (2015) 57:68–72. doi: 10.1097/jom.0000000000000301

PubMed Abstract | Crossref Full Text | Google Scholar

68. London, AJ. Artificial intelligence and black-box medical decisions: accuracy versus explainability. Hast Cent Rep. (2019) 49:15–21. doi: 10.1002/hast.973

PubMed Abstract | Crossref Full Text | Google Scholar

69. Shen, M, Zhang, Y, Zhan, R, Du, T, Shen, P, Lu, X, et al. Predicting the risk of cardiovascular disease in adults exposed to heavy metals: interpretable machine learning. Ecotoxicol Environ Saf. (2025) 290:117570. doi: 10.1016/j.ecoenv.2024.117570

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: occupational heat exposure, hyperuricemia, petrochemical workers, machine learning, renal dysfunction

Citation: Li Q, Wu C, Li M, Zhang Y, Chen Y, Du S, Xu R, Lv Z, Ye W, Zheng W and Xiang J (2025) Association between occupational heat exposure and early renal dysfunction among Chinese petrochemical workers: a combined machine learning and WQS modeling study. Front. Public Health. 13:1648619. doi: 10.3389/fpubh.2025.1648619

Received: 17 June 2025; Accepted: 27 October 2025;
Published: 12 November 2025.

Edited by:

Renata Sisto, National Institute for Insurance against Accidents at Work (INAIL), Italy

Reviewed by:

Xuehai Wang, Karolinska Institutet (KI), Sweden
Laxmikant Borse, Sandip Institute of Pharmaceutical Sciences, India
P. K. Latha, Sri Ramachandra Institute of Higher Education and Research, India

Copyright © 2025 Li, Wu, Li, Zhang, Chen, Du, Xu, Lv, Ye, Zheng and Xiang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jianjun Xiang, amlhbmp1bi54aWFuZ0Bmam11LmVkdS5jbg==; Wei Zheng, emhlbmc3N3dlaUBob3RtYWlsLmNvbQ==; Weimin Ye, eXdtQGZqbXUuZWR1LmNu

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.