Mutational robustness and the role of buffer genes in evolvability.

Organisms rely on mutations to fuel adaptive evolution. However, many mutations impose a negative effect on fitness. Cells may have therefore evolved mechanisms that affect the phenotypic effects of mutations, thus conferring mutational robustness. Specifically, so-called buffer genes are hypothesized to interact directly or indirectly with genetic variation and reduce its effect on fitness. Environmental or genetic perturbations can change the interaction between buffer genes and genetic variation, thereby unmasking the genetic variation's phenotypic effects and thus providing a source of variation for natural selection to act on. This review provides an overview of our understanding of mutational robustness and buffer genes, with the chaperone gene HSP90 as a key example. It discusses whether buffer genes merely affect standing variation or also interact with de novo mutations, how mutational robustness could influence evolution, and whether mutational robustness might be an evolved trait or rather a mere side-effect of complex genetic interactions.

Promoters Constrain Evolution of Expression Levels of Essential Genes in Escherichia coli.

Variability in expression levels in response to random genomic mutations varies among genes, influencing both the facilitation and constraint of phenotypic evolution in organisms. Despite its importance, both the underlying mechanisms and evolutionary origins of this variability remain largely unknown due to the mixed contributions of cis- and trans-acting elements. To address this issue, we focused on the mutational variability of cis-acting elements, that is, promoter regions, in Escherichia coli. Random mutations were introduced into the natural and synthetic promoters to generate mutant promoter libraries. By comparing the variance in promoter activity of these mutant libraries, we found no significant difference in mutational variability in promoter activity between promoter groups, suggesting the absence of a signature of natural selection for mutational robustness. In contrast, the promoters controlling essential genes exhibited a remarkable bias in mutational variability, with mutants displaying higher activities than the wild types being relatively rare compared to those with lower activities. Our evolutionary simulation on a rugged fitness landscape provided a rationale for this vulnerability. These findings suggest that past selection created nonuniform mutational variability in promoters biased toward lower activities of random mutants, which now constrains the future evolution of downstream essential genes toward higher expression levels.

Quantifying the Mutational Robustness of Protein-Coding Genes.

We use large-scale mutagenesis data and computer simulations to quantify the mutational robustness of protein-coding genes by taking into account constraints arising from protein function and the genetic code. Analyses of the distribution of amino acid substitutions from 18 mutagenesis studies revealed an average of 45% of neutral variants; while mutagenesis data of 12 proteins artificially designed under no other constraints but stability, reach an average of 60%. Simulations using a lattice protein model allow us to contrast these estimates to the expected mutational robustness of protein families by generating unbiased samples of foldable sequences, which we find to have 30% of neutral variants. In agreement with mutagenesis data of designed proteins, the model shows that maximally robust protein families might access up to twice the amount of neutral variants observed in the unbiased samples (i.e. 60%). A biophysical model of protein-ligand binding suggests that constraints associated to molecular function have only a moderate impact on robustness of approximately 5 to 10% of neutral variants; and that the direction of this effect depends on the relation between functional performance and thermodynamic stability. Although the genetic code constraints the access of a gene's nucleotide sequence to only 30% of the full distribution of amino acid mutations, it provides an extra 15 to 20% of neutral variants to the estimations above, such that the expected, observed, and maximal robustness of protein-coding genes are approximately 50, 65, and 75%, respectively. We discuss our results in the light of three main hypothesis put forward to explain the existence of mutationally robust genes.

Large-scale in silico mutagenesis experiments reveal optimization of genetic code and codon usage for protein mutational robustness.

BACKGROUND: How, and the extent to which, evolution acts on DNA and protein sequences to ensure mutational robustness and evolvability is a long-standing open question in the field of molecular evolution. We addressed this issue through the first structurome-scale computational investigation, in which we estimated the change in folding free energy upon all possible single-site mutations introduced in more than 20,000 protein structures, as well as through available experimental stability and fitness data. RESULTS: At the amino acid level, we found the protein surface to be more robust against random mutations than the core, this difference being stronger for small proteins. The destabilizing and neutral mutations are more numerous in the core and on the surface, respectively, whereas the stabilizing mutations are about 4% in both regions. At the genetic code level, we observed smallest destabilization for mutations that are due to substitutions of base III in the codon, followed by base I, bases I+III, base II, and other multiple base substitutions. This ranking highly anticorrelates with the codon-anticodon mispairing frequency in the translation process. This suggests that the standard genetic code is optimized to limit the impact of random mutations, but even more so to limit translation errors. At the codon level, both the codon usage and the usage bias appear to optimize mutational robustness and translation accuracy, especially for surface residues. CONCLUSION: Our results highlight the non-universality of mutational robustness and its multiscale dependence on protein features, the structure of the genetic code, and the codon usage. Our analyses and approach are strongly supported by available experimental mutagenesis data.

Moderate Amounts of Epistasis are Not Evolutionarily Stable in Small Populations.

High mutation rates select for the evolution of mutational robustness where populations inhabit flat fitness peaks with little epistasis, protecting them from lethal mutagenesis. Recent evidence suggests that a different effect protects small populations from extinction via the accumulation of deleterious mutations. In drift robustness, populations tend to occupy peaks with steep flanks and positive epistasis between mutations. However, it is not known what happens when mutation rates are high and population sizes are small at the same time. Using a simple fitness model with variable epistasis, we show that the equilibrium fitness has a minimum as a function of the parameter that tunes epistasis, implying that this critical point is an unstable fixed point for evolutionary trajectories. In agent-based simulations of evolution at finite mutation rate, we demonstrate that when mutations can change epistasis, trajectories with a subcritical value of epistasis evolve to decrease epistasis, while those with supercritical initial points evolve towards higher epistasis. These two fixed points can be identified with mutational and drift robustness, respectively.

Chikungunya Virus Vaccine Candidates with Decreased Mutational Robustness Are Attenuated In Vivo and Have Compromised Transmissibility.

Chikungunya virus (CHIKV) is a reemerged arbovirus, a member of the Togaviridae family. It circulates through mosquito vectors mainly of the Aedes family and a mammalian host. CHIKV causes chikungunya fever, a mild to severe disease characterized by arthralgia, with some fatal outcomes described. In the past years, several outbreaks mainly caused by enhanced adaptation of the virus to the vector and ineffective control of the contacts between infected mosquito populations and the human host have been reported. Vaccines represent the best solution for the control of insect-borne viruses, including CHIKV, but are often unavailable. We designed live attenuated CHIKVs by applying a rational genomic design based on multiple replacements of synonymous codons. In doing so, the virus mutational robustness (capacity to maintain phenotype despite introduction of mutations to genotype) is decreased, driving the viral population toward deleterious evolutionary trajectories. When the candidate viruses were tested in the insect and mammalian hosts, we observed overall strong attenuation in both and greatly diminished signs of disease. Moreover, we found that the vaccine candidates elicited protective immunity related to the production of neutralizing antibodies after a single dose. During an experimental transmission cycle between mosquitoes and naive mice, vaccine candidates could be transmitted by mosquito bite, leading to asymptomatic infection in mice with compromised dissemination. Using deep-sequencing technology, we observed an increase in detrimental (stop) codons, which confirmed the effectiveness of this genomic design. Because the approach involves hundreds of synonymous modifications to the genome, the reversion risk is significantly reduced, rendering the viruses promising vaccine candidates.IMPORTANCE Chikungunya fever is a debilitating disease that causes severe pain to the joints, which can compromise the patient's lifestyle for several months and even in some grave cases lead to death. The etiological agent is chikungunya virus, an alphavirus transmitted by mosquito bite. Currently, there are no approved vaccines or treatments against the disease. In our research, we developed novel live attenuated vaccine candidates against chikungunya virus by applying an innovative genomic design. When tested in the insect and mammalian host, the vaccine candidates did not cause disease, elicited strong protection against further infection, and had low risk of reversion to pathogenic phenotypes.

Determining the factors driving selective effects of new nonsynonymous mutations.

The distribution of fitness effects (DFE) of new mutations plays a fundamental role in evolutionary genetics. However, the extent to which the DFE differs across species has yet to be systematically investigated. Furthermore, the biological mechanisms determining the DFE in natural populations remain unclear. Here, we show that theoretical models emphasizing different biological factors at determining the DFE, such as protein stability, back-mutations, species complexity, and mutational robustness make distinct predictions about how the DFE will differ between species. Analyzing amino acid-changing variants from natural populations in a comparative population genomic framework, we find that humans have a higher proportion of strongly deleterious mutations than Drosophila melanogaster. Furthermore, when comparing the DFE across yeast, Drosophila, mice, and humans, the average selection coefficient becomes more deleterious with increasing species complexity. Last, pleiotropic genes have a DFE that is less variable than that of nonpleiotropic genes. Comparing four categories of theoretical models, only Fisher's geometrical model (FGM) is consistent with our findings. FGM assumes that multiple phenotypes are under stabilizing selection, with the number of phenotypes defining the complexity of the organism. Our results suggest that long-term population size and cost of complexity drive the evolution of the DFE, with many implications for evolutionary and medical genomics.

Mapping the Peaks: Fitness Landscapes of the Fittest and the Flattest.

Populations exposed to a high mutation rate harbor abundant deleterious genetic variation, leading to depressed mean fitness. This reduction in mean fitness presents an opportunity for selection to restore fitness through the evolution of mutational robustness. In extreme cases, selection for mutational robustness can lead to flat genotypes (with low fitness but high robustness) outcompeting fit genotypes (with high fitness but low robustness)-a phenomenon known as survival of the flattest. While this effect was previously explored using the digital evolution system Avida, a complete analysis of the local fitness landscapes of fit and flat genotypes has been lacking, leading to uncertainty about the genetic basis of the survival-of-the-flattest effect. Here, we repeated the survival-of-the-flattest study and analyzed the mutational neighborhoods of fit and flat genotypes. We found that the flat genotypes, compared to the fit genotypes, had a reduced likelihood of deleterious mutations as well as an increased likelihood of neutral and, surprisingly, of lethal mutations. This trend holds for mutants one to four substitutions away from the wild-type sequence. We also found that flat genotypes have, on average, no epistasis between mutations, while fit genotypes have, on average, positive epistasis. Our results demonstrate that the genetic causes of mutational robustness on complex fitness landscapes are multifaceted. While the traditional idea of the survival of the flattest emphasized the evolution of increased neutrality, others have argued for increased mutational sensitivity in response to strong mutational loads. Our results show that both increased neutrality and increased lethality can lead to the evolution of mutational robustness. Furthermore, strong negative epistasis is not required for mutational sensitivity to lead to mutational robustness. Overall, these results suggest that mutational robustness is achieved by minimizing heritable deleterious variation.

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.

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.

Empirical analysis of RNA robustness and evolution using high-throughput sequencing of ribozyme reactions.

RNA molecules provide a realistic but tractable model of a genotype to phenotype relationship. This relationship has been extensively investigated computationally using secondary structure prediction algorithms. Enzymatic RNA molecules, or ribozymes, offer access to genotypic and phenotypic information in the laboratory. Advancements in high-throughput sequencing technologies have enabled the analysis of sequences in the lab that now rivals what can be accomplished computationally. This has motivated a resurgence of in vitro selection experiments and opened new doors for the analysis of the distribution of RNA functions in genotype space. A body of computational experiments has investigated the persistence of specific RNA structures despite changes in the primary sequence, and how this mutational robustness can promote adaptations. This article summarizes recent approaches that were designed to investigate the role of mutational robustness during the evolution of RNA molecules in the laboratory, and presents theoretical motivations, experimental methods and approaches to data analysis.

Selection for distinct gene expression properties favours the evolution of mutational robustness in gene regulatory networks.

Mutational robustness is a genotype's tendency to keep a phenotypic trait with little and few changes in the face of mutations. Mutational robustness is both ubiquitous and evolutionarily important as it affects in different ways the probability that new phenotypic variation arises. Understanding the origins of robustness is specially relevant for systems of development that are phylogenetically widespread and that construct phenotypic traits with a strong impact on fitness. Gene regulatory networks are examples of this class of systems. They comprise sets of genes that, through cross-regulation, build the gene activity patterns that define cellular responses, different tissues or distinct cell types. Several empirical observations, such as a greater robustness of wild-type phenotypes, suggest that stabilizing selection underlies the evolution of mutational robustness. However, the role of selection in the evolution of robustness is still under debate. Computer simulations of the dynamics and evolution of gene regulatory networks have shown that selection for any gene activity pattern that is steady and self-sustaining is sufficient to promote the evolution of mutational robustness. Here, I generalize this scenario using a computational model to show that selection for different aspects of a gene activity phenotype increases mutational robustness. Mutational robustness evolves even when selection favours properties that conflict with the stationarity of a gene activity pattern. The results that I present support an important role for stabilizing selection in the evolution of robustness in gene regulatory networks.

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.

Reducing Epistasis and Pleiotropy Can Avoid the Survival of the Flattest Tragedy.

This study investigates whether reducing epistasis and pleiotropy enhances mutational robustness in evolutionary adaptation, utilizing an indirect encoded model within the "survival of the flattest" (SoF) fitness landscape. By simulating genetic variations and their phenotypic consequences, we explore organisms' adaptive mechanisms to maintain positions on higher, narrower evolutionary peaks amidst environmental and genetic pressures. Our results reveal that organisms can indeed sustain their advantageous positions by minimizing the complexity of genetic interactions-specifically, by reducing the levels of epistasis and pleiotropy. This finding suggests a counterintuitive strategy for evolutionary stability: simpler genetic architectures, characterized by fewer gene interactions and multifunctional genes, confer a survival advantage by enhancing mutational robustness. This study contributes to our understanding of the genetic underpinnings of adaptability and robustness, challenging traditional views that equate complexity with fitness in dynamic environments.

Maximum mutational robustness in genotype-phenotype maps follows a self-similar blancmange-like curve.

Phenotype robustness, defined as the average mutational robustness of all the genotypes that map to a given phenotype, plays a key role in facilitating neutral exploration of novel phenotypic variation by an evolving population. By applying results from coding theory, we prove that the maximum phenotype robustness occurs when genotypes are organized as bricklayer's graphs, so-called because they resemble the way in which a bricklayer would fill in a Hamming graph. The value of the maximal robustness is given by a fractal continuous everywhere but differentiable nowhere sums-of-digits function from number theory. Interestingly, genotype-phenotype maps for RNA secondary structure and the hydrophobic-polar (HP) model for protein folding can exhibit phenotype robustness that exactly attains this upper bound. By exploiting properties of the sums-of-digits function, we prove a lower bound on the deviation of the maximum robustness of phenotypes with multiple neutral components from the bricklayer's graph bound, and show that RNA secondary structure phenotypes obey this bound. Finally, we show how robustness changes when phenotypes are coarse-grained and derive a formula and associated bounds for the transition probabilities between such phenotypes.

Analysis of the Contribution of Intrinsic Disorder in Shaping Potyvirus Genetic Diversity.

Intrinsically disordered regions (IDRs) are abundant in the proteome of RNA viruses. The multifunctional properties of these regions are widely documented and their structural flexibility is associated with the low constraint in their amino acid positions. Therefore, from an evolutionary stand point, these regions could have a greater propensity to accumulate non-synonymous mutations (NS) than highly structured regions (ORs, or 'ordered regions'). To address this hypothesis, we compared the distribution of non-synonymous mutations (NS), which we relate here to mutational robustness, in IDRs and ORs in the genome of potyviruses, a major genus of plant viruses. For this purpose, a simulation model was built and used to distinguish a possible selection phenomenon in the biological datasets from randomly generated mutations. We analyzed several short-term experimental evolution datasets. An analysis was also performed on the natural diversity of three different species of potyviruses reflecting their long-term evolution. We observed that the mutational robustness of IDRs is significantly higher than that of ORs. Moreover, the substitutions in the ORs are very constrained by the conservation of the physico-chemical properties of the amino acids. This feature is not found in the IDRs where the substitutions tend to be more random. This reflects the weak structural constraints in these regions, wherein an amino acid polymorphism is naturally conserved. In the course of evolution, potyvirus IDRs and ORs follow different evolutive paths with respect to their mutational robustness. These results have forced the authors to consider the hypothesis that IDRs and their associated amino acid polymorphism could constitute a potential adaptive reservoir.

Robustness and the generalist niche of polyploid species: Genome shock or gradual evolution?

The prevalence of polyploidy in wild and crop species has stimulated debate over its evolutionary advantages and disadvantages. Previous studies have focused on changes occurring at the polyploidization events, including genome-wide changes termed "genome shock," as well as ancient polyploidy. Recent bioinformatics advances and empirical studies of Arabidopsis and wheat relatives are filling a research gap: the functional evolutionary study of polyploid species using RNA-seq, DNA polymorphism, and epigenomics. Polyploid species can become generalists in natura through environmental robustness by inheriting and merging parental stress responses. Their evolvability is enhanced by mutational robustness working on inherited standing variation. The identification of key genes responsible for gradual adaptive evolution will encourage synthetic biological approaches to transfer polyploid advantages to other species.

The evolutionary potential of influenza A virus hemagglutinin is highly constrained by epistatic interactions with neuraminidase.

Antigenic evolution of the influenza A virus (IAV) hemagglutinin (HA) gene limits efforts to effectively control the spread of the virus in the population. Efforts to understand the mechanisms governing HA antigenic evolution typically examine the HA gene in isolation. This can ignore the importance of balancing HA receptor binding activities with the receptor-destroying activities of the viral neuraminidase (NA) to maintain viral fitness. We hypothesize that the need to maintain functional balance with NA significantly constrains the evolutionary potential of the HA. We use deep mutational scanning and show that variation in NA activity significantly reshapes the HA fitness landscape by modulating the overall mutational robustness of HA. Consistent with this, we observe that different NA backgrounds support the emergence of distinct repertoires of HA escape variants under neutralizing antibody pressure. Our results reveal a critical role for intersegment epistasis in influencing the evolutionary potential of the HA gene.

Altering Compositional Properties of Viral Genomes to Design Live-Attenuated Vaccines.

Live-attenuated vaccines have been historically used to successfully prevent numerous diseases caused by a broad variety of RNA viruses due to their ability to elicit strong and perdurable immune-protective responses. In recent years, various strategies have been explored to achieve viral attenuation by rational genetic design rather than using classic and empirical approaches, based on successive passages in cell culture. A deeper understanding of evolutionary implications of distinct viral genomic compositional aspects, as well as substantial advances in synthetic biology technologies, have provided a framework to achieve new viral attenuation strategies. Herein, we will discuss different approaches that are currently applied to modify compositional features of viruses in order to develop novel live-attenuated vaccines.

A Mathematical Analysis of HDV Genotypes: From Molecules to Cells.

Hepatitis D virus (HDV) is classified according to eight genotypes. The various genotypes are included in the HDVdb database, where each HDV sequence is specified by its genotype. In this contribution, a mathematical analysis is performed on RNA sequences in HDVdb. The RNA folding predicted structures of the Genbank HDV genome sequences in HDVdb are classified according to their coarse-grain tree-graph representation. The analysis allows discarding in a simple and efficient way the vast majority of the sequences that exhibit a rod-like structure, which is important for the virus replication, to attempt to discover other biological functions by structure consideration. After the filtering, there remain only a small number of sequences that can be checked for their additional stem-loops besides the main one that is known to be responsible for virus replication. It is found that a few sequences contain an additional stem-loop that is responsible for RNA editing or other possible functions. These few sequences are grouped into two main classes, one that is well-known experimentally belonging to genotype 3 for patients from South America associated with RNA editing, and the other that is not known at present belonging to genotype 7 for patients from Cameroon. The possibility that another function besides virus replication reminiscent of the editing mechanism in HDV genotype 3 exists in HDV genotype 7 has not been explored before and is predicted by eigenvalue analysis. Finally, when comparing native and shuffled sequences, it is shown that HDV sequences belonging to all genotypes are accentuated in their mutational robustness and thermodynamic stability as compared to other viruses that were subjected to such an analysis.