Movement diversity and partial sympatry of coastal and Northeast Arctic cod ecotypes at high latitudes.

Movement diversity within species represent an important but often neglected, component of biodiversity that affects ecological and genetic interactions, as well as the productivity of exploited systems. By combining individual tracking data from acoustic telemetry with novel genetic analyses, we describe the movement diversity of two Atlantic cod Gadus morhua ecotypes in two high-latitude fjord systems: the highly migratory Northeast Arctic cod (NEA cod) that supports the largest cod fishery in the world, and the more sedentary Norwegian coastal cod, which is currently in a depleted state. As predicted, coastal cod displayed a higher level of fjord residency than NEA cod. Of the cod tagged during the spawning season, NEA cod left the fjords permanently to a greater extent and earlier compared to coastal cod, which to a greater extent remained resident and left the fjords temporarily. Despite this overall pattern, horizontal movements atypical for the ecotypes were common with some NEA cod remaining within the fjords year-round and some coastal cod displaying a low fjord fidelity. Fjord residency and exit timing also differed with spawning status and body size, with spawning cod and large individuals tagged during the feeding season more prone to leave the fjords and earlier than non-spawning and smaller individuals. While our results confirm a lower fjord dependency for NEA cod, they highlight a movement diversity within each ecotype and sympatric residency between ecotypes, previously undetected by population-level monitoring. This new knowledge is relevant for the management, which should base their fisheries advice for these interacting ecotypes on their habitat use and seasonal movements.

Contaminants reach everywhere: Fish dietary samples should be surface decontaminated prior to molecular diet analysis.

Knowledge of trophic interaction is necessary to understand the dynamics of ecosystems and develop ecosystem-based management. The key data to measure these interactions should come from large-scale diet analyses with good taxonomic resolution. To that end, molecular methods that analyze prey DNA from guts and feces provide high-resolution dietary taxonomic data. However, molecular diet analysis may also produce unreliable results if the samples are contaminated by external sources of DNA. Employing the freshwater European whitefish (Coregonus lavaretus) as a tracer for sample contamination, we studied the possible route of whitefish in beaked redfish (Sebastes mentella) guts sampled in the Barents Sea. We used whitefish-specific COI primers for diagnostic analysis, and fish-specific 12S and metazoa-specific COI primers for metabarcoding analyses of intestine and stomach contents of fish samples that were either not cleaned, water cleaned, or bleach cleaned after being in contact with whitefish. Both the diagnostic and COI metabarcoding revealed clear positive effects of cleaning samples as whitefish were detected in significantly higher numbers of uncleaned samples compared to water or bleach-cleaned samples. Stomachs were more susceptible to contamination than intestines and bleach cleaning reduced the frequency of whitefish contamination. Also, the metabarcoding approach detected significantly more reads of whitefish in the stomach than in intestine samples. The diagnostic analysis and COI metabarcoding detected contaminants in a higher and comparable number of gut samples than the 12S-based approach. Our study underlines thus the importance of surface decontamination of aquatic samples to obtain reliable diet information from molecular data.

Geographic variation in gene flow from a genetically distinct migratory ecotype drives population genetic structure of coastal Atlantic cod (Gadus morhua L.).

Identifying how physical and biotic factors shape genetic connectivity among populations in time and space is essential to our understanding of the evolutionary trajectory as well as the management of marine species. Atlantic cod is a widespread and commercially important marine species displaying several ecotypes with different life history strategies. Using three sets of SNPs: neutral, informative, and genome-inversion linked, we studied population genetic structure of ~2500 coastal Atlantic cod (CC) from 40 locations along Norway's 2500 km coastline, including nine fjords. We observed: (1) a genetic cline, suggesting a mechanism of isolation by distance, characterized by a declining F ST between CC and North East Arctic Cod (NEAC-genetically distinct migratory ecotype) with increasing latitude, (2) that in the north, samples of CC from outer-fjord areas were genetically more similar to NEAC than were samples of CC from their corresponding inner-fjord areas, (3) greater population genetic differentiation among CC sampled from outer-fjord areas along the coast, than among CC sampled from their corresponding inner-fjord areas, (4) genetic differentiation among samples of CC from both within and among fjords. Collectively, these results permit us to draw two main conclusions. First, that differences in the relative presence of the genetically highly distinct, migratory ecotype NEAC, declining from north to south and from outer to inner fjord, plays the major role in driving population genetic structure of the Norwegian CC. Second, that there is limited connectivity between CC from different fjords. These results suggest that the current management units implemented for this species in Norway should be divided into smaller entities. Furthermore, the situation where introgression from one ecotype drives population genetic structure of another, as is the case here, may exist in other species and geographical regions, thus creating additional challenges for sustainable fisheries management.

Hierarchical genetic structure in an evolving species complex: Insights from genome wide ddRAD data in Sebastes mentella.

The diverse biology and ecology of marine organisms may lead to complex patterns of intraspecific diversity for both neutral and adaptive genetic variation. Sebastes mentella displays a particular life-history as livebearers, for which existence of multiple ecotypes has been suspected to complicate the genetic population structure of the species. Double digest restriction-site associated DNA was used to investigate genetic population structure in S. mentella and to scan for evidence of selection. In total, 42,288 SNPs were detected in 277 fish, and 1,943 neutral and 97 tentatively adaptive loci were selected following stringent filtration. Unprecedented levels of genetic differentiation were found among the previously defined 'shallow pelagic', 'deep pelagic' and 'demersal slope' ecotypes, with overall mean FST = 0.05 and 0.24 in neutral and outlier SNPs, respectively. Bayesian computation estimated a concurrent and historical divergence among these three ecotypes and evidence of local adaptation was found in the S. mentella genome. Overall, these findings imply that the depth-defined habitat divergence of S. mentella has led to reproductive isolation and possibly adaptive radiation among these ecotypes. Additional sub-structuring was detected within the 'shallow' and 'deep' pelagic ecotypes. Population assignment of individual fish showed more than 94% agreement between results based on SNP and previously generated microsatellite data, but the SNP data provided a lower estimate of hybridization among the ecotypes than that by microsatellite data. We identified a SNP panel with only 21 loci to discriminate populations in mixed samples based on a machine-learning algorithm. This first SNP based investigation clarifies the population structure of S. mentella, and provides novel and high-resolution genomic tools for future investigations. The insights and tools provided here can readily be incorporated into the management of S. mentella and serve as a template for other exploited marine species exhibiting similar complex life history traits.

Genomic analysis reveals neutral and adaptive patterns that challenge the current management regime for East Atlantic cod Gadus morhua L.

Challenging long-held perceptions of fish management units can help to protect vulnerable stocks. When a fishery consisting of multiple genetic stocks is managed as a single unit, overexploitation and depletion of minor genetic units can occur. Atlantic cod (Gadus morhua) is an economically and ecologically important marine species across the North Atlantic. The application of new genomic resources, including SNP arrays, allows us to detect and explore novel structure within specific cod management units. In Norwegian waters, coastal cod (i.e. those not undertaking extensive migrations) are divided into two arbitrary management units defined by ICES: one between 62° and 70°N (Norwegian coastal cod; NCC) and one between 58° and 62°N (Norwegian coastal south; NCS). Together, these capture a fishery area of >25,000 km2 containing many spawning grounds. To assess whether these geographic units correctly represent genetic stocks, we analysed spawning cod of NCC and NCS for more than 8,000 SNPs along with samples of Russian White Sea cod, north-east Arctic cod (NEAC: the largest Atlantic stock), and outgroup samples representing the Irish and Faroe Sea's. Our analyses revealed large differences in spatial patterns of genetic differentiation across the genome and revealed a complex biological structure within NCC and NCS. Haplotype maps from four chromosome sets show regional specific SNP indicating a complex genetic structure. The current management plan dividing the coastal cod into only two management units does not accurately reflect the genetic units and needs to be revised. Coastal cod in Norway, while highly heterogenous, is also genetically distinct from neighbouring stocks in the north (NEAC), west (Faroe Island) and the south. The White Sea cod are highly divergent from other cod, possibly yielding support to the earlier notion of subspecies rank.

The pantophysin gene and its relationship with survival in early life stages of Atlantic cod.

Genetic markers are widely used in fisheries management around the world. While the genetic structure and markers selected are usually based on samples from the wild, very few controlled experiments have been carried out to investigate possible differences in influence on traits between markers. Here we examine the bi-allelic gene pantophysin (Pan I), widely used in the management of Atlantic cod, in a series of in vitro crosses under a range of temperatures. It has been proposed that this gene, or another tightly linked gene, may be under strong divergent selection. Resolving this issue is essential in order to interpret results when using this gene marker for stock management. We found no evidence of departure from the expected 1 : 2 : 1 Mendelian ratio for any of the three genotypes during the egg stage, while both the 6 and 12°C temperature regimes in tank experiments favoured the survival of the Pan IAA genotype. No difference in genotype survival was, however, found in a more natural mesocosm environment. Collectively, these results suggest that for the early life stages of Atlantic cod, and under the current experimental conditions, there is no strong consistent influence of Pan I genotype on survival. The results also emphasize the importance of varied experimental studies to verify the importance of environmental factors influencing genotype selection.

Analysis of coastal cod (Gadus morhua L.) sampled on spawning sites reveals a genetic gradient throughout Norway's coastline.

BACKGROUND: Atlantic cod (Gadus morhua L.) has formed the basis of many economically significant fisheries in the North Atlantic, and is one of the best studied marine fishes, but a legacy of overexploitation has depleted populations and collapsed fisheries in several regions. Previous studies have identified considerable population genetic structure for Atlantic cod. However, within Norway, which is the country with the largest remaining catch in the Atlantic, the population genetic structure of coastal cod (NCC) along the entire coastline has not yet been investigated. We sampled > 4000 cod from 55 spawning sites. All fish were genotyped with 6 microsatellite markers and Pan I (Dataset 1). A sub-set of the samples (1295 fish from 17 locations) were also genotyped with an additional 9 microsatellites (Dataset 2). Otoliths were read in order to exclude North East Arctic Cod (NEAC) from the analyses, as and where appropriate. RESULTS: We found no difference in genetic diversity, measured as number of alleles, allelic richness, heterozygosity nor effective population sizes, in the north-south gradient. In both data sets, weak but significant population genetic structure was revealed (Dataset 1: global FST = 0.008, P < 0.0001. Dataset 2: global FST = 0.004, P < 0.0001). While no clear genetic groups were identified, genetic differentiation increased among geographically-distinct samples. Although the locus Gmo132 was identified as a candidate for positive selection, possibly through linkage with a genomic region under selection, overall trends remained when this locus was excluded from the analyses. The most common allele in loci Gmo132 and Gmo34 showed a marked frequency change in the north-south gradient, increasing towards the frequency observed in NEAC in the north. CONCLUSION: We conclude that Norwegian coastal cod displays significant population genetic structure throughout its entire range, that follows a trend of isolation by distance. Furthermore, we suggest that a gradient of genetic introgression between NEAC and NCC contributes to the observed population genetic structure. The current management regime for coastal cod in Norway, dividing it into two stocks at 62°N, represents a simplification of the level of genetic connectivity among coastal cod in Norway, and needs revision.

Geographic extent of introgression in Sebastes mentella and its effect on genetic population structure.

Genetic population structure is often used to identify management units in exploited species, but the extent of genetic differentiation may be inflated by geographic variation in the level of hybridization between species. We identify the genetic population structure of Sebastes mentella and investigate possible introgression within the genus by analyzing 13 microsatellites in 2,562 redfish specimens sampled throughout the North Atlantic. The data support an historical divergence between the "shallow" and "deep" groups, beyond the Irminger Sea where they were described previously. A third group, "slope," has an extended distribution on the East Greenland Shelf, in addition to earlier findings on the Icelandic slope. Furthermore, S. mentella from the Northeast Arctic and Northwest Atlantic waters are genetically different populations. In both areas, interspecific introgression may influence allele frequency differences among populations. Evidence of introgression was found for almost all the identified Sebastes gene pools, but to a much lower extent than suggested earlier. Greenland waters appear to be a sympatric zone for many of the genetically independent Sebastes groups. This study illustrates that the identified groups maintain their genetic integrity in this region despite introgression.

Two adjacent inversions maintain genomic differentiation between migratory and stationary ecotypes of Atlantic cod.

Atlantic cod is composed of multiple migratory and stationary populations widely distributed in the North Atlantic Ocean. The Northeast Arctic cod (NEAC) population in the Barents Sea undertakes annual spawning migrations to the northern Norwegian coast. Although spawning occurs sympatrically with the stationary Norwegian coastal cod (NCC), phenotypic and genetic differences between NEAC and NCC are maintained. In this study, we resolve the enigma by revealing the mechanisms underlying these differences. Extended linkage disequilibrium (LD) and population divergence were demonstrated in a 17.4-Mb region on linkage group 1 (LG1) based on genotypes of 494 SNPs from 192 parents of farmed families of NEAC, NCC or NEACxNCC crosses. Linkage analyses revealed two adjacent inversions within this region that repress meiotic recombination in NEACxNCC crosses. We identified a NEAC-specific haplotype consisting of 186 SNPs that was fixed in NEAC sampled from the Barents Sea, but segregating under Hardy-Weinberg equilibrium in eight NCC stocks. Comparative genomic analyses determine the NEAC configuration of the inversions to be the derived state and date it to ~1.6-2.0 Mya. The haplotype block harbours 763 genes, including candidates regulating swim bladder pressure, haem synthesis and skeletal muscle organization conferring adaptation to long-distance migrations and vertical movements down to large depths. Our results suggest that the migratory ecotype experiences strong directional selection for the two adjacent inversions on LG1. Despite interbreeding between NEAC and NCC, the inversions are maintaining genetic differentiation, and we hypothesize the co-occurrence of multiple adaptive alleles forming a 'supergene' in the NEAC population.

Genetically distinct populations of northern shrimp, Pandalus borealis, in the North Atlantic: adaptation to different temperatures as an isolation factor.

The large-scale population genetic structure of northern shrimp, Pandalus borealis, was investigated over the species' range in the North Atlantic, identifying multiple genetically distinct groups. Genetic divergence among sample localities varied among 10 microsatellite loci (range: FST = -0.0002 to 0.0475) with a highly significant average (FST = 0.0149; P < 0.0001). In contrast, little or no genetic differences were observed among temporal replicates from the same localities (FST = 0.0004; P = 0.33). Spatial genetic patterns were compared to geographic distances, patterns of larval drift obtained through oceanographic modelling, and temperature differences, within a multiple linear regression framework. The best-fit model included all three factors and explained approximately 29% of all spatial genetic divergence. However, geographic distance and larval drift alone had only minor effects (2.5-4.7%) on large-scale genetic differentiation patterns, whereas bottom temperature differences explained most (26%). Larval drift was found to promote genetic homogeneity in parts of the study area with strong currents, but appeared ineffective across large temperature gradients. These findings highlight the breakdown of gene flow in a species with a long pelagic larval phase (up to 3 months) and indicate a role for local adaptation to temperature conditions in promoting evolutionary diversification and speciation in the marine environment.

A genomic island linked to ecotype divergence in Atlantic cod.

The genomic architecture underlying ecological divergence and ecological speciation with gene flow is still largely unknown for most organisms. One central question is whether divergence is genome-wide or localized in 'genomic mosaics' during early stages when gene flow is still pronounced. Empirical work has so far been limited, and the relative impacts of gene flow and natural selection on genomic patterns have not been fully explored. Here, we use ecotypes of Atlantic cod to investigate genomic patterns of diversity and population differentiation in a natural system characterized by high gene flow and large effective population sizes, properties which theoretically could restrict divergence in local genomic regions. We identify a genomic region of strong population differentiation, extending over approximately 20 cM, between pairs of migratory and stationary ecotypes examined at two different localities. Furthermore, the region is characterized by markedly reduced levels of genetic diversity in migratory ecotype samples. The results highlight the genomic region, or 'genomic island', as potentially associated with ecological divergence and suggest the involvement of a selective sweep. Finally, we also confirm earlier findings of localized genomic differentiation in three other linkage groups associated with divergence among eastern Atlantic populations. Thus, although the underlying mechanisms are still unknown, the results suggest that 'genomic mosaics' of differentiation may even be found under high levels of gene flow and that marine fishes may provide insightful model systems for studying and identifying initial targets of selection during ecological divergence.

Genomic signatures of local directional selection in a high gene flow marine organism; the Atlantic cod (Gadus morhua).

BACKGROUND: Marine fishes have been shown to display low levels of genetic structuring and associated high levels of gene flow, suggesting shallow evolutionary trajectories and, possibly, limited or lacking adaptive divergence among local populations. We investigated variation in 98 gene-associated single nucleotide polymorphisms (SNPs) for evidence of selection in local populations of Atlantic cod (Gadus morhua L.) across the species distribution. RESULTS: Our global genome scan analysis identified eight outlier gene loci with very high statistical support, likely to be subject to directional selection in local demes, or closely linked to loci under selection. Likewise, on a regional south/north transect of central and eastern Atlantic populations, seven loci displayed strongly elevated levels of genetic differentiation. Selection patterns among populations appeared to be relatively widespread and complex, i.e. outlier loci were generally not only associated with one of a few divergent local populations. Even on a limited geographical scale between the proximate North Sea and Baltic Sea populations four loci displayed evidence of adaptive evolution. Temporal genome scan analysis applied to DNA from archived otoliths from a Faeroese population demonstrated stability of the intra-population variation over 24 years. An exploratory landscape genetic analysis was used to elucidate potential effects of the most likely environmental factors responsible for the signatures of local adaptation. We found that genetic variation at several of the outlier loci was better correlated with temperature and/or salinity conditions at spawning grounds at spawning time than with geographic distance per se. CONCLUSION: These findings illustrate that adaptive population divergence may indeed be prevalent despite seemingly high levels of gene flow, as found in most marine fishes. Thus, results have important implications for our understanding of the interplay of evolutionary forces in general, and for the conservation of marine biodiversity under rapidly increasing evolutionary pressure from climate and fisheries induced changes in local environments.

Bathymetric barriers promoting genetic structure in the deepwater demersal fish tusk (Brosme brosme).

Population structuring in the North Atlantic deepwater demersal fish tusk (Brosme brosme) was studied with microsatellite DNA analyses. Screening eight samples from across the range of the species for seven loci revealed low but significant genetic heterogeneity (F(ST) = 0.0014). Spatial genetic variability was only weakly related to geographical (Euclidean) distance between study sites or separation of study sites along the path of major ocean currents. Instead, we found a significant effect of habitat, indicated by significant differentiation between relatively closely spaced sites: Rockall, which is surrounded by very deep water (>1000 m), and the Mid-Atlantic Ridge, which is separated from the European slope by a deep ocean basin, were differentiated from relatively homogeneous sites across the Nordic Seas. Limited adult migration across bathymetric barriers in combination with limited intersite exchange of pelagic eggs and larvae due to site-specific circulatory retention or poor survival during drift phases across deep basins may be reducing gene flow. We regard these limitations to gene flow as the most likely mechanisms for the observed population structure in this demersal species. The results underscore the importance of habitat boundaries in marine species.

New haemoglobin genotypes in Atlantic cod, Gadus morhua: possible relation with growth.

In a preliminary study, 121 individually tagged juvenile Atlantic cod (Gadus morhua) were classified according to their haemoglobin genotypes into four groups, i.e., two main haemoglobin genotypes [Hb-I(1/2), Hb-I(2/2)] and two sub-types [Hb-I(1/2b), Hb-I(2/2b)], and reared for 3 months at 10 degrees C, 13 degrees C and T-step (fish reared at 16 degrees C and then subsequently moved to 13 and later to 10 degrees C). Overall growth rates across temperatures were 10% and 19% higher in the Hb-I(2/2b), Hb-I(1/2b) sub-types compared to corresponding Hb-I(2/2) and Hb-I(1/2) main types, respectively. Individual growth rate trajectories varied between the genotypes at all temperatures studied. Our study indicates that under certain environmental conditions higher growth in the two sub-types compared to the main genotypes could be expected. This may indicate difference in other physiological characters not studied here, but seen in previous studies, i.e., oxygen affinity and competitive performance.

Strict conservation of the ITS regions of the ribosomal RNA genes in Atlantic cod (Gadus morhua L.).

The nucleotide sequence of the internal transcribed region (ITS) of ribosomal RNA genes from Atlantic cod (Gadus morhua L.) was determined. The complete ITS region spanned approximately 1113 base pairs. The ITS1 region comprised 532 base pairs, the 5.8S region 159, and the ITS2 region contained 422 base pairs. Sequence data were obtained from a total of 12 samples, one pool from six cod and 11 individuals. The sequencing was carried out in two separate experimental periods employing slightly different methodology. The samples were from two different cod stocks, Norwegian costal cod and North East Arctic cod. The sequence analysis showed that in the 12 samples, the ITS region, including the 5.8S RNA, was identical. The ITS region is thus totally conserved in these two cod stocks. The extreme conservation of the ITS regions in the cod rDNA could reflect the small genome size of cod and/or indicate a specific critical role in the processing of the ribosomal units in cold-adapted species.

Effects of temperature acclimation on lactate dehydrogenase of cod (Gadus morhua): genetic, kinetic and thermodynamic aspects.

The aim of this study was to determine the effects of seasonal temperature variation on the functional properties of lactate dehydrogenase (LDH) from white muscle and liver of Norwegian coastal cod (Gadus morhua) and the possible relevance of LDH allelic variability for thermal acclimation. Two groups of fishes were acclimated to 4 degrees C or 12 degrees C for one year. Polymorphism was observed in only one (Ldh-B) of the three Ldh loci expressed in cod liver and/or muscle. Isozyme expression remained unchanged regardless of acclimation temperature (T(A)). The products of locus Ldh-B comprise only 14-19% (depending on the tissue) of total LDH activities and, consequently, differences between phenotypes are negligible in terms of their effect on LDH total performance. No kinetic (, V(max)) or thermodynamic (E(a), DeltaG) differences were found among Ldh-B phenotypes. Clear kinetic differences were observed between LDH isoforms in the two tissues. However, the Arrhenius activation energy (E(a)) for pyruvate reduction was the same for both tissues (E(a)=47 kJ mol(-1)) at T(A)=12 degrees C. Factors T(A), tissue and phenotype did not reveal a significant effect on the Gibbs free energy change (DeltaG) of the reaction (55.5 kJ mol(-1)). However, at T(A)=4 degrees C, the E(a) was increased (E(a)=53-56 kJ mol(-1)) and the temperature dependence of the constant of substrate inhibition for pyruvate () decreased in both muscle and liver. In conclusion, the strategies of LDH adjustment to seasonal temperature variations in cod involve changes in LDH concentration (quantitative), adjustment of thermodynamic (E(a)) and kinetic () properties of the LDH (modulative) but not the expression of alternative isoforms (qualitative). We assume that the observed increase in E(a) and the decrease of temperature dependence of at low T(A) is the result of structural changes of the LDH molecule (temperature-driven protein folding). We propose a new mechanism of metabolic compensation of seasonal temperature variations - cold acclimation results in changes in the kinetic and thermodynamic properties of LDH in a way that favours aerobic metabolism through reduction of the competition of LDH for pyruvate in normoxic conditions.