Reduced genetic variance among high fitness individuals: inferring stabilizing selection on male sexual displays in Drosophila serrata

Integrating genomics and multivariate evolutionary quantitative genetics: a case study of constraints on sexual selection in Drosophila serrata

In evolutionary quantitative genetics, the genetic variance-covariance matrix, G, and the vector of directional selection gradients, β, are key parameters for predicting multivariate selection responses and genetic constraints. Historically, investigations of G and β have not overlapped with those dissecting the genetic basis of quantitative traits. Thus, it remains unknown whether these parameters reflect pleiotropic effects at individual loci. Here, we integrate multivariate genome-wide association study (GWAS) with G and β estimation in a well-studied system of multivariate constraint: sexual selection on male cuticular hydrocarbons (CHCs) in Drosophila serrata. In a panel of wild-derived re-sequenced lines, we augment genome-based restricted maximum likelihood to estimate G alongside multivariate single nucleotide polymorphism (SNP) effects, detecting 532 significant associations from 1 652 276 SNPs. Constraint was evident, with β lying in a direction of G with low evolvability. Interestingly, minor frequency alleles typically increased male CHC-attractiveness suggesting opposing natural selection on β. SNP effects were significantly misaligned with the major eigenvector of G, gmax, but well aligned to the second and third eigenvectors g2 and g3. We discuss potential factors leading to these varied results including multivariate stabilizing selection and mutational bias. Our framework may be useful as researchers increasingly access genomic methods to study multivariate selection responses in wild populations.

Why does allometry evolve so slowly?

Morphological allometry is striking due to its evolutionary conservatism, making it an example of a certain sort of evolutionary stasis. Organisms that vary in size, whether for developmental, environmental, or evolutionary reasons, adopt shapes that are predictable from that size alone. There are two major hypotheses to explain this. It may be that natural selection strongly favors each allometric pattern, or that organisms lack the development and genetic capacity to produce variant shapes for selection to act on. Using a high-throughput system for measuring the size and shape of Drosophila wings, we documented an allometric pattern that has been virtually unchanged for 40 million years. We performed an artificial selection experiment on the static allometric slope within one species. In just 26 generations, we were able to increase the slope from 1.1 to 1.4, and decrease it to 0.8. Once artificial selection was suspended, the slope rapidly evolved back to a value near the initial static slope. This result decisively rules out the hypothesis that allometry is preserved due to a lack of genetic variation, and provides evidence that natural selection acts to maintain allometric relationships. On the other hand, it seems implausible that selection on allometry in the wing alone could be sufficiently strong to maintain static allometries over millions of years. This suggests that a potential explanation for stasis is selection on a potentially large number of pleiotropic effects. This seems likely in the case of allometry, as the sizes of all parts of the body may be altered when the allometric slope of one body part is changed. Unfortunately, hypotheses about pleiotropy have been very difficult to test. We lay out an approach to begin the systematic study of pleiotropic effects using genetic manipulations and high-throughput phenotyping.

Do slower movers have lower reproductive success and higher mutation load?

Deleterious mutations occur frequently in eukaryotes, resulting in individuals carrying multiple alleles that decrease their fitness. At a population level, if unchecked, accumulation of this mutation load can ultimately lead to extinction. How selection counters the accumulation of mutation load, limiting declines in population fitness, is not well understood. Here, we use manipulative experiments in zebrafish (Danio rerio) to investigate the opportunities for selection on mutation load. Inducing high mutation load through mutagenesis, we applied one generation of within‐family selection on locomotor performance and characterized both the direct response to this selection and the indirect response of reproductive success. Offspring of slow swimming parents exhibited age‐dependent declines in swimming speed, whereas their cousins, with faster swimming parents, did not. This pattern mimics previously documented differences between high and low mutation load populations of zebrafish, suggesting that slow swimming siblings inherited (and transmitted) more mutations than their faster swimming siblings. Crosses among offspring of slow swimming fish had, on average, <75% of the reproductive success of crosses among offspring of fast swimming parents, or crosses of offspring of slow swimmers with offspring of fast swimmers. This evidence of mutationally correlated swimming speed and reproductive success reveals the potential for concordant selection on mutation load through different fitness components. There was no evidence that crosses within families (where parents potentially shared the same mutations inherited from their common ancestor) had lower reproductive success than crosses among families, suggesting that viability selection was not acting predominantly through lethal recessive homozygotes. Rather, patterns of reproductive success are suggestive of effects of mutation number per se on embryo viability. Overall, our results highlight the potential for early life mortality to remove deleterious mutations, and the need to account for this mortality when investigating the evolutionary dynamics of mutation load.

Artificial Selection to Increase the Phenotypic Variance in gmax Fails

Stabilizing selection is important in evolutionary theories of the maintenance of genetic variance and has been invoked as the key process determining macroevolutionary patterns of trait evolution. However, manipulative evidence for the extent of stabilizing selection, particularly on multivariate traits, is lacking. We used artificial disruptive selection in Drosophila serrata as a tool to determine the relative strength of stabilizing selection experienced by multivariate trait combinations with contrasting levels of genetic and mutational variance. Contrary to expectation, when disruptive selection was applied to the major axis of standing genetic variance, g_max, we observed a significant and repeatable decrease in its phenotypic variance. In contrast, the multivariate trait combination predicted to be under strong stabilizing selection showed a significant and repeatable increase in its phenotypic variance. Correlated responses were observed in all selection treatments, and viability selection operating on extreme phenotypes of traits genetically correlated with those directly selected on limited our ability to increase their phenotypic range. Our manipulation revealed that multivariate trait combinations were subject to stabilizing selection; however, we did not observe a direct relationship between the strength of stabilizing selection and the levels of standing genetic variance in multivariate trait combinations. Contrasting patterns of allele frequencies underlying traits with high versus low levels of standing genetic variance may be implicated in determining the response to artificial selection in multivariate trait combinations.

Heritable Micro-environmental Variance Covaries with Fitness in an Outbred Population of Drosophila serrata

The genetic basis of stochastic variation within a defined environment, and the consequences of such micro-environmental variance for fitness are poorly understood . Using a multigenerational breeding design in Drosophila serrata, we demonstrated that the micro-environmental variance in a set of morphological wing traits in a randomly mating population had significant additive genetic variance in most single wing traits. Although heritability was generally low (<1%), coefficients of additive genetic variance were of a magnitude typical of other morphological traits, indicating that the micro-environmental variance is an evolvable trait. Multivariate analyses demonstrated that the micro-environmental variance in wings was genetically correlated among single traits, indicating that common mechanisms of environmental buffering exist for this functionally related set of traits. In addition, through the dominance genetic covariance between the major axes of micro-environmental variance and fitness, we demonstrated that micro-environmental variance shares a genetic basis with fitness, and that the pattern of selection is suggestive of variance-reducing selection acting on micro-environmental variance.

Relaxation of herbivore-mediated selection drives the evolution of genetic covariances between plant competitive and defense traits

Insect herbivores are important mediators of selection on traits that impact plant defense against herbivory and competitive ability. Although recent experiments demonstrate a central role for herbivory in driving rapid evolution of defense and competition-mediating traits, whether and how herbivory shapes heritable variation in these traits remains poorly understood. Here, we evaluate the structure and evolutionary stability of the G matrix for plant metabolites that are involved in defense and allelopathy in the tall goldenrod, Solidago altissima. We show that G has evolutionarily diverged between experimentally replicated populations that evolved in the presence versus the absence of ambient herbivory, providing direct evidence for the evolution of G by natural selection. Specifically, evolution in an herbivore-free habitat altered the orientation of G, revealing a negative genetic covariation between defense- and competition-related metabolites that is typically masked in herbivore-exposed populations. Our results may be explained by predictions of classical quantitative genetic theory, as well as the theory of acquisition-allocation trade-offs. The study provides compelling evidence that herbivory drives the evolution of plant genetic architecture.

Simultaneous Estimation of Additive and Mutational Genetic Variance in an Outbred Population of Drosophila serrata

How new mutations contribute to genetic variation is a key question in biology. Although the evolutionary fate of an allele is largely determined by its heterozygous effect, most estimates of mutational variance and mutational effects derive from highly inbred lines, where new mutations are present in homozygous form. In an attempt to overcome this limitation, middle-class neighborhood (MCN) experiments have been used to assess the fitness effect of new mutations in heterozygous form. However, because MCN populations harbor substantial standing genetic variance, estimates of mutational variance have not typically been available from such experiments. Here we employ a modification of the animal model to analyze data from 22 generations of Drosophila serrata bred in an MCN design. Mutational heritability, measured for eight cuticular hydrocarbons, 10 wing-shape traits, and wing size in this outbred genetic background, ranged from 0.0006 to 0.006 (with one exception), a similar range to that reported from studies employing inbred lines. Simultaneously partitioning the additive and mutational variance in the same outbred population allowed us to quantitatively test the ability of mutation-selection balance models to explain the observed levels of additive and mutational genetic variance. The Gaussian allelic approximation and house-of-cards models, which assume real stabilizing selection on single traits, both overestimated the genetic variance maintained at equilibrium, but the house-of-cards model was a closer fit to the data. This analytical approach has the potential to be broadly applied, expanding our understanding of the dynamics of genetic variance in natural populations.

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.

Time flies: Time of day and social environment affect cuticular hydrocarbon sexual displays in Drosophila serrata

Recent work on Drosophila cuticular hydrocarbons (CHCs) challenges a historical assumption that CHCs in flies are largely invariant. Here, we examine the effect of time of day and social environment on a suite of sexually selected CHCs in Drosophila serrata. We demonstrate that males become more attractive to females during the time of day that flies are most active and when most matings occur, but females become less attractive to males during the same time of day. These opposing temporal changes may reflect differences in selection among the sexes. To evaluate the effect of social environment on male CHC attractiveness, we manipulated male opportunity for mating: male flies were housed either alone, with five females, with five males or with five males and five females. We found that males had the most attractive CHCs when with females, and less attractive CHCs when with competitor males. Social environment mediated how male CHC attractiveness cycled: males housed with females and/or other males showed temporal changes in CHC attractiveness, whereas males housed alone did not. In total, our results demonstrate temporal patterning of male CHCs that is dependent on social environment, and suggest that such changes may be beneficial to males.

Evolutionary constraints in high-dimensional trait sets

Genetic variation for individual traits is typically abundant, but for some multivariate combinations it is very low, suggesting that evolutionary limits might be generated by the geometric distribution of genetic variance. To test this prediction, we artificially selected along all eight genetic eigenvectors of a set of eight quantitative traits in Drosophila serrata. After six generations of 50% truncation selection, at least one replicate population of all treatments responded to selection, allowing us to reject a null genetic subspace as a cause of evolutionary constraint in this system. However, while all three replicate populations of the first five selection treatments displayed a significant response, the remaining three, characterized by low genetic variance in their selection indexes in the base population, displayed inconsistent responses to selection. The observation that only four of the nine replicate populations evolved in response to the direct selection applied to them in these low genetic variance treatments, led us to conclude that a nearly null subspace did limit evolution. Dimensions associated with low genetic variance are often found in multivariate analyses of standing genetic variance in morphological traits, suggesting that the nearly null genetic subspace may be a common mechanism of evolutionary constraint in nature.

A multivariate analysis of genetic variation in the advertisement call of the gray treefrog, Hyla versicolor

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

Comparing G: multivariate analysis of genetic variation in multiple populations

The additive genetic variance-covariance matrix (G) summarizes the multivariate genetic relationships among a set of traits. The geometry of G describes the distribution of multivariate genetic variance, and generates genetic constraints that bias the direction of evolution. Determining if and how the multivariate genetic variance evolves has been limited by a number of analytical challenges in comparing G-matrices. Current methods for the comparison of G typically share several drawbacks: metrics that lack a direct relationship to evolutionary theory, the inability to be applied in conjunction with complex experimental designs, difficulties with determining statistical confidence in inferred differences and an inherently pair-wise focus. Here, we present a cohesive and general analytical framework for the comparative analysis of G that addresses these issues, and that incorporates and extends current methods with a strong geometrical basis. We describe the application of random skewers, common subspace analysis, the 4th-order genetic covariance tensor and the decomposition of the multivariate breeders equation, all within a Bayesian framework. We illustrate these methods using data from an artificial selection experiment on eight traits in Drosophila serrata, where a multi-generational pedigree was available to estimate G in each of six populations. One method, the tensor, elegantly captures all of the variation in genetic variance among populations, and allows the identification of the trait combinations that differ most in genetic variance. The tensor approach is likely to be the most generally applicable method to the comparison of G-matrices from any sampling or experimental design.