Bird evolution: testing the Metaves clade with six new mitochondrial genomes

New view on the organization and evolution of Palaeognathae mitogenomes poses the question on the ancestral gene rearrangement in Aves

BACKGROUND

Bird mitogenomes differ from other vertebrates in gene rearrangement. The most common avian gene order, identified first in Gallus gallus, is considered ancestral for all Aves. However, other rearrangements including a duplicated control region and neighboring genes have been reported in many representatives of avian orders. The repeated regions can be easily overlooked due to inappropriate DNA amplification or genome sequencing. This raises a question about the actual prevalence of mitogenomic duplications and the validity of the current view on the avian mitogenome evolution. In this context, Palaeognathae is especially interesting because is sister to all other living birds, i.e. Neognathae. So far, a unique duplicated region has been found in one palaeognath mitogenome, that of Eudromia elegans.

RESULTS

Therefore, we applied an appropriate PCR strategy to look for omitted duplications in other palaeognaths. The analyses revealed the duplicated control regions with adjacent genes in Crypturellus, Rhea and Struthio as well as ND6 pseudogene in three moas. The copies are very similar and were subjected to concerted evolution. Mapping the presence and absence of duplication onto the Palaeognathae phylogeny indicates that the duplication was an ancestral state for this avian group. This feature was inherited by early diverged lineages and lost two times in others. Comparison of incongruent phylogenetic trees based on mitochondrial and nuclear sequences showed that two variants of mitogenomes could exist in the evolution of palaeognaths. Data collected for other avian mitogenomes revealed that the last common ancestor of all birds and early diverging lineages of Neoaves could also possess the mitogenomic duplication.

CONCLUSIONS

The duplicated control regions with adjacent genes are more common in avian mitochondrial genomes than it was previously thought. These two regions could increase effectiveness of replication and transcription as well as the number of replicating mitogenomes per organelle. In consequence, energy production by mitochondria may be also more efficient. However, further physiological and molecular analyses are necessary to assess the potential selective advantages of the mitogenome duplications.

Characterization of three new mitochondrial genomes of Coraciiformes (Megaceryle lugubris, Alcedo atthis, Halcyon smyrnensis) and insights into their phylogenetics

Coraciiformes contains more than 200 species with great differences on external morphology and life-style. The evolutionary relationships within Coraciiformes and the phylogenetic placement of Coraciiformes in Aves are still questioned. Mitochondrial genome (mitogenome) sequences are popular markers in molecular phylogenetic studies of birds. This study presented the genome characteristics of three new mitogenomes in Coraciiformes and explored the phylogenetic relationships among Coraciiformes and other five related orders with mitogenome data of 30 species. The sizes of three mitogenomes were 17,383 bp (Alcedo atthis), 17,892 bp (Halcyon smyrnensis) and 17,223 bp (Megaceryle lugubris). Each mitogenome contained one control region and 37 genes that were common in vertebrate mitogenomes. The organization of three mitogenomes was identical to the putative ancestral gene order in Aves. Among 13 available Coraciiform mitogenomes, 12 protein coding genes showed indications of negative selection, while the MT-ND6 presented sign of positive selection or relaxed purifying selection. The phylogenetic results supported that Upupidae and Bucerotidae should be separated from Coraciiformes, and that Coraciiformes is more closely related to Piciformes than to Strigiformes, Trogoniformes and Cuculiformes. Our study provide valuable data for further phylogenetic investigation of Coraciiformes.

Deep-Time Demographic Inference Suggests Ecological Release as Driver of Neoavian Adaptive Radiation

Assessing the applicability of theory to major adaptive radiations in deep time represents an extremely difficult problem in evolutionary biology. Neoaves, which includes 95% of living birds, is believed to have undergone a period of rapid diversification roughly coincident with the Cretaceous–Paleogene (K-Pg) boundary. We investigate whether basal neoavian lineages experienced an ecological release in response to ecological opportunity, as evidenced by density compensation. We estimated effective population sizes (Ne) of basal neoavian lineages by combining coalescent branch lengths (CBLs) and the numbers of generations between successive divergences. We used a modified version of Accurate Species TRee Algorithm (ASTRAL) to estimate CBLs directly from insertion–deletion (indel) data, as well as from gene trees using DNA sequence and/or indel data. We found that some divergences near the K-Pg boundary involved unexpectedly high gene tree discordance relative to the estimated number of generations between speciation events. The simplest explanation for this result is an increase in Ne, despite the caveats discussed herein. It appears that at least some early neoavian lineages, similar to the ancestor of the clade comprising doves, mesites, and sandgrouse, experienced ecological release near the time of the K-Pg mass extinction.

Resolving Phylogenetic Relationships within Passeriformes Based on Mitochondrial Genes and Inferring the Evolution of Their Mitogenomes in Terms of Duplications

Mitochondrial genes are placed on one molecule, which implies that they should carry consistent phylogenetic information. Following this advantage, we present a well-supported phylogeny based on mitochondrial genomes from almost 300 representatives of Passeriformes, the most numerous and differentiated Aves order. The analyses resolved the phylogenetic position of paraphyletic Basal and Transitional Oscines. Passerida occurred divided into two groups, one containing Paroidea and Sylvioidea, whereas the other, Passeroidea and Muscicapoidea. Analyses of mitogenomes showed four types of rearrangements including a duplicated control region (CR) with adjacent genes. Mapping the presence and absence of duplications onto the phylogenetic tree revealed that the duplication was the ancestral state for passerines and was maintained in early diverged lineages. Next, the duplication could be lost and occurred independently at least four times according to the most parsimonious scenario. In some lineages, two CR copies have been inherited from an ancient duplication and highly diverged, whereas in others, the second copy became similar to the first one due to concerted evolution. The second CR copies accumulated over twice as many substitutions as the first ones. However, the second CRs were not completely eliminated and were retained for a long time, which suggests that both regions can fulfill an important role in mitogenomes. Phylogenetic analyses based on CR sequences subjected to the complex evolution can produce tree topologies inconsistent with real evolutionary relationships between species. Passerines with two CRs showed a higher metabolic rate in relation to their body mass.

Phylogenetic Signal of Indels and the Neoavian Radiation

The early radiation of Neoaves has been hypothesized to be an intractable “hard polytomy”. We explore the fundamental properties of insertion/deletion alleles (indels), an under-utilized form of genomic data with the potential to help solve this. We scored >5 million indels from >7000 pan-genomic intronic and ultraconserved element (UCE) loci in 48 representatives of all neoavian orders. We found that intronic and UCE indels exhibited less homoplasy than nucleotide (nt) data. Gene trees estimated using indel data were less resolved than those estimated using nt data. Nevertheless, Accurate Species TRee Algorithm (ASTRAL) species trees estimated using indels were generally similar to nt-based ASTRAL trees, albeit with lower support. However, the power of indel gene trees became clear when we combined them with nt gene trees, including a striking result for UCEs. The individual UCE indel and nt ASTRAL trees were incongruent with each other and with the intron ASTRAL trees; however, the combined indel+nt ASTRAL tree was much more congruent with the intronic trees. Finally, combining indel and nt data for both introns and UCEs provided sufficient power to reduce the scope of the polytomy that was previously proposed for several supraordinal lineages of Neoaves.

New Insight into Parrots' Mitogenomes Indicates That Their Ancestor Contained a Duplicated Region

Mitochondrial genomes of vertebrates are generally thought to evolve under strong selection for size reduction and gene order conservation. Therefore, a growing number of mitogenomes with duplicated regions changes our view on the genome evolution. Among Aves, order Psittaciformes (parrots) is especially noteworthy because of its large morphological, ecological, and taxonomical diversity, which offers an opportunity to study genome evolution in various aspects. Former analyses showed that tandem duplications comprising the control region with adjacent genes are restricted to several lineages in which the duplication occurred independently. However, using an appropriate polymerase chain reaction strategy, we demonstrate that early diverged parrot groups contain mitogenomes with the duplicated region. These findings together with mapping duplication data from other mitogenomes onto parrot phylogeny indicate that the duplication was an ancestral state for Psittaciformes. The state was inherited by main parrot groups and was lost several times in some lineages. The duplicated regions were subjected to concerted evolution with a frequency higher than the rate of speciation. The duplicated control regions may provide a selective advantage due to a more efficient initiation of replication or transcription and a larger number of replicating genomes per organelle, which may lead to a more effective energy production by mitochondria. The mitogenomic duplications were associated with phenotypic features and parrots with the duplicated region can live longer, show larger body mass as well as predispositions to a more active flight. The results have wider implications on the presence of duplications and their evolution in mitogenomes of other avian groups.

Cancellous bone and theropod dinosaur locomotion. Part I-an examination of cancellous bone architecture in the hindlimb bones of theropods

This paper is the first of a three-part series that investigates the architecture of cancellous ('spongy') bone in the main hindlimb bones of theropod dinosaurs, and uses cancellous bone architectural patterns to infer locomotor biomechanics in extinct non-avian species. Cancellous bone is widely known to be highly sensitive to its mechanical environment, and has previously been used to infer locomotor biomechanics in extinct tetrapod vertebrates, especially primates. Despite great promise, cancellous bone architecture has remained little utilized for investigating locomotion in many other extinct vertebrate groups, such as dinosaurs. Documentation and quantification of architectural patterns across a whole bone, and across multiple bones, can provide much information on cancellous bone architectural patterns and variation across species. Additionally, this also lends itself to analysis of the musculoskeletal biomechanical factors involved in a direct, mechanistic fashion. On this premise, computed tomographic and image analysis techniques were used to describe and analyse the three-dimensional architecture of cancellous bone in the main hindlimb bones of theropod dinosaurs for the first time. A comprehensive survey across many extant and extinct species is produced, identifying several patterns of similarity and contrast between groups. For instance, more stemward non-avian theropods (e.g. ceratosaurs and tyrannosaurids) exhibit cancellous bone architectures more comparable to that present in humans, whereas species more closely related to birds (e.g. paravians) exhibit architectural patterns bearing greater similarity to those of extant birds. Many of the observed patterns may be linked to particular aspects of locomotor biomechanics, such as the degree of hip or knee flexion during stance and gait. A further important observation is the abundance of markedly oblique trabeculae in the diaphyses of the femur and tibia of birds, which in large species produces spiralling patterns along the endosteal surface. Not only do these observations provide new insight into theropod anatomy and behaviour, they also provide the foundation for mechanistic testing of locomotor hypotheses via musculoskeletal biomechanical modelling.

Closing the gap: Avian lineage splits at a young, narrow seaway imply a protracted history of mixed population response

The evolutionary significance of spatial habitat gaps has been well recognized since Alfred Russel Wallace compared the faunas of Bali and Lombok. Gaps between islands influence population structuring of some species, and flightless birds are expected to show strong partitioning even where habitat gaps are narrow. We examined the population structure of the most numerous living flightless land bird in New Zealand, Weka (Gallirallus australis). We surveyed Weka and their feather lice in native and introduced populations using genetic data gathered from DNA sequences of mitochondrial genes and nuclear β-fibrinogen and five microsatellite loci. We found low genetic diversity among extant Weka population samples. Two genetic clusters were evident in the mtDNA from Weka and their lice, but partitioning at nuclear loci was less abrupt. Many formerly recognized subspecies/species were not supported; instead, we infer one subspecies for each of the two main New Zealand islands. Although currently range restricted, North Island Weka have higher mtDNA diversity than the more wide-ranging southern Weka. Mismatch and neutrality statistics indicate North Island Weka experienced rapid and recent population reduction, while South Island Weka display the signature of recent expansion. Similar haplotype data from a widespread flying relative of Weka and other New Zealand birds revealed instances of North Island-South Island partitioning associated with a narrow habitat gap (Cook Strait). However, contrasting patterns indicate priority effects and other ecological factors have a strong influence on spatial exchange at this scale.

The mitochondrial genome of pin-tailed snipe Gallinago stenura, and its implications for the phylogeny of Charadriiformes

The Charadriiformes, among the most diverse orders of bird, is a good source to research on evolution. The mitochondrial genome sequencing database has rapidly increased in recent years, while Charadriiformes mitogenome has not been well studied. In this research, we determined the complete mitogenome sequence of Gallinago stenura, and comparatively analysed 20 mitogenomes of Charadriiformes. The mitogenomes display moderate size variation, and most of variation due to mutations in the control region. In 13 protein-coding genes, we found: 1. The GC skews are always negative, while the negative AT skews are found in 5 genes, 2. The average uncorrected pairwise distances reveal heterogeneity of evolutionary rate for each gene, 3. The ATG and TAA, respectively, are observed the most commonly start and stop codon. The highest dN/dS is detected for ATP8 (0.16) among Charadriiformes, while the lowest for COI (0.01), indicating that 13 protein-coding genes are evolving under the purifying selection. Predicted secondary structures of tRNAs indicate that the sequences and structures of anticodon, amino acceptor, and TψC arms are highly conserved, and most nucleotide variation is restricted to dihydrouridine arms with obvious indel polymorphisms. A total of 15 conserved sequence boxes were recognized in the control regions, and the 4 bp (5'-AAAC-3') and 7 bp (5'- AAACAAC -3') repeat sequences occurred frequently. Phylogenomic analysis based on the nearly complete mitochondrial genomes strongly supported the monophyly of the order, and the suborder Charadrii is at the basal of Charadriiformes. Moreover, our results well resolved the complexity family-level relationships and clearly depicted the evolutionary processes of Charadriiformes, based on 12 mitochondrial protein-coding genes from 18 families. This study improves our understanding of mitogenomic structure and evolution, which can provide further insights into our understanding of phylogeny and taxonomy in Charadriiformes.

The mitochondrial genome of the sooty swift (Cypseloides fumigatus)

We assembled the mitogenome of Cypseloides fumigatus based on off-target sequences from ultraconserved elements sequencing. We found a total length of 16,850 bp, including 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and one control region, organized in the standard avian gene order. We have built a phylogenetic tree including 26 species of swifts that suggested C. fumigatus as sister species of C. cryptus, and indicated exciting opportunities for biogeographic inferences involving most continents, including Neartic vs Neotropical disjunctions and local radiations across the globe. Finally, we found cases of lack of reciprocal monophyly between named species and high intra-specific divergence, suggesting that population-level studies are warranted.

Evolutionary pathways to convergence in plumage patterns

Background

Avian plumage is ideal for investigating phenotypic convergence because of repeated evolution of the same within-feather patterns. In birds, there are three major types of regular patterns within feathers: scales, bars and spots. Existing models of within-feather pattern development suggest that scales have the simplest developmental mechanism, bars require more stringent regulation than scales, and spots have the strictest developmental parameters. We hypothesized that increasing stringency in the mechanism of pattern formation predicts the evolutionary trajectory of patterns, and hence scales should evolve first, followed by bars and finally spots. Here, using Bayesian phylogenetic modeling we reconstructed pattern evolution in the most spectacularly patterned avian clades – aquatic waterfowl (Anseriformes) and terrestrial gamebirds (Galliformes).

Results

Our analyses suggest that the ancestral state of plumage is an absence of patterns, but with some variability. Independent analyses of seven feather patches reveal that spots evolve after bars and scales. However, both scales and bars evolve frequently from an absence of patterns, contradicting our predictions. Over the whole body, many constraints are conserved from the level of patches, for example the largest number of steps from the ancestral state was required for spots to evolve.

Conclusions

Overall there was remarkable similarity in the inferred evolutionary trajectories of plumage pattern evolution in Galliformes and Anseriformes, suggesting that developmental constraint is similar in these two orders, despite large ecological differences. These evolutionary transitions are largely congruent with a reaction–diffusion based model of pattern formation, but the evolution of bars from an unpatterned ancestor is more common than expected. Our study highlights the promise of testing models of development using comparative methods.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0741-x) contains supplementary material, which is available to authorized users.

Clustering DNA sequences using the out-of-place measure with reduced n-grams

The alignment-free n-gram based method with the out-of-place measures as the distance has been successfully applied to automatic text or natural languages categorization in real time. However, it is not clear about its performance and the selection of n for comparing genome sequences. Here we propose a symmetric version of the out-of-place measure and a new approach for finding the optimal range of n to construct a phylogenetic tree with the symmetric out-of-place measures. Our method is then applied to real genome sequence datasets. The resulting phylogenetic trees are matching with the standard biological classification. It shows that our proposed method is a very powerful tool for phylogenetic analysis in terms of both classification accuracy and computation efficiency.

Application of multigene phylogenetics and site-stripping to resolve intraordinal relationships in the Rhodymeniales (Rhodophyta)

Previous molecular assessments of the red algal order Rhodymeniales have confirmed its monophyly and distinguished the six currently recognized families (viz. Champiaceae, Faucheaceae, Fryeellaceae, Hymenocladiaceae, Lomentariaceae, and Rhodymeniaceae); however, relationships among most of these families have remained unresolved possibly as a result of substitution saturation at deeper phylogenetic nodes. The objective of the current study was to improve rhodymenialean systematics by increasing taxonomic representation and using a more robust multigene dataset of mitochondrial (COB, COI/COI-5P), nuclear (LSU, EF2) and plastid markers (psbA, rbcL). Additionally, we aimed to prevent phylogenetic inference problems associated with substitution saturation (particularly at the interfamilial nodes) by removing fast-evolving sites and analyzing a series of progressively more conservative alignments. The Rhodymeniales was resolved as two major lineages: (i) the Fryeellaceae as sister to the Faucheaceae and Lomentariaceae; and (ii) the Rhodymeniaceae allied to the Champiaceae and Hymenocladiaceae. Support at the interfamilial nodes was highest when 20% of variable sites were removed. Inclusion of Binghamiopsis, Chamaebotrys, and Minium, which were absent in previous phylogenetic investigations, established their phylogenetic affinities while assessment of two genera consistently polyphyletic in phylogenetic analyses, Erythrymenia and Lomentaria, resulted in the proposition of the novel genera Perbella and Fushitsunagia. The taxonomic position of Drouetia was reinvestigated with re-examination of holotype material of D. coalescens to clarify tetrasporangial development in this genus. In addition, we added three novel Australian species to Drouetia as a result of ongoing DNA barcoding assessments-D. aggregata sp. nov., D. scutellata sp. nov., and D. viridescens sp. nov.

Determining the Position of Storks on the Phylogenetic Tree of Waterbirds by Retroposon Insertion Analysis

Despite many studies on avian phylogenetics in recent decades that used morphology, mitochondrial genomes, and/or nuclear genes, the phylogenetic positions of several birds (e.g., storks) remain unsettled. In addition to the aforementioned approaches, analysis of retroposon insertions, which are nearly homoplasy-free phylogenetic markers, has also been used in avian phylogenetics. However, the first step in the analysis of retroposon insertions, that is, isolation of retroposons from genomic libraries, is a costly and time-consuming procedure. Therefore, we developed a high-throughput and cost-effective protocol to collect retroposon insertion information based on next-generation sequencing technology, which we call here the STRONG (Screening of Transposons Obtained by Next Generation Sequencing) method, and applied it to 3 waterbird species, for which we identified 35,470 loci containing chicken repeat 1 retroposons (CR1). Our analysis of the presence/absence of 30 CR1 insertions demonstrated the intra- and interordinal phylogenetic relationships in the waterbird assemblage, namely 1) Loons diverged first among the waterbirds, 2) penguins (Sphenisciformes) and petrels (Procellariiformes) diverged next, and 3) among the remaining families of waterbirds traditionally classified in Ciconiiformes/Pelecaniformes, storks (Ciconiidae) diverged first. Furthermore, our genome-scale, in silico retroposon analysis based on published genome data uncovered a complex divergence history among pelican, heron, and ibis lineages, presumably involving ancient interspecies hybridization between the heron and ibis lineages. Thus, our retroposon-based waterbird phylogeny and the established phylogenetic position of storks will help to understand the evolutionary processes of aquatic adaptation and related morphological convergent evolution.

New Zealand Passerines Help Clarify the Diversification of Major Songbird Lineages during the Oligocene

Passerines are the largest avian order, and the 6,000 species comprise more than half of all extant bird species. This successful radiation probably had its origin in the Australasian region, but dating this origin has been difficult due to a scarce fossil record and poor biogeographic assumptions. Many of New Zealand’s endemic passerines fall within the deeper branches of the passerine radiation, and a well resolved phylogeny for the modern New Zealand element in the deeper branches of the oscine lineage will help us understand both oscine and passerine biogeography. To this end we present complete mitochondrial genomes representing all families of New Zealand passerines in a phylogenetic framework of over 100 passerine species. Dating analyses of this robust phylogeny suggest Passeriformes originated in the early Paleocene, with the major lineages of oscines “escaping” from Australasia about 30 Ma, and radiating throughout the world during the Oligocene. This independently derived conclusion is consistent with the passerine fossil record.

Complete mitochondrial genome of the versicoloured emerald hummingbird Amazilia versicolor, a polymorphic species

The genome of the versicoloured emerald hummingbird (Amazilia versicolor) was partially sequenced in one-sixth of an Illumina HiSeq lane. The mitochondrial genome was assembled using MIRA and MITObim software, yielding a circular molecule of 16,861 bp in length and deposited in GenBank under the accession number KF624601. The mitogenome contained 13 protein-coding genes, 22 transfer tRNAs, 2 ribosomal RNAs and 1 non-coding control region. The molecule was assembled using 21,927 sequencing reads of 100 bp each, resulting in ∼130 × coverage of uniformly distributed reads along the genome. This is the forth mitochondrial genome described for this highly diverse family of birds and may benefit further phylogenetic, phylogeographic, population genetic and species delimitation studies of hummingbirds.

Seven complete mitochondrial genome sequences of bushtits (Passeriformes, Aegithalidae, Aegithalos): the evolution pattern in duplicated control regions

The control region (CR) of the mitochondrial DNA exhibits important functions in replication and transcription, and duplications of the CR have been reported in a wide range of animal groups. In most cases, concerted evolution is expected to explain the high similarity of duplicated CRs. In this paper, we present seven complete mitochondrial genome sequences from the bushtits (genus Aegithalos), in which we discovered two duplicated CRs, and try to survey the evolution pattern of these duplicated CRs. We also found that the duplicated CRs within one individual were almost identical, and variations were concentrated in two sections, one located between a poly-C site and a potential TAS (termination associated sequence) element, the other one located at the 3' end of the duplicated CRs. The phylogenetic analyses of paralogous CRs showed that the tree topology were depending on whether the two high variable regions at the upstream of TAS element and the 3'end of duplicated CRs: when they were concluded, the orthologous copies were closely related; when they were excluded, the paralogous copies in the same lineages were closely related. This may suggest the role of recombination in the evolution of duplicated CRs. Consequently, the recombination was detected, and the breakpoints were found at ∼120 bp (the upstream of the potential TAS element) and ∼1150 bp of the alignment of duplicated CRs. According to these results, we supposed that homologous recombination occurred between paralogous CRs from different mtDNA molecule was proposed as the most suitable mechanism for concerted evolution of the duplicated CRs, and the recombination took place in every replication cycle, so that most part of the duplicated regions remain identical within an individual, while the 5' and 3'end of the duplicated CRs were not involved in recombination, and evolved independently.

Eocene diversification of crown group rails (Aves: Gruiformes: Rallidae)

Central to our understanding of the timing of bird evolution is debate about an apparent conflict between fossil and molecular data. A deep age for higher level taxa within Neoaves is evident from molecular analyses but much remains to be learned about the age of diversification in modern bird families and their evolutionary ecology. In order to better understand the timing and pattern of diversification within the family Rallidae we used a relaxed molecular clock, fossil calibrations, and complete mitochondrial genomes from a range of rallid species analysed in a Bayesian framework. The estimated time of origin of Rallidae is Eocene, about 40.5 Mya, with evidence of intrafamiliar diversification from the Late Eocene to the Miocene. This timing is older than previously suggested for crown group Rallidae, but fossil calibrations, extent of taxon sampling and substantial sequence data give it credence. We note that fossils of Eocene age tentatively assigned to Rallidae are consistent with our findings. Compared to available studies of other bird lineages, the rail clade is old and supports an inference of deep ancestry of ground-dwelling habits among Neoaves.

Phylogeny and forelimb disparity in waterbirds

Previous work has shown that the relative proportions of wing components (i.e., humerus, ulna, carpometacarpus) in birds are related to function and ecology, but these have rarely been investigated in a phylogenetic context. Waterbirds including "Pelecaniformes," Ciconiiformes, Procellariiformes, Sphenisciformes, and Gaviiformes form a highly supported clade and developed a great diversity of wing forms and foraging ecologies. In this study, forelimb disparity in the waterbird clade was assessed in a phylogenetic context. Phylogenetic signal was assessed via Pagel's lambda, Blomberg's K, and permutation tests. We find that different waterbird clades are clearly separated based on forelimb component proportions, which are significantly correlated with phylogeny but not with flight style. Most of the traditional contents of "Pelecaniformes" (e.g., pelicans, cormorants, and boobies) cluster with Ciconiiformes (herons and storks) and occupy a reduced morphospace. These taxa are closely related phylogenetically but exhibit a wide range of ecologies and flight styles. Procellariiformes (e.g., petrels, albatross, and shearwaters) occupy a wide range of morphospace, characterized primarily by variation in the relative length of carpometacarpus and ulna. Gaviiformes (loons) surprisingly occupy a wing morphospace closest to diving petrels and penguins. Whether this result may reflect wing proportions plesiomorphic for the waterbird clade or a functional signal is unclear. A Bayesian approach detecting significant rate shifts across phylogeny recovered two such shifts. At the base of the two sister clades Sphenisciformes + Procellariiformes, a shift to an increase evolutionary rate of change is inferred for the ulna and carpometacarpus. Thus, changes in wing shape begin prior to the loss of flight in the wing-propelled diving clade. Several shifts to slower rate of change are recovered within stem penguins.

Phylogenetic informativeness reconciles ray-finned fish molecular divergence times

Background

Discordance among individual molecular age estimates, or between molecular age estimates and the fossil record, is observed in many clades across the Tree of Life. This discordance is attributed to a variety of variables including calibration age uncertainty, calibration placement, nucleotide substitution rate heterogeneity, or the specified molecular clock model. However, the impact of changes in phylogenetic informativeness of individual genes over time on phylogenetic inferences is rarely analyzed. Using nuclear and mitochondrial sequence data for ray-finned fishes (Actinopterygii) as an example, we extend the utility of phylogenetic informativeness profiles to predict the time intervals when nucleotide substitution saturation results in discordance among molecular ages estimated.

Results

We demonstrate that even with identical calibration regimes and molecular clock methods, mitochondrial based molecular age estimates are systematically older than those estimated from nuclear sequences. This discordance is most severe for highly nested nodes corresponding to more recent (i.e., Jurassic-Recent) divergences. By removing data deemed saturated, we reconcile the competing age estimates and highlight that the older mtDNA based ages were driven by nucleotide saturation.

Conclusions

Homoplasious site patterns in a DNA sequence alignment can systematically bias molecular divergence time estimates. Our study demonstrates that PI profiles can provide a non-arbitrary criterion for data exclusion to mitigate the influence of homoplasy on time calibrated branch length estimates. Analyses of actinopterygian molecular clocks demonstrate that scrutiny of the time scale on which sequence data is informative is a fundamental, but generally overlooked, step in molecular divergence time estimation.

Reweaving the tapestry: a supertree of birds

Our knowledge of the avian tree of life remains uncertain, particularly at deeper levels due to the rapid diversification early in their evolutionary history. They are the most abundant land vertebrate on the planet and have been of great historical interest to systematists. Birds are also economically and ecologically important and as a result are intensively studied, yet despite their importance and interest to humans around 13% of taxa currently on the endangered species list perhaps as a result of human activity. Despite all this no comprehensive phylogeny that includes both extinct and extant species currently exists. Here we present a species-level supertree, constructed using the Matrix Representation with Parsimony method, of Aves containing approximately two thirds of all species from nearly 1000 source phylogenies with a broad taxonomic coverage. The source data for the tree were collected and processed according to a strict protocol to ensure robust and accurate data handling. The resulting tree topology is largely consistent with molecular hypotheses of avian phylogeny. We identify areas that are in broad agreement with current views on avian systematics and also those that require further work. We also highlight the need for leaf-based support measures to enable the identification of rogue taxa in supertrees. This is a first attempt at a supertree of both extinct and extant birds, it is not intended to be utilised in an overhaul of avian systematics or as a basis for taxonomic re-classification but provides a strong basis on which to base further studies on macroevolution, conservation, biodiversity, comparative biology and character evolution, in particular the inclusion of fossils will allow the study of bird evolution and diversification throughout deep time.

Identifying localized biases in large datasets: a case study using the avian tree of life

Large-scale multi-locus studies have become common in molecular phylogenetics, with new studies continually adding to previous datasets in an effort to fully resolve the tree of life. Total evidence analyses that combine existing data with newly collected data are expected to increase the power of phylogenetic analyses to resolve difficult relationships. However, they might be subject to localized biases, with one or a few loci having a strong and potentially misleading influence upon the results. To examine this possibility we combined a newly collected 31-locus dataset that includes representatives of all major avian lineages with a published dataset of 19 loci that has a comparable number of sites (Hackett et al., 2008. Science 320, 1763-1768). This allowed us to explore the advantages of conducting total evidence analyses, and to determine whether it was also important to analyze new datasets independent of published ones. The total evidence analysis yielded results very similar to the published results, with only slightly increased support at a few nodes. However, analyzing the 31- and 19-locus datasets separately highlighted several differences. Two clades received strong support in the published dataset and total evidence analysis, but the support appeared to reflect bias at a single locus (β-fibrinogen [FGB]). The signal in FGB that supported these relationships was sufficient to result in their recovery with bootstrap support, even when combined with 49 loci lacking that signal. FGB did not appear to have a substantial impact upon the results of species tree methods, but another locus (brain-derived neurotrophic factor [BDNF]) did have an impact upon those analyses. These results demonstrated that localized biases can influence large-scale phylogenetic analyses but they also indicated that considering independent evidence and exploring multiple analytical approaches could reveal them.

Beyond phylogeny: pelecaniform and ciconiiform birds, and long-term niche stability

Phylogenetic trees are a starting point for the study of further evolutionary and ecological questions. We show that for avian evolutionary relationships, improved taxon sampling, longer sequences and additional data sets are giving stability to the prediction of the grouping of pelecaniforms and ciconiiforms, thus allowing inferences to be made about long-term niche occupancy. Here we report the phylogeny of the pelecaniform birds and their water-carnivore allies using complete mitochondrial genomes, and show that the basic groupings agree with nuclear sequence phylogenies, even though many short branches are not yet fully resolved. In detail, we show that the Pelecaniformes (minus the tropicbird) and the Ciconiiformes (storks, herons and ibises) form a natural group within a seabird water-carnivore clade. We find pelicans are the closest relatives of the shoebill (in a clade with the hammerkop), and we confirm that tropicbirds are not pelecaniforms. In general, the group appears to be an adaptive radiation into an 'aquatic carnivore' niche that it has occupied for 60-70 million years. From an ecological and life history perspective, the combined pelecaniform-ciconiform group is more informative than focusing on differences in morphology. These findings allow a start to integrating molecular evolution and macroecology.

A phylogeny of birds based on over 1,500 loci collected by target enrichment and high-throughput sequencing

Evolutionary relationships among birds in Neoaves, the clade comprising the vast majority of avian diversity, have vexed systematists due to the ancient, rapid radiation of numerous lineages. We applied a new phylogenomic approach to resolve relationships in Neoaves using target enrichment (sequence capture) and high-throughput sequencing of ultraconserved elements (UCEs) in avian genomes. We collected sequence data from UCE loci for 32 members of Neoaves and one outgroup (chicken) and analyzed data sets that differed in their amount of missing data. An alignment of 1,541 loci that allowed missing data was 87% complete and resulted in a highly resolved phylogeny with broad agreement between the Bayesian and maximum-likelihood (ML) trees. Although results from the 100% complete matrix of 416 UCE loci were similar, the Bayesian and ML trees differed to a greater extent in this analysis, suggesting that increasing from 416 to 1,541 loci led to increased stability and resolution of the tree. Novel results of our study include surprisingly close relationships between phenotypically divergent bird families, such as tropicbirds (Phaethontidae) and the sunbittern (Eurypygidae) as well as between bustards (Otididae) and turacos (Musophagidae). This phylogeny bolsters support for monophyletic waterbird and landbird clades and also strongly supports controversial results from previous studies, including the sister relationship between passerines and parrots and the non-monophyly of raptorial birds in the hawk and falcon families. Although significant challenges remain to fully resolving some of the deep relationships in Neoaves, especially among lineages outside the waterbirds and landbirds, this study suggests that increased data will yield an increasingly resolved avian phylogeny.

Ratite nonmonophyly: independent evidence from 40 novel Loci

Large-scale multilocus studies have become common in molecular phylogenetics, but the best way to interpret these studies when their results strongly conflict with prior information about phylogeny remains unclear. An example of such a conflict is provided by the ratites (the large flightless birds of southern land masses, including ostriches, emus, and rheas). Ratite monophyly is strongly supported by both morphological data and many earlier molecular studies and is used as a textbook example of vicariance biogeography. However, recent studies have indicated that ratites are not monophyletic; instead, the volant tinamous nest inside the ratites rather than forming their sister group within the avian superorder Palaeognathae. Large-scale studies can exhibit biases that reflect a number of factors, including limitations in the fit of the evolutionary models used for analyses and problems with sequence alignment, so the unexpected conclusion that ratites are not monophyletic needs to be rigorously evaluated. A rigorous approach to testing novel hypotheses generated by large-scale studies is to collect independent evidence (i.e., excluding the loci and/or traits used to generate the hypotheses). We used 40 nuclear loci not used in previous studies that investigated the relationship among ratites and tinamous. Our results strongly support the recent molecular studies, revealing that the deepest branch within Palaeognathae separates the ostrich from other members of the clade, rather than the traditional hypothesis that separates the tinamous from the ratites. To ensure these results reflected evolutionary history, we examined potential biases in types of loci used, heterotachy, alignment biases, and discordance between gene trees and the species tree. All analyses consistently supported nonmonophyly of the ratites and no confounding biases were identified. This confirmation that ratites are not monophyletic using independent evidence will hopefully stimulate further comparative research on paleognath development and genetics that might reveal the basis of the morphological convergence in these large, flightless birds.

Podargiform affinities of the enigmatic Fluvioviridavis platyrhamphus and the early diversification of Strisores ("Caprimulgiformes" + Apodiformes)

BACKGROUND

The early Eocene Green River Formation avifauna preserves exceptional exemplars of the earliest unambiguous stem representatives of many extant avian clades. We identify the basal-most member of Podargiformes (extant and fossil stem lineage frogmouths) based on a new specimen of Fluvioviridavis platyrhamphus, a unique neoavian bird from the Fossil Butte Member of the Green River Formation of Wyoming. Extant frogmouths (Podargidae) comprise approximately 13 nocturnal species with an exclusively Australasian distribution.

METHODOLOGY/PRINCIPAL FINDINGS

The new specimen was included in a combined phylogenetic analysis of morphological (osteology and soft tissue) and molecular sequence (cytochrome b, c-myc exon 3, and RAG) data sampling species-level taxa from both extant and extinct members of Steatornithidae, Podargidae, Caprimulgidae, Nyctibiidae, Aegothelidae, and Apodiformes ( = Strisores). New data from F. platyrhamphus help resolve phylogenetic relationships within Strisores, supporting placement of F. platyrhamphus and the European fossil form Masillapodargus longipes as basal parts of Podargiformes and also supporting a sister taxon relationship between Podargiformes and Steatornithiformes (oilbirds) within Strisores. This relationship is recovered only when fossil taxa are included, reaffirming the potential impact of stem fossil taxa on inferences of phylogenetic relationships. The well-preserved mandible and palate of the new specimen demonstrate that many of the unique characteristics of the skull that characterize the crown frogmouth clade Podargidae arose early in the evolutionary history of the clade, over 50 million years ago. Comparisons with the new specimen also indicate that Eurofluvioviridavis and Fluvioviridavis are not closely related.

CONCLUSIONS/SIGNIFICANCE

Together with the European fossil frogmouth Masillapodargus, Fluvioviridavis shows that Podargiformes had a much wider geographic distribution in the past, whereas extant species are restricted to Australasia. The Eocene record of Strisores from the Green River Formation and Messel Formation indicates most major subclade divergences had already occurred by the early-middle Eocene.

Testing hypotheses about the sister group of the passeriformes using an independent 30-locus data set

Although many phylogenetic studies have focused on developing hypotheses about relationships, advances in data collection and computation have increased the feasibility of collecting large independent data sets to rigorously test controversial hypotheses or carefully assess artifacts that may be misleading. One such relationship in need of independent evaluation is the position of Passeriformes (perching birds) in avian phylogeny. This order comprises more than half of all extant birds, and it includes one of the most important avian model systems (the zebra finch). Recent large-scale studies using morphology, mitochondrial, and nuclear sequence data have generated very different hypotheses about the sister group of Passeriformes, and all conflict with an older hypothesis generated using DNA-DNA hybridization. We used novel data from 30 nuclear loci, primarily introns, for 28 taxa to evaluate five major a priori hypotheses regarding the phylogenetic position of Passeriformes. Although previous studies have suggested that nuclear introns are ideal for the resolution of ancient avian relationships, introns have also been criticized because of the potential for alignment ambiguities and the loss of signal due to saturation. To examine these issues, we generated multiple alignments using several alignment programs, varying alignment parameters, and using guide trees that reflected the different a priori hypotheses. Although different alignments and analyses yielded slightly different results, our analyses excluded all but one of the five a priori hypotheses. In many cases, the passerines were sister to the Psittaciformes (parrots), and taxa were members of a larger clade that includes Falconidae (falcons) and Cariamidae (seriemas). However, the position of Coliiformes (mousebirds) was highly unstable in our analyses of 30 loci, and this represented the primary source of incongruence among analyses. Mousebirds were united with passerines or parrots in some analyses, suggesting an additional hypothesis that needs to be considered in future studies. There was no clear evidence that base-compositional convergence, saturation, or long-branch attraction affected our conclusions. These results provide independent evidence excluding four major hypotheses about the position of passerines, allowing the extensive studies on this group to be placed in a more rigorous evolutionary framework.

The evolutionary history of cockatoos (Aves: Psittaciformes: Cacatuidae)

Cockatoos are the distinctive family Cacatuidae, a major lineage of the order of parrots (Psittaciformes) and distributed throughout the Australasian region of the world. However, the evolutionary history of cockatoos is not well understood. We investigated the phylogeny of cockatoos based on three mitochondrial and three nuclear DNA genes obtained from 16 of 21 species of Cacatuidae. In addition, five novel mitochondrial genomes were used to estimate time of divergence and our estimates indicate Cacatuidae diverged from Psittacidae approximately 40.7 million years ago (95% CI 51.6-30.3 Ma) during the Eocene. Our data shows Cacatuidae began to diversify approximately 27.9 Ma (95% CI 38.1-18.3 Ma) during the Oligocene. The early to middle Miocene (20-10 Ma) was a significant period in the evolution of modern Australian environments and vegetation, in which a transformation from mainly mesic to xeric habitats (e.g., fire-adapted sclerophyll vegetation and grasslands) occurred. We hypothesize that this environmental transformation was a driving force behind the diversification of cockatoos. A detailed multi-locus molecular phylogeny enabled us to resolve the phylogenetic placements of the Palm Cockatoo (Probosciger aterrimus), Galah (Eolophus roseicapillus), Gang-gang Cockatoo (Callocephalon fimbriatum) and Cockatiel (Nymphicus hollandicus), which have historically been difficult to place within Cacatuidae. When the molecular evidence is analysed in concert with morphology, it is clear that many of the cockatoo species' diagnostic phenotypic traits such as plumage colour, body size, wing shape and bill morphology have evolved in parallel or convergently across lineages.

Evolution of modern birds revealed by mitogenomics: timing the radiation and origin of major orders

Mitochondrial (mt) genes and genomes are among the major sources of data for evolutionary studies in birds. This places mitogenomic studies in birds at the core of intense debates in avian evolutionary biology. Indeed, complete mt genomes are actively been used to unveil the phylogenetic relationships among major orders, whereas single genes (e.g., cytochrome c oxidase I [COX1]) are considered standard for species identification and defining species boundaries (DNA barcoding). In this investigation, we study the time of origin and evolutionary relationships among Neoaves orders using complete mt genomes. First, we were able to solve polytomies previously observed at the deep nodes of the Neoaves phylogeny by analyzing 80 mt genomes, including 17 new sequences reported in this investigation. As an example, we found evidence indicating that columbiforms and charadriforms are sister groups. Overall, our analyses indicate that by improving the taxonomic sampling, complete mt genomes can solve the evolutionary relationships among major bird groups. Second, we used our phylogenetic hypotheses to estimate the time of origin of major avian orders as a way to test if their diversification took place prior to the Cretaceous/Tertiary (K/T) boundary. Such timetrees were estimated using several molecular dating approaches and conservative calibration points. Whereas we found time estimates slightly younger than those reported by others, most of the major orders originated prior to the K/T boundary. Finally, we used our timetrees to estimate the rate of evolution of each mt gene. We found great variation on the mutation rates among mt genes and within different bird groups. COX1 was the gene with less variation among Neoaves orders and the one with the least amount of rate heterogeneity across lineages. Such findings support the choice of COX 1 among mt genes as target for developing DNA barcoding approaches in birds.

Full mitochondrial genome sequences of two endemic Philippine hornbill species (Aves: Bucerotidae) provide evidence for pervasive mitochondrial DNA recombination

Background

Although nowaday it is broadly accepted that mitochondrial DNA (mtDNA) may undergo recombination, the frequency of such recombination remains controversial. Its estimation is not straightforward, as recombination under homoplasmy (i.e., among identical mt genomes) is likely to be overlooked. In species with tandem duplications of large mtDNA fragments the detection of recombination can be facilitated, as it can lead to gene conversion among duplicates. Although the mechanisms for concerted evolution in mtDNA are not fully understood yet, recombination rates have been estimated from "one per speciation event" down to 850 years or even "during every replication cycle".

Results

Here we present the first complete mt genome of the avian family Bucerotidae, i.e., that of two Philippine hornbills, Aceros waldeni and Penelopides panini. The mt genomes are characterized by a tandemly duplicated region encompassing part of cytochrome b, 3 tRNAs, NADH6, and the control region. The duplicated fragments are identical to each other except for a short section in domain I and for the length of repeat motifs in domain III of the control region. Due to the heteroplasmy with regard to the number of these repeat motifs, there is some size variation in both genomes; with around 21,657 bp (A. waldeni) and 22,737 bp (P. panini), they significantly exceed the hitherto longest known avian mt genomes, that of the albatrosses. We discovered concerted evolution between the duplicated fragments within individuals. The existence of differences between individuals in coding genes as well as in the control region, which are maintained between duplicates, indicates that recombination apparently occurs frequently, i.e., in every generation.

Conclusions

The homogenised duplicates are interspersed by a short fragment which shows no sign of recombination. We hypothesize that this region corresponds to the so-called Replication Fork Barrier (RFB), which has been described from the chicken mitochondrial genome. As this RFB is supposed to halt replication, it offers a potential mechanistic explanation for frequent recombination in mitochondrial genomes.

Characterization of the complete mitochondrial genome of the Rock pigeon, Columba livia (Columbiformes: Columbidae)

The rock pigeon (Columba livia), or Rock dove, is a member of the bird family Columbidae. We mapped the complete mitochondrial genome of the Rock pigeon. The mitochondrial genome of this species is a circular molecule of 17,229 bp in length, encoding a standard set of 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes, plus a putative control region, demonstrating a structure very similar to that of other birds. As found in other vertebrates, most of these genes are coded on the H-strand, except for NADH dehydrogenase subunit 6 (nad6) and eight tRNA genes (Gln, Ala, Asn, Cys, Tyr, Ser(UCN), Pro, Glu). The AT skew and GC skew of the whole genome, protein-coding genes, tRNA, rRNA, and the control region were calculated for the complete mitochondrial genomes of 30 avian species, representing 29 orders. All protein-coding genes initiated with ATG, except for cox1 and nad5, which began with GTG. One extra nucleotide 'C' was present in NADH dehydrogenase subunit 3 (nad3). All tRNA gene sequences have the potential to fold into typical cloverleaf secondary structures. Within the control region, conserved sequences were identified in three domains. Although the conserved blocks, such as ETAS1, ETAS2, CSB1, CSB1-like, and boxes C, D, E, and F, are readily identifiable in the C. livia control region, the typical origin of H-strand replication (O(H)), CSB2 and CSB3 could not be detected. These results provide basic information for phylogenetic analyses of birds, especially Columbiformes species.

Evolution and cytological diversification of the green seaweeds (Ulvophyceae)

The Ulvophyceae, one of the four classes of the Chlorophyta, is of particular evolutionary interest because it features an unrivaled morphological and cytological diversity. Morphological types range from unicells and simple multicellular filaments to sheet-like and complex corticated thalli. Cytological layouts range from typical small cells containing a single nucleus and chloroplast to giant cells containing millions of nuclei and chloroplasts. In order to understand the evolution of these morphological and cytological types, the present paper aims to assess whether the Ulvophyceae are monophyletic and elucidate the ancient relationships among its orders. Our approach consists of phylogenetic analyses (maximum likelihood and Bayesian inference) of seven nuclear genes, small subunit nuclear ribosomal DNA and two plastid markers with carefully chosen partitioning strategies, and models of sequence evolution. We introduce a procedure for fast site removal (site stripping) targeted at improving phylogenetic signal in a particular epoch of interest and evaluate the specificity of fast site removal to retain signal about ancient relationships. From our phylogenetic analyses, we conclude that the ancestral ulvophyte likely was a unicellular uninucleate organism and that macroscopic growth was achieved independently in various lineages involving radically different mechanisms: either by evolving multicellularity with coupled mitosis and cytokinesis (Ulvales-Ulotrichales and Trentepohliales), by obtaining a multinucleate siphonocladous organization where every nucleus provides for its own cytoplasmic domain (Cladophorales and Blastophysa), or by developing a siphonous organization characterized by either one macronucleus or millions of small nuclei and cytoplasmic streaming (Bryopsidales and Dasycladales). We compare different evolutionary scenarios giving rise to siphonous and siphonocladous cytologies and argue that these did not necessarily evolve from a multicellular or even multinucleate state but instead could have evolved independently from a unicellular ancestor.

Complete mitochondrial genome of Otis tarda (Gruiformes: Otididae) and phylogeny of Gruiformes inferred from mitochondrial DNA sequences

The complete nucleotide sequence of mitochondrial genome of the Great bustard (Otis tarda) was determined by using polymerase chain reaction (PCR) method. The genome is 16,849 bp in size, containing 13 protein-coding, 2 ribosomal and 22 transfer RNA genes. Sequences of the tRNA genes can be folded into canonical cloverleaf secondary structure except for tRNA-Cys and tRNA-Ser (AGY), which lose "DHU" arm. Sequence analysis showed that the O. tarda mitochondrial control region (mtCR) contained many elements in common with other avian mtCRs. A microsatellite repeat was found in the 3'-peripheral domain of the O. tarda mtCR. Based on the mitochondrial DNA sequences of 12S rRNA, 16S rRNA and tRNA-Val, a phylogenetic study of Gruiformes was performed. The result showed that Otididae was a sister group to "core Gruiformes" and Charadriiformes with strong support (97% posterior probability values) in Bayesian analysis. The taxonomic status of Rhynochetidae, Mesitornithidae, Pedionomidae and Turnicidae that traditionally belonged to Gruiformes was also discussed in this paper.

The complete mitochondrial genomes of the whistling duck (Dendrocygna javanica) and black swan (Cygnus atratus): dating evolutionary divergence in Galloanserae

Galloanserae is an ancient and diverse avian group, for which comprehensive molecular evidence relevant to phylogenetic analysis in the context of molecular chronology is lacking. In this study, we present two additional mitochondrial genome sequences of Galloanserae (the whistling duck, Dendrocygna javanica, and the black swan, Cygnus atratus) to broaden the scope of molecular phylogenetic reconstruction. The lengths of the whistling duck's and black swan's mitochondrial genomes are 16,753 and 16,748 bases, respectively. Phylogenetic analyses suggest that Dendrocygna is more likely to be in a basal position of the branch consisting of Anatinae and Anserinae, an affiliation that does not conform to its traditional classification. Bayesian approaches were employed to provide a rough timescale for Galloanserae evolution. In general, a narrow range of 95% confidence intervals gave younger estimates than those based on limited genes and estimated that at least two lineages originated before the Coniacian epoch around 90 MYA, well before the Cretaceous-Tertiary boundary. In addition, these results, which were compatible with estimates from fossil evidence, also imply that the origin of numerous genera in Anseriformes took place in the late Oligocene to early Miocene. Taken together, the results presented here provide a working framework for future research on Galloanserae evolution, and they underline the utility of whole mitochondrial genome sequences for the resolution of deep divergence.

The erratic mitochondrial clock: variations of mutation rate, not population size, affect mtDNA diversity across birds and mammals

BACKGROUND

During the last ten years, major advances have been made in characterizing and understanding the evolution of mitochondrial DNA, the most popular marker of molecular biodiversity. Several important results were recently reported using mammals as model organisms, including (i) the absence of relationship between mitochondrial DNA diversity and life-history or ecological variables, (ii) the absence of prominent adaptive selection, contrary to what was found in invertebrates, and (iii) the unexpectedly large variation in neutral substitution rate among lineages, revealing a possible link with species maximal longevity. We propose to challenge these results thanks to the bird/mammal comparison. Direct estimates of population size are available in birds, and this group presents striking life-history trait differences with mammals (higher mass-specific metabolic rate and longevity). These properties make birds the ideal model to directly test for population size effects, and to discriminate between competing hypotheses about the causes of substitution rate variation.

RESULTS

A phylogenetic analysis of cytochrome b third-codon position confirms that the mitochondrial DNA mutation rate is quite variable in birds, passerines being the fastest evolving order. On average, mitochondrial DNA evolves slower in birds than in mammals of similar body size. This result is in agreement with the longevity hypothesis, and contradicts the hypothesis of a metabolic rate-dependent mutation rate. Birds show no footprint of adaptive selection on cytochrome b evolutionary patterns, but no link between direct estimates of population size and cytochrome b diversity. The mutation rate is the best predictor we have of within-species mitochondrial diversity in birds. It partly explains the differences in mitochondrial DNA diversity patterns observed between mammals and birds, previously interpreted as reflecting Hill-Robertson interferences with the W chromosome.

CONCLUSION

Mitochondrial DNA diversity patterns in birds are strongly influenced by the wide, unexpected variation of mutation rate across species. From a fundamental point of view, these results are strongly consistent with a relationship between species maximal longevity and mitochondrial mutation rate, in agreement with the mitochondrial theory of ageing. Form an applied point of view, this study reinforces and extends the message of caution previously expressed for mammals: mitochondrial data tell nothing about species population sizes, and strongly depart the molecular clock assumption.

Toward resolving deep neoaves phylogeny: data, signal enhancement, and priors

We report three developments toward resolving the challenge of the apparent basal polytomy of neoavian birds. First, we describe improved conditional down-weighting techniques to reduce noise relative to signal for deeper divergences and find increased agreement between data sets. Second, we present formulae for calculating the probabilities of finding predefined groupings in the optimal tree. Finally, we report a significant increase in data: nine new mitochondrial (mt) genomes (the dollarbird, New Zealand kingfisher, great potoo, Australian owlet-nightjar, white-tailed trogon, barn owl, a roadrunner [a ground cuckoo], New Zealand long-tailed cuckoo, and the peach-faced lovebird) and together they provide data for each of the six main groups of Neoaves proposed by Cracraft J (2001). We use his six main groups of modern birds as priors for evaluation of results. These include passerines, cuckoos, parrots, and three other groups termed "WoodKing" (woodpeckers/rollers/kingfishers), "SCA" (owls/potoos/owlet-nightjars/hummingbirds/swifts), and "Conglomerati." In general, the support is highly significant with just two exceptions, the owls move from the "SCA" group to the raptors, particularly accipitrids (buzzards/eagles) and the osprey, and the shorebirds may be an independent group from the rest of the "Conglomerati". Molecular dating mt genomes support a major diversification of at least 12 neoavian lineages in the Late Cretaceous. Our results form a basis for further testing with both nuclear-coding sequences and rare genomic changes.