Uncovering the mosaic evolution of the carnivoran skeletal system

The diversity of vertebrate skeletons is often attributed to adaptations to distinct ecological factors such as diet, locomotion, and sensory environment. Although the adaptive evolution of skull, appendicular skeleton, and vertebral column is well studied in vertebrates, comprehensive investigations of all skeletal components simultaneously are rarely performed. Consequently, we know little of how modes of evolution differ among skeletal components. Here, we tested if ecological and phylogenetic effects led to distinct modes of evolution among the cranial, appendicular and vertebral regions in extant carnivoran skeletons. Using multivariate evolutionary models, we found mosaic evolution in which only the mandible, hindlimb and posterior (i.e. last thoracic and lumbar) vertebrae showed evidence of adaptation towards ecological regimes whereas the remaining skeletal components reflect clade-specific evolutionary shifts. We hypothesize that the decoupled evolution of individual skeletal components may have led to the origination of distinct adaptive zones and morphologies among extant carnivoran families that reflect phylogenetic hierarchies. Overall, our work highlights the importance of examining multiple skeletal components simultaneously in ecomorphological analyses. Ongoing work integrating the fossil and palaeoenvironmental record will further clarify deep-time drivers that govern the carnivoran diversity we see today and reveal the complexity of evolutionary processes in multicomponent systems.

Patterns of ontogenetic evolution across extant marsupials reflect different allometric pathways to ecomorphological diversity

The relatively high level of morphological diversity in Australasian marsupials compared to that observed among American marsupials remains poorly understood. We undertake a comprehensive macroevolutionary analysis of ontogenetic allometry of American and Australasian marsupials to examine whether the contrasting levels of morphological diversity in these groups are reflected in their patterns of allometric evolution. We collate ontogenetic series for 62 species and 18 families of marsupials (n = 2091 specimens), spanning across extant marsupial diversity. Our results demonstrate significant lability of ontogenetic allometric trajectories among American and Australasian marsupials, yet a phylogenetically structured pattern of allometric evolution is preserved. Here we show that species diverging more than 65 million years ago converge in their patterns of ontogenetic allometry under animalivorous and herbivorous diets, and that Australasian marsupials do not show significantly greater variation in patterns of ontogenetic allometry than their American counterparts, despite displaying greater magnitudes of extant ecomorphological diversity.

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.

Topographically distinct adaptive landscapes for teeth, skeletons, and size explain the adaptive radiation of Carnivora (Mammalia)

Models of adaptive radiation were originally developed to explain the early, rapid appearance of distinct modes of life within diversifying clades. Phylogenetic tests of this hypothesis have yielded limited support for temporally declining rates of phenotypic evolution across diverse clades, but the concept of an adaptive landscape that links form to fitness, while also crucial to these models, has received more limited attention. Using methods that assess the temporal accumulation of morphological variation and estimate the topography of the underlying adaptive landscape, I found evidence of an early partitioning of mandibulo-dental morphological variation in Carnivora (Mammalia) that occurs on an adaptive landscape with multiple peaks, consistent with classic ideas about adaptive radiation. Although strong support for this mode of adaptive radiation is present in traits related to diet, its signal is not present in body mass data or for traits related to locomotor behavior and substrate use. These findings suggest that adaptive radiations may occur along some axes of ecomorphological variation without leaving a signal in others and that their dynamics are more complex than simple univariate tests might suggest.

Marsupials and Multi-Omics: Establishing New Comparative Models of Neural Crest Patterning and Craniofacial Development

Studies across vertebrates have revealed significant insights into the processes that drive craniofacial morphogenesis, yet we still know little about how distinct facial morphologies are patterned during development. Studies largely point to evolution in GRNs of cranial progenitor cell types such as neural crest cells, as the major driver underlying adaptive cranial shapes. However, this hypothesis requires further validation, particularly within suitable models amenable to manipulation. By utilizing comparative models between related species, we can begin to disentangle complex developmental systems and identify the origin of species-specific patterning. Mammals present excellent evolutionary examples to scrutinize how these differences arise, as sister clades of eutherians and marsupials possess suitable divergence times, conserved cranial anatomies, modular evolutionary patterns, and distinct developmental heterochrony in their NCC behaviours and craniofacial patterning. In this review, I lend perspectives into the current state of mammalian craniofacial biology and discuss the importance of establishing a new marsupial model, the fat-tailed dunnart, for comparative research. Through detailed comparisons with the mouse, we can begin to decipher mammalian conserved, and species-specific processes and their contribution to craniofacial patterning and shape disparity. Recent advances in single-cell multi-omics allow high-resolution investigations into the cellular and molecular basis of key developmental processes. As such, I discuss how comparative evolutionary application of these tools can provide detailed insights into complex cellular behaviours and expression dynamics underlying adaptive craniofacial evolution. Though in its infancy, the field of "comparative evo-devo-omics" presents unparalleled opportunities to precisely uncover how phenotypic differences arise during development.

Sensory adaptations reshaped intrinsic factors underlying morphological diversification in bats

Background

Morphological evolution may be impacted by both intrinsic (developmental, constructional, physiological) and extrinsic (ecological opportunity and release) factors, but can intrinsic factors be altered by adaptive evolution and, if so, do they constrain or facilitate the subsequent diversification of biological form? Bats underwent deep adaptive divergences in skull shape as they evolved different sensory modes; here we investigate the potential impact of this process on two intrinsic factors that underlie morphological variation across organisms, allometry, and modularity.

Results

We use comparative phylogenetic and morphometric approaches to examine patterns of evolutionary allometry and modularity across a 3D geometric morphometric dataset spanning all major bat clades. We show that allometric relationships diverge between echolocators and visually oriented non-echolocators and that the evolution of nasal echolocation reshaped the modularity of the bat cranium.

Conclusions

Shifts in allometry and modularity may have significant consequences on the diversification of anatomical structures, as observed in the bat skull.

Anatomical comparison across heads, fore- and hindlimbs in mammals using network models

Animal body parts evolve with variable degrees of integration that nonetheless yield functional adult phenotypes: but, how? The analysis of modularity with Anatomical Network Analysis (AnNA) is used to quantitatively determine phenotypic modules based on the physical connection among anatomical elements, an approach that is valuable to understand developmental and evolutionary constraints. We created anatomical network models of the head, forelimb, and hindlimb of two taxa considered to represent a 'generalized' eutherian (placental: mouse) and metatherian (marsupial: opossum) anatomical configuration and compared them with our species, which has a derived eutherian configuration. In these models, nodes represent anatomical units and links represent their physical connection. Here, we aimed to identify: (1) the commonalities and differences in modularity between species, (2) whether modules present a potential phylogenetic character, and (3) whether modules preferentially reflect either developmental or functional aspects of anatomy, or a mix of both. We predicted differences between networks of metatherian and eutherian mammals that would best be explained by functional constraints, versus by constraints of development and/or phylogeny. The topology of contacts between bones, muscles, and bones + muscles showed that, among all three species, skeletal networks were more similar than musculoskeletal networks. There was no clear indication that humans and mice are more alike when compared to the opossum overall, even though their musculoskeletal and skeletal networks of fore- and hindlimbs are slightly more similar. Differences were greatest among musculoskeletal networks of heads and next of forelimbs, which showed more variation than hindlimbs, supporting previous anatomical studies indicating that in general the configuration of the hindlimbs changes less across evolutionary history. Most observations regarding the anatomical networks seem to be best explained by function, but an exception is the adult opossum ear ossicles. These ear bones might form an independent module because the incus and malleus are involved in forming a functional primary jaw that enables the neonate to attach to the teat, where this newborn will complete its development. Additionally, the human data show a specialized digit 1 module (thumb/big toe) in both limb types, likely the result of functional and evolutionary pressures, as our ape ancestors had highly movable big toes and thumbs.

Ontogenetic origins of cranial convergence between the extinct marsupial thylacine and placental gray wolf

Phenotypic convergence, describing the independent evolution of similar characteristics, offers unique insights into how natural selection influences developmental and molecular processes to generate shared adaptations. The extinct marsupial thylacine and placental gray wolf represent one of the most extraordinary cases of convergent evolution in mammals, sharing striking cranial similarities despite 160 million years of independent evolution. We digitally reconstructed their cranial ontogeny from birth to adulthood to examine how and when convergence arises through patterns of allometry, mosaicism, modularity, and integration. We find the thylacine and wolf crania develop along nearly parallel growth trajectories, despite lineage-specific constraints and heterochrony in timing of ossification. These constraints were found to enforce distinct cranial modularity and integration patterns during development, which were unable to explain their adult convergence. Instead, we identify a developmental origin for their convergent cranial morphologies through patterns of mosaic evolution, occurring within bone groups sharing conserved embryonic tissue origins. Interestingly, these patterns are accompanied by homoplasy in gene regulatory networks associated with neural crest cells, critical for skull patterning. Together, our findings establish empirical links between adaptive phenotypic and genotypic convergence and provides a digital resource for further investigations into the developmental basis of mammalian evolution.

Phylogeny and foraging behaviour shape modular morphological variation in bat humeri

Bats show a remarkable ecological diversity that is reflected both in dietary and foraging guilds (FGs). Cranial ecomorphological adaptations linked to diet have been widely studied in bats, using a variety of anatomical, computational and mathematical approaches. However, foraging-related ecomorphological adaptations and the concordance between cranial and postcranial morphological adaptations remain unexamined in bats and limited to the interpretation of traditional aerodynamic properties of the wing (e.g. wing loading [WL] and aspect ratio [AR]). For this reason, the postcranial ecomorphological diversity in bats and its drivers remain understudied. Using 3D virtual modelling and geometric morphometrics (GMM), we explored the phylogenetic, ecological and biological drivers of humeral morphology in bats, evaluating the presence and magnitude of modularity and integration. To explore decoupled patterns of variation across the bone, we analysed whole-bone shape, diaphyseal and epiphyseal shape. We also tested whether traditional aerodynamic wing traits correlate with humeral shape. By studying 37 species from 20 families (covering all FGs and 85% of dietary guilds), we found similar patterns of variation in whole-bone and diaphyseal shape and unique variation patterns in epiphyseal shape. Phylogeny, diet and FG significantly correlated with shape variation at all levels, whereas size only had a significant effect on epiphyseal morphology. We found a significant phylogenetic signal in all levels of humeral shape. Epiphyseal shape significantly correlated with wing AR. Statistical support for a diaphyseal-epiphyseal modular partition of the humerus suggests a functional partition of shape variability. Our study is the first to show within-structure modular morphological variation in the appendicular skeleton of any living tetrapod. Our results suggest that diaphyseal shape correlates more with phylogeny, whereas epiphyseal shape correlates with diet and FG.

Evolutionary and ontogenetic changes of the anatomical organization and modularity in the skull of archosaurs

Comparative anatomy studies of the skull of archosaurs provide insights on the mechanisms of evolution for the morphologically and functionally diverse species of crocodiles and birds. One of the key attributes of skull evolution is the anatomical changes associated with the physical arrangement of cranial bones. Here, we compare the changes in anatomical organization and modularity of the skull of extinct and extant archosaurs using an Anatomical Network Analysis approach. We show that the number of bones, their topological arrangement, and modular organization can discriminate birds from non-avian dinosaurs, and crurotarsans. We could also discriminate extant taxa from extinct species when adult birds were included. By comparing within the same framework, juveniles and adults for crown birds and alligator (Alligator mississippiensis), we find that adult and juvenile alligator skulls are topologically similar, whereas juvenile bird skulls have a morphological complexity and anisomerism more similar to those of non-avian dinosaurs and crurotarsans than of their own adult forms. Clade-specific ontogenetic differences in skull organization, such as extensive postnatal fusion of cranial bones in crown birds, can explain this pattern. The fact that juvenile and adult skulls in birds do share a similar anatomical integration suggests the presence of a specific constraint to their ontogenetic growth.

Conflict Resolution for Mesozoic Mammals: Reconciling Phylogenetic Incongruence Among Anatomical Regions

The evolutionary history of Mesozoic mammaliaformes is well studied. Although the backbone of their phylogeny is well resolved, the placement of ecologically specialized groups has remained uncertain. Functional and developmental covariation has long been identified as an important source of phylogenetic error, yet combining incongruent morphological characters altogether is currently a common practice when reconstructing phylogenetic relationships. Ignoring incongruence may inflate the confidence in reconstructing relationships, particularly for the placement of highly derived and ecologically specialized taxa, such as among australosphenidans (particularly, crown monotremes), haramiyidans, and multituberculates. The alternative placement of these highly derived clades can alter the taxonomic constituency and temporal origin of the mammalian crown group. Based on prior hypotheses and correlated homoplasy analyses, we identified cheek teeth and shoulder girdle character complexes as having a high potential to introduce phylogenetic error. We showed that incongruence among mandibulodental, cranial, and postcranial anatomical partitions for the placement of the australosphenidans, haramiyids, and multituberculates could largely be explained by apparently non-phylogenetic covariance from cheek teeth and shoulder girdle characters. Excluding these character complexes brought agreement between anatomical regions and improved the confidence in tree topology. These results emphasize the importance of considering and ameliorating major sources of bias in morphological data, and we anticipate that these will be valuable for confidently integrating morphological and molecular data in phylogenetic and dating analyses.

The embryo of the silky shrew opossum, Caenolestes fuliginosus (Tomes, 1863): First description of the embryo of Paucituberculata

The development of caenolestid marsupials (order Paucituberculata) is virtually unknown. We provide here the first description of Caenolestes fuliginosus embryos collected in the Colombian Central Andes. Our sample of four embryos comes from a single female caught during a fieldtrip at Río Blanco (Manizales, Caldas), in 2014. The sample was processed for macroscopic description using a Standard Event System and for histological descriptions (sectioning and staining). The grade of development of the lumbar flexure and coelomic closure differed between embryos, two of them being more advanced than the others (similar to McCrady's stages 30 and 29, respectively). The pericardial and peritoneal cavities were present, the hepatic anlage was organized in hepatic cords, the heart was in its final position, and the mesonephros was functional. Compared to other Neotropical marsupials, an early appearance of the frontonasal-maxillary fusion and the cervical growth (thickness) was observed; however, absorption of the pharyngeal arches into the body and lung development was delayed. Besides these differences, embryos were similar to equivalent stages in Didelphis virginiana and Monodelphis domestica. Previous proposals of litter size of four for C. fuliginosus are supported.

Evolutionary Integration and Modularity in the Archosaur Cranium

Complex structures, like the vertebrate skull, are composed of numerous elements or traits that must develop and evolve in a coordinated manner to achieve multiple functions. The strength of association among phenotypic traits (i.e., integration), and their organization into highly-correlated, semi-independent subunits termed modules, is a result of the pleiotropic and genetic correlations that generate traits. As such, patterns of integration and modularity are thought to be key factors constraining or facilitating the evolution of phenotypic disparity by influencing the patterns of variation upon which selection can act. It is often hypothesized that selection can reshape patterns of integration, parceling single structures into multiple modules or merging ancestrally semi-independent traits into a strongly correlated unit. However, evolutionary shifts in patterns of trait integration are seldom assessed in a unified quantitative framework. Here, we quantify patterns of evolutionary integration among regions of the archosaur skull to investigate whether patterns of cranial integration are conserved or variable across this diverse group. Using high-dimensional geometric morphometric data from 3D surface scans and computed tomography scans of modern birds (n = 352), fossil non-avian dinosaurs (n = 27), and modern and fossil mesoeucrocodylians (n = 38), we demonstrate that some aspects of cranial integration are conserved across these taxonomic groups, despite their major differences in cranial form, function, and development. All three groups are highly modular and consistently exhibit high integration within the occipital region. However, there are also substantial divergences in correlation patterns. Birds uniquely exhibit high correlation between the pterygoid and quadrate, components of the cranial kinesis apparatus, whereas the non-avian dinosaur quadrate is more closely associated with the jugal and quadratojugal. Mesoeucrocodylians exhibit a slightly more integrated facial skeleton overall than the other grades. Overall, patterns of trait integration are shown to be stable among archosaurs, which is surprising given the cranial diversity exhibited by the clade. At the same time, evolutionary innovations such as cranial kinesis that reorganize the structure and function of complex traits can result in modifications of trait correlations and modularity.

Detecting Mosaic Patterns in Macroevolutionary Disparity

Evolutionary biologists have long sought to understand the full complexity in pattern and process that shapes organismal diversity. Although phylogenetic comparative methods are often used to reconstruct complex evolutionary dynamics, they are typically limited to a single phenotypic trait. Extensions that accommodate multiple traits lack the ability to partition multidimensional data sets into a set of mosaic suites of evolutionarily linked characters. I introduce a comparative framework that identifies heterogeneity in evolutionary patterns across large data sets of continuous traits. Using a model of continuous trait evolution based on the differential accumulation of disparity across lineages in a phylogeny, the approach algorithmically partitions traits into a set of character suites that best explains the data, where each suite displays a distinct pattern in phylogenetic morphological disparity. When applied to empirical data, the approach revealed a mosaic pattern predicted by developmental biology. The evolutionary distinctiveness of individual suites can be investigated in more detail either by fitting conventional comparative models or by directly studying the phylogenetic patterns in disparity recovered during the analysis. This framework can supplement existing comparative approaches by inferring the complex, integrated patterns that shape evolution across the body plan from disparate developmental, morphometric, and environmental sources of phenotypic data.

Marsupial versus placental: assessing the evolutionary changes in the scapula of didelphids and sigmodontines

It is not a new concept that marsupials and placentals are distant and distinct clades among mammals. In South America, these animals coexist, occupy similar niches and, in some cases, are similar in appearance. This is especially true with respect to the locomotor categories of smaller rodents belonging to the family Cricetidae or, more specifically, the subfamily Sigmodontinae, compared with the marsupials of the Didelphidae family. In this study, we have investigated both the similarities and the differences between the two clades by examining locomotion-dependent adaptation, a crucial survival mechanism that has affected the morphology of both clades. We applied geometric morphometrics to quantify the shape of the scapula, which is a very adaptable structure. We found similar morphological adaptations between the clades, especially with respect to adaptation to life in trees. Moreover, Didelphidae are influenced by phylogenetic history to a greater extent than Sigmodontinae with regard to variation of scapula shape and allometry. These differences can be explained by the greater degree of body size variation that exists within the Didelphidae. Didelphidae have an ancient evolutionary history in South America compared with the Sigmodontinae, which have undergone a very successful and rapid diversification more recently.

Modularity increases rate of floral evolution and adaptive success for functionally specialized pollination systems

Angiosperm flowers have diversified in adaptation to pollinators, but are also shaped by developmental and genetic histories. The relative importance of these factors in structuring floral diversity remains unknown. We assess the effects of development, function and evolutionary history by testing competing hypotheses on floral modularity and shape evolution in Merianieae (Melastomataceae). Merianieae are characterized by different pollinator selection regimes and a developmental constraint: tubular anthers adapted to specialized buzz-pollination. Our analyses of tomography-based 3-dimensional flower models show that pollinators selected for functional modules across developmental units and that patterns of floral modularity changed during pollinator shifts. Further, we show that modularity was crucial for Merianieae to overcome the constraint of their tubular anthers through increased rates of evolution in other flower parts. We conclude that modularity may be key to the adaptive success of functionally specialized pollination systems by making flowers flexible (evolvable) for adaptation to changing selection regimes.

Dellinger et al. report the characterization of different pollinator selection regimes on 3D flower models. They show that pollinators selected for functional floral modules, and that this modularity allowed certain species to overcome morphological constraints by increased evolutionary rates in other flower parts.

Individual variation of the masticatory system dominates 3D skull shape in the herbivory-adapted marsupial wombats

Background

Within-species skull shape variation of marsupial mammals is widely considered low and strongly size-dependent (allometric), possibly due to developmental constraints arising from the altricial birth of marsupials. However, species whose skulls are impacted by strong muscular stresses – particularly those produced through mastication of tough food items – may not display such intrinsic patterns very clearly because of the known plastic response of bone to muscle activity of the individual. In such cases, allometry may not dominate within-species shape variation, even if it is a driver of evolutionary shape divergence; ordination of shape in a geometric morphometric context through principal component analysis (PCA) should reveal main variation in areas under masticatory stress (incisor region/zygomatic arches/mandibular ramus); but this main variation should emerge from high individual variability and thus have low eigenvalues.

Results

We assessed the evidence for high individual variation through 3D geometric morphometric shape analysis of crania and mandibles of three species of grazing-specialized wombats, whose diet of tough grasses puts considerable strain on their masticatory system. As expected, we found little allometry and low Principal Component 1 (PC1) eigenvalues within crania and mandibles of all three species. Also as expected, the main variation was in the muzzle, zygomatic arches, and masticatory muscle attachments of the mandibular ramus. We then implemented a new test to ask if the landmark variation reflected on PC1 was reflected in individuals with opposite PC1 scores and with opposite shapes in Procrustes space. This showed that correspondence between individual and ordinated shape variation was limited, indicating high levels of individual variability in the masticatory apparatus.

Discussion

Our results are inconsistent with hypotheses that skull shape variation within marsupial species reflects a constraint pattern. Rather, they support suggestions that individual plasticity can be an important determinant of within-species shape variation in marsupials (and possibly other mammals) with high masticatory stresses, making it difficult to understand the degree to which intrinsic constraints act on shape variation at the within-species level. We conclude that studies that link micro- and macroevolutionary patterns of shape variation might benefit from a focus on species with low-impact mastication, such as carnivorous or frugivorous species.

A tail of evolution: evaluating body length, weight and locomotion as potential drivers of tail length scaling in Australian marsupial mammals

Although mammalian tail length relative to body length is considered indicative of locomotor mode, this association has been difficult to quantify. This could be because the counterweight function of the tail might associate it more with body weight than body length. Alternatively, relative tail length might not be evolutionarily flexible owing to its integration with the remaining skeleton, particularly the spine. Using comparative analyses of morphological means and ranges in Australian marsupials, including the first co-assessment with body weight, our study supports the second hypothesis, i.e. tail length ranges within species, and tail lengths among species are explained better by body length than by body weight. However, all three variables do not differ in phylogenetic signal or rates of evolution. Associations of tail lengths with locomotion are limited, but suggest that scaling slopes, rather than intercepts, are responsible for limited divergence between relative tail lengths at different locomotor modes. This complicates (palaeo-)ecological interpretations of tail length further. We conclude that relative tail length is not a strong predictor of locomotor mode, probably owing to strong integration of tail and body length. The many well-documented bony and soft-tissue adaptations of tails are likely to be better suited to interpretations of locomotor adaptations.

Developmental plasticity associated with early structural integration and evolutionary patterns: Examples of developmental bias and developmental facilitation in the skeletal system

The relation of developmental plasticity to evolutionary diversification is a key component of evolutionary theory involving developmental bias, but the basis of the relationship varies among traits and among taxa. Here I review some scenarios of how structural integration during early organogenesis could influence this relationship. When condensations are highly integrated and dependent on each other during early organogenesis, both plasticity and evolution are restricted, for example size proportions in molar tooth rows and phalanges within a digit. When similar condensations develop and remain separate (in tracheal cartilages and feather buds), they show high levels of variation and diversity in number but not in shape and size, at least at early stages. When non-similar structures form separately and then integrate while still undergoing patterning, high levels of plasticity (in number, size, shape; in rib uncinate processes) or new dimensions of ecologically-significant variation (cusp offset, in mammal teeth) are seen. Although each of these structural integration scenarios is unique, the modulation of evolvability is detectable and informative. Parsing the influence of structural integration at these developmental levels, rather than later-stage structural correlations or only through genetic covariation, may be necessary to advance understanding of evolvability of the phenotype.

Postcranial heterochrony, modularity, integration and disparity in the prenatal ossification in bats (Chiroptera)

BACKGROUND

Self-powered flight is one of the most energy-intensive types of locomotion found in vertebrates. It is also associated with a range of extreme morpho-physiological adaptations that evolved independently in three different vertebrate groups. Considering that development acts as a bridge between the genotype and phenotype on which selection acts, studying the ossification of the postcranium can potentially illuminate our understanding of bat flight evolution. However, the ontogenetic basis of vertebrate flight remains largely understudied. Advances in quantitative analysis of sequence heterochrony and morphogenetic growth have created novel approaches to study the developmental basis of diversification and the evolvability of skeletal morphogenesis. Assessing the presence of ontogenetic disparity, integration and modularity from an evolutionary approach allows assessing whether flight may have resulted in evolutionary differences in the magnitude and mode of development in bats.

RESULTS

We quantitatively compared the prenatal ossification of the postcranium (24 bones) between bats (14 species), non-volant mammals (11 species) and birds (14 species), combining for the first time prenatal sequence heterochrony and developmental growth data. Sequence heterochrony was found across groups, showing that bat postcranial development shares patterns found in other flying vertebrates but also those in non-volant mammals. In bats, modularity was found as an axial-appendicular partition, resembling a mammalian pattern of developmental modularity and suggesting flight did not repattern prenatal postcranial covariance in bats.

CONCLUSIONS

Combining prenatal data from 14 bat species, this study represents the most comprehensive quantitative analysis of chiropteran ossification to date. Heterochrony between the wing and leg in bats could reflect functional needs of the newborn, rather than ecological aspects of the adult. Bats share similarities with birds in the development of structures involved in flight (i.e. handwing and sternum), suggesting that flight altriciality and early ossification of pedal phalanges and sternum are common across flying vertebrates. These results indicate that the developmental modularity found in bats facilitates intramodular phenotypic diversification of the skeleton. Integration and disparity increased across developmental time in bats. We also found a delay in the ossification of highly adaptable and evolvable regions (e.g. handwing and sternum) that are directly associated with flight performance.

A fly in a tube: Macroevolutionary expectations for integrated phenotypes

Phenotypic integration and modularity are ubiquitous features of complex organisms, describing the magnitude and pattern of relationships among biological traits. A key prediction is that these relationships, reflecting genetic, developmental, and functional interactions, shape evolutionary processes by governing evolvability and constraint. Over the last 60 years, a rich literature of research has quantified patterns of integration and modularity across a variety of clades and systems. Only recently has it become possible to contextualize these findings in a phylogenetic framework to understand how trait integration interacts with evolutionary tempo and mode. Here, we review the state of macroevolutionary studies of integration and modularity, synthesizing empirical and theoretical work into a conceptual framework for predicting the effects of integration on evolutionary rate and disparity: a fly in a tube. While magnitude of integration is expected to influence the potential for phenotypic variation to be produced and maintained, thus defining the shape and size of a tube in morphospace, evolutionary rate, or the speed at which a fly moves around the tube, is not necessarily controlled by trait interactions. Finally, we demonstrate this reduced disparity relative to the Brownian expectation for a given rate of evolution with an empirical example: the avian cranium.

Morphological integration in the appendicular skeleton of two domestic taxa: the horse and donkey

Organisms are organized into suites of anatomical structures that typically covary when developmentally or functionally related, and this morphological integration plays a determinant role in evolutionary processes. Artificial selection on domestic species causes strong morphological changes over short time spans, frequently resulting in a wide and exaggerated phenotypic diversity. This raises the question of whether integration constrains the morphological diversification of domestic species and how natural and artificial selection may impact integration patterns. Here, we study the morphological integration in the appendicular skeleton of domestic horses and donkeys, using three-dimensional geometric morphometrics on 75 skeletons. Our results indicate that a strong integration is inherited from developmental mechanisms which interact with functional factors. This strong integration reveals a specialization in the locomotion of domestic equids, partly for running abilities. We show that the integration is stronger in horses than in donkeys, probably because of a greater degree of specialization and predictability of their locomotion. Thus, the constraints imposed by integration are weak enough to allow important morphological changes and the phenotypic diversification of domestic species.

Methods for Studying Morphological Integration and Modularity

Morphological integration and modularity are closely related concepts about how different traits of an organism are correlated. Integration is the overall pattern of intercorrelation; modularity is the partitioning of integration into evolutionarily or developmentally independent blocks of traits. Modularity and integration are usually studied using quantitative phenotypic data, which can be obtained either from extant or fossil organisms. Many methods are now available to study integration and modularity, all of which involve the analysis of patterns found in trait correlation or covariance matrices. We review matrix correlation, random skewers, fluctuating asymmetry, cluster analysis, Euclidean distance matrix analysis (EDMA), graphical modelling, two-block partial least squares, RV coefficients, and theoretical matrix modelling and discuss their similarities and differences. We also review different coefficients that are used to measure correlations. We apply all the methods to cranial landmark data from and ontogenetic series of Japanese macaques,Macaca fuscatato illustrate the methods and their individual strengths and weaknesses. We conclude that the exploratory approaches (cluster analyses of various sorts) were less informative and less consistent with one another than were the results of model testing or comparative approaches. Nevertheless, we found that competing models of modularity and integration are often similar enough that they are not statistically distinguishable; we expect, therefore, that several models will often be significantly correlated with observed data.

Non-decoupled morphological evolution of the fore- and hindlimb of sabretooth predators

Specialized organisms are useful for exploring the combined effects of selection of functional traits and developmental constraints on patterns of phenotypic integration. Sabretooth predators are one of the most interesting examples of specialization among mammals. Their hypertrophied, sabre-shaped upper canines and their powerfully built forelimbs have been interpreted as adaptations to a highly specialized predatory behaviour. Given that the elongated and laterally compressed canines of sabretooths were more vulnerable to fracture than the shorter canines of conical-tooth cats, it has been long hypothesized that the heavily muscled forelimbs of sabretooths were used for immobilizing prey before developing a quick and precise killing bite. However, the effect of this unique adaptation on the covariation between the fore- and the hindlimb has not been explored in a quantitative fashion. In this paper, we investigate if the specialization of sabretooth predators decoupled the morphological variation of their forelimb with respect to their hindlimb or, in contrast, both limbs vary in the same fashion as in conical-tooth cats, which do not show such extreme adaptations in their forelimb. We use 3D geometric morphometrics and different morphological indices to compare the fore- and hindlimb of conical- and sabretooth predators. Our results indicate that the limb bones of sabretooth predators covary following the same trend of conical-tooth cats. Therefore, we show that the predatory specialization of sabretooth predators did not result in a decoupling of the morphological evolution of their fore- and hindlimbs. The role of developmental constraints and natural selection on this coordinate variation between the fore- and the hindlimb is discussed in the light of this new evidence.

Patterns in the bony skull development of marsupials: high variation in onset of ossification and conserved regions of bone contact

Development in marsupials is specialized towards an extremely short gestation and highly altricial newborns. As a result, marsupial neonates display morphological adaptations at birth related to functional constraints. However, little is known about the variability of marsupial skull development and its relation to morphological diversity. We studied bony skull development in five marsupial species. The relative timing of the onset of ossification was compared to literature data and the ossification sequence of the marsupial ancestor was reconstructed using squared-change parsimony. The high range of variation in the onset of ossification meant that no patterns could be observed that differentiate species. This finding challenges traditional studies concentrating on the onset of ossification as a marker for phylogeny or as a functional proxy. Our study presents observations on the developmental timing of cranial bone-to-bone contacts and their evolutionary implications. Although certain bone contacts display high levels of variation, connections of early and late development are quite conserved and informative. Bones that surround the oral cavity are generally the first to connect and the bones of the occipital region are among the last. We conclude that bone contact is preferable over onset of ossification for studying cranial bone development.

Out on a limb: bandicoot limb co-variation suggests complex impacts of development and adaptation on marsupial forelimb evolution

Marsupials display far less forelimb diversity than placentals, possibly because of the laborious forelimb-powered climb to the pouch performed by most marsupial neonates. This is thought to result in stronger morphological integration (i.e., higher co-variance) within the marsupial forelimb skeleton, and lower integration between marsupial fore- and hind limbs, compared to other mammals. Possible mechanisms for this constraint are a fundamental developmental change in marsupial limb patterning, or alternatively more immediate perinatal biomechanical and metabolic requirements. In the latter case, peramelid marsupials (bandicoots), which have neonates that climb very little, should show lower within-limb and higher between-limb integration, compared to other marsupials. We tested this in four peramelid species and the related bilby, using partial correlation analyses of between-landmark linear measurements of limb bones, and Procrustes-based two-block partial least-squares analysis (2B-PLS) of limb bone shapes using the same landmarks. We find extensive between-limb integration in partial correlation analyses of only bone lengths, consistent with a reduction of a short-term biomechanical/allocation constraint in peramelid forelimbs. However, partial correlations of bone proportions and 2B-PLS reveal extensive shape divergence between correlated bone pairs. This result contradicts expectations of developmental constraints or serial homology, instead suggesting a function-driven integration pattern. Comparing visualizations from cross-species principal components analysis and 2B-PLS, we tentatively identify selection for digging and half-bounding as the main driver of bandicoot limb integration patterning. This calls for further assessments of functional versus developmental limb integration in marsupials with a more strenuous neonatal climb to the pouch.

EMMLi: A maximum likelihood approach to the analysis of modularity

Identification of phenotypic modules, semiautonomous sets of highly correlated traits, can be accomplished through exploratory (e.g., cluster analysis) or confirmatory approaches (e.g., RV coefficient analysis). Although statistically more robust, confirmatory approaches are generally unable to compare across different model structures. For example, RV coefficient analysis finds support for both two- and six-module models for the therian mammalian skull. Here, we present a maximum likelihood approach that takes into account model parameterization. We compare model log-likelihoods of trait correlation matrices using the finite-sample corrected Akaike Information Criterion, allowing for comparison of hypotheses across different model structures. Simulations varying model complexity and within- and between-module contrast demonstrate that this method correctly identifies model structure and parameters across a wide range of conditions. We further analyzed a dataset of 3-D data, consisting of 61 landmarks from 181 macaque (Macaca fuscata) skulls, distributed among five age categories, testing 31 models, including no modularity among the landmarks and various partitions of two, three, six, and eight modules. Our results clearly support a complex six-module model, with separate within- and intermodule correlations. Furthermore, this model was selected for all five age categories, demonstrating that this complex pattern of integration in the macaque skull appears early and is highly conserved throughout postnatal ontogeny. Subsampling analyses demonstrate that this method is robust to relatively low sample sizes, as is commonly encountered in rare or extinct taxa. This new approach allows for the direct comparison of models with different parameterizations, providing an important tool for the analysis of modularity across diverse systems.

Do Developmental Constraints and High Integration Limit the Evolution of the Marsupial Oral Apparatus?

Developmental constraints can have significant influence on the magnitude and direction of evolutionary change, and many studies have demonstrated that these effects are manifested on macroevolutionary scales. Phenotypic integration, or the strong interactions among traits, has been similarly invoked as a major influence on morphological variation, and many studies have demonstrated that trait integration changes through ontogeny, in many cases decreasing with age. Here, we unify these perspectives in a case study of the ontogeny of the mammalian cranium, focusing on a comparison between marsupials and placentals. Marsupials are born at an extremely altricial state, requiring, in most cases, the use of the forelimbs to climb to the pouch, and, in all cases, an extended period of continuous suckling, during which most of their development occurs. Previous work has shown that marsupials are less disparate in adult cranial form than are placentals, particularly in the oral apparatus, and in forelimb ontogeny and adult morphology, presumably due to functional selection pressures on these two systems during early postnatal development. Using phenotypic trajectory analysis to quantify prenatal and early postnatal cranial ontogeny in 10 species of therian mammals, we demonstrate that this pattern of limited variation is also apparent in the development of the oral apparatus of marsupials, relative to placentals, but not in the skull more generally. Combined with the observation that marsupials show extremely high integration of the oral apparatus in early postnatal ontogeny, while other cranial regions show similar levels of integration to that observed in placentals, we suggest that high integration may compound the effects of the functional constraints for continuous suckling to ultimately limit the ontogenetic and adult disparity of the marsupial oral apparatus throughout their evolutionary history.

Goldfish morphology as a model for evolutionary developmental biology

Morphological variation of the goldfish is known to have been established by artificial selection for ornamental purposes during the domestication process. Chinese texts that date to the Song dynasty contain descriptions of goldfish breeding for ornamental purposes, indicating that the practice originated over one thousand years ago. Such a well-documented goldfish breeding process, combined with the phylogenetic and embryological proximities of this species with zebrafish, would appear to make the morphologically diverse goldfish strains suitable models for evolutionary developmental (evodevo) studies. However, few modern evodevo studies of goldfish have been conducted. In this review, we provide an overview of the historical background of goldfish breeding, and the differences between this teleost and zebrafish from an evolutionary perspective. We also summarize recent progress in the field of molecular developmental genetics, with a particular focus on the twin-tail goldfish morphology. Furthermore, we discuss unanswered questions relating to the evolution of the genome, developmental robustness, and morphologies in the goldfish lineage, with the goal of blazing a path toward an evodevo study paradigm using this teleost species as a new model species. For further resources related to this article, please visit the WIREs website.

Genetic recapitulation of human pre-eclampsia risk during convergent evolution of reduced placental invasiveness in eutherian mammals

The relationship between phenotypic variation arising through individual development and phenotypic variation arising through diversification of species has long been a central question in evolutionary biology. Among humans, reduced placental invasion into endometrial tissues is associated with diseases of pregnancy, especially pre-eclampsia, and reduced placental invasiveness has also evolved, convergently, in at least 10 lineages of eutherian mammals. We tested the hypothesis that a common genetic basis underlies both reduced placental invasion arising through a developmental process in human placental disease and reduced placental invasion found as a derived trait in the diversification of Euarchontoglires (rodents, lagomorphs, tree shrews, colugos and primates). Based on whole-genome analyses across 18 taxa, we identified 1254 genes as having evolved adaptively across all three lineages exhibiting independent evolutionary transitions towards reduced placental invasion. These genes showed strong evidence of enrichment for associations with pre-eclampsia, based on genetic-association studies, gene-expression analyses and gene ontology. We further used in silico prediction to identify a subset of 199 genes that are likely targets of natural selection during transitions in placental invasiveness and which are predicted to also underlie human placental disorders. Our results indicate that abnormal ontogenies can recapitulate major phylogenetic shifts in mammalian evolution, identify new candidate genes for involvement in pre-eclampsia, imply that study of species with less-invasive placentation will provide useful insights into the regulation of placental invasion and pre-eclampsia, and recommend a novel comparative functional-evolutionary approach to the study of genetically based human disease and mammalian diversification.

The fossil record of phenotypic integration and modularity: A deep-time perspective on developmental and evolutionary dynamics

Variation is the raw material for natural selection, but the factors shaping variation are still poorly understood. Genetic and developmental interactions can direct variation, but there has been little synthesis of these effects with the extrinsic factors that can shape biodiversity over large scales. The study of phenotypic integration and modularity has the capacity to unify these aspects of evolutionary study by estimating genetic and developmental interactions through the quantitative analysis of morphology, allowing for combined assessment of intrinsic and extrinsic effects. Data from the fossil record in particular are central to our understanding of phenotypic integration and modularity because they provide the only information on deep-time developmental and evolutionary dynamics, including trends in trait relationships and their role in shaping organismal diversity. Here, we demonstrate the important perspective on phenotypic integration provided by the fossil record with a study of Smilodon fatalis (saber-toothed cats) and Canis dirus (dire wolves). We quantified temporal trends in size, variance, phenotypic integration, and direct developmental integration (fluctuating asymmetry) through 27,000 y of Late Pleistocene climate change. Both S. fatalis and C. dirus showed a gradual decrease in magnitude of phenotypic integration and an increase in variance and the correlation between fluctuating asymmetry and overall integration through time, suggesting that developmental integration mediated morphological response to environmental change in the later populations of these species. These results are consistent with experimental studies and represent, to our knowledge, the first deep-time validation of the importance of developmental integration in stabilizing morphological evolution through periods of environmental change.

The macroevolutionary consequences of phenotypic integration: from development to deep time

Phenotypic integration is a pervasive characteristic of organisms. Numerous analyses have demonstrated that patterns of phenotypic integration are conserved across large clades, but that significant variation also exists. For example, heterochronic shifts related to different mammalian reproductive strategies are reflected in postcranial skeletal integration and in coordination of bone ossification. Phenotypic integration and modularity have been hypothesized to shape morphological evolution, and we extended simulations to confirm that trait integration can influence both the trajectory and magnitude of response to selection. We further demonstrate that phenotypic integration can produce both more and less disparate organisms than would be expected under random walk models by repartitioning variance in preferred directions. This effect can also be expected to favour homoplasy and convergent evolution. New empirical analyses of the carnivoran cranium show that rates of evolution, in contrast, are not strongly influenced by phenotypic integration and show little relationship to morphological disparity, suggesting that phenotypic integration may shape the direction of evolutionary change, but not necessarily the speed of it. Nonetheless, phenotypic integration is problematic for morphological clocks and should be incorporated more widely into models that seek to accurately reconstruct both trait and organismal evolution.

Quantifying the impact of development on phenotypic variation and evolution

A primary goal of evolutionary biology is to identify the factors that shape phenotypic evolution. According to the theory of natural selection, phenotypic evolution occurs through the differential survival and reproduction of individuals whose traits are selectively advantageous relative to other individuals in the population. This implies that evolution by natural selection is contingent upon the distribution and magnitude of phenotypic variation among individuals, which are in turn the products of developmental processes. Development therefore has the potential to affect the trajectory and rate of phenotypic evolution. Recent research in diverse systems (e.g., mammalian teeth, cichlid skulls, butterfly wings, and marsupial limbs) supports the hypothesis that development biases phenotypic variation and evolution, but suggests that these biases might be system-specific.

Heterochrony, modularity, and the functional evolution of the mechanosensory lateral line canal system of fishes

Background

The canals of the mechanosensory lateral line system are components of the dermatocranium, and demonstrate phenotypic variation in bony fishes. Widened lateral line canals evolved convergently in a limited number of families of teleost fishes and it had been hypothesized that they evolved from narrow canals via heterochrony and explore modularity in the lateral line system. Two species of cichlids with different canal phenotypes were used to test a hypothesis of heterochrony. Histological material prepared from ontogenetic series of Aulonocara stuartgranti (widened canals) and Tramitichromis sp. (narrow canals) was analyzed using ANCOVA to determine rates of increase in canal diameter and neuromast size (length, width) and to compare the timing of onset of critical stages in canal morphogenesis (enclosure, ossification).

Results

A faster rate of increase in canal diameter and neuromast width (but not length), and a delay in onset of canal morphogenesis were found in Aulonocara relative to Tramitichromis. However, rates of increase in canal diameter and neuromast size among canals, among canal portions and among canals segments reveal similar trends within both species.

Conclusion

The evolution of widened lateral line canals is the result of dissociated heterochrony - acceleration in the rate of increase of both canal diameter and neuromast size, and delay in the onset of canal morphogenesis, in Aulonocara (widened canals) relative to Tramitichromis (narrow canals). Common rates of increase in canal diameter and neuromast size among canal portions in different dermatocranial bones and among canal segments reflect the absence of local heterochronies, and suggest modular integration among canals in each species. Thus, canal and neuromast morphology are more strongly influenced by their identities as features of the lateral line system than by the attributes of the dermatocranial bones in which the canals are found. Rate heterochrony manifested during the larval stage ensures that the widened canal phenotype, known to be associated with benthic prey detection in adult Aulonocara, is already present before feeding commences. Heterochrony can likely explain the convergent evolution of widened lateral line canals among diverse taxa. The lateral line system provides a valuable context for novel analyses of the relationship between developmental processes and the evolution of behaviorally and ecologically relevant phenotypes in fishes.

Patterns of postcranial ossification and sequence heterochrony in bats: life histories and developmental trade-offs

The recently increased interest in studies on sequence heterochrony has uncovered developmental variation between species. However, how changes in developmental program are related to shifts in life-history parameters remains largely unsolved. Here we provide the most comprehensive data to date on postcranial ossification sequence of bats and compare them to various boreoeutherian mammals with different locomotive modes. Given that bats are equipped with an elongated manus, we expected to detect characteristic heterochronies particularly related to wing development. Although heterochronies related to wing development were confirmed as predicted, unexpected heterochronies regarding the pedal digits were also found. The timing of ossification onset of pedal phalanges is earlier than other mammals. Particularly, bats deviate from others in that pedal phalanges initiate ossification earlier than manual phalanges. It is known that the foot size of new born bats is close to that of adults, and that it takes several weeks to month until the wing is developed for flight. Given that the foot is required to be firm and stable enough at the time of birth to allow continued attachment to the mother and/or cave walls, we suggest that the accelerated development of the hind foot is linked to their unique life history. Since the forelimb is not mature enough for flight at birth and requires extended postnatal time to be large enough to be fully functional, we postulate that bats invest in earlier development of the hindlimb. We conclud that energy allocation trade-offs can play a significant role in shaping the evolution of development.

Morphological integration in the forelimb of musteloid carnivorans

The forelimb forms a functional unit that allows a variety of behaviours and needs to be mobile, yet at the same time stable. Both mobility and stability are controlled, amongst others, at the level of the elbow joint. This joint is composed of the humero-ulnar articulation, mainly involved during parasagittal movements; and the radio-ulnar articulation, mainly allowing rotation. In contrast, the humero-radial articulation allows both movements of flexion-extension and rotation. Here, we study the morphological integration between each bone of the forelimb at the level of the entire arm, as well as at the elbow joint, in musteloid carnivorans. To do so, we quantitatively test shape co-variation using surface 3D geometric morphometric data. Our results show that morphological integration is stronger for bones that form functional units. Different results are obtained depending on the level of investigation: for the entire arm, results show a greater degree of shape co-variation between long bones of the lower arm than between the humerus and either bone of the lower arm. Thus, at this level the functional unit of the lower arm is comprised of the radius and ulna, permitting rotational movements of the lower arm. At the level of the elbow, results display a stronger shape co-variation between bones allowing flexion and stability (humerus and ulna) than between bones allowing mobility (ulna and radius and humerus and radius). Thus, the critical functional unit appears to be the articulation between the humerus and ulna providing the stability of the joint.

Mammalian skull heterochrony reveals modular evolution and a link between cranial development and brain size

The multiple skeletal components of the skull originate asynchronously and their developmental schedule varies across amniotes. Here we present the embryonic ossification sequence of 134 species, covering all major groups of mammals and their close relatives. This comprehensive data set allows reconstruction of the heterochronic and modular evolution of the skull and the condition of the last common ancestor of mammals. We show that the mode of ossification (dermal or endochondral) unites bones into integrated evolutionary modules of heterochronic changes and imposes evolutionary constraints on cranial heterochrony. However, some skull-roof bones, such as the supraoccipital, exhibit evolutionary degrees of freedom in these constraints. Ossification timing of the neurocranium was considerably accelerated during the origin of mammals. Furthermore, association between developmental timing of the supraoccipital and brain size was identified among amniotes. We argue that cranial heterochrony in mammals has occurred in concert with encephalization but within a conserved modular organization.