Malaria parasites adopt a remarkable variety of morphological life stages as they transition through multiple mammalian host and mosquito vector environments. We profiled the single-cell transcriptomes of thousands of individual parasites, deriving the first high-resolution transcriptional atlas of the entire Plasmodium berghei life cycle. We then used our atlas to precisely define developmental stages of single cells from three different human malaria parasite species, including parasites isolated directly from infected individuals. The Malaria Cell Atlas provides both a comprehensive view of gene usage in a eukaryotic parasite and an open-access reference dataset for the study of malaria parasites.

BackgroundHost defense against pathogenic infection is composed of resistance and tolerance. Resistance is the ability of the host to limit a pathogen burden, whereas tolerance is the ability to limit the deleterious effects of a given pathogen burden. This distinction recognizes that the fittest host does not necessarily have the most aggressive immune system, suggesting that host-pathogen co-evolution involves more than an escalating arms race between pathogen virulence factors and host antimicrobial activity. How a host balances resistance and tolerance and how this balance influences the evolution of host defense remains unanswered. In order to determine how genotype-by-diet interactions and evolutionary costs of each strategy may constrain the evolution of host defense, we measured survival, fecundity, and pathogen burden over five days in ten genotypes of Drosophila melanogaster reared on two diets and infected with the Gram-negative bacterial pathogen Providencia rettgeri.ResultsWe demonstrated two distinct phases of infection: an acute phase that consists of high mortality, low fecundity, and high pathogen loads, and a chronic phase where there was a substantial but stable pathogen load and mortality and fecundity returned to uninfected levels. We demonstrated genetic variation for resistance in both phases of infection, but found genetic variation for tolerance only in the acute phase. We found genotype-by-diet interactions for tolerance, especially in the acute phase, but genotype-by-diet interaction did not significantly shape resistance. We found a diet-dependent positive relationship between resistance and tolerance and a weak evolutionary cost of resistance, but did not detect any costs of tolerance.ConclusionsExisting models of tolerance and resistance are overly simplistic. Multi-phase infections such as that studied here are rarely considered, but we show important differences in determination and evolutionary constraints on tolerance and resistance over the two phases of infection. Our observation of genetic variation for tolerance is inconsistent with simple models that predict evolutionary fixation of tolerance alleles, and instead indicate that genetic variation for resistance and tolerance is likely to be maintained by non-independence between resistance and tolerance, condition-dependent evolutionary costs, and environmental heterogeneity.

Summary Genes of the immune system often evolve rapidly and adaptively, presumably driven by antagonistic interactions with pathogens [1–4]. Those genes encoding secreted antimicrobial peptides (AMPs), however, have failed to exhibit conventional signatures of strong adaptive evolution, especially in arthropods (e.g., [5, 6]) and often segregate for null alleles and gene deletions [3, 4, 7, 8]. Furthermore, quantitative genetic studies have failed to associate naturally occurring polymorphism in AMP genes with variation in resistance to infection [9–11]. Both the lack of signatures of positive selection in AMPs and lack of association between genotype and immune phenotypes have yielded an interpretation that AMP genes evolve under relaxed evolutionary constraint, with enough functional redundancy that variation in, or even loss of, any particular peptide would have little effect on overall resistance [12, 13]. In stark contrast to the current paradigm, we identified a naturally occurring amino acid polymorphism in the antimicrobial peptide, Diptericin, that is highly predictive of resistance to bacterial infection in Drosophila melanogaster [13]. The identical amino acid polymorphism arose in parallel in the sister species D. simulans, by independent mutation with equivalent phenotypic effect. Convergent substitutions to arginine at the same amino acid residue have evolved at least five times across the Drosophila genus. We hypothesize that the alternative alleles are maintained by balancing selection through context-dependent or fluctuating selection. This pattern of evolution appears to be common in antimicrobial peptides, but is invisible to conventional screens for adaptive evolution that are predicated on elevated rates of amino acid divergence.

Understanding the rate of evolutionary change and the genetic architecture that facilitates rapid adaptation is a current challenge in evolutionary biology. Comparative studies show that genes with immune function are among the most rapidly evolving genes across a range of taxa. Here, we use immune defence in natural populations of to understand the rate of evolution in natural populations and the genetics underlying rapid change. We probed the immune system using the natural pathogens and to measure post-infection survival and bacterial load of wild populations collected across seasonal time along a latitudinal transect along eastern North America (Massachusetts, Pennsylvania and Virginia). There are pronounced and repeatable changes in the immune response over the approximately 10 generations between spring and autumn collections, with a significant but less distinct difference observed among geographical locations. Genes with known immune function are not enriched among alleles that cycle with seasonal time, but the immune function of a subset of seasonally cycling alleles in immune genes was tested using reconstructed outbred populations. We find that flies containing seasonal alleles in () have different functional responses to infection and that epistatic interactions among seasonal and () alleles underlie the immune phenotypes observed in natural populations. This rapid, cyclic response to seasonal environmental pressure broadens our understanding of the complex ecological and genetic interactions determining the evolution of immune defence in natural populations.

Malaria parasites have a complex life cycle featuring diverse developmental strategies, each uniquely adapted to navigate specific host environments. Here we use single-cell transcriptomics to illuminate gene usage across the transmission cycle of the most virulent agent of human malaria - Plasmodium falciparum. We reveal developmental trajectories associated with the colonization of the mosquito midgut and salivary glands and elucidate the transcriptional signatures of each transmissible stage. Additionally, we identify both conserved and non-conserved gene usage between human and rodent parasites, which point to both essential mechanisms in malaria transmission and species-specific adaptations potentially linked to host tropism. Together, the data presented here, which are made freely available via an interactive website, provide a fine-grained atlas that enables intensive investigation of the P. falciparum transcriptional journey. As well as providing insights into gene function across the transmission cycle, the atlas opens the door for identification of drug and vaccine targets to stop malaria transmission and thereby prevent disease.