Next-generation small RNA sequencing for microRNAs profiling in the honey bee Apis mellifera

Ame-miR-980-3p participates in autophagy-mediated midgut remodelling in Apis mellifera via targeting Atg2B

Autophagy is a process that serves to degrade damaged proteins and organelles, thereby promoting cell homeostasis, differentiation, development and survival. Many miRNAs have been found to have regulatory roles in autophagy. In insects, it has been shown that autophagy is involved in hormone-regulated programmed cell death during metamorphic midgut remodelling. However, whether this is also true during the remodelling of the honey bee midgut is unclear. In the present study, we explored the relationship between autophagy and midgut remodelling and sought to identify miRNAs involved in this physiological process. We found that autophagy occurred during midgut remodelling and that the inhibition of autophagy resulted in midgut dysplasia in prepupae. Differentially expressed miRNAs enriched in the autophagy signalling pathway during midgut remodelling were identified by small RNA-seq. Ame-miR-980-3p, which targets the autophagy-related gene Atg2B, was screened out. Furthermore, abnormal expression of ame-miR-980-3p in the pupal stage led to the thinning of the midgut wall of newly emerged bees (NE). When ame-miR-980-3p expression was inhibited, the intestinal villi of NE bees became significantly shorter and sparse, and the lipid signal in the peritrophic matrix of Pb almost disappeared, indicating that the adult midgut was underdeveloped and the lipid absorption ability was weakened. Taken together, ame-miR-980-3p targeted Atg2B to participate in the regulation of midgut autophagy in the pupae, and the abnormal expression of ame-miR-980-3p would interfere with cell proliferation and death in the process of midgut remodelling, hinder the formation of adult midgut and eventually lead to adult midgut dysplasia and affect the lipid absorption function of the midgut in Apis mellifera.

MicroRNA profiling of royal jelly extracellular vesicles and their potential role in cell viability and reversing cell apoptosis

MiRNAs are small non-coding RNA molecules that play important regulatory roles in diverse biological processes. Royal jelly, a milky-white substance produced by nurse honeybees (Apis mellifera), is the primary food of queen bees and plays a crucial role in their development. However, little is known about the microRNA (miRNAs) content of royal jelly and their potential functions. In this study, we isolated extracellular vesicles from the royal jelly of 36 samples through sequential centrifugation and targeted nanofiltration and performed high-throughput sequencing to identify and quantify the miRNA content of honeybee royal jelly extracellular vesicles (RJEVs). We found a total of 29 known mature miRNAs and 17 novel miRNAs. Through bioinformatic analysis, we identified several potential target genes of the miRNAs present in royal jelly, including those involved in developmental processes and cell differentiation. To investigate the potential roles of RJEVs in cell viability, RJEVs were supplemented to apoptotic porcine kidney fibroblasts induced by ethanol 6% exposure for 30 min. TUNEL assay showed a significant reduction in the apoptosis percentage after RJEV supplementation when compared with the non-supplemented control group. Moreover, the wound healing assay performed on the apoptotic cells showed a rapid healing capacity of RJEV-supplemented cells compared to the control group. We observed a significant reduction in the expression of the miRNA target genes such as FAM131B, ZEB1, COL5A1, TRIB2, YBX3, MAP2, CTNNA1, and ADAMTS9 suggesting that RJEVs may regulate the target gene expression associated with cellular motility and cell viability. Moreover, RJEVs reduced the expression of apoptotic genes (CASP3, TP53, BAX, and BAK), while significantly increasing the expression of anti-apoptotic genes (BCL2 and BCL-XL). Our findings provide the first comprehensive analysis of the miRNA content of RJEVs and suggest a potential role for these vesicles in the regulation of gene expression and cell survival as well as augmenting cell resurrection or anastasis.

MicroRNAs reshape the immunity of insects in response to bacterial infection

The interaction between bacteria and insects can significantly impact a wide range of different areas because bacteria and insects are widely distributed around the globe. The bacterial-insect interactions have the potential to directly affect human health since insects are vectors for disease transmission, and their interactions can also have economic consequences. In addition, they have been linked to high mortality rates in economically important insects, resulting in substantial economic losses. MicroRNAs (miRNAs) are types of non-coding RNAs involved in regulating gene expression post-transcriptionally. The length of miRNAs ranges from 19 to 22 nucleotides. MiRNAs, in addition to their ability to exhibit dynamic expression patterns, have a diverse range of targets. This enables them to govern various physiological activities in insects, like innate immune responses. Increasing evidence suggests that miRNAs have a crucial biological role in bacterial infection by influencing immune responses and other mechanisms for resistance. This review focuses on some of the most recent and exciting discoveries made in recent years, including the correlation between the dysregulation of miRNA expression in the context of bacterial infection and the progression of the infection. Furthermore, it describes how they profoundly impact the immune responses of the host by targeting the Toll, IMD, and JNK signaling pathways. It also emphasizes the biological function of miRNAs in regulating immune responses in insects. Finally, it also discusses current knowledge gaps about the function of miRNAs in insect immunity, in addition to areas that require more research in the future.

Large-Scale Annotation and Evolution Analysis of MiRNA in Insects

Insects are among the most diverse and successful groups of animals and exhibit great morphological diversity and complexity. The innovation of wings and metamorphosis are some examples of the fascinating biological evolution of insects. Most microRNAs (miRNAs) contribute to canalization by conferring robustness to gene networks and thus increase the heritability of important phenotypes. Though previous studies have demonstrated how miRNAs regulate important phenotypes, little is still known about miRNA evolution in insects. Here, we used both small RNA-seq data and homology searching methods to annotate the miRNA repertoires of 152 arthropod species, including 135 insects and 17 noninsect arthropods. We identified 16,212 miRNA genes, and classified them into highly conserved (62), insect-conserved (90), and lineage-specific (354) miRNA families. The phylogenetic relationship of miRNA binary presence/absence dynamics implies that homoplastic loss of conserved miRNA families tends to occur in far-related morphologically simplified taxa, including scale insects (Coccoidea) and twisted-wing insects (Strepsiptera), leading to inconsistent phylogenetic tree reconstruction. The common ancestor of Insecta shares 62 conserved miRNA families, of which five were rapidly gained in the early winged-insects (Pterygota). We also detected extensive miRNA losses in Paraneoptera that are correlated with morphological reduction, and miRNA gains in early Endopterygota around the time holometabolous metamorphosis appeared. This was followed by abundant miRNA gains in Hymenoptera and Lepidoptera. In summary, we provide a comprehensive data set and a detailed evolutionary analysis of miRNAs in insects. These data will be important for future studies on miRNA functions associated with insect morphological innovation and trait biodiversity.

Elucidation of the microRNA Transcriptome in Western Corn Rootworm Reveals Its Dynamic and Evolutionary Complexity

Diabrotica virgifera virgifera (western corn rootworm, WCR) is one of the most destructive agricultural insect pests in North America. It is highly adaptive to environmental stimuli and crop protection technologies. However, little is known about the underlying genetic basis of WCR behavior and adaptation. More specifically, the involvement of small RNAs (sRNAs), especially microRNAs (miRNAs), a class of endogenous small non-coding RNAs that regulate various biological processes, has not been examined, and the datasets of putative sRNA sequences have not previously been generated for WCR. To achieve a comprehensive collection of sRNA transcriptomes in WCR, we constructed, sequenced, and analyzed sRNA libraries from different life stages of WCR and northern corn rootworm (NCR), and identified 101 conserved precursor miRNAs (pre-miRNAs) in WCR and other Arthropoda. We also identified 277 corn rootworm specific pre-miRNAs. Systematic analyses of sRNA populations in WCR revealed that its sRNA transcriptome, which includes PIWI-interacting RNAs (piRNAs) and miRNAs, undergoes a dynamic change throughout insect development. Phylogenetic analysis of miRNA datasets from model species reveals that a large pool of species-specific miRNAs exists in corn rootworm; these are potentially evolutionarily transient. Comparisons of WCR miRNA clusters to other insect species highlight conserved miRNA-regulated processes that are common to insects. Parallel Analysis of RNA Ends (PARE) also uncovered potential miRNA-guided cleavage sites in WCR. Overall, this study provides a new resource for studying the sRNA transcriptome and miRNA-mediated gene regulation in WCR and other Coleopteran insects.

Identification of anti-inflammatory vesicle-like nanoparticles in honey

Honey has been used as a nutrient, an ointment, and a medicine worldwide for many centuries. Modern research has demonstrated that honey has many medicinal properties, reflected in its anti-microbial, anti-oxidant, and anti-inflammatory bioactivities. Honey is composed of sugars, water and a myriad of minor components, including minerals, vitamins, proteins and polyphenols. Here, we report a new bioactive component‒vesicle-like nanoparticles‒in honey (H-VLNs). These H-VLNs are membrane-bound nano-scale particles that contain lipids, proteins and small-sized RNAs. The presence of plant-originated plasma transmembrane proteins and plasma membrane-associated proteins suggests the potential vesicle-like nature of these particles. H-VLNs impede the formation and activation of the nucleotide-binding domain and leucine-rich repeat related (NLR) family, pyrin domain containing 3 (NLRP3) inflammasome, which is a crucial inflammatory signalling platform in the innate immune system. Intraperitoneal administration of H-VLNs in mice alleviates inflammation and liver damage in the experimentally induced acute liver injury. miR-4057 in H-VLNs was identified in inhibiting NLRP3 inflammasome activation. Together, our studies have identified anti-inflammatory VLNs as a new bioactive agent in honey.

Comparative Identification of MicroRNAs in Apis cerana cerana Workers' Midguts in Responseto Nosema ceranae Invasion

Here, the expression profiles and differentially expressed miRNAs (DEmiRNAs) in the midguts of Apis cerana cerana workers at 7 d and 10 d post-inoculation (dpi) with N. ceranae were investigated via small RNA sequencing and bioinformatics. Five hundred and twenty nine (529) known miRNAs and 25 novel miRNAs were identified in this study, and the expression of 16 predicted miRNAs was confirmed by Stem-loop RT-PCR. A total of 14 DEmiRNAs were detected in the midgut at 7 dpi, including eight up-regulated and six down-regulated miRNAs, while 12 DEmiRNAs were observed in the midgut at 10 dpi, including nine up-regulated and three down-regulated ones. Additionally, five DEmiRNAs were shared, while nine and seven DEmiRNAs were specifically expressed in midguts at 7 dpi and 10 dpi. Gene ontology analysis suggested some DEmiRNAs and corresponding target mRNAs were involved in various functions including immune system processes and response to stimulus. KEGG pathway analysis shed light on the potential functions of some DEmiRNAs in regulating target mRNAs engaged in material and energy metabolisms, cellular immunity and the humoral immune system. Further investigation demonstrated a complex regulation network between DEmiRNAs and their target mRNAs, with miR-598-y, miR-252-y, miR-92-x and miR-3654-y at the center. Our results can facilitate future exploration of the regulatory roles of miRNAs in host responses to N. ceranae, and provide potential candidates for further investigation of the molecular mechanisms underlying eastern honeybee-microsporidian interactions.

The evolution of posttranscriptional regulation

"DNA makes RNA makes protein." After transcription, mRNAs undergo a series of intertwining processes to be finally translated into functional proteins. The "posttranscriptional" regulation (PTR) provides cells an extended option to fine-tune their proteomes. To meet the demands of complex organism development and the appropriate response to environmental stimuli, every step in these processes needs to be finely regulated. Moreover, changes in these regulatory processes are important driving forces underlying the evolution of phenotypic differences across different species. The major PTR mechanisms discussed in this review include the regulation of splicing, polyadenylation, decay, and translation. For alternative splicing and polyadenylation, we mainly discuss their evolutionary dynamics and the genetic changes underlying the regulatory differences in cis-elements versus trans-factors. For mRNA decay and translation, which, together with transcription, determine the cellular RNA or protein abundance, we focus our discussion on how their divergence coordinates with transcriptional changes to shape the evolution of gene expression. Then to highlight the importance of PTR in the evolution of higher complexity, we focus on their roles in two major phenomena during eukaryotic evolution: the evolution of multicellularity and the division of labor between different cell types and tissues; and the emergence of diverse, often highly specialized individual phenotypes, especially those concerning behavior in eusocial insects. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution Translation > Translation Regulation RNA Processing > Splicing Regulation/Alternative Splicing.

Small RNA profiling of low biomass samples: identification and removal of contaminants

Background

Sequencing-based analyses of low-biomass samples are known to be prone to misinterpretation due to the potential presence of contaminating molecules derived from laboratory reagents and environments. DNA contamination has been previously reported, yet contamination with RNA is usually considered to be very unlikely due to its inherent instability. Small RNAs (sRNAs) identified in tissues and bodily fluids, such as blood plasma, have implications for physiology and pathology, and therefore the potential to act as disease biomarkers. Thus, the possibility for RNA contaminants demands careful evaluation.

Results

Herein, we report on the presence of small RNA (sRNA) contaminants in widely used microRNA extraction kits and propose an approach for their depletion. We sequenced sRNAs extracted from human plasma samples and detected important levels of non-human (exogenous) sequences whose source could be traced to the microRNA extraction columns through a careful qPCR-based analysis of several laboratory reagents. Furthermore, we also detected the presence of artefactual sequences related to these contaminants in a range of published datasets, thereby arguing in particular for a re-evaluation of reports suggesting the presence of exogenous RNAs of microbial and dietary origin in blood plasma. To avoid artefacts in future experiments, we also devise several protocols for the removal of contaminant RNAs, define minimal amounts of starting material for artefact-free analyses, and confirm the reduction of contaminant levels for identification of bona fide sequences using ‘ultra-clean’ extraction kits.

Conclusion

This is the first report on the presence of RNA molecules as contaminants in RNA extraction kits. The described protocols should be applied in the future to avoid confounding sRNA studies.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0522-7) contains supplementary material, which is available to authorized users.

The microRNA ame-miR-279a regulates sucrose responsiveness of forager honey bees (Apis mellifera)

Increasing evidence demonstrates that microRNAs (miRNA) play an important role in the regulation of animal behaviours. Honey bees (Apis mellifera) are eusocial insects, with honey bee workers displaying age-dependent behavioural maturation. Many different miRNAs have been implicated in the change of behaviours in honey bees and ame-miR-279a was previously shown to be more highly expressed in nurse bee heads than in those of foragers. However, it was not clear whether this difference in expression was associated with age or task performance. Here we show that ame-miR-279a shows significantly higher expression in the brains of nurse bees relative to forager bees regardless of their ages, and that ame-miR-279a is primarily localized in the Kenyon cells of the mushroom body in both foragers and nurses. Overexpression of ame-miR-279a attenuates the sucrose responsiveness of foragers, while its absence enhances their sucrose responsiveness. Lastly, we determined that ame-miR-279a directly target the mRNA of Mblk-1. These findings suggest that ame-miR-279a plays important roles in regulating honey bee division of labour.

Molecular Evolution of Insect Sociality: An Eco-Evo-Devo Perspective

The evolution of eusociality is a perennial issue in evolutionary biology, and genomic advances have fueled steadily growing interest in the genetic changes underlying social evolution. Along with a recent flurry of research on comparative and evolutionary genomics in different eusocial insect groups (bees, ants, wasps, and termites), several mechanistic explanations have emerged to describe the molecular evolution of eusociality from solitary behavior. These include solitary physiological ground plans, genetic toolkits of deeply conserved genes, evolutionary changes in protein-coding genes, cis regulation, and the structure of gene networks, epigenetics, and novel genes. Despite this proliferation of ideas, there has been little synthesis, even though these ideas are not mutually exclusive and may in fact be complementary. We review available data on molecular evolution of insect sociality and highlight key biotic and abiotic factors influencing social insect genomes. We then suggest both phylogenetic and ecological evolutionary developmental biology (eco-evo-devo) perspectives for a more synthetic view of molecular evolution in insect societies.

Identification and characterization of microRNAs expressed in antennae of Holotrichia parallela Motschulsky and their possible roles in olfactory regulation

MicroRNAs (miRNAs) are small noncoding RNAs that play posttranscriptional, regulatory roles in various biological processes. However, there has been limited investigation into the potential function of miRNAs in olfaction. The coleopteran Holotrichia parallela is an economically important pest, and miRNAs have been identified in only one coleopteran (Tribolium castaneum). Therefore, this study was conducted to identify miRNAs expressed in the antennae of H. parallela and obtain insights into their possible roles in olfaction. By combining deep sequencing and miRDeep2 software, a total of 99 miRNAs, including 76 conserved miRNAs and 23 novel miRNAs, were identified from H. parallela antennae. The 76 conserved miRNAs belong to 63 families and the other 23 may be species specific or tissue specific. The identified miRNAs have many conserved features of miRNAs. Evaluation of the conservation of the identified miRNA families across different species revealed that most of the families are insect specific. The prediction and annotation of targets suggested that 13 of the identified miRNAs participate in olfactory regulation. Gender differences in antennal expression of nine of the olfactory-related miRNAs were confirmed by quantitative real-time PCR.

Shrimp miRNAs regulate innate immune response against white spot syndrome virus infection

MicroRNAs are short noncoding RNAs of RNA interference pathways that regulate gene expression through partial complementary base-pairing to target mRNAs. In this study, miRNAs that are expressed in white spot syndrome virus (WSSV)-infected Penaeus monodon, were identified using next generation sequencing. Forty-six miRNA homologs were identified from WSSV-infected shrimp hemocyte. Stem-loop real-time RT-PCR analysis showed that 11 out of 16 selected miRNAs were differentially expressed upon WSSV infection. Of those, pmo-miR-315 and pmo-miR-750 were highly responsive miRNAs. miRNA target prediction revealed that the miRNAs were targeted at 5'UTR, ORF, and 3'UTR of several immune-related genes such as genes encoding antimicrobial peptides, signaling transduction proteins, heat shock proteins, oxidative stress proteins, proteinases or proteinase inhibitors, proteins in blood clotting system, apoptosis-related proteins, proteins in prophenoloxidase system, pattern recognition proteins and other immune molecules. The highly conserved miRNA homolog, pmo-bantam, was characterized for its function in shrimp. The pmo-bantam was predicted to target the 3'UTR of Kunitz-type serine protease inhibitor (KuSPI). Binding of pmo-bantam to the target sequence of KuSPI gene was analyzed by luciferase reporter assay. Correlation of pmo-bantam and KuSPI expression was observed in lymphoid organ of WSSV-infected shrimp. These results implied that miRNAs might play roles as immune gene regulators in shrimp antiviral response.

MicroRNAs in Honey Bee Caste Determination

The cellular mechanisms employed by some organisms to produce contrasting morphological and reproductive phenotypes from the same genome remains one of the key unresolved issues in biology. Honeybees (Apis mellifera) use differential feeding and a haplodiploid sex determination system to generate three distinct organismal outcomes from the same genome. Here we investigate the honeybee female and male caste-specific microRNA and transcriptomic molecular signatures during a critical time of larval development. Both previously undetected and novel miRNAs have been discovered, expanding the inventory of these genomic regulators in invertebrates. We show significant differences in the microRNA and transcriptional profiles of diploid females relative to haploid drone males as well as between reproductively distinct females (queens and workers). Queens and drones show gene enrichment in physio-metabolic pathways, whereas workers show enrichment in processes associated with neuronal development, cell signalling and caste biased structural differences. Interestingly, predicted miRNA targets are primarily associated with non-physio-metabolic genes, especially neuronal targets, suggesting a mechanistic disjunction from DNA methylation that regulates physio-metabolic processes. Accordingly, miRNA targets are under-represented in methylated genes. Our data show how a common set of genetic elements are differentially harnessed by an organism, which may provide the remarkable level of developmental flexibility required.

Identification and characterization of microRNAs in Eucheuma denticulatum by high-throughput sequencing and bioinformatics analysis

Eucheuma denticulatum, an economically and industrially important red alga, is a valuable marine resource. Although microRNAs (miRNAs) play an essential role in gene post-transcriptional regulation, no research has been conducted to identify and characterize miRNAs in E. denticulatum. In this study, we identified 134 miRNAs (133 conserved miRNAs and one novel miRNA) from 2,997,135 small-RNA reads by high-throughput sequencing combined with bioinformatics analysis. BLAST searching against miRBase uncovered 126 potential miRNA families. A conservation and diversity analysis of predicted miRNA families in different plant species was performed by comparative alignment and homology searching. A total of 4 and 13 randomly selected miRNAs were respectively validated by northern blotting and stem-loop reverse transcription PCR, thereby demonstrating the reliability of the miRNA sequencing data. Altogether, 871 potential target genes were predicted using psRobot and TargetFinder. Target genes classification and enrichment were conducted based on Gene Ontology analysis. The functions of target gene products and associated metabolic pathways were predicted by Kyoto Encyclopedia of Genes and Genomes pathway analysis. A Cytoscape network was constructed to explore the interrelationships of miRNAs, miRNA-target genes and target genes. A large number of miRNAs with diverse target genes will play important roles for further understanding some essential biological processes in E. denticulatum. The uncovered information can serve as an important reference for the protection and utilization of this unique red alga in the future.

Characterisation of the small RNAs in the biomedically important green-bottle blowfly Lucilia sericata

BACKGROUND

The green bottle fly maggot, Lucilia sericata, is a species with importance in medicine, agriculture and forensics. Improved understanding of this species' biology is of great potential benefit to many research communities. MicroRNAs (miRNA) are a short non-protein coding regulatory RNA, which directly regulate a host of protein coding genes at the translational level. They have been shown to have developmental and tissue specific distributions where they impact directly on gene regulation. In order to improve understanding of the biology of L. sericata maggots we have performed small RNA-sequencing of their secretions and tissue at different developmental stages.

RESULTS

We have successfully isolated RNA from the secretions of L. sericata maggots. Illumina small RNA-sequencing of these secretions and the three tissues (crop, salivary gland, gut) revealed that the most common small RNA fragments were derived from ribosomal RNA and transfer RNAs of both insect and bacterial origins. These RNA fragments were highly specific, with the most common tRNAs, such as GlyGCC, predominantly represented by reads derived from the 5' end of the mature maggot tRNA. Each library also had a unique profile of miRNAs with a high abundance of miR-10-5p in the maggot secretions and gut and miR-8 in the food storage organ the crop and salivary glands. The pattern of small RNAs in the bioactive maggot secretions suggests they originate from a combination of saliva, foregut and hindgut tissues. Droplet digital RT-PCR validation of the RNA-sequencing data shows that not only are there differences in the tissue profiles for miRNAs and small RNA fragments but that these are also modulated through developmental stages of the insect.

CONCLUSIONS

We have identified the small-RNAome of the medicinal maggots L. sericata and shown that there are distinct subsets of miRNAs expressed in specific tissues that also alter during the development of the insect. Furthermore there are very specific RNA fragments derived from other non-coding RNAs present in tissues and in the secretions. This new knowledge has applicability in diverse research fields including wound healing, agriculture and forensics.

Developmental regulation of ecdysone receptor (EcR) and EcR-controlled gene expression during pharate-adult development of honeybees (Apis mellifera)

Major developmental transitions in multicellular organisms are driven by steroid hormones. In insects, these, together with juvenile hormone (JH), control development, metamorphosis, reproduction and aging, and are also suggested to play an important role in caste differentiation of social insects. Here, we aimed to determine how EcR transcription and ecdysteroid titers are related during honeybee postembryonic development and what may actually be the role of EcR in caste development of this social insect. In addition, we expected that knocking-down EcR gene expression would give us information on the participation of the respective protein in regulating downstream targets of EcR. We found that in Apis mellifera females, EcR-A is the predominantly expressed variant in postembryonic development, while EcR-B transcript levels are higher in embryos, indicating an early developmental switch in EcR function. During larval and pupal stages, EcR-B expression levels are very low, while EcR-A transcripts are more variable and abundant in workers compared to queens. Strikingly, these transcript levels are opposite to the ecdysteroid titer profile. 20-hydroxyecdysone (20E) application experiments revealed that low 20E levels induce EcR expression during development, whereas high ecdysteroid titers seem to be repressive. By means of RNAi-mediated knockdown (KD) of both EcR transcript variants we detected the differential expression of 234 poly-A⁺ transcripts encoding genes such as CYPs, MRJPs and certain hormone response genes (Kr-h1 and ftz-f1). EcR-KD also promoted the differential expression of 70 miRNAs, including highly conserved ones (e.g., miR-133 and miR-375), as well honeybee-specific ones (e.g., miR-3745 and miR-3761). Our results put in evidence a broad spectrum of EcR-controlled gene expression during postembryonic development of honeybees, revealing new facets of EcR biology in this social insect.

A critical appraisal of the use of microRNA data in phylogenetics

Recent progress in resolving the tree of life continues to expose relationships that resist resolution, which drives the search for novel sources of information to solve these difficult phylogenetic problems. A recent example, the presence and absence of microRNA families, has been vigorously promoted as an ideal source of phylogenetic data and has been applied to several perennial phylogenetic problems. The utility of such data for phylogenetic inference hinges critically both on developing stochastic models that provide a reasonable description of the process that give rise to these data, and also on the careful validation of those models in real inference scenarios. Remarkably, however, the statistical behavior and phylogenetic utility of microRNA data have not yet been rigorously characterized. Here we explore the behavior and performance of microRNA presence/absence data under a variety of evolutionary models and reexamine datasets from several previous studies. We find that highly heterogeneous rates of microRNA gain and loss, pervasive secondary loss, and sampling error collectively render microRNA-based inference of phylogeny difficult. Moreover, our reanalyses fundamentally alter the conclusions for four of the five studies that we reexamined. Our results indicate that the capacity of miRNA data to resolve the tree of life has been overstated, and we urge caution in their application and interpretation.

Identification and characterization of microRNAs in small brown planthopper (Laodephax striatellus) by next-generation sequencing

MicroRNAs (miRNAs) are endogenous non-coding small RNAs that regulate gene expression at the post-transcriptional level and are thought to play critical roles in many metabolic activities in eukaryotes. The small brown planthopper (Laodephax striatellus Fallén), one of the most destructive agricultural pests, causes great damage to crops including rice, wheat, and maize. However, information about the genome of L. striatellus is limited. In this study, a small RNA library was constructed from a mixed L. striatellus population and sequenced by Solexa sequencing technology. A total of 501 mature miRNAs were identified, including 227 conserved and 274 novel miRNAs belonging to 125 and 250 families, respectively. Sixty-nine conserved miRNAs that are included in 38 families are predicted to have an RNA secondary structure typically found in miRNAs. Many miRNAs were validated by stem-loop RT-PCR. Comparison with the miRNAs in 84 animal species from miRBase showed that the conserved miRNA families we identified are highly conserved in the Arthropoda phylum. Furthermore, miRanda predicted 2701 target genes for 378 miRNAs, which could be categorized into 52 functional groups annotated by gene ontology. The function of miRNA target genes was found to be very similar between conserved and novel miRNAs. This study of miRNAs in L. striatellus will provide new information and enhance the understanding of the role of miRNAs in the regulation of L. striatellus metabolism and development.

A virulent strain of deformed wing virus (DWV) of honeybees (Apis mellifera) prevails after Varroa destructor-mediated, or in vitro, transmission

The globally distributed ectoparasite Varroa destructor is a vector for viral pathogens of the Western honeybee (Apis mellifera), in particular the Iflavirus Deformed Wing Virus (DWV). In the absence of Varroa low levels DWV occur, generally causing asymptomatic infections. Conversely, Varroa-infested colonies show markedly elevated virus levels, increased overwintering colony losses, with impairment of pupal development and symptomatic workers. To determine whether changes in the virus population were due Varroa amplifying and introducing virulent virus strains and/or suppressing the host immune responses, we exposed Varroa-naïve larvae to oral and Varroa-transmitted DWV. We monitored virus levels and diversity in developing pupae and associated Varroa, the resulting RNAi response and transcriptome changes in the host. Exposed pupae were stratified by Varroa association (presence/absence) and virus levels (low/high) into three groups. Varroa-free pupae all exhibited low levels of a highly diverse DWV population, with those exposed per os (group NV) exhibiting changes in the population composition. Varroa-associated pupae exhibited either low levels of a diverse DWV population (group VL) or high levels of a near-clonal virulent variant of DWV (group VH). These groups and unexposed controls (C) could be also discriminated by principal component analysis of the transcriptome changes observed, which included several genes involved in development and the immune response. All Varroa tested contained a diverse replicating DWV population implying the virulent variant present in group VH, and predominating in RNA-seq analysis of temporally and geographically separate Varroa-infested colonies, was selected upon transmission from Varroa, a conclusion supported by direct injection of pupae in vitro with mixed virus populations. Identification of a virulent variant of DWV, the role of Varroa in its transmission and the resulting host transcriptome changes furthers our understanding of this important viral pathogen of honeybees.

Honeybees are the most important managed pollinating insect, contributing billions of dollars to annual global agricultural production. Over the last century a parasitic mite, Varroa, has spread worldwide, with significant impacts on honeybee colony health as a consequence of its transmission of a cocktail of viruses while feeding on honeybee ‘blood’. The most important virus for colony health is deformed wing virus (DWV), high levels of which cause developmental deformities and premature ageing resulting in high overwintering colony losses. In experiments on individual Varroa-exposed pupae we demonstrate that a single type of virulent DWV is amplified 1,000–10,000 times in the recipient pupae, despite the mite containing a high diversity of replicating DWV strains. We could recapitulate this by direct injection of pupae with mixed virus populations, showing the virulent strain is advantaged by the route of transmission. In parallel, we detected changes in the immune response and developmental gene expression of the honeybee and propose that these contribute to the characteristic pathogenesis of DWV. Identification of a virulent strain of DWV has implications for therapeutic or prophylactic interventions to improve honeybee colony health, as well as contributing to our understanding of the biology of this important honeybee viral pathogen.

An analysis of the small RNA transcriptome of four developmental stages of the citrus red mite (Panonychus citri)

The citrus red mite (Panonychus citri) can feed on more than 112 plant species around the world. Endogenous small RNAs (sRNAs) have proved to be important components of gene regulation in many eukaryotes. Recently, many sRNAs have been shown to be involved in various biological processes, such as development in many animals, including insects; however, to date, no sRNAs have been reported in the citrus red mite. Using Illumina sequencing, several categories of sRNAs were identified, including 594 known microRNAs (miRNAs) grouped into 206 families and 31 novel miRNAs in the four developmental stages of citrus red mite. In addition, according to bioinformatics analysis and S-Poly(T) miRNA assays, the expression level of many miRNAs varied among the developmental stages. Furthermore, the prediction of miRNAs target genes and their functional annotation indicated that miRNAs are involved in the regulation of multiple pathways in the citrus red mite. As the first report of the sRNA world in citrus red mite, the present study furthers our understanding of the roles played by sRNAs in the development of citrus red mite and the data may help to develop methods of controlling the pests in the field.

Odor memories regulate olfactory receptor expression in the sensory periphery

Odor learning induces structural and functional modifications throughout the olfactory system, but it is currently unknown whether this plasticity extends to the olfactory receptors (Or) in the sensory periphery. Here, we demonstrate that odor learning induces plasticity in olfactory receptor expression in the honeybee, Apis mellifera. Using quantitative RT-PCR analysis, we show that six putative floral scent receptors were differentially expressed in the bee antennae depending on the scent environment that the bees experienced. Or151, which we characterized using an in vitro cell expression system as a broadly tuned receptor binding floral odorants such as linalool, and Or11, the specific receptor for the queen pheromone 9-oxo-decenoic acid, were significantly down-regulated after honeybees were conditioned with the respective odorants in an olfactory learning paradigm. Electroantennogram recordings showed that the neural response of the antenna was similarly reduced after odor learning. Long-term odor memory was essential for inducing these changes, suggesting that the molecular mechanisms involved in olfactory memory also regulate olfactory receptor expression. Our study demonstrates for the first time that olfactory receptor expression is experience-dependent and modulated by scent conditioning, providing novel insight into how molecular regulation at the periphery contributes to plasticity in the olfactory system.

Recipe for a busy bee: microRNAs in Honey Bee caste determination

Social caste determination in the honey bee is assumed to be determined by the dietary status of the young larvae and translated into physiological and epigenetic changes through nutrient-sensing pathways. We have employed Illumina/Solexa sequencing to examine the small RNA content in the bee larval food, and show that worker jelly is enriched in miRNA complexity and abundance relative to royal jelly. The miRNA levels in worker jelly were 7-215 fold higher than in royal jelly, and both jellies showed dynamic changes in miRNA content during the 4(th) to 6(th) day of larval development. Adding specific miRNAs to royal jelly elicited significant changes in queen larval mRNA expression and morphological characters of the emerging adult queen bee. We propose that miRNAs in the nurse bee secretions constitute an additional element in the regulatory control of caste determination in the honey bee.

Sequence and expression characteristics of long noncoding RNAs in honey bee caste development--potential novel regulators for transgressive ovary size

Division of labor in social insect colonies relies on a strong reproductive bias that favors queens. Although the ecological and evolutionary success attained through caste systems is well sketched out in terms of ultimate causes, the molecular and cellular underpinnings driving the development of caste phenotypes are still far from understood. Recent genomics approaches on honey bee developmental biology revealed a set of genes that are differentially expressed genes in larval ovaries and associated with transgressive ovary size in queens and massive cell death in workers. Amongst these, two contigs called special attention, both being over 200 bp in size and lacking apparent coding potential. Herein, we obtained their full cDNA sequences. These and their secondary structure characteristics placed in evidence that they are bona fide long noncoding RNAs (lncRNA) differentially expressed in larval ovaries, thus named lncov1 and lncov2. Genomically, both map within a previously identified QTL on chromosome 11, associated with transgressive ovary size in honey bee workers. As lncov1 was over-expressed in worker ovaries we focused on this gene. Real-time qPCR analysis on larval worker ovaries evidenced an expression peak coinciding with the onset of autophagic cell death. Cellular localization analysis through fluorescence in situ hybridization revealed perinuclear spots resembling omega speckles known to regulate trafficking of RNA-binding proteins. With only four lncRNAs known so far in honey bees, two expressed in the ovaries, these findings open a novel perspective on regulatory factors acting in the fine tuning of developmental processes underlying phenotypic plasticity related to social life histories.

Bacterial infection activates the immune system response and dysregulates microRNA expression in honey bees

In insects, a rapid and massive synthesis of antimicrobial peptides (AMPs) is activated through signaling pathways (Toll and Imd) to combat invading microbial pathogens. However, it is still unclear whether different types of bacteria provoke specific responses. Immune response mechanisms and the activation of specific genes were investigated by challenging Apis mellifera workers with the Gram-negative bacterium Serratia marcescens or the Gram-positive bacterium Micrococcus luteus. The immune system responded by activating most genes of the Toll and Imd pathways, particularly AMP genes. However, genes specifically regulated by M. luteus or S. marcescens were not detected, suggesting an interaction between the signaling pathways that lead to immune effectors synthesis. Despite this finding, kappaB motifs in the 5'-UTRs of selected genes suggest a pathway-specific control of AMP and transferrin-1 gene expression. Regulation by miRNAs was also investigated and revealed a number of candidates for the post-transcriptional regulation of immune genes in bees.

Characterization and comparative profiling of miRNAs in invasive Bemisia tabaci (Gennadius) B and Q

BACKGROUND

MicroRNAs (miRNAs) are small, conserved, non-coding RNAs that post-transcriptionally regulate gene expression. Bemisia tabaci (Gennadius) B and Q are two invasive and dominant whiteflies, and B. tabaci Q has been displacing B in China. Differences in biological traits (fecundity, host range, resistance to insecticides, etc.) as affected by miRNAs might be involved in the displacement. In this study, we performed high-throughput sequencing to identify miRNAs in B. tabaci B and Q.

RESULTS

We identified 170 conserved miRNAs and 15 novel candidates, and found significant differences in the expression of miRNAs between B. tabaci B and Q.

CONCLUSION

Expression levels of miRNAs differ in B. tabaci B vs. Q. Additional research is needed to determine whether these differences are related to differences in the biology of B. tabaci B and Q, and whether these differences help explain why B. tabaci Q is displacing B in China.

Computational identification and characterization of putative miRNAs in Heliothis virescens

Heliothis virescens, a polyphagous pest, is one of the most destructive pests of many crops and vegetables. Various insecticides and pesticides are used by agriculturalists to stop the growth and development of this pest. RNA interference is a new area for the management of pests/insects by inhibiting the growth related RNAs. This involves the miRNAs identification and its characterization. In the present study, computational approach is applied to predict putative miRNA candidates along with their possible target(s) in the Heliothis virescens. A total of 63,662 ESTs were downloaded from dbEST database and processed, trimmed and masked through EGassembler. The H. virescens contigs database obtained after assembly was now used to find the putative miRNA candidates by performing a local BLAST with the miRNAs of insects retrieved from miRBase. We have predicted putative miRNA candidates by homology search against all the reported insect miRNAs. These putative miRNAs candidates were further validated and filtered by different features. In addition, we have also attempted to predict the putative targets of these filtered miRNAs, by making use of 3' untranslated regions of mRNAs from B. mori. These miRNAs and their targets in H. virescens will help in improved understanding of molecular mechanisms of miRNA and development of novel and more precise techniques for better understanding some post transcriptional gene silencing.

Structure, evolution and function of the bi-directionally transcribed iab-4/iab-8 microRNA locus in arthropods

In Drosophila melanogaster, the iab-4/iab-8 locus encodes bi-directionally transcribed microRNAs that regulate the function of flanking Hox transcription factors. We show that bi-directional transcription, temporal and spatial expression patterns and Hox regulatory function of the iab-4/iab-8 locus are conserved between fly and the beetle Tribolium castaneum. Computational predictions suggest iab-4 and iab-8 microRNAs can target common sites, and cell-culture assays confirm that iab-4 and iab-8 function overlaps on Hox target sites in both fly and beetle. However, we observe key differences in the way Hox genes are targeted. For instance, abd-A transcripts are targeted only by iab-8 in Drosophila, whereas both iab-4 and iab-8 bind to Tribolium abd-A. Our evolutionary and functional characterization of a bi-directionally transcribed microRNA establishes the iab-4/iab-8 system as a model for understanding how multiple products from sense and antisense microRNAs target common sites.

Insect microRNAs: biogenesis, expression profiling and biological functions

MicroRNAs (miRNA) are a class of endogenous regulatory RNA molecules 21-24 nucleotides in length that modulate gene expression at the post-transcriptional level via base pairing to target sites within messenger RNAs (mRNA). Typically, the miRNA "seed sequence" (nucleotides 2-8 at the 5' end) binds complementary seed match sites within the 3' untranslated region of mRNAs, resulting in either translational inhibition or mRNA degradation. MicroRNAs were first discovered in Caenorhabditis elegans and were shown to be involved in the timed regulation of developmental events. Since their discovery in the 1990s, thousands of potential miRNAs have since been identified in various organisms through small RNA cloning methods and/or computational prediction, and have been shown to play functionally important roles of gene regulation in invertebrates, vertebrates, plants, fungi and viruses. Numerous functions of miRNAs identified in Drosophila melanogaster have demonstrated a great significance of these regulatory molecules. However, elucidation of miRNA roles in non-drosophilid insects presents a challenging and important task.

The gene vitellogenin affects microRNA regulation in honey bee (Apis mellifera) fat body and brain

In honey bees, Vitellogenin (Vg) is hypothesized to be a major factor affecting hormone signaling, food-related behavior, immunity, stress resistance and lifespan. Likewise microRNAs play important roles in posttranscriptional gene regulation and affect many biological processes. The action of microRNAs and Vg are known to intersect in the context of reproduction; however, the role of these associations on social behavior is unknown. The phenotypic effects of Vg knockdown are best established and studied in the forager stage of workers. Thus, we exploited the well-established RNA interference (RNAi) protocol for Vg knockdown to investigate its downstream effects on microRNA population in honey bee foragers' brain and fat body tissue. To identify microRNAs that are differentially expressed between tissues in control and knockdown foragers, we used µParaflo® microfluidic oligonucleotide microRNA microarrays. Our results show 76 and 74 microRNAs were expressed in the brain of control and knockdown foragers whereas 66 and 69 microRNAs were expressed in the fat body of control and knockdown foragers respectively. Target prediction identified potential seed matches for a differentially expressed subset of microRNAs affected by Vg knockdown. These candidate genes are involved in a broad range of biological processes including insulin signaling, juvenile hormone (JH) and ecdysteroid signaling previously shown to affect foraging behavior. Thus, here we demonstrate a causal link between the Vg knockdown forager phenotype and variation in the abundance of microRNAs in different tissues with possible consequences for regulation of foraging behavior.

Differences in microRNAs and their expressions between foraging and dancing honey bees, Apis mellifera L

Many studies have established that microRNAs (miRNAs) regulate gene expression in various biological processes in mammals and insects including honey bees. Dancing behavior is a form of communication unique to honey bees. However, it remains unclear which miRNAs regulate the dancing behavior in honey bees, and how. In the present study, total small RNAs (sRNAs) in Apis mellifera foragers and dancers were extracted and analyzed by a Solexa Sequencer to determine differentially expressed miRNAs. A small percentage (12.62%) of the unique sRNAs (the number of sequence types) were shared between foragers and dancers, but their expression accounted for 92.92% of the total sRNAs (the number of all sequence reads), and the length of them centered around 22nt. Out of 58 previously identified miRNAs, 54 were present in both foragers and dancers and most of them were down-regulated in dancers. The fold-changes of ame-miR-34, ame-miR-210, ame-miR-278 and ame-miR-282 were higher than 2. 86 and 104 novel miRNAs were detected in foragers and dancers, respectively. Furthermore, two known miRNAs (ame-miR-278 and ame-miR-282) were confirmed, by qPCR, to have lower expressions in dancers. The target genes of ame-miR-278 and ame-miR-282 were associated with kinase, neural function, synaptotagmin and energy. These results indicate that miRNAs are substantially different between the foraging and dancing stages, and suggest that miRNAs might play important roles in regulating dancing behaviors in honey bees.

Epigenetic modification of gene expression in honey bees by heterospecific gland secretions

Background

In the honey bee (Apis mellifera), queen and workers have different behavior and reproductive capacity despite possessing the same genome. The primary substance that leads to this differentiation is royal jelly (RJ), which contains a range of proteins, amino acids, vitamins and nucleic acids. MicroRNA (miRNA) has been found to play an important role in regulating the expression of protein-coding genes and cell biology. In this study, we characterized the miRNAs in RJ from two honey bee sister species and determined their possible effect on transcriptome in one species.

Methodology/Principal Findings

We sequenced the miRNAs in RJ either from A. mellifera (RJM) or A. cerana (RJC). We then determined the global transcriptomes of adult A. mellifera developed from larvae fed either with RJM (mRJM) or RJC (mRJC). Finally we analyzed the target genes of those miRNA that are species specific or differentially expressed in the two honey bee species. We show that there were differences in miRNA between RJM and RJC, and that transcriptomes of adult A. mellifera were affected by the two types of RJ. A high proportion (23.3%) of the affected genes were target genes of differential miRNAs.

Conclusion

We show for the first time that there are differences in miRNAs in RJ between A. mellifera and A. cerana. Further, the differences in transcriptomes of bees reared from these two RJs might be related to miRNA differences of the two species. This study provides the first evidence that heterospecific royal jelly can modify gene expression in honey bees through an epigenetic mechanism.

Next-generation small RNA sequencing for microRNAs profiling in Apis mellifera: comparison between nurses and foragers

MicroRNAs (miRNAs) are endogenous small non-coding RNAs regulating gene expression in animals and plants. To find some differentially expressed miRNAs that may be associated with age-dependent behavioural changes in honey bees (Apis mellifera), we applied next-generation high-throughput sequencing technology to detect small RNAs in nurses and foragers. Our results showed that both nurses and foragers had a complicated small RNA population, and the length of small RNAs varied, 22 nucleotides being the predominant length. Combining deep sequencing and bioinformatic analysis, we discovered that nine known miRNAs were significantly different between nurses and foragers (P < 0.01; absolute value of fold-change ≥ 1). Some of their target genes were related to neural function. Moreover, 67 novel miRNAs were identified in nurses and foragers. Ame-miR-31a and ame-miR-13b were further validated using quantitative reverse-transcription PCR assays. The present study provides new information on the miRNA abundance of honey bees, and enhances our understanding of miRNA function in the regulation of honey bee development.

Identification and developmental profiling of conserved and novel microRNAs in Manduca sexta

MicroRNAs (miRNAs) are a group of small RNAs involved in translation inhibition or mRNA degradation. Due to its large size, Manduca sexta has long been used as a model to study insect physiology and biochemistry. While transcriptome studies have greatly enriched our knowledge on M. sexta structural genes, little is known about posttranscriptional regulation by miRNAs in this lepidopteran species. We constructed four small RNA libraries from embryos, 4th instar feeding larvae, pupae, and adults, obtained 21 million reads of 18-31 nucleotides by Illumina sequencing, and found 163 conserved and 13 novel miRNAs. By searching the M. sexta genome assembly, we identified precursors of 82 conserved miRNAs, 76 of which had mapped reads in one or more of these libraries. After normalization, we compared numbers of miRNA and miRNA-star reads in these libraries and observed abundance changes during development. Interestingly, mse-miR-281-star, mse-miR-31-star, mse-miR-965-star, mse-miR-9a-star, mse-miR-9b-star, mse-miR-2a-star, mse-miR-92b-star and mse-miR-279c-star are either more abundant or maintained at similar levels compared to respective mature miRNA strand. Expression profiling of the first set of miRNAs provided insights to their possible involvement in developmental regulation. This study will aid in the annotation of miRNA genes in the genome.

Behavioral plasticity in honey bees is associated with differences in brain microRNA transcriptome

Small, non-coding microRNAs (miRNAs) have been implicated in many biological processes, including the development of the nervous system. However, the roles of miRNAs in natural behavioral and neuronal plasticity are not well understood. To help address this we characterized the microRNA transcriptome in the adult worker honey bee head and investigated whether changes in microRNA expression levels in the brain are associated with division of labor among honey bees, a well-established model for socially regulated behavior. We determined that several miRNAs were downregulated in bees that specialize on brood care (nurses) relative to foragers. Additional experiments showed that this downregulation is dependent upon social context; it only occurred when nurse bees were in colonies that also contained foragers. Analyses of conservation patterns of brain-expressed miRNAs across Hymenoptera suggest a role for certain miRNAs in the evolution of the Aculeata, which includes all the eusocial hymenopteran species. Our results support the intriguing hypothesis that miRNAs are important regulators of social behavior at both developmental and evolutionary time scales.

Characterization and comparative analysis of small RNAs in three small RNA libraries of the brown planthopper (Nilaparvata lugens)

BACKGROUND

The brown planthopper (BPH), Nilaparvata lugens (Stå;l), which belongs to Homopteran, Delphacidae, is one of the most serious and destructive pests of rice. Feeding BPH with homologous dsRNA in vitro can lead to the death of BPH, which gives a valuable clue to the prevention and control of this pest, however, we know little about its small RNA world.

METHODOLOGY/PRINCIPAL FINDINGS

Small RNA libraries for three developmental stages of BPH (CX-male adult, CC-female adult, CY-last instar female nymph) had been constructed and sequenced. It revealed a prolific small RNA world of BPH. We obtained a final list of 452 (CX), 430 (CC), and 381 (CY) conserved microRNAs (miRNAs), respectively, as well as a total of 71 new miRNAs in the three libraries. All the miRNAs had their own expression profiles in the three libraries. The phylogenic evolution of the miRNA families in BPH was consistent with other species. The new miRNA sequences demonstrated some base biases.

CONCLUSION

Our study discovered a large number of small RNAs through deep sequencing of three small RNA libraries of BPH. Many animal-conserved miRNA families as well as some novel miRNAs have been detected in our libraries. This is the first achievement to discover the small RNA world of BPH. A lot of new valuable information about BPH small RNAs has been revealed which was helpful for studying insect molecular biology and insect resistant research.

Evolution and function of the extended miR-2 microRNA family

MicroRNAs are essential post-transcriptional regulators. Many animal microRNAs are clustered in the genome, and it has been shown that clustered microRNAs may be transcribed as a single transcript. Polycistronic microRNAs are often members of the same family, suggesting a role of tandem duplication in the emergence of clusters. The mir-2 microRNA family is the largest in Drosophila melanogaster, with 8 members that are mostly clustered in the genome. Previous studies suggest that the copy number and genomic distribution of mir-2 family members has been subject to significant change during evolution. The effects of such changes on their function are still unknown. Here we study the evolution of function in the mir-2 family. Our analyses show that, in spite of the change in number and organization among invertebrates, most mir-2 loci produce very similar mature microRNA products. Multiple mature miR-2 sequences are predicted to target genes involved in neural development in Drosophila. These targeting properties are conserved in the distant species Caenorhabditis elegans. Duplication followed by functional diversification is frequent during protein-coding gene evolution. However, our results suggest that the production of microRNA clusters by gene duplication rarely involves functional changes. This pattern of functional redundancy among clustered paralogous microRNAs reflects birth-and-death evolutionary dynamics. However, we identified a small number of mir-2 sequences in Drosophila that may have undergone functional shifts associated with genomic rearrangements. Therefore, redundancy in microRNA families may facilitate the acquisition of novel functional features.

Detection of microRNAs in color space

Motivation: Deep sequencing provides inexpensive opportunities to characterize the transcriptional diversity of known genomes. The AB SOLiD technology generates millions of short sequencing reads in color-space; that is, the raw data is a sequence of colors, where each color represents 2 nt and each nucleotide is represented by two consecutive colors. This strategy is purported to have several advantages, including increased ability to distinguish sequencing errors from polymorphisms. Several programs have been developed to map short reads to genomes in color space. However, a number of previously unexplored technical issues arise when using SOLiD technology to characterize microRNAs.

Results: Here we explore these technical difficulties. First, since the sequenced reads are longer than the biological sequences, every read is expected to contain linker fragments. The color-calling error rate increases toward the 3^′ end of the read such that recognizing the linker sequence for removal becomes problematic. Second, mapping in color space may lead to the loss of the first nucleotide of each read. We propose a sequential trimming and mapping approach to map small RNAs. Using our strategy, we reanalyze three published insect small RNA deep sequencing datasets and characterize 22 new microRNAs.

Availability and implementation: A bash shell script to perform the sequential trimming and mapping procedure, called SeqTrimMap, is available at: http://www.mirbase.org/tools/seqtrimmap/

Contact: antonio.marco@manchester.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.