Diverse historical processes shape deep phylogeographical divergence in the pollinating seed parasite Greya politella

Understanding the historical framework in which species interactions have diversified across landscapes may help to partition the effects of vicariance and geographically variable selection in shaping the geographical mosaic of coevolving species. We used phylogeographical analyses of the pollinating seed parasite Greya politella (Lepidoptera: Prodoxidae) to define the historical processes that may have structured interactions of this species with its host plants across major biogeographical breaks in western North America. Using 648 bp of cytochrome oxidase I and amplified fragment length polymorphisims, we identified deep genetic breaks among some populations consistent with some definitions of cryptic species. A combination of phylogenetic and population genetic approaches indicates that different historical processes may have structured G. politella genetic diversity in four regions: northern Pacific Northwest, southern Oregon, southern Sierra Nevada, and the remainder of California. The northern Pacific Northwest had high genetic diversity likely due to glacial refugia and subsequent spatial expansion, concordant with some other taxa. Populations in southern Oregon possessed unique, closely related haplotypes with restricted gene flow, possibly indicating a long-standing set of populations in this endemic-rich region. Analyses of California populations showed evidence of restricted gene flow and spatial expansion with many closely related haplotypes that occupy a broad geographical range. Southern Sierra Nevada populations were genetically distinct and highly diverse, possibly due to a localized glacial refugium. Together, these results suggest that vicariance and population expansion, possibly in combination with geographically variable selection, have shaped the diversification of G. politella and its interactions with its host plants.     

Phylogenetic relationships of Eurasian Nuthatches (Sitta europaea Linnaeus, 1758) from the Alborz and Zagros Mountains,Iran

The Eurasian Nuthatch (Sitta europaea Linnaeus, 1758) is a resident bird in the Alborz and Zagros deciduous forests. We investigated the phylogenetic relationships and the taxonomic position of the Eurasian Nuthatch among other separated lineages of Eurasia with the help of blood samples collected from 19 individuals belonging to four populations in the Eastern and Western Alborz, as well as in the Northern and Southern Zagros forests. Genetic variation was then analysed using complete ND2 gene sequence (1041?bp) and phylogenetic analysis was done using Bayesian and maximum likelihood inference. Additionally, a median-joining algorithm was used to reveal the relationships among haplotypes. The results of the phylogenetic and haplo-type network analyses indicated that Eurasian Nuthatch haplotypes from the Alborz and Zagros Mountains form lineages distinct from the Asian, Caucasian and European haplotypes. Furthermore, an analysis of molecular variance (AMOVA) detected significant (P<0.001) genetic structure among the lineages. The Asian, European, Caucasian and Alborz lineages diverged from one another by an uncorrected genetic distance ranging from 0.029 to 0.039, while the Zagros lineage showed a slightly lower genetic divergence from the Caucasian lineage (0.006), but it did not share any haplotype with the Caucasian lineage. Thus, we suggest considering five Conservation Significant Units (CSU) for the Eurasian Nuthatches as the result of used dataset.     

Genetic variation of European grayling (Thymallus thymallus) populations in the Western Balkans

In order to elucidate genetic composition of European grayling (Thymallus thymallus) populations in the Western Balkans, the partial mitochondrial DNA (mtDNA) control region was sequenced and 12 microsatellite loci genotyped in 14 populations originating from tributaries of the Adriatic and Danube drainages. Eleven mtDNA haplotypes were found, one confined to the Adriatic clade, one to the Alpine group and the rest to the ‘Balkan’ grayling phylogenetic clade. Haplotypes from the Balkan clade were confined to the Danube drainage and constituted two groups: northern group with haplotypes found in the Slovenian part of the Danube drainage, and southern group, consisting from Bosnia–Herzegovina and Montenegro. Substantial genetic distance between northern and southern groups of haplotypes (0.75–1.8%) and well supported divisions within the northern group indicate very structured grayling population within the studied Danube basin that most probably did not evolve due to vicariance but rather as a consequence of multiple colonization waves that might have occurred during the Pleistocene. Furthermore, genetic distance of ~4% between Adriatic and Danube populations’ haplotypes, suggest that their separation occurred in mid-Pliocene. These findings imply a complex colonization pattern of the Western Balkans drainages. Microsatellite data also confirm high genetic diversity in Western Balkans populations of grayling (on average 7.5 alleles per microsatellite locus and H _exp 0.58). Limited stocking activities were detected based on microsatellites and mtDNA data. Regarding current knowledge of grayling phylogeography appropriate management strategies were proposed to preserve unique, autochthonous grayling populations in Western Balkan.     

Contrasting levels of connectivity and localised persistence characterise the latitudinal distribution of a wind-dispersed rainforest canopy tree

Contrasting signals of genetic divergence due to historic and contemporary gene flow were inferred for Coachwood, Ceratopetalum apetalum (Cunoniaceae), a wind-dispersed canopy tree endemic to eastern Australian warm temperate rainforest. Analysis of nine nuclear microsatellites across 22 localities revealed two clusters between northern and southern regions and with vicariance centred on the wide Hunter River Valley. Within populations diversity was high indicating a relatively high level of pollen dispersal among populations. Genetic variation was correlated to differences in regional biogeography and ecology corresponding to IBRA regions, primary factors being soil type and rainfall. Eleven haplotypes were identified by chloroplast microsatellite analysis from the same 22 localities. A lack of chloroplast diversity within sites demonstrates limited gene flow via seed dispersal. Network representation indicated regional sharing of haplotypes indicative of multiple Pleistocene refugia as well as deep divergences between regional elements of present populations. Chloroplast differentiation between sites in the upper and lower sections of the northern population is reflective of historic vicariance at the Clarence River Corridor. There was no simple vicariance explanation for the distribution of the divergent southern chlorotype, but its distribution may be explained by the effects of drift from a larger initial gene pool. Both the Hunter and Clarence River Valleys represent significant dry breaks within the species range, consistent with this species being rainfall dependent rather than cold-adapted.     

Evolutionary relationship between marble trout of the northern and the southern Adriatic basin

Marble trout (Salmo marmoratus) populate two geographically separated areas in the northern and southern parts of the Adriatic Sea drainage. Although morphologically similar, each population is distinguished by a different set of unrelated mitochondrial haplotypes, suggesting that they have evolved from different ancestors. Due to a possible discordance between mitochondrial and species phylogeny, we performed phylogenetic analysis based on 22 nuclear loci. The results inferred from Maximum-likelihood and Bayesian Inference analysis revealed that northern and southern populations are closely related, forming a monophyletic group. This observation is concordant with the present marble trout classification, which considers both populations as conspecific. On the other hand, our findings are in marked contrast to those of previous mtDNA-based studies and highlight potential dangers of making phylogenetic inferences from mtDNA alone. Reasons for discordance between mtDNA and nDNA phylogeny are discussed with incomplete lineage sorting proposed as the most parsimonious explanation for mtDNA divergence in marble trout.     

Genetic diversity and phylogenetic origin of brown trout Salmo trutta populations in eastern Balkans

The study focuses on the phylogenetic origin and genetic diversity of brown trout in the eastern part of the Balkan Peninsula. It further aims to reveal the impact of human-mediated transfers and stocking with non-indigenous trout on the populations in this area. For these purposes, mtDNA control region and microsatellite variation of 204 individuals from 16 populations were analysed. The results indicate that mtDNA haplotypes from the lower Danube basin and southern Black Sea basins differ substantially from a subclade of the Danubian lineage consisting of haplotypes found so far in the most of the Danube basin and in the Caspian and Aral Sea basins. Considering also the results of demographic analyses, this study evidences a complex evolutionary history of brown trout in the southern and western parts of the Black Sea basin. In the Aegean Sea basin, a high frequency of the central haplotype of Adriatic mtDNA lineage has been found. The other Adriatic lineage haplotypes found in this basin differ from the central haplotype by one mutational step only, indicating a recent evolution of the Adriatic lineage in the Aegean Sea basin. Substantial genetic differentiation among populations and basins was revealed. The hybridization with Atlantic brown trout was indicated in both sea basins, but especially in the Danube basin. Compared to other European regions, it can be inferred that the introgression of exogenous brown trout in the eastern Balkan populations is rather low.     

Mountains too high and valleys too deep drive population structuring and demographics in a Qinghai–Tibetan Plateau frog Nanorana pleskei (Dicroglossidae)

Pleistocene glacial–interglacial climatic oscillations greatly shaped the current genetic structure of many species. However, geographic features may influence the impact of climatic cycling. Distinct geographic and environmental characters between northern and southern parts of the eastern Qinghai–Tibetan Plateau (EQTP) facilitate explorations into the impacts of geographic features on species. The northern parts of EQTP contain large areas of marsh, and the environment is rather homogeneous. In contrast, the southern EQTP harbors complex alpine valleys and a much more heterogeneous setting. We evaluate DNA sequence variation from both the mitochondrial and nuclear genomes in Nanorana pleskei, a species endemic to the EQTP. Hypothesis testing on the evolutionary history of N. pleskei indicates that northern populations can disperse freely, but alpine valleys isolate southern populations. Demographic histories between northern and southern populations also differ. Northern populations appear to have experienced population expansions, while southern frogs exhibit a far more stable demographic history. By combining climatic analyses and species' distribution models, our study suggests that geographic and environmental features drive the differences between the northern and southern EQTP.     

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.     

Molecular Phylogeny and Genetic Divergence of the Introduced Populations of Japanese Guppies, Poecilia reticulata

A phylogenetic tree of guppies from seven introduced Japanese populations was constructed using a mitochondrial gene, and the relationship between these Japanese populations and indigenous South American populations was examined. Phylogenetic trees constructed by maximum parsimony and neighbor joining methods indicated four major groups: Group A: Trinidadian populations; Group B: the northernmost Okinawa populations; Group C: the northern Okinawa populations; Group D: Shimoda populations and the remainder of Trinidadian, Venezuelan, and southern Okinawa populations. Considerable genetic differences were observed among the haplotypes within each Okinawa population similar to that found among haplotypes of different native populations, but not within the Shimoda population, because each Okinawa population included haplotypes belonging to different groups. Genetic differences between Groups B, C, and D might reflect those of the introduced lineages; however, these differences might not result from the divergence after introduction events. These results suggested that multiple introductions of different breeding lineages might cause considerable divergence within and between Okinawa populations.     

The wide but disjunct range of the European mountain plant Androsace lactea L. (Primulaceae) reflects Late Pleistocene range fragmentation and post‐glacial distributional stasis

Aim Our aim was to reconstruct the spatio‐temporal genetic diversification of Androsace lactea, a widely but disjunctly distributed European mountain plant, to test the hypothesis that its distribution is the result of vicariance, in the late Tertiary or during the Pleistocene, or alternatively of long‐distance dispersal. We also addressed the phylogeographic history of the Alps, emphasizing the role of Pleistocene refugia at their northern margin. Location The central and southern European mountain ranges. Methods We gathered amplified fragment length polymorphism (AFLP) data and plastid DNA sequences from one to four individuals of each of 26 populations spanning the entire distribution area. AFLP data were analysed with Bayesian clustering approaches, neighbour‐joining analysis and NeighbourNet. Plastid sequences were used to depict relationships among haplotypes in a statistical parsimony network, to test for population expansions, and to obtain age estimates in a Bayesian framework. Results The AFLP data suggested that many populations were genetically strongly differentiated. The internal structure, however, was weak, and only two major groups of populations, from the north‐western Alps and adjacent regions and from the easternmost Alps, were supported in the neighbour‐joining analysis. One of the Bayesian clustering approaches differentiated three groups of populations: Northern Alps, easternmost Alps and the remaining distribution area. Eleven closely related plastid haplotypes were found, separated by maximally four mutational steps, resulting in a star‐like parsimony network. None of several estimators suggested statistically significant population expansions. The diversification age was inferred to be (mean/median) 0.135/0.08 Ma (95% highest posterior density interval 0.364–0.006 Ma). Main conclusions We found no evidence that long‐distance dispersal shaped the disjunct distribution range; our data rather favoured a vicariance scenario. However, in contrast to the hypothesis that wide but disjunct distributions are old, we conclude that range fragmentation probably happened in the Late Pleistocene, perhaps during the last glaciation. In the Alps, most populations are at least close to formerly unglaciated areas. Our data support distributional stasis and suggest that important refugia were situated at the north‐eastern, but also at the northern and north‐western edges of the Alps, thereby strengthening the evidence for glacial refugia in this strongly glaciated region.     

Reticulated Origin of Domesticated Emmer Wheat Supports a Dynamic Model for the Emergence of Agriculture in the Fertile Crescent

We used supernetworks with datasets of nuclear gene sequences and novel markers detecting retrotransposon insertions in ribosomal DNA loci to reassess the evolutionary relationships among tetraploid wheats. We show that domesticated emmer has a reticulated genetic ancestry, sharing phylogenetic signals with wild populations from all parts of the wild range. The extent of the genetic reticulation cannot be explained by post-domestication gene flow between cultivated emmer and wild plants, and the phylogenetic relationships among tetraploid wheats are incompatible with simple linear descent of the domesticates from a single wild population. A more parsimonious explanation of the data is that domesticated emmer originates from a hybridized population of different wild lineages. The observed diversity and reticulation patterns indicate that wild emmer evolved in the southern Levant, and that the wild emmer populations in south-eastern Turkey and the Zagros Mountains are relatively recent reticulate descendants of a subset of the Levantine wild populations. Based on our results we propose a new model for the emergence of domesticated emmer. During a pre-domestication period, diverse wild populations were collected from a large area west of the Euphrates and cultivated in mixed stands. Within these cultivated stands, hybridization gave rise to lineages displaying reticulated genealogical relationships with their ancestral populations. Gradual movement of early farmers out of the Levant introduced the pre-domesticated reticulated lineages to the northern and eastern parts of the Fertile Crescent, giving rise to the local wild populations but also facilitating fixation of domestication traits. Our model is consistent with the protracted and dispersed transition to agriculture indicated by the archaeobotanical evidence, and also with previous genetic data affiliating domesticated emmer with the wild populations in southeast Turkey. Unlike other protracted models, we assume that humans played an intuitive role throughout the process.     

Distinguishing between primary and secondary intergradation among morphologically differentiated populations of Anolis marmoratus

Distinguishing between primary and secondary intergradation among differentiated populations, and the relative importance of drift and selection, are persistent problems in evolutionary biology. An historical perspective on population interactions can provide insight into the nature of contacts, and thus help resolve these questions. Continuously distributed populations of Anolis marmoratus from the island of Basse Terre in the Guadeloupean archipelago of the Lesser Antilles show a striking degree of geographic variation in morphology. Initial surveys of mtDNA variation from throughout the Guadeloupean Archipelago revealed one case where levels of sequence difference and phylogenetic relationships of alleles from morphologically differentiated populations from the east coast of Basse Terre were consistent with primary intergradation. In this paper, I examine the genetic population structure of a series of populations spanning this north-south cline in morphological variation to test the hypothesis of primary intergradation. Sequences of the mitochondrial cytochrome-b gene from 50 individuals representing five populations spanning the cline were obtained and fourteen unique haplotypes (differing by 2% or less) were detected. Patterns of nucleotide substitution among haplotypes do not deviate from neutral expectation indicating no effect of selection at the level of mtDNA sequences. Estimates of population structure and gene flow were made using both summary statistics for nucleotide diversity (N_at) and cladistic methods. The results are sensitive to the choice of gene flow model, and this is discussed in detail. Mitochondrial variation in the northern populations may not be at equilibrium, and the phylogeny of alleles is consistent with a recent increase in effective population size. Estimates of nucleotide diversity, gene flow, and the phylogenetic relationships of haplotypes indicate that the southern-most population (representing the extreme of morphological variation along this cline) has been relatively isolated from populations to the north and has experienced a reduced effective population size. The apparent clinal variation between the southern population and the others may therefore reflect secondary contact and introgression rather than primary intergradation.     

Molecular phylogeography of endangered sharp-snouted pitviper (Deinagkistrodon acutus; Reptilia, Viperidae) in Mainland China

Using phylogenetic and population genetic approaches, the present study reports the phylogeographic structure of the sharp-snouted pitviper (Deinagkistrodon acutus), a threatened snake species with commercial and medicinal importance in China. The entire mitochondrial ND2 gene (NADH dehydrogenase subunit 2) sequences of 86 individuals of D. acutus from 14 localities across its range in China were determined. Based on the results of phylogenetic analyses, distribution of diagnostic sites, haplotype network, and AMOVA hierarchical analysis, an east-west division of the whole D. acutus population could be observed. Geographically, a line formed by a lake, river, and mountain chain (the Poyang Lake, Gan River to the southern end of the Wuyi Mountains), results in vicariance and approximately vertically splits the range into two and the whole population into two main lineages (western and eastern). The bifurcating tree suggested generally west to east dispersal trend. The data fit the isolation by distance (IBD) model well. Star-like clusters in haplotype network, significantly negative values of Fs statistics, and unimodal mismatch distributions all suggest recent demographic expansions in four areas. The results show that isolation, dispersal, bottleneck, and expansion jointly constitute the history of D. acutus. In a haplotype network, the excessive predominance of central haplotypes, few medium-frequency haplotypes, predominance (73.1%) of the singletons among the derived haplotypes, most of which are connected to the central haplotype by only one mutational step, unsymmetrical campanulate unimodal curve of mismatch distributions and leftwards shift of the peaks, all suggest that the whole D. acutus population is a young population with low genetic diversity. Based on the data, the first priority for conservation action should be given to the Huangshan unit.     

Phylogenetic placement and population structure of Indo‐Pacific bottlenose dolphins (Tursiops aduncus) off Zanzibar,Tanzania, based on mtDNA sequences

Phylogenetic placement of bottlenose dolphins from Zanzibar, East Africa and putative population differentiation between animals found off southern and northern Zanzibar were examined using variation in mtDNA control region sequences. Samples (n= 45) from animals bycaught in fishing gear and skin biopsies collected during boat surveys were compared to published sequences (n= 173) of Indo‐Pacific bottlenose dolphin, Tursiops aduncus, from southeast Australian waters, Chinese/Indonesian waters, and South African waters (which recently was proposed as a new species) and to published sequences of common bottlenose dolphin, Tursiops truncatus. Bayesian and maximum parsimony analyses indicated a close relationship between Zanzibar and South African haplotypes, which are differentiated from both Chinese/Indonesian and Australian T. aduncus haplotypes. Our results suggest that the dolphins found off Zanzibar should be classified as T. aduncus alongside the South African animals. Further, analyses of genetic differentiation showed significant separation between the T. aduncus found off northern and southern Zanzibar despite the relatively short distance (approximately 80 km) between these areas. Much less differentiation was found between southern Zanzibar and South Africa, suggesting a more recent common evolutionary history for these populations than for the northern and southern Zanzibar populations.     

Colonization history of north European field voles (Microtus agrestis) revealed by mitochondrial DNA

The genetic structure of field vole (Microtus agrestis) populations from northern Europe was examined by restriction fragment length polymorphisms of mitochondrial DNA (mtDNA) in 150 individuals from 67 localities. A total of 83 haplotypes was observed, most of which were rare and highly localized geographically. Overall nucleotide diversity was high (134%), but showed a tendency to decrease with higher latitude. Two major mtDNA lineages differing by 2% in nucleotide sequence were identified. A southern mtDNA lineage was observed in field voles from Britain, Denmark and southern and central Sweden, whereas voles from Finland and northern Sweden belonged to a northern lineage. The strict phylogeographic pattern suggests that the present population generic structure in field voles reflects glacial history: the two groups are derived from different glacial refugia, and recolonized Fennoscandia from two directions. A 150–200-km-wide secondary contact zone between the two mtDNA groups was found in northern Sweden. Distinct phylogeographic substructuring was observed within both major mtDNA groups.     

Phylogeography and ecological niche modelling in Eugenia uniflora (Myrtaceae) suggest distinct vegetational responses to climate change between the southern and the northern Atlantic Forest

In this study, we evaluate phylogeographic patterns and predictions of ecological niche modelling (ENM) for Eugenia uniflora (Myrtaceae), a widely distributed taxon in the Atlantic forest domain, to understand the effect of past climatic oscillations on the demographic history of this species. An analysis of phylogeographic population structure and demography was conducted on E. uniflora from 46 localities in natural environments across the distribution range of the species based on three plastid markers. ENM was also performed to predict suitable environments and areas of dramatic decrease in future suitability for the species under distinct representative concentration pathways (RCPs). Eugenia uniflora exhibited higher haplotype and nucleotide diversity in the southern part of its distribution than in the northern part. Two divergent lineages were revealed in the phylogenetic analysis of haplotypes, with an estimated divergence at c. 4.9 Mya. The populations in the northern and central regions of the range probably experienced population growth, whereas populations in the southern region are marked by historical demographic stability. ENM results indicate that the distribution of E. uniflora was fragmented in cool periods and was broader and more connected during warm periods during Pleistocene. The results suggest distinct evolutionary histories in southern to northern populations, indicating region‐specific responses to changes.     

Phylogeography of Podocarpus matudae (Podocarpaceae): pre‐Quaternary relicts in northern Mesoamerican cloud forests

Aim Cloud forests of northern Mesoamerica represent the northern and southern limit of the contact zone between species otherwise characteristic of North or South America. Several phylogeographic studies featuring temperate conifer species have improved our understanding of species responses to environmental changes. In contrast, conifer species that presumably colonized northern Mesoamerica from South America are far less studied. A phylogeographic study of Podocarpus matudae (Podocarpaceae) was conducted to identify any major evolutionary divergences or disjunctions across its range and to determine if its current distribution is associated with pre‐Quaternary climatic and/or long‐distance dispersal events. Location Northern Mesoamerica (Mexico and Guatemala). Methods Sixteen populations (157 individuals) of P. matudae were screened for variation at two plastid DNA markers. The intra‐specific phylogenetic relationships among haplotypes were reconstructed using Bayesian inference. Population genetic analyses were undertaken to gain insight into the evolutionary history of these populations. To test whether genetic divergence among populations occurred at different time‐scales plastid DNA sequence data and fossil‐ and coalescent‐based calibrations were integrated. Results The combination of plastid markers yielded 11 haplotypes. Differentiation among populations based on DNA variation (G_ST) (0.707, SE 0.0807) indicated a clear population structure in P. matudae. Differentiation for ordered alleles (N_ST) (0.811, SE 0.0732) was higher than that for G_ST, indicating phylogeographical structure in P. matudae. Most of the total variation (81.3%, P < 0.0001) was explained by differences among populations. The estimated divergence time between the unique haplotypes from a Guatemalan population and the two most common haplotypes from the Sierra Madre Oriental in Mexico was between 10 and 20 Ma, and further haplotype divergence in the poorly resolved clade of the Sierra Madre Oriental occurred between 3 and 0.5 Ma. Main conclusions Divergence estimations support the hypothesis that extant Podocarpus matudae populations are pre‐Quaternary relicts. This finding is consistent with fossil and pollen data that support a Miocene age for temperate floristic elements in Mesoamerican cloud forests, whereas further haplotype divergence within the Sierra Madre Oriental, Chiapas and Guatemala occurred more recently, coinciding with Pleistocene cloud forest refugia.     

Gene conversion in the mitochondrial genome on interspecific hybridization in voles of the <Emphasis Type="Italic">Clethrionomys</Emphasis> genus

The phenomenon of interspecific hybridization accompanied by transfer of the mitochondrial genome from the northern red-backed vole (Clethrionomys rutilus) to the bank vole (Cl. glareolus) in northeastern Europe is well known already for 25 years. However, the possibility of recombination between homologous segments of maternal and paternal mtDNAs of the voles during fertilization was not previously studied. Analysis of data on variability of nucleotide sequences of the mitochondrial gene for cytochrome b in populations of red-backed and bank voles in the area of their sympatry has shown that as a result of interspecific hybridization, the mitochondrial gene pool of bank voles contains not only mtDNA haplotypes of red-backed vole females, but also mtDNA haplotypes of bank voles bearing short nucleotide tracts of red-backed vole mtDNA. This finding supports the hypothesis that an incomplete elimination of red-backed vole paternal mtDNA during the interspecific hybridization between bank vole females and red-backed vole males leads to the gene conversion of bank vole maternal mtDNA tracts by homologous ones of mtDNA of red-backed vole males.