Genomic Signature of an Avian Lilliput Effect across the K-Pg Extinction

Drivers of avian genomic change revealed by evolutionary rate decomposition

Modern birds have diversified into a striking array of forms, behaviours and ecological roles. Analyses of molecular evolutionary rates can reveal the links between genomic and phenotypic change1-4, but disentangling the drivers of rate variation at the whole-genome scale has been difficult. Using comprehensive estimates of traits and evolutionary rates across a family-level phylogeny of birds5,6, we find that genome-wide mutation rates across lineages are predominantly explained by clutch size and generation length, whereas rate variation across genes is driven by the content of guanine and cytosine. Here, to find the subsets of genes and lineages that dominate evolutionary rate variation in birds, we estimated the influence of individual lineages on decomposed axes of gene-specific evolutionary rates. We find that most of the rate variation occurs along recent branches of the tree, associated with present-day families of birds. Additional tests on axes of rate variation show rapid changes in microchromosomes immediately after the Cretaceous-Palaeogene transition. These apparent pulses of evolution are consistent with major changes in the genetic machineries for meiosis, heart performance, and RNA splicing, surveillance and translation, and correlate with the ecological diversity reflected in increased tarsus length. Collectively, our analyses paint a nuanced picture of avian evolution, revealing that the ancestors of the most diverse lineages of birds underwent major genomic changes related to mutation, gene usage and niche expansion in the early Palaeogene period.

Cretaceous Antarctic bird skull elucidates early avian ecological diversity

Fossils representing Cretaceous lineages of crown clade birds (Aves) are exceptionally rare but are crucial to elucidating major ecological shifts across early avian divergences. Among the earliest known putative crown birds is Vegavis iaai1-5, a foot-propelled diver from the latest Cretaceous (69.2-68.4 million years ago)6 of Antarctica with controversial phylogenetic affinities2,7-10. Initially recovered by phylogenetic analyses as a stem anatid (ducks and closely related species)1,2,11, Vegavis has since been recovered as a stem member of Anseriformes (waterfowl)7-9, or outside Aves altogether10. Here we report a new, nearly complete skull of Vegavis that provides new insight into its feeding ecology and exhibits morphologies that support placement among waterfowl within crown-group birds. Vegavis has an avian beak (absence of teeth and reduced maxilla) and brain shape (hyperinflated cerebrum and ventrally shifted optic lobes). The temporal fossa is well excavated and expansive, indicating that this bird had hypertrophied jaw musculature. The beak is narrow and pointed, and the mandible lacks retroarticular processes. Together, these features comprise a feeding apparatus unlike that of any other known anseriform but like that of other extant birds that capture prey underwater (for example, grebes and loons). The Cretaceous occurrence of Vegavis, with a feeding ecology unique among known Galloanserae (waterfowl and landfowl), is further indication that the earliest anseriform divergences were marked by evolutionary experiments unrepresented in the extant diversity3,11-13.

Shell Constraints on Evolutionary Body Size-Limb Size Allometry Can Explain Morphological Conservatism in the Turtle Body Plan

Turtles are a small clade of vertebrates despite having existed since the Late Triassic. Turtles have a conservative body plan relative to other amniotes, characterized by the presence of a shell and quadrupedality. This morphology is even retained in strong ecological specialists, such as sea turtles, which are secondarily adapted to marine locomotion by strong allometric scaling in their hands. It is possible that the body plan of turtles is strongly influenced by the presence of the shell, acting as a constraint to achieving greater diversity of body forms. Here, we explore the evolutionary allometric relationships of fore- and hindlimb stylopodia (i.e., humerus and femur) with one another as well as their relationship with shell size (carapace length) to assess evidence of constraint. All turtles, including Triassic shelled stem turtles, have near-isometric relationships that do not vary strongly between clades, and evolve at slow evolutionary rates. This indeed indicates that body proportions of turtles are constrained to a narrow range of possibilities. Minor allometric deviations are seen in highly aquatic sea turtles and softshell turtles, which modified their shells by bone losses. Our allometric regressions allow accurate body size estimations for fossils. Several independent sea turtle lineages converged on maximum sizes of 2.2 m of shell length, which may be a biological maximum for the group.

Genome and life-history evolution link bird diversification to the end-Cretaceous mass extinction

Complex patterns of genome evolution associated with the end-Cretaceous [Cretaceous-Paleogene (K-Pg)] mass extinction limit our understanding of the early evolutionary history of modern birds. Here, we analyzed patterns of avian molecular evolution and identified distinct macroevolutionary regimes across exons, introns, untranslated regions, and mitochondrial genomes. Bird clades originating near the K-Pg boundary exhibited numerous shifts in the mode of molecular evolution, suggesting a burst of genomic heterogeneity at this point in Earth's history. These inferred shifts in substitution patterns were closely related to evolutionary shifts in developmental mode, adult body mass, and patterns of metabolic scaling. Our results suggest that the end-Cretaceous mass extinction triggered integrated patterns of evolution across avian genomes, physiology, and life history near the dawn of the modern bird radiation.

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae

BACKGROUND

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems.

RESULTS

We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary.

CONCLUSIONS

Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Paleoneurology of stem palaeognaths clarifies the plesiomorphic condition of the crown bird central nervous system

Lithornithidae, an assemblage of volant Palaeogene fossil birds, provide our clearest insights into the early evolutionary history of Palaeognathae, the clade that today includes the flightless ratites and volant tinamous. The neotype specimen of Lithornis vulturinus, from the early Eocene (approximately 53 million years ago) of Europe, includes a partial neurocranium that has never been thoroughly investigated. Here, we describe these cranial remains including the nearly complete digital endocasts of the brain and bony labyrinth. The telencephalon of Lithornis is expanded and its optic lobes are ventrally shifted, as is typical for crown birds. The foramen magnum is positioned caudally, rather than flexed ventrally as in some crown birds, with the optic lobes, cerebellum, and foramen magnum shifted further ventrally. The overall brain shape is similar to that of tinamous, the only extant clade of flying palaeognaths, suggesting that several aspects of tinamou neuroanatomy may have been evolutionarily conserved since at least the early Cenozoic. The estimated ratio of the optic lobe's surface area relative to the total brain suggests a diurnal ecology. Lithornis may provide the clearest insights to date into the neuroanatomy of the ancestral crown bird, combining an ancestrally unflexed brain with a caudally oriented connection with the spinal cord, a moderately enlarged telencephalon, and ventrally shifted, enlarged optic lobes.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

Complexity of avian evolution revealed by family-level genomes

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

Tip dating and Bayes factors provide insight into the divergences of crown bird clades across the end-Cretaceous mass extinction

The origin of crown birds (Neornithes) remains contentious owing to conflicting divergence time hypotheses obtained from alternative sources of data. The fossil record suggests limited diversification of Neornithes in the Late Mesozoic and a substantial radiation in the aftermath of the Cretaceous-Palaeogene (K-Pg) mass extinction, approximately 66 Ma. Molecular clock studies, however, have yielded estimates for neornithine origins ranging from the Early Cretaceous (130 Ma) to less than 10 Myr before the K-Pg. We use Bayes factors to compare the fit of node ages from different molecular clock studies to an independent morphological dataset. Our results allow us to reject scenarios of crown bird origins deep in the Early Cretaceous, as well as an origin of crown birds within the last 10 Myr of the Cretaceous. The scenario best supported by our analyses is one where Neornithes originated between the Early and Late Cretaceous (ca 100 Ma), while numerous divergences within major neoavian clades either span or postdate the K-Pg. This study affirms the importance of the K-Pg on the diversification of modern birds, and the potential of combined-evidence tip-dating analyses to illuminate recalcitrant 'rocks versus clocks' debates.

A juvenile bird with possible crown-group affinities from a dinosaur-rich Cretaceous ecosystem in North America

BACKGROUND

Living birds comprise the most speciose and anatomically diverse clade of flying vertebrates, but their poor early fossil record and the lack of resolution around the relationships of the major clades have greatly obscured extant avian origins.

RESULTS

Here, I describe a Late Cretaceous bird from North America based on a fragmentary skeleton that includes cranial material and portions of the forelimb, hindlimb, and foot and is identified as a juvenile based on bone surface texture. Several features unite this specimen with crown Aves, but its juvenile status precludes the recognition of a distinct taxon. The North American provenance of the specimen supports a cosmopolitan distribution of early crown birds, clashes with the hypothesized southern hemisphere origins of living birds, and demonstrates that crown birds and their closest relatives coexisted with non-avian dinosaurs that independently converged on avian skeletal anatomy, such as the alvarezsaurids and dromaeosaurids.

CONCLUSIONS

By revealing the ecological and biogeographic context of Cretaceous birds within or near the crown clade, the Lance Formation specimen provides new insights into the contingent nature of crown avian survival through the Cretaceous-Paleogene mass extinction and the subsequent origins of living bird diversity.

Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplications (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome-18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variants (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in testis than ovary. Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

3D atlas of tinamou (Neornithes: Tinamidae) pectoral morphology: Implications for reconstructing the ancestral neornithine flight apparatus

Palaeognathae, the extant avian clade comprising the flightless ratites and flight-capable tinamous (Tinamidae), is the sister group to all other living birds, and recent phylogenetic studies illustrate that tinamous are phylogenetically nested within a paraphyletic assemblage of ratites. As the only extant palaeognaths that have retained the ability to fly, tinamous may provide key information on the nature of the flight apparatus of ancestral crown palaeognaths-and, in turn, crown birds-as well as insight into convergent modifications to the wing apparatus among extant ratite lineages. To reveal new information about the musculoskeletal anatomy of tinamous and facilitate development of computational biomechanical models of tinamou wing function, we generated a three-dimensional musculoskeletal model of the flight apparatus of the extant Andean tinamou (Nothoprocta pentlandii) using diffusible iodine-based contrast-enhanced computed tomography (diceCT). Origins and insertions of the pectoral flight musculature of N. pentlandii are generally consistent with those of other extant volant birds specialized for burst flight, and the entire suite of presumed ancestral neornithine flight muscles are present in N. pentlandii with the exception of the biceps slip. The pectoralis and supracoracoideus muscles are robust, similar to the condition in other extant burst-flying birds such as many extant Galliformes. Contrary to the condition in most extant Neognathae (the sister clade to Palaeognathae), the insertion of the pronator superficialis has a greater distal extent than the pronator profundus, although most other anatomical observations are broadly consistent with the conditions observed in extant neognaths. This work will help form a basis for future comparative studies of the avian musculoskeletal system, with implications for reconstructing the flight apparatus of ancestral crown birds and clarifying musculoskeletal modifications underlying the convergent origins of ratite flightlessness.

Deconstructing age estimates for angiosperms

Estimates of the age of angiosperms from molecular phylogenies vary considerably. As in all estimates of evolutionary timescales from phylogenies, generating these estimates requires assumptions about the rate that molecular sequences are evolving (using clock models) and the time duration of the branches in a phylogeny (using fossil calibrations and branching processes). Often, it is difficult to demonstrate that these assumptions reflect current knowledge of molecular evolution or the fossil record. In this study we re-estimate the age of angiosperms using a minimal set of assumptions, therefore avoiding many of the assumptions inherent to other methods. The age estimates we generate are similar for each of the four datasets analysed, ranging from 130 to 400 Ma, but are far less precise than in previous studies. We demonstrate that this reduction in precision results from making less stringent assumptions about both rate and time, and that the analysed molecular dataset has very little effect on age estimates.

Body size predicts the rate of contemporary morphological change in birds

Variation in evolutionary rates among species is a defining characteristic of the tree of life and may be an important predictor of species' capacities to adapt to rapid environmental change. It is broadly assumed that generation length is an important determinant of microevolutionary rates, and body size is often used as a proxy for generation length. However, body size has myriad biological correlates that could affect evolutionary rates independently from generation length. We leverage two large, independently collected datasets on recent morphological change in birds (52 migratory species breeding in North America and 77 South American resident species) to test how body size and generation length are related to the rates of contemporary morphological change. Both datasets show that birds have declined in body size and increased in wing length over the past 40 y. We found, in both systems, a consistent pattern wherein smaller species declined proportionally faster in body size and increased proportionally faster in wing length. By contrast, generation length explained less variation in evolutionary rates than did body size. Although the mechanisms warrant further investigation, our study demonstrates that body size is an important predictor of contemporary variation in morphological rates of change. Given the correlations between body size and a breadth of morphological, physiological, and ecological traits predicted to mediate phenotypic responses to environmental change, the relationship between body size and rates of phenotypic change should be considered when testing hypotheses about variation in adaptive responses to climate change.

The bacterial origin of mitochondria: Incorrect phylogenies and the importance of metabolic traits

This article provides an updated review on the evolution of mitochondria from bacteria, which were likely related to extant alphaproteobacteria. Particular attention is given to the timeline of oxygen history on Earth and the entwined phases of eukaryotic evolution that produced the animals that still populate our planet. Mitochondria of early-branching unicellular eukaryotes and plants appear to retain partial or vestigial traits that were directly inherited from the alphaproteobacterial ancestors of the organelles. Most of such traits define the current aerobic physiology of mitochondria. Conversely, the anaerobic traits that would be essential in the syntrophic associations postulated for the evolution of eukaryotic cells are scantly present in extant alphaproteobacteria, and therefore cannot help defining from which bacterial lineage the ancestors of mitochondria originated. This question has recently been addressed quantitatively, reaching the novel conclusion that marine bacteria related to Iodidimonas may be the living relatives of protomitochondria. Additional evidence is presented that either support or does not contrast this novel view of the bacterial origin of mitochondria.

Complex macroevolution of pterosaurs

Pterosaurs, the earliest flying tetrapods, are the subject of some recent quantitative macroevolutionary analyses from different perspectives.^1-2 Here, we use an integrative approach involving newly assembled phylogenetic and body size datasets, net diversification rates, morphological rates, and morphological disparity to gain a holistic understanding of the pterosaur macroevolution. The first two parameters are important in quantitative analyses of macroevolution, but they have been rarely used in previous pterosaur studies.^{1,3,4,2,5,6,7,8,9,10,11,12} Our study reveals an ∼115-Ma period-from Early Triassic to Early Cretaceous-of multi-wave increasing net diversification rates and disparity, as well as high morphological rates, followed by an ∼65-Ma period-from Early Cretaceous to the end of the Cretaceous-of mostly negative net diversification rates, decreasing disparity, and relatively low morphological rates in pterosaur evolution. Our study demonstrates the following: (1) body size plays an important role in pterosaur lineage diversification during nearly their whole evolutionary history, and the evolution of locomotion, trophic, and ornamental structures also plays a role in different periods; (2) birds, the other major flying tetrapod group at the time, might have affected pterosaur macroevolution for ∼100 Ma; and (3) different mass extinction events might have affected pterosaur evolution differently. Particularly, the revealed decline in pterosaur biodiversity during the Middle and Late Cretaceous periods provides further support for the possible presence of a biodiversity decline of large-sized terrestrial amniotes starting in the mid-Cretaceous,^13,14 which may have been caused by multiple factors including a global land area decrease during these periods.

Dinosaurs: Comparative Cytogenomics of Their Reptile Cousins and Avian Descendants

Reptiles known as dinosaurs pervade scientific and popular culture, while interest in their genomics has increased since the 1990s. Birds (part of the crown group Reptilia) are living theropod dinosaurs. Chromosome-level genome assemblies cannot be made from long-extinct biological material, but dinosaur genome organization can be inferred through comparative genomics of related extant species. Most reptiles apart from crocodilians have both macro- and microchromosomes; comparative genomics involving molecular cytogenetics and bioinformatics has established chromosomal relationships between many species. The capacity of dinosaurs to survive multiple extinction events is now well established, and birds now have more species in comparison with any other terrestrial vertebrate. This may be due, in part, to their karyotypic features, including a distinctive karyotype of around n = 40 (~10 macro and 30 microchromosomes). Similarity in genome organization in distantly related species suggests that the common avian ancestor had a similar karyotype to e.g., the chicken/emu/zebra finch. The close karyotypic similarity to the soft-shelled turtle (n = 33) suggests that this basic pattern was mostly established before the Testudine-Archosaur divergence, ~255 MYA. That is, dinosaurs most likely had similar karyotypes and their extensive phenotypic variation may have been mediated by increased random chromosome segregation and genetic recombination, which is inherently higher in karyotypes with more and smaller chromosomes.

Cretaceous ornithurine supports a neognathous crown bird ancestor

The bony palate diagnoses the two deepest clades of extant birds: Neognathae and Palaeognathae^1-5. Neognaths exhibit unfused palate bones and generally kinetic skulls, whereas palaeognaths possess comparatively rigid skulls with the pterygoid and palatine fused into a single element, a condition long considered ancestral for crown birds (Neornithes)^3,5-8. However, fossil evidence of palatal remains from taxa close to the origin of Neornithes is scarce, hindering strong inferences regarding the ancestral condition of the neornithine palate. Here we report a new taxon of toothed Late Cretaceous ornithurine bearing a pterygoid that is remarkably similar to those of the extant neognath clade Galloanserae (waterfowl + landfowl). Janavis finalidens, gen. et sp. nov., is generally similar to the well-known Mesozoic ornithurine Ichthyornis in its overall morphology, although Janavis is much larger and exhibits a substantially greater degree of postcranial pneumaticity. We recovered Janavis as the first-known well-represented member of Ichthyornithes other than Ichthyornis, clearly substantiating the persistence of the clade into the latest Cretaceous⁹. Janavis confirms the presence of an anatomically neognathous palate in at least some Mesozoic non-crown ornithurines^10-12, suggesting that pterygoids similar to those of extant Galloanserae may be plesiomorphic for crown birds. Our results, combined with recent evidence on the ichthyornithine palatine¹², overturn longstanding assumptions about the ancestral crown bird palate, and should prompt reevaluation of the purported galloanseran affinities of several bizarre early Cenozoic groups such as the 'pseudotoothed birds' (Pelagornithidae)^13-15.

Orthologous microsatellites, transposable elements, and DNA deletions correlate with generation time and body mass in neoavian birds

The rate of mutation accumulation in germline cells can be affected by cell replication and/or DNA damage, which are further related to life history traits such as generation time and body mass. Leveraging the existing datasets of 233 neoavian bird species, here, we investigated whether generation time and body mass contribute to the interspecific variation of orthologous microsatellite length, transposable element (TE) length, and deletion length and how these genomic attributes affect genome sizes. In nonpasserines, we found that generation time is correlated to both orthologous microsatellite length and TE length, and body mass is negatively correlated to DNA deletions. These patterns are less pronounced in passerines. In all species, we found that DNA deletions relate to genome size similarly as TE length, suggesting a role of body mass dynamics in genome evolution. Our results indicate that generation time and body mass shape the evolution of genomic attributes in neoavian birds.

Comprehensive taxon sampling and vetted fossils help clarify the time tree of shorebirds (Aves, Charadriiformes)

Shorebirds (Charadriiformes) are a globally distributed clade of modern birds and, due to their ecological and morphological disparity, a frequent subject of comparative studies. While molecular phylogenies have been key to establishing the suprafamilial backbone of the charadriiform tree, a number of relationships at both deep and shallow taxonomic levels remain poorly resolved. The timescale of shorebird evolution also remains uncertain as a result of extensive disagreements among the published divergence dating studies, stemming largely from different choices of fossil calibrations. Here, we present the most comprehensive non-supertree phylogeny of shorebirds to date, based on a total-evidence dataset comprising 353 ingroup taxa (90% of all extant or recently extinct species), 27 loci (15 mitochondrial and 12 nuclear), and 69 morphological characters. We further clarify the timeline of charadriiform evolution by time-scaling this phylogeny using a set of 14 up-to-date and thoroughly vetted fossil calibrations. In addition, we assemble a taxonomically restricted 100-locus dataset specifically designed to resolve outstanding problems in higher-level charadriiform phylogeny. In terms of tree topology, our results are largely congruent with previous studies but indicate that some of the conflicts among earlier analyses reflect a genuine signal of pervasive gene tree discordance. Monophyly of the plovers (Charadriidae), the position of the ibisbill (Ibidorhyncha), and the relationships among the five subfamilies of the gulls (Laridae) could not be resolved even with greatly increased locus and taxon sampling. Moreover, several localized regions of uncertainty persist in shallower parts of the tree, including the interrelationships of the true auks (Alcinae) and anarhynchine plovers. Our node-dating and macroevolutionary rate analyses find support for a Paleocene origin of crown-group shorebirds, as well as exceptionally rapid recent radiations of Old World oystercatchers (Haematopodidae) and select genera of gulls. Our study underscores the challenges involved in estimating a comprehensively sampled and carefully calibrated time tree for a diverse avian clade, and highlights areas in need of further research.

The impact of paleoclimatic changes on body size evolution in marine fishes

General rules are useful tools for understanding how organisms evolve. Cope’s rule (tendency to increase in size over evolutionary time) and Bergmann’s rule (tendency to grow to larger sizes in cooler climates) both relate to body size, an important factor that affects the biology, ecology, and physiology of organisms. These rules are well studied in endotherms but remain poorly understood among ectotherms. Here, we show that paleoclimatic changes strongly shaped the trajectory of body size evolution in tetraodontiform fishes. Their body size evolution is explained by both Cope’s and Bergmann’s rules, highlighting the impact of paleoclimatic changes on aquatic organisms, which rely on their environment for temperature regulation and are likely more susceptible than terrestrial vertebrates to climatic changes.

Body size is an important species trait, correlating with life span, fecundity, and other ecological factors. Over Earth’s geological history, climate shifts have occurred, potentially shaping body size evolution in many clades. General rules attempting to summarize body size evolution include Bergmann’s rule, which states that species reach larger sizes in cooler environments and smaller sizes in warmer environments, and Cope’s rule, which poses that lineages tend to increase in size over evolutionary time. Tetraodontiform fishes (including pufferfishes, boxfishes, and ocean sunfishes) provide an extraordinary clade to test these rules in ectotherms owing to their exemplary fossil record and the great disparity in body size observed among extant and fossil species. We examined Bergmann’s and Cope’s rules in this group by combining phylogenomic data (1,103 exon loci from 185 extant species) with 210 anatomical characters coded from both fossil and extant species. We aggregated data layers on paleoclimate and body size from the species examined, and inferred a set of time-calibrated phylogenies using tip-dating approaches for downstream comparative analyses of body size evolution by implementing models that incorporate paleoclimatic information. We found strong support for a temperature-driven model in which increasing body size over time is correlated with decreasing oceanic temperatures. On average, extant tetraodontiforms are two to three times larger than their fossil counterparts, which otherwise evolved during periods of warmer ocean temperatures. These results provide strong support for both Bergmann’s and Cope’s rules, trends that are less studied in marine fishes compared to terrestrial vertebrates and marine invertebrates.

Giant gar from directly above the Cretaceous-Palaeogene boundary suggests healthy freshwater ecosystems existed within thousands of years of the asteroid impact

The Cretaceous-Palaeogene (K-Pg) mass extinction was responsible for the destruction of global ecosystems and loss of approximately three-quarters of species diversity 66 million years ago. Large-bodied land vertebrates suffered high extinction rates, whereas small-bodied vertebrates living in freshwater ecosystems were buffered from the worst effects. Here, we report a new species of large-bodied (1.4-1.5 m) gar based on a complete skeleton from the Williston Basin of North America. The new species was recovered 18 cm above the K-Pg boundary, making it one of the oldest articulated vertebrate fossils from the Cenozoic. The presence of this freshwater macropredator approximately 1.5-2.5 thousand years after the asteroid impact suggests the rapid recovery and reassembly of North American freshwater food webs and ecosystems after the mass extinction.

Macroevolutionary dynamics of climatic niche space

How and why lineages evolve along with niche space as they diversify and adapt to different environments is fundamental to evolution. Progress has been hampered by the difficulties of linking a robust empirical characterization of species niches with flexible evolutionary models that describe their evolution. Consequently, the relative influence of abiotic and biotic factors remains poorly understood. Here, we characterize species' two-dimensional temperature and precipitation niche space occupied (i.e. species niche envelope) as complex geometries and assess their evolution across all Aves using a model that captures heterogeneous evolutionary rates on time-calibrated phylogenies. We find that extant birds coevolved from warm, mesic climatic niches into colder and drier environments and responded to the Cretaceous-Palaeogene (K-Pg) boundary with a dramatic increase in disparity. Contrary to expectations of subsiding rates of niche evolution, our results show that overall rates have increased steadily, with some lineages experiencing exceptionally high evolutionary rates, associated with the colonization of novel niche spaces, and others showing niche stasis. Both competition- and environmental change-driven niche evolution transpire and result in highly heterogeneous rates near the present. Our findings highlight the growing ecological and conservation insights arising from the model-based integration of comprehensive environmental and phylogenetic information.

Investigating the reliability of molecular estimates of evolutionary time when substitution rates and speciation rates vary

Background

An accurate timescale of evolutionary history is essential to testing hypotheses about the influence of historical events and processes, and the timescale for evolution is increasingly derived from analysis of DNA sequences. But variation in the rate of molecular evolution complicates the inference of time from DNA. Evidence is growing for numerous factors, such as life history and habitat, that are linked both to the molecular processes of mutation and fixation and to rates of macroevolutionary diversification. However, the most widely used methods rely on idealised models of rate variation, such as the uncorrelated and autocorrelated clocks, and molecular dating methods are rarely tested against complex models of rate change. One relationship that is not accounted for in molecular dating is the potential for interaction between molecular substitution rates and speciation, a relationship that has been supported by empirical studies in a growing number of taxa. If these relationships are as widespread as current evidence suggests, they may have a significant influence on molecular dates.

Results

We simulate phylogenies and molecular sequences under three different realistic rate variation models—one in which speciation rates and substitution rates both vary but are unlinked, one in which they covary continuously and one punctuated model in which molecular change is concentrated in speciation events, using empirical case studies to parameterise realistic simulations. We test three commonly used “relaxed clock” molecular dating methods against these realistic simulations to explore the degree of error in molecular dates under each model. We find average divergence time inference errors ranging from 12% of node age for the unlinked model when reconstructed under an uncorrelated rate prior using BEAST 2, to up to 91% when sequences evolved under the punctuated model are reconstructed under an autocorrelated prior using PAML.

Conclusions

We demonstrate the potential for substantial errors in molecular dates when both speciation rates and substitution rates vary between lineages. This study highlights the need for tests of molecular dating methods against realistic models of rate variation generated from empirical parameters and known relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02015-8.

Diversification Rate is Associated with Rate of Molecular Evolution in Ray-Finned Fish (Actinopterygii)

Understanding the factors that drive diversification of taxa across the tree of life is a key focus of macroevolutionary research. While the effects of life history, ecology, climate and geography on diversity have been studied for many taxa, the relationship between molecular evolution and diversification has received less attention. However, correlations between rates of molecular evolution and diversification rate have been detected in a range of taxa, including reptiles, plants and birds. A correlation between rates of molecular evolution and diversification rate is a prediction of several evolutionary theories, including the evolutionary speed hypothesis which links variation in mutation rates to differences in speciation rates. If it is widespread, such correlations could also have significant practical impacts, if they are not adequately accounted for in phylogenetic inference of evolutionary rates and timescales. Ray-finned fish (Actinopterygii) offer a prime target to test for this relationship due to their extreme variation in clade size suggesting a wide range of diversification rates. We employ both a sister-pairs approach and a whole-tree approach to test for correlations between substitution rate and net diversification. We also collect life history and ecological trait data and account for potential confounding factors including body size, latitude, max depth and reef association. We find evidence to support a relationship between diversification and synonymous rates of nuclear evolution across two published backbone phylogenies, as well as weak evidence for a relationship between mitochondrial nonsynonymous rates and diversification at the genus level.

exTREEmaTIME: a method for incorporating uncertainty into divergence time estimates

We present a method of divergence time estimation (exTREEmaTIME) that aims to effectively account for uncertainty in divergence time estimates. The method requires a minimal set of assumptions, and, based on these assumptions, estimates the oldest possible divergence times and youngest possible divergence times that are consistent with the assumptions. We use a series of simulations and empirical analyses to illustrate that exTREEmaTIME is effective at representing uncertainty. We then describe how exTREEmaTIME can act as a basis to determine the implications of the more stringent assumptions that are incorporated into other methods of divergence time estimation that produce more precise estimates. This is critically important given that many of the assumptions that are incorporated into these methods are highly complex, difficult to justify biologically, and as such can lead to estimates that are highly inaccurate. This article has an associated First Person interview with the first author of the paper.

Predominantly Eastward Long-Distance Dispersal in Pantropical Ochnaceae Inferred From Ancestral Range Estimation and Phylogenomics

Ochnaceae is a pantropical family with multiple transoceanic disjunctions at deep and shallow levels. Earlier attempts to unravel the processes that led to such biogeographic patterns suffered from insufficient phylogenetic resolution and unclear delimitation of some of the genera. In the present study, we estimated divergence time and ancestral ranges based on a phylogenomic framework with a well-resolved phylogenetic backbone to tackle issues of the timing and direction of dispersal that may explain the modern global distribution of Ochnaceae. The nuclear data provided the more robust framework for divergence time estimation compared to the plastome-scale data, although differences in the inferred clade ages were mostly small. While Ochnaceae most likely originated in West Gondwana during the Late Cretaceous, all crown-group disjunctions are inferred as dispersal-based, most of them as transoceanic long-distance dispersal (LDD) during the Cenozoic. All LDDs occurred in an eastward direction except for the SE Asian clade of Sauvagesieae, which was founded by trans-Pacific dispersal from South America. The most species-rich clade by far, Ochninae, originated from either a widespread neotropical-African ancestor or a solely neotropical ancestor which then dispersed to Africa. The ancestors of this clade then diversified in Africa, followed by subsequent dispersal to the Malagasy region and tropical Asia on multiple instances in three genera during the Miocene-Pliocene. In particular, Ochna might have used the South Arabian land corridor to reach South Asia. Thus, the pantropical distribution of Ochnaceae is the result of LDD either transoceanic or via land bridges/corridors, whereas vicariance might have played a role only along the stem of the family.

Dynamic Patterns of Sex Chromosome Evolution in Neognath Birds: Many Independent Barriers to Recombination at the ATP5F1A Locus

Avian sex chromosomes evolved after the divergence of birds and crocodilians from their common ancestor, so they are younger than the better-studied chromosomes of mammals. It has long been recognized that there may have been several stages to the evolution of avian sex chromosomes. For example, the CHD1 undergoes recombination in paleognaths but not neognaths. Genome assemblies have suggested that there may be variation in the timing of barriers to recombination among Neognathae, but there remains little understanding of the extent of this variability. Here, we look at partial sequences of ATP5F1A, which is on the avian Z and W chromosomes. It is known that recombination of this gene has independently ceased in Galliformes, Anseriformes, and at least five neoavian orders, but whether there are other independent cessations of recombination among Neoaves is not understood. We analyzed a combination of data extracted from published chromosomal-level genomes with data collected using PCR and cloning to identify Z and W copies in 22 orders. Our results suggest that there may be at least 19 independent cessations of recombination within Neognathae, and 3 clades that may still be undergoing recombination (or have only recently ceased recombination). Analyses of ATP5F1A protein sequences revealed an increased amino acid substitution rate for W chromosome gametologs, suggesting relaxed purifying selection on the W chromosome. Supporting this hypothesis, we found that the increased substitution rate was particularly pronounced for buried residues, which are expected to be more strongly constrained by purifying selection. This highlights the dynamic nature of avian sex chromosomes, and that this level of variation among clades means they should be a good system to understand sex chromosome evolution.

Forty new specimens of Ichthyornis provide unprecedented insight into the postcranial morphology of crownward stem group birds

Ichthyornis has long been recognized as a pivotally important fossil taxon for understanding the latest stages of the dinosaur-bird transition, but little significant new postcranial material has been brought to light since initial descriptions of partial skeletons in the 19^th Century. Here, we present new information on the postcranial morphology of Ichthyornis from 40 previously undescribed specimens, providing the most complete morphological assessment of the postcranial skeleton of Ichthyornis to date. The new material includes four partially complete skeletons and numerous well-preserved isolated elements, enabling new anatomical observations such as muscle attachments previously undescribed for Mesozoic euornitheans. Among the elements that were previously unknown or poorly represented for Ichthyornis, the new specimens include an almost-complete axial series, a hypocleideum-bearing furcula, radial carpal bones, fibulae, a complete tarsometatarsus bearing a rudimentary hypotarsus, and one of the first-known nearly complete three-dimensional sterna from a Mesozoic avialan. Several pedal phalanges are preserved, revealing a remarkably enlarged pes presumably related to foot-propelled swimming. Although diagnosable as Ichthyornis, the new specimens exhibit a substantial degree of morphological variation, some of which may relate to ontogenetic changes. Phylogenetic analyses incorporating our new data and employing alternative morphological datasets recover Ichthyornis stemward of Hesperornithes and Iaceornis, in line with some recent hypotheses regarding the topology of the crownward-most portion of the avian stem group, and we establish phylogenetically-defined clade names for relevant avialan subclades to help facilitate consistent discourse in future work. The new information provided by these specimens improves our understanding of morphological evolution among the crownward-most non-neornithine avialans immediately preceding the origin of crown group birds.

Body mass and geographic distribution determined the evolution of the wing flight-feather molt strategy in the Neornithes lineage

The evolutionary history of many organisms is characterized by major changes in morphology and distribution. Specifically, alterations of body mass and geographic distribution may profoundly influence organismal life-history traits. Here, we reconstructed the evolutionary history of flight-feather molt strategy using data from 1,808 Neornithes species. Our analysis suggests that the ancestral molt strategy of first-year birds was partial or entirely absent, and that complete wing flight-feather molt in first-year birds first evolved in the late Eocene and Oligocene (25-40 Ma), at least 30 Myr after birds first evolved. Complete flight-feather molt occurred mainly at equatorial latitudes and in relatively low body mass species, following a diversification of body mass within the lineage. We conclude that both body mass and geographic distribution shaped the evolution of molt strategies and propose that the evolutionary transition towards complete juvenile molt in the Neornithes is a novel, relatively late adaptation.

Causes and Consequences of Apparent Timescaling Across All Estimated Evolutionary Rates

Evolutionary rates play a central role in connecting micro- and macroevolution. All evolutionary rate estimates, including rates of molecular evolution, trait evolution, and lineage diversification, share a similar scaling pattern with time: The highest rates are those measured over the shortest time interval. This creates a disconnect between micro- and macroevolution, although the pattern is the opposite of what some might expect: Patterns of change over short timescales predict that evolution has tremendous potential to create variation and that potential is barely tapped by macroevolution. In this review, we discuss this shared scaling pattern across evolutionary rates. We break down possible explanations for scaling into two categories, estimation error and model misspecification, and discuss how both apply to each type of rate. We also discuss the consequences of this ubiquitous pattern, which can lead to unexpected results when comparing ratesover different timescales. Finally, after addressing purely statistical concerns, we explore a few possibilities for a shared unifying explanation across the three types of rates that results from a failure to fully understand and account for how biological processes scale over time.

Ecological selectivity and the evolution of mammalian substrate preference across the K-Pg boundary

The Cretaceous-Paleogene (K-Pg) mass extinction 66 million years ago was characterized by a worldwide ecological catastrophe and rapid species turnover. Large-scale devastation of forested environments resulting from the Chicxulub asteroid impact likely influenced the evolutionary trajectories of multiple clades in terrestrial environments, and it has been hypothesized to have biased survivorship in favour of nonarboreal lineages across the K-Pg boundary. Here, we evaluate patterns of substrate preferences across the K-Pg boundary among crown group mammals, a group that underwent rapid diversification following the mass extinction. Using Bayesian, likelihood, and parsimony reconstructions, we identify patterns of mammalian ecological selectivity that are broadly similar to those previously hypothesized for birds. Models based on extant taxa indicate predominant K-Pg survivorship among semi- or nonarboreal taxa, followed by numerous independent transitions to arboreality in the early Cenozoic. However, contrary to the predominant signal, some or all members of total-clade Euarchonta (Primates + Dermoptera + Scandentia) appear to have maintained arboreal habits across the K-Pg boundary, suggesting ecological flexibility during an interval of global habitat instability. We further observe a pronounced shift in character state transitions away from plesiomorphic arboreality associated with the K-Pg transition. Our findings are consistent with the hypothesis that predominantly nonarboreal taxa preferentially survived the end-Cretaceous mass extinction, and emphasize the pivotal influence of the K-Pg transition in shaping the early evolutionary trajectories of extant terrestrial vertebrates.

Disentangling the avian altricial-precocial spectrum: Quantitative assessment of developmental mode, phylogenetic signal, and dimensionality

The altricial-precocial spectrum describes patterns of variation in avian developmental mode that greatly influence avian life histories. Appraising a given species' position on this spectrum is therefore fundamental to understanding patterns of avian life history evolution. However, evaluating avian developmental mode remains a relatively subjective task reliant on untested assumptions, including the notion that developmental strategies are distributed along a single dimension of statistical variation. Here, we present a quantitative multivariate framework that objectively discriminates among meaningfully different modes of avian development. We gathered information on seven hatchling and post-hatching traits for up to 4000 extant bird species, and find that most traits related to developmental mode show high phylogenetic signal and little intraclade variation, allowing unknown values to be reliably interpolated. Principal component analyses (PCAs) of these traits illustrate that most variation in hatchling state can be quantified along one dimension of trait space. However, our PCAs also reveal an important second dimension explaining variation in post-hatching behavior, enabling factors related to hatchling state and post-hatching behavior to be disentangled. In order to facilitate future macroevolutionary studies of variation in avian developmental strategies, as well as explorations of covariation between developmental mode and other aspects of avian biology, we present PC scores for 9993 extant avian species.

Evolution and dispersal of snakes across the Cretaceous-Paleogene mass extinction

Mass extinctions have repeatedly shaped global biodiversity. The Cretaceous-Paleogene (K-Pg) mass extinction caused the demise of numerous vertebrate groups, and its aftermath saw the rapid diversification of surviving mammals, birds, frogs, and teleost fishes. However, the effects of the K-Pg extinction on the evolution of snakes-a major clade of predators comprising over 3,700 living species-remains poorly understood. Here, we combine an extensive molecular dataset with phylogenetically and stratigraphically constrained fossil calibrations to infer an evolutionary timescale for Serpentes. We reveal a potential diversification among crown snakes associated with the K-Pg mass extinction, led by the successful colonisation of Asia by the major extant clade Afrophidia. Vertebral morphometrics suggest increasing morphological specialisation among marine snakes through the Paleogene. The dispersal patterns of snakes following the K-Pg underscore the importance of this mass extinction event in shaping Earth's extant vertebrate faunas.

Holosteans contextualize the role of the teleost genome duplication in promoting the rise of evolutionary novelties in the ray-finned fish innate immune system

Over 99% of ray-finned fishes (Actinopterygii) are teleosts, a clade that comprises half of all living vertebrate species that have diversified across virtually all fresh and saltwater ecosystems. This ecological breadth raises the question of how the immunogenetic diversity required to persist under heterogeneous pathogen pressures evolved. The teleost genome duplication (TGD) has been hypothesized as the evolutionary event that provided the substrate for rapid genomic evolution and innovation. However, studies of putative teleost-specific innate immune receptors have been largely limited to comparisons either among teleosts or between teleosts and distantly related vertebrate clades such as tetrapods. Here we describe and characterize the receptor diversity of two clustered innate immune gene families in the teleost sister lineage: Holostei (bowfin and gars). Using genomic and transcriptomic data, we provide a detailed investigation of the phylogenetic history and conserved synteny of gene clusters encoding diverse immunoglobulin domain-containing proteins (DICPs) and novel immune-type receptors (NITRs). These data demonstrate an ancient linkage of DICPs to the major histocompatibility complex (MHC) and reveal an evolutionary origin of NITR variable-joining (VJ) exons that predate the TGD by at least 50 million years. Further characterizing the receptor diversity of Holostean DICPs and NITRs illuminates a sequence diversity that rivals the diversity of these innate immune receptor families in many teleosts. Taken together, our findings provide important historical context for the evolution of these gene families that challenge prevailing expectations concerning the consequences of the TGD during actinopterygiian evolution.

Phylogeny and sex chromosome evolution of palaeognathae

Many paleognaths (ratites and tinamous) have a pair of homomorphic ZW sex chromosomes in contrast to the highly differentiated sex chromosomes of most other birds. To understand the evolutionary causes for the different tempos of sex chromosome evolution, we produced female genomes of 12 paleognathous species and reconstructed the phylogeny and the evolutionary history of paleognathous sex chromosomes. We uncovered that Palaeognathae sex chromosomes had undergone stepwise recombination suppression and formed a pattern of "evolutionary strata". Nine of the 15 studied species' sex chromosomes have maintained homologous recombination in their long pseudoautosomal regions extending more than half of the entire chromosome length. We found that in older strata, the W chromosome suffered more serious functional gene loss. Their homologous Z-linked regions, compared with other genomic regions, have produced an excess of species-specific autosomal duplicated genes that evolved female-specific expression, in contrast to their broadly expressed progenitors. We speculate the "defeminization" of Z chromosome with underrepresentation of female-biased genes and slow divergence of sex chromosomes of paleognaths might be related to their distinctive mode of sexual selection targeting females that evolved in their common ancestors.

Novel Integrative Modeling of Molecules and Morphology across Evolutionary Timescales

Evolutionary models account for either population- or species-level processes but usually not both. We introduce a new model, the FBD-MSC, which makes it possible for the first time to integrate both the genealogical and fossilization phenomena, by means of the multispecies coalescent (MSC) and the fossilized birth–death (FBD) processes. Using this model, we reconstruct the phylogeny representing all extant and many fossil Caninae, recovering both the relative and absolute time of speciation events. We quantify known inaccuracy issues with divergence time estimates using the popular strategy of concatenating molecular alignments and show that the FBD-MSC solves them. Our new integrative method and empirical results advance the paradigm and practice of probabilistic total evidence analyses in evolutionary biology.[Caninae; fossilized birth–death; molecular clock; multispecies coalescent; phylogenetics; species trees.]

Bird neurocranial and body mass evolution across the end-Cretaceous mass extinction: The avian brain shape left other dinosaurs behind

Birds today are the most diverse clade of terrestrial vertebrates, and understanding why extant birds (Aves) alone among dinosaurs survived the Cretaceous-Paleogene mass extinction is crucial to reconstructing the history of life. Hypotheses proposed to explain this pattern demand identification of traits unique to Aves. However, this identification is complicated by a lack of data from non-avian birds. Here, we interrogate survivorship hypotheses using data from a new, nearly complete skull of Late Cretaceous (~70 million years) bird Ichthyornis and reassess shifts in bird body size across the Cretaceous-Paleogene boundary. Ichthyornis exhibited a wulst and segmented palate, previously proposed to have arisen within extant birds. The origin of Aves is marked by larger, reshaped brains indicating selection for relatively large telencephala and eyes but not by uniquely small body size. Sensory system differences, potentially linked to these shifts, may help explain avian survivorship relative to other dinosaurs.

Molecular Evolution of Ecological Specialisation: Genomic Insights from the Diversification of Murine Rodents

Adaptive radiations are characterized by the diversification and ecological differentiation of species, and replicated cases of this process provide natural experiments for understanding the repeatability and pace of molecular evolution. During adaptive radiation, genes related to ecological specialization may be subject to recurrent positive directional selection. However, it is not clear to what extent patterns of lineage-specific ecological specialization (including phenotypic convergence) are correlated with shared signatures of molecular evolution. To test this, we sequenced whole exomes from a phylogenetically dispersed sample of 38 murine rodent species, a group characterized by multiple, nested adaptive radiations comprising extensive ecological and phenotypic diversity. We found that genes associated with immunity, reproduction, diet, digestion, and taste have been subject to pervasive positive selection during the diversification of murine rodents. We also found a significant correlation between genome-wide positive selection and dietary specialization, with a higher proportion of positively selected codon sites in derived dietary forms (i.e., carnivores and herbivores) than in ancestral forms (i.e., omnivores). Despite striking convergent evolution of skull morphology and dentition in two distantly related worm-eating specialists, we did not detect more genes with shared signatures of positive or relaxed selection than in a nonconvergent species comparison. Although a small number of the genes we detected can be incidentally linked to craniofacial morphology or diet, protein-coding regions are unlikely to be the primary genetic basis of this complex convergent phenotype. Our results suggest a link between positive selection and derived ecological phenotypes, and highlight specific genes and general functional categories that may have played an integral role in the extensive and rapid diversification of murine rodents.

The Making of Calibration Sausage Exemplified by Recalibrating the Transcriptomic Timetree of Jawed Vertebrates

Molecular divergence dating has the potential to overcome the incompleteness of the fossil record in inferring when cladogenetic events (splits, divergences) happened, but needs to be calibrated by the fossil record. Ideally but unrealistically, this would require practitioners to be specialists in molecular evolution, in the phylogeny and the fossil record of all sampled taxa, and in the chronostratigraphy of the sites the fossils were found in. Paleontologists have therefore tried to help by publishing compendia of recommended calibrations, and molecular biologists unfamiliar with the fossil record have made heavy use of such works (in addition to using scattered primary sources and copying from each other). Using a recent example of a large node-dated timetree inferred from molecular data, I reevaluate all 30 calibrations in detail, present the current state of knowledge on them with its various uncertainties, rerun the dating analysis, and conclude that calibration dates cannot be taken from published compendia or other secondary or tertiary sources without risking strong distortions to the results, because all such sources become outdated faster than they are published: 50 of the (primary) sources I cite to constrain calibrations were published in 2019, half of the total of 280 after mid-2016, and 90% after mid-2005. It follows that the present work cannot serve as such a compendium either; in the slightly longer term, it can only highlight known and overlooked problems. Future authors will need to solve each of these problems anew through a thorough search of the primary paleobiological and chronostratigraphic literature on each calibration date every time they infer a new timetree, and that literature is not optimized for that task, but largely has other objectives.

The Implications of Lineage-Specific Rates for Divergence Time Estimation

Rate variation adds considerable complexity to divergence time estimation in molecular phylogenies. Here, we evaluate the impact of lineage-specific rates—which we define as among-branch-rate-variation that acts consistently across the entire genome. We compare its impact to residual rates—defined as among-branch-rate-variation that shows a different pattern of rate variation at each sampled locus, and gene-specific rates—defined as variation in the average rate across all branches at each sampled locus. We show that lineage-specific rates lead to erroneous divergence time estimates, regardless of how many loci are sampled. Further, we show that stronger lineage-specific rates lead to increasing error. This contrasts to residual rates and gene-specific rates, where sampling more loci significantly reduces error. If divergence times are inferred in a Bayesian framework, we highlight that error caused by lineage-specific rates significantly reduces the probability that the 95% highest posterior density includes the correct value, and leads to sensitivity to the prior. Use of a more complex rate prior—which has recently been proposed to model rate variation more accurately—does not affect these conclusions. Finally, we show that the scale of lineage-specific rates used in our simulation experiments is comparable to that of an empirical data set for the angiosperm genus Ipomoea. Taken together, our findings demonstrate that lineage-specific rates cause error in divergence time estimates, and that this error is not overcome by analyzing genomic scale multilocus data sets. [Divergence time estimation; error; rate variation.]

Uncertainty in Divergence Time Estimation

Understanding and representing uncertainty is crucial in academic research because it enables studies to build on the conclusions of previous studies, leading to robust advances in a particular field. Here, we evaluate the nature of uncertainty and the manner by which it is represented in divergence time estimation, a field that is fundamental to many aspects of macroevolutionary research, and where there is evidence that uncertainty has been seriously underestimated. We address this issue in the context of methods used in divergence time estimation, and with respect to the manner by which time-calibrated phylogenies are interpreted. With respect to methods, we discuss how the assumptions underlying different methods may not adequately reflect uncertainty about molecular evolution, the fossil record, or diversification rates. Therefore, divergence time estimates may not adequately reflect uncertainty and may be directly contradicted by subsequent findings. For the interpretation of time-calibrated phylogenies, we discuss how the use of time-calibrated phylogenies for reconstructing general evolutionary timescales leads to inferences about macroevolution that are highly sensitive to methodological limitations in how uncertainty is accounted for. By contrast, we discuss how the use of time-calibrated phylogenies to test specific hypotheses leads to inferences about macroevolution that are less sensitive to methodological limitations. Given that many biologists wish to use time-calibrated phylogenies to reconstruct general evolutionary timescales, we conclude that the development of methods of divergence time estimation that adequately account for uncertainty is necessary. [Divergence time estimation; macroevolution; uncertainty.].

The Origin of the Legumes is a Complex Paleopolyploid Phylogenomic Tangle Closely Associated with the Cretaceous-Paleogene (K-Pg) Mass Extinction Event

The consequences of the Cretaceous–Paleogene (K–Pg) boundary (KPB) mass extinction for the evolution of plant diversity remain poorly understood, even though evolutionary turnover of plant lineages at the KPB is central to understanding assembly of the Cenozoic biota. The apparent concentration of whole genome duplication (WGD) events around the KPB may have played a role in survival and subsequent diversification of plant lineages. To gain new insights into the origins of Cenozoic biodiversity, we examine the origin and early evolution of the globally diverse legume family (Leguminosae or Fabaceae). Legumes are ecologically (co-)dominant across many vegetation types, and the fossil record suggests that they rose to such prominence after the KPB in parallel with several well-studied animal clades including Placentalia and Neoaves. Furthermore, multiple WGD events are hypothesized to have occurred early in legume evolution. Using a recently inferred phylogenomic framework, we investigate the placement of WGDs during early legume evolution using gene tree reconciliation methods, gene count data and phylogenetic supernetwork reconstruction. Using 20 fossil calibrations we estimate a revised timeline of legume evolution based on 36 nuclear genes selected as informative and evolving in an approximately clock-like fashion. To establish the timing of WGDs we also date duplication nodes in gene trees. Results suggest either a pan-legume WGD event on the stem lineage of the family, or an allopolyploid event involving (some of) the earliest lineages within the crown group, with additional nested WGDs subtending subfamilies Papilionoideae and Detarioideae. Gene tree reconciliation methods that do not account for allopolyploidy may be misleading in inferring an earlier WGD event at the time of divergence of the two parental lineages of the polyploid, suggesting that the allopolyploid scenario is more likely. We show that the crown age of the legumes dates to the Maastrichtian or early Paleocene and that, apart from the Detarioideae WGD, paleopolyploidy occurred close to the KPB. We conclude that the early evolution of the legumes followed a complex history, in which multiple auto- and/or allopolyploidy events coincided with rapid diversification and in association with the mass extinction event at the KPB, ultimately underpinning the evolutionary success of the Leguminosae in the Cenozoic. [Allopolyploidy; Cretaceous–Paleogene (K–Pg) boundary; Fabaceae, Leguminosae; paleopolyploidy; phylogenomics; whole genome duplication events]

Deep-Time Demographic Inference Suggests Ecological Release as Driver of Neoavian Adaptive Radiation

Assessing the applicability of theory to major adaptive radiations in deep time represents an extremely difficult problem in evolutionary biology. Neoaves, which includes 95% of living birds, is believed to have undergone a period of rapid diversification roughly coincident with the Cretaceous–Paleogene (K-Pg) boundary. We investigate whether basal neoavian lineages experienced an ecological release in response to ecological opportunity, as evidenced by density compensation. We estimated effective population sizes (Ne) of basal neoavian lineages by combining coalescent branch lengths (CBLs) and the numbers of generations between successive divergences. We used a modified version of Accurate Species TRee Algorithm (ASTRAL) to estimate CBLs directly from insertion–deletion (indel) data, as well as from gene trees using DNA sequence and/or indel data. We found that some divergences near the K-Pg boundary involved unexpectedly high gene tree discordance relative to the estimated number of generations between speciation events. The simplest explanation for this result is an increase in Ne, despite the caveats discussed herein. It appears that at least some early neoavian lineages, similar to the ancestor of the clade comprising doves, mesites, and sandgrouse, experienced ecological release near the time of the K-Pg mass extinction.

Noncoding regions underpin avian bill shape diversification at macroevolutionary scales

Recent progress has been made in identifying genomic regions implicated in trait evolution on a microevolutionary scale in many species, but whether these are relevant over macroevolutionary time remains unclear. Here, we directly address this fundamental question using bird beak shape, a key evolutionary innovation linked to patterns of resource use, divergence, and speciation, as a model trait. We integrate class-wide geometric-morphometric analyses with evolutionary sequence analyses of 10,322 protein-coding genes as well as 229,001 genomic regions spanning 72 species. We identify 1434 protein-coding genes and 39,806 noncoding regions for which molecular rates were significantly related to rates of bill shape evolution. We show that homologs of the identified protein-coding genes as well as genes in close proximity to the identified noncoding regions are involved in craniofacial embryo development in mammals. They are associated with embryonic stem cell pathways, including BMP and Wnt signaling, both of which have repeatedly been implicated in the morphological development of avian beaks. This suggests that identifying genotype-phenotype association on a genome-wide scale over macroevolutionary time is feasible. Although the coding and noncoding gene sets are associated with similar pathways, the actual genes are highly distinct, with significantly reduced overlap between them and bill-related phenotype associations specific to noncoding loci. Evidence for signatures of recent diversifying selection on our identified noncoding loci in Darwin finch populations further suggests that regulatory rather than coding changes are major drivers of morphological diversification over macroevolutionary times.

Late Cretaceous neornithine from Europe illuminates the origins of crown birds

Our understanding of the earliest stages of crown bird evolution is hindered by an exceedingly sparse avian fossil record from the Mesozoic era. The most ancient phylogenetic divergences among crown birds are known to have occurred in the Cretaceous period^1-3, but stem-lineage representatives of the deepest subclades of crown birds-Palaeognathae (ostriches and kin), Galloanserae (landfowl and waterfowl) and Neoaves (all other extant birds)-are unknown from the Mesozoic era. As a result, key questions related to the ecology^4,5, biogeography^3,6,7 and divergence times^1,8-10 of ancestral crown birds remain unanswered. Here we report a new Mesozoic fossil that occupies a position close to the last common ancestor of Galloanserae and fills a key phylogenetic gap in the early evolutionary history of crown birds^10,11. Asteriornis maastrichtensis, gen. et sp. nov., from the Maastrichtian age of Belgium (66.8-66.7 million years ago), is represented by a nearly complete, three-dimensionally preserved skull and associated postcranial elements. The fossil represents one of the only well-supported crown birds from the Mesozoic era¹², and is the first Mesozoic crown bird with well-represented cranial remains. Asteriornis maastrichtensis exhibits a previously undocumented combination of galliform (landfowl)-like and anseriform (waterfowl)-like features, and its presence alongside a previously reported Ichthyornis-like taxon from the same locality¹³ provides direct evidence of the co-occurrence of crown birds and avialan stem birds. Its occurrence in the Northern Hemisphere challenges biogeographical hypotheses of a Gondwanan origin of crown birds³, and its relatively small size and possible littoral ecology may corroborate proposed ecological filters^4,5,9 that influenced the persistence of crown birds through the end-Cretaceous mass extinction.

The dynamics of stem and crown groups

The fossil record of the origins of major groups such as animals and birds has generated considerable controversy, especially when it conflicts with timings based on molecular clock estimates. Here, we model the diversity of "stem" (basal) and "crown" (modern) members of groups using a "birth-death model," the results of which qualitatively match many large-scale patterns seen in the fossil record. Typically, the stem group diversifies rapidly until the crown group emerges, at which point its diversity collapses, followed shortly by its extinction. Mass extinctions can disturb this pattern and create long stem groups such as the dinosaurs. Crown groups are unlikely to emerge either cryptically or just before mass extinctions, in contradiction to popular hypotheses such as the "phylogenetic fuse". The patterns revealed provide an essential context for framing ecological and evolutionary explanations for how major groups originate, and strengthen our confidence in the reliability of the fossil record.

Resolving Phylogenetic Relationships within Passeriformes Based on Mitochondrial Genes and Inferring the Evolution of Their Mitogenomes in Terms of Duplications

Mitochondrial genes are placed on one molecule, which implies that they should carry consistent phylogenetic information. Following this advantage, we present a well-supported phylogeny based on mitochondrial genomes from almost 300 representatives of Passeriformes, the most numerous and differentiated Aves order. The analyses resolved the phylogenetic position of paraphyletic Basal and Transitional Oscines. Passerida occurred divided into two groups, one containing Paroidea and Sylvioidea, whereas the other, Passeroidea and Muscicapoidea. Analyses of mitogenomes showed four types of rearrangements including a duplicated control region (CR) with adjacent genes. Mapping the presence and absence of duplications onto the phylogenetic tree revealed that the duplication was the ancestral state for passerines and was maintained in early diverged lineages. Next, the duplication could be lost and occurred independently at least four times according to the most parsimonious scenario. In some lineages, two CR copies have been inherited from an ancient duplication and highly diverged, whereas in others, the second copy became similar to the first one due to concerted evolution. The second CR copies accumulated over twice as many substitutions as the first ones. However, the second CRs were not completely eliminated and were retained for a long time, which suggests that both regions can fulfill an important role in mitogenomes. Phylogenetic analyses based on CR sequences subjected to the complex evolution can produce tree topologies inconsistent with real evolutionary relationships between species. Passerines with two CRs showed a higher metabolic rate in relation to their body mass.