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

Importance of core microRNA pathway genes and microRNAs associated with the defense of Odontotermes formosanus (Shiraki) against Serratia marcescens infection

MicroRNAs (miRNAs) are noncoding small regulatory RNAs involved in diverse biological processes. Odontotermes formosanus (Shiraki) is a polyphagous pest that causes economic damage to agroforestry. Serratia marcescens is a bacterium with great potential for controlling this insect. However, knowledge about the miRNA pathway and the role of miRNAs in O. formosanus defense against SM1 is limited. In this study, OfAgo1, OfDicer1 and OfDrosha were differentially expressed in different castes and tissues. SM1 infection affected the expression of all three genes in O. formosanus. Then, we used specific double-stranded RNAs to silence OfAgo1, OfDicer1 and OfDrosha. Knockdown of these genes enhanced the virulence of SM1 to O. formosanus, suggesting that miRNAs were critical in the defense of O. formosanus against SM1. Furthermore, we sequenced miRNAs from SM1-infected and uninfected O. formosanus. 33 differentially expressed (DE) miRNAs were identified, whereby 22 were upregulated and 11 were downregulated. Finally, the miRNA-mRNA networks were constructed, which further suggested the important role of miRNAs in the defense of O. formosanus against SM1. Totally, O. formosanus miRNA core genes defend against SM1 infection by regulating miRNA expression. This study elucidates the interactions between O. formosanus and SM1 and provides new theories for biological control.

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.

miR-92a-1-p5 Modulated Expression of the flightin Gene Regulates Flight Muscle Formation and Wing Extension in the Pea Aphid, Acyrthosiphon pisum (Hemiptera: Aphidoidea)

Aphids exhibit wing polyphenism. Winged and wingless aphid morphs are produced by parthenogenesis depending on population density and host plant quality. Recent studies showed that microRNAs in alate and apterous individuals have differential expression and are involved in wing dimorphism of Acyrthosiphon pisum. From which miR-92a-1-p5 can target the mRNA of flight muscle gene flightin in vitro, but what effect they have on wing development of aphid is unclear. Here with the nanocarrier-delivered RNA interference (RNAi) method, flightin gene was knocked down in winged nymphs of A. pisum. Results showed that the majority (63.33%) of adults had malformed wings, the shape of dorsal longitudinal muscle (DLM) was deformed severely, the dorsoventral flight muscle (DVM) became wider and looser in aphids with flightin reduction compared with the negative control. Overexpression of miR-92a-1-p5 caused decreased expression of flightin and malformed wings of aphids, with a mutant ratio of 62.50%. Morphological analysis of flight musculature showed the consistent result as that with flightin knockdown. These results suggest that flightin is essential for flight musculature formation and wing extension in A. pisum, which can be modulated by miR-92a-1-p5.

Temporal analysis of microRNAs associated with wing development in the English grain aphid, Sitobion avenae (F.) (Homoptera: Aphidiae)

Molecular mechanisms underlying wing evolution and development have been a point of scientific inquiry for decades. Phloem-feeding aphids are one of the most devastating global insect pests, where dispersal of winged morphs lead to annual movements, migrations, and range expansions. Aphids show a polyphenic wing dimorphism trait, and offer a model to study the role of environment in determining morphological plasticity of a single genotype. Despite recent progresses in the genetic understanding of wing polyphenism, the influence of environmental cues remains unclear. To investigate the involvement of miRNAs in wing development, we sequenced small RNA libraries of the English grain aphid, Sitobion avenae (F.) across six different developmental stages. As a result, we identified 113 conserved and 193 S. avenae-specific miRNAs. Gene Ontology and KEGG pathway analyses of putative target mRNAs for the six differentially expressed miRNAs are enriched for wing development processes. Dietary uptake of miR-263a, miR-316, and miR-184a agomirs and antagomirs led to significantly higher mortality (>70%) and a lower proportion of winged morphs (<5%). On the other hand, wing malformation was observed in miR-2 and miR-306 agomirs and miR-2 and miR-14 antagomirs, respectively, suggesting their involvement in S. avenae wing morphogenesis. These combined results not only shed light on the regulatory role of miRNAs in wing dimorphism, but also provide potential novel targets for the long-term sustainable management of S. avenae, a devastating global grain pest.

An InR/mir-9a/NlUbx regulatory cascade regulates wing diphenism in brown planthoppers

Wing polymorphism significantly contributes to the ecological success of some insect species. For example, the brown planthopper (BPH) Nilaparvata lugens, which is one of the most destructive rice pests in Asia, can develop into either highly mobile long-winged or highly fecund short-winged adult morphs. A recent study reported a highly provocative result that the Hox gene Ultrabithorax (Ubx) is expressed in BPH forewings and showed that this wing development gene is differentially expressed in nymphs that develop into long-winged versus short-winged morphs. Here, we found that Ubx may be a mir-9a target, and used dual luciferase reporter assays and injected micro RNA (miRNA) mimics and inhibitors to confirm the interactions between mir-9a and NlUbx. We measured the mir-9a and NlUbx expression profiles in nymphs and found that the expression of these two biomolecules was negatively correlated. By rearing BPH nymphs on host rice plants with different nutritional status, we were able to characterize a regulatory cascade between insulin receptor genes, mir-9a, and NlUbx that regulate wing length in BPHs. When host quality was low, NlInR1 expression in the nymph terga increased and NlInR2 expression decreased; this led to a higher mir-9a level, which in turn reduced the NlUbx transcript level and ultimately resulted in longer wing lengths. Beyond extending our understanding of the interplay between host plant status and genetic events that modulate polymorphism, we demonstrated both the upstream signal and miRNA-based regulatory mechanism that control Ubx expression in BPH forewings.

Identification of differentially expressed microRNAs under imidacloprid exposure in Sitobion miscanthi

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

Male-Biased microRNA Discovery in the Pea Aphid

Epigenetic mechanisms modulate gene expression levels during development, shaping how a single genome produces a diversity of phenotypes. Here, we begin to explore the epigenetic regulation of sexual dimorphism in pea aphids (Acyrthosiphon pisum) by focusing on microRNAs. Previous analyses of microRNAs in aphids have focused solely on females, so we performed deep sequencing of a sample containing early-stage males. We used this sample, plus samples from Genbank, to find 207 novel pea aphid microRNA coding loci. We localized microRNA loci to a chromosome-level assembly of the pea aphid genome and found that those on the X chromosome have lower overall expression compared to those on autosomes. We then identified a set of 19 putative male-biased microRNAs and found them enriched on the X chromosome. Finally, we performed protein-coding RNA-Seq of first instar female and male pea aphids to identify genes with lower expression in males. 10 of these genes were predicted targets of the 19 male-biased microRNAs. Our study provides the most complete set of microRNAs in the pea aphid to date and serves as foundational work for future studies on the epigenetic control of sexual dimorphism.

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

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

Genome-Wide Identification and Analysis of Nilaparvata lugens microRNAs during Challenge with the Entomopathogenic Fungus Metarhizium anisopliae

The resistance of the notorious rice pest Nilaparvata lugens to many insecticides has caused significant concerns. Our previous study demonstrated that the fungus Metarhizium anisopliae CQMa421 shows great potential for the control of this pest, but the interactions between them are still unclear. Thus, we further investigated fungal infection-related microRNAs (miRNAs) in N. lugens during M. anisopliae CQMa421 challenge using Illumina sequencing. In this study, we constructed twenty-four small RNA libraries over different time courses (i.e., 4 h, 8 h, 16 h, and 24 h). A total of 478.62 M clean reads were collected, with each sample producing more than 13.37 M reads, after the removal of low-quality reads. We identified 2324 miRNAs and their 11,076 target genes within the twenty-four libraries by bioinformatics analysis. Differentially expressed miRNAs (DEmiRNAs), including 58 (32 upregulated vs. 26 downregulated), 62 (30 upregulated vs. 32 downregulated), 126 (71 upregulated vs. 55 downregulated), and 109 (40 upregulated vs. 69 downregulated) DEmiRNAs were identified at 4 h, 8 h, 16 h, and 24 h post-infection, respectively. We further conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis to predict the functions of all target genes of DEmiRNAs. These DEmiRNAs targets identified during 24 h of infection were primarily involved in energy metabolism, lysine degradation, the FoxO signaling pathway, ubiquitin-mediated proteolysis, the mRNA surveillance pathway, and the MAPK signaling pathway. Taken together, our results provide essential information for further study of the interactions between the entomopathogenic fungus M. anisopliae and N. lugens at the posttranscriptional level.

Exploiting somatic piRNAs in Bemisia tabaci enables novel gene silencing through RNA feeding

RNAi usually relies on Dicer-produced siRNAs to induce gene silencing. In many arthropods, another type of RNAi is present in the soma—the piRNA pathway. This work finds exploiting this biology is a viable alternative for gene knockdown.

RNAi promises to reshape pest control by being nontoxic, biodegradable, and species specific. However, due to the plastic nature of RNAi, there is a significant variability in responses. In this study, we investigate small RNA pathways and processing of ingested RNAi trigger molecules in a hemipteran plant pest, the whitefly Bemisia tabaci. Unlike Drosophila, where the paradigm for insect RNAi technology was established, whitefly has abundant somatic piwi-associated RNAs (piRNAs). Long regarded as germline restricted, piRNAs are common in the soma of many invertebrates. We sought to exploit this for a novel gene silencing approach. The main principle of piRNA biogenesis is the recruitment of target RNA fragments into the pathway. As such, we designed synthetic RNAs to possess complementarity to the loci we annotated. Following feeding of these exogenous piRNA triggers knockdown as effective as conventional siRNA-only approaches was observed. These results demonstrate a new approach for RNAi technology that could be applicable to dsRNA-recalcitrant pest species and could be fundamental to realizing insecticidal RNAi against pests.

MicroRNA expression profiling of Plutella xylostella after challenge with B. thuringiensis

The diamondback moth, Plutella xylostella, the main pest of brassica crops, has developed resistance to almost all major classes of insecticides as the farmers rely on insecticides to control this pest. An extensive use of broad-spectrum insecticides against P. xylostella promotes the selection of insecticide resistance, destroy natural enemies, and pollute the environment. In this scenario, it is imperative to use genetic methods such as gene silencing technology as an alternate approach against this pest. Evidence shows that microRNAs play pivotal roles in the regulation of target genes at the post-transcription level and show differential expression under various biological processes. However, the knowledge of their role in insect immunity is still in its infancy. In the present study, we aimed at exploring the response of P. xylostella miRNAs against B. thuringiensis at different time courses (6, 12, 18, 24, and 36 h) by using small RNA sequencing. After data filtration, a combined set of 149 miRNAs was identified from all the libraries. Interestingly, a couple of conserved miRNAs such as miR-1, Let-7, miR-275, miR-184, and miR-10 were listed as abundantly expressed miRNAs after exposure to B. thuringiensis. It is worth mentioning that the differential expression analysis revealed that miR-2, a conserved miRNA, was up-regulated following infection. Furthermore, we experimentally validated the involvement of miR-2b-3p in the regulation of corresponding target trypsin. Our luciferase assay results revealed that miR-2b-3p mimic significantly down-regulated the target gene trypsin indicating that it might play a crucial role in the defense mechanism of P. xylostella against B. thuringiensis infection. On the whole, our findings provide insights into the possible regulatory role of miRNAs in insect immunity in response to microorganisms.

Molecular Mechanisms of Wing Polymorphism in Insects

Many insects are capable of developing into either long-winged or short-winged (or wingless) morphs, which enables them to rapidly match heterogeneous environments. Thus, the wing polymorphism is an adaptation at the root of their ecological success. Wing polymorphism is orchestrated at various levels, starting with the insect's perception of environmental cues, then signal transduction and signal execution, and ultimately the transmitting of signals into physiological adaption in accordance with the particular morph produced. Juvenile hormone and ecdysteroid pathways have long been proposed to regulate wing polymorphism in insects, but rigorous experimental evidence is lacking. The breakthrough findings of ecdysone receptor regulation on transgenerational wing dimorphism in the aphid Acyrthosiphon pisum and of insulin signaling in the planthopper Nilaparvata lugens greatly broaden our understanding of wing polymorphism at the molecular level. Recently, the advent of high-throughput sequencing coupled with functional genomics provides powerful genetic tools for future insights into the molecular bases underlying wing polymorphism in insects.

microRNA regulation in an ancient obligate endosymbiosis

Although many insects are associated with obligate bacterial endosymbionts, the mechanisms by which these host/endosymbiont associations are regulated remain mysterious. While microRNAs (miRNAs) have been recently identified as regulators of host/microbe interactions, including host/pathogen and host/facultative endosymbiont interactions, the role miRNAs may play in mediating host/obligate endosymbiont interactions is virtually unknown. Here, we identified conserved miRNAs that potentially mediate symbiotic interactions between aphids and their obligate endosymbiont, Buchnera aphidicola. Using small RNA sequence data from Myzus persicae and Acyrthosiphon pisum, we annotated 93 M. persicae and 89 A. pisum miRNAs, among which 69 were shared. We found 14 miRNAs that were either highly expressed in aphid bacteriome, the Buchnera-housing tissue, or differentially expressed in bacteriome vs. gut, a non-Buchnera-housing tissue. Strikingly, 10 of these 14 miRNAs have been implicated previously in other host/microbe interaction studies. Investigating the interaction networks of these miRNAs using a custom computational pipeline, we identified 103 miRNA::mRNA interactions shared between M. persicae and A. pisum. Functional annotation of the shared mRNA targets revealed only two over-represented cluster of orthologous group categories: amino acid transport and metabolism, and signal transduction mechanisms. Our work supports a role for miRNAs in mediating host/symbiont interactions between aphids and their obligate endosymbiont Buchnera. In addition, our results highlight the probable importance of signal transduction mechanisms to host/endosymbiont coevolution.

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

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

Decoding the architecture and origins of mechanisms for developmental polyphenism

Developmental polyphenism affords a single genotype multiple solutions to match an organism to its environment. Because polyphenism is the extreme example of how development deviates from a linear genetic blueprint, it demands a genetic explanation for how environmental cues shunt development to hypothetically alternative modules. We highlight several recent advances that have begun to illuminate genetic mechanisms for polyphenism and how this recurring developmental novelty may arise. An emerging genetic knowledge of polyphenism is providing precise targets for testing hypotheses of how switch mechanisms are built-out of olfactory, nutrient-sensing, hormone-reception, and developmental and genetic buffering systems-to accommodate plasticity. Moreover, classic and new model systems are testing the genetic basis of polyphenism's proposed causal roles in evolutionary change.

Zika virus alters the microRNA expression profile and elicits an RNAi response in Aedes aegypti mosquitoes

Zika virus (ZIKV), a flavivirus transmitted primarily by Aedes aegypti, has recently spread globally in an unprecedented fashion, yet we have a poor understanding of host-microbe interactions in this system. To gain insights into the interplay between ZIKV and the mosquito, we sequenced the small RNA profiles in ZIKV-infected and non-infected Ae. aegypti mosquitoes at 2, 7 and 14 days post-infection. ZIKA induced an RNAi response in the mosquito with virus-derived short interfering RNAs and PIWI-interacting RNAs dramatically increased in abundance post-infection. Further, we found 17 host microRNAs (miRNAs) that were modulated by ZIKV infection at all time points. Strikingly, many of these regulated miRNAs have been reported to have their expression altered by dengue and West Nile viruses, while the response was divergent from that induced by the alphavirus Chikungunya virus in mosquitoes. This suggests that conserved miRNA responses occur within mosquitoes in response to flavivirus infection. This study expands our understanding of ZIKV-vector interactions and provides potential avenues to be further investigated to target ZIKV in the mosquito host.