Detection of low-density Plasmodium falciparum infections using amplicon deep sequencing

Application of a new highly multiplexed amplicon sequencing tool to evaluate Plasmodium falciparum antimalarial resistance and relatedness in individual and pooled samples from Dschang, Cameroon

Background

Resistance to antimalarial drugs remains a major obstacle to malaria elimination. Multiplexed, targeted amplicon sequencing is being adopted for surveilling resistance and dissecting the genetics of complex malaria infections. Moreover, genotyping of parasites and detection of molecular markers drug resistance in resource-limited regions requires open-source protocols for processing samples, using accessible reagents, and rapid methods for processing numerous samples including pooled sequencing.

Methods

P lasmodium f alciparum Streamlined Multiplex Antimalarial Resistance and Relatedness Testing (Pf-SMARRT) is a PCR-based amplicon panel consisting of 15 amplicons targeting antimalarial resistance mutations and 9 amplicons targeting hypervariable regions. This assay uses oligonucleotide primers in two pools and a non-proprietary library and barcoding approach.

Results

We evaluated Pf-SMARRT using control mocked dried blood spots (DBS) at varying levels of parasitemia and a mixture of 3D7 and Dd2 strains at known frequencies, showing the ability to genotype at low parasite density and recall within-sample allele frequencies. We then piloted Pf-SMARRT to genotype 100 parasite isolates collected from uncomplicated malaria cases at three health facilities in Dschang, Western Cameroon. Antimalarial resistance genotyping showed high levels of sulfadoxine-pyrimethamine resistance mutations, including 31% prevalence of the DHPS A613S mutation. No K13 candidate or validated artemisinin partial resistance mutations were detected, but one low-level non-synonymous change was observed. Pf-SMARRT's hypervariable targets, used to assess complexity of infections and parasite diversity and relatedness, showed similar levels and patterns compared to molecular inversion probe (MIP) sequencing. While there was strong concordance of antimalarial resistance mutations between individual samples and pools, low-frequency variants in the pooled samples were often missed.

Conclusion

Overall, Pf-SMARRT is a robust tool for assessing parasite relatedness and antimalarial drug resistance markers from both individual and pooled samples. Control samples support that accurate genotyping as low as 1 parasite per microliter is routinely possible.

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.

Tracing the origins of Plasmodium vivax resurgence after malaria elimination on Aneityum Island in Vanuatu

BACKGROUND

Five years after successful malaria elimination, Aneityum Island in Vanuatu experienced an outbreak of Plasmodium vivax of unknown origin in 2002. Epidemiological investigations revealed several potential sources of P. vivax. We aimed to identify the genetic origin of P. vivax responsible for the resurgence.

METHODS

Five P. vivax microsatellite markers were genotyped using DNA extracted from archived blood samples. A total of 69 samples from four P. vivax populations was included: 29 from the outbreak in 2002, seven from Aneityum in 1999 and 2000, 18 from visitors to Aneityum in 2000, and 15 from nearby Tanna Island in 2002. A neighbour-joining phylogenetic tree was constructed to elucidate the relationships among P. vivax isolates. STRUCTURE and principal component analysis were used to assess patterns of genetic structure.

RESULTS

Here we show distinct genetic origins of P. vivax during the outbreak on Aneityum. While the origin of most P. vivax lineages found during the outbreak remains unidentified, limited genetic diversity among these lineages is consistent with a rapid expansion from a recent common ancestor. Contemporaneous P. vivax from neighboring Tanna and potential relapse of P. vivax acquired from other islands in 1999 and 2000 are also identified as minor contributors to the outbreak.

CONCLUSIONS

Multiple reintroductions of P. vivax after elimination highlight the high receptivity and vulnerability to malaria resurgence in island settings of Vanuatu, despite robust surveillance and high community compliance to control measures.

Persistent and multiclonal malaria parasite dynamics despite extended artemether-lumefantrine treatment in children

Standard diagnostics used in longitudinal antimalarial studies are unable to characterize the complexity of submicroscopic parasite dynamics, particularly in high transmission settings. We use molecular markers and amplicon sequencing to characterize post-treatment stage-specific malaria parasite dynamics during a 42 day randomized trial of 3- versus 5 day artemether-lumefantrine in 303 children with and without HIV (ClinicalTrials.gov number NCT03453840). The prevalence of parasite-derived 18S rRNA is >70% in children throughout follow-up, and the ring-stage marker SBP1 is detectable in over 15% of children on day 14 despite effective treatment. We find that the extended regimen significantly lowers the risk of recurrent ring-stage parasitemia compared to the standard 3 day regimen, and that higher day 7 lumefantrine concentrations decrease the probability of ring-stage parasites in the early post-treatment period. Longitudinal amplicon sequencing reveals remarkably dynamic patterns of multiclonal infections that include new and persistent clones in both the early post-treatment and later time periods. Our data indicate that post-treatment parasite dynamics are highly complex despite efficacious therapy, findings that will inform strategies to optimize regimens in the face of emerging partial artemisinin resistance in Africa.

Genotyping Plasmodium falciparum gametocytes using amplicon deep sequencing

BACKGROUND

Understanding the dynamics of gametocyte production in polyclonal Plasmodium falciparum infections requires a genotyping method that detects distinct gametocyte clones and estimates their relative frequencies. Here, a marker was identified and evaluated to genotype P. falciparum mature gametocytes using amplicon deep sequencing.

METHODS

A data set of polymorphic regions of the P. falciparum genome was mined to identify a gametocyte genotyping marker. To assess marker resolution, the number of unique haplotypes in the marker region was estimated from 95 Malawian P. falciparum whole genome sequences. Specificity of the marker for detection of mature gametocytes was evaluated using reverse transcription-polymerase chain reaction of RNA extracted from NF54 mature gametocytes and rings from a non-gametocyte-producing strain of P. falciparum. Amplicon deep sequencing was performed on experimental mixtures of mature gametocytes from two distinct parasite clones, as well as gametocyte-positive P. falciparum field isolates to evaluate the quantitative ability and determine the limit of detection of the genotyping approach.

RESULTS

A 400 bp region of the pfs230 gene was identified as a gametocyte genotyping marker. A larger number of unique haplotypes was observed at the pfs230 marker (34) compared to the sera-2 (18) and ama-1 (14) markers in field isolates from Malawi. RNA and DNA genotyping accurately estimated gametocyte and total parasite clone frequencies when evaluating agreement between expected and observed haplotype frequencies in gametocyte mixtures, with concordance correlation coefficients of 0.97 [95% CI: 0.92-0.99] and 0.92 [95% CI: 0.83-0.97], respectively. The detection limit of the genotyping method for male gametocytes was 0.41 pfmget transcripts/µl [95% CI: 0.28-0.72] and for female gametocytes was 1.98 ccp4 transcripts/µl [95% CI: 1.35-3.68].

CONCLUSIONS

A region of the pfs230 gene was identified as a marker to genotype P. falciparum gametocytes. Amplicon deep sequencing of this marker can be used to estimate the number and relative frequency of parasite clones among mature gametocytes within P. falciparum infections. This gametocyte genotyping marker will be an important tool for studies aimed at understanding dynamics of gametocyte production in polyclonal P. falciparum infections.

Flexible and cost-effective genomic surveillance of P. falciparum malaria with targeted nanopore sequencing

Genomic surveillance of Plasmodium falciparum malaria can provide policy-relevant information about antimalarial drug resistance, diagnostic test failure, and the evolution of vaccine targets. Yet the large and low complexity genome of P. falciparum complicates the development of genomic methods, while resource constraints in malaria endemic regions can limit their deployment. Here, we demonstrate an approach for targeted nanopore sequencing of P. falciparum from dried blood spots (DBS) that enables cost-effective genomic surveillance of malaria in low-resource settings. We release software that facilitates flexible design of amplicon sequencing panels and use this software to design two target panels for P. falciparum. The panels generate 3-4 kbp reads for eight and sixteen targets respectively, covering key drug-resistance associated genes, diagnostic test antigens, polymorphic markers and the vaccine target csp. We validate our approach on mock and field samples, demonstrating robust sequencing coverage, accurate variant calls within coding sequences, the ability to explore P. falciparum within-sample diversity and to detect deletions underlying rapid diagnostic test failure.

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.

Drug resistance and vaccine target surveillance of Plasmodium falciparum using nanopore sequencing in Ghana

Malaria results in over 600,000 deaths annually, with the highest burden of deaths in young children living in sub-Saharan Africa. Molecular surveillance can provide important information for malaria control policies, including detection of antimalarial drug resistance. However, genome sequencing capacity in malaria-endemic countries is limited. We designed and implemented an end-to-end workflow to detect Plasmodium falciparum antimalarial resistance markers and diversity in the vaccine target circumsporozoite protein (csp) using nanopore sequencing in Ghana. We analysed 196 clinical samples and showed that our method is rapid, robust, accurate and straightforward to implement. Importantly, our method could be applied to dried blood spot samples, which are readily collected in endemic settings. We report that P. falciparum parasites in Ghana are mostly susceptible to chloroquine, with persistent sulfadoxine-pyrimethamine resistance and no evidence of artemisinin resistance. Multiple single nucleotide polymorphisms were identified in csp, but their significance is uncertain. Our study demonstrates the feasibility of nanopore sequencing for malaria genomic surveillance in endemic countries.

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.

Population structure and genetic connectivity of Plasmodium falciparum in pre-elimination settings of Southern Africa

To accelerate malaria elimination in the Southern African region by 2030, it is essential to prevent cross-border malaria transmission. However, countries within the region are highly interconnected due to human migration that aids in the movement of the parasite across geographical borders. It is therefore important to better understand Plasmodium falciparum transmission dynamics in the region, and identify major parasite source and sink populations, as well as cross-border blocks of high parasite connectivity. We performed a meta-analysis using collated parasite allelic data generated by microsatellite genotyping of malaria parasites from Namibia, Eswatini, South Africa, and Mozambique (N = 5,314). The overall number of unique alleles was significantly higher (P ≤ 0.01) in Namibia (mean A = 17.3 ± 1.46) compared to South Africa (mean A = 12.2 ± 1.22) and Eswatini (mean A = 13.3 ± 1.27, P ≤ 0.05), whilst the level of heterozygosity was not significantly different between countries. The proportion of polyclonal infections was highest for Namibia (77%), and lowest for Mozambique (64%). A was significant population structure was detected between parasites from the four countries, and patterns of gene flow showed that Mozambique was the major source area and Eswatini the major sink area of parasites between the countries. This study showed strong signals of parasite population structure and genetic connectivity between malaria parasite populations across national borders. This calls for strengthening the harmonization of malaria control and elimination efforts between countries in the southern African region. This data also proves its potential utility as an additional surveillance tool for malaria surveillance on both a national and regional level for the identification of imported cases and/or outbreaks, as well as monitoring for the potential spread of anti-malarial drug resistance as countries work towards malaria elimination.

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.

Plasmodium falciparum gametocyte carriage in longitudinally monitored incident infections is associated with duration of infection and human host factors

Malaria transmission depends on the presence of Plasmodium gametocytes that are the only parasite life stage that can infect mosquitoes. Gametocyte production varies between infections and over the course of infections. Infection duration is highly important for gametocyte production but poorly quantified. Between 2017 and 2019 an all-age cohort of individuals from Tororo, eastern Uganda was followed by continuous passive and routine assessments. We longitudinally monitored 104 incident infections from 98 individuals who were sampled once every 28 days and on any day of symptoms. Among infections that lasted ≥ 3 months, gametocyte appearance was near-universal with 96% of infections having detectable gametocytes prior to clearance. However, most infections were of much shorter duration; 55.7% of asymptomatic infections were detected only once. When considering all asymptomatic infections, regardless of their duration, only 36.3% had detectable gametocytes on at least one time-point prior to parasite clearance. Infections in individuals with sickle-cell trait (HbAS) were more likely to have gametocytes detected (Hazard Rate (HR) = 2.68, 95% CI 1.12, 6.38; p = 0.0231) and had gametocytes detected at higher densities (Density Ratio (DR) = 9.19, 95% CI 2.79, 30.23; p = 0.0002) compared to infections in wildtype (HbAA) individuals. Our findings suggest that a large proportion of incident infections is too short in duration and of too low density to contribute to onward transmission.

WHO antimalarial trial guidelines: good science, bad news?

Estimating antimalarial drug efficacy requires differentiating treatment failures from new infections arising during the several-week follow-up period in drug trials. Genetic profiling of malaria infections can guide this decision but is notoriously difficult in practice. Previous World Health Organisation (WHO) guidelines were based on assumptions with an inherently high risk of underestimating failure rates. A recent update to WHO guidelines recognises a wider range of analyses to overcome these limitations. We discuss these new analyses and their underlying logic. Drug failure rate estimates in moderate to high transmissions areas will become more accurate but will likely rise twofold due to better detection of treatment failures, and the malaria community needs to anticipate and prepare for potentially large increases in estimated failure rates.

Plasmodium falciparum Genetic Diversity in Coincident Human and Mosquito Hosts

Population genetic diversity of Plasmodium falciparum antigenic loci is high despite large bottlenecks in population size during the parasite life cycle. The prevalence of genetically distinct haplotypes at these loci, while well characterized in humans, has not been thoroughly compared between human and mosquito hosts. We assessed parasite haplotype prevalence, diversity, and evenness using human and mosquito P. falciparum infections collected from the same households during a 14-month longitudinal cohort study using amplicon deep sequencing of two antigenic gene fragments (ama1 and csp). To a prior set of infected humans (n = 1,175/2,813; 86.2% sequencing success) and mosquito abdomens (n = 199/1,448; 95.5% sequencing success), we added sequences from infected mosquito heads (n = 134/1,448; 98.5% sequencing success). The overall and sample-level parasite populations were more diverse in mosquitoes than in humans. Additionally, haplotype prevalences were more even in the P. falciparum human population than in the mosquito population, consistent with balancing selection occurring at these loci in humans. In contrast, we observed that infections in humans were more likely to harbor a dominant haplotype than infections in mosquitoes, potentially due to removal of unfit strains by the human immune system. Finally, within a given mosquito, there was little overlap in genetic composition of abdomen and head infections, suggesting that infections may be cleared from the abdomen during a mosquito’s lifespan. Taken together, our observations provide evidence for the mosquito vector acting as a reservoir of sequence diversity in malaria parasite populations.

Concurrent molecular characterization of sand flies and Leishmania parasites by amplicon-based next-generation sequencing

BACKGROUND

Phlebotomine sand flies are vectors of Leishmania parasites, which are the causative agents of leishmaniasis. Herein, we developed an amplicon-based next-generation sequencing (Amp-NGS) to characterize sand flies and Leishmania parasites simultaneously targeting partial fragments of 18S rDNA and ITS1 genes, respectively.

METHODS

Our assay was optimized using reference sand fly (n = 8) and Leishmania spp. (n = 9) samples and validated using wild-caught sand flies from Palestine. The assay was highly specific, and all DNA references were successfully identified to the species level.

RESULTS

Among the wild-caught sand flies (n = 187), Phlebotomus spp. represented 95% of the collected samples (177/187), including Ph. sergenti (147/187, 79%), Ph. papatasi (19/187, 10.2%), Ph. perfiliewi (3/187, 1.6%), Ph. tobbi (2/187, 1.2%) and Ph. syriacus (6/187, 3.2%). Sergentomyia spp. represented only 5% (10/187) of the collected samples and included S. dentata (n = 6), S. fallax (n = 2), S. schwetzi (n = 1) and S. ghesquiere (n = 1). The study observed strong positive correlation between sand fly identification results of the Amp-NGS and morphological identification method (r = 0.84, df = 185, P < 0.001). Some discrepancies between the two methods in the identification of closely related species (i.e. Ph. perfiliewi, Ph. tobbi and Ph. syriacus) were observed. Leishmania DNA was detected and identified as L. tropica in 14 samples (14/187, 7.5%).

CONCLUSIONS

Our assay was sensitive to detect (limit of detection was 0.0016 ng/reaction) and identify Leishmania DNA in sand flies, thus representing a new tool for studying sand flies and their associated Leishmania parasites in endemic areas.

Genetic Diversity of Plasmodium falciparum and Distribution of Antimalarial Drug Resistance Mutations in Symptomatic and Asymptomatic Infections

Malaria control relies on passive case detection, and this strategy fails detecting asymptomatic infections. In addition, infections in endemic areas harbor multiple parasite genotypes that could affect case management and malaria epidemiology.

Parasites and blood-meal hosts of the tsetse fly in Tanzania: a metagenomics study

Background

Tsetse flies can transmit various Trypanosoma spp. that cause trypanosomiasis in humans, wild animals, and domestic animals. Amplicon deep sequencing of the 12S ribosomal RNA (rRNA) gene can be used to detect mammalian tsetse hosts, and the 18S rRNA gene can be used to detect all associated eukaryotic pathogens, including Trypanosoma spp.

Methods

Tsetse flies were collected from the Serengeti National Park (n = 48), Maswa Game Reserve (n = 42), and Tarangire National Park (n = 49) in Tanzania in 2012–13. Amplicon deep sequencing targeting mammal-specific 12S rRNA and 18S rRNA genes was performed to screen the blood-feeding sources of tsetse flies and eukaryotic parasites in tsetse flies, respectively.

Results

12S rRNA gene deep sequencing revealed that various mammals were blood-feeding sources of the tsetse flies, including humans, common warthogs, African buffalos, mice, giraffes, African elephants, waterbucks, and lions. Genes of humans were less frequently detected in Serengeti (P = 0.0024), whereas African buffaloes were detected more frequently as a blood-feeding source (P = 0.0010). 18S rRNA gene deep sequencing showed that six tsetse samples harbored the Trypanosoma gene, which was identified as Trypanosoma godfreyi and Trypanosoma simiae in subsequent ITS1 gene sequencing.

Conclusions

Through amplicon deep sequencing targeting the 12S rRNA and 18S rRNA genes, various mammalian animals were identified as blood-meal sources, and two Trypanosoma species were detected in tsetse flies collected from the Maswa Game Reserve, Serengeti National Park, and Tarangire National Park in Tanzania. This study illustrates the patterns of parasitism of tsetse fly, wild animals targeted by the fly, and Trypanosoma spp. carried by the fly in Tanzania. It may provide essential data for formulating better strategies to control African trypanosomes.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05344-1.

A targeted amplicon sequencing panel to simultaneously identify mosquito species and Plasmodium presence across the entire Anopheles genus

Anopheles is a diverse genus of mosquitoes comprising over 500 described species, including all known human malaria vectors. While a limited number of key vector species have been studied in detail, the goal of malaria elimination calls for surveillance of all potential vector species. Here, we develop a multilocus amplicon sequencing approach that targets 62 highly variable loci in the Anopheles genome and two conserved loci in the Plasmodium mitochondrion, simultaneously revealing both the mosquito species and whether that mosquito carries malaria parasites. We also develop a cheap, nondestructive, and high-throughput DNA extraction workflow that provides template DNA from single mosquitoes for the multiplex PCR, which means specimens producing unexpected results can be returned to for morphological examination. Over 1000 individual mosquitoes can be sequenced in a single MiSeq run, and we demonstrate the panel's power to assign species identity using sequencing data for 40 species from Africa, Southeast Asia, and South America. We also show that the approach can be used to resolve geographic population structure within An. gambiae and An. coluzzii populations, as the population structure determined based on these 62 loci from over 1000 mosquitoes closely mirrors that revealed through whole genome sequencing. The end-to-end approach is quick, inexpensive, robust, and accurate, which makes it a promising technique for very large-scale mosquito genetic surveillance and vector control.

Genetic surveillance for monitoring the impact of drug use on Plasmodium falciparum populations

The use of antimalarial drugs is an effective strategy in the fight against malaria. However, selection of drug resistant parasites is a constant threat to the continued use of this approach. Antimalarial drugs are used not only to treat infections but also as part of population-level strategies to reduce malaria transmission toward elimination. While there is strong evidence that the ongoing use of antimalarial drugs increases the risk of the emergence and spread of drug-resistant parasites, it is less clear how population-level use of drug-based interventions like seasonal malaria chemoprevention (SMC) or mass drug administration (MDA) may contribute to drug resistance or loss of drug efficacy. Critical to sustained use of drug-based strategies for reducing the burden of malaria is the surveillance of population-level signals related to transmission reduction and resistance selection. Here we focus on Plasmodium falciparum and discuss the genetic signatures of a parasite population that are correlated with changes in transmission and related to drug pressure and resistance as a result of drug use. We review the evidence for MDA and SMC contributing to malaria burden reduction and drug resistance selection and examine the use and impact of these interventions in Senegal. Throughout we consider best strategies for ongoing surveillance of both population and resistance signals in the context of different parasite population parameters. Finally, we propose a roadmap for ongoing surveillance during population-level drug-based interventions to reduce the global malaria burden.

Comparative Analysis of Plasmodium falciparum Genotyping via SNP Detection, Microsatellite Profiling, and Whole-Genome Sequencing

Research efforts to combat antimalarial drug resistance rely on quick, robust, and sensitive methods to genetically characterize Plasmodium falciparum parasites. We developed a single-nucleotide polymorphism (SNP)-based genotyping method that can assess 33 drug resistance-conferring SNPs in dhfr, dhps, pfmdr1, pfcrt, and k13 in nine PCRs, performed directly from P. falciparum cultures or infected blood. We also optimized multiplexed fragment analysis and gel electrophoresis-based microsatellite typing methods using a set of five markers that can distinguish 12 laboratory strains of diverse geographical and temporal origin. We demonstrate how these methods can be applied to screen for the multidrug-resistant KEL1/PLA1/PfPailin (KelPP) lineage that has been sweeping across the Greater Mekong Subregion, verify parasite in vitro SNP-editing, identify novel recombinant genetic cross progeny, or cluster strains to infer their geographical origins. Results were compared with Illumina-based whole-genome sequence analysis that provides the most detailed sequence information but is cost-prohibitive. These adaptable, simple, and inexpensive methods can be easily implemented into routine genotyping of P. falciparum parasites in both laboratory and field settings.

Unraveling the Complexity of Imported Malaria Infections by Amplicon Deep Sequencing

Imported malaria and recurrent infections are becoming an emerging issue in many malaria non-endemic countries. This study aimed to determine the molecular patterns of the imported malaria infections and recurrence. Blood samples were collected from patients with imported malaria infections during 2016-2018 in Guangxi Zhuang Autonomous Region, China. Next-generation amplicon deep-sequencing approaches were used to assess parasite genetic diversity, multiplexity of infection, relapse, recrudescence, and antimalarial drug resistance. A total of 44 imported malaria cases were examined during the study, of which 35 (79.5%) had recurrent malaria infections within 1 year. The majority (91.4%) had one recurrent malaria episode, whereas two patients had two recurrences and one patient had three recurrences. A total of 19 recurrence patterns (the species responsible for primary and successive clinical episodes) were found in patients returning from malaria epidemic countries. Four parasite species were detected with a higher than usual proportion (46.2%) of non-falciparum infections or mixed-species infections. An increasing trend of recurrence infections and reduced drug treatment efficacy were observed among the cases of imported malaria. The high recurrence rate and complex patterns of imported malaria from Africa to non-endemic countries have the potential to initiate local transmission, thereby undermining efforts to eliminate locally acquired malaria. Our findings highlight the power of amplicon deep-sequencing applications in molecular epidemiological studies of the imported malaria recurrences.

Identification of Multiple Blastocystis Subtypes in Domestic Animals From Colombia Using Amplicon-Based Next Generation Sequencing

Blastocystis is frequently reported in fecal samples from animals and humans worldwide, and a variety of subtypes (STs) have been observed in wild and domestic animals. In Colombia, few studies have focused on the transmission dynamics and epidemiological importance of Blastocystis in animals. In this study, we characterized the frequency and subtypes of Blastocystis in fecal samples of domestic animals including pigs, minipigs, cows, dogs, horses, goats, sheep, and llama from three departments of Colombia. Of the 118 fecal samples included in this study 81.4% (n = 96) were positive for Blastocystis using a PCR that amplifies a fragment of the small subunit ribosomal RNA (SSU rRNA) gene. PCR positive samples were sequenced by next generation amplicon sequencing (NGS) to determine subtypes. Eleven subtypes were detected, ten previously reported, ST5 (50.7%), ST10 (47.8%), ST25 (34.3%), ST26 (29.8%), ST21 (22.4%), ST23 (22.4%), ST1 (17.9%), ST14 (16.4%), ST24 (14.9%), ST3 (7.5%), and a novel subtype, named ST32 (3.0%). Mixed infection and/or intra -subtype variations were identified in most of the samples. Novel ST32 was observed in two samples from a goat and a cow. To support novel subtype designation, a MinION based sequencing strategy was used to generate the full-length of the SSU rRNA gene. Comparison of full-length nucleotide sequences with those from current valid subtypes supported the designation of ST32. This is the first study in Colombia using NGS to molecularly characterize subtypes of Blastocystis in farm animals. A great diversity of subtypes was observed in domestic animals including subtypes previously identified in humans. Additionally, subtype overlap between the different hosts examined in this study were observed. These findings highlight the presence of Blastocystis subtypes with zoonotic potential in farm animals indicating that farm animals could play a role in transmission to humans.

Zoonotic Malaria: Non-Laverania Plasmodium Biology and Invasion Mechanisms

Malaria, which is caused by Plasmodium parasites through Anopheles mosquito transmission, remains one of the most life-threatening diseases affecting hundreds of millions of people worldwide every year. Plasmodium vivax, which accounts for the majority of cases of recurring malaria caused by the Plasmodium (non-Laverania) subgenus, is an ancient and continuing zoonosis originating from monkey hosts probably outside Africa. The emergence of other zoonotic malarias (P. knowlesi, P. cynomolgi, and P. simium) further highlights the seriousness of the disease. The severity of this epidemic disease is dependent on many factors, including the parasite characteristics, host-parasite interactions, and the pathology of the infection. Successful infection depends on the ability of the parasite to invade the host; however, little is known about the parasite invasion biology and mechanisms. The lack of this information adds to the challenges to malaria control and elimination, hence enhancing the potential for continuation of this zoonosis. Here, we review the literature describing the characteristics, distribution, and genome details of the parasites, as well as host specificity, host-parasite interactions, and parasite pathology. This information will provide the basis of a greater understanding of the epidemiology and pathogenesis of malaria to support future development of strategies for the control and prevention of this zoonotic infection.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Should deep-sequenced amplicons become the new gold-standard for analysing malaria drug clinical trials?

Regulatory clinical trials are required to ensure the continued supply and deployment of effective antimalarial drugs. Patient follow-up in such trials typically lasts several weeks, as the drugs have long half-lives and new infections often occur during this period. “Molecular correction” is therefore used to distinguish drug failures from new infections. The current WHO-recommended method for molecular correction uses length-polymorphic alleles at highly diverse loci but is inherently poor at detecting low-density clones in polyclonal infections. This likely leads to substantial underestimates of failure rates, delaying the replacement of failing drugs with potentially lethal consequences. Deep-sequenced amplicons (AmpSeq) substantially increase the detectability of low-density clones and may offer a new “gold standard” for molecular correction. Pharmacological simulation of clinical trials was used to evaluate the suitability of AmpSeq for molecular correction. We investigated the impact of factors such as the number of amplicon loci analyzed, the informatics criteria used to distinguish genotyping “noise” from real low-density signals, the local epidemiology of malaria transmission, and the potential impact of genetic signals from gametocytes. AmpSeq greatly improved molecular correction and provided accurate drug failure rate estimates. The use of 3 to 5 amplicons was sufficient, and simple, nonstatistical criteria could be used to classify recurrent infections as drug failures or new infections. These results suggest AmpSeq is strongly placed to become the new standard for molecular correction in regulatory trials, with potential extension into routine surveillance once the requisite technical support becomes established.

Navigating Clinical Utilization of Direct-from-Specimen Metagenomic Pathogen Detection: Clinical Applications, Limitations, and Testing Recommendations

BACKGROUND

Metagenomic next generation sequencing (mNGS) is becoming increasingly available for pathogen detection directly from clinical specimens. These tests use target-independent, shotgun sequencing to detect potentially unlimited organisms. The promise of this methodology to aid infection diagnosis is demonstrated through early case reports and clinical studies. However, the optimal role of mNGS in clinical microbiology remains uncertain.

CONTENT

We reviewed studies reporting clinical use of mNGS for pathogen detection from various specimen types, including cerebrospinal fluid, plasma, lower respiratory specimens, and others. Published clinical study data were critically evaluated and summarized to identify promising clinical indications for mNGS-based testing, to assess the clinical impact of mNGS for each indication, and to recognize test limitations. Based on these clinical studies, early testing recommendations are made to guide clinical utilization of mNGS for pathogen detection. Finally, current barriers to routine clinical laboratory implementation of mNGS tests are highlighted.

SUMMARY

The promise of direct-from-specimen mNGS to enable challenging infection diagnoses has been demonstrated through early clinical studies of patients with meningitis or encephalitis, invasive fungal infections, community acquired pneumonia, and other clinical indications. However, the proportion of patient cases with positive clinical impact due to mNGS testing is low in published studies and the cost of testing is high, emphasizing the importance of improving our understanding of 'when to test' and for which patients mNGS testing is appropriate.

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.

Antimalarial Drug Resistance and Implications for the WHO Global Technical Strategy

PURPOSE OF REVIEW

Five years have passed since the World Health Organization released its Global Technical Strategy for Malaria (GTS). In that time, progress against malaria has plateaued. This review focuses on the implications of antimalarial drug resistance for the GTS and how interim progress in parasite genomics and antimalarial pharmacology offer a bulwark against it.

RECENT FINDINGS

For the first time, drug resistance-conferring genes have been identified and validated before their global expansion in malaria parasite populations. More efficient methods for their detection and elaboration have been developed, although low-density infections and polyclonality remain a nuisance to be solved. Clinical trials of alternative regimens for multidrug-resistant malaria have delivered promising results. New agents continue down the development pipeline, while a nascent infrastructure in sub-Saharan Africa for conducting phase I trials and trials of transmission-blocking agents has come to fruition after years of preparation.

SUMMARY

These and other developments can help inform the GTS as the world looks ahead to the next two decades of its implementation. To remain ahead of the threat that drug resistance poses, wider application of genomic-based surveillance and optimization of existing and forthcoming antimalarial drugs are essential.

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.

Sensitive universal detection of blood parasites by selective pathogen-DNA enrichment and deep amplicon sequencing

BACKGROUND

Targeted amplicon deep sequencing (TADS) has enabled characterization of diverse bacterial communities, yet the application of TADS to communities of parasites has been relatively slow to advance. The greatest obstacle to this has been the genetic diversity of parasitic agents, which include helminths, protozoa, arthropods, and some acanthocephalans. Meanwhile, universal amplification of conserved loci from all parasites without amplifying host DNA has proven challenging. Pan-eukaryotic PCRs preferentially amplify the more abundant host DNA, obscuring parasite-derived reads following TADS. Flaherty et al. (2018) described a pan-parasitic TADS method involving amplification of eukaryotic 18S rDNA regions possessing restriction sites only in vertebrates. Using this method, host DNA in total DNA extracts could be selectively digested prior to PCR using restriction enzymes, thereby increasing the number of parasite-derived reads obtained following NGS. This approach showed promise though was only as sensitive as conventional PCR.

RESULTS

Here, we expand on this work by designing a second set of pan-eukaryotic primers flanking the priming sites already described, enabling nested PCR amplification of the established 18S rDNA target. This nested approach facilitated introduction of a second restriction digestion between the first and second PCR, reducing the proportional mass of amplifiable host-derived DNA while increasing the number of PCR amplification cycles. We applied this method to blood specimens containing Babesia, Plasmodium, various kinetoplastids, and filarial nematodes and confirmed its limit of detection (LOD) to be approximately 10-fold lower than previously described, falling within the range of most qPCR methods.

CONCLUSIONS

The assay detects and differentiates the major malaria parasites of humans, along with several other clinically important blood parasites. This represents an important step towards a TADS-based universal parasite diagnostic (UPDx) test with a sufficient LOD for routine applications. Video Abstract.

Sensitive, highly multiplexed sequencing of microhaplotypes from the Plasmodium falciparum heterozygome

Background

Targeted next-generation sequencing offers the potential for consistent, deep coverage of information-rich genomic regions to characterize polyclonal Plasmodium falciparum infections. However, methods to identify and sequence these genomic regions are currently limited.

Methods

A bioinformatic pipeline and multiplex methods were developed to identify and simultaneously sequence 100 targets and applied to dried blood spot (DBS) controls and field isolates from Mozambique. For comparison, whole-genome sequencing data were generated for the same controls.

Results

Using publicly available genomes, 4465 high-diversity genomic regions suited for targeted sequencing were identified, representing the P. falciparum heterozygome. For this study, 93 microhaplotypes with high diversity (median expected heterozygosity = 0.7) were selected along with 7 drug resistance loci. The sequencing method achieved very high coverage (median 99%), specificity (99.8%), and sensitivity (90% for haplotypes with 5% within sample frequency in dried blood spots with 100 parasites/µL). In silico analyses revealed that microhaplotypes provided much higher resolution to discriminate related from unrelated polyclonal infections than biallelic single-nucleotide polymorphism barcodes.

Conclusions

The bioinformatic and laboratory methods outlined here provide a flexible tool for efficient, low-cost, high-throughput interrogation of the P. falciparum genome, and can be tailored to simultaneously address multiple questions of interest in various epidemiological settings.

Implementing parasite genotyping into national surveillance frameworks: feedback from control programmes and researchers in the Asia-Pacific region

The Asia-Pacific region faces formidable challenges in achieving malaria elimination by the proposed target in 2030. Molecular surveillance of Plasmodium parasites can provide important information on malaria transmission and adaptation, which can inform national malaria control programmes (NMCPs) in decision-making processes. In November 2019 a parasite genotyping workshop was held in Jakarta, Indonesia, to review molecular approaches for parasite surveillance and explore ways in which these tools can be integrated into public health systems and inform policy. The meeting was attended by 70 participants from 8 malaria-endemic countries and partners of the Asia Pacific Malaria Elimination Network. The participants acknowledged the utility of multiple use cases for parasite genotyping including: quantifying the prevalence of drug resistant parasites, predicting risks of treatment failure, identifying major routes and reservoirs of infection, monitoring imported malaria and its contribution to local transmission, characterizing the origins and dynamics of malaria outbreaks, and estimating the frequency of Plasmodium vivax relapses. However, the priority of each use case varies with different endemic settings. Although a one-size-fits-all approach to molecular surveillance is unlikely to be applicable across the Asia-Pacific region, consensus on the spectrum of added-value activities will help support data sharing across national boundaries. Knowledge exchange is needed to establish local expertise in different laboratory-based methodologies and bioinformatics processes. Collaborative research involving local and international teams will help maximize the impact of analytical outputs on the operational needs of NMCPs. Research is also needed to explore the cost-effectiveness of genetic epidemiology for different use cases to help to leverage funding for wide-scale implementation. Engagement between NMCPs and local researchers will be critical throughout this process.

Amplicon deep sequencing improves Plasmodium falciparum genotyping in clinical trials of antimalarial drugs

Clinical trials monitoring malaria drug resistance require genotyping of recurrent Plasmodium falciparum parasites to distinguish between treatment failure and new infection occurring during the trial follow up period. Because trial participants usually harbour multi-clonal P. falciparum infections, deep amplicon sequencing (AmpSeq) was employed to improve sensitivity and reliability of minority clone detection. Paired samples from 32 drug trial participants were Illumina deep-sequenced for five molecular markers. Reads were analysed by custom-made software HaplotypR and trial outcomes compared to results from the previous standard genotyping method based on length-polymorphic markers. Diversity of AmpSeq markers in pre-treatment samples was comparable or higher than length-polymorphic markers. AmpSeq was highly reproducible with consistent quantification of co-infecting parasite clones within a host. Outcomes of the three best-performing markers, cpmp, cpp and ama1-D3, agreed in 26/32 (81%) of patients. Discordance between the three markers performed per sample was much lower by AmpSeq (six patients) compared to length-polymorphic markers (eleven patients). Using AmpSeq for discrimination of recrudescence and new infection in antimalarial drug trials provides highly reproducible and robust characterization of clone dynamics during trial follow-up. AmpSeq overcomes limitations inherent to length-polymorphic markers. Regulatory clinical trials of antimalarial drugs will greatly benefit from this unbiased typing method.