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.     

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.     

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.     

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.     

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.     

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.     

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.     

Microgeographical diversification of threespine stickleback: body shape–habitat correlations in a small, ecologically diverse Alaskan drainage

Adaptive radiations are a major source of evolutionary diversity in nature, and understanding how they originate and how organisms diversify during the early stages of adaptive radiation is a major problem in evolutionary biology. The relationship between habitat type and body shape variation was investigated in a postglacial radiation of threespine stickleback in the upper Fish Creek drainage of Cook Inlet, Alaska. Although small, the upper Fish Creek drainage includes ecologically diverse lakes and streams in close proximity to one another that harbour abundant stickleback. Specimens from ancestral anadromous and derived resident freshwater populations differed substantially and could be distinguished by body shape alone, suggesting that the initial stages of adaptation contribute disproportionately to evolutionary divergence. Body shape divergence among resident freshwater populations was also considerable, and phenotypic distances among samples from freshwater populations were associated with habitat type but not geographical distance. As expected, stream stickleback from slow-moving, structurally complex environments tended to have the deepest bodies, stickleback from lakes with a mostly benthic habitat were similar but less extreme, and stickleback from lakes with a mostly limnetic habitat were the most shallow-bodied, elongate fish. Beyond adapting rapidly to conditions in freshwater environments, stickleback can diversify rapidly over small geographical scales in freshwater systems despite opportunities for gene flow. This study highlights the importance of ecological heterogeneity over small geographical scales for evolutionary diversification during the early stages of adaptive radiation, and lays the foundation for future research on this ecologically diverse, postglacial system.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 98 , 139–151.     

Evidence of Adaptive Evolutionary Divergence during Biological Invasion

Rapid phenotypic diversification during biological invasions can either arise by adaptation to alternative environments or by adaptive phenotypic plasticity. Where experimental evidence for adaptive plasticity is common, support for evolutionary diversification is rare. Here, we performed a controlled laboratory experiment using full-sib crosses between ecologically divergent threespine stickleback populations to test for a genetic basis of adaptation. Our populations are from two very different habitats, lake and stream, of a recently invaded range in Switzerland and differ in ecologically relevant morphological traits. We found that in a lake-like food treatment lake fish grow faster than stream fish, resembling the difference among wild type individuals. In contrast, in a stream-like food treatment individuals from both populations grow similarly. Our experimental data suggest that genetically determined diversification has occurred within less than 140 years after the arrival of stickleback in our studied region.     

Allometric relationships in morphological traits associated with foraging,swimming ability,and predator defense reveal adaptations toward brackish and freshwater environments in the threespine stickleback

Freshwater colonization by threespine stickleback has led to divergence in morphology between ancestral marine and derived freshwater populations, making them ideal for studying natural selection on phenotypes. In an open brackish–freshwater system, we previously discovered two genetically distinct stickleback populations that also differ in geometric shape: one mainly found in the brackish water lagoon and one throughout the freshwater system. As shape and size are not perfectly correlated, the aim of this study was to identify the morphological trait(s) that separated the populations in geometric shape. We measured 23 phenotypes likely to be important for foraging, swimming capacity, and defense against predation. The lateral plate morphs in freshwater displayed few significant changes in trait sizes, but the low plated expressed feeding traits more associated with benthic habitats. When comparing the completely plated genetically assigned populations, the freshwater, the hybrids, the migrants and the lagoon fish, many of the linear traits had different slopes and intercepts in trait‐size regressions, precluding our ability to directly compare all traits simultaneously, which most likely results from low variation in body length for the lagoon and migrant population. We found the lagoon stickleback population to be more specialized toward the littoral zone, displaying benthic traits such as large, deep bodies with smaller eyes compared to the freshwater completely plated morph. Further, the lagoon and migrant fish had an overall higher body coverage of lateral plates compared to freshwater fish, and the dorsal and pelvic spines were longer. Evolutionary constraints due to allometric scaling relationships could explain the observed, overall restricted, differences in morphology between the sticklebacks in this study, as most traits have diversified in common allometric trajectories. The observed differences in foraging and antipredation traits between the fish with a lagoon and freshwater genetic signature are likely a result of genetic or plastic adaptations toward brackish and freshwater environments.     

Life-history plasticity in female threespine stickleback

The postglacial adaptive radiation of the threespine stickleback fish (Gasterosteus aculeatus) has been widely used to investigate the roles of both adaptive evolution and plasticity in behavioral and morphological divergence from the ancestral condition represented by present-day oceanic stickleback. These phenotypes tend to exhibit high levels of ecotypic differentiation. Population divergence in life history has also been well studied, but in contrast to behavior and morphology, the extent and importance of plasticity has been much less well studied. In this review, we summarize what is known about life-history plasticity in female threespine stickleback, considering four traits intimately associated with reproductive output: age/size at maturation, level of reproductive effort, egg size and clutch size. We envision life-history plasticity in an iterative, ontogenetic framework, in which females may express plasticity repeatedly across each of several time frames. We contrast the results of laboratory and field studies because, for most traits, these approaches give somewhat different answers. We provide ideas on what the cues might be for observed plasticity in each trait and, when possible, we inquire about the relative costs and benefits to expressed plasticity. We end with an example of how we think plasticity may play out in stickleback life history given what we know of plasticity in the ancestor.     

Cryptic genetic variation and body size evolution in threespine stickleback

The role of environment as a selective agent is well-established. Environment might also influence evolution by altering the expression of genetic variation associated with phenotypes under selection. Far less is known about this phenomenon, particularly its contribution to evolution in novel environments. We investigated how environment affected the evolvability of body size in the threespine stickleback (Gasterosteus aculeatus). Gasterosteus aculeatus is well suited to addressing this question due to the rapid evolution of smaller size in the numerous freshwater populations established following the colonization of new freshwater habitats by an oceanic ancestor. The repeated, rapid evolution of size following colonization contrasts with the general observation of low phenotypic variation in oceanic stickleback. We reared an oceanic population of stickleback under high and low salinity conditions, mimicking a key component of the ancestral environment, and freshwater colonization, respectively. There was low genetic variation for body size under high salinity, but this variance increased significantly when fish were reared under low salinity. We therefore conclude that oceanic populations harbor the standing genetic variation necessary for the evolution of body size, but that this variation only becomes available to selection upon colonization of a new habitat.     

Two morphotypes of lacustrine threespine stickleback,Gasterosteus aculeatus, in Benka Lake, Alaska

Synopsis Within freshwater fishes, a common pattern of diversification of body form and trophic structure has recently been recognized. Two different suites of co-occurring characters appear to allow fish to efficiently forage on either benthic invertebrates or plankton. For threespine stickleback,Gasterosteus aculeatus, these suites of characters have been labeled benthic and limnetic. The forms differ in several morpholoical traits, with the limnetic having a more fusiform shape, larger eyes, longer and more numerous gill rakers, and a smaller more tubular mouth. Benthic and limnetic threespine stickleback are usually found in allopatry or parapatry, and less frequently in sympatry. Within the range of the threespine stickleback, which comprises perhaps tens of thousands of lacustrine populations, the sympatric occurrence of the benthic and limnetic forms has only been established for six lakes within a small region of the Strait of Georgia, British Columbia. We present the first evidence for the presence of sympatric morphotypes of threespine stickleback outside of British Columbia. We examine the nature and extent of this diversification, and present possible explanations for the sympatric occurrence of these morphotypes. We also explore possible reasons for the small number of documented sympatric benthic and limnetic forms of threespine stickleback despite the existence of thousands of apparently suitable lakes.     

Quantitative genetic variation in static allometry in the threespine stickleback

The common pattern of replicated evolution of a consistent shape-environment relationship might reflect selection acting in similar ways within each environment, but divergently among environments. However, phenotypic evolution depends on the availability of additive genetic variation as well as on the direction of selection, implicating a bias in the distribution of genetic variance as a potential contributor to replicated evolution. Allometry, the relationship between shape and size, is a potential source of genetic bias that is poorly understood. The threespine stickleback, Gasterosteus aculeatus, provides an ideal system for exploring the contribution of genetic variance in body shape allometry to evolutionary patterns. The stickleback system comprises marine populations that exhibit limited phenotypic variation, and young freshwater populations which, following independent colonization events, have often evolved similar phenotypes in similar environments. In particular, stickleback diversification has involved changes in both total body size and relative size of body regions (i.e., shape). In a laboratory-reared cohort derived from an oceanic Alaskan population that is phenotypically and genetically representative of the ancestor of the diverse freshwater populations in this region, we determined the phenotypic static allometry, and estimated the additive genetic variation about these population-level allometric functions. We detected significant allometry, with larger fish having relatively smaller heads, a longer base to their second dorsal fin, and longer, shallower caudal peduncles. There was additive genetic variance in body size and in size-independent body shape (i.e., allometric elevation), but typically not in allometric slopes. These results suggest that the parallel evolution of body shape in threespine stickleback is not likely to have been a correlated response to selection on body size, or vice versa. Although allometry is common in fishes, this study highlights the need for additional data on genetic variation in allometric functions to determine how allometry evolves and how it influences phenotypic evolution.     

Effects of genetics and light environment on colour expression in threespine sticklebacks

The genetic basis of traits that are under sexual selection and that are involved in recognizing conspecific mates is poorly known, even in systems in which the phenotypic basis of these traits has been well studied. In the present study, we investigate genetic and environmental influences on nuptial colour, which plays important roles in sexual selection and sexual isolation in species pairs of limnetic and benthic threespine sticklebacks ( Gasterosteus aculeatus species complex). Previous work demonstrated that colour differences among species correlate to differences in the ambient light prevalent in their mating habitat. Red fish are found in clear water and black fish in red-shifted habitats. We used a paternal half-sib split-clutch design to investigate the genetic and environmental basis of nuptial colour. We found genetic differences between a red and a black stickleback population in the expression of both red and black nuptial colour. In addition, the light environment influenced colour expression, and genotype by environment interactions were also present. We found evidence for both phenotypic and genetic correlations between our colour traits; some of these correlations are in opposite directions for our red and black populations. These results suggest that both genetic change and phenotypic plasticity underlie the correlation of male colour with light environment.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 94 , 663–673.     

Independent axes of genetic variation and parallel evolutionary divergence of opercle bone shape in threespine stickleback

Evolution of similar phenotypes in independent populations is often taken as evidence of adaptation to the same fitness optimum. However, the genetic architecture of traits might cause evolution to proceed more often toward particular phenotypes, and less often toward others, independently of the adaptive value of the traits. Freshwater populations of Alaskan threespine stickleback have repeatedly evolved the same distinctive opercle shape after divergence from an oceanic ancestor. Here we demonstrate that this pattern of parallel evolution is widespread, distinguishing oceanic and freshwater populations across the Pacific Coast of North America and Iceland. We test whether this parallel evolution reflects genetic bias by estimating the additive genetic variance-covariance matrix (G) of opercle shape in an Alaskan oceanic (putative ancestral) population. We find significant additive genetic variance for opercle shape and that G has the potential to be biasing, because of the existence of regions of phenotypic space with low additive genetic variation. However, evolution did not occur along major eigenvectors of G, rather it occurred repeatedly in the same directions of high evolvability. We conclude that the parallel opercle evolution is most likely due to selection during adaptation to freshwater habitats, rather than due to biasing effects of opercle genetic architecture.     

Female mate preferences for male body size and shape promote sexual isolation in threespine sticklebacks

Female mate preferences for ecologically relevant traits may enhance natural selection, leading to rapid divergence. They may also forge a link between mate choice within species and sexual isolation between species. Here, we examine female mate preference for two ecologically important traits: body size and body shape. We measured female preferences within and between species of benthic, limnetic, and anadromous threespine sticklebacks (Gasterosteus aculeatus species complex). We found that mate preferences differed between species and between contexts (i.e., within vs. between species). Within species, anadromous females preferred males that were deep bodied for their size, benthic females preferred larger males (as measured by centroid size), and limnetic females preferred males that were more limnetic shaped. In heterospecific mating trials between benthics and limnetics, limnetic females continued to prefer males that were more limnetic like in shape when presented with benthic males. Benthic females showed no preferences for size when presented with limnetic males. These results show that females use ecologically relevant traits to select mates in all three species and that female preference has diverged between species. These results suggest that sexual selection may act in concert with natural selection on stickleback size and shape. Further, our results suggest that female preferences may track adaptation to local environments and contribute to sexual isolation between benthic and limnetic sticklebacks.     

Classification of threespine stickleback along the benthic-limnetic axis

Many species of fish display morphological divergence between individuals feeding on macroinvertebrates associated with littoral habitats (benthic morphotypes) and individuals feeding on zooplankton in the limnetic zone (limnetic morphotypes). Threespine stickleback (Gasterosteus aculeatus L.) have diverged along the benthic-limnetic axis into allopatric morphotypes in thousands of populations and into sympatric species pairs in several lakes. However, only a few well known populations have been studied because identifying additional populations as either benthic or limnetic requires detailed dietary or observational studies. Here we develop a Fisher's linear discriminant function based on the skull morphology of known benthic and limnetic stickleback populations from the Cook Inlet Basin of Alaska and test the feasibility of using this function to identify other morphologically divergent populations. Benthic and limnetic morphotypes were separable using this technique and of 45 populations classified, three were identified as morphologically extreme (two benthic and one limnetic), nine as moderately divergent (three benthic and six limnetic) and the remaining 33 populations as morphologically intermediate. Classification scores were found to correlate with eye size, the depth profile of lakes, and the presence of invasive northern pike (Esox lucius). This type of classification function provides a means of integrating the complex morphological differences between morphotypes into a single score that reflects the position of a population along the benthic-limnetic axis and can be used to relate that position to other aspects of stickleback biology.     

Strong and parallel salinity‐induced phenotypic plasticity in one generation of threespine stickleback

Phenotypic plasticity is a major factor contributing to variation of organisms in nature, yet its evolutionary significance is insufficiently understood. One example system where plasticity might have played an important role in an adaptive radiation is the threespine stickleback (Gasterosteus aculeatus), a fish that has diversified after invading freshwater lakes repeatedly from the marine habitat. The parallel phenotypic changes that occurred in this radiation were extremely rapid. This study evaluates phenotypic plasticity in stickleback body shape in response to salinity in fish stemming from a wild freshwater population. Using a split‐clutch design, we detected surprisingly large phenotypically plastic changes in body shape after one generation. Fish raised in salt water developed shallower bodies and longer jaws, and these changes were consistent and parallel across families. Although this work highlights the effect of phenotypic plasticity, we also find indications that constraints may play a role in biasing the direction of possible phenotypic change. The slopes of the allometric relationship of individual linear traits did not change across treatments, indicating that plastic change does not affect the covariation of traits with overall size. We conclude that stickleback have a large capacity for plastic phenotypic change in response to salinity and that plasticity and evolutionary constraints have likely contributed to the phenotypic diversification of these fish.     

Implications of a fossil stickleback assemblage for Darwinian gradualism

Darwin postulated that a complete fossil record would contain numerous gradual transitions between ancestral and descendant species, but 150 years after publication of The Origin of Species, few such transitions have materialized. The fossil stickleback Gasterosteus doryssus and the deposit in which it occurs provide excellent conditions to detect such transitions. Abundant, well‐preserved fossils occur in a stratigraphic setting with fine temporal resolution. The paleoecology of G. doryssus resembles the ecology of modern lakes that harbour the phenotypically similar three‐spined stickleback Gasterosteus aculeatus. Gasterosteus aculeatus are primitively highly armoured, but G. doryssus comprised two contemporaneous biological species with relatively weak armour, including a near‐shore, benthic feeder (benthic) and an offshore planktivore (limnetic). The benthic species expanded its range into the limnetic zone of the lake, where it apparently switched to planktivory and evolved reduced armour within c. 5000 years in response to directional selection. Although gradual evolution of mean phenotypes occurred, a single major gene caused much of evolutionary change of the pelvic skeleton. Thus, Darwin's expectation that transitions between species in the fossil record would be gradual was met at a fine time scale, but for pelvic structure, a well‐studied trait, his expectation that gradual change would depend entirely on numerous, small, heritable differences among individuals was incorrect.