Integrative Analysis of Differentially Expressed Genes in Time-Course Multi-Omics Data with MINT-DE.

Background: Time-course multi-omics experiments have been highly informative for obtaining a comprehensive understanding of the dynamic relationships between molecules in a biological process, especially if the different profiles are obtained from the same samples. A fundamental step in analyzing time-course multi-omics data involves selecting a short list of genes or gene regions ("sites") that warrant further study. Two important criteria for site selection are the magnitude of change and the temporal dynamic consistency. However, existing methods only consider one of these criteria, while neglecting the other. Results: In our study, we propose a framework called MINT-DE (Multi-omics INtegration of Time-course for Diffierential Expression analysis) to address this limitation. MINT-DE is capable of selecting sites based on summarized measures of both aforementioned aspects. We calculate evidence measures assessing the extent of differential expression for each assay and for the dynamical similarity across assays. Then based on the summary of the evidence assessment measures, sites are ranked. To evaluate the performance of MINT-DE, we apply it to analyze a time-course multi-omics dataset of Drosophila development. We compare the selection obtained from MINT-DE with those obtained from other existing methods. The analysis reveal that MINT-DE is able to identify differentially expressed time-course pairs with the highest correlations. Their corresponding genes are significantly enriched for known biological functions, as measured by gene-gene interaction networks and the Gene Ontology enrichment. Conclusions: These findings suggest the effectiveness of MINT-DE in selecting sites that are both differentially expressed within at least one assay and temporally related across assays. This highlights the potential of MINT-DE to identify biologically important sites for downstream analysis and provide a complementarity of sites that is neglected by existing methods.

Transcriptional programs are activated and microRNAs are repressed within minutes after mating in the Drosophila melanogaster female reproductive tract.

BACKGROUND: The female reproductive tract is exposed directly to the male's ejaculate, making it a hotspot for mating-induced responses. In Drosophila melanogaster, changes in the reproductive tract are essential to optimize fertility. Many changes occur within minutes after mating, but such early timepoints are absent from published RNA-seq studies. We measured transcript abundances using RNA-seq and microRNA-seq of reproductive tracts of unmated and mated females collected at 10-15 min post-mating. We further investigated whether early transcriptome changes in the female reproductive tract are influenced by inhibiting BMPs in secondary cells, a condition that depletes exosomes from the male's ejaculate. RESULTS: We identified 327 differentially expressed genes. These were mostly upregulated post-mating and have roles in tissue morphogenesis, wound healing, and metabolism. Differentially abundant microRNAs were mostly downregulated post-mating. We identified 130 predicted targets of these microRNAs among the differentially expressed genes. We saw no detectable effect of BMP inhibition in secondary cells on transcript levels in the female reproductive tract. CONCLUSIONS: Our results indicate that mating induces early changes in the female reproductive tract primarily through upregulation of target genes, rather than repression. The upregulation of certain target genes might be mediated by the mating-induced downregulation of microRNAs. Male-derived exosomes and other BMP-dependent products were not uniquely essential for this process. Differentially expressed genes and microRNAs provide candidates that can be further examined for their participation in the earliest alterations of the reproductive tract microenvironment.

Time series transcriptome analysis implicates the circadian clock in the Drosophila melanogaster female's response to sex peptide.

Sex peptide (SP), a seminal fluid protein of Drosophila melanogaster males, has been described as driving a virgin-to-mated switch in females, through eliciting an array of responses including increased egg laying, activity, and food intake and a decreased remating rate. While it is known that SP achieves this, at least in part, by altering neuronal signaling in females, the genetic architecture and temporal dynamics of the female's response to SP remain elusive. We used a high-resolution time series RNA-sequencing dataset of female heads at 10 time points within the first 24 h after mating to learn about the genetic architecture, at the gene and exon levels, of the female's response to SP. We find that SP is not essential to trigger early aspects of a virgin-to-mated transcriptional switch, which includes changes in a metabolic gene regulatory network. However, SP is needed to maintain and diversify metabolic changes and to trigger changes in a neuronal gene regulatory network. We further find that SP alters rhythmic gene expression in females and suggests that SP's disruption of the female's circadian rhythm might be key to its widespread effects.

Dense time-course gene expression profiling of the Drosophila melanogaster innate immune response.

BACKGROUND: Immune responses need to be initiated rapidly, and maintained as needed, to prevent establishment and growth of infections. At the same time, resources need to be balanced with other physiological processes. On the level of transcription, studies have shown that this balancing act is reflected in tight control of the initiation kinetics and shutdown dynamics of specific immune genes. RESULTS: To investigate genome-wide expression dynamics and trade-offs after infection at a high temporal resolution, we performed an RNA-seq time course on D. melanogaster with 20 time points post Imd stimulation. A combination of methods, including spline fitting, cluster analysis, and Granger causality inference, allowed detailed dissection of expression profiles, lead-lag interactions, and functional annotation of genes through guilt-by-association. We identified Imd-responsive genes and co-expressed, less well characterized genes, with an immediate-early response and sustained up-regulation up to 5 days after stimulation. In contrast, stress response and Toll-responsive genes, among which were Bomanins, demonstrated early and transient responses. We further observed a strong trade-off with metabolic genes, which strikingly recovered to pre-infection levels before the immune response was fully resolved. CONCLUSIONS: This high-dimensional dataset enabled the comprehensive study of immune response dynamics through the parallel application of multiple temporal data analysis methods. The well annotated data set should also serve as a useful resource for further investigation of the D. melanogaster innate immune response, and for the development of methods for analysis of a post-stress transcriptional response time-series at whole-genome scale.

Identification of a micropeptide and multiple secondary cell genes that modulate Drosophila male reproductive success.

Even in well-characterized genomes, many transcripts are considered noncoding RNAs (ncRNAs) simply due to the absence of large open reading frames (ORFs). However, it is now becoming clear that many small ORFs (smORFs) produce peptides with important biological functions. In the process of characterizing the ribosome-bound transcriptome of an important cell type of the seminal fluid-producing accessory gland of Drosophila melanogaster, we detected an RNA, previously thought to be noncoding, called male-specific abdominal (msa). Notably, msa is nested in the HOX gene cluster of the Bithorax complex and is known to contain a micro-RNA within one of its introns. We find that this RNA encodes a "micropeptide" (9 or 20 amino acids, MSAmiP) that is expressed exclusively in the secondary cells of the male accessory gland, where it seems to accumulate in nuclei. Importantly, loss of function of this micropeptide causes defects in sperm competition. In addition to bringing insights into the biology of a rare cell type, this work underlines the importance of small peptides, a class of molecules that is now emerging as important actors in complex biological processes.

Interactions between the microbiome and mating influence the female's transcriptional profile in Drosophila melanogaster.

Drosophila melanogaster females undergo a variety of post-mating changes that influence their activity, feeding behavior, metabolism, egg production and gene expression. These changes are induced either by mating itself or by sperm or seminal fluid proteins. In addition, studies have shown that axenic females-those lacking a microbiome-have altered fecundity compared to females with a microbiome, and that the microbiome of the female's mate can influence reproductive success. However, the extent to which post-mating changes in transcript abundance are affected by microbiome state is not well-characterized. Here we investigated fecundity and the post-mating transcript abundance profile of axenic or control females after mating with either axenic or control males. We observed interactions between the female's microbiome and her mating status: transcripts of genes involved in reproduction and genes with neuronal functions were differentially abundant depending on the females' microbiome status, but only in mated females. In addition, immunity genes showed varied responses to either the microbiome, mating, or a combination of those two factors. We further observed that the male's microbiome status influences the fecundity of both control and axenic females, while only influencing the transcriptional profile of axenic females. Our results indicate that the microbiome plays a vital role in the post-mating switch of the female's transcriptome.

Female Genetic Contributions to Sperm Competition in Drosophila melanogaster.

In many species, sperm can remain viable in the reproductive tract of a female well beyond the typical interval to remating. This creates an opportunity for sperm from different males to compete for oocyte fertilization inside the female's reproductive tract. In Drosophila melanogaster, sperm characteristics and seminal fluid content affect male success in sperm competition. On the other hand, although genome-wide association studies (GWAS) have demonstrated that female genotype plays a role in sperm competition outcome as well, the biochemical, sensory, and physiological processes by which females detect and selectively use sperm from different males remain elusive. Here, we functionally tested 26 candidate genes implicated via a GWAS for their contribution to the female's role in sperm competition, measured as changes in the relative success of the first male to mate (P1). Of these 26 candidates, we identified eight genes that affect P1 when knocked down in females, and showed that five of them do so when knocked down in the female nervous system. In particular, Rim knockdown in sensory pickpocket (ppk)+ neurons lowered P1, confirming previously published results, and a novel candidate, caup, lowered P1 when knocked down in octopaminergic Tdc2+ neurons. These results demonstrate that specific neurons in the female's nervous system play a functional role in sperm competition and expand our understanding of the genetic, neuronal, and mechanistic basis of female responses to multiple matings. We propose that these neurons in females are used to sense, and integrate, signals from courtship or ejaculates, to modulate sperm competition outcome accordingly.

Roles of Female and Male Genotype in Post-Mating Responses in Drosophila melanogaster.

Mating induces a multitude of changes in female behavior, physiology, and gene expression. Interactions between female and male genotype lead to variation in post-mating phenotypes and reproductive success. So far, few female molecules responsible for these interactions have been identified. Here, we used Drosophila melanogaster from 5 geographically dispersed populations to investigate such female × male genotypic interactions at the female transcriptomic and phenotypic levels. Females from each line were singly-mated to males from the same 5 lines, for a total of 25 combinations. Reproductive output and refractoriness to re-mating were assayed in females from the 25 mating combinations. Female × male genotypic interactions resulted in significant differences in these post-mating phenotypes. To assess whether female × male genotypic interactions affect the female post-mating transcriptome, next-generation RNA sequencing was performed on virgin and mated females at 5 to 6 h post-mating. Seventy-seven genes showed strong variation in mating-induced expression changes in a female × male genotype-dependent manner. These genes were enriched for immune response and odorant-binding functions, and for expression exclusively in the head. Strikingly, variation in post-mating transcript levels of a gene encoding a spermathecal endopeptidase was correlated with short-term egg production. The transcriptional variation found in specific functional classes of genes might be a read-out of female × male compatibility at a molecular level. Understanding the roles these genes play in the female post-mating response will be crucial to better understand the evolution of post-mating responses and related conflicts between the sexes.