- Research article
- Open Access
Asexual expansion of Toxoplasma gondii merozoites is distinct from tachyzoites and entails expression of non-overlapping gene families to attach, invade, and replicate within feline enterocytes
BMC Genomics volume 16, Article number: 66 (2015)
The apicomplexan parasite Toxoplasma gondii is cosmopolitan in nature, largely as a result of its highly flexible life cycle. Felids are its only definitive hosts and a wide range of mammals and birds serve as intermediate hosts. The latent bradyzoite stage is orally infectious in all warm-blooded vertebrates and establishes chronic, transmissible infections. When bradyzoites are ingested by felids, they transform into merozoites in enterocytes and expand asexually as part of their coccidian life cycle. In all other intermediate hosts, however, bradyzoites differentiate exclusively to tachyzoites, and disseminate extraintestinally to many cell types. Both merozoites and tachyzoites undergo rapid asexual population expansion, yet possess different effector fates with respect to the cells and tissues they develop in and the subsequent stages they differentiate into.
To determine whether merozoites utilize distinct suites of genes to attach, invade, and replicate within feline enterocytes, we performed comparative transcriptional profiling on purified tachyzoites and merozoites. We used high-throughput RNA-Seq to compare the merozoite and tachyzoite transcriptomes. 8323 genes were annotated with sequence reads across the two asexually replicating stages of the parasite life cycle. Metabolism was similar between the two replicating stages. However, significant stage-specific expression differences were measured, with 312 transcripts exclusive to merozoites versus 453 exclusive to tachyzoites. Genes coding for 177 predicted secreted proteins and 64 membrane- associated proteins were annotated as merozoite-specific. The vast majority of known dense-granule (GRA), microneme (MIC), and rhoptry (ROP) genes were not expressed in merozoites. In contrast, a large set of surface proteins (SRS) was expressed exclusively in merozoites.
The distinct expression profiles of merozoites and tachyzoites reveal significant additional complexity within the T. gondii life cycle, demonstrating that merozoites are distinct asexual dividing stages which are uniquely adapted to their niche and biological purpose.
Toxoplasma gondii is an intracellular zoonotic parasite that chronically infects 30% of the world’s human population . It has a complex life cycle, infecting a wide range of mammals and birds as intermediate hosts but with felids as the only definitive hosts. Intermediate hosts can become infected through ingestion of oocysts shed into the environment via cat feces or by ingesting tissue cysts in meat or viscera . During acute infection of naïve intermediate hosts, sporozoites from oocysts or bradyzoites from tissue cysts differentiate quickly into tachyzoites, which divide rapidly and disseminate throughout the host’s body. Tachyzoites are, however, controlled efficiently by cell mediated immunity and the parasite reverts to the slowly replicating or quiescent bradyzoite form, which resides within tissue cysts that are particularly abundant in brain and heart muscle but are also present throughout skeletal muscle . When cats ingest tissue cysts, bradyzoites can take a different developmental pathway. They are released from their cysts and invade enterocytes of the small intestine, transforming into schizonts [4,5]. The parasite population that develops in cat enterocytes undergoes a classical coccidian cycle involving several rounds of asexual division and amplification followed by differentiation into macro- and microgamonts, the dimorphic stages of sexual development. Microgametes fertilize macrogametes, producing diploid zygotes that subsequently develop into unsporulated oocysts that are excreted in the feces of the cat. The sexual phase continues within the oocyst as meiosis ensues, followed by mitosis to produce infectious sporozoites, encased within sporocysts inside the oocysts.
The readily culturable, rapidly dividing tachyzoite is the best studied form of T. gondii by far – there is abundant information about cell cycle, metabolism and host parasite interactions for this stage . In contrast, the merozoite, which is the other rapidly dividing asexual form of T. gondii that ultimately generates hundreds of millions of gametes, is the least well studied developmental stage. This is largely because merozoites are not cultivatable in vitro and difficult to access in vivo. Thus, investigation of the molecular mechanisms governing initiation of parasite amplification preceding the development of sexual stages in cats has been severely hampered. The need for a more detailed understanding of the development in the definitive host is underscored by the fact that infected cats shed hundreds of millions of oocysts that can remain infectious for over a year to a wide range of highly susceptible intermediate hosts, including humans . Here, we developed improved protocols for cat infection, parasite isolation, and next generation sequencing to close this knowledge gap by building a transcription profile for the merozoite stage of enteroepithelial development. Using genome-wide comparative transcriptomics, we show that merozoites express distinct gene families in a stage-specific fashion, and fail to express the majority of annotated ROP, GRA and MIC proteins which are upregulated during tachyzoite replication. Among the most highly differentially regulated parasite proteins were several large gene families, including those coding for SRS proteins found on the parasite cell surface. Other key genes expressed by tachyzoites whose products are known to facilitate motility, host seeking, attachment, invasion, and remodeling of the parasitophorous vacuole (PV) within the parasitized host cell were not expressed in merozoites. This strongly suggests that merozoites are biologically distinct and utilize a different suite of genes that are necessary for asexual expansion within feline enterocytes prior to gamont development.
Results and discussion
Purification of merozoites from feline enterocytes for RNA-Seq analysis
T. gondii parasite preparations were generated from enterocytes from an infected cat at onset of patency (5d post infection). The enterocyte cell layer containing replicating parasites as shown by IFA (Figure 1A) from two regions of the rinsed and opened small intestine were selectively harvested by mechanical stripping (Figure 1B). Repeated centrifugation and re-suspension of this material with ice-cold 0.05% Tween80 in PBS yielded microscopically pure fractions of extracellular merozoites (two biological replicates). The released merozoites showed differential staining with Diff-Quick (Figure 1C); no gametocyte stages were observed in this fraction by microscopy. However, even though there is no evidence for this in the RNA-Seq data (see also in Methods) we cannot completely exclude minor contributions from sexual stages (macrogametes, microgametes) to the RNA pool selected for analysis. Nevertheless, based on the overwhelming majority of merozoites in the sample at this early time point of infection we will henceforth refer to this fraction as “merozoites”. Although further Percoll gradient purification yielded merozoites that were completely free of host material, we used the detergent-washed preparations containing minor contamination with host RNA (Figure 1D) for sequencing to avoid any changes to the parasite transcriptome due to additional manipulations.
Merozoites possess highly significant stage-specific mRNA expression differences from tachyzoites
Read mapping revealed a significant number of unique reads for 8323 genes after removal of reads mapping to tRNA and ribosomal RNA sequences, and normalization. Calculation of RPKM values (i.e. mapped reads per kilo base per million reads in dataset) to normalize for dataset size and gene length revealed 7148 genes with threshold RPKMs ≥10 in at least one of the two life-cycle stages. Comparative analysis (tachyzoite versus merozoite mapped reads per gene) was performed using DESeq  to detect differential gene expression. We used a high threshold of ≥ 8-fold difference in mRNA levels (measured as normalized averaged mapped sequence reads per gene from two datasets each for tachyzoites and merozoites) to identify stage-specifically regulated genes . The rigorous analytical approach revealed significant stage-specific expression differences for approximately 10% of genes annotated in ToxoDB (Figure 2A). Applying the ≥ 8-fold difference threshold we identified 453 genes with tachyzoite-specific, versus 312 genes with merozoite-specific expression. This strategy therefore revealed only the most strongly stage-specifically regulated genes; maximal differences of RNA expression values of 54-fold (ORF TGME49_295662, merozoite-specific) and 404-fold (ORF TGME49_215980, tachyzoite-specific) were measured.
RNAs with ≥ 8-fold difference in abundance at each life cycle stage were parsed, and distributions into eleven gene product or functional categories are shown as pie charts (Figure 2B, Additional file 1: Table S1A, B). The global comparison identified hypothetical genes as the most abundant set of stage-specific genes (146/47% in merozoites versus 294/64% in tachyzoites) expressed between the two developmental stages examined. The RNA-Seq data will also serve to annotate, and thereby significantly reduce, the number of hypothetical genes found in both samples. Genes related to metabolism or the cytoskeleton were not differentially regulated, with only a few exceptions. In contrast, the majority of genes implicated in host parasite interactions, such as secreted and surface proteins, were regulated in a stage-specific fashion (Figure 3). For example, significantly more surface (SRS) proteins were expressed in merozoites than in tachyzoites (40 versus 12), and their expression was restricted exclusively to the merozoite developmental stage. In striking contrast, the majority of microneme-, rhoptry-, and dense granule-specific secreted proteins were not expressed in merozoites, suggesting the existence of an unannotated suite of MIC, ROP and GRA proteins that promote merozoite asexual replication.
Merozoite specific gene families
Several large, unannotated gene families containing many predicted secreted and/or membrane-associated proteins, referred to as the T. gondii Family A-E proteins were identified to be largely merozoite stage-regulated (Additional file 2: Figure S1). Family A is predominantly merozoite-specific, with 29 of 33 members expressed above the ≥ 8-fold abundance threshold. Only one Family A gene was expressed exclusively in tachyzoites (Additional file 2: Figure S1A and Additional file 1: Table S2). Members of the remaining four families (with the exception of Family E) were likewise strongly expressed in merozoites (>4-fold change to tachyzoites), but just below the threshold (Additional file 2: Figure S1B). The “miscellaneous” category (Figure 2) also contains members of additional structural protein families many of whose members are strongly stage-regulated. One example is the Lysine-Arginine Rich Unidentified Function (KRUF) family of proteins  (Additional file 1: Table S3). Twelve of the 14 KRUF family members are strongly expressed in tachyzoites but not expressed (≥8-fold lower) in merozoites.
RNA-Seq data provides valuable information about the accuracy of existing gene models in ToxoDB . In this case the newly generated datasets complement existing RNA-Seq data, currently from tachyzoite, bradyzoite, and oocyst stages. To estimate how many gene models have no associated RNA-Seq data we used ToxoDB search strategies. While a detailed revision of ToxoDB gene models based on RNA-Seq data will have to await additional data from cat-derived gametocytes, we used simple ToxoDB in-silico search strategies with stringent criteria to identify merozoite-specific gene models that are currently without any evidence for expression. Several examples could be discovered by collating all predicted gene models for which no or insufficient (<40th percentile) evidence for RNA (n = 42) or protein expression (n = 38) exist at any stage of development and comparing them with the datasets developed in this study. We eliminated any entries with RPKM values >10 in tachyzoites and <10 in merozoites resulting in 10 genes for which significant expression is only documented in merozoites. The group comprises 6 Family A genes (Additional file 2: Figure S1), three hypothetical genes and SRS15B. Although such search strategies are difficult to threshold and require manual follow-up and verification, they represent examples of data mining efforts which are well within reach for experienced ToxoDB users. In addition, fortuitous discovery of large numbers of merozoite-specific RNA-Seq reads mapping to regions without predicted gene models, e.g. in a 4 kb region between TGME49_297820 and TGME49_297830 on chromosome II, indicates that systematic mapping of merozoite data will lead to discovery and annotation of many novel protein coding genes.
Metabolic capacity is similar in merozoites and tachyzoites
Merozoite populations expand rapidly by asexual division whilst remaining strictly confined to the enterocyte monolayer (Figure 1A, B) and presumably are confronted with less variability in terms of bioavailability of metabolites and building blocks than tachyzoites which infect different cell types. In addition, the parasite is ideally positioned to intercept the stream of nutrients flowing from the apical brush border to the basolateral face of the host enterocyte. In the absence of data on enzyme activity levels, mRNA expression can be used effectively to model metabolic flux distributions and stage-specific changes [12,13]. Global comparative analysis of metabolic gene expression was done on the subset of genes used for constructing the iCS382 metabolic model for Toxoplasma . Although we identified 11 “tachyzoite-specific” and 7 “merozoite-specific” genes based on >8-fold difference in mRNA expression within this subset (Additional file 2: Figure S2), there is little evidence in our data that any of these were clustered either on chromosomes or in any particular metabolic pathway. Nevertheless, several genes of the energy metabolism are strongly regulated, with the glycolytic enzyme enolase 1 (ENO1, TGME49_268860) showing the highest difference in expression, albeit with relatively low RPKM values, and is effectively silenced in merozoites (RPKM 80 in tachyzoites, vs 1 in merozoites, adjusted p value: 1.74e-7). ENO1 is a bradyzoite marker and thus not well expressed in tachyzoites to begin with. Indeed, ENO1 appears to be dispensable in tachyzoites: cells with a targeted deletion of the gene grow well in vitro and in vivo, but cyst numbers in chronically infected mice are significantly decreased . Conversely, the ENO2 (a canonical tachyzoite marker) gene is strongly downregulated in bradyzoites. In merozoites ENO2 mRNA is detected at ~2-fold higher levels (RPKM 620) than in tachyzoites, which points to an adaptive response geared towards a higher glycolytic throughput compared to tachyzoites. Although a more detailed analysis of energy metabolism would require RNA-Seq data for these enzymes in all coccidial stages and also bradyzoites integrated into a flux model such as iCS382 , it is tempting to speculate that these adaptations in merozoites help sustain rapid growth in conditions of low oxygen tension in the gut. In addition to their function in glycolysis, nuclear targeted ENO isoenzymes were shown bind to chromosomal DNA and modulate gene expression . The significance of low levels of ENO1 expression in the tachyzoite dataset and the tight silencing of the gene in merozoites in the context of their regulatory functions remain to be investigated. The stage-specific regulation of ENO expression is also consistent with 3.5 fold higher levels of phosphoglycerate kinase PGKI (TGME49_318230) RNA-Seq reads in merozoites. In line with this, differential expression of the two glycolytic lactate dehydrogenase isoenzymes LDH1 and LDH2  in merozoites is very similar to tachyzoites but even more pronounced: virtual silencing of LDH2 (RPKM 2) but strong expression of LDH1 (RPKM 962 vs 472 in tachyzoites). There are also several examples for tachyzoite-specific expression of metabolic factors: both NADPH-generating glucose-6-phosphate 1-dehydrogenase (TGME49_294200), an enzyme of the pentose phosphate pathway, and 8-amino-7-oxononanoate synthase (TGME49_290970), a key enzyme in the biotin synthesis pathway, are not expressed in merozoites (RPKMs 5 and 1, respectively). Biotin provides the prosthetic group for carboxyl transferases and has essential functions in many anabolic and catabolic reactions. Because biotin is contained in food but also supplied by gut microorganisms, a complete shutdown of its synthesis in merozoites suggests that developmental stages in enterocytes rely entirely on scavenging the easily accessible protein.
The most strongly regulated metabolic gene in merozoites relative to tachyzoites is glycerol-3-phosphate dehydrogenase (NAD+) (TGME49_210260), with essentially no expression in the latter (RPKM 12) but a 20-fold higher expression in merozoites (adjusted p-value: 2.7e-21) (Additional file 2: Figure S2). RNA-Seq reads for the homoserine O-acetyltransferase gene coding for an enzyme of the methionine biosynthesis pathway are similarly increased in merozoites. Interestingly, expression of two genes involved in purine metabolism, adenylosuccinate synthetase, (a.k.a., IMP-aspartate ligase, TGME49_279450), and hypoxanthine-xanthine-guanine phosphoribosyl transferase (HXGPRT, TGME49_200320) is > 8-fold higher in merozoites. The latter, a key enzyme in the purine salvage pathway , is expressed at a very high level in merozoites with an RPKM value of 882 compared to 106 in tachyzoites. Another purine metabolism gene, purine nucleoside phosphorylase (PNP, TGME49_307030), is expressed 6-fold higher in merozoites (RPKM 88, adjusted p-value: 1.9e-2), suggesting a possible expansion of purine salvage throughput at least during the asexual phase of parasite population growth during enteroepithelial development.
Altogether, we conclude that in contrast to the marked differences in gene expression of secreted proteins in tachyzoites and merozoites, the expression of genes involved in metabolic processes is largely unaffected. Hence, metabolism is predicted to be broadly similar between tachyzoites and merozoites, although there may be some fine-tuning in the latter for optimal growth in enterocytes. However, no fundamental differences were detectable that would indicate a radically different environment or nutrient availability for growth in enterocytes. A more detailed examination and testing of specific hypotheses using for example the recently developed flux balance analysis models to identify stage-specific metabolic bottlenecks  will be required to confirm or refute this current interpretation as well as identify potential drug targets.
SRS proteins are widely expressed in merozoites
The T. gondii genome encodes several distinct, coccidian-specific surface gene families, including the SRS and SAG-unrelated surface antigens (SUSA) . Arguably the most strikingly regulated set of genes in the merozoite dataset were those coding for the 111 members of the SRS superfamily of proteins annotated in ToxoDB (Version 8.1); expression of more than 52 members of this family was present in merozoites, whereas a separate set of 14 SRS genes were expressed exclusively in tachyzoites (Figure 3). Nearly half of the merozoite expressed genes coding for SRS proteins were present in 5 clusters, SRS12, 15, 22, 26 and 55. SRS proteins are involved in attachment to host cells, but also provoke immune reactions and regulate parasite virulence, which is thought to promote the formation of tissue cysts in intermediate hosts in order to establish persistent, latent infections that facilitate transmission of infection to the definitive host. Previous work has showed that tachyzoites differentially express a number of SRS genes , and these expression differences have been postulated to account for the ability of this stage to invade a broader range of host cells than other coccidians . However, merozoites, which only infect a single cell type (the feline enterocyte) co-dominantly expressed a large repertoire of 52 SRS proteins in a developmental life-cycle stage-dependent manner (Additional file 2: Figure S3). Alternatively, recent work has established that the SRS fold is present in the 10 member Pfs-230-related 6-Cys family of Plasmodium adhesins  that facilitate gamete-gamete recognition and promote gamete recognition and fertilization. This may suggest that the merozoite-restricted SRS genes are less relevant for attachment and invasion within enterocytes, but rather promote gamete development and fertilization. Given the demonstrated ability of SRS proteins to promote immune responses, it is also conceivable that they play a role in stimulating intestinal inflammation and diarrhea, to facilitate the production and dispersal of oocysts.
The striking stage-specificity of expression of SRS genes in merozoites versus tachyzoites (and as previously observed for merozoites versus bradyzoites) raises questions about how their expression is regulated. One clue to this can be obtained from analysis of chromosomal distribution of the various SRS genes [10,18]. However, we find little evidence of concerted clustering of merozoite- or tachyzoite-specific SRS genes on distinct chromosomes (Figure 4), which is not a complete surprise since T. gondii SRS genes are dispersed across all chromosomes and not exclusively in a subtelomeric regional distribution pattern for surface antigen genes that is characteristic for other protozoa [20-23]. Thus, rather than relying on sequential expression of a single gene or promoting rapid ectopic recombination rates to generate surface antigenic variation, T. gondii expresses relatively large sets of non-overlapping, co-dominant SRS transcripts [18,24].
Merozoites express a distinct set of microneme proteins
Micronemes are small organelles clustered at the anterior end of apicomplexan parasites and generally secrete adhesive proteins that are important for motility and invasion. Microneme (MIC) proteins contain one or several adhesive domains such as thrombospondin, von Willebrand factor A, epidermal growth factor (EGF) or PAN-domains (reviewed in ). Comparative analysis of RNA-Seq data shows that the majority of the 26 Toxoplasma MICs and MIC2 associated protein (MIC2AP) accumulate a high number of mapped tachyzoite reads (Additional file 2: Figure S4). Conversely, the range of MIC genes that are specifically expressed in merozoites is considerably reduced and restricted to the highly expressed genes coding for PAN-domain containing proteins MIC17A-C (Additional file 2: Figures S4 and S5), MIC12, a gene annotated as a putative MIC, and a microneme-like protein. Interestingly, both MIC2 and MIC2AP , as well as MIC6-1-4 , which are strongly expressed in tachyzoites and secreted as functional complexes, are virtually silenced in merozoites (Figure 5). Similarly, the RNA-Seq reads for the MIC3 escorter MIC8  are >100-fold lower in merozoites, whilst MIC3 mRNAs are only 6-fold less abundant. Since trafficking to the organelles, as well as deployment on the surface, seems to occur in complexes as a rule , we analyzed the predicted topologies of the six merozoite-specific MIC proteins. None had the predicted transmembrane domain and cytoplasmic tail typical for canonical escorters such as MIC6 . In tachyzoites these adhesin complexes are essential components of the molecular motor and play a key role in providing the machinery for gliding motility .
The number of microneme organelles of different Apicomplexa has been correlated with the requirement for gliding motility , ranging from many in Eimeria sporozoites and merozoites  to none in Theileria zoites . Toxoplasma tachyzoites, which display vigorous gliding motility and are able to cross cell barriers, have many microneme organelles and express a wide range of MIC proteins. Although T. gondii merozoites have a high number of microneme organelles , the complement of expressed MIC and other genes (e.g. AMA1 ) coding for proteins secreted from micronemes appears to be significantly restricted, suggesting that the size of the expressed protein complement is not necessarily correlated to the number of organelles present in the cell. On the other hand, the finding that prominent tachyzoite MIC complexes are not expressed in merozoites (Figures 3 and 5), could be explained by fundamental differences in their mode of gliding motility and cell invasion in contrast to the extraintestinal developmental stages. Because merozoite proliferation is restricted to enterocytes, the most obvious difference is that mature merozoites egress into the gut lumen and invade new host enterocytes from the exposed apical side. This requires some migration through intestinal contents and within the mucus layer rather than through intercellular spaces and on cell surfaces. Because freshly isolated extracellular merozoites display vigorous gliding motility (data not shown) it is unlikely that these stages are without surface adhesins linked to the glideosome complex [31,35]. Rather, as yet unknown MIC adhesins specific for enterocyte surface receptors and containing non-canonical adhesive domains will likely be discovered among hypothetical proteins that are expressed in a merozoite stage-specific manner.
Merozoites only weakly express genes of the tachyzoite moving junction complex
There is a striking difference in expression of genes coding for rhoptry neck proteins (RONs) and AMA1 in tachyzoites versus merozoites (Figure 6). In particular, mRNA levels of genes coding for proteins that combine to form the moving junction as part of the tachyzoite invasion machinery (i.e., AMA1, RON2, RON4, RON 5 and RON8; Figure 6) [36,37] are dramatically lower in merozoites. This supports the contention, raised above, that the mechanism of gliding motility and by extension cell invasion is distinct in merozoites. The dramatically lower expression of the canonical tachyzoite AMA1 (TGME49_255260) in merozoites is not compensated by a significant increase in expression of the sporozoite-specific AMA1 paralog and the associated sporoRON2  is not expressed. A third Toxoplasma AMA1 paralog (TGME49_300130) is weakly expressed in both stages, and its functional role in invasion, if any, remains to be demonstrated. Interestingly, previous claims that AMA1 is essential for tachyzoite invasion  has been revised in light of recent data showing that AMA1 knockout tachyzoites invade host cells normally in vitro . On the one hand, gene expression data alone predicts that invasion via the massively expanded surface of the enterocyte brush border by merozoites is distinct from the well described process observed in extraintestinal stages. On the other, the exact requirements and concepts for formation of a tight junction [36,41-44] which will be moved towards the posterior end of the invading parasite by motor proteins are currently being revisited [40,45]. While gene expression data indicates fundamental differences in cell invasion by merozoites, the question whether merozoites assemble an alternative moving junction machinery for invasion or invade enterocytes by an entirely different mechanism, e.g. based on membrane-zippering such as in Theileria  or Neospora , requires further investigation.
Merozoites express a distinct set of rhoptry proteins
We also investigated RNA expression for a group of 57 genes coding for rhoptry proteins (ROPs) compiled from annotated gene models on ToxoDB (Version 8.1)  and expanded by manual annotation (Additional file 1: Table S4). In total 24 ROP genes were significantly differentially expressed (≥8-fold abundance threshold) between tachyzoites and merozoites (Figure 3 and Additional file 2: Figure S6); 19 genes were considered not expressed in merozoites (RPKM values were <10 in all cases) whereas two were expressed at very low levels (RPKM < 40) in merozoites. Four ROP genes were, however, expressed specifically in merozoites. Interestingly, in a survey of known functions of ROP proteins, we find that the prominent murine virulence factors ROP16 and ROP18/5 are not expressed in merozoites.
ROP16 and 18, which both possess a conserved catalytic triad and phosphoryltransfer activity, have been identified as secreted, strain-specific virulence factors in mice. ROP16 phosphorylates the signal transducers and activators of transcription (STATs) 3 and 6, and regulates interleukin 12 levels [49,50]. The polymorphic ROP18 kinase phosphorylates and thus inhibits loading of immunity-related GTPase (IRG) 6 onto the PV membrane, thereby preventing parasite clearance in infected mice by activated macrophages [51,52]. ROP5, a pseudokinase with an altered catalytic triad, regulates the kinase activity of ROP18  possibly by an allosteric mechanism. Therefore, the absence of IRGs in cats is the simplest explanation for the silencing of the ROP18 (and possibly the ROP5) gene in merozoites rendering expression and secretion of these ROPs as part of a defense strategy redundant during enteroepithelial development.
Merozoites express a distinct set of dense granule proteins
We detected significant differences in expression within the current set of 25 annotated genes encoding dense granule (GRA) proteins (Figure 3, Additional file 2: Figure S7). There was a significant overall bias towards high expression in tachyzoites. Twelve GRA genes were not expressed in merozoites (RPKM < 10), whilst only 2 GRA genes (GRA11 isoforms TGME49_212410 and TGME49_237800, which have, as yet, undefined functions) are >8-fold higher expressed in cat stages (18- and 13-fold, respectively; Figure 3). The high expression of GRA22 in tachyzoites is intriguing, since this protein regulates parasite egress from host cells  and the absence of GRA22 expression in merozoites is therefore puzzling. This might imply that egress from cat enterocytes is less coordinated by the parasite and this may partially explain why parasite development is asynchronous in the cat intestine. GRA genes expressed in tachyzoites fall into two main groups: moderately highly expressed GRA genes (RPKM values in the range of 100) code for proteins that are exported across the PV membrane that function to alter host gene expression (i.e., GRA 15, 16, and 24); versus highly expressed GRA genes (RPKMs > 103) some of whose products participate in the formation of the tubulovesicular network (TVN) within the PV space. Our RNASeq data shows that GRA15 is not expressed in merozoites. This tachyzoite virulence factor has a well-defined role in the activation of several host cell pathways that trigger potent pro-inflammatory responses causally associated with tachyzoite to bradyzoite conversion [49,55-59]. In mice, parasite survival is dependent on combined expression of GRA15 with the STAT-activating form of ROP16 . Both GRA15 and ROP16 genes are not expressed in merozoites, strongly suggesting that JAK-STAT3/6 and NF-κB signaling pathways in enterocytes of the definitive host are not targeted by these secreted factors. GRA24 and 16 are also not expressed in merozoites. GRA24 was recently identified as a parasite soluble effector that traffics to the host cell nucleus, interacts with p38, and is thought to synergize with GRA15 to promote production of the pro-inflammatory cytokine, IL12 . Similarly, GRA16 is exported to the host cell nucleus and affects transcription of host cell genes involved in metabolism and cell cycle processes . Whether these host cellular processes are not targeted by merozoite gene products during their asexual expansion in cat enterocytes, or if as yet unannotated merozoite-specific paralogs exist to perform this function, remain to be established.
The RPKM values for GRA2, 4, and 6 genes, coding for three proteins which are directly associated with the TVN in the PV [62,63], are >8-fold lower in merozoites compared to tachyzoites (Additional file 2: Figure S7). The silencing of these GRAs (RPKMs <10 and log -fold change ≥ 8) in merozoites is consistent with electron microscopy data  showing that type C-E schizonts in cat enterocytes do not elaborate the extensive TVN which is typical of vacuolar spaces of tachyzoites, bradyzoites and type A and B schizonts. This suggests that later stage schizonts dispense with the TVN because of their different mode of cell division ; parasites in type C-E schizonts multiply by endopolygeny rather than by endodyogeny as early stage merozoites, tachyzoites and bradyzoites. Magno et al.  suggested that an important structural function of the TVN might be organization of daughter cells within the PV during multiplication. Involution of the TVN and repression of genes coding for factors associated with these membranes in merozoites is consistent with endopolygeny and attachment of parasites to residual bodies . Conversely, the GRA7 gene, whose 36 kDa product is involved in nutrient acquisition by trapping host endosomal-lysosomal vesicles in the PV , is expressed at ~4-fold lower levels in merozoites (Additional file 2: Figure S7) but is still relatively abundant with an RPKM value of 351 (1837 in tachyzoites). Thus, although not always associated with the PV membrane, GRA7 appears to be universally expressed throughout the life cycle .
Additional factors secreted by dense granules include Kazal-type protease inhibitors PI1, PI2 (TgME49_217430, TgME49_208450) [67-69] (Table S5), which are secreted into the PV and appear to be differentially regulated, despite strong conservation on the structural and sequence level. Whilst read numbers for the PI1 gene, the most abundantly expressed protease inhibitor in tachyzoites , are >2-fold higher in merozoites, PI2 is not expressed in merozoites (RPKMs 50-fold lower). Three of five additional genes annotated as Kazal-type family protease inhibitors (Table S5) show significantly higher expression in merozoites. However, whether the four predicted secreted members of this group also traffic via dense granules is unknown.
Trans and cis regulators of merozoite and tachyzoite biology
Both trans and cis regulators are known to play crucial roles in the transcription of stage specific genes in T. gondii. The trans-acting ApiAP2 transcription factors were discovered in the Apicomplexa almost 10 years ago  and regulate various developmental processes in T. gondii, including bradyzoite conversion, with TgAP2XI-4 , and TgAP2IX-9  promoting and restricting tissue cyst formation, respectively, as well as the tachyzoite cell cycle, with TgAP2XI-5 potentially regulating hundreds of cell cycle-dependent genes . Using stringent criteria we identified three significantly merozoite- or tachyzoite-specifically expressed ApiAP2-coding genes in our study (Additional file 2: Figure S8). One merozoite-specific gene TgAP2VIIa-I, and two tachyzoite-specific genes, TgAP2XII-I and TgAP2VIIa-II showed >8-fold difference in their expression levels, suggesting that these transcription factors may play a role in regulating a suite of stage-specific genes. Relaxing the criterion to include genes differing ≥4-fold in mRNA abundance identified 8 additional tachyzoite-, and 2 additional merozoite-specific genes (a total of 13 regulated stage-specific TgAP2s). Interestingly, several other TgAP2 genes were expressed at very low levels, and in the case of TgAP2XII-3, can be considered silenced in both stages. Overall, differential expression of TgAP2 is not extensive between tachyzoite and merozoites, and this may reflect their similar life cycle fate as rapidly amplifying asexual stages. However, the evenly graded distribution of RPKM values (Additional file 2: Figure S8B) suggests that some differences may become more pronounced when merozoites differentiate to gametocytes. Moreover, the two datasets compared here do not represent synchronized populations, which would average out any dynamics occurring during the merogony cell cycle specifically during development from stage A through stage E schizonts. By contrast, 24 of the 67 predicted TgAP2-coding genes were found to be differentially expressed at different check-points of the tachyzoite cell cycle in a microarray analysis of synchronized parasites . In addition, ApiAP2s have been shown to bind other transcriptional regulators such as histone lysine acetyltransferases , histone deacetlyase  and even other ApiAP2s , producing interactions that may confound correlations between transcription factor activities and their cognate gene transcript levels.
Next, we used the Regulatory Sequence Analysis Tools on-line computational resource (http://www.rsat.eu/) to perform a global analysis of upstream flanking sequences of stage-specifically expressed genes. The analysis returned several enriched six-base DNA motifs within the predicted promoters, i.e. 500 bp genomic regions upstream of the transcriptional start site of either merozoite- or tachyzoite-specific genes. Listed in Figure 7A are the five most significantly enriched 6 bp motifs found. Pattern matching between these enriched motifs produced an amalgamated seven-base DNA motif in both merozoite- and tachyzoite-specific gene promoters, representing putative cis-regulatory elements (Figure 7B and C, respectively). The GAAGAAA motif, present in 21.5% of merozoite promoters, was also present in 19.2% of tachyzoite promoters, indicating that this motif is unlikely to represent a merozoite-specific cis-regulatory element (Figure 7D). However, the GAGACGC motif is clearly enriched in the promoters of tachyzoite-specific genes (present in 20.4% of promoters) over merozoite-specific genes (6.5%). The distribution of this putative cis-regulatory element in the predicted promoters of the 15 most tachyzoite-specific genes is illustrated in Figure 7E. This motif is nearly identical to a motif (A/TGAGACG) previously identified as a functional cis-element in the promoters of dense granule genes, such as TgGRA1, TgGRA2, TgGRA5 and TgGRA6, as well as the tachyzoite marker, TgSRS29B (aka TgSAG1) . A sequence-specific transcription factor capable of activating gene transcription via this cis element has yet to be identified; however, it may play a significant role in promoting and/or maintaining a tachyzoite developmental state.
Taken together, our data support a correlation between promiscuity of invasion and dissemination of the T. gondii parasite in a variety of host cells and tissues. In addition, the data presented here underscore the phenomenon of non-overlapping stage-specific expression of gene sets throughout development (e.g. SRS, ROP, GRA, MIC) as previously proposed  within the Toxoplasma life cycle . Hence, the differentiation of bradyzoites into merozoites, which occurs only in the feline intestine, entails the shut down or downregulation of suites of tachyzoite-specific genes that would normally be activated to support Toxoplama’s promiscuous ability to invade and reproduce asexually in virtually any nucleated cell from essentially any warm-blooded animal. This study provides compelling datasets to begin to explore the suites of genes activated to support merozoite invasion, differentiation, and expansion in cat enterocytes and establishes T. gondii as an excellent genetically-amenable model system to study schizont asexual amplification prior to gamete development, the hallmark of coccidian development in the definitive hosts these coccidian parasites infect. Our data suggest that gliding motility and the formation of a moving junction (or an alternative invasion apparatus) requires machinery that is specifically adapted to breaching the enterocyte surface barrier for every cycle of replication, egress and invasion. Studies of the SRS proteins exclusively expressed during coccidian development will determine if these proteins promote attachment and invasion of feline enterocytes, or whether they participate in gamete-gamete recognition to promote fertilization, analogous to the orthologous role performed by related 6-Cys family proteins during Plasmodium development . Understanding the governance of gene expression and regulation within these families, e.g. via interactions with the TgAP2 family of transcriptional regulators, may prove key to dissecting the full panorama of pathways that regulate the definitive coccidian lifecycle and aid in the complete and accurate annotation of the T. gondii genome database.
Isolation of enteroepithelial parasite stages from a cat infected with bradyzoites from sheep
Oocysts of the CZ T. gondii isolate (type II) were originally isolated from the feces of a captive Siberian tiger (Panthera tigris altaica) at the Dvůr Králové Zoo (Czech Republic). The strain isolation was performed in 2005 by Dr. B. Koudela, and the strain was maintained by passages between mice and cats [81,82]. After obtaining the strain from Dr. Koudela the parasites were passed once through mice and cats to generate ~50 million CZ oocysts, which were used to infect four sheep (104 oocysts per sheep). Equal amounts of brain tissue containing bradyzoites in tissue cysts from three chronically infected sheep was homogenized and fed to cats (55 g per cat). A semi-quantitative analysis of parasite burden in sheep brain tissue was done by magnetic capture real time PCR, to confirm that equal quantities of parasites were inoculated into each cat but no definitive count of cysts was performed. Pilot experiments had established that this protocol results in strong and reproducible infections of cats. All feces from experimentally infected cats were collected and oocysts were enumerated by microscopy every day. For isolation of enteroepithelial stages of T. gondii a cat that was experimentally infected with brain cysts from a sheep was euthanized after onset of patency (day 5 post infection). The small intestine was clamped, removed, and immediately cooled to 0°C. All subsequent procedures for parasite isolation and purification were performed in a cold room (4°C) and on ice. The small intestine (~70 cm) was carefully rinsed with PBS, divided into four sections of equal length, and opened. Enterocytes were selectively removed from villi by gentle scraping with a rubber policeman (Figure 1B). Tissue samples for histology were taken before and after removal of enterocytes. The cells containing numerous parasites were washed off with ice cold PBS and collected by centrifugation (10 min, 900 × g, 2°C). Enterocyte preparations were evaluated microscopically after fixation and staining to identify parasite stages. The two enterocyte samples collected from two different sections of the small intestine were washed with buffer containing Tween80 and syringe-passaged The samples appeared to consist exclusively of free merozoites as no other stages (i.e. gamonts) could be identified by microscopy (Figure 1C).
Histological and fluorescence microscopy examination of small intestinal tissue sections revealed numerous parasite stages that were readily identified as merozoites, schizonts (Figure 1A) and also occasional gamonts due to their very distinctive morphologies. However, the final preparation was clearly highly enriched for merozoites allowing us to focus on interrogating RNA expression datasets specifically for differences between asexually developing merozoites and tachyzoites. The enrichment for merozoites or, more precisely, the paucity of gametocytes in the two biological replicates used for RNA-Seq, was indeed borne out by a lack of known gametocyte markers (e.g. genes coding for flagellar or oocyst wall proteins) in the RNA-Seq datasets. This is in contrast with a preliminary gene expression analysis of stages isolated at or past peak patency (Ramakrishnan et al. unpublished). CZ strain tachyzoites were grown in vitro; confluent human foreskin fibroblast (HFF) cells were infected at an MOI (>5) and harvested from synchronized cultures when 1/3 of the vacuoles had spontaneously lysed.
Experiments involving animals were performed under the direct supervision of a veterinary specialist, and according to Swiss law and guidelines on Animal Welfare and the specific regulations of the Canton of Zurich. Permit number 108/2010 covers all animal experiments presented in this paper and was approved by the Veterinary Office and the Ethics Committee of the Canton of Zurich (Kantonales Veterinäramt Zürich, Zollstrasse 20, 8090 Zürich, Switzerland).
RNA isolation and quality control
Infected tissue and parasite samples were resuspended in 600 μl buffer RLT (QIAGEN) containing 10 μl/ml β-mercaptoethanol. Cell suspensions were then passed through a QIAshredder (QIAGEN) column by centrifugation at ≥ 8000 g for 1 min. The RNA was then purified using the RNeasy Mini Kit (QIAGEN) according to the manufacturer’s protocol (including an on-column DNA digest), and eluted in RNase-free water. The quality of the RNA was analysed using the Agilent RNA 6000 Pico Kit (Agilent) and a Bioanalyzer 2100 (Agilent) (Figure 1D). RNA concentration was determined using a Qubit fluorometer (Invitrogen) together with the RNA assay (Invitrogen).
The genome-wide transcriptome library, was produced using the MicroPolyA Purist Kit (Ambion) and the SOLiD Total RNA-Seq kit (Applied Biosystems). Briefly, approximately 300-500 ng of mRNA were enriched starting from 15-20 μg of total RNA, using MicroPolyA Purist Kit (two rounds of purification to receive only minor ribosomal RNA contamination). The quality and the quantity of the extracted polyA RNA was assessed using a Bioanalyzer (Agilent) Picochip and a Qubit fluorometer (Invitrogen), respectively. The mRNA was then fragmented using RNase III. Ligation of the adaptor mix and reverse transcription were performed following the manufacturer’s protocol. cDNA libraries were selected for fragment sizes between 150 and 250 bp, amplified for 15-18 cycles of PCR using barcoded adaptor primers and purified using the PureLink PCR micro kit (Invitrogen). Library size and concentration were assessed using a Bioanalyzer (Agilent) and a Qubit fluorometer (Invitrogen), respectively.
The poly-A transcriptome libraries were used for emulsion PCR (e-PCR) using a concentration of 0.5 pM. The barcoded libraries were pooled before e-PCR and the resultant beads were loaded on a full SOLiD 5 slide (Applied Biosystems), according to manufacturer’s instructions. SOLiD ToP Sequencing chemistry was used to produce paired-end (50 bp + 35 bp) sequencing reads.
Comparative transcriptomics of merozoites and tachyzoites
We generated RNA-Seq datasets from CZ strain enteroepithelial parasites isolated from mechanically stripped enterocytes and from tachyzoites cultured in human foreskin fibroblasts on a SOLiD5 platform (Applied Biosystems). Two biological samples were sequenced for each stage (two separate parasite and RNA preparations from different sections of the small intestine as well as two different tachyzoite cultures). The samples contained different proportions of parasite and host RNA (Figure 1D). Adjustments to account for these differences were: i) removal of all reads mapping to tRNAs, and ribosomal RNAs prior to DESeq analysis, and ii) normalization and weighting by using transcripts lengths to calculate RPKMs. The reads in colorspace were aligned to the T. gondii ME49 genome from ToxoDB v8.1 using SHRiMP2 mapper [83,84].The count table of reads fitting the gene coordinates have been processed by Bioconductor library DESeq [8,85] to produce the table that includes the mean read number for each gene together with a p-value and fold change. We have applied an 8-fold (log2 ≥ ) change as a cut-off for differential expression.
Immunofluorescence assay and microscopy
Samples from the small intestine of infected cats were embedded in paraffin. Sections of 3 μm were placed on slides with a positively charged coat. After drying over night at 37°C, slides were deparaffinized 3 times 2 min in xylene and then washed twice 1 min in 100% ethanol. Slides were then washed 1 min each in 96%, 70% ethanol, and water. Slides were then submerged in alkaline Target Retrieval Solution pH 9 (Dako), boiled for 20 min in a pressure cooker at 96°C and transferred to water. Permeabilization was done in PBS/0.3% Triton X-100 (permeabilization buffer, PB) for 5 min. Blocking was performed in rabbit serum for at least 1 h at RT. After two washes in PB for 5 min, sheep immune serum against T. gondii diluted in 20% rabbit serum/PB was used to label the parasites by overnight incubation at 4°C. Slides were washed three times with PB for 5 min and samples were incubated with rabbit anti-sheep FITC in 20% rabbit serum/PB for 1 h at RT. Slides were then washed once in PB and the nuclei were counter stained with 1 μg/ml DAPI in PB for 15 min. Slides were washed again twice with PB for 5 min and the samples were mounted using Vectashield (Vector Laboratories). Imaging was performed using a Leica DMI 6000 B microscope and the Leica LAS AF software. Images were processed using ImageJ version 1.47.
Availability of supporting data
The raw sequence data and a complete description of the study are available at the European Nucleotide Archive (Study accession no. PRJEB7935). The data are also available via Toxoplasma Genomics Resource, ToxoDB, [http://www.toxodb.org/toxo/]. All other data sets supporting the results of this article are included within the article and its additional files.
Halonen SK, Weiss LM. Toxoplasmosis. Handb Clin Neurol. 2013;114:125–45.
Jones JL, Dubey JP. Foodborne toxoplasmosis. Clin Infect Dis. 2012;55(6):845–51.
Swierzy IJ, Muhammad M, Kroll J, Abelmann A, Tenter AM, Luder CG. Toxoplasma gondii within skeletal muscle cells: a critical interplay for food-borne parasite transmission. Int J Parasitol. 2014;44(2):91–8.
Ferguson DJ. Use of molecular and ultrastructural markers to evaluate stage conversion of Toxoplasma gondii in both the intermediate and definitive host. Int J Parasitol. 2004;34(3):347–60.
Speer CA, Dubey JP. Ultrastructure of early stages of infections in mice fed Toxoplasma gondii oocysts. Parasitology. 1998;116(Pt 1):35–42.
Kim K, Weiss LM. Toxoplasma: the next 100 years. Microbes Infect. 2008;10(9):978–84.
Torrey EF, Yolken RH. Toxoplasma oocysts as a public health problem. Trends Parasitol. 2013;29(8):380–4.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Fritz HM, Buchholz KR, Chen X, Durbin-Johnson B, Rocke DM, Conrad PA, et al. Transcriptomic analysis of toxoplasma development reveals many novel functions and structures specific to sporozoites and oocysts. PLoS One. 2012;7(2):e29998.
Reid AJ, Vermont SJ, Cotton JA, Harris D, Hill-Cawthorne GA, Konen-Waisman S, et al. Comparative genomics of the apicomplexan parasites Toxoplasma gondii and Neospora caninum: Coccidia differing in host range and transmission strategy. PLoS Pathog. 2012;8(3):e1002567.
Gajria B, Bahl A, Brestelli J, Dommer J, Fischer S, Gao X, et al. ToxoDB: an integrated Toxoplasma gondii database resource. Nucleic Acids Res. 2008;36(Database issue):D553–6.
Colijn C, Brandes A, Zucker J, Lun DS, Weiner B, Farhat MR, et al. Interpreting expression data with metabolic flux models: predicting Mycobacterium tuberculosis mycolic acid production. PLoS Comput Biol. 2009;5(8):e1000489.
Huthmacher C, Hoppe A, Bulik S, Holzhutter HG. Antimalarial drug targets in Plasmodium falciparum predicted by stage-specific metabolic network analysis. BMC Syst Biol. 2010;4:120.
Song C, Chiasson MA, Nursimulu N, Hung SS, Wasmuth J, Grigg ME, et al. Metabolic reconstruction identifies strain-specific regulation of virulence in Toxoplasma gondii. Mol Syst Biol. 2013;9:708.
Mouveaux T, Oria G, Werkmeister E, Slomianny C, Fox BA, Bzik DJ, et al. Nuclear glycolytic enzyme enolase of Toxoplasma gondii functions as a transcriptional regulator. PLoS One. 2014;9(8):e105820.
Yang S, Parmley SF. Toxoplasma gondii expresses two distinct lactate dehydrogenase homologous genes during its life cycle in intermediate hosts. Gene. 1997;184(1):1–12.
Chaudhary K, Darling JA, Fohl LM, Sullivan Jr WJ, Donald RG, Pfefferkorn ER, et al. Purine salvage pathways in the apicomplexan parasite Toxoplasma gondii. J Biol Chem. 2004;279(30):31221–7.
Wasmuth JD, Pszenny V, Haile S, Jansen EM, Gast AT, Sher A, et al. Integrated bioinformatic and targeted deletion analyses of the SRS gene superfamily identify SRS29C as a negative regulator of Toxoplasma virulence. MBio. 2012;3(6):1–13.
Tonkin ML, Arredondo SA, Loveless BC, Serpa JJ, Makepeace KA, Sundar N, et al. Structural and biochemical characterization of Plasmodium falciparum 12 (Pf12) reveals a unique interdomain organization and the potential for an antiparallel arrangement with Pf41. J Biol Chem. 2013;288(18):12805–17.
Myler PJ, Allison J, Agabian N, Stuart K. Antigenic variation in African trypanosomes by gene replacement or activation of alternate telomeres. Cell. 1984;39(1):203–11.
Rubio JP, Thompson JK, Cowman AF. The var genes of Plasmodium falciparum are located in the subtelomeric region of most chromosomes. Embo J. 1996;15(15):4069–77.
Freitas-Junior LH, Bottius E, Pirrit LA, Deitsch KW, Scheidig C, Guinet F, et al. Frequent ectopic recombination of virulence factor genes in telomeric chromosome clusters of P. falciparum. Nature. 2000;407(6807):1018–22.
Marty AJ, Thompson JK, Duffy MF, Voss TS, Cowman AF, Crabb BS. Evidence that Plasmodium falciparum chromosome end clusters are cross-linked by protein and are the sites of both virulence gene silencing and activation. Mol Microbiol. 2006;62(1):72–83.
Jung C, Lee CY, Grigg ME. The SRS superfamily of Toxoplasma surface proteins. Int J Parasitol. 2004;34(3):285–96.
Carruthers VB, Tomley FM. Microneme proteins in apicomplexans. Subcell Biochem. 2008;47:33–45.
Rabenau KE, Sohrabi A, Tripathy A, Reitter C, Ajioka JW, Tomley FM, et al. TgM2AP participates in Toxoplasma gondii invasion of host cells and is tightly associated with the adhesive protein TgMIC2. Mol Microbiol. 2001;41(3):537–47.
Reiss M, Viebig N, Brecht S, Fourmaux MN, Soete M, Di Cristina M, et al. Identification and characterization of an escorter for two secretory adhesins in Toxoplasma gondii. J Cell Biol. 2001;152(3):563–78.
Kessler H, Herm-Gotz A, Hegge S, Rauch M, Soldati-Favre D, Frischknecht F, et al. Microneme protein 8–a new essential invasion factor in Toxoplasma gondii. J Cell Sci. 2008;121(Pt 7):947–56.
Meissner M, Reiss M, Viebig N, Carruthers VB, Toursel C, Tomavo S, et al. A family of transmembrane microneme proteins of Toxoplasma gondii contain EGF-like domains and function as escorters. J Cell Sci. 2002;115(Pt 3):563–74.
Sheiner L, Santos JM, Klages N, Parussini F, Jemmely N, Friedrich N, et al. Toxoplasma gondii transmembrane microneme proteins and their modular design. Mol Microbiol. 2010;77(4):912–29.
Opitz C, Soldati D. ‘The glideosome’: a dynamic complex powering gliding motion and host cell invasion by Toxoplasma gondii. Mol Microbiol. 2002;45(3):597–604.
Scholtyseck E, Mehlhorn H. Ultrastructural study of characteristic organelles (paired organelles, micronemes, micropores) of sporozoa and related organisms. Z Parasitenkd. 1970;34(2):97–127.
Shaw MK. The same but different: the biology of Theileria sporozoite entry into bovine cells. Int J Parasitol. 1997;27(5):457–74.
Hehl AB, Lekutis C, Grigg ME, Bradley PJ, Dubremetz JF, Ortega-Barria E, et al. Toxoplasma gondii homologue of plasmodium apical membrane antigen 1 is involved in invasion of host cells. Infect Immun. 2000;68(12):7078–86.
Meissner M, Schluter D, Soldati D. Role of Toxoplasma gondii myosin A in powering parasite gliding and host cell invasion. Science. 2002;298(5594):837–40.
Alexander DL, Mital J, Ward GE, Bradley P, Boothroyd JC. Identification of the moving junction complex of Toxoplasma gondii: a collaboration between distinct secretory organelles. PLoS Pathog. 2005;1(2):e17.
Beck JR, Chen AL, Kim EW, Bradley PJ. RON5 Is Critical for Organization and Function of the Toxoplasma Moving Junction Complex. PLoS Pathog. 2014;10(3):e1004025.
Poukchanski A, Fritz HM, Tonkin ML, Treeck M, Boulanger MJ, Boothroyd JC. Toxoplasma gondii sporozoites invade host cells using two novel paralogues of RON2 and AMA1. PLoS One. 2013;8(8):e70637.
Mital J, Meissner M, Soldati D, Ward GE. Conditional expression of Toxoplasma gondii apical membrane antigen-1 (TgAMA1) demonstrates that TgAMA1 plays a critical role in host cell invasion. Mol Biol Cell. 2005;16(9):4341–9.
Bargieri DY, Andenmatten N, Lagal V, Thiberge S, Whitelaw JA, Tardieux I, et al. Apical membrane antigen 1 mediates apicomplexan parasite attachment but is dispensable for host cell invasion. Nat Commun. 2013;4:2552.
Tyler JS, Boothroyd JC. The C-terminus of Toxoplasma RON2 provides the crucial link between AMA1 and the host-associated invasion complex. PLoS Pathog. 2011;7(2):e1001282.
Straub KW, Peng ED, Hajagos BE, Tyler JS, Bradley PJ. The moving junction protein RON8 facilitates firm attachment and host cell invasion in Toxoplasma gondii. PLoS Pathog. 2011;7(3):e1002007.
Lamarque M, Besteiro S, Papoin J, Roques M, Vulliez-Le Normand B, Morlon-Guyot J, et al. The RON2-AMA1 interaction is a critical step in moving junction-dependent invasion by apicomplexan parasites. PLoS Pathog. 2011;7(2):e1001276.
Besteiro S, Dubremetz JF, Lebrun M. The moving junction of apicomplexan parasites: a key structure for invasion. Cell Microbiol. 2011;13(6):797–805.
Egarter S, Andenmatten N, Jackson AJ, Whitelaw JA, Pall G, Black JA, et al. The Toxoplasma Acto-MyoA Motor Complex Is Important but Not Essential for Gliding Motility and Host Cell Invasion. PLoS One. 2014;9(3):e91819.
Shaw MK. Cell invasion by Theileria sporozoites. Trends Parasitol. 2003;19(1):2–6.
Hemphill A, Gottstein B, Kaufmann H. Adhesion and invasion of bovine endothelial cells by Neospora caninum. Parasitology. 1996;112(Pt 2):183–97.
Talevich E, Kannan N. Structural and evolutionary adaptation of rhoptry kinases and pseudokinases, a family of coccidian virulence factors. BMC Evol Biol. 2013;13:117.
Jensen KD, Hu K, Whitmarsh RJ, Hassan MA, Julien L, Lu D, et al. Toxoplasma gondii rhoptry 16 kinase promotes host resistance to oral infection and intestinal inflammation only in the context of the dense granule protein GRA15. Infect Immun. 2013;81(6):2156–67.
Yamamoto M, Standley DM, Takashima S, Saiga H, Okuyama M, Kayama H, et al. A single polymorphic amino acid on Toxoplasma gondii kinase ROP16 determines the direct and strain-specific activation of Stat3. J Exp Med. 2009;206(12):2747–60.
Fentress SJ, Behnke MS, Dunay IR, Mashayekhi M, Rommereim LM, Fox BA, et al. Phosphorylation of immunity-related GTPases by a Toxoplasma gondii-secreted kinase promotes macrophage survival and virulence. Cell Host Microbe. 2010;8(6):484–95.
Steinfeldt T, Konen-Waisman S, Tong L, Pawlowski N, Lamkemeyer T, Sibley LD, et al. Phosphorylation of mouse immunity-related GTPase (IRG) resistance proteins is an evasion strategy for virulent Toxoplasma gondii. PLoS Biol. 2010;8(12):e1000576.
Fleckenstein MC, Reese ML, Konen-Waisman S, Boothroyd JC, Howard JC, Steinfeldt T. A Toxoplasma gondii pseudokinase inhibits host IRG resistance proteins. PLoS Biol. 2012;10(7):e1001358.
Okada T, Marmansari D, Li ZM, Adilbish A, Canko S, Ueno A, et al. A novel dense granule protein, GRA22, is involved in regulating parasite egress in Toxoplasma gondii. Mol Biochem Parasitol. 2013;189(1–2):5–13.
Rosowski EE, Lu D, Julien L, Rodda L, Gaiser RA, Jensen KD, et al. Strain-specific activation of the NF-kappaB pathway by GRA15, a novel Toxoplasma gondii dense granule protein. J Exp Med. 2011;208(1):195–212.
Rosowski EE, Saeij JP. Toxoplasma gondii clonal strains all inhibit STAT1 transcriptional activity but polymorphic effectors differentially modulate IFNgamma induced gene expression and STAT1 phosphorylation. PLoS One. 2012;7(12):e51448.
Jensen KD, Wang Y, Wojno ED, Shastri AJ, Hu K, Cornel L, et al. Toxoplasma polymorphic effectors determine macrophage polarization and intestinal inflammation. Cell Host Microbe. 2011;9(6):472–83.
Virreira Winter S, Niedelman W, Jensen KD, Rosowski EE, Julien L, Spooner E, et al. Determinants of GBP recruitment to Toxoplasma gondii vacuoles and the parasitic factors that control it. PLoS One. 2011;6(9):e24434.
Gov L, Karimzadeh A, Ueno N, Lodoen MB. Human innate immunity to Toxoplasma gondii is mediated by host caspase-1 and ASC and parasite GRA15. MBio. 2013;4(4):1–11.
Braun L, Brenier-Pinchart MP, Yogavel M, Curt-Varesano A, Curt-Bertini RL, Hussain T, et al. A Toxoplasma dense granule protein, GRA24, modulates the early immune response to infection by promoting a direct and sustained host p38 MAPK activation. J Exp Med. 2013;210(10):2071–86.
Bougdour A, Durandau E, Brenier-Pinchart MP, Ortet P, Barakat M, Kieffer S, et al. Host cell subversion by Toxoplasma GRA16, an exported dense granule protein that targets the host cell nucleus and alters gene expression. Cell Host Microbe. 2013;13(4):489–500.
Gendrin C, Bittame A, Mercier C, Cesbron-Delauw MF. Post-translational membrane sorting of the Toxoplasma gondii GRA6 protein into the parasite-containing vacuole is driven by its N-terminal domain. Int J Parasitol. 2010;40(11):1325–34.
Mercier C, Adjogble KD, Daubener W, Delauw MF. Dense granules: are they key organelles to help understand the parasitophorous vacuole of all apicomplexa parasites? Int J Parasitol. 2005;35(8):829–49.
Magno RC, Lemgruber L, Vommaro RC, De Souza W, Attias M. Intravacuolar network may act as a mechanical support for Toxoplasma gondii inside the parasitophorous vacuole. Microsc Res Tech. 2005;67(1):45–52.
Coppens I, Dunn JD, Romano JD, Pypaert M, Zhang H, Boothroyd JC, et al. Toxoplasma gondii sequesters lysosomes from mammalian hosts in the vacuolar space. Cell. 2006;125(2):261–74.
Ferguson DJ, Jacobs D, Saman E, Dubremetz JF, Wright SE. In vivo expression and distribution of dense granule protein 7 (GRA7) in the exoenteric (tachyzoite, bradyzoite) and enteric (coccidian) forms of Toxoplasma gondii. Parasitology. 1999;119(Pt 3):259–65.
Morris MT, Carruthers VB. Identification and partial characterization of a second Kazal inhibitor in Toxoplasma gondii. Mol Biochem Parasitol. 2003;128(1):119–22.
Morris MT, Coppin A, Tomavo S, Carruthers VB. Functional analysis of Toxoplasma gondii protease inhibitor 1. J Biol Chem. 2002;277(47):45259–66.
Pszenny V, Davis PH, Zhou XW, Hunter CA, Carruthers VB, Roos DS. Targeted disruption of Toxoplasma gondii serine protease inhibitor 1 increases bradyzoite cyst formation in vitro and parasite tissue burden in mice. Infect Immun. 2012;80(3):1156–65.
Balaji S, Babu MM, Iyer LM, Aravind L. Discovery of the principal specific transcription factors of Apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domains. Nucleic Acids Res. 2005;33(13):3994–4006.
Walker R, Gissot M, Croken MM, Huot L, Hot D, Kim K, et al. The Toxoplasma nuclear factor TgAP2XI-4 controls bradyzoite gene expression and cyst formation. Mol Microbiol. 2013;87(3):641–55.
Radke JB, Lucas O, De Silva EK, Ma Y, Sullivan Jr WJ, Weiss LM, et al. ApiAP2 transcription factor restricts development of the Toxoplasma tissue cyst. Proc Natl Acad Sci U S A. 2013;110(17):6871–6.
Walker R, Gissot M, Huot L, Alayi TD, Hot D, Marot G, et al. Toxoplasma transcription factor TgAP2XI-5 regulates the expression of genes involved in parasite virulence and host invasion. J Biol Chem. 2013;288(43):31127–38.
Behnke MS, Wootton JC, Lehmann MM, Radke JB, Lucas O, Nawas J, et al. Coordinated progression through two subtranscriptomes underlies the tachyzoite cycle of Toxoplasma gondii. PLoS One. 2010;5(8):e12354.
Wang T, Zhou J, Gan X, Wang H, Ding X, Chen L, et al. Toxoplasma gondii induce apoptosis of neural stem cells via endoplasmic reticulum stress pathway. Parasitology. 2014;141(7):988–95.
Saksouk N, Bhatti MM, Kieffer S, Smith AT, Musset K, Garin J, et al. Histone-modifying complexes regulate gene expression pertinent to the differentiation of the protozoan parasite Toxoplasma gondii. Mol Cell Biol. 2005;25(23):10301–14.
Lindner SE, De Silva EK, Keck JL, Llinas M. Structural determinants of DNA binding by a P. falciparum ApiAP2 transcriptional regulator. J Mol Biol. 2010;395(3):558–67.
Mercier C, Lefebvre-Van Hende S, Garber GE, Lecordier L, Capron A, Cesbron-Delauw MF. Common cis-acting elements critical for the expression of several genes of Toxoplasma gondii. Mol Microbiol. 1996;21(2):421–8.
Cowper B, Matthews S, Tomley F. The molecular basis for the distinct host and tissue tropisms of coccidian parasites. Mol Biochem Parasitol. 2012;186(1):1–10.
van Dijk MR, van Schaijk BC, Khan SM, van Dooren MW, Ramesar J, Kaczanowski S, et al. Three members of the 6-cys protein family of Plasmodium play a role in gamete fertility. PLoS Pathog. 2010;6(4):e1000853.
Jurankova J, Opsteegh M, Neumayerova H, Kovarcik K, Frencova A, Balaz V, et al. Quantification of Toxoplasma gondii in tissue samples of experimentally infected goats by magnetic capture and real-time PCR. Vet Parasitol. 2013;193(1–3):95–9.
Basso W, Hartnack S, Pardini L, Maksimov P, Koudela B, Venturini MC, et al. Assessment of diagnostic accuracy of a commercial ELISA for the detection of Toxoplasma gondii infection in pigs compared with IFAT, TgSAG1-ELISA and Western blot, using a Bayesian latent class approach. Int J Parasitol. 2013;43(7):565–70.
Rumble SM, Lacroute P, Dalca AV, Fiume M, Sidow A, Brudno M. SHRiMP: accurate mapping of short color-space reads. PLoS Comput Biol. 2009;5(5):e1000386.
David M, Dzamba M, Lister D, Ilie L, Brudno M. SHRiMP2: sensitive yet practical SHort Read Mapping. Bioinformatics. 2011;27(7):1011–2.
Anders S, McCarthy DJ, Chen Y, Okoniewski M, Smyth GK, Huber W, et al. Count-based differential expression analysis of RNA sequencing data using R and Bioconductor. Nat Protoc. 2013;8(9):1765–86.
We are grateful for the excellent technical assistance by Dr. Manuela Schnyder, Therese Michel and Armin Rüdemann, and the support by the Functional Genomics Center Zurich. This work was supported by fellowships from the University of Zurich to WUB and CL, and in part by the Intramural Research Program of the NIH and NIAID (MEG). NCS is supported by a grant from the Bellberry Foundation. MEG is a scholar of the Canadian Institute for Advanced Research (CIFAR) Program for Integrated Microbial Biodiversity. RAW was supported by a Swiss Government Excellence Scholarship: Postdoctoral Scholarship from the Swiss Confederation.
The authors declare that they have no competing interests.
ABH designed experiments, supervised experimental infections in vivo and in vitro, participated in generation of the biological material for analysis, and in quality controls, performed RNA sequence data analysis, directed and supervised analysis of data, wrote the final draft of the manuscript, and edited the final figures. WUB carried out experimental infections in vivo and in vitro, carried out dissections of host animals, participated in generation of the biological material. CL designed experiments, carried out in vitro infections and dissection of host animals, generated biological material for analysis, performed quality controls, directed and supervised RNA sequencing, participated in analysis of data. CR performed RNA sequence data analysis, directed and supervised analysis of data, designed figures and graphs for the manuscript. MO performed RNA sequence data analysis and quality control. RAW performed RNA sequence data analysis, participated in the preparation of manuscript drafts, designed figures and graphs for the manuscript. MEG performed RNA sequence data analysis participated in the preparation of the final draft of the manuscript, designed figures and graphs for the manuscript. NCS designed experiments, performed RNA sequence data analysis, and participated in the preparation of the final draft of the manuscript. PD designed experiments, supervised experimental infections in vivo and dissections of host animals. All authors read and approved the final manuscript.
Walter U Basso and Christoph Lippuner contributed equally to this work.
Significantly higher expressed genes in merozoites (see also Figure 2). Table S1B. Significantly higher expressed genes in tachyzoites (see also Figure 2). Table S2. List of annotated T. gondii Family A genes. The color key for the expression heat map (log2 fold change) is included. Table S3. Differential expression of KRUF family genes. Table S4. List of definitively and provisionally annotated rhoptry (ROP) genes. The color key for the expression heat map (log2 fold change) is included. Table S5. Protease inhibitors including differentially expressed secreted Kazal-type proteins. The color key for the expression heat map (log2 fold change) is included.
Differential expression of – T. gondii genes and gene families: S1. Differential expression of – T. gondii Family A-E genes. S2. Differential expression of metabolic genes. S3. Differential expression of SRS genes. S4. Differential expression of microneme (MIC) genes annotated in ToxoDB. S5. Differential expression of genes coding for PAN-domain containing proteins. S6. Differential expression of annotated rhoptry (ROP) genes. S7. Differential expression of annotated dense granule (GRA) genes. S8. Differential regulation of TgAP2 genes.
About this article
Cite this article
Hehl, A.B., Basso, W.U., Lippuner, C. et al. Asexual expansion of Toxoplasma gondii merozoites is distinct from tachyzoites and entails expression of non-overlapping gene families to attach, invade, and replicate within feline enterocytes. BMC Genomics 16, 66 (2015). https://doi.org/10.1186/s12864-015-1225-x
- Toxoplasma gondii
- Enteroepithelial development
- Comparative transcriptomics
- Surface antigen
- Stage-specific gene expression