A comparative phylogenomic analysis of birds reveals heterogeneous differentiation processes among Neotropical savannas.

The main objective of this study was to evaluate biogeographical hypotheses of diversification and connection between isolated savannas north (Amazonian savannas) and south (Cerrado core) of the Amazon River. To achieve this, we used genomic markers (genotyping-by-sequencing) to evaluate the genetic structure, population phylogenetic relationships and historical range shifts of four Neotropical passerines with peri-Atlantic distributions: the narrow-billed woodcreeper (Lepidocolaptes angustirostris), the plain-crested elaenia (Elaenia cristata), the grassland sparrow (Ammodramus humeralis) and the white-banded tanager (Neothraupis fasciata). Population genetic analyses indicated that landscape (e.g., geographic distance, landscape resistance and percentage of tree cover) and climate metrics explained divergence among populations in most species, but without indicating a differential role between current and historical factors. Our results did not fully support the hypothesis that isolated populations in Amazonian savannas have been recently derived from the Cerrado core domain. Intraspecific phylogenies and gene flow analyses supported multiple routes of connection between the Cerrado and Amazonian savannas, rejecting the hypothesis that the Atlantic corridor explains the peri-Atlantic distribution. Our results reveal that the biogeographical history of the region is complex and cannot be explained by simple vicariant models.

Comparative phylogenomic patterns in the Baja California avifauna, their conservation implications, and the stages in lineage divergence.

Comparative phylogeography explores the historical congruence of co-distributed species to understand the factors that led to their current genetic and phenotypic structures. Even species that span the same biogeographic barrier can exhibit different phylogeographic structures owing to differences in effective population sizes, genetic marker bias, and dispersal abilities. The Baja California peninsula and adjacent desert regions include several biogeographic barriers, including the Vizcaíno Desert and Sierra de la Laguna (Cape District), that have left phylogeographic patterns in some but not all species. We used genome-wide SNP data to test the hypothesis that the diverse phylogeographic patterns inferred from prior studies were supported. We found that mitochondrial DNA, single nuclear gene, and genome-wide SNP data show that the cactus wren and LeConte's thrasher have a concordant historical division at or near the Vizcaíno Desert in north-central Baja California, the Gila woodpecker is at an intermediate stage of divergence, and the California gnatcatcher lacks phylogeographic structure. None of these four species are classified taxonomically in a way that captures their evolutionary history with the exception of the LeConte's thrasher. We also analyzed mtDNA data on samples of nine other species that span the Vizcaíno Desert, with four showing no apparent division, and six additional species from the Sierra de la Laguna, all but one of which are differentiated. Reasons for contrasting phylogeographic patterns among these species should be explored further with genomic data to test the extent of concordant phylogeographic patterns. The evolutionary division at the Vizcaíno desert is well known in other vertebrate species, and our study further corroborates the extent, profound effect, and importance of this biogeographic boundary. The areas north and south of the Vizcaíno Desert, which contains considerable diversity, should be recognized as historically significant areas for conservation.

Considering the use of the terms strain and adaptation in prion research.

Evolutionary biologists and disease biologists use the terms strain and adaptation in Chronic Wasting Disease (CWD) research in different ways. In evolutionary biology, a strain is a nascent genetic lineage that can be described by a genealogy, and a phylogenetic nomenclature constructed to reflect that genealogy. Prion strains are described as showing distinct host range, clinical presentation, disease progression, and neuropathological and PrP biochemical profiles, and lack information that would permit phylogenetic reconstruction of their history. Prion strains are alternative protein conformations, sometimes derived from the same genotype. I suggest referring to prion strains as ecotypes, because the variant phenotypic conformations ("strains") are a function of the interaction between PRNP amino acid genotype and the host environment. In the case of CWD, a prion ecotype in white-tailed deer would be described by its genotype and the host in which it occurs, such as the H95 + ecotype. However, an evolutionary nomenclature is difficult because not all individuals with the same PRNP genotype show signs of CWD, therefore creating a nomenclature reflecting and one-to-one relationship between PRNP genealogy and CWD presence is difficult. Furthermore, very little information exists on the phylogenetic distribution of CWD ecotypes in wild deer populations. Adaptation has a clear meaning in evolutionary biology, the differential survival and reproduction of individual genotypes. If a new prion ecotype arises in a particular host and kills more hosts or kills at an earlier age, it is the antithesis of the evolutionary definition of adaptation. However, prion strains might be transmitted across generations epigenetically, but whether this represents adaptation depends on the fitness consequences of the strain. Protein phenotypes of PRNP that cause transmissible spongiform encephalopathies (TSEs), and CWD, are maladaptive and would not be propagated genetically or epigenetically via a process consistent with an evolutionary view of adaptation. I suggest terming the process of prion strain origination "phenotypic transformation", and only adaptation if evidence shows they are not maladaptive and persist over evolutionary time periods (e.g., thousands of generations) and across distinct species boundaries (via inheritance). Thus, prion biologists use strain and adaptation, historically evolutionary terms, in quite different ways.

Geographic variation in the PRNP gene and its promoter, and their relationship to chronic wasting disease in North American deer.

PRNP genotypes, number of octarepeats (PHGGGWGQ) and indels in the PRNP promoter can influence the progression of prion disease in mammals. We found no relationship between presence of promoter indels in white-tailed deer and mule deer from Nebraska and CWD presence. White-tailed deer with the 95 H allele and G20D mule deer were more likely to be CWD-free, but unlike other studies white-tailed deer with the 96S allele(s) were equally likely to be CWD-free. We provide the first information on PRNP genotypes and indels in the promoter for Key deer (all homozygous 96SS) and Coues deer (lacked 95 H and 96S alleles, but possessed a uniquely high frequency of 103 T). All deer surveyed were homozygous for three tandem octarepeats.

Genetic and evolutionary considerations of the Chronic Wasting Disease - Human species barrier.

Transmissible spongiform encephalopathies can jump species barriers. In relatively few cases is the possible route of transmission thought to be known, mostly involving humans, cattle and sheep. It is thought that sheep might be the cause of Bovine Spongiform Encephalopathy (BSE) and Chronic Wasting Disease (CWD) in cervids, and that humans might have gotten prion disease (e.g., vCJD) from eating meat from BSE+ cows. A looming societal question is whether humans will acquire a prion disease from ingesting prions from CWD+ deer. On an evolutionary tree of the PRNP gene in mammals, deer, sheep and cow are relatively closely related, whereas these three species are relatively distant from humans. If a prion disease jumped the species barrier from cow to humans, the phylogenetic gap from deer to humans is no greater, and sheer evolutionary distance alone cannot explain a CWD species barrier in humans. Aspects of the PRNP gene were compared among these species to search for genetic differences that might influence the permeability of the species barrier. Human prion disease has been associated with having more than four copies of the octarepeat unit (PHGGGWG), whereas deer, sheep and cow all have three copies. Two amino acid positions in the metal-binding region (96 and 97) have been implicated in species barriers (Breydo and Uversky, 2011), whereas no variation was detected in white-tailed deer and mule deer with and without CWD, or in black-tailed deer, Key deer or Coues deer. Four out of 10 differences between deer and human in the ß2-α2 loop might preclude CWD prions from converting human PrPC to PrPSc because of disruption of a steric zipper. The reasons for a CWD species barrier between deer and humans, if there is one, is still unresolved.

Niche modeling reveals life history shifts in birds at La Brea over the last twenty millennia.

A species presence at a particular site can change over time, resulting in temporally dynamic species pools. Ecological niche models provide estimates of species presence at different time intervals. The avifauna of La Brea includes approximately 120 species dating to approximately 15,000 years ago. Niche models predicted presence at the Last Glacial Maximum for over 90% of 89 landbird species. This confirms that niche modeling produces sensible range estimates at the Last Glacial Maximum. For 97 currently local species that are as yet undocumented at La Brea over 90% were predicted to occur; absence is due to insufficient study, lack of the ecological niche, transient occurrence or a behavioral ability to avoid entrapment. Our 366 niche models provide a prospective checklist of the landbird fauna of La Brea. The models indicate fluidity in life history strategies and a higher proportion of resident birds at the LGM (88% to 60%). We evaluated a subset of 103 species in breeding and winter periods using two climate models (MIROC-ESM, CCSM4) with a variety of differing parameters, finding differences in 5% of the niche models. Niche breadths in bark-foraging birds changed little between the present and LGM, suggesting that greater species diversity at the LGM was due to greater niche availability rather than contractions of niche breadths (i.e., niche partitioning).

Species Limits and Phylogenomic Relationships of Darwin's Finches Remain Unresolved: Potential Consequences of a Volatile Ecological Setting.

Island biotas have become paradigms for illustrating many evolutionary processes. The fauna of the Galapagos Islands includes several taxa that have been focal points for evolutionary studies. Perhaps their most famous inhabitants, Darwin's finches, represent a go-to icon when thinking about how species originate and adapt to the environment. However, unlike other adaptive radiations, past morphological and molecular studies of Darwin's finches have yielded inconsistent hypotheses of species limits and phylogenetic relationships. Expecting that idiosyncrasies of prior data and analytic methods explained different proposed classifications, we were surprised to observe that three new phylogenetic hypotheses derived mostly from the same genomics data were topologically inconsistent. We found that the differences between some of these genomics trees were as great as one would expect between two random trees with the same number of taxa. Thus, the phylogeny of Darwin's finches remains unresolved, as it has for more than a century. A component of phylogenetic uncertainty comes from unclear species limits, under any species concept, in the ground finches (Geospiza) and tree finches (Camarhynchus). We suggest that past authors should have tested the species limits of Lack, rather than uncritically accepting them. In fact, the impressive amount of genomics data do not provide unambiguous hypotheses of the number of species of Geospiza or Camarhynchus, although they imply greater species diversity than Lack's taxonomy. We suggest that insufficient sampling of species populations across islands (35.6% for morphometrics and 20.4% for genomics) prevents accurate diagnoses of species limits. However, it is unknown whether samples from a greater number of islands might result in bridging differences between species, or reveal many new ones. We conclude that attempts to interpret patterns of variation among the finches under standard evolutionary paradigms have obscured some major messages, most specifically the ongoing reciprocal interactions between geographic isolation and lineage divergence, and dispersal and gene flow caused by the volatile ecological conditions in the islands. Although the finches provide textbook examples of natural selection, better understanding of species limits and a robust phylogenetic hypothesis are required to corroborate past hypotheses of speciation and adaptive radiation in the finches of the Galapagos.

Glaciation as a migratory switch.

Migratory behavior in birds is evolutionarily plastic, but it is unclear how this behavior responded during glacial cycles. One view is that at glacial maxima, species simply shifted their breeding ranges south of glacial ice and remained migratory. To test this hypothesis, we constructed ecological niche models for breeding and wintering ranges of 56 species, finding that 70% of currently long-distance North American migrant species likely lacked suitable breeding habitat in North America at the Last Glacial Maximum (LGM), and we hypothesized that they reverted to the ancestral state of being tropical sedentary residents. A smaller percentage of short-distance migrants (27%) experienced altered migratory behavior at the LGM, although the ranges of all species were shifted southward, as traditionally envisioned. We suggest that many species oscillate between sedentary and migratory strategies with each glacial cycle acting as an adaptive switch. Thus, range shifts occur more frequently than speciation events and are inadequately depicted by phylogenetic reconstructions. We suggest that reconstructing the evolutionary history of traits, such as migratory behavior, is best served by using ranges at glacial maxima. A phylogeny of warblers strongly predicted LGM, but not present distributions, and suggested that migration was re-expressed from three tropical centers of warbler diversity. Understanding of evolutionary history will be improved when reconstructions use distributions relevant to the time of character transitions.

The roles of ecology, behaviour and effective population size in the evolution of a community.

Organismal traits such as ecological specialization and migratory behaviour may affect colonization potential, population persistence and degree of isolation, factors that determine the composition and genetic structure of communities. However, studies focusing on community assembly rarely consider these factors jointly. We sequenced 16 nuclear genes and one mitochondrial gene from Caucasian and European populations of 30 forest-dwelling avian species that represent diverse ecological (specialist-generalist) and behavioural (migratory-resident) backgrounds. We tested the effects of organismal traits on population divergence and community assembly in the Caucasus forest, a continental mountain island setting. We found that (i) there is no concordance in divergence times between the Caucasus forest bird populations and their European counterparts, (ii) habitat specialists tend to be more divergent than generalists and (iii) residents tend to be more divergent than migrants. Thus, specialists and residents contribute to the high level of endemism of Caucasus forest avifauna more than do generalists and migrants. Patterns of genetic differentiation are better explained by differences in effective population sizes, an often overlooked factor in comparative studies of phylogeography and speciation, than by divergence times or levels of gene flow. Our results suggest that the Caucasus forest avifauna was assembled through time via dispersal and/or multiple vicariant events, rather than originating simultaneously via a single isolation event. Our study is one of the first multilocus, multispecies analyses revealing how ecological and migratory traits impact the evolutionary history of community formation on a continental island.

How Many Kinds of Birds Are There and Why Does It Matter?

Estimates of global species diversity have varied widely, primarily based on variation in the numbers derived from different inventory methods of arthropods and other small invertebrates. Within vertebrates, current diversity metrics for fishes, amphibians, and reptiles are known to be poor estimators, whereas those for birds and mammals are often assumed to be relatively well established. We show that avian evolutionary diversity is significantly underestimated due to a taxonomic tradition not found in most other taxonomic groups. Using a sample of 200 species taken from a list of 9159 biological species determined primarily by morphological criteria, we applied a diagnostic, evolutionary species concept to a morphological and distributional data set that resulted in an estimate of 18,043 species of birds worldwide, with a 95% confidence interval of 15,845 to 20,470. In a second, independent analysis, we examined intraspecific genetic data from 437 traditional avian species, finding an average of 2.4 evolutionary units per species, which can be considered proxies for phylogenetic species. Comparing recent lists of species to that used in this study (based primarily on morphology) revealed that taxonomic changes in the past 25 years have led to an increase of only 9%, well below what our results predict. Therefore, our molecular and morphological results suggest that the current taxonomy of birds understimates avian species diversity by at least a factor of two. We suggest that a revised taxonomy that better captures avian species diversity will enhance the quantification and analysis of global patterns of diversity and distribution, as well as provide a more appropriate framework for understanding the evolutionary history of birds.

Matching loci surveyed to questions asked in phylogeography.

Although mitochondrial DNA (mtDNA) has long been used for assessing genetic variation within and between populations, its workhorse role in phylogeography has been criticized owing to its single-locus nature. The only choice for testing mtDNA results is to survey nuclear loci, which brings into contrast the difference in locus effective size and coalescence times. Thus, it remains unclear how erroneous mtDNA-based estimates of species history might be, especially for evolutionary events in the recent past. To test the robustness of mtDNA and nuclear sequences in phylogeography, we provide one of the largest paired comparisons of summary statistics and demographic parameters estimated from mitochondrial, five Z-linked and 10 autosomal genes of 30 avian species co-distributed in the Caucasus and Europe. The results suggest that mtDNA is robust in estimating inter-population divergence but not in intra-population diversity, which is sensitive to population size change. Here, we provide empirical evidence showing that mtDNA was more likely to detect population divergence than any other single locus owing to its smaller Ne and thus faster coalescent time. Therefore, at least in birds, numerous studies that have based their inferences of phylogeographic patterns solely on mtDNA should not be readily dismissed.

Sisyphean evolution in Darwin's finches.

The trajectory of speciation involves geographic isolation of ancestral populations followed by divergence by natural selection, genetic drift or sexual selection. Once started, the process may experience fits and starts, as sometimes diverging populations intermittently reconnect. In theory populations might cycle between stages of differentiation and never attain species status, a process we refer to as Sisyphean evolution. We argue that the six putative ground finch species (genus Geospiza) of the Galápagos Islands represent a dramatic example of Sisyphean evolution that has been confused with the standard model of speciation. The dynamic environment of the Galápagos, closely spaced islands, and frequent dispersal and introgression have prevented the completion of the speciation process. We suggest that morphological clusters represent locally adapted ecomorphs, which might mimic, and have been confused with, species, but these ecomorphs do not form separate gene pools and are ephemeral in space and time. Thus the pattern of morphological, behavioural and genetic variation supports recognition of a single species of Geospiza, which we suggest should be recognized as Darwin's ground finch (Geospiza magnirostris). We argue that instead of providing an icon of insular speciation and adaptive radiation, which is featured in nearly every textbook on evolutionary biology, Darwin's ground finch represents a potentially more interesting phenomenon, one of transient morphs trapped in an unpredictable cycle of Sisyphean evolution. Instead of revealing details of the origin of species, the mechanisms underlying the transient occurrence of ecomorphs provide one of the best illustrations of the antagonistic effects of natural selection and introgression.

Fluctuating fire regimes and their historical effects on genetic variation in an endangered shrubland specialist.

The Pleistocene was characterized by worldwide shifts in community compositions. Some of these shifts were a result of changes in fire regimes, which influenced the distribution of species belonging to fire-dependent communities. We studied an endangered juniper-oak shrubland specialist, the black-capped vireo (Vireo atricapilla). This species was locally extirpated in parts of Texas and Oklahoma by the end of the 1980s as a result of habitat change and loss, predation, brood parasitism, and anthropogenic fire suppression. We sequenced multiple nuclear loci and used coalescence methods to obtain a deeper understanding of historical population trends than that typically available from microsatellites or mtDNA. We compared our estimated population history, a long-term history of the fire regime and ecological niche models representing the mid-Holocene, last glacial maximum, and last interglacial. Our Bayesian skyline plots showed a pattern of historical population fluctuation that was consistent with changing fire regimes. Genetic data suggest that the species is genetically unstructured, and that the current population should be orders of magnitude larger than it is at present. We suggest that fire suppression and habitat loss are primary factors contributing to the recent decline of the BCVI, although the role of climate change since the last glacial maximum is unclear at present.

Drastic population fluctuations explain the rapid extinction of the passenger pigeon.

To assess the role of human disturbances in species' extinction requires an understanding of the species population history before human impact. The passenger pigeon was once the most abundant bird in the world, with a population size estimated at 3-5 billion in the 1800s; its abrupt extinction in 1914 raises the question of how such an abundant bird could have been driven to extinction in mere decades. Although human exploitation is often blamed, the role of natural population dynamics in the passenger pigeon's extinction remains unexplored. Applying high-throughput sequencing technologies to obtain sequences from most of the genome, we calculated that the passenger pigeon's effective population size throughout the last million years was persistently about 1/10,000 of the 1800's estimated number of individuals, a ratio 1,000-times lower than typically found. This result suggests that the passenger pigeon was not always super abundant but experienced dramatic population fluctuations, resembling those of an "outbreak" species. Ecological niche models supported inference of drastic changes in the extent of its breeding range over the last glacial-interglacial cycle. An estimate of acorn-based carrying capacity during the past 21,000 y showed great year-to-year variations. Based on our results, we hypothesize that ecological conditions that dramatically reduced population size under natural conditions could have interacted with human exploitation in causing the passenger pigeon's rapid demise. Our study illustrates that even species as abundant as the passenger pigeon can be vulnerable to human threats if they are subject to dramatic population fluctuations, and provides a new perspective on the greatest human-caused extinction in recorded history.

Recent allopatric divergence and niche evolution in a widespread Palearctic bird, the common rosefinch (Carpodacus erythrinus).

A previously published phylogeographic analysis of mtDNA sequences from the widespread Palearctic common rosefinch (Carpodacus erythrinus) suggested the existence of three recently diverged groups, corresponding to the Caucasus, central-western Eurasia, and northeastern Eurasia. We re-evaluated the mtDNA data using coalescence methods and added sequence data from a sex-linked gene. The mtDNA gene tree and SAMOVA supported the distinctiveness of the Caucasian group but not the other two groups. However, UPGMA clustering of mtDNA Φ(ST)-values among populations recovered the three groups. The sex-linked gene tree recovered no phylogeographic signal, which was attributed to recent divergence and insufficient time for sorting of alleles. Overall, coalescence methods indicated a lack of gene flow among the three groups, and population expansion in the central-western and northeastern Eurasia groups. These three groups corresponded to named subspecies, further supporting their validity. A species distribution model revealed potential refugia at the Last Glacial Maximum. These three groups, which we hypothesized are in the early stages of speciation, provided an opportunity for testing tenets of ecological speciation. We showed that the early stages of speciation were not accompanied by ecological niche divergence, consistent with other avian studies.

Multilocus coalescence analyses support a mtDNA-based phylogeographic history for a widespread Palearctic passerine bird, Sitta europaea.

Our understanding of species phylogeography in much of the Palearctic is incomplete. In addition, many existing studies based solely on mitochondrial DNA (mtDNA) can provide a biased view of phylogeographic history because of the effects of lineage sorting, natural selection, or hybridization. We analyzed 13 introns to assess a mtDNA study of the Eurasian nuthatch (Sitta europaea) that suggested a seemingly contemporaneous origin of distinct taxa in the Caucasus, Europe, and Asia. Neutrality tests showed no evidence of selection on either the mtDNA or nuclear sequences. Most nuclear gene trees, except for Z-linked ones, did not recover the three lineages, which we attribute to recent splitting. Analyses of the 13 introns combined revealed the same three groups as did the mtDNA and suggested that nuthatches experienced a trichotomous (or two indistinguishable) split(s) 1-2 million years ago (Mya) and have remained isolated with trifling if not zero gene flow since then, and the Asian group increased in population size. This result demonstrates the usefulness of mtDNA in discovering phylogeographic patterns. The use of multiple nuclear loci facilitated detection of an introgressed individual and improved estimates of process parameters such as divergence time and population expansion. We recommend that phylogeographic studies should be based on both mtDNA and nuclear genes.

Phylogenetic relationships of the mockingbirds and thrashers (Aves: Mimidae).

The mockingbirds, thrashers and allied birds in the family Mimidae are broadly distributed across the Americas. Many aspects of their phylogenetic history are well established, but there has been no previous phylogenetic study that included all species in this radiation. Our reconstructions based on mitochondrial and nuclear DNA sequence markers show that an early bifurcation separated the Mimidae into two clades, the first of which includes North and Middle American taxa (Melanotis, Melanoptila, Dumetella) plus a small radiation that likely occurred largely within the West Indies (Ramphocinclus, Allenia, Margarops, Cinclocerthia). The second and larger radiation includes the Toxostoma thrasher clade, along with the monotypic Sage Thrasher (Oreoscoptes) and the phenotypically diverse and broadly distributed Mimus mockingbirds. This mockingbird group is biogeographically notable for including several lineages that colonized and diverged on isolated islands, including the Socorro Mockingbird (Mimus graysoni, formerly Mimodes) and the diverse and historically important Galapagos mockingbirds (formerly Nesomimus). Our reconstructions support a sister relationship between the Galapagos mockingbird lineage and the Bahama Mockingbird (M. gundlachi) of the West Indies, rather than the Long-tailed Mockingbird (M. longicaudatus) or other species presently found on the South American mainland. Relationships within the genus Toxostoma conflict with traditional arrangements but support a tree based on a preivous mtDNA study. For instance, the southern Mexican endemic Ocellated Thrasher (T. ocellatum) is not an isolated sister species of the Curve-billed thrasher (T. curvirostre).

Collision mortality has no discernible effect on population trends of North American birds.

Avian biodiversity is threatened by numerous anthropogenic factors and migratory species are especially at risk. Migrating birds frequently collide with manmade structures and such losses are believed to represent the majority of anthropogenic mortality for North American birds. However, estimates of total collision mortality range across several orders of magnitude and effects on population dynamics remain unknown. Herein, we develop a novel method to assess relative vulnerability to anthropogenic threats, which we demonstrate using 243,103 collision records from 188 species of eastern North American landbirds. After correcting mortality estimates for variation attributable to population size and geographic overlap with potential collision structures, we found that per capita vulnerability to collision with buildings and towers varied over more than four orders of magnitude among species. Species that migrate long distances or at night were much more likely to be killed by collisions than year-round residents or diurnal migrants. However, there was no correlation between relative collision mortality and long-term population trends for these same species. Thus, although millions of North American birds are killed annually by collisions with manmade structures, this source of mortality has no discernible effect on populations.

The causes of mitochondrial DNA gene tree paraphyly in birds.

Gene tree paraphyly is a potentially serious problem because many phylogenetic and phylogeographic studies assume species are monophyletic. Funk and Omland (Funk, D.J., Omland, K.E., 2003. Species-level paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annu. Rev. Ecol. Evol. Syst. 34, 397-423) found that a seemingly high proportion of bird species (16.7%) were paraphyletic in their mtDNA gene trees. This could imply that mtDNA is an unreliable or even misleading marker for delimiting species. We expand on Funk and Omland's survey and identify the causes of species-level paraphyly in birds. We find that in most cases paraphyly is caused by incorrect taxonomy. In such cases, mtDNA serves systematics by exposing and clarifying taxonomic errors. We find the next most common cause of paraphyly to be incomplete lineage sorting due to recent speciation. Here mtDNA gives a consistent picture of evolution, given the timeframe, but it is not useful for delimiting species and other criteria must be employed. There were relatively few clear instances of paraphyly due to hybridization, though there were more cases where incomplete lineage sorting and hybridization could not be distinguished. We ultimately conclude that, far from a hindrance, mtDNA is generally a useful tool that should continue to facilitate delimitation of avian species.