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

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.

Perchlorate exposure does not induce obesity or non-alcoholic fatty liver disease in zebrafish

Perchlorate is a water-soluble contaminant found throughout the United States and many other countries. Perchlorate competitively inhibits iodide uptake at the sodium/iodide symporter, reducing thyroid hormone synthesis, which can lead to hypothyroidism and metabolic syndromes. Chronic perchlorate exposure induces hepatic steatosis and non-alcoholic fatty liver disease (NAFLD) in developing threespine stickleback (Gasterosteus aculeatus). We hypothesized that perchlorate would also induce zebrafish (Danio rerio) to develop phenotypes consistent with NAFLD and to accumulate lipids throughout the body. We exposed zebrafish embryos to four concentrations of perchlorate treated water (10μg/L, 10mg/L, 30mg/L, and 100mg/L) and a control (0mg/L) over the course of 133 days. Adult zebrafish were euthanized, sectioned, H&E and Oil Red-O stained, and analyzed for liver morphology and whole body lipid accumulation. In a representative section of the liver, we counted the number of lipid droplets and measured the area of each droplet and the total lipid area. For whole body analysis, we calculated the ratio of lipid area to body area within a section. We found that zebrafish exposed to perchlorate did not differ in any measured liver variables or whole body lipid area when compared to controls. In comparison to stickleback, we see a trend that control stickleback accumulate more lipids in their liver than do control zebrafish. Differences between the species indicate that obesogenic effects due to perchlorate exposure are not uniform across fish species, and likely are mediated by evolutionary differences related to geographic location. For example, high latitude fishes such as stickleback evolved to deposit lipid stores for over-winter survival, which may lead to more pronounced obesogenic effects than seen in tropical fish such as zebrafish.

Courtship behavior, nesting microhabitat, and assortative mating in sympatric stickleback species pairs

The maintenance of reproductive isolation in the face of gene flow is a particularly contentious topic, but differences in reproductive behavior may provide the key to explaining this phenomenon. However, we do not yet fully understand how behavior contributes to maintaining species boundaries. How important are behavioral differences during reproduction? To what extent does assortative mating maintain reproductive isolation in recently diverged populations and how important are "magic traits"? Assortative mating can arise as a by-product of accumulated differences between divergent populations as well as an adaptive response to contact between those populations, but this is often overlooked. Here we address these questions using recently described species pairs of three-spined stickleback (Gasterosteus aculeatus), from two separate locations and a phenotypically intermediate allopatric population on the island of North Uist, Scottish Western Isles. We identified stark differences in the preferred nesting substrate and courtship behavior of species pair males. We showed that all males selectively court females of their own ecotype and all females prefer males of the same ecotype, regardless of whether they are from species pairs or allopatric populations. We also showed that mate choice does not appear to be driven by body size differences (a potential "magic trait"). By explicitly comparing the strength of these mating preferences between species pairs and single-ecotype locations, we were able to show that present levels of assortative mating due to direct mate choice are likely a by-product of other adaptations between ecotypes, and not subject to obvious selection in species pairs. Our results suggest that ecological divergence in mating characteristics, particularly nesting microhabitat may be more important than direct mate choice in maintaining reproductive isolation in stickleback species pairs.

Imprecise naming: the anadromous and the sea spawning threespine stickleback should be discriminated by names

Two ecological forms of the threespine stickleback Gasterosteus aculeatus – a strictly marine form and an anadromous form – are often merged in the literature as a single “marine” form. Because we know virtually nothing of the life style of the two oceanic ecotypes in the sea and consequently nothing on reproductive isolation and gene flow I argue for a precise use of the ecological terms “marine” and “anadromous” for these two ecotypes. These terms should be self-describing. The frequent use of terms incorrectly describing intraspecific variation and life style of ecotypes can bias studies on community composition and interactions of populations.

Sexual isolation promotes divergence between parapatric lake and stream stickleback

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

Evolution of reproductive isolation in stickleback fish

To understand how new species form and what causes their collapse, we examined how reproductive isolation evolves during the speciation process, considering species pairs with little to extensive divergence, including a recently collapsed pair. We estimated many reproductive barriers in each of five sets of stickleback fish species pairs using our own data and decades of previous work. We found that the types of barriers important early in the speciation process differ from those important late. Two premating barriers-habitat and sexual isolation-evolve early in divergence and remain two of the strongest barriers throughout speciation. Premating isolation evolves before postmating isolation, and extrinsic isolation is far stronger than intrinsic. Completing speciation, however, may require postmating intrinsic incompatibilities. Reverse speciation in one species pair was characterized by significant loss of sexual isolation. We present estimates of barrier strengths before and after collapse of a species pair; such detail regarding the loss of isolation has never before been documented. Additionally, despite significant asymmetries in individual barriers, which can limit speciation, total isolation was essentially symmetric between species. Our study provides important insight into the order of barrier evolution and the relative importance of isolating barriers during speciation and tests fundamental predictions of ecological speciation.

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.

Fast evolution from precast bricks: genomics of young freshwater populations of threespine stickleback Gasterosteus aculeatus

Adaptation is driven by natural selection; however, many adaptations are caused by weak selection acting over large timescales, complicating its study. Therefore, it is rarely possible to study selection comprehensively in natural environments. The threespine stickleback (Gasterosteus aculeatus) is a well-studied model organism with a short generation time, small genome size, and many genetic and genomic tools available. Within this originally marine species, populations have recurrently adapted to freshwater all over its range. This evolution involved extensive parallelism: pre-existing alleles that adapt sticklebacks to freshwater habitats, but are also present at low frequencies in marine populations, have been recruited repeatedly. While a number of genomic regions responsible for this adaptation have been identified, the details of selection remain poorly understood. Using whole-genome resequencing, we compare pooled genomic samples from marine and freshwater populations of the White Sea basin, and identify 19 short genomic regions that are highly divergent between them, including three known inversions. 17 of these regions overlap protein-coding genes, including a number of genes with predicted functions that are relevant for adaptation to the freshwater environment. We then analyze four additional independently derived young freshwater populations of known ages, two natural and two artificially established, and use the observed shifts of allelic frequencies to estimate the strength of positive selection. Adaptation turns out to be quite rapid, indicating strong selection acting simultaneously at multiple regions of the genome, with selection coefficients of up to 0.27. High divergence between marine and freshwater genotypes, lack of reduction in polymorphism in regions responsible for adaptation, and high frequencies of freshwater alleles observed even in young freshwater populations are all consistent with rapid assembly of G. aculeatus freshwater genotypes from pre-existing genomic regions of adaptive variation, with strong selection that favors this assembly acting simultaneously at multiple loci.

Adaptation to novel environments is a keystone of evolution. There is only a handful of natural and experimental systems in which the process of adaptation has been studied in detail, and each studied system brings its own surprises with regard to the number of loci involved, dynamics of adaptation, extent of interactions between loci and of parallelism between different adapting populations. The threespine stickleback is an excellent model organism for evolutionary studies. Marine-derived freshwater populations of this species have consistently acquired a specific set of morphological, physiological and behavioral traits allowing them to reside in freshwater for their whole lifespan. Previous studies identified several genomic regions responsible for this adaptation. Here, using whole-genome sequencing, we compare the allele frequencies at such regions in four derived freshwater populations of known ages: two natural, and two artificially established in 1978. Knowledge of population ages allows us to infer the strength of selection that acted at these loci. Adaptation of threespine stickleback to freshwater is typically fast, and is driven by strong selection favoring pre-existing alleles that are likely present in the ancestral marine population at low frequencies; however, some of the adaptation may also be due to young population-specific alleles.

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.

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.

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.

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.