The evolutionary history of the extinct ratite moa and New Zealand Neogene paleogeography M. Bunce3-1-2, T. H. Worthybc\ M. J. Phillipsd, R. N. Holdawaye, E. Willerslev', J. Haile'-s, B. Shapiro^, R. P. Scofield1, A. Drummond', P. J. J. Kampk, and A. Cooperb-2 aAncient DNA Laboratory, School of Biological Sciences and Biotechnology, Murdoch University, Perth 6150, Australia; bAustralian Centre for Ancient DNA, University of Adelaide, South Australia 5005, Australia; cSchool of Biological, Earth and Environmental Sciences, University of New South Wales, Sydney, New South Wales 2052, Australia; dCentre for Macroevolution and Macroecology, Research School of Biology, Australian National University, Canberra 2601, Australia; ePalaecol Research Ltd and School of Biological Sciences, University of Canterbury, Christchurch 8041, New Zealand; fCentre for Ancient Genetics, Department of Biology, University of Copenhagen, DK-2100 Copenhagen, Denmark; SDepartment of Zoology, University of Oxford, Oxford 0X13PS, United Kingdom; department of Biology, Pennsylvania State University, University Park, PA 16802; 'Canterbury Museum, Rolleston Avenue, Christchurch 8013, New Zealand; 'Bioinformatics Institute and Department of Computer Sciences, University of Auckland, Auckland 1020, New Zealand; and kDepartment of Earth and Ocean Sciences, University of Waikato, Hamilton 3240, New Zealand Edited by James P. Kennett, University of California, Santa Barbara, CA, and approved September 24, 2009 (received for review June 28, 2009) The ratite moa (Aves: Dinornithiformes) were a speciose group of massive graviportal avian herbivores that dominated the New Zealand (NZ) ecosystem until their extinction =*600 years ago. The phytogeny and evolutionary history of this morphologically diverse order has remained controversial since their initial description in 1839. We synthesize mitochondrial phylogenetic information from 263 subfossil moa specimens from across NZ with morphological, ecological, and new geological data to create the first comprehensive phytogeny, taxonomy, and evolutionary timeframe for all of the species of an extinct order. We also present an important new geological/paleogeographical model of late Cenozoic NZ, which suggests that terrestrial biota on the North and South Island landmasses were isolated for most of the past 20-30 Ma. The data reveal that the patterns of genetic diversity within and between different moa clades reflect a complex history following a major marine transgression in the Oligocene, affected by marine barriers, tectonic activity, and glacial cycles. Surprisingly, the remarkable morphological radiation of moa appears to have occurred much more recently than previous early Miocene (ca. 15 Ma) estimates, and was coincident with the accelerated uplift of the Southern Alps just ca. 5-8.5 Ma. Together with recent fossil evidence, these data suggest that the recent evolutionary history of nearly all of the iconic NZ terrestrial biota occurred principally on just the South Island. ancient DNA | Oligocene Drowning | Dinornithiformes | phylogeny taxonomy The prolonged geographic isolation of New Zealand (80-60 million years) and the paucity of terrestrial mammals created a unique ecosystem dominated by an estimated 245 species of birds (1), providing an unparalleled opportunity to observe evolutionary processes. The most striking of the recent avian radiations is that of the extinct ratite moa (Aves: Dinornithiformes), a speciose order ranging in size from a large turkey to the 3 m tall Dinornis weighing up to 300 kg. Ratites are a basal lineage of birds that are hypothesized to have had a common ancestor ca. 80 million years ago on the Cretaceous southern supercontinent of Gondwana, which subsequently underwent either vicarious speciation as the landmass fragmented (2), and/or flighted dispersal (3). The extant members of the ratite lineage include the ostrich (Africa), emu, cassowary (Australia, New Guinea), rhea (South America), and kiwi (NZ). Extinct ratites include the giant elephant birds (Madagascar) and moa (NZ). Since the first description in 1839 (4) the taxonomy of moa has remained contentious with up to 64 different species and 20 genera assigned at various times (1). The complex geological history of NZ, and significant regional variations in climate, diet, and sexual dimorphism (SD) have resulted in moa being highly variable morphologically, complicating attempts to define species limits. Morphometric and osteological studies have generally recognized two families, Emeidae and Dinornithidae (Table SI), which most recently have been considered to contain eight and three species respectively (1, 5-7). The recovery, amplification, and sequencing of ancient DNA from fossil bone has provided new insights into moa systematics, and revealed extreme cases of SD (8-14). Fig. 1 depicts the current taxonomy of moa (including changes proposed in this article), as well as summarizing species distributions, dimensions, and ecology. In this study we use mitochondrial DNA sequences isolated from 119 specimens, in addition to previously published data (9-13), to study the mode and tempo of moa evolution, phylogeography, and taxonomy, and relate this to the geological and ecological history of NZ. Results and Discussion The accuracy of studies employing molecular clocks to date speciation events is becoming increasingly scrutinized. Problems associated with calibration points, substitution saturation, base composition bias, model selection, and more recently with the time dependency of molecular clocks, can all distort the temporal accuracy of phylogenetic reconstructions (15-18). We used two contrasting approaches in an attempt to establish timeframes for the divergence events within and between moa species: an externally (fossil) calibrated analysis of avian mitochondrial protein-coding sequences, and an internally calibrated Bayesian analysis of radiocarbon-dated moa mitochondrial control region sequences. For the externally calibrated approach, two datasets (designated mtg-10,692 and mtg-2,153) were generated from existing moa GenBank data (9-13, 19) (Dataset SI and Dataset S2). We followed the procedure of Baker et al. (9) and dated the branching events at the base of the moa phylogeny using a 10,692 bp dataset (mtg-10,692) comprising mitochondrial coding region sequences from 26 avian taxa (including 9 ratites), which was externally calibrated from the fossil record (SI Text and Table S2). We used a variety of methods to deal with biases introduced by substitution saturation and phylogenetic artifacts (SI Text), and estimated divergence dates using BEAST v.1.4.8 (20) as described in Methods. Author contributions: A.C. designed research; M.B., E.W., J.H., P.J.J.K., and A.C. performed research; M.B., T.H.W., M.J.P., R.N.H., B.S., R.P.S., A.D., P.J.J.K., and A.C. analyzed data; and M.B., T.H.W., M.J.P., and A.C. wrote the paper. The authors declare no conflict of interest. This article is a PNAS Direct Submission. Data deposition: The sequences reported in this paper have been deposited in the GenBank database (accession nos. GU138993-GU139113). 1M.B. and T.H.W. contributed equally to this work. 2To whom correspondence may be addressed. E-mail: alan.cooper@adelaide.edu.au or m.bunce@murdoch.edu.au. This article contains supporting information online at www.pnas.org/cgi/content/full/ 0906660106/DCSupplemental. www.pnas.org/cgi/doi/10.1073/pnas.0906660106 Family: Dinornithidae _ Dinornis S\ 'laiiiilics: l\\ o species I'l- rohmlus (Soiilli Island, blue) and D. novticzcciUnidinc (Nonh Island. red| Dimensions; 56-249 kg and 90 to 200 cm in height -significant sexual dimorphism wiili females up to three times the mass of males. I labium Hum sing jjtncralisl - lias been found in upland, lev. :.:rii ;ind open loic-l iuiNlaK The lui-ficr forms IVCIipit'li iov. f;li| Family: Megalapterygidae Megalapteryx Systematies: Monotypic. M didinus, (South Island). Dimensions: 28-80 kg and 65 to 95 cm. Pleistocene specimens arc significantly larger than Holocene forms. Habitat: Subalpinc scrub, grassland and high country forests (usually >900m). Systcmatics: Monotypic, A. didiformix. Dimensions: 26-64 kg and 50 to 90 cm Habitat: Non-coastal lowland forests with a continuous canopy. Cook Strait South Island Systematica: Uonolypic, E. i gravis and .goranoitfes). Dimensions: 12-109 kg and 51 to 103 cm. Habitat: Dtiet climates - typically lowland open forest and coastal sites. Pachyornis Emeus Systematies: Mouolypic. /:'. cras.suv. (South Island). Dimensions: ol>-70 kg niul 73 to 00 em in height Habitat: Preference for low land forest (rrsrrally <200m) and swamps. Syslemalics: /' geracod/o (North Island). /( e/e/n'irnno/tio (blue) and P. r»«r,i,'n (green) (South Island). Dimensions: 17-163 kg and 54 (o 1 21 em in height. Habitat: P. rcotm/r, occupied srihalpine grassland, P gerauo/i/es nrrrl /' e/o/Jrooro/uo preferred im land forest edges and sselland vegetation. Fig. 1. Systematies, dimensions, and approximate distributions of moa in the three family, six genera, nine species taxonomy advocated in this study. Estimates of body size and habitat were generated from information contained in Worthy and Holdaway (2002) and are discussed in more detail in the SI Text. A comparison with the generally accepted two family, 11 species classification used from the 1980's to 2002 is presented in Table S1. The resulting posterior distribution of the molecular date estimate for the basal moa branch was then used to calibrate analyses within a shorter dataset of 2153 bp (9) of 29 moa mitochondrial protein coding sequences (mtg-2,153, Table S2) to obtain date estimates for species-level divergences within the moa phytogeny. For the internally calibrated approach, an alignment of 80 Megalapteryx, Dinornis, and Pachyornis control region sequences (389 bp) from radiocarbon-dated moa bones (Table S3) were used as calibration points (or "time-stamped" data) to temporally calibrate a Bayesian framework using BEAST (21). Moa Systematies at Family and Genus Level. Bayesian analysis of the mtg-2,153 mitochondrial coding region dataset (2,153 bp for 29 moa) yielded the maximum a posteriori (MAP) tree shown in Fig. 2. Similar analysis of the control region data from 263 moa control region sequences (389 bp) yielded the higher resolution tree shown in Fig. 3 (presented in detail in Fig. SI), and more detailed analyses with restricted taxonomic coverage shown in Fig. 4. The phylogenetic arrangement in Figs. 2 and 3 has a number of discrepancies with current taxonomic arrangements (Table SI) in which just two families (Dinornithidae and Emei-dae) are recognized, with two subfamilies within Emeidae (Megalapteryginae and Anomalopteryginae). Although these three main lineages are recovered in analyses of the mtg-2153 mitochondrial coding region dataset (Fig. 3), Megalapteryx is clearly the basal moa lineage. This arrangement is consistent with previous genetic studies (9, 11, 12), and recent analyses of morphological characters (1, 22), where Megalapteryx didinus falls as the basal moa lineage, and is outside Emeidae. Conse- ^Anomalopteryx Euryapteryx Marine environments _ Marine environments Marine environments Fig. 2. A spatial and temporal context for the evolution of moa. Molecular phylogeny and date estimates of the moa radiation generated from the mtg-2,153 dataset (see Methods), compared with the new paleogeographic model of Neogene New Zealand. A Bayesian Inference tree is shown with Bayesian posterior probability values (>80%) indicated on the nodes with support, whereas the node bars correspond to the 95% HPD. A series of four paleogeographic maps, based on extensive geological mapping of the area (see Methods and SI Text), show different time horizons during the Neogene and the presence/absence of North and South Island landmasses in central New Zealand. The uplift of the Southern Alps (ca. 5-8.5 million years) and periodic bridging to the North Island in the Pleistocene (<2 million years) appear to be instrumental in moa speciation. The absence of deep (ca. 20 Million years) splits in the moa phylogeny suggest that all recent moa species originated from the southern landmass (see main text), consistent with the phylogeographic distributions (Fig. 3, Fig. 46, and SI Text). quently, we propose that moa taxonomy be revised so that three families are recognized in the Dinornithiformes (Table SI): Dinornithidae (Dinornis); Emeidae (Pachyornis,Emeus, Anoma-lopteryx, Euryapteryx), and Megalapterygidae (new family; type and only genus Megalapteryx). Neither of the presently recognized subfamilies (Anomalopteryginae, Emeinae) is supported by either morphological cladistics (1, 22) or mtDNA analyses (9, 13), and should be discarded. At the genus level, there is strong support (posterior probabilities of 100%) (Fig. 3) for all six moa genera outlined by Cracraft (7) and refined by Worthy and Holdaway (1). The suggestion that species can be recognized solely from genetic distances has previously been identified as problematic (23, 24). Suggestions that this approach can also be applied to the systematic studies of closely related extinct taxa (9) are particularly concerning given the lack of behavioral and ecological data. This problem is exemplified by the strong phylogeographic pattern exhibited within the alpine moaM didinus (Fig. 4 and Fig. SI) where control region genetic distances of 4.7% separate mtDNA "clades". This distance is greater than that between well recognized genera such as Emeus and Euryapteryx (4.5%), and reinforces the idea that such genetic diversity measures are a product of population longevity and niche persistence, as well as other population parameters. Keeping this observation in mind, and integrating recent morphological studies and moa palaeoecological and phylogeographical data, we review the taxonomic status of all moa genera (SI Text) and propose the Bunce et al. Rootofmoa ŕ Family: Dinornithidae lree\| □ North island fossils □ South island fossils Fig. 3. Phylogenetic reconstruction of all moa species. The maximum a posteriori (MAP) tree of 263 mitochondrial control region sequences of each of the moa genera as generated in BEAST (Methods) using a HKY+G + I model on a 389 bp alignment. A more detailed tree in which the tips are visible is presented in Fig. S1. DNA sequences obtained from North and South Island fossils are colored red and blue respectively, and suggest that North Island specimens tend to occupy derived phylogenetic positions (see also Fig. 46 and SI Text). For clarity, posterior probabilities are shown only on basal nodes. Saturation effects in the control region data have resulted in the root of the moa tree being misplaced in this MAP tree, and the mitochondrial protein coding data (Fig. 2) is predicted to give a better estimate of the branching topology at the base of the moa radiation. new taxonomic arrangement shown in Fig. 1 and Table SI. Overall, we have adopted a conservative approach to species designations and regard mtDNA diversity alone as insufficient to define species limits, especially in the context of recognized barriers to gene flow (discussed later in Island Allopatry and Pleistocene Cycles). As new DNA profiling approaches are devised (25) there may be minor future alterations to the taxonomic framework (3 families, 6 genera) proposed here. Establishing a Timeframe for Moa Evolution. Reanalysis of existing mtDNA data. The mean divergence date between Dinornithidae and Emeidae was estimated at just 5.27 Ma [95% highest probability density (HPD): 3.1-9.0] in the BEAST analysis of the mtg-10,692 dataset, with 3rd codon positions RY-coded, and using two ratite and three nonratite calibration priors (SI Text and Table 1). This date is less than half the age of the 15 Ma (95% CI 14.5-15.6 Ma) date inferred using penalized likelihood in r8s used by Baker et al. (9), which we also recovered in r8s analyses of the mtg-10,692 dataset using standard nucleotide coding (15.8 Ma). Importantly although, this estimate drops to 8.1 Ma when RY-coding is used to reduce the impact of saturated 3rd codon positions (SI Text). The estimate drops further, to 7.6 Ma, if the problematic calibration date of 82 Ma for the beginning of rafting between New Zealand and Australia/Antarctica is discarded in favor of other, more secure, avian calibrations (SI Text). Because this geological event would not prevent the dispersal of flying paleognaths, it is now questionable (3). It appears that the differences between our analysis and that of Baker et al. (9) are not primarily with the programs used, but rather the use of the saturated 3rd codon positions and the problematic 82 Ma calibration. When these issues are dealt with, r8s infers a date of 7.6 Ma, well within the 95% HPD of the BEAST analysis. Overall, we prefer the BEAST estimate to those of r8s because the latter assumes autocorrelated rate changes among branches, which 1QQ i—Medi »37 ScArnard JTC Medil?l_alenmufc_í302_Jl O Hadi M5.Hona/comb_l«SS>JJ Midi MS.presideJitsPeiroil.rn (Arthur Midi W7_ttinMiTLiya5U_19D Medi1H0_MdMyttan_Melsom Medi »34_Hane,eoaib_lí60i_68L MediMS.GaribaldLrlelson Medi 553_Char sion_*esrcoa;° Medi pSt.EirdleAirJdrAnriur Hedi«lĚ_Mtí>*Fn_3 250_70 Medi7dl_H0«Blatt^w_S_iaiL5a Medi»5SJ.cwcrBallerHivtr Medi MS_WrvwrerC3aw_WeitCoi.it Madi _6+_Teanau_1632__Cr Medi 353 M-Lutmaie MedlR!7_TaaaiiřVjllev_717«_*Q Medi ..]_M!ta_mar. Madi .55Jumr_Ca«a Madi357_La<*Orball Madi _rCLT_k_lreVall_._77__2_ Madi_AY2___W_rjMa.irjrwv_5erperi HedLCI_th_J:_p_.lľt__E24_5_ Madi 735_Crarnw_ll Medi _ 59_HC-lcb_rnTr__k_H-25 Ma. Furthermore, the endemic North Island biota would presumably descend from taxa that had survived the marine transgression on the limited emergent Northland landmass. Recently detected phylogeographic patterns in NZ skinks are consistent with this model and indicate an early Miocene vertebrate fauna on the North Island (34). Within the moa phylogeny there are no genetic splits that might correspond to a >25 Ma separation (Table 2), indicating that the Quaternary moa lineages appear to have been descended from populations occupying only one island. In light of the molecular dates, the South Island Manuherikia Group fossils (33) suggest that the progenitors of the three moa families probably derive from a single South Island ancestral moa lineage which had survived the Oligocene Drowning event. A South Island ancestral moa lineage is also supported by the phylogeographic distributions of species and specimens in Fig. 3 and Fig. AB (and SI Text) where the majority of North Island taxa and clades occupy derived phylogenetic positions. Discoveries of tuatara (35), acanthisittid wrens (36), and kauri (37) from similar South Island fossil deposits, and a parallel lack of deep genetic splits in these taxa imply that the recent evolutionary history of many iconic NZ endemics was also on the South Island. This latter observation raises important questions about the nature of terrestrial biotic diversity on the North Island from >25 Ma to ca. 1.5-2 Ma, at which point a major biotic interchange is likely to have taken place. The Origin of Diversity within the Dinornithiformes. The rapid morphological radiation of the moa is matched by other island endemics Bunce et al. PNAS I December 8,2009 I vol.106 I no. 49 I 20649 such as the moa-nalos of Hawaii (38), and stands in contrast to the fossil record of many volant species in NZ and Australia (33). The timing of the moa radiation appears to correlate closely with the formation of the Southern Alps, and concomitant increase in habitat diversity across NZ, ca. 5 Ma (39). Although the transpres-sional Alpine Fault became active in the early Miocene (ca. 20 Ma), the uplift resulting in the Southern Alps only began to accelerate ca. 5-8.5 Ma (see discussion in SI Text). The Alps act as a barrier to the circumpolar (westerly) airflow creating wet rainforests to the west, and dryer, warmer climate to the east. Therefore, accelerated uplift generated both upland versus lowland environments, and wet versus dry habitats, forming a mix of forest, shrubland, and grasslands. In this regard, it is important to note that NZ did not experience the pronounced Miocene shift toward aridification detected elsewhere in the world (40, 41). If the niche specialization of the recent species reflects ancestral states, it is possible that the first divergence in the moa tree (Megalapteryx ca. 5.8 Ma) was because of the evolution of upland (Megalapteryx) and lowland (Dinornis fEmeid) ecotypes, whereas the subsequent increase in habitat diversity was the trigger for speciation within Emeidae. Island Allopatry and Pleistocene Cycles. Following the initiation of moa family-level diversity in the early Pliocene (ca. 3.3-5.8 Ma), phylogeographic analyses (Figs. 3 and 4 and Fig. SI) suggest that a further cycle of intrageneric diversity was created by the movement of southern endemics into the North Island ca. 1.5-2 Ma, as the Manawatu Strait started to close and land links first became possible (Fig. 2). The subsequent formation of the seaway now known as Cook Strait is a geologically young event that first formed in a mid-Pleistocene interglacial period (ca. 450 kyr) (42) and prevented gene flow between each island apart from during periods of lowered sea levels during major glaciations (SI Text). An additional consequence of the mountainous topography and rain-shadow zones in both islands is that NZ was able to sustain diverse habitats and refugia throughout the Pleistocene glacial and interglacial periods, despite its comparatively small size (43). The Pleistocene was dominated by long (ca. 100 kyr), cold climatic periods during which grasslands and shrublands were the main vegetation, interspersed with shorter (10 kyr), warm, forest-dominated phases (44). The divergence date estimates in Table 1 suggest that the Pleistocene climate shifts were instrumental in shaping the distribution of moa species and moa clades within NZ. To explore this issue further, we examine three genera in detail to illustrate the role of Pleistocene niche stability in the maintenance of genetic diversity. Emeus. A striking feature of Fig. AB is the lack of significant mtDNA phylogeographic diversity or structure in Emeus crassus, a species that is abundant in Holocene fossil deposits. This observation potentially relates to the drastic reduction of Emeus' favored habitat of wet forests during glaciations, potentially to a single geographic region somewhere in the south-eastern South Island during the last glaciation (Otiran) (ca. 75 to 14 kyr). During the Holocene, a rapid increase in this habitat is likely to have led to a rapid population expansion. This model of a severe bottleneck, followed by a rapid expansion would account for the abundance of fossil remains, but an extremely restricted genetic diversity compared with the other moa genera. Euryapteryx. The phylogeographic patterns seen in Euryapteryx are more complex than observed in other moa species, and probably reflect several dispersal events between the North and South Islands. Euryapteryx curtus is the only moa taxon found on both islands that does not exhibit North/South monophyly, suggesting a relatively recent interchange between the islands. The largely coastal lowland habitat of this species makes it one of the more likely to have utilized Pleistocene land bridges. The date estimates and tree topologies in Table 1 and Figs. 3 and 4 suggest that the common ancestor of Emeus and Euryapteryx diverged in the South Island shortly after the initial island connections (ca. 1.5-2.0 Mya), 20650 | www.pnas.org/cgi/doi/10.1073/pnas.0906660106 with Euryapteryx subsequently diversifying into several major clades. During two different glacial low sea stands (very roughly 450 kyr and 100 kyr), Euryapteryx haplotypes appear to have dispersed to the North Island, with the phylogeographic distribution suggesting that in the penultimate glacial this movement occurred via a land bridge on the western side of Cook Strait. Megalapteryx. The evolutionary history of the South Island upland specialist species Megalapteryx didinus is not well understood, with only two specimens included in a previous genetic study of moa taxonomy (9). Fig. 3 and Fig. AA reveal large amounts of genetic diversity in Megalapteryx relative to other moa genera, which likely relates to the persistence of upland/montane habitat during both glacial and interglacial periods. Despite the deep genetic diversity, there is presently no known parallel morphological diversity, and we consider all four clades members of a single species (see SI Text). The geographic boundaries separating the four clades (A-D, Fig. AA) correspond to the presence of known barriers which effectively split the uplands of the South Island into four regions during the Pleistocene (SI Text). During glacial cycles, extensive ice fields covered the Southern Alps and glaciers flowed outwards from the mountains at numerous points. In central South Island the narrow West Coast was split by glaciers extending to near sea level, while on the eastern side of the Alps major glaciers extended to the lowlands and were accompanied by extensive glacial outwash plains. Interestingly, comparison of the Megalapteryx phytogeny with that of the endemic giant invertebrate alpine scree weta (43) demonstrates striking similarities in the phylogeographic structure of the clades, suggesting common barriers to geneflow between these two upland species. Clearly the complexity of the genetic patterns seen in the different moa taxa provide strong evidence that the use of universal genetic distance(s) to define meaningful taxonomic units is not viable, and moa should not be used as an example of the effectiveness of DNA barcoding (45) (SI Text). Although COI mtDNA signatures clearly discriminate some taxa (as do all mtDNA genes), no single genetic distance can define moa taxa because the genetic diversity of each genus reflects a differing palaeoecological history involving niche persistence, ecology, and population longevity. An accurate model of species diversity is only possible through the integration of a range of data from mitochondrial and nuclear DNA, stable isotopic and ecological signatures, osteology, and morphometries. This is the first detailed genetic study of all members of an entire extinct order, and demonstrates the influential role of tectonic uplift, mountain building, and glacial cycles on the evolution and genetic diversity of NZ terrestrial endemics. The data clarifies moa taxonomic relationships that have remained problematic for >160 years, and suggests that the glacial landscape allowed the rapid diversification of a single moa lineage that survived the Oligocene Drowning. This pattern appears likely to be widespread within other extant NZ endemics. The important new geological model of Neogene NZ emphasizes our current lack of knowledge about the pre-Pliocene landscape, and raises important questions about the role of marine barriers and the biotic diversity of the north and south islands. The combined geological and genetic data suggests that the NZ Neogene terrestrial record is likely to have been marked by the significant loss of terrestrial endemics from a highly unusual environment, which is only just beginning to be characterized. Materials and Methods Samples, DNA Extractions, PCR, and Sequencing. A detailed list of the samples used in this study (location, museum numbers, 14C dates) and extraction methods [after (11)]can be found inTableS3or in previous publications (9-13,19). PCR was used to amplify two overlapping sections of mitochondrial control region (totaling ™389 bp, depending on genus) from the moa samples, and sequencing was performed using an ABI 310 or 3730. Strict ancient DNA (aDNA) guidelines were followed to minimize contamination of samples with exogenous DNA, including Bunce et al. multiple negative extraction and amplification controls to detect contamination. All DNA extractions and PCRs were set-up and performed in a geographically isolated, specialist aDNA laboratory. Three bone subsamples were sent to Copenhagen for independent "bl ind" repl ication, and the resulting sequences matched those previously generated. Sequence Analysis. The 263 aligned moa control region sequences (389 bp, FASTA file in SI Text) were used to jointly estimate a phylogenetic tree and substitution model parameters using Metropolis-Hastings Markov chain Monte Carlo (MCMC) integration. ModelTest 3.06 (46) recommended an HKY+G+I model of substitution for the control region sequences under the Akaike Information Criterion. The tree presented in Fig. 3 (also see Fig. 4and Fig. S1) is the MAP tree (assuming a "strict" molecular clock) obtained from the MCMC analysis (two independent runs each comprising 50 million generations) using BEAST (20). No explicitly informative priors were used in the analysis, so the tree will correspond closely to a maximum likelihood (ML) estimate of the tree under the same model of substitution. Detailed phylogenies for three moa genera shown in Fig. 4 were generated as described above, using just the sequences for these taxa plus an outgroup, as detailed. Timescales for moa evolution were estimated with BEAST v.1.4.8 (20), using 40 million generations. To compare our date estimate for the Dinornithidae/ Emeidae divergence with those from Baker et al. (9), we also inferred dates from the mtg-10,692 dataset using r8s (47) under penalized likelihood (see SI Text). Forthe internally calibrated approach, 80 Megalapteryx, Dinornis, and Pachy-ornis control region sequences (389 bp) with associated radiocarbon ages (Table S3) were used as calibration points (or time-stamped data) and analyzed with BEAST using priors forthe mutation rate (2% to 50% per million years), and the maximum (oldest) age of each genus based on the mtg-2,153 analysis (SI Text). The results are presented in Table 1. ACKNOWLEDGMENTS. We thank the Museum of New Zealand Te Papa Ton-garewa (A.J.D. Tennyson and J.A. Bartle), Canterbury Museum, Otago Museum (S. Heath), the Natural History Museum (S. Chapman), Waitomo Museum (K. Wilde), Otago University (I. Smith), Auckland Museum (B. Gill), and the Oxford Radiocarbon Accelerator Unit (T. Higham) for samples and assistance. Financial support was provided by the Natural Environment Research Council (U.K.), the Wellcome and Leverhulme Trusts, the Arts and Humanities Research Council (U.K.), Australian Research Council, The Marsden Fund, and Foundation for Research, Science and Technology (New Zealand). 1. Worthy TH, Holdaway RN (2002) The Lost World of the Moa (Indiana Univ Press, Bloom-ington, IN), p718. 2. Cracraft J (1974) Continental drift and vertebrate distribution. Ann Rev Ecol Syst 5:189-205. 3. Harshman J, etal. (2008) Phylogenomic evidence for multiple losses of flight in ratite birds. Proc Natl Acad Sci USA 105:13462-13467. 4. Owen R (1840) On the bone of an unknown struthious bird from New Zealand, meeting of November 12, 1839. Proc. Zoo/. Soc. ton. for 1839 Part Vllfrio. 1 xxxiii): 169-171. 5. Archey G (1941) The Moa: A study of the Dinornithiformes. Bulletin of the Auckland Institute and Museum No. 1, pp 116-121. 6. Oliver WRB (1949) The Moas of New Zealand and Australia. Dominion Museum Bulletin No. 15 (Wellington, New Zealand). 7. Cracraft J (1976) Covariation patterns in the post-cranial skeleton of the moas (Aves, Dinornithidae): A factor analytic study. Paleobiol 2:166-173. 8. Cooper A (1997) Studies of avian ancient DNA: From Jurassic park to modern island extinctions. Avian Molecular Evolution and Systematics, ed Mindell DP (Academic, New York), pp 345-373. 9. Baker AJ, Huynen LJ, Haddrath O, Millar CD, Lambert DM (2005) Reconstructing thetempo and mode of evolution in an extinct clade of birds with ancient DNA: The giant moas of New Zealand. Proc Natl Acad Sci USA 102:8257-8262. 10. Bunce M, etal. (2003) Extreme reversed sexual size dimorphism in the extinct New Zealand moa Dinornis. Nature 425:172-175. 11. Cooper A, et al. (2001) Complete mitochondrial genome sequences of two extinct moas clarify ratite evolution. Nature 409:704-707. 12. Huynen L, Millar CD, Scofield RP, Lambert DM (2003) Nuclear DNA sequences detect species limits in ancient moa. Nature 425:175-178. 13. Lambert DM, et al. (2005) Is a large-scale DNA-based inventory of ancient life possible? J Hered 96:279-284. 14. Cooper A, etal. (1992) Independent origins of New Zealand moas and kiwis. Proc Natl Acad Sci USA 89:8741-8744. 15. GlazkoGV, Koonin EV, Rogozin IB (2005) Molecular dating: Ape bones agree with chicken entrails. Trends Gen 21:89-92. 16. Harrison GL, et al. (2004) Four new avian mitochondrial genomes help get to basic evolutionary questions in the late Cretaceous. Mol Biol Evol 21:974-983. 17. Ho SYW, Larson G (2006) Molecular clocks: When times are a-changin'. Trends Gen 22:79-83. 18. Ho SYW, Phillips MJ, Cooper A, Drummond AJ (2005) Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. Mol Biol Evol 22:1561-1568. 19. Haddrath O, Baker AJ (2001) Complete mitochondrial DNA genome sequences of extinct birds: Ratite phylogenetics and thevicariance biogeography hypothesis. Proc Roy Soc Lond 8 268:939-945. 20. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7:214. 21. Shapiro B, et al. (2004) Rise and fall of the Beringian steppe bison. Science 306:1561-1565. 22. Bourdon E, de Ricqles A, Cubo J (2009) A new Transantarctic relationship: Morphological evidence for a Rheidae-Dromaiidae-Casuariidae clade (Aves, Palaeognathae, Ratitae). ZoolJLinn Soc 156:641-663. 23. Moritz C, Cicero C (2004) DNA barcoding: Promise and pitfalls. P/oS Biol 2:1529-1531. 24. Hickerson MJ, Meyer CP, Moritz C (2006) DNA barcoding will often fail to discover new animal species over broad parameter space. Syst Biol 55:729-739. 25. Allentoft ME, et al. (2009) Identification of microsatellites from an extinct moa species using high-throughput (454) sequence data. Biotechniques 46:195-200. 26. Worthy TH, Edwards AR, Millener PR (1991) The fossil record of moas (Aves, Dinornithiformes) older than the Otira (last) glaciation. J Roy Soc NZ 21:101-118. 27. Pulford A, Stern TA (2004) Pliocene exhumation and landscape evolution of central North Island: The role of the upper mantle. J Geophys Res, 10.1029/2003JF000046. 28. Kamp PJJ, et al. (2004) Neogene Stratigraphic architecture and tectonic evolution of Wanganui, King Country, and eastern Taranaki Basins, New Zealand. NZ J Geol Geophys 47:625-644. 29. Kamp PJJ, Furlong KP (2006) Neogene plate tectonic reconstructions and geodynamics of North Island sedimentary basins. New Zealand Petroleum Conference Proceedings, p 16. 30. Field BD, Uruski CI (1997) Cretaceous - Cenozoic geology and petroleum systems of the East Coast Region. NZ Inst Geol Nucl Sci Mono 19:301. 31. Cooper A, Cooper RA (1995) The Oligocene bottleneck and New-Zealand biota - genetic record of a past environmental crisis. Proc Roy Soc Lond B 261:293-302. 32. Jones MEH, Tennyson AJD, Worthy JP, Evans SE, Worthy TH (2009) A sphenodontine (Rhynchocephalia) from the Miocene of New Zealand and palaeobiogeography of the tuatara (Sphenodon). Proc Roy Soc Lond B 276:1385-1390. 33. Worthy TH, Tennyson AJD, Jones C, McNamara J A, Douglas BJ (2007) Miocene waterfowl and other birds from central Otago, New Zealand. J Syst Palaeont 5:1-39. 34. Chappie DG, Ritchie PA, DaughertyCH (2009) Origin, diversification, and systematics of the New Zealand skink fauna (Reptilia: Scincidae). Mol Phylogen Evol 52:470-487. 35. Jones M, Tennyson A, Evans S, Worthy T (2008) The first pre-Pleistocene record of a Tuatara (Sphenodon)-like animal from New Zealand and implications for the Oligocene Drowning. J Vert Paleont 28:97A-98A. 36. Worthy TH, et al. (2009) Biogeographical and Phylogenetic Implications of an Early Miocene Wren (Aves: Passeriformes: Acanthisittidae)from New Zealand. J Vert Paleont, in press. 37. Knapp M, Mudaliar R, Havell D, Wagstaff SJ, Lockhart PJ (2007) The drowning of New Zealand and the problem of Agathis. Syst Biol 56:862-870. 38. Olsen SL, James HF (1991) Descriptions of thirty-two new species of birds from the Hawaiian Islands: Part 1. Non-Passeriformes. Ornithol Mono 45:1-88. 39. Coates G (2002) The Rise and Fall of the Southern Alps (Canterbury Univ Press, Christchurch, New Zealand), p 80. 40. Badgley C, et al. (2008) Ecological changes in Miocene mammalian record show impact of prolonged climatic forcing. Proc Natl Acad Sci USA 105:12145-12149. 41. Shevenell AE, Kennett JP, Lea DW (2004) Middle Miocene southern ocean cooling and Antarctic cryosphere expansion. Science 305:1766-1770. 42. Lewis KB, Carter L, Davey FJ (1994) The opening of Cook Strait - Interglacial tidal scour and aligning basins at a subduction to transform plate edge. Mar Geol 116:293-312. 43. Trewick SA (2001) Scree weta phylogeography: Surviving glaciation and implications for Pleistocene biogeography in New Zealand. New Zealand Journal of Zoology 28:291-298. 44. McGlone MS (1985) Plant biogeography and the late Cenozoic history of New Zealand. NZ J Bot 23:723-749. 45. Waugh J (2007) DNA barcoding in animal species: Progress, potential and pitfalls.fi/ofssays 29:188-197. 46. Posada D, Crandall KA (1998) MODELTEST: Testing the model of DNA substitution. Bioinform 14:817-818. 47. Sanderson MJ (2002) Estimating absolute rates of molecular evolution and divergence times: A penalized likelihood approach. Mol Biol Evol 19:101-109. Bunce et al. PNAS I December 8,2009 I vol.106 I no. 49 I 20651