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

BACKGROUND: Over a decade of vector control by indoor residual spraying (IRS) and long-lasting insecticidal nets (LLINs) distribution on the mainland, and only LLINs on islands had a minimal impact on disease burden in Nchelenge district, northern Zambia. Anopheles funestus and Anopheles gambiae are vectors known only from the mainland. Understanding vector bionomics in the district is necessary for planning and targeting effective vector control. This study aimed to provide information on abundance, seasonality, and Plasmodium falciparum sporozoite infectivity of malaria vectors in Nchelenge, including islands. METHODS: Mosquitoes were collected in 192 CDC indoor light traps set in 56 households between January 2015 and January 2016. Morphological and molecular species identifications and P. falciparum circumsporoites by ELISA were performed. Mosquito counts and relative abundances from the islands and mainland were compared, and household factors associated with vector counts were determined. RESULTS: A total of 5888 anophelines were collected during the study. Of these, 5,704 were female Anopheles funestus sensu lato (s.l.) and 248 female An. gambiae s.l. The highest proportion of An. funestus (n = 4090) was from Chisenga Island and An. gambiae (n = 174) was from Kilwa Island. The highest estimated counts per trap for An. funestus s.l. were from Chisenga island, (89.9, p < 0.001) and from the dry season (78.6, p < 001). For An. gambiae the highest counts per trap were from Kilwa island (3.1, p < 0.001) and the rainy season (2.5, p = 0.007). The highest estimated annual entomological inoculation rate was from Chisenga Island with 91.62 ib/p/y followed by Kilwa Island with 29.77 ib/p/yr, and then Mainland with 19.97 ib/p/yr. CONCLUSIONS: There was varied species abundance and malaria transmission risk across sites and seasons. The risk of malaria transmission was perennial and higher on the islands. The minimal impact of vector control efforts on the mainland was evident, but limited overall. Vector control intervention coverage with effective tools needs to be extended to the islands to effectively control malaria transmission in Nchelenge district.

Temporal and spatial analysis of Plasmodium falciparum genomics reveals patterns of parasite connectivity in a low-transmission district in Southern Province, Zambia.

BACKGROUND: Understanding temporal and spatial dynamics of malaria transmission will help to inform effective interventions and strategies in regions approaching elimination. Parasite genomics are increasingly used to monitor epidemiologic trends, including assessing residual transmission across seasons and importation of malaria into these regions. METHODS: In a low and seasonal transmission setting of southern Zambia, a total of 441 Plasmodium falciparum samples collected from 8 neighbouring health centres between 2012 and 2018 were genotyped using molecular inversion probes (MIPs n = 1793) targeting a total of 1832 neutral and geographically informative SNPs distributed across the parasite genome. After filtering for quality and missingness, 302 samples and 1410 SNPs were retained and used for downstream population genomic analyses. RESULTS: The analyses revealed most (67%, n = 202) infections harboured one clone (monogenomic) with some variation at local level suggesting low, but heterogenous malaria transmission. Relatedness identity-by-descent (IBD) analysis revealed variable distribution of IBD segments across the genome and 6% of pairs were highly-related (IBD ≥ 0.25). Some of the highly-related parasite populations persisted across multiple seasons, suggesting that persistence of malaria in this low-transmission region is fueled by parasites "seeding" across the dry season. For recent years, clusters of clonal parasites were identified that were dissimilar to the general parasite population, suggesting parasite populations were increasingly fragmented at small spatial scales due to intensified control efforts. Clustering analysis using PCA and t-SNE showed a lack of substantial parasite population structure. CONCLUSION: Leveraging both genomic and epidemiological data provided comprehensive picture of fluctuations in parasite populations in this pre-elimination setting of southern Zambia over 7 years.

A Bayesian hierarchical model for signal extraction from protein microarrays.

Protein microarrays are a promising technology that measure protein levels in serum or plasma samples. Due to their high technical variability and high variation in protein levels across serum samples in any population, directly answering biological questions of interest using protein microarray measurements is challenging. Analyzing preprocessed data and within-sample ranks of protein levels can mitigate the impact of between-sample variation. As for any analysis, ranks are sensitive to preprocessing, but loss function based ranks that accommodate major structural relations and components of uncertainty are very effective. Bayesian modeling with full posterior distributions for quantities of interest produce the most effective ranks. Such Bayesian models have been developed for other assays, for example, DNA microarrays, but modeling assumptions for these assays are not appropriate for protein microarrays. Consequently, we develop and evaluate a Bayesian model to extract the full posterior distribution of normalized protein levels and associated ranks for protein microarrays, and show that it fits well to data from two studies that use protein microarrays produced by different manufacturing processes. We validate the model via simulation and demonstrate the downstream impact of using estimates from this model to obtain optimal ranks.

A pre-processing pipeline to quantify, visualize, and reduce technical variation in protein microarray studies.

Technical variation, or variation from non-biological sources, is present in most laboratory assays. Correcting for this variation enables analysts to extract a biological signal that informs questions of interest. However, each assay has different sources and levels of technical variation, and the choice of correction methods can impact downstream analyses. Compared to similar assays such as DNA microarrays, relatively few methods have been developed and evaluated for protein microarrays, a versatile tool for measuring levels of various proteins in serum samples. Here, we propose a pre-processing pipeline to correct for some common sources of technical variation in protein microarrays. The pipeline builds upon an existing normalization method by using controls to reduce technical variation. We evaluate our method using data from two protein microarray studies and by simulation. We demonstrate that pre-processing choices impact the fluorescent-intensity based ranks of proteins, which in turn, impact downstream analysis.

Expanded geographic distribution and host preference of Anopheles gibbinsi (Anopheles species 6) in northern Zambia.

BACKGROUND: Nchelenge District in northern Zambia suffers from holoendemic malaria transmission despite a decade of yearly indoor residual spraying (IRS) and insecticide-treated net (ITN) distributions. One hypothesis for this lack of impact is that some vectors in the area may forage in the early evening or outdoors. Anopheles gibbinsi specimens were identified in early evening mosquito collections performed in this study area, and further insight was gleaned into this taxon, including characterizing its genetic identity, feeding preferences, and potential role as a malaria vector. METHODS: Mosquitoes were collected in July and August 2019 by CDC light traps in Nchelenge District in indoor sitting rooms, outdoor gathering spaces, and animal pens from 16:00-22:00. Host detection by PCR, COI and ITS2 PCR, and circumsporozoite (CSP) ELISA were performed on all samples morphologically identified as An. gibbinsi, and a subset of specimens were selected for COI and ITS2 sequencing. To determine risk factors for increased abundance of An. gibbinsi, a negative binomial generalized linear mixed-effects model was performed with household-level variables of interest. RESULTS: Comparison of COI and ITS2 An. gibbinsi reference sequences to the NCBI database revealed > 99% identity to "Anopheles sp. 6" from Kenya. More than 97% of specimens were morphologically and molecularly consistent with An. gibbinsi. Specimens were primarily collected in animal pen traps (59.2%), followed by traps outdoors near where humans gather (24.3%), and traps set indoors (16.5%). Host DNA detection revealed a high propensity for goats, but 5% of specimens with detected host DNA had fed on humans. No specimens were positive for Plasmodium falciparum sporozoites. Animal pens and inland households > 3 km from Lake Mweru were both associated with increased An. gibbinsi abundance. CONCLUSIONS: This is the first report of An. gibbinsi in Nchelenge District, Zambia. This study provided a species identity for unknown "An. sp. 6" in the NCBI database, which has been implicated in malaria transmission in Kenya. Composite data suggest that this species is largely zoophilic and exophilic, but comes into contact with humans and the malaria parasites they carry. This species should continue to be monitored in Zambia and neighbouring countries as a potential malaria vector.

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

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.

Genetic Evidence of Focal Plasmodium falciparum Transmission in a Pre-elimination Setting in Southern Province, Zambia.

BACKGROUND: Southern Province, Zambia has experienced a dramatic decline in Plasmodium falciparum malaria transmission in the past decade and is targeted for elimination. Zambia's National Malaria Elimination Program recommends reactive case detection (RCD) within 140 m of index households to enhance surveillance and eliminate remaining transmission foci. METHODS: To evaluate whether RCD captures local transmission, we genotyped 26 microsatellites from 106 samples collected from index (n = 27) and secondary (n = 79) cases detected through RCD in the Macha Hospital catchment area between January 2015 and April 2016. RESULTS: Participants from the same RCD event harbored more genetically related parasites than those from different RCD events, suggesting that RCD captures, at least in part, infections related through local transmission. Related parasites clustered in space and time, up to at least 250 m from index households. Spatial analysis identified a putative focal transmission hotspot. CONCLUSIONS: The current RCD strategy detects focal transmission events, although programmatic guidelines to screen within 140 m of index households may fail to capture all secondary cases. This study highlights the utility of parasite genetic data in assessing programmatic interventions, and similar approaches may be useful to malaria elimination programs seeking to tailor intervention strategies to the underlying transmission epidemiology.

The Impact of 3 Years of Targeted Indoor Residual Spraying With Pirimiphos-Methyl on Malaria Parasite Prevalence in a High-Transmission Area of Northern Zambia.

Malaria transmission in northern Zambia has increased in the past decade, despite malaria control activities. Evidence-based intervention strategies are needed to effectively reduce malaria transmission. Zambia's National Malaria Control Centre conducted targeted indoor residual spraying (IRS) in Nchelenge District, Luapula Province, from 2014 to 2016 using the organophosphate insecticide pirimiphos-methyl. An evaluation of the IRS campaign was conducted by the Southern Africa International Centers of Excellence for Malaria Research using actively detected malaria cases in bimonthly household surveys carried out from April 2012 to July 2017. Changes in malaria parasite prevalence after IRS were assessed by season using Poisson regression models with robust standard errors, controlling for clustering of participants in households and demographic, geographical, and climatological covariates. In targeted areas, parasite prevalence declined approximately 25% during the rainy season following IRS with pirimiphos-methyl but did not decline during the dry season or in the overall study area. Within targeted areas, parasite prevalence declined in unsprayed households, suggesting both direct and indirect effects of IRS. The moderate decrease in parasite prevalence within sprayed areas indicates that IRS with pirimiphos-methyl is an effective malaria control measure, but a more comprehensive package of interventions is needed to effectively reduce the malaria burden in this setting.

High Plasmodium falciparum genetic diversity and temporal stability despite control efforts in high transmission settings along the international border between Zambia and the Democratic Republic of the Congo.

BACKGROUND: While the utility of parasite genotyping for malaria elimination has been extensively documented in low to moderate transmission settings, it has been less well-characterized in holoendemic regions. High malaria burden settings have received renewed attention acknowledging their critical role in malaria elimination. Defining the role for parasite genomics in driving these high burden settings towards elimination will enhance future control programme planning. METHODS: Amplicon deep sequencing was used to characterize parasite population genetic diversity at polymorphic Plasmodium falciparum loci, Pfama1 and Pfcsp, at two timepoints in June-July 2016 and January-March 2017 in a high transmission region along the international border between Luapula Province, Zambia and Haut-Katanga Province, the Democratic Republic of the Congo (DRC). RESULTS: High genetic diversity was observed across both seasons and in both countries. No evidence of population structure was observed between parasite populations on either side of the border, suggesting that this region may be one contiguous transmission zone. Despite a decline in parasite prevalence at the sampling locations in Haut-Katanga Province, no genetic signatures of a population bottleneck were detected, suggesting that larger declines in transmission may be required to reduce parasite genetic diversity. Analysing rare variants may be a suitable alternative approach for detecting epidemiologically important genetic signatures in highly diverse populations; however, the challenge is distinguishing true signals from potential artifacts introduced by small sample sizes. CONCLUSIONS: Continuing to explore and document the utility of various parasite genotyping approaches for understanding malaria transmission in holoendemic settings will be valuable to future control and elimination programmes, empowering evidence-based selection of tools and methods to address pertinent questions, thus enabling more efficient resource allocation.

Laboratory diagnosis of Lyme borreliosis: Current state of the art and future perspectives.

This review is directed at physicians and laboratory personnel in private practice and clinics who treat and diagnose Lyme borreliosis (LB) in patients as part of their daily work. A major objective of this paper is to bring together background information on Borrelia (B.) burgdorferi sensu lato (s.l.) and basic clinical knowledge of LB, which is one of the most frequently reported vector-borne diseases in the Northern Hemisphere. The goal is to provide practical guidance for clinicians and for laboratory physicians, and scientists for a better understanding of current achievements and ongoing obstacles in the laboratory diagnosis of LB, an infectious disease that still remains one of the diagnostic chameleons of modern clinical medicine. Moreover, in bringing together current scientific information from guidelines, reviews, and original papers, this review provides recommendations for selecting the appropriate tests in relation to the patient's stage of disease to achieve effective, stage-related application of current direct and indirect laboratory methods for the detection of B. burgdorferi s.l. Additionally, the review aims to discuss the current state of the art concerning the diagnostic potential and limitations of the assays and test methods currently in use to optimize LB patient management and provide insight into the possible future prospects of this rapidly changing area of laboratory medicine.

Controlled release spatial repellent devices (CRDs) as novel tools against malaria transmission: a semi-field study in Macha, Zambia.

BACKGROUND: The emergence of mosquitoes that can avoid indoor-deployed interventions, such as treated bed nets and indoor residual spraying, threatens the mainstay of malaria control in Zambia. Furthermore, the requirement for high coverage of these tools poses operational challenges. Spatial repellents are being assessed to supplement these vector control tools, but limitations exist in the residual effect of the repellent and the need for external power or heat for diffusion of the volatiles. METHODS: A semi-field evaluation of a novel controlled release spatial repellent device (CRD) was conducted in Macha, Zambia. These devices emanate metofluthrin with no need for external power. Devices were deployed in huts within the semi-field system (SFS). Female Anopheles gambiae sensu stricto released within the SFS were trapped overnight by light traps and collected by aspiration the next morning inside and outside of huts to determine the extent of mosquito repellency and the impact on host-seeking and survival. Experiments studied the impact of number of devices as well as the presence of hut occupants. The study was complemented with numerical methods based on computational fluid dynamics to simulate spatial distribution of metofluthrin. RESULTS: Presence of CRDs was associated with significant reductions in indoor counts of mosquitoes, regardless of whether huts were occupied or not. Repellency ranged from 15 to 60% compared to huts with no devices. Reducing the number of devices from 16 to 4 had little impact on repellency. When huts were occupied, indoor mosquito host-seeking was higher in the presence of CRDs, whilst survival was significantly reduced. CONCLUSIONS: This study demonstrated that deployment of as few as four CRDs within a hut was associated with reduced indoor mosquito densities. As would be expected, presence of occupants within huts, resulted in greater indoor catches (both with and without devices). The increased indoor mosquito host-seeking and mortality in huts when devices were present may be explained by the excito-repellency activity of metofluthrin. These semi-field experiments provide preliminary data on the utility of CRD spatial repellents to reduce indoor densities of An. gambiae mosquitoes. Studies will further investigate the impact of CRDs on mosquito behaviour as well as epidemiological protective efficacy.

Malaria knowledge and bed net use in three transmission settings in southern Africa.

BACKGROUND: Insecticide-treated nets (ITNs) reduce malaria morbidity and mortality in endemic areas. Despite increasing availability, the use of ITNs remains limited in some settings. Poor malaria knowledge is a barrier to the widespread use of ITNs. The goal of this study was to assess the levels of malaria knowledge and evaluate factors associated with bed net use among individuals residing in three regions of southern Africa with different levels of malaria transmission and control. METHODS: A cross-sectional study was conducted on a sample of 7535 residents recruited from 2066 households in Mutasa District, Zimbabwe (seasonal malaria transmission), Choma District, Zambia (low transmission) and Nchelenge District, Zambia (high transmission), between March 2012 and March 2017. A standardized questionnaire was used to collect data on demographics, malaria-related knowledge and use of preventive measures. Multivariate logistic regression analyses were used to assess determinants of bed net use. RESULTS: Most of the 3836 adult participants correctly linked mosquito bites to malaria (85.0%), mentioned at least one malaria symptom (95.5%) and knew of the benefit of sleeping under an ITN. Bed net ownership and use were highest in Choma and Nchelenge Districts and lowest in Mutasa District. In multivariate analyses, knowledge of ITNs was associated with a 30-40% increased likelihood of bed net use after adjusting for potential confounders across all sites. Other factors significantly associated with bed net use were age, household size and socioeconomic status, although the direction, strength and size of association varied by study site. Importantly, participants aged 5-14 years had reduced odds of sleeping under a bed net compared to children younger than 5 years. CONCLUSION: Relevant knowledge of ITNs translated into the expected preventive behaviour of sleeping under a bed net, underscoring the need for continued health messaging on malaria prevention. The implementation and delivery of malaria control and elimination interventions needs to consider socioeconomic equity gaps, and target school-age children to ensure access to and improve utilization of ITNs.

Distinct parasite populations infect individuals identified through passive and active case detection in a region of declining malaria transmission in southern Zambia.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

An Operational Framework for Insecticide Resistance Management Planning.

Arthropod vectors transmit organisms that cause many emerging and reemerging diseases, and their control is reliant mainly on the use of chemical insecticides. Only a few classes of insecticides are available for public health use, and the increased spread of insecticide resistance is a major threat to sustainable disease control. The primary strategy for mitigating the detrimental effects of insecticide resistance is the development of an insecticide resistance management plan. However, few examples exist to show how to implement such plans programmatically. We describe the formulation and implementation of a resistance management plan for mosquito vectors of human disease in Zambia. We also discuss challenges, steps taken to address the challenges, and directions for the future.

Underestimation of foraging behaviour by standard field methods in malaria vector mosquitoes in southern Africa.

BACKGROUND: Defining the anopheline mosquito vectors and their foraging behaviour in malaria endemic areas is crucial for disease control and surveillance. The standard protocol for molecular identification of host blood meals in mosquitoes is to morphologically identify fed mosquitoes and then perform polymerase chain reaction (PCR), precipitin tests, or ELISA assays. The purpose of this study was to determine the extent to which the feeding rate and human blood indices (HBIs) of malaria vectors were underestimated when molecular confirmation by PCR was performed on both visually fed and unfed mosquitoes. METHODS: In association with the Southern Africa International Centers of Excellence in Malaria Research (ICEMR), mosquito collections were performed at three sites: Choma district in southern Zambia, Nchelenge district in northern Zambia, and Mutasa district in eastern Zimbabwe. All anophelines were classified visually as fed or unfed, and tested for blood meal species using PCR methods. The HBIs of visually fed mosquitoes were compared to the HBIs of overall PCR confirmed fed mosquitoes by Pearson's Chi-Square Test of Independence. RESULTS: The mosquito collections consisted of Anopheles arabiensis from Choma, Anopheles funestus s.s., Anopheles gambiae s.s. and Anopheles leesoni from Nchelenge, and An. funestus s.s. and An. leesoni from Mutasa. The malaria vectors at all three sites had large human blood indices (HBI) suggesting high anthropophily. When only visually fed mosquitoes tested by PCR for blood meal species were compared to testing those classified as both visually fed and unfed mosquitoes, it was found that the proportion blooded was underestimated by up to 18.7%. For most Anopheles species at each site, there was a statistically significant relationship (P < 0.05) between the HBIs of visually fed mosquitoes and that of the overall PCR confirmed fed mosquitoes. CONCLUSION: The impact on HBI of analysing both visually fed and unfed mosquitoes varied from site to site. This discrepancy may be due to partial blood feeding behaviour by mosquitoes, digestion of blood meals, sample condition, and/or expertise of entomology field staff. It is important to perform molecular testing on all mosquitoes to accurately characterize vector feeding behaviour and develop interventions in malaria endemic areas.

High burden of malaria following scale-up of control interventions in Nchelenge District, Luapula Province, Zambia.

BACKGROUND: Malaria control interventions have been scaled-up in Zambia in conjunction with a malaria surveillance system. Although substantial progress has been achieved in reducing morbidity and mortality, national and local information demonstrated marked heterogeneity in the impact of malaria control across the country. This study reports the high burden of malaria in Nchelenge District, Luapula Province, Zambia from 2006 to 2012 after seven years of control measures. METHODS: Yearly aggregated information on cases of malaria, malaria deaths, use of malaria diagnostics, and malaria control interventions from 2006 to 2012 were obtained from the Nchelenge District Health Office. Trends in the number of malaria cases, methods of diagnosis, malaria positivity rate among pregnant women, and intervention coverage were analysed using descriptive statistics. RESULTS: Malaria prevalence remained high, increasing from 38% in 2006 to 53% in 2012. Increasing numbers of cases of severe malaria were reported until 2010. Intense seasonal malaria transmission was observed with seasonal declines in the number of cases between April and August, although malaria transmission continued throughout the year. Clinical diagnosis without accompanying confirmation declined from 95% in 2006 to 35% in 2012. Intervention coverage with long-lasting insecticide-treated nets and indoor residual spraying increased from 2006 to 2012. CONCLUSIONS: Despite high coverage with vector control interventions, the burden of malaria in Nchelenge District, Zambia remained high. The high parasite prevalence could accurately reflect the true burden, perhaps in part as a consequence of population movement, or improved access to care and case reporting. Quality information at fine spatial scales will be critical for targeting effective interventions and measurement of progress.

A metagenomic analysis of the phase 2 Anopheles gambiae 1000 genomes dataset reveals a wide diversity of cobionts associated with field collected mosquitoes.

The Anopheles gambiae 1000 Genomes (Ag1000G) Consortium previously utilized deep sequencing methods to catalogue genetic diversity across African An. gambiae populations. We analyzed the complete datasets of 1142 individually sequenced mosquitoes through Microsoft Premonition's Bayesian mixture model based (BMM) metagenomics pipeline. All specimens were confirmed as either An. gambiae sensu stricto (s.s.) or An. coluzzii with a high degree of confidence ( > 98% identity to reference). Homo sapiens DNA was identified in all specimens indicating contamination may have occurred either at the time of specimen collection, preparation and/or sequencing. We found evidence of vertebrate hosts in 162 specimens. 59 specimens contained validated Plasmodium falciparum reads. Human hepatitis B and primate erythroparvovirus-1 viral sequences were identified in fifteen and three mosquito specimens, respectively. 478 of the 1,142 specimens were found to contain bacterial reads and bacteriophage-related contigs were detected in 27 specimens. This analysis demonstrates the capacity of metagenomic approaches to elucidate important vector-host-pathogen interactions of epidemiological significance.

Genomics reveals heterogeneous Plasmodium falciparum transmission and population differentiation in Zambia and bordering countries.

Genomic surveillance plays a critical role in monitoring malaria transmission and understanding how the parasite adapts in response to interventions. We conducted genomic surveillance of malaria by sequencing 241 Plasmodium falciparum genomes from regions with varying levels of malaria transmission across Zambia. We found genomic evidence of high levels of within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in the country. We identified country-level clustering of parasites from Zambia and neighboring countries, and distinct clustering of parasites from West Africa. Within Zambia, our identity by descent (IBD) relatedness analysis uncovered spatial clustering of closely related parasite pairs at the local level and rare cases of long-distance sharing. Genomic regions with large shared IBD segments and strong positive selection signatures identified genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Together, our findings enhance our understanding of P. falciparum transmission nationwide in Zambia and highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Genomics reveals heterogeneous Plasmodium falciparum transmission and selection signals in Zambia.

BACKGROUND: Genomic surveillance is crucial for monitoring malaria transmission and understanding parasite adaptation to interventions. Zambia lacks prior nationwide efforts in malaria genomic surveillance among African countries. METHODS: We conducted genomic surveillance of Plasmodium falciparum parasites from the 2018 Malaria Indicator Survey in Zambia, a nationally representative household survey of children under five years of age. We whole-genome sequenced and analyzed 241 P. falciparum genomes from regions with varying levels of malaria transmission across Zambia and estimated genetic metrics that are informative about transmission intensity, genetic relatedness between parasites, and selection. RESULTS: We provide genomic evidence of widespread within-host polygenomic infections, regardless of epidemiological characteristics, underscoring the extensive and ongoing endemic malaria transmission in Zambia. Our analysis reveals country-level clustering of parasites from Zambia and neighboring regions, with distinct separation in West Africa. Within Zambia, identity by descent (IBD) relatedness analysis uncovers local spatial clustering and rare cases of long-distance sharing of closely related parasite pairs. Genomic regions with large shared IBD segments and strong positive selection signatures implicate genes involved in sulfadoxine-pyrimethamine and artemisinin combination therapies drug resistance, but no signature related to chloroquine resistance. Furthermore, differences in selection signatures, including drug resistance loci, are observed between eastern and western Zambian parasite populations, suggesting variable transmission intensity and ongoing drug pressure. CONCLUSIONS: Our findings enhance our understanding of nationwide P. falciparum transmission in Zambia, establishing a baseline for analyzing parasite genetic metrics as they vary over time and space. These insights highlight the urgency of strengthening malaria control programs and surveillance of antimalarial drug resistance.

Malaria is caused by a parasite that is spread to humans via mosquito bites. It is a leading cause of death in children under five years old in sub-Saharan Africa. Analysis of the malaria parasite's complete set of DNA (its genome) can help us to understand transmission of the disease and how this changes in response to different strategies to control the disease. We analyzed the genomes of malaria parasites from children across Zambia. Our study revealed that 77% of children harbored multiple parasite strains, which suggests that local transmission (transmission between people within the same local area) is high. Genetic evidence for long-distance transmission was rarer. Furthermore, our findings suggest parasites are evolving in response to antimalarial drugs. Our study enhances our understanding of malaria dynamics in Zambia and may help to inform strategies for improved surveillance and control.

Temporal genomic analysis of Plasmodium falciparum reveals increased prevalence of mutations associated with delayed clearance following treatment with artemisinin-lumefantrine in Choma District, Southern Province, Zambia.

The emergence of antimalarial drug resistance is an impediment to malaria control and elimination in Africa. Analysis of temporal trends in molecular markers of resistance is critical to inform policy makers and guide malaria treatment guidelines. In a low and seasonal transmission region of southern Zambia, we successfully genotyped 85.5% (389/455) of Plasmodium falciparum samples collected between 2013-2018 from 8 spatially clustered health centres using molecular inversion probes (MIPs) targeting key drug resistance genes. Aside from one sample carrying K13 R622I, none of the isolates carried other World Health Organization-validated or candidate artemisinin partial resistance (ART-R) mutations in K13. However, 13% (CI, 9.6-17.2) of isolates had the AP2MU S160N mutation, which has been associated with delayed clearance following artemisinin combination therapy in Africa. This mutation increased in prevalence between 2015-2018 and bears a genomic signature of selection. During this time period, there was an increase in the MDR1 NFD haplotype that is associated with reduced susceptibility to lumefantrine. Sulfadoxine-pyrimethamine polymorphisms were near fixation. While validated ART-R mutations are rare, a mutation associated with slow parasite clearance in Africa appears to be under selection in southern Zambia.