Phenotypic variation and covariation indicate high evolvability of acoustic communication in crickets

Studying the genetic architecture of sexual traits provides insight into the rate and direction at which traits can respond to selection. Traits associated with few loci and limited genetic and phenotypic constraints tend to evolve at high rates typically observed for secondary sexual characters. Here, we examined the genetic architecture of song traits and female song preferences in the field crickets Gryllus rubens and Gryllus texensis. Song and preference data were collected from both species and interspecific F1 and F2 hybrids. We first analysed phenotypic variation to examine interspecific differentiation and trait distributions in parental and hybrid generations. Then, the relative contribution of additive and additive‐dominance variation was estimated. Finally, phenotypic variance–covariance ( P ) matrices were estimated to evaluate the multivariate phenotype available for selection. Song traits and preferences had unimodal trait distributions, and hybrid offspring were intermediate with respect to the parents. We uncovered additive and dominance variation in song traits and preferences. For two song traits, we found evidence for X‐linked inheritance. On the one hand, the observed genetic architecture does not suggest rapid divergence, although sex linkage may have allowed for somewhat higher evolutionary rates. On the other hand, P matrices revealed that multivariate variation in song traits aligned with major dimensions in song preferences, suggesting a strong selection response. We also found strong covariance between the main traits that are sexually selected and traits that are not directly selected by females, providing an explanation for the striking multivariate divergence in male calling songs despite limited divergence in female preferences.     

Morphometrics and the role of the phenotype in studies of the evolution of developmental mechanisms

Developmental mechanisms are usually assumed to evolve by natural selection of the morphological traits they produce. Therefore, information on phenotypic traits is an important component of comparative studies of development. Morphometrics permits the rigorous quantitative analysis of variation in organismal size and shape, and is increasingly being used in developmental contexts. The new methods of morphometrics combine a geometric concept of shape with the procedures of multivariate statistics, and constitute a powerful and flexible set of tools for analyzing morphological variation. This paper briefly reviews these methods and provides examples of their application in studies of genetic variation and developmental modularity. The results of morphometric analyses can be readily interpreted in relation to the geometry and anatomical structure of the parts under study. Genetic studies of shape in the mouse mandible found two recurrent patterns in environmental and genetic variation from different origins, suggesting that the development system 'channels' the phenotypic expression of variation in similar ways. Moreover, by analyzing the correlations of left-right asymmetries of morphometric traits, it is possible to delimit the spatial extent of developmental modules. These methods complement the experimental approaches of developmental biology and genetics, and can be expected to be especially fruitful in combination with them.     

Genetics of Mandible Form in the Mouse

The underlying determination of phenotypic variability and covariability is described for 14 traits that define the morphological size and shape of the mature mouse mandible. Variability is partitioned into components due to direct additive and dominance genetic effects, indirect maternal additive genetic effects, genetic covariance between direct additive and indirect maternal additive effects and common and residual environmental effects. Multivariate analyses of the dimensionality of genetic variability indicate several complex and independent genetic components underlie the morphological form of the mandible. The multidimensional nature of the genetic components suggests a complex picture with regard to the consequences of selection on mandibular form.     

THE ONTOGENETIC TRAJECTORY OF THE PHENOTYPIC COVARIANCE MATRIX, WITH EXAMPLES FROM CRANIOFACIAL SHAPE IN RATS AND HUMANS

Many classic quantitative genetic theories assume the covariance structure among adult phenotypic traits to be relatively static during evolution. But the cross-sectional covariance matrix arises from the joint variation of a large range of developmental processes and hence is not constant over the period during which a population of developing organisms is actually exposed to selection. To examine how development shapes the phenotypic covariance structure, we ordinate the age-specific covariance matrices of shape coordinates for craniofacial growth in rats and humans. The metric that we use for this purpose is given by the square root of the summed squared log relative eigenvalues. This is the natural metric on the space of positive-definite symmetric matrices, which we introduce and justify in a biometric context. In both species, the covariance matrices appear to change continually throughout the full period of postnatal development. The resulting ontogenetic trajectories alter their direction at major changes of the developmental programs whereas they are fairly straight in between. Consequently, phenotypic covariance matrices—and thus also response to selection—should be expected to vary both over ontogenetic and phylogenetic time scales as different phenotypes are necessarily produced by different developmental pathways.     

Integrating genetic and environmental forces that shape the evolution of geographic variation in a marine snail

Temporal and spatial patterns of phenotypic variation have traditionally been thought to reflect genetic differentiation produced by natural selection. Recently, however, there has been growing interest in how natural selection may shape the genetics of phenotypic plasticity to produce patterns of geographic variation and phenotypic evolution. Because the covariance between genetic and environmental influences can modulate the expression of phenotypic variation, a complete understanding of geographic variation requires determining whether these influences covary in the same (cogradient variation) or in opposing (countergradient variation) directions. We focus on marine snails from rocky intertidal shores as an ideal system to explore how genetic and plastic influences contribute to geographic and historical patterns of phenotypic variation. Phenotypic plasticity in response to predator cues, wave action, and water temperature appear to exert a strong influence on small and large-scale morphological variation in marine snails. In particular, plasticity in snail shell thickness: (i) may contribute to phenotypic evolution, (ii) appears to have evolved across small and large spatial scales, and (iii) may be driven by life history trade-offs tied to architectural constraints imposed by the shell. The plasticity exhibited by these snails represents an important adaptive strategy to the pronounced heterogeneity of the intertidal zone and undoubtedly has played a key role in their evolution.     

Genetic covariance structure of growth in the salamander Ambystoma macrodactylum

The size of an organism at any point during ontogeny often has fitness consequences through either direct selection on size or through selection on size-related morphological, performance, or life history traits. However, the evolutionary response to selection on size across ontogeny (a growth trajectory) may be limited by genetic correlations across ages. Here we characterize the phenotypic and genetic covariance structure of length and mass growth trajectories in a natural population of larval Ambystoma macrodactylum using function-valued quantitative genetic analyses and principal component decomposition. Most of the phenotypic and genetic variation in both growth trajectories appears to be confined to a single principal component describing a pattern of positive covariation among sizes across all ages. Higher order principal components with no significant associated genetic variation were identified for both trajectories, suggesting that evolution towards certain patterns of negative covariation between sizes across ages is constrained. The well-characterized positive relationship between size at metamorphosis and fitness in pond-breeding amphibians predicts that the across-age covariance structure will strongly limit evolution only if there is negative selection on size prior to metamorphosis. The pattern of genetic covariation observed in this study is similar to that observed in other vertebrate taxa, indicating that size may often be highly genetically and phenotypically integrated across ontogeny. Additionally, we find that phenotypic and genetic analyses of growth trajectories can yield qualitatively similar patterns of covariance structure.     

Contemporary and historical evolutionary processes interact to shape patterns of within-lake phenotypic divergences in polyphenic pumpkinseed sunfish, Lepomis gibbosus

Historical and contemporary evolutionary processes can both contribute to patterns of phenotypic variation among populations of a species. Recent studies are revealing how interactions between historical and contemporary processes better explain observed patterns of phenotypic divergence than either process alone. Here, we investigate the roles of evolutionary history and adaptation to current environmental conditions in structuring phenotypic variation among polyphenic populations of sunfish inhabiting 12 postglacial lakes in eastern North America. The pumpkinseed sunfish polyphenism includes sympatric ecomorphs specialized for littoral or pelagic lake habitats. First, we use population genetic methods to test the evolutionary independence of within-lake phenotypic divergences of ecomorphs and to describe patterns of genetic structure among lake populations that clustered into three geographical groupings. We then used multivariate analysis of covariance (MANCOVA) to partition body shape variation (quantified with geometric morphometrics) among the effects of evolutionary history (reflecting phenotypic variation among genetic clusters), the shared phenotypic response of all populations to alternate habitats within lakes (reflecting adaptation to contemporary conditions), and unique phenotypic responses to habitats within lakes nested within genetic clusters. All effects had a significant influence on body form, but the effects of history and the interaction between history and contemporary habitat were larger than contemporary processes in structuring phenotypic variation. This highlights how divergence can be better understood against a known backdrop of evolutionary history.     

Genetic Analysis of Size-Scaling Patterns in the Mouse Mandible

The relationship between multidimensional form of the adult mouse mandible and body size is examined from an ontogenetic perspective. The origin and ontogeny of phenotypic correlations are described in terms of genetic and environmental covariance patterns between adult skeletal morphology and growth in body weight. Different ontogenetic patterns are observed in the genetic correlations, and these can be related to the developmental as well as the functional aspects of mandibular form. The quantitative genetic aspects of craniomandibular growth and morphogenesis are explored, together with an examination of the impact of ontogenetic changes in the genetic variance-covariance structure on morphogenetic integration and evolution by selection.     

QUANTITATIVE GENETICS OF SHAPE IN CRICKET WINGS: DEVELOPMENTAL INTEGRATION IN A FUNCTIONAL STRUCTURE

The role of developmental and genetic integration for evolution is contentious. One hypothesis states that integration acts as a constraint on evolution, whereas an alternative is that developmental and genetic systems evolve to match the functional modularity of organisms. This study examined a morphological structure, the cricket wing, where developmental and functional modules are discordant, making it possible to distinguish the two alternatives. Wing shape was characterized with geometric morphometrics, quantitative genetic information was extracted using a full‐sibling breeding design, and patterns of developmental integration were inferred from fluctuating asymmetry of wing shape. The patterns of genetic, phenotypic, and developmental integration were clearly similar, but not identical. Heritabilities for different shape variables varied widely, but no shape variables were devoid of genetic variation. Simulated selection for specific shape changes produced predicted responses with marked deflections due to the genetic covariance structure. Three hypotheses of modularity according to the wing structures involved in sound production were inconsistent with the genetic, phenotypic, or developmental covariance structure. Instead, there appears to be strong integration throughout the wing. The hypothesis that genetic and developmental integration evolve to match functional modularity can therefore be rejected for this example.     

The role of environment and core‐margin effects on range‐wide phenotypic variation in a montane grasshopper

The integration of genetic information with ecological and phenotypic data constitutes an effective approach to gain insight into the mechanisms determining interpopulation variability and the evolutionary processes underlying local adaptation and incipient speciation. Here, we use the Pyrenean Morales grasshopper (Chorthippus saulcyi moralesi) as study system to (i) analyse the relative role of genetic drift and selection in range‐wide patterns of phenotypic differentiation and (ii) identify the potential selective agents (environment, elevation) responsible for variation. We also test the hypothesis that (iii) the development of dispersal‐related traits is associated with different parameters related to population persistence/turnover, including habitat suitability stability over the last 120 000 years, distance to the species distribution core and population genetic variability. Our results indicate that selection shaped phenotypic differentiation across all the studied morphological traits (body size, forewing length and shape). Subsequent analyses revealed that among‐population differentiation in forewing length was significantly explained by a temperature gradient, suggesting an adaptive response to thermoregulation or flight performance under contrasting temperature regimes. We found support for our hypothesis predicting a positive association between the distance to the species distribution core and the development of dispersal‐related morphology, which suggests an increased dispersal capability in populations located at range edges that, in turn, exhibit lower levels of genetic variability. Overall, our results indicate that range‐wide patterns of phenotypic variation are partially explained by adaptation in response to local environmental conditions and differences in habitat persistence between core and peripheral populations.     

Patterns of phenotypic covariation and correlation in modern humans as viewed from morphological integration

Proportionality of phenotypic and genetic distance is of crucial importance to adequately focus on population history and structure, and it depends on the proportionality of genetic and phenotypic covariance. Constancy of phenotypic covariances is unlikely without constancy of genetic covariation if the latter is a substantial component of the former. If phenotypic patterns are found to be relatively stable, the most probable explanation is that genetic covariance matrices are also stable. Factors like morphological integration account for such stability. Morphological integration can be studied by analyzing the relationships among morphological traits. We present here a comparison of phenotypic correlation and covariance structure among worldwide human populations. Correlation and covariance matrices between 47 cranial traits were obtained for 28 populations, and compared with design matrices representing functional and developmental constraints. Among-population differences in patterns of correlation and covariation were tested for association with matrices of genetic distances (obtained after an examination of 10 Alu-insertions) and with Mahalanobis distances (computed after craniometrical traits). All matrix correlations were estimated by means of Mantel tests. Results indicate that correlation and covariance structure in our species is stable, and that among-group correlation/covariance similarity is not related to genetic or phenotypic distance. Conversely, genetic and morphological distance matrices were highly correlated. Correlation and covariation patterns were largely associated with functional and developmental factors, which probably account for the stability of covariance patterns.     

Evolutionary integration and morphological diversification in complex morphological structures: mandible shape divergence in spiny rats (Rodentia, Echimyidae)

The rodent mandible has become a paradigm for studies on the development and evolution of complex morphological structures. We use a combination of geometric and multivariate morphometric methods in order to assess the correspondence between integration patterns and a priori biological models in the context of evolutionary shape divergence in the mandible of rodents of the family Echimyidae. The correlation of shape distances among operational taxonomic units (individuals, species, genera) in separate morphogenetic components allowed the construction of integration matrices among mandible components for data sets corresponding to varying levels of genetic divergence (intergeneric, interspecific, and intrapopulational). The integration matrices were associated with a priori biological (developmental, genetical, modular) models, and the maximum integration axes (singular warps) were compared with realized axes of maximum interspecific variation (relative warps). The integration pattern and intensity were not stable in data sets with different levels of genetic divergence, and the varying functional demands during the ecological radiation in the family were probably responsible for the differences in observed integration patterns. Developmental and genetic models were significantly associated with the interspecific integration patterns observed, suggesting a role for neutral evolution during the evolutionary divergence of mandible shape. However, directional and stabilizing selection were not discarded as processes responsible for the generation of interspecific integration. The choreography of the morphogenetic components in the mandible is highly flexible and the integrated groups of components can be reorganized depending on functional demands during evolutionary shape changes.     

Insect mating signal and mate preference phenotypes covary among host plant genotypes

Sexual selection acting on small initial differences in mating signals and mate preferences can enhance signal–preference codivergence and reproductive isolation during speciation. However, the origin of initial differences in sexual traits remains unclear. We asked whether biotic environments, a source of variation in sexual traits, may provide a general solution to this problem. Specifically, we asked whether genetic variation in biotic environments provided by host plants can result in signal–preference phenotypic covariance in a host‐specific, plant‐feeding insect. We used a member of the Enchenopa binotata species complex of treehoppers (Hemiptera: Membracidae) to assess patterns of variation in male mating signals and female mate preferences induced by genetic variation in host plants. We employed a novel implementation of a quantitative genetics method, rearing field‐collected treehoppers on a sample of naturally occurring replicated host plant clone lines. We found remarkably high signal–preference covariance among host plant genotypes. Thus, genetic variation in biotic environments influences the sexual phenotypes of organisms living on those environments in a way that promotes assortative mating among environments. This consequence arises from conditions likely to be common in nature (phenotypic plasticity and variation in biotic environments). It therefore offers a general answer to how divergent sexual selection may begin.     

Evolutionary transition between bee pollination and hummingbird pollination in Salvia: Comparing means,variances and covariances of corolla traits

Covariation among traits can modify the evolutionary trajectory of complex structures. This process is thought to operate at a microevolutionary scale, but its long‐term effects remain controversial because trait covariation can itself evolve. Flower morphology, and particularly floral trait (co)variation, has been envisioned as the product of pollinator‐mediated selection. Available evidence suggests that major changes in pollinator assemblages may affect the joint expression of floral traits and their phenotypic integration. We expect species within a monophyletic lineage sharing the same pollinator type will show not only similarity in trait means but also similar phenotypic variance‐covariance structures. Here, we tested this expectation using eighteen Salvia species pollinated either by bees or by hummingbirds. Our findings indicated a nonsignificant multivariate phylogenetic signal and a decoupling between means and variance‐covariance phenotypic matrices of floral traits during the evolution to hummingbird pollination. Mean trait value analyses revealed significant differences between bee‐ and hummingbird‐pollinated Salvia species although fewer differences were detected in the covariance structure between groups. Variance‐covariance matrices were much more similar among bee‐ than hummingbird‐pollinated species. This pattern is consistent with the expectation that, unlike hummingbirds, bees physically manipulate the flower, presumably exerting stronger selection pressures favouring morphological convergence among species. Overall, we conclude that the evolution of hummingbird pollination proceeded through different independent transitions. Thus, although the evolution of hummingbird pollination led to a new phenotypic optimum, the process involved the diversification of the covariance structure.     

Selection in a cycling population: differential response among skeletal traits

Abstract Population density cycles influence phenotypic evolution through both density‐dependent selection during periods of high density and through enhanced genetic drift during periods of low density. We investigated the response of different phenotypic traits to the same density cycles in a population of the yellow‐necked mouse, Apodemus flavicollis, from Bia?owieza National Park in Poland. We examined nonmetric skull traits, skull and mandible size, skull and mandible shape, and transferrin allele frequencies. We found that all of the traits changed significantly over the seven‐year study period. The greatest changes in nonmetric traits and mandible size occurred during periods of increasing density, and the magnitude of changes in skull and mandible shape was correlated with the magnitude of density changes. Frequencies of transferrin alleles changed the most when population density was in decline. Changes among the five phenotypic traits were generally uncorrelated with one another, except for skull and mandible shape. Nonmetric traits were selectively neutral when assessed with Q_ST/F_ST analysis, whereas mandible size, mandible shape, and skull shape showed evidence of fairly strong selection. Selection on skull size was weak or nonexistent. We discuss how different assumptions about the genetic components of variance affect Q_ST estimates when phenotypic variances are substituted for genetic ones. We also found that change in mandible size, mandible shape, skull size, and skull shape were greater than expected under a neutral model given reasonable assumptions about heritability and effective population size.     

Drosophila wing modularity revisited through a quantitative genetic approach

To predict the response of complex morphological structures to selection it is necessary to know how the covariation among its different parts is organized. Two key features of covariation are modularity and integration. The Drosophila wing is currently considered a fully integrated structure. Here, we study the patterns of integration of the Drosophila wing and test the hypothesis of the wing being divided into two modules along the proximo‐distal axis, as suggested by developmental, biomechanical, and evolutionary evidence. To achieve these goals we perform a multilevel analysis of covariation combining the techniques of geometric morphometrics and quantitative genetics. Our results indicate that the Drosophila wing is indeed organized into two main modules, the wing base and the wing blade. The patterns of integration and modularity were highly concordant at the phenotypic, genetic, environmental, and developmental levels. Besides, we found that modularity at the developmental level was considerably higher than modularity at other levels, suggesting that in the Drosophila wing direct developmental interactions are major contributors to total phenotypic shape variation. We propose that the precise time at which covariance‐generating developmental processes occur and/or the magnitude of variation that they produce favor proximo‐distal, rather than anterior‐posterior, modularity in the Drosophila wing.     

Geometric morphometrics and leaf phenotypic plasticity: assessing fluctuating asymmetry and allometry in European white oaks (Quercus)

Because plants are unable to move away from unfavourable habitats and environmental perturbations, leaf phenotypic plasticity facilitates light absorption and gas exchange. Oaks (Quercus spp.) are particularly known for their adaptability and plastic phenotypes, and leaf allometry and developmental instability may represent important mechanisms for their adaptation to environments and evolution. Because of its important role in the adaptation of plant populations to different environments, allometry can be involved in diversifying selection. Developmental instability is related to environmental perturbations and stresses by producing random deviations in structures characterized by bilateral symmetry, such as oak leaves. In addition, developmental instability can also arise from genetic bottlenecks or as a result of hybridization. The splitting of symmetric and asymmetric components of variation and their separate analysis allows the variability in leaf shape traits to be summarized, reducing the variation produced by developmental instability. The geometric morphometric approach is a useful method for the study of leaf asymmetry and allometric patterns. This method provides an important tool for the visualization of shape attributes that characterize species with highly variable leaf phenotypic patterns. In this study, leaf shape and size variability of three white oak species was investigated by means of a two‐dimensional landmark‐based method providing improved knowledge of variance partitioning, species discrimination, fluctuating asymmetry and allometric patterns of variation resulting from the different analyses. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 335–348.     

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype−environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.     

Genetics of body shape and armour variation in threespine sticklebacks

Patterns of genetic variation and covariation can influence the rate and direction of phenotypic evolution. We explored the possibility that the parallel morphological evolution seen in threespine stickleback (Gasterosteus aculeatus) populations colonizing freshwater environments is facilitated by patterns of genetic variation and covariation in the ancestral (marine) population. We estimated the genetic (G) and phenotypic (P) covariance matrices and directions of maximum additive genetic (g(max) ) and phenotypic (p(max) ) covariances of body shape and armour traits. Our results suggest a role for the ancestral G in explaining parallel morphological evolution in freshwater populations. We also found evidence of genetic constraints owing to the lack of variance in the ancestral G. Furthermore, strong genetic covariances and correlations among traits revealed that selective factors responsible for threespine stickleback body shape and armour divergence may be difficult to disentangle. The directions of g(max) and p(max) were correlated, but the correlations were not high enough to imply that phenotypic patterns of trait variation and covariation within populations are very informative of underlying genetic patterns.     

Developmental interactions and the constituents of quantitative variation

Development is the process by which genotypes are transformed into phenotypes. Consequently, development determines the relationship between allelic and phenotypic variation in a population and, therefore, the patterns of quantitative genetic variation and covariation of traits. Understanding the developmental basis of quantitative traits may lead to insights into the origin and evolution of quantitative genetic variation, the evolutionary fate of populations, and, more generally, the relationship between development and evolution. Herein, we assume a hierarchical, modular structure of trait development and consider how epigenetic interactions among modules during ontogeny affect patterns of phenotypic and genetic variation. We explore two developmental models, one in which the epigenetic interactions between modules result in additive effects on character expression and a second model in which these epigenetic interactions produce nonadditive effects. Using a phenotype landscape approach, we show how changes in the developmental processes underlying phenotypic expression can alter the magnitude and pattern of quantitative genetic variation. Additive epigenetic effects influence genetic variances and covariances, but allow trait means to evolve independently of the genetic variances and covariances, so that phenotypic evolution can proceed without changing the genetic covariance structure that determines future evolutionary response. Nonadditive epigenetic effects, however, can lead to evolution of genetic variances and covariances as the mean phenotype evolves. Our model suggests that an understanding of multivariate evolution can be considerably enriched by knowledge of the mechanistic basis of character development.