Skip to main content

Spatial–temporal vector abundance and malaria transmission dynamics in Nchelenge and Lake Mweru islands, a region with a high burden of malaria in northern Zambia



Over a decade of vector control by indoor residual spraying (IRS) and long-lasting insecticidal nets (LLINs) distribution on the mainland, and only LLINs on islands had a minimal impact on disease burden in Nchelenge district, northern Zambia. Anopheles funestus and Anopheles gambiae are vectors known only from the mainland. Understanding vector bionomics in the district is necessary for planning and targeting effective vector control. This study aimed to provide information on abundance, seasonality, and Plasmodium falciparum sporozoite infectivity of malaria vectors in Nchelenge, including islands.


Mosquitoes were collected in 192 CDC indoor light traps set in 56 households between January 2015 and January 2016. Morphological and molecular species identifications and P. falciparum circumsporoites by ELISA were performed. Mosquito counts and relative abundances from the islands and mainland were compared, and household factors associated with vector counts were determined.


A total of 5888 anophelines were collected during the study. Of these, 5,704 were female Anopheles funestus sensu lato (s.l.) and 248 female An. gambiae s.l. The highest proportion of An. funestus (n = 4090) was from Chisenga Island and An. gambiae (n = 174) was from Kilwa Island. The highest estimated counts per trap for An. funestus s.l. were from Chisenga island, (89.9, p < 0.001) and from the dry season (78.6, p < 001). For An. gambiae the highest counts per trap were from Kilwa island (3.1, p < 0.001) and the rainy season (2.5, p = 0.007). The highest estimated annual entomological inoculation rate was from Chisenga Island with 91.62 ib/p/y followed by Kilwa Island with 29.77 ib/p/yr, and then Mainland with 19.97 ib/p/yr.


There was varied species abundance and malaria transmission risk across sites and seasons. The risk of malaria transmission was perennial and higher on the islands. The minimal impact of vector control efforts on the mainland was evident, but limited overall. Vector control intervention coverage with effective tools needs to be extended to the islands to effectively control malaria transmission in Nchelenge district.


The use of long-lasting insecticidal nets (LLINs), indoor residual spraying (IRS), artemisinin-based combination therapy (ACT), and intermittent preventive therapy in pregnancy (IPTp) in the last decade has led to successes against malaria in sub-Saharan Africa [1,2,3,4,5]. Sadly, not all areas receiving control interventions face an easy transition from malaria endemicity to malaria elimination [6]. For some high-transmission areas, the package of control interventions may need to be prescribed differently and specifically. In Zambia, Nchelenge District in the north is an area with high malaria transmission, with a sustained malaria prevalence between 30 and 50%, despite incremental coverage in control interventions comprised of ACTs, IRS, LLINs, IPTp since 2008, and more recently, integrated community case management (ICCM) [7, 8]. Studies from mainland Nchelenge have reported complex vector bionomics of two malaria vector species, Anopheles funestus sensu stricto (s.s.) and An. gambiae s.s. These studies have established that An. funestus is widely distributed and has a population peak during the dry season, while counts of An. gambiae remain relatively and consistently low throughout the year and are more abundant near Lake Mweru on the western border of the district [9,10,11,12]. The few genomics studies completed in Nchelenge have revealed the sympatric existence of Clades I and II of An. funestus s.s. and only S-form An. gambiae, with populations exhibiting varied insecticide susceptibility to carbamates, pyrethroids, and DDT [13,14,15]. It is unknown, however, to what extent these findings can be extended to the island settings within the district, where control measures are much more challenging to implement. The two big islands of the district, Kilwa Island and Chisenga Island, comprise nearly 13% of the district’s total population and, at the time of the study, had received only LLINs as malaria control interventions.

During the study period, cases of malaria reported between January and December 2015 for all age groups from Chisenga and Kilwa Islands rural health centres were 4, 282 and 8, 439, respectively (Nchelenge DHIS). These represented cases per 1000 people of 375 and 746 for Chisenga and Kilwa islands, respectively. Nchelenge Health Centre on the mainland reported 9360 cases with an incidence rate of 476/1000 (Nchelenge DHIS). For each health facility, almost half of the cases (2353 Chisenga Island, 4061 Kilwa Island, and 4477 Nchelenge) were in children 0–5 years old.

Given the enormity of malaria incidence in Nchelenge District, critical data is needed to better understand the underlying factors leading to the lack of success of malaria control in Nchelenge, especially in the context of malaria elimination [16]. Data on the two primary malaria vector species from the islands of Nchelenge is incredibly sparse. Published data was only from Kilwa Island mosquitoes, dating back to a 1926 report on non-malaria vectors [15]. The other known data on vectors came from a recent publication on An. funestus s.s. complete mitogenome [17], which included specimens collected from Kilwa Island during the initial preparatory stages of this study. This study aimed to expand ongoing malaria entomological studies in Nchelenge District to include the main islands in Lake Mweru. The focus was on vector seasonal and spatial abundance, Plasmodium falciparum sporozoite infection rates, and how these data may impact current and future malaria control strategies in the district. The study period was January 2015 to January 2016 as part of the Southern and Central Africa International Centres of Excellence for Malaria Research (ICEMR), a research collaboration set up to study malaria transmission and the impact of control interventions in Zambia, Zimbabwe, and the Democratic Republic of the Congo (DRC).


Study area

Nchelenge District is located in northern Zambia in Luapula Province (Fig. 1) at the mouth of the Luapula River, forming part of the Luapula-Mweru delta. Lake Mweru and the Luapula River form the border with the DRC, and the Luapula-Mweru system interface includes vast swamps and lagoons and a number of islands. Prominent among the islands on the Zambian side of Lake Mweru are Chisenga, Kilwa, and Isokwe (Fig. 1). The mainland along Lake Mweru hosts a network of shallow streams with slow-moving waters and expansive vegetation cover. Mosquitoes are hypothesised to breed along the edges of these streams, which remain inundated for extended periods of the year.

Fig. 1
figure 1

Map of Nchelenge and sampled households across the District (map composed using QGIS 3.22.7)

Nchelenge is largely a rural district with an estimated human population at the time of the survey of 147,927 [18]. The mainland constituted about 77% of the total district population. The island populations by head count from the district health office were estimated at 11409 for Chisenga Island, 3501 for Isokwe Island, and 11305 for Kilwa Island (Nchelenge DHO). The majority of the island's inhabitants are involved in fishing, while also engaging in subsistence farming. The formal housing structures with metal roofing are typical of the district's urban areas. The rural areas comprise houses with mud walls and thatch roofing. Makeshift houses, which expose the inhabitants to mosquito bites, are characteristic of farming fields on the mainland and nearly all houses on Isokwe Island and in fishing camps on the islands. Seasonal human movement between the interior mainland, lakeside, and islands is common. On the mainland, the beginning of the rainy season necessitates the shift to the interior, where there are farming fields. The rainy season is usually between November and May, and the dry season comes between June and September. The collections conducted during the rainy season in 2015 were defined as rainy season 1, and those done in 2016 were designated rainy season 2.

Vector control

Indoor residual spraying and bed net distribution form the main vector control interventions conducted in Nchelenge by the national programme. The IRS programme is normally conducted between September and December of each year. During the study period, IRS with pirimiphos-methyl (Actellic® 300CS) was conducted in 2014 and 2015, and a mass campaign for LLIN (Olyset®) distribution took place in October 2014.

Selection of study households

Satellite imagery of the three islands from Google Maps and an ARCGIS image of mainland Nchelenge District were overlaid with 1 km2 grids [19]. Grids with human habitation were enumerated, and study households were selected by simple random sampling. On the ground, houses within a selected grid (hereafter referred to as a village) were enumerated, and four houses were randomly selected. A total of 56 houses across 14 villages in the study area were subsequently selected and their geographic coordinates recorded (using a Garmin OREGON® 450 GPS gadget). The 14 villages were distributed as follows: Chisenga Island, 4; Isokwe Island, 2; Kilwa Island, 5; and Mainland, 3.

Mosquito collection

The collection of mosquitoes was conducted using CDC light traps hung indoors approximately 150 cm above the floor at the foot end of the sleeping space. A total of 192 traps were set indoors during the entire study period. A total of five rounds of collection were conducted, as follows: First rainy season: 32 households (HH) in January 2015 and 36 HH in March 2015. During the dry season, 37 HH were recorded in July 2015 and 42 HH in October 2015. Finally, the second rainy season 44 HH in January 2016 Each of the selected households was sampled a minimum of one and a maximum of five times (Table 1) over the entire sampling period. Only data from 42 households, a total of 165 trap nights (marked with § in Table 1) sampled 3–5 times, were included in the analysis of malaria vectors’ spatial and temporal abundances.

Table 1 Distribution of households and sampling efforts across sites over the course of the study

Processing of collected mosquitoes

All anopheline mosquitoes collected were morphologically identified by standard identification keys to species groups upon collection in the field [20, 21]. These were individually packed into 0.5-mL centrifuge tubes with desiccant for further laboratory processing. A representative sub-sample of 20% of specimens morphologically identified as female An. funestus sensu lato (s.l.) and An. gambiae s.l. were processed for polymerase chain reaction (PCR) and circumsporozoite ELISA.

Mosquitoes were dissected at the thoracic and abdominal joints into the head, thorax, and abdomen. DNA was extracted from either a leg or abdomen and processed by PCR for species identification [22, 23]. The M and S molecular form determinations of An. gambiae s.s. were completed by PCR, as previously described [24].

TaqMan PCR, as previously described [25], identified the An. funestus clades. The P. falciparum sporozoite infections were determined by the circumsporozoite ELISA assay pf-210 [26, 27]. To remove false positives, the homogenates of all the assays that tested positive on the first run of circumsporozoite ELISA were boiled at 100 °C for ten minutes [28], and the ELISA was repeated. The P. falciparum sporozoite rate was calculated as the number of mosquitoes detected with P. falciparum sporozoites out of the total tested by circumsporozoite ELISA [29]. The biting rate was taken as previously done by Das et al. [10], taking into account the total number of people in the households. Annual entomological inoculation rates were estimated as the product of the biting rate and sporozoite rate multiplied by 365 days.

Data analysis

The mosquito counts were entered in Excel sheets and exported to R and STATA 15 (Release 15, College Station, TX: Stata Corp LLC) for analysis. Morphological data were used for the analyses of species abundance and distribution. Only data from households sampled at least three times were included in the analyses of the seasonal and spatial abundance of vector species. The species PCR data were used to compare vector species infection rates and their spatial and seasonal occurrence.

Variations in total female anophelines collected from each site were evaluated separately for each of the two vector species by multilevel mixed-effects negative binomial regression with household as a random covariate. The covariates of fixed effects were season and site of collection. Marginal effects and means of vector counts were predicted for the season and site of collection. The P. falciparum sporozoite infection rates in Anopheles species were seasonally and spatially compared and predicted by logistic regression analysis. QGIS version 3.10.2 was used to create maps and plot proportions of mosquito counts as well as P. falciparum infection rates.


Proportions of major vectors in the study area

A total of 5888 Anopheles mosquitoes were collected indoors from 56 households in 14 villages (Fig.2). A total of 858 people slept in the houses on the nights of the mosquito sample collection. Of the total collection, 5579 were female anophelines, distributed as follows: Anopheles funestus s.l. = 5323, An. gambiae s.l. = 249, Anopheles coustani s.l. = 4, and unidentified = 3.

Fig. 2
figure 2

Mean composition of the primary malaria vectors An. funestus s.l. and An. gambiae s.l. across the study area (map composed using QGIS3.22.7)

The count data of the two major malaria vectors were highly skewed, with 31 traps recording zero An. funestus s.l. collections and 97 traps recording zero An. gambiae s.l. collections. The maximum count of An. funestus s.l. in a single light trap was 497 from Chisenga Island, and the maximum for An. gambiae s.l. was 30 from Kilwa Island (Table 2).

Table 2 Proportions of the total vectors collected and the maximum count in a single light trap in different sites

Anopheles funestus s.l. count (n = 5240) were dominant over An. gambiae s.l. (n = 238) across all study sites (Table 2). Chisenga Island had the highest proportion of An. funestus s.l. (n = 4090), and Kilwa Island had the highest proportion of An. gambiae s.l. (n = 174).

Factors associated with malaria vector counts

The estimated counts per trap for the two vector species across the study sites ranged from 8.602 to 89.97 for An. funestus s.l. and 0.499 to 3.083 for An. gambiae s.l. (Table 3). On the other hand, seasonal counts per trap for these vectors ranged from 4.869 to 78.606 for An. funestus and 0.818 to 2.548 for An. gambiae (Table 3).

Table 3 Factors associated with the malaria vector counts in CDC light traps set indoors

The An. funestus counts were strongly associated with Chisenga Island (IRR = 9.341, p < 0.001), and those of An. gambiae were so with Kilwa Island (IRR = 4.573, p = 0.004). There was a significant reduction in the estimated An. funestus count during rainy seasons 1 and 2, IRR = 0.189, p < 0.001, and IRR = 0.069, p < 0.001, respectively (Table 3). The An. gambiae count reduced during rainy season 1 (IRR = 0.770, p = 0.494), but the count increased during rainy season 2 (IRR = 1.868, p = 0.115) (Table 3).

Molecular identifications

Molecular species identifications by PCR on a sub-sample (1108 females) successfully amplified 978 (88%) of the specimens processed. DNA from 131 specimens did not amplify.

Identification of species by PCR showed An. funestus s.s. (n = 773) and Anopheles leesoni (n = 10) as members of the An. funestus group. Further mitochondrial DNA analysis of An. funestus s.s. revealed the existence of Clade I (n = 736) and Clade II (n = 37), as previously reported [13]. The An. gambiae s.l. PCR analysis of specimens revealed only An. gambiae s.s. (n = 195). Further processing by TaqMan hydrolysis PCR to determine the M-form (Anopheles coluzzii) and S-forms revealed the S-form of An. gambiae s.s. as the only member of the complex collected.

Plasmodium falciparum infection rates in the vectors

A total of 30 specimens gave a positive P. falciparum sporozoite ELISA result on the first run. When the homogenates were heated to remove false positives, 25 still came out positive. The P. falciparum circumsporozoite ELISA showed positive results in An. funestus s.s. Clade I, 18 (n = 736), and An. gambiae, 7 (n = 195). No sporozoites were detected in An. funestus s.s. Clade II (n = 37) or An. leesoni (n = 10).

In the study sites, P. falciparum infections were found in both primary vectors, An. funestus s.s. and An. gambiae s.s., from Kilwa Island and the mainland. For Isokwe Island and Chisenga Island, only An. funestus s.s. was found infected with P. falciparum. Using the foraging rates based on the count of mosquitoes and household occupants the annual entomological inoculation rates (EIR) were estimated for each collection site. The EIR for Chisenga Island was 91.62 infectious bites per person per year (ib/p/yr), all of which came from An. funestus s.s. Kilwa Island's EIR was 29.77 ib/p/yr, with species contributions of 21.21 ib/p/yr and 8.57 ib/p/yr from An. funestus s.s. and An. gambiae s.s., respectively. Mainland EIR was estimated to be 19.97 ib/p/yr, with An. funestus s.s. contributing 18.12 ib/p/yr and An. gambiae s.s. giving 1.85 ib/p/yr.

Seasonally, P. falciparum sporozoites-infected vectors were detected in both dry and rainy seasons (Fig. 3). In the rainy season, P. falciparum sporozoites were detected among An. funestus s.s. samples from Chisenga Island (Fig. 3A) and in An. gambiae s.s. samples from Kilwa Island and the mainland (Fig. 3C). In the dry season, P. falciparum sporozoite infections were detected in An. funestus s.s. from Isokwe Island, Kilwa Island, and the mainland (Fig. 3B). No dry-season P. falciparum sporozoites were detected in An. gambiae s.s. from any site (Figure 3D).

Fig. 3
figure 3

Spatial and seasonal distribution of P. falciparum sporozoite infections in the primary vectors across study sites. A and B represent infections in An. funestus s.l. in the rainy and dry seasons, respectively. C and D represent infections in An. gambiae s.l. rainy and dry seasons, respectively (maps composed using QGIS 3.22.7)

Seasonal and spatial comparisons in P. falciparum infection rates by logistic regression (Table 4) in the two vector species revealed 54% lower odds of infection in An. gambiae s.s. compared to An. funestus s.s. Clade I (OR = 0.460, 95% CI: 0.184–1.150, p = 0.097). The odds of infection with P. falciparum in vectors were lower on Chisenga Island compared to the mainland (OR = 0.646, 95% CI: 0.198–2.109, p = 0.649). The odds of infection in vectors were higher on Isokwe Island compared to the mainland (OR = 3.564, 95% CI: 0.506–25.100, p = 0.202), and higher on Kilwa Island compared to the mainland (OR = 1.322, 95% CI: 0.388–4.506, p = 0.655).

Table 4 Plasmodium falciparum infection analysis by logistic regression

The estimated probability of An. funestus s.s. Clade I infection with P. falciparum sporozoites during the study was 0.042 (95% CI: 0.16-0.68, p = 0.001). The estimated probability of infection in An. gambiae s.s. was 0.020 (95% CI: 0.004–0.036, p = 0.014). Seasonal predicted probabilities of infection with P. falciparum in vectors were: dry season 0.021 (95% CI: 0.008–0.035, p = 0.002), rainy season 1, 0.015 (95% CI: 0.000–0.034, p = 0.092), and rainy season 2, 0.073 (95% CI: 0.026–0.119, p = 0.002).


This study determined the spatial and seasonal abundances of the two primary vectors of malaria in Nchelenge, including islands. The infection with P. falciparum sporozoites in the two primary vectors of malaria, An. funestus s.s. and An. gambiae s.s., in Nchelenge was determined. There were significant differences in the abundance of An. funestus s.l. on Chisenga Island and of An. gambiae s.l. on Kilwa Island compared to the mainland. The An. funestus s.l. was more abundant during the dry season, and the An. gambiae was abundant during the rainy season. Predicted probabilities for P. falciparum infection were significant for both An. gambiae s.s. and An. funestus s.s., with An. funestus s.s. showing a slight edge over An. gambiae s.s. When the two vectors were considered together, the predicted probability of a malaria vector for infection with P. falciparum during both the dry and rainy seasons was significant. The estimated annual entomological inoculation rate, which is a measure of human exposure to the risk of malaria transmission, was highest on Chisenga Island with 91.62 ib/p/yr, followed by Kilwa Island with 29.77 ib/p/y, and lastly Mainland with 19.97 ib/p/yr.

These results show a dynamic entomological picture driving malaria transmission in the district. The results obtained here both confirmed what was already known about Nchelenge and provided new insights about the malaria vectors on the islands. At the time of this study, the only form of vector control for the islands was net distribution. Anopheles funestus s.l. dominated and was more abundant during the dry season, in agreement with previous findings [10, 12]. This study went further to show that Chisenga Island had disproportionally higher densities of An. funestus s.l. compared to either Kilwa Island or the mainland. The highest counts of An. gambiae s.l. were recorded from Kilwa Island. These findings seem to reflect local differences in the ecologies of these different study sites, which may be important for targeting control interventions [30, 31].

The findings of this study are in support of prior studies that assessed household risk factors for mosquito abundance on the mainland [9, 12]. Mosquito counts found in this study suggest that high bed net usage could be driven by high indoor densities of pyrethroid-resistant mosquitoes [9]. The findings from this study indicate that this scenario would be widespread across Nchelenge, including islands. These scenarios would raise concern as to the efficacy of the use of LLINs as the only vector control intervention for the islands in Nchelenge District. Other interventions, including the use of PBO-synergized LLINs designed to address the high rates of pyrethroid resistance, need to be included. Other household factors, such as thatch roofing and open eaves indicative of poor construction, allowed easier mosquito entry. Improved household structures should be considered a compliment to or integrated into the vector control programme for the district. Studies done in Cameroon, The Gambia, and Kenya have demonstrated that simple and low-cost improvements to house structure would be effective at preventing mosquito entry and reducing indoor densities [32,33,34,35].

The parasite sporozoite rates observed suggested variations in spatial intensity but a sustained, perennial risk of vectors infected with P. falciparum sporozoites. Anopheles funestus s.s. contributed to transmission at all study sites during the dry season and on Chisenga Island during the rainy season. In contrast, An. gambiae appeared to drive transmissions during the rainy season on Kilwa Island and the mainland. The mainland and Kilwa Island were the only sites where both vector species were found infected with P. falciparum. Despite their relatively low numbers, An. gambiae had comparable rates of P. falciparum sporozoite infections with the more abundant An. funestus s.s. on both Kilwa Island and the mainland. This demonstrated that An. gambiae made substantial contribution to malaria transmission in Nchelenge, a fact that needs to be considered when designing and implementing vector control interventions. This study found no P. falciparum-infected Clade II An. funestus s.s. This was in contrast with the prior study by Choi et al. [13] on the mainland that could not attribute any difference in malaria transmission to either Clade, and thus our findings could be attributed to sample size.

Two previous studies, one intensive and another extensive, both conducted in 2012 and 2013, made comparable estimates of the EIRs for the Nchelenge mainland with different calculations of the human biting rate [10, 12]. The EIR estimated for the mainland in the current study was much lower for both An. funestus s.s. and An. gambiae s.s. The EIR is an important indicator of the intensity of malaria transmission in an area. A reduction in EIR is one of the indices used to measure the impact of vector control activity [36, 37]. This study followed a mass LLIN distribution campaign and two years of IRS with pirimiphos-methyl on the mainland. Two studies on the mainland that looked at the effects of three years of IRS with pirimiphos-methyl on parasite prevalence and the number of malaria vectors found that the number of vectors in sprayed homes and targeted areas went down, but parasite prevalence did not change much [38, 39]. The present study findings on the mainland suggest reduced EIR estimates, which seem to agree with the reduced vector abundances and parasite prevalence reported during the study period.

Twelve months of concomitant malaria incidence data were obtained from health centres (HCs) servicing the three individual collection study sites. According to these data, malaria incidence was highest on Kilwa Island HC, then on the mainland Nchelenge HC, and lowest on Chisenga Island HC. In contrast, Chisenga Island had the highest estimated EIR, followed by Kilwa Island, and then the mainland. In an area with high vector counts, the discrepancy between estimated EIR and reported malaria cases may be due to the unreported cases that get treated at home within communities. No prevalence studies were conducted in the present study, nor had they been conducted before on the islands. For the islands in Nchelenge, local factors including ecology, demographics, malaria control use, and mosquito bionomics would need to be considered in the estimation and interpretation of the EIR. Monitoring how EIRs would be affected would be important in the evaluation of the scale-up of the malaria control package on the islands of Nchelenge.


Consistent with prior studies in Nchelenge District, this study reports a high abundance of the primary vector, An. funestus s.s., and a relatively low abundance of An. gambiae, including island sites for the first time. In sites and periods where both vector species were involved in malaria transmission, sporozoite infection rates were comparable, thereby adding to a higher risk for malaria transmission. The combination of IRS and LLINs on the mainland might have a very minimal impact on malaria transmission. Sustained vector control at all the sites should be emphasised with expanded coverage of intervention on the islands. Non-chemical mosquito bite control and prevention options should be considered for Nchelenge. Continued research on the vectors, the parasites, and the human host at different sites is required to fully understand the burden and effectively target control in different areas of Nchelenge.

Availability of data and materials

Raw data on the mosquito collections and entomological indices is available upon reasonable request from the corresponding author.



Artemisinin-based combination therapy


Aeronautical Reconnaissance Coverage Geographic Information System


Centres for Disease Control and prevention




District Health Information System


District Health Office


Deoxyribonucleic Acid


Democratic Republic of Congo


Entomological Inoculation Rate


Enzyme-linked Immunosorbent Assay


Ethics Review Committee


Global Position System


Health Centre




Integrated Community Case Management


International Centres of Excellence for Malaria Research


Intermittent Preventive Therapy in Pregnancy


Indoor Residual Spraying


Insecticide-Treated Nets


Long-lasting Insecticidal Nets

M & S:

Mopti & Savanna strain


Williams Mean


Odds ratio


Polymerase chain reaction


Quantum Geographic Information System


Scientific and Technical Committee


Tropical Diseases Research Centre


  1. Sinka ME, Bangs MJ, Manguin S, Rubio-Palis Y, Chareonviriyaphap T, Coetzee M, et al. A global map of dominant malaria vectors. Parasit Vectors. 2012;5:69.

    Article  PubMed  PubMed Central  Google Scholar 

  2. WHO. World Malaria Report 2020. Geneva, World Health Organization. Accessed 7 Apr 2021.

  3. WHO. The “World Malaria Report 2019” at a glance. Geneva, World Health Organization. Accessed 7 Apr 2021.

  4. Maharaj R, Morris N, Seocharan I, Kruger P, Moonasar D, Mabuza A, et al. The feasibility of malaria elimination in South Africa. Malar J. 2012;11:423.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Moonasar D, Maharaj R, Kunene S, Candrinho B, Saute F, Ntshalintshali N, et al. Towards malaria elimination in the MOSASWA (Mozambique, South Africa and Swaziland) region. Malar J. 2016;15:419.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Dhiman S. Are malaria elimination efforts on right track? An analysis of gains achieved and challenges ahead. Infect Dis Poverty. 2019;8:14.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Mukonka VM, Chanda E, Haque U, Kamuliwo M, Mushinge G, Chileshe J, et al. High burden of malaria following scale-up of control interventions in Nchelenge District, Luapula Province, Zambia. Malar J. 2014;13:153.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Moss WJ, Dorsey G, Mueller I, Laufer MK, Krogstad DJ, Vinetz JM, et al. Malaria epidemiology and control within the international centers of excellence for malaria research. Am J Trop Med Hyg. 2015.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Hast MA, Stevenson JC, Muleba M, Chaponda M, Kabuya J-B, Mulenga M, et al. Risk factors for household vector abundance using indoor CDC light traps in a high malaria transmission area of Northern Zambia. Am J Trop Med Hyg. 2019;101:126–36.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Das S, Muleba M, Stevenson JC, Norris DE. Habitat partitioning of malaria vectors in Nchelenge district, Zambia. Am J Trop Med Hyg. 2016;94:1234–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Das S, Muleba M, Stevenson JC, Pringle JC, Norris DE. Beyond the entomological inoculation rate: Characterising multiple blood feeding behaviour and Plasmodium falciparum multiplicity of infection in Anopheles mosquitoes in northern Zambia. Parasit Vectors. 2017;10:45.

  12. Stevenson JC, Pinchoff J, Muleba M, Lupiya J, Chilusu H, Mwelwa I, et al. Spatio-temporal heterogeneity of malaria vectors in northern Zambia: implications for vector control. Parasit Vectors. 2016;9:510.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Choi KS, Christian R, Nardini L, Wood OR, Agubuzo E, Muleba M, et al. Insecticide resistance and role in malaria transmission of Anopheles funestus populations from Zambia and Zimbabwe. Parasit Vectors. 2014;7:464.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Venter N, Oliver SV, Muleba M, Davies C, Hunt RH, Koekemoer LL, et al. Benchmarking insecticide resistance intensity bioassays for Anopheles malaria vector species against resistance phenotypes of known epidemiological significance. Parasit Vectors. 2017;10:88.

    Article  Google Scholar 

  15. Schwetz J. Contribution à l’étude des moustiques d’Elisabethville et de quelques autres localités du Katanga (Congo belge). Bull Soc Path Exot. 1927;20:170–92.

    Google Scholar 

  16. ZNMEC. National Malaria Elimination Strategic Plan 2015–2020. Zambia National Malaria Elimination Centre, 2015.

  17. Jones CM, Lee Y, Kitchen A, Collier T, Pringle JC, Muleba M, et al. Complete Anopheles funestus mitogenomes reveal an ancient history of mitochondrial lineages and their distribution in southern and central Africa. Sci Rep. 2018;8:9054.

    Article  PubMed  PubMed Central  Google Scholar 

  18. CSO. Population Census Zambia—preliminary report. Lusaka: Central Statistics Office; 2011.

  19. Keating J, Macintyre K, Mbogo C, Githeko A, Regens JL, Swalm C, et al. A geographic sampling strategy for studying relationships between human activity and malaria vectors in urban Africa. Am J Trop Med Hyg. 2003;68:357–65.

    Article  PubMed  Google Scholar 

  20. Gillies MT, Coetzee M. A supplement to the Anophelinae of the South of the Sahara (Afrotropical Region). Publ S Afr Inst Med Res. 1987;55:1–143.

    Google Scholar 

  21. Gillies M, De Meillon B. The Anophelinae of Africa south of the Sahara (Ethiopian zoogeographical region). Publications of the South African Institute for Medical Research. 1968;54:1–343.

  22. Koekemoer LL, Kamau L, Hunt RH, Coetzee M. A cocktail polymerase chain reaction assay to identify members of the Anopheles funestus (Diptera: Culicidae) group. Am J Trop Med Hyg. 2002;6:804–11.

    Article  Google Scholar 

  23. Scott JA, Brogdon WG, Collins FH. Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993;49:520–9.

    Article  CAS  PubMed  Google Scholar 

  24. Favia G, Torre A, Bagayoko M, Lanfrancotti A, Sagnon NF, Toure YT. Molecular identification of sympatric chromosomal forms of Anopheles gambiae and further evidence of their reproductive isolation. Insect Mol Biol. 1997;6:377–83.

    Article  CAS  PubMed  Google Scholar 

  25. Choi KS, Coetzee M, Koekemoer LL. Detection of clade types (clades I and II) within Anopheles funestus sensu stricto by the hydrolysis probe analysis (TaqMan assay). Parasit Vectors. 2013;6:173.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Wirtz RA, Burkot TR, Graves PM, Andre RG. Field evaluation of enzyme-linked immunosorbent assays for Plasmodium falciparum and Plasmodium vivax sporozoites in mosquitoes (Diptera: Culicidae) from Papua New Guinea. J Med Entomol. 1987;24:433–7.

    Article  CAS  PubMed  Google Scholar 

  27. Wirtz RA, Duncan JF, Njelesani EK, Brown AE, Oster CN, Were JBO, et al. ELISA method for detecting Plasmodium falciparum circumsporozoite antibody. Bull World Health Organ. 1989;67:535–42.

  28. Durnez L, Van Bortel W, Denis L, Roelants P, Veracx A, Trung HD, et al. False positive circumsporozoite protein ELISA: a challenge for the estimation of the entomological inoculation rate of malaria and for vector incrimination. Malar J. 2011;10:195.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Mwangangi JM, Muturi EJ, Muriu SM, Nzovu J, Midega JT, Mbogo C. The role of Anopheles arabiensis and Anopheles coustani in indoor and outdoor malaria transmission in Taveta District, Kenya. Parasit Vectors. 2013;6:114.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Nkya TE, Fillinger U, Sangoro OP, Marubu R, Chanda E, Mutero CM. Six decades of malaria vector control in southern Africa: a review of the entomological evidence-base. Malar J. 2022;21:279.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Epopa PS, Collins CM, North A, Millogo AA, Benedict MQ, Tripet F, et al. Seasonal malaria vector and transmission dynamics in western Burkina Faso. Malar J. 2019;18:113.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Atieli H, Menya D, Githeko A, Scott T. House design modifications reduce indoor resting malaria vector densities in rice irrigation scheme area in western Kenya. Malar J. 2009;8:108.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Ogoma SB, Lweitoijera DW, Ngonyani H, Furer B, Russell TL, Mukabana WR, et al. Screening mosquito house entry points as a potential method for integrated control of endophagic filariasis, arbovirus and malaria vectors. PLoS Negl Trop Dis. 2010;4: e773.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Lindsay SW, Jawara M, Paine K, Pinder M, Walraven GEL, Emerson PM. Changes in house design reduce exposure to malaria mosquitoes. Trop Med Int Health. 2003;8:512–7.

    Article  CAS  PubMed  Google Scholar 

  35. Jatta E, Jawara M, Bradley J, Jeffries D, Kandeh B, Knudsen JB, et al. How house design affects malaria mosquito density, temperature, and relative humidity: an experimental study in rural Gambia. Lancet Planet Health. 2018;2:e498-508.

    Article  PubMed  Google Scholar 

  36. Coleman S, Dadzie SK, Seyoum A, Yihdego Y, Mumba P, Dengela D, et al. A reduction in malaria transmission intensity in Northern Ghana after 7 years of indoor residual spraying. Malar J. 1971;16:324.

    Article  Google Scholar 

  37. Salako AS, Dagnon F, Sovi A, Padonou GG, Aïkpon R, Ahogni I, et al. Efficacy of Actellic 300 CS-based indoor residual spraying on key entomological indicators of malaria transmission in Alibori and Donga, two regions of northern Benin. Parasit Vectors. 2019;12:612.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Hast MA, Chaponda M, Muleba M, Kabuya J-B, Lupiya J, Kobayashi T, et al. The impact of 3 years of targeted indoor residual spraying with pirimiphos-methyl on malaria parasite prevalence in a high-transmission area of northern Zambia. Am J Epidemiol. 2019;188:2120–30.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Hast MA, Stevenson JC, Muleba M, Chaponda M, Kabuya JB, Mulenga M, et al. The impact of three years of targeted indoor residual spraying with pirimiphos-methyl on household vector abundance in a high malaria transmission area of northern Zambia. Am J Trop Med Hyg. 2021;104:683–94.

    Article  CAS  Google Scholar 

Download references


We would like to thank the following for their contributions in the mapping, field work, laboratory processing of specimens, DNA sequence and data analyses: Aniset Kamanga, Timothy Shields, James Sichivula, Sebastian Mulenga, David Nkomba, Nelius Venter, Jay Sikalima, Ian Mwelwa, Hunter G. Chilusu, Limonty Simubali and Sophie Berube. We are grateful to the Ministry of Health and the Nchelenge District Health Office for their support. We appreciate the help received through the office of the District Commissioner and the local traditional authority. Special thanks go to the Nchelenge community and to the households that allowed the research team into their houses.


The South Africa-Zambia joint research grant through the Zambian National Science and Technology Council and the South African National Research Foundation as well as a research grant from the NIH funding No. 3U19AI089680 through the International Centres of Excellence for Malaria Research funded this work.

Author information

Authors and Affiliations



MM designed the study, collected and analysed the data, and drafted the first manuscript. KJM, JS, and DN provided guidance on study design and data interpretation. All authors reviewed and approved the manuscript.

Corresponding author

Correspondence to Mbanga Muleba.

Ethics declarations

Ethics approval and consent to participate

This study did not use human specimens, and written consent for indoor collection of mosquitoes was obtained from household heads. This study was approved and cleared by the Tropical Diseases Research Centre, Ndola, Zambia (TDRC, ERC, STC/2015/8).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

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 licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Muleba, M., Mbata, K.J., Stevenson, J.C. et al. Spatial–temporal vector abundance and malaria transmission dynamics in Nchelenge and Lake Mweru islands, a region with a high burden of malaria in northern Zambia. Malar J 22, 327 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: