Detection of novel Plasmodium falciparum haplotypes under treatment pressure in pediatric severe malaria
========================================================================================================

* Balotin Fogang
* Emilie Guillochon
* Claire Kamaliddin
* Gino Agbota
* Sem Ezinmegnon
* Maroufou Jules Alao
* Philippe Deloron
* Gwladys Bertin
* Antoine Claessens

## Abstract

**Background** In Africa, the clearance time for *P. falciparum* severe malaria varies significantly, likely due to the complexity of *P. falciparum* infections and the sequestration phenomenon exhibited by this parasite. This study aims to evaluate different methods to study intra-host dynamics of polygenomic infections during parasite clearance under antimalarial treatment. Additionally, it seeks to determine the association between parasite clearance rate following artesunate or quinine treatment and the genetic complexity of *P. falciparum* in Beninese children with severe malaria.

**Methods** Sixty-five *P. falciparum* severe malaria individuals diagnosed by microscopy and treated with artesunate or quinine were sampled every 8 hours for 24 hours. Using whole genome sequencing (WGS) data, we estimated the multiplicity of infection (MOI) with three algorithms (*Fws*, THE REAL McCOIL, and RoH). We then characterized the *P. falciparum* genetic complexity in WGS-identified polyclonal infections using amplicon sequencing (AmpSeq) on DNA extracted from plasma and from the red blood cells pellet.

**Results** AmpSeq demonstrated greater sensitivity in detecting multiple genomes within isolates compared to WGS methods. The MOI from AmpSeq was significantly higher in RBC pellets compared to plasma (2.4 vs 1.8 distinct microhaplotypes per isolate). However, at parasitaemia over 1000 parasites/uL, the same MOI was detected in both plasma and pellet samples in 85.4% of the isolates. We observed a high variability in parasite clearance rate among participants, but it was not associated with parasite MOI at diagnostic. Interestingly, in 60.9% of participants, previously undetected microhaplotypes appeared in circulation 16 hours after treatment initiation.

**Conclusion** These findings demonstrate that combining different haplotyping techniques effectively determines parasite genetic complexity. Additionally, plasma can be effectively used for parasite genotyping at sufficient parasitaemia levels. The parasite clearance rate of severe malaria is independent of parasite MOI. However, genotyping a single blood sample upon hospital admission does not capture the full spectrum of parasite genotypes present in the infection.

Keywords
*   . *Plasmodium falciparum*
*   genetic complexity
*   parasite clearance
*   severe malaria

## Introduction

*Plasmodium falciparum* is the predominant malaria species and a major contributor to child mortality in Africa. Many children infected with *P. falciparum* succumb to the disease before even reaching a hospital or clinic. Among those who are admitted with severe malaria and receive parenteral antimalarial treatment, approximately 5% do not survive, highlighting the urgent need for effective interventions [1]. Since 2011, the World Health Organisation (WHO) has recommended intravenous artesunate as the first-line antimalarial treatment for severe malaria [2], citing its potential to reduce severe malaria mortality rates by 22.5% compared to quinine [3]. Nevertheless, in cases where artesunate is unavailable, alternative severe malaria treatment such as intramuscular artemether or intravenous quinine are recommended [2].

Despite the widespread recommendation of Artemisinin-based Combination Therapies (ACTs) as the frontline treatment for malaria, reports of *P. falciparum* resistance to ACTs have risen in Asia [4, 5] and more recently in Africa [6, 7]. Indeed, artemisinin resistance is characterized by the slow clearance of parasites *in vivo*, which is a result of reduced drug susceptibility in ring-stage parasites driven by mutations in the *P. falciparum Kelch13* gene [8]. Since the introduction of artemisinin-based combination therapy (ACT) in 2005, no Pfkelch13 mutation associated with artemisinin resistance has been identified in Benin [9, 10].

Although the *PfKelch13* polymorphisms are the strongest predictor of *in vivo* parasite clearance time under treatment with artesunate, this time can be affected by other factors such as the age of the host, host-acquired immunity, initial parasitemia levels, and the developmental stages of parasite [11–13]. Although controversial, studies have shown a link between the genetic complexity of *P. falciparum* and delayed parasite clearance in Africa. Indeed, children infected with multiple strains had nearly a 3-fold increase in treatment failure compared to their age mates infected with a single strain [14]. However, low multiplicity of infection (MOI) at baseline were associated with detectable parasitaemia at 72h post-ACT treatment [15]. The association between antimalarial parasite clearance rate and parasite genetic complexity in individuals with severe malaria (SM) remains to be explored.

In SM, the delay in parasite clearance could also be attributed to the significant sequestration of mature *P. falciparum* parasitized erythrocytes within the tissue capillaries. Mathematical modelling using plasma histidine-rich protein 2 (PfHRP2) levels and parasitaemia have indicated that the number of sequestered parasites is substantially higher than what is detected in the peripheral blood [16–18]. Because the circulating parasitaemia is highly dependent on the parasite developmental stage, the amount of *Plasmodium* DNA in plasma is actually more accurate than parasitemia for diagnostic differentiating severe from uncomplicated malaria [19]. Other studies have also demonstrated that plasma can be used to detect and quantify the *Plasmodium* DNA by PCR [20, 21]. However, whether genomic DNA (gDNA) extracted from plasma could be used as material to reliably estimate the MOI is not known. Although the plasma may contain less parasite DNA than within RBCs, it is unaffected by sequestration and could theoretically contain parasite DNA released from prior schizont ruptures of all the parasite genotypes present within an individual. However, there is currently no evidence regarding the sensitivity of *Plasmodium* genotyping when using plasma samples as the DNA source.

Various sequencing approaches have recently been used to assess the complexity of *P. falciparum*, quantified as the number of parasite genotypes within an isolate [22]. These methods include whole-genome sequencing (WGS)-based techniques using diverse metrics such as *Fws* [23], THE REAL McCOIL [24], and Runs of Homozygosity (RoH) [25], which elucidate intra-host parasite diversity by analysing genome-wide variations across multiple loci. *Fws* is a metric characterizing within-host diversity and its relationship to population-level diversity, with theoretical values ranging from 0 to 1. The RoH indicates genetic relatedness between 2 or more genotypes within an isolate by identifying long blocks of haplotypes that have been inherited from the same parent. The REAL McCOIL is a statistical approach that uses data from all infections to simultaneously estimate allelic frequency and the number of distinct genotypes s in an isolate. Amplicon sequencing (AmpSeq), a technique targeting highly polymorphic loci, has emerged as a novel tool capable of detecting minority microhaplotypes with within-sample frequencies as low as 0.1% [26]. Indeed, the strength of sequencing targeted non-repetitive regions that harbour extensive single nucleotide polymorphisms (SNPs) lies in the fact that all SNPs within an amplicon are linked by a single sequence read, enabling direct microhaplotype identification [26]. Therefore, this study aims to evaluate the association between parasite clearance rate and multiplicity of infection in Beninese children undergoing treatment for severe malaria, using a variety of sequencing methods and biological materials.

## Methods

### Study design

Ethical clearance was obtained from “Comité d’Ethique de la Recherche CER_ISBA Benin” (clearance n°90, 06/06/2016 and clearance n°38, 16/05/2014). During the high malaria transmission season in Southern Benin, children under the age of six diagnosed with severe malaria infection at the “Centre Hospitalier et Universitaire de la Mère et de l’Enfant-Lagune (CHU-MEL) de Cotonou” were enrolled in the study in both 2014 and 2016, as described in Kamaliddin 2019 [27]. Prior to their participation, written informed consent from parents or guardians was obtained. Malaria infection was diagnosed in febrile patients by detecting *P. falciparum* on a Giemsa-stained thick blood smear. Parasite density was estimated by counting the number of parasites per at least 200 white blood cells (WBC) and assuming a standard total blood WBC count of 8000 WBC/µL. Diagnosis and classification of severe malaria were carried out in accordance with WHO guidelines for the management of severe malaria [2]. In brief, cerebral malaria (CM) was characterised by impaired consciousness, indicated by a Blantyre score < 3, after excluding other causes of coma. Severe anaemia (SMA) was identified by a haemoglobin level < 5 g/dL or a haematocrit level < 15%. Severe non-cerebral malaria (SNCM) was defined by the features of severe malaria excluding CM involvement. Furthermore, clinical appreciation was at the discretion of the admitting physician.

Approximately 5 mL of blood were collected from all participants using EDTA-coated tubes before the administration of quinine (in 2014) or artesunate (in 2016) therapy (H0). Subsequently, three additional blood samples were taken at +8 hours (H8), +16 hours (H16), and +24 hours (H24) post-treatment initiation to capture parasite dynamics.

### DNA extraction

Plasma was separated from total blood cells by centrifugation at 1500 rpm for 5 min and immediately stored at −20°C. Red blood cell (RBC) pellet was obtained by depleting white blood cells using a gradient-based separation technique Ficoll (GE Healthcare Life Science). *P. falciparum* genomic DNA (gDNA) was extracted from both 200 µL of red blood cell (RBC) pellets and 200 µL of plasma using the DNEasy Blood kit (QIAGEN) according to the manufacturer’s instructions.

### Complexity of infection using whole genome sequencing (WGS) data

WGS was performed by the Malaria Genomic Epidemiology Network (MalariaGEN) and the samples are part of the *Plasmodium falciparum* Community Project accessible at [https://www.malariagen.net/project/p-falciparum-community-project/](https://www.malariagen.net/project/p-falciparum-community-project/). The raw reads are available on the ENA server under the accession numbers specified in the Supplementary Table 1. Variant Calling Files (VCF) was generated by MalariaGEN as previously described [28]. Briefly, reads were mapped to *P. falciparum* 3D7 v3 reference genome using bwa mem, and the resulting BAM files were subjected to cleaning using Picard tools and GATK. SNPs and indels were called using GATK HaplotypeCaller, with only the core genome considered for analysis.

For downstream analysis bcftools v1.13 [29] and bedtools v2.30 [30] were used for file manipulation. Samples with at least 50% coverage of the genome at 5x were considered. Fraction of within-sample (*Fws*) metrics were computed using the moimix R package, as previously described ([https://github.com/bahlolab/moimix](https://github.com/bahlolab/moimix)). An *Fws* value < 0.95 was indicative of a polyclonal infection. THE REAL McCOIL categorical method was used to estimate the complexity of the infection (MOI) as described by Chang et al. [24] ([https://github.com/Greenhouse-Lab/THEREALMcCOIL](https://github.com/Greenhouse-Lab/THEREALMcCOIL)). Long runs of homozygosity (RoH) and the analysis of heterozygosity in mixed samples were carried out using a custom Python script based on Pearson et al. [25].

### DNA preparation and amplicon sequencing

Three AmpSeq markers, including the genes for conserved plasmodium membrane protein (*cpmp*, PF3D7_0104100), conserved plasmodium protein (*cpp*, PF3D7_1475800) and apical membrane antigen (*ama1-D3*, PF3D7_1133400) were amplified using nested-PCR, as described in [31]. Briefly, primary PCRs were conducted in multiplex for *cpmp/ama1-D3* and monoplex for *cpp*. Nested PCRs were subsequently performed individually for each marker. All amplifications were carried out using the KAPA HiFi Hot Start Ready Mix (Roche) on an Eppendorf Mastercycler Nexus thermocycler. The quality and quantity of nested-PCR products were assessed through gel electrophoresis. Following this, nested-PCR products from each sample were combined in the following proportions: 11 μL of cpmp, 8 μL of cpp, and 4 μL of ama1-D3, for a total of 25 μL. PCR primers and amplification conditions are shown in the Supplementary Table 2.

Libraries preparations and amplicon sequencing were performed by the Genseq platform of the Labex CeMEB (Montpellier). Merged PCR products were purified, sequencing adapters and multiplexing indices were linked by PCR. Indexed products were purified, pooled and quality was verified by electrophoresis using The Fragment Analyzer system (Agilent Biotechnologies). PCR products resulting from RBC pellet or plasma gDNA were separated in two different plates, to prevent the impact of the difference in DNA concentration on the sequencing quality. Finally, amplicon libraries were sequenced on an Illumina MiSeq system in paired-end mode (2 x 300 cycles, 300 bp). The raw AmpSeq reads are available on the Zenedo ([https://zenodo.org/records/13224728](https://zenodo.org/records/13224728)), with corresponding individual accession numbers provided in Supplementary Table 3.

Markers haplotypes were reconstructed using HaplotypR R package [26], accessible at [https://github.com/lerch-a/HaplotypR](https://github.com/lerch-a/HaplotypR). Briefly, the quality of raw reads data was assessed using FastQC. Reads corresponding to each marker were demultiplexed based on primers sequences, with subsequent truncation of primers. Reads were trimmed following sequence quality and then fused. SNPs and microhaplotypes were called with default parameters, requiring a minimum coverage of 3 reads per microhaplotype, at least 25 reads coverage per sample, and a within-host microhaplotype frequency ≥ 0.1%. Microhaplotype sequences that differ by less than two SNPs were deemed identical. The relative frequency of each microhaplotype in an isolate was calculated by dividing the number of reads for each microhaplotype by the total number of reads for that isolate. For downstream analysis, marker showing the highest number of microhaplotypes in an individual was retained. The same marker was chosen for the four time points in an individual. The raw data for AmpSeq are presented in Supplementary Table 3.

### Statistical analysis

All statistical analyses were performed using R software version 4.3.2 [32]. Median comparisons of two quantitative variables were conducted using the Mann–Whitney test, while mean comparisons were assessed using the unpaired t-test. The correlation between two quantitative variables was assessed using the Spearman correlation test. For categorical variables, proportions were analysed using the Chi-squared test. We used linear interpolation to estimate the time required for a 50% reduction in parasitaemia, referred as the extrapolated parasite clearance half-life (PC50%). The clearance curve for each microhaplotype (parasite genotype detected by AmpSeq) was calculated by multiplying the parasitaemia at each time point by the relative frequency of each microhaplotype. Shannon entropy was used to characterize the complexity of microhaplotypes within an isolate obtained from AmpSeq.

## Results

### Characteristics of the study population

During the years 2014 and 2016, a total of 65 children aged 0 to 5 years old suffering from severe malaria and admitted at CHU-MEL in Benin were recruited. Of these participants, 17 received quinine in 2014, while 48 received artesunate in 2016 (Figure 1). These individuals were classified as 26 cerebral malaria (CM), 16 severe non-cerebral malaria (SNCM) and 23 severe malaria anaemia (SMA) (Supplementary Table 1). The median age and sex ratios of participants were comparable between quinine and artemisinin treated groups. However, the geometric mean parasitaemia in 2014 was significantly higher than in 2016 (p< 0.0001). The demographic and parasitological characteristics of the participants at enrolment are provided in Table 1. Blood samples were also collected on arrival at hospital and every 8 hours for 24 hours (Figure 1).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/20/2024.08.20.24312296/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2024/08/20/2024.08.20.24312296/F1)

Figure 1. Study flow diagram.
Sixty-five severe malaria patients were sampled at the start of treatment and every 8 hours for 24 hours afterwards. Out of 230 isolates submitted for whole genome sequencing, 192 yielded high-quality genomes. The multiplicity of infection (MOI) was determined from parasite whole genome sequences. From 34 individuals with at least one polyclonal timepoint, amplicon sequencing (AmpSeq) was performed on plasma and red blood cell pellet samples. A total of 76 RBC pellet isolates and 74 plasma gDNA isolates from 23 individuals were successfully genotyped.

View this table:
[Table 1.](http://medrxiv.org/content/early/2024/08/20/2024.08.20.24312296/T1)

Table 1. Demographic and clinical background of the study participants

### Determining the Multiplicity of Infection (MOI) based on Whole Genome Sequencing (WGS)

We first aimed to characterise the MOI using whole genome sequencing. A total of 230 RBC pellet samples at each timepoint were whole genome sequenced, resulting in 192 good-quality genomes (Figure 1 and Supplementary Table 1). To assess the MOI from these genomes, three algorithms were tested: *FWS*, THE REAL McCOIL and Runs of Homozygosity (RoH). All three algorithms showed robust correlation (Appendix 1). From here onwards any genome with *FWS*> 0.95 was considered monoclonal.

### Determining the MOI based on AmpSeq using gDNA from RBC pellets and plasma

We used the amplicon sequencing approach (AmpSeq) to assess MOI from paired RBC pellets and plasma samples obtained from the same individuals. The rationale was that the RBC pellet represents circulating parasites at a certain timepoint, while plasma could contain parasite DNA from previous schizont ruptures that is independent of sequestration. As described in the methods section, three markers were used (*cpmp*, *cpp*, *ama1*) and the marker showing the highest number of microhaplotypes in an individual was retained (*cpmp, cpp* and *ama1* in 34.8%, 34.8% and 30.4% of individuals, respectively, Supplementary Table 3). AmpSeq was successfully performed at each timepoint sample from 23 individuals with a polyclonal infection, as determined by WGS methods. 69 isolates were successfully genotyped in both RBC pellets and plasma gDNA, 7 were genotyped in RBC pellets gDNA only and 5 in plasma gDNA only (Supplementary Table 1).

In terms of sequence diversity, a total of 65 unique microhaplotypes were detected across all isolates. Among these, 98.5% (64) were identified in pellets, while 70.8% (46) were found in plasma. Furthermore, 19 microhaplotypes (29.2%) were exclusively present in pellets, while only one microhaplotype was only found in plasma. The mean MOI per isolate was significantly higher in RBC pellets compared to plasma (2.4 and 1.8 for pellets and plasma respectively, p = 0.0042) (Figure 2). Among isolates from the same individual, 66.7% (46/69) exhibited an identical number of genotypes detected when using gDNA from both RBC pellets and plasma. Furthermore, at parasitaemia over 1000 parasites/uL, the same MOI was detected in both plasma and pellet samples in 85.4% of the isolates (Supplementary Figure 1). This indicates the influence of parasite load on the sensitivity of AmpSeq when using plasma as the source of DNA. However, these data also show that plasma can be effectively used for parasite genotyping at relatively high parasitaemia levels.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/20/2024.08.20.24312296/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2024/08/20/2024.08.20.24312296/F2)

Figure 2. Comparison of the multiplicity of the infection determined using AmpSeq with gDNA from paired RBC pellet and plasma samples. The MOI per isolate, indicated as a colour scale, is the number of distinct microhaplotypes detected per isolate.

By visualizing the evolution of clonality at each time point (0h, 8h, 16h and 24h) during the 24h after treatment initiation (Figure 3), 7 infections showed the same MOI profile from RBC pellet and plasma over time (Figure 3A), while in 16 infections more microhaplotypes were detected in RBC pellet than in plasma (Figure 3B). Taken together, AmpSeq using gDNA from RBC pellets is more sensitive compared to AmpSeq using gDNA from plasma at low range parasitaemia. From here onwards, the AmpSeq MOI was defined as the number of microhaplotypes detected with gDNA from RBC pellets.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/20/2024.08.20.24312296/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2024/08/20/2024.08.20.24312296/F3)

Figure 3. Relative proportion of microhaplotypes within patients, from AmpSeq using paired gDNA extracted from RBC pellets and plasma samples.
A) Patients with the same microhaplotype profile between plasma and pellet gDNA. B) Patients with more microhaplotypes (higher MOI) detected within RBC pellet gDNA compared to plasma. Each colour indicates a unique microhaplotype within the whole population.

### Comparison of AmpSeq and WGS-based methods for determining MOI from pellet gDNA

A total of 68 isolates had MOI data available for both AmpSeq and WGS. We found a good correlation between AmpSeq and the different WGS-based methods (r= −0.61, r= 0.47 and r= −0.59 for *Fws*, THE REAL McCOIL and RoH, respectively) (Figure 4). In addition, 80.8% and 79.3% of the isolates showed concordant results (either monoclonal or polyclonal) between AmpSeq and *Fws*, or AmpSeq and RoH, respectively (Figure 4B). However, the AmpSeq method was able to detect up to 6 genotypes in an isolate, indicating greater sensitivity of AmpSeq in detecting multiple genomes in an isolate compared to WGS-based methods (Figure 4B). Furthermore, AmpSeq was the only method capable of quantifying the proportion of each genotype identified within an isolate. Together, these findings showed that the combination of different haplotyping techniques can be used to characterise genetic complexity quantitatively and qualitatively within an isolate. The MOI of a polyclonal isolate was then defined from the AmpSeq data.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/20/2024.08.20.24312296/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2024/08/20/2024.08.20.24312296/F4)

Figure 4. Comparison of AmpSeq to different WGS-based algorithms for determining the MOI.
A) Correlation between *Fws* and AmpSeq. B) Proportion of isolates according to clonality in different genotyping techniques. Only individuals with at least one timepoint at *Fws* <0.95 were included in this analysis.

### Variability in parasite clearance rate among individuals and high rate of rebound in parasitaemia during antimalarial treatment

Having established the most sensitive method for MOI determination, we then aimed to characterise the clearance rate within the first 24 hours of treatment. Median parasite clearance time was faster in children treated with artesunate compared to those treated with quinine (Figure 5A), with the mean PC50% of 8.3 and 12.2 hours, respectively (Figure 5B). There was a high variability in parasite clearance half-life within both treatment groups, ranging from 4.1 to 23.3 hours for artesunate and 4.0 to 24.0 hours for quinine (Figure 5B). Of the 48 individuals treated with artesunate, 95.6% (43/45) had a PC50% below 16 hours, compared to 75% (12/16) of individuals treated with quinine (Figure 5B). The PC50% could not be determined for four individuals (023, 014H, 029H, 063H) due to a non-significant decrease in parasitaemia (<50% reduction) between the first and last time points (Supplementary Data 2). Parasite clearance curves are within then same range as previous studies using quinine or artesunate [33] (Supplementary Data 2). These curves indicated substantial variation in initial clearance rates between individuals, highlighting the complex dynamics of parasite response to treatment in the early stage of infection management.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/20/2024.08.20.24312296/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2024/08/20/2024.08.20.24312296/F5)

Figure 5. Parasite clearance curves of severe malaria cases treated with Artesunate or Quinine.
A) Parasite clearance curve according to antimalarial treatment. B) Time for 50% reduction in parasitaemia according to treatment. C) parasite clearance curve highlighting the rebound in parasitaemia post-antimalarial treatment.

We took advantage of the time-series to investigate intra-host competition in polygenomic infections. Interestingly, we observed a ‘rebound’ in the parasitemia determined by microscopy reading. We defined rebound as an increase in parasite densities by at least 100 parasites/µL between two time points, and observed this phenomenon in 34.8% (20/65) of individuals (Figure 6C). Specifically, within the artesunate-treated group, 35.4% (17/48) of individuals experienced a rebound, while in the quinine-treated group, the rebound rate was 17.6% (3/17). These data indicate the potential release of sequestered parasites during anti-malarial treatment.

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/20/2024.08.20.24312296/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2024/08/20/2024.08.20.24312296/F6)

Figure 6. Dynamic of parasite clearance to parasite multiplicity of the infection.
A) Dynamic of Multiplicity of Infection (MOI) determined by AmpSeq throughout the antimalarial treatment. Parasite clearance rate according to clonality at enrolment (H0) for artesunate (B) quinine (C) treatment groups.

### *P. falciparum* clearance rate is not related to MOI at the onset of the treatment

To determine whether the variability in the clearance rate of *P. falciparum* depends on the number of parasite genotypes in the infection, the dynamics of parasite clearance were compared between individuals with polyclonal and monoclonal infections during the first 24h of the severe malaria management under treatment pressure. Independently of the treatment administered, the mean MOI was relatively stable throughout the 24-hour period (Figure 6A). In individuals treated with artesunate, *P. falciparum* clearance showed no significant difference between polyclonal and monoclonal infections (Figure 6B). In contrast, although non-significant, treatment with quinine resulted in faster *P. falciparum* clearance in polyclonal infections (mean PC50% of 4.3) compared to monoclonal infections (mean PC50% of 12.4) (p= 0.1393) (Figure 6C). Part of this difference may be attributed to the higher median parasitaemia at enrolment in individuals treated with quinine and having monoclonal infections compared to those with polyclonal infections (2440000 vs 1520000 p/µL, p= 0.1924). In addition, there was no correlation between MOI at enrolment and parasite clearance rate (r= −0.06, p= 0.6494). These data show that the *P. falciparum* clearance rate does not depend on the initial parasite MOI.

### Sporadic detection of novel microhaplotypes during anti-malarial treatment

To measure the complexity of microhaplotypes within each isolate, we measured the Shannon entropy, a metric dependent on the number and proportion of each microhaplotype within an isolate. Despite the overall parasitaemia decreasing by an average of 335-fold over 24 hours, the microhaplotype entropy only slightly decreased between H0 and H8 (p= 0.0388), and remained stable afterwards (p= 0.6093) (Supplementary Figure 2). This stable MOI is explained by the detection of novel microhaplotypes that were not in circulation at H0. These ‘sporadic’ microhaplotypes were detected in over 60.9% (14/23) of participants, predominantly occurring from 16 hours after the start of treatment (Figure 7), likely indicating the release of sequestered parasites that were previously undetectable. Additionally, half of participants with sporadic microhaplotypes (007H, 008, 012H, 013H, 016H, 024h, 064H, and 066H) had monoclonal infections at H0 but displayed polyclonal infections at later time points (Figure 7A). However, most of the sporadic microhaplotypes were very minor genotypes (Figure 7A).

![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/20/2024.08.20.24312296/F7.medium.gif)

[Figure 7.](http://medrxiv.org/content/early/2024/08/20/2024.08.20.24312296/F7)

Figure 7. Subpopulation clearance curve during anti-malarial treatment, with (A) or without (B) sporadic detection of novel microhaplotypes.
The left panel of each plot shows the relative abundance of microhaplotype subpopulations throughout the antimalarial treatment period. The right panel of each plot shows parasite densities of these subpopulations, calculated by multiplying the total parasitaemia by the relative abundance of the microhaplotypes. The clearance curves for the total parasite density are shown as a background surface (in grey). Each color represents a unique microhaplotype within the population. Sporadic detection is defined as the appearance of a new microhaplotype that was not present at H0.

We inferred the parasitaemia of each circulating genotype based on the total parasitaemia and the proportion of each microhaplotype. Notably, only one individual (067H) harbouring sporadic microhaplotypes showed a rebound in microhaplotype parasitaemia, indicating that the rebound in total parasitaemia observed was unrelated to the release of new genotypes in circulation after the start of the treatment.

Taken together, we here showed that, from blood samples taken on arrival at hospital (H0), it is usually not possible to genotype all *P. falciparum* microhaplotypes from an infection. About two third of the time, one or more distinct strain will appear within 24 hours, presumably as a result of release of sequestered *P. falciparum*.

### Parasite clearance is not correlated with drug resistance allele frequencies

We tested the hypothesis that drug sensitive genomes get cleared faster, by assessing allelic frequencies of resistance markers *crt*, *mdr1, kelch13, dhfr* and *dhps* over the 24 hours period. Of the 45 individuals for whom baseline drug resistance data was available, 95.6% (43/45) carried chloroquine-resistant allele (K76T) (Supplementary Figure 3A). However, the two individuals showing a mixture of sensitive and resistant alleles came from the group treated with artesunate (Supplementary Figure 3A). The frequency of the multidrug resistance gene 1 (*mdr1*) remained stable in the population throughout the 24 hours treatment period (Supplementary Figure 3B). Only one individual with a k13 mutation (V589I) at the time of enrolment was found in the study population (Supplementary Figure 3C). All isolates carried a triple mutation of the *dhfr* gene (N51I, C59R and S108N) (Supplementary Figure 3D). Similarly, 99% of isolates carried the A437G mutation in the *dhps* gene, with 23.2% and 0% a second mutation (Supplementary Figure 3E). Overall, we observed no alteration in genotypes across all analysed resistance markers within the same individual during antimalarial treatment.

## Discussion

In recent years, several approaches have been used to study the intra-host dynamics of polygenomic infections, including WGS which is the most comprehensive approach for genomic epidemiology, providing a complete picture of genetic variation. However, molecular markers, mainly AmpSeq, which combines several highly polymorphic antigen markers, is now the method of choice for high-resolution detection of MOI [22, 34]. We first aim in this study to evaluate different methods to study intra-host dynamics of polygenomic infections during a course of parasite clearance with an antimalarial treatment. Using the WGS data, we observed a robust correlation between the three algorithms (*Fws*, THE REAL McCOIL and RoH) used to determine the clonality of an isolate, in particular higher between *Fws* and RoH with 94.3% consistent results between the two metrics. However, 16.1% of polyclonal isolates according to *FWS* were classified as monoclonal according to THE REAL McCOIL. These data indicate that both the Fws and RoH algorithms exhibit greater sensitivity in detecting multiple infections compared to THE REAL McCOIL. Indeed, THE REAL McCOIL is a more conservative approach, which provides estimates that minimize overestimation of MOI by considering the likelihood of observed genetic variation being due to multiple parasite strains rather than sequencing errors or technical artefacts [24]. In this study, more than 79% of the isolates showed concordant results (either monoclonal or polyclonal) between AmpSeq and Fws or RoH. However, the AmpSeq method was able to detect up to 6 genotypes in an isolate, indicating greater sensitivity of AmpSeq in detecting multiple genomes in an isolate compared to WGS-based methods. Amplicon sequencing targeting highly polymorphic loci has demonstrated the capability to detect minority haplotypes with frequencies as low as 0.1% [26].

Previous studies have demonstrated that plasma can be used to detect and quantify the *Plasmodium* DNA by PCR [20, 21]. However, these studies also indicated that qPCR using whole blood DNA exhibits greater sensitivity than qPCR with plasma DNA. In line with the sequestration phenomenon of *P. falciparum* as shown by quantify the plasma levels of HRP2 or parasite DNA [18, 19], we hypothesized that genotyping parasites using DNA extracted from plasma might detect sequestered *P. falciparum* genotypes. This hypothesis assumes that plasma may contain parasite DNA released from prior schizont ruptures of all the parasite genotypes, which would be unaffected by sequestration at the time of blood sample collection. However, our findings revealed significantly higher microhaplotype complexity in the RBC pellet compared to plasma. These results indicate AmpSeq using gDNA from RBC pellets is more sensitive compared to AmpSeq using gDNA from plasma. To what extent the RBC pellet contains DNA from prior schizont ruptures remains to be determined. Consistent with previous findings where low levels of parasitemia was more effectively detected by qPCR using parasite DNA from whole blood compared to plasma [20], our study also found that parasite load influences the sensitivity of AmpSeq when using plasma as the DNA source. Indeed, the median parasitaemia of isolates with the higher MOI in the RBC pellet compared to plasma was 34.4 fold lower compared to the median parasitaemia of isolates with the same MOI in both biological materials. Interestingly, at parasitaemia levels exceeding 1000 parasites/µL, the same multiplicity of infection (MOI) was detected in both plasma and pellet samples in 85.4% of the isolates. This indicates that plasma can be effectively used for parasite genotyping in retrospective studies.

We observed high variability in parasite clearance half-life within both treatment groups with faster parasite clearance rate in children treated with artesunate compared to those treated with quinine. Artesunate is a fast-acting antimalarial drug with a pronounced effect on ring-stage parasites, whereas quinine primarily targets mature trophozoites [35]. Rapid clearance results from artesunate action on the circulating ring-stage parasites and their subsequent removal predominantly by the spleen, preventing the parasite sequestration [36, 37]. In line with our findings, a recent study of Ghanaian children showed a prolonged delay in parasite clearance among children with severe malaria, even after three days of treatment with artesunate [38]. This increased time to parasite clearance was associated with age, low haemoglobin levels and high number of previous malaria diagnoses [38]. Other studies have noted variations in parasite clearance rates in children treated for uncomplicated falciparum malaria with artemisinin-based combination therapy, both in Kenya and Tanzania [39, 40]. These findings suggest that factors other than drug resistance may contribute to slower parasite clearance in certain infections, as there is no evidence of artemisinin resistance in these regions. It has been shown that, in areas where artemisinin resistance is not present, the parasite clearance time can be affected by factors such as the age of the host, host acquired immunity, initial parasitemia levels, and the developmental stages of parasite [11–13]. Furthermore, although controversial, studies have shown a link between the genetic complexity of *P. falciparum* and delayed parasite clearance in Africa [14, 15]. Our data show that the *P. falciparum* clearance rate in children with severe malaria independent of parasite MOI at the onset of treatment. Notably, *P. falciparum* clearance was faster in polyclonal infections in quinine-treated children. However, this difference was attributed to significantly higher parasitemia at enrolment in individuals treated with quinine.

Interestingly, during antimalarial treatment, a high rate of “rebound” in total parasitemia and specific parasite subpopulations parasitemia was observed. Specifically, an increase in total parasitemia between two time points was detected in 34.8% of individuals. In four individuals (023, 031, 063H, and 064H), minor genotypes present at baseline (H0) increased in relative frequency, becoming the dominant microhaplotype. Additionally, sporadic detection of microhaplotypes occurred in over 60.9% of participants, predominantly from 16 hours after treatment initiation. Half of the participants with sporadic microhaplotypes had monoclonal infections at baseline but exhibited polyclonal infections at subsequent time points. These findings potentially highlight release of sequestered *P. falciparum* genotypes during antimalarial treatment [34] Interestingly, none of these sporadic microhaplotypes were detected in pellet or plasma samples at H0, despite the much higher parasite density at that timepoint that would have facilitated detection of a rare genotype. Either these microhaplotypes recently emerged from the liver, with no or little parasite DNA in the plasma yet, or, probably more likely, these sporadic microhaplotypes are present at extremely low parasitaemia making them undetectable in plasma, and were all sequestered at H0.

Additionally, a study by Marks et al. (2005) reported a parasitological rebound effect in infants treated with a single dose of sulfadoxine-pyrimethamine, attributed to the selection of drug-resistant parasites shortly after drug clearance [41]. However, no mutation associated with artemisinin resistance has yet been identified in Benin [9, 10]. In addition, no alterations were observed in allele frequencies across all analysed resistance markers within the same individual post-antimalarial treatment. Two other studies using AmpSeq detected sporadic microhaplotypes or an increase in the frequency of minor microhaplotypes during the treatment of uncomplicated *P. falciparum* infections in regions with no signs of artemisinin resistance [39, 40]. In *falciparum* malaria, parasitized red blood cells circulate in the peripheral blood for only one-third of the 48-hour asexual cycle [42]; for the remainder of the cycle, they are sequestered in the venules and capillaries. For these infections, the peripheral blood parasitemia at a given time point does not reflect the total parasite burden.

Limitations of this study include the lack of sampling timepoints beyond 24 hours post-treatment, which affected the accurate determination of parasite clearance half-life for some individuals using the standard WWARN method [43]. This was due to the fact that the parasitemia reduction rate was less than 50% of the initial parasitemia at the last recorded timepoint. Therefore, we used linear interpolation to estimate the time required for a 50% reduction in parasitaemia, referred as the extrapolated parasite clearance half-life (PC50%). However, unlike therapeutic efficacy studies typically focusing on day 3, 7 and 21/28, we were able to monitor the dynamics of *P. falciparum* complexity during the initial management of severe malaria.

This study demonstrates that combining different haplotyping techniques can effectively determine genetic complexity. However, AmpSeq is the most sensitive approach and should be the first choice for genotyping to monitor the evolution of genotype clearance following antimalarial treatment. Additionally, plasma can be effectively used for parasite genotyping at sufficient parasitaemia levels. We demonstrated that genotyping a blood sample from a patient on arrival at hospital is typically not sufficient to detect all parasites genotypes present in the infection.

## Conflicts of interest

No conflicts of interest were reported by the authors.

## Financial support

This work was funded by grants from the French National Research Agency (18lJCE15lJ0009lJ01), the ATIP-Avenir programme and by Institut Merieux. The “SPFlJPost doc en France” programme of the “Fondation pour la Recherche Médicale” funded the author “Balotin Fogang” (SPF202209015889).

## Ethical approval

Ethical clearance was obtained from “Comité d’Ethique de la Recherche CER_ISBA Benin” (clearance n°90, 06/06/2016 and clearance n°38, 16/05/2014).

## Supporting information

Supplementary Material [[supplements/312296_file03.docx]](pending:yes)

## Data Availability

All data produced in the present work are contained in the manuscript or are available online at the European Nucleotide Archive (ENA) under the study accession number PRJEB2136 and on Zenodo at https://zenodo.org/records/13224728.

[https://www.ebi.ac.uk/ena/browser/view/PRJEB2136](https://www.ebi.ac.uk/ena/browser/view/PRJEB2136) 

[https://zenodo.org/records/13224728](https://zenodo.org/records/13224728) 

## Author contributions

Conceptualization by A.C and GI.B. Field work and sample processing by C.K, S.E and G.A. Data analysis by B.F, E.G. Writing of original draft by B.F, E.G. Reviewing and editing of manuscript by B.F., E.G, C.K, GI.B, A.C. Funding acquisition by A.C and GI. B.

## Data availability

All whole-genome sequencing (WGS) data generated and analyzed in this study have been deposited in the European Nucleotide Archive (ENA) under the study accession number PRJEB2136, with individual accession numbers in Supplementary Table 1. The raw reads from the Amplicon sequencing (AmpSeq) data are publicly available on Zenodo ([https://zenodo.org/records/13224728](https://zenodo.org/records/13224728)), with corresponding individual accession numbers provided in Supplementary Table 3. The raw data and complete AmpSeq results are included as Supplementary Tables 1 and 2, respectively.

## Description of supplementary Material

**Supplementary Figure 1. Parasitaemia according to the concordance between MOI from paired RBC pellet and plasma samples, determined by AmpSeq.** “Pellet-plasma match” are isolates with identical MOI in plasma and pellet samples. “Low in plasma” are isolates in which the measured MOI was lower in plasma compared to pellet samples. There were only 2 isolates with MOI higher in the plasma than in the pellet, which are not shown on this figure.

**Supplementary Figure 2. Microhaplotype entropy over time during antimalarial treatment.**

**Supplementary Figure 3. Frequency of antimalarial drug resistance markers throughout treatment.**

**Supplementary Table 1. Epidemiological data of study participants**.

**Supplementary Table 2. PCR primers sequences and amplification conditions.** List of forward and reverse primers used for AmpSeq and the amplification program.

**Supplementary Table 3. Amplicon Sequencing results.** AmpSeq read counts and frequencies, for three markers (*cpmp*, *ama1-D3*, *cpp*) used in both plasma and red blood cell pellets. The table identifies the selected marker for each individual, with the marker displaying the highest number of microhaplotypes being retained. The last column indicates ‘sporadic’ microhaplotypes, defined as novel microhaplotypes not detected at H0.

**Appendix 1. Determining the Multiplicity of Infection (MOI) based on Whole Genome Sequencing (WGS).**

**Supplementary Data 1. Non-reference allele frequency (NRAF) distribution across all heterozygous SNPs (left) and heterozygosity calculated in 100-kb bins (right) to highlight RoH.**

**Supplementary Data 2. Individuals parasite clearance curve by timepoints post-antimalarial treatment. Quinine (A) and artesunate (B) treated groups.**

## Acknowledgments

We would like to thank the participants and their parents/guardians for granting approval to participate in this study. We thank the staff of the paediatric wards at the CHUMEL hospital in Cotonou and at Saint Joseph medical centre in Sô-Tchanhoué, Benin, for their help. Our thanks go to MalariaGen for the whole genome sequencing and to the GenSeq platform for the amplicon sequencing.

*   Received August 20, 2024.
*   Revision received August 20, 2024.
*   Accepted August 20, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/)

## References

1.  1.White NJ. Severe malaria. Malar J 2022;21:284.
    
    

2.  2.WHO. Management of severe malaria: a practical handbook. 3rd ed. Geneva: World Health Organization. [https://iris.who.int/handle/10665/79317](https://iris.who.int/handle/10665/79317) (2012, accessed 18 June 2024).
    
    

3.  3.Dondorp AM, Fanello CI, Hendriksen IC, Gomes E, Seni A, et al. Artesunate versus quinine in the treatment of severe falciparum malaria in African children (AQUAMAT): an open-label, randomised trial. Lancet 2010;376:1647–1657.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(10)61924-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21062666&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000284451000030&link_type=ISI) 

4.  4.Phyo AP, Nosten F, Phyo AP, Nosten F. The Artemisinin Resistance in Southeast Asia: An Imminent Global Threat to Malaria Elimination. In: Towards Malaria Elimination - A Leap Forward. IntechOpen. Epub ahead of print 18 July 2018. DOI: 10.5772/intechopen.76519.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.5772/intechopen.76519&link_type=DOI) 

5.  5.Pluijm RW van der, Imwong M, Chau NH, Hoa NT, Thuy-Nhien NT, et al. Determinants of dihydroartemisinin-piperaquine treatment failure in Plasmodium falciparum malaria in Cambodia, Thailand, and Vietnam: a prospective clinical, pharmacological, and genetic study. Lancet Infect Dis 2019;19:952–961.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(19)30391-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31345710&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 

6.  6.Assefa A, Fola AA, Tasew G. Emergence of Plasmodium falciparum strains with artemisinin partial resistance in East Africa and the Horn of Africa: is there a need to panic? Malar J 2024;23:34.
    
    

7.  7.Balikagala B, Fukuda N, Ikeda M, Katuro OT, Tachibana S-I, et al. Evidence of Artemisinin-Resistant Malaria in Africa. N Engl J Med 2021;385:1163–1171.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2101746&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34551228&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 

8.  8.Ariey F, Witkowski B, Amaratunga C, Beghain J, Langlois A-C, et al. A molecular marker of artemisinin-resistant Plasmodium falciparum malaria. Nature 2014;505:50–55.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature12876&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24352242&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000329163300021&link_type=ISI) 

9.  9.Kpemasse A, Dagnon F, Saliou R, Yarou Maye AS, Affoukou CD, et al. Efficacy of Artemether-Lumefantrine for the Treatment of Plasmodium falciparum Malaria in Bohicon and Kandi, Republic of Benin, 2018–2019. Am J Trop Med Hyg 2021;105:670–676.
    
    

10. 10.L’Episcopia M, Doderer-Lang C, Perrotti E, Priuli GB, Cavallari S, et al. Polymorphism analysis of drug resistance markers in Plasmodium falciparum isolates from Benin. Acta Trop 2023;245:106975.
    
    

11. 11.Ataide R, Ashley EA, Powell R, Chan J-A, Malloy MJ, et al. Host immunity to Plasmodium falciparum and the assessment of emerging artemisinin resistance in a multinational cohort. Proc Natl Acad Sci U S A 2017;114:3515–3520.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE0LzEzLzM1MTUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOC8yMC8yMDI0LjA4LjIwLjI0MzEyMjk2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

12. 12.Intharabut B, Kingston HW, Srinamon K, Ashley EA, Imwong M, et al. Artemisinin Resistance and Stage Dependency of Parasite Clearance in Falciparum Malaria. J Infect Dis 2019;219:1483–1489.
    
    

13. 13.WWARN Parasite Clearance Study Group, Abdulla S, Ashley EA, Bassat Q, Bethell D, et al. Baseline data of parasite clearance in patients with falciparum malaria treated with an artemisinin derivative: an individual patient data meta-analysis. Malar J 2015;14:359.
    
    

14. 14.Kyabayinze DJ, Karamagi C, Kiggundu M, Kamya MR, Wabwire-Mangen F, et al. Multiplicity of Plasmodium falciparum infection predicts antimalarial treatment outcome in Ugandan children. Afr Health Sci 2008;8:200–205.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20589125&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 

15. 15.Topazian HM, Moser KA, Ngasala B, Oluoch PO, Forconi CS, et al. Low Complexity of Infection Is Associated With Molecular Persistence of Plasmodium falciparum in Kenya and Tanzania. Front Epidemiol;2. Epub ahead of print 6 June 2022. DOI: 10.3389/fepid.2022.852237.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fepid.2022.852237&link_type=DOI) 

16. 16.Cunnington AJ, Bretscher MT, Nogaro SI, Riley EM, Walther M. Comparison of parasite sequestration in uncomplicated and severe childhood Plasmodium falciparum malaria. J Infect 2013;67:220–230.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jinf.2013.04.013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23623771&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000322737300007&link_type=ISI) 

17. 17.Dondorp AM, Desakorn V, Pongtavornpinyo W, Sahassananda D, Silamut K, et al. Estimation of the total parasite biomass in acute falciparum malaria from plasma PfHRP2. PLoS Med 2005;2:e204.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.0020204&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16104831&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 

18. 18.Fukuda N, Balikagala B, Ueno T, Anywar DA, Kimura E, et al. The Impact of Sequestration on Artemisinin-Induced Parasite Clearance in Plasmodium falciparum Malaria in Africa. Clin Infect Dis Off Publ Infect Dis Soc Am 2023;76:1585–1593.
    
    

19. 19.Imwong M, Stepniewska K, Tripura R, Peto TJ, Lwin KM, et al. The numerical distributions of parasite densities in asymptomatic malaria. J Infect Dis.
    
    

20. 20.Lamikanra AA, Dobaño C, Jiménez A, Nhabomba A, Tsang HP, et al. A direct comparison of real time PCR on plasma and blood to detect Plasmodium falciparum infection in children. Malar J 2012;11:201.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1475-2875-11-201&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22704637&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 

21. 21.Sambe BS, Diagne A, Diatta HAM, Gaba FM, Sarr I, et al. Molecular detection and quantification of Plasmodium vivax DNA in blood pellet and plasma samples from patients in Senegal. Front Parasitol;2. Epub ahead of print 24 April 2023. DOI: 10.3389/fpara.2023.1149738.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fpara.2023.1149738&link_type=DOI) 

22. 22.Ruybal-Pesántez S, McCann K, Vibin J, Siegel S, Auburn S, et al. Molecular markers for malaria genetic epidemiology: progress and pitfalls. Trends Parasitol 2024;40:147–163.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pt.2023.11.006&link_type=DOI) 

23. 23.Manske M, Miotto O, Campino S, Auburn S, Almagro-Garcia J, et al. Analysis of Plasmodium falciparum diversity in natural infections by deep sequencing - Nature. Nature 2012;487:375–379.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11174&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22722859&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000306506500045&link_type=ISI) 

24. 24.Chang H-H, Worby CJ, Yeka A, Nankabirwa J, Kamya MR, et al. THE REAL McCOIL: A method for the concurrent estimation of the complexity of infection and SNP allele frequency for malaria parasites. PLoS Comput Biol 2017;13:e1005348.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1005348&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28125584&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 

25. 25.Pearson RD, Amato R, Auburn S, Miotto O, Almagro-Garcia J, et al. Genomic analysis of local variation and recent evolution in Plasmodium vivax. Nat Genet 2016;48:959–964.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3599&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27348299&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 

26. 26.Lerch A, Koepfli C, Hofmann NE, Messerli C, Wilcox S, et al. Development of amplicon deep sequencing markers and data analysis pipeline for genotyping multi-clonal malaria infections. BMC Genomics 2017;18:864.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12864-017-4260-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29132317&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 

27. 27.Kamaliddin C, Rombaut D, Guillochon E, Royo J, Ezinmegnon S, et al. From genomic to LC-MS/MS evidence: Analysis of PfEMP1 in Benin malaria cases. PLOS ONE 2019;14:e0218012.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0218012&link_type=DOI) 

28. 28.MalariaGEN, Ahouidi A, Ali M, Almagro-Garcia J, Amambua-Ngwa A, et al. An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples. Wellcome Open Res 2021;6:42.
    
    

29. 29.Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, et al. Twelve years of SAMtools and BCFtools. GigaScience 2021;10:giab008.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gigascience/giab008&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33590861&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 

30. 30.Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinforma Oxf Engl 2010;26:841–842.
    
    

31. 31.Gruenberg M, Lerch A, Beck H-P, Felger I. Amplicon deep sequencing improves Plasmodium falciparum genotyping in clinical trials of antimalarial drugs. Sci Rep 2019;9:17790.
    
    

32. 32.R Core Team. R: A Language and Environment for Statistical Computing.
    
    

33. 33.White NJ. Malaria parasite clearance. Malar J 2017;16:88.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-017-1731-1&link_type=DOI) 

34. 34.Jones S, Kay K, Hodel EM, Gruenberg M, Lerch A, et al. Should Deep-Sequenced Amplicons Become the New Gold Standard for Analyzing Malaria Drug Clinical Trials? Antimicrob Agents Chemother;65. Epub ahead of print October 2021. DOI: 10.1128/AAC.00437-21.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/AAC.00437-21&link_type=DOI) 

35. 35.Wilson DW, Langer C, Goodman CD, McFadden GI, Beeson JG. Defining the Timing of Action of Antimalarial Drugs against Plasmodium falciparum. Antimicrob Agents Chemother 2013;57:1455–1467.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjk6IjU3LzMvMTQ1NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzIwLzIwMjQuMDguMjAuMjQzMTIyOTYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

36. 36.Chotivanich K, Udomsangpetch R, McGready R, Proux S, Newton P, et al. Central role of the spleen in malaria parasite clearance. J Infect Dis 2002;185:1538–1541.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/340213&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11992295&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000175314700025&link_type=ISI) 

37. 37.Udomsangpetch R, Pipitaporn B, Krishna S, Angus B, Pukrittayakamee S, et al. Antimalarial drugs reduce cytoadherence and rosetting Plasmodium falciparum. J Infect Dis 1996;173:691–698.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/173.3.691&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8627034&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996TW80600023&link_type=ISI) 

38. 38.Paris L, Tackie RG, Beshir KB, Tampuori J, Awandare GA, et al. Parasite clearance dynamics in children hospitalised with severe malaria in the Ho Teaching Hospital, Volta Region, Ghana. Parasite Epidemiol Control 2022;19:e00276.
    
    

39. 39.Mideo N, Bailey JA, Hathaway NJ, Ngasala B, Saunders DL, et al. A deep sequencing tool for partitioning clearance rates following antimalarial treatment in polyclonal infections. Evol Med Public Health 2016;2016:21–36.
    
    

40. 40.Wamae K, Kimenyi KM, Osoti V, Laurent ZR de, Ndwiga L, et al. Amplicon Sequencing as a Potential Surveillance Tool for Complexity of Infection and Drug Resistance Markers in Plasmodium falciparum Asymptomatic Infections. J Infect Dis 2022;226:920.
    
    

41. 41.Marks F, von Kalckreuth V, Kobbe R, Adjei S, Adjei O, et al. Parasitological rebound effect and emergence of pyrimethamine resistance in Plasmodium falciparum after single-dose sulfadoxine-pyrimethamine. J Infect Dis 2005;192:1962–1965.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/497698&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16267768&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233018200015&link_type=ISI) 

42. 42.White NJ. The parasite clearance curve. Malar J 2011;10:278.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1475-2875-10-278&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21939506&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom) 

43. 43.Flegg JA, Guerin PJ, White NJ, Stepniewska K. Standardizing the measurement of parasite clearance in falciparum malaria: the parasite clearance estimator. Malar J 2011;10:339.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1475-2875-10-339&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22074219&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F20%2F2024.08.20.24312296.atom)