Genetic structure of West Greenland populations of lumpfish Cyclopterus lumpus

Genomic Signatures of Local Adaptation under High Gene Flow in Lumpfish-Implications for Broodstock Provenance Sourcing and Larval Production

Aquaculture of the lumpfish (Cyclopterus lumpus L.) has become a large, lucrative industry owing to the escalating demand for "cleaner fish" to minimise sea lice infestations in Atlantic salmon mariculture farms. We used over 10K genome-wide single nucleotide polymorphisms (SNPs) to investigate the spatial patterns of genomic variation in the lumpfish along the coast of Norway and across the North Atlantic. Moreover, we applied three genome scans for outliers and two genotype-environment association tests to assess the signatures and patterns of local adaptation under extensive gene flow. With our 'global' sampling regime, we found two major genetic groups of lumpfish, i.e., the western and eastern Atlantic. Regionally in Norway, we found marginal evidence of population structure, where the population genomic analysis revealed a small portion of individuals with a different genetic ancestry. Nevertheless, we found strong support for local adaption under high gene flow in the Norwegian lumpfish and identified over 380 high-confidence environment-associated loci linked to gene sets with a key role in biological processes associated with environmental pressures and embryonic development. Our results bridge population genetic/genomics studies with seascape genomics studies and will facilitate genome-enabled monitoring of the genetic impacts of escapees and allow for genetic-informed broodstock selection and management in Norway.

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.

Global, regional, and cryptic population structure in a high gene-flow transatlantic fish

Lumpfish (Cyclopterus lumpus) is a transatlantic marine fish displaying large population sizes and a high potential for dispersal and gene-flow. These features are expected to result in weak population structure. Here, we investigated population genetic structure of lumpfish throughout its natural distribution in the North Atlantic using two approaches: I) 4,393 genome wide SNPs and 95 individuals from 10 locations, and II) 139 discriminatory SNPs and 1,669 individuals from 40 locations. Both approaches identified extensive population genetic structuring with a major split between the East and West Atlantic and a distinct Baltic Sea population, as well as further differentiation of lumpfish from the English Channel, Iceland, and Greenland. The discriminatory loci displayed ~2-5 times higher divergence than the genome wide approach, revealing further evidence of local population substructures. Lumpfish from Isfjorden in Svalbard were highly distinct but resembled most fish from Greenland. The Kattegat area in the Baltic transition zone, formed a previously undescribed distinct genetic group. Also, further subdivision was detected within North America, Iceland, West Greenland, Barents Sea, and Norway. Although lumpfish have considerable potential for dispersal and gene-flow, the observed high levels of population structuring throughout the Atlantic suggests that this species may have a natal homing behavior and local populations with adaptive differences. This fine-scale population structure calls for consideration when defining management units for exploitation of lumpfish stocks and in decisions related to sourcing and moving lumpfish for cleaner fish use in salmonid aquaculture.

Genome- and transcriptome-derived microsatellite loci in lumpfish Cyclopterus lumpus: molecular tools for aquaculture, conservation and fisheries management

The lumpfish Cyclopterus lumpus is commercially exploited in numerous areas of its range in the North Atlantic Ocean, and is important in salmonid aquaculture as a biological agent for controlling sea lice. Despite the economic importance, few genetic resources for downstream applications, such as linkage mapping, parentage analysis, marker-assisted selection (MAS), quantitative trait loci (QTL) analysis, and assessing adaptive genetic diversity are currently available for the species. Here, we identify both genome- and transcriptome-derived microsatellites loci from C. lumpus to facilitate such applications. Across 2,346 genomic contigs, we detected a total of 3,067 microsatellite loci, of which 723 were the most suitable ones for primer design. From 116,555 transcriptomic unigenes, we identified a total of 231,556 microsatellite loci, which may indicate a high coverage of the available STRs. Out of these, primer pairs could only be designed for 6,203 loci. Dinucleotide repeats accounted for 89 percent and 52 percent of the genome- and transcriptome-derived microsatellites, respectively. The genetic composition of the dominant repeat motif types showed differences from other investigated fish species. In the genome-derived microsatellites AC/GT (67.8 percent), followed by AG/CT (15.1 percent) and AT/AT (5.6 percent) were the major motifs. Transcriptome-derived microsatellites showed also most dominantly the AC/GT repeat motif (33 percent), followed by A/T (26.6 percent) and AG/CT (11 percent). Functional annotation of microsatellite-containing transcriptomic sequences showed that the majority of the expressed sequence tags encode proteins involved in cellular and metabolic processes, binding activity and catalytic reactions. Importantly, STRs linked to genes involved in immune system process, growth, locomotion and reproduction were discovered in the present study. The extensive genomic marker information reported here will facilitate molecular ecology studies, conservation initiatives and will benefit many aspects of the breeding programmes of C. lumpus.

Genetic and phenotypic differentiation of lumpfish (Cyclopterus lumpus) across the North Atlantic: implications for conservation and aquaculture

Demand for lumpfish (Cyclopterus lumpus) has increased exponentially over the last decade, both for their roe, which is used as a caviar substitute, and increasingly also as cleaner fish to control sea lice in salmon farming. The species is classified as Near Threatened by the IUCN and there are growing concerns that over-exploitation of wild stocks and translocation of hatchery-reared lumpfish may compromise the genetic diversity of native populations. We carried out a comparative analysis of genetic and phenotypic variation across the species' range to estimate the level of genetic and phenotypic differentiation, and determined patterns of gene flow at spatial scales relevant to management. We found five genetically distinct groups located in the West Atlantic (USA and Canada), Mid Atlantic (Iceland), East Atlantic (Faroe Islands, Ireland, Scotland, Norway and Denmark), English Channel (England) and Baltic Sea (Sweden). Significant phenotypic differences were also found, with Baltic lumpfish growing more slowly, attaining a higher condition factor and maturing at a smaller size than North Atlantic lumpfish. Estimates of effective population size were consistently low across the North East Atlantic (Iceland, Faroe Islands and Norway), the area where most wild lumpfish are fished for their roe, and also for the aquaculture industry. Our study suggests that some lumpfish populations are very small and have low genetic diversity, which makes them particularly vulnerable to over-exploitation and genetic introgression. To protect them we advocate curtailing fishing effort, closing the breeding cycle of the species in captivity to reduce dependence on wild stocks, restricting the translocation of genetically distinct populations, and limiting the risk of farm escapes.