Enriching the gene set analysis of genome-wide data by incorporating directionality of gene expression and combining statistical hypotheses and methods

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

BioMet Toolbox 2.0: genome-wide analysis of metabolism and omics data

Analysis of large data sets using computational and mathematical tools have become a central part of biological sciences. Large amounts of data are being generated each year from different biological research fields leading to a constant development of software and algorithms aimed to deal with the increasing creation of information. The BioMet Toolbox 2.0 integrates a number of functionalities in a user-friendly environment enabling the user to work with biological data in a web interface. The unique and distinguishing feature of the BioMet Toolbox 2.0 is to provide a web user interface to tools for metabolic pathways and omics analysis developed under different platform-dependent environments enabling easy access to these computational tools

Proteome- and Transcriptome-Driven Reconstruction of the Human Myocyte Metabolic Network and Its Use for Identification of Markers for Diabetes

Skeletal myocytes are metabolically active and susceptible to insulin resistance and are thus implicated in type 2 diabetes (T2D). This complex disease involves systemic metabolic changes, and their elucidation at the systems level requires genome-wide data and biological networks. Genome-scale metabolic models (GEMs) provide a network context for the integration of high-throughput data. We generated myocyte-specific RNA-sequencing data and investigated their correlation with proteome data. These data were then used to reconstruct a comprehensive myocyte GEM. Next, we performed a meta-analysis of six studies comparing muscle transcription in T2D versus healthy subjects. Transcriptional changes were mapped on the myocyte GEM, revealing extensive transcriptional regulation in T2D, particularly around pyruvate oxidation, branched-chain amino acid catabolism, and tetrahydrofolate metabolism, connected through the downregulated dihydrolipoamide dehydrogenase. Strikingly, the gene signature underlying this metabolic regulation successfully classifies the disease state of individual samples, suggesting that regulation of these pathways is a ubiquitous feature of myocytes in response to T2D

Six Tissue Transcriptomics Reveals Specific Immune Suppression in Spleen by Dietary Polyunsaturated Fatty Acids

Dietary polyunsaturated fatty acids (PUFA) are suggested to modulate immune function, but the effects of dietary fatty acids composition on gene expression patterns in immune organs have not been fully characterized. In the current study we investigated how dietary fatty acids composition affects the total transcriptome profile, and especially, immune related genes in two immune organs, spleen (SPL) and bone marrow cells (BMC). Four tissues with metabolic function, skeletal muscle (SKM), white adipose tissue (WAT), brown adipose tissue (BAT), and liver (LIV), were investigated as a comparison. Following 8 weeks on low fat diet (LFD), high fat diet (HFD) rich in saturated fatty acids (HFD-S), or HFD rich in PUFA (HFD-P), tissue transcriptomics were analyzed by microarray and metabolic health assessed by fasting blood glucose level, HOMA-IR index, oral glucose tolerance test as well as quantification of crown-like structures in WAT. HFD-P corrected the metabolic phenotype induced by HFD-S. Interestingly, SKM and BMC were relatively inert to the diets, whereas the two adipose tissues (WAT and BAT) were mainly affected by HFD per se (both HFD-S and HFD-P). In particular, WAT gene expression was driven closer to that of the immune organs SPL and BMC by HFDs. The LIV exhibited different responses to both of the HFDs. Surprisingly, the spleen showed a major response to HFD-P (82 genes differed from LFD, mostly immune genes), while it was not affected at all by HFD-S (0 genes differed from LFD). In conclusion, the quantity and composition of dietary fatty acids affected the transcriptome in distinct manners in different organs. Remarkably, dietary PUFA, but not saturated fat, prompted a specific regulation of immune related genes in the spleen, opening the possibility that PUFA can regulate immune function by influencing gene expression in this organ

Diabetes causes marked inhibition of mitochondrial metabolism in pancreatic β-cells

Diabetes is a global health problem caused primarily by the inability of pancreatic β-cells to secrete adequate levels of insulin. The molecular mechanisms underlying the progressive failure of β-cells to respond to glucose in type-2 diabetes remain unresolved. Using a combination of transcriptomics and proteomics, we find significant dysregulation of major metabolic pathways in islets of diabetic βV59M mice, a non-obese, eulipidaemic diabetes model. Multiple genes/proteins involved in glycolysis/gluconeogenesis are upregulated, whereas those involved in oxidative phosphorylation are downregulated. In isolated islets, glucose-induced increases in NADH and ATP are impaired and both oxidative and glycolytic glucose metabolism are reduced. INS-1 β-cells cultured chronically at high glucose show similar changes in protein expression and reduced glucose-stimulated oxygen consumption: targeted metabolomics reveals impaired metabolism. These data indicate hyperglycaemia induces metabolic changes in β-cells that markedly reduce mitochondrial metabolism and ATP synthesis. We propose this underlies the progressive failure of β-cells in diabetes.Peer reviewe

Signatures of cell death and proliferation in perturbation transcriptomics data-from confounding factor to effective prediction

Transcriptional perturbation signatures are valuable data sources for functional genomics. Linking perturbation signatures to screenings opens the possibility to model cellular phenotypes from expression data and to identify efficacious drugs. We linked perturbation transcriptomics data from the LINCS-L1000 project with cell viability information upon genetic (Achilles project) and chemical (CTRP screen) perturbations yielding more than 90 000 signature-viability pairs. An integrated analysis showed that the cell viability signature is a major factor underlying perturbation signatures. The signature is linked to transcription factors regulating cell death, proliferation and division time. We used the cell viability-signature relationship to predict viability from transcriptomics signatures, and identified and validated compounds that induce cell death in tumor cell lines. We showed that cellular toxicity can lead to unexpected similarity of signatures, confounding mechanism of action discovery. Consensus compound signatures predicted cell-specific drug sensitivity, even if the signature is not measured in the same cell line, and outperformed conventional drug-specific features. Our results can help in understanding mechanisms behind cell death and removing confounding factors of transcriptomic perturbation screens. To interactively browse our results and predict cell viability in new gene expression samples, we developed CEVIChE (CEll VIability Calculator from gene Expression; https://saezlab.shinyapps.io/ceviche/)

The fecal metabolome as a functional readout of the gut microbiome

The human gut microbiome plays a key role in human health1, but 16S characterization lacks quantitative functional annotation2. The fecal metabolome provides a functional readout of microbial activity and can be used as an intermediate phenotype mediating host–microbiome interactions3. In this comprehensive description of the fecal metabolome, examining 1,116 metabolites from 786 individuals from a population-based twin study (TwinsUK), the fecal metabolome was found to be only modestly influenced by host genetics (heritability (H2) = 17.9%). One replicated locus at the NAT2 gene was associated with fecal metabolic traits. The fecal metabolome largely reflects gut microbial composition, explaining on average 67.7% (±18.8%) of its variance. It is strongly associated with visceral-fat mass, thereby illustrating potential mechanisms underlying the well-established microbial influence on abdominal obesity. Fecal metabolic profiling thus is a novel tool to explore links among microbiome composition, host phenotypes, and heritable complex traits