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.

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.

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.

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.

The cecal appendix is correlated with greater maximal longevity in mammals.

The cecal appendix had been considered as a useless vestige since Darwin's work, but recent research questioned this idea demonstrating that the cecal appendix appeared among the mammals at least 80 million years ago and has made multiple and independent appearances without any obvious correlation with diet, social life, ecology, or size of the cecum. However, functions and probable selective advantage conferred by this anatomical structure still remain enigmatic. We found, through analyses of data on 258 mammalian species, that cecal appendix presence is correlated with increased maximal observed longevity. This is the first demonstration of a correlation between cecal appendix presence and life history. Interestingly, the classical evolutionary theory of aging that predicts an increased longevity when the extrinsic mortality is reduced has been questioned several times, but recent comparative studies asserted its validity in the taxa, which experience age-dependent and density-dependent mortality, as in mammals. Thus, the cecal appendix may contribute to the increase in longevity through a reduction of extrinsic mortality. A lower risk of fatal infectious diarrhea is one of the most plausible hypotheses that could explain it. However, several hypotheses coexist about the possible functions of the cecal appendix, and our results provide new insights about this much-disputed question. In addition, we show that the cecal appendix arose at least 16 times and was lost only once during the evolutionary history of the considered mammals, an asymmetry that supports the existence of a positive selective of this structure.

A microanatomical and histological study of the scales of the Devonian sarcopterygian Miguashaia bureaui and the evolution of the squamation in coelacanths.

Coelacanths have traditionally been described as morphologically conservative throughout their long evolutionary history, which spans more than 400 million years. After an initial burst during the Devonian, a morphological stasis was long thought to have prevailed since the Carboniferous, as shown by the extant Latimeria. New fossil discoveries have challenged this view, with punctual and sometimes unusual departures from the general coelacanth Bauplan. The dermal skeleton is considered to represent one, if not the main, example of morphological stasis in coelacanth evolution and as a consequence, has remained poorly surveyed. The lack of palaeohistological data on the dermoskeleton has resulted in a poor understanding of the early establishment and evolution of the coelacanth squamation. Here we describe the scales of Miguashaia bureaui from the Upper Devonian of Miguasha, Québec (Canada), revealing histological data for a Palaeozoic coelacanth in great detail and adding to our knowledge on the dermal skeleton of sarcopterygians. Miguashaia displays rounded scales ornamented by tubercules and narrow ridges made of dentine and capped with enamel. At least two generations of superimposed odontodes occur, which is reminiscent of the primitive condition of stem osteichthyans like Andreolepis or Lophosteus, and onychodonts like Selenodus. The middle vascular layer is well developed and shows traces of osteonal remodelling. The basal plate consists of a fully mineralised lamellar bone with a repetitive rotation pattern every five layers indicating a twisted plywood-like arrangement of the collagen plies. Comparisons with the extant Latimeria and other extinct taxa show that these features are consistently conserved across coelacanth evolution with only minute changes in certain taxa. The morphological and histological features displayed in the scales of Miguashaia enable us to draw a comprehensive picture of the onset of the coelacanth squamation and to propose and discuss evolutionary scenarios for the coelacanth dermoskeleton.

Exact Distribution of Divergence Times from Fossil Ages and Tree Topologies.

Being given a phylogenetic tree of both extant and extinct taxa in which the fossil ages are the only temporal information (namely, in which divergence times are considered unknown), we provide a method to compute the exact probability distribution of any divergence time of the tree with regard to any speciation (cladogenesis), extinction, and fossilization rates under the Fossilized Birth-Death model. We use this new method to obtain a probability distribution for the age of Amniota (the synapsid/sauropsid or bird/mammal divergence), one of the most-frequently used dating constraints. Our results suggest an older age (between about 322 and 340 Ma) than has been assumed by most studies that have used this constraint (which typically assumed a best estimate around 310-315 Ma) and provide, for the first time, a method to compute the shape of the probability density for this divergence time. [Divergence times; fossil ages; fossilized birth-death model; probability distribution.].

On trends and patterns in macroevolution: Williston's law and the branchiostegal series of extant and extinct osteichthyans.

BACKGROUND: The branchiostegal series consists of an alignment of bony elements in the posterior portion of the skull of osteichthyan vertebrates. We trace the evolution of the number of elements in a comprehensive survey that includes 440 extant and 66 extinct species. Using a newly updated actinopterygian tree in combination with phylogenetic comparative analyses, we test whether osteichthyan branchiostegals follow an evolutionary trend under 'Williston's law', which postulates that osteichthyan lineages experienced a reduction of bony elements over time. RESULTS: We detected no overall macroevolutionary trend in branchiostegal numbers, providing no support for 'Williston's law'. This result is robust to the subsampling of palaeontological data, but the estimation of the model parameters is much more ambiguous. CONCLUSIONS: We find substantial evidence for a macroevolutionary dynamic favouring an 'early burst' of trait evolution over alternative models. Our study highlights the challenges of accurately reconstructing macroevolutionary dynamics even with large amounts of data about extant and extinct taxa.

Parsimony-based test for identifying changes in evolutionary trends for quantitative characters: implications for the origin of the amniotic egg.

The origin of the amniotic egg was a major event in vertebrate evolution and is thought to have contributed to the spectacular evolutionary radiation of amniotes. We test one of the most popular scenarios proposed by Carroll in 1970 to explain the origin of the amniotic egg using a novel method based on an asymmetric version of linear parsimony (aka Wagner parsimony) for identifying the most parsimonious split of a tree into two parts between which the evolution of the character is allowed to differ. The new method evaluates the cost of splitting a phylogenetic tree at a given node as the integral, over all pairs of asymmetry parameters, of the most parsimonious costs that can be achieved by using the first parameter on the subtree pending from this node and the second parameter elsewhere. By testing all the nodes, we then obtain the most parsimonious split of a tree with regard to the character values at its tips. Among the nine trees and two characters tested, our method yields a total of 517 parsimonious trend changes in Permo-Carboniferous stegocephalians, a single one of which occurs in a part of the tree (among stem-amniotes) where Carroll's scenario predicts that there should have been distinct changes in body size evolutionary trends. This refutes the scenario because the amniote stem does not appear to have elevated rates of evolutionary trend shifts. Our nodal body size estimates offer less discriminating power, but they likewise fail to find strong support for Carroll's scenario.

Likelihood of Tree Topologies with Fossils and Diversification Rate Estimation.

Since the diversification process cannot be directly observed at the human scale, it has to be studied from the information available, namely the extant taxa and the fossil record. In this sense, phylogenetic trees including both extant taxa and fossils are the most complete representations of the diversification process that one can get. Such phylogenetic trees can be reconstructed from molecular and morphological data, to some extent. Among the temporal information of such phylogenetic trees, fossil ages are by far the most precisely known (divergence times are inferences calibrated mostly with fossils). We propose here a method to compute the likelihood of a phylogenetic tree with fossils in which the only considered time information is the fossil ages, and apply it to the estimation of the diversification rates from such data. Since it is required in our computation, we provide a method for determining the probability of a tree topology under the standard diversification model. Testing our approach on simulated data shows that the maximum likelihood rate estimates from the phylogenetic tree topology and the fossil dates are almost as accurate as those obtained by taking into account all the data, including the divergence times. Moreover, they are substantially more accurate than the estimates obtained only from the exact divergence times (without taking into account the fossil record). We also provide an empirical example composed of 50 Permo-Carboniferous eupelycosaur (early synapsid) taxa ranging in age from about 315 Ma (Late Carboniferous) to 270 Ma (shortly after the end of the Early Permian). Our analyses suggest a speciation (cladogenesis, or birth) rate of about 0.1 per lineage and per myr, a marginally lower extinction rate, and a considerable hidden paleobiodiversity of early synapsids. [Extinction rate; fossil ages; maximum likelihood estimation; speciation rate.].

Testing for Depéret's Rule (Body Size Increase) in Mammals using Combined Extinct and Extant Data.

Whether or not evolutionary lineages in general show a tendency to increase in body size has often been discussed. This tendency has been dubbed "Cope's rule" but because Cope never hypothesized it, we suggest renaming it after Depéret, who formulated it clearly in 1907. Depéret's rule has traditionally been studied using fossil data, but more recently a number of studies have used present-day species. While several paleontological studies of Cenozoic placental mammals have found support for increasing body size, most studies of extant placentals have failed to detect such a trend. Here, we present a method to combine information from present-day species with fossil data in a Bayesian phylogenetic framework. We apply the method to body mass estimates of a large number of extant and extinct mammal species, and find strong support for Depéret's rule. The tendency for size increase appears to be driven not by evolution toward larger size in established species, but by processes related to the emergence of new species. Our analysis shows that complementary data from extant and extinct species can greatly improve inference of macroevolutionary processes.

Evolution of organogenesis and the origin of altriciality in mammals.

Mammals feature not only great phenotypic disparity, but also diverse growth and life history patterns, especially in maturity level at birth, ranging from altriciality to precocity. Gestation length, morphology at birth, and other markers of life history are fundamental to our understanding of mammalian evolution. Based on the first synthesis of embryological data and the study of new ontogenetic series, we reconstructed estimates of the ancestral chronology of organogenesis and life-history modes in placental mammals. We found that the ancestor of marsupial and placental mammals was placental-like at birth but had a long, marsupial-like infancy. We hypothesize that mammalian viviparity might have evolved in association with the extension of growth after birth, enabled through lactation, and that mammalian altriciality is inherited from the earliest amniotes. The precocial lifestyle of extant sauropsids and that of many placental mammals were acquired secondarily. We base our conclusions on the best estimates and provide a comprehensive discussion on the methods used and the limitations of our dataset. We provide the most comprehensive embryological dataset ever published, "rescue" old literature sources, and apply available methods and illustrate thus an approach on how to investigate comparatively organogenesis in macroevolution.

Gradual adaptation of bone structure to aquatic lifestyle in extinct sloths from Peru.

Non-pathological densification (osteosclerosis) and swelling (pachyostosis) of bones are the main modifications affecting the skeleton of land vertebrates (tetrapods) that returned to water. However, a precise temporal calibration of the acquisition of such adaptations is still wanting. Here, we assess the timing of such acquisition using the aquatic sloth Thalassocnus, from the Neogene of the Pisco Formation, Peru. This genus is represented by five species occurring in successive vertebrate-bearing horizons of distinct ages. It yields the most detailed data about the gradual acquisition of aquatic adaptations among tetrapods, in displaying increasing osteosclerosis and pachyostosis through time. Such modifications, reflecting a shift in the habitat from terrestrial to aquatic, occurred over a short geological time span (ca 4 Myr). Otherwise, the bones of terrestrial pilosans (sloths and anteaters) are much more compact than the mean mammalian condition, which suggests that the osteosclerosis of Thalassocnus may represent an exaptation.

Assessment of modularity in the urodele skull: an exploratory analysis using ossification sequence data.

The potential presence of developmental modules is studied in the urodele skull using several classical statistical methods that have not previously been used in this context. Principal component analysis (PCA) of ossification sequence data on 21 bones in 21 extant urodele species suggests the presence of up to four developmental modules, but examination of statistically significant correlations using phylogenetic independent contrasts (PIC) and correcting for multiple tests using the false discovery rate suggests the presence of only two modules of uneven size and of two bones that may not be part of these modules. Thus, PCA does not appear to be a reliable method to investigate modularity; direct investigation of statistically significant correlations using PIC or other phylogeny-informed methods is recommended. A binomial test of the distribution of significant correlations between characters shows significant heterogeneity, which suggests that modularity is indeed present in the data. A cluster analysis gives inconsistent results that apparently do not reflect developmental modules. The data include a phylogenetic signal, as shown by a permutation-based test with squared change parsimony, but this is detectable only when the whole matrix is analyzed, and a plot of the tree onto developmental space through Evolutionary PCA shows that homoplasy is pervasive. Evolutionary rates between characters vary about 90-fold. Canonical variates analyses suggest that obligatorily neotenic urodeles may be discriminated from other urodeles on the basis of cranial ossification sequence data.

Incorporating phylogenetic uncertainty on phylogeny-based palaeontological dating and the timing of turtle diversification.

Methods improving the performance of molecular dating of divergence time of clades have improved dramatically in recent years. The calibration of molecular dating using the first appearance of a clade in the fossil record is a crucial step towards inferring the minimal diversification time of various groups and the choice of extinct taxa can strongly influence the molecular dates. Here, we evaluate the uncertainty on the phylogenetic position of extinct taxa through non-parametric bootstrapping. The recognition of phylogenetic uncertainty resulted in the definition of the Bootstrap Uncertainty Range (BUR) for the age of first appearance of a given clade. The BUR is calculated as the interval of geological time in which the diversification of a given clade can be inferred to have occurred, based on the temporal information of the fossil record and the topologies of the bootstrap trees. Divergence times based on BUR analyses were calculated for three clades of turtles: Testudines, Pleurodira and Cryptodira. This resulted in extensive uncertainty ranges of topology-dependent minimal divergence dates for these clades.

Movement in a gravitational field: The question of limb interarticular coordination in terrestrial vertebrates.

In this paper, we demonstrated that interarticular coordination of terrestrial tetrapods emerges from an environment highly constrained by friction and the gravitational field. We briefly review recent works on the jumping behavior in squamates, lemurs and amphibians. We then explore previously published work as well as some unpublished experimental data on human jumping. Finally, we end by inferring locomotion in some of the first limbed vertebrates using a simulation procedure. All these data show that despite changes in shape, structure, and motor controls of taxa, the same spatio-temporal sequence of joint displacements always occurs when the movement is executed in a terrestrial environment. Comparison with aquatic locomotion argues for the hypothesis that this pattern emerged in early terrestrial tetrapods as a response to the gravitational constraint and the terrestrial frictional environment.

Correlation between the presence of a cecal appendix and reduced diarrhea severity in primates: new insights into the presumed function of the appendix.

Increased severity or recurrence risk of some specific infectious diarrhea, such a salmonellosis or Clostridium difficile colitis, have been reported after an appendectomy in human patients. While several other mammals also possess an appendix, the suspected protective function against diarrhea conferred by this structure is known only in humans. From a retrospective collection of veterinary records of 1251 primates attributed to 45 species, including 13 species with an appendix and 32 without, we identified 2855 episodes of diarrhea, 13% of which were classified as severe diarrhea requiring a therapeutic medication or associated with a fatal issue. We identified a lower risk of severe diarrhea among primate species with an appendix, especially in the early part of life when the risk of diarrhea is maximal. Moreover, we observed a delayed onset of diarrhea and of severe diarrhea in species possessing an appendix. Interestingly, none of the primates with an appendix were diagnosed, treated or died of an acute appendicitis during the 20 years of veterinarian follow-up. These results clarify the function of the appendix among primates, as protection against diarrhea. This supports its presumed function in humans and is congruent with the existence of a selective advantage conferred by this structure.

Developmental characters in phylogenetic inference and their absolute timing information.

Despite the recent surge of interest in studying the evolution of development, surprisingly little work has been done to investigate the phylogenetic signal in developmental characters. Yet, both the potential usefulness of developmental characters in phylogenetic reconstruction and the validity of inferences on the evolution of developmental characters depend on the presence of such a phylogenetic signal and on the ability of our coding scheme to capture it. In a recent study, we showed, using simulations, that a new method (called the continuous analysis) using standardized time or ontogenetic sequence data and squared-change parsimony outperformed event pairing and event cracking in analyzing developmental data on a reference phylogeny. Using the same simulated data, we demonstrate that all these coding methods (event pairing and standardized time or ontogenetic sequence data) can be used to produce phylogenetically informative data. Despite some dependence between characters (the position of an event in an ontogenetic sequence is not independent of the position of other events in the same sequence), parsimony analysis of such characters converges on the correct phylogeny as the amount of data increases. In this context, the new coding method (developed for the continuous analysis) outperforms event pairing; it recovers a lower proportion of incorrect clades. This study thus validates the use of ontogenetic data in phylogenetic inference and presents a simple coding scheme that can extract a reliable phylogenetic signal from these data.

The reconstructed evolutionary process with the fossil record.

Using the fossil record yields more detailed reconstructions of the evolutionary process than what is obtained from contemporary lineages only. In this work, we present a stochastic process modeling not only speciation and extinction, but also fossil finds. Next, we derive an explicit formula for the likelihood of a reconstructed phylogeny with fossils, which can be used to estimate the speciation and extinction rates. Finally, we provide a comparative simulation-based evaluation of the accuracy of estimations of these rates from complete phylogenies (including extinct lineages), from reconstructions with contemporary lineages only and from reconstructions with contemporary lineages and the fossil record. Results show that taking the fossil record into account yields more accurate estimates of speciation and extinction rates than considering only contemporary lineages.

Life history and ossification patterns in Miguashaia bureaui reveal the early evolution of osteogenesis in coelacanths.

The study of development is critical for revealing the evolution of major vertebrate lineages. Coelacanths have one of the longest evolutionary histories among osteichthyans, but despite access to extant representatives, the onset of their weakly ossified endoskeleton is still poorly understood. Here we present the first palaeohistological and skeletochronological study of Miguashaia bureaui from the Upper Devonian of Canada, pivotal for exploring the palaeobiology and early evolution of osteogenesis in coelacanths. Cross sections of the caudal fin bones show that the cortex is made of layers of primary bone separated by lines of arrested growth, indicative of a cyclical growth. The medullary cavity displays remnants of calcified cartilage associated with bony trabeculae, characteristic of endochondral ossification. A skeletochronological analysis indicates that rapid growth during a short juvenile period was followed by slower growth in adulthood. Our new analysis highlights the life history and palaeoecology of Miguashaia bureaui and reveals that, despite differences in size and habitat, the poor endoskeletal ossification known in the extant Latimeria chalumnae can be traced back at least 375 million years ago.