Determining the genetic basis of fitness is central to understanding evolution and transmission of microbial pathogens. In human malaria parasites (Plasmodium falciparum), most experimental work on fitness has focused on asexual blood stage parasites, because this stage can be easily cultured, although the transmission of malaria requires both female Anopheles mosquitoes and vertebrate hosts. We explore a powerful approach to identify the genetic determinants of parasite fitness across both invertebrate and vertebrate life-cycle stages of P. falciparum. This combines experimental genetic crosses using humanized mice, with selective whole genome amplification and pooled sequencing to determine genome-wide allele frequencies and identify genomic regions under selection across multiple lifecycle stages. We applied this approach to genetic crosses between artemisinin resistant (ART-R, kelch13-C580Y) and ART-sensitive (ART-S, kelch13-WT) parasites, recently isolated from Southeast Asian patients. Two striking results emerge: we observed (i) a strong genome-wide skew (>80%) towards alleles from the ART-R parent in the mosquito stage, that dropped to ~50% in the blood stage as selfed ART-R parasites were selected against; and (ii) repeatable allele specific skews in blood stage parasites with particularly strong selection (selection coefficient (s) ≤ 0.18/asexual cycle) against alleles from the ART-R parent at loci on chromosome 12 containing MRP2 and chromosome 14 containing ARPS10. This approach robustly identifies selected loci and has strong potential for identifying parasite genes that interact with the mosquito vector or compensatory loci involved in drug resistance.

The emergence and spread of Plasmodium falciparum parasites resistant to front-line antimalarial artemisinin-combination therapies (ACT) threatens to erase the considerable gains against the disease of the last decade. Here, we develop a large-scale phenotypic screening pipeline and use it to carry out a large-scale forward-genetic phenotype screen in P. falciparum to identify genes allowing parasites to survive febrile temperatures. Screening identifies more than 200 P. falciparum mutants with differential responses to increased temperature. These mutants are more likely to be sensitive to artemisinin derivatives as well as to heightened oxidative stress. Major processes critical for P. falciparum tolerance to febrile temperatures and artemisinin include highly essential, conserved pathways associated with protein-folding, heat shock and proteasome-mediated degradation, and unexpectedly, isoprenoid biosynthesis, which originated from the ancestral genome of the parasite’s algal endosymbiont-derived plastid, the apicoplast. Apicoplast-targeted genes in general are upregulated in response to heat shock, as are other Plasmodium genes with orthologs in plant and algal genomes. Plasmodium falciparum parasites appear to exploit their innate febrile-response mechanisms to mediate resistance to artemisinin. Both responses depend on endosymbiont-derived genes in the parasite’s genome, suggesting a link to the evolutionary origins of Plasmodium parasites in free-living ancestors.

The spread of Plasmodium falciparum multidrug resistance highlights the urgency to discover new targets and chemical scaffolds. Unfortunately, lack of experimentally validated functional information about most P. falciparum genes remains a strategic hurdle. Chemogenomic profiling is an established tool for classification of drugs with similar mechanisms of action by comparing drug fitness profiles in a collection of mutants. Inferences of drug mechanisms of action and targets can be obtained by associations between shifts in drug fitness and specific genetic changes in the mutants. In this screen, P. falciparum, piggyBac single insertion mutants were profiled for altered responses to antimalarial drugs and metabolic inhibitors to create chemogenomic profiles. Drugs targeting the same pathway shared similar response profiles and multiple pairwise correlations of the chemogenomic profiles revealed novel insights into drugs’ mechanisms of action. A mutant of the artemisinin resistance candidate gene - “K13-propeller” gene (PF3D7_1343700) exhibited increased susceptibility to artemisinin drugs and identified a cluster of 7 mutants based on similar enhanced responses to the drugs tested. Our approach of chemogenomic profiling reveals artemisinin functional activity, linked by the unexpected drug-gene relationships of these mutants, to signal transduction and cell cycle regulation pathways.

BackgroundPlasmodium falciparum exhibits resistance to the artemisinin component of the frontline antimalarial treatment Artemisinin-based Combination Therapy in South East Asia. Millions of lives will be at risk if artemisinin resistance (ART-R) spreads to Africa. Single non-synonymous mutations in the propeller region of PF3D7_1343700,“K13” are implicated in resistance. In this work, we use transcriptional profiling to characterize a laboratory-generated k13 insertional mutant previously demonstrated to have increased sensitivity to artemisinins to explore the functional role of k13.ResultsA set of RNA-seq and microarray experiments confirmed that the expression profile of k13 is specifically altered during the early ring and early trophozoite stages of the mutant intraerythrocytic development cycle. The down-regulation of k13 transcripts in this mutant during the early ring stage is associated with a transcriptome advance towards a more trophozoite-like state. To discover the specific downstream effect of k13 dysregulation, we developed a new computational method to search for differential gene expression while accounting for the temporal sequence of transcription. We found that the strongest biological signature of the transcriptome shift is an up-regulation of DNA replication and repair genes during the early ring developmental stage and a down-regulation of DNA replication and repair genes during the early trophozoite stage; by contrast, the expressions of housekeeping genes are unchanged. This effect, due to k13 dysregulation, is antagonistic, such that k13 levels are negatively correlated with DNA replication and repair gene expression.ConclusionOur results support a role for k13 as a stress response regulator consistent with the hypothesis that artemisinins mode of action is oxidative stress and k13 as a functional homolog of Keap1 which in humans regulates DNA replication and repair genes in response to oxidative stress.Electronic supplementary materialThe online version of this article (10.1186/s12864-018-5207-7) contains supplementary material, which is available to authorized users.

Scientific research plays a key role in the advancement of human knowledge and pursuit of solutions to important societal challenges. Typically, research occurs within specific institutions where data are generated and subsequently analyzed. Although collaborative science bringing together multiple institutions is now common, in such collaborations the analytical processing of the data is often performed by individual researchers within the team, with only limited internal oversight and critical analysis of the workflow prior to publication. Here, we show how hackathons can be a means of enhancing collaborative science by enabling peer review before results of analyses are published by cross-validating the design of studies or underlying data sets and by driving reproducibility of scientific analyses. Traditionally, in data analysis processes, data generators and bioinformaticians are divided and do not collaborate on analyzing the data. Hackathons are a good strategy to build bridges over the traditional divide and are potentially a great agile extension to the more structured collaborations between multiple investigators and institutions.