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.     

Ancestral Plasticity and Allometry in Threespine Stickleback Fish Reveal Phenotypes Associated with Derived, Freshwater Ecotypes

For over a century, evolutionary biologists have debated whether and how phenotypic plasticity impacts the processes of adaptation and diversification. The empirical tests required to resolve these issues have proven elusive, mainly because it requires documentation of ancestral reaction norms, a difficult prospect where many ancestors are either extinct or have evolved. The threespine stickleback radiation is not limited in this regard, making it an ideal system in which to address general questions regarding the role of plasticity in adaptive evolution. As retreating ice sheets have exposed new habitats, oceanic stickleback founded innumerable freshwater populations, many of which have evolved parallel adaptations to their new environments. Because the founding oceanic population is extant, we can directly evaluate whether specific patterns of ancestral phenotypic expression in the context of novel environments (plasticity), or over ontogeny, predisposed the repeated evolution of "benthic" and "limnetic" ecotypes in shallow and deep lakes, respectively. Consistent with this hypothesis, we found that oceanic stickleback raised in a complex habitat and fed a macroinvertebrate diet expressed traits resembling derived, benthic fish. Alternatively, when reared in a simple environment on a diet of zooplankton, oceanic stickleback developed phenotypes resembling derived, limnetic fish. As fish in both treatments grew, their body depths increased allometrically, as did the size of their mouths, while their eyes became relatively smaller. Allometric trajectories were subtly but significantly impacted by rearing environment. Thus, both environmental and allometric influences on development, along with their interactive effects, produced variation in phenotypes consistent with derived benthic and limnetic fish, which may have predisposed the repeated genetic accommodation of this specific suite of traits. We also found significant shape differences between marine and anadromous stickleback, which has implications for evaluating the ancestral state of stickleback traits.     

Adaptive plasticity and genetic divergence in feeding efficiency during parallel adaptive radiation of whitefish (Coregonus spp.)

Parallel phenotypic divergence in replicated adaptive radiations could either result from parallel genetic divergence in response to similar divergent selection regimes or from equivalent phenotypically plastic response to the repeated occurrence of contrasting environments. In post‐glacial fish, replicated divergence in phenotypes along the benthic‐limnetic habitat axis is commonly observed. Here, we use two benthic‐limnetic species pairs of whitefish from two Swiss lakes, raised in a common garden design, with reciprocal food treatments in one species pair, to experimentally measure whether feeding efficiency on benthic prey has a genetic basis or whether it underlies phenotypic plasticity (or both). To do so, we offered experimental fish mosquito larvae, partially burried in sand, and measured multiple feeding efficiency variables. Our results reveal both, genetic divergence as well as phenotypically plastic divergence in feeding efficiency, with the phenotypically benthic species raised on benthic food being the most efficient forager on benthic prey. This indicates that both, divergent natural selection on genetically heritable traits and adaptive phenotypic plasticity, are likely important mechanisms driving phenotypic divergence in adaptive radiation.     

Responses to simulated winter conditions differ between threespine stickleback ecotypes

Abiotic factors can act as barriers to colonization and drive local adaptation. During colonization, organisms may cope with changes in abiotic factors using existing phenotypic plasticity, but the role of phenotypic plasticity in assisting or hindering the process of local adaptation remains unclear. To address these questions, we explore the role of winter conditions in driving divergence during freshwater colonization and the effects of plasticity on local adaptation in ancestral marine and derived freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We found that freshwater‐resident stickleback had greater tolerance of acute exposure to low temperatures than marine stickleback, but these differences were abolished after acclimation to simulated winter conditions (9L:15D photoperiod at 4 °C). Plasma chloride levels differed between the ecotypes, but showed a similar degree of plasticity between ecotypes. Gene expression of the epithelial calcium channel (ECaC) differed between ecotypes, with the freshwater ecotype demonstrating substantially greater expression than the marine ecotype, but there was no plasticity in this trait under these conditions in either ecotype. In contrast, growth (assessed as final mass) and the expression of an isoform of the electroneutral Na⁺/H⁺ exchanger (NHE3) exhibited substantial change with temperature in the marine ecotype that was not observed in the freshwater ecotype under the conditions tested here, which is consistent with evolution of these traits by a process such as genetic assimilation. These data demonstrate substantial divergence in many of these traits between freshwater and marine stickleback, but also illustrate the complexity of possible relationships between plasticity and local adaptation.     

Adaptation and acclimation of aerobic exercise physiology in Lake Whitefish ecotypes (Coregonus clupeaformis)

The physiological mechanisms underlying local adaptation in natural populations of animals, and whether the same mechanisms contribute to adaptation and acclimation, are largely unknown. Therefore, we tested for evolutionary divergence in aerobic exercise physiology in laboratory bred, size‐matched crosses of ancestral, benthic, normal Lake Whitefish (Coregonus clupeaformis) and derived, limnetic, more actively swimming “dwarf” ecotypes. We acclimated fish to constant swimming (emulating limnetic foraging) and control conditions (emulating normal activity levels) to simultaneously study phenotypic plasticity. We found extensive divergence between ecotypes: dwarf fish generally had constitutively higher values of traits related to oxygen transport (ventricle size) and use by skeletal muscle (percent oxidative muscle, mitochondrial content), and also evolved differential plasticity of mitochondrial function (Complex I activity and flux through Complexes I–IV and IV). The effects of swim training were less pronounced than differences among ecotypes and the traits which had a significant training effect (ventricle protein content, ventricle malate dehydrogenase activity, and muscle Complex V activity) did not differ among ecotypes. Only one trait, ventricle mass, varied in a similar manner with acclimation and adaptation and followed a pattern consistent with genetic accommodation. Overall, the physiological and biochemical mechanisms underlying acclimation and adaptation to swimming activity in Lake Whitefish differ.     

Niche specialization influences adaptive phenotypic plasticity in the threespine stickleback

Phenotypic plasticity may be favored in generalist populations if it increases niche width, even in temporally constant environments. Phenotypic plasticity can increase the frequency of extreme phenotypes in a population and thus allow it to make use of a wide resource spectrum. Here we test the prediction that generalist populations should be more plastic than specialists. In a common-garden experiment, we show that solitary, generalist populations of threespine sticklebacks inhabiting small coastal lakes of British Columbia have a higher degree of morphological plasticity than the more specialized sympatric limnetic and benthic species. The ancestral marine stickleback showed low levels of plasticity similar to those of sympatric sticklebacks, implying that the greater plasticity of the generalist population has evolved recently. Measurements of wild populations show that those with mean trait values intermediate between the benthic and limnetic values indeed have higher morphological variation. Our data indicate that plasticity can evolve rapidly after colonization of a new environment in response to changing niche use.     

Partially repeatable genetic basis of benthic adaptation in threespine sticklebacks

The extent to which convergent adaptation to similar ecological niches occurs by a predictable genetic basis remains a fundamental question in biology. Threespine stickleback fish have undergone an adaptive radiation in which ancestral oceanic populations repeatedly colonized and adapted to freshwater habitats. In multiple lakes in British Columbia, two different freshwater ecotypes have evolved: a deep‐bodied benthic form adapted to forage near the lake substrate, and a narrow‐bodied limnetic form adapted to forage in open water. Here, we use genome‐wide linkage mapping in marine × benthic F2 genetic crosses to test the extent of shared genomic regions underlying benthic adaptation in three benthic populations. We identify at least 100 Quantitative Trait Loci (QTL) harboring genes influencing skeletal morphology. The majority of QTL (57%) are unique to one cross. However, four genomic regions affecting eight craniofacial and armor phenotypes are found in all three benthic populations. We find that QTL are clustered in the genome and overlapping QTL regions are enriched for genomic signatures of natural selection. These findings suggest that benthic adaptation has occurred via both parallel and nonparallel genetic changes.     

Parallel genetic divergence among coastal–marine ecotype pairs of European anchovy explained by differential introgression after secondary contact

Ecophenotypic differentiation among replicate ecotype pairs within a species complex is often attributed to independent outcomes of parallel divergence driven by adaptation to similar environmental contrasts. However, the extent to which parallel phenotypic and genetic divergence patterns have emerged independently is increasingly questioned by population genomic studies. Here, we document the extent of genetic differentiation within and among two geographic replicates of the coastal and marine ecotypes of the European anchovy (Engraulis encrasicolus) gathered from Atlantic and Mediterranean locations. Using a genome‐wide data set of RAD‐derived SNPs, we show that habitat type (marine vs. coastal) is the most important component of genetic differentiation among populations of anchovy. By analysing the joint allele frequency spectrum of each coastal–marine ecotype pair, we show that genomic divergence patterns between ecotypes can be explained by a postglacial secondary contact following a long period of allopatric isolation (c. 300 kyrs). We found strong support for a model including heterogeneous migration among loci, suggesting that secondary gene flow has eroded past differentiation at different rates across the genome. Markers experiencing reduced introgression exhibited strongly correlated differentiation levels among Atlantic and Mediterranean regions. These results support that partial reproductive isolation and parallel genetic differentiation among replicate pairs of anchovy ecotypes are largely due to a common divergence history prior to secondary contact. They moreover provide comprehensive insights into the origin of a surprisingly strong fine‐scale genetic structuring in a high gene flow marine fish, which should improve stock management and conservation actions.     

Parallel evolution of lake whitefish dwarf ecotypes in association with limnological features of their adaptive landscape

Sympatric fish populations observed in many north temperate lakes are among the best models to study the processes of population divergence and adaptive radiation. Despite considerable research on such systems, little is known about the associations between ecological conditions and the extent of ecotypic divergence. In this study, we examined the biotic and abiotic properties of postglacial lakes in which lake whitefish, Coregonus clupeaformis, occur as a derived dwarf ecotype in sympatry with an ancestral normal ecotype. We compared 19 limnological variables between two groups of lakes known from previous studies to harbour sympatric dwarf and normal ecotypes with high and low levels of phenotypic and genetic differentiation respectively. We found clear environmental differences between the two lake groups. Namely, oxygen was the most discriminant variable, where lakes harbouring the most divergent populations were characterized by the greatest hypolimnetic oxygen depletion. These lakes also had lower zooplankton densities and a narrower distribution of zooplantonic prey length. These results suggest that the highest differentiation between sympatric ecotypes occurs in lakes with reduced habitat and prey availability that could increase competition for resources. This in turns supports the hypothesis that parallelism in the extent of phenotypic divergence among sympatric whitefish ecotypes is associated with parallelism in adaptive landscape in terms of differences in limnological characteristics, as well as availability and structure of the zooplanktonic community.     

Phenotypic Plasticity in Response to the Social Environment: Effects of Density and Sex Ratio on Mating Behaviour Following Ecotype Divergence

The ability to express phenotypically plastic responses to environmental cues might be adaptive in changing environments. We studied phenotypic plasticity in mating behaviour as a response to population density and adult sex ratio in a freshwater isopod (Asellus aquaticus). A. aquaticus has recently diverged into two distinct ecotypes, inhabiting different lake habitats (reed Phragmites australis and stonewort Chara tomentosa, respectively). In field surveys, we found that these habitats differ markedly in isopod population densities and adult sex ratios. These spatially and temporally demographic differences are likely to affect mating behaviour. We performed behavioural experiments using animals from both the ancestral ecotype (“reed” isopods) and from the novel ecotype (“stonewort” isopods) population. We found that neither ecotype adjusted their behaviour in response to population density. However, the reed ecotype had a higher intrinsic mating propensity across densities. In contrast to the effects of density, we found ecotype differences in plasticity in response to sex ratio. The stonewort ecotype show pronounced phenotypic plasticity in mating propensity to adult sex ratio, whereas the reed ecotype showed a more canalised behaviour with respect to this demographic factor. We suggest that the lower overall mating propensity and the phenotypic plasticity in response to sex ratio have evolved in the novel stonewort ecotype following invasion of the novel habitat. Plasticity in mating behaviour may in turn have effects on the direction and intensity of sexual selection in the stonewort habitat, which may fuel further ecotype divergence.     

Plastic responses to parents and predators lead to divergent shoaling behaviour in sticklebacks

Population divergence in antipredator defence and behaviour occurs rapidly and repeatedly. Genetic differences, phenotypic plasticity or parental effects may all contribute to divergence, but the relative importance of each of these mechanisms remains unknown. We exposed juveniles to parents and predators to measure how induced changes contribute to shoaling behaviour differences between two threespine stickleback species (benthics and limnetics: Gasterosteus spp). We found that limnetics increased shoaling in response to predator attacks, whereas benthics did not alter their behaviour. Care by limnetic fathers led to increased shoaling in both limnetic and benthic offspring. Shoaling helps limnetics avoid trout and avian predation; our results suggest that this adaptive behaviour is the result of a combination of paternal effects, predator-induced plasticity and genetic differences between species. These results suggest that plasticity substantially contributes to the rapid divergence in shoaling behaviour across the post-Pleistocene radiation of sticklebacks.     

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.     

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.     

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.     

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild‐caught individuals. Here, we investigated the contribution of habitat‐specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake–river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake‐like pond set‐up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set‐up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set‐up were based on a plastic response. Finally, gene expression and bacterial communities of wild‐caught individuals and individuals acclimatized to lake‐like pond conditions showed shifts underlying adaptive phenotypic plasticity.     

Divergence across diet,time and populations rules out parallel evolution in the gut microbiomes of Trinidadian guppies

Diverse microbial consortia profoundly influence animal biology, necessitating an understanding of microbiome variation in studies of animal adaptation. Yet, little is known about such variability among fish, in spite of their importance in aquatic ecosystems. The Trinidadian guppy, Poecilia reticulata, is an intriguing candidate to test microbiome-related hypotheses on the drivers and consequences of animal adaptation, given the recent parallel origins of a similar ecotype across streams. To assess the relationships between the microbiome and host adaptation, we used 16S rRNA amplicon sequencing to characterize gut bacteria of two guppy ecotypes with known divergence in diet, life history, physiology and morphology collected from low-predation (LP) and high-predation (HP) habitats in four Trinidadian streams. Guts were populated by several recurring, core bacteria that are related to other fish associates and rarely detected in the environment. Although gut communities of lab-reared guppies differed from those in the wild, microbiome divergence between ecotypes from the same stream was evident under identical rearing conditions, suggesting host genetic divergence can affect associations with gut bacteria. In the field, gut communities varied over time, across streams and between ecotypes in a stream-specific manner. This latter finding, along with PICRUSt predictions of metagenome function, argues against strong parallelism of the gut microbiome in association with LP ecotype evolution. Thus, bacteria cannot be invoked in facilitating the heightened reliance of LP guppies on lower-quality diets. We argue that the macroevolutionary microbiome convergence seen across animals with similar diets may be a signature of secondary microbial shifts arising some time after host-driven adaptation.     

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.     

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat‐related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model‐based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.     

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.