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.

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.

Extensive phenotypic diversification coexists with little genetic divergence and a lack of population structure in the White Wagtail subspecies complex (Motacilla alba).

Geographically clustered phenotypes often demonstrate consistent patterns in molecular markers, particularly mitochondrial DNA (mtDNA) traditionally used in phylogeographic studies. However, distinct evolutionary trajectories among traits and markers can lead to their discordance. First, geographic structure in phenotypic traits and nuclear molecular markers can be co-aligned but inconsistent with mtDNA (mito-nuclear discordance). Alternatively, phenotypic variation can have little to do with patterns in neither mtDNA nor nuclear markers. Disentangling between these distinct patterns can provide insight into the role of selection, demography and gene flow in population divergence. Here, we examined a previously reported case of strong inconsistency between geographic structure in mtDNA and plumage traits in a widespread polytypic bird species, the White Wagtail (Motacilla alba). We tested whether this pattern is due to mito-nuclear discordance or discrepancy between morphological evolution and both nuclear and mtDNA markers. We analysed population differentiation and structure across six out of nine commonly recognized subspecies using 17 microsatellite loci and a combination of microsatellites and plumage indices in a comprehensively sampled region of a contact between two subspecies. We did not find support for the mito-nuclear discordance hypothesis: nuclear markers indicated a subtle signal of genetic clustering only partially consistent with plumage groups, similar to previous findings that relied on mtDNA. We discuss evolutionary factors that could have shaped the intricate patterns of phenotypic diversification in the White wagtail and the role that repeated selection on plumage 'hotspots' and hybridization may have played.

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.

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.