Genetic and morphometric divergence in threespine stickleback in the Chignik catchment,Alaska

Divergent selection pressures induced by different environmental conditions typically lead to variation in life history, behavior, and morphology. When populations are locally adapted to their current environment, selection may limit movement into novel sites, leading to neutral and adaptive genetic divergence in allopatric populations. Subsequently, divergence can be reinforced by development of pre‐ or postzygotic barriers to gene flow. The threespine stickleback, Gasterosteus aculeatus, is a primarily marine fish that has invaded freshwater repeatedly in postglacial times. After invasion, the established freshwater populations typically show rapid diversification of several traits as they become reproductively isolated from their ancestral marine population. In this study, we examine the genetic and morphometric differentiation between sticklebacks living in an open system comprising a brackish water lagoon, two freshwater lakes, and connecting rivers. By applying a set of microsatellite markers, we disentangled the genetic relationship of the individuals across the diverse environments and identified two genetic populations: one associated with brackish and the other with the freshwater environments. The “brackish” sticklebacks were larger and had a different body shape than those in freshwater. However, we found evidence for upstream migration from the brackish lagoon into the freshwater environments, as fish that were genetically and morphometrically similar to the lagoon fish were found in all freshwater sampling sites. Regardless, few F1‐hybrids were identified, and it therefore appears that some pre‐ and/or postzygotic barriers to gene flow rather than geographic distance are causing the divergence in this system.     

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.     

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.     

Stoichiometric traits of stickleback: Effects of genetic background,rearing environment,and ontogeny

Phenotypes can both evolve in response to, and affect, ecosystem change, but few examples of diverging ecosystem‐effect traits have been investigated. Bony armor traits of fish are good candidates for this because they evolve rapidly in some freshwater fish populations, and bone is phosphorus rich and likely to affect nutrient recycling in aquatic ecosystems. Here, we explore how ontogeny, rearing environment, and bone allocation among body parts affect the stoichiometric phenotype (i.e., stoichiometric composition of bodies and excretion) of threespine stickleback. We use two populations from distinct freshwater lineages with contrasting lateral plating phenotypes (full vs. low plating) and their hybrids, which are mostly fully plated. We found that ontogeny, rearing environment, and body condition were the most important predictors of organismal stoichiometry. Although elemental composition was similar between both populations and their hybrids, we found significant divergence in phosphorus allocation among body parts and in phosphorus excretion rates. Overall, body armor differences did not explain variation in whole body phosphorus, phosphorus allocation, or phosphorus excretion. Evolutionary divergence between these lineages in both allocation and excretion is likely to have important direct consequences for ecosystems, but may be mediated by evolution of multiple morphological or physiological traits beyond plating phenotype.     

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.     

Genomic signatures of the plateless phenotype in the threespine stickleback

Understanding the genetic basis of traits involved in adaptive divergence and speciation is one of the most fundamental objectives in evolutionary biology. Toward that end, we look for signatures of extreme plate loss in the genome of freshwater threespine sticklebacks (Gasterosteus aculeatus). Plateless stickleback have been found in only a few lakes and streams across the world; they represent the far extreme of a phenotypic continuum (plate number) that has been studied for years, although plateless individuals have not yet been the subject of much investigation. We use a dense single nucleotide polymorphism dataset made using RADseq to study fish from three freshwater populations containing plateless and low plated individuals, as well as fish from full plated marine populations. Analyses were performed using FastStructure, sliding windows F_ST, Bayescan and latent factor mixed models to search for genomic differences between the low plated and plateless phenotypes both within and among the three lakes. At least 18 genomic regions which may contribute to within‐morph plate number variation were detected in our low plated stickleback populations. We see no evidence of a selective sweep between low and plateless fish; rather reduction of plate number within the low plated morph seems to be polygenic.     

Variation of telencephalon morphology of the threespine stickleback (Gasterosteus aculeatus) in relation to inferred ecology

We tested the hypothesis that increased telencephalon size has evolved in threespine stickleback fish (Gasterosteus aculeatus) from structurally complex habitats using field‐caught samples from one sea‐run (ancestral) and 18 ecologically diverse freshwater (descendant) populations. Freshwater habitats ranged from shallow, structurally complex lakes with benthic‐foraging stickleback (benthics), to deeper, structurally simple lakes in which stickleback depend more heavily on plankton for prey (generalists). Contrary to our expectations, benthics had smaller telencephala than generalists, but the shape of the telencephalon of the sea‐run and benthic populations were more convex laterally. Convex telencephalon shape may indicate enlargement of the dorsolateral region, which is homologous with the tetrapod hippocampus. Telencephalon morphology is also sexually dimorphic, with larger, less convex telencephala in males. Freshwater stickleback from structurally complex habitats have retained the ancestral telencephalon morphology, but populations that feed more in open habitats on plankton have evolved larger, laterally concave telencephala.     

Selection on the morphology–physiology‐performance nexus: Lessons from freshwater stickleback morphs

Conspecifics inhabiting divergent environments frequently differ in morphology, physiology, and performance, but the interrelationships amongst traits and with Darwinian fitness remains poorly understood. We investigated population differentiation in morphology, metabolic rate, and swimming performance in three‐spined sticklebacks (Gasterosteus aculeatus L.), contrasting a marine/ancestral population with two distinct freshwater morphotypes derived from it: the “typical” low‐plated morph, and a unique “small‐plated” morph. We test the hypothesis that similar to plate loss in other freshwater populations, reduction in lateral plate size also evolved in response to selection. Additionally, we test how morphology, physiology, and performance have evolved in concert as a response to differences in selection between marine and freshwater environments. We raised pure‐bred second‐generation fish originating from three populations and quantified their lateral plate coverage, burst‐ and critical swimming speeds, as well as standard and active metabolic rates. Using a multivariate Q_ST‐F_ST framework, we detected signals of directional selection on metabolic physiology and lateral plate coverage, notably demonstrating that selection is responsible for the reduction in lateral plate coverage in a small‐plated stickleback population. We also uncovered signals of multivariate selection amongst all bivariate trait combinations except the two metrics of swimming performance. Divergence between the freshwater and marine populations exceeded neutral expectation in morphology and in most physiological and performance traits, indicating that adaptation to freshwater habitats has occurred, but through different combinations of traits in different populations. These results highlight both the complex interplay between morphology, physiology and performance in local adaptation, and a framework for their investigation.     

A gene with major phenotypic effects as a target for selection vs. homogenizing gene flow

Genes with major phenotypic effects facilitate quantifying the contribution of genetic vs. plastic effects to adaptive divergence. A classical example is Ectodysplasin (Eda), the major gene controlling lateral plate phenotype in three‐spined stickleback. Completely plated marine stickleback populations evolved repeatedly towards low‐plated freshwater populations, representing a prime example of parallel evolution by natural selection. However, many populations remain polymorphic for lateral plate number. Possible explanations for this polymorphism include relaxation of selection, disruptive selection or a balance between divergent selection and gene flow. We investigated 15 polymorphic stickleback populations from brackish and freshwater habitats in coastal North‐western Europe. At each site, we tracked changes in allele frequency at the Eda gene between subadults in fall, adults in spring and juveniles in summer. Eda genotypes were also compared for body size and reproductive investment. We observed a fitness advantage for the Eda allele for the low morph in freshwater and for the allele for the complete morph in brackish water. Despite these results, the differentiation at the Eda gene was poorly correlated with habitat characteristics. Neutral population structure was the best predictor of spatial variation in lateral plate number, suggestive of a substantial effect of gene flow. A meta‐analysis revealed that the signature of selection at Eda was weak compared to similar studies in stickleback. We conclude that a balance between divergent selection and gene flow can maintain stickleback populations polymorphic for lateral plate number and that ecologically relevant genes may not always contribute much to local adaptation, even when targeted by selection.     

Isolation between sympatric anadromous and resident threespine stickleback species in Mud Lake,Alaska

A sympatric pair of anadromous and resident freshwater threespine stickleback species (Gasterosteus aculeatus species complex) occurs in Mud Lake in the Matanuska-Susitna Valley, Alaska. The two forms differ in an array of morphological traits, including traits associated with predator defense (e.g., spine lengths) and trophic ecology (e.g., number of gill rakers). Mud Lake is only the third lake reported to have anadromous stickleback (which have a complete row of lateral plates) coexisting with low-plated resident stickleback in the absence of intermediate partially plated fish. Microhabitat and seasonal isolation appear to contribute to reproductive isolation between the two forms.     

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.     

The genetics of adaptive shape shift in stickleback: pleiotropy and effect size

The distribution of effect sizes of genes underlying adaptation is unknown ( Orr 2005 ). Are suites of traits that diverged under natural selection controlled by a few pleiotropic genes of large effect (major genes model), by many independently acting genes of small effect (infinitesimal model), or by a combination, with frequency inversely related to effect size (geometric model)? To address this we carried out a quantitative trait loci (QTL) study of a suite of 54 position traits describing body shapes of two threespine stickleback species: an ancestral Pacific marine form and a highly derived benthic species inhabiting a geologically young lake. About half of the 26 detected QTL affected just one coordinate and had small net effects, but several genomic regions affected multiple aspects of shape and had large net effects. The distribution of effect sizes followed the gamma distribution, as predicted by the geometric model of adaptation when detection limits are taken into account. The sex‐determining chromosome region had the largest effect of any QTL. Ancestral sexual dimorphism was similar to the direction of divergence, and was largely eliminated during freshwater adaptation, suggesting that sex differences may provide variation upon which selection can act. Several shape QTL are linked to Eda, a major gene responsible for reduction of lateral body armor in freshwater. Our results are consistent with predictions of the geometric model of adaptation. Shape evolution in stickleback results from a few genes with large and possibly widespread effects and multiple genes of smaller effect.     

Predation-imposed selection on threespine stickleback (Gasterosteus aculeatus) morphology: a test of the refuge use hypothesis

The transition from marine to freshwater life in the threespine stickleback (Gasterosteus aculeatus) is accompanied by complex morphological changes-including reduction in bony armor and change in body shape-but experimental evidence for the selective agents behind these evolutionary transitions is sparse. We investigated whether selection by predatory fish affects threespine stickleback morphology differentially when refuge is absent (pelagic lifestyle-ancestral condition) or present (benthic lifestyle-derived condition). Our results show that selection favors low numbers of lateral plates in habitats with refuge, whereas fully plated individuals have a selective advantage in habitats without refuge. We also found that a decrease in the length of the caudal peduncle increased survival probability, irrespective of habitat. The effect of spine lengths on survival was evident only in a multivariate analysis of selection, implying that it is essential to account for phenotypic and genetic correlations between traits before drawing conclusions about the effects of selection on single traits. Apart from uncovering targets and patterns of predator-induced selection on threespine stickleback morphology, our results provide direct evidence to support the hypothesis that differences in antipredator strategies in pelagic versus benthic sticklebacks could play a role in the repeated, independent cases of plate number reduction following freshwater colonization in this species.     

Growth pattern and survival in populations of Poecilia vivipara (Teleostei; Poeciliidae) inhabiting an environmental gradient: a common garden study

We performed a common garden experiment to assess the existence of genetic differences on growth and body size between two populations of Poecilia vivipara inhabiting extremes of an environmental gradient caused by water salinity in lagoons of Northern Rio de Janeiro State, Brazil: the Campelo lagoon (freshwater) and Açu lagoon (brackish/saltwater). The two populations show extreme differences in average phenotypes for body size, shape and life history (freshwater populations with smaller body size, lower fecundity and larger reproductive allotment). Pregnant females were brought to the lab and the offspring from both groups were kept in a common recirculating system with freshwater. Standard length and survival were measured weekly over a period of 200 days and growth models were fitted and selected with information criteria. The offspring originally from the brackish water lagoon presented larger asymptotic length, higher maximum growth rate but lower survival than the offspring originally from the freshwater lagoon. Potential confounding variables such as density differences due to mortality and maternal effects (offspring size) were included as covariates in comparisons of growth rates between groups. The results are consistent with phenotypic differences among populations having some genetic basis, and with the existence of a trade-off between growth and maintenance due to the high growth/low survival observed in the group that changed from salt to freshwater. Comparisons of captive and natural populations suggest that the influence of environmental factors, such as salinity, food availability, fish density and predation should also be considered relevant to explain phenotypic variation in this system.     

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

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

Intraguild predation leads to genetically based character shifts in the threespine stickleback

Intraguild predation is a common ecological interaction that occurs when a species preys upon another species with which it competes. The interaction is potentially a mechanism of divergence between intraguild prey (IG‐prey) populations, but it is unknown if cases of character shifts in IG‐prey are an environmental or evolutionary response. We investigated the genetic basis and inducibility of character shifts in threespine stickleback from lakes with and without prickly sculpin, a benthic intraguild predator (IG‐predator). Wild populations of stickleback sympatric with sculpin repeatedly show greater defensive armor and water column height preference. We laboratory‐raised stickleback from lakes with and without sculpin, as well as marine stickleback, and found that differences between populations in armor, body shape, and behavior persisted in a common garden. Within the common garden, we raised stickleback half‐families from multiple populations in the presence and absence of sculpin. Although the presence of sculpin induced trait changes in the marine stickleback, we did not observe an induced response in the freshwater stickleback. Behavioral and morphological trait differences between freshwater populations thus have a genetic basis and suggest an evolutionary response to intraguild predation.     

Lateral plate number in low‐plated threespine stickleback: a study of plasticity and heritability

In the threespine stickleback Gasterosteus aculeatus model system, phenotypes are often classified into three morphs according to lateral plate number. Morph identity has been shown to be largely genetically determined, but substantial within‐morph variation in plate number exists. In this study, we test whether plate number has a plastic component in response to salinity in the low‐plated morph using a split‐clutch experiment where families were split in two, one half raised in water at 0 and the other at 30 ppt salt. We find a small salinity‐induced plastic effect on plate number in an unexpected direction, opposite to what we predicted: Fish raised in freshwater on average have slightly more plates than fish raised in saltwater. Our results confirm that heritability of plate number is high. Additionally, we find that variance in plate number at the family level can be predicted from other family level traits, which might indicate that epistatic interactions play a role in creating the observed pattern of lateral plate number variation.     

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.     

Predation's role in repeated phenotypic and genetic divergence of armor in threespine stickleback

Predator-driven divergent selection may cause differentiation in defensive armor in threespine stickleback: (1) predatory fish and birds favor robust armor, whereas (2) predaceous aquatic insects favor armor reduction. Although (1) is well established, no direct experimental evidence exists for (2). I examined the phenotypic and genetic consequences of insect predation using F₂ families from crosses between freshwater and marine stickleback populations. I measured selection on body size, and size-adjusted spine (dorsal and pelvic) and pelvic girdle length, by splitting juvenile F₂ families between control and insect predation treatments, set in pond enclosures. I also examined the effect of insect predation on Ectodysplasin ( Eda ), a gene physically linked to quantitative trait loci for lateral plate number, spine length, and body shape. Insect predation resulted in: (1) significant selection for larger juvenile size, and shorter dorsal spine and pelvic girdle length, (2) higher mortality of individuals missing the pelvic girdle, and (3) selection in favor of the low armor Eda allele. Predatory insects favor less stickleback armor, likely contributing to the widespread reduction of armor in freshwater populations. Because size strongly influences mate choice, predator-driven divergent selection on size may play a substantial role in byproduct reproductive isolation and speciation in threespine stickleback.