A molecular timescale for galliform birds accounting for uncertainty in time estimates and heterogeneity of rates of DNA substitutions across lineages and sites

Mitochondrial and nuclear DNA sequences support a Cretaceous origin of Columbiformes and a dispersal-driven radiation in the Paleocene

Phylogenetic relationships among genera of pigeons and doves (Aves, Columbiformes) have not been fully resolved because of limited sampling of taxa and characters in previous studies. We therefore sequenced multiple nuclear and mitochondrial DNA genes totaling over 9000 bp from 33 of 41 genera plus 8 outgroup taxa, and, together with sequences from 5 other pigeon genera retrieved from GenBank, recovered a strong phylogenetic hypothesis for the Columbiformes. Three major clades were recovered with the combined data set, comprising the basally branching New World pigeons and allies (clade A) that are sister to Neotropical ground doves (clade B), and the Afro-Eurasian and Australasian taxa (clade C). None of these clades supports the monophyly of current families and subfamilies. The extinct, flightless dodo and solitaires (Raphidae) were embedded within pigeons and doves (Columbidae) in clade C, and monophyly of the subfamily Columbinae was refuted because the remaining subfamilies were nested within it. Divergence times estimated using a Bayesian framework suggest that Columbiformes diverged from outgroups such as Apodiformes and Caprimulgiformes in the Cretaceous before the mass extinction that marks the end of this period. Bayesian and maximum likelihood inferences of ancestral areas, accounting for phylogenetic uncertainty and divergence times, respectively, favor an ancient origin of Columbiformes in the Neotropical portion of what was then Gondwana. The radiation of modern genera of Columbiformes started in the Early Eocene to the Middle Miocene, as previously estimated for other avian groups such as ratites, tinamous, galliform birds, penguins, shorebirds, parrots, passerine birds, and toucans. Multiple dispersals of more derived Columbiformes between Australasian and Afro-Eurasian regions are required to explain current distributions.

DNA evidence for a Paleocene origin of the Alcidae (Aves: Charadriiformes) in the Pacific and multiple dispersals across northern oceans

The Alcidae is a group of marine, wing-propelled diving birds known as auks that are distributed along the coasts of the northern oceans. It has been suggested that auks originated in the Pacific coastal shores as early as the Miocene, and dispersed to the Atlantic either through the Arctic coasts of Eurasia and North America (northern dispersal route), or through upwelling zones in the coastal areas of California to Florida (southern dispersal route), before the closure of the Isthmus of Panama in the Pliocene. These hypotheses have not been tested formally because proposed phylogenies failed to recover fully bifurcating, well-supported phylogenetic relationships among and within genera. We therefore constructed a large data set of mitochondrial and nuclear DNA sequences for 21 of the 23 species of extant auks. We also included sequences from two other extant and one extinct species retrieved from GenBank. Our analyses recovered a well-supported phylogenetic hypothesis among and within genera. Aethia is the only genus for which we could not obtain strong support for species relationships, probably due to incomplete lineage sorting. By applying a Bayesian method of molecular dating that allows for rate variation across lineages and genes, we showed that auks became an independent lineage in the Early Paleocene and radiated gradually from the Early Eocene to the Quaternary. Reconstruction of ancestral areas strongly suggests that auks originated in the Pacific during the Paleocene. The southern dispersal route seems to have favored the subsequent colonization of the northern Atlantic Ocean during the Eocene and Oligocene. The northern route across the Arctic Ocean was probably only used more recently after the opening of the Norwegian Sea in the Middle Miocene and the opening of the Bering Strait in the Late Miocene. We postulate that the ancestors of auks lived in a warmer world than that currently occupied by auks, and became gradually adapted to feeding in cool marine currents with high biomass productivity. Hence, warmer tropical waters are now a barrier for the dispersal of auks into the Southern Hemisphere, as it is for penguins in the opposite direction.

Waves of genomic hitchhikers shed light on the evolution of gamebirds (Aves: Galliformes)

BACKGROUND

The phylogenetic tree of Galliformes (gamebirds, including megapodes, currassows, guinea fowl, New and Old World quails, chicken, pheasants, grouse, and turkeys) has been considerably remodeled over the last decades as new data and analytical methods became available. Analyzing presence/absence patterns of retroposed elements avoids the problems of homoplastic characters inherent in other methodologies. In gamebirds, chicken repeats 1 (CR1) are the most prevalent retroposed elements, but little is known about the activity of their various subtypes over time. Ascertaining the fixation patterns of CR1 elements would help unravel the phylogeny of gamebirds and other poorly resolved avian clades.

RESULTS

We analyzed 1,978 nested CR1 elements and developed a multidimensional approach taking advantage of their transposition in transposition character (TinT) to characterize the fixation patterns of all 22 known chicken CR1 subtypes. The presence/absence patterns of those elements that were active at different periods of gamebird evolution provided evidence for a clade (Cracidae + (Numididae + (Odontophoridae + Phasianidae))) not including Megapodiidae; and for Rollulus as the sister taxon of the other analyzed Phasianidae. Genomic trace sequences of the turkey genome further demonstrated that the endangered African Congo Peafowl (Afropavo congensis) is the sister taxon of the Asian Peafowl (Pavo), rejecting other predominantly morphology-based groupings, and that phasianids are monophyletic, including the sister taxa Tetraoninae and Meleagridinae.

CONCLUSION

The TinT information concerning relative fixation times of CR1 subtypes enabled us to efficiently investigate gamebird phylogeny and to reconstruct an unambiguous tree topology. This method should provide a useful tool for investigations in other taxonomic groups as well.

Molecular evolutionary genomics of birds

Insight into the molecular evolution of birds has been offered by the steady accumulation of avian DNA sequence data, recently culminating in the first draft sequence of an avian genome, that of chicken. By studying avian molecular evolution we can learn about adaptations and phenotypic evolution in birds, and also gain an understanding of the similarities and differences between mammalian and avian genomes. In both these lineages, there is pronounced isochore structure with highly variable GC content. However, while mammalian isochores are decaying, they are maintained in the chicken lineage, which is consistent with a biased gene conversion model where the high and variable recombination rate of birds reinforces heterogeneity in GC. In Galliformes, GC is positively correlated with the rate of nucleotide substitution; the mean neutral mutation rate is 0.12-0.15% at each site per million years but this estimate comes with significant local variation in the rate of mutation. Comparative genomics reveals lower d(N)/d(S) ratios on micro- compared to macrochromosomes, possibly due to population genetic effects or a non-random distribution of genes with respect to chromosome size. A non-random genomic distribution is shown by genes with sex-biased expression, with male-biased genes over-represented and female-biased genes under-represented on the Z chromosome. A strong effect of selection is evident on the non-recombining W chromosome with high d(N)/d(S) ratios and limited polymorphism. Nucleotide diversity in chicken is estimated at 4-5 x 10(-3) which might be seen as surprisingly high given presumed bottlenecks during domestication, but is lower than that recently observed in several natural populations of other species. Several important aspects of the molecular evolutionary process of birds remain to be understood and it can be anticipated that the upcoming genome sequence of a second bird species, the zebra finch, as well as the integration of data on gene expression, shall further advance our knowledge of avian evolution.

The genomic landscape of short insertion and deletion polymorphisms in the chicken (Gallus gallus) Genome: a high frequency of deletions in tandem duplicates

It is increasingly recognized that insertions and deletions (indels) are an important source of genetic as well as phenotypic divergence and diversity. We analyzed length polymorphisms identified through partial (0.25x) shotgun sequencing of three breeds of domestic chicken made by the International Chicken Polymorphism Map Consortium. A data set of 140,484 short indel polymorphisms in unique DNA was identified after filtering for microsatellite structures. There was a significant excess of tandem duplicates at indel sites, with deletions of a duplicate motif outnumbering the generation of duplicates through insertion. Indel density was lower in microchromosomes than in macrochromosomes, in the Z chromosome than in autosomes, and in 100 bp of upstream sequence, 5'-UTR, and first introns than in intergenic DNA and in other introns. Indel density was highly correlated with single nucleotide polymorphism (SNP) density. The mean density of indels in pairwise sequence comparisons was 1.9 x 10(-4) indel events/bp, approximately 5% the density of SNPs segregating in the chicken genome. The great majority of indels involved a limited number of nucleotides (median 1 bp), with A-rich motifs being overrepresented at indel sites. The overrepresentation of deletions at tandem duplicates indicates that replication slippage in duplicate sequences is a common mechanism behind indel mutation. The correlation between indel and SNP density indicates common effects of mutation and/or selection on the occurrence of indels and point mutations.

Phylogenetic Position of the New World Quail (Odontophoridae): Eight Nuclear Loci and Three Mitochondrial Regions Contradict Morphology and the Sibley-Ahlquist Tapestry

The evolutionary relationship between the New World quail (Odontophoridae) and other groups of Galliformes has been an area of debate. In particular, the relationship between the New World quail and guineafowl (Numidinae) has been difficult to resolve. We analyzed >8 kb of DNA sequence data from 16 taxa that represent all major lineages of Galliformes to resolve the phylogenetic position of New World quail. A combined data set of eight nuclear loci and three mitochondrial regions analyzed with maximum parsimony, maximum likelihood, and Bayesian methods provide congruent and strong support for New World quail being basal members of a phasianid clade that excludes guineafowl. By contrast, the three mitochondrial regions exhibit modest incongruence with each other. This is reflected in the combined mitochondrial analyses that weakly support the Sibley-Ahlquist topology that placed the New World quail basal in relation to guineafowl and led to the placement of New World quail in its own family, sister to the Phasianidae. However, simulation-based topology tests using the mitochondrial data were unable to reject the topology suggested by our combined (mitochondrial and nuclear) data set. By contrast, similar tests using our most likely topology and our combined nuclear and mitochondrial data allow us to strongly reject the Sibley-Ahlquist topology and a topology based on morphological data that unites Old and New World quail.

Posición Filogenética de las Codornices del Nuevo Mundo (Odontophoridae): Ocho Loci Nucleares y Tres Regiones Mitocondriales Contradicen la Morfología y la Filogenia de Sibley y Ahlquist

Higher-order phylogeny of modern birds (Theropoda, Aves: Neornithes) based on comparative anatomy. II. Analysis and discussion

In recent years, avian systematics has been characterized by a diminished reliance on morphological cladistics of modern taxa, intensive palaeornithogical research stimulated by new discoveries and an inundation by analyses based on DNA sequences. Unfortunately, in contrast to significant insights into basal origins, the broad picture of neornithine phylogeny remains largely unresolved. Morphological studies have emphasized characters of use in palaeontological contexts. Molecular studies, following disillusionment with the pioneering, but non-cladistic, work of Sibley and Ahlquist, have differed markedly from each other and from morphological works in both methods and findings. Consequently, at the turn of the millennium, points of robust agreement among schools concerning higher-order neornithine phylogeny have been limited to the two basalmost and several mid-level, primary groups. This paper describes a phylogenetic (cladistic) analysis of 150 taxa of Neornithes, including exemplars from all non-passeriform families, and subordinal representatives of Passeriformes. Thirty-five outgroup taxa encompassing Crocodylia, predominately theropod Dinosauria, and selected Mesozoic birds were used to root the trees. Based on study of specimens and the literature, 2954 morphological characters were defined; these characters have been described in a companion work, approximately one-third of which were multistate (i.e. comprised at least three states), and states within more than one-half of these multistate characters were ordered for analysis. Complete heuristic searches using 10 000 random-addition replicates recovered a total solution set of 97 well-resolved, most-parsimonious trees (MPTs). The set of MPTs was confirmed by an expanded heuristic search based on 10 000 random-addition replicates and a full ratchet-augmented exploration to ascertain global optima. A strict consensus tree of MPTs included only six trichotomies, i.e. nodes differing topologically among MPTs. Bootstrapping (based on 10 000 replicates) percentages and ratchet-minimized support (Bremer) indices indicated most nodes to be robust. Several fossil Neornithes (e.g. Dinornithiformes, Aepyornithiformes) were placed within the ingroup a posteriori either through unconstrained, heursitic searches based on the complete matrix augmented by these taxa separately or using backbone-constraints. Analysis confirmed the topology among outgroup Theropoda and achieved robust resolution at virtually all levels of the Neornithes. Findings included monophyly of the palaeognathous birds, comprising the sister taxa Tinamiformes and ratites, respectively, and the Anseriformes and Galliformes as monophyletic sister-groups, together forming the sister-group to other Neornithes exclusive of the Palaeognathae (Neoaves). Noteworthy inferences include: (i) the sister-group to remaining Neoaves comprises a diversity of marine and wading birds; (ii) Podicipedidae are the sister-group of Gaviidae, and not closely related to the Phoenicopteridae, as recently suggested; (iii) the traditional Pelecaniformes, including the shoebill (Balaeniceps rex) as sister-taxon to other members, are monophyletic; (iv) traditional Ciconiiformes are monophyletic; (v) Strigiformes and Falconiformes are sister-groups; (vi) Cathartidae is the sister-group of the remaining Falconiformes; (vii) Ralliformes (Rallidae and Heliornithidae) are the sister-group to the monophyletic Charadriiformes, with the traditionally composed Gruiformes and Turniciformes (Turnicidae and Mesitornithidae) sequentially paraphyletic to the entire foregoing clade; (viii) Opisthocomus hoazin is the sister-taxon to the Cuculiformes (including the Musophagidae); (ix) traditional Caprimulgiformes are monophyletic and the sister-group of the Apodiformes; (x) Trogoniformes are the sister-group of Coliiformes; (xi) Coraciiformes, Piciformes and Passeriformes are mutually monophyletic and closely related; and (xii) the Galbulae are retained within the Piciformes. Unresolved portions of the Neornithes (nodes having more than one most-parsimonious solution) comprised three parts of the tree: (a) several interfamilial nodes within the Charadriiformes; (b) a trichotomy comprising the (i) Psittaciformes, (ii) Columbiformes and (iii) Trogonomorphae (Trogoniformes, Coliiformes) + Passerimorphae (Coraciiformes, Piciformes, Passeriformes); and (c) a trichotomy comprising the Coraciiformes, Piciformes and Passeriformes. The remaining polytomies were among outgroups, although several of the highest-order nodes were only marginally supported; however, the majority of nodes were resolved and met or surpassed conventional standards of support. Quantitative comparisons with alternative hypotheses, examination of highly supportive and diagnostic characters for higher taxa, correspondences with prior studies, complementarity and philosophical differences with palaeontological phylogenetics, promises and challenges of palaeogeography and calibration of evolutionary rates of birds, and classes of promising evidence and future directions of study are reviewed. Homology, as applied to avian examples of apparent homologues, is considered in terms of recent theory, and a revised annotated classification of higher-order taxa of Neornithes and other closely related Theropoda is proposed. (c) 2007 The Linnean Society of London, Zoological Journal of the Linnean Society, 2007, 149, 1-95.

Phylogenetic relationships and historical biogeography of neotropical parrots (Psittaciformes: Psittacidae: Arini) inferred from mitochondrial and nuclear DNA sequences

Previous hypotheses of phylogenetic relationships among Neotropical parrots were based on limited taxon sampling and lacked support for most internal nodes. In this study we increased the number of taxa (29 species belonging to 25 of the 30 genera) and gene sequences (6388 base pairs of RAG-1, cyt b, NADH2, ATPase 6, ATPase 8, COIII, 12S rDNA, and 16S rDNA) to obtain a stronger molecular phylogenetic hypothesis for this group of birds. Analyses of the combined gene sequences using maximum likelihood and Bayesian methods resulted in a well-supported phylogeny and indicated that amazons and allies are a sister clade to macaws, conures, and relatives, and these two clades are in turn a sister group to parrotlets. Key morphological and behavioral characters used in previous classifications were mapped on the molecular tree and were phylogenetically uninformative. We estimated divergence times of taxa using the molecular tree and Bayesian and penalized likelihood methods that allow for rate variation in DNA substitutions among sites and taxa. Our estimates suggest that the Neotropical parrots shared a common ancestor with Australian parrots 59 Mya (million of years ago; 95% credibility interval (CrI) 66, 51 Mya), well before Australia separated from Antarctica and South America, implying that ancestral parrots were widespread in Gondwanaland. Thus, the divergence of Australian and Neotropical parrots could be attributed to vicariance. The three major clades of Neotropical parrots originated about 50 Mya (95% CrI 57, 41 Mya), coinciding with periods of higher sea level when both Antarctica and South America were fragmented with transcontinental seaways, and likely isolated the ancestors of modern Neotropical parrots in different regions in these continents. The correspondence between major paleoenvironmental changes in South America and the diversification of genera in the clade of amazons and allies between 46 and 16 Mya suggests they diversified exclusively in South America. Conversely, ancestors of parrotlets and of macaws, conures, and allies may have been isolated in Antarctica and/or the southern cone of South America, and only dispersed out of these southern regions when climate cooled and Antarctica became ice-encrusted about 35 Mya. The subsequent radiation of macaws and their allies in South America beginning about 28 Mya (95% CrI 22, 35 Mya) coincides with the uplift of the Andes and the subsequent formation of dry, open grassland habitats that would have facilitated ecological speciation via niche expansion from forested habitats.

A Mitogenomic Timescale for Birds Detects Variable Phylogenetic Rates of Molecular Evolution and Refutes the Standard Molecular Clock

Current understanding of the diversification of birds is hindered by their incomplete fossil record and uncertainty in phylogenetic relationships and phylogenetic rates of molecular evolution. Here we performed the first comprehensive analysis of mitogenomic data of 48 vertebrates, including 35 birds, to derive a Bayesian timescale for avian evolution and to estimate rates of DNA evolution. Our approach used multiple fossil time constraints scattered throughout the phylogenetic tree and accounts for uncertainties in time constraints, branch lengths, and heterogeneity of rates of DNA evolution. We estimated that the major vertebrate lineages originated in the Permian; the 95% credible intervals of our estimated ages of the origin of archosaurs (258 MYA), the amniote-amphibian split (356 MYA), and the archosaur-lizard divergence (278 MYA) bracket estimates from the fossil record. The origin of modern orders of birds was estimated to have occurred throughout the Cretaceous beginning about 139 MYA, arguing against a cataclysmic extinction of lineages at the Cretaceous/Tertiary boundary. We identified fossils that are useful as time constraints within vertebrates. Our timescale reveals that rates of molecular evolution vary across genes and among taxa through time, thereby refuting the widely used mitogenomic or cytochrome b molecular clock in birds. Moreover, the 5-Myr divergence time assumed between 2 genera of geese (Branta and Anser) to originally calibrate the standard mitochondrial clock rate of 0.01 substitutions per site per lineage per Myr (s/s/l/Myr) in birds was shown to be underestimated by about 9.5 Myr. Phylogenetic rates in birds vary between 0.0009 and 0.012 s/s/l/Myr, indicating that many phylogenetic splits among avian taxa also have been underestimated and need to be revised. We found no support for the hypothesis that the molecular clock in birds "ticks" according to a constant rate of substitution per unit of mass-specific metabolic energy rather than per unit of time, as recently suggested. Our analysis advances knowledge of rates of DNA evolution across birds and other vertebrates and will, therefore, aid comparative biology studies that seek to infer the origin and timing of major adaptive shifts in vertebrates.