Reconciling seascape genetics and fisheries science in three codistributed flatfishes

Uncertainty hampers innovative mixed‐fisheries management by the scales at which connectivity dynamics are relevant to management objectives. The spatial scale of sustainable stock management is species‐specific and depends on ecology, life history and population connectivity. One valuable approach to understand these spatial scales is to determine to what extent population genetic structure correlates with the oceanographic environment. Here, we compare the level of genetic connectivity in three codistributed and commercially exploited demersal flatfish species living in the North East Atlantic Ocean. Population genetic structure was analysed based on 14, 14 and 10 neutral DNA microsatellite markers for turbot, brill and sole, respectively. We then used redundancy analysis (RDA) to attribute the genetic variation to spatial (geographical location), temporal (sampling year) and oceanographic (water column characteristics) components. The genetic structure of turbot was composed of three clusters and correlated with variation in the depth of the pycnocline, in addition to spatial factors. The genetic structure of brill was homogenous, but correlated with average annual stratification and spatial factors. In sole, the genetic structure was composed of three clusters, but was only linked to a temporal factor. We explored whether the management of data poor commercial fisheries, such as in brill and turbot, might benefit from population‐specific information. We conclude that the management of fish stocks has to consider species‐specific genetic structures and may benefit from the documentation of the genetic seascape and life‐history traits.publishedVersionUnit Licence Agreemen

Reconciling seascape genetics and fisheries science in three codistributed flatfishes

Uncertainty hampers innovative mixed‐fisheries management by the scales at which connectivity dynamics are relevant to management objectives. The spatial scale of sustainable stock management is species‐specific and depends on ecology, life history and population connectivity. One valuable approach to understand these spatial scales is to determine to what extent population genetic structure correlates with the oceanographic environment. Here, we compare the level of genetic connectivity in three codistributed and commercially exploited demersal flatfish species living in the North East Atlantic Ocean. Population genetic structure was analysed based on 14, 14 and 10 neutral DNA microsatellite markers for turbot, brill and sole, respectively. We then used redundancy analysis (RDA) to attribute the genetic variation to spatial (geographical location), temporal (sampling year) and oceanographic (water column characteristics) components. The genetic structure of turbot was composed of three clusters and correlated with variation in the depth of the pycnocline, in addition to spatial factors. The genetic structure of brill was homogenous, but correlated with average annual stratification and spatial factors. In sole, the genetic structure was composed of three clusters, but was only linked to a temporal factor. We explored whether the management of data poor commercial fisheries, such as in brill and turbot, might benefit from population‐specific information. We conclude that the management of fish stocks has to consider species‐specific genetic structures and may benefit from the documentation of the genetic seascape and life‐history traits.publishedVersionUnit Licence Agreemen

Evaluating genetic traceability methods for captive bred marine fish and their applications in fisheries management and wildlife forensics

Growing demands for marine fish products is leading to increased pressure on already depleted wild populations and a rise in the aquaculture production. Consequently, more captive bred fish are released into the wild through accidental escape or deliberate restocking, stock enhancement and sea ranching programs. The increased mixing of captive bred fish with wild conspecifics may affect the ecological and/or genetic integrity of wild fish populations. From a fisheries management perspective unambiguous identification tools for captive bred fish will be highly valuable to manage risks. Additionally there is great potential to use these tools in wildlife forensics (i.e. tracing back escapees to their origin and determining mislabelling of seafood products). Using SNP data from captive bred and wild populations of Atlantic cod (Gadus morhua L.) and sole (Solea solea L.), we explored the efficiency of population and parentage assignment techniques for the identification and tracing of captive bred fish. Simulated and empirical data were used to correct for stochastic genetic effects. Overall, parentage assignment performed well when a large effective population size characterizes the broodstock and escapees originate from early generations of captive breeding. Consequently, parentage assignments are particularly useful from a fisheries management perspective to monitor the effects of deliberate releases of captive bred fish on wild populations. Population assignment proved to be more efficient after several generations of captive breeding, which makes it a useful method in forensic applications for well-established aquaculture species. We suggest the implementation of a case by case strategy when choosing the best method

Identification and validation of single nucleotide polymorphisms in growth- and maturation-related candidate genes in sole (Solea solea L.)

Genomic methodologies applied in evolutionary and fisheries research have been of great benefit to understand the marine ecosystem and the management of natural resources. Although single nucleotide polymorphisms (SNPs) are attractive for the study of local adaptation, spatial stock management and traceability, and investigating the effects of fisheries-induced selection, they have rarely been exploited in non-model organisms. This is partly due to difficulties in finding and validating SNPs in species with limited or no genomic resources. Complementary to random genome-scan approaches, a targeted candidate gene approach has the potential to unveil pre-selected functional diversity and provides more in depth information on the action of selection at specific genes. For example genes can be under selective pressure due to climate change and sustained periods of heavy fishing pressure. In this study, we applied a candidate gene approach in sole (Solea solea L.), an important member of the demersal ecosystem. As consumption flatfish it is heavy exploited and has experienced associated life-history changes over the last 60 years. To discover novel genetic polymorphisms in or around genes linked to important life history traits in sole, we screened a total of 76 candidate genes related to growth and maturation using a targeted resequencing approach. We identified in total 86 putative SNPs in 22 genes and validated 29 SNPs using a multiplex single-base extension genotyping assay. We found 22 informative SNPs, of which two represent non-synonymous mutations, potentially of functional relevance. These novel markers should be rapidly and broadly applicable in analyses of natural sole populations, as a measure of the evolutionary signature of overfishing and for initiatives on marker assisted selection

Syntenic relationships between sole and five other (flat)fish.

<p>The chromosomes of four model fish species, namely stickleback (S), tilapia (T), pufferfish (P) and medaka (M), were grouped in A- and B-groups according to their syntenic relationships as described in Sarropoulou et al. (2008), Kai et al. (2011) and Guyon et al. (2012) (left column). The numbers in the grid indicate the number of contigs where sequence homology was found between the sole linkage groups and the chromosomes of the four model species. For each chromosome the sole linkage group with the largest number of homologous sequences is highlighted in grey. Marked with *: the 21 linkage groups that are suggested as chromosome counterpart for sole (or at least part of it). In italics: linkage groups likely to be on the same chromosome as the linkage group marked with * to the left of it. For all 21 putative sole chromosomes (except for LG23) a homologous turbot linkage group is suggested (right column).</p

Number of markers, the corresponding number of distinct contigs and map length for each linkage group of sole.

<p>Number of markers, the corresponding number of distinct contigs and map length for each linkage group of sole.</p

Seascape genetics of a flatfish reveals local selection under high levels of gene flow

Local adaptation is often found to be in a delicate balance with gene flow in marine species with high dispersal potential. Genotyping with mapped transcriptome-derived markers and advanced seascape statistical analyses are proven tools to uncover the genomic basis of biologically relevant traits under environmental selection. Using a panel of 426 gene-linked single nucleotide polymorphisms (SNPs), we scanned 17 samples (n = 539) of sole (Solea solea L.) from the Northeast Atlantic Ocean and applied a node-based seascape analysis. Neutral loci confirmed a clear distinction between the North Sea–Baltic Sea transition zone and the other Eastern Atlantic samples. At a more subtle level, the latter unit split in an English Channel and North Sea group, and a Bay of Biscay and Atlantic Iberian coast group. A fourth group, the Irish and Celtic Sea, was identified with 19 outlier loci. A pattern of isolation by distance (IBD) characterized the latitudinal distribution. Seascape analyses identified winter seawater temperature, food availability and coastal currents to explain a significant component of geographically distributed genetic variation, suggesting that these factors act as drivers of local adaptation. The evidence for local adaptation is in line with the current understanding on the impact of two key ecological factors, the life-history trait winter mortality and the behaviour of inshore/offshore spawning. We conclude that the subtle differentiation between two metapopulations (North Sea and Bay of Biscay) mirrors local adaptation. At least three genomic regions with strong population differentiation point to locally divergent selection. Further functional characterization of these genomic regions should help with formulating adaptive management policies.status: publishe

Sex-averaged linkage map of sole.

<p>Map distances are calculated using the Kosambi mapping function and shown in centimorgans. Combined SNPs are indicated with a ‘C’ at the beginning of their name.</p

Seascape genetics of a flatfish reveals local selection under high levels of gene flow

Local adaptation is often found to be in a delicate balance with gene flow in marine species with high dispersal potential. Genotyping with mapped transcriptome-derived markers and advanced seascape statistical analyses are proven tools to uncover the genomic basis of biologically relevant traits under environmental selection. Using a panel of 426 gene-linked single nucleotide polymorphisms (SNPs), we scanned 17 samples (n¼539) of sole (Solea solea L.) from the Northeast Atlantic Ocean and applied a node-based seascape analysis. Neutral loci confirmed a clear distinction between the North Sea–Baltic Sea transition zone and the other Eastern Atlantic samples. At a more subtle level,the latter unit split in an English Channel and North Sea group, and a Bay of Biscay and Atlantic Iberian coast group. A fourth group, the Irish and Celtic Sea, was identified with 19 outlier loci. A pattern of isolation by distance (IBD) characterized the latitudinal distribution. Seascape analyses identified winter seawater temperature, food availability and coastal currents to explain a significant component of geographically distributed genetic variation, suggesting that these factors act as drivers of local adaptation. The evidence for local adaptation is in line with the current understanding on the impact of two key ecological factors, the life-history trait winter mortality and the behaviour of inshore/offshore spawning. We conclude that the subtle differentiation between two metapopulations (North Sea and Bay of Biscay) mirrors local adaptation.At least three genomic regions with strong population differentiation point to locally divergent selection. Further functional characterization of these genomic regions should help with formulating adaptive management policies