Potential role of phenotypic mutations in the evolution of protein expression and stability

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G·T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R1ρ experiments, we show that dG·rU and dT·rG mismatches in two RNA:DNA hybrids transiently form tautomeric (Genol·T/U $ \mathbin{\lower.3ex\hbox{$\buildrel\textstyle\rightarrow\over {\smash{\leftarrow}\vphantom{_{\vbox to.5ex{\vss}}}}$}}$ G·Tenol/Uenol) and anionic (G·T-/U-) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG·rU- formed with a ten-fold higher propensity relative to dT-·rG and dG·dT- and this could be attributed to the lower pKa (ΔpKa ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G·T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G·U versus G·T, and indicate that anionic Watson-Crick-like G·U- could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

Fitness Effects of Phenotypic Mutations at Proteome-Scale Reveal Optimality of Translation Machinery

Errors in protein translation can lead to non-genetic, phenotypic mutations, including amino acid misincorporations. While phenotypic mutations can increase protein diversity, the systematic characterization of their proteome-wide frequencies and their evolutionary impact has been lacking. Here, we developed a mechanistic model of translation errors to investigate how selection acts on protein populations produced by amino acid misincorporations. We fitted the model to empirical observations of misincorporations obtained from over a hundred mass spectrometry datasets of E. coli and S. cerevisiae. We found that on average 20% to 23% of proteins synthesized in the cell are expected to harbor at least one amino acid misincorporation, and that deleterious misincorporations are less likely to occur. Combining misincorporation probabilities and the estimated fitness effects of amino acid substitutions in a population genetics framework, we found 74% of mistranslation events in E. coli and 94% in S. cerevisiae to be neutral. We further show that the set of available synonymous tRNAs is subject to evolutionary pressure, as the presence of missing tRNAs would increase codon-anticodon cross-reactivity and misincorporation error rates. Overall, we find that the translation machinery is likely optimal in E. coli and S. cerevisiae and that both local solutions at the level of codons and a global solution such as the tRNA pool can mitigate the impact of translation errors. We provide a framework to study the evolutionary impact of codon-specific translation errors and a method for their proteome-wide detection across organisms and conditions.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: Implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G•T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R 1ρ experiments, we show that dG•rU and dT•rG mismatches in two RNA:DNA hybrids transiently form tautomeric (G enol •T/U ⇄G•T enol /U enol ) and anionic (G•T - /U - ) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG•rU - formed with a ten-fold higher propensity relative to dT - •rG and dG•dT - and this could be attributed to the lower pK a (Δ pK a ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G•T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G•U versus G•T, and indicate that anionic Watson-Crick-like G•U - could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

Phenotypic mutations contribute to protein diversity and shape protein evolution

Errors in DNA replication generate genetic mutations, while errors in transcription and translation lead to phenotypic mutations. Phenotypic mutations are orders of magnitude more frequent than genetic ones, yet they are less understood. Here, we review the types of phenotypic mutations, their quantifications, and their role in protein evolution and disease. The diversity generated by phenotypic mutation can facilitate adaptive evolution. Indeed, phenotypic mutations, such as ribosomal frameshift and stop codon readthrough, sometimes serve to regulate protein expression and function. Phenotypic mutations have often been linked to fitness decrease and diseases. Thus, understanding the protein heterogeneity and phenotypic diversity caused by phenotypic mutations will advance our understanding of protein evolution and have implications on human health and diseases.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

Resurrected Rubisco suggests uniform carbon isotope signatures over geologic time

The earliest geochemical indicators of microbes-and the enzymes that powered them-extend back ∼3.8 Ga on Earth. Paleobiologists often attempt to understand these indicators by assuming that the behaviors of extant microbes and enzymes are uniform with those of their predecessors. This consistency in behavior seems at odds with our understanding of the inherent variability of living systems. Here, we examine whether a uniformitarian assumption for an enzyme thought to generate carbon isotope indicators of biological activity, RuBisCO, can be corroborated by independently studying the history of changes recorded within RuBisCO's genetic sequences. We resurrected a Precambrian-age RuBisCO by engineering its ancient DNA inside a cyanobacterium genome and measured the engineered organism's fitness and carbon-isotope-discrimination profile. Results indicate that Precambrian uniformitarian assumptions may be warranted but with important caveats. Experimental studies illuminating early innovations are crucial to explore the molecular foundations of life's earliest traces.

Mistranslation reduces mutation load in evolving proteins through negative epistasis with DNA mutations

Translational errors during protein synthesis cause phenotypic mutations that are several orders of magnitude more frequent than DNA mutations. Such phenotypic mutations may affect adaptive evolution through their interactions with DNA mutations. To study how mistranslation may affect the adaptive evolution of evolving proteins, we evolved populations of green fluorescent protein (GFP) in either high-mistranslation or low-mistranslation Escherichia coli hosts. In both hosts, we first evolved GFP under purifying selection for the ancestral phenotype green fluorescence, and then under directional selection toward the new phenotype yellow fluorescence. High-mistranslation populations evolved modestly higher yellow fluorescence during each generation of evolution than low-mistranslation populations. We demonstrate by high-throughput sequencing that elevated mistranslation reduced the accumulation of deleterious DNA mutations under both purifying and directional selection. It did so by amplifying the fitness effects of deleterious DNA mutations through negative epistasis with phenotypic mutations. In contrast, mistranslation did not affect the incidence of beneficial mutations. Our findings show that phenotypic mutations interact epistatically with DNA mutations. By reducing a population’s mutation load, mistranslation can affect an important determinant of evolvability.

Improve it or lose it: Evolvability cost of competition for expression

Expression level is known to be a strong determinant of a protein's rate of evolution. But the converse can also be true: evolutionary dynamics can affect expression levels of proteins. Having implications in both directions fosters the possibility of an "improve it or lose it" feedback loop, where higher expressed systems are more likely to improve and be expressed even higher, while those that are expressed less are eventually lost to drift. Using a minimal model to study this in the context of a changing environment, we demonstrate that one unexpected consequence of such a feedback loop is that a slow switch to a new environment can allow genotypes to reach higher fitness sooner than a direct exposure to it.

Evidence for a role of phenotypic mutations in virus adaptation

Viruses interact extensively with the host molecular machinery, but the underlying mechanisms are poorly understood. Bacteriophage T7 recruits the small protein thioredoxin of the Escherichia coli host as an essential processivity factor for the viral DNA polymerase. We challenged the phage to propagate in a host in which thioredoxin had been extensively modified to hamper its recruitment. The virus adapted to the engineered host without losing the capability to propagate in the original host, but no genetic mutations were fixed in the thioredoxin binding domain of the viral DNA polymerase. Virus adaptation correlated with mutations in the viral RNA polymerase, supporting that promiscuous thioredoxin recruitment was enabled by phenotypic mutations caused by transcription errors. These results point to a mechanism of virus adaptation that may play a role in cross-species transmission. We propose that phenotypic mutations may generally contribute to the capability of viruses to evade antiviral strategies.

•

Phage adapts to a host modified to hinder the essential recruitment of a host protein

•

No genetic mutations are fixed at the engineered virus-host interaction interface

•

Adaptation is likely linked to phenotypic mutations caused by transcription errors

•

Sub-genomic RNAs may enable this kind of adaptation mechanism in coronaviruses

The Relationship between the Misfolding Avoidance Hypothesis and Protein Evolutionary Rates in the Light of Empirical Evidence

For more than a decade, the misfolding avoidance hypothesis (MAH) and related theories have dominated evolutionary discussions aimed at explaining the variance of the molecular clock across cellular proteins. In this study, we use various experimental data to further investigate the consistency of the MAH predictions with empirical evidence. We also critically discuss experimental results that motivated the MAH development and that are often viewed as evidence of its major contribution to the variability of protein evolutionary rates. We demonstrate, in Escherichia coli and Homo sapiens, the lack of a substantial negative correlation between protein evolutionary rates and Gibbs free energies of unfolding, a direct measure of protein stability. We then analyze multiple new genome-scale data sets characterizing protein aggregation and interaction propensities, the properties that are likely optimized in evolution to alleviate deleterious effects associated with toxic protein misfolding and misinteractions. Our results demonstrate that the propensity of proteins to aggregate, the fraction of charged amino acids, and protein stickiness do correlate with protein abundances. Nevertheless, across multiple organisms and various data sets we do not observe substantial correlations between proteins’ aggregation- and stability-related properties and evolutionary rates. Therefore, diverse empirical data support the conclusion that the MAH and similar hypotheses do not play a major role in mediating a strong negative correlation between protein expression and the molecular clock, and thus in explaining the variability of evolutionary rates across cellular proteins.

The Boggarts of biology: how non-genetic changes influence the genotype

The notion that there is a one-one mapping from genotype to phenotype was overturned a long time ago. Along with genotype and environment, 'non-genetic changes' orchestrated by altered RNA and protein molecules also guide the development of phenotype. The idea that there is a route through which changes in phenotype can lead to changes in genotype impinges on several phenomena of molecular, developmental, evolutionary and applied interest. Phenotypic changes that do not alter the underlying DNA sequence have been studied across model systems (eg: DNA and histone modifications, RNA editing, prion formation) and are known to play an important role in short-term adaptation. However, because of their transient nature and unstable inheritance, the role of such changes in long-term evolution has remained controversial. I classify and review three ways in which non-genetic changes can influence genotype and impact cellular fitness across generations, with an emphasis on the enticing idea that they may act as stepping stones for genetic adaptation. I focus on work from microbial systems and attempt to highlight recent experiments and models that bear on this idea. Overall, I review evidence which suggests that non-genetic changes can impact phenotype via their influence on the genotype, and thus play a role in evolutionary change.

Structures of hyperstable ancestral haloalkane dehalogenases show restricted conformational dynamics

Ancestral sequence reconstruction is a powerful method for inferring ancestors of modern enzymes and for studying structure-function relationships of enzymes. We have previously applied this approach to haloalkane dehalogenases (HLDs) from the subfamily HLD-II and obtained thermodynamically highly stabilized enzymes (ΔT m up to 24 °C), showing improved catalytic properties. Here we combined crystallographic structural analysis and computational molecular dynamics simulations to gain insight into the mechanisms by which ancestral HLDs became more robust enzymes with novel catalytic properties. Reconstructed ancestors exhibited similar structure topology as their descendants with the exception of a few loop deviations. Strikingly, molecular dynamics simulations revealed restricted conformational dynamics of ancestral enzymes, which prefer a single state, in contrast to modern enzymes adopting two different conformational states. The restricted dynamics can potentially be linked to their exceptional stabilization. The study provides molecular insights into protein stabilization due to ancestral sequence reconstruction, which is becoming a widely used approach for obtaining robust protein catalysts.

Protein Structural Information and Evolutionary Landscape by In Vitro Evolution

Protein structure is tightly intertwined with function according to the laws of evolution. Understanding how structure determines function has been the aim of structural biology for decades. Here, we have wondered instead whether it is possible to exploit the function for which a protein was evolutionary selected to gain information on protein structure and on the landscape explored during the early stages of molecular and natural evolution. To answer to this question, we developed a new methodology, which we named CAMELS (Coupling Analysis by Molecular Evolution Library Sequencing), that is able to obtain the in vitro evolution of a protein from an artificial selection based on function. We were able to observe with CAMELS many features of the TEM-1 beta-lactamase local fold exclusively by generating and sequencing large libraries of mutational variants. We demonstrated that we can, whenever a functional phenotypic selection of a protein is available, sketch the structural and evolutionary landscape of a protein without utilizing purified proteins, collecting physical measurements, or relying on the pool of natural protein variants.

How evolution shapes enzyme selectivity - lessons from aminoacyl-tRNA synthetases and other amino acid utilizing enzymes

Aminoacyl-tRNA synthetases (AARSs) charge tRNA with their cognate amino acids. Many other enzymes use amino acids as substrates, yet discrimination against noncognate amino acids that threaten the accuracy of protein translation is a hallmark of AARSs. Comparing AARSs to these other enzymes allowed us to recognize patterns in molecular recognition and strategies used by evolution for exercising selectivity. Overall, AARSs are 2-3 orders of magnitude more selective than most other amino acid utilizing enzymes. AARSs also reveal the physicochemical limits of molecular discrimination. For example, amino acids smaller by a single methyl moiety present a discrimination ceiling of ~200, while larger ones can be discriminated by up to 10⁵ -fold. In contrast, substrates larger by a hydroxyl group challenge AARS selectivity, due to promiscuous H-bonding with polar active site groups. This 'hydroxyl paradox' is resolved by editing. Indeed, when the physicochemical discrimination limits are reached, post-transfer editing - hydrolysis of tRNAs charged with noncognate amino acids, evolved. The editing site often selectively recognizes the edited noncognate substrate using the very same feature that the synthetic site could not efficiently discriminate against. Finally, the comparison to other enzymes also reveals that the selectivity of AARSs is an explicitly evolved trait, showing some clear examples of how selection acted not only to optimize catalytic efficiency with the target substrate, but also to abolish activity with noncognate threat substrates ('negative selection').

Robustness by intrinsically disordered C-termini and translational readthrough

During protein synthesis genetic instructions are passed from DNA via mRNA to the ribosome to assemble a protein chain. Occasionally, stop codons in the mRNA are bypassed and translation continues into the untranslated region (3′-UTR). This process, called translational readthrough (TR), yields a protein chain that becomes longer than would be predicted from the DNA sequence alone. Protein sequences vary in propensity for translational errors, which may yield evolutionary constraints by limiting evolutionary paths. Here we investigated TR in Saccharomyces cerevisiae by analysing ribosome profiling data. We clustered proteins as either prone or non-prone to TR, and conducted comparative analyses. We find that a relatively high frequency (5%) of genes undergo TR, including ribosomal subunit proteins. Our main finding is that proteins undergoing TR are highly expressed and have a higher proportion of intrinsically disordered C-termini. We suggest that highly expressed proteins may compensate for the deleterious effects of TR by having intrinsically disordered C-termini, which may provide conformational flexibility but without distorting native function. Moreover, we discuss whether minimizing deleterious effects of TR is also enabling exploration of the phenotypic landscape of protein isoforms.

The causes of evolvability and their evolution

Evolvability is the ability of a biological system to produce phenotypic variation that is both heritable and adaptive. It has long been the subject of anecdotal observations and theoretical work. In recent years, however, the molecular causes of evolvability have been an increasing focus of experimental work. Here, we review recent experimental progress in areas as different as the evolution of drug resistance in cancer cells and the rewiring of transcriptional regulation circuits in vertebrates. This research reveals the importance of three major themes: multiple genetic and non-genetic mechanisms to generate phenotypic diversity, robustness in genetic systems, and adaptive landscape topography. We also discuss the mounting evidence that evolvability can evolve and the question of whether it evolves adaptively.

What drives phenotypic divergence among coral clonemates of Acropora palmata?

Evolutionary rescue of populations depends on their ability to produce phenotypic variation that is heritable and adaptive. DNA mutations are the best understood mechanisms to create phenotypic variation, but other, less well-studied mechanisms exist. Marine benthic foundation species provide opportunities to study these mechanisms because many are dominated by isogenic stands produced through asexual reproduction. For example, Caribbean acroporid corals are long lived and reproduce asexually via breakage of branches. Fragmentation is often the dominant mode of local population maintenance. Thus, large genets with many ramets (colonies) are common. Here, we observed phenotypic variation in stress responses within genets following the coral bleaching events in 2014 and 2015 caused by high water temperatures. This was not due to genetic variation in their symbiotic dinoflagellates (Symbiodinium "fitti") because each genet of this coral species typically harbours a single strain of S. "fitti". Characterization of the microbiome via 16S tag sequencing correlated the abundance of only two microbiome members (Tepidiphilus, Endozoicomonas) with a bleaching response. Epigenetic changes were significantly correlated with the host's genetic background, the location of the sampled polyps within the colonies (e.g., branch vs. base of colony), and differences in the colonies' condition during the bleaching event. We conclude that long-term microenvironmental differences led to changes in the way the ramets methylated their genomes, contributing to the differential bleaching response. However, most of the variation in differential bleaching response among clonemates of Acropora palmata remains unexplained. This research provides novel data and hypotheses to help understand intragenet variability in stress phenotypes of sessile marine species.

Defining Division of Labor in Microbial Communities

In order to survive and reproduce, organisms must perform a multitude of tasks. However, trade-offs limit their ability to allocate energy and resources to all of these different processes. One strategy to solve this problem is to specialize in some traits and team up with other organisms that can help by providing additional, complementary functions. By reciprocally exchanging metabolites and/or services in this way, both parties benefit from the interaction. This phenomenon, which has been termed functional specialization or division of labor, is very common in nature and exists on all levels of biological organization. Also, microorganisms have evolved different types of synergistic interactions. However, very often, it remains unclear whether or not a given example represents a true case of division of labor. Here we aim at filling this gap by providing a list of criteria that clearly define division of labor in microbial communities. Furthermore, we propose a set of diagnostic experiments to verify whether a given interaction fulfills these conditions. In contrast to the common use of the term, our analysis reveals that both intraspecific and interspecific interactions meet the criteria defining division of labor. Moreover, our analysis identified non-cooperators of intraspecific public goods interactions as growth specialists that divide labor with conspecific producers, rather than being social parasites. By providing a conceptual toolkit, our work will help to unambiguously identify cases of division of labor and stimulate more detailed investigations of this important and widespread type of inter-microbial interaction.

Population Genetics Based Phylogenetics Under Stabilizing Selection for an Optimal Amino Acid Sequence: A Nested Modeling Approach

We present a new phylogenetic approach, selection on amino acids and codons (SelAC), whose substitution rates are based on a nested model linking protein expression to population genetics. Unlike simpler codon models that assume a single substitution matrix for all sites, our model more realistically represents the evolution of protein-coding DNA under the assumption of consistent, stabilizing selection using a cost-benefit approach. This cost-benefit approach allows us to generate a set of 20 optimal amino acid-specific matrix families using just a handful of parameters and naturally links the strength of stabilizing selection to protein synthesis levels, which we can estimate. Using a yeast data set of 100 orthologs for 6 taxa, we find SelAC fits the data much better than popular models by 104-105 Akike information criterion units adjusted for small sample bias. Our results also indicated that nested, mechanistic models better predict observed data patterns highlighting the improvement in biological realism in amino acid sequence evolution that our model provides. Additional parameters estimated by SelAC indicate that a large amount of nonphylogenetic, but biologically meaningful, information can be inferred from existing data. For example, SelAC prediction of gene-specific protein synthesis rates correlates well with both empirical (r=0.33-0.48) and other theoretical predictions (r=0.45-0.64) for multiple yeast species. SelAC also provides estimates of the optimal amino acid at each site. Finally, because SelAC is a nested approach based on clearly stated biological assumptions, future modifications, such as including shifts in the optimal amino acid sequence within or across lineages, are possible.

Identification of Novel Neisseria gonorrhoeae Lineages Harboring Resistance Plasmids in Coastal Kenya

Background

Africa has the highest incidence of gonorrhea in the world. However, little is known about gonococcal populations in this continent or mechanisms of antimicrobial resistance (AMR).

Methods

Whole-genome sequence data were analyzed from 103 Neisseria gonorrhoeae isolates from 73 patients, mainly men who have sex with men, from coastal Kenya. We annotated loci, defined the core genome, defined mechanisms of AMR, and performed phylogenetic analysis. For patients with multiple episodes of gonorrhea, we determined whether infections occurred with related strains.

Results

We identified 3 clusters of isolates that are phylogenetically distinct from isolates found elsewhere. Plasmids were virtually ubiquitous: pTetM and pblaTEM were found in 97%, and 55% of isolates, respectively. This was associated with high doxycycline use for undiagnosed sexually transmitted infections. Twenty-three percent of multiple episodes of gonorrhea in the same individual were caused by a related strain, suggesting inadequate treatment or reinfection.

Conclusions

The prevalence of plasmid-mediated AMR in Kenyan gonococci contrasts with that in wealthy countries, where AMR is largely chromosomally mediated. Antimicrobials have a profound effect on the maintenance of lineages harboring plasmids. Doxycycline can select for tetracycline and penicillin resistance, through plasmid cooperation. Understanding the mechanisms of AMR in high-risk groups is required to inform treatment strategies.

Errors during Gene Expression: Single-Cell Heterogeneity, Stress Resistance, and Microbe-Host Interactions

Gene expression has been considered a highly accurate process, and deviation from such fidelity has been shown previously to be detrimental for the cell. More recently, increasing evidence has supported the notion that the accuracy of gene expression is indeed flexibly variable. The levels of errors during gene expression differ from condition to condition and even from cell to cell within genetically identical populations grown under the same conditions. The different levels of errors resulting from inaccurate gene expression are now known to play key roles in regulating microbial stress responses and host interactions. This minireview summarizes the recent development in understanding the level, regulation, and physiological impact of errors during gene expression.

Gene expression has been considered a highly accurate process, and deviation from such fidelity has been shown previously to be detrimental for the cell. More recently, increasing evidence has supported the notion that the accuracy of gene expression is indeed flexibly variable.

Conformational dynamics and enzyme evolution

Enzymes are dynamic entities, and their dynamic properties are clearly linked to their biological function. It follows that dynamics ought to play an essential role in enzyme evolution. Indeed, a link between conformational diversity and the emergence of new enzyme functionalities has been recognized for many years. However, it is only recently that state-of-the-art computational and experimental approaches are revealing the crucial molecular details of this link. Specifically, evolutionary trajectories leading to functional optimization for a given host environment or to the emergence of a new function typically involve enriching catalytically competent conformations and/or the freezing out of non-competent conformations of an enzyme. In some cases, these evolutionary changes are achieved through distant mutations that shift the protein ensemble towards productive conformations. Multifunctional intermediates in evolutionary trajectories are probably multi-conformational, i.e. able to switch between different overall conformations, each competent for a given function. Conformational diversity can assist the emergence of a completely new active site through a single mutation by facilitating transition-state binding. We propose that this mechanism may have played a role in the emergence of enzymes at the primordial, progenote stage, where it was plausibly promoted by high environmental temperatures and the possibility of additional phenotypic mutations.

Characterization of Reconstructed Ancestral Proteins Suggests a Change in Temperature of the Ancient Biosphere

Understanding the evolution of ancestral life, and especially the ability of some organisms to flourish in the variable environments experienced in Earth’s early biosphere, requires knowledge of the characteristics and the environment of these ancestral organisms. Information about early life and environmental conditions has been obtained from fossil records and geological surveys. Recent advances in phylogenetic analysis, and an increasing number of protein sequences available in public databases, have made it possible to infer ancestral protein sequences possessed by ancient organisms. However, the in silico studies that assess the ancestral base content of ribosomal RNAs, the frequency of each amino acid in ancestral proteins, and estimate the environmental temperatures of ancient organisms, show conflicting results. The characterization of ancestral proteins reconstructed in vitro suggests that ancient organisms had very thermally stable proteins, and therefore were thermophilic or hyperthermophilic. Experimental data supports the idea that only thermophilic ancestors survived the catastrophic increase in temperature of the biosphere that was likely associated with meteorite impacts during the early history of Earth. In addition, by expanding the timescale and including more ancestral proteins for reconstruction, it appears as though the Earth’s surface temperature gradually decreased over time, from Archean to present.

Time-programmable drug dosing allows the manipulation, suppression and reversal of antibiotic drug resistance in vitro

Multi-drug strategies have been attempted to prolong the efficacy of existing antibiotics, but with limited success. Here we show that the evolution of multi-drug-resistant Escherichia coli can be manipulated in vitro by administering pairs of antibiotics and switching between them in ON/OFF manner. Using a multiplexed cell culture system, we find that switching between certain combinations of antibiotics completely suppresses the development of resistance to one of the antibiotics. Using this data, we develop a simple deterministic model, which allows us to predict the fate of multi-drug evolution in this system. Furthermore, we are able to reverse established drug resistance based on the model prediction by modulating antibiotic selection stresses. Our results support the idea that the development of antibiotic resistance may be potentially controlled via continuous switching of drugs.

Variation in Mutational Robustness between Different Proteins and the Predictability of Fitness Effects

Random mutations in genes from disparate protein classes may have different distributions of fitness effects (DFEs) depending on different structural, functional, and evolutionary constraints. We measured the fitness effects of 156 single mutations in the genes encoding AraC (transcription factor), AraD (enzyme), and AraE (transporter) used for bacterial growth on l-arabinose. Despite their different molecular functions these genes all had bimodal DFEs with most mutations either being neutral or strongly deleterious, providing a general expectation for the DFE. This contrasts with the unimodal DFEs previously obtained for ribosomal protein genes where most mutations were slightly deleterious. Based on theoretical considerations, we suggest that the 33-fold higher average mutational robustness of ribosomal proteins is due to stronger selection for reduced costs of translational and transcriptional errors. Whereas the large majority of synonymous mutations were deleterious for ribosomal proteins genes, no fitness effects could be detected for the AraCDE genes. Four mutations in AraC and AraE increased fitness, suggesting that slightly advantageous mutations make up a significant fraction of the DFE, but that they often escape detection due to the limited sensitivity of commonly used fitness assays. We show that the fitness effects of amino acid substitutions can be predicted based on evolutionary conservation, but those weakly deleterious mutations are less reliably detected. This suggests that large-effect mutations and the fraction of highly deleterious mutations can be computationally predicted, but that experiments are required to characterize the DFE close to neutrality, where many mutations ultimately fixed in a population will occur.

Pharmacodynamics of Ceftolozane Combined with Tazobactam against Enterobacteriaceae in a Neutropenic Mouse Thigh Model

Ceftolozane is a new broad-spectrum cephalosporin and is combined with tazobactam to broaden the activity of ceftolozane against strains producing extended-spectrum beta-lactamases (ESBLs). We determined the pharmacodynamics (PD) of the combination in the neutropenic mouse thigh model to determine the optimal exposure of tazobactam. Treatment of CD-1 neutropenic mice was started 2 h after infection with ceftolozane every 2 h (q2h) alone or in combination with tazobactam at different dosing frequencies for 24 h, and the number of CFU in the thighs was determined before and after treatment. The maximum effect model was fit to the dose-response and the pharmacokinetic/PD index (PDI)-response to determine the PDI values for ceftolozane alone and ceftolozane in combination with tazobactam resulting in a static effect and a 1-log kill. The effect of tazobactam was dependent on the percentage of time that the free drug concentration remained above the concentration threshold (percent [Formula: see text]), whereby dosing q2h was more efficacious than dosing every 8 h (q8h), reducing the tazobactam daily dose by a factor 6.9 to 59.0 (n = 3 strains) to obtain a static effect. Using R² as an indicator of the best fit of the percent [Formula: see text]-response relationships, the concentration threshold best correlating with the response varied from 0.5 to 2 mg/liter, depending on the strain. A similar result was obtained when the q2h and q8h regimens were analyzed. For all isolates tested, the mean [Formula: see text] for 0.5 mg/liter tazobactam was 28.2% (range, 17.5 to 45.8%) and 44.4% (range, 26.6 to 54.7%) for a static effect and a 1-log kill, respectively, at ceftolozane exposures that produced a ceftolozane concentration of 4 mg/liter (a concentration greater than the MIC) for 33.9 to 63.3% of a 24-h period under steady-state pharmacokinetic conditions. The main PDI that correlated with the effect of tazobactam was the [Formula: see text] achieved with a C_T of 0.5 mg/liter tazobactam.

On the Potential Origins of the High Stability of Reconstructed Ancestral Proteins

Ancestral reconstruction provides instrumental insights regarding the biochemical and biophysical characteristics of past proteins. A striking observation relates to the remarkably high thermostability of reconstructed ancestors. The latter has been linked to high environmental temperatures in the Precambrian era, the era relating to most reconstructed proteins. We found that inferred ancestors of the serum paraoxonase (PON) enzyme family, including the mammalian ancestor, exhibit dramatically increased thermostabilities compared with the extant, human enzyme (up to 30 °C higher melting temperature). However, the environmental temperature at the time of emergence of mammals is presumed to be similar to the present one. Additionally, the mammalian PON ancestor has superior folding properties (kinetic stability)-unlike the extant mammalian PONs, it expresses in E. coli in a soluble and functional form, and at a high yield. We discuss two potential origins of this unexpectedly high stability. First, ancestral stability may be overestimated by a "consensus effect," whereby replacing amino acids that are rare in contemporary sequences with the amino acid most common in the family increases protein stability. Comparison to other reconstructed ancestors indicates that the consensus effect may bias some but not all reconstructions. Second, we note that high stability may relate to factors other than high environmental temperature such as oxidative stress or high radiation levels. Foremost, intrinsic factors such as high rates of genetic mutations and/or of transcriptional and translational errors, and less efficient protein quality control systems, may underlie the high kinetic and thermodynamic stability of past proteins.

Mistranslation drives the evolution of robustness in TEM-1 β-lactamase

How biological systems such as proteins achieve robustness to ubiquitous perturbations is a fundamental biological question. Such perturbations include errors that introduce phenotypic mutations into nascent proteins during the translation of mRNA. These errors are remarkably frequent. They are also costly, because they reduce protein stability and help create toxic misfolded proteins. Adaptive evolution might reduce these costs of protein mistranslation by two principal mechanisms. The first increases the accuracy of translation via synonymous "high fidelity" codons at especially sensitive sites. The second increases the robustness of proteins to phenotypic errors via amino acids that increase protein stability. To study how these mechanisms are exploited by populations evolving in the laboratory, we evolved the antibiotic resistance gene TEM-1 in Escherichia coli hosts with either normal or high rates of mistranslation. We analyzed TEM-1 populations that evolved under relaxed and stringent selection for antibiotic resistance by single molecule real-time sequencing. Under relaxed selection, mistranslating populations reduce mistranslation costs by reducing TEM-1 expression. Under stringent selection, they efficiently purge destabilizing amino acid changes. More importantly, they accumulate stabilizing amino acid changes rather than synonymous changes that increase translational accuracy. In the large populations we study, and on short evolutionary timescales, the path of least resistance in TEM-1 evolution consists of reducing the consequences of translation errors rather than the errors themselves.

Transcription yield of fully 2'-modified RNA can be increased by the addition of thermostabilizing mutations to T7 RNA polymerase mutants

On average, mutations are deleterious to proteins. Mutations conferring new function to a protein often come at the expense of protein folding or stability, reducing overall activity. Over the years, a panel of T7 RNA polymerases have been designed or evolved to accept nucleotides with modified ribose moieties. These modified RNAs have proven useful, especially in vivo, but the transcriptional yields tend to be quite low. Here we show that mutations previously shown to increase the thermal tolerance of T7 RNA polymerase can increase the activity of mutants with expanded substrate range. The resulting polymerase mutants can be used to generate 2'-O-methyl modified RNA with yields much higher than enzymes currently employed.

The origins of mutational robustness

Biological systems are resistant to genetic changes; a property known as mutational robustness, the origin of which remains an open question. In recent years, researchers have explored emergent properties of biological systems and mechanisms of genetic redundancy to reveal how mutational robustness emerges and persists. Several mechanisms have been proposed to explain the origin of mutational robustness, including molecular chaperones and gene duplication. The latter has received much attention, but its role in robustness remains controversial. Here, I examine recent findings linking genetic redundancy through gene duplication and mutational robustness. Experimental evolution and genome resequencing have made it possible to test the role of gene duplication in tolerating mutations at both the coding and regulatory levels. This evidence as well as previous findings on regulatory reprogramming of duplicates support the role of gene duplication in the origin of robustness.

Visualizing translocation dynamics and nascent transcript errors in paused RNA polymerases in vivo

BACKGROUND

Transcription elongation is frequently interrupted by pausing signals in DNA, with downstream effects on gene expression. Transcription errors also induce prolonged pausing, which can lead to a destabilized genome by interfering with DNA replication. Mechanisms of pausing associated with translocation blocks and misincorporation have been characterized in vitro, but not in vivo.

RESULTS

We investigate the pausing pattern of RNA polymerase (RNAP) in Escherichia coli by a novel approach, combining native elongating transcript sequencing (NET-seq) with RNase footprinting of the transcripts (RNET-seq). We reveal that the G-dC base pair at the 5' end of the RNA-DNA hybrid interferes with RNAP translocation. The distance between the 5' G-dC base pair and the 3' end of RNA fluctuates over a three-nucleotide width. Thus, the G-dC base pair can induce pausing in post-translocated, pre-translocated, and backtracked states of RNAP. Additionally, a CpG sequence of the template DNA strand spanning the active site of RNAP inhibits elongation and induces G-to-A errors, which leads to backtracking of RNAP. Gre factors efficiently proofread the errors and rescue the backtracked complexes. We also find that pausing events are enriched in the 5' untranslated region and antisense transcription of mRNA genes and are reduced in rRNA genes.

CONCLUSIONS

In E. coli, robust transcriptional pausing involves RNAP interaction with G-dC at the upstream end of the RNA-DNA hybrid, which interferes with translocation. CpG DNA sequences induce transcriptional pausing and G-to-A errors.

RNA editing by T7 RNA polymerase bypasses InDel mutations causing unexpected phenotypic changes

DNA-dependent T7 RNA polymerase (T7 RNAP) is the most powerful tool for both gene expression and in vitro transcription. By using a Next Generation Sequencing (NGS) approach we have analyzed the polymorphism of a T7 RNAP-generated mRNA pool of the mboIIM2 gene. We find that the enzyme displays a relatively high level of template-dependent transcriptional infidelity. The nucleotide misincorporations and multiple insertions in A/T-rich tracts of homopolymers in mRNA (0.20 and 0.089%, respectively) cause epigenetic effects with significant impact on gene expression that is disproportionally high to their frequency of appearance. The sequence-dependent rescue of single and even double InDel frameshifting mutants and wild-type phenotype recovery is observed as a result. As a consequence, a heterogeneous pool of functional and non-functional proteins of almost the same molecular mass is produced where the proteins are indistinguishable from each other upon ordinary analysis. We suggest that transcriptional infidelity as a general feature of the most effective RNAPs may serve to repair and/or modify a protein function, thus increasing the repertoire of phenotypic variants, which in turn has a high evolutionary potential.

Lost in transcription: transient errors in information transfer

Errors in information transfer from DNA to RNA to protein are inevitable. Here, we focus on errors that occur in nascent transcripts during transcription, epimutations. Recent approaches using novel cDNA library preparation and next-generation sequencing begin to directly determine the rate of epimutation and allow analysis of the epimutational spectrum of transcription errors, the type and sequence context of the errors produced in a transcript by an RNA polymerase. The phenotypic consequences of transcription errors have been assessed using both forward and reverse epimutation systems. These studies reveal that transient transcription errors can produce a modification of cell phenotype, partial phenotypic suppression of a mutant allele, and a heritable change in cell phenotype, epigenetic switching in a bistable gene network.

On the Nature and Evolutionary Impact of Phenotypic Robustness Mechanisms

Biologists have long observed that physiological and developmental processes are insensitive, or robust, to many genetic and environmental perturbations. A complete understanding of the evolutionary causes and consequences of this robustness is lacking. Recent progress has been made in uncovering the regulatory mechanisms that underlie environmental robustness in particular. Less is known about robustness to the effects of mutations, and indeed the evolution of mutational robustness remains a controversial topic. The controversy has spread to related topics, in particular the evolutionary relevance of cryptic genetic variation. This review aims to synthesize current understanding of robustness mechanisms and to cut through the controversy by shedding light on what is and is not known about mutational robustness. Some studies have confused mutational robustness with non-additive interactions between mutations (epistasis). We conclude that a profitable way forward is to focus investigations (and rhetoric) less on mutational robustness and more on epistasis.

Correlated occurrence and bypass of frame-shifting insertion-deletions (InDels) to give functional proteins

Short insertions and deletions (InDels) comprise an important part of the natural mutational repertoire. InDels are, however, highly deleterious, primarily because two-thirds result in frame-shifts. Bypass through slippage over homonucleotide repeats by transcriptional and/or translational infidelity is known to occur sporadically. However, the overall frequency of bypass and its relation to sequence composition remain unclear. Intriguingly, the occurrence of InDels and the bypass of frame-shifts are mechanistically related - occurring through slippage over repeats by DNA or RNA polymerases, or by the ribosome, respectively. Here, we show that the frequency of frame-shifting InDels, and the frequency by which they are bypassed to give full-length, functional proteins, are indeed highly correlated. Using a laboratory genetic drift, we have exhaustively mapped all InDels that occurred within a single gene. We thus compared the naive InDel repertoire that results from DNA polymerase slippage to the frame-shifting InDels tolerated following selection to maintain protein function. We found that InDels repeatedly occurred, and were bypassed, within homonucleotide repeats of 3–8 bases. The longer the repeat, the higher was the frequency of InDels formation, and the more frequent was their bypass. Besides an expected 8A repeat, other types of repeats, including short ones, and G and C repeats, were bypassed. Although obtained in vitro, our results indicate a direct link between the genetic occurrence of InDels and their phenotypic rescue, thus suggesting a potential role for frame-shifting InDels as bridging evolutionary intermediates.

Homonucleotide repeats are exceptionally prone to insertions and/or deletions of bases (InDels). However, unless they occur in a multiplicity of 3 bases, InDels disrupt the reading frame and are thus expected to be purged from coding regions. Homonucleotide repeats, however, are also vulnerable to slippage by RNA polymerases and the ribosome. Using laboratory evolution techniques, we systematically mapped the occurrence of InDels within a given gene, before and after selection. Our data indicate that frame-shifting InDels were frequently bypassed to give functional proteins at surprisingly high frequencies. Further, we found a strict correlation between the repeat length, the frequency of occurrence of InDels at the DNA level, and the likelihood of bypass by transcriptional/translational slippage. Our results suggest that frame-shifting InDels might comprise functional evolutionary intermediates, and an effective mean of sequence divergence (e.g. when an adjacent InDel restores the frame, resulting in altered sequence and, potentially, in an altered protein structure).

Direct assessment of transcription fidelity by high-resolution RNA sequencing

Cancerous and aging cells have long been thought to be impacted by transcription errors that cause genetic and epigenetic changes. Until now, a lack of methodology for directly assessing such errors hindered evaluation of their impact to the cells. We report a high-resolution Illumina RNA-seq method that can assess noncoded base substitutions in mRNA at 10⁻⁴–10⁻⁵ per base frequencies in vitro and in vivo. Statistically reliable detection of changes in transcription fidelity through ∼10³ nt DNA sites assures that the RNA-seq can analyze the fidelity in a large number of the sites where errors occur. A combination of the RNA-seq and biochemical analyses of the positions for the errors revealed two sequence-specific mechanisms that increase transcription fidelity by Escherichia coli RNA polymerase: (i) enhanced suppression of nucleotide misincorporation that improves selectivity for the cognate substrate, and (ii) increased backtracking of the RNA polymerase that decreases a chance of error propagation to the full-length transcript after misincorporation and provides an opportunity to proofread the error. This method is adoptable to a genome-wide assessment of transcription fidelity.

Experimental evidence for the thermophilicity of ancestral life

Theoretical studies have focused on the environmental temperature of the universal common ancestor of life with conflicting conclusions. Here we provide experimental support for the existence of a thermophilic universal common ancestor. We present the thermal stabilities and catalytic efficiencies of nucleoside diphosphate kinases (NDK), designed using the information contained in predictive phylogenetic trees, that seem to represent the last common ancestors of Archaea and of Bacteria. These enzymes display extreme thermal stabilities, suggesting thermophilic ancestries for Archaea and Bacteria. The results are robust to the uncertainties associated with the sequence predictions and to the tree topologies used to infer the ancestral sequences. Moreover, mutagenesis experiments suggest that the universal ancestor also possessed a very thermostable NDK. Because, as we show, the stability of an NDK is directly related to the environmental temperature of its host organism, our results indicate that the last common ancestor of extant life was a thermophile that flourished at a very high temperature.

Quantifying the adaptive potential of an antibiotic resistance enzyme

For a quantitative understanding of the process of adaptation, we need to understand its "raw material," that is, the frequency and fitness effects of beneficial mutations. At present, most empirical evidence suggests an exponential distribution of fitness effects of beneficial mutations, as predicted for Gumbel-domain distributions by extreme value theory. Here, we study the distribution of mutation effects on cefotaxime (Ctx) resistance and fitness of 48 unique beneficial mutations in the bacterial enzyme TEM-1 β-lactamase, which were obtained by screening the products of random mutagenesis for increased Ctx resistance. Our contributions are threefold. First, based on the frequency of unique mutations among more than 300 sequenced isolates and correcting for mutation bias, we conservatively estimate that the total number of first-step mutations that increase Ctx resistance in this enzyme is 87 [95% CI 75-189], or 3.4% of all 2,583 possible base-pair substitutions. Of the 48 mutations, 10 are synonymous and the majority of the 38 non-synonymous mutations occur in the pocket surrounding the catalytic site. Second, we estimate the effects of the mutations on Ctx resistance by determining survival at various Ctx concentrations, and we derive their fitness effects by modeling reproduction and survival as a branching process. Third, we find that the distribution of both measures follows a Fréchet-type distribution characterized by a broad tail of a few exceptionally fit mutants. Such distributions have fundamental evolutionary implications, including an increased predictability of evolution, and may provide a partial explanation for recent observations of striking parallel evolution of antibiotic resistance.

Comment on "Widespread RNA and DNA sequence differences in the human transcriptome"

Li et al. (Research Articles, 1 July 2011, p. 53; published online 19 May 2011) reported large numbers of differences between DNA and messenger RNA in human cells, indicating unprecedented levels of RNA editing, and including sequence changes not produced by any of the known RNA editing mechanisms. However, common sources of systematic errors in high-throughput sequencing technology, which were not properly accounted for in this study, explain most of the claimed differences.

Level of gene expression is a major determinant of protein evolution in the viral order Mononegavirales

Although the rate at which proteins change is a key parameter in molecular evolution, its determinants are poorly understood in viruses. A variety of factors, including gene length, codon usage bias, protein abundance, protein function, and gene expression level, have been shown to affect the rate of protein evolution in a diverse array of organisms. However, the role of these factors in viral evolution has yet to be addressed. The polar 3'-5' stepwise attenuation of transcription in the Mononegavirales, a group of single-strand negative-sense RNA viruses, provides a unique system to explore the determinants of protein evolution in viruses. We analyzed the relative importance of a variety of factors in shaping patterns of sequence variation in full-length genomes from 13 Mononegavirales species. Our analysis suggests that the level of gene expression, and by extension the relative genomic position of each gene, is a key determinant of the protein evolution in these viruses. This appears to be the consequence of selection for translational robustness, but not for translational accuracy, in highly expressed genes. The small genome size and number of proteins encoded by these viruses allowed us to identify other protein-specific factors that may also play a role in virus evolution, such as host-virus interactions and functional constraints. Finally, we explored the evolutionary pressures acting on noncoding regions in Mononegavirales genomes and observed that, despite being less constrained than coding regions, their evolutionary rates are also associated with genomic position.

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.