Chromosomal Fusions Facilitate Adaptation to Divergent Environments in Threespine Stickleback

Chromosomal fusions are hypothesized to facilitate adaptation to divergent environments, both by bringing together previously unlinked adaptive alleles and by creating regions of low recombination that facilitate the linkage of adaptive alleles; but, there is little empirical evidence to support this hypothesis. Here, we address this knowledge gap by studying threespine stickleback (Gasterosteus aculeatus), in which ancestral marine fish have repeatedly adapted to freshwater across the northern hemisphere. By comparing the threespine and ninespine stickleback (Pungitius pungitius) genomes to a de novo assembly of the fourspine stickleback (Apeltes quadracus) and an outgroup species, we find two chromosomal fusion events involving the same chromosomes have occurred independently in the threespine and ninespine stickleback lineages. On the fused chromosomes in threespine stickleback, we find an enrichment of quantitative trait loci underlying traits that contribute to marine versus freshwater adaptation. By comparing whole-genome sequences of freshwater and marine threespine stickleback populations, we also find an enrichment of regions under divergent selection on these two fused chromosomes. There is elevated genetic diversity within regions under selection in the freshwater population, consistent with a simulation study showing that gene flow can increase diversity in genomic regions associated with local adaptation and our demographic models showing gene flow between the marine and freshwater populations. Integrating our results with previous studies, we propose that these fusions created regions of low recombination that enabled the formation of adaptative clusters, thereby facilitating freshwater adaptation in the face of recurrent gene flow between marine and freshwater threespine sticklebacks.     

Genetic differentiation of freshwater and marine sticklebacks, (Gasterosteus aculeatus) of Eastern Europe

Electrophoretic variation of 13 presumptive enzymatic loci was studied in sticklebacks (Gasterosteus aculeatus) at 6 marine (Baltic) and 13 freshwater sites from Poland. Marine samples are significantly more variable (proportion of polymorphic loci, p # 0.38) than freshwater Samles (p # 0.12). Also the mean number of alleles per locus is significantly higher in the marine samples (n # 1.54) than in freshwater samples (n # 1.14). There is however no significant difference in values of mean heterozygosity between marine and freshwater sticklebacks. Some alleles are found in Baltic sticklebacks only, but all these occur in very low frequencies. Genetic distance between samples is low ranging from D = 0.000053 to 0.031410. However D is significant smaller between marine samples than between freshwater Samles. We suggest than this probaily is the result of larger genetic subdivision of freshwater populations. The low genetic divergence of marine and freshwater populations of the stickleback from poland might be expfained by the recent postglacial colonization of this area from a common refugium.     

Rapid evolution of cold tolerance in stickleback

Climate change is predicted to lead to increased average temperatures and greater intensity and frequency of high and low temperature extremes, but the evolutionary consequences for biological communities are not well understood. Studies of adaptive evolution of temperature tolerance have typically involved correlative analyses of natural populations or artificial selection experiments in the laboratory. Field experiments are required to provide estimates of the timing and strength of natural selection, enhance understanding of the genetics of adaptation and yield insights into the mechanisms driving evolutionary change. Here, we report the experimental evolution of cold tolerance in natural populations of threespine stickleback fish (Gasterosteus aculeatus). We show that freshwater sticklebacks are able to tolerate lower minimum temperatures than marine sticklebacks and that this difference is heritable. We transplanted marine sticklebacks to freshwater ponds and measured the rate of evolution after three generations in this environment. Cold tolerance evolved at a rate of 0.63 haldanes to a value 2.5°C lower than that of the ancestral population, matching values found in wild freshwater populations. Our results suggest that cold tolerance is under strong selection and that marine sticklebacks carry sufficient genetic variation to adapt to changes in temperature over remarkably short time scales.     

Population Structure Limits Parallel Evolution in Sticklebacks

Population genetic theory predicts that small effective population sizes (N_e) and restricted gene flow limit the potential for local adaptation. In particular, the probability of evolving similar phenotypes based on shared genetic mechanisms (i.e., parallel evolution), is expected to be reduced. We tested these predictions in a comparative genomic study of two ecologically similar and geographically codistributed stickleback species (viz. Gasterosteus aculeatus and Pungitius pungitius). We found that P. pungitius harbors less genetic diversity and exhibits higher levels of genetic differentiation and isolation-by-distance than G. aculeatus. Conversely, G. aculeatus exhibits a stronger degree of genetic parallelism across freshwater populations than P. pungitius: 2,996 versus 379 single nucleotide polymorphisms located within 26 versus 9 genomic regions show evidence of selection in multiple freshwater populations of G. aculeatus and P. pungitius, respectively. Most regions involved in parallel evolution in G. aculeatus showed increased levels of divergence, suggestive of selection on ancient haplotypes. In contrast, haplotypes involved in freshwater adaptation in P. pungitius were younger. In accordance with theory, the results suggest that connectivity and genetic drift play crucial roles in determining the levels and geographic distribution of standing genetic variation, providing evidence that population subdivision limits local adaptation and therefore also the likelihood of parallel evolution.     

Non‐parallel divergence across freshwater and marine three‐spined stickleback Gasterosteus aculeatus populations

This work investigated whether multiple freshwater populations of three‐spined stickleback Gasterosteus aculeatus in different freshwater catchments in the Jutland Peninsula, Denmark, derived from the same marine populations show repeated adaptive responses. A total of 327 G. aculeatus collected at 13 sampling locations were screened for genetic variation using a combination of 70 genes putatively under selection and 26 neutral genes along with a marker linked to the ectodysplasin gene (eda), which is strongly correlated with plate armour morphs in the species. A highly significant genetic differentiation was found that was higher among different freshwater samples than between marine–freshwater samples. Tests for selection between marine and freshwater populations showed a very low degree of parallelism and no single nucleotide polymorphism was detected as outlier in all freshwater–marine pairwise comparisons, including the eda. This suggests that G. aculeatus is not necessarily the prime example of parallel local adaptation suggested in much of the literature and that important exceptions exist (i.e. the Jutland Peninsula). While marine populations in the results described here showed a high phenotype–genotype correlation at eda, a low association was found for most of the freshwater populations. The most extreme case was found in the freshwater Lake Hald where all low‐plated phenotypes were either homozygotes for the allele supposed to be associated with completely plated morphs or heterozygotes, but none were homozygotes for the putative low‐plated allele. Re‐examination of data from seven G. aculeatus studies agrees in showing a high but partial association between phenotype–genotype at eda in G. aculeatus freshwater populations and that mismatches occur everywhere in the European regions studied (higher in some areas, i.e. Denmark). This is independent of the eda marker used.     

Partially repeatable genetic basis of benthic adaptation in threespine sticklebacks

The extent to which convergent adaptation to similar ecological niches occurs by a predictable genetic basis remains a fundamental question in biology. Threespine stickleback fish have undergone an adaptive radiation in which ancestral oceanic populations repeatedly colonized and adapted to freshwater habitats. In multiple lakes in British Columbia, two different freshwater ecotypes have evolved: a deep‐bodied benthic form adapted to forage near the lake substrate, and a narrow‐bodied limnetic form adapted to forage in open water. Here, we use genome‐wide linkage mapping in marine × benthic F2 genetic crosses to test the extent of shared genomic regions underlying benthic adaptation in three benthic populations. We identify at least 100 Quantitative Trait Loci (QTL) harboring genes influencing skeletal morphology. The majority of QTL (57%) are unique to one cross. However, four genomic regions affecting eight craniofacial and armor phenotypes are found in all three benthic populations. We find that QTL are clustered in the genome and overlapping QTL regions are enriched for genomic signatures of natural selection. These findings suggest that benthic adaptation has occurred via both parallel and nonparallel genetic changes.     

Rapid Morphological Changes in Threespine Stickleback, Gasterosteus aculeatus, in Freshwater

Synopsis Freshwater and marine threespine stickleback, Gasterosteus aculeatus, differ remarkably in armour plate number and body shape, although differences in other morphological characters are also common. Most freshwater populations have apparently evolved after isolation of marine sticklebacks in freshwater. After colonisation of freshwater habitats, they show rapid morphological changes and associated genetic isolation within as few as eight generations. I transferred fish from marine tide pools to two isolated freshwater ponds, differing in habitat characteristics, at the beginning of the breeding season, when females had ripe ovaries and males had breeding coloration. The first generation fish that I sampled from the ponds had significantly fewer armour plates than their marine ancestors and differed in shape. I also found some significant differences between fish sampled from the larger pond and those from a smaller, adjacent pond. This extremely rapid morphological divergence suggests that either the marine sticklebacks were highly phenotypically plastic or that there was very strong natural selection acting on the first generation within freshwater habitats.     

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.     

Frequency of Ectodysplasin alleles and limited introgression between sympatric threespine stickleback populations

The threespine stickleback (Gasterosteus aculeatus) is primitively an anadromous or resident marine species but has repeatedly colonized fresh water, where predictable phenotypic divergence usually occurs rapidly. A conspicuous element of this divergence is change of the number and position of lateral armor plates from about 33 that cover the entire flank (complete) to <10 anterior plates (low). This difference is caused primarily by variation at the Ectodysplasin (Eda) locus. The low Eda allele appears to be rarer in two geographically adjacent anadromous populations from Cook Inlet, Alaska than in most marine or anadromous populations reported from elsewhere, and there is no evidence of elevated gene flow for Eda between anadromous and resident lake threespine stickleback populations that breed in sympatry. However, the two anadromous populations are divergent for the frequencies of two complete Eda alleles. It is not clear how monomorphic low-plated freshwater populations in Cook Inlet have almost invariably acquired ancestral low Eda alleles from anadromous ancestors in which this allele appears to be extremely rare.     

Heritability of DNA methylation in threespine stickleback (Gasterosteus aculeatus)

Epigenetic mechanisms underlying phenotypic change are hypothesized to contribute to population persistence and adaptation in the face of environmental change. To date, few studies have explored the heritability of intergenerationally stable methylation levels in natural populations, and little is known about the relative contribution of cis- and trans-regulatory changes to methylation variation. Here, we explore the heritability of DNA methylation, and conduct methylation quantitative trait loci (meQTLs) analysis to investigate the genetic architecture underlying methylation variation between marine and freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We quantitatively measured genome-wide DNA methylation in fin tissue using reduced representation bisulfite sequencing of F1 and F2 crosses, and their marine and freshwater source populations. We identified cytosines (CpG sites) that exhibited stable methylation levels across generations. We found that additive genetic variance explained an average of 24–35% of the methylation variance, with a number of CpG sites possibly autonomous from genetic control. We also detected both cis- and trans-meQTLs, with only trans-meQTLs overlapping with previously identified genomic regions of high differentiation between marine and freshwater ecotypes. Finally, we identified the genetic architecture underlying two key CpG sites that were differentially methylated between ecotypes. These findings demonstrate a potential role for DNA methylation in facilitating adaptation to divergent environments and improve our understanding of the heritable basis of population epigenomic variation.     

Natural selection affects multiple aspects of genetic variation at putatively neutral sites across the human genome

A major question in evolutionary biology is how natural selection has shaped patterns of genetic variation across the human genome. Previous work has documented a reduction in genetic diversity in regions of the genome with low recombination rates. However, it is unclear whether other summaries of genetic variation, like allele frequencies, are also correlated with recombination rate and whether these correlations can be explained solely by negative selection against deleterious mutations or whether positive selection acting on favorable alleles is also required. Here we attempt to address these questions by analyzing three different genome-wide resequencing datasets from European individuals. We document several significant correlations between different genomic features. In particular, we find that average minor allele frequency and diversity are reduced in regions of low recombination and that human diversity, human-chimp divergence, and average minor allele frequency are reduced near genes. Population genetic simulations show that either positive natural selection acting on favorable mutations or negative natural selection acting against deleterious mutations can explain these correlations. However, models with strong positive selection on nonsynonymous mutations and little negative selection predict a stronger negative correlation between neutral diversity and nonsynonymous divergence than observed in the actual data, supporting the importance of negative, rather than positive, selection throughout the genome. Further, we show that the widespread presence of weakly deleterious alleles, rather than a small number of strongly positively selected mutations, is responsible for the correlation between neutral genetic diversity and recombination rate. This work suggests that natural selection has affected multiple aspects of linked neutral variation throughout the human genome and that positive selection is not required to explain these observations.     

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.     

Adaptive evolution of lateral plates in three‐spined stickleback Gasterosteus aculeatus: a case study in functional analysis of natural variation

The three‐spined stickleback Gasterosteus aculeatus is a model species for studying questions in ecology and evolution. The rapid diversification of G. aculeatus in post‐glacial freshwater environments, combined with recently developed molecular tools, provides a unique opportunity to study the functional basis of fitness variation in natural populations. In derived freshwater populations, a number of morphological traits have diverged in parallel from the marine ancestral state, including the number of lateral armour plates. Evolution of reduced armour in freshwater populations is due to positive selection from both abiotic and biotic mechanisms. The major effect gene (ectodysplasin‐A or Eda), along with several minor effect genetic regions, has recently been shown to control lateral plate variation. Field experiments have further determined the fitness consequences of allelic variation at the major effect locus. This work helps elucidate the mechanisms connecting genetic variation with phenotypic variation and fitness in the wild, a synthesis that should be applicable to many other phenotypic traits and species of fishes.     

Should I stay or should I go? The Ectodysplasin locus is associated with behavioural differences in threespine stickleback

Adaptive divergence may be facilitated if morphological and behavioural traits associated with local adaptation share the same genetic basis. It is therefore important to determine whether genes underlying adaptive morphological traits are associated with variation in behaviour in natural populations. Positive selection on low-armour alleles at the Ectodysplasin (Eda) locus in threespine stickleback has led to the repeated evolution of reduced armour, following freshwater colonization by fully armoured marine sticklebacks. This adaptive divergence in armour between marine and freshwater populations would be facilitated if the low allele conferred a behavioural preference for freshwater environments. We experimentally tested whether the low allele is associated with preference for freshwater by measuring the preference of each Eda genotype for freshwater versus saltwater after acclimation to either salinity. We found no association between the Eda low allele and preference for freshwater. Instead, the low allele was significantly associated with a reduced preference for the acclimation environment. This behaviour may facilitate the colonization of freshwater habitats from the sea, but could also hinder local adaptation by promoting migration of low alleles between marine and freshwater environments.     

Signatures of selection in the three‐spined stickleback along a small‐scale brackish water – freshwater transition zone

Local adaptation is often obvious when gene flow is impeded, such as observed at large spatial scales and across strong ecological contrasts. However, it becomes less certain at small scales such as between adjacent populations or across weak ecological contrasts, when gene flow is strong. While studies on genomic adaptation tend to focus on the former, less is known about the genomic targets of natural selection in the latter situation. In this study, we investigate genomic adaptation in populations of the three‐spined stickleback Gasterosteus aculeatus L. across a small‐scale ecological transition with salinities ranging from brackish to fresh. Adaptation to salinity has been repeatedly demonstrated in this species. A genome scan based on 87 microsatellite markers revealed only few signatures of selection, likely owing to the constraints that homogenizing gene flow puts on adaptive divergence. However, the detected loci appear repeatedly as targets of selection in similar studies of genomic adaptation in the three‐spined stickleback. We conclude that the signature of genomic selection in the face of strong gene flow is weak, yet detectable. We argue that the range of studies of genomic divergence should be extended to include more systems characterized by limited geographical and ecological isolation, which is often a realistic setting in nature.     

Recurrent selection explains parallel evolution of genomic regions of high relative but low absolute differentiation in a ring species

Recent technological developments allow investigation of the repeatability of evolution at the genomic level. Such investigation is particularly powerful when applied to a ring species, in which spatial variation represents changes during the evolution of two species from one. We examined genomic variation among three subspecies of the greenish warbler ring species, using genotypes at 13 013 950 nucleotide sites along a new greenish warbler consensus genome assembly. Genomic regions of low within‐group variation are remarkably consistent between the three populations. These regions show high relative differentiation but low absolute differentiation between populations. Comparisons with outgroup species show the locations of these peaks of relative differentiation are not well explained by phylogenetically conserved variation in recombination rates or selection. These patterns are consistent with a model in which selection in an ancestral form has reduced variation at some parts of the genome, and those same regions experience recurrent selection that subsequently reduces variation within each subspecies. The degree of heterogeneity in nucleotide diversity is greater than explained by models of background selection, but is consistent with selective sweeps. Given the evidence that greenish warblers have had both population differentiation for a long period of time and periods of gene flow between those populations, we propose that some genomic regions underwent selective sweeps over a broad geographic area followed by within‐population selection‐induced reductions in variation. An important implication of this ‘sweep‐before‐differentiation’ model is that genomic regions of high relative differentiation may have moved among populations more recently than other genomic regions.     

Body armour and lateral‐plate reduction in freshwater three‐spined stickleback Gasterosteus aculeatus: adaptations to a different buoyancy regime?

Several factors related to buoyancy were compared between one marine and two freshwater populations of three‐spined stickleback Gasterosteus aculeatus. Fish from all three populations had buoyancy near to neutral to the ambient water. This showed that neither marine nor freshwater G. aculeatus used swimming and hydrodynamic lift to prevent sinking. Comparing the swimbladder volumes showed that freshwater completely plated G. aculeatus had a significantly larger swimbladder volume than both completely plated marine and low‐plated freshwater G. aculeatus. Furthermore, body tissue density was lower in low‐plated G. aculeatus than in the completely plated marine and freshwater fish. The results show that G. aculeatus either reduce tissue density or increase swimbladder volume to adapt to lower water density. Mass measurements of lateral plates and pelvis showed that loss of body armour in low‐plated G. aculeatus could explain the tissue density difference between low‐plated and completely plated G. aculeatus. This suggests that the common occurrence of plate and armour reduction in freshwater G. aculeatus populations can be an adaptation to a lower water density.     

A low-density SNP array for analyzing differential selection in freshwater and marine populations of threespine stickleback (Gasterosteus aculeatus)

Background

The threespine stickleback (Gasterosteus aculeatus) has become an important model species for studying both contemporary and parallel evolution. In particular, differential adaptation to freshwater and marine environments has led to high differentiation between freshwater and marine stickleback populations at the phenotypic trait of lateral plate morphology and the underlying candidate gene Ectodysplacin (EDA). Many studies have focused on this trait and candidate gene, although other genes involved in marine-freshwater adaptation may be equally important. In order to develop a resource for rapid and cost efficient analysis of genetic divergence between freshwater and marine sticklebacks, we generated a low-density SNP (Single Nucleotide Polymorphism) array encompassing markers of chromosome regions under putative directional selection, along with neutral markers for background.

Results

RAD (Restriction site Associated DNA) sequencing of sixty individuals representing two freshwater and one marine population led to the identification of 33,993 SNP markers. Ninety-six of these were chosen for the low-density SNP array, among which 70 represented SNPs under putatively directional selection in freshwater vs. marine environments, whereas 26 SNPs were assumed to be neutral. Annotation of these regions revealed several genes that are candidates for affecting stickleback phenotypic variation, some of which have been observed in previous studies whereas others are new.

Conclusions

We have developed a cost-efficient low-density SNP array that allows for rapid screening of polymorphisms in threespine stickleback. The array provides a valuable tool for analyzing adaptive divergence between freshwater and marine stickleback populations beyond the well-established candidate gene Ectodysplacin (EDA).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-867) contains supplementary material, which is available to authorized users.     

MULTIPLE EVOLUTIONARY PATHWAYS TO DECREASED LATERAL PLATE COVERAGE IN FRESHWATER THREESPINE STICKLEBACKS

Adaptation to novel environments can be based either on standing genetic variation or variation attributable to new mutations. When standing genetic variation for a functional adaptation is lacking, and variation due to new mutations is not yet available, adaptation is possible only through alternative functional solutions. Reduction in the number of bony lateral plates as an adaptation to freshwater colonization by marine threespine sticklebacks (Gasterosteus aculeatus) has occurred in numerous independent cases through allelic substitution in the ectodysplasin‐a (Eda) gene. Studying the phenotypic and genetic variation in plate number and size in five marine and six freshwater threespine stickleback populations, we found that when variation in Eda was limiting (i.e., alleles associated with the low‐plate morph were missing or in extremely low frequency), plate number reduction did not take place in freshwater populations, but reduced lateral plate coverage was achieved by a reduction in the size of lateral plates. Our results suggest that this phenotypically and genetically discrete "small‐plated" threespine stickleback—which is the dominant form in three northern European freshwater populations—may be functionally equivalent to the low‐plated morph and hence, serve as an example of convergent evolution toward functional similarity in the face of genetic constraints.     

Genome‐wide analysis of allele frequency change in sunflower crop–wild hybrid populations evolving under natural conditions

Crop‐wild hybridization occurs in numerous plant species and could alter the genetic structure and evolutionary dynamics of wild populations. Studying crop‐derived alleles in wild populations is also relevant to assessing/mitigating the risks associated with transgene escape. To date, crop‐wild hybridization has generally been examined via short‐term studies, typically within a single generation, focusing on few traits or genetic markers. Little is known about patterns of selection on crop‐derived alleles over multiple generations, particularly at a genome‐wide scale. Here, we documented patterns of natural selection in an experimental crop × wild sunflower population that was allowed to evolve under natural conditions for two generations at two locations. Allele frequencies at a genome‐wide collection of SNPs were tracked across generations, and a common garden experiment was conducted to compare trait means between generations. These data allowed us to identify instances of selection on crop‐derived alleles/traits and, in concert with QTL mapping results, test for congruence between our genotypic and phenotypic results. We found that natural selection overwhelmingly favours wild alleles and phenotypes. However, crop alleles in certain genomic regions can be favoured, and these changes often occurred in parallel across locations. We did not, however, consistently observe close agreement between our genotypic and phenotypic results. For example, when a trait evolved towards the wild phenotype, wild QTL alleles associated with that trait did not consistently increase in frequency. We discuss these results in the context of crop allele introgression into wild populations and implications for the management of GM crops.