Skip to main content

Spatial-temporal heterogeneity in malaria receptivity is best estimated by vector biting rates in areas nearing elimination

Abstract

Background

Decisions on when vector control can be withdrawn after malaria is eliminated depend on the receptivity or potential of an area to support vector populations. To guide malaria control and elimination programmes, the potential of biting rates, sporozoite rates, entomological inoculation rates and parity rates to estimate malaria receptivity and transmission were compared within and among geographically localised villages of active transmission in the Western Province of the Solomon Islands.

Results

Malaria transmission and transmission potential was heterogeneous in both time and space both among and within villages as defined by anopheline species composition and biting densities. Biting rates during the peak biting period (from 18:00 to 00:00 h) of the primary vector, Anopheles farauti, ranged from less than 0.3 bites per person per half night in low receptivity villages to 26 bites per person in highly receptive villages. Within villages, sites with high anopheline biting rates were significantly clustered. Sporozoite rates provided evidence for continued transmission of Plasmodium falciparum, P. vivax and P. ovale by An. farauti and for incriminating An. hinesorum, as a minor vector, but were unreliable as indicators of transmission intensity.

Conclusions

In the low transmission area studied, sporozoite, entomological inoculation and parity rates could not be measured with the precision required to provide guidance to malaria programmes. Receptivity and potential transmission risk may be most reliably estimated by the vector biting rate. These results support the meaningful design of operational research programmes to ensure that resources are focused on providing information that can be utilised by malaria control programmes to best understand both transmission, transmission risk and receptivity across different areas.

Background

Globally, malaria transmission has fallen significantly with 68% of the reduction in Plasmodium falciparum in Africa attributed to the use of long-lasting insecticide-treated nets (LLINs) [1]. Recently, the World Health Assembly endorsed the Global Technical Strategy for Malaria Control and Elimination (GTS) [2]. The GTS recommends universal access to vector control with LLINs or indoor residual spraying (IRS) to all people at-risk of malaria. Larval source management (including larviciding, insect growth regulators and environmental management) is also recommended in the GTS as a supplemental control measure where larval habitats are few in number, fixed in location and easily accessible. Maintenance of universal access to LLINs or IRS after elimination is recommended in areas both receptive and vulnerable to malaria [2]. Withdrawal of universal access to LLINs or IRS after malaria is eliminated will depend on the risk of resumption of transmission if the malaria parasite is re-introduced, a function of the movement of infected people or mosquitoes (e.g. vulnerability) and vector receptivity (e.g. an environment inherently capable of supporting significant vector populations) [3]. Receptivity is not static and can change with urbanization, alterations in land use patterns and implementation of interventions that permanently reduce the vectorial capacity. Hence, there is a need to better understand how to quantify receptivity to guide malaria elimination programmes.

Although the Solomon Islands has achieved significant reductions in malaria transmission across the past decade, the number of cases has increased since 2015 and the annual parasite incidence was 83.4 cases/1000 population in 2017 (Solomon Islands Ministry of Health and Medical Services, unpublished data). In Western Province, the annual parasite incidence was only 7/1000 population in 2015 and foci of transmission emerged. Concurrently, the malaria parasite ratio changed from predominantly P. falciparum to P. vivax. Stratifying areas by receptivity will be important for targeting resources to where they are required to maintain malaria elimination as well as to respond rapidly to outbreaks [4, 5]. For entomological monitoring of adult anophelines in foci investigations, determining the species composition, receptivity and insecticide resistance are the highest priorities with moderate emphasis on determining the human biting rate, biting time and location [6].

Malaria in the Solomon Islands is transmitted almost exclusively by An. farauti. Following exposure to DDT applied in IRS during the original Malaria Elimination Programme of the 1970s, An. farauti shifted its biting profile from all night with both indoor and outdoor biting to predominantly biting early in the evening and outdoors [7,8,9]. While this behavioural change enables this mosquito to minimize contact with the WHO recommended interventions (LLINs and IRS) that are applied inside houses, these tools still retain efficacy. This is due to An. farauti having a short feeding cycle and individual mosquitoes must complete 5–6 feeding cycles to live long enough to complete the extrinsic incubation period and be infectious. During this time, they are likely to feed indoors at least once, and so, LLINs and IRS can potentially kill a significant proportion of the population [10]. The impacts of LLINs and IRS on vectorial capacity are temporary, and premature withdrawal of these interventions, in the absence of other interventions to permanently reduce vectorial capacity, will leave such areas susceptible to malaria resurgence [3, 11,12,13].

In this context, the utility of different entomological indicators to estimate malaria receptivity, transmission and transmission potential across time and space were directly compared, being biting rates, survivorship (by parity dissections), sporozoite rates and entomological inoculation rates of all human-biting anophelines. The comparison was made across 11 villages that spanned a malaria foci in the Western Province of the Solomon Islands.

Methods

Study sites

The study was conducted in 11 coastal villages located on the volcanic, mountainous, rain-forested islands of Ghizo, Kohinggo, Kolombangara, New Georgia and Ranonnga Islands of Western Province of the Solomon Islands (-8°0'S, 157°0'E) [14]; as well as in Haleta Village on Ngella Sule Island in Central Province (-9°0'S, 159°45'E) (Fig. 1) [10, 14,15,16]. Prior to commencing the entomological survey, the malaria prevalence across the study villages in Western Province was estimated to be 2.1% (measured by PCR in 2013, unpublished data). The area is classified by the Solomon Islands Government as a “near elimination” province. In contrast, Central Province had a malaria prevalence of 13.4% (as measured by PCR in 2012) [17]. Spatial analyses of malaria prevalence in humans from a cross-sectional survey conducted in 2013 identified a significant foci of infection that encompassed the villages of Jack Harbour and Tugivili, but not the villages of Nazareth, Kinamara or Saeragi with New Mala situated on the periphery of the foci (unpublished data).

Fig. 1
figure 1

Map of (a) the Solomon Islands showing (b) the 11 study villages in Western Province (-8°0'S, 157°0'E) and (c) Central Province (-9°0'S, 159°45'E)

Western Province encompasses ≈5000 km2 with a population of 76,649 in 13,762 households [18]. The climate of the region is hot and wet with an annual rainfall of 3725 mm (mean from 1999 to 2010, Munda Airport, New Georgia Island in the Western Province and 2837 mm in the Central Province; Solomon Islands Bureau of Meteorology, unpublished data). Sixty-eight per cent of study village residents self-reported sleeping under a long-lasting insecticide treated net (unpublished data). In both provinces, the mean daily minimum and maximum temperatures were 24 °C and 30 °C, respectively, with an overall mean of 26 °C.

Study period

Unless otherwise specified, anopheline biting densities were estimated in the Western Province between March 2014 and August 2016 for 4 nights per village survey (n = 2064 man-nights), and in Haleta (Central Province) between August 2011 and August 2016 for 5 nights per survey (n = 1534 man-nights). Entomological data reported previously for Haleta in Central Province [10, 14,15,16] is updated here for comparative analyses of vector species composition and heterogeneity in biting rates among and within villages.

Sampling of adult anophelines

Host-seeking (biting) females were sampled with human landing catches (HLC) conducted from 18:00 to 00:00 h by village collectors working outdoors at 10 sites distributed throughout each village. To determine the all-night biting profile, mosquito collections were extended to 06:00 h in Jack Harbour and Saeragi during September 2014. Anophelines landing on the exposed legs and feet of collectors were captured by mouth aspiration and held in individual containers by hour and collection site. Anophelines were identified by morphological criteria [19], prior to dissection for parity determination [20]. Specimens were preserved in 100% ethanol for subsequent PCR-based identification using a DNA sequence of the internal transcribed spacer region 2 of the ribosomal DNA (ITS2) [21] and detection of Plasmodium DNA in heads and thoraces by nested PCR [22].

The annual entomological inoculation rate (EIR) was calculated from the product of the sporozoite rate and the annual biting rate [23, 24]. The sporozoite rate was defined as the proportion of mosquitoes with malaria specific DNA in the head or thorax. The all-night biting rate was calculated by adjusting the estimated biting rate from 18:00–24:00 h to account for the proportion of females estimated to have fed after midnight; based on the all-night collections conducted for An. farauti in Jack Harbour and An. lungae in Saraegi village.

Statistical analysis

Data detailing mosquito surveys and their analyses by dissections and molecular analysis are available from the James Cook University Tropical Data Hub [25, 26]. Differences in the anopheline community composition among villages were analysed by permutational multivariate ANOVA (PERMANOVA; package vegan) [27] and displayed graphically using non-metric multidimensional scaling (nMDS) [28]. Temporal and spatial changes in each species biting rates were compared with a generalized linear model (GLM; package = MASS) with a negative binomial distribution and interacting fixed factors for sample period and village. All analyses were conducted using the R package V3.1.2 [29].

Geographical data were projected in ArcGIS (v10.0) [30] and local spatial clusters of high mosquito densities were detected using FleXScan (v3.1.2) [31] which can identify either circular or irregular shaped clusters [32]. The flexible scan statistic uses a purely spatial Poisson distribution model to identify spatially aggregated clusters with higher than average mosquito densities (“vector foci”) by identifying the spatial window with the greatest ratio of observed to expected cases (Relative Risk). Cluster detection was based on a spatial matrix [33] defined using triangular irregular networks created based on Delaunay Triangulation, with Euclidian distance.

Results

In the Western Province, mosquito surveys were conducted longitudinally in Jack Harbour (n = 10), New Mala (n = 9), Saeragi (n = 9), Kinamara (n = 8), Nazareth (n = 8), Obobulu (n = 5) and Tuguivili (n = 5) and once each in Boboe, Iriri, Kuzi and Koriovuku for a total of 2064 man-nights (each survey consisted of 4 nights of collections with 10 collectors per night; Table 1). Members of both the An. farauti (s.l.) (n = 11,516) and An. lungae (s.l.) (n = 187) complexes were collected outdoors with HLC. PCR analyses estimated that of those morphologically identified to belong to the An. farauti complex; 91% were An. farauti (1376/1520) and 9% An. hinesorum (144/1520; Fig. 2). Of the An. lungae complex specimens, 93% were confirmed by PCR as being An. lungae (178/191) and 7% were An. solomonis (13/191; Fig. 2). The species composition was significantly different among villages, both as a main effect (PERMANOVA, F(1,75) = 2.37, P = 0.001) and as an interaction with time (F(1,75) = 0.88, P = 0.011). Although the species composition was stable over time in some villages, e.g. Jack Harbour, there was strong temporal variation in others, e.g. Kinamara (Fig. 2). The anopheline communities ranged from being exclusively An. farauti in Jack Harbour (and Haleta in the Central Province) to a dominance of An. lungae in Nazareth and Saeragi with other villages (Kinamara) co-dominated by two species. In the ordination plot, the only village that consistently separated on a different gradient was Kinamara, which had varying mixed populations including proportionally elevated numbers of An. hinesorum and An. solomonis (Fig. 3).

Table 1 Timeline of anopheline surveys in Western Province, Solomon Islands
Fig. 2
figure 2

Temporal species composition of Anopheles farauti, Anopheles hinesorum, Anopheles lungae and Anopheles solomonis from villages in Western Province. Total numbers of An. farauti, An. hinesorum, An. lungae and An. solomonis identified by PCR are shown in blue, red, green and orange, respectively. Analyses of the Dec 2014 collection in Jack Harbour confirmed all anophelines were An. farauti (n = 96). Analyses across all collections in Obobulu confirmed 86 % An. farauti (n = 19) and 14 % An. hinesorum (n = 3). Analyses across all collections in Tuguivili confirmed 91 % An. farauti (n = 60), 3 % An. hinesorum (n = 2) and 6 % An. lungae (n = 4). Key: ND, no data; 0, no specimens were caught

Fig. 3
figure 3

Non-metric multidimensional scaling ordination plot comparing the species abundance from different sample periods in study villages in Western (Jack Harbour, Kinamara, Nazareth, Obobulu, Saeragi) and Central (Haleta) Provinces. Each point represents the species composition of one village at one sampling period, and those that are more similar to one another are ordinated closer together. The axis and orientation of the plot is arbitrary

With the exception of Jack Harbour where 10,530 An. farauti were captured by HLC, An. farauti, An. hinesorum, An. lungae and An. solomonis populations were not abundant in the other villages (Fig. 4). Anopheles farauti biting densities from 18:00 to 24:00 h were estimated as 26.3 bites per person per half-night (b/p/h-n) in Jack Harbour and 1.5 b/p/h-n in New Mala, with an average of 0.2 b/p/h-n in all other villages surveyed ≥ 8 times (Kinamara, Nazareth and Saeragi). Large heterogeneities in biting rates were observed for An. farauti over time (β = -0.303, SE = 0.072, P < 0.0001) and among villages (β = -0.982, SE = 0.075, P < 0.0001), noting that there was a significant interaction between time and village (β = 0.063, SE = 0.011, P < 0.0001). For example, the temporal variability of the average biting rate in Jack Harbour across the 10 collection periods ranged from a low of 1.2 b/p/h-n in August 2016 to a high of 73.7 b/p/h-n in June 2014. While the average biting rate of An. hinesorum, An. lungae and An. solomonis among all villages was always very low at 0.1, 0.1 and 0.008 b/p/h-n, respectively, although at times significant increases in the biting rates were observed for An. hinesorum (β = -0.414, SE = 0.162, P = 0108) and An. lungae (β = 0.489, SE = 0.209, P = 0.0194). However, An. solomonis densities remained consistently low (β = -0.112, SE = 0.954, P = 0.906). At such low densities, the GLM model did not detect any influence of village on the landing rates of An. hinesorum (β = -0.225, SE = 0.190, P = 0.237) or An. solomonis (β = -0.471, SE = 1.274, P = 0.711), whereas the densities of An. lungae were influenced by village (β = 1.014, SE = 0.155, P < 0.0001) with a significant interaction with time (β = -0.116, SE = 0.026, P < 0.0001). Seasonal patterns for any of the species or villages were not obvious (due to low densities), except for An. farauti in Jack Harbour, where a peak biting season occurred between April and June.

Fig. 4
figure 4

Longitudinal densities of Anopheles farauti, Anopheles hinesorum, Anopheles lungae and Anopheles solomonis in Jack Harbour, Kinamara, Nazareth, New Mala and Saeragi villages, Western Province, Solomon Islands, estimated by human landing catches from 18:00–00:00 h. Y-axes scale varies by anopheline species

Five entomological surveys were made in Obobulu (March 2014 to August 2015) and Tuguivili (August 2015 to August 2016). In Obobulu, the mean density of An. farauti was 0.12 b/p/h-n and for An. hinesorum was 0.02 b/p/h-n, with neither An. lungae nor An. solomonis captured. In Tuguivili, the mean density of An. farauti was 0.56 b/p/h-n, for An. hinesorum was 0.03 b/p/h-n and for An. lungae was 0.006 b/p/h-n without any An. solomonis caught. Additional single entomological surveys were made in Koriovuku (June 2014), Iriri (October 2014), Kuzi (December 2014) and Boboe (July 2015) villages. In Korivuku, Iriri and Boboe neither An. farauti (s.l.) nor An. lungae (s.l.) were captured from 18:00 to 24:00 h. In Kuzi only An. farauti (s.l.) was collected (1.7 b/p/h-n).

The heads and thoraces of An. farauti (n = 1921), An. hinesorum (n = 39), An. lungae (n = 149) and An. solomonis (n = 13) were analysed by PCR for Plasmodium DNA as an indicator of malaria sporozoites in the salivary glands (Table 2). Hereafter, PCR positives from heads and thoraces will be referred to as sporozoite DNA. Plasmodium falciparum sporozoite DNA positive An. farauti were identified in Jack Harbour (n = 10) and New Mala (n = 6). Sporozoites of P. vivax were identified in An. farauti in Jack Harbour (n = 3) and New Mala (n = 2) with a single sporozoite positive P. ovale detected in An. farauti in Jack Harbour. The head and thorax of one An. hinesorum from Kinamara was positive for P. falciparum DNA. The overall sporozoite rate for An. farauti was 1.1% (n = 22/1921) and for An. hinesorum was 2.6% (n = 1/39; Table 2). Plasmodium DNA was not detected in An. lungae (n = 149) or An. solomonis (n = 13). For An. farauti, the overall annual EIR was estimated to be 26.6 infective bites/person/year (ib/p/y; Table 3) and ranged from 16.5 in 2014 to 55.5 in 2015. For An. hinesorum, the overall estimated EIR was 1.0.

Table 2 Sporozoite-positive Anopheles farauti and Anopheles hinesorum by village in the Western Province
Table 3 Estimated malaria transmission rates by members of the Anopheles farauti complex in Western Province

Peak biting of An. farauti occurred from 19:00 to 20:00 h (Fig. 5). The percentage of overall biting that occurred before 21:00 h was 49% and before midnight was 59%. For An. lungae, peak biting was even earlier (between 18:00–19:00 h) (Fig. 5). The percentage of overall biting by An. lungae that occurred before 21:00 h was 45% with 55% of biting before midnight.

Fig. 5
figure 5

The hourly outdoor biting profiles of An. farauti in Jack Harbour village during September 2014 (top) and An. lungae in Saeragi village during June 2014 (bottom), Western Province, Solomon Islands

Only Jack Harbour village had adequate numbers of An. farauti to allow multiple time-point parity rate estimations. The overall parity rate of An. farauti was 0.50 (393 parous of 785 dissected, pooling 6 surveys between May 2014 and August 2016). Parity varied significantly by time (β = -0.293, SE = 0.064, P < 0.0001), with parity estimates ranging between 0.33 and 0.70. Interestingly, the An. farauti parity rate in New Mala village was also 0.5 from 254 dissections across three survey periods.

Small-scale spatial clustering (foci) of An. farauti densities was investigated in Jack Harbour, New Mala and Haleta; and for An. hinesorum in Kinamara and An. lungae in Saeragi. High spatial variations in biting rates within a village were observed during single surveys. In Jack Harbour village, for example, mean nightly biting catches among the 10 collection sites ranged from 1 to 129 b/p/h-n during the June 2014 survey. In each village, foci of higher than average densities of either An. farauti, An. hinesorum or An. lungae were identified (Fig. 6, Table 4). For An. farauti, the foci contained 4 to 7 sampling sites with the maximum distance across foci ranging from 170 m in Haleta to 558 m in Jack Harbour. Between 47% and 92% of anophelines were captured in foci. The An. hinesorum focus in Kinamara with a maximum size of only 126 m contained 4 sampling sites and accounted for 67% of all An. hinesorum captured. The An. lungae focus in Saeragi with a maximum size of only 52 m contained 2 sampling sites but accounted for 44% of all An. lungae captured.

Fig. 6
figure 6

Spatial distribution and clustering of An. farauti densities in the Western Province villages of Jack Harbour and New Mala and the Central Province village of Haleta, as well as An. hinesorum in Kinamara village and An. lungae in Saeragi village. The panel labels (a-e) refer to the locations of each village on the regional maps. The scales differ by village and represent the total number of female mosquitoes caught in all sampling periods by HLC site

Table 4 Spatial clusters (foci) of An. farauti densities within Jack Harbour, New Mala and Haleta villages as well as An. hinesorum in Kinamara village and An. lungae in Saeragi village

Discussion

Understanding vector species composition and abundance are important considerations for stratifying areas for targeting malaria vector control, particularly as transmission diminishes. The abundance of competent vector species in a suitable climate defines the relative receptivity of a strata for malaria transmission [34]. Receptivity coupled with an influx of parasites (which defines the malaria vulnerability) will guide malaria elimination programmatic decisions on when and where universal access to vector control can be reduced without increasing the risk of re-establishment of transmission after reintroduction of malaria parasites [11]. Malaria transmission in the Solomon Islands with one dominant vector provides a “simple” scenario to estimate transmission, transmission potential and receptivity and the underlying determinants. Anopheline human biting rates and sporozoite infections were analysed from longitudinal surveys among and within villages to define spatial and temporal heterogeneities associated with a malaria focus.

Four anopheline species (An. farauti, An. hinesorum, An. lungae and An. nataliae) were previously found during larval surveys in the Western Province from February to May 2013 [14]. While An. farauti was the predominant anopheline collected in HLCs, only 18% of larvae were identified as An. farauti from the seven villages in which both larval surveys and longitudinal monitoring of adult anophelines was conducted. Anopheles hinesorum was infrequently collected by HLC but was the most commonly collected larval species (66% of specimens) and was found in six of the seven villages, being only absent from Jack Harbour [14]. Particularly noteworthy was the anopheline species composition of Kinamara village in which An. hinesorum was the most frequent human biter (this species was previously reported as a non-human biter in the Solomon Islands [35, 36]). The ecology of Kinamara was unique among the study villages as multiple small fast-flowing steams are found within the village with thick emergent vegetation along the borders of these streams providing harbourage for An. hinesorum larvae.

In contrast to some previous studies, An. lungae and An. hinesorum were collected biting humans in the Solomon Islands, albeit in small numbers [10, 37, 38]. For An. lungae, this behaviour was widespread (though infrequent) in the Western Province, occurring in villages on Ranonngga, Vonavona, Gizho and Kolobangara islands. Anopheles hinesorum was infrequently collected by HLC in the villages of Kinamara, Nazareth, New Mala, Obobulu and Saeragi, despite being the most frequently collected anopheline in larval surveys [14]. Anopheles solomonis (which was not identified during larval surveys) was collected in limited numbers biting humans in the coastal villages of Kinamara, Nazareth and Saeragi. Human-biting An. solomonis was previously reported from an inland village on Santa Isabel Island [39]. The vector status of An. lungae and An. solomonis remains unresolved (sporozoites were not detected, but the sample sizes were inadequate for any definitive conclusions).

Anopheles hinesorum is a significant vector of P. falciparum, P. vivax and P. malariae in Papua New Guinea [40]. In the Solomon Islands, An. hinesorum was previously reported as an animal biter and therefore not a vector of human malarias [39] - until this study found P. falciparum DNA in the head and thorax of a single An. hinesorum, establishing this species as susceptible to infection with P. falciparum. However, its infrequent human-biting habit may limit its role in malaria transmission in the Solomon Islands.

Evidence for ongoing transmission of P. falciparum, P. vivax and P. ovale by An. farauti was found in the Western Province, a low malaria transmission/near elimination province. While 76% of human infections in the Western Province are Plasmodium vivax (unpublished data), P. falciparum sporozoites were detected more than 4-fold more frequently in An. farauti than P. vivax sporozoites. Furthermore, sporozoite rates across cross-sectional surveys varied widely (being frequently zero). Such disparities may be a consequence of low numbers of sporozoite- positive mosquitoes, sampling bias (including serendipitous placement of collection sites and timing of collection periods), measurement errors (including lack of precision and the inherent large fluctuations in mosquito population densities [41]), as well as the likely predominance of P. vivax infections due to relapses [42]. The difficulty and cost associated with both collecting an adequate number of anophelines of a given species, and their analyses for sporozoites makes the programmatic use of sporozoite rates impractical as a surveillance tool to monitor changes in transmission, particularly in low transmission scenarios. Similarly estimates of survivorship (by parity dissection) are often impractical due to the difficulty in capturing sufficient numbers of mosquitoes to track changes in survivorship in low transmission settings.

By extension, the entomological inoculation rate which is the product of the sporozoite rate and the biting rates would suffer from the same implementation constraints as the sporozoite rate in low transmission settings. Significant associations between EIRs and parasite prevalence in children were established for both P. falciparum and P. vivax in a high transmission area in Papua New Guinea [43] and across multiple sites in Africa for P. falciparum [23]. Analyses of EIRs in Africa questioned the use of the EIR as a means to estimate transmission as reductions in EIRs by 95% in some parts of Africa would be required before an impact on parasite rates in humans would be detectable [23]. Similarly, another study reported that the same human malaria prevalence was associated with a wide range of EIRs [44] while in other studies, high malaria rates (> 44%) had very low EIRs (< 0.001) [23].

Estimating the EIR is technically challenging, labor-intensive and costly, and gives only very imprecise estimates with low external validity [44,45,46]. In low transmission areas, EIR measurements as a surveillance tool are not feasible [44]. The discordant associations between infections in humans and EIRs may also be a function of the lack of standardized methods for estimating the EIR that includes not considering the ecological, demographic, and socioeconomic differences across populations [46].

Villages within (Jack Harbour) and near (New Mala) the high malaria focus had higher densities of the primary vector, An. farauti (26.3 and 1.5 b/p/h-n, respectively), than villages outside the foci (average of 0.2 in Kinamara, Nazareth and Saeragi) suggesting that in areas of low transmission intensity, landing rates may better serve as proxies for inoculation rates (David Smith, personal communication). However, estimating potential transmission intensity by vector landing rates is not without its own challenges as landing rates vary widely in time and space being influenced by numerous factors including weather patterns, larval habitats, adult survivorship and flight patterns as well as variations in the attractiveness to and efficiency of individuals to capture mosquitoes. Using the human landing rate (or other proxies for estimating biting rates, e.g. CDC light traps, double net traps, odour (human or animal) baited traps, as measures of receptivity and potential transmission intensity requires longitudinal and representative selection of multiple collection sites to adequately capture the heterogeneity in vector densities within and among villages and the seasonality that characterises vector populations and transmission in many areas. Such small-scale spatial variations and temporal heterogeneity in mosquito densities can have significant consequences for disease transmission [47] and their characterization needs to be adequately described to be programmatically useful [48].

The disparities documented between relative abundances of anopheline species between larval surveys and human biting rates in this study emphasizes the importance of understanding the behaviours of anophelines and their potential as malaria vectors. In this area, larval surveys were not reliable as predictors of relative receptivity for two reasons. Firstly, there is no known method to translate larval survey data to estimates of adult biting numbers. Secondly, a number of the anopheline species collected were infrequent human biters, especially An. hinesorum. Although An. hinesorum was the most commonly collected anopheline in larval surveys and although a single sporozoite-infected An. hinesorum was identified, this species was a very infrequent human biter and unlikely to be able to maintain malaria transmission.

Conclusions

Despite the heterogeneity amongst and within villages, some common characteristics were found for vector bionomics in a malaria focus in the Solomon Islands. Villages within this malaria focus were consistent in that An. farauti, the dominant vector, was the most common anopheline present and was found in higher densities in villages within and near the foci compared to villages outside of the foci. Furthermore, the vector foci within villages in the malaria focus were larger, encompassing a greater proportion of the villages than villages outside the malaria focus. Thus, villages with both higher malaria receptivity and greater transmission were characterized by the presence of the dominant vector with consistently higher biting densities distributed over a larger area compared to villages of lower receptivity. As such, receptivity and potential transmission risk for programmatic decision-making may be most reliably estimated by the vector biting rate. In such low transmission settings, the other entomological indicators (sporozoite rates, entomological inoculation rates and parity rates) were difficult to estimate with precision due to low numbers of mosquitoes, sampling errors and biases.

Abbreviations

b/p/h-n:

bites/person/half-night

EIR:

entomological inoculation rate

GLM:

generalized linear model

GTS:

Global Technical Strategy

HLC:

human landing catch

ib/p/y:

infective bites/person/year

IRS:

indoor residual spraying

LLINs:

long-lasting insecticidal nets

References

  1. Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–11.

    Article  CAS  Google Scholar 

  2. World Health Organization. Global Technical Strategy for Malaria 2016–2030. Geneva: World Health Organization; 2015.

    Google Scholar 

  3. Cohen J, Moonen B, Snow R, Smith D. How absolute is zero? An evaluation of historical and current definitions of malaria elimination. Malar J. 2010;9:213.

    PubMed  Google Scholar 

  4. Noor AM, Alegana VA, Patil AP, Moloney G, Borle M, Yusuf F, et al. Mapping the receptivity of malaria risk to plan the future of control in Somalia. BMJ Open. 2012;2:e001160.

    Article  Google Scholar 

  5. Reiner RC Jr, Le Menach A, Kunene S, Ntshalintshali N, Hsiang MS, Perkins TA, et al. Mapping residual transmission for malaria elimination. eLife. 2015;4:e09520.

    Article  Google Scholar 

  6. World Health Organization. Malaria surveillance, monitoring & evaluation: a reference manual. Geneva: World Health Organization; 2018.

    Google Scholar 

  7. Russell TL, Beebe NW, Cooper RD, Lobo NF, Burkot TR. Successful malaria elimination strategies require interventions that target changing vector behaviours. Malar J. 2013;12:56.

    Article  Google Scholar 

  8. Taylor B. Changes in the feeding behaviour of a malaria vector, Anopheles farauti Lav., following the use of DDT as a residual spray in houses in the British Solomon Islands Protectorate. Trans R Entomol Soc London. 1975;127:227–92.

    Google Scholar 

  9. Paik Y-H, Avery JG. Problem areas in the malaria eradication programme in the British Solomon Islands. P N G Med J. 1973;17:61–7.

    Google Scholar 

  10. Russell TL, Beebe NW, Bugoro H, Apairamo A, Chow WK, Cooper RD, et al. Frequent blood feeding enables insecticide-treated nets to reduce transmission by mosquitoes that bite predominately outdoors. Malar J. 2016;15:1–9.

    Article  Google Scholar 

  11. World Health Organization. From malaria control to malaria elimination: a manual for elimination scenario planning. Geneva: World Health Organization; 2014.

    Google Scholar 

  12. Cohen J, Smith D, Cotter C, Ward A, Yamey G, Sabot O, Moonen B. Malaria resurgence: a systematic review and assessment of its causes. Malar J. 2012;11:122.

    Article  Google Scholar 

  13. World Health Organization. Information note on the risks associated with the scale back of vector control in areas where transmission has been reduced. WHO/HTM/GMP/2015.7. Geneva: World Health Organization; 2015.

  14. Russell TL, Burkot TR, Bugoro H, Apairamo A, Beebe NW, Chow WK, et al. Larval habitats of the Anopheles farauti and Anopheles lungae complexes in the Solomon Islands. Malar J. 2016;15:1–9.

    Article  Google Scholar 

  15. Russell TL, Beebe NW, Bugoro H, Apairamo A, Collins FH, Cooper RD, et al. Anopheles farauti is a homogeneous population that blood feeds early and outdoors in the Solomon Islands. Malar J. 2016;15:1–7.

    Article  Google Scholar 

  16. Russell TL, Beebe NW, Bugoro H, Apairamo A, Cooper RD, Collins FH, et al. Determinants of host feeding success by Anopheles farauti. Malar J. 2016;15:1–9.

    Article  Google Scholar 

  17. Waltmann A, Darcy AW, Harris I, Koepfli C, Lodo J, Vahi V, et al. High rates of asymptomatic, sub-microscopic Plasmodium vivax infection and disappearing Plasmodium falciparum malaria in an area of low transmission in Solomon Islands. PLoS Negl Trop Dis. 2015;9:e0003758.

    Article  Google Scholar 

  18. Solomon Islands National Statistical Office. Solomon Islands 2009 population and housing census. Honiara, Solomon Islands: Solomon Islands Government; 2009.

    Google Scholar 

  19. Belkin JN. The mosquitoes of the South Pacific (Diptera, Culicidae). Berkeley and Los Angeles: University of California Press; 1962.

    Google Scholar 

  20. World Health Organization. Methods and techniques. WHO Offset Publication No. 13. In: Manual on practical entomology in malaria. Part II. Geneva: WHO Division of Malaria and Other Parasitic Diseases; 1975. p. 195.

    Google Scholar 

  21. Beebe NW, Saul A. Discrimination of all members of the Anopheles punctulatus complex by polymerase chain reaction - restriction fragment length polymorphism analysis. Am J Trop Med Hyg. 1995;53:478–81.

    Article  CAS  Google Scholar 

  22. Echeverry DF, Deason N, Makuru V, Davidson J, Xiao H, Nieldbaski J, et al. A fast and robust single PCR for Plasmodium sporozoite detection in mosquitoes using the cytochrome oxidase I gene. Malar J. 2017;16:230.

  23. Beier JC, Killeen GF, Githure JI. Short report: Entomologic inoculation rates and Plasmodium falciparum malaria prevalence in Africa. Am J Trop Med Hyg. 1999;61:109–13.

    Article  CAS  Google Scholar 

  24. Smith DL, Dushoff J, Snow RW, Hay SI. The entomological inoculation rate and Plasmodium falciparum infection in African children. Nature. 2005;438:492–5.

    Article  CAS  Google Scholar 

  25. Burkot TR, Bugoro H, Apairamo A, Cooper RD, Echeverry DF, Odabasi D, et al. Dataset detailing longitudinal entomological surveys across Western Province, Solomon Islands: James Cook University Tropical Data Hub, 2016. doi: https://doi.org/10.25903/5b722d8d5daf3.

  26. Russell TL, Beebe NW, Bugoro H, Apairamo A, Chow WK, Cooper RD, et al. Dataset describing the biting profile, seasonality and feeding cycle of Anopheles farauti in Haleta village. Solomon Islands: James Cook University Tropical Data Hub. 2018. https://doi.org/10.4225/28/56C671268CF73.

  27. Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46.

    Google Scholar 

  28. Clarke KR. Non-parametric multivariate analysis of changes in community structure. Austral Ecol. 1993;18:117–43.

    Article  Google Scholar 

  29. Core R. Team. R: A language and environment for statistical computing [program]. 3.1.2 version. Vienna, Austria: R Foundation for Statistical. Computing. 2013; http://www.R-project.org.

  30. ArcGIS Desktop. Version 10 [program]. Environmental Systems Research Institute: Redlands, CA; 2011.

    Google Scholar 

  31. FleXScan: Software for the flexible scan statistics [program]. 3.1.2 version. Japan; 2013.

  32. Tango T, Takahashi K. A flexibly shaped spatial scan statistic for detecting clusters. Int J Health Geogr. 2005;4:11.

    Article  Google Scholar 

  33. Tango T, Takahashi K. A flexible spatial scan statistic with a restricted likelihood ratio for detecting disease clusters. Stat Med. 2012;31:4207–18.

    Article  Google Scholar 

  34. World Health Organization. WHO malaria terminology. WHO/HTM/GMP/2016.6. Geneva: World Health Organization; 2016.

  35. Foley DH, Meek SR, Bryan JH. The Anopheles punctulatus group of mosquitoes in the Solomon Islands and Vanuatu surveyed by allozyme electrophoresis. Med Vet Entomol. 1994;8:340–50.

    Article  CAS  Google Scholar 

  36. Beebe NW, Bakotee H, Ellis JT, Cooper RD. Differential ecology of Anopheles puntucaltus and three members of the Anopheles farauti complex of mosquitoes on Guadalcanal, Solomon Islands, identified by PCR-RFLP analysis. Med Vet Entomol. 2000;41:308–12.

    Article  Google Scholar 

  37. Bugoro H, Cooper R, Butafa C, Iro'ofa C, Mackenzie D, Chen C-C, Russell T. Bionomics of the malaria vector Anopheles farauti in Temotu Province, Solomon Islands: issues for malaria elimination. Malar J. 2011;10:133.

    Article  Google Scholar 

  38. Bugoro H, Hii J, Butafa C, Iroofa C, Apairamo A, Cooper R, et al. The bionomics of the malaria vector Anopheles farauti in Northern Guadalcanal, Solomon Islands: issues for successful vector control. Malar J. 2014;13:56.

    Article  Google Scholar 

  39. Bugoro H, Iro'ofa C, Mackenzie D, Apairamo A, Hevalao W, Corcoran S, et al. Changes in vector species composition and current vector biology and behaviour will favour malaria elimination in Santa Isabel Province, Solomon Islands. Malar J. 2011;10:287.

    Article  Google Scholar 

  40. Cooper RD, Waterson DGE, Frances SP, Beebe NW, Pluess B, Sweeney AW. Malaria vectors of Papua New Guinea. Int J Parasitol. 2009;39:1495–501.

    Article  CAS  Google Scholar 

  41. Smith DL, McKenzie FE. Statics and dynamics of malaria infection in Anopheles mosquitoes. Malar J. 2004;3:13.

    Article  Google Scholar 

  42. Robinson LJ, Wampfler R, Betuela I, Karl S, White MT, Li Wai Suen CS, et al. Strategies for understanding and reducing the Plasmodium vivax and Plasmodium ovale hypnozoite reservoir in Papua New Guinean children: a randomised placebo-controlled trial and mathematical model. PLoS Med. 2015;12:e1001891.

    Article  Google Scholar 

  43. Burkot TR, Graves PM, Cattan JA, Wirtz RA, Gibson FD. The efficiency of sporozoite transmission in the human malarias, Plasmodium falciparum and P. vivax. Bull World Health Organ. 1987;65:375–80.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Yukich J, Briët O, Bretscher MT, Bennett A, Lemma S, Berhane Y, et al. Estimating Plasmodium falciparum transmission rates in low-endemic settings using a combination of community prevalence and health facility data. PLoS One. 2012;7:e42861.

    Article  CAS  Google Scholar 

  45. Hay SI, Rogers DJ, Toomer JF, Snow RW. Annual Plasmodium falciparum entomological inoculation rates (EIR) across Africa: literature survey, internet access and review. Trans R Soc Trop Med Hyg. 2000;94:113–27.

    Article  CAS  Google Scholar 

  46. Kelly-Hope L, McKenzie FE. The multiplicity of malaria transmission: a review of entomological inoculation rate measurements and methods across sub-Saharan Africa. Malar J. 2009;8:19.

    Article  Google Scholar 

  47. Smith T, Charlwood JD, Takken W, Tanner M, Spiegelhalter DJ. Mapping the densities of malaria vectors within a single village. Acta Trop. 1995;59:1–18.

    Article  CAS  Google Scholar 

  48. Dye C. Vectorial capacity: must we measure all its components? Parasitol Today. 1986;2:203–9.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Technical assistance when conducting the field collections was provided by Jance Oscar of the Vector Borne Disease Control Programme, Solomon Islands. The authors thank the communities of Western Province for their support as well as the support of Albino Bobogare, Director of National Vector Borne Disease Control Programme, Solomon Islands.

Funding

This work was supported by Grant No. 45114 from the Bill and Melinda Gates Foundation to the Malaria Transmission Consortium. In addition, the support of the National Institute of Allergy and Infectious Diseases of the National Institutes of Health for the International Centers of Excellence in Malaria Research in the Southwest Pacific (subaward U19AI08986 to James Cook University) is gratefully acknowledged. The content is solely the responsibility of the authors and does not necessarily represent the official views of the funders or the Australian Defence Force and/or extant Defence Force Policy.

Availability of data and materials

The datasets supporting the conclusions of this article are available in the JCU Tropical Data Hub repository: Western Province at: https://doi.org/10.25903/5b722d8d5daf3; Central Province at: https://doi.org/10.4225/28/56C671268CF73.

Author information

Authors and Affiliations

Authors

Contributions

Study design, manuscript preparation and data analysis: TRB, RDC, JWK, FHC, NFL and TLR. Data collection: TRB, HB, AA, RDC, NAD, DO, HR, NFL and TLR. Laboratory processing: DFE, NAD, VM, JRD, NWB and NFL. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Tanya L. Russell.

Ethics declarations

Ethics approval and consent to participate

Ethical approvals were obtained from the National Health Research & Ethics Committee, Solomon Islands (2011-05-02, HRE002/16), the James Cook University Human Research Ethics Committee, Australia (H4914 and H6488), the University of Notre Dame Institutional Review Board (FWA 00002462) and University Hospitals Case Medical Centre Institutional Review Board for Human Investigation, USA (05–11–11). Meetings were held with community leaders, study participants and village residents, where the aims, the possible risks and potential benefits of the study were explained in Solomon Islands Pidgin. Mosquito collectors were then recruited from village residents and enrolled in the study after the risks were explained and an informed consent agreement signed.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Burkot, T.R., Bugoro, H., Apairamo, A. et al. Spatial-temporal heterogeneity in malaria receptivity is best estimated by vector biting rates in areas nearing elimination. Parasites Vectors 11, 606 (2018). https://doi.org/10.1186/s13071-018-3201-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-018-3201-1

Keywords