Stocking impact and temporal stability of genetic composition in a brackish northern pike population (Esox lucius L.), assessed using microsatellite DNA analysis of historical and contemporary samples

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

Appraising the Genetic Makeup of an Allochthonous Southern Pike Population: An Opportunity to Predict the Evolution of Introgressive Hybridization in Isolated Populations?

Biological invasions are a major threat to the conservation of biodiversity, as invasive species affect native biota through competition, predation, pathogen introduction, habitat alteration, and hybridisation. The present study focuses on a southern pike population, Esox cisalpinus (Teleostei: Esocidae), that has been introduced outside the species' native range. Using microsatellite markers, this study's objective was to gather baseline genetic information and assess the presence of hybrids between this species and E. lucius in the introduced population. The resulting estimates of genetic diversity and effective population size are comparable to those observed in the species' native range. Although different methods yield contrasting and uncertain evidence regarding introgressive hybridization, the presence of late-generation hybrids cannot be completely ruled out. Large numbers of breeders as well as multiple introductions of genetically divergent cohorts and introgressive hybridisation may explain the high genetic diversity of this recently introduced southern pike population. The present study issues a warning that the conservation of southern pike' introgressive hybridisation between northern and southern pike might be underestimated. The genetic information gathered herein may unravel the origin, number of introduction events, and evolutionary trajectory of the introduced population. This information may help us understand the evolution of introgressive hybridisation in the southern pike's native areas.

How well do genetic markers inform about responses to intraspecific admixture? A comparative analysis of microsatellites and RADseq

Background

Fitness consequences of intraspecific genetic admixture can vary from positive to negative depending on the genetic composition of the populations and environmental conditions. Because admixture has potential to influence the success of management and conservation efforts, genetic similarity has been suggested to be used as a proxy to predict the outcome. Studies utilizing microsatellites (a neutral marker) to investigate associations between genetic distance and admixture effects show conflicting results. Marker types that yield information on genome-wide and/or adaptive variation might be more useful for predicting responses to inter-population hybridization. In this study we utilized published data for three populations of pike (Esox lucius) to investigate associations between offspring performance (hatching success) and parental genetic similarity in experimentally purebred and admixed families, based on neutral (microsatellites), genome-wide neutral (RADseq SNPs), and adaptive (SNPs under selection) markers.

Results

Estimated similarity varied among the markers, likely reflecting differences in their inherent properties, but was consistently higher in purebred than admixed families. A significant interaction between marker type and admixture treatment reflected that neutral SNPs yielded higher estimates than adaptive SNPs for admixed families whereas no difference was found for purebred families, which indicates that neutral similarity was not reflective of adaptive similarity. When all samples were pooled, no association between similarity and performance was found for any marker. For microsatellites, similarity was positively correlated with hatching success in purebred families, whereas no association was found in admixed families; however, the direction of the effect differed between the population combinations.

Conclusions

The results strengthen the notion that, as of today, there is no proxy that can reliably predicted the outcome of admixture. This emphasizes the need of further studies to advance knowledge that can shed light on how to safeguard against negative consequences of admixture, and thereby inform management and promote conservation of biological diversity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00974-3.

Catch-related and genetic outcome of adult northern pike Esox lucius stocking in a large river system

Genetic introgression from stocked adult northern pike Esox lucius to a wild self-recruiting population was detected in a large river system and some stocked E. lucius survived up to two spawning seasons and dispersed over several kilometres in the river. Moreover, the catch rate of stocked E. lucius by anglers was low (9.6%), hence suggesting that the efficiency of stocking activity is questionable.

Sex-specific effects of outbreeding on offspring quality in pike (Esox lucius)

Intraspecific genetic admixture occurs when previously separated populations within a species start interbreeding, and it can have either positive, negative, or neutral effects on reproductive performance. As there currently is no reliable predictor for the outcome of admixture, an increased knowledge about admixture effects in different species and populations is important to increase the understanding about what determines the response to admixture. We tested for effects of admixture on F1 offspring quality in three subpopulations of pike (Esox lucius). Gametes were collected in the field, and eggs from each female were experimentally fertilized with milt from a male from each population (one "pure" and two "admixed" treatments). Three offspring quality measures (hatching success, fry survival, and fry length) were determined and compared between (a) pure and admixed population combinations and (b) the sex-specific treatments within each admixed population combination (based on the origin of the male and female, respectively). The results suggested that although there were no overall effects of admixture on offspring quality, the consequences for a given population combination could be sex-specific and thus differ depending on which of the parents originated from one or the other population. All offspring quality traits were influenced by both maternal ID and paternal ID. Sex- and individual-specific effects can have implications for dispersal behavior and gene flow between natural populations, and are important to consider in conservation efforts.

Genomic population structure of freshwater-resident and anadromous ide (Leuciscus idus) in north-western Europe

Climate change experts largely agree that future climate change and associated rises in oceanic water levels over the upcoming decades, will affect marine salinity levels. The subsequent effects on fish communities in estuarine ecosystems however, are less clear. One species that is likely to become increasingly affected by changes in salinity is the ide (Leuciscus idus). The ide is a stenohaline freshwater fish that primarily inhabits rivers, with frequent anadromous behavior when sea salinity does not exceed 15%. Unlike most other anadromous Baltic Sea fish species, the ide has yet to be subjected to large-scale stocking programs, and thus provides an excellent opportunity for studying the natural population structure across the current salinity gradient in the Danish Belts. To explore this, we used Genotyping-by-Sequencing to determine genomic population structure of both freshwater resident and anadromous ide populations in the western Baltic Sea region, and relate the results to the current salinity gradient and the demographic history of ide in the region. The sample sites separate into four clusters, with all anadromous populations in one cluster and the freshwater resident populations in the remaining three. Results demonstrate high level of differentiation between sites hosting freshwater resident populations, but little differentiation among anadromous populations. Thus ide exhibit the genomic population structure of both a typical freshwater species, and a typical anadromous species. In addition to providing a first insight into the population structure of north-western European ide, our data also (1) provide indications of a single illegal introduction by man; (2) suggest limited genetic effects of heavy pollution in the past; and (3) indicate possible historical anadromous behavior in a now isolated freshwater population.

Population differentiation of zander (Sander lucioperca) across native and newly colonized ranges suggests increasing admixture in the course of an invasion

In addition to ecological factors, evolutionary processes can determine the invasion success of a species. In particular, genetic admixture has the potential to induce rapid evolutionary change, which can result from natural or human-assisted secondary contact between differentiated populations. We studied the recent range expansion of zander in Germany focusing on the interplay between invasion and genetic admixture. Historically, the rivers Elbe and Danube harboured the most north-western source populations from which a north-westward range expansion occurred. This was initiated by introducing zander outside its native range into rivers and lakes, and was fostered by migration through artificial canals and stocking from various sources. We analysed zander populations of the native and invaded ranges using nuclear and mitochondrial genetic markers. Three genetic lineages were identified, which were traced to ancestral ranges. Increased genetic diversity and admixture in the invaded region highlighted asymmetric gene flow towards this area. We suppose that the adaptive potential of the invading populations was promoted by genetic admixture, whereas competitive exclusion in the native areas provided a buffer against introgression by novel genotypes. These explanations would be in line with evidence that hybridization can drive evolutionary change under conditions when new niches can be exploited.

Resurrecting an extinct salmon evolutionarily significant unit: archived scales, historical DNA and implications for restoration

Archival scales from 603 sockeye salmon (Oncorhynchus nerka), sampled from May to July 1924 in the lower Columbia River, were analysed for genetic variability at 12 microsatellite loci and compared to 17 present-day O. nerka populations-exhibiting either anadromous (sockeye salmon) or nonanadromous (kokanee) life histories-from throughout the Columbia River Basin, including areas upstream of impassable dams built subsequent to 1924. Statistical analyses identified four major genetic assemblages of sockeye salmon in the 1924 samples. Two of these putative historical groupings were found to be genetically similar to extant evolutionarily significant units (ESUs) in the Okanogan and Wenatchee Rivers (pairwise F(ST)  = 0.004 and 0.002, respectively), and assignment tests were able to allocate 77% of the fish in these two historical groupings to the contemporary Okanogan River and Lake Wenatchee ESUs. A third historical genetic grouping was most closely aligned with contemporary sockeye salmon in Redfish Lake, Idaho, although the association was less robust (pairwise F(ST)  = 0.060). However, a fourth genetic grouping did not appear to be related to any contemporary sockeye salmon or kokanee population, assigned poorly to the O. nerka baseline, and had distinctive early return migration timing, suggesting that this group represents a historical ESU originating in headwater lakes in British Columbia that was probably extirpated sometime after 1924. The lack of a contemporary O. nerka population possessing the genetic legacy of this extinct ESU indicates that efforts to reestablish early-migrating sockeye salmon to the headwater lakes region of the Columbia River will be difficult.

Molecular and phenotypic evidence of a new species of genus Esox (Esocidae, Esociformes, Actinopterygii): the southern pike, Esox flaviae

We address the taxonomic position of the southern European individuals of pike, performing a series of tests and comparisons from morphology, DNA taxonomy and population genetics parameters, in order to support the hypothesis that two species of pike, and not only one, exist in Europe. A strong relationship emerged between a northern genotype supported by COI, Cytb, AFLP and specific fragments, and a phenotype with round spot skin colour pattern and a large number of scales in the lateral line, clearly separated from a southern genotype with other skin colour pattern and a low number of scales in the lateral line. DNA taxonomy, based on a coalescent approach (GMYC) from phylogenetic reconstructions on COI and Cytb together with AFLP admixture analysis, supported the existence of two independently evolving entities. Such differences are not simply due to geographic distances, as northern European samples are more similar to Canadian and Chinese samples than the southern Europe ones. Thus, given that the differences between the two groups of European pike are significant at the phenotypic, genotypic and geographical levels, we propose the identification of two pike species: the already known northern pike (Esox lucius) and the southern pike (E. flaviae n.sp.). The correct identification of these two lineages as independent species should give rise to a ban on the introduction of northern pikes in southern Europe for recreational fishing, due to potential problems of hybridisation.

Comparative estimation of effective population sizes and temporal gene flow in two contrasting population systems

Estimation of effective population sizes (N(e)) and temporal gene flow (N(e)m, m) has many implications for understanding population structure in evolutionary and conservation biology. However, comparative studies that gauge the relative performance of N(e), N(e)m or m methods are few. Using temporal genetic data from two salmonid fish population systems with disparate population structure, we (i) evaluated the congruence in estimates and precision of long- and short-term N(e), N(e)m and m from six methods; (ii) explored the effects of metapopulation structure on N(e) estimation in one system with spatiotemporally linked subpopulations, using three approaches; and (iii) determined to what degree interpopulation gene flow was asymmetric over time. We found that long-term N(e) estimates exceeded short-term N(e) within populations by 2-10 times; the two were correlated in the system with temporally stable structure (Atlantic salmon, Salmo salar) but not in the highly dynamic system (brown trout, Salmo trutta). Four temporal methods yielded short-term N(e) estimates within populations that were strongly correlated, and these were higher but more variable within salmon populations than within trout populations. In trout populations, however, these short-term N(e) estimates were always lower when assuming gene flow than when assuming no gene flow. Linkage disequilibrium data generally yielded short-term N(e) estimates of the same magnitude as temporal methods in both systems, but the two were uncorrelated. Correlations between long- and short-term geneflow estimates were inconsistent between methods, and their relative size varied up to eightfold within systems. While asymmetries in gene flow were common in both systems (58-63% of population-pair comparisons), they were only temporally stable in direction within certain salmon population pairs, suggesting that gene flow between particular populations is often intermittent and/or variable. Exploratory metapopulation N(e) analyses in trout demonstrated both the importance of spatial scale in estimating N(e) and the role of gene flow in maintaining genetic variability within subpopulations. Collectively, our results illustrate the utility of comparatively applying N(e), N(e)m and m to (i) tease apart processes implicated in population structure, (ii) assess the degree of continuity in patterns of connectivity between population pairs and (iii) gauge the relative performance of different approaches, such as the influence of population subdivision and gene flow on N(e) estimation. They further reiterate the importance of temporal sampling replication in population genetics, the value of interpreting N(e)or m in light of species biology, and the need to address long-standing assumptions of current N(e), N(e)m or m models more explicitly in future research.

When invaders meet locally adapted types: rapid moulding of hybrid zones between sculpins (Cottus, Pisces) in the Rhine system

Hybrid zones are commonly studied to dissect the processes that drive divergence among lineages, which have incomplete barriers of reproduction. Most hybrid zones have existed for an extended time making inferences on the initial mode of formation difficult. It is a priori unclear how fast a nascent hybrid zone would form as a response to endogenous and exogenous factors. We have studied several hybrid zones between two lineages of sculpins (Cottus spp.), which emerged due to a recent range expansion of one of the lineages along the river Rhine in the early 1990s. Applying a dense sampling across two contact areas and using a highly informative set of 45 microsatellite markers we found pronounced genetic structure. Steep genetic clines suggest that strong selective forces have shaped the respective hybrid zones from the beginning. We find that the zones are coupled to ecological transitions from small streams to larger rivers. The width of these zones is much smaller than estimates of annual individual dispersal distances, as estimated outside of the hybrid zones. The pattern is apparently not strongly affected by pre- or postzygotic reproductive isolation because numerous backcross hybrids occur within the zones. This suggests that strong natural selection acts against immigrant genotypes. The study exemplifies how local adaptation can play a key role in preventing admixture in dependence of the ecological context.