A proposal to rename the hyperthermophile Pyrococcus woesei as Pyrococcus furiosus subsp. woesei

Pyrococcus species are hyperthermophilic members of the order Thermococcales, with optimal growth temperatures approaching 100 degrees C. All species grow heterotrophically and produce H2 or, in the presence of elemental sulfur (S(o)), H2S. Pyrococcus woesei and P. furiosus were isolated from marine sediments at the same Vulcano Island beach site and share many morphological and physiological characteristics. We report here that the rDNA operons of these strains have identical sequences, including their intergenic spacer regions and part of the 23S rRNA. Both species grow rapidly and produce H2 in the presence of 0.1% maltose and 10-100 microM sodium tungstate in S(o)-free medium. However, P. woesei shows more extensive autolysis than P. furiosus in the stationary phase. Pyrococcus furiosus and P. woesei share three closely related families of insertion sequences (ISs). A Southern blot performed with IS probes showed extensive colinearity between the genomes of P. woesei and P. furiosus. Cloning and sequencing of ISs that were in different contexts in P. woesei and P. furiosus revealed that the napA gene in P. woesei is disrupted by a type III IS element, whereas in P. furiosus, this gene is intact. A type I IS element, closely linked to the napA gene, was observed in the same context in both P. furiosus and P. woesei genomes. Our results suggest that the IS elements are implicated in genomic rearrangements and reshuffling in these closely related strains. We propose to rename P. woesei a subspecies of P. furiosus based on their identical rDNA operon sequences, many common IS elements that are shared genomic markers, and the observation that all P. woesei nucleotide sequences deposited in GenBank to date are > 99% identical to P. furiosus sequences.     

Recurrent high-biomass blooms of <Emphasis Type="Italic">Alexandrium taylorii</Emphasis> (Dinophyceae), a HAB species expanding in the Mediterranean

Summer outbreaks of the dinoflagellate Alexandrium taylorii Balech are recurrent events in nearshore waters of Sicily (Italy)—a central region in the Mediterranean Sea—producing dense yellowish–green patches. Beyond the local phenomenon, the problem covers a broader geographic scale, involving also other European localities, mostly in Spain. Biological, environmental, and molecular data are reported here from a semi-closed bay of Sicily (Vulcano Island, Tyrrhenian Sea, 2000–2003), showing in summer the recurrence of high-biomass blooms and events of water discolouration. Without underestimating the setbacks to the tourism industry, the ecological impact of A. taylorii blooms may be important considering the high levels of biomass produced (West Bay, Vulcano: up to a magnitude order of 10⁷ cells l⁻¹, 50–180 μg-Chla l⁻¹, June 2002 and 2003) and coincident conditions of oxygen supersaturation of the waters (130–170%). Trophic trends in the Tyrrhenian site indicate high amounts of nutrients linked to the increased anthropogenic activity in summer, although recently there has been an apparent shift of the marked eutrophic conditions towards a slighter eutrophy. Genetic data on isolates of A. taylorii from the Mediterranean Sea are also discussed. Molecular analyses implied the sequencing of target rDNA regions (5.8S rDNA and ITS regions) of several isolates from different Mediterranean localities, as well as the application of species-specific PCR assays for rapid species identification in preserved field samples. The confirmation of the specific identity provided new insights into the biogeography of this species and further evidence of the occurrence of A. taylorii in a number of Mediterranean localities, both in the western side (the Catalan coast of Spain) and the eastern area (Greece). Analyses of the molecular diversity of geographically distinct isolates of A. taylorii from Italy, Spain, and Greece based on the 5.8S rDNA-ITS region sequences showed a high level of similarity, indicating the existence of an unique Mediterranean population.     

The population structure of ten Newfoundland outports

Island populations are most informative in the study of the genetic structure of human aggregates. These populations are often of small size, thus violating the Hardy-Weinberg assumption of infinite size. Some geographically isolated island populations are further subdivided by religion, ethnicity, and socioeconomic factors, reducing their effective sizes and facilitating genetic changes due to stochastic processes. Because of extreme geographic and social isolation, fishing communities or outports of Newfoundland have been investigated for genetic microdifferentiation through the founder effect and genetic drift (Crawford et al. 1995). The purpose of this paper is to examine the population structure of 10 Newfoundland outports using the allelic frequencies derived from 12 red cell antigens. To achieve this goal, first we calculated gene frequencies using maximum-likelihood estimation procedures. Second, we used R-matrix methods to explore population differentiation. Third, we regressed mean per-locus heterozygosity on genetic distance from the gene frequency centroid to identify the most isolated populations. On the basis of this information, the three outports of Seal Cove, Island Harbor, and Tilting were found to be genetically differentiated from the other small populations. Moreover, religious and geographic subdivisions appear to explain the observed genetic variation.     

Acquisition of insertion sequences and the GBSi1 intron by Streptococcus agalactiae isolates correlates with the evolution of the species

The prevalence of insertion sequences IS1548, IS861, IS1381, and ISSa4 and of the group II intron GBSi1 within Streptococcus agalactiae human isolates strongly correlates with the genetic lineages obtained by multilocus sequence typing. Our results yielded an evolutionary scheme for the acquisition of these genetic elements linked to the ecosystems from which the isolates were obtained.     

Evolutionary dynamics of insertion sequences in Helicobacter pylori

Prokaryotic insertion sequence (IS) elements behave like parasites in terms of their ability to invade and proliferate in microbial gene pools and like symbionts when they coevolve with their bacterial hosts. Here we investigated the evolutionary history of IS605 and IS607 of Helicobacter pylori, a genetically diverse gastric pathogen. These elements contain unrelated transposase genes (orfA) and also a homolog of the Salmonella virulence gene gipA (orfB). A total of 488 East Asian, Indian, Peruvian, and Spanish isolates were screened, and 18 and 14% of them harbored IS605 and IS607, respectively. IS605 nucleotide sequence analysis (n = 42) revealed geographic subdivisions similar to those of H. pylori; the geographic subdivision was blurred, however, due in part to homologous recombination, as indicated by split decomposition and homoplasy tests (homoplasy ratio, 0.56). In contrast, the IS607 populations (n = 44) showed strong geographic subdivisions with less homologous recombination (homoplasy ratio, 0.2). Diversifying selection (ratio of nonsynonymous change to synonymous change, >1) was evident in approximately 15% of the IS605 orfA codons analyzed but not in the IS607 orfA codons. Diversifying selection was also evident in approximately 2% of the IS605 orfB and approximately 10% of the IS607 orfB codons analyzed. We suggest that the evolution of these elements reflects selection for optimal transposition activity in the case of IS605 orfA and for interactions between the OrfB proteins and other cellular constituents that potentially contribute to bacterial fitness. Taken together, similarities in IS elements and H. pylori population genetic structures and evidence of adaptive evolution in IS elements suggest that there is coevolution between these elements and their bacterial hosts.     

Phylogeography of Asparagopsis taxiformis revisited: Combined mtDNA data provide novel insights into population structure in Japan

Six intraspecific lineages (Lineages 1–6) of Asparagopsis taxiformis have been previously established based on mitochondrial cox2‐cox3 intergenic spacer and a partial cox1 sequences. ‘Lineage 2’ (L2) was suggested to be a recent introduction to the Mediterranean Sea, but its source population has not yet been identified. In order to clarify the nature of northwestern Pacific populations, we performed extensive sampling in Japan (60 individuals from 16 locations) and molecular phylogenetic analyses based on mitochondrial sequences. Sixteen additional individuals, collected from eight locations in the Indo‐Pacific, Caribbean, and Mediterranean regions, were also analyzed. Combined sequence analyses revealed that the Japanese populations only consisted of L2. Out of 19 combined haplotypes identified within L2, two are shared between Japan and the Mediterranean Sea and the Hawaiian Islands, and 12 were identified as endemic to Japan. Genetic analyses of population differentiation suggested that Japanese populations are genetically isolated from the Mediterranean and the Hawaiian populations. A genetic disjunction appears to separate two subpopulations within Japan: one between Toi and Kagoshima and the other between Ojikajima Island and Kagoshima in the Kyushu area.     

Identification of nodule-dominant Rhizobium meliloti strains carrying pRmeGR4b-type plasmid within indigenous soil populations by PCR using primers derived from specific DNA sequences

Abstract: Rhizobium meliloti strain GR4 is a highly infective and competitive bacteria which was isolated in 1975 from a field site in Granada (Spain) and which has a high potential as an inoculant. R. meliloti isolates from alfalfa plants grown in this field site were characterized using polymerase chain reaction. Characterization was based on primers derived from insertion sequence elements (IS Rm3 and IS Rm4 ), plasmid origin of replication (pRmeGR4a repC locus) and plasmid pRmeGR4b specific DNA sequences. Soil isolates harbouring plasmid type pRmeGR4b represented the major infective population in this field site. A direct correlation between the presence of pRmeGR4b-like plasmid and the competitiveness of the strains was found. In addition, four different R. meliloti field populations isolated from Spanish soils were analyzed for the presence of pRmeGR4b related plasmids. Our results indicate that this plasmid type is widespread among R. meliloti field populations and that its frequency within the infective isolates depends on the host plant.     

Living in Vulcan's forge: Algal adaptation to stressful geothermal ponds on Vulcano Island (southern Italy) as a result of pre-selective mutations

Four species of eukaryotic algae proliferate in the sulfureous, acidic (pH 3.1) water of the largest geothermal pond on Vulcano Island (southern Italy). Consequently, this pond constitutes a natural laboratory for analysis of adaptation by phytoplankters to extremely stressful conditions. To distinguish between the pre-selective or post-selective origin of adaptation processes allowing the existence of phytoplankters in the pond, a Luria-Delbrück fluctuation test was carried out with the chlorophycean Dictyosphaerium chlorelloides and the cyanobacterium Microcystis aeruginosa , both isolated from non-extreme waters; natural water from the Vulcano Island pond was used as selective factor. Preselective, resistant D. chlorelloides cells appeared with a frequency of 4.7 × 10⁻⁷ per cell per generation. We propose that the micro-algae inhabiting this stressful pond could be the descendents of chance mutants that arrived in the past or are even arriving at the present. The genetic adaptation of D. chlorelloides to Vulcano waters could help to explain the survival of photosynthesizers in very stressful geothermal waters during the Neoproterozoic 'snowball Earth', a period when primary production collapsed in the biosphere. On the other hand, adaptation to these conditions was not observed in M. aeruginosa , suggesting that cyanobacteria may not be able to develop any kind of adaptation to Vulcano pond water.     

Temporal changes of genetic population structure and diversity in the endangered Blakiston's fish owl (Bubo blakistoni) on Hokkaido Island, Japan, revealed by microsatellite analysis

The Blakiston's fish owl (Bubo blakistoni) population on Hokkaido Island, Japan, decreased to less than one hundred individuals over the last century due to habitat disruption by human activity. Although the ongoing conservation management has slightly restored the population, it remains endangered. In order to assess the genetic variation and population structure of the Blakiston's fish owl in Hokkaido, we genotyped eight microsatellite loci on 120 individuals sampled over the past three decades. The genotype data set showed low levels of genetic variation and gene flow among the geographically isolated five subpopulations. Comparative analysis of past and current populations indicated that some alleles shared by past individuals had been lost, and that genetic variation had declined over the last three decades. The result suggests that the genetic decline may have resulted from inbreeding and/or genetic drift due to bottlenecks in the Hokkaido population. The present study provides invaluable genetic information for the conservation and management of the endangered Blakiston's fish owl in Hokkaido.     

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.     

Contemporary genetic structure and postglacial demographic history of the black scorpionfish,Scorpaena porcus,in the Mediterranean and the Black Seas

Understanding the distribution of genetic diversity in the light of past demographic events linked with climatic shifts will help to forecast evolutionary trajectories of ecosystems within the current context of climate change. In this study, mitochondrial sequences and microsatellite loci were analysed using traditional population genetic approaches together with Bayesian dating and the more recent approximate Bayesian computation scenario testing. The genetic structure and demographic history of a commercial fish, the black scorpionfish, Scorpaena porcus, was investigated throughout the Mediterranean and Black Seas. The results suggest that the species recently underwent population expansions, in both seas, likely concomitant with the warming period following the Last Glacial Maximum, 20 000 years ago. A weak contemporaneous genetic differentiation was identified between the Black Sea and the Mediterranean Sea. However, the genetic diversity was similar for populations of the two seas, suggesting a high number of colonizers entered the Black Sea during the interglacial period and/or the presence of a refugial population in the Black Sea during the glacial period. Finally, within seas, an east/west genetic differentiation in the Adriatic seems to prevail, whereas the Black Sea does not show any structured spatial genetic pattern of its population. Overall, these results suggest that the Black Sea is not that isolated from the Mediterranean, and both seas revealed similar evolutionary patterns related to climate change and changes in sea level.     

Phylogenetic diversity of Synechococcus strains isolated from the East China Sea and the East Sea

Phylogenetic relationships among 33 Synechococcus strains isolated from the East China Sea (ECS) and the East Sea (ES) were studied based on 16S rRNA gene sequences and 16S–23S rRNA gene internal transcribed spacer (ITS) sequences. Pigment patterns of the culture strains were also examined. Based on 16S rRNA gene and ITS sequence phylogenies, the Synechococcus isolates were clustered into 10 clades, among which eight were previously identified and two were novel. Half of the culture strains belonged to clade V or VI. All strains that clustered into novel clades exhibited both phycoerythrobilin and phycourobilin. Interestingly, the pigment compositions of isolates belonging to clades V and VI differed from those reported for other oceanic regions. None of the isolates in clade V showed phycourobilin, whereas strains in clade VI exhibited both phycourobilin and phycoerythrobilin, which is in contrast to previous studies. The presence of novel lineages and the different pigment patterns in the ECS and the ES suggests the possibility that some Synechococcus lineages are distributed only in geographically restricted areas and have evolved in these regions. Therefore, further elucidation of the physiological, ecological, and genetic characteristics of the diverse Synechococcus strains is required to understand their spatial and geographical distribution.     

Evidence of recent lateral gene transfer among hyperthermophilic archaea

A total of 153 nucleotide differences were found over a contiguous 16 kb region between two hyperthermophilic Archaea, Pyrococcus furiosus and Thermococcus litoralis. The 16 kb region in P. furiosus is flanked by insertion sequence (IS) elements with inverted and direct repeats. Both IS elements contain a single open reading frame (ORF) encoding a putative protein of 233 amino acids identified as a transposase. This 16 kb region has the features of a typical bacterial composite transposon and represents a possible mechanism for lateral gene transfer between Archaea or possibly between Archaea and Bacteria. A total of 23 homologous IS elements was found in the genome sequence of P. furiosus, whereas no full-length IS elements were identified in the genomes of Pyrococcus abyssi and Pyrococcus horikoshii. Only one IS element was found in T. litoralis. In P. furiosus and T. litoralis, the 16 kb region contains an ABC transport system for maltose and trehalose that was characterized biochemically for T. litoralis. Regulation of expression studies showed that the malE gene, located on the transposon, and the encoded trehalose/maltose-binding protein (TMBP) are induced in the presence of maltose and trehalose in both P. furiosus and T. litoralis. The implications of transposition as a mechanism for lateral gene transfer among Archaea are discussed.     

Analysis of insertion sequences in thermophilic cyanobacteria: exploring the mechanisms of establishing, maintaining, and withstanding high insertion sequence abundance

Insertion sequences (ISs) are simple mobile genetic elements capable of relocating within a genome. Through this transposition activity, they are known to create mutations which are mostly deleterious to the cell, although occasionally they are beneficial. Two closely related isolates of thermophilic Synechococcus species from hot spring microbial mats are known to harbor a large number of diverse ISs. To explore the mechanism of IS acquisition within natural populations and survival in the face of high IS abundance, we examined IS content and location in natural populations of Synechococcus by comparing metagenomic data to the genomes of fully sequenced cultured isolates. The observed IS distribution in the metagenome was equivalent to the distribution in the isolates, indicating that the cultured isolates are appropriate models for the environmental population. High sequence conservation between IS families shared between the two isolates suggests that ISs are able to move between individuals within populations and between species via lateral gene transfer, consistent with models for IS family accumulation. Most IS families show evidence of recent activity, and interruption of critical genes in some individuals was observed, demonstrating that transposition is an ongoing mutational force in the populations.     

Genotypic characteristics of Mycobacterium avium subsp. hominissuis strains

Mycobacterium avium are typical environmental, non-tuberculosis microorganisms that occasionally cause mycobacteriosis, an infectious disease in wild and domestic animals, birds, and humans. Here, we report the results of the first study on the genetic diversity of the Russian population of M. avium. A total of 85 M. avium subsp. hominissuis (MAH) clinical strains were isolated from patients (including 30 HIV-positive individuals) with mycobacteriosis in St. Petersburg, 2008–2011. The identification of the microorganisms was carried out using biochemical tests and the PCR detection of the mobile elements IS901 and IS900, as well as of the polymorphism of restriction fragments of the hsp65 gene. The genetic diversity of the isolates was evaluated by VNTR typing based on eight variable-number tandem repeats (VNTRs) (292, X3, 25, 47, 3, 7, 10, and 32 [Thibault et al., 2007]). The MAH population studied was characterized by 15 VNTR types, including nine unique patterns and six clusters of isolates with identical eight-digit profiles. The largest clusters (22221128 and 24221128) included 45 (59.2%) and 15 (19.7%) isolates, respectively; the others contained 2–7 strains. The strains of the cluster 2533112’8 possessed a truncated TR10 locus (allele 2′). Taking into account the absence of the epidemiological links between the patients and the fact that the infection was presumably delivered from the environment, the high rate of clustering of MAH isolates can be explained by the low discriminatory power of the eight-locus VNTR-typing scheme (HGDI 0–0.61).     

The structure of allozyme and leaf shape variation in isolated, range-margin populations of the shrub Hippocrepis emerus (Leguminosae)

We investigated allozyme variation and phenotypic variation in leaf shape in 15 populations of the shrub Hippocrepis emerus (Leguminosae) from the three isolated, regional populations in Scandinavia and analysed patterns of differentiation and the hierarchical structuring of diversity on different geographic scales There are pronounced geographic differences between the Norwegian and Swedish isolates of the species and most of the polymorphic allozyme loci show reciprocal fixation in the two isolates The Scandinavian populations of H emerus are not only disjunct on a macrogeographic scale but also show considerable disjunction within regions Within the Oland regional population, a central group of populations shows low levels of Inter-population differentiation the Gotland group of populations is related to this core group of Oland populations The geographically marginal populations on Oland are spatially isolated and show a higher degree of divergence between populations than does the central group of Oland populations We interpret genetic divergence between the marginal populations in terms of genetic drift - as a result of historical fluctuations in habitat availability and population size     

The novel insertion sequences IS1417, IS1418, and IS1419 from Burkholderia glumae and their strain distribution

Three insertion sequences, IS1417, IS1418, and IS1419, were isolated from Burkholderia glumae (formerly Pseudomonas glumae), a gram-negative rice pathogenic bacterium, on the basis of their abilities to activate the expression of the neo gene of the entrap vector pSHI1063. The 1335-bp IS1417 element with 17-bp imperfect terminal inverted repeats was found to be flanked by 5-bp direct repeats of the vector sequence. IS1418 is 865 bp in length and carries 15-bp inverted repeats with a target duplication of 3 bp. The 1215-bp IS1419 sequence is bounded by the 36-bp terminal inverted repeats of the element and 7-bp direct repeats of the vector sequence. IS1417 and IS1418 belong to the IS2 subgroup of the IS3 family and the IS427 subgroup of the IS5 family, respectively, whereas IS1419 does not appear to be a member of any known IS family. Southern blot analysis of DNAs from B. glumae field isolates indicated that those IS elements are widely distributed, but the host range of the three IS elements appears to be limited to B. glumae and some other related species such as B. plantarii. The polymorphisms exhibited in B. glumae isolates suggest that those elements are useful for molecular epidemiological studies of B. glumae infections.     

Novel insertion sequence IS1380 from Acetobacter pasteurianus is involved in loss of ethanol-oxidizing ability.

Acetobacter pasteurianus NCI1380, a thermophilic strain isolated from the surface culture of acetic acid fermentation, showed genetic instability to produce at high frequency spontaneous mutants which were deficient in ethanol oxidation because of the loss of alcohol dehydrogenase activity. Southern hybridization experiments with the cloned alcohol dehydrogenase-cytochrome c gene cluster as the probe showed insertion of an unknown DNA fragment into a specific position in the cytochrome c gene in most of the mutant strains. Cloning and sequencing analyses revealed that the inserted sequence was 1,665 bp in length and had a terminal inverted repeat of 15 bp. In addition, this inserted sequence was found to generate a 4-bp duplication at the inserted site upon transposition. The target site specificity was not very strict, but a TCGA sequence appeared to be preferentially used. The inserted sequence contains two long open reading frames of 461 and 222 amino acids which are overlapped and encoded by different strands. Although these open reading frames showed no homology to any protein registered in the DNA data bases, the longer open reading frame contained many basic amino acids (87 of 461), as was observed with transposases of so-called insertion sequence (IS) elements. All of these characteristics are typical of IS elements, and the sequence was named IS1380. The copy number of IS1380 in a cell of A. pasteurianus NCI1380 was estimated to be about 100. Several strains of acetic acid bacteria also contained IS1380 at high copy numbers. These results suggest that IS1380 is associated with the genetic loss of ethanol-oxidizing ability as well as the genetic instability of acetic acid bacteria in general.     

Acidophiles of saline water at thermal vents of Vulcano,Italy

DNA was extracted from samples taken from close to acidic hydrothermal vents on shore of the Aeolian Island of Vulcano (Italy). RNA gene sequences were amplified by PCR, cloned, and sequenced. A sequence with an origin in samples at 35 degrees and 45 degrees C corresponded to that of a novel Acidithiobacillus species that was isolated from water close to the vents. Novel, iron-oxidizing mesophilic acidophiles were isolated through enrichment cultures with ferrous iron but were not represented in the clone banks of environmental rDNA. These acidophiles were related to Thiobacillus prosperus, which was isolated previously from Vulcano. The archaeal sequences that comprised a clone bank representing a high-temperature sample (75 degrees C) corresponded to those of Acidianus brierleyi and of thermophiles previously isolated from Vulcano, Thermoplasma volcanium and Acidianus infernus.     

Genetic and morphological analyses of tub gurnard Chelidonichthys lucerna populations in Turkish marine waters

Genetic and morphological structure of tub gurnard Chelidonichthys lucerna populations in Turkish marine waters were investigated with mtDNA sequencing of 16S rRNA and morphological characters. C. lucerna samples were collected from the Black Sea, Marmara, Aegean and northeastern Mediterranean coasts of Turkey. The lowest genetic diversity was found in the northeastern Mediterranean (Iskenderun Bay) population, while the highest was in the Marmara population with overall average value of genetic diversity within populations. A total of 14 haplotypes was found, and the highest haplotype diversity was in the Black Sea whereas the lowest was in the northeastern Mediterranean population (Iskenderun Bay). The Black Sea and Iskenderun Bay populations showed the least genetic divergence (0.001081), while the highest was between the Marmara Sea and northeastern Mediterranean (Antalya Bay) populations (0.002067). Pairwise comparisons of genetic distance revealed statistically significant differences (P < 0.05) between the Marmara and both the Aegean and northeastern Mediterranean (Antalya Bay) samples. Neighbour joining tree analyses clustered the northeastern Mediterranean populations (Antalya Bay and Iskenderun Bay) as genetically more interrelated populations, whereas the Aegean Sea population was clustered as most isolated one. Discriminant function analysis of morphological characters showed that only the Black Sea population is differentiated from the other populations.