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

Non-coding RNAs identification and regulatory networks in pathogen-host interaction in the microsporidia congenital infection

BACKGROUND

The interaction networks between coding and non-coding RNAs (ncRNAs) including long non-coding RNA (lncRNA), covalently closed circular RNA (circRNA) and miRNA are significant to elucidate molecular processes of biological activities and interactions between host and pathogen. Congenital infection caused by vertical transmission of microsporidia N. bombycis can result in severe economic losses in the silkworm-feeding industry. However, little is known about ncRNAs that take place in the microsporidia congenital infection. Here we conducted whole-transcriptome RNA-Seq analyses to identify ncRNAs and regulatory networks for both N. bombycis and host including silkworm embryos and larvae during the microsporidia congenital infection.

RESULTS

A total of 4,171 mRNAs, 403 lncRNA, 62 circRNAs, and 284 miRNAs encoded by N. bombycis were identified, among which some differentially expressed genes formed cross-talk and are involved in N. bombycis proliferation and infection. For instance, a lncRNA/circRNA competing endogenous RNA (ceRNA) network including 18 lncRNAs, one circRNA, and 20 miRNAs was constructed to describe 14 key parasites genes regulation, such as polar tube protein 3 (PTP3), ricin-B-lectin, spore wall protein 4 (SWP4), and heat shock protein 90 (HSP90). Regarding host silkworm upon N. bombycis congenital infection, a total of 14,889 mRNAs, 3,038 lncRNAs, 19,039 circRNAs, and 3,413 miRNAs were predicted based on silkworm genome with many differentially expressed coding and non-coding genes during distinct developmental stages. Different species of RNAs form interacting network to modulate silkworm biological processes, such as growth, metamorphosis and immune responses. Furthermore, a lncRNA/circRNA ceRNA network consisting of 140 lncRNAs, five circRNA, and seven miRNAs are constructed hypothetically to describe eight key host genes regulation, such as Toll-6, Serpin-6, inducible nitric oxide synthase (iNOS) and Caspase-8. Notably, cross-species analyses indicate that parasite and host miRNAs play a vital role in pathogen-host interaction in the microsporidia congenital infection.

CONCLUSION

This is the first comprehensive pan-transcriptome study inclusive of both N. bombycis and its host silkworm with a specific focus on the microsporidia congenital infection, and show that ncRNA-mediated regulation plays a vital role in the microsporidia congenital infection, which provides a new insight into understanding the basic biology of microsporidia and pathogen-host interaction.

High Daytime Temperature Responsive MicroRNA Profiles in Developing Grains of Rice Varieties with Contrasting Chalkiness

High temperature impairs starch biosynthesis in developing rice grains and thereby increases chalkiness, affecting the grain quality. Genome encoded microRNAs (miRNAs) fine-tune target transcript abundances in a spatio-temporal specific manner, and this mode of gene regulation is critical for a myriad of developmental processes as well as stress responses. However, the role of miRNAs in maintaining rice grain quality/chalkiness during high daytime temperature (HDT) stress is relatively unknown. To uncover the role of miRNAs in this process, we used five contrasting rice genotypes (low chalky lines Cyp, Ben, and KB and high chalky lines LaGrue and NB) and compared the miRNA profiles in the R6 stage caryopsis samples from plants subjected to prolonged HDT (from the onset of fertilization through R6 stage of caryopsis development). Our small RNA analysis has identified approximately 744 miRNAs that can be grouped into 291 families. Of these, 186 miRNAs belonging to 103 families are differentially regulated under HDT. Only two miRNAs, Osa-miR444f and Osa-miR1866-5p, were upregulated in all genotypes, implying that the regulations greatly varied between the genotypes. Furthermore, not even a single miRNA was commonly up/down regulated specifically in the three tolerant genotypes. However, three miRNAs (Osa-miR1866-3p, Osa-miR5150-3p and canH-miR9774a,b-3p) were commonly upregulated and onemiRNA (Osa-miR393b-5p) was commonly downregulated specifically in the sensitive genotypes (LaGrue and NB). These observations suggest that few similarities exist within the low chalky or high chalky genotypes, possibly due to high genetic variation. Among the five genotypes used, Cypress and LaGrue are genetically closely related, but exhibit contrasting chalkiness under HDT, and thus, a comparison between them is most relevant. This comparison revealed a general tendency for Cypress to display miRNA regulations that could decrease chalkiness under HDT compared with LaGrue. This study suggests that miRNAs could play an important role in maintaining grain quality in HDT-stressed rice.

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.

Exploring miRNA-mRNA regulatory modules responding to tannic acid stress in Micromelalopha troglodyta (Graeser) (Lepidoptera: Notodontidae) via small RNA sequencing

MicroRNAs (miRNAs) are small noncoding RNAs (sRNAs) that regulate gene expression by inhibiting translation or degrading mRNA. Although the functions of miRNAs in many biological processes have been reported, there is currently no research on the possible roles of miRNAs in Micromelalopha troglodyta (Graeser) involved in the response of plant allelochemicals. In this article, six sRNA libraries (three treated with tanic acid and three control) from M. troglodyta were constructed using Illumina sequencing. From the results, 312 known and 43 novel miRNAs were differentially expressed. Notably, some of the most abundant miRNAs, such as miR-432, miR-541-3p, and miR-4448, involved in important physiological processes were also identified. To better understand the function of the targeted genes, we performed Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The results indicated that differentially expressed miRNA targets were involved in metabolism, development, hormone biosynthesis, and immunity. Finally, we visualized a miRNA-mRNA regulatory module that supports the role of miRNAs in host-allelochemical interactions. To our knowledge, this is the first report on miRNAs responding to tannic acid in M. troglodyta. This study provides indispensable information for understanding the potential roles of miRNAs in M. troglodyta and the applications of these miRNAs in M. troglodyta management.

RNA sequencing identifies an ovary-enriched microRNA, miR-311-3p, involved in ovarian development and fecundity by targeting Endophilin B1 in Bactrocera dorsalis

BACKGROUND

The oriental fruit fly, Bactrocera dorsalis, is a highly invasive pest in East Asia and the Pacific. With the development of pesticides resistance, environment-friendly pesticides are urgently needed. MicroRNAs (miRNAs) are critical regulators of numerous biological processes, including reproduction. Thus, it is significant to identify reproductive-related miRNAs in this notorious pest to facilitate its control, such as RNAi-based biopesticides targeting essential miRNAs.

RESULTS

A high-throughput sequencing was carried out to identify miRNAs involved in reproduction from the ovary and fat body at four developmental stages [1 day (d), 5, 9, and 13 days post-eclosion] in female B. dorsalis. Results showed that 98 and 74 miRNAs were differentially expressed in ovary and fat body, respectively, during sexual maturation. Gene ontology analysis showed that target genes involved in oogenesis and lipid particle accounted for 33% and 15% of the total targets, respectively. Among these differentially expressed miRNAs, we found by qPCR that miR-311-3p was enriched in the ovary and down-regulated during sexual maturation. Injection of agomir-miR-311-3p resulted in arrested ovarian development, reduced egg deposition and progeny viability. Endophilin B1 was confirmed to be the target of miR-311-3p, via dual-luciferase assay and expression profiling. Knockdown of Endophilin B1 resulted in reproductive defects similar to those caused by injection of miR-311-3p agomir. Thus, miR-311-3p might play a critical role in female reproduction by targeting Endophilin B1.

CONCLUSION

Our data not only provides knowledge on the abundance of reproductive-related miRNAs and target genes, but also promotes new control strategies for this pest. © 2022 Society of Chemical Industry.

A Genetic Compensation Phenomenon and Global Gene Expression Changes in Sex-miR-2766-3p Knockout Strain of Spodoptera exigua Hübner (Lepidoptera: Noctuidae)

MicroRNAs (miRNAs) drive the post-transcriptional repression of target mRNAs and play important roles in a variety of biological processes. miR-2766-3p is conserved and abundant in Lepidopteran species and may be involved in a variety of biological activities. In this study, Sex-miR-2766-3p was predicted to potentially bind to the 3' untranslated region (UTR) of cap 'n' collar isoform C (CncC) in Spodoptera exigua, and Sex-miR-2766-3p was confirmed to regulate the expression of SeCncC through screening with a luciferase reporter system. Although CRISPR/Cas9 has been extensively utilized to examine insect gene function, studies of miRNA function are still relatively uncommon. Thus, we employed CRISPR/Cas9 to knock out Sex-miR-2766-3p from S. exigua. However, the expression of SeCncC was not significantly altered in the knockout strain (2766-KO) compared with that of the WHS strain. This result suggested that a miRNA knockout might lack phenotypes because of genetic robustness. Additionally, we used transcriptome analysis to examine how the global gene expression patterns of the Sex-miR-2766-3p knockout strain varied. RNA-seq data revealed 1746 upregulated and 2183 downregulated differentially expressed genes (DEGs) in the 2766-KO strain, which might be the result of Sex-miR-2766-3p loss or DNA lesions as the trigger for transcriptional adaptation. GO function classification and KEGG pathway analyses showed that these DEGs were enriched for terms related to binding, catalytic activity, metabolic process, and signal transduction. Our findings demonstrated that S. exigua could compensate for the missing Sex-miR-2766-3p by maintaining the expression of SeCncC by other pathways.

RNA Sequencing Reveals the Potential Adaptation Mechanism to Different Hosts of Grapholita molesta

Grapholita molesta is an important fruit tree worldwide pest which feeds on hosts extensively and does serious harm. In this paper, the growth and development parameters and protease activities of G. molesta fed on different hosts were compared. Using Illumina RNA sequencing technology, 18 midgut samples from five different hosts (apple, pear, plum, peach and peach shoots) and artificial diet were sequenced and compared with the reference genome, resulting in 15269 genes and 2785 predicted new genes. From 15 comparative combinations, DEGs were found from 286 to 4187 in each group, with up-regulated genes from 107 to 2395 and down-regulated genes from 83 to 2665. KEGG pathway analysis showed that DEGs were associated with amino acid metabolism, starch and sucrose metabolism, carbohydrate metabolism, and hydrolase activity. A total of 31 co-expression gene modules of different hosts were identified by WGCNA. qRT-PCR showed that the expression pattern of the trypsin gene was consistent with RNA sequencing. In this study, growth and development parameters, protease activity, DEGs, enrichment analysis and qRT-PCR were combined to reveal the adaptation process to different hosts of G. molesta in many aspects. The results of this study provide a basis for further exploration of the molecular mechanism of host adaptation of G. molesta.

microRNA profile of Hermetia illucens (black soldier fly) and its implications on mass rearing

The growing demands on protein producers and the dwindling available resources have made Hermetia illucens (the black soldier fly, BSF) an economically important species. Insights into the genome of this insect will better allow for robust breeding protocols, and more efficient production to be used as a replacement of animal feed protein. The use of microRNA as a method to understand how gene regulation allows insect species to adapt to changes in their environment, has been established in multiple species. The baseline and life stage expression levels established in this study, allow for insight into the development and sex-linked microRNA regulation in BSF. To accomplish this, microRNA was extracted and sequenced from 15 different libraries with each life stage in triplicate. Of the total 192 microRNAs found, 168 were orthologous to known arthropod microRNAs and 24 microRNAs were unique to BSF. Twenty-six of the 168 microRNAs conserved across arthropods had a statistically significant (p < 0.05) differential expression between Egg to Larval stages. The development from larva to pupa was characterized by 16 statistically significant differentially expressed microRNA. Seven and 9 microRNA were detected as statistically significant between pupa to adult female and pupa to adult male, respectively. All life stages had a nearly equal split between up and down regulated microRNAs. Ten of the unique 24 miRNA were detected exclusively in one life stage. The egg life stage expressed five microRNA (hil-miR-m, hil-miR-p, hil-miR-r, hil-miR-s, and hil-miR-u) not seen in any other life stages. The female adult and pupa life stages expressed one miRNA each hil-miR-h and hil-miR-ac respectively. Both male and female adult life stages expressed hil-miR-a, hil-miR-b, and hil-miR-y. There were no unique microRNAs found only in the larva stage. Twenty-two microRNAs with 56 experimentally validated target genes in the closely related Drosophila melanogaster were identified. Thus, the microRNA found display the unique evolution of BSF, along with the life stages and potential genes to target for robust mass rearing. Understanding of the microRNA expression in BSF will further their use in the crucial search for alternative and sustainable protein sources.

Identification of differentially expressed microRNAs under imidacloprid exposure in Sitobion miscanthi

Imidacloprid is a neonicotinoid that targets sucking pests, such as aphids and the green leaf bug and has been widely applied in wheat fields to control wheat aphids in China. To investigate the involvement of miRNAs in imidacloprid resistance, we sequenced small RNA libraries of Sitobion miscanthi Fabricius, across two different treatments using Illumina short-read sequencing technology. As a result, 265 microRNAs (miRNAs), of which 242 were known and 23 were novel, were identified. Quantitative analysis of miRNA levels showed that 23 miRNAs were significantly up-regulated, and 54 miRNAs were significantly down-regulated in the nymphs of S. miscanthi treated with imidacloprid in comparison with those of the control. Modulation of the abundances of differentially expressed miRNAs, smi-miR-316, smi-miR-1000, and smi-miR-iab-4 by the addition of the corresponding antagomir/inhibitor to the artificial diet significantly changed the susceptibility of S. miscanthi to imidacloprid. Subsequently, the post-transcriptional regulatory mechanism was conducted, smi-miR-278 and smi-miR-316 were confirmed to be participated in the post-transcriptional regulation of nAChRα1A and CYP4CJ6, respectively. The results suggested that miRNAs differentially expressed in response to imidacloprid could play a critical regulatory role in the metabolism of S. miscanthi to imidacloprid.

miRNAs of Aedes aegypti (Linnaeus 1762) conserved in six orders of the class Insecta

Aedes aegypti L. is the most important vector of arboviruses such as dengue, Zika, chikungunya, Mayaro, and yellow fever, which impact millions of people's health per year. MicroRNA profile has been described in some mosquito species as being important for biological processes such as digestion of blood, oviposition, sexual differentiation, insecticide resistance, and pathogens dissemination. We identified the miRNAs of Ae. aegypti females, males and eggs of a reference insecticide susceptible strain New Orleans and compared them with those other insects to determine miRNA fingerprint by new-generation sequencing. The sequences were analyzed using data mining tools and categorization, followed by differential expression analysis and conservation with other insects. A total of 55 conserved miRNAs were identified, of which 34 were of holometabolous insects and 21 shared with hemimetabolous insects. Of these miRNAs, 32 had differential expression within the stages analyzed. Three predominant functions of miRNA were related to embryonic development regulation, metamorphosis, and basal functions. The findings of this research describe new information on Ae. aegypti physiology which could be useful for the development of new control strategies, particularly in mosquito development and metamorphosis processes.

Multiple microRNAs control ecdysone signaling in the midgut of Spodoptera litura

Metamorphosis is one of the most important physiological processes in insects. It is regulated by a serial of ecdysone cascade genes. Recently, lots of microRNAs (miRNAs) were investigated in insects; however, their function in metamorphosis is largely unknown. In the present study, the dynamics of a small RNA population was investigated by RNA sequencing from the midgut of a lepidopteran pest Spodoptera litura during larval-pupal metamorphosis. A total of 101 miRNAs were identified, and 75 miRNAs were differentially expressed during the metamorphic process. The relationship between these differentially expressed miRNAs and 12 ecdysone cascade genes was analyzed by four classical software programs, and a multiple-to-multiple regulatory network was found to exist between these miRNAs and their targets. Among them, miR-14-3p and its two targets (EcR and E75) were chosen for further validation. MiR-14-3p had higher expression level in the 6th instar larvae as compared with either the prepupae or pupae, which was opposite to that of both EcR and E75, two ecdysone cascade genes. Luciferase reporter assay confirmed that both EcR and E75 were regulated by miR-14-3p. Interestingly, the 3' untranslated regions are nearly identical to each other among different transcript variants of the ecdysone cascade genes, including EcR, USP, E75, E74, E78, E93, Hr3, Hr4, Hr39, Krh1 and Ftzf1. Thus, different transcript variants of one ecdysone cascade gene could be regulated by the same miRNA. The above data suggest that the ecdysone signaling pathway is under the tight control of miRNA. These findings expand our understanding of the mechanism of insect metamorphosis and may also provide a novel possibility for the control of pest insects in the future.

High-resolution analysis of differential gene expression during skeletal muscle atrophy and programmed cell death

Skeletal muscles can undergo atrophy and/or programmed cell death (PCD) during development or in response to a wide range of insults, including immobility, cachexia, and spinal cord injury. However, the protracted nature of atrophy and the presence of multiple cell types within the tissue complicate molecular analyses. One model that does not suffer from these limitations is the intersegmental muscle (ISM) of the tobacco hawkmoth Manduca sexta. Three days before the adult eclosion (emergence) at the end of metamorphosis, the ISMs initiate a nonpathological program of atrophy that results in a 40% loss of mass. The ISMs then generate the eclosion behavior and initiate a nonapoptotic PCD during the next 30 h. We have performed a comprehensive transcriptomics analysis of all mRNAs and microRNAs throughout ISM development to better understand the molecular mechanisms that mediate atrophy and death. Atrophy involves enhanced protein catabolism and reduced expression of the genes involved in respiration, adhesion, and the contractile apparatus. In contrast, PCD involves the induction of numerous proteases, DNA methylases, membrane transporters, ribosomes, and anaerobic metabolism. These changes in gene expression are largely repressed when insects are injected with the insect steroid hormone 20-hydroxyecdysone, which delays death. The expression of the death-associated proteins may be greatly enhanced by reductions in specific microRNAs that function to repress translation. This study not only provides fundamental new insights into basic developmental processes, it may also represent a powerful resource for identifying potential diagnostic markers and molecular targets for therapeutic intervention.

Post-transcriptional regulation of insect metamorphosis and oogenesis

Metamorphic transformation from larvae to adults along with the high fecundity is key to insect success. Insect metamorphosis and reproduction are governed by two critical endocrines, juvenile hormone (JH), and 20-hydroxyecdysone (20E). Recent studies have established a crucial role of microRNA (miRNA) in insect metamorphosis and oogenesis. While miRNAs target genes involved in JH and 20E-signaling pathways, these two hormones reciprocally regulate miRNA expression, forming regulatory loops of miRNA with JH and 20E-signaling cascades. Insect metamorphosis and oogenesis rely on the coordination of hormones, cognate genes, and miRNAs for precise regulation. In addition, the alternative splicing of genes in JH and 20E-signaling pathways has distinct functions in insect metamorphosis and oogenesis. We, therefore, focus in this review on recent advances in post-transcriptional regulation, with the emphasis on the regulatory role of miRNA and alternative splicing, in insect metamorphosis and oogenesis. We will highlight important new findings of miRNA interactions with hormonal signaling and alternative splicing of JH receptor heterodimer gene Taiman.

Identification and Integrated Analysis of MicroRNA and mRNA Expression Profiles During Agonistic Behavior in Chinese Mitten Crab (Eriocheir sinensis) Using a Deep Sequencing Approach

As a commercially important species, the Chinese mitten crab (Eriocheir sinensis) has been cultured for a long time in China. Agonistic behavior often causes limb disability and requires much energy, which is harmful to the growth and survival of crabs. In this paper, we divided crabs into a control group (control, no treatment) and an experimental group (fight, agonistic behavior after 1 h) and then collected the thoracic ganglia (TG) to extract RNA. Subsequently, we first used a deep sequencing approach to examine the transcripts of microRNAs (miRNAs) and messenger RNAs (mRNAs) in E. sinensis displaying agonistic behavior. According to the results, we found 29 significant differentially expressed miRNAs (DEMs) and 116 significant differentially expressed unigenes (DEGs). The DEMs esi-miR-199a-5p, esi-let-7d, esi-miR-200a, and esi-miR-200b might participate in the regulation of agonistic behavior by mediating neuroregulation and energy metabolism. Focusing on the transcripts of the mRNAs, the renin–angiotensin system (RAS) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway might be involved in the regulation of agonistic behavior through glucose metabolism as this pathway was significantly enriched with DEGs. Besides, an integrated analysis of the miRNA and mRNA profiles revealed that the retinoid X receptor (RXR) was also involved in visual signal transduction, which was important for agonistic behavior. In addition, four vital agonistic behavior-related metabolic pathways, including the cAMP signaling, MAPK, protein digestion and absorption, and fatty acid metabolism pathways, were significantly enriched with the predicted target unigenes. In conclusion, the findings of this study might provide important insight enhancing our understanding of the underlying molecular mechanisms of agonistic behavior in E. sinensis.

Identification and characterization of miRNAs in an endoparasitoid wasp, Pteromalus puparum

MicroRNAs (miRNAs) are a form of endogenous small noncoding RNAs that regulate protein-coding gene expression at the posttranscriptional level. So far, knowledge of miRNAs in parasitoids remains rudimentary. We investigated miRNAs in Pteromalus puparum, a pupal endoparasitoid wasp with genome and transcriptome sequences completed. In this study, we constructed eight small RNA libraries from selected developmental stages and genders: male embryos, male larvae, male pupae, male adults, mixed-sex embryos, mixed-sex larvae, mixed-sex pupae, and female adults. We identified 254 mature miRNAs with 5p/3p arm features originated from 75 known and 119 novel miRNA genes in P. puparum, 88 of which reside in 26 clusters. The miRNAs in more than half of the clusters exhibit a consistent expression pattern, indicating they were co-transcribed from a long transcript. Comparing miRNA expression in the eight libraries, we found that 84 mature miRNAs were differentially expressed between embryos and larvae, 20 between larvae and pupae, and 26 between pupae and adults. We found some miRNAs were differentially expressed between sexes in embryos (10), larvae (29), pupae (8), and adults (14). Target predictions resulted in 211,571 miRNA-mRNA interactions for 254 different mature miRNAs. These miRNAs may be involved in sexual and developmental regulation of gene expression.

The Role of MicroRNAs in Genome Response to Plant-Lepidoptera Interaction

RNA interference is a known phenomenon of plant immune responses, involving the regulation of gene expression. The key components triggering the silencing of targeted sequences are double-stranded RNA molecules. The regulation of host-pathogen interactions is controlled by miRNA molecules, which regulate the expression of host resistance genes or the genes of the pathogen. The review focused on basic principles of RNA interference as a gene-silencing-based defense mechanism and the role of miRNA molecules in insect genomes. RNA interference as a tool for plant protection management is discussed. The review summarizes current miRNA-based biotechnology approaches for plant protection management.

RNA-sequencing of a citrus bud-feeder, Podagricomela weisei (Coleoptera: Chrysomelidae), reveals xenobiotic metabolism/core RNAi machinery-associated genes and conserved miRNAs

The citrus leaf-mining beetle, Podagricomela weisei Heikertinger, is an important citrus pest that ingests the mesophyll and new shoots. The mechanism underlying the xenobiotic metabolism of P. weisei is not well understood, in part because of a lack of available genomic and transcriptomic data, which has hampered the development of novel pest management approaches [e.g., RNA interference (RNAi)]. In this study, we completed the deep sequencing of the P. weisei transcriptome to identify factors potentially involved in xenobiotic metabolism and the core RNAi machinery. The sequencing of the P. weisei transcriptome generated >27 million clean reads, ultimately yielding 90,410 unigenes with an N50 of 1065 bp. The unigenes were used as queries to search the Nr database, which revealed that 21,847 unigenes were homologous to known genes in various species. Transcripts encoding genes involved in xenobiotic metabolism were identified, including genes encoding cytochrome P450 monooxygenase (P450, 47 unigenes), glutathione S-transferase (GST, 12 unigenes), esterase (EST, 25 unigenes), and the ATP-binding cassette transporter (ABC transporter, 32 unigenes). A parallel sequencing of small RNAs detected 30 conserved miRNAs, with the most abundant being Pwe-miR-1-3p, with an expression level reaching 517,996 reads in the prepared library, followed by Pwe-miR-8-3p (149,402 reads). Genes encoding components of the miRNA, siRNA, and piRNA pathways were also identified, and the results indicated that P. weisei possesses only one of each gene in all three pathways. In summary, this is the first detailed analysis of the transcriptome and small RNAs of P. weisei. The datasets presented herein may form the basis for future molecular characterizations of P. weisei as well as the development of enhanced pest control strategies.

Genome-wide identification, expression profiling, and target gene analysis of microRNAs in the Onion thrips, Thrips tabaci Lindeman (Thysanoptera: Thripidae), vectors of tospoviruses (Bunyaviridae)

Thrips tabaci Lindeman is an important polyphagous insect pest species estimated to cause losses of more than U.S. $1 billion worldwide annually. Chemical insecticides are of limited use in the management of T. tabaci due to the thigmokinetic behavior and development of resistance to insecticides. There is an urgent need to find alternative management strategies. Small noncoding RNAs (sncRNAs) especially microRNAs (miRNAs) hold great promise as key regulators of gene expression in a wide range of organisms. MiRNAs are a group of endogenously originated sncRNA known to regulate gene expression in animals, plants, and protozoans. In this study, we explored these RNAs in T. tabaci using deep sequencing to provide a basis for future studies of their biological and physiological roles in governing gene expression. Apart from snoRNAs and piRNAs, our study identified nine novel and 130 known miRNAs from T. tabaci. Functional classification of the targets for these miRNAs predicted that majority are involved in regulating transcription, translation, signal transduction and genetic information processing. The higher expression of few miRNAs (such as tta-miR-281, tta-miR-184, tta-miR-3533, tta-miR-N1, tta-miR-N7, and tta-miR-N9) in T. tabaci pupal and adult stages reflected their possible role in larval and adult development, metamorphosis, parthenogenesis, and reproduction. This is the first exploration of the miRNAome in T. tabaci, which not only provides insights into their possible role in insect metamorphosis, growth, and development but also offer an important resource for future pest management strategies.

Genome-Wide Identification of Destruxin A-Responsive Immunity-Related MicroRNAs in Diamondback Moth, Plutella xylostella

Plutella xylostella, a global key pest, is one of the major lepidopteran pests of cruciferous vegetables owing to its strong ability of resistance development to a wide range of insecticides. Destruxin A, a mycotoxin of the entomopathogenic fungus, Metarhizium anisopliae, has broad-spectrum insecticidal effects and has been used as an alternative control strategy to reduce harmful effects of insecticides. However, microRNA (miRNA)-regulated reactions against destruxin A have not been elucidated yet. Therefore, here, to identify immunity-related miRNAs, we constructed four small RNA libraries from destruxin A-injected larvae of P. xylostella at three different time courses (2, 4, and 6 h) with a control, and sequenced by Illumina. Our results showed that totally 187 known and 44 novel miRNAs were identified in four libraries by bioinformatic analysis. Interestingly, among differentially expressed known miRNAs, some conserved miRNAs, such as miR-263, miR-279, miR-306, miR-2a, and miR-308, predicted to be involved in regulating immunity-related genes, were also identified. Worthy to mention, miR-306 and miR-279 were also listed as common abundantly expressed miRNA in all treatments. The Kyoto Encyclopedia of Genes and Genomes pathway analysis also indicated that differentially expressed miRNAs were involved in several immunity-related signaling pathways, including toll signaling pathway, IMD signaling pathway, JAK-STAT signaling pathway, and cell adhesion molecules signaling pathway. To the best of our knowledge, this is the first comprehensive report of destruxin A-responsive immunity-related miRNAs in P. xylostella. Our findings will improve in understanding the role of destruxin A-responsive miRNAs in the host immune system and would be useful to develop biological control strategies for controlling P. xylostella.

Genome-Wide Profiling of Plutella xylostella Immunity-Related miRNAs after Isaria fumosorosea Infection

The development of resistance by Plutella xylostella to almost all insecticides is of significant concern all over the world. Entomopathogenic fungi such as Isaria fumosorosea have been used as an alternative to insecticides. However, the knowledge of miRNA-regulated reactions against entomopathogenic fungi is still in its infant stage. In the present study, P. xylostella was challenged with I. fumosorosea at four different time points (12, 18, 24, and 36 h) including a control, to build miRNA libraries by Illumina sequencing. The results of differential expression analysis exhibited that 23 miRNAs were differentially expressed, compared to control, in all treatments. It is worth mentioning, of these, some conserved miRNAs such as miR-2, miR-9a, miR-745, miR-7b, and miR-2767, known to play critical roles in host-pathogen interaction, were also identified. Furthermore, differentially expressed miRNAs were validated by RT-qPCR. Our results provide an essential information for further functional studies of the interaction between I. fumosorosea and P. xylostella at the post-transcriptional level.

An analysis of 67 RNA-seq datasets from various tissues at different stages of a model insect, Manduca sexta

Background

Manduca sexta is a large lepidopteran insect widely used as a model to study biochemistry of insect physiological processes. As a part of its genome project, over 50 cDNA libraries have been analyzed to profile gene expression in different tissues and life stages. While the RNA-seq data were used to study genes related to cuticle structure, chitin metabolism and immunity, a vast amount of the information has not yet been mined for understanding the basic molecular biology of this model insect. In fact, the basic features of these data, such as composition of the RNA-seq reads and lists of library-correlated genes, are unclear. From an extended view of all insects, clear-cut tempospatial expression data are rarely seen in the largest group of animals including Drosophila and mosquitoes, mainly due to their small sizes.

Results

We obtained the transcriptome data, analyzed the raw reads in relation to the assembled genome, and generated heatmaps for clustered genes. Library characteristics (tissues, stages), number of mapped bases, and sequencing methods affected the observed percentages of genome transcription. While up to 40% of the reads were not mapped to the genome in the initial Cufflinks gene modeling, we identified the causes for the mapping failure and reduced the number of non-mappable reads to <8%. Similarities between libraries, measured based on library-correlated genes, clearly identified differences among tissues or life stages. We calculated gene expression levels, analyzed the most abundantly expressed genes in the libraries. Furthermore, we analyzed tissue-specific gene expression and identified 18 groups of genes with distinct expression patterns.

Conclusion

We performed a thorough analysis of the 67 RNA-seq datasets to characterize new genomic features of M. sexta. Integrated knowledge of gene functions and expression features will facilitate future functional studies in this biochemical model insect.

Electronic supplementary material

The online version of this article doi: (10.1186/s12864-017-4147-y) contains supplementary material, which is available to authorized users.

Comparative analysis of miRNA expression during the development of insects of different metamorphosis modes and germ-band types

Background

Do miRNAs contribute to specify the germ-band type and the body structure in the insect embryo? Our goal was to address that issue by studying the changes in miRNA expression along the ontogeny of the German cockroach Blattella germanica, which is a short germ-band and hemimetabolan species.

Results

We sequenced small RNA libraries representing 11 developmental stages of B. germanica ontogeny (with especial emphasis on embryogenesis) and the changes in miRNA expression were examined. Data were compared with equivalent data for two long germ-band holometabolan species Drosophila melanogaster and Drosophila virilis, and the short germ-band holometabolan species Tribolium castaneum. The identification of B. germanica embryo small RNA sequences unveiled miRNAs not detected in previous studies, such as those of the MIR-309 family and 54 novel miRNAs. Four main waves of miRNA expression were recognized (with most miRNA changes occurring during the embryonic stages): the first from day 0 to day 1 of embryogenesis, the second during mid-embryogenesis (days 0–6), the third (with an acute expression peak) on day 2 of embryonic development, and the fourth during post-embryonic development. The second wave defined the boundaries of maternal-to-zygotic transition, with maternal mRNAs being cleared, presumably by Mir-309 and associated scavenger miRNAs.

Conclusion

miRNAs follow well-defined patterns of expression over hemimetabolan ontogeny, patterns that are more diverse during embryonic development than during the nymphal stages. The results suggest that miRNAs play important roles in the developmental transitions between the embryonic stages of development (starting with maternal loading), during which they might influence the germ-band type and metamorphosis mode.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4177-5) contains supplementary material, which is available to authorized users.

Identification of microRNAs and their response to the stress of plant allelochemicals in Aphis gossypii (Hemiptera: Aphididae)

BACKGROUND

MicroRNAs (miRNAs) are a group of short non-coding RNAs involved in the inhibition of protein translation or in mRNA degradation. Although the regulatory roles of miRNAs in various biological processes have been investigated, there is as yet an absence of studies about the regulatory roles of miRNAs involved in the metabolism of plant allelochemicals in insects.

RESULTS

We constructed five small RNA libraries from apterous Aphis gossypii adults that had fed on an artificial diet containing various allelochemicals. Using Illumina sequencing, a total of 73.27 million clean reads was obtained, and 292 miRNAs were identified from A. gossypii. Comparative analysis of read counts indicated that both conserved and novel miRNAs were differently expressed among the five libraries, and the differential expression was validated via qRT-PCR. We found that the transcript levels of several miRNAs were increased or decreased in all of the allelochemical treatment libraries compared to the control. The putative target genes of the miRNAs were predicted using in silico tools, and the target genes of several miRNAs were presumed to be involved in the metabolism of xenobiotic compounds. Furthermore, the target prediction results were confirmed using dual luciferase reporter assay, and Ago-miR-656a-3p was demonstrated to regulate the expression of CYP6J1 post-transcriptionally through binding to the 3' UTR of CYP6J1.

CONCLUSION

Our research results indicate that miRNAs may be involved in the metabolism of plant allelochemicals in A. gossypii, and these results also represent an important new small RNA genomics resource for further studies on this topic.

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.

Identification and characterization of microRNAs in the screwworm flies Cochliomyia hominivorax and Cochliomyia macellaria (Diptera: Calliphoridae)

MicroRNAs (miRNAs) are small noncoding RNAs that modulate gene expression through post-transcriptional regulation. Here, we report the identification and characterization of miRNAs in two closely related screwworm flies with different feeding habits: Cochliomyia hominivorax and Cochliomyia macellaria. The New World screwworm, C. hominivorax, is an obligatory parasite of warm-blooded vertebrates, whereas the secondary screwworm, C. macellaria, is a free-living organism that feeds on decaying organic matter. Here, the small RNA transcriptomes of adults and third-instar larvae of both species were sequenced. A total of 110 evolutionarily conserved miRNAs were identified, and 10 putative precursor miRNAs (pre-miRNAs) were predicted. The relative expression of six selected miRNAs was further investigated, including miRNAs that are related to reproduction and neural processes in other insects. Mature miRNAs were also characterized across an evolutionary time scale, suggesting that the majority of them have been conserved since the emergence of the Arthropoda [540 million years ago (Ma)], Hexapoda (488 Ma) and Brachycera (195 Ma) lineages. This study is the first report of miRNAs for screwworm flies. We also performed a comparative analysis with the hereby predicted miRNAs from the sheep blowfly, Lucilia cuprina. The results presented may advance our understanding of parasitic habits within Calliphoridae and assist further functional studies in blowflies.

MicroRNA-34 directly targets pair-rule genes and cytoskeleton component in the honey bee

MicroRNAs (miRNAs) are key regulators of developmental processes, such as cell fate determination and differentiation. Previous studies showed Dicer knockdown in honeybee embryos disrupt the processing of functional mature miRNAs and impairs embryo patterning. Here we investigated the expression profiles of miRNAs in honeybee embryogenesis and the role of the highly conserved miR-34-5p in the regulation of genes involved in insect segmentation. A total of 221 miRNAs were expressed in honey bee embryogenesis among which 97 mature miRNA sequences have not been observed before. Interestingly, we observed a switch in dominance between the 5-prime and 3-prime arm of some miRNAs in different embryonic stages; however, most miRNAs present one dominant arm across all stages of embryogenesis. Our genome-wide analysis of putative miRNA-target networks and functional pathways indicates miR-34-5p is one of the most conserved and connected miRNAs associated with the regulation of genes involved in embryonic patterning and development. In addition, we experimentally validated that miR-34-5p directly interacts to regulatory elements in the 3'-untranslated regions of pair-rule (even-skipped, hairy, fushi-tarazu transcription factor 1) and cytoskeleton (actin5C) genes. Our study suggests that miR-34-5p may regulate the expression of pair-rule and cytoskeleton genes during early development and control insect segmentation.

Comparative profiling of microRNAs in the winged and wingless English grain aphid, Sitobion avenae (F.) (Homoptera: Aphididae)

MicroRNAs (miRNAs) are short single-stranded non-coding RNAs that regulate gene expression, particularly during development. In this study, 345 miRNAs were identified from the English green aphid, Sitobion avenae (F.), of which 168 were conserved and 177 were S. avenae-specific. Quantitative comparison of miRNA expression levels indicated that 16 and 12 miRNAs were significantly up-regulated in winged and wingless S. avenae small RNA libraries, respectively. Differential expression of these miRNAs was confirmed by real-time quantitative RT-PCR validation. The putative transcript targets for these candidate miRNAs were predicted based on sequences from a model species Drosophila melanogaster and four aphid species Acyrthosiphon pisum, Myzus persicae, Toxoptera citricida, and Aphis gosspii. Gene Ontology and KEGG pathway analyses shed light on the potential functions of these miRNAs in the regulation of genes involved in the metabolism, development and wing polyphenism of S. avenae.

The First Report of miRNAs from a Thysanopteran Insect, Thrips palmi Karny Using High-Throughput Sequencing

Thrips palmi Karny (Thysanoptera: Thripidae) is the sole vector of Watermelon bud necrosis tospovirus, where the crop loss has been estimated to be around USD 50 million annually. Chemical insecticides are of limited use in the management of T. palmi due to the thigmokinetic behaviour and development of high levels of resistance to insecticides. There is an urgent need to find out an effective futuristic management strategy, where the small RNAs especially microRNAs hold great promise as a key player in the growth and development. miRNAs are a class of short non-coding RNAs involved in regulation of gene expression either by mRNA cleavage or by translational repression. We identified and characterized a total of 77 miRNAs from T. palmi using high-throughput deep sequencing. Functional classifications of the targets for these miRNAs revealed that majority of them are involved in the regulation of transcription and translation, nucleotide binding and signal transduction. We have also validated few of these miRNAs employing stem-loop RT-PCR, qRT-PCR and Northern blot. The present study not only provides an in-depth understanding of the biological and physiological roles of miRNAs in governing gene expression but may also lead as an invaluable tool for the management of thysanopteran insects in the future.

Identification and Characterization of Sex-Biased MicroRNAs in Bactrocera dorsalis (Hendel)

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that regulate various biological processes including sexual dimorphism. The oriental fruit fly Bactrocera dorsalis is one of the most destructive agricultural insect pests in many Asian countries. However, no miRNAs have been identified from the separate sex and gonads to elucidate sex gonad differentiation in B. dorsalis. In this study, we constructed four small RNA libraries from whole body of females, males (except ovaries and testes) and ovaries, testes of B. dorsalis for deep sequencing. The data analysis revealed 183 known and 120 novel miRNAs from these libraries. 18 female-biased and 16 male-biased miRNAs that may be involved in sexual differentiation were found by comparing the miRNA expression profiles in the four libraries. Using a bioinformatic approach, we predicted doublesex (dsx) as a target gene of the female-biased miR-989-3p which is considered as the key switch gene in the sex determination of tephritid insects. This study reveals the first miRNA profile related to the sex differentiation and gives a first insight into sex differences in miRNA expression of B. dorsalis which could facilitate studies of the reproductive organ specific roles of miRNAs.

MicroRNA signatures characterizing caste-independent ovarian activity in queen and worker honeybees (Apis mellifera L.)

Queen and worker honeybees differ profoundly in reproductive capacity. The queen of this complex society, with 200 highly active ovarioles in each ovary, is the fertile caste, whereas the workers have approximately 20 ovarioles as a result of receiving a different diet during larval development. In a regular queenright colony, the workers have inactive ovaries and do not reproduce. However, if the queen is sensed to be absent, some of the workers activate their ovaries, producing viable haploid eggs that develop into males. Here, a deep-sequenced ovary transcriptome library of reproductive workers was used as supporting data to assess the dynamic expression of the regulatory molecules and microRNAs (miRNAs) of reproductive and nonreproductive honeybee females. In this library, most of the differentially expressed miRNAs are related to ovary physiology or oogenesis. When we quantified the dynamic expression of 19 miRNAs in the active and inactive worker ovaries and compared their expression in the ovaries of virgin and mated queens, we noted that some miRNAs (miR-1, miR-31a, miR-13b, miR-125, let-7 RNA, miR-100, miR-276, miR-12, miR-263a, miR-306, miR-317, miR-92a and miR-9a) could be used to identify reproductive and nonreproductive statuses independent of caste. Furthermore, integrative gene networks suggested that some candidate miRNAs function in the process of ovary activation in worker bees.

Deep sequencing of small RNA libraries from human prostate epithelial and stromal cells reveal distinct pattern of microRNAs primarily predicted to target growth factors

Complex epithelial and stromal cell interactions are required during the development and progression of prostate cancer. Regulatory small non-coding microRNAs (miRNAs) participate in the spatiotemporal regulation of messenger RNA (mRNA) and regulation of translation affecting a large number of genes involved in prostate carcinogenesis. In this study, through deep-sequencing of size fractionated small RNA libraries we profiled the miRNAs of prostate epithelial (PrEC) and stromal (PrSC) cells. Over 50 million reads were obtained for PrEC in which 860,468 were unique sequences. Similarly, nearly 76 million reads for PrSC were obtained in which over 1 million were unique reads. Expression of many miRNAs of broadly conserved and poorly conserved miRNA families were identified. Sixteen highly expressed miRNAs with significant change in expression in PrSC than PrEC were further analyzed in silico. ConsensusPathDB showed the target genes of these miRNAs were significantly involved in adherence junction, cell adhesion, EGRF, TGF-β and androgen signaling. Let-7 family of tumor-suppressor miRNAs expression was highly pervasive in both, PrEC and PrSC cells. In addition, we have also identified several miRNAs that are unique to PrEC or PrSC cells and their predicted putative targets are a group of transcription factors. This study provides perspective on the miRNA expression in PrEC and PrSC, and reveals a global trend in miRNA interactome. We conclude that the most abundant miRNAs are potential regulators of development and differentiation of the prostate gland by targeting a set of growth factors. Additionally, high level expression of the most members of let-7 family miRNAs suggests their role in the fine tuning of the growth and proliferation of prostate epithelial and stromal cells.

Functional screen for microRNAs of Nilaparvata lugens reveals that targeting of glutamine synthase by miR-4868b regulates fecundity

Insect fecundity is regulated by the interaction of genotypes and the environment. MicroRNAs (miRNAs) also act in insect development and reproduction by regulating genes involved in these physiological processes. Although hundreds of insect miRNAs have been identified, the biological roles of most remain poorly understood. Here, we used a multi-algorithm approach for miRNA target prediction in 3'UTRs of fecundity-related genes in the brown planthopper (BPH) Nilaparvata lugens and identified 38 putative miRNAs targeting 9 fecundity-related genes. High-ranked miRNAs were selected for target validation. Using a dual luciferase reporter assay in S2 cells, we experimentally verified N. lugens glutamine synthetase (NlGS) as an authentic target of microRNA-4868b (miR-4868b). In the females, NlGS protein expression was down-regulated after injection of a miR-4868b mimic but up-regulated after injection of a miR-4868b inhibitor. In addition, overexpression of miR-4868b reduced fecundity, and disrupted ovary development and Vg expression in N. lugens. These findings showed that miR-4868b is involved in regulating N. lugens fecundity by targeting NlGS. Moreover, this study may lead to better understanding of the fecundity of this important agricultural insect pest.

A Diversity of Conserved and Novel Ovarian MicroRNAs in the Speckled Wood (Pararge aegeria)

microRNAs (miRNAs) are important regulators of animal development and other processes, and impart robustness to living systems through post-transcriptional regulation of specific mRNA transcripts. It is postulated that newly emergent miRNAs are generally expressed at low levels and with spatiotemporally restricted expression domains, thus minimising effects of spurious targeting on animal transcriptomes. Here we present ovarian miRNA transcriptome data for two geographically distinct populations of the Speckled Wood butterfly (Pararge aegeria). A total of 74 miRNAs were identified, including 11 newly discovered and evolutionarily-young miRNAs, bringing the total of miRNA genes known from P. aegeria up to 150. We find a positive correlation between miRNA age and expression level. A common set of 55 miRNAs are expressed in both populations. From this set, we identify seven that are consistently either ovary-specific or highly upregulated in ovaries relative to other tissues. This ‘ovary set’ includes miRNAs with known contributions to ovarian function in other insect species with similar ovaries and mode of oogenesis, including miR-989 and miR-2763, plus new candidates for ovarian function. We also note that conserved miRNAs are overrepresented in the ovary relative to the whole body.

Identification of Novel and Conserved microRNAs in Homalodisca vitripennis, the Glassy-Winged Sharpshooter by Expression Profiling

The glassy-winged sharpshooter (GWSS) Homalodisca vitripennis (Hemiptera: Cicadellidae), is a xylem-feeding leafhopper and an important vector of the bacterium Xylella fastidiosa; the causal agent of Pierce’s disease of grapevines. MicroRNAs are a class of small RNAs that play an important role in the functional development of various organisms including insects. In H. vitripennis, we identified microRNAs using high-throughput deep sequencing of adults followed by computational and manual annotation. A total of 14 novel microRNAs that are not found in the miRBase were identified from adult H. vitripennis. Conserved microRNAs were also found in our datasets. By comparison to our previously determined transcriptome sequence of H. vitripennis, we identified the potential targets of the microRNAs in the transcriptome. This microRNA profile information not only provides a more nuanced understanding of the biological and physiological mechanisms that govern gene expression in H. vitripennis, but may also lead to the identification of novel mechanisms for biorationally designed management strategies through the use of microRNAs.

Transcriptomic response of Manduca sexta immune tissues to parasitization by the bracovirus associated wasp Cotesia congregata

During oviposition, Cotesia congregata parasitoid wasps inject into their host, Manduca sexta, some biological factors such as venom, ovarian fluid and a symbiotic polydnavirus (PDV) named Cotesia congregata bracovirus (CcBV). During parasitism, complex interactions occur between wasp-derived factors and host targets that lead to important modifications in host physiology. In particular, the immune response leading to wasp egg encapsulation is inhibited allowing wasp survival. To date, the regulation of host genes during the interaction had only been studied for a limited number of genes. In this study, we analysed the global impact of parasitism on host gene regulation 24 h post oviposition by high throughput 454 transcriptomic analyses of two tissues known to be involved in the host immune response (hemocytes and fat body). To identify specific effects of parasitism on host transcription at this time point, transcriptomes were obtained from non-treated and parasitized larvae, and also from larvae injected with heat-killed bacteria and double stimulated larvae that were parasitized prior to bacterial challenge. Results showed that, immune challenge by bacteria leads to induction of certain antimicrobial peptide (AMP) genes in M. sexta larvae whether they were parasitized or not prior to bacterial challenge. These results show that at 24 h post oviposition pathways leading to expression of AMP genes are not all inactivated suggesting wasps are in an antiseptic environment. In contrast, at this time point genes involved in phenoloxidase activation and cellular immune responses were globally down-regulated after parasitism in accordance with the observed inhibition of wasp egg encapsulation.

Identification and profiling of Manduca sexta microRNAs and their possible roles in regulating specific transcripts in fat body, hemocytes, and midgut

Significance of microRNA-mediated posttranscriptional regulation has been appreciated ever since its discovery. In the tobacco hornworm Manduca sexta, 164 conserved and 16 novel microRNAs have been identified experimentally (Zhang et al., 2012, 2014). To extend the list of microRNAs in this lepidopteran model species and further explore their possible regulatory roles, we constructed and sequenced small RNA libraries of M. sexta fat body, hemocytes and midgut, since transcriptomes of these tissues from the 5th instar larvae had been studied quite extensively. Each library represented a mixture of the same tissues from larvae that were naïve or induced by three different pathogens. From a total of 167 million reads obtained, we identified two new variants of conserved miR-281 and miR-305 and six novel microRNAs. Abundances of all microRNAs were normalized and compared to reveal their differential expression in these three tissues. Star strands of ten microRNAs were present at higher levels than the corresponding mature strands. From a list of tissue-specific transcripts, we predicted target sites in 3'-UTRs using preferentially expressed microRNA groups in each tissue and suggested possible regulatory roles of these microRNAs in energy metabolism, insecticide resistance, and some mitochondrial and immune gene expression. Examining manifold targets, microRNA regulations were suggested of multiple physiological processes. This study has enriched our knowledge of M. sexta microRNAs and how microRNAs potentially coordinate different physiological processes.

Identification and expression profiling of Helicoverpa armigera microRNAs and their possible role in the regulation of digestive protease genes

The present investigation is an effort to determine the possible roles of microRNAs (miRNAs) in the regulation of protease gene expression in Helicoverpa armigera upon exposure to plant protease inhibitors (PIs). Using Illumina platform, deep sequencing of 12 small RNA libraries was performed from H. armigera larvae fed on artificial diet (AD) or recombinant Capsicum annuum PI-7 (rCanPI-7) incorporated diet, at various time intervals (0.5, 2, 6, 12, 24, and 48 h). Sequencing data were analyzed with miRDeep2 software; a total of 186 unique miRNAs were identified from all the 12 libraries, out of which 96 were conserved while 90 were novel. These miRNAs showed all the conserved characteristics of insect miRNAs. Homology analysis revealed that most of the identified miRNAs were insect-specific, and more than 50 miRNAs were Lepidoptera-specific. Several candidate miRNAs (conserved and novel) were differentially expressed in rCanPI-7 fed larvae as compared to the larvae fed on AD. H. armigera miRNAs were found to have target sites in several protease genes as well as in protease regulation related genes such as serine PI and immune reactive PI. As expected, negative correlation in the relative abundance of miRNAs and their target mRNAs was evident from qualitative real time polymerase chain reaction analysis. The investigation revealed potential roles of miRNAs in H. armigera protease gene regulation.

Deep sequencing of small RNA libraries reveals dynamic expression patterns of microRNAs in multiple developmental stages of Bactrocera dorsalis

In eukaryotes, microRNAs (miRNAs) are small, conserved, noncoding RNAs that have emerged as critical regulators of gene expression. The oriental fruit fly Bactrocera dorsalis is one of the most economically important fruit fly pests in East Asia and the Pacific. Although transcriptome analyses have greatly enriched our knowledge of its structural genes, little is known about post-transcriptional regulation by miRNAs in this dipteran species. In this study, small RNA libraries corresponding to four B. dorsalis developmental stages (eggs, larvae, pupae and adults) were constructed and sequenced. Approximately 30.7 million reads of 18-30 nucleotides were obtained, with 123 known miRNAs and 60 novel miRNAs identified amongst these libraries. More than half of the miRNAs were stage-specific during the four developmental stages. A set of miRNAs was found to be up- or down-regulated during development by comparison of their reads at different developmental stages. Moreover, a small part of miRNAs owned both miR-#-3p and miR-#-5p types, with enormously variable miR-#-3p/miR-#-5p ratios in the same library and amongst different developmental stages for each miRNA. Taking these findings together, the current study has uncovered a number of miRNAs and provided insights into their possible involvement in developmental regulation by expression profiling of miRNAs. Further analyses of the expression and function of these miRNAs could increase our understanding of regulatory networks in this insect and lead to novel approaches for its control.

Sex-biased expression of microRNAs in Drosophila melanogaster

Most animals have separate sexes. The differential expression of gene products, in particular that of gene regulators, is underlying sexual dimorphism. Analyses of sex-biased expression have focused mostly on protein-coding genes. Several lines of evidence indicate that microRNAs, a class of major gene regulators, are likely to have a significant role in sexual dimorphism. This role has not been systematically explored so far. Here, I study the sex-biased expression pattern of microRNAs in the model species Drosophila melanogaster. As with protein-coding genes, sex-biased microRNAs are associated with the reproductive function. Strikingly, contrary to protein-coding genes, male-biased microRNAs are enriched in the X chromosome, whereas female microRNAs are mostly autosomal. I propose that the chromosomal distribution is a consequence of high rates of de novo emergence, and a preference for new microRNAs to be expressed in the testis. I also suggest that demasculinization of the X chromosome may not affect microRNAs. Interestingly, female-biased microRNAs are often encoded within protein-coding genes that are also expressed in females. MicroRNAs with sex-biased expression do not preferentially target sex-biased gene transcripts. These results strongly suggest that the sex-biased expression of microRNAs is mainly a consequence of high rates of microRNA emergence in the X chromosome (male bias) or hitchhiked expression by host genes (female bias).

Identification of conserved and novel microRNAs in Manduca sexta and their possible roles in the expression regulation of immunity-related genes

The tobacco hornworm Manduca sexta has served as a model for insect biochemical and physiological research for decades. However, knowledge of the posttranscriptional regulation of gene expression by microRNAs is still rudimentary in this species. Our previous study (Zhang et al., 2012) identified 163 conserved and 13 novel microRNAs in M. sexta, most of which were present at low levels in pupae. To identify additional M. sexta microRNAs and more importantly to examine their possible roles in the expression regulation of immunity-related genes, we constructed four small RNA libraries using fat body and hemocytes from naïve or bacteria-injected larvae and obtained 32.9 million reads of 18-31 nucleotides by Illumina sequencing. Mse-miR-929 and mse-miR-1b (antisense microRNA of mse-miR-1) were predicted in the previous study and now found to be conserved microRNAs in the tissue samples. We also found four novel microRNAs, two of which result from a gene cluster. Mse-miR-281-star, mse-miR-965-star, mse-miR-31-star, and mse-miR-9b-star were present at higher levels than their respective mature strands. Abundance changes of microRNAs were observed after the immune challenge. Based on the quantitative data of mRNA levels in control and induced fat body and hemocytes as well as the results of microRNA target site prediction, we suggest that certain microRNAs and microRNA*s regulate gene expression for pattern recognition, prophenoloxidase activation, cellular responses, antimicrobial peptide synthesis, and conserved intracellular signal transduction (Toll, IMD, JAK-STAT, MAPK-JNK-p38, and small interfering RNA pathways). In summary, this study has enriched our knowledge on M. sexta microRNAs and how some of them may participate in the expression regulation of immunity-related genes.

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.

MiR-182 is up-regulated and targeting Cebpa in hepatocellular carcinoma

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that repress their targets at post transcriptional level. Existing studies have shown that miRNAs are important regulatory genes in hepatocellular carcinoma (HCC), as either tumor suppressors or oncogenes. MiR-122 is normally downregulated in HCC and regarded as a tumor suppressor. Recently miR-122 has been reported to be regulated by CEBPA, which is then involved in a novel pathway to influence proliferation of tumor cells. However it is unknown whether CEBPA is regulated by miRNAs in HCC. In this study, we find that miR-182 is upregulated in HCC model rat, and represses CEBPA in both rat and human. This further improves the current CEBPA/miR-122 pathway that controls the proliferation of tumor cells. These results suggest that miR-182 is a potential oncogene in HCC and could be used as a diagnostic marker and drug target of HCC.

Development associated profiling of chitinase and microRNA of Helicoverpa armigera identified chitinase repressive microRNA

Expression of chitinase is developmentally regulated in insects in consonance with their molting process. During the larval-larval metamorphosis in Helicoverpa armigera, chitinase gene expression varies from high to negligible. In the five-day metamorphic course of fifth-instar larvae, chitinase transcript is least abundant on third day and maximal on fifth day. MicroRNA library prepared from these highest and lowest chitinase-expressing larval stages resulted in isolation of several miRNAs. In silico analysis of sequenced miRNAs revealed three miRNAs having sequence similarity to 3'UTR of chitinase. Gene-targeted specific action of these miRNAs, was investigated by luciferase reporter having 3'UTR of chitinase. Only one of three miRNAs, miR-24, inhibited luciferase expression. Further, a day-wise in vivo quantification of miR-24 in fifth-instar larvae revealed a negative correlation with corresponding chitinase transcript abundance. The force-feeding of synthetic miR-24 induced significant morphological aberrations accompanied with arrest of molting. These miR-24 force-fed larvae revealed significantly reduced chitinase transcript abundance.

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.

Characterization and comparative profiling of MicroRNAs in a sexual dimorphism insect, Eupolyphaga sinensis Walker

BACKGROUND

MicroRNAs are now recognized as key post-transcriptional regulators in animal ontogenesis and phenotypic diversity. Eupolyphaga sinensis Walker (Blattaria) is a sexually dimorphic insect, which is also an important source of material used in traditional Chinese medicine. The male E. sinensis have shorter lifecycles and go through fewer instars than the female. Furthermore, the males have forewings, while the females are totally wingless.

RESULTS

We used the Illumina/Solexa deep sequencing technology to sequence small RNA libraries prepared from the fourth-instar larvae of male and female E. sinensis. 19,097,799 raw reads were yielded in total: 7,817,445 reads from the female library and 11,280,354 from the male, respectively. As a result, we identified 168 known miRNAs belonging to 55 families as well as 204 novel miRNAs. Moreover, 45 miRNAs showed significantly different expression between the female and the male fourth-instar larvae, and we validated 10 of them by Stem-loop qRT-PCR. Some of these differentially expressed miRNAs are related to metamorphosis, development and phenotypic diversity.

CONCLUSIONS/SIGNIFICANCE

This is the first comprehensive description of miRNAs in E. sinensis. The results provide a useful resource for further in-depth study on molecular regulation and evolution of miRNAs. These findings not only enrich miRNAs for hemimetabolans but also lay the foundation for the study of post-transcriptional regulation on the phenomena of sexual dimorphism.