The perils of PCR: can we accurately 'correct' antimalarial trials?

Genotyping methods to distinguish Plasmodium falciparum recrudescence from new infection for the assessment of antimalarial drug efficacy: an observational, single-centre, comparison study

BACKGROUND

Distinguishing Plasmodium falciparum recrudescence from new infections is crucial for the assessment of antimalarial drug efficacy against P falciparum. We aimed to compare the efficacy of different genotyping methods to assess their effect on drug efficacy estimates, particularly in patients from high-transmission settings with polyclonal infections.

METHODS

In this head-to-head comparison study, we compared five different genotyping methods currently used: fast capillary electrophoresis (F-CE) using msp1, msp2, and glurp; high-resolution capillary electrophoresis (H-CE) using msp1, msp2, and glurp; H-CE using microsatellites; targeted amplicon deep sequencing (TADS) using single nucleotide polymorphism (SNP)-rich markers; and high-resolution melting (HRM) analysis using msp1 and msp2. We assessed their sensitivity in detecting minority clones in polyclonal infections, their reproducibility, and the genetic diversity of the markers used. Our study used four well characterised P falciparum laboratory strains mixed in varying ratios, and 20 paired samples collected from an in-vivo clinical trial. The experiments were performed at the Swiss Tropical and Public Health Institute in Basel, Switzerland between May 5, 2020, and Aug 23, 2021.

FINDINGS

H-CE using msp1 and msp2 and TADS revealed the highest sensitivity in detecting minority clones (up to ratios of 1:100 for H-CE and 50:1:1:1 for TADS in the FCB1:HB3 and 3D7:K1:HB3:FCB1 laboratory strain mixtures, respectively), highest reproducibility (intra-assay: 99% and 91% for H-CE and TADS, respectively; inter-assay: 98% and 92% for H-CE and TADS, respectively), and highest genetic diversity in the used markers (up to 36 and 32 unique genotypes in 20 paired samples for H-CE using msp2 and TADS using cpmp, respectively). Microsatellites assessed by H-CE had a lower genetic diversity compared with msp1, msp2, and glurp assessed by H-CE and the SNP-rich markers assessed by TADS, with a maximum of 13 unique genotypes, and some genotypes having allelic frequencies larger than 30%. Markers used by TADS gave the most consistent results in distinguishing recrudescence from new infection across all methods (in 18 of 20 pairs of samples vs 15 of 20 pairs for H-CE).

INTERPRETATION

WHO currently recommends replacing glurp with microsatellites. However, in this study, the replacement of glurp with microsatellites did not change the genotyping outcome, probably due to the lower genetic diversity of microsatellites. More studies with large sample sizes are required to identify the most suitable microsatellites that could replace glurp. Our study indicates that TADS should be considered the gold standard for genotyping to distinguish recrudescence from new infection, and that it should be used to validate other methods.

FUNDING

Swiss Tropical and Public Health Institute.

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.

msp1, msp2, and glurp genotyping to differentiate Plasmodium falciparum recrudescence from reinfections during prevention of reestablishment phase, Sri Lanka, 2014-2019

BACKGROUND

Sri Lanka after eliminating malaria in 2012, is in the prevention of re-establishment (POR) phase. Being a tropical country with high malariogenic potential, maintaining vigilance is important. All malaria cases are investigated epidemiologically and followed up by integrated drug efficacy surveillance (iDES). Occasionally, that alone is not adequate to differentiate Plasmodium falciparum reinfections from recrudescences. This study evaluated the World Health Organization and Medicines for Malaria Venture (MMV) recommended genotyping protocol for the merozoite surface proteins (msp1, msp2) and the glutamate-rich protein (glurp) to discriminate P. falciparum recrudescence from reinfection in POR phase.

METHODS

All P. falciparum patients detected from April 2014 to December 2019 were included in this study. Patients were treated and followed up by iDES up to 28 days and were advised to get tested if they develop fever at any time over the following year. Basic socio-demographic information including history of travel was obtained. Details of the malariogenic potential and reactive entomological and parasitological surveillance carried out by the Anti Malaria Campaign to exclude the possibility of local transmission were also collected. The msp1, msp2, and glurp genotyping was performed for initial and any recurrent infections. Classification of recurrent infections as recrudescence or reinfection was done based on epidemiological findings and was compared with the genotyping outcome.

RESULTS

Among 106 P. falciparum patients, six had recurrent infections. All the initial infections were imported, with a history of travel to malaria endemic countries. In all instances, the reactive entomological and parasitological surveillance had no evidence for local transmission. Five recurrences occurred within 28 days of follow-up and were classified as recrudescence. They have not travelled to malaria endemic countries between the initial and recurrent infections. The other had a recurrent infection after 105 days. It was assumed a reinfection, as he had travelled to the same malaria endemic country in between the two malaria attacks. Genotyping confirmed the recrudescence and the reinfection.

CONCLUSIONS

The msp1, msp2 and glurp genotyping method accurately differentiated reinfections from recrudescence. Since reinfection without a history of travel to a malaria endemic country would mean local transmission, combining genotyping outcome with epidemiological findings will assist classifying malaria cases without any ambiguity.

The temporal dynamics of Plasmodium species infection after artemisinin-based combination therapy (ACT) among asymptomatic children in the Hohoe municipality, Ghana

BACKGROUND

The routine surveillance of asymptomatic malaria using nucleic acid-based amplification tests is essential in obtaining reliable data that would inform malaria policy formulation and the implementation of appropriate control measures.

METHODS

In this study, the prevalence rate and the dynamics of Plasmodium species among asymptomatic children (n = 1697) under 5 years from 30 communities within the Hohoe municipality in Ghana were determined.

RESULTS AND DISCUSSION

The observed prevalence of Plasmodium parasite infection by polymerase chain reaction (PCR) was 33.6% (571/1697), which was significantly higher compared to that obtained by microscopy [26.6% (451/1697)] (P < 0.0001). Based on species-specific analysis by nested PCR, Plasmodium falciparum infection [33.6% (570/1697)] was dominant, with Plasmodium malariae, Plasmodium ovale and Plasmodium vivax infections accounting for 0.1% (1/1697), 0.0% (0/1697), and 0.0% (0/1697), respectively. The prevalence of P. falciparum infection among the 30 communities ranged from 0.0 to 82.5%. Following artesunate-amodiaquine (AS + AQ, 25 mg/kg) treatment of a sub-population of the participants (n = 184), there was a substantial reduction in Plasmodium parasite prevalence by 100% and 79.2% on day 7 based on microscopy and nested PCR analysis, respectively. However, there was an increase in parasite prevalence from day 14 to day 42, with a subsequent decline on day 70 by both microscopy and nested PCR. For parasite clearance rate analysis, we found a significant proportion of the participants harbouring residual Plasmodium parasites or parasite genomic DNA on day 1 [65.0% (13/20)], day 2 [65.0% (13/20)] and day 3 [60.0% (12/20)] after initiating treatment. Of note, gametocyte carriage among participants was low before and after treatment.

CONCLUSION

Taken together, the results indicate that a significant number of individuals could harbour residual Plasmodium parasites or parasite genomic DNA after treatment. The study demonstrates the importance of routine surveillance of asymptomatic malaria using sensitive nucleic acid-based amplification techniques.

WHO antimalarial trial guidelines: good science, bad news?

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

Persistent Submicroscopic Plasmodium falciparum Parasitemia 72 Hours after Treatment with Artemether-Lumefantrine Predicts 42-Day Treatment Failure in Mali and Burkina Faso

A recent randomized controlled trial, the WANECAM (West African Network for Clinical Trials of Antimalarial Drugs) trial, conducted at seven centers in West Africa, found that artemether-lumefantrine, artesunate-amodiaquine, pyronaridine-artesunate, and dihydroartemisinin-piperaquine all displayed good efficacy. However, artemether-lumefantrine was associated with a shorter interval between clinical episodes than the other regimens. In a further comparison of these therapies, we identified cases of persisting submicroscopic parasitemia by quantitative PCR (qPCR) at 72 h posttreatment among WANECAM participants from 5 sites in Mali and Burkina Faso, and we compared treatment outcomes for this group to those with complete parasite clearance by 72 h. Among 552 evaluable patients, 17.7% had qPCR-detectable parasitemia at 72 h during their first treatment episode. This proportion varied among sites, reflecting differences in malaria transmission intensity, but did not differ among pooled drug treatment groups. However, patients who received artemether-lumefantrine and were qPCR positive at 72 h were significantly more likely to have microscopically detectable recurrent Plasmodium falciparum parasitemia by day 42 than those receiving other regimens and experienced, on average, a shorter interval before the next clinical episode. Haplotypes of pfcrt and pfmdr1 were also evaluated in persisting parasites. These data identify a possible threat to the parasitological efficacy of artemether-lumefantrine in West Africa, over a decade since it was first introduced on a large scale.

Increased Sensitivity of Plasmodium falciparum to Artesunate/Amodiaquine Despite 14 Years as First-Line Malaria Treatment, Zanzibar

Artemisinin-based combination therapies (ACTs) are first-line treatments for uncomplicated Plasmodium falciparum malaria. ACT resistance is spreading in Asia but not yet in Africa. Reduced effects of ACT partner drugs have been reported but with little information regarding widely used artesunate/amodiaquine (ASAQ). We studied its efficacy in Zanzibar after 14 years as first-line treatment directly by an in vivo, single-armed trial and indirectly by prevalences of different genotypes in the P. falciparum chloroquine-resistance transporter, multidrug-resistance 1, and Kelch 13 propeller domain genes. In vivo efficacy was higher during 2017 (100%; 95% CI 97.4%-100%) than during 2002-2005 (94.7%; 95% CI 91.9%-96.7%) (p = 0.003). Molecular findings showed no artemisinin resistance-associated genotypes and major increases in genotypes associated with high sensitivity/efficacy for amodiaquine than before ASAQ was introduced. Thus, the efficacy of ASAQ is maintained and appears to be increased after long-term use in contrast to what is observed for other ACTs used in Africa.

Variation in Calculating and Reporting Antimalarial Efficacy against Plasmodium falciparum in Sub-Saharan Africa: A Systematic Review of Published Reports

Antimalarials, in particular artemisinin-based combination therapies (ACTs), are critical tools in reducing the global burden of malaria, which is concentrated in sub-Saharan Africa. Performing and reporting antimalarial efficacy studies in a transparent and standardized fashion permit comparison of efficacy outcomes across countries and time periods. This systematic review summarizes study compliance with WHO laboratory and reporting guidance pertaining to antimalarial therapeutic efficacy studies and evaluates how well studies from sub-Saharan Africa adhered to these guidelines. We included all published studies (January 2020 or before) performed in sub-Saharan Africa where ACT efficacy for treatment of uncomplicated Plasmodium falciparum infection was reported. The primary outcome was a composite indicator for study methodology consistent with WHO guidelines for statistical analysis of corrected efficacy, defined as an article presenting a Kaplan-Meier survival analysis of corrected efficacy or reporting a per-protocol analysis where new infections were excluded from the numerator and denominator. Of 581 articles screened, we identified 279 for the review. Molecular correction was used in 83% (232/279) to distinguish new infections from recrudescences in subjects experiencing recurrent parasitemia. Only 45% (99/221) of articles with therapeutic efficacy as a primary outcome and performing molecular correction reported corrected efficacy outcomes calculated in a way consistent with WHO recommendations. These results indicate a widespread lack of compliance with WHO-recommended methods of analysis, which may result in biases in how antimalarial effectiveness is being measured and reported from sub-Saharan Africa.

Nonparametric Binary Classification to Distinguish Closely Related versus Unrelated Plasmodium falciparum Parasites

Assessing genetic relatedness of Plasmodium falciparum genotypes is a key component of antimalarial efficacy trials. Previous methods have focused on determining a priori definitions of the level of genetic similarity sufficient to classify two infections as sharing the same strain. However, factors such as mixed-strain infections, allelic suppression, imprecise typing methods, and heterozygosity complicate comparisons of apicomplexan genotypes. Here, we introduce a novel method for nonparametric statistical testing of relatedness for P. falciparum parasites. First, the background distribution of genetic distance between unrelated strains is computed. Second, a threshold genetic distance is computed from this empiric distribution of distances to demarcate genetic distances that are unlikely to have arisen by chance. Third, the genetic distance between paired samples is computed, and paired samples with genetic distances below the threshold are classified as related. The method is designed to work with any arbitrary genetic distance definition. We validated this procedure using two independent approaches to calculating genetic distance. We assessed the concordance of the novel nonparametric classification with a gold-standard Bayesian approach for 175 pairs of recurrent P. falciparum episodes from previously published malaria efficacy trials with microsatellite data from five studies in Guinea and Angola. The novel nonparametric approach was 98% sensitive and 84-89% specific in correctly identifying related genotypes compared with a gold-standard Bayesian algorithm. The approach provides a unified and systematic method to statistically assess relatedness of P. falciparum parasites using arbitrary genetic distance methodologies.

Adherence to Mass Drug Administration with Dihydroartemisinin-Piperaquine and Plasmodium falciparum Clearance in Southern Province, Zambia

Mass drug administration (MDA) with artemisinin combination therapy is a potentially useful tool for malaria elimination programs, but its success depends partly on drug effectiveness and treatment coverage in the targeted population. As part of a cluster-randomized controlled trial in Southern Province, Zambia evaluating the impact of MDA and household focal MDA (fMDA) with dihydroartemisinin-piperaquine (DHAp), sub-studies were conducted investigating population drug adherence rates and effectiveness of DHAp as administered in clearing Plasmodium falciparum infections following household mass administration. Adherence information was reported for 181,534 of 336,821 DHAp (53.9%) treatments administered during four rounds of MDA/fMDA, of which 153,197 (84.4%) reported completing the full course of DHAp. The proportion of participants fully adhering to the treatment regimen differed by MDA modality (MDA versus fMDA), RDT status, and whether the first dose was observed by those administering treatments. Among a subset of participants receiving DHAp and selected for longitudinal follow-up, 58 were positive for asexual-stage P. falciparum infection by microscopy at baseline. None of the 45 participants followed up at days 3 and/or 7 were slide positive for asexual-stage parasitemia. For those with longer term follow-up, one participant was positive 47 days after treatment, and two additional participants were positive after 69 days, although these two were determined to be new infections by genotyping. High completion of a 3-day course of DHAp and parasite clearance in the context of household MDA are promising as Zambia's National Malaria Programme continues to weigh appropriate interventions for malaria elimination.

A Computer Modelling Approach To Evaluate the Accuracy of Microsatellite Markers for Classification of Recurrent Infections during Routine Monitoring of Antimalarial Drug Efficacy

Antimalarial drugs have long half-lives, so clinical trials to monitor their efficacy require long periods of follow-up to capture drug failure that may become patent only weeks after treatment. Reinfections often occur during follow-up, so robust methods of distinguishing drug failures (recrudescence) from emerging new infections are needed to produce accurate failure rate estimates. Molecular correction aims to achieve this by comparing the genotype of a patient's pretreatment (initial) blood sample with that of any infection that occurs during follow-up, with matching genotypes indicating drug failure. We use an in silico approach to show that the widely used match-counting method of molecular correction with microsatellite markers is likely to be highly unreliable and may lead to gross under- or overestimates of the true failure rates, depending on the choice of matching criterion. A Bayesian algorithm for molecular correction was previously developed and utilized for analysis of in vivo efficacy trials. We validated this algorithm using in silico data and showed it had high specificity and generated accurate failure rate estimates. This conclusion was robust for multiple drugs, different levels of drug failure rates, different levels of transmission intensity in the study sites, and microsatellite genetic diversity. The Bayesian algorithm was inherently unable to accurately identify low-density recrudescence that occurred in a small number of patients, but this did not appear to compromise its utility as a highly effective molecular correction method for analyzing microsatellite genotypes. Strong consideration should be given to using Bayesian methodology to obtain accurate failure rate estimates during routine monitoring trials of antimalarial efficacy that use microsatellite markers.

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.

Dealing with indeterminate outcomes in antimalarial drug efficacy trials: a comparison between complete case analysis, multiple imputation and inverse probability weighting

BACKGROUND

Antimalarial clinical efficacy studies for uncomplicated Plasmodium falciparum malaria frequently encounter situations in which molecular genotyping is unable to discriminate between parasitic recurrence, either new infection or recrudescence. The current WHO guideline recommends excluding these individuals with indeterminate outcomes in a complete case (CC) analysis. Data from the four artemisinin-based combination (4ABC) trial was used to compare the performance of multiple imputation (MI) and inverse probability weighting (IPW) against the standard CC analysis for dealing with indeterminate recurrences.

METHODS

3369 study participants from the multicentre study (4ABC trial) with molecularly defined parasitic recurrence treated with three artemisinin-based combination therapies were used to represent a complete dataset. A set proportion of recurrent infections (10, 30 and 45%) were reclassified as missing using two mechanisms: a completely random selection (mechanism 1); missingness weakly dependent (mechanism 2a) and strongly dependent (mechanism 2b) on treatment and transmission intensity. The performance of MI, IPW and CC approaches in estimating the Kaplan-Meier (K-M) probability of parasitic recrudescence at day 28 was then compared. In addition, the maximum likelihood estimate of the cured proportion was presented for further comparison (analytical solution). Performance measures (bias, relative bias, standard error and coverage) were reported as an average from 1000 simulation runs.

RESULTS

The CC analyses resulted in absolute underestimation of K-M probability of day 28 recrudescence by up to 1.7% and were associated with reduced precision and poor coverage across all the scenarios studied. Both MI and IPW method performed better (greater consistency and greater efficiency) compared to CC analysis. In the absence of censoring, the analytical solution provided the most consistent and accurate estimate of cured proportion compared to the CC analyses.

CONCLUSIONS

The widely used CC approach underestimates antimalarial failure; IPW and MI procedures provided efficient and consistent estimates and should be considered when reporting the results of antimalarial clinical trials, especially in areas of high transmission, where the proportion of indeterminate outcomes could be large. The analytical solution estimating the cured proportion could provide an alternative approach, in scenarios with minimal censoring due to loss to follow-up or new infections.

Improving Methods for Analyzing Antimalarial Drug Efficacy Trials: Molecular Correction Based on Length-Polymorphic Markers msp-1, msp-2, and glurp

Drug efficacy trials monitor the continued efficacy of front-line drugs against falciparum malaria. Overestimating efficacy results in a country retaining a failing drug as first-line treatment with associated increases in morbidity and mortality, while underestimating drug effectiveness leads to removal of an effective treatment with substantial practical and economic implications.

Drug efficacy trials monitor the continued efficacy of front-line drugs against falciparum malaria. Overestimating efficacy results in a country retaining a failing drug as first-line treatment with associated increases in morbidity and mortality, while underestimating drug effectiveness leads to removal of an effective treatment with substantial practical and economic implications. Trials are challenging: they require long durations of follow-up to detect drug failures, and patients are frequently reinfected during that period. Molecular correction based on parasite genotypes distinguishes reinfections from drug failures to ensure the accuracy of failure rate estimates. Several molecular correction “algorithms” have been proposed, but which is most accurate and/or robust remains unknown. We used pharmacological modeling to simulate parasite dynamics and genetic signals that occur in patients enrolled in malaria drug clinical trials. We compared estimates of treatment failure obtained from a selection of proposed molecular correction algorithms against the known “true” failure rate in the model. Our findings are as follows. (i) Molecular correction is essential to avoid substantial overestimates of drug failure rates. (ii) The current WHO-recommended algorithm consistently underestimates the true failure rate. (iii) Newly proposed algorithms produce more accurate failure rate estimates; the most accurate algorithm depends on the choice of drug, trial follow-up length, and transmission intensity. (iv) Long durations of patient follow-up may be counterproductive; large numbers of new infections accumulate and may be misclassified, overestimating drug failure rate. (v) Our model was highly consistent with existing in vivo data. The current WHO-recommended method for molecular correction and analysis of clinical trials should be reevaluated and updated.

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.

Drug-Resistance and Population Structure of Plasmodium falciparum Across the Democratic Republic of Congo Using High-Throughput Molecular Inversion Probes

A better understanding of the drivers of the spread of malaria parasites and drug resistance across space and time is needed. These drivers can be elucidated using genetic tools. Here, a novel molecular inversion probe (MIP) panel targeting all major drug-resistance mutations and a set of microsatellites was used to genotype Plasmodium falciparum infections of 552 children from the 2013-2014 Demographic and Health Survey conducted in the Democratic Republic of the Congo (DRC). Microsatellite-based analysis of population structure suggests that parasites within the DRC form a homogeneous population. In contrast, sulfadoxine-resistance markers in dihydropteroate synthase show marked spatial structure with ongoing spread of double and triple mutants compared with 2007. These findings suggest that parasites in the DRC remain panmictic despite rapidly spreading antimalarial-resistance mutations. Moreover, highly multiplexed targeted sequencing using MIPs emerges as a cost-effective method for elucidating pathogen genetics in complex infections in large cohorts.

Identifying Recrudescent Plasmodium falciparum in Treated Malaria Patients by Real-time PCR and High Resolution Melt Analysis of Genetic Diversity

Recurrent parasitaemia during follow up of clinical trials of antimalarial drug efficacy results from either recrudescence of parasites surviving treatment or from parasites newly emerging from the hepatic stage of infection. Nested PCR is used to distinguish these two possibilities and the technique is difficult to standardise. There is risk of both false positive and false negative results, leading to misclassification errors. The high-resolution melt (HRM) assay was developed with pairs of conserved primers targeting blocks of merozoite surface protein 1 and 2 (msp1 and msp2) genes, and polymorphisms were compared using sequence-confirmed Plasmodium falciparum DNA samples from laboratory isolates. In this study, the HRM dissociation profiles of msp1 and msp2 amplicons were determined and validated against parasite isolates from malaria patients. The msp1 and msp2 profiles of both laboratory and clinical isolates were reproducibly differentiated by HRM. These rapid assays are performed in a closed-tube system, and so avoid cross-contamination while increasing throughput, which are two major advantages. The HRM assays offer significant gains in simplicity, speed and interpretation of results, and reduced analysis cost, for studies that require discrimination of parasite clones. Assay performance in large-scale studies utilizing DNA samples derived from filter-paper bloodspots should now be evaluated.

Comparison of Two Genotyping Methods for Distinguishing Recrudescence from Reinfection in Antimalarial Drug Efficacy/Effectiveness Trials

Genotyping of allelic variants of Plasmodium falciparum merozoite surface proteins 1 and 2 (msp-1 and msp-2), and the glutamate-rich protein is the gold standard for distinguishing reinfections from recrudescences in antimalarial drug trials. We compared performance of the recently developed 24-single-nucleotide polymorphism (SNP) Barcoding Assay against msp-1 and msp-2 genotyping in a cluster-randomized effectiveness trial of artemether-lumefantrine and dihydroartemisinin-piperaquine in Malawi. Rates of recrudescence and reinfection estimated by the two methods did not differ significantly (Fisher's exact test; P = 0.887 and P = 0.768, respectively). There was a strong agreement between the two methods in predicting treatment outcomes and resolving the genetic complexity of malaria infections in this setting. These results support the use of this SNP assay as an alternative method for correcting antimalarial efficacy/effectiveness data.

Putting evolution in elimination: Winning our ongoing battle with evolving malaria mosquitoes and parasites

Since 2000, the world has made significant progress in reducing malaria morbidity and mortality, and several countries in Africa, South America and South-East Asia are working hard to eliminate the disease. These elimination efforts continue to rely heavily on antimalarial drugs and insecticide-based interventions, which remain the cornerstones of malaria treatment and prevention. However, resistance has emerged against nearly every antimalarial drug and insecticide that is available. In this review we discuss the evolutionary consequences of the way we currently implement antimalarial interventions, which is leading to resistance and may ultimately lead to control failure, but also how evolutionary principles can be applied to extend the lifespan of current and novel interventions. A greater understanding of the general evolutionary principles that are at the core of emerging resistance is urgently needed if we are to develop improved resistance management strategies with the ultimate goal to achieve a malaria-free world.

Species and genotype diversity of Plasmodium in malaria patients from Gabon analysed by next generation sequencing

Background

Six Plasmodium species are known to naturally infect humans. Mixed species infections occur regularly but morphological discrimination by microscopy is difficult and multiplicity of infection (MOI) can only be evaluated by molecular methods. This study investigated the complexity of Plasmodium infections in patients treated for microscopically detected non-falciparum or mixed species malaria in Gabon.

Methods

Ultra-deep sequencing of nucleus (18S rRNA), mitochondrion, and apicoplast encoded genes was used to evaluate Plasmodium species diversity and MOI in 46 symptomatic Gabonese patients with microscopically diagnosed non-falciparum or mixed species malaria.

Results

Deep sequencing revealed a large complexity of confections in patients with uncomplicated malaria, both on species and genotype levels. Mixed infections involved up to four parasite species (Plasmodium falciparum, Plasmodium malariae, Plasmodium ovale curtisi, and P. ovale wallikeri). Multiple genotypes from each species were determined from the asexual 18S rRNA gene. 17 of 46 samples (37%) harboured multiple genotypes of at least one Plasmodium species. The number of genotypes per sample (MOI) was highest in P. malariae (n = 4), followed by P. ovale curtisi (n = 3), P. ovale wallikeri (n = 3), and P. falciparum (n = 2). The highest combined genotype complexity in samples that contained mixed-species infections was seven.

Conclusions

Ultra-deep sequencing showed an unexpected breadth of Plasmodium species and within species diversity in clinical samples. MOI of P. ovale curtisi, P. ovale wallikeri and P. malariae infections were higher than anticipated and contribute significantly to the burden of malaria in Gabon.

Electronic supplementary material

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

Genetic Diversity of Plasmodium falciparum Populations in Malaria Declining Areas of Sabah, East Malaysia

Malaysia has a national goal to eliminate malaria by 2020. Understanding the genetic diversity of malaria parasites in residual transmission foci can provide invaluable information which may inform the intervention strategies used to reach elimination targets. This study was conducted to determine the genetic diversity level of P. falciparum isolates in malaria residual foci areas of Sabah. Malaria active case detection was conducted in Kalabakan and Kota Marudu. All individuals in the study sites were screened for malaria infection by rapid diagnostic test. Blood from P. falciparum-infected individuals were collected on filter paper prior to DNA extraction. Genotyping was performed using merozoite surface protein-1 (MSP-1), merozoite surface protein-2 (MSP-2), glutamate rich protein (GLURP) and 10 neutral microsatellite loci markers. The size of alleles, multiplicity of infection (MOI), mean number of alleles (Na), expected heterozygosity (He), linkage disequilibrium (LD) and genetic differentiation (FST) were determined. In Kalabakan, the MSP-1 and MSP-2 alleles were predominantly K1 and FC27 family types, respectively. The GLURP genotype VI (751-800 bp) was predominant. The MOI for MSP-1 and MSP-2 were 1.65 and 1.20, respectively. The Na per microsatellite locus was 1.70. The He values for MSP-1, MSP-2, GLURP and neutral microsatellites were 0.17, 0.37, 0.70 and 0.33, respectively. In Kota Marudu, the MSP-1 and MSP-2 alleles were predominantly MAD20 and 3D7 family types, respectively. The GLURP genotype IV (651-700 bp) was predominant. The MOI for both MSP-1 and MSP-2 was 1.05. The Na per microsatellite locus was 3.60. The He values for MSP-1, MSP-2, GLURP and neutral microsatellites were 0.24, 0.25, 0.69 and 0.30, respectively. A significant LD was observed in Kalabakan (0.495, p<0.01) and Kota Marudu P. falciparum populations (0.601, p<0.01). High genetic differentiation between Kalabakan and Kota Marudu P. falciparum populations was observed (FST = 0.532). The genetic data from the present study highlighted the limited diversity and contrasting genetic pattern of P. falciparum populations in the malaria declining areas of Sabah.

Cryptosporidium within-host genetic diversity: systematic bibliographical search and narrative overview

Knowledge of the within-host genetic diversity of a pathogen often has broad implications for disease management. Cryptosporidium protozoan parasites are among the most common causative agents of infectious diarrhoea. Current limitations of in vitro culture impose the use of uncultured isolates obtained directly from the hosts as operational units of Cryptosporidium genotyping. The validity of this practice is centred on the assumption of genetic homogeneity of the parasite within the host, and genetic studies often take little account of the within-host genetic diversity of Cryptosporidium. Yet, theory and experimental evidence contemplate genetic diversity of Cryptosporidium at the within-host scale, but this diversity is not easily identified by genotyping methods ill-suited for the resolution of DNA mixtures. We performed a systematic bibliographical search of the occurrence of within-host genetic diversity of Cryptosporidium parasites in epidemiological samples, between 2005 and 2015. Our results indicate that genetic diversity at the within-host scale, in the form of mixed species or intra-species diversity, has been identified in a large number (n=55) of epidemiological surveys of cryptosporidiosis in variable proportions, but has often been treated as a secondary finding and not analysed. As in malaria, there are indications that the scale of this diversity varies between geographical regions, perhaps depending on the prevailing transmission pathways. These results provide a significant knowledge base from which to draw alternative population genetic structure models, some of which are discussed in this paper.

Robust Algorithm for Systematic Classification of Malaria Late Treatment Failures as Recrudescence or Reinfection Using Microsatellite Genotyping

Routine therapeutic efficacy monitoring to measure the response to antimalarial treatment is a cornerstone of malaria control. To correctly measure drug efficacy, therapeutic efficacy studies require genotyping parasites from late treatment failures to differentiate between recrudescent infections and reinfections. However, there is a lack of statistical methods to systematically classify late treatment failures from genotyping data. A Bayesian algorithm was developed to estimate the posterior probability of late treatment failure being the result of a recrudescent infection from microsatellite genotyping data. The algorithm was implemented using a Monte Carlo Markov chain approach and was used to classify late treatment failures using published microsatellite data from therapeutic efficacy studies in Ethiopia and Angola. The algorithm classified 85% of the Ethiopian and 95% of the Angolan late treatment failures as either likely reinfection or likely recrudescence, defined as a posterior probability of recrudescence of <0.1 or >0.9, respectively. The adjusted efficacies calculated using the new algorithm differed from efficacies estimated using commonly used methods for differentiating recrudescence from reinfection. In a high-transmission setting such as Angola, as few as 15 samples needed to be genotyped in order to have enough power to correctly classify treatment failures. Analysis of microsatellite genotyping data for differentiating between recrudescence and reinfection benefits from an approach that both systematically classifies late treatment failures and estimates the uncertainty of these classifications. Researchers analyzing genotyping data from antimalarial therapeutic efficacy monitoring are urged to publish their raw genetic data and to estimate the uncertainty around their classification.

Plasmodium vivax Diversity and Population Structure across Four Continents

Plasmodium vivax is the geographically most widespread human malaria parasite. To analyze patterns of microsatellite diversity and population structure across countries of different transmission intensity, genotyping data from 11 microsatellite markers was either generated or compiled from 841 isolates from four continents collected in 1999–2008. Diversity was highest in South-East Asia (mean allelic richness 10.0–12.8), intermediate in the South Pacific (8.1–9.9) Madagascar and Sudan (7.9–8.4), and lowest in South America and Central Asia (5.5–7.2). A reduced panel of only 3 markers was sufficient to identify approx. 90% of all haplotypes in South Pacific, African and SE-Asian populations, but only 60–80% in Latin American populations, suggesting that typing of 2–6 markers, depending on the level of endemicity, is sufficient for epidemiological studies. Clustering analysis showed distinct clusters in Peru and Brazil, but little sub-structuring was observed within Africa, SE-Asia or the South Pacific. Isolates from Uzbekistan were exceptional, as a near-clonal parasite population was observed that was clearly separated from all other populations (F_ST>0.2). Outside Central Asia F_ST values were highest (0.11–0.16) between South American and all other populations, and lowest (0.04–0.07) between populations from South-East Asia and the South Pacific. These comparisons between P. vivax populations from four continents indicated that not only transmission intensity, but also geographical isolation affect diversity and population structure. However, the high effective population size results in slow changes of these parameters. This persistency must be taken into account when assessing the impact of control programs on the genetic structure of parasite populations.

Plasmodium vivax is the predominant malaria parasite in Latin America, Asia and the South Pacific. Different factors are expected to shape diversity and population structure across continents, e.g. transmission intensity which is much lower in South America as compared to Southeast-Asia and the South Pacific, or geographical isolation of P. vivax populations in the South Pacific. We have compiled data from 841 isolates from South and Central America, Africa, Central Asia, Southeast-Asia and the South Pacific typed with a panel of 11 microsatellite markers. Diversity was highest in Southeast-Asia, where transmission is intermediate-high and migration of infected hosts is high, and lowest in South America and Central Asia where malaria transmission is low and focal. Reducing the panel of microsatellites showed that 2–6 markers are sufficient for genotyping for most drug trials and epidemiological studies, as these markers can identify >90% of all haplotypes. Parasites clustered according to continental origin, with high population differentiation between South American and Central Asian populations and the other populations, and lowest differences between Southeast-Asia and the South Pacific. Current attempts to reduce malaria transmission might change this pattern, but only after transmission is reduced for an extended period of time.

Relevance of undetectably rare resistant malaria parasites in treatment failure: experimental evidence from Plasmodium chabaudi

Resistant malaria parasites are frequently found in mixed infections with drug-sensitive parasites. Particularly early in the evolutionary process, the frequency of these resistant mutants can be extremely low and below the level of molecular detection. We tested whether the rarity of resistance in infections impacted the health outcomes of treatment failure and the potential for onward transmission of resistance. Mixed infections of different ratios of resistant and susceptible Plasmodium chabaudi parasites were inoculated in laboratory mice and dynamics tracked during the course of infection using highly sensitive genotype-specific quantitative polymerase chain reaction (qPCR). Frequencies of resistant parasites ranged from 10% to 0.003% at the onset of treatment. We found that the rarer the resistant parasites were, the lower the likelihood of their onward transmission, but the worse the treatment failure was in terms of parasite numbers and disease severity. Strikingly, drug resistant parasites had the biggest impact on health outcomes when they were too rare to be detected by any molecular methods currently available for field samples. Indeed, in the field, these treatment failures would not even have been attributed to resistance.

Detection of persistent Plasmodium spp. infections in Ugandan children after artemether-lumefantrine treatment

During a longitudinal study investigating the dynamics of malaria in Ugandan lakeshore communities, a consistently high malaria prevalence was observed in young children despite regular treatment. To explore the short-term performance of artemether-lumefantrine (AL), a pilot investigation into parasite carriage after treatment(s) was conducted in Bukoba village. A total of 163 children (aged 2–7 years) with a positive blood film and rapid antigen test were treated with AL; only 8·7% of these had elevated axillary temperatures. On day 7 and then on day 17, 40 children (26·3%) and 33 (22·3%) were positive by microscopy, respectively. Real-time PCR analysis demonstrated that multi-species Plasmodium infections were common at baseline, with 41·1% of children positive for Plasmodium falciparum/Plasmodium malariae, 9·2% for P. falciparum/ Plasmodium ovale spp. and 8·0% for all three species. Moreover, on day 17, 39·9% of children infected with falciparum malaria at baseline were again positive for the same species, and 9·2% of those infected with P. malariae at baseline were positive for P. malariae. Here, chronic multi-species malaria infections persisted in children after AL treatment(s). Better point-of-care diagnostics for non-falciparum infections are needed, as well as further investigation of AL performance in asymptomatic individuals.

Residual Plasmodium falciparum parasitemia in Kenyan children after artemisinin-combination therapy is associated with increased transmission to mosquitoes and parasite recurrence

BACKGROUND

Parasite clearance time after artemisinin-based combination therapy (ACT) may be increasing in Asian and African settings. The association between parasite clearance following ACT and transmissibility is currently unknown.

METHODS

We determined parasite clearance dynamics by duplex quantitative polymerase chain reaction (qPCR) in samples collected in the first 3 days after treatment of uncomplicated malaria with ACT. Gametocyte carriage was determined by Pfs25 quantitative nucleic acid sequence-based amplification assays; infectiousness to mosquitoes by membrane-feeding assays on day 7 after treatment.

RESULTS

Residual parasitemia was detected by qPCR in 31.8% (95% confidence interval [CI], 24.6-39.8) of the children on day 3 after initiation of treatment. Residual parasitemia was associated with a 2-fold longer duration of gametocyte carriage (P = .0007), a higher likelihood of infecting mosquitoes (relative risk, 1.95; 95% CI, 1.17-3.24; P = .015), and a higher parasite burden in mosquitoes (incidence rate ratio, 2.92; 95% CI, 1.61-5.31; P < .001). Children with residual parasitemia were also significantly more likely to experience microscopically detectable parasitemia during follow-up (relative risk, 11.25; 95% CI, 4.08-31.01; P < .001).

CONCLUSIONS

Residual submicroscopic parasitemia is common after ACT and is associated with a higher transmission potential. Residual parasitemia may also have consequences for individual patients because of its higher risk of recurrent parasitemia.

Culture-adapted Plasmodium falciparum isolates from UK travellers: in vitro drug sensitivity, clonality and drug resistance markers

Background

The screening of lead compounds against in vitro parasite cultures is an essential step in the development of novel anti-malarial drugs, but currently relies on laboratory parasite lines established in vitro during the last century. This study sought to establish in continuous culture a series of recent Plasmodium falciparum isolates to represent the current parasite populations in Africa, all of which are now exposed to artemisinin combination therapy.

Methods

Pre-treatment P. falciparum isolates were obtained in EDTA, and placed into continuous culture after sampling of DNA. One post-treatment blood sample was also collected for each donor to monitor parasite clonality during clearance in vivo. IC₅₀ estimates were obtained for 11 anti-malarial compounds for each established parasite line, clonal multiplicity measured in vivo and in vitro, and polymorphic sites implicated in parasite sensitivity to drugs were investigated at the pfmdr1, pfcrt, pfdhfr, pfdhps and pfap2mu loci before and after treatment, and in the cultured lines.

Results

Plasmodium falciparum isolates from seven malaria patients with recent travel to three West African and two East African countries were successfully established in long-term culture. One of these, HL1211, was from a patient with recrudescent parasitaemia 14 days after a full course of artemether-lumefantrine. All established culture lines were shown to be polyclonal, reflecting the in vivo isolates from which they were derived, and at least two lines reliably produce gametocytes in vitro. Two lines displayed high chloroquine IC₅₀ estimates, and carried the CVIET haplotype at codons 72–76, whereas the remaining five lines carried the CVMNK haplotype and were sensitive in vitro. All were sensitive to the endoperoxides dihydroartemisinin and OZ277, but IC₅₀ estimates for lumefantrine varied, with the least sensitive parasites carrying pfmdr1 alleles encoding Asn at codon 86.

Conclusions

This study describes the establishment in continuous culture, in vitro drug sensitivity testing and molecular characterization of a series of multiclonal P. falciparum isolates taken directly from UK malaria patients following recent travel to various malaria-endemic countries in Africa. These “HL” isolates are available as an open resource for studies of drug response, antigenic diversity and other aspects of parasite biology.

Malaria haplotype frequency estimation

We present a Bayesian approach for estimating the relative frequencies of multi-single nucleotide polymorphism (SNP) haplotypes in populations of the malaria parasite Plasmodium falciparum by using microarray SNP data from human blood samples. Each sample comes from a malaria patient and contains one or several parasite clones that may genetically differ. Samples containing multiple parasite clones with different genetic markers pose a special challenge. The situation is comparable with a polyploid organism. The data from each blood sample indicates whether the parasites in the blood carry a mutant or a wildtype allele at various selected genomic positions. If both mutant and wildtype alleles are detected at a given position in a multiply infected sample, the data indicates the presence of both alleles, but the ratio is unknown. Thus, the data only partially reveals which specific combinations of genetic markers (i.e. haplotypes across the examined SNPs) occur in distinct parasite clones. In addition, SNP data may contain errors at non-negligible rates. We use a multinomial mixture model with partially missing observations to represent this data and a Markov chain Monte Carlo method to estimate the haplotype frequencies in a population. Our approach addresses both challenges, multiple infections and data errors.

Transfusion-transmitted malaria in Ghana

BACKGROUND

In sub-Saharan Africa, the prevalence of malaria parasitemia in blood donors varies from 0.6% to 50%. Although the burden of TTM in malaria-endemic countries is unknown, it is recommended that all donated blood is screened for malaria parasites. This study aimed to establish the incidence of TTM and identify a suitable screening test.

METHODS

Pregnant women, children, and immunocompromised malaria-negative transfusion recipients in a teaching hospital in Ghana were recruited over the course of 1 year. Parasites detected in recipients within 14 days of the transfusion were genotyped and compared to parasites in the transfused blood. The presence of genotypically identical parasites in the recipient and the transfused blood confirmed transfusion-transmitted malaria. Four malaria screening tests were compared to assess their usefulness in the context of African blood banks.

RESULTS

Of the 50 patients who received transfusions that were positive for Plasmodium falciparum by polymerase chain reaction (PCR), 7 recipients developed PCR-detectable parasitemia. In only 1 of the 50 recipients (2%) was the parasite identical to that in the transfused blood. The prevalence of P. falciparum malaria in transfused blood was 4.7% (21/445) by microscopy, 13.7% (60/440) by rapid diagnostic test, 18% (78/436) by PCR, and 22.2% (98/442) by enzyme immunoassay.

CONCLUSIONS

Although malaria parasites are commonly detected in blood donors in malaria-endemic areas, transfusion-transmitted malaria occurs infrequently. Policies recommend screening blood donors for malaria, but none of the commonly used methods is sufficiently sensitive to be used by blood banks in malaria-endemic countries.

Long term persistence of clonal malaria parasite Plasmodium falciparum lineages in the Colombian Pacific region

Background

Resistance to chloroquine and antifolate drugs has evolved independently in South America, suggesting that genotype - phenotype studies aimed at understanding the genetic basis of resistance to these and other drugs should be conducted in this continent. This research was conducted to better understand the population structure of Colombian Plasmodium falciparum in preparation for such studies.

Results

A set of 384 SNPs were genotyped in blood spot DNA samples from 447 P. falciparum infected subjects collected over a ten year period from four provinces of the Colombian Pacific coast to evaluate clonality, population structure and linkage disequilibrium (LD). Most infections (81%) contained a single predominant clone. These clustered into 136 multilocus genotypes (MLGs), with 32% of MLGs recovered from multiple (2 – 28) independent subjects. We observed extremely low genotypic richness (R = 0.42) and long persistence of MLGs through time (median = 537 days, range = 1 – 2,997 days). There was a high probability (>5%) of sampling parasites from the same MLG in different subjects within 28 days, suggesting caution is needed when using genotyping methods to assess treatment success in clinical drug trials. Panmixia was rejected as four well differentiated subpopulations (F_ST = 0.084 - 0.279) were identified. These occurred sympatrically but varied in frequency within the four provinces. Linkage disequilibrium (LD) decayed more rapidly (r² = 0.17 for markers <10 kb apart) than observed previously in South American samples.

Conclusions

We conclude that Colombian populations have several advantages for association studies, because multiple clone infections are uncommon and LD decays over the scale of one or a few genes. However, the extensive population structure and low genotype richness will need to be accounted for when designing and analyzing association studies.

Microsatellite analysis of malaria parasites

Microsatellites have been increasingly used to investigate the population structure of malaria parasites, to map genetic loci contributing to phenotypes such as drug resistance and virulence in laboratory crosses and genome-wide association studies and to distinguish between treatment failures and new infections in clinical trials. Here, we provide optimized protocols for genotyping highly polymorphic microsatellites sampled from across the genomes of the human malaria parasites Plasmodium falciparum and P. vivax that have been extensively used in research laboratories worldwide.

Estimating the numbers of malaria infections in blood samples using high-resolution genotyping data

People living in endemic areas often habour several malaria infections at once. High-resolution genotyping can distinguish between infections by detecting the presence of different alleles at a polymorphic locus. However the number of infections may not be accurately counted since parasites from multiple infections may carry the same allele. We use simulation to determine the circumstances under which the number of observed genotypes are likely to be substantially less than the number of infections present and investigate the performance of two methods for estimating the numbers of infections from high-resolution genotyping data. The simulations suggest that the problem is not substantial in most datasets: the disparity between the mean numbers of infections and of observed genotypes was small when there was 20 or more alleles, 20 or more blood samples, a mean number of infections of 6 or less and where the frequency of the most common allele was no greater than 20%. The issue of multiple infections carrying the same allele is unlikely to be a major component of the errors in PCR-based genotyping. Simulations also showed that, with heterogeneity in allele frequencies, the observed frequencies are not a good approximation of the true allele frequencies. The first method that we proposed to estimate the numbers of infections assumes that they are a good approximation and hence did poorly in the presence of heterogeneity. In contrast, the second method by Li et al estimates both the numbers of infections and the true allele frequencies simultaneously and produced accurate estimates of the mean number of infections.

Plasmodium falciparum population dynamics during the early phase of anti-malarial drug treatment in Tanzanian children with acute uncomplicated malaria

Background

This study aimed to explore Plasmodium falciparum population dynamics during the early phase of anti-malarial drug treatment with artemisinin-based combination therapy in children with clinical malaria in a high transmission area in Africa.

Methods

A total of 50 children aged 1-10 years with acute uncomplicated P. falciparum malaria in Bagamoyo District, Tanzania, were enrolled. Participants were hospitalized and received supervised standard treatment with artemether-lumefantrine according to body weight in six doses over 3 days. Blood samples were collected 11 times, i.e. at time of diagnosis (-2 h) and 0, 2, 4, 8, 16, 24, 36, 48, 60 and 72 h after initiation of treatment. Parasite population dynamics were assessed using nested polymerase chain reaction (PCR)-genotyping of merozoite surface protein (msp) 1 and 2.

Results

PCR-analyses from nine sequential blood samples collected after initiation of treatment identified 20 and 21 additional genotypes in 15/50 (30%) and 14/50 (28%) children with msp1 and msp2, respectively, non-detectable in the pre-treatment samples (-2 and 0 h combined). Some 15/20 (75%) and 14/21 (67%) of these genotypes were identified within 24 h, whereas 17/20 (85%) and 19/21 (90%) within 48 h for msp1 and msp2, respectively. The genotype profile was diverse, and varied considerably over time both within and between patients, molecular markers and their respective families.

Conclusion

PCR analyses from multiple blood samples collected during the early treatment phase revealed a complex picture of parasite sub-populations. This underlines the importance of interpreting PCR-outcomes with caution and suggests that the present use of PCR-adjustment from paired blood samples in anti-malarial drug trials may overestimate assessment of drug efficacy in high transmission areas in Africa.

The study is registered at http://www.clinicaltrials.gov with identifier NCT00336375.

Declining responsiveness of Plasmodium falciparum infections to artemisinin-based combination treatments on the Kenyan coast

Background

The emergence of artemisinin-resistant P. falciparum malaria in South-East Asia highlights the need for continued global surveillance of the efficacy of artemisinin-based combination therapies.

Methods

On the Kenyan coast we studied the treatment responses in 474 children 6–59 months old with uncomplicated P. falciparum malaria in a randomized controlled trial of dihydroartemisinin-piperaquine vs. artemether-lumefantrine from 2005 to 2008. (ISRCTN88705995)

Results

The proportion of patients with residual parasitemia on day 1 rose from 55% in 2005–2006 to 87% in 2007–2008 (odds ratio, 5.4, 95%CI, 2.7–11.1; P<0.001) and from 81% to 95% (OR, 4.1, 95%CI, 1.7–9.9; P = 0.002) in the DHA-PPQ and AM-LM groups, respectively. In parallel, Kaplan-Meier estimated risks of apparent recrudescent infection by day 84 increased from 7% to 14% (P = 0.1) and from 6% to 15% (P = 0.05) with DHA-PPQ and AM-LM, respectively. Coinciding with decreasing transmission in the study area, clinical tolerance to parasitemia (defined as absence of fever) declined between 2005–2006 and 2007–2008 (OR body temperature >37.5°C, 2.8, 1.9–4.1; P<0.001). Neither in vitro sensitivity of parasites to DHA nor levels of antibodies against parasite extract accounted for parasite clearance rates or changes thereof.

Conclusions

The significant, albeit small, decline through time of parasitological response rates to treatment with ACTs may be due to the emergence of parasites with reduced drug sensitivity, to the coincident reduction in population-level clinical immunity, or both. Maintaining the efficacy of artemisinin-based therapy in Africa would benefit from a better understanding of the mechanisms underlying reduced parasite clearance rates.

Trial Registration

Controlled-Trials.com ISRCTN88705995

Enhanced understanding of infectious diseases by fusing multiple datasets: a case study on malaria in the Western Brazilian Amazon region

BACKGROUND

A common challenge to the study of several infectious diseases consists in combining limited cross-sectional survey data, collected with a more sensitive detection method, with a more extensive (but biased) syndromic sentinel surveillance data, collected with a less sensitive method. Our article describes a novel modeling framework that overcomes this challenge, resulting in enhanced understanding of malaria in the Western Brazilian Amazon.

METHODOLOGY/PRINCIPAL FINDINGS

A cohort of 486 individuals was monitored using four cross-sectional surveys, where all participants were sampled regardless of symptoms (aggressive-active case detection), resulting in 1,383 microscopy and 1,400 polymerase chain reaction tests. Data on the same individuals were also obtained from the local surveillance facility (i.e., passive and active case detection), totaling 1,694 microscopy tests. Our model accommodates these multiple pathogen and case detection methods. This model is shown to outperform logistic regression in terms of interpretability of its parameters, ability to recover the true parameter values, and predictive performance. We reveal that the main infection determinant was the extent of forest, particularly during the rainy season and in close proximity to water bodies, and participation on forest activities. We find that time residing in Acrelandia (as a proxy for past malaria exposure) decreases infection risk but surprisingly increases the likelihood of reporting symptoms once infected, possibly because non-naïve settlers are only susceptible to more virulent Plasmodium strains. We suggest that the search for asymptomatic carriers should focus on those at greater risk of being infected but lower risk of reporting symptoms once infected.

CONCLUSIONS/SIGNIFICANCE

The modeling framework presented here combines cross-sectional survey data and syndromic sentinel surveillance data to shed light on several aspects of malaria that are critical for public health policy. This framework can be adapted to enhance inference on infectious diseases whenever asymptomatic carriers are important and multiple datasets are available.

Drug-resistant genotypes and multi-clonality in Plasmodium falciparum analysed by direct genome sequencing from peripheral blood of malaria patients

Naturally acquired blood-stage infections of the malaria parasite Plasmodium falciparum typically harbour multiple haploid clones. The apparent number of clones observed in any single infection depends on the diversity of the polymorphic markers used for the analysis, and the relative abundance of rare clones, which frequently fail to be detected among PCR products derived from numerically dominant clones. However, minority clones are of clinical interest as they may harbour genes conferring drug resistance, leading to enhanced survival after treatment and the possibility of subsequent therapeutic failure. We deployed new generation sequencing to derive genome data for five non-propagated parasite isolates taken directly from 4 different patients treated for clinical malaria in a UK hospital. Analysis of depth of coverage and length of sequence intervals between paired reads identified both previously described and novel gene deletions and amplifications. Full-length sequence data was extracted for 6 loci considered to be under selection by antimalarial drugs, and both known and previously unknown amino acid substitutions were identified. Full mitochondrial genomes were extracted from the sequencing data for each isolate, and these are compared against a panel of polymorphic sites derived from published or unpublished but publicly available data. Finally, genome-wide analysis of clone multiplicity was performed, and the number of infecting parasite clones estimated for each isolate. Each patient harboured at least 3 clones of P. falciparum by this analysis, consistent with results obtained with conventional PCR analysis of polymorphic merozoite antigen loci. We conclude that genome sequencing of peripheral blood P. falciparum taken directly from malaria patients provides high quality data useful for drug resistance studies, genomic structural analyses and population genetics, and also robustly represents clonal multiplicity.

How much remains undetected? Probability of molecular detection of human Plasmodia in the field

Background

In malaria endemic areas, most people are simultaneously infected with different parasite clones. Detection of individual clones is hampered when their densities fluctuate around the detection limit and, in case of P. falciparum, by sequestration during part of their life cycle. This has important implications for measures of levels of infection or for the outcome of clinical trials. This study aimed at measuring the detectability of individual P. falciparum and P. vivax parasite clones in consecutive samples of the same patient and at investigating the impact of sampling strategies on basic epidemiological measures such as multiplicity of infection (MOI).

Methods

Samples were obtained in a repeated cross-sectional field survey in 1 to 4.5 years old children from Papua New Guinea, who were followed up in 2-monthly intervals over 16 months. At each follow-up visit, two consecutive blood samples were collected from each child at intervals of 24 hours. Samples were genotyped for the polymorphic markers msp2 for P. falciparum and msp1F3 and MS16 for P. vivax. Observed prevalence and mean MOI estimated from single samples per host were compared to combined data from sampling twice within 24 h.

Findings and Conclusion

Estimated detectability was high in our data set (0.79 [95% CI 0.76–0.82] for P. falciparum and, depending on the marker, 0.61 [0.58–0.63] or 0.73 [0.71–0.75] for P. vivax). When genotyping data from sequential samples, collected 24 hours apart, were combined, the increase in measured prevalence was moderate, 6 to 9% of all infections were missed on a single day. The effect on observed MOI was more pronounced, 18 to 31% of all individual clones were not detected in a single bleed. Repeated sampling revealed little difference between detectability of P. falciparum and P. vivax.

Exposing malaria in-host diversity and estimating population diversity by capture-recapture using massively parallel pyrosequencing

Malaria infections commonly contain multiple genetically distinct variants. Mathematical and animal models suggest that interactions among these variants have a profound impact on the emergence of drug resistance. However, methods currently used for quantifying parasite diversity in individual infections are insensitive to low-abundance variants and are not quantitative for variant population sizes. To more completely describe the in-host complexity and ecology of malaria infections, we used massively parallel pyrosequencing to characterize malaria parasite diversity in the infections of a group of patients. By individually sequencing single strands of DNA in a complex mixture, this technique can quantify uncommon variants in mixed infections. The in-host diversity revealed by this method far exceeded that described by currently recommended genotyping methods, with as many as sixfold more variants per infection. In addition, in paired pre- and posttreatment samples, we show a complex milieu of parasites, including variants likely up-selected and down-selected by drug therapy. As with all surveys of diversity, sampling limitations prevent full discovery and differences in sampling effort can confound comparisons among samples, hosts, and populations. Here, we used ecological approaches of species accumulation curves and capture-recapture to estimate the number of variants we failed to detect in the population, and show that these methods enable comparisons of diversity before and after treatment, as well as between malaria populations. The combination of ecological statistics and massively parallel pyrosequencing provides a powerful tool for studying the evolution of drug resistance and the in-host ecology of malaria infections.