Comparing the metabolic pathways of parasites and their hosts facilitates the identification of new drug targets.

Toxoplasma gondii infects approximately 30% of the world's population, causing disease primarily during pregnancy and in individuals with weakened immune systems. Toxoplasma secretes and exports effector proteins that modulate the host during infection, and several of these proteins are processed by the Golgi-associated aspartyl protease 5 (ASP5). Here, we identify ASP5 substrates by selectively enriching N-terminally derived peptides from wild-type and Δasp5 parasites. We reveal more than 2,000 unique Toxoplasma N-terminal peptides, mapping to both natural N termini and protease cleavage sites. Several of these peptides mapped directly downstream of the characterized ASP5 cleavage site, arginine-arginine-leucine (RRL). We validate candidates as true ASP5 substrates, revealing they are not processed in parasites lacking ASP5 or in wild-type parasites following mutation of the motif from RRL to ARL. All identified ASP5 substrates are dense granule proteins, and interestingly, none appear to be exported, thus differing from the analogous system in related Plasmodium spp. Instead we show that the majority of substrates reside within the parasitophorous vacuole (PV), and its membrane (the PVM), including two kinases and one phosphatase. We show that genetic deletion of WNG2 leads to attenuation in a mouse model, suggesting that this putative kinase is a new virulence factor in Toxoplasma Collectively, these data constitute the first in-depth analyses of ASP5 substrates and shed new light on the role of ASP5 as a maturase of dense granule proteins during the Toxoplasma lytic cycle.IMPORTANCEToxoplasma gondii is one of the most successful human parasites. Central to its success is the arsenal of virulence proteins introduced into the infected host cell. Several of these virulence proteins require direct maturation by the aspartyl protease ASP5, and all require ASP5 for translocation into the host cell, yet the true number of ASP5 substrates and complete repertoire of effectors is currently unknown. Here we selectively enrich N-terminally derived peptides using Terminal Amine Isotopic Labeling of Substrates (TAILS) and use quantitative proteomics to reveal novel ASP5 substrates. We identify, using two different enrichment techniques, new ASP5 substrates and their specific cleavage sites. ASP5 substrates include two kinases and one phosphatase that reside at the host-parasite interface, which are important for infection.

Toxoplasma gondii, an Apicomplexan parasite, causes significant morbidity and mortality, including severe disease in immunocompromised hosts and devastating congenital disease, with no effective treatment for the bradyzoite stage. To address this, we used the Tropical Disease Research database, crystallography, molecular modeling, and antisense to identify and characterize a range of potential therapeutic targets for toxoplasmosis. Phosphoglycerate mutase II (PGMII), nucleoside diphosphate kinase (NDK), ribulose phosphate 3-epimerase (RPE), ribose-5-phosphate isomerase (RPI), and ornithine aminotransferase (OAT) were structurally characterized. Crystallography revealed insights into the overall structure, protein oligomeric states and molecular details of active sites important for ligand recognition. Literature and molecular modeling suggested potential inhibitors and druggability. The targets were further studied with vivoPMO to interrupt enzyme synthesis, identifying the targets as potentially important to parasitic replication and, therefore, of therapeutic interest. Targeted vivoPMO resulted in statistically significant perturbation of parasite replication without concomitant host cell toxicity, consistent with a previous CRISPR/Cas9 screen showing PGM, RPE, and RPI contribute to parasite fitness. PGM, RPE, and RPI have the greatest promise for affecting replication in tachyzoites. These targets are shared between other medically important parasites and may have wider therapeutic potential.

The concept of open data has been gaining traction as a mechanism to increase data use, ensure that data are preserved over time, and accelerate discovery. While epidemiology data sets are increasingly deposited in databases and repositories, barriers to access still remain. ClinEpiDB was constructed as an open-access online resource for clinical and epidemiologic studies by leveraging the extensive web toolkit and infrastructure of the Eukaryotic Pathogen Database Resources (EuPathDB; a collection of databases covering 170+ eukaryotic pathogens, relevant related species, and select hosts) combined with a unified semantic web framework. Here we present an intuitive point-and-click website that allows users to visualize and subset data directly in the ClinEpiDB browser and immediately explore potential associations. Supporting study documentation aids contextualization, and data can be downloaded for advanced analyses. By facilitating access and interrogation of high-quality, large-scale data sets, ClinEpiDB aims to spur collaboration and discovery that improves global health.

The protozoan parasite, Toxoplasma, like many intracellular pathogens, suppresses interferon gamma (IFN-γ)-induced signal transducer and activator of transcription 1 (STAT1) activity. We exploited this well-defined host-pathogen interaction as the basis for a high-throughput screen, identifying nine transcription factors that enhance STAT1 function in the nucleus, including the orphan nuclear hormone receptor TLX. Expression profiling revealed that upon IFN-γ treatment TLX enhances the output of a subset of IFN-γ target genes, which we found is dependent on TLX binding at those loci. Moreover, infection of TLX deficient mice with the intracellular parasite Toxoplasma results in impaired production of the STAT1-dependent cytokine interleukin-12 by dendritic cells and increased parasite burden in the brain during chronic infection. These results demonstrate a previously unrecognized role for this orphan nuclear hormone receptor in regulating STAT1 signaling and host defense and reveal that STAT1 activity can be modulated in a context-specific manner by such "modifiers."

TriTrypDB (http://tritrypdb.org) is an integrated database providing access to genome-scale datasets for kinetoplastid parasites, and supporting a variety of complex queries driven by research and development needs. TriTrypDB is a collaborative project, utilizing the GUS/WDK computational infrastructure developed by the Eukaryotic Pathogen Bioinformatics Resource Center (EuPathDB.org) to integrate genome annotation and analyses from GeneDB and elsewhere with a wide variety of functional genomics datasets made available by members of the global research community, often pre-publication. Currently, TriTrypDB integrates datasets from Leishmania braziliensis, L. infantum, L. major, L. tarentolae, Trypanosoma brucei and T. cruzi. Users may examine individual genes or chromosomal spans in their genomic context, including syntenic alignments with other kinetoplastid organisms. Data within TriTrypDB can be interrogated utilizing a sophisticated search strategy system that enables a user to construct complex queries combining multiple data types. All search strategies are stored, allowing future access and integrated searches. 'User Comments' may be added to any gene page, enhancing available annotation; such comments become immediately searchable via the text search, and are forwarded to curators for incorporation into the reference annotation when appropriate.

Apicomplexan parasites release factors via specialized secretory organelles (rhoptries, micronemes) that are thought to control host cell responses. In order to explore parasite-mediated modulation of host cell signaling pathways, we exploited a phylogenomic approach to characterize the Toxoplasma gondii kinome, defining a 44 member family of coccidian-specific secreted kinases, some of which have been previously implicated in virulence. Comparative genomic analysis suggests that "ROPK" genes are under positive selection, and expression profiling demonstrates that most are differentially expressed between strains and/or during differentiation. Integrating diverse genomic-scale analyses points to ROP38 as likely to be particularly important in parasite biology. Upregulating expression of this previously uncharacterized gene in transgenic parasites dramatically suppresses transcriptional responses in the infected cell. Specifically, parasite ROP38 downregulates host genes associated with MAPK signaling and the control of apoptosis and proliferation. These results highlight the value of integrative genomic approaches in prioritizing candidates for functional validation.

Effective control of the intracellular protozoan parasite Toxoplasma gondii depends on the activation of antigen-specific CD8(+) T-cells that manage acute disease and prevent recrudescence during chronic infection. T-cell activation in turn, requires presentation of parasite antigens by MHC-I molecules on the surface of antigen presenting cells. CD8(+) T-cell epitopes have been defined for several T. gondii proteins, but it is unclear how these antigens enter into the presentation pathway. We have exploited the well-characterized model antigen ovalbumin (OVA) to investigate the ability of parasite proteins to enter the MHC-I presentation pathway, by engineering recombinant expression in various organelles. CD8(+) T-cell activation was assayed using 'B3Z' reporter cells in vitro, or adoptively-transferred OVA-specific 'OT-I' CD8(+) T-cells in vivo. As expected, OVA secreted into the parasitophorous vacuole strongly stimulated antigen-presenting cells. Lower levels of activation were observed using glycophosphatidyl inositol (GPI) anchored OVA associated with (or shed from) the parasite surface. Little CD8(+) T-cell activation was detected using parasites expressing intracellular OVA in the cytosol, mitochondrion, or inner membrane complex (IMC). These results indicate that effective presentation of parasite proteins to CD8(+) T-cells is a consequence of active protein secretion by T. gondii and escape from the parasitophorous vacuole, rather than degradation of phagocytosed parasites or parasite products.

PlasmoDB (http://PlasmoDB.org) is a functional genomic database for Plasmodium spp. that provides a resource for data analysis and visualization in a gene-by-gene or genome-wide scale. PlasmoDB belongs to a family of genomic resources that are housed under the EuPathDB (http://EuPathDB.org) Bioinformatics Resource Center (BRC) umbrella. The latest release, PlasmoDB 5.5, contains numerous new data types from several broad categories--annotated genomes, evidence of transcription, proteomics evidence, protein function evidence, population biology and evolution. Data in PlasmoDB can be queried by selecting the data of interest from a query grid or drop down menus. Various results can then be combined with each other on the query history page. Search results can be downloaded with associated functional data and registered users can store their query history for future retrieval or analysis.

Apicomplexan parasites, including the human pathogens Toxoplasma gondii and Plasmodium falciparum, employ specialized secretory organelles (micronemes, rhoptries, dense granules) to invade and survive within host cells. Because molecules secreted from these organelles function at the host/parasite interface, their identification is important for understanding invasion mechanisms, and central to the development of therapeutic strategies. Using a computational approach based on predicted functional domains, we have identified more than 600 candidate secretory organelle proteins in twelve apicomplexan parasites. Expression in transgenic T. gondii of eight proteins identified in silico confirms that all enter into the secretory pathway, and seven target to apical organelles associated with invasion. An in silico approach intended to identify possible host interacting proteins yields a dataset enriched in secretory/transmembrane proteins, including most of the antigens known to be engaged by apicomplexan parasites during infection. These domain pattern and projected interactome approaches significantly expand the repertoire of proteins that may be involved in host parasite interactions.

The closely related protozoan parasites Toxoplasma gondii and Neospora caninum display similar life cycles, subcellular ultrastructure, invasion mechanisms, metabolic pathways, and genome organization, but differ in their host range and disease pathogenesis. Type II (γ) interferon has long been known to be the major mediator of innate and adaptive immunity to Toxoplasma infection, but genome-wide expression profiling of infected host cells indicates that Neospora is a potent activator of the type I (α/β) interferon pathways typically associated with antiviral responses. Infection of macrophages from mice with targeted deletions in various innate sensing genes demonstrates that host responses to Neospora are dependent on the toll-like receptor Tlr3 and the adapter protein Trif. Consistent with this observation, RNA from Neospora elicits TLR3-dependent type I interferon responses when targeted to the host endo-lysosomal system. Although live Toxoplasma fail to induce type I interferon, heat-killed parasites do trigger this response, albeit much weaker than Neospora, and co-infection studies reveal that T. gondii actively suppresses the production of type I interferon. These findings reveal that eukaryotic pathogens can be potent inducers of type I interferon and that related parasite species interact with this pathway in distinct ways.

Although the genomes of many of the most important human and animal pathogens have now been sequenced, our understanding of the actual proteins expressed by these genomes and how well they predict protein sequence and expression is still deficient. We have used three complementary approaches (two-dimensional electrophoresis, gel-liquid chromatography linked tandem mass spectrometry and MudPIT) to analyze the proteome of Toxoplasma gondii, a parasite of medical and veterinary significance, and have developed a public repository for these data within ToxoDB, making for the first time proteomics data an integral part of this key genome resource.

The Eukaryotic Pathogen, Vector and Host Informatics Resource (VEuPathDB, https://veupathdb.org) represents the 2019 merger of VectorBase with the EuPathDB projects. As a Bioinformatics Resource Center funded by the National Institutes of Health, with additional support from the Welllcome Trust, VEuPathDB supports >500 organisms comprising invertebrate vectors, eukaryotic pathogens (protists and fungi) and relevant free-living or non-pathogenic species or hosts. Designed to empower researchers with access to Omics data and bioinformatic analyses, VEuPathDB projects integrate >1700 pre-analysed datasets (and associated metadata) with advanced search capabilities, visualizations, and analysis tools in a graphic interface. Diverse data types are analysed with standardized workflows including an in-house OrthoMCL algorithm for predicting orthology. Comparisons are easily made across datasets, data types and organisms in this unique data mining platform. A new site-wide search facilitates access for both experienced and novice users. Upgraded infrastructure and workflows support numerous updates to the web interface, tools, searches and strategies, and Galaxy workspace where users can privately analyse their own data. Forthcoming upgrades include cloud-ready application architecture, expanded support for the Galaxy workspace, tools for interrogating host-pathogen interactions, and improved interactions with affiliated databases (ClinEpiDB, MicrobiomeDB) and other scientific resources, and increased interoperability with the Bacterial & Viral BRC.

The Eukaryotic Pathogen, Vector and Host Informatics Resource (VEuPathDB, https://veupathdb.org) is a Bioinformatics Resource Center funded by the National Institutes of Health with additional funding from the Wellcome Trust. VEuPathDB supports >600 organisms that comprise invertebrate vectors, eukaryotic pathogens (protists and fungi) and relevant free-living or non-pathogenic species or hosts. Since 2004, VEuPathDB has analyzed omics data from the public domain using contemporary bioinformatic workflows, including orthology predictions via OrthoMCL, and integrated the analysis results with analysis tools, visualizations, and advanced search capabilities. The unique data mining platform coupled with >3000 pre-analyzed data sets facilitates the exploration of pertinent omics data in support of hypothesis driven research. Comparisons are easily made across data sets, data types and organisms. A Galaxy workspace offers the opportunity for the analysis of private large-scale datasets and for porting to VEuPathDB for comparisons with integrated data. The MapVEu tool provides a platform for exploration of spatially resolved data such as vector surveillance and insecticide resistance monitoring. To address the growing body of omics data and advances in laboratory techniques, VEuPathDB has added several new data types, searches and features, improved the Galaxy workspace environment, redesigned the MapVEu interface and updated the infrastructure to accommodate these changes.

The host immune response has a critical role not only in protection from human leishmaniasis but also in promoting disease severity. Although candidate gene approaches in mouse models of leishmaniasis have been extremely informative, a global understanding of the immune pathways active in lesions from human patients is lacking. To address this issue, genome-wide transcriptional profiling of Leishmania braziliensis-infected cutaneous lesions and normal skin controls was carried out. A signature of the L. braziliensis skin lesion was defined, which includes over 2,000 differentially regulated genes. Pathway-level analysis of this transcriptional response revealed key biological pathways present in cutaneous lesions, generating a testable 'metapathway' model of immunopathology and providing new insights for treatment of human leishmaniasis.

Proteomics data can supplement genome annotation efforts, for example being used to confirm gene models or correct gene annotation errors. Here, we present a large-scale proteogenomics study of two important apicomplexan pathogens: Toxoplasma gondii and Neospora caninum. We queried proteomics data against a panel of official and alternate gene models generated directly from RNASeq data, using several newly generated and some previously published MS datasets for this meta-analysis. We identified a total of 201 996 and 39 953 peptide-spectrum matches for T. gondii and N. caninum, respectively, at a 1% peptide FDR threshold. This equated to the identification of 30 494 distinct peptide sequences and 2921 proteins (matches to official gene models) for T. gondii, and 8911 peptides/1273 proteins for N. caninum following stringent protein-level thresholding. We have also identified 289 and 140 loci for T. gondii and N. caninum, respectively, which mapped to RNA-Seq-derived gene models used in our analysis and apparently absent from the official annotation (release 10 from EuPathDB) of these species. We present several examples in our study where the RNA-Seq evidence can help in correction of the current gene model and can help in discovery of potential new genes. The findings of this study have been integrated into the EuPathDB. The data have been deposited to the ProteomeXchange with identifiers PXD000297and PXD000298.

We report the cloning, expression, and characterization of the first UDP-GalNAc:polypetide N-acetylgalactosaminyltransferase (ppGalNAc-T) from the human disease-causing parasite, Toxoplasma gondii. This enzyme is also the first characterized ppGalNAc-T of protozoan origin. This type of enzyme catalyzes the initial step of mucin-type O-glycosylation, that is, the transfer of GalNAc in O-glycosidic linkage to serine and threonine residues in polypeptides. We used polymerase chain reaction amplification with degenerate primers and hybridization screening of a T. gondii cDNA library to identify this enzyme. The resulting 84-kDa type II membrane protein contains a 49-amino acid N-terminal cytoplasmic domain, a 22-amino acid hydrophobic transmembrane domain, and a 680-amino acid C-terminal lumenal domain. Conceptual translation of the cDNA sequence reveals a relatively long (i.e. 135 amino acids) stem region and the presence of several important sequence motifs. The latter include a glycosyltransferase 1 (GT1) motif containing a DXH sequence, a Gal/GalNAc-T motif, and a region homologous to ricin lectin. Northern blot analysis identified a single 5.5-kb ppGalNAc-T transcript. Comparison of the cDNA and genomic DNA sequences reveals that this transferase is encoded by 10 exons in a 10 kb region. When the recombinant construct was expressed in stably transfected Drosophila melanogaster S2 cells, the purified protein exhibited transferase activity in vitro. The identification of this enzyme in T. gondii demonstrates that this human parasite has its own enzymatic machinery for the O-glycosylation of toxoplasmal proteins.

We wished to produce a single reference gene set for honey bee (Apis mellifera). Our motivation was twofold. First, we wished to obtain an improved set of gene models with increased coverage of known genes, while maintaining gene model quality. Second, we wished to provide a single official gene list that the research community could further utilize for consistent and comparable analyses and functional annotation.

Sporozoites, the invasive form of malaria parasites transmitted by mosquitoes, are quiescent while in the insect salivary glands. Sporozoites only differentiate inside of the hepatocytes of the mammalian host. We show that sporozoite latency is an active process controlled by a eukaryotic initiation factor-2alpha (eIF2alpha) kinase (IK2) and a phosphatase. IK2 activity is dominant in salivary gland sporozoites, leading to an inhibition of translation and accumulation of stalled mRNAs into granules. When sporozoites are injected into the mammalian host, an eIF2alpha phosphatase removes the PO4 from eIF2alpha-P, and the repression of translation is alleviated to permit their transformation into liver stages. In IK2 knockout sporozoites, eIF2alpha is not phosphorylated and the parasites transform prematurely into liver stages and lose their infectivity. Thus, to complete their life cycle, Plasmodium sporozoites exploit the mechanism that regulates stress responses in eukaryotic cells.

GiardiaDB (http://GiardiaDB.org) and TrichDB (http://TrichDB.org) house the genome databases for Giardia lamblia and Trichomonas vaginalis, respectively, and represent the latest additions to the EuPathDB (http://EuPathDB.org) family of functional genomic databases. GiardiaDB and TrichDB employ the same framework as other EuPathDB sites (CryptoDB, PlasmoDB and ToxoDB), supporting fully integrated and searchable databases. Genomic-scale data available via these resources may be queried based on BLAST searches, annotation keywords and gene ID searches, GO terms, sequence motifs and other protein characteristics. Functional queries may also be formulated, based on transcript and protein expression data from a variety of platforms. Phylogenetic relationships may also be interrogated. The ability to combine the results from independent queries, and to store queries and query results for future use facilitates complex, genome-wide mining of functional genomic data.