Empirical and Bayesian approaches to fossil-only divergence times: A study across three reptile clades

Decoupling speciation and extinction reveals both abiotic and biotic drivers shaped 250 million years of diversity in crocodile-line archosaurs

Whereas living representatives of Pseudosuchia, crocodylians, number fewer than 30 species, more than 700 pseudosuchian species are known from their 250-million-year fossil record, displaying far greater ecomorphological diversity than their extant counterparts. With a new time-calibrated tree of >500 species, we use a phylogenetic framework to reveal that pseudosuchian evolutionary history and diversification dynamics were directly shaped by the interplay of abiotic and biotic processes over hundreds of millions of years, supported by information theory analyses. Speciation, but not extinction, is correlated with higher temperatures in terrestrial and marine lineages, with high sea level associated with heightened extinction in non-marine taxa. Low lineage diversity and increased speciation in non-marine species is consistent with opportunities for niche-filling, whereas increased competition may have led to elevated extinction rates. In marine lineages, competition via increased lineage diversity appears to have driven both speciation and extinction. Decoupling speciation and extinction, in combination with ecological partitioning, reveals a more complex picture of pseudosuchian evolution than previously understood. As the number of species threatened with extinction by anthropogenic climate change continues to rise, the fossil record provides a unique window into the drivers that led to clade success and those that may ultimately lead to extinction.

Locomotor and postural diversity among reptiles viewed through the prism of femoral microanatomy: Palaeobiological implications for some Permian and Mesozoic taxa

The water-to-land transition by the first tetrapod vertebrates represents a key stage in their evolution. Selection pressures exerted by this new environment on animals led to the emergence of new locomotor and postural strategies that favoured access to different ecological niches and contributed to their evolutionary success. Today, amniotes show great locomotor and postural diversity, particularly among Reptilia, whose extant representatives include parasagittally locomoting erect and crouched bipeds (birds), sub-parasagittal 'semi-erect' quadrupeds (crocodylians) and sprawling quadrupeds (squamates and turtles). But the different steps leading to such diversity remain enigmatic and the type of locomotion adopted by many extinct species raises questions. This is notably the case of certain Triassic taxa such as Euparkeria and Marasuchus. The exploration of the bone microanatomy in reptiles could help to overcome these uncertainties. Indeed, this locomotor and postural diversity is accompanied by great microanatomical disparity. On land, the bones of the appendicular skeleton support the weight of the body and are subject to multiple constraints that partly shape their external and internal morphology. Here we show how microanatomical parameters measured in cross-section, such as bone compactness or the position of the medullocortical transition, can be related to locomotion. We hypothesised that this could be due to variations in cortical thickness. Using statistical methods that take phylogeny into account (phylogenetic flexible discriminant analyses), we develop different models of locomotion from a sample of femur cross-sections from 51 reptile species. We use these models to infer locomotion and posture in 7 extinct reptile taxa for which they remain debated or not fully clear. Our models produced reliable inferences for taxa that preceded and followed the quadruped/biped and sprawling/erect transitions, notably within the Captorhinidae and Dinosauria. For taxa contemporary with these transitions, such as Terrestrisuchus and Marasuchus, the inferences are more questionable. We use linear models to investigate the effect of body mass and functional ecology on our inference models. We show that body mass seems to significantly impact our model predictions in most cases, unlike the functional ecology. Finally, we illustrate how taphonomic processes can impact certain microanatomical parameters, especially the eccentricity of the section, while addressing some other potential limitations of our methods. Our study provides insight into the evolution of enigmatic locomotion in various early reptiles. Our models and methods could be used by palaeontologists to infer the locomotion and posture in other extinct reptile taxa, especially when considered in combination with other lines of evidence.

Ecomorphological patterns in trigeminal canal branching among sauropsids reveal sensory shift in suchians

The vertebrate trigeminal nerve is the primary mediator of somatosensory information from nerve endings across the face, extending nerve branches through bony canals in the face and mandibles, terminating in sensory receptors. Reptiles evolved several extreme forms of cranial somatosensation in which enhanced trigeminal tissues are present in species engaging in unique mechanosensory behaviors. However, morphology varies by clade and ecology among reptiles. Few lineages approach the extreme degree of tactile somatosensation possessed by crocodylians, the only remaining members of a clade that underwent an ecological transition from the terrestrial to semiaquatic habitat, also evolving a specialized trigeminal system. It remains to be understood how trigeminal osteological correlates inform how adaptations for enhanced cranial sensation evolved in crocodylians. Here we identify an increase in sensory abilities in Early Jurassic crocodylomorphs, preceding the transitions to a semiaquatic habitat. Through quantification of trigeminal neurovascular canal branching patterns in an extant phylogenetic bracket we quantify and identify morphologies associated with sensory behaviors in representative fossil taxa, we find stepwise progression of increasing neurovascular canal density, complexity, and distribution from the primitive archosaurian to the derived crocodilian condition. Model-based inferences of sensory ecologies tested on quantified morphologies of extant taxa with known sensory behaviors indicate a parallel increase in sensory abilities among pseudosuchians. These findings establish patterns of reptile trigeminal ecomorphology, revealing evolutionary patterns of somatosensory ecology.

Total evidence phylogeny of platyrrhine primates and a comparison of undated and tip-dating approaches

There have been multiple published phylogenetic analyses of platyrrhine primates (New World monkeys) using both morphological and molecular data, but relatively few that have integrated both types of data into a total evidence approach. Here, we present phylogenetic analyses of recent and fossil platyrrhines, based on a total evidence data set of 418 morphological characters and 10.2 kilobases of DNA sequence data from 17 nuclear genes taken from previous studies, using undated and tip-dating approaches in a Bayesian framework. We compare the results of these analyses with molecular scaffold analyses using maximum parsimony and Bayesian approaches, and we use a formal information theoretic approach to identify unstable taxa. After a posteriori pruning of unstable taxa, the undated and tip-dating topologies appear congruent with recent molecular analyses and support largely similar relationships, with strong support for Stirtonia as a stem alouattine, Neosaimiri as a stem saimirine, Cebupithecia as a stem pitheciine, and Lagonimico as a stem callitrichid. Both analyses find three Greater Antillean subfossil platyrrhines (Xenothrix, Antillothrix, and Paralouatta) to form a clade that is related to Callicebus, congruent with a single dispersal event by the ancestor of this clade to the Greater Antilles. They also suggest that the fossil Proteropithecia may not be closely related to pitheciines, and that all known platyrrhines older than the Middle Miocene are stem taxa. Notably, the undated analysis found the Early Miocene Panamacebus (currently recognized as the oldest known cebid) to be unstable, and the tip-dating analysis placed it outside crown Platyrrhini. Our tip-dating analysis supports a late Oligocene or earliest Miocene (20.8-27.0 Ma) age for crown Platyrrhini, congruent with recent molecular clock analyses.

Integrative Phylogenetics: Tools for Palaeontologists to Explore the Tree of Life

The modern era of analytical and quantitative palaeobiology has only just begun, integrating methods such as morphological and molecular phylogenetics and divergence time estimation, as well as phenotypic and molecular rates of evolution. Calibrating the tree of life to geological time is at the nexus of many disparate disciplines, from palaeontology to molecular systematics and from geochronology to comparative genomics. Creating an evolutionary time scale of the major events that shaped biodiversity is key to all of these fields and draws from each of them. Different methodological approaches and data employed in various disciplines have traditionally made collaborative research efforts difficult among these disciplines. However, the development of new methods is bridging the historical gap between fields, providing a holistic perspective on organismal evolutionary history, integrating all of the available evidence from living and fossil species. Because phylogenies with only extant taxa do not contain enough information to either calibrate the tree of life or fully infer macroevolutionary dynamics, phylogenies should preferably include both extant and extinct taxa, which can only be achieved through the inclusion of phenotypic data. This integrative phylogenetic approach provides ample and novel opportunities for evolutionary biologists to benefit from palaeontological data to help establish an evolutionary time scale and to test core macroevolutionary hypotheses about the drivers of biological diversification across various dimensions of organisms.

Phylogenomic analyses of echinoid diversification prompt a re-evaluation of their fossil record

Echinoids are key components of modern marine ecosystems. Despite a remarkable fossil record, the emergence of their crown group is documented by few specimens of unclear affinities, rendering their early history uncertain. The origin of sand dollars, one of its most distinctive clades, is also unclear due to an unstable phylogenetic context. We employ 18 novel genomes and transcriptomes to build a phylogenomic dataset with a near-complete sampling of major lineages. With it, we revise the phylogeny and divergence times of echinoids, and place their history within the broader context of echinoderm evolution. We also introduce the concept of a chronospace - a multidimensional representation of node ages - and use it to explore methodological decisions involved in time calibrating phylogenies. We find the choice of clock model to have the strongest impact on divergence times, while the use of site-heterogeneous models and alternative node prior distributions show minimal effects. The choice of loci has an intermediate impact, affecting mostly deep Paleozoic nodes, for which clock-like genes recover dates more congruent with fossil evidence. Our results reveal that crown group echinoids originated in the Permian and diversified rapidly in the Triassic, despite the relative lack of fossil evidence for this early diversification. We also clarify the relationships between sand dollars and their close relatives and confidently date their origins to the Cretaceous, implying ghost ranges spanning approximately 50 million years, a remarkable discrepancy with their rich fossil record.

The impact of molecular data on the phylogenetic position of the putative oldest crown crocodilian and the age of the clade

The use of molecular data for living groups is vital for interpreting fossils, especially when morphology-only analyses retrieve problematic phylogenies for living forms. These topological discrepancies impact on the inferred phylogenetic position of many fossil taxa. In Crocodylia, morphology-based phylogenetic inferences differ fundamentally in placing Gavialis basal to all other living forms, whereas molecular data consistently unite it with crocodylids. The Cenomanian Portugalosuchus azenhae was recently described as the oldest crown crocodilian, with affinities to Gavialis, based on morphology-only analyses, thus representing a potentially important new molecular clock calibration. Here, we performed analyses incorporating DNA data into these morphological datasets, using scaffold and supermatrix (total evidence) approaches, in order to evaluate the position of basal crocodylians, including Portugalosuchus. Our analyses incorporating DNA data robustly recovered Portugalosuchus outside Crocodylia (as well as thoracosaurs, planocraniids and Borealosuchus spp.), questioning the status of Portugalosuchus as crown crocodilian and any future use as a node calibration in molecular clock studies. Finally, we discuss the impact of ambiguous fossil calibration and how, with the increasing size of phylogenomic datasets, the molecular scaffold might be an efficient (though imperfect) approximation of more rigorous but demanding supermatrix analyses.

Genomic diversity and Demographic History of the Dromiciops genus (Marsupialia: Microbiotheriidae)

Three orders represent the South American fauna of marsupials. Of these, Microbiotheria was until recently known as a monotypic genus with the only surviving species Dromiciops gliroides (monito del monte). The recent proposal of a new Dromiciops species (Dromiciops bicinovici), together with new information on the origin and diversification of living microbioterians has changed the prevailing paradigm around the evolutionary history of these emblematic marsupials. Here, we used a RADseq approach to test for evidence of admixture and past or current gene flow among both species of Dromiciops and evaluate the genetic structure within D. gliroides. We analyzed 127 samples of Dromiciops distributed across the known distribution range of both species. We also inferred the joint demographic history of these lineages, thus corroborating the status of D. bozinovici as a distinct species. Demographic history reconstruction indicated that D. bozinovici diverged from D. gliroides around 4my ago and has remained isolated and demographically stable ever since. In contrast, D. gliroides is subdivided into three subclades that experienced recent expansions and moderate gene flow among them (mostly from north to south). Furthermore, genetic distances among populations within D. gliroides were significantly correlated with geographic distances. These results suggest that some of the D. gliroides populations would have survived in glacial refuges, with posterior expansions after ice retreat. Our results have important implications for the systematics of the genus and have profound conservation consequences for the new species, especially considering the fragmentation level of the temperate rainforest.

Fossils improve phylogenetic analyses of morphological characters

Fossils provide our only direct window into evolutionary events in the distant past. Incorporating them into phylogenetic hypotheses of living clades can help time-calibrate divergences, as well as elucidate macroevolutionary dynamics. However, the effect fossils have on phylogenetic reconstruction from morphology remains controversial. The consequences of explicitly incorporating the stratigraphic ages of fossils using tip-dated inference are also unclear. Here, we use simulations to evaluate the performance of inference methods across different levels of fossil sampling and missing data. Our results show that fossil taxa improve phylogenetic analysis of morphological datasets, even when highly fragmentary. Irrespective of inference method, fossils improve the accuracy of phylogenies and increase the number of resolved nodes. They also induce the collapse of ancient and highly uncertain relationships that tend to be incorrectly resolved when sampling only extant taxa. Furthermore, tip-dated analyses under the fossilized birth-death process outperform undated methods of inference, demonstrating that the stratigraphic ages of fossils contain vital phylogenetic information. Fossils help to extract true phylogenetic signals from morphology, an effect that is mediated by both their distinctive morphology and their temporal information, and their incorporation in total-evidence phylogenetics is necessary to faithfully reconstruct evolutionary history.

A phylogenetic study of Haemocystidium parasites and other Haemosporida using complete mitochondrial genome sequences

Haemosporida are diverse vector-borne parasites associated with terrestrial vertebrates. Driven by the interest in species causing malaria (genus Plasmodium), the diversity of avian and mammalian haemosporidian species has been extensively studied, relying mostly on mitochondrial genes, particularly cytochrome b. However, parasites from reptiles have been neglected in biodiversity surveys. Reptilian haemosporidian parasites include Haemocystidium, a genus that shares morphological features with Plasmodium and Haemoproteus. Here, the first complete Haemocystidium mitochondrial DNA (mtDNA) genomes are studied. In particular, three mtDNA genomes from Haemocystidium spp. sampled in Africa, Oceania, and South America, are described. The Haemocystidium mtDNA genomes showed a high A + T content and a gene organization, including an extreme fragmentation of the rRNAs, found in other Haemosporida. These Haemocystidium mtDNA genomes were incorporated in phylogenetic and molecular clock analyses together with a representative sample of haemosporidian parasites from birds, mammals, and reptiles. The recovered phylogeny supported Haemocystidium as a monophyletic group apart from Plasmodium and other Haemosporida. Both the phylogenetic and molecular clock analyses yielded results consistent with a scenario in which haemosporidian parasites radiated with modern birds. Haemocystidium, like mammalian parasite clades, seems to originate from host switches by avian Haemosporida that allowed for the colonization of new vertebrate hosts. This hypothesis can be tested by investigating additional parasite species from all vertebrate hosts, particularly from reptiles. The mtDNA genomes reported here provide baseline data that can be used to scale up studies in haemosporidian parasites of reptiles using barcode approaches.

Bayesian Tip-Dated Phylogenetics in Paleontology: Topological Effects and Stratigraphic Fit

The incorporation of stratigraphic data into phylogenetic analysis has a long history of debate but is not currently standard practice for paleontologists. Bayesian tip-dated (or morphological clock) phylogenetic methods have returned these arguments to the spotlight, but how tip dating affects the recovery of evolutionary relationships has yet to be fully explored. Here I show, through analysis of several data sets with multiple phylogenetic methods, that topologies produced by tip dating are outliers as compared to topologies produced by parsimony and undated Bayesian methods, which retrieve broadly similar trees. Unsurprisingly, trees recovered by tip dating have better fit to stratigraphy than trees recovered by other methods under both the Gap Excess Ratio (GER) and the Stratigraphic Completeness Index (SCI). This is because trees with better stratigraphic fit are assigned a higher likelihood by the fossilized birth-death tree model. However, the degree to which the tree model favors tree topologies with high stratigraphic fit metrics is modulated by the diversification dynamics of the group under investigation. In particular, when net diversification rate is low, the tree model favors trees with a higher GER compared to when net diversification rate is high. Differences in stratigraphic fit and tree topology between tip dating and other methods are concentrated in parts of the tree with weaker character signal, as shown by successive deletion of the most incomplete taxa from two data sets. These results show that tip dating incorporates stratigraphic data in an intuitive way, with good stratigraphic fit an expectation that can be overturned by strong evidence from character data. [fossilized birth-death; fossils; missing data; morphological clock; morphology; parsimony; phylogenetics.].

First evidence for a latitudinal body mass effect in extant Crocodylia and the relationships of their reproductive characters

Relationships between distribution patterns and body size have been documented in many endothermic taxa. However, the evidence for these trends in ectotherms generally is equivocal, and there have been no studies of effects in crocodylians specifically. Here, we examine the relationship between latitudinal distribution and body mass in 20 extant species of crocodylians, as well as the relationships between seven important reproductive variables. Using phylogenetically independent contrasts to inform generalized linear models, we provide the first evidence of a latitudinal effect on adult female body mass in crocodylians. In addition, we explore the relationships between reproductive variables including egg mass, hatchling mass and clutch size. We report no correlation between egg mass and clutch size, upholding previously reported within-species trends. We also find no evidence of a correlation between measures of latitudinal range and incubation temperature, contrasting with the trends found in turtles.

Phylogeny of the Eurasian Wren Nannus troglodytes (Aves: Passeriformes: Troglodytidae) reveals deep and complex diversification patterns of Ibero-Maghrebian and Cyrenaican populations

The Mediterranean Basin represents a Global Biodiversity Hotspot where many organisms show high inter- and intraspecific differentiation. Extant phylogeographic patterns of terrestrial circum-Mediterranean faunas were mainly shaped through Pleistocene range shifts and range fragmentations due to retreat into different glacial refugia. Thus, several extant Mediterranean bird species have diversified by surviving glaciations in different hospitable refugia and subsequently expanded their distribution ranges during the Holocene. Such a scenario was also suggested for the Eurasian Wren (Nannus troglodytes) despite the lack of genetic data for most Mediterranean subspecies. Our phylogenetic multi-locus analysis comprised 18 out of 28 currently accepted subspecies of N. troglodytes, including all but one subspecies which are present in the Mediterranean Basin. The resulting phylogenetic reconstruction dated the onset of the entire Holarctic radiation of three Nannus species to the early Pleistocene. In the Eurasian Wren, two North African subspecies represented separate basal lineages from the Maghreb (N. t. kabylorum) and from the Libyan Cyrenaica (N. t. juniperi), being only distantly related to other Mediterranean populations. Although N. troglodytes appeared to be paraphyletic with respect to the Nearctic Winter Wren (N. hiemalis), respective nodes did not receive strong statistical support. In contrast, paraphyly of the Ibero-Maghrebian taxon N. t. kabylorum was strongly supported. Southern Iberian populations of N. t. kabylorum did not clade with Maghrebian populations of the same subspecies but formed a sister clade to a highly diverse European clade (including nominate N. t. troglodytes and eight further taxa). In accordance with a pattern also found in other birds, Eurasian populations were split into a western clade (Europe, Caucasus) and an eastern clade (Central Asia, Sino-Himalayas, East Asia). This complex phylogeographic pattern revealed cryptic diversification in N. troglodytes, especially in the Iberio-Maghrebian region.

The multi-peak adaptive landscape of crocodylomorph body size evolution

BACKGROUND

Little is known about the long-term patterns of body size evolution in Crocodylomorpha, the > 200-million-year-old group that includes living crocodylians and their extinct relatives. Extant crocodylians are mostly large-bodied (3-7 m) predators. However, extinct crocodylomorphs exhibit a wider range of phenotypes, and many of the earliest taxa were much smaller (< 1.2 m). This suggests a pattern of size increase through time that could be caused by multi-lineage evolutionary trends of size increase or by selective extinction of small-bodied species. Here, we characterise patterns of crocodylomorph body size evolution using a model fitting-approach (with cranial measurements serving as proxies). We also estimate body size disparity through time and quantitatively test hypotheses of biotic and abiotic factors as potential drivers of crocodylomorph body size evolution.

RESULTS

Crocodylomorphs reached an early peak in body size disparity during the Late Jurassic, and underwent an essentially continual decline since then. A multi-peak Ornstein-Uhlenbeck model outperforms all other evolutionary models fitted to our data (including both uniform and non-uniform), indicating that the macroevolutionary dynamics of crocodylomorph body size are better described within the concept of an adaptive landscape, with most body size variation emerging after shifts to new macroevolutionary regimes (analogous to adaptive zones). We did not find support for a consistent evolutionary trend towards larger sizes among lineages (i.e., Cope's rule), or strong correlations of body size with climate. Instead, the intermediate to large body sizes of some crocodylomorphs are better explained by group-specific adaptations. In particular, the evolution of a more aquatic lifestyle (especially marine) correlates with increases in average body size, though not without exceptions.

CONCLUSIONS

Shifts between macroevolutionary regimes provide a better explanation of crocodylomorph body size evolution on large phylogenetic and temporal scales, suggesting a central role for lineage-specific adaptations rather than climatic forcing. Shifts leading to larger body sizes occurred in most aquatic and semi-aquatic groups. This, combined with extinctions of groups occupying smaller body size regimes (particularly during the Late Cretaceous and Cenozoic), gave rise to the upward-shifted body size distribution of extant crocodylomorphs compared to their smaller-bodied terrestrial ancestors.

Evolutionary structure and timing of major habitat shifts in Crocodylomorpha

Extant crocodylomorphs are semiaquatic ambush predators largely restricted to freshwater or estuarine environments, but the group is ancestrally terrestrial and inhabited a variety of ecosystems in the past. Despite its rich ecological history, little effort has focused on elucidating the historical pattern of ecological transitions in the group. Traditional views suggested a single shift from terrestrial to aquatic in the Early Jurassic. However, new fossil discoveries and phylogenetic analyses tend to imply a multiple-shift model. Here we estimate ancestral habitats across a comprehensive phylogeny and show at least three independent shifts from terrestrial to aquatic and numerous other habitat transitions. Neosuchians first invade freshwater habitats in the Jurassic, with up to four subsequent shifts into the marine realm. Thalattosuchians first appear in marine habitats in the Early Jurassic. Freshwater semiaquatic mahajangasuchids are derived from otherwise terrestrial notosuchians. Within nearly all marine groups, some species return to freshwater environments. Only twice have crocodylomorphs reverted from aquatic to terrestrial habitats, both within the crown group. All living non-alligatorid crocodylians have a keratinised tongue with salt-excreting glands, but the lack of osteological correlates for these adaptations complicates pinpointing their evolutionary origin or loss. Based on the pattern of transitions to the marine realm, our analysis suggests at least four independent origins of saltwater tolerance in Crocodylomorpha.

Tip-dating and homoplasy: reconciling the shallow molecular divergences of modern gharials with their long fossil record

Simultaneously analysing morphological, molecular and stratigraphic data suggests a potential resolution to a major remaining inconsistency in crocodylian evolution. The ancient, long-snouted thoracosaurs have always been placed near the Indian gharial Gavialis, but their antiquity (ca 72 Ma) is highly incongruous with genomic evidence for the young age of the Gavialis lineage (ca 40 Ma). We reconcile this contradiction with an updated morphological dataset and novel analysis, and demonstrate that thoracosaurs are an ancient iteration of long-snouted stem crocodylians unrelated to modern gharials. The extensive similarities between thoracosaurs and Gavialis are shown to be an almost 'perfect storm' of homoplasy, combining convergent adaptions to fish-eating, as well resemblances between genuinely primitive traits (thoracosaurs) and atavisms (Gavialis). Phylogenetic methods that ignore stratigraphy (parsimony and undated Bayesian methods) are unable to tease apart these similarities and invariably unite thoracosaurs and Gavialis. However, tip-dated Bayesian approaches additionally consider the large temporal gap separating ancient (thoracosaurs) and modern (Gavialis) iterations of similar long-snouted crocodyliforms. These analyses robustly favour a phylogeny which places thoracosaurs basal to crocodylians, far removed from modern gharials, which accordingly are a very young radiation. This phylogenetic uncoupling of ancient and modern gharial-like crocs is more consistent with molecular clock divergence estimates, and also the bulk of the crocodylian fossil record (e.g. all unequivocal gharial fossils are very young). Provided that the priors and models attribute appropriate relative weights to the morphological and stratigraphic signals-an issue that requires investigation-tip-dating approaches are potentially better able to detect homoplasy and improve inferences about phylogenetic relationships, character evolution and divergence dates.