Sebastian Klie

Metabolomic and transcriptomic stress response of Escherichia coli.

Environmental fluctuations lead to a rapid adjustment of the physiology of Escherichia coli, necessitating changes on every level of the underlying cellular and molecular network. Thus far, the majority of global analyses of E. coli stress responses have been limited to just one level, gene expression. Here, we incorporate the metabolite composition together with gene expression data to provide a more comprehensive insight on system level stress adjustments by describing detailed time‐resolved E. coli response to five different perturbations (cold, heat, oxidative stress, lactose diauxie, and stationary phase). The metabolite response is more specific as compared with the general response observed on the transcript level and is reflected by much higher specificity during the early stress adaptation phase and when comparing the stationary phase response to other perturbations. Despite these differences, the response on both levels still follows the same dynamics and general strategy of energy conservation as reflected by rapid decrease of central carbon metabolism intermediates coinciding with downregulation of genes related to cell growth. Application of co‐clustering and canonical correlation analysis on combined metabolite and transcript data identified a number of significant condition‐dependent associations between metabolites and transcripts. The results confirm and extend existing models about co‐regulation between gene expression and metabolites demonstrating the power of integrated systems oriented analysis.

PlaNet: Combined Sequence and Expression Comparisons across Plant Networks Derived from Seven Species

Genes that are similarly expressed, or coexpressed, are often involved in related biological processes. Such coexpressed relationships also appear to be conserved across species. The PlaNet platform enables comparative analysis of genome-wide coexpression networks across seven plant species, thus enabling prediction of gene function and elucidation of the identity of functional homologs. The model organism Arabidopsis thaliana is readily used in basic research due to resource availability and relative speed of data acquisition. A major goal is to transfer acquired knowledge from Arabidopsis to crop species. However, the identification of functional equivalents of well-characterized Arabidopsis genes in other plants is a nontrivial task. It is well documented that transcriptionally coordinated genes tend to be functionally related and that such relationships may be conserved across different species and even kingdoms. To exploit such relationships, we constructed whole-genome coexpression networks for Arabidopsis and six important plant crop species. The interactive networks, clustered using the HCCA algorithm, are provided under the banner PlaNet (http://aranet.mpimp-golm.mpg.de). We implemented a comparative network algorithm that estimates similarities between network structures. Thus, the platform can be used to swiftly infer similar coexpressed network vicinities within and across species and can predict the identity of functional homologs. We exemplify this using the PSA-D and chalcone synthase-related gene networks. Finally, we assessed how ontology terms are transcriptionally connected in the seven species and provide the corresponding MapMan term coexpression networks. The data support the contention that this platform will considerably improve transfer of knowledge generated in Arabidopsis to valuable crop species.

High-density kinetic analysis of the metabolomic and transcriptomic response of Arabidopsis to eight environmental conditions

The time-resolved response of Arabidopsis thaliana towards changing light and/or temperature at the transcriptome and metabolome level is presented. Plants grown at 21°C with a light intensity of 150 μE m⁻² sec⁻¹ were either kept at this condition or transferred into seven different environments (4°C, darkness; 21°C, darkness; 32°C, darkness; 4°C, 85 μE m⁻² sec⁻¹; 21 °C, 75 μE m⁻² sec⁻¹; 21°C, 300 μE m⁻² sec⁻¹ ; 32°C, 150 μE m⁻² sec⁻¹). Samples were taken before (0 min) and at 22 time points after transfer resulting in (8×) 22 time points covering both a linear and a logarithmic time series totaling 177 states. Hierarchical cluster analysis shows that individual conditions (defined by temperature and light) diverge into distinct trajectories at condition-dependent times and that the metabolome follows different kinetics from the transcriptome. The metabolic responses are initially relatively faster when compared with the transcriptional responses. Gene Ontology over-representation analysis identifies a common response for all changed conditions at the transcriptome level during the early response phase (5-60 min). Metabolic networks reconstructed via metabolite-metabolite correlations reveal extensive environment-specific rewiring. Detailed analysis identifies conditional connections between amino acids and intermediates of the tricarboxylic acid cycle. Parallel analysis of transcriptional changes strongly support a model where in the absence of photosynthesis at normal/high temperatures protein degradation occurs rapidly and subsequent amino acid catabolism serves as the main cellular energy supply. These results thus demonstrate the engagement of the electron transfer flavoprotein system under short-term environmental perturbations.

The Choice between MapMan and Gene Ontology for Automated Gene Function Prediction in Plant Science

Since the introduction of the Gene Ontology (GO), the analysis of high-throughput data has become tightly coupled with the use of ontologies to establish associations between knowledge and data in an automated fashion. Ontologies provide a systematic description of knowledge by a controlled vocabulary of defined structure in which ontological concepts are connected by pre-defined relationships. In plant science, MapMan and GO offer two alternatives for ontology-driven analyses. Unlike GO, initially developed to characterize microbial systems, MapMan was specifically designed to cover plant-specific pathways and processes. While the dependencies between concepts in MapMan are modeled as a tree, in GO these are captured in a directed acyclic graph. Therefore, the difference in ontologies may cause discrepancies in data reduction, visualization, and hypothesis generation. Here provide the first systematic comparative analysis of GO and MapMan for the case of the model plant species Arabidopsis thaliana (Arabidopsis) with respect to their structural properties and difference in distributions of information content. In addition, we investigate the effect of the two ontologies on the specificity and sensitivity of automated gene function prediction via the coupling of co-expression networks and the guilt-by-association principle. Automated gene function prediction is particularly needed for the model plant Arabidopsis in which only half of genes have been functionally annotated based on sequence similarity to known genes. The results highlight the need for structured representation of species-specific biological knowledge, and warrants caution in the design principles employed in future ontologies.

Decreased Nucleotide and Expression Diversity and Modified Coexpression Patterns Characterize Domestication in the Common Bean.

About 60% of the nucleotide diversity was lost during domestication of the common bean. The whole pattern of gene expression has also been affected, with changes in the patterns of coexpression among genes and 18% reduction in the overall diversity of gene expression. About 9% of the genes were selected during domestication, which is associated with further reduced diversity of expression. Using RNA sequencing technology and de novo transcriptome assembly, we compared representative sets of wild and domesticated accessions of common bean (Phaseolus vulgaris) from Mesoamerica. RNA was extracted at the first true-leaf stage, and de novo assembly was used to develop a reference transcriptome; the final data set consists of ∼190,000 single nucleotide polymorphisms from 27,243 contigs in expressed genomic regions. A drastic reduction in nucleotide diversity (∼60%) is evident for the domesticated form, compared with the wild form, and almost 50% of the contigs that are polymorphic were brought to fixation by domestication. In parallel, the effects of domestication decreased the diversity of gene expression (18%). While the coexpression networks for the wild and domesticated accessions demonstrate similar seminal network properties, they show distinct community structures that are enriched for different molecular functions. After simulating the demographic dynamics during domestication, we found that 9% of the genes were actively selected during domestication. We also show that selection induced a further reduction in the diversity of gene expression (26%) and was associated with 5-fold enrichment of differentially expressed genes. While there is substantial evidence of positive selection associated with domestication, in a few cases, this selection has increased the nucleotide diversity in the domesticated pool at target loci associated with abiotic stress responses, flowering time, and morphology.

Analysis of the compartmentalized metabolome – a validation of the non-aqueous fractionation technique

With the development of high-throughput metabolic technologies, a plethora of primary and secondary compounds have been detected in the plant cell. However, there are still major gaps in our understanding of the plant metabolome. This is especially true with regards to the compartmental localization of these identified metabolites. Non-aqueous fractionation (NAF) is a powerful technique for the determination of subcellular metabolite distributions in eukaryotic cells, and it has become the method of choice to analyze the distribution of a large number of metabolites concurrently. However, the NAF technique produces a continuous gradient of metabolite distributions, not discrete assignments. Resolution of these distributions requires computational analyses based on marker molecules to resolve compartmental localizations. In this article we focus on expanding the computational analysis of data derived from NAF. Along with an experimental workflow, we describe the critical steps in NAF experiments and how computational approaches can aid in assessing the quality and robustness of the derived data. For this, we have developed and provide a new version (v1.2) of the BestFit command line tool for calculation and evaluation of subcellular metabolite distributions. Furthermore, using both simulated and experimental data we show the influence on estimated subcellular distributions by modulating important parameters, such as the number of fractions taken or which marker molecule is selected. Finally, we discuss caveats and benefits of NAF analysis in the context of the compartmentalized metabolome.

Conserved Changes in the Dynamics of Metabolic Processes during Fruit Development and Ripening across Species

STATIS analysis of metabolic profiles discerns conserved metabolic processes during development and ripening across fruits. Computational analyses of molecular phenotypes traditionally aim at identifying biochemical components that exhibit differential expression under various scenarios (e.g. environmental and internal perturbations) in a single species. High-throughput metabolomics technologies allow the quantification of (relative) metabolite levels across developmental stages in different tissues, organs, and species. Novel methods for analyzing the resulting multiple data tables could reveal preserved dynamics of metabolic processes across species. The problem we address in this study is 2-fold. (1) We derive a single data table, referred to as a compromise, which captures information common to the investigated set of multiple tables containing data on different fruit development and ripening stages in three climacteric (i.e. peach [Prunus persica] and two tomato [Solanum lycopersicum] cultivars, Ailsa Craig and M82) and two nonclimacteric (i.e. strawberry [Fragaria × ananassa] and pepper [Capsicum chilense]) fruits; in addition, we demonstrate the power of the method to discern similarities and differences between multiple tables by analyzing publicly available metabolomics data from three tomato ripening mutants together with two tomato cultivars. (2) We identify the conserved dynamics of metabolic processes, reflected in the data profiles of the corresponding metabolites that contribute most to the determined compromise. Our analysis is based on an extension to principal component analysis, called STATIS, in combination with pathway overenrichment analysis. Based on publicly available metabolic profiles for the investigated species, we demonstrate that STATIS can be used to identify the metabolic processes whose behavior is similarly affected during fruit development and ripening. These findings ultimately provide insights into the pathways that are essential during fruit development and ripening across species.

Integration of transcriptomics and metabolomics data specifies the metabolic response of Chlamydomonas to rapamycin treatment

Flux phenotypes predicted by constraint-based methods can be refined by the inclusion of heterogeneous data. While recent advances facilitate the integration of transcriptomics and proteomics data, purely stoichiometry-based approaches for the prediction of flux phenotypes by considering metabolomics data are lacking. Here we propose a constraint-based method, termed TREM-Flux, for integrating time-resolved metabolomics and transcriptomics data. We demonstrate the applicability of TREM-Flux in the dissection of the metabolic response of Chlamydomonas reinhardtii to rapamycin treatment by integrating the expression levels of 982 genes and the content of 45 metabolites obtained from two growth conditions. The findings pinpoint cysteine and methionine metabolism to be most affected by the rapamycin treatment. Our study shows that the integration of time-resolved unlabeled metabolomics data in addition to transcriptomics data can specify the metabolic pathways involved in the systems response to a studied treatment.

Principal Components Analysis

Principal components analysis (PCA) is a standard tool in multivariate data analysis to reduce the number of dimensions, while retaining as much as possible of the datas variation. Instead of investigating thousands of original variables, the first few components containing the majority of the datas variation are explored. The visualization and statistical analysis of these new variables, the principal components, can help to find similarities and differences between samples. Important original variables that are the major contributors to the first few components can be discovered as well.This chapter seeks to deliver a conceptual understanding of PCA as well as a mathematical description. We describe how PCA can be used to analyze different datasets, and we include practical code examples. Possible shortcomings of the methodology and ways to overcome these problems are also discussed.

FamNet: A Framework to Identify Multiplied Modules Driving Pathway Expansion in Plants.

Gene module multiplication drives pathway expansion in plants. Gene duplications generate new genes that can acquire similar but often diversified functions. Recent studies of gene coexpression networks have indicated that, not only genes, but also pathways can be multiplied and diversified to perform related functions in different parts of an organism. Identification of such diversified pathways, or modules, is needed to expand our knowledge of biological processes in plants and to understand how biological functions evolve. However, systematic explorations of modules remain scarce, and no user-friendly platform to identify them exists. We have established a statistical framework to identify modules and show that approximately one-third of the genes of a plant’s genome participate in hundreds of multiplied modules. Using this framework as a basis, we implemented a platform that can explore and visualize multiplied modules in coexpression networks of eight plant species. To validate the usefulness of the platform, we identified and functionally characterized pollen- and root-specific cell wall modules that multiplied to confer tip growth in pollen tubes and root hairs, respectively. Furthermore, we identified multiplied modules involved in secondary metabolite synthesis and corroborated them by metabolite profiling of tobacco (Nicotiana tabacum) tissues. The interactive platform, referred to as FamNet, is available at http://www.gene2function.de/famnet.html.