Population structure in Atlantic cod in the eastern North Sea-Skagerrak-Kattegat: early life stage dispersal and adult migration
BMC Research Notes volume 9, Article number: 63 (2016)
In marine fish species, where pelagic egg and larvae drift with ocean currents, population structure has been suggested to be maintained by larval retention due to hydrographic structuring and by homing of adult fish to natal areas. Whilst natal homing of adults has been demonstrated for anadromous and coral reef fishes, there are few documented examples of philopatric migration in temperate marine fish species.
Here, we demonstrate temporally stable genetic differentiation among spawning populations of Atlantic cod (Gadus morhua L.), and present genetic and behavioural evidence for larval drift and philopatric migration in the eastern North Sea-Skagerrak-Kattegat area. We show that juvenile cod collected in the eastern Skagerrak and central Kattegat are genetically similar to cod from offshore spawning areas in the eastern North Sea. Genetic assignment of individual 2–5 year old fish indicates that cod residing at, or migrating towards, spawning areas in Kattegat and the North Sea display philopatric behaviours.
Together these findings suggest a loop between spawning, larval drift and adult return-migrations to spawning areas and underlines that both oceanographic processes and migratory behaviour in the adult phase may be important for stock separation and integrity in marine temperate fishes such as Atlantic cod.
Recognizing temporal and spatial structuring in marine fishes and understanding of the mechanisms responsible for population structure are of vital importance for fisheries management and species conservation [1–4]. Structuring mechanisms include fish migratory behaviour , in particular spawning site fidelity such as natal homing of mature adults , and physical forcing, where hydrodynamic stratification, currents and eddies leads to dispersal, retention or settlement of eggs and larvae [7–10]. Interactions between oceanographic and environmental features, such as bathymetry, temperature and salinity, with egg buoyancy and larval behaviour are thus considered to be of high importance for connectivity and gene flow, as they set the limits for population differentiation .
Because marine fish species often show ontogenetic and seasonal shifts in habitat use, the entire life cycle has to be considered when studying processes shaping population structure [12–14]. Spawning-site fidelity, i.e., adult fish returning repeatedly to spawn at the same location, irrespectively whether they were hatched at this location or not, is widely recognized in many species, e.g., Atlantic Cod, Gadus morhua . Natal homing, where fish actively return and spawn where they were born, on the other hand, has so far been demonstrated only in a few cases, such as coral reef fishes [15–18], salmonids  and Atlantic bluefin tuna, Thunnus thynnus, , although strong circumstantial evidence has been presented for e.g., eels, Anguilla sp. , Atlantic herring, Clupea harengus, , pike, Esox lucius,  and Atlantic cod [24–27]. The distinction between spawning site fidelity and natal homing is fundamental because natal homing behaviour is expected to contribute to population differentiation given sufficient time, whereas spawning-site fidelity may not necessarily lead to reproductive isolation if there is extensive egg and larval drift together with opportunistic and non-philopatric recruitment of juveniles to adult aggregations [13, 28]. Demonstrating natal homing in marine fish species that disperse both passively as planktonic larvae, and by active locomotion as adults, is methodologically challenging and requires that individuals can be tracked from the egg stage to maturation and reproduction, or that reproducing individuals can be assigned to specific spawning populations using distinguishable natal tags [12, 29]. To date, natal homing has been inferred using natural tags such as elemental [6, 15, 20, 23, 30] or genetic fingerprints [17, 31].
Atlantic cod is a marine fish that exhibits population structuring on both large and small spatial scales [32, 33]. In the North Sea-Skagerrak-Kattegat area Atlantic cod comprise a mixture of co-existing resident forms completing their entire life cycle in fjords or sheltered areas , and oceanic populations [25, 30, 35–39]. While adult cod abundance has declined dramatically in the eastern inshore Skagerrak and Kattegat since the 1980′s, juvenile cod show no such trend [36, 40]. Genetic analyses [34, 41] in combination with oceanographic modelling  suggest an extensive drift of cod larvae from the North Sea into coastal Skagerrak, where they settle and possibly mix with juveniles of local coastal origin. The inflow of larvae does not, however, seem to support the diminishing coastal Skagerrak or Kattegat stocks [36, 40]. Rather, it has been hypothesized that the present low abundance of adult cod in eastern Skagerrak is due to return migrations to natal areas in the North Sea [25, 43–45].
Here, we further investigated the population structure and the mechanisms of dispersal and connectivity between cod populations in the eastern North Sea, Skagerrak and Kattegat. First, we investigated the spatial and temporal genetic structure by targeting adult fish at spawning. Second, we examined, by genetic analyses, the temporal variation in larval drift from the eastern North Sea by genetic analyses of juvenile cod in eastern Skagerrak and Kattegat collected in 2005 and 2011. Third, as both a temporally stable population structure was observed, and juveniles were found to be dispersed between areas, it was concluded that alternative structuring mechanisms to physical forcing are likely to be at play for cod in the area. For that reason, we tested whether philopatric migration behaviour, including natal homing, may act as a stock structuring process in cod by matching DNA profiles of individual migrating fish with the genetic information from cod populations on spawning grounds. Using archival tags, Svedäng et al.  showed that 2–5 year old cod released in eastern Skagerrak undertake directional migration towards the North Sea during the spawning season, indicating migration to potential natal spawning grounds . Moreover, Svedäng et al.  identified groups of cod with both directional and resident migratory behaviours. By combining information from behavioural groups identified by data storage tags, with individual genetic analyses, we were able to test specifically the hypothesis that cod return to a likely place of birth at the time of spawning. This combination of methods thus elucidates the links between different life stages, and provides new insights on how population structures in marine systems are maintained and may evolve.
When testing for overall population divergence among the 15 adult and juvenile samples, we found a low but overall significant level of divergence (8 loci: F ST = 0.0027; P = 0001: Additional file 1: Table S1). The level of divergence was fairly similar among loci (cf. Additional file 1: Table S1) and no statistical “outlying” loci were detected (Additional file 2: Figure S1), indicating absence of strong selective forces acting on these loci. Hence, they were considered as suitable for population analyses. The analysis of pairs of samples (Additional file 3: Table S2) show that the two main oceanic populations of cod, in the North Sea and Kattegat, respectively, are clearly genetically divergent (mean pairwise differentiation between the two regions, F ST ± CI0.95 = 0.0036 ± 0.0013; Additional file 4: Table S3). This divergence appears temporally stable, as evidenced by samples collected in different years (Additional file 4: Table S3a), and from comparisons of year-classes (Additional file 4: Table S3b). These findings represent key requirements for statistical assignment of behavioural groups to population of origin (below).
Another important observation from these pairwise comparisons (Additional file 3: Table S2) is that juvenile 0-group cod samples collected along the Eastern Skagerrak coast in 2005 (samples SKJ05a, b, and c) and offshore in 2011 (SKJ11 and KAJ11) were all rather similar to adult eastern North Sea cod (cf. Additional file 3: Table S2), and clustered together with them (Fig. 1).
Genetic assignment of tagged fish
Assignment of adult reference individuals collected in spawning areas showed that about 80 % (83 % for the Kattegat and 79 % for the North Sea) of the fish were assigned back to their original reference population (Table 1: upper rows). The exclusion test showed that all 100 tagged individuals had a probability of 13 % or more (median 70 %) to be encountered in at least one of the two reference populations, and hence all tagged individuals were included in further analyses. The recapture positions and the genetic assignment of the tagged fish (either to the North Sea/western Skagerrak reference or to the Kattegat reference) with five different migratory behaviours are shown in Fig. 2. For the two groups of fish that migrated towards the North Sea, 60 and 63 % of the individuals were assigned to the North Sea reference population (Table 1a). All fish (3 individuals) that migrated from the Skagerrak towards the Kattegat were assigned to the Kattegat reference. A majority (67 %) of the fish that were proposed to be resident in the Kattegat were also assigned genetically to the Kattegat reference (Table 1a). Similarly, 62 % of the resident fish in the Skagerrak were assigned to the North Sea/western Skagerrak reference. Of the 100 tagged fish, 35 individuals showed directional migration in the geolocation data. An exact test showed that the direction of migration was not independent of assignment (Table 1b; P = 0.032), and a higher proportion than random migrated towards the population they were genetically assigned to. This finding is consistent with both philopatry (spawning site fidelity) and natal philopatry (spawning site fidelity to their “birth” site) in the migrating fish.
To assess whether the unequal sizes of the North Sea (n = 201) and Kattegat (n = 435) reference samples affected assignment power, the Kattegat reference was subsampled at n = 201 ten times (equal to the North Sea sample) and the assignment of the tagged fish was accordingly repeated. The results were similar to the original pattern, although somewhat less pronounced (compare Additional file 5: Table S4 with Table 1). Hence, there was no indication of bias towards the larger Kattegat sample.
The genetic assignment of juveniles in this study provides evidence of recurrent transport of Atlantic cod, most likely at the pelagic egg and larval stages, from the North Sea to coastal Skagerrak and Kattegat, in agreement with previous findings [41, 42]. Despite the obvious potential for gene flow caused by this transport, there still remains a weak but temporally stable genetic population structure of cod in the eastern North Sea-Skagerrak-Kattegat region ([34, 46–50], this study). The association between adult individual migratory behaviour and genetic assignments reported herein suggests that natal homing may act as a governing mechanism for maintaining the cod stock structure even on a relatively small spatial scale.
The tagged fish displayed a variety of migratory and resident behaviours, as has been shown previously in Atlantic cod . The numbers of fish in the five different behavioural groups were relatively low and the results should accordingly be interpreted with caution. Nevertheless, our results show that cod with different behaviour mainly act as predicted under the hypothesis of philopatry. The majority of cod migrating towards the western Skagerrak and the North Sea were genetically assigned to the North Sea spawning population, whereas the majority of cod migrating towards the Kattegat or remaining stationary within Kattegat were assigned to the Kattegat spawning population. Cod that remain along the eastern Skagerrak during the course of this study assigned predominantly to the North Sea population, as may be hypothesised from consideration of North Sea larval drift and poor local recruitment in this area. These observations were done despite relatively low genetic differentiation between the Kattegat and the North Sea reference populations, as exemplified by the moderate self-assignment proportions for the two reference samples (cf. Table 1) and the low F ST between them (0.0041). Such levels of genetic differentiation have however been suggested to be of clear biological relevance , are typical of cod also in other regional studies [26, 35, 51], and are indicative of gene flow among populations and/or a recent common historical origin. Gene flow can arise from straying of adults among spawning areas  or from larval drift, if some settled cod enter the local population . Historically, Atlantic cod colonised the Skagerrak-Kattegat area less than 8–10,000 YBP and most populations are evolutionary young [52, 53]. Low differentiation between reference populations reduces the power of the genetic assignment tests , and is a likely cause for the low assignment proportions that we observed for the recaptured cod (typically 80 % self assignment in reference populations and around 65 % for each behaviour group: Table 1a and Additional file 5: Table S4). As a comparison, in a recent study of Atlantic salmon with an order of magnitude higher population divergence (F ST ≈ 0.07), self assignment was 65 % using 14 microsatellite loci and ten source populations (n = 50) . Hence, low assignment proportions of tagged fish should not necessarily be interpreted as evidence for high straying or low natal fidelity in the present study. Instead, the tagged fish are expected to show a somewhat less clear picture than the reference fish, which were specifically targeted at spawning sites at spawning time. As indicated by the exclusion tests, some of the tagged fish may have originated from populations outside the two reference groups, and this would reduce the assignment proportions. Other issues that could affect the results include the possibility that some of the eleven fish that were tagged during spring/summer had not yet started the spawning migration when recaptured in the autumn, and therefore were not classified in the correct behavioural group. Also, some groups consisted of a few individuals only (Table 1). Nevertheless, the proportions of fish in the five behaviour groups that were assigned to the expected population of origin was only slightly lower than the self assignment proportions of individual members in those populations, and in all five cases the majority of fish behaved in accordance with philopatry.
In salmonids, the evolutionary most important factor favouring natal homing has been suggested to be that such behaviour return locally adapted individuals to suitable habitats . For a marine fish such as the Atlantic cod, it appears that certain areas are more suitable for spawning and that these sites are used both by migratory and resident stocks [5, 26, 56–59]. The proximate causes for the natal homing behaviour in cod are still unknown, although recent genome scan studies show that certain genomic regions are associated with migratory behaviour in cod, indicating a genetic basis for migration [39, 60]. The migratory patterns revealed in the present study seem to exclude simple environmental cues, such as prevailing directions of sea currents . Experiments with salmonids have demonstrated the localization of the “natal stream” in anadromous species to be part of a learning process during the juvenile phase ; while small-scale localization rely on chemical cues, large scale orientation is suggested to involve geomagnetic imprinting [63, 64]. A recent experimental study on Sea turtles indicates that the magnetic environment during early development, already at the egg stage, can influence the magnetic orientation behaviour in subsequent life stages . Further, the mechanisms behind natal homing has been extensively studied in tropical fish and recent modelling studies indicate that recruitment to natal reefs are expected to be relatively high and depending on the combination of sensory abilities and oceanographic features [18, 66]. For Atlantic herring, it has been suggested that migration between spawning, wintering and feeding grounds is a socially transferred behaviour, where new year classes adopt the same migratory patterns as older herring cohorts . However, mixing during the juvenile stages and the subsequent migration to different spawning grounds harbouring genetically distinct populations [22, 68], suggests segregation also by natal homing in herring. For marine fishes like Atlantic herring and Atlantic cod, the stock structure may thus be imprinted already in early life stages, before juveniles of different origins intermingle in nursery areas.
Population structure in the North Sea-Skagerrak-Kattegat
The pattern emerging from this and previous studies on cod population structure in the North Sea-Skagerrak-Kattegat region can be summarised as follows: Cod eggs and larvae are transported with currents from the North Sea to coastal Skagerrak. Juvenile cod collected in the coastal Skagerrak in 2001, 2004, 2005 and 2011 were more similar genetically to North Sea spawning aggregations than to local adult cod ([34, 41, this study Fig. 1]). Here we present new evidence that larval transport also protrude into the Kattegat (Fig. 1). The strength and direction of larval drift are apparently governed by an interaction between the size and location of the North Sea cod spawning biomass and the sea-current strength from the North Sea into the Skagerrak/Kattegat during and after spawning, and has a significant effect on the abundance of juvenile cod along the coast [40, 42]. Heath et al.  showed that North Sea cod consists of two genetically distinct units. While our study, and also previous genetic and modelling, have considered the Dogger unit in the central North Sea, recent work using genetic assignment tests on juvenile cod collected in 2014 in Skagerrak and northern Kattegat suggest that both Dogger Bank and the Viking Bank may be important sources of larvae for the area (J. Hemmer Hansen unpubl. data). Nevertheless, as cod larvae are transported into coastal fjord systems, they might be trapped by oceanographic forces, and mixed together with local coastal cod larvae . On the Norwegian Skagerrak coast small local genetically distinct fjord populations persist in spite of this extensive larval drift [34, 47], indicating that foreign larvae are not incorporated into local coastal populations. The distribution pattern of juveniles and adults [36, 40, 44] and tagging studies [25, 43] suggest instead that many juveniles that grow up in coastal Skagerrak migrate back towards the North Sea when reaching maturity, at 2–3 years age; this return migration is corroborated by the genetic assignment of tagged fish presented here. In addition, we find a similar philopatric movement of cod of presumed Kattegat origin. Although conclusive evidence have been lacking, indirect evidence of natal homing have been found in several other cod populations in the eastern [24, 26, 39, 59] and western Atlantic [27, 56, 57].
The distribution of adult fish in the eastern North Sea region can be regarded both as a result of spawning aggregation in certain areas and the return migrations linked to larval drift from these reproduction units. If natal homing is an important mechanisms for cod migration, as our results suggest, an area currently depleted of its local stock components can be repopulated by fish from adjacent areas only slowly, as straying adults aggregate at new locations and give rise to new populations [13, 59, 71]. The disappearance and slow recovery of local stocks of cod in the North Sea region [28, 36, 44, 72] and in the Northwest Atlantic [73, 74] is illustrative in this sense. Even if some inshore enclaves have shown signs of recovery, adjacent areas formerly equally productive, have not been repopulated . Recognition of a possibly slow-changing behavioural stock structuring mechanism is hence of high importance for the assessment and management of marine fisheries [4, 75].
Migratory behaviour of tagged cod
Fish used for tagging were captured either in fishing pots or in trawling sets (maximum duration of 30 min). The fish were retained in tanks supplied with running seawater long enough to determine if they were in suitable condition for tagging. Typically, these were cod that could maintain buoyancy near the bottom of the tank without apparent difficulty and without external injury, such as bloodied fins or net-marks. Healthy cod were then measured to the nearest cm total length. Previous studies have suggested that, while the capture procedure may induce re-equilibration behaviour as a result of the tagging procedure , individuals are likely to return to normal behavioural patterns within 2 weeks of release. All tagged cod were larger than 37 cm and due to the exclusive use of the small archival tag Lotek Model LTD 2410 (http://www.lotek.com), the weight of the tag never exceeded 2.0 % of the fish body weight. Archival tags recording temperature, pressure and light intensity were deployed on 2–5 year old cod between 2003 and 2005 on the Swedish west coast. At release, records on length, weight, sex (if possible) and GPS positions were taken. All tagging was conducted under governmental licence and in adherence with national regulations on the treatment of experimental animals.
In all, we captured, tagged and released 417 individual cod, 142 in the Kattegat and 275 in the eastern Skagerrak. The archival tag was equipped with a real time clock set to UTC, a pressure sensor, external and internal temperature sensors (i.e., inside the capsulation), and a light intensity sensor. The tags were programmed to telescope the data; i.e., once the memory capacity was exhausted, some of the data were overwritten in a specified linear manner. This allows the total storage capacity of the tag to be evenly distributed over the entire mission, without sacrificing the temporal resolution of the data during the logging period. The retained data were saved in blocks of minimum 48 h. The tag also stored several parameter values on a daily basis during the deployment: e.g., estimated longitude and latitude, i.e., onboard processed estimates, sunrise and sunset in UTC, maximum external temperature, and maximum and minimum depths.
All fish were fin-clipped for DNA analyses and released at the location of capture. In total 162 individuals were recaptured between 2003 and 2006, mostly by commercial fishermen, out of which 100 individuals were selected for genetic analyses depending on two conditions: the period of time at liberty exceeded 30 days, and the migratory history of the individual could be unambiguously reconstructed from the archival tag information. Most of the fish finally selected were tagged in the autumn and had been in liberty over the spawning season in the winter, but some (11 out of 100) fish that where tagged in the spring and recaptured during the autumn had not.
Migratory trajectories of individual fish were obtained by retrospective inspection of recorded light intensity data: estimation of local noon (or local midnight) gives records on the longitude, whereas estimation of day length gives the latitude . Three or more consecutive estimates departing more than one longitudinal degree from the previous location was considered as a new, valid geolocated longitude. In addition, the onboard tag algorithm estimated longitude by defining dusk and dawn at civic twilight (zenith equal to 93.44°) as characteristic changes in light intensity, and by recording the times at which they occurred. Due to data storage limitations at the time of deployment, not all days at liberty at sea could be retrospectively inspected. In order to evaluate the automatic onboard geolocation estimates an extended Kalman filter-tracking model was used as for sorting out geographic signals from the tag position time series . The software package KFtrack 0.61 in the R statistical environment  was used for track estimation from day-logged positions. In this way the migration in east–west direction was reconstructed for the entire duration at liberty. These data were subsequently crosschecked and confirmed by the tidal location method , which also give estimates of the latitude. See Svedäng et al.  and Righton et al.  for more details about the tagging procedure and migration reconstruction.
The study material of 100 recaptured and genotyped (below) cod was divided into subsets of fish displaying common migratory patterns, including directional migration to the North Sea from Skagerrak or Kattegat, or to Kattegat from Skagerrak, and nonmigratory within Skagerrak or Kattegat, resulting in five behavioural categories (Table 1). Cod that were recaptured or geolocated west of 10 E° (see Fig. 3) were considered to be migrating towards the North Sea.
Adult reference and juvenile sample collection
All recaptured fish were genotyped and statistically assigned to potential populations of origin using new and existing genetic data on cod [30, 41], collected from different spawning aggregations (Table 2). In the eastern Skagerrak-Kattegat area, the present major spawning grounds and thus the putative populations of origin for the tagged fish, are located in the North Sea/western Skagerrak [36, 41, 42] and in the Kattegat and in the Öresund [82, 83]. Whilst there are no known major spawning areas in eastern Skagerrak, small local fjord populations have been reported on the Norwegian coast in the northern Skagerrak [34, 47], and possibly also on the Swedish coast . These aggregations were not, however, included in the present analysis for several reasons. The populations on the Norwegian coast are small and resident, and located mainly in the northern Skagerrak. On the Swedish coast, the aggregations are depleted and it is unclear if they represent spawning populations. Lastly, it has been shown that the power of the assignment procedure decrease with the number of putative source populations, given that no major contributing source has been left uncovered .
A total of nine samples of adult Atlantic cod representing the two reference groups were collected during 2000–2004 (global n = 636). Group 1, constitutes of spawning aggregations in the eastern North Sea, including the “Dogger unit”  and the western Skagerrak, and group 2 include the Kattegat and Öresund (Fig. 3, 1, Table 2). An additional sample of adult Öresund cod, collected in 2003, was also genotyped for comparison purposes, but was not included in the reference set (Table 2). The fish were collected during the spawning period from January to March by trawling, and care was taken to choose mature fish that were at or close to spawning. Muscle tissue or fin clips for DNA analysis were stored in ethanol.
Juvenile 0-group cod (about 6–12 cm in length) were sampled in June 2005 at three locations in the eastern coastal Skagerrak using a beach seine, and two locations in August 2011 in the Kattegat and Skagerrak by trawling during the IBTS expedition (http://www.ices.dk; Fig. 3, Table 2). All cod samples used were collected in compliance with EU Directive 2010/63/EU, and the national legislations in Sweden (Swedish Board of Agriculture (http://www.jordbruksverket.se,) permit number no. 126-2015 to Swedish University of Agricultural Sciences, Department of Aquatic Resources) and Norway (the Institute of Marine Research have permission to sample cod by the Directorate of Fisheries, Bergen, Norway). No endangered species were used in the present study.
DNA was extracted from fin clips or muscle tissue using the dneasy animal tissue kit (Qiagen Inc.). All fish were genotyped for eight microsatellite DNA loci following published protocols with minor modifications: Gmo2 and Gmo132 ; Gmo3, Gmo8, Gmo19, Gmo34, Gmo35 ; and Tch5 . To increase resolution, the tagged and reference fish were genotyped for additional four loci (i.e., in total twelve): Gmo36, and Gmo37  and Tch13 and Tch22 . The microsatellite DNA fragments were separated on ALF express II (Amersham Pharmacia Biotech), CEQ 8000 Genetic Analysis System (Beckman Coulter) and ABI 3730 Genetic Analyser (Applied Biosystems) automatic sequencers. Repeating two “control” individuals spanning the anticipated allelic ranges on all runs, in addition to internal and external size ladders, ensured scoring consistency among runs and platforms. Two persons scored genotypes independently and any inconsistent scorings were noted and the fish was screened again. We assessed scoring quality of the final data by searching for allele-specific departures from Hardy–Weinberg proportions at each locus, as could arise from allele dropout or null-alleles. This was done by considering one allele at a time, pooling all other alleles at the locus, and calculating F IS for each allele separately based on genotypes thus pooled. We used the Chi square test X2 = n*F IS  with df = 1, where n is the sample size, to assess deviations from HW proportion for that allele. However, no systematic departures were found (data not shown).
Statistical analyses of genetic data
Heterozygosity in adult and juvenile samples was calculated using FSTAT  (Table 2). Population differentiation among samples was estimated using the F ST estimator θ , and statistical significance assessed using the software genepop 4.2. . The software lositan  was used to test for outliers among the eight loci, possibly affected by selection. The pattern of genetic differentiation was explored among the adult and juvenile samples with a Multi Dimensional Scaling plot, based on pair-wise F ST (Additional file 3: Table S2) using cmndscale (r-project.org). The MDS plot (Fig. 1) further indicated that the nine reference samples could be pooled into two regional groups: “North Sea/W Skagerrak” (on the left side in Fig. 1) and “Kattegat” (right side); see also table of pair-wise F ST (Additional file 4: Table S3a). The pooling of reference samples had the benefit of both reducing the number of reference populations and increasing the reference sample sizes without losing significant geographic information. We further assessed the temporal stability of the genetic differentiation between the two regional groups by dividing all fish into year classes based on otolith ageing, and testing for genetic heterogeneity (Additional file 4: Table S3b). This latter test did not indicate any temporal instability in either the North Sea/W Skagerrak or the Kattegat groups and verified the genetic integrity of the two regions. All reference samples were thus pooled into two regions North Sea/W Skagerrak and Kattegat in subsequent statistical analyses; the F ST between the pooled North Sea/W Skagerrak (n = 201) and the pooled Kattegat reference (n = 435) was 0.0041 (P < 0.0001).
Genetic assignment of recaptured tagged individuals to the two pooled reference populations, as well as “self assignment” of reference individuals, was performed with the software geneclass2 . The possibility that the tagged fish originated from unknown spawning populations, not covered in the reference data, was evaluated using the exclusion method in geneclass2, where the likelihood of individual fish belonging to a given reference population was compared with the distribution of likelihoods of 1000 genotypes simulated from each reference population with a Monte Carlo algorithm .
The most likely origin of fish in the five behavioural groups was determined by calculating the proportion of fish that were assigned to each of the two regional reference pools. For the 35 fish that showed directional migration we tested if migration was independent from genetic assignment using an unconditional exact test of independence (Boschloo’s test in the exact R-package, ).
Availability of supporting data
Hauser L, Carvalho GR. Paradigm shifts in marine fisheries genetics: ugly hypotheses slain by beautiful facts. Fish Fish. 2008;9:333–62.
Reiss H, Hoarau G, Dickey-Collas M, Wolff WJ. Genetic population structure of marine fish: mismatch between biological and fisheries management units. Fish Fish. 2009;10:361–95.
Schindler DE, Hilborn R, Chasco B, Boatright CP, Quinn TP, Rogers L, Webster MS. Population diversity and the portfolio effect in an exploited species. Nature. 2010;465:609–12.
Ciannelli L, Fisher JAD, Mauritzen MS, Hunsicker ME, Hidalgo M, Frank KT, Bailey KM. Theory, consequences and evidence of eroding population spatial structure in harvested marine fishes: a review. Mar Ecol Prog Ser. 2013;480:227–43.
Robichaud D, Rose GA. Migratory behaviour and range in Atlantic cod: inference from a century of tagging. Fish Fish. 2004;5:185–214.
Thorrold SR, Latkoczy C, Swart PK, Jones CM. Natal homing in a marine fish metapopulation. Science. 2001;291:297–9.
Cowen RK, Sponaugle S. Larval dispersal and marine population connectivity. Ann Rev Mar Sci. 2009;1:443–66.
Ciannelli L, Knutsen H, Olsen EM, Espeland SH, Asplin L, Jelmert A, Knutsen JA, Stenseth NC. Maintenance of small-scale genetic structure in a marine population in relation to water circulation and egg characteristics. Ecology. 2010;91:2918–30.
Pacariz S, Björk G, Jonsson P, Börjesson P, Svedäng H. A model study of the large scale transport of fish eggs in the Kattegat in relation to egg density. ICES J Mar Sci. 2014;71:345–55.
Myksvoll MS, Jung K-M, Albretsen J, Sundby S. Modelling dispersal of eggs and quantifying connectivity among Norwegian coastal cod subpopulations. ICES J Mar Sci. 2014;71:957–69.
Selkoe KA, Watson JR, White C, Ben-Horin T, Iacchei M, Mitarai S, Siegel DA, Gaines SD, Toonen RJ. Taking the chaos out of genetic patchiness: seascape genetics reveals ecological and oceanographic drivers of genetic patterns in three temperate reef species. Mol Ecol. 2010;19:3708–26.
Bradbury IR, Laurel BJ. Defining ‘natal homing’ in marine fish populations. Mar Ecol Prog Ser. 2007;349:307–8.
Heath MR, Kunzlik PA, Gallego A, Holmes SJ, Wright PJ. A model of meta-population dynamics for North Sea and West of Scotland cod - The dynamic consequences of natal fidelity. Fish Res. 2008;93:92–116.
Frisk MG, Jordaan A, Miller TJ. Moving beyond the current paradigm in marine population connectivity: are adults the missing link? Fish Fish. 2014;15:242–54.
Almany GR, Berumen ML, Thorrold SR, Planes S, Jones GP. Local replenishment of coral reef fish populations in a marine reserve. Science. 2007;316:742–4.
Gerlach G, Atema J, Kingsford MJ, Black KP, Miller-Sims V. Smelling home can prevent dispersal of reef fish larvae. Proc Nat Acad Sci. 2007;104:858–63.
Saenz-Agudelo P, Jones GP, Thorrold SR, Planes S. Patterns and persistence of larval retention and connectivity in a marine fish metapopulation. Mol Ecol. 2012;21:4695–705.
Mouritsen H, Atema J, Kingsford MJ, Gerlash G. Sun compass orientation helps coral reef fish larvae return to their natal reef. PLoS One. 2013;8:e66039. doi:10.1371/journal.pone.0066039.
Hendry AP, Castric V, Kinnison MT, Quinn TP. The evolution of philopatry and dispersal - homing versus straying in salmonids. In: Hendry AP, Stearns SC, editors. Evolution Illuminated: Salmon and their relatives. Oxford: Oxford University Press; 2004.
Rooker JR, Secor DH, DeMetrio G, Schloesser R, Block BA, Neilson JD. Natal homing and connectivity in Atlantic bluefin tuna populations. Science. 2008;322:742–4.
Righton D, Aarestrup K, Jellyman D, Sébert P, van den Thillart G, Tsukamoto K. The Anguilla spp. migration problem: 40 million years of evolution and two millennia of speculation. J Fish Biol. 2012;81:365–86.
Gaggiotti OE, Bekkevold D, Jørgensen HBH, Foll M, Carvalho GR, André C, Ruzzante DE. Disentangling the effects of evolutionary, demographic, and environmental factors influencing genetic structure of natural populations: Atlantic herring as a case study. Evolution. 2009;63:2939–51.
Engstedt O, Engkvist R, Larsson P. Elemental fingerprinting in otoliths reveals natal homing of anadromous Baltic Sea pike (Esox lucius L.). Ecol Fresh Fish. 2013;23:313–21.
Storr-Paulsen M, Wieland K, Hovgaard H, Rätz HJ. Stock structure of Atlantic cod (Gadus morhua) in West Greenland waters: implications of transport and migration. ICES J Mar Res. 2004;61:972–82.
Svedäng H, Righton D, Jonsson P. Migratory behaviour of Atlantic cod Gadus morhua: natal homing is the prime stock separating mechanism. Mar Ecol Prog Ser. 2007;345:1–12.
Nordeide JT, Johansen SD, Jorgensen TE, Karlsen BO, Moum T. Population connectivity among migratory and stationary cod Gadus morhua in the Northeast Atlantic - A review of 80 years of study. Mar Ecol Prog Ser. 2011;435:269–83.
Morris CJ, Green JM, Snelgrove PVR, Pennell CJ, Ollerhead LMN. Temporal and spatial migration of Atlantic cod (Gadus morhua) inside and outside a marine protected area and evidence for the role of prior experience in homing. Can J Fish Aquat Sci. 2014;71:1704–12.
Wright PJ, Neat FC, Gibb FM, Gibb IM, Thordarson H. Evidence for metapopulation structuring in cod from the west of Scotland and North Sea. J Fish Biol. 2006;69:181–99.
Jones GP, Almany GR, Russ GR, Sale PF, Steneck RS, van Oppen MJH, Willis BL. Larval retention and connectivity among populations of corals and reef fishes: history, advances and challenges. Coral Reefs. 2009;28:307–25.
Svedäng H, André C, Jonsson P, Elfman M, Limburg KE. Homing behaviour and otolith chemistry suggest fine-scale sub-population structure within a genetically homogenous Atlantic cod population. Env Biol Fish. 2010;89:383–97.
Jones GP, Planes S, Thorrold SR. Coral fish larvae settle close to home. Curr Biol. 2005;15:1314–8.
O’Leary DB, Coughlan J, Dillane E, McCarthy TV, Cross TF. Microsatellite variation in cod Gadus morhua through-out its geographic range. J Fish Biol. 2007;70:310–35.
Berg P, Jentoft S, Staar B, Ring KH, Knutsen H, Lien S, Jakobsen K, Andre C. Adaptation to low salinity promotes genomic divergence in Atlantic cod (Gadus morhua L.). Genome Biol Evol. 2015;7:1644–63.
Knutsen H, Olsen EM, Jorde PE, Espeland SH, André C, Stenseth NC. Are low but statistically significant levels of genetic differentiation in marine fishes “biologically meaningful”? A case study of coastal Atlantic cod. Mol Ecol. 2011;20:768–83.
Nielsen EE, Wright P, Wright PJ, Hemmer-Hansen J, Poulsen NA, Gibb IM, Meldrup D. Microgeographical population structure of cod Gadus morhua in the North Sea and west of Scotland: the role of sampling loci and individuals. Mar Ecol Prog Ser. 2009;376:213–25.
Svedäng H. The inshore demersal fish community on the Swedish Skagerrak coast: regulation by recruitment from offshore sources. ICES J Mar Res. 2003;60:23–31.
ICES. Report of the working group on the assessment of demersal stocks in the North Sea and Skagerrak (WGNSSK). ICES CM 2012/ACOM:13. 2012a; p.1204.
ICES. Report of the Baltic Fisheries Assessment Working Group 2012 (WGBFAS). ICES CM 2012/ACOM:10. 2012b; p.859.
Hemmer-Hansen J, Nielsen EE, Therkildsen NO, Taylor M, Ogden R, Geffen AJ, Bekkevold D, Helyar S, Pampoulie C, Johansen T. FISHPOPTRACE CONSORTIUM, Carvalho G. A genomic island linked to ecotype divergence in Atlantic cod. Mol Ecol. 2013;22:2653–67.
Cardinale M, Svedäng H. Modelling recruitment and abundance of Atlantic cod, (Gadus morhua), in the eastern Skagerrak-Kattegat (North Sea): evidence of severe depletion due to a prolonged period of high fishing pressure. Fish Res. 2004;69:263–82.
Knutsen H, André C, Jorde PE, Skogen MD, Thuróczy E, Stenseth NC. Transport of North Sea cod larvae into the Skagerrak coastal populations. Proc Roy Soc Biol. 2004;271:1337–44.
Stenseth NC, Jorde PE, Chan KS, Knutsen H, André C, Skogen M, Lekve K. Ecological and genetic impact of Atlantic cod larval drift in the Skagerrak. Proc Roy Soc Biol. 2006;273:1085–92.
Pihl L, Ulmestrand M. Migration pattern of juvenile cod (Gadus morhua) on the Swedish west coast. ICES J Mar Sci. 1993;50:63–70.
Svedäng H, Svenson A. Populations as behavioural units: inference from time series on juvenile cod abundance in the Skagerrak. J Fish Biol. 2006;69:151–64.
Svedäng H, Righton D, Jonsson P. Defining ‘natal homing’ in marine fish populations; need for inference in fishery science: reply to Bradbury & Laurel (2007). Mar Ecol Prog Ser. 2007;347:309–10.
Nielsen EE, Grønskjær P, Meldrup D, Paulsen H. Retention of juveniles within a hybrid zone between North Sea and Baltic Sea Atlantic cod (Gadus morhua). Can J Fish Aquat Sci. 2005;62:2219–25.
Jorde PE, Knutsen H, Espeland SH, Stenseth NC. Spatial scale of genetic structuring in coastal cod Gadus morhua and geographic extent of local populations. Mar Ecol Prog Ser. 2007;343:229–37.
Olsen EM, Knutsen H, Gjøsæter J, Jorde PE, Knutsen JA, Stenseth NC. Small-scale biocomplexity in coastal Atlantic cod supporting a Darwinian perspective on fisheries management. Evol Appl. 2008;1:524–33.
Knutsen H, Jorde PE, André C, Stenseth NC. Fine-scaled geographical population structuring in a highly mobile marine species: the Atlantic cod. Mol Ecol. 2003;12:385–94.
Poulsen NA, Hemmer-Hansen J, Loeschcke V, Carvalho GR, Nielsen EE. Microgeographical population structure and adaptation in Atlantic cod Gadus morhua: spatio-temporal insights from gene-associated DNA markers. Mar Ecol Prog Ser. 2011;436:231–43.
Ruzzante D, Taggart CT, Cook D. A nuclear DNA basis for shelf- and bank-scale population structure in northwest Atlantic cod (Gadus morhua): Labrador to Georges Bank. Mol Ecol. 1998;7:1663–80.
Gyllencreutz R, Backman J, Jakobsson M, Kissel C, Arnold E. Post glacial paleo-oceanography in the Skagerrak. Holocene. 2006;16:975–85.
Bigg GR, Cunningham CW, Ottersen G, Pogson GH, Wadley MR, Williamson P. Ice-age survival of Atlantic cod: agreement between palaeoecology models and genetics. Proc Roy Soc Biol. 2008;275:163–73.
Hansen M, Kenchington E, Nielsen EE. Assigning individual fish to populations using microsatellite DNA markers: methods and applications. Fish Fish. 2001;2:93–112.
Glover KA, Hansen MM, Lien S, Als TD, Høyheim B, Skaala Ø. A comparison of SNP and STR loci for delineating population structure and performing individual genetic assignment. BMC Genet. 2010;11:2.
Bradbury IR, Laurel BJ, Robichaud D, Rose GA, Snelgrove PVR, Gregory RS, Cote D. Windle MJS. Discrete spatial dynamics in a marine broadcast spawner: re-evaluating scales of connectivity and habitat associations in Atlantic cod (Gadus morhua L.) in coastal Newfoundland. Fish Res. 2008;2–3:299–309.
Howell HW, Morin M, Rennels N, Goethel D. Residency of adult Atlantic cod (Gadus morhua) in the western Gulf of Maine. Fish Res. 2008;91:123–32.
Skjaeraasen JE, Meager JJ, Karlsen O, Hutchings JA, Fernö A. Extreme spawning-site fidelity in Atlantic cod. ICES J Mar Sci. 2011;68:1472–7.
Neuenfeldt S, Righton D, Neat F, Wright PJ, Svedäng H, Michalsen K, Subbey S, Steingrund P, Thorsteinsson V, Pampoulie C, Andersen KH, Pedersen MW, Metcalfe J. Analysing migrations of Atlantic cod Gadus morhua in the north-east Atlantic Ocean: then, now and the future. J Fish Biol. 2013;82:741–63.
Karlsen BO, Klingan K, Emblem Å, Jørgensen TE, Jueterbock A, Furmanek T, Hoarau G, Johansen S, Nordeide JT, Moum T. Genomic divergence between the migratory and stationary ecotypes of Atlantic cod. Mol Ecol. 2013;22:5098–111.
Harden-Jones FR. Fish Migrations. London: Arnold; 1968.
Dittman AH, Quinn TP. Pacific salmon homing: mechanisms and ecological basis. J Exp Biol. 1996;199:83–91.
Putman NF, Lohmann KJ, Putman EM, Quinn TP, Klimley AP, Noakes DLG. Evidence for Geomagnetic Imprinting as a Homing Mechanism in Pacific Salmon. Curr Biol. 2013;23:312–6.
Putman NF, Scanian MM, Billman EJ, O’Neill JP, Couture RB, Quinn TP, Lohman KJ, Noakes DLG. An inherited magnetic map guides ocean navigation in juvenile Pacific salmon. Curr Biol. 2014;24:446–50.
Fuxjager MJ, Davidoff KR, Mangiamele LA, Lohmann KJ. The geomagnetic environment in which sea turtle eggs incubate affects subsequent magnetic navigation behaviour of hatchlings. Proc R Soc B. 2014;281:20141218. doi:10.1098/rspb.2014.1218.
Wolanski E, Kingsford MJ. Oceanographic and behavioural assumptions in models of the fate of coral and coral reef fish larvae. J R Soc Interface. 2014;11:20140209. doi:10.1098/rsif.2014.0209.
Corten A. The role of “conservatism” in herring migrations. Rev Fish Biol Fish. 2002;11:339–61.
Ruzzante DE, Mariani S, Bekkevold D, André C, Mosegaard H, Clausen LAW, Dahlgren T, Hutchinson W, Hatfield E, Torstensen E, Brigham J, Simmonds J, Laikre L, Larsson L, Stet R, Ryman N, Carvalho G. Biocomplexity in a highly migratory marine pelagic fish. Proc Roy Soc Biol. 2006;273:1459–64.
Heath MR, Culling MA, Crozier WW, Fox CJ, Gurney WSC, Hutchinson WF, Nielsen EE, O’Sullivan M, Preedy KF, Righton DA, Speirs DC, Taylor MI, Wright PJ, Carvalho GR. Combination of genetics and spatial modelling highlights the sensitivity of cod (Gadus morhua) population diversity in the North Sea to distributions of fishing. ICES J Mar Sci. 2014;71:794–807.
Knutsen H, Olsen EM, Cianelli L, Espeland SH, Knutsen JA, Simonsen SS, Stenseth NC. Egg distribution bottom topography and small-scale population structure in a coastal marine system. Mar Ecol Prog Ser. 2007;333:249–55.
Rose GA, Nelson RJ, Mello LGS. Isolation or metapopulation: whence and whither the Smith Sound cod? Can J Fish Aquat Sci. 2011;68:152–69.
Wright PJ, Millar CP, Gibb FM. Intrastock differences in maturation schedules of Atlantic cod. Gadus morhua. ICES J Mar Sci. 2011;68:1918–27.
Hu L. Wroblewski J Conserving a subpopulation of the northern Atlantic cod metapopulation with a marine protected area. Aquatic Conservation: Mar Fresh Ecosyst. 2009;19:178–93.
Hutchings JA, Rangeley RW. Correlates of recovery for Canadian Atlantic cod (Gadus morhua). Can J Zool. 2011;89:386–400.
Petitgas P, Secor DH, McQuinn I, Huse G, Lo N. Stock collapses and their recovery: mechanisms that establish and maintain life-cycle closure in space and time. ICES J Mar Sci. 2010;67:1841–8.
van der Kooij J, Righton D, Strand E, Michalsen K, Thorsteinsson V, Svedäng H, Neat FC, Neuenfeldt S. Life under pressure: insights from electronic data-storage tags into cod swimbladder function. ICES J Mar Sc. 2007;64:1293–301.
Musyl MK, Brill RW, Curran DS, Gunn JS, Hartog JR, Hill RD, Welch DW, Eveson JP, Boggs CH. Brainard RE. In: Sibert JR, Nielsen JL, editors. Electronic tagging and tracking in marine fisheries reviews: methods and technologies in fish biology and fisheries. Dordrecht: Kluwer Academic Press; 2001. p. 343–68.
Sibert J, Musyl MK, Brill RW. Horizontal movements of bigeye tuna (Thunnus obesus) near Hawaii determined by Kalman filter analysis of archival tagging data. Fish Oceanogr. 2003;12:141–51.
R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0. 2012;http://www.R-project.org/.
Hunter E, Aldridge JN, Metcalfe JD, Arnold GP. Geolocation of free-ranging fish on the European continental shelf as determined from environmental variables. Mar Biol. 2003;142:601–9.
Righton DA, Andersen KH, Neat F, Thorsteinsson V, Steingrund P, Svedäng H, Michalsen K, Hinrichsen H-H, Bendall V, Neuenfeldt S, Wright P, Jonsson P, Huse G, van der Kooij J, Mosegaard H, Hussy K, Metcalfe J. Thermal niche of Atlantic cod, Gadus morhua, limits, tolerance and optima. Mar Ecol Prog Ser. 2010;420:1–13.
Vitale F, Börjesson P, Svedäng H, Casini M. The spatial distribution of cod (Gadus morhua L.) spawning grounds in the Kattegat, eastern North Sea. Fish Res. 2008;90:36–44.
Börjesson P, Jonsson P, Pacariz S, Björk G, Taylor MI, Svedäng H. Spawning of Kattegat cod (Gadus morhua)– mapping spatial distribution by egg surveys. Fish Res. 2013;147:63–71.
Svedäng H, Hagberg J, Börjesson P, Svenson A, Vitale F Bottenfisk i Västerhavet. Fyra studier av beståndens status, utveckling och lekområden vid den svenska västkusten. (Demersal fish in the eastern North Sea. Four studies on stock status, development and spawning areas along the Swedish west coast). Finfo 2004;6. In Swedish with English Summary.
Brooker AL, Cook D, Bentzen P, Wright JM, Doyle RW. Organisation of microsatellites differs between mammals and cold-water teleost fishes. Can J Fish Aquat Sci. 1994;51:1959–66.
Miller KM, Le KD, Beacham TD. Development of tri- and tetranucleotide repeat microsatellite loci in Atlantic cod (Gadus morhua). Mol Ecol. 2000;9:238–9.
O’Reilly PT, Canino MF, Bailey KM, Bentzen P. Isolation of twenty low stutter di- and tetranucleotide microsatellites for population analyses of walleye Pollock and other gadoids. J Fish Biol. 2000;56:1074–86.
Nei M. Molecular Evolutionary Genetics. New York: Columbia University Press; 1987. p. p156.
Goudet J. 2001 http://www.unil.ch/izea/software/fstat.html.
Weir BS. Cockerham CC Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.
Rousset F. Genepop’007: a complete reimplementation of the Genepop software for Windows and Linux. Mol Ecol Res. 2008;8:103–6.
Antao T, Lopes A, Lopes RJ, Beja-Pereira A, Luikart G. LOSITAN: a work bench to detect molecular adaptation based on a F ST -outlier method. BMC Bioinform. 2008;9:323.
Piry SA, Alapetite A, Cornuet JM, Paetkau D, Baudouin L, Estoup A. GeneClass2: a software for genetic assignment and first-generation migrant detection. J Heredity. 2004;95:536–9.
Paetkau D, Slade R, Burden M, Estoup A. Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power. Mol Ecol. 2004;13:55–65.
Calhoun P Exact: Unconditional exact test. R package version 1.4. 2013. http://CRAN.R-project.org/package=Exact.
CA and HS conceived the study, in discussions with HK, AKR and PEJ. HS, PJ and MS performed the tagging and analysed the tagging data. CA, HK, AKR and GD did the genetic analyses. CA and PEJ analysed the genetic data. CA wrote the manuscript in collaboration with all authors. All authors read and approved the final manuscript.
We thank the many fishermen who have reported recaptures and sent back tags and fish from all over the eastern North Sea region, and Anders Svenson, Peter Jacobsson, Ann-Marie Bratt and the crews on R/V Ancylus for their assistance in the tagging practices. Benno Jönsson, Erika Norlinder and Mikael Dahl assisted with genetic analyses. Per Jonsson, Susanne Liljenström, Vidar Øresland and three reviewers gave valuable comments on an earlier version of the manuscript. This project was funded by the Swedish Research Council FORMAS, the GU Centre for marine evolutionary biology (Cemeb), the Swedish Agency for Marine and Water Management (SwAM) agreement 4821-13, the European Fisheries Fund, the Norwegian Research Council and EU Interreg funds.
All authors declare that they have no competing interests.
Additional file 2: Figure S1. Comparison of population differentiation (F ST) and heterozygosity (H E) among 15 samples of cod in eight microsatellite loci to identify outliers and potential candidates for selection using LOSITAN. All loci are candidates to be selectively neutral.
About this article
Cite this article
André, C., Svedäng, H., Knutsen, H. et al. Population structure in Atlantic cod in the eastern North Sea-Skagerrak-Kattegat: early life stage dispersal and adult migration. BMC Res Notes 9, 63 (2016). https://doi.org/10.1186/s13104-016-1878-9