Gene set analysis (GSA) is used to elucidate genome-wide data, in particular transcriptome data. A multitude of methods have been proposed for this step of the analysis, and many of them have been compared and evaluated. Unfortunately, there is no consolidated opinion regarding what methods should be preferred, and the variety of available GSA software and implementations pose a difficulty for the end-user who wants to try out different methods. To address this, we have developed the R package Piano that collects a range of GSA methods into the same system, for the benefit of the end-user. Further on we refine the GSA workflow by using modifications of the gene-level statistics. This enables us to divide the resulting gene set P -values into three classes, describing different aspects of gene expression directionality at gene set level. We use our fully implemented workflow to investigate the impact of the individual components of GSA by using microarray and RNA-seq data. The results show that the evaluated methods are globally similar and the major separation correlates well with our defined directionality classes. As a consequence of this, we suggest to use a consensus scoring approach, based on multiple GSA runs. In combination with the directionality classes, this constitutes a more thorough basis for an enriched biological interpretation.

Background: The intestinal epithelium separates us from the microbiota but also interacts with it and thus affects host immune status and physiology. Previous studies investigated microbiota-induced responses in the gut using intact tissues or unfractionated epithelial cells, thereby limiting conclusions about regional differences in the epithelium. Here, we sought to investigate microbiota-induced transcriptional responses in specific fractions of intestinal epithelial cells. To this end, we used microarray analysis of laser capture microdissection (LCM)-harvested ileal and colonic tip and crypt epithelial fractions from germ-free and conventionally raised mice and from mice during the time course of colonization. Results: We found that about 10% of the host’s transcriptome was microbially regulated, mainly including genes annotated with functions in immunity, cell proliferation, and metabolism. The microbial impact on host gene expression was highly site specific, as epithelial responses to the microbiota differed between cell fractions. Specific transcriptional regulators were enriched in each fraction. In general, the gut microbiota induced a more rapid response in the colon than in the ileum. Conclusions: Our study indicates that the microbiota engage different regulatory networks to alter host gene expression in a particular niche. Understanding host-microbiota interactions on a cellular level may facilitate signaling pathways that contribute to health and disease and thus provide new therapeutic strategies.

Background: The genome-scale metabolic model of Saccharomyces cerevisiae , first presented in 2003, was the first genome-scale network reconstruction for a eukaryotic organism. Since then continuous efforts have been made in order to improve and expand the yeast metabolic network. Results: Here we present iTO977, a comprehensive genome-scale metabolic model that contains more reactions, metabolites and genes than previous models. The model was constructed based on two earlier reconstructions, namely iIN800 and the consensus network, and then improved and expanded using gap-filling methods and by introducing new reactions and pathways based on studies of the literature and databases. The model was shown to perform well both for growth simulations in different media and gene essentiality analysis for single and double knock-outs. Further, the model was used as a scaffold for integrating transcriptomics, and flux data from four different conditions in order to identify transcriptionally controlled reactions, i.e. reactions that change both in flux and transcription between the compared conditions. Conclusion: We present a new yeast model that represents a comprehensive up-to-date collection of knowledge on yeast metabolism. The model was used for simulating the yeast metabolism under four different growth conditions and experimental data from these four conditions was integrated to the model. The model together with experimental data is a useful tool to identify condition-dependent changes of metabolism between different environmental conditions.

Background: Pichia stipitis and Pichia pastoris have long been investigated due to their native abilities to metabolize every sugar from lignocellulose and to modulate methanol consumption, respectively. The latter has been driving the production of several recombinant proteins. As a result, significant advances in their biochemical knowledge, as well as in genetic engineering and fermentation methods have been generated. The release of their genome sequences has allowed systems level research. Results: In this work, genome-scale metabolic models (GEMs) of P. stipitis (iSS884) and P. pastoris (iLC915) were reconstructed. iSS884 includes 1332 reactions, 922 metabolites, and 4 compartments. iLC915 contains 1423 reactions, 899 metabolites, and 7 compartments. Compared with the previous GEMs of P. pastoris, PpaMBEL1254 and iPP668, iLC915 contains more genes and metabolic functions, as well as improved predictive capabilities. Simulations of physiological responses for the growth of both yeasts on selected carbon sources using iSS884 and iLC915 closely reproduced the experimental data. Additionally, the iSS884 model was used to predict ethanol production from xylose at different oxygen uptake rates. Simulations with iLC915 closely reproduced the effect of oxygen uptake rate on physiological states of P. pastoris expressing a recombinant protein. The potential of P. stipitis for the conversion of xylose and glucose into ethanol using reactors in series, and of P. pastoris to produce recombinant proteins using mixtures of methanol and glycerol or sorbitol are also discussed. Conclusions: In conclusion the first GEM of P. stipitis (iSS884) was reconstructed and validated. The expanded version of the P. pastoris GEM, iLC915, is more complete and has improved capabilities over the existing models. Both GEMs are useful frameworks to explore the versatility of these yeasts and to capitalize on their biotechnological potentials.

The gut microbiota plays an important role in human health and disease by acting as a metabolic organ. Metagenomic sequencing has shown how dysbiosis in the gut microbiota is associated with human metabolic diseases such as obesity and diabetes. Modeling may assist to gain insight into the metabolic implication of an altered microbiota. Fast and accurate reconstruction of metabolic models for members of the gut microbiota, as well as methods to simulate a community of microorganisms, are therefore needed. The Integrated Microbial Genomes (IMG) database contains functional annotation for nearly 4,650 bacterial genomes. This tremendous new genomic information adds new opportunities for systems biology to reconstruct accurate genome scale metabolic models (GEMs).