Search

Background: Rapid diagnostic tests based on detection of histidine-rich proteins (HRPs) are widely used for malaria diagnosis, but parasites carrying pfhrp deletions can evade detection and are increasing in frequency in some countries. Models aim to predict conditions under which pfhrp2 and/or pfhrp3 deletions will increase, but a key parameter-the fitness cost of deletions-is unknown., Methods: We removed pfhrp2 and/or pfhrp3 from a Malawian parasite clone using gene editing approaches) and measured fitness costs by conducting pairwise competition experiments., Results: We observed significant fitness costs of 0.087 ± 0.008 (1 standard error) per asexual cycle for pfhrp2 deletion and 0.113 ± 0.008 for the pfhrp2/3 double deletion, relative to the unedited progenitor parasite. Selection against deletions is strong and comparable to that resulting from drug resistance mutations., Conclusions: Prior modeling suggested that diagnostic selection may drive increased frequency of pfhrp deletions only when fitness costs are mild. Our experiments show that costs of pfhrp deletions are higher than these thresholds, but modeling and empirical results can be reconciled if the duration of infection is short. These results may inform future modeling to understand why pfhrp2/3 deletions are increasing in some locations (Ethiopia and Eritrea) but not in others (Mekong region)., Competing Interests: Potential conflicts of interest. All authors: No reported potential conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed, (© The Author(s) 2022. Published by Oxford University Press on behalf of Infectious Diseases Society of America. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Malaria parasites have small extremely AT-rich genomes: microsatellite repeats (1-9 bp) comprise 11% of the genome and genetic variation in natural populations is dominated by repeat changes in microsatellites rather than point mutations. This experiment was designed to quantify microsatellite mutation patterns in Plasmodium falciparum. We established 31 parasite cultures derived from a single parasite cell and maintained these for 114-267 days with frequent reductions to a single cell, so parasites accumulated mutations during ∼13,207 cell divisions. We Illumina sequenced the genomes of both progenitor and end-point mutation accumulation (MA) parasite lines in duplicate to validate stringent calling parameters. Microsatellite calls were 99.89% (GATK), 99.99% (freeBayes), and 99.96% (HipSTR) concordant in duplicate sequence runs from independent sequence libraries, whereas introduction of microsatellite mutations into the reference genome revealed a low false negative calling rate (0.68%). We observed 98 microsatellite mutations. We highlight several conclusions: microsatellite mutation rates (3.12 × 10-7 to 2.16 × 10-8/cell division) are associated with both repeat number and repeat motif like other organisms studied. However, 41% of changes resulted from loss or gain of more than one repeat: this was particularly true for long repeat arrays. Unlike other eukaryotes, we found no insertions or deletions that were not associated with repeats or homology regions. Overall, microsatellite mutation rates are among the lowest recorded and comparable to those in another AT-rich protozoan (Dictyostelium). However, a single infection (>1011 parasites) will still contain over 2.16 × 103 to 3.12 × 104 independent mutations at any single microsatellite locus., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Single-cell genomics is a powerful tool for determining the genetic architecture of complex communities of unicellular organisms. In areas of high transmission, malaria patients are often challenged by the activities of multiple Plasmodium falciparum lineages, which can potentiate pathology, spread drug resistance loci, and also complicate most genetic analysis. Single-cell sequencing of P. falciparum would be key to understanding infection complexity, though efforts are hampered by the extreme nucleotide composition of its genome (∼80% AT-rich). To counter the low coverage achieved in previous studies, we targeted DNA-rich late-stage parasites by Fluorescence-Activated Cell Sorting and whole genome sequencing. Our method routinely generates accurate, near-complete capture of the 23 Mb P. falciparum genome (mean breadth of coverage 90.7%) at high efficiency. Data from 48 single-cell genomes derived from a polyclonal infection sampled in Chikhwawa, Malawi allowed for unambiguous determination of haplotype diversity and recent meiotic events, information that will aid public health efforts., (© The Author(s) 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Reductions in malaria transmission decrease naturally acquired immunity, which may influence the emergence of Plasmodium falciparum artemisinin-resistant phenotypes and genotypes over time., Methods: Antibodies specific for P. falciparum antigens were determined in uncomplicated hyperparasitemic malaria patients over a 10-year period of declining malaria transmission and emerging artemisinin resistance in northwestern Thailand. We investigated the association between antibody levels and both parasite clearance time (PCt½) and artemisinin resistance-associated kelch13 genotypes over time., Results: Immunity to P. falciparum declined prior to 2004, preceding the emergence of artemisinin resistance-associated genotypes and phenotypes (maximum mean change in antibody level per year: anti-MSP142 = -0.17; 95% confidence interval [CI] = -.31 to -.04; P = .01). In this period of declining immunity, and in the absence of kelch13 mutations, PCt½ increased. Between 2007 and 2011, levels of antibodies fluctuated, and higher antibody levels were associated with faster PCt½ (maximum yearly change in PCt½, in hours: EBA140rII = -0.39; 95% CI = -.61 to -.17; P < .001)., Conclusions: Understanding the impact of changing transmission and immunity on the emergence of artemisinin resistance is important particularly as increased malaria control and elimination activities may enhance immunological conditions for the expansion of artemisinin-resistant P. falciparum., (© The Author 2017. Published by Oxford University Press for the Infectious Diseases Society of America.)

Multiple kelch13 alleles conferring artemisinin resistance (ART-R) are currently spreading through Southeast Asian malaria parasite populations, providing a unique opportunity to observe an ongoing soft selective sweep, investigate why resistance alleles have evolved multiple times and determine fundamental population genetic parameters for Plasmodium We sequenced kelch13 (n = 1,876), genotyped 75 flanking SNPs, and measured clearance rate (n = 3,552) in parasite infections from Western Thailand (2001-2014). We describe 32 independent coding mutations including common mutations outside the kelch13 propeller associated with significant reductions in clearance rate. Mutations were first observed in 2003 and rose to 90% by 2014, consistent with a selection coefficient of ∼0.079. ART-R allele diversity rose until 2012 and then dropped as one allele (C580Y) spread to high frequency. The frequency with which adaptive alleles arise is determined by the rate of mutation and the population size. Two factors drive this soft sweep: (1) multiple kelch13 amino-acid mutations confer resistance providing a large mutational target-we estimate the target is 87-163 bp. (2) The population mutation parameter (Θ = 2N e μ) can be estimated from the frequency distribution of ART-R alleles and is ∼5.69, suggesting that short term effective population size is 88 thousand to 1.2 million. This is 52-705 times greater than N e estimated from fluctuation in allele frequencies, suggesting that we have previously underestimated the capacity for adaptive evolution in Plasmodium Our central conclusions are that retrospective studies may underestimate the complexity of selective events and the N e relevant for adaptation for malaria is considerably higher than previously estimated., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Deployment of mefloquine-artesunate (MAS3) on the Thailand-Myanmar border has led to a sustained reduction in falciparum malaria, although antimalarial efficacy has declined substantially in recent years. The role of Plasmodium falciparum K13 mutations (a marker of artemisinin resistance) in reducing treatment efficacy remains controversial., Methods: Between 2003 and 2013, we studied the efficacy of MAS3 in 1005 patients with uncomplicated P. falciparum malaria in relation to molecular markers of resistance., Results: Polymerase chain reaction (PCR)-adjusted cure rates declined from 100% in 2003 to 81.1% in 2013 as the proportions of isolates with multiple Pfmdr1 copies doubled from 32.4% to 64.7% and those with K13 mutations increased from 6.7% to 83.4%. K13 mutations conferring moderate artemisinin resistance (notably E252Q) predominated initially but were later overtaken by propeller mutations associated with slower parasite clearance (notably C580Y). Those infected with both multiple Pfmdr1 copy number and a K13 propeller mutation were 14 times more likely to fail treatment. The PCR-adjusted cure rate was 57.8% (95% confidence interval [CI], 45.4, 68.3) compared with 97.8% (95% CI, 93.3, 99.3) in patients with K13 wild type and Pfmdr1 single copy. K13 propeller mutation alone was a strong risk factor for recrudescence (P = .009). The combined population attributable fraction of recrudescence associated with K13 mutation and Pfmdr1 amplification was 82%., Conclusions: The increasing prevalence of K13 mutations was the decisive factor for the recent and rapid decline in efficacy of artemisinin-based combination (MAS3) on the Thailand-Myanmar border., (© The Author 2016. Published by Oxford University Press for the Infectious Diseases Society of America.)

If copy number variants (CNVs) are predominantly deleterious, we would expect them to be more efficiently purged from populations with a large effective population size (Ne) than from populations with a small Ne. Malaria parasites (Plasmodium falciparum) provide an excellent organism to examine this prediction, because this protozoan shows a broad spectrum of population structures within a single species, with large, stable, outbred populations in Africa, small unstable inbred populations in South America and with intermediate population characteristics in South East Asia. We characterized 122 single-clone parasites, without prior laboratory culture, from malaria-infected patients in seven countries in Africa, South East Asia and South America using a high-density single-nucleotide polymorphism/CNV microarray. We scored 134 high-confidence CNVs across the parasite exome, including 33 deletions and 102 amplifications, which ranged in size from <500 bp to 59 kb, as well as 10,107 flanking, biallelic single-nucleotide polymorphisms. Overall, CNVs were rare, small, and skewed toward low frequency variants, consistent with the deleterious model. Relative to African and South East Asian populations, CNVs were significantly more common in South America, showed significantly less skew in allele frequencies, and were significantly larger. On this background of low frequency CNV, we also identified several high-frequency CNVs under putative positive selection using an FST outlier analysis. These included known adaptive CNVs containing rh2b and pfmdr1, and several other CNVs (e.g., DNA helicase and three conserved proteins) that require further investigation. Our data are consistent with a significant impact of genetic structure on CNV burden in an important human pathogen., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Accurate measurement of malaria parasite clearance rates (CRs) following artemisinin (ART) treatment is critical for resistance surveillance and research, and various CR metrics are currently used. We measured 13 CR metrics in 1472 ART-treated hyperparasitemia infections for which 6-hour parasite counts and parasite genotypes (93 single nucleotide polymorphisms [SNPs]) were available. We used heritability to evaluate the performance of each metric. Heritability ranged from 0.06 ± 0.06 (SD) for 50% parasite clearance times to 0.67 ± 0.04 (SD) for clearance half-lives estimated from 6-hour parasite counts. These results identify the measures that should be avoided and show that reliable clearance measures can be obtained with abbreviated monitoring protocols.

In western Cambodia, malaria parasites clear slowly from the blood after treatment with artemisinin derivatives, but it is unclear whether this results from parasite, host, or other factors specific to this population. We measured heritability of clearance rate by evaluating patients infected with identical or nonidentical parasite genotypes, using methods analogous to human twin studies. A substantial proportion (56%-58%) of the variation in clearance rate is explained by parasite genetics. This has 2 important implications: (1) selection with artemisinin derivatives will tend to drive resistance spread and (2) because heritability is high, the genes underlying parasite clearance rate may be identified by genome-wide association.

Background: Relapses originating from hypnozoites are characteristic of Plasmodium vivax infections. Thus, reappearance of parasitemia after treatment can result from relapse, recrudescence, or reinfection. It has been assumed that parasites causing relapse would be a subset of the parasites that caused the primary infection., Methods: Paired samples were collected before initiation of antimalarial treatment and at recurrence of parasitemia from 149 patients with vivax malaria in Thailand (n=36), where reinfection could be excluded, and during field studies in Myanmar (n=75) and India (n=38)., Results: Combined genetic data from 2 genotyping approaches showed that novel P. vivax populations were present in the majority of patients with recurrent infection (107 [72%] of 149 patients overall [78% of patients in Thailand, 75% of patients in Myanmar {Burma}, and 63% of patients in India]). In 61% of the Thai and Burmese patients and in 55% of the Indian patients, the recurrent infections contained none of the parasite genotypes that caused the acute infection., Conclusions: The P. vivax populations emerging from hypnozoites commonly differ from the populations that caused the acute episode. Activation of heterologous hypnozoite populations is the most common cause of first relapse in patients with vivax malaria.

When selection is strong and beneficial alleles have a single origin, local reductions in genetic diversity are expected. However, when beneficial alleles have multiple origins or were segregating in the population prior to a change in selection regime, the impact on genetic diversity may be less clear. We describe an example of such a "soft" selective sweep in the malaria parasite Plasmodium falciparum that involves adaptive genome rearrangements. Amplification in copy number of genome regions containing the pfmdr1 gene on chromosome 5 confer resistance to mefloquine and spread rapidly in the 1990s. Using flanking microsatellite data and real-time polymerase chain reaction determination of copy number, we show that 5-15 independent amplification events have occurred in parasites on the Thailand/Burma border. The amplified genome regions (amplicons) range in size from 14.7 to 49 kb and contain 2-11 genes, with 2-4 copies arranged in tandem. To examine the impact of drug selection on flanking variation, we genotyped 48 microsatellites on chromosome 5 in 326 parasites from a single Thai location. Diversity was reduced in a 170- to 250-kb (10-15 cM) region of chromosomes containing multiple copies of pfmdr1, consistent with hitchhiking resulting from the rapid recent spread of selected chromosomes. However, diversity immediately flanking pfmdr1 is reduced by only 42% on chromosomes bearing multiple amplicons relative to chromosomes carrying a single copy. We highlight 2 features of these results: 1) All amplicon break points occur in monomeric A/T tracts (9-45 bp). Given the abundance of these tracts in P. falciparum, we expect that duplications will occur frequently at multiple genomic locations and have been underestimated as drivers of phenotypic evolution in this pathogen. 2) The signature left by the spread of amplified genome segments is broad, but results in only limited reduction in diversity. If such "soft" sweeps are common in nature, statistical methods based on diversity reduction may be inefficient at detecting evidence for selection in genome-wide marker screens. This may be particularly likely when mutation rate is high, as appears to be the case for gene duplications, and in pathogen populations where effective population sizes are typically very large.

Background: Our study examined the relative contributions of host, pharmacokinetic, and parasitological factors in determining the therapeutic response to artemether-lumefantrine (AL)., Methods: On the northwest border of Thailand, patients with uncomplicated Plasmodium falciparum malaria were enrolled in prospective studies of AL treatment (4- or 6-dose regimens) and followed up for 42 days. Plasma lumefantrine concentrations were measured by high performance liquid chromatography; malaria parasite pfmdr1 copy number was quantified using a real-time polymerase chain reaction assay (PCR), and in vitro drug susceptibility was tested., Results: All treatments resulted in a rapid clinical response and were well tolerated. PCR-corrected failure rates at day 42 were 13% (95% confidence interval [CI], 9.6%-17%) for the 4-dose regimen and 3.2% (95% CI, 1.8%-4.6%) for the 6-dose regimen. Increased pfmdr1 copy number was associated with a 2-fold (95% CI, 1.8-2.4-fold) increase in lumefantrine inhibitory concentration(50) (P=.001) and an adjusted hazard ratio for risk of treatment failure following completion of a 4-dose regimen, but not a 6-dose regimen, of 4.0 (95% CI, 1.4-11; P=.008). Patients who had lumefantrine levels below 175 ng/mL on day 7 were more likely to experience recrudescence by day 42 (adjusted hazard ratio, 17; 95% CI, 5.5-53), allowing prediction of treatment failure with 75% sensitivity and 84% specificity. The 6-dose regimen ensured that therapeutic levels were achieved in 91% of treated patients., Conclusions: The lumefantrine plasma concentration profile is the main determinant of efficacy of artemether-lumefantrine. Amplification in pfmdr1 determines lumefantrine susceptibility and, therefore, treatment responses when plasma lumefantrine levels are subtherapeutic.

Loci targeted by directional selection are expected to show elevated geographical population structure relative to neutral loci, and a flurry of recent papers have used this rationale to search for genome regions involved in adaptation. Studies of functional mutations that are known to be under selection are particularly useful for assessing the utility of this approach. Antimalarial drug treatment regimes vary considerably between countries in Southeast Asia selecting for local adaptation at parasite loci underlying resistance. We compared the population structure revealed by 10 nonsynonymous mutations (nonsynonymous single-nucleotide polymorphisms [nsSNPs]) in four loci that are known to be involved in antimalarial drug resistance, with patterns revealed by 10 synonymous mutations (synonymous single-nucleotide polymorphisms [sSNPs]) in housekeeping genes or genes of unknown function in 755 Plasmodium falciparum infections collected from 13 populations in six Southeast Asian countries. Allele frequencies at known nsSNPs underlying resistance varied markedly between locations (F(ST) = 0.18-0.66), with the highest frequencies on the Thailand-Burma border and the lowest frequencies in neighboring Lao PDR. In contrast, we found weak but significant geographic structure (F(ST) = 0-0.14) for 8 of 10 sSNPs. Importantly, all 10 nsSNPs showed significantly higher F(ST) (P < 8 x 10(-5)) than simulated neutral expectations based on observed F(ST) values in the putatively neutral sSNPs. This result was unaffected by the methods used to estimate allele frequencies or the number of populations used in the simulations. Given that dense single-nucleotide polymorphism (SNP) maps and rapid SNP assay methods are now available for P. falciparum, comparing genetic differentiation across the genome may provide a valuable aid to identifying parasite loci underlying local adaptation to drug treatment regimes or other selective forces. However, the high proportion of polymorphic sites that appear to be under balancing selection (or linked to selected sites) in the P. falciparum genome violates the central assumption that selected sites are rare, which complicates identification of outlier loci, and suggests that caution is needed when using this approach.

Malaria parasites (Plasmodium falciparum) provide an excellent system in which to study the genomic effects of strong selection in a recombining eukaryote because the rapid spread of resistance to multiple drugs during the last the past 50 years has been well documented, the full genome sequence and a microsatellite map are now available, and haplotype data can be easily generated. We examined microsatellite variation around the dihydrofolate reductase (dhfr) gene on chromosome 4 of P. falciparum. Point mutations in dhfr are known to be responsible for resistance to the antimalarial drug pyrimethamine, and resistance to this drug has spread rapidly in Southeast (SE) Asia after its introduction in 1970s. We genotyped 33 microsatellite markers distributed across chromosome 4 in 61 parasites from a location on the Thailand/Myanmar border. We observed minimal microsatellite length variation in a 12-kb (0.7-cM) region flanking the dhfr gene and diminished variation for approximately 100 kb (6 cM), indicative of a single origin of resistant alleles. Furthermore, we found the same or similar microsatellite haplotypes flanked resistant dhfr alleles sampled from 11 parasite populations in five SE Asian countries indicating recent invasion of a single lineage of resistant dhfr alleles in locations 2000 km apart. Three features of these data are of especially interest. (1). Pyrimethamine resistance is generally assumed to have evolved multiple times because the genetic basis is simple and resistance can be selected easily in the laboratory. Yet our data clearly indicate a single origin of resistant dhfr alleles sampled over a large region of SE Asia. (2). The wide valley ( approximately 6 cM) of reduced variation around dhfr provides "proof-of-principle" that genome-wide association may be an effective way to locate genes under strong recent selection. (3). The width of the selective valley is consistent with predictions based on independent measures of recombination, mutation, and selection intensity, suggesting that we have reasonable estimates of these parameters. We conclude that scanning the malaria parasite genome for evidence of recent selection may prove an extremely effective way to locate genes underlying recently evolved traits such as drug resistance, as well as providing an opportunity to study the dynamics of selective events that have occurred recently or are currently in progress.