Publications

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.

Aquatic snails, the intermediate hosts of schistosomes, harbor a diverse unexplored microbiome. We speculate that this may play a critical role in host-parasite interactions. We summarize our current knowledge of snail microbiomes and highlight future research priorities.

Classical malaria parasite genetic crosses involve isolation, genotyping, and phenotyping of progeny parasites, which is time consuming and laborious. We tested a rapid alternative approach-bulk segregant analysis (BSA)-that utilizes sequencing of bulk progeny populations with and without drug selection for rapid identification of drug resistance loci. We used dihydroartemisinin (DHA) selection in two genetic crosses and investigated how synchronization, cryopreservation, and the drug selection regimen impacted BSA success. We detected a robust quantitative trait locus (QTL) at kelch13 in both crosses but did not detect QTLs at four other candidate loci. QTLs were detected using synchronized, but not unsynchronized progeny pools, consistent with the stage-specific action of DHA. We also successfully applied BSA to cryopreserved progeny pools, expanding the utility of this approach. We conclude that BSA provides a powerful approach for investigating the genetic architecture of drug resistance in Plasmodium falciparum.

Schistosoma mansoni, a snail-borne, blood fluke that infects humans, was introduced into the Americas from Africa during the Trans-Atlantic slave trade. As this parasite shows strong specificity to the snail intermediate host, we expected that adaptation to South American Biomphalaria spp. snails would result in population bottlenecks and strong signatures of selection. We scored 475,081 single nucleotide variants in 143 S. mansoni from the Americas (Brazil, Guadeloupe and Puerto Rico) and Africa (Cameroon, Niger, Senegal, Tanzania, and Uganda), and used these data to ask: (i) Was there a population bottleneck during colonization? (ii) Can we identify signatures of selection associated with colonization? (iii) What were the source populations for colonizing parasites? We found a 2.4- to 2.9-fold reduction in diversity and much slower decay in linkage disequilibrium (LD) in parasites from East to West Africa. However, we observed similar nuclear diversity and LD in West Africa and Brazil, suggesting no strong bottlenecks and limited barriers to colonization. We identified five genome regions showing selection in the Americas, compared with three in West Africa and none in East Africa, which we speculate may reflect adaptation during colonization. Finally, we infer that unsampled populations from central African regions between Benin and Angola, with contributions from Niger, are probably the major source(s) for Brazilian S. mansoni. The absence of a bottleneck suggests that this is a rare case of a serendipitous invasion, where S. mansoni parasites were pre-adapted to the Americas and able to establish with relative ease.

Both theory and experimental data from pathogens suggest that the production of transmission stages should be strongly associated with virulence, but the genetic bases of parasite transmission/virulence traits are poorly understood. The blood fluke Schistosoma mansoni shows extensive variation in numbers of cercariae larvae shed and in their virulence to infected snail hosts, consistent with expected trade-offs between parasite transmission and virulence. We crossed schistosomes from two populations that differ 8-fold in cercarial shedding and in their virulence to Biomphalaria glabrata snail hosts, and determined four-week cercarial shedding profiles in F0 parents, F1 parents and 376 F2 progeny from two independent crosses in inbred snails. Sequencing and linkage analysis revealed that cercarial production is polygenic and controlled by five QTLs (i.e. Quantitative Trait Loci). These QTLs act additively, explaining 28.56% of the phenotypic variation. These results demonstrate that the genetic architecture of key traits relevant to schistosome ecology can be dissected using classical linkage mapping approaches.