Genetic and geographic population structure in the malaria vector, Anopheles farauti, provides a candidate system for pioneering confinable gene-drive releases.

Indoor insecticide applications are the primary tool for reducing malaria transmission in the Solomon Archipelago, a region where Anopheles farauti is the only common malaria vector. Due to the evolution of behavioural resistance in some An. farauti populations, these applications have become less effective. New malaria control interventions are therefore needed in this region, and gene-drives provide a promising new technology. In considering developing a population-specific (local) gene-drive in An. farauti, we detail the species' population genetic structure using microsatellites and whole mitogenomes, finding many spatially confined populations both within and between landmasses. This strong population structure suggests that An. farauti would be a useful system for developing a population-specific, confinable gene-drive for field release, where private alleles can be used as Cas9 targets. Previous work on Anopheles gambiae has used the Cardinal gene for the development of a global population replacement gene-drive. We therefore also analyse the Cardinal gene to assess whether it may be a suitable target to engineer a gene-drive for the modification of local An. farauti populations. Despite the extensive population structure observed in An. farauti for microsatellites, only one remote island population from Vanuatu contained fixed and private alleles at the Cardinal locus. Nonetheless, this study provides an initial framework for further population genomic investigations to discover high-frequency private allele targets in localized An. farauti populations. This would enable the development of gene-drive strains for modifying localised populations with minimal chance of escape and may provide a low-risk route to field trial evaluations.

Dengue and chikungunya virus loads in the mosquito Aedes aegypti are determined by distinct genetic architectures.

Aedes aegypti is the primary vector of the arboviruses dengue (DENV) and chikungunya (CHIKV). These viruses exhibit key differences in their vector interactions, the latter moving more quicky through the mosquito and triggering fewer standard antiviral pathways. As the global footprint of CHIKV continues to expand, we seek to better understand the mosquito's natural response to CHIKV-both to compare it to DENV:vector coevolutionary history and to identify potential targets in the mosquito for genetic modification. We used a modified full-sibling design to estimate the contribution of mosquito genetic variation to viral loads of both DENV and CHIKV. Heritabilities were significant, but higher for DENV (40%) than CHIKV (18%). Interestingly, there was no genetic correlation between DENV and CHIKV loads between siblings. These data suggest Ae. aegypti mosquitoes respond to the two viruses using distinct genetic mechanisms. We also examined genome-wide patterns of gene expression between High and Low CHIKV families representing the phenotypic extremes of viral load. Using RNAseq, we identified only two loci that consistently differentiated High and Low families: a long non-coding RNA that has been identified in mosquito screens post-infection and a distant member of a family of Salivary Gland Specific (SGS) genes. Interestingly, the latter gene is also associated with horizontal gene transfer between mosquitoes and the endosymbiotic bacterium Wolbachia. This work is the first to link the SGS gene to a mosquito phenotype. Understanding the molecular details of how this gene contributes to viral control in mosquitoes may, therefore, also shed light on its role in Wolbachia.

Maintenance of quantitative genetic variance in complex, multitrait phenotypes: the contribution of rare, large effect variants in 2 Drosophila species.

The interaction of evolutionary processes to determine quantitative genetic variation has implications for contemporary and future phenotypic evolution, as well as for our ability to detect causal genetic variants. While theoretical studies have provided robust predictions to discriminate among competing models, empirical assessment of these has been limited. In particular, theory highlights the importance of pleiotropy in resolving observations of selection and mutation, but empirical investigations have typically been limited to few traits. Here, we applied high-dimensional Bayesian Sparse Factor Genetic modeling to gene expression datasets in 2 species, Drosophila melanogaster and Drosophila serrata, to explore the distributions of genetic variance across high-dimensional phenotypic space. Surprisingly, most of the heritable trait covariation was due to few lines (genotypes) with extreme [>3 interquartile ranges (IQR) from the median] values. Intriguingly, while genotypes extreme for a multivariate factor also tended to have a higher proportion of individual traits that were extreme, we also observed genotypes that were extreme for multivariate factors but not for any individual trait. We observed other consistent differences between heritable multivariate factors with outlier lines vs those factors without extreme values, including differences in gene functions. We use these observations to identify further data required to advance our understanding of the evolutionary dynamics and nature of standing genetic variation for quantitative traits.

Natural variation at a single gene generates sexual antagonism across fitness components in Drosophila.

Mutations with conflicting fitness effects in males and females accumulate in sexual populations, reducing their adaptive capacity.1,2 Although quantitative genetic studies indicate that sexually antagonistic polymorphisms are common,3-5 their molecular basis and population genetic properties remain poorly understood.6,7 Here, we show in fruit flies how natural variation at a single gene generates sexual antagonism through phenotypic effects on cuticular hydrocarbon (CHC) traits that function as both mate signals and protectors against abiotic stress8 across a latitudinal gradient. Tropical populations of Drosophila serrata have polymorphic CHCs producing sexual antagonism through opposing but sex-limited effects on these two fitness-related functions. We dissected this polymorphism to a single fatty-acyl CoA reductase gene, DsFAR2-B, that is expressed in oenocyte cells where CHCs are synthesized. RNAi-mediated disruption of the DsFAR2-B ortholog in D. melanogaster oenocytes affected CHCs in a similar way to that seen in D. serrata. Population genomic analysis revealed that balancing selection likely operates at the DsFAR2-B locus in the wild. Our study provides insights into the genetic basis of sexual antagonism in nature and connects sexually varying antagonistic selection on phenotypes with balancing selection on genotypes that maintains molecular variation.

Adaptive divergence in shoot gravitropism creates hybrid sterility in an Australian wildflower.

Natural selection is responsible for much of the diversity we see in nature. Just as it drives the evolution of new traits, it can also lead to new species. However, it is unclear whether natural selection conferring adaptation to local environments can drive speciation through the evolution of hybrid sterility between populations. Here, we show that adaptive divergence in shoot gravitropism, the ability of a plant's shoot to bend upwards in response to the downward pull of gravity, contributes to the evolution of hybrid sterility in an Australian wildflower, Senecio lautus We find that shoot gravitropism has evolved multiple times in association with plant height between adjacent populations inhabiting contrasting environments, suggesting that these traits have evolved by natural selection. We directly tested this prediction using a hybrid population subjected to eight rounds of recombination and three rounds of selection in the field. Our experiments revealed that shoot gravitropism responds to natural selection in the expected direction of the locally adapted population. Using the advanced hybrid population, we discovered that individuals with extreme differences in gravitropism had more sterile crosses than individuals with similar gravitropic responses, which were largely fertile, indicating that this adaptive trait is genetically correlated with hybrid sterility. Our results suggest that natural selection can drive the evolution of locally adaptive traits that also create hybrid sterility, thus revealing an evolutionary connection between local adaptation and the origin of new species.

Highly Replicated Evolution of Parapatric Ecotypes.

Parallel evolution of ecotypes occurs when selection independently drives the evolution of similar traits across similar environments. The multiple origins of ecotypes are often inferred based on a phylogeny that clusters populations according to geographic location and not by the environment they occupy. However, the use of phylogenies to infer parallel evolution in closely related populations is problematic because gene flow and incomplete lineage sorting can uncouple the genetic structure at neutral markers from the colonization history of populations. Here, we demonstrate multiple origins within ecotypes of an Australian wildflower, Senecio lautus. We observed strong genetic structure as well as phylogenetic clustering by geography and show that this is unlikely due to gene flow between parapatric ecotypes, which was surprisingly low. We further confirm this analytically by demonstrating that phylogenetic distortion due to gene flow often requires higher levels of migration than those observed in S. lautus. Our results imply that selection can repeatedly create similar phenotypes despite the perceived homogenizing effects of gene flow.

The impact of artificial selection for Wolbachia-mediated dengue virus blocking on phage WO.

Wolbachia is currently at the forefront of global efforts to control arbovirus transmission from the vector Aedes aegypti. The use of Wolbachia relies on two phenotypes-cytoplasmic incompatibility (CI), conferred by cifA and cifB genes in prophage WO, and Wolbachia-mediated pathogen blocking (WMPB). These traits allow for local, self-sustaining reductions in transmission of dengue (DENV) following release of Wolbachia-infected A. aegypti. Here, aided by previous artificial selection experiment that generated Low and High pathogen blocking lines, we examined the potential link between WMPB and phage WO. We found no evidence that Wolbachia or phage WO relative densities predict DENV blocking strength across selected lines. However, selection resulted in reduced phage WO relative density for the Low WMPB line. The Low blocking line was previously shown to have reduced fitness as a result of selection. Through subsequent genomic analyses, we demonstrate that SNP variation underpinning selection for low blocking led to elevated frequency of potential deleterious SNPs on chromosome 1. The key region on chromosome 1 contains genes relating to cell cycle regulation, oxidative stress, transcriptional pausing, among others, that may have cascading effects on Wolbachia intracellular environment. We hypothesize that reduction in phage WO may be driven by changes in the loci directly under selection for blocking, or by the accumulation of predicted deleterious alleles in linkage disequilibrium with blocking loci resulting from hitchhiking. For the Low line with fewer phage WO, we also detected reduced expression of cifA and cifB CI genes, with patterns of expression varying between somatic and reproductive tissues. In conclusion, we propose that artificial selection for WMPB trait had corresponding impacts on phage WO densities, and also the transcription of CI-causing genes. Future studies may include a more detailed analysis of the regions the A. aegypti chromosome 1's ability to affect WMPB and other Wolbachia-associated intrinsic factors such as phage WO.

Artificial Selection Finds New Hypotheses for the Mechanism of Wolbachia-Mediated Dengue Blocking in Mosquitoes.

Wolbachia is an intracellular bacterium that blocks virus replication in insects and has been introduced into the mosquito, Aedes aegypti for the biocontrol of arboviruses including dengue, Zika, and chikungunya. Despite ongoing research, the mechanism of Wolbachia-mediated virus blocking remains unclear. We recently used experimental evolution to reveal that Wolbachia-mediated dengue blocking could be selected upon in the A. aegypti host and showed evidence that strong levels of blocking could be maintained by natural selection. In this study, we investigate the genetic variation associated with blocking and use these analyses to generate testable hypotheses surrounding the mechanism of Wolbachia-mediated dengue blocking. From our results, we hypothesize that Wolbachia may block virus replication by increasing the regeneration rate of mosquito cells via the Notch signaling pathway. We also propose that Wolbachia modulates the host's transcriptional pausing pathway either to prime the host's anti-viral response or to directly inhibit viral replication.

Selection on Aedes aegypti alters Wolbachia-mediated dengue virus blocking and fitness.

The dengue, Zika and chikungunya viruses are transmitted by the mosquito Aedes aegypti and pose a substantial threat to global public health. Current vaccines and mosquito control strategies have limited efficacy, so novel interventions are needed1,2. Wolbachia are bacteria that inhabit insect cells and have been found to reduce viral infection-a phenotype that is referred to as viral 'blocking'3. Although not naturally found in A. aegypti4, Wolbachia were stably introduced into this mosquito in 20114,5 and were shown to reduce the transmission potential of dengue, Zika and chikungunya6,7. Subsequent field trials showed Wolbachia's ability to spread through A. aegypti populations and reduce the local incidence of dengue fever8. Despite these successes, the evolutionary stability of viral blocking is unknown. Here, we utilized artificial selection to reveal genetic variation in the mosquito that affects Wolbachia-mediated dengue blocking. We found that mosquitoes exhibiting weaker blocking also have reduced fitness, suggesting the potential for natural selection to maintain blocking. We also identified A. aegypti genes that affect blocking strength, shedding light on a possible mechanism for the trait. These results will inform the use of Wolbachia as biocontrol agents against mosquito-borne viruses and direct further research into measuring and improving their efficacy.

The Transcriptional Response of Aedes aegypti with Variable Extrinsic Incubation Periods for Dengue Virus.

Dengue fever is the most prevalent arboviral disease globally. Dengue virus is transmitted primarily by the Aedes aegypti mosquito. One measure of the mosquito's efficiency as a vector is the extrinsic incubation period (EIP), which is the time between the ingestion of viremic blood and the emergence of virions in the saliva. The longer it takes virus to infect the midgut and traverse to the saliva, the fewer opportunities the mosquito will have to transmit the pathogen over its lifetime. We have shown previously that EIP for dengue virus is highly heritable and that it is negatively correlated with vector lifespan. Here, we examined the transcriptional profiles for mosquitoes that varied in their EIP phenotype and identified pathways associated with either short or long EIP. We found that mosquitoes with short EIP have less active immune responses but higher levels of protein translation and calcium ion homeostasis and that mosquitoes with longer EIP may have slower metabolism. These findings indicate a complex interplay between calcium ion distribution, ribosome biogenesis, and metabolism and reveal potential pathways that could be modified to slow the rate of viral progression and hence limit lifetime transmission capability.

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

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.

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.

Family level variation in Wolbachia-mediated dengue virus blocking in Aedes aegypti.

BACKGROUND: The mosquito vector Aedes aegypti is responsible for transmitting a range of arboviruses including dengue (DENV) and Zika (ZIKV). The global reach of these viruses is increasing due to an expansion of the mosquito's geographic range and increasing urbanization and human travel. Vector control remains the primary means for limiting these diseases. Wolbachia pipientis is an endosymbiotic bacterium of insects that has the ability to block the replication of pathogens, including flaviviruses such as DENV or ZIKV, inside the body of the vector. A strain of Wolbachia called wMel is currently being released into wild mosquito populations to test its potential to limit virus transmission to humans. The mechanism that underpins the virus blocking effect, however, remains elusive. METHODS: We used a modified full-sib breeding design in conjunction with vector competence assays in wildtype and wMel-infected Aedes aegypti collected from the field. All individuals were injected with DENV-2 intrathoracically at 5-6 days of age. Tissues were dissected 7 days post-infection to allow quantification of DENV and Wolbachia loads. RESULTS: We show the first evidence of family level variation in Wolbachia-mediated blocking in mosquitoes. This variation may stem from either genetic contributions from the mosquito and Wolbachia genomes or environmental influences on Wolbachia. In these families, we also tested for correlations between strength of blocking and expression level for several insect immunity genes with possible roles in blocking, identifying two genes of interest (AGO2 and SCP-2). CONCLUSIONS: In this study we show variation in Wolbachia-mediated DENV blocking in Aedes aegypti that may arise from genetic contributions and environmental influences on the mosquito-Wolbachia association. This suggests that Wolbachia-mediated blocking may have the ability to evolve through time or be expressed differentially across environments. The long-term efficacy of Wolbachia in the field will be dependent on the stability of blocking. Understanding the mechanism of blocking will be necessary for successful development of strategies that counter the emergence of evolved resistance or variation in its expression under diverse field conditions.

Sexual selection on spontaneous mutations strengthens the between-sex genetic correlation for fitness.

A proposed benefit to sexual selection is that it promotes purging of deleterious mutations from populations. For this benefit to be realized, sexual selection, which is usually stronger on males, must purge mutations deleterious to both sexes. Here, we experimentally test the hypothesis that sexual selection on males purges deleterious mutations that affect both male and female fitness. We measured male and female fitness in two panels of spontaneous mutation-accumulation lines of the fly, Drosophila serrata, each established from a common ancestor. One panel of mutation accumulation lines limited both natural and sexual selection (LS lines), whereas the other panel limited natural selection, but allowed sexual selection to operate (SS lines). Although mutation accumulation caused a significant reduction in male and female fitness in both the LS and SS lines, sexual selection had no detectable effect on the extent of the fitness reduction. Similarly, despite evidence of mutational variance for fitness in males and females of both treatments, sexual selection had no significant impact on the amount of mutational genetic variance for fitness. However, sexual selection did reshape the between-sex correlation for fitness: significantly strengthening it in the SS lines. After 25 generations, the between-sex correlation for fitness was positive but considerably less than one in the LS lines, suggesting that, although most mutations had sexually concordant fitness effects, sex-limited, and/or sex-biased mutations contributed substantially to the mutational variance. In the SS lines this correlation was strong and could not be distinguished from unity. Individual-based simulations that mimick the experimental setup reveal two conditions that may drive our results: (1) a modest-to-large fraction of mutations have sex-limited (or highly sex-biased) fitness effects, and (2) the average fitness effect of sex-limited mutations is larger than the average fitness effect of mutations that affect both sexes similarly.

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.

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.

Evolutionary potential of the extrinsic incubation period of dengue virus in Aedes aegypti.

Dengue fever is the most common arboviral disease worldwide. It is caused by dengue viruses (DENV) and the mosquito Aedes aegypti is its primary vector. One of the most powerful determinants of a mosquito's ability to transmit DENV is the length of the extrinsic incubation period (EIP), the time it takes for a virus to be transmitted by a mosquito after consuming an infected blood meal. Here, we repeatedly measured DENV load in the saliva of individual mosquitoes over their lifetime and used this in combination with a breeding design to determine the extent to which EIP might respond to the evolutionary forces of drift and selection. We demonstrated that genetic variation among mosquitoes contributes significantly to transmission potential and length of EIP. We reveal that shorter EIP is genetically correlated with reduced mosquito lifespan, highlighting negative life-history consequences for virus-infected mosquitoes. This work highlights the capacity for local genetic variation in mosquito populations to evolve and to dramatically affect the nature of human outbreaks. It also provides the impetus for isolating mosquito genes that determine EIP. More broadly, our dual experimental approach offers new opportunities for studying the evolutionary potential of transmission traits in other vector/pathogen systems.

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