Ecomorphological diversification in squamates from conserved pattern of cranial integration

Evolutionary bursts drive morphological novelty in the world's largest skinks

Animal phenotypes evolve and diverge as a result of differing selective pressures and drift. These processes leave unique signatures in patterns of trait evolution, impacting the tempo and mode of morphological macroevolution. While there is a broad understanding of the history of some organismal traits (e.g., body size), there is little consensus about the evolutionary mode of most others. This includes the relative contribution of prolonged (Darwinian gradualist) and episodic (Simpsonian jump) changes toward the evolution of novel morphologies. Here, we use new exon-capture and linear morphological datasets to investigate the tempo and mode of morphological evolution in Australo-Melanesian Tiliquini skinks. We generate a well-supported time-calibrated phylogenomic tree from ∼400 nuclear markers for more than 100 specimens, including undescribed diversity, and provide unprecedented resolution of the rapid Miocene diversification of these lizards. By collecting a morphological dataset that encompasses the lizard body plan (19 traits across the head, body, limb, and tail), we are able to identify that most traits evolve conservatively, but infrequent evolutionary bursts result in morphological novelty. These phenotypic discontinuities occur via rapid rate increases along individual branches, inconsistent with both gradualistic and punctuated equilibrial evolutionary modes. Instead, this "punctuated gradualism" has resulted in the rapid evolution of blue-tongued giants and armored dwarves in the ∼20 million years since colonizing Australia. These results outline the evolutionary pathway toward new morphologies and highlight the heterogeneity of evolutionary tempo and mode, even within individual traits.

The skull of the Turks and Caicos rock iguana, Cyclura carinata (Squamata: Iguanidae)

We provide a detailed and first description of the skull, hyoid apparatus, and trachea of the Turks and Caicos rock iguana, Cyclura carinata (Squamata: Iguanidae). Cyclura is a radiation of iguanas restricted to islands of the Caribbean Sea. Species of Cyclura have high rates of endemism, and all species are severely threatened with extinction. Our anatomical description of this threatened iguana is based on high-resolution computed tomography scans of one adult, one putative adult or near adult, and one juvenile specimen, and includes three-dimensional segmented renderings and visualizations. We discuss some observations of intraspecific and ontogenetic variation, and provide a brief comparison with specimens of another species of Cyclura and published descriptions of other iguanas. Our study provides a cranial osteological framework for Cyclura and augments the body of knowledge on iguana anatomy generally. Finally, we posit that our description and future studies may facilitate identification of fossil Cyclura, which could help understand the paleobiogeography of the genus.

Compound osteoderms preserved in amber reveal the oldest known skink

Scincidae is one of the most species-rich and cosmopolitan clades of squamate reptiles. Abundant disarticulated fossil material has also been attributed to this group, however, no complete pre-Cenozoic crown-scincid specimens have been found. A specimen in Burmite (99 MYA) is the first fossil that can be unambiguously referred to this clade. Our analyses place it as nested within extant skinks, supported by the presence of compound osteoderms formed by articulated small ostedermites. The specimen has a combination of dorsal and ventral compound osteoderms and overlapping cycloid scales that is limited to skinks. We propose that this type of osteoderm evolved as a response to an increased overlap of scales, and to reduced stiffness of the dermal armour. Compound osteoderms could be a key innovation that facilitated diversification in this megadiverse family.

The skull of the semi-aquatic neotropical lizard Echinosaura horrida (Gymnophthalmidae: Cercosaurinae) reveals new synapomorphies within Gymnophthalmoidea

The rough teiid or water cork lizard (Echinosaura horrida) is a small reptile from Colombia and Ecuador placed in a genus that contains eight species and well-known phylogenetic relationships. Here we provide a detailed description and illustrations, bone by bone, of its skull, while we discussed its intraspecific variation by comparing high-resolution computed tomography data from two specimens and the variation within the genus by including previously published data from Echinosaura fischerorum. This allowed to propose putative diagnostic character states for Echinosaura horrida and synapomorphies for Echinosaura. In addition, our discussion includes broader comparisons of new character transformations of the jugal, vomer, orbitosphenoid, and hyoid. These characters are important for diagnosing clades at different levels of the Gymnophthalmoidea phylogeny.

Evolutionary patterns of cat-like carnivorans unveil drivers of the sabertooth morphology

The sabertooth morphology stands as a classic case of convergence, manifesting recurrently across various vertebrate groups, prominently within two carnivorans clades: felids and nimravids. Nonetheless, the evolutionary mechanisms driving these recurring phenotypes remain insufficiently understood, lacking a robust phylogenetic and spatiotemporal framework. We reconstruct the tempo and mode of craniomandibular evolution of Felidae and Nimravidae and evaluate the strength of the dichotomy between conical and saber-toothed species, as well as within saber-toothed morphotypes. To do so, we investigate morphological variation, convergence, phenotypic integration, and evolutionary rates, employing a comprehensive dataset of nearly 200 3D models encompassing mandibles and crania from both extinct and extant feline-like carnivorans, spanning their entire evolutionary timeline. Our results reject the hypothesis of a distinctive sabertooth morphology, revealing instead a continuous spectrum of feline-like phenotypes in both the cranium and mandible, with sporadic instances of unequivocal convergence. Disparity peaked at the end of the Miocene and is usually higher in clades containing taxa with extreme sabertoothed adaptations. We show that taxa with saberteeth exhibit a lower degree of craniomandibular integration, allowing to exhibit a greater range of phenotypes. Those same groups usually show a burst of morphological evolutionary rate at the beginning of their evolutionary history. Consequently, we propose that a reduced degree of integration coupled with rapid evolutionary rates emerge as key components in the development of a sabertooth morphology in multiple clades.

Dynamic evolutionary interplay between ontogenetic skull patterning and whole-head integration

The arrangement and morphology of the vertebrate skull reflect functional and ecological demands, making it a highly adaptable structure. However, the fundamental developmental and macroevolutionary mechanisms leading to different vertebrate skull phenotypes remain unclear. Here we exploit the morphological diversity of squamate reptiles to assess the developmental and evolutionary patterns of skull variation and covariation in the whole head. Our geometric morphometric analysis of a complex squamate ontogenetic dataset (209 specimens, 169 embryos, 44 species), covering stages from craniofacial primordia to fully ossified bones, reveals that morphological differences between snake and lizard skulls arose gradually through changes in spatial relationships (heterotopy) followed by alterations in developmental timing or rate (heterochrony). Along with dynamic spatiotemporal changes in the integration pattern of skull bone shape and topology with surrounding brain tissues and sensory organs, we identify a relatively higher phenotypic integration of the developing snake head compared with lizards. The eye, nasal cavity and Jacobson's organ are pivotal in skull morphogenesis, highlighting the importance of sensory rearrangements in snake evolution. Furthermore, our findings demonstrate the importance of early embryonic, ontogenetic and tissue interactions in shaping craniofacial evolution and ecological diversification in squamates, with implications for the nature of cranio-cerebral relations across vertebrates.

Macroevolutionary drivers of morphological disparity in the avian quadrate

In birds, the quadrate connects the mandible and skull, and plays an important role in cranial kinesis. Avian quadrate morphology may therefore be assumed to have been influenced by selective pressures related to feeding ecology, yet large-scale variation in quadrate morphology and its potential relationship with ecology have never been quantitatively investigated. Here, we used geometric morphometrics and phylogenetic comparative methods to quantify morphological variation of the quadrate and its relationship with key ecological features across a wide phylogenetic sample. We found non-significant associations between quadrate shape and feeding ecology across different scales of phylogenetic comparison; indeed, allometry and phylogeny exhibit stronger relationships with quadrate shape than ecological features. We show that similar quadrate shapes are associated with widely varying dietary ecologies (one-to-many mapping), while divergent quadrate shapes are associated with similar dietary ecologies (many-to-one mapping). Moreover, we show that the avian quadrate evolves as an integrated unit and exhibits strong associations with the morphologies of neighbouring bones. Our results collectively illustrate that quadrate shape has evolved jointly with other elements of the avian kinetic system, with the major crown bird lineages exploring alternative quadrate morphologies, highlighting the potential diagnostic value of quadrate morphology in investigations of bird systematics.

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Drivers of morphological evolution in the toothed whale jaw

Toothed whales (odontocetes) emit high-frequency underwater sounds (echolocate)-an extreme and unique innovation allowing them to sense their prey and environment. Their highly specialized mandible (lower jaw) allows high-frequency sounds to be transmitted back to the inner ear. Echolocation is evident in the earliest toothed whales, but little research has focused on the evolution of mandibular form regarding this unique adaptation. Here, we use a high-density, three-dimensional geometric morphometric analysis of 100 living and extinct cetacean species spanning their ∼50-million-year evolutionary history. Our analyses demonstrate that most shape variation is found in the relative length of the jaw and the mandibular symphysis. The greatest morphological diversity was obtained during two periods of rapid evolution: the initial evolution of archaeocetes (stem whales) in the early to mid-Eocene as they adapted to an aquatic lifestyle, representing one of the most extreme adaptive transitions known, and later on in the mid-Oligocene odontocetes as they became increasingly specialized for a range of diets facilitated by increasingly refined echolocation. Low disparity in the posterior mandible suggests the shape of the acoustic window, which receives sound, has remained conservative since the advent of directional hearing in the aquatic archaeocetes, even as the earliest odontocetes began to receive sounds from echolocation. Diet, echolocation, feeding method, and dentition type strongly influence mandible shape. Unlike in the toothed whale cranium, we found no significant asymmetry in the mandible. We suggest that a combination of refined echolocation and associated dietary specializations have driven morphology and disparity in the toothed whale mandible.

A landmarking protocol for geometric morphometric analysis of squamate endocasts

Landmark-based geometric morphometrics is widely used to study the morphology of the endocast, or internal mold of the braincase, and the diversity associated with this structure across vertebrates. Landmarks, as the basic unit of such methods, are intended to be points of correspondence, selected depending on the question at hand, whose proper definition is essential to guarantee robustness and reproducibility of results. In this study, 20 landmarks are defined to provide a framework to analyze the morphological variability in squamate endocasts. Ten species representing a cross-section of the diversity of Squamata from both phylogenetic and ecological (i.e., habitat) perspectives were considered, to select landmarks replicable throughout the entire clade, regardless of the degree of neuroanatomical resolution of the endocast. To assess the precision, accuracy, and repeatability of these newly defined landmarks, both intraobserver and interobserver error were investigated. Estimates of measurement error show that most of the landmarks established here are highly replicable, and preliminary results suggest that they capture aspects of endocast shape related to both phylogenetic and ecologic signals. This study provides a basis for further examinations of squamate endocast disparity using landmark-based geometric morphometrics.

Endocast, brain, and bones: Correspondences and spatial relationships in squamates

Vertebrate endocasts are widely used in the fields of paleoneurology and comparative neuroanatomy. The validity of endocranial studies is dependent upon the extent to which an endocast reflects brain morphology. Due to the variable neuroanatomical resolution of vertebrate endocasts, direct information about the brain morphology can be sometimes difficult to assess and needs to be investigated across lineages. Here, we employ X-ray computed tomography (CT), including diffusible iodine-based contrast-enhanced CT, to qualitatively compare brains and endocasts in different species of squamates. The relative position of the squamate brain within the skull, as well as the variability that may exist in such spatial relationships, was examined to help clarify the neurological regions evidence on their endocasts. Our results indicate that squamate endocasts provide variable representation of the brain, depending on species and neuroanatomical regions. The olfactory bulbs and peduncles, cerebral hemispheres, as well as the medulla oblongata represent the most easily discernable brain regions from squamate endocasts. In contrast, the position of the optic lobes, the ventral diencephalon and the pituitary may be difficult to determine depending on species. Finally, squamate endocasts provide very limited or no information about the cerebellum. The spatial relationships revealed here between the brain and the surrounding bones may help to identify each of the endocranial region. However, as one-to-one correspondences between a bone and a specific region appear limited, the exact delimitation of these regions may remain challenging according to species. This study provides a basis for further examination and interpretation of squamate endocast disparity.

The affinities of the Late Triassic Cryptovaranoides and the age of crown squamates

Most living reptile diversity is concentrated in Squamata (lizards, including snakes), which have poorly known origins in space and time. Recently, †Cryptovaranoides microlanius from the Late Triassic of the United Kingdom was described as the oldest crown squamate. If true, this result would push back the origin of all major lizard clades by 30-65 Myr and suggest that divergence times for reptile clades estimated using genomic and morphological data are grossly inaccurate. Here, we use computed tomography scans and expanded phylogenetic datasets to re-evaluate the phylogenetic affinities of †Cryptovaranoides and other putative early squamates. We robustly reject the crown squamate affinities of †Cryptovaranoides, and instead resolve †Cryptovaranoides as a potential member of the bird and crocodylian total clade, Archosauromorpha. Bayesian total evidence dating supports a Jurassic origin of crown squamates, not Triassic as recently suggested. We highlight how features traditionally linked to lepidosaurs are in fact widespread across Triassic reptiles. Our study reaffirms the importance of critically choosing and constructing morphological datasets and appropriate taxon sampling to test the phylogenetic affinities of problematic fossils and calibrate the Tree of Life.

Reconstructing the origin and early evolution of the snake brain

Snakes represent one-eighth of terrestrial vertebrate diversity, encompassing various lifestyles, ecologies, and morphologies. However, the ecological origins and early evolution of snakes are controversial topics in biology. To address the paucity of well-preserved fossils and the caveats of osteological traits for reconstructing snake evolution, we applied a different ecomorphological hypothesis based on high-definition brain reconstructions of extant Squamata. Our predictive models revealed a burrowing lifestyle with opportunistic behavior at the origin of crown snakes, reflecting a complex ancestral mosaic brain pattern. These findings emphasize the importance of quantitatively tracking the phenotypic diversification of soft tissues-including the accurate definition of intact brain morphological traits such as the cerebellum-in understanding snake evolution and vertebrate paleobiology. Furthermore, our study highlights the power of combining extant and extinct species, soft tissue reconstructions, and osteological traits in tracing the deep evolution of not only snakes but also other groups where fossil data are scarce.

Detecting natural selection in trait-trait coevolution

No phenotypic trait evolves independently of all other traits, but the cause of trait-trait coevolution is poorly understood. While the coevolution could arise simply from pleiotropic mutations that simultaneously affect the traits concerned, it could also result from multivariate natural selection favoring certain trait relationships. To gain a general mechanistic understanding of trait-trait coevolution, we examine the evolution of 220 cell morphology traits across 16 natural strains of the yeast Saccharomyces cerevisiae and the evolution of 24 wing morphology traits across 110 fly species of the family Drosophilidae, along with the variations of these traits among gene deletion or mutation accumulation lines (a.k.a. mutants). For numerous trait pairs, the phenotypic correlation among evolutionary lineages differs significantly from that among mutants. Specifically, we find hundreds of cases where the evolutionary correlation between traits is strengthened or reversed relative to the mutational correlation, which, according to our population genetic simulation, is likely caused by multivariate selection. Furthermore, we detect selection for enhanced modularity of the yeast traits analyzed. Together, these results demonstrate that trait-trait coevolution is shaped by natural selection and suggest that the pleiotropic structure of mutation is not optimal. Because the morphological traits analyzed here are chosen largely because of their measurability and thereby are not expected to be biased with regard to natural selection, our conclusion is likely general.

Ecological correlates of cranial evolution in the megaradiation of dipsadine snakes

BACKGROUND

Dipsadine snakes represent one of the most spectacular vertebrate radiations that have occurred in any continental setting, with over 800 species in South and Central America. Their species richness is paralleled by stunning ecological diversity, ranging from arboreal snail-eating and aquatic eel-eating specialists to terrestrial generalists. Despite the ecological importance of this clade, little is known about the extent to which ecological specialization shapes broader patterns of phenotypic diversity within the group. Here, we test how habitat use and diet have influenced morphological diversification in skull shape across 160 dipsadine species using micro-CT and 3-D geometric morphometrics, and we use a phylogenetic comparative approach to test the contributions of habitat use and diet composition to variation in skull shape among species.

RESULTS

We demonstrate that while both habitat use and diet are significant predictors of shape in many regions of the skull, habitat use significantly predicts shape in a greater number of skull regions when compared to diet. We also find that across ecological groupings, fossorial and aquatic behaviors result in the strongest deviations in morphospace for several skull regions. We use simulations to address the robustness of our results and describe statistical anomalies that can arise from the application of phylogenetic generalized least squares to complex shape data.

CONCLUSIONS

Both habitat and dietary ecology are significantly correlated with skull shape in dipsadines; the strongest relationships involved skull shape in snakes with aquatic and fossorial lifestyles. This association between skull morphology and multiple ecological axes is consistent with a classic model of adaptive radiation and suggests that ecological factors were an important component in driving morphological diversification in the dipsadine megaradiation.

Bacterial community analysis identifies Klebsiella pneumoniae as a native symbiotic bacterium in the newborn Protobothrops mucrosquamatus

BACKGROUND

The study of the native microbiome of organisms is crucial. The connection between the native microbiome and the host affects the formation of the innate immune system and the organism's growth. However, the native microbiome of newborn venomous snakes has not been reported. Therefore, we aimed to determine the oral and skin microbiomes of newborn Protobothrops mucrosquamatus.

RESULTS

We performed 16 S full-length sequencing on 14 samples collected from 7 newborn P. mucrosquamatus individuals, specifically targeting their oral and skin microbiomes. In terms of the oral and skin microbiome, the main species were Klebsiella pneumoniae lineages. According to subspecies/species analysis, the proportion from highest to lowest was K. quasipneumoniae subsp. similipneumoniae, K. pneumoniae subsp. pneumoniae, and K. pneumoniae subsp. rhinoscleromatis. These three bacteria accounted for 62.5% and 85% of the skin and oral activity, respectively. The oral microbiome of newborn P. mucrosquamatus did not comprise common bacteria found in snakebite wounds or oral cultures in adult snakes. Therefore, the source of other microbiomes in the oral cavities of adult snakes may be the environment or prey. Functional Annotation of the Prokaryotic Taxa analysis showed that the skin/oral native microbiome metabolism was related to fermentation and human infection owing to the dominance of K. pneumoniae lineages. The characteristics of K. pneumoniae may impact the development of venom in venomous snakes.

CONCLUSION

The results of the native microbiome in the oral cavity and skin of newborn P. mucrosquamatus demonstrated that the habitat environment and prey capture may affect the composition of bacteria in adult snakes. We hypothesized that the native microbiome influences newborn venomous snakes and that K. pneumoniae lineages related to citrate fermentation may play a role in venom growth. However, further verification of this is required.

Disentangling Mechanical and Sensory Modules in the Radiation of Noctilionoid Bats

AbstractWith diverse mechanical and sensory functions, the vertebrate cranium is a complex anatomical structure whose shifts between modularity and integration, especially in mechanical function, have been implicated in adaptive diversification. Yet how mechanical and sensory systems and their functions coevolve, as well as how their interrelationship contributes to phenotypic disparity, remain largely unexplored. To examine the modularity, integration, and evolutionary rates of sensory and mechanical structures within the head, we analyzed hard and soft tissue scans from ecologically diverse bats in the superfamily Noctilionoidea, a clade that ranges from insectivores and carnivores to frugivores and nectarivores. We identified eight regions that evolved in a coordinated fashion, thus recognizable as evolutionary modules: five associated with bite force and three linked to olfactory, visual, and auditory systems. Interrelationships among these modules differ between Neotropical leaf-nosed bats (family Phyllostomidae) and other noctilionoids. Consistent with the hypothesis that dietary transitions begin with changes in the capacity to detect novel food items followed by adaptations to process them, peak rates of sensory module evolution predate those of some mechanical modules. We propose that the coevolution of structures influencing bite force, olfaction, vision, and hearing constituted a structural opportunity that allowed the phyllostomid ancestor to take advantage of existing ecological opportunities and contributed to the clade's remarkable radiation.

When Food Fights Back: Skull Morphology and Feeding Behavior of Centipede-Eating Snakes

Feeding is a complex process that involves an integrated response of multiple functional systems. Animals evolve phenotypic integration of complex morphological traits to covary and maximize performance of feeding behaviors. Specialization, such as feeding on dangerous prey, can further shape the integration of behavior and morphology as traits are expected to evolve and maintain function in parallel. Feeding on centipedes, with their powerful forcipules that pinch and inject venom, has evolved multiple times within snakes, including the genus Tantilla. However, the behavioral and morphological adaptations used to consume this dangerous prey are poorly understood. By studying snakes with varying degrees of dietary specialization, we can test the integration of diet, morphology, and behavior to better understand the evolution of consuming difficult prey. We studied the prey preference and feeding behavior of Tantilla using the flat-headed snake (T. gracilis) and the crowned snake (T. coronata), which differ in the percentage of centipedes in their diet. We then quantified cranial anatomy using geometric morphometric data from CT scans. To test prey preference, we offered multiple types of prey and recorded snake behavior. Both species of snakes showed interest in multiple prey types, but only struck or consumed centipedes. To subdue centipedes, crowned snakes used coiling and holding (envenomation) immediately after striking, while flat-headed snakes used the novel behavior of pausing and holding onto centipedes for a prolonged time prior to the completion of swallowing. Each skull element differed in shape after removing the effects of size, position, and orientation. The rear fang was larger in crowned snakes, but the mechanical advantage of the lower jaw was greater in flat-headed snakes. Our results suggest that the integration of behavioral and morphological adaptations is important for the success of subduing and consuming dangerous prey.

A switch in jaw form-function coupling during the evolution of mammals

The evolutionary shift from a single-element ear, multi-element jaw to a multi-element ear, single-element jaw during the transition to crown mammals marks one of the most dramatic structural transformations in vertebrates. Research on this transformation has focused on mammalian middle-ear evolution, but a mandible comprising only the dentary is equally emblematic of this evolutionary radiation. Here, we show that the remarkably diverse jaw shapes of crown mammals are coupled with surprisingly stereotyped jaw stiffness. This strength-based morphofunctional regime has a genetic basis and allowed mammalian jaws to effectively resist deformation as they radiated into highly disparate forms with markedly distinct diets. The main functional consequences for the mandible of decoupling hearing and mastication were a trade-off between higher jaw stiffness versus decreased mechanical efficiency and speed compared with non-mammals. This fundamental and consequential shift in jaw form-function underpins the ecological and taxonomic diversification of crown mammals. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Developmental origin underlies evolutionary rate variation across the placental skull

The placental skull has evolved into myriad forms, from longirostrine whales to globular primates, and with a diverse array of appendages from antlers to tusks. This disparity has recently been studied from the perspective of the whole skull, but the skull is composed of numerous elements that have distinct developmental origins and varied functions. Here, we assess the evolution of the skull's major skeletal elements, decomposed into 17 individual regions. Using a high-dimensional morphometric approach for a dataset of 322 living and extinct eutherians (placental mammals and their stem relatives), we quantify patterns of variation and estimate phylogenetic, allometric and ecological signal across the skull. We further compare rates of evolution across ecological categories and ordinal-level clades and reconstruct rates of evolution along lineages and through time to assess whether developmental origin or function discriminate the evolutionary trajectories of individual cranial elements. Our results demonstrate distinct macroevolutionary patterns across cranial elements that reflect the ecological adaptations of major clades. Elements derived from neural crest show the fastest rates of evolution, but ecological signal is equally pronounced in bones derived from neural crest and paraxial mesoderm, suggesting that developmental origin may influence evolutionary tempo, but not capacity for specialisation. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

High-Density Geometric Morphometric Analysis of Intraspecific Cranial Integration in the Barred Grass Snake (Natrix helvetica) and Green Anole (Anolis carolinensis)

How do phenotypic associations intrinsic to an organism, such as developmental and mechanical processes, direct morphological evolution? Comparisons of intraspecific and clade-wide patterns of phenotypic covariation could inform how population-level trends ultimately dictate macroevolutionary changes. However, most studies have focused on analyzing integration and modularity either at macroevolutionary or intraspecific levels, without a shared analytical framework unifying these temporal scales. In this study, we investigate the intraspecific patterns of cranial integration in two squamate species: Natrix helvetica and Anolis carolinensis. We analyze their cranial integration patterns using the same high-density three-dimensional geometric morphometric approach used in a prior squamate-wide evolutionary study. Our results indicate that Natrix and Anolis exhibit shared intraspecific cranial integration patterns, with some differences, including a more integrated rostrum in the latter. Notably, these differences in intraspecific patterns correspond to their respective interspecific patterns in snakes and lizards, with few exceptions. These results suggest that interspecific patterns of cranial integration reflect intraspecific patterns. Hence, our study suggests that the phenotypic associations that direct morphological variation within species extend across micro- and macroevolutionary levels, bridging these two scales.

Effects of COVID-19 Infection and Vaccination on the Female Reproductive System: A Narrative Review

Several studies and research papers have been published to elucidate and understand the mechanism of the coronavirus disease 2019 (COVID-19) pandemic and its long-term effects on the human body. COVID-19 affects a number of organs, including the female reproductive system. However, less attention has been given to the effects of COVID-19 on the female reproductive system due to their low morbidity. The results of studies investigating the relationship between COVID-19 infection and ovarian function in women of reproductive age have shown the harmless involvement of COVID-19 infection. Several studies have reported the involvement of COVID-19 infection in oocyte quality, ovarian function, and dysfunctions in the uterine endometrium and the menstrual cycle. The findings of these studies indicate that COVID-19 infection negatively affects the follicular microenvironment and dysregulate ovarian function. Although the COVID-19 pandemic and female reproductive health have been studied in humans and animals, very few studies have examined how COVID-19 affects the female reproductive system. The objective of this review is to summarize the current literature and categorize the effects of COVID-19 on the female reproductive system, including the ovaries, uterus, and hormonal profiles. The effects on oocyte maturation, oxidative stress, which causes chromosomal instability and apoptosis in ovaries, in vitro fertilization cycle, high-quality embryos, premature ovarian insufficiency, ovarian vein thrombosis, hypercoagulable state, women’s menstrual cycle, the hypothalamus-pituitary-ovary axis, and sex hormones, including estrogen, progesterone, and the anti-Müllerian hormone, are discussed in particular.

Using pose estimation to identify regions and points on natural history specimens

A key challenge in mobilising growing numbers of digitised biological specimens for scientific research is finding high-throughput methods to extract phenotypic measurements on these datasets. In this paper, we test a pose estimation approach based on Deep Learning capable of accurately placing point labels to identify key locations on specimen images. We then apply the approach to two distinct challenges that each requires identification of key features in a 2D image: (i) identifying body region-specific plumage colouration on avian specimens and (ii) measuring morphometric shape variation in Littorina snail shells. For the avian dataset, 95% of images are correctly labelled and colour measurements derived from these predicted points are highly correlated with human-based measurements. For the Littorina dataset, more than 95% of landmarks were accurately placed relative to expert-labelled landmarks and predicted landmarks reliably captured shape variation between two distinct shell ecotypes ('crab' vs 'wave'). Overall, our study shows that pose estimation based on Deep Learning can generate high-quality and high-throughput point-based measurements for digitised image-based biodiversity datasets and could mark a step change in the mobilisation of such data. We also provide general guidelines for using pose estimation methods on large-scale biological datasets.

Flexible conservatism in the skull modularity of convergently evolved myrmecophagous placental mammals

Background

The skull of placental mammals constitutes one of the best studied systems for phenotypic modularity. Several studies have found strong evidence for the conserved presence of two- and six-module architectures, while the strength of trait correlations (integration) has been associated with major developmental processes such as somatic growth, muscle-bone interactions, and tooth eruption. Among placentals, ant- and termite-eating (myrmecophagy) represents an exemplar case of dietary convergence, accompanied by the selection of several cranial morphofunctional traits such as rostrum elongation, tooth loss, and mastication loss. Despite such drastic functional modifications, the covariance patterns of the skull of convergently evolved myrmecophagous placentals are yet to be studied in order to assess the potential consequences of this dietary shift on cranial modularity.

Results

Here, we performed a landmark-based morphometric analysis of cranial covariance patterns in 13 species of myrmecophagous placentals. Our analyses reveal that most myrmecophagous species present skulls divided into six to seven modules (depending on the confirmatory method used), with architectures similar to those of non-myrmecophagous placentals (therian six modules). Within-module integration is also similar to what was previously described for other placentals, suggesting that most covariance-generating processes are conserved across the clade. Nevertheless, we show that extreme rostrum elongation and tooth loss in myrmecophagid anteaters have resulted in a shift in intermodule correlations in the proximal region of the rostrum. Namely, the naso-frontal and maxillo-palatine regions are strongly correlated with the oro-nasal module, suggesting an integrated rostrum conserved from pre-natal developmental processes. In contrast, the similarly toothless pangolins show a weaker correlation between the anterior rostral modules, resembling the pattern of toothed placentals.

Conclusions

These results reveal that despite some integration shifts related to extreme functional and morphological features of myrmecophagous skulls, cranial modular architectures have conserved the typical mammalian scheme.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02030-9.

Evolutionary origins of the prolonged extant squamate radiation

Squamata is the most diverse clade of terrestrial vertebrates. Although the origin of pan-squamates lies in the Triassic, the oldest undisputed members of extant clades known from nearly complete, uncrushed material come from the Cretaceous. Here, we describe three-dimensionally preserved partial skulls of two new crown lizards from the Late Jurassic of North America. Both species are placed at the base of the skink, girdled, and night lizard clade Pan-Scincoidea, which consistently occupies a position deep inside the squamate crown in both morphological and molecular phylogenies. The new lizards show that several features uniting pan-scincoids with another major lizard clade, the pan-lacertoids, in trees using morphology were convergently acquired as predicted by molecular analyses. Further, the palate of one new lizard bears a handful of ancestral saurian characteristics lost in nearly all extant squamates, revealing an underappreciated degree of complex morphological evolution in the early squamate crown. We find strong evidence for close relationships between the two new species and Cretaceous taxa from Eurasia. Together, these results suggest that early crown squamates had a wide geographic distribution and experienced complicated morphological evolution even while the Rhynchocephalia, now solely represented by the tuatara, was the dominant clade of lepidosaurs.

Evidence for an evo-devo-derived hypothesis on three-dimensional flower shape modularity in a tropical orchid clade

Covarying suites of phenotypic traits, or modules, are increasingly recognized to promote morphological evolution. However, information on how modularity influences flower diversity is rare and lacking for Orchidaceae. Here, we combine high-resolution X-ray computed tomography scanning with three-dimensional geometric morphometrics and phylogenetic comparative methods to test various hypotheses about three-dimensional patterns of flower evolutionary modularity in Malagasy Bulbophyllum orchids and examine rates and modes of module evolution. Based on the four evolutionary modules identified (i.e., sepals, lateral petals, labellum + column-foot, and column-part), our data support the hypothesis that both genetic-developmental and functional adaptive factors shaped evolutionary flower trait covariation in these tropical orchids. In line with "evo-devo" studies, we also find that the labellum evolved independently from the rest of the petal whorl. Finally, we show that modules evolved with different rates, and either in a neutral fashion (only column-part) or under selective constraints, as likely imposed by pollinators. Overall, this study supports current views that modular units can enhance the range and rate of morphological evolution.

Anatomical network analyses reveal evolutionary integration and modularity in the lizards skull

The morphology of lizard skulls is highly diverse, and it is crucial to understand the factors that constrain and promote their evolution to understand how lizards thrive. The results of interactions between cranial bones reflecting these factors can be detected as integration and modularity, and the analysis of integration and modularity allows us to explore the underlying factors. In this study, the integration and modularity of the skulls of lizards and the outgroup tuatara are analyzed using a new method, Anatomical Network Analysis (AnNA), and the factors causing lizards morphological diversity are investigated by comparing them. The comparison of modular structures shows that lizard skulls have high integration and anisomerism, some differences but basically common modular patterns. In contrast, the tuatara shows a different modular pattern from lizards. In addition, the presence of the postorbital bar by jugal and postorbital (postorbitofrontal) also reflect various functional factors by maintaining low integration. The maintenance of basic structures due to basic functional requirements and changes in integration within the modules play a significant role in increasing the morphological diversity of the lizard skull and in the prosperity of the lizards.

Convergence, divergence, and macroevolutionary constraint as revealed by anatomical network analysis of the squamate skull, with an emphasis on snakes

Traditionally considered the earliest-diverging group of snakes, scolecophidians are central to major evolutionary paradigms regarding squamate feeding mechanisms and the ecological origins of snakes. However, quantitative analyses of these phenomena remain scarce. Herein, we therefore assess skull modularity in squamates via anatomical network analysis, focusing on the interplay between ‘microstomy’ (small-gaped feeding), fossoriality, and miniaturization in scolecophidians. Our analyses reveal distinctive patterns of jaw connectivity across purported ‘microstomatans’, thus supporting a more complex scenario of jaw evolution than traditionally portrayed. We also find that fossoriality and miniaturization each define a similar region of topospace (i.e., connectivity-based morphospace), with their combined influence imposing further evolutionary constraint on skull architecture. These results ultimately indicate convergence among scolecophidians, refuting widespread perspectives of these snakes as fundamentally plesiomorphic and morphologically homogeneous. This network-based examination of skull modularity—the first of its kind for snakes, and one of the first to analyze squamates—thus provides key insights into macroevolutionary trends among squamates, with particular implications for snake origins and evolution.

Untangling the relationship between developmental and evolutionary integration

Patterns of integration and modularity among organismal traits are prevalent across the tree of life, and at multiple scales of biological organization. Over the past several decades, researchers have studied these patterns at the developmental, and evolutionary levels. While their work has identified the potential drivers of these patterns at different scales, there appears to be a lack of consensus on the relationship between developmental and evolutionary integration. Here, we review and summarize key studies and build a framework to describe the conceptual relationship between these patterns across organismal scales and illustrate how, and why some of these studies may have yielded seemingly conflicting outcomes. We find that among studies that analyze patterns of integration and modularity using morphological data, the lack of consensus may stem in part from the difficulty of fully disentangling the developmental and functional causes of integration. Nonetheless, in some empirical systems, patterns of evolutionary modularity have been found to coincide with expectations based on developmental processes, suggesting that in some circumstances, developmental modularity may translate to evolutionary modularity. We also advance an extension to Hallgrímsson et al.'s palimpsest model to describe how patterns of trait modularity may shift across different evolutionary scales. Finally, we also propose some directions for future research which will hopefully be useful for investigators interested in these issues.

Fast Likelihood Calculations for Automatic Identification of Macroevolutionary Rate Heterogeneity in Continuous and Discrete Traits

Understanding phenotypic disparity across the tree of life requires identifying where and when evolutionary rates change on phylogeny. A primary methodological challenge in macroevolution is therefore to develop methods for accurate inference of among-lineage variation in rates of phenotypic evolution. Here, we describe a method for inferring among-lineage evolutionary rate heterogeneity in both continuous and discrete traits. The method assumes that the present-day distribution of a trait is shaped by a variable-rate process arising from a mixture of constant-rate processes and uses a single-pass tree traversal algorithm to estimate branch-specific evolutionary rates. By employing dynamic programming optimization techniques and approximate maximum likelihood estimators where appropriate, our method permits rapid exploration of the tempo and mode of phenotypic evolution. Simulations indicate that the method reconstructs rates of trait evolution with high accuracy. Application of the method to datasets on squamate reptile reproduction and turtle body size recovers patterns of rate heterogeneity identified by previous studies but with computational costs reduced by many orders of magnitude. Our results expand the set of tools available for detecting macroevolutionary rate heterogeneity and point to the utility of fast, approximate methods for studying large scale biodiversity dynamics.

The Jurassic rise of squamates as supported by lepidosaur disparity and evolutionary rates

The squamates (lizards, snakes, and relatives) today comprise more than 10,000 species, and yet their sister group, the Rhynchocephalia, is represented by a single species today, the tuatara. The explosion in squamate diversity has been tracked back to the Cretaceous Terrestrial Revolution, 100 million years ago (Ma), the time when flowering plants began their takeover of terrestrial ecosystems, associated with diversification of coevolving insects and insect-eating predators such as lizards, birds, and mammals. Squamates arose much earlier, but their long pre-Cretaceous history of some 150 million years (Myr) is documented by sparse fossils. Here, we provide evidence for an initial radiation of squamate morphology in the Middle and Late Jurassic (174-145 Ma), and show that they established their key ecological roles much earlier than had been assumed, and they have not changed them much since.

Drosophila Wing Integration and Modularity: A Multi-Level Approach to Understand the History of Morphological Structures

Simple Summary

The diverse components of any morphological structure are integrated with respect to each other since they have developed, functioned, and evolved together, a phenomenon known as integration. However, this integration is not absolute but organized in units (i.e., modules) that are relatively independent while participating to generate a structure that acts as a functional whole. Even though most of the studies on modularity and integration have focused on variation among individuals within populations, there are more levels of variation that exhibit modularity and integration, deriving from distinct sources such as genetic variation, phenotypic plasticity, fluctuating asymmetry, evolutionary change, among others. Consequently, the present study focused on analysing the integration and modularity of the wing shape of some of the best-known model organisms, i.e., the genus Drosophila, at the static, developmental, and evolutionary levels to acquire a better insight about how modularity and integration act at different analytical levels. The strong integration and overall similarities observed in the variation pattern at multiple levels suggest a shared mechanism underlying the observed variation in Drosophila’s wing shape and added a new piece of evidence of stasis in the evolutionary history of Drosophila wing.

Abstract

Static, developmental, and evolutionary variation are different sources of morphological variation which can be quantified using morphometrics tools. In the present study we have carried out a comparative multiple level study of integration (i.e., static, developmental, and evolutionary) to acquire insight about the relationships that exist between different integration levels, as well as to better understand their involvement in the evolutionary processes related to the diversification of Drosophila’s wing shape. This approach was applied to analyse wing evolution in 59 species across the whole genus in a large dataset (~10,000 wings were studied). Static integration was analysed using principal component analysis, thus providing an integration measurement for overall wing shape. Developmental integration was studied between wing parts by using a partial least squares method between the anterior and posterior compartments of the wing. Evolutionary integration was analysed using independent contrasts. The present results show that all Drosophila species exhibit strong morphological integration at different levels. The strong integration and overall similarities observed at multiple integration levels suggest a shared mechanism underlying this variation, which could result as consequence of genetic drift acting on the wing shape of Drosophila.

Ontogenetic shift in diet of a large elapid snake is facilitated by allometric change in skull morphology

As snakes are limbless, gape-limited predators, their skull is the main feeding structure involved in prey handling, manipulation and feeding. Ontogenetic changes in prey type and size are likely to be associated with distinct morphological changes in the skull during growth. We investigated ontogenetic variation in diet from stomach contents of 161 Dugite specimens (Pseudonaja affinis, Elapidae) representing the full range of body size for the species, and skull morphology of 46 specimens (range 0.25–1.64 m snout-vent-length; SVL). We hypothesised that changes in prey type throughout postnatal ontogeny would coincide with distinct changes in skull shape. Dugites demonstrate a distinct size-related shift in diet: the smallest individuals ate autotomised reptile tails and reptiles, medium-sized individuals predominantly ate mammals, and the largest individuals had the most diverse diet, including large reptiles. Morphometric analysis revealed that ~40% of the variation in skull shape was associated with body size (SVL). Through ontogeny, skulls changed from a smooth, bulbous cranium with relatively small trophic bones (upper and lower jaws and their attachments), to more rugose bones (as a likely reflection of muscle attachment) and relatively longer trophic bones that would extend gape. Individual shape variation in trophic bone dimensions was greater in larger adults and this likely reflects natural plasticity of individuals feeding on different prey sizes/types. Rather than a distinct morphological shift with diet, the ontogenetic changes were gradual, but positive allometry of individual trophic bones resulted in disproportionate growth of the skull, reflected in increased gape size and mobility of jaw bones in adults to aid the ingestion of larger prey and improve manipulation and processing ability. These results indicate that allometric scaling is an important mechanism by which snakes can change their dietary niche.

Early origins of divergent patterns of morphological evolution on the mammal and reptile stem-lineages

The origin of amniotes 320 million years ago signaled independence from water in vertebrates and was closely followed by divergences within the mammal and reptile stem lineages (Synapsida and Reptilia). Early members of both groups had highly similar morphologies, being superficially “lizard-like” forms with many plesiomorphies. However, the extent to which they might have exhibited divergent patterns of evolutionary change, with the potential to explain the large biological differences between their living members, is unresolved. We use a new, comprehensive phylogenetic dataset to quantify variation in rates and constraints of morphological evolution among Carboniferous–early Permian amniotes. We find evidence for an early burst of evolutionary rates, resulting in the early origins of morphologically distinctive subgroups that mostly persisted through the Cisuralian. Rates declined substantially through time, especially in reptiles. Early reptile evolution was also more constrained compared with early synapsids, exploring a more limited character state space. Postcranial innovation in particular was important in early synapsids, potentially related to their early origins of large body size. In contrast, early reptiles predominantly varied the temporal region, suggesting disparity in skull and jaw kinematics, and foreshadowing the variability of cranial biomechanics seen in reptiles today. Our results demonstrate that synapsids and reptiles underwent an early divergence of macroevolutionary patterns. This laid the foundation for subsequent evolutionary events and may be critical in understanding the substantial differences between mammals and reptiles today. Potential explanations include an early divergence of developmental processes or of ecological factors, warranting cross-disciplinary investigation. [Amniote; body size; constraint; phylogeny; rate.]

X-ray microtomography imaging of craniofacial hard tissues in selected reptile species with different types of dentition

Background

Reptiles exhibit a large heterogeneity in teeth morphology. The main variability comprises the different tooth shape, the type of tooth attachment to the underlying bone, or the ability to replace the teeth.

Findings

Here, we provide full datasets of microtomography scans and 3D models of reptilian dentitions and skulls. We selected representative species for each of 9 reptilian families on the basis of their characteristic dental features. Because there are ≥4 different types of tooth-bone attachments, ranging from the mammalian-like thecodont attachment found in crocodilians to the simple acrodont implantation observed in some lizards, we aimed to evaluate species with different types of tooth-bone attachments. Moreover, another interesting feature varying in reptilian species is the complexity of tooth shape or the number of tooth generations, which can be associated with the type of tooth attachment to the jawbone. Therefore, selected model species also include animals with distinct tooth morphology along the jaw or different number of tooth generations. The development of tooth attachment and relationship of the tooth to the jaw can be further analysed in detail on a large collection of pre-hatching stages of chameleon. Next, we introduce different possibilities for how these datasets can be further used to study tooth-bone relationships or tooth morphology in 3D space. Moreover, these datasets can be valuable for additional morphological and morphometric analyses of reptilian skulls or their individually segmented skeletal elements.

Conclusions

Our collection of microcomputed tomography scans can bring new insight into dental or skeletal research. The broad selection of reptilian species, together with their unique dental features and high quality of these scans including complete series of developmental stages of our model species and provide large opportunities for their reuse. Scans can be further used for virtual reality, 3D printing, or in education.

Evolutionary modularity, integration and disparity in an accretionary skeleton: analysis of venerid Bivalvia

Modular evolution, the relatively independent evolution of body parts, may promote high morphological disparity in a clade. Conversely, integrated evolution via stronger covariation of parts may limit disparity. However, integration can also promote high disparity by channelling morphological evolution along lines of least resistance-a process that may be particularly important in the accumulation of disparity in the many invertebrate systems having accretionary growth. We use a time-calibrated phylogenetic hypothesis and high-density, three-dimensional semilandmarking to analyse the relationship between modularity, integration and disparity in the most diverse extant bivalve family: the Veneridae. In general, venerids have a simple, two-module parcellation of their body that is divided into features of the calcium carbonate shell and features of the internal soft anatomy. This division falls more along developmental than functional lines when placed in the context of bivalve anatomy and biomechanics. The venerid body is tightly integrated in absolute terms, but disparity appears to increase with modularity strength among subclades and ecologies. Thus, shifts towards more mosaic evolution beget higher morphological variance in this speciose family.

Covariation between the cranium and the cervical vertebrae in hominids

The analysis of patterns of integration is crucial for the reconstruction and understanding of how morphological changes occur in a taxonomic group throughout evolution. These patterns are relatively constant; however, both patterns and the magnitudes of integration may vary across species. These differences may indicate morphological diversification, in some cases related to functional adaptations to the biomechanics of organisms. In this study, we analyze patterns of integration between two functional and developmental structures, the cranium and the cervical spine in hominids, and we quantify the amount of divergence of each anatomical element through phylogeny. We applied these methods to three-dimensional data from 168 adult hominid individuals, summing a total of more than 1000 cervical vertebrae. We found the atlas (C1) and axis (C2) display the lowest covariation with the cranium in hominids (Homo sapiens, Pan troglodytes, Pan paniscus, Gorilla gorilla, Gorilla beringei, Pongo pygmaeus). H. sapiens show a relatively different pattern of craniocervical correlation compared with chimpanzees and gorillas, especially in variables implicated in maintaining the balance of the head. Finally, the atlas and axis show lower magnitude of shape change during evolution than the rest of the cervical vertebrae, especially those located in the middle of the subaxial cervical spine. Overall, results suggest that differences in the pattern of craniocervical correlation between humans and gorillas and chimpanzees could reflect the postural differences between these groups. Also, the stronger craniocervical integration and larger magnitude of shape change during evolution shown by the middle cervical vertebrae suggests that they have been selected to play an active role in maintaining head balance.

How to Render Species Comparable Taxonomic Units Through Deep Time: a Case Study on Intraspecific Osteological Variability in Extant and Extinct Lacertid Lizards

Generally, the species is considered to be the only naturally occurring taxon. However, species recognised and defined using different species delimitation criteria cannot readily be compared, impacting studies of biodiversity through Deep Time. This comparability issue is particularly marked when comparing extant with extinct species, because the only available data for species delimitation in fossils is derived from their preserved morphology, which is generally restricted to osteology in vertebrates. Here, we quantify intraspecific, intrageneric, and intergeneric osteological variability in extant species of lacertid lizards using pairwise dissimilarity scores based on a dataset of 253 discrete osteological characters for 99 specimens referred to 24 species. Variability is always significantly lower intraspecifically than between individuals belonging to distinct species of a single genus, which is in turn significantly lower than intergeneric variability. Average values of intraspecific variability and associated standard deviations are consistent (with few exceptions), with an overall average within a species of 0.208 changes per character scored. Application of the same methods to six extinct lacertid species (represented by 40 fossil specimens) revealed that intraspecific osteological variability is inconsistent, which can at least in part be attributed to different researchers having unequal expectations of the skeletal dissimilarity within species units. Such a divergent interpretation of intraspecific and interspecific variability among extant and extinct species reinforces the incomparability of the species unit. Lacertidae is an example where extant species recognised and defined based on a number of delimitation criteria show comparable and consistent intraspecific osteological variability. Here, as well as in equivalent cases, application of those skeletal dissimilarity values to palaeontological species delimitation potentially provides a way to ameliorate inconsistencies created by the use of morphology to define species.

The cichlid oral and pharyngeal jaws are evolutionarily and genetically coupled

Evolutionary constraints may significantly bias phenotypic change, while "breaking" from such constraints can lead to expanded ecological opportunity. Ray-finned fishes have broken functional constraints by developing two jaws (oral-pharyngeal), decoupling prey capture (oral jaw) from processing (pharyngeal jaw). It is hypothesized that the oral and pharyngeal jaws represent independent evolutionary modules and this facilitated diversification in feeding architectures. Here we test this hypothesis in African cichlids. Contrary to our expectation, we find integration between jaws at multiple evolutionary levels. Next, we document integration at the genetic level, and identify a candidate gene, smad7, within a pleiotropic locus for oral and pharyngeal jaw shape that exhibits correlated expression between the two tissues. Collectively, our data show that African cichlid evolutionary success has occurred within the context of a coupled jaw system, an attribute that may be driving adaptive evolution in this iconic group by facilitating rapid shifts between foraging habitats, providing an advantage in a stochastic environment such as the East African Rift-Valley.

What does modularity mean?

Modularity is now generally recognized as a fundamental feature of organisms, one that may have profound consequences for evolution. Modularity has recently become a major focus of research in organismal biology across multiple disciplines including genetics, developmental biology, functional morphology, population and evolutionary biology. While the wealth of new data, and also new theory, has provided exciting and novel insights, the concept of modularity has become increasingly ambiguous. That ambiguity is underlain by diverse intuitions about what modularity means, and the ambiguity is not merely about the meaning of the word-the metrics of modularity are measuring different properties and the methods for delimiting modules delimit them by different, sometimes conflicting criteria. The many definitions, metrics and methods can lead to substantial confusion not just about what modularity means as a word but also about what it means for evolution. Here we review various concepts, using graphical depictions of modules. We then review some of the metrics and methods for analyzing modularity at different levels. To place these in theoretical context, we briefly review theories about the origins and evolutionary consequences of modularity. Finally, we show how mismatches between concepts, metrics and methods can produce theoretical confusion, and how potentially illogical interpretations can be made sensible by a better match between definitions, metrics, and methods.

The limits of convergence: the roles of phylogeny and dietary ecology in shaping non-avian amniote crania

Cranial morphology is remarkably varied in living amniotes and the diversity of shapes is thought to correspond with feeding ecology, a relationship repeatedly demonstrated at smaller phylogenetic scales, but one that remains untested across amniote phylogeny. Using a combination of morphometric methods, we investigate the links between phylogenetic relationships, diet and skull shape in an expansive dataset of extant toothed amniotes: mammals, lepidosaurs and crocodylians. We find that both phylogeny and dietary ecology have statistically significant effects on cranial shape. The three major clades largely partition morphospace with limited overlap. Dietary generalists often occupy clade-specific central regions of morphospace. Some parallel changes in cranial shape occur in clades with distinct evolutionary histories but similar diets. However, members of a given clade often present distinct cranial shape solutions for a given diet, and the vast majority of species retain the unique aspects of their ancestral skull plan, underscoring the limits of morphological convergence due to ecology in amniotes. These data demonstrate that certain cranial shapes may provide functional advantages suited to particular dietary ecologies, but accounting for both phylogenetic history and ecology can provide a more nuanced approach to inferring the ecology and functional morphology of cryptic or extinct amniotes.

Complex macroevolutionary dynamics underly the evolution of the crocodyliform skull

All modern crocodyliforms (alligators, crocodiles and the gharial) are semi-aquatic generalist carnivores that are relatively similar in cranial form and function. However, this homogeneity represents just a fraction of the variation that once existed in the clade, which includes extinct herbivorous and marine forms with divergent skull structure and function. Here, we use high-dimensional three-dimensional geometric morphometrics to quantify whole-skull morphology across modern and fossil crocodyliforms to untangle the factors that shaped the macroevolutionary history and relatively low phenotypic variation of this clade through time. Evolutionary modelling demonstrates that the pace of crocodyliform cranial evolution is initially high, particularly in the extinct Notosuchia, but slows near the base of Neosuchia, with a late burst of rapid evolution in crown-group crocodiles. Surprisingly, modern crocodiles, especially Australian, southeast Asian, Indo-Pacific species, have high rates of evolution, despite exhibiting low variation. Thus, extant lineages are not in evolutionary stasis but rather have rapidly fluctuated within a limited region of morphospace, resulting in significant convergence. The structures related to jaw closing and bite force production (e.g. pterygoid flange and quadrate) are highly variable, reinforcing the importance of function in driving phenotypic variation. Together, these findings illustrate that the apparent conservativeness of crocodyliform skulls betrays unappreciated complexity in their macroevolutionary dynamics.