Single-Cell RNA sequencing reveals cellular heterogeneity, stage transition and antigenic variation during stress adaptation in synchronized Plasmodium falciparum

The malaria parasite has a complex life cycle exhibiting phenotypic and morphogenic variations in two different hosts. Phenotypic cell-to-cell variability can be an important determinant of cellular adaptation, stress tolerance and immune evasion in the host. To investigate cellular heterogeneity, we performed single cell RNA-sequencing (scRNA-seq) of 4949 and 6873 synchronized Plasmodium cells in control and under temperature stress condition. High-resolution clustering of scRNA-seq datasets and a combination of gene signatures allow identification of cellular heterogeneity and stage transition during stress adaptation. We identified a subset of parasites primed for gametogenesis and temperature stress inducted the process of gametogenesis by upregulation of master regulator (AP2-G) of sexual conversion. Interestingly, pseudotime analysis indicated bifurcation for cell-fate decision to gametogenesis at two different stages of intra-erythrocytic cycle. Moreover, we identified a rare population of cells, which is only emerged during the stress condition, showing the reactive state of the pathogen against the temperature stress condition. Interestingly, genes associated with the gametogenesis, chaperon activity and maintenance of cellular homeostasis showed maximum variation under temperature stress condition. Thus, our study suggests that the variability and versatility of the maintenance of cellular homeostasis should enable cells to survive under different stress conditions, and may act as an important stimulator of development of drug-resistance in Plasmodium falciparum