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, gmax, 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.

The properties of spontaneous mutations in the opportunistic pathogen Pseudomonas aeruginosa.

BACKGROUND: Natural genetic variation ultimately arises from the process of mutation. Knowledge of how the raw material for evolution is produced is necessary for a full understanding of several fundamental evolutionary concepts. We performed a mutation accumulation experiment with wild-type and mismatch-repair deficient, mutator lines of the pathogenic bacterium Pseudomonas aeruginosa, and used whole-genome sequencing to reveal the genome-wide rate, spectrum, distribution, leading/lagging bias, and context-dependency of spontaneous mutations. RESULTS: Wild-type base-pair mutation and indel rates were ~10(-10) and ~10(-11) per nucleotide per generation, respectively, and deficiencies in the mismatch-repair system caused rates to increase by over two orders of magnitude. A universal bias towards AT was observed in wild-type lines, but was reversed in mutator lines to a bias towards GC. Biases for which types of mutations occurred during replication of the leading versus lagging strand were detected reciprocally in both replichores. The distribution of mutations along the chromosome was non-random, with peaks near the terminus of replication and at positions intermediate to the replication origin and terminus. A similar distribution bias was observed along the chromosome in natural populations of P. aeruginosa. Site-specific mutation rates were higher when the focal nucleotide was immediately flanked by C:G pairings. CONCLUSIONS: Whole-genome sequencing of mutation accumulation lines allowed the comprehensive identification of mutations and revealed what factors of molecular and genomic architecture affect the mutational process. Our study provides a more complete view of how several mechanisms of mutation, mutation repair, and bias act simultaneously to produce the raw material for evolution.

Regularized regression can improve estimates of multivariate selection in the face of multicollinearity and limited data.

The breeder's equation, Δ z ¯ = G ß   , allows us to understand how genetics (the genetic covariance matrix, G) and the vector of linear selection gradients ß interact to generate evolutionary trajectories. Estimation of ß using multiple regression of trait values on relative fitness revolutionized the way we study selection in laboratory and wild populations. However, multicollinearity, or correlation of predictors, can lead to very high variances of and covariances between elements of ß, posing a challenge for the interpretation of the parameter estimates. This is particularly relevant in the era of big data, where the number of predictors may approach or exceed the number of observations. A common approach to multicollinear predictors is to discard some of them, thereby losing any information that might be gained from those traits. Using simulations, we show how, on the one hand, multicollinearity can result in inaccurate estimates of selection, and, on the other, how the removal of correlated phenotypes from the analyses can provide a misguided view of the targets of selection. We show that regularized regression, which places data-validated constraints on the magnitudes of individual elements of ß, can produce more accurate estimates of the total strength and direction of multivariate selection in the presence of multicollinearity and limited data, and often has little cost when multicollinearity is low. We also compare standard and regularized regression estimates of selection in a reanalysis of three published case studies, showing that regularized regression can improve fitness predictions in independent data. Our results suggest that regularized regression is a valuable tool that can be used as an important complement to traditional least-squares estimates of selection. In some cases, its use can lead to improved predictions of individual fitness, and improved estimates of the total strength and direction of multivariate selection.

Epicuticular compounds of Drosophila subquinaria and D. recens: identification, quantification, and their role in female mate choice.

The epicuticle of various Drosophila species consists of long-chain cuticular hydrocarbons (CHCs) and their derivatives that play a role in waterproofing and a dynamic means of chemical communication. Here, via gas chromatography and mass spectrometry, we identified and quantified the epicuticular composition of D. recens and D. subquinaria, two closely related species that show a pattern of reproductive character displacement in nature. Twenty-four compounds were identified with the most abundant, 11-cis-Vaccenyl acetate, present only in males of each species. Also exclusive to males were five tri-acylglycerides. The 18 remaining compounds were CHCs, all shared between the sexes and species. These CHCs were composed of odd carbon numbers (C29, C31, C33, and C35), with an increase in structural isomers in the C33 and C35 groups. Saturated hydrocarbons comprise only methyl-branched alkanes and were found only in the C29 and C31 groups. Alkenes were the least prevalent, with alkadienes dominating the chromatographic landscape in the longer chain lengths. Sexual dimorphism was extensive with 6/8 of the logcontrast CHCs differing significantly in relative concentration between males and females in D. recens and D. subquinaria, respectively. Males of the two species also differed significantly in relative concentration of six CHCs, while females differed in none. Female-choice mating trials revealed directional sexual selection on male CHCs in a population of each species, consistent with female mate preferences for these traits. The sexual selection vectors differed significantly in multivariate trait space, suggesting that different pheromone blends determine male attractiveness in each species.

Allometry constrains the evolution of sexual dimorphism in Drosophila across 33 million years of divergence.

Sexual dimorphism is widely viewed as adaptive, reflecting the evolution of males and females toward divergent fitness optima. Its evolution, however, may often be constrained by the shared genetic architecture of the sexes, and by allometry. Here, we investigated the evolution of sexual size dimorphism, shape dimorphism, and their allometric relationship, in the wings of 82 taxa in the family Drosophilidae that have been diverging for at least 33 million years. Shape dimorphism among species was remarkably similar, with males characterized by longer, thinner wings than females. There was, however, quantitative variation among species in both size and shape dimorphism, with evidence that they have adapted to different evolutionary optima in different clades on timescales of about 10 million years. Within species, shape dimorphism was predicted by size, and among species, there was a strong relationship between size dimorphism and shape dimorphism. Allometry constrained the evolution of shape dimorphism for the two most variable traits we studied, but dimorphism was evolutionary labile in other traits. The keys for disentangling alternative explanations for dimorphism evolution are studies of natural and sexual selection, together with a deeper understanding of how microevolutionary parameters of evolvability relate to macroevolutionary patterns of divergence.

Digest: Within and between sex covariances can enhance the response to climatic selection.

Do genetic covariances promote or impede rapid adaptation to changing environments? Hangartner et al. found that genetic covariances among traits and between sexes aligned with the inferred direction of selection along a latitudinal cline, suggesting that genetic covariances can augment the evolutionary response to climatic selection.

Cross-sex genetic covariances limit the evolvability of wing-shape within and among species of Drosophila.

The independent evolution of males and females is potentially constrained by both sexes inheriting the same alleles from their parents. This genetic constraint can limit the evolvability of complex traits; however, there are few studies of multivariate evolution that incorporate cross-sex genetic covariances in their predictions. Drosophila wing-shape has emerged as a model high-dimensional phenotype; wing-shape is highly evolvable in contemporary populations, and yet perplexingly stable across phylogenetic timescales. Here, we show that cross-sex covariances in Drosophila melanogaster, given by the B-matrix, may considerably bias wing-shape evolution. Using random skewers, we show that B would constrain the response to antagonistic selection by 90%, on average, but would double the response to concordant selection. Both cross-sex within-trait and cross-sex cross-trait covariances determined the predicted response to antagonistic selection, but only cross-sex within-trait covariances facilitated the predicted response to concordant selection. Similar patterns were observed in the direction of extant sexual dimorphism in D. melanogaster, and in directions of most and least dimorphic variation across the Drosophila phylogeny. Our results highlight the importance of considering between-sex genetic covariances when making predictions about evolution on both macro- and microevolutionary timescales, and may provide one more explanatory piece in the puzzle of stasis.

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.

Accounting for Sampling Error in Genetic Eigenvalues Using Random Matrix Theory.

The distribution of genetic variance in multivariate phenotypes is characterized by the empirical spectral distribution of the eigenvalues of the genetic covariance matrix. Empirical estimates of genetic eigenvalues from random effects linear models are known to be overdispersed by sampling error, where large eigenvalues are biased upward, and small eigenvalues are biased downward. The overdispersion of the leading eigenvalues of sample covariance matrices have been demonstrated to conform to the Tracy-Widom (TW) distribution. Here we show that genetic eigenvalues estimated using restricted maximum likelihood (REML) in a multivariate random effects model with an unconstrained genetic covariance structure will also conform to the TW distribution after empirical scaling and centering. However, where estimation procedures using either REML or MCMC impose boundary constraints, the resulting genetic eigenvalues tend not be TW distributed. We show how using confidence intervals from sampling distributions of genetic eigenvalues without reference to the TW distribution is insufficient protection against mistaking sampling error as genetic variance, particularly when eigenvalues are small. By scaling such sampling distributions to the appropriate TW distribution, the critical value of the TW statistic can be used to determine if the magnitude of a genetic eigenvalue exceeds the sampling error for each eigenvalue in the spectral distribution of a given genetic covariance matrix.

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.

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.

Reproductive character displacement of epicuticular compounds and their contribution to mate choice in Drosophila subquinaria and Drosophila recens.

Interactions between species can alter selection on sexual displays used in mate choice within species. Here we study the epicuticular pheromones of two Drosophila species that overlap partially in geographic range and are incompletely reproductively isolated. Drosophila subquinaria shows a pattern of reproductive character displacement against Drosophila recens, and partial behavioral isolation between conspecific sympatric versus allopatric populations, whereas D. recens shows no such variation in mate choice. First, using manipulative perfuming experiments, we show that females use pheromones as signals for mate discrimination both between species and among populations of D. subquinaria. Second, we show that patterns of variation in epicuticular compounds, both across populations and between species, are consistent with those previously shown for mating probabilities: pheromone compositions differ between populations of D. subquinaria that are allopatric versus sympatric with D. recens, but are similar across populations of D. recens regardless of overlap with D. subquinaria. We also identify differences in pheromone composition among allopatric regions of D. subquinaria. In sum, our results suggest that epicuticular compounds are key signals used by females during mate recognition, and that these traits have diverged among D. subquinaria populations in response to reinforcing selection generated by the presence of D. recens.

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

Directional selection is prevalent in nature, yet phenotypes tend to remain relatively constant, suggesting a limit to trait evolution. However, the genetic basis of this limit is unresolved. Given widespread pleiotropy, opposing selection on a trait may arise from the effects of the underlying alleles on other traits under selection, generating net stabilizing selection on trait genetic variance. These pleiotropic costs of trait exaggeration may arise through any number of other traits, making them hard to detect in phenotypic analyses. Stabilizing selection can be inferred, however, if genetic variance is greater among low- compared to high-fitness individuals. We extend a recently suggested approach to provide a direct test of a difference in genetic variance for a suite of cuticular hydrocarbons (CHCs) in Drosophila serrata. Despite strong directional sexual selection on these traits, genetic variance differed between high- and low-fitness individuals and was greater among the low-fitness males for seven of eight CHCs, significantly more than expected by chance. Univariate tests of a difference in genetic variance were nonsignificant but likely have low power. Our results suggest that further CHC exaggeration in D. serrata in response to sexual selection is limited by pleiotropic costs mediated through other traits.