Jacquelin DeFaveri

Global analysis of genes involved in freshwater adaptation in threespine sticklebacks (Gasterosteus aculeatus).

Examples of parallel evolution of phenotypic traits have been repeatedly demonstrated in threespine sticklebacks (Gasterosteus aculeatus) across their global distribution. Using these as a model, we performed a targeted genome scan—focusing on physiologically important genes potentially related to freshwater adaptation—to identify genetic signatures of parallel physiological evolution on a global scale. To this end, 50 microsatellite loci, including 26 loci within or close to (<6 kb) physiologically important genes, were screened in paired marine and freshwater populations from six locations across the Northern Hemisphere. Signatures of directional selection were detected in 24 loci, including 17 physiologically important genes, in at least one location. Although no loci showed consistent signatures of selection in all divergent population pairs, several outliers were common in multiple locations. In particular, seven physiologically important genes, as well as reference ectodysplasin gene (EDA), showed signatures of selection in three or more locations. Hence, although these results give some evidence for consistent parallel molecular evolution in response to freshwater colonization, they suggest that different evolutionary pathways may underlie physiological adaptation to freshwater habitats within the global distribution of the threespine stickleback.

Characterizing genic and nongenic molecular markers: comparison of microsatellites and SNPs.

The implications of transitioning to single nucleotide polymorphism (SNPs) from microsatellite markers (MSs) have been investigated in a number of population genetics studies, but the effect of genomic location on the amount of information each type of marker reveals has not been explored in detail. We developed novel SNP markers flanking 1 kb regions of 13 genic (within gene or <1 kb away from gene) and 13 nongenic (>10 kb from annotated gene) MSs in the threespine stickleback genome to obtain comparable data for both types of markers. We analysed patterns of genetic diversity and divergence on various geographic scales after converting the SNP loci within each genomic region into haplotypes. Marker type (SNP haplotype or MS) and location (genic or nongenic) significantly affected most estimates of population diversity and divergence. Between‐lineage divergence was significantly higher in SNP haplotypes (genic and nongenic), however, within‐lineage divergence was similar between marker types. Most divergence and diversity measures were uncorrelated between markers, except for population differentiation which was correlated between MSs and SNP haplotypes (both genic and nongenic). Broad‐scale population structure and assignment were similarly resolved by both marker types, however, only the MSs were able to delimit fine‐scale population structuring, particularly when genic and nongenic markers were combined. These results demonstrate that estimates of genetic variability and differentiation among populations can be strongly influenced by marker type, their genomic location in relation to genes and by the interaction of these two factors. This highlights the importance of having an awareness of the inherent strengths and limitations associated with different molecular tools to select the most appropriate methods for accurately addressing various ecological and evolutionary questions.

HETEROGENEOUS GENOMIC DIFFERENTIATION IN MARINE THREESPINE STICKLEBACKS: ADAPTATION ALONG AN ENVIRONMENTAL GRADIENT

Evolutionary divergence among populations occupying ecologically distinct environments can occur even in the face of on‐going gene flow. However, the genetic underpinnings, as well as the scale and magnitude at which this differentiation occurs in marine habitats are not well understood. We investigated the patterns and degree of genomic heterogeneity in threespine sticklebacks (Gasterosteus aculeatus) by assessing genetic variability in 20 nongenic and 20 genic (associated with genes important for freshwater adaptation) microsatellite loci in samples collected from 38 locations spanning the entire Baltic Sea coast to the North Sea boundary. Population divergence (FST ≈ 0.026) and structuring (five genetic clusters) was significantly more pronounced in the genic as compared to nongenic markers (FST ≈ 0.008; no genetic clusters). Patterns of divergence in the genic markers—45% of which were identified as outliers—correlated with local differences in salinity. Yet, a strong positive correlation between divergence in genic and nongenic markers, and their association with environmental factors suggests that adaptive divergence is reducing gene flow across the genome. Apart from providing a clear demonstration of heterogeneous genomic patterns of differentiation in a marine species, the results are indicative of adaptive population structuring across the relatively young Baltic Sea in spite of ample opportunities for gene flow.

Population genomic evidence for adaptive differentiation in Baltic Sea three-spined sticklebacks

BackgroundThe degree of genetic differentiation among populations experiencing high levels of gene flow is expected to be low for neutral genomic sites, but substantial divergence can occur in sites subject to directional selection. Studies of highly mobile marine fish populations provide an opportunity to investigate this kind of heterogeneous genomic differentiation, but most studies to this effect have focused on a relatively low number of genetic markers and/or few populations. Hence, the patterns and extent of genomic divergence in high-gene-flow marine fish populations remain poorly understood.ResultsWe here investigated genome-wide patterns of genetic variability and differentiation in ten marine populations of three-spined stickleback (Gasterosteus aculeatus) distributed across a steep salinity and temperature gradient in the Baltic Sea, by utilizing >30,000 single nucleotide polymorphisms obtained with a pooled RAD-seq approach. We found that genetic diversity and differentiation varied widely across the genome, and identified numerous fairly narrow genomic regions exhibiting signatures of both divergent and balancing selection. Evidence was uncovered for substantial genetic differentiation associated with both salinity and temperature gradients, and many candidate genes associated with local adaptation in the Baltic Sea were identified.ConclusionsThe patterns of genetic diversity and differentiation, as well as candidate genes associated with adaptation, in Baltic Sea sticklebacks were similar to those observed in earlier comparisons between marine and freshwater populations, suggesting that similar processes may be driving adaptation to brackish and freshwater environments. Taken together, our results provide strong evidence for heterogenic genomic divergence driven by local adaptation in the face of gene flow along an environmental gradient in the post-glacially formed Baltic Sea.

Evidence for adaptive phenotypic differentiation in Baltic Sea sticklebacks.

The evidence for adaptive phenotypic differentiation in mobile marine species remains scarce, partly due to the difficulty of obtaining quantitative genetic data to demonstrate the genetic basis of the observed phenotypic differentiation. Using a combination of phenotypic and molecular genetic approaches, we elucidated the relative roles of natural selection and genetic drift in explaining lateral plate number differentiation in threespine sticklebacks (Gasterosteus aculeatus) across the entire Baltic Sea basin (approximately 392 000 km2). We found that phenotypic differentiation (PST = 0.213) in plate number exceeded that in neutral markers (FST = 0.008), suggesting an adaptive basis for the observed differentiation. Because a close correspondence was found between plate phenotype and genotype at a quantitative trait loci (QTL; STN381) tightly linked to the gene (Ectodysplasin) underlying plate variation, the evidence for adaptive differentiation was confirmed by comparison of FST at the QTL (FSTQ = 0.089) with FST at neutral marker loci. Hence, the results provide a comprehensive demonstration of adaptive phenotypic differentiation in a high‐gene‐flow marine environment with direct, rather than inferred, verification for the genetic basis of this differentiation. In general, the results illustrate the utility of PST–FST–FSTQ comparisons in uncovering footprints of natural selection and evolution and add to the growing evidence for adaptive genetic differentiation in high‐gene‐flow marine environments, including that of the relatively young Baltic Sea.

Local adaptation to salinity in the three‐spined stickleback?

Different lines of evidence suggest that the occurrence and extent of local adaptation in high gene flow marine environments – even in mobile and long‐lived vertebrates with complex life cycles – may be more widespread than earlier thought. We conducted a common garden experiment to test for local adaptation to salinity in Baltic Sea sticklebacks (Gasterosteus aculeatus). Fish from three different native salinity regimes (high, mid and low) were subjected to three salinity treatments (high, mid and low) in a full‐factorial experimental design. Irrespective of their origin, fish subjected to low (and mid) salinity treatments exhibited higher juvenile survival, grew to largest sizes and were in better condition than fish subjected to the high salinity treatment. However, a significant interaction between native and treatment salinities – resulting mainly from the poor performance of fish native to low salinity in the high salinity treatment – provided clear cut evidence for adaptation to local variation in salinity. Additional support for this inference was provided by the fact that the results concur with an earlier demonstration of significant differentiation in a number of genes with osmoregulatory functions across the same populations and that the population‐specific responses to salinity treatments exceeded that to be expected by random genetic drift.

Genetic biodiversity in the Baltic Sea: species-specific patterns challenge management

Information on spatial and temporal patterns of genetic diversity is a prerequisite to understanding the demography of populations, and is fundamental to successful management and conservation of species. In the sea, it has been observed that oceanographic and other physical forces can constitute barriers to gene flow that may result in similar population genetic structures in different species. Such similarities among species would greatly simplify management of genetic biodiversity. Here, we tested for shared genetic patterns in a complex marine area, the Baltic Sea. We assessed spatial patterns of intraspecific genetic diversity and differentiation in seven ecologically important species of the Baltic ecosystem—Atlantic herring (Clupea harengus), northern pike (Esox lucius), European whitefish (Coregonus lavaretus), three-spined stickleback (Gasterosteus aculeatus), nine-spined stickleback (Pungitius pungitius), blue mussel (Mytilus spp.), and bladderwrack (Fucus vesiculosus). We used nuclear genetic data of putatively neutral microsatellite and SNP loci from samples collected from seven regions throughout the Baltic Sea, and reference samples from North Atlantic areas. Overall, patterns of genetic diversity and differentiation among sampling regions were unique for each species, although all six species with Atlantic samples indicated strong resistence to Atlantic-Baltic gene-flow. Major genetic barriers were not shared among species within the Baltic Sea; most species show genetic heterogeneity, but significant isolation by distance was only detected in pike and whitefish. These species-specific patterns of genetic structure preclude generalizations and emphasize the need to undertake genetic surveys for species separately, and to design management plans taking into consideration the specific structures of each species.

High degree of genetic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (FST = 0.033 or 0.011, respectively), whereas average FST for directionally selected markers was significantly higher (FST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.

Variation in Age and Size in Fennoscandian Three-Spined Sticklebacks (Gasterosteus aculeatus)

Average age and maximum life span of breeding adult three-spined sticklebacks (Gasterosteus aculeatus) were determined in eight Fennoscandian localities with the aid of skeletochronology. The average age varied from 1.8 to 3.6 years, and maximum life span from three to six years depending on the locality. On average, fish from marine populations were significantly older than those from freshwater populations, but variation within habitat types was large. We also found significant differences in mean body size among different habitat types and populations, but only the population differences remained significant after accounting for variation due to age effects. These results show that generation length and longevity in three-spined sticklebacks can vary significantly from one locality to another, and that population differences in mean body size cannot be explained as a simple consequence of differences in population age structure. We also describe a nanistic population from northern Finland exhibiting long life span and small body size.

Phylogeography of isolated freshwater three-spined stickleback Gasterosteus aculeatus populations in the Adriatic Sea basin.

Analyses of mitochondrial (mt) DNA and microsatellite variation were carried out to examine the relationships between 10 freshwater populations of three-spined sticklebacks Gasterosteus aculeatus along the eastern coast of the Adriatic Sea. Partial sequences of the mtDNA control region and cytochrome b gene, in addition to 15 microsatellite loci, were used to analyse populations from four isolated river catchments. Results uncovered an Adriatic lineage that was clearly divergent from the European lineage, and confirmed that the most divergent and ancient populations are located within the Adriatic lineage as compared with other European populations. Two northern Adriatic populations formed independent clades within the European mitochondrial lineage, suggesting different colonization histories of the different Adriatic populations. Nuclear marker analyses also indicated deep divergence between Adriatic and European populations, albeit with some discordance between the mtDNA phylogeny of the northern Adriatic populations, further highlighting the strong differentiation among the Adriatic populations. The southern populations within the Adriatic lineage were further organized into distinct clades corresponding to respective river catchments and sub-clades corresponding to river tributaries, reflecting a high degree of population structuring within a small geographic region, concurrent with suggestions of existence of several microrefugia within the Balkan Peninsula. The highly divergent clades and haplotypes unique to the southern Adriatic populations further suggest, in accordance with an earlier, more limited survey, that southern Adriatic populations represent an important reservoir for ancient genetic diversity of G. aculeatus.