From May 2005 to April 2020, 11001 were tested for scrub typhus, out of which 2784 were positives. Table 1 describes the demographic features of the study population. The overall positivity rate among the people tested for ST is 25.3 percent. The male to female ratio is 1:1.4. The mean age of scrub typhus positive patients is 44.2 with a standard deviation of 19.2 (Median: 45, IQR: 30, 59). The detailed distribution of cases including gender, occupation and climate is given in supplementary data (Figure S1 to S5).
Table 1: Demographic Characters
Demographic characters
(n = 11001)
|
Total tested
|
Negative
|
Positive
|
Positivity rate (Positive/Total tested)
|
Numbers
|
Percentages
|
Sex
|
Male
|
5570
|
4386
|
1184
|
42.5%
|
21.3%
|
Female
|
5431
|
3831
|
1600
|
57.5%
|
29.5%
|
Total
|
11001
|
8217
|
2784
|
100%
|
25.3%
|
Age group
|
0 - 20 years
|
2036
|
1714
|
322
|
11.6%
|
15.8%
|
21 - 30 years
|
2144
|
1745
|
399
|
14.3%
|
18.6%
|
31 - 40 years
|
1699
|
1254
|
445
|
16%
|
26.2%
|
41 - 50 years
|
1643
|
1144
|
499
|
17.9%
|
30.4%
|
51 - 60 years
|
1588
|
1078
|
510
|
18.3%
|
32.1%
|
>60 years
|
1891
|
1282
|
609
|
21.9%
|
32.2%
|
Total
|
11001
|
8217
|
2784
|
100%
|
25.3%
|
OP/IP
|
IP
|
6902
|
5031
|
1871
|
67.2%
|
27.1%
|
OP
|
4020
|
3186
|
834
|
30%
|
20.8%
|
Missing
|
79
|
0
|
79
|
2.8%
|
NA
|
Total
|
11001
|
8217
|
2784
|
100%
|
25.3%
|
Spatial and Temporal distributions:
The average number of scrub typhus cases reported per year from the entire Vellore district is 186 with the majority of the cases from Katpadi and Vellore blocks with 45 and 44 cases per year respectively (Refer Figure 2). Except for the year May 05 – April 06, Gudiyatham, Katpadi, Vellore and Walajapet blocks have contributed more than 10 scrub typhus cases in all the years.
The mean distance from the patient’s house to CMC hospital ranges from 18.1 km (SD: 19.8 km) in May 19 – April 20 to 31.4 km (SD: 26.5 km) in May 11 – April 12. The overall mean distance in all 15 years from the hospital to scrub typhus patients is 24.8 km (SD: 22.7 km) (Refer Table S3).
Figure 3 (depicted above) shows the time series of scrub typhus positive and the rainfall received, average temperature and humidity of that corresponding month. Increase in cases has occurred only after 2 months of the increased rainfall and it is uniformly observed in all 15 years. The number of ST cases gradually increased with the decrease of the temperature (Max and min) and decreased as the temperature increased (Max and min). This pattern is observed in all 15 years. With relative humidity, the ST cases gradually increased as the maximum and minimum humidity increased and vice versa.
Over the 15 years observed, there is a gradual increase of ST cases from the month of August to February. The ST cases were at peak from September to January and least number of cases were reported from March to July (Refer Figure 4).
In Vellore, the average temperature during the cooler months (August to February) is 26.8°C (Range 22.9 – 31.9°C), the average humidity is 73.2% (Range 63.3 – 83.1%) with a median rainfall of 82.9 mm (IQR: 19.4 mm, 168.5 mm) during this period (Refer Table S1 & S2). The mean temperature in October and November (when scrub typhus cases peak) is 26.9°C (Range 22.6°C – 31.3°C), the average humidity is 78.4% (Range 71.6 – 85.1%) with a median rainfall of 140.5 mm (IQR: 70.9mm, 206.3mm) (Refer Table S1 & S2).
The average temperature during the summer month (March – July) is 31.2°C (Range 24.9 - 37.4°C), the average humidity is 60.1% (Range 48.2 – 71.9%) with a median rainfall of 50.7mm (IQR: 5.7 mm, 97.7 mm) (Refer Table S1 & S2).
Table 2: Month wise correlation between Scrub typhus cases and climate variables
Lag months
|
Rainfall
|
Mean temperature
|
Average humidity
|
r
|
95% CI
|
r
|
95% CI
|
r
|
95% CI
|
L_0
|
0.186*
|
0.047, 0.322
|
-0.698**
|
-0.773, -0.608
|
0.733**
|
0.663, 0.790
|
L_1
|
0.586**
|
0.492, 0.656
|
-0.452**
|
-0.555, -0.326
|
0.700**
|
0.626, 0.755
|
L_2
|
0.751**
|
0.678, 0.809
|
-0.037
|
-0.167, 0.100
|
0.424**
|
0.306, 0.527
|
L_3
|
0.628*
|
0.522, 0.717
|
0.383**
|
0.269, 0.499
|
0.031
|
-0.109, 0.164
|
L_4
|
0.363**
|
0.237, 0.492
|
0.656**
|
0.557, 0.739
|
-0.343**
|
-0.449, -0.225
|
L_5
|
0.051
|
-0.089, 0.193
|
0.732**
|
0.651, 0.792
|
-0.595**
|
-0.673, -0.490
|
**Spearman’s Correlation is significant at the 0.01 level (2-tailed).
*Spearman’s Correlation is significant at the 0.05 level (2-tailed).
For better understanding of the association between scrub typhus and climatic condition, Spearman’s correlation was done to find out the association between the time taken for a spurt in ST cases (lag time) and rainfall, temperature and humidity for up to a lag of 5 months (Refer Table 2). As observed earlier in the month-wise trend for 15 years, rainfall had the maximum positive correlation of 0.751 (95% CI: 0.678, 0.809) with a lag of two months and it is statistically significant (p <0.001). This positive association, however, reduced in the third and fourth month and became insignificant at the lag time of five months.
Mean temperature showed a negative correlation of -0.698 initially and reduced in lag of months and became insignificant in three lag months. However, it showed a positive correlation from 4 lag months.
Average maximum humidity showed a positive correlation with no lag and it continued till 2 lag months. Three lag months showed an insignificant correlation with four and five lag months showing a negative correlation.
Table 3: Spearman correlation coefficient (‘r’) matrix of meteorological variables
|
Rainfall_2monthslag(mm)
|
Mean temperature (°C)
|
Mean relative humidity (%)
|
Rainfall_2monthslag (mm)
|
1
|
-
|
-
|
Mean temperature (°C)
|
-0.612**
|
1
|
-
|
Mean relative humidity (%)
|
0.618**
|
-0.739**
|
1
|
All the p values are significant at 0.001 level
Table 3 describes the correlation coefficient matrix (Spearman) of the meteorological factors such as rainfall, mean temperature and mean relative humidity. Since the rainfall with 2 months lag has the highest correlation, we considered that for computing the Spearman correlation coefficient matrix and negative binomial regression. Whereas for the mean temperature and mean relative humidity, the same month reading with no lag was used. Mean temperature and mean relative humidity seem to be strongly correlated (-0.739) with each other.
Table 4: Negative Binomial regression analysis of climatic factors associated with scrub typhus in Vellore from May 2005 to April 2020
Model
|
Variables
|
β
|
SE
|
p
|
Percent increase = (eβ-1)*100
|
95% CI
|
Lower
|
Higher
|
A
|
Intercept
|
7.756
|
1.0205
|
<0.001
|
|
Rainfall_2monthslag (mm)
|
0.007
|
0.0012
|
<0.001
|
0.70
|
0.40
|
0.90
|
Mean temp (°C)
|
-0.208
|
0.0343
|
<0.001
|
-18.78
|
-24.12
|
-13.15
|
B
|
Intercept
|
-3.059
|
0.675
|
<0.001
|
|
Rainfall_2monthslag (mm)
|
0.005
|
0.0012
|
<0.001
|
0.50
|
0.30
|
0.70
|
Mean RH (%)
|
0.073
|
0.0104
|
<0.001
|
7.57
|
5.44
|
9.86
|
Note: *Negative binomial regression model for monthly scrub typhus without mean relative humidity (A) and model with rainfall, mean humidity and without mean temperature (B). CI = Confidence interval, S.E. = Standard error
To address the multicollinearity issues between mean temperatures and mean relative humidity in regression analysis, two different models (A & B) were developed using negative binomial regression for exploring the relation between mean temperature and mean relative humidity among ST cases (Refer Table 4). Model A was computed using the independent variables like rainfall_2monthslag and mean temperature. Model B with variables like rainfall_2months lag and mean humidity. For an increase of 1°C in mean temperature, the monthly ST case reduces by 18.78% (95% CI: -24.12, -13.15%). On the contrary, for 1 percent increase in mean relative humidity, there is an increase of 7.57% (95% CI: 5.44, 9.86%) of monthly ST cases. Similarly, for an increase of 1mm of rainfall contributes to increase of 0.50 to 0.70% of monthly ST cases after 2 months depending on the variables included in the analysis.
Table 5: Factors associated with Scrub Typhus
Variables
|
ST Positive
|
ST Negative
|
OR (95% CI)
|
p value
|
AOR (95% CI)
|
p value
|
Gender
|
Female
|
1601 (29.5%)
|
3832 (70.5%)
|
1.55 (1.42 – 1.69)
|
<0.001
|
1.16 (1.05 – 1.30)
|
0.006
|
Male
|
1183 (21.2%)
|
4385 (78.8%)
|
Age category
(years)
|
0 – 20 (Ref)
|
322 (15.8%)
|
1714 (84.2%)
|
Reference
|
21 – 30
|
399 (18.6%)
|
1745 (81.4%)
|
1.22 (1.04 – 1.43)
|
0.017
|
0.67 (0.55 – 0.80)
|
<0.001
|
31 – 40
|
445 (26.2%)
|
1254 (73.8%)
|
1.89 (1.61 – 2.22)
|
<0.001
|
0.99 (0.82 – 1.20)
|
0.911
|
41 – 50
|
499 (30.4%)
|
1144 (69.6%)
|
2.32 (1.98 – 2.72)
|
<0.001
|
1.13 (0.93 – 1.37)
|
0.224
|
51 – 60
|
510 (32.1%)
|
1078 (67.9%)
|
2.52 (2.15 – 2.95)
|
<0.001
|
1.23 (1.02 – 1.50)
|
0.034
|
>60
|
609 (32.2%)
|
1282 (67.8%)
|
2.53 (2.17 – 2.95)
|
<0.001
|
1.37 (1.13 – 1.65)
|
0.001
|
Occupation
|
Farmers, daily wages & house wife
|
1580 (33.8%)
|
3100 (66.2%)
|
2.74 (2.47 – 3.04)
|
<0.001
|
2.43 (2.13 – 2.77)
|
<0.001
|
Others
|
665 (15.7%)
|
3577 (84.3%)
|
OP/IP
|
In-patient
|
1871 (27.1%)
|
5031 (72.9%)
|
1.42 (1.30 – 1.56)
|
<0.001
|
1.11 (0.99- 1.24)
|
0.058
|
Out-patient
|
834 (20.7%)
|
3186 (79.3%)
|
Univariate & Multivariate analysis:
Table 5 describes the association of ST with the risk factors analyzed using chi-square and reported with 95% CI of Odds Ratio (OR). In univariate analysis, females, age group of >20 years, farmers, daily wage workers and house wives had the higher odds to be positive for scrub typhus when compared to males, age group 0 - 20 years and other occupations. In-patients also had a higher risk of being ST positive when compared to patients tested as out-patient. However, in multivariate analysis, females, age group of 51 – 60 years and >60 years, farmers, daily wage workers and house wives were at higher odds for being ST positive and the age group of 21 – 30 years showed a protective effect against ST when compared with 0 – 20 years of age.