The phylogeny and divergence timing of the Neoavian radiation remain controversial despite recent progress. We analyzed the genomes of 124 species across all Neoavian orders, using data from 25,460 loci spanning four DNA classes, including 5,756 coding sequences, 12,449 conserved nonexonic elements, 4,871 introns, and 2,384 intergenic segments. We conducted a comprehensive sensitivity analysis to account for the heterogeneity across different DNA classes, leading to an optimal tree of Neoaves with high resolution. This phylogeny features a novel Neoavian dichotomy comprising two monophyletic clades: a previously recognized Telluraves (land birds) and a newly circumscribed Aquaterraves (waterbirds and relatives). Molecular dating analyses with 20 fossil calibrations indicate that the diversification of modern birds began in the Late Cretaceous and underwent a constant and steady radiation across the KPg boundary, concurrent with the rise of angiosperms as well as other major Cenozoic animal groups including placental and multituberculate mammals. The KPg catastrophe had a limited impact on avian evolution compared to the Paleocene-Eocene Thermal Maximum, which triggered a rapid diversification of seabirds. Our findings suggest that the evolution of modern birds followed a slow process of gradualism rather than a rapid process of punctuated equilibrium, with limited interruption by the KPg catastrophe. This study places bird evolution into a new context within vertebrates, with ramifications for the evolution of the Earth's biota.

Melanin pigments play a critical role in physiological processes and shaping animal behaviour. Fossil melanin is a unique resource for understanding the functional evolution of melanin but the impact of fossilisation on molecular signatures for eumelanin and, especially, phaeomelanin is not fully understood. Here we present a model for the chemical taphonomy of fossil eumelanin and phaeomelanin based on thermal maturation experiments using feathers from extant birds. Our results reveal which molecular signatures are authentic signals for thermally matured eumelanin and phaeomelanin, which signatures are artefacts derived from the maturation of non-melanin molecules, and how these chemical data are impacted by sample preparation. Our model correctly predicts the molecular composition of eumelanins in diverse vertebrate fossils from the Miocene and Cretaceous and, critically, identifies direct molecular evidence for phaeomelanin in these fossils. This taphonomic framework adds to the geochemical toolbox that underpins reconstructions of melanin evolution and of melanin-based coloration in fossil vertebrates.

The ancestor of cetaceans underwent a macroevolutionary transition from land to water early in the Eocene Period >50 million years ago. However, little is known about how diverse retroviruses evolved during this shift from terrestrial to aquatic environments. Did retroviruses transition into water accompanying their hosts? Did retroviruses infect cetaceans through cross-species transmission after cetaceans invaded the aquatic environments? Endogenous retroviruses (ERVs) provide important molecular fossils for tracing the evolution of retroviruses during this macroevolutionary transition. Here, we use a phylogenomic approach to study the origin and evolution of ERVs in cetaceans. We identify a total of 8,724 ERVs within the genomes of 25 cetaceans, and phylogenetic analyses suggest these ERVs cluster into 315 independent lineages, each of which represents one or more independent endogenization events. We find that cetacean ERVs originated through two possible routes. 298 ERV lineages may derive from retrovirus endogenization that occurred before or during the transition from land to water of cetaceans, and most of these cetacean ERVs were reaching evolutionary dead-ends. 17 ERV lineages are likely to arise from independent retrovirus endogenization events that occurred after the split of mysticetes and odontocetes, indicating that diverse retroviruses infected cetaceans through cross-species transmission from non-cetacean mammals after the transition to aquatic life of cetaceans. Both integration time and synteny analyses support the recent or ongoing activity of multiple retroviral lineages in cetaceans, some of which proliferated into hundreds of copies within the host genomes. Although ERVs only recorded a proportion of past retroviral infections, our findings illuminate the complex evolution of retroviruses during one of the most marked macroevolutionary transitions in vertebrate history.

Medullary bone is an ephemeral type of bone tissue, today found only in sexually mature female birds, that provides a calcium reservoir for eggshell formation. The presence of medullary bone-like tissues in extant birds, pterosaurs, and dinosaurs distantly related to birds shows that caution must be exercised before concluding that fossils bear medullary bone. Here we describe a new specimen of pengornithid enantiornithine from the Lower Cretaceous Jiufotang Formation. Consisting of an isolated left hindlimb, the three-dimensional preservation contrasts with the crushed preservation characteristic of most Jehol specimens. Histological examinations suggest this resulted from the presence of a thick layer of highly vascular bone spanning the medullary cavities of the femur and tibiotarsus, consistent with expectations for medullary bone in extant birds. Micro-computed tomographic scans reveal small amounts of the same tissue extending into the pedal phalanges. We consider the tissue to be homologous to the medullary bone of Neornithines.

Caviidae is a diverse group of caviomorph rodents that is broadly distributed in South America and is divided into three highly divergent extant lineages: Caviinae (cavies), Dolichotinae (maras), and Hydrochoerinae (capybaras). The fossil record of Caviidae is only abundant and diverse since the late Miocene. Caviids belongs to Cavioidea sensu stricto (Cavioidea s.s.) that also includes a diverse assemblage of extinct taxa recorded from the late Oligocene to the middle Miocene of South America ("eocardiids").

Although stone tools generally co-occur with early members of the genus Homo, they are rarely found in direct association with hominins. We report that both Acheulian and Oldowan artifacts and Homo erectus crania were found in close association at 1.26 million years (Ma) ago at Busidima North (BSN12), and ca. 1.6 to 1.5 Ma ago at Dana Aoule North (DAN5) archaeological sites at Gona, Afar, Ethiopia. The BSN12 partial cranium is robust and large, while the DAN5 cranium is smaller and more gracile, suggesting that H. erectus was probably a sexually dimorphic species. The evidence from Gona shows behavioral diversity and flexibility with a lengthy and concurrent use of both stone technologies by H. erectus, confounding a simple "single species/single technology" view of early Homo.

While pelagornithid or 'bony-toothed' bird fossils representing multiple species are known from Antarctica, a new dentary fragment of a pelagornithid bird from the middle Eocene Submeseta Formation on Seymour Island, Antarctica represents a species with a body size on par with the largest known species in the clade. Measurements from the partial 'toothed' dentary point to a giant body size for the species, although the spacing among the pseudoteeth differs from that published for other pelagornithids. The discrepancy might suggest that previous techniques are not adequate for examination of incomplete material or that another factor such as phylogeny might impact size estimates and comparisons. Combined with a revised stratigraphic position in the early Eocene La Meseta Formation on Seymour Island for the largest pelagornithid tarsometatarsus known, these Antarctic fossils demonstrate the early evolution of giant body size in the clade (by ~ 50 Ma), and they likely represent not only the largest flying birds of the Eocene but also some of the largest volant birds that ever lived (with an estimated 5-6 m wingspan). Furthermore, the distribution of giant-sized pelagornithid fossils across more than 10 million years of Antarctic geological deposits points to a prolonged survival of giant-sized pelagornithids within the southern seas, and their success as a pelagic predatory component of marine and coastal ecosystems alongside early penguins.

The interplay of different virus species in a host cell after infection can affect the adaptation of each virus. Endogenous viral elements, such as endogenous pararetroviruses (PRVs), have arisen from vertical inheritance of viral sequences integrated into host germline genomes. As viral genomic fossils, these sequences can thus serve as valuable paleogenomic data to study the long-term evolutionary dynamics of virus-virus interactions, but they have rarely been applied for this purpose. All extant PRVs have been considered autonomous species in their parasitic life cycle in host cells. Here, we provide evidence for multiple non-autonomous PRV species with structural defects in viral activity that have frequently infected ancient grass hosts and adapted through interplay between viruses. Our paleogenomic analyses using endogenous PRVs in grass genomes revealed that these non-autonomous PRV species have participated in interplay with autonomous PRVs in a possible commensal partnership, or, alternatively, with one another in a possible mutualistic partnership. These partnerships, which have been established by the sharing of noncoding regulatory sequences (NRSs) in intergenic regions between two partner viruses, have been further maintained and altered by the sequence homogenization of NRSs between partners. Strikingly, we found that frequent region-specific recombination, rather than mutation selection, is the main causative mechanism of NRS homogenization. Our results, obtained from ancient DNA records of viruses, suggest that adaptation of PRVs has occurred by concerted evolution of NRSs between different virus species in the same host. Our findings further imply that evaluation of within-host NRS interactions within and between populations of viral pathogens may be important.

The Talbragar Fish Bed is one of Australia's most important Jurassic deposits for freshwater fishes, land plants and aquatic and terrestrial insects. The site has yielded many well preserved fossils, which has led to the formal description of numerous new species and higher taxa. The excellent preservation of many fossils has allowed detailed anatomical studies, e.g. of the early teleost fish Cavenderichthys talbragarensis (Woodward, 1895). Here we report on the fluorescent characteristics and mineral composition of a range of Talbragar fossils. Most specimens fluoresce under ultraviolet, blue and green light. Elemental and mineralogical analyses revealed that the Talbragar fossils consist predominantly of quartz (SiO2), a mineral that is likely to account for the observed fluorescence, with trace kaolinite (Al2Si2O5(OH)4) in some of the fish fossils. Rock matrices are predominantly composed of quartz and goethite (FeO(OH)). Closer inspection of a plant leaf (Pentoxylon australicum White, 1981) establishes fluorescence as a useful tool for the visualisation of anatomical details that are difficult to see under normal light conditions.

Integration of diverse data (molecules, fossils) provides the most robust test of the phylogeny of cetaceans. Positioning key fossils is critical for reconstructing the character change from life on land to life in the water.

Humic substances are breakdown products of decaying organic matter that co-extract with proteins from fossils. These substances are difficult to separate from proteins in solution and interfere with analyses of fossil proteomes. We introduce a method combining multiple recent advances in extraction protocols to both concentrate proteins from fossil specimens with high humic content and remove humics, producing clean samples easily analysed by mass spectrometry (MS). This method includes: (i) a non-demineralizing extraction buffer that eliminates protein loss during the demineralization step in routine methods; (ii) filter-aided sample preparation (FASP) of peptides, which concentrates and digests extracts in one filter, allowing the separation of large humics after digestion; (iii) centrifugal stage tipping, which further clarifies and concentrates samples in a uniform process performed simultaneously on multiple samples. We apply this method to a moa fossil (approx. 800-1000 years) dark with humic content, generating colourless samples and enabling the detection of more proteins with greater sequence coverage than previous MS analyses on this same specimen. This workflow allows analyses of low-abundance proteins in fossils containing humics and thus may widen the range of extinct organisms and regions of their proteomes we can explore with MS.

The Rhinopomatidae, traditionally considered to be one of the most ancient chiropteran clades, remains one of the least known groups of Rhinolophoidea. No relevant fossil record is available for this family. Whereas there have been extensive radiations in related families Rhinolophidae and Hipposideridae, there are only a few species in the Rhinopomatidae and their phylogenetic relationship and status are not fully understood.

In disciplines such as macroevolution that are not amenable to experimentation, scientists usually rely on current observations to test hypotheses about historical events, assuming that "the present is the key to the past." Biogeographers, for example, used this assumption to reconstruct ancestral ranges from the distribution of extant species. Yet, under scenarios of high extinction rates, the biodiversity we observe today might not be representative of the historical diversity and this could result in incorrect biogeographic reconstructions. Here, we introduce a new approach to incorporate into biogeographic inference the temporal, spatial, and environmental information provided by the fossil record, as a direct evidence of the extinct biodiversity fraction. First, inferences of ancestral ranges for those nodes in the phylogeny calibrated with the fossil record are constrained to include the geographic distribution of the fossil. Second, we use fossil distribution and past climate data to reconstruct the climatic preferences and potential distribution of ancestral lineages over time, and use this information to build a biogeographic model that takes into account "ecological connectivity" through time. To show the power of this approach, we reconstruct the biogeographic history of the large angiosperm genus Hypericum, which has a fossil record extending back to the Early Cenozoic. Unlike previous reconstructions based on extant species distributions, our results reveal that Hypericum stem lineages were already distributed in the Holarctic before diversification of its crown-group, and that the geographic distribution of the genus has been relatively stable throughout the climatic oscillations of the Cenozoic. Geographical movement was mediated by the existence of climatic corridors, like Beringia, whereas the equatorial tropical belt acted as a climatic barrier, preventing Hypericum lineages to reach the southern temperate regions. Our study shows that an integrative approach to historical biogeography-that combines sources of evidence as diverse as paleontology, ecology, and phylogenetics-could help us obtain more accurate reconstructions of ancient evolutionary history. It also reveals the confounding effect different rates of extinction across regions have in biogeography, sometimes leading to ancestral areas being erroneously inferred as recent colonization events.

The 28S rRNA genes of several metazoans are interrupted by site-specific targeting non-LTR retrotransposons, such as R2. R2 elements have been deeply analyzed but aspects of their retrotransposition mechanism and the origin of the wide diversity observed are still debated. We characterized six new R2 lineages in four tadpole shrimp species (Notostraca), samples deriving from a parthenogenetic population of Triops cancriformis (R2Tc_it) and from bisexual Lepidurus populations of L. lubbocki (R2Ll), L. couesii (R2LcA, R2LcB, R2LcC) and L. arcticus (R2La). All elements fit the canonical R2 structure but R2Ll which turned out to be a chimera with an additional ORF originating from another R2. Consistently with data on LINEs, R2Ll could be the result of recombination due to reverse transcriptase template jump. The analysis of 28S/R2 5' end junctions further suggests aberrant homologous recombination, as observed in RNA viruses.

We have developed a protocol for reconstructing 3D models of the skulls of extinct species of small mammals. For the first time, the reconstruction uses fragments of fossils from a mixture of different specimens and from related extant species. We use free software and commercial computers to make the process reproducible and usable for the scientific community. We present a semi-quantitative protocol to face the problem of making 3D reconstructions of fossil species that are incomplete in the fossil record and/or represented by a mixture of different individuals, as usually occurs with small vertebrates. Therefore this approach is useful when no complete skull is available. The protocol combines the use of microCT scan technology with a subsequent computer treatment using different software tools for 3D reconstruction from microCT and 3D design and printing (e.g. Fiji, SPIERS, Meshlab, Meshmixer) in a defined order. This kind of free and relatively simple software, plus the detailed description, makes this protocol practicable for researchers who do not necessarily have great deal of experience in working with 3D. As an example, we have performed virtual reconstructions of the skulls of two species of insectivore small mammals (Eulipotyphla): Beremendia fissidens and Dolinasorex glyphodon. The resulting skulls, plus models of the extant shrews Blarina brevicauda, Neomys fodiens, Crocidura russula and Sorex coronatus, make it possible to compare characteristics that can only be observed by means of microCT 3D reconstructions, and given the characteristics of the material, using this protocol. Among the characters we can compare are the position of the mandibles, the spatial relations among all the teeth, the shape of the snout and, in general, all parameters related with the anatomy of the rostrum. Moreover, these reconstructions can be used in different types of context: for anatomical purposes, especially to see internal features or characteristics at whole-skull scale, for bioengineering, animation, or other techniques that need a digital model.

Dental material attributed to Palaeoloxodon huaihoensis from the Middle to Late Pleistocene were recovered over decades from the Penghu Channel during commercial fisheries activities. The National Museum of Nature Science (NMNS) has a collection of such dental material, which differs in size and morphology and likely represents ontogenetic variation and growth trajectory of various age groups of P. huaihoensis. However, little is known regarding age determination. By using length of dental material, enamel thickness (ET), and plate counts, we established the method to distinguish the age of the species, which is directly derived from the extant African forest elephant Loxodonta africana. When measuring signs of allometric growth, we found that in both the upper and lower jaws, tooth width was correlated negatively with lamellar frequency but positively with ET. In the same age group, the number of lamellae was higher in P. huaihoensis than in L. africana. The reconstructed age distribution indicated no difference in the upper or lower jaw. Notably, within our sample, P. huaihoensis is skewed towards adult and older individuals with median age between 33-34.5 years and differed significantly from that of Mammuthus primigenius in the European Kraków Spadzista site. This age distribution pattern is speculated to be related to the harsh environmental conditions and intense intraspecific competition among P. huaihoensis during the last ice age.

Pan-trionychids or softshell turtles are a highly specialized and widespread extant group of aquatic taxa with an evolutionary history that goes back to the Early Cretaceous. The earliest pan-trionychids had already fully developed the "classic" softshell turtle morphology and it has been impossible to resolve whether they are stem members of the family or are within the crown. This has hindered our understanding of the evolution of the two basic body plans of crown-trionychids. Thus it remains unclear whether the more heavily ossified shell of the cyclanorbines or the highly reduced trionychine morphotype is the ancestral condition for softshell turtles. A new pan-trionychid from the Early Cretaceous of Zhejiang, China, Perochelys hengshanensis sp. nov., allows a revision of softshell-turtle phylogeny. Equal character weighting resulted in a topology that is fundamentally inconsistent with molecular divergence date estimates of deeply nested extant species. In contrast, implied weighting retrieved Lower Cretaceous Perochelys spp. and Petrochelys kyrgyzensis as stem trionychids, which is fully consistent with their basal stratigraphic occurrence and an Aptian-Santonian molecular age estimate for crown-trionychids. These results indicate that the primitive morphology for soft-shell turtles is a poorly ossified shell like that of crown-trionychines and that shell re-ossification in cyclanorbines (including re-acquisition of peripheral elements) is secondary.

Catostomidae ("suckers") is a diverse (76 species) and broadly distributed family of Holarctic freshwater fishes with a rich fossil record and a considerable number (∼35%) of threatened and imperiled species. We integrate DNA sequences (three mitochondrial genes, three nuclear genes), morphological data, and fossil information to infer sucker phylogenetic relationships and divergence times using Bayesian "total-evidence" methods, and then test hypotheses about the temporal diversification of the group. Our analyses resolved many nodes within subfamilies and clarified Catostominae relationships to be of the form ((Thoburniini, Moxostomatini), (Erimyzonini, Catostomini)). Patterns of subfamily relationships were incongruent, but mainly supported two placements of the Myxocyprininae; distinguishing these using Bayes factors lent strongest support to a model with Myxocyprininae sister to all remaining sucker lineages. We improved our Bayesian total-evidence dating analysis by excluding problematic characters, using a clock-partitioning scheme identified by Bayesian model selection, and employing a fossilized birth-death tree prior accommodating morphological data and fossils. The resulting chronogram showed that suckers evolved since the Late Cretaceous-Eocene, and that the Catostomini and Moxostomatini clades have accumulated species diversity since the early to mid-Miocene. These results agree with the fossil record and confirm previous hypotheses about dates for the origins of Catostomide and catostomine diversification, but reject previous molecular hypotheses about the timing of divergence of ictiobines, and between Asian-North American lineages. Overall, our findings from a synthesis of multiple data types enhance understanding of the phylogenetic relationships, taxonomic classification, and temporal diversification of suckers, while also highlighting practical methods for improving Bayesian divergence dating models by coupling phylogenetic informativeness profiling with relaxed-clock partitioning.

The fossil record constitutes the primary source of information about the evolutionary history of extant and extinct groups, and many analyses of macroevolution rely on fossils that are accurately placed within phylogenies. To avoid misinterpretation of the fossil record, especially by non-palaeontologists, the proper assessment and communication of uncertainty in fossil placement is crucial. We here use Bayesian morphological phylogenetics to evaluate the classifications of fossil parasitoid wasps (Hymenoptera, Ichneumonidae) and introduce 'RoguePlots' to illustrate placement uncertainty on the phylogeny of extant taxa. Based on an extensive, newly constructed morphological matrix of 222 characters in 24 fossil and 103 extant taxa, we test three different aspects of models of morphological evolution. We find that a model that includes ordered characters, among-character rate variation, and a state-space restricted to observed states achieves the highest marginal likelihoods. The individual RoguePlots reveal large differences in confidence in the placement of the different fossils and allow some refinements to their classification: Polyhelictes bipolarus and Ichninsum appendicrassum are moved from an uncertain subfamily placement to Pimplinae, Plectiscidea lanhami is transferred to Allomacrus in Cylloceriinae (Allomacrus lanhami, comb. nov.), Lithotorus cressoni is moved from Diplazontinae to Orthocentrinae, and we note uncertainty in the generic placements of Rhyssella? vera and Xanthopimpla? messelensis. We discuss potential artefacts that might result in biased posterior probabilities in Bayesian morphological phylogenetic analyses, pertaining to character and taxon sampling, fossilization biases, and model misspecification. Finally, we suggest future directions both in ichneumonid palaeontology and in the way RoguePlots can improve both assessment and representation of placement uncertainty, both in fossils and other rogue taxa.

Old World porcupines (Family: Hystricidae) are the third-largest rodents and inhabit southern Europe, Asia, and most regions of Africa. They are a typical indicator of warm climate and their distribution is restricted to tropical and subtropical zones. In China, porcupines are widely distributed in southern areas of the Yangtze River. However, fossil remains have been identified in a few sites in northern China, among which Tianyuan Cave-near Zhoukoudian site-represents the latest known porcupine fossil record. So far, studies have focused mainly on porcupines' husbandry and domestication but little is known about their intrafamilial phylogenetic relationships and evolutionary history.