From micro- to macroevolution through quantitative genetic variation: positive evidence from field crickets

Quantitative genetics has been introduced to evolutionary biologists with the suggestion that microevolution could be directly linked to macroevolutionary patterns using, among other parameters, the additive genetic variance/ covariance matrix (G) which is a statistical representation of genetic constraints to evolution. However, little is known concerning the rate and pattern of evolution of G in nature, and it is uncertain whether the constraining effect of G is important over evolutionary time scales. To address these issues, seven species of field crickets from the genera Gryllus and Teleogryllus were reared in the laboratory, and quantitative genetic parameters for morphological traits were estimated from each of them using a nested full-sibling family design. We used three statistical approaches (T method, Flury hierarchy, and Mantel test) to compare G matrices or genetic correlation matrices in a phylogenetic framework. Results showed that G matrices were generally similar across species, with occasional differences between some species. We suggest that G has evolved at a low rate, a conclusion strengthened by the consideration that part of the observed across-species variation in G can be explained by the effect of a genotype by environment interaction. The observed pattern of G matrix variation between species could not be predicted by either morphological trait values or phylogeny. The constraint hypothesis was tested by comparing the multivariate orientation of the reconstructed ancestral G matrix to the orientation of the across-species divergence matrix (D matrix, based on mean trait values). The D matrix mainly revealed divergence in size and, to a much smaller extent, in a shape component related to the ovipositor length. This pattern of species divergence was found to be predictable from the ancestral G matrix in agreement with the expectation of the constraint hypothesis. Overall, these results suggest that the G matrix seems to have an influence on species divergence, and that macroevolution can be predicted, at least qualitatively, from quantitative genetic theory. Alternative explanations are discussed.     

DOES GENE FLOW CONSTRAIN ADAPTIVE DIVERGENCE OR VICE VERSA? A TEST USING ECOMORPHOLOGY AND SEXUAL ISOLATION IN TIMEMA CRISTINAE WALKING‐STICKS

Population differentiation often reflects a balance between divergent natural selection and the opportunity for homogenizing gene flow to erode the effects of selection. However, during ecological speciation, trait divergence results in reproductive isolation and becomes a cause, rather than a consequence, of reductions in gene flow. To assess both the causes and the reproductive consequences of morphological differentiation, we examined morphological divergence and sexual isolation among 17 populations of Timema cristinae walking-sticks. Individuals from populations adapted to using Adenostoma as a host plant tended to exhibit smaller overall body size, wide heads, and short legs relative to individuals using Ceonothus as a host. However, there was also significant variation in morphology among populations within host-plant species. Mean trait values for each single population could be reliably predicted based upon host-plant used and the potential for homogenizing gene flow, inferred from the size of the neighboring population using the alternate host and mitochondrial DNA estimates of gene flow. Morphology did not influence the probability of copulation in between-population mating trials. Thus, morphological divergence is facilitated by reductions in gene flow, but does not cause reductions in gene flow via the evolution of sexual isolation. Combined with rearing data indicating that size and shape have a partial genetic basis, evidence for parallel origins of the host-associated forms, and inferences from functional morphology, these results indicate that morphological divergence in T. cristinae reflects a balance between the effects of host-specific natural selection and gene flow. Our findings illustrate how data on mating preferences can help determine the causal associations between trait divergence and levels of gene flow.     

Evolutionary rates for multivariate traits: the role of selection and genetic variation

A fundamental question in evolutionary biology is the relative importance of selection and genetic architecture in determining evolutionary rates. Adaptive evolution can be described by the multivariate breeders'' equation (), which predicts evolutionary change for a suite of phenotypic traits () as a product of directional selection acting on them (β) and the genetic variance–covariance matrix for those traits (G). Despite being empirically challenging to estimate, there are enough published estimates of G and β to allow for synthesis of general patterns across species. We use published estimates to test the hypotheses that there are systematic differences in the rate of evolution among trait types, and that these differences are, in part, due to genetic architecture. We find some evidence that sexually selected traits exhibit faster rates of evolution compared with life-history or morphological traits. This difference does not appear to be related to stronger selection on sexually selected traits. Using numerous proposed approaches to quantifying the shape, size and structure of G, we examine how these parameters relate to one another, and how they vary among taxonomic and trait groupings. Despite considerable variation, they do not explain the observed differences in evolutionary rates.     

Environmental effects on the structure of the G‐matrix

Genetic correlations between traits determine the multivariate response to selection in the short term, and thereby play a causal role in evolutionary change. Although individual studies have documented environmentally induced changes in genetic correlations, the nature and extent of environmental effects on multivariate genetic architecture across species and environments remain largely uncharacterized. We reviewed the literature for estimates of the genetic variance–covariance ( G ) matrix in multiple environments, and compared differences in G between environments to the divergence in G between conspecific populations (measured in a common garden). We found that the predicted evolutionary trajectory differed as strongly between environments as it did between populations. Between‐environment differences in the underlying structure of G (total genetic variance and the relative magnitude and orientation of genetic correlations) were equal to or greater than between‐population differences. Neither environmental novelty, nor the difference in mean phenotype predicted these differences in G . Our results suggest that environmental effects on multivariate genetic architecture may be comparable to the divergence that accumulates over dozens or hundreds of generations between populations. We outline avenues of future research to address the limitations of existing data and characterize the extent to which lability in genetic correlations shapes evolution in changing environments.     

Measuring the contribution of evolution to community trait structure in freshwater zooplankton

There are currently few predictions about when evolutionary processes are likely to play an important role in structuring community features. Determining predictors that indicate when evolution is expected to impact ecological processes in natural landscapes can help researchers identify eco-evolutionary ‘hotspots', where eco-evolutionary interactions are more likely to occur. Using data collected from a survey in freshwater cladoceran communities, landscape population genetic data and phenotypic trait data measured in a common garden, we applied a Bayesian linear model to assess whether the impact of local trait evolution in the keystone species Daphnia magna on cladoceran community trait values could be predicted by population genetic properties (within-population genetic diversity, genetic distance among populations), ecological properties (Simpson's diversity, phenotypic divergence) or environmental divergence. We found that the impact of local trait evolution varied among communities. Moreover, community diversity and phenotypic divergence were found to be better predictors of the contribution of evolution to community trait values than environmental features or genetic properties of the evolving species. Our results thus indicate the importance of ecological context for the impact of evolution on community features. Our study also demonstrates one way to detect signatures of eco-evolutionary interactions in communities inhabiting heterogeneous landscapes using survey data of contemporary ecological and evolutionary structure.     

Sexual selection and population divergence II. Divergence in different sexual traits and signal modalities in field crickets (Teleogryllus oceanicus)

Sexual selection can target many different types of traits. However, the relative influence of different sexually selected traits during evolutionary divergence is poorly understood. We used the field cricket Teleogryllus oceanicus to quantify and compare how five traits from each of three sexual signal modalities and components diverge among allopatric populations: male advertisement song, cuticular hydrocarbon (CHC) profiles and forewing morphology. Population divergence was unexpectedly consistent: we estimated the among‐population (genetic) variance‐covariance matrix, D , for all 15 traits, and D_max explained nearly two‐thirds of its variation. CHC and wing traits were most tightly integrated, whereas song varied more independently. We modeled the dependence of among‐population trait divergence on genetic distance estimated from neutral markers to test for signatures of selection versus neutral divergence. For all three sexual trait types, phenotypic variation among populations was largely explained by a neutral model of divergence. Our findings illustrate how phenotypic integration across different types of sexual traits might impose constraints on the evolution of mating isolation and divergence via sexual selection.     

Rapid divergence of genetic variance-covariance matrix within a natural population

The matrix of genetic variances and covariances (G matrix) represents the genetic architecture of multiple traits sharing developmental and genetic processes and is central for predicting phenotypic evolution. These predictions require that the G matrix be stable. Yet the timescale and conditions promoting G matrix stability in natural populations remain unclear. We studied stability of the G matrix in a 20-year evolution field experiment, where a population of the cosmopolitan parthenogenetic soil nematode Acrobeloides nanus was subjected to drift and divergent selection (benign and stress environments). Selection regime did not influence the level of absolute genetic constraints: under both regimes, two genetic dimensions for three life-history traits were identified. A substantial response to selection in principal components structure and in general matrix pattern was indicated by three statistical methods. G structure was also influenced by drift, with higher divergence under benign conditions. These results show that the G matrix might evolve rapidly in natural populations. The observed high dynamics of G structure probably represents the general feature of asexual species and limits the predictive power of G in phenotypic evolution analyses.     

TYPE I ERROR RATES FOR TESTING GENETIC DRIFT WITH PHENOTYPIC COVARIANCE MATRICES: A SIMULATION STUDY

Studies of evolutionary divergence using quantitative genetic methods are centered on the additive genetic variance–covariance matrix ( G ) of correlated traits. However, estimating G properly requires large samples and complicated experimental designs. Multivariate tests for neutral evolution commonly replace average G by the pooled phenotypic within‐group variance–covariance matrix ( W ) for evolutionary inferences, but this approach has been criticized due to the lack of exact proportionality between genetic and phenotypic matrices. In this study, we examined the consequence, in terms of type I error rates, of replacing average G by W in a test of neutral evolution that measures the regression slope between among‐population variances and within‐population eigenvalues (the Ackermann and Cheverud [AC] test) using a simulation approach to generate random observations under genetic drift. Our results indicate that the type I error rates for the genetic drift test are acceptable when using W instead of average G when the matrix correlation between the ancestral G and P is higher than 0.6, the average character heritability is above 0.7, and the matrices share principal components. For less‐similar G and P matrices, the type I error rates would still be acceptable if the ratio between the number of generations since divergence and the effective population size (t/N_e) is smaller than 0.01 (large populations that diverged recently). When G is not known in real data, a simulation approach to estimate expected slopes for the AC test under genetic drift is discussed.     

Phenotypic divergence of the common toad (Bufo bufo) along an altitudinal gradient: evidence for local adaptation

Variation in the environment can induce different patterns of genetic and phenotypic differentiation among populations. Both neutral processes and selection can influence phenotypic differentiation. Altitudinal phenotypic variation is of particular interest in disentangling the interplay between neutral processes and selection in the dynamics of local adaptation processes but remains little explored. We conducted a common garden experiment to study the phenotypic divergence in larval life-history traits among nine populations of the common toad (Bufo bufo) along an altitudinal gradient in France. We further used correlation among population pairwise estimates of quantitative trait (Q_ST) and neutral genetic divergence (F_ST from neutral microsatellite markers), as well as altitudinal difference, to estimate the relative role of divergent selection and neutral genetic processes in phenotypic divergence. We provided evidence for a neutral genetic differentiation resulting from both isolation by distance and difference in altitude. We found evidence for phenotypic divergence along the altitudinal gradient (faster development, lower growth rate and smaller metamorphic size). The correlation between pairwise Q_STs–F_STs and altitude differences suggested that this phenotypic differentiation was most likely driven by altitude-mediated selection rather than by neutral genetic processes. Moreover, we found different divergence patterns for larval traits, suggesting that different selective agents may act on these traits and/or selection on one trait may constrain the evolution on another through genetic correlation. Our study highlighted the need to design more integrative studies on the common toad to unravel the underlying processes of phenotypic divergence and its selective agents in the context of environmental clines.     

Effects of migration on the genetic covariance matrix

In 1996, Schluter showed that the direction of morphological divergence of closely related species is biased toward the line of least genetic resistance, represented by g_max , the leading eigenvector of the matrix of genetic variance–covariance (the G -matrix). G is used to predict the direction of evolutionary change in natural populations. However, this usage requires that G is sufficiently constant over time to have enough predictive significance. Here, we explore the alternative explanation that G can evolve due to gene flow to conform to the direction of divergence between incipient species. We use computer simulations in a mainland–island migration model with stabilizing selection on two quantitative traits. We show that a high level of gene flow from a mainland population is required to significantly affect the orientation of the G -matrix in an island population. The changes caused by the introgression of the mainland alleles into the island population affect all aspects of the shape of G (size, eccentricity, and orientation) and lead to the alignment of g_max with the line of divergence between the two populations' phenotypic optima. Those changes decrease with increased correlation in mutational effects and with a correlated selection. Our results suggest that high migration rates, such as those often seen at the intraspecific level, will substantially affect the shape and orientation of G , whereas low migration (e.g., at the interspecific level) is unlikely to substantially affect the evolution of G .     

The constancy of the G matrix through species divergence and the effects of quantitative genetic constraints on phenotypic evolution: a case study in crickets

Long-term phenotypic evolution can be modeled using the response-to-selection equation of quantitative genetics, which incorporates information about genetic constraints (the G matrix). However, little is known about the evolution of G and about its long-term importance in constraining phenotypic evolution. We first investigated the degree of conservation of the G matrix across three species of crickets and qualitatively compared the pattern of variation of G to the phylogeny of the group. Second, we investigated the effect of G on phenotypic evolution by comparing the direction of greatest quantitative genetic variation within species (g(max)) to the direction of phenotypic divergence between species (Delta(z)). Each species, Gryllus veletis, G. firmus, and G. pennsylvanicus, was reared in the laboratory using a full-sib breeding design to extract quantitative genetic information. Five morphological traits related to size were measured. G matrices were compared using three statistical approaches: the T method, the Flury hierarchy, and the MANOVA method. Results revealed that the differences between matrices were small and mostly caused by differences in the magnitude of the genetic variation, not by differences in principal component structure. This suggested that the G matrix structure of this group of species was preserved, despite significant phenotypic divergence across species. The small observed differences in G matrices across species were qualitatively consistent with genetic distances, whereas ecological information did not provide a good prediction of G matrix variation. The comparison of g(max) and Delta(z) revealed that the angle between these two vectors was small in two of three species comparisons, whereas the larger angle corresponding to the third species comparison was caused in large part by one of the five traits. This suggests that multivariate phenotypic divergence occurred mostly in a direction predicted by the direction of greatest genetic variation, although it was not possible to demonstrate the causal relationship from G to Delta(z). Overall, this study provided some support for the validity of the predictive power of quantitative genetics over evolutionary time scales.     

THE EVOLUTIONARY STABILITY OF CROSS‐SEX,CROSS‐TRAIT GENETIC COVARIANCES

Although knowledge of the selective agents behind the evolution of sexual dimorphism has advanced considerably in recent years, we still lack a clear understanding of the evolutionary durability of cross‐sex genetic covariances that often constrain its evolution. We tested the relative stability of cross‐sex genetic covariances for a suite of homologous contact pheromones of the fruit fly Drosophila serrata, along a latitudinal gradient where these traits have diverged in mean. Using a Bayesian framework, which allowed us to account for uncertainty in all parameter estimates, we compared divergence in the total amount and orientation of genetic variance across populations, finding divergence in orientation but not total variance. We then statistically compared orientation divergence of within‐sex ( G ) to cross‐sex ( B ) covariance matrices. In line with a previous theoretical prediction, we find that the cross‐sex covariance matrix, B , is more variable than either within‐sex G matrix. Decomposition of B matrices into their symmetrical and nonsymmetrical components revealed that instability is linked to the degree of asymmetry. We also find that the degree of asymmetry correlates with latitude suggesting a role for spatially varying natural selection in shaping genetic constraints on the evolution of sexual dimorphism.     

Genetic covariances promote climatic adaptation in Australian Drosophila*

Evolutionary potential for adaptation hinges upon the orientation of genetic variation for traits under selection, captured by the additive genetic variance-covariance matrix (G), as well as the evolutionary stability of G. Yet studies that assess both the stability of G and its alignment with selection are extraordinarily rare. We evaluated the stability of G in three Drosophila melanogaster populations that have adapted to local climatic conditions along a latitudinal cline. We estimated population- and sex-specific G matrices for wing size and three climatic stress-resistance traits that diverge adaptively along the cline. To determine how G affects evolutionary potential within these populations, we used simulations to quantify how well G aligns with the direction of trait divergence along the cline (as a proxy for the direction of local selection) and how genetic covariances between traits and sexes influence this alignment. We found that G was stable across the cline, showing no significant divergence overall, or in sex-specific subcomponents, among populations. G also aligned well with the direction of clinal divergence, with genetic covariances strongly elevating evolutionary potential for adaptation to climatic extremes. These results suggest that genetic covariances between both traits and sexes should significantly boost evolutionary responses to environmental change.     

Genetic constraints on the evolution of mate recognition under natural selection

Field populations of Drosophila serrata display reproductive character displacement in cuticular hydrocarbons (CHCs) when sympatric with Drosophila birchii. We have previously shown that the naturally occurring pattern of reproductive character displacement can be experimentally replicated by exposing field allopatric populations of D. serrata to experimental sympatry with D. birchii. Here, we tested whether the repeated evolution of reproductive character displacement in natural and experimental populations was a consequence of genetic constraints on the evolution of CHCs. The genetic variance-covariance (G) matrices for CHCs were determined for populations of D. serrata that had evolved in either the presence or absence of D. birchii under field and experimental conditions. Natural selection on mate recognition under both field and experimental sympatric conditions increased the genetic variance in CHCs consistent with a response to selection based on rare alleles. A close association between G eigenstructure and the eigenstructure of the phenotypic divergence (D) matrix in natural and experimental populations suggested that G matrix eigenstructure may have determined the direction in which reproductive character displacement evolved during the reinforcement of mate recognition.     

Sexual dimorphism in the quantitative-genetic architecture of floral, leaf, and allocation traits in Silene latifolia

Many species exhibit sexual dimorphism in a variety of characters, and the underlying genetic architecture of dimorphism potentially involves sex-specific differences in the additive-genetic variance-covariance matrix (G) of dimorphic traits. We investigated the quantitative-genetic structure of dimorphic traits in the dioecious plant Silene latifolia by estimating G (including within-sex matrices, G(m), G(f), and the between-sex variance-covariance matrix, B), and the phenotypic variance-covariance matrix (P) for seven traits. Flower number was the most sexually dimorphic trait, and was significantly genetically correlated with all traits within each sex. Negative genetic correlations between flower size and number suggested a genetic trade-off in investment, but positive environmental correlations between the same traits resulted in no physical evidence for a trade-off in the phenotype. Between-sex genetic covariances for homologous traits were always greater than 0 but smaller than 1, showing that some, but not all, of the variation in traits is caused by genes or alleles with sex-limited expression. Using common principal-components analysis (CPCA), a maximum-likelihood (ML) estimation approach, and element-by-element comparison to compare matrices, we found that G(m) and G(f) differed significantly in eigenstructure because of dissimilarity in covariances involving leaf traits, suggesting the presence of variation in sex-limited genes with pleiotropic effects and/or linkage between sex-limited loci. The sex-specific structure of G is expected to cause differences in the correlated responses to selection within each sex, promoting the further evolution and maintenance of dimorphism.     

Measuring and comparing evolvability and constraint in multivariate characters

The Lande equation forms the basis for our understanding of the short-term evolution of quantitative traits in a multivariate context. It predicts the response to selection as the product of an additive genetic variance matrix and a selection gradient. The selection gradient approximates the force and direction of selection, and the genetic variance matrix quantifies the role of the genetic system in evolution. Attempts to understand the evolutionary significance of the genetic variance matrix are hampered by the fact that the majority of the methods used to characterize and compare variance matrices have not been derived in an explicit theoretical context. We use the Lande equation to derive new measures of the ability of a variance matrix to allow or constrain evolution in any direction in phenotype space. Evolvability captures the ability of a population to evolve in the direction of selection when stabilizing selection is absent. Conditional evolvability captures the ability of a population to respond to directional selection in the presence of stabilizing selection on other trait combinations. We then derive measures of character autonomy and integration from these evolvabilities. We study the properties of these measures and show how they can be used to interpret and compare variance matrices. As an illustration, we show that divergence of wing shape in the dipteran family Drosophilidae has proceeded in directions that have relatively high evolvabilities.     

Population divergence along lines of genetic variance and covariance in the invasive plant Lythrum salicaria in eastern North America

Evolution during biological invasion may occur over contemporary timescales, but the rate of evolutionary change may be inhibited by a lack of standing genetic variation for ecologically relevant traits and by fitness trade-offs among them. The extent to which these genetic constraints limit the evolution of local adaptation during biological invasion has rarely been examined. To investigate genetic constraints on life-history traits, we measured standing genetic variance and covariance in 20 populations of the invasive plant purple loosestrife (Lythrum salicaria) sampled along a latitudinal climatic gradient in eastern North America and grown under uniform conditions in a glasshouse. Genetic variances within and among populations were significant for all traits; however, strong intercorrelations among measurements of seedling growth rate, time to reproductive maturity and adult size suggested that fitness trade-offs have constrained population divergence. Evidence to support this hypothesis was obtained from the genetic variance-covariance matrix (G) and the matrix of (co)variance among population means (D), which were 79.8% (95% C.I. 77.7-82.9%) similar. These results suggest that population divergence during invasive spread of L. salicaria in eastern North America has been constrained by strong genetic correlations among life-history traits, despite large amounts of standing genetic variation for individual traits.     

Phenotypic divergence along lines of genetic variance

Natural populations inhabiting the same environment often independently evolve the same phenotype. Is this replicated evolution a result of genetic constraints imposed by patterns of genetic covariation? We looked for associations between directions of morphological divergence and the orientation of the genetic variance-covariance matrix (G) by using an experimental system of morphological evolution in two allopatric nonsister species of rainbow fish. Replicate populations of both Melanotaenia eachamensis and Melanotaenia duboulayi have independently adapted to lake versus stream hydrodynamic environments. The major axis of divergence (z) among all eight study populations was closely associated with the direction of greatest genetic variance (gmax), suggesting directional genetic constraint on evolution. However, the direction of hydrodynamic adaptation was strongly associated with vectors of G describing relatively small proportions of the total genetic variance, and was only weakly associated with gmax. In contrast, divergence between replicate populations within each habitat was approximately proportional to the level of genetic variance, a result consistent with theoretical predictions for neutral phenotypic divergence. Divergence between the two species was also primarily along major eigenvectors of G. Our results therefore suggest that hydrodynamic adaptation in rainbow fish was not directionally constrained by the dominant eigenvector of G. Without partitioning divergence as a consequence of the adaptation of interest (here, hydrodynamic adaptation) from divergence due to other processes, empirical studies are likely to overestimate the potential for the major eigenvectors of G to directionally constrain adaptive evolution.     

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.     

Evolution and stability of the G-matrix during the colonization of a novel environment

Populations that undergo a process of rapid evolution present excellent opportunities to investigate the mechanisms driving or restraining adaptive divergence. The genetic variance-covariance matrix (G) is often considered to constrain adaptation but little is known about its potential to evolve during phenotypic divergence. We compared the G-matrices of ancestral and recently established ecotype populations of an aquatic isopod (Asellus aquaticus) that have diverged in parallel in two south Swedish lakes. Phenotypic changes after colonization involved a reduction in overall size, lost pigmentation and changes in shape. Comparisons between G-matrices reveal close similarity within the same ecotype from different lakes but some degree of differentiation among ecotypes. Phenotypic divergence has apparently not been much influenced by the orientation of G. Additive genetic variation in the newly colonized habitats has also decreased substantially. This suggests that a process of adaptation from standing genetic variation has occurred and has probably facilitated phenotypic divergence.