Insect-Specific microRNA Involved in the Development of the Silkworm Bombyx mori

Transgenic overexpression of bmo-miR-6498-5p increases resistance to Nosema bombycis in the silkworm, Bombyx mori

Microsporidia are unfriendly microorganisms, and their infections cause considerable damage to economically or environmentally important insects like silkworms and honeybees. Thus, the identification of measures to improve host resistance to microsporidia infections is critically needed. Here, an overexpressed miR-6498-5p transgenic silkworm line was constructed. Importantly, the survival rates and median lethal doses of the transgenic line were clearly higher after infection with Nosema bombycis. H&E staining and RT-qPCR analyses revealed an inhibitory effect on the proliferation of N. bombycis in the transgenic larvae. Metabolomics analysis further revealed the presence of 56 differential metabolites between the two lines. KEGG analysis of these 56 metabolites found that they were involved in various amino acid and vitamin metabolism pathways. Notably, VB6 metabolism was enriched among the metabolites, and the pathway was well known for its involvement in the synthesis, interconversion, and degradation of amino acids. These suggest that miR-6498-5p modifies parasitic environments to inhibit the proliferation of N. bombycis by affecting the host amino acid metabolism. These results demonstrate the potential of microRNAs as biomolecules that can promote resistance to microsporidia and provide new insights and a new approach to generate microsporidia-resistant biological materials.IMPORTANCEMicrosporidia have an extremely wide host range and are capable of infecting a wide variety of insects and vertebrates, including humans, and their lethality to multiple species often poses significant environmental management challenge. Here, we successfully constructed a microsporidium-resistant line in the silkworm, based on the overexpression of miR-6498-5p. Our results strongly support the hypothesis that miR-6498-5p efficiently suppresses the proliferation of Nosema bombycis by regulating the host VB6 metabolism, a key pathway for enzymes involved in amino acid transport and protein metabolism. Our study provides new insights for understanding host anti-pathogen defenses toward microsporidia.

MicroRNAs in the developmental toolbox - a comparative approach to understanding their role in regulating insect development

MicroRNAs are ubiquitous in the genomes of metazoan. Since their discovery during the late 20th century, our understanding of these small, noncoding RNAs has grown rapidly. However, there are still many unknowns about the functional significance of miRNAs - especially in non-model insects. Here I discuss the accumulating evidence that microRNAs are part of gene regulatory networks that determine not only the developmental outcome and/or mediate transitions between stages and alternative developmental pathways. During the last 20 years, researchers have published a multitude of profiling studies that describe changes in miRNAs that may be important for development and catalog potential targets. Proof-of-principle studies document phenotypic changes that occur when candidate genes and/or miRNAs are inhibited or overexpressed. Studies that use both of these approaches along with methods for confirming miRNA-mRNA interaction demonstrate necessary roles for miRNAs withing gene networks. Together all of these types of studies provide essential clues for understanding the function of miRNAs in the developmental toolbox.

Identification and Prediction of Differentially Expressed MicroRNAs Associated with Detoxification Pathways in Larvae of Spodoptera frugiperda

Spodoptera frugiperda poses a severe threat to crops, causing substantial economic losses. The increased use of chemical pesticides has led to resistance in S. frugiperda populations. Micro ribonucleic acids (MicroRNAs or miRNAs) are pivotal in insect growth and development. This study aims to identify miRNAs across different developmental stages of S. frugiperda to explore differential expression and predict target gene functions. High-throughput sequencing of miRNAs was conducted on eggs, 3rd instar larvae, pupae, and adults. Bioinformatics analyses identified differentially expressed miRNAs specifically in larvae, with candidate miRNAs screened to predict target genes, particularly those involved in detoxification pathways. A total of 184 known miRNAs and 209 novel miRNAs were identified across stages. Comparative analysis revealed 54, 15, and 18 miRNAs differentially expressed in larvae, compared to egg, pupa, and adult stages, respectively. Eight miRNAs showed significant differential expression across stages, validated by quantitative reverse transcription PCR (qRT-PCR). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses predicted target genes' functions, identifying eight differentially expressed miRNAs targeting 10 gene families associated with detoxification metabolism, including P450s, glutathione S-transferase (GSTs), ATP-binding cassette (ABC) transporters, and sodium channels. These findings elucidate the species-specific miRNA profiles and regulatory mechanisms of detoxification-related genes in S. frugiperda larvae, offering insights and strategies for effectively managing this pest.

miR-306-5p is involved in chitin metabolism in Aedes albopictus pupae via linc8338-miR-306-5p-XM_019678125.2 axis

Aedes albopictus can transmit several lethal arboviruses. This mosquito has become a sever public health threat due to its rapidly changing global distribution. Chitin, which is the major component of the cuticle and peritrophic membrane (PM), is crucial for the growth and development of insect. microRNAs (miRNAs) play important roles in the posttranscriptional level regulation of gene expression, thereby influencing many biological processes in insects. In this study, an attempt was made to evaluate the role of miR-306-5p in regulating chitin metabolism in Ae. albopictus pupae. Overexpression of miR-306-5p resulted in a significantly reduced survival rate in pupae and an increased malformation rate in adults. Both in vivo and in vitro evidence confirmed the presence of the competing endogenous RNA (ceRNA) regulatory axis (linc8338-miR-306-5p-XM_019678125.2). RNAi of linc8338 and XM_019678125.2 had effects on pupae similar to those of miR-306-5p. The highest expression level of miR-306-5p was found in the midgut, and alteration in the expression of miR-306-5p, XM_019678125.2 and linc8338 induced increased transcript levels of chitin synthase 2 (AaCHS2) and decreased chitinase 10 (AaCht10); as well as increased thickness of the midgut and enlarged midgut epithelial cells. The results of this study highlight the potential of miR-306-5p as a prospective target in mosquito control and confirm that the ceRNA mechanism is involved in chitin metabolism. These findings will provide a basis for further studies to uncover the molecular mechanisms through which ncRNAs regulate chitin metabolism.

Identification of MicroRNAs in the West Nile Virus Vector Culex tarsalis (Diptera: Culicidae)

MicroRNAs (miRNAs) are a group of small noncoding RNAs that regulate gene expression during important biological processes including development and pathogen defense in most living organisms. Presently, no miRNAs have been identified in the mosquito Culex tarsalis (Diptera: Culicidae), one of the most important vectors of West Nile virus (WNV) in North America. We used small RNA sequencing data and in vitro and in vivo experiments to identify and validate a repertoire of miRNAs in Cx. tarsalis mosquitoes. Using bioinformatic approaches we analyzed small RNA sequences from the Cx. tarsalis CT embryonic cell line to discover orthologs for 86 miRNAs. Consistent with other mosquitoes such as Aedes albopictus and Culex quinquefasciatus, miR-184 was found to be the most abundant miRNA in Cx. tarsalis. We also identified 20 novel miRNAs from the recently sequenced Cx. tarsalis genome, for a total of 106 miRNAs identified in this study. The presence of selected miRNAs was biologically validated in both the CT cell line and in adult Cx. tarsalis mosquitoes using RT-qPCR and sequencing. These results will open new avenues of research into the role of miRNAs in Cx. tarsalis biology, including development, metabolism, immunity, and pathogen infection.

Characterization of MicroRNAs Associated with Reproduction in the Brown Planthopper, Nilaparvata lugens

Insects have a robust capacity to produce offspring for propagation, and the reproductive events of female insects have been achieved at the molecular and physiological levels via regulatory gene pathways. However, the roles of MicroRNAs (miRNAs) in the reproductive development of the brown planthopper (BPH), Nilaparvata lugens, remain largely unexplored. To understand the roles of miRNAs in reproductive development, miRNAs were identified by Solexa sequencing in short-winged (SW) female adults of BPH. Small RNA libraries derived from three developmental phases (1 day, 3 days, and 5 days after emergence) were constructed and sequenced. We identified 905 miRNAs, including 263 known and 642 novel miRNAs. Among them, a total of 43 miRNAs were differentially expressed in the three developmental phases, and 14,568 putative targets for 43 differentially expressed miRNAs (DEMs) were predicted by TargetScan and miRanda. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the predicted miRNA targets illustrated the putative roles for these DEMs in reproduction. The progress events were annotated, including oogenesis, lipid biosynthetic process, and related pathways such as apoptosis, ABC transporters, and amino acid metabolism. Four highly abundant DEMs (miR-9a-5p, miR-34-5p, miR-275-3p, and miR-317-3p) were further screened, and miR-34-5p was confirmed to be involved in the regulation of reproduction. Overexpression of miR-34-5p via injecting its mimics reduced fecundity and decreased Vg expression. Moreover, target genes prediction for miR-34-5p showed they might be involved in 20E signaling cascades, apoptosis, and gonadal development, including hormone receptor 4 (HR4), caspase-1 (Cp-1), and spermatogenesis-associated protein 20 (SPATA20). These findings provide a valuable resource for future studies on the role of miRNAs in BPH reproductive development.

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.

Epigenetic Molecular Mechanisms in Insects

Insects are the largest animal group on Earth both in biomass and diversity. Their outstanding success has inspired genetics and developmental research, allowing the discovery of dynamic process explaining extreme phenotypic plasticity and canalization. Epigenetic molecular mechanisms (EMMs) are vital for several housekeeping functions in multicellular organisms, regulating developmental, ontogenetic trajectories and environmental adaptations. In Insecta, EMMs are involved in the development of extreme phenotypic divergences such as polyphenisms and eusocial castes. Here, we review the history of this research field and how the main EMMs found in insects help to understand their biological processes and diversity. EMMs in insects confer them rapid response capacity allowing insect either to change with plastic divergence or to keep constant when facing different stressors or stimuli. EMMs function both at intra as well as transgenerational scales, playing important roles in insect ecology and evolution. We discuss on how EMMs pervasive influences in Insecta require not only the control of gene expression but also the dynamic interplay of EMMs with further regulatory levels, including genetic, physiological, behavioral, and environmental among others, as was earlier proposed by the Probabilistic Epigenesis model and Developmental System Theory.

Analysis of microRNA expression profiles dynamic in different life stages of Haemaphysalis longicornis ticks by deep sequencing of small RNA libraries

The three-host tick Haemaphysalis longicornis is an obligate blood-sucking ectoparasite. In life-stage transitions, microRNAs (miRNAs) show a variety of expression changes. To investigate these changes, deep sequencing technology was applied to identify the conserved and potentially novel miRNAs expressed during the different life stages of H. longicornis. Total RNA from eggs, unfed larvae, unfed nymphs and unfed adults was extracted for deep sequence analysis. Deep sequencing on a Hiseq 4000 generated a total of 111,192,069 reads, grouped into four small RNA (sRNA) libraries, one for each of the four developmental stages of H. longicornis. Among these sequences, 78 conserved and 55 potentially novel miRNAs were identified, including stage-specific and differentially expressed miRNAs. Gene ontology (GO) analysis indicated significantly enriched GO terms related to cell proliferation and differentiation, including specific terms for the processes of development, growth, metabolism, regulation of biological functions, reproduction, and membrane enzyme regular activity. Kyoto Encyclopedia of Gene and Genomes (KEGG) analysis revealed a significant enrichment of the insulin, notch, Hippo, and Wnt signaling pathways for growth and development. Our data highlight the abundance of miRNA changes (conserved and potentially novel) in the different life stages of H. longicornis. In particular, stage-specific miRNAs, as observed, are essential regulators for the development of H. longicornis.

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.

Genome editing in Bombyx mori: New opportunities for silkworm functional genomics and the sericulture industry

In recent years, research in life sciences has been remarkably revolutionized owing to the establishment, development and application of genome editing technologies. Genome editing has not only accelerated fundamental research but has also shown promising applications in agricultural breeding and therapy. In particular, the clustered, regularly interspaced, short palindromic repeat (CRISPR) technology has become an indispensable tool in molecular biology owing to its high efficacy and simplicity. Genome editing tools have also been established in silkworm (Bombyx mori), a model organism of Lepidoptera insects with high economic importance. This has remarkably improved the level and scope of silkworm research and could reveal new mechanisms or targets in basic entomology and pest management studies. In this review, we summarize the progress and potential of genome editing in silkworm and its applications in functional genomic studies for generating novel genetic materials.

Selection of reference genes for expression analysis of plant-derived microRNAs in Plutella xylostella using qRT-PCR and ddPCR

The establishment of an expression quantification system that can be easily applied for the comparison of microRNAs (miRNAs) from biological samples is an important step toward understanding functional mechanisms in organisms. However, there is lack of attention on the selection of reference genes for miRNA expression profiling in insect herbivores. Here, we explored the candidate reference genes in a notorious pest of cruciferous crops, Plutella xylostella, for normalization of miRNA expression in developmental stages and tissues and in response to a change of food source from artificial diet to host plant Arabidopsis thaliana. We first compared the expression levels and stability of eight small RNAs using qRT-PCR, and found that miR11 was the most suitable reference gene for expression quantification of the miRNAs. We then confirmed this finding using digital droplet PCR and further validated with a well-studied cross-kingdom miRNA derived from A. thaliana (ath-miR159a). However, none of the reference genes was applicable for all experimental conditions, and multiple reference genes were sometimes required within the same experiment. Our work provides a method for the selection of reference genes for quantification of plant-derived miRNAs, which paves the way for unveiling their roles in the insect-plant coevolution.

Genome-Wide Screening and Functional Analysis Reveal That the Specific microRNA nlu-miR-173 Regulates Molting by Targeting Ftz-F1 in Nilaparvata lugens

Background: Molting is a crucial physiological behavior during arthropod growth. In the past few years, molting as well as chitin biosynthesis triggered by molting, is subject to regulation by miRNAs. However, how many miRNAs are involved in insect molting at the genome-wide level remains unknown. Results: We deeply sequenced four samples obtained from nymphs at the 2nd-3rd and 4th-5th instars, and then identified 61 miRNAs conserved in the Arthropoda and 326 putative novel miRNAs in the brown planthopper Nilaparvata lugens, a fearful pest of rice. A total of 36 mature miRNAs with significant different expression levels at the genome scale during molting, including 19 conserved and 17 putative novel miRNAs were identified. After comparing the expression profiles, we found that most of the targets of 36 miRNAs showing significantly differential expression were involved in energy and hormone pathways. One of the 17 putative novel miRNAs, nlu-miR-173 was chosen for functional study. nlu-miR-173 acts in 20-hydroxyecdysone signaling through its direct target, N. lugens Ftz-F1(NlFtz-F1), a transcription factor. Furthermore, we found that the transcription of nlu-miR-173 was promoted by Broad-Complex (BR-C), suggesting that its involvement in the 20-hydroxyecdysone pathway contributes to proper molting function. Conclusion: We provided a comprehensive resource of miRNAs associated with insect molting and identified a novel miRNA as a potential target for pest control.

MicroRNA Expression Analysis of Naked Silkworms

The silk gland (SG) is characterized by the synthesis and secretion of silk protein in the economically important silkworm, Bombyx mori L. (Lepidoptera: Bombycidae). Nd and Nd-s are two fibroin-secretion-deficient silkworm mutants. MicroRNA (miRNA) plays an important role in many biological processes, such as cell proliferation, differentiation, and apoptosis. Using the Dazao silkworm as a control, we explored the miRNA expression profiles in the SGs of u02 (Nd) and u05 (Nd-s) to reveal the potential functions of miRNAs in silk protein expression and SG development. Here, we sequenced small RNA libraries made from the whole SGs of three strains. There are 260, 236, and 233 known miRNAs and 20, 18, and 18 potential new miRNAs identified from Dazao, u02, and u05, respectively. Fifty-three miRNAs are differentially expressed between Dazao and u02, 51 between Dazao and u05, and 16 between u02 and u05. Gene ontology/KEGG analyses show that most of the predicted target genes of differentially expressed miRNAs were assigned to functional categories involved in cell proliferation, organ development, and cellular compartment structures. The miRNA expression profile of naked silkworms will pave the way for the understanding of SG development and the regulation of silk protein expression.

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.

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.

Daphnia magna microRNAs respond to nutritional stress and ageing but are not transgenerational

Maternal effects, where the performance of offspring is determined by the condition of their mother, are widespread and may in some cases be adaptive. The crustacean Daphnia magna shows strong maternal effects: offspring size at birth and other proxies for fitness are altered when their mothers are older or when mothers have experienced dietary restriction. The mechanisms for this transgenerational transmission of maternal experience are unknown, but could include changes in epigenetic patterning. MicroRNAs (miRNAs) are regulators of gene expression that have been shown to play roles in intergenerational information transfer, and here, we test whether miRNAs are involved in D. magna maternal effects. We found that miRNAs were differentially expressed in mothers of different ages or nutritional state. We then examined miRNA expression in their eggs, their adult daughters and great granddaughters, which did not experience any treatments. The maternal (treatment) generation exhibited differential expression of miRNAs, as did their eggs, but this was reduced in adult daughters and lost by great granddaughters. Thus, miRNAs are a component of maternal provisioning, but do not appear to be the cause of transgenerational responses under these experimental conditions. MicroRNAs may act in tandem with egg provisioning (e.g., with carbohydrates or fats), and possibly other small RNAs or epigenetic modifications.

Multiple miRNAs jointly regulate the biosynthesis of ecdysteroid in the holometabolous insects, Chilo suppressalis

The accurate rise and fall of active hormones is important for insect development. The ecdysteroids must be cleared in a timely manner. However, the mechanism of suppressing the ecdysteroid biosynthesis at the right time remains unclear. Here, we sequenced a small RNA library of Chilo suppressalis and identified 300 miRNAs in this notorious rice insect pest. Microarray analysis yielded 54 differentially expressed miRNAs during metamorphosis development. Target prediction and in vitro dual-luciferase assays confirmed that seven miRNAs (two conserved and five novel miRNAs) jointly targeted three Halloween genes in the ecdysteroid biosynthesis pathway. Overexpression of these seven miRNAs reduced the titer of 20-hydroxyecdysone (20E), induced mortality, and retarded development, which could be rescued by treatment with 20E. Comparative analysis indicated that the miRNA regulation of metamorphosis development is a conserved process but that the miRNAs involved are highly divergent. In all, we present evidence that both conserved and lineage-specific miRNAs have crucial roles in regulating development in insects by controlling ecdysteroid biosynthesis, which is important for ensuring developmental convergence and evolutionary diversity.

Differentially expressed microRNAs in diapausing versus HCl-treated Bombyx embryos

Differentially expressed microRNAs were detected to explore the molecular mechanisms of diapause termination. The total small RNA of diapause-destined silkworm eggs and HCl-treated eggs was extracted and then sequenced using HiSeq high-throughput method. 44 novel miRNAs were discovered. Compared to those in the diapause-destined eggs, 61 miRNAs showed significant changes in the acid-treated eggs, with 23 being up-regulated and 38 being down-regulated. The potential target genes of differentially expressed miRNAs were predicted by miRanda. Gene Ontology and KEGG pathway enrichment analysis of these potential target genes revealed that they were mainly located within cells and organelles, involved in cellular and metabolic processes, and participated in protein production, processing and transportation. Two differentially expressed genes, Bombyx mori SDH and Bmo-miR-2761-3p, were further analyzed with qRT-PCR. BmSDH was significantly up-regulated in the HCl-treated eggs, while Bmo-miR-2761-3p was down-regulated. These results suggested that these two genes were well coordinated in silkworm eggs. Dual luciferase reporter assay demonstrated that Bmo-miR-2761-3p inhibited the expression of BmSDH.

Small RNA sequencing reveals a role for sugarcane miRNAs and their targets in response to Sporisorium scitamineum infection

BACKGROUND

Sugarcane smut caused by Sporisorium scitamineum leads to a significant reduction in cane yield and sucrose content. MicroRNAs (miRNAs) play an important role in regulating plant responses to biotic stress. The present study was the first to use two sugarcane genotypes, YA05-179 (smut-resistant) and ROC22 (smut-susceptible), to identify differentially expressed miRNAs in sugarcane challenged with S. scitamineum by using high-throughput sequencing.

RESULTS

The predicted target gene number corresponding to known differentially expressed miRNAs in YA05-179 was less than that in ROC22, however most of them were in common. Expression of differential miRNAs under S. scitamineum challenge was mostly downregulated, with similar trends in the two varieties. Gene ontology (GO) analysis showed that the target gene classification of known miRNAs was similar to that of the newly identified miRNAs. These were mainly associated with cellular processes and metabolic processes in the biological process category, as well as combination and catalytic activity in the molecular function category. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed that these predicted target genes involved in a series of physiological and biochemical pathways or disease resistance-related physiological metabolism and signal transduction pathways, suggesting that the molecular interaction mechanism between sugarcane and S. scitamineum was a complex network system. These findings also showed certain predicted target genes of miR5671, miR5054, miR5783, miR5221, and miR6478 play roles in the mitogen-activated protein kinase (MAPK) signaling pathway, plant hormone signal transduction, and plant-pathogen interaction. Quantitative real-time PCR (qRT-PCR) analysis showed that majority of the known miRNAs and its predicted target genes followed a negatively regulated mode. Seven out of eight predicted target genes showed identical expression after 12 h treatment and reached the highest degree of matching at 48 h, indicating that the regulatory role of miRNAs on the target genes in sugarcane was maximized at 48 h after S. scitamineum challenge.

CONCLUSIONS

Taken together, our findings serve as evidence for the association of miRNA expression with the molecular mechanism underlying the pathogenesis of sugarcane smut, particularly on the significance of miRNA levels in relation to the cultivation of smut-resistant sugarcane varieties.

Identification and characterization of microRNAs expressed in the African malaria vector Anopheles funestus life stages using high throughput sequencing

BACKGROUND

Over the past several years, thousands of microRNAs (miRNAs) have been identified in the genomes of various insects through cloning and sequencing or even by computational prediction. However, the number of miRNAs identified in anopheline species is low and little is known about their role. The mosquito Anopheles funestus is one of the dominant malaria vectors in Africa, which infects and kills millions of people every year. Therefore, small RNA molecules isolated from the four life stages (eggs, larvae, pupae and unfed adult females) of An. funestus were sequenced using next generation sequencing technology.

RESULTS

High throughput sequencing of four replicates in combination with computational analysis identified 107 mature miRNA sequences expressed in the An. funestus mosquito. These include 20 novel miRNAs without sequence identity in any organism and eight miRNAs not previously reported in the Anopheles genus but are known in non-anopheles mosquitoes. Finally, the changes in the expression of miRNAs during the mosquito development were determined and the analysis showed that many miRNAs have stage-specific expression, and are co-transcribed and co-regulated during development.

CONCLUSIONS

This study presents the first direct experimental evidence of miRNAs in An. funestus and the first profiling study of miRNA associated with the maturation in this mosquito. Overall, the results indicate that miRNAs play important roles during the growth and development. Silencing such molecules in a specific life stage could decrease the vector population and therefore interrupt malaria transmission.

Expressed microRNA associated with high rate of egg production in chicken ovarian follicles

MicroRNA (miRNA) is a highly conserved class of small noncoding RNA about 19-24 nucleotides in length that function in a specific manner to post-transcriptionally regulate gene expression in organisms. Tissue miRNA expression studies have discovered a myriad of functions for miRNAs in various aspects, but a role for miRNAs in chicken ovarian tissue at 300 days of age has not hitherto been reported. In this study, we performed the first miRNA analysis of ovarian tissues in chickens with low and high rates of egg production using high-throughput sequencing. By comparing low rate of egg production chickens with high rate of egg production chickens, 17 significantly differentially expressed miRNAs were found (P < 0.05), including 11 known and six novel miRNAs. We found that all 11 known miRNAs were involved mainly in pathways of reproduction regulation, such as steroid hormone biosynthesis and dopaminergic synapse. Additionally, expression profiling of six randomly selected differentially regulated miRNAs were validated by quantitative real-time polymerase chain reaction (RT-qPCR). Some miRNAs, such as gga-miR-34b, gga-miR-34c and gga-miR-216b, were reported to regulate processes such as proliferation, cell cycle, apoptosis and metastasis and were expressed differentially in ovaries of chickens with high rates of egg production, suggesting that these miRNAs have an important role in ovary development and reproductive management of chicken. Furthermore, we uncovered that a significantly up-regulated miRNA-gga-miR-200a-3p-is ubiquitous in reproduction-regulation-related pathways. This miRNA may play a special central role in the reproductive management of chicken, and needs to be further studied for confirmation.

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 microRNAs in the white-backed planthopper, Sogatella furcifera

MicroRNAs (miRNAs) are a novel class of small, non-coding endogenous RNAs that play critical regulatory roles in many metabolic activities in eukaryotes. Reports of the identification of miRNAs in Sogatella furcifera (white-backed planthopper), the insect that acts as the only confirmed vector of the southern rice black-streaked dwarf virus (SRBSDV), are limited. In this study, a total of 382 miRNAs were identified in S. furcifera, including 106 conserved and 276 novel miRNAs, using high-throughput sequencing based on two small RNA libraries from viruliferous and non-viruliferous S. furcifera, and these miRNAs belonged to 52 conserved miRNA families and 58 S. furcifera-specific families, respectively. Comparison with miRNAs from 26 insect species and five other species in miRBase showed that more than half of the conserved miRNA families are highly conserved in Hexapoda, while other miRNAs are only conserved in non-dipterans. Furthermore, 4 117 target genes predicted for the 382 identified miRNAs could be categorized into 45 functional groups annotated by Gene Ontology. Compared with non-viruliferous cells, eight up-regulated miRNAs and four down-regulated miRNAs were identified in cells inoculated with SRBSDV, among which miR-14 and miR-n98a may be involved in the immune response to SRBSDV infection. Analyses of the identified miRNAs will provide insights into the roles of these miRNAs in the regulation and expression of genes involved in the metabolism, development and viral infection of S. furcifera.

Identification of Differentially Expressed microRNAs between the Fenpropathrin Resistant and Susceptible Strains in Tetranychus cinnabarinus

The carmine spider mite (Tetranychus cinnabarinus) is one of the most serious pests on crops and its control mainly depends on chemical acaricides. The excessive and improper acaricides use has resulted in mite resistance to many acaricides, including fenpropathrin. Previous studies have indicated fenpropathrin resistance is a complex development process involving many genes, but information on resistance mechanism of post-transcription regulation is rare. Using Illumina sequencing, several categories of sRNAs were identified from susceptible (TS) and fenpropathrin-resistant strains (TR) of T. cinnabarinus, including 75 known microRNAs (miRNAs) and 64 novel miRNAs, whose target genes containing 78592 miRNA-target pairs were predicted by 6 algorithms. Also, 12 significantly differently expressed miRNAs were identified between the TS and TR libraries and RT-qPCR validation also performed a well consistency with sequencing. The targets of significantly differentially expressed miRNAs included 7 glutathione S-transferase, 7 cytochrome P450 and 16 carboxyl/choline esterase genes, their function in fenpropathrin resistance were further analyzed. The present study provides the firstly large-scale characterization of miRNAs in T. cinnabarinus and the comparison between TS and TR strains gives a clue on how miRNA involves in fenpropathrin resistance.

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.

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.

Developmental Changes for the Hemolymph Metabolome of Silkworm (Bombyx mori L.)

Silkworm (Bombyx mori) is a lepidopteran-holometabolic model organism. To understand its developmental biochemistry, we characterized the larval hemolymph metabonome from the third instar to prepupa stage using (1)H NMR spectroscopy whilst hemolymph fatty acid composition using GC-FID/MS. We unambiguously assigned more than 60 metabolites, among which tyrosine-o-β-glucuronide, mesaconate, homocarnosine, and picolinate were reported for the first time from the silkworm hemolymph. Phosphorylcholine was the most abundant metabolite in all developmental stages with exception for the periods before the third and fourth molting. We also found obvious developmental dependence for the hemolymph metabonome involving multiple pathways including protein biosyntheses, glycolysis, TCA cycle, the metabolisms of choline amino acids, fatty acids, purines, and pyrimidines. Most hemolymph amino acids had two elevations during the feeding period of the fourth instar and prepupa stage. Trehalose was the major blood sugar before day 8 of the fifth instar, whereas glucose became the major blood sugar after spinning. C16:0, C18:0 and its unsaturated forms were dominant fatty acids in hemolymph. The developmental changes of hemolymph metabonome were associated with dietary nutrient intakes, biosyntheses of cell membrane, pigments, proteins, and energy metabolism. These findings offered essential biochemistry information in terms of the dynamic metabolic changes during silkworm development.

Identification and characteristics of microRNAs from army worm, Spodoptera frugiperda cell line Sf21

microRNAs play important regulatory role in all intrinsic cellular functions. Amongst lepidopteran insects, miRNAs from only Bombyx mori have been studied extensively with a little focus on Spodoptera sp. In the present study, we identified a total of 226 miRNAs from Spodoptera frugiperda cell line Sf21. Of the total, 116 miRNAs were well conserved within other insects, like B. mori, Drosophila melanogaster and Tribolium castenum while the remaining 110 miRNAs were identified as novel based on comparative analysis with the insect miRNA data set. Landscape distribution analysis based on Sf21 genome assembly revealed clustering of few novel miRNAs. A total of 5 miRNA clusters were identified and the largest one encodes 5 miRNA genes. In addition, 12 miRNAs were validated using northern blot analysis and putative functional role assignment for 6 Sf miRNAs was investigated by examining their relative abundance at different developmental stages of Spodoptera litura and body parts of 6^th instar larvae. Further, we identified a total of 809 potential target genes with GO terms for selected miRNAs, involved in different metabolic and signalling pathways of the insect. The newly identified miRNAs greatly enrich the repertoire of insect miRNAs and analysis of expression profiles reveal their involvement at various steps of biochemical pathways of the army worm.

Chorion genes: a landscape of their evolution, structure, and regulation

Differential regulation at the level of transcription provides a means for controlling gene expression in eukaryotes, especially during development. Insect model systems have been extensively used to decipher the molecular basis of such regulatory cascades, and one of the oldest such model systems is the regulation of chorion gene expression during ovarian follicle maturation. Recent experimental and technological advances have shed new light onto the system, allowing us to revisit it. Thus, in this review we try to summarize almost 40 years' worth of studies on chorion gene regulation while-by comparing Bombyx mori and Drosophila melanogaster models-attempting to present a comprehensive, unified model of the various regulatory aspects of choriogenesis that takes into account the evolutionary conservation and divergence of the underlying mechanisms.

Identification and molecular structure analysis of a new noncoding RNA, a sbRNA homolog, in the silkworm Bombyx mori genome

The small noncoding group of RNAs called stem-bulge RNAs (sbRNAs), first reported in Caenorhabditis elegans, is described as molecules homologous to the Y RNAs, a specific class of noncoding RNAs that is present in vertebrates. This homology indicates the possibility of the existence of sbRNAs in other invertebrate organisms. In this work, we used bioinformatic tools and conserved sequences of sbRNAs from C. Elegans and Y RNAs to search for homologous sbRNA sequences in the Bombyx mori genome. This analysis led to the discovery of one noncoding gene, which was translated into RNA segments and comparatively analysed with segments from human and hamster Y RNAs and C. elegans sbRNAs in molecular dynamic simulations. This gene represents the first evidence for a new sbRNA-like noncoding RNA, the BmsbRNA gene, in this Lepidoptera genome.

MicroRNA Let-7 regulates molting and metamorphosis in the silkworm, Bombyx mori

MicroRNAs (miRNAs) are a class of endogenous, non-coding, regulatory RNA molecules that post-transcriptionally regulate gene expression by binding to the 3'UTRs of mRNA targets and thus cause their degradation or translational inhibition. In insects, important roles of miRNAs in various biological processes have been demonstrated in Drosophila melanogaster. However, biological roles of miRNAs are barely unveiled in the majority of insect species due to limited genetic tools. In the present study, we introduce the transgenic miRNA sponge (miR-SP) technology combining with the binary GAL4/UAS system in the domesticated silkworm, Bombyx mori, to exploit the biological function of an evolutionally conserved miRNA, let-7. We successfully established transgenic silkworm lines in which a miRNA sponge construct targeting BmLet-7 seed region was expressed in a ubiquitous manner directed by A3-GAL4 driver. Transgenic animals showed decreased expression of BmLet-7, leading to developmental arrestment during the larval-larval and larval-pupal transition. Simultaneously, expression levels of the predicted BmLet-7 target genes, FTZ-F1 and Eip74EF (E74), key regulatory factors in the ecdysone pathway, were elevated in transgenic animals. The current study is the first report on application of the transgenic miR-SP technology in non-drosophilid insects, which will not only contribute to better understanding of let-7 biological roles, but also greatly facilitate future miRNA functional analysis in insects.

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.

Identification of conserved and novel microRNAs in the Pacific oyster Crassostrea gigas by deep sequencing

MicroRNAs (miRNAs) play important roles in regulatory processes in various organisms. To date many studies have been performed in the investigation of miRNAs of numerous bilaterians, but limited numbers of miRNAs have been identified in the few species belonging to the clade Lophotrochozoa. In the current study, deep sequencing was conducted to identify the miRNAs of Crassostrea gigas (Lophotrochozoa) at a genomic scale, using 21 libraries that included different developmental stages and adult organs. A total of 100 hairpin precursor loci were predicted to encode miRNAs. Of these, 19 precursors (pre-miRNA) were novel in the oyster. As many as 53 (53%) miRNAs were distributed in clusters and 49 (49%) precursors were intragenic, which suggests two important biogenetic sources of miRNAs. Different developmental stages were characterized with specific miRNA expression patterns that highlighted regulatory variation along a temporal axis. Conserved miRNAs were expressed universally throughout different stages and organs, whereas novel miRNAs tended to be more specific and may be related to the determination of the novel body plan. Furthermore, we developed an index named the miRNA profile age index (miRPAI) to integrate the evolutionary age and expression levels of miRNAs during a particular developmental stage. We found that the swimming stages were characterized by the youngest miRPAIs. Indeed, the large-scale expression of novel miRNAs indicated the importance of these stages during development, particularly from organogenetic and evolutionary perspectives. Some potentially important miRNAs were identified for further study through significant changes between expression patterns in different developmental events, such as metamorphosis. This study broadened the knowledge of miRNAs in animals and indicated the presence of sophisticated miRNA regulatory networks related to the biological processes in lophotrochozoans.

The expression analysis of silk gland-enriched intermediate-size non-coding RNAs in silkworm Bombyx mori

Small non-protein coding RNAs (ncRNAs) play important roles in development, stress response and other cellular processes. Silkworm is an important model for studies on insect genetics and control of Lepidopterous pests. We have previously identified 189 novel intermediate-size ncRNAs in silkworm Bombyx mori, including 40 ncRNAs that showed altered expression in different developmental stages. Here we characterized the functions of these 40 ncRNAs by measuring their expressions in six tissues of the fifth instar larvae using Northern blot and real-time polymerase chain reaction assays. We identified nine ncRNAs (four small nucleolar RNAs and five unclassified ncRNAs) that were enriched in silk gland, including four ncRNAs that showed silk gland-specific expression. We further showed that three of nine silk gland-enriched ncRNAs were predominantly expressed in the anterior silk gland, whereas another three ncRNAs were highly accumulated in the posterior silk gland, suggesting that they may play different roles in fibroin synthesis. Furthermore, an unclassified ncRNA, Bm-152, exhibited converse expression pattern with its antisense host gene gartenzwerg in diverse tissues, and might regulate the expression of gartenzwerg through RNA-protein complex. In addition, two silk gland-enriched ncRNAs Bm-102 and Bm-159 can be found in histone modification complex, which indicated that they might play roles through epigenetic modifications. Taken together, we provided the first expression and preliminary functional analysis of silk gland-enriched ncRNAs, which will help understand the molecular mechanism of silk gland-development and fibroin synthesis.

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.

Differential expression profiling of spleen microRNAs in response to two distinct type II interferons in Tetraodon nigroviridis

MicroRNAs are endogenous, small non-coding RNAs approximately 18-26 nucleotides in length that regulate target gene expression at the post-transcription level. Interferon-γ (IFN-γ) is a Th1 cytokine that is involved in both the innate and adaptive immune responses. We previously identified two IFN-γ genes in green-spotted puffer fish (Tetraodon nigroviridis). To determine whether miRNAs participate in IFN-γ-related immune responses, T. nigroviridis spleen cells were treated with recombinant IFN-γ isoforms, and a Solexa high-throughput sequencing method was used to identify miRNAs. In total, 1,556, 1,538 and 1,573 miRNAs were found in the three samples, and differentially expressed miRNAs were determined. In total, 398 miRNAs were differentially expressed after rIFN-γ1 treatment, and 438 miRNAs were differentially expressed after rIFN-γ2 treatment; additionally, 403 miRNAs were differentially expressed between the treatment groups. Ten differentially expressed miRNAs were chosen for validation using qRT-PCR. Target genes for the differentially expressed miRNAs were predicted, and GO and KEGG analyses were performed. This study provides basic knowledge regarding fish IFN-γ-induced miRNAs and offers clues for further studies into the mechanisms underlying fish IFN-γ-mediated immune responses.

Identification and characterization of microRNAs in the ovaries of multiple and uniparous goats (Capra hircus) during follicular phase

Background

Superior kidding rate is an important economic trait in production of meat goat, and ovulation rate is the precondition of kidding rate. MicroRNAs (miRNAs) play critical roles in almost all ovarian biological processes, including folliculogenesis, follicle development, follicle atresia, luteal development and regression. To find out the different ovarian activity and follicle recruitment with miRNA-mediated posttranscriptional regulation, the small RNAs expressed pattern in the ovarian tissues of multiple and uniparous Anhui White goats during follicular phase was analyzed using Solexa sequencing data.

Results

1008 miRNAs co-expressed, 309 and 433 miRNAs specifically expressed in the ovaries of multiple and uniparous goats during follicular phase were identified. The 10 most highly expressed miRNAs in the multiple library were also the highest expressed in the uniparous library, and there were no significantly different between each other. The highest specific expressed miRNA in the multiple library was miR-29c, and the one in the uniparous library was miR-6406. 35 novel miRNAs were predicted in total. GO annotation and KEGG Pathway analyses were implemented on target genes of all miRNA in two libraries. RT-PCR was applied to detect the expression level of 5 randomly selected miRNAs in multiple and uniparous hircine ovaries, and the results were consistent with the Solexa sequencing data.

Conclusions

In the present study, the different expression of miRNAs in the ovaries of multiple and uniparous goats during follicular phase were characterized and investigated using deep sequencing technology. The result will help to further understand the role of miRNAs in kidding rate regulation and also may help to identify miRNAs which could be potentially used to increase hircine ovulation rate and kidding rate in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-339) contains supplementary material, which is available to authorized users.

Identification and characterization of Cynoglossus semilaevis microRNA response to Vibrio anguillarum infection through high-throughput sequencing

MicroRNAs (miRNA) play key regulatory roles in diverse biological processes. Cynoglossus semilaevis is an important commercial mariculture fish species in China. To identify miRNAs and investigate immune-related miRNAs of C. semilaevis, we performed high-throughput sequencing on three small RNA libraries prepared from C. semilaevis immune tissues (liver, head kidney, spleen, and intestine). One library was prepared under normal conditions (control, CG); two were prepared during Vibrio anguillarum infection, where vibriosis symptoms were obvious and non-obvious (HOSG and NOSG, respectively). We obtained 11,216,875, 12,313,404, and 11,398,695 clean reads per library, respectively. Bioinformatic analysis identified 452 miRNAs, including 24 putative novel miRNAs. We analyzed differentially expressed miRNAs between two libraries using pairwise comparison. For NOSG-CG, there was significant differential expression of 175 (38.72%) miRNAs. There was significant differential expression of 215 (47.57%) miRNAs between HOSG and CG. Compared with CG, The HOSG-NOSG comparison revealed significantly different expression of 122 (26.99%) miRNAs respectively. Real-time quantitative PCR (RT-qPCR) experiments were performed for 10 miRNAs of the three samples, and agreement was found between the sequencing and RT-qPCR data. For miRNAs that were significantly differentially expressed, functional annotation of target genes by Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that a set of miRNAs that were expressed highly abundantly and significantly differentially were might involved in immune system development and immune response. To our understanding, this is the first report of comprehensive identification of C. semilaevis miRNAs being differentially regulated in immune tissues (liver, head kidney, spleen, and intestine) in normal conditions relating to V. anguillarum infection. Many miRNAs were differentially regulated upon pathogen exposure. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in C. semilaevis host-pathogen interactions.

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.

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.

Identification of differentially expressed genes in American cockroach ovaries and testes by suppression subtractive hybridization and the prediction of its miRNAs

Studies on the cockroach have contributed to our understanding of several important developmental processes, especially those that can be easily studied in the embryo. However, our knowledge on late events such as gonad differentiation in the cockroach is still limited. The major aim of the present study was to identify sex-specific genes between adult female and male Periplaneta americana. Two cDNA libraries were constructed using the suppression subtractive hybridization method; a total of 433 and 599 unique sequences were obtained from the forward library and the reverse library, respectively, by cluster assembly, and sequence alignment of 1,032 expressed sequence tags. The analysis of the differentially expressed gene functions allowed these genes to be categorized into three groups: biological process, molecular function, and cellular component. The differentially expressed genes were suggested to be related to the development of the gonads of P. americana. Twelve differentially expressed genes were randomly selected and verified using relative quantitative real-time polymerase chain reaction (qRT-PCR). Meanwhile, by adopting a range of filtering criteria, we predicted two potential microRNA sequences for P. americana, pam-miR100-3p and pam-miR7. To confirm the expression of potential microRNAs (miRNAs) in American cockroach, a qRT-PCR approach was also employed. The data presented here offer the insights into the molecular foundation of sex differences in American cockroach, and the first report for the miRNAs in this species. In addition, the results can be used as a reference for unraveling candidate genes associated with the sex and reproduction of cockroaches.

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.