The utility of cranial ontogeny for phylogenetic inference: a case study in crocodylians using geometric morphometrics

The degree to which the ontogeny of organisms could facilitate our understanding of phylogenetic relationships has long been a subject of contention in evolutionary biology. The famed notion that ‘ontogeny recapitulates phylogeny’ has been largely discredited, but there remains an expectation that closely related organisms undergo similar morphological transformations throughout ontogeny. To test this assumption, we used three‐dimensional geometric morphometric methods to characterize the cranial morphology of 10 extant crocodylian species and construct allometric trajectories that model the post‐natal ontogenetic shape changes. Using time‐calibrated molecular and morphological trees, we employed a suite of comparative phylogenetic methods to assess the extent of phylogenetic signal in these trajectories. All analyses largely demonstrated a lack of significant phylogenetic signal, indicating that ontogenetic shape changes contain little phylogenetic information. Notably, some Mantel tests yielded marginally significant results when analysed with the morphological tree, which suggest that the underlying signal in these trajectories is correlated with similarities in the adult cranial morphology. However, despite these instances, all other analyses, including more powerful tests for phylogenetic signal, recovered statistical and visual evidence against the assumption that similarities in ontogenetic shape changes are commensurate with phylogenetic relatedness and thus bring into question the efficacy of using allometric trajectories for phylogenetic inference.     

How to identify (as opposed to define) a homoplasy: examples from fossil and living great apes

There is much debate on the definitions of homoplasy and homology, and on how to spot them among character states used in a phylogenetic analysis. Many advocate what I call a "processual approach," in which information on genetics, development, function, or other criteria help a priori in identifying two character states as homologous or homoplastic. I argue that the processes represented by these criteria are insufficiently known for most organisms and most characters to be reliably used to identify homoplasies and homologies. Instead, while not foolproof, phylogeny should be the ultimate test for homology. Character states are assumed to be homologous a priori because this is falsifiable and because their initial inclusion in the character-state analysis is based on the assumption that they may be phylogenetically informative. If they fall out as symplesiomorphies or synapomorphies in a phylogenetic analysis, their status as homologies remains unfalsified. If they fall out as homoplasies, having evolved independently in more than one clade, their status as homologous is falsified, and a homoplasy is identified. The character-state transformation series, functional morphology, finer levels of morphological comparison, and the distribution and correlation of characters all help to explain the presence of homoplasies in a given phylogeny. Explaining these homoplasies, and not ignoring them as "noise," should be as much a goal of phylogenetic analysis as the production of a phylogeny. Examples from the fossil record of Miocene hominoids are given to illustrate the advantages of a process-informs-pattern-recognition-after-the-fact approach to understanding the evolution of character states.     

Parallel evolution of larval morphology and habitat in the snail-killing fly genus Tetanocera

In this study, we sequenced one nuclear and three mitochondrial DNA loci to construct a robust estimate of phylogeny for all available species of Tetanocera. Character optimizations suggested that aquatic habitat was the ancestral condition for Tetanocera larvae, and that there were at least three parallel transitions to terrestrial habitat, with one reversal. Maximum likelihood analyses of character state transformations showed significant correlations between habitat transitions and changes in four larval morphological characteristics (cuticular pigmentation and three characters associated with the posterior spiracular disc). We provide evidence that phylogenetic niche conservatism has been responsible for the maintenance of aquatic-associated larval morphological character states, and that concerted convergence and/or gene linkage was responsible for parallel morphological changes that were derived in conjunction with habitat transitions. These habitat-morphology associations were consistent with the action of natural selection in facilitating the morphological changes that occurred during parallel aquatic to terrestrial habitat transitions in Tetanocera.     

Quantification of homoplasy for nucleotide transitions and transversions and a reexamination of assumptions in weighted phylogenetic analysis

Nucleotide transitions are frequently down-weighted relative to transversions in phylogenetic analysis. This is based on the assumption that transitions, by virtue of their greater evolutionary rate, exhibit relatively more homoplasy and are therefore less reliable phylogenetic characters. Relative amounts of homoplastic and consistent transition and transversion changes in mitochondrial protein coding genes were determined from character-state reconstructions on a highly corroborated phylogeny of mammals. We found that although homoplasy was related to evolutionary rates and was greater for transitions, the absolute number of consistent transitions greatly exceeded the number of consistent transversions. Consequently, transitions provided substantially more useful phylogenetic information than transversions. These results suggest that down-weighting transitions may be unwarranted in many cases. This conclusion was supported by the fact that a range of transition: transversion weighting schemes applied to various mitochondrial genes and genomic partitions rarely provided improvement in phylogenetic estimates relative to equal weighting, and in some cases weighting transitions more heavily than transversions was most effective.     

Phylogenetics of the australasian mudfishes: evolution of an eel-like body plan

The mudfish genus Neochanna (Osmeriformes: Galaxiidae) contains six species that exhibit varying degrees of morphological and ecological adaptation to life in swampy conditions. Here, we present the first molecular phylogenetic analysis (16S rRNA+cytochrome b; 1681bp) of the entire genus to (1) test for monophyly of Australian and New Zealand taxa and (2) elucidate morphological character evolution. In addition, we analyse a matrix of 21 morphological characters to test for congruence between mitochondrial DNA and morphology, and to examine total evidence under a Bayesian framework. Molecular data indicate that the diadromous Tasmanian mudfish, N. cleaveri, is sister to a clade of five non-diadromous New Zealand mudfishes. Mapping of morphological characters onto the molecular phylogeny suggests an evolutionary transition from a plesiomorphic "stream" galaxiid morphotype to a more specialised "anguilliform" galaxiid morphotype. Pelvic fins have become increasingly reduced and dorsal, anal, and caudal fins are increasingly confluent. Associated with these changes are elongated nostrils, reduced eyes, and increased anterior cranial ossification. Morphological and total evidence analyses yield similar or identical topologies, respectively. The phylogenetic distribution of diadromy in Neochanna is consistent with a single loss of this character state in New Zealand. However, the strong sister relationship (3.6% divergent; 100% bootstrap support) detected between non-diadromous N. burrowsius (South Island, NZ) and N. rekohua (Chatham Islands) indicates geologically recent dispersal across 850km of ocean. Diadromy may therefore have been retained in the common ancestor of these two mudfish species, and subsequently lost from both lineages.     

Does nasal echolocation influence the modularity of the mammal skull?

In vertebrates, changes in cranial modularity can evolve rapidly in response to selection. However, mammals have apparently maintained their pattern of cranial integration throughout their evolutionary history and across tremendous morphological and ecological diversity. Here, we use phylogenetic, geometric morphometric and comparative analyses to test the hypothesis that the modularity of the mammalian skull has been remodelled in rhinolophid bats due to the novel and critical function of the nasal cavity in echolocation. We predicted that nasal echolocation has resulted in the evolution of a third cranial module, the ‘nasal dome’, in addition to the braincase and rostrum modules, which are conserved across mammals. We also test for similarities in the evolution of skull shape in relation to habitat across rhinolophids. We find that, despite broad variation in the shape of the nasal dome, the integration of the rhinolophid skull is highly consistent with conserved patterns of modularity found in other mammals. Across their broad geographical distribution, cranial shape in rhinolophids follows two major divisions that could reflect adaptations to dietary and environmental differences in African versus South Asian distributions. Our results highlight the potential of a relatively simple modular template to generate broad morphological and functional variation in mammals.     

Morphology and biogeography of Apiaceae subfamily Saniculoideae as inferred by phylogenetic analysis of molecular data

The phylogenetic placements of several African endemic genera at the base of Apiaceae subfamilies Saniculoideae and Apioideae have revolutionized ideas of relationships that affect hypotheses of character evolution and biogeography. Using an explicit phylogeny of subfamily Saniculoideae, we reconstructed the evolutionary history of phenotypic characters traditionally important in classification, identified those characters most useful in supporting relationships, and inferred historical biogeography. The 23 characters examined include those of life history, vegetative morphology, inflorescences, and fruit morphology and anatomy. These characters were optimized over trees derived from maximum parsimony analysis of chloroplast DNA trnQ-trnK sequences from 94 accessions of Apiaceae. The results revealed that many of these characters have undergone considerable modification and that traditional assumptions regarding character-state polarity are often incorrect. Infrasubfamilial relationships inferred by molecular data are supported by one to five morphological characters. However, none of these morphological characters support the monophyly of subfamilies Saniculoideae or Apioideae, the clade of Petagnaea, Eryngium and Sanicula, or the sister-group relationship between Eryngium and Sanicula . Southern African origins of Saniculoideae and of its tribes Steganotaenieae and Saniculeae are supported based on dispersal-vicariance analysis.     

THE RECONSTRUCTION OF ANCESTRAL CHARACTER STATES

The problem of error in the phylogenetic reconstruction of ancestral character states is explored by developing the model of Frumhoff and Reeve (1994). Information about the evolutionary rate of change within a character is inferred from the distribution of its character states on a known phylogeny, and this information is used to impose confidence limits on the error associated with ancestral state inference. Ancestral state inference is found to be remarkably robust under the model assumptions for a wide range of parameter values; however, the probability of error increases when the number of species within a clade is small and/or state-transition probabilities are strongly skewed in favor of the non-ancestral state. The rationale for expecting such a skew, a hypothesis of parallelism, is shown to rely on assumptions of low rates of change in at least two phylogenetically inherited characters: the tendency to occupy a particular ecological niche and the tendency to respond in a particular way to selection. A means for judging the relative likelihoods of parallelism vs. straightforward homology as explanations for a given character-state distribution is suggested. General problems with the model are discussed, as are methods for making it more realistic.     

Phylogeny, diet, and cranial integration in australodelphian marsupials

Studies of morphological integration provide valuable information on the correlated evolution of traits and its relationship to long-term patterns of morphological evolution. Thus far, studies of morphological integration in mammals have focused on placentals and have demonstrated that similarity in integration is broadly correlated with phylogenetic distance and dietary similarity. Detailed studies have also demonstrated a significant correlation between developmental relationships among structures and adult morphological integration. However, these studies have not yet been applied to marsupial taxa, which differ greatly from placentals in reproductive strategy and cranial development and could provide the diversity necessary to assess the relationships among phylogeny, ecology, development, and cranial integration. This study presents analyses of morphological integration in 20 species of australodelphian marsupials, and shows that phylogeny is significantly correlated with similarity of morphological integration in most clades. Size-related correlations have a significant affect on results, particularly in Peramelia, which shows a striking decrease in similarity of integration among species when size is removed. Diet is not significantly correlated with similarity of integration in any marsupial clade. These results show that marsupials differ markedly from placental mammals in the relationships of cranial integration, phylogeny, and diet, which may be related to the accelerated development of the masticatory apparatus in marsupials.     

Modelling rate distributions using character compatibility: implications for morphological evolution among fossil invertebrates

Rate distributions are important considerations when testing hypotheses about morphological evolution or phylogeny. They also have implications about general processes underlying character evolution. Molecular systematists often assume that rates are Poisson processes with gamma distributions. However, morphological change is the product of multiple probabilistic processes and should theoretically be affected by hierarchical integration of characters. Both factors predict lognormal rate distributions. Here, a simple inverse modelling approach assesses the best single-rate, gamma and lognormal models given observed character compatibility for 115 invertebrate groups. Tests reject the single-rate model for nearly all cases. Moreover, the lognormal outperforms the gamma for character change rates and (especially) state derivation rates. The latter in particular is consistent with integration affecting morphological character evolution.     

Head capsule characters in the Hymenoptera and their phylogenetic implications

The head capsule of a taxon sample of three outgroup and 86 ingroup taxa is examined for characters of possible phylogenetic significance within Hymenoptera. 21 morphological characters are illustrated and scored, and their character evolution explored by mapping them onto a phylogeny recently produced from a large morphological data set. Many of the characters are informative and display unambiguous changes. Most of the character support demonstrated is supportive at the superfamily or family level. In contrast, only few characters corroborate deeper nodes in the phylogeny of Hymenoptera.     

The quality of the fossil record and the accuracy of phylogenetic inferences about sampling and diversity

Because phylogenies can be estimated without stratigraphic data and because estimated phylogenies also infer gaps in sampling, some workers have used phylogeny estimates as templates for evaluating sampling from the fossil record and for "correcting" historical diversity patterns. However, it is not known how sampling intensity (the probability of sampling taxa per unit time) and completeness (the proportion of taxa sampled) affect the accuracy of phylogenetic inferences, nor how phylogenetically inferred estimates of sampling and diversity respond to inaccurate estimates of phylogeny. Both issues are addressed with a series of simulations using simple models of character evolution, varying speciation patterns, and various rates of speciation, extinction, character change, and preservation. Parsimony estimates of simulated phylogenies become less accurate as sampling decreases, and inaccurate trees chronically underestimate sampling. Biotic factors such as rates of morphologic change and extinction both affect the accuracy of phylogenetic estimates and thus affect estimated gaps in sampling, indicating that differences in implied sampling need not reflect actual differences in sampling. Errors in inferred diversity are concentrated early in the history of a clade. This, coupled with failure to account for true extinction times (i.e., the Signor-Lipps effect), inflates relative diversity levels early in clade histories. Because factors other than differences in sampling predict differences in the numbers of gaps implied by phylogeny estimates, inferred phylogenies can be misleading templates for evaluating sampling or historical diversity patterns.     

The phylogeny of land plants inferred from 18S rDNA sequences: pushing the limits of rDNA signal?

Previous studies of the phylogeny of land plants based on analysis of 18S ribosomal DNA (rDNA) sequences have generally found weak support for the relationships recovered and at least some obviously spurious relationships, resulting in equivocal inferences of land plant phylogeny. We hypothesized that greater sampling of both characters and taxa would improve inferences of land plant phylogeny based on 18S rDNA sequences. We therefore conducted a phylogenetic analysis of complete (or nearly complete) 18S rDNA sequences for 93 species of land plants and 7 green algal relatives. Parsimony analyses with equal weighting of characters and characters state changes and parsimony analyses weighting (1) stem bases half as much as loop bases and (2) transitions half as much as transversions did not produce substantially different topologies. Although the general structure of the shortest trees is consistent with most hypotheses of land plant phylogeny, several relationships, particularly among major groups of land plants, appear spurious. Increased character and taxon sampling did not substantially improve the performance of 18S rDNA in phylogenetic analyses of land plants, nor did analyses designed to accommodate variation in evolutionary rates among sites. The rate and pattern of 18S rDNA evolution across land plants may limit the usefulness of this gene for phylogeny reconstruction at deep levels of plant phylogeny. We conclude that the mosaic structure of 18S rDNA, consisting of highly conserved and highly variable regions, may contain historical signal at two levels. Rapidly evolving regions are informative for relatively recent divergences (e.g., within angiosperms, seed plants, and ferns), but homoplasy at these sites makes it difficult to resolve relationships among these groups. At deeper levels, changes in the highly conserved regions of small-subunit rDNAs provide signal across all of life. Because constraints imposed by the secondary structure of the rRNA may affect the phylogenetic information content of 18S rDNA, we suggest that 18S rDNA sequences be combined with other data and that methods of analysis be employed to accommodate these differences in evolutionary patterns, particularly across deep divergences in the tree of life.     

Convergence and coevolution in a mutualism: evidence from a molecular phylogeny of Ficus

Abstract.— The interaction between Ficus (Moraceae) and their pollinating wasps (Chalcidoidea: Agaonidae; more than 700 species-specific couples) is one of the most specialized mutualisms found in nature. Both partners of this interaction show extensive variation in their respective biology. Here we investigate Ficus life-history trait evolution and fig/fig wasp coadaptation in the context of a well-resolved molecular phylogeny. Mapping out variations in Ficus life-history traits on an independently derived phylogeny constructed from ribosomal DNA sequences (external and internal transcribed spacer) reveals several parallel transitions in Ficus growth habit and breeding system. Convergent trait evolution might explain the discrepancies between morphological analyses and our molecular reconstruction of the genus. Morphological characters probably correlate with growth habit and breeding system and could therefore be subject to convergent evolution. Furthermore, we reconstruct the evolution of Ficus inflorescence characters that are considered adaptations to the pollinators. Our phylogeny reveals convergences in ostiole shape, stigma morphology, and stamen:ovule ratio. Statistical tests taking into account the phylogenetic relationship of the species show that transitions in ostiole shape are correlated with variation in wasp pollinator head shape, and evolutionary changes in stigma morphology and stamen:ovule ratio correlate with changes in the pollination behavior of the associated wasp. These correlations provide evidence for reciprocal adaptations of morphological characters between these mutualistic partners that have interacted over a long evolutionary time. In light of previous ecological studies on mutualism, we discuss the adaptive significance of these correlations and what they can tell us about the coevolutionary process occurring between figs and their pollinators.     

Insights into the ecology and evolutionary success of crocodilians revealed through bite-force and tooth-pressure experimentation

Background

Crocodilians have dominated predatory niches at the water-land interface for over 85 million years. Like their ancestors, living species show substantial variation in their jaw proportions, dental form and body size. These differences are often assumed to reflect anatomical specialization related to feeding and niche occupation, but quantified data are scant. How these factors relate to biomechanical performance during feeding and their relevance to crocodilian evolutionary success are not known.

Methodology/Principal Findings

We measured adult bite forces and tooth pressures in all 23 extant crocodilian species and analyzed the results in ecological and phylogenetic contexts. We demonstrate that these reptiles generate the highest bite forces and tooth pressures known for any living animals. Bite forces strongly correlate with body size, and size changes are a major mechanism of feeding evolution in this group. Jaw shape demonstrates surprisingly little correlation to bite force and pressures. Bite forces can now be predicted in fossil crocodilians using the regression equations generated in this research.

Conclusions/Significance

Critical to crocodilian long-term success was the evolution of a high bite-force generating musculo-skeletal architecture. Once achieved, the relative force capacities of this system went essentially unmodified throughout subsequent diversification. Rampant changes in body size and concurrent changes in bite force served as a mechanism to allow access to differing prey types and sizes. Further access to the diversity of near-shore prey was gained primarily through changes in tooth pressure via the evolution of dental form and distributions of the teeth within the jaws. Rostral proportions changed substantially throughout crocodilian evolution, but not in correspondence with bite forces. The biomechanical and ecological ramifications of such changes need further examination.     

Stretch Coding and Block Coding: Two New Strategies to Represent Questionably Aligned DNA Sequences

Most coding strategies that address the problem of questionable alignment (elision, case sensitive, missing, polymorphic, gaps as presence/absence matrix) conflict with phylogenetic principles, particularly those relating to the concept of homology (shared similiarity explained by common ancestry). In some cases, the test of conjunction is failed. In other cases, characters that are coded ambiguously can lead to character-state optimization in the terminal taxa that conflicts with the original observations. Only data exclusion and contraction avoid these pitfalls. In highly dissimilar sequences additional character states can represent the available information. Two new methods that accomplish this—block and stretch coding—are introduced here. These two new coding strategies are not in conflict with the test of conjunction and do not contradict the original observations. They are comparable to coding practices with morphological data once the intrinsic differences due to character-state identity and topographical identity have been taken into account. It is suggested that, of the three recoding methods, the one is selected that preserves the maximum potential phylogenetic information as measured with the minimum number of steps required for the particular part of the data matrix. Received: 1 August 2000 / Accepted: 10 July 2001     

Constructional morphology within the head of hammerhead sharks (sphyrnidae)

The study of functional trade‐offs is important if a structure, such as the cranium, serves multiple biological roles, and is, therefore, shaped by multiple selective pressures. The sphyrnid cephalofoil presents an excellent model for investigating potential trade‐offs among sensory, neural, and feeding structures. In this study, hammerhead shark species were chosen to represent differences in head form through phylogeny. A combination of surface‐based geometric morphometrics, computed tomography (CT) volumetric analysis, and phylogenetic analyses were utilized to investigate potential trade‐offs within the head. Hammerhead sharks display a diversity of cranial morphologies where the position of the eyes and nares vary among species, with only minor changes in shape, position, and volume of the feeding apparatus through phylogeny. The basal winghead shark, Eusphyra blochii, has small anteriorly positioned eyes. Through phylogeny, the relative size and position of the eyes change, such that derived species have larger, more medially positioned eyes. The lateral position of the external nares is highly variable, showing no phylogenetic trend. Mouth size and position are conserved, remaining relatively unchanged. Volumetric CT analyses reveal no trade‐offs between the feeding apparatus and the remaining cranial structures. The few trade‐offs were isolated to the nasal capsule volume's inverse correlation with braincase, chondrocranial, and total cephalofoil volume. Eye volume also decreased as cephalofoil width increased. These data indicate that despite considerable changes in head shape, much of the head is morphologically conserved through sphyrnid phylogeny, particularly the jaw cartilages and their associated feeding muscles, with shape change and morphological trade‐offs being primarily confined to the lateral wings of the cephalofoil and their associated sensory structures. J. Morphol. 276:526–539, 2015. © 2015 Wiley Periodicals, Inc.     

Tooth and cranial disparity in the fossil relatives of Sphenodon (Rhynchocephalia) dispute the persistent ‘living fossil’ label

The tuatara (Sphenodon punctatus) is the only living representative of Rhynchocephalia, a group of small vertebrates that originated about 250 million years ago. The tuatara has been referred to as a living fossil; however, the group to which it belongs included a much greater diversity of forms in the Mesozoic. We explore the morphological diversity of Rhynchocephalia and stem lepidosaur relatives (Sphenodon plus 13 fossil relatives) by employing a combination of geometric morphometrics and comparative methods. Geometric morphometrics is used to explore cranium size and shape at interspecific scale, while comparative methods are employed to test association between skull shape and size and tooth number after taking phylogeny into account. Two phylogenetic topologies have been considered to generate a phylomorphospace and quantify the phylogenetic signal in skull shape data, the ancestral state reconstruction as well as morphological disparity using disparity through time plots (DTT). Rhynchocephalia exhibit a significant phylogenetic signal in skull shape that compares well with that computed for other extinct vertebrate groups. A consistent form of allometry has little impact on skull shape evolution while the number of teeth significantly correlates with skull shape also after taking phylogeny into account. The ancestral state reconstruction demonstrates a dramatic shape difference between the skull of Sphenodon and its much larger Cretaceous relative Priosphenodon. Additionally, DTT demonstrates that skull shape disparity is higher between rather than within clades while the opposite applies to skull size and number of teeth. These results were not altered by the use of competing phylogenic hypotheses. Rhynchocephalia evolved as a morphologically diverse group with a dramatic radiation in the Late Triassic and Early Jurassic about 200 million years ago. Differences in size are not marked between species whereas changes in number of teeth are associated with co‐ordinated shape changes in the skull to accommodate larger masticatory muscles. These results show that the tuatara is not the product of evolutionary stasis but that it represents the only survivor of a diverse Mesozoic radiation whose subsequent decline remains to be explained.