Genetic basis of between‐individual and within‐individual variance of docility

Between‐individual variation in phenotypes within a population is the basis of evolution. However, evolutionary and behavioural ecologists have mainly focused on estimating between‐individual variance in mean trait and neglected variation in within‐individual variance, or predictability of a trait. In fact, an important assumption of mixed‐effects models used to estimate between‐individual variance in mean traits is that within‐individual residual variance (predictability) is identical across individuals. Individual heterogeneity in the predictability of behaviours is a potentially important effect but rarely estimated and accounted for. We used 11 389 measures of docility behaviour from 1576 yellow‐bellied marmots (Marmota flaviventris) to estimate between‐individual variation in both mean docility and its predictability. We then implemented a double hierarchical animal model to decompose the variances of both mean trait and predictability into their environmental and genetic components. We found that individuals differed both in their docility and in their predictability of docility with a negative phenotypic covariance. We also found significant genetic variance for both mean docility and its predictability but no genetic covariance between the two. This analysis is one of the first to estimate the genetic basis of both mean trait and within‐individual variance in a wild population. Our results indicate that equal within‐individual variance should not be assumed. We demonstrate the evolutionary importance of the variation in the predictability of docility and illustrate potential bias in models ignoring variation in predictability. We conclude that the variability in the predictability of a trait should not be ignored, and present a coherent approach for its quantification.     

Inclusive heritability: combining genetic and non‐genetic information to study animal behavior and culture

Phenotypic variance results from variation in biological information possessed by individuals. Quantitative geneticists often strive to partition out all environmental variance to measure heritability. Behavioral biologists and ecologists however, require methods to integrate genetic and environmental components of inherited phenotypic variance in order to estimate the evolutionary potential of traits, which encompasses any form of information that is inherited. To help develop this integration, we build on the tools of quantitative genetics and offer the concept of ‘inclusive heritability’ which identifies and unifies the various mechanisms of information transmission across generations. A controversial component of non‐genetic information is animal culture, which is the part of phenotypic variance inherited through social learning. Culture has the unique property of being transmitted horizontally and obliquely, as well as vertically. Accounting for cultural variation would allow us to examine a broader range of evolutionary mechanisms. Culture may, for instance, produce behavioral isolating mechanisms leading to speciation. To advance the study of animal culture, we offer a definition of culture that is rooted in quantitative genetics. We also offer four testable criteria to determine whether a trait is culturally inherited. These criteria may constitute a conceptual tool to study animal culture. We briefly discuss methods to partition out cultural variance. Several authors have recently called for ‘modernizing the modern synthesis’ by including non‐genetic factors such as epigenetics and phenotypic plasticity in order to more fully explain phenotypic evolution. Here, we further propose to broaden the concept of inheritance by incorporating the cultural component of behavior. Applying the concept of inclusive heritability may advance the integration of multiple forms of inheritance into the study of evolution.     

Plasticity and evolution in correlated suites of traits

When organisms are faced with new or changing environments, a central challenge is the coordination of adaptive shifts in many different phenotypic traits. Relationships among traits may facilitate or constrain evolutionary responses to selection, depending on whether the direction of selection is aligned or opposed to the pattern of trait correlations. Attempts to predict evolutionary potential in correlated traits generally assume that correlations are stable across time and space; however, increasing evidence suggests that this may not be the case, and flexibility in trait correlations could bias evolutionary trajectories. We examined genetic and environmental influences on variation and covariation in a suite of behavioural traits to understand if and how flexibility in trait correlations influences adaptation to novel environments. We tested the role of genetic and environmental influences on behavioural trait correlations by comparing Trinidadian guppies (Poecilia reticulata) historically adapted to high‐ and low‐predation environments that were reared under native and non‐native environmental conditions. Both high‐ and low‐predation fish exhibited increased behavioural variance when reared under non‐native vs. native environmental conditions, and rearing in the non‐native environment shifted the major axis of variation among behaviours. Our findings emphasize that trait correlations observed in one population or environment may not predict correlations in another and that environmentally induced plasticity in correlations may bias evolutionary divergence in novel environments.     

A test of Darwin's ‘lop‐eared’ rabbit hypothesis

Integration of evolutionary and developmental biology has stimulated novel insights on the origins and maintenance of phenotypic variation. For instance, phenotypic accommodation predicts that trait covariance originates via a novel developmental input caused by genetic change in one trait, but not the other. Darwin provided a striking example of this process in the ‘lop‐eared’ rabbit by demonstrating that artificial selection for long external ears induced variation in the external auditory meatus. Although this intriguing pattern has been interpreted as evidence of phenotypic accommodation, it is unclear whether it exists and, if it does, whether it is selectively maintained in nature. To address this concern, we examined trait covariance in natural woodrat populations that have likely undergone selection for long ears. We demonstrated a remarkably similar covariance pattern as in the ‘lop‐eared’ rabbit, which was associated with climatic variables along a steep arid‐to‐moist longitudinal gradient. Thus, our results suggest that trait covariance is likely a correlated response to selection. We relate these findings to potential origins of trait covariance owing to altered developmental interactions, such as in phenotypic accommodation. Additional evidence is needed to clarify how phenotypic accommodation and correlated selection promote and maintain trait covariance in natural populations. Nonetheless, our study is the first to support a classic Darwinian example concerning domestication and natural selection.     

Quantifying the genetic component of phenotypic variation in unpedigreed wild plants: tailoring genomic scan for within-population use

The study of adaptive genetic variation in natural populations is central to evolutionary biology. Quantitative genetics methods, however, are hardly applicable to long-lived organisms, and current knowledge on adaptive genetic variation in wild plants mostly refers to annuals and short-lived perennials. Studies on long-lived species are essential to explore possible life-history correlates of genetic variation, selection, and trait heritability. In this paper, we propose a method based on molecular markers to quantify the genetic basis of individual phenotypic differences in wild plants under natural conditions. Rather than focusing on inferring individual relatedness to estimate the heritability of phenotypic traits, we directly estimate the proportion of observed phenotypic variance that is statistically accounted for by genotypic differences between individuals. This is achieved by (i) identifying loci that are correlated across individuals with the phenotypic trait of interest by means of an amplified fragment length polymorphism (AFLP)-based explorative genomic scan, and (ii) fitting multiple regression and linear random effect models to estimate the effects of genotype, environment and genotype × environment on phenotypes. We apply this method to estimate genotypic and environmental effects on cumulative maternal fecundity in a wild population of the long-lived Viola cazorlensis monitored for 20 years. Results show that between 56–63% (depending on estimation method) of phenotypic variance in fecundity is accounted for by genotypic differences in 11 AFLP loci that are significantly related to fecundity. Genotype × environment effects accounted for 38% of fecundity variance, which may help to explain the unexpectedly high levels of genetic variance for fecundity found.     

The correlation between coloration and exploration behaviour varies across hierarchical levels in a wild passerine bird

In vertebrates, darker individuals are often found to be more active and willing to take risks (representing characteristics of a ‘proactive’ coping style), whereas lighter individuals are instead more cautious and less active (representing characteristics of a ‘reactive’ coping style). It is thus generally expected that melanin‐based coloration and proactivity form a suite of positively integrated traits at the among‐individual level. Here, we use a multigenerational pedigree of free‐living great tits (Parus major) to partition variation in, and the correlation between, melanin‐based breast stripe (‘tie’) size and exploration behaviour (a proxy for coping style) into its among‐ and within‐individual components. We show that both traits harbour heritable variation. Against predictions, tie size and speed of exploration were negatively correlated at the among‐individual level due to the combined influences of permanent environmental and additive genetic effects. By contrast, the two traits were weakly positively correlated within individuals (i.e. individuals increasing in tie size after moult tended to become more explorative). The patterns of among‐individual covariance were not caused by correlational selection as we found additive and opposite selection pressures acting on the two traits. These findings imply that testing hypotheses regarding the existence of a ‘syndrome’ at the among‐individual level strictly requires variance partitioning to avoid inappropriate interpretations as the negative ‘unpartitioned’ phenotypic correlation between exploration and tie size resulted from counteracting effects of within‐ and among‐individual correlations. Identifying sources and levels of (co)variation in phenotypic traits is thus critical to our understanding of biological patterns and evolutionary processes.     

Phenotypic and Evolutionary Consequences of Social Behaviours: Interactions among Individuals Affect Direct Genetic Effects

Traditional quantitative genetics assumes that an individual''s phenotype is determined by both genetic and environmental factors. For many animals, part of the environment is social and provided by parents and other interacting partners. When expression of genes in social partners affects trait expression in a focal individual, indirect genetic effects occur. In this study, we explore the effects of indirect genetic effects on the magnitude and range of phenotypic values in a focal individual in a multi-member model analyzing three possible classes of interactions between individuals. We show that social interactions may not only cause indirect genetic effects but can also modify direct genetic effects. Furthermore, we demonstrate that both direct and indirect genetic effects substantially alter the range of phenotypic values, particularly when a focal trait can influence its own expression via interactions with traits in other individuals. We derive a function predicting the relative importance of direct versus indirect genetic effects. Our model reveals that both direct and indirect genetic effects can depend to a large extent on both group size and interaction strength, altering group mean phenotype and variance. This may lead to scenarios where between group variation is much higher than within group variation despite similar underlying genetic properties, potentially affecting the level of selection. Our analysis highlights key properties of indirect genetic effects with important consequences for trait evolution, the level of selection and potentially speciation.     

Trees to treehoppers: genetic variation in host plants contributes to variation in the mating signals of a plant‐feeding insect

Community genetics research has demonstrated ‘bottom‐up’ effects of genetic variation within a plant species in shaping the larger community with which it interacts, such as compositions of arthropod faunas. We demonstrate that such cross‐trophic interactions also influence sexually selected traits. We used a member of the Enchenopa binotata species complex of treehoppers (Hemiptera: Membracidae) to ask whether male mating signals are influenced by host plant genetic variation. We reared a random sample of the treehoppers on potted replicates of a sample of host plant clone lines. We found that treehopper male signals varied according to the clone line on which they developed, showing that genetic variation in host plants affects male treehoppers' behavioural phenotypes. This is the first demonstration of cross‐trophic indirect genetic effects on a sexually selected trait. We discuss how such effects may play an important role in the maintenance of variation and within‐population phenotypic differentiation, thereby promoting evolutionary divergence.     

Life‐history evolution under fluctuating density‐dependent selection and the adaptive alignment of pace‐of‐life syndromes

We present a novel perspective on life‐history evolution that combines recent theoretical advances in fluctuating density‐dependent selection with the notion of pace‐of‐life syndromes (POLSs) in behavioural ecology. These ideas posit phenotypic co‐variation in life‐history, physiological, morphological and behavioural traits as a continuum from the highly fecund, short‐lived, bold, aggressive and highly dispersive ‘fast’ types at one end of the POLS to the less fecund, long‐lived, cautious, shy, plastic and socially responsive ‘slow’ types at the other. We propose that such variation in life histories and the associated individual differences in behaviour can be explained through their eco‐evolutionary dynamics with population density – a single and ubiquitous selective factor that is present in all biological systems. Contrasting regimes of environmental stochasticity are expected to affect population density in time and space and create differing patterns of fluctuating density‐dependent selection, which generates variation in fast versus slow life histories within and among populations. We therefore predict that a major axis of phenotypic co‐variation in life‐history, physiological, morphological and behavioural traits (i.e. the POLS) should align with these stochastic fluctuations in the multivariate fitness landscape created by variation in density‐dependent selection. Phenotypic plasticity and/or genetic (co‐)variation oriented along this major POLS axis are thus expected to facilitate rapid and adaptively integrated changes in various aspects of life histories within and among populations and/or species. The fluctuating density‐dependent selection POLS framework presented here therefore provides a series of clear testable predictions, the investigation of which should further our fundamental understanding of life‐history evolution and thus our ability to predict natural population dynamics.     

Genetic and environmental effects on a condition‐dependent trait: feather growth in Siberian jays

Condition, defined as the amount of ‘internal resources’ an individual can freely allocate, is often assumed to be environmentally determined and to reflect an individual’s health and nutritional status. However, an additive genetic component of condition is possible if it ‘captures’ the genetic variance of many underlying traits as many fitness‐related traits appear to do. Yet, the heritability of condition can be low if selection has eroded much of its additive genetic variance, or if the environmental influences are strong. Here, we tested whether feather growth rate – presumably a condition‐dependent trait – has a heritable component, and whether variation in feather growth rate is related to variation in fitness. To this end, we utilized data from a long‐term population study of Siberian jays (Perisoreus infaustus), and found that feather growth rate, measured as the width of feather growth bars (GB), differed between age‐classes and sexes, but was only weakly related to variation in fitness as measured by annual and life‐time reproductive success. As revealed by animal model analyses, GB width was significantly heritable (h² = 0.10 ± 0.05), showing that this measure of condition is not solely environmentally determined, but reflects at least partly inherited genetic differences among individuals. Consequently, variation in feather growth rates as assessed with ptilochronological methods can provide information about heritable genetic differences in condition.     

A general definition of the heritable variation that determines the potential of a population to respond to selection

Genetic selection is a major force shaping life on earth. In classical genetic theory, response to selection is the product of the strength of selection and the additive genetic variance in a trait. The additive genetic variance reflects a population's intrinsic potential to respond to selection. The ordinary additive genetic variance, however, ignores the social organization of life. With social interactions among individuals, individual trait values may depend on genes in others, a phenomenon known as indirect genetic effects. Models accounting for indirect genetic effects, however, lack a general definition of heritable variation. Here I propose a general definition of the heritable variation that determines the potential of a population to respond to selection. This generalizes the concept of heritable variance to any inheritance model and level of organization. The result shows that heritable variance determining potential response to selection is the variance among individuals in the heritable quantity that determines the population mean trait value, rather than the usual additive genetic component of phenotypic variance. It follows, therefore, that heritable variance may exceed phenotypic variance among individuals, which is impossible in classical theory. This work also provides a measure of the utilization of heritable variation for response to selection and integrates two well-known models of maternal genetic effects. The result shows that relatedness between the focal individual and the individuals affecting its fitness is a key determinant of the utilization of heritable variance for response to selection.     

Age-dependent genetic variance in a life-history trait in the mute swan

Genetic variance in characters under natural selection in natural populations determines the way those populations respond to that selection. Whether populations show temporal and/or spatial constancy in patterns of genetic variance and covariance is regularly considered, as this will determine whether selection responses are constant over space and time. Much less often considered is whether characters show differing amounts of genetic variance over the life-history of individuals. Such age-specific variation, if present, has important potential consequences for the force of natural selection and for understanding the causes of variation in quantitative characters. Using data from a long-term study of the mute swan Cygnus olor, we report the partitioning of phenotypic variance in timing of breeding (subject to strong natural selection) into component parts over 12 different age classes. We show that the additive genetic variance and heritability of this trait are strongly age-dependent, with higher additive genetic variance present in young and, particularly, old birds, but little evidence of any genetic variance for birds of intermediate ages. These results demonstrate that age can have a very important influence on the components of variation of characters in natural populations, and consequently that separate age classes cannot be assumed to be equivalent, either with respect to their evolutionary potential or response.     

Social Selection and Indirect Genetic Effects in Structured populations

Social selection and indirect genetic effects (IGEs) are established concepts in both behavioural ecology and evolutionary genetics. While IGEs describe effects of an individual’s genotype on phenotypes of social partners (and may thus affect their fitness indirectly), the concept of social selection assumes that a given phenotype in one individual affects the fitness of other individuals directly. Although different frameworks, both have been used to investigate the evolution of social traits, such as cooperative behaviour. Despite their similarities (both concepts consider interactions among individuals), they differ in the type of interaction. It remains unclear whether the two concepts make the same predictions about evolutionary trajectories or not. To address this question, we investigate four possible scenarios of social interactions and compare the effects of IGEs and social selection for trait evolution in a multi-trait multi-member model. We show that the two mechanisms can yield similar evolutionary outcomes and that both can create selection pressure at the group level. However, the effect of IGEs can be stronger due to the possibility of feedback loops. Finally, we demonstrate that IGEs, but not social selection gradients, may lead to differences in the direction of evolutionary response between genotypes and phenotypes.     

An evolutionary quantitative genetics model for phenotypic (co)variances under limited dispersal,with an application to socially synergistic traits

Darwinian evolution consists of the gradual transformation of heritable traits due to natural selection and the input of random variation by mutation. Here, we use a quantitative genetics approach to investigate the coevolution of multiple quantitative traits under selection, mutation, and limited dispersal. We track the dynamics of trait means and of variance–covariances between traits that experience frequency‐dependent selection. Assuming a multivariate‐normal trait distribution, we recover classical dynamics of quantitative genetics, as well as stability and evolutionary branching conditions of invasion analyses, except that due to limited dispersal, selection depends on indirect fitness effects and relatedness. In particular, correlational selection that associates different traits within‐individuals depends on the fitness effects of such associations between‐individuals. We find that these kin selection effects can be as relevant as pleiotropy for the evolution of correlation between traits. We illustrate this with an example of the coevolution of two social traits whose association within‐individuals is costly but synergistically beneficial between‐individuals. As dispersal becomes limited and relatedness increases, associations between‐traits between‐individuals become increasingly targeted by correlational selection. Consequently, the trait distribution goes from being bimodal with a negative correlation under panmixia to unimodal with a positive correlation under limited dispersal.     

EXPERIMENTAL EVIDENCE FOR THE EVOLUTION OF INDIRECT GENETIC EFFECTS: CHANGES IN THE INTERACTION EFFECT COEFFICIENT,PSI (Ψ), DUE TO SEXUAL SELECTION

Indirect genetics effects (IGEs)—when the genotype of one individual affects the phenotypic expression of a trait in another—may alter evolutionary trajectories beyond that predicted by standard quantitative genetic theory as a consequence of genotypic evolution of the social environment. For IGEs to occur, the trait of interest must respond to one or more indicator traits in interacting conspecifics. In quantitative genetic models of IGEs, these responses (reaction norms) are termed interaction effect coefficients and are represented by the parameter psi (Ψ). The extent to which Ψ exhibits genetic variation within a population, and may therefore itself evolve, is unknown. Using an experimental evolution approach, we provide evidence for a genetic basis to the phenotypic response caused by IGEs on sexual display traits in Drosophila serrata. We show that evolution of the response is affected by sexual but not natural selection when flies adapt to a novel environment. Our results indicate a further mechanism by which IGEs can alter evolutionary trajectories—the evolution of interaction effects themselves.     

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.     

Exploring patterns of variation in clutch size–density reaction norms in a wild passerine bird

Negative density dependence of clutch size is a ubiquitous characteristic of avian populations and is partly due to within‐individual phenotypic plasticity. Yet, very little is known about the extent to which individuals differ in their degree of phenotypic plasticity, whether such variation has a genetic basis and whether level of plasticity can thus evolve in response to selection. Using 18 years of data of a Dutch great tit population (Parus major), we show that females reduced clutch size with increasing population density (slopes of the reaction norms), differed strongly in their average clutch size (elevations of the reaction norms) at the population‐mean density and that the latter variation was partly heritable. In contrast, we could not detect individual variation in phenotypic plasticity (‘I × E’). Level of plasticity is thus not likely to evolve in response to selection in this population. Observed clutch sizes deviated more from the estimated individual reaction norms in certain years and densities, implying that the within‐individual between‐year variance (so‐called residual variance) of clutch size was heterogeneous with respect to these factors. Given the observational nature of this study, experimental manipulation of density is now warranted to confirm the causality of the observed density effects. Our analyses demonstrate that failure to acknowledge this heterogeneity would have inflated the estimate of ‘I × E’ and led to misinterpretation of the data. This paper thereby emphasizes the fact that heterogeneity in residuals can provide biologically insightful information about the ecological processes underlying the data.     

Characterizing the evolution of genetic variance using genetic covariance tensors

Determining how genetic variance changes under selection in natural populations has proved to be a very resilient problem in evolutionary genetics. In the same way that understanding the availability of genetic variance within populations requires the simultaneous consideration of genetic variance in sets of functionally related traits, determining how genetic variance changes under selection in natural populations will require ascertaining how genetic variance–covariance (G) matrices evolve. Here, we develop a geometric framework using higher-order tensors, which enables the empirical characterization of how G matrices have diverged among populations. We then show how divergence among populations in genetic covariance structure can then be associated with divergence in selection acting on those traits using key equations from evolutionary theory. Using estimates of G matrices of eight male sexually selected traits from nine geographical populations of Drosophila serrata, we show that much of the divergence in genetic variance occurred in a single trait combination, a conclusion that could not have been reached by examining variation among the individual elements of the nine G matrices. Divergence in G was primarily in the direction of the major axes of genetic variance within populations, suggesting that genetic drift may be a major cause of divergence in genetic variance among these populations.     

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 role of common ancestry and gene flow in the evolution of human‐directed play behaviour in dogs

Among‐population variance of phenotypic traits is of high relevance for understanding evolutionary mechanisms that operate in relatively short timescales, but various sources of nonindependence, such as common ancestry and gene flow, can hamper the interpretations. In this comparative analysis of 138 dog breeds, we demonstrate how such confounders can independently shape the evolution of a behavioural trait (human‐directed play behaviour from the Dog Mentality Assessment project). We combined information on genetic relatedness and haplotype sharing to reflect common ancestry and gene flow, respectively, and entered these into a phylogenetic mixed model to partition the among‐breed variance of human‐directed play behaviour while also accounting for within‐breed variance. We found that 75% of the among‐breed variance was explained by overall genetic relatedness among breeds, whereas 15% could be attributed to haplotype sharing that arises from gene flow. Therefore, most of the differences in human‐directed play behaviour among breeds have likely been caused by constraints of common ancestry as a likely consequence of past selection regimes. On the other hand, gene flow caused by crosses among breeds has played a minor, but not negligible role. Our study serves as an example of an analytical approach that can be applied to comparative situations where the effects of shared origin and gene flow require quantification and appropriate statistical control in a within‐species/among‐population framework. Altogether, our results suggest that the evolutionary history of dog breeds has left remarkable signatures on the among‐breed variation of a behavioural phenotype.