What functional genomics has taught us about transcriptional regulation in malaria parasites

Ectopic overexpression of Plasmodium falciparum DNA-/RNA-binding Alba proteins misregulates virulence gene homeostasis during asexual blood development

Alba domain-containing proteins are ubiquitously found in archaea and eukaryotes. By binding to either DNA, RNA, or DNA:RNA hybrids, these proteins function in genome stabilization, chromatin organization, gene regulation, and/or translational modulation. In the malaria parasite Plasmodium falciparum, six Alba domain proteins PfAlba1-6 have been described, of which PfAlba1 has emerged as a "master regulator" of translation during parasite intra-erythrocytic development (IED). Given that a tight control of gene expression is especially important during IED, when malaria pathogenesis manifests, in this study, we focus on three other P. falciparum Albas, PfAlba2-4. Because genetic manipulation of the genomic loci of PfAlba2-4 was unsuccessful, we overexpressed each of these proteins from an episome under a strong constitutive promoter. We observed that PfAlba2 or PfAlba3 overexpression strongly reduced parasite growth and impacted IED stage transitions. In contrast, elevated levels of PfAlba4 were well-tolerated by the parasite. In keeping with this, differential gene expression analysis using RNA-seq of PfAlba2 or PfAlba3 overexpressing strains revealed a significant misregulation of mRNAs encoding virulence factors, such as those related to erythrocyte invasion; a general repression of var gene expression was also apparent. PfAlba4 overexpression, on the other hand, did not significantly perturb the steady-state transcriptome of IED stages and appeared to enhance var mRNA levels. Moreover, distinct sets of genes were targeted by each PfAlba for regulation. Taken together, this study highlights the nonredundant roles of PfAlba proteins in the P. falciparum IED, emphasizing their importance in subtelomeric chromatin biology and RNA regulation.IMPORTANCEThe malaria parasite Plasmodium falciparum tightly controls the expression of its genes at the epigenetic, transcriptional, post-transcriptional, and translational levels to synthesize essential proteins, including virulence factors, in a timely and spatially coordinated manner. A family of six proteins implicated in this process is called PfAlba, characterized by the presence of the DNA-, RNA- or DNA:RNA hybrid-binding Alba domain. To better understand the cellular pathways regulated by this protein family, we overexpressed three PfAlbas during P. falciparum intra-erythrocytic growth and found that high levels of PfAlba2 and PfAlba3 were detrimental to parasite development. This was accompanied by significant changes in the parasite's transcriptome, either with regards to mRNA steady-state levels or expression timing. PfAlba4 overexpression, on the other hand, was well-tolerated by the parasite. Overall, our results delineate specific pathways targeted by individual PfAlbas for regulation and link PfAlba2/PfAlba3 to mutually exclusive expression of the virulence-promoting surface antigen PfEMP1.

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.

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.

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.

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.

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.

A transcription factor helps Plasmodium falciparum gametocytogenesis take shape

The recent study by Campelo Morillo et al. has shown that one of the small number of non-ApiAP2 DNA-binding proteins in the Plasmodium falciparum genome acts as a transcription factor in the gametocytogenesis cascade and is responsible for the gametocyte's distinctive morphology.

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.

Emerging biology of noncoding RNAs in malaria parasites

In eukaryotic organisms, noncoding RNAs (ncRNAs) have been implicated as important regulators of multifaceted biological processes, including transcriptional, posttranscriptional, and epigenetic regulation of gene expression. In recent years, it is becoming clear that protozoan parasites encode diverse ncRNA transcripts; however, little is known about their cellular functions. Recent advances in high-throughput “omic” studies identified many novel long ncRNAs (lncRNAs) in apicomplexan parasites, some of which undergo splicing, polyadenylation, and encode small proteins. To date, only a few of them are characterized, leaving a big gap in our understanding regarding their origin, mode of action, and functions in parasite biology. In this review, we focus on lncRNAs of the human malaria parasite Plasmodium falciparum and highlight their cellular functions and possible mechanisms of action.

Plasmodium vivax pre-erythrocytic stages and the latent hypnozoite

Plasmodium vivax is the most geographically widespread malaria parasite on the planet. This is largely because after mosquito transmission, P. vivax sporozoites can invade hepatocytes and form latent liver stages known as hypnozoites. These persistent liver stages can activate weeks, months or even years after an infected individual suffers a primary clinical infection. Activation then leads to replication and liver stage schizont maturation that ultimately cause relapse of blood stage infection, disease, and onward transmission. Thus, the latent hypnozoite can lie in wait during times when onward transmission is unlikely due to conditions that do not favor the mosquito. For example, in temperate climates where mosquito prevalence is only seasonal. Furthermore, the elimination of hypnozoites is challenging since the hypnozoite reservoir is currently undetectable and not killed by most antimalarial drugs. Here, we review our current knowledge of the pre-erythrocytic stages of the malaria parasite - the sporozoite and liver stages, including the elusive and enigmatic hypnozoite. We focus on our understanding of sporozoite biology, the novel animal models that are available to study the hypnozoite and hypnozoite activation and the ongoing efforts to understand the biological makeup of the hypnozoite that allow for its persistence in the human host.

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 Modular Circuitry of Apicomplexan Cell Division Plasticity

The close-knit group of apicomplexan parasites displays a wide variety of cell division modes, which differ between parasites as well as between different life stages within a single parasite species. The beginning and endpoint of the asexual replication cycles is a 'zoite' harboring the defining apical organelles required for host cell invasion. However, the number of zoites produced per division round varies dramatically and can unfold in several different ways. This plasticity of the cell division cycle originates from a combination of hard-wired developmental programs modulated by environmental triggers. Although the environmental triggers and sensors differ between species and developmental stages, widely conserved secondary messengers mediate the signal transduction pathways. These environmental and genetic input integrate in division-mode specific chromosome organization and chromatin modifications that set the stage for each division mode. Cell cycle progression is conveyed by a smorgasbord of positively and negatively acting transcription factors, often acting in concert with epigenetic reader complexes, that can vary dramatically between species as well as division modes. A unique set of cell cycle regulators with spatially distinct localization patterns insert discrete check points which permit individual control and can uncouple general cell cycle progression from nuclear amplification. Clusters of expressed genes are grouped into four functional modules seen in all division modes: 1. mother cytoskeleton disassembly; 2. DNA replication and segregation (D&S); 3. karyokinesis; 4. zoite assembly. A plug-and-play strategy results in the variety of extant division modes. The timing of mother cytoskeleton disassembly is hard-wired at the species level for asexual division modes: it is either the first step, or it is the last step. In the former scenario zoite assembly occurs at the plasma membrane (external budding), and in the latter scenario zoites are assembled in the cytoplasm (internal budding). The number of times each other module is repeated can vary regardless of this first decision, and defines the modes of cell division: schizogony, binary fission, endodyogeny, endopolygeny.

Direct Nanopore Sequencing of mRNA Reveals Landscape of Transcript Isoforms in Apicomplexan Parasites

Alternative splicing is a widespread phenomenon in metazoans by which single genes are able to produce multiple isoforms of the gene product. However, this has been poorly characterized in apicomplexans, a major phylum of some of the most important global parasites. Efforts have been hampered by atypical transcriptomic features, such as the high AU content of Plasmodium RNA, but also the limitations of short-read sequencing in deciphering complex splicing events. In this study, we utilized the long read direct RNA sequencing platform developed by Oxford Nanopore Technologies to survey the alternative splicing landscape of Toxoplasma gondii and Plasmodium falciparum. We find that while native RNA sequencing has a reduced throughput, it allows us to obtain full-length or nearly full-length transcripts with comparable quantification to Illumina sequencing. By comparing these data with available gene models, we find widespread alternative splicing, particularly intron retention, in these parasites. Most of these transcripts contain premature stop codons, suggesting that in these parasites, alternative splicing represents a pathway to transcriptomic diversity, rather than expanding proteomic diversity. Moreover, alternative splicing rates are comparable between parasites, suggesting a shared splicing machinery, despite notable transcriptomic differences between the parasites. This study highlights a strategy in using long-read sequencing to understand splicing events at the whole-transcript level and has implications in the future interpretation of transcriptome sequencing studies.

IMPORTANCE We have used a novel nanopore sequencing technology to directly analyze parasite transcriptomes. The very long reads of this technology reveal the full-length genes of the parasites that cause malaria and toxoplasmosis. Gene transcripts must be processed in a process called splicing before they can be translated to protein. Our analysis reveals that these parasites very frequently only partially process their gene products, in a manner that departs dramatically from their human hosts.

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.

The highly diverse TATA box-binding proteins among protists: A review

Transcription is the first step of gene expression regulation and is a fundamental mechanism for establishing the viability and development of a cell. The TATA box-binding protein (TBP) interaction with a TATA box in a promoter is one of the best studied mechanisms in transcription initiation. TBP is a transcription factor that is highly conserved from archaea to humans and is essential for the transcription initiated by each of the three RNA polymerases. In addition, the discovery of TBP-related factor 1 (TRF1) and other factors related to TBP shed light on the variability among transcription initiation complexes, thus demonstrating that the compositions of these complexes are, in fact, more complicated than originally believed. Despite these facts, the majority of studies on transcription have been performed on animal, plant and fungal cells, which serve as canonical models, and information regarding protist cells is relatively scarce. The aim of this work is to review the diversity of the TBPs that have been documented in protists and describe some of the specific features that differentiate them from their counterparts in higher eukaryotes.