Genetic variation and differentiation of the two river sculpins,Cottus nozawae andC. amblystomopsis,deduced from allozyme and restriction enzyme-digested mtDNA fragment length polymorphism analyses

Genetic differentiation of the two sibling species,Cottus nozawae andC. amblystomopsis, from the northern part of Japan (Hokkaido Island and the Tohoku District) was investigated using allozyme variations and restriction fragment length polymorphisms of mitochondrial DNA. Although the two species are morphologically very similar, previously being thought to be a single species, they have different life-cycles;C. nozawae has a fluvial life-cycle with a small number of large-sized eggs, whereasC. amblystomopsis is an amphidromous species with a large number of small-sized eggs. Four populations ofC. amblystomopsis from Hokkaido Island and 24 populations ofC. nozawae (22 from Hokkaido Island and 2 from the Tohoku District) were sampled and examined Intrapopulational differentiation in the two species was measured by examining several indexes, including proportion of polymorphic loci (P), mean heterozygosity (H) and nucleotide diversity (π). All measurements were higher in theC. amblystomopsis populations, suggesting that intrapopulational variation inC. nozawae was less than inC. amblystomopsis and reflecting the difference in effective population sizes between them. Cluster analyses were performed using the UPGMA method, based on the data matrices of genetic distance (D) and the net nucleotide difference (δ) between populations. TheC. nozawae andC. amblystomopsis populations from Hokkaido Island composed a large cluster (Hokkaido group), while theC. nozawae populations from the Tohoku District composed a different cluster (Tohoku group). Bootstrap probabilities deduced from 1000 bootstrap replications for presence or absence of restriction sites showed that the mtDNA haplotypes detected within the Tohoku Group occurred in 99.9% of the bootstrap replicates outside the mtDNA haplotypes of the Hokkaido group, while those within the Hokkaido group occurred in 3.5–64.9% of bootstrap replicates. Consequently, the Hokkaido populations of the two species (Hokkaido group) were genetically close to each other, whileC. nozawae from the Tohoku District (Tohoku group) were distant from the Hokkaido group. These results suggest that the ancestral populations of the two species on Hokkaido Island shared the same gene pool, even after becoming geographically isolated from the ancestral population ofC. nozawae in the Tohoku District by the formation of the Tsugaru Straits.     

Patchy distribution of a fluvial sculpin,cottus nozawae, in the gakko river system at the southern margin of its native range

The distribution pattern of the freshwater sculpin,Cottus nozawae, was investigated in the Gakko River system, Tohoku District, northern Japan, that region representing the southern limit of the species distribution.C: nozawae was not found in the lower sections of the river, being instead restricted mainly to the upstream segments of some tributaries. The relationships between habitat variables and density ofC. nozawae analysed at reach scale, revealed that fish density increased with cover area, decreasing stream width and lower maximum water, temperature. Additionally, maximum water temperature was the best predictor of fish density according to stepwise multiple regression analysis. These data indicated that tributary populations ofC. nozawae were isolated from thermally unsuitable environments. Accordingly,C. nozawae was distributed patchily in headwater streams of this river system, in contrast to those on Hokkaido Island, where the species is distributed continuously throughout middle to upper stream courses.     

Genetic variations ofCottus nozawae populations from five tributaries of the Shiribetsu River of southern Hokkaido,Japan

To estimate genetic differentiation and heterogeneity in the landlocked river sculpin,Cottus nozawae, between tributary populations in the same river-system, 107 specimens were captured from 5 tributaries of the Shiribetsu river (course length 128 km), Hokkaido Island and surveyed for allozyme variations and restriction fragment length polymorphisms of mitochondrial DNA (mtDNA). Three and two alleles were seen at theIdh-2 andPgm loci, respectively, but only one locus,Idh-2, out of twenty loci examined was regarded as polymorphic, since the frequency of the most common allele did not exceed 0.95. Three different mtDNA haplotypes were detected, there being replacement of them between the tributary populations. Heterogeneities of allele and haplotype frequencies were significant between some tributary populations, suggesting that genetic differentiation has occurred between them.     

Geographic differentiation of <Emphasis Type="Italic">Cobitis shikokuensis</Emphasis> inferred from mtDNA RFLP analysis

Mitochondrial DNA divergence among populations of the Japanese spinous loach Cobitis shikokuensis, endemic to Shikoku Island, was investigated by restricted fragment length polymorphism analysis. A total of 68 restriction sites on DNA fragments from the cytochrome b to D-loop regions and from the 12S rRNA to 16S rRNA regions, amplified by PCR, were analyzed. A total of 12 haplotypes (plus 6 in outgroups) were detected in 268 specimens collected from 19 localities in seven rivers (and 41 specimens from four localities in three rivers in outgroups). Three of the seven river populations of C. shikokuensis were shown to have unique haplotypes, and four of the seven river populations were monomorphic. The nested structure of the haplotype network for populations of C. shikokuensis exhibited two large clades corresponding to (1) populations from the Shimanto River and its neighbors and (2) two genetically divergent populations in the Shigenobu and Iwamatsu Rivers. The population from the Shimanto River, the largest river inhabited by C. shikokuensis, maintains great haplotype diversity as well as the allozyme diversity previously reported. On the other hand, populations from the Hiji River, the second largest river, which exhibited the highest allozyme diversity, were monomorphic in their mtDNA. The nested clade analysis (NCA) revealed that past fragmentation between the above two clades could occur in the initial distribution process of C. shikokuensis. The large genetic divergence of two river populations from the Shigenobu and Iwamatsu Rivers was inferred to be caused by a process of long distance colonization and fragmentation. MtDNA introgression into the Hiji River population from southern river populations was suggested. Taking genetic divergence into consideration, each river population of C. shikokuensis should be conserved separately as like a distinct species, and conservation programs for the small populations showing less genetic variability should be invoked as soon as possible.     

Geographical and temporal mitochondrial DNA variability in populations of pink salmon

Pink salmon Oncorhynchus gorbuscha from odd and even year generations in rivers of Sakhalin Island, Kuril Island, Kamchatka Peninsula, and Alaska were investigated with five informative restriction endonucleases for mtDNA variation. The odd and even generations from the same rivers of South Sakhalin differed greatly. The time of divergence between the two broodlines was estimated at 0.9-1.1 Myr. The variability of mtDNA in odd year generations was higher than in even year generations and may have been due to' founder' and/or' bottleneck' effects. The differences among river populations within the Sakhalin region in 1991-1993 were not significant and this confirms the highly migratory nature of pink compared with other Pacific salmon. The mtDNA samples revealed statistically significant differences between regions. The northern populations (Kamchatka, Alaska) were less diverse in number and frequency of haplotypes than the southern populations (Sakhalin). This suggests that pink salmon originated in the Sakhalin-Kuril region and that a founder effect during the spread of this species may have restricted the mtDNA variability in other regions.     

Phylogeographic Analysis of Chum Salmon Oncorhynchus keta Walbaum in Asian Populations Based on mtDNA Variation

We used restriction length polymorphism (RFLP) analysis of PCR-amplified fragments of mtDNA to study the genetic structure of chum salmon populations sampled in 1993–2000 during a spawning run in five rivers: Narva (Southern Primorye), Naiba (Sakhalin Island), Sernovodnaya (Kunashir Island, Southern Kuril Islands), Ola (northwestern coast of the Sea of Okhotsk), and Anadyr' (Chukotka Peninsula). In total, 49 haplotypes were identified in 193 fish. Heterogeneity tests showed highly significant (P = 0) differences among all sample pairs. The estimated time of independent divergence of the populations or population groups is in good agreement with the time of Pleistocene glaciations. This result suggests that it is cyclic global changes during this time period that were crucial in determining the within-species divergence in chum salmon. The types of mtDNA genetic variability and mismatch distribution between haplotypes in the populations indicate that the southern regions of the Sea of Okhotsk and Sea of Japan served as refugia for chum salmon during glaciation periods.     

Geographical variation in the cranial and external characters of the little tube-nosed bat,Murina silvatica in the Japanese archipelago

The distribution ofMurina silvatica (Yoshiyuki, 1983) in the Japanese archipelago extends over about 2000 km from north to south. Specimens were obtained from populations in Hokkaido and Yakushima, which are at the northern and southern ends of the range, and from two intermediate locations in Honshu. Measurements of cranial and external morphology were examined for evidence of geographical variation. The results of both multivariate analysis of variance and cluster analysis showed that there was no distinct cline in skull morphology among the Hokkaido, Tohoku and Chubu populations. However, the results of multivariate analysis of variance showed that all measures were significantly greater for the Yakushima population than for the others. Similarly, in a dendrogram of cluster analysis, the Yakushima population formed a cluster that was distinct from the other populations. However, as the difference between the Yakushima population and the other populations was less than the variation found within the Hokkaido, Tohoku and Chubu populations, morphological divergence of the Yakushima population was attributable to intraspecific variation. The island of Yakushima is the most isolated of the four locations and the morphological divergence of this population may be associated with its relative geographic isolation.     

Further Analysis of Intraspecific Sequence Variation of Chloroplast DNA in Primula cuneifolia Ledeb. (Primulaceae): Implications for Biogeography of the Japanese Alpine Flora

Primula cuneifolia Ledeb. (Primulaceae), we analyzed intraspecific variation of the nucleotide sequences of non-coding regions of chloroplast DNA: the intergenic spacers between trnT (UGU) and the trnL (UAA) 5′exon, the trnL (UAA) 3′exon and trnF (GAA), and atpB and rbcL. In 20 populations of P. cuneifolia, 22 nucleotide substitutions and five insertions/deletions were inferred, and their genetic distances ranged from 0.001 to 0.008. Eight distinct haplotypes could be recognized and each haplotype was found to be geographically structured. Three major clades (the Northern, Hokkaido and Southern clades) were revealed in phylogenetic analyses of the haplotypes. The haplotypes of the Northern clade had a wider distribution area in the populations of Mt. Rausu and Rishiri Island of eastern and northern Hokkaido in Japan, northward to Unalaska Island in the Aleutians, and those of the Hokkaido clade were distributed in the populations of central Hokkaido and Mt. Iwaki of the northern Honshu in Japan; in addition, those of the Southern clade were observed only in the populations of the central Honshu. It was shown that the genetic diversifications of the Southern clade were higher than those of the Northern and Hokkaido clades. Furthermore, it was shown that the topology within the Southern clade was hierarchical, and the haplotypes of the Southern populations in the clade were derivative. From these results, we concluded that the cpDNA haplotypes of the three clades in P. cuneifolia arose and assumed the present distribution areas through several cycles of glacial advance and retreat in the Pleistocene. Received 24 June 1998/ Accepted in revised form 28 December 1998     

Genetic variability among endangered Chinese giant salamanders, Andrias davidianus

The endangered Chinese giant salamander (Andrias davidianus) is endemic to mainland China. Genetic divergence among six populations of the species was investigated by means of isozyme electrophoresis and mitochondrial DNA (mtDNA) sequences. Forty allozyme loci were resolved for all populations; the amount of genetic divergence among populations was comparable to that in other amphibians. mtDNA sequences showed a similar level of divergence. The population from Huangshan is distinct from other populations, indicating the existence of localized divergence. Both allozyme and mtDNA data failed to associate the populations into a pattern corresponding to the three Chinese river systems, which may be the consequence of human relocation. Conservation policies should emphasize the protection of localized populations and cessation of human-facilitated introductions. Future studies should focus on investigating the divergence among localized populations from isolated mountain regions, particularly using more fine-grained techniques such as microsatellite DNA.     

Genetic variation and population structure of the Japanese sika deer (Cervus nippon) in the Tohoku District based on mitochondrial D-loop sequences

The sika deer (Cervus nippon) once inhabited the entire Tohoku District, the northeastern part of the main island of Japan. Currently, they are isolated as three discontinuous populations on Mt. Goyo, the Oshika Peninsula, and Kinkazan Island. To assess the genetic diversity and relationships among the sika deer populations in the Tohoku District, we analyzed the mitochondrial DNA D-loop sequences from 177 individuals. We detected a total of five haplotypes. Three haplotypes were present in the population from Mt. Goyo at a haplotype diversity of 0.235 ± 0.061, two haplotypes in the population from the Oshika Peninsula at 0.171 ± 0.064, and only one haplotype was detected in the population from the Kinkazan Island. A significant genetic differentiation was observed among all population pairs. Collectively, our data supports the observed population bottlenecks in the past. Four of the five haplotypes were specific to one of the three populations, whereas only one haplotype was shared between the Mt. Goyo and the Oshika Peninsula populations. This common haplotype may indicate a common ancestral population in the Tohoku District. Conversely, the D-loop haplotypes were completely different among the Kinkazan Island and Oshika Peninsula populations. The lack of a shared haplotype indicates that female gene flow between the two populations is very limited and that the 0.6 km strait acts as a strong barrier.     

Genetic population structure of the fluvial eight-barbel loach <Emphasis Type="Italic">Lefua</Emphasis> sp. 1 in the three river systems in central Honshu,Japan, revealed by microsatellite DNA markers

The fluvial eight-barbel loach Lefua sp. 1 is an undescribed species distributed from the Kinki to Chugoku districts, Honshu, and also on Shikoku Island, Japan. Genetic relationships among local populations are unclear and management units remain undetermined. To aid conservation, we determined genetic population structures from microsatellite loci for 20 populations from three river systems on Honshu. The genetic diversity within populations is relatively low; the majority has experienced genetic bottlenecks. Statistical analysis revealed significant divergence among river systems suggesting that each should be recognized as a management unit. Any conservation program should consider the populations’ genetic uniqueness.     

Differential geographic patterns of mitochondrial DNA variation in two sympatric species of Japanese wood mice, Apodemus speciosus and A. argenteus

We examined the gene sequences of mitochondrial cytochrome b (cyt b) in two Japanese wood mouse species, Apodemus speciosus (n = 89) and A. argenteus (n = 46), which are distributed on the four main islands of Japan (Hokkaido, Honshu, Shikoku, and Kyushu) and on the small islands surrounding them. Apodemus speciosus, the larger of the two species, showed substantial genetic variation, with a maximum of 3% sequence divergence, and remarkable phylogenetic subdivision with two major clades. One clade represents haplotypes from a central region, including Honshu, Shikoku, Kyushu, and their adjacent islands; the other clade includes haplotypes from Hokkaido and the peripheral islands, forming four subclades: a) Hokkaido, b) Sado Island, c) Satsunan Islands, and d) the Izu Islands. Sequence divergence among the four subclades was 1.0 to 1.5%, implying that A. speciosus colonized these geographic regions 0.2 to 0.3 million years ago, assuming a substitution rate of 2.4% per million years. The population on the Izu Islands has preserved haplotypes that are distinct from those in any other region, providing good evidence for the natural colonization of the volcanic islands of the Izu Islands. The cyt b sequence variation had no relation to the karyotypic dimorphism for the eastern (2n = 48) and western (2n = 46) geographic groups, between which a strict border exists at central Honshu. On the other hand, Apodemus argenteus, the smaller of the two species, showed a similar level of sequence divergence (maximum of 3%) but no substantial geographic differentiation: populations in Hokkaido, Sado, and Yakushima shared similar haplotypes with each of the central populations, suggesting that genetic exchanges occurred among the localities in the last 0.15 million years. The apparent genetic structure of the mitochondrial DNA found in the A. speciosus population might be caused solely by long-term existence in insular regions, presumably due to ecological superiority relative to A. argenteus.     

Worldwide migration of parasitic mites as a result of bumblebee commercialization

We investigated the status of infestation by a tracheal mite, Locustacarus buchneri, in natural populations of a Japanese native bumblebee species, Bombus hypocrita, collected on Hokkaido Island and in the Aomori prefecture between 1997 and 2001. We also investigated mite infestation in commercial colonies of the European bumblebee, Bombus terrestris, imported from the Netherlands and Belgium, and the Japanese native species, B. ignitus, imported from the Netherlands, between 1997 and 2001. We detected mites in natural populations of the two B. hypocrita subspecies and in the commercial colonies. Analysis of variations in 535 bp sequences of the mitochondrial cytochrome oxidase subunit 1 (CO1) gene showed that the mite haplotypes in the native populations and in the imported colonies did not overlap in 1997–1999, but in 2000–2001 some mites possessing European CO1 haplotypes were detected in the natural populations of Japanese native bumblebees. In addition, many mites possessing Japanese haplotypes were detected in the imported commercial colonies from Europe. Considering the fact that the Japanese native bumblebees, B. hypocrita, were once exported to Europe for commercialization, these results suggest that bumblebee commercialization has caused overseas migration and cross-infestation of parasitic mites among natural and commercial colonies. However, because the Japanese and European CO1 haplotypes were closely related, there was a possibility that the European haplotypes found in the mites in the Hokkaido Island revealed native variation. To clarify the status of mite invasion, further detailed analysis of genetic variation of the mite, using other genetic markers on additional samples, need to be performed.     

Phylogenetic status of brown trout Salmo trutta populations in five rivers from the southern Caspian Sea and two inland lake basins,Iran: a morphogenetic approach

Interrelationships, origin and phylogenetic affinities of brown trout Salmo trutta populations from the southern Caspian Sea basin, Orumieh and Namak Lake basins in Iran were analysed from complete mtDNA control region sequences, 12 microsatellite loci and morphological characters. Among 129 specimens from six populations, seven haplotypes were observed. Based on mtDNA haplotype data, the Orumieh and southern Caspian populations did not differ significantly, but the Namak basin–Karaj population presented a unique haplotype closely related to the haplotypes of the other populations (0·1% Kimura two‐parameter, K2P divergence). All Iranian haplotypes clustered as a distinct group within the Danube phylogenetic grouping, with an average K2P distance of 0·41% relative to other Danubian haplotypes. The Karaj haplotype in the Namak basin was related to a haplotype (Da26) formerly identified in the Tigris basin in Turkey, to a Salmo trutta oxianus haplotype from the Aral Sea basin, and to haplotype Da1a with two mutational steps, as well as to other Iranian haplotypes with one to two mutational steps, which may indicate a centre of origin in the Caspian basin. In contrast to results of the mtDNA analysis, more pronounced differentiation was observed among the populations studied in the morphological and microsatellite DNA data, except for the two populations from the Orumieh basin, which were similar, possibly due to anthropogenic causes.     

Genetic diversity and historical population structure in the New Zealand mayfly Acanthophlebia cruentata

1. Nucleotide sequences of a 280 base pair region of the cytochrome b gene were used to assess genetic diversity and to infer population histories in the New Zealand mayfly Acanthophlebia cruentata. 2. A hierarchial examination of populations from 19 streams at different spatial scales in the central and northern North Island of New Zealand found 34 haplotypes. A common haplotype was found in all central region streams and unique haplotypes in northern streams. Several central streams had region specific haplotypes with genetically differentiated populations at the 70–100 km scale. 3. Haplotype diversity was high (0.53–0.8) at most sites, but low (0–0.22) in some central sites. amova analyses found significant genetic diversity among regions (69%) and among catchments (58%). Most population pairwise F_ST tests were significant, with non‐significant pairwise tests among sites in the central region and pairs of sites between neighbouring streams. 4. The levels of sequence divergence are interpreted as the result of Pleistocene divergence in multiple refugia, leading to the evolution of regionally unique haplotypes. The low diversity in some central region populations may result from recent colonisation following local extinctions, associated with volcanic events.     

Intraspecific mitochondrial DNA variation in the Colorado potato beetle,Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

Restriction endonuclease analyses of mitochondrial DNA (mtDNA) were used to examine genetic variability and population structure inLeptinotarsa decemlineata (Say). A group of three enzymes, EcoRI, HpaI, and PstI, was used to reveal polymorphism both within and among some of the 10 populations tested, yielding 16 haplotypes in combination. The frequencies of these 16 haplotypes differed significantly across geographic regions, indicating some partitioning of mtDNA haplotypes. Estimates of mtDNA sequence divergence (δ) between haplotypes ranged from 0.016 to 0.135%, suggesting local differentiation of mtDNA in some populations. Analysis of these data suggests that Texas was colonized by more than one mtDNA lineage, most likely originating in Mexico. We hypothesize that a larger founder size for the initial introductions or high levels of variability in the parent population at the edge of the CPB expanding range led to the initial partitioning of haplotypes observed in samples from Texas.     

Multiple origins of Locusta migratoria (Orthoptera: Acrididae) in the Japanese Archipelago and the presence of two major clades in the world: evidence from a molecular approach

Phylogenetic relationships among migratory locust (Locusta migratoria) populations in different climatic regions were analysed by sequencing four mitochondrial DNA regions, with special reference to the origin of Japanese populations. The populations are clearly separated into two clades: one consists of individuals from temperate and cold‐temperate areas of Japan and the Chinese continent, and the other comprises those from subtropical islands of Japan, Hainan Island in China, Timol Leste, Australia, Ethiopia, France, and some individuals from Tsushima Island and Honshu of Japan. The divergence time between the two clades is estimated to be 0.86–1.89 Mya. The phylogenetic analysis revealed that Japanese L. migratoria populations were composed of individuals of six different origins: (1) Hokkaido populations possibly from the Russian continent; (2) Honshu–Kyushu populations from the Chinese continent; (3) Southwest Island populations from Hainan Island or adjacent areas; (4) Ogasawara populations that might have originated from Micronesia; (5) part of the Tsushima population that originated from somewhere in the Asian tropics; and (6) a possible relict population of ancient southern haplotypes that exists in western areas of northern Honshu. The Tsugaru Straits and Tokara Straits have acted as effective geographical barriers, as in other organisms, isolating locust populations for a few thousand years. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 570–581.     

INTRASPECIFIC MOLECULAR PHYLOGENY IN THE YELLOW WARBLER (DENDROICA PETECHIA), AND IMPLICATIONS FOR AVIAN BIOGEOGRAPHY IN THE WEST INDIES

A phylogenetic analysis of mitochondrial DNA (mtDNA) restriction sites was used to examine the evolutionary history of populations of yellow warbler (Dendroica petechia) sampled from North America, Central America, South America, and the West Indies. Thirty-seven haplotypes were identified, and only one was found in more than one of these regions. Estimated sequence divergence among haplotypes ranged from 0.14 to 3.17%, and mtDNAs from North American migratory populations clearly were differentiated from those of most tropical sedentary populations. Parsimony analysis of haplotypes suggested multiple colonizations of the West Indies archipelago and of individual Caribbean islands. The inference of multiple colonizations has important implications for studies of avian ecology and evolution in this region.     

Mitochondrial DNA variability in viviparous and ovoviviparous populations of the urodele Salamandra salamandra

Mitochondrial DNA analysis of 13 populations of S. salamandra along a transect across the North of the Iberian Peninsula showed values of divergence between haplotypes ranging from d = 0.41% to 5.91%. Phenetic and cladistic analysis grouped the isofemale lineages into two main clusters with contrasting phylogeographic patterns. The first group encompasses populations located at each extreme of the Iberian Peninsula. Despite the large geographic distance separating these populations, they exhibit only a minor degree of divergence among haplotypes. In contrast, much higher diversification, in both number of distinct haplotypes, and overall genetic divergence, was observed in the second group of phylogenetically related populations. Surprisingly, this process of radiation and divergence in mtDNA haplotypes occurred in populations in close geographic proximity. All populations sampled in this group are located within a 300 km range, in the central part of our transect across the Northern edge of the Peninsula. Most populations in the central range of our transect exhibit viviparous reproduction — which is derived and highly unusual among urodeles. The genetic distances measured among Asturian (central portion of our transect), viviparous populations are higher than the distances measured between the two main taxonomic clusters. A viviparous population showing an unusual level of mtDNA heterogenetiy is reported and the potential implications of this focus of localized variability are discussed. The dynamics of isofemale lineages among the two reproductive modes was further explored in combination with the previous allozyme data. Several nuclear markers suggest that major mtDNA divergences could be explained by long-term extrinsic barriers to gene flow. Isofemale lineages indicate a narrow secondary contact zone among populations with different reproductive patterns. The existence of viviparous and ovovivparous populations sharing a common haplotype suggests that reproductive transition in S. salamandra could have arisen in absence of genetic mtDNA differentiation. We finally outline a genetic model system where the acquisition of water independence from a primitively aquatic dependent amphibian life cycle can be analyzed from a microevolutionary perspective.     

Genetic diversity of Garra rufa Heckel, 1843 (Teleostei: Cyprinidae) in Anatolia

In the present study, mitochondrial DNA polymerase chain reaction-restriction-fragment length polymorphism (PCR-RFLP) and nuclear DNA inter-simple sequence repeat (ISSR) assays were used to assess the phylogenetic and phylogeographic relationships among Garra rufa samples from Anatolia. The complete mtDNA NADH 3/4 dehydrogenase (ND-3/4) gene amplified by PCR was digested with eight restriction enzymes. These enzymes produced 20 composite haplotypes for G. rufa populations. All the mtDNA haplotypes detected were highly diverged from each other and each lineage had a unique genetic profile. The evaluation of mtDNA PCR-RFLP data coupled with geological history of Anatolia indicated a deep genetic divergence among the mtDNA haplotypes of G. rufa populations from drainages of the Mediterranean and Persian Gulf, suggesting an early isolation of Tigris-Euphrates with Orontes river and other rivers draining into the Mediterranean Sea. In general, data from both mtDNA and nDNA were congruent.