MicroRNAs show diverse and dynamic expression patterns in multiple tissues of Bombyx mori

Integrated analysis of miRNA and mRNA expression profiles in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis

Insect gut microbiota has been reported to participate in regulating host multiple biological processes including metabolism and reproduction. However, the corresponding molecular mechanisms remain largely unknown. Recent studies suggest that microRNAs (miRNAs) are involved in complex interactions between the gut microbiota and the host. Here, we used next-generation sequencing technology to characterize miRNA and mRNA expression profiles and construct the miRNA-gene regulatory network in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis. A total of 3016 differentially expressed genes (DEGs) and 18 differentially expressed miRNAs (DEMs) were identified. Based on the integrated analysis of miRNA and mRNA sequencing data, 229 negatively correlated miRNA-gene pairs were identified from the miRNA-mRNA network. Gene ontology enrichment analysis indicated that DEMs could target several genes involved in the metabolic process, oxidation-reduction process, oogenesis, and insulin signaling pathway. Finally, real-time quantitative polymerase chain reaction further verified the accuracy of RNA sequencing results. In conclusion, our study provides the profiles of miRNA and mRNA expressions under antibiotics treatment and provides an insight into the roles of miRNAs and their target genes in the interaction between the gut microbiota and its host.

miR-275/305 cluster is essential for maintaining energy metabolic homeostasis by the insulin signaling pathway in Bactrocera dorsalis

Increasing evidence indicates that miRNAs play crucial regulatory roles in various physiological processes of insects, including systemic metabolism. However, the molecular mechanisms of how specific miRNAs regulate energy metabolic homeostasis remain largely unknown. In the present study, we found that an evolutionarily conserved miR-275/305 cluster was essential for maintaining energy metabolic homeostasis in response to dietary yeast stimulation in Bactrocera dorsalis. Depletion of miR-275 and miR-305 by the CRISPR/Cas9 system significantly reduced triglyceride and glycogen contents, elevated total sugar levels, and impaired flight capacity. Combined in vivo and in vitro experiments, we demonstrated that miR-275 and miR-305 can bind to the 3'UTR regions of SLC2A1 and GLIS2 to repress their expression, respectively. RNAi-mediated knockdown of these two genes partially rescued metabolic phenotypes caused by inhibiting miR-275 and miR-305. Furthermore, we further illustrated that the miR-275/305 cluster acting as a regulator of the metabolic axis was controlled by the insulin signaling pathway. In conclusion, our work combined genetic and physiological approaches to clarify the molecular mechanism of metabolic homeostasis in response to different dietary stimulations and provided a reference for deciphering the potential targets of physiologically important miRNAs in a non-model organism.

Identification and expression analysis of sex biased miRNAs in chinese hook snout carp Opsariichthys bidens

As an economically important fish, Opsariichthys bidens has obvious sexual dimorphism and strong reproductive capacity, but no epigenetics study can well explain its phenotypic variations. In recent years, many microRNAs involved in the regulation of reproductive development have been explored. In this study, the small RNA libraries of O. bidens on the testis and ovary were constructed and sequenced. A total of 295 known miRNAs were obtained and 100 novel miRNAs were predicted. By comparing testis and ovary libraries, 115 differentially expressed (DE) miRNAs were selected, of which 53 were up-regulated and 62 were down-regulated. A total of 64 GO items (padj < 0.01) and 206 KEGG pathways (padj < 0.01) were enriched in the target gene of miRNA. After that, the expression levels of nine DE miRNAs, including let-7a, miR-146b, miR-18c, miR-202-5p, miR-135c, miR-9-5p, miR-34c-3p, miR-460-5p and miR-338 were verified by qRT-PCR. Furthermore, bidirectional prediction of DE miRNAs and sex-related genes was carried out and the targeting correlation between miR-9-5p and nanos1 was verified by Dual-Luciferase reporter assay. Our findings identified the differentially expressed miRNA and paved the way to new possibilities for the follow-up study on the mechanism of miRNA-mRNA interaction in the gonads of O. bidens.

A dual-luciferase reporter system for characterization of small RNA target genes in both mammalian and insect cells

MicroRNAs (miRNAs) are regulatory RNA molecules that bind to target messenger RNAs (mRNAs) and affect the stability or translational efficiency of the bound mRNAs. Single or dual-luciferase reporter systems have been successfully used to identify miRNA target genes in mammalian cells. These reporter systems, however, are not sensitive enough to verify miRNA-target gene relationships in insect cell lines because the promoters of the target luciferase (usually Renilla) used in these reporter systems are too weak to drive sufficient expression of the target luciferase in insect cells. In this study, we replaced the SV40 promoter in the psiCHECK-2 reporter vector, which is widely used with mammalian cell lines, with the HSV-TK or AC5.1 promoter to yield two new dual-luciferase reporter vectors, designated psiCHECK-2-TK and psiCHECK-2-AC5.1, respectively. Only psiCHECK-2 and psiCHECK-2-AC5.1 had suitable target (Renilla)/reference (firefly) luciferase activity ratios in mammalian (HeLa and HEK293) and insect (Sf9, S2, Helicoverpa zea fat body and ovary) cell lines, while psiCHECK-2-TK had suitable Renilla/firefly luciferase activity ratios regardless of the cell line. Moreover, psiCHECK-2-TK successfully detected the interaction between Helicoverpa armigera miRNA9a and its target, the 3'-untranslated region of heat shock protein 90, in both mammalian and H. zea cell lines, but psiCHECK-2 failed to do so in H. zea cell lines. Furthermore, psiCHECK-2-TK with the target sequence, HzMasc (H. zea Masculinizer), accurately differentiated between H. zea cell lines with or without the negative regulation factor (miRNA or piRNA) of HzMasc. These data demonstrate that psiCHECK-2-TK can be used to functionally characterize small RNA target genes in both mammalian and insect cells.

Identification and characterization of microRNAs in the immature stage of the beneficial predatory bug Arma chinensis Fallou (Hemiptera: Pentatomidae)

MicroRNAs (miRNAs) are a type of small noncoding RNAs that regulate gene expression at the posttranscriptional level and can influence significant biological processes. Arma chinensis (Hemiptera: Pentatomidae) is a predaceous insect species that preys upon a wide variety of insect pests. It is important to explore and understand the molecular mechanisms involving miRNAs in regulating developmental and other gene expression for beneficial insects. However, examination of miRNAs associated with Hemiptera, especially predatory bugs, has been absent or scarce. This study represents the first comprehensive analysis of predatory bug A. chinensis encoded miRNAs through high throughput sequencing and predicts genes and biological processes regulated by the newly identified miRNAs through analyzing their differential expression in and across five nymphal instars. A total of 64 A. chinensis miRNAs, including 46 conserved miRNAs and 18 novel miRNAs, were identified by analysis of high throughput sequence reads mapped to the genome. A total of 2913 potential gene targets for these 64 miRNAs were predicted by comprehensive analyses utilizing miRanda, PITA, and RNAhybrid. Gene Ontology annotation of predicted target genes of A. chinensis suggested the key processes regulated by miRNAs involved biological processes, regulation of cellular processes, and transporter activity. Kyoto Encyclopedia of Genes and Genomes pathway predictions included the Toll and Imd signaling pathway, Valine, leucine and isoleucine degradation, Steroid biosynthesis, the AGE-RAGE signaling pathway in diabetic complications, and Alanine, aspartate and glutamate metabolism. This newly identified miRNAs through analyzing their differential expression, assessment of their predicted functions forms a foundation for further investigation of specific miRNAs.

20-Hydroxyecdysone-responsive microRNAs of insects

20-Hydroxyecdysone (20-HE) plays essential roles in coordinating developmental transitions of insects through responsive protein-coding genes and microRNAs (miRNAs). The involvement of single miRNAs in the ecdysone-signalling pathways has been extensively explored, but the interplay between ecdysone and the majority of miRNAs still remains largely unknown. Here, by small RNA sequencing, we systematically investigated the genome-wide responses of miRNAs to 20-HE in the embryogenic cell lines of Bombyx mori and Drosophila melanogaster. Over 60 and 70 20-HE-responsive miRNAs were identified in the BmE cell line and S2 cell line, respectively. The response of miRNAs to ecdysone exhibited a time-dependent pattern, and the response intensity increased with extending exposure to 20-HE. The relationship between ecdysone and the miRNAs was further explored through knockdown of ecdysone-signalling pathway genes. Specifically, ecdysone regulated the cluster miR-275 and miR-305 through the coordination of BmEcR-B and downstream BmE75B, and the interaction between BmEcR and miR-275 cluster was strengthened by the feedback regulation of BmE75B. Ecdysone induced miR-275-3p and miR-305-5p through the ecdysone response effectors (EcREs) at the upstream of the pre-miR-275 cluster. Overall, the results might help us further understand the relationship between ecdysone signalling pathways and small RNAs in the development and metamorphosis of insects.

A Day in the Life: Identification of Developmentally Regulated MicroRNAs in the Colorado Potato Beetle (Leptinotarsa decemlineata; Coleoptera: Chrysomelidae)

The Colorado potato beetle (Leptinotarsa decemlineata (Say)) is an important pest of the cultivated potato (Solanum tuberosum (L.) [Solanales: Solanaceae]). With its broad resistance toward commonly used insecticides, it is clear that more sophisticated control strategies are needed. Due to their importance in insect development, microRNAs (miRNAs) represent a potential tool to employ in insect control strategies. However, most studies conducted in this area have focused on model species with well-annotated genomes. In this study, next-generation sequencing was used to catalogue the miRNAs produced by L. decemlineata across all eight stages of its development, from eggs to adults. For most stages, the length of miRNAs peaked between 21 and 22 nt, though it was considerably longer for the egg stage (26 nt). Global profiling of miRNAs revealed three distinct developmental clusters: 1) egg stage; 2) early stage (first, second, and third instar); and 3) late stage (fourth instar, prepupae, pupae, and adult). We identified 86 conserved miRNAs and 33 bonafide novel miRNAs, including stage-specific miRNAs and those not previously identified in L. decemlineata. Most of the conserved miRNAs were found in multiple developmental stages, whereas the novel miRNAs were often stage specific with the bulk identified in the egg stage. The identified miRNAs have a myriad of putative functions, including growth, reproduction, and insecticide resistance. We discuss the putative roles of some of the most notable miRNAs in the regulation of L. decemlineata development, as well as the potential applications of this research in Colorado potato beetle management.

Transcriptome-Wide Identification of MicroRNAs and Analysis of Their Potential Roles in Development of Indian Meal Moth (Lepidoptera:Pyralidae)

MicroRNAs (miRNAs) have been reported to play indispensable roles in regulating various developmental processes via the posttranscriptional repression of target genes in insect species. In the present paper, we studied the miRNAs in Indian meal moth (Plodia interpunctella (Hübener)), one of the most economically important stored grains pests around the world. In total, 12 small RNA libraries from four developmental stages of P. interpunctella were constructed, and 178 known and 23 novel miRNAs were identified. In addition, the expression profiles of these miRNAs were assessed across different developmental stages and miRNAs that were highly expressed in eggs, larvae, pupae, and adults were identified. Specifically, 100, 61, and 52 miRNAs were differentially expressed between eggs and larvae, larvae and pupae, and pupae and adults, respectively. The KEGG and GO analysis of the predicted target genes suggested the essential roles of miRNAs in the regulation of complex development of P. interpunctella. Importantly, we also found a set of miRNAs might be involved in the larval metamorphic molting process, with their expressions increasing and then decreasing during the larva-pupa-adult stages of P. interpunctella. In conclusion, the current paper has discovered numerous miRNAs, and some key miRNAs that might be responsible for regulating development in P. interpunctella. To our knowledge, this is the first study to document miRNAs and their expression patterns in interpunctella, and those findings would lay an important molecular foundation for future functional analysis of these miRNAs in P. interpunctella.

Let-7 microRNA is a critical regulator in controlling the growth and function of silk gland in the silkworm

The silk gland is characterized by high protein synthesis. However, the molecular mechanisms controlling silk gland growth and silk protein synthesis remain undetermined. Here we demonstrated that CRISPR/Cas9-based knockdown of let-7 or the whole cluster promoted endoreduplication and enlargement of the silk gland, accompanied by changing silk yield, whereas transgenic overexpression of let-7 led to atrophy and degeneration of the silk gland. Mechanistically, let-7 controls cell growth in the silk gland through coordinating nutrient metabolism processes and energy signalling pathways. Transgenic overexpression of pyruvate carboxylase, a novel target of let-7, resulted in enlargement of the silk glands, which is consistent with the abnormal phenotype of the let-7 knockdown. Overall, our data reveal a previously unknown miRNA-mediated regulation of silk gland growth and physiology and shed light on involvement of let-7 as a critical stabilizer and booster in carbohydrate metabolism, which may have important implications for understanding of the molecular mechanism and physiological function of specialized organs in other species.

Bmo-miR-2780a regulates the expression of the sericin-1 gene of Bombyx mori

Silk production in Bombyx mori L. is largely determined by the expression of genes encoding fibroin and sericin. Here, we examined the regulatory function of a microRNA (miRNA) on silk gene expression using the sericin-1 gene (BmSer-1). First, we downloaded whole mature miRNAs of silkworm from miRBase and identified bmo-miR-2780a as a candidate miRNA for the regulation of BmSer-1 expression. We used semi-quantitative reverse transcription polymerase chain reaction (RT-PCR) with stem-loop primers to investigate the expression profile of bmo-miR-2780a and its predicted target gene BmSer-1 in seven different tissues from 5th instar day-3 larvae, including head, fat body, anterior silk gland (ASG), middle silk gland (MSG), posterior silk gland (PSG), middle gut, and hemolymph. Our results showed that bmo-miR-2780a was specifically expressed in the MSG and that the expression level of BmSer-1 was significantly higher in the MSG than in other tissues. Recombinant plasmids carrying both pri-mir-2780a and Ser1-3'UTR were constructed and then used to cotransfect BmN cells. We further detected the effect of bmo-miR-2780a on Ser-1 in vivo. These results showed that the target gene was significantly decreased by miR-2780a compared with the control group (p < .05), thus indicating that bmo-miR-2780a might negatively regulate the expression of Ser-1.

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.

Genetic Manipulation of MicroRNAs in the Silk Gland of Silkworm, Bombyx Mori

Background

MicroRNAs (miRNAs) are a class of non-coding RNAs with important post-transcriptional regulatory functions. To reveal the function of miRNAs in vivo, the critical step is to change their expression levels in the tissues or organs. In this work, we explored the application of several important genetic techniques in altering the expression of silk gland-specific miR-274 of silkworm (Bombyx mori).

Results

Injection of synthesized microRNA mimics and antagomirs exerted no effect on the expression of miR-274 in the silk gland, miR-274 sponge specifically absorbed miR-274 and down-regulated its expression, transgenic overexpression of miR-274 precursor significantly up-regulated miR-274, and finally tissue-specific CRISPR/Cas9 system achieved deletion of miR-274.

Conclusions

A practical technical system was established for studying the functions of miRNAs in silk gland of Bombyx mori. Our research provides methodological support for the functional study of miRNAs and other noncoding RNAs in the silk gland and more organs in other species.

Bmo-miR-79 downregulates the expression of BmEm4 in the silkworm, Bombyx mori

MicroRNA is an important regulation factor in insect development and metamorphosis. It has been reported that E(spl)m4 is a miRNA-targeted gene, as well as the target of the Notch signaling pathway in Drosophila. The expression of E(spl)m4 can be regulated by microRNA and further affect the neural development of Drosophila. Here, we found that BmEm4, an ortholog of E(spl)m4 from Bombyx mori, was the target gene of bmo-miR-79, with target sites containing the Brd and K boxes of the BmEm4_3'UTR, which was validated by the dual luciferase reporter (DLR) assay. Furthermore, bmo-miR-79 mimics can inhibit the expression of BmEm4 in BmN cells after transfection, and bmo-miR-79 can also inhibit the expression of BmEm4 in different developmental stages of Bombyx mori at a posttranscriptional level, to different degrees. The EMSA test further showed that bmo-miR-79 could bind to BmAGO2, which is the Bombyx mori argonaute2 protein, suggesting that bmo-miR-79 might regulate the expression of BmEm4 by forming miRISC complexes with BmAGO2. Taken together, bmo-miR-79 could regulate the expression of BmEm4 mediated by BmAGO2 and further affect its function in the silkworm Bombyx mori.

Identification and profiling of microRNAs during gonadal development in the giant freshwater prawn Macrobrachium rosenbergii

As post-transcriptional regulators, microRNAs (miRNAs) play an important role in growth and reproductive processes. So far, there is limited information regarding crustacean miRNAs. To explore the potential role of miRNAs in the gonadal development of the prawn Macrobrachium rosenbergii, we constructed seven small RNA libraries from ovarian and testicular tissues at various stages using somatic tissue as the control. A total of 1,954 known and 129 novel miRNAs were retrieved. By comparing differentially expressed miRNAs (DEMs) between testes and ovaries, forty-one miRNAs were identified with sex-biased expression patterns, including 17 ovary-biased and 24 testis-biased patterns. Furthermore, the putative target genes of the sex-biased miRNAs, such as cyclin L1, mitogen-activated protein kinase 7 (MAPK 7), heat shock protein (HSP), and zinc finger protein, were significantly enriched in many reproduction-related pathways including the Gonadotropin-releasing hormone (GnRH) pathway, glycolysis, gluconeogenesis pathway, ovarian steroidogenesis, estrogen signaling pathway, MAPK pathway, Wnt pathway, and insulin signaling pathway, implicating potential regulatory roles of miRNAs in reproduction. These data aid in the further investigation of the mechanism of gonadal development and reproductive regulation mediated by miRNA in M. rosenbergii.

Screening for differentially expressed miRNAs in Aedes albopictus (Diptera: Culicidae) exposed to DENV-2 and their effect on replication of DENV-2 in C6/36 cells

BACKGROUND

The mosquito Aedes albopictus is an important vector for dengue virus (DENV) transmission. The midgut is the first barrier to mosquito infection by DENV, and this barrier is a critical factor affecting the vector competence of the mosquito. However, the molecular mechanism of the interaction between midgut and virus is unknown.

RESULTS

Six small libraries of Ae. albopictus midgut RNAs were constructed, three of which from mosquitoes that were infected with DENV-2 after feeding on infected blood, and another three that remained uninfected with DENV-2 after feeding on same batch of infected blood. A total of 46 differentially expressed miRNAs were identified of which 17 significant differentially expressed miRNAs were selected. Compared to microRNA expression profiles of mosquitoes that were uninfected with DENV-2, 15 microRNAs were upregulated and two were downregulated in mosquitoes that were infected with DENV-2. Among these differentially expressed microRNAs, miR-1767, miR-276-3p, miR-4448 and miR-622 were verified by stem-loop qRT-PCR in samples from seven-day-infected and uninfected midguts and chosen for an in vitro transient transfection assay. miR-1767 and miR-276-3p enhanced dengue virus replication in C6/36 cells, and miR-4448 reduced dengue virus replication.

CONCLUSIONS

To our knowledge, this study is the first to reveal differences in expression levels between mosquitoes infected and uninfected with DENV-2 after feeding on an infected blood meal. It provides useful information on microRNAs expressed in the midgut of Aedes albopictus after exposure to the virus.

microRNAs in Macrobrachium olfersii embryos: Identification, their biogenesis components and potential targets

In embryonic development, microRNAs (miRNAs) regulate the complex gene expression associated with the complexity of embryogenesis. Today, few studies have been conducted on the identification of miRNAs and components of miRNA biogenesis on embryonic development in crustaceans, especially in prawns. In this context, the aim of this study was to identify in silico components of miRNA biogenesis, and miRNAs and potential target genes during embryonic development in the prawn Macrobrachium olfersii through small RNAs and transcriptome analyses. Using the miRDeep2 program, we identified 17 miRNA precursors in M. olfersii, which seven (miR-9, miR-10, miR-92, miR-125, miR-305, miR-1175, and miR-2788) were reported in the miRBase database, indicating high evolutionary conservation of these sequences among animals. The other 10 miRNAs of M. olfersii were novel miRNAs and only similar to Macrobrachium niponnense miRNAs, indicating genus-specific miRNAs. In addition, eight key components of miRNA biogenesis (DROSHA, PASHA/DGCR8, XPO5, RAN, DICER, TRBP2, AGO, and PIWI) were identified in M. olfersii embryos unigenes. In the annotation of miRNA targets, 516 genes were similar to known sequences in the GenBank database. Regarding the conserved miRNAs, we verified that they were differentially expressed during embryonic development in M. olfersii. In conclusion, this is the first study that identifies conserved and novel miRNAs in the prawn M. olfersii with some miRNA target genes involved in embryonic development. Our results will allow further studies on the function of these miRNAs and miRNA biogenesis components during embryonic development in M. olfersii and other prawns of commercial interest.

Expressing MicroRNA Bantam Sponge Drastically Improves the Insecticidal Activity of Baculovirus via Increasing the Level of Ecdysteroid Hormone in Spodoptera exigua Larvae

Bantam is a conserved miRNA highly expressed in insects. We previously showed that the antisense inhibitor (antagomiR) of bantam improved the infection by baculovirus Autographa californica nucleopolyhedrovirus (AcMNPV) in Spodoptera exigua and S. litura larvae. Here, we constructed a recombinant AcMNPV (vPH-banS) expressing bantam sponge, an mRNA containing eight antisense binding sites for bantam. Infection with wild type AcMNPV (WT) or the control recombinant virus vPH resulted in a significant increase of bantam level, whereas infection with vPH-banS led to an approximately 40% reduction of bantam in both Sf9 cells and S. exigua larvae. Although, comparable production of budded virus and polyhedra were detected in vPH-banS-, vPH-, and WT-infected Sf9 cells, vPH-banS showed remarkably increased insecticidal activity in S. exigua larvae. The 50% lethal concentration and the median lethal time of vPH-banS was only 1/40 and 1/2, respectively, of both vPH and WT. Further analysis showed that the level of molting hormone 20-hydroxyecdysone (20E) was significantly higher in larvae infected with vPH-banS than those infected with vPH or WT. This was confirmed by the result that the larvae treated with bantam inhibitor also had a markedly increased 20E level. Moreover, feeding larvae with 20E increased the virus-mediated mortality, whereas feeding with juvenile hormone partially reverted the high insecticidal effect of vPH-banS. Together, our results revealed that vPH-banS infection suppresses the level of bantam, and in turn elevates level of 20E in infected insects, resulting in increased susceptibility to baculovirus infection. Our study provided a novel approach to improve a baculovirus bio-insecticide by interfering with a key homeostasis-regulating miRNA of the host.

microRNA profiling between Bacillus thuringiensis Cry1Ab-susceptible and -resistant European corn borer, Ostrinia nubilalis (Hübner)

Transgenic maize hybrids that express insecticidal Bacillus thuringiensis (Bt) crystalline (Cry) protein toxins effectively protect against the European corn borer, Ostrinia nubilalis, a devastating maize pest. Field monitoring and laboratory selections have detected varying levels of O. nubilalis resistance to Cry1Ab toxin. MicroRNAs (miRNAs) are short noncoding RNAs that are involved in post-transcriptional gene regulation. Their potential roles in the evolution of Bt resistance, however, remain largely unknown. Sequencing of small RNA libraries from the midgut of Cry1Ab-susceptible and resistant O. nubilalis larvae resulted in the discovery of 277 miRNAs, including 248 conserved and 29 novel. Comparative analyses of miRNA expression profiles between the laboratory strains predicted 26 and nine significantly up- and down-regulated transcripts, respectively, in the midgut of Cry1Ab resistant larvae. Amongst 15 differentially regulated miRNAs examined by quantitative real-time PCR, nine (60%) were validated as cosegregating with Cry1Ab resistance in a backcross progeny. Differentially expressed miRNAs were predicted to affect transcripts involved in cell membrane components with functions in metabolism and binding, and the putative Bt-resistance genes aminopeptidase N and cadherin. These results lay the foundation for future investigation of the potential role of miRNAs in the evolution of Bt resistance.

Discovery and functional characterization of microRNAs and their potential roles for gonadal development in spotted knifejaw, Oplegnathus punctatus

The spotted knifejaw (Oplegnathus punctatus) is a newly emerging economical fishery species in China. Studies focused on the regulation of gonadal development and gametogenesis of spotted knifejaw are still insufficient. As a key post-transcriptional regulator, miRNAs have been shown to play important roles in development and reproduction systems. In this study, small RNA deep sequencing in ovary and testis of spotted knifejaw were performed to screen miRNA expression patterns. After sequencing and bioinformatics analysis, a total of 247 conserved known miRNAs and 41 novel miRNAs were identified in spotted knifejaw gonads for the first time. In addition, 36 miRNAs were differentially expressed between testis and ovary. The putative target genes of differentially expressed (DE) miRNAs were significantly enriched in several pathways related to sexual differentiation and gonadal development, such as steroid hormone biosynthesis. Sequencing data was validated through qRT-PCR analysis of selected DE miRNAs. Dual-luciferase reporter analyses of filtered miRNA-target gene pairs confirmed that opu-miR-27b-3p targeted in piwi2 and mov10l1 3' UTRs and down-regulated their expressions in spotted knifejaw. The notion that mov10l1 and piwi2 enhance germ cells proliferation and regulate gonadal development and gametogenesis suggests that opu-miR-27b-3p may attenuated this process in the gonads of spotted knifejaw. These findings provided insights into regulatory roles of gonadal miRNAs and supplied fundamental resources for further studies on miRNA-mediated post-transcriptional regulation in reproductive system of spotted knifejaw.

A comparative analysis of corpora allata-corpora cardiaca microRNA repertoires revealed significant changes during mosquito metamorphosis

The corpora allata (CA) are a pair of endocrine glands with neural connections to the brain and close association with another neuroendocrine organ, the corpora cardiaca (CC). The CA from adult female Aedes aegypti mosquitoes synthesize fluctuating levels of juvenile hormone (JH), which have been linked to the ovarian development and are influenced by nutritional signals. In this study, we investigated the potential involvement of microRNAs (miRNAs), a type of small non-coding RNAs, in the regulation of gene expression in CA-CC complexes during mosquito reproductive development, at stages with distinct JH biosynthesis patterns. We analyzed the miRNA repertoires expressed in the CA-CC of pupae, sugar-fed and blood-fed female Ae. aegypti. In total, 156 mature miRNAs were detected in the CA-CC, with 84 displaying significant differences in expression among the three CA-CC developmental stages. There were more miRNAs that were expressed in pupae, and decreased or were absent after adult emergence, when compared with changes between CA-CC of sugar and blood-fed females. Analysis of the genes identified as potential targets for the CA-CC miRNA repertoires classified them into the broad categories of metabolism, information storage and processing, and cellular processes and signaling; with genes involved in cellular processes and signaling representing the largest portion. Among them, the signal-transduction mechanisms and intracellular trafficking, secretion and vesicular transport contained almost 55% of the genes' targets. A substantial number of miRNAs were differentially abundant in the libraries of the three developmental stages, and those changes were much more notable when pupae and adult stages were compared. We detected putative binding sites for some of the most abundant miRNAs on genes encoding JH biosynthetic enzymes and CC neuropeptides. These studies should help us to gain a better understanding of the regulation of CA-CC activity mediated by miRNAs during major developmental stages in mosquitoes.

MicroRNA profile of silk gland reveals different silk yields of three silkworm strains

Silk proteins are synthesized and secreted by the silk gland. The differential gene expression in it leads to different silk yield among various silkworm strains. As crucial factors, microRNAs (miRNAs) regulate protein synthesis at post-transcriptional level in silk gland. MiRNAs expression level in the silk gland of three silkworm strains (Jingsong, Lan10 and Dazao) was analyzed and 33 differentially expressed miRNAs (DEMs) were discovered between JingSong (JS) and Lan10 (L10), 60 DEMs between JS and Dazao, 54 DEMs between L10 and Dazao respectively. The DEMs target genes were predicted combing with two different methods and their functions were annotated according to gene ontology. Our previous studies showed that a batch of genes related to silk yield were identified in JS and L10 strains by comparative transcriptome and quantitative trait loci (QTL) method. Thirteen DEMs whose target genes are related to protein biosynthesis processes were screened by combining with these researches. Twelve DEMs potentially regulate nineteen genes which exist in our QTL results. Six common DEMs potentially regulate the genes in both of previous results. Finally, five DEMs were selected to verify their expression levels between JS and L10 by qRT-PCR, which showed similar difference as the results of small RNA-sequencing. MiRNAs in the silk gland may directly affect silk protein biosynthesis in different silkworm strains. In current work, we identified a batch of DEMs which potentially regulate the genes related to silk yield. Further functionally study of these miRNAs will contribute to improve varieties and boost the silk yield. Our research provides a basis for studying these miRNAs and their functions in silk production.

Integrated analysis of microRNA and mRNA expression profiles during the sex-differentiation sensitive period in oriental river prawn, Macrobrachium nipponense

Male oriental river prawns (Macrobrachium nipponense) grow faster than females, and therefore, reach larger sizes by harvest time. Histological observations have indicated that the sex-differentiation sensitive period (which includes the formation of the androgenic gland, the testis, and the ovary) is from post-larvae (PL) developmental stage for M. nipponense. In this study, we prepared four microRNA (miRNA) and mRNA libraries using samples collected from sex-differentiation sensitive period (PL7 to PL16) to perform RNA-sequencing for identifying sex-related candidate miRNAs, genes, and metabolic pathways. A total of nine intersection miRNAs were identified, of which three were highly expressed in the androgenic gland, and their expression was verified by quantitative Real-Time PCR (qPCR). These three miRNAs and their 11 predicted target genes may be strong candidates for sex-related miRNAs and sex-related genes in M. nipponense. Five vital sex-related metabolic pathways were also identified that may regulate other sex-differentiation and sex-determination mechanisms. Finding of the study provide important insights to enhance our understanding on sex-differentiation and sex-determination mechanisms for M. nipponense.

De novo characterization of microRNAs in oriental fruit moth Grapholita molesta and selection of reference genes for normalization of microRNA expression

MicroRNAs (miRNAs) are a group of endogenous non-coding small RNAs that have critical regulatory functions in almost all known biological processes at the post-transcriptional level in a variety of organisms. The oriental fruit moth Grapholita molesta is one of the most serious pests in orchards worldwide and threatens the production of Rosacea fruits. In this study, a de novo small RNA library constructed from mixed stages of G. molesta was sequenced through Illumina sequencing platform and a total of 536 mature miRNAs consisting of 291 conserved and 245 novel miRNAs were identified. Most of the conserved and novel miRNAs were detected with moderate abundance. The miRNAs in the same cluster normally showed correlated expressional profiles. A comparative analysis of the 79 conserved miRNA families within 31 arthropod species indicated that these miRNA families were more conserved among insects and within orders of closer phylogenetic relationships. The KEGG pathway analysis and network prediction of target genes indicated that the complex composed of miRNAs, clock genes and developmental regulation genes may play vital roles to regulate the developmental circadian rhythm of G. molesta. Furthermore, based on the sRNA library of G. molesta, suitable reference genes were selected and validated for study of miRNA transcriptional profile in G. molesta under two biotic and six abiotic experimental conditions. This study systematically documented the miRNA profile in G. molesta, which could lay a foundation for further understanding of the regulatory roles of miRNAs in the development and metabolism in this pest and might also suggest clues to the development of genetic-based techniques for agricultural pest control.

The Effect of MicroRNA bantam on Baculovirus AcMNPV Infection in Vitro and in Vivo

The role of microRNA bantam, one of the most abundant microRNAs in Sf9 cells, was studied for its role in baculovirus infection in vitro and in vivo. The expression level of bantam was increased after AcMNPV infection in Sf9 cells and in Spodoptera litura larvae. In Sf9 cells, application of bantam inhibitor or mimic altered the expression of many virus genes, the most affected gene being lef8, gp41 and p10, the expression level of which was increased by 8, 10 and 40 times, respectively, in the presence of bantam inhibitor. Virus DNA replication was decreased in the presence of bantam mimic and increased in the presence of bantam inhibitor in a dose dependent manner. However, the production of budded virus did not change significantly. Feeding the larvae of S. litura and Spodoptera exigua with bantam antagomiR, a more stable form of the inhibitor, resulted in an abnormal larval growth and a decreased pupation rate. In S. litura, larvae died 3.5 days sooner than the control when bantam antagomiR was applied, together with AcMNPV. In infected S. exigua, larval mortality increased from 47% without antagomiR to 80% with it. The results suggest that microRNA bantam plays an important role in insect growth, as well as in baculovirus-insect interaction.

MicroRNAs of the mesothorax in Qinlingacris elaeodes, an alpine grasshopper showing a wing polymorphism with unilateral wing form

MicroRNAs (miRNAs) are now recognized as key post-transcriptional regulators in regulation of phenotypic diversity. Qinlingacris elaeodes is a species of the alpine grasshopper, which is endemic to China. Adult individuals have three wing forms: wingless, unilateral-winged and short-winged. This is an ideal species to investigate the phenotypic plasticity, development and evolution of insect wings because of its case of unilateral wing form in both the sexes. We sequenced a small RNA library prepared from mesothoraxes of the adult grasshoppers using the Illumina deep sequencing technology. Approximately 12,792,458 raw reads were generated, of which the 854,580 high-quality reads were used only for miRNA identification. In this study, we identified 49 conserved miRNAs belonging to 41 families and 69 species-specific miRNAs. Moreover, seven miRNA*s were detected both for conserved miRNAs and species-specific miRNAs, which were supported by hairpin forming precursors based on polymerase chain reaction. This is the first description of miRNAs in alpine grasshoppers. The results provide a useful resource for further studies on molecular regulation and evolution of miRNAs in grasshoppers. These findings not only enrich the miRNAs for insects but also lay the groundwork for the study of post-transcriptional regulation of wing forms.

bmo-miR-0001 and bmo-miR-0015 down-regulate expression of Bombyx mori fibroin light chain gene in vitro

Based on bioinformatic analysis, we selected two novel microRNAs (miRNAs), bmo-miR-0001 and bmo-miR-0015, from high-throughput sequencing of the Bombyx mori larval posterior silk gland (PSG). Firstly, we examined the expression of bmo-miR-0001 and bmo-miR 12 different tissues of the 5th instar Day-3 larvae of the silkworm. The results showed that the expression levels of both bmo-miR-0001 and bmo-miR-0015 were obviously higher in the PSG than in other tissues, implying there is a spatio-temporal condition for bmo-miR-0001 and bmo-miR-0015 to regulate the expression of BmFib-L. To test this hypothesis, we constructed pri-bmo-miR-0001 expressing the plasmid pcDNA3.0 and pri-bmo-miR-0015 expressing the plasmid pcDNA3.0 [ie1-egfp-pri-bmo-miR-0015-SV40]. Finally, the BmN cells were harvested and luciferase activity was detected. The results showed that luciferase activity was reduced significantly (P<0.05) in BmN cells co-transfected by pcDNA3.0 [ie1-egfp-pri-bmo-miR-0001-SV40] or pcDNA3.0 with pGL3.0 [A3-luc-Fib-L-3'UTR-SV40], suggesting that both bmo-miR-0001 and bmo-miR-0015 can down-regulate the expression of BmFib-L in vitro.

Systematic Identification and Characterization of Long Non-Coding RNAs in the Silkworm, Bombyx mori

Long noncoding RNAs (lncRNAs) are emerging as important regulators in various biological processes. However, to date, no systematic characterization of lncRNAs has been reported in the silkworm Bombyx mori. In the present study, we generated eighteen RNA-seq datasets with relatively high depth. Using an in-house designed lncRNA identification pipeline, 11,810 lncRNAs were identified for 5,556 loci. Among these lncRNAs, 474 transcripts were intronic lncRNAs (ilncRNAs), 6,250 transcripts were intergenic lncRNAs (lincRNAs), and 5,086 were natural antisense lncRNAs (lncNATs). Compared with protein-coding mRNAs, silkworm lncRNAs are shorter in terms of full length but longer in terms of exon and intron length. In addition, lncRNAs exhibit a lower level of sequence conservation, more repeat sequences overlapped and higher tissue specificity than protein-coding mRNAs in the silkworm. We found that 69 lncRNA transcripts from 33 gene loci may function as miRNA precursors, and 104 lncRNA transcripts from 72 gene loci may act as competing endogenous RNAs (ceRNAs). In total, 49.47% of all gene loci (2,749/5,556) for which lncRNAs were identified showed sex-biased expression. Co-expression network analysis resulted in 19 modules, 12 of which revealed relatively high tissue specificity. The highlighted darkgoldenrod module was specifically associated with middle and posterior silk glands, and the hub lncRNAs within this module were co-expressed with proteins involved in translation, translocation, and secretory processes, suggesting that these hub lncRNAs may function as regulators of the biosynthesis, translocation, and secretion of silk proteins. This study presents the first comprehensive genome-wide analysis of silkworm lncRNAs and provides an invaluable resource for genetic, evolutionary, and genomic studies of B. mori.

MicroRNA of the fifth-instar posterior silk gland of silkworm identified by Solexa sequencing

No special studies have been focused on the microRNA (miRNA) in the fifth-instar posterior silk gland of Bombyx mori. Here, using next-generation sequencing, we acquired 93.2 million processed reads from 10 small RNA libraries. In this paper, we tried to thoroughly describe how our dataset generated from deep sequencing which was recently published in BMC genomics. Results showed that our findings are largely enriched silkworm miRNA depository and may benefit us to reveal the miRNA functions in the process of silk production.

Sex-biased gene expression and sexual conflict throughout development

Sex-biased gene expression is likely to account for most sexually dimorphic traits because males and females share much of their genome. When fitness optima differ between sexes for a shared trait, sexual dimorphism can allow each sex to express their optimum trait phenotype, and in this way, the evolution of sex-biased gene expression is one mechanism that could help to resolve intralocus sexual conflict. Genome-wide patterns of sex-biased gene expression have been identified in a number of studies, which we review here. However, very little is known about how sex-biased gene expression relates to sex-specific fitness and about how sex-biased gene expression and conflict vary throughout development or across different genotypes, populations, and environments. We discuss the importance of these neglected areas of research and use data from a small-scale experiment on sex-specific expression of genes throughout development to highlight potentially interesting avenues for future research.

MicroRNAs as mediators of insect host-pathogen interactions and immunity

Insects are the most successful group of animals on earth, owing this partly to their very effective immune responses to microbial invasion. These responses mainly include cellular and humoral responses as well as RNA interference (RNAi). Small non-coding RNAs (snRNAs) produced through RNAi are important molecules in the regulation of gene expression in almost all living organisms; contributing to important processes such as development, differentiation, immunity as well as host-microorganism interactions. The main snRNAs produced by the RNAi response include short interfering RNAs, microRNAs and piwi-interacting RNAs. In addition to the host snRNAs, some microorganisms encode snRNAs that affect the dynamics of host-pathogen interactions. In this review, we will discuss the latest developments in regards to the role of microRNA in insect host-pathogen interactions and provide some insights into this rapidly developing area of research.

Development and immunity-related microRNAs of the lepidopteran model host Galleria mellonella

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNAs that act as key players in the post-transcriptional regulation of protein synthesis. Although little is known about their role in complex physiological processes such as development and immunity, our knowledge is expanding rapidly, thanks to the use of model systems. The larvae of the greater wax moth Galleria mellonella are now established as model hosts for pathogens that infect insects or humans. To build on our previously-reported comprehensive G. mellonella transcriptome, here we describe the identification and analysis of development and immunity-related miRNAs, thus providing valuable additional data to promote the use of this model host for the analysis of complex processes.

RESULTS

To screen for miRNAs that are differentially expressed in G. mellonella (1) during metamorphosis or (2) following infection with the entomopathogenic bacterium Serratia entomophila or (3) with the parasitic fungus Metarhizium anisopliae, we designed a microarray containing more than 2000 insect miRNA probe sequences. We identified miRNAs that were significantly expressed in pre-pupae (16), pupae (22) and last-instar larvae infected with M. anisopliae (1) in comparison with untreated last-instar larvae which were used as a reference. We then used our transcriptomic database to identify potential 3' untranslated regions that form miRNA-mRNA duplexes by considering both base pair complementarity and minimum free energy hybridization. We confirmed the co-expression of selected miRNAs (such as miR-71, miR-263a and miR-263b) with their predicted target mRNAs in last-instar larvae, pre-pupae and pupae by RT-PCR. We also identified miRNAs that were expressed in response to infection with bacterial or fungal pathogens, and one miRNA that may act as a candidate mediator of trans-generational immune priming.

CONCLUSIONS

This is the first study to identify miRNAs that are predicted to regulate genes expressed during metamorphosis or in response to infection in the lepidopteran model host G. mellonella.

MicroRNA expression profiling of the fifth-instar posterior silk gland of Bombyx mori

BACKGROUND

The growth and development of the posterior silk gland and the biosynthesis of the silk core protein at the fifth larval instar stage of Bombyx mori are of paramount importance for silk production.

RESULTS

Here, aided by next-generation sequencing and microarry assay, we profile 1,229 microRNAs (miRNAs), including 728 novel miRNAs and 110 miRNA/miRNA* duplexes, of the posterior silk gland at the fifth larval instar. Target gene prediction yields 14,222 unique target genes from 1,195 miRNAs. Functional categorization classifies the targets into complex pathways that include both cellular and metabolic processes, especially protein synthesis and processing.

CONCLUSION

The enrichment of target genes in the ribosome-related pathway indicates that miRNAs may directly regulate translation. Our findings pave a way for further functional elucidation of these miRNAs and their targets in silk production.

Overview of research on Bombyx mori microRNA

MicroRNAs (miRNAs) constitute some of the most significant regulatory factors involved at the post-transcriptional level after gene expression, contributing to the modulation of a large number of physiological processes such as development, metabolism, and disease occurrence. This review comprehensively and retrospectively explores the literature investigating silkworm, Bombyx mori L. (Lepidoptera: Bombicidae), miRNAs published to date, including discovery, identification, expression profiling analysis, target gene prediction, and the functional analysis of both miRNAs and their targets. It may provide experimental considerations and approaches for future study of miRNAs and benefit elucidation of the mechanisms of miRNAs involved in silkworm developmental processes and intracellular activities of other unknown non-coding RNAs.

Identification and developmental profiling of microRNAs in diamondback moth, Plutellaxylostella (L.)

MicroRNAs (miRNAs) are a group of small RNAs involved in various biological processes through negative regulation of mRNAs at the post-transcriptional level. Although miRNA profiles have been documented in over two dozen insect species, few are agricultural pests. In this study, both conserved and novel miRNAs in the diamondback moth, Plutella xylostella L., a devastating insect pest of cruciferous crops worldwide, were documented. High-throughput sequencing of a small RNA library constructed from a mixed life stages of P. xylostella, including eggs, 1st to 4th (last) instar larvae, pupae and adults, identified 384 miRNAs, of which 174 were P. xylostella specific. In addition, temporal expressions of 234 miRNAs at various developmental stages were investigated using a customized microarray analysis. Among the 91 differentially expressed miRNAs, qRT-PCR analysis was used to validate highly expressed miRNAs at each stage. The combined results not only systematically document miRNA profiles in an agriculturally important insect pest, but also provide molecular targets for future functional analysis and, ultimately, genetic-based pest control practice.

RIP-seq of BmAgo2-associated small RNAs reveal various types of small non-coding RNAs in the silkworm, Bombyx mori

Background

Small non-coding RNAs (ncRNAs) are important regulators of gene expression in eukaryotes. Previously, only microRNAs (miRNAs) and piRNAs have been identified in the silkworm, Bombyx mori. Furthermore, only ncRNAs (50-500nt) of intermediate size have been systematically identified in the silkworm.

Results

Here, we performed a systematic identification and analysis of small RNAs (18-50nt) associated with the Bombyx mori argonaute2 (BmAgo2) protein. Using RIP-seq, we identified various types of small ncRNAs associated with BmAGO2. These ncRNAs showed a multimodal length distribution, with three peaks at ~20nt, ~27nt and ~33nt, which included tRNA-, transposable element (TE)-, rRNA-, snoRNA- and snRNA-derived small RNAs as well as miRNAs and piRNAs. The tRNA-derived fragments (tRFs) were found at an extremely high abundance and accounted for 69.90% of the BmAgo2-associated small RNAs. Northern blotting confirmed that many tRFs were expressed or up-regulated only in the BmNPV-infected cells, implying that the tRFs play a prominent role by binding to BmAgo2 during BmNPV infection. Additional evidence suggested that there are potential cleavage sites on the D, anti-codon and TψC loops of the tRNAs. TE-derived small RNAs and piRNAs also accounted for a significant proportion of the BmAgo2-associated small RNAs, suggesting that BmAgo2 could be involved in the maintenance of genome stability by suppressing the activities of transposons guided by these small RNAs. Finally, Northern blotting was also used to confirm the Bombyx 5.8 s rRNA-derived small RNAs, demonstrating that various novel small RNAs exist in the silkworm.

Conclusions

Using an RIP-seq method in combination with Northern blotting, we identified various types of small RNAs associated with the BmAgo2 protein, including tRNA-, TE-, rRNA-, snoRNA- and snRNA-derived small RNAs as well as miRNAs and piRNAs. Our findings provide new clues for future functional studies of the role of small RNAs in insect development and evolution.

MicroRNA-281 regulates the expression of ecdysone receptor (EcR) isoform B in the silkworm, Bombyx mori

Insect development and metamorphosis are regulated by the coordination of ecdysone and juvenile hormones. Insect microRNAs (miRNAs) also act in insect development and metamorphosis by regulating genes in the ecdysone cascade. Although hundreds of insect miRNAs have been identified, the physiological functions of most remain poorly understood. Here, we report that a conserved insect miRNA, microRNA-281 (miR-281), regulates the ecdysone receptor (EcR), in an isoform-specific manner in the silkworm Bombyx mori. The B. mori EcR (BmEcR) gene encodes three isoforms: BmEcR-A, BmEcR-B1 and BmEcR-B2. The 3'UTR regions of A and B genes, which contain multiple potential microRNA targeting sites, are distinct. Target prediction revealed that miR-281 may specifically target the 3'UTR of BmEcR-B. Using a dual luciferase reporter assay in HEK293T cells, we confirmed that miR-281 suppressed transcription of BmEcR-B but not BmEcR-A. The expression of miR-281 and BmEcR-B are well coordinated in the Malpighian tubules from the fourth larval molt to pupation. In the Malpighian tubules of fifth instar larvae, BmEcR-B protein expression was down-regulated after injection of a miR-281 mimic while up-regulated after injection of a miR-281 inhibitor. miR-281 expression was suppressed by 20-hydroxyecdysone treatments but not affected by juvenile hormone treatments. Based on these findings, we propose that miR-281 participates in B. mori developmental regulation in the Malpighian tubules through suppression of BmEcR-B expression.

Involvement of microRNAs in infection of silkworm with bombyx mori cytoplasmic polyhedrosis virus (BmCPV)

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) is one of the most important pathogens of silkworm. MicroRNAs (miRNAs) have been demonstrated to play key roles in regulating host-pathogen interaction. However, there are limited reports on the miRNAs expression profiles during insect pathogen challenges. In this study, four small RNA libraries from BmCPV-infected midgut of silkworm at 72 h post-inoculation and 96 h post-inoculation and their corresponding control midguts were constructed and deep sequenced. A total of 316 known miRNAs (including miRNA*) and 90 novel miRNAs were identified. Fifty-eight miRNAs displayed significant differential expression between the infected and normal midgut (P value < = 0.01 and fold change > = 2.0 or < = 0.5), among which ten differentially expressed miRNA were validated by qRT-PCR method. Further bioinformatics analysis of predicted target genes of differentially expressed miRNAs showed that the miRNA targets were involved in stimulus and immune system process in silkworm.

Identification of microRNAs from Plutella xylostella larvae associated with parasitization by Diadegma semiclausum

MicroRNAs (miRNAs) as small non-coding RNAs play important roles in many biological processes such as development, cell signaling and immune response. Small RNA deep sequencing technology provided an opportunity for a thorough survey of miRNAs in a global key pest Plutella xylostella as well as comparative analysis of miRNA expression profile of the insect in association with parasitization by Diadegma semiclausum. Combining the deep sequencing data and bioinformatics, 235 miRNAs were identified from P. xylostella. Differential expression of host cellular miRNAs in response to parasitism was examined by making small RNA libraries from parasitized and naive second instar larvae of P. xylostella. Bantam, miR-276*, miR-10, miR-31 and miR-184 were detected as five most abundant miRNAs in both libraries and 96 miRNAs were identified that were differentially expressed after parasitization. Bantam*, miR-184 and miR-281* were significantly down-regulated and two miRNAs miR-279b and miR-2944b* were highly induced in parasitized larvae. Interestingly, high copy numbers and differential expression of several miRNA passenger strands (miRNA*) suggest their potential roles in host-parasitoid interaction. In conclusion, expression profiling of miRNAs provided insights into their possible involvement in insect immune response to parasitism and offer an important resource for further studies.

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.

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.

MicroRNAs in metamorphic and non-metamorphic transitions in hemimetabolan insect metamorphosis

Background

Previous work showed that miRNAs play key roles in the regulation of metamorphosis in the hemimetabolan species Blattella germanica. To gain insight about which miRNAs might be important, we have constructed two miRNA libraries, one of the penultimate, pre-metamorphic nymphal instar (N5) and the other of the last, metamorphic nymphal instar (N6).

Results

High throughput sequencing gave 61 canonical miRNAs present in the N5 and N6 libraries, although at different proportions in each. Comparison of both libraries led to the identification of three and 37 miRNAs significantly more expressed in N5 and N6 respectively. Twelve of these 40 miRNAs were then investigated further by qRT-PCR and results indicated that miR-252-3p was well expressed in N5 but not in N6, whereas let-7-5p, miR-100-5p and miR-125-5p showed the reverse pattern. 20-Hydroxyecdysone (20E) tended to stimulate miRNA expression, whereas juvenile hormone (JH) inhibited the 20E stimulatory effect. Expression of let-7, miR-100 and miR-125 was increased by 20E, which has also been observed in D. melanogaster. The only miRNA that was inhibited by 20E was miR-252-3p. The involvement of let-7, miR-100 and miR-125 in metamorphosis has been demonstrated in other insects. Depletion of miR-252-3p caused growth and developmental delays, which suggests that this miRNA is involved in regulating these processes prior to metamorphosis.

Conclusions

The comparative analysis of miRNA libraries from pre-metamorphic (N5) and metamorphic stages (N6) of B. germanica proved to be a useful tool to identify miRNAs with roles in hemimetabolan metamorphosis. Three miRNAs emerged as important factors in the metamorphic stage (N6): let-7-5p, miR-100-5p and miR-125-5p, whereas miR-252-3p appears to be important in the pre-metamorphic stage (N5).

Expression analysis of miRNAs in BmN cells

MicroRNAs (miRNAs) are the family of noncoding single-strand RNA molecules of 21-25 nucleotides in length and play a broad and key regulation role in various physiological and pathological processes including differentiation, apoptosis, proliferation, and tumorigenesis. In Bombyx mori, a total of 487 pre-miRNAs and 562 mature miRNAs were identified by experimental or computational approaches, but their functions remain unknown. To carry out the research of gain-of-function of miRNAs in BmN cells, we firstly identified the endogenous expression of miRNAs in BmN cells by microarray and found that only 73 miRNAs could be detected by miRNA microarray. Then three low abundance or undetected miRNAs, pri-mir-1a, pri-mir-8 and pri-mir-133, were selected to express in BmN cells. The eukaryotic expression vector pIEx-1 harboring baculovirus ie1 promoter and hr5 enhancer was screened and used for expressing miRNA in BmN cells. Three miRNA expression vectors pIEx-1-EGFP-pri-mir-1a/8/133 were constructed, which contained the three corresponding pri-miRNA sequences, respectively. The constructed miRNA vectors were successfully transfected into BmN cells and the qRT-PCR analysis showed that relative abundance of bmo-mir-1a, bmo-mir-8 and bmo-mir-133 in BmN cells transfected with the pIEx-1-EGFP-pri-mir-1a/8/133 is as 32, 4.4 and 904 times as that in BmN cells transfected with the control vector pIEx-1-EGFP, respectively. The present work lays a foundation for the further functional studies of miRNAs in silkworm.

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.

Functional shifts in insect microRNA evolution

MicroRNAs (miRNAs) are short endogenous RNA molecules that regulate gene expression at the posttranscriptional level and have been shown to play critical roles during animal development. The identification and comparison of miRNAs in metazoan species are therefore paramount for our understanding of the evolution of body plans. We have characterized 203 miRNAs from the red flour beetle Tribolium castaneum by deep sequencing of small RNA libraries. We can conclude, from a single study, that the Tribolium miRNA set is at least 15% larger than that in the model insect Drosophila melanogaster (despite tens of high-throughput sequencing experiments in the latter). The rate of birth and death of miRNAs is high in insects. Only one-third of the Tribolium miRNA sequences are conserved in D. melanogaster, and at least 18 Tribolium miRNAs are conserved in vertebrates but lost in Drosophila. More than one-fifth of miRNAs that are conserved between Tribolium and Drosophila exhibit changes in the transcription, genomic organization, and processing patterns that lead to predicted functional shifts. For example, 13% of conserved miRNAs exhibit seed shifting, and we describe arm-switching events in 11% of orthologous pairs. These shifts fundamentally change the predicted targets and therefore function of orthologous miRNAs. In general, Tribolium miRNAs are more representative of the insect ancestor than Drosophila miRNAs and are more conserved in vertebrates.

MicroRNAs of Bombyx mori identified by Solexa sequencing

BACKGROUND

MicroRNA (miRNA) and other small regulatory RNAs contribute to the modulation of a large number of cellular processes. We sequenced three small RNA libraries prepared from the whole body, and the anterior-middle and posterior silk glands of Bombyx mori, with a view to expanding the repertoire of silkworm miRNAs and exploring transcriptional differences in miRNAs between segments of the silk gland.

RESULTS

With the aid of large-scale Solexa sequencing technology, we validated 257 unique miRNA genes, including 202 novel and 55 previously reported genes, corresponding to 324 loci in the silkworm genome. Over 30 known silkworm miRNAs were further corrected in their sequence constitutes and length. A number of reads originated from the loop regions of the precursors of two previously reported miRNAs (bmo-miR-1920 and miR-1921). Interestingly, the majority of the newly identified miRNAs were silkworm-specific, 23 unique miRNAs were widely conserved from invertebrates to vertebrates, 13 unique miRNAs were limited to invertebrates, and 32 were confined to insects. We identified 24 closely positioned clusters and 45 paralogs of miRNAs in the silkworm genome. However, sequence tags showed that paralogs or clusters were not prerequisites for coordinated transcription and accumulation. The majority of silkworm-specific miRNAs were located in transposable elements, and displayed significant differences in abundance between the anterior-middle and posterior silk gland.

CONCLUSIONS

Conservative analysis revealed that miRNAs can serve as phylogenetic markers and function in evolutionary signaling. The newly identified miRNAs greatly enrich the repertoire of insect miRNAs, and provide insights into miRNA evolution, biogenesis, and expression in insects. The differential expression of miRNAs in the anterior-middle and posterior silk glands supports their involvement as new levels in the regulation of the silkworm silk gland.