Abstract
Experimental evolution systems allow the genomic study of
adaptation, and so far this has been done primarily in asexual
systems with small genomes, such as bacteria and yeast 1–3 . Here
we present whole-genome resequencing data from Drosophila
melanogaster populations that have experienced over 600 genera-
tions of laboratory selection for accelerated development. Flies in
these selected populations develop from egg to adult ,20% faster
than flies of ancestral control populations, and have evolved a
number of other correlated phenotypes. On the basis of 688,520
intermediate-frequency, high-quality single nucleotide poly-
morphisms, we identify several dozen genomic regions that show
strong allele frequency differentiation between a pooled sample of
five replicate populations selected for accelerated development and
pooled controls. On the basis of resequencing data from a single
replicate population with accelerated development, as well as single
nucleotide polymorphism data from individual flies from each
replicate population, we infer little allele frequency differentiation
between replicate populations within a selection treatment.
Signatures of selection are qualitatively different than what has
been observed in asexual species; in our sexual populations,
adaptation is not associated with ‘classic’ sweeps whereby newly
arising, unconditionally advantageous mutations become fixed.
More parsimonious explanations include ‘incomplete’ sweep models,
in which mutations have not had enough time to fix, and ‘soft’ sweep
models, in which selection acts on pre-existing, common genetic
variants. We conclude that, at least for life history characters such
as development time, unconditionally advantageous alleles rarely
arise, are associated with small net fitness gains or cannot fix because
selection coefficients change over time.
Experimental evolution uses well-defined selection protocols to force
phenotypic divergence 4,5 . Studies of experimentally evolved popula-
tions have identified mutations responsible for particular adaptations 6
and provided some general insights into the nature of adaptation in
asexually reproducing populations 7 . Adaptation in these populations is
driven by so-called classic selective sweeps, or the fixation of newly
arising beneficial mutations and the genome-wide haplotypes asso-
ciated with them. By contrast, an obligate sexually reproducing system
can harbour a great deal of standing genetic variation on which selection
can act. Standing variation is theoretically predicted to lead to rapid
evolution in novel environments 8 , and case studies of ecologically rel-
evant genes bear out this prediction 9–11 . The idea that short-term evolu-
tion may act primarily on pre-existing intermediate-frequency genetic
variants that are swept the remainder of the way to fixation has been
termed a soft sweep 8,12 model.
We collected genome-wide resequence data for outbred, sexually
reproducing, replicated populations of D. melanogaster selected for
accelerated development and their matched control populations.
Using the Illumina platform, we obtained short-read sequences from
three genomic DNA libraries: a pooled sample of five replicate popula-
tions that have undergone sustained selection for accelerated develop-
ment and early fertility for over 600 generations (ACO); a pooled sample of five replicate ancestral control populations, which experi-
ence no direct selection on development time (CO); and a single ACO
replicate population (ACO 1 ). The ACO treatment has evolved strongly
differentiated life history phenotypes relative to those of the CO treat-
ment 13 (summarized in Fig. 1; see also Supplementary Fig. 1 for the
history of the populations).
To identify single nucleotide polymorphisms (SNPs) significantly
differentiated between the ACO and CO populations, we aligned reads
to the reference genome of Drosophila and considered only those
genomic positions at which there were two observed allelic states.
After quality-filtering, we were left with 688,520 SNPs: approximately
one SNP for every 175 base pairs (bp) on the 120-megabase (Mb)
euchromatic genome (Methods). The average alignment depth at
identified SNPs was ,203 in both the ACO and CO libraries
(Supplementary Fig. 2), and ,103 in the ACO 1 library. For every
SNP, we calculated 2log 10 (P) from a Fisher’s exact test (L 10 FET) for
a difference in allele frequency between the ACO and CO libraries, as
well as the ACO and ACO 1 libraries.
Links and resources
Tags