Mapping Bombina mitochondrial genomes: the conundrum of Carpathian Bombina variegata (Anura: Discoglossidae)

Complete mitochondrial genomes (mtDNA) of five individuals representing two haplotypes of Bombina bombina and three of Bombina variegata were compared using restriction site maps. Phylogenetic analyses reveal three ancient mitochondrial lineages: (1) two very similar haplotypes A and B of B. bombina ; (2) almost identical haplotypes D and E of B. variegata ; and (3) haplotype C of B. variegata . Haplotype C is as different from haplotypes D/E as from A/B. These data are strikingly discordant with relationships based on morphology and allozymes. Haplotypes C and D/E represent a pre-Pleistocene mitochondrial divergence within B. variegata , nearly coincident with speciation between B. variegata and B. bombina . Geographical partitioning of the two divergent B. variegata mitochondrial lineages indicates repeated localization of the lineages in separate glacial refugia during the Pleistocene. That nuclear genes do not show a similar divergence, but rather indicate relatively free genetic exchange between populations with divergent mtDNAs, suggests that males dispersed much more widely than females during expansions from glacial refugia. Comparison of Bombina mtDNA maps with a restriction site map of Xenopus laevis mtDNA revealed 16 homologous sites; 12 of these may be nearly invariant across primitive anuran mtDNAs. Two distinct regions of heteroplasmy, representing two regions with variable numbers of sequence repeats [length variable (LV) regions], were characterized. Comparison with the Xenopus map places LV1, present in all five haplotypes, near the 5'-end of the control region, and LV2 present only in B. variegata , near the 3'-end. Although phylogenetic analyses did not group the two major B. variegata lineages together, presence of LV2 in both lineages supports placement of both within B. variegata .     

Phylogeography of the fire-bellied toads Bombina: independent Pleistocene histories inferred from mitochondrial genomes

The fire-bellied toads Bombina bombina and Bombina variegata, interbreed in a long, narrow zone maintained by a balance between selection and dispersal. Hybridization takes place between local, genetically differentiated groups. To quantify divergence between these groups and reconstruct their history and demography, we analysed nucleotide variation at the mitochondrial cytochrome b gene (1096 bp) in 364 individuals from 156 sites representing the entire range of both species. Three distinct clades with high sequence divergence (K2P = 8-11%) were distinguished. One clade grouped B. bombina haplotypes; the two other clades grouped B. variegata haplotypes. One B. variegata clade included only Carpathian individuals; the other represented B. variegata from the southwestern parts of its distribution: Southern and Western Europe (Balkano-Western lineage), Apennines, and the Rhodope Mountains. Differentiation between the Carpathian and Balkano-Western lineages, K2P approximately 8%, approached interspecific divergence. Deep divergence among European Bombina lineages suggests their preglacial origin, and implies long and largely independent evolutionary histories of the species. Multiple glacial refugia were identified in the lowlands adjoining the Black Sea, in the Carpathians, in the Balkans, and in the Apennines. The results of the nested clade and demographic analyses suggest drastic reductions of population sizes during the last glacial period, and significant demographic growth related to postglacial colonization. Inferred history, supported by fossil evidence, demonstrates that Bombina ranges underwent repeated contractions and expansions. Geographical concordance between morphology, allozymes, and mtDNA shows that previous episodes of interspecific hybridization have left no detectable mtDNA introgression. Either the admixed populations went extinct, or selection against hybrids hindered mtDNA gene flow in ancient hybrid zones.     

Comparative Analysis of Mitochondrial Genomes in <Emphasis Type="Italic">Bombina</Emphasis> (Anura; Bombinatoridae)

The complete mitochondrial genomes of two basal anurans, Bombina bombina and B. variegata (Anura; Bombinatoridae), were sequenced. The gene order of their mitochondrial DNA (mtDNA) is identical to that of canonical vertebrate mtDNA. In contrast, we show that there are structural differences in regulatory regions and protein coding genes between the mtDNA of these two closely related species. Corrected sequence divergence between the mtDNA of B. bombina and B. variegata amounts to 8.7% (2.3% divergence in amino acids). Comparisons with two East Asian congeners show that the control region contains two repeat regions, LV1 and LV2, present in all species except for B. bombina, in which LV2 has been secondarily lost. The rRNAs and tRNAs are characterized by low nucleotide divergence. The protein coding genes are considerably more disparate, although functional constraint is high but variable among genes, as evidenced by dN/dS ratios. A mtDNA phylogeny established the distribution of autapomorphic nonsynonomous substitutions in the mitogenomes of B. bombina and B. variegata. Nine of 98 nonsynonomous substitutions led to radical amino acid replacements that may alter mitochondrial protein function. Most radical substitutions were found in ND2, ND4, or ND5, encoding mitochondrial subunits of complex I of the electron transport system. The extensive divergence between the mitogenomes of B. bombina and B. variegata is discussed in terms of its possible role in impeding gene flow in natural hybrid zones between these two species.     

Evolution of Bombina bombina and Bombina variegata (Anura: Discoglossidae) in the Carpathian Basin: a history of repeated mt-DNA introgression across species

The structure and geographic location of hybrid zones change through time. Current patterns result from present and historical population-environment interactions that act on each of the hybridizing taxa. This is particularly evident for species involved in complex hybrid zones, such as that formed by the toad species Bombina bombina and Bombina variegata (Anura: Discoglossidae), which interact along extensive areas in Central Europe. We used data on external morphology and partial sequences of the cytochrome oxidase I (cox1) and nicotinamide adenine dinucleotid dehydrogenase subunit 4 (nad4) mitochondrial DNA (mt-DNA) genes to analyze the current patterns of genetic structure shown by both species of Bombina along their contact zone in Hungary. Phylogenetic, phylogeographic, and historical demography analyses were applied to 1.5kb mt-DNA obtained from 119 individuals representing 24 populations from Hungary and additional specimens from Slovakia, Albania, and Bosnia-Herzegovina. We use these data to infer the evolutionary history of the isolated populations of B. variegata in Hungary and to discriminate between competing biogeographic scenarios accounting for the historical interactions between species in this region. Results from the inferred phylogenetic branching pattern and sequence divergence among species and populations support the following: (i) recent population expansion has occurred in Hungarian populations of B. bombina, which are genetically very homogeneous; (ii) the Hungarian populations of B. variegata correspond to two distinct mitochondrial lineages (Carpathian and Alpine, respectively); average maximum-likelihood-corrected sequence divergence between these lineages is 8.96% for cox1 and 10.85% for nad4; (iii) mt-DNA divergence among the three isolated western populations of B. variegata from Transdanubia is low, with four closely related haplotypes, which suggests that the isolation between these populations is the result of a recent process, possibly mediated by the invasion of B. bombina; and (iv) we have detected discordances between morphology and mt-DNA data in the Transdanubia region (Bakony Mountains, Mecsek Mountains, Orség area), suggesting mt-DNA introgression across species in this regions. These results are discussed with reference to previous biogeographic hypotheses.     

Limited mitochondrial DNA introgression in a Bombina hybrid zone

We analysed the variation at mitochondrial DNA (mtDNA) and six allozyme loci in two transects across a hybrid zone between the fire-bellied toads Bombina bombina and Bombina variegata in southern Poland. The mtDNA cline was narrower than allozyme clines in one transect (Przemyśl, β = 0.435, w  = 3.42 km) and shifted to the B .  bombina side in both (near Kraków by 1.32 km, α = 0.437, and near Przemyśl by 0.97 km, α = 0.319). Cytonuclear associations were weak. Narrowed mtDNA clines could be a by-product of female demography and lowered effective population size of mtDNA. Alternatively, restricted mtDNA introgression is a likely result of negative epistasis in recombinants augmented by environment-dependent selection because divergence of Bombina mitochondrial and nuclear genomes is large. The shift of the mtDNA cline, contrary to expectations from initial mating preferences and fecundity differences between the species, suggests that on the B .  bombina side of the zone hybrid females with B .  variegata mtDNA have a higher chance of leaving progeny.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 295–306.     

Quantitative immunological studies of the albumins of several species of fire bellied toads, genus Bombina.

1. Rabbit antisera against purified serum albumin of Bombina bombina were used to study relationships between B. bombina, B. variegata, and B. orientalis. 2. Quantitative micro-complement fixation tests indicated the albumins of B. bombina from central Poland and Bulgaria were indistinguishable. The albumins from several populations of B. variegata differed very slightly from that of B. bombina. The albumin of B. orientalis was quite distinct from that of B. bombina. 3. Using albumin as a molecular clock, we estimated B. bombina and B. variegata diverged within the last million years, whereas the B. orientalis lineage diverged roughly 10-12 mil yr ago.     

Reproductive biology of the fire-bellied toads, Bombina bombina and B. variegata (Anura: Discoglossidae): egg size, clutch size and larval period length differences

Clutch size, egg size, larval size and duration of larval period were estimated for samples of Bombina bombina and B. variegata from several populations. Bombina variegata lays fewer eggs per batch (mean 17.4) than B. bombina (mean 32.5). The volume of B. variegata eggs is 2.4 times greater and the dry weight 2.3 times greater than that of B. bombina. As a consequence, the volume of B. variegata hatchlings and their dry weight are respectively 1.9 times and 2.5 times greater than that of B. bombina. Analysis of variance showed that variation of egg size in B. variegata is significantly greater than in B. bombina. When reared under laboratory conditions freshly metamorphosed individuals of both species showed no difference in body size or weight. Although growth rates are equal, tadpoles of B. variegata reached metamorphosis 14 days earlier than B. bombina tadpoles because of differences in size of hatchlings. Differences in reproductive parameters found for the two European Bombina species can be interpreted as adaptations to contrasting environments in which they breed.     

Differentiation among Austrian populations of Norway spruce [Picea abies (L.) Karst.] assayed by mitochondrial DNA markers

Thirty-seven populations of Norway spruce [Picea abies (L.) Karst.] across the Austrian Alps and Bohemian Massif were sampled to elucidate the geographical pattern of genetic differentiation. Three polymorphic mitochondrial DNA (mtDNA) loci were surveyed. Two or three alleles were detected at each locus, resulting in seven multilocus mtDNA haplotypes (A–F). Western populations proved to be monomorphic, whereas eastern and central Austrian populations were slightly to highly polymorphic, respectively. As revealed by spatial analysis of molecular variance and Monmonier’s analysis, the two main haplotypes A and B are not randomly distributed. Haplotype A was restricted to central and eastern Austria, whereas haplotype B occurred in all Austrian populations but was the only haplotype identified in western populations. This pattern may be explained by different glacial refugia located in the Dinaric Alps and the Carpathian mountains.     

A rangewide phylogeography of Hermann's tortoise, Testudo hermanni (Reptilia: Testudines: Testudinidae): implications for taxonomy

Hermann's tortoise (Testudo hermanni), the best-known western Palaearctic tortoise species, has a rare natural distribution pattern comprising the Mediterranean areas of the Iberian, Apennine, and Balkan Peninsulas, as well as Sicily, Corsica and Sardinia. The western part of this range is traditionally considered habitat for T. h. hermanni, while T. h. boettgeri occurs in the Balkans. Taxonomy of this tortoise has been challenged in recent years, with the two subspecies being considered full species and the central Dalmatian populations of T. h. boettgeri being considered a third species, T. hercegovinensis. Using an mtDNA fragment approximately 1150 bp long (cytochrome b gene and adjacent portion of tRNA-Thr gene), we investigated mtDNA diversity with regard to contrasting concepts of two subspecies or three species. Seven closely related haplotypes were identified from the western Mediterranean and 15 different, in part much-differentiated, haplotypes from the Balkans. Western Mediterranean haplotypes differ from Balkan haplotypes in 16–42 mutation steps. One to seven mutation steps occur within western Mediterranean populations. Balkan haplotypes, differing in 1−37 nucleotides, group in parsimony network analysis into three major assemblages that display, in part, a similar degree of differentiation to that of western Mediterranean haplotypes relative to Balkan haplotypes. Rates of sequence evolution are different in both regions, and low divergence, palaeogeography and the fossil record suggest a slower molecular clock in the western Mediterranean. While monophyly in western Mediterranean haplotypes is well-supported, conflicting evidence is obtained for Balkan haplotypes; maximum parsimony supports monophyly of Balkan haplotypes, but other phylogenetic analyses (Bayesian, ML, ME) indicate Balkan haplotypes could be paraphyletic with respect to the western Mediterranean clade. These results imply a process of differentiation not yet complete despite allopatry in the western Mediterranean and the Balkans, and suggest all populations of T. hermanni are conspecific. In the western Mediterranean no clear geographical pattern in haplotype distribution is found. Distribution of Balkan haplotypes is more structured. One group of similar haplotypes occurs in the eastern Balkans (Bulgaria, Republic of Macedonia, Romania and the Greek regions Evvia, Macedonia, Peloponnese, Thessaly and Thrace). Two distinct haplotypes, differing in eight to nine mutation steps from the most common haplotype of the first group, are confined to the western slope of the Taygetos Mts. in the Peloponnese. Yet another group, connected over between four and 23 mutation steps with haplotypes of the eastern Balkan group, occurs along the western slope of the Dinarid and Pindos Mts. (Istria, Dalmatia, western Greece). Taygetos haplotypes are nested within other haplotypes in all phylogenetic analyses and support for monophyly of the other Balkan groups is at best weak. We conclude that using the traditional two subspecies model should be continued for T. hermanni. Phylogeographies of T. hermanni and Emys orbicularis, another codistributed chelonian, are markedly different, but share a few similarities. Both were forced to retreat to southern refuges during Pleistocene glaciations. With the advent of Holocene warming, E. orbicularis underwent rapid range expansion and temperate regions of Europe and adjacent Asia were recolonized from refuges in the Balkans and the northern Black Sea Region. By contrast, T. hermanni remained more or less confined to refuges and nearby regions, resulting in a much smaller range, and allopatric and parapatric distribution of haplotype groups and clades. MtDNA lineages are more diverse in E. orbicularis than they are in T. hermanni on southern European peninsulas, indicating several distinct glacial refuges in close proximity and extensive intergradation during Holocene range expansion for E. orbicularis. In T. hermanni it is likely that only on the Balkan Peninsula was more than one refuge located, corresponding to the parapatric ranges of haplotype groups currently there. On the old western Mediterranean islands Corsica and Sardinia no differentiated (E. orbicularis) or only weakly differentiated haplotypes (T. hermanni) occur, even though there is evidence for the presence of both species on Corsica since at least the Middle Pleistocene. High mountain chains constitute major barriers separating distinct mtDNA clades or groups in each species.     

Turkish honeybees: genetic variation and evidence for a fourth lineage of Apis mellifera mtDNA

The mtDNA of bees from 84 colonies of Turkish honeybees (Apis mellifera) was surveyed for variation at four diagnostic restriction sites and the sequence of a noncoding intergenic region. These colonies came from 16 locations, ranging from European Turkey and the western Mediterranean coast to the Caucasus Mountains along the Georgian border, the eastern Lake Van region, and the extreme south. Combined restriction site and sequence data revealed four haplotypes. Three haplotypes belonged to the eastern Mediterranean mtDNA lineage. The fourth haplotype, which had a novel restriction site pattern and noncoding sequence, was found in samples from the extreme south, near the Syrian border. We found two different noncoding sequences among the eastern Mediterranean haplotypes. The "Caucasian" sequence matches that described from A. m. caucasica, and the "Anatolian" sequence matches that of A. m. carnica. The frequency of the "Caucasian" sequence was highest (98-100%) in sites near the Georgian border and decreased steeply to the south and west. Elsewhere the Anatolian sequence was found. In European Turkey (Thrace) a restriction site polymorphism previously reported from A. m. carnica in Austria and the Balkans was present at high frequency. A novel mtDNA haplotype with a unique restriction site pattern and noncoding sequence was found among bees from Hatay, in the extreme south near the Syrian border. This haplotype differed from the three previously known lineages of honeybee mtDNA--African, western European, and eastern Mediterranean-and may represent a fourth mitochondrial lineage.     

mtDNA variation of aboriginal Siberians reveals distinct genetic affinities with Native Americans.

The mtDNA variation of 411 individuals from 10 aboriginal Siberian populations was analyzed in an effort to delineate the relationships between Siberian and Native American populations. All mtDNAs were characterized by PCR amplification and restriction analysis, and a subset of them was characterized by control region sequencing. The resulting data were then compiled with previous mtDNA data from Native Americans and Asians and were used for phylogenetic analyses and sequence divergence estimations. Aboriginal Siberian populations exhibited mtDNAs from three (A, C, and D) of the four haplogroups observed in Native Americans. However, none of the Siberian populations showed mtDNAs from the fourth haplogroup, group B. The presence of group B deletion haplotypes in East Asian and Native American populations but their absence in Siberians raises the possibility that haplogroup B could represent a migratory event distinct from the one(s) which brought group A, C, and D mtDNAs to the Americas. Our findings support the hypothesis that the first humans to move from Siberia to the Americas carried with them a limited number of founding mtDNAs and that the initial migration occurred between 17,000-34,000 years before present.     

Low mitochondrial DNA variation and recent colonization of Scandinavia by the wood lemming Myopus schisticolor

The population genetic structure of wood lemmings ( Myopus schisticolor ) from Scandinavia, Finland and western Siberia was examined by restriction fragment length polymorphism of mitochondrial DNA (mtDNA) in 45 individuals from six localities. The 12 observed mtDNA haplotypes demonstrated a distinct phylogeographic pattern, suggesting that the postglacial colonization of Scandinavia by wood lemmings occurred from north-east. However, a very low level of haplotype and nucleotide diversity and a lack of geographical structure were found within Scandinavia. The limited mtDNA diversity in the Scandinavian populations probably reflects recent divergence in situ after colonization by a limited number of founders. Allozyme data support this scenario.     

Comparative study of carbohydrate chains released from the oviducal mucins of the two very closely related amphibian species Bombina bombina and Bombina variegata

The eggs of amphibians are surrounded by an extracellular matrix, termed jelly coat, which is mainly composed of hydrated mucin-type glycoproteins. These highly glycosylated molecules are synthesized by the oviduct and play an important role in the fertilization process. From a structural and chemical point of view, these oviducal mucins are very different from one species to another and they could be involved in the species-specificity of gamete interactions or could influence the parasite tropism. Bombina bombina and Bombina variegata are the two most closely related species within the genus, which hybridize readily in nature. Divergence occurred during geographic isolation estimated at 2-7 million years ago. The oviducal mucins of these species have been studied at the carbohydrate level, and the primary structures of 28 compounds have been established by NMR spectroscopy. The carbohydrate chains released from the oviducal mucins of the two species were similar and characterized by the common sequences GlcNAc(beta 1-3)[Fuc(alpha 1-4)]GlcNAc(beta 1-6) and GlcNAc(alpha 1-4)Gal(beta 1-4)Gal(beta 1-3) attached to GalNAc-ol (core 2). Nevertheless, some differences confirmed the strict species-specificity of amphibian oviducal carbohydrate chains observed previously. On the one hand, the presence of beta Gal 1,4-linked to beta GlcNAc in B. bombina, but not in B. variegata, can indicate that beta 4GalT: beta GlcNAc and beta 4GalT: beta Gal are two distinct glycosyltransferases. On the other hand, deaminoneuraminic acid (Kdn) is present in B. bombina, and N -glycolylneuraminic acid (NeuGc) in B. variegata. Although the enzymes involved in the biosynthesis of Kdn are not as well characterized, it can be suggested that at least one step of the biosynthetic pathway of NeuAc has been disrupted, leading the B. bombina oviducal NeuAc-9-synthase to use Man-6-P as a substrate, instead of ManNAc-6-P.     

MITOCHONDRIAL DNA VARIABILITY OF STRIPED DOLPHINS (STENELLA COERULEOALBA) IN THE SPANISH MEDITERRANEAN WATERS

Frozen muscle samples from 44 striped dolphins stranded on the Spanish Mediterranean coasts from 1990 to 1993 have been studied by means of mitochondrial DNA (mtDNA) restriction site analysis. Thirty-five of these dolphins were affected by a die-off occurring during this time in the western Mediterranean Sea. The mtDNA from each dolphin was digested with 15 restriction endonucleases that recognized 61 different restriction sites. The specific location of these sites on the mitochondrial gene map allowed us to determine the distribution of variability along this molecule. From the restriction analysis, a total of 15 different composite patterns or haplotypes was obtained and their phylogenetic relationships established both by cladistic and phenetic methods.
Estimates of sequence divergence among the 15 haplotypes ranged from 0.148% to 0.919%. Haplotype diversity (average heterozygosity) obtained was 0.789, a high value compared to other available data on cetaceans, especially considering that only a small area of the distribution of the species has been sampled. Neither population subdivision nor any evidence of correlation between haplotype frequencies and temporal or geographic distribution of dolphins has been detected.     

Phylogeny and biogeography of the alpine newt Mesotriton alpestris (Salamandridae, Caudata), inferred from mtDNA sequences

In this paper, we performed phylogenetic analyses of Mesotriton alpestris populations from the entire range of species distribution, using fragments of two mtDNA genes, cytochrome b (309bp) and 16S rRNA ( approximately 500bp). Sequence diversity patterns and phylogenetic analyses reveal the existence of a relict lineage (Clade A) of late Miocene origin, comprising populations from south-eastern Serbia. This lineage is proposed to be ancestor to a western and an eastern lineage, which diverged during the middle Pliocene. The western lineage is further divided in two clades (Clades B, C) of middle Pliocene origin that represent populations from Italy (B) and populations from central Europe and Iberia (C). Further subdivision, dated back to the middle-late Pliocene, was found within the eastern lineage, representing southern (Clade D) and central-northern (Clade E) Balkan populations, respectively. Extensive sequence divergence, implying greater isolation in multiple refugia, is found within eastern clades, while the western clades seem to have been involved in the colonization of central, western and north-eastern Europe from a hypothetical refugium in central Europe. The extent of divergence does not support the current taxonomy indicating cryptic speciation in the Balkans, while paedomorphic lineages were found to have been evolved during early-middle Pleistocene probably as a response to the ongoing dramatic climatic oscillations.     

Mitochondrial DNA phylogeography in northern grasshopper mice (Onychomys leucogaster) — the influence of Quaternary climatic oscillations on population dispersion and divergence

The details of mitochondrial DNA (mtDNA) phylogenetic structure of the northern grasshopper mouse Onychomys leucogaster were examined using populations from a postulated area of endemism that includes three arid regions (Colorado Plateaus, Interior Plains, and Wyoming Basins) in western North America. Fifteen tetra- and heptanucleotide restriction enzymes were used to assay restriction-site variation in a 2150-bp PCR-amplified fragment of mtDNA representing the ND2 and part of the COI gene regions. A total of 18 mtDNA haplotypes were detected. Although overall genetic divergence among these haplotypes was low (average = 1.1%), phylogeographic structuring was apparent. Notably, a clear phylogenetic split separated one group of haplotypes restricted to the Wyoming Basins from all others. This phylogenetic split was further corroborated by examination of nucleotide sequence variation from a 270-bp stretch of the mtDNA cytochrome b gene. Overall geographic and phylogenetic patterns suggest a complex history of geographic structuring and subsequent mixing of populations of grasshopper mice throughout the late Pleistocene. These patterns of variation are evaluated relative to alternative hypotheses about biotic responses to Quaternary climatic oscillations in western North American arid regions.     

Hybridization of Bombina bombina and B. variegata (Anura, Discoglossidae) at a sharp ecotone in western Ukraine: comparisons across transects and over time

Abstract Bombina bombina and B. variegata are two anciently diverged toad taxa that have adapted to different breeding habitats yet hybridize freely in zones of overlap where their parapatric distributions meet. Here, we report on a joint genetic and ecological analysis of a hybrid zone in the vicinity of Stryi in western Ukraine. We used five unlinked allozyme loci, two nuclear single nucleotide polymorphisms and a mitochondrial DNA haplotype as genetic markers. Parallel allele frequency clines with a sharp central step occur across a sharp ecotone, where transitions in aquatic habitat, elevation, and terrestrial vegetation coincide. The width of the hybrid zone, estimated as the inverse of the maximum gradient in allele frequency, is 2.3 km. This is the smallest of four estimates derived from different clinal transects across Europe. We argue that the narrow cline near Stryi is mainly due to a combination of habitat distribution and habitat preference. Adult toads show a preference for either ponds ( B. bombina ) or puddles ( B. variegata ), which is known to affect the distribution of genotypes within the hybrid zones. At Stryi, it should cause a reduction of the dispersal rate across the ecotone and thus narrow the cline. A detailed comparison of all five intensively studied Bombina transects lends support to the hypothesis that habitat distribution plus habitat preference can jointly affect the structure of hybrid zones and, ultimately, the resulting barriers to gene flow between differentiated gene pools. This study also represents a resampling of an area that was last studied more than 70 years ago. Our allele-frequency clines largely coincide with those that were described then on the basis of morphological variation. However, we found asymmetrical introgression of B. variegata genes into B. bombina territory along the bank of a river.     

A new cryptic species of pond turtle from southern Italy, the hottest spot in the range of the genus Emys (Reptilia, Testudines, Emydidae)

Geographic variation in the mtDNA haplotypes (cytochrome b gene) of 127 European pond turtles from Italy was investigated. Thirty‐eight of the Italian samples were also studied by nuclear fingerprinting (ISSR PCR) and compared with samples from other parts of the range representing all nine currently known mtDNA lineages of Emys orbicularis. Our genetic findings were compared against morphological data sets (measurements, colour pattern) for 109 adult turtles from southern Italy. Italy is displaying on a small geographical scale the most complicated variation known over the entire distributional area of Emys (North Africa over Europe and Asia Minor to the Caspian and Aral Seas). The Tyrrhenic coast of the Apennine Peninsula, the Mt. Pollino area and Basilicata are inhabited by Emys orbicularis galloitalica, a subspecies harbouring a distinct mtDNA lineage. The same lineage is also found in Sardinia. Along the Adriatic coast of Italy and on the Salentine Peninsula (Apulia, southern Italy), another morphologically distinctive subspecies (Emys orbicularis hellenica) occurs, which also bears a different mtDNA lineage. A higher diversity of mtDNA haplotypes in the south of the Apennine Peninsula suggests that the glacial refugia of E. o. galloitalica and E. o. hellenica were located here. A further refuge of E. o. hellenica probably existed in the southern Balkans. The west coasts of the Balkans and Corfu have probably been colonized from Italy and not from the geographically closer southern Balkanic refuge. In Sicily, a third mtDNA lineage is distributed, which is sister to all other known lineages of Emys. Morphologically, Sicilian pond turtles resemble E. o. galloitalica. However, nuclear fingerprinting revealed a clear distinctiveness of the Sicilian taxon, whereas no significant divergence was detected between representatives of the other eight mtDNA lineages of Emys. Furthermore, nuclear fingerprinting provided no evidence for current or past gene flow between the Sicilian taxon and the mainland subspecies of E. orbicularis. Therefore, Sicilian pond turtles are described here as a species new to science. Some populations in Calabria and on the Salentine Peninsula comprise individuals of different mtDNA lineages. We interpret this as a natural contact. However, we cannot exclude that these syntopic occurrences are the result of human activity. For example, in other parts of Italy, the natural distribution pattern of Emys is obscured by allochthonous turtles. This could also be true for southern Italy. The discovery of the complex taxonomic differentiation in southern Italy requires reconsidering conservation strategies.     

DIFFERING LEVELS OF AMONG-POPULATION DIVERGENCE IN THE MITOCHONDRIAL DNA OF PERIODICAL CICADAS RELATED TO HISTORICAL BIOGEOGRAPHY

Mitochondrial DNA (mtDNA) haplotypes were determined for 118 individuals of 13-and 17-year periodical cicadas (genus Magicicada) collected from 16 localities throughout the Midwest and eastern United States. Two distinct mtDNA lineages, identified as A and B, differ by 2.5% based on analysis of fragment patterns and restriction maps. Observed levels of mtDNA diversity within each lineage are low compared to estimates for other taxa. The two lineages are regionally segregated, with the boundary line occurring at a latitude of approximately 33° North. The levels of mtDNA diversity and population genetic structure differ within the two lineages. There is a remarkably low level of mean mtDNA divergence and no genetic structure in lineage A, whereas lineage B exhibits an order of magnitude higher level of mtDNA diversity and significant genetic structure among sampled populations. The low level of mtDNA diversity in cicadas may be attributed to (1) a population bottleneck that most likely occurred during the Pleistocene, (2) recent colonization following the retreat of the glaciers and the expansion of deciduous forests, and/or (3) high among-family reproductive variance (as a consequence of large population size, high fecundity, aggregative behavior of adults, and clumping of eggs). The difference in mtDNA diversity and population genetic structure between the lineages suggests that they experienced different biogeographic histories; we relate this to Pleistocene changes.     

The divergence of the dolly varden char Salvelinus malma in Asian Northern Pacific populations inferred from the PCR-RFLP analysis of the mitochondrial DNA

Genetic differentiation of the dolly varden char Salvelinus malma Walbaum was studied in five populations from the western part of the Northern Pacific. Using restriction analysis (RFLP), we examined polymorphism of three mitochondrial DNA (mtDNA) fragments amplified in polymerase chain reaction (PCR). MtDNA haplotypes were shown to fall into two phylogenetic groups, which probably reflect the existence of two previously described subspecies of Asian dolly varden, S. malma malma and S. malma krascheninnikovi. The divergence of mtDNA nucleotide sequences in the dolly varden subspecies (about 4%) corresponds to the differences between the valid char species from the genus Salvelinus.