Environmental change exposes beneficial epistatic interactions in a catalytic RNA

Natural selection drives populations of individuals towards local peaks in a fitness landscape. These peaks are created by the interactions between individual mutations. Fitness landscapes may change as an environment changes. In a previous contribution, we discovered a variant of the Azoarcus group I ribozyme that represents a local peak in the RNA fitness landscape. The genotype at this peak is distinguished from the wild-type by four point mutations. We here report ribozyme fitness data derived from constructing all possible combinations of these point mutations. We find that these mutations interact epistatically. Importantly, we show that these epistatic interactions change qualitatively in the three different environments that we studied. We find examples where the relative fitness of a ribozyme can change from neutral or negative in one environment, to positive in another. We also show that the fitness effect of a specific GC-AU base pair switch is dependent on both the environment and the genetic context. Moreover, the mutations that we study improve activity at the cost of decreased structural stability. Environmental change is ubiquitous in nature. Our results suggest that such change can facilitate adaptive evolution by exposing new peaks of a fitness landscape. They highlight a prominent role for genotype-environment interactions in doing so.     

Landscape configuration determines gene flow and phenotype in a flightless forest-edge ground-dwelling bush-cricket, Pholidoptera griseoaptera

Spatial configuration of habitats influences genetic structure and population fitness whereas it affects mainly species with limited dispersal ability. To reveal how habitat fragmentation determines dispersal and dispersal-related morphology in a ground-dispersing insect species we used a bush-cricket (Pholidoptera griseoaptera) which is associated with forest-edge habitat. We analysed spatial genetic patterns together with variability of the phenotype in two forested landscapes with different levels of fragmentation. While spatial configuration of forest habitats did not negatively affect genetic characteristics related to the fitness of sampled populations, genetic differentiation was found higher among populations from an extensive forest. Compared to an agricultural matrix between forest patches, the matrix of extensive forest had lower permeability and posed barriers for the dispersal of this species. Landscape configuration significantly affected also morphological traits that are supposed to account for species dispersal potential; individuals from fragmented forest patches had longer hind femurs and a higher femur to pronotum ratio. This result suggests that selection pressure act differently on populations from both landscape types since dispersal-related morphology was related to the level of habitat fragmentation. Thus observed patterns may be explained as plastic according to the level of landscape configuration; while anthropogenic fragmentation of habitats for this species can lead to homogenization of spatial genetic structure.     

The evolution of life histories in spatially heterogeneous environments: Optimal reaction norms revisited

Summary Natural populations live in heterogeneous environments, where habitat variation drives the evolution of phenotypic plasticity. The key feature of population structure addressed in this paper is the net flow of individuals from source (good) to sink (poor) habitats. These movements make it necessary to calculate fitness across the full range of habitats encountered by the population, rather than independently for each habitat. As a consequence, the optimal phenotype in a given habitat not only depends on conditions there but is linked to the performance of individuals in other habitats. We generalize the Euler-Lotka equation to define fitness in a spatially heterogeneous environment in which individuals disperse among habitats as newborn and then stay in a given habitat for life. In this case, maximizing fitness (the rate of increase over all habitats) is equivalent to maximizing the reproductive value of newborn in each habitat but not to maximizing the rate of increase that would result if individuals in each habitat were an isolated population. The new equation can be used to find optimal reaction norms for life history traits, and examples are calculated for age at maturity and clutch size. In contrast to previous results, the optimal reaction norm differs from the line connecting local adaptations of isolated populations each living in only one habitat. Selection pressure is higher in good and frequent habitats than in poor and rare ones. A formula for the relative importance of these two factors allows predictions of the habitat in which the genetic variance about the optimal reaction norm should be smallest.     

Ginkgo: spatially-explicit simulator of complex phylogeographic histories

We present Ginkgo, a software package for agent-based, forward-time simulations of genealogies of multiple unlinked loci from diploid populations. Ginkgo simulates the evolution of one or more species on a spatially explicit landscape of cells. The user of the software can specify the geographical and environmental characteristics of the landscape, and these properties can change according to a prespecified schedule. The geographical elements modelled include the arrangement of cells and movement rates between particular cells. Each species has a function that can calculate a fitness score for any combination of an individual organism's phenotype and environmental characteristics. The user can control the number of fitness factors (the dimensionality of the cell-specific fitness factors and the individuals phenotypic vectors) and the weighting of each of these dimensions in the fitness calculation. Cell-specific fitness trait optima can be specified across the landscape to mimic differences in habitat. In addition to their differing fitness functions, species can differ in terms of their vagility and fecundity. Genealogies and occurrence data can be produced at any time during the simulation in NEXUS and ESRI Ascii Grid formats, respectively.     

Frequency-dependent reproductive success in female common lizards: a real-life hawk–dove–bully game?

Alternative strategies are characterised by context-dependent fitness payoffs, which means that their fitness depends on the frequency and the nature of their interactions with one or more strategies. The analysis of the variation of the fitness of each strategy in different social environments can elucidate the evolutionary dynamics played by the strategies. In the common lizard, three female colour types (yellow, orange and mixed) are associated with alternative reaction norms in reproduction and social behaviour that signal alternative strategies. To clarify the nature of colour-specific interactions and their influence on female fitness, we analysed the response of female reproductive success to an experimental manipulation of colour frequencies in natural populations. We found that juvenile body condition at birth for all colour types was negatively affected by the local frequency of yellow females. In addition, we found that mixed females had higher clutch hatching success in the populations where orange females were frequent. These results prove that female reproduction is sensitive to the social environment, and are consistent with a scenario of a hawk–dove–bully game, in which yellow females are aggressive hawks, orange females non-aggressive doves, and mixed females have a context-dependent bully strategy. In this system, the plastic bully strategy would confer a reproductive advantage to putative heterozygotes in some social environments, which could allow the maintenance of the system through context-dependent overdominance effects.     

GENERALISTS,SPECIALISTS, AND THE EVOLUTION OF PHENOTYPIC PLASTICITY IN SYMPATRIC POPULATIONS OF DISTINCT SPECIES

The evolution of phenotypic plasticity is studied in a model with two reproductively isolated “species” in a coarse-grained environment, consisting of two types of habitats. A quantitative genetic model for selection was constructed, in which habitats differ in the optimal value for a focal trait, and with random dispersal among habitats. The main interest was to study the effects of different selection regimes. Three cases were investigated: (1) without any limits to plasticity; (2) without genetic variation for plasticity; and (3) with a fitness cost for phenotypically plastic reactions. In almost all cases a generalist strategy to exploit both habitats emerged. Without any limits to plasticity, optimal adaptive reactions evolved. Without any genetic variation for plasticity, a compromise strategy with an intermediate, fixed phenotype evolved, whereas in the presence of costs a plastic compromise between the demands of the habitats and the costs associated with plasticity was found. Specialization and phenotypic differentiation was only found when selection within habitats was severe and optimal phenotypes for different habitats were widely different. Under soft selection (local regulation of population numbers in each habitat) the specialists coexisted; under hard selection (global regulation of population numbers) one specialist outcompeted the other. The prevalent evolutionary outcome of compromises rather than specialization implies that costs or constraints are not necessarily detectable as local adaptation in transplantation or translocation experiments.     

Status-dependent selection in the dimorphic beetle Onthophagus taurus.

The occurrence of alternative reproductive phenotypes is widespread in most animal taxa. The majority of known examples best fit the notion of alternative tactics within a conditional strategy where the fitness pay-offs depend on an individual's competitive ability or status. Individuals are proposed as "choosing" the tactic that maximizes their fitness, given their status relative to others in the population. Theoretically, status-dependent selection should determine when an animal should switch between alternative tactics. While a number of studies have demonstrated unequal fitness pay-offs associated with alternative tactics, none, to our knowledge, have examined the fitness functions necessary for predicting when individuals should switch between tactics. Here, we use a dimorphic male beetle in order to provide the first empirically derived fitness function across alternative reproductive phenotypes. Our data provide empirical support for a game-theoretic conditional strategy that has evolved under status-dependent selection.     

Evolutionary dynamics of habitat use

I examine the evolution of alternate genotypes that use two habitats that differ in vegetative cover, focusing on the interplay between ecological dynamics of the community and changes in selective advantage. Facultative habitat choice can stabilize a predator population that would cycle if isolated in the more open habitat. This has important implications for the evolution of habitat use strategies. Local stability arising from facultative habitat use allows any number of behavioural genotypes to co-exist: selective use of the open habitat, selective use of the dense habitat, opportunistic use of both habitats in proportion to availability, and facultative switching between habitats to maximize energy gain. Co-existence occurs because the fitness landscape is flat at the ecological equilibrium imposed by the facultative genotype. In contrast, ecological instability favours the evolution of genotypes with behavioural flexibility to avoid being in the wrong place at the wrong time or selective exploitation of one of the habitats. Uncertain information about habitat quality erodes the adaptive advantage of otherwise optimal behaviours, favouring a bet-hedging behavioural strategy synonymous with partial habitat preferences. These results suggest that ecological dynamics could have a strong influence on behavioural heterogeneity within forager populations and that a mixed ESS for habitat use should predominate.     

Consequences of divergent temperature optima in a host–parasite system

It was suggested that parasite infections become more severe with rising temperature, as expected during global warming. In ectothermic systems, the growth of a parasite and therefore its reproductive capacity is expected to increase with temperature. However, the outcome of the interaction depends on the temperature optima of both host and parasite. Here we used experimental infections of three‐spined stickleback fish Gasterosteus aculeatus with its specific tapeworm parasite Schistocephalus solidus to investigate in detail the temperature optima for both host and parasite. We analyzed the fitness consequences thereof, focusing on growth and immunity of the host, and growth and offspring production of the parasite as fitness correlates. We checked for potential differences among populations, using the offspring of hosts and parasites derived from four study sites in Iceland, Germany and Spain that differ in average annual temperature ranging between 4.8°C and 18.4°C. We found differences in temperature optima of host and parasites that were quite consistent across the populations: while sticklebacks grew faster and had higher immune activity at low temperatures, the parasites did not even grow fast enough to reach sexual maturity in these conditions. By contrast, with increasing temperatures, parasite growth, egg production and offspring hatching increased strongly while host immunity and growth were impaired. Our results show that divergent temperature optima of hosts and parasites can have drastic fitness consequences and support the expectation that some parasites will benefit from global warming.     

Influence of the rate of introduction on the fitness of restored populations

Most studies using demographic PVA models in a context of species restoration have concluded that rather than the rate of introduction, the total number of individuals released had the most important significant influence on the chance of success. In this article we use a genetic simulation model including deleterious and adaptive alleles to assess the impact of the method of release on the change in population mean fitness. We systematically compare a strategy that consists in releasing all individuals at the same time with a strategy that consists in staggering releases over a long period of time. Our results show that the former strategy is more beneficial for long-term fitness when considering advantageous genes only, while the latter is better when considering deleterious genes only. If deleterious and adaptive alleles are considered together, the best strategy depends then essentially on which of these types of alleles has the stronger influence on the change in total fitness. Although the relative contributions of the variance in total fitness due to adaptive and deleterious alleles may vary with the initial frequencies and the selective and dominance effects of these alleles, our results show that the optimal rate of release is mostly dependant on the expected long-term population size. Thus from a genetic view-point, the temporal release strategy of reintroduced populations should be considered with respect to their environment's carrying capacity.     

Plasticity of plant defense and its evolutionary implications in wild populations of Boechera stricta

Phenotypic plasticity is thought to impact evolutionary trajectories by shifting trait values in a direction that is either favored by natural selection (“adaptive” plasticity) or disfavored (“nonadaptive” plasticity). However, it is unclear how commonly each of these types of plasticity occurs in natural populations. To answer this question, we measured glucosinolate defensive chemistry and reproductive fitness in over 1500 individuals of the wild perennial mustard Boechera stricta, planted in four common gardens across central Idaho, United States. Glucosinolate profiles—including total glucosinolate concentration as well as the relative abundances and overall diversity of different compounds—were strongly plastic both among habitats and within habitats. Patterns of glucosinolate plasticity varied greatly among genotypes. Plasticity among sites was predicted to affect fitness in 27.1% of cases; more often than expected by chance, glucosinolate plasticity increased rather than decreased relative fitness. In contrast, we found no evidence for within‐habitat selection on glucosinolate reaction norm slopes (i.e., plasticity along a continuous environmental gradient). Together, our results indicate that glucosinolate plasticity may improve the ability of B. stricta populations to persist after migration to new habitats.     

Evolution in heterogeneous environments: Effects of migration on habitat specialization

Summary Richard Levins introduced fitness sets as a tool for investigating evolution within heterogeneous environments. Evolutionary game theory permits a synthesis and generalization of this approach by considering the evolutionary response of organisms to any scale of habitat heterogeneity. As scales of heterogeneity increase from fine to coarse, the evolutionary stable strategy (ESS) switches from a single generalist species to several species that become increasingly specialized on distinct habitats. Depending upon the organisms' ecology, the switch from one to two species may occur at high migration rates (relatively fine-grained environment), or may only occur at very low migration rates (coarse-grained environment). At the ESS, the evolutionary context of a species is the entire landscape, while its ecological context may be a single habitat.Evolution towards the ESS can be represented with adaptive landscapes. In the absence of frequency-dependence, shifting from a single strategy ESS to a two strategy ESS poses the problem of evolving across valleys in the adaptive surface to occupy new peaks (hence, Sewell Wright's shifting balance theory). Frequency-dependent processes facilitate evolution across valleys. If a system with a two strategy ESS is constrained to possess a single strategy, the population may actually evolve a strategy that minimizes fitness. Because the population now rests at the bottom of a valley, evolution by natural selection can drive populations to occupy both peaks.     

Demographic variation in cycad populations inhabiting contrasting forest fragments

Reduced habitat quality after fragmentation can significantly affect population viability, but the effects of differing quality of the remaining habitat on population fitness are rarely evaluated. Here, I compared fragmented populations of the cycad Zamia melanorrhachis from habitats with different history and subject to contrasting levels of disturbance to explore potential demographic differences in populations across habitat patches that could differ in habitat quality. Secondary-forest fragments had a lower canopy cover and soil moisture than remnant-forest fragments, which may represent a harsh environment for this cycad. A smaller average plant size and lower population density in the secondary-forest fragments support the hypothesis that these fragments may be of lower quality, e.g., if plants have reduced survival and/or fecundity in these habitats. However, variation in the stage-structure of populations (i.e., the relative proportions of non-reproductive and reproductive plants) was associated with the area of the forest fragments rather than the type of habitat (remnant versus secondary forest). These results suggest that different demographic parameters may respond differently to habitat fragmentation, which may be explained if processes like adult survival and recruitment depend on different characteristics of the habitat, e.g., average light/water availability versus suitable area for plant establishment. This study shows that forest fragments may differ drastically in environmental conditions and can sustain populations that can vary in their demography. Understanding how forest fragments may represent different habitat types is relevant for evaluating population viability in a heterogeneous landscape and for designing conservation programs that account for this heterogeneity.     

Interactions between soil habitat and geographic range location affect plant fitness

Populations are often found on different habitats at different geographic locations. This habitat shift may be due to biased dispersal, physiological tolerances or biotic interactions. To explore how fitness of the native plant Chamaecrista fasciculata depends on habitat within, at and beyond its range edge, we planted seeds from five populations in two soil substrates at these geographic locations. We found that with reduced competition, lifetime fitness was always greater or equivalent in one habitat type, loam soils, though early-season survival was greater on sand soils. At the range edge, natural populations are typically found on sand soil habitats, which are also less competitive environments. Early-season survival and fitness differed among source populations, and when transplanted beyond the range edge, range edge populations had greater fitness than interior populations. Our results indicate that even when the optimal soil substrate for a species does not change with geographic range location, the realized niche of a species may be restricted to sub-optimal habitats at the range edge because of the combined effects of differences in abiotic and biotic effects (e.g. competitors) between substrates.     

Stochasticity, complex spatial structure, and the feasibility of the shifting balance theory

Sewall Wright's shifting balance theory of evolution posits a mechanism by which a structured population may escape local fitness optima and find a global optimum. We examine a one-locus, two-allele model of underdominance in populations with differing spatial arrangements of demes, both analytically and with Monte Carlo simulations. We find that inclusion of variance in interpatch connectivities can significantly reduce the number of generations required for fixation of the more favorable allele relative to island and stepping-stone models. Although time to fixation increases with migration rate in all cases, the presence of one or two relatively isolated demes may reduce the number of generations by 80% or more. These results suggest that the shifting balance process may operate under less restrictive conditions than those found with a simple spatial arrangement of demes.     

Genetic algorithms and evolution

The genetic algorithm (GA) as developed by Holland (1975, Adaptation in Natural and Artificial Systems. Ann Arbor: University of Michigan Press) is an optimization technique based on natural selection. We use a modified version of this technique to investigate which aspects of natural selection make it an efficient search procedure. Our main modification to Holland's GA is the subdividing of the population into semi-isolated demes. We consider two examples. One is a fitness landscape with many local optima. The other is a model of singing in birds that has been previously analysed using dynamic programming. Both examples have epistatic interactions. In the first example we show that the GA can find the global optimum and that its success is improved by subdividing the population. In the second example we show that GAs can evolve to the optimal policy found by dynamic programming.     

The Fitness Consequences of Aneuploidy Are Driven by Condition-Dependent Gene Effects

Aneuploidy is a hallmark of tumor cells, and yet the precise relationship between aneuploidy and a cell’s proliferative ability, or cellular fitness, has remained elusive. In this study, we have combined a detailed analysis of aneuploid clones isolated from laboratory-evolved populations of Saccharomyces cerevisiae with a systematic, genome-wide screen for the fitness effects of telomeric amplifications to address the relationship between aneuploidy and cellular fitness. We found that aneuploid clones rise to high population frequencies in nutrient-limited evolution experiments and show increased fitness relative to wild type. Direct competition experiments confirmed that three out of four aneuploid events isolated from evolved populations were themselves sufficient to improve fitness. To expand the scope beyond this small number of exemplars, we created a genome-wide collection of >1,800 diploid yeast strains, each containing a different telomeric amplicon (Tamp), ranging in size from 0.4 to 1,000 kb. Using pooled competition experiments in nutrient-limited chemostats followed by high-throughput sequencing of strain-identifying barcodes, we determined the fitness effects of these >1,800 Tamps under three different conditions. Our data revealed that the fitness landscape explored by telomeric amplifications is much broader than that explored by single-gene amplifications. As also observed in the evolved clones, we found the fitness effects of most Tamps to be condition specific, with a minority showing common effects in all three conditions. By integrating our data with previous work that examined the fitness effects of single-gene amplifications genome-wide, we found that a small number of genes within each Tamp are centrally responsible for each Tamp’s fitness effects. Our genome-wide Tamp screen confirmed that telomeric amplifications identified in laboratory-evolved populations generally increased fitness. Our results show that Tamps are mutations that produce large, typically condition-dependent changes in fitness that are important drivers of increased fitness in asexually evolving populations.     

Competition between phenotypes

Summary We present two models for phenotypic-dependent interspecific competition. In both cases the survivorship of individuals of one population depends on the entire phenotypic distribution of the other species. The first model considers a continuously varying metric trait, with assortative or random mating; the second model examines a character controlled by two alleles at a single locus. Pursuing the notion that each population maximizes its mean fitness we define a vector-optimum strategy using the concepts of cooperative and competitive optima. It is found that the dynamical constraints placed on the equations of motion by Mendelian genetics often prevent a population from evolving to a strategic optimum. However, for the single locus case with complete dominance, the competitive optimum always coincides with some dynamical equilibrium on the Hardy-Weinberg manifold.     

Individual dispersal,landscape connectivity and ecological networks

Connectivity is classically considered an emergent property of landscapes encapsulating individuals' flows across space. However, its operational use requires a precise understanding of why and how organisms disperse. Such movements, and hence landscape connectivity, will obviously vary according to both organism properties and landscape features. We review whether landscape connectivity estimates could gain in both precision and generality by incorporating three fundamental outcomes of dispersal theory. Firstly, dispersal is a multi‐causal process; its restriction to an ‘escape reaction’ to environmental unsuitability is an oversimplification, as dispersing individuals can leave excellent quality habitat patches or stay in poor‐quality habitats according to the relative costs and benefits of dispersal and philopatry. Secondly, species, populations and individuals do not always react similarly to those cues that trigger dispersal, which sometimes results in contrasting dispersal strategies. Finally, dispersal is a major component of fitness and is thus under strong selective pressures, which could generate rapid adaptations of dispersal strategies. Such evolutionary responses will entail spatiotemporal variation in landscape connectivity. We thus strongly recommend the use of genetic tools to: (i) assess gene flow intensity and direction among populations in a given landscape; and (ii) accurately estimate landscape features impacting gene flow, and hence landscape connectivity. Such approaches will provide the basic data for planning corridors or stepping stones aiming at (re)connecting local populations of a given species in a given landscape. This strategy is clearly species‐ and landscape‐specific. But we suggest that the ecological network in a given landscape could be designed by stacking up such linkages designed for several species living in different ecosystems. This procedure relies on the use of umbrella species that are representative of other species living in the same ecosystem.     

THE SELECTION BARRIER BETWEEN POPULATIONS SUBJECT TO STABILIZING SELECTION

The introgression of genes carried by a small group of immigrants is studied. The recipient and the donor populations differ at several autosomal loci subject to weak selection, and two allelic forms of each gene are considered. Fitness variation is determined by additive allelic effects, by dominance effects, and by two-locus additive-by-additive epistatic interaction of the effects of the alleles. The fate of the group of immigrants is quantified by the selection barrier that describes the cumulative mean fitness of the hybrids and hybrid descendants relative to the fitness of the resident population. The monomorphic and the polymorphic loci of the recipient population contribute differently to the selection barrier. If the genetic difference between recipient and donor population is small, then the contribution of the monomorphic loci is dominated by a positive term dependent on the difference in gene frequencies. The contribution of the polymorphic loci depends only on the difference of the leading order in the pairwise linkage disequilibria between the two populations. This contribution may be positive or negative; and, thus, polymorphic loci may either contribute to the barrier or inflate the introgression.