GroEL dependency affects codon usage--support for a critical role of misfolding in gene evolution

Gene product diversity: adaptive or not?

One gene does not equal one RNA or protein. The genomic revolution has revealed numerous different RNA and protein molecules that can be produced from one gene, such as circular RNAs generated by back-splicing, proteins with residues mismatching the genomic encoding because of RNA editing, and proteins extended in the C terminus via stop codon readthrough in translation. Are these diverse products results of exquisite gene regulations or imprecise biological processes? While there are cases where the gene product diversity appears beneficial, genome-scale patterns suggest that much of this diversity arises from nonadaptive, molecular errors. This finding has important implications for studying the functions of diverse gene products and for understanding the fundamental properties and evolution of cellular life.

Chaperone requirements for de novo folding of Saccharomyces cerevisiae septins

Polymers of septin protein complexes play cytoskeletal roles in eukaryotic cells. The specific subunit composition within complexes controls functions and higher-order structural properties. All septins have globular GTPase domains. The other eukaryotic cytoskeletal NTPases strictly require assistance from molecular chaperones of the cytosol, particularly the cage-like chaperonins, to fold into oligomerization-competent conformations. We previously identified cytosolic chaperones that bind septins and influence the oligomerization ability of septins carrying mutations linked to human disease, but it was unknown to what extent wild-type septins require chaperone assistance for their native folding. Here we use a combination of in vivo and in vitro approaches to demonstrate chaperone requirements for de novo folding and complex assembly by budding yeast septins. Individually purified septins adopted nonnative conformations and formed nonnative homodimers. In chaperonin- or Hsp70-deficient cells, septins folded slower and were unable to assemble posttranslationally into native complexes. One septin, Cdc12, was so dependent on cotranslational chaperonin assistance that translation failed without it. Our findings point to distinct translation elongation rates for different septins as a possible mechanism to direct a stepwise, cotranslational assembly pathway in which general cytosolic chaperones act as key intermediaries.

Network analysis of synonymous codon usage

MOTIVATION

Most amino acids are encoded by multiple synonymous codons, some of which are used more rarely than others. Analyses of positions of such rare codons in protein sequences revealed that rare codons can impact co-translational protein folding and that positions of some rare codons are evolutionarily conserved. Analyses of their positions in protein 3-dimensional structures, which are richer in biochemical information than sequences alone, might further explain the role of rare codons in protein folding.

RESULTS

We model protein structures as networks and use network centrality to measure the structural position of an amino acid. We first validate that amino acids buried within the structural core are network-central, and those on the surface are not. Then, we study potential differences between network centralities and thus structural positions of amino acids encoded by conserved rare, non-conserved rare and commonly used codons. We find that in 84% of proteins, the three codon categories occupy significantly different structural positions. We examine protein groups showing different codon centrality trends, i.e. different relationships between structural positions of the three codon categories. We see several cases of all proteins from our data with some structural or functional property being in the same group. Also, we see a case of all proteins in some group having the same property. Our work shows that codon usage is linked to the final protein structure and thus possibly to co-translational protein folding.

AVAILABILITY AND IMPLEMENTATION

https://nd.edu/∼cone/CodonUsage/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Widespread non-modular overlapping codes in the coding regions

Messenger RNAs (mRNAs) consist of a coding region (open reading frame (ORF)) and two untranslated regions (UTRs), 5'UTR and 3'UTR. Ribosomes travel along the coding region, translating nucleotide triplets (called codons) to a chain of amino acids. The coding region was long believed to mainly encode the amino acid content of proteins, whereas regulatory signals reside in the UTRs and in other genomic regions. However, in recent years we have learned that the ORF is expansively populated with various regulatory signals, or codes, which are related to all gene expression steps and additional intracellular aspects. In this paper, we review the current knowledge related to overlapping codes inside the coding regions, such as the influence of synonymous codon usage on translation speed (and, in turn, the effect of translation speed on protein folding), ribosomal frameshifting, mRNA stability, methylation, splicing, transcription and more. All these codes come together and overlap in the ORF sequence, ensuring production of the right protein at the right time.

The stop-and-go traffic regulating protein biogenesis: How translation kinetics controls proteostasis

Generating a functional proteome requires the ribosome to carefully regulate disparate co-translational processes that determine the fate of nascent polypeptides. With protein synthesis being energetically expensive, the ribosome must balance the costs of efficiently making a protein with those of properly folding it. Emerging as a primary means of regulating this trade-off is the nonuniform rate of translation elongation that defines translation kinetics. The varying speeds with which the ribosome progresses along a transcript have been implicated in several aspects of protein biogenesis, including co-translational protein folding and translational fidelity, as well as gene expression by mediating mRNA decay and protein quality control pathways. The optimal translation kinetics required to efficiently execute these processes can be distinct. Thus, the ribosome is tasked with tightly regulating translation kinetics to balance these processes while maintaining adaptability for changing cellular conditions. In this review, we first discuss the regulatory role of translation elongation in protein biogenesis and what factors influence elongation kinetics. We then describe how changes in translation kinetics signal downstream pathways that dictate the fate of nascent polypeptides. By regulating these pathways, the kinetics of translation elongation has emerged as a critical tool for driving gene expression and maintaining proteostasis through varied mechanisms, including nascent chain folding and binding different ribosome-associated machinery. Indeed, a growing number of examples demonstrate the important role of local changes in elongation kinetics in modulating the pathophysiology of human disease.

Codon usage regulates protein structure and function by affecting translation elongation speed in Drosophila cells

Codon usage biases are found in all eukaryotic and prokaryotic genomes and have been proposed to regulate different aspects of translation process. Codon optimality has been shown to regulate translation elongation speed in fungal systems, but its effect on translation elongation speed in animal systems is not clear. In this study, we used a Drosophila cell-free translation system to directly compare the velocity of mRNA translation elongation. Our results demonstrate that optimal synonymous codons speed up translation elongation while non-optimal codons slow down translation. In addition, codon usage regulates ribosome movement and stalling on mRNA during translation. Finally, we show that codon usage affects protein structure and function in vitro and in Drosophila cells. Together, these results suggest that the effect of codon usage on translation elongation speed is a conserved mechanism from fungi to animals that can affect protein folding in eukaryotic organisms.

Evolution of Chaperonin Gene Duplication in Stigonematalean Cyanobacteria (Subsection V)

Chaperonins promote protein folding and are known to play a role in the maintenance of cellular stability under stress conditions. The group I bacterial chaperonin complex comprises GroEL, that forms a barrel-like oligomer, and GroES that forms the lid. In most eubacteria the GroES/GroEL chaperonin is encoded by a single-copy bicistronic operon, whereas in cyanobacteria up to three groES/groEL paralogs have been documented. Here we study the evolution and functional diversification of chaperonin paralogs in the heterocystous, multi-seriate filament forming cyanobacterium Chlorogloeopsis fritschii PCC 6912. The genome of C. fritschii encodes two groES/groEL operons (groESL1, groESL1.2) and a monocistronic groEL gene (groEL2). A phylogenetic reconstruction reveals that the groEL2 duplication is as ancient as cyanobacteria, whereas the groESL1.2 duplication occurred at the ancestor of heterocystous cyanobacteria. A comparison of the groEL paralogs transcription levels under different growth conditions shows that they have adapted distinct transcriptional regulation. Our results reveal that groEL1 and groEL1.2 are upregulated during diazotrophic conditions and the localization of their promoter activity points towards a role in heterocyst differentiation. Furthermore, protein–protein interaction assays suggest that paralogs encoded in the two operons assemble into hybrid complexes. The monocistronic encoded GroEL2 is not forming oligomers nor does it interact with the co-chaperonins. Interaction between GroES1.2 and GroEL1.2 could not be documented, suggesting that the groESL1.2 operon does not encode a functional chaperonin complex. Functional complementation experiments in Escherichia coli show that only GroES1/GroEL1 and GroES1/GroEL1.2 can substitute the native operon. In summary, the evolutionary consequences of chaperonin duplication in cyanobacteria include the retention of groESL1 as a housekeeping gene, subfunctionalization of groESL1.2 and neofunctionalization of the monocistronic groEL2 paralog.

Codon optimality, bias and usage in translation and mRNA decay

The advent of ribosome profiling and other tools to probe mRNA translation has revealed that codon bias - the uneven use of synonymous codons in the transcriptome - serves as a secondary genetic code: a code that guides the efficiency of protein production, the fidelity of translation and the metabolism of mRNAs. Recent advancements in our understanding of mRNA decay have revealed a tight coupling between ribosome dynamics and the stability of mRNA transcripts; this coupling integrates codon bias into the concept of codon optimality, or the effects that specific codons and tRNA concentrations have on the efficiency and fidelity of the translation machinery. In this Review, we first discuss the evidence for codon-dependent effects on translation, beginning with the basic mechanisms through which translation perturbation can affect translation efficiency, protein folding and transcript stability. We then discuss how codon effects are leveraged by the cell to tailor the proteome to maintain homeostasis, execute specific gene expression programmes of growth or differentiation and optimize the efficiency of protein production.

The Molecular Chaperone DnaK Is a Source of Mutational Robustness

Molecular chaperones, also known as heat-shock proteins, refold misfolded proteins and help other proteins reach their native conformation. Thanks to these abilities, some chaperones, such as the Hsp90 protein or the chaperonin GroEL, can buffer the deleterious phenotypic effects of mutations that alter protein structure and function. Hsp70 chaperones use a chaperoning mechanism different from that of Hsp90 and GroEL, and it is not known whether they can also buffer mutations. Here, we show that they can. To this end, we performed a mutation accumulation experiment in Escherichia coli, followed by whole-genome resequencing. Overexpression of the Hsp70 chaperone DnaK helps cells cope with mutational load and completely avoid the extinctions we observe in lineages evolving without chaperone overproduction. Additionally, our sequence data show that DnaK overexpression increases mutational robustness, the tolerance of its clients to nonsynonymous nucleotide substitutions. We also show that this elevated mutational buffering translates into differences in evolutionary rates on intermediate and long evolutionary time scales. Specifically, we studied the evolutionary rates of DnaK clients using the genomes of E. coli, Salmonella enterica, and 83 other gamma-proteobacteria. We find that clients that interact strongly with DnaK evolve faster than weakly interacting clients. Our results imply that all three major chaperone classes can buffer mutations and affect protein evolution. They illustrate how an individual protein like a chaperone can have a disproportionate effect on the evolution of a proteome.

Trm9-Catalyzed tRNA Modifications Regulate Global Protein Expression by Codon-Biased Translation

Post-transcriptional modifications of transfer RNAs (tRNAs) have long been recognized to play crucial roles in regulating the rate and fidelity of translation. However, the extent to which they determine global protein production remains poorly understood. Here we use quantitative proteomics to show a direct link between wobble uridine 5-methoxycarbonylmethyl (mcm⁵) and 5-methoxy-carbonyl-methyl-2-thio (mcm⁵s²) modifications catalyzed by tRNA methyltransferase 9 (Trm9) in tRNA^Arg(UCU) and tRNA^Glu(UUC) and selective translation of proteins from genes enriched with their cognate codons. Controlling for bias in protein expression and alternations in mRNA expression, we find that loss of Trm9 selectively impairs expression of proteins from genes enriched with AGA and GAA codons under both normal and stress conditions. Moreover, we show that AGA and GAA codons occur with high frequency in clusters along the transcripts, which may play a role in modulating translation. Consistent with these results, proteins subject to enhanced ribosome pausing in yeast lacking mcm⁵U and mcm⁵s²U are more likely to be down-regulated and contain a larger number of AGA/GAA clusters. Together, these results suggest that Trm9-catalyzed tRNA modifications play a significant role in regulating protein expression within the cell.

Here we present evidence for a more complicated role for transfer RNAs (tRNAs) than as mere adapters that link the genetic code in messenger RNA (mRNA) to the amino acid sequence of a protein during translation. tRNAs have long been known to be modified with dozens of different chemical structures other than the 4 canonical ribonucleosides, though the role of these modifications in controlling translation is poorly understood. By quantifying the expression of thousands of proteins in the yeast S. cerevisiae, we identified a mechanistic link between modified ribonucleosides located at the wobble position of two tRNAs, tRNA^Arg(UCU) and tRNA^Glu(UUC), and the translation of proteins derived from genes enriched with codons read by these tRNAs: AGA and GAA. In cells lacking the enzyme that inserts these modifications, tRNA methyltransferase 9 (Trm9), we found a significant reduction in proteins from genes enriched with AGA and GAA codons and with runs of these codons. Also, mRNAs enriched with runs of AGA and GAA codons are subject to stalled translation on ribosomes in yeast lacking mcm⁵U and mcm⁵s²U. Together, these results reveal a distinct role for Trm9-catalyzed tRNA modifications in selectively regulating the expression of proteins enriched with AGA and GAA codons.

Cold adaptation of a psychrophilic chaperonin from Psychrobacter sp. and its application for heterologous protein expression

OBJECTIVES

A chaperonin, PsyGroELS, from the Antarctic psychrophilic bacterium Psychrobacter sp. PAMC21119, was examined for its role in cold adaptation when expressed in a mesophilic Escherichia coli strain.

RESULTS

Growth of E. coli harboring PsyGroELS at 10 °C was increased compared to the control strain. A co-expression system using PsyGroELS was developed to increase productivity of the psychrophilic enzyme PsyEst9. PsyEst9 was cloned and expressed using three E. coli variants that co-expressed GroELS from PAMC21119, E. coli, or Oleispira antarctica RB8(T). Co-expression with PsyGroELS was more effective for the production of PsyEst9 compared tothe other chaperonins.

CONCLUSION

PsyGroELS confers cold tolerance to E. coli, and shows potential as an effective co-expression system for the stable production of psychrophilic proteins.

Biophysics of protein evolution and evolutionary protein biophysics

The study of molecular evolution at the level of protein-coding genes often entails comparing large datasets of sequences to infer their evolutionary relationships. Despite the importance of a protein's structure and conformational dynamics to its function and thus its fitness, common phylogenetic methods embody minimal biophysical knowledge of proteins. To underscore the biophysical constraints on natural selection, we survey effects of protein mutations, highlighting the physical basis for marginal stability of natural globular proteins and how requirement for kinetic stability and avoidance of misfolding and misinteractions might have affected protein evolution. The biophysical underpinnings of these effects have been addressed by models with an explicit coarse-grained spatial representation of the polypeptide chain. Sequence-structure mappings based on such models are powerful conceptual tools that rationalize mutational robustness, evolvability, epistasis, promiscuous function performed by 'hidden' conformational states, resolution of adaptive conflicts and conformational switches in the evolution from one protein fold to another. Recently, protein biophysics has been applied to derive more accurate evolutionary accounts of sequence data. Methods have also been developed to exploit sequence-based evolutionary information to predict biophysical behaviours of proteins. The success of these approaches demonstrates a deep synergy between the fields of protein biophysics and protein evolution.

RNA chaperones buffer deleterious mutations in E. coli

Both proteins and RNAs can misfold into non-functional conformations. Protein chaperones promote native folding of nascent polypeptides and refolding of misfolded species, thereby buffering mutations that compromise protein structure and function. Here, we show that RNA chaperones can also act as mutation buffers that enhance organismal fitness. Using competition assays, we demonstrate that overexpression of select RNA chaperones, including three DEAD box RNA helicases (DBRHs) (CsdA, SrmB, RhlB) and the cold shock protein CspA, improves fitness of two independently evolved Escherichia coli mutator strains that have accumulated deleterious mutations during short- and long-term laboratory evolution. We identify strain-specific mutations that are deleterious and subject to buffering when introduced individually into the ancestral genotype. For DBRHs, we show that buffering requires helicase activity, implicating RNA structural remodelling in the buffering process. Our results suggest that RNA chaperones might play a fundamental role in RNA evolution and evolvability.

Hsp90 promotes kinase evolution

Heat-shock protein 90 (Hsp90) promotes the maturation and stability of its client proteins, including many kinases. In doing so, Hsp90 may allow its clients to accumulate mutations as previously proposed by the capacitor hypothesis. If true, Hsp90 clients should show increased evolutionary rate compared with nonclients; however, other factors, such as gene expression and protein connectivity, may confound or obscure the chaperone’s putative contribution. Here, we compared the evolutionary rates of many Hsp90 clients and nonclients in the human protein kinase superfamily. We show that Hsp90 client status promotes evolutionary rate independently of, but in a small magnitude similar to that of gene expression and protein connectivity. Hsp90’s effect on kinase evolutionary rate was detected across mammals, specifically relaxing purifying selection. Hsp90 clients also showed increased nucleotide diversity and harbored more damaging variation than nonclient kinases across humans. These results are consistent with the central argument of the capacitor hypothesis that interaction with the chaperone allows its clients to harbor genetic variation. Hsp90 client status is thought to be highly dynamic with as few as one amino acid change rendering a protein dependent on the chaperone. Contrary to this expectation, we found that across protein kinase phylogeny Hsp90 client status tends to be gained, maintained, and shared among closely related kinases. We also infer that the ancestral protein kinase was not an Hsp90 client. Taken together, our results suggest that Hsp90 played an important role in shaping the kinase superfamily.

Exposing synonymous mutations

Synonymous codon changes, which do not alter protein sequence, were previously thought to have no functional consequence. Although this concept has been overturned in recent years, there is no unique mechanism by which these changes exert biological effects. A large repertoire of both experimental and bioinformatic methods has been developed to understand the effects of synonymous variants. Results from this body of work have provided global insights into how biological systems exploit the degeneracy of the genetic code to control gene expression, protein folding efficiency, and the coordinated expression of functionally related gene families. Although it is now clear that synonymous variants are important in a variety of contexts, from human disease to the safety and efficacy of therapeutic proteins, there is no clear consensus on the approaches to identify and validate these changes. Here, we review the diverse methods to understand the effects of synonymous mutations.

Entrainment to periodic initiation and transition rates in a computational model for gene translation

Periodic oscillations play an important role in many biomedical systems. Proper functioning of biological systems that respond to periodic signals requires the ability to synchronize with the periodic excitation. For example, the sleep/wake cycle is a manifestation of an internal timing system that synchronizes to the solar day. In the terminology of systems theory, the biological system must entrain or phase-lock to the periodic excitation. Entrainment is also important in synthetic biology. For example, connecting several artificial biological systems that entrain to a common clock may lead to a well-functioning modular system. The cell-cycle is a periodic program that regulates DNA synthesis and cell division. Recent biological studies suggest that cell-cycle related genes entrain to this periodic program at the gene translation level, leading to periodically-varying protein levels of these genes. The ribosome flow model (RFM) is a deterministic model obtained via a mean-field approximation of a stochastic model from statistical physics that has been used to model numerous processes including ribosome flow along the mRNA. Here we analyze the RFM under the assumption that the initiation and/or transition rates vary periodically with a common period . We show that the ribosome distribution profile in the RFM entrains to this periodic excitation. In particular, the protein synthesis pattern converges to a unique periodic solution with period . To the best of our knowledge, this is the first proof of entrainment in a mathematical model for translation that encapsulates aspects such as initiation and termination rates, ribosomal movement and interactions, and non-homogeneous elongation speeds along the mRNA. Our results support the conjecture that periodic oscillations in tRNA levels and other factors related to the translation process can induce periodic oscillations in protein levels, and may suggest a new approach for re-engineering genetic systems to obtain a desired, periodic, protein synthesis rate.

Cyanobacterial heat-shock response: role and regulation of molecular chaperones

Cyanobacteria constitute a morphologically diverse group of oxygenic photoautotrophic microbes which range from unicellular to multicellular, and non-nitrogen-fixing to nitrogen-fixing types. Sustained long-term exposure to changing environmental conditions, during their three billion years of evolution, has presumably led to their adaptation to diverse ecological niches. The ability to maintain protein conformational homeostasis (folding–misfolding–refolding or aggregation–degradation) by molecular chaperones holds the key to the stress adaptability of cyanobacteria. Although cyanobacteria possess several genes encoding DnaK and DnaJ family proteins, these are not the most abundant heat-shock proteins (Hsps), as is the case in other bacteria. Instead, the Hsp60 family of proteins, comprising two phylogenetically conserved proteins, and small Hsps are more abundant during heat stress. The contribution of the Hsp100 (ClpB) family of proteins and of small Hsps in the unicellular cyanobacteria (Synechocystis and Synechococcus) as well as that of Hsp60 proteins in the filamentous cyanobacteria (Anabaena) to thermotolerance has been elucidated. The regulation of chaperone genes by several cis-elements and trans-acting factors has also been well documented. Recent studies have demonstrated novel transcriptional and translational (mRNA secondary structure) regulatory mechanisms in unicellular cyanobacteria. This article provides an insight into the heat-shock response: its organization, and ecophysiological regulation and role of molecular chaperones, in unicellular and filamentous nitrogen-fixing cyanobacterial strains.

Flagellated algae protein evolution suggests the prevalence of lineage-specific rules governing evolutionary rates of eukaryotic proteins

Understanding the general rules governing the rate of protein evolution is fundamental to evolutionary biology. However, attempts to address this issue in yeasts and mammals have revealed considerable differences in the relative importance of determinants for protein evolutionary rates. This phenomenon was previously explained by the fact that yeasts and mammals are different in many cellular and genomic properties. Flagellated algae species have several cellular and genomic characteristics that are intermediate between yeasts and mammals. Using partial correlation analyses on the evolution of 6,921 orthologous proteins from Chlamydomonas reinhardtii and Volvox carteri, we examined factors influencing evolutionary rates of proteins in flagellated algae. Previous studies have shown that mRNA abundance and gene compactness are strong determinants for protein evolutionary rates in yeasts and mammals, respectively. We show that both factors also influence algae protein evolution with mRNA abundance having a larger impact than gene compactness on the rates of algae protein evolution. More importantly, among all the factors examined, coding sequence (CDS) length has the strongest (positive) correlation with protein evolutionary rates. This correlation between CDS length and the rates of protein evolution is not due to alignment-related issues or domain density. These results suggest no simple and universal rules governing protein evolutionary rates across different eukaryotic lineages. Instead, gene properties influence the rate of protein evolution in a lineage-specific manner.

New universal rules of eukaryotic translation initiation fidelity

The accepted model of eukaryotic translation initiation begins with the scanning of the transcript by the pre-initiation complex from the 5′end until an ATG codon with a specific nucleotide (nt) context surrounding it is recognized (Kozak rule). According to this model, ATG codons upstream to the beginning of the ORF should affect translation. We perform for the first time, a genome-wide statistical analysis, uncovering a new, more comprehensive and quantitative, set of initiation rules for improving the cost of translation and its efficiency. Analyzing dozens of eukaryotic genomes, we find that in all frames there is a universal trend of selection for low numbers of ATG codons; specifically, 16–27 codons upstream, but also 5–11 codons downstream of the START ATG, include less ATG codons than expected. We further suggest that there is selection for anti optimal ATG contexts in the vicinity of the START ATG. Thus, the efficiency and fidelity of translation initiation is encoded in the 5′UTR as required by the scanning model, but also at the beginning of the ORF.

The observed nt patterns suggest that in all the analyzed organisms the pre-initiation complex often misses the START ATG of the ORF, and may start translation from an alternative initiation start-site. Thus, to prevent the translation of undesired proteins, there is selection for nucleotide sequences with low affinity to the pre-initiation complex near the beginning of the ORF. With the new suggested rules we were able to obtain a twice higher correlation with ribosomal density and protein levels in comparison to the Kozak rule alone (e.g. for protein levels r = 0.7 vs. r = 0.31; p<10⁻¹²).

Gene translation is an important step of the intra-cellular protein synthesis, which is a central process in all living organisms. Thus, understanding how translation efficiency is encoded in transcripts has ramifications to every biomedical discipline. The aim of the current study is to decipher the way translation initiation fidelity is encoded in eukaryotic transcripts, and how evolution shapes the beginning of transcripts. Based on the genomes of dozens of organisms we were able to derive a new, more precise, set of rules related to this process, facilitating a high resolution view of the mechanisms aiding translation initiation fidelity. Among others, we show that there is a universal trend of selection for low numbers of ATG codons upstream, but also in the 5–11 codons downstream of the START ATG, presumably to prevent translation of alternative ORFs over the main one. With the new suggested rules we were able to obtain a twice higher correlation with ribosomal density and protein levels in comparison to the previous translation initiation efficiency rule.

Efficient translation initiation dictates codon usage at gene start

Rare codons are enriched at gene start in many genomes. Genome analysis and experimental testing show that this enrichment evolved to keep the ribosome binding site free from stable mRNA structures, in order to facilitate efficient translation initiation.

The use of rare codons coincides with suppression of mRNA structures at the ribosome binding site across genomes.

There is preferential selection for synonymous codons that reduce GC-content at the beginning of genes and a stronger pressure for rare codon usage in GC-rich organisms.

Amino acids encoded by AU-rich codons are preferred at gene start.

Experimental results show that mRNA structure at translation start strongly influences protein yield.

The genetic code is degenerate; thus, protein evolution does not uniquely determine the coding sequence. One of the puzzles in evolutionary genetics is therefore to uncover evolutionary driving forces that result in specific codon choice. In many bacteria, the first 5–10 codons of protein-coding genes are often codons that are less frequently used in the rest of the genome, an effect that has been argued to arise from selection for slowed early elongation to reduce ribosome traffic jams. However, genome analysis across many species has demonstrated that the region shows reduced mRNA folding consistent with pressure for efficient translation initiation. This raises the possibility that unusual codon usage is a side effect of selection for reduced mRNA structure. Here we discriminate between these two competing hypotheses, and show that in bacteria selection favours codons that reduce mRNA folding around the translation start, regardless of whether these codons are frequent or rare. Experiments confirm that primarily mRNA structure, and not codon usage, at the beginning of genes determines the translation rate.

Positively charged residues are the major determinants of ribosomal velocity

Both for understanding mechanisms of disease and for the design of transgenes, it is important to understand the determinants of ribosome velocity, as changes in the rate of translation are important for protein folding, error attenuation, and localization. While there is great variation in ribosomal occupancy along even a single transcript, what determines a ribosome's occupancy is unclear. We examine this issue using data from a ribosomal footprinting assay in yeast. While codon usage is classically considered a major determinant, we find no evidence for this. By contrast, we find that positively charged amino acids greatly retard ribosomes downstream from where they are encoded, consistent with the suggestion that positively charged residues interact with the negatively charged ribosomal exit tunnel. Such slowing is independent of and greater than the average effect owing to mRNA folding. The effect of charged amino acids is additive, with ribosomal occupancy well-predicted by a linear fit to the density of positively charged residues. We thus expect that a translated poly-A tail, encoding for positively charged lysines regardless of the reading frame, would act as a sandtrap for the ribosome, consistent with experimental data.

Synonymous codon changes in the oncogenes of the cottontail rabbit papillomavirus lead to increased oncogenicity and immunogenicity of the virus

Papillomaviruses use rare codons with respect to the host. The reasons for this are incompletely understood but among the hypotheses is the concept that rare codons result in low protein production and this allows the virus to escape immune surveillance. We changed rare codons in the oncogenes E6 and E7 of the cottontail rabbit papillomavirus to make them more mammalian-like and tested the mutant genomes in our in vivo animal model. While the amino acid sequences of the proteins remained unchanged, the oncogenic potential of some of the altered genomes increased dramatically. In addition, increased immunogenicity, as measured by spontaneous regression, was observed as the numbers of codon changes increased. This work suggests that codon usage may modify protein production in ways that influence disease outcome and that evaluation of synonymous codons should be included in the analysis of genetic variants of infectious agents and their association with disease.

Evidence for deep phylogenetic conservation of exonic splice-related constraints: splice-related skews at exonic ends in the brown alga Ectocarpus are common and resemble those seen in humans

The control of RNA splicing is often modulated by exonic motifs near splice sites. Chief among these are exonic splice enhancers (ESEs). Well-described ESEs in mammals are purine rich and cause predictable skews in codon and amino acid usage toward exonic ends. Looking across species, those with relatively abundant intronic sequence are those with the more profound end of exon skews, indicative of exonization of splice site recognition. To date, the only intron-rich species that have been analyzed are mammals, precluding any conclusions about the likely ancestral condition. Here, we examine the patterns of codon and amino acid usage in the vicinity of exon-intron junctions in the brown alga Ectocarpus siliculosus, a species with abundant large introns, known SR proteins, and classical splice sites. We find that amino acids and codons preferred/avoided at both 3' and 5' ends in Ectocarpus, of which there are many, tend, on average, to also be preferred/avoided at the same exon ends in humans. Moreover, the preferences observed at the 5' ends of exons are largely the same as those at the 3' ends, a symmetry trend only previously observed in animals. We predict putative hexameric ESEs in Ectocarpus and show that these are purine rich and that there are many more of these identified as functional ESEs in humans than expected by chance. These results are consistent with deep phylogenetic conservation of SR protein binding motifs. Assuming codons preferred near boundaries are "splice optimal" codons, in Ectocarpus, unlike Drosophila, splice optimal and translationally optimal codons are not mutually exclusive. The exclusivity of translationally optimal and splice optimal codon sets is thus not universal.

Evolutionary conservation of codon optimality reveals hidden signatures of cotranslational folding

The choice of codons can influence local translation kinetics during protein synthesis. Whether codon preference is linked to co-translational regulation of polypeptide folding remains unclear. Here, we derive a revised translational efficiency scale that incorporates the competition between tRNA supply and demand. Applying this scale to ten closely related yeasts, we uncover the evolutionary conservation of codon optimality in eukaryotes. This analysis reveals universal patterns of conserved optimal and nonoptimal codons, often in clusters, which associate with the secondary structure of the translated polypeptides independent of the levels of expression. Our analysis suggests an evolved function for codon optimality in regulating the rhythm of elongation to facilitate co-translational polypeptide folding, beyond its previously proposed role of adapting to the cost of expression. These findings establish how mRNA sequences are generally under selection to optimize the co-translational folding of corresponding polypeptides.

Cumulative impact of chaperone-mediated folding on genome evolution

Molecular chaperones support protein folding and unfolding along with assembly and translocation of protein complexes. Chaperones have been recognized as important mediators between an organismal genotype and phenotype as well as important maintainers of cellular fitness under environmental conditions that induce high mutational loads. Here we review recent studies revealing that the folding assistance supplied by chaperones is evident in genomic sequences implicating chaperone-mediated folding as an influential factor during protein evolution. Interaction of protein with chaperones ensures a proper folding and function, yet an adaptation to obligatory dependence on such assistance may be irreversible, representing an evolutionary trap. A correlation between the requirement for a chaperone and protein expression level indicates that the evolution of substrate-chaperone interaction is bounded by the required substrate abundance within the cell. Accumulating evidence suggests that the utility of chaperones is governed by a delicate balance between their help in mitigating the risks of protein misfolding and aggregate formation on one hand and the slower rate of protein maturation and the energetic cost of chaperone synthesis on the other.

On the steady-state distribution in the homogeneous ribosome flow model

A central biological process in all living organisms is gene translation. Developing a deeper understanding of this complex process may have ramifications to almost every biomedical discipline. Reuveni et al. recently proposed a new computational model of gene translation called the Ribosome Flow Model (RFM). In this paper, we consider a particular case of this model, called the Homogeneous Ribosome Flow Model (HRFM). From a biological viewpoint, this corresponds to the case where the transition rates of all the coding sequence codons are identical. This regime has been suggested recently based on experiments in mouse embryonic cells. We consider the steady-state distribution of the HRFM. We provide formulas that relate the different parameters of the model in steady state. We prove the following properties: 1) the ribosomal density profile is monotonically decreasing along the coding sequence; 2) the ribosomal density at each codon monotonically increases with the initiation rate; and 3) for a constant initiation rate, the translation rate monotonically decreases with the length of the coding sequence. In addition, we analyze the translation rate of the HRFM at the limit of very high and very low initiation rate, and provide explicit formulas for the translation rate in these two cases. We discuss the relationship between these theoretical results and biological findings on the translation process.

Evolutionary layering and the limits to cellular perfection

Although observations from biochemistry and cell biology seemingly illustrate hundreds of examples of exquisite molecular adaptations, the fact that experimental manipulation can often result in improvements in cellular infrastructure raises the question as to what ultimately limits the level of molecular perfection achievable by natural selection. Here, it is argued that random genetic drift can impose a strong barrier to the advancement of molecular refinements by adaptive processes. Moreover, although substantial improvements in fitness may sometimes be accomplished via the emergence of novel cellular features that improve on previously established mechanisms, such advances are expected to often be transient, with overall fitness eventually returning to the level before incorporation of the genetic novelty. As a consequence of such changes, increased molecular/cellular complexity can arise by Darwinian processes, while yielding no long-term increase in adaptation and imposing increased energetic and mutational costs.

Stability analysis of the ribosome flow model

Gene translation is a central process in all living organisms. Developing a better understanding of this complex process may have ramifications to almost every biomedical discipline. Recently, Reuveni et al. proposed a new computational model of this process called the ribosome flow model (RFM). In this study, we show that the dynamical behavior of the RFM is relatively simple. There exists a unique equilibrium point e and every trajectory converges to e. Furthermore, convergence is monotone in the sense that the distance to e can never increase. This qualitative behavior is maintained for any feasible set of parameter values, suggesting that the RFM is highly robust. Our analysis is based on a contraction principle and the theory of monotone dynamical systems. These analysis tools may prove useful in studying other properties of the RFM as well as additional intracellular biological processes.

Chaperones divide yeast proteins into classes of expression level and evolutionary rate

It has long been known that many proteins require folding via molecular chaperones for their function. Although it has become apparent that folding imposes constraints on protein sequence evolution, the effects exerted by different chaperone classes are so far unknown. We have analyzed data of protein interaction with the chaperones in Saccharomycescerevisiae using network methods. The results reveal a distinct community structure within the network that was hitherto undetectable with standard statistical tools. Sixty-four yeast chaperones comprise ten distinct modules that are defined by interaction specificity for their 2,691 interacting proteins. The classes of interacting proteins that are in turn defined by their dedicated chaperone modules are distinguished by various physiochemical protein properties and are characterized by significantly different protein expression levels, codon usage, and amino acid substitution rates. Correlations between substitution rate, codon bias, and gene expression level that have long been known for yeast are apparent at the level of the chaperone-defined modules. This indicates that correlated expression, conservation, and codon bias levels for yeast genes are attributable to previously unrecognized effects of protein folding. Proteome-wide categories of chaperone–substrate specificity uncover novel hubs of functional constraint in protein evolution that are conserved across 20 fungal genomes.

Error prevention and mitigation as forces in the evolution of genes and genomes

Why are short introns rarely a multiple of three nucleotides long? Why do essential genes cluster? Why are genes in operons often lined up in the order in which they are needed in the encoded pathway? In this Opinion article, we argue that these and many other - ostensibly disparate - observations are all pieces of an emerging picture in which multiple aspects of gene anatomy and genome architecture have evolved in response to error-prone gene expression.

Composite effects of gene determinants on the translation speed and density of ribosomes

Background

Translation is a central process of life, and its regulation is crucial for cell growth. In this article, focusing on two model organisms, Escherichia coli and Saccharomyces cerevisiae, we study how three major local features of a gene's coding sequence (its adaptation to the tRNA pool, its amino acid charge, and its mRNA folding energy) affect its translation elongation.

Results

We find that each of these three different features has a non-negligible distinct correlation with the speed of translation elongation. In addition, each of these features might contribute independently to slowing down ribosomal speed at the beginning of genes, which was suggested in previous studies to improve ribosomal allocation and the cost of translation, and to decrease ribosomal jamming. Remarkably, a model of ribosomal translation based on these three basic features highly correlated with the genomic profile of ribosomal density. The robustness to transcription errors in terms of the values of these features is higher at the beginnings of genes, suggesting that this region is important for translation.

Conclusions

The reported results support the conjecture that translation elongation speed is affected by the three coding sequence determinants mentioned above, and not only by adaptation to the tRNA pool; thus, evolution shapes all these determinants along the coding sequences and across genes to improve the organism's translation efficiency.

Preventing dangerous nonsense: selection for robustness to transcriptional error in human genes

Nonsense Mediated Decay (NMD) degrades transcripts that contain a premature STOP codon resulting from mistranscription or missplicing. However NMD's surveillance of gene expression varies in efficiency both among and within human genes. Previous work has shown that the intron content of human genes is influenced by missplicing events invisible to NMD. Given the high rate of transcriptional errors in eukaryotes, we hypothesized that natural selection has promoted a dual strategy of “prevention and cure” to alleviate the problem of nonsense transcriptional errors. A prediction of this hypothesis is that NMD's inefficiency should leave a signature of “transcriptional robustness” in human gene sequences that reduces the frequency of nonsense transcriptional errors. For human genes we determined the usage of “fragile” codons, prone to mistranscription into STOP codons, relative to the usage of “robust” codons that do not generate nonsense errors. We observe that single-exon genes have evolved to become robust to mistranscription, because they show a significant tendency to avoid fragile codons relative to robust codons when compared to multi-exon genes. A similar depletion is evident in last exons of multi-exon genes. Histone genes are particularly depleted of fragile codons and thus highly robust to transcriptional errors. Finally, the protein products of single-exon genes show a strong tendency to avoid those amino acids that can only be encoded using fragile codons. Each of these observations can be attributed to NMD deficiency. Thus, in the human genome, wherever the “cure” for nonsense (i.e. NMD) is inefficient, there is increased reliance on the strategy of nonsense “prevention” (i.e. transcriptional robustness). This study shows that human genes are exposed to the deleterious influence of transcriptional errors. Moreover, it suggests that gene expression errors are an underestimated phenomenon, in molecular evolution in general and in selection for genomic robustness in particular.

In biological systems mistakes are made constantly because the cellular machinery is complex and error-prone. Mistakes are made in copying DNA for transmission to offspring (“genetic mutations”) but are much more frequent in the day-to-day task of gene expression. Genetic mutations are heritable and therefore have long been the almost exclusive focus of evolutionary genetics research. In contrast, gene expression errors are not inherited and have tended to be disregarded in evolutionary studies. Here we show how human genes have evolved a mechanism to reduce the occurrence of a specific type of gene expression error—transcriptional errors that create premature STOP codons (so-called “nonsense errors”). Nonsense errors are potentially highly toxic for the cell, so natural selection has evolved a strategy called Nonsense Mediated Decay (NMD) to “cure” such errors. However this cure is inefficient. Here we describe how a preventative strategy of “transcriptional robustness” has evolved to decrease the frequency of nonsense errors. Moreover, these “prevention and cure” strategies are used interchangeably—the most transcriptionally robust genes are those for which NMD is most inefficient. Our work implies that gene expression errors play an important role as supporting actors to genetic mutations in molecular evolution, particularly in the evolution of robustness.

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.

Determinants of translation efficiency and accuracy

A given protein sequence can be encoded by an astronomical number of alternative nucleotide sequences. Recent research has revealed that this flexibility provides evolution with multiple ways to tune the efficiency and fidelity of protein translation and folding.

Proper functioning of biological cells requires that the process of protein expression be carried out with high efficiency and fidelity. Given an amino-acid sequence of a protein, multiple degrees of freedom still remain that may allow evolution to tune efficiency and fidelity for each gene under various conditions and cell types. Particularly, the redundancy of the genetic code allows the choice between alternative codons for the same amino acid, which, although ‘synonymous,' may exert dramatic effects on the process of translation. Here we review modern developments in genomics and systems biology that have revolutionized our understanding of the multiple means by which translation is regulated. We suggest new means to model the process of translation in a richer framework that will incorporate information about gene sequences, the tRNA pool of the organism and the thermodynamic stability of the mRNA transcripts. A practical demonstration of a better understanding of the process would be a more accurate prediction of the proteome, given the transcriptome at a diversity of biological conditions.

Simultaneous analysis of foodborne pathogenic bacteria by an oligonucleotide microarray assay

A rapid and accurate method for simultaneous identification of foodborne infectious pathogens was developed based on oligonucleotide microarray technology. The proposed identification method is based on PCR amplification of the target region of the groEL genes with degenerate primers, followed by the PCR products hybridization with oligonucleotide probes specific for species. The groEL gene amplification products of seventeen species of pathogenic bacteria were hybridized to the oligonucleotide array. Hybridization results were analyzed with digoxigenin-linked enzyme reaction. Results indicated that fifteen species of pathogenic bacteria showed high sensitivity and specificity for the oligonucleotide array, while two other species gave cross-reaction with the E. coli. Our results suggested that microarray analysis of foodborne infectious pathogens might be very useful for simultaneous identification of bacterial pathogens. The oligonucleotide array can also be applied to samples collected in clinical settings of foodborne infections. The superiority of oligonucleotide array over other tests lies on its rapidity, accuracy and efficiency in the diagnosis, treatment and control of foodborne infections.

Functional analysis of archaeal MBF1 by complementation studies in yeast

Background

Multiprotein-bridging factor 1 (MBF1) is a transcriptional co-activator that bridges a sequence-specific activator (basic-leucine zipper (bZIP) like proteins (e.g. Gcn4 in yeast) or steroid/nuclear-hormone receptor family (e.g. FTZ-F1 in insect)) and the TATA-box binding protein (TBP) in Eukaryotes. MBF1 is absent in Bacteria, but is well- conserved in Eukaryotes and Archaea and harbors a C-terminal Cro-like Helix Turn Helix (HTH) domain, which is the only highly conserved, classical HTH domain that is vertically inherited in all Eukaryotes and Archaea. The main structural difference between archaeal MBF1 (aMBF1) and eukaryotic MBF1 is the presence of a Zn ribbon motif in aMBF1. In addition MBF1 interacting activators are absent in the archaeal domain. To study the function and therefore the evolutionary conservation of MBF1 and its single domains complementation studies in yeast (mbf1Δ) as well as domain swap experiments between aMBF1 and yMbf1 were performed.

Results

In contrast to previous reports for eukaryotic MBF1 (i.e. Arabidopsis thaliana, insect and human) the two archaeal MBF1 orthologs, TMBF1 from the hyperthermophile Thermoproteus tenax and MMBF1 from the mesophile Methanosarcina mazei were not functional for complementation of an Saccharomyces cerevisiae mutant lacking Mbf1 (mbf1Δ). Of twelve chimeric proteins representing different combinations of the N-terminal, core domain, and the C-terminal extension from yeast and aMBF1, only the chimeric MBF1 comprising the yeast N-terminal and core domain fused to the archaeal C-terminal part was able to restore full wild-type activity of MBF1.

However, as reported previously for Bombyx mori, the C-terminal part of yeast Mbf1 was shown to be not essential for function. In addition phylogenetic analyses revealed a common distribution of MBF1 in all Archaea with available genome sequence, except of two of the three Thaumarchaeota; Cenarchaeum symbiosum A and Nitrosopumilus maritimus SCM1.

Conclusions

The absence of MBF1-interacting activators in the archaeal domain, the presence of a Zn ribbon motif in the divergent N-terminal domain of aMBF1 and the complementation experiments using archaeal- yeast chimeric proteins presented here suggests that archaeal MBF1 is not able to functionally interact with the transcription machinery and/or Gcn4 of S. cerevisiae. Based on modeling and structural prediction it is tempting to speculate that aMBF1 might act as a single regulator or non-essential transcription factor, which directly interacts with DNA via the positive charged linker or the basal transcription machinery via its Zn ribbon motif and the HTH domain. However, also alternative functions in ribosome biosynthesis and/or functionality have been discussed and therefore further experiments are required to unravel the function of MBF1 in Archaea.

Reviewers

This article was reviewed by William Martin, Patrick Forterre, John van der Oost and Fabian Blombach (nominated by Eugene V Koonin (United States)). For the full reviews, please go to the Reviewer's Reports section.

A gene's ability to buffer variation is predicted by its fitness contribution and genetic interactions

Background

Many single-gene knockouts result in increased phenotypic (e.g., morphological) variability among the mutant's offspring. This has been interpreted as an intrinsic ability of genes to buffer genetic and environmental variation. A phenotypic capacitor is a gene that appears to mask phenotypic variation: when knocked out, the offspring shows more variability than the wild type. Theory predicts that this phenotypic potential should be correlated with a gene's knockout fitness and its number of negative genetic interactions. Based on experimentally measured phenotypic capacity, it was suggested that knockout fitness was unimportant, but that phenotypic capacitors tend to be hubs in genetic and physical interaction networks.

Methodology/Principal Findings

We re-analyse the available experimental data in a combined model, which includes knockout fitness and network parameters as well as expression level and protein length as predictors of phenotypic potential. Contrary to previous conclusions, we find that the strongest predictor is in fact haploid knockout fitness (responsible for 9% of the variation in phenotypic potential), with an additional contribution from the genetic interaction network (5%); once these two factors are taken into account, protein-protein interactions do not make any additional contribution to the variation in phenotypic potential.

Conclusions/Significance

We conclude that phenotypic potential is not a mysterious “emergent” property of cellular networks. Instead, it is very simply determined by the overall fitness reduction of the organism (which in its compromised state can no longer compensate for multiple factors that contribute to phenotypic variation), and by the number (and presumably nature) of genetic interactions of the knocked-out gene. In this light, Hsp90, the prototypical phenotypic capacitor, may not be representative: typical phenotypic capacitors are not direct “buffers” of variation, but are simply genes encoding central cellular functions.

Impact of translational error-induced and error-free misfolding on the rate of protein evolution

Theoretical calculations suggest that, in addition to translational error-induced protein misfolding, a non-negligible fraction of misfolded proteins are error free.

We propose that the anticorrelation between the expression level of a protein and its rate of sequence evolution be explained by an overarching protein-misfolding-avoidance hypothesis that includes selection against both error-induced and error-free protein misfolding, and verify this model by a molecular-level evolutionary simulation.

We provide strong empirical evidence for the protein-misfolding-avoidance hypothesis, including a positive correlation between protein expression level and stability, enrichment of misfolding-minimizing codons and amino acids in highly expressed genes, and stronger evolutionary conservation of residues in which nonsynonymous changes are more likely to increase protein misfolding.

The rate of protein sequence evolution has long been of central interest to molecular evolutionists. Different proteins of the same species evolve at vastly different rates, which is commonly explained by a variation in functional constraint among different proteins (Kimura and Ohta, 1974). However, it is unclear how to quantify the functional constraint of a protein from the knowledge of its function. In the past decade, various types of genomic data from model organisms have been examined to look for the determinants of the rate of protein sequence evolution. The most unexpected discovery was a very strong anticorrelation between the expression level and evolutionary rate of a protein (E–R anticorrelation) (Pal et al, 2001). The prevailing explanation of the E–R anticorrelation is the translational robustness hypothesis (Drummond et al, 2005). This hypothesis posits that mistranslation induces protein misfolding, which is toxic to cells (Figure 1). Consequently, highly expressed proteins are under stronger pressures to be translationally robust and thus are more constrained in sequence evolution. However, the impact of the other source of misfolded proteins, translational error-free proteins (Figure 1), has not been evaluated. By theoretical calculation, computer simulation, and empirical data analysis, we examined the role of selection against both error-induced and error-free protein misfolding in creating the E–R correlation.

Our theoretical calculations suggested that a non-negligible fraction of misfolded proteins are error free. We estimated that when a protein is not very stable, on average ∼20% of misfolded molecules are error free. However, when a protein is very stable, this fraction reduces to ∼5%, which is probably a result of natural selection against protein misfolding.

We conducted a molecular-level evolutionary simulation (Figure 2A) using three different schemes: error-induced misfolding only, error-free misfolding only, and both types of misfolding. As expected, results from the first simulation are similar to those from a previous study that considers only error-induced misfolding (Drummond and Wilke, 2008). Interestingly, the second and third simulations can also generate the same patterns, including a positive correlation between the protein expression level and the unfolding energy (ΔG) of the error-free protein (Figure 2B), a negative correlation between the expression level and the fraction of protein molecules that misfold after being mistranslated (Figure 2C), a negative correlation between ΔG and the evolutionary rate (Figure 2D), and a negative correlation between the expression level and the evolutionary rate (i.e., the E–R anticorrelation) (Figure 2E). Furthermore, we found that selection against protein misfolding is more effective in reducing error-free misfolding than error-induced misfolding.

Based on these results, we propose that an overarching protein-misfolding-avoidance hypothesis that includes both sources of misfolding is superior to the prevailing translational robustness hypothesis, which considers only error-induced misfolding. We tested three key predictions of the protein-misfolding-avoidance hypotheses using yeast data. First, we showed that, consistent with our prediction, a positive correlation exists between the protein expression level and stability, which is measured by the unfolding energy or melting temperature. In addition, protein expression level is negatively correlated with protein aggregation propensity. Second, we found that codons minimizing protein misfolding are used more frequently in highly expressed proteins than in lowly expressed ones. Third, we showed that, within the same protein, amino acid residues in which random nonsynonymous mutations are more likely to increase protein misfolding are evolutionarily more conserved.

Together, these results provide unambiguous evidence that avoidance of both error-induced and error-free protein misfolding is a major source of the E–R anticorrelation and that protein stability and mistranslation have important roles in protein evolution.

What determines the rate of protein evolution is a fundamental question in biology. Recent genomic studies revealed a surprisingly strong anticorrelation between the expression level of a protein and its rate of sequence evolution. This observation is currently explained by the translational robustness hypothesis in which the toxicity of translational error-induced protein misfolding selects for higher translational robustness of more abundant proteins, which constrains sequence evolution. However, the impact of error-free protein misfolding has not been evaluated. We estimate that a non-negligible fraction of misfolded proteins are error free and demonstrate by a molecular-level evolutionary simulation that selection against protein misfolding results in a greater reduction of error-free misfolding than error-induced misfolding. Thus, an overarching protein-misfolding-avoidance hypothesis that includes both sources of misfolding is superior to the translational robustness hypothesis. We show that misfolding-minimizing amino acids are preferentially used in highly abundant yeast proteins and that these residues are evolutionarily more conserved than other residues of the same proteins. These findings provide unambiguous support to the role of protein-misfolding-avoidance in determining the rate of protein sequence evolution.

Evolution of molecular error rates and the consequences for evolvability

Making genes into gene products is subject to predictable errors, each with a phenotypic effect that depends on a normally cryptic sequence. Many cryptic sequences have strongly deleterious effects, for example when they cause protein misfolding. Strongly deleterious effects can be avoided globally by avoiding making errors (e.g., via proofreading machinery) or locally by ensuring that each error has a relatively benign effect. The local solution requires powerful selection acting on every cryptic site and so evolves only in large populations. Small populations with less effective selection evolve global solutions. Here we show that for a large range of realistic intermediate population sizes, the evolutionary dynamics are bistable and either solution may result. The local solution facilitates the genetic assimilation of cryptic genetic variation and therefore substantially increases evolvability.

Transition from positive to neutral in mutation fixation along with continuing rising fitness in thermal adaptive evolution

It remains to be determined experimentally whether increasing fitness is related to positive selection, while stationary fitness is related to neutral evolution. Long-term laboratory evolution in Escherichia coli was performed under conditions of thermal stress under defined laboratory conditions. The complete cell growth data showed common continuous fitness recovery to every 2°C or 4°C stepwise temperature upshift, finally resulting in an evolved E. coli strain with an improved upper temperature limit as high as 45.9°C after 523 days of serial transfer, equivalent to 7,560 generations, in minimal medium. Two-phase fitness dynamics, a rapid growth recovery phase followed by a gradual increasing growth phase, was clearly observed at diverse temperatures throughout the entire evolutionary process. Whole-genome sequence analysis revealed the transition from positive to neutral in mutation fixation, accompanied with a considerable escalation of spontaneous substitution rate in the late fitness recovery phase. It suggested that continually increasing fitness not always resulted in the reduction of genetic diversity due to the sequential takeovers by fit mutants, but caused the accumulation of a considerable number of mutations that facilitated the neutral evolution.

The detailed results of a two-year in vitro thermal adaptive evolution experiment are described. A laboratory-evolved E. coli strain with an improved upper temperature limit, as high as 45.9°C, was acquired after 523 days of serial transfer, equivalent to 7,560 generations, in nutrient-limited medium. The complete daily records of cell growth exhibited universal two-phase fitness dynamics, a rapid growth recovery phase followed by a gradual increasing growth phase, throughout the entire evolutionary process. Genome-sequence analysis not only showed considerable escalation of the spontaneous substitution rate, but also revealed the transition from positive to nearly neutral in mutation fixation. Particularly, even with the rising fitness of bacterial cells, neutrality was observed in molecular evolution. These observations suggested that a discrete evolutionary mode occurred in the continuous evolutionary route, linking Darwinian adaptive selection with Kimura's neutral evolution. Such transition from beneficial to neutral path may be adopted as an evolutionary strategy robust to rigorous environmental changes.

Chaperonin-dependent accelerated substitution rates in prokaryotes

Many proteins require the assistance of molecular chaperones in order to fold efficiently. Chaperones are known to mask the effects of mutations that induce misfolding because they can compensate for the deficiency in spontaneous folding. One of the best studied chaperones is the eubacterial GroEL/GroES system. In Escherichia coli, three classes of proteins have been distinguished based on their degree of dependency on GroEL for folding: 1) those that do not require GroEL, 2) those that require GroEL in a temperature-dependent manner, and 3) those that obligately require GroEL for proper folding. The buffering effects of GroEL have so far been observed in experimental regimens, but their effect on genomes during evolution has not been examined. Using 446 sequenced proteobacterial genomes, we have compared the frequency of amino acid replacements among orthologs of 236 proteins corresponding to the three categories of GroEL dependency determined for E. coli. Evolutionary rates are significantly correlated with GroEL dependency upon folding with GroEL dependency class accounting for up to 84% of the variation in amino acid substitution rates. Greater GroEL dependency entails increased evolutionary rates with GroEL obligatory proteins (Class III) evolving on average up to 15% faster than GroEL partially dependent proteins (Class II) and 35% faster than GroEL-independent proteins (Class I). Moreover, GroEL dependency class correlations are strictly conserved throughout all proteobacteria surveyed, as is a significant correlation between folding class and codon bias. The results suggest that during evolution, GroEL-dependent folding increases evolutionary rate by buffering the deleterious effects of misfolding-related mutations.

The effect of chaperonin buffering on protein evolution

Molecular chaperones are highly conserved and ubiquitous proteins that help other proteins in the cell to fold. Pioneering work by Rutherford and Lindquist suggested that the chaperone Hsp90 could buffer (i.e., suppress) phenotypic variation in its client proteins and that alternate periods of buffering and expression of these variants might be important in adaptive evolution. More recently, Tokuriki and Tawfik presented an explicit mechanism for chaperone-dependent evolution, in which the Escherichia coli chaperonin GroEL facilitated the folding of clients that had accumulated structurally destabilizing but neofunctionalizing mutations in the protein core. But how important an evolutionary force is chaperonin-mediated buffering in nature? Here, we address this question by modeling the per-residue evolutionary rate of the crystallized E. coli proteome, evaluating the relative contributions of chaperonin buffering, functional importance, and structural features such as residue contact density. Previous findings suggest an interaction between codon bias and GroEL in limiting the effects of misfolding errors. Our results suggest that the buffering of deleterious mutations by GroEL increases the evolutionary rate of client proteins. We then examine the evolutionary fate of GroEL clients in the Mycoplasmas, a group of bacteria containing the only known organisms that lack chaperonins. We show that GroEL was lost once in the common ancestor of a monophyletic subgroup of Mycoplasmas, and we evaluate the effect of this loss on the subsequent evolution of client proteins, providing evidence that client homologs in 11 Mycoplasma species have lost their obligate dependency on GroEL for folding. Our analyses indicate that individual molecules such as chaperonins can have significant effects on proteome evolution through their modulation of protein folding.

Signatures of protein biophysics in coding sequence evolution

Since the early days of molecular evolution, the conventional wisdom has been that the evolution of protein-coding genes is primarily determined by functional constraints. Yet recent evidence indicates that the evolution of these genes is strongly shaped by the biophysical processes of protein synthesis, protein folding, and specific as well as nonspecific protein-protein interactions. Selection pressures related to these biophysical processes affect primarily the amino-acid sequence of genes, but they also leave their mark on synonymous sites at the nucleotide level. While evidence for specific selection pressures related to protein biophysics is strong, there is currently no unifying framework that integrates the various selection pressures on coding sequences and disentangles their relative importance.