Phylogenetic divergence in leatherside chub (Gila copei) inferred from mitochondrial cytochrome b sequences

We examined intra-specific phylogenetic relationships in leatherside chub, Gila copei. The complete mitochondrial (mt) cytochrome b gene (1140 bp) was sequenced for 30 individuals from 10 populations that span the geographical distribution of this species. Traditional phylogenetic analyses revealed two deeply divergent and evolutionarily distinct mtDNA clades that are geographically separated in northern and southern drainage basins. Interpopulation sequence variation between clades ranged from 7.7 to 8.1%. The northern clade was genetically more similar and phylogenetically more closely related to the selected out-group Lepidomeda m. mollispinus than to the southern clade, suggesting that the taxonomy of this species may require revision. Sequence variation among populations within clades ranged from 0 to 0.3% in the north and from 0 to 0.7% in the south. Statistical parsimony was used to construct phylogenetic networks of haplotypes within clades. Nested clade analysis revealed that geographical fragmentation has played an important role in genetic structuring within northern and southern clades.     

Phylogeography of two European newt species--discordance between mtDNA and morphology

The newts Triturus vulgaris and Triturus montandoni are sister species that exhibit contrasting levels of intraspecific morphological variation. Triturus vulgaris has a broad Eurasiatic distribution encompassing both formerly glaciated and unglaciated areas and shows substantial morphological differentiation in the southern part of its range, while T. montandoni, confined to the Carpathians, is morphologically uniform. We analysed sequence variation of two mtDNA fragments of the total length of c. 1850 bp in 285 individuals of both species collected from 103 localities. Phylogenetic analysis of 200 unique haplotypes defined 12 major clades, their age estimated at c. 4.5-1.0 million years (Myr). Most of the older clades were found in the southern part of the range, and also in central Europe, mainly in Romania. The distribution of mtDNA clades points to the existence of several glacial refugia, located in the Caucasus region, Anatolia, the Balkan Peninsula, Italy, and more to the north in central Europe. The concordance between mtDNA based phylogeny and the distribution of T. vulgaris subspecies was weak. Triturus montandoni haplotypes did not form a monophyletic group. Instead they were found in six clades, in five of them mixed with T. vulgaris haplotypes, most likely as a result of past or ongoing hybridization and multiple introgression of mtDNA from T. vulgaris to T. montandoni. Patterns of sequence variation within clades suggested long-term demographic stability in the southern groups, moderate and relatively old demographic growth in the populations inhabiting central Europe, and high growth in some of the groups that colonized northern parts of Europe after the last glacial maximum.     

Genetic and geographical differentiation of <Emphasis Type="Italic">Pandaka</Emphasis> gobies in Japan

The mitochondrial DNA (mtDNA) phylogeny of Japanese Pandaka species (Perciformes: Gobiidae) was inferred from partial nucleotide sequences of the mitochondrial 12S and 16S rRNA genes (1083bp). The resultant mtDNA tree showed two major clades (clade I and clade II), which were inconsistent with the present taxonomic classification. One of the major clades was further divided into two geographical groups, distributed on the Japanese Major Islands (clade I-A) and from Amami-oshima Island to Iriomote Island (clade I-B). The mtDNA haplotypes in clade II were found only on Iriomote Island. The mtDNA divergences in clade I indicated that the Japanese Major Island (clade I-A) and Ryukyu (clade I-B) groups have been geographically isolated from each other for millions of years, based on the putative molecular divergence rate. The geographical distributions of mtDNA haplotypes in clade I-A and clade I-B also suggested that Pandaka gobies had not dispersed to distant offshore islands, indicating that their geographical differentiation may be closely associated with the geological history of the Japanese and Ryukyu Archipelagos.This revised version was published online in January 2005 with corrections to the repetition of the 1st authors name.     

Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers

We compared mitochondrial DNA and gill-raker number variation in populations of the European whitefish Coregonus lavaretus (L.) species complex to illuminate their evolutionary history, and discuss mechanisms behind diversification. Using single-strand conformation polymorphism (SSCP) and sequencing 528 bp of combined parts of the cytochrome oxidase b (cyt b) and NADH dehydrogenase subunit 3 (ND3) mithochondrial DNA (mtDNA) regions, we documented phylogeographic relationships among populations and phylogeny of mtDNA haplotypes. Demographic events behind geographical distribution of haplotypes were inferred using nested clade analysis (NCA) and mismatch distribution. Concordance between operational taxonomical groups, based on gill-raker numbers, and mtDNA patterns was tested. Three major mtDNA clades were resolved in Europe: a North European clade from northwest Russia to Denmark, a Siberian clade from the Arctic Sea to southwest Norway, and a South European clade from Denmark to the European Alps, reflecting occupation in different glacial refugia. Demographic events inferred from NCA were isolation by distance, range expansion, and fragmentation. Mismatch analysis suggested that clades which colonized Fennoscandia and the Alps expanded in population size 24 500-5800 years before present, with minute female effective population sizes, implying small founder populations during colonization. Gill-raker counts did not commensurate with hierarchical mtDNA clades, and poorly with haplotypes, suggesting recent origin of gill-raker variation. Whitefish designations based on gill-raker numbers were not associated with ancient clades. Lack of congruence in morphology and evolutionary lineages implies that the taxonomy of this species complex should be reconsidered.     

Molecular phylogeography of common garter snakes (Thamnophis sirtalis) in western North America: implications for regional historical forces

Complete ND2 and partial ND4 and cytochrome b mitochondrial DNA (mtDNA) sequences were analysed to evaluate the phylogeographic patterns of common garter snakes (Thamnophis sirtalis) in western North America. This species is widely distributed throughout North America, and exhibits extensive phenotypic variation in the westernmost part of its range. The overall phylogeographic pattern based on mtDNA sequences is concordant with results from studies of other species in this region, implicating historical vicariant processes during the Pleistocene and indicating bottleneck effects of recent dispersal into postglacial habitat. Indeed, the topology is statistically consistent with the hypothesis of both southern (Great Basin and California) and northern (Haida Gwaii) refugia. Specifically, we identified genetic breaks among three major clades: Northwest Coastal populations, Intermountain populations, and all California populations. The California clade contained the only other well-supported branching patterns detected; relationships among populations within the two northern clades were indistinguishable. These molecular splits contrast sharply with all prior geographical analyses of phenotypic variation in T. sirtalis in this region. Our results suggest that the extensive phenotypic variation in western T. sirtalis has been shaped more by local evolutionary forces than by shared common ancestry. Consequently, we consider all morphologically based subspecies designations of T. sirtalis in this region invalid because they do not reflect reciprocal monophyly of the mtDNA sequences.     

Phylogeography of Sylvirana latouchii (Anura, Ranidae) in Taiwan

Biogeographic studies are important for understanding the natural history of faunas. To comprehend the geographical patterns of genetic variation in anurans in Taiwan, we investigated the genetic structure of Sylvirana latouchii (Anura, Ranidae) from 31 populations by using mitochondrial DNA (mtDNA) cytochrome b sequences. A neighbor-joining tree of 38 haplotypes revealed three major divergences in Taiwanese S. latouchii: the northern, western, and eastern-and-southern clades. Each clade was restricted to a single geographical district and showed obvious differentiation. The patterns of geographical divergence in this species reflect common historical events experienced by other native animals distributed in Taiwan. The order of divergence times between clades was inferred using a molecular clock test. The population relationship of S. latouchii between Taiwan and mainland China is discussed. Further study employing more populations of S. latouchii from mainland China is necessary to clarify the original geographical patterns and migratory history of this species.     

Combining molecular, morphological and ecological data to infer species boundaries in a cryptic tropical pitviper

Few operational methods exist for delimiting species boundaries, and these usually require sampling strategies that are unrealistic for widespread organisms that occur at low densities. Here we apply molecular, morphological and ecological species delimitation criteria to a wide-ranging, fragmented group of Asian green pitvipers, the Popeia popeiorum complex. A mitochondrial DNA phylogeny for the group indicates two well-differentiated clades, corresponding mainly to northern and southern parts of its range. Strong phylogeographical structure within each clade suggests isolation in forest refugia during the Pliocene and a southward colonization of the Sunda islands during the Pleistocene. Multivariate analysis of morphological characters reveals a generally conserved pattern of geographical variation, incongruent with the recovered phylogenetic history. We compare groups delineated by mtDNA variation to morphological and ecological divisions in the complex, and discuss the implications of these for the taxonomy of the group. Discordance between species boundaries inferred from different criteria suggests that combining independent sources of data provides the most reliable estimation of species boundaries in organisms that are difficult to sample in large numbers.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 87 , 343–364.     

Phylogeography of European moose (Alces alces) based on contemporary mtDNA data and archaeological records

Phylogeography can help to determine LGM refugia and postglacial migration routes. However, the locations of LGM refugial areas in eastern Europe are not clear. Moose (Alces alces) is presently a common species in central and north-eastern Europe, but there are no studies showing its phylogenetic pattern and genetic diversity across its whole continuous range. Moose never became extinct in the eastern part of its range, and the eastern mtDNA lineage has the largest effective population size. The present study shows the phylogeographic pattern and genetic diversity of European moose and compares the results of mtDNA analyses with the archaeological record of the species to identify its LGM refugia and postglacial migration routes. I combined the mtDNA control region sequences obtained in all studies of moose in Europe and western Asia. The genetic data were then compared with the archaeological records of the species dated to the LGM. I found that the European moose lineage inhabits Europe and western Asia. It is composed of two clades: the eastern and the central-western, consisting of a total of six discrete haplogroups. The most complex, the eastern clade, has the largest range. Some of the haplogroups have narrow or scattered distributions and two are common in almost the whole range. Genetic diversity hotspots were detected in contact zones of different mtDNA haplogroups rather than in the LGM refugial areas of moose. Archaeological records dated to the LGM were found in several localities in central, southern and eastern Europe as well as in western Asia. The range of the moose during the LGM was much larger than previously thought. The eastern clade survived the LGM in western Siberia, the Ural Mountains and Russian plain. LGM refugia of moose were also located in the Caucasus, Carpathians, Balkans and northern Italy.     

Genetic divergence and phylogeographic relationships among european perch (Perca fluviatilis) populations reflect glacial refugia and postglacial colonization

We used the widely distributed freshwater fish, perch (Perca fluviatilis), to investigate the postglacial colonization routes of freshwater fishes in Europe. Genetic variability within and among drainages was assessed using mitochondrial DNA (mtDNA) D-loop sequencing and RAPD markers from 55 populations all over Europe as well as one Siberian population. High level of structuring for both markers was observed among drainages and regions, while little differentiation was seen within drainages and regions. Phylogeographic relationships among European perch were determined from the distribution of 35 mtDNA haplotypes detected in the samples. In addition to a distinct southern European group, which includes a Greek and a southern Danubian population, three major groups of perch are observed: the western European drainages, the eastern European drainages including the Siberian population, and Norwegian populations from northern Norway, and western side of Oslofjord. Our data suggest that present perch populations in western and northern Europe were colonized from three main refugia, located in southeastern, northeastern and western Europe. In support of this, nested cladistic analysis of mtDNA clade and nested clade distances suggested historical range expansion as the main factor determining geographical distribution of haplotypes. The Baltic Sea has been colonized from all three refugia, and northeastern Europe harbours descendants from both eastern European refugia. In the upper part of the Danube lineages from the western European and the southern European refugia meet. The southern European refugium probably did not contribute to the recolonization of other western and northern European drainages after the last glaciation. However, phylogenetic analyses suggest that the southern European mtDNA lineage is the most ancient, and therefore likely to be the founder of all present perch lineages. The colonization routes used by perch probably also apply to other freshwater species with similar distribution patterns.     

Coalescence methods reveal the impact of vicariance on the spatial genetic structure of Elephantulus edwardii (Afrotheria, Macroscelidea)

Within the Macroscelidea 15 species of elephant-shrews are recognized, of which nine occur in the southern African subregion. The Cape rock elephant-shrew (Elephantulus edwardii) is the only strictly endemic South African elephant-shrew species. Recent distribution data suggest that E. edwardii is continuously distributed from Namaqualand in the Western Cape Province to Port Elizabeth in the Eastern Cape Province. Molecular sequences from the mitochondrial cytochrome b gene and variable control region indicate significant substructure within the Cape rock elephant-shrew across its distribution. Our data unequivocally showed the presence of a northern Namaqua and central Fynbos clade with four evolutionary lineages identified within the latter. The geographical delimitation of the northern and central clades corresponds closely with patterns reported for other rock-dwelling vertebrate species, indicating a shared biogeographical history for these taxa in South Africa. A coalescent method revealed the effects of ancestral polymorphism in shaping the Namaqua and Fynbos populations since their divergence ~1.7 million years ago. Furthermore, our analyses uncovered a distinct Karoo lineage(s) that does not correspond to any of the previously described and/or currently recognized species, and we therefore argue for the possible recognition of a new sister taxon to E. edwardii. The taxonomic affinities of this clade were examined by sequencing corresponding regions from the type specimens of species described in the past, but which presently are synonimized within E. edwardii. Our results reveal the morphological misidentification of one of these types, accentuating the problems of field identification.     

Phylogeography and Pleistocene refugia of the Little Owl Athene noctua inferred from mtDNA sequence data

Pleistocene glaciations greatly affected the distribution of genetic diversity in animal populations. The Little Owl is widely distributed in temperate regions and could have survived the last glaciations in southern refugia. To describe the phylogeographical structure of European populations, we sequenced the mitochondrial cytochrome c oxidase I (COI) and control region (CR1) in 326 individuals sampled from 22 locations. Phylogenetic analyses of COI identified two deeply divergent clades: a western haplogroup distributed in western and northwestern Europe, and an eastern haplogroup distributed in southeastern Europe. Faster evolving CR1 sequences supported the divergence between these two main clades, and identified three subgroups within the eastern clade: Balkan, southern Italian and Sardinian. Divergence times estimated from COI with fossil calibrations indicate that the western and eastern haplogroups split 2.01–1.71 Mya. Slightly different times for splits were found using the standard 2% rate and 7.3% mtDNA neutral substitution rate. CR1 sequences dated the origin of endemic Sardinian haplotypes at 1.04–0.26 Mya and the split between southern Italian and Balkan haplogroups at 0.72–0.21 Mya, coincident with the onset of two Pleistocene glaciations. Admixture of mtDNA haplotypes was detected in northern Italy and in central Europe. These findings support a model of southern Mediterranean and Balkan refugia, with postglacial expansion and secondary contacts for Little Owl populations. Central and northern Europe was predominantly recolonized by Little Owls from Iberia, whereas expansion out of the Balkans was more limited. Northward expansion of the Italian haplogroup was probably prevented by the Alps, and the Sardinian haplotypes remained confined to the island. Results showed a clear genetic pattern differentiating putative subspecies. Genetic distances between haplogroups were comparable with those recorded between different avian species.     

Trans-Arctic dispersals and the evolution of a circumpolar marine fish species complex, the capelin (Mallotus villosus)

Trans-Arctic dispersals and population and range expansions during the Pleistocene enhanced opportunities for evolutionary diversification and contributed to the process of speciation within the capelin, a northern marine-fish complex exhibiting a circumpolar distribution. Capelin is composed of four highly divergent and geographically discrete mitochondrial DNA (mtDNA) clades (609 bp; cytochrome b). Two clades occur in the North Atlantic, one associated with Canadian Atlantic waters, including Hudson Bay, and the second distributed from West Greenland to the Barents Sea. Two additional clades occur in the Arctic and northeast Pacific Oceans, representing the most recent divergence within the capelin phylogenetic tree. Judged from mtDNA diversity, capelin populations comprising all clades experienced at least one demographic and spatial reduction-expansion episode during recent Pleistocene glaciations that imprinted their molecular architecture. The large contemporary populations in the northeast Pacific and Arctic Oceans exhibited significant genetic structure whereas no such structure was detected in the equally extensive North Atlantic clades. All clades are characterized by one or two prevalent mtDNA haplotypes distributed over the entire range of the clade. Assuming a Pacific ancestor for capelin, we infer that capelin dispersed on two separate occasions to the North Atlantic. A more recent event resulted in the isolation of eastern Pacific and Arctic clades, with the Arctic clade positioned for a potential third Atlantic invasion, as revealed by the presence of this clade in the Labrador Sea. The Labrador Sea is a potential contact zone for three of the four capelin clades.     

Molecular evidence that Langeronia macrocirra and Langeronia cf. parva (Trematoda: Pleurogenidae) parasites of anurans from Mexico are conspecific

The genus Langeronia parasitizing the intestine of several species of anurans is distributed from North to Central America. We identified Langeronia macrocirra and Langeronia cf. parva from the same host and localities, and present here new data not applicable about their tegumental surface by scanning electron microscopy. We compared sequences of the rDNA ITS2 region and mtDNA cox1 gene for the two morphotypes. ITS2 exhibited a high degree of conservation. Phylogenetic reconstruction using cox1 revealed three clades (I, II, and III), which did not correspond to a previous identification or host. Little divergence was found within clades: sequences were identical in clade I, whereas clade II had 0.27% and clade III had 1.08%. Inter-clade divergence reached 8.69% (I vs. III). This pattern of genetic divergence indicated that both taxa probably belong to the same species, so we posit that the morphological changes could be correlated with development. Increasing sample size and geographical coverage will contribute to the taxonomy of the genus based on morphological and molecular evidence, and will open tracks toward the use of DNA barcodes to the genus in Mexico.     

A phylogenetic delimitation of the "Sphagnum subsecundum complex" (Sphagnaceae, Bryophyta)

A seemingly obvious but sometimes overlooked premise of any evolutionary analysis is delineating the group of taxa under study. This is especially problematic in some bryophyte groups because of morphological simplicity and convergence. This research applies information from nucleotide sequences for eight plastid and nuclear loci to delineate a group of northern hemisphere peat moss species, the so-called Sphagnum subsecundum complex, which includes species known to be gametophytically haploid or diploid (i.e., sporophytically diploid-tetraploid). Despite the fact that S. subsecundum and several species in the complex have been attributed disjunct ranges that include all major continents, phylogenetic analyses suggest that the group is actually restricted to Europe and eastern North America. Plants from western North America, from California to Alaska, which are morphologically similar to species of the S. subsecundum complex in eastern N. America and Europe, actually belong to a different deep clade within Sphagnum section Subsecunda. One species often considered part of the S. subsecundum complex, S. contortum, likely has a reticulate history involving species in the two deepest clades within section Subsecunda. Nucleotide sequences have a strong geographic structure across the section Subsecunda, but shallow tip clades suggest repeated long-distance dispersal in the section as well.     

Phylogeography of a mountain lizard species: an ancient fragmentation process mediated by riverine barriers in the Liolaemus monticola complex (Sauria: Liolaemidae)

Liolaemus monticola is a mountain lizard species, with a widespread distribution from central Chile that displays several highly polymorphic chromosomal races. Our study determined the phylogeographic structuring and relationships among three chromosomal races of L. monticola in Chile. Mitochondrial DNA (mtDNA) sequences of the cytochrome b gene were examined using the following phylogenetic methods: maximum parsimony, maximum likelihood, Bayesian inference and nested clade phylogeographic analyses (NCPAs). These methods revealed two major monophyletic clades (north and south) in the L. monticola species, with non-overlapping geographical locations separated by the Maipo and Yeso rivers (except one hybrid, from a zone of secondary contact). The NCPA showed that a past fragmentation process likely resulted in the separation of the two clades. The southern clade includes all samples of the 'Southern, 2 n  = 34' race; the northern clade is comprised of all remaining derived chromosomal races: the 'Northern, 2 n  = 38–40 and the Multiple Fission, 2 n  = 42–44' races. Our results support the hypothesis of a geographical and genetic split resulting from allopatric processes caused by riparian barriers acting over a long time period. The inferred biogeographical scenario shows that populations have moved from the south to the north using the Andean mountains as the primary corridor for dispersal.     

Genetic differentiation of an endangered capercaillie (<Emphasis Type="Italic">Tetrao urogallus</Emphasis>) population at the Southern edge of the species range

The low-latitude limits of species ranges are thought to be particularly important as long-term stores of genetic diversity and hot spots for speciation. The Iberian Peninsula, one of the main glacial refugia in Europe, houses the southern distribution limits of a number of boreal species. The capercaillie is one such species with a range extending northwards to cover most of Europe from Iberia to Scandinavia and East to Siberia. The Cantabrian Range, in North Spain, constitutes the contemporary south-western distribution limit of the species. In contrast to all other populations, which live in pure or mixed coniferous forests, the Cantabrian population is unique in inhabiting pure deciduous forests. We have assessed the existence of genetic differentiation between this and other European populations using microsatellite and mitochondrial DNA (mtDNA) extracted from capercaillie feathers. Samples were collected between 2001 and 2004 across most of the current distribution of the Cantabrian population. Mitochondrial DNA analysis showed that the Cantabrian birds form a distinct clade with respect to all the other European populations analysed, including the Alps, Black Forest, Scandinavia and Russia, which are all members of a discrete clade. Microsatellite DNA from Cantabrian birds reveals the lowest genetic variation within the species in Europe. The existence of birds from both mtDNA clades in the Pyrenees and evidence from microsatellite frequencies for two different groups, points to the existence of a Pyrenean contact zone between European and Cantabrian type birds. The ecological and genetic differences of the Cantabrian capercaillies qualify them as an Evolutionarily Significant Unit and support the idea of the importance of the rear edge for speciation. Implications for capercaillie taxonomy and conservation are discussed.     

Exploring factors shaping population genetic structure of the freshwater fish Sinocyclocheilus grahami (Teleostei, Cyprinidae)

Phylogeographical analyses on Sinocyclocheilus grahami samples from seven localities within the Lake Dianchi Basin in China were conducted to explore the main factors shaping population structure within this species. Phylogenetic and network analyses revealed two major clades in 24 mtDNA haplotypes. One clade included three haplotypes exclusively from samples of the lower basin and another clade encompassed other haplotypes from samples of the upper basin. The estimated divergence time between the two clades predated the river capture event connecting the lower and upper lake basin and thus supported geographical isolation as the main factor shaping genetic divergence between these two clades. Furthermore, analysis of molecular variance and pair-wise Φ_ST distances revealed significant genetic differentiation within the upper basin. Mantel tests clearly supported patterns of differentiation arose purely as a result of isolation by distance. These results further highlight the importance of geographical isolation in shaping differentiation within this species.     

Following the rivers: historical reconstruction of California voles Microtus californicus (Rodentia: Cricetidae) in the deserts of eastern California

The California vole, Microtus californicus, restricted to habitat patches where water is available nearly year‐round, is a remnant of the mesic history of the southern Great Basin and Mojave deserts of eastern California. The history of voles in this region is a model for species‐edge population dynamics through periods of climatic change. We sampled voles from the eastern deserts of California and examined variation in the mitochondrial cytb gene, three nuclear intron regions, and across 12 nuclear microsatellite markers. Samples are allocated to two mitochondrial clades: one associated with southern California and the other with central and northern California. The limited mtDNA structure largely recovers the geographical distribution, replicated by both nuclear introns and microsatellites. The most remote population, Microtus californicus scirpensis at Tecopa near Death Valley, was the most distinct. This population shares microsatellite alleles with both mtDNA clades, and both its northern clade nuclear introns and southern clade mtDNA sequences support a hybrid origin for this endangered population. The overall patterns support two major invasions into the desert through an ancient system of riparian corridors along streams and lake margins during the latter part of the Pleistocene followed by local in situ divergence subsequent to late Pleistocene and Holocene drying events. Changes in current water resource use could easily remove California voles from parts of the desert landscape.     

Phylogeography of white-spotted Charr (Salvelinus leucomaenis) inferred from mitochondrial DNA sequences

The white-spotted charr (Salvelinus leucomaenis) is a coldwater-adapted fish distributed in far-eastern Asia. To assess phylogeographic patterns of this species over most of its range in the Japanese archipelago and Sakhalin Island, Russia, we examined nucleotide sequences of the mitochondrial DNA (mtDNA) cytochrome b region (557 bp) in 141 individuals from 50 populations. A total of 33 (5.5%) nucleotide positions were polymorphic and defined 29 haplotypes. Phylogenetic analysis assigned the observed haplotypes to four main clades, which were characterized by the idiosyncrasies and discontinuity of geographic distributions. The nested clade analyses revealed that the geographical distribution patterns of some haplotypes and clades were explained by historical event such as past fragmentation. Although substantial genetic differentiation was found among the four main clades, their geographic distributions overlapped extensively in several regions. Since white-spotted charr can potentially use both freshwater and marine environments, coexistence among different lineages can be attributed to secondary contact through range expansion by migratory individuals during multiple glacial periods after interglacial isolation. Finally, our data demonstrate that the current subspecies designation does not reflect the phylogeography of this species based on mtDNA analysis. Hierarchical analysis (AMOVA) also showed that genetic variation was far more pronounced within subspecies than among subspecies (i.e., among discrete regions). These results suggest that each population, rather than each subspecies, must be treated as an evolutionarily significant unit.     

Selective neutrality of mitochondrial ND2 sequences, phylogeography and species limits in Sitta europaea

Variation and geographic differentiation in mitochondrial DNA (mtDNA) was studied in the widespread and phenotypically variable Eurasian nuthatch (Sitta europaea). To assess whether sequences were evolving in a selectively neutral fashion, we used McDonald–Kreitman [Nature 351 (1991) 652] tests and a tree-based method, which suggested that although ND2 sequences are affected by natural selection against slightly deleterious alleles, the effects do not compromise phylogeographic inferences. Three phylogenetic species-level clades of nuthatches were discovered, corresponding to the Caucasus, southern Europe, and northern Europe plus Asia. Unimodal mismatch distributions within each clade suggest that populations have undergone recent growth. A westward range expansion was inferred from the geographic pattern in nucleotide diversity. Although samples were insufficient, it is possible that nuthatches in England and Japan are recently differentiated. Two specimens of the subspecies S. e. arctica formed a sister group to all other S. europaea, differing by ca. 10% uncorrected sequence divergence, pointing the need for additional study of this phenotypically distinct taxon. As with other species, mtDNA data support major phenotypic distinctions, but not subspecies.