Hierarchical transcriptional control regulates Plasmodium falciparum sexual differentiation

A single cell atlas of sexual development in Plasmodium falciparum

The developmental decision made by malaria parasites to become sexual underlies all malaria transmission. Here, we describe a rich atlas of short- and long-read single-cell transcriptomes of over 37,000 Plasmodium falciparum cells across intraerythrocytic asexual and sexual development. We used the atlas to explore transcriptional modules and exon usage along sexual development and expanded it to include malaria parasites collected from four Malian individuals naturally infected with multiple P. falciparum strains. We investigated genotypic and transcriptional heterogeneity within and among these wild strains at the single-cell level, finding differential expression between different strains even within the same host. These data are a key addition to the Malaria Cell Atlas interactive data resource, enabling a deeper understanding of the biology and diversity of transmission stages.

Regulation of sexual commitment in malaria parasites - a complex affair

Malaria blood stage parasites commit to either one of two distinct cellular fates while developing within erythrocytes of their mammalian host: they either undergo another round of asexual replication or they differentiate into nonreplicative transmissible gametocytes. Depending on the state of infection, either path may support or impair the ultimate goal of human-to-human transmission via the mosquito vector. Malaria parasites therefore evolved strategies to control investments into asexual proliferation versus gametocyte formation. Recent work provided fascinating molecular insight into shared and unique mechanisms underlying the control and environmental modulation of sexual commitment in the two most widely studied malaria parasite species, Plasmodium falciparum and P. berghei. With this review, we aim at placing these findings into a comparative mechanistic context.

PbARID-associated chromatin remodeling events are essential for gametocyte development in Plasmodium

Gametocyte development of the Plasmodium parasite is a key step for transmission of the parasite. Male and female gametocytes are produced from a subpopulation of asexual blood-stage parasites, but the mechanisms that regulate the differentiation of sexual stages are still under investigation. In this study, we investigated the role of PbARID, a putative subunit of a SWI/SNF chromatin remodeling complex, in transcriptional regulation during the gametocyte development of P. berghei. PbARID expression starts in early gametocytes before the manifestation of male and female-specific features, and disruption of its gene results in the complete loss of gametocytes with detectable male features and the production of abnormal female gametocytes. ChIP-seq analysis of PbARID showed that it forms a complex with gSNF2, an ATPase subunit of the SWI/SNF chromatin remodeling complex, associating with the male cis-regulatory element, TGTCT. Further ChIP-seq of PbARID in gsnf2-knockout parasites revealed an association of PbARID with another cis-regulatory element, TGCACA. RIME and DNA-binding assays suggested that HDP1 is the transcription factor that recruits PbARID to the TGCACA motif. Our results indicated that PbARID could function in two chromatin remodeling events and paly essential roles in both male and female gametocyte development.

Epigenetic regulation as a therapeutic target in the malaria parasite Plasmodium falciparum

Over the past thirty years, epigenetic regulation of gene expression has gained increasing interest as it was shown to be implicated in illnesses ranging from cancers to parasitic diseases. In the malaria parasite, epigenetics was shown to be involved in several key steps of the complex life cycle of Plasmodium, among which asexual development and sexual commitment, but also in major biological processes like immune evasion, response to environmental changes or DNA repair. Because epigenetics plays such paramount roles in the Plasmodium parasite, enzymes involved in these regulating pathways represent a reservoir of potential therapeutic targets. This review focuses on epigenetic regulatory processes and their effectors in the malaria parasite, as well as the inhibitors of epigenetic pathways and their potential as new anti-malarial drugs. Such types of drugs could be formidable tools that may contribute to malaria eradication in a context of widespread resistance to conventional anti-malarials.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. However, in low transmission settings where most mosquitoes become infected with only a single parasite clone, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. To investigate whether this clonality was potentially associated with the persistence and spatial spread of the mutation, we performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n=1,409) through estimation of identity by descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally report polymorphisms exhibiting evidence of selection for drug resistance or other phenotypes and reported a novel pfk13 mutation (G718S) as well as 61 nonsynonymous substitutions that increased markedly in frequency. However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Promising antimalarial hits from phenotypic screens: a review of recently-described multi-stage actives and their modes of action

Over the last two decades, global malaria cases caused by Plasmodium falciparum have declined due to the implementation of effective treatments and the use of insecticides. However, the COVID-19 pandemic caused major disruption in the timely delivery of medical goods and diverted public health resources, impairing malaria control. The emergence of resistance to all existing frontline antimalarials underpins an urgent need for new antimalarials with novel mechanisms of action. Furthermore, the need to reduce malaria transmission and/or prevent malaria infection has shifted the focus of antimalarial research towards the discovery of compounds that act beyond the symptomatic blood stage and also impact other parasite life cycle stages. Phenotypic screening has been responsible for the majority of new antimalarial lead compounds discovered over the past 10 years. This review describes recently reported novel antimalarial hits that target multiple parasite stages and were discovered by phenotypic screening during the COVID-19 pandemic. Their modes of action and targets in blood stage parasites are also discussed.

Machine Learning Approaches Identify Chemical Features for Stage-Specific Antimalarial Compounds

Efficacy data from diverse chemical libraries, screened against the various stages of the malaria parasite Plasmodium falciparum, including asexual blood stage (ABS) parasites and transmissible gametocytes, serve as a valuable reservoir of information on the chemical space of compounds that are either active (or not) against the parasite. We postulated that this data can be mined to define chemical features associated with the sole ABS activity and/or those that provide additional life cycle activity profiles like gametocytocidal activity. Additionally, this information could provide chemical features associated with inactive compounds, which could eliminate any future unnecessary screening of similar chemical analogs. Therefore, we aimed to use machine learning to identify the chemical space associated with stage-specific antimalarial activity. We collected data from various chemical libraries that were screened against the asexual (126 374 compounds) and sexual (gametocyte) stages of the parasite (93 941 compounds), calculated the compounds' molecular fingerprints, and trained machine learning models to recognize stage-specific active and inactive compounds. We were able to build several models that predict compound activity against ABS and dual activity against ABS and gametocytes, with Support Vector Machines (SVM) showing superior abilities with high recall (90 and 66%) and low false-positive predictions (15 and 1%). This allowed the identification of chemical features enriched in active and inactive populations, an important outcome that could be mined for essential chemical features to streamline hit-to-lead optimization strategies of antimalarial candidates. The predictive capabilities of the models held true in diverse chemical spaces, indicating that the ML models are therefore robust and can serve as a prioritization tool to drive and guide phenotypic screening and medicinal chemistry programs.

The global transcriptome of Plasmodium falciparum mid-stage gametocytes (stages II-IV) appears largely conserved and gametocyte-specific gene expression patterns vary in clinical isolates

Our overall understanding of the developmental biology of malaria parasites has been greatly enhanced by recent advances in transcriptomic analysis. However, most of these investigations rely on laboratory strains (LS) that were adapted into in vitro culture many years ago, and the transcriptomes of clinical isolates (CI) circulating in human populations have not been assessed. In this study, RNA-seq was used to compare the global transcriptome of mid-stage gametocytes derived from three short-term cultured CI, with gametocytes derived from the NF54 reference laboratory strain. The core transcriptome appeared to be consistent between CI- and LS-derived gametocyte preparations, but some important differences were also observed. A majority of gametocyte-specific genes (43/53) appear to have relatively higher expression in CI-derived gametocytes than in LS-derived gametocytes, but a K-means clustering analysis showed that genes involved in flagellum- and microtubule-based processes (movement/motility) were more abundant in both groups, albeit with some differences between them. In addition, gametocytes from one CI described as CI group II gametocytes (CI:GGII) showed gene expression variation in the form of reduced gametocyte-specific gene expression compared to the other two CI-derived gametocytes (CI gametocyte group I, CI:GGI), although the mixed developmental stages used in our study is a potential confounder, only partially mitigated by the inclusion of multiple replicates for each CI. Overall, our study suggests that there may be subtle differences in the gene expression profiles of mid-stage gametocytes from CI relative to the NF54 reference strain of Plasmodium falciparum. Thus, it is necessary to deploy gametocyte-producing clinical parasite isolates to fully understand the diversity of gene expression strategies that may occur during the sequestered development of parasite sexual stages. IMPORTANCE Maturing gametocytes of Plasmodium falciparum are known to sequester away from peripheral circulation into the bone marrow until they are mature. Blocking gametocyte sequestration can prevent malaria transmission from humans to mosquitoes, but most studies aim to understand gametocyte development utilizing long-term adapted laboratory lines instead of clinical isolates. This is a particular issue for our understanding of the sexual stages, which are known to decrease rapidly during adaptation to long-term culture, meaning that many LS are unable to produce transmissible gametocytes. Using RNA-seq, we investigated the global transcriptome of mid-stage gametocytes derived from three clinical isolates and a reference strain (NF54). This identified important differences in gene expression profiles between immature gametocytes of CI and the NF54 reference strain of P. falciparum, suggesting increased investment in gametocytogenesis in clinical isolates. Our transcriptomic data highlight the use of clinical isolates in studying the morphological, cellular features and molecular biology of gametocytes.

The role of long noncoding RNAs in malaria parasites

The human malaria parasites, including Plasmodium falciparum, persist as a major cause of global morbidity and mortality. The recent stalling of progress toward malaria elimination substantiates a need for novel interventions. Controlled gene expression is central to the parasite's numerous life cycle transformations and adaptation. With few specific transcription factors (TFs) identified, crucial roles for chromatin states and epigenetics in parasite transcription have become evident. Although many chromatin-modifying enzymes are known, less is known about which factors mediate their impacts on transcriptional variation. Like those of higher eukaryotes, long noncoding RNAs (lncRNAs) have recently been shown to have integral roles in parasite gene regulation. This review aims to summarize recent developments and key findings on the role of lncRNAs in P. falciparum.

Phenotypic Screens Identify Genetic Factors Associated with Gametocyte Development in the Human Malaria Parasite Plasmodium falciparum

Transmission of the deadly malaria parasite Plasmodium falciparum from humans to mosquitoes is achieved by specialized intraerythrocytic sexual forms called gametocytes. Though the crucial regulatory mechanisms leading to gametocyte commitment have recently come to light, networks of genes that control sexual development remain to be elucidated. Here, we report a pooled-mutant screen to identify genes associated with gametocyte development in P. falciparum. Our results categorized genes that modulate gametocyte progression as hypoproducers or hyperproducers of gametocytes, and the in-depth analysis of individual clones confirmed phenotypes in sexual commitment rates and putative functions in gametocyte development. We present a new set of genes that have not been implicated in gametocytogenesis before and demonstrate the potential of forward genetic screens in isolating genes impacting parasite sexual biology, an exciting step toward the discovery of new antimalarials for a globally significant pathogen. IMPORTANCE Blocking human-to-vector transmission is an essential step toward malaria elimination. Gametocytes are solely responsible for achieving this transmission and represent an opportunity for therapeutic intervention. While these falciform-shaped parasite stages were first discovered in the 1880s, our understanding of the genetic determinants responsible for their formation and molecular mechanisms that drive their development is limited. In this work, we developed a scalable screening methodology with piggyBac mutants to identify genes that influence the development of gametocytes in the most lethal human malaria parasite, P. falciparum. By doing so, we lay the foundation for large-scale functional genomic studies specifically designed to address remaining questions about sexual commitment, maturation, and mosquito infection in P. falciparum. Such functional genetic screens will serve to expedite the identification of essential pathways and processes for the development of novel transmission-blocking agents.

Heat-shock responses: systemic and essential ways of malaria parasite survival

Fever is a part of the human innate immune response that contributes to limiting microbial growth and development in many infectious diseases. For the parasite Plasmodium falciparum, survival of febrile temperatures is crucial for its successful propagation in human populations as well as a fundamental aspect of malaria pathogenesis. This review discusses recent insights into the biological complexity of the malaria parasite's heat-shock response, which involves many cellular compartments and essential metabolic processes to alleviate oxidative stress and accumulation of damaged and unfolded proteins. We highlight the overlap between heat-shock and artemisinin resistance responses, while also explaining how the malaria parasite adapts its fever response to fight artemisinin treatment. Additionally, we discuss how this systemic and essential fight for survival can also contribute to parasite transmission to mosquitoes.

Isoliensinine from Cissampelos pariera rhizomes exhibits potential gametocytocidal and anti-malarial activities against Plasmodium falciparum clinical isolates

BACKGROUND

The unmet demand for effective malaria transmission-blocking agents targeting the transmissible stages of Plasmodium necessitates intensive discovery efforts. In this study, a bioactive bisbenzylisoquinoline (BBIQ), isoliensinine, from Cissampelos pariera (Menispermaceae) rhizomes was identified and characterized for its anti-malarial activity.

METHODS

Malaria SYBR Green I fluorescence assay was performed to evaluate the in vitro antimalarial activity against D6, Dd2, and F32-ART5 clones, and immediate ex vivo (IEV) susceptibility for 10 freshly collected P. falciparum isolates. To determine the speed- and stage-of-action of isoliensinine, an IC₅₀ speed assay and morphological analyses were performed using synchronized Dd2 asexuals. Gametocytocidal activity against two culture-adapted gametocyte-producing clinical isolates was determined using microscopy readouts, with possible molecular targets and their binding affinities deduced in silico.

RESULTS

Isoliensinine displayed a potent in vitro gametocytocidal activity at mean IC₅₀^gam values ranging between 0.41 and 0.69 µM for Plasmodium falciparum clinical isolates. The BBIQ compound also inhibited asexual replication at mean IC₅₀^Asexual of 2.17 µM, 2.22 µM, and 2.39 µM for D6, Dd2 and F32-ART5 respectively, targeting the late-trophozoite to schizont transition. Further characterization demonstrated a considerable immediate ex vivo potency against human clinical isolates at a geometric mean IC₅₀^IEV = 1.433 µM (95% CI 0.917-2.242). In silico analyses postulated a probable anti-malarial mechanism of action by high binding affinities for four mitotic division protein kinases; Pfnek1, Pfmap2, Pfclk1, and Pfclk4. Additionally, isoliensinine was predicted to possess an optimal pharmacokinetics profile and drug-likeness properties.

CONCLUSION

These findings highlight considerable grounds for further exploration of isoliensinine as an amenable scaffold for malaria transmission-blocking chemistry and target validation.

Identification of Potential Drug Targets in Erythrocyte Invasion Pathway of Plasmodium falciparum

The erythrocyte invasion phase plays a critical role in multiplication, sexual determination, and drug resistance in Plasmodium falciparum. In order to identify the critical genes and pathways in the erythrocyte invasion phase, the gene set (GSE129949) and the RNA-Seq count data for the W2mef strain were used for further analysis. An integrative bioinformatics study was performed to scrutinize genes as potential drug targets. 487 differentially expressed genes (DEGs) with an adjusted P value < 0.001 enriched 47 Gene Ontology (GO) terms that were over-represented based on hyper-geometric analysis P value < 0.01. Protein-Protein interaction network analysis was done using DEGs with higher confidence interactions (PPI score threshold = 0.7). MCODE and cytoHubba apps were utilized to define the hub proteins and rank them based on multiple topological analyses and MCODE scores. Furthermore, Gene Set Enrichment Analysis (GSEA) was carried out by using 322 gene sets from the MPMP database. The genes involved in multiple significant gene sets were determined by leading-edge analysis. Our study identified six genes encoding proteins that could be potential drug targets involved in the erythrocyte invasion phase related to merozoites motility, cell-cycle regulation, G-dependent protein kinase phosphorylation in schizonts, control of microtubule assembly, and sexual commitment. The druggability of those proteins was calculated based on the DCI (Drug Confidence Index) and predicted binding pockets' values. The protein that showed the best binding pocket value was subjected to deep learning-based virtual screening. The study identified the best small molecule inhibitors in terms of drug-binding score against the proteins for inhibitor identification.

Plasmodium falciparum gametocytes display global chromatin remodelling during sexual differentiation

BACKGROUND

The protozoan malaria parasite Plasmodium falciparum has a complex life cycle during which it needs to differentiate into multiple morphologically distinct life forms. A key process for transmission of the disease is the development of male and female gametocytes in the human blood, yet the mechanisms determining sexual dimorphism in these haploid, genetically identical sexual precursor cells remain largely unknown. To understand the epigenetic program underlying the differentiation of male and female gametocytes, we separated the two sexual forms by flow cytometry and performed RNAseq as well as comprehensive ChIPseq profiling of several histone variants and modifications.

RESULTS

We show that in female gametocytes the chromatin landscape is globally remodelled with respect to genome-wide patterns and combinatorial usage of histone variants and histone modifications. We identified sex specific differences in heterochromatin distribution, implicating exported proteins and ncRNAs in sex determination. Specifically in female gametocytes, the histone variants H2A.Z/H2B.Z were highly enriched in H3K9me3-associated heterochromatin. H3K27ac occupancy correlated with stage-specific gene expression, but in contrast to asexual parasites this was unlinked to H3K4me3 co-occupancy at promoters in female gametocytes.

CONCLUSIONS

Collectively, we defined novel combinatorial chromatin states differentially organising the genome in gametocytes and asexual parasites and unravelled fundamental, sex-specific differences in the epigenetic code. Our chromatin maps represent an important resource for future understanding of the mechanisms driving sexual differentiation in P. falciparum.

Regulators of male and female sexual development are critical for the transmission of a malaria parasite

Malaria transmission to mosquitoes requires a developmental switch in asexually dividing blood-stage parasites to sexual reproduction. In Plasmodium berghei, the transcription factor AP2-G is required and sufficient for this switch, but how a particular sex is determined in a haploid parasite remains unknown. Using a global screen of barcoded mutants, we here identify genes essential for the formation of either male or female sexual forms and validate their importance for transmission. High-resolution single-cell transcriptomics of ten mutant parasites portrays the developmental bifurcation and reveals a regulatory cascade of putative gene functions in the determination and subsequent differentiation of each sex. A male-determining gene with a LOTUS/OST-HTH domain as well as the protein interactors of a female-determining zinc-finger protein indicate that germ-granule-like ribonucleoprotein complexes complement transcriptional processes in the regulation of both male and female development of a malaria parasite.

Host-Erythrocytic Sphingosine-1-Phosphate Regulates Plasmodium Histone Deacetylase Activity and Exhibits Epigenetic Control over Cell Death and Differentiation

The evolution of resistance to practically all antimalarial drugs poses a challenge to the current malaria elimination and eradication efforts. Given that the epigenome of Plasmodium falciparum governs several crucial parasite functions, pharmaceutical interventions with transmission-blocking potential that target epigenetic molecular markers and regulatory mechanisms are likely to encounter drug resistance. In the malaria parasite, histone deacetylases (HDACs) are essential epigenetic modulators that regulate cellular transcriptional rearrangements, notably the molecular mechanisms underlying parasite proliferation and differentiation. We establish "lipid sequestration" as a mechanism by which sphingolipids, specifically Sphingosine-1-Phosphate (S1P) (a metabolic product of Sphingosine Kinase 1 [SphK-1]), regulate epigenetic reprogramming in the parasite by interacting with, and modulating, the histone-deacetylation activity of PfHDAC-1, thereby regulating Plasmodium pathogenesis. Furthermore, we demonstrate that altering host S1P levels with PF-543, a potent and selective Sphk-1 inhibitor, dysregulates PfHDAC-1 activity, resulting in a significant increase in the global histone acetylation signals and, consequently, transcriptional modulation of genes associated with gametocytogenesis, virulence, and proliferation. Our findings point to a hitherto unrecognized functional role for host S1P-mediated sphingolipid signaling in modulating PfHDAC-1's enzymatic activity and, as a result, the parasite's dynamic genome-wide transcriptional patterns. The epigenetic regulation of parasite proliferation and sexual differentiation offers a novel approach for developing host-targeted therapeutics to combat malaria resistance to conventional regimens. IMPORTANCE Sphingolipid is an 18-carbon amino-alcohol-containing lipid with a sphingosine backbone, which when phosphorylated by sphingosine kinase 1 (SphK-1), generates sphingosine-1-phosphate (S1P), an essential lipid signaling molecule. Dysregulation of S1P function has been observed in a variety of pathologies, including severe malaria. The malaria parasite Plasmodium acquires a host S1P pool for its growth and survival. Here, we describe the molecular attuning of histone deacetylase-1 (PfHDAC-1), a crucial epigenetic modulator that contributes to the establishment of epigenetic chromatin states and parasite survival, in response to S1P binding. Our findings highlight the host lipid-mediated epigenetic regulation of malaria parasite key genes.

High-throughput analysis of the transcriptional patterns of sexual genes in malaria

BACKGROUND

Plasmodium falciparum (Pf) is the leading protozoan causing malaria, the most devastating parasitic disease. To ensure transmission, a small subset of Pf parasites differentiate into the sexual forms (gametocytes). Since the abundance of these essential parasitic forms is extremely low within the human host, little is currently known about the molecular regulation of their sexual differentiation, highlighting the need to develop tools to investigate Pf gene expression during this fundamental mechanism.

METHODS

We developed a high-throughput quantitative Reverse-Transcription PCR (RT-qPCR) platform to robustly monitor Pf transcriptional patterns, in particular, systematically profiling the transcriptional pattern of a large panel of gametocyte-related genes (GRG). Initially, we evaluated the technical performance of the systematic RT-qPCR platform to ensure it complies with the accepted quality standards for: (i) RNA extraction, (ii) cDNA synthesis and (iii) evaluation of gene expression through RT-qPCR. We then used this approach to monitor alterations in gene expression of a panel of GRG upon treatment with gametocytogenesis regulators.

RESULTS

We thoroughly elucidated GRG expression profiles under treatment with the antimalarial drug dihydroartemisinin (DHA) or the metabolite choline over the course of a Pf blood cycle (48 h). We demonstrate that both significantly alter the expression pattern of PfAP2-G, the gametocytogenesis master regulator. However, they also markedly modify the developmental rate of the parasites and thus might bias the mRNA expression. Additionally, we screened the effect of the metabolites lactate and kynurenic acid, abundant in severe malaria, as potential regulators of gametocytogenesis.

CONCLUSIONS

Our data demonstrate that the high-throughput RT-qPCR method enables studying the immediate transcriptional response initiating gametocytogenesis of the parasites from a very low volume of malaria-infected RBC samples. The obtained data expand the current knowledge of the initial alterations in mRNA profiles of GRG upon treatment with reported regulators. In addition, using this method emphasizes that asexual parasite stage composition is a crucial element that must be considered when interpreting changes in GRG expression by RT-qPCR, specifically when screening for novel compounds that could regulate Pf sexual differentiation.

Molecular mechanisms of cellular quiescence in apicomplexan parasites

Quiescence is a reversible nonproliferative cellular state that allows organisms to persist through unfavorable conditions. Quiescence can be stimulated by a wide range of external or intrinsic factors. Cells undergo a coordinated molecular program to enter and exit from the quiescent state, which is governed by signaling, transcriptional and translational changes, epigenetic mechanisms, metabolic switches, and changes in cellular architecture. These mechanisms have been extensively studied in model organisms, and a growing number of studies have identified conserved mechanisms in apicomplexan parasites. Quiescence in the context of a parasitic infection has significant clinical impact: quiescent forms may underlie treatment failures, relapsing infections, and stress tolerance. Here, we review the latest understanding of quiescence in apicomplexa, synthesizing these studies to highlight conserved mechanisms, and identifying technologies to assist in further characterization of quiescence. Understanding conserved mechanisms of quiescence in apicomplexans will provide avenues for transmission prevention and radical cure of infections.

PfHMGB2 has a role in malaria parasite mosquito infection

Differentiation of asexually replicating parasites into gametocytes is critical for successful completion of the sexual phase of the malaria parasite life cycle. Gametes generated from gametocytes fuse to form a zygote which differentiates into ookinetes and oocysts. The sporozoites are formed inside oocysts which migrate to the salivary glands for next cycle of human infection. These morphologically and functionally distinct stages require stage-specific gene expression via specific transcriptional regulators. The capacity of high mobility group box (HMGB) proteins to interact with DNA in a sequence independent manner enables them to regulate higher order chromosome organization and regulation of gene expression. Plasmodium falciparum HMGB2 (PfHMGB2) shows a typical L- shaped predicted structure which is similar to mammalian HMG box proteins and shows very high protein sequence similarity to PyHMGB2 and PbHMGB2. Functional characterization of PfHMGB2 by gene deletion (Pfhmgb2¯) showed that knockout parasites develop normally as asexual stages and undergo gametocytogenesis. Transmission experiments revealed that Pfhmgb2¯ can infect mosquitoes and develop as oocyst stages. However, transmission was reduced compared to wild type (WT) parasites and as a consequence, the salivary gland sporozoites were reduced in number. In summary, we demonstrate that PfHMGB2 has no role in asexual growth and a modest role in sexual phase development and parasite transmission to the mosquito.

Comparative single-cell transcriptional atlases of Babesia species reveal conserved and species-specific expression profiles

Babesia is a genus of apicomplexan parasites that infect red blood cells in vertebrate hosts. Pathology occurs during rapid replication cycles in the asexual blood stage of infection. Current knowledge of Babesia replication cycle progression and regulation is limited and relies mostly on comparative studies with related parasites. Due to limitations in synchronizing Babesia parasites, fine-scale time-course transcriptomic resources are not readily available. Single-cell transcriptomics provides a powerful unbiased alternative for profiling asynchronous cell populations. Here, we applied single-cell RNA sequencing to 3 Babesia species (B. divergens, B. bovis, and B. bigemina). We used analytical approaches and algorithms to map the replication cycle and construct pseudo-synchronized time-course gene expression profiles. We identify clusters of co-expressed genes showing "just-in-time" expression profiles, with gradually cascading peaks throughout asexual development. Moreover, clustering analysis of reconstructed gene curves reveals coordinated timing of peak expression in epigenetic markers and transcription factors. Using a regularized Gaussian graphical model, we reconstructed co-expression networks and identified conserved and species-specific nodes. Motif analysis of a co-expression interactome of AP2 transcription factors identified specific motifs previously reported to play a role in DNA replication in Plasmodium species. Finally, we present an interactive web application to visualize and interactively explore the datasets.

The Plasmodium falciparum Nuclear Protein Phosphatase NIF4 Is Required for Efficient Merozoite Invasion and Regulates Artemisinin Sensitivity

Artemisinin resistance in Plasmodium falciparum has been associated with a mutation in the NLI-interacting factor-like phosphatase PfNIF4, in addition to the mutations in the Kelch13 protein as the major determinant. We found that PfNIF4 was predominantly expressed at the schizont stage and localized in the nuclei of the parasite. To elucidate the functions of PfNIF4 in P. falciparum, we performed PfNIF4 knockdown (KD) using the inducible ribozyme system. PfNIF4 KD attenuated merozoite invasion and affected gametocytogenesis. PfNIF4 KD parasites also showed significantly increased in vitro susceptibility to artemisinins. Transcriptomic and proteomic analysis revealed that PfNIF4 KD led to the downregulation of gene categories involved in invasion and artemisinin resistance (e.g., mitochondrial function, membrane, and Kelch13 interactome) at the trophozoite and/or schizont stage. Consistent with PfNIF4 being a protein phosphatase, PfNIF4 KD resulted in an overall upregulation of the phosphoproteome of infected erythrocytes. Quantitative phosphoproteomic profiling identified a set of PfNIF4-regulated phosphoproteins with functional similarity to FCP1 substrates, particularly proteins involved in chromatin organization and transcriptional regulation. Specifically, we observed increased phosphorylation of Ser2/5 of the tandem repeats in the C-terminal domain (CTD) of RNA polymerase II (RNAPII) upon PfNIF4 KD. Furthermore, using the TurboID-based proteomic approach, we identified that PfNIF4 interacted with the RNAPII components, AP2-domain transcription factors, and chromatin-modifiers and binders. These findings suggest that PfNIF4 may act as the RNAPII CTD phosphatase, regulating the expression of general and parasite-specific cellular pathways during the blood-stage development.

Distinct Histone Post-translational Modifications during Plasmodium falciparum Gametocyte Development

Histones are the building units of nucleosomes, which constitute chromatin. Histone post-translational modifications (PTMs) play an essential role in epigenetic gene regulation. The Plasmodium falciparum genome encodes canonical and variant histones and a collection of conserved enzymes for histone PTMs and chromatin remodeling. Herein, we profiled the P. falciparum histone PTMs during the development of gametocytes, the obligatory stage for parasite transmission. Mass spectrometric analysis of histones extracted from the early, middle, and late stages of gametocytes identified 457 unique histone peptides with 90 PTMs, of which 50% were novel. The gametocyte histone PTMs display distinct patterns from asexual stages, with many new methylation sites in histones H3 and H3.3 (e.g., K14, K18, and K37). Quantitative analyses revealed a high abundance of acetylation in H3 and H4, mono-methylation of H3/H3.3 K37, and ubiquitination of H3BK112, suggesting that these PTMs play critical roles in gametocytes. Gametocyte histones also showed extensive and unique combinations of PTMs. These data indicate that the parasite harbors distinct transcription regulation mechanisms during gametocyte development and lay the foundation for further characterization of epigenetic regulation in the life cycle of the malaria parasite.

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

Like most human pathogens, the malaria parasite Plasmodium falciparum experiences strong selection pressure from public health interventions such as drug treatment. While most commonly studied in the context of drug targets and related pathways, parasite adaptation to control measures likely extends to phenotypes beyond drug resistance. Here, we use modeling to explore how control measures can reduce levels of within-host competition between P. falciparum genotypes and favor higher rates of sexual investment. We validate these predictions with longitudinally sampled genomic data from French Guiana during a period of malaria decline and find that the most strongly selected genes are enriched for transcription factors involved in commitment to and development of the parasite’s sexual gametocyte form.

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.

ApiAP2 Gene-Network Regulates Gametocytogenesis in Plasmodium Parasites

Malaria is a mosquito-borne infectious disease, caused by unicellular Apicomplexan protozoa of the genus Plasmodium. The sexual stage of Plasmodium is one of the most fascinating aspects of the Plasmodium life cycle, yet relatively less explored until now. The production of sexually fit gametocytes through gametocytogenesis is essential to the transmission of the Plasmodium parasite into an anopheline mosquito vector. Understanding how gametocytogenesis is regulated promotes the identification of novel drug targets and also the development of transmission-blocking vaccines that would help reduce the disease burden in endemic areas. Transcriptional regulation in Plasmodium parasites is primarily controlled by a family of twenty-seven Apicomplexan Apetela 2 (ApiAP2) genes which act in a cascade to enable the parasite to progress through its asexual replication as well as gametocytogenesis. Here, we review the latest progress made on members of the ApiAP2 family characterized as key players of the transcriptional machinery of gametocytes. Further, we will highlight the transcriptional regulation network of ApiAP2 genes at each stage of gametocytogenesis.

Streamlined and Robust Stage-Specific Profiling of Gametocytocidal Compounds Against Plasmodium falciparum

Malaria elimination is dependent on the ability to target both the pathogenic and transmissible stages of the human malaria parasite, Plasmodium falciparum. These forms of the parasite are differentiated by unique developmental stages, each with their own biological mechanisms and processes. These individual stages therefore also respond differently to inhibitory compounds, and this complicates the discovery of multistage active antimalarial agents. The search for compounds with transmission-blocking activity has focused on screening for activity on mature gametocytes, with only limited descriptions available for the activity of such compounds on immature stage gametocytes. This therefore poses a gap in the profiling of antimalarial agents for pan-reactive, multistage activity to antimalarial leads. Here, we optimized an effective and robust strategy for the simple and cost-effective description of the stage-specific action of gametocytocidal antimalarial compounds.

Adapt or Die: Targeting Unique Transmission-Stage Biology for Malaria Elimination

Plasmodium parasites have a complex life cycle that includes development in the human host as well as the Anopheles vector. Successful transmission of the parasite between its host and vector therefore requires the parasite to balance its investments in asexual replication and sexual reproduction, varying the frequency of sexual commitment to persist within the human host and generate future opportunities for transmission. The transmission window is extended further by the ability of stage V gametocytes to circulate in peripheral blood for weeks, whereas immature stage I to IV gametocytes sequester in the bone marrow and spleen until final maturation. Due to the low gametocyte numbers in blood circulation and with the ease of targeting such life cycle bottlenecks, transmission represents an efficient target for therapeutic intervention. The biological process of Plasmodium transmission is a multistage, multifaceted process and the past decade has seen a much deeper understanding of the molecular mechanisms and regulators involved. Clearly, specific and divergent processes are used during transmission compared to asexual proliferation, which both poses challenges but also opportunities for discovery of transmission-blocking antimalarials. This review therefore presents an update of our molecular understanding of gametocyte and gamete biology as well as the status of transmission-blocking activities of current antimalarials and lead development compounds. By defining the biological components associated with transmission, considerations for the development of new transmission-blocking drugs to target such untapped but unique biology is suggested as an important, main driver for transmission-blocking drug discovery.

PfARID Regulates P. falciparum Malaria Parasite Male Gametogenesis and Female Fertility and Is Critical for Parasite Transmission to the Mosquito Vector

Sexual reproduction of Plasmodium falciparum parasites is critical to the spread of malaria in the human population. The factors that regulate gene expression underlying formation of fertilization-competent gametes, however, remain unknown. Here, we report that P. falciparum expresses a protein with an AT-rich interaction domain (ARID) which, in other organisms, is part of chromatin remodeling complexes. P. falciparum ARID (PfARID) localized to the parasite nucleus and is critical for the formation of male gametes and fertility of female gametes. PfARID gene deletion (Pfarid^–) gametocytes showed downregulation of gene expression important for gametogenesis, antigenic variation, and cell signaling and for parasite development in the mosquito. Our study identifies PfARID as a critical nuclear protein involved in regulating the gene expression landscape of mature gametocytes. This establishes fertility and also prepares the parasite for postfertilization events that are essential for infection of the mosquito vector.

Identification of Co-Existing Mutations and Gene Expression Trends Associated With K13-Mediated Artemisinin Resistance in Plasmodium falciparum

Plasmodium falciparum infects millions and kills thousands of people annually the world over. With the emergence of artemisinin and/or multidrug resistant strains of the pathogen, it has become even more challenging to control and eliminate the disease. Multiomics studies of the parasite have started to provide a glimpse into the confounding genetics and mechanisms of artemisinin resistance and identified mutations in Kelch13 (K13) as a molecular marker of resistance. Over the years, thousands of genomes and transcriptomes of artemisinin-resistant/sensitive isolates have been documented, supplementing the search for new genes/pathways to target artemisinin-resistant isolates. This meta-analysis seeks to recap the genetic landscape and the transcriptional deregulation that demarcate artemisinin resistance in the field. To explore the genetic territory of artemisinin resistance, we use genomic single-nucleotide polymorphism (SNP) datasets from 2,517 isolates from 15 countries from the MalariaGEN Network (The Pf3K project, pilot data release 4, 2015) to dissect the prevalence, geographical distribution, and co-existing patterns of genetic markers associated with/enabling artemisinin resistance. We have identified several mutations which co-exist with the established markers of artemisinin resistance. Interestingly, K13-resistant parasites harbor α-ß hydrolase and putative HECT domain-containing protein genes with the maximum number of SNPs. We have also explored the multiple, publicly available transcriptomic datasets to identify genes from key biological pathways whose consistent deregulation may be contributing to the biology of resistant parasites. Surprisingly, glycolytic and pentose phosphate pathways were consistently downregulated in artemisinin-resistant parasites. Thus, this meta-analysis highlights the genetic and transcriptomic features of resistant parasites to propel further exploratory studies in the community to tackle artemisinin resistance.

A Putative Plasmodium RNA-Binding Protein Plays a Critical Role in Female Gamete Fertility and Parasite Transmission to the Mosquito Vector

Plasmodium falciparum sexual stage gametocytes are critical for parasite transmission from the human host to the mosquito vector. Mature gametocytes generate fertile male (micro-) or female (macro-) gametes upon activation inside the mosquito midgut. While a number of parasite genes have been described that are critical for P. falciparum gametogenesis and fertility, no parasite gene has been shown to have a unique function in macrogametes. The genome of P. falciparum encodes numerous RNA-binding proteins. We identified a novel protein containing a putative RNA-binding domain, which we named Macrogamete-Contributed Factor Essential for Transmission (MaCFET). This protein is expressed in the asexual and sexual stages. Parasites that carry a deletion of MaCFET (Pfmacfet¯), developed normally as asexual stages, indicating that its function is not essential for the asexual proliferation of the parasite in vitro. Furthermore, Pfmacfet¯ male and female gametocytes developed normally and underwent activation to form microgametes and macrogametes. However, by utilizing genetic crosses, we demonstrate that Pfmacfet¯ parasites suffer a complete female-specific defect in successful fertilization. Therefore, PfMaCFET is a critical female-contributed factor for parasite transmission to the mosquito. Based on its putative RNA-binding properties, PfMaCFET might be in involved in the regulation of mRNAs that encode female-specific functions for fertilization or female-contributed factors needed post fertilization.

The transcriptome from asexual to sexual in vitro development of Cystoisospora suis (Apicomplexa: Coccidia)

The apicomplexan parasite Cystoisospora suis is an enteropathogen of suckling piglets with woldwide distribution. As with all coccidian parasites, its lifecycle is characterized by asexual multiplication followed by sexual development with two morphologically distinct cell types that presumably fuse to form a zygote from which the oocyst arises. However, knowledge of the sexual development of C. suis is still limited. To complement previous in vitro studies, we analysed transcriptional profiles at three different time points of development (corresponding to asexual, immature and mature sexual stages) in vitro via RNASeq. Overall, transcription of genes encoding proteins with important roles in gametes biology, oocyst wall biosynthesis, DNA replication and axonema formation as well as proteins with important roles in merozoite biology was identified. A homologue of an oocyst wall tyrosine rich protein of Toxoplasma gondii was expressed in macrogametes and oocysts of C. suis. We evaluated inhibition of sexual development in a host-free culture for C. suis by antiserum specific to this protein to evaluate whether it could be exploited as a candidate for control strategies against C. suis. Based on these data, targets can be defined for future strategies to interrupt parasite transmission during sexual development.

Histone Modification Landscapes as a Roadmap for Malaria Parasite Development

Plasmodium falciparum remains the deadliest parasite species in the world, responsible for 229 million cases of human malaria in 2019. The ability of the P. falciparum parasite to progress through multiple life cycle stages and thrive in diverse host and vector species hinges on sophisticated mechanisms of epigenetic regulation of gene expression. Emerging evidence indicates such epigenetic control exists in concentric layers, revolving around core histone post-translational modification (PTM) landscapes. Here, we provide a necessary update of recent epigenome research in malaria parasites, focusing specifically on the ability of dynamic histone PTM landscapes to orchestrate the divergent development and differentiation pathways in P. falciparum parasites. In addition to individual histone PTMs, we discuss recent findings that imply functional importance for combinatorial PTMs in P. falciparum parasites, representing an operational histone code. Finally, this review highlights the remaining gaps and provides strategies to address these to obtain a more thorough understanding of the histone modification landscapes that are at the center of epigenetic regulation in human malaria parasites.

The transcriptional regulator HDP1 controls expansion of the inner membrane complex during early sexual differentiation of malaria parasites

Transmission of Plasmodium falciparum and other malaria parasites requires their differentiation from asexual blood stages into gametocytes, the non-replicative sexual stage necessary to infect the mosquito vector. This transition involves changes in gene expression and chromatin reorganization that result in the activation and silencing of stage-specific genes. However, the genomes of malaria parasites have been noted for their limited number of transcriptional and chromatin regulators, and the molecular mediators of these changes remain largely unknown. We recently identified homeodomain protein 1 (HDP1) as a DNA-binding protein, first expressed in gametocytes, that enhances the expression of key genes critical for early sexual differentiation. The discovery of HDP1 marks a new class of transcriptional regulator in malaria parasites outside of the better-characterized ApiAP2 family. Here, using molecular biology, biochemistry and microscopy techniques, we show that HDP1 is essential for gametocyte maturation, facilitating the necessary upregulation of inner membrane complex components during early gametocytogenesis that gives P. falciparum gametocytes their characteristic shape.

A dynamic and combinatorial histone code drives malaria parasite asexual and sexual development

Histone posttranslational modifications (PTMs) frequently co-occur on the same chromatin domains or even in the same molecule. It is now established that these “histone codes” are the result of cross talk between enzymes that catalyze multiple PTMs with univocal readout as compared with these PTMs in isolation. Here, we performed a comprehensive identification and quantification of histone codes of the malaria parasite, Plasmodium falciparum. We used advanced quantitative middle-down proteomics to identify combinations of PTMs in both the proliferative, asexual stages and transmissible, sexual gametocyte stages of P. falciparum. We provide an updated, high-resolution compendium of 77 PTMs on H3 and H3.3, of which 34 are newly identified in P. falciparum. Coexisting PTMs with unique stage distinctions were identified, indicating that many of these combinatorial PTMs are associated with specific stages of the parasite life cycle. We focused on the code H3R17me2K18acK23ac for its unique presence in mature gametocytes; chromatin proteomics identified a gametocyte-specific SAGA-like effector complex including the transcription factor AP2-G2, which we tied to this specific histone code, as involved in regulating gene expression in mature gametocytes. Ultimately, this study unveils previously undiscovered histone PTMs and their functional relationship with coexisting partners. These results highlight that investigating chromatin regulation in the parasite using single histone PTM assays might overlook higher-order gene regulation for distinct proliferation and differentiation processes.

•

First middle-down chromatin proteomics compendium of the malaria parasite, Plasmodium falciparum.

•

Novel histone PTMs (including arginine methylation) in both asexual parasites and transmissible gametocytes.

•

Histone PTM cross talk is dynamic life cycle stage stratified.

•

Gametocytes rely on histone PTM connectivity to allow onward transmission.

•

AP2-G2 is an important effector of H3K18acK23ac in mature gametocytes.

A transcriptional switch controls sex determination in Plasmodium falciparum

Sexual reproduction and meiotic sex are deeply rooted in the eukaryotic tree of life, but mechanisms determining sex or mating types are extremely varied and are only well characterized in a few model organisms¹. In malaria parasites, sexual reproduction coincides with transmission to the vector host. Sex determination is non-genetic, with each haploid parasite capable of producing either a male or a female gametocyte in the human host². The hierarchy of events and molecular mechanisms that trigger sex determination and maintenance of sexual identity are yet to be elucidated. Here we show that the male development 1 (md1) gene is both necessary and sufficient for male fate determination in the human malaria parasite Plasmodium falciparum. We show that Md1 has a dual function stemming from two separate domains: in sex determination through its N terminus and in male development from its conserved C-terminal LOTUS/OST-HTH domain. We further identify a bistable switch at the md1 locus, which is coupled with sex determination and ensures that the male-determining gene is not expressed in the female lineage. We describe one of only a few known non-genetic mechanisms of sex determination in a eukaryote and highlight Md1 as a potential target for interventions that block malaria transmission.

A Detailed Gene Expression Map of Giardia Encystation

Giardia intestinalis is an intestinal protozoan parasite that causes diarrheal infections worldwide. A key process to sustain its chain of transmission is the formation of infectious cysts in the encystation process. We combined deep RNAseq of a broad range of encystation timepoints to produce a high-resolution gene expression map of Giardia encystation. This detailed transcriptomic map of encystation confirmed a gradual change of gene expression along the time course of encystation, showing the most significant gene expression changes during late encystation. Few genes are differentially expressed early in encystation, but the major cyst wall proteins CWP-1 and -2 are highly up-regulated already after 3.5 h encystation. Several transcription factors are sequentially up-regulated throughout the process, but many up-regulated genes at 7, 10, and 14 h post-induction of encystation have binding sites in the upstream regions for the Myb2 transcription factor, suggesting that Myb2 is a master regulator of encystation. We observed major changes in gene expression of several meiotic-related genes from 10.5 h of encystation to the cyst stage, and at 17.5 h encystation, there are changes in many different metabolic pathways and protein synthesis. Late encystation, 21 h to cysts, show extensive gene expression changes, most of all in VSP and HCMP genes, which are involved in antigenic variation, and genes involved in chromatin modifications. This high-resolution gene expression map of Giardia encystation will be an important tool in further studies of this important differentiation process.

Global diversity of the gene encoding the Pfs25 protein-a Plasmodium falciparum transmission-blocking vaccine candidate

Background

Vaccines against the sexual stages of the malarial parasite Plasmodium falciparum are indispensable for controlling malaria and abrogating the spread of drug-resistant parasites. Pfs25, a surface antigen of the sexual stage of P. falciparum, is a leading candidate for transmission-blocking vaccine development. While clinical trials have reported that Pfs25-based vaccines are safe and effective in inducing transmission-blocking antibodies, the extent of the genetic diversity of Pfs25 in malaria endemic populations has rarely been studied. Thus, this study aimed to investigate the global diversity of Pfs25 in P. falciparum populations.

Methods

A database of 307 Pfs25 sequences of P. falciparum was established. Population genetic analyses were performed to evaluate haplotype and nucleotide diversity, analyze haplotypic distribution patterns of Pfs25 in different geographical populations, and construct a haplotype network. Neutrality tests were conducted to determine evidence of natural selection. Homology models of the Pfs25 haplotypes were constructed, subjected to molecular dynamics (MD), and analyzed in terms of flexibility and percentages of secondary structures.

Results

The Pfs25 gene of P. falciparum was found to have 11 unique haplotypes. Of these, haplotype 1 (H1) and H2, the major haplotypes, represented 70% and 22% of the population, respectively, and were dominant in Asia, whereas only H1 was dominant in Africa, Central America, and South America. Other haplotypes were rare and region-specific, resulting in unique distribution patterns in different geographical populations. The diversity in Pfs25 originated from ten single-nucleotide polymorphism (SNP) loci located in the epidermal growth factor (EGF)-like domains and anchor domain. Of these, an SNP at position 392 (GGA/GCA), resulting in amino acid substitution 131 (Gly/Ala), defined the two major haplotypes. The MD results showed that the structures of H1 and H2 variants were relatively similar. Limited polymorphism in Pfs25 could likely be due to negative selection.

Conclusions

The study successfully established a Pfs25 sequence database that can become an essential tool for monitoring vaccine efficacy, designing assays for detecting malaria carriers, and conducting epidemiological studies of P. falciparum. The discovery of the two major haplotypes, H1 and H2, and their conserved structures suggests that the current Pfs25-based vaccines could be used globally for malaria control.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05078-6.

Chemogenomic Fingerprints Associated with Stage-Specific Gametocytocidal Compound Action against Human Malaria Parasites

Kinase-focused inhibitors previously revealed compounds with differential activity against different stages of Plasmodium falciparum gametocytes. MMV666810, a 2-aminopyrazine, is more active on late-stage gametocytes, while a pyrazolopyridine, MMV674850, preferentially targets early-stage gametocytes. Here, we probe the biological mechanisms underpinning this differential stage-specific killing using in-depth transcriptome fingerprinting. Compound-specific chemogenomic profiles were observed with MMV674850 treatment associated with biological processes shared between asexual blood stage parasites and early-stage gametocytes but not late-stage gametocytes. MMV666810 has a distinct profile with clustered gene sets associated primarily with late-stage gametocyte development, including Ca²⁺-dependent protein kinases (CDPK1 and 5) and serine/threonine protein kinases (FIKK). Chemogenomic profiling therefore highlights essential processes in late-stage gametocytes, on a transcriptional level. This information is important to prioritize compounds that preferentially compromise late-stage gametocytes for further development as transmission-blocking antimalarials.

CRISPR/Cas9-engineered inducible gametocyte producer lines as a valuable tool for Plasmodium falciparum malaria transmission research

The malaria parasite Plasmodium falciparum replicates inside erythrocytes in the blood of infected humans. During each replication cycle, a small proportion of parasites commits to sexual development and differentiates into gametocytes, which are essential for parasite transmission via the mosquito vector. Detailed molecular investigation of gametocyte biology and transmission has been hampered by difficulties in generating large numbers of these highly specialised cells. Here, we engineer P. falciparum NF54 inducible gametocyte producer (iGP) lines for the routine mass production of synchronous gametocytes via conditional overexpression of the sexual commitment factor GDV1. NF54/iGP lines consistently achieve sexual commitment rates of 75% and produce viable gametocytes that are transmissible by mosquitoes. We also demonstrate that further genetic engineering of NF54/iGP parasites is a valuable tool for the targeted exploration of gametocyte biology. In summary, we believe the iGP approach developed here will greatly expedite basic and applied malaria transmission stage research.

A cascade of transcriptional repression determines sexual commitment and development in Plasmodium falciparum

Gametocytogenesis, the process by which malaria parasites produce sexual forms that can infect mosquitoes, is essential for the transmission of malaria. A transcriptional switch of the pfap2-g gene triggers sexual commitment, but how the complex multi-step process is precisely programed remains largely unknown. Here, by systematic functional screening of a panel of ApiAP2 transcription factors, we identify six new ApiAP2 members associated with gametocytogenesis in Plasmodium falciparum. Among these, PfAP2-G5 (PF3D7_1139300) was found to be indispensable for gametocytogenesis. This factor suppresses the transcriptional activity of the pfap2-g gene via binding to both the upstream region and exonic gene body, the latter is linked to the maintenance of local heterochromatin structure, thereby preventing initiation of sexual commitment. Removal of this repressive effect through pfap2-g5 knockout disrupts the asexual replication cycle and promotes sexual commitment accompanied by upregulation of pfap2-g expression. However, the gametocytes produced fail to mature fully. Further analyses show that PfAP2-G5 is essential for gametocyte maturation, and causes the down-regulation of pfap2-g and a set of early gametocyte genes activated by PfAP2-G prior to gametocyte development. Collectively, our findings reveal a regulation cascade of gametocyte production in malaria parasites, and provide a new target for transmission blocking interventions.

MALBoost: a web-based application for gene regulatory network analysis in Plasmodium falciparum

Background

Gene Regulatory Networks (GRN) produce powerful insights into transcriptional regulation in cells. The power of GRNs has been underutilized in malaria research. The Arboreto library was incorporated into a user-friendly web-based application for malaria researchers (http://malboost.bi.up.ac.za). This application will assist researchers with gaining an in depth understanding of transcriptomic datasets.

Methods

The web application for MALBoost was built in Python-Flask with Redis and Celery workers for queue submission handling, which execute the Arboreto suite algorithms. A submission of 5–50 regulators and total expression set of 5200 genes is permitted. The program runs in a point-and-click web user interface built using Bootstrap4 templates. Post-analysis submission, users are redirected to a status page with run time estimates and ultimately a download button upon completion. Result updates or failure updates will be emailed to the users.

Results

A web-based application with an easy-to-use interface is presented with a use case validation of AP2-G and AP2-I. The validation set incorporates cross-referencing with ChIP-seq and transcriptome datasets. For AP2-G, 5 ChIP-seq targets were significantly enriched with seven more targets presenting with strong evidence of validated targets.

Conclusion

The MALBoost application provides the first tool for easy interfacing and efficiently allows gene regulatory network construction for Plasmodium. Additionally, access is provided to a pre-compiled network for use as reference framework. Validation for sexually committed ring-stage parasite targets of AP2-G, suggests the algorithm was effective in resolving “traditionally” low-level signatures even in bulk RNA datasets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03848-2.

Benzimidazole Derivatives Are Potent against Multiple Life Cycle Stages of Plasmodium falciparum Malaria Parasites

The continued emergence of resistance to front-line antimalarial treatments is of great concern. Therefore, new compounds that potentially have a novel target in various developmental stages of Plasmodium parasites are needed to treat patients and halt the spread of malaria. Here, several benzimidazole derivatives were screened for activity against the symptom-causing intraerythrocytic asexual blood stages and the transmissible gametocyte stages of P. falciparum. Submicromolar activity was obtained for 54 compounds against asexual blood stage parasites with 6 potent at IC₅₀ < 100 nM while not displaying any marked toxicity against mammalian cells. Nanomolar potency was also observed against gametocytes with two compounds active against early stage gametocytes and two compounds active against late-stage gametocytes. The transmission-blocking potential of the latter was confirmed as they could prevent male gamete exflagellation and the lead compound reduced transmission by 72% in an in vivo mosquito feeding model. These compounds therefore have activity against multiple stages of Plasmodium parasites with potential for differential targets.

Machine Learning Uses Chemo-Transcriptomic Profiles to Stratify Antimalarial Compounds With Similar Mode of Action

The rapid development of antimalarial resistance motivates the continued search for novel compounds with a mode of action (MoA) different to current antimalarials. Phenotypic screening has delivered thousands of promising hit compounds without prior knowledge of the compounds’ exact target or MoA. Whilst the latter is not initially required to progress a compound in a medicinal chemistry program, identifying the MoA early can accelerate hit prioritization, hit-to-lead optimization and preclinical combination studies in malaria research. The effects of drug treatment on a cell can be observed on systems level in changes in the transcriptome, proteome and metabolome. Machine learning (ML) algorithms are powerful tools able to deconvolute such complex chemically-induced transcriptional signatures to identify pathways on which a compound act and in this manner provide an indication of the MoA of a compound. In this study, we assessed different ML approaches for their ability to stratify antimalarial compounds based on varied chemically-induced transcriptional responses. We developed a rational gene selection approach that could identify predictive features for MoA to train and generate ML models. The best performing model could stratify compounds with similar MoA with a classification accuracy of 76.6 ± 6.4%. Moreover, only a limited set of 50 biomarkers was required to stratify compounds with similar MoA and define chemo-transcriptomic fingerprints for each compound. These fingerprints were unique for each compound and compounds with similar targets/MoA clustered together. The ML model was specific and sensitive enough to group new compounds into MoAs associated with their predicted target and was robust enough to be extended to also generate chemo-transcriptomic fingerprints for additional life cycle stages like immature gametocytes. This work therefore contributes a new strategy to rapidly, specifically and sensitively indicate the MoA of compounds based on chemo-transcriptomic fingerprints and holds promise to accelerate antimalarial drug discovery programs.

Recent achievements and doors opened for coccidian parasite research and development through transcriptomics of enteric sexual stages

The Coccidia is the largest group of parasites within the Apicomplexa, a phylum of unicellular, obligate parasites characterized by the possession of an apical complex of organelles and structures in the asexual stages of their life cycles, as well as by a sexual reproductive phase that occurs enterically in host animals. Coccidian sexual reproduction involves morphologically distinct microgametes and macrogametes that combine to form a diploid zygote and, ultimately, following meiosis and mitosis, haploid, infectious sporozoites, inside sporocysts within an oocyst. Recent transcriptomic analyses have identified genes involved in coccidian sexual stage development and reproduction, including genes encoding for microgamete- and macrogamete-specific proteins with roles in gamete motility, fusion and fertilization, and in the formation of the resilient oocyst wall that allows coccidians to persist for long periods in the environment. Transcriptomics has also provided important clues about the regulation of gene expression in the transformation of parasites from one developmental stage to the next, a complex sequence of events that may involve transcription factors such as the apicomplexan Apetala2 (ApiAP2) family, alternative splicing, regulatory RNAs and MORC (a microrchida homologue and regulator of sexual stage development in Toxoplasma gondii). The molecular dissection of coccidian sexual development and reproduction by transcriptomic analyses may lead to the development of novel transmission-blocking strategies.

The PfAP2-G2 transcription factor is a critical regulator of gametocyte maturation

Differentiation from asexual blood stages to mature sexual gametocytes is required for the transmission of malaria parasites. Here, we report that the ApiAP2 transcription factor, PfAP2-G2 (PF3D7_1408200) plays a critical role in the maturation of Plasmodium falciparum gametocytes. PfAP2-G2 binds to the promoters of a wide array of genes that are expressed at many stages of the parasite life cycle. Interestingly, we also find binding of PfAP2-G2 within the gene body of almost 3,000 genes, which strongly correlates with the location of H3K36me3 and several other histone modifications as well as Heterochromatin Protein 1 (HP1), suggesting that occupancy of PfAP2-G2 in gene bodies may serve as an alternative regulatory mechanism. Disruption of pfap2-g2 does not impact asexual development, but the majority of sexual parasites are unable to mature beyond stage III gametocytes. The absence of pfap2-g2 leads to overexpression of 28% of the genes bound by PfAP2-G2 and none of the PfAP2-G2 bound genes are downregulated, suggesting that it is a repressor. We also find that PfAP2-G2 interacts with chromatin remodeling proteins, a microrchidia (MORC) protein, and another ApiAP2 protein (PF3D7_1139300). Overall our data demonstrate that PfAP2-G2 establishes an essential gametocyte maturation program in association with other chromatin-related proteins.

H3K36 methylation reprograms gene expression to drive early gametocyte development in Plasmodium falciparum

Background

The Plasmodium sexual gametocyte stages are the only transmissible form of the malaria parasite and are thus responsible for the continued transmission of the disease. Gametocytes undergo extensive functional and morphological changes from commitment to maturity, directed by an equally extensive control program. However, the processes that drive the differentiation and development of the gametocyte post-commitment, remain largely unexplored. A previous study reported enrichment of H3K36 di- and tri-methylated (H3K36me2&3) histones in early-stage gametocytes. Using chromatin immunoprecipitation followed by high-throughput sequencing, we identify a stage-specific association between these repressive histone modifications and transcriptional reprogramming that define a stage II gametocyte transition point.

Results

Here, we show that H3K36me2 and H3K36me3 from stage II gametocytes are associated with repression of genes involved in asexual proliferation and sexual commitment, indicating that H3K36me2&3-mediated repression of such genes is essential to the transition from early gametocyte differentiation to intermediate development. Importantly, we show that the gene encoding the transcription factor AP2-G as commitment master regulator is enriched with H3K36me2&3 and actively repressed in stage II gametocytes, providing the first evidence of ap2-g gene repression in post-commitment gametocytes. Lastly, we associate the enhanced potency of the pan-selective Jumonji inhibitor JIB-04 in gametocytes with the inhibition of histone demethylation including H3K36me2&3 and a disruption of normal transcriptional programs.

Conclusions

Taken together, our results provide the first description of an association between global gene expression reprogramming and histone post-translational modifications during P. falciparum early sexual development. The stage II gametocyte-specific abundance of H3K36me2&3 manifests predominantly as an independent regulatory mechanism targeted towards genes that are repressed post-commitment. H3K36me2&3-associated repression of genes is therefore involved in key transcriptional shifts that accompany the transition from early gametocyte differentiation to intermediate development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00393-9.

Targeting Gametocytes of the Malaria Parasite Plasmodium falciparum in a Functional Genomics Era: Next Steps

Mosquito transmission of the deadly malaria parasite Plasmodium falciparum is mediated by mature sexual forms (gametocytes). Circulating in the vertebrate host, relatively few intraerythrocytic gametocytes are picked up during a bloodmeal to continue sexual development in the mosquito vector. Human-to-vector transmission thus represents an infection bottleneck in the parasite's life cycle for therapeutic interventions to prevent malaria. Even though recent progress has been made in the identification of genetic factors linked to gametocytogenesis, a plethora of genes essential for sexual-stage development are yet to be unraveled. In this review, we revisit P. falciparum transmission biology by discussing targetable features of gametocytes and provide a perspective on a forward-genetic approach for identification of novel transmission-blocking candidates in the future.

Multistage and transmission-blocking targeted antimalarials discovered from the open-source MMV Pandemic Response Box

Chemical matter is needed to target the divergent biology associated with the different life cycle stages of Plasmodium. Here, we report the parallel de novo screening of the Medicines for Malaria Venture (MMV) Pandemic Response Box against Plasmodium asexual and liver stage parasites, stage IV/V gametocytes, gametes, oocysts and as endectocides. Unique chemotypes were identified with both multistage activity or stage-specific activity, including structurally diverse gametocyte-targeted compounds with potent transmission-blocking activity, such as the JmjC inhibitor ML324 and the antitubercular clinical candidate SQ109. Mechanistic investigations prove that ML324 prevents histone demethylation, resulting in aberrant gene expression and death in gametocytes. Moreover, the selection of parasites resistant to SQ109 implicates the druggable V-type H⁺-ATPase for the reduced sensitivity. Our data therefore provides an expansive dataset of compounds that could be redirected for antimalarial development and also point towards proteins that can be targeted in multiple parasite life cycle stages.

Here, Reader et al. screen the Medicines for Malaria Venture Pandemic Response Box in parallel against Plasmodiumasexual and liver stage parasites, stage IV/V gametocytes, gametes, oocysts and as endectocides. They identify two potent transmission-blocking drugs: a histone demethylase inhibitor ML324 and the antitubercular SQ109.

PfMAP-2 is essential for male gametogenesis in the malaria parasite Plasmodium falciparum

In malaria parasites, male gametogenesis is a proliferative stage essential for parasite transmission to the mosquito vector. It is a rapid process involving three rounds of genome replication alternating with closed endomitoses, and assembly of axonemes to produce eight flagellated motile microgametes. Studies in Plasmodium berghei have highlighted tight regulation of gametogenesis by a network of kinases. The P. berghei MAPK homologue PbMAP-2 is dispensable for asexual development but important at the induction of axoneme motility. However, in P. falciparum, causing the most severe form of human malaria, PfMAP-2 was suggested to be essential for asexual proliferation indicating distinct functions for MAP-2 in these two Plasmodium species. We here show that PfMAP-2 is dispensable for asexual growth but important for male gametogenesis in vitro. Similar to PbMAP-2, PfMAP-2 is required for initiating axonemal beating but not for prior DNA replication or axoneme formation. In addition, single and double null mutants of PfMAP-2 and the second P. falciparum MAPK homologue PfMAP-1 show no defect in asexual proliferation, sexual commitment or gametocytogenesis. Our results suggest that MAPK activity plays no major role in the biology of both asexual and sexual blood stage parasites up until the point of male gametogenesis.