Distinct parasite populations infect individuals identified through passive and active case detection in a region of declining malaria transmission in southern Zambia
Malaria Journal volume 16, Article number: 154 (2017)
Substantial reductions in the burden of malaria have been documented in parts of sub-Saharan Africa, with elimination strategies and goals being formulated in some regions. Within this context, understanding the epidemiology of low-level malaria transmission is crucial to achieving and sustaining elimination. A 24 single-nucleotide-polymorphism Plasmodium falciparum molecular barcode was used to characterize parasite populations from infected individuals identified through passive and active case detection in an area approaching malaria elimination in southern Zambia.
The study was conducted in the catchment area of Macha Hospital in Choma District, Southern Province, Zambia, where the parasite prevalence declined over the past decade, from 9.2% in 2008 to less than 1% in 2013. Parasite haplotypes from actively detected, P. falciparum-infected participants enrolled in a serial cross-sectional, community-based cohort study from 2008 to 2013 and from passively detected, P. falciparum-infected individuals enrolled at five rural health centres from 2012 to 2015 were compared. Changes in P. falciparum genetic relatedness, diversity and complexity were analysed as malaria transmission declined.
Actively detected cases identified in the community were most commonly rapid diagnostic test negative, asymptomatic and had submicroscopic parasitaemia. Phylogenetic reconstruction using concatenated 24 SNP barcode revealed a separation of parasite haplotypes from passively and actively detected infections, consistent with two genetically distinct parasite populations. For passively detected infections identified at health centres, the proportion of detectable polyclonal infections was consistently low in all seasons, in contrast with actively detected infections in which the proportion of polyclonal infections was high. The mean genetic divergence for passively detected infections was 34.5% for the 2012–2013 transmission season, 37.8% for the 2013–2014 season, and 30.8% for the 2014–2015 season. The mean genetic divergence for actively detected infections was 22.3% in the 2008 season and 29.0% in the 2008–2009 season and 9.9% across the 2012–2014 seasons.
Distinct parasite populations were identified among infected individuals identified through active and passive surveillance, suggesting that infected individuals detected through active surveillance may not have contributed substantially to ongoing transmission. As parasite prevalence and diversity within these individuals declined, resource-intensive efforts to identify the chronically infected reservoir may not be necessary to eliminate malaria in this setting.
Substantial reductions in the burden of malaria have been documented in parts of sub-Saharan Africa  and malaria elimination goals have been proposed at regional, national, and subnational levels [2,3,4]. As areas transition from high or moderate to low malaria transmission and approach elimination, understanding the epidemiology of low-level transmission is crucial to achieving and sustaining elimination [2, 3, 5]. The relative magnitude of infection attributed to locally acquired and imported cases is of interest, specifically the role of the chronically infected, asymptomatic reservoir in maintaining local malaria transmission [6,7,8,9,10]. One way to address this question is to determine the genetic relatedness of parasites between actively detected, predominantly chronically infected, asymptomatic individuals and passively detected, predominantly acutely infected, symptomatic malaria cases . By comparing the genetic overlap between parasites infecting these populations, inferences can be made regarding the relative contribution of the asymptomatic reservoir to symptomatic malaria cases. Understanding the role of these two populations in local malaria transmission dynamics can guide interventions, particularly reactive case detection and focal or mass drug administration strategies that target the chronically infected, asymptomatic reservoir.
A Plasmodium falciparum molecular barcode assay consisting of 24 unlinked, single nucleotide polymorphisms (SNPs) has been used to characterize unique genetic signatures and track circulating P. falciparum parasite populations . The barcode was developed to elucidate malaria transmission dynamics by tracking the genetic diversity and complexity of the parasite over time and space [12,13,14], and specifically for use in resource-limited settings as the highest level of technology required is a polymerase chain reaction (PCR) assay .
This study was conducted in an area of southern Zambia that experienced a dramatic decrease in malaria transmission over the past decade . The molecular barcode was used to determine the relatedness between parasites over time in infected individuals identified through passive and active surveillance from the same geographic area, as well as changes in parasite genetic complexity and diversity.
Study site and population
The study was conducted in the catchment area of Macha Hospital in Choma District, Southern Province, Zambia, where there is a single rainy season which lasts from November through April, followed by a cool dry season from April until August and a hot dry season from August through November. Malaria transmission peaks during the rainy season  and the primary vector is Anopheles arabiensis [16, 17]. The hospital catchment area is approximately 1200 km2 with roughly 30,000 residents consisting primarily of villagers living in small, scattered homesteads. The prevalence of P. falciparum infection declined in this area over the past decade, from 9.2% in 2008 to less than 1% in 2013 . Artemisinin-based combination therapy (ACT) with artemether–lumefantrine was introduced as first-line anti-malarial therapy in Zambia in 2002  and into the study area in 2004. In Zambia, long-lasting insecticide-treated nets (LLINs) are distributed through antenatal care clinics and additional mass distribution campaigns . LLINs were widely distributed in the study area in 2007  and more than 11,000 LLINs were distributed from nine health posts in the catchment area of Macha Hospital in June 2012, with additional LLINs distributed in 2014 according to the Office of the Macha Hospital Environmental Health Technician. ITN ownership was estimated at 83% in the study area in 2013 .
Active community-based malaria surveillance
Satellite images were used to develop a sampling frame for the random selection of households to enroll participants in longitudinal and cross-sectional malaria surveys . The identification and enumeration of households was done manually to delineate household and non-household structures (kraals, schools, and larger structures) . Households selected from the sampling frame were enrolled in either a cross-sectional or longitudinal cohort. Households enrolled in the cross-sectional cohort were surveyed once, whereas households enrolled in the longitudinal cohort were repeatedly surveyed every 2 months. Surveys were conducted each month, alternating between the cross-sectional cohort and the longitudinal cohort, from February 2008 through October 2013. For each study visit, a questionnaire was administered to collect demographic data, history of recent malaria symptoms and treatment, healthcare-seeking behaviour, knowledge of malaria risk and prevention, and long-lasting insecticidal net use. Recent malaria symptoms were defined as having a documented fever greater than 38 °C by tympanic temperature, or reporting a fever and chills within the previous 48 h. A blood sample was collected by finger prick for a rapid diagnostic test (RDT), microscopy and blood was spotted on Whatman 903™ Protein Saver cards. These dried blood spots (DBS) were used to detect Plasmodium parasites by a nested PCR as described below and molecular barcoding was attempted for all samples positive for Plasmodium parasites by nested PCR (Additional file 1: Figure S1). Participants positive by RDT were offered treatment with artemether-lumefantrine (Coartem®).
Passive malaria surveillance at rural health centres
Fourteen rural health centres (RHCs) that serve the catchment population surrounding Macha Hospital sent a weekly text message report of the number of RDTs used, number of positive RDTs, and the number of people treated for malaria to study staff at Macha Research Trust . Five of these RHCs within the study area where active surveillance was conducted were included in a sub-study from 2012 through 2015 that collected DBS from a convenience sample of individuals with a positive RDT and recorded their demographic data. Molecular barcoding was attempted for all samples collected at the RHCs (Additional file 1: Figure S1).
The DBS were stored at −20 °C in individual plastic bags containing desiccant until DNA extraction. DBS collected from February to September 2008 were initially stored at room temperature and subsequently at −20 °C. Parasite DNA was extracted using the Chelex© method from one dried blood spot . Briefly, the DBS were placed in 1.5 mL microcentrifuge tubes, 1 mL of 0.1% weight by volume saponin in 1× phosphate buffered saline (PBS) was added and the mixture was incubated for 10 min at room temperature. The tubes were centrifuged for 2 min at 14,000 rpm, the supernatant discarded and 1 mL of 1× PBS was added. The tubes were again centrifuged for 2 min at 14,000 rpm, the supernatant discarded and 150 μL of 2% weight by volume Chelex© solution and 50 μL of DNase free water were added and the tubes were boiled for 8 min. The tubes were then centrifuged for 1 min at 14,000 rpm and approximately 150 μL of DNA was stored at −20 °C.
Infection was confirmed using a Plasmodium multi-species nested PCR assay as previously described [18, 23] targeting a 815 base-pair segment of the mitochondrial cytochrome b gene (cytb). Briefly, in the primary PCR step, 6 μL of DNA extract was added to 0.2 mL tubes containing a 19 μL reaction mix made up of DNase free water and final concentrations of 200 μM dNTPs, 1× ThermoPol PCR buffer (New England Biolabs, Ipswich, MA, USA), 1 μM forward and 1 μM reverse primers and DNA Taq polymerase (New England Biolabs, Ipswich, MA) in 25 μL reaction volume. In the nested PCR step, 3 μL of the primary PCR product was added to 0.2 mL PCR tubes containing 22 μL of reaction mix with DNase free water and the final concentrations of dNTPs, Thermopol buffer, forward and reverse primers, and Taq DNA polymerase in a 25 μL reaction mix as described above. No template controls were included in each experiment and reactions were run in a Techne™ TC-412 thermo cycler (Bibby Scientific Limited, Staffordshire, UK). Amplified product was detected by electrophoresis on 1% agarose gel and viewed under UV light.
The 24 SNP molecular barcode assays were run using a TaqMan protocol at the Macha Research Trust laboratory in Zambia . DNA was extracted from a second dried blood spot for samples positive by nested PCR using the Chelex© method. Due to low parasite DNA concentrations, samples were pre-amplified prior to performing the 24 SNP molecular barcode assay . The pre-amplification step was done by adding 5 μL of DNA from each sample to 10 μL of TaqMan pre-amplification master mix, and 0.2× pooled assay mixture, containing of forward and reverse primers for each of the 24 SNPs as previously described .
Pre-amplified samples were diluted 1:20 with TE buffer prior to running the 24 SNP molecular barcoding assay as previously described . For each of the 24 SNP assays, 2.0 μL of pre-amplified sample DNA was added to 10 μL of TaqMan master mix, 7.5 μL distilled water, and 0.5 μL TaqMan commercially available primer and probe assay mixture. For each assay, three known positive controls and two negative, no-template controls were run. Positive controls consisted of DNA samples from P. falciparum strains obtained from MR4 (now known as BEI Resources, Manassas, VA, USA) with known haplotypes for all 24 SNPs. Typically, 12 SNP assays were run for five samples at a time with controls on a 96-well plate. The assays were run on the Applied Biosystems StepOnePlus™ (Thermo Scientific, Waltham, MA, USA), and the Roche LightCycler 480™ (Roche Diagnostics Corporation, Indianapolis, IN, USA) real time PCR systems.
Parasites sampled from peripheral blood are those in the haploid intra-erythrocytic stage of their lifecycle. . Thus, SNP calls were made automatically based on allelic discrimination plots using software programs accompanying the real time PCR systems as one of the two alleles or mixed [25, 26]. In cases where SNP calls could not be made automatically, determination was made manually by the study investigators. Otherwise the SNP call was classified as failed. Samples with failed SNP calls were repeated up to three times. If a sample failed on all repeated assays, the SNP was treated as missing data.
Haplotypes, composed of the concatenated 24 SNPs, were aligned using the ClustalW method in MEGA6 software, with mixed calls represented as a third possible allele . Phylogenetic analyses were conducted using the maximum likelihood methods for actively and passively detected malaria infections separately, grouped by malaria transmission season, to explore temporal phylogenetic clustering of haplotypes. Haplotypes for the combined dataset were aligned using the same method and a phylogenetic tree was constructed using the maximum likelihood method.
Parasite genetic complexity and diversity
Parasite genetic complexity was determined by the number of mixed calls at each of the 24 SNPs. Samples with four or more mixed calls were categorized as polyclonal infections . Temporal trends were graphed and variation in the proportion of polyclonal infections was determined using the Wilcoxon rank-sum test. Samples with more than half missing data were excluded (eight samples from active surveillance and seven samples from passive surveillance).
Parasite genetic diversity was evaluated by determining divergence from the most common barcode for each transmission season stratified by whether the parasite was identified through passive or active surveillance. The most common barcode for each season was determined by allelic frequency at each SNP to determine the nucleotide diversity. A modified ‘SNP π’, was developed to account for missing data and mixed allele calls, and was used to measure the seasonal parasite genetic divergence . Details of the methods for calculating the modified ‘SNP π’ are described in the Additional file 2: Appendix.
Data on age and sex analysed for passively detected cases. Age, sex, symptoms of malaria, and RDT and microscopy results were analysed for actively detected cases. Median ages between actively and passively detected infected individuals and across transmission seasons were compared using a Kruskal–Wallis test. Sex differences between groups and across transmission seasons were compared using a chi2 test. RDT and microscopy results across transmission seasons were compared for actively detected cases using a chi2 test.
The 24 SNP molecular barcode was performed on 72 samples from actively detected malaria cases collected from 2008 to 2013 and 46 samples from passively detected malaria cases collected from 2012 to 2015. Households of actively detected cases and the RHC from which the passively detected cases were identified came from the same underlying source population (Fig. 1). The median age was 13 years (IQR = 7–30) for infected individuals identified through passive case detection and 14 years (IQR = 11–21) for those identified through active case detection (p = 0.63) (Table 1). There were no statistically significant differences in sex between the two study populations (p = 0.34), and no statistically significant differences in median age or sex were observed in passively (p = 0.15, p = 0.63 respectively) or actively (p = 0.79, p = 0.96 respectively) detected cases across the malaria transmission seasons (Table 1). While all actively detected cases were identified by PCR, no statistically significant differences were observed in microscopy results across the malaria transmission seasons (p = 0.07), with most infected individuals detected by active case detection negative by microscopy (Table 1). Statistically significant differences were observed in the RDT results across malaria transmission seasons (p = 0.001), with the proportion of RDT positive infections variable, but all seasons had a higher proportion of RDT negative infections among individuals identified through active surveillance (Table 1). Statistically significant differences were observed in the proportion of individuals reporting symptoms of malaria for each season (p = 0.02), with trends in symptomatic infections decreasing from the 2008 season through the 2010–2011 season, then subsequently increasing (Table 1). For both actively and passively detected malaria cases, no identical molecular barcodes were detected and all haplotypes were unique.
The phylogenetic tree showed evidence of clustering of haplotypes from passively detected cases during the 2012–2013 annual malaria transmission season along with some mixing during the other two seasons (Fig. 2a). However, there was no evidence of phylogenetic clustering during the 2013–2014 and 2014–2015 annual malaria transmission seasons (Fig. 2a). For actively detected cases, the phylogenetic tree showed no evidence of phylogenetic clustering during any of the malaria transmission seasons (Fig. 2b). In a combined phylogenetic tree, parasite haplotypes from passively and actively detected cases were largely divergent (Fig. 3), suggesting two genetically distinct parasite populations (Fig. 3).
Genetic complexity and diversity
Genetic complexity was approximated by calculating the proportion of polyclonal infections in a transmission season. The proportion of polyclonal infections was consistently low for passively detected cases in all seasons, with no polyclonal infections identified during the 2012–2013 and 2014–2015 annual transmission seasons (Table 2; Fig. 4a). In contrast, the proportion of polyclonal infections among actively detected cases was consistently high for all malaria transmission seasons (Fig. 4b). The last two seasons were combined due to the low numbers of cases but all infections were polyclonal (Table 2; Fig. 4b). No seasonal trends in genetic complexity were identified (Fig. 4b).
The mean genetic divergence for passively detected cases was 34.5% for the 2012–2013 season, 37.8% for the 2013–2014 season, and 30.8% for the 2014–2015 season (Table 2; Fig. 5a). The mean genetic divergence for actively detected cases was 22.3% in the 2008 season and 29.0% in the 2008–2009 season (Table 2; Fig. 5b). This decreased to 9.9% in the combined 2012–2014 seasons (Table 2; Fig. 5b). Overall, the genetic divergence remained high among passively detected cases but decreased for actively detected cases (Fig. 5b).
Distinct parasite populations were found in individuals identified through passive and active surveillance in a region of declining malaria transmission in southern Zambia using a SNP-based molecular barcode for P. falciparum. These results suggest that the population of mostly asymptomatic, actively detected infections may not significantly contribute to ongoing transmission in this setting and may not be the primary source of infection leading to clinical malaria cases.
The molecular barcode can be used to document parasite diversity over time by providing information on the number of unique barcodes present each year, the persistence of unique barcodes between years, and the genetic divergence measured by nucleotide diversity within and between years. The molecular barcode determines whether the major allele, minor allele, or a mixture of major and minor alleles is present in the haploid blood stage of the parasite infection for each of the 24 SNPs. The resulting barcode represents the parasite’s haplotype and can be used to compare the genetic relatedness between individual (monoclonal and polyclonal) infections and populations of parasites infecting different communities . Variation arises as a result of genetic recombination and outcrossing during the sexual stage of infection in the mosquito vector when gametocytes combine . SNPs with a mixture of major and minor alleles are referred to as mixed (polyclonal) infections, with higher numbers of mixed infections representing higher genetic complexity. The frequency of mixed infections approximates the level of genetic complexity and provides information about the burden of infection due to genetically unique parasites .
Typically, as malaria transmission declines, a genetic bottleneck is created, reducing opportunities for outcrossing in the mosquito midgut and leading to reduced parasite diversity among passively detected, symptomatic malaria cases [13, 28]. This decline in diversity is accompanied by a decline in the complexity of infection, as fewer unique parasites circulate in the population [13, 28]. However, parasite diversity among passively detected, symptomatic infections did not decline in southern Zambia and the complexity of infection was low. As malaria transmission was consistently low between the 2012 and 2015 seasons, the sustained, high parasite diversity in symptomatic infections may be evidence of imported malaria with subsequent local transmission among individuals susceptible to clinical malaria.
Passively detected cases also had fewer complex infections, with a low proportion of polyclonal infections, consistent with recent infection with single parasite clones. In contrast, actively detected cases had more complex infections, with a high proportion of polyclonal infections. Between the 2008 and 2009–2010 transmission seasons there was an indication of a non-statistically significant, decreasing trend in the proportion of polyclonal infections; however, this proportion subsequently increased and remained high, indicating that the actively detected population of mostly asymptomatic and sub-patent infections was composed of many different parasite clones. While transmission declined, these individuals may have harboured parasites acquired over time and maintained at low levels of parasitaemia. Chronic infections have been shown to persist for up to a decade [29,30,31,32].
Parasite genetic diversity remained relatively high and constant for passively detected cases throughout the observed malaria transmission seasons. This finding, along with the phylogenetic separation from actively detected cases and the paucity of detectable polyclonal infections, suggests these passively detected cases may represented recent infections with parasites imported into the study area. In contrast, parasite genetic diversity was lower among actively detected cases, consistent with a chronically infected population. The decreasing parasite diversity among the actively detected cases over the study period is consistent with loss of parasites from this chronically infected reservoir as malaria transmission declined.
Previous studies using this genotyping approach were reported from Senegal and Malawi. The molecular barcode was used in Senegal to determine the clonal and epidemic expansion of passively detected clinical P. falciparum infections . After enhanced deployment of ITNs and use of RDTs and ACT between 2005 and 2011, reductions in parasite genetic diversity and complexity were observed . Between 2006 and 2013, the molecular barcode was used to determine the decline in malaria transmission among passively detected infections, with re-introduction in 2012 . In Malawi, the molecular barcode was used to determine differences in the complexity of infection between severe clinical malaria cases and cerebral malaria cases among children younger than 5 years of age . Children with cerebral malaria had less complex infections than those with severe malaria . To the authors’ knowledge, this is the first time this molecular barcode was used to determine the genetic relatedness between actively and passively detected infections.
Analyses of spatial hotspots of asymptomatic parasitemia and clinical malaria cases have been previously published from low and moderate transmission settings in Kenya and Mali [32,33,34,35]. In the Kenya studies, no spatial overlap was identified between hotspots of asymptomatic malaria infections and clinical malaria infections [33, 34]. This spatial disaggregation supports the hypothesis that the asymptomatic reservoir may not be the major source of infection leading to clinical malaria in low transmission settings. The parasite prevalence in the Macha study area is lower than reported in the Kenya studies. The genetic separation of parasites infecting the asymptomatic and symptomatic populations identified in the current study provides additional support that as areas approach elimination the chronically infected population may not contribute substantially to ongoing transmission and a high proportion of clinical cases are likely due to seasonal importation. This has important implications for the need for resource-intensive efforts to identify and treat the chronically-infected, asymptomatic reservoir through reactive test-and-treat strategies, focal or mass drug administration, or using highly sensitivity diagnostic tests.
There were several major limitations of these analyses. The first was that passive and active case detection did not overlap in time except for one malaria transmission season. However, given the absence of seasonal phylogenetic clustering for both passively and actively detected cases, and because all samples came from the same underlying source population, the combined phylogenetic tree is likely an accurate representation of the phylogenetic relationships between the two parasite populations. Additionally, when constructing the combined phylogenetic tree, samples were indexed by season so within season relationships between the two populations would not be masked. The second limitation was the high number of mixed infections from which single haplotypes were difficult to resolve. To account for this, a third allele was used to indicate mixed infections in calculating the phylogenetic relatedness and mixed infections were accounted for when calculating genetic divergence [13, 14]. The third limitation was that the sample size was relatively small, although this is to be expected in a pre-elimination setting. The fourth limitation was the low level of parasitaemia, also an expected finding in a pre-elimination setting. Pre-amplification of parasite DNA was performed to increase the molecular barcode yield but potentially introduced bias in the allele frequency. However, pre-amplification sites were barcode specific to reduce amplification in non-specific parasite DNA, and all samples were pre-amplified regardless of level of parasitaemia . After pre-amplification, assays with greater than three failures were treated as missing data and samples with over half of data missing were dropped from analyses. However, missing data may have been informative as these samples had low levels of parasitaemia.
Distinct parasite populations were found among infected individuals identified through active and passive surveillance in a low transmission setting in southern Zambia. The actively detected population represented a chronically infected, reservoir that may not contribute substantially to ongoing transmission. As parasite prevalence and diversity within the actively detected cases declined over time, resource-intensive efforts to identify and treat the chronically infected reservoir through reactive case detection, or focal or mass drug administration, may not be necessary to eliminate malaria in this setting.
O’Meara WP, Mangeni JN, Steketee R, Greenwood B. Changes in the burden of malaria in sub-Saharan Africa. Lancet Infect Dis. 2010;10:545–55.
Cotter C, Sturrock HJ, 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.
World Health Organization. Disease surveillance for malaria elimination: an operational manual. Geneva: World Health Organization; 2012.
malERA Consultative Group on Monitoring Evaluation, and Surveillance. A research agenda for malaria eradication: monitoring, evaluation, and surveillance. PLoS Med. 2011;8:1000400.
Hay SI, Smith DL, Snow RW. Measuring malaria endemicity from intense to interrupted transmission. Lancet Infect Dis. 2008;8:369–78.
Bousema T, Okell L, Felger I, Drakeley C. Asymptomatic malaria infections: detectability, transmissibility and public health relevance. Nat Rev Microbiol. 2014;12:833–40.
Imwong M, Nguyen TN, Tripura R, Peto TJ, Lee SJ, Lwin KM, et al. The epidemiology of subclinical malaria infections in South–East Asia: findings from cross-sectional surveys in Thailand–Myanmar border areas, Cambodia, and Vietnam. Malar J. 2015;14:381.
Harris I, Sharrock WW, Bain LM, Gray KA, Bobogare A, Boaz L, et al. A large proportion of asymptomatic Plasmodium infections with low and sub-microscopic parasite densities in the low transmission setting of Temotu Province, Solomon Islands: challenges for malaria diagnostics in an elimination setting. Malar J. 2010;9:254.
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.
Karl S, Gurarie D, Zimmerman PA, King CH, St Pierre TG, Davis TM. A sub-microscopic gametocyte reservoir can sustain malaria transmission. PLoS ONE. 2011;6:e20805.
Daniels R, Volkman SK, Milner DA, Mahesh N, Neafsey DE, Park DJ, et al. A general SNP-based molecular barcode for Plasmodium falciparum identification and tracking. Malar J. 2008;7:223.
Obaldia N 3rd, Baro NK, Calzada JE, Santamaria AM, Daniels R, Wong W, et al. Clonal outbreak of Plasmodium falciparum infection in eastern Panama. J Infect Dis. 2015;211:1087–96.
Daniels R, Chang HH, Sene PD, Park DC, Neafsey DE, Schaffner SF, et al. Genetic surveillance detects both clonal and epidemic transmission of malaria following enhanced intervention in Senegal. PLoS ONE. 2013;8:e60780.
Neafsey DE, Schaffner SF, Volkman SK, Park D, Montgomery P, Milner DA Jr, et al. Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence. Genome Biol. 2008;9:R171.
Pinchoff J, Hamapumbu H, Kobayashi T, Simubali L, Stevenson JC, Norris DE, et al. Factors associated with sustained use of long-lasting insecticide-treated nets following a reduction in malaria transmission in Southern Zambia. Am J Trop Med Hyg. 2015;93:954–60.
Kent R, Thuma P, Mharakurwa S, Norris D. 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.
Fornadel CM, Norris LC, Glass GE, Norris DE. Analysis of Anopheles arabiensis blood feeding behavior in southern Zambia during the two years after introduction of insecticide-treated bed nets. Am J Trop Med Hyg. 2010;83:848–53.
Laban NM, Kobayashi T, Hamapumbu 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.
Steketee R, Sipilanyambi N, Chimumbwa J, Banda J, Mohamed A, Miller J, et al. National malaria control and scaling up for impact: the Zambian experience through 2006. Am J Trop Med Hyg. 2008;79:45–52.
MOH. Zambia National Malaria Indicator Survey 2012. Lusaka: Ministry of Health Government of the Republic of Zambia; 2012. p. 2012.
Moss W, Hamapumbu 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. Malar J. 2011;10:163.
Shiff CJ, Stoyanov C, Choobwe C, Kamanga A, Mukonka VM. Measuring malaria by passive case detection: a new perspective based on Zambian experience. Malar J. 2013;12:120.
Steenkeste N, Incardona S, Chy S, Duval L, Ekala MT, Lim P, et al. Towards high-throughput molecular detection of Plasmodium: new approaches and molecular markers. Malar J. 2009;8:86.
Mharakurwa S, Daniels R, Scott A, Wirth DF, Thuma P, Volkman SK. Pre-amplification methods for tracking low-grade Plasmodium falciparum populations during scaled-up interventions in Southern Zambia. Malar J. 2014;13:89.
Tojo Y, Asahina J, Miyashita Y, Takahashi M, Matsumoto N, Hasegawa S, et al. Development of an automation system for single nucleotide polymorphisms genotyping using bio-strand, a new three-dimensional microarray. J Biosci Bioeng. 2005;99:120–4.
Gut IG. Automation in genotyping of single nucleotide polymorphisms. Hum Mutat. 2001;17:475–92.
Milner DA Jr, Vareta J, Valim C, Montgomery J, Daniels RF, Volkman SK, et al. Human cerebral malaria and Plasmodium falciparum genotypes in Malawi. Malar J. 2012;11:35.
Daniels RF, Schaffner SF, Wenger EA, Proctor JL, Chang HH, Wong W, et al. Modeling malaria genomics reveals transmission decline and rebound in Senegal. Proc Natl Acad Sci USA. 2015;112:7067–72.
Kitchen AD, Chiodini PL, Tossell J. Detection of malarial DNA in blood donors—evidence of persistent infection. Vox Sang. 2014;107:123–31.
Assennato SM, Berzuini A, Foglieni B, Spreafico M, Allain JP, Prati D. Plasmodium genome in blood donors at risk for malaria after several years of residence in Italy. Transfusion. 2014;54:2419–24.
D’Ortenzio E, Godineau N, Fontanet A, Houze S, Bouchaud O, Matheron S, et al. Prolonged Plasmodium falciparum infection in immigrants, Paris. Emerg Infect Dis. 2008;14:323–6.
Ashley EA, White NJ. The duration of Plasmodium falciparum infections. Malar J. 2014;13:500.
Kangoye DT, Noor A, Midega J, Mwongeli J, Mkabili D, Mogeni P, et al. Malaria hotspots defined by clinical malaria, asymptomatic carriage, PCR and vector numbers in a low transmission area on the Kenyan Coast. Malar J. 2016;15:213.
Zhou G, Afrane YA, Malla S, Githeko AK, Yan G. Active case surveillance, passive case surveillance and asymptomatic malaria parasite screening illustrate different age distribution, spatial clustering and seasonality in western Kenya. Malar J. 2015;14:41.
Sissoko MS, van den Hoogen LL, Samake Y, Tapily A, Diarra AZ, Coulibaly M, et al. Spatial patterns of Plasmodium falciparum clinical incidence, asymptomatic parasite carriage and anopheles density in two villages in Mali. Am J Trop Med Hyg. 2015;93:790–7.
KMS conceived of the study, coordinated and carried out the analyses, and drafted the manuscript. BK conducted the laboratory assays and participated in manuscript preparation. TK participated in the design, analyses and manuscript preparation. MNSS conducted the laboratory assays and participated in manuscript preparation. SM participated in the study design and analysis. GC participated in the analysis and manuscript preparation. DEN participated in the design, analysis, and manuscript preparation. JCS participated in the design of the study and manuscript preparation. PET participated in the design of the study and manuscript preparation. WJM conceived of the study, participated in its design, and drafted the manuscript. All authors read and approved the final manuscript.
We thank the study participants for their participation and commitment to this study, and the health workers that collected samples and surveys as part of passive case detection. We also thank the Macha Research Trust field team for their work on study activities and data collection. We thank Julia Pringle for advice and input on population genetics analytic methods. Lastly, we thank Andre Hackman for his work on data management for this study.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated during and/or analysed during the current study are not publicly available to respect the wishes of the Ministry of Health of the Government of Zambia, but are available from the corresponding author on reasonable request.
The study was approved by the Institutional Review Board at the Johns Hopkins Bloomberg School of Public Health and the Review Board and the Tropical Diseases Research Centre Ethics Review Committee in Ndola, Zambia. Informed consent was obtained from all participants.
This work was supported by the Johns Hopkins Malaria Research Institute, the Bloomberg Family Foundation, and the Division of Microbiology and Infectious Diseases, National Institutes of Allergy and Infectious Diseases, National Institutes of Health as part of the International Centers of Excellence for Malaria Research (U19 AI089680).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 2. The modified ‘SNP π’ used to account for missing data and mixed allele calls was calculated as follows.
About this article
Cite this article
Searle, K.M., Katowa, B., Kobayashi, T. 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 16, 154 (2017). https://doi.org/10.1186/s12936-017-1810-3