The Apicomplexan AP2 family: integral factors regulating Plasmodium development

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.

Size-dependent enhancement of gene expression by Plasmodium 5'UTR introns

BACKGROUND

Eukaryotic genes contain introns that are removed by the spliceosomal machinery during mRNA maturation. Introns impose a huge energetic burden on a cell; therefore, they must play an essential role in maintaining genome stability and/or regulating gene expression. Many genes (> 50%) in Plasmodium parasites contain predicted introns, including introns in 5' and 3' untranslated regions (UTR). However, the roles of UTR introns in the gene expression of malaria parasites remain unknown.

METHODS

In this study, an episomal dual-luciferase assay was developed to evaluate gene expression driven by promoters with or without a 5'UTR intron from four Plasmodium yoelii genes. To investigate the effect of the 5'UTR intron on endogenous gene expression, the pytctp gene was tagged with 3xHA at the N-terminal of the coding region, and parasites with or without the 5'UTR intron were generated using the CRISPR/Cas9 system.

RESULTS

We showed that promoters with 5'UTR introns had higher activities in driving gene expression than those without 5'UTR introns. The results were confirmed in recombinant parasites expressing an HA-tagged gene (pytctp) driven by promoter with or without 5'UTR intron. The enhancement of gene expression was intron size dependent, but not the DNA sequence, e.g. the longer the intron, the higher levels of expression. Similar results were observed when a promoter from one strain of P. yoelii was introduced into different parasite strains. Finally, the 5'UTR introns were alternatively spliced in different parasite development stages, suggesting an active mechanism employed by the parasites to regulate gene expression in various developmental stages.

CONCLUSIONS

Plasmodium 5'UTR introns enhance gene expression in a size-dependent manner; the presence of alternatively spliced mRNAs in different parasite developmental stages suggests that alternative slicing of 5'UTR introns is one of the key mechanisms in regulating parasite gene expression and differentiation.

Telomere maintenance in African trypanosomes

Telomere maintenance is essential for genome integrity and chromosome stability in eukaryotic cells harboring linear chromosomes, as telomere forms a specialized structure to mask the natural chromosome ends from DNA damage repair machineries and to prevent nucleolytic degradation of the telomeric DNA. In Trypanosoma brucei and several other microbial pathogens, virulence genes involved in antigenic variation, a key pathogenesis mechanism essential for host immune evasion and long-term infections, are located at subtelomeres, and expression and switching of these major surface antigens are regulated by telomere proteins and the telomere structure. Therefore, understanding telomere maintenance mechanisms and how these pathogens achieve a balance between stability and plasticity at telomere/subtelomere will help develop better means to eradicate human diseases caused by these pathogens. Telomere replication faces several challenges, and the "end replication problem" is a key obstacle that can cause progressive telomere shortening in proliferating cells. To overcome this challenge, most eukaryotes use telomerase to extend the G-rich telomere strand. In addition, a number of telomere proteins use sophisticated mechanisms to coordinate the telomerase-mediated de novo telomere G-strand synthesis and the telomere C-strand fill-in, which has been extensively studied in mammalian cells. However, we recently discovered that trypanosomes lack many telomere proteins identified in its mammalian host that are critical for telomere end processing. Rather, T. brucei uses a unique DNA polymerase, PolIE that belongs to the DNA polymerase A family (E. coli DNA PolI family), to coordinate the telomere G- and C-strand syntheses. In this review, I will first briefly summarize current understanding of telomere end processing in mammals. Subsequently, I will describe PolIE-mediated coordination of telomere G- and C-strand synthesis in T. brucei and implication of this recent discovery.

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.

Unravelling the sexual developmental biology of Cystoisospora suis, a model for comparative coccidian parasite studies

Introduction

The apicomplexan parasite Cystoisospora suis has global significance as an enteropathogen of suckling piglets. Its intricate life cycle entails a transition from an asexual phase to sexual development, ultimately leading to the formation of transmissible oocysts.

Methods

To advance our understanding of the parasite's cellular development, we complemented previous transcriptome studies by delving into the proteome profiles at five distinct time points of in vitro cultivation through LC/MS-MS analysis.

Results

A total of 1,324 proteins were identified in the in vitro developmental stages of C. suis, and 1,082 proteins were identified as significantly differentially expressed. Data are available via ProteomeXchange with identifier PXD045050. We performed BLAST, GO enrichment, and KEGG pathway analyses on the up- and downregulated proteins to elucidate correlated events in the C. suis life cycle. Our analyses revealed intriguing metabolic patterns in macromolecule metabolism, DNA- and RNA-related processes, proteins associated with sexual stages, and those involved in cell invasion, reflecting the adaptation of sexual stages to a nutrient-poor and potentially stressful extracellular environment, with a focus on enzymes involved in metabolism and energy production.

Discussion

These findings have important implications for understanding the developmental biology of C. suis as well as other, related coccidian parasites, such as Eimeria spp. and Toxoplasma gondii. They also support the role of C. suis as a new model for the comparative biology of coccidian tissue cyst 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.

PbAP2-FG2 and PbAP2R-2 function together as a transcriptional repressor complex essential for Plasmodium female development

Gametocyte development is a critical step in the life cycle of Plasmodium. Despite the number of studies on gametocyte development that have been conducted, the molecular mechanisms regulating this process remain to be fully understood. This study investigates the functional roles of two female-specific transcriptional regulators, PbAP2-FG2 and PbAP2R-2, in P. berghei. Knockout of pbap2-fg2 or pbap2r-2 impairs female gametocyte development, resulting in developmental arrest during ookinete development. ChIP-seq analyses of these two factors indicated their colocalization on the genome, suggesting that they function as a complex. These analyses also revealed that their target genes contained a variety of genes, including both male and female-enriched genes. Moreover, differential expression analyses showed that these target genes were upregulated through the disruption of pbap2-fg2 or pbap2r-2, indicating that these two factors function as a transcriptional repressor complex in female gametocytes. Formation of a complex between PbAP2-FG2 and PbAP2R-2 was confirmed by RIME, a method that combines ChIP and MS analysis. In addition, the analysis identified a chromatin regulator PbMORC as an interaction partner of PbAP2-FG2. Comparative target analysis between PbAP2-FG2 and PbAP2-G demonstrated a significant overlap between their target genes, suggesting that repression of early gametocyte genes activated by PbAP2-G is one of the key roles for this female transcriptional repressor complex. Our results indicate that the PbAP2-FG2-PbAP2R-2 complex-mediated repression of the target genes supports the female differentiation from early gametocytes.

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.

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.

Identification of a novel AP2 transcription factor in zygotes with an essential role in Plasmodium ookinete development

The sexual phase of Plasmodium represents a crucial step in malaria transmission, during which these parasites fertilize and form ookinetes to infect mosquitoes. Plasmodium development after fertilization is thought to proceed with female-stored mRNAs until the formation of a retort-form ookinete; thus, transcriptional activity in zygotes has previously been considered quiescent. In this study, we reveal the essential role of transcriptional activity in zygotes by investigating the function of a newly identified AP2 transcription factor, AP2-Z, in P. berghei. ap2-z was previously reported as a female transcriptional regulator gene whose disruption resulted in developmental arrest at the retort stage of ookinetes. In this study, although ap2-z was transcribed in females, we show that it was translationally repressed by the DOZI complex and translated after fertilization with peak expression at the zygote stage. ChIP-seq analysis of AP2-Z shows that it binds on specific DNA motifs, targeting the majority of genes known as an essential component of ookinetes, which largely overlap with the AP2-O targets, as well as genes that are unique among the targets of other sexual transcription factors. The results of this study also indicate the existence of a cascade of transcription factors, beginning with AP2-G, that proceeds from gametocytogenesis to ookinete formation.

Sexual development in Plasmodium parasites, a causative agent of malaria, is essential for their transmission from vertebrate hosts to mosquitoes. This important developmental process proceeds as follows: formation of a gametocyte/gamete, fertilization and conversion of the zygote into the mosquito midgut invasive stage, called the ookinete. As a target of transmission blocking strategies, it is important to understand the mechanisms regulating Plasmodium sexual development. In this study, we assessed transcriptional regulation after fertilization by investigating the function of a novel transcription factor, AP2-Z. The results revealed the essential role of de novo transcription activated by AP2-Z in zygotes for promoting ookinete development. As transcriptional activity during the zygote stage has previously been considered silent in Plasmodium, novel genes important for ookinete formation can now be explored in the target genes of AP2-Z. Investigating the functions of these genes can help us understand the mechanisms of Plasmodium zygote/ookinete development and identify new targets for transmission blocking vaccines.

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.

Paving the Way: Contributions of Big Data to Apicomplexan and Kinetoplastid Research

In the age of big data an important question is how to ensure we make the most out of the resources we generate. In this review, we discuss the major methods used in Apicomplexan and Kinetoplastid research to produce big datasets and advance our understanding of Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania biology. We debate the benefits and limitations of the current technologies, and propose future advancements that may be key to improving our use of these techniques. Finally, we consider the difficulties the field faces when trying to make the most of the abundance of data that has already been, and will continue to be, generated.

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

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

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

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

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.

Targeting HDACs of apicomplexans: structural insights for a better treatment

Aetiologic agents of diseases such as malaria and toxoplasmosis are found in representatives of the phylum Apicomplexa. Therefore, apicomplexan parasites are known to have a significant impact on public health. Epigenetic factors such as histone acetylation/deacetylation are among the main mechanisms of gene regulation in these parasites. Histone deacetylases (HDACs) have aroused a great deal of interest over the past 20 years for being promising targets in the development of drugs for treating several diseases such as cancer. In addition, they have also been shown to be effective for parasitic diseases. However, little is known about the structure of these proteins, as well as their interactions with specific ligands. In this paper, we modelled 14 HDACs from different apicomplexan parasites and performed molecular docking with 12 ligands analogous to the HDAC inhibitors FR235222 and apicidin, which had previously been tested againstToxoplasma gondiiandPlasmodium falciparum. In thisin silicostudy, we were able to gather relevant structural data regarding these proteins as well as insights into protein–ligand interactions for testing and developing drugs for these diseases.

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.

Chromosome splitting of Plasmodium berghei using the CRISPR/Cas9 system

Spatial arrangement of chromosomes is responsible for gene expression in Plasmodium parasites. However, methods for rearranging chromosomes have not been established, which makes it difficult to investigate its role in detail. Here, we report a method for splitting chromosome in rodent malaria parasite by CRISPR/Cas9 system using fragments in which a telomere and a centromere were incorporated. The resultant split chromosomes segregated accurately into daughter parasites by the centromere. In addition, elongation of de novo telomeres were observed, indicating its proper function. Furthermore, chromosome splitting had no effect on development of parasites. Splitting of the chromosome is expected to alter its spatial arrangement, and our method will thus be useful for investigating its biological role related with gene expression.

Contribution of Transcriptome to Elucidate the Biology of Plasmodium spp

In the present review, we discuss some of the new technologies that have been applied to elucidate how Plasmodium spp escape from the immune system and subvert the host physiology to orchestrate the regulation of its biological pathways. Our manuscript describes how techniques such as microarray approaches, RNA-Seq and single-cell RNA sequencing have contributed to the discovery of transcripts and changed the concept of gene expression regulation in closely related malaria parasite species. Moreover, the text highlights the contributions of high-throughput RNA sequencing for the current knowledge of malaria parasite biology, physiology, vaccine target and the revelation of new players in parasite signaling.

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.

Pb103 Regulates Zygote/Ookinete Development in Plasmodium berghei via Double Zinc Finger Domains

Sexual reproduction of Plasmodium parasites takes place in anopheline mosquitoes, where male and female gametes fuse to form zygotes and then ookinetes. These processes are orchestrated by stage-specific protein expression, which is mediated in part by translational repression. Accumulating evidence shows that RNA binding proteins (RBPs) play crucial roles in these processes. Here, we report the characterization of P. berghei 103 (Pb103), which encodes a protein possessing double zinc finger domains (ZFs), an RBP. Reporter parasites expressing azami green fluorescent protein (AGFP) under the endogenous Pb103 gene promoter (Pb103-AGFP reporter) showed that the AGFP fluorescent signal was detected from gametes to ookinetes, while AGFP mRNA was translationally repressed in female gametocytes. The Pb103-disrupted parasites (Pb103(−)) grew and produced gametocytes with similar efficiencies to those of wild-type parasites. However, no oocysts were formed in mosquitoes fed Pb103(−). An in vitro fertilization assay showed abortion at the zygote stage in Pb103(−), suggesting that Pb103 plays a critical role in zygote/ookinete development. Cross-fertilization assays with Pb103(−) and male- or female-sterile parasites revealed that Pb103 was essential exclusively for female gametes. To identify the domains critical for zygote/ookinete development, transgenic parasites expressing partially deleted Pb103 were generated and assayed for ookinete maturation. As a result, deleting either of two ZFs but not the C-terminal region abolished zygote/ookinete development, highlighting the indispensable roles of ZFs in parasite sexual development, most likely via translational repression.

Dynamically expressed genes provide candidate viability biomarkers in a model coccidian

Eimeria parasites cause enteric disease in livestock and the closely related Cyclosporacayetanensis causes human disease. Oocysts of these coccidian parasites undergo maturation (sporulation) before becoming infectious. Here, we assessed transcription in maturing oocysts of Eimeria acervulina, a widespread chicken parasite, predicted gene functions, and determined which of these genes also occur in C. cayetanensis. RNA-Sequencing yielded ~2 billion paired-end reads, 92% of which mapped to the E. acervulina genome. The ~6,900 annotated genes underwent temporally-coordinated patterns of gene expression. Fifty-three genes each contributed >1,000 transcripts per million (TPM) throughout the study interval, including cation-transporting ATPases, an oocyst wall protein, a palmitoyltransferase, membrane proteins, and hypothetical proteins. These genes were enriched for 285 gene ontology (GO) terms and 13 genes were ascribed to 17 KEGG pathways, defining housekeeping processes and functions important throughout sporulation. Expression differed in mature and immature oocysts for 40% (2,928) of all genes; of these, nearly two-thirds (1,843) increased their expression over time. Eight genes expressed most in immature oocysts, encoding proteins promoting oocyst maturation and development, were assigned to 37 GO terms and 5 KEGG pathways. Fifty-six genes underwent significant upregulation in mature oocysts, each contributing at least 1,000 TPM. Of these, 40 were annotated by 215 GO assignments and 9 were associated with 18 KEGG pathways, encoding products involved in respiration, carbon fixation, energy utilization, invasion, motility, and stress and detoxification responses. Sporulation orchestrates coordinated changes in the expression of many genes, most especially those governing metabolic activity. Establishing the long-term fate of these transcripts in sporulated oocysts and in senescent and deceased oocysts will further elucidate the biology of coccidian development, and may provide tools to assay infectiousness of parasite cohorts. Moreover, because many of these genes have homologues in C. cayetanensis, they may prove useful as biomarkers for risk.

Single-Cell RNA Sequencing Reveals Cellular Heterogeneity and Stage Transition under Temperature Stress in Synchronized Plasmodium falciparum Cells

The malaria parasite has a complex life cycle exhibiting phenotypic and morphogenic variations in two different hosts by existing in heterogeneous developmental states. To investigate this cellular heterogeneity of the parasite within the human host, we performed single-cell RNA sequencing of synchronized Plasmodium cells under control and temperature treatment conditions. Using the Malaria Cell Atlas (https://www.sanger.ac.uk/science/tools/mca) as a guide, we identified 9 subtypes of the parasite distributed across known intraerythrocytic stages. Interestingly, temperature treatment results in the upregulation of the AP2-G gene, the master regulator of sexual development in a small subpopulation of the parasites. Moreover, we identified a heterogeneous stress-responsive subpopulation (clusters 5, 6, and 7 [∼10% of the total population]) that exhibits upregulation of stress response pathways under normal growth conditions. We also developed an online exploratory tool that will provide new insights into gene function under normal and temperature stress conditions. Thus, our study reveals important insights into cell-to-cell heterogeneity in the parasite population under temperature treatment that will be instrumental toward a mechanistic understanding of cellular adaptation and population dynamics in Plasmodium falciparum. IMPORTANCE The malaria parasite has a complex life cycle exhibiting phenotypic variations in two different hosts accompanied by cell-to-cell variability that is important for stress tolerance, immune evasion, and drug resistance. To investigate cellular heterogeneity determined by gene expression, we performed single-cell RNA sequencing (scRNA-seq) of about 12,000 synchronized Plasmodium cells under physiologically relevant normal (37°C) and temperature stress (40°C) conditions phenocopying the cyclic bouts of fever experienced during malarial infection. In this study, we found that parasites exhibit transcriptional heterogeneity in an otherwise morphologically synchronized culture. Also, a subset of parasites is continually committed to gametocytogenesis and stress-responsive pathways. These observations have important implications for understanding the mechanisms of drug resistance generation and vaccine development against the malaria parasite.

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.

Gene regulation of the avian malaria parasite Plasmodium relictum, during the different stages within the mosquito vector

The malaria parasite Plasmodium relictum is one of the most widespread species of avian malaria. As in the case of its human counterparts, bird Plasmodium undergoes a complex life cycle infecting two hosts: the arthropod vector and the vertebrate host. In this study, we examined transcriptomes of P. relictum (SGS1) during crucial timepoints within its vector, Culex pipiens quinquefasciatus. Differential gene-expression analyses identified genes linked to the parasites life-stages at: i) a few minutes after the blood meal is ingested, ii) during peak oocyst production phase, iii) during peak sporozoite phase and iv) during the late-stages of the infection. A large amount of genes coding for functions linked to host-immune invasion and multifunctional genes was active throughout the infection cycle. One gene associated with a conserved Plasmodium membrane protein with unknown function was upregulated throughout the parasite development in the vector, suggesting an important role in the successful completion of the sporogonic cycle. Gene expression analysis further identified genes, with unknown functions to be significantly differentially expressed during the infection in the vector as well as upregulation of reticulocyte-binding proteins, which raises the possibility of the multifunctionality of these RBPs. We establish the existence of highly stage-specific pathways being overexpressed during the infection. This first study of gene-expression of a non-human Plasmodium species in its vector provides a comprehensive insight into the molecular mechanisms of the common avian malaria parasite P. relictum and provides essential information on the evolutionary diversity in gene regulation of the Plasmodium's vector stages.

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

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.

Pleiotropic Roles for the Plasmodium berghei RNA Binding Protein UIS12 in Transmission and Oocyst Maturation

Colonization of the mosquito host by Plasmodium parasites is achieved by sexually differentiated gametocytes. Gametocytogenesis, gamete formation and fertilization are tightly regulated processes, and translational repression is a major regulatory mechanism for stage conversion. Here, we present a characterization of a Plasmodium berghei RNA binding protein, UIS12, that contains two conserved eukaryotic RNA recognition motifs (RRM). Targeted gene deletion resulted in viable parasites that replicate normally during blood infection, but form fewer gametocytes. Upon transmission to Anopheles stephensi mosquitoes, both numbers and size of midgut-associated oocysts were reduced and their development stopped at an early time point. As a consequence, no salivary gland sporozoites were formed indicative of a complete life cycle arrest in the mosquito vector. Comparative transcript profiling in mutant and wild-type infected red blood cells revealed a decrease in transcript abundance of mRNAs coding for signature gamete-, ookinete-, and oocyst-specific proteins in uis12(-) parasites. Together, our findings indicate multiple roles for UIS12 in regulation of gene expression after blood infection in good agreement with the pleiotropic defects that terminate successful sporogony and onward transmission to a new vertebrate host.

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 Comparison of Stage Conversion in the Coccidian Apicomplexans Toxoplasma gondii, Hammondia hammondi, and Neospora caninum

Stage conversion is a critical life cycle feature for several Apicomplexan parasites as the ability to switch between life forms is critical for replication, dissemination, pathogenesis and ultimately, transmission to a new host. In order for these developmental transitions to occur, the parasite must first sense changes in their environment, such as the presence of stressors or other environmental signals, and then respond to these signals by initiating global alterations in gene expression. As our understanding of the genetic components required for stage conversion continues to broaden, we can better understand the conserved mechanisms for this process and unique components and their contribution to pathogenesis by comparing stage conversion in multiple closely related species. In this review, we will discuss what is currently known about the mechanisms driving stage conversion in Toxoplasma gondii and its closest relatives Hammondia hammondi and Neospora caninum. Work by us and others has shown that these species have some important differences in the way that they (1) progress through their life cycle and (2) respond to stage conversion initiating stressors. To provide a specific example of species-specific complexities associated with stage conversion, we will discuss our recent published and unpublished work comparing stress responses in T. gondii and H. hammondi.

Comparative analysis of gene expression between Babesia bovis blood stages and kinetes allowed by improved genome annotation

Throughout their life cycle, Babesia parasites alternate between a mammalian host, where they cause babesiosis, and the tick vector. Transition between hosts results in distinct environmental signals that influence patterns of gene expression, consistent with the morphological and functional changes operating in the parasites during their life stages. In addition, comparing differential patterns of gene expression among mammalian and tick parasite stages can provide clues for developing improved methods of control. Hereby, we upgraded the genome assembly of Babesia bovis, a bovine hemoparasite, closing a 139 kbp gap, and used RNA-Seq datasets derived from mammalian blood and tick kinete stages to update the genome annotation. Of the originally annotated genes, 1,254 required structural changes, and 326 new genes were identified, leading to a different predicted proteome compared to the original annotation. Next, the RNA-Seq data was used to identify B. bovis genes that were differentially expressed in the vertebrate and arthropod hosts. In blood stages, 28% of the genes were upregulated up to 300 fold, whereas 26% of the genes in kinetes, a tick stage, were upregulated up to >19,000 fold. We thus discovered differentially expressed genes that may play key biological roles and serve as suitable targets for the development of vaccines to control bovine babesiosis.

Analysis of Plasmodium vivax schizont transcriptomes from field isolates reveals heterogeneity of expression of genes involved in host-parasite interactions

Plasmodium vivax gene regulation remains difficult to study due to the lack of a robust in vitro culture method, low parasite densities in peripheral circulation and asynchronous parasite development. We adapted an RNA-seq protocol “DAFT-seq” to sequence the transcriptome of four P. vivax field isolates that were cultured for a short period ex vivo before using a density gradient for schizont enrichment. Transcription was detected from 78% of the PvP01 reference genome, despite being schizont-enriched samples. This extensive data was used to define thousands of 5′ and 3′ untranslated regions, some of which overlapped with neighbouring transcripts, and to improve the gene models of 352 genes, including identifying 20 novel gene transcripts. This dataset has also significantly increased the known amount of heterogeneity between P. vivax schizont transcriptomes from individual patients. The majority of genes found to be differentially expressed between the isolates lack Plasmodium falciparum homologs and are predicted to be involved in host-parasite interactions, with an enrichment in reticulocyte binding proteins, merozoite surface proteins and exported proteins with unknown function. An improved understanding of the diversity within P. vivax transcriptomes will be essential for the prioritisation of novel vaccine targets.

Genome wide distribution of G-quadruplexes and their impact on gene expression in malaria parasites

Mechanisms of transcriptional control in malaria parasites are still not fully understood. The positioning patterns of G-quadruplex (G4) DNA motifs in the parasite's AT-rich genome, especially within the var gene family which encodes virulence factors, and in the vicinity of recombination hotspots, points towards a possible regulatory role of G4 in gene expression and genome stability. Here, we carried out the most comprehensive genome-wide survey, to date, of G4s in the Plasmodium falciparum genome using G4Hunter, which identifies G4 forming sequences (G4FS) considering their G-richness and G-skewness. We show an enrichment of G4FS in nucleosome-depleted regions and in the first exon of var genes, a pattern that is conserved within the closely related Laverania Plasmodium parasites. Under G4-stabilizing conditions, i.e., following treatment with pyridostatin (a high affinity G4 ligand), we show that a bona fide G4 found in the non-coding strand of var promoters modulates reporter gene expression. Furthermore, transcriptional profiling of pyridostatin-treated parasites, shows large scale perturbations, with deregulation affecting for instance the ApiAP2 family of transcription factors and genes involved in ribosome biogenesis. Overall, our study highlights G4s as important DNA secondary structures with a role in Plasmodium gene expression regulation, sub-telomeric recombination and var gene biology.

An intrinsic oscillator drives the blood stage cycle of the malaria parasite Plasmodium falciparum

The blood stage of the infection of the malaria parasite Plasmodium falciparum exhibits a 48-hour developmental cycle that culminates in the synchronous release of parasites from red blood cells, which triggers 48-hour fever cycles in the host. This cycle could be driven extrinsically by host circadian processes or by a parasite-intrinsic oscillator. To distinguish between these hypotheses, we examine the P. falciparum cycle in an in vitro culture system and show that the parasite has molecular signatures associated with circadian and cell cycle oscillators. Each of the four strains examined has a different period, which indicates strain-intrinsic period control. Finally, we demonstrate that parasites have low cell-to-cell variance in cycle period, on par with a circadian oscillator. We conclude that an intrinsic oscillator maintains Plasmodium's rhythmic life cycle.

Refining the transcriptome of the human malaria parasite Plasmodium falciparum using amplification-free RNA-seq

Background

Plasmodium parasites undergo several major developmental transitions during their complex lifecycle, which are enabled by precisely ordered gene expression programs. Transcriptomes from the 48-h blood stages of the major human malaria parasite Plasmodium falciparum have been described using cDNA microarrays and RNA-seq, but these assays have not always performed well within non-coding regions, where the AT-content is often 90–95%.

Results

We developed a directional, amplification-free RNA-seq protocol (DAFT-seq) to reduce bias against AT-rich cDNA, which we have applied to three strains of P. falciparum (3D7, HB3 and IT). While strain-specific differences were detected, overall there is strong conservation between the transcriptional profiles. For the 3D7 reference strain, transcription was detected from 89% of the genome, with over 78% of the genome transcribed into mRNAs. We also find that transcription from bidirectional promoters frequently results in non-coding, antisense transcripts. These datasets allowed us to refine the 5′ and 3′ untranslated regions (UTRs), which can be variable, long (> 1000 nt), and often overlap those of adjacent transcripts.

Conclusions

The approaches applied in this study allow a refined description of the transcriptional landscape of P. falciparum and demonstrate that very little of the densely packed P. falciparum genome is inactive or redundant. By capturing the 5′ and 3′ ends of mRNAs, we reveal both constant and dynamic use of transcriptional start sites across the intraerythrocytic developmental cycle that will be useful in guiding the definition of regulatory regions for use in future experimental gene expression studies.

Messenger RNAs with large numbers of upstream open reading frames are translated via leaky scanning and reinitiation in the asexual stages of Plasmodium falciparum

The genome of Plasmodium falciparum has one of the most skewed base-pair compositions of any eukaryote, with an AT content of 80–90%. As start and stop codons are AT-rich, the probability of finding upstream open reading frames (uORFs) in messenger RNAs (mRNAs) is high and parasite mRNAs have an average of 11 uORFs in their leader sequences. Similar to other eukaryotes, uORFs repress the translation of the downstream open reading frame (dORF) in P. falciparum, yet the parasite translation machinery is able to bypass these uORFs and reach the dORF to initiate translation. This can happen by leaky scanning and/or reinitiation.

In this report, we assessed leaky scanning and reinitiation by studying the effect of uORFs on the translation of a dORF, in this case, the luciferase reporter gene, and showed that both mechanisms are employed in the asexual blood stages of P. falciparum. Furthermore, in addition to the codon usage of the uORF, translation of the dORF is governed by the Kozak sequence and length of the uORF, and inter-cistronic distance between the uORF and dORF. Based on these features whole-genome data was analysed to uncover classes of genes that might be regulated by uORFs. This study indicates that leaky scanning and reinitiation appear to be widespread in asexual stages of P. falciparum, which may require modifications of existing factors that are involved in translation initiation in addition to novel, parasite-specific proteins.

Single-cell transcription analysis of Plasmodium vivax blood-stage parasites identifies stage- and species-specific profiles of expression

Plasmodium vivax and P. falciparum, the parasites responsible for most human malaria worldwide, exhibit striking biological differences, which have important clinical consequences. Unfortunately, P. vivax, unlike P. falciparum, cannot be cultivated continuously in vitro, which limits our understanding of its biology and, consequently, our ability to effectively control vivax malaria. Here, we describe single-cell gene expression profiles of 9,215 P. vivax parasites from bloodstream infections of Aotus and Saimiri monkeys. Our results show that transcription of most P. vivax genes occurs during short periods of the intraerythrocytic cycle and that this pattern of gene expression is conserved in other Plasmodium species. However, we also identify a strikingly high proportion of species-specific transcripts in late schizonts, possibly associated with the specificity of erythrocyte invasion. Our findings provide new and robust markers of blood-stage parasites, including some that are specific to the elusive P. vivax male gametocytes, and will be useful for analyzing gene expression data from laboratory and field samples.

Analysis of individual Plasmodium vivax parasites reveals the tight control of the expression of most genes during the intra-erythrocytic cycle and the differentiation of male and female gametocytes, and highlights differences between the development of P. vivax and P. falciparum.

The second life of Plasmodium in the mosquito host: gene regulation on the move

Malaria parasites face dynamically changing environments and strong selective constraints within human and mosquito hosts. To survive such hostile and shifting conditions, Plasmodium switches transcriptional programs during development and has evolved mechanisms to adjust its phenotype through heterogeneous patterns of gene expression. In vitro studies on culture-adapted isolates have served to set the link between chromatin structure and functional gene expression. Yet, experimental evidence is limited to certain stages of the parasite in the vertebrate, i.e. blood, while the precise mechanisms underlying the dynamic regulatory landscapes during development and in the adaptation to within-host conditions remain poorly understood. In this review, we discuss available data on transcriptional and epigenetic regulation in Plasmodium mosquito stages in the context of sporogonic development and phenotypic variation, including both bet-hedging and environmentally triggered direct transcriptional responses. With this, we advocate the mosquito offers an in vivo biological model to investigate the regulatory networks, transcription factors and chromatin-modifying enzymes and their modes of interaction with regulatory sequences, which might be responsible for the plasticity of the Plasmodium genome that dictates stage- and cell type-specific blueprints of gene expression.

Efficient Single-Gene and Gene Family Editing in the Apicomplexan Parasite Eimeria tenella Using CRISPR-Cas9

Eimeria species are pathogenic protozoa with a wide range of hosts and the cause of poultry coccidiosis, which results in huge economic losses to the poultry industry. These parasites encode a genome of ∼8000 genes that control a highly coordinated life cycle of asexual replication and sexual differentiation, transmission, and virulence. However, the function and physiological importance of the large majority of these genes remain unknown mostly due to the lack of tools for systematic analysis of gene functions. Here, we report the first application of CRISPR-Cas9 gene editing technology in Eimeria tenella for analysis of gene function at a single gene level as well as for systematic functional analysis of an entire gene family. Using a transgenic line constitutively expressing Cas9, we demonstrated successful and efficient loss of function through non-homologous end joining as well as guided homologous recombination. Application of this approach to the study of the localization of EtGRA9 revealed that the gene encodes a secreted protein whose cellular distribution varied during the life cycle. Systematic disruption of the ApiAp2 transcription factor gene family using this approach revealed that 23 of the 33 factors expressed by this parasite are essential for development and survival in the host. Our data thus establish CRISPR-Cas9 as a powerful technology for gene editing in Eimeria and will set the stage for systematic functional analysis of its genome to understand its biology and pathogenesis, and will make it possible to identify and validate new targets for coccidiosis therapy.

A single-nucleotide polymorphism in a Plasmodium berghei ApiAP2 transcription factor alters the development of host immunity

The acquisition of malaria immunity is both remarkably slow and unpredictable. At present, we know little about the malaria parasite genes that influence the host's ability to mount a protective immune response. Here, we show that a single-nucleotide polymorphism (SNP) resulting in a single amino acid change (S to F) in an ApiAP2 transcription factor in the rodent malaria parasite Plasmodium berghei (Pb) NK65 allowed infected mice to mount a T helper cell 1 (T_H1)-type immune response that controlled subsequent infections. As compared to PbNK65^S, PbNK65^F parasites differentially expressed 46 genes, most of which are predicted to play roles in immune evasion. PbNK65^F infections resulted in an early interferon-γ response and a later expansion of germinal centers, resulting in high levels of infected red blood cell-specific T_H1-type immunoglobulin G2b (IgG2b) and IgG2c antibodies. Thus, the Pb ApiAP2 transcription factor functions as a critical parasite virulence factor in malaria infections.

Identification of a Master Regulator of Differentiation in Toxoplasma

Toxoplasma gondii chronically infects a quarter of the world's population, and its recrudescence can cause life-threatening disease in immunocompromised individuals and recurrent ocular lesions in the immunocompetent. Acute-stage tachyzoites differentiate into chronic-stage bradyzoites, which form intracellular cysts resistant to immune clearance and existing therapies. The molecular basis of this differentiation is unknown, despite being efficiently triggered by stresses in culture. Through Cas9-mediated screening and single-cell profiling, we identify a Myb-like transcription factor (BFD1) necessary for differentiation in cell culture and in mice. BFD1 accumulates during stress and its synthetic expression is sufficient to drive differentiation. Consistent with its function as a transcription factor, BFD1 binds the promoters of many stage-specific genes and represents a counterpoint to the ApiAP2 factors that dominate our current view of parasite gene regulation. BFD1 provides a genetic switch to study and control Toxoplasma differentiation and will inform prevention and treatment of chronic infections.

G-strand binding protein 2 is involved in asexual and sexual development of Plasmodium berghei

G-strand binding protein 2 (GBP2) is a Ser/Arg-rich (SR) protein involved in mRNA surveillance and nuclear mRNA quality control in yeast. However, the roles of GBP2 in virulence and sexual development in Plasmodium parasites are unclear, although GBP2 is involved in the asexual development of Plasmodium berghei, the rodent malaria parasite. In this study, we investigated the role of GBP2 in virulence and sexual development of P. berghei using gbp2-deleted P. berghei (Δgbp2 parasites). Then, to identify factors affected by gbp2 deletion, we performed a comparative proteomic analysis of the Δgbp2 parasites. We found that GBP2 was not associated with the development of experimental cerebral malaria during infection with P. berghei, but asexual development of the parasite was delayed with deletion of gbp2. However, the development of P. berghei gametocytes was significantly reduced with deletion of gbp2. Comparative proteomic analysis revealed that the levels of adenosine deaminase (ADA), purine nucleoside phosphorylase (PNP), and hypoxanthine-guanine phosphoribosyltransferase (HGPRT) in Δgbp2 parasites were significantly higher than those in wild-type (WT) parasites, suggesting that biosynthesis of purine nucleotides may be involved in function of GBP2. Therefore, we investigated the effect of purine starvation on the sexual development and proteome. In nt1-deleted P. berghei (Δnt1 parasites), the production of male and female gametocytes was significantly reduced compared to that in WT parasites. Moreover, we found that protein levels of GBP2 in Δnt1 parasites were markedly lower than in WT parasites. These findings suggest that GBP2 is primarily involved in the sexual development of malaria parasites, and its function may be suppressed by purine starvation.

Female-specific gene regulation in malaria parasites by an AP2-family transcription factor

The malaria gametocyte, the gamete precursor, is the essential stage for malaria transmission to the mosquito vector. In the vertebrate host's blood, it develops into a mature male or female capable of transforming into a gamete in the mosquito blood meal. Despite the importance of this stage in the malaria life cycle, the genetic regulation of gametocyte development is poorly understood. In particular, transcription factors involved in sex-specific gene expression have not been identified. In this paper, we report that an AP2-family transcription factor, AP2-FG, is responsible for female-specific gene regulation. AP2-FG expression in Plasmodium berghei was observed exclusively in female gametocytes, in the beginning of 4-6 h before sexual dimorphism manifests in developing gametocytes. AP2-FG disruption resulted in the arrest of female maturation, but did not affect the development of males. Chromatin immunoprecipitation sequencing analysis suggested that AP2-FG directly regulates over 700 genes. Its targets include genes for female gametocyte-specific functions, such as gametogenesis, fertilization and zygote development. AP2-FG binding to target gene promoters was associated with a 10 bp sequence motif. These results indicate that AP2-FG plays a role in the differentiation of early gametocytes to mature females by governing a female-specific gene expression repertoire.

Hierarchical transcriptional control regulates Plasmodium falciparum sexual differentiation

BACKGROUND

Malaria pathogenesis relies on sexual gametocyte forms of the malaria parasite to be transmitted between the infected human and the mosquito host but the molecular mechanisms controlling gametocytogenesis remains poorly understood. Here we provide a high-resolution transcriptome of Plasmodium falciparum as it commits to and develops through gametocytogenesis.

RESULTS

The gametocyte-associated transcriptome is significantly different from that of the asexual parasites, with dynamic gene expression shifts characterizing early, intermediate and late-stage gametocyte development and results in differential timing for sex-specific transcripts. The transcriptional dynamics suggest strict transcriptional control during gametocytogenesis in P. falciparum, which we propose is mediated by putative regulators including epigenetic mechanisms (driving active repression of proliferation-associated processes) and a cascade-like expression of ApiAP2 transcription factors.

CONCLUSIONS

The gametocyte transcriptome serves as the blueprint for sexual differentiation and will be a rich resource for future functional studies on this critical stage of Plasmodium development, as the intraerythrocytic transcriptome has been for our understanding of the asexual cycle.

Protozoan persister-like cells and drug treatment failure

Antimicrobial treatment failure threatens our ability to control infections. In addition to antimicrobial resistance, treatment failures are increasingly understood to derive from cells that survive drug treatment without selection of genetically heritable mutations. Parasitic protozoa, such as Plasmodium species that cause malaria, Toxoplasma gondii and kinetoplastid protozoa, including Trypanosoma cruzi and Leishmania spp., cause millions of deaths globally. These organisms can evolve drug resistance and they also exhibit phenotypic diversity, including the formation of quiescent or dormant forms that contribute to the establishment of long-term infections that are refractory to drug treatment, which we refer to as 'persister-like cells'. In this Review, we discuss protozoan persister-like cells that have been linked to persistent infections and discuss their impact on therapeutic outcomes following drug treatment.

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

Malaria parasites are characterized by a complex life cycle that is accompanied by dynamic gene expression patterns. The factors and mechanisms that regulate gene expression in these parasites have been searched for even before the advent of next generation sequencing technologies. Functional genomics approaches have substantially boosted this area of research and have yielded significant insights into the interplay between epigenetic, transcriptional and post-transcriptional mechanisms. Recently, considerable progress has been made in identifying sequence-specific transcription factors and DNA-encoded regulatory elements. Here, we review the insights obtained from these efforts including the characterization of core promoters, the involvement of sequence-specific transcription factors in life cycle progression and the mapping of gene regulatory elements. Furthermore, we discuss recent developments in the field of functional genomics and how they might contribute to further characterization of this complex gene regulatory network.

Predicting gene expression in the human malaria parasite Plasmodium falciparum using histone modification, nucleosome positioning, and 3D localization features

Empirical evidence suggests that the malaria parasite Plasmodium falciparum employs a broad range of mechanisms to regulate gene transcription throughout the organism's complex life cycle. To better understand this regulatory machinery, we assembled a rich collection of genomic and epigenomic data sets, including information about transcription factor (TF) binding motifs, patterns of covalent histone modifications, nucleosome occupancy, GC content, and global 3D genome architecture. We used these data to train machine learning models to discriminate between high-expression and low-expression genes, focusing on three distinct stages of the red blood cell phase of the Plasmodium life cycle. Our results highlight the importance of histone modifications and 3D chromatin architecture in Plasmodium transcriptional regulation and suggest that AP2 transcription factors may play a limited regulatory role, perhaps operating in conjunction with epigenetic factors.