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

Transcriptome-Wide Evaluation Characterization of microRNAs and Assessment of Their Functional Roles as Regulators of Diapause in Ostrinia furnacalis Larvae (Lepidoptera: Crambidae)

microRNAs (miRNAs) function as vital regulators of diapause in insects through their ability to post-transcriptionally suppress target gene expression. In this study, the miRNA of Ostrinia furnacalis, an economically important global crop pest species, was characterized. For the included analyses, 9 small RNA libraries were constructed using O. furnacalis larvae in different diapause states (non-diapause, ND; diapause, D; diapause-termination, DT). The results identified 583 total miRNAs, of which 256 had previously been identified, whereas 327 were novel. Furthermore, comparison analysis revealed that 119 and 27 miRNAs were differentially expressed in the D vs. ND and DT vs. D, respectively. Moreover, the expression patterns of their miRNAs were also analyzed. GO and KEGG analysis of the target genes of differentially expressed miRNAs highlighted the importance of these miRNAs as diapause regulators in O. furnacalis, especially through metabolic processes, endocrine processes, 20-hydroxyecdysone, and circadian clock signaling pathways. In summary, this study highlighted the involvement of specific miRNAs in the control of diapause in O. furnacalis. To the best of our knowledge, this is the first study to identify miRNA expression patterns in O. furnacalis, thereby providing reference and novel evidence enhancing our current understanding of how small RNAs influence insect diapause.

Conserved microRNAs miR-8-3p and miR-2a-3 targeting chitin biosynthesis to regulate the molting process of Sogatella furcifera (Horváth)(Hemiptera: Delphacidae)

Molting is a key solution to growth restriction in insects. The periodic synthesis and degradation of chitin, one of the major components of the insect epidermis, is necessary for insect growth. MicroRNA (miRNA) have been implicated in molting regulation, yet their involvement in the interplay interaction between the chitin synthesis pathway and 20-hydroxyecdysone signaling remains poorly understood. In this study, soluble trehalase (Tre1) and phosphoacetylglucosamine mutase (PAGM) were identified as targets of conserved miR-8-3p and miR-2a-3, respectively. The expression profiles of miR-8-3p-SfTre1 and miR-2a-3-SfPAGM exhibited an opposite pattern during the different developmental stages, indicating a negative regulatory relationship between them. This relationship was confirmed by an in vitro dual-luciferase reporter system. Overexpression of miR-8-3p and miR-2a-3 by injection of mimics inhibited the expression of their respective target genes and increased mortality, leading to death in the pre-molting, and molting death phenomena. They also caused a decrease in chitin content and expression levels of key genes in the chitin synthesis pathway (SfTre1, SfTre2, SfHK, SfG6PI, SfGFAT, SfGNA, SfPAGM, SfUAP, SfCHS1, SfCHS1a, and SfCHS1b). Conversely, the injection of miRNA inhibitors resulted in the upregulation of the expression levels of these genes. Following 20E treatment, the expression levels of miR-8-3p and miR-2a-3 decreased significantly, while their corresponding target genes increased significantly. These results indicate that miR-8-3p and miR-2a-3 play a regulatory role in the molting of Sogatella furcifera by targeting SfTre1 and SfPAGM, respectively. These findings provide new potential targets for the development of subsequent new control strategies.

Impact of Ae-GRD on Ivermectin Resistance and Its Regulation by miR-71-5p in Aedes aegypti

iGABAR, a member of the Cys-loop ligand-gated ion channel superfamily, is a significant target of the insecticide ivermectin (IVM). GRD is the potential subunit of the insect iGABAR. However, little information about GRD in Ae. aegypti has been reported. In this study, we involved cloning and characterizing the iGABAR subunit GRD of Ae. aegypti (Ae-GRD). Sequence analysis indicated that Ae-GRD, as part of the cysteine-loop ligand-gated ion channel family, is similar to other insect GRD. RNA interference (RNAi) was employed to explore IVM resistance in Ae. aegypti, resulting in a significant reduction in Ae-GRD expression (p < 0.05), and the mortality of Ae. aegypti adults with Ae-GRD knockdown was significantly decreased after exposure to ivermectin. Bioinformatics prediction identified miR-71-5p as a potential regulator of Ae-GRD. In vitro, dual-luciferase reporter assays confirmed that Ae-GRD expression was regulated by miR-71-5p. Microinjection of miR-71-5p mimics upregulated miR-71-5p expression and downregulated Ae-GRD gene expression, reducing mortality by 34.52% following IVM treatment. Conversely, microinjection of a miR-71-5p inhibitor decreased miR-71-5p expression but did not affect the susceptibility to IVM despite increased Ae-GRD expression (p < 0.05). In conclusion, Ae-GRD, as one of the iGABA receptor subunits, is a potential target of ivermectin. It may influence ivermectin resistance by modulating the GABA signaling pathway. The inhibition of Ae-GRD expression by miR-71-5p decreased ivermectin resistance and consequently lowered the mortality rate of Ae. aegypti mosquitoes. This finding provides empirical evidence of the relationship between Ae-GRD and its miRNA in modulating insecticide resistance, offering novel perspectives for mosquito control strategies.

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

Elucidating the ecophysiology of soybean pod-sucking stinkbug Riptortus pedestris (Hemiptera: Alydidae) based on de novo genome assembly and transcriptome analysis

Food security is important for the ever-growing global population. Soybean, Glycine max (L.) Merr., is cultivated worldwide providing a key source of food, protein and oil. Hence, it is imperative to maintain or to increase its yield under different conditions including challenges caused by abiotic and biotic stresses. In recent years, the soybean pod-sucking stinkbug Riptortus pedestris has emerged as an important agricultural insect pest in East, South and Southeast Asia. Here, we present a genomics resource for R. pedestris including its genome assembly, messenger RNA (mRNA) and microRNA (miRNA) transcriptomes at different developmental stages and from different organs. As insect hormone biosynthesis genes (genes involved in metamorphosis) and their regulators such as miRNAs are potential targets for pest control, we analyzed the sesquiterpenoid (juvenile) and ecdysteroid (molting) hormone biosynthesis pathway genes including their miRNAs and relevant neuropeptides. Temporal gene expression changes of these insect hormone biosynthesis pathways were observed at different developmental stages. Similarly, a diet-specific response in gene expression was also observed in both head and salivary glands. Furthermore, we observed that microRNAs (bantam, miR-14, miR-316, and miR-263) of R. pedestris fed with different types of soybeans were differentially expressed in the salivary glands indicating a diet-specific response. Interestingly, the opposite arms of miR-281 (-5p and -3p), a miRNA involved in regulating development, were predicted to target Hmgs genes of R. pedestris and soybean, respectively. These observations among others highlight stinkbug's responses as a function of its interaction with soybean. In brief, the results of this study not only present salient findings that could be of potential use in pest management and mitigation but also provide an invaluable resource for R. pedestris as an insect model to facilitate studies on plant-pest interactions.

MicroR-9c-5p and novel-mir50 co-target Akt to regulate Lasioderma serricorne reproduction

High fecundity is a common characteristic of insect pests which increases the difficulty of population control. Serine/threonine kinase Akt is an indispensable component of the insulin signaling pathway. Silencing of LsAkt severely hinders reproduction in Lasioderma serricorne, a stored product insect pest. However, the post-transcriptional pathway of LsAkt in L. serricorne remains unknown. This study identified 2 binding sites of miR-9c-5p and novel-mir50 in the coding sequences of LsAkt. The expression profiles of 2 microRNAs (miRNAs) and LsAkt displayed an opposite pattern during the adult stages. Luciferase reporter assay showed that novel-mir50 and miR-9c-5p could downregulate the expression of LsAkt. Overexpression of miR-9c-5p and novel-mir50 by injection of mimics inhibited the expression of LsAkt and reduced oviposition, decreased egg hatchability, and blocked ovarian development. It also decreased the expression of genes involved in ovarian development (LsVg and LsVgR) and the nutritional signaling pathway (LsTOR, LsS6K, and Ls4EBP), and reduced the phosphorylation of Akt. Conversely, injection of miR-9c-5p and novel-mir50 inhibitors induced the expressions of LsAkt, LsVg, LsVgR, LsTOR, LsS6K, and Ls4EBP, enhanced Akt phosphorylation level, and accelerated ovarian development. Injection of bovine insulin downregulated the expression of miR-9c-5p and novel-mir50 and upregulated the LsAkt expression. It also rescued the reproductive development defects associated with miR-9c-5p/novel-mir50 overexpression, forming a positive regulatory loop of insulin signaling. These results indicate that miR-9c-5p/novel-mir50 regulates the female reproduction of L. serricorne by targeting Akt in response to insulin signaling. The data also demonstrate the effects of the insulin/miRNA/Akt regulatory axis in insect reproduction.

MicroRNA Targets PAP1 to Mediate Melanization in Plutella xylostella (Linnaeus) Infected by Metarhizium anisopliae

MicroRNAs (miRNAs) play a pivotal role in important biological processes by regulating post-transcriptional gene expression and exhibit differential expression patterns during development, immune responses, and stress challenges. The diamondback moth causes significant economic damage to crops worldwide. Despite substantial advancements in understanding the molecular biology of this pest, our knowledge regarding the role of miRNAs in regulating key immunity-related genes remains limited. In this study, we leveraged whole transcriptome resequencing data from Plutella xylostella infected with Metarhizium anisopliae to identify specific miRNAs targeting the prophenoloxidase-activating protease1 (PAP1) gene and regulate phenoloxidase (PO) cascade during melanization. Seven miRNAs (pxy-miR-375-5p, pxy-miR-4448-3p, pxy-miR-279a-3p, pxy-miR-3286-3p, pxy-miR-965-5p, pxy-miR-8799-3p, and pxy-miR-14b-5p) were screened. Luciferase reporter assays confirmed that pxy-miR-279a-3p binds to the open reading frame (ORF) and pxy-miR-965-5p to the 3' untranslated region (3' UTR) of PAP1. Our experiments demonstrated that a pxy-miR-965-5p mimic significantly reduced PAP1 expression in P. xylostella larvae, suppressed PO activity, and increased larval mortality rate. Conversely, the injection of pxy-miR-965-5p inhibitor could increase PAP1 expression and PO activity while decreasing larval mortality rate. Furthermore, we identified four LncRNAs (MSTRG.32910.1, MSTRG.7100.1, MSTRG.6802.1, and MSTRG.22113.1) that potentially interact with pxy-miR-965-5p. Interference assays using antisense oligonucleotides (ASOs) revealed that silencing MSTRG.7100.1 and MSTRG.22113.1 increased the expression of pxy-miR-965-5p. These findings shed light on the potential role of pxy-miR-965-5p in the immune response of P. xylostella to M. anisopliae infection and provide a theoretical basis for biological control strategies targeting the immune system of this pest.

The miRNAs let-7b and miR-141 Coordinately Regulate Vitellogenesis by Modulating Methyl Farnesoate Degradation in the Swimming Crab Portunus trituberculatus

Methyl farnesoate (MF), a crucial sesquiterpenoid hormone, plays a pivotal role in the reproduction of female crustaceans, particularly in the vitellogenesis process. Despite extensive research on its functions, the molecular mechanisms that regulate MF levels during the vitellogenic phase remain largely elusive. This study investigates the roles of microRNAs (miRNAs), significant post-transcriptional regulators of gene expression, in controlling MF levels in the swimming crab Portunus trituberculatus. Through bioinformatic analysis, four miRNAs were identified as potential regulators targeting two genes encoding Carboxylesterases (CXEs), which are key enzymes in MF degradation. Dual luciferase reporter assays revealed that let-7b and miR-141 suppress CXE1 and CXE2 expression by directly binding to their 3' UTRs. In vivo overexpression of let-7b and miR-141 significantly diminished CXE1 and CXE2 levels, consequently elevating hemolymph MF and enhancing vitellogenin expression. Spatiotemporal expression profile analysis showed that these two miRNAs and their targets exhibited generally opposite patterns during ovarian development. These findings demonstrate that let-7b and miR-141 collaboratively modulate MF levels by targeting CXEs, thus influencing vitellogenesis in P. trituberculatus. Additionally, we found that the expression of let-7b and miR-141 were suppressed by MF, constituting a regulatory loop for the regulation of MF levels. The findings contribute novel insights into miRNA-mediated ovarian development regulation in crustaceans and offer valuable information for developing innovative reproduction manipulation techniques for P. trituberculatus.

The RNase III enzyme Dicer1 is essential for larval development in Bombyx mori

MicroRNAs (miRNAs) are important regulators of nearly all aspects of biological processes in eukaryotes. During the biogenesis of miRNAs, the RNase III enzyme Dicer processes double-strand precursor miRNAs into mature miRNAs and promotes the assembly of RNA-induced silencing complexes (RISCs). Dicer has been reported to participate in a wide range of physiological processes, including development and immunity, in some insect species. However, the physiological roles of Dicer in lepidopterans remain poorly understood. In this study, we investigated the function of Bombyx mori Dicer1. We first performed sequence alignment and found that the sequence of functional domains of Dicer1 are varied among Lepidoptera, Diptera, Coleoptera, Blattaria, and Orthoptera. Using a binary clustered regularly interspaced palindromic repeats (CRISPR) / CRISPR-associated protein 9 genome editing approach, we showed that BmDicer1 mutants have arrested development from the 3rd instar into the 4th instar. RNA sequencing analysis indicated that the defects in BmDicer1 mutants are due to dysregulation of genes that encode proteins involved in metabolism, protein degradation, absorption, and renin-angiotensin pathways. Analysis using quantitative real-time polymerase chain reaction showed that mutation of BmDicer1 altered expression of miRNAs and their target genes. Therefore, our study demonstrates the critical roles of BmDicer1 in miRNA biogenesis and larval development in silkworm.

Identification and Characterization of Development-Related microRNAs in the Red Flour Beetle, Tribolium castaneum

MicroRNAs (miRNAs) play important roles in insect growth and development, but they were poorly studied in insects. In this study, a total of 883 miRNAs were detected from the early embryo (EE), late larva (LL), early pupa (EP), late pupa (LP), and early adult (EA) of Tribolium castaneum by microarray assay. Further analysis identified 179 differentially expressed unique miRNAs (DEmiRNAs) during these developmental stages. Of the DEmiRNAs, 102 DEmiRNAs exhibited stage-specific expression patterns during development, including 53 specifically highly expressed miRNAs and 20 lowly expressed miRNAs in EE, 19 highly expressed miRNAs in LL, 5 weakly expressed miRNAs in EP, and 5 abundantly expressed miRNAs in EA. These miRNAs were predicted to target 747, 265, 472, 234, and 121 genes, respectively. GO enrichment analysis indicates that the targets were enriched by protein phosphorylation, calcium ion binding, sequence-specific DNA binding transcription factor activity, and cytoplasm. An RNA interference-mediated knockdown of the DEmiRNAs tca-miR-6-3p, tca-miR-9a-3p, tca-miR-9d-3p, tca-miR-11-3p, and tca-miR-13a-3p led to defects in metamorphosis and wing development of T. castaneum. This study has completed the identification and characterization of development-related miRNAs in T. castaneum, and will enable us to investigate their roles in the growth and development of insect.

MicroRNA miR-2765-3p regulates reproductive diapause by targeting FoxO in Galeruca daurica

The forkhead box O (FoxO), as a conserved transcription factor, plays an indispensable role in regulating insect diapause. However, how FoxO is regulated to control diapause in insects remains unknown. In this study, we discovered functional binding sites for miR-2765-3p in the 3' untranslated region of FoxO in Galeruca daurica. The luciferase reporter assay showed that miR-2765-3p targeted FoxO and suppressed its expression. The expression profiles of miR-2765-3p and FoxO displayed opposite patterns during the female developmental process. Overexpression of miR-2765-3p by the injection of the miR-2765-3p agomir into adult females reduced FoxO expression, leading to the suppression of lipid accumulation, promotion of ovarian development, and inhibition of reproductive diapause. This is similar to the phenotype that results from the depletion of FoxO by injecting dsFoxO into adult females. In addition, the repression of miR-2765-3p by injecting the miR-2765-3p antagomir increased the FoxO transcript level, leading to the stimulation of lipid accumulation, depression of ovarian development, and induction of reproductive diapause. A hormone injection assay showed that the juvenile hormone (JH) agonist (methoprene) upregulated miR-2765-3p and downregulated FoxO. Notably, injecting methoprene rescued ovarian development defects associated with miR-2765-3p inhibition. These findings indicate that the JH/miR-2765-3p/FoxO axis plays a vital role in the regulation of reproductive diapause in G. daurica.

The MicroRNA Ame-Bantam-3p Controls Larval Pupal Development by Targeting the Multiple Epidermal Growth Factor-like Domains 8 Gene (megf8) in the Honeybee, Apis mellifera

20-Hydroxyecdysone (20E) plays an essential role in coordinating developmental transitions in insects through responsive protein-coding genes and microRNAs (miRNAs). However, the interplay between 20E and miRNAs during insect metamorphosis is unknown. In this study, using small RNA sequencing, a comparative miRNA transcriptomic analysis in different development stages, and 20E treatment, we identified ame-bantam-3p as a key candidate miRNA involved in honeybee metamorphosis. Target prediction and in vitro dual-luciferase assays confirmed that ame-bantam-3p interacts with the coding region of the megf8 gene and promotes its expression. Meanwhile, temporal expression analysis revealed that the expression of ame-bantam-3p is higher in the larval stage than in prepupal and pupal stages, and that this expression pattern is similar to that of megf8. In vivo, we found that the mRNA level of megf8 was significantly increased after the injection of ame-bantam-3p agomir. A 20E feeding assay showed that 20E downregulated the expression of both ame-bantam-3p and its target gene megf8 on larval days five, six, and seven. Meanwhile, the injection of ame-bantam-3p agomir also reduced the 20E titer, as well as the transcript levels of essential ecdysteroid synthesis genes, including Dib, Phm, Sad, and Nvd. The transcript levels of 20E cascade genes, including EcRA, ECRB1, USP, E75, E93, and Br-c, were also significantly decreased after ame-bantam-3p agomir injection. However, ame-bantam-3p antagomir injection and dsmegf8 injection showed the opposite effect to ame-bantam-3p agomir injection. Ame-bantam-3p agomir treatment ultimately led to mortality and the failure of larval pupation by inhibiting ecdysteroid synthesis and the 20E signaling pathway. However, the expression of 20E signaling-related genes was significantly increased after megf8 knockdown, and larvae injected with dsmegf8 showed early pupation. Combined, our results indicate that ame-bantam-3p is involved in the 20E signaling pathway through positively regulating its target gene megf8 and is indispensable for larval–pupal development in the honeybee. These findings may enhance our understanding of the relationship between 20E signaling and small RNAs during honeybee development.

The variation of a cytochrome P450 gene, CYP6G4, drives the evolution of Musca domestica L. (Diptera: Muscidae) resistance to insecticides in China

Long term and excessive insecticide use has resulted in some environmental problems and especially, insecticide resistance evolution in insect pests. The variation of cytochrome P450 monooxygenases (P450s), associated with the metabolic detoxification of toxic xenobiotics, is often involved in insecticide resistance. Here, we found that the variation in a P450 gene, CYP6G4, is the most important driver of carbamates resistance in the house fly (Musca domestica). Deciphering the detailed molecular mechanisms of the insecticide resistance is critical for performing suitable insecticide resistance management strategies. Our research results revealed that the combination of amino acid mutations (110C-330E-360N/S, 110C-330E-360S) of CYP6G4 could improve the resistance to propoxur. The nucleotide variations in the promoter region of CYP6G4 significantly increased the luciferase activity by the reporter gene assays. Additionally, miR-281-1-5p was confirmed to post-transcriptionally down-regulate the expression of CYP6G4. These findings suggest that three independent mechanisms; amino acid mutations of the P450 protein, mutations in the promoter region and low expression of post-trans-regulatory factors, as the powerful strategies for the insect resistance to toxic compounds, play a crucial role in the evolutionary processes of insecticide resistance.

Responses of microRNA in digestive glands of mussel Mytilus galloprovincialis exposed to polystyrene nanoplastics

Polystyrene nanoplastics (PS-NPs) are typical accumulated nanoplastics in the marine environment and organisms, and have strong potential risks to marine ecological environment and human health. MiRNAs could respond to and participate in the response process of environmental stressors. However, the response of miRNAs to nanoplastics has not been fully explored. In this study, miRNA responses of digestive glands in mussels Mytilus galloprovincialis treated by 200 nm PS-NPs (20, 200, 2000 μg/L) for 7 days were characterized by BGISEQ-500 deep sequencing and bioinformatics analysis, along with histopathological quantification with planimetric parameters on hematoxylin and eosin (H&E) staining. Results showed that one novel miRNA (novel_mir63) and seven known miRNAs (miR-34_2, miR-34_5, miR-281_8, let-7-5p_6, miR-10, miR-124, miR-29b-3p) were significantly (adjusted P-value < 0.05) differentially expressed after PS-NPs treatments, and most of them were down-regulated expect for novel_mir63 and miR-34_2. Function analysis of target genes corresponding to these differentially expressed miRNAs indicated that PS-NPs disturbed the process related to metabolism, aging, cardiac function, neural excitation, and repairment. Among them, acetyl-CoA C-acetyltransferase and purine metabolism pathway played vital connection roles. Meanwhile, significantly morphology changes of digestive tubes obtained from H&E stained sections also implied severely disrupted metabolic capability in digestive glands, reflected by significantly increased mean diverticular radius (MDR) and mean luminal radius (MLR) values and the ratio of MLR to mean epithelial thickness (MET), and significantly decreased MET value and MET/MDR. Overall, these findings have revealed new characterization of miRNAs and their target genes in mussel M. galloprovincialis under PS-NPs stress, and provide important clues to further elucidate the toxicity mechanisms of PS-NPs.

miR-309a is a regulator of ovarian development in the oriental fruit fly Bactrocera dorsalis

Fecundity is arguably one of the most important life history traits, as it is closely tied to fitness. Most arthropods are recognized for their extreme reproductive capacity. For example, a single female of the oriental fruit fly Bactrocera dorsalis, a highly invasive species that is one of the most destructive agricultural pests worldwide, can lay more than 3000 eggs during its life span. The ovary is crucial for insect reproduction and its development requires further investigation at the molecular level. We report here that miR-309a is a regulator of ovarian development in B. dorsalis. Our bioinformatics and molecular studies have revealed that miR-309a binds the transcription factor pannier (GATA-binding factor A/pnr), and this activates yolk vitellogenin 2 (Vg 2) and vitellogenin receptor (VgR) advancing ovarian development. We further show that miR-309a is under the control of juvenile hormone (JH) and independent from 20-hydroxyecdysone. Thus, we identified a JH-controlled miR-309a/pnr axis that regulates Vg2 and VgR to control the ovarian development. This study has further enhanced our understanding of molecular mechanisms governing ovarian development and insect reproduction. It provides a background for identifying targets for controlling important Dipteran pests.

The ovary is a very critical organ for insect reproduction. Especially, many insect pests are famous for their large reproductive capacity. Therefore, understanding the molecular mechanisms involved in ovarian development could significantly contribute in the development of new insect pest control strategies. In this study, we report that miR-309a regulates the development of the ovary in an important dipteran pest, B. dorsalis, through a transcriptional factor, pannier (GATA-binding factor A/pnr), which in turn directly mediates the expression of yolk vitellogenin 2 (Vg 2) and vitellogenin receptor (VgR). Moreover, miR-309a is under the upstream control of juvenile hormone (JH). Here, in Dipterans, a novel JH-miR-309a-pnr-Vg-related genes regulatory pathway was found in ovarian development. This finding advances our understanding of a mechanism regulating insect ovarian development and provides new insights for potential targets to control dipteran pests through the reproductive strategy.

Identification and profiling of stable microRNAs in hemolymph of young and old Locusta migratoria fifth instars

Since the discovery of the first microRNA (miRNA) in the nematode Caenorhabditis elegans, numerous novel miRNAs have been identified which can regulate presumably every biological process in a wide range of metazoan species. In accordance, several insect miRNAs have been identified and functionally characterized. While regulatory RNA pathways are traditionally described at an intracellular level, studies reporting on the presence and potential role of extracellular (small) sRNAs have been emerging in the last decade, mainly in mammalian systems. Interestingly, evidence in several species indicates the functional transfer of extracellular RNAs between donor and recipient cells, illustrating RNA-based intercellular communication. In insects, however, reports on extracellular small RNAs are emerging but the number of detailed studies is still very limited. Here, we demonstrate the presence of stable sRNAs in the hemolymph of the migratory locust, Locusta migratoria. Moreover, the levels of several extracellular miRNAs (ex-miRNAs) present in locust hemolymph differed significantly between young and old fifth nymphal instars. In addition, we performed a 'proof of principle' experiment which suggested that extracellularly delivered miRNA molecules are capable of affecting the locusts' development.

A transcriptomic atlas underlying developmental plasticity of seasonal forms of Bicyclus anynana butterflies

Organisms residing in regions with alternating seasons often develop different phenotypes, or forms, in each season. These forms are often adaptations to each season and result from an altered developmental response to specific environmental cues such as temperature. While multiple studies have examined form-specific gene expression profiles in a diversity of species, little is known about how environments and developmental transitions, cued by hormone pulses, alter post-transcriptional patterns. In this study, we examine how gene expression, alternative splicing, and miRNA-mediated gene silencing in Bicyclus anynana butterfly hindwing tissue, varies across two rearing temperatures at four developmental timepoints. These timepoints flank two temperature-sensitive periods that coincide with two pulses of the insect hormone 20E. Our results suggest that developmental transitions, coincident with 20E pulses, elicit a greater impact on all these transcriptomic patterns than rearing temperatures per se. More similar transcriptomic patterns are observed pre-20E pulses than those observed post-20E pulses. We also found functionally distinct sets of differentially expressed and differentially spliced genes in the seasonal forms. Furthermore, around 10% of differentially expressed genes are predicted to be direct targets of, and regulated by, differentially expressed miRNAs between the seasonal forms. Many differentially expressed genes, miRNAs, or differentially spliced genes potentially regulate eyespot size plasticity, and we validated the differential splicing pattern of one such gene, daughterless. We present a comprehensive and interactive transcriptomic atlas of the hindwing tissue of both seasonal forms of B. anynana throughout development, a model organism of seasonal plasticity.

Transcriptome-wide identification of microRNAs in response to 20-hydroxyecdysone in Galeruca daurica

Both 20-hydroxyecdysone (20E) and miRNAs have multiple functions in the regulation of various physiological processes in insects. However, little is known about the interaction between 20E and miRNAs. In this study, six small RNA libraries were constructed from the adult Galeruca daurica treated with 20E and dimethyl sulfoxide (DMSO), respectively. Using small RNA sequencing, a total of 183 miRNAs, including 140 known and 43 novel miRNAs, were identified. Compared with the control (DMSO), 52 miRNAs (21 up-regulated and 31 down-regulated) were significantly differentially expressed after 20E treatment. The KEGG and GO analysis of the predicted genes targeted by 20E-responsive miRNAs indicate that 20E may influence the metabolic change during reproductive diapause in G. daurica via regulating miRNAs.

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.

Long Noncoding RNA lncR17454 Regulates Metamorphosis of Silkworm Through let-7 miRNA Cluster

A number of long noncoding RNAs (lncRNAs) have been identified in silkworm, but little is known about their functions. Recent study showed that the let-7 miRNA cluster (contains let-7, miR-2795, and miR-100) was transcribed from the last exon of lncRNA lncR17454 in silkworm. To investigate the functional role of lncR17454, dsRNAs of lncR17454 were injected into the hemolymph of 1-d-old third-instar larvae of Bombyx mori, repression of lncR17454 led to molting arrestment during the larval-larval and larval-pupal transition of silkworm, which was consistent to the result as let-7 knockdown in other studies. The expression level of mature let-7, miR-100, and miR-2795 decreased 40%, 36%, and 40%, respectively, while the mRNA level of two predicted target genes of let-7, the Broad Complex isoform 2 (BR-C-Z2) and the BTB-Zinc finger transcription repression factor gene Abrupt (Ab), increased significantly after lncR17454 knockdown. In contrast, when adding the 20-Hydroxyecdysone (20E) to silkworm BmN4 cell lines, the expression level of lncR17454 and let-7 cluster all increased significantly, but the expression of Abrupt, the predicted target gene of let-7, was repressed. Dual-luciferase reporter assays confirmed Abrupt was the real target of let-7. Here we found that the lncRNA lncR17454 can play regulator roles in the metamorphosis of silkworm through let-7 miRNA cluster and the ecdysone signaling pathway, which will provide new clues for lepidopteran pest control.

The miR-282-5p regulates larval moulting process by targeting chitinase 5 in Bombyx mori

Moulting is critical for growth, development and survival in insects. As the main components of cuticle, dynamic change of chitin is consistent with the moulting process. Chitinase is the main enzyme to mediate chitin metabolism in the old cuticle. To avoid over-degrading chitin from the new cuticle, the expression of chitinase must be precisely regulated. In this study, we performed microRNA-sequencing to investigate expression change of microRNAs in silkworm epidermis during the moulting process. A comparative microRNA transcriptomic analysis from different moulting stages and 20-hydroxyecdysone (20E) treatment identified bmo-miR-282-5p as a candidate. By the bioinformatic analysis, chitinase 5 (BmCht5) was predicted to be a target of bmo-miR-282-5p. Meanwhile, a temporal expression analysis revealed that BmCht5 only expressed at moulting D3 stage, whereas bmo-miR-282-5p showed a converse pattern, in which its transcript signal disappeared at this time point. Furthermore, a luciferase assay and agomir treatment demonstrated that bmo-miR-282-5p suppressed transcript of BmCht5 in vivo. As a result, injection of 282-5p agomir triggered 40% death due to moulting failure. In addition, RNA interference (RNAi)-mediated silencing of BmCht5 caused 30% developmental defect. Taken together, our data demonstrate the coordinated regulation of chitinase 5 by conserved miR-282-5p, and the 20E signalling pathway is essential for the normal moulting process in the domesticated silkworm.

De novo RNA-Seq and Annotation of Sesquiterpenoid and Ecdysteroid Biosynthesis Genes and MicroRNAs in a Spider Mite Eotetranychus kankitus

Eotetranychus kankitus is an important mite pest in citrus, but molecular data on the developmental processes of E. kankitus are lacking. The different development stages mix of E. kankitus was used to sequence for transcriptome and small RNAs to identify genes and predict miRNAs associated with sesquiterpenoid and ecdysteroid biosynthesis and signaling pathways. More than 36 million clean reads were assembled and 67,927 unigenes were generated. Of the unigenes, 19,300 were successfully annotated through annotation databases NR, SwissProt, COG, GO, KEGG, PFAM, and KOG. The transcripts were involved in sesquiterpenoid biosynthesis (11 genes) and ecdysteroid biosynthesis and signaling pathway (13 genes). Another, small RNA library was obtained and 31 conserved miRNAs were identified. Five most abundant miRNAs were Ek-miR-5735, Ek-miR-1, Ek-miR-263a, Ek-miR-184, and Ek-miR-8. The target genes related to sesquiterpenoid and ecdysteroid showed that 10 of the conserved miRNAs could potentially target the sesquiterpenoid and ecdysteroid pathway according to four-prediction software, sRNAT, miRanda, RNAhybrid, and Risearch2. Thus, the results of this study will provide bioinformatics information for further molecular studies of E. kankitus which may facilitate improved pest control strategies.

Spatio-Temporal Profiling of Metarhizium anisopliae-Responsive microRNAs Involved in Modulation of Plutella xylostella Immunity and Development

Metarhizium anisopliae, a ubiquitous pathogenic fungus, regulates a wide array of the insect pest population. The fungus has been employed to control Plutella xylostella, an insecticide-resistant destructive lepidopteran pest, which causes substantial economic losses in crops worldwide. Integration of modern gene-silencing technologies in pest control strategies has become more crucial to counter pesticide-resistant insects. MicroRNAs (miRNA) play essential roles in the various biological process via post-transcriptional gene regulation. In the present study, RNA-seq analysis of control (CK36h, CK72h) and fungal-infected (T36h, T72h) midguts was performed to reveal underlying molecular mechanisms occurring in larval midgut at different time courses. We aimed at exploring M. anisopliae-responsive miRNAs and their target genes involved in development and immunity. After data filtration, a combined set of 170 miRNAs were identified from all libraries. Interestingly, miR-281, miR-263, miR-1, miR-6094 and miR-8 were listed among the most abundantly expressed conserved miRNAs. Furthermore, we experimentally studied the role of differentially expressed miR-11912-5p in regulating corresponding target trypsin-like serine proteinase (Px_TLSP). The luciferase assay (in vitro) revealed that miRNA-11912-5p significantly downregulated its target gene, suggesting it might play a crucial role in defense mechanism of P. xylostella against M.+ anisopliae infection. We used synthetic miRNA mimic/inhibitor (in vivo), to overexpress/silence miRNA, which showed harmful effects on larval duration, survival and adult fecundity. Additionally, fungal application in the presence of mimics revealed enhanced sensitivity of P. xylostella to infection. Our finding provides an insight into the relatively obscure molecular mechanisms involved in insect midgut during the fungal infection.

Transcriptome dynamics during metamorphosis of imaginal discs into wings and thoracic dorsum in Apis mellifera castes

BACKGROUND

Much of the complex anatomy of a holometabolous insect is built from disc-shaped epithelial structures found inside the larva, i.e., the imaginal discs, which undergo a rapid differentiation during metamorphosis. Imaginal discs-derived structures, like wings, are built through the action of genes under precise regulation.

RESULTS

We analyzed 30 honeybee transcriptomes in the search for the gene expression needed for wings and thoracic dorsum construction from the larval wing discs primordia. Analyses were carried out before, during, and after the metamorphic molt and using worker and queen castes. Our RNA-seq libraries revealed 13,202 genes, representing 86.2% of the honeybee annotated genes. Gene Ontology analysis revealed functional terms that were caste-specific or shared by workers and queens. Genes expressed in wing discs and descendant structures showed differential expression profiles dynamics in premetamorphic, metamorphic and postmetamorphic developmental phases, and also between castes. At the metamorphic molt, when ecdysteroids peak, the wing buds of workers showed maximal gene upregulation comparatively to queens, thus underscoring differences in gene expression between castes at the height of the larval-pupal transition. Analysis of small RNA libraries of wing buds allowed us to build miRNA-mRNA interaction networks to predict the regulation of genes expressed during wing discs development.

CONCLUSION

Together, these data reveal gene expression dynamics leading to wings and thoracic dorsum formation from the wing discs, besides highlighting caste-specific differences during wing discs metamorphosis.

Genome-wide identification and characterization of long non-coding RNAs in Tribolium castaneum

Long non-coding RNAs (lncRNAs) are poorly understood in insects. In this study, we performed genome-wide analysis of lncRNAs in Tribolium castaneum by RNA-seq. In total, 4516 lncRNA transcripts corresponding to 3917 genes were identified from late embryos, early larvae, late larvae, early pupae, late pupae and early adults of T. castaneum, including 3152 novel lncRNAs and 1364 known lncRNAs. These lncRNAs have few exons and transcripts, and are short in length. During development, they exhibited nine different expression patterns. Functionally, they can act either by targeting messenger RNAs (1813 lncRNAs) and lncRNAs (45 lncRNAs) or as micro RNA (miRNA) precursors (46 lncRNAs). LncRNAs were observed to target the metabolic enzymes of glycolysis, TCA cycle and amino acids, demonstrating that lncRNAs control metabolism by regulating metabolic enzymes. Moreover, lncRNAs were shown to participate in cell differentiation and development via their targets. As miRNA precursors, lncRNAs could participate in the ecdysone signaling pathway. This study provides comprehensive information for lncRNAs of T. castaneum, and will promote functional analysis and target identification of lncRNAs in the insect.

Characterization and target gene analysis of microRNAs in the antennae of the parasitoid wasp Microplitis mediator

MicroRNAs (miRNAs), a class of non-coding single-strand RNA molecules encoded by endogenous genes, are about 21-24 nucleotides long and are involved in the post-transcriptional regulation of gene expression in plants and animals. Generally, the types and quantities of miRNAs in the different tissues of an organism are diverse, and these divergences may be related to their specific functions. Here we have identified 296 known miRNAs and 46 novel miRNAs in the antennae of the parasitoid wasp Microplitis mediator by high-throughput sequencing. Thirty-three miRNAs were predicted to target olfactory-associated genes, including odorant binding proteins (OBPs), chemosensory proteins, odorant receptors (ORs), ionotropic receptors (IRs) and gustatory receptors. Among these, 17 miRNAs were significantly highly expressed in the antennae, four miRNAs were highly expressed both in the antennae and head or wings, while the remaining 12 miRNAs were mainly expressed in the head, thorax, abdomen, legs and wings. Notably, miR-9a-5p and miR-2525-3p were highly expressed in male antennae, whereas miR-1000-5p and novel-miR-13 were enriched in female antennae. The 17 miRNAs highly expressed in antennae are likely to be associated with olfaction, and were predicted to target one OBP (targeted by miR-3751-3p), one IR (targeted by miR-7-5p) and 14 ORs (targeted by 15 miRNAs including miR-6-3p, miR-9a-5p, miR-9b-5p, miR-29-5p, miR-71-5p, miR-275-3p, miR-1000-5p, miR-1000-3p, miR-2525-3p, miR-6012-3p, miR-9719-3p, novel-miR-10, novel-miR-13, novel-miR-14 and novel-miR-28). These candidate olfactory-associated miRNAs are all likely to be involved in chemoreception through the regulation of chemosensory gene expression in the antennae of M. mediator.

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

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

Differential expression profiles of miRNAs and their putative targets in Schistosoma mansoni during its life cycle

BACKGROUND

Schistosomiasis is a disease caused by Schistosoma. Due to its complex life cycle, evolutionary position and sexual dimorphism, schistosomes have several mechanisms of gene regulation. MicroRNAs (miRNAs) are short endogenous RNAs that regulate gene expression at the post-transcriptional level by targeting mRNA transcripts.

OBJECTIVES

Here, we tested 12 miRNAs and identified their putative targets using a computational approach.

METHODS

We performed the expression profiles of a set of miRNAs and their putative targets during the parasite's life cycle by quantitative reverse transcription-polymerase chain reaction (qRT-PCR).

FINDINGS

Our results showed differential expression patterns of the mature miRNAs sma-miR-250; sma-miR-92a; sma-miR-new_4-3p; sma-miR-new_4-5p; sma-miR-new_5-5p; sma-miR-new_12-5p; sma-miR-new_13-3p and sma-miR-new_13-5p. Interestingly, many of the putative target genes are linked to oxidative phosphorylation and are up-regulated in adult-worms, which led us to suggest that miRNAs might play important roles in the post-transcriptional regulation of genes related to energetic metabolism inversion during parasite development. It is noteworthy that the expression of sma-miR-new_13-3p exhibited a negative correlation on SmNADH:ubiquinone oxidoreductase complex I.

MAIN CONCLUSIONS

Our analysis revealed putative miRNA genes related to important biological processes, such as transforming growth factor beta (TGF-β) signaling, proteasome regulation, glucose and lipid metabolism, immune system evasion and transcriptional regulation.

Expression profiling of Echinococcus multilocularis miRNAs throughout metacestode development in vitro

The neglected zoonotic disease alveolar echinococcosis (AE) is caused by the metacestode stage of the tapeworm parasite Echinococcus multilocularis. MicroRNAs (miRNAs) are small non-coding RNAs with a major role in regulating gene expression in key biological processes. We analyzed the expression profile of E. multilocularis miRNAs throughout metacestode development in vitro, determined the spatial expression of miR-71 in metacestodes cultured in vitro and predicted miRNA targets. Small cDNA libraries from different samples of E. multilocularis were sequenced. We confirmed the expression of 37 miRNAs in E. multilocularis being some of them absent in the host, such as miR-71. We found a few miRNAs highly expressed in all life cycle stages and conditions analyzed, whereas most miRNAs showed very low expression. The most expressed miRNAs were miR-71, miR-9, let-7, miR-10, miR-4989 and miR-1. The high expression of these miRNAs was conserved in other tapeworms, suggesting essential roles in development, survival, or host-parasite interaction. We found highly regulated miRNAs during the different transitions or cultured conditions analyzed, which might suggest a role in the regulation of developmental timing, host-parasite interaction, and/or in maintaining the unique developmental features of each developmental stage or condition. We determined that miR-71 is expressed in germinative cells and in other cell types of the germinal layer in E. multilocularis metacestodes cultured in vitro. MiRNA target prediction of the most highly expressed miRNAs and in silico functional analysis suggested conserved and essential roles for these miRNAs in parasite biology. We found relevant targets potentially involved in development, cell growth and death, lifespan regulation, transcription, signal transduction and cell motility. The evolutionary conservation and expression analyses of E. multilocularis miRNAs throughout metacestode development along with the in silico functional analyses of their predicted targets might help to identify selective therapeutic targets for treatment and control of AE.

Alveolar echinococcosis (AE) is a zoonotic disease caused by the metacestode stage of the helminth parasite Echinococcus multilocularis. Current treatment requires surgery and/or prolonged drug therapy. Thus, novel strategies for the treatment of AE are needed. MicroRNAs (miRNAs), a class of small ~22-nucleotide (nt) non-coding RNAs with a major role in regulating gene expression, have been suggested as potential therapeutic targets for treatment and control of helminth parasite infections. In this work, we analyzed the expression profile of E. multilocularis miRNAs throughout metacestode development in vitro. We predicted functional roles for highly expressed miRNAs and found that they could be involved in essential roles for survival and development in the host. We determined that E. multilocularis miR-71, a highly expressed miRNA that is absent in the human host, is expressed in germinative cells and in other cell types of the germinal layer in E. multilocularis metacestodes cultured in vitro. Germinative cells are a relevant cell type to target for anti-echinococcosis drug development. MiRNAs that are absent in the human host, involved in essential functions, highly expressed and/or expressed in germinative cells in E. multilocularis metacestodes may represent selective therapeutic targets for treatment and control of AE.

Hormonal regulation of microRNA expression dynamics in the gut of the yellow fever mosquito Aedes aegypti

The yellow fever mosquito Aedes aegypti is an obligatory blood feeder and a major arboviral disease vector, evoking severe public health concerns worldwide. In adult female mosquitoes, the gut is critical for blood digestion and pathogen entry. We aimed for a systematic exploration of microRNA expression dynamics in the gut during the gonadotrophic cycle. Small RNA libraries were constructed from female mosquito gut tissues at five time points. Unsupervised hierarchical clustering revealed three expression clusters (early, mid and late) peaking at sequential time points - 24, 48 and 72 h posteclosion. Differentially expressed miRNAs were identified at 24 h post-blood meal (PBM). Depletions of Methoprene-tolerant [Met; the juvenile hormone (JH) receptor] and Ecdysone receptor [EcR; the receptor to 20-hydroxyecdysone (20E)] were performed using dsRNA to these genes to investigate impacts on microRNA expressions. Our results suggest that Met-mediated signalling downregulates miRNA expression from the early cluster and upregulates that from the late cluster. EcR signalling either up- or downregulated miRNA levels at 24 h PBM, indicating a differential effect of this receptor in miRNA gene expression. Furthermore, miR-281, which is the most abundant miRNA in the gut tissue, is induced and repressed by Met- and EcR-mediated signalling, respectively. Systematic depletion using synthetic antagomir and phenotype examinations indicate that miR-281 is obligatory for the normal progression of blood digestion, ovarian development and reproduction. Collectively, this study unveils expression dynamics of microRNAs in the female gut tissue during the gonadotrophic cycle and demonstrates that they are affected by JH and 20E signalling.

Identification and profiling of Bactrocera dorsalis microRNAs and their potential roles in regulating the developmental transitions of egg hatching, molting, pupation and adult eclosion

MicroRNAs (miRNAs) are endogenous small noncoding RNAs (18-25 nt) that are involved in many physiological processes including development, cancer, immunity, apoptosis and host-microbe interactions through post-transcriptional regulation of gene expression. In this study, we measured the profile of small RNAs over the developmental transitions of the oriental fruit fly Bactrocera dorsalis from egg hatching, molting, and pupation to adult eclosion. We identified 250 miRNAs, including 83 known and 167 novel miRNAs, and 47 isomiRNAs. In addition, we identified the miRNAs differentially expressed over the developmental transitions. Interestingly, the miR-309 cluster, the miR-2 cluster/family and the let-7 cluster were among these differentially expressed miRNAs, suggesting a role in the regulation of egg hatching, molting and pupation/adult eclosion, respectively. Moreover, a detailed analysis of the temporal expression patterns of 14 highly expressed miRNAs in the pupal stage revealed three types of expression profiles. Furthermore, injection of a miR-100 mimic in the 3rd instar larvae resulted in a significant decrease in pupation and adult eclosion rates, whereas injection of a miR-317 antagomir resulted in a significant decrease in the pupation rate and a decrease in the pupation time, indicating that miR-100 and miR-317 are involved in the process of pupation. Finally, injection of a miR-100/miR-285 mimic or antagomir in pupae resulted in a significant decrease in the eclosion rate and a significant increase in the prevalence of a partial eclosion phenotype, implying the involvement of miR-100 and miR-285 in the process of adult eclosion. This study identified critical miRNAs involved in the transitions of this important holometabolic model and pest insect B. dorsalis from egg hatching to adult eclosion, thus providing a useful resource for exploring the regulatory role of miRNAs during insect post-embryonic development.

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.

lincRNA_Tc13743.2-miR-133-5p-TcGSTm02 regulation pathway mediates cyflumetofen resistance in Tetranychus cinnabarinus

Differential expression of metabolic detoxification enzymes is an important mechanism involved in pesticide/acaricide resistance of mite pests. The competing endogenous RNA hypothesis offers a new opportunity to investigate post-transcriptional regulation of those genes. In this study, 4454 long non-coding RNAs were identified in the carmine spider mite Tetranychus cinnabarinus by transcriptome sequencing. Software-based predictions indicated that a long intergenic non-coding RNA, (lincRNA)_Tc13743.2 and a detoxification enzyme gene, TcGSTm02, both contained a microRNA (miR-133-5p) response element. Over-expression of lincRNA_Tc13743.2 and TcGSTm02 were detected in a cyflumetofen-resistant T. cinnabarinus strain (CyR), whereas down-regulation of miR-133-5p was observed in this strain. Conversely, up-regulation of miR-133-5p could inhibit TcGSTm02 expression levels, and both lincRNA_Tc13743.2 and TcGSTm02 were significantly enriched in miR-133-5p biotin-avidin pull-down assays. RNA-binding protein immunoprecipitation assay showed that lincRNA_Tc13743.2 and TcGSTm02 bound to a silencing complex containing miR-133-5p. Moreover, a luciferase reporter assay based on a human cell line revealed that over-expression of lincRNA_Tc13743.2 could significantly reduce the inhibition exerted by miR-133-5p through the TcGSTm02 3'UTR. In addition, co-localization of lincRNA_Tc13743.2 and miR-133-5p was detected using fluorescence in situ hybridization, suggesting that lincRNA_Tc13743.2 interacts directly with miR-133-5p in spider mites. More importantly, silencing the expression of lincRNA_Tc13743.2 significantly reduced the expression levels of TcGSTm02 and increased the sensitivity of spider mites to cyflumetofen. Our data show that lincRNA_Tc13743.2 up-regulates TcGSTm02 expression by competing for miR-133-5p binding, demonstrating that a lincRNA_Tc13743.2-miR-133-5p-TcGSTm02 pathway mediates cyflumetofen resistance in mites.

MicroRNA-2738 regulates gene expression in the sex determination pathway in Bombyx mori

MicroRNAs (miRNAs) are a class of short, non-coding transcripts that bind to 3'-untranslated regions to trigger messenger RNA degradation or translational inhibition. Here we explored how miRNAs regulate sex determination in Bombyx mori, a lepidopteran model insect. Genes known to be involved in sex determination, BmPSI, Bmdsx, and BmMasc, are predicted targets of the species-specific miR-2738. Using a dual luciferase reporter assay in HEK293T cells, we confirmed that miR-2738 suppressed transcription of BmPSI, Bmdsx, and BmMasc. The levels of BmPSI and BmMasc were significantly down-regulated in B. mori miR-2738 overexpression. In contrast, the genetic disruption of miR-2738 using the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 transgenic system increased the levels of BmPSI and BmMasc transcripts, whereas splicing of Bmdsx was unaltered by miR-2738 depletion or overexpression. Taken together, this study implicates miR-2738 as a minor regulator of sex determination genes in the silkworm.

20-Hydroxyecdysone-responsive microRNAs of insects

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

miR-34 regulates larval growth and wing morphogenesis by directly modulating ecdysone signalling and cuticle protein in Bombyx mori

microRNAs (miRNA) are small non-coding RNAs that modulate the myriad biological activities by targeting genes, and many studies showed that miRNAs played a pivotal role in insect development. Here, we find that Bm-miRNA (miR-34) controls larval growth and wing morphology by targeting BmE74 and BmCPG4. Overexpression of miR-34 in the whole body caused a smaller body size, partially displays deformed wings and venation defects in adults. Ablation of miR-34 by transgenic CRISPR/Cas9 technology resulted in a severe developmental delay during the larval stage. Moreover, we confirmed that miR-34 directly targeted BmE74 and BmCPG4 by using a dual luciferase reporter assay in HEK293T cells. Remarkably, loss-of-function of BmCPG4 caused wing defects, which was similar to the phenotype of miR-34 overexpression in animals. In addition, our analysis revealed that ecdysone strongly inhibited miR-34 expression in vivo. Taken together, our study identifies miR-34 as a modulator that regulates larval growth and wing morphogenesis by directly modulating ecdysone signalling and cuticle protein in Bombyx mori.

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.

MiR-189942 regulates fufenozide susceptibility by modulating ecdysone receptor isoform B in Plutella xylostella (L.)

Although dibenzoylhydrazine-type non-steroidal ecdysone agonists, such as fufenozide, have an excellent performance record, the emergence of resistance could severely compromise the efficacy of these compounds in integrated pest management programs. To investigate possible mechanisms of resistance, we investigated the regulation of the expression of the PxEcR-B gene encoding the ecdysone receptor isoform B (PxEcR-B), which is the specific target of fufenozide in P. xylostella. Bioinformatics analysis revealed a putative miR-189942 binding site in the 3'-UTR of PxEcR-B mRNA. In a PxEcR-B 3'-UTR luciferase reporter system, miR-189942 downregulated the luciferase activity, and these effects were abolished by a deletion mutation in the putative miR-189942 binding site. Moreover, at 96 h after treatment with an agomir (mimic) or antagomir (inhibitor) of miR-189942, PxEcR-B expression was decreased by 71 ± 4% and increased by 4.19- fold respectively. Furthermore, overexpression or knockdown of miR-189942 changed the sensitivity of P. xylostella to fufenozide in vivo but had no influence on the sensitivity to chlorantraniliprole, which does not target PxEcR-B. These data indicate that miR-189942 suppressed PxEcR-B expression via binding at the 3'-UTR of PxEcR-B, thus increasing the tolerance of P. xylostella to fufenozide. These findings provide empirical evidence of the involvement of miRNAs in the regulation of insecticide resistance.

Analysis of miRNAs in the Heads of Different Castes of the Bumblebee Bombus lantschouensis (Hymenoptera: Apidae)

Bumblebees are important insect pollinators for many wildflowers and crops. MicroRNAs (miRNAs) are endogenous non-coding small RNAs that regulate different biological functions in insects. In this study, the miRNAs in the heads of the three castes of the bumblebee Bombus lantschouensis were identified and characterized by small RNA deep sequencing. The significant differences in the expression of miRNAs and their target genes were analyzed. The results showed that the length of the small RNA reads from males, queens, and workers was distributed between 18 and 30 nt, with a peak at 22 nt. A total of 364 known and 89 novel miRNAs were identified from the heads of the three castes. The eight miRNAs with the highest expressed levels in males, queens, and workers were identical, although the order of these miRNAs based on expression differed. The male vs. queen, male vs. worker, and worker vs. queen comparisons identified nine, fourteen, and four miRNAs with significant differences in expression, respectively. The different castes were clustered based on the differentially expressed miRNAs (DE miRNAs), and the expression levels of the DE miRNAs obtained by RT-qPCR were consistent with the read counts obtained through Solexa sequencing. The putative target genes of these DE miRNAs were enriched in 29 Gene Ontology (GO) terms, and catalytic activity was the most enriched GO term, as demonstrated by its association with 2837 target genes in the male vs. queen comparison, 3535 target genes in the male vs. worker comparison, and 2185 target genes in the worker vs. queen comparison. This study highlights the characteristics of the miRNAs in the three B. lantschouensis castes and will aid further studies on the functions of miRNAs in bumblebees.

Comparative analysis of shrimp (Penaeus vannamei) miRNAs expression profiles during WSSV infection under experimental conditions and in pond culture

In recent years, the importance of viral and host microRNAs (miRNAs) in mediating viral replication and control of host cellular machinery, has been realised and increasing efforts have been taken in order to understand the interactions of miRNAs from host and pathogen during infection. However, all existing studies has thus far been conducted in controlled experimental conditions and the veracity of these data for field conditions are yet to be established. In this framework, small RNA sequencing was performed to identify the miRNAs involved in shrimp (Penaeus vannamei) immune responses under two different WSSV infection conditions of natural infection and experimentally challenged conditions. The expression profiles of miRNAs of shrimp infected with WSSV under two contrasting conditions were compared and as a result, 23365 known miRNAs and 481 novel miRNAs were identified. Amongst the most abundantly expressed miRNAs, the hypoxia related miR-210 and immune pathway related miR-29b were expressed only in infected shrimps of both conditions. miR-8-5p, having a functional role in modulation of chitin biosynthesis was exclusively represented in higher numbers in the WSSV -infected shrimps under natural conditions whilst four of the miRNAs (mja-miR-6493-5p, mja-miR-6492, mmu-miR-3968, tcf-miR-9b-5p) identified from shrimps collected from pond culture targeted chitinase, an important enzyme involved in growth and moulting in shrimps, indicating an interaction between WSSV infection and moult cycle under culture conditions. Some of the miRNAs (tca-miR-87b-3p, cte-miR-277a) and miRNAs belonging to class miR-9, miR-981 that were identified only in WSSV infected shrimps under experimental conditions, are known to respond against WSSV infection in shrimps. Moreover, the miRNA target prediction revealed several immune-related gene targets such as cathepsin, c-type lectin, haemocyanin and ubiquitin protein ligase were commonly identified under both the conditions. However, the miRNAs identified from challenge experiment had wide number of gene targets as compared to the miRNAs of natural infection. The shrimp miRNA mja-miR-6489-3p, was also found to target early virus gene wsv001 of WSSV. Our study, therefore, provides the comparative analysis of miRNA expression from shrimp during WSSV infection in two different conditions.

microRNA-14 as an efficient suppressor to switch off ecdysone production after ecdysis in insects

The precise increase and decrease of hormone ecdysone are critical for accurate development in insects. Most previous works focus on transcriptional activation of ecdysone production; however, little is known about the mechanism of switching off ecdysone biosynthesis after ecdysis. Here, we showed that the precursor microRNA-14 (pre-miR-14) encodes two mature miRNAs in silkworm; both of these two mature miRNAs regulate various genes in the ecdysone-signalling pathway. Bmo-miR-14-5p targets on nine genes whereas Bmo-miR-14-3p targets on two genes in the same pathway. These two mature miRNAs increased immediately after the ecdysis, efficiently suppressing the 20-hydroxyecdysone (20E) biosynthesis, the upstream regulation, and the downstream response genes. Knocking down either of two mature miRNAs or both of them delays moult development, impairing development synchrony in antagomir-treated groups. In addition, overexpressing Bmo-miR-14-5p but not Bmo-miR-14-3p significantly affected the 20E titer and increased the moulting time variation, suggesting that Bmo-miR-14-5p, though it is less abundant, has more potent effects in development regulation than Bmo-miR-14-3p. In summary, we present evidence that a pre-miRNA encodes two mature miRNAs targeting on the same pathway, which significantly improves miRNA regulation efficiencies to programmatically switch off ecdysone biosynthesis.

MicroRNA Let-7 targets the ecdysone signaling pathway E75 gene to control larval-pupal development in Bactrocera dorsalis

MicroRNAs (miRNAs) regulate various biological processes during insect development; however, their role in larval-pupal development in oriental fruit fly, Bactrocera dorsalis (Hendel) remains unknown. In the current study, we address the biological function of a conserved miRNA, Bdo-Let-7 in the regulation of BdE75 gene, which belongs to the ecdysone signaling pathway and participates in the larval-pupal development in B. dorsalis. Using dual luciferase reporter assay in HEK293T cells we show that Bdo-Let-7 miRNA interacts with the 3' untranslated region of BdE75 gene and suppresses its expression. The Bdo-Let-7 and BdE75 are also co-expressed in the larval-pupal stages and in different tissues of B. dorsalis. In in vivo experiments, the injection of Bdo-Let-7 agomir and antagomir in third instar larvae down- and up-regulated the expression of BdE75, respectively. The 20-hydroxyecdysone (20E) injection assay shows that 20E up-regulated the expression of Bdo-Let-7 on the 5th day of the larvae. Moreover, abnormal pupation and eclosion were observed after larval Bdo-Let-7 antagomir injection. Based on these results, we show that Bdo-Let-7 regulates the ecdysone signaling pathway through the exact dose of BdE75 gene, and is indispensable for normal larval-pupal development in B. dorsalis.

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

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

Clustered miR-2, miR-13a, miR-13b and miR-71 coordinately target Notch gene to regulate oogenesis of the migratory locust Locusta migratoria

MicroRNAs (miRNAs), ∼22-nt small noncoding RNAs with a crucial role in various biological processes of organisms, are usually clustered in the genome. However, little is known about the miRNA clusters involved in insect reproduction. By small RNA sequencing and quantification followed by qRT-PCR, we found that the expression of invertebrate-specific miR-2/13/71 cluster including miR-2, miR-13a, miR-13b and miR-71 significantly decreased after adult ecdysis of the migratory locust, Locusta migratoria. Luciferase reporter assay and RNA immunoprecipitation demonstrated that miR-2/13/71 bound to the protein coding sequence of Notch and downregulated its expression. Injection of miR-2/13/71 agomiRs led to significant decrease of Notch expression as well as markedly reduced levels of Vitellogenin mRNA, suppressed oocyte maturation and impaired ovarian growth. Moreover, the expression of miR-2/13/71 was repressed by juvenile hormone (JH). Our results thus point to a previously unidentified mechanism by which JH-repressed miR-2/13/71 coordinately downregulates Notch to modulate insect reproduction. The increase of JH and decrease of miR-2/13/71 expression in both previtellogenic and vitellogenic stages of adult females ensure a high level of Notch expression, critically contributing to JH-dependent vitellogenesis and oogenesis.

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.

The microRNAs let-7 and miR-278 regulate insect metamorphosis and oogenesis by targeting the juvenile hormone early-response gene Krüppel-homolog 1

Krüppel-homolog 1 (Kr-h1), a zinc-finger transcription factor, inhibits larval metamorphosis and promotes adult reproduction by transducing juvenile hormone (JH). Although the transcriptional regulation of Kr-h1 has been extensively studied, little is known about its regulation at the post-transcriptional level. Using the migratory locust Locusta migratoria as a model system, we report here that the microRNAs let-7 and miR-278 bound to the Kr-h1 coding sequence and downregulated its expression. Application of let-7 and miR-278 mimics (agomiRs) significantly reduced the level of Kr-h1 transcripts, resulting in partially precocious metamorphosis in nymphs as well as markedly decreased yolk protein precursors, arrested ovarian development and blocked oocyte maturation in adults. Moreover, the expression of let-7 and miR-278 was repressed by JH, constituting a regulatory loop of JH signaling. This study thus reveals a previously unknown regulatory mechanism whereby JH suppresses the expression of let-7 and miR-278, which, together with JH induction of Kr-h1 transcription, prevents the precocious metamorphosis of nymphs and stimulates the reproduction of adult females. These results advance our understanding of the coordination of JH and miRNA regulation in insect development.

Juvenile hormone and sesquiterpenoids in arthropods: Biosynthesis, signaling, and role of MicroRNA

Arthropod molting and reproduction are precisely controlled by the levels of sesquiterpenoids, a class of C15 hormones derived from three isoprene units. The two major functional arthropod sesquiterpenoids are juvenile hormone (JH) and methyl farnesoate (MF). In hemimetabolous insects (such as the aphids, bugs, and cockroaches) and holometabolous insects (such as beetles, bees, butterflies, and flies), dramatic decrease in the titers of JH and/or MF promote metamorphosis from larvae to adults either directly or through an intermediate pupal stage, respectively. JH is absent in crustaceans (lobster, shrimp, crab) and other arthropods (chelicerates such as ticks, mites, spiders, scorpions and myriapods such as millipede and centipedes). In some crustaceans, molting and reproduction is dependent on changing levels of MF. The regulation of sesquiterpenoid production is thus crucial in the life cycle of arthropods. Dynamic and complex mechanisms have evolved to regulate sesquiterpenoid production. Noncoding RNAs such as the microRNAs are primary regulators. This article provides an overview of microRNAs that are known to regulate sesquiterpenoid production in arthropods.

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.

Parasitic insect-derived miRNAs modulate host development

Parasitic wasps produce several factors including venom, polydnaviruses (PDVs) and specialized wasp cells named teratocytes that benefit the survival of offspring by altering the physiology of hosts. However, the underlying molecular mechanisms for the alterations remain unclear. Here we find that the teratocytes of Cotesia vestalis, an endoparasitoid of the diamondback moth Plutella xylostella, and its associated bracovirus (CvBV) can produce miRNAs and deliver the products into the host via different ways. Certain miRNAs in the parasitized host are mainly produced by teratocytes, while the expression level of miRNAs encoded by CvBV can be 100-fold greater in parasitized hosts than non-parasitized ones. We further show that one teratocyte-produced miRNA (Cve-miR-281-3p) and one CvBV-produced miRNA (Cve-miR-novel22-5p-1) arrest host growth by modulating expression of the host ecdysone receptor (EcR). Altogether, our results show the first evidence of cross-species regulation by miRNAs in animal parasitism and their possible function in the alteration of host physiology during parasitism.

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

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