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.     

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.     

Stabilizing selection and adaptive evolution in a combination of two traits in an arctic ungulate

Stabilizing selection is thought to be common in wild populations and act as one of the main evolutionary mechanisms, which constrain phenotypic variation. When multiple traits interact to create a combined phenotype, correlational selection may be an important process driving adaptive evolution. Here, we report on phenotypic selection and evolutionary changes in two natal traits in a semidomestic population of reindeer (Rangifer tarandus) in northern Finland. The population has been closely monitored since 1969, and detailed data have been collected on individuals since they were born. Over the length of the study period (1969–2015), we found directional and stabilizing selection toward a combination of earlier birth date and heavier birth mass with an intermediate optimum along the major axis of the selection surface. In addition, we demonstrate significant changes in mean traits toward earlier birth date and heavier birth mass, with corresponding genetic changes in breeding values during the study period. Our results demonstrate evolutionary changes in a combination of two traits, which agree closely with estimated patterns of phenotypic selection. Knowledge of the selective surface for combinations of genetically correlated traits are vital to predict how population mean phenotypes and fitness are affected when environments change.     

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.     

Quantification and decomposition of environment‐selection relationships

In nature, selection varies across time in most environments, but we lack an understanding of how specific ecological changes drive this variation. Ecological factors can alter phenotypic selection coefficients through changes in trait distributions or individual mean fitness, even when the trait‐absolute fitness relationship remains constant. We apply and extend a regression‐based approach in a population of Soay sheep (Ovis aries) and suggest metrics of environment‐selection relationships that can be compared across studies. We then introduce a novel method that constructs an environmentally structured fitness function. This allows calculation of full (as in existing approaches) and partial (acting separately through the absolute fitness function slope, mean fitness, and phenotype distribution) sensitivities of selection to an ecological variable. Both approaches show positive overall effects of density on viability selection of lamb mass. However, the second approach demonstrates that this relationship is largely driven by effects of density on mean fitness, rather than on the trait‐fitness relationship slope. If such mechanisms of environmental dependence of selection are common, this could have important implications regarding the frequency of fluctuating selection, and how previous selection inferences relate to longer term evolutionary dynamics.     

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.     

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.     

The evolution of genetic canalization under fluctuating selection

Abstract.— If the direction of selection changes from generation to generation, the ability to respond to selection is maladaptive: the response to selection in one generation leads to reduced fitness in the next. Because the response is determined by the amount of genetic variance expressed at the phenotypic level, rapidly fluctuating selection should favor modifier genes that reduce the phenotypic effect of alleles segregating at structural loci underlying the trait. Such reduction in phenotypic expression of genetic variation has been named "genetic canalization." I support this argument with a series of two- and multilocus models with alternating linear selection and Gaussian selection with fluctuating optimum. A canalizing modifier gene affects the fitness of its carriers in three ways: (1) it reduces the phenotypic consequences of genetic response to previous selection; (2) it reduces the genetic response to selection, which is manifested as linkage disequilibrium between the modifier and structural loci; and (3) it reduces the phenotypic variance. The first two effects reduce fitness under directional selection sustained for several generations, but improve fitness when the direction of selection has just been reversed. The net effect tends to favor a canalizing modifier under rapidly fluctuating selection regimes (period of eight generations or less). The third effect improves fitness of the modifier allele if the fitness function is convex and reduces it if the function is concave. Under fluctuating Gaussian selection, the population is more likely to experience the concave portion of the fitness function when selection is stronger. Therefore, only weak to moderately strong fluctuating Gaussian selection favors genetic canalization. This paper considerably broadens the conditions that favor genetic canalization, which so far has only been postulated to evolve under long-term stabilizing selection.     

The evolution of trade-offs under directional and correlational selection

Using quantitative genetic theory, we develop predictions for the evolution of trade-offs in response to directional and correlational selection. We predict that directional selection favoring an increase in one trait in a trade-off will result in change in the intercept but not the slope of the trade-off function, with the mean value of the selected trait increasing and that of the correlated trait decreasing. Natural selection will generally favor an increase in some combination of trait values, which can be represented as directional selection on an index value. Such selection induces both directional and correlational selection on the component traits. Theory predicts that selection on an index value will also change the intercept but not the slope of the trade-off function but because of correlational selection, the direction of change in component traits may be in the same or opposite directions. We test these predictions using artificial selection on the well-established trade-off between fecundity and flight capability in the cricket, Gryllus firmus and compare the empirical results with a priori predictions made using genetic parameters from a separate half-sibling experiment. Our results support the predictions and illustrate the complexity of trade-off evolution when component traits are subject to both directional and correlational selection.     

Domestication and fitness in the wild: A multivariate view

Domesticated species continually escaping and interbreeding with wild relatives impose a migration load on wild populations. As domesticated stocks become increasingly different as a result of artificial and natural selection in captivity, fitness of escapees in the wild is expected to decline, reducing the effective rate of migration into wild populations. Recent theory suggest that this may alleviate and eventually eliminate the resulting migration load. I develop a multivariate model of trait and wild fitness evolution resulting from the joint effects of artificial and natural selection in the captive environment. Initially, the evolutionary trajectory is dominated by the effects of artificial selection causing a fast initial decline in fitness of escapees in the wild. In later phases, through the counteracting effects of correlational multivariate natural selection in captivity, the mean phenotype is pushed in directions of weak stabilizing selection, allowing a sustained response in the trait subject to artificial selection. Provided that there is some alignment between the adaptive landscapes in the wild and in captivity, these phases are associated with slower rates of decline in wild fitness of the domesticated stock, suggesting that detrimental effects on wild populations are likely to remain a concern in the foreseeable future.     

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.     

Apparent directional selection by biased pleiotropic mutation

Pleiotropic effects of deleterious mutations are considered to be among the factors responsible for genetic constraints on evolution by long-term directional selection acting on a quantitative trait. If pleiotropic phenotypic effects are biased in a particular direction, mutations generate apparent directional selection, which refers to the covariance between fitness and the trait owing to a linear association between the number of mutations possessed by individuals and the genotypic values of the trait. The present analysis has shown how the equilibrium mean value of the trait is determined by a balance between directional selection and biased pleiotropic mutations. Assuming that genes act additively both on the trait and on fitness, the total variance-standardized directional selection gradient was decomposed into apparent and true components. Experimental data on mutation bias from the bristle traits of Drosophila and life history traits of Daphnia suggest that apparent selection explains a small but significant fraction of directional selection pressure that is observed in nature; the data suggest that changes induced in a trait by biased pleiotropic mutation (i.e., by apparent directional selection) are easily compensated for by (true) directional selection.     

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.     

Evolutionary aspects and sensitivity studies of some major gene models

Extensive biometrical and statistically oriented studies in segregation and pedigree analyses reflect current efforts to demonstrate major gene factors playing a significant role for a whole hierarchy of multifactorial diseases and related risk factors exhibiting continuous variation. The evolutionary aspects of the changes in gene frequencies of some major gene one locus models admitting a broad range of genotype-phenotype associations and different forms of selection functions are investigated. The flexibility of differences among the genotypic-phenotypic distribution can take account of variable penetrance expressivity, complex multifarious heterogeneous background effects, or partial dominance concepts. The phenotype distribution and selection function are assumed to be time invariant such that the environments with which the population interacts do not depend on either the phenotypes or the genotypes present in the population of any particular generation. Viability selection optimizing or directional acts on the phenotypic level. We consider random mating, and concentrate mostly on evaluating the nature of the equilibrium structure for the cases of “strong” and “weak” selection. For weak stabilizing selection the determinants of superior genotypic fitness in the class of phenotypic symmetric distributions reside in minimizing a combination of the phenotypic variance and the deviation of the phenotypic mean from the optimal phenotype. With equal means of central phenotype values, a canalizing selection effect signifying fitness superiority for the genotype with minimal variance is in force. For strong stabilizing selection the genotype-phenotype density at the optimal value determines the relative genotype fitness value. For directional selection the determinants of the selection realizations depend on a “standardized” deviation of the mean phenotype distributional value relative to its total variance. The effects of symmetry as against asymmetry in the genotype distributions with prescribed means and variances were investigated by numerical computations.     

GENE FLOW FROM DOMESTICATED SPECIES TO WILD RELATIVES: MIGRATION LOAD IN A MODEL OF MULTIVARIATE SELECTION

Domesticated species frequently spread their genes into populations of wild relatives through interbreeding. The domestication process often involves artificial selection for economically desirable traits. This can lead to an indirect response in unknown correlated traits and a reduction in fitness of domesticated individuals in the wild. Previous models for the effect of gene flow from domesticated species to wild relatives have assumed that evolution occurs in one dimension. Here, I develop a quantitative genetic model for the balance between migration and multivariate stabilizing selection. Different forms of correlational selection consistent with a given observed ratio between average fitness of domesticated and wild individuals offsets the phenotypic means at migration–selection balance away from predictions based on simpler one-dimensional models. For almost all parameter values, correlational selection leads to a reduction in the migration load. For ridge selection, this reduction arises because the distance the immigrants deviates from the local optimum in effect is reduced. For realistic parameter values, however, the effect of correlational selection on the load is small, suggesting that simpler one-dimensional models may still be adequate in terms of predicting mean population fitness and viability.     

Dominance Genetic Variance for Traits Under Directional Selection in Drosophila serrata

In contrast to our growing understanding of patterns of additive genetic variance in single- and multi-trait combinations, the relative contribution of nonadditive genetic variance, particularly dominance variance, to multivariate phenotypes is largely unknown. While mechanisms for the evolution of dominance genetic variance have been, and to some degree remain, subject to debate, the pervasiveness of dominance is widely recognized and may play a key role in several evolutionary processes. Theoretical and empirical evidence suggests that the contribution of dominance variance to phenotypic variance may increase with the correlation between a trait and fitness; however, direct tests of this hypothesis are few. Using a multigenerational breeding design in an unmanipulated population of Drosophila serrata, we estimated additive and dominance genetic covariance matrices for multivariate wing-shape phenotypes, together with a comprehensive measure of fitness, to determine whether there is an association between directional selection and dominance variance. Fitness, a trait unequivocally under directional selection, had no detectable additive genetic variance, but significant dominance genetic variance contributing 32% of the phenotypic variance. For single and multivariate morphological traits, however, no relationship was observed between trait–fitness correlations and dominance variance. A similar proportion of additive and dominance variance was found to contribute to phenotypic variance for single traits, and double the amount of additive compared to dominance variance was found for the multivariate trait combination under directional selection. These data suggest that for many fitness components a positive association between directional selection and dominance genetic variance may not be expected.     

Strong artificial selection in the wild results in predicted small evolutionary change

Estimates of genetic variation and selection allow for quantitative predictions of evolutionary change, at least in controlled laboratory experiments. Natural populations are, however, different in many ways, and natural selection on heritable traits does not always result in phenotypic change. To test whether we were able to predict the evolutionary dynamics of a complex trait measured in a natural, heterogeneous environment, we performed, over an 8-year period, a two-way selection experiment on clutch size in a subdivided island population of great tits (Parus major). Despite strong artificial selection, there was no clear evidence for evolutionary change at the phenotypic level. Environmentally induced differences in clutch size among years are, however, large and can mask evolutionary changes. Indeed, genetic changes in clutch size, inferred from a statistical model, did not deviate systematically from those predicted. Although this shows that estimates of genetic variation and selection can indeed provide quantitative predictions of evolutionary change, also in the wild, it also emphasizes that demonstrating evolution in wild populations is difficult, and that the interpretation of phenotypic trends requires great care.     

Protected polymorphism and evolutionary stability in pleiotropic models with trait-specific dominance

When alleles have pleiotropic effects on a number of quantitative traits, the degree of dominance between a pair of alleles can be different for each trait. Such trait-specific dominance has been studied previously in models for the maintenance of genetic variation by antagonistic effects of an allele on two fitness components. By generalizing these models to an arbitrary number of fitness components or other phenotypic traits with different degrees of dominance, I show that genetic polymorphism is generally impossible without antagonistic fitness effects of different traits and without trait-specific dominance. I also investigate dominance and pleiotropy from a more long-term evolutionary perspective, allowing for the study of general ecological scenarios, and I discuss the effects of trait-specific dominance on evolutionary stability criteria. When selection is mainly directional and only trait-specific dominance and antagonism cause the emergence of polymorphism, then these polymorphisms can be overtaken by single mutants again, such that they are probably short-lived on an evolutionary time scale. Near evolutionarily singular points where directional selection is absent, trait-specific dominance and overdominance facilitate the emergence of polymorphism and cause evolutionary divergence in some cases. An important outcome of these models is that trait-specific dominance allows for the emergence of genetic polymorphisms without a selective disadvantage for heterozygotes. This removes the scope for the evolution of assortative mate choice and affects dominance modification. Sympatric speciation by disruptive ecological selection requires this heterozygote disadvantage in order to evolve, and therefore it becomes less plausible if the emergence of genetic polymorphism usually occurs via trait-specific dominance and antagonistic effects.     

The effect of trait type and strength of selection on heritability and evolvability in an island bird population

The heritability (h²) of fitness traits is often low. Although this has been attributed to directional selection having eroded genetic variation in direct proportion to the strength of selection, heritability does not necessarily reflect a trait's additive genetic variance and evolutionary potential (“evolvability”). Recent studies suggest that the low h² of fitness traits in wild populations is caused not by a paucity of additive genetic variance (V_A) but by greater environmental or nonadditive genetic variance (V_R). We examined the relationship between h² and variance‐standardized selection intensities (i or β_σ), and between evolvability (I_A:V_A divided by squared phenotypic trait mean) and mean‐standardized selection gradients (β_μ). Using 24 years of data from an island population of Savannah sparrows, we show that, across diverse traits, h² declines with the strength of selection, whereas I_A and I_R (V_R divided by squared trait mean) are independent of the strength of selection. Within trait types (morphological, reproductive, life‐history), h², I_A, and I_R are all independent of the strength of selection. This indicates that certain traits have low heritability because of increased residual variance due to the age at which they are expressed or the multiple factors influencing their expression, rather than their association with fitness.     

Is photomorphogenic shade avoidance adaptive? Perspectives from population biology

Photomorphogenic shade avoidance responses provide an ideal model system for integrating genetic, physiological and population biology approaches to the study of adaptive plasticity. The adaptive plasticity hypothesis predicts that shade avoidance phenotypes induced by low ratios of red to far-red light (R:FR) will have high relative fitness in dense stands, but will suffer a fitness disadvantage at low density. Experiments with transgenic and mutant plants in which photomorphogenic genes are disabled, as well as phenotype manipulation by means of altered R:FR, strongly support the shade avoidance hypothesis. The observation of photomorphogenic ecotypes in different selective environments also suggests that the shade avoidance response has undergone adaptive evolution. Quantitative genetic variation in R:FR sensitivity has been detected in wild populations, indicating that the evolutionary potential exists for response to natural selection. However, evolutionary response may be constrained by genetic correlations among developmentally linked traits. Therefore it cannot be assumed that an observed suite of photomorphogenic responses represents an adaptive optimum for every trait.