Testing extinction events and temporal shifts in diversification and fossilization rates through the skyline Fossilized Birth-Death (FBD) model: The example of some mid-Permian synapsid extinctions

In the last decade, the Fossilized Birth-Death (FBD) process has yielded interesting clues about the evolution of biodiversity through time. To facilitate such studies, we extend our method to compute the probability density of phylogenetic trees of extant and extinct taxa in which the only temporal information is provided by the fossil ages (i.e. without the divergence times) in order to deal with the piecewise constant FBD process, known as the "skyline FBD", which allows rates to change between pre-defined time intervals, as well as modelling extinction events at the bounds of these intervals. We develop approaches based on this method to assess hypotheses about the diversification process and to answer questions such as "Does a mass extinction occur at this time?" or "Is there a change in the fossilization rate between two given periods?". Our software can also yield Bayesian and maximum-likelihood estimates of the parameters of the skyline FBD model under various constraints. These approaches are applied to a simulated dataset in order to test their ability to answer the questions above. Finally, we study an updated dataset of Permo-Carboniferous synapsids to get additional insights into the dynamics of biodiversity change in three clades (Ophiacodontidae, Edaphosauridae and Sphenacodontidae) in the Pennsylvanian (Late Carboniferous) and Cisuralian (Early Permian), and to assess support for end-Sakmarian (or Artinskian) and end-Cisuralian mass extinction events discussed in previous studies.

Evolution of posture in amniotes-Diving into the trabecular architecture of the femoral head

Extant amniotes show remarkable postural diversity. Broadly speaking, limbs with erect (strongly adducted, more vertically oriented) posture are found in mammals that are particularly heavy (graviportal) or show good running skills (cursorial), while crouched (highly flexed) limbs are found in taxa with more generalized locomotion. In Reptilia, crocodylians have a "semi-erect" (somewhat adducted) posture, birds have more crouched limbs and lepidosaurs have sprawling (well-abducted) limbs. Both synapsids and reptiles underwent a postural transition from sprawling to more erect limbs during the Mesozoic Era. In Reptilia, this postural change is prominent among archosauriforms in the Triassic Period. However, limb posture in many key Triassic taxa remains poorly known. In Synapsida, the chronology of this transition is less clear, and competing hypotheses exist. On land, the limb bones are subject to various stresses related to body support that partly shape their external and internal morphology. Indeed, bone trabeculae (lattice-like bony struts that form the spongy bone tissue) tend to orient themselves along lines of force. Here, we study the link between femoral posture and the femoral trabecular architecture using phylogenetic generalized least squares. We show that microanatomical parameters measured on bone cubes extracted from the femoral head of a sample of amniote femora depend strongly on body mass, but not on femoral posture or lifestyle. We reconstruct ancestral states of femoral posture and various microanatomical parameters to study the "sprawling-to-erect" transition in reptiles and synapsids, and obtain conflicting results. We tentatively infer femoral posture in several hypothetical ancestors using phylogenetic flexible discriminant analysis from maximum likelihood estimates of the microanatomical parameters. In general, the trabecular network of the femoral head is not a good indicator of femoral posture. However, ancestral state reconstruction methods hold great promise for advancing our understanding of the evolution of posture in amniotes.

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.

Distributions of extinction times from fossil ages and tree topologies: the example of mid-Permian synapsid extinctions

Given a phylogenetic tree that includes only extinct, or a mix of extinct and extant taxa, where at least some fossil data are available, we present a method to compute the distribution of the extinction time of a given set of taxa under the Fossilized-Birth-Death model. Our approach differs from the previous ones in that it takes into account (i) the possibility that the taxa or the clade considered may diversify before going extinct and (ii) the whole phylogenetic tree to estimate extinction times, whilst previous methods do not consider the diversification process and deal with each branch independently. Because of this, our method can estimate extinction times of lineages represented by a single fossil, provided that they belong to a clade that includes other fossil occurrences. We assess and compare our new approach with a standard previous one using simulated data. Results show that our method provides more accurate confidence intervals. This new approach is applied to the study of the extinction time of three Permo-Carboniferous synapsid taxa (Ophiacodontidae, Edaphosauridae, and Sphenacodontidae) that are thought to have disappeared toward the end of the Cisuralian (early Permian), or possibly shortly thereafter. The timing of extinctions of these three taxa and of their component lineages supports the idea that the biological crisis in the late Kungurian/early Roadian consisted of a progressive decline in biodiversity throughout the Kungurian.

Evolution of metabolic scaling among the tetrapod: Effect of phylogeny, the geologic time of class formation and uniformity of species within a class

The metabolic scaling in the animal has been discussed for over 90 years, but no consensus has been reached. Our analysis of 2,126 species of vertebrates reveals a significant allometric exponent heterogeneity. We show that classes of terrestrial vertebrates exhibit the evolution of metabolic scaling. Both the allometric coefficient "a" and the allometric exponent "b" change naturally, but differently depending on the geological time of group formation. The allometric coefficient "a" shows the measure of the evolutionary development of systems that forms resting metabolism in animals. Endothermic classes, such as birds and mammals, have a metabolic rate that is in an order of magnitude higher than that in ectothermic classes, including amphibians and reptiles. In the terrestrial vertebrate phylogeny, we find that the metabolic scaling is characterized by three main allometric exponent values: b = 3/4 (mammals), b > 3/4 (ectotherms, such as amphibians and reptiles), and b < 3/4 (birds). The heterogeneity of the allometric exponent is a natural phenomenon associated with the general evolution of vertebrates. The scaling factor decreases depending on both the external design and the size (birds vs mammals) of the animal. The metabolic rate and uniformity of species within a class increase as the geological start date of formation of the class approaches the present time. The higher the mass-specific standard metabolic rate in the class, the slower metabolic rate grows with increasing body size in this class. Our results lay the groundwork for further exploration of the evolutionary and ecological aspects of the development of metabolic scaling in animals. This article is protected by copyright. All rights reserved.

Tracking the Origin and Early Evolution of Reptiles

The origin of Reptilia and the biostratigraphic and palaeobiogeographic distribution of its early representatives are still poorly understood. An independent source of information may come from the extensive Carboniferous footprint record of reptiles, which is arguably richer and more complete than the skeletal record. Nevertheless, previous studies often failed to provide useful information because they were based on poorly preserved material and/or characters non-exclusive of reptile tracks. In fact, a large part of the supposed early reptile tracks can be assigned to the anamniote ichnotaxon Hylopus hardingi. Here, we revise the ichnotaxon Hylopus hardingi based on anatomy-consistent material, attribute it to anamniote reptiliomorphs, and distinguish it from Notalacerta missouriensis, the earliest ichnotaxon that can be attributed to reptiles, and the somewhat younger Varanopus microdactylus (attributed to parareptiles, such as bolosaurians) and Dromopus lacertoides (attributed to araeoscelid reptiles and non-varanodontine varanopids). These attributions are based on correlating morphofunctional features of tracks and skeletons. Multivariate analysis of trackway parameters indicates that the late Bashkirian Notalacerta missouriensis and Hylopus hardingi differ markedly in their trackway patterns from Late Mississippian Hylopus hardingi and Late Pennsylvanian reptile tracks, which appear to share a derived amniote-like type of gait. While the first occurrence/appearance of reptile tracks in the tetrapod footprint record during the late Bashkirian corresponds to the first occurrence/appearance of reptiles in the skeletal record, footprints significantly enlarge the paleobiogeographic distribution of the group, suggesting an earlier radiation of reptiles during the Bashkirian throughout North America and possibly North Africa. Dromopus appeared in the Kasimovian together with the diapsid group Araeoscelidia, but footprints from Western-European occurrences enlarge the paleobiogeographic distribution of diapsids and varanopids. Varanopus and bolosaurian parareptiles appear in the Gzhelian of North America. Older parareptiles are, however, known from the late Moscovian. In all, the footprint record of early reptiles supplements the skeletal record, suggesting possible future lines of research.

Evolutionary Change in Locomotion Close to the Origin of Amniotes Inferred From Trackway Data in an Ancestral State Reconstruction Approach

Among amniote and non-amniote tetrapod trackways from late Carboniferous to early Permian deposits, certain trackway measures vary notably. Some of this variability can be attributed to evolutionary changes in trackmaker anatomy and locomotion style close to the origin of amniotes. Here we demonstrate that steps in early amniote locomotion evolution can be addressed by applying methods of ancestral state reconstruction on trackway data – a novel approach in tetrapod ichnology. Based on (a) measurements of 186 trackways referred to the Carboniferous and early Permian ichnogenera Batrachichnus, Limnopus, Hylopus, Amphisauropus, Matthewichnus, Ichniotherium, Dimetropus, Tambachichnium, Erpetopus, Varanopus, Hyloidichnus, Notalacerta and Dromopus, (b) correlation of these ichnotaxa with specific groups of amphibian, reptiliomorph, synapsid, and reptilian trackmakers based on imprint morphology and (c) known skeletal-morphology-based phylogenies of the supposed trackmakers, we infer ancestral states for functionally controlled trackway measures in a maximum likelihood approach. The most notable finding of our analysis is a concordant change in trackway parameters within a series of ancestral amniote trackmakers, which reflects an evolutionary change in locomotion: In the ancestors of amniotes and diadectomorphs, an increase in body size was accompanied by a decrease in (normalized) gauge width and glenoacetabular length and by a change in imprint orientation toward a more trackway-parallel and forward-pointing condition. In the subsequent evolution of diadectomorph, synapsid and reptilian trackmakers after the diversification of the clades Cotylosauria (Amniota + Diadectomorpha) and Amniota, stride length increased whereas gauges decreased further or remained relatively narrow within most lineages. In accordance with this conspicuous pattern of evolutionary change in trackway measures, we interpret the body size increase as an underlying factor that triggered the reorganization of the locomotion apparatus. The secondary increase in stride length, which occurred convergently within distinct groups, is interpreted as an increase in locomotion capability when the benefits of reorganization came into effect. The track-trackmaker pair of Ichniotherium sphaerodactylum and Orobates pabsti from the early Permian Bromacker locality of the Thuringian Forest, proposed in earlier studies as a suitable ancestral amniote track-trackmaker model, fits relatively well with our modeled last common ancestor of amniotes – with the caveat that the Bromacker material is younger and some of the similarities appear to be due to convergence.

An Evolutionary Microcircuit Approach to the Neural Basis of High Dimensional Sensory Processing in Olfaction

Odor stimuli consist of thousands of possible molecules, each molecule with many different properties, each property a dimension of the stimulus. Processing these high dimensional stimuli would appear to require many stages in the brain to reach odor perception, yet, in mammals, after the sensory receptors this is accomplished through only two regions, the olfactory bulb and olfactory cortex. We take a first step toward a fundamental understanding by identifying the sequence of local operations carried out by microcircuits in the pathway. Parallel research provided strong evidence that processed odor information is spatial representations of odor molecules that constitute odor images in the olfactory bulb and odor objects in olfactory cortex. Paleontology provides a unique advantage with evolutionary insights providing evidence that the basic architecture of the olfactory pathway almost from the start ∼330 million years ago (mya) has included an overwhelming input from olfactory sensory neurons combined with a large olfactory bulb and olfactory cortex to process that input, driven by olfactory receptor gene duplications. We identify a sequence of over 20 microcircuits that are involved, and expand on results of research on several microcircuits that give the best insights thus far into the nature of the high dimensional processing.

Raising names from the dead: A time-calibrated phylogeny of frog shells (Bursidae, Tonnoidea, Gastropoda) using mitogenomic data

With 59 Recent species, Bursidae, known as «frog shells», are a small but widely distributed group of tropical and subtropical gastropods that are most diverse in the Indo-West Pacific. The present study is aimed at reconstructing phylogenetic relationships of bursid gastropods based on extensive and representative taxon sampling. Five genetic markers (cytochrome c oxidase subunit I (cox1), 16 s and 12 s rRNA mitochondrial genes, 28 s rRNA and Histone H3 nuclear gene) were sequenced for over 30 species in every known genus but Crossata. Furthermore, we sequenced the complete mt-genome of 9 species (10 specimens) (Aspa marginata, Marsupina bufo, Korrigania quirihorai, Korrigania fijiensis, Tutufa rubeta, Bursa lamarckii, Lampasopsis rhodostoma (twice), Bufonaria perelegans and Bursa aff. tuberosissima). Our analysis recovered Bursidae as a monophyletic group, whereas the genus Bursa was found to be polyphyletic. The genera Talisman and Dulcerana are resurrected and the genera Alanbeuella gen. nov. and Korrigania gen. nov. are described. Dating analysis using 21 extinct taxa for node and simplified tip calibrations was performed, showing a diversification of the group in two phases. Diversification may be linked to tectonic events leading to biodiversity relocation from the western Tethys toward the Indo-Pacific.

Rates and Rocks: Strengths and Weaknesses of Molecular Dating Methods

I present here an in-depth, although non-exhaustive, review of two topics in molecular dating. Clock models, which describe the evolution of the rate of evolution, are considered first. Some of the shortcomings of popular approaches-uncorrelated clock models in particular-are presented and discussed. Autocorrelated models are shown to be more reasonable from a biological perspective. Some of the most recent autocorrelated models also rely on a coherent treatment of instantaneous and average substitution rates while previous models are based on implicit approximations. Second, I provide a brief overview of the processes involved in collecting and preparing fossil data. I then review the main techniques that use this data for calibrating the molecular clock. I argue that, in its current form, the fossilized birth-death process relies on assumptions about the mechanisms underlying fossilization and the data collection process that may negatively impact the date estimates. Node-dating approaches make better use of the data available, even though they rest on paleontologists' intervention to prepare raw fossil data. Altogether, this study provides indications that may help practitioners in selecting appropriate methods for molecular dating. It will also hopefully participate in defining the contour of future methodological developments in the field.

Quantifying the Error of Secondary vs. Distant Primary Calibrations in a Simulated Environment

Using calibrations to obtain absolute divergence times is standard practice in molecular clock studies. While the use of primary (e.g., fossil) calibrations is preferred, this approach can be limiting because of their rarity in fast-growing datasets. Thus, alternatives need to be explored, such as the use of secondary (molecularly-derived) calibrations that can anchor a timetree in a larger number of nodes. However, the use of secondary calibrations has been discouraged in the past because of concerns in the error rates of the node estimates they produce with an apparent high precision. Here, we quantify the amount of errors in estimates produced by the use of secondary calibrations relative to true times and primary calibrations placed on distant nodes. We find that, overall, the inaccuracies in estimates based on secondary calibrations are predictable and mirror errors associated with primary calibrations and their confidence intervals. Additionally, we find comparable error rates in estimated times from secondary calibrations and distant primary calibrations, although the precision of estimates derived from distant primary calibrations is roughly twice as good as that of estimates derived from secondary calibrations. This suggests that increasing dataset size to include primary calibrations may produce divergence times that are about as accurate as those from secondary calibrations, albeit with a higher precision. Overall, our results suggest that secondary calibrations may be useful to explore the parameter space of plausible evolutionary scenarios when compared to time estimates obtained with distant primary calibrations.