A deep sequencing tool for partitioning clearance rates following antimalarial treatment in polyclonal 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.

Targeted amplicon deep sequencing of ama1 and mdr1 to track within-host P. falciparum diversity throughout treatment in a clinical drug trial

Introduction

Antimalarial therapeutic efficacy studies are routinely conducted in malaria-endemic countries to assess the effectiveness of antimalarial treatment strategies. Targeted amplicon sequencing (AmpSeq) uniquely identifies and quantifies genetically distinct parasites within an infection. In this study, AmpSeq of Plasmodium falciparum apical membrane antigen 1 ( ama1), and multidrug resistance gene 1 ( mdr1), were used to characterise the complexity of infection (COI) and drug-resistance genotypes, respectively.

Methods

P. falciparum-positive samples were obtained from a triple artemisinin combination therapy clinical trial conducted in 30 children under 13 years of age between 2018 and 2019 in Kilifi, Kenya. Nine of the 30 participants presented with recurrent parasitemia from day 26 (624h) onwards. The ama1 and mdr1 genes were amplified and sequenced, while msp1, msp2 and glurp data were obtained from the original clinical study.

Results

The COI was comparable between ama1 and msp1, msp2 and glurp; overall, ama1 detected more microhaplotypes. Based on ama1, a stable number of microhaplotypes were detected throughout treatment until day 3. Additionally, a recrudescent infection was identified with an ama1 microhaplotype initially observed at 30h and later in an unscheduled follow-up visit. Using the relative frequencies of ama1 microhaplotypes and parasitemia, we identified a fast (<1h) and slow (>5h) clearing microhaplotype. As expected, only two mdr1 microhaplotypes (NF and NY) were identified based on the combination of amino acid polymorphisms at codons 86 and 184.

Conclusions

This study highlights AmpSeq as a tool for highly-resolution tracking of parasite microhaplotypes throughout treatment and can detect variation in microhaplotype clearance estimates. AmpSeq can also identify slow-clearing microhaplotypes, a potential early sign of selection during treatment. Consequently, AmpSeq has the capability of improving the discriminatory power to distinguish recrudescences from reinfections accurately.

Dynamics of parasite growth in genetically diverse Plasmodium falciparum isolates

Multiple parasite lineages with different proliferation rates or fitness may coexist within a clinical malaria isolate, resulting in complex growth interactions and variations in phenotype. To elucidate the dynamics of parasite growth in multiclonal isolates, we measured growth rates (GRs) of three Plasmodium falciparum Cambodian isolates, including IPC_3445 (MRA-1236), IPC_5202 (MRA-1240), IPC_6403 (MRA-1285), and parasite lineages previously cloned from each of these isolates by limiting dilution. Following synchronization, in vitro cultures of each parasite line were maintained over four consecutive asexual cycles (192 h), with thin smears prepared at each 48-h cycle to estimate GR and fold change in parasitemia (FCP). Cell cycle time (CCT), the duration it takes for ring-stage parasites to develop into mature schizonts, was measured by monitoring the development of 0-3-h post-invasion rings for up to 52 h post-incubation. Laboratory lines 3D7 (MRA-102) and Dd2 (MRA-150) were used as controls. Significant differences in GR, FCP, and CCT were observed between parasite isolates and clonal lineages from each isolate. The parasite lines studied here have well-defined growth phenotypes and will facilitate basic malaria research and development of novel malaria interventions. These lines are available to malaria researchers through the MR4 collection of NIAID's BEI Resources Program.

Low Complexity of Infection Is Associated With Molecular Persistence of Plasmodium falciparum in Kenya and Tanzania

Background

Plasmodium falciparum resistance to artemisinin-based combination therapies (ACTs) is a threat to malaria elimination. ACT-resistance in Asia raises concerns for emergence of resistance in Africa. While most data show high efficacy of ACT regimens in Africa, there have been reports describing declining efficacy, as measured by both clinical failure and prolonged parasite clearance times.

Methods

Three hundred children aged 2-10 years with uncomplicated P. falciparum infection were enrolled in Kenya and Tanzania after receiving treatment with artemether-lumefantrine. Blood samples were taken at 0, 24, 48, and 72 h, and weekly thereafter until 28 days post-treatment. Parasite and host genetics were assessed, as well as clinical, behavioral, and environmental characteristics, and host anti-malarial serologic response.

Results

While there was a broad range of clearance rates at both sites, 85% and 96% of Kenyan and Tanzanian samples, respectively, were qPCR-positive but microscopy-negative at 72 h post-treatment. A greater complexity of infection (COI) was negatively associated with qPCR-detectable parasitemia at 72 h (OR: 0.70, 95% CI: 0.53-0.94), and a greater baseline parasitemia was marginally associated with qPCR-detectable parasitemia (1,000 parasites/uL change, OR: 1.02, 95% CI: 1.01-1.03). Demographic, serological, and host genotyping characteristics showed no association with qPCR-detectable parasitemia at 72 h. Parasite haplotype-specific clearance slopes were grouped around the mean with no association detected between specific haplotypes and slower clearance rates.

Conclusions

Identifying risk factors for slow clearing P. falciparum infections, such as COI, are essential for ongoing surveillance of ACT treatment failure in Kenya, Tanzania, and more broadly in sub-Saharan Africa.

Targeted Amplicon deep sequencing of ama1 and mdr1 to track within-host P. falciparum diversity throughout treatment in a clinical drug trial

Antimalarial therapeutic efficacy studies are routinely conducted in malaria-endemic countries to assess the effectiveness of antimalarial treatment strategies. Targeted amplicon deep sequencing (TADS) uniquely identifies and quantifies genetically distinct parasites within an infection. In this study, TADS Plasmodium falciparum apical membrane antigen 1 (ama1), and multidrug resistance gene 1 (mdr1), were used to characterize the complexity of infection (COI) and drug-resistance genotypes, respectively. P. falciparum positive samples were obtained from a triple artemisinin combination therapy clinical trial conducted in 30 children under 13 years of age between 2018 and 2019 in Kilifi, Kenya. Of the 30 participants, 9 presented with recurrent parasitemia from day 26 (624h) onwards. The ama1 and mdr1 genes were amplified and sequenced, while msp1, msp2 and glurp data were obtained from the original clinical study. The COI was comparable between ama1 and msp1, msp2 and glurp, however, overall ama1 detected more haplotypes. Based on ama1, a stable number of haplotypes were detected throughout treatment up until day 3. Additionally, a recrudescent infection was identified with an ama1 haplotype initially observed at 30h and later in an unscheduled follow-up visit. Using the relative frequencies of ama1 haplotypes and parasitaemia, we identified a fast (<1h) and slow (>5h) clearing haplotype. As expected, only two mdr1 haplotypes (NF and NY) were identified based on the combination of amino acid polymorphisms at codons 86 and 184. This study highlights TADS as a sensitive tool for tracking parasite haplotypes throughout treatment and can detect variation in haplotype clearance estimates. TADS can also identify slow clearing haplotypes, a potential early sign of selection during treatment. Consequently, TADS has the capability of improving the discriminatory power to accurately distinguish recrudescences from reinfections.

Potential Emergence of Plasmodium Resistance to Artemisinin Induced by the Use of Artemisia annua for Malaria and COVID-19 Prevention in Sub-African Region

Plasmodium resistance to antimalarial drugs is an obstacle to the elimination of malaria in endemic areas. This situation is particularly dramatic for Africa, which accounts for nearly 92% of malaria cases worldwide. Drug pressure has been identified as a key factor in the emergence of antimalarial drug resistance. Indeed, this pressure is favoured by several factors, including the use of counterfeit forms of antimalarials, inadequate prescription controls, poor adherence to treatment regimens, dosing errors, and the increasing use of other forms of unapproved antimalarials. This resistance has led to the replacement of chloroquine (CQ) by artemisinin-based combination therapies (ACTs) which are likely to become ineffective in the coming years due to the uncontrolled use of Artemisia annua in the sub-Saharan African region for malaria prevention and COVID-19. The use of Artemisia annua for the prevention of malaria and COVID-19 could be an important factor in the emergence of resistance to Artemisinin-based combination therapies.

Safety and efficacy of intermittent presumptive treatment with sulfadoxine-pyrimethamine using rapid diagnostic test screening and treatment with dihydroartemisinin-piperaquine at the first antenatal care visit (IPTp-SP+): study protocol for a randomized controlled trial

Background

Intermittent preventive treatment in pregnancy (IPTp) with sulfadoxine-pyrimethamine (SP) is recommended by the World Health Organization for the prevention of malaria in pregnancy (MIP)-associated adverse outcomes in high burden areas. However, the efficacy of IPTp-SP has decreased in step with increasing parasite drug resistance. Suitable alternative strategies are needed.

Methods

This is a protocol for a phase IIIb open-label, two-armed randomized controlled superiority trial to assess the safety and efficacy of a hybrid approach to IPTp combining screening and treatment with dihydroartemisinin-piperaquine (DP) to the current IPTp-SP regimen at the first antenatal care clinic visit. Pregnant women without HIV infection and without signs or symptoms of malaria will be randomized to either standard IPTp-SP or hybrid IPTp-SP plus screening and treatment (IPTp-SP+). In the IPTp-SP+ arm, participants who screen positive by rapid diagnostic test for P. falciparum will be treated with DP at the first antenatal visit while those who screen negative will receive SP per current guidelines. All participants will be administered SP on days 35 and 63 and will be actively followed biweekly up to day 63 and then monthly until delivery. Infants will be followed until 1 year after delivery. The primary endpoint is incident PCR-confirmed MIP at day 42. Secondary endpoints include incident MIP at other time points, placental malaria, congenital malaria, hemoglobin trends, birth outcomes, and incidence of adverse events in infants up to the first birthday.

Discussion

A hybrid approach to IPTp that combines screening and treatment with an artemisinin-based combination therapy at the first visit with standard IPTp-SP is hypothesized to confer added benefit over IPTp-SP alone in a high malaria transmission area with prevalent SP resistant parasites.

Trial registration

Pan African Clinical Trials Registry 201905721140808. Registered retrospectively on 11 May 2019

Supplementary Information

The online version contains supplementary material available at 10.1186/s13063-021-05745-0.

Immune selection suppresses the emergence of drug resistance in malaria parasites but facilitates its spread

Although drug resistance in Plasmodium falciparum typically evolves in regions of low transmission, resistance spreads readily following introduction to regions with a heavier disease burden. This suggests that the origin and the spread of resistance are governed by different processes, and that high transmission intensity specifically impedes the origin. Factors associated with high transmission, such as highly immune hosts and competition within genetically diverse infections, are associated with suppression of resistant lineages within hosts. However, interactions between these factors have rarely been investigated and the specific relationship between adaptive immunity and selection for resistance has not been explored. Here, we developed a multiscale, agent-based model of Plasmodium parasites, hosts, and vectors to examine how host and parasite dynamics shape the evolution of resistance in populations with different transmission intensities. We found that selection for antigenic novelty (“immune selection”) suppressed the evolution of resistance in high transmission settings. We show that high levels of population immunity increased the strength of immune selection relative to selection for resistance. As a result, immune selection delayed the evolution of resistance in high transmission populations by allowing novel, sensitive lineages to remain in circulation at the expense of the spread of a resistant lineage.

In contrast, in low transmission settings, we observed that resistant strains were able to sweep to high population prevalence without interference. Additionally, we found that the relationship between immune selection and resistance changed when resistance was widespread. Once resistance was common enough to be found on many antigenic backgrounds, immune selection stably maintained resistant parasites in the population by allowing them to proliferate, even in untreated hosts, when resistance was linked to a novel epitope. Our results suggest that immune selection plays a role in the global pattern of resistance evolution.

Drug resistance in the malaria parasite, Plasmodium falciparum, presents an ongoing public health challenge, but aspects of its evolution are poorly understood. Although antimalarial resistance is common worldwide, it can typically be traced to just a handful of evolutionary origins. Counterintuitively, although Sub Saharan Africa bears 90% of the global malaria burden, resistance typically originates in regions where transmission intensity is low. In high transmission regions, infections are genetically diverse, and hosts have significant standing adaptive immunity, both of which are known to suppress the frequency of resistance within infections. However, interactions between immune-driven selection, transmission intensity, and resistance have not been investigated. Using a multiscale, agent-based model, we found that high transmission intensity slowed the evolution of resistance via its effect on host population immunity. High host immunity strengthened selection for antigenic novelty, interfering with selection for resistance and allowing sensitive lineages to suppress resistant lineages in untreated hosts. However, once resistance was common in the circulating parasite population, immune selection maintained it in the population at a high prevalence. Our findings provide a novel explanation for observations about the origin of resistance and suggest that adaptive immunity is a critical component of selection.

Dissection of haplotype-specific drug response phenotypes in multiclonal malaria isolates

Natural infections of Plasmodium falciparum, the parasite responsible for the deadliest form of human malaria, often comprise multiple parasite lineages (haplotypes). Multiclonal parasite isolates may exhibit variable phenotypes including different drug susceptibility profiles over time due to the presence of multiple haplotypes. To test this hypothesis, three P. falciparum Cambodian isolates IPC_3445 (MRA-1236), IPC_5202 (MRA-1240) and IPC_6403 (MRA-1285) suspected to be multiclonal were cloned by limiting dilution, and the resulting clones genotyped at 24 highly polymorphic single nucleotide polymorphisms (SNPs). Isolates harbored up to three constituent haplotypes, and exhibited significant variability (p < 0.05) in susceptibility to chloroquine, mefloquine, artemisinin and piperaquine as measured by half maximal drug inhibitory concentration (IC₅₀) assays and parasite survival assays, which measure viability following exposure to pharmacologically relevant concentrations of antimalarial drugs. The IC₅₀ of the most abundant haplotype frequently reflected that of the uncloned parental isolate, suggesting that a single haplotype dominates the antimalarial susceptibility profile and masks the effect of minor frequency haplotypes. These results indicate that phenotypic variability in parasite isolates is often due to the presence of multiple haplotypes. Depending on intended end-use, clinical isolates should be cloned to yield single parasite lineages with well-defined phenotypes and genotypes. The availability of such standardized clonal parasite lineages through NIAID's BEI Resources program will aid research directed towards the development of diagnostics and interventions including drugs against malaria.

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.

Parasite clearance, cure rate, post-treatment prophylaxis and safety of standard 3-day versus an extended 6-day treatment of artemether-lumefantrine and a single low-dose primaquine for uncomplicated Plasmodium falciparum malaria in Bagamoyo district, Tanzania: a randomized controlled trial

Background

Artemisinin-based combination therapy (ACT) resistant Plasmodium falciparum represents an increasing threat to Africa. Extended ACT regimens from standard 3 to 6 days may represent a means to prevent its development and potential spread in Africa.

Methods

Standard 3-day treatment with artemether–lumefantrine (control) was compared to extended 6-day treatment and single low-dose primaquine (intervention); in a randomized controlled, parallel group, superiority clinical trial of patients aged 1–65 years with microscopy confirmed uncomplicated P. falciparum malaria, enrolled in Bagamoyo district, Tanzania. The study evaluated parasite clearance, including proportion of PCR detectable P. falciparum on days 5 and 7 (primary endpoint), cure rate, post-treatment prophylaxis, safety and tolerability. Clinical, and laboratory assessments, including ECG were conducted during 42 days of follow-up. Blood samples were collected for parasite detection (by microscopy and PCR), molecular genotyping and pharmacokinetic analyses. Kaplan–Meier survival analyses were done for both parasite clearance and recurrence.

Results

A total of 280 patients were enrolled, 141 and 139 in the control and intervention arm, respectively, of whom 121 completed 42 days follow-up in each arm. There was no difference in proportion of PCR positivity across the arms at day 5 (80/130 (61.5%) vs 89/134 (66.4%), p = 0.44), or day 7 (71/129 (55.0%) vs 70/134 (52.2%), p = 0.71). Day 42 microscopy determined cure rates (PCR adjusted) were 97.4% (100/103) and 98.3% (110/112), p = 0.65, in the control and intervention arm, respectively. Microscopy determined crude recurrent parasitaemia during follow-up was 21/121 (17.4%) in the control and 14/121 (11.6%) in the intervention arm, p = 0.20, and it took 34 days and 42 days in the respective arms for 90% of the patients to remain without recurrent parasitaemia. Lumefantrine exposure was significantly higher in intervention arm from D3 to D42, but cardiac, biochemical and haematological safety was high and similar in both arms.

Conclusion

Extended 6-day artemether–lumefantrine treatment and a single low-dose of primaquine was not superior to standard 3-day treatment for ACT sensitive P. falciparum infections but, importantly, equally efficacious and safe. Thus, extended artemether–lumefantrine treatment may be considered as a future treatment regimen for ACT resistant P. falciparum, to prolong the therapeutic lifespan of ACT in Africa.

Trial registration ClinicalTrials.gov, NCT03241901. Registered July 27, 2017 https://clinicaltrials.gov/show/NCT03241901

Evaluation of residual submicroscopic Plasmodium falciparum parasites 3 days after initiation of treatment with artemisinin-based combination therapy

Plasmodium falciparum resistance against artemisinin has not emerged in Africa; however, there are reports of the presence of polymerase chain reaction-determined residual submicroscopic parasitaemia detected on day 3 after artemisinin-based combination therapy (ACT). These residual submicroscopic parasites are thought to represent tolerant/resistant parasites against artemisinin, the fast-acting component of the combination. This review focused on residual submicroscopic parasitaemia, what it represents, and its significance on the emergence and spread of artemisinin resistance in Africa. Presence of residual submicroscopic parasitemia on day 3 after treatment initiation leaves question on whether successful treatment is attained with ACT. Thus there is a need to determine the potential public health implication of the PCR-determined residual submicroscopic parasitaemia observed on day 3 after ACT. Robust techniques, such as in vitro cultivation, should be used to evaluate if the residual submicroscopic parasites detected on day 3 after ACT are viable asexual parasites, or gametocytes, or the DNA of the dead parasites waiting to be cleared from the circulation. Such techniques would also evaluate the transmissibility of these residual parasites.

PCR correction strategies for malaria drug trials: updates and clarifications

Malaria drug trials conducted in endemic areas face a major challenge in their analysis because it is difficult to establish whether parasitaemia in blood samples collected after treatment indicate drug failure or a new infection acquired after treatment. It is therefore vital to reliably distinguish drug failures from new infections in order to obtain accurate estimates of drug failure rates. This distinction can be achieved for Plasmodium falciparum by comparing parasite genotypes obtained at the time of treatment (the baseline) and on the day of recurring parasitaemia. Such PCR correction is required to obtain accurate failure rates, even for new effective drugs. Despite the routine use of PCR correction in surveillance of drug resistance and in clinical drug trials, limitations inherent to the molecular genotyping methods have led some researchers to question the validity of current PCR correction strategies. Here we describe and discuss recent developments in these genotyping approaches, with a particular focus on method validation and limitations of the genotyping strategies. Our aim is to update scientists from public and private bodies who are working on the development, deployment, and surveillance of new malaria drugs. We aim to promote discussion around these issues and argue for the adoption of improved standardised PCR correction methodologies.

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.

Recrudescence, Reinfection, or Relapse? A More Rigorous Framework to Assess Chloroquine Efficacy for Plasmodium vivax Malaria

Background

Plasmodium vivax resistance to chloroquine (CQ) has been reported worldwide, although the World Health Organization clinical drug efficacy studies protocol does not permit classification of patient outcomes.

Methods

We enrolled 40 patients with P. vivax malaria in northeastern Cambodia, where >17% treatment failures were previously reported. Patients were treated with CQ (30 mg/kg) and followed for 2 months, with frequent clinical examination and capillary blood sample collection for microscopy, molecular parasite detection and genotyping, and drug concentration measurements. Reinfections were prevented by relocating patients to a transmission-free area.

Results

P. vivax parasites were eliminated in all patients by day 3. Genomic analyses revealed that all clones in polyclonal infections were cleared at the same rate, indicating their equal susceptibility to CQ. CQ blood concentrations were below the therapeutic level in all recurrent infections (24 of 40 patients), which were efficiently cleared by a second course of CQ treatment. Genotyping (128 SNPs barcode) and sequences of entire parasite genome (Whole-Genome Sequencing, Illumina) indicated that two thirds (6 of 8) of the recurrent parasites resulted from heterologous relapses whose 50% are from by sibling/recombinant clones.

Conclusions

No evidence of CQ resistance was observed. Our data suggest that P. vivax antimalarial drug resistance is likely overestimated and that the current guidelines for clinical drug studies of P. vivax malaria need to be revised.

Mosquitoes as a feasible sentinel group for anti-malarial resistance surveillance by Next Generation Sequencing of Plasmodium falciparum

Background

Plasmodium falciparum drug resistance surveillance is key to successful disease control and eradication. Contemporary methods that only allow determination of prevalence of resistance are expensive, time consuming and require ethical considerations. A newer method involving Next Generation Sequencing (NGS) permits obtaining frequency of resistance while allowing to detect minority variants in mixed infections. Here, NGS was tested for P. falciparum resistance marker detection in mosquito samples as a feasible and suitable alternative for molecular resistance surveillance. Anopheles funestus were collected in southern Mozambique using CDC light traps and manual collections. DNA was extracted from either whole mosquito, head-thorax and abdomen separately or pools of five mosquitoes. These samples were screened for P. falciparum and if positive for k13, pfcrt, pfmdr1, pfdhps and pfdhfr mutations related to anti-malarial drug resistance with Sanger sequencing and NGS.

Results

Among the 846 samples screened for P. falciparum, 122 were positive by 18S ssrDNA qPCR with an infection rate of 23.6%. No mutations were observed for k13 and pfcrt72-76 and almost zero for pfmdr86, but quintuple pfdhfr/pfdhps mutations were near fixation and about half of the isolates contained the pfmdr184F polymorphism. Similar allele frequencies of resistance markers were estimated with NGS in comparison with the prevalence of markers obtained with the gold standard Sanger sequencing.

Conclusions

Pooled deep sequencing of P. falciparum isolates extracted from mosquitoes is a promising, efficient and cost-effective method to quantify allele frequencies at population level which allows to detect known and unknown markers of resistance in single and mixed infections in a timelier manner. Using mosquitoes as sentinel group and focusing on allele frequency opposed to prevalence, permits active surveillance across a more homogeneous geographical range.

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

BACKGROUND

Deep sequencing of targeted genomic regions is becoming a common tool for understanding the dynamics and complexity of Plasmodium infections, but its lower limit of detection is currently unknown. Here, a new amplicon analysis tool, the Parallel Amplicon Sequencing Error Correction (PASEC) pipeline, is used to evaluate the performance of amplicon sequencing on low-density Plasmodium DNA samples. Illumina-based sequencing of two Plasmodium falciparum genomic regions (CSP and SERA2) was performed on two types of samples: in vitro DNA mixtures mimicking low-density infections (1-200 genomes/μl) and extracted blood spots from a combination of symptomatic and asymptomatic individuals (44-653,080 parasites/μl). Three additional analysis tools-DADA2, HaplotypR, and SeekDeep-were applied to both datasets and the precision and sensitivity of each tool were evaluated.

RESULTS

Amplicon sequencing can contend with low-density samples, showing reasonable detection accuracy down to a concentration of 5 Plasmodium genomes/μl. Due to increased stochasticity and background noise, however, all four tools showed reduced sensitivity and precision on samples with very low parasitaemia (< 5 copies/μl) or low read count (< 100 reads per amplicon). PASEC could distinguish major from minor haplotypes with an accuracy of 90% in samples with at least 30 Plasmodium genomes/μl, but only 61% at low Plasmodium concentrations (< 5 genomes/μl) and 46% at very low read counts (< 25 reads per amplicon). The four tools were additionally used on a panel of extracted parasite-positive blood spots from natural malaria infections. While all four identified concordant patterns of complexity of infection (COI) across four sub-Saharan African countries, COI values obtained for individual samples differed in some cases.

CONCLUSIONS

Amplicon deep sequencing can be used to determine the complexity and diversity of low-density Plasmodium infections. Despite differences in their approach, four state-of-the-art tools resolved known haplotype mixtures with similar sensitivity and precision. Researchers can therefore choose from multiple robust approaches for analysing amplicon data, however, error filtration approaches should not be uniformly applied across samples of varying parasitaemia. Samples with very low parasitaemia and very low read count have higher false positive rates and call for read count thresholds that are higher than current default recommendations.

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

Longitudinal tracking of individual Plasmodium falciparum strains in multi-clonal infections is essential for investigating infection dynamics of malaria. The traditional genotyping techniques did not permit tracking changes in individual clone density during persistent natural infections. Amplicon deep sequencing (Amp-Seq) offers a tool to address this knowledge gap. The sensitivity of Amp-Seq for relative quantification of clones was investigated using three molecular markers, ama1-D2, ama1-D3, and cpmp. Amp-Seq and length-polymorphism based genotyping were compared for their performance in following minority clones in longitudinal samples from Papua New Guinea. Amp-Seq markers were superior to length-polymorphic marker msp2 in detecting minority clones (sensitivity Amp-Seq: 95%, msp2: 85%). Multiplicity of infection (MOI) by Amp-Seq was 2.32 versus 1.73 for msp2. The higher sensitivity had no effect on estimates of force of infection because missed minority clones were detected in preceding or succeeding bleeds. Individual clone densities were tracked longitudinally by Amp-Seq despite MOI > 1, thus providing an additional parameter for investigating malaria infection dynamics. Amp-Seq based genotyping of longitudinal samples improves detection of minority clones and estimates of MOI. Amp-Seq permits tracking of clone density over time to study clone competition or the dynamics of specific, i.e. resistance-associated genotypes.

Intra-host dynamics of co-infecting parasite genotypes in asymptomatic malaria patients

Malaria-infected individuals often harbor mixtures of genetically distinct parasite genotypes. We studied intra-host dynamics of parasite genotypes co-infecting asymptomatic adults in an area of intense malaria transmission in Chikhwawa, Malawi. Serial blood samples (5 ml) were collected over seven consecutive days from 25 adults with asymptomatic Plasmodium falciparum malaria and analyzed to determine whether a single peripheral blood sample accurately captures within-host parasite diversity. Blood samples from three of the participants were also analyzed by limiting dilution cloning and SNP genotyping of the parasite clones isolated to examine both the number and relatedness of co-infecting parasite haplotypes. We observed rapid turnover of co-infecting parasite genotypes in 88% of the individuals sampled (n = 22) such that the genetic composition of parasites infecting these individuals changed dramatically over the course of seven days of follow up. Nineteen of the 25 individuals sampled (76%) carried multiple parasite genotypes at baseline. Analysis of serial blood samples from three of the individuals revealed that they harbored 6, 12 and 17 distinct parasite haplotypes respectively. Approximately 70% of parasite haplotypes recovered from the three extensively sampled individuals were unrelated (proportion of shared alleles <83.3%) and were deemed to have primarily arisen from superinfection (inoculation of unrelated parasite haplotypes through multiple mosquito bites). The rest were related at the half-sib level or greater and were deemed to have been inoculated into individual human hosts via parasite co-transmission from single mosquito bites. These findings add further to the growing weight of evidence indicating that a single blood sample poorly captures within-host parasite diversity and underscore the importance of repeated blood sampling to accurately capture within-host parasite ecology. Our data also demonstrate a more pronounced role for parasite co-transmission in generating within-host parasite diversity in high transmission settings than previously assumed. Taken together, these findings have important implications for understanding the evolution of drug resistance, malaria transmission, parasite virulence, allocation of gametocyte sex ratios and acquisition of malaria immunity.

Molecular assays for antimalarial drug resistance surveillance: A target product profile

Antimalarial drug resistance is a major constraint for malaria control and elimination efforts. Artemisinin-based combination therapy is now the mainstay for malaria treatment. However, delayed parasite clearance following treatment with artemisinin derivatives has now spread in the Greater Mekong Sub region and may emerge or spread to other malaria endemic regions. This spread is of great concern for malaria control programmes, as no alternatives to artemisinin-based combination therapies are expected to be available in the near future. There is a need to strengthen surveillance systems for early detection and response to the antimalarial drug resistance threat. Current surveillance is mainly done through therapeutic efficacy studies; however these studies are complex and both time- and resource-intensive. For multiple common antimalarials, parasite drug resistance has been correlated with specific genetic mutations, and the molecular markers associated with antimalarial drug resistance offer a simple and powerful tool to monitor the emergence and spread of resistant parasites. Different techniques to analyse molecular markers associated with antimalarial drug resistance are available, each with advantages and disadvantages. However, procedures are not adequately harmonized to facilitate comparisons between sites. Here we describe the target product profiles for tests to analyse molecular markers associated with antimalarial drug resistance, discuss how use of current techniques can be standardised, and identify the requirements for an ideal product that would allow malaria endemic countries to provide useful spatial and temporal information on the spread of resistance.

Fitness Costs and the Rapid Spread of kelch13-C580Y Substitutions Conferring Artemisinin Resistance

Fitness costs are key determinants of whether drug resistance alleles establish and how fast they spread within populations. More than 125 different kelch13 alleles, each containing a different amino acid substitution, have arisen in Southeast Asian malaria parasite (Plasmodium falciparum) populations under artemisinin selection over the past 15 years in a dramatic example of a soft selective event.

Fitness costs are key determinants of whether drug resistance alleles establish and how fast they spread within populations. More than 125 different kelch13 alleles, each containing a different amino acid substitution, have arisen in Southeast Asian malaria parasite (Plasmodium falciparum) populations under artemisinin selection over the past 15 years in a dramatic example of a soft selective event. However, just one of these alleles (C580Y) is now outcompeting other alleles in multiple different countries and is spreading toward fixation. Here we examine the fitness consequences of C580Y, relative to another less successful kelch13 mutation (R561H), to try to explain the distinctive dynamics of C580Y. We hypothesized that C580Y will show lower fitness costs than other kelch13 substitutions in the absence of artemisinin treatment. We used CRISPR/Cas9 methods to introduce single mutations (C580Y or R561H) or synonymous control edits into a wild-type parasite isolated on the Thailand-Myanmar border, conducted replicated head-to-head competition assays, and determined the outcome of competition using deep sequencing of kelch13 amplicons. Contrary to our predictions, these experiments reveal that C580Y carries higher fitness costs (s [selection coefficient] = 0.15 ± 0.008 [1 standard error {SE}]) than R561H (s = 0.084 ± 0.005). Furthermore, R561H outcompetes C580Y in direct competition (s = 0.065 ± 0.004). We conclude that fitness costs of C580Y in isolation are unlikely to explain the rapid spread of this substitution.

The impact of within-host ecology on the fitness of a drug-resistant parasite

Background and objectives

The rate of evolution of drug resistance depends on the fitness of resistant pathogens. The fitness of resistant pathogens is reduced by competition with sensitive pathogens in untreated hosts and so enhanced by competitive release in drug-treated hosts. We set out to estimate the magnitude of those effects on a variety of fitness measures, hypothesizing that competitive suppression and competitive release would have larger impacts when resistance was rarer to begin with.

Methodology

We infected mice with varying densities of drug-resistant Plasmodium chabaudi malaria parasites in a fixed density of drug-sensitive parasites and followed infection dynamics using strain-specific quantitative PCR.

Results

Competition with susceptible parasites reduced the absolute fitness of resistant parasites by 50–100%. Drug treatment increased the absolute fitness from 2- to >10 000-fold. The ecological context and choice of fitness measure was responsible for the wide variation in those estimates. Initial population growth rates poorly predicted parasite abundance and transmission probabilities.

Conclusions and implications

(i) The sensitivity of estimates of pathogen fitness to ecological context and choice of fitness measure make it difficult to derive field-relevant estimates of the fitness costs and benefits of resistance from experimental settings. (ii) Competitive suppression can be a key force preventing resistance from emerging when it is rare, as it is when it first arises. (iii) Drug treatment profoundly affects the fitness of resistance. Resistance evolution could be slowed by developing drug use policies that consider in-host competition.

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.

Altered life history strategies protect malaria parasites against drugs

Drug resistance has been reported against all antimalarial drugs, and while parasites can evolve classical resistance mechanisms (e.g., efflux pumps), it is also possible that changes in life history traits could help parasites evade the effects of treatment. The life history of malaria parasites is governed by an intrinsic resource allocation problem: specialized stages are required for transmission, but producing these stages comes at the cost of producing fewer of the forms required for within-host survival. Drug treatment, by design, alters the probability of within-host survival, and so should alter the costs and benefits of investing in transmission. Here, we use a within-host model of malaria infection to predict optimal patterns of investment in transmission in the face of different drug treatment regimes and determine the extent to which alternative patterns of investment can buffer the fitness loss due to drugs. We show that over a range of drug doses, parasites are predicted to adopt "reproductive restraint" (investing more in asexual replication and less in transmission) to maximize fitness. By doing so, parasites recoup some of the fitness loss imposed by drugs, though as may be expected, increasing dose reduces the extent to which altered patterns of transmission investment can benefit parasites. We show that adaptation to drug-treated infections could result in more virulent infections in untreated hosts. This work emphasizes that in addition to classical resistance mechanisms, drug treatment generates selection for altered parasite life history. Understanding how any shifts in life history will alter the efficacy of drugs, as well as any limitations on such shifts, is important for evaluating and predicting the consequences of drug treatment.

Effects of multiple sources of genetic drift on pathogen variation within hosts

Changes in pathogen genetic variation within hosts alter the severity and spread of infectious diseases, with important implications for clinical disease and public health. Genetic drift may play a strong role in shaping pathogen variation, but analyses of drift in pathogens have oversimplified pathogen population dynamics, either by considering dynamics only at a single scale-such as within hosts or between hosts-or by making drastic simplifying assumptions, for example, that host immune systems can be ignored or that transmission bottlenecks are complete. Moreover, previous studies have used genetic data to infer the strength of genetic drift, whereas we test whether the genetic drift imposed by pathogen population processes can be used to explain genetic data. We first constructed and parameterized a mathematical model of gypsy moth baculovirus dynamics that allows genetic drift to act within and between hosts. We then quantified the genome-wide diversity of baculovirus populations within each of 143 field-collected gypsy moth larvae using Illumina sequencing. Finally, we determined whether the genetic drift imposed by host-pathogen population dynamics in our model explains the levels of pathogen diversity in our data. We found that when the model allows drift to act at multiple scales-including within hosts, between hosts, and between years-it can accurately reproduce the data, but when the effects of drift are simplified by neglecting transmission bottlenecks and stochastic variation in virus replication within hosts, the model fails. A de novo mutation model and a purifying selection model similarly fail to explain the data. Our results show that genetic drift can play a strong role in determining pathogen variation and that mathematical models that account for pathogen population growth at multiple scales of biological organization can be used to explain this variation.

SeekDeep: single-base resolution de novo clustering for amplicon deep sequencing

PCR amplicon deep sequencing continues to transform the investigation of genetic diversity in viral, bacterial, and eukaryotic populations. In eukaryotic populations such as Plasmodium falciparum infections, it is important to discriminate sequences differing by a single nucleotide polymorphism. In bacterial populations, single-base resolution can provide improved resolution towards species and strains. Here, we introduce the SeekDeep suite built around the qluster algorithm, which is capable of accurately building de novo clusters representing true, biological local haplotypes differing by just a single base. It outperforms current software, particularly at low frequencies and at low input read depths, whether resolving single-base differences or traditional OTUs. SeekDeep is open source and works with all major sequencing technologies, making it broadly useful in a wide variety of applications of amplicon deep sequencing to extract accurate and maximal biologic information.

A deep sequencing approach to estimate Plasmodium falciparum complexity of infection (COI) and explore apical membrane antigen 1 diversity

BACKGROUND

Humans living in regions with high falciparum malaria transmission intensity harbour multi-strain infections comprised of several genetically distinct malaria haplotypes. The number of distinct malaria parasite haplotypes identified from an infected human host at a given time is referred to as the complexity of infection (COI). In this study, an amplicon-based deep sequencing method targeting the Plasmodium falciparum apical membrane antigen 1 (pfama1) was utilized to (1) investigate the relationship between P. falciparum prevalence and COI, (2) to explore the population genetic structure of P. falciparum parasites from malaria asymptomatic individuals participating in the 2007 Demographic and Health Survey (DHS) in the Democratic Republic of Congo (DRC), and (3) to explore selection pressures on geospatially divergent parasite populations by comparing AMA1 amino acid frequencies in the DRC and Mali.

RESULTS

A total of 900 P. falciparum infections across 11 DRC provinces were examined. Deep sequencing of both individuals, for COI analysis, and pools of individuals, to examine population structure, identified 77 unique pfama1 haplotypes. The majority of individual infections (64.5%) contained polyclonal (COI > 1) malaria infections based on the presence of genetically distinct pfama1 haplotypes. A minimal correlation between COI and malaria prevalence as determined by sensitive real-time PCR was identified. Population genetic analyses revealed extensive haplotype diversity, the vast majority of which was shared across the sites. AMA1 amino acid frequencies were similar between parasite populations in the DRC and Mali.

CONCLUSIONS

Amplicon-based deep sequencing is a useful tool for the detection of multi-strain infections that can aid in the understanding of antigen heterogeneity of potential malaria vaccine candidates, population genetics of malaria parasites, and factors that influence complex, polyclonal malaria infections. While AMA1 and other diverse markers under balancing selection may perform well for understanding COI, they may offer little geographic or temporal discrimination between parasite populations.

Development of amplicon deep sequencing markers and data analysis pipeline for genotyping multi-clonal malaria infections

Background

Amplicon deep sequencing permits sensitive detection of minority clones and improves discriminatory power for genotyping multi-clone Plasmodium falciparum infections. New amplicon sequencing and data analysis protocols are needed for genotyping in epidemiological studies and drug efficacy trials of P. falciparum.

Methods

Targeted sequencing of molecular marker csp and novel marker cpmp was conducted in duplicate on mixtures of parasite culture strains and 37 field samples. A protocol allowing to multiplex up to 384 samples in a single sequencing run was applied. Software “HaplotypR” was developed for data analysis.

Results

Cpmp was highly diverse (H_e = 0.96) in contrast to csp (H_e = 0.57). Minority clones were robustly detected if their frequency was >1%. False haplotype calls owing to sequencing errors were observed below that threshold.

Conclusions

To reliably detect haplotypes at very low frequencies, experiments are best performed in duplicate and should aim for coverage of >10′000 reads/amplicon. When compared to length polymorphic marker msp2, highly multiplexed amplicon sequencing displayed greater sensitivity in detecting minority clones.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4260-y) contains supplementary material, which is available to authorized users.

Malaria Epidemiology at the Clone Level

Genotyping to distinguish between parasite clones is nowadays a standard in many molecular epidemiological studies of malaria. It has become crucial in drug trials and to follow individual clones in epidemiological studies, and to understand how drug resistance emerges and spreads. Here, we review the applications of the increasingly available genotyping tools and whole-genome sequencing data, and argue for a better integration of population genetics findings into malaria-control strategies.

Malaria Genomics in the Era of Eradication

The first reference genome assembly for the Plasmodium falciparum malaria parasite was completed over a decade ago, and the impact of this and other genomic resources on malaria research has been significant. Genomic resources for other malaria parasites are being established, even as P. falciparum continues to be the focus of development of new genomic methods and applications. Here we review the impact and applications of genomic data on malaria research, and discuss future needs and directions as genomic data generation becomes less expensive and more decentralized. Specifically, we focus on how population genomic strategies can be utilized to advance the malaria eradication agenda.

What is a vector?

Many important and rapidly emerging pathogens of humans, livestock and wildlife are 'vector-borne'. However, the term 'vector' has been applied to diverse agents in a broad range of epidemiological systems. In this perspective, we briefly review some common definitions, identify the strengths and weaknesses of each and consider the functional differences between vectors and other hosts from a range of ecological, evolutionary and public health perspectives. We then consider how the use of designations can afford insights into our understanding of epidemiological and evolutionary processes that are not otherwise apparent. We conclude that from a medical and veterinary perspective, a combination of the 'haematophagous arthropod' and 'mobility' definitions is most useful because it offers important insights into contact structure and control and emphasizes the opportunities for pathogen shifts among taxonomically similar species with similar feeding modes and internal environments. From a population dynamics and evolutionary perspective, we suggest that a combination of the 'micropredator' and 'sequential' definition is most appropriate because it captures the key aspects of transmission biology and fitness consequences for the pathogen and vector itself. However, we explicitly recognize that the value of a definition always depends on the research question under study.This article is part of the themed issue 'Opening the black box: re-examining the ecology and evolution of parasite transmission'.

Selective sweep suggests transcriptional regulation may underlie Plasmodium vivax resilience to malaria control measures in Cambodia

Cambodia, in which both Plasmodium vivax and Plasmodium falciparum are endemic, has been the focus of numerous malaria-control interventions, resulting in a marked decline in overall malaria incidence. Despite this decline, the number of P vivax cases has actually increased. To understand better the factors underlying this resilience, we compared the genetic responses of the two species to recent selective pressures. We sequenced and studied the genomes of 70 P vivax and 80 P falciparum isolates collected between 2009 and 2013. We found that although P falciparum has undergone population fracturing, the coendemic P vivax population has grown undisrupted, resulting in a larger effective population size, no discernable population structure, and frequent multiclonal infections. Signatures of selection suggest recent, species-specific evolutionary differences. Particularly, in contrast to P falciparum, P vivax transcription factors, chromatin modifiers, and histone deacetylases have undergone strong directional selection, including a particularly strong selective sweep at an AP2 transcription factor. Together, our findings point to different population-level adaptive mechanisms used by P vivax and P falciparum parasites. Although population substructuring in P falciparum has resulted in clonal outgrowths of resistant parasites, P vivax may use a nuanced transcriptional regulatory approach to population maintenance, enabling it to preserve a larger, more diverse population better suited to facing selective threats. We conclude that transcriptional control may underlie P vivax's resilience to malaria control measures. Novel strategies to target such processes are likely required to eradicate P vivax and achieve malaria elimination.

Dissecting malaria biology and epidemiology using population genetics and genomics

Molecular approaches have an increasingly recognized utility in surveillance of malaria parasite populations, not only in defining prevalence and incidence with higher sensitivity than traditional methods, but also in monitoring local and regional parasite transmission patterns. In this review, we provide an overview of population genetic and genomic studies of human-infecting Plasmodium species, highlighting recent advances in the field. In accordance with the renewed impetus for malaria eradication, many studies are now using genetic and genomic epidemiology to support local evidence-based intervention strategies. Microsatellite genotyping remains a popular approach for both Plasmodium falciparum and Plasmodium vivax. However, with the increasing availability of whole genome sequencing data enabling effective single nucleotide polymorphism-based panels tailored to a given study question and setting, this approach is gaining popularity. The availability of new reference genomes for Plasmodium malariae and Plasmodium ovale should see a surge in similar molecular studies on these currently neglected species. Genomic studies are revealing new insights into important adaptive mechanisms of the parasite including antimalarial drug resistance. The advent of new methodologies such as selective whole genome amplification for dealing with extensive human DNA in low density field isolates should see genome-wide approaches becoming routine for parasite surveillance once the economic costs outweigh the current cost benefits of targeted approaches.