Fine-scale environmentally associated spatial structure of lumpfish (Cyclopterus lumpus) across the Northwest Atlantic.

Lumpfish, Cyclopterus lumpus, have historically been harvested throughout Atlantic Canada and are increasingly in demand as a solution to controlling sea lice in Atlantic salmon farms-a process which involves both the domestication and the transfer of lumpfish between geographic regions. At present, little is known regarding population structure and diversity of wild lumpfish in Atlantic Canada, limiting attempts to assess the potential impacts of escaped lumpfish individuals from salmon pens on currently at-risk wild populations. Here, we characterize the spatial population structure and genomic-environmental associations of wild populations of lumpfish throughout the Northwest Atlantic using both 70K SNP array data and whole-genome re-sequencing data (WGS). At broad spatial scales, our results reveal a large environmentally associated genetic break between the southern populations (Gulf of Maine and Bay of Fundy) and northern populations (Newfoundland and the Gulf of St. Lawrence), linked to variation in ocean temperature and ice cover. At finer spatial scales, evidence of population structure was also evident in a distinct coastal group in Newfoundland and significant isolation by distance across the northern region. Both evidence of consistent environmental associations and elevated genome-wide variation in F ST values among these three regional groups supports their biological relevance. This study represents the first extensive description of population structure of lumpfish in Atlantic Canada, revealing evidence of broad and fine geographic scale environmentally associated genomic diversity. Our results will facilitate the commercial use of lumpfish as a cleaner fish in Atlantic salmon aquaculture, the identification of lumpfish escapees, and the delineation of conservation units of this at-risk species throughout Atlantic Canada.

Whole-genome sequencing reveals fine-scale environment-associated divergence near the range limits of a temperate reef fish.

Environmental variation is increasingly recognized as an important driver of diversity in marine species despite the lack of physical barriers to dispersal and the presence of pelagic stages in many taxa. A robust understanding of the genomic and ecological processes involved in structuring populations is lacking for most marine species, often hindering management and conservation action. Cunner (Tautogolabrus adspersus) is a temperate reef fish with both pelagic early life-history stages and strong site-associated homing as adults; the species is also of interest for use as a cleaner fish in salmonid aquaculture in Atlantic Canada. We aimed to characterize genomic and geographic differentiation of cunner in the Northwest Atlantic. To achieve this, a chromosome-level genome assembly for cunner was produced and used to characterize spatial population structure throughout Atlantic Canada using whole-genome sequencing. The genome assembly spanned 0.72 Gbp and 24 chromosomes; whole-genome sequencing of 803 individuals from 20 locations from Newfoundland to New Jersey identified approximately 11 million genetic variants. Principal component analysis revealed four regional Atlantic Canadian groups. Pairwise FST and selection scans revealed signals of differentiation and selection at discrete genomic regions, including adjacent peaks on chromosome 10 across multiple pairwise comparisons (i.e. FST 0.5-0.75). Redundancy analysis suggested association of environmental variables related to benthic temperature and oxygen range with genomic structure. Results suggest regional scale diversity in this temperate reef fish and can directly inform the collection and translocation of cunner for aquaculture applications and the conservation of wild populations throughout the Northwest Atlantic.

Genomic and machine learning-based screening of aquaculture-associated introgression into at-risk wild North American Atlantic salmon (Salmo salar) populations.

The negative genetic impacts of gene flow from domestic to wild populations can be dependent on the degree of domestication and exacerbated by the magnitude of pre-existing genetic differences between wild populations and the domestication source. Recent evidence of European ancestry within North American aquaculture Atlantic salmon (Salmo salar) has elevated the potential impact of escaped farmed salmon on often at-risk wild North American salmon populations. Here, we compare the ability of single nucleotide polymorphism (SNP) and microsatellite (SSR) marker panels of different sizes (7-SSR, 100-SSR and 220K-SNP) to detect introgression of European genetic information into North American wild and aquaculture populations. Linear regression comparing admixture predictions for a set of individuals common to the three datasets showed that the 100-SSR panel and 7-SSR panels replicated the full 220K-SNP-based admixture estimates with low accuracy (r2 of .64 and .49, respectively). Additional tests explored the effects of individual sample size and marker number, which revealed that ~300 randomly selected SNPs could replicate the 220K-SNP admixture predictions with greater than 95% fidelity. We designed a custom SNP panel (301-SNP) for European admixture detection in future monitoring work and then developed and tested a python package, salmoneuadmix (https://github.com/CNuge/SalmonEuAdmix), which uses a deep neural network to make de novo estimates of individuals' European admixture proportion without the need to conduct complete admixture analysis utilizing baseline samples. The results demonstrate the mobilization of targeted SNP panels and machine learning in support of at-risk species conservation and management.

Resolving fine-scale population structure and fishery exploitation using sequenced microsatellites in a northern fish.

The resiliency of populations and species to environmental change is dependent on the maintenance of genetic diversity, and as such, quantifying diversity is central to combating ongoing widespread reductions in biodiversity. With the advent of next-generation sequencing, several methods now exist for resolving fine-scale population structure, but the comparative performance of these methods for genetic assignment has rarely been tested. Here, we evaluate the performance of sequenced microsatellites and a single nucleotide polymorphism (SNP) array to resolve fine-scale population structure in a critically important salmonid in north eastern Canada, Arctic Charr (Salvelinus alpinus). We also assess the utility of sequenced microsatellites for fisheries applications by quantifying the spatial scales of movement and exploitation through genetic assignment of fishery samples to rivers of origin and comparing these results with a 29-year tagging dataset. Self-assignment and simulation-based analyses of 111 genome-wide microsatellite loci and 500 informative SNPs from 28 populations of Arctic Charr in north-eastern Canada identified largely river-specific genetic structure. Despite large differences (~4X) in the number of loci surveyed between panels, mean self-assignment accuracy was similar with the microsatellite loci and the SNP panel (>90%). Subsequent analysis of 996 fishery-collected samples using the microsatellite panel revealed that larger rivers contribute greater numbers of individuals to the fishery and that coastal fisheries largely exploit individuals originating from nearby rivers, corroborating results from traditional tagging experiments. Our results demonstrate the efficacy of sequence-based microsatellite genotyping to advance understanding of fine-scale population structure and harvest composition in northern and understudied species.

Estimating the relative fitness of escaped farmed salmon offspring in the wild and modelling the consequences of invasion for wild populations.

Throughout their native range, wild Atlantic salmon populations are threatened by hybridization and introgression with escapees from net-pen salmon aquaculture. Although domestic-wild hybrid offspring have shown reduced fitness in laboratory and field experiments, consequential impacts on population abundance and genetic integrity remain difficult to predict in the field, in part because the strength of selection against domestic offspring is often unknown and context-dependent. Here, we follow a single large escape event of farmed Atlantic salmon in southern Newfoundland and monitor changes in the in-river proportions of hybrids and feral individuals over time using genetically based hybrid identification. Over a three-year period following the escape, the overall proportion of wild parr increased consistently (total wild proportion of 71.6%, 75.1% and 87.5% each year, respectively), with subsequent declines in feral (genetically pure farmed individuals originating from escaped, farmed adults) and hybrid parr. We quantify the strength of selection against parr of aquaculture ancestry and explore the genetic and demographic consequences for populations in the region. Within-cohort changes in the relative proportions of feral and F1 parr suggest reduced relative survival compared to wild individuals over the first (0.15 and 0.81 for feral and F1, respectively) and second years of life (0.26, 0.83). These relative survivorship estimates were used to inform an individual-based salmon eco-genetic model to project changes in adult abundance and overall allele frequency across three invasion scenarios ranging from short-term to long-term invasion and three relative survival scenarios. Modelling results indicate that total population abundance and time to recovery were greatly affected by relative survivorship and predict significant declines in wild population abundance under continued large escape events and calculated survivorship. Overall, this work demonstrates the importance of estimating the strength of selection against domestic offspring in the wild to predict the long-term impact of farmed salmon escape events on wild populations.

Extensive hybridization following a large escape of domesticated Atlantic salmon in the Northwest Atlantic.

Domestication is rife with episodes of interbreeding between cultured and wild populations, potentially challenging adaptive variation in the wild. In Atlantic salmon, Salmo salar, the number of domesticated individuals far exceeds wild individuals, and escape events occur regularly, yet evidence of the magnitude and geographic scale of interbreeding resulting from individual escape events is lacking. We screened juvenile Atlantic salmon using 95 single nucleotide polymorphisms following a single, large aquaculture escape in the Northwest Atlantic and report the landscape-scale detection of hybrid and feral salmon (27.1%, 17/18 rivers). Hybrids were reproductively viable, and observed at higher frequency in smaller wild populations. Repeated annual sampling of this cohort revealed decreases in the presence of hybrid and feral offspring over time. These results link previous observations of escaped salmon in rivers with reports of population genetic change, and demonstrate the potential negative consequences of escapes from net-pen aquaculture on wild populations.