Research article | Open | Published:
A North American stem turaco, and the complex biogeographic history of modern birds
BMC Evolutionary Biologyvolume 18, Article number: 102 (2018)
Earth’s lower latitudes boast the majority of extant avian species-level and higher-order diversity, with many deeply diverging clades restricted to vestiges of Gondwana. However, palaeontological analyses reveal that many avian crown clades with restricted extant distributions had stem group relatives in very different parts of the world.
Our phylogenetic analyses support the enigmatic fossil bird Foro panarium Olson 1992 from the early Eocene (Wasatchian) of Wyoming as a stem turaco (Neornithes: Pan-Musophagidae), a clade that is presently endemic to sub-Saharan Africa. Our analyses offer the first well-supported evidence for a stem musophagid (and therefore a useful fossil calibration for avian molecular divergence analyses), and reveal surprising new information on the early morphology and biogeography of this clade. Total-clade Musophagidae is identified as a potential participant in dispersal via the recently proposed ‘North American Gateway’ during the Palaeogene, and new biogeographic analyses illustrate the importance of the fossil record in revealing the complex historical biogeography of crown birds across geological timescales.
In the Palaeogene, total-clade Musophagidae was distributed well outside the range of crown Musophagidae in the present day. This observation is consistent with similar biogeographic observations for numerous other modern bird clades, illustrating shortcomings of historical biogeographic analyses that do not incorporate information from the avian fossil record.
Living birds are among the world’s most diverse and widely distributed tetrapods; they inhabit a myriad of different environments, and exhibit enormous disparity in their forms and lifestyles . However, this striking diversity is not distributed evenly across the globe. Today, Earth’s lower latitudes support the majority of avian species-level and higher-order diversity, and many deeply diverging clades are restricted to present-day vestiges of Gondwana, including Africa, Australasia, and South America [2,3,4]. This common distributional pattern has prompted the proposal of a Gondwanan origin for living birds, a hypothesis that has been used to corroborate arguments for an ancient Mesozoic diversification of the avian crown group [5,6,7]. Such arguments, which have historically favored a vicariant Gondwanan origin for crown birds (Neornithes), have often ignored data gained from the Palaeogene fossil bird record, which has improved substantially in recent years thanks to new discoveries and diagnoses based on rigorous phylogenetic analyses . In fact, phylogenetic hypotheses for many Northern Hemisphere Palaeogene bird fossils may cast doubt on the hypothesis of a Mesozoic Gondwanan origin of Neornithes, as many crown-clades with restricted extant distributions appear to have stem-group relatives in very different parts of the world [8,9,10]. For example, palaeontological analyses have suggested that taxa as diverse as total group seriemas (Cariamidae) [11,12,13,14,15,16], mousebirds (Coliidae) [17,18,19,20,21,22], and courols (Leptosomidae) [11, 23, 24] – all of which are currently restricted to formerly Gondwanan landmasses – have early stem-group representatives in the Palaeogene of the Northern Hemisphere. Furthermore, stem group representatives of some extant clades currently restricted to the New World (e.g., hummingbirds, Trochilidae) were formerly distributed in the Old World [25,26,27,28,29,30,31], while the opposite is true for certain extant taxa endemic to the Old World (e.g., the roller + ground roller clade, Coracioidea [32,33,34,35]). A recent attempt to address the biogeographic origin of Neornithes, incorporating the Cenozoic avian fossil record, suggests a Mesozoic origin in West Gondwana (comprising what is now South America, West Antarctica, and portions of East Antarctica), followed by subsequent expansion into North America via an early Palaeogene land mass linking South America and North America .
Earlier work on neornithine historical biogeography , which did not incorporate the Cenozoic avian fossil record, drew three major conclusions: 1) that Neornithes originated in the Southern Hemisphere; 2) that the distributions of major groups of crown birds were influenced by the breakup of Gondwana (which took place almost entirely in the Mesozoic); and 3), as a corollary of 2), that crown birds were not substantially affected by the end-Cretaceous mass extinction event. However, several recent divergence time analyses advocate a largely Cenozoic adaptive radiation of higher-level clades within the avian crown group [36, 37]. Furthermore, the Late Cretaceous fossil record suggests a devastating impact of the K-Pg mass extinction event on avian diversity ~ 66 million years ago (making it unlikely that many lineages of crown birds survived this extinction event) . As a result, the origins of today’s pervasive ‘trans-Antarctic’ neornithine biogeographic distributions are in need of additional study . Notably, analytical reconstructions of avian historical biogeography have only recently begun to incorporate the early fossil record of the avian crown group [4, 39, 40], placing a premium on robustly supported phylogenetic hypotheses for early crown bird fossils.
Despite the importance of addressing the phylogenetic position of Palaeogene fossil birds, numerous significant specimens are in need of reevaluation. One notable example is Foro panarium Olson 1992, from the Early Eocene (Wasatchian) Fossil Butte Member of the Green River Formation, Wyoming. The holotype and only known specimen of F. panarium (USNM 336261) is represented by a well-preserved, nearly complete, semi-articulated skeleton . In the initial description of this specimen, Olson  noted several osteological features shared with extant Musophagidae (including similarities in the scleral ossicles, ectethmoid, and the presence and shape of the pectineal process of the pelvis), and cautiously referred F. panarium to Order Cuculiformes Wagler 1830. Despite noting that this grouping (comprising Opisthocomidae, Musophagidae, and Cuculidae) probably did not represent a monophyletic group, this taxonomic decision was made on the basis of greater perceived overall similarity between F. panarium and the Hoatzin (Opisthocomus hoazin), turacos (Musophagidae), and cuckoos (Cuculidae) than to any other clade of extant birds. Despite the thorough initial description of F. panarium , and the near-complete nature of the holotype, the position of F. panarium has not been reassessed in an explicitly phylogenetic context.
Here we show that new phylogenetic analyses support Foro panarium as a stem turaco (Pan-Musophagidae). We thereby extend the list of Palaeogene fossil birds exhibiting geographic distributions that are extremely dissimilar to those of their closest living relatives, as all extant turacos (comprising a clade of ~ 24 species) are endemic to sub-Saharan Africa . Although the higher order phylogenetic placement of crown Musophagidae has historically been labile (e.g., contrast [36, 37, 43,44,45,46,47,48,49,50,51,52]), several recent studies uphold a comparatively recent common ancestor of crown turacos and crown cuckoos (Cuculidae), possibly as part of a broader clade including crown bustards (Otididae) [36, 37, 46, 53].
We tested the phylogenetic position of Foro panarium (Fig. 1) by performing a suite of phylogenetic analyses under different topological constraints and optimality criteria. The holotype of F. panarium (which has never previously been included in a phylogenetic analysis) was coded from direct observation. Parsimony and Bayesian phylogenetic analyses were performed in PAUP* v.4.0b10  and MrBayes v3.2.2 , respectively, with Palaeognathae specified as the outgroup. The character/taxon matrix (Additional file 1) consisted of 46 taxa and 153 morphological characters (of which 65 could be coded for F. panarium), and is based on a revised version of the dataset published by Mayr et al.  (itself a modification of Mayr and Clarke ). Two characters were newly added to the Mayr et al. dataset: 152. bill short and stout with broad processus maxillaris of the os nasale: no (0), yes (1); and 153. furcula unfused at midline: no (0), yes (1).
Given pervasive incongruities between the morphological phylogenetic topology inferred by Mayr and Clarke  and recent phylogenomic analyses of neornithine interrelationships (e.g., [36, 37, 45, 46, 57]), we analyzed the morphological matrix under a series of hierarchical topological constraints (‘scaffolds’ sensu Lee ) informed by the Hackett et al.  and Prum et al.  topologies. The same system of constraints was adhered to under parsimony and Bayesian optimality criteria. First, phylogenetic analyses were performed on the unconstrained morphological dataset (Fig. 2a). Subsequent analyses applied variations of the 50% ML majority-rule topology from Hackett et al.  as a topological scaffold (Fig. 2b, c, and d). These constrained analyses fixed the phylogenetic interrelationships of all taxa in the Hackett et al.  majority-rule topology, except those for which Foro was considered a potential fossil total-group representative by Olson  (i.e., Opisthocomidae, Cuculidae, and Musophagidae). The first scaffold analysis (Backbone 1) did not fix the phylogenetic position of these three taxa, and subsequent analyses (Backbone 2–4) sequentially fixed the position of Opisthocomidae, Opisthocomidae + Cuculidae, and finally Opisthocomidae + Cuculidae + Musophagidae. Additionally, constrained analyses were performed following the recent phylogenomic topology of Prum et al. , applying the same methodology. Constraint trees are provided as separate Additional file 2.
Several osteological similarities between F. panarium and the Hoatzin (Opisthocomus hoazin) were noted in the original description of F. panarium . Given the historical difficulty of identifying the Hoatzin’s extant sister taxon (the lineage leading to O. hoazin may represent the single longest branch in the neornithine tree of life ), and its strikingly unusual skeleton, an effort was made to shorten the phylogenetic branch leading to O. hoazin by including two stem opisthocomids in the phylogenetic analysis: Hoazinavis lacustris Mayr et al. 2011, from the late Oligocene/early Miocene of Brazil; and Namibiavis senutae Mayr 2014, from the late early Miocene of Namibia [56, 59]. The position of these fossils was left unconstrained in all phylogenetic analyses. Results of phylogenetic analyses are presented in Fig. 2.
For the parsimony analyses, heuristic searches were performed under tree-bisection-reconnection (TBR) branch swapping with 1000 random stepwise sequence addition replicates. Minimum branch lengths were set to collapse. Node support was calculated using bootstrap frequencies with 1000 bootstrap replicates and 10 random sequence addition replicates. Characters 55, 71, and 91 were treated as ordered (following the original Mayr and Clarke dataset ).
Bayesian phylogenetic analyses were run using MrBayes  on the CIPRES Science Gateway , using the Mk model of morphological evolution  with gamma-distributed rate variation and variable coding. All analyses were performed with two concurrent runs, sampling frequency of 1000, and four Metropolis-coupled chains (T = 0.1) for 10 million generations. Characters 55, 71, and 91 were again treated as ordered. Analyses were checked for convergence using standard MrBayes diagnostics (e.g., PRSF< 0.01, mixing between chains > 20%) and Tracer (v1.5)  (e.g., ESS > 200 for all estimated parameters). For all summary statistics, a relative burn-in of 25% was applied. Additional file 3: Table S1 provides the full character/taxon matrix. Although variations of the core anatomical dataset used to score F. panarium in this study have been subjected to numerous phylogenetic analyses under parsimony (e.g., [50, 63,64,65]), to our knowledge, the present study represents the first time this matrix has been analyzed within a Bayesian phylogenetic framework.
Historical biogeographic analyses
To obtain a time-scaled phylogeny, we used a majority-rule consensus tree identified by Mesquite  based on a sample of 1000 trees from the posterior distribution of Jetz et al.  (http://www.birdtree.org). This tree also applied the Hackett et al.  topology as a higher-order phylogenetic scaffold, allowing for direct comparisons with the constrained analyses described above. Following Cracraft , we restricted our sampling to crown neornithine family-level clades whose present-day biogeographic distributions are either exclusively Gondwanan or exclusively Laurasian.
The biogeographic history of birds in the Northern Hemisphere has been complex throughout the Cenozoic. In particular, the European fossil record indicates that several bird groups restricted to lower-latitudes today, including total group Struthionidae, Bucorvidae, and Opisthocomidae, may have migrated into higher latitude environments in the Northern Hemisphere during a period of warmer global temperatures around the Miocene climatic optimum ~ 15 MYA, after having arisen in lower-latitude settings [15, 68, 69]. In order to avoid conflating Palaeogene biogeographic patterns with subsequent overprinting of Northern Hemisphere avifaunas by dispersal from the tropics during the Miocene , we elected to restrict our fossil sampling to Palaeogene localities.
To evaluate the complexity of avian biogeographic history in light of the Palaeogene bird fossil record, we performed two separate analytical reconstructions of the biogeographic origin of the avian crown clade using multiple techniques. These analyses were designed to test the capacity of fossils to illustrate profound fluctuations in avian biogeographic distributions that would otherwise be unknowable. The first was a quantitative approximation of the scenario put forth by Cracraft , incorporating every extant non-passerine clade traditionally ranked at the family level that exhibits an exclusively ‘Gondwanan’ or ‘Laurasian’ biogeographic distribution, and no fossils (Fig. 3a). Although certain clades (e.g., Nyctibiidae) extend into southernmost North America, or even slightly beyond (e.g., Cracidae), the northward extensions of their ranges are almost certainly attributable to contiguous dispersal following uplift of the Panamanian isthmus. Accordingly, these clades were scored as Gondwanan, as they were by Cracraft . Other clades, such as the Podargidae (frogmouth) crown group, are found exclusively at lower latitudes, but extend from Australia into Southeast Asia (not a vestige of Gondwana). These clades were excluded from the Gondwanan-only coding scheme. It should be noted, however, that probable stem group podargids are known from the Palaeogene of both Europe and North America; as such, the Palaeogene distribution of the Podargidae total-group greatly exceeded the apparently relictual range occupied by the crown clade today (e.g., [63, 70]). The same is true of several other extant pantropical groups not included in this analysis, including trogons (Trogoniformes) and parrots (Psittaciformes), both of which have extensive Northern Hemisphere total group fossil records .
The second analysis (Fig. 3b) builds upon the first by incorporating all well-supported extinct Palaeogene sisters to the extant clades in Fig. 3a. For a complete list of the fossil taxa included, see Table 1. We illustrate the ability of the avian fossil record to introduce important new biogeographic information by performing ancestral area reconstructions in a Bayesian framework using the R (v. 2.15.1)  package phytools (function make.simmap) , under the equal rates (ER) model with an estimated stationary prior distribution on the root, empirical rate matrix, and 500 simulation replicates. Extant and extinct taxa were coded as either ‘Gondwanan’ or ‘Laurasian’, and the results of these reconstructions are presented at the nodes of Fig. 3a and b.
In addition to the ‘simplistic’ biogeographic reconstruction described above, we performed biogeographical ancestral state reconstructions using RASP (v.3.2 build 20,160,719)  under two models: S-DIVA (Statistical-Dispersal Vicariance Analysis; ) and BayArea . Both analyses were conducted using the phylogeny built from the full dataset containing both extant and fossil taxa, with Todidae coded as Laurasian (i.e., equivalent to the tree shown in Fig. 3b). As we use only a single phylogenetic tree and do not exclude any possible combined geographic ranges, our S-DIVA analysis is equivalent to using the standard DIVA model . The maximum number of areas at each node was set to two, extinctions were allowed, and ancestral node ages were entered manually.
For the BayArea analysis, we represented the relative positions of Laurasia and Gondwana by using the centroid latitude/longitude coordinates of the two landmasses, as required by the BayArea algorithm. To do this, we used the program GPlates (v.1.5.0 build 16,091)  to reconstruct the tectonic plate positions of the two supercontinents at 118 MYA (the approximate root age of the Jetz et al. 2012 tree). The rotation model, coastline data, model gridmarks, global isochron data, continent-ocean boundary data, and spreading ridge data were sourced from Seton et al. . The GPlates polygon digitization tool was used to trace approximate outlines of Laurasia and Gondwana by hand, resulting in 56 and 35 modern longitude-latitude coordinates, respectively. These coordinates were then exported as ESRI shapefiles, and the centroid coordinates of each polygonal area were then calculated using the gCentroid function in the R  package rgdal . Because the polygon enclosing the Laurasian landmass overlaps with the North Pole, additional processing was needed to obtain its centroid due to the limitations of the gCentroid algorithm. Namely, the original longitude-latitude coordinates, which were specified under the WGS84 geographic coordinate system, were first converted to a planar coordinate reference system (CRS) via EPSG projection 3575 (North Pole LAEA [Lambert azimuthal equal-area] Europe) using the rgdal function spTransform. The centroid was then calculated using the converted planar coordinates using the gCentroid function. Finally, the planar centroid coordinates were reconverted into longitude-latitude coordinates under the WGS84 CRS. The final coordinates were 79.65 lat and 27.97 long for Laurasia and − 43.32 lat and 9.41 long for Gondwana. The BayArea analysis was run using these coordinates for 50 million generations, sampling every 1000 generations, using default parameter and prior settings. A 25% burn-in of 12.5 million generations was implemented after the run before summarization of results.
Recently, the use of fossil tip-dating  to investigate the influence of fossils on ancestral state reconstructions has increased, allowing for more accurate reconstructions (e.g., ). However, due to the lack of a sufficiently taxonomically-inclusive anatomical phylogenetic dataset for the avian crown, a ‘total-evidence’ phylogenetic analysis that explicitly incorporates nucleotide sequences and anatomical data across crown bird diversity, along with fossil data from Palaeogene representatives, was beyond the scope of this investigation (although such an analysis has recently been applied to the avian subclade Galloanserae ). Instead, we incorporated Palaeogene fossils into the Jetz et al.  time-calibrated consensus tree by grafting them into their phylogenetic positions as inferred by independent phylogenetic analyses (Table 1). The ages of the fossils grafted into the consensus tree followed the age of the earliest well-supported representative of the fossil species (Table 1). These age estimates follow best practices for justifying minimum age constraints ; thus, we apply the youngest possible age for each fossil, inclusive of error. The node subtending each grafted fossil and its extant sister clade was placed 20 million years prior to the age of the fossil, except for the nodes within total-clade Spheniscidae and Anseranatidae, because this would conflict with the age for the node subtending crown Spheniscidae and Anseranatidae and their extant sisters from the Jetz et al.  analysis. For these clades, a five million-year offset was applied instead. See Additional file 2 for details regarding alternative biogeographic parameterizations for taxa with biogeographic distributions that were challenging to characterize.
Digital calipers sensitive to 0.01 mm were used to measure the total length of the femur, tibiotarsus, and tarsometatarsus for adult representatives of Otidimorphae (Musophagidae, Cuculidae, and Otididae). Measurements are presented in Additional file 4: Table S2 and plotted in Fig. 4. In total, six species of Musophagidae were examined, representing every major musophagid subclade. Twenty-eight species of Cuculidae were examined from across the extant diversity of cuckoos, and four species of Otididae were measured. Measurements from extant taxa were compared to measurements from Foro panarium presented in Olson  and to the possible stem cuculid Eocuculus cherpinae [84, 85]. A hindlimb index was computed as the ratio between the sum of tibiotarsus and tarsometatarsus length, and femur length. Results of morphometric comparisons are displayed in Fig. 4.
Despite exhibiting an enigmatic combination of osteological features not seen in any extant birds , Foro panarium was supported in our analyses as the extinct sister to crown Musophagidae under all optimality criteria and topological constraints (Fig. 2). In our Bayesian analyses, posterior probabilities for the exclusive Foro-Musophagidae clade ranged from 0.80 (Hackett Backbone with cuckoos, hoatzin, and turacos unconstrained) to 0.96 (Prum Constraint). Under parsimony, the Foro-Musophagidae node was supported with bootstrap scores ranging from 0.57 (Hackett Backbone with cuckoos, hoatzin, and turacos unconstrained) to 0.95 (Hackett Backbone with all extant taxa constrained). In the unconstrained parsimony analyses, the most parsimonious tree (which resolved Foro as sister to Musophagidae) was found to be 215 steps; 217 were necessary to resolve Foro as sister to Cuculidae, and 220 were necessary to resolve Foro as sister to Opisthocomidae.
We found substantial incongruence between our unconstrained morphological analyses and recently proposed genetic sequence-based hypotheses for neornithine interrelationships (e.g., [36, 37, 45]) under both parsimony and Bayesian optimality criteria, underscoring the extant avian radiation as an exemplar of discordance among alternative phylogenetic datasets [86, 87]. Although considerable progress has been made toward resolving the deepest nodes within the neoavian radiation, complete consilience has been elusive even among recent next-generation phylogenomic datasets [36, 37], perhaps as a result of widespread incomplete lineage sorting during the early stages of the neoavian radiation , or long-branch attraction artifacts in datasets with sparse taxonomic samples . From an anatomical standpoint, achieving congruence between phenotypic and genotypic analyses of neornithine interrelationships demands a renewed focus on expanding and refining existing neornithine character/taxon datasets, and additional effort to integrate phenotypic and genotypic data in phylogenetic studies [80,81,82, 89, 90].
Foro panarium exhibits a furcula unfused at its midline (Character 153:1), which is recovered as an unambiguously optimized synapomorphy of a Foro + Musophagidae clade. Although the majority of the left furcular ramus is separated from the right ramus and rotated perpendicular to its original position (Fig. 1), the medial portion of the left ramus is in its original orientation. Together with the complete and unbroken right furcular ramus, the medial portion of the left ramus bounds an unfused midline symphysis. Together, the medial and lateral components of the left furcular ramus match the length and shape of the intact right ramus. Monophyly of crown Musophagidae to the exclusion of F. panarium is supported by several character states, including: a largely ossified septum nasale (character 8:1); a pygostyle perforated by a foramen at its caudoventral end (character 61:0); an extremitas omalis of the furcula with strongly developed, laterally protruding facies articularis acrocoracoidea (character 62:1); and crista deltopectoralis of humerus extending distal to crista bicipitalis and proximodorsal portion of bone with sigmoidally curved margin (character 151:1). F. panarium exhibits the same states for these four characters that bustards (Otididae) do, which presumably reflects the plesiomorphic condition for Otidimorphae.
Although character state optimization varied across our alternative phylogenetic analyses, character states consistently optimizing as unambiguous local synapomorphies of a Foro panarium + Musophagidae clade included 93:0 (large tubercula praeacetabularia of the pelvis); 86:1 (carpometacarpus, proximal end of os metacarpale minus dorsoventally wide and strongly deflected ventrally); 87:1 (os carpi ulnare with crus longum greatly abbreviated); 152:1 (bill short and stout with broad processus maxillaris of the os nasale); and 153:1 (furcula unfused at midline). Character state 87:1 was posited as a potential synapomorphy of a Musophagidae + Cuculidae clade by Hughes —accordingly, this character state did not optimize as an unambiguous F. panarium + Musophagidae synapomorphy in analyses imposing backbone constraints where Musophagidae and Cuculidae are closely related (e.g. ), and may in fact represent a synapomorphy for Otidimorphae more generally.
Despite the inclusion of the branch-shortening stem hoatzins Hoazinavis lacustris and Namibiavis senutae [56, 59], an exclusive Foro + Opisthocomidae clade was never recovered. In our Bayesian analyses, N. senutae and H. lacustris never formed an exclusive clade with crown Opisthocomidae under majority rules consensus (instead emerging in a polytomy with Opisthocomus), likely due to missing data arising from the incomplete nature of those fossil specimens. Osteological similarities between F. panarium and Opisthocomus hoazin, which mostly relate to the overall shape of the skull and mandible (including the ‘short, hook-like retroarticular process’ ), are therefore best interpreted as the result of convergence.
Ancestral area reconstructions
Our historical biogeographic reconstructions revealed marked differences between analyses including and excluding Palaeogene fossil birds. In the extant-only biogeographic reconstruction (Fig. 3a), the hypothesis that pervasive ‘trans-Antarctic’ biogeographic patterns across crown birds indicate a Southern Hemisphere origin for Neornithes was strongly supported (pro ). Indeed, the maximum likelihood ancestral area reconstruction for the neornithine root is inferred to have been Gondwanan with 100% probability. However, inclusion of Palaeogene bird fossils obfuscates this inference (Fig. 3b). All nine variations of the extant+fossil dataset yield ambiguous reconstructions, with the probability of an exclusively Gondwanan root node varying between 45.2 and 52.2% (Additional file 5: Figure S1, contra ). This analysis illustrates the unique potential of the fossil record to illustrate where representatives of major avian subclades were formerly distributed, and underscores the profound historical fluctuations in avian biogeography revealed by the fossil record.
The results of our S-DIVA analysis are shown in Fig. 5. This method reconstructs the root as being 100% Gondwanan despite the inclusion of Laurasian fossils, similar to our naïve likelihood extant-taxa-only results. Under the simple DIVA model, branch lengths are ignored; we believe that this limitation of the method explains the simplistic 100% Gondwanan reconstruction, due to the ubiquity of Gondwanan distributions among the extant taxa in our analysis. In contrast, the BayArea analysis reconstructs a far more ambiguous result, with the root node reconstructed as 81.69% Laurasian/Gondwanan, 11.33% Gondwanan, and 6.98% Laurasian (Fig. 6).
Hindlimb ratios (HR) for extant musophagids (which are arboreal) are uniformly low, with the highest measured musophagid ratio exhibited by the White-bellied Go-away Bird, Corythaixoides leucogaster (HR = 2.32) (Fig. 4; Additional file 4: Table S2). Musophagid hindlimb ratios are exceeded by most Cuculidae, although the lowest overall hindlimb ratio was exhibited by the generally arboreal Common Cuckoo, Cuculus canorus (HR = 2.13). The exclusively terrestrial Otididae exhibit uniformly high hindlimb ratios, with exceptionally high values seen in Ardeotis (A. arabs, HR = 3.63; A. kori, HR = 3.53). Foro panarium exhibits a hindlimb ratio (HR = 2.77) intermediate between that of the Little Bustard (Tetrax tetrax, HR = 2.66) and the Great Bustard (Otis tarda, = 2.85), which also closely corresponds to ratios exhibited by some (terrestrial) neomorphine ground cuckoos, including Dromococcyx phasianellus (HR = 2.78) and Geococcyx californianus (HR = 2.70). The hindlimb ratio of F. panarium greatly exceeds those of any extant musophagids, although the femur length (54.1 mm) falls within the range of variation exhibited by crown turacos, suggesting a live body mass between 636 g–714 g . This estimated body mass falls between the mean body size of the large extant turacos Crinifer zonurus and Corythaeola cristata, and within the known range of variation for small extant bustards such as Eupodotis ruficrista .
Discussion and conclusions
The evolutionary history of total group turacos (Pan-Musophagidae)
Crown turacos (Musophagidae, Musophagiformes) comprise ~ 24 traditionally recognized species, many of which exhibit vivid plumage coloration conferred by unique porphyrin pigments . Until recently, the phylogenetic position of crown turacos had been among the most recalcitrant phylogenetic problems in neornithine systematics, with turacos having been alternatively allied with a number of distantly related clades (e.g., cuckoos [93, 94], the hoatzin , mousebirds , and waterbirds ). Recently, however, the phylogenetic position of turacos within Neoaves has begun to come into focus, with congruent results emerging from next-generation phylogenomic datasets supporting a close relationship between turacos, bustards (Otididae), and cuckoos (Cuculidae) [36, 37, 46] (although McCormack et al.  did not sample Cuculidae). Both McCormack et al.  and Jarvis et al.  inferred a sister-group relationship between turacos and bustards, whereas Prum et al.  inferred a sister group relationship between turacos and a bustard + cuckoo clade. The most exclusive clade uniting Musophagidae, Cuculidae, and Otididae, with variable ingroup interrelationships, has been named Otidimorphae [36, 37].
The known fossil record of early turacos is sparse. Although molecular divergence dating analyses conflict markedly on the estimated antiquity of the turaco total group (e.g., contrast Mesozoic estimates from Jetz et al.  and Palaeogene estimates from Prum et al. ), no convincing total group turaco fossils have previously been identified in sediments older than the early Oligocene [8, 96, 97]. Other putative Oligocene records of turacos  have been reinterpreted as belonging to other clades . Our inference that Foro panarium represents an early Eocene stem turaco provides a well-supported calibration point for divergence time analyses as both the earliest representative of the turaco total group and the earliest fossil representative of Otidimorphae, thereby filling a major gap in the Cenozoic neornithine fossil record.
Assuming a post-Cretaceous radiation of crown Neoaves [36,37,38, 44, 99, 100], an early Eocene (~ 52.1MYA) stem turaco might be expected to exhibit a mix of crown turaco apomorphies as well as plesiomorphies that have been overwritten in extant Musophagidae. Indeed, Foro panarium appears intermediate, possessing some crown turaco apomorphies (e.g., unfused furcula), while still exhibiting numerous osteological plesiomorphies (e.g., unossified septum nasale; lack of prominent quill knobs on ulna).
Among the most surprising features of the skeleton of Foro panarium is the presence of long, gracile legs . Extant turacos are highly arboreal, and exhibit fairly short legs (Fig. 4; Additional file 4: Table S2). The emerging hypothesis that turacos share a recent common ancestor with bustards (large-bodied, terrestrial birds with long legs) raises questions regarding the evolution of divergent hindlimb proportions in extant otidimorphs [36, 37]. Among crown group Otidimorphae, the hindlimb proportions of F. panarium are much more similar to those of bustards and certain New World ground cuckoos than they are to extant turacos (Fig. 4). The elongated hindlimbs of F. panarium may therefore approximate the ancestral condition for Otidimorphae . The shortened hindlimbs of many cuckoos and turacos may have arisen either on their shared stem to the exclusion of bustards (as advocated by the Jarvis et al.  TENT topology), or independently along both the crownward portion of the turaco stem and the cuckoo stem (following the Prum et al.  topology). Unfortunately, the early evolutionary history of Cuculidae is poorly understood, with no well-supported stem group representatives . As a result, the total group cuckoo fossil record may have little to contribute to the discussion of the pattern and timing of otidimorph hindlimb reduction – and the evolution of arboreal habits – at present. However, the late Eocene-early Oligocene taxon Eocuculus cherpinae Chandler 1999 (known from the Florissant Beds of Colorado and the Lubéron area of France) may represent a stem cuculid, although this hypothesis awaits further testing on the basis of phylogenetic analysis and the discovery of additional fossils [84, 85, 102]. Based on measurements of the holotype , the hindlimb index of E. cherpinae is considerably lower than any extant otidimorph (Fig. 4; Additional file 4: Table S2); if this taxon does in fact represent a stem cuculid, it may suggest that cuckoos arose from taxa with extremely short hindlimbs. However, a clade comprising the relatively long-legged ground cuckoos (Neomorphinae), together with the new world Crotophaginae, represents the extant sister taxon to all other extant cuckoos . The phylogenetic position of Neomorphinae is consistent with a long-legged plesiomorphic condition for Otidimorphae, which would suggest that the short legs of E. cherpinae are autapomorphic (if E. cherpinae is indeed a stem cuckoo). Hindlimb proportions of F. panarium and other otidimorphs, in combination with exclusively terrestrial habits in bustards and widespread terrestrial ecologies within the extant cuckoo radiation, suggest that turacos may be descended from ground-dwelling antecedents, and that the arboreal habits of extant turacos arose relatively recently in their evolutionary history .
Claramunt and Cracraft  infer the probability of a South American vs. African origin for Pan-Musophagidae as equally parsimonious. If a South American origin for total-clade turacos is accurate, then the presence of an early stem turaco in the Eocene of North America is consistent with their interpretation of an early Cenozoic dispersal event from South America through North America (the ‘North American Gateway’ hypothesis). This scenario is congruent with prior inferences of North American origins for the (presently) Old World clades Coliidae and Leptosomidae, as well as the cosmopolitan Coraciiformes, corroborating the argument that North America played a pivotal role in the early evolutionary history of clades that subsequently radiated in the Old World .
Given the striking similarity of North American and European avifaunas in the early Eocene [32, 104, 105], it would be somewhat surprising if stem group turacos are not eventually recovered from Messel, or other early Cenozoic European localities. Europe was connected to North America across the Greenland-Scotland ridge in the early Eocene [106,107,108], which has been causally attributed to similarities in mammalian, squamate, and avian faunas from that time [109,110,111]. However, given that only one specimen referable to Foro panarium has yet been identified, it is possible that early stem turacos were comparatively rare in habitats surrounding lakes such as those of the Green River system. If so, they may have been similarly uncommon in the lakeside habitats surrounding Messel . Numerous avian taxa from the Palaeogene of the Northern Hemisphere, whose present-day distributions are restricted to the Afrotropical zoogeographic realm (e.g., Coliiformes), are known to have persisted in Europe until the Miocene . Therefore, as-yet undetected early stem group musophagids may have survived longer in Europe than they did in North America (if they were present in Europe in the first place), and may be expected to be recovered in stratigraphically younger Palaeogene European localities.
Implications for molecular divergence time analyses
The early Eocene Fossil Butte Member of the Green River Formation has been the subject of detailed dating analyses. Smith et al.  reported a radiometric date of 51.97 ± 0.16 Ma from a K-feldspar tuff just above the middle unit of the Fossil Butte Member . As a result, Foro panarium easily represents the oldest representative of Otidimorphae yet known. Although undescribed fossils purportedly belonging to Otididae have been recovered from the Oligocene of Kazakhstan  and from an unknown horizon within the Oligocene Quercy fissure fillings of France , no pre-Oligocene remains of Otididae have yet been reported. Similarly, although extant Cuculidae are distributed on every continent except Antarctica , potential fossil representatives of total-clade Cuculidae are unknown from sediments older than the late Eocene Florissant Fossil Beds of Colorado, the roughly contemporaneous Cypress Hills Formation of Saskatchewan, and the early Oligocene Pichovet locality in France [8, 84, 85, 102, 116].
Our inferred phylogenetic position for Foro panarium as sister to Musophagidae hints not only at the antiquity of the musophagid total group, but also indicates that additional phylogenetic divergences within Otidimorphae had taken place by the early Eocene. If Musophagidae represents the extant sister taxon to a Cuculidae + Otididae clade (sensu ), then at least the stem group of this latter clade would have been present in the early Eocene. If, instead, Cuculidae is sister to a Musophagidae + Otididae clade (sensu ), then stem representatives of all three of these extant clades must have been present in the early Eocene.
Foro panarium therefore represents a valuable calibration point for neornithine molecular divergence dating analyses. Otidimorphae are among the deepest diverging clades within living birds (diverging near the base of Neoaves; [36, 37]). Thus, the presence of a stem musophagid in the early Eocene of North America has important implications for understanding the tempo of the extant neornithine radiation—one of the most contentious topics in contemporary ornithology [38, 117,118,119,120,121,122]. Recent summaries of useful criteria for well-justified fossil calibrations (e.g., [83, 123,124,125]) identify six desirable features for a potential calibrating fossil; F. panarium satisfies all of these criteria, and therefore should help shed new light on the tempo and mode of the extant avian radiation.
Neornithine historical biogeography and the importance of fossils
The importance of incorporating fossils into macroevolutionary analyses has been well documented. Including fossils can meaningfully influence phylogenetic analyses (e.g., [81, 126,127,128,129,130]), ancestral state reconstructions (e.g., [81, 120, 131]), and models of trait evolution (e.g., [132, 133]). Similarly, synthesizing biogeographic data for extant clades and their fossil relatives has potential to provide a more holistic understanding of the dynamics of geographic range evolution than do studies focused exclusively on either the current distributions of clades or those of their fossil antecedents [134,135,136,137]. Indeed, the Palaeogene record of fossil birds has been discussed at length with respect to the historical biogeography of crown birds [8,9,10, 15, 32, 104, 105, 138, 139]. However, the results of these studies have rarely been incorporated into historical biogeographic analyses of living birds groups (however, see ). As a result, the Palaeogene fossil record of birds has often been dismissed for having been inadequately examined from a phylogenetic perspective (e.g., ). Although much of the early work on Palaeogene bird fossils was conducted outside of an explicitly phylogenetic paradigm (including the initial description of Foro panarium ), a host of Palaeogene birds have been subjected to phylogenetic analyses over the last 20 years, providing strong evidence for the presence of the stem lineages of many higher clades (‘Orders’ under the traditional Linnean hierarchy) in the Palaeogene, often well outside the geographical limits of their respective crown clades .
The present study represents one of the first efforts to explicitly incorporate the Palaeogene neornithine fossil record into analytical reconstructions of higher-order avian historical biogeography. The results underscore the importance of the fossil record for informing our understanding of both the pattern by which crown neornithine biogeography has evolved and the timing of the extant neornithine radiation. Cracraft  identified numerous examples of apparent trans-Antarctic distribution patterns across the neornithine tree of life – that is, avian higher clades restricted to landmasses that were formerly part of the Gondwanan supercontinent (i.e., Africa, Madagascar, Antarctica, Asia, Australia/New Zealand, and South America). Indeed, an ancestral area reconstruction based only on the current geographic distributions of these clades recovers a maximally-supported result in which both the most recent common ancestor of all living birds, and most internal nodes, favor an origin somewhere within Gondwana or its derivatives (Fig. 3a). This result has been interpreted to suggest that vicariance related to the breakup of the Gondwanan supercontinent, which was largely complete by the end of the Mesozoic era, was responsible for establishing this striking biogeographic pattern . Because this biogeographic scenario necessarily invokes numerous neornithine divergences deep in the Cretaceous, this scenario has also been taken as support for a major Mesozoic radiation of crown birds . Despite apparent corroboration of this temporal scenario by numerous molecular divergence time analyses [67, 140,141,142], the hypothesis of a major radiation of the neornithine crown in the Mesozoic has elicited mounting criticism in recent years. Recent molecular divergence dating analyses have begun to corroborate a rapid post-Cretaceous radiation of the major groups of crown birds [36, 37, 120,121,122]. If the results of these more recent analyses reflect the approximate age of the avian crown, any model invoking Mesozoic Gondwanan vicariance as a driver of crown avian biogeographic patterns can be soundly rejected.
Large-scale analytical reconstructions of the biogeographic origin of crown birds favor a West Gondwanan origin of Neornithes (comprised of the landmasses now belonging to South America and portions of Antarctica), a hypothesis that is robust to the inclusion of numerous Palaeogene fossils from the Northern Hemisphere . Although including Palaeogene fossils in our biogeographic analysis rendered the ‘Gondwanan vs. Laurasian’ reconstruction ambiguous (Fig. 3), that result is not inconsistent with the West Gondwanan origin of crown birds advocated by Claramunt and Cracraft . Regardless, it is clear that the Palaeogene fossil bird record from the Northern Hemisphere will be instrumental for unraveling the biogeographic history of Neornithes. Indeed, the presence of Foro panarium in the early Palaeogene of North America, and its identification as a stem turaco, is consistent with the ‘North American Gateway’ hypothesis – the idea that, after initially diverging in West Gondwana, North America served as a stepping stone for many avian clades on their way to colonizing the Old World via high latitude land bridges . The continued identification of Laurasian fossils belonging to the early stem lineages of clades presently restricted to the Southern Hemisphere, such as F. panarium, may fundamentally alter our understanding of the early biogeographic origins of living bird clades. Continued efforts to assess the phylogenetic affinities of bird fossils from the Northern and Southern Hemispheres will be a crucial step towards fully understanding the biogeographic and temporal origins of modern bird diversity.
The ambiguity of the BayArea result mirrors that obtained from our naïve likelihood method with fossils included, and underscores the important influence of branch length information in historical biogeographic reconstructions (and comparative analyses more broadly). Although most of the extant taxa in our analyses are found on modern landmasses that comprised the Gondwanan paleo-supercontinent, most stem fossil representatives of these groups are found instead on Laurasia, and are subtended by relatively short branches – particularly in comparison to the length of subtending branches leading to extant taxa, along which there is comparatively ample opportunity for biogeographic events such as dispersal and vicariance to occur. In conjunction with our naïve likelihood analyses, the results of our explicitly biogeographic reconstructions of the distribution of early birds provide an important illustration of the influence of fossils in such analyses. When informative branch length data are incorporated, the fossil record provides indispensable data on evolutionary and biogeographic history, leading to reconstructions that may be unexpected when one considers extant data alone. This is true even when analyses containing only extant data result in seemingly unambiguous and unequivocal reconstructions (e.g., a 100% Gondwanan reconstruction for the root node). Indeed, our results serve as a cautionary warning against overreliance on extant data and high statistical support values when studying evolutionary processes that are fundamentally historical in nature.
Million years ago
Gill FB. Ornithology, Second Edition. New York: W. H. Freeman and Co.; 1995.
Cracraft J. Continental drift, paleoclimatology, and the evolution and biogeography of birds. J Zool. 1973;169(4):455–543.
Cracraft J. Avian evolution, Gondwana biogeography and the cretaceous-tertiary mass extinction event. Proc R Soc B Biol Sci. 2001;268(1466):459–69.
Claramunt S, Cracraft J. A new time tree reveals earth history’s imprint on the evolution of modern birds. Sci Adv. 2015;1(11):e1501005.
Hedges SB, Parker PH, Sibley CG, Kumar S. Continental breakup and the ordinal diversification of birds and mammals. Nature. 1996;381(6579):226–9.
Rambaut A, Bromham L. Estimating divergence dates from molecular sequences. Mol Biol Evol. 1998;15(4):442–8.
Cooper A, Penny D. Mass survival of birds across the cretaceous-tertiary boundary: molecular evidence. Science. 1997;275:1109–13.
Mayr G. Paleogene Fossil Birds. Berlin: Springer; 2009.
Mourer-Chauviré C. Les oiseaux fossiles des phosphorites du Quercy (Éocène supérieur a Oligocène supérieur): Implications paléobiogéographiques. Geobios. 1982;15(Supplement 1):413–26.
Olson SL. Aspects of global avifaunal dynamics during the Cenozoic. In: Ouellet H, editor. Acta XIX congressus internationalis ornithologici. Ottawa: University of Ottawa Press; 1989. p. 2023–9.
Mayr G. Avian Remains from the Middle Eocene of the Geiseltal (Sachsen-Anhalt, Germany). In: Zhou Z, Zhang F, editors. Proceedings of the 5th symposium of the Society of Avian Paleontology and Evolution Beijing, 1–4 June 2000. Beijing: Science Press; 2002. p. 77–96.
Mayr G. New or previously unrecorded avian taxa from the middle Eocene of Messel (Hessen, Germany). Mitt Mus Naturkd Berl Geowiss Reihe. 2000;3:207–19.
Peters DS. Idiornis tuberculata n. spec., ein weiterer ungewöhnlicher Vogel aus der Grube Messel (Aves: Gruiformes: Cariamidae: Idiornithinae). In: Peters DS, editor. Acta Palaeornithologica, vol. 181. Frankfurt: Cour Forsch Inst Senckenberg; 1995. p. 107–19.
Mourer-Chauviré C. Les Gruiformes (Aves) des Phosphorites du Quercy (France). 1. Sous-ordre Cariamae (Cariamidae et Phorusrhacidae). Palaeovertebrata. 1983;13(4):83–143.
Mourer-Chauviré C. Les relations entre les avifaunes du Tertiaire inférieur d'Europe et d'Amérique du Sud. Bull Soc Geol Fr. 1999;170(1):85–90.
Mourer-Chauviré C. The avifauna of the Eocene and Oligocene phosphorites du Quercy (France): an updated list. Strata Ser. 2006;1(13):135–49.
Mayr G, Peters D. The mousebirds (Aves: Coliiformes) from the middle Eocene of Grube Messel (Hessen, Germany). Senckenb Lethaea. 1998;78(1–2):179–97.
Mayr G. New specimens of the middle Eocene fossil mousebird Selmes absurdipes Peters 1999. Ibis. 2001;143(4):427–34.
Houde P, Olson SL. A radiation of coly-like birds from the early Eocene of North America (Aves: Sandcoleiformes new order). In: Campbell KE, editor. Papers in avian paleontology honoring Pierce Brodkorb, vol. 36. Los Angeles City: Natural History Museum Science Series; 1992. p. 137–60.
Ksepka DT, Clarke JA. Affinities of Palaeospiza bella and the phylogeny and biogeography of mousebirds (Coliiformes). Auk. 2009;126(2):245–59.
Ksepka DT, Clarke JA. New fossil mousebird (Aves: Coliiformes) with feather preservation provides insight into the ecological diversity of an Eocene north American avifauna. Zool J Linnean Soc. 2010;160:684–706.
Zelenkov NV, Dyke G. The fossil record and evolution of mousebirds (Aves: Coliiformes). Palaeontology. 2008;51:1403–18.
Mayr G. The Madagascan “cuckoo-roller” (Aves: Leptosomidae) is not a roller–notes on the phylogenetic affinities and evolutionary history of a “living fossil”. Acta Ornithologica. 2008;43(2):226–30.
Mayr G. A new species of Plesiocathartes (Aves:? Leptosomidae) from the middle Eocene of Messel, Germany. PaleoBios. 2002;22(1):10–20.
Karhu A. A new family of swift-like birds from the Paleogene of Europe. Paleontol J. 1988;3:78–88.
Karhu A: Morphological divergence within the order Apodiformes as revealed by the structure of the humerus. In: Papers in avian paleontology honoring Pierce Brodkorb. Edited by Campbell KE, Los Angeles: Nat Hist Mus Los Angeles Cty Sci Ser; 1992,vol.36 379–384.
Karhu A: A new genus and species of the family Jungornithidae (Apodiformes) from the late Eocene of the Northern Caucasus, with comments on the ancestry of hummingbirds. In: Avian paleontology at the close of the 20th century Proceedings of the 4th international meeting of the Society of Avian Paleontology and Evolution. Edited by Olson SL, Washington, D.C.: Smithson Contrib Paleobiol; 1999:vol. 89 207–216.
Mayr G. Phylogeny of early tertiary swifts and hummingbirds (Aves: Apodiformes). Auk. 2003;120(1):145.
Mayr G. Old world fossil record of modern-type hummingbirds. Science. 2004;304(5672):861.
Bochenski Z, Bochenski Z. An old world hummingbird from the Oligocene: a new fossil from polish Carpathians. J Ornithol. 2008;149(2):211–6.
Louchart A, Tourment N, Carrier J, Roux T, Mourer-Chauviré C. Hummingbird with modern feathering: an exceptionally well-preserved Oligocene fossil from southern France. Naturwissenschaften. 2008;95(2):171–5.
Houde P, Olson SL. Small arboreal nonpasserine birds from the early Tertiary of western North America. In: Ouellet H, editor. Acta XIX Congressus Internationalis Ornithologici. Ottawa, University of Ottawa Press; 1989. p. 2030–6.
Mayr G, Mourer-Chauviré C, Weidig I. Osteology and systematic position of the Eocene Primobucconidae (Aves, Coraciiformes sensu stricto), with first records from Europe. J Syst Palaeontol. 2004;2:1–12.
Clarke JA, Ksepka DT, Smith NA, Norell MA. Combined phylogenetic analysis of a new north American fossil species confirms widespread Eocene distribution for stem rollers (Aves, Coracii). Zool J Linnean Soc. 2009;157(3):586–611.
Ksepka DT, Clarke JA. Primobucco mcgrewi (Aves: Coracii) from the Eocene Green River formation: new anatomical data from the earliest constrained record of stem rollers. J Vertebr Paleontol. 2010;30(1):215–25.
Jarvis ED, Mirarab S, Aberer AJ, Li B, Houde P, Li C, Ho SYW, Faircloth BC, Nabholz B, Howard JT, et al. Whole-genome analyses resolve early branches in the tree of life of modern birds. Science. 2014;346(6215):1320–31.
Prum RO, Berv JS, Dornburg A, Field DJ, Townsend JP, Lemmon EM, Lemmon AR. A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing. Nature. 2015;526(7574):569–73.
Longrich NR, Tokaryk T, Field DJ. Mass extinction of birds at the cretaceous-Paleogene (K-Pg) boundary. Proc Natl Acad Sci. 2011;108(37):15253–7.
Cracraft J, Claramunt S. Conceptual and analytical worldviews shape differences about global avian biogeography. J Biogeogr. 2017;44(4):958-60.
Mayr G. Avian higher level biogeography: southern hemispheric origins or southern hemispheric relicts? J Biogeogr. 2017;44(4):956–8.
Olson SL. A new family of primitive landbirds from the lower Eocene Green River formation of Wyoming. In: Campbell KE, editor. Papers in avian paleontology honoring Pierce Brodkorb, vol. vol. 36: Natural History Museum of Los Angeles County science series; 1992. p. 137–60.
Turner DA. Turacos (Musophagidae). In: del Hoyo J, Elliott A, Sargatal J, editors. Handbook of the Birds of the World Sandgrouse to Cuckoos, vol. vol 4. Lynx Edicions: Barcelona; 1997.
van Tuinen M, Sibley CG, Hedges SB. The early history of modern birds inferred from DNA sequences of nuclear and mitochondrial ribosomal genes. Mol Biol Evol. 2000;17(3):451–7.
Ericson PGP, Anderson CL, Britton T, Elzanowksi A, Johannson US, Källersjö M, Ohlson JI, Parsons TJ, Zuccon D, Mayr G. Diversification of Neoaves: integration of molecular sequence data and fossils. Biol Lett. 2006;2:543–7.
Hackett SJ, Kimball RT, Reddy S, Bowie RCK, Braun EL, Braun MJ, Chojnowski JL, Cox WA, Han K-L, Harshman J, et al. A phylogenomic study of birds reveals their evolutionary history. Science. 2008;320:1763–8.
McCormack JE, Harvey MG, Faircloth BC, Crawford NG, Glenn TC, Brumfield RT. A phylogeny of birds based on over 1,500 loci collected by target enrichment and high-throughput sequencing. PLoS One. 2013;8(1):e54848.
Sibley CG, Ahlquist JE. Phylogeny and classification of birds: a study in molecular evolution. New Haven: Yale University Press; 1990.
Mayr G, Manegold A, Johansson US. Monophyletic groups within ‘higher land birds’ - comparison of morphological and molecular data. J Zool Syst Evol Res. 2003;41:233–48.
Mindell DP, Sorenson MD, Huddleston CJ, Miranda HC, Knight A, Swachuk SJ, Yuri T. Phylogenetic relationships among and within select avian orders based on mitochondrial DNA. In: Mindell DP, editor. Avian molecular evolution and systematics. San Diego, CA: Academic Press; 1997. p. 213–47.
Mayr G, Clarke J. The deep divergences of neornithine birds: a phylogenetic analysis of morphological characters. Cladistics. 2003;19(6):527–53.
Johansson US, Parsons TJ, Irestedt M, Ericson PGP. Clades within the ‘higher land birds’ evaluated by nuclear DNA sequences. J Zool Syst Evol Res. 2001;39(1–2):37–51.
Sorenson MD, Oneal E, García-Moreno J, Mindell DP. More taxa, more characters: the hoatzin problem is still unresolved. Mol Biol Evol. 2003;20(9):1484–99.
Hughes JM. Monophyly and phylogeny of cuckoos (Aves, Cuculidae) inferred from osteological characters. Zool J Linnean Soc. 2000;130:263–307.
Swofford DL. PAUP* Phylogenetic Analysis Using Parsimony (*and Other Methods) v.4.0b10 for Macintosh. Sunderland: Sinauer Associates; 2002.
Ronquist F, Huelsenbeck JP. MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.
Mayr G, Alvarenga H, Mourer-Chauviré C. Out of Africa: fossils shed light on the origin of the hoatzin, an iconic Neotropic bird. Naturwissenschaften. 2011;98(11):961–6.
Ericson PGP, Zuccon D, Ohlson JI, Johansson US, Alvarenga H, Prum RO. Higher level phylogeny and morphological evolution of tyrant flycatchers, cotingas, manakins and their allies (Aves: Tyrannida). Mol Phylogenet Evol. 2006;40:471–83.
Lee MSY. Turtle origins: insights from phylogenetic retrofitting and molecular scaffolds. J Evol Biol. 2013;26:2729–38.
Mayr G, De Pietri VL. Earliest and first northern hemisphere hoatzin fossils substantiate old world origin of a “Neotropic endemic”. Naturwissenschaften. 2014;101(2):143–8.
Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES science gateway for inference of large phylogenetic trees. Gateway Comp Environ Workshop (GCE), 2010. 2010;(Ieee):1–8.
Lewis PO. A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol. 2001;50:913–25.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.
Nesbitt SJ, Ksepka DT, Clarke JA. Podargiform affinities of the enigmatic Fluvioviridavis platyrhamphus and the early diversification of Strisores (“Caprimulgiformes” + Apodiformes). PLoS One. 2011;6(11):e26350.
Bertelli S, Lindow BEK, Dyke GJ, Mayr G. Another charadriiform-like bird from the lower Eocene of Denmark. Paleontol J. 2013;47(11):1282–301.
Ksepka DT, Clarke JA, Nesbitt SJ, Kulp FB, Grande L. Fossil evidence of wing shape in a stem relative of swifts and hummingbirds (Aves, pan-Apodiformes). Proc R Soc B Biol Sci. 2013;280(1761):20130580.
Mesquite: a modular system for evolutionary analysis. Version 2.75.
Jetz W, Thomas GH, Joy JB, Hartmann K, Mooers AO. The global diversity of birds in space and time. Nature. 2012;491:444–8.
Brunet J. Oiseaux de l’Éocène supérieur du bassin de Paris. Ann Paleontol. 1970;56:3–57.
Mayr G. Two-phase extinction of “southern hemispheric” birds in the Cenozoic of Europe and the origin of the Neotropic avifauna. Palaeobiodiversity and Palaeoenvironments. 2011;91:325–33.
Mayr G. Eocene fossils and the early evolution of frogmouths (Podargiformes): further specimens of Masillapodargus and a comparison with Fluvioviridavis. Palaeobio Palaeoenviron. 2015;95(4):587–96.
Team RDC. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2003.
Revell LJ. Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3:217–23.
Yu Y, Harris AJ, Blair C, He X. RASP (reconstruct ancestral state in phylogenies): a tool for historical biogeography. Mol Phylogenet Evol. 2015;87:46–9.
Yu Y, Harris A, He X. S-DIVA (statistical dispersal-Vicariance analysis): a tool for inferring biogeographic histories. Mol Phylogenet Evol. 2010;56(2):848–50.
Landis MJ, Matzke NJ, Moore BR, Huelsenbeck JP. Bayesian analysis of biogeography when the number of areas is large. Syst Biol. 2013;62(6):789–804.
Ronquist F, Cannatella D. Dispersal-vicariance analysis: a new approach to the quantification of historical biogeography. Syst Biol. 1997;46(1):195–203.
GPlates [Comput Softw] In., v.1.5.0 build 16091 EDN Retrieved from https://www.gplates.org/. 2015.
Seton M, Müller R, Zahirovic S, Gaina C, Torsvik T, Shephard G, Talsma A, Gurnis M, Turner M, Maus S. Global continental and ocean basin reconstructions since 200Ma. Earth Sci Rev. 2012;113(3):212–70.
Bivand R, Keitt T, Rowlingson B. rgdal: bindings for the geospatial data abstraction library. R package version. 2014:08–16.
Ronquist F, Klopfstein S, Vilhelmsen L, Schulmeister S, Murray DL, Rasnitsyn AP. A total-evidence approach to dating with fossils, applied to the early radiation of the hymenoptera. Syst Biol. 2012;61:973–99.
Hsiang AY, Field DJ, Webster TH, Behlke ADB, Davis MB, Racicot RA, Gauthier JA. The origin of snakes: revealing the ecology, behavior, and evolutionary history of early snakes using genomics, phenomics, and the fossil record. BMC Evol Biol. 2015;15(1):1–22.
Worthy TH, Degrange FJ, Handley WD, Lee MSY. The evolution of giant flightless birds and novel phylogenetic relationships for extinct fowl (Aves, Galloanseres). Royal Soc Open Sci. 2017;4(10):170975.
Parham JF, Donoghue PCJ, Bell CJ, Calway TD, Head JJ, Holroyd PA, Inoue J, Irmis RB, Joyce WG, Ksepka DT, et al. Best practices for justifying fossil calibrations. Syst Biol. 2012;61:346–59.
Mayr G. A specimen of Eocuculus Chandler, 1999 (Aves, ? Cuculidae) from the early Oligocene of France. Geobios. 2006;39(6):865–72.
Chandler RM. Fossil birds of Florissant, Colorado: with a description of a new genus and species of cuckoo. Geol Res Div Tech Rep. 1999;99:49–52.
Smith ND, Grande L, Clarke JA. A new species of Threskiornithidae-like bird (Aves, Ciconiiformes) from the Green River formation (Eocene) of Wyoming. J Vertebr Paleontol. 2013;33(2):363–81.
Mayr G. The origins of crown group birds: molecules and fossils. Palaeontology. 2014;57(2):231–42.
Suh A, Smeds L, Ellegren H. The dynamics of incomplete lineage sorting across the ancient adaptive radiation of neoavian birds. PLoS Biol. 2015;13(8):e1002224.
Ksepka DT. Broken gears in the avian molecular clock: new phylogenetic analyses support stem galliform status for Gallinuloides wyomingensis and rallid affinities for Amitabha urbsinterdictensis. Cladistics. 2009;25(2):173–97.
Pyron RA. Divergence time estimation using fossils as terminal taxa and the origins of Lissamphibia. Syst Biol. 2011;60(4):466–81.
Field DJ, Lynner C, Brown C, Darroch SAF. Skeletal correlates for body mass estimation in modern and fossil flying birds. PLoS One. 2013;8(11):e82000.
Dunning Jr, JB. CRC handbook of avian body masses. CRC Press; 2007.
Sibley CG, Ahlquist JE. The relationships of the hoatzin. Auk. 1973;90(1):1–13.
de Queiroz K, Good DA. The scleral ossicles of Opisthocomus and their phylogenetic significance. Auk. 1988;105(1):29–35.
Brush AH. Comparison of egg-white proteins: effect of electrophoretic conditions. Biochem Syst Ecol. 1979;7(2):155–65.
Rasmussen DT, Olson SL, Simons EL. Fossil birds from the Oligocene Jebel Qatrani formation, Fayum Province, Egypt. Smithson Contrib Paleobiol. 1987;62:1–20.
Field DJ. Preliminary paleoecological insights from the Pliocene avifauna of Kanapoi, Kenya: implications for the ecology of Australopithecus anamensis. J Hum Evol. 2017; https://www.sciencedirect.com/science/article/pii/S0047248417303494.
Ballmann P. Ein neuer vertreter der Musophagidae (Aves) aus dem Chattium von Gaimersheim bei Ingolstadt (Bayern). Mitteilungen der Bayerischen Staatssammlung für Paläontologie und Histor. Geologie. 1970;10:271–5.
Feduccia A. Explosive evolution in tertiary birds and mammals. Science. 1995;267:637–8.
Ksepka DT, Phillips MJ. Avian diversification patterns across the K-Pg boundary: influence of calibrations, datasets, and model misspecification. Ann Mo Bot Gard. 2015;100(4):300–28.
Field DJ, Bercovici A, Berv JS, Dunn RE, Fastovsky DE, Lyson TR, Vajda V, Gauthier JA. Early evolution of modern birds structured by global forest collapse at the end-cretaceous mass extinction. Current Biol. 2018;28(11):1825–31.
Mayr G. Pumiliornis tessellatus Mayr, 1999 revisited — new data on the osteology and possible phylogenetic affinities of an enigmatic middle Eocene bird. Paläontol Z. 2008;82(3):247–53.
Sorenson MD, Payne RB. Molecular systematics: cuckoo phylogeny inferred from mitochondrial DNA sequences. Pp. 68–94 in R.B. Payne. Bird Families of the World: Cuckoos. Oxford University Press; 2005.
Peters DS. Zoogeographical relationships of the Eocene avifauna from Messel (Germany). In: Bell BD, Cossee RO, JEC F, Heather BD, Hitchmough RA, CJR R, Williams MJ, editors. Acta XX congressus internationalis ornithologici. Christchurch: New Zealand: Ornithological Congress Trust Board; 1991. p. 572–7.
Blondel J, Mourer-Chauviré C. Evolution and history of the western Palaearctic avifauna. Trends Ecol Evol. 13(12):488–92.
McKenna MC. Eocene paleolatitude, climate and mammals of Ellesmere Island. Palaeogeogr Palaeoclimatol Palaeoecol. 1980;30:349–62.
McKenna MC. Cenozoic paleogeography of North Atlantic land bridges. In Structure and development of the Greenland-Scotland Ridge (pp. 351-399). Boston: Springer; 1983.
McKenna MC, Bott MHP, Saxov S, Talwani M, Thiede J. Structure and development of the Greenland-Scotland ridge: new methods and concepts. In: Bott M, Saxov S, Talwani M, Thiede J, editors. Structure and development of the Greenland-Scotland ridge: new methods and concepts, vol. vol. 8: Springer Science & Business Media; 1983. p. 351.
Prothero DR. The Eocene–Oligocene transition: paradise lost. New York: Columbia University Press; 1994.
Lindow BEK, Dyke GJ. Bird evolution in the Eocene: climate change in Europe and a Danish fossil fauna. Biol Rev. 2006;81(04):483–99.
Smith KT. A new lizard assemblage from the earliest Eocene (zone Wa0) of the Bighorn Basin, Wyoming, USA: biogeography during the warmest interval of the Cenozoic. J Syst Palaeontol. 2009;7(3):299–358.
Schaal S, Ziegler W. Messel – ein Schaufenster in die Geschichte der Erde und des Lebens. Frankfurt: Kramer; 1988.
Smith ME, Chamberlain KR, Singer BS, Carroll AR. Eocene clocks agree: coeval 40Ar/39Ar, U-Pb, and astronomical ages from the Green River formation. Geology. 2010;38(6):527–30.
Kurochkin EN. A survey of the Paleogene birds of Asia. Smithson Contrib Paleobiol. 1976;27:75–86.
Payne RB: Family Cuculidae (cuckoos). In: Handbook of the Birds of the World, Sandgrouse to Cuckoos. del Hoyo J, Elliott A, Sargatal J. Barcelona: Lynx Edicions; 1997:vol 4 508–607.
Weigel RD. Oligocene birds from Saskatchewan. Q J Florida Acad Sci. 1963;26(3):257–62.
Feduccia A. Avian extinction at the end of the cretaceous: assessing the magnitude and subsequent explosive radiation. Cretac Res. 2014;50:1–15.
Lindow BE. Bird Evolution Across the K–Pg Boundary and the Basal Neornithine Diversification. In Living Dinosaurs: The Evolutionary History of Modern Birds (pp. 338–354). Chichester: Wiley; 2011.
Mayr G. The age of the crown group of passerine birds and its evolutionary significance - molecular calibrations versus the fossil record. Syst Biodivers. 2013;11:7–13.
Berv JS, Field DJ. Genomic signature of an avian Lilliput effect across the K-Pg extinction. Syst Biol. 2018;67(1):1–13.
Ksepka DT, Stidham T, Williamson T. Early Paleocene landbird supports rapid phylogenetic and morphological diversification of crown birds after the K-Pg mass extinction. Proc Natl Acad Sci. 2017;114(30):8047–52.
Field DJ. Big-time insights from a tiny bird fossil. Proc Natl Acad Sci. 2017;114(30):7750–2.
Müller J, Reisz RR. Four well-constrained calibration points from the vertebrate fossil record for molecular clock estimates. BioEssays. 2005;27(10):1069–75.
Parham J, Irmis RB: Caveats on the use of fossil calibrations for molecular dating: a comment on near. et al. Am Nat 2008, 171(1):132–136.
Ksepka DT, Benton MJ, Carrano MT, Gandolfo MA, Head JJ, Hermsen EJ, Joyce WG, Lamm KS, Patane JSL, Phillips MJ, et al. Synthesizing and databasing fossil calibrations: divergence dating and beyond. Biol Lett. 2011;7(6):801–3.
Doyle JA, Donoghue MJ. Paleobotanical perspectives on vascular plant evolution: the importance of fossils in elucidating seed plant phylogeny and macroevolution. Rev Palaeobot Palynol. 1987;50(1):63–95.
Gauthier J, Kluge A, Rowe T. Amniote phylogeny and the importance of fossils. Cladistics. 1988;4(2):105–209.
Donoghue MJ, Doyle JA, Gauthier J, Kluge AG, Rowe T. The importance of fossils in phylogeny reconstruction. Annu Rev Ecol Syst. 1989;20:431–60.
Lee MSY. Uninformative characters and apparent conflict between molecules and morphology. Mol Biol Evol. 2001;18(4):676–80.
Smith ND, Turner AH. Morphology's role in phylogeny reconstruction: perspectives from paleontology. Syst Biol. 2005;54(1):166–73.
Finarelli JA, Flynn JJ. Ancestral state reconstruction of body size in the Caniformia (Carnivora, Mammalia): the effects of incorporating data from the fossil record. Syst Biol. 2006;55(2):301–13.
Slater GJ, Harmon LJ, Alfaro ME. Integrating fossils with molecular phylogenies improves inference of trait evolution. Evolution. 2012;66(12):3931–44.
Mitchell JS. Extant-only comparative methods fail to recover the disparity preserved in the bird fossil record. Evolution. 2015;69(9):2414–24.
Lieberman BS. Paleobiogeography: using fossils to study global change, plate tectonics, and evolution. Top Geobiol. 2000;16:1–208.
Lieberman BS. Phylogenetic biogeography with and without the fossil record: gauging the effects of extinction and paleontological incompleteness. Palaeogeogr Palaeoclimatol Palaeoecol. 2002;178(1–2):39–52.
Lieberman BS. Paleobiogeography: the relevance of fossils to biogeography. Annu Rev Ecol Evol Syst. 2003;34:51–69.
Murray AM. Eocene cichlid fishes from Tanzania, East Africa. J Vertebr Paleontol. 2001;20(4):651–64.
James HF. Paleogene fossils and the radiation of modern birds. Auk. 2005;122(4):1049–54.
Mayr G. The Paleogene fossil record of birds in Europe. Biol Rev. 2005;80(04):515.
Baker AJ, Pereira SL, Paton TA. Phylogenetic relationships and divergence times of Charadriiformes genera: multigene evidence for the cretaceous origin of at least 14 clades of shorebirds. Biol Lett. 2007;3(2):205–10.
Barker FK, Cibois A, Schikler P, Felsenstein J, Cracraft J. Phylogeny and diversification of the largest avian radiation. Proc Natl Acad Sci U S A. 2004;101:11040–5.
Mitchell KJ, Cooper A, Phillips MJ. Comment on “whole-genome analyses resolve early branches in the tree of life of modern birds”. Science. 2015;349(6255):1460.
Nesbitt SJ, Clarke JA. The anatomy and taxonomy of the exquisitely preserved Green River formation (early Eocene) Lithornithids (Aves) and the relationships of Lithornithidae. Bull Am Mus Nat Hist. 2016;1–91.
We gratefully acknowledge collections assistance from Helen James and Mark Florence for making USNM 336261 available for study, Brian Schmidt and Chris Milenski for providing access to extant comparative material at the USNM, and Kristof Zyskowski for providing access to extant bird skeletons at the YPM. Discussion with Helen James, Jacques Gauthier, Bhart-Anjan Bhullar, Richard Prum, Michael Donoghue, Daniel Ksepka, and Krister Smith improved the quality of the manuscript. We are especially grateful to Gerald Mayr for discussion on character codings and for providing the character/taxon matrix from Mayr et al. 2011, to Sterling Nesbitt for discussion of Palaeogene stem tinamous, and to Lawrence M. Witmer and Ryan C. Ridgely of the Ohio University MicroCT Scanning Facility for use of the Musophaga rossae μCT scan. We thank the International Turaco Society for providing illustrations.
DJF was supported by NSERC (PGS-D), a Smithsonian Predoctoral Fellowship, and a 50th Anniversary Prize Fellowship from the University of Bath.
Availability of data and materials
All supporting data are included in the additional files.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional information on materials and methods. (DOCX 142 kb)
Nexus file for phylogenetic analysis. (ZIP 7 kb)
Table S1. Character/taxon matrix used for the phylogenetic analyses. Character descriptions follow Mayr 2011 for characters 6, 15, 22, 31, 32, 44, 57, 58, 65, 67, 85, 94, 100, 105, 106, 149, 150, and 151. For the remaining characters, descriptions follow Mayr & Clarke 2003, except newly added characters 152 and 153 as described in the main text. (PDF 123 kb)
Table S2. Hindlimb measurements for fossil and recent Otidimorphae. (PDF 70 kb)
Figure S1. Relative likelihood of a Gondwanan vs. Laurasian avian common ancestor based on alternative parameterization of the historical biogeographic analyses. (PNG 211 kb)