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.
Abstract
Environmental DNA (eDNA) has attracted interest in relation to fisheries, with its possibilities for species identification and promises for species quantification. In the context of ...fisheries catches, eDNA can be most useful for the estimation of bycatch proportions. The assessment of species mixtures in large catches (>1000 t) is challenging, especially when morphologically similar species are to be differentiated. We used an experimental set-up to simulate industrial sprat fishery catches, and tested two types of water (blood water and discharge water) derived from this simulated fishery for their suitability in reliable species quantification. We analysed nine mixtures of sprat and herring—the main bycatch species. Species-specific quantitative PCR was used for species identification and quantification. Species-to-species weight fractions and eDNA fractions in mixtures showed a strong correlation. Accounting for size-based differences in DNA abundance (allometrically scaled weight) reduced the estimated standard error on weight fraction prediction from 0.064 to 0.054 in blood water, and from 0.080 to 0.075 in discharge water when comparing the weight-based model with the allometrically scaled weight model, respectively. Accounting for allometric scalling in genetic analyses of fisheries process water can serve as a more precise method for the assessment of bycatch, thus in a wider sense improve the quality of fisheries-dependent data.
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.
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.
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.
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.
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.
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.
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).
Estimates of Atlantic cod Gadus morhua reproductive success, determined using experimental spawning groups and genetic paternity assignment of offspring, showed that within‐group variance in male ...size correlated positively with the degree of male mating skew, predicting a decrease in male reproductive skew with decreasing size variation among males under natural conditions.
Determining spatio-temporal distributions of fish populations is of interest to marine ecology, in general, and to fisheries science in particular. Genetic mixed-stock analysis is routinely applied ...in several anadromous fishes for determining migratory routes and timing but has rarely been used for marine fishes, for which population differentiation is commonly weak and the method presumably less powerful. We used microsatellite information for Northeast Atlantic herringClupea harengusL. populations and mixed stocks to address 2 questions. We used simulated mixture samples and 3 different statistical approaches to determine whether mixed stock composition could be determined with accuracy. Simulations showed that the applied approaches and mixture samples of 100 individuals enabled detailed composition analyses on a regional level, with resolution for tracing the ecologically dominant Rügen (Greifswalder Bodden) herring population. We then estimated spatio-temporal variation in herring migratory behaviour in the Skagerrak from 17 mixed samples collected over 2 seasons and 2 yr, and identified hitherto undescribed differences in distributions among populations that feed and winter in the area.
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.