Twenty years on from Developmental Plasticity and Evolution: middle-range theories and how to test them

In Developmental Plasticity and Evolution, Mary-Jane West-Eberhard argued that the developmental mechanisms that enable organisms to respond to their environment are fundamental causes of adaptation and diversification. Twenty years after publication of this book, this once so highly controversial claim appears to have been assimilated by a wealth of studies on 'plasticity-led' evolution. However, we suggest that the role of development in explanations for adaptive evolution remains underappreciated in this body of work. By combining concepts of evolvability from evolutionary developmental biology and quantitative genetics, we outline a framework that is more appropriate to identify developmental causes of adaptive evolution. This framework demonstrates how experimental and comparative developmental biology and physiology can be leveraged to put the role of plasticity in evolution to the test.

Adaptive diversification and niche packing on rugged fitness landscapes

Explaining the emergence of diversity and the coexistence of competing types has long been one of the main goals of ecological theory. Rugged fitness landscapes have often been used to explain diversity through the presence of local peaks, or adaptive zones, in the fitness landscape acting as available niches for different species. Alternatively, niche-packing and theories based on limiting similarity describe frequency-dependent selection leading to the organic differentiation of a continuous phenotype space into multiple coexisting types. By combining rugged carrying capacity landscapes with frequency-dependent selection, here we investigate the effects of ruggedness on adaptive diversification and stably maintained diversity. We show that while increased ruggedness often leads to a decreased opportunity for adaptive diversification, it is the shape of the global carrying capacity function, not the local ruggedness, that determines the diversity of the ESS and the total diversity a system can stably maintain.

Evolution of enzyme levels in metabolic pathways: A theoretical approach. Part 2

Metabolism is essential for cell function and adaptation. Because of their central role in metabolism, kinetic parameters and enzyme concentrations are under constant selective pressure to adapt the fluxes of the metabolic networks to the needs of the organism. In line with various studies dealing with enzyme evolution, we recently developed a model of the evolution of enzyme concentrations under selection for increased flux, considered as a proxy for fitness (Coton et al., 2022). With this model, taking into account two realistic cellular constraints, competition for resources and co-regulation, we determined the evolutionary equilibria and range of neutral variations of enzyme concentrations. In this article, we expanded this model by considering that the enzymes in a pathway can belong to different co-regulation groups. We determined the equilibria and showed that the constraints modify the adaptive landscape by limiting the number of independent dimensions. We also showed that any trade-off between enzyme concentrations is sufficient to limit the flux and relax selection for increasing the concentration of other enzymes. Even though this model is based on simplifying assumptions, the complexity of the relationship between enzyme concentrations prevents the formal analysis of the range of neutral variation of enzyme concentrations. However, we could show that selection for maximizing the flux results in selective neutrality for all enzymes regardless the constraints applied, giving generality to the prediction of Hartl et al. (1985).

A common mechanism drives the alignment between the micro- and macroevolution of primate molars

A central challenge for biology is to reveal how different levels of biological variation interact and shape diversity. However, recent experimental studies have indicated that prevailing models of evolution cannot readily explain the link between micro- and macroevolution at deep timescales. Here, we suggest that this paradox could be the result of a common mechanism driving a correlated pattern of evolution. We examine the proportionality between genetic variance and patterns of trait evolution in a system whose developmental processes are well understood to gain insight into how such alignment between morphological divergence and genetic variation might be maintained over macroevolutionary time. Primate molars present a model system by which to link developmental processes to evolutionary dynamics because of the biased pattern of variation that results from the developmental architecture regulating their formation. We consider how this biased variation is expressed at the population level, and how it manifests through evolution across primates. There is a strong correspondence between the macroevolutionary rates of primate molar divergence and their genetic variation. This suggests a model of evolution in which selection is closely aligned with the direction of genetic variance, phenotypic variance, and the underlying developmental architecture of anatomical traits.

Evolution of enzyme levels in metabolic pathways: A theoretical approach

The central role of metabolism in cell functioning and adaptation has given rise to countless studies on the evolution of enzyme-coding genes and network topology. However, very few studies have addressed the question of how enzyme concentrations change in response to positive selective pressure on the flux, considered a proxy of fitness. In particular, the way cellular constraints, such as resource limitations and co-regulation, affect the adaptive landscape of a pathway under selection has never been analyzed theoretically. To fill this gap, we developed a model of the evolution of enzyme concentrations that combines metabolic control theory and an adaptive dynamics approach, and integrates possible dependencies between enzyme concentrations. We determined the evolutionary equilibria of enzyme concentrations and their range of neutral variation, and showed that they differ with the properties of the enzymes, the constraints applied to the system and the initial enzyme concentrations. Simulations of long-term evolution confirmed all analytical and numerical predictions, even though we relaxed the simplifying assumptions used in the analytical treatment.

Trophic niche shifts and phenotypic trait evolution are largely decoupled in Australasian parrots

BACKGROUND

Trophic shifts from one dietary niche to another have played major roles in reshaping the evolutionary trajectories of a wide range of vertebrate groups, yet their consequences for morphological disparity and species diversity differ among groups.

METHODS

Here, we use phylogenetic comparative methods to examine whether the evolution of nectarivory and other trophic shifts have driven predictable evolutionary pathways in Australasian psittaculid parrots in terms of ecological traits such as body size, beak shape, and dispersal capacity.

RESULTS

We found no evidence for an 'early-burst' scenario of lineage or morphological diversification. The best-fitting models indicate that trait evolution in this group is characterized by abrupt phenotypic shifts (evolutionary jumps), with no sign of multiple phenotypic optima correlating with different trophic strategies. Thus, our results point to the existence of weak directional selection and suggest that lineages may be evolving randomly or slowly toward adaptive peaks they have not yet reached.

CONCLUSIONS

This study adds to a growing body of evidence indicating that the relationship between avian morphology and feeding ecology may be more complex than usually assumed and highlights the importance of adding more flexible models to the macroevolutionary toolbox.

Multidimensional divergent selection, local adaptation, and speciation

Divergent selection applied to one or more traits drives local adaptation and may lead to ecological speciation. Divergent selection on many traits might be termed "multidimensional" divergent selection. There is a commonly held view that multidimensional divergent selection is likely to promote local adaptation and speciation to a greater extent than unidimensional divergent selection. We disentangle the core concepts underlying dimensionality as a property of the environment, phenotypes, and genome. In particular, we identify a need to separate the overall strength of selection and the number of loci affected from dimensionality per se, and to distinguish divergence dimensionality from dimensionality of stabilizing selection. We then critically scrutinize this commonly held view that multidimensional selection promotes speciation, re-examining the evidence base from theory, experiments, and nature. We conclude that the evidence base is currently weak and generally suffers from confounding of possible causal effects. Finally, we propose several mechanisms by which multidimensional divergent selection and related processes might influence divergence, both as a driver and as a barrier.

Exploring adaptive landscapes across deep time: A case study using echinoid body size

Adaptive landscapes are a common way of conceptualizing the phenotypic evolution of lineages across deep time. Although multiple approaches exist to implement this concept into operational models of trait evolution, inferring adaptive landscapes from comparative datasets remains challenging. Here, I explore the macroevolutionary dynamics of echinoid body size using data from over 5000 specimens and a phylogenetic framework incorporating a dense fossil sampling and spanning approximately 270 million years. Furthermore, I implement a novel approach of exploring alternative parameterizations of adaptive landscapes that succeeds in finding simpler, yet better-fitting models. Echinoid body size has been constrained to evolve within a single adaptive optimum for much of the clade's history. However, most of the morphological disparity of echinoids was generated by multiple regime shifts that drove the repeated evolution of miniaturized and gigantic forms. Events of body size innovation occurred predominantly in the Late Cretaceous and were followed by a drastic slowdown following the Cretaceous-Paleogene mass extinction. The discovery of these patterns is contingent upon directly sampling fossil taxa. The macroevolution of echinoid body size is therefore characterized by a late increase in disparity (likely linked to an expansion of ecospace), generated by active processes driving lineages toward extreme morphologies.

Loss of consumers constrains phenotypic evolution in the resulting food web

The loss of biodiversity is altering the structure of ecological networks; however, we are currently in a poor position to predict how these altered communities will affect the evolution of remaining populations. Theory on fitness landscapes provides a framework for predicting how selection alters the evolutionary trajectory and adaptive potential of populations, but often treats the network of interacting populations as a “black box.” Here, we integrate ecological networks and fitness landscapes to examine how changes in food‐web structure shape phenotypic evolution. We conducted a field experiment that removed a guild of larval parasitoids that imposed direct and indirect selection pressures on an insect herbivore. We then measured herbivore survival as a function of three key phenotypic traits to estimate directional, quadratic, and correlational selection gradients in each treatment. We used these selection gradients to characterize the slope and curvature of the fitness landscape to understand the direct and indirect effects of consumer loss on phenotypic evolution. We found that the number of traits under directional selection increased with the removal of larval parasitoids, indicating evolution was more constrained toward a specific combination of traits. Similarly, we found that the removal of larval parasitoids altered the curvature of the fitness landscape in such a way that tended to decrease the evolvability of the traits we measured in the next generation. Our results suggest that the loss of trophic interactions can impose greater constraints on phenotypic evolution. This indicates that the simplification of ecological communities may constrain the adaptive potential of remaining populations to future environmental change.

Testing eco-evolutionary predictions using fossil data: Phyletic evolution following ecological opportunity

Fossil sequences provide observations of phenotypes within a lineage over time and represent essential data for increasing our understanding of phyletic evolution beyond microevolutionary timescales. I investigate if fossil time series of the diatom Stephanodiscus niagarae/yellowstonensis follow evolutionary dynamics compatible with hypotheses for how the adaptive landscape changes when a population enters a new environment. The lineage-which has a remarkably detailed stratigraphic record-invaded Yellowstone Lake immediately after recession of ice from the basin 14,000 years ago. Several phyletic models portraying different types of evolutionary dynamics-both compatible and not compatible with changes in the adaptive landscape following ecological opportunity-were fitted to the fossil times-series of S. niagarae/yellowstonensis. Different models best describe the three analyzed traits. Two of the models (a new model of decelerated evolution and an Ornstein-Uhlenbeck model) capture trait dynamics compatible with an event of ecological opportunity, whereas the third model (random walk) does not. Entering a new environment may accordingly affect trait dynamics for thousands of years, but the effects can vary across phenotypes. However, tests of model adequacy reveal shortcomings in all three models explaining the trait dynamics, suggesting model development is needed to more fully understand the phyletic evolution in S. niagarae/yellowstonensis.

The multi-peak adaptive landscape of crocodylomorph body size evolution

BACKGROUND

Little is known about the long-term patterns of body size evolution in Crocodylomorpha, the > 200-million-year-old group that includes living crocodylians and their extinct relatives. Extant crocodylians are mostly large-bodied (3-7 m) predators. However, extinct crocodylomorphs exhibit a wider range of phenotypes, and many of the earliest taxa were much smaller (< 1.2 m). This suggests a pattern of size increase through time that could be caused by multi-lineage evolutionary trends of size increase or by selective extinction of small-bodied species. Here, we characterise patterns of crocodylomorph body size evolution using a model fitting-approach (with cranial measurements serving as proxies). We also estimate body size disparity through time and quantitatively test hypotheses of biotic and abiotic factors as potential drivers of crocodylomorph body size evolution.

RESULTS

Crocodylomorphs reached an early peak in body size disparity during the Late Jurassic, and underwent an essentially continual decline since then. A multi-peak Ornstein-Uhlenbeck model outperforms all other evolutionary models fitted to our data (including both uniform and non-uniform), indicating that the macroevolutionary dynamics of crocodylomorph body size are better described within the concept of an adaptive landscape, with most body size variation emerging after shifts to new macroevolutionary regimes (analogous to adaptive zones). We did not find support for a consistent evolutionary trend towards larger sizes among lineages (i.e., Cope's rule), or strong correlations of body size with climate. Instead, the intermediate to large body sizes of some crocodylomorphs are better explained by group-specific adaptations. In particular, the evolution of a more aquatic lifestyle (especially marine) correlates with increases in average body size, though not without exceptions.

CONCLUSIONS

Shifts between macroevolutionary regimes provide a better explanation of crocodylomorph body size evolution on large phylogenetic and temporal scales, suggesting a central role for lineage-specific adaptations rather than climatic forcing. Shifts leading to larger body sizes occurred in most aquatic and semi-aquatic groups. This, combined with extinctions of groups occupying smaller body size regimes (particularly during the Late Cretaceous and Cenozoic), gave rise to the upward-shifted body size distribution of extant crocodylomorphs compared to their smaller-bodied terrestrial ancestors.

Feather evolution exemplifies sexually selected bridges across the adaptive landscape

Over the last two decades, paleontologists have pieced together the early evolutionary history of feathers. Simple hair-like feathers served as insulating pelage, but the first feathers with complex branching structures and a plainer form evolved for the purpose of sexual display. The evolution of these complex display feathers was essential to the later evolution of flight. Feathers illustrate how sexual selection can generate complex novel phenotypes, which are then available for natural selection to modify and direct toward novel functions. In the longstanding metaphor of the adaptive landscape, sexual selection is a means by which lineages resting on one adaptive peak may gradually bridge a gap to another peak, without the landscape itself being first altered by environmental changes.

Performance in three shell functions predicts the phenotypic distribution of hard-shelled turtles

Adaptive landscapes have served as fruitful guides to evolutionary research for nearly a century. Current methods guided by landscape frameworks mostly utilize evolutionary modeling (e.g., fitting data to Ornstein-Uhlenbeck models) to make inferences about adaptive peaks. Recent alternative methods utilize known relationships between phenotypes and functional performance to derive information about adaptive landscapes; this information can then help explain the distribution of species in phenotypic space and help infer the relative importance of various functions for guiding diversification. Here, data on performance for three turtle shell functions-strength, hydrodynamic efficiency, and self-righting ability-are used to develop a set of predicted performance optima in shell shape space. The distribution of performance optima shows significant similarity to the distribution of existing turtle species and helps explain the absence of shells in otherwise anomalously empty regions of morphospace. The method outperforms a modeling-based approach in inferring the location of reasonable adaptive peaks and in explaining the shape of the phenotypic distributions of turtle shells. Performance surface-based methods allow researchers to more directly connect functional performance with macroevolutionary diversification, and to explain the distribution of species (including presences and absences) across phenotypic space.

Endogenous Molecular-Cellular Network Cancer Theory: A Systems Biology Approach

In light of ever apparent limitation of the current dominant cancer mutation theory, a quantitative hypothesis for cancer genesis and progression, endogenous molecular-cellular network hypothesis has been proposed from the systems biology perspective, now for more than 10 years. It was intended to include both the genetic and epigenetic causes to understand cancer. Its development enters the stage of meaningful interaction with experimental and clinical data and the limitation of the traditional cancer mutation theory becomes more evident. Under this endogenous network hypothesis, we established a core working network of hepatocellular carcinoma (HCC) according to the hypothesis and quantified the working network by a nonlinear dynamical system. We showed that the two stable states of the working network reproduce the main known features of normal liver and HCC at both the modular and molecular levels. Using endogenous network hypothesis and validated working network, we explored genetic mutation pattern in cancer and potential strategies to cure or relieve HCC from a totally new perspective. Patterns of genetic mutations have been traditionally analyzed by posteriori statistical association approaches in light of traditional cancer mutation theory. One may wonder the possibility of a priori determination of any mutation regularity. Here, we found that based on the endogenous network theory the features of genetic mutations in cancers may be predicted without any prior knowledge of mutation propensities. Normal hepatocyte and cancerous hepatocyte stable states, specified by distinct patterns of expressions or activities of proteins in the network, provide means to directly identify a set of most probable genetic mutations and their effects in HCC. As the key proteins and main interactions in the network are conserved through cell types in an organism, similar mutational features may also be found in other cancers. This analysis yielded straightforward and testable predictions on an accumulated and preferred mutation spectrum in normal tissue. The validation of predicted cancer state mutation patterns demonstrates the usefulness and potential of a causal dynamical framework to understand and predict genetic mutations in cancer. We also obtained the following implication related to HCC therapy, (1) specific positive feedback loops are responsible for the maintenance of normal liver and HCC; (2) inhibiting proliferation and inflammation-related positive feedback loops, and simultaneously inducing liver-specific positive feedback loop is predicated as the potential strategy to cure or relieve HCC; (3) the genesis and regression of HCC is asymmetric. In light of the characteristic property of the nonlinear dynamical system, we demonstrate that positive feedback loops must be existed as a simple and general molecular basis for the maintenance of phenotypes such as normal liver and HCC, and regulating the positive feedback loops directly or indirectly provides potential strategies to cure or relieve HCC.

With a little help from my friends: cooperation can accelerate the rate of adaptive valley crossing

Background

Natural selection favors changes that lead to genotypes possessing high fitness. A conflict arises when several mutations are required for adaptation, but each mutation is separately deleterious. The process of a population evolving from a genotype encoding for a local fitness maximum to a higher fitness genotype is termed an adaptive peak shift.

Results

Here we suggest cooperative behavior as a factor that can facilitate adaptive peak shifts. We model cooperation in a public goods scenario, wherein each individual contributes resources that are later equally redistributed among all cooperating individuals. We use mathematical modeling and stochastic simulations to study the effect of cooperation on peak shifts in both panmictic and structured populations. Our results show that cooperation can substantially affect the rate of complex adaptation. Furthermore, we show that cooperation increases the population diversity throughout the peak shift process, thus increasing the robustness of the population to sudden environmental changes.

Conclusions

We provide a new explanation to adaptive valley crossing in natural populations and suggest that the long term evolution of a species depends on its social behavior.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0983-2) contains supplementary material, which is available to authorized users.

Exploring the hidden landscape of female preferences for complex signals

A major challenge in evolutionary biology is explaining the origins of complex phenotypic diversity. In animal communication, complex signals may evolve from simpler signals because novel signal elements exploit preexisting biases in receivers' sensory systems. Investigating the shape of female preference functions for novel signal characteristics is a powerful, but underutilized, method to describe the adaptive landscape potentially guiding complex signal evolution. We measured female preference functions for characteristics of acoustic appendages added to male calling songs in the grasshopper Chorthippus biguttulus, which naturally produces only simple songs. We discovered both hidden preferences for and biases against novel complex songs, and identified rules governing song attractiveness based on multiple characteristics of both the base song and appendage. The appendage's temporal position and duration were especially important: long appendages preceding the song often made songs less attractive, while following appendages were neutral or weakly attractive. Appendages had stronger effects on songs of shorter duration, but did not restore the attractiveness of very unattractive songs. We conclude that sensory biases favor, within predictable limits, the evolution of complex songs in grasshoppers. The function-valued approach is an important tool in determining the generality of these limits in other taxa and signaling modalities.

Adaptive multiscapes: an up-to-date metaphor to visualize molecular adaptation

Background

Wright’s metaphor of the fitness landscape has shaped and conditioned our view of the adaptation of populations for almost a century. Since its inception, and including criticism raised by Wright himself, the concept has been surrounded by controversy. Among others, the debate stems from the intrinsic difficulty to capture important features of the space of genotypes, such as its high dimensionality or the existence of abundant ridges, in a visually appealing two-dimensional picture. Two additional currently widespread observations come to further constrain the applicability of the original metaphor: the very skewed distribution of phenotype sizes (which may actively prevent, due to entropic effects, the achievement of fitness maxima), and functional promiscuity (i.e. the existence of secondary functions which entail partial adaptation to environments never encountered before by the population).

Results

Here we revise some of the shortcomings of the fitness landscape metaphor and propose a new “scape” formed by interconnected layers, each layer containing the phenotypes viable in a given environment. Different phenotypes within a layer are accessible through mutations with selective value, while neutral mutations cause displacements of populations within a phenotype. A different environment is represented as a separated layer, where phenotypes may have new fitness values, other phenotypes may be viable, and the same genotype may yield a different phenotype, representing genotypic promiscuity. This scenario explicitly includes the many-to-many structure of the genotype-to-phenotype map. A number of empirical observations regarding the adaptation of populations in the light of adaptive multiscapes are reviewed.

Conclusions

Several shortcomings of Wright’s visualization of fitness landscapes can be overcome through adaptive multiscapes. Relevant aspects of population adaptation, such as neutral drift, functional promiscuity or environment-dependent fitness, as well as entropic trapping and the concomitant impossibility to reach fitness peaks are visualized at once. Adaptive multiscapes should aid in the qualitative understanding of the multiple pathways involved in evolutionary dynamics.

Reviewers

This article was reviewed by Eugene Koonin and Ricard Solé.

Context dependence in complex adaptive landscapes: frequency and trait-dependent selection surfaces within an adaptive radiation of Caribbean pupfishes

The adaptive landscape provides the foundational bridge between micro- and macroevolution. One well-known caveat to this perspective is that fitness surfaces depend on ecological context, including competitor frequency, traits measured, and resource abundance. However, this view is based largely on intraspecific studies. It is still unknown how context-dependence affects the larger features of peaks and valleys on the landscape which ultimately drive speciation and adaptive radiation. Here, I explore this question using one of the most complex fitness landscapes measured in the wild in a sympatric pupfish radiation endemic to San Salvador Island, Bahamas by tracking survival and growth of laboratory-reared F2 hybrids. I present new analyses of the effects of competitor frequency, dietary isotopes, and trait subsets on this fitness landscape. Contrary to expectations, decreasing competitor frequency increased survival only among very common phenotypes, whereas less common phenotypes rarely survived despite few competitors, suggesting that performance, not competitor frequency, shapes large-scale features of the fitness landscape. Dietary isotopes were weakly correlated with phenotype and growth, but did not explain additional survival variation. Nonlinear fitness surfaces varied substantially among trait subsets, revealing one-, two-, and three-peak landscapes, demonstrating the complexity of selection in the wild, even among similar functional traits.

Empirically founded genotype-phenotype maps from mammalian cyclic nucleotide-gated ion channels

A major barrier between evolutionary and functional biology is the difficulty of determining appropriate genotype-phenotype-fitness maps, particularly in metazoans. Concrete perspectives towards unifying these approaches are offered by studies on the physiological systems that depend on ion channel dynamics. I focus on the cyclic nucleotide-gated (CNG) channels implicated in the photoreceptor's response to light. From an evolutionary standpoint, sensory systems offers interpretative advantages, as the relation between the sensory response and environment is relatively straightforward. For CNG and other ion channels, extensive data are available about the physiological consequences of scanning mutagenesis on sensitive protein domains, such as the conduction pore. Mutant ion channels can be easily studied in living cells, so that the relation between genotypes and phenotypes is less speculative than usual. By relying on relatively simple theoretical frameworks, I used these data to relate the sequence space with phenotypes at increasing hierarchical levels. These empirical genotype-phenotype and phenotype-phenotype landscapes became smoother at higher integration levels, especially in heterozygous condition. The epistatic interaction between sites was analyzed from double mutant constructs. Magnitude epistasis was common. Moreover, evidence of reciprocal sign epistasis and the presence of permissive mutations were also observed, which suggest how adaptive regions can be connected across maladaptive valleys. The approach I describe suggests a way to better relate the evolutionary dynamics with the underlying physiology.

Coarse dark patterning functionally constrains adaptive shifts from aposematism to crypsis in strawberry poison frogs

Ecological specialization often requires tight coevolution of several traits, which may constrain future evolutionary pathways and make species more prone to extinction. Aposematism and crypsis represent two specialized adaptations to avoid predation. We tested whether the combined effects of color and pattern on prey conspicuousness functionally constrain or facilitate shifts between these two adaptations. We combined data from 17 natural populations of strawberry poison frogs, Oophaga pumilio with an experimental approach using digitalized images of frogs and chickens as predators. We show that bright coloration often co-occurs with coarse patterning among the natural populations. Dull green frogs with coarse patterning are rare in nature but in the experiment they were as easily detected as bright red frogs suggesting that this trait combination represents a transient evolutionary state toward aposematism. Hence, a gain of either bright color or coarse patterning leads to conspicuousness, but a transition back to crypsis would be functionally constrained in populations with both bright color and coarse patterning by requiring simultaneous changes in two traits. Thus, populations (or species) signaling aposematism by conspicuous color should be less likely to face an evolutionary dead end and more likely to radiate than populations with both conspicuous color and coarse patterning.

Phenotypic and genotypic convergences are influenced by historical contingency and environment in yeast

Different organisms have independently and recurrently evolved similar phenotypic traits at different points throughout history. This phenotypic convergence may be caused by genotypic convergence and in addition, constrained by historical contingency. To investigate how convergence may be driven by selection in a particular environment and constrained by history, we analyzed nine life-history traits and four metabolic traits during an experimental evolution of six yeast strains in four different environments. In each of the environments, the population converged toward a different multivariate phenotype. However, the evolution of most traits, including fitness components, was constrained by history. Phenotypic convergence was partly associated with the selection of mutations in genes involved in the same pathway. By further investigating the convergence in one gene, BMH1, mutated in 20% of the evolved populations, we show that both the history and the environment influenced the types of mutations (missense/nonsense), their location within the gene itself, as well as their effects on multiple traits. However, these effects could not be easily predicted from ancestors' phylogeny or past selection. Combined, our data highlight the role of pleiotropy and epistasis in shaping a rugged fitness landscape.

Epistasis, pleiotropy, and the mutation load in sexual and asexual populations

Mutation may impose a substantial load on populations, which varies according to the reproductive mode of organisms. Over the past years, various authors used adaptive landscape models to predict the long-term effect of mutation on mean fitness; however, many of these studies assumed very weak mutation rates, so that at most one mutation segregates in the population. In this article, we derive several simple approximations (confirmed by simulations) for the mutation load at high mutation rate (U), using a general model that allows us to play with the number of selected traits (n), the degree of pleiotropy of mutations, and the shape of the fitness function (which affects the average sign and magnitude of epistasis among mutations). When mutations have strong fitness effects, the equilibrium fitness W¯ of sexuals and asexuals is close to e(-U); under weaker mutational effects, sexuals reach a different regime where W¯ is a simple function of U and of a parameter describing the shape of the fitness function. Contrarily to weak mutation results showing that W¯ is an increasing function of population size and a decreasing function of n, these parameters may have opposite effects in sexual populations at high mutation rate.