Global proteome turnover analyses of the Yeasts S. cerevisiae and S. pombe

Quantitating translational control: mRNA abundance-dependent and independent contributions and the mRNA sequences that specify them

Translation rate per mRNA molecule correlates positively with mRNA abundance. As a result, protein levels do not scale linearly with mRNA levels, but instead scale with the abundance of mRNA raised to the power of an ‘amplification exponent’. Here we show that to quantitate translational control, the translation rate must be decomposed into two components. One, TR_mD, depends on the mRNA level and defines the amplification exponent. The other, TR_mIND, is independent of mRNA amount and impacts the correlation coefficient between protein and mRNA levels. We show that in Saccharomyces cerevisiae TR_mD represents ∼20% of the variance in translation and directs an amplification exponent of 1.20 with a 95% confidence interval [1.14, 1.26]. TR_mIND constitutes the remaining ∼80% of the variance in translation and explains ∼5% of the variance in protein expression. We also find that TR_mD and TR_mIND are preferentially determined by different mRNA sequence features: TR_mIND by the length of the open reading frame and TR_mD both by a ∼60 nucleotide element that spans the initiating AUG and by codon and amino acid frequency. Our work provides more appropriate estimates of translational control and implies that TR_mIND is under different evolutionary selective pressures than TR_mD.

Curated protein information in the Saccharomyces genome database

Due to recent advancements in the production of experimental proteomic data, the Saccharomyces genome database (SGD; www.yeastgenome.org) has been expanding our protein curation activities to make new data types available to our users. Because of broad interest in post-translational modifications (PTM) and their importance to protein function and regulation, we have recently started incorporating expertly curated PTM information on individual protein pages. Here we also present the inclusion of new abundance and protein half-life data obtained from high-throughput proteome studies. These new data types have been included with the aim to facilitate cellular biology research.

Database URL: www.yeastgenome.org

CRISPR/Cas9-based genome editing for simultaneous interference with gene expression and protein stability

Interference with genes is the foundation of reverse genetics and is key to manipulation of living cells for biomedical and biotechnological applications. However, classical genetic knockout and transcriptional knockdown technologies have different drawbacks and offer no control over existing protein levels. Here, we describe an efficient genome editing approach that affects specific protein abundances by changing the rates of both RNA synthesis and protein degradation, based on the two cross-kingdom control mechanisms CRISPRi and the N-end rule for protein stability. In addition, our approach demonstrates that CRISPRi efficiency is dependent on endogenous gene expression levels. The method has broad applications in e.g. study of essential genes and antibiotics discovery.

A reference-based protein degradation assay without global translation inhibitors

Although it is widely appreciated that the use of global translation inhibitors, such as cycloheximide, in protein degradation assays may result in artefacts, these inhibitors continue to be employed, owing to the absence of robust alternatives. We describe here the promoter reference technique (PRT), an assay for protein degradation with two advantageous features: a reference protein and a gene-specific inhibition of translation. In PRT assays, one measures, during a chase, the ratio of a test protein to a long-lived reference protein, a dihydrofolate reductase (DHFR). The test protein and DHFR are coexpressed, in the yeast Saccharomyces cerevisiae, on a low-copy plasmid from two identical P _TDH3 promoters containing additional, previously developed DNA elements. Once transcribed, these elements form 5'-RNA aptamers that bind to the added tetracycline, which represses translation of aptamer-containing mRNAs. The selectivity of repression avoids a global inhibition of translation. This selectivity is particularly important if a component of a relevant proteolytic pathway (e.g. a specific ubiquitin ligase) is itself short-lived. We applied PRT to the Pro/N-end rule pathway, whose substrates include the short-lived Mdh2 malate dehydrogenase. Mdh2 is targeted for degradation by the Gid4 subunit of the GID ubiquitin ligase. Gid4 is also a metabolically unstable protein. Through analyses of short-lived Mdh2 as a target of short-lived Gid4, we illustrate the advantages of PRT over degradation assays that lack a reference and/or involve cycloheximide. In sum, PRT avoids the use of global translation inhibitors during a chase and also provides a "built-in" reference protein.

Perspective on architecture and assembly of membrane contact sites

Membrane contact sites (MCS) are platforms of physical contact between different organelles. They are formed through interactions involving lipids and proteins, and function in processes such as calcium and lipid exchange, metabolism and organelle biogenesis. In this article, we discuss emerging questions regarding the architecture, organisation and assembly of MCS, such as: What is the contribution of different components to the interaction between organelles? How is the specific composition of different types of membrane contacts sites established and maintained? How are proteins and lipids spatially organised at MCS and how does that influence their function? How dynamic are MCS on the molecular and ultrastructural level? We highlight current state of research and point out experimental approaches that promise to contribute to a spatiomechanistic understanding of MCS functions.

Cell-wide analysis of protein thermal unfolding reveals determinants of thermostability

Temperature-induced cell death is thought to be due to protein denaturation, but the determinants of thermal sensitivity of proteomes remain largely uncharacterized. We developed a structural proteomic strategy to measure protein thermostability on a proteome-wide scale and with domain-level resolution. We applied it to Escherichia coli, Saccharomyces cerevisiae, Thermus thermophilus, and human cells, yielding thermostability data for more than 8000 proteins. Our results (i) indicate that temperature-induced cellular collapse is due to the loss of a subset of proteins with key functions, (ii) shed light on the evolutionary conservation of protein and domain stability, and (iii) suggest that natively disordered proteins in a cell are less prevalent than predicted and (iv) that highly expressed proteins are stable because they are designed to tolerate translational errors that would lead to the accumulation of toxic misfolded species.

Determinants and Regulation of Protein Turnover in Yeast

Protein turnover maintains the recycling needs of the proteome, and its malfunction has been linked to aging and age-related diseases. However, not all proteins turnover equally, and the factors that contribute to accelerate or slow down turnover are mostly unknown. We measured turnover rates for 3,160 proteins in exponentially growing yeast and analyzed their dependence on physical, functional, and genetic properties. We found that functional characteristics, including protein localization, complex membership, and connectivity, have greater effect on turnover than sequence elements. We also found that protein turnover and mRNA turnover are correlated. Analysis under nutrient perturbation and osmotic stress revealed that protein turnover highly depends on cellular state and is faster when proteins are being actively used. Finally, stress-induced changes in protein and transcript abundance correlated with changes in protein turnover. This study provides a resource of protein turnover rates and principles to understand the recycling needs of the proteome under basal conditions and perturbation.

Tracking the evolution of 3D gene organization demonstrates its connection to phenotypic divergence

It has recently been shown that the organization of genes in eukaryotic genomes, and specifically in 3D, is strongly related to gene expression and function and partially conserved between organisms. However, previous studies of 3D genomic organization analyzed each organism independently from others. Here, we propose an approach for unified inter-organismal analysis of gene organization based on a network representation of Hi-C data. We define and detect four classes of spatially co-evolving orthologous modules (SCOMs), i.e. gene families that co-evolve in their 3D organization, based on patterns of divergence and conservation of distances. We demonstrate our methodology on Hi-C data from Saccharomyces cerevisiae and Schizosaccharomyces pombe, and identify, among others, modules relating to RNA splicing machinery and chromatin silencing by small RNA which are central to S. pombe's lifestyle. Our results emphasize the importance of 3D genomic organization in eukaryotes and suggest that the evolutionary mechanisms that shape gene organization affect the organism fitness and phenotypes. The proposed algorithms can be utilized in future studies of genome evolution and comparative analysis of spatial genomic organization in different tissues, conditions and single cells.

Regulation of Rab5 isoforms by transcriptional and post-transcriptional mechanisms in yeast

Rab5 GTPases are master regulators of early endosome biogenesis and transport. The genome of Saccharomyces cerevisiae encodes three Rab5 proteins: Vps21, the major isoform, Ypt52 and Ypt53. Here, we show that Vps21 is the most abundant Rab5 protein and Ypt53 is the least abundant. In stressed cells, Ypt53 levels increase but never exceed that of Vps21. Its induction requires the transcription factors Crz1 and Gis1. In growing cells, the expression of Ypt53 is suppressed by post-transcriptional mechanisms mediated by the untranslated regions of the YPT53 mRNA. Based on genetic experiments, these sequences appear to stimulate deadenylation, Pat1-activated decapping and Xrn1-mediated mRNA degradation. Once this regulation is bypassed, Ypt53 protein levels surpass Vps21, and Ypt53 is sufficient to maintain endosomal function and cell growth.

ACBD3 functions as a scaffold to organize the Golgi stacking proteins and a Rab33b-GAP

Golgin45 plays important roles in Golgi stack assembly and is known to bind both the Golgi stacking protein GRASP55 and Rab2 in the medial-Golgi cisternae. In this study, we sought to further characterize the cisternal adhesion complex using a proteomics approach. We report here that Acyl-CoA binding domain containing 3 (ACBD3) is likely to be a novel binding partner of Golgin45. ACBD3 interacts with Golgin45 via its GOLD domain, while its co-expression significantly increases Golgin45 targeting to the Golgi. Furthermore, ACBD3 recruits TBC1D22, a Rab33b GTPase activating protein (GAP), to a large multi-protein complex containing Golgin45 and GRASP55. These results suggest that ACBD3 may provide a scaffolding to organize the Golgi stacking proteins and a Rab33b-GAP at the medial-Golgi.

A Simple Light Isotope Metabolic Labeling (SLIM-labeling) Strategy: A Powerful Tool to Address the Dynamics of Proteome Variations In Vivo

Many quantitative proteomics strategies rely on in vivo metabolic incorporation of amino acids with modified stable isotope profiles into proteins. These methods give rise to multiple ions for each peptide, with possible distortion of the isotopolog distribution, making the overall analytical process complex. We validated an alternative strategy, simple light isotope metabolic labeling (SLIM-labeling), which alleviates many of these problems. SLIM-labeling is based on the in vivo reduction of the isotopic composition of proteins using metabolic precursors with a unique light isotope composition to label all amino acids. This brings a new dimension to in-depth, high resolution MS-based quantitative proteomics. Here, we describe a ¹²C-based SLIM-labeling strategy using U-[¹²C]-glucose as the metabolic precursor of all amino acids in the pathogenic yeast Candida albicans Monoisotopic ion intensity increased exponentially following ¹²C enrichment, substantially improving peptide identification scores and protein sequence coverage in bottom-up analyses. Multiplexing samples of ¹²C composition varying from natural abundance (98.93%) to 100% makes it possible to address relative quantification issues, keeping all the critical information for each peptide within a single isotopolog cluster. We applied this method to measure, for the first time, protein turnover at the proteome scale in Candida albicans and its modulation by inhibitors of the proteasome and vacuolar protein degradation systems.

High-throughput Screening for Protein-based Inheritance in S. cerevisiae

The encoding of biological information that is accessible to future generations is generally achieved via changes to the DNA sequence. Long-lived inheritance encoded in protein conformation (rather than sequence) has long been viewed as paradigm-shifting but rare. The best characterized examples of such epigenetic elements are prions, which possess a self-assembling behavior that can drive the heritable manifestation of new phenotypes. Many archetypal prions display a striking N/Q-rich sequence bias and assemble into an amyloid fold. These unusual features have informed most screening efforts to identify new prion proteins. However, at least three known prions (including the founding prion, PrP^Sc) do not harbor these biochemical characteristics. We therefore developed an alternative method to probe the scope of protein-based inheritance based on a property of mass action: the transient overexpression of prion proteins increases the frequency at which they acquire a self-templating conformation. This paper describes a method for analyzing the capacity of the yeast ORFeome to elicit protein-based inheritance. Using this strategy, we previously found that >1% of yeast proteins could fuel the emergence of biological traits that were long-lived, stable, and arose more frequently than genetic mutation. This approach can be employed in high throughput across entire ORFeomes or as a targeted screening paradigm for specific genetic networks or environmental stimuli. Just as forward genetic screens define numerous developmental and signaling pathways, these techniques provide a methodology to investigate the influence of protein-based inheritance in biological processes.

Quantitative characterization of the auxin-inducible degron: a guide for dynamic protein depletion in single yeast cells

Perturbations are essential for the interrogation of biological systems. The auxin-inducible degron harbors great potential for dynamic protein depletion in yeast. Here, we thoroughly and quantitatively characterize the auxin-inducible degron in single yeast cells. We show that an auxin concentration of 0.25 mM is necessary for fast and uniform protein depletion between single cells, and that in mother cells proteins are depleted faster than their daughters. Although, protein recovery starts immediately after removal of auxin, it takes multiple generations before equilibrium is reached between protein synthesis and dilution, which is when the original protein levels are restored. Further, we found that blue light, used for GFP excitation, together with auxin results in growth defects, caused by the photo-destruction of auxin to its toxic derivatives, which can be avoided if indole-free auxin substitutes are used. Our work provides guidelines for the successful combination of microscopy, microfluidics and the auxin-inducible degron, offering the yeast community an unprecedented tool for dynamic perturbations on the single cell level.

Multiplexed gene control reveals rapid mRNA turnover

The rates of mRNA synthesis and decay determine the mRNA expression level. The two processes are under coordinated control, which makes the measurements of these rates challenging, as evidenced by the low correlation among the methods of measurement of RNA half-lives. We developed a minimally invasive method, multiplexed gene control, to shut off expression of genes with controllable synthetic promoters. The method was validated by measuring the ratios of the nascent to mature mRNA molecules and by measuring the half-life with endogenous promoters that can be controlled naturally or through inserting short sequences that impart repressibility. The measured mRNA half-lives correlated highly with those obtained with the metabolic pulse-labeling method in yeast. However, mRNA degradation was considerably faster in comparison to previous estimates, with a median half-life of around 2 min. The half-life permits the estimation of promoter-dependent and promoter-independent transcription rates. The dynamical range of the promoter-independent transcription rates was larger than that of the mRNA half-lives. The rapid mRNA turnover and the broad adjustability of promoter-independent transcription rates are expected to have a major impact on stochastic gene expression and gene network behavior.

Independent regulation of gene expression level and noise by histone modifications

The inherent stochasticity generates substantial gene expression variation among isogenic cells under identical conditions, which is frequently referred to as gene expression noise or cell-to-cell expression variability. Similar to (average) expression level, expression noise is also subject to natural selection. Yet it has been observed that noise is negatively correlated with expression level, which manifests as a potential constraint for simultaneous optimization of both. Here, we studied expression noise in human embryonic cells with computational analysis on single-cell RNA-seq data and in yeast with flow cytometry experiments. We showed that this coupling is overcome, to a certain degree, by a histone modification strategy in multiple embryonic developmental stages in human, as well as in yeast. Importantly, this epigenetic strategy could fit into a burst-like gene expression model: promoter-localized histone modifications (such as H3K4 methylation) are associated with both burst size and burst frequency, which together influence expression level, while gene-body-localized ones (such as H3K79 methylation) are more associated with burst frequency, which influences both expression level and noise. We further knocked out the only “writer” of H3K79 methylation in yeast, and observed that expression noise is indeed increased. Consistently, dosage sensitive genes, such as genes in the Wnt signaling pathway, tend to be marked with gene-body-localized histone modifications, while stress responding genes, such as genes regulating autophagy, tend to be marked with promoter-localized ones. Our findings elucidate that the “division of labor” among histone modifications facilitates the independent regulation of expression level and noise, extend the “histone code” hypothesis to include expression noise, and shed light on the optimization of transcriptome in evolution.

Gene expression noise, or cell-to-cell expression variability, has been a topic of intense interest for more than a decade. The prevailing model of “burst-like transcription” mediated by the promoter transitions between on and off states explains the formation of noise in eukaryotes. Albeit widely accepted, the cis- elements that determine the burst frequency and burst size remain largely unknown. Here we systematically examined the relationship between transcriptional burst frequency/size and all major histone modifications in various cell types, including human embryonic cells, mouse embryonic stem cells, and yeast, and found that histone markers can be divided into two groups based on their associations with burst frequency/ size. Coincidently, promoter-localized histone markers are associated with both burst size and burst frequency whereas gene-body-localized ones are more associated with burst frequency. We further knocked out a gene that is responsible for “writing” a gene-body histone mark in yeast, and found that burst frequency is indeed reduced. Our findings reveal a new mechanism of transcriptional burst regulation and shed light on the simultaneous optimization of gene expression level and noise in evolution.

ER stress causes widespread protein aggregation and prion formation

ER stress results in widespread aggregation of proteins that are not localized to the ER or are part of the secretory system. Hamdan et al. demonstrate that amorphous and amyloidogenic protein aggregation is an indirect consequence of perturbing ER homeostasis.

Disturbances in endoplasmic reticulum (ER) homeostasis create a condition termed ER stress. This activates the unfolded protein response (UPR), which alters the expression of many genes involved in ER quality control. We show here that ER stress causes the aggregation of proteins, most of which are not ER or secretory pathway proteins. Proteomic analysis of the aggregated proteins revealed enrichment for intrinsically aggregation-prone proteins rather than proteins which are affected in a stress-specific manner. Aggregation does not arise because of overwhelming proteasome-mediated degradation but because of a general disruption of cellular protein homeostasis. We further show that overexpression of certain chaperones abrogates protein aggregation and protects a UPR mutant against ER stress conditions. The onset of ER stress is known to correlate with various disease processes, and our data indicate that widespread amorphous and amyloid protein aggregation is an unanticipated outcome of such stress.

Longevity of major coenzymes allows minimal de novo synthesis in microorganisms

Coenzymes are vital for cellular metabolism and act on the full spectrum of enzymatic reactions. Intrinsic chemical reactivity, enzyme promiscuity and high flux through their catalytic cycles make coenzymes prone to damage. To counteract such compromising factors and ensure stable levels of functional coenzymes, cells use a complex interplay between de novo synthesis, salvage, repair and degradation. However, the relative contribution of these factors is currently unknown, as is the overall stability of coenzymes in the cell. Here, we use dynamic ¹³C-labelling experiments to determine the half-life of major coenzymes of Escherichia coli. We find that coenzymes such as pyridoxal 5-phosphate, flavins, nicotinamide adenine dinucleotide (phosphate) and coenzyme A are remarkably stable in vivo and allow biosynthesis close to the minimal necessary rate. In consequence, they are essentially produced to compensate for dilution by growth and passed on over generations of cells. Exceptions are antioxidants, which are short-lived, suggesting an inherent requirement for increased renewal. Although the growth-driven turnover of stable coenzymes is apparently subject to highly efficient end-product homeostasis, we exemplify that coenzyme pools are propagated in excess in relation to actual growth requirements. Additional testing of Bacillus subtilis and Saccharomyces cerevisiae suggests that coenzyme longevity is a conserved feature in biology.

The link between adjacent codon pairs and mRNA stability

Background

Evidence in diverse organisms suggests that codon optimality is a major determinant of mRNA translation and degradation. Codon optimality is thought to act by modulating the efficiency of ribosome elongation. In Saccharomyces cerevisiae, a recent study has identified 17 adjacent codon pairs that mediate strong inhibition of translation elongation. However, relationships between the inhibitory codon pairs and other aspects of gene expression are unknown.

Results

To gain insights into how the inhibitory codon pairs may affect aspects of gene expression, we utilized existing datasets to conduct genome-scale analyses in S. cerevisiae. Our analysis revealed the following points. First, the inhibitory codon pairs are significantly associated with faster mRNA decay. The association is not solely due to the content of nucleotides, individual codons, or dipeptides encoded by the inhibitory codon pairs. Second, the inhibitory codon pairs cannot fully explain the previously known relationship of codon optimality with mRNA stability, suggesting that optimality of individual codons and properties of adjacent codon pairs both contribute to gene regulation. Finally, although the inhibitory codon pairs are associated with slower mRNA synthesis and protein instability, the associations can be attributed to usage bias in individual codons.

Conclusions

This study suggests an association of inhibitory codon pairs with mRNA stability and thus another layer of complexity in the codon-mediated gene regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3749-8) contains supplementary material, which is available to authorized users.

The Pivotal Role of Protein Phosphorylation in the Control of Yeast Central Metabolism

Protein phosphorylation is the most frequent eukaryotic post-translational modification and can act as either a molecular switch or rheostat for protein functions. The deliberate manipulation of protein phosphorylation has great potential for regulating specific protein functions with surgical precision, rather than the gross effects gained by the over/underexpression or complete deletion of a protein-encoding gene. In order to assess the impact of phosphorylation on central metabolism, and thus its potential for biotechnological and medical exploitation, a compendium of highly confident protein phosphorylation sites (p-sites) for the model organism Saccharomyces cerevisiae has been analyzed together with two more datasets from the fungal pathogen Candida albicans. Our analysis highlights the global properties of the regulation of yeast central metabolism by protein phosphorylation, where almost half of the enzymes involved are subject to this sort of post-translational modification. These phosphorylated enzymes, compared to the nonphosphorylated ones, are more abundant, regulate more reactions, have more protein–protein interactions, and a higher fraction of them are ubiquitinated. The p-sites of metabolic enzymes are also more conserved than the background p-sites, and hundreds of them have the potential for regulating metabolite production. All this integrated information has allowed us to prioritize thousands of p-sites in terms of their potential phenotypic impact. This multi-source compendium should enable the design of future high-throughput (HTP) mutation studies to identify key molecular switches/rheostats for the manipulation of not only the metabolism of yeast, but also that of many other biotechnologically and medically important fungi and eukaryotes.

Protein Degradation Rate in Arabidopsis thaliana Leaf Growth and Development

We applied ¹⁵N labeling approaches to leaves of the Arabidopsis thaliana rosette to characterize their protein degradation rate and understand its determinants. The progressive labeling of new peptides with ¹⁵N and measuring the decrease in the abundance of >60,000 existing peptides over time allowed us to define the degradation rate of 1228 proteins in vivo. We show that Arabidopsis protein half-lives vary from several hours to several months based on the exponential constant of the decay rate for each protein. This rate was calculated from the relative isotope abundance of each peptide and the fold change in protein abundance during growth. Protein complex membership and specific protein domains were found to be strong predictors of degradation rate, while N-end amino acid, hydrophobicity, or aggregation propensity of proteins were not. We discovered rapidly degrading subunits in a variety of protein complexes in plastids and identified the set of plant proteins whose degradation rate changed in different leaves of the rosette and correlated with leaf growth rate. From this information, we have calculated the protein turnover energy costs in different leaves and their key determinants within the proteome.

Protein turnover measurement using selected reaction monitoring-mass spectrometry (SRM-MS)

Protein turnover represents an important mechanism in the functioning of cells, with deregulated synthesis and degradation of proteins implicated in many diseased states. Therefore, proteomics strategies to measure turnover rates with high confidence are of vital importance to understanding many biological processes. In this study, the more widely used approach of non-targeted precursor ion signal intensity (MS1) quantification is compared with selected reaction monitoring (SRM), a data acquisition strategy that records data for specific peptides, to determine if improved quantitative data would be obtained using a targeted quantification approach. Using mouse liver as a model system, turnover measurement of four tricarboxylic acid cycle proteins was performed using both MS1 and SRM quantification strategies. SRM outperformed MS1 in terms of sensitivity and selectivity of measurement, allowing more confident determination of protein turnover rates. SRM data are acquired using cheaper and more widely available tandem quadrupole mass spectrometers, making the approach accessible to a larger number of researchers than MS1 quantification, which is best performed on high mass resolution instruments. SRM acquisition is ideally suited to focused studies where the turnover of tens of proteins is measured, making it applicable in determining the dynamics of proteins complexes and complete metabolic pathways.This article is part of the themed issue 'Quantitative mass spectrometry'.

Chm7 and Heh1 collaborate to link nuclear pore complex quality control with nuclear envelope sealing

The integrity of the nuclear envelope barrier relies on membrane remodeling by the ESCRTs, which seal nuclear envelope holes and contribute to the quality control of nuclear pore complexes (NPCs); whether these processes are mechanistically related remains poorly defined. Here, we show that the ESCRT-II/III chimera, Chm7, is recruited to a nuclear envelope subdomain that expands upon inhibition of NPC assembly and is required for the formation of the storage of improperly assembled NPCs (SINC) compartment. Recruitment to sites of NPC assembly is mediated by its ESCRT-II domain and the LAP2-emerin-MAN1 (LEM) family of integral inner nuclear membrane proteins, Heh1 and Heh2. We establish direct binding between Heh2 and the "open" forms of both Chm7 and the ESCRT-III, Snf7, and between Chm7 and Snf7. Interestingly, Chm7 is required for the viability of yeast strains where double membrane seals have been observed over defective NPCs; deletion of CHM7 in these strains leads to a loss of nuclear compartmentalization suggesting that the sealing of defective NPCs and nuclear envelope ruptures could proceed through similar mechanisms.

Analysis of the Yeast Peptidome and Comparison with the Human Peptidome

Peptides function as signaling molecules in species as diverse as humans and yeast. Mass spectrometry-based peptidomics techniques provide a relatively unbiased method to assess the peptidome of biological samples. In the present study, we used a quantitative peptidomic technique to characterize the peptidome of the yeast Saccharomyces cerevisiae and compare it to the peptidomes of mammalian cell lines and tissues. Altogether, 297 yeast peptides derived from 75 proteins were identified. The yeast peptides are similar to those of the human peptidome in average size and amino acid composition. Inhibition of proteasome activity with either bortezomib or epoxomicin led to decreased levels of some yeast peptides, suggesting that these peptides are generated by the proteasome. Approximately 30% of the yeast peptides correspond to the N- or C-terminus of the protein; the human peptidome is also highly represented in N- or C-terminal protein fragments. Most yeast and humans peptides are derived from a subset of abundant proteins, many with functions involving cellular metabolism or protein synthesis and folding. Of the 75 yeast proteins that give rise to peptides, 24 have orthologs that give rise to human and/or mouse peptides and for some, the same region of the proteins are found in the human, mouse, and yeast peptidomes. Taken together, these results support the hypothesis that intracellular peptides may have specific and conserved biological functions.

Organelle acidification negatively regulates vacuole membrane fusion in vivo

The V-ATPase is a proton pump consisting of a membrane-integral V0 sector and a peripheral V1 sector, which carries the ATPase activity. In vitro studies of yeast vacuole fusion and evidence from worms, flies, zebrafish and mice suggested that V0 interacts with the SNARE machinery for membrane fusion, that it promotes the induction of hemifusion and that this activity requires physical presence of V0 rather than its proton pump activity. A recent in vivo study in yeast has challenged these interpretations, concluding that fusion required solely lumenal acidification but not the V0 sector itself. Here, we identify the reasons for this discrepancy and reconcile it. We find that acute pharmacological or physiological inhibition of V-ATPase pump activity de-acidifies the vacuole lumen in living yeast cells within minutes. Time-lapse microscopy revealed that de-acidification induces vacuole fusion rather than inhibiting it. Cells expressing mutated V0 subunits that maintain vacuolar acidity were blocked in this fusion. Thus, proton pump activity of the V-ATPase negatively regulates vacuole fusion in vivo. Vacuole fusion in vivo does, however, require physical presence of a fusion-competent V0 sector.

Maintenance-energy requirements and robustness of Saccharomyces cerevisiae at aerobic near-zero specific growth rates

Background

Saccharomyces cerevisiae is an established microbial platform for production of native and non-native compounds. When product pathways compete with growth for precursors and energy, uncoupling of growth and product formation could increase product yields and decrease formation of biomass as a by-product. Studying non-growing, metabolically active yeast cultures is a first step towards developing S. cerevisiae as a robust, non-growing cell factory. Microbial physiology at near-zero growth rates can be studied in retentostats, which are continuous-cultivation systems with full biomass retention. Hitherto, retentostat studies on S. cerevisiae have focused on anaerobic conditions, which bear limited relevance for aerobic industrial processes. The present study uses aerobic, glucose-limited retentostats to explore the physiology of non-dividing, respiring S. cerevisiae cultures, with a focus on industrially relevant features.

Results

Retentostat feeding regimes for smooth transition from exponential growth in glucose-limited chemostat cultures to near-zero growth rates were obtained by model-aided experimental design. During 20 days of retentostats cultivation, the specific growth rate gradually decreased from 0.025 h⁻¹ to below 0.001 h⁻¹, while culture viability remained above 80 %. The maintenance requirement for ATP (m_ATP) was estimated at 0.63 ± 0.04 mmol ATP (g biomass)⁻¹ h⁻¹, which is ca. 35 % lower than previously estimated for anaerobic retentostats. Concomitant with decreasing growth rate in aerobic retentostats, transcriptional down-regulation of genes involved in biosynthesis and up-regulation of stress-responsive genes resembled transcriptional regulation patterns observed for anaerobic retentostats. The heat-shock tolerance in aerobic retentostats far exceeded previously reported levels in stationary-phase batch cultures. While in situ metabolic fluxes in retentostats were intentionally low due to extreme caloric restriction, off-line measurements revealed that cultures retained a high metabolic capacity.

Conclusions

This study provides the most accurate estimation yet of the maintenance-energy coefficient in aerobic cultures of S. cerevisiae, which is a key parameter for modelling of industrial aerobic, glucose-limited fed-batch processes. The observed extreme heat-shock tolerance and high metabolic capacity at near-zero growth rates demonstrate the intrinsic potential of S. cerevisiae as a robust, non-dividing microbial cell factory for energy-intensive products.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0501-z) contains supplementary material, which is available to authorized users.

Degradation Parameters from Pulse-Chase Experiments

Pulse-chase experiments are often used to study the degradation of macromolecules such as proteins or mRNA. Considerations for the choice of pulse length include the toxicity of the pulse to the cell and maximization of labeling. In the general case of non-exponential decay, varying the length of the pulse results in decay patterns that look different. Analysis of these patterns without consideration to pulse length would yield incorrect degradation parameters. Here we propose a method that constructively includes pulse length in the analysis of decay patterns and extracts the parameters of the underlying degradation process. We also show how to extract decay parameters reliably from measurements taken during the pulse phase.

A Minimalistic Resource Allocation Model to Explain Ubiquitous Increase in Protein Expression with Growth Rate

Most proteins show changes in level across growth conditions. Many of these changes seem to be coordinated with the specific growth rate rather than the growth environment or the protein function. Although cellular growth rates, gene expression levels and gene regulation have been at the center of biological research for decades, there are only a few models giving a base line prediction of the dependence of the proteome fraction occupied by a gene with the specific growth rate. We present a simple model that predicts a widely coordinated increase in the fraction of many proteins out of the proteome, proportionally with the growth rate. The model reveals how passive redistribution of resources, due to active regulation of only a few proteins, can have proteome wide effects that are quantitatively predictable. Our model provides a potential explanation for why and how such a coordinated response of a large fraction of the proteome to the specific growth rate arises under different environmental conditions. The simplicity of our model can also be useful by serving as a baseline null hypothesis in the search for active regulation. We exemplify the usage of the model by analyzing the relationship between growth rate and proteome composition for the model microorganism E.coli as reflected in recent proteomics data sets spanning various growth conditions. We find that the fraction out of the proteome of a large number of proteins, and from different cellular processes, increases proportionally with the growth rate. Notably, ribosomal proteins, which have been previously reported to increase in fraction with growth rate, are only a small part of this group of proteins. We suggest that, although the fractions of many proteins change with the growth rate, such changes may be partially driven by a global effect, not necessarily requiring specific cellular control mechanisms.

Protein turnover analysis in Salmonella Typhimurium during infection by dynamic SILAC, Topograph, and quantitative proteomics

Protein turnover affects protein abundance and phenotypes. Comprehensive investigation of protein turnover dynamics has the potential to provide substantial information about gene expression. Here we report a large-scale protein turnover study in Salmonella Typhimurium during infection by quantitative proteomics. Murine macrophage-like RAW 264.7 cells were infected with SILAC labeled Salmonella. Bacterial cells were extracted after 0, 30, 60, 120, and 240 min. Mass spectrometry analyses yielded information about Salmonella protein turnover dynamics and a software program named Topograph was used for the calculation of protein half lives. The half lives of 311 proteins from intracellular Salmonella were obtained. For bacteria cultured in control medium (DMEM), the half lives for 870 proteins were obtained. The calculated median of protein half lives was 69.13 and 99.30 min for the infection group and the DMEM group, respectively, indicating an elevated protein turnover at the initial stage of infection. Gene ontology analyses revealed that a number of protein functional groups were significantly regulated by infection, including proteins involved in ribosome, periplasmic space, cellular amino acid metabolic process, ion binding, and catalytic activity. The half lives of proteins involved in purine metabolism pathway were found to be significantly shortened during infection.

AnGeLi: A Tool for the Analysis of Gene Lists from Fission Yeast

Genome-wide assays and screens typically result in large lists of genes or proteins. Enrichments of functional or other biological properties within such lists can provide valuable insights and testable hypotheses. To systematically detect these enrichments can be challenging and time-consuming, because relevant data to compare against query gene lists are spread over many different sources. We have developed AnGeLi (Analysis of Gene Lists), an intuitive, integrated web-tool for comprehensive and customized interrogation of gene lists from the fission yeast, Schizosaccharomyces pombe. AnGeLi searches for significant enrichments among multiple qualitative and quantitative information sources, including gene and phenotype ontologies, genetic and protein interactions, numerous features of genes, transcripts, translation, and proteins such as copy numbers, chromosomal positions, genetic diversity, RNA polymerase II and ribosome occupancy, localization, conservation, half-lives, domains, and molecular weight among others, as well as diverse sets of genes that are co-regulated or lead to the same phenotypes when mutated. AnGeLi uses robust statistics which can be tailored to specific needs. It also provides the option to upload user-defined gene sets to compare against the query list. Through an integrated data submission form, AnGeLi encourages the community to contribute additional curated gene lists to further increase the usefulness of this resource and to get the most from the ever increasing large-scale experiments. AnGeLi offers a rigorous yet flexible statistical analysis platform for rich insights into functional enrichments and biological context for query gene lists, thus providing a powerful exploratory tool through which S. pombe researchers can uncover fresh perspectives and unexpected connections from genomic data. AnGeLi is freely available at: www.bahlerlab.info/AnGeLi.

Yeast Interspecies Comparative Proteomics Reveals Divergence in Expression Profiles and Provides Insights into Proteome Resource Allocation and Evolutionary Roles of Gene Duplication

Omics analysis is a versatile approach for understanding the conservation and diversity of molecular systems across multiple taxa. In this study, we compared the proteome expression profiles of four yeast species (Saccharomyces cerevisiae, Saccharomyces mikatae, Kluyveromyces waltii, and Kluyveromyces lactis) grown on glucose- or glycerol-containing media. Conserved expression changes across all species were observed only for a small proportion of all proteins differentially expressed between the two growth conditions. Two Kluyveromyces species, both of which exhibited a high growth rate on glycerol, a nonfermentative carbon source, showed distinct species-specific expression profiles. In K. waltii grown on glycerol, proteins involved in the glyoxylate cycle and gluconeogenesis were expressed in high abundance. In K. lactis grown on glycerol, the expression of glycolytic and ethanol metabolic enzymes was unexpectedly low, whereas proteins involved in cytoplasmic translation, including ribosomal proteins and elongation factors, were highly expressed. These marked differences in the types of predominantly expressed proteins suggest that K. lactis optimizes the balance of proteome resource allocation between metabolism and protein synthesis giving priority to cellular growth. In S. cerevisiae, about 450 duplicate gene pairs were retained after whole-genome duplication. Intriguingly, we found that in the case of duplicates with conserved sequences, the total abundance of proteins encoded by a duplicate pair in S. cerevisiae was similar to that of protein encoded by nonduplicated ortholog in Kluyveromyces yeast. Given the frequency of haploinsufficiency, this observation suggests that conserved duplicate genes, even though minor cases of retained duplicates, do not exhibit a dosage effect in yeast, except for ribosomal proteins. Thus, comparative proteomic analyses across multiple species may reveal not only species-specific characteristics of metabolic processes under nonoptimal culture conditions but also provide valuable insights into intriguing biological principles, including the balance of proteome resource allocation and the role of gene duplication in evolutionary history.

Next-generation analysis of gene expression regulation--comparing the roles of synthesis and degradation

Technological advances now enable routine measurement of mRNA and protein abundances, and estimates of their rates of synthesis and degradation that inform on their values and the degree of change in response to stimuli. Importantly, more and more data on time-series experiments are emerging, e.g. of cells responding to stress, enabling first insights into a new dimension of gene expression regulation - its dynamics and how it allows for very different response signals across genes. This review discusses recently published methods and datasets, their impact on what we now know about the relationships between concentrations and synthesis rates of mRNAs and proteins in yeast and mammalian cells, their evolution, and new hypotheses on translation regulatory mechanisms generated by approaches that involve ribosome footprinting.

Proteolytic regulation of metabolic enzymes by E3 ubiquitin ligase complexes: lessons from yeast

Eukaryotic organisms use diverse mechanisms to control metabolic rates in response to changes in the internal and/or external environment. Fine metabolic control is a highly responsive, energy-saving process that is mediated by allosteric inhibition/activation and/or reversible modification of preexisting metabolic enzymes. In contrast, coarse metabolic control is a relatively long-term and expensive process that involves modulating the level of metabolic enzymes. Coarse metabolic control can be achieved through the degradation of metabolic enzymes by the ubiquitin-proteasome system (UPS), in which substrates are specifically ubiquitinated by an E3 ubiquitin ligase and targeted for proteasomal degradation. Here, we review select multi-protein E3 ligase complexes that directly regulate metabolic enzymes in Saccharomyces cerevisiae. The first part of the review focuses on the endoplasmic reticulum (ER) membrane-associated Hrd1 and Doa10 E3 ligase complexes. In addition to their primary roles in the ER-associated degradation pathway that eliminates misfolded proteins, recent quantitative proteomic analyses identified native substrates of Hrd1 and Doa10 in the sterol synthesis pathway. The second part focuses on the SCF (Skp1-Cul1-F-box protein) complex, an abundant prototypical multi-protein E3 ligase complex. While the best-known roles of the SCF complex are in the regulation of the cell cycle and transcription, accumulating evidence indicates that the SCF complex also modulates carbon metabolism pathways. The increasing number of metabolic enzymes whose stability is directly regulated by the UPS underscores the importance of the proteolytic regulation of metabolic processes for the acclimation of cells to environmental changes.

Feasibility of protein turnover studies in prototroph Saccharomyces cerevisiae strains

Quantitative proteomics studies of yeast that use metabolic labeling with amino acids rely on auxotrophic mutations of one or more genes on the amino acid biosynthesis pathways. These mutations affect yeast metabolism and preclude the study of some biological processes. Overcoming this limitation, it has recently been described that proteins in a yeast prototrophic strain can also be metabolically labeled with heavy amino acids. However, the temporal profiles of label incorporation under the different phases of the prototroph's growth have not been examined. Labeling trajectories are important in the study of protein turnover and dynamics, in which label incorporation into proteins is monitored across many time points. Here we monitored protein labeling trajectories for 48 h after a pulse with heavy lysine in a yeast prototrophic strain and compared them with those of a lysine auxotrophic yeast. Labeling was successful in prototroph yeast during exponential growth phase but not in stationary phase. Furthermore, we were able to determine the half-lives of more than 1700 proteins during exponential phase of growth with high accuracy and reproducibility. We found a median half-life of 2 h in both strains, which corresponds with the cellular doubling time. Nucleolar and ribosomal proteins showed short half-lives, whereas mitochondrial proteins and other energy production enzymes presented longer half-lives. Except for some proteins involved in lysine biosynthesis, we observed a high correlation in protein half-lives between prototroph and auxotroph strains. Overall, our results demonstrate the feasibility of using prototrophs for proteomic turnover studies and provide a reliable data set of protein half-lives in exponentially growing yeast.

Genome-wide screen identifies a novel p97/CDC-48-dependent pathway regulating ER-stress-induced gene transcription

The accumulation of misfolded proteins in the endoplasmic reticulum (ER) activates the Unfolded Protein Response (UPR(ER)) to restore ER homeostasis. The AAA(+) ATPase p97/CDC-48 plays key roles in ER stress by promoting both ER protein degradation and transcription of UPR(ER) genes. Although the mechanisms associated with protein degradation are now well established, the molecular events involved in the regulation of gene transcription by p97/CDC-48 remain unclear. Using a reporter-based genome-wide RNAi screen in combination with quantitative proteomic analysis in Caenorhabditis elegans, we have identified RUVB-2, a AAA(+) ATPase, as a novel repressor of a subset of UPR(ER) genes. We show that degradation of RUVB-2 by CDC-48 enhances expression of ER stress response genes through an XBP1-dependent mechanism. The functional interplay between CDC-48 and RUVB-2 in controlling transcription of select UPR(ER) genes appears conserved in human cells. Together, these results describe a novel role for p97/CDC-48, whereby its role in protein degradation is integrated with its role in regulating expression of ER stress response genes.