Molecules and fossils tell distinct yet complementary stories of mammal diversification

Imbalanced speciation pulses sustain the radiation of mammals

The evolutionary histories of major clades, including mammals, often comprise changes in their diversification dynamics, but how these changes occur remains debated. We combined comprehensive phylogenetic and fossil information in a new "birth-death diffusion" model that provides a detailed characterization of variation in diversification rates in mammals. We found an early rising and sustained diversification scenario, wherein speciation rates increased before and during the Cretaceous-Paleogene (K-Pg) boundary. The K-Pg mass extinction event filtered out more slowly speciating lineages and was followed by a subsequent slowing in speciation rates rather than rebounds. These dynamics arose from an imbalanced speciation process, with separate lineages giving rise to many, less speciation-prone descendants. Diversity seems to have been brought about by these isolated, fast-speciating lineages, rather than by a few punctuated innovations.

A large therian mammal from the Late Cretaceous of South America

Theria represent an extant clade that comprises placental and marsupial mammals. Here we report on the discovery of a new Late Cretaceous mammal from southern Patagonia, Patagomaia chainko gen. et sp. nov., represented by hindlimb and pelvic elements with unambiguous therian features. We estimate Patagomaia chainko attained a body mass of 14 kg, which is considerably greater than the 5 kg maximum body mass of coeval Laurasian therians. This new discovery demonstrates that Gondwanan therian mammals acquired large body size by the Late Cretaceous, preceding their Laurasian relatives, which remained small-bodied until the beginning of the Cenozoic. Patagomaia supports the view that the Southern Hemisphere was a cradle for the evolution of modern mammalian clades, alongside non-therian extinct groups such as meridiolestidans, gondwanatherians and monotremes.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

The relationship between geographic range size and rates of species diversification

Range size is a universal characteristic of every biological species, and is often assumed to affect diversification rate. There are strong theoretical arguments that large-ranged species should have higher rates of diversification. On the other hand, the observation that small-ranged species are often phylogenetically clustered might indicate high diversification of small-ranged species. This discrepancy between theory and the data may be caused by the fact that typical methods of data analysis do not account for range size changes during speciation. Here we use a cladogenetic state-dependent diversification model applied to mammals to show that range size changes during speciation are ubiquitous and small-ranged species indeed diversify generally slower, as theoretically expected. However, both range size and diversification are strongly influenced by idiosyncratic and spatially localized events, such as colonization of an archipelago or a mountain system, which often override the general pattern of range size evolution.

Towards understanding paleoclimate impacts on primate de novo genes

De novo genes are genes that emerge as new genes in some species, such as primate de novo genes that emerge in certain primate species. Over the past decade, a great deal of research has been conducted regarding their emergence, origins, functions, and various attributes in different species, some of which have involved estimating the ages of de novo genes. However, limited by the number of species available for whole-genome sequencing, relatively few studies have focused specifically on the emergence time of primate de novo genes. Among those, even fewer investigate the association between primate gene emergence with environmental factors, such as paleoclimate (ancient climate) conditions. This study investigates the relationship between paleoclimate and human gene emergence at primate species divergence. Based on 32 available primate genome sequences, this study has revealed possible associations between temperature changes and the emergence of de novo primate genes. Overall, findings in this study are that de novo genes tended to emerge in the recent 13 MY when the temperature continues cooling, which is consistent with past findings. Furthermore, in the context of an overall trend of cooling temperature, new primate genes were more likely to emerge during local warming periods, where the warm temperature more closely resembled the environmental condition that preceded the cooling trend. Results also indicate that both primate de novo genes and human cancer-associated genes have later origins in comparison to random human genes. Future studies can be in-depth on understanding human de novo gene emergence from an environmental perspective as well as understanding species divergence from a gene emergence perspective.

A timescale for placental mammal diversification based on Bayesian modeling of the fossil record

The timing of the placental mammal radiation has been the focus of debate over the efficacy of competing methods for establishing evolutionary timescales. Molecular clock analyses estimate that placental mammals originated before the Cretaceous-Paleogene (K-Pg) mass extinction, anywhere from the Late Cretaceous to the Jurassic. However, the absence of definitive fossils of placentals before the K-Pg boundary is compatible with a post-Cretaceous origin. Nevertheless, lineage divergence must occur before it can be manifest phenotypically in descendent lineages. This, combined with the non-uniformity of the rock and fossil records, requires the fossil record to be interpreted rather than read literally. To achieve this, we introduce an extended Bayesian Brownian bridge model that estimates the age of origination and, where applicable, extinction through a probabilistic interpretation of the fossil record. The model estimates the origination of placentals in the Late Cretaceous, with ordinal crown groups originating at or after the K-Pg boundary. The results reduce the plausible interval for placental mammal origination to the younger range of molecular clock estimates. Our findings support both the Long Fuse and Soft Explosive models of placental mammal diversification, indicating that the placentals originated shortly prior to the K-Pg mass extinction. The origination of many modern mammal lineages overlapped with and followed the K-Pg mass extinction.

Chemical signaling glands are unlinked to species diversification in lizards

Sexual selection has long been thought to increase species diversification. Sexually selected traits, such as sexual signals that contribute to reproductive isolation, were thought to promote diversification. However, studies exploring links between sexually selected traits and species diversification have thus far primarily focused on visual or acoustic signals. Many animals often employ chemical signals (i.e., pheromones) for sexual communications, but large-scale analyses on the role of chemical communications in driving species diversification have been missing. Here, for the first time, we investigate whether traits associated with chemical communications-the presence of follicular epidermal glands-promote diversification across 6,672 lizard species. In most analyses, we found no strong association between the presence of follicular epidermal glands and species diversification rates, either across all lizard species or at lower phylogenetic scales. Previous studies suggest that follicular gland secretions act as species recognition signals that prevent hybridization during speciation in lizards. However, we show that geographic range overlap was no different in sibling species pairs with and without follicular epidermal glands. Together, these results imply that either follicular epidermal glands do not primarily function in sexual communications or sexually selected traits in general (here chemical communication) have a limited effect on species diversification. In our additional analysis accounting for sex-specific differences in glands, we again found no detectable effect of follicular epidermal glands on species diversification rates. Thus, our study challenges the general role of sexually selected traits in broad-scale species diversification patterns.

A new species of the highly polytypic South American rodent Ctenomys increases the diversity of the magellanicus clade

The subterranean rodent Ctenomys is the most polytypic South American mammal genus and one of the most speciose and rapidly diversifying mammal genera in the world. Its systematics is unstable due to the underlying accelerated diversification processes that give rise to evolutionary lineages at different stages of differentiation and to remarkable morphological homogeneity even among long-differentiated species. As a result, species boundaries are often difficult to define. Diversity of this genus in the coastal area of central Argentina has been extensively studied, with two independent lineages currently recognized while a distinct third population had not been previously detected. Through a phylogenetic analysis based on combined morphological and molecular evidence, Bayesian estimates of divergence times, and morphometric and morphological assessments, we recognize this third population as an independently evolving lineage. The new species, Ctenomys pulcersp. nov., is here described for both the living fauna and the fossil record of the Pampean region of central Argentina. According to phylogenetic results, Ctenomys pulcersp. nov. belongs to the essentially Patagonian magellanicus clade, and would have diverged from its sister species, Ctenomys bidaui, during the middle Pleistocene (ca. 0.4 Ma). Its current distribution in the fixed and semifixed dunes of the coastal Pampean region is assumed to represent a relict of a wider and continuous distribution of potentially suitable environments during the late Pleistocene. Ctenomys pulcersp. nov. occurs in a particularly fragile natural system subjected to profound disturbances caused by diverse anthropic actions and therefore measures for the conservation of its habitat will be indispensable.

Cross-disciplinary information for understanding macroevolution

Many different macroevolutionary models can produce the same observations. Despite efforts in building more complex and realistic models, it may still be difficult to distinguish the processes that have generated the biodiversity we observe. In this opinion we argue that we can make new progress by reaching out across disciplines, relying on independent data and theory to constrain macroevolutionary inference. Using mainly paleontological insights and data, we illustrate how we can eliminate less plausible or implausible models, and/or parts of parameter space, while applying comparative phylogenetic approaches. We emphasize that such cross-disciplinary insights and data can be drawn between many other disciplines relevant to macroevolution. We urge cross-disciplinary training, and collaboration using common-use databases as a platform for increasing our understanding.

Phylogenetics and the Cenozoic radiation of lampreys

The development of a movable jaw is one of the most important transitions in the evolutionary history of animals.¹ Jawed vertebrates rapidly diversified after appearing approximately 470 million years ago. Today, only lampreys and hagfishes represent the once dominant jawless grade^2,3,4 and comprise less than 1% of living vertebrate species. Their relationship to other vertebrates ranks among the more contentious problems in animal phylogenetics.^{5,6,7,8,9,10,11,12} Further, the phylogenetic relationships within lampreys and hagfishes remain unclear,^13,14,15 and the ages of their living lineages are largely unexplored.^16,17 Because of their importance for the genomic and developmental changes that prefigured jawed vertebrate diversity,^18,19,20,21 the evolutionary history of lampreys and hagfishes is a major frontier of organismal biology. Of these two clades, lampreys²² are more ecologically diverse, exhibiting freshwater, anadromous, and fully marine forms, as well as parasitic and nonparasitic species.^23,24 Here, we present a new phylogeny and historical biogeographic reconstruction of all living lampreys. Whereas the early diversification of this clade tracks Pangaean fragmentation, lampreys also rapidly radiated in the northern hemisphere during the mid-Cretaceous and directly after the Cretaceous-Paleogene extinction. These radiations mirrored concurrent ones in other animals and plants and coincided with changes to lamprey ecology and feeding behavior. Our results suggest that 80% of living lamprey clades appeared in the last 20 million years of Earth history. Rather than gradually accumulating since the oldest stem-group forms appeared in the early Paleozoic, living lamprey biodiversity results from diversifications extending from the Cretaceous to present.

Bursts in skull evolution weakened with time

The skull shapes of mammals diversified more rapidly early in their history

Alternating regimes of shallow and deep-sea diversification explain a species-richness paradox in marine fishes

The deep sea contains a surprising diversity of life, including iconic fish groups such as anglerfishes and lanternfishes. Still, >65% of marine teleost fish species are restricted to the photic zone <200 m, which comprises less than 10% of the ocean's total volume. From a macroevolutionary perspective, this paradox may be explained by three hypotheses: 1) shallow water lineages have had more time to diversify than deep-sea lineages, 2) shallow water lineages have faster rates of speciation than deep-sea lineages, or 3) shallow-to-deep sea transition rates limit deep-sea richness. Here we use phylogenetic comparative methods to test among these three non-mutually exclusive hypotheses. While we found support for all hypotheses, the disparity in species richness is better described as the uneven outcome of alternating phases that favored shallow or deep diversification over the past 200 million y. Shallow marine teleosts became incredibly diverse 100 million y ago during a period of warm temperatures and high sea level, suggesting the importance of reefs and epicontinental settings. Conversely, deep-sea colonization and speciation was favored during brief episodes when cooling temperatures increased the efficiency of the ocean's carbon pump. Finally, time-variable ecological filters limited shallow-to-deep colonization for much of teleost history, which helped maintain higher shallow richness. A pelagic lifestyle and large jaws were associated with early deep-sea colonists, while a demersal lifestyle and a tapered body plan were typical of later colonists. Therefore, we also suggest that some hallmark characteristics of deep-sea fishes evolved prior to colonizing the deep sea.

Evolution of the ancestral mammalian karyotype and syntenic regions

Computational reconstruction of ancestral mammalian karyotypes revealed a comprehensive picture of the chromosome rearrangements that occurred over the evolutionary history of mammals. Ancient gene order, in some cases extending to full chromosomes, was found conserved for more than 300 My, demonstrating strong evolutionary constraint against rearrangements in some regions. Conserved segments of chromosomes are enriched for genes that control developmental processes. Therefore, Darwinian selection likely maintains ancient gene combinations while allowing for genomic innovations within or near chromosomal sites that break and rearrange over evolutionary time. The revealed relationship between the three-dimensional structure of chromosomes and the evolutionary stability of chromosome segments provides additional insights into the mechanisms of chromosome evolution and diseases associated with genome rearrangements.

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

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.

Predicting long pendant edges in model phylogenies, with applications to biodiversity and tree inference

In the simplest phylogenetic diversification model (the pure-birth Yule process), lineages split independently at a constant rate λ for time t. The length of a randomly chosen edge (either interior or pendant) in the resulting tree has an expected value that rapidly converges to 12λ as t grows, and thus is essentially independent of t. However, the behaviour of the length L of the longest pendant edge reveals remarkably different behaviour: L converges to t/2 as the expected number of leaves grows. Extending this model to allow an extinction rate μ (where μ < λ), we also establish a similar result for birth-death trees, except that t/2 is replaced by t/2 ⋅ (1 - μ/λ). This 'complete' tree may contain subtrees that have died out before time t; for the 'reduced tree' that just involves the leaves present at time t and their direct ancestors, the longest pendant edge length L again converges to t/2. Thus, there is likely to be at least one extant species whose associated pendant branch attaches to the tree approximately half-way back in time to the origin of the entire clade. We also briefly consider the length of the shortest edges. Our results are relevant to phylogenetic diversity indices in biodiversity conservation, and to quantifying the length of aligned sequences required to correctly infer a tree. We compare our theoretical results with simulations, and with the branch lengths from a recent phylogenetic tree of all mammals.

Multiple evolutionary origins and losses of tooth complexity in squamates

Teeth act as tools for acquiring and processing food, thus holding a prominent role in vertebrate evolution. In mammals, dental-dietary adaptations rely on tooth complexity variations controlled by cusp number and pattern. Complexity increase through cusp addition has dominated the diversification of mammals. However, studies of Mammalia alone cannot reveal patterns of tooth complexity conserved throughout vertebrate evolution. Here, we use morphometric and phylogenetic comparative methods across fossil and extant squamates to show they also repeatedly evolved increasingly complex teeth, but with more flexibility than mammals. Since the Late Jurassic, multiple-cusped teeth evolved over 20 times independently from a single-cusped common ancestor. Squamates frequently lost cusps and evolved varied multiple-cusped morphologies at heterogeneous rates. Tooth complexity evolved in correlation with changes in plant consumption, resulting in several major increases in speciation. Complex teeth played a critical role in vertebrate evolution outside Mammalia, with squamates exemplifying a more labile system of dental-dietary evolution.