- Research article
- Open Access
Complete mitochondrial genomes of living and extinct pigeons revise the timing of the columbiform radiation
BMC Evolutionary Biologyvolume 16, Article number: 230 (2016)
Pigeons and doves (Columbiformes) are one of the oldest and most diverse extant lineages of birds. However, the nature and timing of the group’s evolutionary radiation remains poorly resolved, despite recent advances in DNA sequencing and assembly and the growing database of pigeon mitochondrial genomes. One challenge has been to generate comparative data from the large number of extinct pigeon lineages, some of which are morphologically unique and therefore difficult to place in a phylogenetic context.
We used ancient DNA and next generation sequencing approaches to assemble complete mitochondrial genomes for eleven pigeons, including the extinct Ryukyu wood pigeon (Columba jouyi), the thick-billed ground dove (Alopecoenas salamonis), the spotted green pigeon (Caloenas maculata), the Rodrigues solitaire (Pezophaps solitaria), and the dodo (Raphus cucullatus). We used a Bayesian approach to infer the evolutionary relationships among 24 species of living and extinct pigeons and doves.
Our analyses indicate that the earliest radiation of the Columbidae crown group most likely occurred during the Oligocene, with continued divergence of major clades into the Miocene, suggesting that diversification within the Columbidae occurred more recently than has been reported previously.
The lineage of pigeons and doves, Columbiformes, is one of the most diverse non-passerine orders of birds. Columbiformes are the sixth most speciose order among the 40 traditionally recognized orders of living birds, according to species counts by the International Ornithologist’s Committee World Birdlist . Pigeons and doves inhabit every land area outside the Arctic and Antarctic, and display a wide range of variation in their ecological adaptations, although their relatively conserved anatomy and morphology has obscured phylogenetic relationships within the family. Recent whole genome analyses resolved the placement of pigeons and doves as sister to sandgrouses (Pterocliformes) and mesites (Mesitornithiformes) [2, 3].
Previous genetic analyses have helped to clarify cryptic relationships among some branches of the Columbiformes, and have provided insights into the timing of the diversification of this group [4–11]. For example, ancient DNA extracted from the remains of two large flightless pigeons, the extinct dodo (Raphus cucullatus) and its sister species, the solitaire (Pezophaps solitaria), suggested that the closest living relative of these species is the Nicobar pigeon (Caloenas nicobarica) . This work also suggested that the dodo and solitaire lineages diverged 18–36 Million years ago (Mya), during the late Oligocene . This date was biogeographically interesting because it was prior to the emergence of the two islands to which the flightless species were endemic . Similarly, old divergence estimates were obtained in a later study by Pereira et al. , who used a more taxonomically comprehensive phylogeny of both mitochondrial and nuclear DNA to infer that the dodo and solitaire diverged from the Caloenas lineage 33.5–50 Mya and from each other 15–30 Mya. This study also concluded that the entire columbiform lineage probably originated during the Cretaceous, and that the main period of diversification among columbiformes occurred at the Paleocene/Eocene boundary .
Here, we revisit the timing of the origin of and diversification within the columbiformes using a data set of complete mitochondrial genomes from a taxonomically broad selection of pigeons and doves. We assemble complete mitochondrial genomes from eleven pigeon species and from the yellow-throated sandgrouse, Pterocles gutturalis. Our new genomes include those of the extinct Ryukyu wood pigeon (Columba jouyi), the extremely rare thick-billed ground dove (Alopecoenas salamonis), which is known from only two specimens [14, 15], the spotted green pigeon (Caloenas maculata), which is also known as Liverpool pigeon and is represented by only one surviving museum specimen , the Rodrigues solitaire (Pezophaps solitaria), and the dodo (Raphus cucullatus). Using these newly assembled mitochondrial genomes and available published mitochondrial genomes, we estimate a phylogeny to infer the major evolutionary relationships among the pigeons and doves using both a Bayesian and maximum likelihood approach. Further, we use a molecular clock approach, calibrated using whole-genome data [2, 3], to infer the timing of divergence between Columbiformes, Pterocliformes, and the Galloansera.
Results and discussion
A new phylogeny for pigeons and doves
Both the ML and Bayesian approaches to inferring a mitochondrial phylogeny result in the same overall topology, with strong statistical support for most nodes (Fig. 1). The phylogeny supports two major clades, an Indo-Pacific clade (Fig. 1, yellow), and a Holarctic clade (Fig. 1, blue) that also includes New World pigeons. This or similar structure has been observed previously in more taxonomically focused data sets [5, 17, 18].
Within the Indo-Pacific clade, our results corroborate several previously supported relationships. Similarly to Jønsson et al.  and Moyle et al. , we find that Alopecoenas is more closely related to the Zebra dove (Geopelia striata) than it is to the Luzon bleeding-heart dove (Gallicolumba luzonica). Our results also support the previously identified close relationship between the dodo (Raphus cucullatus), solitaire (Pezophaps solitaria), and Nicobar pigeon (Caloenas nicobarica) , and the sister relationship between the Nicobar pigeon and the spotted green pigeon (Caloenas maculata) . We also find that this group is sister to the crowned pigeons Goura, supporting the result of Shapiro et al. , but in contrast to Pereira et al. , who placed the Caloenas/Raphus lineage within the Didunculinae clade.
Also within the Indo-Pacific clade, our results suggest that the genus Otidiphaps, instead of belonging to the Didunculinae clade, is sister to the highly diverse fruit pigeon clade , which is represented in our phylogeny by Hemiphaga. Further investigation that includes other species that have been suggested to be closely related to these genera, such as the thick-billed ground pigeon (Trugon terrestris) , will help to further disentangle these evolutionary relationships.
Within the Holarctic and New World clade, we identify a strongly supported subclade that includes the genera Leptotila, Zenaida and Geotrygon. These species occur from North to South America and diverged from each other during the early Miocene. The close relationship between Zenaida and Leptotila conflicts with previous results  that placed the white-tipped dove, Leptotila verreauxi, closer to the violaceous quail-dove, Geotrygon violacea, than to Zenaida. The difference between our and previous phylogenies may be attributable to including the mourning dove, Zenaida macroura, in our analysis. The Leptotila/Zenaida/Geotrygon clade has strong support in the Bayesian analysis, but weak support in the Maximum Likelihood tree, highlighting the challenge of inferring and interpreting phylogenetic relationships from taxonomically limited data sets.
Inferring the timing of diversification within pigeons and doves
The combination of recent genome-scale analyses of avian evolution  and our new data set of complete mitochondrial genomes provides an opportunity to recalibrate the timing of the origin of and radiation within the Columbiformes. When inferring time-calibrated phylogenies, careful consideration is required with respect to each fossil or type of calibration employed . Theoretical and empirical work have shown that using multiple calibration points generally leads to more robust estimates of evolutionary rates [21, 22]. Unfortunately, no fossils are known from within the family of pigeons and doves that could be used as calibration [23, 24]. Therefore, we used the time-scales estimated by Jarvis et al.  and Prum et al. , which agree with each other with respect to the timing of diversification of Columbiformes.
We find that pigeons and doves most likely began to diversify during the late Oligocene, and continued to diversify into the Miocene (Fig. 1). Specifically, we find that the Holarctic and Indo-Pacific clades diverge around 24.7 Mya (95 % credibility interval, CI, 18.9–31.3 Mya), similar to , which place the divergence of two columbiform clades during the mid-Oligocene. This timing is in contrast to previous studies, which suggested that the Columbiformes radiated much earlier and more slowly, over the course of the Eocene and Oligocene . The transition from the Eocene to the Oligocene corresponds to a period of widespread global cooling [26, 27] and associated geological changes, including the opening of the Drake Passage and the formation of the Wallacea region . These changing global conditions may in part explain the timing of the rapid diversification within the Columbidae, which contains many highly dispersive, “supertramp” species .
We estimate that the dodo and solitaire diverged around 13.1 Mya (95 % CI 9.5–17.3 Mya), during the Early/Middle Miocene transition, rather than around the Oligocene/Miocene transition (22.8 Mya  and 25.6 Mya ), as previously proposed. We also find that the common ancestor of the dodo and solitaire diverged from Caloenas around 18 Mya (95 % CI 13.6–23.2 Mya), rather than during the Middle or latest Eocene (33.6 Mya  and 42.6 Mya ). Although our estimated divergence dates are more recent than those proposed previously, these dates indicate that both flightless pigeons diverged from their closest flying relative at least 10 Mya prior to the emergence of Mauritius and Rodrigues Islands, to which the flightless birds were endemic [12, 30, 31]. This finding corroborates previous claims that these lineages must have originated elsewhere . The passenger pigeon, Ectopistes migratorius, is known to be closely related to the lineage of large New Word pigeons, represented in our study by the band-tailed pigeon, Patagioenas fasciata [10, 11]. However, the timing of divergence between these lineages has been unknown. Our phylogeny indicates that passenger pigeons and band-tailed pigeons share a common ancestor around 12.4 Mya (95 % CI 9.0–16.3 Mya). This common ancestor diverged from other Old World pigeons during the transition between the Oligocene and the Miocene, around 16.2 Mya (95 % CI 11.7–20.5 Mya). This result contrasts with the results of Pereira et al. , which placed the split of Patagioenas/Ectopistes and the remaining Columbids around 35 Mya.
Rapid diversification, such as that identified here for the Columbiformes, may lead to variation among gene trees due to the effects of incomplete lineage sorting [32, 33], which can lead to inference of different phylogenies for different loci. In future, therefore, it will be important to confirm the evolutionary hypotheses presented here using multiple, independently inherited markers. The use of additional calibration points will also likely increase the precision of the nodes age estimates. Nonetheless, the strong support for the branching order of our phylogeny provides new insights into many of the cryptic evolutionary relationships among pigeons and doves, and attests to the resolving power of complete mitochondrial genomes, at least for inference of the evolutionary history of this locus. Broader taxonomic sampling and the addition of a greater diversity of extinct lineages and calibration points may further resolve the timing and nature of evolutionary diversification within this very diverse group of birds.
We present a new phylogeny of the pigeons and doves (Columbiformes) based on complete mitochondrial genomes from 24 pigeon species including several extinct or extremely rare species. The branching order in the phylogenetic tree is strongly statistically supported. By including a molecular rate calibration from recent genome-scale analyses, we infer that the lineage of pigeons and doves began to diversify during the Oligocene/Miocene transition, which is a more recent diversification than has been suggested previously. We interpret the phylogenetic results in the context of previous research, and support the recognition of the genus Alopecoenas.
DNA extraction and sequencing
We obtained recent or historic tissues for 16 pigeon and dove specimens for the purposes of generating mitochondrial genomes. This included bone powder for the dodo and solitaire, feather for the spotted green pigeon and tooth-billed pigeon, and toe pads for all other samples (Table 1). For modern samples we used muscle tissue for mourning dove, muscle and blood tissue for band-tailed pigeon, and liver tissue for the Nicobar pigeon.
We processed modern tissue samples at two institutions: the UCSC Paleogenomics Lab Modern Facility (N2009-0024, BTP2013), and the Church Lab, at Harvard University (AMNH DOT 14025, Zm1). For all the recent samples, we extracted DNA using the Qiagen DNeasy Blood & Tissues Kit (Qiagen, USA) according to the manufacturer’s instructions and sheared the resulting DNA into fragments <1,000 base pairs (bp) long. We size-selected DNA from the Nicobar Pigeon and band-tailed pigeon prior to sequencing. Mourning dove DNA was sequenced without size selection.
For historic samples, we extracted DNA and prepared genomic libraries in isolated, purpose-built ancient DNA facilities at four institutions: the UCSC Paleogenomics ab (samples AMNH 612456, OMNH 1764, AMNH 224546, OMNH 1762, AMNH 616460, Zm1, S1B1, OUMNH 1759), the University of Copenhagen (sample ZMUC AVES-105485), the McMaster University Ancient DNA Centre (samples FMNH 47395, FMNH 47396, and FMNH 47397), and Griffith University (sample D3538). At UCSC, we extracted DNA following , in which we digested tissues in buffer modified from the Qiagen Blood & Tissue Kit that comprised 150 μL Buffer ATL, 30 μL proteinase K solution, and 20 μL of 1 M dithiothreitol (DTT), in a rotation incubator at 56 °C for 48 h, and then purified DNA using the Qiagen Nucleotide Removal Kit according to the manufacturers protocol. At McMaster University, we extracted DNA using a phenol:chloroform:isoamyl alcohol and chloroform based solution, or “in-house” silica columns with an extraction to binding buffer ratio of 1:2 and 30 μL silica beads [34, 35]. At Copenhagen University, we first drilled 0.01 g of bone powder through the Foramen Magnum from the inside of the braincase of the “Gottorp” Dodo Specimen (ZMUC AVES-105485), not damaging the exterior of the skull, and using appropriate anti-contamination precautions. We then extracted DNA from the bone powder as in , and purified the extract following . We then constructed an Illumina sequencing library using a blunt-end protocol . The library was indexed by amplification with Accuprime under the following conditions: 95 °C for 1 min, then 12 cycles of 95 °C for 15 s, 60 °C for 30 s, 68 °C for 30 s, followed by 68 °C for seven minutes. At Griffith University, we extracted DNA from a feather of the single surviving spotted green pigeon specimen following . The resulting DNA was treated with Uracil-DNA Glycosylase (Thermo), the NEBNext End Repair Module (NEB), the NEBNext Quick Ligation Module (NEB) and the Bst DNA Polymerase Large Fragment (NEB), respectively, each time purifying with the MinElute PCR Purification Kit (Qiagen). We prepared sequencing libraries according to  and cleaned the resulting libraries with the AxyPrep Mag PCR Clean-Up Kit (Axygen). The libraries were sequenced as 50 bp single-end reads on an Illumina HiSeq 2500 sequencing system at the Macrogen Inc, South Korea. For samples AMNH DOT 14025 and Zm1 DNA extracts were prepared for sequencing using the Illumina TruSeq kit following manufacturer’s protocols.
For samples extracted at McMaster and UCSC, we prepared uniquely-barcoded Illumina sequencing libraries following  and cleaned the libraries using Sera-Mag SPRI SpeedBeads (ThermoScientific) in 18 % PEG-8000. We generated paired-end sequence data from pooled libraries using both an Illumina MiSeq (1 × 75 bp) at UCSC and Illumina HiSeq2000 (1 × 100 bp) at the Vincent J. Coates Genomic Sequencing Center at UC Berkeley, the University of Copenhagen, McMaster University and the University of Toronto.
Assembly of mitochondrial genomes
To assemble a taxonomically diverse data set of pigeons and doves, we downloaded short read files (SRA) from the online database Sequence Read Archive + (http://sra.dnanexus.com) for Columba rupestris (accession SRS346866 ), Pterocles gutturalis (accession SRP029347 [40, 41]) and Ectopistes migratorius (accession SRS391366). We assembled the mitochondrial genomes of passenger pigeon specimens FMNH 47395 and FMNH 47397 (SRA accession SRS391366 for both specimens, GenBank accession JQ692598 for FMNH 47397). We assembled a new mitochondrial genome for FMNH 47397 because the previous assembly contains improperly assembled regions (16 s rRNA, ND6, and D-loop).
We processed our sequencing data and the previously published SRAs by removing adapters using SeqPrep (https://github.com/jstjohn/SeqPrep). For the extinct species, sequence fragments were sufficiently short that we also merged the paired reads also using SeqPrep, enforcing a minimum overlap of 10 base-pairs between forward and reverse reads. We mapped the processed reads to the reference mitochondrial genome of Columba livia (GenBank accession NC_013978.1 ), and other published pigeon mitochondrial genomes using MIA (https://github.com/udo-stenzel/mapping-iterative-assembler), which is an iterative, reference-based, short-fragment assembler designed for circular genomes . For each genome, we required a minimum of three unique molecules (3× coverage) to call a consensus base at each site; otherwise, bases were called as “N.” Finally, we inspected the assemblies by eye using Geneious R8.1 and corrected poorly assembled regions with a second set of iterative mapping assemblies (http://www.geneious.com, ).
In addition to the sequences described above, we downloaded previously published mitochondrial data for 12 species of pigeon and dove and chicken (Table 1). We also concatenated mitochondrial genes 16S, 12S, ATP8, ATP6, CYTB, COIII, COII, COI, ND5, ND4, ND4L, ND3, ND2, ND1 from the Namaqua sandgrouse, creating a partial mitochondrial genome (GenBank accessions DQ385063.1, DQ385080.1, DQ385097.1, DQ385114.1, DQ385131.1, DQ385148.1, DQ385165.1, DQ385182.1, DQ385199.1, DQ385216.1, DQ385233.1, DQ385250.1, DQ385267.1, DQ385284 ). We aligned all the full mitochondrial genomes and the concatenated sandgrouse mitochondrial genes using MUSCLE , and visually checked the alignment using SeaView v.4.5.4 .
We partitioned the alignment in order to account for different evolutionary rates along the different regions of the mitochondrial genome. We split the alignment into six distinct partitions: 12S and 16S ribosomal RNA genes, all tRNA genes, the hypervariable region, and three partitions for the protein coding genes, according to their codon position. We used PartitionFinder  to select the models of evolutionary evolution for each partition and estimated the phylogenetic relationships among the mitochondrial genomes in our data set using BEAST 1.8.1  and RAxML v.8.2.0 . For BEAST we assumed a lognormal uncorrelated relaxed clock  for each partition, and a Birth-Death speciation process  for the tree prior. To calibrate the molecular clock, we placed a normal prior on the age of the divergence between Columbiformes and Pterocliformes of 55 Mya ± 15 Mya [2, 3]. We ran six MCMC chains for 30 million states, sampling trees and model parameters every 3000 states. We discarded the first 30 % as burn-in, and visually inspected the remainder for convergence using Tracer v1.6 .
Since 3rd codon positions evolve more rapidly than 1st and 2nd codon positions, we tested for evidence of saturation at these sites using the method described in . We found no evidence of saturation (Additional file 1: Figure S3). Nevertheless, since changes in substitution biases can result in systematic error , and this is more likely to affect synonymous substitutions, we repeated the BEAST analysis using only 1st and 2nd codon positions of coding genes. We obtained the same topology entirely consistent with the use of the entire alignment, suggesting that our estimation of the topology using the whole alignment is not driven by systematic error.
If rates are too variable over the history of a particular group, dating analyses will be unreliable. The absence of multiple calibration points limits our ability to estimate the extent to which rates may have varied over the history of this clade. However, we found that the standard deviation (SD) in the rate variation over the tree, as obtained by BEAST, is sufficiently small to suggest a good fit of a molecular clock model (8.627E-4 SD). Despite this, the determinants of rate variation might be heritable , and may, therefore, be better reflected by an autocorrelated clock model. Therefore, in addition to the uncorrelated relaxed clock method from BEAST, we also ran a dating analysis with an autocorrelated rate model using MCMCTREE . We used the same calibration as used for the BEAST analysis, and the same partitioned dataset. The 95 % HPD (highest posterior density) for estimates of node age obtained from the two Bayesian approaches all overlap and are reported in Additional file 2: Table S1 and Additional file 3: Figure S2.
We performed a maximum likelihood (ML) phylogenetic analysis on the same data set using RAxML v.8.2.0 . We assumed the GTRGAMMA model for each partition, and performed 1000 bootstrap replicates.
To evaluate the robustness of our molecular clock approach, we performed an additional analysis in which we estimated divergence times at the nodes in our phylogeny using Reltime . Reltime uses a maximum likelihood approach to calculate branch lengths in substitutions per site and computes branch-specific relative rates without calibrations. We used the topology generated by BEAST as input tree and a GTR evolutionary model with four gamma categories. We then converted the relative divergence time into absolute times using the same calibration as used in the BEAST analysis, except that Reltime treats the calibration information as a flat time range. All Reltime calculations were done in MEGA7 . The ages estimated using this approach fell within the 95 % HPD of estimates from the BEAST and MCMCTREE analyses (Additional file 1: Table S1, Additional file 2: Figure S2), with slightly older mean results compared to BEAST. Based on these results, we infer that while additional calibration points may increase precision in node age estimates, their inclusion is not likely to alter the shape of the tree topology significantly.
12S ribosomal RNA
16S ribosomal RNA
ATP synthase subunit 6
ATP synthase subunit 8
Cytochrome oxidase subunit I
Cytochrome oxidase subunit II
Cytochrome oxidase subunit III
Million years ago
- ND1 NADH:
Ubiquinone oxidoreductase core subunit 1
Ubiquinone oxidoreductase core subunit 2
- ND3 NADH:
Ubiquinone oxidoreductase core subunit 3
- ND4 NADH:
Ubiquinone oxidoreductase core subunit 4
- ND4L NADH:
Ubiquinone oxidoreductase core subunit 4 L
- ND5 NADH:
Ubiquinone oxidoreductase core subunit 5
Sequence read archive
Gill FB, Donsker DB. In: Gill FB, Donsker DB, editors. IOC World Bird List (v 5.2) [Internet]. 2015. Available from: http://www.worldbirdnames.org/ioc-lists/crossref/.
Jarvis ED, Mirarab S, Aberer AJ, Li B, Houde P, Li C, et al. Whole-genome analyses resolve early branches in the tree of life of modern birds. 2014. Available from: www.sciencemag.org/content/346/6215/1320/suppl/DC1.
Prum RO, Berv JS, Dornburg A, Field DJ, Townsend JP, Lemmon EM, et al. A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing. Nature [Internet]. 2015;526:569–73. Nature Publishing Group, a division of Macmillan Publishers Limited. Available from: http://dx.doi.org/10.1038/nature15697.
Jønsson KA, Irestedt M, Bowie RCK, Christidis L, Fjeldså J. Systematics and biogeography of Indo-Pacific ground-doves. Mol Phylogenet Evol. 2011;59(2):538–43.
Pereira SL, Johnson KP, Clayton DH, Baker AJ. Mitochondrial and nuclear DNA sequences support a Cretaceous origin of Columbiformes and a dispersal-driven radiation in the Paleocene. Syst Biol. 2007;56:656–72.
Johnson KP, Clayton DH. A molecular phylogeny of the dove genus zenaida: mitochondrial and nuclear dna sequences. Condor [Internet]. 2000;102:864. Available from: http://aoucospubs.org/doi/abs/10.1650/0010-5422(2000)102%5B0864%3AAMPOTD%5D2.0.CO%3B2 [cited 2015 Jun 18].
Johnson KP, de Kort S, Dinwoodey K, Mateman AC, ten Cate C, Lessells CM, et al. A molecular phylogeny of the dove Genera Streptopelia and Columba. Auk [Internet]. 2001;118:874. Available from: http://aoucospubs.org/doi/abs/10.1642/0004-8038(2001)118%5B0874%3AAMPOTD%5D2.0.CO%3B2 [cited 2015 Jun 18].
Cibois A, Thibault JC, Bonillo C, Filardi CE, Watling D, Pasquet E. Phylogeny and biogeography of the fruit doves (Aves: Columbidae). Mol Phylogenet Evol. 2014;70:442–53.
Moyle RG, Jones RM, Andersen MJ. A reconsideration of Gallicolumba (Aves: Columbidae) relationships using fresh source material reveals pseudogenes, chimeras, and a novel phylogenetic hypothesis. Mol Phylogenet Evol. 2013;66(3):1060–6.
Johnson KP, Clayton DH, Dumbacher JP, Fleischer RC. The flight of the passenger pigeon: phylogenetics and biogeographic history of an extinct species. Mol Phylogenet Evol. 2010;57(1):455–8.
Fulton TL, Wagner SM, Fisher C, Shapiro B. Nuclear DNA from the extinct Passenger Pigeon (Ectopistes migratorius) confirms a single origin of New World pigeons. Ann Anat [Internet]. 2012;194:52–7. Elsevier GmbH. Available from: http://dx.doi.org/10.1016/j.aanat.2011.02.017.
McDougall I, Chamalaun FH. Isotopic Dating and Geomagnetic Polarity Studies on Volcanic Rocks from Mauritius, Indian Ocean. Geol Soc Am Bull [Internet]. 1969;80:1419–42. Available from: http://gsabulletin.gsapubs.org/content/80/8/1419.short.
Shapiro B, Sibthorpe D, Rambaut A, Austin J, Wragg GM, Bininda-Emonds ORP, et al. Flight of the dodo. Science. 2002.
Gibbs D, Barnes E, Cox J. Pigeons and Doves: A guide to the pigeons and doves of the world [Internet]. 1st ed. New Haven: Yale University Press; 2001.
Danielsen F, Filardi CE, Jønsson KA, Kohaia V, Krabbe N, Kristensen JB, et al. Endemic avifaunal biodiversity and tropical forest loss in Makira, a mountainous Pacific island. Singap. J Trop Geogr. 2010;30:100–14.
Heupink TH, van Grouw H, Lambert DM. The mysterious Spotted Green Pigeon and its relation to the Dodo and its kindred. BMC Evol Biol [Internet]. 2014;14:136. Available from: http://bmcevolbiol.biomedcentral.com/articles/10.1186/1471-2148-14-136 [cited 2 May 2016].
Gibb GC, Penny D. Two aspects along the continuum of pigeon evolution: a South-Pacific radiation and the relationship of pigeons within Neoaves. Mol Phylogenet Evol [Internet]. 2010;56:698–706. Available from: http://dx.doi.org/10.1016/j.ympev.2010.04.016.
Fulton TL, Wagner SM, Fisher C, Shapiro B. Nuclear DNA from the extinct passenger pigeon (Ectopistes migratorius) confirms a single origin of new world pigeons. Ann Anat. 2012;194(1):52–7.
Pacheco MA, Battistuzzi FU, Lentino M, Aguilar RF, Kumar S, Escalante AA. Evolution of modern birds revealed by mitogenomics: timing the radiation and origin of major orders. Mol Biol Evol. 2011;28(6):1927–42.
Warnock RCM, Parham JF, Joyce WG, Lyson TR, Donoghue PCJ. Calibration uncertainty in molecular dating analyses: there is no substitute for the prior evaluation of time priors. Proc Biol Sci [Internet]. 2015;282:20141013. Available from: http://rspb.royalsocietypublishing.org/content/282/1798/20141013 [cited 4 May 2016].
Near TJ, Sanderson MJ. Assessing the quality of molecular divergence time estimates by fossil calibrations and fossil-based model selection. Philos Trans R Soc Lond B Biol Sci [Internet]. 2004;359:1477–83. Available from: http://rstb.royalsocietypublishing.org/content/359/1450/1477 [cited 4 May 2016].
Graur D, Martin W. Reading the entrails of chickens: molecular timescales of evolution and the illusion of precision. Trends Genet [Internet]. 2004;20:80–6. Elsevier. Available from: http://www.cell.com/article/S0168952503003421/fulltext [cited 30 Mar 2016].
Worthy TH, Hand SJ, Worthy JP, Tennyson AJD, Scofield RP. A large fruit pigeon (Columbidae) from the early Miocene of New Zealand. Auk [Internet]. 2009;126:649–56. Available from: http://www.bioone.org/doi/abs/10.1525/auk.2009.08244 [cited 4 May 2016].
Mayr G. The paleogene fossil record of birds in Europe. Biol Rev Camb Philos Soc [Internet]. 2005;80:515–42. Available from: http://doi.wiley.com/10.1017/S1464793105006779 [cited 18 Sep 2015].
Claramunt S, Cracraft J. A new time tree reveals earth historys imprint on the evolution of modern birds. Sci Adv [Internet]. 2015;1:e1501005. American Association for the Advancement of Science. Available from: http://advances.sciencemag.org/content/1/11/e1501005.abstract [cited 12 Dec 2015].
Zachos J, Pagani M, Sloan L, Thomas E, Billups K. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science [Internet]. 2001;292:686–93. Available from: http://www.sciencemag.org/content/292/5517/686.abstract [cited 2 Feb 2015].
Bohaty SM, Zachos JC, Delaney ML. Foraminiferal Mg/Ca evidence for southern ocean cooling across the eocene–oligocene transition. Earth Planet Sci Lett [Internet]. 2012;317–318:251–61. Available from: http://www.sciencedirect.com/science/article/pii/S0012821X11007059 [cited 18 Sep 2015].
Hall R. Southeast Asia’s changing palaeogeography. Blumea - biodiversity. Evol Biogeogr Plants [Internet]. 2009;54:148–61. Nationaal Herbarium Nederland. Available from: http://www.ingentaconnect.com/content/nhn/blumea/2009/00000054/F0030001/art00026?token=004f12494d7e2a46762c6b665d58663f257023796d42673f7b2f267738703375686f49bee7cea2a cited [19 Sep 2015].
Whittaker RJ, Fernández-Palacios JM. Island biogeography: ecology, evolution, and conservation. Oxford: Oxford University Press; 2007.
Saddul P. Mauritius - a geomorphological analysis. Mauritius: Mahatma Gandhi Institute; 1995.
Cheke A, Hume JP. Lost Land of the Dodo: The Ecological History of Mauritius, Réunion, and Rodrigues. New Haven: Yale University Press; 2008.
Rokas A, Krüger D, Carroll SB. Animal evolution and the molecular signature of radiations compressed in time. Science [Internet]. 2005;310:1933–8. Available from: http://www.sciencemag.org/content/310/5756/1933.abstract [cited 19 Jun 2015].
Glor RE. Phylogenetic insights on adaptive radiation. Annu Rev Ecol Evol Syst [Internet]. 2010;41:251–70. Annual Reviews. Available from: http://www.annualreviews.org/doi/abs/10.1146/annurev.ecolsys.39.110707.173447 [cited 19 Jun 2015].
Enk J, Devault A, Debruyne R, King CE, Treangen T, O’Rourke D, et al. Complete Columbian mammoth mitogenome suggests interbreeding with woolly mammoths. Genome Biol [Internet]. 2011;12:R51. Available from: http://genomebiology.com/2011/12/5/R51 cited 27 May 2015].
Rohland N, Siedel H, Hofreiter M. A rapid column-based ancient DNA extraction method for increased sample throughput. Mol Ecol Resour [Internet]. 2010;10:677–83. http://www.ncbi.nlm.nih.gov/pubmed/21565072 [cited 27 May 2015].
Shapiro B, Drummond AJ, Rambaut A, Wilson MC, Matheus PE, Sher AV, et al. Rise and fall of the Beringian steppe bison. Science [Internet]. 2004;306:1561–5. Available from: http://www.sciencemag.org/content/306/5701/1561.full [cited 13 Jul 2015].
Yang DY, Eng B, Waye JS, Dudar JC, Saunders SR. Technical note: improved DNA extraction from ancient bones using silica-based spin columns. Am J Phys Anthropol [Internet]. 1998;105:539–43. Available from: http://www.ncbi.nlm.nih.gov/pubmed/9584894 [cited 5 Aug 2015].
Meyer M, Kircher M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb Protoc [Internet]. 2010;2010:pdb.prot5448. Available from: http://cshprotocols.cshlp.org/content/2010/6/pdb.prot5448.abstract [cited 11 Jul 2014].
Cui J, Zhao W, Huang Z, Jarvis ED, Gilbert MTP, Walker PJ, et al. Low frequency of paleoviral infiltration across the avian phylogeny. Genome Biol [Internet]. 2014;15:539. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=4272516&tool=pmcentrez&rendertype=abstract [cited 22 May 2015].
Zhang G, Li B, Li C, Gilbert MTP, Jarvis ED, Wang J. Comparative genomic data of the avian phylogenomics project. Gigascience [Internet]. 2014;3:26. BioMed Central. Available from: http://gigascience.biomedcentral.com/articles/10.1186/2047-217X-3-26 [cited 3 May 2016].
Zhang G, Li C, Li Q, Li B, Larkin DM, Lee C, et al. Comparative genomics reveals insights into avian genome evolution and adaptation. Science (80-) [Internet]. 2014;346:1311–20. Available from: http://www.sciencemag.org/content/346/6215/1311 [cited 11 Dec 2014].
Kan XZ, Li XF, Zhang LQ, Chen L, Qian CJ, Zhang XW, et al. Characterization of the complete mitochondrial genome of the Rock pigeon, Columba livia (Columbiformes: Columbidae). Genet Mol Res [Internet]. 2010;9:1234–49. Available from: http://www.ncbi.nlm.nih.gov/pubmed/20603809 [cited 5 Aug 2015].
Green RE, Malaspinas A-S, Krause J, Briggs AW, Johnson PLF, Uhler C, et al. A complete neandertal mitochondrial genome sequence determined by high-throughput sequencing. Cell [Internet]. 2008;134:416–26. Available from: http://linkinghub.elsevier.com/retrieve/pii/S0092867408007733.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics [Internet]. 2012;28:1647–9. Available from: http://bioinformatics.oxfordjournals.org/content/28/12/1647.long [cited 10 Jul 2014].
Paton TA, Baker AJ. Sequences from 14 mitochondrial genes provide a well-supported phylogeny of the Charadriiform birds congruent with the nuclear RAG-1 tree. Mol Phylogenet Evol. 2006;39(3):657–67.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res [Internet]. 2004;32:1792–7. Available from: http://nar.oxfordjournals.org/content/32/5/1792.long [cited 11 Jul 2014].
Gouy M, Guindon S, Gascuel O. SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol [Internet]. 2010;27:221–4. Available from: http://mbe.oxfordjournals.org/content/27/2/221.full [cited 12 Jul 2014].
Lanfear R, Calcott B, Ho SYW, Guindon S. Partitionfinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol [Internet]. 2012;29:1695–701. Available from: http://mbe.oxfordjournals.org/content/29/6/1695 [cited 4 Jun 2015].
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol [Internet]. 2012;29:1969–73. Available from: http://mbe.oxfordjournals.org/cgi/doi/10.1093/molbev/mss075.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics [Internet]. 2014;30:1312–3. Available from: http://bioinformatics.oxfordjournals.org/content/30/9/1312 [cited 15 Jul 2014].
Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol [Internet]. 2006;4:e88. Available from: http://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.0040088 [cited 9 Jul 2014].
Gernhard T. The conditioned reconstructed process. J Theor Biol [Internet]. 2008;253:769–78. Available from: http://www.sciencedirect.com/science/article/pii/S0022519308001811.
Rambaut A, Drummond AJ. Tracer V1.6 [Internet]. 2013. Available from: http://tree.bio.ed.ac.uk/software/tracer/.
Philippe H, Sörhannus U, Baroin A, Perasso R, Gasse F, Adoutte A. Comparison of molecular and paleontological data in diatoms suggests a major gap in the fossil record. J Evol Biol. 1994;7:247–65.
Philippe H. Opinion: long branch attraction and protist phylogeny. Protist [Internet]. 2000;151:307–16. Urban & Fischer. Available from: http://linkinghub.elsevier.com/retrieve/pii/S1434461004700292.
Nabholz B, Glémin S, Galtier N, Ballard J, Whitlock M, Lane N, et al. The erratic mitochondrial clock: variations of mutation rate, not population size, affect mtDNA diversity across birds and mammals. BMC Evol Biol [Internet]. 2009;9:54. BioMed Central. Available from: http://bmcevolbiol.biomedcentral.com/articles/10.1186/1471-2148-9-54 [cited 15 Jul 2016].
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.
Tamura K, Battistuzzi FU, Billing-Ross P, Murillo O, Filipski A, Kumar S. Estimating divergence times in large molecular phylogenies. Proc Natl Acad Sci U S A [Internet]. 2012;109:19333–8. Available from: http://www.pnas.org/content/109/47/19333.abstract [cited 3May 2016].
Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. [Internet]. 2016;msw054. Available from: http://mbe.oxfordjournals.org/content/33/7/1870 [cited 23 Mar 2016].
Jang KH, Ryu SH, Kang S-G, Hwang UW. Complete mitochondrial genome of the Japanese wood pigeon, Columba janthina janthina (Columbiformes, Columbidae). Mitochondrial DNA [Internet]. 2014; 1–2. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25431823 [cited 5 Aug 2015].
Wu H, Liu B-G, Hu G-Z, Liu J-H, Yuan L, Pan Y-S. The complete mitochondrial genome of Archangel pigeon. Mitochondrial DNA [Internet]. 2014; 1–2. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25007353 [cited 5 Aug 2015].
Nishibori M, Shimogiri T, Hayashi T, Yasue H. Molecular evidence for hybridization of species in the genus Gallus except for Gallus varius. Anim Genet [Internet]. 2005;36:367–75. Available from: http://www.ncbi.nlm.nih.gov/pubmed/16167978 [cited 5 Aug 2015].
Yan SQ, Guo PC, Li YM, Qi SM, Bai CY, Zhao ZH, et al. Complete mitochondrial genome of the Spotted dove (Streptopelia chinensis). Mitochondrial DNA [Internet]. 2015; 1–2. Available from: http://www.ncbi.nlm.nih.gov/pubmed/25600734 [cited 5 Aug 2015].
The authors acknowledge the contributions of the ornithology collections of the American Museum of Natural History, courtesy of Paul Sweet and Thomas J. Trombone, and the Field Museum of Natural History, courtesy of Dave Willard, for providing Ectopistes, Columba, Gallicolumba, Patagioenas, and Alopecoenas archival tissue samples. We thank Malgosia Nowak-Kemp and the Oxford University Museum of Natural History for providing access to Goura and Didunculus archival tissues. We thank World Museum (National Museums Liverpool), Clemency Fisher, Hein van Grouw, and Tony Parker for providing samples from the Caloenas maculata specimen. We thank Sal Alvarez of Exotic Wings International aviaries for supplying a band-tailed pigeon blood sample. We thank Lily Shiue, Steven Salzberg, Joshua Kapp, Steven Weber, Daniela Puiu, Beatriz Mello, and Gemma Murray for their contributions in data production and technical support. We also thank the reviewers and editor for the comments and suggestions.
AERS was funded by the Ciência sem Fronteiras Fellowship (CAPES, Brazil). BJN was funded by the Founding Funders, Angel Funders, and many supporters of Revive & Restore (http://longnow.org/revive/our-supporters/). This work was supported in part by the Gordon and Betty Moore Foundation and the Packard Foundation. Additional private funds towards sequencing were donated by James Sartor, Sarahí Avelar Aguiñaga, Janette and Anton J. Novak Jr., Walta and Timothy J. Novak, Marietta Koppang, Carmen and Nathon Maxwell.
Availability of data and materials
DNA sequences GenBank accessions KX902235–KX902250.
AERS, BJN, JH, THH, and BS carried out the molecular lab work, AERS, BJN, and BS. designed the study, AERS and BJN performed the analyses. AERS, BJN, and BS wrote the paper with contributions from JH, THH and JF, and all authors contributed to the final draft. THH, MTPG, HP, GMC, and BS also contributed reagents. BJN painted the images for Raphus cucullatus and Pezophaps solitaria on Fig. 1. All authors contributed to the discussion of results and gave their final approval for publication.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Time tree obtained using the BEAST method. Each branch is colored according to the gradient on the top left of the figure, from lower to higher rate estimates for the coding genes. Each branch has been labelled with its rate. (PDF 185 kb)
Divergence times estimates from BEAST, MCMCTREE and Reltime, including minimum and maximum 95 % CI. Node numbers refer to Additional file 1: Figure S1. (XLSX 10 kb)
A) Time tree obtained using the Reltime method. Each node received a number. B) Dated nodes, including 95 % CI. Brown circles denotes BEAST results, yellow squares MCMCTREE results, and blue circles Reltime results. Node ID relates to node numbers in panel A, and can be seen in table format in Additional file 1: Table S1. (PDF 250 kb)