Sequence heterochrony and the evolution of development

Integration, heterochrony, and adaptation in pedal digits of syndactylous marsupials

Background

Marsupial syndactyly is a curious morphology of the foot found in all species of diprotodontian and peramelemorph marsupials. It is traditionally defined as a condition in which digits II and III of the foot are bound by skin and are reduced. Past treatments of marsupial syndactyly have not considered the implications of this unique morphology for broader issues of digit development and evolution, and the ongoing debate regarding its phylogenetic meaning lacks a broad empirical basis. This study undertakes the first interdisciplinary characterisation of syndactyly, using variance/covariance matrix comparisons of morphometric measurements, locomotor indices, ossification sequences, and re-assessment of the largely anecdotal data on the phylogenetic distribution of tarsal/metatarsal articulations and "incipient syndactyly".

Results

Syndactylous digits have virtually identical variance/covariance matrices and display heterochronic ossification timing with respect to digits IV/V. However, this does not impact on overall locomotor adaptation patterns in the syndactylous foot as determined by analysis of locomotor predictor ratios. Reports of incipient syndactyly in some marsupial clades could not be confirmed; contrary to previous claims, syndactyly does not appear to impact on tarsal bone arrangement.

Conclusion

The results suggest that marsupial syndactyly originates from a constraint that is rooted in early digit ontogeny and results in evolution of the syndactylous digits as a highly integrated unit. Although convergent evolution appears likely, syndactyly in Diprotodontia and Peramelemorpha may occur through homologous developmental processes. We argue that the term "syndactyly" is a misnomer because the marsupial condition only superficially resembles its name-giving human soft-tissue syndactyly.

FUNCTIONAL EQUIVALENCE OF MORPHOLOGIES ENABLES MORPHOLOGICAL AND ECOLOGICAL DIVERSITY

Diversity in organismal forms among taxa is thought to reflect distinct selection pressures across environments. The central assumption underlying this expectation is that taxa experiencing similar selection have similar response to that selection. However, because selection acts on trait function, taxa similarity in selection response depends crucially on the relationship between function and morphology. Further, when a trait consists of multiple parts, changes in function in response to selection can result from modification of different parts, and adaptation to the same environment might result in functional but not morphological similarity. Here, we address the extent to which functional and morphological diversity in masticatory apparatus of soricid shrews reflects a shared ecological characteristic of their diet type. We examine the factors limiting morphological variation across shrew species by assessing the relative contribution of trait function (biomechanics of the jaw), ecology, and phylogeny to species similarity in mandibular traits. We found that species that shared diet type were functionally but not morphologically similar. The presence of multiple semi-independently varying traits enabled functional equivalence of composite foraging morphologies and resulted in variable response to selection exerted by similar diet. We show that functional equivalence of multiple morphologies enabled persistence of differences in habitat use (e.g., habitat moisture and coverage) among species that specialize on the same diet. We discuss the importance of developmental and functional integration among traits for evolutionary diversification of morphological structures that generate equivalent functions.

An integrative approach identifies developmental sequence heterochronies in freshwater basommatophoran snails

Adopting an integrative approach to the study of sequence heterochrony, we compared the timing of developmental events encompassing a mixture of developmental stages and functional traits in the embryos of 12 species of basommatophoran snails in an explicit phylogenetic framework. PARSIMOV analysis demonstrated clear functional heterochronies associated both with basal branches within the phylogeny and with terminal speciation events. A consensus of changes inferred under both accelerated transformation and delayed transformation optimizations identified four heterochronies where the direction of movement was known plus six twin heterochronies where the relative movements of the two events could not be assigned. On average, 0.5 and 0.58 events were inferred to have changed their position in the developmental sequence on internal and terminal branches of the phylogeny, respectively; these values are comparable with frequencies of sequence heterochrony reported in mammals. Directional heterochronies such as the early occurrence of body flexing in relation to the ontogeny of the eye spots, heart beat, and free swimming events occurred convergently and/or at different levels (i.e., familial, generic, and species) within the phylogeny. Such a functional approach to the study of developmental sequences has highlighted the possibility that heterochrony may have played a prominent role in the evolution of this group of invertebrates.

Cranial modularity and sequence heterochrony in mammals

Heterochrony, the temporal shifting of developmental events relative to each other, requires a degree of autonomy among those processes or structures. Modularity, the division of larger structures or processes into autonomous sets of internally integrated units, is often discussed in relation to the concept of heterochrony. However, the relationship between the developmental modules derived from studies of heterochrony and evolutionary modules, which should be of adaptive importance and relate to the genotype-phenotype map, has not been explicitly studied. I analyzed a series of sectioned and whole cleared-and-stained embryological and neonatal specimens, supplemented with published ontogenetic data, to test the hypothesis that bones within the same phenotypic modules, as determined by morphometric analysis, are developmentally integrated and will display coordinated heterochronic shifts across taxa. Modularity was analyzed in cranial bone ossification sequences of 12 therian mammals. A dataset of 12-18 developmental events was used to assess if modularity in developmental sequences corresponds to six phenotypic modules, derived from a recent morphometric analysis of cranial modularity in mammals. Kendall's tau was used to measure rank correlations, with randomization tests for significance. If modularity in developmental sequences corresponds to observed phenotypic modules, bones within a single phenotypic module should show integration of developmental timing, maintaining the same timing of ossification relative to each other, despite differences in overall ossification sequences across taxa. Analyses did not find any significant conservation of developmental timing within the six phenotypic modules, meaning that bones that are highly integrated in adult morphology are not significantly integrated in developmental timing.

Heterochrony, cannibalism, and the evolution of viviparity in Salamandra salamandra

The way in which novelties that lead to macroevolutionary events originate is a major question in evolutionary biology, and one that can be addressed using the fire salamander (Salamandra salamandra) as a model system. It is exceptional among amphibians in displaying intraspecific diversity of reproductive strategies. In S. salamandra, two distinct modes of reproduction co-occur: the common mode, ovoviviparity (females giving birth to many small larvae), and a phylogenetically derived reproductive strategy, viviparity (females producing only a few large, fully metamorphosed juveniles, which are nourished maternally). We examine the relationship between heterochronic modifications of the ontogeny and the evolution of the new reproductive mode in the fire salamander. The in vitro development of embryos of ovoviviparous and viviparous salamanders from fertilization to metamorphosis is compared, highlighting the key events that distinguish the two modes of reproduction. We identify the heterochronic events that, together with the intrauterine cannibalistic behavior, characterize the derived viviparous reproductive strategy. The ways in which evolutionary novelties can arise by modification of developmental programs can be studied in S. salamandra. Moreover, the variation in reproductive modes and the associated variation of sequences of development occur in neighboring, conspecific populations. Thus, S. salamandra is a unique biological system in which evolutionary developmental research questions can be addressed at the level of populations.

Skull ontogeny: developmental patterns of fishes conserved across major tetrapod clades

In vertebrates, the ontogeny of the bony skull forms a particularly complex part of embryonic development. Although this area used to be restricted to neontology, recent discoveries of fossil ontogenies provide an additional source of data. One of the most detailed ossification sequences is known from Permo-Carboniferous amphibians, the branchiosaurids. These temnospondyls form a near-perfect link between the piscine osteichthyans and the various clades of extant tetrapods, retaining a full complement of dermal bones in the skull. For the first time, the broader evolutionary significance of these event sequences is analyzed, focusing on the identification of sequence heterochronies. A set of 120 event pairs was analyzed by event pair cracking, which helped identify active movers. A cladistic analysis of the event pair data was also carried out, highlighting some shared patterns between widely divergent clades of tetrapods. The analyses revealed an unexpected degree of similarity between the widely divergent taxa. Most interesting is the apparently modular composition of the cranial sequence: five clusters of bones were discovered in each of which the elements form in the same time window: (1) jaw bones, (2) marginal palatal elements, (3) circumorbital bones, (4) skull roof elements, and (5) neurocranial ossifications. In the studied taxa, these "modules" have in most cases been shifted fore and back on the trajectory relative to the Amia sequence, but did not disintegrate. Such "modules" might indicate a high degree of evolutionary limitation (constraint).

Larval ectoderm, organizational homology, and the origins of evolutionary novelty

Comprehending the origin of marine invertebrate larvae remains a key domain of research for evolutionary biologists, including the repeated origin of direct developmental modes in echinoids. In order to address the latter question, we surveyed existing evidence on relationships of homology between the ectoderm territories of two closely related sea urchin species in the genus Heliocidaris that differ in their developmental mode. Additionally, we explored a recently articulated idea about homology called 'organizational homology' (Müller 2003. In: Müller GB, Newman SA, editors. Origination of organismal form: beyond the gene in developmental and evolutionary biology. Cambridge, MA: A Bradford Book, The MIT Press. p 51-69. ) in the context of this specific empirical case study. Applying the perspective of organizational homology to our experimental system of congeneric echinoids has led us to a new hypothesis concerning the ectoderm evolution in these species. The extravestibular ectoderm of the direct developer Heliocidaris erythrogramma is a novel developmental territory that arose as a fusion of the oral and aboral ectoderm territories found in indirect developing echinoids such as Heliocidaris tuberculata. This hypothesis instantiates a theoretical principle concerning the origin of developmental modules, 'integration', which has been neglected because the opposite theoretical principle, 'parcellation', is more readily observable in events such as gene duplication and divergence (Wagner 1996. Am Zool 36:36-43).

Craniofacial development in marsupial mammals: Developmental origins of evolutionary change

Biologists have long studied the evolutionary consequences of the differences in reproductive and life history strategies of marsupial and eutherian mammals. Over the past few decades, the impact of these strategies on the development of the marsupial embryo and neonate has received attention. In this review, the differences in development in the craniofacial region in marsupial and eutherian mammals will be discussed. The review will highlight differences at the organogenic and cellular levels, and discuss hypotheses for shifts in the expression of important regulatory genes. The major difference in the organogenic period is a whole-scale shift in the relative timing of central nervous system structures, in particular those of the forebrain, which are delayed in marsupials, relative to the structures of the oral-facial apparatus. Correlated with the delay in development of nervous system structures, the ossification of the bones of the neurocranium are delayed, while those of the face are accelerated. This study will also review work showing that the neural crest, which provides much of the cellular material to the facial skeleton and may also carry important patterning information, is notably accelerated in its development in marsupials. Potential consequences of these observations for hypotheses on constraint, evolutionary integration, and the existence of developmental modules is discussed. Finally, the implications of these results for hypotheses on the genetic modulation of craniofacial patterning are presented.

Status and applications of genomic resources for the gray, short-tailed opossum, Monodelphis domestica, an American marsupial model for comparative biology

Owing to its small size, favourable reproductive characteristics, and simple husbandry, the gray, short-tailed opossum, Monodelphis domestica, has become the most widely distributed and intensively utilised laboratory-bred research marsupial in the world today. This article provides an overview of the current state and future projections of genomic resources for this species and discusses the potential impact of this growing resource base on active research areas that use M. domestica as a model system. The resources discussed include: fully arrayed, bacterial artificial chromosome (BAC) libraries; an expanding linkage map; developing full-genome BAC-contig and chromosomal fluorescence in situ hybridisation maps; public websites providing access to the M. domestica whole-genome-shotgun sequence trace database and the whole-genome sequence assembly; and a new project underway to create an expressed-sequence database and microchip expression arrays for functional genomics applications. Major research areas discussed span a variety of genetic, evolutionary, physiologic, reproductive, developmental, and behavioural topics, including: comparative immunogenetics; genomic imprinting; reproductive biology; neurobiology; photobiology and carcinogenesis; genetics of lipoprotein metabolism; developmental and behavioural endocrinology; sexual differentiation and development; embryonic and fetal development; meiotic recombination; genome evolution; molecular evolution and phylogenetics; and more.

Developmental mechanisms facilitating the evolution of bills and quills

Beaks and feathers epitomize inimitable avian traits. Within individuals and across species there exists astounding diversity in the size, shape, arrangement, and colour of beaks and feathers in association with various functional adaptations. What has enabled the concomitantly divergent evolution of beaks and feathers? The common denominator may lie in their developmental programmes. As revealed through recent transplant experiments using quail and duck embryos, the developmental programme for each structure utilizes mesenchyme as a dominant source of species-specific patterning information, acts as a module of closely coupled molecular and histogenic events, and operates with a high degree of spatial and temporal plasticity. By synergizing these three features, the developmental programmes underlying beaks and feathers likely have the essential potential to react spontaneously to novel conditions and new gene functions, and as a consequence are well equipped to generate and accommodate innovative phenotypes during the course of evolution.

Molar scaling in strepsirrhine primates

We examined how maxillary molar dimensions change with body and skull size estimates among 54 species of living and subfossil strepsirrhine primates. Strepsirrhine maxillary molar areas tend to scale with negative allometry, or possibly isometry, relative to body mass. This observation supports several previous scaling analyses showing that primate molar areas scale at or slightly below geometric similarity relative to body mass. Strepsirrhine molar areas do not change relative to body mass(0.75), as predicted by the metabolic scaling hypothesis. Relative to basicranial length, maxillary molar areas tend to scale with positive allometry. Previous claims that primate molar areas scale with positive allometry relative to body mass appear to rest on the incorrect assumption that skull dimensions scale isometrically with body mass. We identified specific factors that help us to better understand these observed scaling patterns. Lorisiform and lemuriform maxillary molar scaling patterns did not differ significantly, suggesting that the two infraorders had little independent influence on strepsirrhine scaling patterns. Contrary to many previous studies of primate dental allometry, we found little evidence for significant differences in molar area scaling patterns among frugivorous, folivorous, and insectivorous groups. We were able to distinguish folivorous species from frugivorous and insectivorous taxa by comparing M1 lengths and widths. Folivores tend to have a mesiodistally elongated M1 for a given buccolingual M1 width when compared to the other two dietary groups. It has recently been shown that brain mass has a strong influence on primate dental eruption rates. We extended this comparison to relative maxillary molar sizes, but found that brain mass appears to have little influence on the size of strepsirrhine molars. Alternatively, we observed a strong correlation between the relative size of the facial skull and relative molar areas among strepsirrhines. We hypothesize that this association may be underlain by a partial sharing of the patterning of development between molar and facial skull elements.

A comparative study of embryonic development of some bird species with different patterns of postnatal growth

Some studies show that birds with high postnatal growth rates (e.g. altricial species) are characterized by a rapid early development of "supply" organs, such as digestive organs. Birds with low postnatal growth rates (e.g. precocial species) exhibit a slower early development of these organs and a more rapid early development of other "demand" organs, such as brain, muscles, skeleton and feathers. To test whether these differences can be traced back to early embryonic development and whether they can be associated with changes in developmental timing, i.e. heterochrony, we compared embryos of the precocial quail and the altricial fieldfare, two bird species with low and high postnatal growth rates, respectively. We used classical staging techniques that use developmental landmarks to categorize embryonic maturity as well as morphological measurements. These techniques were combined with immune detection of muscle specific proteins in the somites. Our data showed that the anlagen of the head, brain and eyes develop earlier in the quail than in the fieldfare in contrast to the gut which develops earlier in the fieldfare than in the quail. Our data also showed that the quail and the fieldfare displayed different rates of myotome formation in the somites which contribute to muscle formation in the limbs and thorax. We believe these observations are connected with important differences in neonatal characteristics, such as the size of the brain, eyes, organs for locomotion and digestion. This leads us to the conclusion that selection for late ontogenetic characteristics can alter early embryonic development and that growth rate is of fundamental importance for the patterning of avian embryonic development. It also appears that this comparative system offers excellent opportunities to test hypotheses about heterochrony.

Heterochrony and the evolution of poecilogony: generating larval diversity

Poecilogony is the production of more than one type of young within a single species of marine invertebrate. We chose a poecilogonous polychaete to investigate potential differences in morphogenesis among offspring that are polymorphic in dispersal potentials (planktonic, benthic) and trophic modes (planktotrophy, adelphophagy). Differences in morphogenesis occur and are strongly influenced by maternal type. Females that provide extra-embryonic nutrition (as nurse eggs; type III females) also produce offspring with an accelerated onset of juvenile traits, relative to planktotrophic offspring of females that do not provide extra-embryonic nutrition (type I females). Thus, progeny of some females appear morphologically preadapted for a benthic lifestyle. Surprisingly, differences in phenotype among offspring do not parallel offspring ecotype, as offspring with early onset of juvenile traits (III) are ecologically bimodal. Some Type III offspring eat the nurse eggs (adelphophagy), have accelerated development, and hatch as benthic juveniles. In contrast, their siblings hatch as small, planktotrophic, dispersive larvae that are morphologically similar to their type III siblings, but ecologically similar to Type I planktotrophic larvae. We propose that poecilogony evolved through sequence heterochrony in morphogenesis with accelerated onset of juvenile traits in type III offspring. In addition, we suggest that heterochrony in life-history events (hatching, metamorphosis) also occurs, thereby generating offspring that are dimorphic in both phenotype and ecotype. Over time, selection acting on different levels of ontogeny (morphogenesis vs. dispersal) may balance this polymorphism and allow poecilogony to persist.

COMPARATIVE DEVELOPMENTAL PHYSIOLOGY:An Interdisciplinary Convergence

Comparative developmental physiology spans genomics to physiological ecology and evolution. Although not a new discipline, comparative developmental physiology's position at the convergence of development, physiology and evolution gives it prominent new significance. The contributions of this discipline may be particularly influential as physiologists expand beyond genomics to a true systems synthesis, integrating molecular through organ function in multiple organ systems. This review considers how developing physiological systems are directed by genes yet respond to environment and how these characteristics both constrain and enable evolution of physiological characters. Experimental approaches and methodologies of comparative developmental physiology include studying event sequences (heterochrony and heterokairy), describing the onset and progression of physiological regulation, exploiting scaling, expanding the list of animal models, using genetic engineering, and capitalizing on new miniaturized technologies for physiological investigation down to the embryonic level. A synthesis of these approaches is likely to generate a more complete understanding of how physiological systems and, indeed, whole animals develop and how populations evolve.

Schultz’s Unruly Rule: Dental Developmental Sequences and Schedules in Small-Bodied, Folivorous Lemurs

Schultz's rule (as reconstructed by Smith) states that there is a relationship between the pattern (or relative order) of eruption of molar versus secondary (replacement) teeth and the overall pace (or absolute timing) of growth and maturation. Species with 'fast' life histories (rapid dental development, rapid growth, early sexual maturation, short life spans) are said to exhibit relatively early eruption of the molars and late eruption of the secondary replacement teeth (premolars, canines, incisors), whereas species with 'slow' life histories are said to exhibit relatively late eruption of the molars and early eruption of the secondary dentition. In a recent review, B.H. Smith noted that primates with tooth combs might violate this rule because tooth combs tend to erupt early, regardless of the pace of life history. We show that exceptions to Schultz's rule among lemurs are not limited to the relative timing of eruption of the tooth comb. Rather, among lemurs, some species with extremely accelerated dental development exhibit a pattern of eruption of molars and of secondary teeth in direct opposition to the expectations of Schultz's rule. We focus particularly on the pattern (order) and pace (absolute timing) of dental development and eruption in Avahi and Lepilemur - two relatively small, nocturnal folivores with rapid dental development. These taxa differ markedly in their eruption sequences (the premolars erupt after M2 and M3 in Lepilemur but not Avahi ). We offer an explanation for the failure of Schultz's rule to predict these differences. Schultz's rule presumes that eruption timing is dependent on the size of the jaw and that, therefore, molar crown formation and eruption will be delayed in species with slow-growing jaws. We show that a variety of processes (including developmental imbrication) allows the crowns of permanent teeth to form and to erupt into jaws that might appear to be too small to accommodate them.

Regionalization of axial skeleton in the lungfishNeoceratodus forsteri (Dipnoi)

Differentiation of the axial skeleton into distinct regions, once thought to be characteristic of the Tetrapoda, also occurs in the actinopterygian Danio rerio. In these taxa, the boundary between the cervical-thoracic regions correlates with Hoxc6 expression and morphological features such as position of the pectoral fin and associated nerves, and the absence of ribs. In the lungfish Neoceratodus, a member of the extant sister taxon to the Tetrapoda, the first vertebral element to chondrify is situated well posterior to the skull, developing from somites 6 and 7 (6/7) and associated with an enlarged cranial rib and nerves innervating the pectoral fin. Two vertebral elements develop later and more anteriorly, associated with somites 4/5 and 5/6. These three elements become incorporated into the occipital region of the skull during Neoceratodus ontogeny, until the cranial rib itself articulates to the rear of the skull. These features of early development indicate a regionalization of the Neoceratodus vertebral column: the cranial rib marks the boundary between the cervical and thoracic regions, the two more anterior vertebrae lacking ribs represent the cervical region, while somites 1-4 (cranial half), lacking any vertebral development, represent the occipital region. However, the cervical region of the vertebral column is effectively lost during ontogeny of Neoceratodus. A recognizable cervical region in the tetrapod vertebral column, as in zebrafish, suggests that cervical vertebrae are not incorporated into the skull but maintained as distinct elements of the column, representing an important shift in relative developmental timing and the influence of heterochrony in this region during the fish-tetrapod transition.

Regional dissociated heterochrony in multivariate analysis

Heterochrony, the classic framework to study ontogeny and phylogeny, in essence relies on a univariate concept of shape. Though principal component plots of multivariate shape data seem to resemble classical bivariate allometric plots, the language of heterochrony cannot be translated directly into general multivariate methodology. We simulate idealized multivariate ontogenetic trajectories and demonstrate their behavior in principal component plots in shape space and in size-shape space. The concept of "dissociation", which is conventionally regarded as a change in the relationship between shape change and size change, appears to be algebraically the same as regional dissociation - the variation of apparent heterochrony by region. Only if the trajectories of two related species lie along exactly the same path in shape space can the classic terminology of heterochrony apply so that pure dissociation of size change against shape change can be detected. We demonstrate a geometric morphometric approach to these issues using adult and subadult crania of 48 Pan paniscus and 47 P. troglodytes. On each specimen we digitized 47 landmarks and 144 semilandmarks on ridge curves and the external neurocranial surface. The relation between these two species' growth trajectories is too complex for a simple summary in terms of global heterochrony.

Comparative and phylogenetic analysis of developmental sequences

Event pairing has been proposed for the optimization of developmental sequences (event sequences) on a given phylogenetic hypothesis (cladogram) to determine instances of sequence heterochrony. Here, we show that event pairing is faulty, leading to the optimization of impossible hypothetical ancestors, the underestimation of the lengths of the developmental sequences on the tree, and the proposition of synapomorphies that are not supported by the data. When used for phylogenetic analysis, event pairing can even produce cladograms that are inconsistent with the data. These errors are caused by the fact that event pairing treats dependent features as if they were independent. We present a new method for comparative and phylogenetic analysis of developmental sequences that does not exhibit these errors. Our method applies Search-based character optimization and treats the entire developmental sequence as a single character that is then analyzed by using an edit cost function, which specifies the transformation cost between pairs of observed and unobserved character states, and dynamic programming. In other words, the developmental sequence is directly optimized on the tree. We used event pairing as an edit cost function, but others are possible.

Development of physiological regulatory systems: altering the timing of crucial events

There is currently tremendous interest in how the physiology of individual animals changes and develops during ontogeny. One of the key areas is the extent to which the timing and/or rate of physiological development is fixed within an individual and to what extent can it be altered. We propose that plasticity in the timing of the onset of a particular physiological regulatory system during an individuals development be referred to as physiological heterokairy (to clearly distinguish this phenomenon from physiological heterochrony, which is an evolutionary pattern), and we marshal evidence for three different patterns of heterokairy: 1. altering relative position in the physiological itinerary; 2. altering overall rate of development per se and; 3. a combination of 1 and 2. Using these patterns as a starting point, we develop a framework for investigating physiological heterokairy which takes cognizance of the facts that multiple components of each regulatory system could appear at different times and multiple regulatory systems could come 'on-line' at different times. We finish by placing physiological heterokairy in the wider context of its ecological and evolutionary implications and its relationship to physiological genomics and heterochrony.

The evolution of development: two portraits of skull ossification in pipoid frogs

Development creates morphology, and the study of developmental processes has repeatedly shed light on patterns of morphological evolution. However, development itself evolves as well, often concomitantly with changes in life history or in morphology. In this paper, two approaches are used to examine the evolution of skull development in pipoid frogs. Pipoids have highly unusual morphologies and life histories compared to other frogs, and their development also proves to be remarkable. First, a phylogenetic examination of skull bone ossification sequences reveals that jaw ossification occurs significantly earlier in pipoids than in other frogs; this represents a reversal to the primitive vertebrate condition. Early jaw ossification in pipoids is hypothesized to result from the absence of certain larval specializations possessed by other frogs, combined with unusual larval feeding behaviors. Second, thin-plate spline morphometric studies of ontogenetic shape change reveal important differences between pipoid skull development and that of other frogs. In the course of frog evolution, there has been a shift away from salamander-like patterns of ontogenetic shape change. The pipoids represent the culmination of this trend, and their morphologies are highly derived in numerous respects. This study represents the first detailed examination of the evolution of skull development in a diverse vertebrate clade within a phylogenetic framework. It is also the first study to examine ossification sequences across vertebrates, and the first to use thin-plate spline morphometrics to quantitatively describe ontogenetic trajectories.