Genetic analyses verify sexually mature escaped farmed Atlantic cod and farmed cod eggs in the natural environment.

Elucidating the effects of domesticated organisms escaping into the natural environment represents a topic of importance in both evolutionary and conservation biology. However, when excluding the abundant data on salmonids, there is a lack of knowledge on this topic for marine fish aquaculture, which continues to expand globally. In order to bridge this empirical gap, we investigated a suspected escape of sexually mature domesticated Atlantic cod from a commercial marine fish farm in northern Norway. This involved genotyping samples of fish from cages on the farm, putatively identified escapees and wild cod captured in the region and samples of recently spawned eggs collected in the sea. Genetic analyses confirmed a farmed ancestry of the suspected escapees, and significantly, 27% of the sampled cod eggs. Furthermore, statistical analyses revealed a strong reduction in genetic variation in all samples of the farmed cod, including low effective population size and high degree of siblingship. These results thus document the escape of sexually mature adult cod and the release of fertilized domesticated cod eggs into the natural environment. Although it is possible that some of the mature escapees spawned post-escape, the fact that only a single egg of potential hybrid farmed × wild origin was identified, together with the high number of mature cod in the farm, points to within cage spawning as the primary source of these eggs. This suggestion is supported by oceanic particle-drift modelling, verifying that transport of eggs between the farm and the egg sampling locations was plausible. This study represents a rare documentation of interaction between domesticated and wild populations for a marine fish, pointing towards potential impacts on the local wild population.

Contrasting management regimes indicative of mesopredator release in temperate coastal fish assemblages.

The absence of functional top predators has been proposed as a mechanism acting to shape fish assemblages in temperate marine ecosystems, with cascading effects on lower trophic levels. We explore this scenario by comparing the trophic and functional status of fish assemblages in Norwegian marine national parks, open to fishing, to a nearby coastal seascape that harbors a system of marine protected areas (MPAs) including a no-take zone. Demersal fish assemblages were sampled using fyke nets over three consecutive seasons. Atlantic cod (Gadus morhua) is potentially a dominant top predator in this ecosystem, and historically, this and other gadids have been targeted by the full range of former and present fisheries. In the present study, we find that average body size of the Atlantic cod was significantly larger in the zoned seascape compared to the unprotected areas (mean ± SD: 36.6 cm ± 14.38 vs. 23.4 ± 7.50; p < .001) and that the unprotected seascape was characterized by a higher abundance of mesopredator fish species. These observations are consistent with the hypothesis that the protection of top predators within MPAs aids to control the mesopredator populations and provides empirical support to the notion that the present state of many coastal fish assemblages is driven by mesopredator release linked to functional depletion of large top predators.

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.

Combining population genomics with demographic analyses highlights habitat patchiness and larval dispersal as determinants of connectivity in coastal fish species.

Gene flow shapes spatial genetic structure and the potential for local adaptation. Among marine animals with nonmigratory adults, the presence or absence of a pelagic larval stage is thought to be a key determinant in shaping gene flow and the genetic structure of populations. In addition, the spatial distribution of suitable habitats is expected to influence the distribution of biological populations and their connectivity patterns. We used whole genome sequencing to study demographic history and reduced representation (double-digest restriction associated DNA) sequencing data to analyse spatial genetic structure in broadnosed pipefish (Syngnathus typhle). Its main habitat is eelgrass beds, which are patchily distributed along the study area in southern Norway. Demographic connectivity among populations was inferred from long-term (~30-year) population counts that uncovered a rapid decline in spatial correlations in abundance with distance as short as ~2 km. These findings were contrasted with data for two other fish species that have a pelagic larval stage (corkwing wrasse, Symphodus melops; black goby, Gobius niger). For these latter species, we found wider spatial scales of connectivity and weaker genetic isolation-by-distance patterns, except where both species experienced a strong barrier to gene flow, seemingly due to lack of suitable habitat. Our findings verify expectations that a fragmented habitat and absence of a pelagic larval stage promote genetic structure, while presence of a pelagic larvae stage increases demographic connectivity and gene flow, except perhaps over extensive habitat gaps.

Stabilizing selection on Atlantic cod supergenes through a millennium of extensive exploitation.

Life on Earth has been characterized by recurring cycles of ecological stasis and disruption, relating biological eras to geological and climatic transitions through the history of our planet. Due to the increasing degree of ecological abruption caused by human influences many advocate that we now have entered the geological era of the Anthropocene, or "the age of man." Considering the ongoing mass extinction and ecosystem reshuffling observed worldwide, a better understanding of the drivers of ecological stasis will be a requisite for identifying routes of intervention and mitigation. Ecosystem stability may rely on one or a few keystone species, and the loss of such species could potentially have detrimental effects. The Atlantic cod (Gadus morhua) has historically been highly abundant and is considered a keystone species in ecosystems of the northern Atlantic Ocean. Collapses of cod stocks have been observed on both sides of the Atlantic and reported to have detrimental effects that include vast ecosystem reshuffling. By whole-genome resequencing we demonstrate that stabilizing selection maintains three extensive "supergenes" in Atlantic cod, linking these genes to species persistence and ecological stasis. Genomic inference of historic effective population sizes shows continued declines for cod in the North Sea-Skagerrak-Kattegat system through the past millennia, consistent with an early onset of the marine Anthropocene through industrialization and commercialization of fisheries throughout the medieval period.

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.

Demographic history has shaped the strongly differentiated corkwing wrasse populations in Northern Europe.

Understanding the biological processes involved in genetic differentiation and divergence between populations within species is a pivotal aim in evolutionary biology. One particular phenomenon that requires clarification is the maintenance of genetic barriers despite the high potential for gene flow in the marine environment. Such patterns have been attributed to limited dispersal or local adaptation, and to a lesser extent to the demographic history of the species. The corkwing wrasse (Symphodus melops) is an example of a marine fish species where regions of particular strong divergence are observed. One such genetic break occurred at a surprisingly small spatial scale (FST ~0.1), over a short coastline (<60 km) in the North Sea-Skagerrak transition area in southwestern Norway. Here, we investigate the observed divergence and purported reproductive isolation using genome resequencing. Our results suggest that historical events during the post-glacial recolonization route can explain the present population structure of the corkwing wrasse in the northeast Atlantic. While the divergence across the break is strong, we detected ongoing gene flow between populations over the break suggesting recent contact or negative selection against hybrids. Moreover, we found few outlier loci and no clear genomic regions potentially being under selection. We concluded that neutral processes and random genetic drift e.g., due to founder events during colonization have shaped the population structure in this species in Northern Europe. Our findings underline the need to take into account the demographic process in studies of divergence processes.

Interbreeding between local and translocated populations of a cleaner fish in an experimental mesocosm predicts risk of disrupted local adaptation.

Translocation of organisms within or outside its native range carries the risk of modifying the community of the recipient ecosystems and induces gene flow between locally adapted populations or closely related species. In this study, we evaluated the genetic consequences of large-scale translocation of cleaner wrasses that has become a common practice within the salmon aquaculture industry in northern Europe to combat sea lice infestation. A major concern with this practice is the potential for hybridization of escaped organisms with the local, recipient wrasse population, and thus potentially introduce exogenous alleles and breaking down coadapted gene complexes in local populations. We investigated the potential threat for such genetic introgressions in a large seminatural mesocosm basin. The experimental setting represented a simulated translocation of corkwing wrasse (Symphodus melops) that occurs on a large scale in the Norwegian salmon industry. Parentage assignment analysis of mesocosm's offspring revealed 30% (195 out of 651 offspring) interbreeding between the two populations, despite their being genetically (F ST = 0.094, p < 0.05) and phenotypically differentiated. Moreover, our results suggest that reproductive fitness of the translocated western population doubled that of the local southern population. Our results confirm that human translocations may overcome the impediments imposed by natural habitat discontinuities and urge for immediate action to manage the genetic resources of these small benthic wrasses.

Stable coexistence of genetically divergent Atlantic cod ecotypes at multiple spatial scales.

Coexistence in the same habitat of closely related yet genetically different populations is a phenomenon that challenges our understanding of local population structure and adaptation. Identifying the underlying mechanisms for such coexistence can yield new insight into adaptive evolution, diversification and the potential for organisms to adapt and persist in response to a changing environment. Recent studies have documented cryptic, sympatric populations of Atlantic cod (Gadus morhua) in coastal areas. We analysed genetic origin of 6,483 individual cod sampled annually over 14 years from 125 locations along the Norwegian Skagerrak coast and document stable coexistence of two genetically divergent Atlantic cod ecotypes throughout the study area and study period. A "fjord" ecotype dominated in numbers deep inside fjords while a "North Sea" ecotype was the only type found in offshore North Sea. Both ecotypes coexisted in similar proportions throughout coastal habitats at all spatial scales. The size-at-age of the North Sea ecotype on average exceeded that of the fjord ecotype by 20% in length and 80% in weight across all habitats. Different growth and size among individuals of the two types might be one of several ecologically significant variables that allow for stable coexistence of closely related populations within the same habitat. Management plans, biodiversity initiatives and other mitigation strategies that do not account for the mixture of species ecotypes are unlikely to meet objectives related to the sustainability of fish and fisheries.

Are we underestimating the occurrence of sympatric populations?

Sympatric populations are conspecific populations that coexist spatially. They are of interest in evolutionary biology by representing the potential first steps of sympatric speciation and are important to identify and monitor in conservation management. Reviewing the literature pertaining to sympatric populations, we find that most cases of sympatry appear coupled to phenotypic divergence, implying ease of detection. In comparison, phenotypically cryptic, sympatric populations seem rarely documented. We explore the statistical power for detecting population mixtures from genetic marker data, using commonly applied tests for heterozygote deficiency (i.e., Wahlund effect) and the structure software, through computer simulations. We find that both tests are efficient at detecting population mixture only when genetic differentiation is high, sample size and number of genetic markers are reasonable and the sympatric populations happen to occur in similar proportions in the sample. We present an approximate expression based on these experimental factors for the lower limit of FST , beyond which power for structure collapses and only the heterozygote-deficiency tests retain some, although low, power. The findings suggest that cases of cryptic sympatry may have passed unnoticed in population genetic screenings using number of loci typical of the pre-genomics era. Hence, cryptic sympatric populations may be more common than hitherto thought, and we urge more attention being diverted to their detection and characterization.

Inferring genetic connectivity in real populations, exemplified by coastal and oceanic Atlantic cod.

Genetic data are commonly used to estimate connectivity between putative populations, but translating them to demographic dispersal rates is complicated. Theoretical equations that infer a migration rate based on the genetic estimator FST , such as Wright's equation, FST ≈ 1/(4Nem + 1), make assumptions that do not apply to most real populations. How complexities inherent to real populations affect migration was exemplified by Atlantic cod in the North Sea and Skagerrak and was examined within an age-structured model that incorporated genetic markers. Migration was determined under various scenarios by varying the number of simulated migrants until the mean simulated level of genetic differentiation matched a fixed level of genetic differentiation equal to empirical estimates. Parameters that decreased the Ne /Nt ratio (where Ne is the effective and Nt is the total population size), such as high fishing mortality and high fishing gear selectivity, increased the number of migrants required to achieve empirical levels of genetic differentiation. Higher maturity-at-age and lower selectivity increased Ne /Nt and decreased migration when genetic differentiation was fixed. Changes in natural mortality, fishing gear selectivity, and maturity-at-age within expected limits had a moderate effect on migration when genetic differentiation was held constant. Changes in population size had the greatest effect on the number of migrants to achieve fixed levels of FST , particularly when genetic differentiation was low, FST ≈ 10-3 Highly variable migration patterns, compared with constant migration, resulted in higher variance in genetic differentiation and higher extreme values. Results are compared with and provide insight into the use of theoretical equations to estimate migration among real populations.

Can we rely on selected genetic markers for population identification? Evidence from coastal Atlantic cod.

The use of genetic markers under putative selection in population studies carries the potential for erroneous identification of populations and misassignment of individuals to population of origin. Selected markers are nevertheless attractive, especially in marine organisms that are characterized by weak population structure at neutral loci. Highly fecund species may tolerate the cost of strong selective mortality during early life stages, potentially leading to a shift in offspring genotypes away from the parental proportions. In Atlantic cod, recent genetic studies have uncovered different genotype clusters apparently representing phenotypically cryptic populations that coexist in coastal waters. Here, we tested if a high-graded SNP panel specifically designed to classify individual cod to population of origin may be unreliable because of natural selection acting on the SNPs or their linked background. Temporal samples of cod were collected from two fjords, starting at the earliest life stage (pelagic eggs) and carried on until late autumn (bottom-settled juveniles), covering the period during summer of high natural mortality. Despite the potential for selective mortality during the study period, we found no evidence for selection, as both cod types occurred throughout the season, already in the earliest egg samples, and there was no evidence for a shift during the season in the proportions of one or the other type. We conclude that high-graded marker panels under putative natural selection represent a valid and useful tool for identifying biological population structure in this highly fecund species and presumably in others.

Genome architecture enables local adaptation of Atlantic cod despite high connectivity.

Adaptation to local conditions is a fundamental process in evolution; however, mechanisms maintaining local adaptation despite high gene flow are still poorly understood. Marine ecosystems provide a wide array of diverse habitats that frequently promote ecological adaptation even in species characterized by strong levels of gene flow. As one example, populations of the marine fish Atlantic cod (Gadus morhua) are highly connected due to immense dispersal capabilities but nevertheless show local adaptation in several key traits. By combining population genomic analyses based on 12K single nucleotide polymorphisms with larval dispersal patterns inferred using a biophysical ocean model, we show that Atlantic cod individuals residing in sheltered estuarine habitats of Scandinavian fjords mainly belong to offshore oceanic populations with considerable connectivity between these diverse ecosystems. Nevertheless, we also find evidence for discrete fjord populations that are genetically differentiated from offshore populations, indicative of local adaptation, the degree of which appears to be influenced by connectivity. Analyses of the genomic architecture reveal a significant overrepresentation of a large ~5 Mb chromosomal rearrangement in fjord cod, previously proposed to comprise genes critical for the survival at low salinities. This suggests that despite considerable connectivity with offshore populations, local adaptation to fjord environments may be enabled by suppression of recombination in the rearranged region. Our study provides new insights into the potential of local adaptation in high gene flow species within fine geographical scales and highlights the importance of genome architecture in analyses of ecological adaptation.

The role of the Strait of Gibraltar in shaping the genetic structure of the Mediterranean Grenadier, Coryphaenoides mediterraneus, between the Atlantic and Mediterranean Sea.

Population genetic studies of species inhabiting the deepest parts of the oceans are still scarce and only until recently we started to understand how oceanographic processes and topography affect dispersal and gene flow patterns. The aim of this study was to investigate the spatial population genetic structure of the bathyal bony fish Coryphaenoides mediterraneus, with a focus on the Atlantic-Mediterranean transition. We used nine nuclear microsatellites and the mitochondrial cytochrome c oxidase I gene from 6 different sampling areas. No population genetic structure was found within Mediterranean with both marker types (mean ΦST = 0.0960, FST = -0.0003, for both P > 0.05). However, within the Atlantic a contrasting pattern of genetic structure was found for the mtDNA and nuclear markers (mean ΦST = 0.2479, P < 0.001; FST = -0.0001, P > 0.05). When comparing samples from Atlantic and Mediterranean they exhibited high and significant levels of genetic divergence (mean ΦST = 0.7171, FST = 0.0245, for both P < 0.001) regardless the genetic marker used. Furthermore, no shared haplotypes were found between Atlantic and Mediterranean populations. These results suggest very limited genetic exchange between Atlantic and Mediterranean populations of C. mediterraneus, likely due to the shallow bathymetry of the Strait of Gibraltar acting as a barrier to gene flow. This physical barrier not only prevents the direct interactions between the deep-living adults, but also must prevent interchange of pelagic early life stages between the two basins. According to Bayesian simulations it is likely that Atlantic and Mediterranean populations of C. mediterraneus were separated during the late Pleistocene, which is congruent with results for other deep-sea fish from the same region.

Habitat Discontinuities Separate Genetically Divergent Populations of a Rocky Shore Marine Fish.

Habitat fragmentation has been suggested to be responsible for major genetic differentiations in a range of marine organisms. In this study, we combined genetic data and environmental information to unravel the relative role of geography and habitat heterogeneity on patterns of genetic population structure of corkwing wrasse (Symphodus melops), a rocky shore species at the northern limit of its distribution range in Scandinavia. Our results revealed a major genetic break separating populations inhabiting the western and southern coasts of Norway. This genetic break coincides with the longest stretch of sand in the whole study area, suggesting habitat fragmentation as a major driver of genetic differentiation of this obligate rocky shore benthic fish in Scandinavia. The complex fjords systems extending along the western coast of Norway appeared responsible for further regional genetic structuring. Our findings indicate that habitat discontinuities may lead to significant genetic fragmentation over short geographical distances, even for marine species with a pelagic larval phase, as for this rocky shore fish.

"Islands of Divergence" in the Atlantic Cod Genome Represent Polymorphic Chromosomal Rearrangements.

In several species genetic differentiation across environmental gradients or between geographically separate populations has been reported to center at "genomic islands of divergence," resulting in heterogeneous differentiation patterns across genomes. Here, genomic regions of elevated divergence were observed on three chromosomes of the highly mobile fish Atlantic cod (Gadus morhua) within geographically fine-scaled coastal areas. The "genomic islands" extended at least 5, 9.5, and 13 megabases on linkage groups 2, 7, and 12, respectively, and coincided with large blocks of linkage disequilibrium. For each of these three chromosomes, pairs of segregating, highly divergent alleles were identified, with little or no gene exchange between them. These patterns of recombination and divergence mirror genomic signatures previously described for large polymorphic inversions, which have been shown to repress recombination across extensive chromosomal segments. The lack of genetic exchange permits divergence between noninverted and inverted chromosomes in spite of gene flow. For the rearrangements on linkage groups 2 and 12, allelic frequency shifts between coastal and oceanic environments suggest a role in ecological adaptation, in agreement with recently reported associations between molecular variation within these genomic regions and temperature, oxygen, and salinity levels. Elevated genetic differentiation in these genomic regions has previously been described on both sides of the Atlantic Ocean, and we therefore suggest that these polymorphisms are involved in adaptive divergence across the species distributional range.

Population structure in Atlantic cod in the eastern North Sea-Skagerrak-Kattegat: early life stage dispersal and adult migration.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

The Pillars of Hercules as a bathymetric barrier to gene flow promoting isolation in a global deep-sea shark (Centroscymnus coelolepis).

Knowledge of the mechanisms limiting connectivity and gene flow in deep-sea ecosystems is scarce, especially for deep-sea sharks. The Portuguese dogfish (Centroscymnus coelolepis) is a globally distributed and near threatened deep-sea shark. C. coelolepis population structure was studied using 11 nuclear microsatellite markers and a 497-bp fragment from the mtDNA control region. High levels of genetic homogeneity across the Atlantic (Φ(ST) = -0.0091, F(ST) = 0.0024, P > 0.05) were found suggesting one large population unit at this basin. The low levels of genetic divergence between Atlantic and Australia (Φ(ST) = 0.0744, P < 0.01; F(ST) = 0.0015, P > 0.05) further suggested that this species may be able to maintain some degree of genetic connectivity even across ocean basins. In contrast, sharks from the Mediterranean Sea exhibited marked genetic differentiation from all other localities studied (Φ(ST) = 0.3808, F(ST) = 0.1149, P < 0.001). This finding suggests that the shallow depth of the Strait of Gibraltar acts as a barrier to dispersal and that isolation and genetic drift may have had an important role shaping the Mediterranean shark population over time. Analyses of life history traits allowed the direct comparison among regions providing a complete characterization of this shark's populations. Sharks from the Mediterranean had markedly smaller adult body size and size at maturity compared to Atlantic and Pacific individuals. Together, these results suggest the existence of an isolated and unique population of C. coelolepis inhabiting the Mediterranean that most likely became separated from the Atlantic in the late Pleistocene.

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.

Samples from subdivided populations yield biased estimates of effective size that overestimate the rate of loss of genetic variation.

Many empirical studies estimating effective population size apply the temporal method that provides an estimate of the variance effective size through the amount of temporal allele frequency change under the assumption that the study population is completely isolated. This assumption is frequently violated, and the magnitude of the resulting bias is generally unknown. We studied how gene flow affects estimates of effective size obtained by the temporal method when sampling from a population system and provide analytical expressions for the expected estimate under an island model of migration. We show that the temporal method tends to systematically underestimate both local and global effective size when populations are connected by gene flow, and the bias is sometimes dramatic. The problem is particularly likely to occur when sampling from a subdivided population where high levels of gene flow obscure identification of subpopulation boundaries. In such situations, sampling in a manner that prevents biased estimates can be difficult. This phenomenon might partially explain the frequently reported unexpectedly low effective population sizes of marine populations that have raised concern regarding the genetic vulnerability of even exceptionally large populations.