Measuring differential gene expression by short read sequencing: quantitative comparison to 2-channel gene expression microarrays

<p>Abstract</p> <p>Background</p> <p>High-throughput cDNA synthesis and sequencing of poly(A)-enriched RNA is rapidly emerging as a technology competing to replace microarrays as a quantitative platform for measuring gene expression.</p> <p>Results</p> <p>Consequently, we compared full length cDNA sequencing to 2-channel gene expression microarrays in the context of measuring differential gene expression. Because of its comparable cost to a gene expression microarray, our study focused on the data obtainable from a single lane of an Illumina 1 G sequencer. We compared sequencing data to a highly replicated microarray experiment profiling two divergent strains of <it>S. cerevisiae</it>.</p> <p>Conclusion</p> <p>Using a large number of quantitative PCR (qPCR) assays, more than previous studies, we found that neither technology is decisively better at measuring differential gene expression. Further, we report sequencing results from a diploid hybrid of two strains of <it>S. cerevisiae </it>that indicate full length cDNA sequencing can discover heterozygosity and measure quantitative allele-specific expression simultaneously.</p

Accurate proteome-wide protein quantification from high-resolution 15N mass spectra

In quantitative mass spectrometry-based proteomics, the metabolic incorporation of a single source of 15N-labeled nitrogen has many advantages over using stable isotope-labeled amino acids. However, the lack of a robust computational framework for analyzing the resulting spectra has impeded wide use of this approach. We have addressed this challenge by introducing a new computational methodology for analyzing 15N spectra in which quantification is integrated with identification. Application of this method to an Escherichia coli growth transition reveals significant improvement in quantification accuracy over previous methods

Riboneogenesis in Yeast

SummaryGlucose is catabolized in yeast via two fundamental routes, glycolysis and the oxidative pentose phosphate pathway, which produces NADPH and the essential nucleotide component ribose-5-phosphate. Here, we describe riboneogenesis, a thermodynamically driven pathway that converts glycolytic intermediates into ribose-5-phosphate without production of NADPH. Riboneogenesis begins with synthesis, by the combined action of transketolase and aldolase, of the seven-carbon bisphosphorylated sugar sedoheptulose-1,7-bisphosphate. In the pathway's committed step, sedoheptulose bisphosphate is hydrolyzed to sedoheptulose-7-phosphate by the enzyme sedoheptulose-1,7-bisphosphatase (SHB17), whose activity we identified based on metabolomic analysis of the corresponding knockout strain. The crystal structure of Shb17 in complex with sedoheptulose-1,7-bisphosphate reveals that the substrate binds in the closed furan form in the active site. Sedoheptulose-7-phosphate is ultimately converted by known enzymes of the nonoxidative pentose phosphate pathway to ribose-5-phosphate. Flux through SHB17 increases when ribose demand is high relative to demand for NADPH, including during ribosome biogenesis in metabolically synchronized yeast cells

Directing Experimental Biology: A Case Study in Mitochondrial Biogenesis

Computational approaches have promised to organize collections of functional genomics data into testable predictions of gene and protein involvement in biological processes and pathways. However, few such predictions have been experimentally validated on a large scale, leaving many bioinformatic methods unproven and underutilized in the biology community. Further, it remains unclear what biological concerns should be taken into account when using computational methods to drive real-world experimental efforts. To investigate these concerns and to establish the utility of computational predictions of gene function, we experimentally tested hundreds of predictions generated from an ensemble of three complementary methods for the process of mitochondrial organization and biogenesis in Saccharomyces cerevisiae. The biological data with respect to the mitochondria are presented in a companion manuscript published in PLoS Genetics (doi:10.1371/journal.pgen.1000407). Here we analyze and explore the results of this study that are broadly applicable for computationalists applying gene function prediction techniques, including a new experimental comparison with 48 genes representing the genomic background. Our study leads to several conclusions that are important to consider when driving laboratory investigations using computational prediction approaches. While most genes in yeast are already known to participate in at least one biological process, we confirm that genes with known functions can still be strong candidates for annotation of additional gene functions. We find that different analysis techniques and different underlying data can both greatly affect the types of functional predictions produced by computational methods. This diversity allows an ensemble of techniques to substantially broaden the biological scope and breadth of predictions. We also find that performing prediction and validation steps iteratively allows us to more completely characterize a biological area of interest. While this study focused on a specific functional area in yeast, many of these observations may be useful in the contexts of other processes and organisms

Predicting Cellular Growth from Gene Expression Signatures

Maintaining balanced growth in a changing environment is a fundamental systems-level challenge for cellular physiology, particularly in microorganisms. While the complete set of regulatory and functional pathways supporting growth and cellular proliferation are not yet known, portions of them are well understood. In particular, cellular proliferation is governed by mechanisms that are highly conserved from unicellular to multicellular organisms, and the disruption of these processes in metazoans is a major factor in the development of cancer. In this paper, we develop statistical methodology to identify quantitative aspects of the regulatory mechanisms underlying cellular proliferation in Saccharomyces cerevisiae. We find that the expression levels of a small set of genes can be exploited to predict the instantaneous growth rate of any cellular culture with high accuracy. The predictions obtained in this fashion are robust to changing biological conditions, experimental methods, and technological platforms. The proposed model is also effective in predicting growth rates for the related yeast Saccharomyces bayanus and the highly diverged yeast Schizosaccharomyces pombe, suggesting that the underlying regulatory signature is conserved across a wide range of unicellular evolution. We investigate the biological significance of the gene expression signature that the predictions are based upon from multiple perspectives: by perturbing the regulatory network through the Ras/PKA pathway, observing strong upregulation of growth rate even in the absence of appropriate nutrients, and discovering putative transcription factor binding sites, observing enrichment in growth-correlated genes. More broadly, the proposed methodology enables biological insights about growth at an instantaneous time scale, inaccessible by direct experimental methods. Data and tools enabling others to apply our methods are available at http://function.princeton.edu/growthrate

Uniform nomenclature for the mitochondrial contact site and cristae organizing system

The mitochondrial inner membrane contains a large protein complex that functions in inner membrane organization and formation of membrane contact sites. The complex was variably named the mitochondrial contact site complex, mitochondrial inner membrane organizing system, mitochondrial organizing structure, or Mitofilin/Fcj1 complex. To facilitate future studies, we propose to unify the nomenclature and term the complex "mitochondrial contact site and cristae organizing system" and its subunits Mic10 to Mic60