Identification and genome-wide prediction of DNA binding specificities for the ApiAP2 family of regulators from the malaria parasite

DNA sequence and chromatin differentiate sequence-specific transcription factor binding in the human malaria parasite Plasmodium falciparum

Development of the malaria parasite, Plasmodium falciparum, is regulated by a limited number of sequence-specific transcription factors (TFs). However, the mechanisms by which these TFs recognize genome-wide binding sites is largely unknown. To address TF specificity, we investigated the binding of two TF subsets that either bind CACACA or GTGCAC DNA sequence motifs and further characterized two additional ApiAP2 TFs, PfAP2-G and PfAP2-EXP, which bind unique DNA motifs (GTAC and TGCATGCA). We also interrogated the impact of DNA sequence and chromatin context on P. falciparum TF binding by integrating high-throughput in vitro and in vivo binding assays, DNA shape predictions, epigenetic post-translational modifications, and chromatin accessibility. We found that DNA sequence context minimally impacts binding site selection for paralogous CACACA-binding TFs, while chromatin accessibility, epigenetic patterns, co-factor recruitment, and dimerization correlate with differential binding. In contrast, GTGCAC-binding TFs prefer different DNA sequence context in addition to chromatin dynamics. Finally, we determined that TFs that preferentially bind divergent DNA motifs may bind overlapping genomic regions due to low-affinity binding to other sequence motifs. Our results demonstrate that TF binding site selection relies on a combination of DNA sequence and chromatin features, thereby contributing to the complexity of P. falciparum gene regulatory mechanisms.

An ApiAp2 Transcription Factor with a Dispensable Role in Plasmodium berghei Life Cycle

Malaria parasites have a complex life cycle and undergo replication and population expansion within vertebrate hosts and mosquito vectors. These developmental transitions rely on changes in gene expression and chromatin reorganization that result in the activation and silencing of stage-specific genes. The ApiAp2 family of DNA-binding proteins plays an important role in regulating gene expression in malaria parasites. Here, we characterized the ApiAp2 protein in Plasmodium berghei, which we termed Ap2-D. In silico analysis revealed that Ap2-D has three beta-sheets followed by a helix at the C-terminus for DNA binding. Using gene tagging with 3XHA-mCherry, we found that Ap2-D is expressed in Plasmodium blood stages and is present in the parasite cytoplasm and nucleus. Surprisingly, our gene deletion study revealed a completely dispensable role for Ap2-D in the entirety of the P. berghei life cycle. Ap2-D KO parasites were found to grow in the blood successfully and progress through the mosquito midgut and salivary glands. Sporozoites isolated from mosquito salivary glands were infective for hepatocytes and achieved similar patency as WT in mice. We emphasize the importance of genetic validation of antimalarial drug targets before progressing them to drug discovery.

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.

Ready for renascence in mosquito: The regulation of gene expression in Plasmodium sexual development

Malaria remains one of the most perilous vector-borne infectious diseases for humans globally. Sexual gametocyte represents the exclusive stage at which malaria parasites are transmitted from the vertebrate to the Anopheles host. The feasible and effective approach to prevent malaria transmission is by addressing the sexual developmental processes, that is, gametocytogenesis and gametogenesis. Thus, this review will comprehensively cover advances in the regulation of gene expression surrounding the transmissible stages, including epigenetic, transcriptional, and post-transcriptional control.

Coordinated regulation of gene expression in Plasmodium female gametocytes by two transcription factors

Gametocytes play key roles in the Plasmodium lifecycle. They are essential for sexual reproduction as precursors of the gametes. They also play an essential role in parasite transmission to mosquitoes. Elucidation of the gene regulation at this stage is essential for understanding these two processes at the molecular level and for developing new strategies to break the parasite lifecycle. We identified a novel Plasmodium transcription factor (TF), designated as a partner of AP2-FG or PFG. In this article, we report that this TF regulates the gene expression in female gametocytes in concert with another female-specific TF AP2-FG. Upon the disruption of PFG, majority of female-specific genes were significantly downregulated, and female gametocyte lost the ability to produce ookinetes. ChIP-seq analysis showed that it was located in the same position as AP2-FG, indicating that these two TFs form a complex. ChIP-seq analysis of PFG in AP2-FG-disrupted parasites and ChIP-seq analysis of AP2-FG in PFG-disrupted parasites demonstrated that PFG mediates the binding of AP2-FG to a ten-base motif and that AP2-FG binds another motif, GCTCA, in the absence of PFG. In promoter assays, this five-base motif was identified as another female-specific cis-acting element. Genes under the control of the two forms of AP2-FG, with or without PFG, partly overlapped; however, each form had target preferences. These results suggested that combinations of these two forms generate various expression patterns among the extensive genes expressed in female gametocytes.

Localization in vivo and in vitro confirms EnApiAP2 protein encoded by ENH_00027130 as a nuclear protein in Eimeria necatrix

Introduction

Apicomplexan AP2 family of proteins (ApiAP2) are transcription factors (TFs) that regulate parasite growth and development, but little is known about the ApiAP2 TFs in Eimeria spp. ENH_00027130 sequence is predicted to encode a Eimeria necatrix ApiAP2 protein (EnApiAP2).

Methods

The cDNAs encoding full-length and truncated EnApiAP2 protein were cloned and sequenced, respectively. Then, the two cDNAs were cloned into the pET28a(+) expression vector and expressed expressed in Escherichia coli BL21. The mouse polyclonal antibody (pAb) and monoclonal antibody (mAb) against recombinant EnApiAP2 (rEnApiAP2) and EnApiAP2tr (rEnApiAP2tr) were prepared and used to localize the native EnApiAP2 protein in E. necatrix, respectively. Finally, the recombinant pEGFP-C1-ΔNLS-EnApiAP2s (knockout of a nuclear localization sequence, NLS) and pEGFP-C1-EnApiAP2 plasmid were constructed and transfected into DF-1 cells, respectively, to further observe subcellular localization of EnApiAP2 protein.

Results

The EnApiAP2 gene had a size of 5019 bp and encoded 1672 amino acids, containing a conserved AP2 domain with a secondary structure consisting of an α-helix and three antiparallel β-strands. The rEnApiAP2 and rEnApiAP2tr were predominantly expressed in the form of inclusion bodies, and could be recognized by the 6×His tag mAb and the serum of convalescent chickens after infection with E. necatrix, respectively. The native EnApiAP2 protein was detected in sporozoites (SZ) and second generation merozoites (MZ-2) extracts, with a size of approximately 210 kDa. A quantitative real-time PCR (qPCR) analysis showed that the transcription level of EnApiAP2 was significantly higher in SZ than in MZ-2, third generation merozoites (MZ-3) and gametocytes (P<0.01). EnApiAP2 protein was localized in the nuclei of SZ, MZ-2 and MZ-3 of E. necatrix. The protein of EnApiAP2 was localized in the nucleus of the DF-1 cells, whereas the ΔNLS-EnApiAP2 was expressed in the cytoplasm, which further confirmed that EnApiAP2 is nucleoprotein.

Discussion

EnApiAP2 protein encoded by ENH_00027130 sequence was localized in the nucleus of E. necatrix parasites, and relied on the NLS for migration to DF-1 cell nucleus. The function of EnApiAP2 need further study.

Identifying unique exposure-specific transgenerational differentially DNA methylated region epimutations in the genome using hybrid deep learning prediction models

Exposure to environmental toxicants can lead to epimutations in the genome and an increase in differential DNA methylated regions (DMRs) that have been linked to increased susceptibility to various diseases. However, the unique effect of particular toxicants on the genome in terms of leading to unique DMRs for the toxicants has been less studied. One hurdle to such studies is the low number of observed DMRs per toxicants. To address this hurdle, a previously validated hybrid deep-learning cross-exposure prediction model is trained per exposure and used to predict exposure-specific DMRs in the genome. Given these predicted exposure-specific DMRs, a set of unique DMRs per exposure can be identified. Analysis of these unique DMRs through visualization, DNA sequence motif matching, and gene association reveals known and unknown links between individual exposures and their unique effects on the genome. The results indicate the potential ability to define exposure-specific epigenetic markers in the genome and the potential relative impact of different exposures. Therefore, a computational approach to predict exposure-specific transgenerational epimutations was developed, which supported the exposure specificity of ancestral toxicant actions and provided epigenome information on the DMR sites predicted.

DNA-binding protein PfAP2-P regulates parasite pathogenesis during malaria parasite blood stages

Malaria-associated pathogenesis such as parasite invasion, egress, host cell remodelling and antigenic variation requires concerted action by many proteins, but the molecular regulation is poorly understood. Here we have characterized an essential Plasmodium-specific Apicomplexan AP2 transcription factor in Plasmodium falciparum (PfAP2-P; pathogenesis) during the blood-stage development with two peaks of expression. An inducible knockout of gene function showed that PfAP2-P is essential for trophozoite development, and critical for var gene regulation, merozoite development and parasite egress. Chromatin immunoprecipitation sequencing data collected at timepoints matching the two peaks of pfap2-p expression demonstrate PfAP2-P binding to promoters of genes controlling trophozoite development, host cell remodelling, antigenic variation and pathogenicity. Single-cell RNA sequencing and fluorescence-activated cell sorting revealed de-repression of most var genes in Δpfap2-p parasites. Δpfap2-p parasites also overexpress early gametocyte marker genes, indicating a regulatory role in sexual stage conversion. We conclude that PfAP2-P is an essential upstream transcriptional regulator at two distinct stages of the intra-erythrocytic development cycle.

MalariaSED: a deep learning framework to decipher the regulatory contributions of noncoding variants in malaria parasites

Malaria remains one of the deadliest infectious diseases. Transcriptional regulation effects of noncoding variants in this unusual genome of malaria parasites remain elusive. We developed a sequence-based, ab initio deep learning framework, MalariaSED, for predicting chromatin profiles in malaria parasites. The MalariaSED performance was validated by published ChIP-qPCR and TF motifs results. Applying MalariaSED to ~ 1.3 million variants shows that geographically differentiated noncoding variants are associated with parasite invasion and drug resistance. Further analysis reveals chromatin accessibility changes at Plasmodium falciparum rings are partly associated with artemisinin resistance. MalariaSED illuminates the potential functional roles of noncoding variants in malaria parasites.

Cohesin contributes to transcriptional repression of stage-specific genes in the human malaria parasite

The complex life cycle of the human malaria parasite, Plasmodium falciparum, is driven by specific transcriptional programs, but it is unclear how most genes are activated or silenced at specific times. There is an association between transcription and spatial organization; however, the molecular mechanisms behind genome organization are unclear. While P. falciparum lacks key genome-organizing proteins found in metazoans, it has all core components of the cohesin complex. To investigate the role of cohesin in P. falciparum, we functionally characterize the cohesin subunit Structural Maintenance of Chromosomes protein 3 (SMC3). SMC3 knockdown during early stages of the intraerythrocytic developmental cycle (IDC) upregulates a subset of genes involved in erythrocyte egress and invasion, which are normally expressed at later stages. ChIP-seq analyses reveal that during the IDC, SMC3 enrichment at the promoter regions of these genes inversely correlates with gene expression and chromatin accessibility. These data suggest that SMC3 binding contributes to the repression of specific genes until their appropriate time of expression, revealing a new mode of stage-specific gene repression in P. falciparum.

Epigenetic Regulation and Chromatin Remodeling in Malaria Parasites

Plasmodium falciparum, the human malaria parasite, infects two hosts and various cell types, inducing distinct morphological and physiological changes in the parasite in response to different environmental conditions. These variations required the parasite to adapt and develop elaborate molecular mechanisms to ensure its spread and transmission. Recent findings have significantly improved our understanding of the regulation of gene expression in P. falciparum. Here, we provide an up-to-date overview of technologies used to highlight the transcriptomic adjustments occurring in the parasite throughout its life cycle. We also emphasize the complementary and complex epigenetic mechanisms regulating gene expression in malaria parasites. This review concludes with an outlook on the chromatin architecture, the remodeling systems, and how this 3D genome organization is critical in various biological processes.

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.

Plasmodium falciparum Development from Gametocyte to Oocyst: Insight from Functional Studies

Malaria elimination may never succeed without the implementation of transmission-blocking strategies. The transmission of Plasmodium spp. parasites from the human host to the mosquito vector depends on circulating gametocytes in the peripheral blood of the vertebrate host. Once ingested by the mosquito during blood meals, these sexual forms undergo a series of radical morphological and metabolic changes to survive and progress from the gut to the salivary glands, where they will be waiting to be injected into the vertebrate host. The design of effective transmission-blocking strategies requires a thorough understanding of all the mechanisms that drive the development of gametocytes, gametes, sexual reproduction, and subsequent differentiation within the mosquito. The drastic changes in Plasmodium falciparum shape and function throughout its life cycle rely on the tight regulation of stage-specific gene expression. This review outlines the mechanisms involved in Plasmodium falciparum sexual stage development in both the human and mosquito vector, and zygote to oocyst differentiation. Functional studies unravel mechanisms employed by P. falciparum to orchestrate the expression of stage-specific functional products required to succeed in its complex life cycle, thus providing us with potential targets for developing new therapeutics. These mechanisms are based on studies conducted with various Plasmodium species, including predominantly P. falciparum and the rodent malaria parasites P. berghei. However, the great potential of epigenetics, genomics, transcriptomics, proteomics, and functional genetic studies to improve the understanding of malaria as a disease remains partly untapped because of limitations in studies using human malaria parasites and field isolates.

Histone modification analysis reveals common regulators of gene expression in liver and blood stage merozoites of Plasmodium parasites

Gene expression in malaria parasites is subject to various layers of regulation, including histone post-translational modifications (PTMs). Gene regulatory mechanisms have been extensively studied during the main developmental stages of Plasmodium parasites inside erythrocytes, from the ring stage following invasion to the schizont stage leading up to egress. However, gene regulation in merozoites that mediate the transition from one host cell to the next is an understudied area of parasite biology. Here, we sought to characterize gene expression and the corresponding histone PTM landscape during this stage of the parasite lifecycle through RNA-seq and ChIP-seq on P. falciparum blood stage schizonts, merozoites, and rings, as well as P. berghei liver stage merozoites. In both hepatic and erythrocytic merozoites, we identified a subset of genes with a unique histone PTM profile characterized by a region of H3K4me3 depletion in their promoter. These genes were upregulated in hepatic and erythrocytic merozoites and rings, had roles in protein export, translation, and host cell remodeling, and shared a DNA motif. These results indicate that similar regulatory mechanisms may underlie merozoite formation in the liver and blood stages. We also observed that H3K4me2 was deposited in gene bodies of gene families encoding variant surface antigens in erythrocytic merozoites, which may facilitate switching of gene expression between different members of these families. Finally, H3K18me and H2K27me were uncoupled from gene expression and were enriched around the centromeres in erythrocytic schizonts and merozoites, suggesting potential roles in the maintenance of chromosomal organization during schizogony. Together, our results demonstrate that extensive changes in gene expression and histone landscape occur during the schizont-to-ring transition to facilitate productive erythrocyte infection. The dynamic remodeling of the transcriptional program in hepatic and erythrocytic merozoites makes this stage attractive as a target for novel anti-malarial drugs that may have activity against both the liver and blood stages.

PfAP2-MRP DNA-binding protein is a master regulator of parasite pathogenesis during malaria parasite blood stages

Malaria pathogenicity results from the parasite's ability to invade, multiply within and then egress from the host red blood cell (RBC). Infected RBCs are remodeled, expressing antigenic variant proteins (such as PfEMP1, coded by the var gene family) for immune evasion and survival. These processes require the concerted actions of many proteins, but the molecular regulation is poorly understood. We have characterized an essential Plasmodium specific Apicomplexan AP2 (ApiAP2) transcription factor in Plasmodium falciparum (PfAP2-MRP; Master Regulator of Pathogenesis) during the intraerythrocytic developmental cycle (IDC). An inducible gene knockout approach showed that PfAP2-MRP is essential for development during the trophozoite stage, and critical for var gene regulation, merozoite development and parasite egress. ChIP-seq experiments performed at 16 hour post invasion (h.p.i.) and 40 h.p.i. matching the two peaks of PfAP2-MRP expression, demonstrate binding of PfAP2-MRP to the promoters of genes controlling trophozoite development and host cell remodeling at 16 h.p.i. and antigenic variation and pathogenicity at 40 h.p.i. Using single-cell RNA-seq and fluorescence-activated cell sorting, we show de-repression of most var genes in Δpfap2-mrp parasites that express multiple PfEMP1 proteins on the surface of infected RBCs. In addition, the Δpfap2-mrp parasites overexpress several early gametocyte marker genes at both 16 and 40 h.p.i., indicating a regulatory role in the sexual stage conversion. Using the Chromosomes Conformation Capture experiment (Hi-C), we demonstrate that deletion of PfAP2-MRP results in significant reduction of both intra-chromosomal and inter-chromosomal interactions in heterochromatin clusters. We conclude that PfAP2-MRP is a vital upstream transcriptional regulator controlling essential processes in two distinct developmental stages during the IDC that include parasite growth, chromatin structure and var gene expression.

Toxoplasma gondii AP2XII-2 Contributes to Transcriptional Repression for Sexual Commitment

Toxoplasma gondii is a widespread protozoan parasite that has a significant impact on human and veterinary health. The parasite undergoes a complex life cycle involving multiple hosts and developmental stages. How Toxoplasma transitions between life cycle stages is poorly understood yet central to controlling transmission. Of particular neglect are the factors that contribute to its sexual development, which takes place exclusively in feline intestines. While epigenetic repressors have been shown to play an important role in silencing the spurious gene expression of sexually committed parasites, the specific factors that recruit this generalized machinery to the appropriate genes remain largely unexplored. Here, we establish that a member of the AP2 transcription factor family, AP2XII-2, is targeted to genomic loci associated with sexually committed parasites along with epigenetic regulators of transcriptional silencing, HDAC3 and MORC. Despite its widespread association with gene promoters, AP2XII-2 is required for the silencing of relatively few genes. Using the CUT&Tag (cleavage under targets and tagmentation) methodology, we identify two major genes associated with sexual development downstream of AP2XII-2 control, AP2X-10 and the amino acid hydroxylase AAH1. Our findings show that AP2XII-2 is a key contributor to the gene regulatory pathways modulating Toxoplasma sexual development. IMPORTANCE Toxoplasma gondii is a parasite that undergoes its sexual stage exclusively in feline intestines, making cats a major source of transmission. A better understanding of the proteins controlling the parasite's life cycle stage transitions is needed for the development of new therapies aimed at treating toxoplasmosis and the transmission of the infection. Genes that regulate the sexual stages need to be turned on and off at the appropriate times, activities that are mediated by specific transcription factors that recruit general machinery to silence or activate gene expression. In this study, we identify a transcription factor called AP2XII-2 as being important for the repression of a subset of sexual stage genes, including a sexual stage-specific AP2 factor (AP2X-10) and a protein (AAH1) required to construct the infectious oocysts expelled from infected cats.

In silico Structure Prediction, Molecular Docking, and Dynamic Simulation of Plasmodium falciparum AP2-I Transcription Factor

Plasmodium falciparum Apicomplexan Apetala 2 Invasion (PfAP2-I) transcription factor (TF) is a protein that regulates the expression of a subset of gene families involved in P. falciparum red blood cell (RBC) invasion. Inhibiting PfAP2-I TF with small molecules represents a potential new antimalarial therapeutic target to combat drug resistance, which this study aims to achieve. The 3D model structure of PfAP2-I was predicted ab initio using ROBETTA prediction tool and was validated using Save server 6.0 and MolProbity. Computed Atlas of Surface Topography of proteins (CASTp) 3.0 was used to predict the active sites of the PfAP2-I modeled structure. Pharmacophore modeling of the control ligand and PfAP2-I modeled structure was carried out using the Pharmit server to obtain several compounds used for molecular docking analysis. Molecular docking and postdocking studies were conducted using AutoDock vina and Discovery studio. The designed ligands' toxicity predictions and in silico drug-likeness were performed using the SwissADME predictor and OSIRIS Property Explorer. The modeled protein structure from the ROBETTA showed a validation result of 96.827 for ERRAT, 90.2% of the amino acid residues in the most favored region for the Ramachandran plot, and MolProbity score of 1.30 in the 98th percentile. Five (5) best hit compounds from molecular docking analysis were selected based on their binding affinity (between -8.9 and -11.7 Kcal/mol) to the active site of PfAP2-I and were considered for postdocking studies. For the absorption, distribution, metabolism, elimination, and toxicity (ADMET) properties, compound MCULE-7146940834 had the highest drug score (0.63) and drug-likeness (6.76). MCULE-7146940834 maintained a stable conformation within the flexible protein's active site during simulation. The good, estimated binding energies, drug-likeness, drug score, and molecular dynamics simulation interaction observed for MCULE-7146940834 against PfAP2-I show that MCULE-7146940834 can be considered a lead candidate for PfAP2-I inhibition. Experimental validations should be carried out to ascertain the efficacy of these predicted best hit compounds.

Inhibitors of ApiAP2 protein DNA binding exhibit multistage activity against Plasmodium parasites

Plasmodium parasites are reliant on the Apicomplexan AP2 (ApiAP2) transcription factor family to regulate gene expression programs. AP2 DNA binding domains have no homologs in the human or mosquito host genomes, making them potential antimalarial drug targets. Using an in-silico screen to dock thousands of small molecules into the crystal structure of the AP2-EXP (Pf3D7_1466400) AP2 domain (PDB:3IGM), we identified putative AP2-EXP interacting compounds. Four compounds were found to block DNA binding by AP2-EXP and at least one additional ApiAP2 protein. Our top ApiAP2 competitor compound perturbs the transcriptome of P. falciparum trophozoites and results in a decrease in abundance of log₂ fold change > 2 for 50% (46/93) of AP2-EXP target genes. Additionally, two ApiAP2 competitor compounds have multi-stage anti-Plasmodium activity against blood and mosquito stage parasites. In summary, we describe a novel set of antimalarial compounds that interact with AP2 DNA binding domains. These compounds may be used for future chemical genetic interrogation of ApiAP2 proteins or serve as starting points for a new class of antimalarial therapeutics.

Plasmodium parasites are the causative agent of malaria, which resulted in over 600,000 deaths in 2021. Due to resistance arising for every antimalarial therapeutic deployed to date, new drug targets and druggable pathways must be explored. To address this concern, we used a molecular docking screen to predict competitors of DNA binding by the parasite specific family of Apicomplexan AP2 (ApiAP2) transcription factor proteins for testing in vitro and in vivo. We find that ApiAP2 competing compounds have antimalarial activity consistent with the disruption of gene regulation. This work will further our understanding of both the biological role and targetability of parasite transcriptional regulation.

Epigenetic and Epitranscriptomic Gene Regulation in Plasmodium falciparum and How We Can Use It against Malaria

Malaria, caused by Plasmodium parasites, is still one of the biggest global health challenges. P. falciparum is the deadliest species to humans. In this review, we discuss how this parasite develops and adapts to the complex and heterogenous environments of its two hosts thanks to varied chromatin-associated and epigenetic mechanisms. First, one small family of transcription factors, the ApiAP2 proteins, functions as master regulators of spatio-temporal patterns of gene expression through the parasite life cycle. In addition, chromatin plasticity determines variable parasite cell phenotypes that link to parasite growth, virulence and transmission, enabling parasite adaptation within host conditions. In recent years, epitranscriptomics is emerging as a new regulatory layer of gene expression. We present evidence of the variety of tRNA and mRNA modifications that are being characterized in Plasmodium spp., and the dynamic changes in their abundance during parasite development and cell fate. We end up outlining that new biological systems, like the mosquito model, to decipher the unknowns about epigenetic mechanisms in vivo; and novel methodologies, to study the function of RNA modifications; are needed to discover the Achilles heel of the parasite. With this new knowledge, future strategies manipulating the epigenetics and epitranscriptomic machinery of the parasite have the potential of providing new weapons against malaria.

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.

Division and Transmission: Malaria Parasite Development in the Mosquito

The malaria parasite life cycle alternates between two hosts: a vertebrate and the female Anopheles mosquito vector. Cell division, proliferation, and invasion are essential for parasite development, transmission, and survival. Most research has focused on Plasmodium development in the vertebrate, which causes disease; however, knowledge of malaria parasite development in the mosquito (the sexual and transmission stages) is now rapidly accumulating, gathered largely through investigation of the rodent malaria model, with Plasmodium berghei. In this review, we discuss the seminal genome-wide screens that have uncovered key regulators of cell proliferation, invasion, and transmission during Plasmodium sexual development. Our focus is on the roles of transcription factors, reversible protein phosphorylation, and molecular motors. We also emphasize the still-unanswered important questions around key pathways in cell division during the vector transmission stages and how they may be targeted in future studies.

Nuclear Factor AP2X-4 Governs the Expression of Cell Cycle- and Life Stage-Regulated Genes and is Critical for Toxoplasma Growth

Toxoplasma gondii is a ubiquitous pathogen infecting one third of the world's population and diverse animals. It has a complex life cycle alternating among different developmental stages, which contributes to its transmission and pathogenesis. The parasite has a sophisticated gene regulation network that enables timely expression of genes at designated stages. However, little is known about the underlying regulatory mechanisms. Here, we identified an AP2 family transcription factor named TgAP2X-4, which was crucial for parasite growth during the acute infection stage. TgAP2X-4 deletion leads to reduced expression of many genes that are normally upregulated during the M phase of the cell cycle. These include genes that encode rhoptry neck proteins that are key for parasite invasion. As a result, the Δap2X-4 mutant displayed significantly decreased efficiency of host cell invasion. Transcriptomic analyses suggested that TgAP2X-4 also regulates a large group of genes that are typically induced during chronic infection, such as BAG1 and LDH2. Given the diverse impacts on gene expression, TgAP2X-4 inactivation results in severely impaired parasite growth, as well as drastic attenuation of parasite virulence and complete inability to form chronic infection. Therefore, TgAP2X-4 represents a candidate for antitoxoplasmic drug and vaccine designs. IMPORTANCE Toxoplasma gondii has a complicated gene regulation network that allows "just in time" expression of genes to cope with the physiological needs at each stage during the complex life cycle. However, how such regulation is achieved is largely unknown. Here, we identified a transcription factor named TgAP2X-4 that is critical for the growth and life cycle progression of the parasite. Detailed analyses found that TgAP2X-4 regulated the expression of many cell cycle-regulated genes, including a subset of rhoptry genes that were essential for the parasites to enter host cells. It also regulated the expression of many genes involved in the development of chronic infection. Because of the diverse impacts on gene expression, TgAP2X-4 inactivation caused reduced parasite growth in vitro and attenuated virulence in vivo. Therefore, it is a potential target for drug or vaccine designs against Toxoplasma infections.

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.

Plasmodium falciparum S-Adenosylmethionine Synthetase Is Essential for Parasite Survival through a Complex Interaction Network with Cytoplasmic and Nuclear Proteins

S-adenosylmethionine synthetase (SAMS) is a key enzyme for the synthesis of the lone methyl donor S-adenosyl methionine (SAM), which is involved in transmethylation reactions and hence required for cellular processes such as DNA, RNA, and histone methylation, but also polyamine biosynthesis and proteostasis. In the human malaria parasite Plasmodium falciparum, PfSAMS is encoded by a single gene and has been suggested to be crucial for malaria pathogenesis and transmission; however, to date, PfSAMS has not been fully characterized. To gain deeper insight into the function of PfSAMS, we generated a conditional gene knockdown (KD) using the glmS ribozyme system. We show that PfSAMS localizes to the cytoplasm and the nucleus of blood-stage parasites. PfSAMS-KD results in reduced histone methylation and leads to impaired intraerythrocytic growth and gametocyte development. To further determine the interaction network of PfSAMS, we performed a proximity-dependent biotin identification analysis. We identified a complex network of 1114 proteins involved in biological processes such as cell cycle control and DNA replication, or transcription, but also in phosphatidylcholine and polyamine biosynthesis and proteasome regulation. Our findings highlight the diverse roles of PfSAMS during intraerythrocytic growth and sexual stage development and emphasize that PfSAMS is a potential drug target.

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.

Sources of transcription variation in Plasmodium falciparum

Variation in transcript abundance can contribute to both short-term environmental response and long-term evolutionary adaptation. Most studies are designed to assess differences in mean transcription levels and do not consider other potentially important and confounding sources of transcriptional variation. Detailed quantification of variation sources will improve our ability to detect and identify the mechanisms that contribute to genome-wide transcription changes that underpin adaptive responses. To quantify innate levels of expression variation, we measured mRNA levels for more than 5000 genes in the malaria parasite, Plasmodium falciparum, among clones derived from two parasite strains across biologically and experimentally replicated batches. Using a mixed effects model, we partitioned the total variation among four sources - between strain, within strain, environmental batch effects, and stochastic noise. We found 646 genes with significant variation attributable to at least one of these sources. These genes were categorized by their predominant variation source and further examined using gene ontology enrichment analysis to associate function with each source of variation. Genes with environmental batch effect and within strain transcript variation may contribute to phenotypic plasticity, while genes with between strain variation may contribute to adaptive responses and processes that lead to parasite strain-specific survival under varied conditions.

Malaria parasites do respond to heat

The capacity of malaria parasites to respond to changes in their environment at the transcriptional level has been the subject of debate, but recent evidence has unambiguously demonstrated that Plasmodium spp. can produce adaptive transcriptional responses when exposed to some specific types of stress. These include metabolic conditions and febrile temperature. The Plasmodium falciparum protective response to thermal stress is similar to the response in other organisms, but it is regulated by a transcription factor evolutionarily unrelated to the conserved transcription factor that drives the heat shock (HS) response in most eukaryotes. Of the many genes that change expression during HS, only a subset constitutes an authentic response that contributes to parasite survival.

PfAP2-EXP2, an Essential Transcription Factor for the Intraerythrocytic Development of Plasmodium falciparum

Plasmodium falciparum undergoes a series of asexual replications in human erythrocytes after infection, which are effective targets for combatting malaria. Here, we report roles of an ApiAP2 transcription factor PfAP2-EXP2 (PF3D7_0611200) in the intraerythrocytic developmental cycle of P. falciparum. PfAP2-EXP2 conditional knockdown resulted in an asexual growth defect but without an appreciable effect on parasite morphology. Further ChIP-seq analysis revealed that PfAP2-EXP2 targeted genes related to virulence and interaction between erythrocytes and parasites. Especially, PfAP2-EXP2 regulation of euchromatic genes does not depend on recognizing specific DNA sequences, while a CCCTAAACCC motif is found in its heterochromatic binding sites. Combined with transcriptome profiling, we suggest that PfAP2-EXP2 is participated in the intraerythrocytic development by affecting the expression of genes related to cell remodeling at the schizont stage. In summary, this study explores an ApiAP2 member plays an important role for the P. falciparum blood-stage replication, which suggests a new perspective for malaria elimination.

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.

OUP accepted manuscript

Heterochromatin-associated gene silencing controls multiple physiological processes in malaria parasites, however, little is known concerning the regulatory network and cis-acting sequences involved in the organization of heterochromatin and how they modulate heterochromatic gene expression. Based on systematic profiling of genome-wide occupancy of eighteen Apicomplexan AP2 transcription factors by ChIP-seq analysis, we identify and characterize eight heterochromatin-associated factors (PfAP2-HFs), which exhibit preferential enrichment within heterochromatic regions but with differential coverage profiles. Although these ApiAP2s target euchromatic gene loci via specific DNA motifs, they are likely integral components of heterochromatin independent of DNA motif recognition. Systematic knockout screenings of ApiAP2 factors coupled with RNA-seq transcriptomic profiling revealed three activators and three repressors of heterochromatic gene expression including four PfAP2-HFs. Notably, expression of virulence genes is either completely silenced or significantly reduced upon the depletion of PfAP2-HC. Integrated multi-omics analyses reveal autoregulation and feed-forward loops to be common features of the ApiAP2 regulatory network, in addition to the occurrence of dynamic interplay between local chromatin structure and ApiAP2s in transcriptional control. Collectively, this study provides a valuable resource describing the genome-wide landscape of the ApiAP2 family and insights into functional divergence and cooperation within this family during the blood-stage development of malaria parasites.

Decrypting the complexity of the human malaria parasite biology through systems biology approaches

The human malaria parasite, Plasmodium falciparum, is a unicellular protozoan responsible for over half a million deaths annually. With a complex life cycle alternating between human and invertebrate hosts, this apicomplexan is notoriously adept at evading host immune responses and developing resistance to all clinically administered treatments. Advances in omics-based technologies, increased sensitivity of sequencing platforms and enhanced CRISPR based gene editing tools, have given researchers access to more in-depth and untapped information about this enigmatic micro-organism, a feat thought to be infeasible in the past decade. Here we discuss some of the most important scientific achievements made over the past few years with a focus on novel technologies and platforms that set the stage for subsequent discoveries. We also describe some of the systems-based methods applied to uncover gaps of knowledge left through single-omics applications with the hope that we will soon be able to overcome the spread of this life-threatening disease.

Essential role of a Plasmodium berghei heat shock protein (PBANKA_0938300) in gametocyte development

The continued existence of Plasmodium parasites in physiologically distinct environments during their transmission in mosquitoes and vertebrate hosts requires effector proteins encoded by parasite genes to provide adaptability. Parasites utilize their robust stress response system involving heat shock proteins for their survival. Molecular chaperones are involved in maintaining protein homeostasis within a cell during stress, protein biogenesis and the formation of protein complexes. Due to their critical role in parasite virulence, they are considered targets for therapeutic interventions. Our results identified a putative P. berghei heat shock protein (HSP) belonging to the HSP40 family (HspJ62), which is abundantly induced upon heat stress and expressed during all parasite stages. To determine the role HspJ62, a gene-disrupted P. berghei transgenic line was developed (ΔHspJ62), which resulted in disruption of gametocyte formation. Such parasites were unable to form subsequent sexual stages because of disrupted gametogenesis, indicating the essential role of HspJ62 in gametocyte formation. Transcriptomic analysis of the transgenic line showed downregulation of a number of genes, most of which were specific to male or female gametocytes. The transcription factor ApiAP2 was also downregulated in ΔHspJ62 parasites. Our findings suggest that the downregulation of ApiAP2 likely disrupts the transcriptional regulation of sexual stage genes, leading to impaired gametogenesis. This finding also highlights the critical role that HspJ62 indirectly plays in the development of P. berghei sexual stages and in facilitating the conversion from the asexual blood stage to the sexual stage. This study characterizes the HspJ62 protein as a fertility factor because parasites lacking it are unable to transmit to mosquitoes. This study adds an important contribution to ongoing research aimed at understanding gametocyte differentiation and formation in parasites. The molecule adds to the list of potential drug targets that can be targeted to inhibit parasite sexual development and consequently parasite transmission.

Full-length transcriptome analysis and identification of transcript structures in Eimeria necatrix from different developmental stages by single-molecule real-time sequencing

Background

Eimeria necatrix is one of the most pathogenic parasites, causing high mortality in chickens. Although its genome sequence has been published, the sequences and complete structures of its mRNA transcripts remain unclear, limiting exploration of novel biomarkers, drug targets and genetic functions in E. necatrix.

Methods

Second-generation merozoites (MZ-2) of E. necatrix were collected using Percoll density gradients, and high-quality RNA was extracted from them. Single-molecule real-time (SMRT) sequencing and Illumina sequencing were combined to generate the transcripts of MZ-2. Combined with the SMRT sequencing data of sporozoites (SZ) collected in our previous study, the transcriptome and transcript structures of E. necatrix were studied.

Results

SMRT sequencing yielded 21,923 consensus isoforms in MZ-2. A total of 17,151 novel isoforms of known genes and 3918 isoforms of novel genes were successfully identified. We also identified 2752 (SZ) and 3255 (MZ-2) alternative splicing (AS) events, 1705 (SZ) and 1874 (MZ-2) genes with alternative polyadenylation (APA) sites, 4019 (SZ) and 2588 (MZ-2) fusion transcripts, 159 (SZ) and 84 (MZ-2) putative transcription factors (TFs) and 3581 (SZ) and 2039 (MZ-2) long non-coding RNAs (lncRNAs). To validate fusion transcripts, reverse transcription-PCR was performed on 16 candidates, with an accuracy reaching up to 87.5%. Sanger sequencing of the PCR products further confirmed the authenticity of chimeric transcripts. Comparative analysis of transcript structures revealed a total of 3710 consensus isoforms, 815 AS events, 1139 genes with APA sites, 20 putative TFs and 352 lncRNAs in both SZ and MZ-2.

Conclusions

We obtained many long-read isoforms in E. necatrix SZ and MZ-2, from which a series of lncRNAs, AS events, APA events and fusion transcripts were identified. Information on TFs will improve understanding of transcriptional regulation, and fusion event data will greatly improve draft versions of gene models in E. necatrix. This information offers insights into the mechanisms governing the development of E. necatrix and will aid in the development of novel strategies for coccidiosis control.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05015-7.

A heat-shock response regulated by the PfAP2-HS transcription factor protects human malaria parasites from febrile temperatures

Periodic fever is a characteristic clinical feature of human malaria, but how parasites survive febrile episodes is not known. Although the genomes of Plasmodium species encode a full set of chaperones, they lack the conserved eukaryotic transcription factor HSF1, which activates the expression of chaperones following heat shock. Here, we show that PfAP2-HS, a transcription factor in the ApiAP2 family, regulates the protective heat-shock response in Plasmodium falciparum. PfAP2-HS activates the transcription of hsp70-1 and hsp90 at elevated temperatures. The main binding site of PfAP2-HS in the entire genome coincides with a tandem G-box DNA motif in the hsp70-1 promoter. Engineered parasites lacking PfAP2-HS have reduced heat-shock survival and severe growth defects at 37 °C but not at 35 °C. Parasites lacking PfAP2-HS also have increased sensitivity to imbalances in protein homeostasis (proteostasis) produced by artemisinin, the frontline antimalarial drug, or the proteasome inhibitor epoxomicin. We propose that PfAP2-HS contributes to the maintenance of proteostasis under basal conditions and upregulates specific chaperone-encoding genes at febrile temperatures to protect the parasite against protein damage.

A unique GCN5 histone acetyltransferase complex controls erythrocyte invasion and virulence in the malaria parasite Plasmodium falciparum

The histone acetyltransferase GCN5-associated SAGA complex is evolutionarily conserved from yeast to human and functions as a general transcription co-activator in global gene regulation. In this study, we identified a divergent GCN5 complex in Plasmodium falciparum, which contains two plant homeodomain (PHD) proteins (PfPHD1 and PfPHD2) and a plant apetela2 (AP2)-domain transcription factor (PfAP2-LT). To dissect the functions of the PfGCN5 complex, we generated parasite lines with either the bromodomain in PfGCN5 or the PHD domain in PfPHD1 deleted. The two deletion mutants closely phenocopied each other, exhibiting significantly reduced merozoite invasion of erythrocytes and elevated sexual conversion. These domain deletions caused dramatic decreases not only in histone H3K9 acetylation but also in H3K4 trimethylation, indicating synergistic crosstalk between the two euchromatin marks. Domain deletion in either PfGCN5 or PfPHD1 profoundly disturbed the global transcription pattern, causing altered expression of more than 60% of the genes. At the schizont stage, these domain deletions were linked to specific down-regulation of merozoite genes involved in erythrocyte invasion, many of which contain the AP2-LT binding motif and are also regulated by AP2-I and BDP1, suggesting targeted recruitment of the PfGCN5 complex to the invasion genes by these specific factors. Conversely, at the ring stage, PfGCN5 or PfPHD1 domain deletions disrupted the mutually exclusive expression pattern of the entire var gene family, which encodes the virulent factor PfEMP1. Correlation analysis between the chromatin state and alteration of gene expression demonstrated that up- and down-regulated genes in these mutants are highly correlated with the silent and active chromatin states in the wild-type parasite, respectively. Collectively, the PfGCN5 complex represents a novel HAT complex with a unique subunit composition including an AP2 transcription factor, which signifies a new paradigm for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

Epigenetic regulation of gene expression plays essential roles in orchestrating the general and parasite-specific cellular pathways in the malaria parasite Plasmodium falciparum. To better understand the epigenetic mechanisms in this parasite, we characterized the histone acetyltransferase GCN5-mediated transcription regulation during intraerythrocytic development of the parasite. Using tandem affinity purification and proteomic characterization, we identified that the PfGCN5-associated complex contains nine core components, including two PHD domain proteins (PfPHD1 and PfPHD2) and an AP2-domain transcription factor, which is divergent from the canonical GCN5 complexes evolutionarily conserved from yeast to human. To understand the functions of the PfGCN5 complex, we performed domain deletions in two subunits of this complex, PfGCN5 and PfPHD1. We found that the two deletion mutants displayed very similar growth phenotypes, including significantly reduced merozoite invasion rates and elevated sexual conversion. These two mutants were associated with dramatic decreases in histone H3K9 acetylation and H3K4 trimethylation, which led to global changes in chromatin states and gene expression. Consistent with the phenotypes, genes significantly affected by the PfGCN5 and PfPHD1 gene disruption include those participating in parasite-specific pathways such as invasion, virulence, and sexual development. In conclusion, this study presents a new model of the PfGCN5 complex for targeting the co-activator complex to regulate general and parasite-specific cellular processes in this low-branching parasitic protist.

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.

The Transcription Factor PfAP2-O Influences Virulence Gene Transcription and Sexual Development in Plasmodium falciparum

The human malaria parasite Plasmodium falciparum expresses variant PfEMP1 proteins on the infected erythrocyte, which function as ligands for endothelial receptors in capillary vessels, leading to erythrocyte sequestration and severe malaria. The factors that orchestrate the mono-allelic expression of the 45–90 PfEMP1-encoding var genes within each parasite genome are still not fully identified. Here, we show that the transcription factor PfAP2-O influences the transcription of var genes. The temporary knockdown of PfAP2-O leads to a complete loss of var transcriptional memory and a decrease in cytoadherence in CD36 adherent parasites. AP2-O-knocked-down parasites exhibited also significant reductions in transmission through Anopheles mosquitoes. We propose that PfAP2-O is, beside its role in transmission stages, also one of the virulence gene transcriptional regulators and may therefore be exploited as an important target to disrupt severe malaria and block parasite transmission.

Mechanisms of triggering malaria gametocytogenesis by AP2-G

The transcription factor (TF) AP2-G is essential for gametocytogenesis in the malaria parasite; however, it remains unclear if AP2-G determines commitment to sexual stage development fate in the schizont stage, or whether AP2-G directly initiates sexual stage differentiation and development beginning in the late-trophozoite stage. In this study, we addressed this issue by investigating the expression profile of AP2-G and determining genome-wide target genes in Plasmodium berghei. Fluorescence microscopy showed that AP2-G expression was first observed in the parasite 12 h after erythrocyte invasion and peaked at 18 h when sexual features were first manifested in early gametocytes. Expression of AP2-G decreased with manifestation of sex-specific features. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) was performed at peak AP2-G expression and identified over 1000 binding sites in the genome. The main binding motif of the TF predicted from the binding sites was GTACNY. Predicted targets contained a number of genes related to protein biogenesis, suggesting that AP2-G plays a role in establishing a cellular basis required for sexual differentiation. AP2-G binding sites also existed upstream of gametocyte-specific TFs, namely AP2-G2, AP2-FG, and AP2-G itself. Furthermore, the target contained two AP2 TF-related genes. Disruption of these genes resulted in the arrest of ookinete development. These results suggest another role of AP2-G: activating a transcriptional cascade to promote conversion into early gametocytes. Taken together, AP2-G is involved not in establishing sexual commitment of schizonts, but rather in triggering the initiation of differentiation and the early development of gametocytes in the late trophozoite stage.

Global transcriptome landscape of the rabbit protozoan parasite Eimeria stiedae

Background

Coccidiosis caused by Eimeria stiedae is a widespread and economically significant disease of rabbits. The lack of studies on the life-cycle development and host interactions of E. stiedae at the molecular level has hampered our understanding of its pathogenesis.

Methods

In this study, we present a comprehensive transcriptome landscape of E. stiedae to illustrate its dynamic development from unsporulated oocysts to sporulated oocysts, merozoites, and gametocytes, and to identify genes related to parasite-host interactions during parasitism using combined PacBio single-molecule real-time and Illumina RNA sequencing followed by bioinformatics analysis and qRT-PCR validation.

Results

In total, 12,582 non-redundant full-length transcripts were generated with an average length of 1808 bp from the life-cycle stages of E. stiedae. Pairwise comparisons between stages revealed 8775 differentially expressed genes (DEGs) showing highly significant description changes, which compiled a snapshot of the mechanisms underlining asexual and sexual biology of E. stiedae including oocyst sporulation between unsporulated and sporulated oocysts; merozoite replication between sporulated oocysts and merozoites; and gametophyte development and gamete generation between merozoites and gametocytes. Further, 248 DEGs were grouped into nine series clusters and five groups by expression patterns, and showed that parasite–host interaction-related genes predominated in merozoites and gametocytes and were mostly involved in steroid biosynthesis and lipid metabolism and carboxylic acid. Additionally, co-expression analyses identified genes associated with development and host invasion in unsporulated and sporulated oocysts and immune interactions during gametocyte parasitism.

Conclusions

This is the first study, to our knowledge, to use the global transcriptome profiles to decipher molecular changes across the E. stiedae life cycle, and these results not only provide important information for the molecular characterization of E. stiedae, but also offer valuable resources to study other apicomplexan parasites with veterinary and public significance.

Graphic Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04811-5.

A single-cell atlas of Plasmodium falciparum transmission through the mosquito

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.

Here the authors use single-cell RNA-seq to profile the transmission stages of the human malaria parasite Plasmodium falciparum as it progresses through the Anopheles mosquito. They highlight unique patterns of gene usage throughout this development and identify potential pleiotropic genes that function at multiple life cycle stages.

Peculiarities of Plasmodium falciparum Gene Regulation and Chromatin Structure

The highly complex life cycle of the human malaria parasite, Plasmodium falciparum, is based on an orchestrated and tightly regulated gene expression program. In general, eukaryotic transcription regulation is determined by a combination of sequence-specific transcription factors binding to regulatory DNA elements and the packaging of DNA into chromatin as an additional layer. The accessibility of regulatory DNA elements is controlled by the nucleosome occupancy and changes of their positions by an active process called nucleosome remodeling. These epigenetic mechanisms are poorly explored in P. falciparum. The parasite genome is characterized by an extraordinarily high AT-content and the distinct architecture of functional elements, and chromatin-related proteins also exhibit high sequence divergence compared to other eukaryotes. Together with the distinct biochemical properties of nucleosomes, these features suggest substantial differences in chromatin-dependent regulation. Here, we highlight the peculiarities of epigenetic mechanisms in P. falciparum, addressing chromatin structure and dynamics with respect to their impact on transcriptional control. We focus on the specialized chromatin remodeling enzymes and discuss their essential function in P. falciparum gene regulation.

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.

The ApiAP2 factor PfAP2-HC is an integral component of heterochromatin in the malaria parasite Plasmodium falciparum

Malaria parasites undergo a complex life cycle in the human host and the mosquito vector. The ApiAP2 family of DNA-binding proteins plays a dominant role in parasite development and life cycle progression. Most ApiAP2 factors studied to date act as transcription factors regulating stage-specific gene expression. Here, we characterized an ApiAP2 factor in Plasmodium falciparum that we termed PfAP2-HC. We demonstrate that PfAP2-HC specifically binds to heterochromatin throughout the genome. Intriguingly, PfAP2-HC does not bind DNA in vivo and recruitment of PfAP2-HC to heterochromatin is independent of its DNA-binding domain but strictly dependent on heterochromatin protein 1. Furthermore, our results suggest that PfAP2-HC functions neither in the regulation of gene expression nor in heterochromatin formation or maintenance. In summary, our findings reveal PfAP2-HC as a core component of heterochromatin in malaria parasites and identify unexpected properties and substantial functional divergence among the members of the ApiAP2 family of regulatory proteins.

•

The ApiAP2 factor AP2-HC is a core component of heterochromatin in malaria parasites

•

Binding of AP2-HC to heterochromatin strictly depends on heterochromatin protein 1

•

The AP2 DNA-binding domain of AP2-HC is dispensable for heterochromatin association

Plasmodium Reproduction, Cell Size, and Transcription: How to Cope With Increasing DNA Content?

Plasmodium, the unicellular parasite that causes malaria, evolved a highly unusual mode of reproduction. During its complex life cycle, invasive or transmissive stages alternate with proliferating stages, where a single parasite can produce tens of thousands of progeny. In the clinically relevant blood stage of infection, the parasite replicates its genome up to thirty times and forms a multinucleated cell before daughter cells are assembled. Thus, within a single cell cycle, Plasmodium develops from a haploid to a polypoid cell, harboring multiple copies of its genome. Polyploidy creates several biological challenges, such as imbalances in genome output, and cells can respond to this by changing their size and/or alter the production of RNA species and protein to achieve expression homeostasis. However, the effects and possible adaptations of Plasmodium to the massively increasing DNA content are unknown. Here, we revisit and embed current Plasmodium literature in the context of polyploidy and propose potential mechanisms of the parasite to cope with the increasing gene dosage.

Identification of long regulatory elements in the genome of Plasmodium falciparum and other eukaryotes

Long regulatory elements (LREs), such as CpG islands, polydA:dT tracts or AU-rich elements, are thought to play key roles in gene regulation but, as opposed to conventional binding sites of transcription factors, few methods have been proposed to formally and automatically characterize them. We present here a computational approach named DExTER (Domain Exploration To Explain gene Regulation) dedicated to the identification of candidate LREs (cLREs) and apply it to the analysis of the genomes of P. falciparum and other eukaryotes. Our analyses show that all tested genomes contain several cLREs that are somewhat conserved along evolution, and that gene expression can be predicted with surprising accuracy on the basis of these long regions only. Regulation by cLREs exhibits very different behaviours depending on species and conditions. In P. falciparum and other Apicomplexan organisms as well as in Dictyostelium discoideum, the process appears highly dynamic, with different cLREs involved at different phases of the life cycle. For multicellular organisms, the same cLREs are involved in all tissues, but a dynamic behavior is observed along embryonic development stages. In P. falciparum, whose genome is known to be strongly depleted of transcription factors, cLREs are predictive of expression with an accuracy above 70%, and our analyses show that they are associated with both transcriptional and post-transcriptional regulation signals. Moreover, we assessed the biological relevance of one LRE discovered by DExTER in P. falciparum using an in vivo reporter assay. The source code (python) of DExTER is available at https://gite.lirmm.fr/menichelli/DExTER.

Plasmodium transcription repressor AP2-O3 regulates sex-specific identity of gene expression in female gametocytes

Male and female gametocytes are sexual precursor cells essential for mosquito transmission of malaria parasite. Differentiation of gametocytes into fertile gametes (known as gametogenesis) relies on the gender‐specific transcription program. How the parasites establish distinct repertoires of transcription in the male and female gametocytes remains largely unknown. Here, we report that an Apetala2 family transcription factor AP2‐O3 operates as a transcription repressor in the female gametocytes. AP2‐O3 is specifically expressed in the female gametocytes. AP2‐O3‐deficient parasites produce apparently normal female gametocytes. Nevertheless, these gametocytes fail to differentiate into fully fertile female gametes, leading to developmental arrest in fertilization and early development post‐fertilization. AP2‐O3 disruption causes massive upregulation of transcriptionally dormant male genes and simultaneously downregulation of highly transcribed female genes in the female gametocytes. AP2‐O3 targets a substantial proportion of the male genes by recognizing an 8‐base DNA motif. In addition, the maternal AP2‐O3 is removed after fertilization, which is required for the zygote to ookinete development. Therefore, the global transcriptional repression of the male genes in the female gametocytes is required for safeguarding female‐specific transcriptome and essential for the mosquito transmission of Plasmodium.

PfAP2-G2 Is Associated to Production and Maturation of Gametocytes in Plasmodium falciparum via Regulating the Expression of PfMDV-1

Gametocyte is the sole form of the Plasmodium falciparum which is transmissible to the mosquito vector. Here, we report that an Apicomplexan Apetala2 (ApiAP2) family transcription factor, PfAP2-G2 (Pf3D7_1408200), plays a role in the development of gametocytes in P. falciparum by regulating the expression of PfMDV-1 (Pf3D7_1216500). Reverse transcriptase-quantitative PCR (RT-qPCR) analysis showed that PfAP2-G2 was highly expressed in the ring stage. Indirect immunofluorescence assay showed nuclear localization of PfAP2-G2 in asexual stages. The knockout of PfAP2-G2 led to a ~95% decrease in the number of mature gametocytes with a more substantial influence on the production and maturation of the male gametocytes, resulting in a higher female/male gametocyte ratio. To test the mechanism of this phenotype, RNA-seq and RT-qPCR showed that disruption of PfAP2-G2 led to the down-regulation of male development gene-1 (PfMDV-1) in asexual stages. We further found that PfAP2-G2 was enriched at the transcriptional start site (TSS) of PfMDV-1 by chromatin immunoprecipitation and qPCR assay in both ring stage and schizont stage, which demonstrated that PfMDV-1 is one of the targets of PfAP2-G2. In addition, RT-qPCR also showed that PfAP2-G (Pf3D7_1222600), the master regulator for sexual commitment, was also down-regulated in the PfAP2-G2 knockout parasites in the schizont stage, but no change in the ring stage. This phenomenon suggested that PfAP2-G2 played a role at the asexual stage for the development of parasite gametocytes and warrants further investigations in regulatory pathways of PfAP2-G2.

Preparing for Transmission: Gene Regulation in Plasmodium Sporozoites

Plasmodium sporozoites are transmitted to mammals by anopheline mosquitoes and first infect the liver, where they transform into replicative exoerythrocytic forms, which subsequently release thousands of merozoites that invade erythrocytes and initiate the malaria disease. In some species, sporozoites can transform into dormant hypnozoites in the liver, which cause malaria relapses upon reactivation. Transmission from the insect vector to a mammalian host is a critical step of the parasite life cycle, and requires tightly regulated gene expression. Sporozoites are formed inside oocysts in the mosquito midgut and become fully infectious after colonization of the insect salivary glands, where they remain quiescent until transmission. Parasite maturation into infectious sporozoites is associated with reprogramming of the sporozoite transcriptome and proteome, which depends on multiple layers of transcriptional and post-transcriptional regulatory mechanisms. An emerging scheme is that gene expression in Plasmodium sporozoites is controlled by alternating waves of transcription activity and translational repression, which shape the parasite RNA and protein repertoires for successful transition from the mosquito vector to the mammalian host.

A systems-level gene regulatory network model for Plasmodium falciparum

Many of the gene regulatory processes of Plasmodium falciparum, the deadliest malaria parasite, remain poorly understood. To develop a comprehensive guide for exploring this organism's gene regulatory network, we generated a systems-level model of P. falciparum gene regulation using a well-validated, machine-learning approach for predicting interactions between transcription regulators and their targets. The resulting network accurately predicts expression levels of transcriptionally coherent gene regulatory programs in independent transcriptomic data sets from parasites collected by different research groups in diverse laboratory and field settings. Thus, our results indicate that our gene regulatory model has predictive power and utility as a hypothesis-generating tool for illuminating clinically relevant gene regulatory mechanisms within P. falciparum. Using the set of regulatory programs we identified, we also investigated correlates of artemisinin resistance based on gene expression coherence. We report that resistance is associated with incoherent expression across many regulatory programs, including those controlling genes associated with erythrocyte-host engagement. These results suggest that parasite populations with reduced artemisinin sensitivity are more transcriptionally heterogenous. This pattern is consistent with a model where the parasite utilizes bet-hedging strategies to diversify the population, rendering a subpopulation more able to navigate drug treatment.