Relative role of life-history traits and historical factors in shaping genetic population structure of sardines (Sardina pilchardus)

Genomic analysis of NE Atlantic sardine (Sardina pilchardus) reveals reduced variation in a recently established North Sea population and directs reconsideration of management units

The European sardine (Sardina pilchardus) is under intense fishing pressure and exhibits distributional/abundance shifts linked to environmental change. The current understanding of population demographics needed for sustainable management is uncertain due to concerns that previous genetic studies lacked resolution and limited sampling of sardine north of the Bay of Biscay. To address these issues, we performed mtDNA sequencing and genome wide SNP analysis of samples collected across the Bay of Biscay, Celtic Sea, English Channel and North Sea. The complete SNP data reported a lack of structure throughout the sampled area compatible with high gene flow. A consensus suite of positive outlier SNPs was identified which reported a significant correlation with geographical distance with the largest differentiation between the southern Bay of Biscay and North Sea samples which also reported a significant mtDNA ΦST. While the roles of dispersal limitation and environmental heterogeneity underpinning this require further study, this adds to growing evidence that selection is influencing sardine population structure against a background of high gene flow. The results indicate that while there may be a level of demographic independence between North Sea and South Biscay sardine, the current delimitation of central (Biscay) and northern (Channel and Celtic Sea) operational stocks may misrepresent connectivity between the Biscay and Channel. The North Sea sample exhibited markedly lower mtDNA and nuclear variation than other samples. As sardine have only recently invaded the North Sea such reduced genetic variation is compatible with predictions for peripheral leading-edge populations but contrasts with patterns for other small pelagic species and emphasises the need to consider species-specific genetic structure in ecosystem-based management. Nascent management of the North Sea sardine fishery must ensure that current low levels of genetic diversity are not eroded further as this may decrease the species adaptive potential and inhibit its expansion.

Santamaría MT, Faria R, Machado AM, Fonseca MM, Froufe E, C. Castro LF. Population Genomics Reveals the Underlying Structure of the Small Pelagic European Sardine and Suggests Low Connectivity within Macaronesia

The European sardine (Sardina pilchardus, Walbaum 1792) is indisputably a commercially important species. Previous studies using uneven sampling or a limited number of makers have presented sometimes conflicting evidence of the genetic structure of S. pilchardus populations. Here, we show that whole genome data from 108 individuals from 16 sampling areas across 5000 km of the species' distribution range (from the Eastern Mediterranean to the archipelago of Azores) support at least three genetic clusters. One includes individuals from Azores and Madeira, with evidence of substructure separating these two archipelagos in the Atlantic. Another cluster broadly corresponds to the center of the distribution, including the sampling sites around Iberia, separated by the Almeria-Oran front from the third cluster that includes all of the Mediterranean samples, except those from the Alboran Sea. Individuals from the Canary Islands appear to belong to the Mediterranean cluster. This suggests at least two important geographical barriers to gene flow, even though these do not seem complete, with many individuals from around Iberia and the Mediterranean showing some patterns compatible with admixture with other genetic clusters. Genomic regions corresponding to the top outliers of genetic differentiation are located in areas of low recombination indicative that genetic architecture also has a role in shaping population structure. These regions include genes related to otolith formation, a calcium carbonate structure in the inner ear previously used to distinguish S. pilchardus populations. Our results provide a baseline for further characterization of physical and genetic barriers that divide European sardine populations, and information for transnational stock management of this highly exploited species towards sustainable fisheries.

Phylogeography of the veined squid, Loligo forbesii, in European waters

The veined squid, Loligo forbesii Steenstrup, 1856, occurs at the European Shelf areas including the Azores and represents a valuable resource for the European commercial fishery in the North East Atlantic. However, very little is known about its population structure and phylogeography. This lack of knowledge also impedes the development of sustainable fishery management for this species. The present study combined the use of two types of markers that retrieve patterns of gene flow in different time spans; the analysis of 16 nuclear microsatellites and sequencing of the mitochondrial cytochrome oxidase subunit I (COI). Whereas the high mutation rate of microsatellites allows the description of recent patterns of connectivity in species, the lower mutation rate of COI provides phylogeographic patterns on a longer timescale. A total of 347 individuals of L. forbesii were investigated from nearly the entire distribution range of the species, including the North East Atlantic Shelf, the Azores and the Mediterranean. Individuals from the Western and Eastern Mediterranean Sea have never been included in a genetic study before. We were able to analyse COI sequences from all 12 sampling areas and define three clades of L. forbesii. Due to our large sampling area, we are presenting 13 COI-haplotypes that were previously unknown. The microsatellite analysis does not include the Azores but three main clades could be identified at the remaining 11 sampling sites. Low F_ST values indicate gene flow over large geographical distances. However, the genetically significant differences and an additional slight grouping in the microsatellite structure reveal that geographical barriers seem to influence the population structure and reduce gene flow. Furthermore, both markers provide strong evidence that the observed phylogeographic pattern reflects the geographical history of the Azores and the Mediterranean Sea.

Population Genetic Diversity of Two Marine Gobies (Gobiiformes: Gobiidae) from the North-Eastern Atlantic and the Mediterranean Sea

Gobies (Gobiiformes: Gobiidae) are the most species-rich family of fishes in general, and the most abundant fish group in the European seas. Nonetheless, our knowledge on many aspects of their biology, including the population genetic diversity, is poor. Although barriers to gene flow are less apparent in the marine environment, the ocean is not a continuous habitat, as has been shown by studies on population genetics of various marine biota. For the first time, European marine goby species which cannot be collected by common fishery techniques were studied. The population genetic structure of two epibenthic species, Gobius geniporus and Gobius cruentatus, from seven localities across their distribution ranges was assessed, using one mitochondrial (cytochrome b) and one nuclear gene (first intron of ribosomal protein gene S7). Our results showed that there is a great diversity of haplotypes of mitochondrial gene cytochrome b in both species at all localities. Global fixation indices (FST) indicated a great differentiation of populations in both studied gobies. Our results did not show a geographic subdivision to individual populations. Instead, the data correspond with the model of migration which allow divergence and recurrent migration from the ancestral population. The estimated migration routes coincide with the main currents in the studied area. This matches well the biology of the studied species, with adults exhibiting only short-distance movements and planktonic larval stages.

Connectivity and population structure of albacore tuna across southeast Atlantic and southwest Indian Oceans inferred from multidisciplinary methodology

Albacore tuna (Thunnus alalunga) is an important target of tuna fisheries in the Atlantic and Indian Oceans. The commercial catch of albacore is the highest globally among all temperate tuna species, contributing around 6% in weight to global tuna catches over the last decade. The accurate assessment and management of this heavily exploited resource requires a robust understanding of the species’ biology and of the pattern of connectivity among oceanic regions, yet Indian Ocean albacore population dynamics remain poorly understood and its level of connectivity with the Atlantic Ocean population is uncertain. We analysed morphometrics and genetics of albacore (n = 1,874) in the southwest Indian (SWIO) and southeast Atlantic (SEAO) Oceans to investigate the connectivity and population structure. Furthermore, we examined the species’ dispersal potential by modelling particle drift through major oceanographic features. Males appear larger than females, except in South African waters, yet the length–weight relationship only showed significant male–female difference in one region (east of Madagascar and Reunion waters). The present study produced a genetic differentiation between the southeast Atlantic and southwest Indian Oceans, supporting their demographic independence. The particle drift models suggested dispersal potential of early life stages from SWIO to SEAO and adult or sub-adult migration from SEAO to SWIO.

Genetic analysis redraws the management boundaries for the European sprat

Sustainable fisheries management requires detailed knowledge of population genetic structure. The European sprat is an important commercial fish distributed from Morocco to the Arctic circle, Baltic, Mediterranean, and Black seas. Prior to 2018, annual catch advice on sprat from the International Council for the Exploration of the Sea (ICES) was based on five putative stocks: (a) North Sea, (b) Kattegat-Skagerrak and Norwegian fjords, (c) Baltic Sea, (d) West of Scotland-southern Celtic Seas, and (e) English Channel. However, there were concerns that the sprat advice on stock size estimates management plan inadequately reflected the underlying biological units. Here, we used ddRAD sequencing to develop 91 SNPs that were thereafter used to genotype approximately 2,500 fish from 40 locations. Three highly distinct and relatively homogenous genetic groups were identified: (a) Norwegian fjords; (b) Northeast Atlantic including the North Sea, Kattegat-Skagerrak, Celtic Sea, and Bay of Biscay; and (c) Baltic Sea. Evidence of genetic admixture and possibly physical mixing was detected in samples collected from the transition zone between the North and Baltic seas, but not between any of the other groups. These results have already been implemented by ICES with the decision to merge the North Sea and the Kattegat-Skagerrak sprat to be assessed as a single unit, thus demonstrating that genetic data can be rapidly absorbed to align harvest regimes and biological units.

Weak population structure and recent demographic expansion of the monogenean parasite Kapentagyrus spp. infecting clupeid fishes of Lake Tanganyika, East Africa

Lake Tanganyika, East Africa, is the oldest and deepest African Great Lake and harbours one of the most diverse fish assemblages on earth. Two clupeid fishes, Limnothrissa miodon and Stolothrissa tanganicae, constitute a major part of the total fish catch, making them indispensable for local food security. Parasites have been proposed as indicators of stock structure in highly mobile pelagic hosts. We examined the monogeneans Kapentagyrus limnotrissae and Kapentagyrus tanganicanus (Dactylogyridae) infecting these clupeids to explore the parasites' lake-wide population structure and patterns of demographic history. Samples were collected at seven sites distributed across three sub-basins of the lake. Intraspecific morphological variation of the monogeneans (n = 380) was analysed using morphometrics and geomorphometrics of sclerotised structures. Genetic population structure of both parasite species (n = 246) was assessed based on a 415 bp fragment of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Overall, we observed a lack of clear geographical morphological differentiation in both parasites along a north-south axis. This lack of geographical population structure was also reflected by a large proportion of shared haplotypes, and a pattern of seemingly unrestricted gene flow between populations. Significant morphological and genetic differentiation between some populations might reflect temporal differentiation rather than geographical isolation. Overall, the shallow population structure of both species of Kapentagyrus reflects the near-panmictic population structure of both host species as previously reported. Morphological differences related to host species identity of K. tanganicanus were consistent with incipient speciation at the genetic level. Both parasite species experienced a recent demographic expansion, which might be linked to paleohydrological events. Finally, interspecific hybridisation was found in Kapentagyrus, representing the first case in dactylogyrid monogeneans.

Genetic homogeneity of the Pacific thread herring (Opisthonema libertate) (Günther, 1867) in the Eastern Pacific, inferred from mtDNA sequences

In the present study, the population genetic structure of the Pacific thread herring (Opisthonema libertate) was analyzed through mitochondrial DNA (mtDNA) control region sequences. Organisms were collected from June 2015 to July 2015 from four commercial landing sites (Bahia Magdalena, Guaymas and Mazatlan, Mexico, and Puntarenas, Costa Rica) and one artisanal sampling site (Puerto La Libertad, El Salvador). A total of 125 sequences were analyzed. High levels of haplotype (h = 0.990) and nucleotide (π = 0.030) diversity were found. Pairwise Φst comparisons indicated differences attributed mainly to the organisms from El Salvador. However, Bayesian inferences did not support the existence of different populations. The haplotype distribution between locations did not show a clear phylogeographic pattern. Mismatched distribution showed a unimodal pattern for the five sampled areas, indicative of sudden demographic expansion. These results were supported by Bayesian skyline plot. Our results do not support the hypothesis that the Pacific thread herring presents a population genetic structure. Future genetic comparisons should include a larger number of samples as well as more polymorphic molecular markers to further support our results.

First genomic study on Lake Tanganyika sprat Stolothrissa tanganicae: a lack of population structure calls for integrated management of this important fisheries target species

BACKGROUND

Clupeid fisheries in Lake Tanganyika (East Africa) provide food for millions of people in one of the world's poorest regions. Due to climate change and overfishing, the clupeid stocks of Lake Tanganyika are declining. We investigate the population structure of the Lake Tanganyika sprat Stolothrissa tanganicae, using for the first time a genomic approach on this species. This is an important step towards knowing if the species should be managed separately or as a single stock. Population structure is important for fisheries management, yet understudied for many African freshwater species. We hypothesize that distinct stocks of S. tanganicae could be present due to the large size of the lake (isolation by distance), limnological variation (adaptive evolution), or past separation of the lake (historical subdivision). On the other hand, high mobility of the species and lack of obvious migration barriers might have resulted in a homogenous population.

RESULTS

We performed a population genetic study on wild-caught S. tanganicae through a combination of mitochondrial genotyping (96 individuals) and RAD sequencing (83 individuals). Samples were collected at five locations along a north-south axis of Lake Tanganyika. The mtDNA data had low global FST and, visualised in a haplotype network, did not show phylogeographic structure. RAD sequencing yielded a panel of 3504 SNPs, with low genetic differentiation (FST = 0.0054; 95% CI: 0.0046-0.0066). PCoA, fineRADstructure and global FST suggest a near-panmictic population. Two distinct groups are apparent in these analyses (FST = 0.1338 95% CI: 0.1239,0.1445), which do not correspond to sampling locations. Autocorrelation analysis showed a slight increase in genetic difference with increasing distance. No outlier loci were detected in the RADseq data.

CONCLUSION

Our results show at most very weak geographical structuring of the stock and do not provide evidence for genetic adaptation to historical or environmental differences over a north-south axis. Based on these results, we advise to manage the stock as one population, integrating one management strategy over the four riparian countries. These results are a first comprehensive study on the population structure of these important fisheries target species, and can guide fisheries management.

Restricted Gene Flow for Gadus macrocephalus from Yellow Sea Based on Microsatellite Markers: Geographic Block of Tsushima Current

The Pacific cod Gadus macrocephalus is a demersal, economically important fish in the family Gadidae. Population genetic differentiation of Pacific cod was examined across its northwestern Pacific range by screening variation of eight microsatellite loci in the present study. All four populations exhibited high genetic diversity. Pairwise fixation index (Fst) suggested a moderate to high level of genetic differentiation among populations. Population of the Yellow Sea (YS) showed higher genetic difference compared to the other three populations based on the results of pairwise Fst, three-dimensional factorial correspondence analysis (3D-FCA) and STRUCTURE, which implied restricted gene flow among them. Wilcoxon signed rank tests suggested no significant heterozygosity excess and no recent genetic bottleneck events were detected. Microsatellite DNA is an effective molecular marker for detecting the phylogeographic pattern of Pacific cod, and these Pacific cod populations should be three management units.

Scale morphometry allows discrimination of European sardine Sardina pilchardus and round sardinella Sardinella aurita and among their local populations

Landmark-based geometric morphometric analysis revealed differences in scale shape between European sardine Sardina pilchardus and round sardinella Sardinella aurita as well as among the local populations of each species. Fish scale measurements from four different areas in the central and eastern Mediterranean Sea showed that the mean scale shape of the two species using landmark data could be differentiated with high certainty. Populations of S. aurita from the central and eastern Mediterranean Sea could be separated reliably (P < 0·001) with an average discrimination rate of 91%, whereas the average discrimination of the S. pilchardus populations was lower (80%), albeit still high.

The Story of a Hitchhiker: Population Genetic Patterns in the Invasive Barnacle Balanus(Amphibalanus) improvisus Darwin 1854

Understanding the ecological and evolutionary forces that determine the genetic structure and spread of invasive species is a key component of invasion biology. The bay barnacle, Balanus improvisus (= Amphibalanus improvisus), is one of the most successful aquatic invaders worldwide, and is characterised by broad environmental tolerance. Although the species can spread through natural larval dispersal, human-mediated transport through (primarily) shipping has almost certainly contributed to the current global distribution of this species. Despite its worldwide distribution, little is known about the phylogeography of this species. Here, we characterize the population genetic structure and model dispersal dynamics of the barnacle B. improvisus, and describe how human-mediated spreading via shipping as well as natural larval dispersal may have contributed to observed genetic variation. We used both mitochondrial DNA (cytochrome c oxidase subunit I: COI) and nuclear microsatellites to characterize the genetic structure in 14 populations of B. improvisus on a global and regional scale (Baltic Sea). Genetic diversity was high in most populations, and many haplotypes were shared among populations on a global scale, indicating that long-distance dispersal (presumably through shipping and other anthropogenic activities) has played an important role in shaping the population genetic structure of this cosmopolitan species. We could not clearly confirm prior claims that B. improvisus originates from the western margins of the Atlantic coasts; although there were indications that Argentina could be part of a native region. In addition to dispersal via shipping, we show that natural larval dispersal may play an important role for further colonisation following initial introduction.

Discovery of Genome-Wide Microsatellite Markers in Scombridae: A Pilot Study on Albacore Tuna

Recent developments in sequencing technologies and bioinformatics analysis provide a greater amount of DNA sequencing reads at a low cost. Microsatellites are the markers of choice for a variety of population genetic studies, and high quality markers can be discovered in non-model organisms, such as tuna, with these recent developments. Here, we use a high-throughput method to isolate microsatellite markers in albacore tuna, Thunnus alalunga, based on coupling multiplex enrichment and next-generation sequencing on 454 GS-FLX Titanium pyrosequencing. The crucial minimum number of polymorphic markers to infer evolutionary and ecological processes for this species has been described for the first time. We provide 1670 microsatellite design primer pairs, and technical and molecular genetics selection resulting in 43 polymorphic microsatellite markers. On this panel, we characterized 34 random and selectively neutral markers («neutral») and 9 «non-neutral» markers. The variability of «neutral» markers was screened with 136 individuals of albacore tuna from southwest Indian Ocean (42), northwest Indian Ocean (31), South Africa (31), and southeast Atlantic Ocean (32). Power analysis demonstrated that the panel of genetic markers can be applied in diversity and population genetics studies. Global genetic diversity for albacore was high with a mean number of alleles at 16.94; observed heterozygosity 66% and expected heterozygosity 77%. The number of individuals was insufficient to provide accurate results on differentiation. Of the 9 «non-neutral» markers, 3 were linked to a sequence of known function. The one is located to a sequence having an immunity function (ThuAla-Tcell-01) and the other to a sequence having energy allocation function (ThuAla-Hki-01). These two markers were genotyped on the 136 individuals and presented different diversity levels. ThuAla-Tcell-01 has a high number of alleles (20), heterozygosity (87–90%), and assignment index. ThuAla-Hki-01 has a lower number of alleles (9), low heterozygosity (24–27%), low assignment index and significant inbreeding. Finally, the 34 «neutral» and 3 «non-neutral» microsatellites markers were tested on four economically important Scombridae species—Thunnus albacares, Thunnus thynnus, Thunnus obesus, and Acanthocybium solandri.

Comparative pattern of genetic structure in two Mediterranean killifishes Aphanius fasciatus and Aphanius iberus inferred from both mitochondrial and nuclear data

In this study, genetic variation was assessed in Aphanius fasciatus and Aphanius iberus characterized by similar ecological traits but with very different distribution ranges in the Mediterranean area. Five populations of A. iberus and five of A. fasciatus were analysed using five polymorphic microsatellite loci and partial mitochondrial control region (D-loop) sequences. Congruent results were found with both nuclear and mitochondrial molecular markers. The results showed that similar levels of genetic divergence, based on mitochondrial control region sequences, are present among populations of A. iberus and among populations of A. fasciatus despite the very different geographic distance existing among the examined populations of the two species (low geographic distance in A. iberus and high in A. fasciatus). A possible explanation could be that the populations of A. iberus were isolated for a longer time than the populations of A. fasciatus supporting the hypothesis that the split in the lineage leading to A. iberus is older than the split in the lineage leading to A. fasciatus. The possibility that the wide circum-Mediterranean distribution of A. fasciatus ensures the high connectivity of its populations, preventing, in some cases, local differentiation, however, cannot be ruled out.

Testing the consistency of connectivity patterns for a widely dispersing marine species

Connectivity is widely recognized as an important component in developing effective management and conservation strategies. Although managers are generally most interested in demographic, rather than genetic connectivity, new analytic approaches are able to provide estimates of both demographic and genetic connectivity measures from genetic data. Combining such genetic data with mathematical models represents a powerful approach for accurately determining patterns of population connectivity. Here, we use microsatellite markers to investigate the genetic population structure of the New Zealand Rock Lobster, Jasus edwardsii, which has one of the longest known larval durations of all marine species (>2 years), a very large geographic range (>5500 km), and has been the subject of extensive dispersal modeling. Despite earlier mitochondrial DNA studies finding homogeneous genetic structure, the mathematical model suggests that there are source-sink dynamics for this species. We found evidence of genetic structure in J. edwardsii populations with three distinct genetic groups across New Zealand and a further Australian group; these groups and patterns of gene flow were generally congruent with the earlier mathematical model. Of particular interest was the consistent identification of a self-recruiting population/region from both modeling and genetic approaches. Although there is the potential for selection and harvesting to influence the patterns we observed, we believe oceanographic processes are most likely responsible for the genetic structure observed in J. edwardsii. Our results, using a species at the extreme end of the dispersal spectrum, demonstrate that source-sink population dynamics may still exist for such species.

Imprints from genetic drift and mutation imply relative divergence times across marine transition zones in a pan-European small pelagic fish (Sprattus sprattus)

Geographic distributions of most temperate marine fishes are affected by postglacial recolonisation events, which have left complex genetic imprints on populations of marine species. This study investigated population structure and demographic history of European sprat (Sprattus sprattus L.) by combining inference from both mtDNA and microsatellite genetic markers throughout the species' distribution. We compared effects from genetic drift and mutation for both genetic markers in shaping genetic differentiation across four transition zones. Microsatellite markers revealed significant isolation by distance and a complex population structure across the species' distribution (overall θ(ST)=0.038, P<0.01). Across transition zones markers indicated larger effects of genetic drift over mutations in the northern distribution of sprat contrasting a stronger relative impact of mutation in the species' southern distribution in the Mediterranean region. These results were interpreted to reflect more recent divergence times between northern populations in accordance with previous findings. This study demonstrates the usefulness of comparing inference from different markers and estimators of divergence for phylogeographic and population genetic studies in species with weak genetic structure, as is the case in many marine species.

Genetic assessment of population structure and connectivity in the threatened Mediterranean coral Astroides calycularis (Scleractinia, Dendrophylliidae) at different spatial scales

Understanding dispersal patterns, population structure and connectivity among populations is helpful in the management and conservation of threatened species. Molecular markers are useful tools as indirect estimators of these characteristics. In this study, we assess the population genetic structure of the orange coral Astroides calycularis in the Alboran Sea at local and regional scale, and at three localities outside of this basin. Bayesian clustering methods, traditional F-statistics and D(est) statistics were used to determine the patterns of genetic structure. Likelihood and coalescence approaches were used to infer migration patterns and effective population sizes. The results obtained reveal a high level of connectivity among localities separated by as much as 1 km and moderate levels of genetic differentiation among more distant localities, somewhat corresponding with a stepping-stone model of gene flow and connectivity. These data suggest that connectivity among populations of this coral is mainly driven by the biology of the species, with low dispersal abilities; in addition, hydrodynamic processes, oceanographic fronts and the distribution of rocky substrate along the coastline may influence larval dispersal.

Characterization of a Bayesian genetic clustering algorithm based on a Dirichlet process prior and comparison among Bayesian clustering methods

Background

A Bayesian approach based on a Dirichlet process (DP) prior is useful for inferring genetic population structures because it can infer the number of populations and the assignment of individuals simultaneously. However, the properties of the DP prior method are not well understood, and therefore, the use of this method is relatively uncommon. We characterized the DP prior method to increase its practical use.

Results

First, we evaluated the usefulness of the sequentially-allocated merge-split (SAMS) sampler, which is a technique for improving the mixing of Markov chain Monte Carlo algorithms. Although this sampler has been implemented in a preceding program, HWLER, its effectiveness has not been investigated. We showed that this sampler was effective for population structure analysis. Implementation of this sampler was useful with regard to the accuracy of inference and computational time. Second, we examined the effect of a hyperparameter for the prior distribution of allele frequencies and showed that the specification of this parameter was important and could be resolved by considering the parameter as a variable. Third, we compared the DP prior method with other Bayesian clustering methods and showed that the DP prior method was suitable for data sets with unbalanced sample sizes among populations. In contrast, although current popular algorithms for population structure analysis, such as those implemented in STRUCTURE, were suitable for data sets with uniform sample sizes, inferences with these algorithms for unbalanced sample sizes tended to be less accurate than those with the DP prior method.

Conclusions

The clustering method based on the DP prior was found to be useful because it can infer the number of populations and simultaneously assign individuals into populations, and it is suitable for data sets with unbalanced sample sizes among populations. Here we presented a novel program, DPART, that implements the SAMS sampler and can consider the hyperparameter for the prior distribution of allele frequencies to be a variable.

Mitochondrial phylogeography of the European sprat (Sprattus sprattus L., Clupeidae) reveals isolated climatically vulnerable populations in the Mediterranean Sea and range expansion in the northeast Atlantic

We examined the genetic structure of the European sprat (Sprattus sprattus) by means of a 530-bp sequence of the mitochondrial control region from 210 fish originating from seven sampling localities of its distributional range. Phylogeographical analysis of 128 haplotypes showed a phylogenetic separation into two major clades with the Strait of Sicily acting as a barrier to gene flow between them. While no population differentiation was observed based on analysis of molecular variance and net nucleotide differences between samples of the Baltic Sea, the North Sea and the Bay of Biscay nor between the Black Sea and the Bosporus, a strong population differentiation between these samples and two samples from the Mediterranean Sea was found. Further, the biggest genetic distance was observed within the Mediterranean Sea between the populations of the Gulf of Lyon and the Adriatic Sea, indicating genetic isolation of these regions. Low genetic diversities and star-like haplotype networks of both Mediterranean Sea populations point towards recent demographic expansion scenarios after low population size, which is further supported by negative F(S) values and unimodal mismatch distributions with a low mean. Along the northeast Atlantic coast, a northwards range expansion of a large and stable population can be assumed. The history of a diverse but differentiated Black Sea population remains unknown due to uncertainties in the palaeo-oceanography of this sea. Our genetic data did not confirm the presently used classification into subspecies but are only preliminary in the absence of nuclear genetic analyses.