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

Dynamics of Plasmodium species and genotype diversity in adults with asymptomatic infection in Gabon

OBJECTIVES

We investigated the diversity and dynamics of Plasmodium infection in serially collected samples from asymptomatic participants of a clinical trial assessing the efficacy and safety of ivermectin in Gabon. We checked whether the baseline sample reflected the P. falciparum genotype and Plasmodium species diversity seen over 7 days of follow-up.

METHODS

Blood samples were collected at inclusion, every 8 hours until hour 72, daily until day 7, and on day 14. Plasmodium species was determined by qPCR and pfmsp1 length polymorphism was assessed for P. falciparum genotyping.

RESULTS

In 17/48 (35%) individuals, all pfmsp1 genotypes identified during the assessed period were detected at baseline; in 31/48 (65%), new genotypes were found during follow-up. Additional sampling at hour 24 allowed the identification of all genotypes seen over 7 days in 50% of the individuals. Ivermectin did not impact the genotype dynamics. Mixed Plasmodium spp. infections were detected in 28/49 (57%) individuals at baseline, and detection of non-falciparum infections during follow-up varied.

CONCLUSIONS

Our results reveal complex intra-host dynamics of P. falciparum genotypes and Plasmodium species and underscore the importance of serial sampling in clinical trials for antimalarial drugs with asymptomatically P. falciparum-infected individuals. This might allow a more accurate identification of genotypes in multiple infections, impacting the assessment of drug efficacy.

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.

PHARE: a bioinformatics pipeline for compositional profiling of multiclonal Plasmodium falciparum infections from long-read Nanopore sequencing data

BACKGROUND

The emergence of drug-resistant clones of Plasmodium falciparum is a major public health concern, and the ability to detect and track the spread of these clones is crucial for effective malaria control and treatment. However, in endemic settings, malaria infected people often carry multiple P. falciparum clones simultaneously making it likely to miss drug-resistant clones using traditional molecular typing methods.

OBJECTIVES

Our goal was to develop a bioinformatics pipeline for compositional profiling in multiclonal P. falciparum samples, sequenced using the Oxford Nanopore Technologies MinION platform.

METHODS

We developed the 'Finding P. falciparum haplotypes with resistance mutations in polyclonal infections' (PHARE) pipeline using existing bioinformatics tools and custom scripts written in python. PHARE was validated on three control datasets containing P. falciparum DNA of four laboratory strains at varying mixing ratios. Additionally, the pipeline was tested on clinical samples from children admitted to a paediatric hospital in the Central African Republic.

RESULTS

The PHARE pipeline achieved high recall and accuracy rates in all control datasets. The pipeline can be used on any gene and was tested with amplicons of the P. falciparum drug resistance marker genes pfdhps, pfdhfr and pfK13.

CONCLUSIONS

The PHARE pipeline helps to provide a more complete picture of drug resistance in the circulating P. falciparum population and can help to guide treatment recommendations. PHARE is freely available under the GNU Lesser General Public License v.3.0 on GitHub: https://github.com/Fippu/PHARE.

Genotyping Plasmodium falciparum gametocytes using amplicon deep sequencing

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Amplicon Deep Sequencing Reveals Multiple Genetic Events Lead to Treatment Failure with Atovaquone-Proguanil in Plasmodium falciparum

Atovaquone-proguanil (AP) is used as treatment for uncomplicated malaria, and as a chemoprophylactic agent against Plasmodium falciparum. Imported malaria remains one of the top causes of fever in Canadian returning travelers. Twelve sequential whole-blood samples before and after AP treatment failure were obtained from a patient diagnosed with P. falciparum malaria upon their return from Uganda and Sudan. Ultradeep sequencing was performed on the cytb, dhfr, and dhps markers of treatment resistance before and during the episode of recrudescence. Haplotyping profiles were generated using three different approaches: msp2-3D7 agarose and capillary electrophoresis, and cpmp using amplicon deep sequencing (ADS). A complexity of infection (COI) analysis was conducted. De novo cytb Y268C mutants strains were observed during an episode of recrudescence 17 days and 16 h after the initial malaria diagnosis and AP treatment initiation. No Y268C mutant reads were observed in any of the samples prior to the recrudescence. SNPs in the dhfr and dhps genes were observed upon initial presentation. The haplotyping profiles suggest multiple clones mutating under AP selection pressure (COI > 3). Significant differences in COI were observed by capillary electrophoresis and ADS compared to the agarose gel results. ADS using cpmp revealed the lowest haplotype variation across the longitudinal analysis. Our findings highlight the value of ultra-deep sequencing methods in the understanding of P. falciparum haplotype infection dynamics. Longitudinal samples should be analyzed in genotyping studies to increase the analytical sensitivity.

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.