Integrated pathway modules using time-course metabolic profiles and EST data from Milnesium tardigradum

BackgroundTardigrades are multicellular organisms, resistant to extreme environmental changes such as heat, drought, radiation and freezing. They outlast these conditions in an inactive form (tun) to escape damage to cellular structures and cell death. Tardigrades are apparently able to prevent or repair such damage and are therefore a crucial model organism for stress tolerance. Cultures of the tardigrade Milnesium tardigradum were dehydrated by removing the surrounding water to induce tun formation. During this process and the subsequent rehydration, metabolites were measured in a time series by GC-MS. Additionally expressed sequence tags are available, especially libraries generated from the active and inactive state. The aim of this integrated analysis is to trace changes in tardigrade metabolism and identify pathways responsible for their extreme resistance against physical stress.ResultsIn this study we propose a novel integrative approach for the analysis of metabolic networks to identify modules of joint shifts on the transcriptomic and metabolic levels. We derive a tardigrade-specific metabolic network represented as an undirected graph with 3,658 nodes (metabolites) and 4,378 edges (reactions). Time course metabolite profiles are used to score the network nodes showing a significant change over time. The edges are scored according to information on enzymes from the EST data. Using this combined information, we identify a key subnetwork (functional module) of concerted changes in metabolic pathways, specific for de- and rehydration. The module is enriched in reactions showing significant changes in metabolite levels and enzyme abundance during the transition. It resembles the cessation of a measurable metabolism (e.g. glycolysis and amino acid anabolism) during the tun formation, the production of storage metabolites and bioprotectants, such as DNA stabilizers, and the generation of amino acids and cellular components from monosaccharides as carbon and energy source during rehydration.ConclusionsThe functional module identifies relationships among changed metabolites (e.g. spermidine) and reactions and provides first insights into important altered metabolic pathways. With sparse and diverse data available, the presented integrated metabolite network approach is suitable to integrate all existing data and analyse it in a combined manner.

[1]  J. Clegg,et al.  Cryptobiosis--a peculiar state of biological organization. , 2001, Comparative biochemistry and physiology. Part B, Biochemistry & molecular biology.

[2]  M. Clark,et al.  Dormancy and Resistance in Harsh Environments , 2010 .

[3]  H. Köhler,et al.  Stress gene (hsp70) sequences and quantitative expression in Milnesium tardigradum (Tardigrada) during active and cryptobiotic stages , 2004, Journal of Experimental Biology.

[4]  Stan Pounds,et al.  Estimating the Occurrence of False Positives and False Negatives in Microarray Studies by Approximating and Partitioning the Empirical Distribution of P-values , 2003, Bioinform..

[5]  Jason A. Papin,et al.  Applications of genome-scale metabolic reconstructions , 2009, Molecular systems biology.

[6]  J. Hoheisel,et al.  Correspondence analysis applied to microarray data , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Ralph O. Schill,et al.  Anhydrobiotic Abilities of Tardigrades , 2010 .

[8]  Roberto Guidetti,et al.  Heat shock proteins in encysted and anhydrobiotic eutardigrades , 2012 .

[9]  Jan Hummel,et al.  Retention index thresholds for compound matching in GC-MS metabolite profiling. , 2008, Journal of chromatography. B, Analytical technologies in the biomedical and life sciences.

[10]  Joachim Selbig,et al.  Decision tree supported substructure prediction of metabolites from GC-MS profiles , 2010, Metabolomics.

[11]  P Green,et al.  Base-calling of automated sequencer traces using phred. II. Error probabilities. , 1998, Genome research.

[12]  Adam Godzik,et al.  Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences , 2006, Bioinform..

[13]  Steffen Hengherr,et al.  Trehalose and anhydrobiosis in tardigrades – evidence for divergence in responses to dehydration , 2008, The FEBS journal.

[14]  Benno Schwikowski,et al.  Discovering regulatory and signalling circuits in molecular interaction networks , 2002, ISMB.

[15]  Ron Shamir,et al.  Identification of functional modules using network topology and high-throughput data , 2007, BMC Systems Biology.

[16]  R. O. Schill,et al.  Induction of Hsp70 by desiccation, ionising radiation and heat-shock in the eutardigrade Richtersius coronifer. , 2007, Comparative biochemistry and physiology. Part B, Biochemistry & molecular biology.

[17]  J. Nielsen,et al.  Uncovering transcriptional regulation of metabolism by using metabolic network topology. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Frank Förster,et al.  Transcriptome survey of the anhydrobiotic tardigrade Milnesium tardigradum in comparison with Hypsibius dujardini and Richtersius coronifer , 2010, BMC Genomics.

[19]  G. A. Mack,et al.  K-Sample Rank Tests for Umbrella Alternatives , 1981 .

[20]  Tobias Müller,et al.  Identifying functional modules in protein–protein interaction networks: an integrated exact approach , 2008, ISMB.

[21]  Jean YH Yang,et al.  Bioconductor: open software development for computational biology and bioinformatics , 2004, Genome Biology.

[22]  Stefan Wiemann,et al.  KEGGgraph: a graph approach to KEGG PATHWAY in R and bioconductor , 2009, Bioinform..

[23]  Olga G. Troyanskaya,et al.  Integrated functional networks of process, tissue, and developmental stage specific interactions in Arabidopsis thaliana , 2010, BMC Systems Biology.

[24]  S. Panke,et al.  Putative regulatory sites unraveled by network-embedded thermodynamic analysis of metabolome data , 2006, Molecular systems biology.

[25]  A. R. Jonckheere,et al.  A DISTRIBUTION-FREE k-SAMPLE TEST AGAINST ORDERED ALTERNATIVES , 1954 .

[26]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[27]  J. Nielsen,et al.  Integration of metabolome data with metabolic networks reveals reporter reactions , 2006, Molecular systems biology.

[28]  Roded Sharan,et al.  Efficient Algorithms for Detecting Signaling Pathways in Protein Interaction Networks , 2006, J. Comput. Biol..

[29]  Alexander Erban,et al.  TagFinder for the quantitative analysis of gas chromatography - mass spectrometry (GC-MS)-based metabolite profiling experiments , 2008, Bioinform..

[30]  P. Dixon.  VEGAN, a package of R functions for community ecology , 2003 .

[31]  Jörg Bernhardt,et al.  Modeling antibiotic and cytotoxic effects of the dimeric isoquinoline IQ-143 on metabolism and its regulation in Staphylococcus aureus, Staphylococcus epidermidis and human cells , 2011, Genome Biology.

[32]  Frank Förster,et al.  Stress response in tardigrades: differential gene expression of molecular chaperones , 2010, Cell Stress and Chaperones.

[33]  T. J. Terpstra,et al.  The asymptotic normality and consistency of kendall's test against trend, when ties are present in one ranking , 1952 .

[34]  Luke Hunter,et al.  Interpreting Metabolomic Profiles using Unbiased Pathway Models , 2010, PLoS Comput. Biol..

[35]  O. Fiehn,et al.  Metabolite profiling for plant functional genomics , 2000, Nature Biotechnology.

[36]  Gilbert Saporta,et al.  L'analyse des données , 1981 .

[37]  B. Palsson,et al.  Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. , 2003, Genome research.

[38]  Olga G. Troyanskaya,et al.  Coordinated Concentration Changes of Transcripts and Metabolites in Saccharomyces cerevisiae , 2009, PLoS Comput. Biol..

[39]  Gábor Csárdi,et al.  The igraph software package for complex network research , 2006 .

[40]  C. Daub,et al.  BMC Systems Biology , 2007 .

[41]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[42]  P. Green,et al.  Base-calling of automated sequencer traces using phred. I. Accuracy assessment. , 1998, Genome research.

[43]  Tobias Müller,et al.  Robustness and accuracy of functional modules in integrated network analysis , 2012, Bioinform..

[44]  Markus J. Herrgård,et al.  Network-based prediction of human tissue-specific metabolism , 2008, Nature Biotechnology.

[45]  Peter Westh,et al.  Trehalose accumulation in the tardigrade Adorybiotus coronifer during anhydrobiosis , 1991 .

[46]  Akiyasu C. Yoshizawa,et al.  KAAS: an automatic genome annotation and pathway reconstruction server , 2007, Environmental health perspectives.

[47]  Robert Gentleman,et al.  Using GOstats to test gene lists for GO term association , 2007, Bioinform..

[48]  Jens Nielsen,et al.  Global Transcriptional and Physiological Responses of Saccharomyces cerevisiae to Ammonium, l-Alanine, or l-Glutamine Limitation , 2006, Applied and Environmental Microbiology.

[49]  Tobias Müller,et al.  Bioinformatics Applications Note Systems Biology Bionet: an R-package for the Functional Analysis of Biological Networks , 2022 .

[50]  Peter Westh,et al.  Cryptobiosis in the Eutardigrade Adorybiotus (Richtersius) coronifer: Tolerance to Alcohols, Temperature and de novo Protein Synthesis , 2001 .

[51]  Marcus Frohme,et al.  Investigating heat shock proteins of tardigrades in active versus anhydrobiotic state using shotgun proteomics , 2011 .

[52]  Frank Förster,et al.  Tardigrade workbench: comparing stress-related proteins, sequence-similar and functional protein clusters as well as RNA elements in tardigrades , 2009, BMC Genomics.

[53]  Frank Förster,et al.  Bioinformatics identifies tardigrade molecular adaptations including the DNA‐j family and first steps towards dynamical modelling , 2011 .

[54]  Joachim Kopka,et al.  Nonsupervised construction and application of mass spectral and retention time index libraries from time-of-flight gas chromatography-mass spectrometry metabolite profiles. , 2007, Methods in molecular biology.

[55]  K. Ingemar Jönsson,et al.  Trehalose in three species of desiccation tolerant tardigrades , 2010 .

[56]  Lazzaro Spallanzani.  Opuscoli de fisica animale e vegetabile , 1826 .