Malaria Epidemiology at the Clone Level

Extensive genetic diversity in Plasmodium vivax from Sudan and its genetic relationships with other geographical isolates

Plasmodium vivax, traditionally overlooked has experienced a notable increase in cases in East Africa. This study investigated the geographical origin and genetic diversity of P. vivax in Sudan using 14 microsatellite markers. A total of 113 clinical P. vivax samples were collected from two different ecogeographical zones, New Halfa and Khartoum, in Sudan. Additionally, 841 geographical samples from the database were incorporated for a global genetic analysis to discern genetic relationships among P. vivax isolates on regional and worldwide scales. On the regional scale, our findings revealed 91 unique and 8 shared haplotypes among the Sudan samples, showcasing a remarkable genetic diversity compared to other geographical isolates and supporting the hypothesis that P. vivax originated from Africa. On a global scale, distinct genetic clustering of P. vivax isolates from Africa, South America, and Asia (including Papua New Guinea and Solomon Island) was observed, with limited admixture among the three clusters. Principal component analysis emphasized the substantial contribution of African isolates to the observed global genetic variation. The Sudanese populations displayed extensive genetic diversity, marked by significant multi-locus linkage disequilibrium, suggesting an ancestral source of P. vivax variation globally and frequent recombination among the isolates. Notably, the East African P. vivax exhibited similarity with some Asian isolates, indicating potential recent introductions. Overall, our results underscore the effectiveness of utilizing microsatellite markers for implementing robust control measures, given their ability to capture extensive genetic diversity and linkage disequilibrium patterns.

Evaluating the performance of Plasmodium falciparum genetic metrics for inferring National Malaria Control Programme reported incidence in Senegal

BACKGROUND

Genetic surveillance of the Plasmodium falciparum parasite shows great promise for helping National Malaria Control Programmes (NMCPs) assess parasite transmission. Genetic metrics such as the frequency of polygenomic (multiple strain) infections, genetic clones, and the complexity of infection (COI, number of strains per infection) are correlated with transmission intensity. However, despite these correlations, it is unclear whether genetic metrics alone are sufficient to estimate clinical incidence.

METHODS

This study examined parasites from 3147 clinical infections sampled between the years 2012-2020 through passive case detection (PCD) across 16 clinic sites spread throughout Senegal. Samples were genotyped with a 24 single nucleotide polymorphism (SNP) molecular barcode that detects parasite strains, distinguishes polygenomic (multiple strain) from monogenomic (single strain) infections, and identifies clonal infections. To determine whether genetic signals can predict incidence, a series of Poisson generalized linear mixed-effects models were constructed to predict the incidence level at each clinical site from a set of genetic metrics designed to measure parasite clonality, superinfection, and co-transmission rates.

RESULTS

Model-predicted incidence was compared with the reported standard incidence data determined by the NMCP for each clinic and found that parasite genetic metrics generally correlated with reported incidence, with departures from expected values at very low annual incidence (< 10/1000/annual [‰]).

CONCLUSIONS

When transmission is greater than 10 cases per 1000 annual parasite incidence (annual incidence > 10‰), parasite genetics can be used to accurately infer incidence and is consistent with superinfection-based hypotheses of malaria transmission. When transmission was < 10‰, many of the correlations between parasite genetics and incidence were reversed, which may reflect the disproportionate impact of importation and focal transmission on parasite genetics when local transmission levels are low.

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.

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.

Can incorporating genotyping data into efficacy estimators improve efficiency of early phase malaria vaccine trials?

BACKGROUND

Early phase malaria vaccine field trials typically measure malaria infection by PCR or thick blood smear microscopy performed on serially sampled blood. Vaccine efficacy (VE) is the proportion reduction in an endpoint due to vaccination and is often calculated as VEHR = 1-hazard ratio or VERR = 1-risk ratio. Genotyping information can distinguish different clones and distinguish multiple infections over time, potentially increasing statistical power. This paper investigates two alternative VE endpoints incorporating genotyping information: VEmolFOI, the vaccine-induced proportion reduction in incidence of new clones acquired over time, and VEC, the vaccine-induced proportion reduction in mean number of infecting clones per exposure.

METHODS

Power of VEmolFOI and VEC was compared to that of VEHR and VERR by simulations and analytic derivations, and the four VE methods were applied to three data sets: a Phase 3 trial of RTS,S malaria vaccine in 6912 African infants, a Phase 2 trial of PfSPZ Vaccine in 80 Burkina Faso adults, and a trial comparing Plasmodium vivax incidence in 466 Papua New Guinean children after receiving chloroquine + artemether lumefantrine with or without primaquine (as these VE methods can also quantify effects of other prevention measures). By destroying hibernating liver-stage P. vivax, primaquine reduces subsequent reactivations after treatment completion.

RESULTS

In the trial of RTS,S vaccine, a significantly reduced number of clones at first infection was observed, but this was not the case in trials of PfSPZ Vaccine or primaquine, although the PfSPZ trial lacked power to show a reduction. Resampling smaller data sets from the large RTS,S trial to simulate phase 2 trials showed modest power gains from VEC compared to VEHR for data like those from RTS,S, but VEC is less powerful than VEHR for trials in which the number of clones at first infection is not reduced. VEmolFOI was most powerful in model-based simulations, but only the primaquine trial collected enough serial samples to precisely estimate VEmolFOI. The primaquine VEmolFOI estimate decreased after most control arm liver-stage infections reactivated (which mathematically resembles a waning vaccine), preventing VEmolFOI from improving power.

CONCLUSIONS

The power gain from the genotyping methods depends on the context. Because input parameters for early phase power calculations are often uncertain, these estimators are not recommended as primary endpoints for small trials unless supported by targeted data analysis.

TRIAL REGISTRATIONS

NCT00866619, NCT02663700, NCT02143934.

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.

Evaluating the performance of Plasmodium falciparum genetics for inferring National Malaria Control Program reported incidence in Senegal

Genetic surveillance of the Plasmodium falciparum parasite shows great promise for helping National Malaria Control Programs (NMCPs) assess parasite transmission. Genetic metrics such as the frequency of polygenomic (multiple strain) infections, genetic clones, and the complexity of infection (COI, number of strains per infection) are correlated with transmission intensity. However, despite these correlations, it is unclear whether genetic metrics alone are sufficient to estimate clinical incidence. Here, we examined parasites from 3,147 clinical infections sampled between the years 2012-2020 through passive case detection (PCD) across 16 clinic sites spread throughout Senegal. Samples were genotyped with a 24 single nucleotide polymorphism (SNP) molecular barcode that detects parasite strains, distinguishes polygenomic (multiple strain) from monogenomic (single strain) infections, and identifies clonal infections. To determine whether genetic signals can predict incidence, we constructed a series of Poisson generalized linear mixed-effects models to predict the incidence level at each clinical site from a set of genetic metrics designed to measure parasite clonality, superinfection, and co-transmission rates. We compared the model-predicted incidence with the reported standard incidence data determined by the NMCP for each clinic and found that parasite genetic metrics generally correlated with reported incidence, with departures from expected values at very low annual incidence (<10/1000/annual [‰]). When transmission is greater than 10 cases per 1000 annual parasite incidence (annual incidence >10 ‰), parasite genetics can be used to accurately infer incidence and is consistent with superinfection-based hypotheses of malaria transmission. When transmission was <10 ‰, we found that many of the correlations between parasite genetics and incidence were reversed, which we hypothesize reflects the disproportionate impact of importation and focal transmission on parasite genetics when local transmission levels are low.

Can incorporating genotyping data into efficacy estimators improve efficiency of early phase malaria vaccine trials?

Background

Early phase malaria vaccine field trials typically measure malaria infection by PCR or thick blood smear microscopy performed on serially sampled blood. Vaccine efficacy (VE) is the proportion reduction in an endpoint due to vaccination and is often calculated as V E H R = 1 - hazard ratio or V E R R = 1 - risk ratio. Genotyping information can distinguish different clones and distinguish multiple infections over time, potentially increasing statistical power. This paper investigates two alternative VE endpoints incorporating genotyping information: V E m o l F O I , the vaccine-induced proportion reduction in incidence of new clones acquired over time, and V E C , the vaccine-induced proportion reduction in mean number of infecting clones per exposure.

Methods

We used simulations and analytic derivations to compare power of these methods to V E H R and V E R R and applied them to three data sets: a Phase 3 trial of RTS,S malaria vaccine in 6912 African infants, a Phase 2 trial of PfSPZ Vaccine in 80 Burkina Faso adults, and a trial comparing Plasmodium vivax incidence in 466 Papua New Guinean children after receiving chloroquine + artemether lumefantrine with or without primaquine (as these VE methods can also quantify effects of other prevention measures). By destroying hibernating liver-stage P. vivax, primaquine reduces subsequent reactivations after treatment completion.

Results

The RTS,S vaccine significantly reduced the number of clones at first infection, but PfSPZ vaccine and primaquine did not. Resampling smaller data sets from the large RTS,S trial to simulate phase 2 trials showed modest power gains from V E C compared to V E H R for data like RTS,S, but V E C is less powerful than V E H R for vaccines which do not reduce the number of clones at first infection. V E m o l F O I was most powerful in model-based simulations, but only the primaquine trial collected enough serial samples to precisely estimate V E m o l F O I . The primaquine V E m o l F O I estimate decreased after most control arm liver-stage infections reactivated (which mathematically resembles a waning vaccine), preventing V E m o l F O I from improving power.

Conclusions

The power gain from the genotyping methods depends on the context. Because input parameters for early phase power calculations are often uncertain, we recommend against these estimators as primary endpoints for small trials unless supported by targeted data analysis.

Trial registrations

NCT00866619, NCT02663700, NCT02143934.

Plasmodium falciparum population structure inferred by msp1 amplicon sequencing of parasites collected from febrile patients in Kenya

BACKGROUND

Multiplicity of infection (MOI) is an important measure of Plasmodium falciparum diversity, usually derived from the highly polymorphic genes, such as msp1, msp2 and glurp as well as microsatellites. Conventional methods of deriving MOI lack fine resolution needed to discriminate minor clones. This study used amplicon sequencing (AmpliSeq) of P. falciparum msp1 (Pfmsp1) to measure spatial and temporal genetic diversity of P. falciparum.

METHODS

264 P. falciparum positive blood samples collected from areas of differing malaria endemicities between 2010 and 2019 were used. Pfmsp1 gene was amplified and amplicon libraries sequenced on Illumina MiSeq. Sequences were aligned against a reference sequence (NC_004330.2) and clustered to detect fragment length polymorphism and amino acid variations.

RESULTS

Children < 5 years had higher parasitaemia (median = 23.5 ± 5 SD, p = 0.03) than the > 5-14 (= 25.3 ± 5 SD), and those > 15 (= 25.1 ± 6 SD). Of the alleles detected, 553 (54.5%) were K1, 250 (24.7%) MAD20 and 211 (20.8%) RO33 that grouped into 19 K1 allelic families (108-270 bp), 14 MAD20 (108-216 bp) and one RO33 (153 bp). AmpliSeq revealed nucleotide polymorphisms in alleles that had similar sizes, thus increasing the K1 to 104, 58 for MAD20 and 14 for RO33. By AmpliSeq, the mean MOI was 4.8 (± 0.78, 95% CI) for the malaria endemic Lake Victoria region, 4.4 (± 1.03, 95% CI) for the epidemic prone Kisii Highland and 3.4 (± 0.62, 95% CI) for the seasonal malaria Semi-Arid region. MOI decreased with age: 4.5 (± 0.76, 95% CI) for children < 5 years, compared to 3.9 (± 0.70, 95% CI) for ages 5 to 14 and 2.7 (± 0.90, 95% CI) for those > 15. Females' MOI (4.2 ± 0.66, 95% CI) was not different from males 4.0 (± 0.61, 95% CI). In all regions, the number of alleles were high in the 2014-2015 period, more so in the Lake Victoria and the seasonal transmission arid regions.

CONCLUSION

These findings highlight the added advantages of AmpliSeq in haplotype discrimination and the associated improvement in unravelling complexity of P. falciparum population structure.

Multiplexed ddPCR-amplicon sequencing reveals isolated Plasmodium falciparum populations amenable to local elimination in Zanzibar, Tanzania

Zanzibar has made significant progress toward malaria elimination, but recent stagnation requires novel approaches. We developed a highly multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug-resistance loci, and successfully sequenced 290 samples from five districts covering both main islands. Here, we elucidate fine-scale Plasmodium falciparum population structure and infer relatedness and connectivity of infections using an identity-by-descent (IBD) approach. Despite high genetic diversity, we observe pronounced fine-scale spatial and temporal parasite genetic structure. Clusters of near-clonal infections on Pemba indicate persistent local transmission with limited parasite importation, presenting an opportunity for local elimination efforts. Furthermore, we observe an admixed parasite population on Unguja and detect a substantial fraction (2.9%) of significantly related infection pairs between Zanzibar and the mainland, suggesting recent importation. Our study provides a high-resolution view of parasite genetic structure across the Zanzibar archipelago and provides actionable insights for prioritizing malaria elimination efforts.

Biology and epidemiology of Plasmodium falciparum and Plasmodium vivax gametocyte carriage: Implication for malaria control and elimination

Malaria is among the leading public health problems worldwide. Female anopheles mosquito orchestrates the transmission of malaria by taking gametocytes and introducing sporozoite while taking blood meals. Interrupting transmission is the major strategy for malaria elimination. The gametocyte stage is essential for the onward transmission of malaria. Thus, understanding its basic biology and epidemiology is key to malaria control and elimination. Therefore, the current review focuses on revealing the biology, prevalence, and determinants of gametocyte carriage as well as its implication on mitigation of malaria. It also illustrates the role of asymptomatic and sub-microscopic Plasmodium infections and G-6-PD deficiency in gametocyte carriage and hence malaria transmission. Gametocytogenesis is initiated at committed merozoites and gives rise to the development of gametocytes. The trigger for gametocytogenesis depends on the host, parasite, and intervention factors. Gametocytes pass through five developmental stages identifiable by molecular markers. A considerable number of malaria patients carry gametocytes at a sub-microscopic level, thereby serving as a potential infectious reservoir of transmission. Factors involving the human host, Plasmodium parasite, and intervention parameters play a critical role in gametocyte biology and prevalence. The contribution of asymptomatic and sub-microscopic infections to malaria transmission is unknown. The clear impact of G-6-PD deficiency on malaria control and elimination remains unclear. Lack of clarity on such issues might impede the success of interventions. Basic science and epidemiological studies should continue to overcome the challenges and cope with the ever-evolving parasite and guide interventions.

Molecular Surveillance of Malaria Using the PF AmpliSeq Custom Assay for Plasmodium falciparum Parasites from Dried Blood Spot DNA Isolates from Peru

Malaria molecular surveillance has great potential to support national malaria control programs (NMCPs), informing policy for its control and elimination. Here, we present a new three-day workflow for targeted resequencing of markers in 13 resistance-associated genes, histidine rich protein 2 and 3 (hrp2&3), a country (Peru)-specific 28 SNP-barcode for population genetic analysis, and apical membrane antigen 1 (ama1), using Illumina short-read sequencing technology. The assay applies a multiplex PCR approach to amplify all genomic regions of interest in a rapid and easily standardizable procedure and allows simultaneous amplification of a high number of targets at once, therefore having great potential for implementation into routine surveillance practice by NMCPs. The assay can be performed on routinely collected filter paper blood spots and can be easily adapted to different regions to investigate either regional trends or in-country epidemiological changes.

Drought correlates with reduced infection complexity and possibly prevalence in a decades-long study of the lizard malaria parasite Plasmodium mexicanum

Microparasites often exist as a collection of genetic 'clones' within a single host (termed multi-clonal, or complex, infections). Malaria parasites are no exception, with complex infections playing key roles in parasite ecology. Even so, we know little about what factors govern the distribution and abundance of complex infections in natural settings. Utilizing a natural dataset that spans more than 20 years, we examined the effects of drought conditions on infection complexity and prevalence in the lizard malaria parasite Plasmodium mexicanum and its vertebrate host, the western fence lizard, Sceloporus occidentalis. We analyzed data for 14,011 lizards sampled from ten sites over 34 years with an average infection rate of 16.2%. Infection complexity was assessed for 546 infected lizards sampled during the most recent 20 years. Our data illustrate significant, negative effects of drought-like conditions on infection complexity, with infection complexity expected to increase by a factor of 2.27 from the lowest to highest rainfall years. The relationship between rainfall and parasite prevalence is somewhat more ambiguous; when prevalence is modeled over the full range in years, a 50% increase in prevalence is predicted between the lowest and highest rainfall years, but this trend is not apparent or is reversed when data are analyzed over a shorter timeframe. To our knowledge, this is the first reported evidence for drought affecting the abundance of multi-clonal infections in malaria parasites. It is not yet clear what mechanism might connect drought with infection complexity, but the correlation we observed suggests that additional research on how drought influences parasite features like infection complexity, transmission rates and within-host competition may be worthwhile.

Malaria Molecular Surveillance in the Peruvian Amazon with a Novel Highly Multiplexed Plasmodium falciparum AmpliSeq Assay

Molecular surveillance for malaria has great potential to support national malaria control programs (NMCPs). To bridge the gap between research and implementation, several applications (use cases) have been identified to align research, technology development, and public health efforts. For implementation at NMCPs, there is an urgent need for feasible and cost-effective tools. We designed a new highly multiplexed deep sequencing assay (Pf AmpliSeq), which is compatible with benchtop sequencers, that allows high-accuracy sequencing with higher coverage and lower cost than whole-genome sequencing (WGS), targeting genomic regions of interest. The novelty of the assay is its high number of targets multiplexed into one easy workflow, combining population genetic markers with 13 nearly full-length resistance genes, which is applicable for many different use cases. We provide the first proof of principle for hrp2 and hrp3 deletion detection using amplicon sequencing. Initial sequence data processing can be performed automatically, and subsequent variant analysis requires minimal bioinformatic skills using any tabulated data analysis program. The assay was validated using a retrospective sample collection (n = 254) from the Peruvian Amazon between 2003 and 2018. By combining phenotypic markers and a within-country 28-single-nucleotide-polymorphism (SNP) barcode, we were able to distinguish different lineages with multiple resistance haplotypes (in dhfr, dhps, crt and mdr1) and hrp2 and hrp3 deletions, which have been increasing in recent years. We found no evidence to suggest the emergence of artemisinin (ART) resistance in Peru. These findings indicate a parasite population that is under drug pressure but is susceptible to current antimalarials and demonstrate the added value of a highly multiplexed molecular tool to inform malaria strategies and surveillance systems. IMPORTANCE While the power of next-generation sequencing technologies to inform and guide malaria control programs has become broadly recognized, the integration of genomic data for operational incorporation into malaria surveillance remains a challenge in most countries where malaria is endemic. The main obstacles include limited infrastructure, limited access to high-throughput sequencing facilities, and the need for local capacity to run an in-country analysis of genomes at a large-enough scale to be informative for surveillance. In addition, there is a lack of standardized laboratory protocols and automated analysis pipelines to generate reproducible and timely results useful for relevant stakeholders. With our standardized laboratory and bioinformatic workflow, malaria genetic surveillance data can be readily generated by surveillance researchers and malaria control programs in countries of endemicity, increasing ownership and ensuring timely results for informed decision- and policy-making.

Recent Advances in Imported Malaria Pathogenesis, Diagnosis, and Management

Purpose of Review

Malaria is an important human parasitic disease affecting the population of tropical, subtropical regions as well as travelers to these areas.The purpose of this article is to provide clinicians practicing in non-endemic areas with a comprehensive overview of the recent data on microbiologic and pathophysiologic features of five Plasmodium parasites, clinical presentation of uncomplicated and severe cases, modern diagnostic methods, and treatment of malaria.

Recent Findings

Employment of robust surveillance programs, rapid diagnostic tests, highly active artemisinin-based therapy, and the first malaria vaccine have led to decline in malaria incidence; however, emerging drug resistance, disruptions due to the COVID-19 pandemic, and other socio-economic factors have stalled the progress.

Summary

Clinicians practicing in non-endemic areas such as the United States should consider a diagnosis of malaria in returning travelers presenting with fever, utilize rapid diagnostic tests if available at their practice locations in addition to microscopy, and timely initiate guideline-directed management as delays in treatment can lead to poor clinical outcomes.

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

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

Heterogeneity in prevalence of subclinical Plasmodium falciparum and Plasmodium vivax infections but no parasite genomic clustering in the Chittagong Hill Tracts, Bangladesh

BACKGROUND

Malaria remains endemic in Bangladesh, with the majority of cases occurring in forested, mountainous region in the Chittagong Hill Tracts (CHT). This area is home to Bengali and diverse groups of indigenous people (Pahari) residing largely in mono-ethnic villages.

METHODS

1002 individuals of the 9 most prominent Pahari and the Bengali population were randomly selected and screened by RDT and qPCR. Parasites were genotyped by msp2 and deep sequencing of 5 amplicons (ama1-D3, cpmp, cpp, csp, and msp7) for Plasmodium falciparum (n = 20), and by microsatellite (MS) typing of ten loci and amplicon sequencing of msp1 for Plasmodium vivax (n = 21). Population structure was analysed using STRUCTURE software. Identity-by-state (IBS) was calculated as a measure of parasite relatedness and used to generate relatedness networks.

RESULTS

The prevalence of P. falciparum and P. vivax infection was 0.7% by RDT (P. falciparum 6/1002; P. vivax 0/1002, mixed: 1/1002) and 4% by qPCR (P. falciparum 21/1002; P. vivax 16/1002, mixed: 5/1002). Infections were highly clustered, with 64% (27/42) of infections occurring in only two Pahari groups, the Khumi and Mro. Diversity was high; expected heterozygosity was 0.93 for P. falciparum and 0.81 for P. vivax. 85.7% (18/21) of P. vivax and 25% (5/20) of P. falciparum infections were polyclonal. No population structure was evident for either species, suggesting high transmission and gene flow among Pahari groups.

CONCLUSIONS

High subclinical infection prevalence and genetic diversity mirror ongoing transmission. Control activities should be specifically directed to Pahari groups at greatest risk.

Tenfold difference in DNA recovery rate: systematic comparison of whole blood vs. dried blood spot sample collection for malaria molecular surveillance

BACKGROUND

Molecular and genomic surveillance is becoming increasingly used to track malaria control and elimination efforts. Blood samples can be collected as whole blood and stored at - 20 °C until DNA extraction, or as dried blood spots (DBS), circumventing the need for a cold chain. Despite the wide use of either method, systematic comparisons of how the method of blood sample preservation affects the limit of detection (LOD) of molecular diagnosis and the proportion of DNA recovered for downstream applications are lacking.

METHODS

Extractions based on spin columns, magnetic beads, Tween-Chelex, and direct PCR without prior extraction were compared for whole blood and dried blood spots (DBS) using dilution series of Plasmodium falciparum culture samples. Extracted DNA was quantified by qPCR and droplet digital PCR (ddPCR).

RESULTS

DNA recovery was 5- to 10-fold higher for whole blood compared to DBS, resulting in a 2- to 3-fold lower LOD for both extraction methods compared to DBS. For whole blood, a magnetic bead-based method resulted in a DNA recovery rate of 88-98% when extracting from whole blood compared to 17-33% for a spin-column based method. For extractions from DBS, the magnetic bead-based method resulted in 8-20% DNA recovery, while the spin-column based method resulted in only 2% DNA recovery. The Tween-Chelex method was superior to other methods with 15-21% DNA recovery, and even more sensitive than extractions from whole blood samples. The direct PCR method was found to have the lowest LOD overall for both, whole blood and DBS.

CONCLUSIONS

Pronounced differences in LOD and DNA yield need to be considered when comparing prevalence estimates based on molecular methods and when selecting sampling protocols for other molecular surveillance applications.

Plasmodium falciparum genomic surveillance reveals spatial and temporal trends, association of genetic and physical distance, and household clustering

Molecular epidemiology using genomic data can help identify relationships between malaria parasite population structure, malaria transmission intensity, and ultimately help generate actionable data to assess the effectiveness of malaria control strategies. Genomic data, coupled with geographic information systems data, can further identify clusters or hotspots of malaria transmission, parasite genetic and spatial connectivity, and parasite movement by human or mosquito mobility over time and space. In this study, we performed longitudinal genomic surveillance in a cohort of 70 participants over four years from different neighborhoods and households in Thiès, Senegal—a region of exceptionally low malaria transmission (entomological inoculation rate less than 1). Genetic identity (identity by state, IBS) was established using a 24-single nucleotide polymorphism molecular barcode, identity by descent was calculated from whole genome sequence data, and a hierarchical Bayesian regression model was used to establish genetic and spatial relationships. Our results show clustering of genetically similar parasites within households and a decline in genetic similarity of parasites with increasing distance. One household showed extremely high diversity and warrants further investigation as to the source of these diverse genetic types. This study illustrates the utility of genomic data with traditional epidemiological approaches for surveillance and detection of trends and patterns in malaria transmission not only by neighborhood but also by household. This approach can be implemented regionally and countrywide to strengthen and support malaria control and elimination efforts.

Persistent Plasmodium falciparum infections enhance transmission-reducing immunity development

Subclinical infections that serve as reservoir populations to drive transmission remain a hurdle to malaria control. Data on infection dynamics in a geographical area is required to strategically design and implement malaria interventions. In a longitudinal cohort, we monitored Plasmodium falciparum infection prevalence and persistence, and anti-parasite immunity to gametocyte and asexual antigens for 10 weeks. Of the 100 participants, only 11 were never infected, whilst 16 had persistent infections detected by reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), and one participant had microscopic parasites at all visits. Over 70% of the participants were infected three or more times, and submicroscopic gametocyte prevalence was high, ≥ 48% of the parasite carriers. Naturally induced responses against recombinant Pfs48/45.6C, Pfs230proC, and EBA175RIII-V antigens were not associated with either infection status or gametocyte carriage, but the antigen-specific IgG titers inversely correlated with parasite and gametocyte densities consistent with partial immunity. Longitudinal analysis of gametocyte diversity indicated at least four distinct clones circulated throughout the study period. The high prevalence of children infected with distinct gametocyte clones coupled with marked variation in infection status at the individual level suggests ongoing transmission and should be targeted in malaria control programs.

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.

Diagnosing the drug resistance signature in Plasmodium falciparum: a review from contemporary methods to novel approaches

The genome sequence project of the human malaria parasite Plasmodium falciparum reveal variations in the parasite DNA sequence. Most of these variations are single nucleotide polymorphism (SNP). A high frequency of single nucleotide polymorphism (SNP) in the Plasmodium falciparum population is usually a benchmark for anti-malarial resistance which allows parasites to be elusive to the chemotherapeutic agents, vaccine and vector control strategies, resulting in the leading cause of morbidity and mortality globally. The high density of drug resistance signature markers such as pfcrt,pfmdr1, pfdhps, pfdhfr, pfkelch13, pfatpase6 and pfmrp1 in the genome opens up a scope for the study of the genetic basis of this elusive parasite. The precise and prompt diagnosis of resistance strains of parasite plays vital role in malaria elimination program.This review probably shed light on contemporary SNP diagnostic tools used in molecular surveillance of Plasmodium falciparum drug resistance in terms of mechanism, reaction modalities, and development with their virtues and shortcomings.

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.

Monitoring Plasmodium falciparum and Plasmodium vivax using microsatellite markers indicates limited changes in population structure after substantial transmission decline in Papua New Guinea

Monitoring the genetic structure of pathogen populations may be an economical and sensitive approach to quantify the impact of control on transmission dynamics, highlighting the need for a better understanding of changes in population genetic parameters as transmission declines. Here we describe the first population genetic analysis of two major human malaria parasites, Plasmodium falciparum (Pf) and Plasmodium vivax (Pv), following nationwide distribution of long-lasting insecticide-treated nets (LLINs) in Papua New Guinea (PNG). Parasite isolates from pre- (2005-2006) and post-LLIN (2010-2014) were genotyped using microsatellite markers. Despite parasite prevalence declining substantially (East Sepik Province: Pf = 54.9%-8.5%, Pv = 35.7%-5.6%, Madang Province: Pf = 38.0%-9.0%, Pv: 31.8%-19.7%), genetically diverse and intermixing parasite populations remained. Pf diversity declined modestly post-LLIN relative to pre-LLIN (East Sepik: R_s  = 7.1-6.4, H_E  = 0.77-0.71; Madang: R_s  = 8.2-6.1, H_E  = 0.79-0.71). Unexpectedly, population structure present in pre-LLIN populations was lost post-LLIN, suggesting that more frequent human movement between provinces may have contributed to higher gene flow. Pv prevalence initially declined but increased again in one province, yet diversity remained high throughout the study period (East Sepik: R_s  = 11.4-9.3, H_E  = 0.83-0.80; Madang: R_s  = 12.2-14.5, H_E  = 0.85-0.88). Although genetic differentiation values increased between provinces over time, no significant population structure was observed at any time point. For both species, a decline in multiple infections and increasing clonal transmission and significant multilocus linkage disequilibrium post-LLIN were positive indicators of impact on the parasite population using microsatellite markers. These parameters may be useful adjuncts to traditional epidemiological tools in the early stages of transmission reduction.

Molecular epidemiology of Plasmodium falciparum by multiplexed amplicon deep sequencing in Senegal

BACKGROUND

Molecular epidemiology can provide important information regarding the genetic diversity and transmission of Plasmodium falciparum, which can assist in designing and monitoring elimination efforts. However, malaria molecular epidemiology including understanding the genetic diversity of the parasite and performing molecular surveillance of transmission has been poorly documented in Senegal. Next Generation Sequencing (NGS) offers a practical, fast and high-throughput approach to understand malaria population genetics. This study aims to unravel the population structure of P. falciparum and to estimate the allelic diversity, multiplicity of infection (MOI), and evolutionary patterns of the malaria parasite using the NGS platform.

METHODS

Multiplex amplicon deep sequencing of merozoite surface protein 1 (PfMSP1) and merozoite surface protein 2 (PfMSP2) in fifty-three P. falciparum isolates from two epidemiologically different areas in the South and North of Senegal, was carried out.

RESULTS

A total of 76 Pfmsp1 and 116 Pfmsp2 clones were identified and 135 different alleles were found, 56 and 79 belonged to the pfmsp1 and pfmsp2 genes, respectively. K1 and IC3D7 allelic families were most predominant in both sites. The local haplotype diversity (Hd) and nucleotide diversity (π) were higher in the South than in the North for both genes. For pfmsp1, a high positive Tajima's D (TD) value was observed in the South (D = 2.0453) while negative TD value was recorded in the North (D = - 1.46045) and F-Statistic (Fst) was 0.19505. For pfmsp2, non-directional selection was found with a highly positive TD test in both areas and Fst was 0.02111. The mean MOI for both genes was 3.07 and 1.76 for the South and the North, respectively, with a statistically significant difference between areas (p = 0.001).

CONCLUSION

This study revealed a high genetic diversity of pfmsp1 and pfmsp2 genes and low genetic differentiation in P. falciparum population in Senegal. The MOI means were significantly different between the Southern and Northern areas. Findings also showed that multiplexed amplicon deep sequencing is a useful technique to investigate genetic diversity and molecular epidemiology of P. falciparum infections.

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.

qRT-PCR versus IFA-based Quantification of Male and Female Gametocytes in Low-Density Plasmodium falciparum Infections and Their Relevance for Transmission

Background

Accurate quantification of female and male gametocytes and sex ratios in asymptomatic low-density malaria infections are important for assessing their transmission potential. Gametocytes often escape detection even by molecular methods, therefore ultralow gametocyte densities were quantified in large blood volumes.

Methods

Female and male gametocytes were quantified in 161 PCR-positive Plasmodium falciparum infections from a cross-sectional survey in Papua New Guinea. Ten-fold concentrated RNA from 800 µL blood was analyzed using female-specific pfs25 and male-specific pfmget or mssp qRT-PCR. Gametocyte sex ratios from qRT-PCR were compared with those from immunofluorescence assays (IFA).

Results

Gametocytes were identified in 58% (93/161) P. falciparum-positive individuals. Mean gametocyte densities were frequently below 1 female and 1 male gametocyte/µL by qRT-PCR. The mean proportion of males was 0.39 (95% confidence interval, 0.33–0.44) by pfs25/pfmget qRT-PCR; this correlated well with IFA results (Pearsons r² = 0.91; P < .001). A Poisson model fitted to our data predicted 16% P. falciparum-positive individuals that are likely to transmit, assuming at least 1 female and 1 male gametocyte per 2.5 µL mosquito bloodmeal.

Conclusions

Based on model estimates of female and male gametocytes per 2.5 µL blood, P. falciparum-positive individuals detected exclusively by ultrasensitive diagnostics are negligible for human-to-mosquito transmission.

Estimating the transmission potential of ultralow-density malaria infections informs interventions. Almost all infections with ≥1 female and male gametocyte per 2.5 µL mosquito bloodmeal, and thus with highest likelihood of contributing to human-to-mosquito transmission, were detectable by standard molecular diagnostics.

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.

A Targeted "Next-Generation" Sequencing-Informatic Approach to Define Genetic Diversity in Theileria orientalis Populations within Individual Cattle: Proof-of-Principle

Oriental theileriosis is an economically important tickborne disease of bovines, caused by some members of the Theileria orientalis complex. Currently, 11 distinct operational taxonomic units (OTUs), or genotypes, are recognized based on their major piroplasm surface protein (MPSP) gene sequences. Two of these genotypes (i.e., chitose and ikeda) are recognized as pathogenic in cattle, causing significant disease in countries of the Asia-Pacific region. However, the true extent of genetic variation and associated virulence/pathogenicity within this complex is unknown. Here, we undertook a proof-of-principle study of a small panel of genomic DNAs (n = 13) from blood samples originating from individual cattle known to harbor T. orientalis, in order to assess the performance of a targeted “next-generation” sequencing-informatic approach to identify genotypes. Five genotypes (chitose, ikeda, buffeli, type 4, and type 5) were defined; multiple genotypes were found within individual samples, with dominant and minor sequence types representing most genotypes. This study indicates that this sequencing-informatic workflow could be useful to assess the nature and extent of genetic variation within and among populations of T. orientalis on a large scale, and to potentially employ panels of distinct gene markers for expanded molecular epidemiological investigations of socioeconomically important protistan pathogens more generally.

High Genetic Diversity of Plasmodium falciparum in the Low-Transmission Setting of the Kingdom of Eswatini

BACKGROUND

To better understand transmission dynamics, we characterized Plasmodium falciparum genetic diversity in Eswatini, where transmission is low and sustained by importation.

METHODS

Twenty-six P. falciparum microsatellites were genotyped in 66% of confirmed cases (2014-2016; N = 582). Population and within-host diversity were used to characterize differences between imported and locally acquired infections. Logistic regression was used to assess the added value of diversity metrics to classify imported and local infections beyond epidemiology data alone.

RESULTS

Parasite population in Eswatini was highly diverse (expected heterozygosity [HE] = 0.75) and complex: 67% polyclonal infections, mean multiplicity of infection (MOI) 2.2, and mean within-host infection fixation index (FWS) 0.84. Imported cases had comparable diversity to local cases but exhibited higher MOI (2.4 vs 2.0; P = .004) and lower mean FWS (0.82 vs 0.85; P = .03). Addition of MOI and FWS to multivariate analyses did not increase discrimination between imported and local infections.

CONCLUSIONS

In contrast to the common perception that P. falciparum diversity declines with decreasing transmission intensity, Eswatini isolates exhibited high parasite diversity consistent with high rates of malaria importation and limited local transmission. Estimates of malaria transmission intensity from genetic data need to consider the effect of importation, especially as countries near elimination.

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.

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.

Long-distance transmission patterns modelled from SNP barcodes of Plasmodium falciparum infections in The Gambia

Malaria has declined significantly in The Gambia and determining transmission dynamics of Plasmodium falciparum can help targeting control interventions towards elimination. This can be inferred from genetic similarity between parasite isolates from different sites and timepoints. Here, we imposed a P. falciparum life cycle time on a genetic distance likelihood model to determine transmission paths from a 54 SNP barcode of 355 isolates. Samples were collected monthly during the 2013 malaria season from six pairs of villages spanning 300 km from western to eastern Gambia. There was spatial and temporal hierarchy in pairwise genetic relatedness, with the most similar barcodes from isolates within the same households and village. Constrained by travel data, the model detected 60 directional transmission events, with 27% paths linking persons from different regions. We identified 13 infected individuals (4.2% of those genotyped) responsible for 2 to 8 subsequent infections within their communities. These super-infectors were mostly from high transmission villages. When considering paths between isolates from the most distant regions (west vs east) and travel history, there were 3 transmission paths from eastern to western Gambia, all at the peak (October) of the malaria transmission season. No paths with known travel originated from the extreme west to east. Although more than half of all paths were within-village, parasite flow from east to west may contribute to maintain transmission in western Gambia, where malaria transmission is already low. Therefore, interrupting malaria transmission in western Gambia would require targeting eastern Gambia, where malaria prevalence is substantially higher, with intensified malaria interventions.

Applying next-generation sequencing to track falciparum malaria in sub-Saharan Africa

Next-generation sequencing (NGS) technologies are increasingly being used to address a diverse range of biological and epidemiological questions. The current understanding of malaria transmission dynamics and parasite movement mainly relies on the analyses of epidemiologic data, e.g. case counts and self-reported travel history data. However, travel history data are often not routinely collected or are incomplete, lacking the necessary level of accuracy. Although genetic data from routinely collected field samples provides an unprecedented opportunity to track the spread of malaria parasites, it remains an underutilized resource for surveillance due to lack of local awareness and capacity, limited access to sensitive laboratory methods and associated computational tools and difficulty in interpreting genetic epidemiology data. In this review, the potential roles of NGS in better understanding of transmission patterns, accurately tracking parasite movement and addressing the emerging challenges of imported malaria in low transmission settings of sub-Saharan Africa are discussed. Furthermore, this review highlights the insights gained from malaria genomic research and challenges associated with integrating malaria genomics into existing surveillance tools to inform control and elimination strategies.

A pilot randomized controlled trial to compare the effectiveness of two 14-day primaquine regimens for the radical cure of vivax malaria in South India

BACKGROUND

Radical cure of Plasmodium vivax malaria requires treatment with a blood schizonticide and a hypnozoitocide (primaquine) to eradicate the dormant liver stages. There has been uncertainty about the operational effectiveness and optimum dosing of the currently recommended 14-day primaquine (PQ) course.

METHODS

A two centre, randomized, open-label, two arm study was conducted in South India. Patients were randomized to receive either high dose (0.5 mg base/kg body weight) or conventional dose (0.25 mg/kg) PQ for 14 days. Plasma concentrations of PQ and carboxyprimaquine (CPQ) on the 7th day of treatment were measured by reverse phase high performance liquid chromatography. Study subjects were followed up for 6 months. Recurrent infections were genotyped using capillary fragment length polymorphism of two PCR-amplified microsatellite markers (MS07 and MS 10).

RESULTS

Fifty patients were enrolled. Baseline characteristics and laboratory features did not differ significantly between the groups. Mean age of the study population was 42 ± 16.0 years. Recurrences 80-105 days later occurred in 4 (8%) patients, two in each the groups. All recurrences had the same microsatellite genotype as that causing the index infection suggesting all were relapses. One relapse was associated with low CPQ concentrations suggesting poor adherence.

CONCLUSIONS

This small pilot trial supports the effectiveness of the currently recommended lower dose (0.25 mg/kg/day) 14 day PQ regimen for the radical cure of vivax malaria in South India. Trial registration Clinical Trials Registry-India, CTRI/2017/03/007999. Registered 3 March 2017, http://ctri.nic.in/Clinicaltrials/regtrial.php?modid=1&compid=19&EncHid=82755.86366 .

Multiplicity and molecular epidemiology of Plasmodium vivax and Plasmodium falciparum infections in East Africa

Background

Parasite genetic diversity and multiplicity of infection (MOI) affect clinical outcomes, response to drug treatment and naturally-acquired or vaccine-induced immunity. Traditional methods often underestimate the frequency and diversity of multiclonal infections due to technical sensitivity and specificity. Next-generation sequencing techniques provide a novel opportunity to study complexity of parasite populations and molecular epidemiology.

Methods

Symptomatic and asymptomatic Plasmodium vivax samples were collected from health centres/hospitals and schools, respectively, from 2011 to 2015 in Ethiopia. Similarly, both symptomatic and asymptomatic Plasmodium falciparum samples were collected, respectively, from hospitals and schools in 2005 and 2015 in Kenya. Finger-pricked blood samples were collected and dried on filter paper. Long amplicon (> 400 bp) deep sequencing of merozoite surface protein 1 (msp1) gene was conducted to determine multiplicity and molecular epidemiology of P. vivax and P. falciparum infections. The results were compared with those based on short amplicon (117 bp) deep sequencing.

Results

A total of 139 P. vivax and 222 P. falciparum samples were pyro-sequenced for pvmsp1 and pfmsp1, yielding a total of 21 P. vivax and 99 P. falciparum predominant haplotypes. The average MOI for P. vivax and P. falciparum were 2.16 and 2.68, respectively, which were significantly higher than that of microsatellite markers and short amplicon (117 bp) deep sequencing. Multiclonal infections were detected in 62.2% of the samples for P. vivax and 74.8% of the samples for P. falciparum. Four out of the five subjects with recurrent P. vivax malaria were found to be a relapse 44–65 days after clearance of parasites. No difference was observed in MOI among P. vivax patients of different symptoms, ages and genders. Similar patterns were also observed in P. falciparum except for one study site in Kenyan lowland areas with significantly higher MOI.

Conclusions

The study used a novel method to evaluate Plasmodium MOI and molecular epidemiological patterns by long amplicon ultra-deep sequencing. The complexity of infections were similar among age groups, symptoms, genders, transmission settings (spatial heterogeneity), as well as over years (pre- vs. post-scale-up interventions). This study demonstrated that long amplicon deep sequencing is a useful tool to investigate multiplicity and molecular epidemiology of Plasmodium parasite infections.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2337-y) contains supplementary material, which is available to authorized users.

Plasmodium vivax and Plasmodium falciparum infection dynamics: re-infections, recrudescences and relapses

Background

In malaria endemic populations, complex patterns of Plasmodium vivax and Plasmodium falciparum blood-stage infection dynamics may be observed. Genotyping samples from longitudinal cohort studies for merozoite surface protein (msp) variants increases the information available in the data, allowing multiple infecting parasite clones in a single individual to be identified. msp genotyped samples from two longitudinal cohorts in Papua New Guinea (PNG) and Thailand were analysed using a statistical model where the times of acquisition and clearance of each clone in every individual were estimated using a process of data augmentation.

Results

For the populations analysed, the duration of blood-stage P. falciparum infection was estimated as 36 (95% Credible Interval (CrI): 29, 44) days in PNG, and 135 (95% CrI 94, 191) days in Thailand. Experiments on simulated data indicated that it was not possible to accurately estimate the duration of blood-stage P. vivax infections due to the lack of identifiability between a single blood-stage infection and multiple, sequential blood-stage infections caused by relapses. Despite this limitation, the method and data point towards short duration of blood-stage P. vivax infection with a lower bound of 24 days in PNG, and 29 days in Thailand. On an individual level, P. vivax recurrences cannot be definitively classified into re-infections, recrudescences or relapses, but a probabilistic relapse phenotype can be assigned to each P. vivax sample, allowing investigation of the association between epidemiological covariates and the incidence of relapses.

Conclusion

The statistical model developed here provides a useful new tool for in-depth analysis of malaria data from longitudinal cohort studies, and future application to data sets with multi-locus genotyping will allow more detailed investigation of infection dynamics.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2318-1) contains supplementary material, which is available to authorized users.

Biological concepts in recurrent Plasmodium vivax malaria

A curious aspect of the evolution of the hypnozoite theory of malarial relapse is its transmogrification from theory into 'fact', this being of historical, linguistic, scientific and sociological interest. As far as it goes, the hypnozoite explanation for relapse is almost certainly correct. I contend, however, that many of the genotypically homologous, non-reinfection, relapse-like Plasmodium vivax recurrences that researchers ascribe to hypnozoite activation are probably hypnozoite-independent. Indeed, some malariologists are starting to recognize that homologous P. vivax recurrences have most likely been overattributed to activation of hypnozoites. Hitherto identified, non-hypnozoite, possible plasmodial sources of recurrence that must be considered, besides circulating erythrocytic stages, include parasites in splenic dendritic cells, other cells in the spleen (in addition to infected erythrocytes there), bone marrow (importantly) and the skin. I argue that we need to take into account the possibility of a dual or multiple extra-vascular origin of P. vivax non-reinfection recurrences, not arbitrarily discount it. The existence of a P. vivax reservoir(s) is a topical subject and one of practical importance for malaria eradication. Pertinent drug-associated matters are also discussed, as is the dormancy-related significance of clues provided by blood-stage-induced malarial infection.