Henning Redestig

Extension of the visualization tool mapman to allow statistical analysis of arrays, display of coresponding genes, and comparison with known responses

MapMan is a user-driven tool that displays large genomics datasets onto diagrams of metabolic pathways or other processes. Here, we present new developments, including improvements of the gene assignments and the user interface, a strategy to visualize multilayered datasets, the incorporation of statistics packages, and extensions of the software to incorporate more biological information including visualization of coresponding genes and horizontal searches for similar global responses across large numbers of arrays.

pcaMethods—a bioconductor package providing PCA methods for incomplete data

UNLABELLED pcaMethods is a Bioconductor compliant library for computing principal component analysis (PCA) on incomplete data sets. The results can be analyzed directly or used to estimate missing values to enable the use of missing value sensitive statistical methods. The package was mainly developed with microarray and metabolite data sets in mind, but can be applied to any other incomplete data set as well. AVAILABILITY http://www.bioconductor.org

PageMan: An interactive ontology tool to generate, display, and annotate overview graphs for profiling experiments

BackgroundMicroarray technology has become a widely accepted and standardized tool in biology. The first microarray data analysis programs were developed to support pair-wise comparison. However, as microarray experiments have become more routine, large scale experiments have become more common, which investigate multiple time points or sets of mutants or transgenics. To extract biological information from such high-throughput expression data, it is necessary to develop efficient analytical platforms, which combine manually curated gene ontologies with efficient visualization and navigation tools. Currently, most tools focus on a few limited biological aspects, rather than offering a holistic, integrated analysis.ResultsHere we introduce PageMan, a multiplatform, user-friendly, and stand-alone software tool that annotates, investigates, and condenses high-throughput microarray data in the context of functional ontologies. It includes a GUI tool to transform different ontologies into a suitable format, enabling the user to compare and choose between different ontologies. It is equipped with several statistical modules for data analysis, including over-representation analysis and Wilcoxon statistical testing. Results are exported in a graphical format for direct use, or for further editing in graphics programs.PageMan provides a fast overview of single treatments, allows genome-level responses to be compared across several microarray experiments covering, for example, stress responses at multiple time points. This aids in searching for trait-specific changes in pathways using mutants or transgenics, analyzing development time-courses, and comparison between species. In a case study, we analyze the results of publicly available microarrays of multiple cold stress experiments using PageMan, and compare the results to a previously published meta-analysis.PageMan offers a complete users guide, a web-based over-representation analysis as well as a tutorial, and is freely available at http://mapman.mpimp-golm.mpg.de/pageman/.ConclusionPageMan allows multiple microarray experiments to be efficiently condensed into a single page graphical display. The flexible interface allows data to be quickly and easily visualized, facilitating comparisons within experiments and to published experiments, thus enabling researchers to gain a rapid overview of the biological responses in the experiments.

The Metabolic Response of Heterotrophic Arabidopsis Cells to Oxidative Stress

To cope with oxidative stress, the metabolic network of plant cells must be reconfigured either to bypass damaged enzymes or to support adaptive responses. To characterize the dynamics of metabolic change during oxidative stress, heterotrophic Arabidopsis (Arabidopsis thaliana) cells were treated with menadione and changes in metabolite abundance and 13C-labeling kinetics were quantified in a time series of samples taken over a 6 h period. Oxidative stress had a profound effect on the central metabolic pathways with extensive metabolic inhibition radiating from the tricarboxylic acid cycle and including large sectors of amino acid metabolism. Sequential accumulation of metabolites in specific pathways indicated a subsequent backing up of glycolysis and a diversion of carbon into the oxidative pentose phosphate pathway. Microarray analysis revealed a coordinated transcriptomic response that represents an emergency coping strategy allowing the cell to survive the metabolic hiatus. Rather than attempt to replace inhibited enzymes, transcripts encoding these enzymes are in fact down-regulated while an antioxidant defense response is mounted. In addition, a major switch from anabolic to catabolic metabolism is signaled. Metabolism is also reconfigured to bypass damaged steps (e.g. induction of an external NADH dehydrogenase of the mitochondrial respiratory chain). The overall metabolic response of Arabidopsis cells to oxidative stress is remarkably similar to the superoxide and hydrogen peroxide stimulons of bacteria and yeast (Saccharomyces cerevisiae), suggesting that the stress regulatory and signaling pathways of plants and microbes may share common elements.

Identification of metabolic and biomass QTL in Arabidopsis thaliana in a parallel analysis of RIL and IL populations.

Plant growth and development are tightly linked to primary metabolism and are subject to natural variation. In order to obtain an insight into the genetic factors controlling biomass and primary metabolism and to determine their relationships, two Arabidopsis thaliana populations [429 recombinant inbred lines (RIL) and 97 introgression lines (IL), derived from accessions Col-0 and C24] were analyzed with respect to biomass and metabolic composition using a mass spectrometry-based metabolic profiling approach. Six and 157 quantitative trait loci (QTL) were identified for biomass and metabolic content, respectively. Two biomass QTL coincide with significantly more metabolic QTL (mQTL) than statistically expected, supporting the notion that the metabolic profile and biomass accumulation of a plant are linked. On the same basis, three out the six biomass QTL can be simulated purely on the basis of metabolic composition. QTL based on analysis of the introgression lines were in substantial agreement with the RIL-based results: five of six biomass QTL and 55% of the mQTL found in the RIL population were also found in the IL population at a significance level of P ≤ 0.05, with >80% agreement on the allele effects. Some of the differences could be attributed to epistatic interactions. Depending on the search conditions, metabolic pathway-derived candidate genes were found for 24–67% of all tested mQTL in the database AraCyc 3.5. This dataset thus provides a comprehensive basis for the detection of functionally relevant variation in known genes with metabolic function and for identification of genes with hitherto unknown roles in the control of metabolism.

The OMA orthology database in 2015: function predictions, better plant support, synteny view and other improvements

The Orthologous Matrix (OMA) project is a method and associated database inferring evolutionary relationships amongst currently 1706 complete proteomes (i.e. the protein sequence associated for every protein-coding gene in all genomes). In this update article, we present six major new developments in OMA: (i) a new web interface; (ii) Gene Ontology function predictions as part of the OMA pipeline; (iii) better support for plant genomes and in particular homeologs in the wheat genome; (iv) a new synteny viewer providing the genomic context of orthologs; (v) statically computed hierarchical orthologous groups subsets downloadable in OrthoXML format; and (vi) possibility to export parts of the all-against-all computations and to combine them with custom data for ‘client-side’ orthology prediction. OMA can be accessed through the OMA Browser and various programmatic interfaces at http://omabrowser.org.

TargetSearch - a Bioconductor package for the efficient preprocessing of GC-MS metabolite profiling data

BackgroundMetabolite profiling, the simultaneous quantification of multiple metabolites in an experiment, is becoming increasingly popular, particularly with the rise of systems-level biology. The workhorse in this field is gas-chromatography hyphenated with mass spectrometry (GC-MS). The high-throughput of this technology coupled with a demand for large experiments has led to data pre-processing, i.e. the quantification of metabolites across samples, becoming a major bottleneck. Existing software has several limitations, including restricted maximum sample size, systematic errors and low flexibility. However, the biggest limitation is that the resulting data usually require extensive hand-curation, which is subjective and can typically take several days to weeks.ResultsWe introduce the TargetSearch package, an open source tool which is a flexible and accurate method for pre-processing even very large numbers of GC-MS samples within hours. We developed a novel strategy to iteratively correct and update retention time indices for searching and identifying metabolites. The package is written in the R programming language with computationally intensive functions written in C for speed and performance. The package includes a graphical user interface to allow easy use by those unfamiliar with R.ConclusionsTargetSearch allows fast and accurate data pre-processing for GC-MS experiments and overcomes the sample number limitations and manual curation requirements of existing software. We validate our method by carrying out an analysis against both a set of known chemical standard mixtures and of a biological experiment. In addition we demonstrate its capabilities and speed by comparing it with other GC-MS pre-processing tools. We believe this package will greatly ease current bottlenecks and facilitate the analysis of metabolic profiling data.

Metabolomic correlation-network modules in Arabidopsis based on a graph-clustering approach

BackgroundDeciphering the metabolome is essential for a better understanding of the cellular metabolism as a system. Typical metabolomics data show a few but significant correlations among metabolite levels when data sampling is repeated across individuals grown under strictly controlled conditions. Although several studies have assessed topologies in metabolomic correlation networks, it remains unclear whether highly connected metabolites in these networks have specific functions in known tissue- and/or genotype-dependent biochemical pathways.ResultsIn our study of metabolite profiles we subjected root tissues to gas chromatography-time-of-flight/mass spectrometry (GC-TOF/MS) and used published information on the aerial parts of 3 Arabidopsis genotypes, Col-0 wild-type, methionine over-accumulation 1 (mto1), and transparent testa4 (tt4) to compare systematically the metabolomic correlations in samples of roots and aerial parts. We then applied graph clustering to the constructed correlation networks to extract densely connected metabolites and evaluated the clusters by biochemical-pathway enrichment analysis. We found that the number of significant correlations varied by tissue and genotype and that the obtained clusters were significantly enriched for metabolites included in biochemical pathways.ConclusionsWe demonstrate that the graph-clustering approach identifies tissue- and/or genotype-dependent metabolomic clusters related to the biochemical pathway. Metabolomic correlations complement information about changes in mean metabolite levels and may help to elucidate the organization of metabolically functional modules.

Metabolomic approaches toward understanding nitrogen metabolism in plants

Plants can assimilate inorganic nitrogen (N) sources to organic N such as amino acids. N is the most important of the mineral nutrients required by plants and its metabolism is tightly coordinated with carbon (C) metabolism in the fundamental processes that permit plant growth. Increased understanding of N regulation may provide important insights for plant growth and improvement of quality of crops and vegetables because N as well as C metabolism are fundamental components of plant life. Metabolomics is a global biochemical approach useful to study N metabolism because metabolites not only reflect the ultimate phenotypes (traits), but can mediate transcript levels as well as protein levels directly and/or indirectly under different N conditions. This review outlines analytical and bioinformatic techniques particularly used to perform metabolomics for studying N metabolism in higher plants. Examples are used to illustrate the application of metabolomic techniques to the model plants Arabidopsis and rice, as well as other crops and vegetables.

Compensation for systematic cross-contribution improves normalization of mass spectrometry based metabolomics data.

Most mass spectrometry based metabolomics studies are semiquantitative and depend on efficient normalization techniques to suppress systematic error. A common approach is to include isotope-labeled internal standards (ISs) and then express the estimated metabolite abundances relative to the IS. Because of problems such as insufficient chromatographic resolution, however, the analytes may directly influence estimates of the IS, a phenomenon known as cross-contribution (CC). Normalization using ISs that suffer from CC effects will cause significant loss of information if the interfering analytes are associated with the studied factors. We present a novel normalization algorithm, which compensates for systematic CC effects that can be traced back to a linear association with the experimental design. The proposed method was found to be superior at purifying the signal of interest compared to current normalization methods when applied to two biological data sets and a multicomponent dilution mixture. Our method is applicable to data from randomized and designed experiments that use ISs to monitor the systematic error.