Evolution of quantitative traits in the wild: mind the ecology

Recent advances in the quantitative genetics of traits in wild animal populations have created new interest in whether natural selection, and genetic response to it, can be detected within long-term ecological studies. However, such studies have re-emphasized the fact that ecological heterogeneity can confound our ability to infer selection on genetic variation and detect a population''s response to selection by conventional quantitative genetics approaches. Here, I highlight three manifestations of this issue: counter gradient variation, environmentally induced covariance between traits and the correlated effects of a fluctuating environment. These effects are symptomatic of the oversimplifications and strong assumptions of the breeder''s equation when it is applied to natural populations. In addition, methods to assay genetic change in quantitative traits have overestimated the precision with which change can be measured. In the future, a more conservative approach to inferring quantitative genetic response to selection, or genomic approaches allowing the estimation of selection intensity and responses to selection at known quantitative trait loci, will provide a more precise view of evolution in ecological time.     

THE CONTRIBUTION OF MATERNAL EFFECTS TO SELECTION RESPONSE: AN EMPIRICAL TEST OF COMPETING MODELS

Maternal effects can dramatically influence the evolutionary process, in some cases facilitating and in others hindering adaptive evolution. Maternal effects have been incorporated into quantitative genetic models using two theoretical frameworks: the variance‐components approach, which partitions variance into direct and maternal components, and the trait‐based approach, which assumes that maternal effects are mediated by specific maternal traits. Here, we demonstrate parallels between these models and test their ability to predict evolutionary change. First, we show that the two approaches predict equivalent responses to selection in the absence of maternal effects mediated by traits that are themselves maternally influenced. We also introduce a correction factor that may be applied when such cascading maternal effects are present. Second, we use several maternal effect models, as well as the standard breeder's equation, to predict evolution in response to artificial selection on flowering time in American bellflower, Campanulastrum americanum. Models that included maternal effects made much more accurate predictions of selection response than the breeder's equation. Maternal effect models differed somewhat in their fit, with a version of the trait‐based model providing the best fit. We recommend fitting such trait‐based models when possible and appropriate to make the most accurate evolutionary predictions.     

Evolutionary potential of morphological traits across different life‐history stages

Despite accumulating examples of selection acting on heritable traits in the wild, predicted evolutionary responses are often different from observed phenotypic trends. Various explanations have been suggested for these mismatches. These include within‐individual changes across lifespan that can create important variation in genetic architecture of traits and selection acting on them, but also potential problems with the methodological approach used to predict evolutionary responses of traits. Here, we used an 8‐year data set on tree swallow (Tachycineta bicolor) to first assess the effects of differences among three nestling life‐history stages on the genetic (co)variances of two morphological traits (body mass and primary feather length) and the selection acting on them over three generations. We then estimated the evolutionary potential of these traits by predicting their evolutionary responses using the breeder's equation and the secondary theorem of selection approaches. Our results showed variation in strength and direction of selection and slight changes in trait variance across ages. Predicted evolutionary responses differed importantly between both approaches for half of the trait–age combinations we studied, suggesting the presence of environmentally induced correlations between focal traits and fitness possibly biasing breeder's equation predictions. Our results emphasize that predictions of evolutionary potential for morphological traits are likely to be highly variable, both in strength and direction, depending on the life stage and method used, thus mitigating our capacity to predict adaptation and persistence of wild populations.     

ESTIMATING UNCERTAINTY IN MULTIVARIATE RESPONSES TO SELECTION

Predicting the responses to natural selection is one of the key goals of evolutionary biology. Two of the challenges in fulfilling this goal have been the realization that many estimates of natural selection might be highly biased by environmentally induced covariances between traits and fitness, and that many estimated responses to selection do not incorporate or report uncertainty in the estimates. Here we describe the application of a framework that blends the merits of the Robertson–Price Identity approach and the multivariate breeder's equation to address these challenges. The approach allows genetic covariance matrices, selection differentials, selection gradients, and responses to selection to be estimated without environmentally induced bias, direct and indirect selection and responses to selection to be distinguished, and if implemented in a Bayesian‐MCMC framework, statistically robust estimates of uncertainty on all of these parameters to be made. We illustrate our approach with a worked example of previously published data. More generally, we suggest that applying both the Robertson–Price Identity and the multivariate breeder's equation will facilitate hypothesis testing about natural selection, genetic constraints, and evolutionary responses.     

Predicting evolutionary responses to selection on polyandry in the wild: additive genetic covariances with female extra-pair reproduction

The evolutionary forces that underlie polyandry, including extra-pair reproduction (EPR) by socially monogamous females, remain unclear. Selection on EPR and resulting evolution have rarely been explicitly estimated or predicted in wild populations, and evolutionary predictions are vulnerable to bias due to environmental covariances and correlated selection through unmeasured traits. However, evolutionary responses to (correlated) selection on any trait can be directly predicted as additive genetic covariances (cov_A) with appropriate components of relative fitness. I used comprehensive life-history, paternity and pedigree data from song sparrows (Melospiza melodia) to estimate cov_A between a female''s liability to produce extra-pair offspring and two specific fitness components: relative annual reproductive success (ARS) and survival to recruitment. All three traits showed non-zero additive genetic variance. Estimates of cov_A were positive, predicting evolution towards increased EPR, but 95% credible intervals overlapped zero. There was therefore no conclusive prediction of evolutionary change in EPR due to (correlated) selection through female ARS or recruitment. Negative environmental covariance between EPR and ARS would have impeded evolutionary prediction from phenotypic selection differentials. These analyses demonstrate an explicit quantitative genetic approach to predicting evolutionary responses to components of (correlated) selection on EPR that should be unbiased by environmental covariances and unmeasured traits.     

COMPARING STRENGTHS OF DIRECTIONAL SELECTION: HOW STRONG IS STRONG?

The fundamental equation in evolutionary quantitative genetics, the Lande equation, describes the response to directional selection as a product of the additive genetic variance and the selection gradient of trait value on relative fitness. Comparisons of both genetic variances and selection gradients across traits or populations require standardization, as both are scale dependent. The Lande equation can be standardized in two ways. Standardizing by the variance of the selected trait yields the response in units of standard deviation as the product of the heritability and the variance-standardized selection gradient. This standardization conflates selection and variation because the phenotypic variance is a function of the genetic variance. Alternatively, one can standardize the Lande equation using the trait mean, yielding the proportional response to selection as the product of the squared coefficient of additive genetic variance and the mean-standardized selection gradient. Mean-standardized selection gradients are particularly useful for summarizing the strength of selection because the mean-standardized gradient for fitness itself is one, a convenient benchmark for strong selection. We review published estimates of directional selection in natural populations using mean-standardized selection gradients. Only 38 published studies provided all the necessary information for calculation of mean-standardized gradients. The median absolute value of multivariate mean-standardized gradients shows that selection is on average 54% as strong as selection on fitness. Correcting for the upward bias introduced by taking absolute values lowers the median to 31%, still very strong selection. Such large estimates clearly cannot be representative of selection on all traits. Some possible sources of overestimation of the strength of selection include confounding environmental and genotypic effects on fitness, the use of fitness components as proxies for fitness, and biases in publication or choice of traits to study.     

INTERACTING PHENOTYPES AND THE EVOLUTIONARY PROCESS. III. SOCIAL EVOLUTION

Interactions among conspecifics influence social evolution through two distinct but intimately related paths. First, they provide the opportunity for indirect genetic effects (IGEs), where genes expressed in one individual influence the expression of traits in others. Second, interactions can generate social selection when traits expressed in one individual influence the fitness of others. Here, we present a quantitative genetic model of multivariate trait evolution that integrates the effects of both IGEs and social selection, which have previously been modeled independently. We show that social selection affects evolutionary change whenever the breeding value of one individual covaries with the phenotype of its social partners. This covariance can be created by both relatedness and IGEs, which are shown to have parallel roles in determining evolutionary response. We show that social selection is central to the estimation of inclusive fitness and derive a version of Hamilton's rule showing the symmetrical effects of relatedness and IGEs on the evolution of altruism. We illustrate the utility of our approach using altruism, greenbeards, aggression, and weapons as examples. Our model provides a general predictive equation for the evolution of social phenotypes that encompasses specific cases such as kin selection and reciprocity. The parameters can be measured empirically, and we emphasize the importance of considering both IGEs and social selection, in addition to relatedness, when testing hypotheses about social evolution.     

Testing for biases in selection on avian reproductive traits and partitioning direct and indirect selection using quantitative genetic models

Key life history traits such as breeding time and clutch size are frequently both heritable and under directional selection, yet many studies fail to document microevolutionary responses. One general explanation is that selection estimates are biased by the omission of correlated traits that have causal effects on fitness, but few valid tests of this exist. Here, we show, using a quantitative genetic framework and six decades of life‐history data on two free‐living populations of great tits Parus major, that selection estimates for egg‐laying date and clutch size are relatively unbiased. Predicted responses to selection based on the Robertson–Price Identity were similar to those based on the multivariate breeder's equation (MVBE), indicating that unmeasured covarying traits were not missing from the analysis. Changing patterns of phenotypic selection on these traits (for laying date, linked to climate change) therefore reflect changing selection on breeding values, and genetic constraints appear not to limit their independent evolution. Quantitative genetic analysis of correlational data from pedigreed populations can be a valuable complement to experimental approaches to help identify whether apparent associations between traits and fitness are biased by missing traits, and to parse the roles of direct versus indirect selection across a range of environments.     

Evolutionary stasis of a heritable morphological trait in a wild fish population despite apparent directional selection

Comparing observed versus theoretically expected evolutionary responses is important for our understanding of the evolutionary process, and for assessing how species may cope with anthropogenic change. Here, we document directional selection for larger female size in Atlantic salmon, using pedigree‐derived estimates of lifetime reproductive success as a fitness measure. We show the trait is heritable and, thus, capable of responding to selection. The Breeder's Equation, which predicts microevolution as the product of phenotypic selection and heritability, predicted evolution of larger size. This was at odds, however, with the observed lack of either phenotypic or genetic temporal trends in body size, a so‐called “paradox of stasis.” To investigate this paradox, we estimated the additive genetic covariance between trait and fitness, which provides a prediction of evolutionary change according to Robertson's secondary theorem of selection (STS) that is unbiased by missing variables. The STS prediction was consistent with the observed stasis. Decomposition of phenotypic selection gradients into genetic and environmental components revealed a potential upward bias, implying unmeasured factors that covary with trait and fitness. These results showcase the power of pedigreed, wild population studies—which have largely been limited to birds and mammals—to study evolutionary processes on contemporary timescales.     

Selection on skewed characters and the paradox of stasis

Observed phenotypic responses to selection in the wild often differ from predictions based on measurements of selection and genetic variance. An overlooked hypothesis to explain this paradox of stasis is that a skewed phenotypic distribution affects natural selection and evolution. We show through mathematical modeling that, when a trait selected for an optimum phenotype has a skewed distribution, directional selection is detected even at evolutionary equilibrium, where it causes no change in the mean phenotype. When environmental effects are skewed, Lande and Arnold's (1983) directional gradient is in the direction opposite to the skew. In contrast, skewed breeding values can displace the mean phenotype from the optimum, causing directional selection in the direction of the skew. These effects can be partitioned out using alternative selection estimates based on average derivatives of individual relative fitness, or additive genetic covariances between relative fitness and trait (Robertson–Price identity). We assess the validity of these predictions using simulations of selection estimation under moderate sample sizes. Ecologically relevant traits may commonly have skewed distributions, as we here exemplify with avian laying date — repeatedly described as more evolutionarily stable than expected — so this skewness should be accounted for when investigating evolutionary dynamics in the wild.     

Natural selection and quantitative genetics of life-history traits in Western women: a twin study

Whether contemporary human populations are still evolving as a result of natural selection has been hotly debated. For natural selection to cause evolutionary change in a trait, variation in the trait must be correlated with fitness and be genetically heritable and there must be no genetic constraints to evolution. These conditions have rarely been tested in human populations. In this study, data from a large twin cohort were used to assess whether selection will cause a change among women in a contemporary Western population for three life-history traits: age at menarche, age at first reproduction, and age at menopause. We control for temporal variation in fecundity (the "baby boom" phenomenon) and differences between women in educational background and religious affiliation. University-educated women have 35% lower fitness than those with less than seven years education, and Roman Catholic women have about 20% higher fitness than those of other religions. Although these differences were significant, education and religion only accounted for 2% and 1% of variance in fitness, respectively. Using structural equation modeling, we reveal significant genetic influences for all three life-history traits, with heritability estimates of 0.50, 0.23, and 0.45, respectively. However, strong genetic covariation with reproductive fitness could only be demonstrated for age at first reproduction, with much weaker covariation for age at menopause and no significant covariation for age at menarche. Selection may, therefore, lead to the evolution of earlier age at first reproduction in this population. We also estimate substantial heritable variation in fitness itself, with approximately 39% of the variance attributable to additive genetic effects, the remainder consisting of unique environmental effects and small effects from education and religion. We discuss mechanisms that could be maintaining such a high heritability for fitness. Most likely is that selection is now acting on different traits from which it did in pre-industrial human populations.     

DIRECT AND INDIRECT SEXUAL SELECTION AND QUANTITATIVE GENETICS OF MALE TRAITS IN GUPPIES (POECILIA RETICULATA)

Abstract.— The ornamentation and displays on which sexual attractiveness and thus mating success are based may be complex and comprise several traits. Predicting the outcome of sexual selection on such complex phenotypes requires an understanding of both the direct operation of selection on each trait and the indirect consequences of selection operating directly on genetically correlated traits. Here we report the results of a quantitative genetic analysis of the ornamentation, sexual attractiveness, and mating success of male guppies (Poecilia reticulata). We analyze male ornamentation both from the point of view of single ornamental traits (e.g., the area of each color) and of composite measures of the way the entire pattern is likely to be perceived by females (e.g., the mean and contrast in chroma). We demonstrate that there is substantial additive genetic variation in almost all measures of male ornamentation and that much of this variation may be Y linked. Attractiveness and mating success are positively correlated at the phenotypic and genetic level. Orange area and chroma, the area of a male's tail, and the color contrast of his pattern overall are positively correlated with attractiveness and/or mating success at the phenotypic and genetic levels. Using attractiveness and mating success as measures of fitness, we estimate gradients of linear directional sexual selection operating on each male trait and use equations of multivariate evolutionary change to predict the response of male ornamentation to this sexual selection. From these analyses, we predict that indirect selection may have important effects on the evolution of male guppy color patterns.     

The maintenance of heritable variation through social competition

The paradoxical persistence of heritable variation for fitness-related traits is an evolutionary conundrum that remains a preeminent problem in evolutionary biology. Here we describe a simple mechanism in which social competition results in the evolutionary maintenance of heritable variation for fitness related traits. We demonstrate this mechanism using a genetic model with two primary assumptions: the expression of a trait depends upon success in social competition for limited resources; and competitive success of a genotype depends on the genotypes that it competes against. We find that such social competition generates heritable (additive) genetic variation for "competition-dependent" traits. This heritable variation is not eroded by continuous directional selection because, rather than leading to fixation of favored alleles, selection leads instead to allele frequency cycling due to the concerted coevolution of the social environment with the effects of alleles. Our results provide a mechanism for the maintenance of heritable variation in natural populations and suggest an area for research into the importance of competition in the genetic architecture of fitness related traits.     

SELECTION AND EVOLUTION OF CAUSALLY COVARYING TRAITS

When traits cause variation in fitness, the distribution of phenotype, weighted by fitness, necessarily changes. The degree to which traits cause fitness variation is therefore of central importance to evolutionary biology. Multivariate selection gradients are the main quantity used to describe components of trait‐fitness covariation, but they quantify the direct effects of traits on (relative) fitness, which are not necessarily the total effects of traits on fitness. Despite considerable use in evolutionary ecology, path analytic characterizations of the total effects of traits on fitness have not been formally incorporated into quantitative genetic theory. By formally defining “extended” selection gradients, which are the total effects of traits on fitness, as opposed to the existing definition of selection gradients, a more intuitive scheme for characterizing selection is obtained. Extended selection gradients are distinct quantities, differing from the standard definition of selection gradients not only in the statistical means by which they may be assessed and the assumptions required for their estimation from observational data, but also in their fundamental biological meaning. Like direct selection gradients, extended selection gradients can be combined with genetic inference of multivariate phenotypic variation to provide quantitative prediction of microevolutionary trajectories.     

Natural selection, evolvability and bias due to environmental covariance in the field in an annual plant

Estimates of the form and magnitude of natural selection based on phenotypic relationships between traits and fitness measures can be biased when environmental factors influence both relative fitness and phenotypic trait values. I quantified genetic variances and covariances, and estimated linear and quadratic selection coefficients, for seven traits of an annual plant grown in the field. For replicates of 50 paternal half-sib families, coefficients of selection were calculated both for individual phenotypic values of the traits and for half-sib family mean values. The potential for evolutionary response was supported by significant heritability and phenotypic directional selection for several traits but contradicted by the absence of significant genetic variation for fitness estimates and evidence of bias in phenotypic selection coefficients due to environmental covariance for at least two of the traits analysed. Only studies of a much wider range of organisms and traits will reveal the frequency and extent of such bias.     

Selection and evolutionary potential of spring arrival phenology in males and females of a migratory songbird

The timing of annual life‐history events affects survival and reproduction of all organisms. A changing environment can perturb phenological adaptations and an important question is if populations can evolve fast enough to track the environmental changes. Yet, little is known about selection and evolutionary potential of traits determining the timing of crucial annual events. Migratory species, which travel between different climatic regions, are particularly affected by global environmental changes. To increase our understanding of evolutionary potential and selection of timing traits, we investigated the quantitative genetics of arrival date at the breeding ground using a multigenerational pedigree of a natural great reed warbler (Acrocephalus arundinaceus) population. We found significant heritability of 16.4% for arrival date and directional selection for earlier arrival in both sexes acting through reproductive success, but not through lifespan. Mean arrival date advanced with 6 days over 20 years, which is in exact accordance with our predicted evolutionary response based on the breeder's equation. However, this phenotypic change is unlikely to be caused by microevolution, because selection seems mainly to act on the nongenetic component of the trait. Furthermore, demographical changes could also not account for the advancing arrival date. Instead, a strong correlation between spring temperatures and population mean arrival date suggests that phenotypic plasticity best explains the advancement of arrival date in our study population. Our study dissects the evolutionary and environmental forces that shape timing traits and thereby increases knowledge of how populations cope with rapidly changing environments.     

Cultural transmission and the evolution of human behaviour: a general approach based on the Price equation

Transmitted culture can be viewed as an inheritance system somewhat independent of genes that is subject to processes of descent with modification in its own right. Although many authors have conceptualized cultural change as a Darwinian process, there is no generally agreed formal framework for defining key concepts such as natural selection, fitness, relatedness and altruism for the cultural case. Here, we present and explore such a framework using the Price equation. Assuming an isolated, independently measurable culturally transmitted trait, we show that cultural natural selection maximizes cultural fitness, a distinct quantity from genetic fitness, and also that cultural relatedness and cultural altruism are not reducible to or necessarily related to their genetic counterparts. We show that antagonistic coevolution will occur between genes and culture whenever cultural fitness is not perfectly aligned with genetic fitness, as genetic selection will shape psychological mechanisms to avoid susceptibility to cultural traits that bear a genetic fitness cost. We discuss the difficulties with conceptualizing cultural change using the framework of evolutionary theory, the degree to which cultural evolution is autonomous from genetic evolution, and the extent to which cultural change should be seen as a Darwinian process. We argue that the nonselection components of evolutionary change are much more important for culture than for genes, and that this and other important differences from the genetic case mean that different approaches and emphases are needed for cultural than genetic processes.     

Two sexes,one genome: the evolutionary dynamics of intralocus sexual conflict

As the evolutionary interests of males and females are frequently divergent, a trait value that is optimal for the fitness of one sex is often not optimal for the other. A shared genome also means that the same genes may underlie the same trait in both sexes. This can give rise to a form of sexual antagonism, known as intralocus sexual conflict (IASC). Here, a tug‐of‐war over allelic expression can occur, preventing the sexes from reaching optimal trait values, thereby causing sex‐specific reductions in fitness. For some traits, it appears that IASC can be resolved via sex‐specific regulation of genes that subsequently permits sexual dimorphism; however, it seems that whole‐genome resolution may be impossible, due to the genetic architecture of certain traits, and possibly due to the changing dynamics of selection. In this review, we explore the evolutionary mechanisms of, and barriers to, IASC resolution. We also address the broader consequences of this evolutionary feud, the possible interactions between intra‐ and interlocus sexual conflict (IRSC: a form of sexual antagonism involving different loci in each sex), and draw attention to issues that arise from using proxies as measurements of conflict. In particular, it is clear that the sex‐specific fitness consequences of sexual dimorphism require characterization before making assumptions concerning how this relates to IASC. Although empirical data have shown consistent evidence of the fitness effects of IASC, it is essential that we identify the alleles mediating these effects in order to show IASC in its true sense, which is a “conflict over shared genes.”     

Variation,selection and evolution of function-valued traits

We describe an emerging framework for understanding variation, selection and evolution of phenotypic traits that are mathematical functions. We use one specific empirical example – thermal performance curves (TPCs) for growth rates of caterpillars – to demonstrate how models for function-valued traits are natural extensions of more familiar, multivariate models for correlated, quantitative traits. We emphasize three main points. First, because function-valued traits are continuous functions, there are important constraints on their patterns of variation that are not captured by multivariate models. Phenotypic and genetic variation in function-valued traits can be quantified in terms of variance-covariance functions and their associated eigenfunctions: we illustrate how these are estimated as well as their biological interpretations for TPCs. Second, selection on a function-valued trait is itself a function, defined in terms of selection gradient functions. For TPCs, the selection gradient describes how the relationship between an organism's performance and its fitness varies as a function of its temperature. We show how the form of the selection gradient function for TPCs relates to the frequency distribution of environmental states (caterpillar temperatures) during selection. Third, we can predict evolutionary responses of function-valued traits in terms of the genetic variance-covariance and the selection gradient functions. We illustrate how non-linear evolutionary responses of TPCs may occur even when the mean phenotype and the selection gradient are themselves linear functions of temperature. Finally, we discuss some of the methodological and empirical challenges for future studies of the evolution of function-valued traits.     

Quantitative genetic versions of Hamilton's rule with empirical applications

Hamilton''s theory of inclusive fitness revolutionized our understanding of the evolution of social interactions. Surprisingly, an incorporation of Hamilton''s perspective into the quantitative genetic theory of phenotypic evolution has been slow, despite the popularity of quantitative genetics in evolutionary studies. Here, we discuss several versions of Hamilton''s rule for social evolution from a quantitative genetic perspective, emphasizing its utility in empirical applications. Although evolutionary quantitative genetics offers methods to measure each of the critical parameters of Hamilton''s rule, empirical work has lagged behind theory. In particular, we lack studies of selection on altruistic traits in the wild. Fitness costs and benefits of altruism can be estimated using a simple extension of phenotypic selection analysis that incorporates the traits of social interactants. We also discuss the importance of considering the genetic influence of the social environment, or indirect genetic effects (IGEs), in the context of Hamilton''s rule. Research in social evolution has generated an extensive body of empirical work focusing—with good reason—almost solely on relatedness. We argue that quantifying the roles of social and non-social components of selection and IGEs, in addition to relatedness, is now timely and should provide unique additional insights into social evolution.