Genetic independence between traits separated by metamorphosis is widespread but varies with biological function.

Why is metamorphosis so pervasive? Does it facilitate the independent (micro)evolution of quantitative traits in distinct life stages, similarly to how it enables some limbs and organs to develop at specific life stages? We tested this hypothesis by measuring the expression of 6400 genes in 41 Drosophila melanogaster inbred lines at larval and adult stages. Only 30% of the genes showed significant genetic correlations between larval and adult expression. By contrast, 46% of the traits showed some level of genetic independence between stages. Gene ontology terms enrichment revealed that across stages correlated traits were often involved in proteins synthesis, insecticide resistance and innate immunity, while a vast number of genes expression traits associated with energy metabolism were independent between life stages. We compared our results to a similar case: genetic constraints between males and females in gonochoric species (i.e. sexual antagonism). We expected selection for the separation between males and females to be higher than between juvenile and adult functions, as gonochorism is a more common strategy in the animal kingdom than metamorphosis. Surprisingly, we found that inter-stage constraints were lower than inter-sexual genetic constraints. Overall, our results show that metamorphosis enables a large part of the transcriptome to evolve independently at different life stages.

The phenome-wide distribution of genetic variance.

A general observation emerging from estimates of additive genetic variance in sets of functionally or developmentally related traits is that much of the genetic variance is restricted to few trait combinations as a consequence of genetic covariance among traits. While this biased distribution of genetic variance among functionally related traits is now well documented, how it translates to the broader phenome and therefore any trait combination under selection in a given environment is unknown. We show that 8,750 gene expression traits measured in adult male Drosophila serrata exhibit widespread genetic covariance among random sets of five traits, implying that pleiotropy is common. Ultimately, to understand the phenome-wide distribution of genetic variance, very large additive genetic variance-covariance matrices (G) are required to be estimated. We draw upon recent advances in matrix theory for completing high-dimensional matrices to estimate the 8,750-trait G and show that large numbers of gene expression traits genetically covary as a consequence of a single genetic factor. Using gene ontology term enrichment analysis, we show that the major axis of genetic variance among expression traits successfully identified genetic covariance among genes involved in multiple modes of transcriptional regulation. Our approach provides a practical empirical framework for the genetic analysis of high-dimensional phenome-wide trait sets and for the investigation of the extent of high-dimensional genetic constraint.

The measure and significance of Bateman's principles.

Bateman's principles explain sex roles and sexual dimorphism through sex-specific variance in mating success, reproductive success and their relationships within sexes (Bateman gradients). Empirical tests of these principles, however, have come under intense scrutiny. Here, we experimentally show that in replicate groups of red junglefowl, Gallus gallus, mating and reproductive successes were more variable in males than in females, resulting in a steeper male Bateman gradient, consistent with Bateman's principles. However, we use novel quantitative techniques to reveal that current methods typically overestimate Bateman's principles because they (i) infer mating success indirectly from offspring parentage, and thus miss matings that fail to result in fertilization, and (ii) measure Bateman gradients through the univariate regression of reproductive over mating success, without considering the substantial influence of other components of male reproductive success, namely female fecundity and paternity share. We also find a significant female Bateman gradient but show that this likely emerges as spurious consequences of male preference for fecund females, emphasizing the need for experimental approaches to establish the causal relationship between reproductive and mating success. While providing qualitative support for Bateman's principles, our study demonstrates how current approaches can generate a misleading view of sex differences and roles.

High-throughput phenotyping to characterise range use behaviour in broiler chickens.

A key characteristic of free-range chicken farming is to enable chickens to spend time outdoors. However, each chicken may use the available areas for roaming in variable ways. To check if, and how, broilers use their outdoor range at an individual level, we need to reliably characterise range use behaviour. Traditional methods relying on visual scans require significant time investment and only provide discontinuous information. Passive RFID (Radio Frequency Identification) systems enable tracking individually tagged chickens' when they go through pop-holes; hence, they only provide partial information on the movements of individual chickens. Here, we describe a new method to measure chickens' range use and test its reliability on three ranges each containing a different breed. We used an active RFID system to localise chickens in their barn, or in one of nine zones of their range, every 30 seconds and assessed range-use behaviour in 600 chickens belonging to three breeds of slow- or medium-growing broilers used for outdoor production (all <40 g daily weight gain). From those real-time locations, we determined five measures to describe daily range use: time spent in the barn, number of outdoor accesses, number of zones visited in a day, gregariousness (an index that increases when birds spend time in zones where other birds are), and numbers of zone changes. Principal Component Analyses (PCAs) were performed on those measures, in each production system, to create two synthetic indicators of chickens' range use behaviour. The first two PCA axes represented over 90% of the variance of the five measures and were both consistent over time and correlated with independent visual scans. Contributions of the five measures to the PCAs were similar among breeds, except for the correlation between the number of outdoor accesses and the four other measures. PC1 correlated with time spent inside the barn and zone changes frequency, whilst PC2 was explained by exploration of the range. Taken together, PC1 and PC2 indicators showed that range use increased with age, outdoor temperature (in spring), and did not differ between males and females. Importantly, daily scores for both indicators were repeatable among individuals - particularly in PC1 - showing inter-individual variability on range-use. The characterisation of broiler behaviour around their range with these reliable and repeatable indicators provides novel tools to help understand individual variations of range-use in free-range farming.

Mutational Pleiotropy and the Strength of Stabilizing Selection Within and Between Functional Modules of Gene Expression.

Variational modules, sets of pleiotropically covarying traits, affect phenotypic evolution, and therefore are predicted to reflect functional modules, such that traits within a variational module also share a common function. Such an alignment of function and pleiotropy is expected to facilitate adaptation by reducing the deleterious effects of mutations, and by allowing coordinated evolution of functionally related sets of traits. Here, we adopt a high-dimensional quantitative genetic approach using a large number of gene expression traits in Drosophila serrata to test whether functional grouping, defined by gene ontology (GO terms), predicts variational modules. Mutational or standing genetic covariance was significantly greater than among randomly grouped sets of genes for 38% of our functional groups, indicating that GO terms can predict variational modularity to some extent. We estimated stabilizing selection acting on mutational covariance to test the prediction that functional pleiotropy would result in reduced deleterious effects of mutations within functional modules. Stabilizing selection within functional modules was weaker than that acting on randomly grouped sets of genes in only 23% of functional groups, indicating that functional alignment can reduce deleterious effects of pleiotropic mutation but typically does not. Our analyses also revealed the presence of variational modules that spanned multiple functions.

Rapid evolution of the intersexual genetic correlation for fitness in Drosophila melanogaster.

Sexual antagonism (SA) arises when male and female phenotypes are under opposing selection, yet genetically correlated. Until resolved, antagonism limits evolution toward optimal sex-specific phenotypes. Despite its importance for sex-specific adaptation and existing theory, the dynamics of SA resolution are not well understood empirically. Here, we present data from Drosophila melanogaster, compatible with a resolution of SA. We compared two independent replicates of the "LHM " population in which SA had previously been described. Both had been maintained under identical, controlled conditions, and separated for around 200 generations. Although heritabilities of male and female fitness were similar, the intersexual genetic correlation differed significantly, being negative in one replicate (indicating SA) but close to zero in the other. Using population sequencing, we show that phenotypic differences were associated with population divergence in allele frequencies at nonrandom loci across the genome. Large frequency changes were more prevalent in the population without SA and were enriched at loci mapping to genes previously shown to have sexually antagonistic relationships between expression and fitness. Our data suggest that rapid evolution toward SA resolution has occurred in one of the populations and open avenues toward studying the genetics of SA and its resolution.

Transcriptome-wide effects of sexual selection on the fate of new mutations.

Sexual selection on males is predicted to have widespread effects on genetic variation as a consequence of the pleiotropic allelic effects on sexual and nonsexual traits. We manipulated the opportunity for sexual selection on males during 27 generations of mutation accumulation in inbred lines of Drosophila serrata, and used a microarray platform to investigate the effect of sexual selection on the expression of 2689 genes. While gene expression signal was, on average, higher in the absence of sexual selection, this difference was small (0.1%). In contrast, sexual selection impacted substantially on the mutational variance in gene expression. Over all genes, mutational variance in gene expression was, on average, 42% higher when sexual selection operated than when it was absent. Our results indicate that sexual selection on males can generate widespread effects across the genome. An increase in mutational variance without a corresponding change in mean suggested that most expression traits were unlikely to be under direct sexual selection. Instead, the mutational variance in gene expression traits is consistent with divergence generated by widespread pleiotropic associations with traits affecting male mating success.

Pleiotropic mutations are subject to strong stabilizing selection.

The assumption that pleiotropic mutations are more deleterious than mutations with more restricted phenotypic effects is an important premise in models of evolution. However, empirical evidence supporting this assumption is limited. Here, we estimated the strength of stabilizing selection on mutations affecting gene expression in male Drosophila serrata. We estimated the mutational variance (VM) and the standing genetic variance (VG) from two well-matched panels of inbred lines: a panel of mutation accumulation (MA) lines derived from a single inbred ancestral line and a panel of inbred lines derived from an outbred population. For 855 gene-expression traits, we estimated the strength of stabilizing selection as s = VM/VG. Selection was observed to be relatively strong, with 17% of traits having s > 0.02, a magnitude typically associated with life-history traits. Randomly assigning expression traits to five-trait sets, we used factor analytic mixed modeling in the MA data set to identify covarying traits that shared pleiotropic mutations. By assigning traits to the same trait sets in the outbred line data set, we then estimated s for the combination of traits affected by pleiotropic mutation. For these pleiotropic combinations, the median s was three times greater than s acting on the individual component traits, and 46% of the pleiotropic trait combinations had s > 0.02. Although our analytical approach was biased toward detecting mutations with relatively large effects, likely overestimating the average strength of selection, our results provide widespread support for the prediction that stronger selection can act against mutations with pleiotropic effects.

The nature and extent of mutational pleiotropy in gene expression of male Drosophila serrata.

The nature and extent of mutational pleiotropy remain largely unknown, despite the central role that pleiotropy plays in many areas of biology, including human disease, agricultural production, and evolution. Here, we investigate the variation in 11,604 gene expression traits among 41 mutation accumulation (MA) lines of Drosophila serrata. We first confirmed that these expression phenotypes were heritable, detecting genetic variation in 96% of them in an outbred, natural population of D. serrata. Among the MA lines, 3385 (29%) of expression traits were variable, with a mean mutational heritability of 0.0005. In most traits, variation was generated by mutations of relatively small phenotypic effect; putative mutations with effects of greater than one phenotypic standard deviation were observed for only 8% of traits. With most (71%) traits unaffected by any mutation, our data provide no support for universal pleiotropy. We further characterized mutational pleiotropy in the 3385 variable traits, using sets of 5, randomly assigned, traits. Covariance among traits chosen at random with respect to their biological function is expected only if pleiotropy is extensive. Taking an analytical approach in which the variance unique to each trait in the random 5-trait sets was partitioned from variance shared among traits, we detected significant (at 5% false discovery rate) mutational covariance in 21% of sets. This frequency of statistically supported covariance implied that at least some mutations must pleiotropically affect a substantial number of traits (>70; 0.6% of all measured traits).