- 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.
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 (9–11). 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. 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. 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. 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. 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. 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. 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. 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 (48–50). 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 (64–66). 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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
28. China. MoHotPsRo, GBZT189.7–2007-measurement of physical factors in the workplace part 7: High temperature. Beijing: China Standards Press (2007).
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
Edited by:
Renata Sisto, National Institute for Insurance against Accidents at Work (INAIL), ItalyReviewed by:
Xuehai Wang, Karolinska Institutet (KI), SwedenLaxmikant 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
Rong Xu4