The Cryptosporidium parvum ApiAP2 gene family: insights into the evolution of apicomplexan AP2 regulatory systems

Single cell expression and chromatin accessibility of the Toxoplasma gondii lytic cycle identifies AP2XII-8 as an essential ribosome regulon driver

Sequential lytic cycles driven by cascading transcriptional waves underlie pathogenesis in the apicomplexan parasite Toxoplasma gondii. This parasite's unique division by internal budding, short cell cycle, and jumbled up classically defined cell cycle stages have restrained in-depth transcriptional program analysis. Here, unbiased transcriptome and chromatin accessibility maps throughout the lytic cell cycle are established at the single-cell level. Correlated pseudo-timeline assemblies of expression and chromatin profiles maps transcriptional versus chromatin level transition points promoting the cell division cycle. Sequential clustering analysis identifies functionally related gene groups promoting cell cycle progression. Promoter DNA motif mapping reveals patterns of combinatorial regulation. Pseudo-time trajectory analysis reveals transcriptional bursts at different cell cycle points. The dominant burst in G1 is driven largely by transcription factor AP2XII-8, which engages a conserved DNA motif, and promotes the expression of 44 ribosomal proteins encoding regulon. Overall, the study provides integrated, multi-level insights into apicomplexan transcriptional regulation.

Cryptosporidium life cycle small molecule probing implicates translational repression and an Apetala 2 transcription factor in macrogamont differentiation

The apicomplexan parasite Cryptosporidium is a leading cause of childhood diarrhea in developing countries. Current treatment options are inadequate and multiple preclinical compounds are being actively pursued as potential drugs for cryptosporidiosis. Unlike most apicomplexans, Cryptosporidium spp. sequentially replicate asexually and then sexually within a single host to complete their lifecycles. Anti-cryptosporidial compounds are generally identified or tested through in vitro phenotypic assays that only assess the asexual stages. Therefore, compounds that specifically target the sexual stages remain unexplored. In this study, we leveraged the ReFRAME drug repurposing library against a newly devised multi-readout imaging assay to identify small-molecule compounds that modulate macrogamont differentiation and maturation. RNA-seq studies confirmed selective modulation of macrogamont differentiation for 10 identified compounds (9 inhibitors and 1 accelerator). The collective transcriptomic profiles of these compounds indicates that translational repression accompanies Cryptosporidium sexual differentiation, which we validated experimentally. Additionally, cross comparison of the RNA-seq data with promoter sequence analysis for stage-specific genes converged on a key role for an Apetala 2 (AP2) transcription factor (cgd2_3490) in differentiation into macrogamonts. Finally, drug annotation for the ReFRAME hits indicates that an elevated supply of energy equivalence in the host cell is critical for macrogamont formation.

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

Introduction

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

Methods

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

Results

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

Discussion

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

Cell cycle-regulated ApiAP2s and parasite development: the Toxoplasma paradigm

The cell division cycle of T. gondii is driven by cyclically expressed ApiAP2 transcription factors (AP2s) that promote gene sets (regulons) associated with specific biological functions. AP2s drive other AP2s, thereby propelling the progressive gene expression waves defining the lytic cycle. AP2s can act as dimers by themselves, in combination with other AP2s (constitutive or cyclical) or in complexes with epigenetic factors. Exit from the cell cycle into either the extracellular state or differentiation into bradyzoites results in major changes in gene expression. Surprisingly, both transitions lead to expression of a shared set of unique AP2s that suggest a shared stress response that, governed by the specific conditions, leads to different outcomes.

Babesia gibsoni Whole-Genome Sequencing, Assembling, Annotation, and Comparative Analysis

The intracellular protozoan parasite Babesia gibsoni infects canine erythrocytes and causes babesiosis. The hazards to animal health have increased due to the rise of B. gibsoni infections and medication resistance. However, the lack of high-quality full-genome sequencing sets has expanded the obstacles to the development of pathogeneses, drugs, and vaccines. In this study, the whole genome of B. gibsoni was sequenced, assembled, and annotated. The genomic size of B. gibsoni was 7.94 Mbp in total. Four chromosomes with the size of 0.69 Mb, 2.10 Mb, 2.77 Mb, and 2.38 Mb, respectively, 1 apicoplast (28.4 Kb), and 1 mitochondrion (5.9 Kb) were confirmed. KEGG analysis revealed 2,641 putative proteins enriched on 316 pathways, and GO analysis showed 7,571 annotations of the nuclear genome in total. Synteny analysis showed a high correlation between B. gibsoni and B. bovis. A new divergent point of B. gibsoni occurred around 297.7 million years ago, which was earlier than that of B. bovis, B. ovata, and B. bigemina. Orthology analysis revealed 22 and 32 unique genes compared to several Babesia spp. and apicomplexan species. The metabolic pathways of B.gibsoni were characterized, pointing to a minimal size of the genome. A species-specific secretory protein SA1 and 19 homologous genes were identified. Selected specific proteins, including apetala 2 (AP2) factor, invasion-related proteins BgAMA-1 and BgRON2, and rhoptry function proteins BgWH_04g00700 were predicted, visualized, and modeled. Overall, whole-genome sequencing provided molecular-level support for the diagnosis, prevention, clinical treatment, and further research of B. gibsoni. IMPORTANCE The whole genome of B. gibsoni was first sequenced, annotated, and disclosed. The key part of genome composition, four chromosomes, was comparatively analyzed for the first time. A full-scale phylogeny evolution analysis based on the whole-genome-wide data of B. gibsoni was performed, and a new divergent point on the evolutionary path was revealed. In previous reports, molecular studies were often limited by incomplete genomic data, especially in key areas like life cycle regulation, metabolism, and host-pathogen interaction. With the whole-genome sequencing of B. gibsoni, we provide useful genetic data to encourage the exploration of new terrain and make it feasible to resolve the theoretical and practical problems of babesiosis.

Genetic Ablation of a Female-Specific Apetala 2 Transcription Factor Blocks Oocyst Shedding in Cryptosporidium parvum

The apicomplexan parasite Cryptosporidium is a leading global cause of diarrheal disease, and the infection poses a particularly grave threat to young children and those with weakened immune function. Infection occurs by ingestion of meiotic spores called oocysts, and transmission relies on fecal shedding of new oocysts. The entire life cycle thus occurs in a single host and features asexual as well as sexual forms of replication. Here, we identify and locus tag two Apetala 2-type (AP2) transcription factors and demonstrate that they are exclusively expressed in male and female gametes, respectively. To enable functional studies of essential genes in Cryptosporidium parvum, we develop and validate a small-molecule-inducible gene excision system, which we apply to the female factor AP2-F to achieve conditional gene knockout. Analyzing this mutant, we find the factor to be dispensable for asexual growth and early female fate determination in vitro but to be required for oocyst shedding in infected animals in vivo. Transcriptional analyses conducted in the presence or absence of AP2-F revealed that the factor controls the transcription of genes encoding crystalloid body proteins, which are exclusively expressed in female gametes. In C. parvum, the organelle is restricted to sporozoites, and its loss in other apicomplexan parasites leads to blocked transmission. Overall, our development of conditional gene ablation in C. parvum provides a robust method for genetic analysis in this parasite that enabled us to identify AP2-F as an essential regulator of transcription required for oocyst shedding and transmission. IMPORTANCE The parasite Cryptosporidium infects millions of people worldwide each year, leading to life-threatening diarrheal disease in young children and immunosuppressed individuals. There is no vaccine and only limited treatment. Transmission occurs via the fecal-oral route by an environmentally resilient spore-like oocyst. Infection takes place in the intestinal epithelium, where parasites initially propagate asexually before transitioning to male and female gametes, with sex leading to the formation of new oocysts. The essential role of sexual development for continuous infection and transmission makes it an attractive target for therapy and prevention. To study essential genes and potential drug targets across the life cycle, we established inducible gene excision for C. parvum. We determined that the female-specific transcription factor AP2-F is not required for asexual growth and early female development in vitro but is necessary for oocyst shedding in vivo. This work enhances the genetic tools available to study Cryptosporidium gene function.

Putative SET-domain methyltransferases in Cryptosporidium parvum and histone methylation during infection

Cryptosporidium parvum is a leading cause of diarrhoeal illness worldwide being a significant threat to young children and immunocompromised patients, but the pathogenesis caused by this parasite remains poorly understood. C. parvum was recently linked with oncogenesis. Notably, the mechanisms of gene expression regulation are unexplored in Cryptosporidium and little is known about how the parasite impact host genome regulation. Here, we investigated potential histone lysine methylation, a dynamic epigenetic modification, during the life cycle of the parasite. We identified SET-domain containing proteins, putative lysine methyltransferases (KMTs), in the C. parvum genome and classified them phylogenetically into distinct subfamilies (namely CpSET1, CpSET2, CpSET8, CpKMTox and CpAKMT). Our structural analysis further characterized CpSET1, CpSET2 and CpSET8 as histone lysine methyltransferases (HKMTs). The expression of the CpSET genes varies considerably during the parasite life cycle and specific methyl-lysine antibodies showed dynamic changes in parasite histone methylation during development (CpSET1:H3K4; CpSET2:H3K36; CpSET8:H4K20). We investigated the impact of C. parvum infection on the host histone lysine methylation. Remarkably, parasite infection led to a considerable decrease in host H3K36me3 and H3K27me3 levels, highlighting the potential of the parasite to exploit the host epigenetic regulation to its advantage. This is the first study to describe epigenetic mechanisms occurring throughout the parasite life cycle and during the host–parasite interaction. A better understanding of histone methylation in both parasite and host genomes may highlight novel infection control strategies.

Division and Transmission: Malaria Parasite Development in the Mosquito

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

Evolution of the DEHYDRATION-RESPONSIVE ELEMENT-BINDING PROTEIN subfamily in green plants

Adapting to unfavorable environments is a necessary step in plant terrestrialization and radiation. The dehydration-responsive element-binding (DREB) protein subfamily plays a pivotal role in plant abiotic stress regulation. However, relationships between the origin and expansion of the DREB subfamily and adaptive evolution of land plants are still being elucidated. Here, we constructed the evolutionary history of the DREB subfamily by compiling APETALA2/ethylene-responsive element-binding protein superfamily genes from 169 representative species of green plants. Through extensive phylogenetic analyses and comparative genomic analysis, our results revealed that the DREB subfamily diverged from the ethylene-responsive factor (ERF) subfamily in the common ancestor of Zygnemophyceae and Embryophyta during the colonization of land by plants, followed by expansions to form three different ancient archetypal genes in Zygnemophyceae species, designated as groups archetype-I, archetype-II/III, and archetype-IV. Four large-scale expansions paralleling the evolution of land plants led to the nine-subgroup divergence of group archetype-II/III in angiosperms, and five whole-genome duplications during Brassicaceae and Poaceae radiation shaped the diversity of subgroup IIb-1. We identified a Poaceae-specific gene in subgroup IIb-1, ERF014, remaining in a Poaceae-specific microsynteny block and co-evolving with a small heat shock protein cluster. Expression analyses demonstrated that heat acclimation may have driven the neofunctionalization of ERF014s in Pooideae by engaging in the conserved heat-responsive module in Poaceae. This study provides insights into lineage-specific expansion and neofunctionalization in the DREB subfamily, together with evolutionary information valuable for future functional studies of plant stress biology.

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.

The developmental trajectories of Toxoplasma stem from an elaborate epigenetic rewiring

Toxoplasma gondii is considered to be one of the most successful parasitic pathogens. It owes this success to its flexibility in responding to signals emanating from the different environments it encounters during its multihost life cycle. The adaptability of this unicellular organism relies on highly coordinated and evolutionarily optimized developmental abilities that allow it to adopt the forms best suited to the requirements of each environment. Here we discuss recent outstanding studies that have uncovered how master regulators epigenetically regulate the cryptic process of sexual development and the transition to chronicity. We also highlight the molecular and technical advances that allow the field to embark on a new journey of epigenetic reprogramming of T. gondii development.

Global transcriptome landscape of the rabbit protozoan parasite Eimeria stiedae

Background

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

Methods

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

Results

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

Conclusions

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

Graphic Abstract

Supplementary Information

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

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.

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

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

Analysis of Long Non-Coding RNA in Cryptosporidium parvum Reveals Significant Stage-Specific Antisense Transcription

Cryptosporidium is a protist parasite that has been identified as the second leading cause of moderate to severe diarrhea in children younger than two and a significant cause of mortality worldwide. Cryptosporidium has a complex, obligate, intracellular but extra cytoplasmic lifecycle in a single host. How genes are regulated in this parasite remains largely unknown. Long non-coding RNAs (lncRNAs) play critical regulatory roles, including gene expression across a broad range of organisms. Cryptosporidium lncRNAs have been reported to enter the host cell nucleus and affect the host response. However, no systematic study of lncRNAs in Cryptosporidium has been conducted to identify additional lncRNAs. In this study, we analyzed a C. parvum in vitro strand-specific RNA-seq developmental time series covering both asexual and sexual stages to identify lncRNAs associated with parasite development. In total, we identified 396 novel lncRNAs, mostly antisense, with 86% being differentially expressed. Surprisingly, nearly 10% of annotated mRNAs have an antisense transcript. lncRNAs occur most often at the 3' end of their corresponding sense mRNA. Putative lncRNA regulatory regions were identified and many appear to encode bidirectional promoters. A positive correlation between lncRNA and upstream mRNA expression was observed. Evolutionary conservation and expression of lncRNA candidates was observed between C. parvum, C. hominis and C. baileyi. Ten C. parvum protein-encoding genes with antisense transcripts have P. falciparum orthologs that also have antisense transcripts. Three C. parvum lncRNAs with exceptional properties (e.g., intron splicing) were experimentally validated using RT-PCR and RT-qPCR. This initial characterization of the C. parvum non-coding transcriptome facilitates further investigations into the roles of lncRNAs in parasite development and host-pathogen interactions.

Essential Genes of the Parasitic Apicomplexa

Genome-scale mutagenesis screens for genes essential for apicomplexan parasite survival have been completed in three species: Plasmodium falciparum, the major human malaria parasite, Plasmodium berghei, a model rodent malaria parasite, and the more distantly related Toxoplasma gondii, the causative agent of toxoplasmosis. These three species share 2606 single-copy orthologs, 1500 of which have essentiality data in all three screens. In this review, we explore the overlap between these datasets to define the core essential genes of the phylum Apicomplexa. We further discuss the implications of these groundbreaking studies for understanding apicomplexan parasite biology, and we identify promising areas of focus for developing new pan-apicomplexan parasite interventions.

Internal daughter formation of Toxoplasma gondii tachyzoites is coordinated by transcription factor TgAP2IX-5

Toxoplasma gondii rapidly propagates through endodyogeny of tachyzoites, a process in which daughter parasites divide within the cell of the mother parasite. Recent studies have revealed that transcription factors with AP2-domain participate in the process of cell division in T. gondii. However, the concise regulation of the division cycles by AP2 proteins is poorly understood. In this study, we evaluated the effect of the transcription factor TgAP2IX-5 on the daughter cell formation in T. gondii. TgAP2IX-5 is a nuclear protein and is highly expressed during the S phase of the cell cycle of tachyzoites. TgAP2IX-5-disrupted strain showed a severe defect in replication and completely blocked lytic parasite growth. Following 3-indoleacetic acid treatment or without treatment of AP2IX-5-AID-3HA tagged strain for 30 min, 1 and 2 hr, the differentially expressed genes were 8, 54 and 202, respectively. Among these genes, the significantly downregulated ones were AP2 proteins, inner membrane complex (IMC) proteins and SAG-related proteins. Interestingly, loss of TgAP2IX-5 leads to a defect in internal daughter IMC formation and abnormalities in the morphology of organelles during cell division. Together, our study suggests that TgAP2IX-5 is crucial in regulating IMC formation of daughter cells in T. gondii.

Cryptosporidium: host and parasite transcriptome in infection

Cryptosporidium is a waterborne gastrointestinal parasite that causes outbreaks of diarrheal disease worldwide. Despite the impact of this parasite on human health there are no effective drugs or vaccines. Transcriptomic data can provide insights into host-parasite interactions that lead to identification of targets for therapeutic interventions. However, for Cryptosporidium, interpreting transcriptomes has been challenging, in part due to the presence of multiple life cycle stages, the lack of appropriate host cells and the inability to culture the parasite through its complete life cycle. The recent improvements in cell culture and the ability to tag and isolate specific life cycle stages will radically improve transcriptomic data and advance our understanding of Cryptosporidium host-parasite interactions.

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.

Life cycle progression and sexual development of the apicomplexan parasite Cryptosporidium parvum

The apicomplexan parasite Cryptosporidium is a leading global cause of severe diarrhoeal disease and an important contributor to early childhood mortality. Currently, there are no fully effective treatments or vaccines available. Parasite transmission occurs through ingestion of oocysts, through either direct contact or consumption of contaminated water or food. Oocysts are meiotic spores and the product of parasite sex. Cryptosporidium has a single-host life cycle in which both asexual and sexual processes occur in the intestine of infected hosts. Here, we genetically engineered strains of Cryptosporidium to make life cycle progression and parasite sex tractable. We derive reporter strains to follow parasite development in culture and in infected mice and define the genes that orchestrate sex and oocyst formation through mRNA sequencing of sorted cells. After 2 d, parasites in cell culture show pronounced sexualization, but productive fertilization does not occur and infection falters. By contrast, in infected mice, male gametes successfully fertilize female parasites, which leads to meiotic division and sporulation. To rigorously test for fertilization, we devised a two-component genetic-crossing assay using a reporter that is activated by Cre recombinase. Our findings suggest obligate developmental progression towards sex in Cryptosporidium, which has important implications for the treatment and prevention of the infection.

Infection with Cryptosporidium parvum is a leading cause of severe diarrhoeal disease and childhood mortality worldwide. Using tools they recently developed to genetically engineer Cryptosporidium, the authors define life cycle stage-specific markers and generate reporter parasites, making life cycle progression and parasite sex tractable.

ApiAP2 Transcription Factors in Apicomplexan Parasites

Apicomplexan parasites are protozoan organisms that are characterised by complex life cycles and they include medically important species, such as the malaria parasite Plasmodium and the causative agents of toxoplasmosis (Toxoplasma gondii) and cryptosporidiosis (Cryptosporidium spp.). Apicomplexan parasites can infect one or more hosts, in which they differentiate into several morphologically and metabolically distinct life cycle stages. These developmental transitions rely on changes in gene expression. In the last few years, the important roles of different members of the ApiAP2 transcription factor family in regulating life cycle transitions and other aspects of parasite biology have become apparent. Here, we review recent progress in our understanding of the different members of the ApiAP2 transcription factor family in apicomplexan parasites.

The Epigenome, Cell Cycle, and Development in Toxoplasma

Toxoplasma gondii is a common veterinary and human pathogen that persists as latent bradyzoite forms within infected hosts. The ability of the parasite to interconvert between tachyzoite and bradyzoite is key for pathogenesis of toxoplasmosis, particularly in immunocompromised individuals. The transition between tachyzoites and bradyzoites is epigenetically regulated and coupled to the cell cycle. Recent epigenomic studies have begun to elucidate the chromatin states associated with developmental switches in T. gondii. Evidence is also emerging that AP2 transcription factors both activate and repress the bradyzoite developmental program. Further studies are needed to understand the mechanisms by which T. gondii transduces environmental signals to coordinate the epigenetic and transcriptional machinery that are responsible for tachyzoite-bradyzoite interconversion.

A Knockout Screen of ApiAP2 Genes Reveals Networks of Interacting Transcriptional Regulators Controlling the Plasmodium Life Cycle

A family of apicomplexa-specific proteins containing AP2 DNA-binding domains (ApiAP2s) was identified in malaria parasites. This family includes sequence-specific transcription factors that are key regulators of development. However, functions for the majority of ApiAP2 genes remain unknown. Here, a systematic knockout screen in Plasmodium berghei identified ten ApiAP2 genes that were essential for mosquito transmission: four were critical for the formation of infectious ookinetes, and three were required for sporogony. We describe non-essential functions for AP2-O and AP2-SP proteins in blood stages, and identify AP2-G2 as a repressor active in both asexual and sexual stages. Comparative transcriptomics across mutants and developmental stages revealed clusters of co-regulated genes with shared cis promoter elements, whose expression can be controlled positively or negatively by different ApiAP2 factors. We propose that stage-specific interactions between ApiAP2 proteins on partly overlapping sets of target genes generate the complex transcriptional network that controls the Plasmodium life cycle.

•

Mutants in 11 of 26 apiAP2 genes reveal gene functions in mosquito transmission

•

Co-expression clustering across mutants and stages reveals molecular phenotypes

•

Multifunctional apiAP2 genes create complex regulatory networks in Plasmodium

•

Ap2-g2 is a transcriptional repressor in both asexual and sexual blood stages

Comparative transcriptome analysis of second- and third-generation merozoites of Eimeria necatrix

Background

Eimeria is a common genus of apicomplexan parasites that infect diverse vertebrates, most notably poultry, causing serious disease and economic losses. Eimeria species have complex life-cycles consisting of three developmental stages. However, the molecular basis of the Eimeria reproductive mode switch remains an enigma.

Methods

Total RNA extracted from second- (MZ-2) and third-generation merozoites (MZ-3) of Eimeria necatrix was subjected to transcriptome analysis using RNA sequencing (RNA-seq) followed by qRT-PCR validation.

Results

A total of 6977 and 6901 unigenes were obtained from MZ-2 and MZ-3, respectively. Approximately 2053 genes were differentially expressed genes (DEGs) between MZ-2 and MZ-3. Compared with MZ-2, 837 genes were upregulated and 1216 genes were downregulated in MZ-3. Approximately 95 genes in MZ-2 and 48 genes in MZ-3 were further identified to have stage-specific expression. Gene ontology category and KEGG analysis suggested that 216 upregulated genes in MZ-2 were annotated by 70 GO assignments, 242 upregulated genes were associated with 188 signal pathways, while 321 upregulated genes in MZ-3 were annotated by 56 GO assignments, 322 upregulated genes were associated with 168 signal pathways. The molecular functions of upregulated genes in MZ-2 were mainly enriched for protein degradation and amino acid metabolism. The molecular functions of upregulated genes in MZ-3 were mainly enriched for transcriptional activity, cell proliferation and cell differentiation.

Conclusions

To the best of our knowledge, this is the first RNA-seq data study of the MZ-2 and MZ-3 stages of E. necatrix; it demonstrates a high number of differentially expressed genes between the MZ-2 and MZ-3 of E. necatrix. This study forms a basis for deciphering the molecular mechanisms underlying the shift from the second to third generation schizogony in Eimeria. It also provides valuable resources for future studies on Eimeria, and provides insight into the understanding of reproductive mode plasticity in different Eimeria species.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2325-z) contains supplementary material, which is available to authorized users.

Genomics of apicomplexan parasites

The increasing prevalence of infections involving intracellular apicomplexan parasites such as Plasmodium, Toxoplasma, and Cryptosporidium (the causative agents of malaria, toxoplasmosis, and cryptosporidiosis, respectively) represent a significant global healthcare burden. Despite their significance, few treatments are available; a situation that is likely to deteriorate with the emergence of new resistant strains of parasites. To lay the foundation for programs of drug discovery and vaccine development, genome sequences for many of these organisms have been generated, together with large-scale expression and proteomic datasets. Comparative analyses of these datasets are beginning to identify the molecular innovations supporting both conserved processes mediating fundamental roles in parasite survival and persistence, as well as lineage-specific adaptations associated with divergent life-cycle strategies. The challenge is how best to exploit these data to derive insights into parasite virulence and identify those genes representing the most amenable targets. In this review, we outline genomic datasets currently available for apicomplexans and discuss biological insights that have emerged as a consequence of their analysis. Of particular interest are systems-based resources, focusing on areas of metabolism and host invasion that are opening up opportunities for discovering new therapeutic targets.

Toxoplasma gondii AP2IX-4 Regulates Gene Expression during Bradyzoite Development

Toxoplasma gondii is a single-celled parasite that persists in its host as a transmissible tissue cyst. How the parasite converts from its replicative form to the bradyzoites housed in tissue cysts is not well understood, but the process clearly involves changes in gene expression. Here we report that parasites lacking a cell cycle-regulated transcription factor called AP2IX-4 display reduced frequencies of tissue cyst formation in culture and in a mouse model of infection. Parasites missing AP2IX-4 lose the ability to regulate bradyzoite genes during tissue cyst development. Expressed in developing bradyzoites still undergoing division, AP2IX-4 may serve as a useful marker in the study of transitional forms of the parasite.

Toxoplasma gondii is a protozoan parasite of great importance to human and animal health. In the host, this obligate intracellular parasite persists as a tissue cyst that is imperceptible to the immune response and unaffected by current therapies. The tissue cysts facilitate transmission through predation and give rise to chronic cycles of toxoplasmosis in immunocompromised patients. Transcriptional changes accompany conversion of the rapidly replicating tachyzoites into the encysted bradyzoites, and yet the mechanisms underlying these alterations in gene expression are not well defined. Here we show that AP2IX-4 is a nuclear protein exclusively expressed in tachyzoites and bradyzoites undergoing division. Knockout of AP2IX-4 had no discernible effect on tachyzoite replication but resulted in a reduced frequency of tissue cyst formation following alkaline stress induction—a defect that is reversible by complementation. AP2IX-4 has a complex role in regulating bradyzoite gene expression, as the levels of many bradyzoite mRNAs dramatically increased beyond those seen under conditions of normal stress induction in AP2IX-4 knockout parasites exposed to alkaline media. The loss of AP2IX-4 also resulted in a modest virulence defect and reduced cyst burden in chronically infected mice, which was reversed by complementation. These findings illustrate that the transcriptional mechanisms responsible for tissue cyst development operate across the intermediate life cycle from the dividing tachyzoite to the dormant bradyzoite.

IMPORTANCEToxoplasma gondii is a single-celled parasite that persists in its host as a transmissible tissue cyst. How the parasite converts from its replicative form to the bradyzoites housed in tissue cysts is not well understood, but the process clearly involves changes in gene expression. Here we report that parasites lacking a cell cycle-regulated transcription factor called AP2IX-4 display reduced frequencies of tissue cyst formation in culture and in a mouse model of infection. Parasites missing AP2IX-4 lose the ability to regulate bradyzoite genes during tissue cyst development. Expressed in developing bradyzoites still undergoing division, AP2IX-4 may serve as a useful marker in the study of transitional forms of the parasite.

Opposing Transcriptional Mechanisms Regulate Toxoplasma Development

Toxoplasma infections are lifelong because of the development of the bradyzoite tissue cyst, which is effectively invisible to the immune system. Despite the important clinical consequences of this developmental pathway, the molecular basis of the switch mechanisms that control tissue cyst formation is still poorly understood. Significant changes in gene expression are associated with tissue cyst development, and ApiAP2 transcription factors are an important mechanism regulating this developmental transcriptome. However, the molecular composition of these ApiAP2 complexes and the operating principles of ApiAP2 mechanisms are not well defined. Here we establish that competing ApiAP2 transcriptional mechanisms operate to regulate this clinically important developmental pathway.

The Toxoplasma biology that underlies human chronic infection is developmental conversion of the acute tachyzoite stage into the latent bradyzoite stage. We investigated the roles of two alkaline-stress-induced ApiAP2 transcription factors, AP2IV-3 and AP2IX-9, in bradyzoite development. These factors were expressed in two overlapping waves during bradyzoite development, with AP2IX-9 increasing expression earlier than AP2IV-3, which peaked as AP2IX-9 expression was declining. Disruption of the AP2IX-9 gene enhanced, while deletion of AP2IV-3 gene decreased, tissue cyst formation, demonstrating that these factors have opposite functions in bradyzoite development. Conversely, conditional overexpression of FKBP-modified AP2IX-9 or AP2IV-3 with the small molecule Shield 1 had a reciprocal effect on tissue cyst formation, confirming the conclusions of the knockout experiments. The AP2IX-9 repressor and AP2IV-3 activator tissue cyst phenotypes were borne out in gene expression studies that determined that many of the same bradyzoite genes were regulated in an opposite manner by these transcription factors. A common gene target was the canonical bradyzoite marker BAG1, and mechanistic experiments determined that, like AP2IX-9, AP2IV-3 regulates a BAG1 promoter-luciferase reporter and specifically binds the BAG1 promoter in parasite chromatin. Altogether, these results suggest that the AP2IX-9 transcriptional repressor and the AP2IV-3 transcriptional activator likely compete to control bradyzoite gene expression, which may permit Toxoplasma to better adapt to different tissue environments and select a suitable host cell for long-term survival of the dormant tissue cyst.

IMPORTANCEToxoplasma infections are lifelong because of the development of the bradyzoite tissue cyst, which is effectively invisible to the immune system. Despite the important clinical consequences of this developmental pathway, the molecular basis of the switch mechanisms that control tissue cyst formation is still poorly understood. Significant changes in gene expression are associated with tissue cyst development, and ApiAP2 transcription factors are an important mechanism regulating this developmental transcriptome. However, the molecular composition of these ApiAP2 complexes and the operating principles of ApiAP2 mechanisms are not well defined. Here we establish that competing ApiAP2 transcriptional mechanisms operate to regulate this clinically important developmental pathway.

Comparative Bioinformatics Analysis of Transcription Factor Genes Indicates Conservation of Key Regulatory Domains among Babesia bovis, Babesia microti, and Theileria equi

Apicomplexa tick-borne hemoparasites, including Babesia bovis, Babesia microti, and Theileria equi are responsible for bovine and human babesiosis and equine theileriosis, respectively. These parasites of vast medical, epidemiological, and economic impact have complex life cycles in their vertebrate and tick hosts. Large gaps in knowledge concerning the mechanisms used by these parasites for gene regulation remain. Regulatory genes coding for DNA binding proteins such as members of the Api-AP2, HMG, and Myb families are known to play crucial roles as transcription factors. Although the repertoire of Api-AP2 has been defined and a HMG gene was previously identified in the B. bovis genome, these regulatory genes have not been described in detail in B. microti and T. equi. In this study, comparative bioinformatics was used to: (i) identify and map genes encoding for these transcription factors among three parasites' genomes; (ii) identify a previously unreported HMG gene in B. microti; (iii) define a repertoire of eight conserved Myb genes; and (iv) identify AP2 correlates among B. bovis and the better-studied Plasmodium parasites. Searching the available transcriptome of B. bovis defined patterns of transcription of these three gene families in B. bovis erythrocyte stage parasites. Sequence comparisons show conservation of functional domains and general architecture in the AP2, Myb, and HMG proteins, which may be significant for the regulation of common critical parasite life cycle transitions in B. bovis, B. microti, and T. equi. A detailed understanding of the role of gene families encoding DNA binding proteins will provide new tools for unraveling regulatory mechanisms involved in B. bovis, B. microti, and T. equi life cycles and environmental adaptive responses and potentially contributes to the development of novel convergent strategies for improved control of babesiosis and equine piroplasmosis.

Comparative genomics reveals Cyclospora cayetanensis possesses coccidia-like metabolism and invasion components but unique surface antigens

Background

Cyclospora cayetanensis is an apicomplexan that causes diarrhea in humans. The investigation of foodborne outbreaks of cyclosporiasis has been hampered by a lack of genetic data and poor understanding of pathogen biology. In this study we sequenced the genome of C. cayetanensis and inferred its metabolism and invasion components based on comparative genomic analysis.

Results

The genome organization, metabolic capabilities and potential invasion mechanism of C. cayetanensis are very similar to those of Eimeria tenella. Propanoyl-CoA degradation, GPI anchor biosynthesis, and N-glycosylation are some apparent metabolic differences between C. cayetanensis and E. tenella. Unlike Eimeria spp., there are no active LTR-retrotransposons identified in C. cayetanensis. The similar repertoire of host cell invasion-related proteins possessed by all coccidia suggests that C. cayetanensis has an invasion process similar to the one in T. gondii and E. tenella. However, the significant reduction in the number of identifiable rhoptry protein kinases, phosphatases and serine protease inhibitors indicates that monoxenous coccidia, especially C. cayetanensis, have limited capabilities or use a different system to regulate host cell nuclear activities. C. cayetanensis does not possess any cluster of genes encoding the TA4-type SAG surface antigens seen in E. tenella, and may use a different family of surface antigens in initial host cell interactions.

Conclusions

Our findings indicate that C. cayetanensis possesses coccidia-like metabolism and invasion components but unique surface antigens. Amino acid metabolism and post-translation modifications of proteins are some major differences between C. cayetanensis and other apicomplexans. The whole genome sequence data of C. cayetanensis improve our understanding of the biology and evolution of this major foodborne pathogen and facilitate the development of intervention measures and advanced diagnostic tools.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2632-3) contains supplementary material, which is available to authorized users.

Cis regulatory motifs and antisense transcriptional control in the apicomplexan Theileria parva

Background

Theileria parva is an intracellular parasite that causes a lymphoproliferative disease in cattle. It does so by inducing cancer-like phenotypes in the host cells it infects, although the molecular and regulatory mechanisms involved remain poorly understood. RNAseq data, and the resulting updated genome annotation now available for this parasite, offer an unprecedented opportunity to characterize the genomic features associated with gene regulation in this species. Our previous analyses revealed a T. parva genome even more gene-dense than previously thought, with many adjacent loci overlapping each other, not only at the level of untranslated sequences (UTRs) but even in coding sequences.

Results

Despite this compactness, Theileria intergenic regions show a pattern of size distribution indicative of monocistronic gene transcription. Three previously described motifs are conserved among Theileria species and highly prevalent in promoter regions near or at the transcription start sites. We found novel motifs at many transcription termination sites, as well as upstream of parasite genes thought to be critical for host transformation. Adjacent genes that could be regulated by antisense transcription from an overlapping transcriptional unit are syntenic between T. parva and P. falciparum at a frequency higher than expected by chance, suggesting the presence of common, and evolutionary old, regulatory mechanisms in the phylum Apicomplexa.

Conclusions

We propose a model of transcription with conserved sense and antisense transcription from a few taxonomically ubiquitous and several species-specific promoter motifs. Interestingly, the gene networks regulated by conserved promoters are themselves, in most cases, not conserved between species or genera.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2444-5) contains supplementary material, which is available to authorized users.

ApiAP2 Factors as Candidate Regulators of Stochastic Commitment to Merozoite Production in Theileria annulata

Background

Differentiation of one life-cycle stage to the next is critical for survival and transmission of apicomplexan parasites. A number of studies have shown that stage differentiation is a stochastic process and is associated with a point that commits the cell to a change over in the pattern of gene expression. Studies on differentiation to merozoite production (merogony) in T. annulata postulated that commitment involves a concentration threshold of DNA binding proteins and an auto-regulatory loop.

Principal Findings

In this study ApiAP2 DNA binding proteins that show changes in expression level during merogony of T. annulata have been identified. DNA motifs bound by orthologous domains in Plasmodium were found to be enriched in upstream regions of stage-regulated T. annulata genes and validated as targets for the T. annulata AP2 domains by electrophoretic mobility shift assay (EMSA). Two findings were of particular note: the gene in T. annulata encoding the orthologue of the ApiAP2 domain in the AP2-G factor that commits Plasmodium to gametocyte production, has an expression profile indicating involvement in transmission of T. annulata to the tick vector; genes encoding related domains that bind, or are predicted to bind, sequence motifs of the type 5'-(A)CACAC(A) are implicated in differential regulation of gene expression, with one gene (TA11145) likely to be preferentially up-regulated via auto-regulation as the cell progresses to merogony.

Conclusions

We postulate that the Theileria factor possessing the AP2 domain orthologous to that of Plasmodium AP2-G may regulate gametocytogenesis in a similar manner to AP2-G. In addition, paralogous ApiAP2 factors that recognise 5'-(A)CACAC(A) type motifs could operate in a competitive manner to promote reversible progression towards the point that commits the cell to undergo merogony. Factors possessing AP2 domains that bind (or are predicted to bind) this motif are present in the vector-borne genera Theileria, Babesia and Plasmodium, and other Apicomplexa; leading to the proposal that the mechanisms that control stage differentiation will show a degree of conservation.

The ability of vector-borne Apicomplexan parasites (Babesia, Plasmodium and Theileria) to change from one life-cycle stage to the next is critical for establishment of infection and transmission to new hosts. Stage differentiation steps of both Plasmodium and Theileria are known to involve stochastic transition through an intermediate form to a point that commits the cell to generate the next stage in the life-cycle. In this study we have identified genes encoding ApiAP2 DNA binding proteins in Theileria annulata that are differentially expressed during differentiation from the macroschizont stage, through merozoite production (merogony) to the piroplasm stage. The results provide evidence that the ApiAp2 factor in Theileria that possesses the orthologue of the Plasmodium AP2-G domain may also operate to regulate gametocytogenesis, and that progression to merogony is promoted by the ability of a merozoite DNA binding protein to preferentially up-regulate its own production. In addition, identification of multiple ApiAP2 DNA binding domains that bind related motifs within and across vector-borne Apicomplexan genera lead to the proposal that the mechanisms that promote the transition from asexual to sexual replication will show a degree of conservation.

UniPROBE, update 2015: new tools and content for the online database of protein-binding microarray data on protein-DNA interactions

The Universal PBM Resource for Oligonucleotide Binding Evaluation (UniPROBE) serves as a convenient source of information on published data generated using universal protein-binding microarray (PBM) technology, which provides in vitro data about the relative DNA-binding preferences of transcription factors for all possible sequence variants of a length k (‘k-mers’). The database displays important information about the proteins and displays their DNA-binding specificity data in terms of k-mers, position weight matrices and graphical sequence logos. This update to the database documents the growth of UniPROBE since the last update 4 years ago, and introduces a variety of new features and tools, including a new streamlined pipeline that facilitates data deposition by universal PBM data generators in the research community, a tool that generates putative nonbinding (i.e. negative control) DNA sequences for one or more proteins and novel motifs obtained by analyzing the PBM data using the BEEML-PBM algorithm for motif inference. The UniPROBE database is available at http://uniprobe.org.