Combining experimental evolution with next-generation sequencing: a powerful tool to study adaptation from standing genetic variation.

Evolve and resequence (E&R) is a new approach to investigate the genomic responses to selection during experimental evolution. By using whole genome sequencing of pools of individuals (Pool-Seq), this method can identify selected variants in controlled and replicable experimental settings. Reviewing the current state of the field, we show that E&R can be powerful enough to identify causative genes and possibly even single-nucleotide polymorphisms. We also discuss how the experimental design and the complexity of the trait could result in a large number of false positive candidates. We suggest experimental and analytical strategies to maximize the power of E&R to uncover the genotype-phenotype link and serve as an important research tool for a broad range of evolutionary questions.

Experimental evolution has a long tradition in biology (Garland and Rose, 2009). By exposing an evolving population to conditions chosen by the researcher, it is possible to study the response to this selection regime. A recent review highlighted the broad range of applications that have been investigated with this methodology and concluded that the breadth of research questions is only limited by the creativity of the experimenter (Kawecki ). In addition to the great diversity of experimental designs, experimental evolution provides a unique advantage compared with other evolutionary analyses: the ability to replicate an experiment under identical conditions. Through this replication, experimenters are able to distinguish between stochastic and deterministic effects. Until recently, experimental evolution has mainly focused on phenotypes, sometimes combined with the analysis of a small number of markers (see, for example, Nuzhdin ; Teotonio ). In the wake of the latest sequencing technologies and the ongoing drop in DNA sequencing costs, however, the ultimate goal to connect the phenotypic response to the underlying genetic changes during an experimental evolution study has now come within reach.

Depending on the starting population, two conceptually different approaches of experimental evolution can be distinguished. Either the experiment starts from a genetically homogeneous (invariable) population or from a polymorphic population. In the first approach, adaptation occurs through the accumulation of new beneficial mutations during the experiment (Elena and Lenski, 2003). These experiments therefore require very large population sizes and many generations to ensure a sufficient mutation supply and are thus largely restricted to microorganisms. Alternatively, experiments starting with a polymorphic population do not require novel mutations as selection can act on beneficial alleles that are already present at the beginning of the experiment. Given the massive genetic variation that is present in the starting population, the key challenge for this approach is distinguishing between selected and neutral variants. Neither randomly selected markers nor whole genome sequencing of a few representative individuals can provide sufficient information about the true target(s) of selection. Rather, genome-wide polymorphism data are needed.

As whole genome sequencing is still not feasible for large numbers of individuals, experimental evolution studies starting from polymorphic base populations rely on a modified next-generation sequencing approach. Rather than sequencing individuals separately, DNA of multiple individuals from a population are sequenced together (Pool-Seq). This method is more cost effective than sequencing of individuals (Futschik and Schlötterer, 2010) and yields highly accurate genome-wide allele frequency estimates (reviewed in Rellstab ; Schlötterer ). The combination of experimental evolution with Pool-Seq is also known as Evolve and Resequence (E&R; Turner ; Figure 1). Here, we review the state of the art of whole genome polymorphism analysis in experimental evolution studies relying primarily on segregating variation in the starting population.

Figure 1

Overview of E&R studies. (a) A population of flies is exposed for 60 generations to ultraviolet (UV) radiation (purple arrows). We assume here, for the sake of illustration, that darker pigmentation is beneficial in high UV environments, whereby darker flies will increase in frequency. (b) At the genotypic level, the allele frequency of the causative allele (dark brown) will increase, more so than hitchhiking variants (dark gray background) that will be recombined onto other backgrounds (breaks between dark and light gray background). (c) The allele frequencies of the starting population and the selected population are measured with Pool-Seq. (d) Causative variants can be identified by contrasting the allele frequencies between base and selected population and visualized with Manhattan plots. A full color version of this figure is available at the Heredity journal online.

In many experimental evolution studies, researchers select for a well-defined trait in a controlled environment. This assures that both the phenotypic and the underlying genomic response are triggered either directly or indirectly by the selection regime applied during the experiment. Thus, E&R studies provide a complementary approach to genome-wide association studies (GWASs) and linkage mapping experiments as strategies to connect genotype and phenotype.

Support for experimental evolution results by independent methods

Despite its conceptual appeal, E&R studies face a lot of uncharted territory. For example, guidelines for experimental design and data analysis (Box 1) were not available for the first studies. Therefore, several E&R studies used additional techniques to provide independent evidence to support the E&R results. In the following, we will highlight the results of E&R studies based on truncating selection and accompanying validation strategies adopted by some of these studies.

Hypoxia tolerance

Having selected a Drosophila melanogaster population for an increased ability to tolerate low oxygen concentrations over 200 generations, Zhou identified 188 candidate genes located in genomic regions that responded to the selection regime. Of these genes, 28 were previously implicated in hypoxia or similar phenotypes and 12 were linked to the Notch pathway. As previous gene expression studies had associated the Notch pathway with hypoxia, the authors concluded that they had successfully identified genes involved in this trait. We note, however, that the genomic regions reported to be responding to selection were rather large, probably because of a small number of founder haplotypes (see below). Hence, it is not clear how many of the candidate genes were actually selected in the population during the study.

Genetic basis of aging

Experimental evolution has a long tradition in the study of aging and other life history traits. Remolina selected long-lived flies for 50 generations and compared them with unselected controls propagated in a similar manner. Contrary to many other studies that inferred selection on the basis of individual single-nucleotide polymorphisms (SNPs), this study searched for regions of reduced variability in 1 kb windows. In total, 156 genes were identified to show the signature of adaptive divergence between selected and unselected lines. To validate these candidates, the authors measured gene expression divergence between both groups. Using a false discovery rate of <0.1, 25 candidate genes were found to be differentially expressed.

Parasitoid resistance

Following just five generations of selection for resistance against the parasitoid Asobara tabida, a consistent level of resistance was seen across all four replicates of a D. melanogaster population (Jalvingh ). Contrasting allele frequencies in selected and control flies using a similar window-based approach as Remolina , the authors concluded that <5% of the genome was influenced by selection. Among 345 genes located in the selected regions, 91 could be linked to pathways associated with immune response. Using two independent expression analyses related to A. tabida resistance, the authors found that some differentially expressed genes were located in selected regions, but no significant overlap between their data set and either of the expression analyses could be detected.

Courtship song in Drosophila

The vibration of wings is an important courtship signal in fruit flies, and has species-specific features. Male courtship song is characterized by several features, including the length of the interpulse interval (IPI). After demonstrating variation for this trait among D. melanogaster strains, Turner and Miller (2012) selected from a polymorphic base population for short and long IPI. After 14 generations, the experimental populations had diverged for IPI and were subjected to Pool-Seq. Despite a low empirical false discovery rate of 0.42%, >13 000 variants were significantly differentiated between short and long IPI flies. A significant under-representation of candidate SNPs on the X chromosome was also observed. In a subsequent study, the authors performed a GWAS based on the same set of lines that were used to generate the starting population for the experimental evolution study (Turner ). Although no SNP was significant after correction for multiple testing, SNPs with a high differentiation in the E&R study also tended to have low P-values in the GWAS. Conversely, none of the most significant SNPs in the GWAS were found among the 13 000 most differentiated ones in the experimental evolution study. This was taken as evidence for variation in IPI being caused by many loci, rather than a small number of large effect genes. Validation of two candidate genes by quantitative complementation tests revealed that one of them, Syntropin-like 1, had a small, but significant, effect on IPI.

Drosophila C virus resistance

The Drosophila C virus (DCV) is one of the best-studied pathogens of D. melanogaster. Exposing a natural population sample for 20 generations to DCV resulted in an increased survival after infection, suggesting that resistance alleles increased in frequency (Martins ). By applying Pool-Seq in four replicate populations, Martins identified two genomic regions where a variant increased in frequency in the selected populations. Interestingly, one of these regions had also been identified in a previous GWAS (Magwire ) and in both studies the same SNP in the gene pastrel was identified as the most likely target of selection. In addition, the involvement of Ubc-E2H in the second candidate region was validated with RNA interference. With the same SNP being indentified with E&R and GWAS, the DCV resistance is probably the most convincing example of E&R having identified the causative link between genotype and phenotype.

Evolutionary inferences of selection trajectories

The E&R method can potentially offer much more than measuring differentiated allele frequencies between two selection regimes or between selected and control populations. By sampling evolving populations at multiple time points, it is also possible to study the trajectories of the selected alleles and thus elucidate their evolutionary dynamics. Such trajectories have been largely studied in a modified experimental evolution design termed laboratory natural selection (Garland and Rose, 2009). Rather than selecting for a specific phenotype, this approach exposes populations to a defined environment where, as in nature, better adapted individuals have a higher reproductive success. Surprisingly, all studies of allele frequency dynamics detected a similar behavior of selected alleles.

Parts generated a polymorphic population of recombinant cells from two diverged yeast strains and then subjected it to high temperature for up to 12 generations. By following allele frequency changes during adaptation, the authors were not only able to pinpoint 21 selected genomic regions, but they also captured interesting dynamics for these loci that were not compatible with classic directional selection. Rather than increasing in frequency until becoming fixed, most favored alleles plateaued at intermediate frequencies. This reduction in the selection coefficient was later confirmed via an elegant population genetic model (Illingworth ).

A similar pattern has also been reported for two D. melanogaster E&R studies. The first study exposed a D. melanogaster population to a novel high-temperature environment and sampled allele frequency changes at two different time points, specifically after 15 and 37 generations (Orozco-terWengel ). Among all SNPs with allele frequency changes greater than expected under genetic drift during the experiment, the authors focused on the 2000 SNPs showing the most significant change across all three replicates. Although in the first 15 generations the majority of the candidate SNPs experienced a frequency increase of ~28%, in subsequent generations the allele frequencies had plateaued without becoming fixed. Most importantly, the authors also ruled out the possibility that this pattern was an analytical artifact (Orozco-terWengel ). The second D. melanogaster study, which did not analyze time series data, reported a large genomic response in flies selected for accelerated development over 600 generations, but found little support for selective sweeps resulting in their fixation (Burke ).

The reason for these puzzling dynamics is not yet understood. In experiments with changing environments, such as the fluctuating temperature used by Orozco-terWengel , marginal overdominance may explain the plateauing of selected alleles. In the other two experiments where selection was constant, thus marginal overdominance cannot serve as a universal explanation. Alternatively, recessive deleterious alleles or heterozygous advantage may explain the plateaus in the Drosophila data, but not in the haploid yeast strains. Finally, it has been proposed that the observed pattern could be explained by selection on a complex trait with several contributing loci: pronounced allele frequency changes are expected as long as the trait is far away from the fitness optimum, but slows down as the optimum is being approached (Chevin and Hospital, 2008). Further empirical testing is required to distinguish between the different explanations for the plateauing of putatively selected SNPs.

The mystery of the large number of candidate SNPs

One common observation in all E&R studies is that a massive number of candidate SNPs are identified, even after rigorous correction for multiple testing. Importantly, such large numbers of selected SNPs are not compatible with the observed large frequency changes, wherein 30% increases for selected alleles are not uncommon (Smith, 1968; Nuzhdin and Turner, 2013). One apparent explanation for the large number of candidate SNPs is that selection may act on a moderate number of loci that drag along many linked neutral variants, a phenomenon known as hitchhiking. In particular, studies that rely on either small experimental population sizes or have starting populations with high levels of linkage disequilibrium are expected to show in a selection signature comprising broad regions of adjacent SNPs.

Although this pattern can be clearly seen in some E&R studies, it is not sufficient to explain the patterns observed for populations with large population sizes. For example, two studies on flies exposed to new thermal environments (Orozco-terWengel ; Tobler ) tested explicitly whether narrow-range linkage could explain the excess of significant SNPs, but found that only SNPs within ±200 bp of the focal SNPs were affected, ruling this out as a general explanation. Tobler further investigated this question by comparing independent sets of replicates that were started from the same base population and had subsequently evolved independently from each other under the same selection regime. The authors found a very good general concordance between replicates for SNPs with allele frequency changes deviating from neutral expectations. However, this concordance was also apparent when only short introns were analyzed. Because short introns are, to a good approximation, evolving neutrally in D. melanogaster, no significant concordance is expected among sites located in these regions. The authors concluded that instead long-range linkage disequilibrium with selected sites may be responsible for the correlated response of SNPs located in short introns. Such long-range linkage disequilibrium could result from segregating chromosomal inversions that are common in D. melanogaster. In fact, using inversion-specific SNP markers, it has been shown that some inversion frequency changes in these experimental populations were probably driven by selection (Kapun ). Inversions are not the only cause of long-range linkage disequilibrium, however. Beneficial alleles occurring at a low frequency in the starting population will, by chance, have an association with all SNPs private to the haplotypes upon which the beneficial variant occurs. The lower the starting frequency of the beneficial allele, the more spurious long-range associations will be generated. By using individual-based computer simulations, Tobler demonstrated that such long-range linkage disequilibrium does result in many false positives because of linkage extending over several megabases. Consistent with this idea, a 1-Mb genomic region on chromosome 3R was found to harbor a large number of candidate SNPs with many putatively selected alleles in this region having risen from low frequencies in the starting population to high frequencies in the evolved populations (Orozco-terWengel ; Tobler ).

The future of E&R

Based on the results of recent E&R studies, it is apparent that E&R could be a powerful method to complement ongoing linkage mapping and GWAS approaches (Table 1). This has been demonstrated by a recent E&R study that identified the causative SNP for at least one gene determining a trait with a simple genetic basis (DCV resistance) (Martins ), whereas complementary results were obtained for a GWAS and an E&R study on a more complex trait (see, for example, Turner ). Nevertheless, is has also become clear that E&R faces its own specific challenges that need to be considered when interpreting the data, some of them are discussed below.

Table 1

Features of different approaches aiming to link genotype and phenotype

	E&R	Classic GWAS	GWAS in reference panel	Pool-GWAS	Linkage mapping
Analysis of heterozygous individuals	+	+	−	+	+
Repeated phenotyping	Every generation	−	+	−	−
Sensitivity to environmental noise	Low because of repeated phenotyping in every generation and replication	High	Low because of repeated phenotyping of identical genotypes	Moderate because of replication	High
Well-established analysis strategies	−	+	+	−	+
Mapping resolution	High	High	High	High	Moderate (cost effective), high (expensive)
Genetic diversity analyzed	High	High	Moderate–high	High	Limited to parental genotypes
Inference of effect size	Selection, coefficient	+	+	−	+
Randomized genetic background	+, in starting population and repeated mixing by sexual reproduction during the experiment	+, but sensitive to population structure that can be accounted for in analysis	+, but sensitive to population structure that can be accounted for in analysis given a sufficient sample size	+, but sensitive to population structure that can be accounted for in analysis	+
Genotyping/sequencing costs	Low because of Pool-Seq	High for establishment, no costs for follow-up experiments	High for establishment, no costs for follow-up experiments	Low because of Pool-Seq	Low because of the use of genetic markers (Rad-Tag sequencing)
Sampling effort	High because of maintenance of experimental populations	Depends on species	High for establishment, low later on	Moderate	Moderate
Analysis of multiple traits from the same genotypes	−	+	+	−	+
Replication	Yes, is common practice	Only across different populations	Requires an independent reference panel	Yes, it is common practice. Easy to expand to multiple populations	Requires independent mapping families
Influence of allele frequency (conditional on presence in the sample)	Low power for high-frequency alleles, low-frequency alleles are often lost	Yes	Yes	Yes	No
Trajectories of selected variants	+	−	−	−	−
Identification of adaptive variants in a defined environment	+	−	−	−	−

Abbreviations: E&R, Evolve and resequence; GWAS, genome-wide association study; Pool-Seq, Pool-sequencing.

Towards improved experimental designs

Current E&R studies employ a diverse array of experimental designs, but until very recently no guidance was available on how to optimize the power of these designs to detect selected loci. Three forward simulation studies (Baldwin-Brown ; Kofler and Schlötterer, 2014; Kessner and Novembre, 2014b) have explored the most important factors for an optimal experimental design. All three studies showed that increasing the number of replicates and experimental population size resulted in a higher power to detect selected loci. The strength of selection was also found to have a major impact, with both very strong and very weak selection being problematic (Kofler and Schlötterer, 2014). Although weakly selected sites failed to show a detectable allele frequency change, strong selection caused the fixation of many linked neutral variants, precluding the identification of the causative SNP. Furthermore, the detection of selected alleles becomes more difficult as the experiment continues, as causative SNPs eventually become fixed in the population while drift gradually reduces the signal-to-noise ratio.

The history of the starting population also has an important influence on the power of the study. The results showed that the amount of variation in the starting population is key. Experiments with starting populations using as many independent lines as possible had the highest power because of the low level of linkage disequilibrium. Another approach to reduce linkage disequilibrium is pre-experiment cultivation of the starting population in the laboratory. Computer simulations showed that this experimental approach resulted in the loss of favorable alleles and that the increase in power was moderate with strongly selected alleles benefitting most from this strategy (Kofler and Schlötterer, 2014). Importantly, laboratory adaptation during the pre-experiment cultivation probably does not confound the subsequent analysis: in recent E&R study, the same starting population was selected in two different environmental conditions, and very few SNPs appeared to be selected in both treatments, and this is unexpected if laboratory adaptation is important (Tobler ).

One further factor influencing the power of a study is the sequence coverage. Although for strongly selected sites a coverage of 50 × is fully sufficient, weakly selected sites require a substantially higher coverage (up to 200 × ) to estimate the allele frequencies to a level of precision that permits the reliable detection of small frequency differences (Kofler and Schlötterer, 2014) (Box 1). Nevertheless, even when a large number of loci are selected, almost 60% of the target sites can be identified using an appropriate experimental design (Kofler and Schlötterer, 2014). This aptly demonstrates the enormous potential of experimental evolution to identify the target(s) of selection. One further strategy to improve the performance of E&R studies is to include haplotype information (Kessner ). Current methods require knowledge of the haplotypes in the starting population, however, that will become increasingly difficult with larger number of founder chromosomes.

Because of its compact genome, high-quality reference genomic sequence, short generation time and ease of cultivation, D. melanogaster has been frequently used for E&R studies. However, D. melanogaster harbors many segregating inversions that could negatively affect the power of experimental evolution. Therefore, we suggest that future experiments make use of D. simulans, a close relative of D. melanogaster, that is almost free of segregating inversions (Aulard ). Adding to the attraction of D. simulans, it has a substantially improved reference genome (Hu ; Palmieri ), and with latest advances in genome editing (Liu ; Terns and Terns, 2014) rigorous functional testing can also be applied outside of genetic model organisms. In addition to changing the focal species, we strongly recommend increasing the number of replicates, number of founder chromosomes and the experimental population size. Although population sizes of ∼1000 individuals are currently at the upper end of Drosophila E&R studies, it is possible to increase this number by an order of magnitude. Not only is this expected to have an impact on the dynamics of phenotypic change (Weber, 1996), but also on the accuracy of the identification of targets of selection (Baldwin-Brown ; Kofler and Schlötterer, 2014). In order to identify causative variants with E&R, several test statistics have been developed (Box 2), some of which show remarkable differences in statistical power under a given evolutionary scenario (Figure 2). Furthermore, new statistical approaches that take full advantage of trajectories from multiple time points and across several replicates have the potential to increase the power of E&R studies substantially (Terhorst and Song, 2014; Topa ).

Figure 2

Performance of different test statistics used in E&R studies. Receiver operator characteristic (ROC) curves that contrast the true positive rate with the false positive rate. We extended the results of Kofler and Schlötterer (2014) by including the pooled Hs/D test, the pooled Hs/Hc test and the pooled FST test (Remolina ). Briefly, Kofler and Schlötterer (2014) simulated E&R with a base population that captures the pattern of polymorphism in a natural D. melanogaster population. They simulated 60 generations of selection with a population size of 1000 and 3 replicates. Results are shown for SNP-based analysis (top graphs) and for a window-based (bottom graphs) analyses using either 150 strongly (left graphs) or 150 weakly (right graphs) selected loci. The behavior of the Cochran–Mantel–Haenszel (CMH) and pooled FST tests are very similar, resulting in largely overlapping curves. We note that this comparison is mainly for illustrative purpose and it may be that different evolutionary scenarios change the behavior of the test statistics. A full color version of this figure is available at the Heredity journal online.

For the long-term success of experimental evolution studies of adaptation from standing variation, it would be helpful to introduce other models that have short generation times, can be cultivated at large effective population sizes and have high recombination rates to uncouple linked sites. Notably, the widely used model organisms yeast and Caenorhabditis elegans are not optimal for this purpose, as recombination is possible only under restricted conditions for these species (see, for example, Parts ; Teotonio ). One possible model, however, may be Caenorhabditis remanei that is obligate sexually reproducing and can be cultivated at large population sizes. Furthermore, natural populations appear to harbor substantial levels of natural variation (Cutter ).

Validating candidates from E&R studies

It is important to distinguish between validating allele frequency estimates obtained from Pool-Seq and validating candidate loci identified in E&R studies. Because Pool-Seq has been shown to obtain reliable allele frequency estimates when some minimum quality criteria are met (Rellstab ; Schlötterer ), we will focus here on the second aspect of validation.

Traits with a simple genetic basis, such as DCV resistance, are best validated by functional analysis of the identified genes and variants. Whether the preferred approach is knockdown of the identified genes by RNA interference, quantitative complementation tests or allelic replacements of candidate SNPs depends on the trait of interest. The validation of candidates for complex traits, however, is a notoriously challenging enterprise as the effect sizes of individual mutations tend to be very small. Confirming the predicted effects of candidate variants by another method such as GWAS provides another feasible strategy (Turner ). Lack of replication does not necessarily indicate lack of an effect, however; for example, if different populations or samples are used for GWAS and E&R, the validation is complicated by allele frequency variation and possible epistatic interactions. Reversing the selection regime may be a particularly appealing validation approach for some E&R experimental designs. Populations that have been selected to drive a trait in one direction could subsequently be selected in the opposite direction, for example, by moving a population from a high to a low temperature regime. Previously, it has been shown that these reverse selection schemes can change the phenotype and allele frequencies at SNP markers in the opposite direction (Teotonio ). Nevertheless, reverse selection will address the problem of linkage between selected and neutral sites only to a moderate extent. One further possibility to validate candidate loci with small effects is via experimental evolution with competing genotypes that differ only in the allele(s) of interest. In these competition assays, continued selection over multiple generations may validate even small functional differences between alleles.

Although gene expression analysis has also been used previously to validate E&R results, we caution that in the absence of a good understanding of how expression differences could affect phenotypes, the interpretation of expression data may be too complex to serve as a stringent validation of candidate genes/SNPs.

Finally, we end this section with a cautionary note on the ability to functionally validate candidate SNPs identified in E&R studies. One implicit assumption is that a larger number of generations will increase the power to detect functionally important loci, but qualitatively similar results are obtained independently of the generation at which the tests are being performed (ignoring variation in selection coefficients among loci). Nevertheless, a recent trajectory analysis found almost entirely different sets of candidates depending on whether generation 15 or 37 was compared with the starting population (Orozco-terWengel ). The reason for this surprising observation is that some alleles increased rapidly early on, but then did not change thereafter, resulting in a frequency plateau between generations 15 and 37. In contrast, many other alleles increased more slowly, but continuously, achieving a higher frequency change by generation 37 than the plateauing alleles. Hence, we suggest that functional validations may benefit from the inclusion of the selection trajectories of candidate loci when comparing E&R results with GWAS or linkage mapping studies.

E&R unlimited?

We anticipate that as our ability to reliably interpret E&R results continues to improve, there will be an increasing number of studies that will apply this method to a substantially broader range of taxa and species. So far, most E&R studies focused on Drosophila, but this approach has also been successfully applied to species with longer generation times and smaller population sizes (Johansson ; Rubin ; Kelly ; Beissinger ). Until now, E&R has been used for rather simple research questions, largely concerned with linking genotypes with phenotypes. But, given the inherent flexibility of experimental evolution framework (Kawecki ), we anticipate a broader use of E&R studies in the future. This could include investigating the impact of migration, different combinations of selective environmental conditions, fluctuating environments or the influence of genetic composition of the starting population, among many others.

E&R studies may be further expanded to study not only the response of the host genome, but also the dynamics of pathogens or endosymbionts during the experiment. One nice example for the potential of this approach comes from the analysis of Wolbachia strains in an experiment that was designed to identify the genomic response of adaptation of D. melanogaster to novel environments (Orozco-terWengel ; Tobler ). After exposing D. melanogaster and its Wolbachia endosymbiont to two different temperature regimes, temperature-dependent differences in the dynamics of Wolbachia strains were uncovered. Although the frequencies of three different Wolbachia clades remained stable in the hot environment, one clade increased from ~25% to ~80% in <15 generations in the cold environment (Versace ).

The E&R approach can also be extended beyond DNA polymorphism to incorporate gene expression levels (Yampolsky ; Hollis ). Contrasting gene expression levels of differentially selected populations provides a powerful and complementary approach to elucidate the selective response. A good example comes from a recent study that compared gene expression response in experimental D. melanogaster populations where either a monogamous or polygamous mating system was enforced over 65 generations (Hollis ). Theory suggests that genes showing sex-biased expression levels due to sexually antagonistic selection—that is, genes affecting traits with different fitness optima between the sexes—should evolve female-like expression levels in a monogamous system, where selection on male traits is relaxed. The experimental results fit this expectation, showing that genes known to have sex-biased expression had feminized their expression in the monogamous, but not the polygamous, population by the end of the study.

Finally, we anticipate that the combination of allele frequency changes with gene expression dynamics will not only assist in the identification of causative variants, but will also provide a more complete picture of the selected trait, in particular when combined with time series analyses.

Data archiving

There were no data to deposit.