Weather Indices. The average temperature was 28.7°C and average humidity was 80.4% across the two study days in November with the maximum temperature reaching 31.8°C (Table 1). In April, the average temperature was 29.0°C and average humidity was 76.5% across the two study collection days, with the maximum temperature reaching 34.0°C.
Table 1
Summarized statistics of worker characteristics, biomarkers, and weather indices.
| November (baseline), n = 40 | April (end of season), n = 40 | P-value |
Demographics and Clinical | | | |
Age, years, mean ± SD | 33 ± 11 | - | - |
HbA1c, %, mean ± SD | 5.6 ± 0.4 | Not measured | - |
SBP, mmHg, mean ± SD | 107 ± 9 | 116 ± 12 | < 0.01* |
DBP, mmHg, mean ± SD | 72 ± 8 | 68 ± 10 | 0.01* |
- HTN, n (%) | 0 (0%) | 3 (8%) | - |
Tons of cane cut, median (IQR) | 5.6 (5.0, 6.2) | 5.4 (4.0, 6.2) | 0.25 |
Biomarkers, median (IQR) | | | |
Serum Uric Acid, mg/dL | 5.84 (4.85, 6.57) | 5.40 (4.80, 6.29) | 0.20 |
Serum Creatine Kinase, units/L | 260.0 (217.0, 379.5) | 321.0 (248.0, 463.0) | 0.03* |
Serum Creatinine, µmol/L | 77.81 (69.85, 93.73) | 76.93 (68.97, 91.96) | 0.64 |
Serum Cystatin C, mg/L | 0.75 (0.53, 0.94) | 0.59 (0.48, 0.79) | 0.03* |
eGFR-Creatinine, ml/min/1.73m2 | 110.62 (88.82, 121.58) | 108.41 (97.35,118.08) | 0.77 |
- < 90, n (%) | 10 (25%) | 6 (15%) | < 0.01* |
eGFR-Cystatin C, ml/min/1.73m2 | 118.13 (92.35, 146.21) | 136.43 (113.2, 150.41) | 0.03* |
- < 90, n (%) | 9 (23%) | 5 (13%) | 0.25 |
Weather Indices, 2-day average [day 1, day 2] | | | |
Mean Humidity, % | 80.4 [84.2, 76.6] | 76.5 [81.9, 71.1] | - |
Maximum Humidity, % | 100.0 [100.0, 100.0] | 100.0 [100.0, 100.0] | - |
Mean Temperature, °C | 28.7 [28.2, 29.1] | 29.0 [28.0, 30.0] | - |
Maximum Temperature, °C | 31.6 [31.4, 31.8] | 32.9 [31.7, 34.0] | - |
SD, standard deviation. IQR, interquartile range. SBP, systolic blood pressure. DBP, diastolic blood pressure. HTN, Hypertension was defined as systolic blood pressure ≥ 140 mmHg or diastolic blood pressure ≥ 90 mmHg. eGFR, estimated glomerular filtration rate. P-values were obtained by paired Wilcoxon signed-rank test. |
* Indicates significance at p < 0.05. |
Clinical and Biomarker Data. We included 40 participants in this study with data collected at the start and the end of the harvest, as presented in Table 1. The mean age was 33 years. We observed more hypertension in April as compared to November based on higher systolic blood pressure readings. Work intensity, measured by tons cut, was slightly higher in November compared to April but was not significantly different. We observed significantly higher CK levels in April compared to November.
{Insert Table 1}
HSP70 and anti-HSP70. At the start of the harvest in November (baseline), the median serum level of eHSP70 was 2.24 (interquartile range (IQR) 1.7, 3.16) ng/mL and, in April, was 3.33 (IQR 2.31, 4.32) ng/mL (Table 2). For serum anti-HSP70, the median level was 86.90 (IQR 47.34, 120.65) µg/mL in November and 98.31 (IQR 60.93-145.13) µg/mL in April. On average, serum eHSP70 levels in April were 1.1 ng/ml higher than baseline levels in November (p < 0.01). Similarly, the serum anti-HSP70 levels were 30.6 µg/mL higher in April as compared to November (p < 0.01). Distributions of HSP70 and anti-HSP70 and log transformed HSP70 and anti-HSP70 are shown in Figs. 1 and 2.
Table 2
Summarized statistics for extracellular heat shock protein 70 (eHSP70) and antibodies to HSP70 (anti-HSP70) levels.
| November (baseline), n = 40 | April (end of season), n = 40 | P-value | Change from November to April* |
eHSP70, ng/mL, median (IQR) | 2.24 (1.71, 3.16) | 3.33 (2.31, 4.32) | 0.01 | Mean ± SD: 1.14 ± 2.70 Median (IQR): 0.98 (-0.11, 2.29) |
Anti-HSP70, µg/mL, median (IQR) | 86.90 (47.34, 120.65) | 98.31 (60.93, 145.13) | 0.01 | Mean ± SD: 30.56 ± 67.96 Median (IQR): 16.04 (-12.33, 53.83) |
IQR, interquartile range. SD, standard deviation. P-values between November and April values were obtained by paired Wilcoxon signed-rank test. |
*Positive change represents higher values in April. |
{Insert FIGURES 1 and 2}
November serum eHSP70 levels were moderately positively correlated with April eHSP70 levels (r: 0.38; p = 0.02) and November serum anti-HSP70 levels were strongly positively correlated with April anti-HSP70 levels (r: 0.71; p < 0.01). We observed that cross-harvest change in eHSP70 was positively correlated with cross-harvest change in anti-HSP70 (r = 0.40; p = 0.01).
{Insert Table 2}
-
Heat Exposure and HSP70 and anti-HSP70 levels
We conducted a repeated-measure analysis to investigate the acute associations between log-transformed serum eHSP70 and anti-HSP70 and meteorological conditions. In the univariate models, increasing mean and maximum temperatures were both associated with increasing serum eHSP70 (Table 3). Decreasing average humidity was associated with increasing serum eHSP70. In the multivariable model, maximum temperature remained associated with serum eHSP70 (β: 0.21, 95% CI: 0.09, 0.33) after controlling for age, systolic and diastolic blood pressure (Table 4). Decreasing uric acid levels were associated with increasing serum eHSP70 in the univariate model but this relationship was no longer significant in the multivariable model (β: -0.09, 95% CI -0.21, 0.03). We included maximum temperature as the one meteorological variable in the multivariable model, due to multicollinearity issues and consistency with the serum anti-HSP70 model (see below).
Table 3
Results from univariate linear mixed-effect models examining associations between serum extracellular Heat Shock Protein 70 (eHSP70), and antibodies to HSP70, (anti-HSP70) and weather indices and markers of kidney function.
| eHSP70 Models | | Anti-HSP70 Models | |
| EstimateA (95% CI) | p-value | EstimateA (95% CI) | p-value |
Clinical | | | | |
Age | 0.01 (-0.01, 0.02) | 0.47 | 0.001 (-0.02, 0.03) | 0.92 |
HbA1c | 0.24 (-0.16, 0.64) | 0.24 | -0.628 (-1.27, 0.01) | 0.05 |
Systolic | 0.01 (-0.004, 0.02) | 0.20 | 0.015 (0.01, 0.03) | < 0.01* |
Diastolic | -0.01 (-0.02, 0.01) | 0.47 | 0.0001 (-0.01, 0.01) | 0.99 |
Biomarkers | | | | |
Uric Acid | -0.12 (-0.23, -0.001) | 0.048* | -0.138 (-0.28, 0.01) | 0.07 |
Creatine Kinase | 0.0003 (-0.0003, 0.001) | 0.36 | 0.00038 (-0.0002, 0.001) | 0.21 |
eGFR-Creatinine | -0.001 (-0.01, 0.01) | 0.83 | 0.001 (-0.01, 0.01) | 0.88 |
eGFR-Cystatin C | 0.002 (-0.002, 0.01) | 0.28 | -0.002 (-0.01, 0.0) | 0.20 |
eHSP70 | - | | 0.085 (0.03, 0.14) | < 0.01* |
Weather Indices | | | | |
Mean Humidity | -0.05 (-0.08, -0.02) | < 0.01* | -0.032 (-0.06, 0) | 0.04* |
Mean Temperature, °C | 0.31 (0.12, 0.5) | < 0.01* | 0.16 (-0.02, 0.34) | 0.08 |
Maximum Temperature, °C | 0.18 (0.08, 0.27) | < 0.01* | 0.091 (0.01, 0.17) | 0.03* |
A Coefficient of log-transformed HSP70 and anti-HSP70 concentrations. |
* Indicates significance at p < 0.05. |
Table 4
Results from multivariable linear mixed-effect models examining associations between serum extracellular Heat Shock Protein 70, (eHSP70), and antibodies to HSP70, (anti-HSP70), and weather indices and markers of kidney function.
| eHSP70 Model | | Anti-HSP70 Model | |
| Estimate A (95% CI) | p-value | Estimate A (95% CI) | p-value |
Age, years | 0.01 (-0.01, 0.02) | 0.48 | 0.004 (-0.02, 0.03) | 0.75 |
Systolic BP, mmHg | -0.002 (-0.02, 0.01) | 0.76 | 0.01 (-0.01, 0.02) | 0.16 |
Diastolic BP, mmHg | 0.01 (-0.01, 0.02) | 0.44 | -0.004 (-0.02, 0.01) | 0.69 |
Maximum Temperature, °C | 0.21 (0.09, 0.33) | < 0.01* | 0.009 (-0.13, 0.11) | 0.99 |
Uric Acid, mg/dL | -0.09 (-0.21, 0.03) | 0.12 | - | |
eHSP70, ng/mL | - | | 0.065 (0.01, 0.13) | 0.03* |
A Coefficient of log-transformed HSP70 and anti-HSP70 level. |
* Indicates significance at p < 0.05. |
In the univariate models for anti-HSP70, increasing systolic blood pressure, maximum temperature, and eHSP70 were significantly associated with increases in anti-HSP70 levels. Decreasing average humidity was associated with increasing anti-HSP70 levels. In the multivariable model, eHSP70 (β: 0.067, 95% CI 0.01, 0.13) remained significant after controlling for age, systolic and diastolic blood pressure. HbA1c (β: -0.664, 95% CI: -1.34, 0.012) and maximum temperature (β: 0.001, 95% CI: -0.12, 0.12) were no longer significant. Of note, the hottest of the four days also had the lowest average humidity, with a maximum temperature of 34.0°C and a humidity of 71%, possibly explaining the observed opposite directions of the β estimates for temperature and humidity.
{Insert Table 3}
{Insert Table 4}
-
Kidney function and eHSP70 and anti-HSP70 levels
Having reduced kidney function (eGFR-Cystatin C < 90 mL/min/1.73 m2) was more common in November than April (23% and 13%, respectively) with a similar trend with eGFR-Creatinine (Table 1). We did not see acute relationships in the repeated measures analysis between serum eHSP70 or anti-HSP70 levels with markers of muscle breakdown (CK) or markers of kidney function (eGFR-Creatinine and eGFR-Cystatin) at the end of the work shift (Table 3).
We examined correlations between eHSP70 and anti-HSP70 levels and markers of kidney function across the harvest. There was a moderate negative correlation between cross-harvest change in eGFR-Cystatin C and April anti-HSP70 levels (r = -0.49, p < 0.01) and cross-harvest change in anti-HSP70 levels (r= -0.38, p = 0.02). Negative cross-harvest change in eGFR-Cystatin C represents worsening in kidney function. Thus, worsening kidney function from November to April was correlated with higher levels of anti-HSP70 levels in April. Additionally, worsening kidney function from November to April was also correlated with increases in anti-HSP70 from November to April. We did not see significant correlations with cross-harvest change in eGFR-Creatinine. We performed post-hoc analyses based on these results to assess differences in anti-HSP70 levels between workers who had declines in eGFR-Cystatin C across the harvest (n = 15, change ranged from − 96.0 to -1.8) and workers who maintained or had increases in eGFR-Cystatin C across the harvest (n = 25, change ranged from 5.4 to 79.0). The 15 workers with declines across the harvest had significantly higher levels of April anti-HSP70 levels (median: 135.18 µg/mL; IQR: 105.48, 183.63) compared to the 25 workers who maintained or had increases in eGFR-Cystatin C (median: 72.3 µg/mL; IQR: 47.86, 109.30), Wilcoxon rank-sum test (p < 0.01). We also observed that those workers with declines in eGFR-Cystatin C across the harvest had greater increases in anti-HSP70 levels across the harvest (median change: 48.72; IQR: 20.87, 74.74) compared to workers who maintained or had increases in eGFR-Cystatin C (median change: 5.27; IQR: -25.61, 38.12), Wilcoxon rank-sum test (p = 0.04).