- Research article
- Open Access
The germline of the malaria mosquito produces abundant miRNAs, endo-siRNAs, piRNAs and 29-nt small RNAs
© Castellano et al.; licensee BioMed Central. 2015
Received: 4 August 2014
Accepted: 19 January 2015
Published: 19 February 2015
Small RNAs include different classes essential for endogenous gene regulation and cellular defence against genomic parasites. However, a comprehensive analysis of the small RNA pathways in the germline of the mosquito Anopheles gambiae has never been performed despite their potential relevance to reproductive capacity in this malaria vector.
We performed small RNA deep sequencing during larval and adult gonadogenesis and find that they predominantly express four classes of regulatory small RNAs. We identified 45 novel miRNA precursors some of which were sex-biased and gonad-enriched , nearly doubling the number of previously known miRNA loci. We also determine multiple genomic clusters of 24-30 nt Piwi-interacting RNAs (piRNAs) that map to transposable elements (TEs) and 3’UTR of protein coding genes. Unusually, many TEs and the 3’UTR of some endogenous genes produce an abundant peak of 29-nt small RNAs with piRNA-like characteristics. Moreover, both sense and antisense piRNAs from TEs in both Anopheles gambiae and Drosophila melanogaster reveal novel features of piRNA sequence bias. We also discovered endogenous small interfering RNAs (endo-siRNAs) that map to overlapping transcripts and TEs.
This is the first description of the germline miRNome in a mosquito species and should prove a valuable resource for understanding gene regulation that underlies gametogenesis and reproductive capacity. We also provide the first evidence of a piRNA pathway that is active against transposons in the germline and our findings suggest novel piRNA sequence bias. The contribution of small RNA pathways to germline TE regulation and genome defence in general is an important finding for approaches aimed at manipulating mosquito populations through the use of selfish genetic elements.
In recent years additional layers of complexity have been revealed in the regulation of gene expression following the discovery in animals of several classes of small RNA molecules that can act at both the transcriptional and post-transcriptional level. Many of these small RNAs themselves show tissue-specific expression and have been shown to be essential for correct organogenesis and developmental progression. Anopheles gambiae mosquitoes are the major vectors of Plasmodium malaria parasites. Successful malaria control initiatives in the past have all relied on reducing the reproductive capacity of mosquito populations. Therefore a better understanding of the processes that regulate sexual development and, in particular, gonadogenesis and gametogenesis could provide novel targets for vector control. Small regulatory RNAs such as microRNAs (miRNAs) and other classes of small non-coding RNAs play a role in the germline of many organisms in germline stem cell maintenance and in restricting the expression of transposable elements yet little is known about the diversity of small RNAs and their contribution in the malaria mosquito [1-3].
Each of the different small RNA classes are characterised by their ability to interact with Argonaute (AGO) proteins, all of which are involved in gene silencing mechanisms . Studies, principally performed with the fruitfly Drosophila melanogaster, indicate that in the fly AGO proteins can be divided into two different clades. One clade contains AGO1 and AGO2 that are expressed ubiquitously and function in gene silencing through binding with microRNAs (miRNAs) and small interfering RNAs (siRNAs), respectively . The second clade contains the Piwi proteins, specifically expressed in the germline, composed of AGO3, Piwi and Aubergine (Aub) that bind to Piwi interacting RNAs (piRNAs).
miRNAs are a large class of ~21-24-nt small RNAs produced by DICER1 processing of endogenously expressed RNA hairpin structures [6,7] that are involved in post-transcriptional gene repression [5,8]. After DICER1 processing, the derived mature duplexes are unwound, loaded onto AGO1 to guide it on the gene targets that are recognized through incomplete base pairing with the loaded single stranded, miRNA, enabling AGO1 to repress protein translation and/or destabilise the mRNA transcript [5,8]. On the other hand, siRNAs are exactly 21-nt long RNA molecules and produced by sequential DICER2 cleavage of long double strand RNAs (dsRNAs) . They are then unwound and loaded onto AGO2 as a single stranded guide siRNA . The complete complementarity between the loaded siRNA and the target permits AGO2 to mediate the cleavage of the target, that occurs opposite to the 10th and the 11th nt of the annealed siRNA . piRNAs are, instead, 24-29–nt long and are particularly expressed in the germline, from discrete genomic loci and have been shown in several organisms to be involved in the silencing of genomic repeats and active transposable elements (TEs) [2,11,12]. Their biogenesis does not depend on DICER enzymes, but they are generated by a primary biogenesis pathway that is not completely understood and by a secondary biogenesis pathway called the ping-pong mechanism: piRNAs derived from one genomic strand generate the 5’ end of new piRNAs from the opposite strand due to the endonucleolytic activity of the Piwi proteins [2,13]. So far miRNAs have been exclusively involved in silencing of endogenous gene targets while siRNAs are involved in the repression of host genes, TEs and viruses. piRNAs are predominantly involved in the silencing of TEs, although a few examples of control of non-TE elements may exist (reviewed in ). Indeed in the mosquito Aedes aegypti, in addition to canonical RNAi-mediated silencing of viruses [15-17], piRNA-like molecules have also been implicated in this process.
A few studies have attempted to identify computationally the microRNA pool in Anopheline mosquitoes yet only a small fraction of these microRNAs have been confirmed by sequencing or other approaches and the microRNA complement for A. gambiae is much smaller than that described for Drosophila, suggesting that it is far from completely described [20-25]. Moreover, two Dicer enzymes, DCR1 and DCR2, and 5 AGO proteins AGO1-5 have been identified in A. gambiae suggesting that all the genetic machinery is there to process the full range of small regulatory RNAs . Due to this lack of knowledge we set out to clone and analyse the small RNA populations present during the formation of the gonads in each sex to evaluate the contribution of each pathway to this essential process.
We identify a large number of novel miRNAs, some of which are sex-biased and/or enriched during gametogenesis, some of which represent new miRNA gene clusters, and others an expansion of existing clusters. We also identify endo-siRNAs derived from overlapping, convergently transcribed protein coding genes (cis-NAT-siRNAs)  and TEs and 5’-half-tRNAs that are 32-nt small RNAs formed by the processing of tRNA hairpin [28,29] whose expression was significantly downregulated in pre-vitellogenic ovaries. We finally identified in the gonads a large class of piRNAs, predominantly derived from TEs, for which we described novel sequence bias that may be relevant for piRNA recognition, loading and/or biogenesis.
This is the first report describing the complexity of the small RNA transcriptome during the essential processes of gonadogenesis and gametogonesis in the mosquito and the predominant role of piRNAs in limiting transposon proliferation in the germline is relevant to vector control approaches that propose these elements as agents to modify mosquito populations.
Results and discussion
Characteristics of the Anopheles gambiae small RNA population
Novel miRNAs discovered in the mosquito
In the analysed tissues we confirmed the expression of all the 67 precursors producing miRNAs from Anopheles gambiae annotated on the last release of miRBase (Release 20) database (Additional file 2). Apart from aga-miR-929 and aga-miR-219 that presented 2 and 5 total reads respectively, all the other miRNAs ranged from 28 to 974,600 reads and were expressed in multiple tissues (Additional files 3 and 4). The most expressed miRNA corresponded to aga-miR-263a (Additional file 4). In addition the top 20 most expressed miRNAs represent 96.3% of all the miRNA reads. Mapping the reads derived from the analysed samples we could also identify 62 previously unidentified “star” miRNA molecules that map on the opposite arm and with 2-nt 3’ overhang of the miRNA precursor respect to the known miRNA (Additional files 2 and 3).
In mammals and plants a significant fraction of miRNAs derive from repetitive sequences [38-42]. To evaluate the extent of miRNAs that originate from repetitive sequences in the mosquito we used RepeatMasker (http://www.repeatmasker.org) on all the miRNA precursors, including both known and novel genes identified in this study. Surprisingly, miRNAs in the mosquito rarely derive from repetitive sequences, with only one novel miRNA derived from a low complexity (CGT)n repeat.
Overall the novel miRNAs identified were prevalently derived from intergenic regions (35%) or the sense strand of introns (62%) (Additional file 5). Just 2 novel miRNAs seem to derive from exonic regions.
A recent study using a similar sequencing approach applied to A. gambiae mosquitoes pre- and post-bloodmeal revealed 58 novel miRNAs, showing the utility of a targeted deep sequencing approach . 18 of these miRNAs were also revealed by our analysis, ultimately reducing the number of truly novel miRNAs in our study to 45 (Additional file 5).
Several miRNA genes are upregulated in the gonads and specifically during gametogenesis
We next performed a differential expression analysis of both the novel and known miRNAs between the different analysed tissues. The read counts for each miRNA were compared between conditions using the Bioconductor DESeq package, which uses a negative binomial distribution model to test for differential expression in deep sequencing datasets (Additional file 7) [44,45]. Quantification of relative expression by qPCR was used to confirm tissue-specific expression for a subset of these microRNAs that showed varying but significant differential expression according to DESeq (Additional file 8b). Hierarchical clustering according to expression profile identified 3 main classes according to their breadth and intensity of expression (Figure 2C). Unsurprisingly Class I, comprising miRNAs abundantly expressed across all samples, was predominantly made up of previously annotated miRNA, whereas classes II and III, comprising lowly-expressed miRNAs and tissue-restricted miRNAs, respectively, were enriched for novel miRNAs, reflecting the greater sensitivity afforded by dissecting gonadal tissues at various stages of development. Among the tissue restricted miRNAs we found 41 miRNAs that were differentially expressed (adjusted P- value < 0.01) between the testes and pre-vitellogenic ovaries, including two (aga-mir-2944a-2 and aga-mir-286b) that were also upregulated both in the testes and during oogenesis, suggesting a general role in gametogenesis (Additional files 7 and 8a). Validating our approach, microRNAs with known roles in oogenesis such as miR-989 were heavily ovary-enriched in our analysis [20,46]. Larval stages of the mosquito show virtually no sexual dimorphism in external features and consistent with this we failed to find microRNAs that were significantly differentially expressed between the two sexes. However mIR-989 and aga-mir-10361b, that were ovary and testes-enriched in the adult respectively, were similarly enriched in the immature larval gonad of each sex (Additional file 7: Table S4), indicating that each microRNA must play an early role in the formation of the gonad. Interestingly, testes-biased miRNAs showed a non-random chromosomal distribution – of the 31 miRNAs that showed strong testes-bias 20 (64.5%) mapped exclusively to the X chromosome, compared to only 23% for non-testes biased miRNAs (Additional file 7: Table S4; Fisher’s exact test, p < 2E-05). Protein coding genes on the X chromosome are usually silenced due to meiotic sex chromosome inactivation (MSCI) during spermatogenesis , however our findings are consistent with recent evidence that X-linked microRNAs can escape this inactivation . Moreover the high number of testes-biased novel microRNAs we identified on the X chromosome suggests that for microRNA genes at least this chromosome represents a favourable environment for male-biased microRNAs to evolve, as others have suggested .
Endo-siRNAs are preferentially expressed in adult testes of Anopheles gambiae
AGAP003387 locus produced high abundance of one 29-nt piRNA from its 3’UTR
The 29-nt fraction derived from mRNAs actually derived almost exclusively from the sense strand of the 3’UTR of the gene AGAP003387 and was dominated by a unique sequence that is the most abundant read in all the samples (about 1-4% of all the reads) (Figure 4D). All the sample tissues presented this identical pattern. The nuclear size distribution (26-29 nt) of each unique small RNA sequence (Figure 4D) that derives from this locus and the presence of a bias of uridine (U) as the first nucleotide in these sequences, suggests that they have some piRNA-like characteristics (Figure 4E).
piRNAs map to discrete genomic loci and are produced abundantly from transposable elements
In order to define piRNA expression clusters in the Anopheles genome we selected only reads containing a U1 or A10 that should represent a piRNA-enriched fraction from all samples. We than mapped these reads on the genome with perfect match at a unique position. This strategy has been successfully used to find the genomic origin of piRNAs in other organisms . We considered loci as producing piRNAs if the average length of reads was 26-27-nt, abundance was more than 50 mean count of reads and loci larger than 100-bp. Using these criteria we revealed that Anopheles piRNAs map to more than 1500 discrete genomic loci (Additional file 13).
A subgroup of TEs produce predominantly 29-nt sized small RNAs
Anopheles piRNAs show novel sequence bias in addition to the classical characteristics of A10 and U1 bias
Further analysis of the sequence characteristics of putative Anopheles piRNAs that derive from both the sense and the antisense strand of TEs revealed other nucleotide biases and further complexity at positions of the piRNA in addition to the described A10 bias for sense piRNAs and U1 bias for antisense piRNAs. In particular piRNAs derived from the sense strand, in addition to showing a predominance of A10, are also highly enriched for U1. To exclude the possibility that the sense piRNA population was a mixture of primary (U1) and secondary (A10) piRNAs, we analyzed just the reads derived from the TEs that contain A10 and we observed that they also preferentially contained a U at position 1, confirming that these sequence biases frequently manifest in the same piRNA molecule (data not shown). Interestingly, the ratio A10 vs U1 decreased with increasing piRNA size in the various tissues, suggesting qualitative differences between the species of piRNA that may be related to their mode of biogenesis (Figure 5B and D). To our knowledge this is the first report of such sequence bias and was revealed by the additional sensitivity in detection afforded by separating piRNAs into different size classes. To confirm if this phenomenon held true across different organisms we re-analysed the dataset of Drosophila piRNAs that were specifically bound by the Piwi proteins PIWI, AUB and AGO3 , split these by size class and used pLogo to check for sequence bias (Additional file 16). As expected all populations displayed the known signatures yet strikingly the population of piRNAs derived from AGO3, that are heavily enriched for sense strand piRNAs, additionally contained a sequence bias of U (40%) at the 5’ terminal position that was statistically significant and similar to what we observed in the mosquito (Additional file 16C). This U1 enrichment also for sense piRNAs suggests that AGO3 preferentially loads U1 enriched small RNAs as for other argonaute family members in worm, fly and fly and human [58-61]. We also noticed a significant enrichment of cytosines in the last 2 nt of the piRNAs of each size class that was present in each of the two species regardless of piRNA size (Figure 5 and Additional file 16).
Anopheles gambiae produces an abundant class of 32-nt half tRNA that is down-regulated in pre-vitellogenic ovaries
A distinct population of longer small RNAs, 32 nucleotides in length, were among the most abundant class of RNA revealed in all samples that mapped exclusively to tRNA genes, and specifically the 5’ end. These 5’ half tRNAs have been previously described in other organisms and have recently been shown to have wide ranging yet fully elucidated roles in a wide range of host processes ranging from translational inhibition, stress response and signalling [29,62,63]. Though we cannot speculate a role for these half tRNAs in the mosquito it is interesting to note that they were abundantly expressed in all samples but showed a 30-fold upregulation in ovaries specifically during vitellogenesis (Additional file 7).
In this study we have greatly increased the known complement of small regulatory RNAs in the mosquito A. gambiae, identifying a large number of novel microRNAs, revealing the extent of endogenous siRNA production and describing for the first time the presence of piRNAs with likely roles in transposon control in the germline and possibly also roles in the control of a limited number of endogenous genes.
By focusing on the gonads at various stages of development we were able to sample the microRNAs expressed during the development and maintenance of these tissues at a much higher sensitivity than would have been afforded by examining the whole animal. This is the case even more so for the testes, given the relatively small size of this tissue and may in part explain the higher proportion of novel microRNAs recovered from the testis-enriched pool. The microRNA profile of both the ovary and testes through development and at different stages of gametogenesis described here should provide a valuable resource for better understanding the regulation of gene expression during this crucial process. In the future further functional analysis through mutagenesis or mis-expression of those microRNAs specifically enriched during the process of male or female gametogenesis should reveal candidate genes whose disruption could block mosquito reproduction.
The piRNA pathway has been shown in several organisms to be essential in the germline for repressing the expression of transposable elements, whose unregulated transposition would cause deleterious effects and loss of reproductive fitness. This report is the first description of the piRNA landscape in the reproductive tissues of the mosquito. We detected abundant piRNAs expressed from numerous clusters around the genome with sequence homology to all transposable element families identified in the A. gambiae genome. Different TE species produce distinct piRNA pools that were indicative of primary piRNA silencing only or a combination of primary and ping-pong produced piRNAs. The three Piwi proteins PIWI, AGO3 and Aub have different roles in the two mechanisms of piRNA generation and mutations in any of them lead to a de-repression of TE activity in the model insect Drosophila, indicating that they have non-redundant roles [64,65]. Interestingly our results show that different transposable elements can produce piRNA pools with markedly different characteristics, consistent with a route of biogenesis by specific members of the Piwi family. In Drosophila, members of the gypsy family of retrotransposons are expressed as virion-like particles in the somatic follicle cells of the gonad that surround the oocyte. In the follicle cells only PIWI is expressed of the 3 Piwi proteins, leading to silencing by primary piRNAs generated from transcripts of flamenco, a locus rich in degenerate copies of gypsy . Similarly we saw that mosquito piRNAs mapping to the various gypsy elements were almost exclusively anti-sense, consistent with a primary biogenesis by PIWI, though as yet we do not know if there is a functional equivalent of flamenco as a piRNA master locus in the mosquito, or if in the somatic follicle only PIWI is expressed. Orthologues of each of the 3 Piwi genes are present in the A. gambiae genome and expression analysis in the closely related mosquito A. stephensi reveals that each is ovary-enriched and upregulated during the process of oogenesis  As the genetic tools available in the mosquito improve it should be possible to dissect the role of various AGO and Piwi proteins in small RNA biogenesis and activity, and to confirm Piwi-association of some of the putative piRNA classes that we have revealed.
In terms of the characteristics of the piRNA populations we discovered, many features are conserved with other organisms such as size (24-30 nt), the sequence bias at positions 1 and 10 that are signatures of a ping pong mechanism of biogenesis. On the other hand, we reveal some novelties in the generation of piRNAs such as a population of piRNAs recognizing the 3’ UTR of endogenous genes, previously undocumented sequence biases in piRNAs corresponding to the sense strand of TEs and an abundant 29 nt class of piRNA.
Despite the wide interest in piRNAs in recent years, many questions remain unresolved. For example, although the mechanism by which transposable elements are repressed in most cases is likely a combination of mRNA degradation through the Slicer activity of piRNA-loaded Piwi members and piRNA-guided heterochromatic silencing of expression, there are also reports of piRNA-directed translational repression [67,68]. Moreover, how the initial trigger RNA from the piRNA cluster is produced, recognized and loaded onto a Piwi protein is still not clearly resolved. The novel characteristics of the piRNA pool in the mosquito adds clues to their origin may help in resolving the mechanism of piRNA biogenesis in general. Furthermore, there are several instances of piRNAs and other small RNAs derived from repetitive elements on the sex chromosomes that have been co-opted to mediate interactions between the chromosomes and can play a role in sex determination [69,70]. It remains to be seen whether such a role exists in the mosquito, where much of the Y chromosome is made of repetitive sequences including several transposon relics .
Certainly the extensive characterization of the piRNA pool provided here and the demonstration of extensive TE-derived piRNA pools that are abundantly expressed in the germline provides an answer as to how the mosquito manages to control the proliferation of TEs in its genome. Attempts to introduce into a mosquito population anti-pathogen or otherwise beneficial constructs through the use of TEs designed to spread at super-Mendelian frequency, as has been proposed , will have to make contingency for this genome defence mechanism.
All animal work was conducted according to UK Home Office Regulations and approved under Home Office License PPL 70/6453.
Gonads were dissected in PBS from 3–4 day old adult male and female mosquitoes (G3 strain) that had been reared under standard insectary conditions (adults reared at 27°C, 80% humidity, fed ad libitum on 5% glucose). To obtain vitellogenic ovaries females were dissected 24 hours after a bloodmeal from an anaesthetised mouse. In order to sex larval stages that are otherwise morphologically indistinguishable we used a mosquito line containing a sex-linked insertion of the visible marker gene RFP (Nolan, unpublished) that could allow us to unambiguously separate male and female larvae in the progeny of a cross between RFP-positive males and wild type females. Each sample contained tissue from a minimum of 10 individuals. The small size and delicate nature of the larval gonads prevented their removal intact from the larval body. In order to enrich for the gonadal tissue we dissected segments 5 to 7 from L4 larvae using fine 30 gauge needles. We separately confirmed under microscopy that this section consistently contained both male and female developing gonads. RNA was then isolated from this tissues using TRIzol Reagent (Life Technologies), following the manufactured instructions. RNA quality and integrity was assessed using a Bioanalyzer instrument (Agilent Technologies Genomics).
Small RNA libraries preparation and sequencing
One μg of purified RNA from two biological replicate per condition was used to prepare small RNA libraries according to the TruSeq Small RNA Sample preparation kit (Illumina) instructions. Fifty base pair (bp), single end sequencing was performed using the HiSeq 2000 instrument (Illumina).
Reads preprocessing and mapping
We clipped out the 3’ adaptors form the reads obtained at the end of the sequencing run using the FASX-toolkit from Hannon lab (http://hannonlab.cshl.edu/fastx_toolkit/) before further analyses. The processed reads from each sample were then mapped to the AgamP3 genome assembly using Bowtie version 0.12.7 allowing for 0 mismatches. Annotation of small RNAs was usually done mapping the reads on data downloaded from VectorBase (https://www.vectorbase.org/downloads) with the exception of transposable elements that were instead obtained from the University of California at Santa Cruz (USCS) Genome Browser (http://genome.ucsc.edu).
miRNA discovery and miRNA expression profile
To discover novel miRNAs we used both miRCat from the UEA small RNA Workbench (https://srna-workbench.cmp.uea.ac.uk/) and miRDeep2 algorithms [31,32] maintaining default setting and filtering reads by size ≥17. Among the novel candidates discovered using this approaches, high-confidence miRNAs containing both mature and star sequences complementary with 2-nt 3′ overhang detected in multiple samples were considered. To discover mirtrons we mapped all the reads coming from the various samples on introns smaller than 100-nt and manually inspected them. Newly discovered miRNAs were then quantified among the various tissue samples using the quantifier module from the miRDeep2 package. Differential expression profile of both novel and known miRNAs across the samples was performed using DESeq from the Bioconductor project (http://www.bioconductor.org).
Normalised expression values from the mIRDeep analysis were used to create a heatmap with the software GENE-E (www.broadinstitute.org) using log10 values. Only mIRs with a read count >20 in at least one condition were included in the heatmap and were hierarchically clustered by Euclidean distance.
Evaluation of miRNAs located on genomic repeats
To evaluate whether annotated miRNAs are located on genomic repeats and to discover the nature of these repeats, novel and downloaded miRNA precursors from the miRBase release 20 were analyzed using the RepeatMasker script, version 3.2.8 (http://www.repeatmasker.org/).
qPCR validation of miRNA expression
Reverse transcription was primed using a stem-loop primer with 8 nucleotides of complementarity to the target mIR of interest and cDNA synthesized using Superscript III reverse transcriptase, followed by PCR amplification using a mIR-specific primer and a primer specific to the stem loop . Reactions were performed in a StepOne Plus RT-PCR machine (Invitrogen) and PCR product was quantified by measuring SYBR green incorporation. The broadly expressed and non-sex-specific microRNA bantam was used as an internal control and the comparative Ct method was used to compare miRNA amounts between samples. 2 biological replicates were performed for each sample (TE, BF or OV) and a mimimum of 10 individuals were included per sample.
piRNA cluster analysis
To discover piRNA clusters, few clusters of small RNA reads that map to TEs were first selected and verified that the reads that map on the antisense strand predominantly contained an uridine (U) as the first nucleotide whereas the reads that map on the sense strand contained an adenosine (A) at the 10th position. This indicated that also piRNAs that derive from TE in mosquito are produced through a secondary ping-pong pathway as in Drosophila . Next, only reads containing a U as first nucleotide or an adenosine A at the 10th position were bioinformatically selected from each sample in order to enrich the reads of piRNAs. Clusters of piRNAs were than identified using the SiLoCo implementation from the UEA small RNA Workbench , considering chromosomal regions no larger than 100 base pairs (bp), loci containing average read size that range from 26 to 27 and containing more than 50 reads.
piRNA motif analysis
Analysis of Drosophila publicly available samples
Analysis of Drosophila derived sequence reads was performed on fastq files downloaded from NCBI Gene Expression Omnibus (GEO), (GSE6734, GSE15378, GSE11086) using the same tools and procedures described above for Anopheles.
Availability of supporting data
The datasets supporting the results of this article have been deposited at the European Nucleotide Archive with submission number PRJEB7896.
We are grateful to José Afonso Guerra-Assunção for initial attempts at microRNA prediction and Dan Lawson and Nikolai Windbichler for advice. This work was supported by grants from the European Commission FP7 projects INFRAVEC (grant agreement no. 228421) and the Foundation for the National Institutes of Health through the Vector-Based Control of Transmission: Discovery Research (VCTR) program of the Grand Challenges in Global Health initiative.
- Girard A, Sachidanandam R, Hannon GJ, Carmell MA. A germline-specific class of small RNAs binds mammalian Piwi proteins. Nature. 2006;442(7099):199–202.PubMedGoogle Scholar
- Brennecke J, Aravin AA, Stark A, Dus M, Kellis M, Sachidanandam R, et al. Discrete small RNA-generating loci as master regulators of transposon activity in Drosophila. Cell. 2007;128(6):1089–103.View ArticlePubMedGoogle Scholar
- Forstemann K, Tomari Y, Du T, Vagin VV, Denli AM, Bratu DP, et al. Normal microRNA maturation and germ-line stem cell maintenance requires Loquacious, a double-stranded RNA-binding domain protein. PLoS Biol. 2005;3(7):e236.View ArticlePubMed CentralPubMedGoogle Scholar
- Meister G. Argonaute proteins: functional insights and emerging roles. Nat Rev Genet. 2013;14(7):447–59.View ArticlePubMedGoogle Scholar
- Okamura K, Ishizuka A, Siomi H, Siomi MC. Distinct roles for Argonaute proteins in small RNA-directed RNA cleavage pathways. Genes Dev. 2004;18(14):1655–66.View ArticlePubMed CentralPubMedGoogle Scholar
- Saito K, Ishizuka A, Siomi H, Siomi MC. Processing of pre-microRNAs by the Dicer-1-Loquacious complex in Drosophila cells. PLoS Biol. 2005;3(7):e235.View ArticlePubMed CentralPubMedGoogle Scholar
- Jiang F, Ye X, Liu X, Fincher L, McKearin D, Liu Q. Dicer-1 and R3D1-L catalyze microRNA maturation in Drosophila. Genes Dev. 2005;19(14):1674–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu Q, Rand TA, Kalidas S, Du F, Kim HE, Smith DP, et al. R2D2, a bridge between the initiation and effector steps of the Drosophila RNAi pathway. Science. 2003;301(5641):1921–5.View ArticlePubMedGoogle Scholar
- Elbashir SM, Lendeckel W, Tuschl T. RNA interference is mediated by 21- and 22-nucleotide RNAs. Genes Dev. 2001;15(2):188–200.View ArticlePubMed CentralPubMedGoogle Scholar
- Saito K, Nishida KM, Mori T, Kawamura Y, Miyoshi K, Nagami T, et al. Specific association of Piwi with rasiRNAs derived from retrotransposon and heterochromatic regions in the Drosophila genome. Genes Dev. 2006;20(16):2214–22.View ArticlePubMed CentralPubMedGoogle Scholar
- Vagin VV, Sigova A, Li C, Seitz H, Gvozdev V, Zamore PD. A distinct small RNA pathway silences selfish genetic elements in the germline. Science. 2006;313(5785):320–4.View ArticlePubMedGoogle Scholar
- Gunawardane LS, Saito K, Nishida KM, Miyoshi K, Kawamura Y, Nagami T, et al. A slicer-mediated mechanism for repeat-associated siRNA 5’ end formation in Drosophila. Science. 2007;315(5818):1587–90.View ArticlePubMedGoogle Scholar
- Siomi MC, Sato K, Pezic D, Aravin AA. PIWI-interacting small RNAs: the vanguard of genome defence. Nat Rev Mol Cell Biol. 2011;12(4):246–58.View ArticlePubMedGoogle Scholar
- Campbell CL, Keene KM, Brackney DE, Olson KE, Blair CD, Wilusz J, et al. Aedes aegypti uses RNA interference in defense against Sindbis virus infection. BMC Microbiol. 2008;8:47.View ArticlePubMed CentralPubMedGoogle Scholar
- Myles KM, Morazzani EM, Adelman ZN. Origins of alphavirus-derived small RNAs in mosquitoes. RNA Biol. 2009;6(4):387–91.View ArticlePubMed CentralPubMedGoogle Scholar
- Myles KM, Wiley MR, Morazzani EM, Adelman ZN. Alphavirus-derived small RNAs modulate pathogenesis in disease vector mosquitoes. Proc Natl Acad Sci U S A. 2008;105(50):19938–43.View ArticlePubMed CentralPubMedGoogle Scholar
- Morazzani EM, Wiley MR, Murreddu MG, Adelman ZN, Myles KM. Production of virus-derived ping-pong-dependent piRNA-like small RNAs in the mosquito soma. PLoS Pathog. 2012;8(1):e1002470.View ArticlePubMed CentralPubMedGoogle Scholar
- Hess AM, Prasad AN, Ptitsyn A, Ebel GD, Olson KE, Barbacioru C, et al. Small RNA profiling of Dengue virus-mosquito interactions implicates the PIWI RNA pathway in anti-viral defense. BMC Microbiol. 2011;11:45.View ArticlePubMed CentralPubMedGoogle Scholar
- Mead EA, Tu Z. Cloning, characterization, and expression of microRNAs from the Asian malaria mosquito. Anopheles stephensi. BMC Genomics. 2008;9:244.View ArticlePubMed CentralPubMedGoogle Scholar
- Mead EA, Li M, Tu Z, Zhu J. Translational regulation of Anopheles gambiae mRNAs in the midgut during Plasmodium falciparum infection. BMC Genomics. 2012;13:366.View ArticlePubMed CentralPubMedGoogle Scholar
- Thirugnanasambantham K, Hairul-Islam VI, Saravanan S, Subasri S, Subastri A. Computational approach for identification of Anopheles gambiae miRNA involved in modulation of host immune response. Appl Biochem Biotechnol. 2013;170(2):281–91.View ArticlePubMedGoogle Scholar
- Winter F, Edaye S, Huttenhofer A, Brunel C. Anopheles gambiae miRNAs as actors of defence reaction against Plasmodium invasion. Nucleic Acids Res. 2007;35(20):6953–62.View ArticlePubMed CentralPubMedGoogle Scholar
- Chatterjee R, Chaudhuri K. An approach for the identification of microRNA with an application to Anopheles gambiae. Acta Biochim Pol. 2006;53(2):303–9.PubMedGoogle Scholar
- Jain S, Rana V, Shrinet J, Sharma A, Tridibes A, Sunil S, et al. Blood feeding and Plasmodium infection alters the miRNome of Anopheles stephensi. PLoS One. 2014;9(5):e98402.View ArticlePubMed CentralPubMedGoogle Scholar
- Hoa NT, Keene KM, Olson KE, Zheng L. Characterization of RNA interference in an Anopheles gambiae cell line. Insect Biochem Mol Biol. 2003;33(9):949–57.View ArticlePubMedGoogle Scholar
- Borsani O, Zhu J, Verslues PE, Sunkar R, Zhu JK. Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell. 2005;123(7):1279–91.View ArticlePubMed CentralPubMedGoogle Scholar
- Elbarbary RA, Takaku H, Uchiumi N, Tamiya H, Abe M, Takahashi M, et al. Modulation of gene expression by human cytosolic tRNase Z(L) through 5’-half-tRNA. PLoS One. 2009;4(6):e5908.View ArticlePubMed CentralPubMedGoogle Scholar
- Lee YS, Shibata Y, Malhotra A, Dutta A. A novel class of small RNAs: tRNA-derived RNA fragments (tRFs). Genes Dev. 2009;23(22):2639–49.View ArticlePubMed CentralPubMedGoogle Scholar
- Papathanos PA, Windbichler N, Menichelli M, Burt A, Crisanti A. The vasa regulatory region mediates germline expression and maternal transmission of proteins in the malaria mosquito Anopheles gambiae: a versatile tool for genetic control strategies. BMC Mol Biol. 2009;10:65.View ArticlePubMed CentralPubMedGoogle Scholar
- Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40(1):37–52.View ArticlePubMed CentralPubMedGoogle Scholar
- Stocks MB, Moxon S, Mapleson D, Woolfenden HC, Mohorianu I, Folkes L, et al. The UEA sRNA workbench: a suite of tools for analysing and visualizing next generation sequencing microRNA and small RNA datasets. Bioinformatics. 2012;28(15):2059–61.View ArticlePubMed CentralPubMedGoogle Scholar
- Castellano L, Stebbing J. Deep sequencing of small RNAs identifies canonical and non-canonical miRNA and endogenous siRNAs in mammalian somatic tissues. Nucleic Acids Res. 2013;41(5):3339–51.View ArticlePubMed CentralPubMedGoogle Scholar
- Ladewig E, Okamura K, Flynt AS, Westholm JO, Lai EC. Discovery of hundreds of mirtrons in mouse and human small RNA data. Genome Res. 2012;22(9):1634–45.View ArticlePubMed CentralPubMedGoogle Scholar
- Chiang HR, Schoenfeld LW, Ruby JG, Auyeung VC, Spies N, Baek D, et al. Mammalian microRNAs: experimental evaluation of novel and previously annotated genes. Genes Dev. 2010;24(10):992–1009.View ArticlePubMed CentralPubMedGoogle Scholar
- Ruby JG, Jan CH, Bartel DP. Intronic microRNA precursors that bypass Drosha processing. Nature. 2007;448(7149):83–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Westholm JO, Ladewig E, Okamura K, Robine N, Lai EC. Common and distinct patterns of terminal modifications to mirtrons and canonical microRNAs. RNA. 2012;18(2):177–92.View ArticlePubMed CentralPubMedGoogle Scholar
- Piriyapongsa J, Jordan IK. A family of human microRNA genes from miniature inverted-repeat transposable elements. PLoS One. 2007;2(2):e203.View ArticlePubMed CentralPubMedGoogle Scholar
- Piriyapongsa J, Jordan IK. Dual coding of siRNAs and miRNAs by plant transposable elements. RNA. 2008;14(5):814–21.View ArticlePubMed CentralPubMedGoogle Scholar
- Sun J, Zhou M, Mao Z, Li C. Characterization and evolution of microRNA genes derived from repetitive elements and duplication events in plants. PLoS One. 2012;7(4):e34092.View ArticlePubMed CentralPubMedGoogle Scholar
- Smalheiser NR, Torvik VI. Mammalian microRNAs derived from genomic repeats. Trends Genet. 2005;21(6):322–6.View ArticlePubMedGoogle Scholar
- Yuan Z, Sun X, Liu H, Xie J. MicroRNA genes derived from repetitive elements and expanded by segmental duplication events in mammalian genomes. PLoS One. 2011;6(3):e17666.View ArticlePubMed CentralPubMedGoogle Scholar
- Biryukova I, Ye T, Levashina E. Transcriptome-wide analysis of microRNA expression in the malaria mosquito Anopheles gambiae. BMC Genomics. 2014;15:557.View ArticlePubMed CentralPubMedGoogle Scholar
- Anders S, McCarthy DJ, Chen Y, Okoniewski M, Smyth GK, Huber W, et al. Count-based differential expression analysis of RNA sequencing data using R and Bioconductor. Nat Protoc. 2013;8(9):1765–86.View ArticlePubMedGoogle Scholar
- Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.View ArticlePubMed CentralPubMedGoogle Scholar
- Kugler JM, Verma P, Chen YW, Weng R, Cohen SM. miR-989 is required for border cell migration in the Drosophila ovary. PLoS One. 2013;8(7):e67075.View ArticlePubMed CentralPubMedGoogle Scholar
- Turner JM. Meiotic sex chromosome inactivation. Development. 2007;134(10):1823–31.View ArticlePubMedGoogle Scholar
- Song R, Ro S, Michaels JD, Park C, McCarrey JR, Yan W. Many X-linked microRNAs escape meiotic sex chromosome inactivation. Nat Genet. 2009;41(4):488–93.View ArticlePubMed CentralPubMedGoogle Scholar
- Marco A. Sex-biased expression of microRNAs in Drosophila melanogaster. Open Biol. 2014;4:140024.View ArticlePubMed CentralPubMedGoogle Scholar
- Babiarz JE, Ruby JG, Wang Y, Bartel DP, Blelloch R. Mouse ES cells express endogenous shRNAs, siRNAs, and other Microprocessor-independent, Dicer-dependent small RNAs. Genes Dev. 2008;22(20):2773–85.View ArticlePubMed CentralPubMedGoogle Scholar
- Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP. The impact of microRNAs on protein output. Nature. 2008;455(7209):64–71.View ArticlePubMed CentralPubMedGoogle Scholar
- Czech B, Malone CD, Zhou R, Stark A, Schlingeheyde C, Dus M, et al. An endogenous small interfering RNA pathway in Drosophila. Nature. 2008;453(7196):798–802.View ArticlePubMed CentralPubMedGoogle Scholar
- Ghildiyal M, Seitz H, Horwich MD, Li C, Du T, Lee S, et al. Endogenous siRNAs derived from transposons and mRNAs in Drosophila somatic cells. Science. 2008;320(5879):1077–81.View ArticlePubMed CentralPubMedGoogle Scholar
- Robine N, Lau NC, Balla S, Jin Z, Okamura K, Kuramochi-Miyagawa S, et al. A broadly conserved pathway generates 3’UTR-directed primary piRNAs. Curr Biol. 2009;19(24):2066–76.View ArticlePubMed CentralPubMedGoogle Scholar
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J. Repbase update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005;110(1–4):462–7.View ArticlePubMedGoogle Scholar
- O’Shea JP, Chou MF, Quader SA, Ryan JK, Church GM, Schwartz D. pLogo: a probabilistic approach to visualizing sequence motifs. Nat Methods. 2013;10(12):1211–2.View ArticlePubMedGoogle Scholar
- Lau NC, Robine N, Martin R, Chung WJ, Niki Y, Berezikov E, et al. Abundant primary piRNAs, endo-siRNAs, and microRNAs in a Drosophila ovary cell line. Genome Res. 2009;19(10):1776–85.View ArticlePubMed CentralPubMedGoogle Scholar
- Frank F, Sonenberg N, Nagar B. Structural basis for 5’-nucleotide base-specific recognition of guide RNA by human AGO2. Nature. 2010;465(7299):818–22.View ArticlePubMedGoogle Scholar
- Ghildiyal M, Xu J, Seitz H, Weng Z, Zamore PD. Sorting of Drosophila small silencing RNAs partitions microRNA* strands into the RNA interference pathway. RNA. 2010;16(1):43–56.View ArticlePubMed CentralPubMedGoogle Scholar
- Hu HY, Yan Z, Xu Y, Hu H, Menzel C, Zhou YH, et al. Sequence features associated with microRNA strand selection in humans and flies. BMC Genomics. 2009;10:413.View ArticlePubMed CentralPubMedGoogle Scholar
- Lau NC, Lim LP, Weinstein EG, Bartel DP. An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans. Science. 2001;294(5543):858–62.View ArticlePubMedGoogle Scholar
- Ivanov P, Emara MM, Villen J, Gygi SP, Anderson P. Angiogenin-Induced tRNA Fragments Inhibit Translation Initiation. Mol Cell. 2011;43(4):613–23.View ArticlePubMed CentralPubMedGoogle Scholar
- Haussecker D, Huang Y, Lau A, Parameswaran P, Fire AZ, Kay MA. Human tRNA-derived small RNAs in the global regulation of RNA silencing. Rna-Publication Rna Soc. 2010;16(4):673–95.View ArticleGoogle Scholar
- Malone CD, Brennecke J, Dus M, Stark A, McCombie WR, Sachidanandam R, et al. Specialized piRNA pathways act in germline and somatic tissues of the Drosophila ovary. Cell. 2009;137(3):522–35.View ArticlePubMed CentralPubMedGoogle Scholar
- Li C, Vagin VV, Lee S, Xu J, Ma S, Xi H, et al. Collapse of germline piRNAs in the absence of Argonaute3 reveals somatic piRNAs in flies. Cell. 2009;137(3):509–21.View ArticlePubMed CentralPubMedGoogle Scholar
- Macias V, Coleman J, Bonizzoni M, James AA. piRNA pathway gene expression in the malaria vector mosquito Anopheles stephensi. Insect Mol Biol. 2014;23(5):579–86.View ArticlePubMed CentralPubMedGoogle Scholar
- Gou LT, Dai P, Yang JH, Xue YC, Hu YP, Zhou Y, et al. Pachytene piRNAs instruct massive mRNA elimination during late spermiogenesis. Cell Res. 2014;24(6):680–700.View ArticlePubMed CentralPubMedGoogle Scholar
- Rouget C, Papin C, Boureux A, Meunier AC, Franco B, Robine N, et al. Maternal mRNA deadenylation and decay by the piRNA pathway in the early Drosophila embryo. Nature. 2010;467(7319):1128–32.View ArticlePubMedGoogle Scholar
- Aravin AA, Klenov MS, Vagin VV, Bantignies F, Cavalli G, Gvozdev VA. Dissection of a natural RNA silencing process in the Drosophila melanogaster germ line. Mol Cell Biol. 2004;24(15):6742–50.View ArticlePubMed CentralPubMedGoogle Scholar
- Kiuchi T, Koga H, Kawamoto M, Shoji K, Sakai H, Arai Y, et al. A single female-specific piRNA is the primary determiner of sex in the silkworm. Nature. 2014;509(7502):633–6.View ArticlePubMedGoogle Scholar
- Krzywinski J, Nusskern DR, Kern MK, Besansky NJ. Isolation and characterization of Y chromosome sequences from the African malaria mosquito Anopheles gambiae. Genetics. 2004;166(3):1291–302.View ArticlePubMed CentralPubMedGoogle Scholar
- James AA. Gene drive systems in mosquitoes: rules of the road. Trends Parasitol. 2005;21(2):64–7.View ArticlePubMedGoogle Scholar
- Chen CF, Ridzon DA, Broomer AJ, Zhou ZH, Lee DH, Nguyen JT, et al. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005;33(20):e179.View ArticlePubMed CentralPubMedGoogle Scholar
- Crooks GE, Hon G, Chandonia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004;14(6):1188–90.View ArticlePubMed CentralPubMedGoogle Scholar
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.