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.

Genomic insights into host and parasite interactions during intracellular infection by Toxoplasma gondii.

To gain insights into the molecular interactions of an intracellular pathogen and its host cell, we studied the gene expression and chromatin states of human fibroblasts infected with the Apicomplexan parasite Toxoplasma gondii. We show a striking activation of host cell genes that regulate a number of cellular processes, some of which are protective of the host cell, others likely to be advantageous to the pathogen. The simultaneous capture of host and parasite genomic information allowed us to gain insights into the regulation of the T. gondii genome. We show how chromatin accessibility and transcriptional profiling together permit novel annotation of the parasite's genome, including more accurate mapping of known genes and the identification of new genes and cis-regulatory elements. Motif analysis reveals not only the known T. gondii AP2 transcription factor-binding site but also a previously-undiscovered candidate TATA box-containing motif at one-quarter of promoters. By inferring the transcription factor and upstream cell signaling responses involved in the host cell, we can use genomic information to gain insights into T. gondii's perturbation of host cell physiology. Our resulting model builds on previously-described human host cell signalling responses to T. gondii infection, linked to induction of specific transcription factors, some of which appear to be solely protective of the host cell, others of which appear to be co-opted by the pathogen to enhance its own survival.

toxoMine: an integrated omics data warehouse for Toxoplasma gondii systems biology research.

Toxoplasma gondii (T. gondii) is an obligate intracellular parasite that must monitor for changes in the host environment and respond accordingly; however, it is still not fully known which genetic or epigenetic factors are involved in regulating virulence traits of T. gondii. There are on-going efforts to elucidate the mechanisms regulating the stage transition process via the application of high-throughput epigenomics, genomics and proteomics techniques. Given the range of experimental conditions and the typical yield from such high-throughput techniques, a new challenge arises: how to effectively collect, organize and disseminate the generated data for subsequent data analysis. Here, we describe toxoMine, which provides a powerful interface to support sophisticated integrative exploration of high-throughput experimental data and metadata, providing researchers with a more tractable means toward understanding how genetic and/or epigenetic factors play a coordinated role in determining pathogenicity of T. gondii. As a data warehouse, toxoMine allows integration of high-throughput data sets with public T. gondii data. toxoMine is also able to execute complex queries involving multiple data sets with straightforward user interaction. Furthermore, toxoMine allows users to define their own parameters during the search process that gives users near-limitless search and query capabilities. The interoperability feature also allows users to query and examine data available in other InterMine systems, which would effectively augment the search scope beyond what is available to toxoMine. toxoMine complements the major community database ToxoDB by providing a data warehouse that enables more extensive integrative studies for T. gondii. Given all these factors, we believe it will become an indispensable resource to the greater infectious disease research community.

DNA repair in species with extreme lifespan differences.

Differences in DNA repair capacity have been hypothesized to underlie the great range of maximum lifespans among mammals. However, measurements of individual DNA repair activities in cells and animals have not substantiated such a relationship because utilization of repair pathways among animals--depending on habitats, anatomical characteristics, and life styles--varies greatly between mammalian species. Recent advances in high-throughput genomics, in combination with increased knowledge of the genetic pathways involved in genome maintenance, now enable a comprehensive comparison of DNA repair transcriptomes in animal species with extreme lifespan differences. Here we compare transcriptomes of liver, an organ with high oxidative metabolism and abundant spontaneous DNA damage, from humans, naked mole rats, and mice, with maximum lifespans of ~120, 30, and 3 years, respectively, with a focus on genes involved in DNA repair. The results show that the longer-lived species, human and naked mole rat, share higher expression of DNA repair genes, including core genes in several DNA repair pathways. A more systematic approach of signaling pathway analysis indicates statistically significant upregulation of several DNA repair signaling pathways in human and naked mole rat compared with mouse. The results of this present work indicate, for the first time, that DNA repair is upregulated in a major metabolic organ in long-lived humans and naked mole rats compared with short-lived mice. These results strongly suggest that DNA repair can be considered a genuine longevity assurance system.

Distinct strains of Toxoplasma gondii feature divergent transcriptomes regardless of developmental stage.

Using high through-put RNA sequencing, we assayed the transcriptomes of three different strains of Toxoplasma gondii representing three common genotypes under both in vitro tachyzoite and in vitro bradyzoite-inducing alkaline stress culture conditions. Strikingly, the differences in transcriptional profiles between the strains, RH, PLK, and CTG, is much greater than differences between tachyzoites and alkaline stressed in vitro bradyzoites. With an FDR of 10%, we identified 241 genes differentially expressed between CTG tachyzoites and in vitro bradyzoites, including 5 putative AP2 transcription factors. We also observed a close association between cell cycle regulated genes and differentiation. By Gene Set Enrichment Analysis (GSEA), there are a number of KEGG pathways associated with the in vitro bradyzoite transcriptomes of PLK and CTG, including pyrimidine metabolism and DNA replication. These functions are likely associated with cell-cycle arrest. When comparing mRNA levels between strains, we identified 1,526 genes that were differentially expressed regardless of culture-condition as well as 846 differentially expressed only in bradyzoites and 542 differentially expressed only in tachyzoites between at least two strains. Using GSEA, we identified that ribosomal proteins were expressed at significantly higher levels in the CTG strain than in either the RH or PLK strains. This association holds true regardless of life cycle stage.