EvolQG - An R package for evolutionary quantitative genetics

Exact expressions and numerical evaluation of average evolvability measures for characterizing and comparing [Formula: see text] matrices

Theory predicts that the additive genetic covariance ([Formula: see text]) matrix determines a population's short-term (in)ability to respond to directional selection-evolvability in the Hansen-Houle sense-which is typically quantified and compared via certain scalar indices called evolvability measures. Often, interest is in obtaining the averages of these measures across all possible selection gradients, but explicit formulae for most of these average measures have not been known. Previous authors relied either on approximations by the delta method, whose accuracy is generally unknown, or Monte Carlo evaluations (including the random skewers analysis), which necessarily involve random fluctuations. This study presents new, exact expressions for the average conditional evolvability, average autonomy, average respondability, average flexibility, average response difference, and average response correlation, utilizing their mathematical structures as ratios of quadratic forms. The new expressions are infinite series involving top-order zonal and invariant polynomials of matrix arguments, and can be numerically evaluated as their partial sums with, for some measures, known error bounds. Whenever these partial sums numerically converge within reasonable computational time and memory, they will replace the previous approximate methods. In addition, new expressions are derived for the average measures under a general normal distribution for the selection gradient, extending the applicability of these measures into a substantially broader class of selection regimes.

Mammalian forelimb evolution is driven by uneven proximal-to-distal morphological diversity

Vertebrate limb morphology often reflects the environment due to variation in locomotor requirements. However, proximal and distal limb segments may evolve differently from one another, reflecting an anatomical gradient of functional specialization that has been suggested to be impacted by the timing of development. Here, we explore whether the temporal sequence of bone condensation predicts variation in the capacity of evolution to generate morphological diversity in proximal and distal forelimb segments across more than 600 species of mammals. Distal elements not only exhibit greater shape diversity, but also show stronger within-element integration and, on average, faster evolutionary responses than intermediate and upper limb segments. Results are consistent with the hypothesis that late developing distal bones display greater morphological variation than more proximal limb elements. However, the higher integration observed within the autopod deviates from such developmental predictions, suggesting that functional specialization plays an important role in driving within-element covariation. Proximal and distal limb segments also show different macroevolutionary patterns, albeit not showing a perfect proximo-distal gradient. The high disparity of the mammalian autopod, reported here, is consistent with the higher potential of development to generate variation in more distal limb structures, as well as functional specialization of the distal elements.

Secondary contact rather than coexistence-Erebia butterflies in the Alps

Secondary contact zones are ideal systems to study the processes that govern the evolution of reproductive barriers, especially at advanced stages of the speciation process. An increase in reproductive isolation resulting from selection against maladaptive hybrids is thought to contribute to reproductive barrier buildup in secondary contact zones. Although such processes have been invoked for many systems, it remains unclear to which extent they influence contact zone dynamics in nature. Here, we study a very narrow contact zone between the butterfly species Erebia cassioides and Erebia tyndarus in the Swiss Alps. We quantified phenotypic traits related to wing shape and reproduction as well as ecology to compare the degree of intra- and interspecific differentiation. Even though only very few first-generation hybrids occur, we find no strong indications for current reinforcing selection, suggesting that if reinforcement occurred in our system, it likely operated in the past. Additionally, we show that both species differ less in their ecological niche at the contact zone than elsewhere, which could explain why coexistence between these butterflies may currently not be possible.

Low genetic variation of foliar traits among Prosopis chilensis (Leguminosae) provenances

Prosopis chilensis (Molina) Stuntz (Leguminosae) is a valuable native species in Argentina that has been proposed to be used in reforestation, afforestation and restoration programmes. Natural provenances show important differentiation in height, shape, spine size, fruits and foliar traits throughout their distribution in the semiarid Monte ecoregion. The goal of this work was to characterize the genetic basis of the leaf variation in P. chilensis aiming to contribute to the improvement management program. We analyzed morphological variation and estimate narrow sense heritability for ten quantitative traits from a provenance-progeny trial founded from open pollinated families. We assessed the variance components by a generalized linear mixed model. Differences among provenances were quantified through univariate Q_ST statistics and multivariate discriminant analysis of principal components. Finally, univariate and multivariate neutrality test were conducted to unveil the evolutionary forces that shape the variation. Univariate and multivariate analysis showed low genetic variation in foliar traits among provenances grown in the common garden. Consistently, the Q_ST estimates for each trait were low. Both, the univariate (Q_ST-F_ST comparison) and the multivariate neutrality test suggest that the leaf variation among provenances may be shaped by genetic drift rather than selective forces. Heritability estimates were significant only for leaflet apex and leaflet apex/leaflet area. Since genetic variation for most foliar traits among provenances estimated under controlled environmental conditions were very low or absent, the variation described in the wild would be explained merely by plastic response to varying environments. These results are discussed in terms of adaptive strategies and the use of different provenances as seed sources within the framework of the improvement program. It is expected that P. chilensis seeds or seedlings from trees selected under economical criteria will be able to develop in different areas thanks to the phenotypic plasticity of leaf traits.

Development and function explain the modular evolution of phalanges in gecko lizards

Selective regimes favouring the evolution of functional specialization probably affect covariation among phenotypic traits. Phalanges of most tetrapods develop from a conserved module that constrains their relative proportions. In geckos, however, biomechanical specializations associated with adhesive toepads involve morphological variation in the autopodium and might reorganize such modular structures. We tested two hypotheses to explain the modular architecture of hand bones in geckos, one based on developmental interactions and another incorporating functional associations related to locomotion, and compared the empirical support for each hypothetical module between padded and padless lineages. We found strong evidence for developmental modules in most species, which probably reflects embryological constraints during phalangeal formation. Although padded geckos exhibit a functional specialization involving the hyperextension of the distal phalanges that is absent in padless species, the padless species are the ones that show a distal functional module with high integration. Some ancestrally padless geckos apparently deviate from developmental predictions and present a relatively weak developmental module of phalanges and a strongly integrated distal module, which may reflect selective regimes involving incipient frictional adhesion in digit morphology. Modularity of digit elements seems dynamic along the evolutionary history of geckos, being associated with the presence/absence of adhesive toepads.

Cranial shape diversification in horses: variation and covariation patterns under the impact of artificial selection

The potential of artificial selection to dramatically impact phenotypic diversity is well known. Large-scale morphological changes in domestic species, emerging over short timescales, offer an accelerated perspective on evolutionary processes. The domestic horse (Equus caballus) provides a striking example of rapid evolution, with major changes in morphology and size likely stemming from artificial selection. However, the microevolutionary mechanisms allowing to generate this variation in a short time interval remain little known. Here, we use 3D geometric morphometrics to quantify skull morphological diversity in the horse, and investigate modularity and integration patterns to understand how morphological associations contribute to cranial evolvability in this taxon. We find that changes in the magnitude of cranial integration contribute to the diversification of the skull morphology in horse breeds. Our results demonstrate that a conserved pattern of modularity does not constrain large-scale morphological variations in horses and that artificial selection has impacted mechanisms underlying phenotypic diversity to facilitate rapid shape changes. More broadly, this study demonstrates that studying microevolutionary processes in domestic species produces important insights into extant phenotypic diversity.

Patterns of skeletal integration in birds reveal that adaptation of element shapes enables coordinated evolution between anatomical modules

Birds show tremendous ecological disparity in spite of strong biomechanical constraints imposed by flight. Modular skeletal evolution is generally accepted to have facilitated this, with distinct body regions showing semi-independent evolutionary trajectories. However, this hypothesis has received little scrutiny. We analyse evolutionary modularity and ecomorphology using three-dimensional data from across the entire skeleton in a phylogenetically broad sample of extant birds. We find strongly modular evolution of skeletal element sizes within body regions (head, trunk, forelimb and hindlimb). However, element shapes show substantially less modularity, have stronger relationships to ecology, and provide evidence that ecological adaptation involves coordinated evolution of elements across different body regions. This complicates the straightforward paradigm in which modular evolution facilitated the ecological diversification of birds. Our findings suggest the potential for undetected patterns of morphological evolution in even well-studied groups, and advance the understanding of the interface between evolutionary integration and ecomorphology.

Intraspecific trait variability and community assembly in hawkmoths (Lepidoptera: Sphingidae) across an elevational gradient in the eastern Himalayas, India

We investigated some aspects of hawkmoth community assembly at 13 elevations along a 200- to 2770-m transect in the eastern Himalayas, a little studied biodiversity hot spot of global importance. We measured the morphological traits of body mass, wing loading, and wing aspect ratio of 3,301 free-ranging individuals of 76 species without having to collect or even constrain them. We used these trait measurements and T-statistic metrics to assess the strength of intracommunity ("internal") and extra-community ("external") filters which determine the composition of communities vis-a-vis the regional pool of species.The trait distribution of constituent species turned out to be nonrandom subsets of the community-trait distribution, providing strong evidence for internal filtering in all elevational communities. The external filter metric was more ambiguous. However, the elevational dependence of many metrics including that of the internal filter provided evidence for external (i.e., environmental) filtering. On average, a species occupied as much as 50%-75% of the total community-trait space, yet the T-statistic metric for internal filter was sufficiently sensitive to detect a strong nonrandom structure in the trait distribution.We suggest that the change in T-statistic metrics along the environmental gradient may provide more clues to the process of community assembly than previously envisaged. A large, smoothly varying and well-sampled environmental span would make it easier to discern them. Developing T-statistics for combined analysis of multiple traits will perhaps provide a more accurate picture of internal/filtering and niche complementarity. Moths are a hyperdiverse taxon and a very important component of many ecosystems. Our technique for accurately measuring body and wing dimensions of free-ranging moths can generate trait database for a large number of individuals in a time- and resource-efficient manner for a variety of community assembly studies using this important taxon.

Quercus species divergence is driven by natural selection on evolutionarily less integrated traits

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism's potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea (Matt.) Liebl. and Quercus robur L.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species' evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module "leaf size gradient" was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under "global regulation" in vegetative bud establishment and development. Lamina basal shape and intercalary veins' number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits' ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.

Variation in brown rat cranial shape shows directional selection over 120 years in New York City

Urbanization exposes species to novel environments and selection pressures that may change morphological traits within a population. We investigated how the shape and size of crania and mandibles changed over time within a population of brown rats (Rattus norvegicus) living in Manhattan, New York, USA, a highly urbanized environment. We measured 3D landmarks on the cranium and mandible of 62 adult individuals sampled in the 1890s and 2010s. Static allometry explained approximately 22% of shape variation in crania and mandible datasets, while time accounted for approximately 14% of variation. We did not observe significant changes in skull size through time or between the sexes. Estimating the P-matrix revealed that directional selection explained temporal change of the crania but not the mandible. Specifically, rats from the 2010s had longer noses and shorter upper molar tooth rows, traits identified as adaptive to colder environments and higher quality or softer diets, respectively. Our results highlight the continual evolution to selection pressures. We acknowledge that urban selection pressures impacting cranial shape likely began in Europe prior to the introduction of rats to Manhattan. Yet, our study period spanned changes in intensity of artificial lighting, human population density, and human diet, thereby altering various aspects of rat ecology and hence pressures on the skull.

Loss of consumers constrains phenotypic evolution in the resulting food web

The loss of biodiversity is altering the structure of ecological networks; however, we are currently in a poor position to predict how these altered communities will affect the evolution of remaining populations. Theory on fitness landscapes provides a framework for predicting how selection alters the evolutionary trajectory and adaptive potential of populations, but often treats the network of interacting populations as a “black box.” Here, we integrate ecological networks and fitness landscapes to examine how changes in food‐web structure shape phenotypic evolution. We conducted a field experiment that removed a guild of larval parasitoids that imposed direct and indirect selection pressures on an insect herbivore. We then measured herbivore survival as a function of three key phenotypic traits to estimate directional, quadratic, and correlational selection gradients in each treatment. We used these selection gradients to characterize the slope and curvature of the fitness landscape to understand the direct and indirect effects of consumer loss on phenotypic evolution. We found that the number of traits under directional selection increased with the removal of larval parasitoids, indicating evolution was more constrained toward a specific combination of traits. Similarly, we found that the removal of larval parasitoids altered the curvature of the fitness landscape in such a way that tended to decrease the evolvability of the traits we measured in the next generation. Our results suggest that the loss of trophic interactions can impose greater constraints on phenotypic evolution. This indicates that the simplification of ecological communities may constrain the adaptive potential of remaining populations to future environmental change.

Modelling of population structure through contemporary groups in genetic evaluation

Background

Forest trees can occupy extensive geography and environmentally highly variable areas which result in high genetic variability in the direction of pressure from natural selection. At the same time, the majority of conifer species are wind-pollinated from both short and long distances, resulting in wide-spread gene flow, which can lead to maladaptation to local conditions. Quantitative analyses of provenance/progeny tests correct for genetic differences between populations to ensure unbiased genetic parameters are obtained. Commonly, the provenance effect is fitted as a fixed term or can be implemented as a contemporary group in the pedigree.

Results

The use of a provenance effect, either as a fixed term or as the same contemporary groups in both maternal and paternal sides of the pedigree, resulted in fairly similar precision of genetic parameters in our case. However, when we developed a phantom contemporary group for the paternal side of the pedigree that considered a different genetic quality of pollen compared with the maternal contribution from trees in the local environment, the model fit and accuracy of breeding values increased.

Conclusion

Consideration of the mating dynamics and the vector of gene flow are important factors in modelling contemporary genetic groups, particularly when implementing pedigrees within a mixed model framework to obtain unbiased estimates of genetic parameters. This approach is especially important in traits involved in local adaptation.

Mouse Skull Mean Shape and Shape Robustness Rely on Different Genetic Architectures and Different Loci

The genetic architecture of skull shape has been extensively studied in mice and the results suggest a highly polygenic and additive basis. In contrast few studies have explored the genetic basis of the skull variability. Canalization and developmental stability are the two components of phenotypic robustness. They have been proposed to be emergent properties of the genetic networks underlying the development of the trait itself, but this hypothesis has been rarely tested empirically. Here we use outbred mice to investigate the genetic architecture of canalization of the skull shape by implementing a genome-wide marginal epistatic test on 3D geometric morphometric data. The same data set had been used previously to explore the genetic architecture of the skull mean shape and its developmental stability. Here, we address two questions: (1) Are changes in mean shape and changes in shape variance associated with the same genomic regions? and (2) Do canalization and developmental stability rely on the same loci and genetic architecture and do they involve the same patterns of shape variation? We found that unlike skull mean shape, among-individual shape variance and fluctuating asymmetry (FA) show a total lack of additive effects. They are both associated with complex networks of epistatic interactions involving many genes (protein-coding and regulatory elements). Remarkably, none of the genomic loci affecting mean shape contribute these networks despite their enrichment for genes involved in craniofacial variation and diseases. We also found that the patterns of shape FA and individual variation are largely similar and rely on similar multilocus epistatic genetic networks, suggesting that the processes channeling variation within and among individuals are largely common. However, the loci involved in these two networks are completely different. This in turn underlines the difference in the origin of the variation at these two levels, and points at buffering processes that may be specific to each level.

The modular organization of roe deer (Capreolus capreolus) body during ontogeny: the effects of sex and habitat

Background

As a small artiodactyl, the roe deer (Capreolus capreolus L.) is characterized by biological plasticity and great adaptability demonstrated by their survival under a wide variety of environmental conditions. In order to depict patterns of phenotypic variation of roe deer body this study aims to quantify variation during ontogenetic development and determine how sex-specific reproductive investment and non-uniform habitat differences relate to phenotypic variation and do these differential investments mold the patterns of phenotypic variation through modular organisation.

Results

Patterns of phenotypic correlation among body traits change during the ontogeny of roe deer, with differential influence of sex and habitat type. Modularity was found to be a feature of closed habitats with trunk+forelimbs+hindlimbs as the best supported integration/modularity hypothesis for both sexes. The indices of integration and evolvability vary with habitat type, age and sex where increased integration is followed by decreased evolvability.

Conclusion

This is the first study that quantifies patterns of correlation in the roe deer body and finds pronounced changes in correlation structure during ontogeny affected by sex and habitat type. The correlation structure of the roe deer body is developmentally written over the course of ontogeny but we do not exclude the influence of function on ontogenetic changes. Modularity arises with the onset of reproduction (subadults not being modular) and is differentially expressed in males and females from different habitats. Both adult males and females show modularity in primordial, closed habitats. Overall, all these findings are important as they provide support to the idea that modularity can evolve at the population level and change fast within a species.

Adaptive radiation along a deeply conserved genetic line of least resistance in Anolis lizards

On microevolutionary timescales, adaptive evolution depends upon both natural selection and the underlying genetic architecture of traits under selection, which may constrain evolutionary outcomes. Whether such genetic constraints shape phenotypic diversity over macroevolutionary timescales is more controversial, however. One key prediction is that genetic constraints should bias the early stages of species divergence along “genetic lines of least resistance” defined by the genetic (co)variance matrix, G. This bias is expected to erode over time as species means and G matrices diverge, allowing phenotypes to evolve away from the major axis of variation. We tested for evidence of this signal in West Indian Anolis lizards, an iconic example of adaptive radiation. We found that the major axis of morphological evolution was well aligned with a major axis of genetic variance shared by all species despite separation times of 20–40 million years, suggesting that divergence occurred along a conserved genetic line of least resistance. Further, this signal persisted even as G itself evolved, apparently because the largest evolutionary changes in G were themselves aligned with the line of genetic least resistance. Our results demonstrate that the signature of genetic constraint may persist over much longer timescales than previously appreciated, even in the presence of evolving genetic architecture. This pattern may have arisen either because pervasive constraints have biased the course of adaptive evolution or because the G matrix itself has been shaped by selection to conform to the adaptive landscape.

Testing the stability of behavioural coping style across stress contexts in the Trinidadian guppy

Within populations, individuals can vary in stress response, a multivariate phenomenon comprising neuroendocrine, physiological and behavioural traits.Verbal models of individual stress "coping style" have proposed that the behavioural component of this variation can be described as a single axis, with each individual's coping style being consistent across time and stress contexts.Focusing on this behavioural component of stress response and combining repeated measures of multiple traits with a novel multivariate modelling framework, we test for the existence of coping style variation and assess its stability across contexts in the Trinidadian guppy (Poecilia reticulata).Specifically, we test the following hypotheses: (1) there exists repeatable among-individual behavioural (co)variation ("personality") within a mild stress context consistent with a risk-averse-risk-prone continuum of behavioural coping style, (2) there is population-level plasticity in behaviour as a function of stressor severity, (3) there is among-individual variation in plasticity (i.e. IxE), and (4) the presence of IxE reduces cross-context stability of behavioural coping style.We found significant repeatable among-individual behavioural (co)variation in the mild stress context (open field trial), represented as an I matrix. However, I was not readily described by a simple risk-averse-risk-prone continuum as posited by the original coping style model. We also found strong evidence for population-level changes in mean behaviour with increasing stressor severity (simulated avian and piscine predation risks).Single-trait analyses did show the presence of individual-by-environment interactions (IxE), as among-individual cross-context correlations were significantly less than +1. However, multitrait analysis revealed the consequences of this plasticity variation were minimal. Specifically, we found little evidence for changes in the structure of I between mild and moderate stress contexts overall, and only minor changes between the two moderate contexts (avian vs. piscine predator).We show that a multivariate approach to assessing changes in among-individual (co)variance across contexts can prevent the over-interpretation of statistically significant, but small, individual-by-environment effects. While behavioural flexibility enables populations (and individuals) to respond rapidly to changes in the environment, multivariate personality structure can be conserved strongly across such contexts. A plain language summary is available for this article.

Shape Covariation (or the Lack Thereof) Between Vertebrae and Other Skeletal Traits in Felids: The Whole is Not Always Greater than the Sum of Parts

Within carnivorans, cats show comparatively little disparity in overall morphology, with species differing mainly in body size. However, detailed shape analyses of individual osteological structures, such as limbs or skulls, have shown that felids display significant morphological differences that correlate with their observed ecological and behavioural ranges. Recently, these shape analyses have been extended to the felid axial skeleton. Results demonstrate a functionally-partitioned vertebral column, with regions varying greatly in level of correlation between shape and ecology. Moreover, a clear distinction is evident between a phylogenetically-constrained neck region and a selection-responsive posterior spine. Here, we test whether this regionalisation of function reflected in vertebral column shape is also translated into varying levels of phenotypic integration between this structure and most other skeletal elements. We accomplish this comparison by performing pairwise tests of integration between vertebral and other osteological units, quantified with 3D geometric morphometric data and analysed both with and without phylogenetic correction. To our knowledge, this is the first study to test for integration across a comprehensive sample of whole-skeleton elements. Our results show that, prior to corrections, strong covariation is present between vertebrae across the vertebral column and all other elements, with the exception of the femur. However, most of these significant correlations disappear after correcting for phylogeny, which is a significant influence on cranial and limb morphology of felids and other carnivorans. Our results thus suggest that the vertebral column of cats displays relative independence from other skeletal elements and may represent several distinct evolutionary morphological modules.

Evolution of a complex phenotype with biphasic ontogeny: Contribution of development versus function and climatic variation to skull modularity in toads

The theory of morphological integration and modularity predicts that if functional correlations among traits are relevant to mean population fitness, the genetic basis of development will be molded by stabilizing selection to match functional patterns. Yet, how much functional interactions actually shape the fitness landscape is still an open question. We used the anuran skull as a model of a complex phenotype for which we can separate developmental and functional modularity. We hypothesized that functional modularity associated to functional demands of the adult skull would overcome developmental modularity associated to bone origin at the larval phase because metamorphosis would erase the developmental signal. We tested this hypothesis in toad species of the Rhinella granulosa complex using species phenotypic correlation pattern (P-matrices). Given that the toad species are distributed in very distinct habitats and the skull has important functions related to climatic conditions, we also hypothesized that differences in skull trait covariance pattern are associated to differences in climatic variables among species. Functional and hormonal-regulated modules are more conspicuous than developmental modules only when size variation is retained on species P-matrices. Without size variation, there is a clear modularity signal of developmental units, but most species have the functional model as the best supported by empirical data without allometric size variation. Closely related toad species have more similar climatic niches and P-matrices than distantly related species, suggesting phylogenetic niche conservatism. We infer that the modularity signal due to embryonic origin of bones, which happens early in ontogeny, is blurred by the process of growth that occurs later in ontogeny. We suggest that the species differing in the preferred modularity model have different demands on the orbital functional unit and that species contrasting in climate are subjected to divergent patterns of natural selection associated to neurocranial allometry and T3 hormone regulation.

Intense natural selection preceded the invasion of new adaptive zones during the radiation of New World leaf-nosed bats

The family Phyllostomidae, which evolved in the New World during the last 30 million years, represents one of the largest and most morphologically diverse mammal families. Due to its uniquely diverse functional morphology, the phyllostomid skull is presumed to have evolved under strong directional selection; however, quantitative estimation of the strength of selection in this extraordinary lineage has not been reported. Here, we used comparative quantitative genetics approaches to elucidate the processes that drove cranial evolution in phyllostomids. We also quantified the strength of selection and explored its association with dietary transitions and specialization along the phyllostomid phylogeny. Our results suggest that natural selection was the evolutionary process responsible for cranial diversification in phyllostomid bats. Remarkably, the strongest selection in the phyllostomid phylogeny was associated with dietary specialization and the origination of novel feeding habits, suggesting that the adaptive diversification of phyllostomid bats was triggered by ecological opportunities. These findings are consistent with Simpson’s quantum evolutionary model of transitions between adaptive zones. The multivariate analyses used in this study provides a powerful tool for understanding the role of evolutionary processes in shaping phenotypic diversity in any group on both micro- and macroevolutionary scales.

The evolution of phenotypic integration: How directional selection reshapes covariation in mice

Variation is the basis for evolution, and understanding how variation can evolve is a central question in biology. In complex phenotypes, covariation plays an even more important role, as genetic associations between traits can bias and alter evolutionary change. Covariation can be shaped by complex interactions between loci, and this genetic architecture can also change during evolution. In this article, we analyzed mouse lines experimentally selected for changes in size to address the question of how multivariate covariation changes under directional selection, as well as to identify the consequences of these changes to evolution. Selected lines showed a clear restructuring of covariation in their cranium and, instead of depleting their size variation, these lines increased their magnitude of integration and the proportion of variation associated with the direction of selection. This result is compatible with recent theoretical works on the evolution of covariation that take the complexities of genetic architecture into account. This result also contradicts the traditional view of the effects of selection on available covariation and suggests a much more complex view of how populations respond to selection.

Morphological modularity in the vertebral column of Felidae (Mammalia, Carnivora)

BACKGROUND

Previous studies have demonstrated that the clear morphological differences among vertebrae across the presacral column are accompanied by heterogeneous functional signals in vertebral shape. Further, several lines of evidence suggest that the mammalian axial skeleton is a highly modular structure. These include its composition of serial units, a trade-off between high shape variance and strong conservation of vertebral count, and direct association of regions with anterior expression sites of Hox genes. Here we investigate the modular organisation of the presacral vertebral column of modern cats (Felidae, Carnivora, Mammalia) with pairwise comparisons of vertebral shape covariation (i.e. integration) and evaluate our results against hypotheses of developmental and functional modularity. We used three-dimensional geometric morphometrics to quantify vertebral shape and then assessed integration between pairs of vertebrae with phylogenetic two-block partial least square analysis (PLS).

RESULTS

Six modules were identified in the pairwise analyses (vertebrae included are designated as 'C' for cervical, 'T' for thoracic, and 'L' for lumbar): an anterior module (C1 to T1); a transitional module situated between the last cervicals and first thoracics (C6 to T2); an anterior to middle thoracic set (T4 to T8); an anticlinal module (T10 and T11); a posterior set composed of the last two thoracics and lumbars (T12 to L7); and a module showing covariation between the cervicals and the posterior set (T12 to L7). These modules reflect shared developmental pathways, ossification timing, and observed ecological shape diversification in living species of felids.

CONCLUSIONS

We show here that patterns of shape integration reflect modular organisation of the vertebral column of felids. Whereas this pattern corresponds with hypotheses of developmental and functional regionalisation in the axial skeleton, it does not simply reflect major vertebral regions. This modularity may also have permitted vertebral partitions, specifically in the posterior vertebral column, to be more responsive to selection and achieve higher morphological disparity than other vertebral regions.

Fluctuations in Evolutionary Integration Allow for Big Brains and Disparate Faces

In theory, evolutionary modularity allows anatomical structures to respond differently to selective regimes, thus promoting morphological diversification. These differences can then influence the rate and direction of phenotypic evolution among structures. Here we use geometric morphometrics and phenotypic matrix statistics to compare rates of craniofacial evolution and estimate evolvability in the face and braincase modules of a clade of teleost fishes (Gymnotiformes) and a clade of mammals (Carnivora), both of which exhibit substantial craniofacial diversity. We find that the face and braincase regions of both clades display different degrees of integration. We find that the face and braincase evolve at similar rates in Gymnotiformes and the reverse in Carnivora with the braincase evolving twice as fast as the face. Estimates of evolvability and constraints in these modules suggest differential responses to selection arising from fluctuations in phylogenetic integration, thus influencing differential rates of skull-shape evolution in these two clades.

Concordance between stabilizing sexual selection, intraspecific variation, and interspecific divergence in Phymata

Empirical studies show that lineages typically exhibit long periods of evolutionary stasis and that relative levels of within-species trait covariance often correlate with the extent of between-species trait divergence. These observations have been interpreted by some as evidence of genetic constraints persisting for long periods of time. However, an alternative explanation is that both intra- and interspecific variation are shaped by the features of the adaptive landscape (e.g., stabilizing selection). Employing a genus of insects that are diverse with respect to a suite of secondary sex traits, we related data describing nonlinear phenotypic (sexual) selection to intraspecific trait covariances and macroevolutionary divergence. We found support for two key predictions (1) that intraspecific trait covariation would be aligned with stabilizing selection and (2) that there would be restricted macroevolutionary divergence in the direction of stabilizing selection. The observed alignment of all three matrices offers a point of caution in interpreting standing variability as metrics of evolutionary constraint. Our results also illustrate the power of sexual selection for determining variation observed at both short and long timescales and account for the apparently slow evolution of some secondary sex characters in this lineage.

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.

Are more diverse parts of the mammalian skull more labile?

Morphological variation is unevenly distributed within the mammalian skull; some of its parts have diversified more than others. It is commonly thought that this pattern of variation is mainly the result of the structural organization of the skull, as defined by the pattern and magnitude of trait covariation. Patterns of trait covariation can facilitate morphological diversification if they are aligned in the direction of selection, or these patterns can constrain diversification if oriented in a different direction. Within this theoretical framework, it is thought that more variable parts possess patterns of trait covariation that made them more capable of evolutionary change, that is, are more labile. However, differences in the degree of morphological variation among skull traits could arise despite variation in trait lability if, for example, some traits have evolved at a different rate and/or undergone stabilizing selection. Here, we test these hypotheses in the mammalian skull using 2D geometric morphometrics to quantify skull shape and estimating constraint, rates of evolution, and lability. Contrary to the expectations, more variable parts of the skull across mammalian species are less capable of evolutionary change than are less variable skull parts. Our results suggest that patterns of morphological variation in the skull could result from differences in rate of evolution and stabilizing selection.