Colour polymorphism and genetic variation in Idotea baltica populations from the Adriatic Sea and Baltic Sea

Phenotypic and genetic variation was studied in two of the four European subspecies of the marine isopod Idotea baltica; the Mediterranean I. b. basteri and the Baltic I. b. baltica. Spatial and temporal patterns of colour polymorphism were analysed in northern Adriatic and western Baltic Sea populations. Pronounced differences in phenotype composition were observed between populations of both subspecies as seen in the distribution of various colour variants bilineata, lineata, flavafusca and several combined forms). Compared with Adriatic samples, western Baltic Sea populations show higher phenotypic diversity. To obtain an estimate of the degree of genetic divergence between the subspecies, 12 gene-enzyme systems were investigated electrophoretically. The results obtained indicate a relatively high level of genetic variation; I. b. basteri from the nothern Adriatic tends to be more polymorphic and more heterozygous than I. b. baltica from the western Baltic. Both subspecies share identical electrophoretic mobilities of the homologous enzyme proteins examined; however, in allelic composition they exhibit significant differences at approximately half the number of loci scored. The genetic distance (Nei's D) measured at the subspecific level was 0.04. Amounts and geographical patterns of variation, observed both in colour phenotype and electrophoretic variation, are considered.     

Biocomplexity in Populations of European Anchovy in the Adriatic Sea

The sustained exploitation of marine populations requires an understanding of a species'' adaptive seascape so that populations can track environmental changes from short- and long-term climate cycles and from human development. The analysis of the distributions of genetic markers among populations, together with correlates of life-history and environmental variability, can provide insights into the extent of adaptive variation. Here, we examined genetic variability among populations of mature European anchovies (n = 531) in the Adriatic (13 samples) and Tyrrhenian seas (2 samples) with neutral and putative non-neutral microsatellite loci. These genetic markers failed to confirm the occurrence of two anchovy species in the Adriatic Sea, as previously postulated. However, we found fine-scale population structure in the Adriatic, especially in northern areas, that was associated with four of the 13 environmental variables tested. Geographic gradients in sea temperature, salinity and dissolved oxygen appear to drive adaptive differences in spawning time and early larval development among populations. Resolving adaptive seascapes in Adriatic anchovies provides a means to understand mechanisms underpinning local adaptation and a basis for optimizing exploitation strategies for sustainable harvests.     

Phylogeography of the red coral (<Emphasis Type="Italic">Corallium rubrum</Emphasis>): inferences on the evolutionary history of a temperate gorgonian

The red coral Corallium rubrum (Cnidaria, Octocorallia) is an exploited, long-lived sessile species from the Mediterranean Sea and the adjacent coastline in the Atlantic Ocean. Surveys of genetic variation using microsatellites have shown that populations of C. rubrum are characterized by strong differentiation at the local scale but a study of the phylogeography of this species was still lacking. Here, we used seven polymorphic microsatellite loci, together with sequence data from an intron of the elongation factor 1 (EF1) gene, to investigate the genetic structure of C. rubrum across its geographical range in the western Mediterranean Sea and in the Adriatic Sea. The EF1 sequences were also used to analyse the consequences of demographic fluctuations linked with past environmental change. Clustering analysis with microsatellite loci highlighted three to seven genetic groups with the distinction of North African and Adriatic populations; this distinction appeared significant with AMOVA and differentiation tests. Microsatellite and EF1 data extended the isolation by distance pattern previously observed for this species at the western Mediterranean scale. EF1 sequences confirmed the genetic differentiation observed between most samples with microsatellites. A statistical parsimony network of EF1 haplotypes provided no evidence of high sequence divergence among regions, suggesting no long-term isolation. Selective neutrality tests on microsatellites and EF1 were not significant but should be interpreted with caution in the case of EF1 because of the low sample sizes for this locus. Our results suggest that recent Quaternary environmental fluctuations had a limited impact on the genetic structure of C. rubrum.     

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

Little is known about local adaptations in marine fishes since population genetic surveys in these species have typically not applied genetic markers subject to selection. In this study, we used a candidate gene approach to investigate adaptive population divergence in the European flounder (Platichthys flesus L.) throughout the northeastern Atlantic. We contrasted patterns of genetic variation in a presumably neutral microsatellite baseline to patterns from a heat-shock cognate protein gene, Hsc70. Using two different neutrality tests we found that the microsatellite data set most likely represented a neutral baseline. In contrast, Hsc70 strongly deviated from neutral expectations. Importantly, when estimating standardized levels of population divergence (F(ST)'), we also found a large discrepancy in the patterns of structuring in the two data sets. Thus, samples grouped according to geographical or historical proximity with regards to microsatellites, but according to environmental similarities with regards to Hsc70. The differences between the data sets were particularly pronounced in pairwise comparisons involving populations in the western and central Baltic Sea. For instance, the genetic differentiation between geographically close Baltic Sea and North Sea populations was found to be 0.02 and 0.45 for microsatellites and Hsc70 respectively. Our results strongly suggest adaptive population divergence and indicate local adaptations at the DNA level in a background of high levels of gene flow, typically found in many marine fish species. Furthermore, this study highlights the usefulness of the candidate gene approach for demonstrating local selection in non-model organisms such as most marine fishes.     

Marine landscapes and population genetic structure of herring (Clupea harengus L.) in the Baltic Sea

Numerically small but statistically significant genetic differentiation has been found in many marine fish species despite very large census population sizes and absence of obvious barriers to migrating individuals. Analyses of morphological traits have previously identified local spawning groups of herring (Clupea harengus L.) in the environmentally heterogeneous Baltic Sea, whereas allozyme markers have not revealed differentiation. We analysed variation at nine microsatellite loci in 24 samples of spring-spawning herring collected at 11 spawning locations throughout the Baltic Sea. Significant temporal differentiation was observed at two locations, which we ascribe to sympatrically spawning but genetically divergent 'spawning waves'. Significant differentiation was also present on a geographical scale, though pairwise F(ST) values were generally low, not exceeding 0.027. Partial Mantel tests showed no isolation by geographical distance, but significant associations were observed between genetic differentiation and environmental parameters (salinity and surface temperature) (0.001 < P < or = 0.099), though these outcomes were driven mainly by populations in the southwestern Baltic Sea, which also exhibits the steepest environmental gradients. Application of a novel method for detecting barriers to gene flow by combining geographical coordinates and genetic differentiation allowed us to identify two zones of lowered gene flow. These zones were concordant with the separation of the Baltic Sea into major basins, with environmental gradients and with differences in migration behaviour. We suggest that similar use of landscape genetics approaches may increase the understanding of the biological significance of genetic differentiation in other marine fishes.     

Managing polyploidy in ex situ conservation genetics: the case of the critically endangered Adriatic sturgeon (Acipenser naccarii)

While the current expansion of conservation genetics enables to address more efficiently the management of threatened species, alternative methods for genetic relatedness data analysis in polyploid species are necessary. Within this framework, we present a standardized and simple protocol specifically designed for polyploid species that can facilitate management of genetic diversity, as exemplified by the ex situ conservation program for the tetraploid Adriatic sturgeon Acipenser naccarii. A critically endangered endemic species of the Adriatic Sea tributaries, its persistence is strictly linked to the ex situ conservation of a single captive broodstock currently decimated to about 25 individuals, which represents the last remaining population of Adriatic sturgeon of certain wild origin. The genetic variability of three F1 broodstocks available as future breeders was estimated based on mitochondrial and microsatellite information and compared with the variability of the parental generation. Genetic data showed that the F1 stocks have only retained part of the genetic variation present in the original stock due to the few parent pairs used as founders. This prompts for the urgent improvement of the current F1 stocks by incorporating new founders that better represent the genetic diversity available. Following parental allocation based on band sharing values, we set up a user-friendly tool for selection of candidate breeders according to relatedness between all possible parent-pairs that secures the use of non-related individuals. The approach developed here could also be applied to other endangered tetraploid sturgeon species overexploited for caviar production, particularly in regions lacking proper infrastructure and/or expertise.     

Reduced genetic variation and strong genetic population structure in the freshwater killifish Valencia letourneuxi (Valenciidae) based on nuclear and mitochondrial markers

The genetic variation of the critically endangered Corfu killifish (Valencia letourneuxi), an endemic freshwater fish species of the western Balkans, was assessed for nine populations sampled in eight water systems in western continental Greece, the Peloponnese and the Ionian Island of Corfu, using mitochondrial and microsatellite markers. The analyses were based on data from three mtDNA regions (D‐loop, COI and 16S rRNA sequences) and 14 microsatellite loci. Samples from the congeneric species Valencia hispanica and the phylogenetically closely related species Aphanius fasciatus were also used in the study as outgroups. Both the mitochondrial and the microsatellite analyses revealed three distinct population groupings associated with the geographical distribution of the populations: one southern group occupying rivers draining to the Patraikos Gulf, the second one including the populations flowing into the Amvrakikos Gulf and the third, more northern group, including the other populations from rivers in Corfu Island and Epirus flowing into the Ionian Sea. Within these groupings there is limited genetic differentiation between populations; in addition, there is reduced intrapopulation genetic variation, evidenced by low heterozygosity values, number of alleles and haplotype diversity. In terms of taxonomic implications and appropriate management actions for conservation, our data suggest that the major population groups should be regarded at least as three distinct conservation units (CUs), with translocation and restocking actions to take place only within the geographical range of the CU concerned. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 334–349.     

Population genetic structure of wild and hatchery black rockfish Sebastes inermis in Korea, assessed using cross-species microsatellite markers

The population structure of the black rockfish, Sebastes inermis (Sebastidae), was estimated using 10 microsatellite loci developed for S. schlegeli on samples of 174 individuals collected from three wild and three hatchery populations in Korea. Reduced genetic variation was detected in hatchery strains [overall number of alleles (N(A)) = 8.07; allelic richness (A(R)) = 7.37; observed heterozygosity (H(O)) = 0.641] compared with the wild samples (overall N(A) = 8.43; A(R) = 7.83; H(O) = 0.670), but the difference was not significant. Genetic differentiation among the populations was significant (overall F(ST) = 0.0237, P < 0.05). Pairwise F(ST) tests, neighbor-joining tree, and principal component analyses showed significant genetic heterogeneity among the hatchery strains and between wild and hatchery strains, but not among the wild populations, indicating high levels of gene flow along the southern coast of Korea, even though the black rockfish is a benthic, non-migratory marine species. Genetic differentiation among the hatchery strains could reflect genetic drift due to intensive breeding practices. Thus, in the interests of optimal resource management, genetic variation should be monitored and inbreeding controlled within stocks in commercial breeding programs. Information on genetic population structure based on cross-species microsatellite markers can aid in the proper management of S. inermis populations.     

Patterns of genetic diversity and biogeographical history of the tropical wetland tree, Pterocarpus officinalis (Jacq.), in the Caribbean basin

Studies examining intraspecific variation in plant species with widespread distributions and disjunct populations have mainly concentrated on temperate species. Here, we determined the genetic structure of a broadly distributed wetland tropical tree, Pterocarpus officinalis (Jacq.), from eight Neotropical populations using amplified length fragment polymorphisms (AFLP). AFLPs proved highly variable with almost half (48%) of the genetic variation at these loci occurring among individuals within populations. Nonetheless, there was a strong geographical pattern in the distribution of AFLP variation within P. officinalis. Caribbean and continental populations fell into two well-defined genetic clusters supported by the presence of a number of unique AFLP bands. Within these two regions, there were also strong genetic differences among populations, caused mainly by frequency differences in AFLP bands, making it difficult to determine the evolutionary relationships among populations. In addition, our analysis of P. officinalis revealed striking differences in the levels of AFLP variation among the eight populations sampled. In general, Caribbean populations had lower genetic diversity than continental populations. Moreover, there was a clear loss in AFLP diversity with distance from the continent among Caribbean populations. The overall genetic pattern within P. officinalis suggests that past colonization history, coupled with genetic drift within local populations, rather than contemporary gene flow are the major forces shaping variation within this species.     

Phylogeographic evidence for the postglacial colonization of the North and Baltic Sea coasts from inland glacial refugia by Triglochin maritima L.

We investigated the geographical distribution of genetic variation in 67 individuals of Triglochin maritima from 38 localities across Europe using AFLP markers. Analysis of genetic variation resulted in the recognition of two major genetic groups. Apart from few geographical outliers, these are distributed (1) along the Atlantic coasts of Portugal, Spain and France and (2) in the North Sea area, the Baltic Sea area, at central European inland localities, the northern Adriatic Sea coast and the Mediterranean coast of southwest France. Considering possible range shifts of T. maritima in reaction to Quaternary climatic changes as deduced from the present-day northern temperature limit of the species, Quaternary changes of coastline in the North Sea area and the very recent origin of the Baltic Sea, we conclude that the coastal populations of T. maritima in the North Sea and Baltic Sea areas originated from inland populations.     

Control region haplotype variation in the central Mediterranean common sole indicates geographical isolation and population structuring in Italian stocks

Nine common sole Solea vulgaris samples from five fishery Management Units (MUs) of the central Mediterranean exhibited differences in a control region sequence marker. Parsimony network showed two low-divergent haplotype phylogroups. The former predominated in samples from the Ligurian, Tyrrhenian, Adriatic and north-western Ionian Seas (MUs 9, 10, 17, 18 and northern part of 19, respectively), whereas the latter was abundant in the south-western Ionian (southern part of MU 19). The geographical clustering of maternal lineages accounts for population structuring and indicates geographical isolation of common sole stocks in the Mediterranean. Several life-history traits of common sole and hydro-geographical features of the basin might support this pattern of differentiation. Haplotype frequency differences were detected among samples either from different MUs or within the same MU. This situation indicates the presence of partially subdivided or nearly panmictic population units, whose exploitation should be based on appropriate bio-ecological features. The usefulness of control region sequence marker enhances the routinely use of the genetic stock structure analysis in low-dispersal demersal marine resources.     

Life history and genetic characterisation of sea trout Salmo trutta in the Adriatic Sea

1. The brown trout Salmo trutta is characterised by both anadromous (sea trout) and resident populations, naturally occurring in Atlantic and Ponto-Caspian rivers. Sea trout are currently considered absent from rivers of the Mediterranean area, probably because of the non-optimal chemical–physical characteristics of the Mediterranean Sea. However, the occasional bycatch of smoltified S. trutta in the Adriatic Sea is well known among fishermen and the biological explanation of this phenomenon is still controversial. The aim of this study was to compare the genetic diversity of freshwater and marine brown trout to try to understand the factors underlying the presence of putative anadromous brown trout in the Adriatic Sea.
2. In this study, we analysed the genetic diversity of: (1) wild brown trout collected from the Esino River (central Italy); (2) a domestic strain of brown trout used for stocking the study area; and (3) a sample of Adriatic sea trout collected near the outlet of the Esino River. Together with genetic analysis, we carried out scale analysis in order to track the freshwater/marine stages of the life cycle in the sea trout samples. The genetic characterisation was carried out by polymerase chain reaction–restriction fragment length polymorphism analysis of the mtDNA fragment ND-5/6 and the nuclear locus LDH-C1* and by genotyping 15 microsatellite loci. The genetic polymorphism obtained was used to investigate intra- and inter-population genetic diversity, rates of genetic introgression between wild and domestic samples and the origin of sea trout specimens by using assignment tests.
3. Our genetic analyses demonstrated that the sea trout analysed in this study are from the domestic strain of Atlantic origin used in central Italy for stocking activities. The level of genetic introgression between native and domestic samples is high in the Esino River. The populations more resilient to introgressive hybridisation appeared to be those living in the portion of the river network dominated by carbonate rocks. Assignment tests (GeneClass) suggest the existence of a link between stocking efforts and the freshwater origin of the sea trout. In addition, data obtained from the analysis of scales, size measurement, and sex determination showed a pattern of smolt age, size, and sex ratio very similar to those observed in other anadromous populations.
4. In conclusion, the present study highlighted that sea trout from the central Adriatic Sea originated from brown trout of Atlantic origin inhabiting the Esino River. Their seaward migratory behaviour could represent a consequence of an active migration instead of a passive displacement by water flow. Our results also showed that traditional stocking practices represent a negative activity for the conservation of the last Mediterranean native S. trutta populations.

     

Allozyme, mitochondrial-DNA, and morphometrie variability indicate cryptic species of anchovy (Engraulis encrasicolus)

Previous surveys of population structure in the Atlantic-Mediterranean anchovy Engraulis encrasicolus L. have reported heterogeneity in morphology, allozyme frequencies, and mitochondrial DNA haplotype frequencies at a regional scale. In particular, two stocks of anchovy have been detected in the Adriatic Sea. In this paper, the available data is reviewed with the aim to relate genetic variation to geography at the widest possible geographical scale, for investigating the evolutionary mechanisms underlying stock structure in anchovy. Correspondence analysis of allozyme frequencies (24 samples, three polymorphic loci) compiled from the literature indicates three distinct entities in the Mediterranean Sea. Open-sea or oceanic anchovy populations are genetically different from inshore-water populations within a region (Nei's ^ G _ST = 0.035–0.067), while broadscale geographical variation is weak for each of these two habitat-specific forms (^ G _ST = 0.005–0.006). Mitochondrial-DNA haplotype frequencies support the distinction between an inshore form and an oceanic form (^ G _ST = 0.067–0.107), with virtually no genetic differences among oceanic populations across the Gulf of Biscay, the western Mediterranean and the Ionian Sea (^ G _ST = −0.001). If natural selection on marker loci is unimportant, these results indicate the occurrence of two parapatric, genetically distinct, habitat-specific forms that are widely distributed throughout the Mediterranean Sea. Persistent allele and haplotype-frequency differences between these forms indicate reproductive isolation and the presence of an E. encrasicolus species complex in the Mediterranean. © 2002 The Linnean Society of London, Biological Journal of the Linnean Society 75 : 261–269.     

Phylogeography of isolated freshwater three-spined stickleback Gasterosteus aculeatus populations in the Adriatic Sea basin

Analyses of mitochondrial (mt) DNA and microsatellite variation were carried out to examine the relationships between 10 freshwater populations of three-spined sticklebacks Gasterosteus aculeatus along the eastern coast of the Adriatic Sea. Partial sequences of the mtDNA control region and cytochrome b gene, in addition to 15 microsatellite loci, were used to analyse populations from four isolated river catchments. Results uncovered an Adriatic lineage that was clearly divergent from the European lineage, and confirmed that the most divergent and ancient populations are located within the Adriatic lineage as compared with other European populations. Two northern Adriatic populations formed independent clades within the European mitochondrial lineage, suggesting different colonization histories of the different Adriatic populations. Nuclear marker analyses also indicated deep divergence between Adriatic and European populations, albeit with some discordance between the mtDNA phylogeny of the northern Adriatic populations, further highlighting the strong differentiation among the Adriatic populations. The southern populations within the Adriatic lineage were further organized into distinct clades corresponding to respective river catchments and sub-clades corresponding to river tributaries, reflecting a high degree of population structuring within a small geographic region, concurrent with suggestions of existence of several microrefugia within the Balkan Peninsula. The highly divergent clades and haplotypes unique to the southern Adriatic populations further suggest, in accordance with an earlier, more limited survey, that southern Adriatic populations represent an important reservoir for ancient genetic diversity of G. aculeatus.     

Successful worldwide invasion of the veined rapa whelk, <Emphasis Type="Italic">Rapana venosa</Emphasis>, despite a dramatic genetic bottleneck

Biological invasions represent an important component of global change, with potentially huge detrimental effects on native biological biodiversity and ecosystems. Knowledge about invasion history provides information about the invasion process and the origin and genetic composition of invading populations. To clarify the source and invasive routes of a successful world-wide invader, the veined rapa whelk, Rapana venosa, genetic variability of samples from five representative native populations from coasts of Japan and China and 13 worldwide invasive populations was analyzed using 11 nuclear microsatellite loci. A dramatic decrease of genetic variation was detected in the invasive populations compared with the native populations. The results demonstrated that R. venosa was capable of establishing itself in many areas despite a dramatic genetic bottleneck, suggesting that a remarkable reduction of genetic diversity is not a limiting factor for short-term success of this invasive species. Considering the lack of mitochondrial variation previously observed in the invasive populations, the dramatic genetic bottleneck and the allele distribution detected using microsatellites suggested that the original introduced Black Sea population could have been founded by very few individuals, perhaps only a single female and a single male. The initial invasive Black Sea population was likely an accidental introduction from Japan, and then invaded the Adriatic Sea by range expansion, which served as a source for subsequent invasive populations in Europe and America by various transport vectors. In addition, microsatellite alleles in the invasive populations showed a tendency to mutate with the addition or deletion of a single repeat, which is consistent with the stepwise mutation model. Our findings provide a good example of how an aquatic invader with a drastic genetic bottleneck and very low genetic diversity rapidly expands its geographical range.     

First Microsatellite Loci of Red Mullet (<Emphasis Type="Italic">Mullus barbatus</Emphasis>) and Their Application to Genetic Structure Analysis of Adriatic Shared Stock

In order to study the genetic structure of the Adriatic shared stock of red mullet (Mullus barbatus), we developed a set of dinucleotide microsatellite markers. A dinucleotide-enriched genomic library was obtained, and 6 polymorphic dinucleotide loci were successfully optimized. The markers showed high expected heterozygosity (from 0.68 to 0.92) and allele number (from 12 to 33); thus they appear to be suitable for detecting genetic differences in the population of red mullet. Four Adriatic samples were subsequently analyzed for microsatellite variation, and the results showed subtle but statistically significant genetic differentiation, indicating that the Adriatic red mullet may group into local, genetically isolated populations. No correlation between geographic distance and genetic differentiation was observed. In addition, the evidence of recent bottlenecks in the Adriatic samples indicates that the observed population subdivision might reflect random local allelic variations, generated by reproductive success, survival rates, or fishing pressure.     

Mitochondrial DNA Diversity of the Asian Moon Scallop, <Emphasis Type="Italic">Amusium pleuronectes</Emphasis> (Pectinidae), in Thailand

Sequence variation of the mitochondrial DNA 16S rRNA region of the Asian moon scallop, Amusium pleuronectes, was surveyed in seven populations along the coast of Thailand. A total of 16 unique haplotypes were detected among 174 individuals with a total 27 variable sites out of 534 bp sequenced. The mitochondrial haplotypes grouped into two distinct arrays (estimated to differ by about 2.62% to 2.99% nucleotide divergence) that characterized samples collected from the Gulf of Thailand versus the Andaman Sea. Low levels of intrapopulation variation were observed, while in contrast, significant divergence was observed between populations from the Gulf of Thailand and Andaman Sea. Results of AMOVA reveal a high F _ST value (0.765) and showed that the majority of the total genetic variance (76.03%) occurred among groups (i.e., Andaman Sea and the Gulf of Thailand) and little among populations within the group (0.52%) and within populations (23.45%). The genetic differentiation between the populations recorded in the present study is similar to that observed in a variety of marine species in the Indo-Pacific. The implications of the findings for management of A. pleuronectes genetic resources in Thailand are discussed.     

Mitochondrial DNA Reveals Genetic Structuring of Pinna nobilis across the Mediterranean Sea

Pinna nobilis is the largest endemic Mediterranean marine bivalve. During past centuries, various human activities have promoted the regression of its populations. As a consequence of stringent standards of protection, demographic expansions are currently reported in many sites. The aim of this study was to provide the first large broad-scale insight into the genetic variability of P. nobilis in the area that encompasses the western Mediterranean, Ionian Sea, and Adriatic Sea marine ecoregions. To accomplish this objective twenty-five populations from this area were surveyed using two mitochondrial DNA markers (COI and 16S). Our dataset was then merged with those obtained in other studies for the Aegean and Tunisian populations (eastern Mediterranean), and statistical analyses (Bayesian model-based clustering, median-joining network, AMOVA, mismatch distribution, Tajima’s and Fu’s neutrality tests and Bayesian skyline plots) were performed. The results revealed genetic divergence among three distinguishable areas: (1) western Mediterranean and Ionian Sea; (2) Adriatic Sea; and (3) Aegean Sea and Tunisian coastal areas. From a conservational point of view, populations from the three genetically divergent groups found may be considered as different management units.     

Phylogeography of the planktonic chaetognath Sagitta setosa reveals isolation in European seas

Numerous planktonic species have disjunct distribution patterns in the world's oceans. However, it is unclear whether these are truly unconnected by gene flow, or whether they are composed of morphologically cryptic species. The marine planktonic chaetognath Sagitta setosa Müller has a discontinuous geographic distribution over the continental shelf in the northeastern Atlantic, Mediterranean Sea, and Black Sea. Morphological variation between these populations has been described, but overlaps and is therefore unsuitable to determine the degree of isolation between populations. To test whether disjunct populations are also genetically disjunct, we sequenced a 504-bp fragment of mitochondrial DNA comprising the cytochrome oxidase II region of 86 individuals. Sequences were highly variable; each represented a different haplotype. Within S. setosa, sequence divergence ranged from 0.2 to 8.1% and strong phylogeographic structure was found, with four main groups corresponding to the northeastern Atlantic, Mediterranean Sea (including Ligurian Sea, Tyrrhenian Sea and Gulf of Gabes), Adriatic Sea, and Black Sea. Two of these (Atlantic and Black Sea) were resolved as monophyletic clades, thus gene flow between disjunct populations of S. setosa has been extremely limited and lineage sorting has taken place. The deepest divergence was between Atlantic and Mediterranean/Black Sea populations followed by a split between Mediterranean and Black Sea populations. The Mediterranean/Black Sea clade comprised three groups, with the Adriatic Sea as the most likely sister clade of the Black Sea. These data are consistent with a colonization of the Black Sea from the Mediterranean. Furthermore, a possible cryptic species was found in the Black Sea with 23.1% sequence divergence from S. setosa. Two possibilities for the evolutionary origin of this species are proposed, namely, that it represents a relict species from the ancient Paratethys, or that it represents another chaetognath species that colonized the Black Sea more recently. Even though the exact timing of disjunction of S. setosa populations remains unclear, on the basis of the geological and paleoclimatic history of the European basins and our estimates of net nucleotide divergence, we suggest that disjunct populations arose through vicariance resulting from the cyclical changes in temperature and sea levels during the Pleistocene. We conclude that these populations have remained disjunct, not because of limited dispersal ability, but because of the inability to maintain viable populations in suboptimal, geographically intermediate areas.