Red blood cell tension protects against severe malaria in the Dantu blood group.

Malaria has had a major effect on the human genome, with many protective polymorphisms-such as the sickle-cell trait-having been selected to high frequencies in malaria-endemic regions1,2. The blood group variant Dantu provides 74% protection against all forms of severe malaria in homozygous individuals3-5, a similar degree of protection to that afforded by the sickle-cell trait and considerably greater than that offered by the best malaria vaccine. Until now, however, the protective mechanism has been unknown. Here we demonstrate the effect of Dantu on the ability of the merozoite form of the malaria parasite Plasmodium falciparum to invade red blood cells (RBCs). We find that Dantu is associated with extensive changes to the repertoire of proteins found on the RBC surface, but, unexpectedly, inhibition of invasion does not correlate with specific RBC-parasite receptor-ligand interactions. By following invasion using video microscopy, we find a strong link between RBC tension and merozoite invasion, and identify a tension threshold above which invasion rarely occurs, even in non-Dantu RBCs. Dantu RBCs have higher average tension than non-Dantu RBCs, meaning that a greater proportion resist invasion. These findings provide both an explanation for the protective effect of Dantu, and fresh insight into why the efficiency of P. falciparum invasion might vary across the heterogenous populations of RBCs found both within and between individuals.

Declines in prevalence alter the optimal level of sexual investment for the malaria parasite Plasmodium falciparum.

Successful infectious disease interventions can result in large reductions in parasite prevalence. Such demographic change has fitness implications for individual parasites and may shift the parasite's optimal life history strategy. Here, we explore whether declining infection rates can alter Plasmodium falciparum's investment in sexual versus asexual growth. Using a multiscale mathematical model, we demonstrate how the proportion of polyclonal infections, which decreases as parasite prevalence declines, affects the optimal sexual development strategy: Within-host competition in multiclone infections favors a greater investment in asexual growth whereas single-clone infections benefit from higher conversion to sexual forms. At the same time, drug treatment also imposes selection pressure on sexual development by shortening infection length and reducing within-host competition. We assess these models using 148 P. falciparum parasite genomes sampled in French Guiana over an 18-y period of intensive intervention (1998 to 2015). During this time frame, multiple public health measures, including the introduction of new drugs and expanded rapid diagnostic testing, were implemented, reducing P. falciparum malaria cases by an order of magnitude. Consistent with this prevalence decline, we see an increase in the relatedness among parasites, but no single clonal background grew to dominate the population. Analyzing individual allele frequency trajectories, we identify genes that likely experienced selective sweeps. Supporting our model predictions, genes showing the strongest signatures of selection include transcription factors involved in the development of P. falciparum's sexual gametocyte form. These results highlight how public health interventions impose wide-ranging selection pressures that affect basic parasite life history traits.

Resolving drug selection and migration in an inbred South American Plasmodium falciparum population with identity-by-descent analysis.

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.

Capturing in vivo RNA transcriptional dynamics from the malaria parasite Plasmodium falciparum.

To capture the transcriptional dynamics within proliferating cells, methods to differentiate nascent transcription from preexisting mRNAs are desired. One approach is to label newly synthesized mRNA transcripts in vivo through the incorporation of modified pyrimidines. However, the human malaria parasite, Plasmodium falciparum, is incapable of pyrimidine salvage for mRNA biogenesis. To capture cellular mRNA dynamics during Plasmodium development, we engineered parasites that can salvage pyrimidines through the expression of a single bifunctional yeast fusion gene, cytosine deaminase/uracil phosphoribosyltransferase (FCU). We show that expression of FCU allows for the direct incorporation of thiol-modified pyrimidines into nascent mRNAs. Using developmental stage-specific promoters to express FCU-GFP enables the biosynthetic capture and in-depth analysis of mRNA dynamics from subpopulations of cells undergoing differentiation. We demonstrate the utility of this method by examining the transcriptional dynamics of the sexual gametocyte stage transition, a process that is essential to malaria transmission between hosts. Using the pfs16 gametocyte-specific promoter to express FCU-GFP in 3D7 parasites, we found that sexual stage commitment is governed by transcriptional reprogramming and stabilization of a subset of essential gametocyte transcripts. We also measured mRNA dynamics in F12 gametocyte-deficient parasites and demonstrate that the transcriptional program required for sexual commitment and maturation is initiated but likely aborted due to the absence of the PfAP2-G transcriptional regulator and a lack of gametocyte-specific mRNA stabilization. Biosynthetic labeling of Plasmodium mRNAs is incredibly versatile, can be used to measure transcriptional dynamics at any stage of parasite development, and will allow for future applications to comprehensively measure RNA-protein interactions in the malaria parasite.

Open Source Drug Discovery with the Malaria Box Compound Collection for Neglected Diseases and Beyond.

A major cause of the paucity of new starting points for drug discovery is the lack of interaction between academia and industry. Much of the global resource in biology is present in universities, whereas the focus of medicinal chemistry is still largely within industry. Open source drug discovery, with sharing of information, is clearly a first step towards overcoming this gap. But the interface could especially be bridged through a scale-up of open sharing of physical compounds, which would accelerate the finding of new starting points for drug discovery. The Medicines for Malaria Venture Malaria Box is a collection of over 400 compounds representing families of structures identified in phenotypic screens of pharmaceutical and academic libraries against the Plasmodium falciparum malaria parasite. The set has now been distributed to almost 200 research groups globally in the last two years, with the only stipulation that information from the screens is deposited in the public domain. This paper reports for the first time on 236 screens that have been carried out against the Malaria Box and compares these results with 55 assays that were previously published, in a format that allows a meta-analysis of the combined dataset. The combined biochemical and cellular assays presented here suggest mechanisms of action for 135 (34%) of the compounds active in killing multiple life-cycle stages of the malaria parasite, including asexual blood, liver, gametocyte, gametes and insect ookinete stages. In addition, many compounds demonstrated activity against other pathogens, showing hits in assays with 16 protozoa, 7 helminths, 9 bacterial and mycobacterial species, the dengue fever mosquito vector, and the NCI60 human cancer cell line panel of 60 human tumor cell lines. Toxicological, pharmacokinetic and metabolic properties were collected on all the compounds, assisting in the selection of the most promising candidates for murine proof-of-concept experiments and medicinal chemistry programs. The data for all of these assays are presented and analyzed to show how outstanding leads for many indications can be selected. These results reveal the immense potential for translating the dispersed expertise in biological assays involving human pathogens into drug discovery starting points, by providing open access to new families of molecules, and emphasize how a small additional investment made to help acquire and distribute compounds, and sharing the data, can catalyze drug discovery for dozens of different indications. Another lesson is that when multiple screens from different groups are run on the same library, results can be integrated quickly to select the most valuable starting points for subsequent medicinal chemistry efforts.

Metabolomic Profiling of the Malaria Box Reveals Antimalarial Target Pathways.

The threat of widespread drug resistance to frontline antimalarials has renewed the urgency for identifying inexpensive chemotherapeutic compounds that are effective against Plasmodium falciparum, the parasite species responsible for the greatest number of malaria-related deaths worldwide. To aid in the fight against malaria, a recent extensive screening campaign has generated thousands of lead compounds with low micromolar activity against blood stage parasites. A subset of these leads has been compiled by the Medicines for Malaria Venture (MMV) into a collection of structurally diverse compounds known as the MMV Malaria Box. Currently, little is known regarding the activity of these Malaria Box compounds on parasite metabolism during intraerythrocytic development, and a majority of the targets for these drugs have yet to be defined. Here we interrogated the in vitro metabolic effects of 189 drugs (including 169 of the drug-like compounds from the Malaria Box) using ultra-high-performance liquid chromatography-mass spectrometry (UHPLC-MS). The resulting metabolic fingerprints provide information on the parasite biochemical pathways affected by pharmacologic intervention and offer a critical blueprint for selecting and advancing lead compounds as next-generation antimalarial drugs. Our results reveal several major classes of metabolic disruption, which allow us to predict the mode of action (MoA) for many of the Malaria Box compounds. We anticipate that future combination therapies will be greatly informed by these results, allowing for the selection of appropriate drug combinations that simultaneously target multiple metabolic pathways, with the aim of eliminating malaria and forestalling the expansion of drug-resistant parasites in the field.

Flexible and cost-effective genomic surveillance of P. falciparum malaria with targeted nanopore sequencing.

Genomic surveillance of Plasmodium falciparum malaria can provide policy-relevant information about antimalarial drug resistance, diagnostic test failure, and the evolution of vaccine targets. Yet the large and low complexity genome of P. falciparum complicates the development of genomic methods, while resource constraints in malaria endemic regions can limit their deployment. Here, we demonstrate an approach for targeted nanopore sequencing of P. falciparum from dried blood spots (DBS) that enables cost-effective genomic surveillance of malaria in low-resource settings. We release software that facilitates flexible design of amplicon sequencing panels and use this software to design two target panels for P. falciparum. The panels generate 3-4 kbp reads for eight and sixteen targets respectively, covering key drug-resistance associated genes, diagnostic test antigens, polymorphic markers and the vaccine target csp. We validate our approach on mock and field samples, demonstrating robust sequencing coverage, accurate variant calls within coding sequences, the ability to explore P. falciparum within-sample diversity and to detect deletions underlying rapid diagnostic test failure.

Barcoding Genetically Distinct Plasmodium falciparum Strains for Comparative Assessment of Fitness and Antimalarial Drug Resistance.

The repeated emergence of antimalarial drug resistance in Plasmodium falciparum, including to the current frontline antimalarial artemisinin, is a perennial problem for malaria control. Next-generation sequencing has greatly accelerated the identification of polymorphisms in resistance-associated genes but has also highlighted the need for more sensitive and accurate laboratory tools to profile current and future antimalarials and to quantify the impact of drug resistance acquisition on parasite fitness. The interplay of fitness and drug response is of fundamental importance in understanding why particular genetic backgrounds are better at driving the evolution of drug resistance in natural populations, but the impact of parasite fitness landscapes on the epidemiology of drug resistance has typically been laborious to accurately quantify in the lab, with assays being limited in accuracy and throughput. Here we present a scalable method to profile fitness and drug response of genetically distinct P. falciparum strains with well-described sensitivities to several antimalarials. We leverage CRISPR/Cas9 genome-editing and barcode sequencing to track unique barcodes integrated into a nonessential gene (pfrh3). We validate this approach in multiplex competitive growth assays of three strains with distinct geographical origins. Furthermore, we demonstrate that this method can be a powerful approach for tracking artemisinin response as it can identify an artemisinin resistant strain within a mix of multiple parasite lines, suggesting an approach for scaling the laborious ring-stage survival assay across libraries of barcoded parasite lines. Overall, we present a novel high-throughput method for multiplexed competitive growth assays to evaluate parasite fitness and drug response. IMPORTANCE The complex interplay between antimalarial resistance and parasite fitness has important implications for understanding the development and spread of drug resistance alleles and the impact of genetic background on transmission. One limitation with current methodologies to measure parasite fitness is the ability to scale this beyond simple head-to-head competition experiments between a wildtype control line and test line, with a need for a scalable approach that allows tracking of parasite growth in complex mixtures. In our study, we have used CRISPR editing to insert unique DNA barcodes into a safe-harbor genomic locus to tag multiple parasite strains and use next-generation sequencing to read out strain dynamics. We observe inherent fitness differences between the strains, as well as sensitive modulation of responses to challenge with clinically relevant antimalarials, including artemisinin.

Design and implementation of multiplexed amplicon sequencing panels to serve genomic epidemiology of infectious disease: A malaria case study.

Multiplexed PCR amplicon sequencing (AmpSeq) is an increasingly popular application for cost-effective monitoring of threatened species and managed wildlife populations, and shows strong potential for the genomic epidemiology of infectious disease. AmpSeq data from infectious microbes can inform disease control in multiple ways, such as by measuring drug resistance marker prevalence, distinguishing imported from local cases, and determining the effectiveness of therapeutics. We describe the design and comparative evaluation of two new AmpSeq assays for Plasmodium falciparum malaria parasites: a four-locus panel ("4CAST") composed of highly diverse antigens, and a 129-locus panel ("AMPLseq") composed of drug resistance markers, highly diverse loci for inferring relatedness, and a locus to detect Plasmodium vivax co-infection. We explore the performance of each panel in various public health use cases with in silico simulations as well as empirical experiments. The 4CAST panel appears highly suitable for evaluating the number of distinct parasite strains within samples (complexity of infection), showing strong performance across a wide range of parasitaemia levels without a DNA pre-amplification step. For relatedness inference, the larger AMPLseq panel performs similarly to two existing panels of comparable size, despite differences in the data and approach used for designing each panel. Finally, we describe an R package (paneljudge) that facilitates the design and comparative evaluation of genetic panels for relatedness estimation, and we provide general guidance on the design and implementation of AmpSeq panels for the genomic epidemiology of infectious disease.

Defining multiplicity of vector uptake in transfected Plasmodium parasites.

The recurrent emergence of drug resistance in Plasmodium falciparum increases the urgency to genetically validate drug resistance mechanisms and identify new targets. Reverse genetics have facilitated genome-scale knockout screens in Plasmodium berghei and Toxoplasma gondii, in which pooled transfections of multiple vectors were critical to increasing scale and throughput. These approaches have not yet been implemented in human malaria species such as P. falciparum and P. knowlesi, in part because the extent to which pooled transfections can be performed in these species remains to be evaluated. Here we use next-generation sequencing to quantitate uptake of a pool of 94 barcoded vectors. The distribution of vector acquisition allowed us to estimate the number of barcodes and DNA molecules taken up by the parasite population. Dilution cloning of P. falciparum transfectants showed that individual clones possess as many as seven episomal barcodes, revealing that an intake of multiple vectors is a frequent event despite the inefficient transfection efficiency. Transfection of three spectrally-distinct fluorescent reporters allowed us to evaluate different transfection methods and revealed that schizont-stage transfection limited the tendency for parasites to take up multiple vectors. In contrast to P. falciparum, we observed that the higher transfection efficiency of P. knowlesi resulted in near complete representation of the library. These findings have important implications for how reverse genetics can be scaled in culturable Plasmodium species.

Spatio-temporal dynamics of Plasmodium falciparum transmission within a spatial unit on the Colombian Pacific Coast.

As malaria control programmes concentrate their efforts towards malaria elimination a better understanding of malaria transmission patterns at fine spatial resolution units becomes necessary. Defining spatial units that consider transmission heterogeneity, human movement and migration will help to set up achievable malaria elimination milestones and guide the creation of efficient operational administrative control units. Using a combination of genetic and epidemiological data we defined a malaria transmission unit as the area contributing 95% of malaria cases diagnosed at the catchment facility located in the town of Guapi in the South Pacific Coast of Colombia. We provide data showing that P. falciparum malaria transmission is heterogeneous in time and space and analysed, using topological data analysis, the spatial connectivity, at the micro epidemiological level, between parasite populations circulating within the unit. To illustrate the necessity to evaluate the efficacy of malaria control measures within the transmission unit in order to increase the efficiency of the malaria control effort, we provide information on the size of the asymptomatic reservoir, the nature of parasite genotypes associated with drug resistance as well as the frequency of the Pfhrp2/3 deletion associated with false negatives when using Rapid Diagnostic Tests.