A phylogenetic hypothesis for passerine birds: taxonomic and biogeographic implications of an analysis of nuclear DNA sequence data

Passerine birds comprise over half of avian diversity, but have proved difficult to classify. Despite a long history of work on this group, no comprehensive hypothesis of passerine family-level relationships was available until recent analyses of DNA-DNA hybridization data. Unfortunately, given the value of such a hypothesis in comparative studies of passerine ecology and behaviour, the DNA-hybridization results have not been well tested using independent data and analytical approaches. Therefore, we analysed nucleotide sequence variation at the nuclear RAG-1 and c-mos genes from 69 passerine taxa, including representatives of most currently recognized families. In contradiction to previous DNA-hybridization studies, our analyses suggest paraphyly of suboscine passerines because the suboscine New Zealand wren Acanthisitta was found to be sister to all other passerines. Additionally, we reconstructed the parvorder Corvida as a basal paraphyletic grade within the oscine passerines. Finally, we found strong evidence that several family-level taxa are misplaced in the hybridization results, including the Alaudidae, Irenidae, and Melanocharitidae. The hypothesis of relationships we present here suggests that the oscine passerines arose on the Australian continental plate while it was isolated by oceanic barriers and that a major northern radiation of oscines (i.e. the parvorder Passerida) originated subsequent to dispersal from the south.     

A phylogenetic supertree of oscine passerine birds (Aves: Passeri)

Oscine passerine birds make up almost half of all avian diversity. Relationships within the group, and its classification, have long been controversial. Over the last 10 years numerous molecular phylogenies have been published. We compiled source phylogenies from 99 published studies to construct an oscine supertree. We aimed to illustrate weak and strong parts of the phylogeny and set targets for future phylogenetic work and therefore preferred a heuristic approach where we judged the adequacy of taxon sampling and molecular method of each source tree instead of using matrices and automated tree-building programs. We present an estimate of the phylogenetic relationships of 1723 extant and one extinct species of oscine passerine birds (Aves: Passeri) — more than 37% of the total. We included 34/35 (97%) families, 38/39 (97%) subfamilies and 40/43 (93%) tribes. Overall resolution is 83% of a fully bifurcating tree. The basal lineages are all distributed in the Australo-Papuan region, but several more distal lineages dispersed out of this region and radiated in other parts of the world. However, taxa of the Australian region suffer from larger evolutionary gaps and the deep branches of the Sylvioidea and nine South American primaried oscines are still poorly resolved.     

Multiple copies of coding as well as pseudogene c-mos sequence exist in three lacertid species

The analysis of a 581 bp section of the nuclear gene c-mos revealed multiple copies of putative functional sequences as well as pseudogenes in three closely related lacertid species Lacerta laevis, L. kulzeri and L. cyanisparsa. A phylogenetic analysis of c-mos in comparison with a molecular phylogeny based on the mitochondrial cytochrome b gene supports our findings. The study also provides new insights into the phylogenetic relationships of L. cyanisparsa and L. laevis.Pseudogenes of the three species share 11 single-nucleotide substitutions, a 1 bp deletion and a premature stop codon but differ by group-specific mutations. This result suggests that the c-mos gene has become duplicated and subsequently silenced already in the common ancestor of the three species. Sequence divergence suggests that the duplication and the loss of function occurred in the late Miocene/early Pliocene, i.e., about 5 million years ago. Indications of gene conversion are discussed.We suggest that future studies using c-mos for phylogenetic studies should provide evidence for the orthology of the sequences compared.     

Four new avian mitochondrial genomes help get to basic evolutionary questions in the late cretaceous

Good phylogenetic trees are required to test hypotheses about evolutionary processes. We report four new avian mitochondrial genomes, which together with an improved method of phylogenetic analysis for vertebrate mt genomes give results for three questions in avian evolution. The new mt genomes are: magpie goose (Anseranas semipalmata), an owl (morepork, Ninox novaeseelandiae); a basal passerine (rifleman, or New Zealand wren, Acanthisitta chloris); and a parrot (kakapo or owl-parrot, Strigops habroptilus). The magpie goose provides an important new calibration point for avian evolution because the well-studied Presbyornis fossils are on the lineage to ducks and geese, after the separation of the magpie goose. We find, as with other animal mitochondrial genomes, that RY-coding is helpful in adjusting for biases between pyrimidines and between purines. When RY-coding is used at third positions of the codon, the root occurs between paleognath and neognath birds (as expected from morphological and nuclear data). In addition, passerines form a relatively old group in Neoaves, and many modern avian lineages diverged during the Cretaceous. Although many aspects of the avian tree are stable, additional taxon sampling is required.     

Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction

Both mitochondrial and nuclear gene sequences have been employed in efforts to reconstruct deep-level phylogenetic relationships. A fundamental question in molecular systematics concerns the efficacy of different types of sequences in recovering clades at different taxonomic levels. We compared the performance of four mitochondrial data sets (cytochrome b, cytochrome oxidase II, NADH dehydrogenase subunit I, 12S rRNA-tRNA-16S rRNA) and eight nuclear data sets (exonic regions of alpha-2B adrenergic receptor, aquaporin, ss-casein, gamma-fibrinogen, interphotoreceptor retinoid binding protein, kappa-casein, protamine, von Willebrand Factor) in recovering deep-level mammalian clades. We employed parsimony and minimum-evolution with a variety of distance corrections for superimposed substitutions. In 32 different pairwise comparisons between these mitochondrial and nuclear data sets, we used the maximum set of overlapping taxa. In each case, the variable-length bootstrap was used to resample at the size of the smaller data set. The nuclear exons consistently performed better than mitochondrial protein and rRNA-tRNA coding genes on a per-residue basis in recovering benchmark clades. We also concatenated nuclear genes for overlapping taxa and made comparisons with concatenated mitochondrial protein-coding genes from complete mitochondrial genomes. The variable-length bootstrap was used to score the recovery of benchmark clades as a function of the number of resampled base pairs. In every case, the nuclear concatenations were more efficient than the mitochondrial concatenations in recovering benchmark clades. Among genes included in our study, the nuclear genes were much less affected by superimposed substitutions. Nuclear genes having appropriate rates of substitution should receive strong consideration in efforts to reconstruct deep-level phylogenetic relationships.     

The relationships and origins of the New Zealand wattlebirds (Passeriformes, Callaeatidae) from DNA sequence analyses

The monophyly of the endemic New Zealand wattlebirds (Callaeatidae) was examined through the sequencing of nuclear RAG-1 and c-mos genes and comparison to other passerine sequences. The New Zealand wattlebirds were strongly supported to be monophyletic and were nested within Corvida. An estimate for the time of divergence of the New Zealand wattlebirds indicated that the ancestors of this family arrived via transoceanic dispersal after the separation of New Zealand from Gondwana. Long branches separated the three New Zealand wattlebird genera from one another and relationships among them were unresolved, even in analyses including a further 1.5 kb of mitochondrial DNA sequences. However, most of the analyses supported either a basally diverging huia or kokako.     

Resolution of phylogenetic relationships among recently evolved species as a function of amount of DNA sequence: an empirical study based on woodpeckers (Aves: Picidae)

Synonymous substitutions in the 13 mitochondrial encoded protein genes form a large pool of characters that should approach the ideal for phylogenetic analysis of being independently and identically distributed. Pooling sequences from multiple mitochondrial protein-coding genes should result in statistically more powerful estimates of relationships among species that diverged sufficiently recently that most nucleotide substitutions are synonymous. Cytochrome oxidase I (COI) was sequenced for woodpecker species for which cytochrome b (cyt b) sequences were available. A pairing-design test based on the normal distribution indicated that cyt b evolves more rapidly than COI when all nucleotides are compared but their rates are equal for synonymous substitutions. Nearly all of the phylogenetically informative substitutions among woodpeckers are synonymous. Statistical support for relationships, as measured by bootstrap proportions, increased as the number of nucleotides increased from 1047 (cyt b) to 1512 (COI) to 2559 nucleotides (aggregate data set). Pseudo-bootstrap replicates showed the same trend and increasing the amount of sequence beyond the actual length of 2559 nucleotides to 5120 (2x) resulted in stronger bootstrap support, even though the amount of phylogenetic information was the same. However, the amount of sequence required to resolve an internode depends on the length of the internode and its depth in the phylogeny.     

Molecular phylogenetics of finches and sparrows: consequences of character state removal in cytochrome b sequences

The complete mitochondrial cytochrome b genes of 53 genera of oscine passerine birds representing the major groups of finches and some allies were compared. Phylogenetic trees resulting from three levels of character partition removal (no data removed, transitions at third positions of codons removed, and all transitions removed [transversion parsimony]) were generally concordant, and all supported several basic statements regarding relationships of finches and finch-like birds, including: (1) larks (Alaudidae) show no close relationship to any finch group; (2) Peucedramus (olive warbler) is phylogenetically far removed from true wood warblers; (3) a clade consisting of fringillids, passerids, motacillids, and emberizids is supported, and this clade is characterized by evolution of a vestigial 10th wing primary; and (4) Hawaiian honeycreepers are derived from within the cardueline finches. Excluding transition substitutions at third positions of codons resulted in phylogenetic trees similar to, but with greater bootstrap nodal support than, trees derived using either all data (equally weighted) or transversion parsimony. Relative to the shortest trees obtained using all data, the topologies obtained after elimination of third-position transitions showed only slight increases in realized treelength and homoplasy. These increases were negligable compared to increases in overall nodal support; therefore, this partition removal scheme may enhance recovery of deep phylogenetic signal in protein-coding DNA datasets.     

Patterns of song evolution and sexual selection in the oropendolas and caciques

Although oscine bird song is widely thought to have evolvedunder the influence of sexual selection, few studies have usedphylogenetic comparative methods to investigate how these vocalizationshave changed historically. In the present study, we use a molecularphylogeny based on mitochondrial sequence data to reconstructvocal evolution in the oropendolas and caciques, an oscine groupwith diverse taxon-specific song patterns and a wide range inlevels of sexual size dimorphism. Our reconstructions show thatlarge changes in song organization and structure have occurredon branches of the phylogeny with relatively high levels ofsize dimorphism. The particular vocal components that changed,however, often differed in different phylogenetic lineages.These patterns indicate that sexual selection has had importantinfluences on song evolution in these birds, but has targeteddifferent components of song in different taxa. Our resultsprovide insight into how sexual selection influences bird songand suggest directions for future research to uncover the behavioralmechanisms driving vocal evolution.     

Molecular Phylogenetics of Finches and Sparrows: Consequences of Character State Removal in CytochromebSequences

The complete mitochondrial cytochromebgenes of 53 genera of oscine passerine birds representing the major groups of finches and some allies were compared. Phylogenetic trees resulting from three levels of character partition removal (no data removed, transitions at third positions of codons removed, and all transitions removed [transversion parsimony]) were generally concordant, and all supported several basic statements regarding relationships of finches and finch-like birds, including: (1) larks (Alaudidae) show no close relationship to any finch group; (2)Peucedramus(olive warbler) is phylogenetically far removed from true wood warblers; (3) a clade consisting of fringillids, passerids, motacillids, and emberizids is supported, and this clade is characterized by evolution of a vestigial 10th wing primary; and (4) Hawaiian honeycreepers are derived from within the cardueline finches. Excluding transition substitutions at third positions of codons resulted in phylogenetic trees similar to, but with greater bootstrap nodal support than, trees derived using either all data (equally weighted) or transversion parsimony. Relative to the shortest trees obtained using all data, the topologies obtained after elimination of third-position transitions showed only slight increases in realized treelength and homoplasy. These increases were negligable compared to increases in overall nodal support; therefore, this partition removal scheme may enhance recovery of deep phylogenetic signal in protein-coding DNA datasets.     

Rates and patterns of mitochondrial DNA sequence evolution in fringilline finches (fringilla spp.) and the greenfinch (Carduelis chloris)

Rates and patterns of evolution in partial sequences of five mitochondrial genes (cytochrome b, ATPase 6, NADH dehydrogenase subunit 5, tRNA(Glu), and the control region) were compared among taxa in the passerine bird genera Fringilla and Carduelis. Rates of divergence do not vary significantly among genes, even in comparisons with the control region. Rate variation among lineages is significant only for the control region and NADH dehydrogenase subunit 5, and patterns of variation are consistent with the expectations of neutral theory. Base composition is biased in all genes but is stationary among lineages, and there is evidence for directional mutation pressure only in the control region. Despite these similarities, patterns of substitution differ among genes, consistent with alternative regimes of selective constraint. Rates of nonsynonymous substitution are higher in NADH dehydrogenase subunit 5 than in other protein-coding genes, and transitions exist in elevated proportions relative to transversions. Transitions appear to accumulate linearly with time in tRNA(Glu), and despite exhibiting the highest overall rate of divergence among species, there are no transversional changes in this gene. Finally, for resolving phylogenetic relationships among Fringilla taxa, the combined protein-coding data are broadly similar to those of the control region in terms of phylogenetic informativeness and statistical support.      

Ultraviolet visual sensitivity in three avian lineages: paleognaths, parrots, and passerines

Ultraviolet (UV) light-transmitted signals play a major role in avian foraging and communication, subserving functional roles in feeding, mate choice, egg recognition, and nestling discrimination. Sequencing functionally relevant regions of the short wavelength sensitive type 1 (SWS1) opsin gene that is responsible for modulating the extent of SWS1 UV sensitivity in birds allows predictions to be made about the visual system's UV sensitivity in species where direct physiological or behavioral measures would be impractical or unethical. Here, we present SWS1 segment sequence data from representative species of three avian lineages for which visually based cues for foraging and communication have been investigated to varying extents. We also present a preliminary phylogenetic analysis and ancestral character state reconstructions of key spectral tuning sites along the SWS1 opsin based on our sequence data. The results suggest ubiquitous ultraviolet SWS1 sensitivity (UVS) in both paleognaths, including extinct moa (Emeidae), and parrots, including the nocturnal and flightless kakapo (Strigops habroptilus), and in most, but not all, songbird (oscine) lineages, and confirmed violet sensitivity (VS) in two suboscine families. Passerine hosts of avian brood parasites were included both UVS and VS taxa, but sensitivity did not co-vary with egg rejection behaviors. The results should stimulate future research into the functional parallels between the roles of visual signals and the genetic basis of visual sensitivity in birds and other taxa.     

Nuclear corroboration of DNA-DNA hybridization in deep phylogenies of hummingbirds, swifts, and passerines: the phylogenetic utility of ZENK (ii)

This paper documents the phylogenetic utility of ZENK at the avian intra-ordinal level using hummingbirds, swifts, and passerines as case studies. ZENK sequences (1.7 kb) were used to reconstruct separate gene trees containing the major lineages of each group, and the three trees were examined for congruence with existing DNA-DNA hybridization trees. The results indicate both that ZENK is an appropriate nuclear marker for resolving relationships deep in the avian tree, and that many relationships within these three particular groups are congruent among the different datasets. Specifically, within hummingbirds there was topological agreement that the major hummingbird lineages diverged in a graded manner from the "hermits," to the "mangoes," to the "coquettes," to the "emeralds," and finally to a sister relationship between the "mountain-gems" and the "bees." Concerning swifts, the deepest divergences were congruent: treeswifts (Hemiprocnidae) were sister to the typical swifts (Apodidae), and the subfamily Apodinae was monophyletic relative to Cypseloidinae. Within Apodinae, however, were short, unresolved branches among the swiftlets, spinetails, and more typical swifts; a finding which coincides with other datasets. Within passerine birds, there was congruent support for monophyly of sub-oscines and oscines, and within sub-oscines, for monophyly of New World groups relative to the Old World lineages. New World sub-oscines split into superfamilies Furnaroidea and Tyrannoidea, with the Tyrannoid relationships completely congruent among ZENK and DNA-DNA hybridization trees. Within Furnaroidea, however, there was some incongruence regarding the positions of Thamnophilidae and Formicariidae. Concerning oscine passerines, both datasets showed a split between Corvida and Passerida and confirmed the traditional membership of passerid superfamilies Muscicapoidea and Passeroidea. Monophyly of Sylvioidea, however, remained uncertain, as did the relationships among the superfamiles themselves. These results are strikingly similar to other recent findings and indicative of continuing uncertainty about the higher level relationships of oscine passerines.     

Two new avian mitochondrial genomes (penguin and goose) and a summary of bird and reptile mitogenomic features

We report complete mitochondrial (mt) genomes for a penguin (little blue, Eudyptula minor) and a goose (greater white-fronted, Anser albifrons). A revised annotation of avian and reptile mt genomes has been carried out, which improves consistency of labeling gene start and stop positions. In conjunction with this, a summary of mt gene features is presented and a number of conserved patterns and interesting differences identified. The protein-coding genes from the two new genomes were analysed together with those from 17 other birds plus outgroup (reptile) taxa. The unrooted amino acid tree from 19 avian genomes was locally stable with many high bootstrap values using several maximum likelihood methods. In particular, Anseriformes (goose and duck) grouped strongly with Galliformes (chicken) to form Gallianseres, while the penguin paired firmly with the stork. The position where the outgroup joined the avian tree varied with the combination of outgroup taxa used. The three best supported positions of the root were passerine, but the traditional rooting position between paleognaths and neognaths could not be excluded.     

Phylogeny and biogeography of Oriolidae (Aves: Passeriformes)

Understanding oscine passerine dispersal patterns out of their Australian area of origin is hampered by a paucity of robust phylogenies. We constructed a molecular phylogeny of the oscine family, Oriolidae, which is distributed from Australia through to the Old World. We used the phylogeny to assess direction and timing of dispersal and whether dispersal can be linked with the well‐documented movements of geological terranes in the Indonesian Archipelago. We sampled 29 of 33 species of Oriolidae from fresh tissue and from toe pads from museum specimens, and examined two nuclear introns and two mitochondrial genes. Model‐based phylogenetic analyses yielded strong support for clades that generally mirrored classical systematics. Biogeographical analyses and divergence time estimates demonstrated that the family originated in the Australo‐Papuan region from where it dispersed first to Asia and then onwards to Africa and the Philippines before back‐colonising Asia and the Indonesian archipelago. Thus, contrary to several other avian families in the region, Oriolidae represents a sequential dispersal pattern from Australia to Africa via Asia. However, it is noteworthy that the Pacific islands and archipelagos remain uncolonised and that members inhabiting Wallacea are recent colonisers suggesting that Oriolidae are poorly adapted to island life.     

Conserved sequence motifs, alignment, and secondary structure for the third domain of animal 12S rRNA

Secondary structure models are an important step for aligning sequences, understanding probabilities of nucleotide substitutions, and evaluating the reliability of phylogenetic reconstructions. A set of conserved sequence motifs is derived from comparative sequence analysis of 184 invertebrate and vertebrate taxa (including many taxa from the same genera, families, and orders) with reference to a secondary structure model for domain III of animal mitochondrial small subunit (12S) ribosomal RNA. A template is presented to assist with secondary structure drawing. Our model is similar to previous models but is more specific to mitochondrial DNA, fitting both invertebrate and vertebrate groups, including taxa with markedly different nucleotide compositions. The second half of the domain III sequence can be difficult to align precisely, even when secondary structure information is considered. This is especially true for comparisons of anciently diverged taxa, but well-conserved motifs assist in determining biologically meaningful alignments. Patterns of conservation and variability in both paired and unpaired regions make differential phylogenetic weighting in terms of "stems" and "loops" unsatisfactory. We emphasize looking carefully at the sequence data before and during analyses, and advocate the use of conserved motifs and other secondary structure information for assessing sequencing fidelity.      

Use of mitogenomic information in teleostean molecular phylogenetics: a tree-based exploration under the maximum-parsimony optimality criterion

We explored the phylogenetic utility and limits of the individual and concatenated mitochondrial genes for reconstructing the higher-level relationships of teleosts, using the complete (or nearly complete) mitochondrial DNA sequences of eight teleosts (including three newly determined sequences), whose relative phylogenetic positions were noncontroversial. Maximum-parsimony analyses of the nucleotide and amino acid sequences of 13 protein-coding genes from the above eight teleosts, plus two outgroups (bichir and shark), indicated that all of the individual protein-coding genes, with the exception of ND5, failed to recover the expected phylogeny, although unambiguously aligned sequences from 22 concatenated transfer RNA (tRNA) genes (stem regions only) recovered the expected phylogeny successfully with moderate statistical support. The phylogenetic performance of the 13 protein-coding genes in recovering the expected phylogeny was roughly classified into five groups, viz. very good (ND5, ND4, COIII, COI), good (COII, cyt b), medium (ND3, ND2), poor (ND1, ATPase 6), and very poor (ND4L, ND6, ATPase 8). Although the universality of this observation was unclear, analysis of successive concatenation of the 13 protein-coding genes in the same ranking order revealed that the combined data sets comprising nucleotide sequences from the several top-ranked protein-coding genes (no 3rd codon positions) plus the 22 concatenated tRNA genes (stem regions only) best recovered the expected phylogeny, with all internal branches being supported by bootstrap values >90%. We conclude that judicious choice of mitochondrial genes and appropriate data weighting, in conjunction with purposeful taxonomic sampling, are prerequisites for resolving higher-level relationships in teleosts under the maximum-parsimony optimality criterion.     

Congruent Avian Phylogenies Inferred from Mitochondrial and Nuclear DNA Sequences

Recent molecular studies addressing the phylogenetic relationships of avian orders have had conflicting results. While studies using nuclear DNA sequences tend to support traditional taxonomic views, also supported by morphological data [(paleognaths (galloanseres (all other birds)))], with songbirds forming a clade within Neoaves (all other birds), analyses with complete mtDNA genomes have resulted in topologies that place songbirds as one of the earliest-diverging avian lineages. Considering that over half of the extant bird species are songbirds, these different results have very different implications for our understanding of avian evolution. We analyzed data sets comprising nearly 4 kb of mitochondrial DNA (mtDNA) (complete 12S, ND1, ND2, and cytochrome b) plus 600 bp of the nuclear gene c-mos for 15 birds that were chosen to represent all major avian clades and to minimize potential long-branch attraction problems; we used a partition-specific maximum likelihood approach. Our results show congruence with respect to the ingroup among phylogenies obtained with mtDNA and the nuclear gene c-mos, separately or combined. The data sets support a traditional avian taxonomy, with paleognaths (ratites and tinamous) occupying a basal position and with songbirds more derived and forming a monophyletic group. We also show that, for mtDNA studies, turtles may be a better outgroup for birds than crocodilians because of their slower rate of sequence evolution.     

Genetic evidence supports song learning in the three-wattled bellbird Procnias tricarunculata (Cotingidae)

Vocal learning is thought to have evolved in three clades of birds (parrots, hummingbirds, and oscine passerines), and three clades of mammals (whales, bats, and primates). Behavioural data indicate that, unlike other suboscine passerines, the three-wattled bellbird Procnias tricarunculata (Cotingidae) is capable of vocal learning. Procnias tricarunculata shows conspicuous vocal ontogeny, striking geographical variation in song, and rapid temporal change in song within a population. Deprivation studies of vocal development in P. tricarunculata are impractical. Here, we report evidence from mitochondrial DNA sequences and nuclear microsatellite loci that genetic variation within and among the four allopatric breeding populations of P. tricarunculata is not congruent with variation in vocal behaviour. Sequences of the mitochondrial DNA control region document extensive haplotype sharing among localities and song types, and no phylogenetic resolution of geographical populations or behavioural groups. The vocally differentiated, allopatric breeding populations of P. tricarunculata are only weakly genetically differentiated populations, and are not distinct taxa. Mitochondrial DNA and microsatellite variation show small (2.9% and 13.5%, respectively) but significant correlation with geographical distance, but no significant residual variation by song type. Estimates of the strength of selection that would be needed to maintain the observed geographical pattern in vocal differentiation if songs were genetically based are unreasonably high, further discrediting the hypothesis of a genetic origin of vocal variation. These data support a fourth, phylogenetically independent origin of avian vocal learning in Procnias. Geographical variations in P. tricarunculata vocal behaviour are likely culturally evolved dialects.