A global characterization of the translational and transcriptional programs induced by methionine restriction through ribosome profiling and RNA-seq

Background: Among twenty amino acids, methionine has a special role as it is coded by the translation initiation codon and methionyl-tRNAi (Met-tRNAi) is required for the assembly of the translation initiation complex. Thus methionine may play a special role in global gene regulation. Methionine has also been known to play important roles in cell growth, development, cancer, and aging. In this work, we characterize the translational and transcriptional programs induced by methionine restriction (MetR) and investigate the potential mechanisms through which methionine regulates gene expression, using the budding yeast S. cerevisiae as the model organism. Results: Using ribosomal profiling and RNA-seq, we observed a broad spectrum of gene expression changes in response to MetR and identified hundreds of genes whose transcript level and/or translational efficiency changed significantly. These genes show clear functional themes, suggesting that cell slows down its growth and cell cycle progression and increases its stress resistance and maintenance in response to MetR. Interestingly, under MetR cell also decreases glycolysis and increases respiration, and increased respiration was linked to lifespan extension caused by caloric restriction. Analysis of genes whose translational efficiency changed significantly under MetR revealed different modes of translational regulation: 1) Ribosome loading patterns in the 5'UTR and coding regions of genes with increased translational efficiency suggested mechanisms both similar and different from that for the translational regulation of Gcn4 under general amino acid starvation condition; 2) Genes with decreased translational efficiency showed strong enrichment of lysine, glutamine, and glutamate codons, supporting the model that methionine can regulate translation by controlling tRNA thiolation. Conclusions: MetR induced a broad spectrum of gene expression changes at both the transcriptional and translational levels, with clear functional themes indicative of the physiological state of the cell under MetR. Different modes of translational regulation were induced by MetR, including the regulation of the ribosome loading at 5'UTR and regulation by tRNA thiolation. Since MetR extends the lifespan of many species, the list of genes we identified in this study can be good candidates for studying the mechanisms of lifespan extension.National Institutes of Health [AG043080]SCI(E)ARTICLE1

Mammalian microRNAs predominantly act to decrease target mRNA levels

MicroRNAs (miRNAs) are endogenous ~22-nucleotide RNAs that mediate important gene-regulatory events by pairing to the mRNAs of protein-coding genes to direct their repression. Repression of these regulatory targets leads to decreased translational efficiency and/or decreased mRNA levels, but the relative contributions of these two outcomes have been largely unknown, particularly for endogenous targets expressed at low-to-moderate levels. Here, we use ribosome profiling to measure the overall effects on protein production and compare these to simultaneously measured effects on mRNA levels. For both ectopic and endogenous miRNA regulatory interactions, lowered mRNA levels account for most (≥84%) of the decreased protein production. These results show that changes in mRNA levels closely reflect the impact of miRNAs on gene expression and indicate that destabilization of target mRNAs is the predominant reason for reduced protein output.National Institutes of Health (U.S.

Shape, Size, and Robustness: Feasible Regions in the Parameter Space of Biochemical Networks

The concept of robustness of regulatory networks has received much attention in the last decade. One measure of robustness has been associated with the volume of the feasible region, namely, the region in the parameter space in which the system is functional. In this paper, we show that, in addition to volume, the geometry of this region has important consequences for the robustness and the fragility of a network. We develop an approximation within which we could algebraically specify the feasible region. We analyze the segment polarity gene network to illustrate our approach. The study of random walks in the parameter space and how they exit the feasible region provide us with a rich perspective on the different modes of failure of this network model. In particular, we found that, between two alternative ways of activating Wingless, one is more robust than the other. Our method provides a more complete measure of robustness to parameter variation. As a general modeling strategy, our approach is an interesting alternative to Boolean representation of biochemical networks

Identification of a Topological Characteristic Responsible for the Biological Robustness of Regulatory Networks

Attribution of biological robustness to the specific structural properties of a regulatory network is an important yet unsolved problem in systems biology. It is widely believed that the topological characteristics of a biological control network largely determine its dynamic behavior, yet the actual mechanism is still poorly understood. Here, we define a novel structural feature of biological networks, termed ‘regulation entropy’, to quantitatively assess the influence of network topology on the robustness of the systems. Using the cell-cycle control networks of the budding yeast (Saccharomyces cerevisiae) and the fission yeast (Schizosaccharomyces pombe) as examples, we first demonstrate the correlation of this quantity with the dynamic stability of biological control networks, and then we establish a significant association between this quantity and the structural stability of the networks. And we further substantiate the generality of this approach with a broad spectrum of biological and random networks. We conclude that the regulation entropy is an effective order parameter in evaluating the robustness of biological control networks. Our work suggests a novel connection between the topological feature and the dynamic property of biological regulatory networks

riboviz: analysis and visualization of ribosome profiling datasets

Abstract Background Using high-throughput sequencing to monitor translation in vivo, ribosome profiling can provide critical insights into the dynamics and regulation of protein synthesis in a cell. Since its introduction in 2009, this technique has played a key role in driving biological discovery, and yet it requires a rigorous computational toolkit for widespread adoption. Description We have developed a database and a browser-based visualization tool, riboviz, that enables exploration and analysis of riboseq datasets. In implementation, riboviz consists of a comprehensive and flexible computational pipeline that allows the user to analyze private, unpublished datasets, along with a web application for comparison with published yeast datasets. Source code and detailed documentation are freely available from https://github.com/shahpr/RiboViz . The web-application is live at www.riboviz.org. Conclusions riboviz provides a comprehensive database and analysis and visualization tool to enable comparative analyses of ribosome-profiling datasets. This toolkit will enable both the community of systems biologists who study genome-wide ribosome profiling data and also research groups focused on individual genes to identify patterns of transcriptional and translational regulation across different organisms and conditions

Comprehensive comparative analysis of strand-specific RNA sequencing methods

Strand-specific, massively parallel cDNA sequencing (RNA-seq) is a powerful tool for transcript discovery, genome annotation and expression profiling. There are multiple published methods for strand-specific RNA-seq, but no consensus exists as to how to choose between them. Here we developed a comprehensive computational pipeline to compare library quality metrics from any RNA-seq method. Using the well-annotated Saccharomyces cerevisiae transcriptome as a benchmark, we compared seven library-construction protocols, including both published and our own methods. We found marked differences in strand specificity, library complexity, evenness and continuity of coverage, agreement with known annotations and accuracy for expression profiling. Weighing each method's performance and ease, we identified the dUTP second-strand marking and the Illumina RNA ligation methods as the leading protocols, with the former benefitting from the current availability of paired-end sequencing. Our analysis provides a comprehensive benchmark, and our computational pipeline is applicable for assessment of future protocols in other organisms.Howard Hughes Medical InstituteUnited States-Israel Binational Science Foundatio

Robustness under Functional Constraint: The Genetic Network for Temporal Expression in Drosophila Neurogenesis

Precise temporal coordination of gene expression is crucial for many developmental processes. One central question in developmental biology is how such coordinated expression patterns are robustly controlled. During embryonic development of the Drosophila central nervous system, neural stem cells called neuroblasts express a group of genes in a definite order, which leads to the diversity of cell types. We produced all possible regulatory networks of these genes and examined their expression dynamics numerically. From the analysis, we identified requisite regulations and predicted an unknown factor to reproduce known expression profiles caused by loss-of-function or overexpression of the genes in vivo, as well as in the wild type. Following this, we evaluated the stability of the actual Drosophila network for sequential expression. This network shows the highest robustness against parameter variations and gene expression fluctuations among the possible networks that reproduce the expression profiles. We propose a regulatory module composed of three types of regulations that is responsible for precise sequential expression. This study suggests that the Drosophila network for sequential expression has evolved to generate the robust temporal expression for neuronal specification

Multiple Aggregates and Aggresomes of C-Terminal Truncated Human αA-Crystallins in Mammalian Cells and Protection by αB-Crystallin

Cleavage of 11 (αA162), 5 (αA168) and 1 (αA172) residues from the C-terminus of αA-crystallin creates structurally and functionally different proteins. The formation of these post-translationally modified αA-crystallins is enhanced in diabetes. In the present study, the fate of the truncated αA-crystallins expressed in living mammalian cells in the presence and absence of native αA- or αB-crystallin has been studied by laser scanning confocal microscopy (LSM).YFP tagged αAwt, αA162, αA168 and αA172, were individually transfected or co-transfected with CFP tagged αAwt or αBwt, expressed in HeLa cells and studied by LSM. Difference in protein aggregation was not caused by different level of α-crystallin expression because Western blotting results showed nearly same level of expression of the various α-crystallins. The FRET-acceptor photo-bleaching protocol was followed to study in situ protein-protein interaction. αA172 interacted with αAwt and αBwt better than αA168 and αA162, interaction of αBwt being two-fold stronger than that of αAwt. Furthermore, aggresomes were detected in cells individually expressing αA162 and αA168 constructs and co-expression with αBwt significantly sequestered the aggresomes. There was no sequestration of aggresomes with αAwt co-expression with the truncated constructs, αA162 and αA168. Double immunocytochemistry technique was used for co-localization of γ-tubulin with αA-crystallin to demonstrate the perinuclear aggregates were aggresomes.αA172 showed the strongest interaction with both αAwt and αBwt. Native αB-crystallin provided protection to partially unfolded truncated αA-crystallins whereas native αA-crystallin did not. Aggresomes were detected in cells expressing αA162 and αA168 and αBwt co-expression with these constructs diminished the aggresome formation. Co-localization of γ-tubulin in perinuclear aggregates validates for aggresomes