ToxoDB: an integrated Toxoplasma gondii database resource

Yeast three-hybrid screen identifies TgBRADIN/GRA24 as a negative regulator of Toxoplasma gondii bradyzoite differentiation

Differentiation of the protozoan parasite Toxoplasma gondii into its latent bradyzoite stage is a key event in the parasite's life cycle. Compound 2 is an imidazopyridine that was previously shown to inhibit the parasite lytic cycle, in part through inhibition of parasite cGMP-dependent protein kinase. We show here that Compound 2 can also enhance parasite differentiation, and we use yeast three-hybrid analysis to identify TgBRADIN/GRA24 as a parasite protein that interacts directly or indirectly with the compound. Disruption of the TgBRADIN/GRA24 gene leads to enhanced differentiation of the parasite, and the TgBRADIN/GRA24 knockout parasites show decreased susceptibility to the differentiation-enhancing effects of Compound 2. This study represents the first use of yeast three-hybrid analysis to study small-molecule mechanism of action in any pathogenic microorganism, and it identifies a previously unrecognized inhibitor of differentiation in T. gondii. A better understanding of the proteins and mechanisms regulating T. gondii differentiation will enable new approaches to preventing the establishment of chronic infection in this important human pathogen.

Forward genetics screens using macrophages to identify Toxoplasma gondii genes important for resistance to IFN-γ-dependent cell autonomous immunity

Toxoplasma gondii, the causative agent of toxoplasmosis, is an obligate intracellular protozoan pathogen. The parasite invades and replicates within virtually any warm blooded vertebrate cell type. During parasite invasion of a host cell, the parasite creates a parasitophorous vacuole (PV) that originates from the host cell membrane independent of phagocytosis within which the parasite replicates. While IFN-dependent-innate and cell mediated immunity is important for eventual control of infection, innate immune cells, including neutrophils, monocytes and dendritic cells, can also serve as vehicles for systemic dissemination of the parasite early in infection. An approach is described that utilizes the host innate immune response, in this case macrophages, in a forward genetic screen to identify parasite mutants with a fitness defect in infected macrophages following activation but normal invasion and replication in naïve macrophages. Thus, the screen isolates parasite mutants that have a specific defect in their ability to resist the effects of macrophage activation. The paper describes two broad phenotypes of mutant parasites following activation of infected macrophages: parasite stasis versus parasite degradation, often in amorphous vacuoles. The parasite mutants are then analyzed to identify the responsible parasite genes specifically important for resistance to induced mediators of cell autonomous immunity. The paper presents a general approach for the forward genetics screen that, in theory, can be modified to target parasite genes important for resistance to specific antimicrobial mediators. It also describes an approach to evaluate the specific macrophage antimicrobial mediators to which the parasite mutant is susceptible. Activation of infected macrophages can also promote parasite differentiation from the tachyzoite to bradyzoite stage that maintains chronic infection. Therefore, methodology is presented to evaluate the importance of the identified parasite gene to establishment of chronic infection.

Phylogenomic analysis of Emiliania huxleyi provides evidence for haptophyte-stramenopile association and a chimeric haptophyte nuclear genome

Emiliania huxleyi is a haptophyte alga of uncertain phylogenetic affinity containing a secondarily derived, chlorophyll c containing plastid. We sought to characterize its relationships with other taxa by quantifying the bipartitions in which it was included from a group of single protein phylogenetic trees in a way that allowed for variation in taxonomic content and accounted for paralogous sequences. The largest number of sequences supported a phylogenetic relationship of E. huxleyi with the stramenopiles, in particular Aureococcus anophagefferens. Far fewer nuclear sequences gave strong support to the placement of this coccolithophorid with the cryptophyte, Guillardia theta. The majority of the sequences that did support this relationship did not have plastid related functions. These results suggest that the haptophytes may be more closely allied with the heterokonts than with the cryptophytes. Another small set of genes associated E. huxleyi with the Viridiplantae with high support. While these genes could have been acquired with a plastid, the lack of plastid related functions among the proteins for which they code and the lack of other organisms with chlorophyll c containing plastids within these bipartitions suggests other explanations may be possible. This study also identified several genes that may have been transferred from the haptophyte lineage to the dinoflagellates Karenia brevis and Karlodinium veneficum as a result of their haptophyte derived plastid, including some with non-photosynthetic functions.

Improving the gene structure annotation of the apicomplexan parasite Neospora caninum fulfils a vital requirement towards an in silico-derived vaccine

Neospora caninum is an apicomplexan parasite which can cause abortion in cattle, instigating major economic burden. Vaccination has been proposed as the most cost-effective control measure to alleviate this burden. Consequently the overriding aspiration for N. caninum research is the identification and subsequent evaluation of vaccine candidates in animal models. To save time, cost and effort, it is now feasible to use an in silico approach for vaccine candidate prediction. Precise protein sequences, derived from the correct open reading frame, are paramount and arguably the most important factor determining the success or failure of this approach. The challenge is that publicly available N. caninum sequences are mostly derived from gene predictions. Annotated inaccuracies can lead to erroneously predicted vaccine candidates by bioinformatics programs. This study evaluates the current N. caninum annotation for potential inaccuracies. Comparisons with annotation from a closely related pathogen, Toxoplasma gondii, are also made to distinguish patterns of inconsistency. More importantly, a mRNA sequencing (RNA-Seq) experiment is used to validate the annotation. Potential discrepancies originating from a questionable start codon context and exon boundaries were identified in 1943 protein coding sequences. We conclude, where experimental data were available, that the majority of N. caninum gene sequences were reliably predicted. Nevertheless, almost 28% of genes were identified as questionable. Given the limitations of RNA-Seq, the intention of this study was not to replace the existing annotation but to support or oppose particular aspects of it. Ideally, many studies aimed at improving the annotation are required to build a consensus. We believe this study, in providing a new resource on gene structure and annotation, is a worthy contributor to this endeavour.

The structure of bradyzoite-specific enolase from Toxoplasma gondii reveals insights into its dual cytoplasmic and nuclear functions

In addition to catalyzing a central step in glycolysis, enolase assumes a remarkably diverse set of secondary functions in different organisms, including transcription regulation as documented for the oncogene c-Myc promoter-binding protein 1. The apicomplexan parasite Toxoplasma gondii differentially expresses two nuclear-localized, plant-like enolases: enolase 1 (TgENO1) in the latent bradyzoite cyst stage and enolase 2 (TgENO2) in the rapidly replicative tachyzoite stage. A 2.75 Å resolution crystal structure of bradyzoite enolase 1, the second structure to be reported of a bradyzoite-specific protein in Toxoplasma, captures an open conformational state and reveals that distinctive plant-like insertions are located on surface loops. The enolase 1 structure reveals that a unique residue, Glu164, in catalytic loop 2 may account for the lower activity of this cyst-stage isozyme. Recombinant TgENO1 specifically binds to a TTTTCT DNA motif present in the cyst matrix antigen 1 (TgMAG1) gene promoter as demonstrated by gel retardation. Furthermore, direct physical interactions of both nuclear TgENO1 and TgENO2 with the TgMAG1 gene promoter are demonstrated in vivo using chromatin immunoprecipitation (ChIP) assays. Structural and biochemical studies reveal that T. gondii enolase functions are multifaceted, including the coordination of gene regulation in parasitic stage development. Enolase 1 provides a potential lead in the design of drugs against Toxoplasma brain cysts.

Recombinant expression, purification, and crystallization of the glutaminyl-tRNA synthetase from Toxoplasma gondii

Aminoacyl tRNA synthetases play a critical role in protein synthesis by providing precursor transfer-RNA molecules correctly charged with their cognate amino-acids. The essential nature of these enzymes make them attractive targets for designing new drugs against important pathogenic protozoans like Toxoplasma. Because no structural data currently exists for a protozoan glutaminyl-tRNA synthetase (QRS), an understanding of its potential as a drug target and its function in the assembly of the Toxoplasma multi-aminoacyl tRNA (MARS) complex is therefore lacking. Here we describe the optimization of expression and purification conditions that permitted the recovery and crystallization of both domains of the Toxoplasma QRS enzyme from a heterologous Escherichia coli expression system. Expression of full-length QRS was only achieved after the addition of an N-terminal histidine affinity tag and the isolated protein was active on both cellular and in vitro produced Toxoplasma tRNA. Taking advantage of the proteolytic susceptibility of QRS to cleavage into component domains, N-terminal glutathione S-transferase (GST) motif-containing domain fragments were isolated and crystallization conditions discovered. Isolation of the C-terminal catalytic domain was accomplished after subcloning the domain and optimizing expression conditions. Purified catalytic domain survived cryogenic storage and yielded large diffraction-quality crystals over-night after optimization of screening conditions. This work will form the basis of future structural studies into structural-functional relationships of both domains including potential targeted drug-design studies and investigations into the assembly of the Toxoplasma MARS complex.

Secreted Toxoplasma gondii molecules interfere with expression of MHC-II in interferon gamma-activated macrophages

The obligate intracellular protozoan parasite Toxoplasma gondii interferes with major histocompatibility complex class II antigen presentation to dampen host CD4(+) T cell responses. While it is known that T. gondii inhibits major histocompatibility complex class II gene transcription and expression in infected host cells, the mechanism of this host manipulation is unknown. Here, we show that soluble parasite proteins inhibit IFNγ-induced expression of major histocompatibility complex class II on the surface of the infected cell in a dose-dependent response that was abolished by protease treatment. Subcellular fractionation of T. gondii tachyzoites revealed that the major histocompatibility complex class II inhibitory activity co-partitioned with rhoptries and/or dense granules. However, parasite mutants deleted for single rhoptries or dense granules genes (ROP1, 4/7, 14, 16 and 18 or GRA 2-9 and 12 knock-out strains) retained the ability to inhibit expression of major histocompatibility complex class II. In addition, excreted/secreted antigens released by extracellular tachyzoites displayed immunomodulatory activity characterized by an inhibition of major histocompatibility complex class II expression, and reduced expression and release of TNFα by macrophages. Tandem MS analysis of parasite excreted/secreted antigens generated a list of T. gondii secreted proteins that may participate in major histocompatibility complex class II inhibition and the modulation of host immune functions.

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

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

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

Widespread occurrence of organelle genome-encoded 5S rRNAs including permuted molecules

5S Ribosomal RNA (5S rRNA) is a universal component of ribosomes, and the corresponding gene is easily identified in archaeal, bacterial and nuclear genome sequences. However, organelle gene homologs (rrn5) appear to be absent from most mitochondrial and several chloroplast genomes. Here, we re-examine the distribution of organelle rrn5 by building mitochondrion- and plastid-specific covariance models (CMs) with which we screened organelle genome sequences. We not only recover all organelle rrn5 genes annotated in GenBank records, but also identify more than 50 previously unrecognized homologs in mitochondrial genomes of various stramenopiles, red algae, cryptomonads, malawimonads and apusozoans, and surprisingly, in the apicoplast (highly derived plastid) genomes of the coccidian pathogens Toxoplasma gondii and Eimeria tenella. Comparative modeling of RNA secondary structure reveals that mitochondrial 5S rRNAs from brown algae adopt a permuted triskelion shape that has not been seen elsewhere. Expression of the newly predicted rrn5 genes is confirmed experimentally in 10 instances, based on our own and published RNA-Seq data. This study establishes that particularly mitochondrial 5S rRNA has a much broader taxonomic distribution and a much larger structural variability than previously thought. The newly developed CMs will be made available via the Rfam database and the MFannot organelle genome annotator.

Genetic mapping reveals that sinefungin resistance in Toxoplasma gondii is controlled by a putative amino acid transporter locus that can be used as a negative selectable marker

Quantitative trait locus (QTL) mapping studies have been integral in identifying and understanding virulence mechanisms in the parasite Toxoplasma gondii. In this study, we interrogated a different phenotype by mapping sinefungin (SNF) resistance in the genetic cross between type 2 ME49-FUDR(r) and type 10 VAND-SNF(r). The genetic map of this cross was generated by whole-genome sequencing of the progeny and subsequent identification of single nucleotide polymorphisms (SNPs) inherited from the parents. Based on this high-density genetic map, we were able to pinpoint the sinefungin resistance phenotype to one significant locus on chromosome IX. Within this locus, a single nonsynonymous SNP (nsSNP) resulting in an early stop codon in the TGVAND_290860 gene was identified, occurring only in the sinefungin-resistant progeny. Using CRISPR/CAS9, we were able to confirm that targeted disruption of TGVAND_290860 renders parasites sinefungin resistant. Because disruption of the SNR1 gene confers resistance, we also show that it can be used as a negative selectable marker to insert either a positive drug selection cassette or a heterologous reporter. These data demonstrate the power of combining classical genetic mapping, whole-genome sequencing, and CRISPR-mediated gene disruption for combined forward and reverse genetic strategies in T. gondii.

DB-AT: a 2015 update to the Full-parasites database brings a multitude of new transcriptomic data for apicomplexan parasites

The previous release of our Full-parasites database (http://fullmal.hgc.jp/) brought enhanced functionality, an expanded full-length cDNA content, and new RNA-Seq datasets from several important apicomplexan parasites. The 2015 update witnesses the major shift in the databases content with focus on diverse transcriptomes of the apicomplexan parasites. The content of the database was substantially enriched with transcriptome information for new apicomplexan parasites. The latest version covers a total of 17 species, with addition of our newly generated RNA-Seq data of a total of 909 150 388 tags. Moreover, we have generated and included two novel and unique datasets, which represent diverse nature of transcriptomes in individual parasites in vivo and in vitro. One is the data collected from 116 Indonesian patients infected with Plasmodium falciparum. The other is a series of transcriptome data collected from a total of 38 single cells of P. falciparum cultured in vitro. We believe that with the recent advances our database becomes an even better resource and a unique platform in the analysis of apicomplexan parasites and their interaction with their hosts. To adequately reflect the recent modifications and the current content we have changed the database name to DB-AT—DataBase of Apicomplexa Transcriptomes.

Systematic identification of the lysine succinylation in the protozoan parasite Toxoplasma gondii

Lysine succinylation is a new posttranslational modification identified in histone proteins of Toxoplasma gondii, an obligate intracellular parasite of the phylum Apicomplexa. However, very little is known about their scope and cellular distribution. Here, using LC-MS/MS to identify parasite peptides enriched by immunopurification with succinyl lysine antibody, we produced the first lysine succinylome in this parasite. Overall, a total of 425 lysine succinylation sites that occurred on 147 succinylated proteins were identified in extracellular Toxoplasma tachyzoites, which is a proliferative stage that results in acute toxoplasmosis. With the bioinformatics analysis, it is shown that these succinylated proteins are evolutionarily conserved and involved in a wide variety of cellular functions such as metabolism and epigenetic gene regulation and exhibit diverse subcellular localizations. Moreover, we defined five types of definitively conserved succinylation site motifs, and the results imply that lysine residue of a polypeptide with lysine on the +3 position and without lysine at the -1 to +2 position is a preferred substrate of lysine succinyltransferase. In conclusion, our findings suggest that lysine succinylation in Toxoplasma involves a diverse array of cellular functions, although the succinylation occurs at a low level.

Plasticity between MyoC- and MyoA-glideosomes: an example of functional compensation in Toxoplasma gondii invasion

The glideosome is an actomyosin-based machinery that powers motility in Apicomplexa and participates in host cell invasion and egress from infected cells. The central component of the glideosome, myosin A (MyoA), is a motor recruited at the pellicle by the acylated gliding-associated protein GAP45. In Toxoplasma gondii, GAP45 also contributes to the cohesion of the pellicle, composed of the inner membrane complex (IMC) and the plasma membrane, during motor traction. GAP70 was previously identified as a paralog of GAP45 that is tailored to recruit MyoA at the apical cap in the coccidian subgroup of the Apicomplexa. A third member of this family, GAP80, is demonstrated here to assemble a new glideosome, which recruits the class XIV myosin C (MyoC) at the basal polar ring. MyoC shares the same myosin light chains as MyoA and also interacts with the integral IMC proteins GAP50 and GAP40. Moreover, a central component of this complex, the IMC-associated protein 1 (IAP1), acts as the key determinant for the restricted localization of MyoC to the posterior pole. Deletion of specific components of the MyoC-glideosome underscores the installation of compensatory mechanisms with components of the MyoA-glideosome. Conversely, removal of MyoA leads to the relocalization of MyoC along the pellicle and at the apical cap that accounts for residual invasion. The two glideosomes exhibit a considerable level of plasticity to ensure parasite survival.

Toxoplasma gondii can infect most warm-blooded animals, and is an important opportunistic pathogen for humans. This obligate intracellular parasite is able to invade virtually all nucleated cells, and as with most parasites of the Apicomplexa phylum, relies on a substrate-dependent gliding motility to actively penetrate into host cells and egress from infected cells. The conserved molecular machine (named glideosome) powering motility is located at the periphery of the parasite and involves the molecular motor, myosin A (MyoA). The glideosome exists in three flavors, exhibiting the same overall organization and sharing some common components while being spatially restricted to the central IMC, the apical cap and the basal pole of the parasite, respectively. The central and apical glideosomes are associated with MyoA (MyoA-glideosome) whereas the basal complex recruits myosin C (MyoC). Deleting components of the MyoC-glideosome uncovers the existence of complementary and compensatory mechanisms that ensure successful establishment of infection. This study highlights a higher degree of complexity and plasticity of the gliding machinery.

Structural and functional divergence of the aldolase fold in Toxoplasma gondii

Parasites of the phylum Apicomplexa are highly successful pathogens of humans and animals worldwide. As obligate intracellular parasites, they have significant energy requirements for invasion and gliding motility that are supplied by various metabolic pathways. Aldolases have emerged as key enzymes involved in these pathways, and all apicomplexans express one or both of fructose 1,6-bisphosphate (F16BP) aldolase and 2-deoxyribose 5-phosphate (dR5P) aldolase (DERA). Intriguingly, Toxoplasma gondii, a highly successful apicomplexan parasite, expresses F16BP aldolase (TgALD1), d5RP aldolase (TgDERA), and a divergent dR5P aldolase-like protein (TgDPA) exclusively in the latent bradyzoite stage. While the importance of TgALD1 in glycolysis is well established and TgDERA is also likely to be involved in parasite metabolism, the detailed function of TgDPA remains elusive. To gain mechanistic insight into the function of different T. gondii aldolases, we first determined the crystal structures of TgALD1 and TgDPA. Structural analysis revealed that both aldolases adopt a TIM barrel fold accessorized with divergent secondary structure elements. Structural comparison of TgALD1 and TgDPA with members of their respective enzyme families revealed that, while the active-site residues are conserved in TgALD1, key catalytic residues are absent in TgDPA. Consistent with this observation, biochemical assays showed that, while TgALD1 was active on F16BP, TgDPA was inactive on dR5P. Intriguingly, both aldolases are competent to bind polymerized actin in vitro. Altogether, structural and biochemical analyses of T. gondii aldolase and aldolase-like proteins reveal diverse functionalization of the classic TIM barrel aldolase fold.

Large, rapidly evolving gene families are at the forefront of host-parasite interactions in Apicomplexa

The Apicomplexa is a phylum of parasitic protozoa, which includes the malaria parasite Plasmodium, amongst other species that can devastate human and animal health. The past decade has seen the release of genome sequences for many of the most important apicomplexan species, providing an excellent basis for improving our understanding of their biology. One of the key features of each genome is a unique set of large, variant gene families. Although closely related species share the same families, even different types of malaria parasite have distinct families. In some species they tend to be found at the ends of chromosomes, which may facilitate aspects of gene expression regulation and generation of sequence diversity. In others they are scattered apparently randomly across chromosomes. For some families there is evidence they are involved in antigenic variation, immune regulation and immune evasion. For others there are no known functions. Even where function is unknown these families are most often predicted to be exposed to the host, contain much sequence diversity and evolve rapidly. Based on these properties it is clear that they are at the forefront of host–parasite interactions. In this review I compare and contrast the genomic context, gene structure, gene expression, protein localization and function of these families across different species.

Proteomics reveals differences in protein abundance and highly similar antigenic profiles between Besnoitia besnoiti and Besnoitia tarandi

Besnoitia besnoiti and Besnoitia tarandi are two cyst-forming apicomplexan parasites of the genus Besnoitia. B. besnoiti uses cattle as an intermediate host, in which it causes a disease that progresses in two sequential phases: the acute anasarca stage and the chronic scleroderma stage. Reindeer and caribou act as intermediate hosts for B. tarandi, which causes clinical signs similar to those caused by B. besnoiti. Previous studies demonstrated high molecular similarity, as determined by 18S and ITS-1 RNA sequences, between these Besnoitia spp., and strong serological cross-reactivity between these species has recently been demonstrated. Thus, a difference gel electrophoresis approach and mass spectrometry analysis were used to describe the proteomes and explore differences in protein abundance between B. besnoiti and B. tarandi in tachyzoite extracts. Immunoproteomes were also compared using 2-DE immunoblotting with polyclonal sera from experimentally infected rabbits. From approximately 1400 spots detected in DIGE-gels, 28 and 29 spots were differentially abundant in B. besnoiti and B. tarandi tachyzoites, respectively (± 1.5-fold, p<0.05). Four and 13 spots were exclusively detected in B. besnoiti and B. tarandi, respectively. Of the 32 differentially abundant spots analyzed by MALDI-TOF/MS, 6 up-regulated B. besnoiti proteins (LDH; HSP90; purine nucleoside phosphorylase and 3 hypothetical proteins) and 6 up-regulated B. tarandi proteins (G3PDH; LDH; PDI; mRNA decapping protein and 2 hypothetical proteins) were identified. Interestingly, no specific antigen spots were recognized by sera on any of the Besnoitia species studied and a similar antigen profile has been observed for B. tarandi and B. besnoiti sera when cross reactions were studied. This fact corroborates the difficulty in discerning Besnoitia infections using current serological assays. The present study underscores the importance of sequencing the B. besnoiti genome for species diversity studies of the genus Besnoitia.

kdetrees: Non-parametric estimation of phylogenetic tree distributions

MOTIVATION

Although the majority of gene histories found in a clade of organisms are expected to be generated by a common process (e.g. the coalescent process), it is well known that numerous other coexisting processes (e.g. horizontal gene transfers, gene duplication and subsequent neofunctionalization) will cause some genes to exhibit a history distinct from those of the majority of genes. Such 'outlying' gene trees are considered to be biologically interesting, and identifying these genes has become an important problem in phylogenetics.

RESULTS

We propose and implement kdetrees, a non-parametric method for estimating distributions of phylogenetic trees, with the goal of identifying trees that are significantly different from the rest of the trees in the sample. Our method compares favorably with a similar recently published method, featuring an improvement of one polynomial order of computational complexity (to quadratic in the number of trees analyzed), with simulation studies suggesting only a small penalty to classification accuracy. Application of kdetrees to a set of Apicomplexa genes identified several unreliable sequence alignments that had escaped previous detection, as well as a gene independently reported as a possible case of horizontal gene transfer. We also analyze a set of Epichloë genes, fungi symbiotic with grasses, successfully identifying a contrived instance of paralogy.

AVAILABILITY AND IMPLEMENTATION

Our method for estimating tree distributions and identifying outlying trees is implemented as the R package kdetrees and is available for download from CRAN.

Horizontal gene transfer is a significant driver of gene innovation in dinoflagellates

The dinoflagellates are an evolutionarily and ecologically important group of microbial eukaryotes. Previous work suggests that horizontal gene transfer (HGT) is an important source of gene innovation in these organisms. However, dinoflagellate genomes are notoriously large and complex, making genomic investigation of this phenomenon impractical with currently available sequencing technology. Fortunately, de novo transcriptome sequencing and assembly provides an alternative approach for investigating HGT. We sequenced the transcriptome of the dinoflagellate Alexandrium tamarense Group IV to investigate how HGT has contributed to gene innovation in this group. Our comprehensive A. tamarense Group IV gene set was compared with those of 16 other eukaryotic genomes. Ancestral gene content reconstruction of ortholog groups shows that A. tamarense Group IV has the largest number of gene families gained (314-1,563 depending on inference method) relative to all other organisms in the analysis (0-782). Phylogenomic analysis indicates that genes horizontally acquired from bacteria are a significant proportion of this gene influx, as are genes transferred from other eukaryotes either through HGT or endosymbiosis. The dinoflagellates also display curious cases of gene loss associated with mitochondrial metabolism including the entire Complex I of oxidative phosphorylation. Some of these missing genes have been functionally replaced by bacterial and eukaryotic xenologs. The transcriptome of A. tamarense Group IV lends strong support to a growing body of evidence that dinoflagellate genomes are extraordinarily impacted by HGT.

Auranofin is highly efficacious against Toxoplasma gondii in vitro and in an in vivo experimental model of acute toxoplasmosis

BACKGROUND

The mainstay of toxoplasmosis treatment targets the folate biosynthetic pathways and has not changed for the last 50 years. The activity of these chemotherapeutic agents is restricted to one lifecycle stage of Toxoplasma gondii, they have significant toxicity, and the impending threat of emerging resistance to these agents makes the discovery of new therapies a priority. We now demonstrate that auranofin, an orally administered gold containing compound that was FDA approved for treatment of rheumatoid arthritis, has activity against Toxoplasma gondii in vitro (IC50 = 0.28 µM) and in vivo (1 mg/kg).

METHODS/PRINCIPAL FINDINGS

Replication within human foreskin fibroblasts of RH tachyzoites was inhibited by auranofin. At 0.4 µM, auranofin inhibited replication, as measured by percent infected fibroblasts at 24 hrs, (10.94% vs. 24.66% of controls; p = 0.0003) with no effect on parasite invasion (16.95% vs. 12.91% p = 0.4331). After 18 hrs, 62% of extracellular parasites treated with auranofin were non-viable compared to control using an ATP viability assay (p = 0.0003). In vivo, a previously standardized chicken embryo model of acute toxoplasmosis was used. Fourteen day old chicken embryos were injected through the chorioallantoic vein with 1×104 tachyzoites of the virulent RH strain. The treatment group received one dose of auranofin at the time of inoculation (1 mg/kg estimated body weight). On day 5, auranofin-treated chicken embryos were 100% protected against death (p = 0.0002) and had a significantly reduced parasite load as determined by histopathology, immunohistochemistry and by the number of parasites quantified by real-time PCR.

CONCLUSIONS

These results reveal in vitro and in vivo activity of auranofin against T. gondii, suggesting that it may be an effective alternative treatment for toxoplasmosis.

Multisystemic toxoplasmosis associated with a type II-like Toxoplasma gondii strain in a New Zealand fur seal (Arctocephalus forsteri) from New South Wales, Australia

We report the first confirmed case of toxoplasmosis in an Australian pinniped. Presence of Toxoplasma gondii DNA was detected in the brain of a free-ranging subadult New Zealand fur seal (Arctocephalus forsteri) with nonsuppurative meningoencephalitis, hypophysitis, posterior uveitis, retrobulbar cellulitis, and myocarditis associated with protozoan cysts and tachyzoites. The emaciated seal stranded moribund on a beach in northern Sydney in New South Wales. Histopathology coupled with specific immunohistochemistry and PCR assays confirmed the presence of T. gondii. The T. gondii sample (NZfs8825) identified in this study has an identical genotype as the type II (ToxoDB PCR-RFLP genotype #1) based on the direct sequencing and virtual RFLP of multilocus DNA markers including SAG1, 5'- and 3'-SAG2, alt.SAG2, SAG3, BTUB, GRA6, c22-8, c29-2, L358, PK1, and Apico. Direct sequencing of T. gondii B1 DNA marker from the T. gondii sample (NZfs8825) identified a type II-like strain, based on presence of non-archetypal B1 gene polymorphisms previously reported as unique to Australia. This study suggests that T. gondii oocysts originating from mainland Australia, which has a large population of feral cats, may act as a disease threat to native marine fauna. Therefore, emerging toxoplasmosis in the Arctic has a relevant parallel in the Southern Ocean within Australian waters with yet unknown relevance to Antarctica.

An atomic-resolution view of neofunctionalization in the evolution of apicomplexan lactate dehydrogenases

Malate and lactate dehydrogenases (MDH and LDH) are homologous, core metabolic enzymes that share a fold and catalytic mechanism yet possess strict specificity for their substrates. In the Apicomplexa, convergent evolution of an unusual LDH from MDH produced a difference in specificity exceeding 12 orders of magnitude. The mechanisms responsible for this extraordinary functional shift are currently unknown. Using ancestral protein resurrection, we find that specificity evolved in apicomplexan LDHs by classic neofunctionalization characterized by long-range epistasis, a promiscuous intermediate, and few gain-of-function mutations of large effect. In canonical MDHs and LDHs, a single residue in the active-site loop governs substrate specificity: Arg102 in MDHs and Gln102 in LDHs. During the evolution of the apicomplexan LDH, however, specificity switched via an insertion that shifted the position and identity of this 'specificity residue' to Trp107f. Residues far from the active site also determine specificity, as shown by the crystal structures of three ancestral proteins bracketing the key duplication event. This work provides an unprecedented atomic-resolution view of evolutionary trajectories creating a nascent enzymatic function.

Gene Set Enrichment Analysis (GSEA) of Toxoplasma gondii expression datasets links cell cycle progression and the bradyzoite developmental program

Background

Large amounts of microarray expression data have been generated for the Apicomplexan parasite Toxoplasma gondii in an effort to identify genes critical for virulence or developmental transitions. However, researchers’ ability to analyze this data is limited by the large number of unannotated genes, including many that appear to be conserved hypothetical proteins restricted to Apicomplexa. Further, differential expression of individual genes is not always informative and often relies on investigators to draw big-picture inferences without the benefit of context. We hypothesized that customization of gene set enrichment analysis (GSEA) to T. gondii would enable us to rigorously test whether groups of genes serving a common biological function are co-regulated during the developmental transition to the latent bradyzoite form.

Results

Using publicly available T. gondii expression microarray data, we created Toxoplasma gene sets related to bradyzoite differentiation, oocyst sporulation, and the cell cycle. We supplemented these with lists of genes derived from community annotation efforts that identified contents of the parasite-specific organelles, rhoptries, micronemes, dense granules, and the apicoplast. Finally, we created gene sets based on metabolic pathways annotated in the KEGG database and Gene Ontology terms associated with gene annotations available at http://www.toxodb.org. These gene sets were used to perform GSEA analysis using two sets of published T. gondii expression data that characterized T. gondii stress response and differentiation to the latent bradyzoite form.

Conclusions

GSEA provides evidence that cell cycle regulation and bradyzoite differentiation are coupled. Δgcn5A mutants unable to induce bradyzoite-associated genes in response to alkaline stress have different patterns of cell cycle and bradyzoite gene expression from stressed wild-type parasites. Extracellular tachyzoites resemble a transitional state that differs in gene expression from both replicating intracellular tachyzoites and in vitro bradyzoites by expressing genes that are enriched in bradyzoites as well as genes that are associated with the G1 phase of the cell cycle. The gene sets we have created are readily modified to reflect ongoing research and will aid researchers’ ability to use a knowledge-based approach to data analysis facilitating the development of new insights into the intricate biology of Toxoplasma gondii.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-515) contains supplementary material, which is available to authorized users.

A survey of innovation through duplication in the reduced genomes of twelve parasites

We characterize the prevalence, distribution, divergence, and putative functions of detectable two-copy paralogs and segmental duplications in the Apicomplexa, a phylum of parasitic protists. Apicomplexans are mostly obligate intracellular parasites responsible for human and animal diseases (e.g. malaria and toxoplasmosis). Gene loss is a major force in the phylum. Genomes are small and protein-encoding gene repertoires are reduced. Despite this genomic streamlining, duplications and gene family amplifications are present. The potential for innovation introduced by duplications is of particular interest. We compared genomes of twelve apicomplexans across four lineages and used orthology and genome cartography to map distributions of duplications against genome architectures. Segmental duplications appear limited to five species. Where present, they correspond to regions enriched for multi-copy and species-specific genes, pointing toward roles in adaptation and innovation. We found a phylum-wide association of duplications with dynamic chromosome regions and syntenic breakpoints. Trends in the distribution of duplicated genes indicate that recent, species-specific duplicates are often tandem while most others have been dispersed by genome rearrangements. These trends show a relationship between genome architecture and gene duplication. Functional analysis reveals: proteases, which are vital to a parasitic lifecycle, to be prominent in putative recent duplications; a pair of paralogous genes in Toxoplasma gondii previously shown to produce the rate-limiting step in dopamine synthesis in mammalian cells, a possible link to the modification of host behavior; and phylum-wide differences in expression and subcellular localization, indicative of modes of divergence. We have uncovered trends in multiple modes of duplicate divergence including sequence, intron content, expression, subcellular localization, and functions of putative recent duplicates that highlight the role of duplications in the continuum of forces that have shaped these genomes.

Identification of T. gondii myosin light chain-1 as a direct target of TachypleginA-2, a small-molecule inhibitor of parasite motility and invasion

Motility of the protozoan parasite Toxoplasma gondii plays an important role in the parasite’s life cycle and virulence within animal and human hosts. Motility is driven by a myosin motor complex that is highly conserved across the Phylum Apicomplexa. Two key components of this complex are the class XIV unconventional myosin, TgMyoA, and its associated light chain, TgMLC1. We previously showed that treatment of parasites with a small-molecule inhibitor of T. gondii invasion and motility, tachypleginA, induces an electrophoretic mobility shift of TgMLC1 that is associated with decreased myosin motor activity. However, the direct target(s) of tachypleginA and the molecular basis of the compound-induced TgMLC1 modification were unknown. We show here by “click” chemistry labelling that TgMLC1 is a direct and covalent target of an alkyne-derivatized analogue of tachypleginA. We also show that this analogue can covalently bind to model thiol substrates. The electrophoretic mobility shift induced by another structural analogue, tachypleginA-2, was associated with the formation of a 225.118 Da adduct on S57 and/or C58, and treatment with deuterated tachypleginA-2 confirmed that the adduct was derived from the compound itself. Recombinant TgMLC1 containing a C58S mutation (but not S57A) was refractory to click labelling and no longer exhibited a mobility shift in response to compound treatment, identifying C58 as the site of compound binding on TgMLC1. Finally, a knock-in parasite line expressing the C58S mutation showed decreased sensitivity to compound treatment in a quantitative 3D motility assay. These data strongly support a model in which tachypleginA and its analogues inhibit the motility of T. gondii by binding directly and covalently to C58 of TgMLC1, thereby causing a decrease in the activity of the parasite’s myosin motor.

Whole genome profiling of spontaneous and chemically induced mutations in Toxoplasma gondii

Background

Next generation sequencing is helping to overcome limitations in organisms less accessible to classical or reverse genetic methods by facilitating whole genome mutational analysis studies. One traditionally intractable group, the Apicomplexa, contains several important pathogenic protozoan parasites, including the Plasmodium species that cause malaria.

Here we apply whole genome analysis methods to the relatively accessible model apicomplexan, Toxoplasma gondii, to optimize forward genetic methods for chemical mutagenesis using N-ethyl-N-nitrosourea (ENU) and ethylmethane sulfonate (EMS) at varying dosages.

Results

By comparing three different lab-strains we show that spontaneously generated mutations reflect genome composition, without nucleotide bias. However, the single nucleotide variations (SNVs) are not distributed randomly over the genome; most of these mutations reside either in non-coding sequence or are silent with respect to protein coding. This is in contrast to the random genomic distribution of mutations induced by chemical mutagenesis. Additionally, we report a genome wide transition vs transversion ratio (ti/tv) of 0.91 for spontaneous mutations in Toxoplasma, with a slightly higher rate of 1.20 and 1.06 for variants induced by ENU and EMS respectively. We also show that in the Toxoplasma system, surprisingly, both ENU and EMS have a proclivity for inducing mutations at A/T base pairs (78.6% and 69.6%, respectively).

Conclusions

The number of SNVs between related laboratory strains is relatively low and managed by purifying selection away from changes to amino acid sequence. From an experimental mutagenesis point of view, both ENU (24.7%) and EMS (29.1%) are more likely to generate variation within exons than would naturally accumulate over time in culture (19.1%), demonstrating the utility of these approaches for yielding proportionally greater changes to the amino acid sequence. These results will not only direct the methods of future chemical mutagenesis in Toxoplasma, but also aid in designing forward genetic approaches in less accessible pathogenic protozoa as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-354) contains supplementary material, which is available to authorized users.

Toxoplasma gondii merozoite gene expression analysis with comparison to the life cycle discloses a unique expression state during enteric development

BACKGROUND

Considerable work has been carried out to understand the biology of tachyzoites and bradyzoites of Toxoplasma gondii in large part due to in vitro culture methods for these stages. However, culturing methods for stages that normally develop in the gut of the definitive felid host, including the merozoite and sexual stages, have not been developed hindering the ability to study a large portion of the parasite's life cycle. Here, we begin to unravel the molecular aspects of enteric stages by providing new data on merozoite stage gene expression.

RESULTS

To profile gene expression differences in enteric stages we harvested merozoites from the intestine of infected cats and hybridized mRNA to the Affymetrix Toxoplasma GeneChip. We analyzed the merozoite data in context of the life cycle by comparing it to previously published data for the oocyst, tachyzoite, and bradyzoite stages. Principal component analysis highlighted the unique profile of merozoites, placing them approximately half-way on a continuum between the tachyzoite/bradyzoite and oocyst samples. Prior studies have shown that antibodies to surface antigen one (SAG1) and many dense granule proteins do not label merozoites: our microarray data confirms that these genes were not expressed at this stage. Also, the expression for many rhoptry and microneme proteins was drastically reduced while the expression for many surface antigens was increased at the merozoite stage. Gene Ontology and KEGG analysis revealed that genes involved in transcription/translation and many metabolic pathways were upregulated at the merozoite stage, highlighting unique growth requirements of this stage. To functionally test these predictions, we demonstrated that an upstream promoter region of a merozoite specific gene was sufficient to control expression in merozoites in vivo.

CONCLUSIONS

Merozoites are the first developmental stage in the coccidian cycle that takes place within the gut of the definitive host. The data presented here describe the global gene expression profile of the merozoite stage and the creation of transgenic parasite strains that show stage-specific expression of reporter genes in the cat intestine. These data and reagents will be useful in unlocking how the parasite senses and responds to the felid gut environment to initiate enteric development.

The apical complex provides a regulated gateway for secretion of invasion factors in Toxoplasma

The apical complex is the definitive cell structure of phylum Apicomplexa, and is the focus of the events of host cell penetration and the establishment of intracellular parasitism. Despite the importance of this structure, its molecular composition is relatively poorly known and few studies have experimentally tested its functions. We have characterized a novel Toxoplasma gondii protein, RNG2, that is located at the apical polar ring--the common structural element of apical complexes. During cell division, RNG2 is first recruited to centrosomes immediately after their duplication, confirming that assembly of the new apical complex commences as one of the earliest events of cell replication. RNG2 subsequently forms a ring, with the carboxy- and amino-termini anchored to the apical polar ring and mobile conoid, respectively, linking these two structures. Super-resolution microscopy resolves these two termini, and reveals that RNG2 orientation flips during invasion when the conoid is extruded. Inducible knockdown of RNG2 strongly inhibits host cell invasion. Consistent with this, secretion of micronemes is prevented in the absence of RNG2. This block, however, can be fully or partially overcome by exogenous stimulation of calcium or cGMP signaling pathways, respectively, implicating the apical complex directly in these signaling events. RNG2 demonstrates for the first time a role for the apical complex in controlling secretion of invasion factors in this important group of parasites.

Characterization of a homolog of DEAD-box RNA helicases in Toxoplasma gondii as a marker of cytoplasmic mRNP stress granules

Toxoplasma gondii is an obligate intracellular protozoan which infects one-third of the human population. Due to its high infection prevalence, Toxoplasma offers an ideal system for the study of host-parasite interaction. Similar to other eukaryotes, Toxoplasma maintains levels and localization of cytoplasmic mRNAs throughout its life cycle as part of a gene regulation network to meet all cellular and biochemical requirements. More recently, it was reported that the presence of cytoplasmic mRNA granules could contribute to the parasite pathogenesis and viability. Here we identified a novel Toxoplasma DEAD-box RNA helicase, referred to as Toxoplasma gondiiHomolog of DOZI (TgHoDI), because of its high homology (81%) to Plasmodium DOZI. TgHoDI is the functional ortholog of yeast DHH1, and its function was authenticated by complementation studies in Δdhh1 yeast strain. We demonstrated that TgHoDI is a marker of cytoplasmic RNA stress granules, which assemble when the parasites experience cellular stresses and translational arrest.

Assembly of the novel five-component apicomplexan multi-aminoacyl-tRNA synthetase complex is driven by the hybrid scaffold protein Tg-p43

In Toxoplasma gondii, as in other eukaryotes, a subset of the amino-acyl-tRNA synthetases are arranged into an abundant cytoplasmic multi-aminoacyl-tRNA synthetase (MARS) complex. Through a series of genetic pull-down assays, we have identified the enzymes of this complex as: methionyl-, glutaminyl-, glutamyl-, and tyrosyl-tRNA synthetases, and we show that the N-terminal GST-like domain of a partially disordered hybrid scaffold protein, Tg-p43, is sufficient for assembly of the intact complex. Our gel filtration studies revealed significant heterogeneity in the size and composition of isolated MARS complexes. By targeting the tyrosyl-tRNA synthetases subunit, which was found exclusively in the complete 1 MDa complex, we were able to directly visualize MARS particles in the electron microscope. Image analyses of the negative stain data revealed the observed heterogeneity and instability of these complexes to be driven by the intrinsic flexibility of the domain arrangements within the MARS complex. These studies provide unique insights into the assembly of these ubiquitous but poorly understood eukaryotic complexes.

Molecular cloning and characterization of NcROP2Fam-1, a member of the ROP2 family of rhoptry proteins in Neospora caninum that is targeted by antibodies neutralizing host cell invasion in vitro

Recent publications demonstrated that a fragment of a Neospora caninum ROP2 family member antigen represents a promising vaccine candidate. We here report on the cloning of the cDNA encoding this protein, N. caninum ROP2 family member 1 (NcROP2Fam-1), its molecular characterization and localization. The protein possesses the hallmarks of ROP2 family members and is apparently devoid of catalytic activity. NcROP2Fam-1 is synthesized as a pre-pro-protein that is matured to 2 proteins of 49 and 55 kDa that localize to rhoptry bulbs. Upon invasion the protein is associated with the nascent parasitophorous vacuole membrane (PVM), evacuoles surrounding the host cell nucleus and, in some instances, the surface of intracellular parasites. Staining was also observed within the cyst wall of 'cysts' produced in vitro. Interestingly, NcROP2Fam-1 was also detected on the surface of extracellular parasites entering the host cells and antibodies directed against NcROP2Fam-1-specific peptides partially neutralized invasion in vitro. We conclude that, in spite of the general belief that ROP2 family proteins are intracellular antigens, NcROP2Fam-1 can also be considered as an extracellular antigen, a property that should be taken into account in further experiments employing ROP2 family proteins as vaccines.

Phosphatidylethanolamine synthesis in the parasite mitochondrion is required for efficient growth but dispensable for survival of Toxoplasma gondii

Toxoplasma gondii is a highly prevalent obligate intracellular parasite of the phylum Apicomplexa, which also includes other parasites of clinical and/or veterinary importance, such as Plasmodium, Cryptosporidium, and Eimeria. Acute infection by Toxoplasma is hallmarked by rapid proliferation in its host cells and requires a significant synthesis of parasite membranes. Phosphatidylethanolamine (PtdEtn) is the second major phospholipid class in T. gondii. Here, we reveal that PtdEtn is produced in the parasite mitochondrion and parasitophorous vacuole by decarboxylation of phosphatidylserine (PtdSer) and in the endoplasmic reticulum by fusion of CDP-ethanolamine and diacylglycerol. PtdEtn in the mitochondrion is synthesized by a phosphatidylserine decarboxylase (TgPSD1mt) of the type I class. TgPSD1mt harbors a targeting peptide at its N terminus that is required for the mitochondrial localization but not for the catalytic activity. Ablation of TgPSD1mt expression caused up to 45% growth impairment in the parasite mutant. The PtdEtn content of the mutant was unaffected, however, suggesting the presence of compensatory mechanisms. Indeed, metabolic labeling revealed an increased usage of ethanolamine for PtdEtn synthesis by the mutant. Likewise, depletion of nutrients exacerbated the growth defect (∼56%), which was partially restored by ethanolamine. Besides, the survival and residual growth of the TgPSD1mt mutant in the nutrient-depleted medium also indicated additional routes of PtdEtn biogenesis, such as acquisition of host-derived lipid. Collectively, the work demonstrates a metabolic cooperativity between the parasite organelles, which ensures a sustained lipid synthesis, survival and growth of T. gondii in varying nutritional milieus.

The histone code of Toxoplasma gondii comprises conserved and unique posttranslational modifications

Epigenetic gene regulation has emerged as a major mechanism for gene regulation in all eukaryotes. Histones are small, basic proteins that constitute the major protein component of chromatin, and posttranslational modifications (PTM) of histones are essential for epigenetic gene regulation. The different combinations of histone PTM form the histone code for an organism, marking functional units of chromatin that recruit macromolecular complexes that govern chromatin structure and regulate gene expression. To characterize the repertoire of Toxoplasma gondii histone PTM, we enriched histones using standard acid extraction protocols and analyzed them with several complementary middle-down and bottom-up proteomic approaches with the high-resolution Orbitrap mass spectrometer using collision-induced dissociation (CID), higher-energy collisional dissociation (HCD), and/or electron transfer dissociation (ETD) fragmentation. We identified 249 peptides with unique combinations of PTM that comprise the T. gondii histone code. T. gondii histones share a high degree of sequence conservation with human histones, and many modifications are conserved between these species. In addition, T. gondii histones have unique modifications not previously identified in other species. Finally, T. gondii histones are modified by succinylation, propionylation, and formylation, recently described histone PTM that have not previously been identified in parasitic protozoa. The characterization of the T. gondii histone code will facilitate in-depth analysis of how epigenetic regulation affects gene expression in pathogenic apicomplexan parasites and identify a new model system for elucidating the biological functions of novel histone PTM.

IMPORTANCE

Toxoplasma gondii is among the most common parasitic infections in humans. The transition between the different stages of the T. gondii life cycle are essential for parasite virulence and survival. These differentiation events are accompanied by significant changes in gene expression, and the control mechanisms for these transitions have not been elucidated. Important mechanisms that are involved in the control of gene expression are the epigenetic modifications that have been identified in several eukaryotes. T. gondii has a full complement of histone-modifying enzymes, histones, and variants. In this paper, we identify over a hundred PTM and a full repertoire of PTM combinations for T. gondii histones, providing the first large-scale characterization of the T. gondii histone code and an essential initial step for understanding how epigenetic modifications affect gene expression and other processes in this organism.

The histone code of Toxoplasma gondii comprises conserved and unique posttranslational modifications

Epigenetic gene regulation has emerged as a major mechanism for gene regulation in all eukaryotes. Histones are small, basic proteins that constitute the major protein component of chromatin, and posttranslational modifications (PTM) of histones are essential for epigenetic gene regulation. The different combinations of histone PTM form the histone code for an organism, marking functional units of chromatin that recruit macromolecular complexes that govern chromatin structure and regulate gene expression. To characterize the repertoire of Toxoplasma gondii histone PTM, we enriched histones using standard acid extraction protocols and analyzed them with several complementary middle-down and bottom-up proteomic approaches with the high-resolution Orbitrap mass spectrometer using collision-induced dissociation (CID), higher-energy collisional dissociation (HCD), and/or electron transfer dissociation (ETD) fragmentation. We identified 249 peptides with unique combinations of PTM that comprise the T. gondii histone code. T. gondii histones share a high degree of sequence conservation with human histones, and many modifications are conserved between these species. In addition, T. gondii histones have unique modifications not previously identified in other species. Finally, T. gondii histones are modified by succinylation, propionylation, and formylation, recently described histone PTM that have not previously been identified in parasitic protozoa. The characterization of the T. gondii histone code will facilitate in-depth analysis of how epigenetic regulation affects gene expression in pathogenic apicomplexan parasites and identify a new model system for elucidating the biological functions of novel histone PTM.

Toxoplasma gondii is among the most common parasitic infections in humans. The transition between the different stages of the T. gondii life cycle are essential for parasite virulence and survival. These differentiation events are accompanied by significant changes in gene expression, and the control mechanisms for these transitions have not been elucidated. Important mechanisms that are involved in the control of gene expression are the epigenetic modifications that have been identified in several eukaryotes. T. gondii has a full complement of histone-modifying enzymes, histones, and variants. In this paper, we identify over a hundred PTM and a full repertoire of PTM combinations for T. gondii histones, providing the first large-scale characterization of the T. gondii histone code and an essential initial step for understanding how epigenetic modifications affect gene expression and other processes in this organism.

Structure of the Toxoplasma gondii ROP18 kinase domain reveals a second ligand binding pocket required for acute virulence

At least a third of the human population is infected with the intracellular parasite Toxoplasma gondii, which contributes significantly to the disease burden in immunocompromised and neutropenic hosts and causes serious congenital complications when vertically transmitted to the fetus. Genetic analyses have identified the Toxoplasma ROP18 Ser/Thr protein kinase as a major factor mediating acute virulence in mice. ROP18 is secreted into the host cell during the invasion process, and its catalytic activity is required for the acute virulence phenotype. However, its precise molecular function and regulation are not fully understood. We have determined the crystal structure of the ROP18 kinase domain, which is inconsistent with a previously proposed autoinhibitory mechanism of regulation. Furthermore, a sucrose molecule bound to our structure identifies an additional ligand-binding pocket outside of the active site cleft. Mutational analysis confirms an important role for this pocket in virulence.

Apoptosis-like cell death pathways in the unicellular parasite Toxoplasma gondii following treatment with apoptosis inducers and chemotherapeutic agents: a proof-of-concept study

Ancient pathways of an apoptosis-like cell death have been identified in unicellular eukaryotes including protozoan parasites. Here, we examined programmed cell death in the apicomplexan Toxoplasma gondii which is a common intracellular pathogen of humans and warm-blooded animals. Treatment of extracellular T. gondii with various pro-apoptotic stimuli significantly induced DNA strand breaks as revealed by TUNEL and flow cytometry. Using staurosporine or miltefosine as pro-apoptotic stimuli, parasites also presented a reduced cell size, i.e. pyknosis and externalized phosphatidylserine while the plasma membrane remained intact. Importantly, staurosporine also induced DNA strand breaks in intracellular T. gondii. Data mining of the Toxoplasma genome resource identified 17 putative cell death-associated genes encoding proteases, a nuclease and several apoptosis regulators. Staurosporine-treated parasites but not controls strongly up-regulated several of these genes in a time-dependent fashion with a putative PDCD2 protein being more than 100-fold up-regulated. However, the mitochondrial membrane potential (ΔΨ_m) remained intact and caspase-like activity increased only slightly during staurosporine-triggered cell death. As compared to staurosporine, the transcriptional response of parasites to miltefosine was more restricted but PDCD2 was again strongly induced. Furthermore, T. gondii lost their ΔΨ_m and rapidly presented strong caspase-like activity during miltefosine treatment. Consequently, protease inhibitors abrogated miltefosine-induced but not staurosporine-induced Toxoplasma cell death. Finally, toxoplasmacidal drugs triggered DNA strand breaks in extracellular T. gondii. Interestingly, clindamycin also induced markers of an apoptosis-like cell death in intracellular parasites. Together, the data indicate that T. gondii possesses ancient apoptosis-like cell death machinery which can be triggered by chemotherapeutic agents.

Unravelling the Neospora caninum secretome through the secreted fraction (ESA) and quantification of the discharged tachyzoite using high-resolution mass spectrometry-based proteomics

BACKGROUND

The apicomplexan parasite Neospora caninum causes neosporosis, a disease that leads to abortion or stillbirth in cattle, generating an economic impact on the dairy and beef cattle trade. As an obligatory intracellular parasite, N. caninum needs to invade the host cell in an active manner to survive. The increase in parasite cytosolic Ca2+ upon contact with the host cell mediates critical events, including the exocytosis of phylum-specific secretory organelles and the activation of the parasite invasion motor. Because invasion is considered a requirement for pathogen survival and replication within the host, the identification of secreted proteins (secretome) involved in invasion may be useful to reveal interesting targets for therapeutic intervention.

METHODS

To chart the currently missing N. caninum secretome, we employed mass spectrometry-based proteomics to identify proteins present in the N. caninum tachyzoite using two different approaches. The first approach was identifying the proteins present in the tachyzoite-secreted fraction (ESA). The second approach was determining the relative quantification through peptide stable isotope labelling of the tachyzoites submitted to an ethanol secretion stimulus (discharged tachyzoite), expecting to identify the secreted proteins among the down-regulated group.

RESULTS

As a result, 615 proteins were identified at ESA and 2,011 proteins quantified at the discharged tachyzoite. We have analysed the connection between the secreted and the down-regulated proteins and searched for putative regulators of the secretion process among the up-regulated proteins. An interaction network was built by computational prediction involving the up- and down-regulated proteins. The mass spectrometry proteomics data have been deposited to the ProteomeXchange with identifier PXD000424.

CONCLUSIONS

The comparison between the protein abundances in ESA and their measure in the discharged tachyzoite allowed for a more precise identification of the most likely secreted proteins. Information from the network interaction and up-regulated proteins was important to recognise key proteins potentially involved in the metabolic regulation of secretion. Our results may be helpful to guide the selection of targets to be investigated against Neospora caninum and other Apicomplexan organisms.

Global analysis of apicomplexan protein S-acyl transferases reveals an enzyme essential for invasion

The advent of techniques to study palmitoylation on a whole proteome scale has revealed that it is an important reversible modification that plays a role in regulating multiple biological processes. Palmitoylation can control the affinity of a protein for lipid membranes, which allows it to impact protein trafficking, stability, folding, signalling and interactions. The publication of the palmitome of the schizont stage of Plasmodium falciparum implicated a role for palmitoylation in host cell invasion, protein export and organelle biogenesis. However, nothing is known so far about the repertoire of protein S-acyl transferases (PATs) that catalyse this modification in Apicomplexa. We undertook a comprehensive analysis of the repertoire of Asp-His-His-Cys cysteine-rich domain (DHHC-CRD) PAT family in Toxoplasma gondii and Plasmodium berghei by assessing their localization and essentiality. Unlike functional redundancies reported in other eukaryotes, some apicomplexan-specific DHHCs are essential for parasite growth, and several are targeted to organelles unique to this phylum. Of particular interest is DHHC7, which localizes to rhoptry organelles in all parasites tested, including the major human pathogen P. falciparum. TgDHHC7 interferes with the localization of the rhoptry palmitoylated protein TgARO and affects the apical positioning of the rhoptry organelles. This PAT has a major impact on T. gondii host cell invasion, but not on the parasite's ability to egress.

Gene fusion analysis in the battle against the African endemic sleeping sickness

The protozoan Trypanosoma brucei causes African Trypanosomiasis or sleeping sickness in humans, which can be lethal if untreated. Most available pharmacological treatments for the disease have severe side-effects. The purpose of this analysis was to detect novel protein-protein interactions (PPIs), vital for the parasite, which could lead to the development of drugs against this disease to block the specific interactions. In this work, the Domain Fusion Analysis (Rosetta Stone method) was used to identify novel PPIs, by comparing T. brucei to 19 organisms covering all major lineages of the tree of life. Overall, 49 possible protein-protein interactions were detected, and classified based on (a) statistical significance (BLAST e-value, domain length etc.), (b) their involvement in crucial metabolic pathways, and (c) their evolutionary history, particularly focusing on whether a protein pair is split in T. brucei and fused in the human host. We also evaluated fusion events including hypothetical proteins, and suggest a possible molecular function or involvement in a certain biological process. This work has produced valuable results which could be further studied through structural biology or other experimental approaches so as to validate the protein-protein interactions proposed here. The evolutionary analysis of the proteins involved showed that, gene fusion or gene fission events can happen in all organisms, while some protein domains are more prone to fusion and fission events and present complex evolutionary patterns.

Genetic manipulation in Δku80 strains for functional genomic analysis of Toxoplasma gondii

Targeted genetic manipulation using homologous recombination is the method of choice for functional genomic analysis to obtain a detailed view of gene function and phenotype(s). The development of mutant strains with targeted gene deletions, targeted mutations, complemented gene function, and/or tagged genes provides powerful strategies to address gene function, particularly if these genetic manipulations can be efficiently targeted to the gene locus of interest using integration mediated by double cross over homologous recombination.

Due to very high rates of nonhomologous recombination, functional genomic analysis of Toxoplasma gondii has been previously limited by the absence of efficient methods for targeting gene deletions and gene replacements to specific genetic loci. Recently, we abolished the major pathway of nonhomologous recombination in type I and type II strains of T. gondii by deleting the gene encoding the KU80 protein^1,2. The Δku80 strains behave normally during tachyzoite (acute) and bradyzoite (chronic) stages in vitro and in vivo and exhibit essentially a 100% frequency of homologous recombination. The Δku80 strains make functional genomic studies feasible on the single gene as well as on the genome scale^1-4.

Here, we report methods for using type I and type II Δku80Δhxgprt strains to advance gene targeting approaches in T. gondii. We outline efficient methods for generating gene deletions, gene replacements, and tagged genes by targeted insertion or deletion of the hypoxanthine-xanthine-guanine phosphoribosyltransferase (HXGPRT) selectable marker. The described gene targeting protocol can be used in a variety of ways in Δku80 strains to advance functional analysis of the parasite genome and to develop single strains that carry multiple targeted genetic manipulations. The application of this genetic method and subsequent phenotypic assays will reveal fundamental and unique aspects of the biology of T. gondii and related significant human pathogens that cause malaria (Plasmodium sp.) and cryptosporidiosis (Cryptosporidium).

Genetic basis for phenotypic differences between different Toxoplasma gondii type I strains

BACKGROUND

Toxoplasma gondii has a largely clonal population in North America and Europe, with types I, II and III clonal lineages accounting for the majority of strains isolated from patients. RH, a particular type I strain, is most frequently used to characterize Toxoplasma biology. However, compared to other type I strains, RH has unique characteristics such as faster growth, increased extracellular survival rate and inability to form orally infectious cysts. Thus, to identify candidate genes that could account for these parasite phenotypic differences, we determined genetic differences and differential parasite gene expression between RH and another type I strain, GT1. Moreover, as differences in host cell modulation could affect Toxoplasma replication in the host, we determined differentially modulated host processes among the type I strains through host transcriptional profiling.

RESULTS

Through whole genome sequencing, we identified 1,394 single nucleotide polymorphisms (SNPs) and insertions/deletions (indels) between RH and GT1. These SNPs/indels together with parasite gene expression differences between RH and GT1 were used to identify candidate genes that could account for type I phenotypic differences. A polymorphism in dense granule protein, GRA2, determined RH and GT1 differences in the evasion of the interferon gamma response. In addition, host transcriptional profiling identified that genes regulated by NF-ĸB, such as interleukin (IL)-12p40, were differentially modulated by the different type I strains. We subsequently showed that this difference in NF-ĸB activation was due to polymorphisms in GRA15. Furthermore, we observed that RH, but not other type I strains, recruited phosphorylated IĸBα (a component of the NF-ĸB complex) to the parasitophorous vacuole membrane and this recruitment of p- IĸBα was partially dependent on GRA2.

CONCLUSIONS

We identified candidate parasite genes that could be responsible for phenotypic variation among the type I strains through comparative genomics and transcriptomics. We also identified differentially modulated host pathways among the type I strains, and these can serve as a guideline for future studies in examining the phenotypic differences among type I strains.

Regulation of ATG8 membrane association by ATG4 in the parasitic protist Toxoplasma gondii

In the process of autophagy, the Atg8 protein is conjugated, through a ubiquitin-like system, to the lipid phosphatidylethanolamine (PE) to associate with the membrane of forming autophagosomes. There, it plays a crucial role in the genesis of these organelles and in autophagy in general. In most eukaryotes, the cysteine peptidase Atg4 processes the C terminus of cytosolic Atg8 to regulate its association with autophagosomal membranes and also delipidates Atg8 to release this protein from membranes. The parasitic protist Toxoplasma gondii contains a functional, yet apparently reduced, autophagic machinery. T. gondii Atg8 homolog, in addition to a cytosolic and occasionally autophagosomal localization, also localizes to the apicoplast, a nonphotosynthetic plastid bounded by four membranes. Our attempts to interfere with TgATG8 function showed that it appears to be essential for parasite multiplication inside its host cell. This protein also displays a peculiar C terminus that does not seem to necessitate processing prior to membrane association and yet an unusually large Toxoplasma homolog of ATG4 is predicted in the parasite genome. A TgATG4 conditional expression mutant that we have generated is severely affected in growth, and displays significant alterations at the organellar level, noticeably with a fragmentation of the mitochondrial network and a loss of the apicoplast. TgATG4-depleted parasites appear to be defective in the recycling of membrane-bound TgATG8. Overall, our data highlight a role for the TgATG8 conjugation pathway in maintaining the homeostasis of the parasite's organelles and suggest that Toxoplasma has evolved a specialized autophagic machinery with original regulation.

Structural and evolutionary adaptation of rhoptry kinases and pseudokinases, a family of coccidian virulence factors

Background

The widespread protozoan parasite Toxoplasma gondii interferes with host cell functions by exporting the contents of a unique apical organelle, the rhoptry. Among the mix of secreted proteins are an expanded, lineage-specific family of protein kinases termed rhoptry kinases (ROPKs), several of which have been shown to be key virulence factors, including the pseudokinase ROP5. The extent and details of the diversification of this protein family are poorly understood.

Results

In this study, we comprehensively catalogued the ROPK family in the genomes of Toxoplasma gondii, Neospora caninum and Eimeria tenella, as well as portions of the unfinished genome of Sarcocystis neurona, and classified the identified genes into 42 distinct subfamilies. We systematically compared the rhoptry kinase protein sequences and structures to each other and to the broader superfamily of eukaryotic protein kinases to study the patterns of diversification and neofunctionalization in the ROPK family and its subfamilies. We identified three ROPK sub-clades of particular interest: those bearing a structurally conserved N-terminal extension to the kinase domain (NTE), an E. tenella-specific expansion, and a basal cluster including ROP35 and BPK1 that we term ROPKL. Structural analysis in light of the solved structures ROP2, ROP5, ROP8 and in comparison to typical eukaryotic protein kinases revealed ROPK-specific conservation patterns in two key regions of the kinase domain, surrounding a ROPK-conserved insert in the kinase hinge region and a disulfide bridge in the kinase substrate-binding lobe. We also examined conservation patterns specific to the NTE-bearing clade. We discuss the possible functional consequences of each.

Conclusions

Our work sheds light on several important but previously unrecognized features shared among rhoptry kinases, as well as the essential differences between active and degenerate protein kinases. We identify the most distinctive ROPK-specific features conserved across both active kinases and pseudokinases, and discuss these in terms of sequence motifs, evolutionary context, structural impact and potential functional relevance.

By characterizing the proteins that enable these parasites to invade the host cell and co-opt its signaling mechanisms, we provide guidance on potential therapeutic targets for the diseases caused by coccidian parasites.

The unfolded protein response in the protozoan parasite Toxoplasma gondii features translational and transcriptional control

The unfolded protein response (UPR) is an important regulatory network that responds to perturbations in protein homeostasis in the endoplasmic reticulum (ER). In mammalian cells, the UPR features translational and transcriptional mechanisms of gene expression aimed at restoring proteostatic control. A central feature of the UPR is phosphorylation of the α subunit of eukaryotic initiation factor-2 (eIF2) by PERK (EIF2AK3/PEK), which reduces the influx of nascent proteins into the ER by lowering global protein synthesis, coincident with preferential translation of key transcription activators of genes that function to expand the processing capacity of this secretory organelle. Upon ER stress, the apicomplexan parasite Toxoplasma gondii is known to induce phosphorylation of Toxoplasma eIF2α and lower translation initiation. To characterize the nature of the ensuing UPR in this parasite, we carried out microarray analyses to measure the changes in the transcriptome and in translational control during ER stress. We determined that a collection of transcripts linked with the secretory process are induced in response to ER stress, supporting the idea that a transcriptional induction phase of the UPR occurs in Toxoplasma. Furthermore, we determined that about 500 gene transcripts showed enhanced association with translating ribosomes during ER stress. Many of these target genes are suggested to be involved in gene expression, including JmjC5, which continues to be actively translated during ER stress. This study indicates that Toxoplasma triggers a UPR during ER stress that features both translational and transcriptional regulatory mechanisms, which is likely to be important for parasite invasion and development.

Identification of Besnoitia besnoiti proteins that showed differences in abundance between tachyzoite and bradyzoite stages by difference gel electrophoresis

Bovine besnoitiosis is a chronic and debilitating disease, caused by the apicomplexan parasite Besnoitia besnoiti. Infection of cattle by B. besnoiti is governed by the tachyzoite stage, which is related to acute infection, and the bradyzoite stage gathered into macroscopic cysts located in subcutaneous tissue in the skin, mucosal membranes and sclera conjunctiva and related to persistence and chronic infection. However, the entire life cycle of this parasite and the molecular mechanisms underlying tachyzoite-to-bradyzoite conversion remain unknown. In this context, a different antigenic pattern has been observed between tachyzoite and bradyzoite extracts. Thus, to identify stage-specific proteins, a difference gel electrophoresis (DIGE) approach was used on tachyzoite and bradyzoite extracts followed by mass spectrometry (MS) analysis. A total of 130 and 132 spots were differentially expressed in bradyzoites and tachyzoites, respectively (average ratio ± 1.5, P<0.05 in t-test). Furthermore, 25 differentially expressed spots were selected and analysed by MALDI-TOF/MS. As a result, 5 up-regulated bradyzoite proteins (GAPDH, ENO1, LDH, SOD and RNA polymerase) and 5 up-regulated tachyzoite proteins (ENO2; LDH; ATP synthase; HSP70 and PDI) were identified. The present results set the basis for the identification of new proteins as drug targets. Moreover, the role of these proteins in tachyzoite-to-bradyzoite conversion and the role of the host cell environment should be a subject of further research.