Identification and characterization of extensive intra-molecular associations between 3′-UTRs and their ORFs

During eukaryotic translation, mRNAs may form intra-molecular interactions between distant domains. The 5′-cap and the polyA tail were shown to interact through their associated proteins, and this can induce physical compaction of the mRNA in vitro. However, the stability of this intra-molecular association in translating mRNAs and whether additional contacts exist in vivo are largely unknown. To explore this, we applied a novel approach in which several endogenous polysomal mRNAs from Saccharomyces cerevisiae were cleaved near their stop codon and the resulting 3′-UTR fragments were tested either for co-sedimentation or co-immunoprecipitation (co-IP) with their ORFs. In all cases a significant fraction of the 3′-UTR fragments sedimented similarly to their ORF-containing fragments, yet the extent of co-sedimentation differed between mRNAs. Similar observations were obtained by a co-IP assay. Interestingly, various treatments that are expected to interfere with the cap to polyA interactions had no effect on the co-sedimentation pattern. Moreover, the 3′-UTR appeared to co-sediment with different regions from within the ORF. Taken together, these results indicate extensive physical associations between 3′-UTRs and their ORFs that vary between genes. This implies that polyribosomal mRNAs are in a compact configuration in vivo

Transcriptome and Proteome Exploration to Model Translation Efficiency and Protein Stability in Lactococcus lactis

This genome-scale study analysed the various parameters influencing protein levels in cells. To achieve this goal, the model bacterium Lactococcus lactis was grown at steady state in continuous cultures at different growth rates, and proteomic and transcriptomic data were thoroughly compared. Ratios of mRNA to protein were highly variable among proteins but also, for a given gene, between the different growth conditions. The modeling of cellular processes combined with a data fitting modeling approach allowed both translation efficiencies and degradation rates to be estimated for each protein in each growth condition. Estimated translational efficiencies and degradation rates strongly differed between proteins and were tested for their biological significance through statistical correlations with relevant parameters such as codon or amino acid bias. These efficiencies and degradation rates were not constant in all growth conditions and were inversely proportional to the growth rate, indicating a more efficient translation at low growth rate but an antagonistic higher rate of protein degradation. Estimated protein median half-lives ranged from 23 to 224 min, underlying the importance of protein degradation notably at low growth rates. The regulation of intracellular protein level was analysed through regulatory coefficient calculations, revealing a complex control depending on protein and growth conditions. The modeling approach enabled translational efficiencies and protein degradation rates to be estimated, two biological parameters extremely difficult to determine experimentally and generally lacking in bacteria. This method is generic and can now be extended to other environments and/or other micro-organisms

Transient Phenomena in Gene Expression after Induction of Transcription

When transcription of a gene is induced by a stimulus, the number of its mRNA molecules changes with time. Here we discuss how this time evolution depends on the shape of the mRNA lifetime distribution. Analysis of the statistical properties of this change reveals transient effects on polysomes, ribosomal profiles, and rate of protein synthesis. Our studies reveal that transient phenomena in gene expression strongly depend on the specific form of the mRNA lifetime distribution

The 3-Base Periodicity and Codon Usage of Coding Sequences Are Correlated with Gene Expression at the Level of Transcription Elongation

Background: Gene transcription is regulated by DNA transcriptional regulatory elements, promoters and enhancers that are located outside the coding regions. Here, we examine the characteristic 3-base periodicity of the coding sequences and analyse its correlation with the genome-wide transcriptional profile of yeast. Principal Findings: The analysis of coding sequences by a new class of indices proposed here identified two different sources of 3-base periodicity: the codon frequency and the codon sequence. In exponentially growing yeast cells, the codon-frequency component of periodicity accounts for 71.9 % of the variability of the cellular mRNA by a strong association with the density of elongating mRNA polymerase II complexes. The mRNA abundance explains most of the correlation between the codon-frequency component of periodicity and protein levels. Furthermore, pyrimidine-ending codons of the four-fold degenerate small amino acids alanine, glycine and valine are associated with genes with double the transcription rate of those associated with purine-ending codons. Conclusions: We demonstrate that the 3-base periodicity of coding sequences is higher than expected by the codon usage frequency (CUF) and that its components, associated with codon bias and amino acid composition, are correlated with gene expression, principally at the level of transcription elongation. This indicates a role of codon sequences in maximising the transcription efficiency in exponentially growing yeast cells. Moreover, the results contrast with the common Darwinia

Pichia pastoris regulates its gene-specific response to different carbon sources at the transcriptional, rather than the translational, level

Background: The methylotrophic, Crabtree-negative yeast Pichia pastoris is widely used as a heterologous protein production host. Strong inducible promoters derived from methanol utilization genes or constitutive glycolytic promoters are typically used to drive gene expression. Notably, genes involved in methanol utilization are not only repressed by the presence of glucose, but also by glycerol. This unusual regulatory behavior prompted us to study the regulation of carbon substrate utilization in different bioprocess conditions on a genome wide scale. Results: We performed microarray analysis on the total mRNA population as well as mRNA that had been fractionated according to ribosome occupancy. Translationally quiescent mRNAs were defined as being associated with single ribosomes (monosomes) and highly-translated mRNAs with multiple ribosomes (polysomes). We found that despite their lower growth rates, global translation was most active in methanol-grown P. pastoris cells, followed by excess glycerol- or glucose-grown cells. Transcript-specific translational responses were found to be minimal, while extensive transcriptional regulation was observed for cells grown on different carbon sources. Due to their respiratory metabolism, cells grown in excess glucose or glycerol had very similar expression profiles. Genes subject to glucose repression were mainly involved in the metabolism of alternative carbon sources including the control of glycerol uptake and metabolism. Peroxisomal and methanol utilization genes were confirmed to be subject to carbon substrate repression in excess glucose or glycerol, but were found to be strongly de-repressed in limiting glucose-conditions (as are often applied in fed batch cultivations) in addition to induction by methanol. Conclusions: P. pastoris cells grown in excess glycerol or glucose have similar transcript profiles in contrast to S. cerevisiae cells, in which the transcriptional response to these carbon sources is very different. The main response to different growth conditions in P. pastoris is transcriptional; translational regulation was not transcript-specific. The high proportion of mRNAs associated with polysomes in methanol-grown cells is a major finding of this study; it reveals that high productivity during methanol induction is directly linked to the growth condition and not only to promoter strength

Genome-Scale Analysis of Translation Elongation with a Ribosome Flow Model

We describe the first large scale analysis of gene translation that is based on a model that takes into account the physical and dynamical nature of this process. The Ribosomal Flow Model (RFM) predicts fundamental features of the translation process, including translation rates, protein abundance levels, ribosomal densities and the relation between all these variables, better than alternative (‘non-physical’) approaches. In addition, we show that the RFM can be used for accurate inference of various other quantities including genes' initiation rates and translation costs. These quantities could not be inferred by previous predictors. We find that increasing the number of available ribosomes (or equivalently the initiation rate) increases the genomic translation rate and the mean ribosome density only up to a certain point, beyond which both saturate. Strikingly, assuming that the translation system is tuned to work at the pre-saturation point maximizes the predictive power of the model with respect to experimental data. This result suggests that in all organisms that were analyzed (from bacteria to Human), the global initiation rate is optimized to attain the pre-saturation point. The fact that similar results were not observed for heterologous genes indicates that this feature is under selection. Remarkably, the gap between the performance of the RFM and alternative predictors is strikingly large in the case of heterologous genes, testifying to the model's promising biotechnological value in predicting the abundance of heterologous proteins before expressing them in the desired host

Asc1 Supports Cell-Wall Integrity Near Bud Sites by a Pkc1 Independent Mechanism

Background: The yeast ribosomal protein Asc1 is a WD-protein family member. Its mammalian ortholog, RACK1 was initially discovered as a receptor for activated protein C kinase (PKC) that functions to maintain the active conformation of PKC and to support its movement to target sites. In the budding yeast though, a connection between Asc1p and the PKC signaling pathway has never been reported. Methodology/Principal Findings: In the present study we found that asc1-deletion mutant (asc1D) presents some of the hallmarks of PKC signaling mutants. These include an increased sensitivity to staurosporine, a specific Pkc1p inhibitor, and susceptibility to cell-wall perturbing treatments such as hypotonic- and heat shock conditions and zymolase treatment. Microscopic analysis of asc1D cells revealed cell-wall invaginations near bud sites after exposure to hypotonic conditions, and the dynamic of cells ’ survival after this stress further supports the involvement of Asc1p in maintaining the cell-wall integrity during the mid-to late stages of bud formation. Genetic interactions between asc1 and pkc1 reveal synergistic sensitivities of a double-knock out mutant (asc1D/pkc1D) to cell-wall stress conditions, and high basal level of PKC signaling in asc1D. Furthermore, Asc1p has no effect on the cellular distribution or redistribution of Pkc1p at optimal or at cell-wall stress conditions. Conclusions/Significance: Taken together, our data support the idea that unlike its mammalian orthologs, Asc1p act

Global quantification of mammalian gene expression control

Gene expression is a multistep process that involves the transcription, translation and turnover of messenger RNAs and proteins. Although it is one of the most fundamental processes of life, the entire cascade has never been quantified on a genome-wide scale. Here we simultaneously measured absolute mRNA and protein abundance and turnover by parallel metabolic pulse labelling for more than 5,000 genes in mammalian cells. Whereas mRNA and protein levels correlated better than previously thought, corresponding half-lives showed no correlation. Using a quantitative model we have obtained the first genome-scale prediction of synthesis rates of mRNAs and proteins. We find that the cellular abundance of proteins is predominantly controlled at the level of translation. Genes with similar combinations of mRNA and protein stability shared functional properties, indicating that half-lives evolved under energetic and dynamic constraints. Quantitative information about all stages of gene expression provides a rich resource and helps to provide a greater understanding of the underlying design principles

Minimization of Biosynthetic Costs in Adaptive Gene Expression Responses of Yeast to Environmental Changes

Yeast successfully adapts to an environmental stress by altering physiology and fine-tuning metabolism. This fine-tuning is achieved through regulation of both gene expression and protein activity, and it is shaped by various physiological requirements. Such requirements impose a sustained evolutionary pressure that ultimately selects a specific gene expression profile, generating a suitable adaptive response to each environmental change. Although some of the requirements are stress specific, it is likely that others are common to various situations. We hypothesize that an evolutionary pressure for minimizing biosynthetic costs might have left signatures in the physicochemical properties of proteins whose gene expression is fine-tuned during adaptive responses. To test this hypothesis we analyze existing yeast transcriptomic data for such responses and investigate how several properties of proteins correlate to changes in gene expression. Our results reveal signatures that are consistent with a selective pressure for economy in protein synthesis during adaptive response of yeast to various types of stress. These signatures differentiate two groups of adaptive responses with respect to how cells manage expenditure in protein biosynthesis. In one group, significant trends towards downregulation of large proteins and upregulation of small ones are observed. In the other group we find no such trends. These results are consistent with resource limitation being important in the evolution of the first group of stress responses