Rapid ecological speciation in three-spined stickleback Gasterosteus aculeatus from Middleton Island, Alaska: the roles of selection and geographic isolation

Dimensionality and Modularity of Adaptive Variation: Divergence in Threespine Stickleback from Diverse Environments

AbstractPopulations are subjected to diverse environmental conditions that affect fitness and induce evolutionary or plastic responses, resulting in phenotypic divergence. Some authors contend that such divergence is concentrated along a single major axis of trait covariance even if that axis does not lead populations directly toward a fitness optimum. Other authors argue that divergence can occur readily along many phenotype axes at the same time. We use populations of threespine stickleback (Gasterosteus aculeatus) from 14 lakes with contrasting ecological conditions to find some resolution along the continuum between these two extremes. Unlike many previous studies, we included several functional suites of traits (defensive, swimming, trophic) potentially subject to different sources of selection. We find that populations exhibit dimensionality of divergence that is high enough to preclude a history of constraint along a single axis-both for divergence in multivariate mean trait values and for the structure of trait covariances. Dimensionality varied among trait suites and were strongly influenced by the inclusion of specific traits, and integration of trait suites varied between populations. We leverage this variation into new insights about the process of divergence and suggest that similar analyses could increase understanding of other adaptive radiations.

Freshwater Colonization, Adaptation, and Genomic Divergence in Threespine Stickleback

The Threespine Stickleback is ancestrally a marine fish, but many marine populations breed in fresh water (i.e., are anadromous), facilitating their colonization of isolated freshwater habitats a few years after they form. Repeated adaptation to fresh water during at least 10 My and continuing today has led to Threespine Stickleback becoming a premier system to study rapid adaptation. Anadromous and freshwater stickleback breed in sympatry and may hybridize, resulting in introgression of freshwater-adaptive alleles into anadromous populations, where they are maintained at low frequencies as ancient standing genetic variation. Anadromous stickleback have accumulated hundreds of freshwater-adaptive alleles that are disbursed as few loci per marine individual and provide the basis for adaptation when they colonize fresh water. Recent whole-lake experiments in lakes around Cook Inlet, Alaska have revealed how astonishingly rapid and repeatable this process is, with the frequency of 40% of the identified freshwater-adaptive alleles increasing from negligible (∼1%) in the marine founder to ≥50% within ten generations in fresh water, and freshwater phenotypes evolving accordingly. These high rates of genomic and phenotypic evolution imply very intense directional selection on phenotypes of heterozygotes. Sexual recombination rapidly assembles freshwater-adaptive alleles that originated in different founders into multilocus freshwater haplotypes, and regions important for adaptation to freshwater have suppressed recombination that keeps advantageous alleles linked within large haploblocks. These large haploblocks are also older and appear to have accumulated linked advantageous mutations. The contemporary evolution of Threespine Stickleback has provided broadly applicable insights into the mechanisms that facilitate rapid adaptation.

Phosphorus limitation does not drive loss of bony lateral plates in freshwater stickleback (Gasterosteus aculeatus)

Connecting the selective forces that drive the evolution of phenotypes to their underlying genotypes is key to understanding adaptation, but such connections are rarely tested experimentally. Threespine stickleback (Gasterosteus aculeatus) are a powerful model for such tests because genotypes that underlie putatively adaptive traits have been identified. For example, a regulatory mutation in the Ectodysplasin (Eda) gene causes a reduction in the number of bony armor plates, which occurs rapidly and repeatedly when marine sticklebacks invade freshwater. However, the source of selection on plate loss in freshwater is unknown. Here, we tested whether dietary reduction of phosphorus can account for selection on plate loss due to a growth advantage of low-plated fish in freshwater. We crossed marine fish heterozygous for the 16 kilobase freshwater Eda haplotype and compared the growth of offspring with different genotypes under contrasting levels of dietary phosphorus in both saltwater and freshwater. Eda genotype was not associated with growth differences in any treatment, or with mechanisms that could mitigate the impacts of phosphorus limitation, such as differential phosphorus deposition, phosphorus excretion, or intestine length. This study highlights the importance of experimentally testing the putative selective forces acting on phenotypes and their underlying genotypes in the wild.

Adaptation via pleiotropy and linkage: Association mapping reveals a complex genetic architecture within the stickleback Eda locus

Genomic mapping of the loci associated with phenotypic evolution has revealed genomic "hotspots," or regions of the genome that control multiple phenotypic traits. This clustering of loci has important implications for the speed and maintenance of adaptation and could be due to pleiotropic effects of a single mutation or tight genetic linkage of multiple causative mutations affecting different traits. The threespine stickleback (Gasterosteus aculeatus) is a powerful model for the study of adaptive evolution because the marine ecotype has repeatedly adapted to freshwater environments across the northern hemisphere in the last 12,000 years. Freshwater ecotypes have repeatedly fixed a 16 kilobase haplotype on chromosome IV that contains Ectodysplasin (Eda), a gene known to affect multiple traits, including defensive armor plates, lateral line sensory hair cells, and schooling behavior. Many additional traits have previously been mapped to a larger region of chromosome IV that encompasses the Eda freshwater haplotype. To identify which of these traits specifically map to this adaptive haplotype, we made crosses of rare marine fish heterozygous for the freshwater haplotype in an otherwise marine genetic background. Further, we performed fine-scale association mapping in a fully interbreeding, polymorphic population of freshwater stickleback to disentangle the effects of pleiotropy and linkage on the phenotypes affected by this haplotype. Although we find evidence that linked mutations have small effects on a few phenotypes, a small 1.4-kb region within the first intron of Eda has large effects on three phenotypic traits: lateral plate count, and both the number and patterning of the posterior lateral line neuromasts. Thus, the Eda haplotype is a hotspot of adaptation in stickleback due to both a small, pleiotropic region affecting multiple traits as well as multiple linked mutations affecting additional traits.

Genetic divergence outpaces phenotypic evolution among threespine stickleback populations in old freshwater habitats

Species such as threespine stickleback (Gasterosteus aculeatus) that inhabit divergent selective environments and that have diversified on different time scales can be of value for understanding evolutionary processes. Here we synthesize high-resolution genotypic and phenotypic data to explore a largely unstudied distribution of threespine stickleback populations living in oceanic and freshwater habitats along coastal and inland regions of Oregon. Many inland aquatic habitats of Oregon remained unglaciated during the last ice age, meaning that some extant Oregon lake and river stickleback may have descended from freshwater populations established long before more well-studied, post-glacial freshwater populations. To address the degree of congruence between genetic and phenotypic divergence, we directly compared Oregon stickleback to much younger (post-glacial) Alaskan populations. We found phenotypic variation in Oregon stickleback to be primarily partitioned between oceanic and freshwater habitats, as has been documented in other stickleback systems. However, the main axis of genetic divergence was between coastal and inland regions regardless of habitat type. Furthermore, when comparing patterns between Oregon and Alaska we found similar levels of phenotypic divergence, but much greater genetic divergence among Oregon's populations. The Oregon stickleback system therefore appears well suited for future studies linking genotypic and phenotypic change, further extending the utility of this small fish to provide general insights into evolutionary processes.

Evidence for contemporary and historical gene flow between guppy populations in different watersheds, with a test for associations with adaptive traits

In dendritic river systems, gene flow is expected to occur primarily within watersheds. Yet, rare cross-watershed transfers can also occur, whether mediated by (often historical) geological events or (often contemporary) human activities. We explored these events and their potential evolutionary consequences by analyzing patterns of neutral genetic variation (microsatellites) and adaptive phenotypic variation (male color) in wild guppies (Poecilia reticulata) distributed across two watersheds in northern Trinidad. We found the expected signatures of within-watershed gene flow; yet we also inferred at least two instances of cross-watershed gene flow-one in the upstream reaches and one further downstream. The upstream cross-watershed event appears to be very recent (41 ± 13 years), suggesting dispersal via recent flooding or undocumented human-mediated transport. The downstream cross-watershed event appears to be considerably older (577 ± 265 years), suggesting a role for rare geological or climatological events. Alongside these strong signatures of both contemporary and historical gene flow, we found little evidence of impacts on presumably adaptive phenotypic differentiation, except perhaps in the one instance of very recent cross-watershed gene flow. Selection in this system seems to overpower gene flow-at least on the spatiotemporal scales investigated here.

Repeated Selection of Alternatively Adapted Haplotypes Creates Sweeping Genomic Remodeling in Stickleback

Heterogeneous genetic divergence can accumulate across the genome when populations adapt to different habitats while still exchanging alleles. How long does diversification take and how much of the genome is affected? When divergence occurs in parallel from standing genetic variation, how often are the same haplotypes involved? We explore these questions using restriction site-associated DNA sequencing genotyping data and show that broad-scale genomic repatterning, fueled by copious standing variation, can emerge in just dozens of generations in replicate natural populations of threespine stickleback fish (Gasterosteus aculeatus). After the catastrophic 1964 Alaskan earthquake, marine stickleback colonized newly created ponds on seismically uplifted islands. We find that freshwater fish in these young ponds differ from their marine ancestors across the same genomic segments previously shown to have diverged in much older lake populations. Outside of these core divergent regions the genome shows no population structure across the ocean-freshwater divide, consistent with strong local selection acting in alternative environments on stickleback populations still connected by significant gene flow. Reinforcing this inference, a majority of divergent haplotypes that are at high frequency in ponds are detectable in the sea, even across great geographic distances. Building upon previous population genomics work in this model species, our data suggest that a long history of divergent selection and gene flow among stickleback populations in oceanic and freshwater habitats has maintained polymorphisms of alternatively adapted DNA sequences that facilitate parallel evolution.

Impact of the huge 2011 Tohoku-oki tsunami on the phenotypes and genotypes of Japanese coastal threespine stickleback populations

On March 11, 2011, a large earthquake occurred, causing a tsunami which struck the Pacific coast of northeast Japan. We investigated the ecological and genetic effects of the large tsunami on the threespine stickleback (genus Gasterosteus) populations in Otsuchi Town, which was one of the most severely damaged areas after the tsunami. Our environmental surveys showed that spring water may have contributed to the habitat recovery. Morphological analysis of the stickleback before and after the tsunami showed morphological shifts in the gill raker number, which is a foraging trait. Genetic analyses revealed that the allelic richness of one population was maintained after the tsunami, whereas that of another decreased in 2012 and then started to recover in 2013. Additionally, we found that the large tsunami and ground subsidence created new spring water-fed pools with sticklebacks, suggesting that the tsunami brought sticklebacks into these pools. Genetic analysis of this population showed that this population might be derived from hybridization between freshwater Gasterosteus aculeatus and anadromous G. nipponicus. Overall, our data indicate that tsunamis can influence morphologies and genetic structures of freshwater fishes. Furthermore, spring water may play important roles in the maintenance and creation of fish habitats, faced with environmental disturbance.

Evolution of stickleback in 50 years on earthquake-uplifted islands

How rapidly can animal populations in the wild evolve when faced with sudden environmental shifts? Uplift during the 1964 Great Alaska Earthquake abruptly created freshwater ponds on multiple islands in Prince William Sound and the Gulf of Alaska. In the short time since the earthquake, the phenotypes of resident freshwater threespine stickleback fish on at least three of these islands have changed dramatically from their oceanic ancestors. To test the hypothesis that these freshwater populations were derived from oceanic ancestors only 50 y ago, we generated over 130,000 single-nucleotide polymorphism genotypes from more than 1,000 individuals using restriction site-associated DNA sequencing (RAD-seq). Population genomic analyses of these data support the hypothesis of recent and repeated, independent colonization of freshwater habitats by oceanic ancestors. We find evidence of recurrent gene flow between oceanic and freshwater ecotypes where they co-occur. Our data implicate natural selection in phenotypic diversification and support the hypothesis that the metapopulation organization of this species helps maintain a large pool of genetic variation that can be redeployed rapidly when oceanic stickleback colonize freshwater environments. We find that the freshwater populations, despite population genetic analyses clearly supporting their young age, have diverged phenotypically from oceanic ancestors to nearly the same extent as populations that were likely founded thousands of years ago. Our results support the intriguing hypothesis that most stickleback evolution in fresh water occurs within the first few decades after invasion of a novel environment.

Phenotypic plasticity and epigenetic marking: an assessment of evidence for genetic accommodation

The relationship between genotype (which is inherited) and phenotype (the target of selection) is mediated by environmental inputs on gene expression, trait development, and phenotypic integration. Phenotypic plasticity or epigenetic modification might influence evolution in two general ways: (1) by stimulating evolutionary responses to environmental change via population persistence or by revealing cryptic genetic variation to selection, and (2) through the process of genetic accommodation, whereby natural selection acts to improve the form, regulation, and phenotypic integration of novel phenotypic variants. We provide an overview of models and mechanisms for how such evolutionary influences may be manifested both for plasticity and epigenetic marking. We point to promising avenues of research, identifying systems that can best be used to address the role of plasticity in evolution, as well as the need to apply our expanding knowledge of genetic and epigenetic mechanisms to our understanding of how genetic accommodation occurs in nature. Our review of a wide variety of studies finds widespread evidence for evolution by genetic accommodation.

Genetic divergence of a sympatric lake-resident-anadromous three-spined stickleback Gasterosteus aculeatus species pair

The genetic relationship between sympatric, morphologically divergent populations of anadromous and lake-resident three-spined stickleback Gasterosteus aculeatus in the Jim Creek drainage of Cook Inlet, Alaska, was examined using microsatellite loci and mitochondrial d-loop sequence data. Resident samples differed substantially from sympatric anadromous samples in the Jim Creek drainage with the magnitude of the genetic divergence being similar to that between allopatric resident and anadromous populations in other areas. Resident samples were genetically similar within the Jim Creek drainage, as were the anadromous samples surveyed. Neighbour-joining and Structure cluster analysis grouped the samples into four genetic clusters by ecomorph (anadromous v. all resident) and geographic location of the resident samples (Jim Creek, Mat-Su and Kenai). There was no evidence of hybridization between resident and anadromous G. aculeatus in the Jim Creek drainage, which thus appear to be reproductively isolated.

The population structure and recent colonization history of Oregon threespine stickleback determined using restriction-site associated DNA-sequencing

Understanding how genetic variation is partitioned across genomes within and among populations is a fundamental problem in ecological and evolutionary genetics. To address this problem, we studied the threespine stickleback fish, which has repeatedly undergone parallel phenotypic and genetic differentiation when oceanic fish have invaded freshwater habitats. While significant evolutionary genetic research has been performed using stickleback from geographic regions that have been deglaciated in the last 20 000 years, less research has focused on freshwater populations that predate the last glacial maximum. We performed restriction-site associated DNA-sequencing (RAD-seq) based population genomic analyses on stickleback from across Oregon, which was not glaciated during the last maximum. We sampled stickleback from coastal, Willamette Basin and central Oregon sites, analysed their genetic diversity using RAD-seq, performed structure analyses, reconstructed their phylogeographic history and tested the hypothesis of recent stickleback introduction into central Oregon, where incidence of this species was only recently documented. Our results showed a clear phylogeographic break between coastal and inland populations, with oceanic populations exhibiting the lowest levels of divergence from one another. Willamette Basin and central Oregon populations formed a clade of closely related populations, a finding consistent with a recent introduction of stickleback into central Oregon. Finally, genome-wide analysis of genetic diversity (π) and correlations of alleles within individuals in subpopulations (FIS) supported a role for introgressive hybridization in coastal populations and a recent expansion in central Oregon. Our results exhibit the power of next-generation sequencing genomic approaches such as RAD-seq to identify both historical population structure and recent colonization history.

Repeated lake-stream divergence in stickleback life history within a Central European lake basin

Life history divergence between populations inhabiting ecologically distinct habitats might be a potent source of reproductive isolation, but has received little attention in the context of speciation. We here test for life history divergence between threespine stickleback inhabiting Lake Constance (Central Europe) and multiple tributary streams. Otolith analysis shows that lake fish generally reproduce at two years of age, while their conspecifics in all streams have shifted to a primarily annual life cycle. This divergence is paralleled by a striking and consistent reduction in body size and fecundity in stream fish relative to lake fish. Stomach content analysis suggests that life history divergence might reflect a genetic or plastic response to pelagic versus benthic foraging modes in the lake and the streams. Microsatellite and mitochondrial markers further reveal that life history shifts in the different streams have occurred independently following the colonization by Lake Constance stickleback, and indicate the presence of strong barriers to gene flow across at least some of the lake-stream habitat transitions. Given that body size is known to strongly influence stickleback mating behavior, these barriers might well be related to life history divergence.

Population genetic dynamics of three-spined sticklebacks (Gasterosteus aculeatus) in anthropogenic altered habitats

In industrialized and/or agriculturally used landscapes, inhabiting species are exposed to a variety of anthropogenic changes in their environments. Genetic diversity may be reduced if populations encounter founder events, bottlenecks, or isolation. Conversely, genetic diversity may increase if populations adapt to changes in selective regimes in newly created habitats. With the present study, genetic variability of 918 sticklebacks from 43 samplings (21.3 ± 3.8 per sample) at 36 locations from cultivated landscapes in Northwest Germany was analyzed at nine neutral microsatellite loci. To test if differentiation is influenced by habitat alterations, sticklebacks were collected from ancient running waters and adjacent artificial stagnant waters, from brooks with salt water inflow of anthropogenic and natural origin and adjacent freshwater sites. Overall population structure was dominated by isolation by distance (IBD), which was significant across all populations, and analysis of molecular variance (AMOVA) revealed that 10.6% of the variation was explained by river catchment area. Populations in anthropogenic modified habitats deviated from the general IBD structure and in the AMOVA, grouping by habitat type running/stagnant water explained 4.9% of variation and 1.4% of the variation was explained by salt-/freshwater habitat. Sticklebacks in salt-polluted water systems seem to exhibit elevated migratory activity between fresh- and saltwater habitats, reducing IBD. In other situations, populations showed distinct signs of genetic isolation, which in some locations was attributed to mechanical migration barriers, but in others to potential anthropogenic induced bottleneck or founder effects. The present study shows that anthropogenic habitat alterations may have diverse effects on the population genetic structure of inhabiting species. Depending on the type of habitat change, increased genetic differentiation, diversification, or isolation are possible consequences.

Partial reproductive isolation of a recently derived resident-freshwater population of threespine stickleback (Gasterosteus aculeatus) from its putative anadromous ancestor

We used no-choice mating trials to test for assortative mating between a newly derived resident-freshwater population (8-22 generations since founding) of threespine stickleback (Gasterosteus aculeatus) in Loberg Lake, Alaska and its putative anadromous ancestor as well as a morphologically convergent but distantly related resident-freshwater population. Partial reproductive isolation has evolved between the Loberg Lake population and its ancestor within a remarkably short time period. However, Loberg stickleback readily mate with morphologically similar, but distantly related resident-freshwater stickleback. Partial premating isolation is asymmetrical; anadromous females and smaller resident-freshwater males from Loberg Lake readily mate, but the anadromous males and smaller Loberg females do not. Our results indicate that premating isolation can begin to evolve in allopatry within a few generations after isolation as a correlated effect of evolution of reduced body size.

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.

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.