Stabilizing selection has been predicted to change genetic variances and covariances so that the orientation of the genetic variance-covariance matrix (G) becomes aligned with the orientation of the fitness surface, but it is less clear how directional selection may change G. Here we develop statistical approaches to the comparison of G with vectors of linear and nonlinear selection. We apply these approaches to a set of male sexually selected cuticular hydrocarbons (CHCs) of Drosophila serrata. Even though male CHCs displayed substantial additive genetic variance, more than 99% of the genetic variance was orientated 74.9 degrees away from the vector of linear sexual selection, suggesting that open-ended female preferences may greatly reduce genetic variation in male display traits. Although the orientation of G and the fitness surface were found to differ significantly, the similarity present in eigenstructure was a consequence of traits under weak linear selection and strong nonlinear (convex) selection. Associating the eigenstructure of G with vectors of linear and nonlinear selection may provide a way of determining what long-term changes in G may be generated by the processes of natural and sexual selection.

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.

Determining the dimensionality of G provides an important perspective on the genetic basis of a multivariate suite of traits. Since the introduction of Fisher's geometric model, the number of genetically independent traits underlying a set of functionally related phenotypic traits has been recognized as an important factor influencing the response to selection. Here, we show how the effective dimensionality of G can be established, using a method for the determination of the dimensionality of the effect space from a multivariate general linear model introduced by Amemiya (1985). We compare this approach with two other available methods, factor-analytic modeling and bootstrapping, using a half-sib experiment that estimated G for eight cuticular hydrocarbons of Drosophila serrata. In our example, eight pheromone traits were shown to be adequately represented by only two underlying genetic dimensions by Amemiya's approach and factor-analytic modeling of the covariance structure at the sire level. In contrast, bootstrapping identified four dimensions with significant genetic variance. A simulation study indicated that while the performance of Amemiya's method was more sensitive to power constraints, it performed as well or better than factor-analytic modeling in correctly identifying the original genetic dimensions at moderate to high levels of heritability. The bootstrap approach consistently overestimated the number of dimensions in all cases and performed less well than Amemiya's method at subspace recovery.

In many species, females display preferences for extreme male signal traits, but it has not been determined if such preferences evolve as a consequence of females gaining genetic benefits from exercising choice. If females prefer extreme male traits because they indicate male genetic quality that will enhance the fitness of offspring, a genetic correlation will evolve between female preference genes and genes that confer offspring fitness. We show that females of Drosophila serrata prefer extreme male cuticular hydrocarbon (CHC) blends, and that this preference affects offspring fitness. Female preference is positively genetically correlated with offspring fitness, indicating that females have gained genetic benefits from their choice of males. Despite male CHCs experiencing strong sexual selection, the genes underlying attractive CHCs also conferred lower offspring fitness, suggesting a balance between sexual selection and natural selection may have been reached in this population.

Sexual selection in natural populations acts on highly heritable traits and tends to be relatively strong, implicating sexual selection as a causal agent in many phenotypic radiations. Sexual selection appears to be ineffectual in promoting phenotypic divergence among contemporary natural populations, however, and there is little evidence from artificial selection experiments that sexual fitness can evolve. Here, we demonstrate that a multivariate male trait preferred by Drosophila serrata females can respond to selection and results in the maintenance of male mating success. The response to selection was associated with a gene of major effect increasing in frequency from 12 to 35% in seven generations. No further response to selection, or increase in frequency of the major gene, was observed between generations 7 and 11, indicating an evolutionary limit had been reached. Genetic analyses excluded both depletion of genetic variation and overdominance as causes of the evolutionary limit. Relaxing artificial selection resulted in the loss of 52% of the selection response after a further five generations, demonstrating that the response under artificial sexual selection was opposed by antagonistic natural selection. We conclude that male D. serrata sexually selected traits, and attractiveness to D. serrata females conferred by these traits, were held at an evolutionary limit by the lack of genetic variation that would allow an increase in sexual fitness while simultaneously maintaining nonsexual fitness. Our results suggest that sexual selection is unlikely to cause divergence among natural populations without a concomitant change in natural selection, a conclusion consistent with observational evidence from natural populations.genetic variance | genetic constraint