Next Article in Journal
Features Linked to Treatment Outcomes in Behavioral Addictions and Related Disorders
Previous Article in Journal
Which Aspects of Work Safety Satisfaction Are Important to Mental Health of Healthcare Workers during COVID-19 Pandemic in Poland?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Latent Trajectories of Haematological, Hepatic, and Renal Profiles after Oil Spill Exposure: A Longitudinal Analysis

by
Benjamin Atta Owusu
1,2,3,4,
Apiradee Lim
1,3,
Nitinun Pongsiri
1,2,
Chanthip Intawong
5,
Sunthorn Rheanpumikankit
6,
Saijit Suksri
6 and
Thammasin Ingviya
3,4,7,*
1
Department of Mathematics and Computer Science, Faculty of Science and Technology, Prince of Songkla University, Pattani Campus, Pattani 94000, Thailand
2
Multidisciplinary Research and Innovation Centre, Kumasi AOK569, Ghana
3
Air Pollution and Health Effect Research Center, Prince of Songkla University, Hat Yai Campus, Songkhla 90110, Thailand
4
Division of Digital Innovation and Data Analytics, Faculty of Medicine, Prince of Songkla University, Hat Yai Campus, Songkhla 90110, Thailand
5
Occupational Medicine Department, Rayong Hospital, Rayong 21000, Thailand
6
Rayong Provincial Public Health Office, Rayong 21000, Thailand
7
Department of Family and Preventive Medicine, Faculty of Medicine, Prince of Songkla University, Hat Yai Campus, Songkhla 90110, Thailand
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(4), 2871; https://doi.org/10.3390/ijerph20042871
Submission received: 25 November 2022 / Revised: 15 January 2023 / Accepted: 28 January 2023 / Published: 6 February 2023
(This article belongs to the Section Environmental Health)

Abstract

:
Exposure to polycyclic aromatic hydrocarbons (PAHs) and volatile organic compounds (VOCs) in crude oil has carcinogenic effects on various organ systems. This longitudinal cohort study examined the effects of oil spill exposure on the haematological, hepatic, and renal profiles of Rayong oil spill clean-up workers. The sample included 869 clean-up workers from the Rayong oil spill. Latent class mixture models were used to investigate and classify the longitudinal trajectories and trends of the haematological, hepatic, and renal indices. Subgroup analysis was used to evaluate the association between the urinary metabolites of PAHs and VOCs and haematological, hepatic, and renal parameters. Most clean-up workers (97.6%) had increasing levels of white blood cells (WBCs) (0.03 × 103 cells/µL), 94.90% of the workers had a significantly increasing trend of blood urea nitrogen (0.31 mg/dL per year), and 87.20% had a significantly increasing trend of serum creatinine (0.01 mg/dL per year). A high–decreasing trend of WBCs was seen in 2.42% (−0.73 × 103 per year). Post-exposure changes in haematological, renal, and hepatic profiles are present in workers exposed to the Rayong oil spill. This indicates possible long-term health complications and worsening renal function after exposure to PAHs and VOCs in crude oil.

1. Introduction

Oil spill disasters have released toxic chemicals, such as polycyclic aromatic hydrocarbons (PAHs) and volatile organic compounds (VOCs), into marine and terrestrial ecosystems for decades following the invention of liquid propellants [1,2]. In most oil spills, local citizens voluntarily participate in clean-up activities and become exposed to various levels of toxicants. PAHs and VOCs in crude oil, such as benzo[a]pyrene, are classified as Group I human carcinogens by the International Agency for Research on Cancer [3,4]. The intake of these compounds, primarily by inhalation, may induce severe adverse health effects, including alterations in haematological, hepatic, and renal functions [5,6].
The quantification of urinary metabolites of PAH and VOCs to estimate exposure levels has been investigated previously [7,8]. The levels of 1-hydroxypyrene-glucuronide (1-OHPG) and trans,trans-muconic acid (t,t-MA) have been widely used to assess environmental exposure to PAHs [7,8] and VOCs (e.g., benzene), respectively [9,10,11]. The levels of urine metabolites represent quantifiable levels of exposure from multiple routes and sources, for instance, cigarette smoke and oil spills [12,13].
The complete blood count (CBC) and serum levels of renal and liver biomarkers are commonly used to assess the effects of crude oil exposure in humans and other organisms. Several studies have used haematological, hepatic, and renal parameters to assess the effects of exposure to PAH and VOC on human health [6,14,15,16]. However, the effects on these parameters vary according to the level and duration of exposure, as well as by the influence of other environmental factors.
The Rayong oil spill of 2013 was caused by a ruptured pipeline, whereby more than 50,000 L of crude oil spilled into the Gulf of Thailand. Until January 2022, this was the most recent marine disaster in Thailand. An observational baseline study conducted by Ingviya et al. [17] assessed urinary 1-OHPG and t,t-MA concentrations to determine the internal exposure to PAHs and VOCs among clean-up workers of the Rayong oil spill on different clean-up days. The study showed that workers who participated during the initial days of clean-up had significantly higher exposure, as indicated by higher levels of 1-OHPG, while the t,t-MA concentration was undetectable in the majority of clean-up workers throughout the clean-up period [17]. Similarly, in an observational study, Owusu et al. [18] assessed endpoint changes in haematological, hepatic, and renal indices among 570 Rayong oil spill clean-up workers. The exposure levels were differentiated based on their work on different clean-up days.
In this study, we hypothesised that exposure to the Rayong oil spill induced different nonlinear trajectories of haematological, hepatic, and renal profiles. We aimed to examine the long-term effects of toxicants (PAHs and VOCs) by assessing the 5-year trajectories of haematological, hepatic, and renal indices after the 2013 Rayong oil spill. Primarily, we examined the presence of distinct longitudinal trajectories of these indices among clean-up workers and their associated factors, including age. Finally, long-term changes in the haematological, hepatic, and renal indices by the exposure levels of PAHs and VOCs were assessed using urinary concentrations of the toxicant metabolites. Assessing these trajectories can guide and encourage the implementation of health follow-up protocols for workers who clean oil spills.

2. Materials and Methods

2.1. Study Design and Setting

In this retrospective longitudinal cohort study, we analysed haematological, liver, and kidney function parameters, measured at baseline and during a 5-year follow-up after the 2013 Rayong oil spill. The exposure data of the cohort were retrieved from the baseline post-shift urinary concentrations of 1-OHPG and t,t-MA, quantified by Ingviya et al. [17]. Data were retrieved from the record forms initially collected by Rayong Hospital as part of health services.

2.2. Ethics Approval and Informed Consent

Clean-up workers verbally consented to the use of their data for further analysis to improve the quality of health services and follow-up protocols.

2.3. Study Population and Measurements

During and after the Rayong oil spill clean-up activities, the Rayong Provincial Health Office collaborated with Rayong Hospital to undertake a 5-year health surveillance of clean-up workers. After the oil spill, the workers were invited to undergo annual health assessments, including laboratory tests, CBCs, and hepatic and renal function parameters. The dataset included 869 Rayong oil spill clean-up workers that attended at least one follow-up visit between 2014 and 2018. A subgroup analysis was conducted on 169 workers who provided urine samples to determine 1-OHPG and t,t-MA urinary concentrations. Demographic information including age, sex, and background occupation was collected. Urinary cotinine levels were measured, and a cut-off of 50 ng/mL used to differentiate between non-smokers and smokers [19]. Data on specific clean-up activities were also collected. Haematological, hepatic, and renal parameters at baseline and follow-up were retrieved from the electronic records of the Rayong Provincial Health Office and Rayong Hospital.

2.4. Quantification of 1-OHPG and t,t-MA Concentrations

The concentrations of 1-OHPG and t,t-MA were quantified in urine samples collected during clean-up in 2013. Immunoaffinity chromatography and synchronous fluorescence spectroscopy were used to measure the 1-OHPG concentration at Paul Strickland laboratory at John Hopkins University (limit of detection 0.04 pmol/mL; coefficient of variation 5.6%) [17]. Baseline 1-OHPG concentration was categorised as high (>5.0 pmol/mL), moderate (1.0–5.0 pmol/mL), or low (<1.0 pmol/mL), based on a study by Kang et al. [20].
Urinary t,t-MA concentration was analysed using high-performance liquid chromatography with fluorescence detection in three laboratories, including Rayong Hospital and two other private laboratories. The agreement between the t,t-MA measurements reported by the three laboratories was 99.99% [21]. The t,t-MA concentration was classified as detectable or undetectable because a large proportion of the samples had levels below the limit of detection (0.01 mg/dL) of the assay.

2.5. Measurement of Haematological, Hepatic, and Renal Parameters

Approved medical protocols in a standardised laboratory at the Royal Hospital were used for analysis. The CBC parameters, including haemoglobin, haematocrit, red blood cell count (RBC), white blood cell count (WBC), absolute neutrophil count (ANC), and platelet count, were measured using an automatic CBC analyser. Hepatic [serum aspartate transaminase (AST) and alanine transaminase (ALT)] and renal indices [blood urea nitrogen (BUN) and serum creatinine (Cr)] were determined using enzymatic assay techniques.

2.6. Statistical Analysis

Descriptive statistics (mean, standard deviation, and percentages) were used for demographic data of the oil spill clean-up workers and to summarise the haematological and hepatic profiles. Generalised mixed-effects models accounted for the expected within-subject correlation between the repeated measures. Latent class mixture models (LCMMs) were used to investigate and classify the longitudinal trajectories and trends of the haematological, hepatic, and renal indices, while adjusting for age [22,23,24]. The “lcmm” package in R was used for the latent class trajectory analysis [25]. Polynomial specifications as functions of time, with a number of latent classes ranging from one to three, were used in the LCMMs to assess the different trends. The optimal number of latent classes for each model was determined using the Bayesian information criterion (BIC) and posterior probabilities above a 0.7 threshold [25]. The latent classes for each haematological, hepatic, and renal index were labelled using the observed mean at baseline (high, low, and normal) and their adjusted trends during the study (stable, increasing, and decreasing). The definitions of high, low, and normal concentrations were based on the standard medical reference ranges for each index at baseline [26]. We examined the yearly class-specific observed and predicted means of each haematological, hepatic, and renal index. Furthermore, we examined the effects of age and the latent longitudinal trajectories of the indices by assessing age-specific trajectories.
A subgroup analysis of 169 oil spill clean-up workers was performed to determine the effects of baseline urinary concentrations of 1-OHPG and t,t-MA and smoking status on the overall trajectories and trends of haematological and hepatorenal profiles. Age-centred generalised linear mixed models were used for the subgroup analysis. The characteristics of the subgroup, including age, gender, and background occupation, were compared to the characteristics of the full cohort (869 workers) to assess the possibility of selection bias. A two-sample t-test was used to test the hypothesis that the average age of the full cohort was equal to the average age of the subgroup. Similarly, a Chi-square test with Yates’ continuity correction was used to test the proportion of gender and background occupation of the subgroup and the full cohort. All statistical analyses and graphical displays were performed using R version 4.0.5 [27].

3. Results

3.1. Demographics of All Oil Spill Clean-Up Workers (n = 869)

During the clean-up, 2118 workers engaged in various activities. The demographic characteristics of these workers have been reported previously [9,17]. During the 5-year follow-up surveillance, 869 clean-up workers reported to Rayong Hospital at least once. The demographics of the 869 participants are summarised in Table 1.
The mean (±SD) age of the clean-up workers at baseline was 39.47 ± 10.23 years, and 86.30% of the study sample was male. Civilian volunteers and PTT global chemical (PTTGC) staff accounted for 29.80% and 38.80% of all clean-up workers, respectively. Based on the cotinine levels, 109 (12.50%) workers were active smokers (>50 ng/mL). Demographic characteristics and baseline concentrations of 1-OHPG and t,t-MA are summarised in Table 1.

3.2. Characteristics of Subgroup of Oil Spill Clean-Up Workers (n = 169)

A total of 169 oil spill clean-up workers provided post-shift urine samples and reported at least one follow-up visit. The demographic characteristics of the 169 workers are summarised in Table 2. The mean (±SD) age was 39.72 ± 9.91 years and 94.70% were men. Civilian volunteers accounted for 44.40% of the 169 workers. At baseline post-shift, 17 (10.10%) clean-up workers had high levels of 1-OHPG (>5.0 pmol/mol), 55 (32.50%) had moderate levels (1.00–5.00 pmol/mL), and 97 (57.40%) had low levels (<1.0 pmol/mol). Overall, 48 (28.40%) workers had detectable urinary t,t-MA levels. The demographic characteristics and baseline concentrations of 1-OHPG and t,t-MA in the 169 workers are summarised in Table 2.
Comparing between the subgroup and full cohort of 869 workers, the mean ages were not significantly different (Supplementary Table S1). However, the characteristics of the subgroup and the whole cohort were different in terms of gender and background occupation (Supplementary Table S1).

3.3. Latent Trajectories and Trends of Haematological, Hepatic, and Renal Parameters (n = 869)

The LCMMs were used to classify the trends of haematological, hepatic, and renal parameters of all 869 participants. The analyses showed that Rayong oil spill clean-up workers had between one and four latent trends for each index.
The best LCMM (based on BIC and posterior probabilities) for haemoglobin had three latent groups (high-, low-, and normal-stable). There were 33 clean-up workers (3.80%) with low haemoglobin levels at baseline and a stable trend over time (low-stable). Four workers (0.46%) had high haemoglobin levels at baseline, but with a significantly decreasing trend of 1.50 g/dL per year. At the baseline, the low-stable latent group had observed and predicted mean haemoglobin levels of 10.9 g/dL and 11.1 g/dL, respectively. The observed and predicted mean haemoglobin levels for the low-stable class were relatively stable over the period of the study, with no significant differences over time. The observed mean haemoglobin level for the high-decreasing group ranged between 16.8 g/dL at the baseline and 9.22 g/dL in 2018, while the predicted level ranged from 16.9 g/dL to 9.74 g/dL. The clean-up workers had three latent groups based on the levels and trends of haematocrit. Average haematocrit percentage was high among three (0.35%) clean-up workers, with a decreasing trend (4.84% per year). The remaining 866 (99.65) workers had normal haematocrit percentage at baseline, of which 6.10% had relatively lower average haematocrit percentage (indicated with asterisks), compared to 90.60%. Among the workers with high decreasing trend, the observed mean haematocrit percent was 46.9% at baseline, while the mean percent at the fifth year was 17%. The observed and predicted mean haematocrit levels for most of the workers were almost the same for every year. The WBC count of 2.42% of the clean-up workers demonstrated a high–decreasing trend (−0.73 × 103 cells/µL per year). Among these workers, the observed mean WBC count reduced from 12.79 × 103 cells/µL at baseline to 8.00 × 103 cells/µL in 2018. The mean WBC count for 97.5% of the workers increased slightly from 7.05 × 103 cells/µL to 7.11 × 103 cells/µL. All oil spill clean-up workers had normal RBC, platelet, and ANC levels at baseline. The majority of the workers (95.00%) had increasing trends of RBC from an average of 5.08 cells/µL to 5.12 cells/µL (trend = 0.01 cells/µL per year), and 8.75% had decreasing trends of ANC over time (−0.08 × 103 cells/µL per year) (Figure 1).
The average BUN level grouped the Rayong oil clean-up workers into two latent classes (normal–increasing and normal–stable), whereas three latent groups were observed for Cr levels (normal–decreasing, normal–stable, and normal-increasing). Although the average BUN and Cr levels were normal at baseline for all clean-up workers, 90.20% of the workers had a significantly increasing trend of BUN (0.31 mg/dL per year), and 85.70% had a significantly increasing trend of Cr (0.01 mg/dL per year). The observed mean BUN level for the normal–increasing group was 10.97 mg/dL at the baseline, and 12.65 mg/dL at the fifth year of the follow-up. Meanwhile, the Cr level increased from 0.98 mg/dL to 1.02 mg/dL among the normal increasing group. The latent trajectories and trends of haematological, hepatic, and renal parameters for all oil spill clean-up workers are shown in Figure 1, while the observed and predicted class-specific means for each haematological, hepatic, and renal index are shown in Figure 2.

3.4. Longitudinal Trajectories by Age Group

The trajectories of haemoglobin, haematocrit WBC, RBC, and platelets were similar for all age groups; however, the levels of these parameters at baseline were lower for the age group of 50 years and older. These parameters showed a slightly increasing trend over time, which did not differ by age.
The average BUN at baseline ranged between 10.5–12.5 mg/dL among all the workers. Average BUN and Cr levels increased throughout the study period. These trajectories were consistent for all age groups, except for those aged 50 years and above, whose BUN and Cr increased sharply throughout the study period.
AST and ALT levels showed stable trajectories among clean-up workers aged <50 years. Among workers older than 50 years, AST and ALT had an N-shaped trajectory, reaching peak levels three years after the clean-up work. The longitudinal trajectories of haematological, hepatic, and renal parameters according to age group are shown in Figure 3.
Although the RBC, WBC, and ANC show multiple slope trajectories, the variations between the baseline and the final follow-up were still within the normal range. These trajectories could be within normal biological variation.

3.5. Longitudinal Trajectories by Smoking Status (n = 169)

The longitudinal trajectories of haematological, hepatic, and renal parameters varied significantly according to smoking status. The average haemoglobin and haematocrit levels in non-smokers exhibited a decreasing trend over the study period, whereas the levels fluctuated in smokers.
The average WBC and ANC levels among non-smokers showed an N-shaped trajectory during the 5 years after the clean-up. Meanwhile, the WBC and ANC counts from smokers had a decreasing trajectory, particularly between the third and fourth years. The trajectories of BUN were stable during the 2 years and showed an increasing trend at 3–5 years after the clean-up for both smokers and non-smokers. The Cr, AST, and ALT levels of smokers and non-smokers fluctuated during 5 years of follow-up. The longitudinal trajectories of haematological, hepatic, and renal parameters by smoking status of the 169 oil spill clean-up workers are shown in Figure 4.

3.6. Longitudinal Trajectories by Concentrations of Urinary 1-OHPG and t,t-MA (n = 169)

Most haematological indices showed varying trajectories. Haemoglobin and haematocrit levels among workers with high urinary 1-OHPG concentrations exhibited a sigmoid pattern with decreasing levels over time compared with the other exposure groups. Workers with low urinary 1-OHPG concentrations had stable trajectories of haemoglobin and haematocrit levels over time. The trajectories of WBC, ANC, and platelets were almost linearly stable during the five years, except workers with high urinary 1-OHPG concentrations exhibited a relatively lower platelet count than those in the other exposure groups. Figure 5 depicts the trajectories of all the indices based on urinary 1-OHPG concentrations.
Urinary t,t-MA was not detectable in 71.6% of the workers. The longitudinal trajectories of workers with detectable and undetectable concentrations of t,t-MA differed for most indices. The haematocrit level for a detectable t,t-MA concentration showed a U-shaped trajectory and decreased gradually until the fourth follow-up. WBC and ANC counts were generally higher in the detectable group than in the undetectable group. The longitudinal trajectories of the platelet counts were similar for both t,t-MA groups. The trends of BUN and Cr levels did not visibly differ between the detectable and undetectable groups. The longitudinal trajectories of the haematological, hepatic, and renal indices according to the baseline t,t-MA concentrations are shown in Figure 6.

4. Discussion

This longitudinal cohort study identified different latent longitudinal trajectories and trends in the haematological, hepatic, and renal profiles of 869 Rayong oil spill clean-up workers. Linking exposure to PAHs and VOCs, subgroup analyses assessed the longitudinal trend of the profiles among 169 cleaners with available data for urinary concentrations of 1-OHPG and t,t-MA.
Generally, no significant changes in haemoglobin and haematocrit levels were found among most clean-up workers, except for an observable decreased haemoglobin (<10 g/dL) and haematocrit (<24%) among a small group of workers. These workers might develop anaemic symptoms, such as fatigue, syncope, or lower work tolerance [28]. Age and cigarette smoking at the time of the oil spill clean-up showed minimal to no effects on haemoglobin and haematocrit levels. Subgroup analysis indicated that high 1-OHPG and detectable t,t-MA levels might be associated with reduced haemoglobin and haematocrit levels at baseline and with decreasing trajectories over time. Acutely elevated haemoglobin has been reported among a Canadian sub-population exposed to VOCs [28], which is consistent with the high haemoglobin level observed among a fraction of Rayong oil spill clean-up workers. Some studies conducted 2–7 years after oil spill exposure reported no significant changes in haemoglobin levels [5,6]. Studies have found that PAHs and VOCs can induce haemolytic anaemia [29,30], which may result in reduced haemoglobin and haematocrit levels. Another study by Kamal et al., also reported that exposure to PAHs is associated with reduced haematocrit levels in humans [15]. Consistently, our subgroup analysis revealed that urinary PAH and VOC metabolite levels were associated with reduced haemoglobin and haematocrit levels. Although chronic smoking is known to increase haemoglobin and haematocrit levels in human and animal models in response to hypoxaemia [31,32], this study cohort of smokers exhibited decreased or stable haemoglobin and haematocrit levels during the follow-up period. These findings suggest that the longitudinal changes in haemoglobin and haematocrit among clean-up workers are likely due to exposure to PAH and VOCs from other sources, including crude oil, rather than cigarette smoking.
Two remarkable trends in WBC, RBC, and platelet counts should be noted in workers. These haematological indices exhibited reduced levels and decreasing trends over time in a small group of workers, while some workers exhibited increasing trajectories in these indices. From the trajectory analysis, both trends were not associated with age and were only partially associated with smoking. Workers with high and detectable concentrations of 1-OHPG and t,t-MA at baseline had higher trajectories of WBC, RBC, and platelet counts than those with low or undetectable concentrations.
Increased WBC, RBC, and platelet counts due to oil spill exposure have been reported in previous studies [6,33]. A population-based study reported that high concentrations of VOC and PAH metabolites were associated with elevated levels of haematological indices, including WBC and RBC [28]. Similarly, the concentrations of 1-OHPG and t,t-MA were associated with high WBC and RBC counts in Rayong oil spill workers. Following exposure to oil spills and xenobiotic metabolism, benzene-active metabolites, including tt-muconaldehyde, 1,4-benzoquinone, hydroquinone, and catechol, can suppress the bone marrow [34,35,36]. Such haematopoietic damage can manifest in many ways, including alterations in the haematological profiles. Therefore, the alterations in haematological functions identified among the Rayong oil spill clean-up workers may be attributable to the suppression of bone marrow function by the effects of PAHs and VOCs in crude oil.
All oil-spill clean-up workers had normal BUN and Cr levels at baseline, but most had corresponding increasing trends. The longitudinal trajectories of the kidney indices were age-related, as higher BUN and Cr levels were observed among workers in the older age group. Higher and moderate 1-OHPG and detectable t,t-MA concentrations correlated with increasing BUN levels. Despite the observable association between BUN, Cr, and age [37], the possible effects of oil spill exposure on renal indices cannot be ruled out [5,6]. PAH metabolites, such as 1-OHPG, can cause oxidative stress and induce kidney damage [38,39]. AST and ALT levels are surrogates for hepatocellular injury and are highly sensitive to PAH and VOC exposure, as reported in various studies [40,41]. Although studies have reported possible liver damage from exposure to PAHs and VOCs, longitudinal trends in our research indicated no significant changes.

Strengths and Limitations of the Study

This study found possible alterations in the haematopoietic system of oil spill clean-up workers. The findings of this study should be interpreted with consideration of its limitations. Most clean-up workers who attended follow-up visits did not provide urine samples at the baseline. The information on the concentrations of 1-OHPG and t,t-MA to assess exposure levels was only available for a relatively small group of workers. Nonetheless, the characteristics of this subgroup are similar to those of the entire cohort. Therefore, our results may be generalisable to the entire cohort. Additionally, selection (volunteering) bias was possible because the follow-up visit was voluntary. The characteristics of the workers who provided urine samples were different from that of the full cohort of 869 workers. Therefore, results from analysis based on urine biomarker concentration should be interpreted cautiously. Although most of the significant trends are shown in the small groups, with relatively small sample size, the trends are clinically significant. For instance, for haematocrit, the level decreases until possible anaemic symptoms could occur. Despite these limitations, this study provides insight into the possible long-term effects of oil spill exposure.
Previous studies involving subjects exposed to oil spill exposure have reported general findings. However, the latent class trajectory analysis used in this study identified small unique groups of oil spill clean-up workers with significant changes in haematological, hepatic, or renal profiles. Our findings suggest possible long-term haematopoietic function alterations from oil spill exposure, indicated by non-linear trajectories and significant long-term changes in haematological, hepatic, and renal profiles in workers after the oil spill. However, to assess whether these results were incidental or due to oil spill exposure, more studies including prospective cohort studies with a rigorous methodologic design and with periodic biomarker measurements for the oil spill clean-up workers are needed.

5. Conclusions

More than 5 years after exposure to PAHs and VOCs through oil spill clean-ups, alterations in the haematological, hepatic, and renal systems were still observable among the workers. Therefore, crude oil-related activities such as oil spill clean-up should be classified as potentially dangerous jobs exposing workers to carcinogens and haematotoxins. Personal protective equipment should be used in crude oil-related activities such as oil spill clean-up to reduce the risks of the exposure to carcinogens and haematotoxins. Health monitoring including blood sampling for haematological, hepatic, and renal indices should be planned before the workers participate in oil spill clean-ups to include assessments throughout the pre-placement of jobs and immediately post-shift, and then periodically for years following the clean-up. Future studies should assess the changes of biological indices adding oxidative stress biomarkers such as malondialdehyde and epigenetic changes, starting from pre-placement examination to further study the possibility of the acute and long-term effects from the oil spill.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph20042871/s1, Table S1: Comparison of demographic between the cohort (N = 869) and the subgroup (N = 169).

Author Contributions

Conceptualization, B.A.O., A.L., C.I., S.R. and T.I.; Data curation, C.I., and S.S.; Formal analysis, B.A.O., N.P. and T.I.; Investigation, B.A.O., N.P. and T.I.; Methodology, B.A.O., N.P. and T.I.; Project administration, A.L. and T.I.; Resources, C.I., S.R., S.S. and T.I.; Supervision, A.L. and T.I.; Validation, B.A.O., A.L. and T.I.; Visualization, B.A.O. and N.P.; Writing—original draft, B.A.O., N.P. and T.I.; Writing—review & editing, B.A.O., A.L., N.P., C.I., S.R., S.S. and T.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. However, B.A.O. received partial tuition fee funding from the Faculty of Science and Technology, Prince of Songkla University Pattani Campus and Air Pollution and Health Effect Research Center, Prince of Songkla University, Hat Yai Campus, Thailand.

Institutional Review Board Statement

The Research Ethics Committee of the Prince of Songkla University approved this study (Institutional Review Board [IRB] approval number psu.pn.1-001/64).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data, codes, and materials generated and analysed in this published article are available from the corresponding author upon reasonable request.

Acknowledgments

The authors are grateful to the Multidisciplinary Research and Innovation Centre, Kumasi, Ghana, for assisting with the data analysis. We thank Rayong Hospital for providing the oil spill data. Special gratitude is given to Don McNeil for his comments and valuable suggestions throughout this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Beyer, J.; Jonsson, G.; Porte, C.; Krahn, M.M.; Ariese, F. Analytical Methods for Determining Metabolites of Polycyclic Aromatic Hydrocarbon (PAH) Pollutants in Fish Bile: A Review. Environ. Toxicol. Pharmacol. 2010, 30, 224–244. [Google Scholar] [CrossRef] [PubMed]
  2. De Hoop, L.; Schipper, A.M.; Leuven, R.S.E.W.; Huijbregts, M.A.J.; Olsen, G.H.; Smit, M.G.D.; Hendriks, A.J. Sensitivity of Polar and Temperate Marine Organisms to Oil Components. Environ. Sci. Technol. 2011, 45, 9017–9023. [Google Scholar] [CrossRef] [PubMed]
  3. Baan, R.; Grosse, Y.; Straif, K.; Secretan, B.; El Ghissassi, F.; Bouvard, V.; Benbrahim-Tallaa, L.; Guha, N.; Freeman, C.; Galichet, L.; et al. A Review of Human Carcinogens--Part F: Chemical Agents and Related Occupations. Lancet Oncol. 2009, 10, 1143–1144. [Google Scholar] [CrossRef] [PubMed]
  4. International Agency for Research on Cancer. Occupational Exposures in Petroleum Refining; Crude on and Major Petroleum Fuels. Int. Agency Res. Cancer 1988, 45, 331. [Google Scholar]
  5. Choi, Y.H.; Hong, J.Y.; Lee, M.S. A Retrospective Mid- and Long-Term Follow-up Study on the Changes in Hematologic Parameters in the Highly Exposed Residents of the Hebei Spirit Oil Spill in Taean, South Korea. Osong Public Health Res. Perspect. 2017, 8, 358–366. [Google Scholar] [CrossRef]
  6. D’Andrea, M.A.; Reddy, G.K. The Development of Long-Term Adverse Health Effects in Oil Spill Cleanup Workers of the Deepwater Horizon Offshore Drilling Rig Disaster. Front. Public Health 2018, 6, 117. [Google Scholar] [CrossRef]
  7. Jongeneelen, F.J. Benchmark Guideline for Urinary 1-Hydroxypyrene as Biomarker of Occupational Exposure to Polycyclic Aromatic Hydrocarbons. Ann. Occup. Hyg. 2001, 45, 3–13. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Ding, J.; Shen, G.; Zhong, J.; Wang, C.; Wei, S.; Chen, C.; Chen, Y.; Lu, Y.; Shen, H.; et al. Dietary and Inhalation Exposure to Polycyclic Aromatic Hydrocarbons and Urinary Excretion of Monohydroxy Metabolites—A Controlled Case Study in Beijing, China. Environ. Pollut. 2014, 184, 515–522. [Google Scholar] [CrossRef]
  9. Sithisarankul, P.; Intawong, C. Preliminary Report of Health Effects Among Oil Spill Cleanup Workers and Volunteers, Thailand, 2013. J. Health Res. 2015, 29, 197–201. [Google Scholar]
  10. Wiwanitkit, V.; Suwansaksri, J.; Nasuan, P. Urine Trans,Trans-Muconic Acid as a Biomarker for Benzene Exposure in Gas Station Attendants in Bangkok, Thailand. Ann. Clin. Lab. Sci. 2001, 31, 399–401. [Google Scholar]
  11. Bechtold, W.E.; Lucier, G.; Birnbaum, L.S.; Yin, S.N.; Li, G.L.; Henderson, R.F. Muconic Acid Determinations in Urine as a Biological Exposure Index for Workers Occupationally Exposed to Benzene. Am. Ind. Hyg. Assoc. J. 1991, 52, 473–478. [Google Scholar] [CrossRef] [PubMed]
  12. Jain, R.B. Levels of Selected Urinary Metabolites of Volatile Organic Compounds among Children Aged 6-11 Years. Environ. Res. 2015, 142, 461–470. [Google Scholar] [CrossRef] [PubMed]
  13. Strickland, P.; Kang, D.; Sithisarankul, P. Polycyclic Aromatic Hydrocarbon Metabolites in Urine as Biomarkers of Exposure and Effect. Environ. Health Perspect. 1996, 104, 927–932. [Google Scholar] [CrossRef] [PubMed]
  14. Abou-ElWafa, H.S.; Albadry, A.A.; El-Gilany, A.H.; Bazeed, F.B. Some Biochemical and Hematological Parameters among Petrol Station Attendants: A Comparative Study. BioMed Res. Int. 2015, 2015, 418724. [Google Scholar] [CrossRef] [PubMed]
  15. Kamal, A.; Cincinelli, A.; Martellini, T.; Malik, R.N. Biomarkers of PAH Exposure and Hematologic Effects in Subjects Exposed to Combustion Emission during Residential (and Professional) Cooking Practices in Pakistan. Environ. Sci. Pollut. Res. 2016, 23, 1284–1299. [Google Scholar] [CrossRef]
  16. Samadi, M.T.; Shakerkhatibi, M.; Poorolajal, J.; Rahmani, A.; Rafieemehr, H.; Hesam, M. Association of Long Term Exposure to Outdoor Volatile Organic Compounds (BTXS)with pro-Inflammatory Biomarkers and Hematologic Parameters in Urban Adults: A Cross-Sectional Study in Tabriz, Iran. Ecotoxicol. Environ. Saf. 2019, 180, 152–159. [Google Scholar] [CrossRef]
  17. Ingviya, T.; Intawong, C.; Abubaker, S.; Strickland, P.T. Exposure Assessment of Rayong Oil Spill Cleanup Workers. Expo. Health 2020, 12, 617–628. [Google Scholar] [CrossRef]
  18. Owusu, B.A.; Lim, A.; Intawong, C.; Rheanpumikankit, S.; Suksri, S.; Ingviya, T. Haematological, Renal, and Hepatic Function Changes among Rayong Oil Spill Clean-up Workers: A Longitudinal Study. Int. Arch. Occup. Environ. Health 2022, 95, 1481–1489. [Google Scholar] [CrossRef]
  19. Zielinska-Danch, W.; Wardas, W.; Sobczak, A.; Szołtysek-Bołdys, I. Estimation of Urinary Cotinine Cut-off Points Distinguishing Non-Smokers, Passive and Active Smokers. Biomarkers 2007, 12, 484–496. [Google Scholar] [CrossRef]
  20. Kang, D.; Rothman, N.; Cho, S.H.; Lim, H.S.; Kwon, H.J.; Kim, S.M.; Schwartz, B.; Strickland, P.T. Association of Exposure to Polycyclic Aromatic Hydrocarbons (Estimated from Job Category) with Concentration of 1-Hydroxypyrene Glucuronide in Urine from Workers at a Steel Plant. Occup. Environ. Med. 1995, 52, 593–599. [Google Scholar] [CrossRef]
  21. Intawong, C.; Khunkitti, R.; Sripaung, N.; Yapun, C.; Buachoom, P.; Tupwongse, V.; Sithisarankul, P. Interlaboratory Comparison for Urinary Trans, Trans-Muconic Acid Testing in Rayong Province. Thammasat Med. J. 2015, 15, 363–375. [Google Scholar]
  22. Dayimu, A.; Qian, W.; Fan, B.; Wang, C.; Li, J.; Wang, S.; Ji, X.; Zhou, G.; Zhang, T.; Xue, F. Trajectories of Haemoglobin and Incident Stroke Risk: A Longitudinal Cohort Study. BMC Public Health 2019, 19, 1395. [Google Scholar] [CrossRef] [PubMed]
  23. Farooqui, W.A.; Uddin, M.; Qadeer, R.; Shafique, K. Trajectories of Vital Status Parameters and Risk of Mortality among Acute Organophosphorus Poisoning Patients-A Latent Class Growth Analysis. BMC Public Health 2020, 20, 1538. [Google Scholar] [CrossRef] [PubMed]
  24. Randall, D.A.; Patterson, J.A.; Gallimore, F.; Morris, J.M.; Simpson, J.M.; McGee, T.M.; Ford, J.B. Haemoglobin Trajectories during Pregnancy and Associated Outcomes Using Pooled Maternity and Hospitalization Data from Two Tertiary Hospitals. Vox Sang. 2019, 114, 842–852. [Google Scholar] [CrossRef] [PubMed]
  25. Proust-Lima, C.; Philipps, V.; Liquet, B. Estimation of Extended Mixed Models Using Latent Classes and Latent Processes: The R Package Lcmm. J. Stat. Softw. 2017, 78. [Google Scholar] [CrossRef]
  26. Kasper, D.; Fauci, A.; Hauser, S.; Longo, D.; Jameson, L.J.; Loscalzo, J. Harrison’s Principles of Internal Medicine Vol 1, 20th ed.; Mcgraw-Hill: New York, NY, USA, 2015. [Google Scholar]
  27. R Core Team. The R Project for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  28. Cakmak, S.; Cole, C.; Hebbern, C.; Andrade, J.; Dales, R. Associations between Blood Volatile Organic Compounds, and Changes in Hematologic and Biochemical Profiles, in a Population-Based Study. Environ. Int. 2020, 145, 106121. [Google Scholar] [CrossRef]
  29. Kponee, K.Z.; Chiger, A.; Kakulu, I.I.; Vorhees, D.; Heiger-Bernays, W. Petroleum Contaminated Water and Health Symptoms: A Cross-Sectional Pilot Study in a Rural Nigerian Community. Environ. Health Glob. Access Sci. Source 2015, 14, 86. [Google Scholar] [CrossRef]
  30. McLoone, P.; Dyussupov, O.; Nurtlessov, Z.; Kenessariyev, U.; Kenessary, D. The Effect of Exposure to Crude Oil on the Immune System. Health Implications for People Living near Oil Exploration Activities. Int. J. Environ. Health Res. 2021, 31, 762–787. [Google Scholar] [CrossRef]
  31. Malenica, M.; Prnjavorac, B.; Bego, T.; Dujic, T.; Semiz, S.; Skrbo, S.; Gusic, A.; Hadzic, A.; Causevic, A. Effect of Cigarette Smoking on Haematological Parameters in Healthy Population. Med. Arch. 2017, 71, 132–136. [Google Scholar] [CrossRef]
  32. Pedersen, K.M.; Çolak, Y.; Ellervik, C.; Hasselbalch, H.C.; Bojesen, S.E.; Nordestgaard, B.G. Smoking and Increased White and Red Blood Cells: A Mendelian Randomization Approach in the Copenhagen General Population Study. Arterioscler. Thromb. Vasc. Biol. 2019, 39, 965–977. [Google Scholar] [CrossRef]
  33. Ramírez, M.I.; Arévalo-Jaramillo, A.P.; Espinosa, C.I.; Bailon-Moscoso, N. Is the Anemia in Men an Effect of the Risk of Crude Oil Contamination? Toxicol. Rep. 2022, 9, 480–486. [Google Scholar] [CrossRef] [PubMed]
  34. Shahsavani, S.; Fararouei, M.; Soveid, M.; Hoseini, M.; Dehghani, M. The Association between the Urinary Biomarkers of Polycyclic Aromatic Hydrocarbons and Risk of Metabolic Syndromes and Blood Cell Levels in Adults in a Middle Eastern Area. J. Environ. Health Sci. Eng. 2021, 19, 1667–1680. [Google Scholar] [CrossRef] [PubMed]
  35. Ross, D. The Role of Metabolism and Specific Metabolites in Benzene-Induced Toxicity: Evidence and Issues. J. Toxicol. Environ. Health Part A 2000, 61, 357–372. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, L.; He, X.; Bi, Y.; Ma, Q. Stem Cell and Benzene-Induced Malignancy and Hematotoxicity. Chem. Res. Toxicol. 2012, 25, 1303–1315. [Google Scholar] [CrossRef] [PubMed]
  37. Poli, D.; Mozzoni, P.; Pinelli, S.; Cavallo, D.; Papaleo, B.; Caporossi, L. Sex Difference and Benzene Exposure: Does It Matter? Int. J. Environ. Res. Public Health 2022, 19, 2339. [Google Scholar] [CrossRef]
  38. Palackal, N.T.; Lee, S.H.; Harvey, R.G.; Blair, I.A.; Penning, T.M. Activation of Polycyclic Aromatic Hydrocarbon Trans-Dihydrodiol Proximate Carcinogens by Human Aldo-Keto Reductase (AKR1C) Enzymes and Their Functional Overexpression in Human Lung Carcinoma (A549) Cells. J. Biol. Chem. 2002, 277, 24799–24808. [Google Scholar] [CrossRef]
  39. Sun, S.; Mao, W.; Tao, S.; Zou, X.; Tian, S.; Qian, S.; Yao, C.; Zhang, G.; Chen, M. Polycyclic Aromatic Hydrocarbons and the Risk of Kidney Stones in US Adults: An Exposure-Response Analysis of NHANES 2007–2012. Int. J. Gen. Med. 2021, 14, 2665–2676. [Google Scholar] [CrossRef]
  40. Zhu, L.; Fang, Z.; Jin, Y.; Chang, W.; Huang, M.; He, L.; Chen, Y.; Yao, Y. Association between Serum Alanine and Aspartate Aminotransferase and Blood Pressure: A Cross-Sectional Study of Chinese Freshmen. BMC Cardiovasc. Disord. 2021, 21, 472. [Google Scholar] [CrossRef]
  41. Wu, M.-T.; Mao, I.-F.; Wypij, D.; Ho, C.-K.; Chen, J.-R.; Christiani, D.C. Serum Liver Function Profiles in Coking Workers. Am. J. Ind. Med. 1997, 32, 478–486. [Google Scholar] [CrossRef]
Figure 1. Latent trends of the haematological and hepatic parameters. (A) Haemoblogin: Only 3 significant latent groups: Low-stable, Normal-stable, high-decreasing. (B) Haematocrit: Only 3 significant latent groups: Normal-low-stable, Normal-high-stable, High-decreasing. (C) White blood cell: Only 2 significant latent groups: Normal-increasing, Normal-decreasing. (D) Red blood cell: Only 3 significant latent groups: Normal-decreasing, Normal-stable, Normal-increasing. (E) Platelet counts: Only 2 significant latent groups: Normal-increasing slow, Normal-increasing sharp. (F) Absolute Neutrophil count: Two significant latent groups: Normal-stable, Normal-decreasing. (G) Blood urea nitrogen: Only 2 significant latent groups: Normal-stable, Normal-increasing. (H) Serum creatinine: Only 3 significant latent groups: Normal-decreasing, Normal-stable, Normal-increasing. (I) Aspartate aminotransferase: Only 1 latent group: Normal-decreasing. (J) aminotransferase: Only 2 significant latent groups: Normal-decreasing, Normal-stable Note: (*) represents significant trends.
Figure 1. Latent trends of the haematological and hepatic parameters. (A) Haemoblogin: Only 3 significant latent groups: Low-stable, Normal-stable, high-decreasing. (B) Haematocrit: Only 3 significant latent groups: Normal-low-stable, Normal-high-stable, High-decreasing. (C) White blood cell: Only 2 significant latent groups: Normal-increasing, Normal-decreasing. (D) Red blood cell: Only 3 significant latent groups: Normal-decreasing, Normal-stable, Normal-increasing. (E) Platelet counts: Only 2 significant latent groups: Normal-increasing slow, Normal-increasing sharp. (F) Absolute Neutrophil count: Two significant latent groups: Normal-stable, Normal-decreasing. (G) Blood urea nitrogen: Only 2 significant latent groups: Normal-stable, Normal-increasing. (H) Serum creatinine: Only 3 significant latent groups: Normal-decreasing, Normal-stable, Normal-increasing. (I) Aspartate aminotransferase: Only 1 latent group: Normal-decreasing. (J) aminotransferase: Only 2 significant latent groups: Normal-decreasing, Normal-stable Note: (*) represents significant trends.
Ijerph 20 02871 g001
Figure 2. Class-specific observed and predicted means of haematological, hepatic, and renal profiles.
Figure 2. Class-specific observed and predicted means of haematological, hepatic, and renal profiles.
Ijerph 20 02871 g002
Figure 3. Longitudinal trajectories by age group at baseline.
Figure 3. Longitudinal trajectories by age group at baseline.
Ijerph 20 02871 g003
Figure 4. Longitudinal trajectories by smoking status.
Figure 4. Longitudinal trajectories by smoking status.
Ijerph 20 02871 g004
Figure 5. Longitudinal trajectories of the haematological, hepatic, and renal indices by the baseline concentration of 1-OHPG.
Figure 5. Longitudinal trajectories of the haematological, hepatic, and renal indices by the baseline concentration of 1-OHPG.
Ijerph 20 02871 g005
Figure 6. Longitudinal trajectories of the haematological, hepatic, and renal indices by the baseline concentration of t,t-MA.
Figure 6. Longitudinal trajectories of the haematological, hepatic, and renal indices by the baseline concentration of t,t-MA.
Ijerph 20 02871 g006
Table 1. Description of the participants’ demographics.
Table 1. Description of the participants’ demographics.
Demographic FactorDescriptionNumber of WorkersPercentage
SexMen75086.30%
Women11913.70%
Age group at baseline20–2914016.10%
30–3916719.20%
40–4927631.80%
50+9310.70%
Not stated19322.20%
Background occupation* PTTGC staff33738.80%
Civilian25929.80%
Military27331.40%
Smoking status
Non-smokersCotinine < 50 ng/mL606.90%
SmokersCotinine > 50 ng/mL10912.50%
UnknownCotinine levels are not available.70080.60%
* PTT global chemical (PTTGC).
Table 2. Characteristics of the subgroup with available urine samples.
Table 2. Characteristics of the subgroup with available urine samples.
Demographic FactorsDescriptionNumber of WorkersPercentage
SexMen16094.70%
Women95.30%
Age group at baseline (years)20–292213.00%
30–394023.90%
40–497343.10%
50+3420.00%
Background occupationPTTGC staff2917.20%
Civilian7544.40%
Military6538.40%
Cotinine level (ng/mL)Median (1st–3rd quartiles)7.51 (3.11–1040.79)
Non-smokersCotinine < 50 ng/mL10964.50%
SmokersCotinine > 50 ng/mL6035.50%
* 1-OHPG (n = 169)Median (1st–3rd quartiles)0.76 (0.31–2.27)
High1-OHPG > 5.0 pmol/mL1710.10%
Low1-OHPG < 1.0 pmol/mL9757.40%
Moderate1.0–5.0 pmol/mL5532.50%
** t,t-MA (n = 169)Median (1st–3rd quartiles)0. 00 (0.00–36.40)
Detectable4828.40%
Undetectable12171.60%
* 1-hydroxypyrene-glucuronide (1-OHPG); ** trans,trans-muconic acid (t,t-MA).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Owusu, B.A.; Lim, A.; Pongsiri, N.; Intawong, C.; Rheanpumikankit, S.; Suksri, S.; Ingviya, T. Latent Trajectories of Haematological, Hepatic, and Renal Profiles after Oil Spill Exposure: A Longitudinal Analysis. Int. J. Environ. Res. Public Health 2023, 20, 2871. https://doi.org/10.3390/ijerph20042871

AMA Style

Owusu BA, Lim A, Pongsiri N, Intawong C, Rheanpumikankit S, Suksri S, Ingviya T. Latent Trajectories of Haematological, Hepatic, and Renal Profiles after Oil Spill Exposure: A Longitudinal Analysis. International Journal of Environmental Research and Public Health. 2023; 20(4):2871. https://doi.org/10.3390/ijerph20042871

Chicago/Turabian Style

Owusu, Benjamin Atta, Apiradee Lim, Nitinun Pongsiri, Chanthip Intawong, Sunthorn Rheanpumikankit, Saijit Suksri, and Thammasin Ingviya. 2023. "Latent Trajectories of Haematological, Hepatic, and Renal Profiles after Oil Spill Exposure: A Longitudinal Analysis" International Journal of Environmental Research and Public Health 20, no. 4: 2871. https://doi.org/10.3390/ijerph20042871

APA Style

Owusu, B. A., Lim, A., Pongsiri, N., Intawong, C., Rheanpumikankit, S., Suksri, S., & Ingviya, T. (2023). Latent Trajectories of Haematological, Hepatic, and Renal Profiles after Oil Spill Exposure: A Longitudinal Analysis. International Journal of Environmental Research and Public Health, 20(4), 2871. https://doi.org/10.3390/ijerph20042871

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop