Coordinated progression through two subtranscriptomes underlies the tachyzoite cycle of Toxoplasma gondii

BCC0 collaborates with IMC32 and IMC43 to form the Toxoplasma gondii essential daughter bud assembly complex

Toxoplasma gondii divides by endodyogeny, in which two daughter buds are formed within the cytoplasm of the maternal cell using the inner membrane complex (IMC) as a scaffold. During endodyogeny, components of the IMC are synthesized and added sequentially to the nascent daughter buds in a tightly regulated manner. We previously showed that the early recruiting proteins IMC32 and IMC43 form an essential daughter bud assembly complex which lays the foundation of the daughter cell scaffold in T. gondii. In this study, we identify the essential, early recruiting IMC protein BCC0 as a third member of this complex by using IMC32 as bait in both proximity labeling and yeast two-hybrid screens. We demonstrate that BCC0's localization to daughter buds depends on the presence of both IMC32 and IMC43. Deletion analyses and functional complementation studies reveal that residues 701-877 of BCC0 are essential for both its localization and function and that residues 1-899 are sufficient for function despite minor mislocalization. Pairwise yeast two-hybrid assays additionally demonstrate that BCC0's essential domain binds to the coiled-coil region of IMC32 and that BCC0 and IMC43 do not directly interact. This data supports a model for complex assembly in which an IMC32-BCC0 subcomplex initially recruits to nascent buds via palmitoylation of IMC32 and is locked into the scaffold once bud elongation begins by IMC32 binding to IMC43. Together, this study dissects the organization and function of a complex of three early recruiting daughter proteins which are essential for the proper assembly of the IMC during endodyogeny.

Integrating patients in time series clinical transcriptomics data

MOTIVATION

Analysis of time series transcriptomics data from clinical trials is challenging. Such studies usually profile very few time points from several individuals with varying response patterns and dynamics. Current methods for these datasets are mainly based on linear, global orderings using visit times which do not account for the varying response rates and subgroups within a patient cohort.

RESULTS

We developed a new method that utilizes multi-commodity flow algorithms for trajectory inference in large scale clinical studies. Recovered trajectories satisfy individual-based timing restrictions while integrating data from multiple patients. Testing the method on multiple drug datasets demonstrated an improved performance compared to prior approaches suggested for this task, while identifying novel disease subtypes that correspond to heterogeneous patient response patterns.

AVAILABILITY AND IMPLEMENTATION

The source code and instructions to download the data have been deposited on GitHub at https://github.com/euxhenh/Truffle.

Temporal gene expression during asexual development of the apicomplexan Sarcocystis neurona

Asexual replication in the apicomplexan Sarcocystis neurona involves two main developmental stages: the motile extracellular merozoite and the sessile intracellular schizont. Merozoites invade host cells and transform into schizonts that undergo replication via endopolygeny to form multiple (64) daughter merozoites that are invasive to new host cells. Given that the capabilities of the merozoite vary significantly from the schizont, the patterns of transcript levels throughout the asexual lifecycle were determined and compared in this study. RNA-Seq data were generated from extracellular merozoites and four intracellular schizont development time points. Of the 6,938 genes annotated in the S. neurona genome, 6,784 were identified in the transcriptome. Of these, 4,111 genes exhibited significant differential expression between the merozoite and at least one schizont development time point. Transcript levels were significantly higher for 2,338 genes in the merozoite and 1,773 genes in the schizont stages. Included in this list were genes encoding the secretory pathogenesis determinants (SPDs), which encompass the surface antigen and SAG-related sequence (SAG/SRS) and the secretory organelle proteins of the invasive zoite stage (micronemes, rhoptries, and dense granules). As anticipated, many of the S. neurona SPD gene transcripts were abundant in merozoites. However, several SPD transcripts were elevated in intracellular schizonts, suggesting roles unrelated to host cell invasion and the initial establishment of the intracellular niche. The hypothetical genes that are potentially unique to the genus Sarcocystis are of particular interest. Their conserved expression patterns are instructive for future investigations into the possible functions of these putative Sarcocystis-unique genes.

IMPORTANCE

The genus Sarcocystis is an expansive clade within the Apicomplexa, with the species S. neurona being an important cause of neurological disease in horses. Research to decipher the biology of S. neurona and its host-pathogen interactions can be enhanced by gene expression data. This study has identified conserved apicomplexan orthologs in S. neurona, putative Sarcocystis-unique genes, and gene transcripts abundant in the merozoite and schizont stages. Importantly, we have identified distinct clusters of genes with transcript levels peaking during different intracellular schizont development time points, reflecting active gene expression changes across endopolygeny. Each cluster also has subsets of transcripts with unknown functions, and investigation of these seemingly Sarcocystis-unique transcripts will provide insights into the interesting biology of this parasite genus.

The Crk4-Cyc4 complex regulates G2/M transition in Toxoplasma gondii

A versatile division of apicomplexan parasites and a dearth of conserved regulators have hindered the progress of apicomplexan cell cycle studies. While most apicomplexans divide in a multinuclear fashion, Toxoplasma gondii tachyzoites divide in the traditional binary mode. We previously identified five Toxoplasma CDK-related kinases (Crk). Here, we investigated TgCrk4 and its cyclin partner TgCyc4. We demonstrated that TgCrk4 regulates conventional G2 phase processes, such as repression of chromosome rereplication and centrosome reduplication, and acts upstream of the spindle assembly checkpoint. The spatial TgCyc4 dynamics supported the TgCrk4-TgCyc4 complex role in the coordination of chromosome and centrosome cycles. We also identified a dominant TgCrk4-TgCyc4 complex interactor, TgiRD1 protein, related to DNA replication licensing factor CDT1 but played no role in licensing DNA replication in the G1 phase. Our results showed that TgiRD1 also plays a role in controlling chromosome and centrosome reduplication. Global phosphoproteome analyses identified TgCrk4 substrates, including TgORC4, TgCdc20, TgGCP2, and TgPP2ACA. Importantly, the phylogenetic and structural studies suggest the Crk4-Cyc4 complex is limited to a minor group of the binary dividing apicomplexans.

AP2XII-1 and AP2XI-2 Suppress Schizogony Gene Expression in Toxoplasma gondii

Toxoplasma gondii is an intracellular parasite that is important in medicine and veterinary science and undergoes distinct developmental transitions in its intermediate and definitive hosts. The switch between stages of T. gondii is meticulously regulated by a variety of factors. Previous studies have explored the role of the microrchidia (MORC) protein complex as a transcriptional suppressor of sexual commitment. By utilizing immunoprecipitation and mass spectrometry, constituents of this protein complex have been identified, including MORC, Histone Deacetylase 3 (HDAC3), and several ApiAP2 transcription factors. Conditional knockout of MORC or inhibition of HDAC3 results in upregulation of a set of genes associated with schizogony and sexual stages in T. gondii tachyzoites. Here, our focus extends to two primary ApiAP2s (AP2XII-1 and AP2XI-2), demonstrating their significant impact on the fitness of asexual tachyzoites and their target genes. Notably, the targeted disruption of AP2XII-1 and AP2XI-2 resulted in a profound alteration in merozoite-specific genes targeted by the MORC-HDAC3 complex. Additionally, considerable overlap was observed in downstream gene profiles between AP2XII-1 and AP2XI-2, with AP2XII-1 specifically binding to a subset of ApiAP2 transcription factors, including AP2XI-2. These findings reveal an intricate cascade of ApiAP2 regulatory networks involved in T. gondii schizogony development, orchestrated by AP2XII-1 and AP2XI-2. This study provides valuable insights into the transcriptional regulation of T. gondii growth and development, shedding light on the intricate life cycle of this parasitic pathogen.

The Toxoplasma gondii F-Box Protein L2 Functions as a Repressor of Stage Specific Gene Expression

Toxoplasma gondii is a foodborne pathogen that can cause severe and life-threatening infections in fetuses and immunocompromised patients. Felids are its only definitive hosts, and a wide range of animals, including humans, serve as intermediate hosts. When the transmissible bradyzoite stage is orally ingested by felids, they transform into merozoites that expand asexually, ultimately generating millions of gametes for the parasite sexual cycle. However, bradyzoites in intermediate hosts differentiate exclusively to disease-causing tachyzoites, which rapidly disseminate throughout the host. Though tachyzoites are well-studied, the molecular mechanisms governing transitioning between developmental stages are poorly understood. Each parasite stage can be distinguished by a characteristic transcriptional signature, with one signature being repressed during the other stages. Switching between stages require substantial changes in the proteome, which is achieved in part by ubiquitination. F-box proteins mediate protein poly-ubiquitination by recruiting substrates to SKP1, Cullin-1, F-Box protein E3 ubiquitin ligase (SCF-E3) complexes. We have identified an F-box protein named Toxoplasma gondii F-Box Protein L2 (TgFBXL2), which localizes to distinct perinucleolar sites. TgFBXL2 is stably engaged in an SCF-E3 complex that is surprisingly also associated with a COP9 signalosome complex that negatively regulates SCF-E3 function. At the cellular level, TgFBXL2-depleted parasites are severely defective in centrosome replication and daughter cell development. Most remarkable, RNAseq data show that TgFBXL2 conditional depletion induces the expression of stage-specific genes including a large cohort of genes necessary for sexual commitment. Together, these data suggest that TgFBXL2 is a latent guardian of stage specific gene expression in Toxoplasma and poised to remove conflicting proteins in response to an unknown trigger of development.

Systematic characterization of all Toxoplasma gondii TBC domain-containing proteins identifies an essential regulator of Rab2 in the secretory pathway

Toxoplasma gondii resides in its intracellular niche by employing a series of specialized secretory organelles that play roles in invasion, host cell manipulation, and parasite replication. Rab GTPases are major regulators of the parasite's secretory traffic that function as nucleotide-dependent molecular switches to control vesicle trafficking. While many of the Rab proteins have been characterized in T. gondii, precisely how these Rabs are regulated remains poorly understood. To better understand the parasite's secretory traffic, we investigated the entire family of Tre2-Bub2-Cdc16 (TBC) domain-containing proteins, which are known to be involved in vesicle fusion and secretory protein trafficking. We first determined the localization of all 18 TBC domain-containing proteins to discrete regions of the secretory pathway or other vesicles in the parasite. Second, we use an auxin-inducible degron approach to demonstrate that the protozoan-specific TgTBC9 protein, which localizes to the endoplasmic reticulum (ER), is essential for parasite survival. Knockdown of TgTBC9 results in parasite growth arrest and affects the organization of the ER and mitochondrial morphology. TgTBC9 knockdown also results in the formation of large lipid droplets (LDs) and multi-membranous structures surrounded by ER membranes, further indicating a disruption of ER functions. We show that the conserved dual-finger active site in the TBC domain of the protein is critical for its GTPase-activating protein (GAP) function and that the Plasmodium falciparum orthologue of TgTBC9 can rescue the lethal knockdown. We additionally show by immunoprecipitation and yeast 2 hybrid analyses that TgTBC9 preferentially binds Rab2, indicating that the TBC9-Rab2 pair controls ER morphology and vesicular trafficking in the parasite. Together, these studies identify the first essential TBC protein described in any protozoan and provide new insight into intracellular vesicle trafficking in T. gondii.

An apicoplast-localized deubiquitinase contributes to the cell growth and apicoplast homeostasis of Toxoplasma gondii

Toxoplasma gondii is among the most important parasites worldwide. The apicoplast is a unique organelle shared by all Apicomplexan protozoa. Increasing lines of evidence suggest that the apicoplast possesses its own ubiquitination system. Deubiquitination is a crucial step executed by deubiquitinase (DUB) during protein ubiquitination. While multiple components of ubiquitination have been identified in T. gondii, the deubiquitinases involved remain unknown. The aim of the current study was to delineate the localization of TgOTU7 and elucidate its functions. TgOTU7 was specifically localized at the apicoplast, and its expression was largely regulated during the cell cycle. Additionally, TgOTU7 efficiently breaks down ubiquitin chains, exhibits linkage-nonspecific deubiquitinating activity and is critical for the lytic cycle and apicoplast biogenesis, similar to the transcription of the apicoplast genome and the nuclear genes encoding apicoplast-targeted proteins. Taken together, the results indicate that the newly described deubiquitinase TgOTU7 specifically localizes to the apicoplast and affects the cell growth and apicoplast homeostasis of T. gondii.

Toxoplasma gondii F-Box Protein L2 Silences Feline-Restricted Genes Necessary for Sexual Commitment

Toxoplasma gondii is a foodborne pathogen that can cause severe and life-threatening infections in fetuses and immunocompromised patients. Felids are its only definitive hosts, and a wide range of animals, including humans, serve as intermediate hosts. When the transmissible bradyzoite stage is orally ingested by felids, they transform into merozoites that expand asexually, ultimately generating millions of gametes for the parasite sexual cycle. However, bradyzoites in intermediate hosts differentiate exclusively to disease-causing tachyzoites, which rapidly disseminate throughout the host. Though tachyzoites are well-studied, the molecular mechanisms governing transitioning between developmental stages are poorly understood. Each parasite stage can be distinguished by a characteristic transcriptional signature, with one signature being repressed during the other stages. Switching between stages requires substantial changes in the proteome, which is achieved in part by ubiquitination. F-box proteins mediate protein poly-ubiquitination by recruiting substrates to SKP1, Cullin-1, F-Box protein E3 ubiquitin ligase (SCF-E3) complexes. We have identified an F-box protein named Toxoplasma gondii F-Box Protein L2 (TgFBXL2), which localizes to distinct nuclear sites. TgFBXL2 is stably engaged in an SCF-E3 complex that is surprisingly also associated with a COP9 signalosome complex that negatively regulates SCF-E3 function. At the cellular level, TgFBXL2-depleted parasites are severely defective in centrosome replication and daughter cell development. Most remarkable, RNA seq data show that TgFBXL2 conditional depletion induces the expression of genes necessary for sexual commitment. We suggest that TgFBXL2 is a latent guardian of sexual stage development in Toxoplasma and poised to remove conflicting proteins in response to an unknown trigger of sexual development.

The multifaceted roles of Myb domain-containing proteins in apicomplexan parasites

Apicomplexan parasites are a large and diverse clade of protists responsible for significant diseases of humans and animals. Central to the ability of these parasites to colonize their host and evade immune responses is an expanded repertoire of gene-expression programs that requires the coordinated action of complex transcriptional networks. DNA-binding proteins and chromatin regulators are essential orchestrators of apicomplexan gene expression that often act in concert. Although apicomplexan genomes encode various families of putative DNA-binding proteins, most remain functionally and mechanistically unexplored. This review highlights the versatile role of myeloblastosis (Myb) domain-containing proteins in apicomplexan parasites as transcription factors and chromatin regulators. We explore the diversity of Myb domain structure and use phylogenetic analysis to identify common features across the phylum. This provides a framework to discuss functional heterogeneity and regulation of Myb domain-containing proteins particularly emphasizing their role in parasite differentiation.

The Toxoplasma protein phosphatase 6 catalytic subunit (TgPP6C) is essential for cell cycle progression and virulence

Protein phosphatases are post-translational regulators of Toxoplasma gondii proliferation, tachyzoite-bradyzoite differentiation and pathogenesis. Here, we identify the putative protein phosphatase 6 (TgPP6) subunits of T. gondii and elucidate their role in the parasite lytic cycle. The putative catalytic subunit TgPP6C and regulatory subunit TgPP6R likely form a complex whereas the predicted structural subunit TgPP6S, with low homology to the human PP6 structural subunit, does not coassemble with TgPP6C and TgPP6R. Functional studies showed that TgPP6C and TgPP6R are essential for parasite growth and replication. The ablation of TgPP6C significantly reduced the synchronous division of the parasite's daughter cells during endodyogeny, resulting in disordered rosettes. Moreover, the six conserved motifs of TgPP6C were required for efficient endodyogeny. Phosphoproteomic analysis revealed that ablation of TgPP6C predominately altered the phosphorylation status of proteins involved in the regulation of the parasite cell cycle. Deletion of TgPP6C significantly attenuated the parasite virulence in mice. Immunization of mice with TgPP6C-deficient type I RH strain induced protective immunity against challenge with a lethal dose of RH or PYS tachyzoites and Pru cysts. Taken together, the results show that TgPP6C contributes to the cell division, replication and pathogenicity in T. gondii.

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.

AP2XII-1 is a negative regulator of merogony and presexual commitment in Toxoplasma gondii

IMPORTANCE

Sexual development is vital for the transmission, genetic hybridization, and population evolution of apicomplexan pathogens, which include several clinically relevant parasites, such as Plasmodium, Eimeria, and Toxoplasma gondii. Previous studies have demonstrated different morphological characteristics and division patterns between asexual and sexual stages of the parasites. However, the primary regulation is poorly understood. A transition from the asexual to the sexual stage is supposedly triggered/accompanied by rewiring of gene expression and controlled by transcription factors and chromatin modulators. Herein, we discovered a tachyzoite-specific transcriptional factor AP2XII-1, which represses the presexual development in the asexual tachyzoite stage of T. gondii. Conditional knockdown of AP2XII-1 perturbs tachyzoite proliferation by endodyogeny and drives a transition to a morphologically and transcriptionally distinct merozoite stage. The results also suggest a hierarchical transcriptional regulation of sexual development by AP2 factors and provide a path to culturing merozoites and controlling inter-host transmission of T. gondii.

An ex vivo model of Toxoplasma recrudescence reveals developmental plasticity of the bradyzoite stage

IMPORTANCE

The classical depiction of the Toxoplasma lifecycle is bradyzoite excystation conversion to tachyzoites, cell lysis, and immune control, followed by the reestablishment of bradyzoites and cysts. In contrast, we show that tachyzoite growth slows independent of the host immune response at a predictable time point following excystation. Furthermore, we demonstrate a host cell-dependent pathway of continuous amplification of the cyst-forming bradyzoite population. The developmental plasticity of the excysted bradyzoites further underlines the critical role the cyst plays in the flexibility of the lifecycle of this ubiquitous parasite. This revised model of Toxoplasma recrudescence uncovers previously unknown complexity in the clinically important bradyzoite stage of the parasite, which opens the door to further study these novel developmental features of the Toxoplasma intermediate life cycle.

High-throughput identification of Toxoplasma gondii effector proteins that target host cell transcription

Intracellular pathogens and other endosymbionts reprogram host cell transcription to suppress immune responses and recalibrate biosynthetic pathways. This reprogramming is critical in determining the outcome of infection or colonization. We combine pooled CRISPR knockout screening with dual host-microbe single-cell RNA sequencing, a method we term dual perturb-seq, to identify the molecular mediators of these transcriptional interactions. Applying dual perturb-seq to the intracellular pathogen Toxoplasma gondii, we are able to identify previously uncharacterized effector proteins and directly infer their function from the transcriptomic data. We show that TgGRA59 contributes to the export of other effector proteins from the parasite into the host cell and identify an effector, TgSOS1, that is necessary for sustained host STAT6 signaling and thereby contributes to parasite immune evasion and persistence. Together, this work demonstrates a tool that can be broadly adapted to interrogate host-microbe transcriptional interactions and reveal mechanisms of infection and immune evasion.

Identification of IMC43, a novel IMC protein that collaborates with IMC32 to form an essential daughter bud assembly complex in Toxoplasma gondii

The inner membrane complex (IMC) of Toxoplasma gondii is essential for all phases of the parasite's life cycle. One of its most critical roles is to act as a scaffold for the assembly of daughter buds during replication by endodyogeny. While many daughter IMC proteins have been identified, most are recruited after bud initiation and are not essential for parasite fitness. Here, we report the identification of IMC43, a novel daughter IMC protein that is recruited at the earliest stages of daughter bud initiation. Using an auxin-inducible degron system we show that depletion of IMC43 results in aberrant morphology, dysregulation of endodyogeny, and an extreme defect in replication. Deletion analyses reveal a region of IMC43 that plays a role in localization and a C-terminal domain that is essential for the protein's function. TurboID proximity labelling and a yeast two-hybrid screen using IMC43 as bait identify 30 candidate IMC43 binding partners. We investigate two of these: the essential daughter protein IMC32 and a novel daughter IMC protein we named IMC44. We show that IMC43 is responsible for regulating the localization of both IMC32 and IMC44 at specific stages of endodyogeny and that this regulation is dependent on the essential C-terminal domain of IMC43. Using pairwise yeast two-hybrid assays, we determine that this region is also sufficient for binding to both IMC32 and IMC44. As IMC43 and IMC32 are both essential proteins, this work reveals the existence of a bud assembly complex that forms the foundation of the daughter IMC during endodyogeny.

An apicomplexan bromodomain protein, TgBDP1, associates with diverse epigenetic factors to regulate essential transcriptional processes in Toxoplasma gondii

The protozoan pathogen Toxoplasma gondii relies on tight regulation of gene expression to invade and establish infection in its host. The divergent gene regulatory mechanisms of Toxoplasma and related apicomplexan pathogens rely heavily on regulators of chromatin structure and histone modifications. The important contribution of histone acetylation for Toxoplasma in both acute and chronic infection has been demonstrated, where histone acetylation increases at active gene loci. However, the direct consequences of specific histone acetylation marks and the chromatin pathway that influences transcriptional regulation in response to the modification are unclear. As a reader of lysine acetylation, the bromodomain serves as a mediator between the acetylated histone and transcriptional regulators. Here we show that the bromodomain protein, TgBDP1, which is conserved among Apicomplexa and within the Alveolata superphylum, is essential for Toxoplasma asexual proliferation. Using cleavage under targets and tagmentation, we demonstrate that TgBDP1 is recruited to transcriptional start sites of a large proportion of parasite genes. Transcriptional profiling during TgBDP1 knockdown revealed that loss of TgBDP1 leads to major dysregulation of gene expression, implying multiple roles for TgBDP1 in both gene activation and repression. This is supported by interactome analysis of TgBDP1 demonstrating that TgBDP1 forms a core complex with two other bromodomain proteins and an ApiAP2 factor. This core complex appears to interact with other epigenetic factors such as nucleosome remodeling complexes. We conclude that TgBDP1 interacts with diverse epigenetic regulators to exert opposing influences on gene expression in the Toxoplasma tachyzoite. IMPORTANCE Histone acetylation is critical for proper regulation of gene expression in the single-celled eukaryotic pathogen Toxoplasma gondii. Bromodomain proteins are "readers" of histone acetylation and may link the modified chromatin to transcription factors. Here, we show that the bromodomain protein TgBDP1 is essential for parasite survival and that loss of TgBDP1 results in global dysregulation of gene expression. TgBDP1 is recruited to the promoter region of a large proportion of parasite genes, forms a core complex with two other bromodomain proteins, and interacts with different transcriptional regulatory complexes. We conclude that TgBDP1 is a key factor for sensing specific histone modifications to influence multiple facets of transcriptional regulation in Toxoplasma gondii.

Functional screening reveals Toxoplasma prenylated proteins required for endocytic trafficking and rhoptry protein sorting

In the apicomplexans, endocytosed cargos (e.g., hemoglobin) are trafficked to a specialized organelle for digestion. This follows a unique endocytotic process at the micropore/cytostome in these parasites. However, the mechanism underlying endocytic trafficking remains elusive, due to the repurposing of classical endocytic proteins for the biogenesis of apical organelles. To resolve this issue, we have exploited the genetic tractability of the model apicomplexan Toxoplasma gondii, which ingests host cytosolic materials (e.g., green fluorescent protein[GFP]). We determined an association between protein prenylation and endocytic trafficking, and using an alkyne-labeled click chemistry approach, the prenylated proteome was characterized. Genome editing, using clustered regularly interspaced short palindromic repaet/CRISPR-associated nuclease 9 (CRISPR/Cas9), was efficiently utilized to generate genetically modified lines for the functional screening of 23 prenylated candidates. This identified four of these proteins that regulate the trafficking of endocytosed GFP vesicles. Among these proteins, Rab1B and YKT6.1 are highly conserved but are non-classical endocytic proteins in eukaryotes. Confocal imaging analysis showed that Rab1B and Ras are substantially localized to both the trans-Golgi network and the endosome-like compartments in the parasite. Conditional knockdown of Rab1B caused a rapid defect in secretory trafficking to the rhoptry bulb, suggesting a trafficking intersection role for the key regulator Rab1B. Further experiments confirmed a critical role for protein prenylation in regulating the stability/activity of these proteins (i.e., Rab1B and YKT6.1) in the parasite. Our findings define the molecular basis of endocytic trafficking and reveal a potential intersection function of Rab1B on membrane trafficking in T. gondii. This might extend to other related protists, including the malarial parasites. IMPORTANCE The protozoan Toxoplasma gondii establishes a permissive niche, in host cells, that allows parasites to acquire large molecules such as proteins. Numerous studies have demonstrated that the parasite repurposes the classical endocytic components for secretory sorting to the apical organelles, leaving the question of endocytic transport to the lysosome-like compartment unclear. Recent studies indicated that endocytic trafficking is likely to associate with protein prenylation in malarial parasites. This information promoted us to examine this association in the model apicomplexan T. gondii and to identify the key components of the prenylated proteome that are involved. By exploiting the genetic tractability of T. gondii and a host GFP acquisition assay, we reveal four non-classical endocytic proteins that regulate the transport of endocytosed cargos (e.g., GFP) in T. gondii. Thus, we extend the principle that protein prenylation regulates endocytic trafficking and elucidate the process of non-classical endocytosis in T. gondii and potentially in other related protists.

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

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

scDual-Seq of Toxoplasma gondii-infected mouse BMDCs reveals heterogeneity and differential infection dynamics

Dendritic cells and macrophages are integral parts of the innate immune system and gatekeepers against infection. The protozoan pathogen, Toxoplasma gondii, is known to hijack host immune cells and modulate their immune response, making it a compelling model to study host-pathogen interactions. Here we utilize single cell Dual RNA-seq to parse out heterogeneous transcription of mouse bone marrow-derived dendritic cells (BMDCs) infected with two distinct genotypes of T. gondii parasites, over multiple time points post infection. We show that the BMDCs elicit differential responses towards T. gondii infection and that the two parasite lineages distinctly manipulate subpopulations of infected BMDCs. Co-expression networks define host and parasite genes, with implications for modulation of host immunity. Integrative analysis validates previously established immune pathways and additionally, suggests novel candidate genes involved in host-pathogen interactions. Altogether, this study provides a comprehensive resource for characterizing host-pathogen interplay at high-resolution.

Toxoplasma ERK7 protects the apical complex from premature degradation

Accurate cellular replication balances the biogenesis and turnover of complex structures. In the apicomplexan parasite Toxoplasma gondii, daughter cells form within an intact mother cell, creating additional challenges to ensuring fidelity of division. The apical complex is critical to parasite infectivity and consists of apical secretory organelles and specialized cytoskeletal structures. We previously identified the kinase ERK7 as required for maturation of the apical complex in Toxoplasma. Here, we define the Toxoplasma ERK7 interactome, including a putative E3 ligase, CSAR1. Genetic disruption of CSAR1 fully suppresses loss of the apical complex upon ERK7 knockdown. Furthermore, we show that CSAR1 is normally responsible for turnover of maternal cytoskeleton during cytokinesis, and that its aberrant function is driven by mislocalization from the parasite residual body to the apical complex. These data identify a protein homeostasis pathway critical for Toxoplasma replication and fitness and suggest an unappreciated role for the parasite residual body in compartmentalizing processes that threaten the fidelity of parasite development.

Toxoplasma gondii encodes an array of TBC-domain containing proteins including an essential regulator that targets Rab2 in the secretory pathway

Toxoplasma gondii resides in its intracellular niche by employing a series of specialized secretory organelles that play roles in invasion, host-cell manipulation and parasite replication. Rab GTPases are major regulators of the parasite's secretory traffic that function as nucleotide dependent molecular switches to control vesicle trafficking. While many of the Rab proteins have been characterized in T. gondii , precisely how these Rabs are regulated remains poorly understood. To better understand the parasite's secretory traffic, we investigated the entire family of Tre2-Bub2-Cdc16 (TBC)-domain containing proteins, which are known to be involved in vesicle fusion and secretory protein trafficking. We first determined the localization of all 18 TBC-domain containing proteins to discrete regions of the secretory pathway or other vesicles in the parasite. We then use an auxin-inducible degron approach to demonstrate that the protozoan-specific TgTBC9 protein that localizes to the ER is essential for parasite survival. Knockdown of TgTBC9 results in parasite growth arrest and affects the organization of the ER and Golgi apparatus. We show that the conserved dual-finger active site in the TBC-domain of the protein is critical for its GTPase-activating protein (GAP) function and that the P. falciparum orthologue of TgTBC9 can rescue the lethal knockdown. We additionally show by immunoprecipitation and yeast two hybrid analyses that TgTBC9 directly binds Rab2, indicating that this TBC-Rab pair controls ER to Golgi traffic in the parasite. Together, these studies identify the first essential TBC protein described in any protozoan, provide new insight into intracellular vesicle trafficking in T. gondii , and reveal promising targets for the design of novel therapeutics that can specifically target apicomplexan parasites.

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.

Toxoplasma gondii HOOK-FTS-HIP Complex is Critical for Secretory Organelle Discharge during Motility, Invasion, and Egress

Members of the Apicomplexa phylum possess specialized secretory organelles that discharge, apically and in a timely regulated manner, key factors implicated in parasite motility, host cell invasion, egress and subversion of host cellular functions. The mechanisms regulating trafficking and apical docking of these secretory organelles are only partially elucidated. Here, we characterized two conserved endosomal trafficking regulators known to promote vesicle transport and/or fusion, HOOK and Fused Toes (FTS), in the context of organelle discharge in Toxoplasma gondii. TgHOOK and TgFTS form a complex with a coccidian-specific partner, named HOOK interacting partner (HIP). TgHOOK displays an apically enriched vesicular pattern and concentrates at the parasite apical tip where it colocalizes with TgFTS and TgHIP. Functional investigations revealed that TgHOOK is dispensable but fitness conferring. The protein regulates the apical positioning and secretion of micronemes and contributes to egress, motility, host cell attachment, and invasion. Conditional depletion of TgFTS or TgHIP impacted on the same processes but led to more severe phenotypes. This study provides evidence of endosomal trafficking regulators involved in the apical exocytosis of micronemes and possibly as a consequence or directly on the discharge of the rhoptries. IMPORTANCE Toxoplasma gondii affects between 30 and 80% of the human population, poses a life-threatening risk to immunocompromised individuals, and is a cause of abortion and birth defects following congenital transmission. T. gondii belongs to the phylum of Apicomplexa characterized by a set of unique apical secretory organelles called the micronemes and rhoptries. Upon host cell recognition, this obligatory intracellular parasite secretes specific effectors contained in micronemes and rhoptries to promote parasite invasion of host cells and subsequent persistence. Here, we identified novel T. gondii endosomal trafficking regulators and demonstrated that they regulate microneme organelle apical positioning and exocytosis, thereby strongly contributing to host cell invasion and parasite virulence.

IMC29 Plays an Important Role in Toxoplasma Endodyogeny and Reveals New Components of the Daughter-Enriched IMC Proteome

The Toxoplasma inner membrane complex (IMC) is a unique organelle that plays critical roles in parasite motility, invasion, egress, and replication. The IMC is delineated into the apical, body, and basal regions, defined by proteins that localize to these distinct subcompartments. The IMC can be further segregated by proteins that localize specifically to the maternal IMC, the daughter bud IMC, or both. While the function of the maternal IMC has been better characterized, the precise roles of most daughter IMC components remain poorly understood. Here, we demonstrate that the daughter protein IMC29 plays an important role in parasite replication. We show that Δimc29 parasites exhibit severe replication defects, resulting in substantial growth defects and loss of virulence. Deletion analyses revealed that IMC29 localization is largely dependent on the N-terminal half of the protein containing four predicted coiled-coil domains while IMC29 function requires a short C-terminal helical region. Using proximity labeling, we identify eight novel IMC proteins enriched in daughter buds, significantly expanding the daughter IMC proteome. We additionally report four novel proteins with unique localizations to the interface between two parasites or to the outer face of the IMC, exposing new subregions of the organelle. Together, this work establishes IMC29 as an important early daughter bud component of replication and uncovers an array of new IMC proteins that provides important insights into this organelle. IMPORTANCE The inner membrane complex (IMC) is a conserved structure across the Apicomplexa phylum, which includes obligate intracellular parasites that cause toxoplasmosis, malaria, and cryptosporidiosis. The IMC is critical for the parasite to maintain its intracellular lifestyle, particularly in providing a scaffold for daughter bud formation during parasite replication. While many IMC proteins in the later stages of division have been identified, components of the early stages of division remain unknown. Here, we focus on the early daughter protein IMC29, demonstrating that it is crucial for faithful parasite replication and identifying specific regions of the protein that are important for its localization and function. We additionally use proximity labeling to reveal a suite of daughter-enriched IMC proteins, which represent promising candidates to further explore this IMC subcompartment.

Babesia divergens egress from host cells is orchestrated by essential and druggable kinases and proteases

Apicomplexan egress from host cells is fundamental to the spread of infection and is poorly characterized in Babesia spp., parasites of veterinary importance and emerging zoonoses. Through the use of video microscopy, transcriptomics and chemical genetics, we have implicated signaling, proteases and gliding motility as key drivers of egress by Babesia divergens. We developed reverse genetics to perform a knockdown screen of putative mediators of egress, identifying kinases and proteases involved in distinct steps of egress (ASP3, PKG and CDPK4) and invasion (ASP2, ASP3 and PKG). Inhibition of egress leads to continued intracellular replication, indicating exit from the replication cycle is uncoupled from egress. Chemical genetics validated PKG, ASP2 and ASP3 as druggable targets in Babesia spp. All taken together, egress in B. divergens more closely resembles T. gondii than the more evolutionarily-related Plasmodium spp. We have established a molecular framework for biological and translational studies of B. divergens egress.

Functional Characterization of 15 Novel Dense Granule Proteins in Toxoplasma gondii Using the CRISPR-Cas9 System

The analysis of the subcellular localization and function of dense granule proteins (GRAs) is of central importance for the understanding of host-parasite interaction and pathogenesis of Toxoplasma gondii infection. Here, we identified 15 novel GRAs and used C-terminal endogenous gene tagging to determine their localization at the intravacuolar network (IVN), parasitophorous vacuole (PV), or PV membrane (PVM) in the tachyzoites and at the periphery of the bradyzoites-containing cysts. The functions of the 15 gra genes were examined in type I RH strain and 5 of these gra genes were also evaluated in the cyst-forming type II Pru strain. The 15 novel gra genes were successfully disrupted by using CRISPR-Cas9 mediated homologous recombination and the results showed that 13 gra genes were not individually essential for T. gondii replication in vitro or virulence in mice during acute and chronic infection. Intriguingly, deletion of TGME49_266410 and TGME49_315910 in both RH and Pru strains decreased the parasite replication in vitro and attenuated its virulence, and also reduced the cyst-forming ability of the Pru strain in mice during chronic infection. Comparison of the transcriptomic profiles of the 15 gra genes suggests that they may play roles in other life cycle stages and genotypes of T. gondii. Taken together, our findings improve the understanding of T. gondii pathogenesis and demonstrate the involvement of two novel GRAs, TGME49_266410 and TGME49_315910, in the parasite replication and virulence. IMPORTANCE Dense granule proteins (GRAs) play important roles in Toxoplasma gondii pathogenicity. However, the functions of many putative GRAs have not been elucidated. Here, we found that 15 novel GRAs are secreted into intravacuolar network (IVN), parasitophorous vacuole (PV), or PV membrane (PVM) in tachyzoites and are located at the periphery of the bradyzoite-containing cysts. TGME49_266410 and TGME49_315910 were crucial to the growth of RH and Pru strains in vitro. Deletion of TGME49_266410 and TGME49_315910 attenuated the parasite virulence in mice. However, disruption of other 13 gra genes did not have a significant impact on the proliferation and pathogenicity of T. gondii in vitro or in vivo. The marked effects of the two novel GRAs (TGME49_266410 and TGME49_315910) on the in vitro growth and virulence of T. gondii are notable and warrant further elucidation of the temporal and spatial dynamics of translocation of these two novel GRAs and how do they interfere with host cell functions.

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

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

TgTKL4 Is a Novel Kinase That Plays an Important Role in Toxoplasma Morphology and Fitness

Protein kinases of the protozoan parasite Toxoplasma gondii have been shown to play key roles in regulating parasite motility, invasion, replication, egress and survival within the host. The tyrosine kinase-like (TKL) kinase family of proteins are a set of poorly studied kinases that our recent studies have indicated play a critical role in Toxoplasma biology. In this study, we focused on TgTKL4, another member of the TKL family that is predicted to confer parasite fitness. Endogenous tagging of TgTKL4 identified it as a temporally oscillating kinase with dynamic localization in the parasite. Gene disruption experiments suggested that TgTKL4 is important for Toxoplasma propagation in vitro, and loss of this kinase resulted in replication and invasion defects. During parasite division, TgTKL4 expression was limited to the synthesis and mitosis-cytokinesis phases and, interestingly, loss of TgTKL4 led to defects in Toxoplasma morphology. Further analysis of the parasite cytoskeleton indicated that the subpellicular microtubules are shorter and more widely spaced in parasites lacking TgTKL4. Although loss of TgTKL4 caused only moderate changes in the gene expression profile, TgTKL4 null mutants exhibited significant changes in their global phospho-proteome, including in proteins that constitute the parasite cytoskeleton. Additionally, mice inoculated intraperitoneally with TgTKL4 knockout parasites showed increased survival rates, suggesting that TgTKL4 plays an important role in acute toxoplasmosis. Together, these findings suggest that TgTKL4 mediates a signaling pathway that regulates parasite morphology and is an important factor required for parasite fitness in vitro and in vivo. IMPORTANCE Toxoplasma gondii is a protozoan parasite that can cause life-threatening disease in mammals; hence, identifying key factors required for parasite growth and pathogenesis is important to develop novel therapeutics. In this study, we identified and characterized another member of the newly described TKL family, TgTKL4, a cell cycle-regulated kinase. By disrupting TgTKL4, we determined that this kinase is required for normal parasite growth in vitro and that loss of this kinase results in parasites with reduced competence in replication and invasion processes. Specifically, Toxoplasma parasites lacking TgTKL4 had defects in cytoskeletal arrangement, resulting in parasites with abnormal morphology. Phospho-proteome studies provided further clues that decreased phosphorylation of proteins that constitute the Toxoplasma cytoskeleton could be responsible for altered morphology in TgTKL4-deficient parasites. Additionally, loss of TgTKL4 resulted in attenuation of virulence in the animal model, suggesting that TgTKL4 is an important virulence factor. Hence, this study provides a novel insight into the importance of a TgTKL4 as a fitness-determining factor for Toxoplasma propagation in vitro and pathogenesis in vivo.

Effect of deleting four Toxoplasma gondii calcium-binding EGF domain-containing proteins on parasite replication and virulence

Several calcium-binding proteins including calcium-dependent protein kinases play important roles in several facets of the intracellular infection cycle of the apicomplexan protozoan parasite Toxoplasma gondii. However, the role of the calcium-binding epidermal growth factor (EGF) domain-containing proteins (CBDPs) remains poorly understood. In this study, we examined the functions of four CBDP genes in T. gondii RH strain of type I by generating knock-out strains using CRISPR-Cas9 system. We investigated the ability of mutant strains deficient in CBDP1, CBDP2, CBDP3, or CBDP4 to form plaques, replicate intracellularly, and egress from the host cells. The results showed that no definite differences between any of these four CBDP mutant strains and the wild-type strain in terms of their ability to form plaques, intracellular replication, and egress. Additionally, CBDP mutants did not exhibit any significant attenuated virulence compared to the wild-type strain in mice. The expression profiles of CBDP2-4 genes were conserved among T. gondii strains of different genotypes, life cycle stages, and developmental forms. Whether other CBDP genes play any roles in the pathogenicity of T. gondii strains of different genotypes remains to be elucidated.

Toxoplasma Shelph, a Phosphatase Located in the Parasite Endoplasmic Reticulum, Is Required for Parasite Virulence

Toxoplasma gondii is a single-celled parasitic eukaryote that evolved to successfully propagate in any nucleated cell. As with any other eukaryote, its life cycle is regulated by signaling pathways controlled by kinases and phosphatases. T. gondii encodes an atypical bacterial-like phosphatase absent from mammalian genomes, named Shelph, after its first identification in the psychrophilic bacterium Schewanella sp. Here, we demonstrate that Toxoplasma Shelph is an active phosphatase localized in the parasite endoplasmic reticulum. The phenotyping of a shelph knockout (KO) line showed a minor impairment in invasion on human fibroblasts, while the other steps of the parasite lytic cycle were not affected. In contrast with Plasmodium ortholog Shelph1, this invasion deficiency was not correlated with any default in the biogenesis of secretory organelles. However, Shelph-KO parasites displayed a much-pronounced defect in virulence in vivo. These phenotypes could be rescued by genetic complementation, thus supporting an important function for Shelph in the context of a natural infection. IMPORTANCE Toxoplasma gondii belongs to the Apicomplexa phylum, which comprises more than 5,000 species, among which is Plasmodium falciparum, the notorious agent of human malaria. Intriguingly, the Apicomplexa genomes encode at least one phosphatase closely related to the bacterial Schewanella phosphatase, or Shelph. To better understand the importance of these atypical bacterial enzymes in eukaryotic parasites, we undertook the functional characterization of T. gondii Shelph. Our results uncovered its subcellular localization and its enzymatic activity, revealed its subtle involvement during the tachyzoite invasion step of the lytic cycle, and more importantly, highlighted a critical requirement of this phosphatase for parasite propagation in mice. Overall, this study revealed an unexpected role for T. gondii Shelph in the maintenance of parasite virulence in vivo.

IMC10 and LMF1 mediate mitochondrial morphology through mitochondrion-pellicle contact sites in Toxoplasma gondii

The single mitochondrion of Toxoplasma gondii is highly dynamic, being predominantly in a peripherally distributed lasso-shape in intracellular parasites and collapsed in extracellular parasites. The peripheral positioning of the mitochondrion is associated with apparent contacts between the mitochondrion membrane and the parasite pellicle. The outer mitochondrial membrane-associated protein LMF1 is critical for the correct positioning of the mitochondrion. Intracellular parasites lacking LMF1 fail to form the lasso-shaped mitochondrion. To identify other proteins that tether the mitochondrion of the parasite to the pellicle, we performed a yeast two-hybrid screen for LMF1 interactors. We identified 70 putative interactors localized in different cellular compartments, such as the apical end of the parasite, mitochondrial membrane and the inner membrane complex (IMC), including with the pellicle protein IMC10. Using protein-protein interaction assays, we confirmed the interaction of LMF1 with IMC10. Conditional knockdown of IMC10 does not affect parasite viability but severely affects mitochondrial morphology in intracellular parasites and mitochondrial distribution to the daughter cells during division. In effect, IMC10 knockdown phenocopies disruption of LMF1, suggesting that these two proteins define a novel membrane tether between the mitochondrion and the IMC in Toxoplasma. This article has an associated First Person interview with the first author of the paper.

Identification of new components of the basal pole of Toxoplasma gondii provides novel insights into its molecular organization and functions

The Toxoplasma gondii tachyzoite is a singled-cell obligate intracellular parasite responsible for the acute phase of toxoplasmosis. This polarized cell exhibits an apical complex, a hallmark of the phylum Apicomplexa, essential for motility, invasion, and egress from the host cell. Located on the opposite end of the cell is the basal complex, an elaborated cytoskeletal structure that also plays critical roles in the lytic cycle of the parasite, being involved in motility, cell division, constriction and cytokinesis, as well as intravacuolar cell-cell communication. Nevertheless, only a few proteins of this structure have been described and functionally assessed. In this study, we used spatial proteomics to identify new basal complex components (BCC), and in situ imaging, including ultrastructure expansion microscopy, to position them. We thus confirmed the localization of nine BCCs out of the 12 selected candidates and assigned them to different sub-compartments of the basal complex, including two new domains located above the basal ring and below the posterior cup. Their functional investigation revealed that none of these BCCs are essential for parasite growth in vitro. However, one BCC is critical for constricting of the basal complex, likely through direct interaction with the class VI myosin heavy chain J (MyoJ), and for gliding motility. Four other BCCs, including a phosphatase and a guanylate-binding protein, are involved in the formation and/or maintenance of the intravacuolar parasite connection, which is required for the rosette organization and synchronicity of cell division.

Not your Mother's MAPKs: Apicomplexan MAPK function in daughter cell budding

Reversible phosphorylation by protein kinases is one of the core mechanisms by which biological signals are propagated and processed. Mitogen-activated protein kinases, or MAPKs, are conserved throughout eukaryotes where they regulate cell cycle, development, and stress response. Here, we review advances in our understanding of the function and biochemistry of MAPK signaling in apicomplexan parasites. As expected for well-conserved signaling modules, MAPKs have been found to have multiple essential roles regulating both Toxoplasma tachyzoite replication and sexual differentiation in Plasmodium. However, apicomplexan MAPK signaling is notable for the lack of the canonical kinase cascade that normally regulates the networks, and therefore must be regulated by a distinct mechanism. We highlight what few regulatory relationships have been established to date, and discuss the challenges to the field in elucidating the complete MAPK signaling networks in these parasites.

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

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

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.

The first apicoplast tRNA thiouridylase plays a vital role in the growth of Toxoplasma gondii

Toxoplasmosis caused by the protozoan Toxoplasma gondii is one of the most common parasitic diseases in humans and almost all warm-blooded animals. Lys, Glu, and Gln-specific tRNAs contain a super-modified 2-thiourea (s²U) derivatives at the position 34, which is essential for all living organisms by maintaining the structural stability and aminoacylation of tRNA, and the precision and efficiency of codon recognition during protein translation. However, the enzyme(s) involved in this modification in T. gondii remains elusive. In this report, three putative tRNA-specific 2-thiolation enzymes were identified, of which two were involved in the s²U34 modification of tRNA^Lys, tRNA^Glu, and tRNA^Gln. One was named TgMnmA, an apicoplast-located tRNA-specific 2-thiolation enzyme in T. gondii. Knockout of TgMnmA showed that this enzyme is important for the lytic cycle of tachyzoites. Loss of TgMnmA also led to abnormities in apicoplast biogenesis and severely disturbed apicoplast genomic transcription. Notably, mice survived from the infection with 10 TgMnmA-KO RH tachyzoites. These findings provide new insights into s²U34 tRNA modification in Apicomplexa, and suggest TgMnmA, the first apicoplast tRNA thiouridylase identified in all apicomplexans, as a potential drug target.

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.

The BCC7 Protein Contributes to the Toxoplasma Basal Pole by Interfacing between the MyoC Motor and the IMC Membrane Network

T. gondii is a eukaryotic parasite that has evolved a stage called tachyzoite which multiplies in host cells by producing two daughter cells internally. These nascent tachyzoites bud off their mother and repeat the division process until the expanding progenies escape to settle and multiply in other host cells. Over these intra- and extra-cellular phases, the tachyzoite maintains an essential apicobasal polarity that emerges through a unique bidirectional budding process of the elongating cells. This process requires the assembly of several molecular complexes that, at the nascent pole, encompass structural and myosin motor elements. To characterize a recently identified basal pole marker named BCC7 with respect to the posterior myosin J and myosin C motors, we used conventional biochemistry as well as advanced proteomic and in silico analysis in conjunction with live and super resolution microscopy of transgenic fluorescent tachyzoites. We document that BCC7 forms a ribbed ring below which myosin C motor entities distribute regularly. In addition, we identified—among 13 BCC7 putative partners—two novel and five known members of the inner membrane complex (IMC) family which ends at the apical side of the ring. Therefore, BCC7 could assist the stabilization of the IMC plaques and contribute to the parasite biomechanical properties.

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.

TSS-seq of Toxoplasma gondii sporozoites revealed a novel motif in stage-specific promoters

Toxoplasma gondii is one of the most common zoonotic protozoan parasites. It has three major infectious stages: rapidly multiplying tachyzoites (Tz), slowly replicating bradyzoites (Bz) and a resting/free-living stage, sporozoites (Sz). The regulatory mechanisms governing stage-specific gene expression are not fully understood. Few transcriptional start sites (TSS) are known for Sz. In this study, we obtained TSS of Sz using an oligo-capping method and RNA-seq analysis. We identified 1,043,503 TSS in the Sz transcriptome. These defined 38,973 TSS clusters, of which, 11,925 were expressed in Sz and 1535 TSS differentially expressed in Sz. Based on these data, we defined promoter regions and novel sporozoite stage-specific motifs using MEME. TGTANNTACA was distributed around -55 to -75 regions from each TSS. Interestingly, the same motif was reported in another apicomplexan, Plasmodium berghei, as a cis-element of female-specific gametocyte genes, implying the presence of common regulatory machinery. Further comparative analysis should better define the distribution and function of these elements in other members of this important parasitic phylum.

Transcriptional signatures of clonally derived Toxoplasma tachyzoites reveal novel insights into the expression of a family of surface proteins

Toxoplasma gondii has numerous, large, paralogous gene families that are likely critical for supporting its unparalleled host range: nearly any nucleated cell in almost any warm-blooded animal. The SRS (SAG1-related sequence) gene family encodes over 100 proteins, the most abundant of which are thought to be involved in parasite attachment and, based on their stage-specific expression, evading the host immune response. For most SRS proteins, however, little is understood about their function and expression profile. Single-parasite RNA-sequencing previously demonstrated that across an entire population of lab-grown tachyzoites, transcripts for over 70 SRS genes were detected in at least one parasite. In any one parasite, however, transcripts for an average of only 7 SRS genes were detected, two of which, SAG1 and SAG2A, were extremely abundant and detected in virtually all. These data do not address whether this pattern of sporadic SRS gene expression is consistently inherited among the progeny of a given parasite or arises independently of lineage. We hypothesized that if SRS expression signatures are stably inherited by progeny, subclones isolated from a cloned parent would be more alike in their expression signatures than they are to the offspring of another clone. In this report, we compare transcriptomes of clonally derived parasites to determine the degree to which expression of the SRS family is stably inherited in individual parasites. Our data indicate that in RH tachyzoites, SRS genes are variably expressed even between parasite samples subcloned from the same parent within approximately 10 parasite divisions (72 hours). This suggests that the pattern of sporadically expressed SRS genes is highly variable and not driven by inheritance mechanisms, at least under our conditions.

Genome-wide localization of histone variants in Toxoplasma gondii implicates variant exchange in stage-specific gene expression

BACKGROUND

Toxoplasma gondii is a protozoan parasite that differentiates from acute tachyzoite stages to latent bradyzoite forms in response to environmental cues that modify the epigenome. We studied the distribution of the histone variants CenH3, H3.3, H2A.X, H2A.Z and H2B.Z, by genome-wide chromatin immunoprecipitation to understand the role of variant histones in developmental transitions of T. gondii parasites.

RESULTS

H3.3 and H2A.X were detected in telomere and telomere associated sequences, whereas H3.3, H2A.X and CenH3 were enriched in centromeres. Histones H2A.Z and H2B.Z colocalize with the transcriptional activation mark H3K4me3 in promoter regions surrounding the nucleosome-free region upstream of the transcription start site. The H2B.Z/H2A.Z histone pair also localizes to the gene bodies of genes that are silent but poised for activation, including bradyzoite stage-specific genes. The majority of H2A.X and H2A.Z/H2B.Z loci do not overlap, consistent with variant histones demarcating specific functional regions of chromatin. The extent of enrichment of H2A.Z/H2B.Z (and H3.3 and H2A.X) within the entire gene (5'UTR and gene body) reflects the timing of gene expression during the cell cycle, suggesting that dynamic turnover of H2B.Z/H2A.Z occurs during the tachyzoite cell cycle. Thus, the distribution of the variant histone H2A.Z/H2B.Z dimer defines active and developmentally silenced regions of the T. gondii epigenome including genes that are poised for expression.

CONCLUSIONS

Histone variants mark functional regions of parasite genomes with the dynamic placement of the H2A.Z/H2B.Z dimer implicated as an evolutionarily conserved regulator of parasite and eukaryotic differentiation.

Functional Characterization of 17 Protein Serine/Threonine Phosphatases in Toxoplasma gondii Using CRISPR-Cas9 System

Protein serine/threonine phosphatases (PSPs), found in various plants and protozoa, are involved in the regulation of various biological processes. However, very little is known about the role of PSPs in the pathogenicity of the apicomplexan protozoan Toxoplasma gondii. Herein, the subcellular localization of 17 PSPs (PP5, PP7, EFPP, SLP, PPM3F, PPM4, PPM5A, PPM5B, PPM6, PPM8, PPM9, PPM12, PPM14, PPM18, CTD1, CTD2, and CTD3) was examined by 6× HA tagging of endogenous genes in C-terminal. The PSPs were detected in the cytoplasm (PP5, EFPP, PPM8, and CTD2), dense granules (SLP), nucleus (PPM4 and PPM9), inner membrane complex (PPM12), basal complex (CTD3), and apical pole (PP7). The remaining PSPs exhibited low or undetectable level of expression. To characterize the contribution of these genes to the infectivity of T. gondii, knock-out (KO) strains of type I RH strain deficient in the 17 psp genes and KO type II Pru strain deficient in pp7 and slp genes were constructed. The pathogenicity of individual RHΔpsp mutants was characterized in vitro using plaque, egress, and intracellular replication assays, and mouse infection, while pathogenicity of PruΔpp7 and PruΔslp mutant strains was evaluated by examining the parasite lytic cycle in vitro and assessment of brain cyst burden in mice. No significant differences were observed between 16 RHΔpsp strains and wild-type (WT) RH strain. However, RHΔpp7 exhibited significantly lower invasion efficiency and parasitophorous vacuole formation in vitro, and less virulence in mice compared with other RHΔpsp and WT strains. In addition, PruΔpp7 exhibited marked attenuation of virulence and significant reduction in the brain cyst burden in mice compared with PruΔslp and WT strains, suggesting the key role of PP7 in the virulence of T. gondii. Comparative transcriptomic profiling of the 17 psp genes showed that they may play different roles in the pathogenesis of different genotypes or life cycle stages of T. gondii. These findings provide new insight into the role of PSPs in the pathogenesis of T. gondii.

The Extracellular Milieu of Toxoplasma's Lytic Cycle Drives Lab Adaptation, Primarily by Transcriptional Reprogramming

Evolve and resequencing (E&R) was applied to lab adaptation of Toxoplasma gondii for over 1,500 generations with the goal of mapping host-independent in vitro virulence traits. Phenotypic assessments of steps across the lytic cycle revealed that only traits needed in the extracellular milieu evolved. Nonsynonymous single-nucleotide polymorphisms (SNPs) in only one gene, a P4 flippase, fixated across two different evolving populations, whereas dramatic changes in the transcriptional signature of extracellular parasites were identified. Newly developed computational tools correlated phenotypes evolving at different rates with specific transcriptomic changes. A set of 300 phenotype-associated genes was mapped, of which nearly 50% is annotated as hypothetical. Validation of a select number of genes by knockouts confirmed their role in lab adaptation and highlights novel mechanisms underlying in vitro virulence traits. Further analyses of differentially expressed genes revealed the development of a “pro-tachyzoite” profile as well as the upregulation of the fatty acid biosynthesis (FASII) pathway. The latter aligned with the P4 flippase SNP and aligned with a low abundance of medium-chain fatty acids at low passage, indicating this is a limiting factor in extracellular parasites. In addition, partial overlap with the bradyzoite differentiation transcriptome in extracellular parasites indicated that stress pathways are involved in both situations. This was reflected in the partial overlap between the assembled ApiAP2 and Myb transcription factor network underlying the adapting extracellular state with the bradyzoite differentiation program. Overall, E&R is a new genomic tool successfully applied to map the development of polygenic traits underlying in vitro virulence of T. gondii.

IMPORTANCE It has been well established that prolonged in vitro cultivation of Toxoplasma gondii augments progression of the lytic cycle. This lab adaptation results in increased capacities to divide, migrate, and survive outside a host cell, all of which are considered host-independent virulence factors. However, the mechanistic basis underlying these enhanced virulence features is unknown. Here, E&R was utilized to empirically characterize the phenotypic, genomic, and transcriptomic changes in the non-lab-adapted strain, GT1, during 2.5 years of lab adaptation. This identified the shutdown of stage differentiation and upregulation of lipid biosynthetic pathways as the key processes being modulated. Furthermore, lab adaptation was primarily driven by transcriptional reprogramming, which rejected the starting hypothesis that genetic mutations would drive lab adaptation. Overall, the work empirically shows that lab adaptation augments T. gondii’s in vitro virulence by transcriptional reprogramming and that E&R is a powerful new tool to map multigenic traits.

Revisiting the Role of Toxoplasma gondii ERK7 in the Maintenance and Stability of the Apical Complex

Toxoplasma gondii extracellular signal-regulated kinase 7 (ERK7) is known to contribute to the integrity of the apical complex and to participate in the final step of conoid biogenesis. In the absence of ERK7, mature parasites lose their conoid complex and are unable to glide, invade, or egress from host cells. In contrast to a previous report, we show here that the depletion of ERK7 phenocopies the depletion of the apical cap protein AC9 or AC10. The absence of ERK7 leads to the loss of the apical polar ring (APR), the disorganization of the basket of subpellicular microtubules (SPMTs), and a severe impairment in microneme secretion. Ultrastructure expansion microscopy (U-ExM), coupled to N-hydroxysuccinimide ester (NHS-ester) staining on intracellular parasites, offers an unprecedented level of resolution and highlights the disorganization of the rhoptries as well as the dilated plasma membrane at the apical pole in the absence of ERK7. Comparative proteomics analysis of wild-type and ERK7-depleted parasites confirmed the disappearance of known apical complex proteins, including markers of the apical polar ring and a new apical cap named AC11. Concomitantly, the absence of ERK7 led to an accumulation of microneme proteins, resulting from the defect in the exocytosis of the organelles. AC9-depleted parasites were included as controls and exhibited an increase in inner membrane complex proteins, with two new proteins assigned to this compartment, namely, IMC33 and IMC34. IMPORTANCE The conoid is an enigmatic, dynamic organelle positioned at the apical tip of the coccidian subgroup of the Apicomplexa, close to the apical polar ring (APR) from which the subpellicular microtubules (SPMTs) emerge and through which the secretory organelles (micronemes and rhoptries) reach the plasma membrane for exocytosis. In Toxoplasma gondii, the conoid protrudes concomitantly with microneme secretion, during egress, motility, and invasion. The conditional depletion of the apical cap structural protein AC9 or AC10 leads to a disorganization of SPMTs as well as the loss of the APR and conoid, resulting in a microneme secretion defect and a block in motility, invasion, and egress. We show here that the depletion of the kinase ERK7 phenocopies AC9 and AC10 mutants. The combination of ultrastructure expansion microscopy and NHS-ester staining revealed that ERK7-depleted parasites exhibit a dilated apical plasma membrane and an altered positioning of the rhoptries, while electron microscopy images unambiguously highlight the loss of the APR.

A Novel Toxoplasma Inner Membrane Complex Suture-Associated Protein Regulates Suture Protein Targeting and Colocalizes with Membrane Trafficking Machinery

The cytoskeleton of Toxoplasma gondii is composed of the inner membrane complex (IMC) and an array of underlying microtubules that provide support at the periphery of the parasite. Specific subregions of the IMC carry out distinct roles in replication, motility, and host cell invasion. Building on our previous in vivo biotinylation (BioID) experiments of the IMC, we identified here a novel protein that localizes to discrete puncta that are embedded in the parasite’s cytoskeleton along the IMC sutures. Gene knockout analysis showed that loss of the protein results in defects in cytoskeletal suture protein targeting, cytoskeletal integrity, parasite morphology, and host cell invasion. We then used deletion analyses to identify a domain in the N terminus of the protein that is critical for both localization and function. Finally, we used the protein as bait for in vivo biotinylation, which identified several other proteins that colocalize in similar spot-like patterns. These putative interactors include several proteins that are implicated in membrane trafficking and are also associated with the cytoskeleton. Together, these data reveal an unexpected link between the IMC sutures and membrane trafficking elements of the parasite and suggest that the suture puncta are likely a portal for trafficking cargo across the IMC.

The Role of Type II Fatty Acid Synthesis Enzymes FabZ, ODSCI, and ODSCII in the Pathogenesis of Toxoplasma gondii Infection

Toxoplasma gondii is an obligate intracellular protozoan parasite, which has a worldwide distribution and can infect a large number of warm-blooded animals and humans. T. gondii must colonize and proliferate inside the host cells in order to maintain its own survival by securing essential nutrients for the development of the newly generated tachyzoites. The type II fatty acid biosynthesis pathway (FASII) in the apicoplast is essential for the growth and survival of T. gondii. We investigated whether deletion of genes in the FASII pathway influences the in vitro growth and in vivo virulence of T. gondii. We focused on beta-hydroxyacyl-acyl carrier protein dehydratase (FabZ) and oxidoreductase, short chain dehydrogenase/reductase family proteins ODSCI and ODSCII. We constructed T. gondii strains deficient in FabZ, ODSCI, and ODSCII using CRISPR-Cas9 gene editing technology. The results of immunofluorescence assay, plaque assay, proliferation assay and egress assay showed that in RHΔFabZ strain the apicoplast was partly lost and the growth ability of the parasite in vitro was significantly inhibited, while for RHΔODSCI and RHΔODSCII mutant strains no similar changes were detected. RHΔFabZ exhibited reduced virulence for mice compared with RHΔODSCI and RHΔODSCII, as shown by the improved survival rate. Deletion of FabZ in the PRU strain significantly decreased the brain cyst burden in mice compared with PRUΔODSCI and PRUΔODSCII. Collectively, these findings suggest that FabZ contributes to the growth and virulence of T. gondii, while ODSCI and ODSCII do not contribute to these traits.

Cryptosporidium rhoptry effector protein ROP1 injected during invasion targets the host cytoskeletal modulator LMO7

The parasite Cryptosporidium invades and replicates in intestinal epithelial cells and is a leading cause of diarrheal disease and early childhood mortality. The molecular mechanisms that underlie infection and pathogenesis are largely unknown. Here, we delineate the events of host cell invasion and uncover a mechanism unique to Cryptosporidium. We developed a screen to identify parasite effectors, finding the injection of multiple parasite proteins into the host from the rhoptry organelle. These factors are targeted to diverse locations within the host cell and its interface with the parasite. One identified effector, rhoptry protein 1 (ROP1), accumulates in the terminal web of enterocytes through direct interaction with the host protein LIM domain only 7 (LMO7) an organizer of epithelial cell polarity and cell-cell adhesion. Genetic ablation of LMO7 or ROP1 in mice or parasites, respectively, impacts parasite burden in vivo in opposite ways. Taken together, these data provide molecular insight into how Cryptosporidium manipulates its intestinal host niche.

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.