Study population
A total of 1,977 individuals (215 from Honiara, 221 from Guadalcanal, 392 from Western Malaita, and 996 from Isabel) participated in the study. Participants had a median age of 30 years, with 62% being female (Table 1). The average tympanic temperature of participants during the survey was 37.1°C, and 34 people had a temperature exceeding 38°C at the time of the survey, with the maximum temperature recorded being 40.6°C.
Table 1
Study population summary characteristics.
Characteristic
|
Summary
|
Survey dates
|
Apr 2018
|
Number of participants
|
1,977
|
Number of samples analysed
|
1,914
|
Age – Range
|
5–86 years
|
Age – Median
|
30 years
|
Percentage female
|
62% (n = 1,229)
|
P. falciparum prevalence
|
1.2% (n = 23/1914)
|
P. vivax prevalence
|
1.5% (n = 28/1914)
|
Percentage P. vivax
|
61% (n = 28/46)
|
Vector control and arbovirus protection
Insecticide-treated bed net use (predominately long-lasting insecticide treated nets (LLINs)) the night prior to the survey was reported by 63% of participants and significantly differed by province (χ2 = 141.49 df = 4, p < 0.0001; Fig. 2), ranging from 48% in Malaita Province to 75% in Isabel Province. The use of ITNs also varied by villages from 29% in Lilisiana, Malaita to 95% in Hovukoilo Village in Isabel (Table S1). There was limited use of other mosquito protective measures; 10% of study participants used mosquito coils (usually a volatile pyrethroid) and 5% of survey participants’ houses had window screens. Topical repellents were used by only 0.8% of survey participants. The use of other mosquito control measures was significantly related to the province (χ2 = 301.27 df = 4, p < 0.0001), with higher mosquito coil use in Honiara and Guadalcanal (39% and 29%, respectively), and 16% of participant houses in Western Malaita having window-screens (Fig. 2).
Plasmodium prevalence by PCR
Blood samples from 1,914 participants were analysed by PCR for Plasmodium spp .DNA (henceforth, DNA malaria positive participants are referred to as malaria-positive and used to determine malaria prevalences): 46 participants were malaria positive with 17 individuals positive for P. falciparum, 22 with P. vivax, 6 with both P. falciparum and P. vivax, and a single individual with P. ovale (Table 2). Of Plasmodium PCR positive individuals, seven were febrile (temperature ≥ 38°C): 2 were positive for P. falciparum, 2 were positive for P. vivax and 3 were PCR positive with both P. falciparum and P. vivax; Table 3). Thus, the percentage of febrile or symptomatic participants for P. falciparum was 22% (n = 5/23) and for P. vivax was 17% (n = 5/28).
Table 2
The prevalence of Plasmodium DNA-positives across Provinces in the Solomon Islands among residents of all ages above 5 years.
|
|
Number positive
|
Prevalence
|
Province
|
Participants
|
Mixed*
|
Pf
|
Pv
|
Po
|
Pf
|
Pv
|
Overall
|
Honiara
|
211
|
0
|
0
|
4
|
0
|
0.0%
|
1.9%
|
1.9%
|
Guadalcanal
|
369
|
1
|
0
|
3
|
0
|
0.3%
|
1.1%
|
1.1%
|
Isabel
|
946
|
0
|
0
|
0
|
0
|
0.0%
|
0.0%
|
0.0%
|
Malaita
|
388
|
5
|
17
|
15
|
1
|
5.7%
|
5.2%
|
9.8%
|
Overall
|
1,914
|
6
|
17
|
22
|
1
|
1.2%
|
1.5%
|
2.4%
|
*The mixed infections contained both P. falciparum and P. vivax |
Table 3
The number and percentage of participants that were positive for either P. falciparum, or P. vivax summarised by the various explanatory variables.
|
|
P. falciparum
|
P. vivax
|
Parameter
|
Total
|
n
|
%
|
n
|
%
|
Gender
|
|
|
|
|
|
Female
|
1197
|
14
|
1.2%
|
14
|
1.2%
|
Male
|
717
|
9
|
1.3%
|
14
|
1.9%
|
Fever
|
|
|
|
|
|
Yes
|
34
|
5
|
14.7%
|
5
|
14.7%
|
No
|
1880
|
18
|
0.9%
|
23
|
1.2%
|
Domestic travel history
|
|
|
|
Yes
|
79
|
0
|
0.0%
|
2
|
2.5%
|
No
|
1835
|
23
|
1.3%
|
26
|
1.4%
|
International travel history
|
|
|
Yes
|
49
|
0
|
0.0%
|
1
|
2.0%
|
No
|
1865
|
23
|
1.2%
|
27
|
1.4%
|
Bednet use
|
|
|
|
|
|
Yes
|
1275
|
11
|
0.9%
|
6
|
0.5%
|
No
|
639
|
12
|
1.9%
|
22
|
3.4%
|
Malaria history
|
|
|
|
|
Yes
|
884
|
14
|
1.6%
|
15
|
1.7%
|
No
|
1016
|
9
|
0.9%
|
13
|
1.3%
|
Medicine use
|
|
|
|
|
Yes
|
33
|
0
|
0.0%
|
1
|
3.0%
|
No
|
1881
|
23
|
1.2%
|
27
|
1.4%
|
For P. falciparium, the base GLM model was most substantially improved by adding village (100% wAIC support). Sequentially the model was improved by adding temperature (93% wAIC support, Fig. 3). These explanatory variables of village and temperature were significant (log-likelihood ratio test) and were included in the most parsimonious model (Table 4). None of the other remaining candidate factors were able to further improve model fit. Of note is that although the prevalence of P. falciparum was reduced almost by half from 1.9–0.9% by using an ITN (Table 3), this factor did not explain sufficient variation to be included in the final model.
Table 4
Final set of nested models evaluated to determine which best predicted the prevalence of Plasmodium DNA-positives.
Model
|
df
|
AIC
|
∆AIC
|
wAIC
|
χ2
|
p value
|
P. falciparum
|
|
|
|
|
|
|
Village
|
19
|
183.91
|
7.31
|
0.0147
|
|
|
Village + Temperature
|
20
|
176.59
|
0.00
|
0.5694
|
9.31
|
0.0022*
|
Village + Temperature +
Bednet use
|
21
|
177.22
|
0.63
|
0.4157
|
1.37
|
0.241
|
P. vivax
|
|
|
|
|
|
|
Village
|
19
|
253.64
|
23.87
|
< 0.0001
|
|
|
Village + Temperature
|
20
|
241.25
|
11.48
|
0.0023
|
14.39
|
0.0001*
|
Village + Temperature +
Bednet use
|
21
|
231.74
|
1.97
|
0.2714
|
11.51
|
0.0007*
|
Village + Temperature +
Bednet use + Gender
|
22
|
229.77
|
0.00
|
0.7261
|
3.98
|
0.0463*
|
Model comparison was made on the basis ∆AIC, wAIC and goodness-of-fit using maximum likelihood estimation. The full list of explanatory variables included village, gender, temperature, age, bednet use and malaria history. |
Much of the variation in P. falciparum prevalence was explained by village, and this species was extremely heterogeneous across the provinces, with only 2 small foci of transmission identified among 19 villages surveyed (Fig. 1). Plasmodium falciparum DNA positive samples were extremely localised, there was 1 mixed positive sample from Guadalcanal (n = 119) and the remaining 22 positive participants were from only 2 villages in Malaita. Overall, 69% of the P. falciparum positives were from a single village, where the prevalence was 15.5% (Fig. 1). When compared at the provincial level, P. falciparium prevalence was highest in Malaita (Table 2).
For P. vivax, the base GLM model was most substantially improved by adding village (99% wAIC support). Sequentially the model was improved by adding temperature (87% wAIC support, Fig. 3), bednet use (94% wAIC support) and gender (51% wAIC support). These explanatory variables of village and temperature were significant (log-likelihood ratio test), and were included in the most parsimonious model (Table 4). Adding age or malaria history was unable to further improve the model fit. While age was not included in the most parsimonious P. vivax model, a univariate GLM did pick up an influence of age on infection (χ2 = 9.63, df = 1, p = 0.0019; Fig. 3). For gender, males were more likely to be positive for P. vivax (Table 3).
The P. vivax positive samples were heterogeneous across the villages, but they were not as extremely localised as the P. falciparum infections. Plasmodium vivax DNA was detected in participants from Guadalcanal (n = 4), Honiara (n = 4) and Malaita (n = 20), across 7 villages. Overall 63% of the P. vivax positive individuals were found in 2 villages, where the prevalence was 8.0% and 9.7% (Fig. 1). Plasmodium spp. DNA was not detected from Isabel province participants. In the current study, statistical hotspots were not defined by geospatial statistics due to the difficulty to show statistical significance at such very low transmission intensities where only isolated foci remain.