Metabolic reconstruction identifies strain‐specific regulation of virulence in Toxoplasma gondii

Metabolic changes in Toxoplasma gondii-infected host cells measured by autofluorescence imaging

Toxoplasma gondii, the causative agent of toxoplasmosis, is an obligate intracellular parasite that infects warm-blooded vertebrates across the world. In humans, seropositivity rates of T. gondii range from 10% to 90% across communities. Despite its prevalence, few studies address how T. gondii infection changes the metabolism of host cells. In this study, we investigate how T. gondii manipulates the host cell metabolic environment by monitoring the metabolic response over time using noninvasive autofluorescence lifetime imaging of single cells, metabolite analysis, extracellular flux analysis, and reactive oxygen species (ROS) production. Autofluorescence lifetime imaging indicates that infected host cells become more oxidized and have an increased proportion of bound NAD(P)H compared to uninfected controls. Over time, infected cells also show decreases in levels of intracellular glucose and lactate, increases in oxygen consumption, and variability in ROS production. We further examined changes associated with the pre-invasion "kiss and spit" process using autofluorescence lifetime imaging, which also showed a more oxidized host cell with an increased proportion of bound NAD(P)H over 48 hours compared to uninfected controls, suggesting that metabolic changes in host cells are induced by T. gondii kiss and spit even without invasion.IMPORTANCEThis study sheds light on previously unexplored changes in host cell metabolism induced by T. gondii infection using noninvasive, label-free autofluorescence imaging. In this study, we use optical metabolic imaging (OMI) to measure the optical redox ratio (ORR) in conjunction with fluorescence lifetime imaging microscopy (FLIM) to noninvasively monitor single host cell response to T. gondii infection over 48 hours. Collectively, our results affirm the value of using autofluorescence lifetime imaging to noninvasively monitor metabolic changes in host cells over the time course of a microbial infection. Understanding this metabolic relationship between the host cell and the parasite could uncover new treatment and prevention options for T. gondii infections worldwide.

ToxoNet: A high confidence map of protein-protein interactions in Toxoplasma gondii

The apicomplexan intracellular parasite Toxoplasma gondii is a major food borne pathogen that is highly prevalent in the global population. The majority of the T. gondii proteome remains uncharacterized and the organization of proteins into complexes is unclear. To overcome this knowledge gap, we used a biochemical fractionation strategy to predict interactions by correlation profiling. To overcome the deficit of high-quality training data in non-model organisms, we complemented a supervised machine learning strategy, with an unsupervised approach, based on similarity network fusion. The resulting combined high confidence network, ToxoNet, comprises 2,063 interactions connecting 652 proteins. Clustering identifies 93 protein complexes. We identified clusters enriched in mitochondrial machinery that include previously uncharacterized proteins that likely represent novel adaptations to oxidative phosphorylation. Furthermore, complexes enriched in proteins localized to secretory organelles and the inner membrane complex, predict additional novel components representing novel targets for detailed functional characterization. We present ToxoNet as a publicly available resource with the expectation that it will help drive future hypotheses within the research community.

Architect: A tool for aiding the reconstruction of high-quality metabolic models through improved enzyme annotation

Constraint-based modeling is a powerful framework for studying cellular metabolism, with applications ranging from predicting growth rates and optimizing production of high value metabolites to identifying enzymes in pathogens that may be targeted for therapeutic interventions. Results from modeling experiments can be affected at least in part by the quality of the metabolic models used. Reconstructing a metabolic network manually can produce a high-quality metabolic model but is a time-consuming task. At the same time, current methods for automating the process typically transfer metabolic function based on sequence similarity, a process known to produce many false positives. We created Architect, a pipeline for automatic metabolic model reconstruction from protein sequences. First, it performs enzyme annotation through an ensemble approach, whereby a likelihood score is computed for an EC prediction based on predictions from existing tools; for this step, our method shows both increased precision and recall compared to individual tools. Next, Architect uses these annotations to construct a high-quality metabolic network which is then gap-filled based on likelihood scores from the ensemble approach. The resulting metabolic model is output in SBML format, suitable for constraints-based analyses. Through comparisons of enzyme annotations and curated metabolic models, we demonstrate improved performance of Architect over other state-of-the-art tools, notably with higher precision and recall on the eukaryote C. elegans and when compared to UniProt annotations in two bacterial species. Code for Architect is available at https://github.com/ParkinsonLab/Architect. For ease-of-use, Architect can be readily set up and utilized using its Docker image, maintained on Docker Hub.

Pierre J, Xiong X, Howland JG, Josselyn SA, Frankland PW, Parkinson J. Toxoplasma infection in male mice alters dopamine-sensitive behaviors and host gene expression patterns associated with neuropsychiatric disease

During chronic infection, the single celled parasite, Toxoplasma gondii, can migrate to the brain where it has been associated with altered dopamine function and the capacity to modulate host behavior, increasing risk of neurocognitive disorders. Here we explore alterations in dopamine-related behavior in a new mouse model based on stimulant (cocaine)-induced hyperactivity. In combination with cocaine, infection resulted in heightened sensorimotor deficits and impairment in prepulse inhibition response, which are commonly disrupted in neuropsychiatric conditions. To identify molecular pathways in the brain affected by chronic T. gondii infection, we investigated patterns of gene expression. As expected, infection was associated with an enrichment of genes associated with general immune response pathways, that otherwise limits statistical power to identify more informative pathways. To overcome this limitation and focus on pathways of neurological relevance, we developed a novel context enrichment approach that relies on a customized ontology. Applying this approach, we identified genes that exhibited unexpected patterns of expression arising from the combination of cocaine exposure and infection. These include sets of genes which exhibited dampened response to cocaine in infected mice, suggesting a possible mechanism for some observed behaviors and a neuroprotective effect that may be advantageous to parasite persistence. This model offers a powerful new approach to dissect the molecular pathways by which T. gondii infection contributes to neurocognitive disorders.

Cryptosporidium felis differs from other Cryptosporidium spp. in codon usage

Cryptosporidium spp. are important enteric pathogens in a wide range of vertebrates including humans. Previous comparative analysis revealed conservation in genome composition, gene content, and gene organization among Cryptosporidium spp., with a progressive reductive evolution in metabolic pathways and invasion-related proteins. In this study, we sequenced the genome of zoonotic pathogen Cryptosporidium felis and conducted a comparative genomic analysis. While most intestinal Cryptosporidium species have similar genomic characteristics and almost complete genome synteny, fewer protein-coding genes and some sequence inversions and translocations were found in the C. felis genome. The C. felis genome exhibits much higher GC content (39.6 %) than other Cryptosporidium species (24.3–32.9 %), especially at the third codon position (GC3) of protein-coding genes. Thus, C. felis has a different codon usage, which increases the use of less energy costly amino acids (Gly and Ala) encoded by GC-rich codons. While the tRNA usage is conserved among Cryptosporidium species, consistent with its higher GC content, C. felis uses a unique tRNA for GTG for valine instead of GTA in other Cryptosporidium species. Both mutational pressures and natural selection are associated with the evolution of the codon usage in Cryptosporidium spp., while natural selection seems to drive the codon usage in C. felis. Other unique features of the C. felis genome include the loss of the entire traditional and alternative electron transport systems and several invasion-related proteins. Thus, the preference for the use of some less energy costly amino acids in C. felis may lead to a more harmonious parasite–host interaction, and the strengthened host-adaptation is reflected by the further reductive evolution of metabolism and host invasion-related proteins.

Recent achievements and doors opened for coccidian parasite research and development through transcriptomics of enteric sexual stages

The Coccidia is the largest group of parasites within the Apicomplexa, a phylum of unicellular, obligate parasites characterized by the possession of an apical complex of organelles and structures in the asexual stages of their life cycles, as well as by a sexual reproductive phase that occurs enterically in host animals. Coccidian sexual reproduction involves morphologically distinct microgametes and macrogametes that combine to form a diploid zygote and, ultimately, following meiosis and mitosis, haploid, infectious sporozoites, inside sporocysts within an oocyst. Recent transcriptomic analyses have identified genes involved in coccidian sexual stage development and reproduction, including genes encoding for microgamete- and macrogamete-specific proteins with roles in gamete motility, fusion and fertilization, and in the formation of the resilient oocyst wall that allows coccidians to persist for long periods in the environment. Transcriptomics has also provided important clues about the regulation of gene expression in the transformation of parasites from one developmental stage to the next, a complex sequence of events that may involve transcription factors such as the apicomplexan Apetala2 (ApiAP2) family, alternative splicing, regulatory RNAs and MORC (a microrchida homologue and regulator of sexual stage development in Toxoplasma gondii). The molecular dissection of coccidian sexual development and reproduction by transcriptomic analyses may lead to the development of novel transmission-blocking strategies.

Amino Acid Metabolism in Apicomplexan Parasites

Obligate intracellular pathogens have coevolved with their host, leading to clever strategies to access nutrients, to combat the host's immune response, and to establish a safe niche for intracellular replication. The host, on the other hand, has also developed ways to restrict the replication of invaders by limiting access to nutrients required for pathogen survival. In this review, we describe the recent advancements in both computational methods and high-throughput -omics techniques that have been used to study and interrogate metabolic functions in the context of intracellular parasitism. Specifically, we cover the current knowledge on the presence of amino acid biosynthesis and uptake within the Apicomplexa phylum, focusing on human-infecting pathogens: Toxoplasma gondii and Plasmodium falciparum. Given the complex multi-host lifecycle of these pathogens, we hypothesize that amino acids are made, rather than acquired, depending on the host niche. We summarize the stage specificities of enzymes revealed through transcriptomics data, the relevance of amino acids for parasite pathogenesis in vivo, and the role of their transporters. Targeting one or more of these pathways may lead to a deeper understanding of the specific contributions of biosynthesis versus acquisition of amino acids and to design better intervention strategies against the apicomplexan parasites.

Proteomic Characterization of Host-Pathogen Interactions during Bovine Trophoblast Cell Line Infection by Neospora caninum

Despite the importance of bovine neosporosis, relevant knowledge gaps remain concerning the pathogenic mechanisms of Neospora caninum. Infection of the placenta is a crucial event in the pathogenesis of the disease; however, very little is known about the relation of the parasite with this target organ. Recent studies have shown that isolates with important variations in virulence also show different interactions with the bovine trophoblast cell line F3 in terms of proliferative capacity and transcriptome host cell modulation. Herein, we used the same model of infection to study the interaction of Neospora with these target cells at the proteomic level using LC-MS/MS over the course of the parasite lytic cycle. We also analysed the proteome differences between high- (Nc-Spain7) and low-virulence (Nc-Spain1H) isolates. The results showed that mitochondrial processes and metabolism were the main points of Neospora-host interactions. Interestingly, Nc-Spain1H infection showed a higher level of influence on the host cell proteome than Nc-Spain7 infection.

Modeling the metabolic interplay between a parasitic worm and its bacterial endosymbiont allows the identification of novel drug targets

The filarial nematode Brugia malayi represents a leading cause of disability in the developing world, causing lymphatic filariasis in nearly 40 million people. Currently available drugs are not well-suited to mass drug administration efforts, so new treatments are urgently required. One potential vulnerability is the endosymbiotic bacteria Wolbachia-present in many filariae-which is vital to the worm. Genome scale metabolic networks have been used to study prokaryotes and protists and have proven valuable in identifying therapeutic targets, but have only been applied to multicellular eukaryotic organisms more recently. Here, we present iDC625, the first compartmentalized metabolic model of a parasitic worm. We used this model to show how metabolic pathway usage allows the worm to adapt to different environments, and predict a set of 102 reactions essential to the survival of B. malayi. We validated three of those reactions with drug tests and demonstrated novel antifilarial properties for all three compounds.

Comparative genomic analysis of three intestinal species reveals reductions in secreted pathogenesis determinants in bovine-specific and non-pathogenic Cryptosporidium species

The three common intestinal Cryptosporidium species in cattle differ significantly in host range, pathogenicity and public health significance. While Cryptosporidium parvum is pathogenic in pre-weaned calves and has a broad host range, C. bovis and C. ryanae are largely non-pathogenic and bovine-specific species in post-weaned calves. Thus far, only the genome of C. parvum has been sequenced. To improve our understanding of the genetic determinants of biological differences among Cryptosporidium spcies, we sequenced the genomes of C. bovis and C. ryanae and conducted a comparative genomics analysis. The genome of C. bovis has a gene content and organization more similar to C. ryanae than to other Cryptosporidium species sequenced to date; the level of similarity in amino acid and nucleotide sequences between the two species is 75.2 and 69.4 %, respectively. A total of 3723 and 3711 putative protein-encoding genes were identified in the genomes of C. bovis and C. ryanae, respectively, which are fewer than the 3981 in C. parvum. Metabolism is similar among the three species, although energy production pathways are further reduced in C. bovis and C. ryanae. Compared with C. parvum, C. bovis and C. ryanae have lost 14 genes encoding mucin-type glycoproteins and three for insulinase-like proteases. Other gene gains and losses in the two bovine-specific and non-pathogenic species also involve the secretory pathogenesis determinants (SPDs); they have lost all genes encoding MEDLE, FLGN and SKSR proteins, and two of the three genes for NFDQ proteins, but have more genes encoding secreted WYLE proteins, secreted leucine-rich proteins and GPI-anchored adhesin PGA18. The only major difference between C. bovis and C. ryanae is in nucleotide metabolism. In addition, half of the highly divergent genes between C. bovis and C. ryanae encode secreted or membrane-bound proteins. Therefore, C. bovis and C. ryanae have gene organization and metabolic pathways similar to C. parvum, but have lost some invasion-associated mucin glycoproteins, insulinase-like proteases, MEDLE secretory proteins and other SPDs. The multiple gene families under positive selection, such as helicase-associated domains, AMP-binding domains, protein kinases, mucins, insulinases and TRAPs could contribute to differences in host specificity and pathogenicity between C. parvum and C. bovis. Biological studies should be conducted to assess the contribution of these copy number variations to the narrow host range and reduced pathogenicity of C. bovis and C. ryanae.

Comparative tachyzoite proteome analyses among six Neospora caninum isolates with different virulence

The biological variability among Neospora caninum isolates has been widely shown, however, the molecular basis that determines this diversity has not been thoroughly elucidated to date. The latest studies have focused on a limited number of isolates. Therefore, the goal of the present study was to compare the proteome of a larger number of N. caninum isolates with different origins and virulence. Label-free LC-MS/MS was used to investigate the tachyzoite proteomic differences among Nc-Bahia, Nc-Spain4H and Nc-Spain7, representing high virulence isolates and Nc-Ger6, Nc-Spain2H and Nc-Spain1H, representing low virulence isolates. Pairwise comparisons between all isolates and between high virulence and low virulence groups identified a subset of proteins with higher abundance in high virulence isolates. These proteins were involved in energy and redox metabolism, and DNA/RNA processing, which might determine the faster growth rates and parasite survival of the high virulence isolates. Highlighted proteins included a predicted member of the rhoptry kinase family ROP20 specific for N. caninum, Bradyzoite pseudokinase 1 and several dense granule proteins. DNA polymerase, which was more abundant in all high virulence isolates in all comparisons, might also be implicated in virulence. These results reveal insights into possible mechanisms involved in specific phenotypic traits and virulence in N. caninum, and the relevance of these candidate proteins for N. caninum virulence deserves further investigation.

Functional and Computational Genomics Reveal Unprecedented Flexibility in Stage-Specific Toxoplasma Metabolism

To survive and proliferate in diverse host environments with varying nutrient availability, the obligate intracellular parasite Toxoplasma gondii reprograms its metabolism. We have generated and curated a genome-scale metabolic model (iTgo) for the fast-replicating tachyzoite stage, harmonized with experimentally observed phenotypes. To validate the importance of four metabolic pathways predicted by the model, we have performed in-depth in vitro and in vivo phenotyping of mutant parasites including targeted metabolomics and CRISPR-Cas9 fitness screening of all known metabolic genes. This led to unexpected insights into the remarkable flexibility of the parasite, addressing the dependency on biosynthesis or salvage of fatty acids (FAs), purine nucleotides (AMP and GMP), a vitamin (pyridoxal-5P), and a cofactor (heme) in both the acute and latent stages of infection. Taken together, our experimentally validated metabolic network leads to a deeper understanding of the parasite's biology, opening avenues for the development of therapeutic intervention against apicomplexans.

Vitamin and cofactor acquisition in apicomplexans: Synthesis versus salvage

The Apicomplexa phylum comprises diverse parasitic organisms that have evolved from a free-living ancestor. These obligate intracellular parasites exhibit versatile metabolic capabilities reflecting their capacity to survive and grow in different hosts and varying niches. Determined by nutrient availability, they either use their biosynthesis machineries or largely depend on their host for metabolite acquisition. Because vitamins cannot be synthesized by the mammalian host, the enzymes required for their synthesis in apicomplexan parasites represent a large repertoire of potential therapeutic targets. Here, we review recent advances in metabolic reconstruction and functional studies coupled to metabolomics that unravel the interplay between biosynthesis and salvage of vitamins and cofactors in apicomplexans. A particular emphasis is placed on Toxoplasma gondii, during both its acute and latent stages of infection.

Metabolic interactions between Toxoplasma gondii and its host

Toxoplasma gondii is an obligate intracellular parasite belonging to the phylum Apicomplexa that infects all warm-blooded animals, including humans. T. gondii can replicate in every nucleated host cell by orchestrating metabolic interactions to derive crucial nutrients. In this review, we summarize the current status of known metabolic interactions of T. gondii with its host cell and discuss open questions and promising experimental approaches that will allow further dissection of the host-parasite interface and discovery of ways to efficiently target both tachyzoite and bradyzoite forms of T. gondii, which are associated with acute and chronic infection, respectively.

A Genome-Scale Metabolic Reconstruction of Phytophthora infestans With the Integration of Transcriptional Data Reveals the Key Metabolic Patterns Involved in the Interaction of Its Host

Phytophthora infestans, the causal agent of late blight disease, affects potatoes and tomatoes worldwide. This plant pathogen has a hemibiotrophic lifestyle, having an initial biotrophic infection phase during which the pathogen spreads within the host tissue, followed by a necrotrophic phase in which host cell death is induced. Although increasing information is available on the molecular mechanisms, underlying the distinct phases of the hemibiotrophic lifestyle, studies that consider the entire metabolic processes in the pathogen while undergoing the biotrophic, transition to necrotrophic, and necrotrophic phases have not been conducted. In this study, the genome-scale metabolic reconstruction of P. infestans was achieved. Subsequently, transcriptional data (microarrays, RNA-seq) was integrated into the metabolic reconstruction to obtain context-specific (metabolic) models (CSMs) of the infection process, using constraint-based reconstruction and analysis. The goal was to identify specific metabolic markers for distinct stages of the pathogen's life cycle. Results indicate that the overall metabolism show significant changes during infection. The most significant changes in metabolism were observed at the latest time points of infection. Metabolic activity associated with purine, pyrimidine, fatty acid, fructose and mannose, arginine, glycine, serine, and threonine amino acids appeared to be the most important metabolisms of the pathogen during the course of the infection, showing high number of reactions associated with them and expression switches at important stages of the life cycle. This study provides a framework for future throughput studies of the metabolic changes during the hemibiotrophic life cycle of this important plant pathogen.

Biomedical applications of genome-scale metabolic network reconstructions of human pathogens

The growing global threat of antibiotic resistant human pathogens has coincided with improved methods for developing and using genome-scale metabolic network reconstructions. Consequently, there has been an increase in the number of high-quality reconstructions of relevant human and zoonotic pathogens. Novel biomedical applications of pathogen reconstructions focus on three key aspects of pathogen behavior: the evolution of antibiotic resistance, virulence factor production, and host-pathogen interactions. New methods using these reconstructions aim to improve understanding of microbe pathogenicity and guide the development of new therapeutic strategies. This review summarizes the latest ways that genome-scale metabolic network reconstructions have been used to study human pathogens and suggests future applications with the potential to mitigate infectious disease.

Modulation of cis- and trans- Golgi and the Rab9A-GTPase during infection by Besnoitia besnoiti, Toxoplasma gondii and Neospora caninum

Like most intracellular pathogens, the apicomplexan parasites Besnoitia besnoiti, Toxoplasma gondii and Neospora caninum scavenge metabolites from their host cells. Recruitment of the Golgi complex to the vicinity of the parasitophorous vacuole (PV) is likely to aid in this process. In this work, we comparatively assessed B. besnoiti, T. gondii and N. caninum infected human retinal pigmented epithelial (hTERT-RPE-1) cells at 24 h post-infection and used antibodies to confirm Golgi ribbon compaction in B. besnoiti, and Golgi ribbon dispersion in T. gondii, while no alteration in Golgi morphology was seen in N. caninum infected cells. In either case, the Golgi stacks of infected cells contained both cis- (GM130) and trans- (TGN46) Golgi proteins. The localization of Rab9A, an important regulator of endosomal trafficking, was also studied. GFP-tagged Rab9A was recruited to the vicinity of the PV of all three parasites. Toxoplasma-infected cells exhibited increased expression of Rab9A in comparison to non-infected cells. However, Rab9A expression levels remained unaltered upon infection with N. caninum and B. besnoiti tachyzoites. In contrast to Rab9A, a GFP-tagged dominant negative mutant form of Rab9A (Rab9A DN), was not recruited to the PV, and the expression of Rab9A DN did not affect host cell invasion nor replication by all three parasites. Thus, B. besnoiti, T. gondii and N. caninum show similarities but also differences in how they affect constituents of the endosomal/secretory pathways.

Integrative transcriptome and proteome analyses define marked differences between Neospora caninum isolates throughout the tachyzoite lytic cycle

Neospora caninum is one of the main causes of transmissible abortion in cattle. Intraspecific variations in virulence have been widely shown among N. caninum isolates. However, the molecular basis governing such variability have not been elucidated to date. In this study label free LC-MS/MS was used to investigate proteome differences between the high virulence isolate Nc-Spain7 and the low virulence isolate Nc-Spain1H throughout the tachyzoite lytic cycle. The results showed greater differences in the abundance of proteins at invasion and egress with 77 and 62 proteins, respectively. During parasite replication, only 19 proteins were differentially abundant between isolates. The microneme protein repertoire involved in parasite invasion and egress was more abundant in the Nc-Spain1H isolate, which displays a lower invasion rate. Rhoptry and dense granule proteins, proteins related to metabolism and stress responses also showed differential abundances between isolates. Comparative RNA-Seq analyses during tachyzoite egress were also performed, revealing an expression profile of genes associated with the bradyzoite stage in the low virulence Nc-Spain1H isolate. The differences in proteome and RNA expression profiles between these two isolates reveal interesting insights into likely mechanisms involved in specific phenotypic traits and virulence in N. caninum.

SIGNIFICANCE

The molecular basis that governs biological variability in N. caninum and the pathogenesis of neosporosis has not been well-established yet. This is the first study in which high throughput technology of LC-MS/MS and RNA-Seq is used to investigate differences in the proteome and transcriptome between two well-characterized isolates. Both isolates displayed different proteomes throughout the lytic cycle and the transcriptomes also showed marked variations but were inconsistent with the proteome results. However, both datasets identified a pre-bradyzoite status of the low virulence isolate Nc-Spain1H. This study reveals interesting insights into likely mechanisms involved in virulence in N. caninum and shed light on a subset of proteins that are potentially involved in the pathogenesis of this parasite.

Genomics of apicomplexan parasites

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

Toxoplasma growth in vitro is dependent on exogenous tyrosine and is independent of AAH2 even in tyrosine-limiting conditions

Toxoplasma gondii is an obligate intracellular parasite capable of infecting virtually all nucleated cell types in almost all warm-blooded animals. Interestingly, Toxoplasma has a relatively full repertoire of amino acid biosynthetic machinery, perhaps reflecting its broad host range and, consequently, its need to adapt to a wide array of amino acid resources. Although Toxoplasma has been shown to be auxotrophic for tryptophan and arginine, it has not previously been determined if Toxoplasma is also auxotrophic for tyrosine. Toxoplasma tachyzoites and bradyzoites were recently found to express an amino acid hydroxylase (AAH2) that is capable of synthesizing tyrosine and dihydroxyphenylalanine (DOPA) from phenylalanine; however, the role of AAH2 in tachyzoite and bradyzoite infection has not yet been identified. To determine if Toxoplasma requires exogenous tyrosine for growth, we performed growth assays on tachyzoites and bradyzoites in nutrient-rich media titrated with varying amounts of tyrosine. We found that Toxoplasma tachyzoites form significantly smaller plaques in tyrosine-limiting media in a dose-dependent manner and that this phenotype is not affected by deletion of TgAAH2. To determine if bradyzoites require exogenous tyrosine for growth, we induced differentiation from tachyzoites in vitro in tyrosine-limiting media and found that replication and vacuole number are all decreased in tyrosine-deficient media. Importantly, culture of confluent human fibroblasts in tyrosine-deficient media does not affect their viability, indicating that, at least in vitro, the need for tyrosine is at the level of Toxoplasma, not the host cell supporting its growth.

Lessons from the genomes and transcriptomes of filarial nematodes

Human filarial infections are a leading cause of morbidity in the developing world. While a small arsenal of drugs exists to treat these infections, there remains a tremendous need for the development of additional interventions. Recent genome sequences and transcriptome analyses of filarial nematodes have provided novel biological insight and allowed for the prediction of novel drug targets as well as potential vaccine candidates. In this review, we discuss the currently available data, insights gained into the metabolism of these organisms, and how the filaria field can move forward by leveraging these data.

The genome of Onchocerca volvulus, agent of river blindness

Human onchocerciasis is a serious neglected tropical disease caused by the filarial nematode Onchocerca volvulus that can lead to blindness and chronic disability. Control of the disease relies largely on mass administration of a single drug, and the development of new drugs and vaccines depends on a better knowledge of parasite biology. Here, we describe the chromosomes of O. volvulus and its Wolbachia endosymbiont. We provide the highest-quality sequence assembly for any parasitic nematode to date, giving a glimpse into the evolution of filarial parasite chromosomes and proteomes. This resource was used to investigate gene families with key functions that could be potentially exploited as targets for future drugs. Using metabolic reconstruction of the nematode and its endosymbiont, we identified enzymes that are likely to be essential for O. volvulus viability. In addition, we have generated a list of proteins that could be targeted by Federal-Drug-Agency-approved but repurposed drugs, providing starting points for anti-onchocerciasis drug development.

Latent toxoplasmosis is associated with neurocognitive impairment in young adults with and without chronic HIV infection

We evaluated the impact of latent toxoplasmosis (LT) on neurocognitive (NC) and neurobehavioural functioning in young adults with and without chronic HIV infection, using a standardised NC test battery, self-reported Beck Depression Inventory, Frontal System Behavior Scale, MINI-International Neuropsychiatric Interview and risk-assessment battery. 194 young adults (median age 24years, 48.2% males) with chronic HIV infection (HIV+) since childhood and 51 HIV seronegative (HIV-) participants were included. HIV+ individuals had good current immunological status (median CD4: 479 cells/μl) despite a low CD4 nadir (median: 93 cells/μl). LT (positive anti-Toxoplasma IgG antibodies) was present in one third of participants. The impairment rates in the HIV- with and without Toxo were not significantly different (p=0.17). However, we observed an increasing trend (p<0.001) in impairment rates with HIV and LT status: HIV-/LT- (6.1%); HIV-/LT+ (22%), HIV+/LT- (31%), HIV+/LT+ (49%). In a multivariable analysis using the entire study group there were main effects on cognition for HIV and also for LT. Within the HIV+ group LT was associated with worse performance globally (p=0.006), in memory (p=0.009), speed of information processing (p=0.01), verbal (p=0.02) and learning (p=0.02) domains. LT was not associated with depressive symptoms, frontal systems dysfunction or risk behaviors in any of the groups. HIV participants with lower Toxoplasma antibody concentration had worse NC performance, with higher GDS values (p=0.03) and worse learning (p=0.002), memory (p=0.006), speed of information processing (p=0.01) T scores. Latent Toxoplasmosis may contribute to NC impairment in young adults, including those with and without chronic HIV infection.

No association between latent toxoplasmosis and multiple body measures in U.S. adults

Toxoplasma gondii (Nicolle et Manceaux, 1908) is an intracellular parasite that can cause ongoing latent infection persisting for the duration of a non-definitive host's life. Affecting approximately one-third of the world's population, latent toxoplasmosis has been associated with neuropsychological outcomes and a previous report suggested an association between latent toxoplasmosis and adult height. Given the large number of people with latent toxoplasmosis and its potential associations with human height, we sought to better understand the association between latent toxoplasmosis and human morphology by evaluating seropositivity for T. gondii and multiple body measures reported in the National Health and Nutrition Examination Survey III (NHANES III) and in the more recent continuous NHANES data sets from the United States Centers for Disease Control and Prevention for which data on T. gondii are available. In these analyses, latent toxoplasmosis was not associated with any of the body measures assessed in the NHANES datasets even after taking into account interactions between latent toxoplasmosis and testosterone suggesting that in these samples, latent toxoplasmosis is not associated with adult morphology including height.

Metabolic flux analyses of Pseudomonas aeruginosa cystic fibrosis isolates

Pseudomonas aeruginosa is a metabolically versatile wide-ranging opportunistic pathogen. In humans P. aeruginosa causes infections of the skin, urinary tract, blood, and the lungs of Cystic Fibrosis patients. In addition, P. aeruginosa's broad environmental distribution, relatedness to biotechnologically useful species, and ability to form biofilms have made it the focus of considerable interest. We used ¹³C metabolic flux analysis (MFA) and flux balance analysis to understand energy and redox production and consumption and to explore the metabolic phenotypes of one reference strain and five strains isolated from the lungs of cystic fibrosis patients. Our results highlight the importance of the oxidative pentose phosphate and Entner-Doudoroff pathways in P. aeruginosa growth. Among clinical strains we report two divergent metabolic strategies and identify changes between genetically related strains that have emerged during a chronic infection of the same patient. MFA revealed that the magnitude of fluxes through the glyoxylate cycle correlates with growth rates.

Determination of Ancylostoma caninum ova viability using metabolic profiling

Differentiation between viable and non-viable hookworm ova in environmental samples is necessary in order to implement strategies to mitigate re-infections in endemic regions. In this study, an untargeted metabolic profiling method was developed that utilised gas chromatography-mass spectrometry (GC-MS) in order to investigate hookworm ova viability. Ancylostoma caninum was used to investigate the metabolites within viable and non-viable ova. Univariate and multivariate statistical analyses of the data resulted in the identification of 53 significant metabolites across all hookworm ova samples. The major compounds observed in viable and non-viable hookworm ova were tetradecanoic acid, commonly known as myristic acid [fold change (FC) = 0.4], and dodecanoic acid, commonly known as lauric acid (FC = 0.388). Additionally, the viable ova had self-protecting metabolites such as prostaglandins, a typical feature absent in non-viable ova. The results of this study demonstrate that metabolic profiling using GC-MS methods can be used to determine the viability of canine hookworm ova. Further studies are needed to assess the applicability of metabolic profiling using GC-MS to detect viable hookworm ova in the mixed (viable and non-viable) populations from environmental samples and identify the metabolites specific to human hookworm species.

Local admixture of amplified and diversified secreted pathogenesis determinants shapes mosaic Toxoplasma gondii genomes

Toxoplasma gondii is among the most prevalent parasites worldwide, infecting many wild and domestic animals and causing zoonotic infections in humans. T. gondii differs substantially in its broad distribution from closely related parasites that typically have narrow, specialized host ranges. To elucidate the genetic basis for these differences, we compared the genomes of 62 globally distributed T. gondii isolates to several closely related coccidian parasites. Our findings reveal that tandem amplification and diversification of secretory pathogenesis determinants is the primary feature that distinguishes the closely related genomes of these biologically diverse parasites. We further show that the unusual population structure of T. gondii is characterized by clade-specific inheritance of large conserved haploblocks that are significantly enriched in tandemly clustered secretory pathogenesis determinants. The shared inheritance of these conserved haploblocks, which show a different ancestry than the genome as a whole, may thus influence transmission, host range and pathogenicity.

Cancer cell metabolism as new targets for novel designed therapies

Metabolic processes are altered in cancer cells, which obtain advantages from this metabolic reprogramming in terms of energy production and synthesis of biomolecules that sustain their uncontrolled proliferation. Due to the conceptual progresses in the last decade, metabolic reprogramming was recently included as one of the new hallmarks of cancer. The advent of high-throughput technologies to amass an abundance of omic data, together with the development of new computational methods that allow the integration and analysis of omic data by using genome-scale reconstructions of human metabolism, have increased and accelerated the discovery and development of anticancer drugs and tumor-specific metabolic biomarkers. Here we review and discuss the latest advances in the context of metabolic reprogramming and the future in cancer research.

Pan-phylum Comparison of Nematode Metabolic Potential

Nematodes are among the most important causative pathogens of neglected tropical diseases. The increased availability of genomic and transcriptomic data for many understudied nematode species provides a great opportunity to investigate different aspects of their biology. Increasingly, metabolic potential of pathogens is recognized as a critical determinant governing their development, growth and pathogenicity. Comparing metabolic potential among species with distinct trophic ecologies can provide insights on overall biology or molecular adaptations. Furthermore, ascertaining gene expression at pathway level can help in understanding metabolic dynamics over development. Comparison of biochemical pathways (or subpathways, i.e. pathway modules) among related species can also retrospectively indicate potential mistakes in gene-calling and functional annotation. We show with numerous illustrative case studies that comparisons at the level of pathway modules have the potential to uncover biological insights while remaining computationally tractable. Here, we reconstruct and compare metabolic modules found in the deduced proteomes of 13 nematodes and 10 non-nematode species (including hosts of the parasitic nematode species). We observed that the metabolic potential is, in general, concomitant with phylogenetic and/or ecological similarity. Varied metabolic strategies are required among the nematodes, with only 8 out of 51 pathway modules being completely conserved. Enzyme comparison based on topology of metabolic modules uncovered diversification between parasite and host that can potentially guide therapeutic intervention. Gene expression data from 4 nematode species were used to study metabolic dynamics over their life cycles. We report unexpected differential metabolism between immature and mature microfilariae of the human filarial parasite Brugia malayi. A set of genes potentially important for parasitism is also reported, based on an analysis of gene expression in C. elegans and the human hookworm Necator americanus. We illustrate how analyzing and comparing metabolism at the level of pathway modules can improve existing knowledge of nematode metabolic potential and can provide parasitism related insights. Our reconstruction and comparison of nematode metabolic pathways at a pan-phylum and inter-phylum level enabled determination of phylogenetic restrictions and differential expression of pathways. A visualization of our results is available at http://nematode.net and the program for identification of module completeness (modDFS) is freely available at SourceForge. The methods reported will help biologists to predict biochemical potential of any organism with available deduced proteome, to direct experiments and test hypotheses.

We reconstructed metabolic pathways of 23 organisms including 13 nematode species, using their complete deduced protein coding sequences and compared them to 10 non-nematodes. We observed that metabolic potential availability is concomitant with phylogenetic and/or ecological similarity, with the exceptions providing interesting case studies. We also studied changes in metabolic profiles under different developmental stages of 4 nematode species using stage-specific transcriptomic data. A comparison of the variation patterns in these profiles led to recognition of modules that share metabolic profiles at various life-cycle stages or during development. The undertaken analysis improved genome annotation and the obtained results provided insight into parasitism, resulting in identification of taxonomically-restricted pathways and enzymes that may provide new mechanisms for control of nematode infections.

Metabolic Needs and Capabilities of Toxoplasma gondii through Combined Computational and Experimental Analysis

Toxoplasma gondii is a human pathogen prevalent worldwide that poses a challenging and unmet need for novel treatment of toxoplasmosis. Using a semi-automated reconstruction algorithm, we reconstructed a genome-scale metabolic model, ToxoNet1. The reconstruction process and flux-balance analysis of the model offer a systematic overview of the metabolic capabilities of this parasite. Using ToxoNet1 we have identified significant gaps in the current knowledge of Toxoplasma metabolic pathways and have clarified its minimal nutritional requirements for replication. By probing the model via metabolic tasks, we have further defined sets of alternative precursors necessary for parasite growth. Within a human host cell environment, ToxoNet1 predicts a minimal set of 53 enzyme-coding genes and 76 reactions to be essential for parasite replication. Double-gene-essentiality analysis identified 20 pairs of genes for which simultaneous deletion is deleterious. To validate several predictions of ToxoNet1 we have performed experimental analyses of cytosolic acetyl-CoA biosynthesis. ATP-citrate lyase and acetyl-CoA synthase were localised and their corresponding genes disrupted, establishing that each of these enzymes is dispensable for the growth of T. gondii, however together they make a synthetic lethal pair.

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.

Systems-based analysis of the Sarcocystis neurona genome identifies pathways that contribute to a heteroxenous life cycle

Sarcocystis neurona is a member of the coccidia, a clade of single-celled parasites of medical and veterinary importance including Eimeria, Sarcocystis, Neospora, and Toxoplasma. Unlike Eimeria, a single-host enteric pathogen, Sarcocystis, Neospora, and Toxoplasma are two-host parasites that infect and produce infectious tissue cysts in a wide range of intermediate hosts. As a genus, Sarcocystis is one of the most successful protozoan parasites; all vertebrates, including birds, reptiles, fish, and mammals are hosts to at least one Sarcocystis species. Here we sequenced Sarcocystis neurona, the causal agent of fatal equine protozoal myeloencephalitis. The S. neurona genome is 127 Mbp, more than twice the size of other sequenced coccidian genomes. Comparative analyses identified conservation of the invasion machinery among the coccidia. However, many dense-granule and rhoptry kinase genes, responsible for altering host effector pathways in Toxoplasma and Neospora, are absent from S. neurona. Further, S. neurona has a divergent repertoire of SRS proteins, previously implicated in tissue cyst formation in Toxoplasma. Systems-based analyses identified a series of metabolic innovations, including the ability to exploit alternative sources of energy. Finally, we present an S. neurona model detailing conserved molecular innovations that promote the transition from a purely enteric lifestyle (Eimeria) to a heteroxenous parasite capable of infecting a wide range of intermediate hosts.

Sarcocystis neurona is a member of the coccidia, a clade of single-celled apicomplexan parasites responsible for major economic and health care burdens worldwide. A cousin of Plasmodium, Cryptosporidium, Theileria, and Eimeria, Sarcocystis is one of the most successful parasite genera; it is capable of infecting all vertebrates (fish, reptiles, birds, and mammals—including humans). The past decade has witnessed an increasing number of human outbreaks of clinical significance associated with acute sarcocystosis. Among Sarcocystis species, S. neurona has a wide host range and causes fatal encephalitis in horses, marine mammals, and several other mammals. To provide insights into the transition from a purely enteric parasite (e.g., Eimeria) to one that forms tissue cysts (Toxoplasma), we present the first genome sequence of S. neurona. Comparisons with other coccidian genomes highlight the molecular innovations that drive its distinct life cycle strategies.

Metabolic sensor governing bacterial virulence in Staphylococcus aureus

An effective metabolism is essential to all living organisms, including the important human pathogen Staphylococcus aureus. To establish successful infection, S. aureus must scavenge nutrients and coordinate its metabolism for proliferation. Meanwhile, it also must produce an array of virulence factors to interfere with host defenses. However, the ways in which S. aureus ties its metabolic state to its virulence regulation remain largely unknown. Here we show that citrate, the first intermediate of the tricarboxylic acid (TCA) cycle, binds to and activates the catabolite control protein E (CcpE) of S. aureus. Using structural and site-directed mutagenesis studies, we demonstrate that two arginine residues (Arg145 and Arg256) within the putative inducer-binding cavity of CcpE are important for its allosteric activation by citrate. Microarray analysis reveals that CcpE tunes the expression of 126 genes that comprise about 4.7% of the S. aureus genome. Intriguingly, although CcpE is a major positive regulator of the TCA-cycle activity, its regulon consists predominantly of genes involved in the pathogenesis of S. aureus. Moreover, inactivation of CcpE results in increased staphyloxanthin production, improved ability to acquire iron, increased resistance to whole-blood-mediated killing, and enhanced bacterial virulence in a mouse model of systemic infection. This study reveals CcpE as an important metabolic sensor that allows S. aureus to sense and adjust its metabolic state and subsequently to coordinate the expression of virulence factors and bacterial virulence.