High evolutionary constraints limited adaptive responses to past climate changes in toad skulls

Interactions among traits that build a complex structure may be represented as genetic covariation and correlation. Genetic correlations may act as constraints, deflecting the evolutionary response from the direction of natural selection. We investigated the relative importance of drift, selection, and constraints in driving skull divergence in a group of related toad species. The distributional range of these species encompasses very distinct habitats with important climatic differences and the species are primarily distinguished by differences in their skulls. Some parts of the toad skull, such as the snout, may have functional relevance in reproductive ecology, detecting water cues. Thus, we hypothesized that the species skull divergence was driven by natural selection associated with climatic variation. However, given that all species present high correlations among skull traits, our second prediction was of high constraints deflecting the response to selection. We first extracted the main morphological direction that is expected to be subjected to selection by using within- and between-species covariance matrices. We then used evolutionary regressions to investigate whether divergence along this direction is explained by climatic variation between species. We also used quantitative genetics models to test for a role of random drift versus natural selection in skull divergence and to reconstruct selection gradients along species phylogeny. Climatic variables explained high proportions of between-species variation in the most selected axis. However, most evolutionary responses were not in the direction of selection, but aligned with the direction of allometric size, the dimension of highest phenotypic variance in the ancestral population. We conclude that toad species have responded to selection related to climate in their skulls, yet high evolutionary constraints dominated species divergence and may limit species responses to future climate change.     

Cranial modularity shifts during mammalian evolution

The mammalian skull has been studied as several separate functional components for decades, but the study of modularity is a more recent, integrative approach toward quantitative examination of independent subsets of highly correlated traits, or modules. Although most studies of modularity focus on developmental and genetic systems, phenotypic modules have been noted in many diverse morphological structures. However, few studies have provided empirical data for comparing modules across higher taxonomic levels, limiting the ability to assess the broader evolutionary significance of modularity. This study uses 18-32 three-dimensional cranial landmarks to analyze phenotypic modularity in 106 mammalian species and demonstrates that cranial modularity is generally conserved in the evolution of therian mammals (marsupials and placentals) but differs between therians and monotremes, the two extant subclasses of Mammalia. Within therians, cluster analyses identify six distinct modules, but only three modules display significant integration in all species. Monotremes display only two highly integrated modules. Specific hypotheses of functional and developmental influences on cranial bones were tested. Theoretical correlation matrices for bones were constructed on the basis of shared function, tissue origin, or mode of ossification, and all three of these models are significantly correlated with observed correlation matrices for the mammalian cranium.     

SIZE VARIATION,GROWTH STRATEGIES,AND THE EVOLUTION OF MODULARITY IN THE MAMMALIAN SKULL

Allometry is a major determinant of within‐population patterns of association among traits and, therefore, a major component of morphological integration studies. Even so, the influence of size variation over evolutionary change has been largely unappreciated. Here, we explore the interplay between allometric size variation, modularity, and life‐history strategies in the skull from representatives of 35 mammalian families. We start by removing size variation from within‐species data and analyzing its influence on integration magnitudes, modularity patterns, and responses to selection. We also carry out a simulation in which we artificially alter the influence of size variation in within‐taxa matrices. Finally, we explore the relationship between size variation and different growth strategies. We demonstrate that a large portion of the evolution of modularity in the mammalian skull is associated to the evolution of growth strategies. Lineages with highly altricial neonates have adult variation patterns dominated by size variation, leading to high correlations among traits regardless of any underlying modular process and impacting directly their potential to respond to selection. Greater influence of size variation is associated to larger intermodule correlations, less individualized modules, and less flexible responses to natural selection.     

Why the short face? Developmental disintegration of the neurocranium drives convergent evolution in neotropical electric fishes

Convergent evolution is widely viewed as strong evidence for the influence of natural selection on the origin of phenotypic design. However, the emerging evo‐devo synthesis has highlighted other processes that may bias and direct phenotypic evolution in the presence of environmental and genetic variation. Developmental biases on the production of phenotypic variation may channel the evolution of convergent forms by limiting the range of phenotypes produced during ontogeny. Here, we study the evolution and convergence of brachycephalic and dolichocephalic skull shapes among 133 species of Neotropical electric fishes (Gymnotiformes: Teleostei) and identify potential developmental biases on phenotypic evolution. We plot the ontogenetic trajectories of neurocranial phenotypes in 17 species and document developmental modularity between the face and braincase regions of the skull. We recover a significant relationship between developmental covariation and relative skull length and a significant relationship between developmental covariation and ontogenetic disparity. We demonstrate that modularity and integration bias the production of phenotypes along the brachycephalic and dolichocephalic skull axis and contribute to multiple, independent evolutionary transformations to highly brachycephalic and dolichocephalic skull morphologies.     

Roles for modularity and constraint in the evolution of cranial diversity among Anolis lizards

Complex organismal structures are organized into modules, suites of traits that develop, function, and vary in a coordinated fashion. By limiting or directing covariation among component traits, modules are expected to represent evolutionary building blocks and to play an important role in morphological diversification. But how stable are patterns of modularity over macroevolutionary timescales? Comparative analyses are needed to address the macroevolutionary effect of modularity, but to date few have been conducted. We describe patterns of skull diversity and modularity in Caribbean Anolis lizards. We first diagnose the primary axes of variation in skull shape and then examine whether diversification of skull shape is concentrated to changes within modules or whether changes arose across the structure as a whole. We find no support for the hypothesis that cranial modules are conserved as species diversify in overall skull shape. Instead we find that anole skull shape and modularity patterns independently converge. In anoles, skull modularity is evolutionarily labile and may reflect the functional demands of unique skull shapes. Our results suggest that constraints have played little role in limiting or directing the diversification of head shape in Anolis lizards.     

Anatomical Network Analysis Shows Decoupling of Modular Lability and Complexity in the Evolution of the Primate Skull

Modularity and complexity go hand in hand in the evolution of the skull of primates. Because analyses of these two parameters often use different approaches, we do not know yet how modularity evolves within, or as a consequence of, an also-evolving complex organization. Here we use a novel network theory-based approach (Anatomical Network Analysis) to assess how the organization of skull bones constrains the co-evolution of modularity and complexity among primates. We used the pattern of bone contacts modeled as networks to identify connectivity modules and quantify morphological complexity. We analyzed whether modularity and complexity evolved coordinately in the skull of primates. Specifically, we tested Herbert Simon’s general theory of near-decomposability, which states that modularity promotes the evolution of complexity. We found that the skulls of extant primates divide into one conserved cranial module and up to three labile facial modules, whose composition varies among primates. Despite changes in modularity, statistical analyses reject a positive feedback between modularity and complexity. Our results suggest a decoupling of complexity and modularity that translates to varying levels of constraint on the morphological evolvability of the primate skull. This study has methodological and conceptual implications for grasping the constraints that underlie the developmental and functional integration of the skull of humans and other primates.     

The Evolution of Modularity in the Mammalian Skull II: Evolutionary Consequences

Changes in patterns and magnitudes of integration may influence the ability of a species to respond to selection. Consequently, modularity has often been linked to the concept of evolvability, but their relationship has rarely been tested empirically. One possible explanation is the lack of analytical tools to compare patterns and magnitudes of integration among diverse groups that explicitly relate these aspects to the quantitative genetics framework. We apply such framework here using the multivariate response to selection equation to simulate the evolutionary behavior of several mammalian orders in terms of their flexibility, evolvability and constraints in the skull. We interpreted these simulation results in light of the integration patterns and magnitudes of the same mammalian groups, described in a companion paper. We found that larger magnitudes of integration were associated with a blur of the modules in the skull and to larger portions of the total variation explained by size variation, which in turn can exert a strong evolutionary constraint, thus decreasing the evolutionary flexibility. Conversely, lower overall magnitudes of integration were associated with distinct modules in the skull, to smaller fraction of the total variation associated with size and, consequently, to weaker constraints and more evolutionary flexibility. Flexibility and constraints are, therefore, two sides of the same coin and we found them to be quite variable among mammals. Neither the overall magnitude of morphological integration, the modularity itself, nor its consequences in terms of constraints and flexibility, were associated with absolute size of the organisms, but were strongly associated with the proportion of the total variation in skull morphology captured by size. Therefore, the history of the mammalian skull is marked by a trade-off between modularity and evolvability. Our data provide evidence that, despite the stasis in integration patterns, the plasticity in the magnitude of integration in the skull had important consequences in terms of evolutionary flexibility of the mammalian lineages.     

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.     

Skull ontogeny and modularity in two species of Lagenorhynchus: Morphological and ecological implications

Comparisons of skull shape between closely related species can provide information on the role that phylogeny and function play in cranial evolution. We used 3D‐anatomical landmarks in order to study the skull ontogeny of two closely related species, Lagenorhynchus obscurus and Lagenorhynchus australis , with a total sample of 52 skulls. We found shared trends between species, such as the relative compression of the neurocranium and the enlargement of the rostrum during ontogeny. However, these are common mammalian features, associated with prenatal brain development and sensory capsules. Moreover, we found a posterior displacement of the external nares and infraorbital foramina, and a strong development of the rostrum in an anteroposterior direction. Such trends are associated with the process of telescoping and have been observed in postnatal ontogeny of other odontocetes, suggesting a constraint in the pattern. Interspecific differences related to the deepness of facial region, robustness of the feeding apparatus and rostrum orientation may be related with the specific lifestyles of L. obscurus and L. australis . We also tested the presence of three different modules in the skull (basicranium, neurocranium, rostrum), all of which presented strong integration. Only the rostrum showed a different ontogenetic trajectory between species. Even though we detected directional asymmetry, changes in this feature along ontogeny were not detectable. Because asymmetry may be related to echolocation, our results suggest a functional importance of directional asymmetry from the beginning of postnatal life. J. Morphol. 278:203–214, 2017. © 2016 Wiley Periodicals,Inc.     

A MODULAR FRAMEWORK CHARACTERIZES MICRO‐ AND MACROEVOLUTION OF OLD WORLD MONKEY DENTITIONS

The study of modularity can provide a foundation for integrating development into studies of phenotypic evolution. The dentition is an ideal phenotype for this as it is developmentally relatively simple, adaptively highly significant, and evolutionarily tractable through the fossil record. Here, we use phenotypic variation in the dentition to test a hypothesis about genetic modularity. Quantitative genetic analysis of size variation in the baboon dentition indicates a genetic modular framework corresponding to tooth type categories. We analyzed covariation within the dentitions of six species of Old World monkeys (OWMs) to assess the macroevolutionary extent of this framework: first by estimating variance–covariance matrices of linear tooth size, and second by performing a geometric morphometric (GM) analysis of tooth row shape. For both size and shape, we observe across OWMs a framework of anterior and postcanine modules, as well as submodularity between the molars and premolars. Our results of modularity by tooth type suggest that adult variation in the OWM dentition is influenced by early developmental processes such as odontogenesis and jaw patterning. This study presents a comparison of genotypic modules to phenotypic modules, which can be used to better understand their action across evolutionary time scales.     

Ecological,historical and evolutionary determinants of modularity in weighted seed‐dispersal networks

Modularity is a recurrent and important property of bipartite ecological networks. Although well‐resolved ecological networks describe interaction frequencies between species pairs, modularity of bipartite networks has been analysed only on the basis of binary presence–absence data. We employ a new algorithm to detect modularity in weighted bipartite networks in a global analysis of avian seed‐dispersal networks. We define roles of species, such as connector values, for weighted and binary networks and associate them with avian species traits and phylogeny. The weighted, but not binary, analysis identified a positive relationship between climatic seasonality and modularity, whereas past climate stability and phylogenetic signal were only weakly related to modularity. Connector values were associated with foraging behaviour and were phylogenetically conserved. The weighted modularity analysis demonstrates the dominating impact of ecological factors on the structure of seed‐dispersal networks, but also underscores the relevance of evolutionary history in shaping species roles in ecological communities.     

Morphological integration and serial homology: A case study of the cranium and anterior vertebrae in salamanders

Serial homology or the repetition of equivalent developmental units and their derivatives is a phenomenon encountered in a variety of organisms, with the vertebrate axial skeleton as one of the most notable examples. Serially homologous structures can be viewed as an appropriate model system for studying morphological integration and modularity, due to the strong impact of development on their covariation. Here, we explored the pattern of morphological integration of the cranium and the first three serially homologous structures (atlas, first, and second trunk vertebrae) in salamandrid salamanders, using micro-CT scanning and three-dimensional geometric morphometrics. We explored the integration between structures at static and evolutionary levels. Effects of allometry on patterns of modularity were also taken into account. At the static level (within species), we analyzed inter-individual variation in shape to detect functional modules and intra-individual variation to detect developmental modules. Significant integration (based on inter-individual variation) among all structures was detected and allometry is shown to be an important integrating factor. The pattern of intra-individual, asymmetric variation indicates statistically significant developmental integration between the cranium and the atlas and between the first two trunk vertebrae. At the evolutionary level (among species), the cranium, atlas, and trunk vertebrae separate as different modules. Our results show that morphological integration at the evolutionary level coincides with morphological and functional differentiation of the axial skeleton, allowing the more or less independent evolutionary changes of the cranial skeleton and the vertebral column, regardless of the relatively strong integration at the static level. The observed patterns of morphological integration differ across levels, indicating different impacts of developmental and phylogenetic constraints and functional demands.     

Multilevel assessment of the Lacertid lizard cranial modularity

Different factors and processes that produce phenotypic variation at the individual, population, or interspecific level can influence or alter the covariance structure among morphological traits. Therefore, studies of the patterns of integration and modularity at multiple levels—static, ontogenetic, and evolutionary, can provide invaluable data on underlying factors and processes that structured morphological variation, directed, or constrained evolutionary changes. Our dataset, consisting of cranium shape data for 14 lizard species from the family Lacertidae, with substantial samples of hatchlings and adults along with their inferred evolutionary relationships, enabled us to assess modularity and morphological integration at all three levels. Five, not mutually exclusive modularity hypotheses of lizard cranium, were tested, and the effects of allometry on intensity and the pattern of integration and modularity were estimated. We used geometric morphometrics to extract symmetric and asymmetric, as well as allometric and nonallometric, components of shape variation. At the static level, firm confirmation of cranial modularity was found for hypotheses which separate anterior and posterior functional compartments of the skull. At the ontogenetic level, two alternative hypotheses (the “anteroposterior” and “neurodermatocranial” hypotheses) of ventral cranial modularity were confirmed. At the evolutionary level, the “neurodermatocranial” hypothesis was confirmed for the ventral cranium, which is in accordance with the pattern observed at the ontogenetic level. The observed pattern of static modularity could be driven by functional demands and can be regarded as adaptive. Ontogenetic modularity and evolutionary modularity show the same developmental origin, indicating conservatism of modularity patterns driven by developmental constraints.     

Mosaic heterochrony and evolutionary modularity: the trilobite genus Zacanthopsis as a case study

Logical connections exist between evolutionary modularity and heterochrony, two unifying and structuring themes in the expanding field of evolutionary developmental biology. The former sees complex phenotypes as being made up of semi-independent units of evolutionary transformation; the latter requires such a modular organization of phenotypes to occur in a localized or mosaic fashion. This conceptual relationship is illustrated here by analyzing the evolutionary changes in the cranidial ontogeny of two related species of Cambrian trilobites. With arguments from comparative developmental genetics and functional morphology, we delineate putative evolutionary modules within the cranidium and examine patterns of evolutionary changes in ontogeny at both global and local scales. Results support a case of mosaic heterochrony, that is, a combination of local heterochronies affecting the different parts individuated in the cranidium, leading to the complex pattern of allometric repatterning observed at the global scale. Through this example, we show that recasting morphological analyses of complex phenotypes with a priori knowledge or hypotheses about their organizational and variational properties can significantly improve our interpretation and understanding of evolutionary changes among related taxa, fossil and extant. Such considerations open avenues to investigate the large-scale dynamics of modularity and its role in phenotypic evolution.     

Functional decoupling between flowers and leaves in the Ameroglossum pernambucense complex can facilitate local adaptation across a pollinator and climatic heterogeneous landscape

Decoupling between floral and leaf traits is expected in plants with specialized pollination systems to assure a precise flower–pollinator fit, irrespective of leaf variation associated with environmental heterogeneity (functional modularity). Nonetheless, developmental interactions among floral traits also decouple flowers from leaves regardless of selection pressures (developmental modularity). We tested functional modularity in the hummingbird‐pollinated flowers of the Ameroglossum pernambucense complex while controlling for developmental modularity. Using two functional traits responsible for flower–pollinator fit [floral tube length (TL) and anther–nectary distance (AN)], one floral trait not linked to pollination [sepal length (SL), control for developmental modularity] and one leaf trait [leaf length (LL)], we found evidence of flower functional modularity. Covariation between TL and AN was ca. two‐fold higher than the covariation of either of these traits with sepal and leaf lengths, and variations in TL and AN, important for a precise flower–pollinator fit, were smaller than SL and LL variations. Furthermore, we show that previously reported among‐population variation of flowers associated with local pollinator phenotypes was independent from SL and LL variations. These results suggest that TL and AN are functionally linked to fit pollinators and sufficiently decoupled from developmentally related floral traits (SL) and vegetative traits (LL). These results support previous evidences of population differentiation due to local adaptation in the A. pernambucense complex and shed light on the role of flower–leaf decoupling for local adaptation in species distributed across biotic and abiotic heterogeneous landscapes.     

Association between integration structure and functional evolution in the opercular four‐bar apparatus of the threespine stickleback,Gasterosteus aculeatus (Pisces: Gasterosteidae)

Phenotypes may evolve to become integrated in response to functional demands. Once evolved, integrated phenotypes, often modular, can also influence the trajectory of subsequent responses to selection. Clearly, connecting modularity and functionally adaptive evolution has been challenging. The teleost skull and jaw structures are useful for understanding this connection because of the key roles that these structures play in feeding in novel environments with different prey resources. In the present study, we examined such a structure in the threespine stickleback: the opercular four‐bar lever that functions in jaw opening. Comparing oceanic and two fresh‐water populations, we find marked phenotypic divergence in the skull opercular region, and the major axes of morphological and functional variation of the lever are found to be highly correlated. All three populations share the same global skull integration structure, and a conserved, strongly‐supported modular organization is evident in the region encompassing the lever. Importantly, a boundary between two modules that subdivides the lever apparatus corresponds to the region of most prominent morphological evolution. The matched modular phenotypic and functional architecture of head and jaw structures of stickleback therefore may be important for facilitating their rapid adaptive transitions between highly divergent habitats. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 375–390.     

Ontogenetic scaling of the human nose in a longitudinal sample: Implications for genus Homo facial evolution

Researchers have hypothesized that nasal morphology, both in archaic Homo and in recent humans, is influenced by body mass and associated oxygen consumption demands required for tissue maintenance. Similarly, recent studies of the adult human nasal region have documented key differences in nasal form between males and females that are potentially linked to sexual dimorphism in body size, composition, and energetics. To better understand this potential developmental and functional dynamic, we first assessed sexual dimorphism in the nasal cavity in recent humans to determine when during ontogeny male‐female differences in nasal cavity size appear. Next, we assessed whether there are significant differences in nasal/body size scaling relationships in males and females during ontogeny. Using a mixed longitudinal sample we collected cephalometric and anthropometric measurements from n = 20 males and n = 18 females from 3.0 to 20.0+ years of age totaling n = 290 observations. We found that males and females exhibit similar nasal size values early in ontogeny and that sexual dimorphism in nasal size appears during adolescence. Moreover, when scaled to body size, males exhibit greater positive allometry in nasal size compared to females. This differs from patterns of sexual dimorphism in overall facial size, which are already present in our earliest age groups. Sexually dimorphic differences in nasal development and scaling mirror patterns of ontogenetic variation in variables associated with oxygen consumption and tissue maintenance. This underscores the importance of considering broader systemic factors in craniofacial development and may have important implications for the study of patters craniofacial evolution in the genus Homo. Am J Phys Anthropol 153:52–60, 2014. © 2013 Wiley Periodicals, Inc.     

Effects of environmental perturbations during postnatal development on the phenotypic integration of the skull

Integration and modularity are fundamental determinants of how natural selection effects evolutionary change in complex multivariate traits. Interest in the study of the specific developmental basis of integration through experimental approaches is fairly recent and it has mainly focused on its genetic determinants. In this study, we present evidence that postnatal environmental perturbations can modify the covariance structure by influencing the variance of some developmental processes relative to the variances of other processes that contribute to such structure. We analyzed the effects of the reduction of nutrient supply in different ontogenetic stages (i.e. before and after weaning, and from birth to adulthood) in Rattus norvegicus. Our results show that this environmental perturbation alters the phenotypic variation/covariation structure of the principal modules of the skull (base, vault, and face). The covariance matrices of different treatment groups exhibit low correlations and are significantly different, indicating that the treatments influence covariance structure. Postnatal nutrient restriction also increases the variance of somatic growth. This increased variance drives an increase in overall integration of cranial morphology through the correlated allometric effects of size variation. The extent of this increase in integration depends on the time and duration of the nutritional restriction. These results support the conclusion that environmental perturbations can influence integration and thus covariance structure via developmental plasticity.     

ADAPTIVE RADIATIONS,ECOLOGICAL SPECIALIZATION,AND THE EVOLUTIONARY INTEGRATION OF COMPLEX MORPHOLOGICAL STRUCTURES

The evolutionary integration of complex morphological structures is a macroevolutionary pattern in which morphogenetic components evolve in a coordinated fashion, which can result from the interplay among processes of developmental, genetic integration, and different types of selection. We tested hypotheses of ecological versus developmental factors underlying patterns of within‐species and evolutionary integration in the mandible of phyllostomid bats, during the most impressive ecological and morphological radiation among mammals. Shape variation of mandibular morphogenetic components was associated with diet, and the transition of integration patterns from developmental to within‐species to evolutionary was examined. Within‐species (as a proxy to genetic) integration in different lineages resembled developmental integration regardless of diet specialization, however, evolutionary integration patterns reflected selection in different mandibular components. For dietary specializations requiring extensive functional changes in mastication patterns or biting, such as frugivores and sanguivores, the evolutionary integration pattern was not associated with expected within‐species or developmental integration. On the other hand, specializations with lower mastication demands or without major functional reorganization (such as nectarivores and carnivores), presented evolutionary integration patterns similar to the expected developmental pattern. These results show that evolutionary integration patterns are largely a result of independent selection on specific components regardless of developmental modules.     

Functional features and ontogenetic changes in reproductive allocation and partitioning strategies of plant modules

In the concept of modularity, plant modules are considered as iterative units and their changes are analyzed in terms of number or size. This paper, however, analyses changes with respect to the reproductive functional performance of modules and individual plant age. Patterns of resource allocation and partitioning in reproductive modules (fruits) are compared between two different age groups of a bushy perennial plant, Cistus ladanifer.Although modules do not differ in their allocation strategies (young plant modules produce the same seed and packing/protective structure biomass as old plant modules), their partitioning strategies change with plant age: young plant modules produce a larger number of lighter seeds than old plants. These differences have a direct consequence on the plant pre-dispersal fitness, which is not counteracted by insect predation on reproductive modules. These results are empirical evidence of a differentiation in the performance of reproductive modules with the ontogenetic development of this plant species. We think that the consideration of such kind of changes in module features is essential in the analysis of the iterative construction of plants.