Metabolic Gene Card for ORF: SKI2/YLR398C

This Metabolic Gene Card gives an overview of the metabolic profile upon deletion of the gene SKI2/YLR398C. It further lists other genes that have a similar metabolic impact, and hence, likely have a similar function. Further, we provide an enrichment analysis to dissect the obtained functional clusters, and compare the obtained information with physical and genetic interaction networks.


Ref. Functional Metabolomics Describes the Yeast Biosynthetic Regulome
Michael Mülleder, Enrica Calvani, Mohammad Tauqeer Alam, Richard Kangda Wang, Florian Eckerstorfer, Aleksej Zelezniak and Markus Ralser
Cell. 2016          Ralser laboratory

The gene card provides a summary about the response of the prototrophic Saccharomyces cerevisiae deletion strain SKI2/YLR398C under exponential growth in minimum synthetic medium - a condition were biosynthetic metabolism is expected to be highly active. The concentration of the free amino acid profile is a direct phenotypic readout for the gene deletion and allows us to estimate the genetic contribution of each of 4678 tested genes towards amino acid homeostasis.
Genes with similar function are expected to cause a similar metabolic signature. We assign metabolically similar strains to 280 distinct groups using hierarchical clustering and test them for overrepresentation of gene ontology (GO) functional terms, phenotypes and co-citation in literature. Thereby, we can suggest a gene function for not well characterised genes.

Free amino acid concentration profile

Amino acid levels are given in relation to the average of 4678 deletion strains. Robust estimates for mean, standard deviation, and covariance for each amino acid, and the Mahalanobis distance as metric for the overall change on the amino acid profile, were calculated using the Minimum Covariance Determinant (MCD) (R package: robustbase). Accordingly, we determined significance of change in the univariate case with a Z-test and in the multivariate statistic with a χ2-test statistic. All p values were adjusted for multiple testing with the FDR method from Benjamini & Hochberg (p.adjust function in R).

plot of chunk unnamed-chunk-4

Identification of 1519 altered amino acid profiles, 1459 by multivariate (χ2-test) and 1069 by univariate statistics (Z-test). Deletion strain SKI2/YLR398C is highlighted.

plot of chunk unnamed-chunk-5

The amino acids are variance scaled (z-scores) and ranked from low to high for easier interpretation. Significance of change is indicated in color and values given in table below.

to the top

Functional information derived and recapitulated from similarity of metabolic profiles

Similar metabolic profiles are indicative for a common genetic role in the regulation of amino acid metabolism. Furthermore, they are informative about functional relationship and therefore valuable for hypothesis generation and for functional characterisation of not well characterised open reading frames.

Groups of strains with similar amino acid profiles were identified by consensus clustering using 500 repetitions of hierarchical clustering based on pairwise Mahalanobis distances upon perturbation of the data. The obtained consensus matrix, which records the number of times two strains were assigned to the same cluster, was substracted from 1 and used as distance metric to identify 280 metabolically distinct and robust clusters (R package dynamicTreeCut).

Deletion mutants are assigned to groups with similar metabolic profiles

plot of chunk unnamed-chunk-7

Illustrated are the concentration values in a heatmap for deletion strains with a similar amino acid profile to SKI2/YLR398C. The dendrogram was calculated using Wards Criterion as agglomeration method. Concentration values are variance scaled (z-scores) and magnitude of change indicated in color.

Support for functional similarity by available genetic and physical interaction networks

To summarise the available data for the set of metabolically similar strains we construct a network based on genetic and physical interactions obtained from the BIOGRID database using the freely available esyn toolkit (2014 Bean et al). The online tool organises the information from BIOGRID and adds several more functionalities to analyse the network. It allows to retrieve information for each interaction, add common interactors, manually add or delete genes, perform network analysis tasks, or integrate data from Yeastmine.

Genes of similar metabolic signature: Enrichment analysis for Gene Ontology, KEGG, phenotypic similarity or co-citations

Enrichment analysis is a further tool to put new findings into the context of established knowledge. In order to find predictive enrichments, we remove the strain of interest from the cluster and the background set before the enrichment analysis. We use the gene ontology definition (GO Consortium), KEGG annotation and the phenotype and literature database from SGD. Enrichment analysis for GO annotations was performed in R with the GOstats package. In this package p values were not adjusted for multiple testing.

Enrichment analysis for phenotypes and literature annotations was done using hypergeometric testing. Obtained p values were corrected with the FDR method from Benjamini & Hochberg (p.adjust function in R).

Enrichment analysis based on metabolic signature for:

Gene ontology enrichment based on metabolic signature: biological process

Gene ontology enrichment based on metabolic signature: molecular function

Gene ontology enrichment based on metabolic signature: cellular component

to the top

Enrichment analysis based on metabolic signature: KEGG functional terms

Enrichment analysis based on metabolic signature: phenotypic similarity

Enrichment analysis based on metabolic signature: co-citations

to the top


Mülleder, M., Calvani, E., Alam, M.T., Wang, R.K., Eckerstorfer, F., Zelezniak, A., and Ralser, M. (2016). Functional Metabolomics Describes the Yeast Biosynthetic Regulome. Cell. dx.doi.org/10.1016/j.cell.2016.09.007

Information on gene ontology terms, KEGG terms, phenotype and co-citations was retrieved on 2016-08-10.