Longitudinal tracking and quantification of individual Plasmodium falciparum clones in complex infections

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.

Molecular tools are crucial for malaria elimination

The eradication of Plasmodium parasites, responsible for malaria, is a daunting global public health task. It requires a comprehensive approach that addresses symptomatic, asymptomatic, and submicroscopic cases. Overcoming this challenge relies on harnessing the power of molecular diagnostic tools, as traditional methods like microscopy and rapid diagnostic tests fall short in detecting low parasitaemia, contributing to the persistence of malaria transmission. By precisely identifying patients of all types and effectively characterizing malaria parasites, molecular tools may emerge as indispensable allies in the pursuit of malaria elimination. Furthermore, molecular tools can also provide valuable insights into parasite diversity, drug resistance patterns, and transmission dynamics, aiding in the implementation of targeted interventions and surveillance strategies. In this review, we explore the significance of molecular tools in the pursuit of malaria elimination, shedding light on their key contributions and potential impact on public health.

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.

Microsatellites reveal high polymorphism and high potential for use in anti-malarial efficacy studies in areas with different transmission intensities in mainland Tanzania

BACKGROUND

Tanzania is currently implementing therapeutic efficacy studies (TES) in areas of varying malaria transmission intensities as per the World Health Organization (WHO) recommendations. In TES, distinguishing reinfection from recrudescence is critical for the determination of anti-malarial efficacy. Recently, the WHO recommended genotyping polymorphic coding genes, merozoite surface proteins 1 and 2 (msp1 and msp2), and replacing the glutamate-rich protein (glurp) gene with one of the highly polymorphic microsatellites in Plasmodium falciparum to adjust the efficacy of antimalarials in TES. This study assessed the polymorphisms of six neutral microsatellite markers and their potential use in TES, which is routinely performed in Tanzania.

METHODS

Plasmodium falciparum samples were obtained from four TES sentinel sites, Kibaha (Pwani), Mkuzi (Tanga), Mlimba (Morogoro) and Ujiji (Kigoma), between April and September 2016. Parasite genomic DNA was extracted from dried blood spots on filter papers using commercial kits. Genotyping was done using six microsatellites (Poly-α, PfPK2, TA1, C3M69, C2M34 and M2490) by capillary method, and the data were analysed to determine the extent of their polymorphisms and genetic diversity at the four sites.

RESULTS

Overall, 83 (88.3%) of the 94 samples were successfully genotyped (with positive results for ≥ 50.0% of the markers), and > 50.0% of the samples (range = 47.6-59.1%) were polyclonal, with a mean multiplicity of infection (MOI) ranging from 1.68 to 1.88 among the four sites. There was high genetic diversity but limited variability among the four sites based on mean allelic richness (RS = 7.48, range = 7.27-8.03, for an adjusted minimum sample size of 18 per site) and mean expected heterozygosity (He = 0.83, range = 0.80-0.85). Cluster analysis of haplotypes using STRUCTURE, principal component analysis, and pairwise genetic differentiation (FST) did not reveal population structure or clustering of parasites according to geographic origin. Of the six markers, Poly-α was the most polymorphic, followed by C2M34, TA1 and C3M69, while M2490 was the least polymorphic.

CONCLUSION

Microsatellite genotyping revealed high polyclonality and genetic diversity but no significant population structure. Poly-α, C2M34, TA1 and C3M69 were the most polymorphic markers, and Poly-α alone or with any of the other three markers could be adopted for use in TES in Tanzania.

Plasmodium falciparum transmission in the highlands of Ethiopia is driven by closely related and clonal parasites

Malaria cases are frequently recorded in the Ethiopian highlands even at altitudes above 2000 m. The epidemiology of malaria in the Ethiopian highlands, and, in particular, the role of importation by human migration from the highly endemic lowlands is not well understood. We sequenced 187 Plasmodium falciparum samples from two sites in the Ethiopian highlands, Gondar (n = 159) and Ziway (n = 28), using a multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug resistance loci. Here, we characterize the parasite population structure and genetic relatedness. We identify moderate parasite diversity (mean HE : 0.54) and low infection complexity (74.9% monoclonal). A significant percentage of infections share microhaplotypes, even across transmission seasons and sites, indicating persistent local transmission. We identify multiple clusters of clonal or near-clonal infections, highlighting high genetic relatedness. Only 6.3% of individuals diagnosed with P. falciparum reported recent travel. Yet, in clonal or near-clonal clusters, infections of travellers were frequently observed first in time, suggesting that parasites may have been imported and then transmitted locally. 31.1% of infections are pfhrp2-deleted and 84.4% pfhrp3-deleted, and 28.7% have pfhrp2/3 double deletions. Parasites with pfhrp2/3 deletions and wild-type parasites are genetically distinct. Mutations associated with resistance to sulphadoxine-pyrimethamine or suggested to reduce sensitivity to lumefantrine are observed at near-fixation. In conclusion, genomic data corroborate local transmission and the importance of intensified control in the Ethiopian highlands.

Amplicon deep sequencing of five highly polymorphic markers of Plasmodium falciparum reveals high parasite genetic diversity and moderate population structure in Ethiopia

BACKGROUND

Plasmodium falciparum genetic diversity can add information on transmission intensity and can be used to track control and elimination interventions.

METHODS

Dried blood spots (DBS) were collected from patients who were recruited for a P. falciparum malaria therapeutic efficacy trial in three malaria endemic sites in Ethiopia from October to December 2015, and November to December 2019. qPCR-confirmed infections were subject to amplicon sequencing of polymorphic markers ama1-D3, csp, cpp, cpmp, msp7. Genetic diversity, the proportion of multiclonal infections, multiplicity of infection, and population structure were analysed.

RESULTS

Among 198 samples selected for sequencing, data was obtained for 181 samples. Mean MOI was 1.38 (95% CI 1.24-1.53) and 17% (31/181) of infections were polyclonal. Mean He across all markers was 0.730. Population structure was moderate; populations from Metema and Metehara 2015 were very similar to each other, but distinct from Wondogent 2015 and Metehara 2019.

CONCLUSION

The high level of parasite genetic diversity and moderate population structure in this study suggests frequent gene flow of parasites among sites. The results obtained can be used as a baseline for additional parasite genetic diversity and structure studies, aiding in the formulation of appropriate control strategies in Ethiopia.

A genomic platform for surveillance and antigen discovery in Plasmodium spp. using long-read amplicon sequencing

Many vaccine candidate proteins in the malaria parasite Plasmodium falciparum are under strong immunological pressure and confer antigenic diversity. We present a sequencing and data analysis platform for the genomic surveillance of the insertion or deletion (indel)-rich antigens merozoite surface protein 1 (MSP1), MSP2, glutamate-rich protein (GLURP), and CSP from P. falciparum using long-read circular consensus sequencing (CCS) in multiclonal malaria isolates. Our platform uses 40 PCR primers per gene to asymmetrically barcode and identify multiclonal infections in pools of up to 384 samples. With msp2, we validated the method using 235 mock infections combining 10 synthetic variants at different concentrations and infection complexities. We applied this strategy to P. falciparum isolates from a longitudinal cohort in Tanzania. Finally, we constructed an analysis pipeline that streamlines the processing and interpretation of epidemiological and antigenic diversity data from demultiplexed FASTQ files. This platform can be easily adapted to other polymorphic antigens of interest in Plasmodium or any other human pathogen.

Reconstructing multi-strain pathogen interactions from cross-sectional survey data via statistical network inference

Infectious diseases often involve multiple pathogen species or multiple strains of the same pathogen. As such, knowledge of how different pathogens interact is key to understand and predict the outcome of interventions targeting only a subset of species or strains involved in disease. Population-level data may be useful to infer pathogen strain interactions, but most previously used inference methods only consider uniform interactions between all strains or focus on marginal pairwise interactions. As such, these methods are prone to bias induced by indirect interactions through other strains. Here, we evaluated statistical network inference for reconstructing heterogeneous interactions from cross-sectional surveys detecting joint presence/absence patterns of pathogen strains within hosts. We applied various network models to simulated survey data, representing endemic infection states of multiple pathogen strains with potential interactions in acquisition or clearance of infection. Satisfactory performance was demonstrated by the estimators converging to the true interactions. Accurate reconstruction of interaction networks was achieved by regularization or penalization for sample size. Although performance deteriorated in the presence of host heterogeneity, this was overcome by correcting for individual-level risk factors. Our work demonstrates how statistical network inference could prove useful for detecting multi-strain pathogen interactions and may have applications beyond epidemiology.

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.

Irrigation-Induced Environmental Changes Sustain Malaria Transmission and Compromise Intervention Effectiveness

BACKGROUND

Irrigated agriculture enhances food security, but it potentially promotes mosquito-borne disease transmission and affects vector intervention effectiveness. This study was conducted in the irrigated and nonirrigated areas of rural Homa Bay and Kisumu Counties, Kenya.

METHODS

We performed cross-sectional and longitudinal surveys to determine Plasmodium infection prevalence, clinical malaria incidence, molecular force of infection (molFOI), and multiplicity of infection. We examined the impact of irrigation on the effectiveness of the new interventions.

RESULTS

We found that irrigation was associated with >2-fold higher Plasmodium infection prevalence and 3-fold higher clinical malaria incidence compared to the nonirrigated area. Residents in the irrigated area experienced persistent, low-density parasite infections and higher molFOI. Addition of indoor residual spraying was effective in reducing malaria burden, but the reduction was more pronounced in the nonirrigated area than in the irrigated area.

CONCLUSIONS

Our findings collectively suggest that irrigation may sustain and enhance Plasmodium transmission and affects intervention effectiveness.

Genomic approaches for monitoring transmission dynamics of malaria: A case for malaria molecular surveillance in Sub-Saharan Africa

Transmission dynamics is an important indicator for malaria control and elimination. As we move closer to eliminating malaria in Sub-Saharan Africa (sSA), transmission indices with higher resolution (genomic approaches) will complement our current measurements of transmission. Most of the present programmatic knowledge of malaria transmission patterns are derived from assessments of epidemiologic and clinical data, such as case counts, parasitological estimates of parasite prevalence, and Entomological Inoculation Rates (EIR). However, to eliminate malaria from endemic areas, we need to track changes in the parasite population and how they will impact transmission. This is made possible through the evolving field of genomics and genetics, as well as the development of tools for more in-depth studies on the diversity of parasites and the complexity of infections, among other topics. If malaria elimination is to be achieved globally, country-specific elimination activities should be supported by parasite genomic data from regularly collected blood samples for diagnosis, surveillance and possibly from other programmatic interventions. This presents a unique opportunity to track the spread of malaria parasites and shed additional light on intervention efficacy. In this review, various genetic techniques are highlighted along with their significance for an enhanced understanding of transmission patterns in distinct topological settings throughout Sub-Saharan Africa. The importance of these methods and their limitations in malaria surveillance to guide control and elimination strategies, are explored.

Asia-Pacific ICEMR: Understanding Malaria Transmission to Accelerate Malaria Elimination in the Asia Pacific Region

Gaining an in-depth understanding of malaria transmission requires integrated, multifaceted research approaches. The Asia-Pacific International Center of Excellence in Malaria Research (ICEMR) is applying specifically developed molecular and immunological assays, in-depth entomological assessments, and advanced statistical and mathematical modeling approaches to a rich series of longitudinal cohort and cross-sectional studies in Papua New Guinea and Cambodia. This is revealing both the essential contribution of forest-based transmission and the particular challenges posed by Plasmodium vivax to malaria elimination in Cambodia. In Papua New Guinea, these studies document the complex host-vector-parasite interactions that are underlying both the stunning reductions in malaria burden from 2006 to 2014 and the significant resurgence in transmission in 2016 to 2018. Here we describe the novel analytical, surveillance, molecular, and immunological tools that are being applied in our ongoing Asia-Pacific ICEMR research program.

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.

Amplicon sequencing as a potential surveillance tool for complexity of infection and drug resistance markers in Plasmodium falciparum asymptomatic infections

Background

Genotyping Plasmodium falciparum subpopulations in malaria infections is an important aspect of malaria molecular epidemiology to understand within-host diversity and the frequency of drug resistance markers.

Methods

We characterized P. falciparum genetic diversity in asymptomatic infections and subsequent first febrile infections using amplicon sequencing (AmpSeq) of ama1 in Coastal Kenya. We also examined temporal changes in haplotype frequencies of mdr1, a drug-resistant marker.

Results

We found >60% of the infections were polyclonal (complexity of infection [COI] >1) and there was a reduction in COI over time. Asymptomatic infections had a significantly higher mean COI than febrile infections based on ama1 sequences (2.7 [95% confidence interval {CI}, 2.65–2.77] vs 2.22 [95% CI, 2.17–2.29], respectively). Moreover, an analysis of 30 paired asymptomatic and first febrile infections revealed that many first febrile infections (91%) were due to the presence of new ama1 haplotypes. The mdr1-YY haplotype, associated with chloroquine and amodiaquine resistance, decreased over time, while the NY (wild type) and the NF (modulates response to lumefantrine) haplotypes increased.

Conclusions

This study emphasizes the utility of AmpSeq in characterizing parasite diversity as it can determine relative proportions of clones and detect minority clones. The usefulness of AmpSeq in antimalarial drug resistance surveillance is also highlighted.

Comparative analysis of targeted next-generation sequencing for Plasmodium falciparum drug resistance markers

Well-defined molecular resistance markers are available for a range of antimalarial drugs, and molecular surveillance is increasingly important for monitoring antimalarial drug resistance. Different genotyping platforms are available, but these have not been compared in detail. We compared Targeted Amplicon Deep sequencing (TADs) using Ion Torrent PGM with Illumina MiSeq for the typing of antimalarial drug resistance genes. We developed and validated protocols to type the molecular resistance markers pfcrt, pfdhfr, pfdhps, pfmdr1, pfkelch, and pfcytochrome b, in Plasmodium falciparum for the Ion Torrent PGM and Illumina MiSeq sequencing platforms. With P. falciparum 3D7 and K1 as reference strains, whole blood samples (N = 20) and blood spots from Rapid Diagnostic Test (RDT) samples (N = 5) from patients with uncomplicated falciparum malaria from Ubon Ratchathani were assessed on both platforms and compared for coverage (average reads per amplicon), sequencing accuracy, variant accuracy, false positive rate, false negative rate, and alternative allele detection, with conventional Sanger sequencing as the reference method for SNP calling. Both whole blood and RDT samples could be successfully sequenced using the Ion Torrent PGM and Illumina MiSeq platforms. Coverage of reads per amplicon was higher with Illumina MiSeq (28,886 reads) than with Ion Torrent PGM (1754 reads). In laboratory generated artificial mixed infections, the two platforms could detect the minor allele down to 1% density at 500X coverage. SNPs calls from both platforms were in complete agreement with conventional Sanger sequencing. The methods can be multiplexed with up to 96 samples per run, which reduces cost by 86% compared to conventional Sanger sequencing. Both platforms, using the developed TAD protocols, provide an accurate method for molecular surveillance of drug resistance markers in P. falciparum, but Illumina MiSeq provides higher coverage than Ion Torrent PGM.

SARS-CoV-2 variant detection with ADSSpike

The SARS-CoV-2 coronavirus pandemic has been an unprecedented challenge to global pandemic response and preparedness. With the continuous appearance of new SARS-CoV-2 variants, it is imperative to implement tools for genomic surveillance and diagnosis in order to decrease viral transmission and prevalence. The ADSSpike workflow was developed with the goal of identifying signature SNPs from the S gene associated with SARS-CoV-2 variants through amplicon deep sequencing. Seventy-two samples were sequenced, and 30 mutations were identified. Among those, signature SNPs were linked to 2 Zeta-VOI (P.2) samples and one to the Alpha-VOC (B.1.17). An average depth of 700 reads was found to properlycorrectly identify all SNPs and deletions pertinent to SARS-CoV-2 mutants. ADSSpike is the first workflow to provide a practical, cost-effective, and scalable solution to diagnose SARS-CoV-2 VOC/VOI in the clinical laboratory, adding a valuable tool to public health measures to fight the COVID-19 pandemic for approximately $41.85 USD/reaction.

Approximate likelihood-based estimation method of multiple-type pathogen interactions: An application to longitudinal pneumococcal carriage data

While the serotypes of Streptococcus pneumoniae are known to compete during colonization in human hosts, our knowledge of how competition occurs is still incomplete. New insights of pneumococcal between‐type competition could be generated from carriage data obtained by molecular‐based detection methods, which record more complete sets of serotypes involved in co‐carriage than when detection is done by culture. Here, we develop a Bayesian estimation method for inferring between‐type interactions from longitudinal data recording the presence/absence of the types at discrete observation times. It allows inference from data containing co‐carriage of two or more serotypes, which is often the case when pneumococcal presence is determined by molecular‐based methods. The computational burden posed by the increased number of types detected in co‐carriage is addressed by approximating the likelihood under a multi‐state model with the likelihood of only those trajectories with minimum number of acquisition and clearance events between observation times. The proposed method's performance was validated on simulated data. The estimates of the interaction parameters of acquisition and clearance were unbiased in settings with short sampling intervals between observation times. With less frequent sampling, the estimates of the interaction parameters became more biased, but their ratio, which summarizes the total interaction, remained unbiased. Confounding due to unobserved heterogeneity in exposure could be corrected by including individual‐level random effects. In an application to empirical data about pneumococcal carriage in infants, we found new evidence for between‐serotype competition in clearance, although the effect size was small.

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.

Systematic review of Plasmodium falciparum and Plasmodium vivax polyclonal infections: Impact of prevalence, study population characteristics, and laboratory procedures

Multiple infections of genetically distinct clones of the same Plasmodium species are common in many malaria endemic settings. Mean multiplicity of infection (MOI) and the proportion of polyclonal infections are often reported as surrogate marker of transmission intensity, yet the relationship with traditional measures such as parasite prevalence is not well understood. We have searched Pubmed for articles on P. falciparum and P. vivax multiplicity, and compared the proportion of polyclonal infections and mean MOI to population prevalence. The impact of the genotyping method, number of genotyping markers, method for diagnosis (microscopy/RDT vs. PCR), presence of clinical symptoms, age, geographic region, and year of sample collection on multiplicity indices were assessed. For P. falciparum, 153 studies met inclusion criteria, yielding 275 individual data points and 33,526 genotyped individuals. The proportion of polyclonal infections ranged from 0–96%, and mean MOI from 1–6.1. For P. vivax, 54 studies met inclusion criteria, yielding 115 data points and 13,325 genotyped individuals. The proportion of polyclonal infections ranged from 0–100%, and mean MOI from 1–3.8. For both species, the proportion of polyclonal infections ranged from very low to close to 100% at low prevalence, while at high prevalence it was always high. Each percentage point increase in prevalence resulted in a 0.34% increase in the proportion of polyclonal P. falciparum infections (P<0.001), and a 0.78% increase in the proportion of polyclonal P. vivax infections (P<0.001). In multivariable analysis, higher prevalence, typing multiple markers, diagnosis of infections by PCR, and sampling in Africa were found to result in a higher proportion of P. falciparum polyclonal infections. For P. vivax, prevalence, year of study, typing multiple markers, and geographic region were significant predictors. In conclusion, polyclonal infections are frequently present in all settings, but the association between multiplicity and prevalence is weak.

Development of a new barcode-based, multiplex-PCR, next-generation-sequencing assay and data processing and analytical pipeline for multiplicity of infection detection of Plasmodium falciparum

BACKGROUND

Simultaneous infection with multiple malaria parasite strains is common in high transmission areas. Quantifying the number of strains per host, or the multiplicity of infection (MOI), provides additional parasite indices for assessing transmission levels but it is challenging to measure accurately with current tools. This paper presents new laboratory and analytical methods for estimating the MOI of Plasmodium falciparum.

METHODS

Based on 24 single nucleotide polymorphisms (SNPs) previously identified as stable, unlinked targets across 12 of the 14 chromosomes within P. falciparum genome, three multiplex PCRs of short target regions and subsequent next generation sequencing (NGS) of the amplicons were developed. A bioinformatics pipeline including B4Screening pathway removed spurious amplicons to ensure consistent frequency calls at each SNP location, compiled amplicons by SNP site diversity, and performed algorithmic haplotype and strain reconstruction. The pipeline was validated by 108 samples generated from cultured-laboratory strain mixtures in different proportions and concentrations, with and without pre-amplification, and using whole blood and dried blood spots (DBS). The pipeline was applied to 273 smear-positive samples from surveys conducted in western Kenya, then providing results into StrainRecon Thresholding for Infection Multiplicity (STIM), a novel MOI estimator.

RESULTS

The 24 barcode SNPs were successfully identified uniformly across the 12 chromosomes of P. falciparum in a sample using the pipeline. Pre-amplification and parasite concentration, while non-linearly associated with SNP read depth, did not influence the SNP frequency calls. Based on consistent SNP frequency calls at targeted locations, the algorithmic strain reconstruction for each laboratory-mixed sample had 98.5% accuracy in dominant strains. STIM detected up to 5 strains in field samples from western Kenya and showed declining MOI over time (q < 0.02), from 4.32 strains per infected person in 1996 to 4.01, 3.56 and 3.35 in 2001, 2007 and 2012, and a reduction in the proportion of samples with 5 strains from 57% in 1996 to 18% in 2012.

CONCLUSION

The combined approach of new multiplex PCRs and NGS, the unique bioinformatics pipeline and STIM could identify 24 barcode SNPs of P. falciparum correctly and consistently. The methodology could be applied to field samples to reliably measure temporal changes in MOI.

Sex-based differences in clearance of chronic Plasmodium falciparum infection

Multiple studies have reported a male bias in incidence and/or prevalence of malaria infection in males compared to females. To test the hypothesis that sex-based differences in host-parasite interactions affect the epidemiology of malaria, we intensively followed Plasmodium falciparum infections in a cohort in a malaria endemic area of eastern Uganda and estimated both force of infection (FOI) and rate of clearance using amplicon deep-sequencing. We found no evidence of differences in behavioral risk factors, incidence of malaria, or FOI by sex. In contrast, females cleared asymptomatic infections at a faster rate than males (hazard ratio [HR]=1.82, 95% CI 1.20 to 2.75 by clone and HR = 2.07, 95% CI 1.24 to 3.47 by infection event) in multivariate models adjusted for age, timing of infection onset, and parasite density. These findings implicate biological sex-based differences as an important factor in the host response to this globally important pathogen.

Microgeographic Epidemiology of Malaria Parasites in an Irrigated Area of Western Kenya by Deep Amplicon Sequencing

To improve food security, investments in irrigated agriculture are anticipated to increase throughout Africa. However, the extent to which environmental changes from water resource development will impact malaria epidemiology remains unclear. This study was designed to compare the sensitivity of molecular markers used in deep amplicon sequencing for evaluating malaria transmission intensities and to assess malaria transmission intensity at various proximities to an irrigation scheme. Compared to ama1, csp, and msp1 amplicons, cpmp required the smallest sample size to detect differences in infection complexity between transmission risk zones. Transmission intensity was highest within 5 km of the irrigation scheme by polymerase chain reaction positivity rate, infection complexity, and linkage disequilibrium. The irrigated area provided a source of parasite infections for the surrounding 2- to 10-km area. This study highlights the suitability of the cpmp amplicon as a measure for transmission intensities and the impact of irrigation on microgeographic epidemiology of malaria parasites.

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.

Strengthening Surveillance Systems for Malaria Elimination by Integrating Molecular and Genomic Data

Unprecedented efforts in malaria control over the last 15 years have led to a substantial decrease in both morbidity and mortality in most endemic settings. However, these progresses have stalled over recent years, and resurgence may cause dramatic impact on both morbidity and mortality. Nevertheless, elimination efforts are currently going on with the objective of reducing malaria morbidity and mortality by 90% and malaria elimination in at least 35 countries by 2030. Strengthening surveillance systems is of paramount importance to reach those targets, and the integration of molecular and genomic techniques into routine surveillance could substantially improve the quality and robustness of data. Techniques such as polymerase chain reaction (PCR) and quantitative PCR (qPCR) are increasingly available in malaria endemic countries, whereas others such as sequencing are already available in a few laboratories. However, sequencing, especially next-generation sequencing (NGS), requires sophisticated infrastructure with adequate computing power and highly trained personnel for data analysis that require substantial investment. Different techniques will be required for different applications, and cost-effective planning must ensure the appropriate use of available resources. The development of national and sub-regional reference laboratories could help in minimizing the resources required in terms of equipment and trained staff. Concerted efforts from different stakeholders at national, sub-regional, and global level are needed to develop the required framework to establish and maintain these reference laboratories.

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.

Asymptomatic malaria in the clinical and public health context

Introduction: Historically, the global community has focused on the control of symptomatic malaria. However, interest in asymptomatic malaria has been growing, particularly in the context of malaria elimination.Areas covered: We undertook a comprehensive PubMed literature review on asymptomatic malaria as it relates to detection and elimination with emphasis between 2014 and 2019. Diagnostic tools with a low limit of detection (LOD) have allowed us to develop a more detailed understanding of asymptomatic malaria and its impact. These highly sensitive diagnostics have demonstrated that the prevalence of asymptomatic malaria is greater than previously thought. In addition, it is now possible to detect the malaria reservoir in the community, something that was previously not feasible. Asymptomatic malaria has previously not been treated, but research has begun to examine whether treating individuals with asymptomatic malaria may lead to health benefits. Finally, we have begun to understand the importance of asymptomatic malaria in ongoing transmission.Expert opinion: Therefore, with malaria elimination back on the agenda, asymptomatic malaria can no longer be ignored, especially in light of new ultra-sensitive diagnostic tools.

Understanding P. falciparum Asymptomatic Infections: A Proposition for a Transcriptomic Approach

Malaria is still a significant public health burden in the tropics. Infection with malaria causing parasites results in a wide range of clinical disease presentations, from severe to uncomplicated or mild, and in the poorly understood asymptomatic infections. The complexity of asymptomatic infections is due to the intricate interplay between factors derived from the human host, parasite, and environment. Asymptomatic infections often go undetected and provide a silent natural reservoir that sustains malaria transmission. This creates a major obstacle for malaria control and elimination efforts. Numerous studies have tried to characterize asymptomatic infections, unanimously revealing that host immunity is the underlying factor in the maintenance of these infections and in the risk of developing febrile malaria infections. An in-depth understanding of how host immunity and parasite factors interact to cause malaria disease tolerance is thus required. This review primarily focuses on understanding anti-inflammatory and pro-inflammatory responses to asymptomatic infections in malaria endemic areas, to present the view that it is potentially the shift in host immunity toward an anti-inflammatory profile that maintains asymptomatic infections after multiple exposures to malaria. Conversely, symptomatic infections are skewed toward a pro-inflammatory immune profile. Moreover, we propose that these infections can be better interrogated using next generation sequencing technologies, in particular RNA sequencing (RNA-seq), to investigate the immune system using the transcriptome sampled during a clearly defined asymptomatic infection.

Impact of Multiplicity of Plasmodium falciparum Infection on Clinical Disease in Malawi

Multiplicity of infection (MOI), the number of unique Plasmodium falciparum parasite genotypes found in one infected individual, may contribute to the development of clinical malaria disease. However, the independent contribution of MOI and parasite density to clinical disease has not been well characterized. We conducted a two-year longitudinal cohort study of adults and children in a high-transmission setting in Malawi to test the hypothesis that increased MOI was independently associated with clinical disease, after accounting for parasite density. Of 1,062 episodes of infection, 477 (44.9%) were associated with symptoms. After controlling for repeated measures within an individual, key demographic factors, and parasite density, there was no association between MOI and clinical disease (OR = 1.02, 95% CI: 0.70-1.51). Although the limited ability to discern MOI in low-density asymptomatic infections may have impacted our results, we conclude that MOI is not an independent risk factor for clinical disease.

Use cases for genetic epidemiology in malaria elimination

Background

While traditional epidemiological approaches have supported significant reductions in malaria incidence across many countries, higher resolution information about local and regional malaria epidemiology will be needed to efficiently target interventions for elimination. The application of genetic epidemiological methods for the analysis of parasite genetics has, thus far, primarily been confined to research settings. To illustrate how these technical methods can be used to advance programmatic and operational needs of National Malaria Control Programmes (NMCPs), and accelerate global progress to eradication, this manuscript presents seven use cases for which genetic epidemiology approaches to parasite genetic data are informative to the decision-making of NMCPs.

Methods

The use cases were developed through a highly iterative process that included an extensive review of the literature and global guidance documents, including the 2017 World Health Organization’s Framework for Malaria Elimination, and collection of stakeholder input. Semi-structured interviews were conducted with programmatic and technical experts about the needs and opportunities for genetic epidemiology methods in malaria elimination.

Results

Seven use cases were developed: Detect resistance, Assess drug resistance gene flow, Assess transmission intensity, Identify foci, Determine connectivity of parasite populations, Identify imported cases, and Characterize local transmission chains. The method currently used to provide the information sought, population unit for implementation, the pre-conditions for using these approaches, and post-conditions intended as a product of the use case were identified for each use case.

Discussion

This framework of use cases will prioritize research and development of genetic epidemiology methods that best achieve the goals of NMCPs, and ultimately, inform the establishment of normative policy guidance for their uses. With significant engagement of stakeholders from malaria endemic countries and collaboration with local programme experts to ensure strategic implementation, genetic epidemiological approaches have tremendous potential to accelerate global malaria elimination efforts.

Electronic supplementary material

The online version of this article (10.1186/s12936-019-2784-0) contains supplementary material, which is available to authorized users.