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.

Genetic structure of West Greenland populations of lumpfish Cyclopterus lumpus.

In this study, 11 microsatellite markers were used to determine the structure of West Greenlandic lumpfish Cyclopterus lumpus populations across six spawning locations spanning >1500 km and compared with neighbouring populations in Canada and Iceland. To evaluate whether data allow for identification of origin of C. lumpus in Greenlandic waters, genetic assignment analysis was performed for 86 C. lumpus sampled on a feeding migration. Significant structuring with isolation by distance was observed in the West Greenland samples and two major subpopulations, north and south, were suggested. Based on FST values, closer relationships were observed between Greenland and Canada, than Greenland and Iceland. Surprisingly, the North Greenland population showed more similarities with Canadian samples, than did the geographically closer south-west Greenland population. Origin could be assigned for a high proportion of non-spawning fish and demonstrated a marked east-west spatial separation of fish of Greenlandic and Icelandic genotypes.

Signatures of natural selection between life cycle stages separated by metamorphosis in European eel.

BACKGROUND: Species showing complex life cycles provide excellent opportunities to study the genetic associations between life cycle stages, as selective pressures may differ before and after metamorphosis. The European eel presents a complex life cycle with two metamorphoses, a first metamorphosis from larvae into glass eels (juvenile stage) and a second metamorphosis into silver eels (adult stage). We tested the hypothesis that different genes and gene pathways will be under selection at different life stages when comparing the genetic associations between glass eels and silver eels. RESULTS: We used two sets of markers to test for selection: first, we genotyped individuals using a panel of 80 coding-gene single nucleotide polymorphisms (SNPs) developed in American eel; second, we investigated selection at the genome level using a total of 153,423 RAD-sequencing generated SNPs widely distributed across the genome. Using the RAD approach, outlier tests identified a total of 2413 (1.57%) potentially selected SNPs. Functional annotation analysis identified signal transduction pathways as the most over-represented group of genes, including MAPK/Erk signalling, calcium signalling and GnRH (gonadotropin-releasing hormone) signalling. Many of the over-represented pathways were related to growth, while others could result from the different conditions that eels inhabit during their life cycle. CONCLUSIONS: The observation of different genes and gene pathways under selection when comparing glass eels vs. silver eels supports the adaptive decoupling hypothesis for the benefits of metamorphosis. Partitioning the life cycle into discrete morphological phases may be overall beneficial since it allows the different life stages to respond independently to their unique selection pressures. This might translate into a more effective use of food and niche resources and/or performance of phase-specific tasks (e.g. feeding in the case of glass eels, migrating and reproducing in the case of silver eels).

Assessing patterns of hybridization between North Atlantic eels using diagnostic single-nucleotide polymorphisms.

The two North Atlantic eel species, the European eel (Anguilla anguilla) and the American eel (Anguilla rostrata), spawn in partial sympatry in the Sargasso Sea, providing ample opportunity to interbreed. In this study, we used a RAD (Restriction site Associated DNA) sequencing approach to identify species-specific diagnostic single-nucleotide polymorphisms (SNPs) and design a low-density array that combined with screening of a diagnostic mitochondrial DNA marker. Eels from Iceland (N=159) and from the neighboring Faroe Islands (N=29) were genotyped, along with 94 larvae (49 European and 45 American eel) collected in the Sargasso Sea. Our SNP survey showed that the majority of Icelandic eels are pure European eels but there is also an important contribution of individuals of admixed ancestry (10.7%). Although most of the hybrids were identified as F1 hybrids from European eel female × American eel male crosses, backcrosses were also detected, including a first-generation backcross (F1 hybrid × pure European eel) and three individuals identified as second-generation backcrosses originating from American eel × F1 hybrid backcrosses interbreeding with pure European eels. In comparison, no hybrids were observed in the Faroe Islands, the closest bodies of land to Iceland. It is possible that hybrids show an intermediate migratory behaviour between the two parental species that ultimately brings hybrid larvae to the shores of Iceland, situated roughly halfway between the Sargasso Sea and Europe. Only two hybrids were observed among Sargasso Sea larvae, both backcrosses, but no F1 hybrids, that points to temporal variation in the occurrence of hybridization.

Imprints from genetic drift and mutation imply relative divergence times across marine transition zones in a pan-European small pelagic fish (Sprattus sprattus).

Geographic distributions of most temperate marine fishes are affected by postglacial recolonisation events, which have left complex genetic imprints on populations of marine species. This study investigated population structure and demographic history of European sprat (Sprattus sprattus L.) by combining inference from both mtDNA and microsatellite genetic markers throughout the species' distribution. We compared effects from genetic drift and mutation for both genetic markers in shaping genetic differentiation across four transition zones. Microsatellite markers revealed significant isolation by distance and a complex population structure across the species' distribution (overall θ(ST)=0.038, P<0.01). Across transition zones markers indicated larger effects of genetic drift over mutations in the northern distribution of sprat contrasting a stronger relative impact of mutation in the species' southern distribution in the Mediterranean region. These results were interpreted to reflect more recent divergence times between northern populations in accordance with previous findings. This study demonstrates the usefulness of comparing inference from different markers and estimators of divergence for phylogeographic and population genetic studies in species with weak genetic structure, as is the case in many marine species.

Application of SNPs for population genetics of nonmodel organisms: new opportunities and challenges.

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.

An assessment of the spatial scale of local adaptation in brown trout (Salmo trutta L.): footprints of selection at microsatellite DNA loci.

Local adaptation is considered a paradigm in studies of salmonid fish populations. Yet, little is known about the geographical scale of local adaptation. Is adaptive divergence primarily evident at the scale of regions or individual populations? Also, many salmonid populations are subject to spawning intrusion by farmed conspecifics that experience selection regimes fundamentally different from wild populations. This prompts the question if adaptive differences between wild populations and hatchery strains are more pronounced than between different wild populations? We addressed these issues by analyzing variation at 74 microsatellite loci (including anonymous and expressed sequence tag- and quantitative trait locus-linked markers) in 15 anadromous wild brown trout (Salmo trutta L.) populations, representing five geographical regions, along with two lake populations and two hatchery strains used for stocking some of the populations. F(ST)-based outlier tests revealed more outlier loci between different geographical regions separated by 522 ± 228 km (mean ± s.d.) than between populations within regions separated by 117 ± 79 km (mean ± s.d.). A significant association between geographical distance and number of outliers between regions was evident. There was no evidence for more outliers in comparisons involving hatchery trout, but the loci under putative selection generally were not the same as those found to be outliers between wild populations. Our study supports the notion of local adaption being increasingly important at the scale of regions as compared with individual populations, and suggests that loci involved in adaptation to captive environments are not necessarily the same as those involved in adaptive divergence among wild populations.

Detecting population structure in a high gene-flow species, Atlantic herring (Clupea harengus): direct, simultaneous evaluation of neutral vs putatively selected loci.

In many marine fish species, genetic population structure is typically weak because populations are large, evolutionarily young and have a high potential for gene flow. We tested whether genetic markers influenced by natural selection are more efficient than the presumed neutral genetic markers to detect population structure in Atlantic herring (Clupea harengus), a migratory pelagic species with large effective population sizes. We compared the spatial and temporal patterns of divergence and statistical power of three traditional genetic marker types, microsatellites, allozymes and mitochondrial DNA, with one microsatellite locus, Cpa112, previously shown to be influenced by divergent selection associated with salinity, and one locus located in the major histocompatibility complex class IIA (MHC-IIA) gene, using the same individuals across analyses. Samples were collected in 2002 and 2003 at two locations in the North Sea, one location in the Skagerrak and one location in the low-saline Baltic Sea. Levels of divergence for putatively neutral markers were generally low, with the exception of single outlier locus/sample combinations; microsatellites were the most statistically powerful markers under neutral expectations. We found no evidence of selection acting on the MHC locus. Cpa112, however, was highly divergent in the Baltic samples. Simulations addressing the statistical power for detecting population divergence showed that when using Cpa112 alone, compared with using eight presumed neutral microsatellite loci, sample sizes could be reduced by up to a tenth while still retaining high statistical power. Our results show that the loci influenced by selection can serve as powerful markers for detecting population structure in high gene-flow marine fish species.

Characterization of nine polymorphic microsatellite markers in sprat (Sprattus sprattus L.).

Nine polymorphic microsatellites were isolated from sprat (Sprattus sprattus) using a microsatellite enrichment protocol and selective hybridization with a biotinylated (AC)(12) probe. The loci showed different variation patterns in a Baltic Sea population (44 individuals) with mean number of alleles at 12.7 and mean observed heterozygosity at 0.78. These microsatellite loci are expected to be used for taxonomic considerations in sprat, stock differentiation and population genetic analysis.

Male reproductive competition in spawning aggregations of cod (Gadus morhua, L.).

Reproductive competition may lead to a large skew in reproductive success among individuals. Very few studies have analysed the paternity contribution of individual males in spawning aggregations of fish species with huge census population sizes. We quantified the variance in male reproductive success in spawning aggregations of cod under experimental conditions over an entire spawning season. Male reproductive success was estimated by microsatellite-based parentage analysis of offspring produced in six separate groups of spawning cod. In total, 1340 offspring and 102 spawnings distributed across a spawning season were analysed. Our results show that multiple males contributed sperm to most spawnings but that paternity frequencies were highly skewed among males, with larger males on average siring higher proportions of offspring. It was further indicated that male reproductive success was dependent on the magnitude of the size difference between a female and a male. We discuss our results in relation to the cod mating system. Finally, we suggest that the highly skewed distribution of paternity success observed in cod may be a factor contributing to the low effective population size/census population size ratios observed in many marine organisms.