Do birds of a feather flock together?

bird

Editor's Introduction

Ancient DNA reveals elephant birds and kiwi are sister taxa and clarifies ratite bird evolution

The species we see on Earth today have been shaped by millions of years of evolution. Over great periods of time, species often change and adapt in response to environmental cues and constraints. One of the great mysteries of science has been to determine the origins of the species alive today. The history of a group of large flightless birds called ratites has puzzled evolutionary biologists for years. In the current study, scientists use DNA recovered from fossils of extinct elephant birds, an ancient ratite, to clarify the mysterious origin of these flightless birds.

Paper Details

Original title
Ancient DNA reveals elephant birds and kiwi are sister taxa and clarifies ratite bird evolution
Authors
Kieren Mitchell
Original publication date
Reference
Vol. 344 no. 6186 pp. 898-900
Issue name
Science
DOI
DOI: 10.1126/science.1251981

Abstract

The evolution of the ratite birds has been widely attributed to vicariant speciation, driven by the Cretaceous breakup of the supercontinent Gondwana. The early isolation of Africa and Madagascar implies that the ostrich and extinct Madagascan elephant birds (Aepyornithidae) should be the oldest ratite lineages. We sequenced the mitochondrial genomes of two elephant birds and performed phylogenetic analyses, which revealed that these birds are the closest relatives of the New Zealand kiwi and are distant from the basal ratite lineage of ostriches. This unexpected result strongly contradicts continental vicariance and instead supports flighted dispersal in all major ratite lineages. We suggest that convergence toward gigantism and flightlessness was facilitated by early Tertiary expansion into the diurnal herbivory niche after the extinction of the dinosaurs.

Report

Despite extensive studies, the evolutionary history of the giant flightless ratite birds of the Southern Hemisphere landmasses and the related flighted tinamous of South America has remained a major unresolved question. The ratites and tinamous, termed “palaeognaths” due to their shared basal palate structure, form the sister taxon to all other living birds (neognaths). The living ratites are one of the few bird groups composed largely of giant terrestrial herbivores and include: the emu and cassowary in Australia and New Guinea, the kiwi in New Zealand, the ostrich in Africa, and the rhea in South America. In addition, two recently extinct groups included the largest birds known: the moa from New Zealand (height up to 2 to 3 m, 250 kg in weight) (1) and elephant birds from Madagascar (2 to 3 m in height, up to 275 kg in weight) (23). Ratites have been believed to have originated through vicariant speciation driven by the continental breakup of the supercontinent Gondwana on the basis of congruence between the sequence of continental rifting and the presumed order of lineage divergence and distribution of ratites (45).

New Zealand is the only landmass to have supported two major ratite lineages: the giant herbivorous moa and the chicken-sized, nocturnal, omnivorous kiwi. Morphological phylogenetic analyses initially suggested that these two groups were each other’s closest relatives (67), presumably diverging after the isolation of an ancestral form following the separation of New Zealand and Australia in the late Cretaceous ~80 to 60 million years ago (Ma) (8). However, subsequent studies suggest that kiwi are more closely related to the Australasian emu and cassowaries (910), whereas the closest living relatives of the giant moa are the flighted South American tinamous (1114). The latter relationship was completely unexpected on morphological grounds and suggests a more complex evolutionary history than predicted by a model of strict vicariant speciation. By rendering ratites paraphyletic, the relationship between the moa and tinamous also strongly suggests that gigantism and flightlessness have evolved multiple times among palaeognaths (1213).

Perhaps the most enigmatic of the modern palaeognaths are the recently extinct giant Madagascan elephant birds. Africa and Madagascar were the first continental fragments to rift from the supercontinent Gondwana, separating from the other continents (and each other) completely during the Early Cretaceous (~130 to 100 Ma) (15). Consequently, the continental vicariance model predicts that elephant birds and ostriches should be the basal palaeognath lineages (16). Most molecular analyses recover the ostrich in a basal position, consistent with a vicariant model. However, the phylogenetic position of the elephant birds remains unresolved, as cladistic studies of ratite morphology are sensitive to character choice and may be confounded by convergence (17), whereas DNA studies have been hampered by the generally poor molecular preservation of elephant bird remains (18).

We used hybridization enrichment with in-solution RNA arrays of palaeognath mitochondrial genome sequences and high-throughput sequencing to sequence near-complete mitochondrial genomes from both elephant bird genera: Aepyornis and Mullerornis. Phylogenetic analyses placed the two taxa, Aepyornis hildebrandti (15,547 base pairs) and Mullerornis agilis (15,731 base pairs), unequivocally as the sister taxa to the kiwi (Fig. 1 and fig. S1). This result was consistently retrieved, regardless of phylogenetic method or taxon sampling, and was strongly supported by topological tests (19). To our knowledge, no previous study has suggested this relationship, probably because of the disparate morphology, ecology, and distribution of the two groups. Elephant birds were herbivorous, almost certainly diurnal, and among the largest birds known, whereas kiwi are highly derived omnivores, nocturnal, and about two orders of magnitude smaller. Elephant birds more closely resemble the moa, and analyses of morphology have suggested a close relationship between these taxa (17). However, adding morphological characters to our molecular data set increased support for the relationship between elephant birds and kiwi (figs. S2 and S3) and allowed for identification of several distinctive character states that diagnose this clade [see list in (19)].

140523-b-evo-f1.jpg

Fig 1.  Phylogenetic position of the elephant birds from mitochondrial sequence data. Bayesian posterior probabilities and maximum likelihood bootstrap are presented in black below each branch; asterisks denote branches that received maximum possible support (bootstrap = 100%, Bayesian posterior probability = 1.0). Divergence dates [blue numbers above branches; blue bars represent 95% highest posterior density (HPD) intervals] were inferred with six well-supported node age constraints (table S5). Blue arrows mark the minimum date for the evolution of flightlessness in lineages for which fossil evidence is available (2122). The scale is given in millions of years before the present. Silhouettes indicate the relative size of representative taxa. Species diversity for each major clade is presented in parentheses, with extinct groups shown in red. The dagger symbol (†) indicates that the number of elephant bird species is uncertain.

Obtaining mitochondrial DNA sequences

Before the authors could infer information about the phylogenetic history of the ratites, they had to sequence mitochondrial DNA from their two elephant bird fossils.

To separate the mitochondrial DNA from the nuclear DNA in their samples, they used a technology called hybridization enrichment. Basically, this technique uses RNA “baits” to “fish” out only the mitochondrial DNA while excluding other types of DNA. 

The link below gives a more detailed description of how the technology works.

http://www.mycroarray.com/mybaits/MYbaits+sequence+capture+target+enrichment+kit.html

DNA Sequencing

To sequence the mitochondrial DNA from two elephant bird fossils, the authors used a technology called semiconductor sequencing using the Ion Torrent system from Life Technologies.  Essentially, this technique uses a semiconductor chip to read changes in pH to sequence DNA.

Developing the phylogenetic tree

After comparing their mitochondrial DNA sequences to those of several other palaeognaths, the authors could infer a phylogenetic tree based on DNA sequence similarity.

The branches give us information about how each of the different groups of species, or clades, evolved from one common ancestor. 

Based on sequence similarity, elephant birds and kiwis are each other’s closest relatives on the tree.

How long ago did these species diverge?

From the sequence information, the authors can also estimate how long ago each of the branches on the phylogenetic tree diverged from the common ancestor.  

The blue numbers in this figure represent how many millions of years ago each branch split into two.  For example, the kiwi and the elephant bird appear to have diverged from a common ancestor about 50 million years ago.  

The blue arrows give us an estimate of how long ago the trait of flightlessness evolved for each branch.  For example, the authors estimate that flightlessness for the ostrich evolved approximately 20 million years ago.   

Speciation by continental vicariance provides a poor explanation of the close relationship between elephant birds and kiwi. Madagascar and New Zealand have never been directly connected, and molecular dates calculated from the genetic data suggest that kiwi and elephant birds diverged after the breakup of Gondwana (Fig. 1 and fig. S4). However, mean node age estimates among palaeognath lineages are sensitive to taxon sampling (Fig. 2), so molecular dating provides limited power for testing hypotheses about ratite biogeography. Depending on taxon sampling, estimates for the basal divergence among palaeognaths are equally consistent with the separation of Africa ~100 Ma (15) and the Cretaceous-Tertiary (KPg) boundary (~65 Ma) (Fig. 2). Thus, topological comparisons may be a more robust tool to test hypotheses of vicariance and connection.

140523-b-evo-f2.jpg

Fig 2. Sensitivity of palaeognath age estimation to taxon sampling and genetic loci used. Mean and 95% HPD intervals are displayed for the age (basal divergence) of crown palaeognaths, as inferred under several data set permutations. The y axis represents time before the present in millions of years, whereas age estimates for individual data sets are arrayed on the x axis for both mitochondrial (circles) and nuclear loci (triangles). Results are presented for the following taxon sets: A, a full taxon set including ratites, tinamous, neognaths, and crocodilians; B, ratites, tinamous, and neognaths only; C, ratites, neognaths, and crocodilians only; and D, ratites and neognaths only. Taxon sets A and C are each calibrated with six fossil node constraints, whereas taxon sets B and D are calibrated with a subset of four relevant constraints (table S5). Taxon sets are represented visually with silhouettes of an ostrich (ratites), a flying tinamou (tinamous), a duck (neognaths), and an alligator (crocodilians). Analyses excluding the rate-anomalous tinamous (taxon sets C and D) retrieve a young age near the KPg boundary.

What is molecular dating?

Molecular dating uses differences in DNA sequence from one species to another to estimate how long ago the species diverged.  

This technique assumes that genetic changes, or mutations, accumulate over time at roughly the same rate.  So, the greater the difference in genetic sequence, the farther back the divergence date.  

The authors used this technique with their mitochondrial DNA sequences to estimate how long ago different ratite taxa diverged.

How does taxon sampling influence divergence date estimations?

Different studies of palaeognath phylogeny have reported very different estimations of divergence date.  But why is that?

To find out, the authors chose four different groups of taxa, each one containing one palaeognath, and used molecular dating to estimate the time of divergence of the palaeognath from the other taxa.

Group A had the most diversity while group D had the least diversity.

As this figure shows, each of the four taxon groups have drastically different estimations of divergence date.      

How does the data source influence divergence data estimation?

The authors also compared divergence date estimations based on nuclear DNA sequence versus mitochondrial DNA sequence for each of their four taxon groups.  

Taxon Sets B and D had roughly the same date estimation regardless of sequence source.  However, taxon sets A and C had considerably different date estimations depending on the source of the DNA sequence.

A cautionary tale 

This analysis of taxon sampling and DNA sequence source emphasizes the importance of prudence when selecting samples for any scientific study.  

It is crucial that the sample is truly representative of the entire group.  Unfortunately, when dealing with ancient phylogenies, it is not always possible to predict how well a specific sample will represent a group.  

The phylogenetic placement of the elephant bird as sister to the kiwi creates a marked discordance between the order of continental breakup (Fig. 3, A and B) and the sequence of palaeognath divergences (Fig. 3C). Instead, it appears that the common ancestor of elephant birds and kiwis was probably flighted and capable of long-distance dispersal, which is supported by a small, possibly flighted kiwi relative from the Early Miocene of New Zealand (20). Together, the phylogenetic position of the flighted tinamous and apparent flighted ancestor of the kiwi and elephant bird imply that every major ratite lineage independently lost flight (Fig. 1). We suggest that flighted dispersal was the primary driver of the distribution of palaeognath lineages and that the discordance between distribution and phylogeny is more consistent with lineage turnover in a phylogenetically diverse, flighted, and widespread clade. Early Tertiary palaeognaths were capable of long-distance dispersal, with remains found well outside the range of modern ratites, including the flighted lithornithids in North America and Europe and the flightless Palaeotis and Remiornis in Europe (Fig. 3A) (21). Rapid diversification through flighted dispersal also provides an explanation for the short and often poorly supported internodes amongst basal extant ratite lineages (1314).

140523-b-evo-f3.jpg

Fig 3. Conflict between inferred palaeognath phylogeny and the topology predicted by continental vicariance.(A) Relative position of continents during the Late Cretaceous and Tertiary. Continental landmasses are colored according to order of severance from the remaining Gondwanan landmass: Africa and Madagascar first (dark gray; 100 to 130 Ma), followed by New Zealand (red; 60 to 80 Ma), then finally Australia, Antarctica, and South America (green; 30 to 50 Ma). Palaeognath-bearing fossil localities from the late Palaeocene and Eocene (2122) are represented by circles (flighted taxa) and triangles (flightless taxa). (B) Predicted phylogeny of ratites under a model of speciation governed solely by continental vicariance. (C) Palaeognath phylogeny as inferred in the present study (see Fig. 1).

The Breakup of Gondwana

Gondwana is a one of the two landmasses formed after the breakup of supercontinent Pangaea.  

The sequence of the breakup of this landmass greatly influenced the ratite phylogeny developed before the completion of the current study.

https://www.bbc.com/news/av/science-environment-16945234/animation-shows-rise-and-fall-of-last-supercontinent

Panel A

This panel summarizes the progression of the breakup of Gondwana, beginning 160 million years ago.  

The authors have colored each landmass based on the order in which they separated from the supercontinent.  

As the continents broke apart and became separated by bodies of water, the palaeognath species became farther and farther apart from one another.  

The bottom map in the panel indicates the location of several palaeognath fossil remains, showing that the species were widespread after the break-up of the super continent. 

Panel B

Panels B and C both depict palaeognath phylogenies.  Each species silhouette represents a different taxon.  They are colored to correspond in Panel A to the continent on which they thrived. 

The tree in this panel is the inferred phylogeny based on the continental vicariance model of speciation.  

This means that the drawn relationships were based on the order of the breakup of Gondwana, summarized in Panel A.  

Beginning from the top, the taxa are tinamou, rhea, emu, moa, kiwi, elephant bird, and ostrich.   

Panel C

This phylogeny was inferred based on the data collected in this study.  

The taxa were ordered based on mitochondrial DNA sequence similarity and not by proximity to one another after continental breakup.  

From top to bottom, the taxa are moa, tinamou, elephant bird, kiwi, emu, rhea, and ostrich.

Early ratite evolution appears to have been dominated by flighted dispersal and parallel evolution, with flightlessness evolving a minimum of six times and gigantism a minimum of five (Fig. 1) (22), suggesting that adaptations for cursoriality may have confounded phylogenetic inference. Elsewhere, avian gigantism and flightlessness are almost exclusively observed in island environments in the absence of mammalian predators and competitors (e.g., the dodo). However, each of the landmasses occupied by ratites (excluding New Zealand) is now home to a diverse mammalian fauna. We suggest that the initial evolution of flightless ratites began in the ecological vacuum after the KPg mass-extinction event and the extinction of the dinosaurs (1221). Most mammals appear to have remained relatively small and unspecialized for up to 10 Ma after the KPg extinction (23), potentially providing a window of opportunity for the evolution of large flightless herbivores in continental bird lineages. The early Tertiary fossil record supports this interpretation, with geographically widespread flighted palaeognath fossils (Fig. 3A) (22) and the appearance of other flightless avian herbivores such as gastornithids in Europe and North America, dromornithids in Australia, and Brontornis from South America (21). After the early Tertiary, the increasing prevalence of morphologically diverse mammalian competitors is likely to have prevented flightlessness from developing in other continental bird lineages.

The kiwi and tinamous are the only recent palaeognath lineages to not exhibit gigantism, and both taxa co-occur with a second palaeognath lineage (moa and rhea, respectively) that is both much larger and not their closest relative. We suggest that the disparity in size between co-occurring lineages may be a result of the relative timing of arrival of ancestral flighted palaeognaths coupled with competitive exclusion: The first palaeognath to arrive on each landmass monopolized the available niche space for large flightless herbivores and omnivores, forcing subsequent arrivals to adopt an alternative role and remain much smaller. For example, the South American ancestors of the rhea lineage (Diogenornis) were already large and flightless at 55 Ma (21) when the tinamou lineage originated. The absence of sympatric lineages of small palaeognaths on other landmasses in the recent past may reflect unavailability of alternative niches upon arrival (e.g., due to diversification of herbivorous mammals during the early Tertiary) or subsequent competition with mammals and/or neognathous birds. It is presumably the latter that has necessitated the maintenance of flight in the tinamous.

References and Notes

  1. M. Bunce, T. H. Worthy, M. J. Phillips, R. N. Holdaway, E. Willerslev, J. Haile, B. Shapiro, R. P. Scofield, A. Drummond, P. J. Kamp, A. Cooper, The evolutionary history of the extinct ratite moa and New Zealand Neogene paleogeography. Proc. Natl. Acad. Sci. U.S.A. 106, 20646–20651 (2009).

  2. T. H. Worthy, R. N. Holdaway, The Lost World of the Moa (Indiana Univ. Press, Bloomington, IN, 2002).

  3. J. P. Hume, M. Walters, Extinct Birds (Bloomsbury, London, 2012).

  4. J. Cracraft, Continental drift, paleoclimatology, and the evolution and biogeography of birds. J. Zool. 169, 455–543 (1973).

  5. K. Lee, J. Feinstein, J. Cracraft, in Avian Molecular Evolution and Systematics, D. Mindell, Ed. (Academic Press, New York, 1997), pp. 173–208.

  6. T. J. Parker, On the cranial osteology, classification and phylogeny of Dinornithidae. Trans. Zool. Soc. London 13, 373–431 (1895). 

  7. J. Cracraft, Phylogeny and evolution of the ratite birds. Ibis 116, 494–521 (1974). 

  8. W. P. Schellart, G. S. Lister, V. G. Toy, A Late Cretaceous and Cenozoic reconstruction of the Southwest Pacific region: Tectonics controlled by subduction and slab rollback processes. Earth Sci. Rev. 76, 191–233 (2006). 

  9. A. Cooper, C. Mourer-Chauviré, G. K. Chambers, A. von Haeseler, A. C. Wilson, S. Pääbo, Independent origins of New Zealand moas and kiwis. Proc. Natl. Acad. Sci. U.S.A. 89, 8741–8744 (1992).

  10. A. H. Bledsoe, A phylogenetic analysis of postcranial skeletal characters of the ratite birds. Ann. Carnegie Mus. 57, 73 (1988). 

  11. O. Haddrath, A. J. Baker, Multiple nuclear genes and retroposons support vicariance and dispersal of the palaeognaths, and an Early Cretaceous origin of modern birds. Proc. Biol. Sci. 279, 4617–4625 (2012). 

  12. M. J. Phillips, G. C. Gibb, E. A. Crimp, D. Penny, Tinamous and moa flock together: Mitochondrial genome sequence analysis reveals independent losses of flight among ratites. Syst. Biol. 59, 90–107 (2010). 

  13. J. Harshman, E. L. Braun, M. J. Braun, C. J. Huddleston, R. C. Bowie, J. L. Chojnowski, S. J. Hackett, K. L. Han, R. T. Kimball, B. D. Marks, K. J. Miglia, W. S. Moore, S. Reddy, F. H. Sheldon, D. W. Steadman, S. J. Steppan, C. C. Witt, T. Yuri, Phylogenomic evidence for multiple losses of flight in ratite birds. Proc. Natl. Acad. Sci. U.S.A. 105, 13462–13467 (2008). 

  14. J. V. Smith, E. L. Braun, R. T. Kimball, Ratite nonmonophyly: Independent evidence from 40 novel Loci. Syst. Biol. 62, 35–49 (2013). 

  15. J. R. Ali, D. W. Krause, Late Cretaceous bioconnections between Indo- Madagascar and Antarctica: Refutation of the Gunnerus Ridge causeway hypothesis. J. Biogeogr. 38, 1855–1872 (2011). 

  16. P. Johnston, New morphological evidence supports congruent phylogenies and Gondwana vicariance for palaeognathous birds. Zool. J. Linn. Soc. 163, 959–982 (2011). 

  17. T. H. Worthy, R. P. Scofield, Twenty-first century advances in knowledge of the biology of moa (Aves: Dinornithiformes): A new morphological analysis and moa diagnoses revised. N.Z. J. Zool. 39, 87–153 (2012). 

  18. A. Cooper, C. Lalueza-Fox, S. Anderson, A. Rambaut, J. Austin, R. Ward, Complete mitochondrial genome sequences of two extinct moas clarify ratite evolution. Nature 409, 704–707 (2001). 

  19. For details, see materials and methods on Science Online.

  20. T. H. Worthy, J. P. Worthy, A. J. D. Tennyson, S. W. Salisbury, S. J. Hand, R. P. Scofield, in Paleornithological Research 2013 - Proceedings of the 8th International Meeting of the Society of Avian Paleontology and Evolution, U. B. Göhlich, A. Kroh, Eds. (Natural History Museum Vienna, Vienna, Austria, 2013), pp. 63–80.

  21. G. Mayr, in Paleogene Fossil Birds (Springer, Berlin, 2009), pp. 25–34.

  22. See supplementary text on Science Online.

  23. K. Black, M. Archer, S. Hand, H. Godthelp, in Earth and Life, J. Talent, Ed. (Springer Netherlands, 2012), pp. 983–1078.

  24. M. Meyer, M. Kircher, Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harbor Protoc. 2010, (2010). 

  25. M. Martin, Cutadapt removes adapter sequences from high-throughput sequencing reads. Bioinformatics Action 17, 10 (2012). 

  26. C. Hahn, L. Bachmann, B. Chevreux, Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads—A baiting and iterative mapping approach. Nucleic Acids Res. 41, e129 (2013). 

  27. R. E. Green, A. S. Malaspinas, J. Krause, A. W. Briggs, P. L. Johnson, C. Uhler, M. Meyer, J. M. Good, T. Maricic, U. Stenzel, K. Prüfer, M. Siebauer, H. A. Burbano, M. Ronan, J. M. Rothberg, M. Egholm, P. Rudan, D. Brajkovi?, Z. Ku?an, I. Gusi?, M. Wikström, L. Laakkonen, J. Kelso, M. Slatkin, S. Pääbo, A complete Neandertal mitochondrial genome sequence determined by high-throughput sequencing. Cell 134, 416–426 (2008). 

  28. H. Li, B. Handsaker, A. Wysoker, T. Fennell, J. Ruan, N. Homer, G. Marth, G. Abecasis, R. Durbin, The Sequence Alignment/Map (SAM) format and SAMtools. Bioinformatics 25, 2078–2079 (2009). 

  29. A. Ginolhac, M. Rasmussen, M. T. P. Gilbert, E. Willerslev, L. Orlando, mapDamage: Testing for damage patterns in ancient DNA sequences. Bioinformatics 27, 2153–2155 (2011). 

  30. MEGAN 5 (Metagenome Analyzer), written by D. H. Huson, original design by D. H. Huson and S. C. Schuster, with contributions from S. Mitra, D. C. Richter, P. Rupek, H.-J. Ruscheweyh, R. Tappu, and N. Weber; http://ab.inf.uni-tuebingen.de/software/megan5/.

  31. H. Li, R. Durbin, Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009). 

  32. M. Kircher, in Ancient DNA: Methods and Protocols, B. Shapiro, M. Hofreiter, Eds. (Springer, New York, 2012), vol. 840, pp. 197–228.

  33. M. Schubert, A. Ginolhac, S. Lindgreen, J. F. Thompson, K. A. Al-Rasheid, E. Willerslev, A. Krogh, L. Orlando, Improving ancient DNA read mapping against modern reference genomes. BMC Genomics 13, 178 (2012). 

  34. R. Lanfear, B. Calcott, S. Y. W. Ho, S. Guindon, PartitionFinder: Combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol. Biol. Evol. 29, 1695–1701 (2012). 

  35. A. Stamatakis, P. Hoover, J. Rougemont, A rapid bootstrap algorithm for the RAxML Web servers. Syst. Biol. 57, 758–771 (2008). 

  36. A. Stamatakis, RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22, 2688–2690 (2006). 

  37. F. Ronquist, M. Teslenko, P. van der Mark, D. L. Ayres, A. Darling, S. Höhna, B. Larget, L. Liu, M. A. Suchard, J. P. Huelsenbeck, MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539–542 (2012).

  38. A. J. Drummond, A. Rambaut, BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol.7, 214 (2007). 

  39. A. Rambaut, A. J. Drummond, Tracer v1.4, available from http://tree.bio.ed.ac.uk/software/tracer/ (2007).

  40. X. Xia, Z. Xie, M. Salemi, L. Chen, Y. Wang, An index of substitution saturation and its application. Mol. Phylogenet. Evol. 26, 1–7 (2003). 

  41. X. Xia, P. Lemey, in The Phylogenetic Handbook: A Practical Approach to Phylogenetic Analysis and Hypothesis Testing, P. Lemey, M. Salemi, A. Vandamme, Eds. (Cambridge Univ. Press, Cambridge, 2009).

  42. X. Xia, Z. Xie, DAMBE: Software package for data analysis in molecular biology and evolution. J. Hered. 92, 371–373 (2001). 

  43. D. L. Swofford, PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods) (Sinauer Associates, Sunderland, MA, 2002).

  44. M. D. Sorenson, E. A. Franzosa, TreeRot, version 3 (Boston Univ., Boston, 2007).

  45. J. Gatesy, G. Amato, M. Norell, R. DeSalle, C. Hayashi, Combined support for wholesale taxic atavism in gavialine crocodylians. Syst. Biol. 52, 403–422 (2003). 

  46. P. O. Lewis, A likelihood approach to estimating phylogeny from discrete morphological character data. Syst. Biol. 50, 913–925 (2001). 

  47. P. Houde, S. L. Olson, Paleognathous carinate birds from the early tertiary of North America. Science 214, 1236–1237 (1981). 

  48. P. W. Houde, Palaeognathous birds from the early Tertiary of the Northern Hemisphere. Publ. Nuttall Ornithol. Club 22, 1 (1988). 

  49. G. Grellet-Tinner, G. J. Dyke, The eggshell of the Eocene bird Lithornis. Acta Palaeontol. Pol. 50, 831 (2005). 

  50. P. G. Foster, D. A. Hickey, Compositional bias may affect both DNA-based and protein-based phylogenetic reconstructions. J. Mol. Evol. 48, 284–290 (1999). 

  51. M. J. Phillips, D. Penny, The root of the mammalian tree inferred from whole mitochondrial genomes. Mol. Phylogenet. Evol. 28, 171–185 (2003). 

  52. P. J. Lockhart, M. A. Steel, M. D. Hendy, D. Penny, Recovering evolutionary trees under a more realistic model of sequence evolution. Mol. Biol. Evol. 11, 605–612 (1994). 

  53. S. J. Gould, Of Kiwi eggs and the Liberty Bell. Nat. Hist. 95, 20 (1986). 

  54. S. J. Gould, in Bully for Brontosaurus, S. J. Gould, Ed. (W. Norton, New York, 1991), pp. 109–123.

  55. H. M. F. Alvarenga, Uma ave ratitae do Paleoceno Brasileiro: Bacia calcária de Itaboraí, Estado do Rio de Janeiro, Brasil. Bol. Mus. Nac. Geol. 41, 1 (1983). 

  56. C. W. Andrew, On the pelvis and hind-limb of Mullerornis betsilei M.-Edw. & Grand.; with a note on the occurrence of a ratite bird in the Upper Eocene beds of the Fayum, Egypt. Proc. Zool. Soc. Lond. 74, 163–171 (1904). 

  57. L. D. Martin, in Papers in Avian Paleontology Honoring Pierce Brodkorb, in Natural History Museum of Los Angeles County Science Series 36, K. E. Campbell, Ed. (Natural History Museum of Los Angeles County, Los Angeles, 1992), pp. 97–125.

  58. P. W. Houde, H. Haubold, Palaeotis weigelti restudied: A small middle Eocene ostrich (Aves: Struthioniformes). Palaeovertebrata 17, 27 (1987). 

  59. C. P. Tambussi, J. I. Noriega, A. Ga?dzicki, A. Tatur, M. A. Reguero, S. F. Vizcaino, Ratite bird from the Paleogene La Meseta Formation, Seymour Island, Antarctica. Polish Polar Res. 15, 15–20 (1994).

  60. M. J. Benton, in The Fossil Record 2, M. J. Benton, Ed. (Chapman & Hall, London, 1993), pp. 681–715.

  61. M. J. Benton, P. C. J. Donoghue, Paleontological evidence to date the tree of life. Mol. Biol. Evol. 24, 26–53 (2007). 

  62. J. Gauthier, Saurischian monophyly and the origin of birds. Mem. Calif. Acad. Sci. 8, 1 (1986). 

  63. P. C. Sereno, Basal archosaurs: Phylogenetic relationships and functional implications. Soc. Vertebr. Paleontol. Mem. 2, 1–53 (1991). 

  64. M. J. Benton, Scleromochlus taylori and the origin of dinosaurs and pterosaurs. Philos. Trans. R. Soc. London Ser. B 354, 1423–1446 (1999). 

  65. C. A. Brochu, Patterns of calibration age sensitivity with quartet dating methods. J. Paleontol. 78, 7–30 (2004). 

  66. T. E. Williamson, Brachychampsa sealeyi, sp nov., (Crocodylia, Alligatoroidea) from the Upper Cretaceous (lower Campanian) Menefee Formation, northwestern New Mexico. J. Vertebr. Paleontol. 16, 421–431 (1996). 

  67. A. D. Buscaloni, F. Ortega, D. Vasse, New crocodiles (Eusuchia: Alligatoroidea) from the Upper Cretaceous of southern Europe. Comptes Rendus l’Acad. Sci. Paris. Sci. Terre Plandtes 325, 525 (1997). 

  68. X.-C. Wu, A. P. Russell, D. B. Brinkman, A review of Leidyosuchus canadensis Lambe, 1907 (Archosauria: Crocodylia) and an assessment of cranial variation based upon new material. Can. J. Earth Sci. 38, 1665–1687 (2001). 

  69. D. R. Schwimmer, King of the Crocodylians: The Paleobiology of Deinosuchus. (Indiana Univ. Press, Bloomington, IN, 2002).

  70. J. A. Clarke, C. P. Tambussi, J. I. Noriega, G. M. Erickson, R. A. Ketcham, Definitive fossil evidence for the extant avian radiation in the Cretaceous. Nature 433, 305–308 (2005). 

  71. Z. Zhonghe, Evolutionary radiation of the Jehol Biota: Chronological and ecological perspectives. Geol. J. 41, 377–393 (2006).

  72. K. E. Slack, C. M. Jones, T. Ando, G. L. Harrison, R. E. Fordyce, U. Arnason, D. Penny, Early penguin fossils, plus mitochondrial genomes, calibrate avian evolution. Mol. Biol. Evol. 23, 1144–1155 (2006). 

  73. C. Patterson, P. V. Rich, The fossil history of the emus, Dromaius (Aves: Dromaiinae). Rec. South Austral. Mus. (Adelaide) 21, 85 (1987). 

  74. T. H. Worthy, S. J. Hand, M. Archer, Phylogenetic relationships of the Australian Oligo-Miocene ratite Emuarius gidju Casuariidae. Integr. Zool. 9, 148–166 (2014).

  75. W. E. Boles, “Revision of Dromaius gidju Patterson and Rich 1987 from Riversleigh Northwestern Queensland Australia with a Reassessment of its generic position.” in Natural History Museum of Los Angeles City Science Series 195 (Natural History Museum of Los Angeles City, Los Angeles, 1992).

  76. O. Haddrath, A. J. Baker, Complete mitochondrial DNA genome sequences of extinct birds: Ratite phylogenetics and the vicariance biogeography hypothesis. Proc. R. Soc. London Ser. B 268, 939–945 (2001). 

Acknowledgments: This study was funded by the New Zealand Marsden Fund and the Australian Research Council. Grid computing facilities were provided by eRSA (e-research South Australia) and CIPRES (Cyberinfrastructure for Phylogenetic Research). We thank the Museum of New Zealand Te Papa Tongarewa (P. Millener, S. Bartle, A. Tennyson), Natural History Museum Oslo (N. Heintz), and National Museum of Natural History Paris (C. Lefèvre) for samples. We acknowledge D. Penny, M. Phillips, D. Burney and R. Ward for valuable advice and assistance. All data associated with this study are available on GenBank (accession nos. KJ749824 and KJ749825) and DRYAD (doi:10.5061/dryad.2727k).