- Research article
- Open Access
Bioinformatic prediction, deep sequencing of microRNAs and expression analysis during phenotypic plasticity in the pea aphid, Acyrthosiphon pisum
© Legeai et al; licensee BioMed Central Ltd. 2010
- Received: 7 July 2009
- Accepted: 5 May 2010
- Published: 5 May 2010
Post-transcriptional regulation in eukaryotes can be operated through microRNA (miRNAs) mediated gene silencing. MiRNAs are small (18-25 nucleotides) non-coding RNAs that play crucial role in regulation of gene expression in eukaryotes. In insects, miRNAs have been shown to be involved in multiple mechanisms such as embryonic development, tissue differentiation, metamorphosis or circadian rhythm. Insect miRNAs have been identified in different species belonging to five orders: Coleoptera, Diptera, Hymenoptera, Lepidoptera and Orthoptera.
We developed high throughput Solexa sequencing and bioinformatic analyses of the genome of the pea aphid Acyrthosiphon pisum in order to identify the first miRNAs from a hemipteran insect. By combining these methods we identified 149 miRNAs including 55 conserved and 94 new miRNAs. Moreover, we investigated the regulation of these miRNAs in different alternative morphs of the pea aphid by analysing the expression of miRNAs across the switch of reproduction mode. Pea aphid microRNA sequences have been posted to miRBase: http://microrna.sanger.ac.uk/sequences/
Our study has identified candidates as putative regulators involved in reproductive polyphenism in aphids and opens new avenues for further functional analyses.
- Deep Sequencing
- Microfluidic Chip
- Mature miRNAs
- miRNA Precursor
- Parthenogenetic Female
MicroRNAs (miRNAs) are small (18-24 nucleotides) non-coding RNAs (ncRNAs) that regulate gene expression in eukaryotes. Numerous miRNA genes have been found in animal genomes. They are located either within intronic sequences of mRNA-coding genes, or in intergenic regions. Many miRNAs are highly conserved throughout evolution (reviewed in ). However, there are also taxa-specific miRNAs [2–5]. One such miRNA is iab-4 which is only described in insect species and is involved in wing formation . The description of insect miRNAs in miRBase  remains largely restricted to Diptera (Drosophila melanogaster, D. pseudoobscura, Anopheles gambiae), Hymenoptera (Apis mellifera, Nasonia vitripennis), Coleoptera (Tribolium castaneum), Orthoptera (Locusta migratoria) and Lepidoptera (Bombyx morii), which all diverged about 280 million-years ago . The insect species with the greatest number of miRNAs (147) represented in miRBase is D. melanogaster, which until recently was also the only insect for which deep sequencing results have been combined with a thorough bioinformatics analysis . This number has recently been surpassed in a study on the silkworm, B. mori  and the locust Locusta migratoria . Functional analyses of insect miRNAs has been mainly restricted to D. melanogaster [9, 12] where they have been shown to affect multiple biological processes such as embryo development and tissue differentiation, cell proliferation or morphological changes. Until now, no miRNAs have been described from the Hemiptera, a group of insects that includes many of the world's most damaging insect pests such as aphids, whiteflies, and scales.
MiRNAs are generated from genome-encoded precursors that form hairpin structures with imperfect base-paired segments. MiRNAs processing occurs in several steps (reviewed in [13, 14]). The primary transcripts (pri-miRNA) are essentially synthesised by RNA polymerase II and cleaved in the nucleus by the RNAse III Drosha. The resulting 70 nucleotide-hairpin pre-miRNA is exported in the cytoplasm and the final miRNA maturation step is mediated by the RNAse III Dicer that produces 18-24 nucleotide long miRNA duplexes. One strand of this duplex is incorporated into the RNA-induced silencing complex (miRISC) and guides this miRISC to degrade, destabilize or translationally inhibit the mRNA targets (reviewed in ). Post-transcriptional regulation of gene expression by miRNAs is highly complex since a single miRNA can regulate hundred of genes, and a single gene may be regulated by multiple miRNAs.
The newly assembled genome of the pea aphid Acyrthosiphon pisum recently released by the International Aphid Genomics Consortium  has enabled the identification of miRNAs in a hemipteran insect for the first time. Aphids are herbivorous insects characterised by the unique ability to feed exclusively on phloem sap. The serious damage caused by aphids is partly due to their amazing ability to adapt to environmental variations . This adaptive capacity is largely explained by their phenotypic plasticity that allows the production of distinct phenotypic morphs in response to environmental changes. This is illustrated by the switch of reproductive mode in response to seasonal changes: parthenogenetic females are produced in spring and summer whereas sexual females and males are produced in autumn [reviewed in ]. The pea aphid also shows a unique duplication of the miRNA processing machinery suggesting that miRNAs play a significant role in the life cycle of these insects . Here, we present a list of 149 A. pisum miRNAs identified by combining in vivo and bioinformatic approaches. Moreover, the regulation of these miRNAs in phenotypic plasticity has been investigated and miRNAs known to be regulated by insect endocrine pathways were identified as differentially expressed in different morphs involved in sexual and asexual reproduction.
Identification of A. pisum miRNAs
The recently sequenced and assembled A. pisum genome  was used to identify miRNAs in the pea aphid using three complementary strategies. First, miRNAs listed in miRBase were used to identify pea aphid miRNAs by sequence homology. Second, small RNAs extracted from parthenogenetic females of the pea aphid were sequenced by deep sequencing technology and analysed using the miRDeep program . Finally, the pea aphid genome was screened for putative miRNA precursors using a new algorithm, GR4500.
In the first approach we blasted 1275 insect miRNAs (D. melanogaster and 11 other Drosophila species, A. gambiae, B. mori, T. castaneum, A. melifera, Locusta migratoria) deposited in miRBase release 14 against the A. pisum genome. The 200 nucleotide-long genomic sequences surrounding each hit were investigated for their ability to fold into typical miRNA precursor hairpins. In total 43 different putative aphid miRNAs showed both sequence homology and a potential hairpin structure (Additional file 1). These 43 A. pisum miRNAs candidates corresponded to 44 precursors. The expression of 33 of them was verified by RT-PCR (data not shown).
In the second approach we undertook deep sequencing of small RNAs from a mixed generation sample of parthenogenetic female pea aphids. Approximately 3 million sequences were generated, corresponding to approximately 850,000 unique sequences. MiRNAs are known to show obvious size preference and tend to be between 18 and 24 nucleotides in length. Analysis of the length distribution of the pea aphid mature miRNAs showed a peak at 22nt as previously reported in D. melanogaster. Moreover, similar to other organisms, many of the predicted pea aphid miRNA mature sequences (59%) start with a uridine residue at their 5'end [20–22].
MiRDeep is a selective program for identifying miRNAs from deep sequencing data that minimises the false positive rate. It selects reads that align exactly to the genome, at less than 5 different locations, and selects putative precursors of miRNA that match a high number of reads on the hairpin structure where the mature miRNA is supposed to be. Using a pool of several 454 sequencing runs, Friedländer et al.  showed that miRDeep was able to retrieve 89% of the known miRNA from Caenorhabditis elegans, and recovered 73% of the known miRNA incorporated in a set of 106 Solexa sequences from a HeLa cell small RNA sample sequenced by Solexa technology. In addition to miRDeep, in order to increase the number of predicted pea aphid miRNAs, we implemented a complementary approach, called GR4500. This algorithm first screened the pea aphid genome for miRNA putative precursors (hairpin structure, see methods) using a classifier trained with a set of 30 validated pea aphid miRNAs. GR4500 then compared the selected hairpins to the Solexa reads for biological validation and for identifying the mature sequence. The classifier selected 4402 candidates from 2.5 million genomic hairpins of more than 63 nucleotides. The comparison of these candidates with the 850,000 sequences of small RNAs previously obtained by deep sequencing confirmed the expression of 116 of them. These 116 precursors corresponded to 98 different mature miRNAs (Additional file 1). GR4500 identified 38 new mature miRNAs not found in miRBase or identified by miRDeep from 45 genomic precursors. None of these 38 miRNAs have orthologues in other organisms.
Altogether, the combination of homology search, miRDeep analysis of deep sequencing and GR4500 scan of the genome allowed the identification of 149 mature miRNAs candidates in the pea aphid (Fig 1). These 149 mature miRNAs candidates correspond to 176 precursors of miRNAs (pre-miRNAs) (Additional file 1). This is similar to the number of annotated miRNAs on miRBase for D. melanogaster (157 miRNAs). For all these pea aphid miRNA candidates, we have evidence of their expression since mature sequences have been detected either by RT-PCR and/or deep sequencing. It is admitted that microRNAs identified by deep sequencing, have to be supported by multiple reads and/or by the existence of its complementary strand mir* (Meyers et al. 2008). Applying this rule to the pea aphid predicted miRNA by miRDeep and GR4500 lead to the subdivision between 98 A. pisum miRNAs (≥ 5 reads and/or mir* detected) and 51 miRNA candidates (< 5 reads without mir*) (Additional file 2). These miRNA candidates still require confirmation that will be provided by future extensive deep sequencing, such as in other morphs than parthenogenetic individuals used for this study.
Genomic organisation of miRNA precursors
As expected, the majority of miRNA precursors (86%) were located in intergenic regions (97 precursors) or introns (55 precursors) With the exception of two miRNAs located in unassembled reads (raw sequences from the genome of the pea aphid not assembled in the first release of the genome, see IAGC et al. 2009) (Additional file 1), the remainders of the miRNA precursors were located in regions annotated as exons from protein coding genes (12.5%, 22). However, only two of these correspond to predicted genes whose annotation is supported with biological evidence (ESTs), encoding a heterogeneous nuclear ribonucleoprotein K and a vesicular mannose-binding lectin. Exonic miRNAs have been previously described in mammals [27, 28]. However, with a few rare exceptions, exons that encode miRNAs do not code for proteins and are "mRNA-like non coding RNAs". The remaining 20 exonic miRNAs of the pea aphid are found in predicted genes whose annotation is not yet supported by biological evidence. Moreover, among the 20 miRNAs identified within exon, 15 were well supported either by abundant Solexa reads and/or the existence of mir*. This could indicate false prediction in the protein coding gene set, or that these predicted genes correspond to mRNA-like non coding RNAs.
We investigated the physical distribution of the 176 precursors of miRNAs along the different assembled genomic scaffolds of the pea aphid genome. 71.6% (126) of miRNA precursors were identified as singletons (1 miRNA locus per scaffold) while 28.4% (50) of miRNAs precursors were distributed in 17 clusters composed of up to 6 miRNAs per cluster (Additional file 1). Among the 17 pea aphid miRNAs clusters, 3 clusters are composed of the same identical mature miRNAs present as multicopy (Ap-mir-X40, Ap-mir-971, Ap-mir-X5) with slight differences within precursor sequences, suggesting a very recent duplication of these miRNAs. We identified 4 clusters composed of closely related miRNAs (e.g. Ap-mir92a and Ap-mir-92b) with slight differences in their precursor or mature sequences, and 10 clusters composed of different mature miRNAs, suggesting more ancient duplication events. The length of these clusters varies from 160 nucleotides (in that case, the miRNA precursors overlap) to 30 kb.
Two aphid miRNA clusters are unusually long: cluster EQ127026-cl1 (Ap-mir-277/Ap-mir-317/Ap-mir-34) at 15 kb and cluster EQ127560-cl1 (Ap-mir-263b/Ap-mir-228) at 30. This EQ127026-cl1 cluster is also conserved in D. melanogaster and A. mellifera. Moreover, long length miRNA clusters up to 62 kb have been identified in the mouse genome . Cluster EQ119865-cl1 (Ap-mir-2a-1a/Ap-mir-2a-1b/Ap-mir-2b/Ap-mir-13a/Ap-mir-13b/Ap-mir-71), and cluster EQ127026-cl1 also show similar organisation between the honey bee and the pea aphid genomes. However, some clusters conserved in D. melanogaster and honey bee have diverged in aphids. The duplication of Api-mir-307 observed in the pea aphid genome has not been reported in other insect genomes. The cluster including let-7, mir-100 and mir-125 conserved in D. melanogaster and in A. mellifera is incomplete in A. pisum, since Api-mir-125 was not found. Sequence analysis of the scaffold EQ112277-cl1 containing Ap-mir-100 and Ap-let-7 showed a sequence partially homologous to mir-125 (18 bases identical on 22). However, the genomic sequence surrounding this putative mir-125 sequence did not show any hairpin structure. These variations in clusters organisation could indicate an adaptation for each insect order of their miRNA gene contents possibly related to life history traits.
Functional annotation of pea aphid miRNAs
Expression of miRNAs in alternative pea aphid morphs
Aphids have a complicated life cycle characterized by a phase of viviparous clonal reproduction (parthenogenesis) that alternates with a phase of sexual reproduction. This switch is triggered by a sensing of the decrease of the day length in autumn (reviewed in ). The different alternative morphs produced are: parthenogenetic females named virginoparae producing viviparous parthenogenetic females, oviparae sexual females, sexual males, and parthenogenetic females named sexuparae producing sexual males and females. The mechanisms by which aphids detect and respond to differences in day length are not known. MiRNAs have been implicated in modulating circadian rhythm responses in insects . On this basis, it is reasonable to hypothesize that miRNAs may also contribute to the regulation of aphid polyphenism.
Interestingly, among them, Ap-let-7, Ap-mir2a-1 and Ap-mir100 showed differential expression between oviparae and the two other parthenogenetic morphs (virginoparae and sexuparae). Ap-let-7 and Ap-mir-100 were up-regulated in oviparae, and Ap-mir2a-1 was down-regulated in oviparae. Ap-let-7 and Ap-mir100 belong to the same cluster and follow similar expression patterns. Mir-2 is involved in apoptosis regulation during development  and let-7 and mir-100 have been reported to be involved in metamorphosis  and are up-regulated in response to ecdysone , a hormone involved in insect development. Three miRNAs (Ap-mir-34, Ap-mir-X47 and Ap-mir-X103) and two mir* (Ap-mir307* and Ap-mirX-52*) showed different expression between the two parthenogenetic morphs: sexuparae and virginoparae, that differ by the type of embryos they contain (sexual vs asexual), and by the conditions of rearing (short-day vs long-day). Four of these miRNAs were specific to the virginoparae/sexuparae comparison, while Ap-mir-X103 was also up-regulated in oviparae. These 5 miRNAs are particularly interesting candidates for the switch of reproduction mode from parthenogenesis to sexual reproduction. Three of these miRNAs are newly-identified aphid miRNAs and their function remains to be determined. Interestingly, mir-34 has been shown to be regulated in D. melanogaster by ecdysone as well as by juvenile hormones . Juvenile hormones are known to be involved in aphids in the transduction of the photoperiodic signal from the brain to the ovaries during the switch of reproductive mode . Thus, our study has identified a strong candidate (Ap-mir-34) as a putative regulator involved in reproductive polyphenism in aphids and opens new avenues for further functional analyses.
This work has established a catalog of miRNA genes in the pea aphid that represents an essential base of knowledge-base for investigating the miRNA post-transcriptional regulation of key biological traits for an organism whose adaptation is shaped by phenotypic plasticity. Deciphering the gene regulation network between miRNAs and their mRNA targets in the pea aphid remains an objective that this study opens.
The LSR1 clone  of the pea aphid Acyrthosiphon pisum was reared and maintained as clonal individuals (parthenogenesis) on the plant Vicia fabae at 18°C under a 16 h photoperiod for long-day condition (parthenogenesis). Sexual individuals and sexuparae were produced by rearing the pea aphid at 18°C under a 12 h photoperiod (short-day condition) as described in .
Homology identification of A. pisum miRNAs
For homology prediction of A. pisum miRNAs, insect miRNAs (D. melanogaster and the 11 other Drosophila species, A. mellifera, A. gambiae, B. mori, T. castaneum and Locusta migratoria) were retrieved from miRBase release 14  and blasted against the pea aphid genome. Putative miRNAs were selected as sequences identical or slightly identical (with 1 or 2 nt different) with the original homologous mature miRNA sequence retrieved from miRBase. A sliding window of 150 bases of genomic sequences surrounding each putative miRNAs were retrieved and its ability to fold into a potential miRNA hairpin precursor was investigated using the RNA fold program  and/or the MiRAlign program ). Secondary structure should have free energy change less than or equal to -18 kcal/mole.
To investigate the conservation of mature aphid miRNAs with other miRNAs listed in miRBase was analysed using the Similarity SEARCH program (SSEARCH) available on the miRBase website with default parameters. A minimum SSEARCH score of 135 was selected to be significant.
Small RNA sequencing
Total RNA was extracted from a mixed culture of parthenogenetic females of the LSR1 clone  of pea aphid using PureZOL RNA Isolation Reagent (Bio-Rad, Hayward CA). The RNA concentration and purity were determined photometrically by measuring absorbance at 260 nm and A260/A280 ratio using the NanoDrop ND-1000 spectrophotometer (Nanodrop Technologies). 40 μg of ethanol precipitated RNA was sent in duplicate to Illumina Inc. (Hayward CA) for size fractionation (<50 bp) and deep sequencing. Results of small RNA sequencing have been deposited in Gene Expression Omnibus (Additional file 4).
Mapping of the short reads Illumina sequences on the pea aphid genome
Approximately 3 million total sequences were processed to remove linker sequences and quantify the resulting unique sequences. The 851,979 unique sequences (mean: 25.2 bp, median: 26 bp, max: 33 bp, min: 1 bp) were compared to the A. pisum genome using NCBI megablast (wod size: 12, no filter). For further analysis, we used an extremely selective process by conserving only exact match on the complete short read. Sequences that mapped to more than 5 different places on the genome were removed. It resulted in 305,055 alignments, covering only 205,500 (24%) of the unique reads, but 44% (1,333,398) of all reads.
MiRDeep analysis was done according to the miRDeep protocol, and a data set of putative precursor sequences was obtained from the genome, based on the predicted mature miRNA alignments. The secondary structures of the putative precursors were predicted by RNAfold . Short reads were aligned to their putative genomic precursor to create signatures of the precursors. Finally, we run miRDeep with the RNAfold stability test using those structures and signatures inputs.
The pea aphid genome was scanned to find miRNA-like hairpins. For that, we first used RNAfold  on 120 nucleotides windows with an overlap of 100 nucleotides, and selected in each window all hairpins longer than 63 bases with at least 70% arm base pairing. This led to approximately 2.5 million genomic hairpins. We then developed a set of features to discriminate between miRNA and non miRNA hairpins: folding energy, total and maximum internal loop size and symmetry, arms and terminal loop sizes; GC%, complexity score (in terms of repetitions of small motifs) and the score delivered by the Microprocessor SVM  which characterizes the presence of a Drosha recognition site. Using 38 pea aphid miRNAs already found by homology search as a learning set, we constructed a classifier using all these features. We first used a decision tree implemented in the R tree package to build a prototype classifier and then hand-tuned it to improve its selectivity. Additional file 5 lists the features and their corresponding cut-off thresholds. 4402 hairpins were selected and used for short reads alignment to create signatures of the precursors. Finally, as a last step, MiRDeep parametered with the sensitive option (-x) was used to extract mature sequences using those signatures and structures.
MiRNAs profiling in three different morphs of the pea aphid
Based on the sequence of aphid mature miRNAs, we designed probes. The custom μparaflo™ microfluidic chip (LC Sciences Houston, USA) contained one probe for each mature aphid miRNA and its corresponding mir*. Each probe was repeated 10 times on the chip to ensure assay reproducibility. Multiple control probes were included on each chip for quality controls, sample labelling and assay conditions. Among the control probes, PUC2PM-20B and PUC2 MM-20B are the perfect match and single-based match detection probes, respectively, of a 20-mer RNA positive control sequence that is spiked into the RNA samples before labelling. The detection probes were made by in situ synthesis using PGR (photogenerated reagent) chemistry. The expression of miRNAs was analysed in three different adult morphs: parthenogenetic females reared under long day condition that produce parthenogenetic embryos, parthenogenetic females reared under short day condition that produce sexual embryos, and sexual females reared under short day condition that produce eggs. Total RNAs were extracted from 10 adult females 48 h after the final moult. Three independent samples were analysed for each reproductive morph. Total RNAs were extracted using miRNeasy purification kit (QIAGEN) according to manufacturer's instructions. Quality of RNAs was checked with the Bioanalyser (Agilent). The assay started from 5 μg total RNA sample, which was size fractionated using a YM-100 Microcon centrifugal filter (Millipore) and the small RNAs (< 300 bases) isolated were 3'-extended with a poly(A) tail using poly(A) polymerase. An oligonucleotide tag was then ligated to the poly(A) tail for later fluorescent dye staining; two different tags were used for the two RNA samples in dual-sample experiments. Hybridization was performed overnight on a μParaflo microfluidic chip using a micro-circulation pump (Atactic Technologies). On the microfluidic chip, each detection probe consisted of a chemically modified nucleotide coding segment complementary to target microRNA or other control RNAs and a spacer segment of polyethylene glycol to extend the coding segment away from the substrate. The hybridization melting temperatures were balanced by chemical modifications of the detection probes. Hybridization used 100 μL 6× SSPE buffer (0.90 M NaCl, 60 mM Na2HPO4, 6 mM EDTA, pH 6.8) containing 25% formamide at 34°C. After RNA hybridization, tag-conjugating Cy3 were circulated through the microfluidic chip for dye staining. Fluorescence images were collected using a laser scanner (GenePix 4000B, Molecular Device) and digitized using Array-Pro image analysis software (Media Cybernetics).
For each chip and each probe, the average signal value and its standard deviation were quantified. Data were analyzed by first subtracting the background, then integrating all the signals corresponding for the same probe for one given chip. A transcript to be listed as detectable must meets at least two conditions: signal intensity higher than 3× (background standard deviation) and spot CV < 0.5. CV. When repeating probes are present on an array, a transcript is listed as detectable only if the signals from at least 50% of the repeating probes are above detection level. Then normalization the signals from all arrays were performed using a LOWESS filter (Locally-weighted Regression). Results obtained in the different reproductive morphs were compared by comparing the ratio of the two sets of detected signals (log2) and p-values of the t-test were calculated. Differentially detected signals were those with p-value < 0.05. The results of microfluidic experiments have been deposited in Gene Expression Omnibus (Additional file 4).
This work has been supported by INRA-SPE (AIP Séquençage, Programme innovant, AIP INRA/INRIA), ANR Holocentrism and the CSIRO OCE postdoctoral fellowship program. We thank the sequencing platform of Biogenouest (Roscoff, France) for their help at the beginning of the project. The authors would also like to thank Dr Shuji Shigenonu (Okazaki Institute for Integrative Bioscience, Japan) for his RNA extraction for Illumina/Solexa sequencing and Andrew Spriggs (CSIRO Plant Industry, Canberra, Australia) for his bioinformatic advice. Beverley Osborne is acknowledged for his writing assistance.
- Kim VN, Han J, Siomi MC: Biogenesis of small RNAs in animals. Nat Rev Mol Cell Biol. 2009, 10: 126-39. 10.1038/nrm2632.PubMedView ArticleGoogle Scholar
- Bentwich I, Avniel A, Karov Y, Aharonov R, Gilad S, Barad O, Barzilai A, Einat P, Einav U, Meiri E, Sharon E, Spector Y, Bentwich Z: Identification of hundreds of conserved and nonconserved human microRNAs. Nat Genet. 2005, 37: 766-770. 10.1038/ng1590.PubMedView ArticleGoogle Scholar
- Berezikov E, Guryev V, Belt van de J, Wienholds E, Plasterk RH, Cuppen E: Phylogenetic shadowing and computational identification of human microRNA genes. Cell. 2005, 120: 21-24. 10.1016/j.cell.2004.12.031.PubMedView ArticleGoogle Scholar
- Weaver DB, Anzola JM, Evans JD, Reid JG, Reese JT, Childs KL, Zdobnov EM, Samanta MP, Miller J, Elsik CG: Computational and transcriptional evidence for microRNAs in the honey bee genome. Genome Biol. 2007, 8: R97-10.1186/gb-2007-8-6-r97.PubMed CentralPubMedView ArticleGoogle Scholar
- Winter F, Edaye S, Huettenhofer A, Brunel C: Anopheles gambiae miRNAs as actors of defence reaction against Plasmodium invasion. Nucl Acid Res. 2007, 35: 6953-6962. 10.1093/nar/gkm686.View ArticleGoogle Scholar
- Ronshaugen M, Biemar , Piel J, Levine M, Lai EC: The Drosophila microRNA iab-4 causes a dominant homeotic transformation of halteres to wings. Genes Dev. 2005, 19: 2947-2952. 10.1101/gad.1372505.PubMed CentralPubMedView ArticleGoogle Scholar
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucl Acids Res. 2008, 36: D154-D158. 10.1093/nar/gkm952.PubMed CentralPubMedView ArticleGoogle Scholar
- Grimaldi D, Engels MS: Evolution of the Insects. 2005, Cambridge University Press: New YorkGoogle Scholar
- Lu J, Shen Y, Wu Q, Kumar S, He B, Shi S, Carthew RW, Wang SM, Wu CI: The birth and death of microRNA genes in Drosophila . Nat Genet. 2008, 40: 351-355. 10.1038/ng.73.PubMedView ArticleGoogle Scholar
- Zhang Y, Zhou X, Ge X, Jiang J, Li M, Jia S, Yang X, Kan Y, Miao X, Zhao G, Li F, Huang Y: Insect-specific microRNA involved in the development of the silkworm Bombyx mori . PLoS One. 2009, 4: e4677-10.1371/journal.pone.0004677.PubMed CentralPubMedView ArticleGoogle Scholar
- Wei Y, Chen S, Yang P, Ma Z, Kang L: Characterization and comparative profiling of the small RNA transcriptomes in two phases of locust. Genome Biol. 2009, 10: R6-10.1186/gb-2009-10-1-r6.PubMed CentralPubMedView ArticleGoogle Scholar
- Jaubert S, Méreau A, Antoniewski C, Tagu D: MicroRNAs in Drosophila : the magic wand to enter the chamber of secrets?. Biochimie. 2007, 89: 1211-1220. 10.1016/j.biochi.2007.05.012.PubMedView ArticleGoogle Scholar
- Carthew RW, Sontheimer EJ: Origins and Mechanisms of miRNAs and siRNAs. Cell. 2009, 136: 642-655. 10.1016/j.cell.2009.01.035.PubMed CentralPubMedView ArticleGoogle Scholar
- Winter J, Jung S, Keller S, Gregory RI, Diederichs S: Many roads to maturity: microRNA biogenesis pathways and their regulation. Nat Cell Biol. 2009, 11: 228-234. 10.1038/ncb0309-228.PubMedView ArticleGoogle Scholar
- International Aphid Genomics Consortium: Genome sequence of the pea aphid Acyrthosiphon pisum. PLoS Biol. 2010, 8: e1000313-10.1371/journal.pbio.1000313.View ArticleGoogle Scholar
- Tagu D, Klingler JP, Moya A, Simon JC: Early progress in aphid genomics and consequences for plant-aphid interactions studies. Mol Plant Microb Interac. 2008, 6: 701-708. 10.1094/MPMI-21-6-0701.View ArticleGoogle Scholar
- Le Trionnaire G, Hardie J, Jaubert-Possamai S, Simon JC, Tagu D: Shifting from asexual to sexual reproduction in aphids: physiological and developmental aspects. Biol Cell. 2008, 100: 441-451. 10.1042/BC20070135.PubMedView ArticleGoogle Scholar
- Jaubert-Possamai S, Rispe C, Tanguy S, Gordon K, Walsh T, Edwards O, Tagu D: Expansion of the miRNA pathway in the hemipteran insect Acyrthosiphon pisum . Mol Biol Evol. 2010, 27 (5): 979-87. 10.1093/molbev/msp256.PubMed CentralPubMedView ArticleGoogle Scholar
- Friedländer MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, Rajewsky N: Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotech. 2008, 26: 407-415. 10.1038/nbt1394.View ArticleGoogle Scholar
- Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, Tuschl T: Identification of tissue-specific microRNAs from mouse. Curr Biol. 2002, 12: 735-739. 10.1016/S0960-9822(02)00809-6.PubMedView ArticleGoogle Scholar
- Reinhart BJ, Bartel DP: Small RNAs correspond to centromere heterochromatic repeats. Science. 2002, 297: 1831-1831. 10.1126/science.1077183.PubMedView ArticleGoogle Scholar
- Singh J, Nagaraju J: In silico prediction and characterization of microRNAs from red flour beetle (Tribolium castaneum). Ins Mol Biol. 2008, 17: 427-436. 10.1111/j.1365-2583.2008.00816.x.View ArticleGoogle Scholar
- Soares AR, Pereira PM, Santos B, Egas C, Gomes AC, Arrais J, Oliveira JL, Moura GR, Santos MA: Parallel DNA pyrosequencing unveils new zebrafish microRNAs. BMC Genomics. 2009, 10: 195-10.1186/1471-2164-10-195.PubMed CentralPubMedView ArticleGoogle Scholar
- Glazov EA, Cottee PA, Barris WC, Moore RJ, Dalrymple BP, Tizard ML: A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach. Genome Res. 2008, 18: 957-964. 10.1101/gr.074740.107.PubMed CentralPubMedView ArticleGoogle Scholar
- Okamura K, Phillips MD, Tyler DM, Duan H, Chou Y, Lai E: The regulatory activity of microRNA* species has substantial influence on microRNA and 3'UTR evolution. Nat Struct Mol Biol. 2008, 15: 354-363. 10.1038/nsmb.1409.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu N, Okamura K, Phillips MD, Chung W, Lai E: The evolution and functional diversification of animal microRNA genes. Cell Res. 2008, 18: 985-996. 10.1038/cr.2008.278.PubMed CentralPubMedView ArticleGoogle Scholar
- Rodriguez A, Griffiths-Jones S, Ashurst JL, Bradley A: Identification of mammalian microRNA host genes and transcription units. Genome Res. 2004, 14: 1902-1910. 10.1101/gr.2722704.PubMed CentralPubMedView ArticleGoogle Scholar
- Weber MJ: New human and mouse microRNA genes found by homology search. FEBS J. 2005, 272: 59-73. 10.1111/j.1432-1033.2004.04389.x.PubMedView ArticleGoogle Scholar
- Ro S, Park C, Sanders KM, McCarrey JR, Yan W: Cloning and expression profiling of testis-expressed microRNAs. Dev Biol. 2007, 311: 592-602. 10.1016/j.ydbio.2007.09.009.PubMed CentralPubMedView ArticleGoogle Scholar
- Cheng HYM, Obrietan K: Revealing a role of microRNAs in the regulation of the biological clock. Cell Cycle. 2007, 6: 3034-3038.PubMedView ArticleGoogle Scholar
- Leaman D, Chen PY, Fak J, Yalcin A, Pearce M, Unnerstall U, Marks DS, Sander C, Tuschl T, Gaul U: Antisense-mediated depletion reveals essential and specific functions of microRNAs in Drosophila development. Cell. 2005, 121: 1097-1108. 10.1016/j.cell.2005.04.016.PubMedView ArticleGoogle Scholar
- Sempere LF, Dubrovsky EB, Dubrovskaya VA, Berger EM, Ambros V, Ros A: The expression of the let-7 small regulatory RNA is controlled by ecdysone during metamorphosis in Drosophila melanogaster . Dev Biol. 2002, 244: 170-179. 10.1006/dbio.2002.0594.PubMedView ArticleGoogle Scholar
- Corbitt TS, Hardie J: Juvenile hormone effects on polymorphism in the pea aphid, Acyrthosiphon pisum. Entomol Exp Appl. 1985, 38: 131-135. 10.1007/BF00361923.View ArticleGoogle Scholar
- Le Trionnaire G, Jaubert S, Sabater-Munoz B, Benedetto A, Bonhomme J, Prunier-Leterme N, Martinez-Torres D, Simon JC, Tagu D: Seasonal photoperiodism regulates the expression of cuticular and signalling protein genes in the pea aphid. Ins Biochem Mol Biol. 2007, 37: 1094-1102. 10.1016/j.ibmb.2007.06.008.View ArticleGoogle Scholar
- Hofacker IL, Fontana W, Stadler PF, Bonhoeffer LS, Tacker M, Shuster P: Fast folding and comparison of RNA secondary structures. Monatshefte. 1994, 125: 167-188. 10.1007/BF00818163.View ArticleGoogle Scholar
- Wang X, Zhang J, Li F, Gu J, He T, Zhang X, Li Y: MicroRNA identification based on sequence and structure alignment. Bioinformatics. 2005, 21: 3610-3614. 10.1093/bioinformatics/bti562.PubMedView ArticleGoogle Scholar
- Helvik SA, Snøve O, Saetrom P: Reliable prediction of Drosha processing improves microRNA gene prediction. Bioinformatics. 2007, 23: 142-149. 10.1093/bioinformatics/btl570.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.