Rapid evolution in response to introduced predators II: the contribution of adaptive plasticity

Background  

Introductions of non-native species can significantly alter the selective environment for populations of native species, which can respond through phenotypic plasticity or genetic adaptation. We examined phenotypic and genetic responses of Daphnia populations to recent introductions of non-native fish to assess the relative roles of phenotypic plasticity versus genetic change in causing the observed patterns. The Daphnia community in alpine lakes throughout the Sierra Nevada of California (USA) is ideally suited for investigation of rapid adaptive evolution because there are multiple lakes with and without introduced fish predators. We conducted common-garden experiments involving presence or absence of chemical cues produced by fish and measured morphological and life-history traits in Daphnia melanica populations collected from lakes with contrasting fish stocking histories. The experiment allowed us to assess the degree of population differentiation due to fish predation and examine the contribution of adaptive plasticity in the response to predator introduction.     

Adaptive divergence in body size overrides the effects of plasticity across natural habitats in the brown trout

The evolution of life-history traits is characterized by trade-offs between different selection pressures, as well as plasticity across environmental conditions. Yet, studies on local adaptation are often performed under artificial conditions, leaving two issues unexplored: (i) how consistent are laboratory inferred local adaptations under natural conditions and (ii) how much phenotypic variation is attributed to phenotypic plasticity and to adaptive evolution, respectively, across environmental conditions? We reared fish from six locally adapted (domesticated and wild) populations of anadromous brown trout (Salmo trutta) in one semi-natural and three natural streams and recorded a key life-history trait (body size at the end of first growth season). We found that population-specific reaction norms were close to parallel across different streams and Q_ST was similar – and larger than F_ST – within all streams, indicating a consistency of local adaptation in body size across natural environments. The amount of variation explained by population origin exceeded the variation across stream environments, indicating that genetic effects derived from adaptive processes have a stronger effect on phenotypic variation than plasticity induced by environmental conditions. These results suggest that plasticity does not “swamp” the phenotypic variation, and that selection may thus be efficient in generating genetic change.     

Rapid evolution of cold tolerance in stickleback

Climate change is predicted to lead to increased average temperatures and greater intensity and frequency of high and low temperature extremes, but the evolutionary consequences for biological communities are not well understood. Studies of adaptive evolution of temperature tolerance have typically involved correlative analyses of natural populations or artificial selection experiments in the laboratory. Field experiments are required to provide estimates of the timing and strength of natural selection, enhance understanding of the genetics of adaptation and yield insights into the mechanisms driving evolutionary change. Here, we report the experimental evolution of cold tolerance in natural populations of threespine stickleback fish (Gasterosteus aculeatus). We show that freshwater sticklebacks are able to tolerate lower minimum temperatures than marine sticklebacks and that this difference is heritable. We transplanted marine sticklebacks to freshwater ponds and measured the rate of evolution after three generations in this environment. Cold tolerance evolved at a rate of 0.63 haldanes to a value 2.5°C lower than that of the ancestral population, matching values found in wild freshwater populations. Our results suggest that cold tolerance is under strong selection and that marine sticklebacks carry sufficient genetic variation to adapt to changes in temperature over remarkably short time scales.     

Genomic signatures of experimental adaptive radiation in Drosophila

Abiotic environmental factors play a fundamental role in determining the distribution, abundance and adaptive diversification of species. Empowered by new technologies enabling rapid and increasingly accurate examination of genomic variation in populations, researchers may gain new insights into the genomic background of adaptive radiation and stress resistance. We investigated genomic variation across generations of large‐scale experimental selection regimes originating from a single founder population of Drosophila melanogaster, diverging in response to ecologically relevant environmental stressors: heat shock, heat knock down, cold shock, desiccation and starvation. When compared to the founder population, and to parallel unselected controls, there were more than 100,000 single nucleotide polymorphisms (SNPs) displaying consistent allelic changes in response to selective pressures across generations. These SNPs were found in both coding and noncoding sequences, with the highest density in promoter regions, and involved a broad range of functionalities, including molecular chaperoning by heat‐shock proteins. The SNP patterns were highly stressor‐specific despite considerable variation among line replicates within each selection regime, as reflected by a principal component analysis, and co‐occurred with selective sweep regions. Only ~15% of SNPs with putatively adaptive changes were shared by at least two selective regimes, while less than 1% of SNPs diverged in opposite directions. Divergent stressors driving evolution in the experimental system of adaptive radiation left distinct genomic signatures, most pronounced in starvation and heat‐shock selection regimes.     

Stabilizing selection and adaptive evolution in a combination of two traits in an arctic ungulate

Stabilizing selection is thought to be common in wild populations and act as one of the main evolutionary mechanisms, which constrain phenotypic variation. When multiple traits interact to create a combined phenotype, correlational selection may be an important process driving adaptive evolution. Here, we report on phenotypic selection and evolutionary changes in two natal traits in a semidomestic population of reindeer (Rangifer tarandus) in northern Finland. The population has been closely monitored since 1969, and detailed data have been collected on individuals since they were born. Over the length of the study period (1969–2015), we found directional and stabilizing selection toward a combination of earlier birth date and heavier birth mass with an intermediate optimum along the major axis of the selection surface. In addition, we demonstrate significant changes in mean traits toward earlier birth date and heavier birth mass, with corresponding genetic changes in breeding values during the study period. Our results demonstrate evolutionary changes in a combination of two traits, which agree closely with estimated patterns of phenotypic selection. Knowledge of the selective surface for combinations of genetically correlated traits are vital to predict how population mean phenotypes and fitness are affected when environments change.     

Natural selection drives patterns of lake-stream divergence in stickleback foraging morphology

To what extent are patterns of biological diversification determined by natural selection? We addressed this question by exploring divergence in foraging morphology of threespine stickleback fish inhabiting lake and stream habitats within eight independent watersheds. We found that lake fish generally displayed more developed gill structures and had more streamlined bodies than did stream fish. Diet analysis revealed that these morphological differences were associated with limnetic vs. benthic foraging modes, and that the extent of morphological divergence within watersheds reflected differences in prey resources utilized by lake and stream fish. We also found that patterns of divergence were unrelated to patterns of phenotypic trait (co)variance within populations (i.e. the ‘line of least resistance’). Instead, phenotypic (co)variances were more likely to have been shaped by adaptation to lake vs. stream habitats. Our study thus implicates natural selection as a strong deterministic force driving morphological diversification in lake–stream stickleback. The strength of this inference was obtained by complementing a standard analysis of parallel divergence in means between discrete habitat categories (lake vs. stream) with quantitative estimates of selective forces and information on trait (co)variances.     

The adaptive dynamics of niche constructing traits in spatially subdivided populations: evolving posthumous extended phenotypes

Niche construction, by which organisms modify the environment in which they live, has been proposed to affect the evolution of many phenotypic traits. But what about the evolution of a niche constructing trait itself, whose expression changes the pattern of natural selection to which the trait is exposed in subsequent generations? This article provides an inclusive fitness analysis of selection on niche constructing phenotypes, which can affect their environment from local to global scales in arbitrarily spatially subdivided populations. The model shows that phenotypic effects of genes extending far beyond the life span of the actor can be affected by natural selection, provided they modify the fitness of those individuals living in the future that are likely to have inherited the niche construction lineage of the actor. Present benefits of behaviors are thus traded off against future indirect costs. The future costs will generally result from a complicated interplay of phenotypic effects, population demography and environmental dynamics. To illustrate these points, I derive the adaptive dynamics of a trait involved in the consumption of an abiotic resource, where resource abundance in future generations feeds back to the evolutionary dynamics of the trait.     

Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae)

Natural selection drives local adaptation, potentially even at small temporal and spatial scales. As a result, adaptive genetic and phenotypic divergence can occur among populations living in different habitats. We investigated patterns of differentiation between contrasting lake and stream habitats in the cyprinid fish European minnow (Phoxinus phoxinus) at both the morphological and genomic levels using geometric morphometrics and AFLP markers, respectively. We also used a spatial correlative approach to identify AFLP loci associated with environmental variables representing potential selective forces responsible for adaptation to divergent habitats. Our results identified different morphologies between lakes and streams, with lake fish presenting a deeper body and caudal peduncle compared to stream fish. Body shape variation conformed to a priori predictions concerning biomechanics and swimming performance in lakes vs. streams. Moreover, morphological differentiation was found to be associated with several environmental variables, which could impose selection on body and caudal peduncle shape. We found adaptive genetic divergence between these contrasting habitats in the form of 'outlier' loci (2.9%) whose genetic divergence exceeded neutral expectations. We also detected additional loci (6.6%) not associated with habitat type (lake vs. stream), but contributing to genetic divergence between populations. Specific environmental variables related to trophic dynamics, landscape topography and geography were associated with several neutral and outlier loci. These results provide new insights into the morphological divergence and genetic basis of adaptation to differentiated habitats.     

Reversal of evolutionary downsizing caused by selective harvest of large fish

Evolutionary responses to the long-term exploitation of individuals from a population may include reduced growth rate, age at maturation, body size and productivity. Theoretical models suggest that these genetic changes may be slow or impossible to reverse but rigorous empirical evidence is lacking. Here, we provide the first empirical demonstration of a genetically based reversal of fishing-induced evolution. We subjected six populations of silverside fish (Menidia menidia) to three forms of size-selective fishing for five generations, thereby generating twofold differences among populations in mean weight and yield (biomass) at harvest. This was followed by an additional five generations during which size-selective harvest was halted. We found that evolutionary changes were reversible. Populations evolving smaller body size when subjected to size-selective fishing displayed a slow but significant increase in size when fishing ceased. Neither phenotypic variance in size nor juvenile survival was reduced by the initial period of selective fishing, suggesting that sufficient genetic variation remained to allow recovery. By linear extrapolation, we predict full recovery in about 12 generations, although the rate of recovery may taper off near convergence. The recovery rate in any given wild population will also depend on other agents of selection determined by the specifics of life history and environment. By contrast, populations that in the first five generations evolved larger size and yield showed little evidence of reversal. These results show that populations have an intrinsic capacity to recover genetically from harmful evolutionary changes caused by fishing, even without extrinsic factors that reverse the selection gradient. However, harvested species typically have generation times of 3–7 years, so recovery may take decades. Hence, the need to account for evolution in managing fisheries remains.     

Drainage-structuring of ancestral variation and a common functional pathway shape limited genomic convergence in natural high- and low-predation guppies

Studies of convergence in wild populations have been instrumental in understanding adaptation by providing strong evidence for natural selection. At the genetic level, we are beginning to appreciate that the re-use of the same genes in adaptation occurs through different mechanisms and can be constrained by underlying trait architectures and demographic characteristics of natural populations. Here, we explore these processes in naturally adapted high- (HP) and low-predation (LP) populations of the Trinidadian guppy, Poecilia reticulata. As a model for phenotypic change this system provided some of the earliest evidence of rapid and repeatable evolution in vertebrates; the genetic basis of which has yet to be studied at the whole-genome level. We collected whole-genome sequencing data from ten populations (176 individuals) representing five independent HP-LP river pairs across the three main drainages in Northern Trinidad. We evaluate population structure, uncovering several LP bottlenecks and variable between-river introgression that can lead to constraints on the sharing of adaptive variation between populations. Consequently, we found limited selection on common genes or loci across all drainages. Using a pathway type analysis, however, we find evidence of repeated selection on different genes involved in cadherin signaling. Finally, we found a large repeatedly selected haplotype on chromosome 20 in three rivers from the same drainage. Taken together, despite limited sharing of adaptive variation among rivers, we found evidence of convergent evolution associated with HP-LP environments in pathways across divergent drainages and at a previously unreported candidate haplotype within a drainage.     

A test for the genetic basis of natural selection: an individual-based longitudinal study in a stream-dwelling fish

In addition to the well-studied evolutionary parameters of (1) phenotype-fitness covariance and (2) the genetic basis of phenotypic variation, adaptive evolution by natural selection requires that (3) fitness variation is effected by heritable genetic differences among individuals and (4) phenotype-fitness covariances must be, at least in part, underlain by genetic covariances. These latter two requirements for adaptive evolutionary change are relatively unstudied in natural populations. Absence of the latter requirements could explain stasis of apparently directionally selected heritable traits. We provide complementary analyses of selection and variation at phenotypic and genetic levels for juvenile growth rate in brook charr Salvelinus fontinalis in Freshwater River, Newfoundland, Canada. Contrary to the vast majority of reports in fish, we found very little viability selection of juvenile body size. Large body size appears nonetheless to be selectively advantageous via a relationship with early maturity. Genetic patterns in evolutionary parameters largely reflected phenotypic patterns. We have provided inference of selection based on longitudinal data, which are uncommon in high fecundity organisms. Furthermore we have provided a practicable framework for further studies of the genetic basis of natural selection.     

Relaxed predation results in reduced phenotypic integration in a suite of dragonflies

Although changes in magnitude of single traits responding to selective agents have been studied intensively, little is known about selection shaping networks of traits and their patterns of covariation. However, this is central for our understanding of phenotypic evolution as traits are embedded in a multivariate environment with selection affecting a multitude of traits simultaneously rather than individually. Here, we investigate inter‐ and intraspecific patterns of trait integration (trait correlations) in the larval abdomen of dragonflies as a response to a change in predator selection. Species of the dragonfly genus Leucorrhinia underwent a larval habitat shift from predatory fish to predatory dragonfly‐dominated lakes with an associated relaxation in selection pressure from fish predation. Our results indicate that the habitat‐shift‐induced relaxed selection pressure caused phenotypic integration of abdominal traits to be reduced. Intraspecific findings matched patterns comparing species from both habitats with higher abdominal integration in response to predatory fish. This higher integration is probably a result of faster burst swimming speed. The abdomen holds the necessary morphological machinery to successfully evade predatory fish via burst swimming. Hence, abdominal traits have to function in a tight coordinated manner, as maladaptive variation and consequently nonoptimal burst swimming would cause increased mortality. In predatory dragonfly‐dominated lakes, no such strong link between burst swimming and mortality is present. Our findings highlight the importance of studying multivariate trait relationships as a response to selection for understanding patterns of phenotypic diversification.     

Contrasting post-settlement selection results in many-to-one mapping of high performance phenotypes in the Hawaiian waterfall-climbing goby <Emphasis Type="Italic">Sicyopterus stimpsoni</Emphasis>

Natural selection drives adaptive evolution, but contrasting environmental pressures may lead to trade-offs between phenotypes that confer different performances. Such trade-offs may weaken the strength of selection and/or generate complex fitness surfaces with multiple local optima that correspond to different selection regimes. We evaluated how differences in patterns of phenotypic selection might promote morphological differences between subpopulations of the amphidromous Hawaiian waterfall-climbing goby, Sicyopterus stimpsoni. We conducted laboratory experiments on fish from the islands of Kaua‘i and Hawai‘i (the “Big Island”) to compare patterns of linear and nonlinear selection, and the opportunity for selection, that result from two contrasting pressures, predator evasion and waterfall climbing, which vary in intensity between islands. We found directional and nonlinear selection were strongest when individuals were exposed to their primary selective pressures (predator evasion on Kaua‘i, waterfall climbing on the Big Island). However, the opportunity for selection was greater for the non-primary pressure: climbing on Kaua‘i, predator evasion on the Big Island. Canonical rotation of the nonlinear gamma matrix demonstrated that individuals from Kaua‘i and the Big Island occupy regions near their local fitness peaks for some traits. Therefore, selection for predator evasion on Kaua‘i and climbing on the Big Island may be less effective in promoting morphological changes in this species, because variation of functionally important traits in their respective environments may have been reduced by directional or stabilizing selection. These results demonstrate that despite constraints on the opportunities for selection, population differences in phenotypic traits can arise due to differences in selective regimes. For S. stimpsoni, sufficient variation exists in other locomotor traits, allowing for necessary levels of performance in the contrasting selective regime (i.e., climbing on Kaua‘i and predator evasion on the Big Island) through many-to-one-mapping, which may be essential for the survival of local populations in an evanescent island environment.     

Ecological causes of morphological evolution in the three‐spined stickleback

The central assumption of evolutionary theory is that natural selection drives the adaptation of populations to local environmental conditions, resulting in the evolution of adaptive phenotypes. The three‐spined stickleback (Gasterosteus aculeatus) displays remarkable phenotypic variation, offering an unusually tractable model for understanding the ecological mechanisms underpinning adaptive evolutionary change. Using populations on North Uist, Scotland we investigated the role of predation pressure and calcium limitation on the adaptive evolution of stickleback morphology and behavior. Dissolved calcium was a significant predictor of plate and spine morph, while predator abundance was not. Stickleback latency to emerge from a refuge varied with morph, with populations with highly reduced plates and spines and high predation risk less bold. Our findings support strong directional selection in three‐spined stickleback evolution, driven by multiple selective agents.     

The repeatability of adaptive radiation during long-term experimental evolution of Escherichia coli in a multiple nutrient environment

Adaptive radiations occur when a species diversifies into different ecological specialists due to competition for resources and trade-offs associated with the specialization. The evolutionary outcome of an instance of adaptive radiation cannot generally be predicted because chance (stochastic events) and necessity (deterministic events) contribute to the evolution of diversity. With increasing contributions of chance, the degree of parallelism among different instances of adaptive radiations and the predictability of an outcome will decrease. To assess the relative contributions of chance and necessity during adaptive radiation, we performed a selection experiment by evolving twelve independent microcosms of Escherichia coli for 1000 generations in an environment that contained two distinct resources. Specialization to either of these resources involves strong trade-offs in the ability to use the other resource. After selection, we measured three phenotypic traits: 1) fitness, 2) mean colony size, and 3) colony size diversity. We used fitness relative to the ancestor as a measure of adaptation to the selective environment; changes in colony size as a measure of the evolution of new resource specialists because colony size has been shown to correlate with resource specialization; and colony size diversity as a measure of the evolved ecological diversity. Resource competition led to the rapid evolution of phenotypic diversity within microcosms. Measurements of fitness, colony size, and colony size diversity within and among microcosms showed that the repeatability of adaptive radiation was high, despite the evolution of genetic variation within microcosms. Consistent with the observation of parallel evolution, we show that the relative contributions of chance are far smaller and less important than effects due to adaptation for the traits investigated. The two-resource environment imposed similar selection pressures in independent populations and promoted parallel phenotypic adaptive radiations in all independently evolved microcosms.     

Rapid,broad‐scale gene expression evolution in experimentally harvested fish populations

Gene expression changes potentially play an important role in adaptive evolution under human‐induced selection pressures, but this has been challenging to demonstrate in natural populations. Fishing exhibits strong selection pressure against large body size, thus potentially inducing evolutionary changes in life history and other traits that may be slowly reversible once fishing ceases. However, there is a lack of convincing examples regarding the speed and magnitude of fisheries‐induced evolution, and thus, the relevant underlying molecular‐level effects remain elusive. We use wild‐origin zebrafish (Danio rerio) as a model for harvest‐induced evolution. We experimentally demonstrate broad‐scale gene expression changes induced by just five generations of size‐selective harvesting, and limited genetic convergence following the cessation of harvesting. We also demonstrate significant allele frequency changes in genes that were differentially expressed after five generations of size‐selective harvesting. We further show that nine generations of captive breeding induced substantial gene expression changes in control stocks likely due to inadvertent selection in the captive environment. The large extent and rapid pace of the gene expression changes caused by both harvest‐induced selection and captive breeding emphasizes the need for evolutionary enlightened management towards sustainable fisheries.     

Evolution of the additive genetic variance–covariance matrix under continuous directional selection on a complex behavioural phenotype

Given the pace at which human-induced environmental changes occur, a pressing challenge is to determine the speed with which selection can drive evolutionary change. A key determinant of adaptive response to multivariate phenotypic selection is the additive genetic variance–covariance matrix (G). Yet knowledge of G in a population experiencing new or altered selection is not sufficient to predict selection response because G itself evolves in ways that are poorly understood. We experimentally evaluated changes in G when closely related behavioural traits experience continuous directional selection. We applied the genetic covariance tensor approach to a large dataset (n = 17 328 individuals) from a replicated, 31-generation artificial selection experiment that bred mice for voluntary wheel running on days 5 and 6 of a 6-day test. Selection on this subset of G induced proportional changes across the matrix for all 6 days of running behaviour within the first four generations. The changes in G induced by selection resulted in a fourfold slower-than-predicted rate of response to selection. Thus, selection exacerbated constraints within G and limited future adaptive response, a phenomenon that could have profound consequences for populations facing rapid environmental change.     

Monitoring adaptive genetic responses to environmental change

Widespread environmental changes including climate change, selective harvesting and landscape alterations now greatly affect selection regimes for most organisms. How animals and plants can adapt to these altered environments via contemporary evolution is thus of strong interest. We discuss how to use genetic monitoring to study adaptive responses via repeated analysis of the same populations over time, distinguishing between phenotypic and molecular genetics approaches. After describing monitoring designs, we develop explicit criteria for demonstrating adaptive responses, which include testing for selection and establishing clear links between genetic and environmental change. We then review a few exemplary studies that explore adaptive responses to climate change in Drosophila, selective responses to hunting and fishing, and contemporary evolution in Daphnia using resurrected resting eggs. We further review a broader set of 44 studies to assess how well they meet the proposed criteria, and conclude that only 23% fulfill all criteria. Approximately half (43%) of these studies failed to rule out the alternative hypothesis of replacement by a different, better-adapted population. Likewise, 34% of the studies based on phenotypic variation did not test for selection as opposed to drift. These shortcomings can be addressed via improved experimental designs and statistical testing. We foresee monitoring of adaptive responses as a future valuable tool in conservation biology, for identifying populations unable to evolve at sufficiently high rates and for identifying possible donor populations for genetic rescue. Technological advances will further augment the realization of this potential, especially next-generation sequencing technologies that allow for monitoring at the level of whole genomes.     

Secondary contact seeds phenotypic novelty in cichlid fishes

Theory proposes that genomic admixture between formerly reproductively isolated populations can generate phenotypic novelty for selection to act upon. Secondary contact may therefore be a significant promoter of phenotypic novelty that allows species to overcome environmental challenges and adapt to novel environments, including during adaptive radiation. To date, this has largely been considered from the perspective of interspecific hybridization at contact zones. However, it is also possible that this process occurs more commonly between natural populations of a single species, and thus its importance in adaptive evolution may have been underestimated. In this study, we tested the consequences of genomic introgression during apparent secondary contact between phenotypically similar lineages of the riverine cichlid fish Astatotilapia calliptera. We provide population genetic evidence of a secondary contact zone in the wild, and then demonstrate using mate-choice experiments that both lineages can reproduce together successfully in laboratory conditions. Finally, we show that genomically admixed individuals display extreme phenotypes not observed in the parental lineages. Collectively, the evidence shows that secondary contact can drive the evolution of phenotypic novelty, suggesting that pulses of secondary contact may repeatedly seed genetic novelty, which when coupled with ecological opportunity could promote rapid adaptive evolution in natural circumstances.