Malaria parasites break down host haemoglobin into peptides and amino acids in the digestive vacuole for export to the parasite cytoplasm for growth: interrupting this process is central to the mode of action of several antimalarial drugs. Mutations in the chloroquine (CQ) resistance transporter, pfcrt, located in the digestive vacuole membrane, confer CQ resistance in Plasmodium falciparum, and typically also affect parasite fitness. However, the role of other parasite loci in the evolution of CQ resistance is unclear. Here we use a combination of population genomics, genetic crosses and gene editing to demonstrate that a second vacuolar transporter plays a key role in both resistance and compensatory evolution. Longitudinal genomic analyses of the Gambian parasites revealed temporal signatures of selection on a putative amino acid transporter (pfaat1) variant S258L, which increased from 0% to 97% in frequency between 1984 and 2014 in parallel with the pfcrt1 K76T variant. Parasite genetic crosses then identified a chromosome 6 quantitative trait locus containing pfaat1 that is selected by CQ treatment. Gene editing demonstrated that pfaat1 S258L potentiates CQ resistance but at a cost of reduced fitness, while pfaat1 F313S, a common southeast Asian polymorphism, reduces CQ resistance while restoring fitness. Our analyses reveal hidden complexity in CQ resistance evolution, suggesting that pfaat1 may underlie regional differences in the dynamics of resistance evolution, and modulate parasite resistance or fitness by manipulating the balance between both amino acid and drug transport.
Publications
2023
BACKGROUND: The trematode parasite Schistosoma mansoni uses an aquatic snail intermediate and a vertebrate definitive host to complete its life cycle. We previously showed that a key transmission trait-the number of cercariae larvae shed from infected Biomphalaria spp. snails-varies significantly within and between different parasite populations and is genetically controlled by five loci. We investigated the hypothesis that the success of parasite genotypes showing high propagative fitness in the intermediate snail host may be offset by lower reproductive fitness in the definitive vertebrate host.
METHODS: We investigated this trade-off hypothesis by selecting parasite progeny producing high or low number of larvae in the snail and then comparing fitness parameters and virulence in the rodent host. We infected inbred BALB/c mice using two Schistosoma mansoni parasite lines [high shedder (HS) and low shedder (LS) lines] isolated from F2 progeny generated by genetic crosses between SmLE (HS parent) and SmBRE (LS parent) parasites. We used the F3 progeny to infect two populations of inbred Biomphalaria glabrata snails. We then compared life history traits and virulence of these two selected parasite lines in the rodent host to understand pleiotropic effects of genes determining cercarial shedding in parasites infecting the definitive host.
RESULTS: HS parasites shed high numbers of cercariae, which had a detrimental impact on snail physiology (measured by laccase-like activity and hemoglobin rate), regardless of the snail genetic background. In contrast, selected LS parasites shed fewer cercariae and had a lower impact on snail physiology. Similarly, HS worms have a higher reproductive fitness and produced more viable F3 miracidia larvae than LS parasites. This increase in transmission is correlated with an increase in virulence toward the rodent host, characterized by stronger hepato-splenomegaly and hepatic fibrosis.
CONCLUSIONS: These experiments revealed that schistosome parasite propagative and reproductive fitness was positively correlated in intermediate and definitive host (positive pleiotropy). Therefore, we rejected our trade-off hypothesis. We also showed that our selected schistosome lines exhibited low and high shedding phenotype regardless of the intermediate snail host genetic background. .
We describe the MalariaGEN Pf7 data resource, the seventh release of Plasmodium falciparum genome variation data from the MalariaGEN network. It comprises over 20,000 samples from 82 partner studies in 33 countries, including several malaria endemic regions that were previously underrepresented. For the first time we include dried blood spot samples that were sequenced after selective whole genome amplification, necessitating new methods to genotype copy number variations. We identify a large number of newly emerging crt mutations in parts of Southeast Asia, and show examples of heterogeneities in patterns of drug resistance within Africa and within the Indian subcontinent. We describe the profile of variations in the C-terminal of the csp gene and relate this to the sequence used in the RTS,S and R21 malaria vaccines. Pf7 provides high-quality data on genotype calls for 6 million SNPs and short indels, analysis of large deletions that cause failure of rapid diagnostic tests, and systematic characterisation of six major drug resistance loci, all of which can be freely downloaded from the MalariaGEN website.
2022
The human malaria parasite Plasmodium falciparum is globally widespread, but its prevalence varies significantly between and even within countries. Most population genetic studies in P. falciparum focus on regions of high transmission where parasite populations are large and genetically diverse, such as sub-Saharan Africa. Understanding population dynamics in low transmission settings, however, is of particular importance as these are often where drug resistance first evolves. Here, we use the Pacific Coast of Colombia and Ecuador as a model for understanding the population structure and evolution of Plasmodium parasites in small populations harboring less genetic diversity. The combination of low transmission and a high proportion of monoclonal infections means there are few outcrossing events and clonal lineages persist for long periods of time. Yet despite this, the population is evolutionarily labile and has successfully adapted to changes in drug regime. Using newly sequenced whole genomes, we measure relatedness between 166 parasites, calculated as identity by descent (IBD), and find 17 distinct but highly related clonal lineages, six of which have persisted in the region for at least a decade. This inbred population structure is captured in more detail with IBD than with other common population structure analyses like PCA, ADMIXTURE, and distance-based trees. We additionally use patterns of intra-chromosomal IBD and an analysis of haplotypic variation to explore past selection events in the region. Two genes associated with chloroquine resistance, crt and aat1, show evidence of hard selective sweeps, while selection appears soft and/or incomplete at three other key resistance loci (dhps, mdr1, and dhfr). Overall, this work highlights the strength of IBD analyses for studying parasite population structure and resistance evolution in regions of low transmission, and emphasizes that drug resistance can evolve and spread in small populations, as will occur in any region nearing malaria elimination.
In 2005 Richard Carter's group surprised the malaria genetics community with an elegant approach to rapidly mapping the genetic basis of phenotypic traits in rodent malaria parasites. This approach, which he termed "linkage group selection", utilized bulk pools of progeny, rather than individual clones, and exploited simple selection schemes to identify genome regions underlying resistance to drug treatment (or other phenotypes). This work was the first application of "bulk segregant" methodologies for genetic mapping in microbes: this approach is now widely used in yeast, and across multiple recombining pathogens ranging from Aspergillus fungi to Schistosome parasites. Genetic crosses of human malaria parasites (for which Richard Carter was also a pioneer) can now be conducted in humanized mice, providing new opportunities for exploiting bulk segregant approaches for a wide variety of malaria parasite traits. We review the application of bulk segregant approaches to mapping malaria parasite traits and suggest additional developments that may further expand the utility of this powerful approach.
What genes determine in vitro growth and nutrient utilization in asexual blood-stage malaria parasites? Competition experiments between NF54, clone 3D7, a lab-adapted African parasite, and a recently isolated Asian parasite (NHP4026) reveal contrasting outcomes in different media: 3D7 outcompetes NHP4026 in media containing human serum, while NHP4026 outcompetes 3D7 in media containing AlbuMAX, a commercial lipid-rich bovine serum formulation. To determine the basis for this polymorphism, we conducted parasite genetic crosses using humanized mice and compared genome-wide allele frequency changes in three independent progeny populations cultured in media containing human serum or AlbuMAX. This bulk segregant analysis detected three quantitative trait loci (QTL) regions [on chromosome (chr) 2 containing aspartate transaminase AST; chr 13 containing EBA-140; and chr 14 containing cysteine protease ATG4] linked with differential growth in serum or AlbuMAX in each of the three independent progeny pools. Selection driving differential growth was strong (s = 0.10 - 0.23 per 48-hour lifecycle). We conducted validation experiments for the strongest QTL on chr 13: competition experiments between ΔEBA-140 and 3D7 wildtype parasites showed fitness reversals in the two medium types as seen in the parental parasites, validating this locus as the causative gene. These results (i) demonstrate the effectiveness of bulk segregant analysis for dissecting fitness traits in P. falciparum genetic crosses, and (ii) reveal intimate links between red blood cell invasion and nutrient composition of growth media. Use of parasite crosses combined with bulk segregant analysis will allow systematic dissection of key nutrient acquisition/metabolism and red blood cell invasion pathways in P. falciparum.
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).
Oxamniquine (OXA) is a prodrug activated by a sulfotransferase (SULT) that was only active against Schistosoma mansoni. We have reengineered OXA to be effective against S. haematobium and S. japonicum. Three derivatives stand out, CIDD-0066790, CIDD-0072229, and CIDD-0149830 as they kill all three major human schistosome species. However, questions remain. Is the OXA mode of action conserved in derivatives? RNA-interference experiments demonstrate that knockdown of the SmSULT, ShSULT, and SjSULT results in resistance to CIDD-0066790. Confirming that the OXA-derivative mode of action is conserved. Next is the level of expression of the schistosome SULTs in each species, as well as changes in SULT expression throughout development in S. mansoni. Using multiple tools, our data show that SmSULT has higher expression compared to ShSULT and SjSULT. Third, is the localization of SULT in the adult, multicellular eucaryotic schistosome species. We utilized fluorescence in situ hybridization and uptake of radiolabeled OXA to determine that multiple cell types throughout the adult schistosome worm express SULT. Thus, we hypothesize the ability of many cells to express the sulfotransferase accounts for the ability of the OXA derivatives to kill adult worms. Our studies demonstrate that the OXA derivatives are able to kill all three human schistosome species and thus will be a useful complement to PZQ.
Drug resistance mutations tend to disrupt key physiological processes and frequently carry fitness costs, which are a central determinant of the rate of spread of these mutations in natural populations. Head-to-head competition assays provide a standard approach to measuring fitness for malaria parasites. These assays typically use a standardized culture medium containing RPMI 1640, which has a 1.4- to 5.5-fold higher concentration of amino acids than human blood. In this rich medium, we predict that fitness costs will be underestimated because resource competition is weak. We tested this prediction using an artemisinin-sensitive parasite edited to contain kelch-C580Y or R561H mutations conferring resistance to artemisinin or synonymous control mutations. We examined the impact of these single amino acid mutations on fitness, using replicated head-to-head competition experiments conducted in media containing (i) normal RPMI, (ii) modified RPMI with reduced amino acid concentration, (iii) RPMI containing only isoleucine, or (iv) 3-fold diluted RPMI. We found a significant 1.3- to 1.4-fold increase in fitness costs measured in modified and isoleucine-only media relative to normal media, while fitness costs were 2.5-fold higher in diluted media. We conclude that fitness costs are strongly affected by media composition and will be significantly underestimated in normal RPMI. Several components differed between media, including pABA and sodium bicarbonate concentrations, so we cannot directly determine which is responsible. Elevated fitness costs in nature will limit spread of artemisinin (ART) resistance but will also promote evolution of compensatory mutations that restore fitness and can be exploited to maximize selection in laboratory experiments.