Transcriptional and epigenetic responses to mating and aging in Drosophila melanogaster
© Zhou et al.; licensee BioMed Central Ltd. 2014
Received: 24 July 2014
Accepted: 13 October 2014
Published: 23 October 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.
KeywordsmicroRNA Histone modification ChipSeq Environmental plasticity Systems genetics
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.
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 .
Number of genes associated with changes in histone modifications*
Post-mating changes in young flies
Post-mating changes in aged flies
Changes associated with aging in virgin flies
Changes associated with aging in mated flies
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.
Numbers of candidate histone regulated genes with altered expression after mating or upon aging*
Virgin vs. Mated
Young vs. Aged
Young vs. Aged
Virgin vs. Mated
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.
Effects of micro-RNA overexpression on target gene expression
Real time PCR confirmation of interactions between miRNAs and target genes*
Expression in virgin vs. mated female flies ( P-Value)
miRNA overexpression cross
Effect on target gene expression
Effects of micro-RNA overexpression on egg laying
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 (w 1118 : 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.
We thank Rachel Lyman, Sarah McAdam and Genevieve Evans for technical assistance and Dr. Wen Huang for valuable discussions.
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Bushati N, Cohen SM: MicroRNA functions. Annu Rev Cell Dev Biol. 2007, 23: 175-205. 10.1146/annurev.cellbio.23.090506.123406.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Wolfner MF: Battle and ballet: molecular interactions between the sexes in Drosophila. J Hered. 2009, 100 (4): 399-410. 10.1093/jhered/esp013.PubMed CentralPubMedView ArticleGoogle Scholar
- Chapman T: Seminal fluid-mediated fitness traits in Drosophila. Heredity (Edinb). 2001, 87 (Pt 5): 511-521.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- Tadros W, Lipshitz HD: The maternal-to-zygotic transition: a play in two acts. Development. 2009, 136 (18): 3033-3042. 10.1242/dev.033183.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Lyko F, Ramsahoye BH, Jaenisch R: DNA methylation in Drosophila melanogaster. Nature. 2000, 408 (6812): 538-540. 10.1038/35046205.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Kouzarides T: SnapShot: Histone-modifying enzymes. Cell. 2007, 131 (4): 822-10.1016/j.cell.2007.11.005.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Brand AH, Perrimon N: Targeted gene expression as a means of altering cell fates and generating dominant phenotypes. Development. 1993, 118 (2): 401-415.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9 (4): 357-359. 10.1038/nmeth.1923.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.