High-resolution micro-epidemiology of parasite spatial and temporal dynamics in a high malaria transmission setting in Kenya

Analytic optimization of Plasmodium falciparum marker gene haplotype recovery from amplicon deep sequencing of complex mixtures

Molecular epidemiologic studies of malaria parasites and other pathogens commonly employ amplicon deep sequencing (AmpSeq) of marker genes derived from dried blood spots (DBS) to answer public health questions related to topics such as transmission and drug resistance. As these methods are increasingly employed to inform direct public health action, it is important to rigorously evaluate the risk of false positive and false negative haplotypes derived from clinically-relevant sample types. We performed a control experiment evaluating haplotype recovery from AmpSeq of 5 marker genes (ama1, csp, msp7, sera2, and trap) from DBS containing mixtures of DNA from 1 to 10 known P. falciparum reference strains across 3 parasite densities in triplicate (n = 270 samples). While false positive haplotypes were present across all parasite densities and mixtures, we optimized censoring criteria to remove 83% (148/179) of false positives while removing only 8% (67/859) of true positives. Post-censoring, the median pairwise Jaccard distance between replicates was 0.83. We failed to recover 35% (477/1365) of haplotypes expected to be present in the sample. Haplotypes were more likely to be missed in low-density samples with <1.5 genomes/μL (OR: 3.88, CI: 1.82-8.27, vs. high-density samples with ≥75 genomes/μL) and in samples with lower read depth (OR per 10,000 reads: 0.61, CI: 0.54-0.69). Furthermore, minority haplotypes within a sample were more likely to be missed than dominant haplotypes (OR per 0.01 increase in proportion: 0.96, CI: 0.96-0.97). Finally, in clinical samples the percent concordance across markers for multiplicity of infection ranged from 40%-80%. Taken together, our observations indicate that, with sufficient read depth, the majority of haplotypes can be successfully recovered from DBS while limiting the false positive rate.

Genomic malaria surveillance of antenatal care users detects reduced transmission following elimination interventions in Mozambique

Routine sampling of pregnant women at first antenatal care (ANC) visits could make Plasmodium falciparum genomic surveillance more cost-efficient and convenient in sub-Saharan Africa. We compare the genetic structure of parasite populations sampled from 289 first ANC users and 93 children from the community in Mozambique between 2015 and 2019. Samples are amplicon sequenced targeting 165 microhaplotypes and 15 drug resistance genes. Metrics of genetic diversity and relatedness, as well as the prevalence of drug resistance markers, are consistent between the two populations. In an area targeted for elimination, intra-host genetic diversity declines in both populations (p = 0.002-0.007), while for the ANC population, population genetic diversity is also lower (p = 0.0004), and genetic relatedness between infections is higher (p = 0.002) than control areas, indicating a recent reduction in the parasite population size. These results highlight the added value of genomic surveillance at ANC clinics to inform about changes in transmission beyond epidemiological data.

Molecular markers for malaria genetic epidemiology: progress and pitfalls

Over recent years, progress in molecular markers for genotyping malaria parasites has enabled informative studies of epidemiology and transmission dynamics. Results have highlighted the value of these tools for surveillance to support malaria control and elimination strategies. There are many different types and panels of markers available for malaria parasite genotyping, and for end users, the nuances of these markers with respect to 'use case', resolution, and accuracy, are not well defined. This review clarifies issues surrounding different molecular markers and their application to malaria control and elimination. We describe available marker panels, use cases, implications for different transmission settings, limitations, access, cost, and data accuracy. The information provided can be used as a guide for molecular epidemiology and surveillance of malaria.

Genomic malaria surveillance of antenatal care users detects reduced transmission following elimination interventions in Mozambique

Routine sampling of pregnant women at first antenatal care (ANC) visits could make Plasmodium falciparum genomic surveillance more cost-efficient and convenient in sub-Saharan Africa. We compared the genetic structure of parasite populations sampled from 289 first ANC attendees and 93 children from the community in Mozambique between 2015 and 2019. Samples were amplicon sequenced targeting 165 microhaplotypes and 15 drug resistance genes. Metrics of genetic diversity and relatedness, as well as the prevalence of drug resistance markers, were consistent between the two populations. In an area targeted for elimination, intra-host genetic diversity declined in both populations (p=0.002-0.007), while for the ANC population, population genetic diversity was also lower (p=0.0004), and genetic relatedness between infections were higher (p=0.002) than control areas, indicating a recent reduction in the parasite population size. These results highlight the added value of genomic surveillance at ANC clinics to inform about changes in transmission beyond epidemiological data.

Sampling for malaria molecular surveillance

Strategic use of Plasmodium falciparum genetic variation has great potential to inform public health actions for malaria control and elimination. Malaria molecular surveillance (MMS) begins with a strategy to identify and collect parasite samples, guided by public-health priorities. In this review we discuss sampling design practices for MMS and point out epidemiological, biological, and statistical factors that need to be considered. We present examples for different use cases, including detecting emergence and spread of rare variants, establishing transmission sources and inferring changes in malaria transmission intensity. This review will potentially guide the collection of samples and data, serve as a starting point for further methodological innovation, and enhance utilization of MMS to support malaria elimination.

Analytic optimization of Plasmodium falciparum marker gene haplotype recovery from amplicon deep sequencing of complex mixtures

Molecular epidemiologic studies of malaria parasites commonly employ amplicon deep sequencing (AmpSeq) of marker genes derived from dried blood spots (DBS) to answer public health questions related to topics such as transmission and drug resistance. As these methods are increasingly employed to inform direct public health action, it is important to rigorously evaluate the risk of false positive and false negative haplotypes derived from clinically-relevant sample types. We performed a control experiment evaluating haplotype recovery from AmpSeq of 5 marker genes (ama1, csp, msp7, sera2, and trap) from DBS containing mixtures of DNA from 1 to 10 known P. falciparum reference strains across 3 parasite densities in triplicate (n=270 samples). While false positive haplotypes were present across all parasite densities and mixtures, we optimized censoring criteria to remove 83% (148/179) of false positives while removing only 8% (67/859) of true positives. Post-censoring, the median pairwise Jaccard distance between replicates was 0.83. We failed to recover 35% (477/1365) of haplotypes expected to be present in the sample. Haplotypes were more likely to be missed in low-density samples with <1.5 genomes/μL (OR: 3.88, CI: 1.82-8.27, vs. high-density samples with ≥75 genomes/μL) and in samples with lower read depth (OR per 10,000 reads: 0.61, CI: 0.54-0.69). Furthermore, minority haplotypes within a sample were more likely to be missed than dominant haplotypes (OR per 0.01 increase in proportion: 0.96, CI: 0.96-0.97). Finally, in clinical samples the percent concordance across markers for multiplicity of infection ranged from 40%-80%. Taken together, our observations indicate that, with sufficient read depth, haplotypes can be successfully recovered from DBS while limiting the false positive rate.

Amplicon Deep Sequencing Reveals Multiple Genetic Events Lead to Treatment Failure with Atovaquone-Proguanil in Plasmodium falciparum

Atovaquone-proguanil (AP) is used as treatment for uncomplicated malaria, and as a chemoprophylactic agent against Plasmodium falciparum. Imported malaria remains one of the top causes of fever in Canadian returning travelers. Twelve sequential whole-blood samples before and after AP treatment failure were obtained from a patient diagnosed with P. falciparum malaria upon their return from Uganda and Sudan. Ultradeep sequencing was performed on the cytb, dhfr, and dhps markers of treatment resistance before and during the episode of recrudescence. Haplotyping profiles were generated using three different approaches: msp2-3D7 agarose and capillary electrophoresis, and cpmp using amplicon deep sequencing (ADS). A complexity of infection (COI) analysis was conducted. De novo cytb Y268C mutants strains were observed during an episode of recrudescence 17 days and 16 h after the initial malaria diagnosis and AP treatment initiation. No Y268C mutant reads were observed in any of the samples prior to the recrudescence. SNPs in the dhfr and dhps genes were observed upon initial presentation. The haplotyping profiles suggest multiple clones mutating under AP selection pressure (COI > 3). Significant differences in COI were observed by capillary electrophoresis and ADS compared to the agarose gel results. ADS using cpmp revealed the lowest haplotype variation across the longitudinal analysis. Our findings highlight the value of ultra-deep sequencing methods in the understanding of P. falciparum haplotype infection dynamics. Longitudinal samples should be analyzed in genotyping studies to increase the analytical sensitivity.

How can the complex epidemiology of malaria in India impact its elimination?

Malaria is a human health hazard in the tropical and subtropical zones of the globe and is poised to be eliminated by the year 2030. Despite a decrease in incidence in the past two decades, many endemic countries, including India, report cases regularly. The epidemiology of malaria in India is unique owing to several features of the Plasmodium parasites, Anopheles vectors, ecoepidemiological situations conducive to disease transmission, and susceptible humans living in rural and forested areas. Limitations in public health reach, and poor health-seeking behaviour of vulnerable populations living in hard-to-reach areas, add to the problem. We bring all of these factors together in a comprehensive framework and opine that, in spite of complexities, targeted elimination of malaria in India is achievable with planned programmatic approaches.

Performance of SNP barcodes to determine genetic diversity and population structure of Plasmodium falciparum in Africa

Panels of informative biallelic single nucleotide polymorphisms (SNPs) have been proposed to be an economical method to fast-track the population genetic analysis of Plasmodium falciparum in malaria-endemic areas. Whilst used successfully in low-transmission areas where infections are monoclonal and highly related, we present the first study to evaluate the performance of these 24- and 96-SNP molecular barcodes in African countries, characterised by moderate-to-high transmission, where multiclonal infections are prevalent. For SNP barcodes it is generally recommended that the SNPs chosen i) are biallelic, ii) have a minor allele frequency greater than 0.10, and iii) are independently segregating, to minimise bias in the analysis of genetic diversity and population structure. Further, to be standardised and used in many population genetic studies, these barcodes should maintain characteristics i) to iii) across various iv) geographies and v) time points. Using haplotypes generated from the MalariaGEN P. falciparum Community Project version six database, we investigated the ability of these two barcodes to fulfil these criteria in moderate-to-high transmission African populations in 25 sites across 10 countries. Predominantly clinical infections were analysed, with 52.3% found to be multiclonal, generating high proportions of mixed-allele calls (MACs) per isolate thereby impeding haplotype construction. Of the 24- and 96-SNPs, loci were removed if they were not biallelic and had low minor allele frequencies in all study populations, resulting in 20- and 75-SNP barcodes respectively for downstream population genetics analysis. Both SNP barcodes had low expected heterozygosity estimates in these African settings and consequently biased analyses of similarity. Both minor and major allele frequencies were temporally unstable. These SNP barcodes were also shown to identify weak genetic differentiation across large geographic distances based on Mantel Test and DAPC. These results demonstrate that these SNP barcodes are vulnerable to ascertainment bias and as such cannot be used as a standardised approach for malaria surveillance in moderate-to-high transmission areas in Africa, where the greatest genomic diversity of P. falciparum exists at local, regional and country levels.

R H: a genetic metric for measuring intrahost Plasmodium falciparum relatedness and distinguishing cotransmission from superinfection

Multiple-strain (polygenomic) infections are a ubiquitous feature of Plasmodium falciparum parasite population genetics. Under simple assumptions of superinfection, polygenomic infections are hypothesized to be the result of multiple infectious bites. As a result, polygenomic infections have been used as evidence of repeat exposure and used to derive genetic metrics associated with high transmission intensity. However, not all polygenomic infections are the result of multiple infectious bites. Some result from the transmission of multiple, genetically related strains during a single infectious bite (cotransmission). Superinfection and cotransmission represent two distinct transmission processes, and distinguishing between the two could improve inferences regarding parasite transmission intensity. Here, we describe a new metric, R H, that utilizes the correlation in allelic state (heterozygosity) within polygenomic infections to estimate the likelihood that the observed complexity resulted from either superinfection or cotransmission. R H is flexible and can be applied to any type of genetic data. As a proof of concept, we used R H to quantify polygenomic relatedness and estimate cotransmission and superinfection rates from a set of 1,758 malaria infections genotyped with a 24 single nucleotide polymorphism (SNP) molecular barcode. Contrary to expectation, we found that cotransmission was responsible for a significant fraction of 43% to 53% of the polygenomic infections collected in three distinct epidemiological regions in Senegal. The prediction that polygenomic infections frequently result from cotransmission stresses the need to incorporate estimates of relatedness within polygenomic infections to ensure the accuracy of genomic epidemiology surveillance data for informing public health activities.

Hotspots and super-spreaders: Modelling fine-scale malaria parasite transmission using mosquito flight behaviour

Malaria hotspots have been the focus of public health managers for several years due to the potential elimination gains that can be obtained from targeting them. The identification of hotspots must be accompanied by the description of the overall network of stable and unstable hotspots of malaria, especially in medium and low transmission settings where malaria elimination is targeted. Targeting hotspots with malaria control interventions has, so far, not produced expected benefits. In this work we have employed a mechanistic-stochastic algorithm to identify clusters of super-spreader houses and their related stable hotspots by accounting for mosquito flight capabilities and the spatial configuration of malaria infections at the house level. Our results show that the number of super-spreading houses and hotspots is dependent on the spatial configuration of the villages. In addition, super-spreaders are also associated to house characteristics such as livestock and family composition. We found that most of the transmission is associated with winds between 6pm and 10pm although later hours are also important. Mixed mosquito flight (downwind and upwind both with random components) were the most likely movements causing the spread of malaria in two out of the three study areas. Finally, our algorithm (named MALSWOTS) provided an estimate of the speed of malaria infection progression from house to house which was around 200-400 meters per day, a figure coherent with mark-release-recapture studies of Anopheles dispersion. Cross validation using an out-of-sample procedure showed accurate identification of hotspots. Our findings provide a significant contribution towards the identification and development of optimal tools for efficient and effective spatio-temporal targeted malaria interventions over potential hotspot areas.

Plasmodium falciparum importation does not sustain malaria transmission in a semi-arid region of Kenya

Human movement impacts the spread and transmission of infectious diseases. Recently, a large reservoir of Plasmodium falciparum malaria was identified in a semi-arid region of northwestern Kenya historically considered unsuitable for malaria transmission. Understanding the sources and patterns of transmission attributable to human movement would aid in designing and targeting interventions to decrease the unexpectedly high malaria burden in the region. Toward this goal, polymorphic parasite genes (ama1, csp) in residents and passengers traveling to Central Turkana were genotyped by amplicon deep sequencing. Genotyping and epidemiological data were combined to assess parasite importation. The contribution of travel to malaria transmission was estimated by modelling case reproductive numbers inclusive and exclusive of travelers. P. falciparum was detected in 6.7% (127/1891) of inbound passengers, including new haplotypes which were later detected in locally-transmitted infections. Case reproductive numbers approximated 1 and did not change when travelers were removed from transmission networks, suggesting that transmission is not fueled by travel to the region but locally endemic. Thus, malaria is not only prevalent in Central Turkana but also sustained by local transmission. As such, interrupting importation is unlikely to be an effective malaria control strategy on its own, but targeting interventions locally has the potential to drive down transmission.

Exposure to diverse Plasmodium falciparum genotypes shapes the risk of symptomatic malaria in incident and persistent infections: A longitudinal molecular epidemiologic study in Kenya

BACKGROUND

Repeated exposure to malaria infections could protect against symptomatic progression, as people develop adaptive immunity to infections acquired over time.

METHODS

We investigated how new, recurrent, and persistent Plasmodium falciparum infections were associated with the odds of developing symptomatic compared to asymptomatic malaria. Using a 14-month longitudinal cohort in Western Kenya, we used amplicon deep sequencing of two polymorphic genes (pfama1 and pfcsp) to assess overlap of parasite genotypes (represented by haplotypes) acquired within an individual's successive infections. We hypothesized infections with novel haplotypes would increase the odds of symptomatic malaria.

RESULTS

After excluding initial infections, we observed 534 asymptomatic and 88 symptomatic infections across 186 people. We detected 109 pfcsp haplotypes, and each infection was classified as harboring novel, recurrent or persistent haplotypes. Incident infections with only new haplotypes had higher odds of symptomatic malaria when compared to infections with only recurrent haplotypes [odds ratio (OR): 3.24, 95% confidence interval (CI): 1.20 to 8.78], but infections with both new and recurrent haplotypes [OR: 0.64, 95% CI: 0.15 to 2.65] did not. Assessing persistent infections, those with mixed (persistent with new or recurrent) haplotypes [OR: 0.77, 95% CI: 0.21 to 2.75] had no association with symptomatic malaria compared to infections with only persistent haplotypes. Results were similar for pfama1.

CONCLUSIONS

These results confirm that incident infections with only novel haplotypes were associated with increased odds of symptomatic malaria compared to infections with only recurrent haplotypes but this relationship was not seen when haplotypes persisted over time in consecutive infections.

Genotyping cognate Plasmodium falciparum in humans and mosquitoes to estimate onward transmission of asymptomatic infections

Malaria control may be enhanced by targeting reservoirs of Plasmodium falciparum transmission. One putative reservoir is asymptomatic malaria infections and the scale of their contribution to transmission in natural settings is not known. We assess the contribution of asymptomatic malaria to onward transmission using a 14-month longitudinal cohort of 239 participants in a high transmission site in Western Kenya. We identify P. falciparum in asymptomatically- and symptomatically-infected participants and naturally-fed mosquitoes from their households, genotype all parasites using deep sequencing of the parasite genes pfama1 and pfcsp, and use haplotypes to infer participant-to-mosquito transmission through a probabilistic model. In 1,242 infections (1,039 in people and 203 in mosquitoes), we observe 229 (pfcsp) and 348 (pfama1) unique parasite haplotypes. Using these to link human and mosquito infections, compared with symptomatic infections, asymptomatic infections more than double the odds of transmission to a mosquito among people with both infection types (Odds Ratio: 2.56; 95% Confidence Interval (CI): 1.36-4.81) and among all participants (OR 2.66; 95% CI: 2.05-3.47). Overall, 94.6% (95% CI: 93.1-95.8%) of mosquito infections likely resulted from asymptomatic infections. In high transmission areas, asymptomatic infections are the major contributor to mosquito infections and may be targeted as a component of transmission reduction.

Within-household clustering of genetically related Plasmodium falciparum infections in a moderate transmission area of Uganda

BACKGROUND

Evaluation of genetic relatedness of malaria parasites is a useful tool for understanding transmission patterns, but patterns are not easily detectable in areas with moderate to high malaria transmission. To evaluate the feasibility of detecting genetic relatedness in a moderate malaria transmission setting, relatedness of Plasmodium falciparum infections was measured in cohort participants from randomly selected households in the Kihihi sub-county of Uganda (annual entomological inoculation rate of 27 infectious bites per person).

METHODS

All infections detected via microscopy or Plasmodium-specific loop mediated isothermal amplification from passive and active case detection during August 2011-March 2012 were genotyped at 26 microsatellite loci, providing data for 349 samples from 230 participants living in 80 households. Pairwise genetic relatedness was calculated using identity by state (IBS).

RESULTS

As expected, genetic diversity was high (mean heterozygosity [He] = 0.73), and the majority (76.5 %) of samples were polyclonal. Despite the high genetic diversity, fine-scale population structure was detectable, with significant spatiotemporal clustering of highly related infections. Although the difference in malaria incidence between households at higher (mean 1127 metres) versus lower elevation (mean 1015 metres) was modest (1.4 malaria cases per person-year vs. 1.9 per person-year, respectively), there was a significant difference in multiplicity of infection (2.2 vs. 2.6, p = 0.008) and, more strikingly, a higher proportion of highly related infections within households (6.3 % vs. 0.9 %, p = 0.0005) at higher elevation compared to lower elevation.

CONCLUSIONS

Genetic data from a relatively small number of diverse, multiallelic loci reflected fine scale patterns of malaria transmission. Given the increasing interest in applying genetic data to augment malaria surveillance, this study provides evidence that genetic data can be used to inform transmission patterns at local spatial scales even in moderate transmission areas.

Household prevalence of P. falciparum, P. vivax, and P. ovale in the Democratic Republic of the Congo, 2013-2014

In a cross-sectional molecular study in the Democratic Republic of Congo, 78% of households had at least one member infected with Plasmodium falciparum, vivax, and/or ovale spp. 47% of children and 33% of adults tested positive for at least one species. Risk factors varied by species and age group.