Skip to main content
  • Research article
  • Open access
  • Published:

Female responses to experimental removal of sexual selection components in Drosophila melanogaster

Abstract

Background

Despite the common assumption that multiple mating should in general be favored in males, but not in females, to date there is no consensus on the general impact of multiple mating on female fitness. Notably, very little is known about the genetic and physiological features underlying the female response to sexual selection pressures. By combining an experimental evolution approach with genomic techniques, we investigated the effects of single and multiple matings on female fecundity and gene expression. We experimentally manipulated the opportunity for mating in replicate populations of Drosophila melanogaster by removing components of sexual selection, with the aim of testing differences in short term post-mating effects of females evolved under different mating strategies.

Results

We show that monogamous females suffer decreased fecundity, a decrease that was partially recovered by experimentally reversing the selection pressure back to the ancestral state. The post-mating gene expression profiles of monogamous females differ significantly from promiscuous females, involving 9% of the genes tested (approximately 6% of total genes in D. melanogaster). These transcripts are active in several tissues, mainly ovaries, neural tissues and midgut, and are involved in metabolic processes, reproduction and signaling pathways.

Conclusions

Our results demonstrate how the female post-mating response can evolve under different mating systems, and provide novel insights into the genes targeted by sexual selection in females, by identifying a list of candidate genes responsible for the decrease in female fecundity in the absence of promiscuity.

Background

The evolution of mating strategies, and in particular of multiple mating by females (polyandry), has attracted a great deal of attention in recent decades [1]-[3]. The debate on the adaptive significance of female multiple mating stems from the common assumption that males exhibit a stronger positive covariance between promiscuity and reproductive success than females. In other words, males may gain more offspring by repeated matings than females, even though both sexes have the same average numbers of matings, mates and offspring [3]. Polyandry is also assumed to carry costs in terms of time and energy for additional matings [4],[5] or physical injury [6],[7], as well as an increased risk of predation and infection during copulation [7].

Empirical studies however, show that in the vast majority of species, females often mate with more than one male [1],[8],[9]. Theoretically, female polyandry can be promoted by selection if males provide resource benefits, through the ejaculate [10]-[12] or through additional paternal care [13],[14], or if some males do not provide viable sperm or insufficient ejaculate to fertilize the ova [15],[16]. It has also be proposed as a strategy to reduce sexual harassment [17]. Moreover, there could be indirect genetic benefits by acquiring `good genes’, compatible genes or in producing genetically diverse progeny or promoting sperm competition [18]. Finally, multiple mating can be non-adaptive for females in the presence of strong selection for multiple mating on males coupled with a strong intersexual genetic correlation for mating propensity [19].

Recently, experimental evolution studies have been an increasingly popular approach to evaluate the effect of different mating systems on male or female fitness [20], predominantly within the framework of sexually antagonistic co-evolution. Under a promiscuous (multiple) mating system, where the fitness values of an individual and its mate are not perfectly correlated, pre-copulatory and post-copulatory intrasexual competition are expected to result in the evolution of traits that increase the reproductive success of members of one sex at the expense of the other, in a co-evolutionary arms race called interlocus sexual conflict [21]. An eminent example of harm induced by males to females in an attempt to maximize their mating rate and fertilization success is represented by Drosophila melanogaster, in which courtship and transfer of seminal fluid are known to increase female mortality rate and decrease lifetime reproductive success while increasing male competitive abilities [22].

Holland and Rice [23] were the first to investigate the change in female reproductive success in populations of D. melanogaster using an experimental evolution design where sexual selection was removed by enforcing single partner (monogamy) and random mating assignment. They found that monogamous populations had greater net reproductive rate than (promiscuous) controls, while fecundity of monogamous females was reduced after mating with ancestral (promiscuous) males [23]. After this seminal paper, several other studies employed a similar methodology in different taxa [24]-[34], briefly reviewed by Edward et al. [20], with some degree of variation in experimental design and outcome.

Regardless of the adaptive significance of female polyandry, the genetic basis of the fitness components that depend on different mating strategies is a key aspect, which has so far received little attention. In other words, very little information is available about the characteristics and identity of the genes that respond to an alteration of sexual selection but see [35],[36]. With modern genomic techniques, it is possible to scan whole genomes and transcriptomes and associate them with the corresponding phenotypes. Coupling experimental evolution with genome sequencing or transcriptome profiling is a very recent and successful approach [37],[38], in that it experimentally magnifies the variation in the trait of interest and produces a greater resolving power in identifying structure of molecular networks and adaptive processes [39]. However, these methods to our knowledge have not been widely applied to sexual selection studies so far.

Conversely, other aspects of the fruit fly reproductive biology are much better known. In recent years, considerable quantities of data have been collected on the female physiological changes associated with the shifts in female mating status. Molecular and genomics techniques have been employed to investigate the effects of mating in D. melanogaster [40]-[44]. In particular, several detailed studies have focused on seminal fluid components on female post-mating physiology, leading to the identification of several seminal fluid proteins (SFP) and to the isolation of their effect in females reviewed in [45],[46], including the characterization of sex-peptide and its receptor [47],[48].

Here, we integrate these approaches to investigate the evolutionary response of populations experiencing differing sexual selection pressures, at both a phenotypic and genomic level, allowing a direct comparison between the two. We begin by using experimental evolution to evaluate the effects of the removal of components of sexual selection in a laboratory-adapted population of D. melanogaster. The effects of enforced monogamy are then investigated both in terms of differences in female reproductive output and in female post-mating response, measured as genome-wide gene expression profiles. In addition, for populations that have evolved under enforced monogamy we subsequently reverse the selection pressure back to the ancestral promiscuous state and again investigate how reproductive output is affected, demonstrating that our results are unlikely to be due to inbreeding. We take an exploratory approach to investigate the characteristics and biology of those transcripts identified as being influenced by the experimental selection regimes, with the ultimate aim of understanding in more detail which biological processes in females are associated with evolutionary changes in number of matings.

Methods

Fly stocks

All flies used to constitute the experimental evolution lines were derived from a large outbred wild-type population of D. melanogaster (LHM) that had been maintained under the same rearing protocol for over 400 non-overlapping generations for a detailed description, see [49]. The population is maintained in a set of 56 vials at a large size (1792 adults) under competitive conditions and at moderate larval density in standard rearing environment: 25°C, cornmeal/molasses/yeast/agar medium, 12 h:12 h light/dark cycle, 16 individuals of each sex per vial (25 mm × 95 mm) with a 14 day generation cycle. We applied the same culturing condition to our experimental lines, unless otherwise specified.

Experimental evolution lines

An overview of the entire experiment with timings of the various assays is given in Figure 1. In March 2008, a replicate of the ancestral LHM population was obtained by allowing females to lay eggs for 18 h (Day 0). On the day of emergence, Day 10, we collected 384 virgin males and 384 virgin females from the base population, and randomly assigned them to 2 treatments, each constituted by 4 replicate populations of 96 individuals, and stored separately by sex. On Day 13, males and females were placed together in fresh vials (16 pairs per vial, three vials per population) with 6 mg fresh yeast and allowed to mate. In one treatment (hereafter referred to as “monogamous treatment”; M), males were removed after 1 h under brief CO2 anesthesia and discarded. During this window of time, in our LHM population virtually all the sexually mature and healthy females mate, but none of the females mate twice, due to their refractory period. We performed a preliminary study using time-lapse photography see [50] to confirm this pattern: we placed 25 vials containing 10 virgin males and 10 virgin females in an incubator under standard conditions and monitored their activity for 12 h. The results show a peak in mating activity (often 10 pairs simultaneously) between 10 and 30 minutes, followed by a long (>1 h) refractory period, after which females start re-mating (Additional file 1: Figure S8). In the other treatment (hereafter “promiscuous treatment”; P), males were left in the vials with the females and allowed to mate further. On day 14 (i.e. Day 0 of the following generation), the flies were transferred to fresh vials to oviposit for 18 h. The following day (Day 1), eggs were counted and those exceeding 150 were removed to ensure a uniform larval environment. On Day 10, 96 individuals per replicate population were collected as virgin (48 males and 48 females) and the same culturing conditions described above were applied every generation. It should be born in mind then that while frequency of mating differs between the two treatments, they also show differences in adult sex ratio and density following removal of males in the M treatment.

Figure 1
figure 1

Overview of experiment . Experimental populations evolving under monogamous (M, in blue) and promiscuous (P, in red) mating systems were derived from the LHM base population at generation GO, the reverse selection lines (MP, in black) occurred at G95. The intervals between each of the 5 female fecundity trials, as well as body size and microarray assays are also shown.

Body size

After 30 generations of experimental evolution we harvested 40 males and 40 females from each replicate population to assess whether there had been a change in body size, which one may expect if the two selection regimes experienced different levels of inbreeding. A single wing was removed from each individual, mounted on a slide using transparent tape and photographed using bright field-illumination (×40 magnification). Length was measured using the straight line tool in ImageJ [51], from the intersection of the anterior cross vein and longitudinal vein 3 (L3) to the intersection of L3 with the distal wing margin [52].

Female fecundity

We obtained measures of fecundity of females in our experimental lines under four different conditions: after a single mating and after being continuously exposed to males, during the whole period in which they were allowed to oviposit in our selection regime, or during a longer timeframe, to account for potential shifts in the resource that females allocate to eggs over time.

The effect of the treatment on female fecundity was assayed with a factorial design in four different trials, after 30, 31, 50 and 58 generations of experimental evolution, with slightly different experimental designs (trials 1–4 respectively, see below). For each trial, the following protocol was applied: on day 14 of the chosen generation, a replicate of the experimental lines were obtained, by allowing flies to oviposit for an additional 24 h in fresh vials. The populations obtained were cultured with standard protocol (i.e. as in the promiscuous treatment) for a generation to remove the majority of any parental effects. On Day 10 of the following generation, 160 females and 160 males for each treatment and replicate were collected as virgins and stored separately (10 vials of males and 10 vials of females for each of the 8 experimental lines). On Day 13, half of the females from each experimental line (5 vials) were crossed to males from the same experimental line, and the other half (5 vials) were crossed to males from a single replicate population of the other treatment (females from replicate 1 of the monogamous treatment were mated to males from replicate 1 of the promiscuous treatment, replicate 2 of the monogamous treatment was paired with replicate 2 of the promiscuous treatment, and so on), and allowed to mate in fresh vials containing 6 mg live yeast. At this stage, the trials differed in their design, as follows.

In trial 1, males were removed and discarded the following day (after 30 h, Day 14), while 14 individual females were transferred in oviposition test tubes, and allowed to oviposit for 18 hours, corresponding to the window of time in which eggs laid by females in the experimental populations were retained for the next generation. Females were then discarded, the tubes refrigerated for 24 hours and the eggs counted. In trial 2, the protocol employed was identical to the one described for trial 1, except the males were removed and discarded after 1 h, allowing females to mate only once.

In trial 3, after 1 h, all the males were removed and discarded. Groups of 16 females were allowed to oviposit, and were transferred every 12 h (at 9:00 and at 21:00) to a fresh vial for four days (6 times, 7 time-points) to avoid excessive larval density, then discarded. When the new generation emerged, progeny were counted. In trial 4, the protocol employed was very similar to the one described for trial 3, with the following differences: after crossing target males and females, males were not removed from the vials; also, during the four days of oviposition the flies were transferred every 6 h during the daylight hours (9:00, 15:00 and 21:00; 9 times, 10 time-points).

Reversed selection lines

After 95 generations of experimental evolution a third treatment was established using surplus flies harvested from each of the 4 monogamous populations. In this new treatment, the rearing protocol was identical to that for the promiscuous treatment, and therefore flies in these populations experienced a reversal of the selection pressure from a monogamous to a promiscuous mating system (hereafter refereed to as the MP treatment). After a further 25 generations of experimental evolution (generation 120 in total) all populations were cultured with standard protocol for a single generation to remove parental effects, then an assay of female fecundity from all replicate populations and treatments was performed (trial 5; n = 53-64 individual females per population), using the same protocol employed for trial 1 (see above).

Microarray data

After 46 generations, on day 14, replicates of the experimental lines were obtained, by allowing flies to oviposit for additional 24 h in fresh vials. The populations obtained were cultured with standard protocol (i.e. as the promiscuous treatment) for a generation to remove parental effects. On Day 10 of the second generation, 64 females and 64 males for each treatment and replicate were collected as virgins and stored separately (4 vials of males and 4 vials of females for each of the 8 experimental lines). On Day 13, half of the females from each experimental line (2 vials) were crossed to males from the same experimental line, and the other half (2 vials) were crossed to males from a single replicate population of the other treatment, and allowed to mate in fresh vials containing 6 mg live yeast. After 1 h, all the males were removed and discarded, while the females were randomly divided in two groups of 8 flies under brief CO2 anesthesia, to be used as a main sample and its backup. A single mating treatment was used in order to allow a more direct comparison with previously published data (see below). After 6 hours, the females (whole body) were flash-frozen in liquid nitrogen and stored at −80°C for no more than four days until RNA extraction. Hence, for each of the replicate population we collected 4 independent samples of eight females, two samples of females mated to males of the same replicate population, and two samples of females mated to males of a single replicate population in the other treatment, giving a total of 32 samples. Total RNA was extracted independently from each sample using Trizol (Invitrogen, Carlsbad, CA, USA) and purified with an RNeasy Mini Kit (Qiagen, Hilden, Germany). RNA quality and quantity was assessed with an Agilent Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The RNA samples were prepared and hybridized to Affymetrix Drosophila GeneChip 2.0 microarrays at the Uppsala Array Platform (Uppsala, Sweden) following manufacturer’s instructions. The arrays were scanned in two batches of 16, balanced for replicate population of origin and replicate population of origin of the males to which they were mated [GEO:GSE48385].

Statistical analysis

All statistical analyses were run in the R environment version 2.11.1 for most analyses, version 3.0.1 for body size and reversed experimental evolution assay, available at http://www.r-project.org [53].

Male and female body size was analyzed using a full factorial linear model (lm function; mating system and sex as fixed effects; no interaction term) using within replicate means to avoid pseudoreplication (n = 36).

Female fecundity data for each trial were analyzed using linear models (lm function). In all cases, amount of eggs or progeny produced was averaged within replicate population and summed across time points (for trials 3 and 4), to avoid pseudoreplication. We fitted the following model to each dataset:

y i j k = f i + m j + I i j + e i j k

with i = {1,2}; j = {1,2}; k = {1,…,4}; where y is the number of progeny/eggs produced by females after each cross, f is the treatment of origin of the females (monogamous or promiscuous, fixed effect), m is the treatment of origin of the males to which females were mated (monogamous or promiscuous, fixed effect) and I is their interaction. The interaction term was subsequently dropped, because it was not significant in any trial and did not improve the fit of the models (P > 0.25 for all models).

Microarray data were analysed using the BioConductor suite of packages [54] in R. To pre-process the raw expression data, we used the standard RMA (Robust Multichip Average) algorithm [55] implemented in the affy package [56]. After pre-processing the resulting dataset was filtered to exclude features according the following criteria: (i) probe sets without an Entrez Gene ID annotation, (ii) Affymetrix quality control probe sets, (iii) if multiple probe sets mapped to the same Entrez Gene ID, only the probe set with the highest coefficient of variation was retained. Out of the original 18952 features, the filtering step removed 6380 probe sets, while 12572 probe sets, corresponding to as many known genes, were retained for the statistical analyses.

Significance of differential expression was assessed using the package limma Linear Models for Microarray Data; [57]. A model matrix was designed to fit a parameter for every combination of replicate population of origin of females (n = 8) and population of origin of males to which females were mated (n = 8), for a total of 16 parameters. An additional random effect with two levels was fitted to control for the batch effect, and estimated borrowing information between features, by constraining the within-block correlations to be equal across features and by using empirical Bayes methods to moderate the standard deviations [57]. A contrast matrix was designed to obtain the contrasts of interest: the main effect of treatment of origin of females, the main effect of treatment of origin of males to which females were mated, and their interaction. All the resulting P values were corrected for multiple testing to obtain a maximum false discovery rate of 5% FDR; [58]; corrected P <0.05.

We used a mean-rank gene set enrichment test MR-GSE test, implemented in LIMMA; [59] to test whether the sets of up-regulated or down-regulated significant transcripts showed a tendency to be up- or down-regulated after mating, using the t-values from a contrast between virgin and singly mated females from the same population, from a previously published study [44]. Genes in these sets that were found to be up-regulated by mating were classified as `Virgin-like’ if they were also lower in monogamous females relative to promiscuous females, or if they were down-regulated by mating and higher in monogamous females relative to promiscuous ones. We therefore defined as `Virgin-like’ the subset of transcripts for which the expression in the monogamous female is lower than the promiscuous females if the gene is up-regulated by mating, or higher if the gene is down-regulated by mating; in other words, genes whose profile is more similar to a virgin fly. `Mated-like’ genes are the complimentary set of genes. A MR-GSE test was also used to test whether the set of significant transcripts showed a tendency to be associated with female fitness, using the t-values from a previously published study [60] on the same population.

Among the genes found to be differentially expressed, we identified transcriptional modules of correlated expression across-tissues using the hopach package Hierarchical Ordererd Partitioning and Collapsing Hybrid; [61]. We computed a distance matrix using the pairwise correlations r ij between the expression of the significant transcripts across different tissues of D. melanogaster. The tissue-specific expression data were produced by the FlyAtlas team [62], available on the Gene Expression Omnibus [GEO:GSE7763], and normalized according to a method described elsewhere [60]. The clustering algorithm built a hierarchical tree by recursively partitioning or collapsing clusters at each level, using MSS (Median Split Silhouette) criteria to identify the level of the tree with maximally homogeneous clusters [61].

We selected the modules containing more than 50 genes (the clusters more likely to provide biologically meaningful within-module summary statistics) and analyzed them to identify whether they showed: (i) association with genes involved in female post-mating response (see above); (ii) association with female fitness (see above); (iii) non random chromosomal distribution (based on genes present on the microarray); (iv) over-represented Gene Ontology categories; (v) tissue enrichment or specificity. Non-random chromosomal distribution was assessed with a Fisher’s exact test on the expected and observed number of genes on each chromosome (P < 0.01). To identify GO categories enriched for particular subsets of transcripts, we used a hypergeometric test for over-representation P < 0.01, GOSTATS package; [63]. Allometric differences in tissue size between treatments may contribute to patterns of differential expression in whole body extractions. By exploring tissue enrichment or specificity within the set of differentially expressed genes we are therefore able to provide a list of candidate tissues that may have been influenced by the mating system treatment. We identified tissue-enriched or tissue-specific transcripts using data from the FlyAtlas database (described above; 2-fold cut-off). The tissue specific expression levels for the list of transcripts in each module were obtained, and the modules were tested for over-abundance of genes of interest in a target tissue using a one-tailed Fisher’s exact test. All the reported P values were Bonferroni corrected for testing on multiple tissues (P < 0.01, n = 17).

Results

Body size

After 30 generations of experimental evolution the body sizes of individuals from the two selection regimes remained virtually unchanged, with males and females of promiscuous lines having approximately 1% smaller wings than those in monogamous lines (mating system effect: F 3,12 = 1.84, P = 0.200). There is no evidence that males evolved to become larger in the promiscuous lines relative to the ancestral population see Additional file 1: Figures S2 and S3 in [64].

Female fecundity

The reproductive output of promiscuous females was greater than those of monogamous females, regardless of the males they were mated with (Figure 2; see Methods). This difference was significant when measured as number of eggs laid in a 18 h period, corresponding to the oviposition period in every generation of experimental evolution (as well as the ancestral population), both after a single mating with a male (F 1,13 = 4.89, P = 0.046), or being continuously exposed to males (F 1,13 = 11.32, P = 0.005). When measured as number of adult progeny emerging from eggs laid during a period of 4 days, this difference was significant only after a single mating (F 1,13 = 6.06, P = 0.029), again regardless of the male partner, but not when females where continuously exposed to males (F 1,13 = 2.93, P = 0.110), although the effect sizes were comparable for direction and magnitude (Table 1).

Figure 2
figure 2

Female fecundity . Reproductive output of females evolved under monogamous and promiscuous selection regimes after mating once (panels B and C) or being continuously exposed to males (panels A and D), during the normal reproductive window (panel A and B; 14 individual females) or during a longer interval (4 days; panel C and D; 16 individual females). Closed circles represent median values, boxes limited to interquartile range, and whiskers at minimum and maximum values.

Table 1 Linear model results on female fecundity

Reversed selection lines

As described in the Methods, the reversed selection lines (MP) were established after the monogamous and promiscuous populations had already undergone 95 generations of selection. All three treatments were then run for a further 25 generations prior to the final assays of fecundity being performed (Figure 1, trail 5). Despite this substantial additional period of experimental evolution, relative to the first round of assays (a further 90 generations), the patterns of reproductive output between monogamous and promiscuous females was remarkably similar at these two time points (M vs P % differences: G30 45%; G120 36%), indicating that the majority of phenotypic evolution had occurred within the first 30 generations (posthoc Tukey HSD: M vs P, t =2.716, P = 0.0279; Figure 3). The phenotypic change in reproductive output of females from the MP populations following 25 generations of reversed selection was smaller and non-significant. It did however indicate a reversal in reproductive output had occurred; posthoc tests showed that the reproductive output of females from the MP treatment was intermediate to both monogamous and promiscuous treatments (Tukey HSD: M vs MP, t = 0.913, P = 0.6365; MP vs P, t =1.804, P = 0.1850; Figure 3).

Figure 3
figure 3

Reversed selection . Individual level reproductive output of females evolved under monogamous (M), monogamous then promiscuous (MP) and promiscuous (P) selection regimes. Results of the posthoc analysis are given above the plotting frame, letters not shared indicate treatments that show statistical significant differences (see Results for details). Closed circles represent median values, boxes limited to interquartile range, and whiskers at minimum and maximum values.

Gene expression profiles

After 46 generations of selection, we tested the difference in female genome-wide post-mating response, by measuring gene expression in adult D. melanogaster females evolved under different sexual selection regimes (monogamous and promiscuous). After multiple testing correction, females from the monogamous and promiscuous populations showed a significant difference in the expression of 1141 transcripts (≈9% of the transcripts analyzed, at 5% F.D.R.), while male type and the interaction of female type and male type did not significantly affect the post-mating expression patterns. Among the differentially expressed transcripts, 438 were up-regulated and 703 down-regulated in monogamous females (Binomial test: ratio = 0.38, P < 0.0001).

We compared the expression profile of these transcripts with the female post-mating response characteristic to the ancestral population [44], and found that the expression level of 728 genes is altered in monogamous females to a lesser extent after mating, compared to promiscuous females (hereafter `virgin-like’, see Methods), while the post-mating reaction of the remaining 413 genes is altered to a higher extent in monogamous females compared to promiscuous females (hereafter `mated-like’), and their proportion is higher than expected by chance (Binomial test: ratio = 0.64, P < 0.0001). In general, genes that are down-regulated in monogamous vs. promiscuous mated females tend to be switched on (up-regulated) by mating (two-tailed MR-GSE test: P < 0.0001, Figure 4A) and genes that are up-regulated in monogamous vs. promiscuous mated females tend to be switched off (down-regulated) by mating (two-tailed MR-GSE test: P < 0.0001, Figure 4A). We also used previously published and independently derived data [60] to test the relationship between the significant transcripts identified in this study and female fecundity, and found them to be over-represented among genes strongly associated with female fitness (irrespective of the sign of the association; MR-GSE test, P < 0.0001, Figure 4B).

Figure 4
figure 4

Association with post-mating response and female fitness . (A) Density distribution of significant up-regulated (blue) and down-regulated (red) transcripts in monogamous (versus promiscuous) females along all the tested genes, ranked according to their post-mating reaction (data from a previously published study on the same population; Innocenti and Morrow [44]). (B) Density distribution of the significant transcripts along all the tested genes, ranked by the t-value of their association with female fitness (data from a previously published study on the same population; Innocenti and Morrow [60]).

In order to identify clusters of transcripts co-expressed in one or more tissues, and hence possibly involved in similar biological function, we calculated modules of correlated expression among the significant transcripts using data from the FlyAtlas database [62]. Among them, we selected the 7 clusters containing more than 50 genes (Figure 5), which represented about 75% of the significant transcripts, and evaluated their post-mating expression profile in comparison to the ancestral population, their tissue specificity, chromosomal distribution and over-presentation among Gene Ontology categories (see Additional file 1).

Figure 5
figure 5

Transcriptional modules . Level-plot representing the matrix of pair-wise correlation for the expression of the 1141 significant transcripts across tissues of D. melanogaster data from [62]. The correlation matrix has been used to compute modules of correlated expression (separated by grey lines). The 7 modules containing more than 50 genes are labeled.

Module 1 contains genes highly specific for the male gonads, showing little or no expression in other tissues (Additional file 1: Figures S1D,E). Overall, they do not tend to be perturbed by mating (Additional file 1: Figures S1A). Over-represented Gene Ontology (GO) terms indicate that the activity of a portion of these genes is linked to mitochondrial cellular components (cellular respiration, electron transport, Additional file 1: Table S1).

Module 2 is a large cluster of genes active in the majority of the tissue types (but generally not in the gonads, Additional file 1: Figures S2E), and significantly over-expressed in the head, eyes, carcass, fat body, heart and spermatheca. These transcripts are subject to changes in expression levels after mating (Table 2), with mated monogamous flies showing a more virgin-like expression profile for these genes (Additional file 1: Figures S2A). They are chiefly involved in enzymatic metabolic activity (oxidation reduction, proteolysis, Additional file 1: Table S2). The left arm of chromosome 2 is enriched for this set of genes (Fisher exact test: Odds-ratio = 1.49, P = 0.005).

Table 2 Summary description of the main modules

The activity of genes clustered in module 3 is very similar to those of module 2: these transcripts are active ubiquitously in the fruit fly tissues (Additional file 1: Figures S3D,E) and the up-regulated subset is significantly enriched among the set of genes which respond to mating (Table 2). Although not significant under our cut-off, a higher than expected proportion of these genes lies on chromosome 2 L (Odds-ratio = 1.47, P = 0.045). GO terms associated with these genes include, again, strong cytoplasmic enzymatic activity (oxidation reduction, catalytic activity, Additional file 1: Table S3).

Module 4 presents the most distinctive and peculiar patterns. The majority of these genes are down-regulated in the monogamous treatment (91 out of 124, Table 2) and tend to be strongly affected by mating and distinctly more virgin-like in monogamous females (Additional file 1: Figures S4A). These transcripts are consistently highly expressed in the midgut, but relatively silent in all the other tissues (Additional file 1: Figures S4D,E), and most of their activity is linked to metabolic processes, mainly peptidase and hydrolase activity (Additional file 1: Table S4). The distribution on the chromosomes is significantly skewed towards the right arm of chromosome 2 (Odds-ratio = 1.49, P = 0.009).

Modules 5 and 6 show highest relative expression levels in the ovaries, although the transcripts are also active at slightly lower levels in every other tissue. These genes tend to be overall weakly down-regulated after mating (Additional file 1: Figures S5A and S6A). Module 5 showed relative virgin-like expression in monogamous females compared to promiscuous females (Table 2). Perhaps unsurprisingly, sexual reproduction and female gamete generation were among the most enriched biological processes, while the same sets of genes were linked to nucleic acid and protein binding molecular functions (Additional file 1: Table S5). Module 6, also significantly over-expressed in the brain, showed enrichment for biological processes such as behaviour and signaling processes (Additional file 1: Table S6).

Module 7 contains genes significantly more active in neural tissues: brain, thoracic ganglion, head and eyes (Additional file 1: Figures S7D,E). They are mostly down-regulated (Table 2) in monogamous females, but tend to be up-regulated after mating (virgin-like in monogamous females, Additional file 1: Figures S7A).

All modules tend to be associated with female fitness, with transcripts in modules 2, 3 and 4 showing the strongest association (Table 2).

Discussion

In our study, we experimentally manipulated the mating system in replicate populations of D. melanogaster, by removing components of post-copulatory sexual selection, with the aim of testing differences in short term post-mating reaction of females evolved under different mating strategies. We showed that monogamous females suffer decreased fecundity, regardless of the type of male they were mated with, or whether mated once or continuously exposed to males. We also showed that monogamous females could recover some of this loss in fecundity if the selection pressure was reversed experimentally. Previously, Holland and Rice [23] removed sexual selection in experimental lines from the same population (LHM) by manipulating sex ratio, and found that (i) monogamous females showed higher `net reproductive rate’ (female fecundity and offspring survival) than controls when mated with males from their own populations, and (ii) monogamous females showed lower fecundity than controls when mated once to ancestral (promiscuous) males [23]. Monogamous males, in turn, evolved decreased courtship rate. Our experiment employed a different design, which allowed mass mating (mate choice and pre-copulatory intra-sexual competition) but a single mating event in the monogamous treatment, in order to leave selection on courtship rate unaffected. Our results showed no effect of male type on female fecundity (and consequently no interaction between male and female type), which suggests there is no evidence that males evolved decreased courtship intensity or a less harmful ejaculate. It is thus unlikely that the decrease in fecundity of monogamous females reflects a selective pressure towards less `resistant’ females. On the other hand, the experimental treatment removed continuous male harassment in the monogamous environment and decreased population density during selection. Relaxed selection on resistance to male harassment may have allowed the accumulation of deleterious mutations or recombination of extant genetic variation with sub-optimal epistatic effects which resulted in overall decrease in mean female fitness in the monogamous environment; suggesting a potential indirect benefit to females of multiple mating. The decline in female fitness in monogamous lines is not likely to be due to simple differences in population size and subsequent inbreeding, since not only does theoretical and previous empirical work indicate that n = 96 is above the threshold for drift decay see [65], our reversed experimental evolution treatment showed that the significant differences in fecundity between monogamous and promiscuous females disappeared when monogamous females experienced a reintroduction of a promiscuous mating system. Such a response would not occur if monogamous populations had simply become bottlenecked. This is further supported by the minimal differences in body size seen between monogamous and promiscuous treatments, a trait that is sensitive to inbreeding [66].

The results of our genome-wide expression analysis confirmed a significant difference between post-mating reaction between monogamous and promiscuous females, while the evolutionary history of the males to which they were mated did not influence their expression profiles. The genes that evolved to respond differently to mating accounted for around 9% of the transcriptome tested. When comparing transcriptional changes which occur when a female switches between the virgin and mated status with differential expression between mated monogamous and promiscuous females, it is clear that genes which are up-regulated by mating tend to be down-regulated in monogamous vs. promiscuous females (and those down-regulated by mating tend to be up-regulated in monogamous vs promiscuous females; Figure 4A), i.e. monogamous females generally show a more `virgin-like’ expression profile. Similarly, significant genes tend to be over-represented among candidate genes known to be associated with female fitness. Taken together, these two lines of evidence can be interpreted as strong support for the phenotypic results showing decreased female fecundity: monogamous females seem to exhibit a weaker post-mating reaction, both in terms of the extent to which genes are expressed and how many eggs they lay.

When analysing and partitioning these genes according to the tissue where they are predominantly active, we can identify 4 broad categories: transcripts active in i) the midgut, ii) the ovaries, iii) neural tissues and iv) a wide range of tissues. The midgut (module 3, Table 2) provides the strongest and clearest signal, (Additional file 1: Figures S4E), and genes active in this tissue are mainly linked to enzymatic activity (Additional file 1: Table S4). Such genes are usually activated by mating and show a decreased response in monogamous females (Additional file 1: Figures S4A). Significant genes in the ovaries (module 5 and 6) are involved in gamete production, while in the neural tissues they regulate signaling processes and transmembrane transport activity (Additional file 1: Table S5, S6). The last category (module 2 and 3) contains genes expressed in a diverse array of tissues (Table 2), although known to show overall very high correlation for expression [60]. The Gene Ontology categories involved (oxidoreductase activity, lipid and sugars storage/metabolism, Additional file 1: Table S2, S3) seem to indicate a predominant function in energy production and resource consumption. An additional, small category points to transcripts mostly active in the testes, and its interpretation is problematic, given the sex-limited nature of this tissue. This set of genes, however, which are not involved in a normal post-mating reaction (Additional file 1: Figures S1A), could be selected due to pleiotropic activity in other tissues, or exhibit non-random patterns of coexpression with selected transcripts.

More recently, two studies in D. melanogaster [67] and D. pseudoobscura [36] have also utilized experimental evolution to investigate the effects of enforced monogamy (while allowing multiple mating) on expression profiles. Hollis et al. [67] examined expression profiles from unmated male and female flies and found a general pattern of feminization of the sex-biased genes in monogamous populations, consistent with predictions from sexual conflict theory that sex-specific selection may favour and maintain patterns of sex-specific trait expression. Similar processes are likely to have occurred in our populations although a direct comparison of the results is difficult since the factors explored in the two overlap only partially. Nonetheless it was noted that feminization extends beyond sex-limited tissues. Immonen et al. [36] enforced monogamy for 100 generations in replicate populations of D. pseudoobscura and also examined expression profiles from unmated females. They found that females evolving under enforced monogamy may have elevated investment in somatic rather than reproductive processes, but in contrast to Hollis et al. expression profiles of those females were less feminized, which is likely due to qualitative differences between the two species in the costs of mating for females. Taking these three studies together, it is clear that we still have much to learn before we have a complete picture of how sexual selection pressures influence the transcriptomes of males and females throughout their reproductive lives.

Conclusions

This study, in combination with the independent characterization of post-mating expression profiles in females and the relationship between transcript abundance and female fitness in the ancestral population, provides a list of genes associated with changes in female fecundity caused by evolution under different sexual selective pressures. Given the general agreement between what is already known about the effects of the male ejaculate on females and their fitness in D. melanogaster, with the patterns of tissue specificity and biological processes identified here, these data provide an indication as to which genes may be targeted by the process of post-mating sexual selection in this promiscuously mating population.

Availability of supporting data

The data sets supporting the results of this article are available in the Dryad repository, doi:http://dx.doi.org/10.5061/dryad.2dp4q [68] and the Gene Expression Omnibus [GSE48385]: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE48385 [69].

Authors’ contributions

PI and EM conceived of the study, carried out the forward experimental evolution and drafted the manuscript. IF carried out the reverse selection and time lapse of mating behaviour PI carried out the fitness assays and the microarray experiment. All authors read and approved the final manuscript.

Additional file

Abbreviations

F.D.R.:

False discovery rate

GO:

Gene ontology

HSD:

Honest significant difference

LHM:

Base population

M:

Population evolved under a monogamous mating system

MP:

Population derived from M but then evolved under a promiscuous mating system

MR-GSE:

Mean-rank gene set enrichment

MSS:

Median split silhouette

P:

Population evolved under a promiscuous mating system

RNA:

Ribosomal nucleic acid

SFP:

Seminal fluid proteins

References

  1. Andersson M: Sexual Selection. 1994, Princeton University Press, Princeton

    Google Scholar 

  2. Birkhead T: Promiscuity: An Evolutionary History of Sperm Competition. 2000, Harvard University Press, Faber and Faber, London

    Google Scholar 

  3. Shuster SM, Wade MJ: Mating Systems and Strategies. 2003, Princeton University Press, Princeton

    Google Scholar 

  4. Jormalainen V, Merilaita S, Riihimäki J: Costs of intersexual conflict in the isopod Idotea baltica. J Evol Biol. 2001, 14: 763-772. 10.1046/j.1420-9101.2001.00325.x.

    Article  Google Scholar 

  5. Wedell N, Kvarnemo C, Lessells CM, Tregenza T: Sexual conflict and life histories. Anim Behav. 2006, 71: 999-1011. 10.1016/j.anbehav.2005.06.023.

    Article  Google Scholar 

  6. Boeuf BJL, Mesnick S: Sexual behavior of male Northern Elephant Seals: I. Lethal injuries to adult females. Behav. 1991, 116: 143-162. 10.1163/156853990X00400.

    Article  Google Scholar 

  7. Hurst GD, Sharpe RG, Broomfield AH, Walker LE, Majerus TM, Zakharov IA, Majerus ME: Sexually transmitted disease in a promiscuous insect, Adalia bipunctata. Ecol Entomol. 1995, 20: 230-236. 10.1111/j.1365-2311.1995.tb00452.x.

    Article  Google Scholar 

  8. Ridley M: Mating frequency and fecundity in insects. Biol Rev. 1988, 63: 509-549. 10.1111/j.1469-185X.1988.tb00669.x.

    Article  Google Scholar 

  9. Simmons LW: Sperm Competition and Its Evolutionary Consequences in the Insects. 2001, Princeton University Press, Princeton

    Google Scholar 

  10. Thornhill R: Sexual selection and nuptial feeding behavior in Bittacus apicalis (Insecta: Mecoptera). Am Nat. 1976, 110: 529-548. 10.1086/283089.

    Article  Google Scholar 

  11. Boggs CL, Gilbert LE: Male contribution to egg production in butterflies: evidence for transfer of nutrients at mating. Science. 1979, 206: 83-84. 10.1126/science.206.4414.83.

    Article  PubMed  CAS  Google Scholar 

  12. Turner ME, Anderson WW: Multiple mating and female fitness in Drosophila pseudoobscura. Evolution. 1983, 37: 714-723. 10.2307/2407913.

    Article  Google Scholar 

  13. Stacey PB: Female promiscuity and male reproductive success in social birds and mammals. Am Nat. 1982, 120: 51-64. 10.1086/283969.

    Article  Google Scholar 

  14. Davies NB: Cooperation and conflict among dunnocks, Prunella modularis, in a variable mating system. Anim Behav. 1985, 33: 628-648. 10.1016/S0003-3472(85)80087-7.

    Article  Google Scholar 

  15. Gibson R, Jewell P: Semen quality, female choice and multiple mating in domestic sheep - a test of Trivers sexual competence hypothesis. Behav. 1982, 80: 9-31. 10.1163/156853982X00418.

    Article  Google Scholar 

  16. Gromko MH, Newport ME, Kortier MG: Sperm dependence of female receptivity to remating in Drosophila melanogaster. Evolution. 1984, 38: 1273-1282. 10.2307/2408634.

    Article  Google Scholar 

  17. Svärd L, Wiklund C: Different ejaculate delivery strategies in 1st versus subsequent matings in the Swallowtail Butterfly Papilio machaon. Behav Ecol Sociobiol. 1986, 18: 325-330. 10.1007/BF00299663.

    Article  Google Scholar 

  18. Yasui Y: The `genetic benefits’ of female multiple mating reconsidered. Trends Ecol Evol. 1998, 13: 246-250. 10.1016/S0169-5347(98)01383-4.

    Article  PubMed  CAS  Google Scholar 

  19. Halliday T, Arnold SJ: Multiple mating by females: a perspective from quantitative genetics. Anim Behav. 1987, 35: 939-941. 10.1016/S0003-3472(87)80138-0.

    Article  Google Scholar 

  20. Edward DA, Fricke C, Chapman T: Adaptations to sexual selection and sexual conflict: insights from experimental evolution and artificial selection. Phil Trans R Soc B. 2010, 365: 2541-2548. 10.1098/rstb.2010.0027.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Parker GA: Sexual Selection and Sexual Conflict. Sex Sel Reprod Compet Insects. Edited by: Blum MS, Blum NA. 1979, Academic Press, London, 123-166. 10.1016/B978-0-12-108750-0.50010-0.

    Google Scholar 

  22. Chapman T, Liddle LF, Kalb JM, Wolfner MF, Partridge L: Cost of mating in Drosophila melanogaster females is mediated by male accessory gland products. Nature. 1995, 373: 241-244. 10.1038/373241a0.

    Article  PubMed  CAS  Google Scholar 

  23. Holland B, Rice WR: Experimental removal of sexual selection reverses intersexual antagonistic coevolution and removes a reproductive load. Proc Natl Acad Sci U S A. 1999, 96: 5083-5088. 10.1073/pnas.96.9.5083.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  24. Hosken D, Garner TW, Ward P: Sexual conflict selects for male and female reproductive characters. Curr Biol. 2001, 11: 489-493. 10.1016/S0960-9822(01)00146-4.

    Article  PubMed  CAS  Google Scholar 

  25. Hosken DJ, Martin OY, Wigby S, Chapman T, Hodgson DJ: Sexual conflict and reproductive isolation in flies. Biol Lett. 2009, 5: 697-699. 10.1098/rsbl.2009.0066.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. Martin O, Hosken D: Costs and benefits of evolving under experimentally enforced polyandry or monogamy. Evolution. 2003, 57: 2765-2772. 10.1111/j.0014-3820.2003.tb01518.x.

    Article  PubMed  Google Scholar 

  27. Crudgington HS, Beckerman AP, Brüstle L, Green K, Snook RR: Experimental removal and elevation of sexual selection: does sexual selection generate manipulative males and resistant females?. Am Nat. 2005, 165: S72-S87. 10.1086/429353.

    Article  PubMed  Google Scholar 

  28. Tilszer M, Antoszczyk K, Ssaek N, Zajaoc E, Radwan J: Evolution under relaxed sexual conflict in the Bulb Mite Rhizoglyphus robini. Evolution. 2006, 60: 1868-1873. 10.1111/j.0014-3820.2006.tb00530.x.

    Article  PubMed  Google Scholar 

  29. Bacigalupe LD, Crudgington HS, Hunter F, Moore AJ, Snook RR: Sexual conflict does not drive reproductive isolation in experimental populations of Drosophila pseudoobscura. J Evol Biol. 2007, 20: 1763-1771. 10.1111/j.1420-9101.2007.01389.x.

    Article  PubMed  CAS  Google Scholar 

  30. Fricke C, Arnqvist G: Rapid adaptation to a novel host in a Seed Beetle (Callosobrucus maculatus): the role of sexual selection. Evolution. 2007, 61: 440-454. 10.1111/j.1558-5646.2007.00038.x.

    Article  PubMed  Google Scholar 

  31. LaMunyon CW, Bouban O, Cutter AD: Postcopulatory sexual selection reduces genetic diversity in experimental populations of Caenorhabditis elegans. J Hered. 2007, 98: 67-72. 10.1093/jhered/esl052.

    Article  PubMed  CAS  Google Scholar 

  32. Simmons LW, Garcia Gonzalez F: Evolutionary reduction in testes size and competitive fertilization success in response to the experimental removal of sexual selection in Dung Beetles. Evolution. 2008, 62: 2580-2591. 10.1111/j.1558-5646.2008.00479.x.

    Article  PubMed  Google Scholar 

  33. Gay L, Eady PE, Vasudev R, Hosken DJ, Tregenza T: Does reproductive isolation evolve faster in larger populations via sexually antagonistic coevolution?. Biol Lett. 2009, 5: 693-696. 10.1098/rsbl.2009.0072.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Maklakov A, Bonduriansky R, Brooks R: Sex differences, sexual selection, and ageing: and experimental evolution approach. Evolution. 2009, 63: 2491-2503. 10.1111/j.1558-5646.2009.00750.x.

    Article  PubMed  Google Scholar 

  35. Immonen E, Ritchie MG: The genomic response to courtship song stimulation in female Drosophila melanogaster. Proc R Soc B Biol Sci. 2011, 279: 1359-1365. 10.1098/rspb.2011.1644.

    Article  Google Scholar 

  36. Immonen E, Snook RR, Ritchie MG: Mating system variation drives rapid evolution of the female transcriptome in Drosophila pseudoobscura. Ecol Evol. 2014, 4: 2186-2201.

    PubMed  PubMed Central  Google Scholar 

  37. Burke MK, Dunham JP, Shahrestani P, Thornton KR, Rose MR, Long AD: Genome-wide analysis of a long-term evolution experiment with Drosophila. Nature. 2010, 467: 587-590. 10.1038/nature09352.

    Article  PubMed  CAS  Google Scholar 

  38. Zhou D, Udpa N, Gersten M, Visk DW, Bashir A, Xue J, Frazer KA, Posakony JW, Subramaniam S, Bafna V, Haddad GG: Experimental selection of hypoxia-tolerant Drosophila melanogaster. Proc Natl Acad Sci U S A. 2011, 108: 2349-2354. 10.1073/pnas.1010643108.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  39. Turner TL, Stewart AD, Fields AT, Rice WR, Tarone AM: Population-based resequencing of experimentally evolved populations reveals the genetic basis of body size variation in Drosophila melanogaster. PLoS Genet. 2011, 7: e1001336-10.1371/journal.pgen.1001336.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  40. Lawniczak MK, Begun DJ: A genome-wide analysis of courting and mating responses in Drosophila melanogaster females. Genome. 2004, 47: 900-910. 10.1139/g04-050.

    Article  PubMed  CAS  Google Scholar 

  41. McGraw LA, Gibson G, Clark AG, Wolfner MF: Genes regulated by mating, sperm, or seminal proteins in mated female Drosophila melanogaster. Curr Biol. 2004, 14: 1509-1514. 10.1016/j.cub.2004.08.028.

    Article  PubMed  CAS  Google Scholar 

  42. McGraw LA, Clark AG, Wolfner MF: Post-mating gene expression profiles of female Drosophila melanogaster in response to time and to four male accessory gland proteins. Genetics. 2008, 179: 1395-1408. 10.1534/genetics.108.086934.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  43. Mack PD, Kapelnikov A, Heifetz Y, Bender M: Mating-responsive genes in reproductive tissues of female Drosophila melanogaster. Proc Natl Acad Sci U S A. 2006, 103: 10358-10363. 10.1073/pnas.0604046103.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  44. Innocenti P, Morrow EH: Immunogenic males: a genome-wide analysis of reproduction and the cost of mating in Drosophila melanogaster females. J Evol Biol. 2009, 22: 964-973. 10.1111/j.1420-9101.2009.01708.x.

    Article  PubMed  CAS  Google Scholar 

  45. Ravi Ram K, Wolfner MF: Seminal influences: Drosophila Acps and the molecular interplay between males and females during reproduction. Integr Comp Biol. 2007, 47: 427-445. 10.1093/icb/icm046.

    Article  PubMed  CAS  Google Scholar 

  46. Avila FW, Sirot LK, Laflamme BA, Rubinstein CD, Wolfner MF: Insect seminal fluid proteins: identification and function. Ann Rev Entomol. 2011, 56: 21-40. 10.1146/annurev-ento-120709-144823.

    Article  CAS  Google Scholar 

  47. Domanitskaya EV, Liu H, Chen S, Kubli E: The hydroxyproline motif of male sex peptide elicits the innate immune response in Drosophila females. FEBS J. 2007, 274: 5659-5668. 10.1111/j.1742-4658.2007.06088.x.

    Article  PubMed  CAS  Google Scholar 

  48. Yapici N, Kim Y-J, Ribeiro C, Dickson BJ: A receptor that mediates the post-mating switch in Drosophila reproductive behaviour. Nature. 2008, 451: 33-38. 10.1038/nature06483.

    Article  PubMed  Google Scholar 

  49. Rice WR, Linder JE, Friberg U, Lew TA, Morrow EH, Stewart AD: Inter-locus antagonistic coevolution as an engine of speciation: assessment with hemiclonal analysis. Proc Natl Acad Sci U S A. 2005, 102: 6527-6534. 10.1073/pnas.0501889102.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  50. Kuijper B, Morrow EH: Direct observation of female mating frequency using time-lapse photography. Fly (Austin). 2009, 3: 118-120. 10.4161/fly.8053.

    Article  Google Scholar 

  51. Rasband WS: ImageJ. 1997, U. S. National Institutes of Health, Bethesda, Maryland, U.S.A.

    Google Scholar 

  52. Partridge L, Hoffman AA, Jones JS: Male size and mating success in Drosophila melanogaster and D. pseudoobscura under field conditions. Anim Behav. 1987, 35: 468-476. 10.1016/S0003-3472(87)80272-5.

    Article  Google Scholar 

  53. R: A Language and Environment for Statistical Computing. 2009, R Foundation for Statistical Computing, Vienna, Austria

  54. Gentleman R, Carey V, Bates D, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini A, Sawitzki G, Smith C, Smyth G, Tierney L, Yang J, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31: e15-10.1093/nar/gng015.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Gautier L, Cope L, Bolstad BM: Irizarry RA: affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004, 20: 307-315. 10.1093/bioinformatics/btg405.

    Article  PubMed  CAS  Google Scholar 

  57. Smyth G: Limma: linear models for microarray data. Bioinforma Comput Biol Solut Using R Bioconductor. 2005, 397-420. 10.1007/0-387-29362-0_23.

  58. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995, 57: 289-300.

    Google Scholar 

  59. Michaud J, Simpson K, Escher R, Buchet-Poyau K, Beissbarth T, Carmichael C, Ritchie M, Schutz F, Cannon P, Liu M, Shen X, Ito Y, Raskind W, Horwitz M, Osato M, Turner D, Speed T, Kavallaris M, Smyth G, Scott H: Integrative analysis of RUNX1 downstream pathways and target genes. BMC Genomics. 2008, 9: 363-10.1186/1471-2164-9-363.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Innocenti P, Morrow EH: The sexually antagonistic genes of Drosophila melanogaster. PLoS Biol. 2010, 8: e1000335-10.1371/journal.pbio.1000335.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Van der Laan MJ, Pollard KS: A new algorithm for hybrid hierarchical clustering with visualization and the bootstrap. J Stat Plan Inference. 2003, 117: 275-303. 10.1016/S0378-3758(02)00388-9.

    Article  Google Scholar 

  62. Chintapalli VR, Wang J, Dow JAT: Using FlyAtlas to identify better Drosophila melanogaster models of human disease. Nat Genet. 2007, 39: 715-720. 10.1038/ng2049.

    Article  PubMed  CAS  Google Scholar 

  63. Falcon S, Gentleman R: Using GOstats to test gene lists for GO term association. Bioinforma Oxf Engl. 2007, 23: 257-258. 10.1093/bioinformatics/btl567.

    Article  CAS  Google Scholar 

  64. Morrow EH, Leijon A, Meerupati A: Hemiclonal analysis reveals significant genetic, environmental and genotype x environment effects on sperm size in Drosophila melanogaster. J Evol Biol. 2008, 21: 1692-1702. 10.1111/j.1420-9101.2008.01585.x.

    Article  PubMed  CAS  Google Scholar 

  65. Morrow EH, Stewart AD, Rice WR: Assessing the extent of genome-wide intralocus sexual conflict via experimentally enforced gender-limited selection. J Evol Biol. 2008, 21: 1046-1054. 10.1111/j.1420-9101.2008.01542.x.

    Article  PubMed  CAS  Google Scholar 

  66. Radwan J, Drewniak M: Inbreeding depression for size but not for symmetry in Drosophila melanogaster. Hereditas. 2001, 134: 85-89.

    PubMed  CAS  Google Scholar 

  67. Hollis B, Houle D, Yan Z, Kawecki TJ, Keller L: Evolution under monogamy feminizes gene expression in Drosophila melanogaster. Nat Commun. 2014, 5: 3482-10.1038/ncomms4482.

    Article  PubMed  Google Scholar 

  68. Innocenti P, Flis I, Morrow EH: Data from: female responses to experimental removal of sexual selection components in Drosophila melanogaster. Dryad Digital Repository 2014. doi:10.5061/dryad.2dp4q.,

  69. Innocenti P, Morrow EH Data from: Female responses to experimental removal of sexual selection components in Drosophila melanogaster. Gene Expression Omnibus 2014. , [http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE48385]

Download references

Acknowledgments

We thank Jessica Abbott, Tim Connallon, Björn Rogell and two anonymous referees for comments on the manuscript and the Uppsala Array Platform. Funding was provided by the Swedish Research Council (Grant #2008-5533 and #2011-3701), the European Research Council (Grant #280632) and a Royal Society University Research Fellowship to EM.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Edward H Morrow.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

12862_2014_239_MOESM1_ESM.pdf

Additional file 1: Additional analyses outlining module characteristics, mating frequency data in base population and GO analyses for each module. (PDF 2 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Innocenti, P., Flis, I. & Morrow, E.H. Female responses to experimental removal of sexual selection components in Drosophila melanogaster. BMC Evol Biol 14, 239 (2014). https://doi.org/10.1186/s12862-014-0239-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-014-0239-3

Keywords