Transcriptional and epigenetic responses to mating and aging in Drosophila melanogaster
BMC Genomics volume 15, Article number: 927 (2014)
Phenotypic plasticity allows organisms to respond rapidly to changing environmental circumstances, and understanding its genomic basis can yield insights regarding the underlying genes and genetic networks affecting complex phenotypes. Female Drosophila melanogaster undergo dramatic physiological changes mediated by seminal fluid components transferred upon mating, including decreased longevity. Their physiological and behavioral effects have been well characterized, but little is known about resulting changes in regulation of gene expression or the extent to which mating-induced changes in gene expression are the same as those occurring during aging.
We assessed genome-wide mRNA, microRNA, and three common histone modifications implicated in gene activation for young and aged virgin and mated female D. melanogaster in a factorial design. We identified phenotypically plastic transcripts and epigenetic modifications associated with mating and aging. We used these data to derive phenotypically plastic regulatory networks associated with mating of young flies, and aging of virgin and mated flies. Many of the mRNAs, microRNAs and epigenetic modifications associated with mating of young flies also occur with age in virgin flies, which may reflect mating-induced accelerated aging. We functionally tested the plastic regulatory networks by overexpressing environmentally sensitive microRNAs. Overexpression resulted in altered expression of ~70% of candidate target genes, and in all cases affected oviposition.
Our results implicate microRNAs as mediators of phenotypic plasticity associated with mating and provide a comprehensive documentation of the genomic and epigenomic changes that accompany mating- and aging-induced physiological changes in female D. melanogaster.
Phenotypic plasticity is the phenomenon whereby organisms with the same genotype can achieve different phenotypic endpoints in response to changes in external or internal environments. Phenotypic plasticity thus allows organisms to respond rapidly to changing environmental circumstances. Understanding the genomic basis of phenotypic plasticity can yield insights regarding the underlying genes and genetic networks affecting complex phenotypes, since mutations in genes whose expression changes in different environments also affect the phenotype [1–5].
Most previous studies of genome wide transcriptional responses to environmental perturbations have utilized microarray technology, which has a lower dynamic range than direct sequencing of mRNAs and can be influenced by DNA polymorphisms between the target genotype and the reference genotype used to design the array. We used whole genome RNA sequencing to alleviate these limitations and included a comprehensive analysis of microRNAs (miRNAs). MicroRNAs bind to the 3′ untranslated regions of mRNA, preventing its translation, polyadenylation or stability, and play essential roles in signaling pathways that regulate development and differentiation [6, 7]. Although miRNAs have been implicated in lifespan determination in Caenorhabditis elegans and Drosophila[8, 9], little is known regarding their role in phenotypic plasticity. Finally, epigenetic modifications have been proposed as an important aspect of phenotypic plasticity through their regulation of gene expression [10–12]. However, we know little of the complex relationships between genome-wide transcriptional responses to environmental perturbations and their regulation by miRNAs and epigenetic modifications.
Here, we utilized next generation sequencing technology to assess the plastic responses of mRNAs, miRNAs and three major histone modifications (H3K4me1, H3K4me3, and H3K9ac) to mating and aging in D. melanogaster females, in a factorial design. Mating and aging have previously been shown to elicit transcriptional responses in Drosophila[4, 9, 13–20]. During mating, male D. melanogaster transfer a complex mixture of accessory gland derived peptides and proteins together with sperm into the female reproductive tract. These peptides exert profound physiological and behavioral changes in the female which are thought to enhance the male’s competitive reproductive success [21–23]. They activate oocyte development and accelerate egg laying, reduce receptivity of the female to additional matings, and enhance the female’s immune competence [21–23]. In addition, mated females have a reduced lifespan compared to virgin females [4, 24]. The effects of several accessory gland peptides have been well characterized, including ovulin  and the sex peptide . The latter exerts its effect by binding to receptors that are expressed in the reproductive tract and in the central nervous system . While the abundances of many hundreds of transcripts are known to change with age [2, 8, 9, 14–16, 20, 27], it is not known to what extent these patterns vary between virgin and mated flies and whether the reduced longevity of mated females reflects an acceleration of normal senescence. We used our data to construct interaction networks between phenotypically plastic miRNAs, their target genes and associated histone 3 markers. We functionally validated the effects of miRNAs on oviposition in mated females to demonstrate that miRNAs are instrumental in facilitating post-mating physiological and behavioral changes in female D. melanogaster.
We assessed changes in genome wide mRNA, miRNA and epigenetic modifications for 3–5 day old and 4-week old virgin Canton S (B) D. melanogaster females and females which were maintained in the continuous presence of males with multiple opportunities for mating. We refer to the latter as “mated” females. We evaluated overall differences in transcriptional and epigenetic responses to mating and to aging, as well as differences in the molecular signatures of plasticity between young and old virgin and mated flies, and mated and virgin young and old flies. This design enables us to assess whether or to what extent mating induces changes in molecular response profiles that would otherwise occur later in life.
We obtained mRNA expression profiles for each of the four physiological conditions by RNA sequencing. An average 95.7% of sequences had quality scores higher than 37.65, with over 65× coverage of the estimated 30 Mb Drosophila transcriptome. We detected 26,151 mRNA sequences, including 15,610 annotated transcripts and 10,541 previously unannotated transcripts. Some of the unannotated transcripts could be artifacts of mapping short reads corresponding to premature unspliced mRNAs back to the genome; others represent novel intergenic long non-coding RNAs. We identified 647 unique genes that were differentially expressed post-mating or as a function of age. These include 344 genes with predicted transcripts of unknown function, 23 annotated non-coding transcripts and two novel (i.e., previously unannotated) genes. The remaining 278 genes included multiple members of gene families, such as the IM, Jonah, Cp, Cpr, Twdl, Osi, Vm, Mal and Obp families (Additional file 1). Among the differentially expressed genes, 474 showed altered transcript levels in young flies after mating; 160 showed altered transcript levels between young and old virgin flies; and 279 genes showed altered expression between young and old mated flies (Additional file 1). There was modest overlap among these categories (Additional file 2: Figure S1). Only 18 differentially expressed transcripts were identified between aged virgin and mated flies. This suggests that the largest post-mating differences in transcript profiles are observed in young females and these differences fade as flies age.
We performed Gene Ontology (GO) enrichment analyses for the two treatments with the largest numbers of environmentally responsive transcripts: the plastic response to mating for young flies (Figure 1a) and the plastic response to aging of virgin flies (Figure 1b). Transcripts with altered expression levels following mating show strong enrichment of GO terms associated with egg development and proteolysis, in agreement with a previous study ; as well as GO terms associated with bacterial defense and immune responses (Figure 1a). GO categories associated with immune/defense response, egg shell formation and metabolism were enriched in the comparison between young and old virgin flies (Figure 1b), also in agreement with previous studies [14–16]. There is partial concordance between GO categories between the transcript profiles of young flies before and after mating and young and aged virgin flies, indicating that some biological processes undergo specific post-mating or aging-related changes in transcript abundance levels, whereas others reflect post-mating changes that would occur later in the lifespan in the absence of mating.
We extracted and sequenced small RNAs (<200 nt, including miRNAs) from the same samples used to assess genomic responses of mRNAs to mating and aging. Examination of the distribution profiles of small RNAs showed prominent peaks at 22 bp and 30 bp, corresponding to miRNAs and rasiRNAs (the Drosophila equivalent of piRNAs), respectively (Figure 2). A prominent peak of small RNA fragments around 15 bp is evident in replicate samples of flies following mating. These fragments are transient and are not observed for older mated flies, nor are they observed in young or old virgin flies. These small mRNA fragments are likely due to degradation of maternal RNA in embryos during the maternal to zygotic transition following fertilization [29, 30]. We mapped the 15-18 bp RNA fragments across all samples to the genome sequence and identified 147 genes that give rise to degradation fragments exclusively in young flies after mating (Additional file 3: Table S1). Among these genes, 64 encode predicted transcripts of unknown function. These genes also comprise genes associated with developmental processes and, intriguingly, four genes encoding chemoreceptors (Gr39b, Gr93d, Or13a, Or45a).
We detected 149 mature miRNAs from 105 precursors, of which 41 miRNAs from 39 precursors with 28 conserved seeds showed significant differences across the four conditions (Additional file 4). Similar numbers of miRNAs show plastic responses to the four treatments. We find 27 differentially regulated miRNAs after mating in young flies and 22 in old flies; and 17 miRNAs differentially regulated with aging in virgin flies and 28 in mated flies, with considerable overlap in plastic responses of miRNAs among the four physiological conditions (Additional file 5: Figure S2). Of particular interest is the miRNA 309 cluster, which, in addition to mir-309, comprises mir-286, mir-3, mir-4, mir-5 and mir-6. The latter gives rise to three alternative stem-loop configurations with similar conserved seed regions [31, 32]. Zygotic expression of these miRNAs has been associated with degradation of maternal mRNA in embryos [29, 32]. Indeed, these miRNAs are only expressed in mated flies and are virtually absent in virgins (Figure 3). Among the 147 genes that are uniquely implicated in maternal degradation, 45 contain 21 conserved binding sites for miRNAs that are differentially expressed after mating, and 18 genes are targets of the miRNA 309 cluster.
Epigenetic regulation in D. melanogaster occurs primarily through association of modified histones with target DNA sequences . Numerous post-translational modifications of histone proteins have been identified, including methylation, acetylation, phosphorylation, ubiquitylation, SUMOylation, biotinylation and ADP-ribosylation [11, 34, 35]. Some histone modifications act in cis to change chromatin structure directly, whereas others act in trans to affect the recruitment of a protein complex that regulates gene expression . Modified histone protein residues can serve as docking sites for transcription factors and other chromatin-modifying enzymes to regulate transcription or induce further chromatin remodeling [11, 12, 34, 35]. We examined three common histone 3 modifications (H3K4me1, H3K4me3, H3K9ac) previously implicated in gene activation .
We used immuno-coprecipitation to identify DNA segments differentially associated with each of the histone 3 modifications across the genome. We identified 3,484, 5,467, and 4,986 peaks of co-precipitated DNA fragments for H3K4me1, H3K4me3 and H3K9ac, respectively; of which 710 (H3K4me1), 1,513 (H3K4me3) and 1,819 (H3K9ac) co-precipitated fragments were differentially expressed among the four physiological conditions (Additional file 6: Figure S3). We found similar numbers of up- and down-regulated genes for each variable histone modification (Table 1; Additional files 7, 8 and 9). The number of histone modifications associated with aging was larger than those associated with post-mating changes, especially for H3K4me3 (Table 1). We identified candidate histone regulated genes associated with co-precipitated fragments by applying a 2 kb window downstream and upstream of each gene. On average, 3.3 candidate genes were associated with each histone peak (Table 1), of which ~24% are novel transcripts.
About 10% of differentially modified fragments after mating in young flies were shared between H3K4me1 and H3K4me3, and between H3K4me1 and H3K9ac, while about 15% were shared among H3K4me3 and H3K9ac, and only about 1% of these fragments were associated with all three modified histone marks (Additional file 10: Figure S4). Similarly, about 17% of differentially modified fragments identified as a function of aging of virgin females were shared between H3K4me1 and H3K4me3, and between H3K4me1 and H3K9ac, while as much as 36% of precipitated DNA sequences were shared among H3K4me3 and H3K9ac, and about 5% of sequences were associated with all three modified histone marks (Additional file 10: Figure S4).
We asked to what extent histone modifications induced by mating are permanent throughout the lifespan of reproductively active females. We found that 18%, 9% and 11% of H3K4me1, H3K4me3 and H3K9ac modifications, respectively, that occurred post-mating in young flies persisted at 4 weeks of age when females were maintained in the presence of males. Next, we asked how many of these histone modifications also occur as a result of aging in the absence of mating. We found that 65%, 73% and 77% of these H3K4me1, H3K4me3 and H3K9ac modifications, respectively, ultimately occur in virgin flies, but only at later age. These markers, therefore, represent epigenetic mating-induced accelerated aging changes and are associated with up to 539 unique candidate genes, of which 421 are annotated (Additional file 11). These include genes associated with oocyte differentiation, including stonewall, 14-3-3zeta, bazooka, missing oocyte, and egalitarian. Histone modifications that affect these candidate genes may reflect persistent post-mating stimulation of oogenesis.
To assess to what extent environmentally plastic transcripts might be causally associated with epigenetic regulation we asked which putative histone regulated genes showed altered transcriptional regulation across the four experimental conditions (Table 2). The majority of genes with altered transcriptional regulation corresponded with candidate target genes for the three histone marks when young virgin flies are compared to young mated flies and when young mated flies are compared to aged mated flies (Table 2). Many fewer genes with altered transcriptional regulation corresponded with candidate target genes when young virgin flies were compared to aged virgins and when old virgin flies were compared to old mated flies (Table 2). Although these proportions roughly parallel differences in the numbers of transcripts with altered abundance levels as a consequence of mating or aging, they illustrate that different transcripts interact with histone 3 upon mating and as a function of age. Among genes with altered transcriptional regulation that interact with modified histone 3, only a few are associated with more than one histone mark, indicating target specific interactions for each histone modification with these genes (Table 2).
Phenotypically plastic regulatory networks
We identified target genes that undergo changes in expression for all miRNAs with altered transcriptional regulation, and recruited genes known to interact with these target genes from the Drosophila Interactions Database  to obtain integrative networks that provide a comprehensive visualization of post-mating or aging-associated changes in whole genome transcriptional profiles, and superimposed histone marks on the networks. We only included genes with more than one interaction with miRNA target genes. Since very few transcripts were altered in the comparison of old virgin and mated flies, we derived these plastic regulatory networks for three conditions: young virgin and mated flies, young and old mated flies, and young and old virgin flies.
The network associated with post-mating changes in young flies comprised 46 target genes for 20 miRNAs, and these target genes have 92 interacting partners (Figure 4; Additional file 12). In addition to genes that encode predicted transcripts of unknown function, developmental genes feature prominently among the latter. The network reveals several notable features. First, there are shared targets for multiple miRNAs. Second, these target genes interact with multiple partners and each interacting partner interacts with multiple miRNA targets, thus forming a collage of interacting ensembles that are embedded in the network. Third, 13% of miRNA target genes and 37% of their interacting partners are targets for regulation through histone modifications. In addition, two miRNAs, mir-193 and mir-1000, are themselves potential targets for H3K4me3. It is reasonable to infer from the network structure that regulatory modifications at a limited number of focal genes can affect the expression of a suite of interacting partners.
A similar analysis of transcriptional changes between young and old virgin flies gives rise to a smaller network, comprised of only seven miRNAs, five target genes, and 14 genes interacting with these targets (Figure 5; Additional file 12). Of the latter, 10 are potential targets for regulation by modified histone 3, as are two of the miRNA target genes. Among the seven miRNAs, five overlap with the network associated with post-mating changes, including mir-193. It should be noted that in this case the composition of differentially expressed miRNA target genes restricts connectivity, which limits the size of the network. Nevertheless, it is clear that the genetic architectures that underlie changes in gene expression are distinct with regard to post-mating changes and aging, as might be expected.
The network associated with aging in mated flies comprises 15 miRNAs, 27 miRNA target genes and 64 interacting partners (Figure 6; Additional file 12). Two miRNAs, mir-193 and mir-279 are potential targets of H3K4me3 and H3K4me1, respectively. Eight miRNA target genes (~29%) and 39 (~60%) interacting partners are potentially regulated by histone modifications.
Effects of micro-RNA overexpression on target gene expression
To test the connectivity of the computational networks (Figures 4, 5 and 6) and evaluate causality, we overexpressed six miRNAs for which overexpression lines were available (mir-281-1, mir-286, mir-34, mir-92b, mir-310 and mir-988) using the binary GAL4-UAS expression system  under a universal ubiquitin promoter. Evidence that the driver was effective comes from the observation that no viable offspring could be obtained when mir-281-1 and mir-310 were overexpressed. The remaining overexpression lines generated F1 progeny, in which we used quantitative RT-PCR to monitor expression of target genes in virgin and mated flies contemporaneously with F1 progeny obtained by crossing the control progenitor strain to the ubiquitin-GAL4 driver line. We assessed whether changes in their corresponding target genes that were differentially expressed between young virgin and mated females were significantly different after mating in control lines compared to the miRNA overexpression lines (Student’s t-test, P <0.05). Despite differences in genetic backgrounds between the GAL4-UAS offspring and the original Canton S (B) line used to derive the networks, we observed significant effects for 13 of the 20 target genes tested (70%; Table 3). Thus, perturbation of miRNAs results in disruption of target gene expression that occurs during changes in physiological state and directly implicates miRNAs in the regulation of post-mating and aging-dependent gene expression.
Effects of micro-RNA overexpression on egg laying
Stimulation of oogenesis and oviposition is the most profound change induced in females as a consequence of mating . To assess causality between post-mating modulation of miRNA expression and oviposition, we measured egg laying after mating by control females and females overexpressing mir-286, mir-34, mir-92b, or mir-988. Mating causes a reduction in expression of mir-34, mir-92b and mir-988 and an increase in expression of mir-286 (Figure 3; Additional file 4). Females that overexpress mir-34 and mir-92b lay 34% and 37% fewer eggs after mating, whereas females that overexpress mir-286 and mir-988 lay 16% and 33% more eggs than controls (Figure 7). Thus, post-mating changes in miRNA expression levels are causally associated with physiological changes.
We identified genomic and epigenomic changes in gene expression that are associated with mating and aging in female D. melanogaster. We found that 26 miRNAs change expression levels after mating, including transient expression of members of the mir-309 cluster associated with the maternal to zygotic transition that results in maternal mRNA breakdown [29, 30]. The mir-309 cluster is analogous to mir-430, which mediates the maternal to zygotic transition in zebrafish . Up-regulation of the mir-309 cluster after mating is accompanied by the transient accumulation of small RNA fragments, likely representing breakdown products of degraded maternal RNAs in embryos. In D. melanogaster, loci required for the destabilization of distinct subsets of maternal mRNAs during the maternal to zygotic transition have been mapped to different chromosome arms , indicating that different mechanisms may mediate maternal RNA degeneration in the zygote. The appearance of RNA fragments only in young flies after mating could be due to subsequent changes that may include epigenetic modifications, which increase the efficiency of the removal of such fragments during repeated matings.
We identified target genes for miRNAs that change transcript abundance levels after mating. We examined similar changes in genome-wide expression levels of miRNAs and their associated target genes that occur as females age. In addition, we analyzed epigenetic modifications mediated via histone 3 that are associated with mating and aging. We were able to construct integrative genetic networks from a subset of genes with altered expression due to mating or aging to visualize part of the genetic architectures that underlie these physiological changes. We found that the largest differences in transcript profiles occur in young females after mating. Although the genetic architectures that underlie changes in gene expression are distinct with regard to post-mating changes and aging, we found epigenetic modifications that occur post-mating in young flies and at later age in virgin flies, consistent with mating-induced accelerated aging and consequent shortening of the lifespan of mated females  (Additional file 3).
The power of our analysis is limited by the small number of replicate samples and histone modifications analyzed, and possible differences in genomic responses in different genetic backgrounds. Thus, the complement of genes that undergo altered regulation as a consequence of physiological state and epigenetic modifications that could regulate their expression provides only a partial view of the genetic architectures that regulate physiological state-induced regulation of gene expression. Furthermore, genes associated with histone marks are commonly assigned based on chromosomal location [36, 41, 42]. Demonstrating causal relationships with candidate genes within such chromosomal regions remains challenging, and it is difficult to unambiguously identify false positives or exclude long-distance effects of modified histones on gene expression. We alleviated this concern to some extent by focusing on histone modifications associated with co-localized genes that undergo altered expression. We also note that our analyses were performed on whole flies and we cannot exclude that some changes in miRNA expression or histone modifications may be compartmentalized to certain tissues.
A large proportion of genes with changes in transcript abundance are members of gene families. Members of the Jonah, Cp, Twdl, and Mal gene families show antagonistic changes in transcript abundances after mating and as a result of aging, whereas Vm family members show changes in the same direction. Similarly, some genes with altered expression that belong to gene ontology categories shared between mating and aging conditions also show antagonistic effects. We speculate that these antagonistic changes may be part of the molecular machinery that mediates trade-offs between reproduction and lifespan . In addition to restricted expression in a limited number of cells, antagonistic regulation in different tissues could also be a limiting factor in detecting genes that undergo altered transcriptional regulation (e.g. we did not detect altered transcript abundance of smaug, which triggers the maternal to zygotic degradation [30, 43]).We observed many histone modifications after mating or upon aging that are not accompanied by changes in corresponding transcript abundances. Here, the target gene might be remote from the histone-DNA interaction site. Alternatively, interaction between the histone and the target gene might prevent activation of this gene at a later stage or inactivation of a gene that is actively expressed. Furthermore, the networks centered on miRNAs and their target genes (Figures 4, 5 and 6) do not include all transcripts with altered expression levels, nor do they include all potential target genes for epigenetic regulation. Nevertheless, the complex structure of these networks raises the possibility that epigenetic regulation of target genes for histone 3 could also indirectly affect the expression of interconnected genes.
In addition to epigenetic modifications that uniquely occur after mating, we identified identical epigenetic markers that occur also in virgin females at advanced age. These mark post-mating changes that relate to the shortening of female lifespan after mating, possibly due to the increased energy expenditure required for egg production [21, 24].
Overexpression of miRNAs confirmed the biological relevance of the derived environmentally plastic regulatory networks. Despite differences in genetic backgrounds, altered expression of miRNAs correlated with altered expression of ~70% of their candidate target genes, and in all cases miRNA overexpression affected post-mating oviposition. Although we analyzed only a subset of miRNAs for which viable overexpression lines could be obtained, it is reasonable to assume that these observations would extend to the larger ensemble of miRNAs.
Since this study focused on whole flies, it remains to be determined to what extent the genomic and epigenomic changes we have observed are compartmentalized to reproductive organs or the central nervous system. Previous studies have implicated steroid hormones in regulating expression of miRNAs [44, 45]. Odorant binding proteins expressed in the male accessory gland could conceivably serve as carriers for such ligands . However, the precise mechanisms by which signals transmitted during mating regulate gene expression remain to be explored. Nevertheless, this study presents to the best of our knowledge the first demonstration that miRNAs play a critical role in mediating post-mating changes in female flies and a comprehensive documentation of the genomic and epigenomic changes that accompany mating-induced physiological changes as well as aging in female D. melanogaster.
To gain insight regarding the genes and regulatory networks underpinning plastic changes after mating and during aging, we identified mRNAs, microRNAs and epigenetic modifications with significantly different expression after mating or as a function of aging. We used these data to derive phenotypically plastic regulatory networks centered on environmentally sensitive microRNAs associated with aging and mating, and identified several biomarkers of mating-induced accelerated aging. Overexpression of several plastic microRNAs resulted in altered expression of candidate target genes and affected oviposition. MicroRNAs are thus critical in mediating post-mating changes in gene expression. These data provide comprehensive documentation of the genomic and epigenomic changes that accompany mating- and aging-induced physiological changes in female D. melanogaster.
Isogenic Canton S (B) flies  were reared on cornmeal-molasses-yeast medium at 25°C under a 12 hour light–dark cycle. Adults were collected immediately after eclosion and subjected to two treatments. In the first treatment, adults were allowed to mate at a density of 25 females and 25 males per vial. In the second treatment, 50 females were placed in each vial and maintained as virgins. Flies were transferred to new vials every 2 days. Female flies were collected at 3–5 days (young flies) or at 4 weeks (aged flies) at the same time of day (1 pm-3 pm) and flash frozen on dry ice.
mRNA and miRNA sequencing
Two independent samples of 25 female flies were used for each condition (young virgin flies, young mated flies, aged virgin flies and aged mated flies). RNA was extracted and separated into a small RNA (<200 nts) and large RNA fraction (>200 nts) using the Qiagen miRNeasy kit. rRNA was depleted from the large RNA fraction using Ribo-Zero rRNA Removal Kits (Epicentre, Inc.) to enrich for mRNA. We used enriched mRNA to prepare bar-coded libraries. Four mRNA samples were pooled in equal molarity from each replicate and sequenced in one 50 bp single read Illumina HiSeq2000 lane. Small RNA libraries were prepared using the NEXflex Small RNA Sequencing kit (Bioo, Inc.) and enriched for miRNA, siRNA and piRNA following gel electrophoresis. All 8 bar-coded samples were pooled in equal molarities and sequenced in one 50 bp single read Illumina HiSeq2000 lane.
mRNA sequences were aligned and assembled using Bowtie2-2.0.6, Tophat-2.0.7 and Cufflinks-2.0.2 [46, 47]. We used Cuffdiff to analyze pairwise differential expression between the conditions of transcripts, genes, splicing and promoter uses . We first trimmed small RNA sequences to remove the adapter sequence before aligning them to the Drosophila miRBase  using miRExpress . We then used quantile normalized counts for differential analysis implemented by the R package (lmPerm) using permutation with a one-way ANOVA model, Y = μ + C + ϵ, where μ is the overall mean, C designates condition and ϵ the error term. We used FDR to correct for multiple-testing under dependency . One pooled standard deviation was used for post hoc pairwise comparisons. We identified putative miRNA targets using the TargetScan Fly  and MinoTar data bases . Raw data can be accessed in the National Center for Biotechnology Information Sequence Read Archive (SRA) under accession number SRP048388.
Chromatin immunoprecipitation and ChIP sequencing
Two independent samples of 400 female flies were used for each condition. Flies were ground in liquid nitrogen. Samples were then cross-linked with 1% formaldehyde for 10 min and lysed, and chromatin was digested into mono- and di-nucleosomes with micrococcal nuclease. Aliquots of digested chromatin samples were then immunoprecipitated with antibodies against H3K4me1, H3K4me3, and H3K9ac. Antibodies of which the specificities have been ascertained by the modENCODE project [36, 54] were obtained from Abcam, Inc., and Active Motif, Inc.
Enriched chromatin extracts were used to prepare bar-coded sequencing libraries. ChIP libraries of the four conditions from each replicate were pooled and sequenced in one 50 bp single read Illumina HiSeq2000 lane. ChIP sequences were aligned to the reference genome using Bowtie2-2.0.6 followed by peak calling using MACS1.4.2 (Model-based Analysis for ChIP-Seq) [47, 55]. Peak calls were merged across samples if the peak was present in more than one sample and the peak width bigger than 146 bp. Sequence counts were then summed for each peak, and normalized across samples by calculating Fragment Per Kilobase per Million reads mapped (FPKM), corrected for background using median FPKM of non-peak intervals in each sample and followed by quantile normalization across samples. We used background corrected normalized values of FPKM for each peak for differential modification analysis, as described above, using permutation with one-way ANOVA and FDR to correct for multiple-testing under dependency . We used one pooled standard deviation for post hoc pairwise comparisons. Raw data can be accessed in the National Center for Biotechnology Information Sequence Read Archive (SRA) under accession number SRP048404.
Gene ontology analysis
Gene ontology enrichment analyses for transcripts that changed expression after mating or aging was done using DAVID  with a cut-off for nominal statistical significance of P < 0.01.
Confirmation of interactions between miRNAs and target genes
We obtained UAS-mir-281-1, UAS-mir-286, UAS-mir-34, UAS-mir-92b, UAS-mir-988, mir-310 and the co-isogenic control lines from the Bloomington Stock Center. Females from each homozygous UAS-mir line and their progenitor control were crossed with males from a ubiquitin-GAL4 driver line (w1118: Ubi-Gal4). We did not obtain viable offspring when mir-281-1 or mir-310 was overexpressed under the ubiquitin promoter. For the other crosses 4–6 day old virgin and mated F1 females were flash frozen for RNA extraction. We extracted total RNA from 10 female flies of each replicate with three biological replicates of each cross. RNAs were converted to cDNA using the iScript cDNA Synthesis Kit (Bio-Rad, Inc.). We used real time quantitative PCR on an ABI 7900HT qRT-PCR thermal cycler with 3 technical replicates for each sample and ABsolute qPCR SYBR Green Mixes (Thermo Scientific, Inc.) to evaluate the difference in expression of miRNA target genes between virgin and mated F1 females in controls and between virgin and mated F1 females that overexpress miRNAs (Table 3) using the following primers: CCCTGCTCCTTCTCCTTCTT and TGAACTTGAACCCCACCTTG for Mlc1, GGTCATCGTCTGTCCGATCT and CTAACTTCCTGGGCGACAAC for dyl, CACCCCAATTCACTCTGGTT and CCTTCAGCCGAAGATTGAAA for Egfr, GGAGACACAGCTCCGACAAT and ATTCTGGCTACCCGGCTTAT for Vha68-1, TTTCGGTGTCTTGTGTCTGG and AGGCAGTCGCTGTAGATGCT for Osi7, GTCCCTGATCGGCGATAATA and AACCCAGGAGCTTCTTCAACT for CG15212, TCCGAGTACTCTCCCTCCAG and TCTTGCTAGTCCAGCCACCT for CG4000, CAGTAGTGGAGCAGGCAGTG and ATGATAGGCCTCCTCCTTGG for Pkcdelta, GGGATAGGCTGGATTGGATT and CGTCGAGAGTAGCTCCGATT for CG31869, AGTAGTTGCGATGCCAATCC and AGGCGCTATAATTTCGATGC for CG8303, TGATCTTGGCAGTGGACTTG and AAGATTCTGCGCAAGAAGGA for Prm, CACTGGCATACCTTTGGTTG and TGGCAGTGACATTGATTCGT for fln, AGAATTGCCTTTTGCCACA and CCGACTACGATTCGAAAAGC for Tusp, ATTGCAGTGTGGTTGTTCCA and CCAGTTTGGTTGGAGACGAT for CG34040, GCGAAAGAGAGTCCGAGAGA and GGTTGTGTGCTGGTCAGTTG for dan, CATCGGCATCAGCATCTATC and GCATAGATGGATGGGTGGAC for CG16867, GTAATGACGGTGGCCTTTTC and CTTTGACCATCTTGGCCACT for Gbeta76C, AGGTCTTGGGCCTAGGTGTT and GAGGCCGATCCTGATGAATA for prc, CAAATCCTGCGCCTCATAGT and CACATCGAGTGCCTTGGATA for CG15571, GGTGGAGCTGGAGAACTACG and AGGTCTCTCCGTCGACATTC for CG31086, TCACCAGATTCCCGCTTATC and GTAACCACCGGATGATGAGG for Ela and ACTGAGAGCCAGGACTGGAA and CTTCCGGTTTCACAGATGGT for vkg. Gapdh expression (Primers: AGGCGTTTGTGACTTCTGGAA and TCTGGCCGTTGAGCATTTC) in each sample was used to normalize gene expression by mean normalization.
We collected 5–7 day old-mated F1 control females and females overexpressing miRNAs under the ubiquitin-GAL4 promoter and counted the number of eggs laid by 5 females during 18 hours (4:00 pm-11:00 am). Ten replicates of 5 F1 females were tested from each cross. We used Dunnett’s test to assess variation of egg laying between the control and miRNA overexpression crosses.
Lakowski B, Hekimi S: The genetics of caloric restriction in Caenorhabditis elegans. Proc Natl Acad Sci U S A. 1998, 95 (22): 13091-13096. 10.1073/pnas.95.22.13091.
Magwire MM, Yamamoto A, Carbone MA, Roshina NV, Symonenko AV, Pasyukova EG, Morozova TV, Mackay TFC: Quantitative and molecular genetic analyses of mutations increasing Drosophila life span. PLoS Genet. 2010, 6 (7): e1001037-10.1371/journal.pgen.1001037.
Weindruch R, Kayo T, Lee CK, Prolla TA: Microarray profiling of gene expression in aging and its alteration by caloric restriction in mice. J Nutr. 2001, 131 (3): 918S-923S.
Zhou S, Campbell TG, Stone EA, Mackay TFC, Anholt RRH: Phenotypic plasticity of the Drosophila transcriptome. PLoS Genet. 2012, 8 (3): e1002593-10.1371/journal.pgen.1002593.
Sambandan D, Carbone MA, Anholt RRH, Mackay TFC: Phenotypic plasticity and genotype by environment interaction for olfactory behavior in Drosophila melanogaster. Genetics. 2008, 179 (2): 1079-1088. 10.1534/genetics.108.086769.
Bourc’his D, Voinnet O: A small-RNA perspective on gametogenesis, fertilization, and early zygotic development. Science. 2010, 330 (6004): 617-622. 10.1126/science.1194776.
Bushati N, Cohen SM: MicroRNA functions. Annu Rev Cell Dev Biol. 2007, 23: 175-205. 10.1146/annurev.cellbio.23.090506.123406.
Kato M, Slack FJ: Ageing and the small, non-coding RNA world. Ageing Res Rev. 2013, 12 (1): 429-435. 10.1016/j.arr.2012.03.012.
Liu N, Landreh M, Cao K, Abe M, Hendriks GJ, Kennerdell JR, Zhu Y, Wang LS, Bonini NM: The microRNA miR-34 modulates ageing and neurodegeneration in Drosophila. Nature. 2012, 482 (7386): 519-523. 10.1038/nature10810.
Duncan EJ, Gluckman PD, Dearden PK: Epigenetics, plasticity, and evolution: How do we link epigenetic change to phenotype?. J Exp Zool B Mol Dev Evol. 2014, 322 (4): 208-220. 10.1002/jez.b.22571.
Gardner KE, Allis CD, Strahl BD: Operating on chromatin, a colorful language where context matters. J Mol Biol. 2011, 409 (1): 36-46. 10.1016/j.jmb.2011.01.040.
Spencer VA, Davie JR: Role of covalent modifications of histones in regulating gene expression. Gene. 1999, 240 (1): 1-12. 10.1016/S0378-1119(99)00405-9.
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 (3): 1395-1408. 10.1534/genetics.108.086934.
Pletcher SD, Macdonald SJ, Marguerie R, Certa U, Stearns SC, Goldstein DB, Partridge L: Genome-wide transcript profiles in aging and calorically restricted Drosophila melanogaster. Curr Biol. 2002, 12 (9): 712-723. 10.1016/S0960-9822(02)00808-4.
McCarroll SA, Murphy CT, Zou S, Pletcher SD, Chin CS, Jan YN, Kenyon C, Bargmann CI, Li H: Comparing genomic expression patterns across species identifies shared transcriptional profile in aging. Nat Genet. 2004, 36 (2): 197-204. 10.1038/ng1291.
Lai CQ, Parnell LD, Lyman RF, Ordovas JM, Mackay TFC: Candidate genes affecting Drosophila life span identified by integrating microarray gene expression analysis and QTL mapping. Mech Ageing Dev. 2007, 128 (3): 237-249. 10.1016/j.mad.2006.12.003.
Smith G, Fang Y, Liu X, Kenny J, Cossins AR, de Oliveira CC, Etges WJ, Ritchie MG: Transcriptome-wide expression variation associated with environmental plasticity and mating success in cactophilic Drosophila mojavensis. Evolution. 2013, 67 (7): 1950-1963. 10.1111/evo.12082.
Dalton JE, Kacheria TS, Knott SR, Lebo MS, Nishitani A, Sanders LE, Stirling EJ, Winbush A, Arbeitman MN: Dynamic, mating-induced gene expression changes in female head and brain tissues of Drosophila melanogaster. BMC Genomics. 2010, 11: 541-10.1186/1471-2164-11-541.
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 (16): 1509-1514. 10.1016/j.cub.2004.08.028.
Lai CQ, Leips J, Zou W, Roberts JF, Wollenberg KR, Parnell LD, Zeng ZB, Ordovas JM, Mackay TFC: Speed-mapping quantitative trait loci using microarrays. Nat Methods. 2007, 4 (10): 839-841. 10.1038/nmeth1084.
Wolfner MF: Battle and ballet: molecular interactions between the sexes in Drosophila. J Hered. 2009, 100 (4): 399-410. 10.1093/jhered/esp013.
Chapman T: Seminal fluid-mediated fitness traits in Drosophila. Heredity (Edinb). 2001, 87 (Pt 5): 511-521.
Wolfner MF: The gifts that keep on giving: physiological functions and evolutionary dynamics of male seminal proteins in Drosophila. Heredity (Edinb). 2002, 88 (2): 85-93. 10.1038/sj.hdy.6800017.
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 (6511): 241-244. 10.1038/373241a0.
Rubinstein CD, Wolfner MF: Drosophila seminal protein ovulin mediates ovulation through female octopamine neuronal signaling. Proc Natl Acad Sci U S A. 2013, 110 (43): 17420-17425. 10.1073/pnas.1220018110.
Yapici N, Kim YJ, Ribeiro C, Dickson BJ: A receptor that mediates the post-mating switch in Drosophila reproductive behaviour. Nature. 2008, 451 (7174): 33-37. 10.1038/nature06483.
Wilson RH, Lai CQ, Lyman RF, Mackay TFC: Genomic response to selection for postponed senescence in Drosophila. Mech Ageing Dev. 2013, 134 (3–4): 79-88.
Norga KK, Gurganus MC, Dilda CL, Yamamoto A, Lyman RF, Patel PH, Rubin GM, Hoskins RA, Mackay TFC, Bellen HJ: Quantitative analysis of bristle number in Drosophila mutants identifies genes involved in neural development. Curr Biol. 2003, 13 (16): 1388-1396. 10.1016/S0960-9822(03)00546-3.
Barckmann B, Simonelig M: Control of maternal mRNA stability in germ cells and early embryos. BBA-Gene Regul Mech. 2013, 1829 (6–7): 714-724.
Tadros W, Lipshitz HD: The maternal-to-zygotic transition: a play in two acts. Development. 2009, 136 (18): 3033-3042. 10.1242/dev.033183.
Ruby JG, Stark A, Johnston WK, Kellis M, Bartel DP, Lai EC: Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs. Genome Res. 2007, 17 (12): 1850-1864. 10.1101/gr.6597907.
Bushati N, Stark A, Brennecke J, Cohen SM: Temporal reciprocity of miRNAs and their targets during the maternal-to-zygotic transition in Drosophila. Curr Biol. 2008, 18 (7): 501-506. 10.1016/j.cub.2008.02.081.
Lyko F, Ramsahoye BH, Jaenisch R: DNA methylation in Drosophila melanogaster. Nature. 2000, 408 (6812): 538-540. 10.1038/35046205.
Allis CD, Berger SL, Cote J, Dent S, Jenuwien T, Kouzarides T, Pillus L, Reinberg D, Shi Y, Shiekhattar R, Shilatifard A, Workman J, Zhang Y: New nomenclature for chromatin-modifying enzymes. Cell. 2007, 131 (4): 633-636. 10.1016/j.cell.2007.10.039.
Kouzarides T: SnapShot: Histone-modifying enzymes. Cell. 2007, 131 (4): 822-10.1016/j.cell.2007.11.005.
Nègre N, Brown CD, Ma L, Bristow CA, Miller SW, Wagner U, Kheradpour P, Eaton ML, Loriaux P, Sealfon R, Li Z, Ishii H, Spokony RF, Chen J, Hwang L, Cheng C, Auburn RP, Davis MB, Domanus M, Shah PK, Morrison CA, Zieba J, Suchy S, Senderowicz L, Victorsen A, Bild NA, Grundstad AJ, Hanley D, MacAlpine DM, Mannervik M, et al: A cis-regulatory map of the Drosophila genome. Nature. 2011, 471 (7339): 527-531. 10.1038/nature09990.
Murali T, Pacifico S, Yu J, Guest S, Roberts GG, Finley RL: DroID 2011: a comprehensive, integrated resource for protein, transcription factor, RNA and gene interactions for Drosophila. Nucleic Acids Res. 2011, 39 (Database issue): D736-D743.
Brand AH, Perrimon N: Targeted gene expression as a means of altering cell fates and generating dominant phenotypes. Development. 1993, 118 (2): 401-415.
Giraldez AJ, Mishima Y, Rihel J, Grocock RJ, Van Dongen S, Inoue K, Enright AJ, Schier AF: Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science. 2006, 312 (5770): 75-79. 10.1126/science.1122689.
De Renzis S, Elemento O, Tavazoie S, Wieschaus EF: Unmasking activation of the zygotic genome using chromosomal deletions in the Drosophila embryo. PLoS Biol. 2007, 5 (5): e117-10.1371/journal.pbio.0050117.
Yin H, Sweeney S, Raha D, Snyder M, Lin H: A high-resolution whole-genome map of key chromatin modifications in the adult Drosophila melanogaster. PLoS Genet. 2011, 7 (12): e1002380-10.1371/journal.pgen.1002380.
Nagarajan M, Veyrieras JB, de Dieuleveult M, Bottin H, Fehrmann S, Abraham AL, Croze S, Steinmetz LM, Gidrol X, Yvert G: Natural single-nucleosome epi-polymorphisms in yeast. PLoS Genet. 2010, 6 (4): e1000913-10.1371/journal.pgen.1000913.
Benoit B, He CH, Zhang F, Votruba SM, Tadros W, Westwood JT, Smibert CA, Lipshitz HD, Theurkauf WE: An essential role for the RNA-binding protein Smaug during the Drosophila maternal-to-zygotic transition. Development. 2009, 136 (6): 923-932. 10.1242/dev.031815.
Jin H, Kim VN, Hyun S: Conserved microRNA miR-8 controls body size in response to steroid signaling in Drosophila. Genes Dev. 2012, 26 (13): 1427-1432. 10.1101/gad.192872.112.
Kucherenko MM, Barth J, Fiala A, Shcherbata HR: Steroid-induced microRNA let-7 acts as a spatio-temporal code for neuronal cell fate in the developing Drosophila brain. EMBO J. 2012, 31 (24): 4511-4523. 10.1038/emboj.2012.298.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL: TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14 (4): R36-10.1186/gb-2013-14-4-r36.
Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9 (4): 357-359. 10.1038/nmeth.1923.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L: Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012, 7 (3): 562-578. 10.1038/nprot.2012.016.
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34: D140-D144. 10.1093/nar/gkj112.
Wang WC, Lin FM, Chang WC, Lin KY, Huang HD, Lin NS: miRExpress: Analyzing high-throughput sequencing data for profiling microRNA expression. BMC Bioinformatics. 2009, 10: 328-10.1186/1471-2105-10-328.
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I: Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001, 125 (1–2): 279-284.
Kheradpour P, Stark A, Roy S, Kellis M: Reliable prediction of regulator targets using 12 Drosophila genomes. Genome Res. 2007, 17 (12): 1919-1931. 10.1101/gr.7090407.
Schnall-Levin M, Zhao Y, Perrimon N, Berger B: Conserved microRNA targeting in Drosophila is as widespread in coding regions as in 3ʹUTRs. Proc Natl Acad Sci U S A. 2010, 107 (36): 15751-15756. 10.1073/pnas.1006172107.
Egelhofer TA, Minoda A, Klugman S, Lee K, Kolasinska-Zwierz P, Alekseyenko AA, Cheung MS, Day DS, Gadel S, Gorchakov AA, Gu T, Kharchenko PV, Kuan S, Latorre I, Linder-Basso D, Luu Y, Ngo Q, Perry M, Rechtsteiner A, Riddle NC, Schwartz YB, Shanower GA, Vielle A, Ahringer J, Elgin SC, Kuroda MI, Pirrotta V, Ren B, Strome S, Park PJ, et al: An assessment of histone-modification antibody quality. Nat Struct Mol Biol. 2011, 18 (1): 91-93. 10.1038/nsmb.1972.
Feng J, Liu T, Qin B, Zhang Y, Liu XS: Identifying ChIP-seq enrichment using MACS. Nat Protoc. 2012, 7 (9): 1728-1740. 10.1038/nprot.2012.101.
da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.
We thank Rachel Lyman, Sarah McAdam and Genevieve Evans for technical assistance and Dr. Wen Huang for valuable discussions.
The authors declare that they have no competing interests.
SZ, TFCM and RRHA designed the experiments, SZ performed the experiments and analyzed the data, SZ, TFCM and RRHA wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Differences in gene expression levels. Table S1. Differences in gene expression levels between young and aged female flies. Table S2. Differences in gene expression levels between virgin and mated female flies. (XLSX 119 KB)
Additional file 2: Figure S1: Four-way Venn diagram that shows overlap among transcripts that change in abundance level after mating of young flies (green circle) and aged flies (red circle), and as a result of aging in virgin flies (blue circle) and mated flies (yellow circle). (TIFF 2 MB)
Additional file 5: Figure S2: Four-way Venn diagram that shows overlap among miRNAs that change in abundance level after mating of young flies (green circle) and aged flies (red circle), and as a result of aging in virgin flies (blue circle) and mated flies (yellow circle). (TIFF 2 MB)
Additional file 6: Figure S3: Four-way Venn diagrams that show overlap among histone modifications that change in modification intensity after mating of young flies (green circle) and aged flies (red circle), and as a result of aging in virgin flies (blue circle) and mated flies (yellow circle). (a) H3K4me1 modification. (b) H3K4me3 modification. (c) H3K9ac modification. (TIFF 998 KB)
Additional file 7: Changes in H3K4me1 modification. Table S4. Changes in H3K4me1 modification between young and aged female flies. Table S5. Changes in H3K4me1 modification between virgin and mated young female flies. (XLSX 298 KB)
Additional file 8: Changes in H3K4me3 modification. Table S6. Changes in H3K4me3 modification between young and aged female flies. Table S7. Changes in H3K4me3 modification between virgin and mated female flies. (XLSX 610 KB)
Additional file 9: Changes in H3K9ac modification. Table S8. (a) Changes in H3K9ac modification between young and aged virgin female flies. (b) Changes in H3K9ac modification between young and aged mated female flies. Table S9. (a) Changes in H3K9ac modification between virgin and mated young female flies. (b) Changes in H3K9ac modification between virgin and mated aged female flies. (XLSX 715 KB)
Additional file 10: Figure S4: Three-way Venn diagrams that show overlap among H3K4me1 (blue circle), H3K4me3 (green circle) and H3K9ac (yellow circle) modifications that change in modification intensity after mating of young flies (a), and aged flies (b), and as a result of aging in virgin flies (c) and mated flies (d). (TIFF 832 KB)
Additional file 11: Histone modifications and associated genes that are response for mating-induced accelerated aging. Table S10. Genes associated with H3K4me1 modification changes. Table S11. Genes associated with H3K4me3 modification changes. Table S12. Genes associated with H3K9ac modification changes. (XLSX 49 KB)
Additional file 12: Genes identified from the DroID database that can interact with miRNA target genes which change expression after mating or upon aging. Table S13. Genes that can interact with miRNA target genes which change expression in young flies after mating. Table S14. Genes that can interact with miRNA target genes which change expression in virgin flies upon aging. Table S15. Genes that can interact with miRNA target genes which change expression in mated flies upon aging. (XLSX 76 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Zhou, S., Mackay, T.F. & Anholt, R.R. Transcriptional and epigenetic responses to mating and aging in Drosophila melanogaster. BMC Genomics 15, 927 (2014). https://doi.org/10.1186/1471-2164-15-927