GLOBAL SURVEY OF MITOCHONDRIAL DNA SEQUENCES IN THE THREESPINE STICKLEBACK: EVIDENCE FOR RECENT MIGRATIONS

Environmental DNA unveils deep phylogeographic structure of a freshwater fish

Phylogeography bears an important part in ecology and evolution. However, current phylogeographic studies are largely constrained by limited numbers of individual samples. Using an environmental DNA (eDNA) assay for phylogeographic analyses, this study provides detailed information regarding the history of Siberian stone loach Barbatula toni, a primary freshwater fish across the whole range of Hokkaido, Japan. Based on an eDNA metabarcoding on 293 river water samples, we detected eDNA from B. toni in 189 rivers. A total of 51 samples, representing the entire island, were then selected from the B. toni eDNA-positive sample set for the subsequent analyses. To elucidate the phylogeographic structure of B. toni, newly developed eDNA metabarcoding primers (Barba-cytb-F/R) were applied to these samples, specifically targeting their haplotypic variation in cytochrome b. After a bioinformatic processing to mitigate haplotypic false positives, a total of 50 eDNA haplotypes were identified. Two regionally restricted, genetically distinct lineages of the species were revealed as a result of phylogeographic analyses on the haplotypes and tissue-derived DNA from B. toni. According to a molecular clock analysis, they have been genetically isolated for at least 1.5 million years, suggesting their ancient origin and colonisation of Hokkaido, presumably in the glacial periods. These results demonstrate how freshwater fishes can alter their distributions over evolutionary timescales and how eDNA assay can deepen our understanding of phylogeography.

Characterizing phenotypic diversity in marine populations of the threespine stickleback

The threespine stickleback (Gasterosteus aculeatus) is an important model for studying the evolution of vertebrate morphology. Sticklebacks inhabit freshwater, brackish, and marine northern hemisphere waters. Anadromous and marine populations (hereafter marine) are assumed to have remained unchanged morphologically from ancestral marine sticklebacks, despite marine environments varying on regional and local scales. Recent studies suggest that genetic and phenotypic structure exists in marine populations, yet the scale of this variation, and its ecological causes remain unclear. Our goal was to assess morphological trait variation in marine stickleback populations around Southern British Columbia (BC) and determine if oceanographic and habitat characteristics were associated with this variation. Between May-July 2019, we sampled 534 sticklebacks from 15 sites around Vancouver Island, a region characterized by a large diversity of oceanographic and habitat features. We characterized trait variation using two-dimensional (2D) geometric morphometric analysis, comparing individuals between oceanographic regions and habitats. We focused on head and body shape. We found that marine sticklebacks varied morphologically among and between regions and habitats, but the variation did not appear to be related to environmental variation. Sexual dimorphism was the largest source of variation, but oceanographic and habitat variables influenced differences between sexes. We concluded that marine sticklebacks offer abundant opportunities for expanding our knowledge of drivers of morphology.

Complex population history affects admixture analyses in nine-spined sticklebacks

Introgressive hybridization is an important process in evolution but challenging to identify, undermining the efforts to understand its role and significance. On the contrary, many analytical methods assume direct descent from a single common ancestor, and admixture among populations can violate their assumptions and lead to seriously biased results. A detailed analysis of 888 whole-genome sequences of nine-spined sticklebacks (Pungitius pungitius) revealed a complex pattern of population ancestry involving multiple waves of gene flow and introgression across northern Europe. The two recognized lineages were found to have drastically different histories, and their secondary contact zone was wider than anticipated, displaying a smooth gradient of foreign ancestry with some curious deviations from the expected pattern. Interestingly, the freshwater isolates provided peeks into the past and helped to understand the intermediate states of evolutionary processes. Our analyses and findings paint a detailed picture of the complex colonization history of northern Europe and provide backdrop against which introgression and its role in evolution can be investigated. However, they also expose the challenges in analyses of admixed populations and demonstrate how hidden admixture and colonization history misleads the estimation of admixture proportions and population split times.

Extreme intraspecific divergence in mitochondrial haplotypes makes the threespine stickleback fish an emerging evolutionary mutant model for mito-nuclear interactions

Mitochondrial DNA is primarily maternally inherited in most animals and evolves about 10 times faster than biparentally inherited nuclear DNA. Mitochondrial dysfunction (mt-dys) arises when interactions between the co-evolving mitochondrial and nuclear genomes are perturbed in essential processes like oxidative phosphorylation (OXPHOS). Over time mt-dys can lead to mitochondrial diseases (mt-diseases), which are surprisingly prevalent and include common diseases such as Alzheimer's, Parkinson's, and diabetes. Unfortunately, the strong impact that intraspecific mitochondrial and nuclear genetic variation has on mt-disease complicates its study and the development of effective treatments. Animal models have advanced our understanding of mt-disease but their relevance to human conditions is often limited by their relatively low nuclear genetic diversity. Many traditional laboratory models also typically have a single mitochondrial haplotype (mitotype), in stark contrast to over 5,000 mitotypes in humans worldwide. The threespine stickleback fish has an evolutionary history that has made it a favorable evolutionary mutant model (EMM) for studying mito-nuclear interactions and possibly mt-diseases. EMMs are species with naturally evolved states that mimic maladaptive human diseases. In threespine stickleback, a period of isolation followed by introgression of the mitochondrial genome from a sister species resulted in the maintenance of two distinct mitochondrial haplotypes which continue to segregate within many populations of wild stickleback. The existence of two mitogenomes segregating in numerous genetically diverse populations provides a unique system for exploring complex mito-nuclear dynamics. Here we provide the first complete coding region analysis of the two threespine stickleback mitotypes, whose mitogenomic divergence exceeds that of other mammalian models for mitochondrial disease and even that between ancient and modern humans. We find that divergence is not uniform across the mitogenome, but primarily impacts protein coding genes, and significantly impacts proteins in Complex I of OXPHOS. The full characterization of these highly divergent intraspecific mitotypes provides a foundation for the development of threespine stickleback as an EMM for mito-nuclear interactions.

Genomic changes underlying repeated niche shifts in an adaptive radiation

In adaptive radiations, single lineages rapidly diversify by adapting to many new niches. Little is known yet about the genomic mechanisms involved, that is, the source of genetic variation or genomic architecture facilitating or constraining adaptive radiation. Here, we investigate genomic changes associated with repeated invasion of many different freshwater niches by threespine stickleback in the Haida Gwaii archipelago, Canada, by resequencing single genomes from one marine and 28 freshwater populations. We find 89 likely targets of parallel selection in the genome that are enriched for old standing genetic variation. In contrast to theoretical expectations, their genomic architecture is highly dispersed with little clustering. Candidate genes and genotype-environment correlations match the three major environmental axes predation regime, light environment, and ecosystem size. In a niche space with these three dimensions, we find that the more divergent a new niche from the ancestral marine habitat, the more loci show signatures of parallel selection. Our findings suggest that the genomic architecture of parallel adaptation in adaptive radiation depends on the steepness of ecological gradients and the dimensionality of the niche space.

Threespine Stickleback in Lake Constance: The Ecology and Genomic Substrate of a Recent Invasion

Invasive species can be powerful models for studying contemporary evolution in natural environments. As invading organisms often encounter new habitats during colonization, they will experience novel selection pressures. Threespine stickleback (Gasterosteus aculeatus complex) have recently colonized large parts of Switzerland and are invasive in Lake Constance. Introduced to several watersheds roughly 150 years ago, they spread across the Swiss Plateau (400–800 m a.s.l.), bringing three divergent hitherto allopatric lineages into secondary contact. As stickleback have colonized a variety of different habitat types during this recent range expansion, the Swiss system is a useful model for studying contemporary evolution with and without secondary contact. For example, in the Lake Constance region there has been rapid phenotypic and genetic divergence between a lake population and some stream populations. There is considerable phenotypic variation within the lake population, with individuals foraging in and occupying littoral, offshore pelagic, and profundal waters, the latter of which is a very unusual habitat for stickleback. Furthermore, adults from the lake population can reach up to three times the size of adults from the surrounding stream populations, and are large by comparison to populations globally. Here, we review the historical origins of the threespine stickleback in Switzerland, and the ecomorphological variation and genomic basis of its invasion in Lake Constance. We also outline the potential ecological impacts of this invasion, and highlight the interest for contemporary evolution studies.

A key metabolic gene for recurrent freshwater colonization and radiation in fishes

Colonization of new ecological niches has triggered large adaptive radiations. Although some lineages have made use of such opportunities, not all do so. The factors causing this variation among lineages are largely unknown. Here, we show that deficiency in docosahexaenoic acid (DHA), an essential ω-3 fatty acid, can constrain freshwater colonization by marine fishes. Our genomic analyses revealed multiple independent duplications of the fatty acid desaturase gene Fads2 in stickleback lineages that subsequently colonized and radiated in freshwater habitats, but not in close relatives that failed to colonize. Transgenic manipulation of Fads2 in marine stickleback increased their ability to synthesize DHA and survive on DHA-deficient diets. Multiple freshwater ray-finned fishes also show a convergent increase in Fads2 copies, indicating its key role in freshwater colonization.

Specificity of resistance and geographic patterns of virulence in a vertebrate host-parasite system

BACKGROUND

Host genotype - parasite genotype co-evolutionary dynamics are influenced by local biotic and abiotic environmental conditions. This results in spatially heterogeneous selection among host populations. How such heterogeneous selection influences host resistance, parasite infectivity and virulence remains largely unknown. We hypothesized that different co-evolutionary trajectories of a vertebrate host-parasite association result in specific virulence patterns when assessed on a large geographic scale. We used two reference host populations of three-spined sticklebacks and nine strains of their specific cestode parasite Schistocephalus solidus from across the Northern Hemisphere for controlled infection experiments. Host and parasite effects on infection phenotypes including host immune gene expression were determined.

RESULTS

S. solidus strains grew generally larger in hosts coming from a population with high parasite diversity and low S. solidus prevalence (DE hosts). Hosts from a population with low parasite diversity and high S. solidus prevalence (NO hosts) were better able to control the parasite's growth, regardless of the origin of the parasite. Host condition and immunological parameters converged upon infection and parasite growth showed the same geographic pattern in both host types.

CONCLUSION

Our results suggest that NO sticklebacks evolved resistance against a variety of S. solidus strains, whereas DE sticklebacks are less resistant against S. solidus. Our data provide evidence that differences in parasite prevalence can cause immunological heterogeneity and that parasite size, a proxy for virulence and resistance, is, on a geographic scale, determined by main effects of the host and the parasite and less by an interaction of both genotypes.

Contemporary ancestor? Adaptive divergence from standing genetic variation in Pacific marine threespine stickleback

Background

Populations that have repeatedly colonized novel environments are useful for studying the role of ecology in adaptive divergence – particularly if some individuals persist in the ancestral habitat. Such “contemporary ancestors” can be used to demonstrate the effects of selection by comparing phenotypic and genetic divergence between the derived population and their extant ancestors. However, evolution and demography in these “contemporary ancestors” can complicate inferences about the source (standing genetic variation, de novo mutation) and pace of adaptive divergence. Marine threespine stickleback (Gasterosteus aculeatus) have colonized freshwater environments along the Pacific coast of North America, but have also persisted in the marine environment. To what extent are marine stickleback good proxies of the ancestral condition?

Results

We sequenced > 5800 variant loci in over 250 marine stickleback from eight locations extending from Alaska to California, and phenotyped them for platedness and body shape. Pairwise F_ST varied from 0.02 to 0.18. Stickleback were divided into five genetic clusters, with a single cluster comprising stickleback from Washington to Alaska. Plate number, Eda, body shape, and candidate loci showed evidence of being under selection in the marine environment. Comparisons to a freshwater population demonstrated that candidate loci for freshwater adaptation varied depending on the choice of marine populations.

Conclusions

Marine stickleback are structured into phenotypically and genetically distinct populations that have been evolving as freshwater stickleback evolved. This variation complicates their usefulness as proxies of the ancestors of freshwater populations. Lessons from stickleback may be applied to other “contemporary ancestor”-derived population studies.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1228-8) contains supplementary material, which is available to authorized users.

The genomic landscape at a late stage of stickleback speciation: High genomic divergence interspersed by small localized regions of introgression

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68-1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow.

Rapid adaptive evolution of colour vision in the threespine stickleback radiation

Vision is a sensory modality of fundamental importance for many animals, aiding in foraging, detection of predators and mate choice. Adaptation to local ambient light conditions is thought to be commonplace, and a match between spectral sensitivity and light spectrum is predicted. We use opsin gene expression to test for local adaptation and matching of spectral sensitivity in multiple independent lake populations of threespine stickleback populations derived since the last ice age from an ancestral marine form. We show that sensitivity across the visual spectrum is shifted repeatedly towards longer wavelengths in freshwater compared with the ancestral marine form. Laboratory rearing suggests that this shift is largely genetically based. Using a new metric, we found that the magnitude of shift in spectral sensitivity in each population corresponds strongly to the transition in the availability of different wavelengths of light between the marine and lake environments. We also found evidence of local adaptation by sympatric benthic and limnetic ecotypes to different light environments within lakes. Our findings indicate rapid parallel evolution of the visual system to altered light conditions. The changes have not, however, yielded a close matching of spectrum-wide sensitivity to wavelength availability, for reasons we discuss.

Fish Parasite Dinoflagellates Haidadinium ichthyophilum and Piscinoodinium Share a Recent Common Ancestor

The dinoflagellate Haidadinium ichthyophilum Buckland-Nicks, Reimchen and Garbary 1997 is an ectoparasite of the spine-deficient, three-spine stickleback Gasterosteus aculeatus L. Reimchen 1984, a fish endemic to Rouge Lake, Haida Gwaii. Haidadinium ichthyophilum proved difficult to assign taxonomically because its morphology and complex life cycle exhibited defining characteristics of both autotrophic and heterotrophic dinoflagellates, and was tentatively assigned to the Phytodiniales. Here, we characterized a 492 bp fragment of the small subunit ribosomal RNA (SSU rRNA) from preserved H. ichthyophilum cysts. In SSU phylogeny, H. ichthyophilum branches with the fish parasites, Piscinoodinium sp., strongly supporting the inclusion of H. ichthyophilum within the Suessiales.

Phylogeography and Conservation Genetics of the Ibero-Balearic Three-Spined Stickleback (Gasterosteus aculeatus)

Genetic isolation and drift may imperil peripheral populations of wide-ranging species more than central ones. Therefore, information about species genetic variability and population structure is invaluable for conservation managers. The Iberian populations of three-spined stickleback lie at the southwestern periphery of the European distribution of Gasterosteus aculeatus. This teleost is a protected species in Portugal and Spain and local extinctions have been reported in both countries during the last decades. Our objectives were (i) to determine whether the Iberian populations of G. aculeatus are unique or composed of any of the major evolutionary lineages previously identified and (ii) to assess the evolutionary potential of these peripheral populations. We genotyped 478 individuals from 17 sites at 10 polymorphic microsatellite loci to evaluate the genetic variability and differentiation of the Ibero-Balearic populations. We also sequenced 1,165 bp of the mitochondrial genome in 331 of those individuals in order to complement the estimates of genetic diversity in the Ibero-Balearic region. We predicted the evolutionary potential of the different sites analysed based on the contribution of each of them to total allelic/mitochondrial diversity. An intraspecific phylogeny at European level was reconstructed using our data from the mitochondrial cytochrome b gene (755 bp) and published sequences. The so-called Transatlantic, European and Mediterranean mitochondrial lineages were found to be present in the Ibero-Balearic region. Their phylogeography suggests a history of multiple colonisations. The nuclear results show, however, a strong correlation between population structure and drainage system. The following basins should be prioritised by conservation policies in order to preserve those populations with the highest evolutionary potential: the Portuguese Vouga and Tagus as well as the Spanish Majorca and Limia. Maintenance of their connectivity, control of exotic species and monitoring of habitat properties are strongly recommended in those areas. Genetic variation alone cannot, however, ensure the persistence of these peripheral southern populations of G. aculeatus. On the one hand, the analysis of a historical sample from Eastern Spain (Penyscola) revealed no genetic erosion, which suggests a fairly sudden extinction of that population. On the other hand, the reintroduction program implemented in the Valencian Community has mostly failed despite our finding of similar level of genetic diversity between the wild source and the captive-bred released individuals.

Region-wide and ecotype-specific differences in demographic histories of threespine stickleback populations, estimated from whole genome sequences

We analysed 81 whole genome sequences of threespine sticklebacks from Pacific North America, Greenland and Northern Europe, representing 16 populations. Principal component analysis of nuclear SNPs grouped populations according to geographical location, with Pacific populations being more divergent from each other relative to European and Greenlandic populations. Analysis of mitogenome sequences showed Northern European populations to represent a single phylogeographical lineage, whereas Greenlandic and particularly Pacific populations showed admixture between lineages. We estimated demographic history using a genomewide coalescence with recombination approach. The Pacific populations showed gradual population expansion starting >100 Kya, possibly reflecting persistence in cryptic refuges near the present distributional range, although we do not rule out possible influence of ancient admixture. Sharp population declines ca. 14-15 Kya were suggested to reflect founding of freshwater populations by marine ancestors. In Greenland and Northern Europe, demographic expansion started ca. 20-25 Kya coinciding with the end of the Last Glacial Maximum. In both regions, marine and freshwater populations started to show different demographic trajectories ca. 8-9 Kya, suggesting that this was the time of recolonization. In Northern Europe, this estimate was surprisingly late, but found support in subfossil evidence for presence of several freshwater fish species but not sticklebacks 12 Kya. The results demonstrate distinctly different demographic histories across geographical regions with potential consequences for adaptive processes. They also provide empirical support for previous assumptions about freshwater populations being founded independently from large, coherent marine populations, a key element in the Transporter Hypothesis invoked to explain the widespread occurrence of parallel evolution across freshwater stickleback populations.

Displaced phylogeographic signals from Gyrodactylus arcuatus, a parasite of the three-spined stickleback Gasterosteus aculeatus, suggest freshwater glacial refugia in Europe

We examined the global mitochondrial phylogeography of Gyrodactylus arcuatus, a flatworm ectoparasite of three-spined stickleback Gasterosteus aculeatus. In accordance with the suggested high divergence rate of 13%/million years, the genetic variation of the parasite was high: haplotype diversity h=0.985 and nucleotide diversity π=0.0161. The differentiation among the parasite populations was substantial (Φst=0.759), with two main allopatric clades (here termed Euro and North) accounting for 54% of the total genetic variation. The diversity center of the Euro clade was in the Baltic Sea, while the North clade was spread across the Barents and White Seas. A single haplotype within the North clade was found in the western and eastern Pacific Ocean. Divergence of main clades was estimated to be circa 200 thousand years ago. Each main clade was further divided into six distinct subclades, estimated to have diverged in isolation since 135 thousand years ago. This second division corresponds approximately to the Eemian interglacial predating the last glacial maximum. A demographic expansion of the subclades is associated with colonisation of northern Europe since the last glacial maximum, circa 15-40 thousand years ago. The parasite phylogeny is most likely explained by sequential isolated bottlenecks and expansions in numerous allopatric refugia. The postglacial intermingling and high variation in the marine parasite populations, separately in the Baltic and Barents Seas, suggest low competition of divergent parasite matrilines, coupled with a large population size and high rate of dispersal of hosts. The genetic contribution of the assumed refugial fish populations maintaining the parasite during the last glacial maximum was not detected among the marine sticklebacks, which perhaps were infected after range expansion.

The aggregate site frequency spectrum for comparative population genomic inference

Understanding how assemblages of species responded to past climate change is a central goal of comparative phylogeography and comparative population genomics, an endeavour that has increasing potential to integrate with community ecology. New sequencing technology now provides the potential to perform complex demographic inference at unprecedented resolution across assemblages of nonmodel species. To this end, we introduce the aggregate site frequency spectrum (aSFS), an expansion of the site frequency spectrum to use single nucleotide polymorphism (SNP) data sets collected from multiple, co-distributed species for assemblage-level demographic inference. We describe how the aSFS is constructed over an arbitrary number of independent population samples and then demonstrate how the aSFS can differentiate various multispecies demographic histories under a wide range of sampling configurations while allowing effective population sizes and expansion magnitudes to vary independently. We subsequently couple the aSFS with a hierarchical approximate Bayesian computation (hABC) framework to estimate degree of temporal synchronicity in expansion times across taxa, including an empirical demonstration with a data set consisting of five populations of the threespine stickleback (Gasterosteus aculeatus). Corroborating what is generally understood about the recent postglacial origins of these populations, the joint aSFS/hABC analysis strongly suggests that the stickleback data are most consistent with synchronous expansion after the Last Glacial Maximum (posterior probability = 0.99). The aSFS will have general application for multilevel statistical frameworks to test models involving assemblages and/or communities, and as large-scale SNP data from nonmodel species become routine, the aSFS expands the potential for powerful next-generation comparative population genomic inference.

Iterative development and the scope for plasticity: contrasts among trait categories in an adaptive radiation

Phenotypic plasticity can influence evolutionary change in a lineage, ranging from facilitation of population persistence in a novel environment to directing the patterns of evolutionary change. As the specific nature of plasticity can impact evolutionary consequences, it is essential to consider how plasticity is manifested if we are to understand the contribution of plasticity to phenotypic evolution. Most morphological traits are developmentally plastic, irreversible, and generally considered to be costly, at least when the resultant phenotype is mis-matched to the environment. At the other extreme, behavioral phenotypes are typically activational (modifiable on very short time scales), and not immediately costly as they are produced by constitutive neural networks. Although patterns of morphological and behavioral plasticity are often compared, patterns of plasticity of life history phenotypes are rarely considered. Here we review patterns of plasticity in these trait categories within and among populations, comprising the adaptive radiation of the threespine stickleback fish Gasterosteus aculeatus. We immediately found it necessary to consider the possibility of iterated development, the concept that behavioral and life history trajectories can be repeatedly reset on activational (usually behavior) or developmental (usually life history) time frames, offering fine tuning of the response to environmental context. Morphology in stickleback is primarily reset only in that developmental trajectories can be altered as environments change over the course of development. As anticipated, the boundaries between the trait categories are not clear and are likely to be linked by shared, underlying physiological and genetic systems.

Evolutionary Influences of Plastic Behavioral Responses Upon Environmental Challenges in an Adaptive Radiation

At the end of the 19th century, the suggestion was made by several scientists, including J. M. Baldwin, that behavioral responses to environmental change could both rescue populations from extinction (Baldwin Effect) and influence the course of subsequent evolution. Here we provide the historical and theoretical background for this argument and offer evidence of the importance of these ideas for understanding how animals (and other organisms that exhibit behavior) will respond to the rapid environmental changes caused by human activity. We offer examples from long-term research on the evolution of behavioral and other phenotypes in the adaptive radiation of the threespine stickleback fish (Gasterosteus aculeatus), a radiation in which it is possible to infer ancestral patterns of behavioral plasticity relative to the post-glacial freshwater radiation in northwestern North America, and to use patterns of parallelism and contemporary evolution to understand adaptive causes of responses to environmental modification. Our work offers insights into the complexity of cognitive responses to environmental change, and into the importance of examining multiple aspects of the phenotype simultaneously, if we are to understand how behavioral shifts contribute to the persistence of populations and to subsequent evolution. We conclude by discussing the origins of apparent novelties induced by environmental shifts, and the importance of accounting for geographic variation within species if we are to accurately anticipate the effects of anthropogenic environmental modification on the persistence and evolution of animals.

The evolution of antipredator behaviour following relaxed and reversed selection in Alaskan threespine stickleback fish

Changing environments, whether through natural or anthropogenic causes, can lead to the loss of some selective pressures ('relaxed selection') and possibly even the reinstatement of selective agents not encountered for many generations ('reversed selection'). We examined the outcome of relaxed and reversed selection in the adaptive radiation of the threespine stickleback fish, Gasterostues aculeatus L., in which isolated populations encounter a variety of predation regimes. Oceanic stickleback, which represent the ancestral founders of the freshwater radiation, encounter many piscivorous fish. Derived, freshwater populations, on the other hand, vary with respect to the presence of predators. Some populations encounter native salmonids, whereas others have not experienced predation by large fish in thousands of generations (relax-selected populations). Some relax-selected populations have had sport fish, including rainbow trout, Oncorhynchus mykiss, introduced within the past several decades (reverse-selected). We examined the behavioural responses of stickleback from three populations of each type to simulated attacks by trout and birds to determine whether relaxed and reversed selection has led to divergence in behaviour, and whether this divergence was predator specific. Fish from trout-free populations showed weak responses to trout, as predicted, but these responses were similar to those of oceanic (ancestral) populations. Fish from populations that co-occur with trout, whether native or introduced, showed elevated antipredator responses, indicating that in freshwater, trout predation selects for enhanced antipredator responses, which can evolve extremely rapidly. Comparison of laboratory-reared and wild-caught individuals suggests a combination of learned and genetic components to this variation. Responses to a model bird flyover were weakly linked to predation environment, indicating that the loss of predation by trout may partially influence the evolution of responses to birds. Our results reject the hypothesis that the consistent presence of predatory birds has been sufficient to maintain responses to piscivorous fish under periods of relaxed selection.

Landscape genetics of Schistocephalus solidus parasites in threespine stickleback (Gasterosteus aculeatus) from Alaska

The nature of gene flow in parasites with complex life cycles is poorly understood, particularly when intermediate and definitive hosts have contrasting movement potential. We examined whether the fine-scale population genetic structure of the diphyllobothriidean cestode Schistocephalus solidus reflects the habits of intermediate threespine stickleback hosts or those of its definitive hosts, semi-aquatic piscivorous birds, to better understand complex host-parasite interactions. Seventeen lakes in the Cook Inlet region of south-central Alaska were sampled, including ten in the Matanuska-Susitna Valley, five on the Kenai Peninsula, and two in the Bristol Bay drainage. We analyzed sequence variation across a 759 bp region of the mitochondrial DNA (mtDNA) cytochrome oxidase I region for 1,026 S. solidus individuals sampled from 2009-2012. We also analyzed allelic variation at 8 microsatellite loci for 1,243 individuals. Analysis of mtDNA haplotype and microsatellite genotype variation recovered evidence of significant population genetic structure within S. solidus. Host, location, and year were factors in structuring observed genetic variation. Pairwise measures revealed significant differentiation among lakes, including a pattern of isolation-by-distance. Bayesian analysis identified three distinct genotypic clusters in the study region, little admixture within hosts and lakes, and a shift in genotype frequencies over time. Evidence of fine-scale population structure in S. solidus indicates that movement of its vagile, definitive avian hosts has less influence on gene flow than expected based solely on movement potential. Observed patterns of genetic variation may reflect genetic drift, behaviors of definitive hosts that constrain dispersal, life history of intermediate hosts, and adaptive specificity of S. solidus to intermediate host genotype.

Extensive copy-number variation of young genes across stickleback populations

Duplicate genes emerge as copy-number variations (CNVs) at the population level, and remain copy-number polymorphic until they are fixed or lost. The successful establishment of such structural polymorphisms in the genome plays an important role in evolution by promoting genetic diversity, complexity and innovation. To characterize the early evolutionary stages of duplicate genes and their potential adaptive benefits, we combine comparative genomics with population genomics analyses to evaluate the distribution and impact of CNVs across natural populations of an eco-genomic model, the three-spined stickleback. With whole genome sequences of 66 individuals from populations inhabiting three distinct habitats, we find that CNVs generally occur at low frequencies and are often only found in one of the 11 populations surveyed. A subset of CNVs, however, displays copy-number differentiation between populations, showing elevated within-population frequencies consistent with local adaptation. By comparing teleost genomes to identify lineage-specific genes and duplications in sticklebacks, we highlight rampant gene content differences among individuals in which over 30% of young duplicate genes are CNVs. These CNV genes are evolving rapidly at the molecular level and are enriched with functional categories associated with environmental interactions, depicting the dynamic early copy-number polymorphic stage of genes during population differentiation.

After a locus is duplicated in a genome, individuals from a population instantaneously differ in the number of copies of this locus producing a copy-number variation (CNV). Over time, the joint effects of selection and other evolutionary forces will act to either eliminate the extra genetic copy or retain it. Depending on this evolutionary interplay, young duplications, including newly duplicated genes, can persist for millions of years as CNVs. CNVs may especially be prevalent between populations that have colonized and adapted to disparate environments in which selective pressures differ. Using whole genome sequences from several populations of three-spined sticklebacks that inhabit different environments, we find that a third of young duplicated genes are CNVs. These young CNV genes are enriched with environmental response functions and evolving rapidly at the molecular level, making them promising candidates for a role in the rapid ecological adaptation to novel environments.

Comparative analysis of Japanese three-spined stickleback clades reveals the Pacific Ocean lineage has adapted to freshwater environments while the Japan Sea has not

Divergent selection and adaptive divergence can increase phenotypic diversification amongst populations and lineages. Yet adaptive divergence between different environments, habitats or niches does not occur in all lineages. For example, the colonization of freshwater environments by ancestral marine species has triggered adaptive radiation and phenotypic diversification in some taxa but not in others. Studying closely related lineages differing in their ability to diversify is an excellent means of understanding the factors promoting and constraining adaptive evolution. A well-known example of the evolution of increased phenotypic diversification following freshwater colonization is the three-spined stickleback. Two closely related stickleback lineages, the Pacific Ocean and the Japan Sea occur in Japan. However, Japanese freshwater stickleback populations are derived from the Pacific Ocean lineage only, suggesting the Japan Sea lineage is unable to colonize freshwater. Using stable isotope data and trophic morphology, we first show higher rates of phenotypic and ecological diversification between marine and freshwater populations within the Pacific Ocean lineage, confirming adaptive divergence has occurred between the two lineages and within the Pacific Ocean lineage but not in the Japan Sea lineage. We further identified consistent divergence in diet and foraging behaviour between marine forms from each lineage, confirming Pacific Ocean marine sticklebacks, from which all Japanese freshwater populations are derived, are better adapted to freshwater environments than Japan Sea sticklebacks. We suggest adaptive divergence between ancestral marine populations may have played a role in constraining phenotypic diversification and adaptive evolution in Japanese sticklebacks.

Reversed brain size sexual dimorphism accompanies loss of parental care in white sticklebacks

Uncovering factors that shape variation in brain morphology remains a major challenge in evolutionary biology. Recently, it has been shown that brain size is positively associated with level of parental care behavior in various taxa. One explanation for this pattern is that the cognitive demands of performing complex parental care may require increased brain size. This idea is known as the parental brain hypothesis (PBH). We set out to test the predictions of this hypothesis in wild populations of threespine stickleback (Gasterosteus aculeatus). These fish are commonly known to exhibit (1) uniparental male care and (2) sexual dimorphism in brain size (males>females). To test the PBH, we took advantage of the existence of closely related populations of stickleback that display variation in parental care behavior: common marine threespine sticklebacks (uniparental male care) and white threespine sticklebacks (no care). To begin, we quantified genetic differentiation among two common populations and three white populations from Nova Scotia. We found overall low differentiation among populations, although F_ST was increased in between-type comparisons. We then measured the brain weights of males and females from all five populations along with two additional common populations from British Columbia. We found that sexual dimorphism in brain size is reversed in white stickleback populations: males have smaller brains than females. Thus, while several alternatives need to be ruled out, the PBH appears to be a reasonable explanation for sexual dimorphism in brain size in threespine sticklebacks.

Unique mitochondrial DNA lineages in Irish stickleback populations: cryptic refugium or rapid recolonization?

Repeated recolonization of freshwater environments following Pleistocene glaciations has played a major role in the evolution and adaptation of anadromous taxa. Located at the western fringe of Europe, Ireland and Britain were likely recolonized rapidly by anadromous fishes from the North Atlantic following the last glacial maximum (LGM). While the presence of unique mitochondrial haplotypes in Ireland suggests that a cryptic northern refugium may have played a role in recolonization, no explicit test of this hypothesis has been conducted. The three-spined stickleback is native and ubiquitous to aquatic ecosystems throughout Ireland, making it an excellent model species with which to examine the biogeographical history of anadromous fishes in the region. We used mitochondrial and microsatellite markers to examine the presence of divergent evolutionary lineages and to assess broad-scale patterns of geographical clustering among postglacially isolated populations. Our results confirm that Ireland is a region of secondary contact for divergent mitochondrial lineages and that endemic haplotypes occur in populations in Central and Southern Ireland. To test whether a putative Irish lineage arose from a cryptic Irish refugium, we used approximate Bayesian computation (ABC). However, we found no support for this hypothesis. Instead, the Irish lineage likely diverged from the European lineage as a result of postglacial isolation of freshwater populations by rising sea levels. These findings emphasize the need to rigorously test biogeographical hypothesis and contribute further evidence that postglacial processes may have shaped genetic diversity in temperate fauna.

Genetic architecture of parallel pelvic reduction in ninespine sticklebacks

Teleost fish genomes are known to be evolving faster than those of other vertebrate taxa. Thus, fish are suited to address the extent to which the same vs. different genes are responsible for similar phenotypic changes in rapidly evolving genomes of evolutionary independent lineages. To gain insights into the genetic basis and evolutionary processes behind parallel phenotypic changes within and between species, we identified the genomic regions involved in pelvic reduction in Northern European ninespine sticklebacks (Pungitius pungitius) and compared them to those of North American ninespine and threespine sticklebacks (Gasterosteus aculeatus). To this end, we conducted quantitative trait locus (QTL) mapping using 283 F₂ progeny from an interpopulation cross. Phenotypic analyses indicated that pelvic reduction is a recessive trait and is inherited in a simple Mendelian fashion. Significant QTL for pelvic spine and girdle lengths were identified in the region of the Pituitary homeobox transcription factor 1 (Pitx1) gene, also responsible for pelvic reduction in threespine sticklebacks. The fact that no QTL was observed in the region identified in the mapping study of North American ninespine sticklebacks suggests that an alternative QTL for pelvic reduction has emerged in this species within the past 1.6 million years after the split between Northern European and North American populations. In general, our study provides empirical support for the view that alternative genetic mechanisms that lead to similar phenotypes can evolve over short evolutionary time scales.

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.

Phylogeography and adaptation genetics of stickleback from the Haida Gwaii archipelago revealed using genome-wide single nucleotide polymorphism genotyping

Threespine stickleback populations are model systems for studying adaptive evolution and the underlying genetics. In lakes on the Haida Gwaii archipelago (off western Canada), stickleback have undergone a remarkable local radiation and show phenotypic diversity matching that seen throughout the species distribution. To provide a historical context for this radiation, we surveyed genetic variation at >1000 single nucleotide polymorphism (SNP) loci in stickleback from over 100 populations. SNPs included markers evenly distributed throughout genome and candidate SNPs tagging adaptive genomic regions. Based on evenly distributed SNPs, the phylogeographic pattern differs substantially from the disjunct pattern previously observed between two highly divergent mtDNA lineages. The SNP tree instead shows extensive within watershed population clustering and different watersheds separated by short branches deep in the tree. These data are consistent with separate colonizations of most watersheds, despite underlying genetic connections between some independent drainages. This supports previous suppositions that morphological diversity observed between watersheds has been shaped independently, with populations exhibiting complete loss of lateral plates and giant size each occurring in several distinct clades. Throughout the archipelago, we see repeated selection of SNPs tagging candidate freshwater adaptive variants at several genomic regions differentiated between marine-freshwater populations on a global scale (e.g. EDA, Na/K ATPase). In estuarine sites, both marine and freshwater allelic variants were commonly detected. We also found typically marine alleles present in a few freshwater lakes, especially those with completely plated morphology. These results provide a general model for postglacial colonization of freshwater habitat by sticklebacks and illustrate the tremendous potential of genome-wide SNP data sets hold for resolving patterns and processes underlying recent adaptive divergences.

Progressive recombination suppression and differentiation in recently evolved neo-sex chromosomes

Recombination suppression leads to the structural and functional differentiation of sex chromosomes and is thus a crucial step in the process of sex chromosome evolution. Despite extensive theoretical work, the exact processes and mechanisms of recombination suppression and differentiation are not well understood. In threespine sticklebacks (Gasterosteus aculeatus), a different sex chromosome system has recently evolved by a fusion between the Y chromosome and an autosome in the Japan Sea lineage, which diverged from the ancestor of other lineages approximately 2 Ma. We investigated the evolutionary dynamics and differentiation processes of sex chromosomes based on comparative analyses of these divergent lineages using 63 microsatellite loci. Both chromosome-wide differentiation patterns and phylogenetic inferences with X and Y alleles indicated that the ancestral sex chromosomes were extensively differentiated before the divergence of these lineages. In contrast, genetic differentiation appeared to have proceeded only in a small region of the neo-sex chromosomes. The recombination maps constructed for the Japan Sea lineage indicated that recombination has been suppressed or reduced over a large region spanning the ancestral and neo-sex chromosomes. Chromosomal regions exhibiting genetic differentiation and suppressed or reduced recombination were detected continuously and sequentially in the neo-sex chromosomes, suggesting that differentiation has gradually spread from the fusion point following the extension of recombination suppression. Our study illustrates an ongoing process of sex chromosome differentiation, providing empirical support for the theoretical model postulating that recombination suppression and differentiation proceed in a gradual manner in the very early stage of sex chromosome evolution.

Parallel and nonparallel ecological, morphological and genetic divergence in lake-stream stickleback from a single catchment

Parallel phenotypic evolution in similar environments has been well studied in evolutionary biology; however, comparatively little is known about the influence of determinism and historical contingency on the nature, extent and generality of this divergence. Taking advantage of a novel system containing multiple lake-stream stickleback populations, we examined the extent of ecological, morphological and genetic divergence between three-spined stickleback present in parapatric environments. Consistent with other lake-stream studies, we found a shift towards a deeper body and shorter gill rakers in stream fish. Morphological shifts were concurrent with changes in diet, indicated by both stable isotope and stomach contents analysis. Performing a multivariate test for shared and unique components of evolutionary response to the distance gradient from the lake, we found a strong signature of parallel adaptation. Nonparallel divergence was also present, attributable mainly to differences between river locations. We additionally found evidence of genetic substructuring across five lake-stream transitions, indicating that some level of reproductive isolation occurs between populations in these habitats. Strong correlations between pairwise measures of morphological, ecological and genetic distance between lake and stream populations supports the hypothesis that divergent natural selection between habitats drives adaptive divergence and reproductive isolation. Lake-stream stickleback divergence in Lough Neagh provides evidence for the deterministic role of selection and supports the hypothesis that parallel selection in similar environments may initiate parallel speciation.

Shifts in morphology and diet of non-native sticklebacks introduced into Japanese crater lakes

An increasing number of exotic animals are causing ecological problems. Therefore, for better ecosystem management, it is important to understand how exotic species colonize and adapt to novel environments. The threespine sticklebacks (Gasterosteus aculeatus) can be a good vertebrate model system to explore the ecological and genetic mechanisms of adaptation not only in natural populations, but also in non-native populations. Although morphological changes have been documented in several introduced populations of stickleback, little is known about the dietary changes during colonization into novel environments. Here, we investigated the morphological and dietary changes of exotic threespine stickleback populations introduced into three Japanese crater lakes (Lake Towada, Lake Kussharo, and Lake Shikotsu). Sticklebacks were introduced into the crater lakes likely along with salmonids transplanted for aquaculture. The stickleback population in Lake Kussharo had multiple mitochondrial haplotypes and had larger phenotypic variances than other crater lake stickleback populations that had only one mitochondrial haplotype. Compilation of historical data on the morphology and stomach contents of the Lake Towada stickleback population showed that substantial shifts in body size and stomach contents occurred after colonization. Some of these changes may be related to an outbreak of the Schistocephalus parasite. These results suggest that sticklebacks can change their morphology and trophic ecology when they colonize novel environments. Therefore, extreme care should be taken when salmonids are transported between watersheds for aquaculture and that long-term monitoring of exotic species is essential for ecosystem management. In addition, further genetic studies on phenotypic changes in crater lake sticklebacks would help elucidate the genetic mechanisms underlying the adaptation of exotic fishes to novel environments.

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.

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.

Distinct lineages of Schistocephalus parasites in threespine and ninespine stickleback hosts revealed by DNA sequence analysis

Parasitic interactions are often part of complex networks of interspecific relationships that have evolved in biological communities. Despite many years of work on the evolution of parasitism, the likelihood that sister taxa of parasites can co-evolve with their hosts to specifically infect two related lineages, even when those hosts occur sympatrically, is still unclear. Furthermore, when these specific interactions occur, the molecular and physiological basis of this specificity is still largely unknown. The presence of these specific parasitic relationships can now be tested using molecular markers such as DNA sequence variation. Here we test for specific parasitic relationships in an emerging host-parasite model, the stickleback-Schistocephalus system. Threespine and ninespine stickleback fish are intermediate hosts for Schistocephalus cestode parasites that are phenotypically very similar and have nearly identical life cycles through plankton, stickleback, and avian hosts. We analyzed over 2000 base pairs of COX1 and NADH1 mitochondrial DNA sequences in 48 Schistocephalus individuals collected from threespine and ninespine stickleback hosts from disparate geographic regions distributed across the Northern Hemisphere. Our data strongly support the presence of two distinct clades of Schistocephalus, each of which exclusively infects either threespine or ninespine stickleback. These clades most likely represent different species that diverged soon after the speciation of their stickleback hosts. In addition, genetic structuring exists among Schistocephalus taken from threespine stickleback hosts from Alaska, Oregon and Wales, although it is much less than the divergence between hosts. Our findings emphasize that biological communities may be even more complex than they first appear, and beg the question of what are the ecological, physiological, and genetic factors that maintain the specificity of the Schistocephalus parasites and their stickleback hosts.

Phylogeography of ninespine sticklebacks (Pungitius pungitius) in North America: glacial refugia and the origins of adaptive traits

The current geographical distribution of the ninespine stickleback (Pungitius pungitius) was shaped in large part by the glaciation events of the Pleistocene epoch (2.6 Mya-1 Kya). Previous efforts to elucidate the phylogeographical history of the ninespine stickleback in North America have focused on a limited set of morphological traits, some of which are likely subject to widespread convergent evolution, thereby potentially obscuring relationships among populations. In this study, we used genetic information from both mitochondrial DNA (mtDNA) sequences and nuclear microsatellite markers to determine the phylogenetic relationships among ninespine stickleback populations. We found that ninespine sticklebacks in North America probably dispersed from at least three glacial refugia-the Mississippi, Bering, and Atlantic refugia-not two as previously thought. However, by applying a molecular clock to our mtDNA data, we found that these three groups diverged long before the most recent glacial period. Our new phylogeny serves as a critical framework for examining the evolution of derived traits in this species, including adaptive phenotypes that evolved multiple times in different lineages. In particular, we inferred that loss of the pelvic (hind fin) skeleton probably evolved independently in populations descended from each of the three putative North American refugia.

Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags

Next-generation sequencing technology provides novel opportunities for gathering genome-scale sequence data in natural populations, laying the empirical foundation for the evolving field of population genomics. Here we conducted a genome scan of nucleotide diversity and differentiation in natural populations of threespine stickleback (Gasterosteus aculeatus). We used Illumina-sequenced RAD tags to identify and type over 45,000 single nucleotide polymorphisms (SNPs) in each of 100 individuals from two oceanic and three freshwater populations. Overall estimates of genetic diversity and differentiation among populations confirm the biogeographic hypothesis that large panmictic oceanic populations have repeatedly given rise to phenotypically divergent freshwater populations. Genomic regions exhibiting signatures of both balancing and divergent selection were remarkably consistent across multiple, independently derived populations, indicating that replicate parallel phenotypic evolution in stickleback may be occurring through extensive, parallel genetic evolution at a genome-wide scale. Some of these genomic regions co-localize with previously identified QTL for stickleback phenotypic variation identified using laboratory mapping crosses. In addition, we have identified several novel regions showing parallel differentiation across independent populations. Annotation of these regions revealed numerous genes that are candidates for stickleback phenotypic evolution and will form the basis of future genetic analyses in this and other organisms. This study represents the first high-density SNP-based genome scan of genetic diversity and differentiation for populations of threespine stickleback in the wild. These data illustrate the complementary nature of laboratory crosses and population genomic scans by confirming the adaptive significance of previously identified genomic regions, elucidating the particular evolutionary and demographic history of such regions in natural populations, and identifying new genomic regions and candidate genes of evolutionary significance.

Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags

Next-generation sequencing technology provides novel opportunities for gathering genome-scale sequence data in natural populations, laying the empirical foundation for the evolving field of population genomics. Here we conducted a genome scan of nucleotide diversity and differentiation in natural populations of threespine stickleback (Gasterosteus aculeatus). We used Illumina-sequenced RAD tags to identify and type over 45,000 single nucleotide polymorphisms (SNPs) in each of 100 individuals from two oceanic and three freshwater populations. Overall estimates of genetic diversity and differentiation among populations confirm the biogeographic hypothesis that large panmictic oceanic populations have repeatedly given rise to phenotypically divergent freshwater populations. Genomic regions exhibiting signatures of both balancing and divergent selection were remarkably consistent across multiple, independently derived populations, indicating that replicate parallel phenotypic evolution in stickleback may be occurring through extensive, parallel genetic evolution at a genome-wide scale. Some of these genomic regions co-localize with previously identified QTL for stickleback phenotypic variation identified using laboratory mapping crosses. In addition, we have identified several novel regions showing parallel differentiation across independent populations. Annotation of these regions revealed numerous genes that are candidates for stickleback phenotypic evolution and will form the basis of future genetic analyses in this and other organisms. This study represents the first high-density SNP-based genome scan of genetic diversity and differentiation for populations of threespine stickleback in the wild. These data illustrate the complementary nature of laboratory crosses and population genomic scans by confirming the adaptive significance of previously identified genomic regions, elucidating the particular evolutionary and demographic history of such regions in natural populations, and identifying new genomic regions and candidate genes of evolutionary significance.