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.

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.

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.

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.

Hypersalinity drives convergent bone mass increases in Miocene marine mammals from the Paratethys.

Pachyosteosclerosis-a condition that creates dense, bulky bones-often characterizes the early evolution of secondarily aquatic tetrapods like whales and dolphins1-3 but then usually fades away as swimming efficiency increases.4 Here, we document a remarkable reversal of this pattern, namely the convergent re-emergence of bone densification in Miocene seals, dolphins, and whales from the epicontinental Paratethys Sea of eastern Europe and central Asia. This phenomenon was driven by imbalanced remodeling and inhibited resorption of primary trabeculae and coincided with hypersaline conditions-the Badenian salinity crisis-that affected the Central Paratethys between 13.8 and 13.4 Ma.5 Dense bones acting as ballast would have facilitated efficient swimming in the denser and more buoyant water and hence were likely adaptive in this setting. From the Central Paratethys, pachyosteosclerosis subsequently spread eastward, where it became a defining feature of the endemic late Miocene whale assemblage.6,7.

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.

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.

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.

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.].

The relationship between genome size and metabolic rate in extant vertebrates.

Genome size has long been hypothesized to affect the metabolic rate in various groups of animals. The mechanism behind this proposed association is the nucleotypic effect, in which large nucleus and cell sizes influence cellular metabolism through surface area-to-volume ratios. Here, we provide a review of the recent literature on the relationship between genome size and metabolic rate. We also conduct an analysis using phylogenetic comparative methods and a large sample of extant vertebrates. We find no evidence that the effect of genome size improves upon models in explaining metabolic rate variation. Not surprisingly, our results show a strong positive relationship between metabolic rate and body mass, as well as a substantial difference in metabolic rate between endothermic and ectothermic vertebrates, controlling for body mass. The presence of endothermy can also explain elevated rate shifts in metabolic rate whereas genome size cannot. We further find no evidence for a punctuated model of evolution for metabolic rate. Our results do not rule out the possibility that genome size affects cellular physiology in some tissues, but they are consistent with previous research suggesting little support for a direct functional connection between genome size and basal metabolic rate in extant vertebrates. This article is part of the theme issue 'Vertebrate palaeophysiology'.

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.

Phylogeny of Paleozoic limbed vertebrates reassessed through revision and expansion of the largest published relevant data matrix.

The largest published phylogenetic analysis of early limbed vertebrates (Ruta M, Coates MI. 2007. Journal of Systematic Palaeontology 5:69-122) recovered, for example, Seymouriamorpha, Diadectomorpha and (in some trees) Caudata as paraphyletic and found the "temnospondyl hypothesis" on the origin of Lissamphibia (TH) to be more parsimonious than the "lepospondyl hypothesis" (LH)-though only, as we show, by one step. We report 4,200 misscored cells, over half of them due to typographic and similar accidental errors. Further, some characters were duplicated; some had only one described state; for one, most taxa were scored after presumed relatives. Even potentially continuous characters were unordered, the effects of ontogeny were not sufficiently taken into account, and data published after 2001 were mostly excluded. After these issues are improved-we document and justify all changes to the matrix-but no characters are added, we find (Analysis R1) much longer trees with, for example, monophyletic Caudata, Diadectomorpha and (in some trees) Seymouriamorpha; Ichthyostega either crownward or rootward of Acanthostega; and Anthracosauria either crownward or rootward of Temnospondyli. The LH is nine steps shorter than the TH (R2; constrained) and 12 steps shorter than the "polyphyly hypothesis" (PH-R3; constrained). Brachydectes (Lysorophia) is not found next to Lissamphibia; instead, a large clade that includes the adelogyrinids, urocordylid "nectrideans" and aïstopods occupies that position. As expected from the taxon/character ratio, most bootstrap values are low. Adding 56 terminal taxa to the original 102 increases the resolution (and decreases most bootstrap values). The added taxa range in completeness from complete articulated skeletons to an incomplete lower jaw. Even though the lissamphibian-like temnospondyls Gerobatrachus, Micropholis and Tungussogyrinus and the extremely peramorphic salamander Chelotriton are added, the difference between LH (R4; unconstrained) and TH (R5) rises to 10 steps, that between LH and PH (R6) to 15; the TH also requires several more regains of lost bones than the LH. Casineria, in which we tentatively identify a postbranchial lamina, emerges rather far from amniote origins in a gephyrostegid-chroniosuchian grade. Bayesian inference (Analysis EB, settings as in R4) mostly agrees with R4. High posterior probabilities are found for Lissamphibia (1.00) and the LH (0.92); however, many branches remain weakly supported, and most are short, as expected from the small character sample. We discuss phylogeny, approaches to coding, methods of phylogenetics (Bayesian inference vs. equally weighted vs. reweighted parsimony), some character complexes (e.g. preaxial/postaxial polarity in limb development), and prospects for further improvement of this matrix. Even in its revised state, the matrix cannot provide a robust assessment of the phylogeny of early limbed vertebrates. Sufficient improvement will be laborious-but not difficult.

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.

Exceptional soft tissues preservation in a mummified frog-eating Eocene salamander.

Fossils are almost always represented by hard tissues but we present here the exceptional case of a three-dimensionally preserved specimen that was 'mummified' (likely between 40 and 34 million years ago) in a terrestrial karstic environment. This fossil is the incomplete body of a salamander, Phosphotriton sigei, whose skeleton and external morphology are well preserved, as revealed by phase-contrast synchrotron X-ray microtomography. In addition, internal structures composed of soft tissues preserved in three dimensions are now identified: a lung, the spinal cord, a lumbosacral plexus, the digestive tract, muscles and urogenital organs that may be cloacal glands. These are among the oldest known cases of three-dimensional preservation of these organs in vertebrates and shed light on the ecology of this salamander. Indeed, the digestive tract contains remains of a frog, which represents the only known case of an extinct salamander that fed on a frog, an extremely rare type of predation in extant salamanders. These new data improve our scarce knowledge on soft tissue anatomy of early urodeles and should prove useful for future biologists and palaeontologists working on urodele evolutionary biology. We also suggest that the presence of bat guano and carcasses represented a close source of phosphorus, favouring preservation of soft tissues. Bone microanatomy indicates that P. sigei was likely amphibious or terrestrial, and was probably not neotenic.

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.].

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.