Deciphering genetic structure and inferring connectivity in marine species have been challenging due to weak genetic differentiation and limited resolution offered by traditional genotypic methods. The main goal of this study was to assess how a population genomics framework could help delineate the genetic structure of the American lobster (Homarus americanus) throughout much of the species' range and increase the assignment success of individuals to their location of origin. We genotyped 10 156 filtered SNPs using RAD sequencing to delineate genetic structure and perform population assignment for 586 American lobsters collected in 17 locations distributed across a large portion of the species' natural distribution range. Our results revealed the existence of a hierarchical genetic structure, first separating lobsters from the northern and southern part of the range (FCT  = 0.0011; P-value = 0.0002) and then revealing a total of 11 genetically distinguishable populations (mean FST  = 0.00185; CI: 0.0007-0.0021, P-value < 0.0002), providing strong evidence for weak, albeit fine-scale population structuring within each region. A resampling procedure showed that assignment success was highest with a subset of 3000 SNPs having the highest FST . Applying Anderson's (Molecular Ecology Resources, 2010, 10, 701) method to avoid 'high-grading bias', 94.2% and 80.8% of individuals were correctly assigned to their region and location of origin, respectively. Lastly, we showed that assignment success was positively associated with sample size. These results demonstrate that using a large number of SNPs improves fine-scale population structure delineation and population assignment success in a context of weak genetic structure. We discuss the implications of these findings for the conservation and management of highly connected marine species, particularly regarding the geographic scale of demographic independence.

Understanding the evolutionary mechanisms generating parallel genomic divergence patterns among replicate ecotype pairs remains an important challenge in speciation research. We investigated the genomic divergence between the anadromous parasitic river lamprey (Lampetra fluviatilis) and the freshwater-resident nonparasitic brook lamprey (Lampetra planeri) in nine population pairs displaying variable levels of geographic connectivity. We genotyped 338 individuals with RAD sequencing and inferred the demographic divergence history of each population pair using a diffusion approximation method. Divergence patterns in geographically connected population pairs were better explained by introgression after secondary contact, whereas disconnected population pairs have retained a signal of ancient migration. In all ecotype pairs, models accounting for differential introgression among loci outperformed homogeneous migration models. Generating neutral predictions from the inferred divergence scenarios to detect highly differentiated markers identified greater proportions of outliers in disconnected population pairs than in connected pairs. However, increased similarity in the most divergent genomic regions was found among connected ecotype pairs, indicating that gene flow was instrumental in generating parallelism at the molecular level. These results suggest that heterogeneous genomic differentiation and parallelism among replicate ecotype pairs have partly emerged through restricted introgression in genomic islands.

Salmonids rank among the most socioeconomically valuable fishes and the most targeted species by stocking with hatchery-reared individuals. Here, we used molecular parentage analysis to assess the reproductive success of wild-and hatchery-born Atlantic salmon over three consecutive years in a small river in Québec. Yearly restocking in this river follows a single generation of captive breeding. Among the adults returning to the river to spawn, between 11% and 41% each year were born in hatchery. Their relative reproductive success (RRS) was nearly half that of wild-born fish (0.55). RRS varied with life stage, being 0.71 for fish released at the fry stage and 0.42 for fish released as smolt. The lower reproductive success of salmon released as smolt was partly mediated by the modification of the proportion of single-sea-winter/multi-sea-winter fish. Overall, our results suggest that modifications in survival and growth rates alter the life-history strategies of these fish at the cost of their reproductive success. Our results underline the potential fitness decrease, warn on long-term evolutionary consequences for the population of repeated stocking and support the adoption of more natural rearing conditions for captive juveniles and their release at a younger stage, such as unfed fry.

Stocking represents the most important management tool worldwide to increase and sustain commercial and recreational fisheries in a context of overexploitation. Genetic impacts of this practice have been investigated in many studies, which examined population and individual admixture, but few have investigated determinants of these processes. Here, we addressed these questions from the genotyping at 19 microsatellite loci of 3341 adult lake trout (Salvelinus namaycush) from 72 unstocked and stocked lakes. Results showed an increase in genetic diversity and a twofold decrease in the extent of genetic differentiation among stocked populations when compared to unstocked. Stocked populations were characterized by significant admixture at both population and individual levels. Moreover, levels of admixture in stocked populations were strongly correlated with stocking intensity and a threshold value of total homogenization between source and stocked populations was identified. Our results also suggest that under certain scenarios, the genetic impacts of stocking could be of short duration. Overall, our study emphasizes the important alteration of the genetic integrity of stocked populations and the need to better understand determinants of admixture to optimize stocking strategies and to conserve the genetic integrity of wild populations.

Disentangling the effects of natural environmental features and anthropogenic factors on the genetic structure of endangered populations is an important challenge for conservation biology. Here, we investigated the combined influences of major environmental features and stocking with non-native fish on the genetic structure and local adaptation of Atlantic salmon (Salmo salar) populations. We used 17 microsatellite loci to genotype 975 individuals originating from 34 French rivers. Bayesian analyses revealed a hierarchical genetic structure into five geographically distinct clusters. Coastal distance, geological substrate and river length were strong predictors of population structure. Gene flow was higher among rivers with similar geologies, suggesting local adaptation to geological substrate. The effect of river length was mainly owing to one highly differentiated population that has the farthest spawning grounds off the river mouth (up to 900km) and the largest fish, suggesting local adaptation to river length. We detected high levels of admixture in stocked populations but also in neighbouring ones, implying large-scale impacts of stocking through dispersal of non-native individuals. However, we found relatively few admixed individuals suggesting a lower fitness of stocked fish and/or some reproductive isolation between wild and stocked individuals. When excluding stocked populations, genetic structure increased as did its correlation with environmental factors. This study overall indicates that geological substrate and river length are major environmental factors influencing gene flow and potential local adaptation among Atlantic salmon populations but that stocking with non-native individuals may ultimately disrupt these natural patterns of gene flow among locally adapted populations.

Little is known about the genetic basis differentiating resident and anadromous forms found in many salmonid species. Using a medium-density SNP array, we documented genomic diversity and divergence at 2336 genetically mapped loci among three pairs of North American anadromous and freshwater Atlantic salmon populations. Our results show that across the genome, freshwater populations have lower diversity and a smaller proportion of private polymorphism relative to anadromous populations. Moreover, differentiation was more pronounced among freshwater than among anadromous populations at multiple spatial scales, suggesting a large effect of genetic drift in these isolated freshwater populations. Using nonhierarchical and hierarchical genome scans, we identified hundreds of markers spread across the genome that are potentially under divergent selection between anadromous and freshwater populations, but few outlier loci were repeatedly found in all three freshwater–anadromous comparisons. Similarly, a sliding window analysis revealed numerous regions of high divergence that were nonparallel among the three comparisons. These last results show little evidence for the parallel evolution of alleles selected for in freshwater populations, but suggest nonparallel adaptive divergence at many loci of small effects distributed through the genome. Overall, this study emphasizes the important role of genetic drift in driving genome-wide reduction in diversity and divergence in freshwater Atlantic salmon populations and suggests a complex multigenic basis of adaptation to resident and anadromous strategies with little parallelism.

While the stocking of captive-bred fish has been occurring for decades and has had substantial immediate genetic and evolutionary impacts on wild populations, its long-term consequences have only been weakly investigated. Here, we conducted a spatiotemporal analysis of 1428 Atlantic salmon sampled from 1965 to 2006 in 25 populations throughout France to investigate the influence of stocking on the neutral genetic structure in wild Atlantic salmon (Salmo salar) populations. On the basis of the analysis of 11 microsatellite loci, we found that the overall genetic structure among populations dramatically decreased over the period studied. Admixture rates among populations were highly variable, ranging from a nearly undetectable contribution from donor stocks to total replacement of the native gene pool, suggesting extremely variable impacts of stocking. Depending on population, admixture rates either increased, remained stable, or decreased in samples collected between 1998 and 2006 compared to samples from 1965 to 1987, suggesting either rising, long-lasting or short-term impacts of stocking. We discuss the potential mechanisms contributing to this variability, including the reduced fitness of stocked fish and persistence of wild locally adapted individuals.

Genome scans represent powerful approaches to investigate the action of natural selection on the genetic variation of natural populations and to better understand local adaptation. This is very useful, for example, in the field of conservation biology and evolutionary biology. Thanks to Next Generation Sequencing, genomic resources are growing exponentially, improving genome scan analyses in non-model species. Thousands of SNPs called using Reduced Representation Sequencing are increasingly used in genome scans. Besides, genome sequences are also becoming increasingly available, allowing better processing of short-read data, offering physical localization of variants, and improving haplotype reconstruction and data imputation. Ultimately, genome sequences are also becoming the raw material for selection inferences. Here, we discuss how the increasing availability of such genomic resources, notably genome sequences, influences the detection of signals of selection. Mainly, increasing data density and having the information of physical linkage data expand genome scans by (i) improving the overall quality of the data, (ii) helping the reconstruction of demographic history for the population studied to decrease false-positive rates and (iii) improving the statistical power of methods to detect the signal of selection. Of particular importance, the availability of a high-quality reference genome can improve the detection of the signal of selection by (i) allowing matching the potential candidate loci to linked coding regions under selection, (ii) rapidly moving the investigation to the gene and function and (iii) ensuring that the highly variable regions of the genomes that include functional genes are also investigated. For all those reasons, using reference genomes in genome scan analyses is highly recommended.