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.     

ACROSS‐ENVIRONMENT GENETIC CORRELATIONS AND THE FREQUENCY OF SELECTIVE ENVIRONMENTS SHAPE THE EVOLUTIONARY DYNAMICS OF GROWTH RATE IN IMPATIENS CAPENSIS

Trade‐offs can exist within and across environments, and constrain evolutionary trajectories. To examine the effects of competition and resource availability on trade‐offs, we grew individuals of recombinant inbred lines of Impatiens capensis in a factorial combination of five densities with two light environments (full light and neutral shade) and used a Bayesian logistic growth analysis to estimate intrinsic growth rates. To estimate across‐environment constraints, we developed a variance decomposition approach to principal components analysis, which accounted for sample size, model‐fitting, and within‐RIL variation prior to eigenanalysis. We detected negative across‐environment genetic covariances in intrinsic growth rates, although only under full‐light. To evaluate the potential importance of these covariances, we surveyed natural populations of I. capensis to measure the frequency of different density environments across space and time. We combined our empirical estimates of across‐environment genetic variance–covariance matrices and frequency of selective environments with hypothetical (yet realistic) selection gradients to project evolutionary responses in multiple density environments. Selection in common environments can lead to correlated responses to selection in rare environments that oppose and counteract direct selection in those rare environments. Our results highlight the importance of considering both the frequency of selective environments and the across‐environment genetic covariances in traits simultaneously.     

RESISTANCE OF GENETIC CORRELATION STRUCTURE TO DIRECTIONAL SELECTION IN DROSOPHILA MELANOGASTER

The genetic covariance and correlation matrices for five morphological traits were estimated from four populations of fruit flies, Drosophila melanogaster, to measure the extent of change in genetic covariances as a result of directional selection. Two of the populations were derived from lines that had undergone selection for large or small thorax length over the preceding 23 generations. A third population was constituted using flies from control lines that were maintained with equivalent population sizes as the selected lines. The fourth population contained flies from the original cage population from which the selected and control lines had been started. Tests of the homogeneity of covariance matrices using maximum likelihood techniques revealed significant changes in covariance structure among the selected lines. Prediction of base population trait means from selected line means under the assumption of constant genetic covariances indicated that genetic covariances for the small population differed more from the base population than did the covariances for the large population. The predicted small population means diverged farther from the expected means because the additive genetic variance associated with several traits increased in value and most of the genetic covariances associated with one trait changed in sign. These results illustrate that genetic covariances may remain nearly constant in some situations while changing markedly in others. Possible developmental reasons for the genetic changes are discussed.     

A MULTIVARIATE ANALYSIS OF GENETIC VARIATION IN THE ADVERTISEMENT CALL OF THE GRAY TREEFROG,HYLA VERSICOLOR

Genetic variation in sexual displays is crucial for an evolutionary response to sexual selection, but can be eroded by strong selection. Identifying the magnitude and sources of additive genetic variance underlying sexually selected traits is thus an important issue in evolutionary biology. We conducted a quantitative genetics experiment with gray treefrogs (Hyla versicolor) to investigate genetic variances and covariances among features of the male advertisement call. Two energetically expensive traits showed significant genetic variation: call duration, expressed as number of pulses per call, and call rate, represented by its inverse, call period. These two properties also showed significant genetic covariance, consistent with an energetic constraint to call production. Combining the genetic variance–covariance matrix with previous estimates of directional sexual selection imposed by female preferences predicts a limited increase in call duration but no change in call rate despite significant selection on both traits. In addition to constraints imposed by the genetic covariance structure, an evolutionary response to sexual selection may also be limited by high energetic costs of long‐duration calls and by preferences that act most strongly against very short‐duration calls. Meanwhile, the persistence of these preferences could be explained by costs of mating with males with especially unattractive calls.     

HOMOGENEITY OF THE GENETIC VARIANCE-COVARIANCE MATRIX FOR ANTIPREDATOR TRAITS IN TWO NATURAL POPULATIONS OF THE GARTER SNAKE THAMNOPHIS ORDINOIDES

Quantitative genetic models of evolution rely on the genetic variance-covariance matrix to predict the phenotypic response to selection. Both prospective and retrospective studies of phenotypic evolution across generations rely on assumptions about the constancy of patterns of genetic covariance through time. In the absence of robust theoretical predictions about the stability of genetic covariances, this assumption must be tested with empirical comparisons of genetic parameters among populations and species. Genetic variance-covariance matrices were estimated for a suite of antipredator traits in two populations of the northwestern garter snake, Thamnophis ordinoides. The characters studied include color pattern and antipredator behaviors that interact to facilitate escape from predators. Significant heritabilities for all traits were detected in both populations. Genetic correlations and covariances were found among behaviors in both populations and between color pattern and behavior in one of the populations. Phenotypic means differed among populations, but pairwise comparisons revealed no heterogeneity of genetic parameters between the populations. The structure of the genetic variance-covariance matrix has apparently not changed significantly during the divergence of these two populations.     

Evolutionary quantitative genetics of nonlinear developmental systems

In quantitative genetics, the effects of developmental relationships among traits on microevolution are generally represented by the contribution of pleiotropy to additive genetic covariances. Pleiotropic additive genetic covariances arise only from the average effects of alleles on multiple traits, and therefore the evolutionary importance of nonlinearities in development is generally neglected in quantitative genetic views on evolution. However, nonlinearities in relationships among traits at the level of whole organisms are undeniably important to biology in general, and therefore critical to understanding evolution. I outline a system for characterizing key quantitative parameters in nonlinear developmental systems, which yields expressions for quantities such as trait means and phenotypic and genetic covariance matrices. I then develop a system for quantitative prediction of evolution in nonlinear developmental systems. I apply the system to generating a new hypothesis for why direct stabilizing selection is rarely observed. Other uses will include separation of purely correlative from direct and indirect causal effects in studying mechanisms of selection, generation of predictions of medium‐term evolutionary trajectories rather than immediate predictions of evolutionary change over single generation time‐steps, and the development of efficient and biologically motivated models for separating additive from epistatic genetic variances and covariances.     

QUANTITATIVE GENETIC DIVERGENCE AND STANDING GENETIC (CO)VARIANCE IN THERMAL REACTION NORMS ALONG LATITUDE

Although the potential to adapt to warmer climate is constrained by genetic trade‐offs, our understanding of how selection and mutation shape genetic (co)variances in thermal reaction norms is poor. Using 71 isofemale lines of the fly Sepsis punctum, originating from northern, central, and southern European climates, we tested for divergence in juvenile development rate across latitude at five experimental temperatures. To investigate effects of evolutionary history in different climates on standing genetic variation in reaction norms, we further compared genetic (co)variances between regions. Flies were reared on either high or low food resources to explore the role of energy acquisition in determining genetic trade‐offs between different temperatures. Although the latter had only weak effects on the strength and sign of genetic correlations, genetic architecture differed significantly between climatic regions, implying that evolution of reaction norms proceeds via different trajectories at high latitude versus low latitude in this system. Accordingly, regional genetic architecture was correlated to region‐specific differentiation. Moreover, hot development temperatures were associated with low genetic variance and stronger genetic correlations compared to cooler temperatures. We discuss the evolutionary potential of thermal reaction norms in light of their underlying genetic architectures, evolutionary histories, and the materialization of trade‐offs in natural environments.     

PHENOTYPIC PLASTICITY AND EPIGENETIC MARKING: AN ASSESSMENT OF EVIDENCE FOR GENETIC ACCOMMODATION

The relationship between genotype (which is inherited) and phenotype (the target of selection) is mediated by environmental inputs on gene expression, trait development, and phenotypic integration. Phenotypic plasticity or epigenetic modification might influence evolution in two general ways: (1) by stimulating evolutionary responses to environmental change via population persistence or by revealing cryptic genetic variation to selection, and (2) through the process of genetic accommodation, whereby natural selection acts to improve the form, regulation, and phenotypic integration of novel phenotypic variants. We provide an overview of models and mechanisms for how such evolutionary influences may be manifested both for plasticity and epigenetic marking. We point to promising avenues of research, identifying systems that can best be used to address the role of plasticity in evolution, as well as the need to apply our expanding knowledge of genetic and epigenetic mechanisms to our understanding of how genetic accommodation occurs in nature. Our review of a wide variety of studies finds widespread evidence for evolution by genetic accommodation.     

MUTATIONAL CONTRIBUTIONS TO GENETIC VARIANCE-COVARIANCE MATRICES: AN EXPERIMENTAL APPROACH USING INDUCED MUTATIONS IN ARABIDOPSIS THALIANA

Genetic potential for evolutionary change and covariational constraints are typically summarized as the genetic variance-covariance matrix G , and there is currently debate over the extent to which G remains effectively constant during the course of adaptive evolution. However, G provides only a temporally restricted view of constraints that ignores possible biases in how new mutations affect multivariate phenotypes. We used chemical mutagenesis to study the effect of mutations as summarized by the mutational covariance matrix, M , in Arabidopsis thaliana. By introducing mutations into three isogenic strains of A. thaliana, we were able to quantify M directly as the genetic variance-covariance matrix of mutagenized lines. Induced mutations generally did not alter the means of the six morphology and life-history traits we measured, but they did affect the levels of available genetic variation and the covariances among traits. However, these effects were not consistent among the three isogenic lines; that is, there were significant differences among the lines in both the number of mutations produced by ethyl-methane-sulfonate treatment and the M matrices they induced. The evolutionary implications of the dependence of M on the number of mutations, the particular genetic background, and the mutagenic sampling of loci in the genome are discussed in light of commonly applied models of multivariate evolution and the potential for the genetic architecture itself to change in ways that facilitate the coordinated evolution of complex phenotypes.     

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.     

Variation in pleiotropy and the mutational underpinnings of the G-matrix

The pattern and extent of pleiotropic gene action can contribute substantially to the internal structure and shape of the additive genetic variance-covariance matrix (G)--a key determinant of evolutionary trajectories. We use data from our study (Estes et al. 2004) on the univariate effects of mutation in a mismatch-repair-defective strain, msh-2, of Caenorhabditis elegans to address the impact of increasing levels of selection on the magnitude and pattern of genetic covariance due to new mutations. Mutational covariances between three life-history traits are shown to exhibit a weak pattern of decline with increasing population size (increasing selection), while the orientation of mutational matrices remains reasonably constant. This suggests that mutations with smaller effects on fitness may tend to be slightly more confined in their influence than large-effect mutations (i.e., small-effect mutations reduce the magnitude of covariation between characters), but do not change the direction of this covariation.     

POPULATION STRUCTURE OF MORPHOLOGICAL TRAITS IN CLARKIA DUDLEYANA. II. CONSTANCY OF WITHIN-POPULATION GENETIC VARIANCE

Recent quantitative genetic studies have attempted to infer long-term selection responsible for differences in observed phenotypes. These analyses are greatly simplified by the assumption that the within-population genetic variance remains constant through time and over space, or for the multivariate case, that the matrix of additive genetic variances and covariances (G matrix) is constant. We examined differences in G matrices and the association of these differences with differences in multivariate means (Mahalanobis D²) among 11 populations of the California endemic annual plant, Clarkia dudleyana. Based on nine continuous morphological traits, the relationship between Mahalanobis D² and a distance measure summarizing differences in G matrices reflected no concomitant change in (co)variances with changes in means. Based on both broad- and narrow-sense analyses, we found little evidence that G matrices differed between populations. These results suggest that both the additive and nonadditive (co)variances for traits have remained relatively constant despite changes in means.     

EVOLUTION OF RESISTANCE TO A MULTIPLE‐HERBIVORE COMMUNITY: GENETIC CORRELATIONS,DIFFUSE COEVOLUTION,AND CONSTRAINTS ON THE PLANT'S RESPONSE TO SELECTION

Although plants are generally attacked by a community of several species of herbivores, relatively little is known about the strength of natural selection for resistance in multiple‐herbivore communities—particularly how the strength of selection differs among herbivores that feed on different plant organs or how strongly genetic correlations in resistance affect the evolutionary responses of the plant. Here, we report on a field study measuring natural selection for resistance in a diverse community of herbivores of Solanum carolinense. Using linear phenotypic‐selection analyses, we found that directional selection acted to increase resistance to seven species. Selection was strongest to increase resistance to fruit feeders, followed by flower feeders, then leaf feeders. Selection favored a decrease in resistance to a stem borer. Bootstrapping analyses showed that the plant population contained significant genetic variation for each of 14 measured resistance traits and significant covariances in one‐third of the pairwise combinations of resistance traits. These genetic covariances reduced the plant's overall predicted evolutionary response for resistance against the herbivore community by about 60%. Diffuse (co)evolution was widespread in this community, and the diffuse interactions had an overwhelmingly constraining (rather than facilitative) effect on the plant's evolution of resistance.     

ADAPTIVE RADIATION ALONG GENETIC LINES OF LEAST RESISTANCE

Are measurements of quantitative genetic variation useful for predicting long-term adaptive evolution? To answer this question, I focus on g_max, the multivariate direction of greatest additive genetic variance within populations. Original data on threespine sticklebacks, together with published genetic measurements from other vertebrates, show that morphological differentiation between species has been biased in the direction of g_max for at least four million years, despite evidence that natural selection is the cause of differentiation. This bias toward the direction of evolution tends to decay with time. Rate of morphological divergence between species is inversely proportional to θ, the angle between the direction of divergence and the direction of greatest genetic variation. The direction of greatest phenotypic variance is not identical with g_max, but for these data is nearly as successful at predicting the direction of species divergence. I interpret the findings to mean that genetic variances and covariances constrain adaptive change in quantitative traits for reasonably long spans of time. An alternative hypothesis, however, cannot be ruled out: that morphological differentiation is biased in the direction g_max because divergence and g_max are both shaped by the same natural selection pressures. Either way, the results reveal that adaptive differentiation occurs principally along “genetic lines of least resistance.”     

The effect of temperature and wing morphology on quantitative genetic variation in the cricket Gryllus firmus, with an appendix examining the statistical properties of the Jackknife-MANOVA method of matrix comparison

We investigated the effect of temperature and wing morphology on the quantitative genetic variances and covariances of five size-related traits in the sand cricket, Gryllus firmus. Micropterous and macropterous crickets were reared in the laboratory at 24, 28 and 32 degrees C. Quantitative genetic parameters were estimated using a nested full-sib family design, and (co)variance matrices were compared using the T method, Flury hierarchy and Jackknife-manova method. The results revealed that the mean phenotypic value of each trait varied significantly among temperatures and wing morphs, but temperature reaction norms were not similar across all traits. Micropterous individuals were always smaller than macropterous individuals while expressing more phenotypic variation, a finding discussed in terms of canalization and life-history trade-offs. We observed little variation between the matrices of among-family (co)variation corresponding to each combination of temperature and wing morphology, with only one matrix of six differing in structure from the others. The implications of this result are discussed with respect to the prediction of evolutionary trajectories.     

PATERNAL CARE: DIRECT AND INDIRECT GENETIC EFFECTS OF FATHERS ON OFFSPRING PERFORMANCE

Knowledge of how genetic effects arising from parental care influence the evolution of offspring traits comes almost exclusively from studies of maternal care. However, males provide care in some taxa, and often this care differs from females in quality or quantity. If variation in paternal care is genetically based then, like maternal care and maternal effects, paternal effects may have important consequences for the evolution of offspring traits via indirect genetic effects (IGEs). IGEs and direct–indirect genetic covariances associated with parental care can contribute substantially to total heritability and influence predictions about how traits respond to selection. It is unknown, however, if the magnitude and sign of parental effects arising from fathers are the same as those arising from mothers. We used a reciprocal cross‐fostering experiment to quantify environmental and genetic effects of paternal care on offspring performance in the burying beetle, Nicrophorus vespilloides. We found that IGEs were substantial and direct–indirect genetic covariances were negative. Combined, these patterns led to low total heritabilities for offspring performance traits. Thus, under paternal care, offspring performance traits are unlikely to evolve in response to selection, and variation in these traits will be maintained in the population despite potentially strong selection on these traits. These patterns are similar to those generated by maternal care, indicating that the genetic effects of care on offspring performance are independent of the caregiver's sex.     

VARIATION IN GENETIC ARCHITECTURE OF CALLING SONG AMONG POPULATIONS OF ALLONEMOBIUS SOCIUS,A. FASCIATUS,AND A HYBRID POPULATION: DRIFT OR SELECTION?

Predictions using quantitative genetic models generally assume that the variance-covariance matrices remain constant over time. This assumption is based on the supposition that selection is generally weak and hence variation lost through selection can be replaced by new mutations. Whether this is generally true can only be ascertained from empirical studies. Ideally for such a study we should be able to make a prediction concerning the relative strength of selection versus genetic drift. If the latter force is prevalent then the variance-covariances matrices should be proportional to each other. Previous studies have indicated that females in the two sibling cricket species Allonemobius socius and A. fasciatus do not discriminate between males of the two species by their calling song. Therefore, differences between the calling song of the two males most likely result from drift rather than sexual selection. We test this hypothesis by comparing the genetic architecture of calling song of three populations of A. fasciatus with two populations of A. socius. We found no differences among populations within species, but significant differences in the G (genetic) and P (phenotypic) matrices between species, with the matrices being proportional as predicted under the hypothesis of genetic drift. Because of the proportional change in the (co)variances no differences between species are evident in the heritabilities or genetic correlations. Comparison of the two species with a hybrid population from a zone of overlap showed highly significant nonproportional variation in genetic architecture. This variation is consistent with a general mixture of two separate genomes or selection. Qualitative conclusions reached using the phenotypic matrices are the same as those reached using the genetic matrices supporting the hypothesis that the former may be used as surrogate measures of the latter.     

A test for the genetic basis of natural selection: an individual-based longitudinal study in a stream-dwelling fish

In addition to the well-studied evolutionary parameters of (1) phenotype-fitness covariance and (2) the genetic basis of phenotypic variation, adaptive evolution by natural selection requires that (3) fitness variation is effected by heritable genetic differences among individuals and (4) phenotype-fitness covariances must be, at least in part, underlain by genetic covariances. These latter two requirements for adaptive evolutionary change are relatively unstudied in natural populations. Absence of the latter requirements could explain stasis of apparently directionally selected heritable traits. We provide complementary analyses of selection and variation at phenotypic and genetic levels for juvenile growth rate in brook charr Salvelinus fontinalis in Freshwater River, Newfoundland, Canada. Contrary to the vast majority of reports in fish, we found very little viability selection of juvenile body size. Large body size appears nonetheless to be selectively advantageous via a relationship with early maturity. Genetic patterns in evolutionary parameters largely reflected phenotypic patterns. We have provided inference of selection based on longitudinal data, which are uncommon in high fecundity organisms. Furthermore we have provided a practicable framework for further studies of the genetic basis of natural selection.     

A TEST OF THE CONJECTURE THAT G‐MATRICES ARE MORE STABLE THAN B‐MATRICES

The G‐matrix occupies an important position in evolutionary biology both as a summary of the inheritance of quantitative traits and as an ingredient in predicting how those traits will respond to selection and drift. Consequently, the stability of G has an important bearing on the accuracy of predicted evolutionary trajectories. Furthermore, G should evolve in response to stable features of the adaptive landscape and their trajectories through time. Although the stability and evolution of G might be predicted from knowledge of selection in natural populations, most empirical comparisons of G‐matrices have been made in the absence of such a priori predictions. We present a theoretical argument that within‐sex G‐matrices should be more stable than between‐sex B‐matrices because they are more powerfully exposed to multivariate stabilizing selection. We tested this conjecture by comparing estimates of B‐ and within‐sex G‐matrices among three populations of the garter snake Thamnophis elegans. Matrix comparisons using Flury's hierarchical approach revealed that within‐sex G‐matrices had four principal components in common (full CPC), whereas B‐matrices had only a single principal component in common and eigenvalues that were more variable among populations. These results suggest that within‐sex G is more stable than B , as predicted by our theoretical argument.     

MEASURING SELECTION AND CONSTRAINT IN THE EVOLUTION OF GROWTH

We present a quantitative genetic model for the evolution of growth trajectories that makes no assumptions about the shapes of growth trajectories that are possible. Evolution of a population's mean growth trajectory is governed by the selection gradient function and the additive genetic covariance function. The selection gradient function is determined by the impact of changes in size on the birth and death rates at different ages, and can be estimated for natural populations. The additive genetic covariance function can also be estimated empirically, as we demonstrate with four vertebrate populations. Using the genetic data from mice, a computer simulation shows that evolution of a growth trajectory can be constrained by the absence of genetic variation for certain changes in the trajectory's shape. These constraints can be visualized with an analysis of the covariance function. Results from four vertebrate populations show that while each has substantial genetic variation for some evolutionary changes in its growth trajectory, most types of changes have little or no variation available. This suggests that constraints may often play an important role in the evolution of growth.