The population genetic theory of hidden variation and genetic robustness

Evolution of coadaptation in a subdivided population

The interplay between population subdivision and epistasis is investigated by studying the fixation probability of a coadapted haplotype in a subdivided population. Analytical and simulation models are developed to study the evolutionary fate of two conditionally neutral mutations that interact epistatically to enhance fitness. We find that the fixation probability of a coadapted haplotype shows a marked increase when the population is genetically subdivided and subpopulations are loosely connected by migration. Moderate migration and isolation allow the propagation of the mutant alleles across subpopulations, while at the same time preserving the favorable allelic combination established within each subpopulation. Together they create the condition most favorable for the ultimate fixation of the coadapted haplotype. On the basis of the analytical and simulation results, we discuss the fundamental role of population subdivision and restricted gene flow in promoting the evolution of functionally integrated systems, with some implications for the shifting-balance theory of evolution.

Timescales of genetic and epigenetic inheritance

According to classical evolutionary theory, phenotypic variation originates from random mutations that are independent of selective pressure. However, recent findings suggest that organisms have evolved mechanisms to influence the timing or genomic location of heritable variability. Hypervariable contingency loci and epigenetic switches increase the variability of specific phenotypes; error-prone DNA replicases produce bursts of variability in times of stress. Interestingly, these mechanisms seem to tune the variability of a given phenotype to match the variability of the acting selective pressure. Although these observations do not undermine Darwin's theory, they suggest that selection and variability are less independent than once thought.

The Cost of Keeping Eggs Fresh: Quantitative Genetic Variation in Females that Mate Late Relative to Sexual Maturation

In many species, females abandon mate choice to ensure that eggs are fertilized before they are lost. But why do females not just maintain oocytes longer if there is a benefit to mate choice? We conducted a quantitative genetic study in the cockroach Nauphoeta cinerea to test whether genetic constraints prevent the evolution of oocyte maintenance or selection against oocyte loss is weak when females mate late relative to sexual maturity. We found standing genetic variation within the population and no evidence for genetic constraints. Levels of genetic variation are of the magnitude found for life-history traits in general, suggesting that this trait has been exposed to selection. We unexpectedly found two categories of females: those that delay reproduction and those that reproduce at a normal time when mating late, which could indicate alternative strategies. However, frequency-dependent selection does not maintain this variation as females that delay always reproduce less well. Given these findings, we suggest that there may be advantages to egg degradation. The evolution of maintenance of fertilizable oocytes over time would then be constrained by the need to maintain the mechanism by which females control the distribution of resources between current and future reproductive events.

Mutations leading to loss of sporulation ability in Bacillus subtilis are sufficiently frequent to favor genetic canalization

We measured the rate of mutations impairing sporulation ability in Bacillus subtilis as 0.003 in a mutator population, following 6000 generations of strong selection for sporulation that have previously been described. This means that the product of the population size and the functional mutation rate is approximately 10(5), well within the parameter range for which genetic canalization of sporulation ability is expected.

Control of canalization and evolvability by Hsp90

Partial reduction of Hsp90 increases expression of morphological novelty in qualitative traits of Drosophila and Arabidopsis, but the extent to which the Hsp90 chaperone also controls smaller and more likely adaptive changes in natural quantitative traits has been unclear. To determine the effect of Hsp90 on quantitative trait variability we deconstructed genetic, stochastic and environmental components of variation in Drosophila wing and bristle traits of genetically matched flies, differing only by Hsp90 loss-of-function or wild-type alleles. Unexpectedly, Hsp90 buffering was remarkably specific to certain normally invariant and highly discrete quantitative traits. Like the qualitative trait phenotypes controlled by Hsp90, highly discrete quantitative traits such as scutellor and thoracic bristle number are threshold traits. When tested across genotypes sampled from a wild population or in laboratory strains, the sensitivity of these traits to many types of variation was coordinately controlled, while continuously variable bristle types and wing size, and critically invariant left-right wing asymmetry, remained relatively unaffected. Although increased environmental variation and developmental noise would impede many types of selection response, in replicate populations in which Hsp90 was specifically impaired, heritability and ‘extrinsic evolvability’, the expected response to selection, were also markedly increased. However, despite the overall buffering effect of Hsp90 on variation in populations, for any particular individual or genotype in which Hsp90 was impaired, the size and direction of its effects were unpredictable. The trait and genetic-background dependence of Hsp90 effects and its remarkable bias toward invariant or canalized traits support the idea that traits evolve independent and trait-specific mechanisms of canalization and evolvability through their evolution of non-linearity and thresholds. Highly non-linear responses would buffer variation in Hsp90-dependent signaling over a wide range, while over a narrow range of signaling near trait thresholds become more variable with increasing probability of triggering all-or-none developmental responses.

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.

Wagner's canalization model

Wagner (1996, Does evolutionary plasticity evolve? Evolution 50, 1008-1023.) and Siegal and Bergman, 2002 and Azevedo et al., 2006 have studied a simple model of the evolution of a network of N genes, in order to explain the observed phenomenon that systems evolve to be robust. These authors primarily considered the case N=10 and used simulations to reach their conclusions. Here we investigate this model in more detail, considering systems of different sizes with and without recombination, and with selection for convergence instead of to a specified limit. For the simpler evolutionary model lacking recombination, we analyze the system as a neutral network. This allows us to describe the equilibrium distribution networks within genotype space. Our results show that, given a sufficiently large population size, the qualitative observation that systems evolve to be robust, is itself robust, as it does not depend on the details of the model. In simple terms, robust systems have more viable offspring, so the evolution of robustness is merely selection for increased fecundity, an observation that is well known in the theory of neutral networks.

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.

Response of fluctuating and directional asymmetry to selection on wing shape in Drosophila melanogaster

We tested whether directional selection on an index-based wing character in Drosophila melanogaster affected developmental stability and patterns of directional asymmetry. We selected for both an increase (up selection) and a decrease (down selection) of the index value on the left wing and compared patterns of fluctuating and directional asymmetry in the selection index and other wing traits across selection lines. Changes in fluctuating asymmetry across selection lines were predominantly small, but we observed a tendency for fluctuating asymmetry to decrease in the up-selected lines in both replicates. Because changes in fluctuating asymmetry depended on the direction of selection, and were not related to changes in trait size, these results fail to support existing hypotheses linking directional selection and developmental stability. Selection also produced a pattern of directional asymmetry that was similar in all selected lines whatever the direction of selection. This result may be interpreted as a release of genetic variance in directional asymmetry under selection.

Effect of E(sev) and Su(Raf) Hsp83 mutants and trans-heterozygotes on bristle trait means and variation in Drosophila melanogaster

The Hsp90 protein encoded by the Hsp83 gene is required for the development of many traits in Drosophila. Hsp83 is also thought to play a role in the expression of phenotypic and genetic variability for subsequent selection and evolutionary change. Here we examine the impact of different E(sev) and Su(Raf) Hsp83 mutants on means and phenotypic variances of invariant and variable bristle traits. One of the mutants influenced the normally invariant thoracic bristle number, while none affected invariant scutellar bristle number. E(sev) alleles consistently influenced variable bristle traits while there were fewer effects of the Su(Raf) alleles. For the variable traits, none of the Hsp83 alleles had any effect on phenotypic variance, environmental variance, or developmental stability of the bristle traits. When alleles were combined in trans-heterozygotes, there were both cumulative and complementary effects on thoracic and variable bristle trait numbers, depending on the allelic combination. Overall, the results suggest that Hsp83 mutants do not have detectable effects on the phenotypic or environmental variance of bristle traits and that complementation of E(sev) and Su(Raf) Hsp83 mutants can extend to thoracic bristles as well as previously reported effects on viability. Some allelic combinations lead to more severe effects on variable bristle trait means than do single Hsp83 mutations.

The brachymorph mouse and the developmental-genetic basis for canalization and morphological integration

Although it is well known that many mutations influence phenotypic variability as well as the mean, the underlying mechanisms for variability effects are very poorly understood. The brachymorph (bm) phenotype results from an autosomal recessive mutation in the phosphoadenosine-phosphosulfate synthetase 2 gene (Papps2). A major cranial manifestation is a dramatic reduction in the growth of the chondrocranium which results from undersulfation of glycosaminoglycans (GAGs) in the cartilage matrix. We found that this reduction in the growth of the chondrocranium is associated with an altered pattern of craniofacial shape variation, a significant increase in phenotypic variance and a dramatic increase in morphological integration for craniofacial shape. Both effects are largest in the basicranium. The altered variation pattern indicates that the mutation produces developmental influences on shape that are not present in the wildtype. As the mutation dramatically reduces sulfation of GAGs, we infer that this influence is variation among individuals in the degree of sulfation, or variable expressivity of the mutation. This variation may be because of genetic variation at other loci that influence sulfation, environmental effects, or intrinsic effects. We infer that chondrocranial development exhibits greater sensitivity to variation in the sulfation of chondroitin sulfate when the degree of sulfation is low. At normal levels, sulfation probably contributes minimally to phenotypic variation. This case illustrates canalization in a particular developmental-genetic context.

Loss of dispensable genes is not adaptive in yeast

A substantial share of genes identified in yeast can be deleted without visible phenotypic effects. Current debate concentrates on the possible roles of seemingly dispensable genes. The costs of maintaining unnecessary functions has attracted little attention. The hypothesis of antagonistic pleiotropy postulates that adaptations to different constituents of the environment are likely to interfere with each other, and therefore loss of unnecessary functions is potentially advantageous. We tested an entire collection of nonessential yeast gene deletions in a benign and nutritionally rich environment in which the number of dispensable genes was particularly high. We applied a series of competition experiments that could detect differences in relative fitness of approximately 0.005. No beneficial deletions were found, except perhaps for the deletion of about a dozen genes that slightly improved competitive ability; however, a functional explanation of the fitness advantage is lacking. The paucity of beneficial gene deletions is striking because genetic adaptations to laboratory conditions are regularly observed in yeast. However, it accords with the finding that the gene contents of four species of Saccharomyces are nearly identical, despite up to 20 million years of independent evolution and extensive DNA sequence divergence. Such extreme conservation of functions would be improbable if there were periods of selection promoting the loss of temporarily dispensable genes. The evident cohesion of the yeast genomes may be their evolved feature or an intrinsic property of complex genetic systems.

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.

The Evolutionary Genetics of Canalization

Evolutionary genetics has recently made enormous progress in understanding how genetic variation maps into phenotypic variation. However why some traits are phenotypically invariant despite apparent genetic and environmental changes has remained a major puzzle. In the 1940s, Conrad Hal Waddington coined the concept and term "canalization" to describe the robustness of phenotypes to perturbation; a similar concept was proposed by Waddington's contemporary Ivan Ivanovich Schmalhausen. This paper reviews what has been learned about canalization since Waddington. Canalization implies that a genotype's phenotype remains relatively invariant when individuals of a particular genotype are exposed to different environments (environmental canalization) or when individuals of the same single- or multilocus genotype differ in their genetic background (genetic canalization). Consequently, genetic canalization can be viewed as a particular kind of epistasis, and environmental canalization and phenotypic plasticity are two aspects of the same phenomenon. Canalization results in the accumulation of phenotypically cryptic genetic variation, which can be released after a "decanalizing" event. Thus, canalized genotypes maintain a cryptic potential for expressing particular phenotypes, which are only uncovered under particular decanalizing environmental or genetic conditions. Selection may then act on this newly released genetic variation. The accumulation of cryptic genetic variation by canalization may therefore increase evolvability at the population level by leading to phenotypic diversification under decanalizing conditions. On the other hand, under canalizing conditions, a major part of the segregating genetic variation may remain phenotypically cryptic; canalization may therefore, at least temporarily, constrain phenotypic evolution. Mechanistically, canalization can be understood in terms of transmission patterns, such as epistasis, pleiotropy, and genotype by environment interactions, and in terms of genetic redundancy, modularity, and emergent properties of gene networks and biochemical pathways. While different forms of selection can favor canalization, the requirements for its evolution are typically rather restrictive. Although there are several methods to detect canalization, there are still serious problems with unambiguously demonstrating canalization, particularly its adaptive value.

The role of epistatic gene interactions in the response to selection and the evolution of evolvability

It has been argued that the architecture of the genotype-phenotype map determines evolvability, but few studies have attempted to quantify these effects. In this article we use the multilinear epistatic model to study the effects of different forms of epistasis on the response to directional selection. We derive an analytical prediction for the change in the additive genetic variance, and use individual-based simulations to understand the dynamics of evolvability and the evolution of genetic architecture. This shows that the major determinant for the evolution of the additive variance, and thus the evolvability, is directional epistasis. Positive directional epistasis leads to an acceleration of evolvability, while negative directional epistasis leads to canalization. In contrast, pure non-directional epistasis has little effect on the response to selection. One consequence of this is that the classical epistatic variance components, which do not distinguish directional and non-directional effects, are useless as predictors of evolutionary dynamics. The build-up of linkage disequilibrium also has negligible effects. We argue that directional epistasis is likely to have major effects on evolutionary dynamics and should be the focus of empirical studies of epistasis.

Stress-induced variation in evolution: from behavioural plasticity to genetic assimilation

Extreme environments are closely associated with phenotypic evolution, yet the mechanisms behind this relationship are poorly understood. Several themes and approaches in recent studies significantly further our understanding of the importance that stress-induced variation plays in evolution. First, stressful environments modify (and often reduce) the integration of neuroendocrinological, morphological and behavioural regulatory systems. Second, such reduced integration and subsequent accommodation of stress-induced variation by developmental systems enables organismal 'memory' of a stressful event as well as phenotypic and genetic assimilation of the response to a stressor. Third, in complex functional systems, a stress-induced increase in phenotypic and genetic variance is often directional, channelled by existing ontogenetic pathways. This accounts for similarity among individuals in stress-induced changes and thus significantly facilitates the rate of adaptive evolution. Fourth, accumulation of phenotypically neutral genetic variation might be a common property of locally adapted and complex organismal systems, and extreme environments facilitate the phenotypic expression of this variance. Finally, stress-induced effects and stress-resistance strategies often persist for several generations through maternal, ecological and cultural inheritance. These transgenerational effects, along with both the complexity of developmental systems and stressor recurrence, might facilitate genetic assimilation of stress-induced effects. Accumulation of phenotypically neutral genetic variance by developmental systems and phenotypic accommodation of stress-induced effects, together with the inheritance of stress-induced modifications, ensure the evolutionary persistence of stress-response strategies and provide a link between individual adaptability and evolutionary adaptation.

Evidence for canalization of Distal-less function in the leg of Drosophila melanogaster

A considerable body of theory pertaining to the evolution of canalization has emerged recently, yet there have been few empirical investigations of their predictions. To address this, patterns of canalization and trait correlation were investigated under the individual and joint effects of the introgression of a loss-of-function allele of the Distal-less gene and high-temperature stress on a panel of iso-female lines. Variation was examined for number of sex comb teeth and the length of the basi-tarsus on the pro-thoracic leg of male Drosophila melanogaster. I demonstrate that whereas there is evidence for trait canalization, there is no evidence to support the hypothesis of the evolution of genetic canalization as a response to microenvironmental canalization. Furthermore, I demonstrate that although there are genetic correlations between these traits, there is no association between their measures of canalization. I discuss the prospects of the evolutionary lability of the Distal-less gene within the context of changes in genetic variation and covariation.

Flexibility in a gene network affecting a simple behavior in Drosophila melanogaster

Gene interactions are emerging as central to understanding the realization of any phenotype. To probe the flexibility of interactions in a defined gene network, we isolated a set of 16 interacting genes in Drosophila, on the basis of their alteration of a quantitative behavioral phenotype-the loss of coordination in a temperature-sensitive allele of Syntaxin1A. The interactions inter se of this set of genes were then assayed in the presence and in the absence of the original Syntaxin1A mutation to ask whether the relationships among the 16 genes remain stable or differ after a change in genetic context. The pattern of epistatic interactions that occurs within this set of variants is dramatically altered in the two different genetic contexts. The results imply considerable flexibility in the network interactions of genes.