An analysis of the small RNA transcriptome of four developmental stages of the citrus red mite (Panonychus citri)

Genome-wide analysis of Panonychus citri microRNAs with a focus on potential insecticidal activity of 4 microRNAs to eggs and nymphs

Panonychus citri McGregor (Acari: Tetranychidae), a destructive citrus pest, causes considerable annual economic losses due to its short lifespan and rapid resistance development. MicroRNA (miRNA)-induced RNA interference is a promising approach for pest control because of endogenous regulation of pest growth and development. To search for miRNAs with potential insecticidal activity in P. citri, genome-wide analysis of miRNAs at different developmental stages was conducted, resulting in the identification of 136 miRNAs, including 73 known and 63 novel miRNAs. A total of 17 isomiRNAs and 12 duplicated miRNAs were characterized. MiR-1 and miR-252-5p were identified as reference miRNAs for P. citri and Tetranychus urticae. Based on differential expression analysis, treatments with miR-let-7a and miR-315 mimics and the miR-let-7a antagomir significantly reduced the egg hatch rate and resulted in abnormal egg development. Overexpression or downregulation of miR-34-5p and miR-305-5p through feeding significantly decreased the adult eclosion rate and caused molting defects. The 4 miRNAs, miR-let-7a, miR-315, miR-34-5p, and miR-305-5p, had important regulatory functions and insecticidal properties in egg hatching and adult eclosion. In general, these data advance our understanding of miRNAs in mite biology, which can assist future studies on insect-specific miRNA-based green pest control technology.

Regulatory roles of three miRNAs on allergen mRNA expression in Tyrophagus putrescentiae

BACKGROUND

Tyrophagus putresecentiae is an important mite species in rural and urban environments, causing sensitization and allergic disease. While evidence suggests that microRNAs (miRNAs) may regulate the expression of allergen-encoding genes, no study has directly investigated this possibility. Here, this gap was addressed by profiling miRNAs and elucidating their target allergen messenger RNAs (mRNAs) in this mite species.

METHODS

Small RNA and transcriptome libraries were constructed for eggs, larvae, nymphs, and adults. After deep miRNA and whole-transcriptome sequencing were performed, the miRNA and allergen-encoding mRNA regulatory networks were explored.

RESULTS

A total of 540 miRNAs were identified, including 155 with expression levels differing significantly across the four mite developmental stages (p < .01), 59 of which were novel. The mRNA expression for allergens was higher for Tyr p 1 in adults than in other developmental stages; Tyr p 2-5, 7, 10, 13, 33, and 34 in immature stages; and Tyr p 28, 35, and 36 in eggs and adults. A combined miRNA and transcriptome bioinformatics analysis showed that allergen Tyr p 3 was regulated by miRNA PC-5p-5698441_1, Tyr p 4 was regulated by PC-5p-7050653_1, and Tyr p 34 was regulated by PC-5p-5534223_1 and PC-5p-5698441_1. These three allergen mRNA and three miRNAs were identified using qRT-PCR, and their regulatory roles were confirmed by double-fluorescent reporter gene system and site-directed mutagenesis technology.

CONCLUSIONS

For the first time, allergen mRNA expression and miRNAs were profiled throughout the life cycle for an allergen-producing mite, and the results showed that miRNAs bind to target allergen mRNAs to regulate their expression.

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.

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

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

Analysis of the microRNA transcriptome of Daphnia pulex during aging

Daphnia pulex is an important food organism that exhibits a particular mode of reproduction known as cyclical parthenogenesis (asexual) and sexual reproduction. Regulation of the aging process by microRNAs (miRNAs) is a research hotspot in miRNA studies. To investigate a possible role of miRNAs in regulating aging and senescence, we used Illumina HiSeq to sequence two miRNA libraries from 1-day-old (1d) and 25-day-old (25d) D. pulex specimens. In total, we obtained 11,218,097 clean reads and 28,569 unique miRNAs from 1d specimens and 11,819,106 clean reads and 44,709 unique miRNAs from 25d specimens. Bioinformatic analyses was used to identify 1335 differentially expressed miRNAs from known miRNAs, including 127 miRNAs that exhibited statistically significant differences (P < 0.01); 92 miRNAs were upregulated and 35 were downregulated. Quantitative real-time (qRT)-PCR experiments were performed for nine miRNAs from five samples (1d, 5d, 10d, 15d, 20d and 25d) during the aging process, and the sequencing and qRT-PCR data were found to be consistent. Ninety-four miRNAs were predicted to correspond to 2014 target genes in known miRNAs with 4032 target gene sites. Sixteen pathways changed significantly (P < 0.05) at different developmental stages, revealing many important principles of the miRNA regulatory aging network of D. pulex. Overall, the difference in miRNA expression profile during aging of D. pulex forms a basis for further studies aimed at understanding the role of miRNAs in regulating aging, reproductive transformation, senescence, and longevity.

Characterization of transcriptome in the Indian meal moth Plodia interpunctella (Lepidoptera: Pyralidae) and gene expression analysis during developmental stages

The Indian meal moth Plodia interpunctella (Lepidoptera: Pyralidae) is a worldwide pest that causes serious damage to stored foods. Although many efforts have been conducted on this species due to its economic importance, the study of genetic basis of development, behavior and insecticide resistance has been greatly hampered due to lack of genomic information. In this study, we used high throughput sequencing platform to perform a de novo transcriptome assembly and tag-based digital gene expression profiling (DGE) analyses across four different developmental stages of P. interpunctella (egg, third-instar larvae, pupae and adult). We obtained approximate 9gigabyte (GB) of clean data and recovered 84,938 unigenes, including 37,602 clusters and 47,336 singletons. These unigenes were annotated using BLAST against the non-redundant protein databases and then functionally classified based on Gene Ontology (GO), Clusters of Orthologous Groups (COG), and Kyoto Encyclopedia of Genes and Genomes databases (KEGG). A large number of differentially expressed genes were identified by pairwise comparisons among different developmental stages. Gene expression profiles dramatically changed between developmental stage transitions. Some of these differentially expressed genes were related to digestion and cuticularization. Quantitative real-time PCR results of six randomly selected genes conformed the findings in the DGEs. Furthermore, we identified over 8000 microsatellite markers and 97,648 single nucleotide polymorphisms which will be useful for population genetics studies of P. interpunctella. This transcriptomic information provided insight into the developmental basis of P. interpunctella and will be helpful for establishing integrated management strategies and developing new targets of insecticides for this serious pest.

Pervasive microRNA Duplication in Chelicerates: Insights from the Embryonic microRNA Repertoire of the Spider Parasteatoda tepidariorum

MicroRNAs are small (∼22 nt) noncoding RNAs that repress translation and therefore regulate the production of proteins from specific target mRNAs. microRNAs have been found to function in diverse aspects of gene regulation within animal development and many other processes. Among invertebrates, both conserved and novel, lineage specific, microRNAs have been extensively studied predominantly in holometabolous insects such as Drosophila melanogaster However little is known about microRNA repertoires in other arthropod lineages such as the chelicerates. To understand the evolution of microRNAs in this poorly sampled subphylum, we characterized the microRNA repertoire expressed during embryogenesis of the common house spider Parasteatoda tepidariorum We identified a total of 148 microRNAs in P. tepidariorum representing 66 families. Approximately half of these microRNA families are conserved in other metazoans, while the remainder are specific to this spider. Of the 35 conserved microRNAs families 15 had at least two copies in the P. tepidariorum genome. A BLAST-based approach revealed a similar pattern of duplication in other spiders and a scorpion, but not among other chelicerates and arthropods, with the exception of a horseshoe crab. Among the duplicated microRNAs we found examples of lineage-specific tandem duplications, and the duplication of entire microRNA clusters in three spiders, a scorpion, and in a horseshoe crab. Furthermore, we found that paralogs of many P. tepidariorum microRNA families exhibit arm switching, which suggests that duplication was often followed by sub- or neofunctionalization. Our work shows that understanding the evolution of microRNAs in the chelicerates has great potential to provide insights into the process of microRNA duplication and divergence and the evolution of animal development.

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

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

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

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

Transcriptome Analysis of the Carmine Spider Mite, Tetranychus cinnabarinus (Boisduval, 1867) (Acari: Tetranychidae), and Its Response to β-Sitosterol

Tetranychus cinnabarinus (Acari: Tetranychidae) is a worldwide polyphagous agricultural pest that has the title of resistance champion among arthropods. We reported previously the identification of the acaricidal compound β-sitosterol from Mentha piperita and Inula japonica. However, the acaricidal mechanism of β-sitosterol is unclear. Due to the limited genetic research carried out, we de novo assembled the transcriptome of T. cinnabarinus using Illumina sequencing and conducted a differential expression analysis of control and β-sitosterol-treated mites. In total, we obtained >5.4 G high-quality bases for each sample with unprecedented sequencing depth and assembled them into 22,941 unigenes. We identified 617 xenobiotic metabolism-related genes involved in detoxification, binding, and transporting of xenobiotics. A highly expanded xenobiotic metabolic system was found in mites. T. cinnabarinus detoxification genes-including carboxyl/cholinesterase and ABC transporter class C-were upregulated after β-sitosterol treatment. Defense-related proteins, such as Toll-like receptor, legumain, and serine proteases, were also activated. Furthermore, other important genes-such as the chloride channel protein, cytochrome b, carboxypeptidase, peritrophic membrane chitin binding protein, and calphostin-may also play important roles in mites' response to β-sitosterol. Our results demonstrate that high-throughput-omics tool facilitates identification of xenobiotic metabolism-related genes and illustration of the acaricidal mechanisms of β-sitosterol.

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

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