BackgroundCaecilians (Gymnophiona) are the least speciose extant lissamphibian order, yet living forms capture approximately 250 million years of evolution since their earliest divergences. This long history is reflected in the broad range of skull morphologies exhibited by this largely fossorial, but developmentally diverse, clade. However, this diversity of form makes quantification of caecilian cranial morphology challenging, with highly variable presence or absence of many structures. Consequently, few studies have examined morphological evolution across caecilians. This extensive variation also raises the question of degree of conservation of cranial modules (semi-autonomous subsets of highly-integrated traits) within this clade, allowing us to assess the importance of modular organisation in shaping morphological evolution. We used an intensive surface geometric morphometric approach to quantify cranial morphological variation across all 32 extant caecilian genera. We defined 16 cranial regions using 53 landmarks and 687 curve and 729 surface sliding semilandmarks. With these unprecedented high-dimensional data, we analysed cranial shape and modularity across caecilians assessing phylogenetic, allometric and ecological influences on cranial evolution, as well as investigating the relationships among integration, evolutionary rate, and morphological disparity.ResultsWe found highest support for a ten-module model, with greater integration of the posterior skull. Phylogenetic signal was significant (Kmult = 0.87, p < 0.01), but stronger in anterior modules, while allometric influences were also significant (R2 = 0.16, p < 0.01), but stronger posteriorly. Reproductive strategy and degree of fossoriality were small but significant influences on cranial morphology (R2 = 0.03–0.05), after phylogenetic (p < 0.03) and multiple-test (p < 0.05) corrections. The quadrate-squamosal ‘cheek’ module was the fastest evolving module, perhaps due to its pivotal role in the unique dual jaw-closing mechanism of caecilians. Highly integrated modules exhibited both high and low disparities, and no relationship was evident between integration and evolutionary rate.ConclusionsOur high-dimensional approach robustly characterises caecilian cranial evolution and demonstrates that caecilian crania are highly modular and that cranial modules are shaped by differential phylogenetic, allometric, and ecological effects. More broadly, and in contrast to recent studies, this work suggests that there is no simple relationship between integration and evolutionary rate or disparity.Electronic supplementary materialThe online version of this article (10.1186/s12862-018-1342-7) contains supplementary material, which is available to authorized users.

Synopsis Advances in imaging technologies, such as computed tomography (CT) and surface scanning, have facilitated the rapid generation of large datasets of high-resolution three-dimensional (3D) specimen reconstructions in recent years. The wealth of phenotypic information available from these datasets has the potential to inform our understanding of morphological variation and evolution. However, the ever-increasing ease of compiling 3D datasets has created an urgent need for sophisticated methods of capturing high-density shape data that reflect the biological complexity in form. Landmarks often do not take full advantage of the rich shape information available from high-resolution 3D specimen reconstructions, as they are typically restricted to sutures or processes that can be reliably identified across specimens and exclude most of the surface morphology. The development of sliding and surface semilandmark techniques has greatly enhanced the quantification of shape, but their application to diverse datasets can be challenging, especially when dealing with the variable absence of some regions within a structure. Using comprehensive 3D datasets of crania that span the entire clades of birds, squamates and caecilians, we demonstrate methods for capturing morphology across incredibly diverse shapes. We detail many of the difficulties associated with applying semilandmarks to comparable regions across highly disparate structures, and provide solutions to some of these challenges, while considering the consequences of decisions one makes in applying these approaches. Finally, we analyze the benefits of high-density sliding semilandmark approaches over landmark-only studies for capturing shape across diverse organisms and discuss the promise of these approaches for the study of organismal form.

The field of comparative morphology has entered a new phase with the rapid generation of high-resolution three-dimensional (3D) data. With freely available 3D data of thousands of species, methods for quantifying morphology that harness this rich phenotypic information are quickly emerging. Among these techniques, high-density geometric morphometric approaches provide a powerful and versatile framework to robustly characterize shape and phenotypic integration, the covariances among morphological traits. These methods are particularly useful for analyses of complex structures and across disparate taxa, which may share few landmarks of unambiguous homology. However, high-density geometric morphometrics also brings challenges, for example, with statistical, but not biological, covariances imposed by placement and sliding of semilandmarks and registration methods such as Procrustes superimposition. Here, we present simulations and case studies of high-density datasets for squamates, birds, and caecilians that exemplify the promise and challenges of high-dimensional analyses of phenotypic integration and modularity. We assess: (1) the relative merits of “big” high-density geometric morphometrics data over traditional shape data; (2) the impact of Procrustes superimposition on analyses of integration and modularity; and (3) differences in patterns of integration between analyses using high-density geometric morphometrics and those using discrete landmarks. We demonstrate that for many skull regions, 20–30 landmarks and/or semilandmarks are needed to accurately characterize their shape variation, and landmark-only analyses do a particularly poor job of capturing shape variation in vault and rostrum bones. Procrustes superimposition can mask modularity, especially when landmarks covary in parallel directions, but this effect decreases with more biologically complex covariance patterns. The directional effect of landmark variation on the position of the centroid affects recovery of covariance patterns more than landmark number does. Landmark-only and landmark-plus-sliding-semilandmark analyses of integration are generally congruent in overall pattern of integration, but landmark-only analyses tend to show higher integration between adjacent bones, especially when landmarks placed on the sutures between bones introduces a boundary bias. Allometry may be a stronger influence on patterns of integration in landmark-only analyses, which show stronger integration prior to removal of allometric effects compared to analyses including semilandmarks. High-density geometric morphometrics has its challenges and drawbacks, but our analyses of simulated and empirical datasets demonstrate that these potential issues are unlikely to obscure genuine biological signal. Rather, high-density geometric morphometric data exceed traditional landmark-based methods in characterization of morphology and allow more nuanced comparisons across disparate taxa. Combined with the rapid increases in 3D data availability, high-density morphometric approaches have immense potential to propel a new class of studies of comparative morphology and phenotypic integration.

Metamorphosis is widespread across the animal kingdom and induces fundamental changes in the morphology, habitat, and resources used by an organism during its lifetime. Metamorphic species are likely to experience more dynamic selective pressures through ontogeny compared to those with single-phase life cycles, which may drive divergent evolutionary dynamics. Here, we reconstruct the cranial evolution of the salamander using geometric morphometric data from 148 species spanning their full phylogenetic, developmental, and ecological diversity. We demonstrate that life cycle influences cranial shape diversity and rate of evolution. Shifts in rate of cranial evolution are consistently associated with transitions from biphasic to either direct-developing or paedomorphic life cycle strategies. Direct-developers exhibit the slowest rates of evolution and lowest disparity, and paedomorphic species the highest. Species undergoing complete metamorphosis (biphasic and direct-developing) exhibit greater cranial modularity (evolutionary independence among regions) than do paedomorphic species, which undergo differential metamorphosis. Biphasic and direct-developing species also display elevated disparity relative to evolutionary rate for bones associated with feeding, whereas this is not the case for paedomorphic species. Metamorphosis has profoundly influenced salamander cranial evolution, requiring greater autonomy of cranial elements and facilitating the rapid evolution of regions that are remodelled through ontogeny. Rather than compounding functional constraints on variation, metamorphosis appears to have promoted salamander morphological evolution over 180 million years, which may explain the ubiquity of this complex life cycle strategy across disparate organisms. MainDevelopmental processes play a fundamental role in structuring the morphological diversity of organisms 1-3 , being both a driver of, and a constraint on, phenotypic change 1,4,5 . As such, shifts in development and life history can have profound impacts on the evolutionary trajectories of lineages. Early attempts to delineate these effects resulted in the recapitulationist doctrine, stating that ontogeny replicates phylogeny 6 . However, studies of groups such as amphibians have shown that shifts in developmental strategies (e.g., biphasic, direct-development, paedomorphy, and viviparity) have occurred many times. In some cases, metamorphic species can even eliminate later ontogenetic stages and mature with larval