Next Article in Journal
Small RNAs of Haloferax mediterranei: Identification and Potential Involvement in Nitrogen Metabolism
Previous Article in Journal
Pathogenicity Islands Distribution in Non-O157 Shiga Toxin-Producing Escherichia coli (STEC)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 9
Issue 2
10.3390/genes9020082
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Re-Arrangements in the Cytoplasmic Distribution of Small RNAs Following the Maternal-to-Zygotic Transition in Drosophila Embryos

by
Mehmet Ilyas Cosacak
1,2,3,
Hatice Yiğit
1,†,
Caghan Kizil
2,3 and
Bünyamin Akgül
1,*
1
Department of Molecular Biology and Genetics, İzmir Institute of Technology, Gülbahçeköyü, 35430 İzmir, Turkey
2
German Center for Neurodegenerative Diseases (DZNE) Dresden, Helmholtz Association, Arnoldstr. 18, 01307 Dresden, Germany
3
Center for Regenerative Therapies Dresden (CRTD), Technische Universität Dresden, Fetscherstr. 105, 01307 Dresden, Germany
*
Author to whom correspondence should be addressed.
†
Present address: Institute for Developmental Biology and Neurobiology, Johannes Gutenberg-University Mainz, Bentzelweg 3, 55128 Mainz, Germany.
Genes 2018, 9(2), 82; https://doi.org/10.3390/genes9020082
Submission received: 30 November 2017 / Revised: 31 January 2018 / Accepted: 2 February 2018 / Published: 10 February 2018
(This article belongs to the Section Molecular Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Small ribonucleic acids (RNAs) are known to regulate gene expression during early development. However, the dynamics of interaction between small RNAs and polysomes during this process is largely unknown. To investigate this phenomenon, 0–1 h and 7–8 h Drosophila melanogaster embryos were fractionated on sucrose density gradients into four fractions based on A254 reading (1) translationally inactive messenger ribonucleoprotein (mRNP), (2) 60S, (3) monosome, and (4) polysome. Comparative analysis of deep-sequencing reads from fractionated and un-fractionated 0–1 h and 7–8 h embryos revealed development-specific co-sedimentation pattern of small RNAs with the cellular translation machinery. Although most micro RNAs (miRNAs) did not have a specific preference for any state of the translational machinery, we detected fraction-specific enrichment of a few miRNAs such as dme-miR-1-3p, -184-3p, 5-5p and 263-5p. More interestingly, we observed changes in the subcellular location of a subset of miRNAs in fractionated embryos despite no measurable difference in their amount in unfractionated embryos. Transposon-derived endo small interfering RNAs (siRNAs) were over-expressed in 7–8 h embryos and associated mainly with the mRNP fraction. In contrast, transposon-derived PIWI-interacting RNAs (piRNA), which were more abundant in 0–1 h embryos, co-sedimented primarily with the polysome fractions. These results suggest that there appears to be a complex interplay among the small RNAs with respect to their polysome-cosedimentation pattern during early development in Drosophila melanogaster.

1. Introduction

In Drosophila melanogaster, early developmental processes until the maternal-to-zygotic transition (MZT) require the function of nearly 7000 maternally-loaded RNAs [1,2]. The MZT transfers the developmental control to the zygotic genome and induces the degradation of maternal messenger RNAs (mRNAs). Accumulating evidence suggests that small non-coding RNAs (ncRNAs) play an essential role in regulating the MZT both in insects and mammals [3,4,5,6,7].
There are at least three types of well-documented ncRNAs, miRNAs, endogenous silencing RNAs (endo-siRNAs) and piRNAs [8,9,10]. Micro RNAs are ubiquitously expressed small RNAs of ~22 nucleotide in length that regulate gene expression post-transcriptionally by decreasing RNA stability or repressing translation at the initiation, 60S joining or post-initiation steps [11,12,13,14]. PIWI-interacting RNAs appear to stabilize the germ cell genome and induce deadenylation of maternal mRNAs in the early Drosophila melanogaster embryo [7,15]. An abundance of sense and antisense endo-siRNAs are generated during transcription from different sites around the promoter and termination regions or from sense mRNAs and longer antisense transcripts [16]. Endo-siRNAs can have various regulatory functions ranging from heterochromatin formation to post-transcriptional gene regulation [8].
There appears to be a complex interplay among the three major classes of small ncRNAs in regulating the MZT. The piRNA pathway has been reported to modulate maternal mRNA deadenylation and decay in the early Drosophila melanogaster embryos [7,17]. Both maternal and sperm-borne endo-siRNAs regulate fertilization and early developmental gene regulatory processes [18,19,20]. Substantial changes have been observed in the miRNA expression profiles of embryos at different developmental stages in Drosophila melanogaster [15,21,22] and targeted knockout mutations result in various phenotypes [23]. The similarity in the miRNA expression pattern in mature oocyte and zygote suggests the presence of maternally inherited miRNAs [24]. The significance of maternal miRNAs pre-MZT is largely unknown, except for dme-miR-34, which is maternally inherited and important for neurogenesis in Drosophila melanogaster [25]. It is proposed that endo-siRNAs and miRNAs may be responsible for regulating early developmental processes before and after the MZT, respectively [26].
Coordinated regulation of mRNAs and miRNAs through the MZT stage is crucial for proper development. In Drosophila melanogaster, most maternal RNA clearance is mediated by the RNA-binding protein Smaug [2]. Dme-miR-309 plays a vital role in post-MZT elimination of maternal RNAs [5]. Interestingly, maternal and zygotic miRNA expressions are coordinated as well. For example, Zelda, in conjunction with maternal morphogens and other zygotic transcription factors, trans-activates zygotic miRNA transcription [27] while Wispy is responsible for adenylation-dependent degradation of maternal miRNAs [28].
Despite a great progress in small RNA-mediated gene regulation, the control of small RNAs themselves has only been tackled recently [29,30]. For instance, miRNA biogenesis can be regulated at multiple steps, including transcription, post-transcriptional processing, editing and intracytoplasmic localization. Apparently, there is a correlation between intracytoplasmic location of (miRNA-induced silencing complex (miRISC)) complexes and their function. Whereas, P-body-associated miRNAs may be involved in mRNA degradation [31], polysome-association of miRNAs may be correlated with miRNA-mediated translational regulation [32]. Interestingly, polysome occupancy of miRNAs was shown to differ between human embryonic stem cells and foreskin fibroblast cells [33]. Differential polysome occupancy, which is apparently influenced by the choice of seed not the abundance, is correlated with the target sequence. Even more interestingly, differential polysome association results in the formation of diverse miRNA effector complexes that are regulated by extracellular signals [34].
We exploited the presence of extensive post-transcriptional gene regulatory networks that occur in Drosophila melanogaster embryos to investigate intracellular dynamics of small RNAs pre- and post-MZT by deep-sequencing. Comparative analysis of fractionated and unfractionated 0–1 h and 7–8 h embryo cytosolic extracts showed that polysomal and non-polysomal fraction-associated small RNA profiles change dramatically following the MZT. Transposon-derived and intergenic-region-derived small RNAs are more abundant in 0–1 h embryos while miRNAs are more abundant in 7–8 h embryos in unfractionated embryos. Micro RNAs appear to interact with the translational machinery at all states, suggesting that each miRISC resides in distinct cytoplasmic reservoirs. Of the two types of transposon-derived small RNAs, siRNAs are expressed at 7–8 h embryos and are primarily associated with the non-polysomal messenger ribonucleoprotein (mRNP) complexes. PIWI-interacting RNAs, on the other hand, are detected more abundantly at 0–1 h embryos and relatively co-localize with polysomal complexes. Altogether these results suggest that each type of small RNAs possesses a different polysome association profile and their profile appears to be regulated during the MZT stage in Drosophila melanogaster.

2. Materials and Methods

2.1. Small RNA Deep-Sequencing Data Analysis

Embryo collection, polysome profiling and small RNA deep-sequencing data have been previously described ([35], GEO accession number GSE35443). We used the small RNA-seq data to investigate the small RNA dynamics during early development in Drosophila melanogaster. To this extent, the 36-bp multiplex sequences were split into their corresponding samples according to their barcodes. Only sequences without any ambiguity were used for further analyses.
After barcode selection, the 3′ adapter sequence was trimmed from the raw reads in 4 steps by using the 3′ adapter sequence 5′-ATCTCGTATGCCGTCTTCTGCTTGT-3′: (1) Trimming full adapter sequences identified inserts with a minimum length of 7 nt, (2) If no adapter sequences were found, the last base of the 3′ adapter was removed in a stepwise manner down to a minimum adapter size of 4 nt, which permit identifying inserts up to 8–28 nt, (3) Finally, the remaining reads were searched for adapter sequences with mismatches, and (4) Inserts with poly(A) and/or multiple unidentifiable nucleotides (Ns) were removed from the data. The resulting 15–29-bp inserts were collapsed into a fasta format and used for subsequent analyses.
For quality control analyses, nexalign program [36] was used to align all sequences to the genome (dmel-r6.18, flybase.org) [37] and CCA-appended mature transfer RNAs (tRNAs) (flybase.org). The sequences were categorized in an order as those with exact match, one (insertion or deletion or mismatch), two or three mismatches. The remaining sequences were grouped as unmapped.
All known RNA sequences were downloaded from flybase (dmel-r6.18) except (1) hairpin and miRNAs sequences from mirbase (www.mirbase.org) (Release 21) [38], (2) ribosomal RNA (rRNA) sequences (5.8S, 18S and 28S) from GenBank (M21017.1) [39], (3) Repbase collection (RepBase22.12) from Jurka et al. 2005 [40], and (4) piRNA cluster genomic coordinates from Brennecke et al. in 2007 [41]. The piRNA clusters from sense strand were extracted by an in-house-algorithm from dmel-r5.39. The 15–29-bp deep-sequencing reads were aligned to these known Drosophila melanogaster RNAs in the following order: rRNA, miRNA hairpin, tRNA, miscellaneous RNA (miscRNA), ncRNA, transposon, transcript, intron, pseudogene and intergenic region. The initial alignment was carried out for exact matches followed by alignment with one, two or three mismatches. In order to analyse the piRNA cluster and Repbase collection, the transposon-mapped sequences were remapped for exact matches to the repbase collection and the intergenic-region- and transposon-mapped reads were realigned to the piRNA clusters.
To calculate miRNA expression levels, we used the exact-matched mature miRNAs sequences and hairpin sequences in miRBase (Release 21). The fold-change for each miRNA was calculated using the formula (fold change = ((7–8h_RPM + 10)/(0–1h_RPM + 10)) in which 10 reads were added to each read to eliminate overestimations in lower reads. The resulting miRNAs were clustered using Cluster 3.0 [42] and visualized by java Treeview [43]. The transposon sequences from Flybase were used to identify the transposon-derived piRNA expression levels. We aligned all transposon-derived 23–29-nt reads to the Repbase collection to obtain the log2 ratio of 0–1 h and 7–8 h RPM values. In order to validate that the 23–29-nt reads are indeed piRNAs, we searched for two features in these reads (1) a 10-nt 5′-5′ complementarity and (2) a preference for U and A at the nucleotide positions 1 and 10, respectively. To analyse the piRNA clusters, the previously published piRNA cluster coordinates were used [41]. The 21-nt reads were not used for piRNA analysis, as they represent siRNA population. All scripts for small RNA sequence data analyses were written in Python and available upon request.
We downloaded 10 fastq files for 0–2 h and 11 fastq files for 6–8 h embryonic transcripts from modENCODE [44]. All fastq files were aligned to Drosophila melanogaster genome (dmel-r6.18) using GSNAP [45]. featureCounts [46,47] was used to assign reads to each gene using GTF file (dmel-all-r6.18.gtf) downloaded from flybase.org. We used R package DESeq2 [48] to normalize reads, calculate fold changes and p-adjusted values. Differentially expressed genes were selected using 2-fold change and p adjusted value <0.1. The pipeline for RNA-seq of mRNA data analyses were written in R and available upon request.

2.2. Quantitative PCR Analysis

RNAs smaller than 200 nt were isolated from unfractionated 0–1 h and 7–8 h embryos using the miRVana miRNA isolation kit according to the manufacturer’s instruction (Thermo Fisher Scientific, Waltham, MA, USA). Complementary DNA was prepared from 200 ng small RNAs using RT2 miRNA cDNA kit (QIAGEN, Hilden, Germany). qPCR was performed in duplicates of three biological replicates (Lightcycler 480, Roche, Basel, Switzerland) using the 3′ standard reverse primer (5′TAGGT22-3′) and the 5′ forward primer (mature miRNA sequence). Relative miRNA expression was calculated after normalizing the data with U6 ncRNA.

3. Results

Blastoderm cellularization, which takes place between 2.5 and 3 h after fertilization in Drosophila melanogaster, is the first developmental process that requires post-MZT zygotic transcription [2]. Small RNAs play an important role in the regulation of the MZT to set the stage for post-MZT embryonic events [3,4,10]. To gain insight into the expression levels of small RNAs during MZT, we analyzed the small RNA-seq data that was previously deposited to GEO ([35], GSE35443). This data set includes small RNA-seq of 0–1 h and 7–8 h embryos. Half of these embryos were used to isolate total RNAs from unfractionated embryos (Figure 1A, un-fractionated (UF)) as a reference. The other halves were then fractionated, on sucrose density gradients (SDGs), into 4 major sub-fractions based on A254 absorbance: (1) mRNPs devoid of rRNAs, (2) 60S, containing 28S rRNA, (3) monosome, and (4) polysome. Here we report the comparative analysis of dynamic changes that we have observed in the intracytoplasmic localization of all small RNAs during early development in Drosophila melanogaster.
One major concern in polysome profiling is the potential for a greater RNA degradation during sample processing prior to downstream events. Thus, we first checked the RNA quality of total RNAs and small RNAs using Agilent RNA 6000 Nano kit and Agilent Small RNA kit, respectively (Agilent, Santa Clara, CA, USA) (Figure S1). Total RNA capillary electrophoresis verified high quality of RNAs phenol-extracted from fractions [35]. Interestingly, each fraction appears to possess a unique small RNA profile especially based on small RNAs that are smaller than tRNAs. When we calculated the frequency of deep-sequencing reads based on the read size, a great majority (88.8–99.5%) of adaptor-containing reads contained inserts of 15–29 bp in size, further verifying a high quality of input RNA. The low percentage of 0–14-bp inserts in fractionated RNAs, compared to that in unfractionated total RNAs, was consistent with the high quality of fractionated RNAs (Table S1). Interestingly, 1 h total RNA contained the greatest percentage of 0–14 bp inserts, congruous with the destabilization of maternal mRNAs [2]. Of the 12,553,921 total reads, 71.95% of the fragments matched perfectly to the Drosophila melanogaster genome with a number of unique sequences in each sample ranging from 71,830 to 757,803. Based on the relative length distribution, we observed enrichment in 22–28-bp RNAs (Figure S2).

3.1. Small RNA Populations Localize in Distinct Cytoplasmic Reservoirs

We used the nexalign program to align sense or antisense sequences to the known RNAs [36]. A pre-alignment to all known RNAs prompted us to align sequences in an order as described in the Methods section. We obtained a perfect match ratio of 71–82% whereas more than 90% of sequences mapped with exact or a single mismatch. Quite interestingly, each small RNA population appeared to co-sediment with polysomal fractions to a different extent (Figure 1B). In the non-polysomal mRNP fraction, the majority of small RNAs contained tRNA-derived fragments (tRFs), with 31% and 37% in 0–1 h and 7–8 h, respectively. Based on northern blot results, we detected differentially expressed tRFs in accordance with the deep-sequencing data, which was published previously [35]. In the polysomal fraction, the majority of small RNAs was derived from rRNAs and transposons. The pattern of enrichment with respect to the developmental stage (1 h versus 8 h) was similar in fractionated and unfractionated samples, indicating the potential biological significance rather than random degradation introduced during fractionation.
A relatively high percentage of rRNA-derived fragments in 80S and polysomal fractions of both 0–1 h and 7–8 h embryos is conceivable considering the abundance of ribosomes in these fractions (Figure 1B). We noticed a much higher rRNA-derived fragment in the total RNAs of 7–8 h embryos compared to 0–1 h embryos (Figure 1B) despite quite similar RNA quality based on Bioanalyzer capillary electrophoresis (Figure S1A). Interestingly, while the percentage of rRNA-derived fragments are similar in mRNP, 60S and 80S fractions, it is higher in polysomal fractions in 7–8 h embryos. Normally, degraded fragments would be expected to enrich in mRNP fractions, however these small RNAs co-sediment mainly with the polysomal fraction, suggesting a biological degradation or processing from rRNAs rather than random degradation. Additionally, the transcript-derived fragments are more abundant in 0–1 h embryos, which probably indicates either clearance of maternal mRNAs or biogenesis of novel zygotic small RNAs. We then focused on miRNAs and transposon-derived small RNAs as they constitute the bulk of small RNAs in our samples.

3.2. Micro RNAs Interact with Cellular Translational Machinery at All States

Although the temporal expression of miRNAs during early development is well documented [15], the extent to which each dysregulated miRNA is associated with polysomes is unknown. To this end, we firstly aligned the RNA-seq reads to the miRNA hairpin with perfect matches. A total of 263 mature miRNAs were detected with at least 1 read in at least one sample. Setting the minimal threshold at 50 RPM filtered out 169 miRNAs, yielding a total of 94 mature miRNAs for further analyses.
To identify fraction-specific miRNAs, we selected ten most abundant hairpin sequences from each sample (10 hairpins from each sample of total, mRNP, 60S, monosome and polysome, yielding a total of 50 hairpins) and identified 18 common hairpin sequences that constituted 82.5–93.6% of all hairpin sequences with perfect matches in each sample (Table 1). Thus, these 18 miRNAs represent the most abundant miRNAs in each sample. We observed dynamic changes in the expression of many miRNAs in unfractionated total RNAs in agreement with the published results (Table 2 [15]). Analysis of miRNA expression in fractionated samples yielded several interesting points (Table 1). Firstly, for the majority of miRNAs, the miRNA levels in fractionated samples were comparable to those in the unfractionated samples, probably suggesting non-selective distribution of these miRNAs throughout the cytoplasm. Secondly, although most miRNAs are distributed nearly evenly throughout the fractions, a few miRNAs were enriched in particular fractions. For instance, dme-miR-5-5p is enriched in the mRNP fraction both in 0–1 h and 7–8 h embryos whereas dme-miR-1-3p is enriched in the 80S fraction. Thirdly and more importantly, the extent to which a miRNA is associated with a particular fraction appears to be relatively similar for most miRNAs. However, for a few miRNAs, the degree of association is developmentally regulated. For instance, dme-miR-9c-5p makes up of 11.5 and 4.1% of miRNAs in 0–1 h and 7–8 h embryos, respectively.
It is well known that maternal miRNAs are degraded after the MZT is completed and zygotic miRNA transcription modulates zygotic gene regulatory networks [28]. It is unknown however whether the polysome status of zygotic miRNAs follow that of maternal miRNAs. There are potentially two possibilities (1) if the target sequence determines the extent of association to different components present in each fraction, we would expect differential polysome association pre- and post-MZT, assuming that the targets of miRNAs pre- and post-MZT stage are different, (2) if the association of miRNAs with different components is independent from the target sequences, then we would expect to see a similar polysome-association pattern. To differentiate between these two hypotheses, we first calculated the miRNA expression ratios (0–1 h/7–8 h) in the unfractionated and fractionated samples. We then checked how the ratio in the unfractionated samples is reflected upon that in the fractionated samples. We assumed that if the cytoplasmic fate of a miRNA does not change pre- and post-MZT, the ratios obtained from the unfractionated and fractionated samples should be similar between 0–1 h and 7–8 h embryos. In contrast, if there are any changes in the localization of miRNAs, the difference in the expression level obtained from the unfractionated samples should manifest in a particular sub-cellular fraction. This approach resulted in identification of 4 different miRNA groups with each having a unique expression behaviour.
The first group includes 41 miRNAs that are over-expressed in 7–8 h embryos (Table 2, G1; Figure 2). While 17 miRNAs co-sediment with the complexes in the mRNP fraction, some co-sediment with the 60S fraction. Interestingly, for 9 miRNAs, we did not see a proportional increase in their expression in fractionated RNAs although their expression was increased in un-fractionated embryos, suggesting the involvement of nuclear retention or other unknown mechanisms. The relative number of miRNAs classified in the second group does not change in 0–1 h and 7–8 h un-fractionated embryos, indicating a similar transcriptional activity and/or miRNA stability. However, we detected dynamic changes in their subcellular locations following the MZT (Table 2, G2; Figure 2). For instance, despite no difference in the total dme-miR-9a-5p amount in un-fractionated embryos, this specific miRNA becomes more polysome-associated in 7–8 h embryos (Log2 fold = 3.3). The third group includes 29 miRNAs, whose expression decreases in 7–8 h embryos (Table 2, G3; Figure 2). Interestingly, the decrease in the miRNA expression is not distributed evenly throughout the 4 fractions, suggesting a preference for a specific fraction. For instance, 20 out of these 29 miRNAs sediment specifically with non-polysomal fractions, particularly in the mRNP fraction. This indicates that the majority of small RNAs acts at mRNP complexes in early embryo compared to the 7–8 h ones. For some miRNAs, we detected equal distribution in the cytoplasm despite a decrease in the unfractionated embryos, suggesting perhaps more efficient nucleo-cytoplasmic transport. The fourth group includes the miRNAs whose relative abundance in the un-fractionated and fractionated embryos are different. These miRNAs behave similarly to those in the second group in that the transcriptional output from these miRNAs does not lead to proportional changes in their cytoplasmic location, suggesting differential localization pathways pre- and post-MZT (Table 2, G4; Figure 2). We then selected six miRNAs, representing at least one miRNA per group, to validate RNA-seq data by quantitative PCR (qPCR) analyses of total RNAs from unfractionated embryos (Table 3). Results from both methods were in parallel to each other.

3.3. Micro RNA Cluster Behaviours and miRNA Editing

Some miRNAs are known to be transcribed as members of gene clusters [50]. Expression studies revealed that miRNA clusters are co-expressed [51] and cluster members are coordinated during target regulation [52]. Assuming such coordination during target RNA regulation, we hypothesized that the cytoplasmic fate of the members should be similar. Thus, we compared the extent to which each cluster member is associated with four different fractions in our experimental design, each representing a different translational state of the cell. Our cluster analysis showed that miRNA cluster members behave similarly with respect to their cytoplasmic localization (Table S2).
We also checked the frequency of post-transcriptional miRNA editing events likely to occur during the early development in Drosophila melanogaster as miRNA editing is commonly used in eukaryotes to modulate the targets of miRNAs [53]. We first aligned our sequences to the known miRNA sequences and looked for the sequences that align to the known sequences with a single mismatch. Based on this approach, we identified one candidate editing event (dme-miR-986, C→T at 11th position) (Figure S3). The PCR-amplification and sequencing of dme-miR-986 from P2 strain embryos and S2 cells showed that this particular difference in the sequence stems from a single nucleotide polymorphism (SNP) not an editing process.

3.4. Cytoplasmic Micro RNA Re-Arrangement Does Not Correlate with Target Abundance

Our observation in the cytoplasmic distribution of miRNAs is interesting in the sense that localized miRISC complexes may differentially regulate their target mRNAs. To test this hypothesis, we examined a potential correlation between the four miRNA groups that we identified (Table 2) and the target mRNA abundance in 0–2 h and 6–8 h embryos [44]. Thus, we first checked the quality of 0–2 h and 6–8 h embryo RNA-seq data (Figure S4A,B [44]) and then identified the differentially expressed transcripts (Figure S4C). Of these mRNAs differentially expressed in 0–2 h and 6–8 h embryos, we selected the experimentally validated targets of miRNAs included in our four miRNA groups. When we categorized the target transcripts based on their abundance (i.e., up-regulated, down-regulated, no-change), we did not detect a correlation between the cytoplasmic distribution of miRNAs and their effect on target mRNA abundance (Figure S5), suggesting perhaps the presence of some other mechanisms.

3.5. Transposon-Derived Small Interfering RNAs and PIWI-Interacting RNAs Interact with Different Complexes

Some piRNAs are maternally deposited into the oocyte and they might be involved in the degradation of maternal mRNAs from the embryo [7]. Also, MILI/MIWI is associated with polysomal fractions [54]. However intra-cytoplasmic distribution and re-arrangements, if any, of piRNAs during early development are unknown. We used the Repbase [40] collection to calculate the small RNAs generated from transposons. The expression of transposon-derived small RNAs decreased towards 7–8 h, suggesting the significance of these small RNAs probably before the MZT. To ensure that we selected the transposon-derived piRNAs, we looked for two main features associated with piRNAs. The first feature is the 10-nucleotide complementarity at their 5′ ends as generated by the ping pong cycle (Figure 3A). The second feature in Drosophila melanogaster is the presence of a U nucleotide at the 1st position and an A nucleotide at the 10th position (Figure 3B; reference [55]). Furthermore, we only used the reads in length of 23–29 bp to select the sequences that match the aforementioned two criteria.
We divided transposon-derived small RNA transcripts into two groups based on their sizes: transposon-derived piRNAs of 23–29 nt and transposon-derived siRNAs of 21 nt [55]. To globally compare the polysome association of transposon-derived piRNAs to that of transposon-derived siRNAs in embryos, we calculated the read frequency of transposon-derived small RNAs. We noticed a much higher peak at the 21 bp reads in 7–8 h embryo total RNAs (Figure 3D) compared with 0–1 h (Figure 3C). Moreover, the 21-bp read frequency is higher in mRNP-associated complexes decreasing towards polysomal fractions. This data suggest that transposon-derived siRNA expression is more abundant in later developmental stages although piRNA expression appears to be slightly higher in 0–1 h embryos. Additionally, the intracellular localization of transposon-derived siRNAs appears to be different from that of transposon-derived piRNAs (Figure 3C,D). Although transposon-derived piRNAs are distributed throughout four fractions nearly equally, transposon-derived siRNAs are predominantly associated with the mRNP fraction. We then checked the frequency of reads from the 42AB cluster as it is one of the best characterized piRNA clusters in Drosophila melanogaster [56]. This analysis showed, in parallel to previous findings, production of transcripts from both strands: Interestingly, there appears to be a strand-bias in polysome-associated transcripts especially in 0–1 h embryos (Figure 4B).
To substantiate our observation that transposon-derived piRNAs are highly expressed earlier compared to the siRNAs, we collected previously published deep-sequencing data from different developmental stages of Drosophila melanogaster [15,41,57]. We used the same strategy as in Figure 4 to trace the temporal expression of siRNAs and piRNAs. Our analysis showed that the siRNA expression is relatively low in early developmental stages but increases significantly during later developmental stages (Figure 5A). Using the same data sets, we also checked the relative abundance of Drosophila melanogaster small RNAs in various developmental stages. This analysis validated the notion that the transposon-derived small RNA abundance drops in later developmental stages while miRNA expression level increases (Figure 5B). Moreover, the 21-nt siRNAs derived from transposon increase while the piRNA levels (23–29 nt) decrease.

4. Discussion

We present here an in-depth analysis of cytoplasmic distribution of small RNAs in early development with respect to their association with the cellular translational machinery. The comparative analysis of the most abundant small RNAs reveals that the spatial location of each type of small RNA is quite different (Figure 1 and Table S3). Interestingly, miRNAs appear to utilize numerous molecular mechanisms as they interact with the translational machinery at all states. piRNAs, which are expressed more abundantly at the pre-MZT stage, are associated with polysomal complexes. In contrast to piRNAs, transposon-derived siRNAs are more abundantly expressed at the post-MZT stage and are primarily found in the mRNP fraction.
During the maternal-to-zygotic transition two crucial events take place: (1) clearance of maternal mRNAs and (2) transcriptional activation of the zygotic genome. The clearance of maternal mRNAs is orchestrated by numerous RNA-Binding proteins, including Smaug [1,2] and small RNAs such as miRNAs, siRNAs and piRNAs [4,7,59,60,61]. Differential expression and specific functions of miRNAs in early development were previously reported [15,22] and our findings are in agreement with these results. For instance, dme-miR-8 and 10 are up-regulated in 7–8 h embryos in total RNAs while dme-miR-310 and 311 are down-regulated.
Although the differential expression of miRNAs during the MZT stage is well documented, the polysome profile of miRNAs during this switch is unknown. Our findings suggest that the majority of miRNAs appear to be associated with distinct fractions and probably with distinct complexes as a result (e.g., mRNP or polysomal complexes) (Table 2). For instance, dme-miR-263-5p, 5-5p and 9c-5p are primarily found in the mRNP fraction whereas dme-miR-184-3p is mainly part of polysomal fractions. Bantam-3p and dme-miR-1-3p are primarily enriched in 60S and 80S fractions, respectively (Table 2). Although we cannot attribute any mechanistic insight into the fraction-specific functions of miRNAs yet, it might be noteworthy to speculate that this specific localization might be associated with association of miRNAs with different components in each fraction. For instance, since mRNP complexes are known to include translationally inactive mRNAs [62], the miRNAs in this fraction might be involved in sequestration of mRNAs (or storage) away from the translational machinery. The miRNAs in the 60S and 80S fractions are likely to interfere with the translation initiation while the polysome-associated miRNAs probably interfere with translation elongation.
One other interesting observation from our findings is that the change in the number of miRNAs in total RNAs (transcriptional or post-transcriptional increases/decreases) is not distributed equally throughout each fraction in the embryos. Some of differentially expressed miRNAs are directed towards the translationally inactive mRNP complexes whereas some others are destined for the translationally active polysomal complexes (Table 2). Although the polysome association of miRNAs has been documented before [63], the switch in their state suggests an interesting layer of regulation. We interpret all these observations to mean that miRNAs have a pre-defined cytoplasmic fate following their nucleo-cytoplasmic transport. When we looked at the validated target abundance of miRNAs categorized under four different groups, we did not see any correlation between the cytoplasmic localization pattern of miRNAs and their target abundance. We might not have detected any correlation for at least three reasons; (1) other regulatory mechanisms, e.g., transcriptional regulation, could be in effect, (2) target RNA stability may not change upon miRNA binding (i.e., translational regulation), (3) effect on target mRNAs could be residual, making it difficult to detect. Thus, further studies with reporter target sequences would be required to precisely trace the actual effect of cytoplasmic miRNA distribution on its target.
For a few miRNAs (Table 2, dme-miR-1-3p, dme-miR-2b-3p, dme-miR-92a-3p, dme-miR-92b-39 and dme-miR-1012-3p) miRISC complexes appear to switch from one translational state to another. For example, mRNP-fraction-associated miR-1-3p switches to polysomal fractions following the MZT. We cannot, however, conclusively claim that miRNAs switch their intracellular location as some of these post-MZT miRNAs could be zygotically transcribed miRNAs. In this situation, there should be additional factor(s), e.g., miRNA- binding proteins that specify the polysome status of zygotically transcribed miRNAs. In either case, our data suggests that the activity of miRNP complexes may be spatially modulated through the differential cytoplasmic localization of these complexes pre- and post-MZT. More direct evidence is required, however, to demonstrate whether the eukaryotic cells utilize intracellular re-localization of miRNPs as a means to modulate miRNA function and target mRNA(s) as a result.
We also detected cytoplasmically localized transposon-derived small RNAs, which are classified as siRNAs (21 nt) or piRNAs (23–29 nt). PIWI-interacting RNAs protect the zygotic genome against infection by retroviruses that can come from the surrounding follicle cells or endogenous transposons within the female germline [4]. PIWI-interacting RNAs function together with the well-characterized Smaug to direct the CCR4 deadenylase to specific mRNAs, thus facilitating maternal mRNA deadenylation and decay [7]. The cytoplasmic localization of piRNAs in early development is consistent with their role in deadenylation of maternal mRNAs. The exact role of endogenous siRNAs in early development is not well-defined. Based on the observation that the mRNA profiles of wild type and Dgcr8 null mouse oocytes were identical, it was proposed that endo-siRNAs, rather than miRNAs, are responsible for the Dicer knockout phenotype observed in mice [64]. Thus, endo-siRNAs and piRNAs are primarily involved in gametogenesis and very early development whereas miRNAs appear to get involved in later developmental stages [65,66]. Our data is consistent with this view in that transposon-derived transcripts are more abundant in 0–1 h embryos (both fractionated and unfractionated, Figure 1B) whereas miRNA expression is up-regulated in 7–8 h embryos. Endogenous siRNAs are well-known for their role in heterochromatin formation in the nucleus [67]. The identification of endo-siRNAs in the cytoplasm points to two possibilities (1) the process of maturation as they are produced from double-stranded RNAs (dsRNAs) in cytoplasm, (2) a potential role for these small RNAs in post-transcriptional gene regulation.
Transposon-derived siRNAs and piRNAs differ from each other in two ways: (1) 21-nt transposon-derived siRNAs are highly expressed in 7–8 h embryos whereas piRNAs abound in 0–1 h embryos (Figure 2 and Figure 3); and (2) siRNAs are mainly associated with the complexes in the mRNP fraction while piRNAs are largely associated with the polysomal complexes. It is quite interesting that both siRNAs and piRNAs are produced from the same transposons, yet they interact with the cellular translational machinery quite differently. Another interesting point is that we detected a strand-bias in the reads obtained from the dual-strand transposon 42AB in 0–1 h embryos especially in the polysome fraction (Figure 5B). It remains to be investigated whether there is a functional relationship between the biased-production of piRNAs and development.
We previously reported that tRNAs also serve as templates for small RNAs, i.e., ~28-nt tRFs [35]. These small RNAs are expressed at both stages (Figure 1B) and in fact throughout the embryonic development and in mature flies as well. They are mainly associated with non-polysomal fractions, resembling transposon-derived siRNAs. Studies on tRF function mainly focused on miRNA-like regulatory functions. In Drosophila melanogaster, tRFs are immunoprecipitated with anti-AGO1 antibody [68]. Interestingly, they inhibit cap-dependent translation initiation [69]. It is unknown however whether tRFs are specifically involved in modulation of translation at the MZT stage.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/9/2/82/s1. Figure S1: Total RNA quality control analysis, Figure S2: The insert length distribution of small RNAs, Figure S3: SNP on dme-miR-986-5p, Figure S4: Analysis of 0-2 and 6-8 hour old embryonic RNA-Seq data, Figure S5: Correlation between miRNA and their targets [67], Table S1: The percentage of hairpin and mature miRNA mapped reads in all reads mapped to miRNAs, Table S2: miRNA cluster expression analysis, Table S3: Alignment of the reads to known RNAs.

Acknowledgments

We thank C.-P. David Tu of Pennsylvania State University for providing the Drosophila melanogaster P2 strain. We also thank the İYTE Biotechnology Research Center for the instrumental support. This work was supported by the Scientific and Technical Research Council of Turkey [104T144 to BA].

Author Contributions

B.A. conceived and designed the experiments; H.Y. performed the original experiments; M.I.C. performed small RNA-seq and modENCODE RNA-seq data analyses in B.A. and C.K. labs, respectively; B.A. and M.I.C. wrote and revised the paper; H.Y. contributed to the editing of the paper; C.K. read and confirmed.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Renzis, S.; Elemento, O.; Tavazoie, S.; Wieschaus, E.F. Unmasking activation of the zygotic genome using chromosomal deletions in the Drosophila embryo. PLoS Biol. 2007, 5, e117. [Google Scholar] [CrossRef]
  2. Tadros, W.; Goldman, A.L.; Babak, T.; Menzies, F.; Vardy, L.; Orr-Weaver, T.; Hughes, T.R.; Westwood, J.T.; Smibert, C.A.; Lipshitz, H.D. Smaug is a major regulator of maternal mRNA destabilization in Drosophila and its translation is activated by the PAN GU kinase. Dev. Cell 2007, 12, 143–155. [Google Scholar] [CrossRef] [PubMed]
  3. Alberti, C.; Cochella, L. A framework for understanding the roles of miRNAs in animal development. Development 2017, 144, 2548–2559. [Google Scholar] [CrossRef] [PubMed]
  4. Bourc’his, D.; Voinnet, O. A small-RNA perspective on gametogenesis, fertilization, and early zygotic development. Science 2010, 330, 617–622. [Google Scholar] [CrossRef] [PubMed]
  5. Bushati, N.; Stark, A.; Brennecke, J.; Cohen, S.M. Temporal reciprocity of miRNAs and their targets during the maternal-to-zygotic transition in Drosophila. Curr. Biol. 2008, 18, 501–506. [Google Scholar] [CrossRef] [PubMed]
  6. Ma, J.; Flemr, M.; Stein, P.; Berninger, P.; Malik, R.; Zavolan, M.; Svoboda, P.; Schultz, R.M. MicroRNA activity is suppressed in mouse oocytes. Curr. Biol. 2010, 20, 265–270. [Google Scholar] [CrossRef] [PubMed]
  7. Rouget, C.; Papin, C.; Boureux, A.; Meunier, A.C.; Franco, B.; Robine, N.; Lai, E.C.; Pelisson, A.; Simonelig, M. Maternal mRNA deadenylation and decay by the piRNA pathway in the early Drosophila embryo. Nature 2010, 467, 1128–1132. [Google Scholar] [CrossRef] [PubMed]
  8. Carthew, R.W.; Sontheimer, E.J. Origins and mechanisms of miRNAs and siRNAs. Cell 2009, 136, 642–655. [Google Scholar] [CrossRef] [PubMed]
  9. Ghildiyal, M.; Zamore, P.D. Small silencing RNAs: An expanding universe. Nat. Rev. Genet. 2009, 10, 94–108. [Google Scholar] [CrossRef] [PubMed]
  10. Weick, E.M.; Miska, E.A. PiRNAs: From biogenesis to function. Development 2014, 141, 3458–3471. [Google Scholar] [CrossRef] [PubMed]
  11. Chendrimada, T.P.; Finn, K.J.; Ji, X.; Baillat, D.; Gregory, R.I.; Liebhaber, S.A.; Pasquinelli, A.E.; Shiekhattar, R. MicroRNA silencing through risc recruitment of Eif6. Nature 2007, 447, 823–828. [Google Scholar] [CrossRef] [PubMed]
  12. Guo, H.; Ingolia, N.T.; Weissman, J.S.; Bartel, D.P. Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature 2010, 466, 835–840. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Kiriakidou, M.; Tan, G.S.; Lamprinaki, S.; De Planell-Saguer, M.; Nelson, P.T.; Mourelatos, Z. An mRNA m7G cap binding-like motif within human Ago2 represses translation. Cell 2007, 129, 1141–1151. [Google Scholar] [CrossRef] [PubMed]
  14. Petersen, C.P.; Bordeleau, M.E.; Pelletier, J.; Sharp, P.A. Short RNAs repress translation after initiation in mammalian cells. Mol. Cell 2006, 21, 533–542. [Google Scholar] [CrossRef] [PubMed]
  15. Aravin, A.A.; Lagos-Quintana, M.; Yalcin, A.; Zavolan, M.; Marks, D.; Snyder, B.; Gaasterland, T.; Meyer, J.; Tuschl, T. The small RNA profile during Drosophila melanogaster development. Dev. Cell 2003, 5, 337–350. [Google Scholar] [CrossRef]
  16. Okamura, K.; Lai, E.C. Endogenous small interfering RNAs in animals. Nat. Rev. Mol. Cell Biol. 2008, 9, 673–678. [Google Scholar] [CrossRef] [PubMed]
  17. Barckmann, B.; Pierson, S.; Dufourt, J.; Papin, C.; Armenise, C.; Port, F.; Grentzinger, T.; Chambeyron, S.; Baronian, G.; Desvignes, J.P.; et al. Aubergine iCLIP reveals piRNA-dependent decay of mRNAs involved in germ cell development in the early embryo. Cell Rep. 2015, 12, 1205–1216. [Google Scholar] [CrossRef] [PubMed]
  18. Claycomb, J.M. Ancient endo-siRNA pathways reveal new tricks. Curr. Biol. 2014, 24, R703–R715. [Google Scholar] [CrossRef] [PubMed]
  19. Yuan, S.; Schuster, A.; Tang, C.; Yu, T.; Ortogero, N.; Bao, J.; Zheng, H.; Yan, W. Sperm-borne miRNAs and endo-siRNAs are important for fertilization and preimplantation embryonic development. Development 2016, 143, 635–647. [Google Scholar] [CrossRef] [PubMed]
  20. Zhang, H.; Liu, J.; Tai, Y.; Zhang, X.; Zhang, J.; Liu, S.; Lv, J.; Liu, Z.; Kong, Q. Identification and characterization of l1-specific endo-siRNAs essential for early embryonic development in pig. Oncotarget 2017, 8, 23167–23176. [Google Scholar] [CrossRef] [PubMed]
  21. Aboobaker, A.A.; Tomancak, P.; Patel, N.; Rubin, G.M.; Lai, E.C. Drosophila microRNAs exhibit diverse spatial expression patterns during embryonic development. Proc. Natl. Acad. Sci. USA 2005, 102, 18017–18022. [Google Scholar] [CrossRef] [PubMed]
  22. Leaman, D.; Chen, P.Y.; Fak, J.; Yalcin, A.; Pearce, M.; Unnerstall, U.; Marks, D.S.; Sander, C.; Tuschl, T.; Gaul, U. Antisense-mediated depletion reveals essential and specific functions of microRNAs in Drosophila development. Cell 2005, 121, 1097–1108. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, Y.W.; Song, S.; Weng, R.; Verma, P.; Kugler, J.M.; Buescher, M.; Rouam, S.; Cohen, S.M. Systematic study of Drosophila microRNA functions using a collection of targeted knockout mutations. Dev. Cell 2014, 31, 784–800. [Google Scholar] [CrossRef] [PubMed]
  24. Tang, F.; Kaneda, M.; O’Carroll, D.; Hajkova, P.; Barton, S.C.; Sun, Y.A.; Lee, C.; Tarakhovsky, A.; Lao, K.; Surani, M.A. Maternal microRNAs are essential for mouse zygotic development. Genes Dev. 2007, 21, 644–648. [Google Scholar] [CrossRef] [PubMed]
  25. Soni, K.; Choudhary, A.; Patowary, A.; Singh, A.R.; Bhatia, S.; Sivasubbu, S.; Chandrasekaran, S.; Pillai, B. Mir-34 is maternally inherited in Drosophila melanogaster and Danio rerio. Nucleic Acids Res. 2013, 41, 4470–4480. [Google Scholar] [CrossRef] [PubMed]
  26. Svoboda, P.; Flemr, M. The role of miRNAs and endogenous siRNAs in maternal-to-zygotic reprogramming and the establishment of pluripotency. EMBO Rep. 2010, 11, 590–597. [Google Scholar] [CrossRef] [PubMed]
  27. Fu, S.; Nien, C.Y.; Liang, H.L.; Rushlow, C. Co-activation of microRNAs by Zelda is essential for early Drosophila development. Development 2014, 141, 2108–2118. [Google Scholar] [CrossRef] [PubMed]
  28. Lee, M.; Choi, Y.; Kim, K.; Jin, H.; Lim, J.; Nguyen, T.A.; Yang, J.; Jeong, M.; Giraldez, A.J.; Yang, H.; et al. Adenylation of maternally inherited microRNAs by Wispy. Mol. Cell 2014, 56, 696–707. [Google Scholar] [CrossRef] [PubMed]
  29. Ha, M.; Kim, V.N. Regulation of microRNA biogenesis. Nat. Rev. Mol. Cell Biol. 2014, 15, 509–524. [Google Scholar] [CrossRef] [PubMed]
  30. Watanabe, T.; Lin, H. Posttranscriptional regulation of gene expression by PIWI proteins and piRNAs. Mol. Cell 2014, 56, 18–27. [Google Scholar] [CrossRef] [PubMed]
  31. Liu, J.; Valencia-Sanchez, M.A.; Hannon, G.J.; Parker, R. MicroRNA-dependent localization of targeted mRNAs to mammalian P-bodies. Nat. Cell Biol. 2005, 7, 719–723. [Google Scholar] [CrossRef] [PubMed]
  32. Maroney, P.A.; Yu, Y.; Fisher, J.; Nilsen, T.W. Evidence that microRNAs are associated with translating messenger RNAs in human cells. Nat. Struct. Mol. Biol. 2006, 13, 1102–1107. [Google Scholar] [CrossRef] [PubMed]
  33. Molotski, N.; Soen, Y. Differential association of microRNAs with polysomes reflects distinct strengths of interactions with their mRNA targets. RNA 2012, 18, 1612–1623. [Google Scholar] [CrossRef] [PubMed]
  34. Wu, P.H.; Isaji, M.; Carthew, R.W. Functionally diverse microRNA effector complexes are regulated by extracellular signaling. Mol. Cell 2013, 52, 113–123. [Google Scholar] [CrossRef] [PubMed]
  35. Goktas, C.; Yigit, H.; Cosacak, M.I.; Akgul, B. Differentially expressed tRNA-derived small RNAs co-sediment primarily with non-polysomal fractions in Drosophila. Genes 2017, 8, 333. [Google Scholar] [CrossRef] [PubMed]
  36. De Hoon, M.J.; Taft, R.J.; Hashimoto, T.; Kanamori-Katayama, M.; Kawaji, H.; Kawano, M.; Kishima, M.; Lassmann, T.; Faulkner, G.J.; Mattick, J.S.; et al. Cross-mapping and the identification of editing sites in mature microRNAs in high-throughput sequencing libraries. Genome Res. 2010, 20, 257–264. [Google Scholar] [CrossRef] [PubMed]
  37. Gramates, L.S.; Marygold, S.J.; Santos, G.D.; Urbano, J.M.; Antonazzo, G.; Matthews, B.B.; Rey, A.J.; Tabone, C.J.; Crosby, M.A.; Emmert, D.B.; et al. Flybase at 25: Looking to the future. Nucleic Acids Res. 2017, 45, D663–D671. [Google Scholar] [CrossRef] [PubMed]
  38. Griffiths-Jones, S.; Grocock, R.J.; van Dongen, S.; Bateman, A.; Enright, A.J. miRBase: MicroRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34, D140–D144. [Google Scholar] [CrossRef] [PubMed]
  39. Benson, D.A.; Karsch-Mizrachi, I.; Lipman, D.J.; Ostell, J.; Wheeler, D.L. Genbank. Nucleic Acids Res 2008, 36, D25–D30. [Google Scholar] [CrossRef] [PubMed]
  40. Jurka, J.; Kapitonov, V.V.; Pavlicek, A.; Klonowski, P.; Kohany, O.; Walichiewicz, J. Repbase update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110, 462–467. [Google Scholar] [CrossRef] [PubMed]
  41. Brennecke, J.; Aravin, A.A.; Stark, A.; Dus, M.; Kellis, M.; Sachidanandam, R.; Hannon, G.J. Discrete small RNA-generating loci as master regulators of transposon activity in Drosophila. Cell 2007, 128, 1089–1103. [Google Scholar] [CrossRef] [PubMed]
  42. De Hoon, M.J.; Imoto, S.; Nolan, J.; Miyano, S. Open source clustering software. Bioinformatics 2004, 20, 1453–1454. [Google Scholar] [CrossRef] [PubMed]
  43. Saldanha, A.J. Java Treeview—Extensible visualization of microarray data. Bioinformatics 2004, 20, 3246–3248. [Google Scholar] [CrossRef] [PubMed]
  44. Celniker, S.E.; Dillon, L.A.; Gerstein, M.B.; Gunsalus, K.C.; Henikoff, S.; Karpen, G.H.; Kellis, M.; Lai, E.C.; Lieb, J.D.; MacAlpine, D.M.; et al. Unlocking the secrets of the genome. Nature 2009, 459, 927–930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Wu, T.D.; Nacu, S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics 2010, 26, 873–881. [Google Scholar] [CrossRef] [PubMed]
  46. Liao, Y.; Smyth, G.K.; Shi, W. The subread aligner: Fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013, 41, e108. [Google Scholar] [CrossRef] [PubMed]
  47. Liao, Y.; Smyth, G.K.; Shi, W. featureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef] [PubMed]
  48. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  49. Bor, Y.C.; Swartz, J.; Li, Y.; Coyle, J.; Rekosh, D.; Hammarskjold, M.-L. Northern blot analysis of mRNA from mammalian polyribosomes. Protocol Exchange 2006. [Google Scholar] [CrossRef]
  50. Altuvia, Y.; Landgraf, P.; Lithwick, G.; Elefant, N.; Pfeffer, S.; Aravin, A.; Brownstein, M.J.; Tuschl, T.; Margalit, H. Clustering and conservation patterns of human microRNAs. Nucleic Acids Res. 2005, 33, 2697–2706. [Google Scholar] [CrossRef] [PubMed]
  51. Baskerville, S.; Bartel, D.P. Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. RNA 2005, 11, 241–247. [Google Scholar] [CrossRef] [PubMed]
  52. Wang, J.; Haubrock, M.; Cao, K.M.; Hua, X.; Zhang, C.Y.; Wingender, E.; Li, J. Regulatory coordination of clustered microRNAs based on microRNA-transcription factor regulatory network. BMC Syst. Biol. 2011, 5, 199. [Google Scholar] [CrossRef] [PubMed]
  53. Kawahara, Y.; Zinshteyn, B.; Sethupathy, P.; Iizasa, H.; Hatzigeorgiou, A.G.; Nishikura, K. Redirection of silencing targets by adenosine-to-inosine editing of miRNAs. Science 2007, 315, 1137–1140. [Google Scholar] [CrossRef] [PubMed]
  54. Grivna, S.T.; Pyhtila, B.; Lin, H. MIWI associates with translational machinery and PIWI-interacting RNAs (piRNAs) in regulating spermatogenesis. Proc. Natl. Acad. Sci. USA 2006, 103, 13415–13420. [Google Scholar] [CrossRef] [PubMed]
  55. Ghildiyal, M.; Seitz, H.; Horwich, M.D.; Li, C.; Du, T.; Lee, S.; Xu, J.; Kittler, E.L.; Zapp, M.L.; Weng, Z.; et al. Endogenous siRNAs derived from transposons and mRNAs in Drosophila somatic cells. Science 2008, 320, 1077–1081. [Google Scholar] [CrossRef] [PubMed]
  56. Aravin, A.; Gaidatzis, D.; Pfeffer, S.; Lagos-Quintana, M.; Landgraf, P.; Iovino, N.; Morris, P.; Brownstein, M.J.; Kuramochi-Miyagawa, S.; Nakano, T.; et al. A novel class of small RNAs bind to MILI protein in mouse testes. Nature 2006, 442, 203–207. [Google Scholar] [CrossRef] [PubMed]
  57. Klattenhoff, C.; Xi, H.; Li, C.; Lee, S.; Xu, J.; Khurana, J.S.; Zhang, F.; Schultz, N.; Koppetsch, B.S.; Nowosielska, A.; et al. The Drosophila Hp1 homolog rhino is required for transposon silencing and piRNA production by dual-strand clusters. Cell 2009, 138, 1137–1149. [Google Scholar] [CrossRef] [PubMed]
  58. Crooks, G.E.; Hon, G.; Chandonia, J.M.; Brenner, S.E. Weblogo: A sequence logo generator. Genome Res. 2004, 14, 1188–1190. [Google Scholar] [CrossRef] [PubMed]
  59. Chung, W.J.; Okamura, K.; Martin, R.; Lai, E.C. Endogenous RNA interference provides a somatic defense against Drosophila transposons. Curr. Biol. 2008, 18, 795–802. [Google Scholar] [CrossRef] [PubMed]
  60. Giraldez, A.J.; Mishima, Y.; Rihel, J.; Grocock, R.J.; Van Dongen, S.; Inoue, K.; Enright, A.J.; Schier, A.F. Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science 2006, 312, 75–79. [Google Scholar] [CrossRef] [PubMed]
  61. Tchurikov, N.A.; Kretova, O.V. Both piRNA and siRNA pathways are silencing transcripts of the suffix element in the Drosophila melanogaster germline and somatic cells. PLoS ONE 2011, 6, e21882. [Google Scholar] [CrossRef] [PubMed]
  62. Zong, Q.; Schummer, M.; Hood, L.; Morris, D.R. Messenger RNA translation state: The second dimension of high-throughput expression screening. Proc. Natl. Acad. Sci. USA 1999, 96, 10632–10636. [Google Scholar] [CrossRef] [PubMed]
  63. Krichevsky, A.M.; King, K.S.; Donahue, C.P.; Khrapko, K.; Kosik, K.S. A microRNA array reveals extensive regulation of microRNAs during brain development. RNA 2003, 9, 1274–1281. [Google Scholar] [CrossRef] [PubMed]
  64. Suh, N.; Baehner, L.; Moltzahn, F.; Melton, C.; Shenoy, A.; Chen, J.; Blelloch, R. MicroRNA function is globally suppressed in mouse oocytes and early embryos. Curr. Biol. 2010, 20, 271–277. [Google Scholar] [CrossRef] [PubMed]
  65. Dallaire, A.; Simard, M.J. The implication of microRNAs and endo-siRNAs in animal germline and early development. Dev. Biol. 2016, 416, 18–25. [Google Scholar] [CrossRef] [PubMed]
  66. Suh, N.; Blelloch, R. Small RNAs in early mammalian development: From gametes to gastrulation. Development 2011, 138, 1653–1661. [Google Scholar] [CrossRef] [PubMed]
  67. Fagegaltier, D.; Bouge, A.L.; Berry, B.; Poisot, E.; Sismeiro, O.; Coppee, J.Y.; Theodore, L.; Voinnet, O.; Antoniewski, C. The endogenous siRNA pathway is involved in heterochromatin formation in Drosophila. Proc. Natl. Acad. Sci. USA 2009, 106, 21258–21263. [Google Scholar] [CrossRef] [PubMed]
  68. Karaiskos, S.; Naqvi, A.S.; Swanson, K.E.; Grigoriev, A. Age-driven modulation of tRNA-derived fragments in Drosophila and their potential targets. Biol. Direct. 2015, 10, 51. [Google Scholar] [CrossRef] [PubMed]
  69. Sobala, A.; Hutvagner, G. Small RNAs derived from the 5’ end of tRNA can inhibit protein translation in human cells. RNA Biol. 2013, 10, 553–563. [Google Scholar] [CrossRef] [PubMed]
Genes 09 00082 g001 550
Figure 1. Alignment and categorization of small RNAs based on their origin. (A) Experimental design for total RNA extraction from un-fractionated (UF) and fractionated embryos has been published [35]. The profile is modified from Bor et al. 2006 [49] (B) Percentage of small RNAs in un-fractionated and fractionated embryos. The sequences were aligned in an order from down to top as shown in the figure legend. The aligned sequences included not only the perfect matches but also indel and up to 3 mismatches. mRNP: messenger ribonucleoprotein; ncRNA: non-coding RNA; miscRNA: miscellaneous RNA.
Figure 1. Alignment and categorization of small RNAs based on their origin. (A) Experimental design for total RNA extraction from un-fractionated (UF) and fractionated embryos has been published [35]. The profile is modified from Bor et al. 2006 [49] (B) Percentage of small RNAs in un-fractionated and fractionated embryos. The sequences were aligned in an order from down to top as shown in the figure legend. The aligned sequences included not only the perfect matches but also indel and up to 3 mismatches. mRNP: messenger ribonucleoprotein; ncRNA: non-coding RNA; miscRNA: miscellaneous RNA.
Genes 09 00082 g001
Genes 09 00082 g002 550
Figure 2. The Clustering of miRNA expression level. RNA-seq reads were first aligned to 430 mature miRNA sequences with perfect matches, which yielded a total of 263 mature miRNAs with at least 1 read in at least one sample. Setting the minimal threshold at 50 RPM filtered out 169 miRNAs, resulting in a total of 94 mature miRNAs. The log2 fold changes of these 263 miRNAs were clustered using Gene Cluster 3.0 [42] and visualized by Java Treeview [43] (A); Relative distribution of miRNA in fractions. In order calculate the relative abundance of miRNAs in each fraction, we calculated the percentage of each mature miRNA in the sum of all fractions of the developmental stage (B). All clustering was carried out using “Euclidian Distance” and “Complete Linkage” in Cluster 3.0.
Figure 2. The Clustering of miRNA expression level. RNA-seq reads were first aligned to 430 mature miRNA sequences with perfect matches, which yielded a total of 263 mature miRNAs with at least 1 read in at least one sample. Setting the minimal threshold at 50 RPM filtered out 169 miRNAs, resulting in a total of 94 mature miRNAs. The log2 fold changes of these 263 miRNAs were clustered using Gene Cluster 3.0 [42] and visualized by Java Treeview [43] (A); Relative distribution of miRNA in fractions. In order calculate the relative abundance of miRNAs in each fraction, we calculated the percentage of each mature miRNA in the sum of all fractions of the developmental stage (B). All clustering was carried out using “Euclidian Distance” and “Complete Linkage” in Cluster 3.0.
Genes 09 00082 g002
Genes 09 00082 g003a 550Genes 09 00082 g003b 550
Figure 3. Analysis of transposon-derived transcripts. (A) The 5′-5′ complementarity among the transposon matched small RNAs. The complementarity of reads derived from transposons were calculated by looking at 5’-to-5’ complementarity. The frequency at each length of complementation was plotted for all reads and reads normalized by the number of mapping sites in the genome. X-axis is the 5’-to-5’ distance and y-axis indicates read number. (-: minus) values indicate antisense reads; (B) The nucleotide bias at position 1 and 10 of the transposon-derived small RNAs. Sequence logos were generated using WebLogo [58]. The small RNAs that have a 10-nt 5’-to-5’ complementation from Panel A. The nucleotide distribution of transposon-derived small RNAs in 0–1 h (C) and 7–8 h embryos (D).
Figure 3. Analysis of transposon-derived transcripts. (A) The 5′-5′ complementarity among the transposon matched small RNAs. The complementarity of reads derived from transposons were calculated by looking at 5’-to-5’ complementarity. The frequency at each length of complementation was plotted for all reads and reads normalized by the number of mapping sites in the genome. X-axis is the 5’-to-5’ distance and y-axis indicates read number. (-: minus) values indicate antisense reads; (B) The nucleotide bias at position 1 and 10 of the transposon-derived small RNAs. Sequence logos were generated using WebLogo [58]. The small RNAs that have a 10-nt 5’-to-5’ complementation from Panel A. The nucleotide distribution of transposon-derived small RNAs in 0–1 h (C) and 7–8 h embryos (D).
Genes 09 00082 g003aGenes 09 00082 g003b
Genes 09 00082 g004 550
Figure 4. The Clustering of transposon-derived transcripts. (A) The reads were first aligned to the Repbase collection [40], all reads and uniquely mapped reads were used to calculate the transcript expression levels. Then the log2 fold change was calculated and clustered as in Figure 2A; (B) The reads that are mapped to one of the most abundant transposon cluster in Drosophila melanogaster, 42AB.
Figure 4. The Clustering of transposon-derived transcripts. (A) The reads were first aligned to the Repbase collection [40], all reads and uniquely mapped reads were used to calculate the transcript expression levels. Then the log2 fold change was calculated and clustered as in Figure 2A; (B) The reads that are mapped to one of the most abundant transposon cluster in Drosophila melanogaster, 42AB.
Genes 09 00082 g004
Genes 09 00082 g005a 550Genes 09 00082 g005b 550
Figure 5. Analysis of small RNA abundance in various developmental stages. (A) The percentage of transposon-derived small RNAs; (B) The small RNAs aligned to know RNAs as in Figure 1B for all datasets downloaded from GEO with the following GEO numbers; OSSC_#1:GSM385744, OSSC_#2:GSM385748, OSSC_#3:GSM385821, OSSC_#4:GSM385822, Unfertilzed_Eggs:GSM1550732, embryo_0–1h_#1:GSM180330, embryo_0–1h_#2:GSM286604, embryo_0–1h_#3:GSM286613, embryo_0–1h_#4:GSM868349, embryo_2–6h_#1:GSM180331, embryo_2–6h_#2:GSM286606, embryo_7–8h:GSM868350, embryo_6–10h:GSM180332, embryo_12–24h_#1:GSM180333, embryo_12–24h_#2:GSM364902, Larvae_1st_Ins:GSM360256, Larvae_2nd_Ins:GSM360257, Pupa_0–1day:GSM360260, Pupa_0–2day:GSM360262, Imaginal_Discs_#1:GSM180335, Imaginal_Discs_#2:GSM275691, Male_Head_#1:GSM180328, Male_Head_#2:GSM286601, Male_Body_#2:GSM180329, Male_Body_#2:GSM286602, S2 Cells:GSM180337, S2KC_Cells:GSM272651, S2GMR_Cells:GSM272652, KC_Cells:GSM272653. OSSC: Ovary Somatic Sheet (OSS) cell line. We numbered datasets (#) if they are the same time points. We only used 15–29 nt inserts.
Figure 5. Analysis of small RNA abundance in various developmental stages. (A) The percentage of transposon-derived small RNAs; (B) The small RNAs aligned to know RNAs as in Figure 1B for all datasets downloaded from GEO with the following GEO numbers; OSSC_#1:GSM385744, OSSC_#2:GSM385748, OSSC_#3:GSM385821, OSSC_#4:GSM385822, Unfertilzed_Eggs:GSM1550732, embryo_0–1h_#1:GSM180330, embryo_0–1h_#2:GSM286604, embryo_0–1h_#3:GSM286613, embryo_0–1h_#4:GSM868349, embryo_2–6h_#1:GSM180331, embryo_2–6h_#2:GSM286606, embryo_7–8h:GSM868350, embryo_6–10h:GSM180332, embryo_12–24h_#1:GSM180333, embryo_12–24h_#2:GSM364902, Larvae_1st_Ins:GSM360256, Larvae_2nd_Ins:GSM360257, Pupa_0–1day:GSM360260, Pupa_0–2day:GSM360262, Imaginal_Discs_#1:GSM180335, Imaginal_Discs_#2:GSM275691, Male_Head_#1:GSM180328, Male_Head_#2:GSM286601, Male_Body_#2:GSM180329, Male_Body_#2:GSM286602, S2 Cells:GSM180337, S2KC_Cells:GSM272651, S2GMR_Cells:GSM272652, KC_Cells:GSM272653. OSSC: Ovary Somatic Sheet (OSS) cell line. We numbered datasets (#) if they are the same time points. We only used 15–29 nt inserts.
Genes 09 00082 g005aGenes 09 00082 g005b
Table
Table 1. The percentage of the most abundant 10 micro RNAs. The sequences that perfectly matched to hairpin and mature miRNAs were used to calculate read frequency. The most abundant 10 miRNAs were selected from each sample and then pooled to generate a panel of 18 miRNAs that appear to be the most abundant miRNAs in 0–1 h and 7–8 h Drosophila melanogaster embryos. mRNP, messenger ribonucleoprotein.
Table 1. The percentage of the most abundant 10 micro RNAs. The sequences that perfectly matched to hairpin and mature miRNAs were used to calculate read frequency. The most abundant 10 miRNAs were selected from each sample and then pooled to generate a panel of 18 miRNAs that appear to be the most abundant miRNAs in 0–1 h and 7–8 h Drosophila melanogaster embryos. mRNP, messenger ribonucleoprotein.
HairpinMature_3bp
mRNP60S80SPolysomeTotalmRNP60S80SPolysomeTotal
0–1 h7–8 h0–1 h7–8 h0–1 h7–8 h0–1 h7–8 h0–1 h7–8 h0–1 h7–8 h0–1 h7–8 h0–1 h7–8 h0–1 h7–8 h0–1 h7–8 h
dme-bantam-3p0.33.436.327.335.926.42.91.01.530.30.33.335.726.235.726.12.60.91.530.1
dme-miR-1-3p7.75.53.17.914.819.91.07.543.913.77.75.53.17.914.819.91.07.443.913.7
dme-miR-184-3p17.59.67.64.610.46.127.025.911.13.617.59.67.64.610.46.127.025.911.13.6
dme-miR-263a-5p18.024.71.54.23.45.133.122.23.59.018.024.71.54.23.45.133.122.13.59.0
dme-miR-283-5p0.20.30.91.31.31.91.41.80.21.40.20.30.91.21.31.81.41.80.11.4
dme-miR-286-3p2.010.87.118.17.518.96.810.32.211.62.010.87.118.17.518.96.810.32.211.5
dme-miR-305-5p1.00.01.50.30.40.10.30.01.70.11.00.01.30.20.30.00.30.01.70.1
dme-miR-306-5p2.70.00.10.00.10.00.10.02.40.12.70.00.10.00.10.00.10.02.40.1
dme-miR-311-3p2.10.00.60.11.00.10.20.12.60.12.10.00.60.11.00.10.20.12.60.1
dme-miR-315-5p0.20.80.60.80.30.51.40.80.22.00.20.80.60.80.30.51.40.80.22.0
dme-miR-5-5p16.427.23.44.32.43.55.97.03.34.516.427.23.13.92.33.35.96.93.24.2
dme-miR-7-5p0.11.50.51.20.61.61.11.40.31.00.11.50.51.20.61.61.11.40.31.0
dme-miR-8-3p0.91.23.02.93.33.01.21.71.32.50.51.12.72.53.12.71.21.61.22.2
dme-miR-92b-3p0.30.01.61.02.01.20.20.30.80.10.30.01.60.92.01.20.20.30.80.1
dme-miR-958-3p0.00.14.25.10.50.60.10.20.00.80.00.14.15.00.50.60.00.10.00.8
dme-miR-9a-5p6.05.31.63.51.72.91.78.64.16.46.05.31.63.41.72.91.78.64.16.4
dme-miR-9c-5p11.54.18.82.46.01.89.14.49.42.611.54.18.82.46.01.89.04.49.42.5
dme-miR-iab-4-5p1.20.10.10.10.00.00.10.00.50.11.60.10.00.10.00.00.20.10.70.1
Table
Table 2. Categorization of micro RNAs based on their association with translational machinery. Based on fold of induction between 0–1 h and 7–8 h embryos, miRNAs were divided into 4 major miRNA groups (Group 1-overexpressed in 7–8 h unfractionated embryos; Group 2-equally expressed in unfractionated embryos; Group 3-downregulated in 7–8 h unfractionated embryos; Group 4; Others). The fold changes in red, green and light blue indicate increase, decrease and no significant change in 7–8 h embryos, respectively.
Table 2. Categorization of micro RNAs based on their association with translational machinery. Based on fold of induction between 0–1 h and 7–8 h embryos, miRNAs were divided into 4 major miRNA groups (Group 1-overexpressed in 7–8 h unfractionated embryos; Group 2-equally expressed in unfractionated embryos; Group 3-downregulated in 7–8 h unfractionated embryos; Group 4; Others). The fold changes in red, green and light blue indicate increase, decrease and no significant change in 7–8 h embryos, respectively.
GmiRNATotmRNP60S80SPolyGmiRNATotmRNP60S80SPoly
1dme-miR-286-3p2.63.12.22.11.62dme-miR-11-3p0.81.30.80.52.4
1dme-miR-283-5p3.61.51.31.31.32dme-miR-9b-5p0.40.8−0.4−0.21.7
1dme-miR-7-5p2.04.52.02.11.32dme-miR-1012-3p−0.8−0.8−0.5−1.01.4
1dme-miR-957-3p2.61.31.21.20.82dme-miR-1010-3p0.2−0.60.40.51.4
1dme-miR-263a-5p1.61.12.31.30.42dme-miR-5-5p0.61.41.21.31.2
1dme-miR-315-5p3.92.61.31.30.22dme-miR-12-5p0.30.40.61.21.1
1dme-miR-314-3p5.62.11.31.00.02dme-miR-13a-3p0.20.61.30.20.0
1dme-bantam-5p2.61.71.51.2−0.82dme-miR-31a-3p0.20.01.00.7−0.3
1dme-miR-958-3p7.03.61.10.81.82dme-miR-124-3p0.0−0.31.20.2−0.3
1dme-miR-956-3p3.31.21.00.60.72dme-miR-965-5p0.9−0.40.70.80.8
1dme-miR-998-3p1.31.21.50.50.42dme-miR-1003-3p−0.8−0.50.10.00.0
1dme-miR-2c-5p1.41.21.60.20.22dme-miR-190-5p0.40.10.00.1−0.7
1dme-miR-304-5p6.01.50.90.50.12dme-miR-282-3p−0.1−1.9−1.7−0.7−0.7
1dme-miR-8-3p1.11.80.80.51.32dme-miR-2a-3p0.60.20.80.5−1.5
1dme-miR-983-5p2.01.20.70.30.02dme-miR-33-5p−0.3−2.4−1.6−1.1−1.9
1dme-bantam-3p4.64.30.40.3−0.52dme-miR-281-2-5p−0.5−1.30.3−0.10.6
1dme-miR-2b-3p1.61.20.40.8−1.12dme-miR-2b-2-5p−0.8−1.60.51.00.5
1dme-miR-276a-3p3.30.91.41.00.82dme-miR-375-3p0.1−2.30.8−0.20.1
1dme-miR-5-3p2.70.81.71.02.23dme-miR-3-5p−2.80.0−0.1−0.10.0
1dme-miR-263b-5p1.50.71.91.21.53dme-miR-184-3p−1.4−0.20.10.00.9
1dme-miR-137-3p1.3−0.61.61.31.43dme-miR-312-5p−3.8−0.2−1.1−0.6−0.6
1dme-miR-1002-5p3.1−0.21.21.91.43dme-miR-6-1-5p−4.2−0.3−0.7−0.80.0
1dme-miR-10-5p1.90.21.20.72.23dme-miR-3-3p−1.8−0.40.40.50.7
1dme-miR-252-5p1.50.81.70.70.53dme-miR-999-3p−1.7−0.4−0.2−0.10.6
1dme-miR-316-5p3.40.01.30.9−0.13dme-miR-305-3p−1.0−0.7−0.2−0.60.0
1dme-miR-1006-3p2.50.51.10.50.03dme-miR-9c-5p−1.7−0.8−1.0−1.00.0
1dme-miR-968-5p4.20.61.10.30.03dme-miR-13b-2-5p−1.4−0.80.70.60.1
1dme-miR-1000-5p1.90.11.10.60.73dme-miR-9b-3p−2.5−1.10.0−0.1−0.3
1dme-miR-987-5p2.10.70.20.11.13dme-miR-79-3p−2.3−1.1−0.60.00.7
1dme-miR-31a-5p1.30.90.50.10.83dme-miR-92a-3p−2.3−1.30.3−0.31.4
1dme-miR-13b-3p1.10.10.90.70.53dme-miR-124-5p−1.5−1.30.50.50.4
1dme-miR-927-3p1.80.60.60.60.53dme-miR-1002-3p−1.6−1.60.30.20.0
1dme-miR-996-3p1.0−0.2−0.20.20.13dme-miR-308-3p−1.4−1.7−0.8−0.1−0.5
1dme-miR-284-5p2.40.40.60.10.03dme-miR-310-3p−3.0−1.7−1.4−1.7−0.3
1dme-miR-927-5p2.1−0.60.50.40.03dme-miR-92b-3p−2.5−1.90.10.01.2
1dme-miR-87-3p2.6−0.30.70.0−0.13dme-miR-968-3p−2.1−2.20.30.50.7
1dme-miR-281-3p1.6−0.2−0.2−0.7−0.13dme-miR-995-3p−1.4−2.4−0.1−0.40.4
1dme-miR-1008-3p1.4−0.20.50.3−0.53dme-miR-312-3p−5.2−2.5−1.2−0.9−0.2
1dme-miR-8-5p1.1−1.61.01.02.63dme-miR-279-3p−3.9−2.60.30.00.8
1dme-miR-10-3p2.2−1.00.61.22.03dme-miR-6-2-5p−4.9−2.7−1.0−0.10.4
1dme-miR-14-3p1.0−1.41.5−0.31.13dme-miR-iab-8-3p−2.3−3.20.80.5−0.3
4dme-miR-1-3p−1.50.22.21.23.83dme-miR-iab-4-5p−2.3−3.31.00.7−0.3
4dme-miR-988-3p−1.0−0.11.41.20.03dme-miR-275-3p−2.9−3.6−1.2−1.40.2
4dme-miR-4-3p−1.1−0.52.20.71.43dme-miR-282-5p−2.5−3.7−1.7−1.2−1.4
4dme-miR-31b-5p−1.2−0.90.1−0.91.73dme-miR-311-3p−5.2−4.8−1.2−2.0−0.2
4dme-miR-981-3p1.0−1.9−0.10.60.43dme-miR-305-5p−4.0−4.9−1.8−1.8−1.3
2dme-miR-9a-5p0.80.51.91.53.33dme-miR-306-5p−5.1−5.2−0.6−0.8−1.1
Table
Table 3. Quantitative PCR validation of miRNA expression levels in unfractionated embryos. Quantitative PCR analyses were carried out with total RNAs from 0–1 h and 7–8 h embryos in which U6 ncRNA was used for normalization. Relative miRNA expression was then calculated in Log2 scale (7–8 h embryos/0–1 h embryos).
Table 3. Quantitative PCR validation of miRNA expression levels in unfractionated embryos. Quantitative PCR analyses were carried out with total RNAs from 0–1 h and 7–8 h embryos in which U6 ncRNA was used for normalization. Relative miRNA expression was then calculated in Log2 scale (7–8 h embryos/0–1 h embryos).
miRNAqPCRRNA-Seq
dme-miR-1-3p−1.1−1.5
dme-miR-3-3p−1.5−1.8
dme-miR-5-5p0.240.6
dme-miR-7-5p1.72.0
dme-miR-8-3p1.31.1
dme-miR-1002-5p3.93.1

Share and Cite

MDPI and ACS Style

Cosacak, M.I.; Yiğit, H.; Kizil, C.; Akgül, B. Re-Arrangements in the Cytoplasmic Distribution of Small RNAs Following the Maternal-to-Zygotic Transition in Drosophila Embryos. Genes 2018, 9, 82. https://doi.org/10.3390/genes9020082

AMA Style

Cosacak MI, Yiğit H, Kizil C, Akgül B. Re-Arrangements in the Cytoplasmic Distribution of Small RNAs Following the Maternal-to-Zygotic Transition in Drosophila Embryos. Genes. 2018; 9(2):82. https://doi.org/10.3390/genes9020082

Chicago/Turabian Style

Cosacak, Mehmet Ilyas, Hatice Yiğit, Caghan Kizil, and Bünyamin Akgül. 2018. "Re-Arrangements in the Cytoplasmic Distribution of Small RNAs Following the Maternal-to-Zygotic Transition in Drosophila Embryos" Genes 9, no. 2: 82. https://doi.org/10.3390/genes9020082

APA Style

Cosacak, M. I., Yiğit, H., Kizil, C., & Akgül, B. (2018). Re-Arrangements in the Cytoplasmic Distribution of Small RNAs Following the Maternal-to-Zygotic Transition in Drosophila Embryos. Genes, 9(2), 82. https://doi.org/10.3390/genes9020082

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop