Analysis of plant-derived miRNAs in animal small RNA datasets

ABSTRACT: BACKGROUND: Plants contain significant quantities of small RNAs (sRNAs) derived from various sRNA biogenesis pathways. Many of these sRNAs play regulatory roles in plants. Previous analysis revealed that numerous sRNAs in corn, rice and soybean seeds have high sequence similarity to animal genes. However, exogenous RNA is considered to be unstable within the gastrointestinal tract of many animals, thus limiting potential for any adverse effects from consumption of dietary RNA. A recent paper reported that putative plant miRNAs were detected in animal plasma and serum, presumably acquired through ingestion, and may have a functional impact in the consuming organisms. RESULTS: To address the question of how common this phenomenon could be, we searched for plant miRNAs sequences in public sRNA datasets from various tissues of mammals, chicken and insects. Our analyses revealed that plant miRNAs were present in the animal sRNA datasets, and significantly miR168 was extremely over-represented. Furthermore, all or nearly all (>96%) miR168 sequences were monocot derived for most datasets, including datasets for two insects reared on dicot plants in their respective experiments. To investigate if plant-derived miRNAs, including miR168, could accumulate and move systemically in insects, we conducted insect feeding studies for three insects including corn rootworm, which has been shown to be responsive to plant-produced long double-stranded RNAs. CONCLUSIONS: Our analyses suggest that the observed plant miRNAs in animal sRNA datasets can originate in the process of sequencing, and that accumulation of plant miRNAs via dietary exposure is not universal in animals.     

A selective insecticidal protein from Pseudomonas mosselii for corn rootworm control

The coleopteran insect western corn rootworm (WCR, Diabrotica virgifera virgifera) is an economically important pest in North America and Europe. Transgenic corn plants producing Bacillus thuringiensis (Bt) insecticidal proteins have been useful against this devastating pest, but evolution of resistance has reduced their efficacy. Here, we report the discovery of a novel insecticidal protein, PIP‐47Aa, from an isolate of Pseudomonas mosselii. PIP‐47Aa sequence shows no shared motifs, domains or signatures with other known proteins. Recombinant PIP‐47Aa kills WCR, two other corn rootworm pests (Diabrotica barberi and Diabrotica undecimpunctata howardi) and two other beetle species (Diabrotica speciosa and Phyllotreta cruciferae), but it was not toxic to the spotted lady beetle (Coleomegilla maculata) or seven species of Lepidoptera and Hemiptera. Transgenic corn plants expressing PIP‐47Aa show significant protection from root damage by WCR. PIP‐47Aa kills a WCR strain resistant to mCry3A and does not share rootworm midgut binding sites with mCry3A or AfIP‐1A/1B from Alcaligenes that acts like Cry34Ab1/Cry35Ab1. Our results indicate that PIP‐47Aa is a novel insecticidal protein for controlling the corn rootworm pests.     

crw1 - A Novel Maize Mutant Highly Susceptible to Foliar Damage by the Western Corn Rootworm Beetle

Western corn rootworm (WCR), Diabrotica virgifera virgifera LeConte (Coleoptera: Chrysomelidae), is the most destructive insect pest of corn (Zea mays L.) in the United States. The adult WCR beetles derive their nourishment from multiple sources including corn pollen and silks as well as the pollen of alternate hosts. Conversely, the corn foliage is largely neglected as a food source by WCR beetles, leading to a perception of a passive interaction between the two. We report here a novel recessive mutation of corn that was identified and named after its foliar susceptibility to corn rootworm beetles (crw1). The crw1 mutant under field conditions was exceptionally susceptible to foliar damage by WCR beetles in an age-specific manner. It exhibits pleiotropic defects on cell wall biochemistry, morphology of leaf epidermal cells and lower structural integrity via differential accumulation of cell wall bound phenolic acids. These findings indicate that crw1 is perturbed in a pathway that was not previously ascribed to WCR susceptibility, as well as implying the presence of an active mechanism(s) deterring WCR beetles from devouring corn foliage. The discovery and characterization of this mutant provides a unique opportunity for genetic analysis of interactions between maize and adult WCR beetles and identify new strategies to control the spread and invasion of this destructive pest.     

Control of coleopteran insect pests through RNA interference

Commercial biotechnology solutions for controlling lepidopteran and coleopteran insect pests on crops depend on the expression of Bacillus thuringiensis insecticidal proteins, most of which permeabilize the membranes of gut epithelial cells of susceptible insects. However, insect control strategies involving a different mode of action would be valuable for managing the emergence of insect resistance. Toward this end, we demonstrate that ingestion of double-stranded (ds)RNAs supplied in an artificial diet triggers RNA interference in several coleopteran species, most notably the western corn rootworm (WCR) Diabrotica virgifera virgifera LeConte. This may result in larval stunting and mortality. Transgenic corn plants engineered to express WCR dsRNAs show a significant reduction in WCR feeding damage in a growth chamber assay, suggesting that the RNAi pathway can be exploited to control insect pests via in planta expression of a dsRNA.     

Systemic RNAi in western corn rootworm,Diabrotica virgifera virgifera,does not involve transitive pathways

Western corn rootworm (WCR, Diabrotica virgifera virgifera LeConte) is highly sensitive to orally delivered double‐stranded RNA (dsRNA). RNAi in WCR is systemic and spreads throughout the insect body. This raises the question whether transitive RNAi is a mechanism that functions in WCR to amplify the RNAi response via production of secondary siRNA. Secondary siRNA production is achieved through RNA‐dependent RNA polymerase (RdRP) activity in other eukaryotic organisms, but RdRP has not been identified in WCR and any other insects. This study visualized the spread of the RNAi‐mediated knockdown of Dv v‐ATPase C mRNA throughout the WCR gut and other tissues using high‐sensitivity branched DNA in situ hybridization. Furthermore, we did not detect either secondary siRNA production or transitive RNAi in WCR through siRNA sequence profile analysis. Nucleotide mismatched sequences introduced into either the sense or antisense strand of v‐ATPase C dsRNAs were maintained in siRNAs derived from WCR fed with the mismatched dsRNAs in a strand specific manner. The distribution of all siRNAs was restricted to within the original target sequence regions, which may indicate the lack of new dsRNA synthesis leading to production of secondary siRNA. Thus, the systemic spread of RNAi in WCR may be derived from the original dsRNA molecules taken up from the gut lumen. These results indicate that the initial dsRNA dose is important for a lethal systemic RNAi response in WCR and have implications in developing effective dsRNA traits to control WCR and in resistance management to prolong the durability of RNAi trait technology.     

Identification of novel microRNA-like-coding sites on the long-stem microRNA precursors in Arabidopsis

Plant microRNA (miRNA) is a crucial regulator of gene expression. It has been reported that more than one miRNA/miRNA* duplex could be produced from a microRNA precursor (pre-miRNA). In this study, we performed a comprehensive search for the novel miRNA candidates on the pre-miRNAs of Arabidopsis. AGO1 enrichment, co-existence of the miRNA*-like coordinates, and unique genome-wide match sites were taken into consideration for candidate screening. As a result, 43 miRNA-like candidates derived from 25 pre-miRNAs were identified. Among these candidates, 31 strong candidates from 22 pre-miRNAs passed all the filtering steps. Interestingly, some of these miRNA-like candidates showed organ-specific expression patterns. After target prediction and degradome sequencing data-based validation, five miRNA candidate–target pairs (ath-miR863-5p.2–AT1G76550.1, ath-miR822.2–AT5G03552.1, ath-miR822.3–AT5G02350.1, sRNA4–AT1G66290.1 and sRNA6–AT1G66310.1) were identified, providing a basis for in-depth functional analysis of these miRNA candidates. These results could update the current understanding of the biogenesis and the action of the plant miRNAs.     

SPORTS1.0: A Tool for Annotating and Profiling Non-coding RNAs Optimized for rRNA- and tRNA-derived Small RNAs

High-throughput RNA-seq has revolutionized the process of small RNA (sRNA) discovery, leading to a rapid expansion of sRNA categories. In addition to the previously well-characterized sRNAs such as microRNAs (miRNAs), piwi-interacting RNAs (piRNAs), and small nucleolar RNA (snoRNAs), recent emerging studies have spotlighted on tRNA-derived sRNAs (tsRNAs) and rRNA-derived sRNAs (rsRNAs) as new categories of sRNAs that bear versatile functions. Since existing software and pipelines for sRNA annotation are mostly focused on analyzing miRNAs or piRNAs, here we developed the sRNA annotation pipelineoptimized for rRNA- and tRNA-derived sRNAs (SPORTS1.0). SPORTS1.0 is optimized for analyzing tsRNAs and rsRNAs from sRNA-seq data, in addition to its capacity to annotate canonical sRNAs such as miRNAs and piRNAs. Moreover, SPORTS1.0 can predict potential RNA modification sites based on nucleotide mismatches within sRNAs. SPORTS1.0 is precompiled to annotate sRNAs for a wide range of 68 species across bacteria, yeast, plant, and animal kingdoms, while additional species for analyses could be readily expanded upon end users’ input. For demonstration, by analyzing sRNA datasets using SPORTS1.0, we reveal that distinct signatures are present in tsRNAs and rsRNAs from different mouse cell types. We also find that compared to other sRNA species, tsRNAs bear the highest mismatch rate, which is consistent with their highly modified nature. SPORTS1.0 is an open-source software and can be publically accessed at https://github.com/junchaoshi/sports1.0.     

Abnormally high digestive enzyme activity and gene expression explain the contemporary evolution of a Diabrotica biotype able to feed on soybeans

Western corn rootworm (Diabrotica virgifera) (WCR) depends on the continuous availability of corn. Broad adoption of annual crop rotation between corn (Zea mays) and nonhost soybean (Glycine max) exploited WCR biology to provide excellent WCR control, but this practice dramatically reduced landscape heterogeneity in East-central Illinois and imposed intense selection pressure. This selection resulted in behavioral changes and "rotation-resistant" (RR) WCR adults. Although soybeans are well defended against Coleopteran insects by cysteine protease inhibitors, RR-WCR feed on soybean foliage and remain long enough to deposit eggs that will hatch the following spring and larvae will feed on roots of planted corn. Other than documenting changes in insect mobility and egg laying behavior, 15 years of research have failed to identify any diagnostic differences between wild-type (WT)- and RR-WCR or a mechanism that allows for prolonged RR-WCR feeding and survival in soybean fields. We documented differences in behavior, physiology, digestive protease activity (threefold to fourfold increases), and protease gene expression in the gut of RR-WCR adults. Our data suggest that higher constitutive activity levels of cathepsin L are part of the mechanism that enables populations of WCR to circumvent soybean defenses, and thus, crop rotation. These new insights into the mechanism of WCR tolerance of soybean herbivory transcend the issue of RR-WCR diagnostics and management to link changes in insect gut proteolytic activity and behavior with landscape heterogeneity. The RR-WCR illustrates how agro-ecological factors can affect the evolution of insects in human-altered ecosystems.     

Identification of MicroRNAs in Helicoverpa armigera and Spodoptera litura Based on Deep Sequencing and Homology Analysis

The current identification of microRNAs (miRNAs) in insects is largely dependent on genome sequences. However, the lack of available genome sequences inhibits the identification of miRNAs in various insect species. In this study, we used a miRNA database of the silkworm Bombyx mori as a reference to identify miRNAs in Helicoverpa armigera and Spodoptera litura using deep sequencing and homology analysis. Because all three species belong to the Lepidoptera, the experiment produced reliable results. Our study identified 97 and 91 conserved miRNAs in H. armigera and S. litura, respectively. Using the genome of B. mori and BAC sequences of H. armigera as references, 1 novel miRNA and 8 novel miRNA candidates were identified in H. armigera, and 4 novel miRNA candidates were identified in S. litura. An evolutionary analysis revealed that most of the identified miRNAs were insect-specific, and more than 20 miRNAs were Lepidoptera-specific. The investigation of the expression patterns of miR-2a, miR-34, miR-2796-3p and miR-11 revealed their potential roles in insect development. miRNA target prediction revealed that conserved miRNA target sites exist in various genes in the 3 species. Conserved miRNA target sites for the Hsp90 gene among the 3 species were validated in the mammalian 293T cell line using a dual-luciferase reporter assay. Our study provides a new approach with which to identify miRNAs in insects lacking genome information and contributes to the functional analysis of insect miRNAs.     

Identification and Characteristics of microRNAs from Army Worm,Spodoptera frugiperda Cell Line Sf21

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

Ultrastructural Changes Caused by Snf7 RNAi in Larval Enterocytes of Western Corn Rootworm (Diabrotica virgifera virgifera Le Conte)

The high sensitivity to oral RNA interference (RNAi) of western corn rootworm (WCR, Diabrotica virgifera virgifera Le Conte) provides a novel tool for pest control. Previous studies have shown that RNAi of DvSnf7, an essential cellular component of endosomal sorting complex required for transport (ESCRT), caused deficiencies in protein de-ubiquitination and autophagy, leading to WCR death. Here we investigated the detailed mechanism leading to larval death by analyzing the ultrastructural changes in midgut enterocytes of WCR treated with double-stranded RNA (ds-DvSnf7). The progressive phases of pathological symptoms caused by DvSnf7-RNAi in enterocytes include: 1) the appearance of irregularly shaped macroautophagic complexes consisting of relatively large lysosomes and multi-lamellar bodies, indicative of failure in autolysosome formation; 2) cell sloughing and loss of apical microvilli, and eventually, 3) massive loss of cellular contents indicating loss of membrane integrity. These data suggest that the critical functions of Snf7 in insect midgut cells demonstrated by the ultrastructural changes in DvSnf7 larval enterocytes underlies the conserved essential function of the ESCRT pathway in autophagy and membrane stability in other organisms.