Gene expression plasticity evolves in response to colonization of freshwater lakes in threespine stickleback

Phenotypic plasticity is predicted to facilitate individual survival and/or evolve in response to novel environments. Plasticity that facilitates survival should both permit colonization and act as a buffer against further evolution, with contemporary and derived forms predicted to be similarly plastic for a suite of traits. On the other hand, given the importance of plasticity in maintaining internal homeostasis, derived populations that encounter greater environmental heterogeneity should evolve greater plasticity. We tested the evolutionary significance of phenotypic plasticity in coastal British Columbian postglacial populations of threespine stickleback (Gasterosteus aculeatus) that evolved under greater seasonal extremes in temperature after invading freshwater lakes from the sea. Two ancestral (contemporary marine) and two derived (contemporary freshwater) populations of stickleback were raised near their thermal tolerance extremes, 7 and 22 °C. Gene expression plasticity was estimated for more than 14 000 genes. Over five thousand genes were similarly plastic in marine and freshwater stickleback, but freshwater populations exhibited significantly more genes with plastic expression than marine populations. Furthermore, several of the loci shown to exhibit gene expression plasticity have been previously implicated in the adaptive evolution of freshwater populations, including a gene involved in mitochondrial regulation (PPARAa). Collectively, these data provide molecular evidence that highlights the importance of plasticity in colonization and adaptation to new environments.     

DISENTANGLING THE ROLE OF PHENOTYPIC PLASTICITY AND GENETIC DIVERGENCE IN CONTEMPORARY ECOTYPE FORMATION DURING A BIOLOGICAL INVASION

The occurrence of contemporary ecotype formation through adaptive divergence of populations within the range of an invasive species typically requires standing genetic variation but can be facilitated by phenotypic plasticity. The relative contributions of both of these to adaptive trait differentiation have rarely been simultaneously quantified in recently diverging vertebrate populations. Here we study a case of intraspecific divergence into distinct lake and stream ecotypes of threespine stickleback that evolved in the past 140 years within the invasive range in Switzerland. Using a controlled laboratory experiment with full‐sib crosses and treatments mimicking a key feature of ecotypic niche divergence, we test if the phenotypic divergence that we observe in the wild results from phenotypic plasticity or divergent genetic predisposition. Our experimental groups show qualitatively similar phenotypic divergence as those observed among wild adults. The relative contribution of plasticity and divergent genetic predisposition differs among the traits studied, with traits related to the biomechanics of feeding showing a stronger genetic predisposition, whereas traits related to locomotion are mainly plastic. These results implicate that phenotypic plasticity and standing genetic variation interacted during contemporary ecotype formation in this case.     

A test of the "flexible stem" model of evolution: ancestral plasticity, genetic accommodation, and morphological divergence in the threespine stickleback radiation

If an ancestral stem group repeatedly colonizes similar environments, developmental plasticity specific to that group should consistently give rise to similar phenotypes. Parallel selection on those similar phenotypes could lead to the repeated evolution of characteristic ecotypes, a property common to many adaptive radiations. A key prediction of this "flexible stem" model of adaptive radiation is that patterns of phenotypic divergence in derived groups should mirror patterns of developmental plasticity in their common ancestor. The threespine stickleback radiation provides an excellent opportunity to test this prediction because the marine form is representative of the ancestral stem group, which has repeatedly given rise to several characteristic ecotypes. We examined plasticity of several aspects of shape and trophic morphology in response to diets characteristic of either the derived benthic ecotype or the limnetic ecotype. When marine fish were reared on alternative diets, plasticity of head and mouth shape paralleled phenotypic divergence between the derived ecotypes, supporting the flexible stem model. Benthic and limnetic fish exhibited patterns of plasticity similar to those of the marine population; however, some differences in population means were present, as well as subtle differences in shape plasticity in the benthic population, indicating a role for genetic accommodation in this system.     

Should I stay or should I go? The Ectodysplasin locus is associated with behavioural differences in threespine stickleback

Adaptive divergence may be facilitated if morphological and behavioural traits associated with local adaptation share the same genetic basis. It is therefore important to determine whether genes underlying adaptive morphological traits are associated with variation in behaviour in natural populations. Positive selection on low-armour alleles at the Ectodysplasin (Eda) locus in threespine stickleback has led to the repeated evolution of reduced armour, following freshwater colonization by fully armoured marine sticklebacks. This adaptive divergence in armour between marine and freshwater populations would be facilitated if the low allele conferred a behavioural preference for freshwater environments. We experimentally tested whether the low allele is associated with preference for freshwater by measuring the preference of each Eda genotype for freshwater versus saltwater after acclimation to either salinity. We found no association between the Eda low allele and preference for freshwater. Instead, the low allele was significantly associated with a reduced preference for the acclimation environment. This behaviour may facilitate the colonization of freshwater habitats from the sea, but could also hinder local adaptation by promoting migration of low alleles between marine and freshwater environments.     

Comparative Analysis of Japanese Three-Spined Stickleback Clades Reveals the Pacific Ocean Lineage Has Adapted to Freshwater Environments while the Japan Sea Has Not

Divergent selection and adaptive divergence can increase phenotypic diversification amongst populations and lineages. Yet adaptive divergence between different environments, habitats or niches does not occur in all lineages. For example, the colonization of freshwater environments by ancestral marine species has triggered adaptive radiation and phenotypic diversification in some taxa but not in others. Studying closely related lineages differing in their ability to diversify is an excellent means of understanding the factors promoting and constraining adaptive evolution. A well-known example of the evolution of increased phenotypic diversification following freshwater colonization is the three-spined stickleback. Two closely related stickleback lineages, the Pacific Ocean and the Japan Sea occur in Japan. However, Japanese freshwater stickleback populations are derived from the Pacific Ocean lineage only, suggesting the Japan Sea lineage is unable to colonize freshwater. Using stable isotope data and trophic morphology, we first show higher rates of phenotypic and ecological diversification between marine and freshwater populations within the Pacific Ocean lineage, confirming adaptive divergence has occurred between the two lineages and within the Pacific Ocean lineage but not in the Japan Sea lineage. We further identified consistent divergence in diet and foraging behaviour between marine forms from each lineage, confirming Pacific Ocean marine sticklebacks, from which all Japanese freshwater populations are derived, are better adapted to freshwater environments than Japan Sea sticklebacks. We suggest adaptive divergence between ancestral marine populations may have played a role in constraining phenotypic diversification and adaptive evolution in Japanese sticklebacks.     

Contemporary ecotypic divergence during a recent range expansion was facilitated by adaptive introgression

Although rapid phenotypic evolution during range expansion associated with colonization of contrasting habitats has been documented in several taxa, the evolutionary mechanisms that underlie such phenotypic divergence have less often been investigated. A strong candidate for rapid ecotype formation within an invaded range is the three‐spine stickleback in the Lake Geneva region of central Europe. Since its introduction only about 140 years ago, it has undergone a significant expansion of its range and its niche, now forming phenotypically differentiated parapatric ecotypes that occupy either the pelagic zone of the large lake or small inlet streams, respectively. By comparing museum collections from different times with contemporary population samples, we here reconstruct the evolution of parapatric phenotypic divergence through time. Using genetic data from modern samples, we infer the underlying invasion history. We find that parapatric habitat‐dependent phenotypic divergence between the lake and stream was already present in the first half of the twentieth century, but the magnitude of differentiation increased through time, particularly in antipredator defence traits. This suggests that divergent selection between the habitats occurred and was stable through much of the time since colonization. Recently, increased phenotypic differentiation in antipredator defence traits likely results from habitat‐dependent selection on alleles that arrived through introgression from a distantly related lineage from outside the Lake Geneva region. This illustrates how hybridization can quickly promote phenotypic divergence in a system where adaptation from standing genetic variation was constrained.     

ADAPTIVE DIVERGENCE BETWEEN FRESHWATER AND MARINE STICKLEBACKS: INSIGHTS INTO THE ROLE OF PHENOTYPIC PLASTICITY FROM AN INTEGRATED ANALYSIS OF CANDIDATE GENE EXPRESSION

Debate surrounding the integration of phenotypic plasticity within the neo‐Darwinian paradigm has recently intensified, but is largely dominated by conceptual abstractions. Advances in our capacities to identify candidate genes, and quantify their levels of expression, now facilitate the study of natural variation in inherently plastic traits, and may lead to a more concrete understanding of plasticity's role in adaptive evolution. We present data from parapatric threespine stickleback (Gasterosteus aculeatus) demes inhabiting geologically recent, freshwater and saltwater zones of a large estuary. Reaction norms for survival confirm adaptation to local salinity conditions. Analysis of osmoregulatory candidate gene expression within an ecological quantitative genetics framework suggests putative mechanisms underlying adaptive variation, and provides insights into the role of ancestral trait plasticity in this divergence. A sodium–potassium ATPase (ATP1A1) is identified as a candidate gene for freshwater adaptation. In addition to heritable variation for gene expression, we infer significant correlation between measures of expression and individual fitness. Overall results indicate a loss of plasticity in the freshwater deme. We discuss how this is consistent with adaptation facilitated by ancestral plasticity as a heuristic example that may prove useful for future, explicit tests of the genetic assimilation hypothesis.     

Phenotypic divergence among threespine stickleback that differ in nuptial coloration

By studying systems in their earliest stages of differentiation, we can learn about the evolutionary forces acting within and among populations and how those forces could contribute to reproductive isolation. Such an understanding would help us to better discern and predict how selection leads to the maintenance of multiple morphs within a species, rather than speciation. The postglacial adaptive radiation of the threespine stickleback (Gasterosteus aculeatus) is one of the best‐studied cases of evolutionary diversification and rapid, repeated speciation. Following deglaciation, marine stickleback have continually invaded freshwater habitats across the northern hemisphere and established resident populations that diverged innumerable times from their oceanic ancestors. Independent freshwater colonization events have yielded broadly parallel patterns of morphological differences in freshwater and marine stickleback. However, there is also much phenotypic diversity within and among freshwater populations. We studied a lesser‐known freshwater “species pair” found in southwest Washington, where male stickleback in numerous locations have lost the ancestral red sexual signal and instead develop black nuptial coloration. We measured phenotypic variation in a suite of traits across sites where red and black stickleback do not overlap in distribution and at one site where they historically co‐occurred. We found substantial phenotypic divergence between red and black morphs in noncolor traits including shape and lateral plating, and additionally find evidence that supports the hypothesis of sensory drive as the mechanism responsible for the evolutionary switch in color from red to black. A newly described third “mixed” morph in Connor Creek, Washington, differs in head shape and size from the red and black morphs, and we suggest that their characteristics are most consistent with hybridization between anadromous and freshwater stickleback. These results lay the foundation for future investigation of the underlying genetic basis of this phenotypic divergence as well as the evolutionary processes that may drive, maintain, or limit divergence among morphs.     

Phenotypic Plasticity and the Ecotypic Differentiation of Aggressive Behavior in Threespine Stickleback

The threespine stickleback fish, Gasterosteus aculeatus, has undergone a remarkable postglacial adaptive radiation in which an ancient oceanic ancestor has given rise to uncountable freshwater populations. The radiation is characterized by repeated, independent evolution of similar derived phenotypes under similar environmental conditions. A common pattern of divergence is caused by differences in habitat that favor morphological and behavioral features that enhance efficiency of feeding on plankton (limnetic ecotypes) vs. those that enhance efficiency of feeding on benthic invertebrates (benthic ecotypes). These two ecotypes exhibit consistently different patterns of courtship and of foraging and cannibalistic behavior (divergent behavioral syndromes). Here, we demonstrate that there also exist differences in aggression toward conspecifics that are likely to be characteristic of the ecotypes. We report differences in patterns of aggression toward rivals between the ecotypes and offer evidence of differences in the patterns of phenotypic plasticity (norms of reaction) for these traits across population types, and of differences in the incorporation of aggressive elements of behavior in courtship. These data support an earlier suggestion that differences in aggressive tendencies could have facilitated assortative mating between the four benthic–limnetic species pairs found in British Columbia lakes, and they demonstrate the need to evaluate divergent behavioral phenotypes in this radiation as phenotypic norms of reaction rather than as fixed traits.     

Do floral and ecogeographic isolation allow the co‐occurrence of two ecotypes of Anacamptis papilionacea (Orchidaceae)?

Ecotypes are relatively frequent in flowering plants and considered central in ecological speciation as local adaptation can promote the insurgence of reproductive isolation. Without geographic isolation, gene flow usually homogenizes the allopatrically generated phenotypic and ecological divergences, unless other forms of reproductive isolation keep them separated. Here, we investigated two orchid ecotypes with marked phenotypic floral divergence that coexist in contact zones. We found that the two ecotypes show different ecological habitat preferences with one being more climatically restricted than the other. The ecotypes remain clearly morphologically differentiated both in allopatry and in sympatry and differed in diverse floral traits. Despite only slightly different flowering times, the two ecotypes achieved floral isolation thanks to different pollination strategies. We found that both ecotypes attract a wide range of insects, but the ratio of male/female attracted by the two ecotypes was significantly different, with one ecotype mainly attracts male pollinators, while the other mainly attracts female pollinators. As a potential consequence, the two ecotypes show different pollen transfer efficiency. Experimental plots with pollen staining showed a higher proportion of intra‐ than interecotype movements confirming floral isolation between ecotypes in sympatry while crossing experiments excluded evident postmating barriers. Even if not completely halting the interecotypes pollen flow in sympatry, such incipient switch in pollination strategy between ecotypes may represent a first step on the path toward evolution of sexual mimicry in Orchidinae.     

LOCAL ADAPTATION AND THE EVOLUTION OF PHENOTYPIC PLASTICITY IN TRINIDADIAN GUPPIES (POECILIA RETICULATA)

Divergent selection pressures across environments can result in phenotypic differentiation that is due to local adaptation, phenotypic plasticity, or both. Trinidadian guppies exhibit local adaptation to the presence or absence of predators, but the degree to which predator‐induced plasticity contributes to population differentiation is less clear. We conducted common garden experiments on guppies obtained from two drainages containing populations adapted to high‐ and low‐predation environments. We reared full‐siblings from all populations in treatments simulating the presumed ancestral (predator cues present) and derived (predator cues absent) conditions and measured water column use, head morphology, and size at maturity. When reared in presence of predator cues, all populations had phenotypes that were typical of a high‐predation ecotype. However, when reared in the absence of predator cues, guppies from high‐ and low‐predation regimes differed in head morphology and size at maturity; the qualitative nature of these differences corresponded to those that characterize adaptive phenotypes in high‐ versus low‐predation environments. Thus, divergence in plasticity is due to phenotypic differences between high‐ and low‐predation populations when reared in the absence of predator cues. These results suggest that plasticity might initially play an important role during colonization of novel environments, and then evolve as a by‐product of adaptation to the derived environment.     

Quick divergence but slow convergence during ecotype formation in lake and stream stickleback pairs of variable age

When genetic constraints restrict phenotypic evolution, diversification can be predicted to evolve along so‐called lines of least resistance. To address the importance of such constraints and their resolution, studies of parallel phenotypic divergence that differ in their age are valuable. Here, we investigate the parapatric evolution of six lake and stream threespine stickleback systems from Iceland and Switzerland, ranging in age from a few decades to several millennia. Using phenotypic data, we test for parallelism in ecotypic divergence between parapatric lake and stream populations and compare the observed patterns to an ancestral‐like marine population. We find strong and consistent phenotypic divergence, both among lake and stream populations and between our freshwater populations and the marine population. Interestingly, ecotypic divergence in low‐dimensional phenotype space (i.e. single traits) is rapid and seems to be often completed within 100 years. Yet, the dimensionality of ecotypic divergence was highest in our oldest systems and only there parallel evolution of unrelated ecotypes was strong enough to overwrite phylogenetic contingency. Moreover, the dimensionality of divergence in different systems varies between trait complexes, suggesting different constraints and evolutionary pathways to their resolution among freshwater systems.     

Lack of parallel genetic patterns underlying the repeated ecological divergence of beach and stream‐spawning kokanee salmon

Recent progress in methods for detecting adaptive population divergence in situ shows promise for elucidating the conditions under which selection acts to generate intraspecific diversity. Rapid ecological diversification is common in fishes; however, the role of phenotypic plasticity and adaptation to local environments is poorly understood. It is now possible to investigate genetic patterns to make inferences regarding phenotypic traits under selection and possible mechanisms underlying ecotype divergence, particularly where similar novel phenotypes have arisen in multiple independent populations. Here, we employed a bottom‐up approach to test for signatures of directional selection associated with divergence of beach‐ and stream‐spawning kokanee, the obligate freshwater form of sockeye salmon (Oncorhynchus nerka). Beach‐ and stream‐spawners co‐exist in many post‐glacial lakes and exhibit distinct reproductive behaviours, life‐history traits and spawning habitat preferences. Replicate ecotype pairs across five lakes in British Columbia, Canada were genotyped at 57 expressed sequence tag‐linked and anonymous microsatellite loci identified in a previous genome scan. Fifteen loci exhibited signatures of directional selection (high F_ST outliers), four of which were identified in multiple lakes. However, the lack of parallel genetic patterns across all lakes may be a result of: 1) an inability to detect loci truly under selection; 2) alternative genetic pathways underlying ecotype divergence in this system; and/or 3) phenotypic plasticity playing a formative role in driving kokanee spawning habitat differences. Gene annotations for detected outliers suggest pathogen resistance and energy metabolism as potential mechanisms contributing to the divergence of beach‐ and stream‐spawning kokanee, but further study is required.     

Sexual isolation promotes divergence between parapatric lake and stream stickleback

Speciation can be initiated by adaptive divergence between populations in ecologically different habitats, but how sexually based reproductive barriers contribute to this process is less well understood. We here test for sexual isolation between ecotypes of threespine stickleback fish residing in adjacent lake and stream habitats in the Lake Constance basin, Central Europe. Mating trials exposing females to pairings of territorial lake and stream males in outdoor mesocosms allowing for natural reproductive behaviour reveal that mating occurs preferentially between partners of the same ecotype. Compared to random mating, this sexual barrier reduces gene flow between the ecotypes by some 36%. This relatively modest strength of sexual isolation is surprising because comparing the males between the two ecotypes shows striking differentiation in traits generally considered relevant to reproductive behaviour (body size, breeding coloration, nest size). Analysing size differences among the individuals in the mating trials further indicates that assortative mating is not related to ecotype differences in body size. Overall, we demonstrate that sexually based reproductive isolation promotes divergence in lake–stream stickleback along with other known reproductive barriers, but we also caution against inferring strong sexual isolation from the observation of strong population divergence in sexually relevant traits.     

Detecting Holocene divergence in the anadromous-freshwater three-spined stickleback (Gasterosteus aculeatus) system

The anadromous-freshwater three-spined stickleback (Gasterosteus aculeatus) system allows for inferring the role of adaptation in speciation with a high level of accuracy because the freshwater ecotype has evolved multiple times from a uniform anadromous ancestor. A cause for concern is that independent evolution among drainages is not guaranteed in areas with a poorly resolved glacial history. This is the case for the west European great rivers, whose downstream valleys flanked the southern limit of the late Pleistocene ice sheet. We tested for independent and postglacial colonization of these valleys hypothesizing that the relationships among anadromous and freshwater sticklebacks correspond to a raceme structure. We compared the reduction in plate number accompanying this colonization to the genetic differentiation using 13 allozyme and five microsatellite loci in 350 individuals. Overall microsatellite differentiation (F(ST) = 0.147) was twice as large as allozyme differentiation (F(ST) = 0.066). Although habitat-specific gene flow may mask the ancestral relationships among both ecotypes, levels of microsatellite differentiation supported the hypothesis of raceme-like divergence, reflecting independent colonizations rather than the presence of two distinct evolutionary clades. Under an infinite alleles model and in the absence of gene flow, the observed freshwater divergence might be reached after 440 (microsatellites) to 4500 (allozymes) generations. Hence, the anadromous-freshwater stickleback system most likely diverged postglacially. We conclude that the reduction in plate number in two freshwater basins probably occurred independently, and that its considerable variation among populations is not in agreement with the time since divergence.     

Evolution and stability of the G-matrix during the colonization of a novel environment

Populations that undergo a process of rapid evolution present excellent opportunities to investigate the mechanisms driving or restraining adaptive divergence. The genetic variance-covariance matrix (G) is often considered to constrain adaptation but little is known about its potential to evolve during phenotypic divergence. We compared the G-matrices of ancestral and recently established ecotype populations of an aquatic isopod (Asellus aquaticus) that have diverged in parallel in two south Swedish lakes. Phenotypic changes after colonization involved a reduction in overall size, lost pigmentation and changes in shape. Comparisons between G-matrices reveal close similarity within the same ecotype from different lakes but some degree of differentiation among ecotypes. Phenotypic divergence has apparently not been much influenced by the orientation of G. Additive genetic variation in the newly colonized habitats has also decreased substantially. This suggests that a process of adaptation from standing genetic variation has occurred and has probably facilitated phenotypic divergence.     

Chromosomal Fusions Facilitate Adaptation to Divergent Environments in Threespine Stickleback

Chromosomal fusions are hypothesized to facilitate adaptation to divergent environments, both by bringing together previously unlinked adaptive alleles and by creating regions of low recombination that facilitate the linkage of adaptive alleles; but, there is little empirical evidence to support this hypothesis. Here, we address this knowledge gap by studying threespine stickleback (Gasterosteus aculeatus), in which ancestral marine fish have repeatedly adapted to freshwater across the northern hemisphere. By comparing the threespine and ninespine stickleback (Pungitius pungitius) genomes to a de novo assembly of the fourspine stickleback (Apeltes quadracus) and an outgroup species, we find two chromosomal fusion events involving the same chromosomes have occurred independently in the threespine and ninespine stickleback lineages. On the fused chromosomes in threespine stickleback, we find an enrichment of quantitative trait loci underlying traits that contribute to marine versus freshwater adaptation. By comparing whole-genome sequences of freshwater and marine threespine stickleback populations, we also find an enrichment of regions under divergent selection on these two fused chromosomes. There is elevated genetic diversity within regions under selection in the freshwater population, consistent with a simulation study showing that gene flow can increase diversity in genomic regions associated with local adaptation and our demographic models showing gene flow between the marine and freshwater populations. Integrating our results with previous studies, we propose that these fusions created regions of low recombination that enabled the formation of adaptative clusters, thereby facilitating freshwater adaptation in the face of recurrent gene flow between marine and freshwater threespine sticklebacks.