The past,present and future of reproductive skew theory and experiments

A major evolutionary question is how reproductive sharing arises in cooperatively breeding species despite the inherent reproductive conflicts in social groups. Reproductive skew theory offers one potential solution: each group member gains or is allotted inclusive fitness equal to or exceeding their expectation from reproducing on their own. Unfortunately, a multitude of skew models with conflicting predictions has led to confusion in both testing and evaluating skew theory. The confusion arises partly because one set of models (the ‘transactional’ type) answer the ultimate evolutionary question of what ranges of reproductive skew can yield fitness‐enhancing solutions for all group members. The second set of models (‘compromise’) give an evolutionarily proximate, game‐theoretic evolutionarily stable state (ESS) solution that determines reproductive shares based on relative competitive abilities. However, several predictions arising from compromise models require a linear payoff to increased competition and do not hold with non‐linear payoffs. Given that for most species it may be very difficult or impossible to determine the true relationship between effort devoted to competition and reproductive share gained, compromise models are much less predictive than previously appreciated. Almost all skew models make one quantitative prediction (e.g. realized skew must fall within ranges predicted by transactional models), and two qualitative predictions (e.g. variation in relatedness or competitive ability across groups affects skew). A thorough review of the data finds that these three predictions are relatively rarely supported. As a general rule, therefore, the evolution of cooperative breeding appears not to be dependent on the ability of group members to monitor relatedness or competitive ability in order to adjust their behaviour dynamically to gain reproductive share. Although reproductive skew theory fails to predict within‐group dynamics consistently, it does better at predicting quantitative differences in skew across populations or species. This suggests that kin selection can play a significant role in the evolution of sociality. To advance our understanding of reproductive skew will require focusing on a broader array of factors, such as the frequency of mistaken identity, delayed fitness payoffs, and selection pressures arising from across‐group competition. We furthermore suggest a novel approach to investigate the sharing of reproduction that focuses on the underlying genetics of skew. A quantitative genetics approach allows the partitioning of variance in reproductive share itself or that of traits closely associated with skew into genetic and non‐genetic sources. Thus, we can determine the heritability of reproductive share and infer whether it actually is the focus of natural selection. We view the ‘animal model’ as the most promising empirical method where the genetics of reproductive share can be directly analyzed in wild populations. In the quest to assess whether skew theory can provide a framework for understanding the evolution of sociality, quantitative genetics will be a central tool in future research.     

Evolution and stability of the G-matrix on a landscape with a moving optimum

In quantitative genetics, the genetic architecture of traits, described in terms of variances and covariances, plays a major role in determining the trajectory of evolutionary change. Hence, the genetic variance-covariance matrix (G-matrix) is a critical component of modern quantitative genetics theory. Considerable debate has surrounded the issue of G-matrix constancy because unstable G-matrices provide major difficulties for evolutionary inference. Empirical studies and analytical theory have not resolved the debate. Here we present the results of stochastic models of G-matrix evolution in a population responding to an adaptive landscape with an optimum that moves at a constant rate. This study builds on the previous results of stochastic simulations of G-matrix stability under stabilizing selection arising from a stationary optimum. The addition of a moving optimum leads to several important new insights. First, evolution along genetic lines of least resistance increases stability of the orientation of the G-matrix relative to stabilizing selection alone. Evolution across genetic lines of least resistance decreases G-matrix stability. Second, evolution in response to a continuously changing optimum can produce persistent maladaptation for a correlated trait, even if its optimum does not change. Third, the retrospective analysis of selection performs very well when the mean G-matrix (G) is known with certainty, indicating that covariance between G and the directional selection gradient beta is usually small enough in magnitude that it introduces only a small bias in estimates of the net selection gradient. Our results also show, however, that the contemporary G-matrix only serves as a rough guide to G. The most promising approach for the estimation of G is probably through comparative phylogenetic analysis. Overall, our results show that directional selection actually can increase stability of the G-matrix and that retrospective analysis of selection is inherently feasible. One major remaining challenge is to gain a sufficient understanding of the G-matrix to allow the confident estimation of G.     

数量性状发育遗传模型及其分析方法的研究进展

发育遗传模型是同时反映性状遗传和发育本质、提供影响遗传变异及调整发育进程的有关因素的信息的模型。建立在群体遗传学基础上的直接效应模型适用于单一基因控制的简单性状。渐成模型将遗传变异分解成直接分量和渐成分量(母体效应和互作效应),能更好地反映有机体遗传和发育的生物学机制。生长轨迹模型有效地综合了复杂性状各分量的发育动态,可获得连续的、综合的、详细的、动态的发育信息。条件遗传分析方法不仅可以估算特定时间段的净效应,且可将净效应分解为不同遗传分量,了解各效应分量的相对贡献。 Abstract：Developmental genetic models and analysis methods for quantitative traits are presented.Developmental genetic models should reflect the genetic and developmental essence,and provide the information of the factors influencing the genetic variation and the developmental process.Direct effect models,which based on the population genetics,may be suitable to analyze simple traits with single gene.Epigenetic models can decompose the whole genetic variation into direct and epigenetic components (maternal effects and epigenetic interaction effects),so that biological mechanism can be better understood.Growth trace models effectively synthesize the developmental dynamics of components of complex traits.With them,continuous,compositive,detailed,and dynamic information of development is available.Conditional analysis method can not only estimate the net effects in a specific time interval,but also depose them into genetic components and help to appreciate the contributions of different effects.     

Genetic variability under mutation selection balance

A fundamental problem in evolutionary genetics is understanding how high levels of genetic variation in quantitative traits are maintained in natural populations. Variation is removed by the natural selection of individuals with optimal phenotypes and is recovered by mutation; however, previous analyses had indicated that a mutation-selection balance was insufficient to maintain observed levels of genetic variation in these traits. Using more general models, however, it has recently been shown that it is indeed a sufficient mechanism. These models can be used to explore other phenomena in evolutionary biology.     

Genomic Heritability: What Is It?

Whole-genome regression methods are being increasingly used for the analysis and prediction of complex traits and diseases. In human genetics, these methods are commonly used for inferences about genetic parameters, such as the amount of genetic variance among individuals or the proportion of phenotypic variance that can be explained by regression on molecular markers. This is so even though some of the assumptions commonly adopted for data analysis are at odds with important quantitative genetic concepts. In this article we develop theory that leads to a precise definition of parameters arising in high dimensional genomic regressions; we focus on the so-called genomic heritability: the proportion of variance of a trait that can be explained (in the population) by a linear regression on a set of markers. We propose a definition of this parameter that is framed within the classical quantitative genetics theory and show that the genomic heritability and the trait heritability parameters are equal only when all causal variants are typed. Further, we discuss how the genomic variance and genomic heritability, defined as quantitative genetic parameters, relate to parameters of statistical models commonly used for inferences, and indicate potential inferential problems that are assessed further using simulations. When a large proportion of the markers used in the analysis are in LE with QTL the likelihood function can be misspecified. This can induce a sizable finite-sample bias and, possibly, lack of consistency of likelihood (or Bayesian) estimates. This situation can be encountered if the individuals in the sample are distantly related and linkage disequilibrium spans over short regions. This bias does not negate the use of whole-genome regression models as predictive machines; however, our results indicate that caution is needed when using marker-based regressions for inferences about population parameters such as the genomic heritability.     

The ecology and evolution of host-plant resistance to insects

Genetic techniques have yielded new insights into plant-herbivore coevolution. Quantitative genetic tests of herbivory theory reveal that in some cases insect herbivores impose selection on resistance traits. Also, some resistance traits are costly while others appear not to be, and genetic models can explain these results. Genetic variation in plant resistance influences insect community structure by modifying interactions of herbivores with competitors and natural enemies. Therefore, models of multispecies coevolution are more realistic than pairwise coevolutionary models. Ecological genetics will facilitate further theoretical and empirical exploration of multispecies coevolution of plants and herbivores.     

Stability of the G-matrix in a population experiencing pleiotropic mutation, stabilizing selection, and genetic drift

Abstract. Quantitative genetics theory provides a framework that predicts the effects of selection on a phenotype consisting of a suite of complex traits. However, the ability of existing theory to reconstruct the history of selection or to predict the future trajectory of evolution depends upon the evolutionary dynamics of the genetic variance-covariance matrix (G-matrix). Thus, the central focus of the emerging field of comparative quantitative genetics is the evolution of the G-matrix. Existing analytical theory reveals little about the dynamics of G, because the problem is too complex to be mathematically tractable. As a first step toward a predictive theory of G-matrix evolution, our goal was to use stochastic computer models to investigate factors that might contribute to the stability of G over evolutionary time. We were concerned with the relatively simple case of two quantitative traits in a population experiencing stabilizing selection, pleiotropic mutation, and random genetic drift. Our results show that G-matrix stability is enhanced by strong correlational selection and large effective population size. In addition, the nature of mutations at pleiotropic loci can dramatically influence stability of G. In particular, when a mutation at a single locus simultaneously changes the value of the two traits (due to pleiotropy) and these effects are correlated, mutation can generate extreme stability of G. Thus, the central message of our study is that the empirical question regarding G-matrix stability is not necessarily a general question of whether G is stable across various taxonomic levels. Rather, we should expect the G-matrix to be extremely stable for some suites of characters and unstable for others over similar spans of evolutionary time.     

Converging phenomics and genomics to study natural variation in plant photosynthetic efficiency

In recent years developments in plant phenomic approaches and facilities have gradually caught up with genomic approaches. An opportunity lies ahead to dissect complex, quantitative traits when both genotype and phenotype can be assessed at a high level of detail. This is especially true for the study of natural variation in photosynthetic efficiency, for which forward genetics studies have yielded only a little progress in our understanding of the genetic layout of the trait. High‐throughput phenotyping, primarily from chlorophyll fluorescence imaging, should help to dissect the genetics of photosynthesis at the different levels of both plant physiology and development. Specific emphasis should be directed towards understanding the acclimation of the photosynthetic machinery in fluctuating environments, which may be crucial for the identification of genetic variation for relevant traits in food crops. Facilities should preferably be designed to accommodate phenotyping of photosynthesis‐related traits in such environments. The use of forward genetics to study the genetic architecture of photosynthesis is likely to lead to the discovery of novel traits and/or genes that may be targeted in breeding or bio‐engineering approaches to improve crop photosynthetic efficiency. In the near future, big data approaches will play a pivotal role in data processing and streamlining the phenotype‐to‐gene identification pipeline.     

The influence of nonrandom extra‐pair paternity on heritability estimates derived from wild pedigrees

Quantitative genetic analysis is often fundamental for understanding evolutionary processes in wild populations. Avian populations provide a model system due to the relative ease of inferring relatedness among individuals through observation. However, extra‐pair paternity (EPP) creates erroneous links within the social pedigree. Previous work has suggested this causes minor underestimation of heritability if paternal misassignment is random and hence not influenced by the trait being studied. Nevertheless, much literature suggests numerous traits are associated with EPP and the accuracy of heritability estimates for such traits remains unexplored. We show analytically how nonrandom pedigree errors can influence heritability estimates. Then, combining empirical data from a large great tit (Parus major) pedigree with simulations, we assess how heritability estimates derived from social pedigrees change depending on the mode of the relationship between EPP and the focal trait. We show that the magnitude of the underestimation is typically small (<15%). Hence, our analyses suggest that quantitative genetic inference from pedigrees derived from observations of social relationships is relatively robust; our approach also provides a widely applicable method for assessing the consequences of nonrandom EPP.     

Estimation and testing of parent-of-origin effects for quantitative traits

Recent progress in developing family-based association methods has extended their use to the analysis of quantitative traits in the offspring and to the estimation, for dichotomous traits, of the relative contribution of genetic and environmental mechanisms for parent-of-origin effects. However, many traits of interest are not naturally measured on a binary scale yet are suspected or known to be influenced by imprinted genes, and there is consequent interest in seeking evidence for parent-of-origin effects at these loci. Here we show how simple linear models can be used to estimate these parent-of-origin effects for a broad class of phenotypes; in particular, normally distributed quantitative traits are easily dealt with.     

Phylogenetic analysis and comparative data: a test and review of evidence

The question is often raised whether it is statistically necessary to control for phylogenetic associations in comparative studies. To investigate this question, we explore the use of a measure of phylogenetic correlation, lambda, introduced by Pagel (1999), that normally varies between 0 (phylogenetic independence) and 1 (species' traits covary in direct proportion to their shared evolutionary history). Simulations show lambda to be a statistically powerful index for measuring whether data exhibit phylogenetic dependence or not and whether it has low rates of Type I error. Moreover, lambda is robust to incomplete phylogenetic information, which demonstrates that even partial information on phylogeny will improve the accuracy of phylogenetic analyses. To assess whether traits generally show phylogenetic associations, we present a quantitative review of 26 published phylogenetic comparative data sets. The data sets include 103 traits and were chosen from the ecological literature in which debate about the need for phylogenetic correction has been most acute. Eighty-eight percent of data sets contained at least one character that displayed significant phylogenetic dependence, and 60% of characters overall (pooled across studies) showed significant evidence of phylogenetic association. In 16% of tests, phylogenetic correlation could be neither supported nor rejected. However, most of these equivocal results were found in small phylogenies and probably reflect a lack of power. We suggest that the parameter lambda be routinely estimated when analyzing comparative data, since it can also be used simultaneously to adjust the phylogenetic correction in a manner that is optimal for the data set, and we present an example of how this may be done.     

Social traits,social networks and evolutionary biology

The social environment is both an important agent of selection for most organisms, and an emergent property of their interactions. As an aggregation of interactions among members of a population, the social environment is a product of many sets of relationships and so can be represented as a network or matrix. Social network analysis in animals has focused on why these networks possess the structure they do, and whether individuals’ network traits, representing some aspect of their social phenotype, relate to their fitness. Meanwhile, quantitative geneticists have demonstrated that traits expressed in a social context can depend on the phenotypes and genotypes of interacting partners, leading to influences of the social environment on the traits and fitness of individuals and the evolutionary trajectories of populations. Therefore, both fields are investigating similar topics, yet have arrived at these points relatively independently. We review how these approaches are diverged, and yet how they retain clear parallelism and so strong potential for complementarity. This demonstrates that, despite separate bodies of theory, advances in one might inform the other. Techniques in network analysis for quantifying social phenotypes, and for identifying community structure, should be useful for those studying the relationship between individual behaviour and group‐level phenotypes. Entering social association matrices into quantitative genetic models may also reduce bias in heritability estimates, and allow the estimation of the influence of social connectedness on trait expression. Current methods for measuring natural selection in a social context explicitly account for the fact that a trait is not necessarily the property of a single individual, something the network approaches have not yet considered when relating network metrics to individual fitness. Harnessing evolutionary models that consider traits affected by genes in other individuals (i.e. indirect genetic effects) provides the potential to understand how entire networks of social interactions in populations influence phenotypes and predict how these traits may evolve. By theoretical integration of social network analysis and quantitative genetics, we hope to identify areas of compatibility and incompatibility and to direct research efforts towards the most promising areas. Continuing this synthesis could provide important insights into the evolution of traits expressed in a social context and the evolutionary consequences of complex and nuanced social phenotypes.     

Effects of larval crowding on quantitative variation for development time and viability in Drosophila melanogaster

Competition between individuals belonging to the same species is a universal feature of natural populations and is the process underpinning organismal adaptation. Despite its importance, still comparatively little is known about the genetic variation responsible for competitive traits. Here, we measured the phenotypic variation and quantitative genetics parameters for two fitness‐related traits—egg‐to‐adult viability and development time—across a panel of Drosophila strains under varying larval densities. Both traits exhibited substantial genetic variation at all larval densities, as well as significant genotype‐by‐environment interactions (GEIs). GEI was attributable to changes in the rank order of reaction norms for both traits, and additionally to differences in the between‐line variance for development time. The coefficient of genetic variation increased under stress conditions for development time, while it was higher at both high and low densities for viability. While development time also correlated negatively with fitness at high larval densities—meaning that fast developers have high fitness—there was no correlation with fitness at low density. This result suggests that GEI may be a common feature of fitness‐related genetic variation and, further, that trait values under noncompetitive conditions could be poor indicators of individual fitness. The latter point could have significant implications for animal and plant breeding programs, as well as for conservation genetics.     

Estimating evolutionary parameters when viability selection is operating

Some individuals die before a trait is measured or expressed (the invisible fraction), and some relevant traits are not measured in any individual (missing traits). This paper discusses how these concepts can be cast in terms of missing data problems from statistics. Using missing data theory, I show formally the conditions under which a valid evolutionary inference is possible when the invisible fraction and/or missing traits are ignored. These conditions are restrictive and unlikely to be met in even the most comprehensive long-term studies. When these conditions are not met, many selection and quantitative genetic parameters cannot be estimated accurately unless the missing data process is explicitly modelled. Surprisingly, this does not seem to have been attempted in evolutionary biology. In the case of the invisible fraction, viability selection and the missing data process are often intimately linked. In such cases, models used in survival analysis can be extended to provide a flexible and justified model of the missing data mechanism. Although missing traits pose a more difficult problem, important biological parameters can still be estimated without bias when appropriate techniques are used. This is in contrast to current methods which have large biases and poor precision. Generally, the quantitative genetic approach is shown to be superior to phenotypic studies of selection when invisible fractions or missing traits exist because part of the missing information can be recovered from relatives.     

Quantitative trait locus analysis of longitudinal quantitative trait data in complex pedigrees

There is currently considerable interest in genetic analysis of quantitative traits such as blood pressure and body mass index. Despite the fact that these traits change throughout life they are commonly analyzed only at a single time point. The genetic basis of such traits can be better understood by collecting and effectively analyzing longitudinal data. Analyses of these data are complicated by the need to incorporate information from complex pedigree structures and genetic markers. We propose conducting longitudinal quantitative trait locus (QTL) analyses on such data sets by using a flexible random regression estimation technique. The relationship between genetic effects at different ages is efficiently modeled using covariance functions (CFs). Using simulated data we show that the change in genetic effects over time can be well characterized using CFs and that including parameters to model the change in effect with age can provide substantial increases in power to detect QTL compared with repeated measure or univariate techniques. The asymptotic distributions of the methods used are investigated and methods for overcoming the practical difficulties in fitting CFs are discussed. The CF-based techniques should allow efficient multivariate analyses of many data sets in human and natural population genetics.     

Evolutionary quantitative genetics of nonlinear developmental systems

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

TRANSLATING BETWEEN MICROEVOLUTIONARY PROCESS AND MACROEVOLUTIONARY PATTERNS: THE CORRELATION STRUCTURE OF INTERSPECIFIC DATA

As species evolve along a phylogenetic tree, we expect closely related species to retain some phenotypic similarities due to their shared evolutionary histories. The amount of expected similarity depends both on the hierarchical phylogenetic structure, and on the specific magnitude and types of evolutionary changes that accumulate during each generation. In this study, we show how models of microevolutionary change can be translated into the resulting macroevolutionary patterns. We illustrate how the structure of phenotypic covariances expected in interspecific measurements can be derived, and how this structure depends on the microevolutionary forces guiding phenotypic change at each generation. We then explore the covariance structure expected from several simple microevolutionary models of phenotypic evolution, including various combinations of random genetic drift, directional selection, stabilizing selection, and environmental change, as well as models of punctuated or burst-like evolution. We find that stabilizing selection leads to patterns of exponential decrease of between species covariance with phylogenetic distance. This is different from the usual linear patterns of decrease assumed in most comparative and systematic methods. Nevertheless, linear patterns of decrease can result from many processes in addition to random genetic drift, such as directional and fluctuating selection as well as modes of punctuated change. Our framework can be used to develop methods for (1) phylogenetic reconstruction; (2) inference of the evolutionary process from comparative data; and (3) conducting or evaluating statistical analyses of comparative data while taking phylogenetic history into account.