Scott L. Allen

Genomic Evidence that Sexual Selection Impedes Adaptation to a Novel Environment

Sexual selection is widely appreciated for generating remarkable phenotypic diversity, but its contribution to adaptation and the purging of deleterious mutations is unresolved. To provide insight into the impact of sexual selection on naturally segregating polymorphisms across the genome, we previously evolved 12 populations of Drosophila serrata in a novel environment employing a factorial manipulation of the opportunities for natural and sexual selection. Here, we genotype more than 1,400 SNPs in the evolved populations and reveal that sexual selection affected many of the same genomic regions as natural selection, aligning with it as often as opposing it. Intriguingly, more than half of the 80 SNPs showing treatment effects revealed an interaction between natural and sexual selection. For these SNPs, while sexual selection alone often caused a change in allele frequency in the same direction as natural selection alone, when natural and sexual selection occurred together, changes in allele frequency were greatly reduced or even reversed. This suggests an antagonism between natural and sexual selection arising from male-induced harm to females. Behavioral experiments showed that males preferentially courted and mated with high-fitness females, and that the harm associated with this increased male attention eliminated the female fitness advantage. During our experiment, females carrying otherwise adaptive alleles may therefore have disproportionally suffered male-induced harm due to their increased sexual attractiveness. These results suggest that a class of otherwise adaptive mutations may not contribute to adaptation when mating systems involve sexual conflict and male mate preferences.

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).

The Genomic Distribution of Sex-Biased Genes in Drosophila serrata: X Chromosome Demasculinization, Feminization, and Hyperexpression in Both Sexes

The chromosomal distribution of genes with sex-biased expression is often nonrandom, and in species with XY sex chromosome systems, it is common to observe a deficit of X-linked male-biased genes and an excess of X-linked female-biased genes. One explanation for this pattern is that sex-specific selection has shaped the gene content of the X. Alternatively, the deficit of male-biased and excess of female-biased genes could be an artifact of differences between the sexes in the global expression level of their X chromosome(s), perhaps brought about by a lack of dosage compensation in males and hyperexpression in females. In the montium fruit fly, Drosophila serrata, both these explanations can account for a deficit of male-biased and excess of female-biased X-linked genes. Using genome-wide expression data from multiple male and female tissues (n = 176 hybridizations), we found that testis- and accessory gland-specific genes are underrepresented whereas female ovary-specific genes are overrepresented on the X chromosome, suggesting that X-linkage is disfavored for male function genes but favored for female function genes. However, genes with such sex-specific functions did not fully account for the deficit of male-biased and excess of female-biased X-linked genes. We did, however, observe sex differences in the global expression level of the X chromosome and autosomes. Surprisingly, and in contrast to other species where a lack of dosage compensation in males is responsible, we found that hyperexpression of X-linked genes in both sexes leads to this imbalance in D. serrata. Our results highlight how common genomic distributions of sex-biased genes, even among closely related species, may arise via quite different evolutionary processes.

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.

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.

Sex-biased transcriptome divergence along a latitudinal gradient

Sex‐dependent gene expression is likely an important genomic mechanism that allows sex‐specific adaptation to environmental changes. Among Drosophila species, sex‐biased genes display remarkably consistent evolutionary patterns; male‐biased genes evolve faster than unbiased genes in both coding sequence and expression level, suggesting sex differences in selection through time. However, comparatively little is known of the evolutionary process shaping sex‐biased expression within species. Latitudinal clines offer an opportunity to examine how changes in key ecological parameters also influence sex‐specific selection and the evolution of sex‐biased gene expression. We assayed male and female gene expression in Drosophila serrata along a latitudinal gradient in eastern Australia spanning most of its endemic distribution. Analysis of 11 631 genes across eight populations revealed strong sex differences in the frequency, mode and strength of divergence. Divergence was far stronger in males than females and while latitudinal clines were evident in both sexes, male divergence was often population specific, suggesting responses to localized selection pressures that do not covary predictably with latitude. While divergence was enriched for male‐biased genes, there was no overrepresentation of X‐linked genes in males. By contrast, X‐linked divergence was elevated in females, especially for female‐biased genes. Many genes that diverged in D. serrata have homologs also showing latitudinal divergence in Drosophila simulans and Drosophila melanogaster on other continents, likely indicating parallel adaptation in these distantly related species. Our results suggest that sex differences in selection play an important role in shaping the evolution of gene expression over macro‐ and micro‐ecological spatial scales.

Single-molecule sequencing of the Drosophila serrata genome

Long-read sequencing technology promises to greatly enhance de novo assembly of genomes for nonmodel species. Although the error rates of long reads have been a stumbling block, sequencing at high coverage permits the self-correction of many errors. Here, we sequence and de novo assemble the genome of Drosophila serrata, a species from the montium subgroup that has been well-studied for latitudinal clines, sexual selection, and gene expression, but which lacks a reference genome. Using 11 PacBio single-molecule real-time (SMRT cells), we generated 12 Gbp of raw sequence data comprising ∼65 × whole-genome coverage. Read lengths averaged 8940 bp (NRead50 12,200) with the longest read at 53 kbp. We self-corrected reads using the PBDagCon algorithm and assembled the genome using the MHAP algorithm within the PBcR assembler. Total genome length was 198 Mbp with an N50 just under 1 Mbp. Contigs displayed a high degree of chromosome arm-level conservation with the D. melanogaster genome and many could be sensibly placed on the D. serrata physical map. We also provide an initial annotation for this genome using in silico gene predictions that were supported by RNA-seq data.

Genetic constraints on microevolutionary divergence of sex-biased gene expression

The evolution of sex-specific phenotypes is an important dimension of diversification and local adaptation. The sex-dependent regulation of gene expression is considered a key genomic mechanism facilitating sex-dependent adaptation. In many species, genes with male-biased expression evolve faster in DNA sequence and expression level than genes with female-biased or sexually monomorphic expression. While positive selection may be responsible for rapid DNA sequence evolution, why expression of male-biased genes also evolves rapidly remains unclear. Beyond sex differences in selection, some aspects of the genetic architecture of gene expression could contribute to the rapid evolution of male-biased gene expression. First, male-biased genes might simply have greater standing genetic variance than female-biased genes. Second, male-biased genes could be less constrained by pleiotropy, either within or between sexes. Here, we evaluate these alternative explanations on an intraspecific scale using a series of quantitative genetic experiments conducted on natural variation in male and female gene expression in the fly Drosophila serrata. Male-biased genes had significantly higher genetic variance than female-biased genes and were generally more narrowly expressed across tissues, suggesting lower within-individual pleiotropy. However, consistent with stronger constraints due to between-sex pleiotropy, their between-sex genetic correlations, rMF, were higher than for female-biased genes and more strongly negatively associated with sex bias. Using an extensive clinal dataset, we tested whether sex differences in gene expression divergence among populations have been shaped by pleiotropy. Here too, male-biased gene divergence was more strongly associated with between-sex pleiotropy than was female-biased gene divergence. Systematic differences in genetic variance and pleiotropy may be important factors influencing sex-specific adaptation arising through changes in gene expression. This article is part of the theme issue ‘Linking local adaptation with the evolution of sex differences’.

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

Collet et al. adopt a high-dimensional quantitative genetic approach using gene expression traits to test for the presence of modularity of the genotype-phenotype map, where traits contributing to the same function (functional modularity)... 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.

A Genomic Reference Panel for Drosophila serrata

Here we describe a collection of re-sequenced inbred lines of Drosophila serrata, sampled from a natural population situated deep within the species endemic distribution in Brisbane, Australia. D. serrata is a member of the speciose montium group whose members inhabit much of south east Asia and has been well studied for aspects of climatic adaptation, sexual selection, sexual dimorphism, and mate recognition. We sequenced 110 lines that were inbred via 17-20 generations of full-sib mating at an average coverage of 23.5x with paired-end Illumina reads. 15,228,692 biallelic SNPs passed quality control after being called using the Joint Genotyper for Inbred Lines (JGIL). Inbreeding was highly effective and the average levels of residual heterozygosity (0.86%) were well below theoretical expectations. As expected, linkage disequilibrium decayed rapidly, with r2 dropping below 0.1 within 100 base pairs. With the exception of four closely related pairs of lines which may have been due to technical errors, there was no statistical support for population substructure. Consistent with other endemic populations of other Drosophila species, preliminary population genetic analyses revealed high nucleotide diversity and, on average, negative Tajima’s D values. A preliminary GWAS was performed on a cuticular hydrocarbon trait, 2-Me-C28 revealing 4 SNPs passing Bonferroni significance residing in or near genes. One gene Cht9 may be involved in the transport of CHCs from the site of production (oenocytes) to the cuticle. Our panel will facilitate broader population genomic and quantitative genetic studies of this species and serve as an important complement to existing D. melanogaster panels that can be used to test for the conservation of genetic architectures across the Drosophila genus.