A test of the European Pleistocene refugial paradigm, using a Western Palaearctic endemic bird species.

Hewitt's paradigm for effects of Pleistocene glaciations on European populations assumes their isolation in peninsular refugia during glacial maxima, followed by re-colonization of broader Europe during interstadials. This paradigm is well supported by studies of poorly dispersing taxa, but highly dispersive birds have not been included. To test this paradigm, we use the dunnock (Prunella modularis), a Western Palaearctic endemic whose range includes all major European refugia. MtDNA gene tree, multilocus species tree and species delimitation analyses indicate the presence of three distinct lineages: one in the Iberian refugium, one in the Caucasus refugium, and one comprising the Italian and Balkan refugia and broader Europe. Our gene flow analysis suggests isolation of both the Iberian and Caucasus lineages but extensive exchange between Italy, the Balkans and broader Europe. Demographic stability could not be rejected for any refugial population, except the very recent expansion in the Caucasus. By contrast, northern European populations may have experienced two expansion periods. Iberia and Caucasus had much smaller historical populations than other populations. Although our results support the paradigm, in general, they also suggest that in highly dispersive taxa, isolation of neighbouring refugia was incomplete, resulting in large super-refugial populations.

Multi-locus reassessment of a striking discord between mtDNA gene trees and taxonomy across two congeneric species complexes.

Resolving relationships among members of the yellow and citrine wagtail species complexes is among the greatest challenges in avian systematics due to arguably the most dramatic disagreements between traditional taxonomy and mtDNA phylogeny. Each species complex is divided into three geographically cohesive mtDNA clades. Each clade from one species complex has a sister from the other complex. Furthermore, one cross-complex pair is more distantly related to the remaining two pairs than are several other wagtail species. To test mtDNA gene tree topology, we sequenced the mtDNA ND2 gene and 11 nuclear introns for seven wagtail species. Our mtDNA gene tree reconstruction supported the results of previous studies, thereby confirming the disagreement between mtDNA phylogeny and taxonomy. However, our multi-locus species tree which used mtDNA clades as "taxa" was consistent with traditional taxonomy regardless of whether mtDNA was included in the analysis or not. Our multi-locus data suggest that despite the presence of strongly supported, geographically structured mtDNA variation, the mtDNA gene tree misrepresents the evolutionary history of the yellow and citrine wagtail complexes. This mito-nuclear discord results from mtDNA representing the biogeographic, but not evolutionary history of these recently radiated Palearctic wagtails.

Geographic patterns of mtDNA and Z-linked sequence variation in the Common Chiffchaff and the 'chiffchaff complex'.

The Common Chiffchaff Phylloscopus collybita is an abundant, polytypic Palearctic bird. Validity of some of its subspecies is controversial and birds from some parts of the species range remain unclassified taxonomically. The relationships among populations from different geographic areas have not been sufficiently explored with molecular data. In this study we analyzed the relationships among the four species in the 'chiffchaff complex' (Common Chiffchaff, Iberian Chiffchaff P. ibericus, Canary Islands Chiffchaff P. canariensis and Mountain Chiffchaff P. sindianus), and the patterns of intraspecific geographic variation in the mtDNA ND2 gene and intron 9 of the Z-linked aconitase gene (ACO1I9) across the Common Chiffchaff range, including a recently discovered population breeding on Mt. Hermon (Anti-Lebanon mountains). Our data supported the monophyly of the chiffchaff complex and its current systematics at the species level. Within the Common Chiffchaff, the Siberian race P. c. tristis was the most differentiated subspecies and may represent a separate or incipient species. Other Common Chiffchaff subspecies also were differentiated in their mtDNA, however, lineages of neighboring subspecies formed wide zones of introgression. The Mt. Hermon population was of mixed genetic origin but contained some birds with novel unique lineage that could not be assigned to known subspecies. All Common Chiffchaff lineages diverged at the end of the Ionian stage of Pleistocene. Lineage sorting of ACO1I9 alleles was not as complete as that of mtDNA. Chiffchaff species were mostly distinct at ACO1I9, except the Common and Canary Islands Chiffchaffs that shared multiple alleles. An AMOVA identified geographic structure in Common Chiffchaff ACO1I9 variation that was broadly consistent with that of mtDNA ND2 gene. The genetic and other data suggest the chiffchaff complex to be a group of evolutionarily young taxa that represent a paradigm of 'species evolution in action' from intergrading subspecies through to apparently complete biological speciation.

Effects of asymmetric nuclear introgression, introgressive mitochondrial sweep, and purifying selection on phylogenetic reconstruction and divergence estimates in the Pacific clade of Locustella warblers.

When isolated but reproductively compatible populations expand geographically and meet, simulations predict asymmetric introgression of neutral loci from a local to invading taxon. Genetic introgression may affect phylogenetic reconstruction by obscuring topology and divergence estimates. We combined phylogenetic analysis of sequences from one mtDNA and 12 nuDNA loci with analysis of gene flow among 5 species of Pacific Locustella warblers to test for presence of genetic introgression and its effects on tree topology and divergence estimates. Our data showed that nuDNA introgression was substantial and asymmetrical among all members of superspecies groups whereas mtDNA showed no introgression except a single species pair where the invader's mtDNA was swept by mtDNA of the local species. This introgressive sweep of mtDNA had the opposite direction of the nuDNA introgression and resulted in the paraphyly of the local species' mtDNA haplotypes with respect to those of the invader. Тhe multilocus nuDNA species tree resolved all inter- and intraspecific relationships despite substantial introgression. However, the node ages on the species tree may be underestimated as suggested by the differences in node age estimates based on non-introgressing mtDNA and introgressing nuDNA. In turn, the introgressive sweep and strong purifying selection appear to elongate internal branches in the mtDNA gene tree.

Complex biogeographic history of a Holarctic passerine.

Our analysis of the ND2 sequences revealed six clades within winter wrens (Troglodytes troglodytes). These clades corresponded to six geographical regions: western Nearctic, eastern Nearctic, eastern Asia, Nepal, Caucasus and Europe, and differed by 3-8.8% of sequence divergence. Differences among regions explained 96% of the sequence variation in winter wren. Differences among individuals within localities explained 3% of the sequence variation, and differences among localities within regions explained 1%. Grouping sequences into subspecies instead of localities did not change these proportions. Proliferation of the six clades coincided with Early and Middle Pleistocene glaciations. The distribution of winter wren clades can be explained by a series of five consecutive vicariant events. Western Nearctic wrens diverged from the Holarctic ancestor 1.6 Myr before the present time (MYBP). Eastern Nearctic and Palaearctic wrens diverged 1 MYBP. Eastern and western Palaearctic birds diverged 0.83 MYBP. Nepalese and east Asian wrens diverged 0.67 MYBP, and Caucasian birds diverged from European wrens 0.54 MYBP. The winter wren has a much greater degree of inter- and intracontinental differentiation than the three other Holarctic birds studied to date--dunlin (Calidris alpina), common raven (Corvus corax) and three-toed woodpecker (Picoides trydactylus)--and represents an example of cryptic speciation that has been overlooked.

Limited phylogeographic signal in sex-linked and autosomal loci despite geographically, ecologically, and phenotypically concordant structure of mtDNA variation in the Holarctic avian genus Eremophila.

Phylogeographic studies of Holarctic birds are challenging because they involve vast geographic scale, complex glacial history, extensive phenotypic variation, and heterogeneous taxonomic treatment across countries, all of which require large sample sizes. Knowledge about the quality of phylogeographic information provided by different loci is crucial for study design. We use sequences of one mtDNA gene, one sex-linked intron, and one autosomal intron to elucidate large scale phylogeographic patterns in the Holarctic lark genus Eremophila. The mtDNA ND2 gene identified six geographically, ecologically, and phenotypically concordant clades in the Palearctic that diverged in the Early-Middle Pleistocene and suggested paraphyly of the horned lark (E. alpestris) with respect to the Temminck's lark (E. bilopha). In the Nearctic, ND2 identified five subclades which diverged in the Late Pleistocene. They overlapped geographically and were not concordant phenotypically or ecologically. Nuclear alleles provided little information on geographic structuring of genetic variation in horned larks beyond supporting the monophyly of Eremophila and paraphyly of the horned lark. Multilocus species trees based on two nuclear or all three loci provided poor support for haplogroups identified by mtDNA. The node ages calculated using mtDNA were consistent with the available paleontological data, whereas individual nuclear loci and multilocus species trees appeared to underestimate node ages. We argue that mtDNA is capable of discovering independent evolutionary units within avian taxa and can provide a reasonable phylogeographic hypothesis when geographic scale, geologic history, and phenotypic variation in the study system are too complex for proposing reasonable a priori hypotheses required for multilocus methods. Finally, we suggest splitting the currently recognized horned lark into five Palearctic and one Nearctic species.

Geographic mode of speciation in a mountain specialist Avian family endemic to the Palearctic.

Mountains host greater avian diversity than lowlands at the same latitude due to their greater diversity of habitats stratified along an elevation gradient. Here we test whether this greater ecological heterogeneity promotes sympatric speciation. We selected accentors (Prunellidae), an avian family associated with mountains of the Palearctic, as a model system. Accentors differ in their habitat/elevation preferences and south-central Siberia and Himalayan regions each host 6 of the 13 species in the family. We used sequences of the mtDNA ND2 gene and the intron 9 of the Z chromosome specific ACO1 gene to reconstruct a complete species-level phylogeny of Prunellidae. The tree based on joint analysis of both loci was used to reconstruct the family's biogeographic history and to date the diversification events. We also analyzed the relationship between the node age and sympatry, to determine the geographic mode of speciation in Prunellidae. Our data suggest a Miocene origin of Prunellidae in the Himalayan region. The major division between alpine species (subgenus Laiscopus) and species associated with shrubs (subgenus Prunella) and initial diversification events within the latter happened within the Himalayan region in the Miocene and Pliocene. Accentors colonized other parts of the Palearctic during the Pliocene-Pleistocene transition. This spread across the Palearctic resulted in rapid diversification of accentors. With only a single exception dating to 0.91 Ma, lineages younger than 1.5 Ma are allopatric. In contrast, sympatry values for older nodes are >0. There was no relationship between node age and range symmetry. Allopatric speciation (not to include peripatric) is the predominant geographic mode of speciation in Prunellidae despite the favorable conditions for ecological diversification in the mountains and range overlaps among species.

Recent evolutionary history of the bluethroat (Luscinia svecica) across Eurasia.

We analysed mitochondrial DNA (mtDNA) sequences from 154 bluethroats (Luscinia svecica) sampled at 21 sites throughout much of their Eurasian range. A previously reported, single base-pair mtDNA difference between L. s. svecica and L. s. namnetum was inconsistent upon expanded geographical sampling. A significant FST value (0.29) and an isolation-by-distance effect show the existence of geographical differentiation. Phylogenetic analysis of haplotypes revealed northern and southern groups, although lineage sorting is incomplete. There was no geographical structure to the haplotype tree within groups, and currently recognized subspecies were not supported. A minimum evolution tree based on pairwise mtDNA genetic distances among average samples showed the same two broadly distributed northern and southern groups. These groups abut in the centre of the latitudinal range, and were possibly isolated by forest that developed and spread westward over the last 15 000 years. Pairwise FST values averaged 0.16 in the southern group, 0.04 in the northern group, and 0.42 between groups. Mismatch distributions suggested population growth in each group, with that in the south being more recent. In the northern group, the geographical pattern in tau suggested northward and eastward expansion. Analysis of nucleotide diversity suggested westward expansion in the southern group. The northern group had higher nucleotide diversity than the southern group, consistent with a larger current population size in the north. Given the significant FST, incompletely sorted haplotype tree, and broadly patterned minimum evolution tree, L. svevica appears to represent a species at an intermediate stage of differentiation between panmixia and reciprocal monophyly.