Adaptive divergence in a high gene flow environment: Hsc70 variation in the European flounder (Platichthys flesus L.)

Population Genetic Study on the European Flounder (Platichthys flesus) from the Southern Baltic Sea Using SNPs and Microsatellite Markers

The European flounder (Platichthys flesus), which is closely related to the recently discovered Baltic flounder (Platichthys solemdali), is currently the third most commercially fished species in the Baltic Sea. According to the available data from the Polish Fisheries Monitoring Center and fishermen's observations, the body condition indices of the species in the Baltic Sea have declined in recent years. The aim of the present study was to obtain information on the current patterns of genetic variability and the population structure of the European flounder and to verify whether the Baltic flounder is present in the southern Baltic Sea. Moreover, we aimed to verify whether the observed decline in the body condition indices of the species in the Baltic Sea might be associated with adaptive alterations in its gene pool due to increased fishing pressure. For this purpose, 190 fish were collected from four locations along the central coastline of Poland, i.e., Mechelinki, Władysławowo, the Vistula Lagoon in 2018, and the Słupsk Bank in 2020. The fish were morphologically analyzed and then genetically screened by the application of nineteen microsatellite DNA and two diagnostic SNP markers. The examined European flounder specimens displayed a high level of genetic diversity (PIC = 0.832-0.903, I = 2.579-2.768). A lack of significant genetic differentiation (Fst = 0.004, p > 0.05) was observed in all the examined fish, indicating that the European flounder in the sampled area constitutes a single genetic cluster. A significant deficiency in heterozygotes (Fis = 0.093, p < 0.05) and overall deviations from Hardy-Weinberg expectations (H-WE) were only detected in fish sampled from the Słupsk Bank. The estimated effective population size (Ne) among the sampled fish groups varied from 712 (Słupsk Bank) to 10,115 (Władysławowo and Mechelinki). However, the recorded values of the Garza-Williamson indicator (M = 0.574-0.600) and the lack of significant (p > 0.05) differences in Heq > He under the SMM model did not support the species' population size changes in the past. The applied SNP markers did not detect the presence of the Baltic flounder among the fish sampled from the studied area. The analysis of an association between biological traits and patterns of genetic diversity did not detect any signs of directional selection or density-dependent adaptive changes in the gene pool of the examined fish that might be caused by increased fishing pressure.

The role of genomic signatures of directional selection and demographic history in the population structure of a marine teleost with high gene flow

Recent studies have uncovered patterns of genomic divergence in marine teleosts where panmixia due to high gene flow has been the general paradigm. These signatures of divergent selection are often impacted by structural variants, acting as "supergenes" facilitating local adaptation. The highly dispersing European plaice (Pleuronectes platessa)-in which putative structural variants (i.e., inversions) have been identified-has successfully colonized the brackish water ecosystem of the Baltic Sea. Thus, the species represents an ideal opportunity to investigate how the interplay of gene flow, structural variants, natural selection, past demographic history, and gene flow impacts on population (sub)structuring in marine systems. Here, we report on the generation of an annotated draft plaice genome assembly in combination with population sequencing data-following the salinity gradient from the Baltic Sea into the North Sea together with samples from Icelandic waters-to illuminate genome-wide patterns of divergence. Neutral markers pointed at large-scale panmixia across the European continental shelf associated with high gene flow and a common postglacial colonization history of shelf populations. However, based on genome-wide outlier loci, we uncovered signatures of population substructuring among the European continental shelf populations, i.e., suggesting signs of ongoing selection. Genome-wide selection analyses (xp-EHH) and the identification of genes within genomic regions of recent selective sweeps-overlapping with the outlier loci-suggest that these represent the signs of divergent selection. Our findings provide support for genomic divergence driven by local adaptation in the face of high gene flow and elucidate the relative importance of demographic history versus adaptive divergence in shaping the contemporary population genetic structure of a marine teleost. The role of the putative inversion(s) in the substructuring-and potentially ongoing adaptation-was seemingly not substantial.

Genome-wide analysis reveals the genetic stock structure of hoki (Macruronus novaezelandiae)

The assessment of the genetic structuring of biodiversity is crucial for management and conservation. This is particularly critical for widely distributed and highly mobile deep-water teleosts, such as hoki (Macruronus novaezelandiae). This species is significant to Māori people and supports the largest commercial fishery in New Zealand, but uncertainty about its stock structure presents a challenge for management. Here, we apply a comprehensive genomic analysis to shed light on the demographic structure of this species by (1) assembling the genome, (2) generating a catalogue of genome-wide SNPs to infer the stock structure and (3) identifying regions of the genome under selection. The final genome assembly used short and long reads and is near complete, representing 93.8% of BUSCO genes, and consisting of 566 contigs totalling 501 Mb. Whole-genome re-sequencing of 510 hoki sampled from 14 locations around New Zealand and Australia, at a read depth greater than 10×, produced 227,490 filtered SNPs. Analyses of these SNPs were able to resolve the stock structure of hoki into two genetically and geographically distinct clusters, one including the Australian and the other one all New Zealand locations, indicating genetic exchange between these regions is limited. Location differences within New Zealand samples were much more subtle (global F ST  = 0.0006), and while small and significant differences could be detected, they did not conclusively identify additional substructures. Ten putative adaptive SNPs were detected within the New Zealand samples, but these also did not geographically partition the dataset further. Contemporary and historical N e estimation suggest the current New Zealand population of hoki is large yet declining. Overall, our study provides the first genomic resources for hoki and provides detailed insights into the fine-scale population genetic structure to inform the management of this species.

High gene flow in polar cod (Boreogadus saida) from West-Svalbard and the Eurasian Basin

The current and projected environmental change of the Arctic Ocean contrasts sharply with the limited knowledge of its genetic biodiversity. Polar cod Boreogadus saida (Lepechin, 1774) is an abundant circumpolar marine fish and ecological key species. The central role of polar cod in the Arctic marine food web warrants a better understanding of its population structure and connectivity. In this study, the genetic population structure of 171 juveniles, collected from several fjords off West-Svalbard (Billefjorden, Hornsund and Kongsfjorden), the northern Sophia Basin and the Eurasian Basin of the Arctic Ocean, was analysed using nine DNA microsatellite loci. Genetic analyses indicated moderate to high genetic diversity, but absence of spatial population structure and isolation-by-distance, suggesting ongoing gene flow between the studied sampling regions. High levels of connectivity may be key for polar cod to maintain populations across wide spatial scales. The adaptive capacity of the species will be increasingly important to face challenges such as habitat fragmentation, ocean warming and changes in prey composition. In view of a limited understanding of the population dynamics and evolution of polar cod, a valuable next step to predict future developments should be an integrated biological evaluation, including population genomics, a life-history approach, and habitat and biophysical dispersal modelling.

Evolution at two time frames: ancient structural variants involved in post-glacial divergence of the European plaice (Pleuronectes platessa)

Changing environmental conditions can lead to population diversification through differential selection on standing genetic variation. Structural variant (SV) polymorphisms provide examples of ancient alleles that in time become associated with novel environmental gradients. The European plaice (Pleuronectes platessa) is a marine flatfish showing large allele-frequency differences at two putative SVs associated with environmental variation. In this study, we explored the contribution of these SVs to population structure across the North East Atlantic. We compared genome-wide population structure using sets of RAD-sequencing SNPs with the spatial structure of the SVs. We found that in contrast to the rest of the genome, the SVs were only weakly associated with an isolation-by-distance pattern. Indeed, both SVs showed important variation in haplogroup frequencies, with the same haplogroup increasing both along the salinity gradient of the Baltic Sea, and found in high frequency in the northern-range margin of the Atlantic. Phylogenetic analyses suggested that the SV alleles are much older than the age of the Baltic Sea itself. These results suggest that the SVs are older than the age of the environmental gradients with which they currently co-vary. Altogether, our results suggest that the plaice SVs were shaped by evolutionary processes occurring at two time frames, firstly following their origin, ancient spread and maintenance in the ancestral populations, and secondly related to their current association with more recently formed environmental gradients such as those found in the North Sea-Baltic Sea transition zone.

The hemoglobin Gly16β1Asp polymorphism in turbot (Scophthalmus maximus) is differentially distributed across European populations

Turbot is an important flatfish widely distributed along the European coasts, whose fishery is centered in the North Sea. The commercial value of the species has boosted a successful aquaculture sector in Europe and China. Body growth is the main target of turbot breeding programs and is also a key trait related to local adaptation to temperature and salinity. Differences in growth rate and optimal growth temperature in turbot have been shown to be associated with a hemoglobin polymorphism reported more than 50 years ago. Here, we identified a Gly16Asp amino acid substitution in the β1 globin subunit by searching for genetic variation in the five functional globin genes within the whole annotated turbot genome. We predicted increased stability of the turbot hemoglobin by the replacement of the conserved Gly with the negative charged Asp residue that is consistent with the higher rate of αβ dimer assembly in the human J-Baltimore Gly16β->Asp mutant than in normal HbA. The turbot Hbβ1-Gly16 variant dominated in the northern populations examined, particularly in the Baltic Sea, while the Asp allele showed elevated frequencies in southern populations and was the prevalent variant in the Adriatic Sea. Body weight did not associate with the Hbβ1 genotypes at farming conditions (i.e., high oxygen levels, feeding ad libitum) after analyzing 90 fish with high growth dispersal from nine turbot families. Nevertheless, all data at hand suggest that the turbot hemoglobin polymorphism has an adaptive significance in the variable wild conditions regarding temperature and oxygen availability.

Spatial and environmental influences on selection in a clock gene coding trinucleotide repeat in Canada lynx (Lynx canadensis)

Clock genes exhibit substantial control over gene expression and ultimately life-histories using external cues such as photoperiod, and are thus likely to be critical for adaptation to shifting seasonal conditions and novel environments as species redistribute their ranges under climate change. Coding trinucleotide repeats (cTNRs) are found within several clock genes, and may be interesting targets of selection due to their containment within exonic regions and elevated mutation rates. Here, we conduct inter-specific characterization of the NR1D1 cTNR between Canada lynx and bobcat, and intra-specific spatial and environmental association analyses of neutral microsatellites and our functional cTNR marker, to investigate the role of selection on this locus in Canada lynx. We report signatures of divergent selection between lynx and bobcat, with the potential for hybrid-mediated gene flow in the area of range overlap. We also provide evidence that this locus is under selection across Canada lynx in eastern Canada, with both spatial and environmental variables significantly contributing to the explained variation, after controlling for neutral population structure. These results suggest that cTNRs may play an important role in the generation of functional diversity within some mammal species, and allow for contemporary rates of adaptation in wild populations in response to environmental change. We encourage continued investment into the study of cTNR markers to better understand their broader relevance to the evolution and adaptation of mammals.

Skim-Sequencing Based Genotyping Reveals Genetic Divergence of the Wild and Domesticated Population of Black Tiger Shrimp (Penaeus monodon) in the Indo-Pacific Region

The domestication of a wild-caught aquatic animal is an evolutionary process, which results in genetic discrimination at the genomic level in response to strong artificial selection. Although black tiger shrimp (Penaeus monodon) is one of the most commercially important aquaculture species, a systematic assessment of genetic divergence and structure of wild-caught and domesticated broodstock populations of the species is yet to be documented. Therefore, we used skim sequencing (SkimSeq) based genotyping approach to investigate the genetic structure of 50 broodstock individuals of P. monodon species, collected from five sampling sites (n = 10 in each site) across their distribution in Indo-Pacific regions. The wild-caught P. monodon broodstock population were collected from Malaysia (MS) and Japan (MJ), while domesticated broodstock populations were collected from Madagascar (MMD), Hawaii, HI, USA (MMO), and Thailand (MT). After various filtering process, a total of 194,259 single nucleotide polymorphism (SNP) loci were identified, in which 4983 SNP loci were identified as putatively adaptive by the pcadapt approach. In both datasets, pairwise F_ST estimates high genetic divergence between wild and domesticated broodstock populations. Consistently, different spatial clustering analyses in both datasets categorized divergent genetic structure into two clusters: (1) wild-caught populations (MS and MJ), and (2) domesticated populations (MMD, MMO and MT). Among 4983 putatively adaptive SNP loci, only 50 loci were observed to be in the coding region. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that non-synonymous mutated genes might be associated with the energy production, metabolic functions, respiration regulation and developmental rates, which likely act to promote adaptation to the strong artificial selection during the domestication process. This study has demonstrated the applicability of SkimSeq in a highly duplicated genome of P. monodon specifically, across a range of genetic backgrounds and geographical distributions, and would be useful for future genetic improvement program of this species in aquaculture.

Population Genomics of an Anadromous Hilsa Shad Tenualosa ilisha Species across Its Diverse Migratory Habitats: Discrimination by Fine-Scale Local Adaptation

The migration of anadromous fish in heterogenic environments unceasingly imposes a selective pressure that results in genetic variation for local adaptation. However, discrimination of anadromous fish populations by fine-scale local adaptation is challenging because of their high rate of gene flow, highly connected divergent population, and large population size. Recent advances in next-generation sequencing (NGS) have expanded the prospects of defining the weakly structured population of anadromous fish. Therefore, we used NGS-based restriction site-associated DNA (NextRAD) techniques on 300 individuals of an anadromous Hilsa shad (Tenualosa ilisha) species, collected from nine strategic habitats, across their diverse migratory habitats, which include sea, estuary, and different freshwater rivers. The NextRAD technique successfully identified 15,453 single nucleotide polymorphism (SNP) loci. Outlier tests using the FST OutFLANK and pcadapt approaches identified 74 and 449 SNPs (49 SNPs being common), respectively, as putative adaptive loci under a divergent selection process. Our results, based on the different cluster analyses of these putatively adaptive loci, suggested that local adaptation has divided the Hilsa shad population into two genetically structured clusters, in which marine and estuarine collection sites were dominated by individuals of one genetic cluster and different riverine collection sites were dominated by individuals of another genetic cluster. The phylogenetic analysis revealed that all the riverine populations of Hilsa shad were further subdivided into the north-western riverine (turbid freshwater) and the north-eastern riverine (clear freshwater) ecotypes. Among all of the putatively adaptive loci, only 36 loci were observed to be in the coding region, and the encoded genes might be associated with important biological functions related to the local adaptation of Hilsa shad. In summary, our study provides both neutral and adaptive contexts for the observed genetic divergence of Hilsa shad and, consequently, resolves the previous inconclusive findings on their population genetic structure across their diverse migratory habitats. Moreover, the study has clearly demonstrated that NextRAD sequencing is an innovative approach to explore how dispersal and local adaptation can shape genetic divergence of non-model anadromous fish that intersect diverse migratory habitats during their life-history stages.

Evidence for contemporary and historical gene flow between guppy populations in different watersheds, with a test for associations with adaptive traits

In dendritic river systems, gene flow is expected to occur primarily within watersheds. Yet, rare cross-watershed transfers can also occur, whether mediated by (often historical) geological events or (often contemporary) human activities. We explored these events and their potential evolutionary consequences by analyzing patterns of neutral genetic variation (microsatellites) and adaptive phenotypic variation (male color) in wild guppies (Poecilia reticulata) distributed across two watersheds in northern Trinidad. We found the expected signatures of within-watershed gene flow; yet we also inferred at least two instances of cross-watershed gene flow-one in the upstream reaches and one further downstream. The upstream cross-watershed event appears to be very recent (41 ± 13 years), suggesting dispersal via recent flooding or undocumented human-mediated transport. The downstream cross-watershed event appears to be considerably older (577 ± 265 years), suggesting a role for rare geological or climatological events. Alongside these strong signatures of both contemporary and historical gene flow, we found little evidence of impacts on presumably adaptive phenotypic differentiation, except perhaps in the one instance of very recent cross-watershed gene flow. Selection in this system seems to overpower gene flow-at least on the spatiotemporal scales investigated here.

Do genetic drift and accumulation of deleterious mutations preclude adaptation? Empirical investigation using RADseq in a northern lacustrine fish

Understanding genomic signatures of divergent selection underlying long-term adaptation in populations located in heterogeneous environments is a key goal in evolutionary biology. In this study, we investigated neutral, adaptive and deleterious genetic variation using 7,192 SNPs in 31 Lake Trout (Salvelinus namaycush) populations (n = 673) from Québec, Canada. Average genetic diversity was low, weakly shared among lakes, and positively correlated with lake size, indicating a major role for genetic drift subsequent to lake isolation. Putatively deleterious mutations were on average at lower frequencies than the other SNPs, and their abundance relative to the entire polymorphism in each population was positively correlated with inbreeding, suggesting that the effectiveness of purifying selection was negatively correlated with inbreeding, as predicted from theory. Despite evidence for pronounced genetic drift and inbreeding, several outlier loci were associated with temperature and found in or close to genes with biologically relevant functions notably related to heat stress and immune responses. Outcomes of gene-temperature associations were influenced by the inclusion of the most inbred populations, in which allele frequencies deviated the most from model predictions. This result illustrates challenge in identifying gene-environment associations in cases of high genetic drift and restricted gene flow and suggests limited adaptation in populations experiencing higher inbreeding. We discuss the relevance of these findings for the conservation and management, notably regarding stocking and genetic rescue, of Lake Trout populations and other species inhabiting highly fragmented habitats.

Immigrant reproductive dysfunction facilitates ecological speciation

The distributions of species are not only determined by where they can survive - they must also be able to reproduce. Although immigrant inviability is a well-established concept, the fact that immigrants also need to be able to effectively reproduce in foreign environments has not been fully appreciated in the study of adaptive divergence and speciation. Fertilization and reproduction are sensitive life-history stages that could be detrimentally affected for immigrants in non-native habitats. We propose that "immigrant reproductive dysfunction" is a hitherto overlooked aspect of reproductive isolation caused by natural selection on immigrants. This idea is supported by results from experiments on an externally fertilizing fish (sand goby, Pomatoschistus minutus). Growth and condition of adults were not affected by non-native salinity whereas males spawning as immigrants had lower sperm motility and hatching success than residents. We interpret these results as evidence for local adaptation or acclimation of sperm, and possibly also components of paternal care. The resulting loss in fitness, which we call "immigrant reproductive dysfunction," has the potential to reduce gene flow between populations with locally adapted reproduction, and it may play a role in species distributions and speciation.

Signatures of selection in mammalian clock genes with coding trinucleotide repeats: Implications for studying the genomics of high-pace adaptation

Climate change is predicted to affect the reproductive ecology of wildlife; however, we have yet to understand if and how species can adapt to the rapid pace of change. Clock genes are functional genes likely critical for adaptation to shifting seasonal conditions through shifts in timing cues. Many of these genes contain coding trinucleotide repeats, which offer the potential for higher rates of change than single nucleotide polymorphisms (SNPs) at coding sites, and, thus, may translate to faster rates of adaptation in changing environments. We characterized repeats in 22 clock genes across all annotated mammal species and evaluated the potential for selection on repeat motifs in three clock genes (NR1D1,CLOCK, and PER1) in three congeneric species pairs with different latitudinal range limits: Canada lynx and bobcat (Lynx canadensis and L. rufus), northern and southern flying squirrels (Glaucomys sabrinus and G. volans), and white‐footed and deer mouse (Peromyscus leucopus and P. maniculatus). Signatures of positive selection were found in both the interspecific comparison of Canada lynx and bobcat, and intraspecific analyses in Canada lynx. Northern and southern flying squirrels showed differing frequencies at common CLOCK alleles and a signature of balancing selection. Regional excess homozygosity was found in the deer mouse at PER1 suggesting disruptive selection, and further analyses suggested balancing selection in the white‐footed mouse. These preliminary signatures of selection and the presence of trinucleotide repeats within many clock genes warrant further consideration of the importance of candidate gene motifs for adaptation to climate change.

Proteomic responses of European flounder to temperature and hypoxia as interacting stressors: Differential sensitivities of populations

In the context of global change, ectotherms are increasingly impacted by abiotic perturbations. Along the distribution area of a species, the populations at low latitudes are particularly exposed to temperature increase and hypoxic events. In this study, we have compared the proteomic responses in the liver of European flounder populations, by using 2-D electrophoresis. One southern peripheral population from Portugal vs two northern core populations from France, were reared in a common garden experiment. Most of the proteomic differences were observed between the two experimental conditions, a cold vs a warm and hypoxic conditions. Consistent differentiations between populations were observed in accumulation of proteins involved in the bioenergetics- and methionine-metabolisms, fatty acids transport, and amino-acid catabolism. The specific regulation of crucial enzymes like ATP-synthase and G6PDH, in the liver of the southern population, could be related to a possible local adaptation. This southern peripheral population is spatially distant from northern core populations and has experienced dissimilar ecological conditions; thus it may contain genotypes that confer resilience to climate changes.

Adaptation and evolution of deep-sea scale worms (Annelida: Polynoidae): insights from transcriptome comparison with a shallow-water species

Polynoid scale worms (Polynoidae, Annelida) invaded deep-sea chemosynthesis-based ecosystems approximately 60 million years ago, but little is known about their genetic adaptation to the extreme deep-sea environment. In this study, we reported the first two transcriptomes of deep-sea polynoids (Branchipolynoe pettiboneae, Lepidonotopodium sp.) and compared them with the transcriptome of a shallow-water polynoid (Harmothoe imbricata). We determined codon and amino acid usage, positive selected genes, highly expressed genes and putative duplicated genes. Transcriptome assembly produced 98,806 to 225,709 contigs in the three species. There were more positively charged amino acids (i.e., histidine and arginine) and less negatively charged amino acids (i.e., aspartic acid and glutamic acid) in the deep-sea species. There were 120 genes showing clear evidence of positive selection. Among the 10% most highly expressed genes, there were more hemoglobin genes with high expression levels in both deep-sea species. The duplicated genes related to DNA recombination and metabolism, and gene expression were only enriched in deep-sea species. Deep-sea scale worms adopted two strategies of adaptation to hypoxia in the chemosynthesis-based habitats (i.e., rapid evolution of tetra-domain hemoglobin in Branchipolynoe or high expression of single-domain hemoglobin in Lepidonotopodium sp.).

Population genomics of an endemic Mediterranean fish: differentiation by fine scale dispersal and adaptation

The assessment of the genetic structuring of biodiversity is crucial for management and conservation. For species with large effective population sizes a low number of markers may fail to identify population structure. A solution of this shortcoming can be high-throughput sequencing that allows genotyping thousands of markers on a genome-wide approach while facilitating the detection of genetic structuring shaped by selection. We used Genotyping-by-Sequencing (GBS) on 176 individuals of the endemic East Atlantic peacock wrasse (Symphodus tinca), from 6 locations in the Adriatic and Ionian seas. We obtained a total of 4,155 polymorphic SNPs and we observed two strong barriers to gene flow. The first one differentiated Tremiti Islands, in the northwest, from all the other locations while the second one separated east and south-west localities. Outlier SNPs potentially under positive selection and neutral SNPs both showed similar patterns of structuring, although finer scale differentiation was unveiled with outlier loci. Our results reflect the complexity of population genetic structure and demonstrate that both habitat fragmentation and positive selection are on play. This complexity should be considered in biodiversity assessments of different taxa, including non-model yet ecologically relevant organisms.

Discrete phenotypes are not underpinned by genome-wide genetic differentiation in the squat lobster Munida gregaria (Crustacea: Decapoda: Munididae): a multi-marker study covering the Patagonian shelf

BACKGROUND

DNA barcoding has demonstrated that many discrete phenotypes are in fact genetically distinct (pseudo)cryptic species. Genetically identical, isogenic individuals, however, can also express similarly different phenotypes in response to a trigger condition, e.g. in the environment. This alternative explanation to cryptic speciation often remains untested because it requires considerable effort to reject the hypothesis that the observed underlying genetic homogeneity of the different phenotypes may be trivially caused by too slowly evolving molecular markers. The widespread squat lobster Munida gregaria comprises two discrete ecotypes, gregaria s. str. and subrugosa, which were long regarded as different species due to marked differences in morphological, ecological and behavioral traits. We studied the morphometry and genetics of M. gregaria s. l. and tested (1) whether the phenotypic differences remain stable after continental-scale sampling and inclusion of different life stages, (2) and whether each phenotype is underpinned by a specific genotype.

RESULTS

A total number of 219 gregaria s. str. and subrugosa individuals from 25 stations encompassing almost entire range in South America were included in morphological and genetic analyses using nine unlinked hypervariable microsatellites and new COI sequences. Results from the PCA and using discriminant functions demonstrated that the morphology of the two forms remains discrete. The mitochondrial data showed a shallow, star-like haplotype network and complete overlap of genetic distances within and among ecotypes. Coalescent-based species delimitation methods, PTP and GMYC, coherently suggested that haplotypes of both ecotypes forms a single species. Although all microsatellite markers possess sufficient genetic variation, AMOVA, PCoA and Bayesian clustering approaches revealed no genetic clusters corresponding to ecotypes or geographic units across the entire South-American distribution. No evidence of isolation-by-distance could be detected for this species in South America.

CONCLUSIONS

Despite their pronounced bimodal morphologies and different lifestyles, the gregaria s. str. and subrugosa ecotypes form a single, dimorphic species M. gregaria s. l.. Based on adequate geographic coverage and multiple independent polymorphic loci, there is no indication that each phenotype may have a unique genetic basis, leaving phenotypic plasticity or localized genomic islands of speciation as possible explanations.

Temporal stability of genetic variability and differentiation in the three-spined stickleback (Gasterosteus aculeatus)

Temporal variation in allele frequencies, whether caused by deterministic or stochastic forces, can inform us about interesting demographic and evolutionary phenomena occurring in wild populations. In spite of the continued surge of interest in the genetics of three-spined stickleback (Gasterosteus aculeatus) populations, little attention has been paid towards the temporal stability of allele frequency distributions, and whether there are consistent differences in effective size (Ne) of local populations. We investigated temporal stability of genetic variability and differentiation in 15 microsatellite loci within and among eight collection sites of varying habitat type, surveyed twice over a six-year time period. In addition, Nes were estimated with the expectation that they would be lowest in isolated ponds, intermediate in larger lakes and largest in open marine sites. In spite of the marked differences in genetic variability and differentiation among the study sites, the temporal differences in allele frequencies, as well as measures of genetic diversity and differentiation, were negligible. Accordingly, the Ne estimates were temporally stable, but tended to be lower in ponds than in lake or marine habitats. Hence, we conclude that allele frequencies in putatively neutral markers in three-spined sticklebacks seem to be temporally stable - at least over periods of few generations - across a wide range of habitat types differing markedly in levels of genetic variability, effective population size and gene flow.

Population genomic evidence for adaptive differentiation in Baltic Sea three-spined sticklebacks

Background

The degree of genetic differentiation among populations experiencing high levels of gene flow is expected to be low for neutral genomic sites, but substantial divergence can occur in sites subject to directional selection. Studies of highly mobile marine fish populations provide an opportunity to investigate this kind of heterogeneous genomic differentiation, but most studies to this effect have focused on a relatively low number of genetic markers and/or few populations. Hence, the patterns and extent of genomic divergence in high-gene-flow marine fish populations remain poorly understood.

Results

We here investigated genome-wide patterns of genetic variability and differentiation in ten marine populations of three-spined stickleback (Gasterosteus aculeatus) distributed across a steep salinity and temperature gradient in the Baltic Sea, by utilizing >30,000 single nucleotide polymorphisms obtained with a pooled RAD-seq approach. We found that genetic diversity and differentiation varied widely across the genome, and identified numerous fairly narrow genomic regions exhibiting signatures of both divergent and balancing selection. Evidence was uncovered for substantial genetic differentiation associated with both salinity and temperature gradients, and many candidate genes associated with local adaptation in the Baltic Sea were identified.

Conclusions

The patterns of genetic diversity and differentiation, as well as candidate genes associated with adaptation, in Baltic Sea sticklebacks were similar to those observed in earlier comparisons between marine and freshwater populations, suggesting that similar processes may be driving adaptation to brackish and freshwater environments. Taken together, our results provide strong evidence for heterogenic genomic divergence driven by local adaptation in the face of gene flow along an environmental gradient in the post-glacially formed Baltic Sea.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0130-8) contains supplementary material, which is available to authorized users.

A tale of two seas: contrasting patterns of population structure in the small-spotted catshark across Europe

Elasmobranchs represent important components of marine ecosystems, but they can be vulnerable to overexploitation. This has driven investigations into the population genetic structure of large-bodied pelagic sharks, but relatively little is known of population structure in smaller demersal taxa, which are perhaps more representative of the biodiversity of the group. This study explores spatial population genetic structure of the small-spotted catshark (Scyliorhinus canicula), across European seas. The results show significant genetic differences among most of the Mediterranean sample collections, but no significant structure among Atlantic shelf areas. The data suggest the Mediterranean populations are likely to have persisted in a stable and structured environment during Pleistocene sea-level changes. Conversely, the Northeast Atlantic populations would have experienced major changes in habitat availability during glacial cycles, driving patterns of population reduction and expansion. The data also provide evidence of male-biased dispersal and female philopatry over large spatial scales, implying complex sex-determined differences in the behaviour of elasmobranchs. On the basis of this evidence, we suggest that patterns of connectivity are determined by trends of past habitat stability that provides opportunity for local adaptation in species exhibiting philopatric behaviour, implying that resilience of populations to fisheries and other stressors may differ across the range of species.

Population genomics of marine fishes: next-generation prospects and challenges

Over the past few years, technological advances have facilitated giant leaps forward in our ability to generate genome-wide molecular data, offering exciting opportunities for gaining new insights into the ecology and evolution of species where genomic information is still limited. Marine fishes are valuable organisms for advancing our understanding of evolution on historical and contemporary time scales, and here we highlight areas in which research on these species is likely to be particularly important in the near future. These include possibilities for gaining insights into processes on ecological time scales, identifying genomic signatures associated with population divergence under gene flow, and determining the genetic basis of phenotypic traits. We also consider future challenges pertaining to the implementation of genome-wide coverage through next-generation sequencing and genotyping methods in marine fishes. Complications associated with fast decay of linkage disequilibrium, as expected for species with large effective population sizes, and the possibility that adaptation is associated with both soft selective sweeps and polygenic selection, leaving complex genomic signatures in natural populations, are likely to challenge future studies. However, the combination of high genome coverage and new statistical developments offers promising solutions. Thus, the next generation of studies is likely to truly facilitate the transition from population genetics to population genomics in marine fishes. This transition will advance our understanding of basic evolutionary processes and will offer new possibilities for conservation and management of valuable marine resources.

Importance of plasticity and local adaptation for coping with changing salinity in coastal areas: a test case with barnacles in the Baltic Sea

Background

Salinity plays an important role in shaping coastal marine communities. Near-future climate predictions indicate that salinity will decrease in many shallow coastal areas due to increased precipitation; however, few studies have addressed this issue. The ability of ecosystems to cope with future changes will depend on species’ capacities to acclimatise or adapt to new environmental conditions. Here, we investigated the effects of a strong salinity gradient (the Baltic Sea system – Baltic, Kattegat, Skagerrak) on plasticity and adaptations in the euryhaline barnacle Balanus improvisus. We used a common-garden approach, where multiple batches of newly settled barnacles from each of three different geographical areas along the Skagerrak-Baltic salinity gradient were exposed to corresponding native salinities (6, 15 and 30 PSU), and phenotypic traits including mortality, growth, shell strength, condition index and reproductive maturity were recorded.

Results

We found that B. improvisus was highly euryhaline, but had highest growth and reproductive maturity at intermediate salinities. We also found that low salinity had negative effects on other fitness-related traits including initial growth and shell strength, although mortality was also lowest in low salinity. Overall, differences between populations in most measured traits were weak, indicating little local adaptation to salinity. Nonetheless, we observed some population-specific responses – notably that populations from high salinity grew stronger shells in their native salinity compared to the other populations, possibly indicating adaptation to differences in local predation pressure.

Conclusions

Our study shows that B. improvisus is an example of a true brackish-water species, and that plastic responses are more likely than evolutionary tracking in coping with future changes in coastal salinity.

Fin whale MDH-1 and MPI allozyme variation is not reflected in the corresponding DNA sequences

The appeal of genetic inference methods to assess population genetic structure and guide management efforts is grounded in the correlation between the genetic similarity and gene flow among populations. Effects of such gene flow are typically genomewide; however, some loci may appear as outliers, displaying above or below average genetic divergence relative to the genomewide level. Above average population, genetic divergence may be due to divergent selection as a result of local adaptation. Consequently, substantial efforts have been directed toward such outlying loci in order to identify traits subject to local adaptation. Here, we report the results of an investigation into the molecular basis of the substantial degree of genetic divergence previously reported at allozyme loci among North Atlantic fin whale (Balaenoptera physalus) populations. We sequenced the exons encoding for the two most divergent allozyme loci (MDH-1 and MPI) and failed to detect any nonsynonymous substitutions. Following extensive error checking and analysis of additional bioinformatic and morphological data, we hypothesize that the observed allozyme polymorphisms may reflect phenotypic plasticity at the cellular level, perhaps as a response to nutritional stress. While such plasticity is intriguing in itself, and of fundamental evolutionary interest, our key finding is that the observed allozyme variation does not appear to be a result of genetic drift, migration, or selection on the MDH-1 and MPI exons themselves, stressing the importance of interpreting allozyme data with caution. As for North Atlantic fin whale population structure, our findings support the low levels of differentiation found in previous analyses of DNA nucleotide loci.

A review of molecular approaches for investigating patterns of coevolution in marine host-parasite relationships

Parasites and their relationships with hosts play a crucial role in the evolutionary pathways of every living organism. One method of investigating host-parasite systems is using a molecular approach. This is particularly important as analyses based solely on morphology or laboratory studies of parasites and their hosts do not take into account historical evolutionary interactions that can shape the distribution, abundance and population structure of parasites and their hosts. However, the predominant host-parasite coevolution literature has focused on terrestrial hosts and their parasites, and there still is a lack of studies in marine environments. Given that marine systems are generally more open than terrestrial ones, they provide fascinating opportunities for large-scale (as well as small-scale) geographic studies. Further, patterns and processes of genetic structuring and systematics are becoming more available across many different taxa (but especially fishes) in many marine systems, providing an excellent basis for examining whether parasites follow host population/species structure. In this chapter, we first highlight the factors and processes that challenge our ability to interpret evolutionary patterns of coevolution of hosts and their parasites in marine systems at different spatial, temporal and taxonomic scales. We then review the use of the most commonly utilized genetic markers in studying marine host-parasite systems. We give an overview and discuss which molecular methodologies resolve evolutionary relationships best and also discuss the applicability of new approaches, such as next-generation sequencing and studies utilizing functional markers to gain insights into more contemporary processes shaping host-parasite relationships.

Local adaptation to salinity in the three-spined stickleback?

Different lines of evidence suggest that the occurrence and extent of local adaptation in high gene flow marine environments - even in mobile and long-lived vertebrates with complex life cycles - may be more widespread than earlier thought. We conducted a common garden experiment to test for local adaptation to salinity in Baltic Sea sticklebacks (Gasterosteus aculeatus). Fish from three different native salinity regimes (high, mid and low) were subjected to three salinity treatments (high, mid and low) in a full-factorial experimental design. Irrespective of their origin, fish subjected to low (and mid) salinity treatments exhibited higher juvenile survival, grew to largest sizes and were in better condition than fish subjected to the high salinity treatment. However, a significant interaction between native and treatment salinities - resulting mainly from the poor performance of fish native to low salinity in the high salinity treatment - provided clear cut evidence for adaptation to local variation in salinity. Additional support for this inference was provided by the fact that the results concur with an earlier demonstration of significant differentiation in a number of genes with osmoregulatory functions across the same populations and that the population-specific responses to salinity treatments exceeded that to be expected by random genetic drift.

Population genetic studies revealed local adaptation in a high gene-flow marine fish, the small yellow croaker (Larimichthys polyactis)

The genetic differentiation of many marine fish species is low. Yet local adaptation may be common in marine fish species as the vast and changing marine environment provides more chances for natural selection. Here, we used anonymous as well as known protein gene linked microsatellites and mitochondrial DNA to detect the population structure of the small yellow croaker (Larimichthys polyactis) in the Northwest Pacific marginal seas. Among these loci, we detected at least two microsatellites, anonymous H16 and HSP27 to be clearly under diversifying selection in outlier tests. Sequence cloning and analysis revealed that H16 was located in the intron of BAHCC1 gene. Landscape genetic analysis showed that H16 mutations were significantly associated with temperature, which further supported the diversifying selection at this locus. These marker types presented different patterns of population structure: (i) mitochondrial DNA phylogeny showed no evidence of genetic divergence and demonstrated only one glacial linage; (ii) population differentiation using putatively neutral microsatellites presented a pattern of high gene flow in the L. polyactis. In addition, several genetic barriers were identified; (iii) the population differentiation pattern revealed by loci under diversifying selection was rather different from that revealed by putatively neutral loci. The results above suggest local adaptation in the small yellow croaker. In summary, population genetic studies based on different marker types disentangle the effects of demographic history, migration, genetic drift and local adaptation on population structure and also provide valuable new insights for the design of management strategies in L. polyactis.

High degree of genetic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene-based genome-scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection-as determined with several outlier detection methods-was low (F(ST) = 0.033 or 0.011, respectively), whereas average FST for directionally selected markers was significantly higher (F(ST) = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.

Population genetics provides new insights into biomarker prevalence in dab (Limanda limanda L.): a key marine biomonitoring species

Bioindicators are species for which some quantifiable aspect of its biology, a biomarker, is assumed to be sensitive to ecosystem health. However, there is frequently a lack of information on the underlying genetic and environmental drivers shaping the spatiotemporal variance in prevalence of the biomarkers employed. Here, we explore the relative role of potential variables influencing the spatiotemporal prevalence of biomarkers in dab, Limanda limanda, a species used as a bioindicator of marine contaminants. Firstly, the spatiotemporal genetic structure of dab around UK waters (39 samples across 15 sites for four years: 2005–2008) is evaluated with 16 microsatellites. Two temporally stable groups are identified corresponding to the North and Irish Seas (average between basin = 0.007; = 0.022). Secondly, we examine the association between biomarker prevalence and several variables, including genetic structuring, age and contaminant exposure. Genetic structure had significant interactive effects, together with age and some contaminants, in the prevalence of some of the biomarkers considered, namely hyperpigmentation and liver lesions. The integration of these data sets enhanced our understanding of the relationship between biomarker prevalence, exposure to contaminants and population-specific response, thereby yielding more informative predictive models of response and prospects for environmental remediation.

Detection of molecular signatures of selection at microsatellite loci in the South African abalone (Haliotis midae) using a population genomic approach

Identifying genomic regions that may be under selection is important for elucidating the genetic architecture of complex phenotypes underlying adaptation to heterogeneous environments. A population genomic approach, using a classical neutrality test and various Fst-outlier detection methods was employed to evaluate genome-wide polymorphism data in order to identify loci that may be candidates for selection amongst six populations (three cultured and three wild) of the South African abalone, Haliotis midae. Approximately 9% of the genome-wide microsatellite markers were putatively subject to directional selection, whilst 6-18% of the genome is thought to be influenced by balancing selection. Genetic diversity estimates for candidate loci under directional selection was significantly reduced in comparison to candidate neutral loci, whilst candidate balancing selection loci demonstrated significantly higher levels of genetic diversity (Kruskal-Wallis test, P<0.05). Pairwise Fst estimates based on candidate directional selection loci also demonstrated increased levels of differentiation between study populations. Various candidate loci under selection showed significant inter-chromosomal linkage disequilibrium, suggesting possible gene-networks underling adaptive phenotypes. Furthermore, several loci had significant hits to known genes when performing BLAST searches to NCBI's non-redundant databases, whilst others are known to be derived from expressed sequences even though homology to a known gene could not be established. A number of loci also demonstrated relatively high similarity to transposable elements. The association of these loci to functional and genomically active sequences could in part explain the observed signatures of selection.

Heterogeneous Genomic Differentiation in marine threespine sticklebacks: adaptation along an environmental gradient

Evolutionary divergence among populations occupying ecologically distinct environments can occur even in the face of on-going gene flow. However, the genetic underpinnings, as well as the scale and magnitude at which this differentiation occurs in marine habitats are not well understood. We investigated the patterns and degree of genomic heterogeneity in threespine sticklebacks (Gasterosteus aculeatus) by assessing genetic variability in 20 nongenic and 20 genic (associated with genes important for freshwater adaptation) microsatellite loci in samples collected from 38 locations spanning the entire Baltic Sea coast to the North Sea boundary. Population divergence (F(ST) ≈ 0.026) and structuring (five genetic clusters) was significantly more pronounced in the genic as compared to nongenic markers (F(ST) ≈ 0.008; no genetic clusters). Patterns of divergence in the genic markers--45% of which were identified as outliers--correlated with local differences in salinity. Yet, a strong positive correlation between divergence in genic and nongenic markers, and their association with environmental factors suggests that adaptive divergence is reducing gene flow across the genome. Apart from providing a clear demonstration of heterogeneous genomic patterns of differentiation in a marine species, the results are indicative of adaptive population structuring across the relatively young Baltic Sea in spite of ample opportunities for gene flow.

Molecular signatures of lineage-specific adaptive evolution in a unique sea basin: the example of an anadromous goby Leucopsarion petersii

Climate changes on various time scales often shape genetic novelty and adaptive variation in many biotas. We explored molecular signatures of directional selection in populations of the ice goby Leucopsarion petersii inhabiting a unique sea basin, the Sea of Japan, where a wide variety of environments existed in the Pleistocene in relation to shifts in sea level by repeated glaciations. This species consisted of two historically allopatric lineages, the Japan Sea (JS) and Pacific Ocean (PO) lineages, and these have lived under contrasting marine environments that are expected to have imposed different selection regimes caused by past climatic and current oceanographic factors. We applied a limited genome-scan approach using seven candidate genes for phenotypic differences between two lineages in combination with 100 anonymous microsatellite loci. Neuropeptide Y (NPY) gene, which is an important regulator of food intake and potent orexigenic agent, and three anonymous microsatellites were identified as robust outliers, that is, candidate loci potentially under directional selection, by multiple divergence- and diversity-based outlier tests in comparisons focused on multiple populations of the JS vs. PO lineages. For these outlier loci, populations of the JS lineage had putative signals of selective sweeps. Additionally, real-time quantitative PCR analysis using fish reared in a common environment showed a higher expression level for NPY gene in the JS lineage. Thus, this study succeeded in identifying candidate genomic regions under selection across populations of the JS lineage and provided evidence for lineage-specific adaptive evolution in this unique sea basin.

Microsatellite length variation in candidate genes correlates with habitat in the gilthead sea bream Sparus aurata

The genetic basis and evolutionary implications of local adaptation in high gene flow marine organisms are still poorly understood. In several Mediterranean fish species, alternative migration patterns exist between individuals entering coastal lagoons that offer favourable conditions for growth and those staying in the sea where environmental conditions are less subject to rapid and stressful change. Whether these coexisting strategies are phenotypically plastic or include a role for local adaptation through differential survival needs to be determined. Here, we explore the genetic basis of alternate habitat use in western Mediterranean populations of the gilthead sea bream (Sparus aurata). Samples from lagoonal and open-sea habitats were typed for three candidate gene microsatellite loci, seven anonymous microsatellites and 44 amplified fragment length polymorphism markers to test for genotype-environment associations. While anonymous markers globally indicated high levels of gene flow across geographic locations and habitats, non-neutral differentiation patterns correlated with habitat type were found at two candidate microsatellite loci located in the promoter region of the growth hormone and prolactin genes. Further analysis of these two genes revealed that a mechanism based on habitat choice alone could not explain the distribution of genotype frequencies at a regional scale, thus implying a role for differential survival between habitats. We also found an association between allele size and habitat type, which, in the light of previous studies, suggests that polymorphisms in the proximal promoter region could influence gene expression by modulating transcription factor binding, thus providing a potential explanatory link between genotype and growth phenotype in nature.

Exploring neutral and adaptive processes in expanding populations of gilthead sea bream, Sparus aurata L., in the North-East Atlantic

Recent studies in empirical population genetics have highlighted the importance of taking into account both neutral and adaptive genetic variation in characterizing microevolutionary dynamics. Here, we explore the genetic population structure and the footprints of selection in four populations of the warm-temperate coastal fish, the gilthead sea bream (Sparus aurata), whose recent northward expansion has been linked to climate change. Samples were collected at four Atlantic locations, including Spain, Portugal, France and the South of Ireland, and genetically assayed using a suite of species-specific markers, including 15 putatively neutral microsatellites and 23 expressed sequence tag-linked markers, as well as a portion of the mitochondrial DNA (mtDNA) control region. Two of the putatively neutral markers, Bld-10 and Ad-10, bore signatures of strong directional selection, particularly in the newly established Irish population, although the potential 'surfing effect' of rare alleles at the edge of the expansion front was also considered. Analyses after the removal of these loci suggest low but significant population structure likely affected by some degree of gene flow counteracting random genetic drift. No signal of historic divergence was detected at mtDNA. BLAST searches conducted with all 38 markers used failed to identify specific genomic regions associated to adaptive functions. However, the availability of genomic resources for this commercially valuable species is rapidly increasing, bringing us closer to the understanding of the interplay between selective and neutral evolutionary forces, shaping population divergence of an expanding species in a heterogeneous milieu.

Comparative population genetic analysis of bocaccio rockfish Sebastes paucispinis using anonymous and gene-associated simple sequence repeat loci

Comparative population genetic analyses of traditional and emergent molecular markers aid in determining appropriate use of new technologies. The bocaccio rockfish Sebastes paucispinis is a high gene-flow marine species off the west coast of North America that experienced strong population decline over the past 3 decades. We used 18 anonymous and 13 gene-associated simple sequence repeat (SSR) loci (expressed sequence tag [EST]-SSRs) to characterize range-wide population structure with temporal replicates. No F(ST)-outliers were detected using the LOSITAN program, suggesting that neither balancing nor divergent selection affected the loci surveyed. Consistent hierarchical structuring of populations by geography or year class was not detected regardless of marker class. The EST-SSRs were less variable than the anonymous SSRs, but no correlation between F(ST) and variation or marker class was observed. General linear model analysis showed that low EST-SSR variation was attributable to low mean repeat number. Comparative genomic analysis with Gasterosteus aculeatus, Takifugu rubripes, and Oryzias latipes showed consistently lower repeat number in EST-SSRs than SSR loci that were not in ESTs. Purifying selection likely imposed functional constraints on EST-SSRs resulting in low repeat numbers that affected diversity estimates but did not affect the observed pattern of population structure.

A review of the application of molecular genetics for fisheries management and conservation of sharks and rays

Since the first investigation 25 years ago, the application of genetic tools to address ecological and evolutionary questions in elasmobranch studies has greatly expanded. Major developments in genetic theory as well as in the availability, cost effectiveness and resolution of genetic markers were instrumental for particularly rapid progress over the last 10 years. Genetic studies of elasmobranchs are of direct importance and have application to fisheries management and conservation issues such as the definition of management units and identification of species from fins. In the future, increased application of the most recent and emerging technologies will enable accelerated genetic data production and the development of new markers at reduced costs, paving the way for a paradigm shift from gene to genome-scale research, and more focus on adaptive rather than just neutral variation. Current literature is reviewed in six fields of elasmobranch molecular genetics relevant to fisheries and conservation management (species identification, phylogeography, philopatry, genetic effective population size, molecular evolutionary rate and emerging methods). Where possible, examples from the Indo-Pacific region, which has been underrepresented in previous reviews, are emphasized within a global perspective.

Matching genetics with oceanography: directional gene flow in a Mediterranean fish species

Genetic connectivity and geographic fragmentation are two opposing mechanisms determining the population structure of species. While the first homogenizes the genetic background across populations the second one allows their differentiation. Therefore, knowledge of processes affecting dispersal of marine organisms is crucial to understand their genetic distribution patterns and for the effective management of their populations. In this study, we use genetic analyses of eleven microsatellites in combination with oceanographic satellite and dispersal simulation data to determine distribution patterns for Serranus cabrilla, a ubiquitous demersal broadcast spawner, in the Mediterranean Sea. Pairwise population F(ST) values ranged between -0.003 and 0.135. Two genetically distinct clusters were identified, with a clear division located between the oceanographic discontinuities at the Ibiza Channel (IC) and the Almeria-Oran Front (AOF), revealing an admixed population in between. The Balearic Front (BF) also appeared to dictate population structure. Directional gene flow on the Spanish coast was observed as S. cabrilla dispersed from west to east over the AOF, from north to south on the IC and from south of the IC towards the Balearic Islands. Correlations between genetic and oceanographic data were highly significant. Seasonal changes in current patterns and the relationship between ocean circulation patterns and spawning season may also play an important role in population structure around oceanographic fronts.