Skip to main content

Improving the efficiency of reactive case detection for malaria elimination in southern Zambia: a cross-sectional study

Abstract

Background

Reactive case detection (RCD) seeks to enhance malaria surveillance and control by identifying and treating parasitaemic individuals residing near index cases. In Zambia, this strategy starts with passive detection of symptomatic incident malaria cases at local health facilities or by community health workers, with subsequent home visits to screen-and-treat residents in the index case and neighbouring (secondary) households within a 140-m radius using rapid diagnostic tests (RDTs). However, a small circular radius may not be the most efficient strategy to identify parasitaemic individuals in low-endemic areas with hotspots of malaria transmission. To evaluate if RCD efficiency could be improved by increasing the probability of identifying parasitaemic residents, environmental risk factors and a larger screening radius (250 m) were assessed in a region of low malaria endemicity.

Methods

Between January 12, 2015 and July 26, 2017, 4170 individuals residing in 158 index and 531 secondary households were enrolled and completed a baseline questionnaire in the catchment area of Macha Hospital in Choma District, Southern Province, Zambia. Plasmodium falciparum prevalence was measured using PfHRP2 RDTs and quantitative PCR (qPCR). A Quickbird™ high-resolution satellite image of the catchment area was used to create environmental risk factors in ArcGIS, and generalized estimating equations were used to evaluate associations between risk factors and secondary households with parasitaemic individuals.

Results

The parasite prevalence in secondary (non-index case) households was 0.7% by RDT and 1.8% by qPCR. Overall, 8.5% (n = 45) of secondary households had at least one resident with parasitaemia by qPCR or RDT. The risk of a secondary household having a parasitaemic resident was significantly increased in proximity to higher order streams and marginally with increasing distance from index households. The adjusted OR for proximity to third- and fifth-order streams were 2.97 (95% CI 1.04–8.42) and 2.30 (95% CI 1.04–5.09), respectively, and that for distance to index households for each 50 m was 1.24 (95% CI 0.98–1.58).

Conclusion

Applying proximity to streams as a screening tool, 16% (n = 3) more malaria-positive secondary households were identified compared to using a 140-m circular screening radius. This analysis highlights the potential use of environmental risk factors as a screening strategy to increase RCD efficiency.

Background

In countries where a state of low-endemicity for malaria has been established and maintained, strategies and policies geared toward elimination of both vector and parasites have begun to take form. An essential component that has been incorporated in many of these elimination strategies is malaria control using environmental risk factors [1]. These strategies take advantage of the spatial distribution of malaria, which varies depending on the ecology and population, but in regions with low endemicity is often concentrated in small, isolated areas or “hot spots” comprised mostly of asymptomatic individuals still infectious to mosquitoes [2, 3]. Many malaria endemic countries have a surveillance system in place for identifying symptomatic cases in real-time (passive case detection or PCD); however, this system fails to reach asymptomatic individuals [1]. Active case detection (ACD) is a surveillance method recommended by the World Health Organization (WHO) in low transmission settings in which symptomatic and asymptomatic individuals are screened and treated for malaria [4]. Reactive case detection (RCD) is a form of ACD that was designed to take advantage of the spatial and temporal clustering of asymptomatic individuals within “hot spots” by using passively detected cases as triggers to initiate screening and treatment of individuals living in proximity to those cases [5, 6]. RCD is implemented in many countries working towards malaria elimination, including Zambia [6], South Africa [7], Brazil [8], Cambodia [9], and India [10].

The application of RCD in many of these settings differs in important features, such as the optimal screening radius and the number of households investigated [10]. In each instance, however, RCD is labour and time-intensive, requiring significant human resources, many rapid diagnostic tests (RDTs), and often long travel times between households [10, 11]. The utility of RCD in low transmission settings has been debated in part due to operational challenges during implementation and the use of less sensitive diagnostic tools such as RDTs and microscopy, which miss low density infections [1, 2, 10, 12]. Other limitations of RCD are its inability to reach populations with poor access, as well as the reliance on incident symptomatic cases seeking care to find “hot spots” comprised of asymptomatic individuals [10, 12].

Residual transmission in “hot spots” is driven by many local environmental factors such as vegetation and availability of aquatic habitats that determine vector density and heterogeneity [13]. For example, Anopheles larval sites contract and cluster around permanent aquatic habitats during the dry season, and expand during the wet season [14]. Various topographical features can also predict incident cases [13]. To increase efforts towards elimination, RCD may be improved by including environmental risk factors into the screening process, leveraging the heterogeneous nature of malaria transmission as a function of environmental features to guide asymptomatic case detection [13].

The Government of Zambia launched their RCD strategy in 2011 to enhance malaria surveillance and engage health systems at the community level to identify and treat individuals infected with Plasmodium falciparum who did not seek care or had minimal or no symptoms [15,16,17]. This RCD strategy is part of the National Malaria Elimination Strategic Plan (NMESP) to eliminate malaria in Zambia by 2021 and is employed in communities where parasite prevalence is approximately 1% and ten or fewer cases are passively detected at health facilities [15, 16, 18]. In Zambia, RCD starts with passive detection of a symptomatic malaria index case using P. falciparum histidine-rich protein-2 (PfHRP2) RDTs at a rural health clinic or by community health workers (CHWs) at rural health posts. CHWs then perform household visits to screen-and-treat residents within the index household as well as neighbouring or secondary households within a 140-m radius [15].

Studies have shown that environmental risk factors can be used to identify households likely to have parasitaemic residents [19,20,21,22,23,24]; however, the use of such environmental risk factors have not been explored in southern Zambia in the context of RCD. Building on prior work that assessed the efficiency of RCD in southern Zambia, the predictive ability of environmental risk factors was evaluated at varied spatial scales to identify parasitaemic residents of households located within a larger radius of 250 m from an index household [18, 25, 26].

Methods

Study site

Households were enrolled into the RCD study in the catchment area of Macha Hospital in Choma District, Southern Province, Zambia between January 12, 2015 and July 26, 2017 [18, 25, 27, 28]. The region has a tropical savannah climate with the rainy season occurring from December to April, followed by a cool dry season from May to August, and a hot dry season from September to November as previously described [15, 18, 22, 26]. Malaria transmission is propagated by the primary vector Anopheles arabiensis, which peaks during the rainy season. Infections are almost exclusively due to P. falciparum [18, 25, 29]. The major malaria control interventions are case management with artemisinin-based combination therapy (ACT) introduced in 2004, long-lasting insecticide-treated nets (LLINs) that were introduced in 2007 and redistributed approximately every 3 years with the most recent being in November 2017, and targeted mass drug administration (MDA) and indoor residual spraying (IRS) largely outside the study area [18, 26].

Reactive case detection

RCD eligibility and enrollment started at thirteen health centres and 23 health posts within the study catchment area where symptomatic individuals positive for malaria by PfHRP2-based RDT (index cases) triggered follow-up visits by a CHW and study team from Macha Research Trust [15, 18]. The study field team received notifications of an index case through SMS text messages from the health centre staff, after which they visited the household of the index case as well as secondary households located within a 250 m radius of the index case within 1 week of notification [15]. The RCD radius was expanded from the government recommended 140 to 250 m for the study. If the index case travelled outside their home district in the previous month and stayed overnight, they were not eligible for RCD screening through the government program. The field team was trained to administer consent, perform RDT testing, provide ACT for uncomplicated malaria, collect finger prick blood on filter paper, administer surveys, collect data using electronic tablets, and educate participants on malaria transmission and prevention [18].

Study population

The study population consisted of residents in an index case household and secondary households within 250 m of an index case. Households were single or multiple houses belonging to the main and extended family [30]. When index and secondary households were screened, all members of a household were eligible for enrollment. After written informed consent, including parental permission and assent for older children, a questionnaire was administered to obtain demographic information, knowledge of malaria transmission, malaria symptoms, travel history, and malaria prevention methods [18]. Parents or guardians completed surveys on behalf of participants younger than 16 years. A blood sample was collected for a PfHRP2-based RDT (SD Bioline, Standard Diagnostics Inc, Gyeonggi-do, Republic of Korea) and as dried blood spots (DBS) on Whatman 903™ Protein Saver cards (GE Healthcare Bio-Sciences Corporation, Piscataway, NJ) for quantitative PCR (qPCR) [18, 19]. Household residents found to be RDT positive were offered artemether/lumefantrine (Coartem®), the standard treatment for uncomplicated malaria in Zambia. Global positioning system (GPS) coordinates were obtained at each household using hand-held GPSMAP® 62 devices (Garmin Ltd, Olathe, Kansas) and mapped using ArcGIS version 10.5 (Environmental Systems Research Institute, Redlands, California) on a high resolution Quickbird™ satellite image of the catchment area [25, 28].

Parasite prevalence

Parasite prevalence was determined using the PfHRP2-based RDT results and detection of P. falciparum mitochondrial cytochrome b gene (cytb) by qPCR. PfHRP2 RDT readings were performed according to the manufacturer’s instructions [31]. DBS samples for qPCR were stored in plastic bags with desiccants and transported to the laboratory at Macha Research Trust for further drying. Samples were re-packed and stored at − 20 °C until parasite DNA extraction was performed using the Chelex® method [28, 32]. qPCR was performed with the Applied Biosystems StepOnePlus™ Real-Time PCR System (Thermo Scientific, Waltham, MA, USA). Primers specific to P. falciparum cytb were used to amplify, detect and quantify P. falciparum DNA using SYBR® Green fluorescence [30, 33, 34]. Filter paper spotted with laboratory-cultured parasites (NF54) and dilutions of 3D7 genomic DNA were used as standards [32]. The limit of detection was established as one parasite/uL [18]. The qPCR reaction consisted of 5 µL DNA template, 5 µL SYBR® Green PCR Master Mix (ThermoFisher), 200 nM forward primer (5′ CCT GAT AAT GCT ATC GTA 3′), and 200 nM reverse primer (5′ TAA TAC AAT TAC TAA ACC AGC 3′) [18]. Amplification with correct melting temperature was considered positive and the amplicon was further confirmed on a 4% agarose gel [18].

Environmental risk factors

A Quickbird™ satellite image of the 1200 km2 catchment area provided by the GeoEye-1 satellite (DigitalGlobe Services, Inc., Denver, Colorado) in 2017 and comprised of four-band 1.64-m spatial resolution and 0.41-m resolution panchromatic single-band imagery. Remote sensing data was imported into ArcGIS version 10.5 to geocode index and secondary households [25, 28]. All data layers were projected onto the Universal Transverse Mercator (UTM), Southern Hemisphere, Zone 35, WGS1984. A digital elevation model (DEM) with 90-m resolution was obtained from Shuttle Radar Topography Mission (SRTM) version 3, processed in ERDAS Imagine 2011 software (Hexagon Geospatial, Madison, Al) and imported into ArcGIS [18, 26]. The ArcHydro Tools module of ArcGIS was used to build a stream network according to the Strahler stream classification that assigned order values of 1, 2, 3, etc. based on a hierarchy of tributaries, such that two small first-order streams join to form a second-order stream [19, 24, 35]. A shapefile for roads was created by digitizing roads in ArcGIS based on a 1:50,000 topographic map of Zambia and the satellite image. Households with one or more malaria positive individuals by RDT or qPCR (excluding the index case) were classified as positive households [19, 25]. Index and secondary households were compared based on the following baseline characteristics: median age of residents per household, number of individuals per household, number of individuals 5 years and younger per household, number of parasitaemic individuals, insecticide-treated bed net ownership, floor material, and cooking energy source [21].

Secondary households with and without parasitaemic individuals were also compared using the same variables. In addition to household descriptive variables, environmental risk factors characterizing the surroundings of secondary households and previously shown to be associated with malaria risk were evaluated on the following levels: (1) household-level, defined as environmental risk factors within 100-m radius of a household; (2) cluster-level, defined as environmental risk factors within 250-m radius of a household; and (3) neighbourhood-level, defined as environmental risk factors outside the 250-m screening radius.

Household-level risk factors included the number of animal pens within a 100-m radius of the main house structure and distance to nearest animal pens [36, 37]. Cluster-level risk factors included distance to index households and elevation difference with index households [20, 23, 24]. Neighbourhood-level risk factors included distance to the main road and distance to streams [19, 23]. If the distance between the index and secondary households was more than 300 m, coordinates were manually cross-referenced and re-mapped by the field team when necessary. Elevation differences were generated by taking the difference in elevation as recorded by GPS devices from each secondary household and its corresponding index household. Missing elevation coordinates were extracted from the DEM. Animal pens were manually digitized in ArcGIS and were defined as enclosed, dark- or light-brown, oblong, circular, or rectangular roofless structures of any size within a 100-m radius of the main house structure. Animal pens that were visually problematic to identify in ArcGIS were cross-referenced with Google Earth Images captured in 2017. From the individual stream order distances, the closest stream to the secondary household was identified. All risk factors generated using ArcGIS were imported into STATA 14.2 for statistical analyses. Figure 1 illustrates secondary households within the 140- and 250-m screening radius of the positive index household and the proximity of animal pens and streams.

Fig. 1
figure 1

High-resolution Quickbird™ satellite image for catchment area of Macha Hospital in Choma District

Statistical analysis

The Chi square test for proportions and Wilcoxon rank-sum for means were used to compare household descriptive variables between index and secondary households, as well as secondary households with and without parasitaemic individuals. The analysis was restricted to all participants in secondary households who provided consent, completed the survey, and had RDT or PCR results. Generalized linear models with inference based on generalized estimating equations (GEE) were used to estimate the cross-sectional population average effect for the difference in odds of a positive vs. negative secondary households for each environmental risk factor. The GEE model was chosen for its ability to account for the clustering of secondary households around the index household and to more accurately estimate standard errors. The outcome was a binary variable, distinguishing secondary households with parasitaemic individuals and those without. The household-, cluster-, and neighbourhood-level environmental risk factors were assessed for collinearity using variance inflation factor values. Variables included in the model were: distance to index household (per 50 m), distance to main road (per 50 m), elevation difference with index household (per 10 m), number of animal pens within 100 m, presence of animal pen (yes vs. no), season, and a categorical variable identifying nearest streams order 1 through 5. Model fit was evaluated using the Hosmer–Lemeshow goodness of fit test and significance was evaluated at a p value of 0.05 (Additional file 1: Table S1).

Results

Characteristics of the study population

Between January 12, 2015 and July 26, 2017, 4170 individuals (excluding index cases) in 689 households received a screening visit, completed a survey, and had RDT or qPCR result available. Of the 689 households, 77% (531) were classified as secondary households comprised of 2926 residents. The median number of individuals per households was 5 (interquartile range [IQR]: 3–7) and 7 (IQR: 5–10) (P < 0.001) for secondary and index households, respectively, with index households having more children 5 years and younger (median: 2.0, IQR 1.0–2.0) (P < 0.001). Almost half (46.8%) of index households had at least one parasitaemic individual other than the index case by qPCR or RDT, while only 8.5% of secondary households had at least one parasitaemic individual (P < 0.001). There were no differences in median age per household, insecticide-treated bed net ownership, household floor material, and preferred cooking energy between index and secondary households (Table 1).

Table 1 Household characteristics comparing index and secondary households enrolled from January 2015–July 2017

The median age within negative secondary households was 17.8 years (IQR: 13.7–25.4 years) and that within positive secondary households was 15.9 years (IQR: 13.1–20.7 years) (P = 0.16) (Table 2). The median number of individuals per households was 5 (IQR: 3–7) and 7 (IQR: 5–8) (P < 0.001) for negative and positive secondary households, respectively, with more individuals 5 years and younger residing in positive secondary households (P < 0.01) (Table 2). The composition of parasitaemic individuals in secondary households ranged from one (78% of households) to three individuals (2% of households) (P < 0.001). No differences were found in insecticide-treated bed net ownership, household floor material, and preferred cooking energy between positive and negative secondary households.

Table 2 Household characteristics comparing negative and positive secondary households enrolled January 2015–July 2017

Malaria prevalence

Excluding index cases, 153 participants were positive for malaria by RDT or qPCR, 37% (56) of whom resided in 45 secondary households. The parasite prevalence for residents in secondary households was 0.7% using RDT and 1.8% using qPCR, while that for residents in index households was 2.7% and 7.3%, respectively. In secondary households, the parasite prevalence remained below 2.0% during the rainy, cool dry, and hot dry season, while a pattern of seasonal transmission was observed among index households with RDT parasite prevalence of 2.8%, 3.3%, and 1.2% during the rainy, cool dry, and hot dry season, respectively, and qPCR parasite prevalence of 8.5%, 6.1%, and 4.1%, for the corresponding seasons for index households (Fig. 2).

Fig. 2
figure 2

Parasite prevalence (%) for index and secondary households, January 2015–July 2017 in southern Zambia

Environmental household risk factors

Secondary households with parasitaemic individuals were not significantly further away from index households (median: 179.8 m, IQR: 121.7–226.9 m) compared to those without parasitaemic individuals (median: 164.7 m, IQR: 104.1–210.1 m) (P = 0.17) but they were in closer proximity to streams (median: 335.3 m, IQR: 242.2–539.1 m) (P < 0.01) (Table 3). Although no statistically significant differences were detected for the other environmental risk factors, positive secondary households exhibited a trend previously associated with increased malaria risk [19,20,21, 26, 36,37,38]: they were at lower elevation from the index household (median: − 0.7, IQR: − 8.9 to 9.4 m) (P = 0.52); further away from the main road (median: 9266.6 m, IQR: 4130.4–17,353.7 m) (P = 0.28); and more likely to have an animal pen (64.4%) (P = 0.69). If an animal pen was present, they were also closer to it than negative secondary households (median: 36.9 m, IQR: 20.6–50.5 m) (P = 0.43) (Table 3).

Table 3 Environmental household risk factors comparing negative and positive secondary households enrolled January 2015–July 2017

Efficiency of reactive case detection

The association between positive secondary households and environmental risk factors at the household-, cluster-, and neighbourhood-level, was evaluated for 45 positive and 485 negative secondary households. One negative secondary household was excluded as it was > 300 m from the index household. Cluster- and neighbourhood-level risk factors were associated with positive secondary households in the multivariate GEE model, while no significant association was observed for household-level risk factors (Table 4). The cluster-level risk factor that was marginally associated with a positive secondary household was increasing distance to the index household (P = 0.07). As the distance to index households increased by 50 m, the odds of a positive secondary household increased by 24% (adjusted odds ratio [OR]: 1.24, 95% confidence interval [CI] 0.98–1.58). No association was observed for elevation. The neighbourhood-level risk factor associated with positive secondary households were third- and fifth-order streams. Positive secondary households were 3.14 times more likely to be located near a third-order stream (OR: 2.97, 95% CI 1.04–8.42) (P = 0.041) and 3.20 times more likely to be located near a fifth-order stream (OR: 2.30, 95% CI 1.04–5.09) (P = 0.040), compared to a negative secondary household. Fourth-order streams (OR: 1.62, 95% CI 0.21–12.65) also exhibited a trend of increased risk; however, the association was not statistically significant (P = 0.64). No association was observed for the distance to the main road. The household-level risk factor of number of animal pens within 100-m was not associated with positive secondary households; however, if an animal pen was present, the odds of being a positive household increased by 60% (OR: 1.60, 95% CI 0.57–4.47), although the association was not statistically significant (P = 0.37).

Table 4 Crude and adjusted OR for the association between environmental risk factors and positive secondary households

An increased risk of identifying a positive secondary household near streams suggests that environmental features could potentially guide RCD screening strategies. To evaluate if this increased probability of identifying secondary household with parasitaemic individuals would require fewer households to be screened, streams closest to an index household were used to find secondary households within the 250-m radius. Using only streams, a total of 137 secondary households, 14% (n = 19) that had parasitaemic residents, representing 42% of all positive secondary households in the study sample, would have been identified. The current RCD screening method using only the 140-m radius identified 210 secondary households; however, only 8% (n = 16) of these secondary households had parasitaemic residents, representing 36% of all positive secondary households in the study sample. These results indicate that incorporating environmental risk factors such as streams, within a larger screening radius, would allow for more parasitaemic individuals to be identified while screening fewer number of households, consequently increasing the efficiency of the RCD programme [25]. If streams nearest to index households were used to find secondary households within a 250-m radius, 16% (n = 3) more secondary households with parasitaemic residents would have been screened and treated compared to the national RCD strategy with 140-m screening radius.

Discussion

Environmental risk factors were associated with the probability of finding households with parasitaemic residents using RCD as demonstrated in other studies in Zambia [13]. In the low transmission setting of Choma District, Zambia, identifying streams located near index households to guide and direct screening has the potential to improve RCD and affect transmission by identifying households with asymptomatic infections. These findings are in line with a previous study conducted within the same study area in 2008 where it was shown that households within 1.98 km from a third-order stream were 2.8 (95% CI 1.2–6.9) times more likely to have an RDT positive resident than those within 6 km [26].

Although no associations were found with the other environmental risk factors such as distance to a main road, elevation, season, and number and presence of animal pens, non-parametric comparisons between positive and negative secondary households exhibited a trend of increased malaria risk for these risk factors [19,20,21, 26, 36,37,38]. The risk associated with animal pens varies in the literature depending on vector behaviour. Anopheles arabiensis has been reported to be anthropophilic in southern Zambia but also displays zoophilic habits by feeding opportunistically on non-human blood sources [36]. Other vectors besides An. arabiensis might also have an important role in transmission as P. falciparum-infected Anopheles squamosus exhibiting outdoor zoophagic feeding behaviour were recently identified in the area. Early studies in Choma District, Zambia found that ownership of cattle reduced the risk of P. falciparum infection by 87% while others have found less conclusive evidence [36, 37]. For elevation, however, it has been clearly shown that increased elevation offers protection against malaria infection [13, 20, 24, 26, 38]. Since index and secondary households in this study were located only < 300 m from each other and variation in elevation was minimal, it is unlikely that the elevation would influence malaria risk at this spatial scale. Distance from the index household marginally increased the probability of finding positive secondary households (OR: 1.24, 95% CI 0.98–1.58), in contrast to other studies. Larsen et al. observed a decreased risk for neighbouring households located further away, and Bulterys et al. found an adjusted OR of 0.26 (95% CI 0.07–0.98) as distance between households increased. Finally, distance to the main road has often been treated as an indicator of increased malaria risk. In Chongwe District, Zambia, the odds of RDT positive households increased by 5% for every 500-m increase in distance from the road [39]. As only proximity to the main road was examined, it is possible that constant use from vehicles, animal carts, and people prevented mosquito breeding sites from developing undisturbed, reducing this as a risk factor. Less frequently used subsidiary roads and rural paths (not included in the analysis) could provide more opportune mosquito breeding sites closer to residences as their composition allows for easier accumulation of aquatic breeding sites compared to the tarmac and concrete main road.

The use of environmental risk factors for malaria risk prediction is a common approach to malaria control and has been employed in various transmission settings around the world. For example, in Chimoio, Mozambique a GIS-based spatial model was designed to estimate areas of risk using temperature, precipitation, altitude, slope, distance to water bodies, distance to roads, normalized difference vegetation index (NDVI), land use, and land cover [40]. The model identified that 4% of Chimoio was at high risk for contracting malaria, with precipitation as a key risk factor for the entire area studied [40]. Another study in south Sumatra, Indonesia used ordinary least square and geographically weighted regression to show that altitude, distance to forest, and rainfall determined overall malaria incidence with considerable heterogeneity at the village level [41]. These findings were consistent with other studies in Cambodia, Addis Ababa, Ethiopia, and Rondôia, Brazil [41]. Despite the extensive literature on environmental risk factors for malaria, their application within the context of RCD has been limited.

Many studies evaluating the efficiency of RCD highlight its inability to halt infections in areas of low transmission due to the use of less sensitive RDTs, travel-related infections, and large budgetary requirements [2, 18]. A major concern for RCD-based strategies is that asymptomatic individuals will be missed if no clinical cases report to CHWs [42]. A survey in coastal Kenya found that asymptomatic and symptomatic infections do not necessarily overlap spatially, and that clusters of symptomatic infections have greater temporal stability over more than 10 years [42]. Another issue often highlighted is the different criteria and screening radii employed by countries to define and recruit neighbouring households [42] For example, RCD data from four villages in the Myanmar-Thailand border determined that RCD would only be successful at a radius of 150 m, and any screening occurring beyond this radius would not perform better than random screening [2]. Another study in Pailin Province in western Cambodia screened the nearest five households for every fifteenth index case and ten nearest households for every 30th index case. Using this approach, they predicted only 40% of infections and concluded that RCD was not recommended in a setting targeting elimination [43]. However, with the shortcomings of a circular radius and the various implementation challenges, for RCD to be an effective method for malaria elimination in these low-endemic countries a tailored approach adapted to the local parasite epidemiology, vector biology, and living/working environment of the community must be considered key for it to succeed.

This study used environmental risk factors for malaria to characterize the low transmission setting to improve RCD efficiency. Previous work on enhancing RCD efficiency in Southern Province, Zambia has also shown that time-invariant measures of the environment, such as the topographic position index (TPI; measure of an area’s relative elevation to find slopes, valleys, and ridges), the convergence index (CI; measure of an area’s propensity to pool water), median enhanced vegetation index (EVI; measure of vegetation density), and the topographical wetness index (TWI; measure of water flow) were stronger predictors for identifying parasitaemic individuals than demographics of incident symptomatic cases [13, 26]. They showed that ridges and upper slopes (at a TPI scale of 270 m) and wetter regions (TWI > 10.2) were associated with finding more parasitaemic individuals during RCD [13]. These findings, along with the current study, support the significance of water bodies in improving the efficiency of RCD. Third through fifth-order streams are mid-level streams that may not always be suitable for larval development; however, larvae have been collected from water at the edges of these streams (unpublished findings). During the dry season, as water accumulates into smaller pools, they become ideal larval development sites. These streams can also serve as important markers for nearby areas with similar high water table harbouring larvae [13, 26, 38, 44]. And as these streams can be challenging to locate depending on size and season, spatial risk maps with topographical measures, such as CI and TWI, can offer guidance to CHWs to possibly reach clusters of asymptomatic carriers otherwise missed during regular RCD screening. Other water sources such as dams, are also important determinants for malaria transmission as was shown in Ethiopia, where reservoir water level management suppressed larval development [45].

In addition to the use of streams as a screening tool, RCD efficiency could benefit from the combined use of RDTs and highly sensitive qPCR. For this study region, the overall parasite prevalence (3.7%) was mostly driven by qPCR as parasite prevalence by RDT was only 1.3%. Although costly, sensitive molecular methods such as qPCR are critical in low endemic settings to detect potential parasite-transmitting asymptomatic carriers. Even ultra-sensitive RDTs (uRDTs), such as the new Alere™ Malaria Ag P.f uRDT which was designed for low transmission settings, may not be sufficiently sensitive alternatives to SD Bioline PfHRP2 RDTs [46]. The Alere™ us-RDT has a ten-fold lower limit of detection for PfHRP2 compared to regular RDTs but missed 56% of PCR-detectable P. falciparum infections in a low-endemic setting in Myanmar, and in Papua New Guinea the test missed 50% of P. falciparum infections otherwise detectable by qPCR [47].

There were several limitations to this study. Restricting environmental variables within set radii raises concerns for edge-effect associations. For example, animal pens located just outside the 100-m radius were not counted as belonging to neighbouring households, thus potentially underestimating the number of animal pens associated with a household. Not all environmental risk factors important for malaria transmission were evaluated. Vegetation cover, an important indicator of available mosquito habitat, could also be a useful screening tool and has yet to be evaluated for RCD strategies [13, 48, 49]. Finally, the risk factors shown to be associated with positive households in this low transmission setting of Choma District, Zambia may not be applicable in other endemic regions.

The effectiveness of RCD ultimately depends on the number of cases found and treated in a timely manner and the resources allocated during implementation [2]. However, it is important to consider the added value of a tailored RCD approach based on demographic and ecological risk factors and more sensitive diagnostic tools to fully reap the benefits of this screening method to achieve malaria elimination [1]. In Cambodia, where infection is linked to occupation and mobility, an expanded RCD approach was implemented in which individuals who were coworkers of a symptomatic index case in settings of high malaria infection, such as forests and plantations, were also screened [1]. The expanded RCD had a detection rate of 3.9% compared to 0–2% using the classic RCD approach [1]. Through this adapted RCD design, Cambodia’s National Malaria Control Programme sought to identify and treat asymptomatic individuals within a discrete population whose members shared a common malaria risk through occupations such as logging, mining, and migrant labour [1]. The environmental risk factors identified in this study demonstrate that, even in low transmission settings, a tailored approach is possible; however, further work is needed to fully understand how these risk factors vary across district and season and how they can be modified to guide RCD strategies nationally in Zambia.

Conclusion

This study identified higher order streams as risk factors for parasitaemia in households neighbouring an index case as part of RCD in rural southern Zambia. These risk factors have the potential to improve the efficiency of RCD in a low transmission setting by not only identifying parasitaemic individuals more efficiently but also potentially reducing the number of households needed to be screened. Combined with other strategies for malaria elimination, tailored RCD approaches can help realize the goal of malaria elimination in Zambia.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request and through ClinEpiDB (https://clinepidb.org/ce/app/).

Abbreviations

ACD:

Active case detection

ACT:

Artemisinin-based combination therapy

aOR:

Adjusted odds ratio

CHW(s):

Community health worker(s)

CI:

Convergence index

Cytb :

Cytochrome b gene

DBS:

Dried blood spots

DEM:

Digital elevation model

EVI:

Enhanced vegetation index

GEE:

Generalized estimating equations

GIS:

Geographic information system

GPS:

Global positioning system

IRS:

Indoor residual spraying

IQR:

Inter quartile range

LLINs:

Long lasting insecticide treated nets

MDA:

Mass drug administration

NDVI:

Normalized difference vegetation index

NMESP:

National Malaria Elimination Strategic Plan

OR:

Odds ratio

PCD:

Passive case detection

PfHRP2:

P. falciparum histidine-rich protein-2

qPCR:

Quantitative polymerase chain reaction

RCD:

Reactive case detection

RDT(s):

Rapid diagnostic test(s)

SRTM:

Shuttle Radar Topography Mission

TPI:

Topographic position index

TWI:

Topographical wetness index

uRDT(s):

Ultra-sensitive RDT(s)

UTM:

Universal transverse mercator

WHO:

World Health Organization

References

  1. Rossi G, Van den Bergh R, Nguon C, Debackere M, Vernaeve L, Khim N, et al. Adapting reactive case detection strategies for falciparum malaria in a low-transmission area in Cambodia. Clin Infect Dis. 2017;66:296–8.

    Article  Google Scholar 

  2. Parker DM, Landier J, Von Seidlein L, Dondorp A, White L, Hanboonkunupakarn B, et al. Limitations of malaria reactive case detection in an area of low and unstable transmission on the Myanmar-Thailand border. Malar J. 2016;15:571.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Cotter C, Sturrock HJW, Hsiang MS, Liu J, Phillips AA, Hwang J, et al. The changing epidemiology of malaria elimination: new strategies for new challenges. Lancet. 2013;382:900–11.

    Article  PubMed  Google Scholar 

  4. WHO. Disease surveillance for malaria elimination: An operational manual. Geneva, World Health Organization, 2012. https://apps.who.int/iris/bitstream/handle/10665/44852/9789241503334_eng.pdf;jsessionid=051421372DD46E2264D749EA8463FD5B?sequence=1. Accessed 01 Mar 2019.

  5. Moonen B, Cohen JM, Snow RW, Slutsker L, Drakeley C, Smith DL, et al. Operational strategies to achieve and maintain malaria elimination. Lancet. 2010;376:1592–603.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Stresman GH, Kamanga A, Moono P, Hamapumbu H, Mharakurwa S, Kobayashi T, et al. A method of active case detection to target reservoirs of asymptomatic malaria and gametocyte carriers in a rural area in Southern Province, Zambia. Malar J. 2010;9:265.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Gerardin J, Bever CA, Hamainza B, Miller JM, Eckhoff PA, Wenger EA. Optimal population-level infection detection strategies for malaria control and elimination in a spatial model of malaria transmission. PLoS Comput Biol. 2016;12:e1004707.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Branch O, Casapia WM, Gamboa DV, Hernandez JN, Alava FF, Roncal N, et al. Clustered local transmission and asymptomatic Plasmodium falciparum and Plasmodium vivax malaria infections in a recently emerged, hypoendemic Peruvian Amazon community. Malar J. 2005;4:27.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Parker DM, Matthews SA, Yan G, Zhou G, Lee MC, Sirichaisinthop J, et al. Microgeography and molecular epidemiology of malaria at the Thailand-Myanmar border in the malaria pre-elimination phase. Malar J. 2015;14:198.

    Article  PubMed  PubMed Central  Google Scholar 

  10. van Eijk AM, Ramanathapuram L, Sutton PL, Kanagaraj D, Priya GSL, Ravishankaran S, et al. What is the value of reactive case detection in malaria control? A case-study in India and a systematic review. Malar J. 2016;15:67.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. Zelman BW, Baral R, Zarlinda I, Coutrier FN, Sanders KC, Cotter C, et al. Costs and cost-effectiveness of malaria reactive case detection using loop-mediated isothermal amplification compared to microscopy in the low transmission setting of Aceh Province, Indonesia. Malar J. 2018;17:220.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Sturrock HJW, Novotny JM, Kunene S, Dlamini S, Zulu Z, Cohen JM, et al. Reactive case detection for malaria elimination: real-life experience from an ongoing program in Swaziland. PLoS ONE. 2013;8:e63830.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Larsen DA, Bridges DJ, Miller J, Cheelo S, Winters A, Ngwenya-Kangombe T, et al. Location, location, location: environmental factors better predict malaria-positive individuals during reactive case detection than index case demographics in Southern Province, Zambia. Malar J. 2017;16:18.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Bousema T, Griffin JT, Sauerwein RW, Smith DL, Churcher TS, Takken W, et al. Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Med. 2012;9:e1001165.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Larsen DA, Chisha Z, Winters B, Mwanza M, Kamuliwo M, Mbwili C, et al. Malaria surveillance in low-transmission areas of Zambia using reactive case detection. Malar J. 2015;14:465.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. USAID. President’s Malaria Initiative. Zambia. Malaria Operational Plan FY; 2019. https://www.pmi.gov/docs/default-source/default-document-library/malaria-operational-plans/fy19/fy-2019-zambia-malaria-operational-plan.pdf?sfvrsn=3. Accessed 19 Apr 2020.

  17. National Malaria Control Programme. National Malaria Elimination Strategic Plan 2017–2021. 2015. https://www.nmec.org.zm/malaria-elimination-strategic-plan. Accessed 30 Mar 2018.

  18. Deutsch-Feldman M, Hamapumbu H, Lubinda J, Musonda M, Katowa B, et al. Efficiency of a malaria reactive test-and-treat program in Southern Zambia: a prospective, observational study. Am J Trop Med Hyg. 2018;98:1382–8.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Pinchoff J, Chaponda M, Shields T, Lupiya J, Kobayashi T, Mulenga M, et al. Predictive malaria risk and uncertainty mapping in Nchelenge District, Zambia: evidence of widespread, persistent risk and implications for targeted interventions. Am J Trop Med Hyg. 2015;93:1260–7.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Chirombo J, Lowe R, Kazembe L. Using structured additive regression models to estimate risk factors of malaria: analysis of 2010 Malawi malaria indicator survey data. PLoS ONE. 2014;9:e101116.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Graves PM, Richards FO, Ngondi J, Emerson PM, Shargie EB, Endeshaw T, et al. Individual, household and environmental risk factors for malaria infection in Amhara, Oromia and SNNP regions of Ethiopia. Trans R Soc Trop Med Hyg. 2009;103:1211–20.

    Article  PubMed  Google Scholar 

  22. 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;93:5–15.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Clennon JA, Kamanga A, Musapa M, Shiff C, Glass GE. Identifying malaria vector breeding habitats with remote sensing data and terrain-based landscape indices in Zambia. Int J Health Geogr. 2010;9:58.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Pinchoff J, Chaponda M, Shields TM, Sichivula J, Muleba M, Mulenga M, et al. Individual and household level risk factors associated with malaria in Nchelenge District, a region with perennial transmission: a serial cross-sectional study from 2012 to 2015. PLoS ONE. 2016;11:e0156717.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Searle KM, Shields T, Hamapumba H, Kobayashi T, Mharakurwa S, Thuma PE, et al. Efficiency of household reactive case detection for malaria in rural Southern Zambia: simulations based on cross-sectional surveys from two epidemiological settings. PLoS ONE. 2013;8:e70972.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Moss WJ, Hamapumba H, Kobayashi T, Shields T, Kamanga A, Clennon J, et al. Use of remote sensing to identify spatial risk factors for malaria in a region of declining transmission: a cross-sectional and longitudinal community survey. Malar J. 2011;10:163.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Mharakurwa S, Thuma PE, Norris DE, Mulenga M, Chalwe V, Chipeta J, et al. Malaria epidemiology and control in Southern Africa. Acta Trop. 2012;121:202–6.

    Article  PubMed  Google Scholar 

  28. Searle KM, Katowa B, Kobayashi T, Siame MNS, Mharakurwa S, Carpi G, et al. Distinct parasite populations infect individuals identified through passive and active case detection in a region of declining malaria transmission in Southern Zambia. Malar J. 2017;16:154.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Kent RJ, Thuma PE, Mharakurwa S, Douglas EN. Seasonality, blood feeding behavior, and transmission of Plasmodium falciparum by Anopheles arabiensis after an extended drought in Southern Zambia. Am J Trop Med Hyg. 2007;76:267–74.

    Article  PubMed  Google Scholar 

  30. WHO. World Malaria Report. Geneva, World Health Organization; 2016. http://apps.who.int/iris/bitstream/handle/10665/252038/9789241511711-eng.pdf;jsessionid=8AB072B00140FD1DFB1331B2C69C43C1?sequence=1. Accessed 05 Feb 2018.

  31. Boyce MR, O’Meara WP. Use of malaria RDTs in various health contexts across Sub-Saharan Africa: a systematic review. BMC Public Health. 2017;17:470.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Laban NM, Kobayashi T, Hamapumba H, Sullivan D, Mharakurwa S, Thuma PE, et al. Comparison of a PfHRP2-based rapid diagnostic test and PCR for malaria in a low prevalence setting in rural Southern Zambia: implications for elimination. Malar J. 2015;14:25.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Xu W, Morris U, Aydin-Schmidt B, Msellem M, Shakely D, Petzold M, et al. SYBR green real-time PCR-RFLP assay targeting the Plasmodium cytochrome B gene—a highly sensitive molecular tool for malaria parasite detection and species determination. PLoS ONE. 2015;10:e0120210.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Ranadive N, Kunene S, Darteh S, Ntshalintshali N, Nhlabathi N, Dlamini N, et al. Limitations of rapid diagnostic testing in patients with suspected malaria: a diagnostic accuracy evaluation from Swaziland, a low-endemicity country aiming for malaria elimination. Clin Infect Dis. 2017;64:1221–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Strahler A. Dynamic basis of geomorphology. Geol Soc Am Bull. 1952;63:923–38.

    Article  Google Scholar 

  36. Bulterys PL, Mharakurwa S, Thuma PE. Cattle, other domestic animal ownership, and distance between dwelling structures are associated with reduced risk of recurrent Plasmodium falciparum infection in Southern Zambia. Trop Med Int Health. 2009;14:522–8.

    Article  PubMed  Google Scholar 

  37. Mayagaya VS, Nkwengulila G, Lyimo IN, Kihonda J, Mtambala H, Ngonyani H, et al. The Impact of livestock on the abundance, resting behaviour and sporozoite rate of malaria vectors in Southern Tanzania. Malar J. 2015;14:17.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Gómez-Barroso D, García-Carrasco E, Herrador Z, Ncogo P, Romay-Barja M, Mangue MEO, et al. Spatial clustering and risk factors of malaria infections in Bata District, Equatorial Guinea. Malar J. 2017;16:146.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Pinchoff J, Henostroza G, Carter BS, Roberts ST, Hatwiinda S, Hamainza B, et al. Spatial patterns of incident malaria cases and their household contacts in a single clinic catchment area of Chongwe District, Zambia. Malar J. 2015;4:305.

    Article  Google Scholar 

  40. Ferrao JL, Niquisse S, Mendes JM, Painho M. Mapping and modelling malaria risk areas using climate, socio-demographic and clinical variables in Chimoio, Mozambique. Int J Environ Res Public Health. 2018;15:E79.

    Article  Google Scholar 

  41. Hasyim H, Nursafingi A, Haque U, Montag D, Groneberg DA, Dhimal M, et al. Spatial modelling of malaria cases associated with environmental factors in South Sumatra, Indonesia. Malar J. 2018;17:87.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Fontoura PS, Finco BF, Lima NF, de Carvalho JF, Vinetz JM, Castro MC, et al. Reactive case detection for Plasmodium vivax malaria elimination in rural Amazonia. PLoS Negl Trop Dis. 2016;10:e0005221.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Hustedt J, Canavati SE, Rang C, Ashton RA, Khim N, Berne L, et al. Reactive case-detection of malaria in Pailin Province, Western Cambodia: lessons from a year-long evaluation in a pre-elimination setting. Malar J. 2016;15:132.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Texas Aquatic Science. Streams and Rivers; 2019. https://texasaquaticscience.org/streams-rivers-aquatic-science-texas/. Accessed 24 Mar 2019.

  45. Kibret S, Wilson GG, Ryder D, Tekie H, Petros B. Environmental and meteorological factors linked to malaria transmission around large dams at three ecological settings in Ethiopia. Malar J. 2019;18:54.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Das S, Jang IK, Barney B, Peck R, Rek JC, Arinaitwe E, et al. Performance of a high-sensitivity rapid diagnostic test for Plasmodium falciparum malaria in asymptomatic individuals from Uganda and Myanmar and naive human challenge infections. Am J Trop Med Hyg. 2017;97:1540–50.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Hofmann NE, Greunberg M, Nate E, Ura A, Rodriguez-Rodriguez D, Salib M, et al. Assessment of ultra-sensitive malaria diagnosis versus standard molecular diagnostics for malaria elimination: an in-depth molecular community cross-sectional study. Lancet Infect Dis. 2018;18:1108–16.

    Article  PubMed  Google Scholar 

  48. Rubio-Palis Y, Bevilacqua M, Medina DA, Moreno JE, Cardenas L, Sanchez V, et al. Malaria entomological risk factors in relation to land cover in the lower Caura river basin, Venezuela. Mem Inst Oswaldo Cruz. 2013;108:220–8.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Ricotta EE, Frese SA, Choobwe C, Louis TA, Shiff CJ. Evaluating local vegetation cover as a risk factor for malaria transmission: a new analytical approach using ImageJ. Malar J. 2014;13:94.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We acknowledge the Zambian National Malaria Elimination Programme and local health centre staff for their support. We thank the members of the community in Macha, Zambia for their participation in the surveys and the field and laboratory teams at the Macha Research Trust for their work and commitment towards the Southern and Central Africa International Centers of Excellence for Malaria Research.

Funding

This analysis was supported by the Johns Hopkins Malaria Research Institute, the Bloomberg Family Foundation, and the Division of Microbiology and Infectious Diseases, National Institute of Allergy and Infectious Diseases as part of the International Centers of Excellence for Malaria Research (U19AI089680).

Author information

Authors and Affiliations

Authors

Consortia

Contributions

Original idea and study design contributed by WJM. HH, oversaw all field operations and managed data collection along with ML, and JM under guidance from JCS and PET. MM and BK performed DNA extractions and qPCR under guidance of TK. Cleaning and creation of analytical dataset in STATA and environmental predictors in ArcGIS was performed by FRPB, advised by TMS, KMS, and WJM. Figures and tables for manuscript were created by FRPB, in addition to analysing and interpreting the data. TMS, KMS, and ML created the spatial risk map. KMS helped with choosing statistical models. FRPB wrote manuscript and WJM was a major contributor in writing the manuscript. All authors read and approved the final manuscript

Corresponding author

Correspondence to Fiona R. P. Bhondoekhan.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Johns Hopkins Bloomberg School of Public Health Institutional Review Board and the Tropical Diseases Research Center Ethics Review Committee.

Consent for publication

The publication of this manuscript was approved by the National Health Research Authority—Zambia prior to publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary information

Additional file 1: Table S1.

QIC goodness-of-fit statistics for the GEE models.

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 http://creativecommons.org/licenses/by/4.0/. 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 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

Bhondoekhan, F.R.P., Searle, K.M., Hamapumbu, H. et al. Improving the efficiency of reactive case detection for malaria elimination in southern Zambia: a cross-sectional study. Malar J 19, 175 (2020). https://doi.org/10.1186/s12936-020-03245-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-020-03245-1

Keywords