Perspective: Evolution and detection of genetic robustness

Robustness promotes evolvability of thermotolerance in an RNA virus

BACKGROUND

The ability for an evolving population to adapt to a novel environment is achieved through a balance of robustness and evolvability. Robustness is the invariance of phenotype in the face of perturbation and evolvability is the capacity to adapt in response to selection. Genetic robustness has been posited, depending on the underlying mechanism, to either decrease the efficacy of selection, or increase the possibility of future adaptation. However, the true effect of genetic robustness on evolvability in biological systems remains uncertain.

RESULTS

Here we demonstrate that genetic robustness increases evolvability of thermotolerance in laboratory populations of the RNA virus phi6. We observed that populations founded by robust clones evolved greater resistance to heat shock, relative to populations founded by brittle (less-robust) clones. Thus, we provide empirical evidence for the idea that robustness can promote evolvability in this environment, and further suggest that evolvability can arise indirectly via selection for robustness, rather than through direct selective action.

CONCLUSION

Our data imply that greater tolerance of mutational change is associated with virus adaptability in a new niche, a finding generally relevant to evolutionary biology, and informative for elucidating how viruses might evolve to emerge in new habitats and/or overcome novel therapies.

Architecture and inherent robustness of a bacterial cell-cycle control system

A closed-loop control system drives progression of the coupled stalked and swarmer cell cycles of the bacterium Caulobacter crescentus in a near-mechanical step-like fashion. The cell-cycle control has a cyclical genetic circuit composed of four regulatory proteins with tight coupling to processive chromosome replication and cell division subsystems. We report a hybrid simulation of the coupled cell-cycle control system, including asymmetric cell division and responses to external starvation signals, that replicates mRNA and protein concentration patterns and is consistent with observed mutant phenotypes. An asynchronous sequential digital circuit model equivalent to the validated simulation model was created. Formal model-checking analysis of the digital circuit showed that the cell-cycle control is robust to intrinsic stochastic variations in reaction rates and nutrient supply, and that it reliably stops and restarts to accommodate nutrient starvation. Model checking also showed that mechanisms involving methylation-state changes in regulatory promoter regions during DNA replication increase the robustness of the cell-cycle control. The hybrid cell-cycle simulation implementation is inherently extensible and provides a promising approach for development of whole-cell behavioral models that can replicate the observed functionality of the cell and its responses to changing environmental conditions.

Quantifying the decanalizing effects of spontaneous mutations in rhabditid nematodes

The evolution of canalization, the robustness of the phenotype to environmental or genetic perturbation, has attracted considerable recent interest. A key step toward understanding the evolution of any phenotype is characterizing the rate at which mutation introduces genetic variation for the trait (the mutational variance, V(M)) and the average directional effects of mutations on the trait mean (DeltaM). In this study, the mutational parameters for canalization of productivity and body volume are quantified in two sets of mutation accumulation lines of nematodes in the genus Caenorhabditis and are compared with the mutational parameters for the traits themselves. Four results emerge: (1) spontaneous mutations consistently decanalize the phenotype; (2) the mutational parameters for decanalization, V(M) (quantified as mutational heritability) and DeltaM, are of the same order of magnitude as the same parameters for the traits themselves; (3) the mutational parameters for canalization are roughly correlated with the parameters for the traits themselves across taxa; and (4) there is no evidence that residual segregating overdominant loci contribute to the decay of canalization. These results suggest that canalization is readily evolvable and that any evolutionary factor that causes mutations to accumulate will, on average, decanalize the phenotype.

Parasites lead to evolution of robustness against gene loss in host signaling networks

Many biological networks can maintain their function against single gene loss. However, the evolutionary mechanisms responsible for such robustness remain unclear. Here, we demonstrate that antagonistic host–parasite interactions can act as a selective pressure driving the emergence of robustness against gene loss. Using a model of host signaling networks and simulating their coevolution with parasites that interfere with network function, we find that networks evolve both redundancy and specific architectures that allow them to maintain their response despite removal of proteins. We show that when the parasite pressure is removed, subsequent evolution can lead to loss of redundancy while architecture-based robustness is retained. Contrary to intuition, increased parasite virulence hampers evolution of robustness by limiting the generation of population level diversity in the host. However, when robustness emerges under high virulence, it tends to be stronger. These findings predict an increased presence of robustness mechanisms in biological networks operating under parasite interference. Conversely, the presence of such mechanisms could indicate current or past parasite interference.

Critical dynamics in genetic regulatory networks: examples from four kingdoms

The coordinated expression of the different genes in an organism is essential to sustain functionality under the random external perturbations to which the organism might be subjected. To cope with such external variability, the global dynamics of the genetic network must possess two central properties. (a) It must be robust enough as to guarantee stability under a broad range of external conditions, and (b) it must be flexible enough to recognize and integrate specific external signals that may help the organism to change and adapt to different environments. This compromise between robustness and adaptability has been observed in dynamical systems operating at the brink of a phase transition between order and chaos. Such systems are termed critical. Thus, criticality, a precise, measurable, and well characterized property of dynamical systems, makes it possible for robustness and adaptability to coexist in living organisms. In this work we investigate the dynamical properties of the gene transcription networks reported for S. cerevisiae, E. coli, and B. subtilis, as well as the network of segment polarity genes of D. melanogaster, and the network of flower development of A. thaliana. We use hundreds of microarray experiments to infer the nature of the regulatory interactions among genes, and implement these data into the Boolean models of the genetic networks. Our results show that, to the best of the current experimental data available, the five networks under study indeed operate close to criticality. The generality of this result suggests that criticality at the genetic level might constitute a fundamental evolutionary mechanism that generates the great diversity of dynamically robust living forms that we observe around us.

Molecular chaperones and selection against mutations

Background

Molecular chaperones help to restore the native states of proteins after their destabilization by external stress. It has been proposed that another function of chaperones is to maintain the activity of proteins destabilized by mutation, weakening the selection against suboptimal protein variants. This would allow for the accumulation of genetic variation which could then be exposed during environmental perturbation and facilitate rapid adaptation.

Results

We focus on studies describing interactions of chaperones with mutated polypeptides. There are some examples that chaperones can alleviate the deleterious effects of mutations through increased assistance of destabilized proteins. These experiments are restricted to bacteria and typically involve overexpression of chaperones. In eukaryotes, it was found that the malfunctioning of chaperones aggravated phenotypic aberrations associated with mutations. This effect could not be linked to chaperone-mediated stabilization of mutated proteins. More likely, the insufficient activity of chaperones inflicted a deregulation of multiple cellular systems, including those responsible for signaling and therefore important in development. As to why the assistance of mutated proteins by chaperones seems difficult to demonstrate, we note that chaperone-assisted folding can often co-exist with chaperone-assisted degradation. There is growing evidence that some chaperones, including those dependent on Hsp90, can detect potentially functional but excessively unstable proteins and direct them towards degradation instead of folding. This implies that at least some mutations are exposed rather than masked by the activity of molecular chaperones.

Conclusion

It is at present impossible to determine whether molecular chaperones are mostly helpers or examiners of mutated proteins because experiments showing either of these roles are very few. Depending on whether assistance or disposal prevails, molecular chaperones could speed up or slow down evolution of protein sequences. Similar uncertainties arise when the concept of chaperones (mostly Hsp90) as general regulators of evolvability is considered. If the two roles of chaperones are antagonistic, then any (even small) modification of the chaperone activities to save mutated polypeptides could lead to increased misfolding and aggregation of other proteins. This would be a permanent burden, different from the stochastic cost arising from indiscriminate buffering of random mutations of which many are maladaptive.

Reviewers

This article was reviewed by A. S. Kondrashov, J. Höhfeld (nominated by A. Eyre-Walker) and D. A. Drummond (nominated by C. Adami). For the full reviews, please go to the Reviewers' comments section.

The evolution of developmental patterning under genetic duplication constraints

Of considerable interest are the evolutionary and developmental origins of complex, adaptive structures and the mechanisms that stabilize these structures. We consider the relationship between the evolutionary process of gene duplication and deletion and the stability of morphogenetic patterns produced by interacting activators and inhibitors. We compare the relative stability of patterns with a single activator and inhibitor (two-dimensional system) against a 'redundant' system with two activators or two inhibitors (three-dimensional system). We find that duplication events can both expand and contract the space of patterns. We study developmental robustness in terms of stochastic escape times from this space, also known as a 'canalization potential'. We embed the output of pattern formation into an explicit evolutionary model of gene duplication, gene loss and variation in the steepness of the canalization potential. We find that under all constant conditions, the system evolves towards a preference for steep potentials associated with low phenotypic variability and longer lifespans. This preference leads to an overall decrease in the density of redundant genotypes as developmental robustness neutralizes the advantages of genetic robustness.

Quasi-species evolution in subdivided populations favours maximally deleterious mutations

Most models of quasi-species evolution predict that populations will evolve to occupy areas of sequence space with the greatest concentration of neutral sequences, thus minimizing the deleterious mutation rate and creating mutationally 'robust' genomes. In contrast, empirical studies of the principal model of quasi-species evolution, RNA viruses, suggest that the effects of deleterious mutations are more severe than in similar DNA-based microbes. We demonstrate that populations divided into discrete patches connected by dispersal may favour genotypes where the deleterious effect of non-neutral mutations is maximized. This effect is especially strong in the absence of back mutation and when the amount of time spent in hosts prior to dispersal is intermediate. Our results indicate that RNA viruses that produce acute infections initiated by a small number of virions are expected to evolve fragile genetic architectures when compared with other RNA viruses.

In silico genetic robustness analysis of microRNA secondary structures: potential evidence of congruent evolution in microRNA

Background

Robustness is a fundamental property of biological systems and is defined as the ability to maintain stable functioning in the face of various perturbations. Understanding how robustness has evolved has become one of the most attractive areas of research for evolutionary biologists, as it is still unclear whether genetic robustness evolved as a direct consequence of natural selection, as an intrinsic property of adaptations, or as congruent correlate of environment robustness. Recent studies have demonstrated that the stem-loop structures of microRNA (miRNA) are tolerant to some structural changes and show thermodynamic stability. We therefore hypothesize that genetic robustness may evolve as a correlated side effect of the evolution for environmental robustness.

Results

We examine the robustness of 1,082 miRNA genes covering six species. Our data suggest the stem-loop structures of miRNA precursors exhibit a significantly higher level of genetic robustness, which goes beyond the intrinsic robustness of the stem-loop structure and is not a byproduct of the base composition bias. Furthermore, we demonstrate that the phenotype of miRNA buffers against genetic perturbations, and at the same time is also insensitive to environmental perturbations.

Conclusion

The results suggest that the increased robustness of miRNA stem-loops may result from congruent evolution for environment robustness. Potential applications of our findings are also discussed.

Maladaptation and the paradox of robustness in evolution

Background

Organisms use a variety of mechanisms to protect themselves against perturbations. For example, repair mechanisms fix damage, feedback loops keep homeostatic systems at their setpoints, and biochemical filters distinguish signal from noise. Such buffering mechanisms are often discussed in terms of robustness, which may be measured by reduced sensitivity of performance to perturbations.

Methodology/Principal Findings

I use a mathematical model to analyze the evolutionary dynamics of robustness in order to understand aspects of organismal design by natural selection. I focus on two characters: one character performs an adaptive task; the other character buffers the performance of the first character against perturbations. Increased perturbations favor enhanced buffering and robustness, which in turn decreases sensitivity and reduces the intensity of natural selection on the adaptive character. Reduced selective pressure on the adaptive character often leads to a less costly, lower performance trait.

Conclusions/Significance

The paradox of robustness arises from evolutionary dynamics: enhanced robustness causes an evolutionary reduction in the adaptive performance of the target character, leading to a degree of maladaptation compared to what could be achieved by natural selection in the absence of robustness mechanisms. Over evolutionary time, buffering traits may become layered on top of each other, while the underlying adaptive traits become replaced by cheaper, lower performance components. The paradox of robustness has widespread implications for understanding organismal design.

Direct selection on genetic robustness revealed in the yeast transcriptome

BACKGROUND

Evolutionary theory predicts that organisms should evolve the ability to produce high fitness phenotypes in the face of environmental disturbances (environmental robustness) or genetic mutations (genetic robustness). While several studies have uncovered mechanisms that lead to both environmental and genetic robustness, we have yet to understand why some components of the genome are more robust than others. According to evolutionary theory, environmental and genetic robustness will have different responses to selective forces. Selection on environmental robustness for a trait is expected to be strong and related to the fitness costs of altering that trait. In contrast to environmental robustness, selection on genetic robustness for a trait is expected to be largely independent of the fitness cost of altering the trait and instead should correlate with the standing genetic variation for the trait that can potentially be buffered. Several mechanisms that provide both environmental and genetic robustness have been described, and this correlation could be explained by direct selection on both forms of robustness (direct selection hypothesis), or through selection on environmental robustness and a correlated response in genetic robustness (congruence hypothesis).

METHODOLOGY/PRINCIPAL FINDINGS

Using both published and novel data on gene expression in the yeast Saccharomyces cerevisiae, we find that genetic robustness is correlated with environmental robustness across the yeast genome as predicted by the congruence hypothesis. However, we also show that environmental robustness, but not genetic robustness, is related to per-gene fitness effects. In contrast, genetic robustness is significantly correlated with network position, suggesting that genetic robustness has been under direct selection.

CONCLUSIONS/SIGNIFICANCE

We observed a significant correlation between our measures of genetic and environmental robustness, in agreement with the congruence hypothesis. However, this correlation alone cannot explain the co-variance of genetic robustness with position in the protein interaction network. We therefore conclude that direct selection on robustness has played a role in the evolution of genetic robustness in the transcriptome.

RSRE: RNA structural robustness evaluator

Biological robustness, defined as the ability to maintain stable functioning in the face of various perturbations, is an important and fundamental topic in current biology, and has become a focus of numerous studies in recent years. Although structural robustness has been explored in several types of RNA molecules, the origins of robustness are still controversial. Computational analysis results are needed to make up for the lack of evidence of robustness in natural biological systems. The RNA structural robustness evaluator (RSRE) web server presented here provides a freely available online tool to quantitatively evaluate the structural robustness of RNA based on the widely accepted definition of neutrality. Several classical structure comparison methods are employed; five randomization methods are implemented to generate control sequences; sub-optimal predicted structures can be optionally utilized to mitigate the uncertainty of secondary structure prediction. With a user-friendly interface, the web application is easy to use. Intuitive illustrations are provided along with the original computational results to facilitate analysis. The RSRE will be helpful in the wide exploration of RNA structural robustness and will catalyze our understanding of RNA evolution. The RSRE web server is freely available at http://biosrv1.bmi.ac.cn/RSRE/ or http://biotech.bmi.ac.cn/RSRE/.

Selection for robustness in mutagenized RNA viruses

Mutational robustness is defined as the constancy of a phenotype in the face of deleterious mutations. Whether robustness can be directly favored by natural selection remains controversial. Theory and in silico experiments predict that, at high mutation rates, slow-replicating genotypes can potentially outcompete faster counterparts if they benefit from a higher robustness. Here, we experimentally validate this hypothesis, dubbed the "survival of the flattest," using two populations of the vesicular stomatitis RNA virus. Characterization of fitness distributions and genetic variability indicated that one population showed a higher replication rate, whereas the other was more robust to mutation. The faster replicator outgrew its robust counterpart in standard competition assays, but the outcome was reversed in the presence of chemical mutagens. These results show that selection can directly favor mutational robustness and reveal a novel viral resistance mechanism against treatment by lethal mutagenesis.

The theory of biological robustness and its implication in cancer

One of the essential issues in systems biology is to identify fundamental principles that govern living organisms at the system level. In this chapter, I argue that robustness is a fundamental feature of living systems where its relationship with evolution-trade-offs among robustness, fragility, resource demands, and performance-provides a possible framework for how biological systems have evolved and been organized. In addition, diseases can be con- sidered as a manifestation of fragility of the system. In some cases, such as cancer, the disease state establishes its own robustness against therapeutic interventions. Understanding robustness and its intrinsic properties will provide us with a more profound understanding of biological systems, their anomalies, and countermeasures.

Evolution of robustness to noise and mutation in gene expression dynamics

Phenotype of biological systems needs to be robust against mutation in order to sustain themselves between generations. On the other hand, phenotype of an individual also needs to be robust against fluctuations of both internal and external origins that are encountered during growth and development. Is there a relationship between these two types of robustness, one during a single generation and the other during evolution? Could stochasticity in gene expression have any relevance to the evolution of these types of robustness? Robustness can be defined by the sharpness of the distribution of phenotype; the variance of phenotype distribution due to genetic variation gives a measure of 'genetic robustness', while that of isogenic individuals gives a measure of 'developmental robustness'. Through simulations of a simple stochastic gene expression network that undergoes mutation and selection, we show that in order for the network to acquire both types of robustness, the phenotypic variance induced by mutations must be smaller than that observed in an isogenic population. As the latter originates from noise in gene expression, this signifies that the genetic robustness evolves only when the noise strength in gene expression is larger than some threshold. In such a case, the two variances decrease throughout the evolutionary time course, indicating increase in robustness. The results reveal how noise that cells encounter during growth and development shapes networks' robustness to stochasticity in gene expression, which in turn shapes networks' robustness to mutation. The necessary condition for evolution of robustness, as well as the relationship between genetic and developmental robustness, is derived quantitatively through the variance of phenotypic fluctuations, which are directly measurable experimentally.

Adaptive walks toward a moving optimum

We investigate how the dynamics and outcomes of adaptation by natural selection are affected by environmental stability by simulating adaptive walks in response to an environmental change of fixed magnitude but variable speed. Here we consider monomorphic lineages that adapt by the sequential fixation of beneficial mutations. This is modeled by selecting short RNA sequences for folding stability and secondary structure conservation at increasing temperatures. Using short RNA sequences allows us to describe adaptive outcomes in terms of genotype (sequence) and phenotype (secondary structure) and to follow the dynamics of fitness increase. We find that slower rates of environmental change affect the dynamics of adaptive walks by reducing the fitness effect of fixed beneficial mutations, as well as by increasing the range of time in which the substitutions of largest effect are likely to occur. In addition, adaptation to slower rates of environmental change results in fitter endpoints with fewer possible end phenotypes relative to lineages that adapt to a sudden change. This suggests that care should be taken when experiments using sudden environmental changes are used to make predictions about adaptive responses to gradual change.

Empirical support for optimal virulence in a castrating parasite

The trade-off hypothesis for the evolution of virulence predicts that parasite transmission stage production and host exploitation are balanced such that lifetime transmission success (LTS) is maximised. However, the experimental evidence for this prediction is weak, mainly because LTS, which indicates parasite fitness, has been difficult to measure. For castrating parasites, this simple model has been modified to take into account that parasites convert host reproductive resources into transmission stages. Parasites that kill the host too early will hardly benefit from these resources, while postponing the killing of the host results in diminished returns. As predicted from optimality models, a parasite inducing castration should therefore castrate early, but show intermediate levels of virulence, where virulence is measured as time to host killing. We studied virulence in an experimental system where a bacterial parasite castrates its host and produces spores that are not released until after host death. This permits estimating the LTS of the parasite, which can then be related to its virulence. We exposed replicate individual Daphnia magna (Crustacea) of one host clone to the same amount of bacterial spores and followed individuals until their death. We found that the parasite shows strong variation in the time to kill its host and that transmission stage production peaks at an intermediate level of virulence. A further experiment tested for the genetic basis of variation in virulence by comparing survival curves of daphniids infected with parasite spores obtained from early killing versus late killing infections. Hosts infected with early killer spores had a significantly higher death rate as compared to those infected with late killers, indicating that variation in time to death was at least in part caused by genetic differences among parasites. We speculate that the clear peak in lifetime reproductive success at intermediate killing times may be caused by the exceptionally strong physiological trade-off between host and parasite reproduction. This is the first experimental study to demonstrate that the production of propagules is highest at intermediate levels of virulence and that parasite genetic variability is available to drive the evolution of virulence in this system.

Balancing robustness and evolvability

Can a single unifying mathematical framework help to explain robustness - the ability of organisms to persist in the face of changing conditions - at all biological scales, from biochemical to ecological?

Evolution of complexity in signaling pathways

It is not clear how biological pathways evolve to mediate a certain physiological response and why they show a level of complexity that is generally above the minimum required to achieve such a response. One possibility is that pathway complexity increases due to the nature of evolutionary mechanisms. Here, we analyze this possibility by using mathematical models of biological pathways and evolutionary simulations. Starting with a population of small pathways of three proteins, we let the population evolve with mutations that affect pathway structure through duplication or deletion of existing proteins, deletion or creation of interactions among them, or addition of new proteins. Our simulations show that such mutational events, coupled with a selective pressure, leads to growth of pathways. These results indicate that pathways could be driven toward complexity via simple evolutionary mechanisms and that complexity can arise without any specific selective pressure for it. Furthermore, we find that the level of complexity that pathways evolve toward depends on the selection criteria. In general, we find that final pathway size tends to be lower when pathways evolve under stringent selection criteria. This leads to the counterintuitive conclusion that simple response requirements on a pathway would facilitate its evolution toward higher complexity.

Why molecular chaperones buffer mutational damage: a case study with a yeast Hsp40/70 system

The malfunctioning of molecular chaperones may result in uncovering genetic variation. The molecular basis of this phenomenon remains largely unknown. Chaperones rescue proteins unfolded by environmental stresses and therefore they might also help to stabilize mutated proteins and thus mask damages. To test this hypothesis, we carried out a genomewide mutagenesis followed by a screen for mutations that were synthetically harmful when the RAC-Ssb1/2 cytosolic chaperones were inactive. Mutants with such a phenotype were found and mapped to single nucleotide substitutions. However, neither the genes identified nor the nature of genetic lesions implied that folding of the mutated proteins was being supported by the chaperones. In a second screen, we identified temperature-sensitive (ts) mutants, a phenotype indicative of structural instability of proteins. We tested these for an association with sensitivity to loss of chaperone activity but found no such correlation as might have been expected if the chaperones assisted the folding of mutant proteins. Thus, molecular chaperones can mask the negative effects of mutations but the mechanism of such buffering need not be direct. A plausible role of chaperones is to stabilize genetic networks, thus making them more tolerant to malfunctioning of their constituents.

Subsets of SNPs define rare genotype classes that predict ischemic heart disease

Single nucleotide polymorphisms (SNPs) are hypothesized to explain the genetic predisposition to ischemic heart disease (IHD) in the general population. Lack of evidence for a role of such variation is fostering pessimism about the utility of genetic information in the practice of medicine. In this study we determined the utility of exonic and 5' SNPs in apolipoprotein E (APOE) and lipoprotein lipase (LPL) when considered singly and in combination for predicting incidence of IHD in 8,456 individuals from the general population during 24 years of follow-up. In men, LPL D9N improved prediction of IHD (P = 0.03) beyond smoking, diabetes and hypertension. The group of men heterozygous and homozygous for the rare D9N variant had a hazard ratio (HR) of 1.69 (95% confidence interval = 1.10-2.58) relative to the most common genotype. Pairwise combinations of D9N with -219G > T in APOE and N291S and S447X in LPL significantly improved the prediction of IHD (P = 0.05 in women, P = 0.04 in men, P = 0.03 in men, respectively) beyond smoking, diabetes and hypertension, and identified subgroups of individuals (n = 6-94) with highly significant HRs of 1.92-4.35. These results were validated in a case-control study (n = 8,806). In conclusion, we present evidence that combinations of SNPs in APOE and LPL identify subgroups of individuals at substantially increased risk of IHD beyond that associated with smoking, diabetes and hypertension.

Effect of random and hub gene disruptions on environmental and mutational robustness in Escherichia coli

Background

Genome-wide profiling has allowed the regulatory interaction networks of many organisms to be visualised and the pattern of connections between genes to be studied. These networks are non-random, following a power-law distribution with a small number of well-connected 'hubs' and many genes with only one or a few connections. Theoretical work predicts that power-law networks display several unique properties. One of the most biologically interesting of these is an intrinsic robustness to disturbance such that removal of a random gene will have little effect on network function. Conversely, targeted removal of a hub gene is expected to have a large effect.

Results

We compared the response of Escherichia coli to environmental and mutational stress following disruption of random or hub genes. We found that disruption of random genes had less effect on robustness to environmental stress than did the targeted disruption of hub genes. In contrast, random disruption strains were slightly less robust to the effect of mutational stress than were hub disruption strains. When we compared the effect of each disruption on environmental and mutational stress, we found a negative relationship, such that strains that were more environmentally robust tended to be less robust to mutational stress.

Conclusion

Our results demonstrate that mutant strains of E. coli respond differently to stress, depending on whether random or hub genes are disrupted. This difference indicates that the power-law distribution of regulatory interactions has biological significance, making random disruptions less deleterious to organisms facing environmental stress. That E. coli can reduce the effect of environmental stress without reducing the phenotypic effect of additional mutations, indicates that robustness and evolvability need not be antagonistic.

Effects of random mutations in the human immunodeficiency virus type 1 transcriptional promoter on viral fitness in different host cell environments

A mutation's effect on fitness or phenotype may in part depend on the interaction of the mutation with the environment. The resulting phenotype or fitness is important, since it determines the adaptive potential of a species. To date, most studies have focused on alterations to protein-coding regions of the genome and their consequential fitness effects. Non-protein-coding regulatory regions have been largely neglected, although they make up a large and important part of an organism's genome. Here, we use human immunodeficiency virus type 1 as a model system to investigate fitness effects of random mutations in noncoding DNA-binding sites of the transcriptional promoter. We determined 242 fitness values for 35 viral promoter mutants with one, two, or three mutations across seven distinct cellular environments and identified that (i) all mutants have an effect in at least one cellular environment; (ii) fitness effects are highly dependent on the cellular environment; (iii) disadvantageous and advantageous mutations occur at high and similar frequencies; and (iv) epistatic effects of multiple mutations are rare. Our results underline the evolutionary potential of regulatory regions and indicate that DNA-binding sites evolve under strong selection, while at the same time, they are very plastic to environmental change.

Geometric robustness theory and biological networks

We provide a geometric framework for investigating the robustness of information flows over biological networks. We use information measures to quantify the impact of knockout perturbations on simple networks. Robustness has two components, a measure of the causal contribution of a node or nodes, and a measure of the change or exclusion dependence, of the network following node removal. Causality is measured as statistical contribution of a node to network function, whereas exclusion dependence measures a distance between unperturbed network and reconfigured network function. We explore the role that redundancy plays in increasing robustness, and how redundacy can be exploited through error-correcting codes implemented by networks. We provide examples of the robustness measure when applied to familiar boolean functions such as the AND, OR and XOR functions. We discuss the relationship between robustness measures and related measures of complexity and how robustness always implies a minimal level of complexity.

Mechanisms of genetic robustness in RNA viruses

Two key features of RNA viruses are their compacted genomes and their high mutation rate. Accordingly, deleterious mutations are common and have an enormous impact on viral fitness. In their multicellular hosts, robustness can be achieved by genomic redundancy, including gene duplication, diploidy, alternative metabolic pathways and biochemical buffering mechanisms. However, here we review evidence suggesting that during RNA virus evolution, alternative robustness mechanisms may have been selected. After briefly describing how genetic robustness can be quantified, we discuss mechanisms of intrinsic robustness arising as consequences of RNA-genome architecture, replication peculiarities and quasi-species population dynamics. These intrinsic robustness mechanisms operate efficiently at the population level, despite the mutational sensitivity shown by individual genomes. Finally, we discuss the possibility that viruses might exploit cellular buffering mechanisms for their own benefit, producing a sort of extrinsic robustness.

Direct evolution of genetic robustness in microRNA

Genetic robustness, the invariance of the phenotype in the face of genetic perturbations, can endow the organism with reduced susceptibility to mutations. A large body of work in recent years has focused on the origins, mechanisms, and consequences of robustness in a wide range of biological systems. Despite the apparent prevalence of mutational robustness in nature, however, its evolutionary origins are still unclear. Does robustness evolve directly by natural selection or is it merely a correlated byproduct of other phenotypic traits? By examining microRNA (miRNA) genes of several eukaryotic species, we show that the structure of miRNA precursor stem-loops exhibits a significantly high level of mutational robustness in comparison with random RNA sequences with similar stem-loop structures. Hence, this excess robustness of miRNA goes beyond the intrinsic robustness of the stem-loop hairpin structure. Furthermore, we show that it is not the byproduct of a base composition bias or of thermodynamic stability. These findings suggest that the excess robustness of miRNA stem-loops is the result of direct evolutionary pressure toward increased robustness. We further demonstrate that this adaptive robustness evolves to compensate for structures with low intrinsic robustness.

Evolution of mutational robustness in an RNA virus

Mutational (genetic) robustness is phenotypic constancy in the face of mutational changes to the genome. Robustness is critical to the understanding of evolution because phenotypically expressed genetic variation is the fuel of natural selection. Nonetheless, the evidence for adaptive evolution of mutational robustness in biological populations is controversial. Robustness should be selectively favored when mutation rates are high, a common feature of RNA viruses. However, selection for robustness may be relaxed under virus co-infection because complementation between virus genotypes can buffer mutational effects. We therefore hypothesized that selection for genetic robustness in viruses will be weakened with increasing frequency of co-infection. To test this idea, we used populations of RNA phage φ6 that were experimentally evolved at low and high levels of co-infection and subjected lineages of these viruses to mutation accumulation through population bottlenecking. The data demonstrate that viruses evolved under high co-infection show relatively greater mean magnitude and variance in the fitness changes generated by addition of random mutations, confirming our hypothesis that they experience weakened selection for robustness. Our study further suggests that co-infection of host cells may be advantageous to RNA viruses only in the short term. In addition, we observed higher mutation frequencies in the more robust viruses, indicating that evolution of robustness might foster less-accurate genome replication in RNA viruses.

RNA phage viruses evolved under high co-infection of host cells are less robust to mutations than those propagated under low co-infection, suggesting that co-infection may be advantageous to RNA viruses only in the short term.

Evolution of genetic potential

Organisms employ a multitude of strategies to cope with the dynamical environments in which they live. Homeostasis and physiological plasticity buffer changes within the lifetime of an organism, while stochastic developmental programs and hypermutability track changes on longer time-scales. An alternative long-term mechanism is "genetic potential"--a heightened sensitivity to the effects of mutation that facilitates rapid evolution to novel states. Using a transparent mathematical model, we illustrate the concept of genetic potential and show that as environmental variability decreases, the evolving population reaches three distinct steady state conditions: (1) organismal flexibility, (2) genetic potential, and (3) genetic robustness. As a specific example of this concept we examine fluctuating selection for hydrophobicity in a single amino acid. We see the same three stages, suggesting that environmental fluctuations can produce allele distributions that are distinct not only from those found under constant conditions, but also from the transient allele distributions that arise under isolated selective sweeps.

The population genetic theory of hidden variation and genetic robustness

One of the most solid generalizations of transmission genetics is that the phenotypic variance of populations carrying a major mutation is increased relative to the wild type. At least some part of this higher variance is genetic and due to release of previously hidden variation. Similarly, stressful environments also lead to the expression of hidden variation. These two observations have been considered as evidence that the wild type has evolved robustness against genetic variation, i.e., genetic canalization. In this article we present a general model for the interaction of a major mutation or a novel environment with the additive genetic basis of a quantitative character under stabilizing selection. We introduce an approximation to the genetic variance in mutation-selection-drift balance that includes the previously used stochastic Gaussian and house-of-cards approximations as limiting cases. We then show that the release of hidden genetic variation is a generic property of models with epistasis or genotype-environment interaction, regardless of whether the wild-type genotype is canalized or not. As a consequence, the additive genetic variance increases upon a change in the environment or the genetic background even if the mutant character state is as robust as the wild-type character. Estimates show that this predicted increase can be considerable, in particular in large populations and if there are conditionally neutral alleles at the loci underlying the trait. A brief review of the relevant literature suggests that the assumptions of this model are likely to be generic for polygenic traits. We conclude that the release of hidden genetic variance due to a major mutation or environmental stress does not demonstrate canalization of the wild-type genotype.

Evolutionary capacitance may be favored by natural selection

Evolutionary capacitors phenotypically reveal a stock of cryptic genetic variation in a reversible fashion. The sudden and reversible revelation of a range of variation is fundamentally different from the gradual introduction of variation by mutation. Here I study the invasion dynamics of modifiers of revelation. A modifier with the optimal rate of revelation mopt has a higher probability of invading any other population than of being counterinvaded. mopt varies with the population size N and the rate theta at which environmental change makes revelation adaptive. For small populations less than a minimum cutoff Nmin, all revelation is selected against. Nmin is typically quite small and increases only weakly, with theta-1/2. For large populations with N>1/theta, mopt is approximately 1/N. Selection for the optimum is highly effective and increases in effectiveness with larger N>>1/theta. For intermediate values of N, mopt is typically a little less than theta and is only weakly favored over less frequent revelation. The model is analogous to a two-locus model for the evolution of a mutator allele. It is a fully stochastic model and so is able to show that selection for revelation can be strong enough to overcome random drift.

Evolution of dominance in metabolic pathways

Dominance is a form of phenotypic robustness to mutations. Understanding how such robustness can evolve provides a window into how the relation between genotype and phenotype can evolve. As such, the issue of dominance evolution is a question about the evolution of inheritance systems. Attempts at explaining the evolution of dominance have run into two problems. One is that selection for dominance is sensitive to the frequency of heterozygotes. Accordingly, dominance cannot evolve unless special conditions lead to the presence of a high frequency of mutant alleles in the population. Second, on the basis of theoretical results in metabolic control analysis, it has been proposed that metabolic systems possess inherent constraints. These hypothetical constraints imply the default manifestation of dominance of the wild type with respect to the effects of mutations at most loci. Hence, some biologists have maintained that an evolutionary explanation is not relevant to dominance. In this article, we put into question the hypothetical assumption of default metabolic constraints. We show that this assumption is based on an exclusion of important nonlinear interactions that can occur between enzymes in a pathway. With an a priori exclusion of such interactions, the possibility of epistasis and hence dominance modification is eliminated. We present a theoretical model that integrates enzyme kinetics and population genetics to address dominance evolution in metabolic pathways. In the case of mutations that decrease enzyme concentrations, and given the mechanistic constraints of Michaelis-Menten-type catalysis, it is shown that dominance of the wild type can be extensively modified in a two-enzyme pathway. Moreover, we discuss analytical results indicating that the conclusions from the two-enzyme case can be generalized to any number of enzymes. Dominance modification is achieved chiefly through changes in enzyme concentrations or kinetic parameters such as k(cat), both of which can alter saturation levels. Low saturation translates into higher levels of dominance with respect to mutations that decrease enzyme concentrations. Furthermore, it is shown that in the two-enzyme example, dominance evolves as a by-product of selection in a manner that is insensitive to the frequency of heterozygotes. Using variation in k(cat) as an example of modifier mutations, it is shown that the latter can have direct fitness effects in addition to dominance modification effects. Dominance evolution can occur in a frequency-insensitive manner as a result of selection for such dual-effects alleles. This type of selection may prove to be a common pattern for the evolution of phenotypic robustness to mutations.

Molecular Population Genetics of Redundant Floral-Regulatory Genes in Arabidopsis thaliana

Functional redundancy between duplicated genes is predicted to be transitory, as one gene either loses its function or gains a new function, or both genes accrue degenerative, yet complimentary mutations. Yet there are many examples where functional redundancy has been maintained between gene duplicates. To determine whether selection is acting on functionally redundant gene duplicates, we performed molecular evolution and population genetic analyses between two pairs of functionally redundant MADS-box genes from the model plant Arabidopsis thaliana: SEPALLATA1 (SEP1) and SEPALLATA2 (SEP2), involved in floral organ identity, and SHATTERPROOF1 (SHP1) and SHATTERPROOF2 (SHP2), involved in seed shattering. We found evidence for purifying selection acting to constrain functional divergence between paralogous genes. The protein evolution of both pairs of duplicate genes is functionally constrained, as evidenced by Ka/Ks ratios of 0.16 between paralogs. This functional constraint is stronger in the highly conserved DNA-binding and protein-binding MIK region than in the C-terminal region. We also assayed the evolutionary forces acting between orthologs of the SEP and SHP genes in A. thaliana and the closely related species, Arabidopsis lyrata. Heterogeneity analyses of the polymorphism-to-divergence ratio indicate selective sweeps have occurred within the transcriptional unit of SHP1 and the promoter of SHP2 in the A. thaliana lineage. Similar analyses identified a significant reduction in polymorphism within the SEP1 locus, spanning the 3' region of intron 1 to exon 3, that may represent an intragenic sweep within the SEP1 locus. We discuss whether the evolutionary forces acting on SEP1 and SEP2 versus SHP1 and SHP2 vary according to their position in the floral developmental pathway, as found with other floral-regulatory genes.

Genetic variation in cancer predisposition: mutational decay of a robust genetic control network

A computational model of cancer progression is used to study how mutations in genes that control tumor initiation and progression accumulate in populations. The model assumes that cancer occurs only after a cell lineage has progressed through a series of stages. The greater the number of stages, the more strongly the individual is protected against cancer. It is shown that an extra stage initially improves the survival of individuals by decreasing mortality from cancer. However, the additional buffering by an extra stage reduces the impact of any single hereditary mutation and therefore allows the accumulation of more nonlethal mutations in the population. Extra stages thereby lead to the evolution of partially decreased cancer mortality and significantly increased genetic predisposition to disease in the population as a whole. In general, the model illustrates how all robust control networks allow the accumulation of deleterious mutations. An increase in the number of buffering components leads to significant mutational decay in the protection provided by each buffering component and increased genetic predisposition to disease. An extra buffering component's net contribution to survival and reproduction is often small.