Phylogenetic differentiation between Atlantic <Emphasis Type="Italic">Scomber colias</Emphasis> and Pacific <Emphasis Type="Italic">Scomber japonicus</Emphasis> based on nuclear DNA sequences

In the classical taxonomy, three Scomber species are distinguished: S. scombrus, S. australasicus, and S. japonicus. Yet, some fish taxonomists have recently recognized Scomber colias, inhabiting the Atlantic Ocean, as a separate species from S. japonicus, distributed in the Pacific Ocean. Such proposal was based on significant mitochondrial DNA divergence as well as great phenotypic variation among individuals from these two ocean basins. However, in the absence of nuclear DNA data this issue remains still controversial. In this study, a phylogenetic analysis of nuclear 5S rDNA sequences was performed. A total of 30 individuals of S. colias collected in the Atlantic and 34 specimens of S. japonicus from the Pacific were characterized. Moreover, nine individuals of Pacific S. australasicus and eight of Atlantic S. scombrus were included. Maximum likelihood, maximum parsimony, and neighbor-joining analyses revealed the presence of two well-supported distinct clades corresponding to S. colias and S. japonicus, respectively. Altogether, morphologic and genetic data are in agreement with the recognition of two different species, S. colias in the Atlantic, and S. japonicus in the Pacific.     

Differential population structuring of two closely related fish species, the mackerel (Scomber scombrus) and the chub mackerel (Scomber japonicus), in the Mediterranean Sea

Population genetic structures of the mackerel (Scomber scombrus) and chub mackerel (Scomber japonicus) were studied in the Mediterranean Sea. Fragments of 272 bp (S. scomber) and 387 bp (S. japonicus) of the 5'-end of the mitochondrial control region were sequenced from spawning individuals collected off the coasts of Greece, Italy, Spain, and Portugal. High levels of mitochondrial control region haplotypic diversity (> 0.98) were found for both Scomber species. Nucleotide diversity was higher in the mackerel (0.022) than in the chub mackerel (0.017). Global F(ST) values were also higher and significant in the mackerel (0.024, P < 0.0001) as opposed to the chub mackerel (0.003, P > 0.05). Molecular variance analyses showed differential genetic structuring for these two closely related species. There is extensive gene flow between Mediterranean Sea and Atlantic Ocean populations of chub mackerel, which are organized into a larger panmictic unit. In contrast, Mediterranean Sea populations of mackerel show some degree of genetic differentiation and are structured along an east-west axis. The analysed eastern Mediterranean Sea mackerel populations (Greece, Italy) are clearly separated from that of the western Mediterranean Sea (Barcelona), which forms a panmictic unit with eastern Atlantic Ocean populations. The genetic structures of both species showed asymmetric migration patterns and indicated population expansion.     

Phylogeography of the green turtle, Chelonia mydas, in the Southwest Indian Ocean

Patterns of mitochondrial DNA (mtDNA) variation were used to analyse the population genetic structure of southwestern Indian Ocean green turtle (Chelonia mydas) populations. Analysis of sequence variation over 396 bp of the mtDNA control region revealed seven haplotypes among 288 individuals from 10 nesting sites in the Southwest Indian Ocean. This is the first time that Atlantic Ocean haplotypes have been recorded among any Indo-Pacific nesting populations. Previous studies indicated that the Cape of Good Hope was a major biogeographical barrier between the Atlantic and Indian Oceans because evidence for gene flow in the last 1.5 million years has yet to emerge. This study, by sampling localities adjacent to this barrier, demonstrates that recent gene flow has occurred from the Atlantic Ocean into the Indian Ocean via the Cape of Good Hope. We also found compelling genetic evidence that green turtles nesting at the rookeries of the South Mozambique Channel (SMC) and those nesting in the North Mozambique Channel (NMC) belong to separate genetic stocks. Furthermore, the SMC could be subdivided in two different genetic stocks, one in Europa and the other one in Juan de Nova. We suggest that this particular genetic pattern along the Mozambique Channel is attributable to a recent colonization from the Atlantic Ocean and is maintained by oceanic conditions in the northern and southern Mozambique Channel that influence early stages in the green turtle life cycle.     

Hierarchical analyses of genetic variation of samples from breeding and feeding grounds confirm the genetic partitioning of northwest Atlantic and South Atlantic populations of swordfish (Xiphias gladius L.)

In species with high migratory potential, the genetic signal revealing population differentiation is often obscured by population admixture. To our knowledge, the explicit comparison of genetic samples from known spawning and feeding areas has not been conducted for any highly migratory pelagic fish species. This study examines the geographic heterogeneity of swordfish mitochondrial DNA (mtDNA) lineages within the Atlantic Ocean using 330 base pairs of sequence of the control region from 480 individuals. Hierarchical analyses of sequence variation were conducted to test whether samples from areas identified as the corresponding spawning and feeding grounds for the northwest (NW) Atlantic (Caribbean and Georges Banks-US northeast) and the South Atlantic (Brazil-Uruguay and Gulf of Guinea), were more closely related to each other than to samples from any other region, including the Mediterranean Sea, the Indian Ocean, and the Pacific Ocean. Phylogeographic analyses reveal that swordfish mtDNA phylogeny is characterized by incomplete lineage sorting and secondary contact of two highly divergent clades. However, despite this complex phylogenetic signature, results from an analysis of nucleotide diversity and from an analysis of molecular variance (AMOVA) were for the most part concordant and indicate that NW Atlantic and South Atlantic swordfish belong to separate populations. The mtDNA distinctiveness of NW Atlantic and South Atlantic swordfish populations is indicative of philopatric behavior in swordfish towards breeding and feeding areas.     

The curious case of Hermodice carunculata (Annelida: Amphinomidae): evidence for genetic homogeneity throughout the Atlantic Ocean and adjacent basins

Over the last few decades, advances in molecular techniques have led to the detection of strong geographic population structure and cryptic speciation in many benthic marine taxa, even those with long‐lived pelagic larval stages. Polychaete annelids, in particular, generally show a high degree of population divergence, especially in mitochondrial genes. Rarely have molecular studies confirmed the presence of ‘cosmopolitan’ species. The amphinomid polychaete Hermodice carunculata was long considered the sole species within its genus, with a reported distribution throughout the Atlantic and adjacent basins. However, recent studies have indicated morphological differences, primarily in the number of branchial filaments, between the East and West Atlantic populations; these differences were invoked to re‐instate Hermodice nigrolineata, formerly considered a junior synonym of H. carunculata. We utilized sequence data from two mitochondrial (cytochrome c oxidase subunit I, 16S rDNA) markers and one nuclear (internal transcribed spacer) marker to examine the genetic diversity of Hermodice throughout its distribution range in the Atlantic Ocean, including the Mediterranean Sea, the Caribbean Sea, the Gulf of Mexico and the Gulf of Guinea. Our analyses revealed generally low genetic divergences among collecting localities and between the East and West Atlantic, although phylogenetic trees based on mitochondrial data indicate the presence of a private lineage in the Mediterranean Sea. A re‐evaluation of the number of branchial filaments confirmed differences between East and West Atlantic populations; however, the differences were not diagnostic and did not reflect the observed genetic population structure. Rather, we suspect that the number of branchial filaments is a function of oxygen saturation in the environment. Our results do not support the distinction between H. carunculata in the West Atlantic and H. nigrolineata in the East Atlantic. Instead, they re‐affirm the older notion that H. carunculata is a cohesive species with a broad distribution across the Atlantic Ocean.     

Fine Spatial Structure of Atlantic Hake (<Emphasis Type="Italic">Merluccius merluccius</Emphasis>) Stocks Revealed by Variation at Microsatellite Loci

Genetic variation at 5 microsatellite loci was analyzed for European hake Merluccius merluccius sampled from 9 different regions in the Atlantic Ocean and the Mediterranean Sea. Significant genetic differentiation was found between samples, suggesting a fine subdivision of Atlantic and Mediterranean hake stocks. These results are discussed in the context of the decline of demersal fish species, probably due to overfishing.     

短沟对虾两个野生群体遗传多样性的RAPD分析

利用RAPD标记技术检测了厦门和汕头沿海2个短沟对虾群体基因组DNA的多态性,并对其遗传多样性进行了分析。从40条随机引物中筛选出13个10bp引物。共扩增出65条清晰可重复的DNA片段,片断长度为100—2200bp,在2个群体间没有检测到特异的片段。厦门和汕头群体的多态片段比例分别为87．69％和89．23％,杂合度分别为0．212和0．218,遗传多样性指数分别为0．2847和0．2913,两群体间的遗传距离为0．018,FST值为0．004。可见两野生群体种质资源仍然维持在良好水平,遗传分化程度很低,可能是同一种群,具有进一步开发的潜力。     

New Nuclear SNP Markers Unravel the Genetic Structure and Effective Population Size of Albacore Tuna (Thunnus alalunga)

In the present study we have investigated the population genetic structure of albacore (Thunnus alalunga, Bonnaterre 1788) and assessed the loss of genetic diversity, likely due to overfishing, of albacore population in the North Atlantic Ocean. For this purpose, 1,331 individuals from 26 worldwide locations were analyzed by genotyping 75 novel nuclear SNPs. Our results indicated the existence of four genetically homogeneous populations delimited within the Mediterranean Sea, the Atlantic Ocean, the Indian Ocean and the Pacific Ocean. Current definition of stocks allows the sustainable management of albacore since no stock includes more than one genetic entity. In addition, short- and long-term effective population sizes were estimated for the North Atlantic Ocean albacore population, and results showed no historical decline for this population. Therefore, the genetic diversity and, consequently, the adaptive potential of this population have not been significantly affected by overfishing.     

Genetic identification of Anisakis larvae in European hake from Atlantic and Mediterranean waters for stock recognition

The occurrence of seven species of the larval parasitic nematode Anisakis , which can be used as a biological tag for hake Merluccius merluccius stocks throughout their geographical range, is reported. Hake were collected from 14 localities in the Mediterranean Sea and the Atlantic Ocean. Anisakis larvae ( n  = 1950), which were recovered, were identified to species by means of genetic markers (allozymes). Within Anisakis type I, the larvae of A. pegreffii , A. simplex s.s ., A. typica and A. ziphidarum were detected, while within Anisakis type II, A. physeteris , A. brevispiculata and Anisakis sp. were identified. There were significant differences in the relative proportions of the various Anisakis species identified in hake samples from the Mediterranean Sea and Atlantic Ocean, suggesting the existence of different stocks of M. merluccius in European waters.     

The genetic population structure of Thunnus thynnus (Linnaeus, 1758) in the Mediterranean Sea,a controversial issue

The Atlantic bluefin tuna (ABFT), Thunnus thynnus (Linnaeus, 1758), is an important commercial species managed as two different stocks, western and eastern Atlantic, with their spawning grounds in the Gulf of Mexico and in the Mediterranean Sea, respectively. The eastern Atlantic stock has been overexploited in the last decades, leading to the application of specific management measures introduced by the International Commission for the Atlantic Tuna (ICCAT). A clear understanding of the genetic structure of ABFT Mediterranean population should be pursued in order to support management decisions. To date the genetic studies on the Mediterranean ABFT, carried out with different molecular markers and sampling procedures, have produced unclear results. Here, we analysed ABFT samples from central and western Mediterranean Sea with mitochondrial sequences and 11 microsatellite loci to investigate, among the others, the area of the Strait of Messina, where environmental conditions seem to support a resident population of ABFT. Furthermore, genetic analyses of mitochondrial sequences were carried out including nucleotide sequences of Adriatic ABFT wild larvae retrieved from GenBank. Among the investigated areas a genetic differentiation was detected between the Strait of Messina and the Tyrrhenian Sea with microsatellite loci according to the exact G test, but not to the Bayesian analyses carried out with STRUCTURE. The analyses with mitochondrial sequences do not reveal any differentiation among sampled areas, however, a highly significant genetic divergence was observed between the Adriatic mitochondrial sequences retrieved from GenBank and the central‐western Mediterranean sequences obtained in the present work. Our results provide some evidence of population structure of Mediterranean ABFT adding pieces to a still unclear picture.     

Genetic structure and environmental heterogeneity in the European hake (Merluccius merluccius)

This study aimed at assessing the genetic structure and the state of the stocks of the European hake (Merluccius merluccius). To this end, 15 samples were taken from the whole range of the species and analysed using allozymes. Since 11 samples were taken from the poorly studied Mediterranean Sea, the results obtained provided a complete picture of the hake's genetic structure and an initial insight into its relationships with environmental features. Atlantic and Mediterranean hake populations are separated by the Almeria-Oran front. This area has been proved to be the boundary between Atlantic and Mediterranean stocks of many marine organisms, but some doubt exists concerning the efficaciousness of the local gyres as barriers to the gene flow. Our data have evidenced a latitudinal cline at loci Gapdh and Gpi-2 within the Mediterranean Sea, with a further steep change across the Almeria-Oran front. The genetic pattern showed a strong correlation with the values of the salinity both at the surface and at -320 m and of the salinity + temperature at the surface, suggesting a role for these parameters in maintaining the genetic differentiation among the two population groups through selective processes. Finally, the levels of genetic variability were found to be slightly lower in the depleted Atlantic stock than in the Mediterranean one.     

Review of albacore tuna, <Emphasis Type="Italic">Thunnus alalunga</Emphasis>, biology,fisheries and management

Albacore is one of the most important commercially harvested species in the world’s oceans. Despite a long history of scientific research, there is no global review or synthesis of knowledge about the species across all oceanic regions. We analysed 613 published studies that report on albacore and summarize the current state of knowledge on biology, stock structure, fisheries and management. To describe the trends in albacore fisheries, we examined the catch and effort databases of Regional Fisheries Management Organisations. The stocks of albacore are generally largest in the Pacific Ocean and smallest in the Mediterranean Sea. The biology of Atlantic and Pacific Ocean stocks are well documented, while the Indian Ocean and the Mediterranean Sea stocks are more data deficient. These two latter areas should be considered as priorities for future research and data collection in order to better understand the state of global stocks of albacore tuna. Improved information would also assist with delineating stock boundaries needed for sustainable management of this species.     

Population genetic structure of North Atlantic, Mediterranean Sea and Sea of Cortez fin whales, Balaenoptera physalus (Linnaeus 1758): analysis of mitochondrial and nuclear loci

Samples were collected from 407 fin whales, Balaenoptera physalus , at four North Atlantic and one Mediterranean Sea summer feeding area as well as the Sea of Cortez in the Pacific Ocean. For each sample, the sex, the sequence of the first 288 nucleotides of the mitochondrial (mt) control region and the genotype at six microsatellite loci were determined. A significant degree of divergence was detected at all nuclear and mt loci between North Atlantic/Mediterranean Sea and the Sea of Cortez. However, the divergence time estimated from the mt sequences was substantially lower than the time elapsed since the rise of the Panama Isthmus, suggesting occasional gene flow between the North Pacific and North Atlantic ocean after the separation of the two oceans. Within the North Atlantic and Mediterranean Sea, significant levels of heterogeneity were observed in the mtDNA between the Mediterranean Sea, the eastern (Spain) and the western (the Gulf of Maine and the Gulf of St Lawrence) North Atlantic. Samples collected off West Greenland and Iceland could not be unequivocally assigned to either of the two areas. The homogeneity tests performed using the nuclear data revealed significant levels of divergence only between the Mediterranean Sea and the Gulf of St Lawrence or West Greenland. In conclusion, our results suggest the existence of several recently diverged populations in the North Atlantic and Mediterranean Sea, possibly with some limited gene flow between adjacent populations, a population structure which is consistent with earlier population models proposed by Kellogg, Ingebrigtsen, and Sergeant.     

Phylogeography and population structure of the Atlantic and Mediterranean green turtle Chelonia mydas: a mitochondrial DNA control region sequence assessment

Mitochondrial (mt) DNA sequences were analysed to resolve the phylogeography and population genetic structure of Atlantic and Mediterranean populations of green turtles ( Chelonia mydas ). Analysis of sequence variation over 487 base pairs of the control (D-loop) region identified 18 haplotypes among 147 individuals from nine nesting populations. Pairwise comparisons of haplotype frequencies distinguished most nesting colonies, indicating significant genetic differentiation among rookeries and a strong propensity for natal homing behaviour by nesting females. Comparison of control region sequence data to earlier restriction fragment length polymorphism (RFLP) data for the same individuals demonstrates approximately a sixfold higher substitution rate in the 5' end of the control region. The sequence data provide higher resolution both in terms of the number of mtDNA genotype variants and the phylogeographic relationships detected within the Atlantic region, and reveal a gene genealogy that distinguishes two groups of haplotypes corresponding to (i) the western Caribbean and Mediterranean, and (ii) eastern Caribbean, South Atlantic and West Africa. The data suggest that phylogeographic patterns in the Atlantic Ocean may be interpreted in terms of female nest site fidelity and episodic dispersal events. The distribution of mtDNA haplotypes within the region is thus explained by the geological and climatic alternations (glacial and interglacial) over the last million years.     

Phylogeography of the Calonectris shearwaters using molecular and morphometric data

We investigated phylogenetic relationships and the biogeographic history of the Calonectris species complex, using both molecular and biometric data from one population of the Cape Verde shearwater Calonectris edwardsii (Cape Verde Islands), one from the streaked shearwater C. leucomelas (western Pacific Ocean) and 26 from Cory's shearwater populations distributed across the Atlantic (C. d. borealis) and the Mediterranean (C. d. diomedea). The streaked shearwater appeared as the most basal and distant clades, whereas the genetic divergences among the three main clades within the Palearctic were similar. Clock calibrations match the first speciation event within Calonectris to the Panama Isthmus formation, suggesting a vicariant scenario for the divergence of the Pacific and the Palearctic clades. The separation between the Atlantic and Mediterranean clades would have occurred in allopatry by range contraction followed by local adaptation during the major biogeographic events of the Pleistocene. The endemic form from Cape Verde probably evolved as a result of ecological divergence from the Mediterranean subspecies. Finally, one Mediterranean population (Almeria) was unexpectedly grouped into the Atlantic subspecies clade, both by genetic and by morphometric analyses, pointing out the Almeria-Oran oceanographic front (AOOF) as the actual divide between the two Cory's shearwater subspecies. Our results highlight the importance of oceanographic boundaries as potentially effective barriers shaping population and species phylogeographical structure in pelagic seabirds.     

Octopus vulgaris (Cuvier, 1797) in the Mediterranean Sea: Genetic Diversity and Population Structure

The common octopus, Octopus vulgaris Cuvier 1797, is a largely exploited cephalopod species in the Mediterranean Sea and the Atlantic Ocean, as well as along the coasts of Africa, Brazil and Japan, where its taxonomic identity is still debated. The assessment of its genetic structure is a pressing need to correctly manage the resource and to avoid overfishing and collapsing of local stocks. Here we analysed genetic variation and population structure of O. vulgaris using thirteen microsatellite loci in seven sampling localities from the Mediterranean Sea and one from the Atlantic Ocean. We also used a DNA barcoding approach by COI gene fragment to understand the phylogenetic relationships among the specimens here investigated and the ones whose sequences are available in literature. Our results reveal high levels of allelic richness and moderate heterozygosity in all samples investigated, and a pronounced differentiation of the Atlantic and Sicilian specimens. This latter aspect seems to support the isolation of the biota within the Strait of Messina. A certain degree of differentiation was detected among the other geographic samples within the Mediterranean Sea, which is more compatible with an island model than isolation by distance. The occurrence of null alleles affected more genetic diversity indices than population structure estimations. This study provides new insights about the genetic diversity and structure of O. vulgaris in the area of interest, which can be used as guidelines for a fisheries management perspective.     

Range‐wide population structure of European sea bass Dicentrarchus labrax

The euryhaline European sea bass Dicentrarchus labrax L., inhabiting the coasts of the eastern Atlantic Ocean and Mediterranean Sea, has had many opportunities for differentiation throughout its large natural range. However, evidence for this has been incompletely documented geographically and with an insufficient number of markers. Therefore, its full range was sampled at 22 sites and individuals were genotyped with a suite of mapped markers, including 14 microsatellite loci (N = 536) and 46 neutral or gene‐linked single nucleotide polymorphisms (SNPs; N = 644). We confirm that the Atlantic and Mediterranean basins harbour two distinct lineages. Within the Atlantic Ocean no pattern was obvious based on the microsatellite and SNP genotypes, except for a subtle difference between South‐eastern and North‐eastern Atlantic sea bass attributed to limited introgression of alleles of Mediterranean origin. SNP genotypes of the Mediterranean lineage differentiated into three groups, probably under the influence of geographical isolation. The Western Mediterranean group showed genetic homogeneity without evidence for outlier loci. The Adriatic group appeared as a distinct unit. The Eastern Mediterranean group showed a longitudinal gradient of genotypes and most interestingly an outlier locus linked to the somatolactin gene. Overall, the spatial pattern fits those observed with other taxa of between‐basin segregation and within‐basin connectivity, which concurs well with the swimming capabilities of European sea bass. Evidence from a few outlier loci in this and other studies encourages further exploration of its regional connectivity and adaptive evolution.     

Genetic differentiation and secondary contact zone in the seagrass Cymodocea nodosa across the Mediterranean–Atlantic transition region

Aim A central question in evolutionary ecology is the nature of environmental barriers that can limit gene flow and induce population genetic divergence, a first step towards speciation. Here we study the geographical barrier constituted by the transition zone between the Atlantic Ocean and the Mediterranean Sea, using as our model Cymodocea nodosa, a seagrass distributed throughout the Mediterranean and in the Atlantic, from central Portugal to Mauritania. We also test predictions about the genetic footprints of Pleistocene glaciations. Location The Atlantic–Mediterranean transition region and adjacent areas in the Atlantic (Mauritania to south‐west Portugal) and the Mediterranean. Methods We used eight microsatellite markers to compare 20 seagrass meadows in the Atlantic and 27 meadows in the Mediterranean, focusing on the transition between these basins. Results Populations from these two regions form coherent groups containing several unique, high‐frequency alleles for the Atlantic and for the Mediterranean, with some admixture west of the Almeria–Oran Front (Portugal, south‐west Spain and Morocco). These are populations where only one or a few genotypes were found, for all but Cadiz, but remarkably still show the footprint of a contact zone. This extremely low genotypic richness at the Atlantic northern edge contrasts with the high values (low clonality) at the Atlantic southern edge and in most of the Mediterranean. The most divergent populations are those at the higher temperature range limits: the southernmost Atlantic populations and the easternmost Mediterranean, both potential footprints of vicariance. Main conclusions A biogeographical transition region occurs close to the Almeria–Oran front. A secondary contact zone in Atlantic Iberia and Morocco results from two distinct dispersal sources: the Mediterranean and southernmost Atlantic populations, possibly during warmer interglacial or post‐glacial periods. The presence of high‐frequency diagnostic alleles in present‐day disjunct populations from the southernmost Atlantic region indicates that their separation from all remaining populations is ancient, and suggests an old, stable rear edge.     

Genetic divergence between Atlantic and Indo-Pacific stocks of bigeye tuna (Thunnus obesus) and admixture around South Africa

Two mitochondrial DNA segments of the bigeye tuna (Thunnus obesus) were amplified by polymerase chain reaction (PCR), and restriction fragment length polymorphism (RFLP) analyses of these segments were used for the genetic stock study. The variation in a segment flanking the ATPase and COIII genes was low; only two genotypes (α and β) were detected by RsaI digestion. Yet a large difference in the genotype distribution was observed between ocean basin samples. The α type predominated in four Atlantic samples, where 178 of 244 individuals were the α type. In contrast, only one of 195 individuals collected in the Indo‐Pacific was the α type? The frequency of the α type varied considerably from 0 to 80% among seven samples collected off the Cape of Good Hope. The variation found in the other segment, containing the D‐loop region, was much higher; two endonucleases (DpnII and RsaI) detected five genotypes each and 15 composite genotypes. A highly significant difference in genotype frequencies was observed between the Atlantic and Indo‐Pacific samples, but no heterogeneity was observed among the four Atlantic or among four Indo‐Pacific samples. These results clearly indicate that not only gene flow, but also fish migration, between the Atlantic and Indian Oceans are severely restricted, and that fishes from these distinct stocks are intermingling around South Africa. The simple and diagnostic genetic marker found in this study can be used to estimate mixing ratios between Atlantic and Indian stocks around South Africa.     

Multilocus phylogenetic analysis of the genus Atherina (Pisces: Atherinidae)

Sand-smelts are small fishes inhabiting inshore, brackish and freshwater environments and with a distribution in the eastern Atlantic and Mediterranean Sea, extending south into the Indian Ocean. Here, we present a broad phylogenetic analysis of the genus Atherina using three mitochondrial (control region, 12S and 16S) and two nuclear markers (rhodopsin and 2nd intron of S7). Phylogenetic analyses fully support the monophyly of the genus. Two anti-tropical clades were identified, separating the South African Atherina breviceps from the north-eastern Atlantic and Mediterranean Atherina' species. In European waters, two groups were found. The first clade formed by a well supported species-pair: Atherina presbyter (eastern Atlantic) and Atherina hepsetus (Mediterranean), both living in marine waters; a second clade included Atherina boyeri (brackish and freshwater environments) and two independent lineages of marine punctated and non-punctated fishes, recently proposed as separate species. Sequence divergence values strongly suggest multiple species within the A. boyeri complex.