Morphological evolution and modularity of the caecilian skull

Caecilians (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.

We found highest support for a ten-module model, with greater integration of the posterior skull. Phylogenetic signal was significant (K_mult = 0.87, p < 0.01), but stronger in anterior modules, while allometric influences were also significant (R² = 0.16, p < 0.01), but stronger posteriorly. Reproductive strategy and degree of fossoriality were small but significant influences on cranial morphology (R² = 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.

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