Next Article in Journal
Trafficking and Function of the Voltage-Gated Sodium Channel β2 Subunit
Previous Article in Journal
Hinge Region in DNA Packaging Terminase pUL15 of Herpes Simplex Virus: A Potential Allosteric Target for Antiviral Drugs
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomolecules
Volume 9
Issue 10
10.3390/biom9100602
biomolecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Host Plant-Derived miRNAs Potentially Modulate the Development of a Cosmopolitan Insect Pest, Plutella xylostella

by
Ling-ling Zhang
1,2,
Xiao-dong Jing
1,2,
Wei Chen
1,2,
Yue Wang
1,2,
Jun-han Lin
1,2,
Ling Zheng
1,2,
Yu-hong Dong
1,2,
Li Zhou
1,2,
Fei-fei Li
1,2,
Fei-ying Yang
1,2,
Lu Peng
1,2,
Liette Vasseur
1,2,3,
Wei-yi He
1,2,* and
Min-sheng You
1,2,*
1
State Key Laboratory for Ecological Pest Control of Fujian and Taiwan Crops, Institute of Applied Ecology, Fujian Agriculture and Forestry University, Fuzhou 350002, China
2
International Joint Research Laboratory of Ecological Pest Control, Ministry of Education, Fujian Agriculture and Forestry University, Fuzhou 350002, China
3
Department of Biological Sciences, Brock University, St. Catharines, ON L2S3A1, Canada
*
Authors to whom correspondence should be addressed.
Biomolecules 2019, 9(10), 602; https://0-doi-org.brum.beds.ac.uk/10.3390/biom9100602
Submission received: 20 August 2019 / Revised: 26 September 2019 / Accepted: 10 October 2019 / Published: 12 October 2019
(This article belongs to the Section Molecular Biology)
Browse Figures
Versions Notes

Abstract

:
Plant microRNAs (miRNAs) have recently been reported to be involved in the cross-kingdom regulation of specific cellular and physiological processes in animals. However, little of this phenomenon is known for the communication between host plant and insect herbivore. In this study, the plant-derived miRNAs in the hemolymph of a cruciferous specialist Plutella xylostella were identified by small RNAs sequencing. A total of 39 miRNAs with typical characteristics of plant miRNAs were detected, of which 24 had read counts ≥ 2 in each library. Three plant-derived miRNAs with the highest read counts were validated, and all of them were predicted to target the hemocyanin domains-containing genes of P. xylostella. The luciferase assays in the Drosophila S2 cell demonstrated that miR159a and novel-7703-5p could target BJHSP1 and PPO2 respectively, possibly in an incomplete complementary pairing mode. We further found that treatment with agomir-7703-5p significantly influenced the pupal development and egg-hatching rate when reared on the artificial diet. The developments of both pupae and adults were severely affected upon their transfer to Arabidopsis thaliana, but this might be independent of the cross-kingdom regulation of the three plant-derived miRNAs on their target genes in P. xylostella, based on expression analysis. Taken together, our work reveals that the plant-derived miRNAs could break the barrier of the insect mid-gut to enter the circulatory system, and potentially regulate the development of P. xylostella. Our findings provide new insights into the co-evolution of insect herbivore and host plant, and novel direction for pest control using plant-derived miRNAs.

1. Introduction

MicroRNAs (miRNAs) are small non-coding RNAs of 19–24 nucleotides [1]. In both animals and plants, miRNAs play an important role in post-transcriptional regulation by interacting with target mRNAs involved in cell differentiation, apoptosis, metabolism and response to environmental stress [2]. However, the miRNAs in the two kingdoms are quite different with the plant precursor miRNA molecules being longer and more variable than those of animals [3]. Proteins participating in the biogenesis of miRNAs also differ between plants and animals [4,5] with 2′-O-methylation at 3′ end by the methyltransferase Hua enhancer 1 (HEN1) being involved in plants [6]. Animal miRNAs regulate gene expression through the 3′-UTR regions of their target mRNAs to repress translation, while most plant miRNAs mediate cleavage of target mRNAs in the coding regions [7].
It has been reported that plant miRNAs are able to pass through the strict gastrointestinal barrier of animals when ingested, circulating in the blood, and eventually exerting biological functions by regulating gene expression in various tissues or cells [8,9,10,11]. MiRNA-mediated cross-kingdom communication has been first reported in mammalian liver, where rice-derived miR168a alters the physiological status of mice via its regulating of the low-density lipoprotein receptor adapter protein 1 (LDLRAP1) [12]. Plant-derived miRNAs have been closely associated with cancer development [13,14] and viral infection [15,16]. These findings suggest that the miRNAs of plants when ingested may contribute to the physiological changes of animals [17].
In insects, plant-derived miRNAs have been found to accumulate and move systemically in lepidopteran (Helicoverpa zea and Spodoptera frugiperda) and coleopteran (Diabrotica virgifera virgifera) insects under the controlled feeding experiments [18]. According to small RNA (sRNA) sequencing and digital droplet polymerase chain reaction (PCR), several mulberry miRNAs can enter different tissues of silkworm with no effect on potential target genes and phenotype [19]. In honeybees, pollen-derived miR162a enriched in beebread targets mTOR prevents the differentiation of larvae into queens, which develop instead as workers [20]. The miR162a-mTOR network involved in caste dimorphism is also conservative in the non-social insect, Drosophila melanogaster (common fruit fly) [20]. In two cereal aphids, Schizaphis graminum and Sipha flava, 13 sorghum miRNAs (eight known and five novel) and three novel barley miRNAs have been detected and are predicted to target aphid genes involved in detoxification, and starch and sucrose metabolism [21].
However, the concept has been controversial since the first identification of the cross-kingdom regulation of plant-derived miRNAs [12]. Other findings have shown that inadequate or little general uptake of dietary miRNAs can be found in different tissues of adult bees [22] and in the plasma of mammalian blood cells, as an example, with routine plant ingestion [23,24]. Further, two independent studies through the bioinformatic approach and diet feeding experiment revealed the non-universality of dietary uptake and an absence of cross-kingdom gene regulation by plant miRNAs [18,25]. Further research is therefore essential to investigate the possibility of such a concept in other species.
Currently, miRNAs studies related to Plutella xylostella, a worldwide destructive pest of cruciferous vegetables [26,27,28], are mainly targeting insecticide resistance [29,30,31] and immune response [32,33]. As an oligophagous insect herbivore, P. xylostella may have long been evolving with miRNAs derived from host plants. Therefore, in this study, we explored whether plant-derived miRNAs could be detected in the circulatory system and exert functions in P. xylostella using high-throughput sequencing, molecular, cellular and biological assays.

2. Materials and Methods

2.1. Insects and Sampling

The insecticide-susceptible G88 strain of P. xylostella was maintained under laboratory conditions and on a commercially-available artificial diet (#F9772-DBM, Frontier Scientific Services, Newark, DE, USA) for three years without exposure to any host plant. The GC strain of P. xylostella was established by transferring the newly emerged 1st instar G88 onto Col-0 type A. thaliana until 4th instar for sampling. Insects were reared at 25 ± 1 °C, an RH of 60 ± 5% and a photoperiod of 16:8 h (L:D).

2.2. Plants

Arabidopsis thaliana (Col-0 ecotype) was used in this study. Seeds were surface-sterilized for 5 min with 75% ethanol and washed five times using sterile water. Seeds were suspended in sterile water and then vernalized at 4 °C in the dark for three days before germination.
Seedlings were grown at 25 °C, under a 16:8 h photoperiod, on vermiculite, and were provided with a commercial nutrient solution (Coolaber, Beijing, China) every day. The intensity of light was set to 2,000 lux and RH at 65 ± 5%.

2.3. Small RNA Sequencing

To build two small RNA (sRNA) libraries of the GC strain, for each library, the hemolymph of 30 first-day 4th instar larvae were dissected and pooled. Hemolymph samples were collected in the lysis buffer of a Quick-RNA Tissue/Insect Microprep kit (#R2030, ZymoResearch, Irvine, CA, USA) using a glass needle, and centrifuged at 13,000× g for 1 min at room temperature for RNA extraction, using the same kit according to the manufacturer’s instructions. Sequencing libraries were constructed using the NEBNext® UltraTM small RNA Sample Library Prep Kit for Illumina, and the sequencing procedure was performed by Biomarker Technologies Corporation (Illumina HiSeq2500, Beijing, China). Clean reads were compared to the known plant microRNA (miRNA) precursors deposited in the online miRBase (version 21, http://www.mirbase.org/). Unknown plant miRNAs were inferred by miRdeep2 software using the settings of “-g -1 -l 250 -b 0” based on the information of precursors in the genome location and their structural energies [34]. The identified miRNAs that were supposed to be plant-derived were further used to map the P. xylostella genome using Bowtie [35] with a default parameter to exclude the possibility of the presence of these miRNAs.

2.4. RNA Extraction and qRT-PCR

For miRNA profiling, total RNA was extracted using the Quick-RNA Tissue/Insect Microprep kit (#R2030, ZymoResearch, Irvine, CA, USA) according to the manufacturer’s instructions. The quality and concentrations of total RNA were estimated by electrophoresis and a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Shanghai, China). The first-strand cDNA of miRNA was synthesized using the stem-loop strategy and a GoScriptTM Reverse Transcription System kit (#TM316, Promega, Madison, WI, USA). The miRNA-specific reverse primer was synthesized by adding the reverse complementary sequence of the last 8 bp of mature miRNA to the 3′ terminal of the common stem-loop structure (Table S1). The remainder of the sequence of mature miRNA was added to a universal adapter at the 5′ end as the forward primer paired with a universal reverse primer. The relative levels of miRNAs were normalized using the U6 snRNA gene [36].
For the RNA extraction of each of the three biological replicates, six larvae of the 3rd instar or four larvae of the 4th instar were collected and pooled. The pooled samples were frozen in liquid nitrogen immediately and stored at −80 °C until use. For mRNA profiling, total RNA was isolated using TRIzol reagent (#15596018, Invitrogen, Waltham, MA, USA). according to the manufacturer’s instruction. The quality and concentrations of total RNAs were estimated by electrophoresis and the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Shanghai, China). The first-strand cDNA was synthesized using a GoScriptTM Reverse Transcription System kit (#TM316, Promega, Madison, WI, USA). The gene specific primer pairs were presented in Table S1. EF1 was used as the internal control gene for normalizing the mRNA expression levels [37].
Quantitative real-time polymerase chain reaction (qRT-PCR) assays for miRNA and mRNA were conducted on the Real time polymerase chain reaction (PCR) system (CFX96, Bio-rad, Hercules, CA, USA). For each biological replicate, the qRT-PCR tests were conducted using three technical replicates. The reaction condition was as follows: 95 °C for 3 min, 40 cycles at 95 °C for 30 s and 58 °C for 30 s. At the end of the cycling, melting curves of the resulting PCR products were obtained by gradually increasing the temperature from 55 °C to 95 °C with an increment of 0.5 °C /5 s to check the specificity of the primers. The qRT-PCR data were analyzed using methods of 2−∆∆Ct [38] for comparing the miRNA/mRNA expression levels between the control and treatments, or 2−∆Ct [39] for comparing expression levels between two strains.

2.5. Target Prediction and Transcriptome Sequencing

Based on the CDSs of the P. xylostella genome and a UTR dataset for 3,050 genes obtained from DBM-DB [36], the potential target binding sites of three miRNAs were predicted by miRanda [40] and RNAhybrid [41] using the built-in default parameter and the parameter of Drosophila melanogaster, respectively, with a setting of minimal free energy (mfe) ≤ −15 kcal/mol.
For transcriptome sequencing, total RNA was isolated from G88 larvae at 24 h post-injection of each of three miRNA agomirs. Larvae treated with agomir-NC were used as the control. Six larvae were pooled for each replicate, and two replicates were prepared for each treatment. A total of 8 RNA samples were sequenced and analyzed at BGI Genomics Company (BGISeq-500, Shenzhen, China). Genome information of P. xylostella were downloaded from DBM-DB [42]. Differentially-expressed genes between the miRNA agomir treatment and the control were annotated and classified with GO terms and the KEGG pathway, of which down-regulated genes were used for generating intersection genes with in silico results. The enriched GO term/KEGG pathway was defined as a set of genes, which were annotated to a specific GO term/KEGG pathway and overrepresented in the differentially-expressed genes when compared with the number of genes in the genome annotated to this specific GO term/KEGG pathway [43].

2.6. Cell Line and Luciferase Assays

The Drosophila S2 cell line was maintained at 24 °C in Grace medium plus 10% fetal bovine serum (FBS) (#16140, Gibco, Waltham, MA, USA). The potential miRNA target sites (Table S1) were added with the restriction enzyme sites of SacI and SalII, which were complementary to multiple cloning sites of pmirGLO vector (#E1330, Promega, Madison, WI, USA), and ligated into the vector in the original orientation (wild-type). Mutations were introduced at the binding region of the miRNA target sites, which were used as the negative control (mutant). The S2 cells plated in 96-well plates (#3988, Costar, New York, NY, USA) were co-transfected with 100 ng of each plasmid, 100 nM agomir or agomir-NC (negative control with no similarity to the miRNAs in the miRBase), and 2 µL CellfectinII reagent (#10362100, Invitrogen, Waltham, MA, USA). The S2 cells were analyzed for luciferase activity using the Dual-Glo® Luciferase Assay System and GloMax® 20/20 Luminometer (#E2920, Promega, Madison, WI, USA) at 48 h post-transfection. The relative luciferase activity (firefly luciferase/renilla luciferase) of the group containing agomir-NC and the wild-type pmirGLO vector was used as our control. For each transfection, the relative luciferase activity was averaged from five biological replicates.

2.7. miRNA Agomir Treatments

The agomirs (double-strand mature miRNAs) of miR-159a, miR-166a-3p and novel-7703-5p were synthesized by Sangon Biotech (Shanghai, China) (Table S1). The anti-sense strand of agomir was modified with cholesterol and four thioskeletons at the 3′-terminal, four thioskeletons at the 5′-terminal, and 2′-methoxy groups for the entire strand. First-day 3rd instar larvae of the G88 strain were selected for injection (69 nL of each miRNA agomir at 100 µM per larva) using nanoliter injection 2000 (World Precision Instruments, USA). Larvae injected with agomir-NC were used as the control group. For phenotypic investigation, three biological replicates were used, and each replicate contained 20 larvae treated with miRNA agomir or agomir-NC.

2.8. Western Blots

Total proteins were extracted using a phenol/SDS extraction method [44] and quantified by a BCA Protein Quantification Kit (Yeasen, Shanghai, China). The heat-denatured proteins were subjected to 10% SDS-PAGE electrophoresis and then transferred to a PVDF membrane (Millipore, Burlington, MA, USA). After blocking in BSA (1:1000, Yeasen, Shanghai, China) at room temperature for 1 h, the membrane was incubated in anti-PPO2 polyclonal antibody (1:1000) or anti-tubulin antibody (1:3,000) at 4 °C overnight. After that, one of the secondary antibodies was added, and the reactant was incubated at room temperature for 1 h.
The immunological blots were detected using ECL Western Blot Kit (#1705060S, Bio-rad, Hercules, CA, USA) and photographed using a Fusion FX imaging system (version 7.0, Vilber, Marne-la-Vallée, France). The density of the bands was analyzed using ImageJ (version 1.52) [45].

2.9. Statistical Analysis

Data were presented as mean ± SD and analyzed by IBM SPSS Statistics version 21 software (Columbia, Carolina, USA). Student’s t test or one-way Analysis of Variance (ANOVA) followed by Tukey’s multiple comparisons was used to determine the significant difference at p-value < 0.05.

2.10. Data Availability

The raw sequence data of sRNA-seq and RNA-seq have been deposited in the Genome Sequence Archive (GSA) [46] in the BIG Data Center [47] operated by the Beijing Institute of Genomics (BIG), Chinese Academy of Sciences, under accession numbers CRA001879 and CRA001896, respectively, which are publicly accessible at http://bigd.big.ac.cn/gsa.

3. Results

3.1. Characterization of Plant-Derived miRNAs in Larval Hemolymph of P. xylostella

The two sRNA libraries of the larval hemolymph of the GC strain shared 89.57% common reads (Figure 1A). After filtering rRNAs, tRNAs, scRNAs, snRNAs, snoRNAs and repeats-associated reads, the residual unannotated reads that included miRNAs accounted for more than 50% (Figure 1B,C). A total of 39 miRNAs, including 33 previously-reported and 6 novel sequences from the larval hemolymph of P. xylostella were well aligned to the genome of Arabidopsis thaliana (Table 1). After filtering those miRNAs with read counts < 2 in each library, the final miRNA set contained 21 known and three novel sequences. The Bowtie-based mapping result showed that none of the 39 miRNAs could be mapped to our P. xylostella genome, indicating that they are all plant-derived. The length distribution of these miRNAs was dominant at 21 nt (Figure 1D) and the first base at 5′-end had a bias for U (Figure 1E), which were consistent with the specific recognition and cleavage mode of Dicer-like enzymes (DCL1) on plant miRNAs precursor.
To confirm that these plant-derived miRNAs were not from contamination during sRNA-seq or sample preparation, stem-loop PCR was used to clone three plant miRNAs with top high read counts (>100), including miR159a, miR166a-3p and novel-7703-5p (Table 1). All of these miRNAs were validated using Sanger sequencing and were consistent with sequences from sRNA-seq. The relative levels of three miRNAs in the larval hemolymph of the GC strain were significantly higher than in the G88 strain, although all of them were also detected in the G88 strain (Figure 1F), which may be due to the plant-derived additive in the artificial diet (#F9772-DBM, Frontier Scientific Services, USA), including linseed oil, soy flour and wheat germ.
By searching mature sequences of miR159a, miR166a-3p and novel-7703-5p with that of other species in the miRBase (version 22) [48], we found that both miR159a and miR166a-3p widely existed in monocotyledons, Brassicaceae, Fabaceae, and Rosaceae.
The mature sequence of miR159a varied at both ends, while miR166a-3p was highly conserved (Figure S1A,B). No existing sequences of other plants in the miRBase were similar to novel-7703-5p, except for A. thaliana. To further test the conservative nature of miR159a, miR166a-3p and novel-7703-5p in the host plant, we used stem loop PCR to clone these miRNAs in cabbage (Brassica oleracea), oilseed rape (Brassica napus), and radish (Raphanus sativus). Their mature sequences were identical in the tested host plants (Figure S1C–E), which might imply a universal connection between these miRNAs and P. xylostella.

3.2. Target Prediction of Plant-Derived miRNAs in P. xylostella

Based on the predicted CDSs of the P. xylostella genome [49] and a UTR dataset for 3,050 genes generated by a previously reported transcriptome data [43], two types of miRNA target prediction software, miRanda [40] and RNAHybrid [41], were used to predict the target genes and their corresponding binding sites using the built-in default parameter and the parameter of Drosophila melanogaster, respectively. The predicted target genes from both software were integrated for each of the miRNAs for further analysis (Tables S2–S4).
To further improve the accuracy of in silico prediction, we used RNA-seq to screen the differentially expressed genes (DEGs) in the 3rd instar larvae of the G88 strain in response to the injection of each miRNA agomir. The relative levels of the three miRNAs showed a peak at 24 h post-injection and decreased rapidly by 48 h (Figure S2). Therefore, RNA-seq was performed to capture the global change of gene expression in the larvae of P. xylostella treated with a miRNA agomir at 24 h. Compared with the negative control, the total numbers of DEGs were 1,969, 1,242 and 2,280 for treatments of agomir-159a, agomir-166a-3p and agomir-7703-5p, respectively. Most of the DEGs were up-regulated with few of them down-regulated (Figures S3–S5). Annotation through GO terms and KEGG pathways for gene enrichment analysis of the DEGs between each treatment and the control suggested that miR159a, miR166a-3p or novel-7703-5p might potentially influence cellular and metabolism processes in P. xylostella (Figures S3–S5).
Given that miRNAs are known to mostly mediate the down-regulation of the target genes, the down-regulated DEGs in each agomir treatment were considered to be the corresponding target genes. Based on the intersection of in silico prediction and down-regulated genes from RNA-seq data, we identified 12, 1 and 8 target genes for miR159a, miR166a-3p and novel-7703-5p, respectively (Table 2). Three of these target genes belonged to the arthropod hemocyanin superfamily, members of which participate in multiple important functions in the hemolymph of arthropods [50]. They were predicted to encode basic juvenile hormone-suppressible protein 1 (BJHSP1, Px006820), basic juvenile hormone-suppressible protein 2 (BJHSP2, Px007031) and phenoloxidase subunit 2 (PPO2, Px002274) targeted by miR159a, miR166a-3p and novel-7703-5p, respectively (Figure 2A–C). All three conserved hemocyanin domains (C, M and N) were contained but differently organized within each of the three target genes (Figure S6).

3.3. Suppression and Binding of Plant-Derived miRNAs to Target Genes in P. xylostella

To understand whether the expression of target genes could be regulated by three plant-derived miRNAs in vivo, we injected the 3rd instar larvae with each miRNA agomir and measured the expression using qRT-PCR. The results showed that the expression levels of the target genes were significantly suppressed when compared with the control (Figure 2D–F). We further examined the protein levels of PPO2 in response to the agomir-7703-5p treatment using Western blot and found that PPO2 levels were decreased at 24 h post-injection and soon recovered by 48 h (Figure 3), showing a same pattern as the gene expression.
The potential binding sites for miR159a, miR166a-3p and novel-7703-5p, in an incomplete complementary pairing mode, were all located in the open reading frame of their target genes (Figure 2A–C). To test whether miR159a, miR166a-3p and novel-7703-5p could bind to their target genes in vitro, the binding sites were inserted into downstream of the luciferase gene in the pmirGLO vector and transfected into S2 cells. The normalized luciferase activity was significantly increased when binding sites of miR159a and novel-7703-5p were mutated, while there was no obvious change using the vector containing the mutated binding site of miR166a-3p (Figure 4). These results indicated that miR159a and novel-7703-5p had a direct relationship with BJHSP1 and PPO2, respectively. Further, we noticed that there might be an unexpected interaction between agomir-NC and the potential binding site of novel-7703-5p on PPO2. However, the degree of binding by agomir-NC on PPO2 might be less than that of novel-7703-5p (Figure 4C).

3.4. Effects of Plant-Derived miRNAs on the Development of P. xylostella

To further explore the effects of plant-derived miRNAs on the development of P. xylostella, biological parameters reflecting the fitness were compared between the miRNA agomir treatments and agomir-NC control. We found that there were no significant differences in the survival rate of larvae between the control and treatments at 24 h and 48 h post-injection, as well as for the weight of the pupae, rate of pupation, number of laid eggs and the hatching rate (Figure S7). However, treatments of agomir-166a-3p and agomir-7703-5p significantly shortened the duration of the larval stage compared with the control (Figure 5A). Further, the agomir-7703-5p treatment also significantly resulted in the prolongation of the pupal stage, a higher percentage of abnormalities of the pupae and a decrease in the adult emergence rate (Figure 5B–D). Abnormal pupae died soon after emergence, with their body turning black and wizened (Figure 5E).

3.5. Link of Decreased Fitness of the Larvae on Host Plant with Plant-Derived miRNAs

After the transfer of newly-emerged 1st instar larvae of the G88 strain onto A. thaliana, there was little impact on the duration and the survival rate of larvae, pupal duration and adult emergence rate (Figure S8). However, the pupal weight, pupation rate, number of laid eggs and hatching rate were significantly reduced when compared with the G88 individuals feeding on artificial diet (Figure 6), indicating influence on the quality of pupae and adults rather than the larvae of P. xylostella.
To test whether plant-derived miRNAs were involved in a plant defense to an insect herbivore, we examined the relative levels of the three miRNAs in intact leaves of A. thaliana and the leaves fed by the 3rd instar larvae of the G88 strain for 8 h, based on when most defense-related transcripts are up-regulated [51]. The results showed that the levels of miR159a, miR166a-3p and novel-7703-5p in leaves damaged by P. xylostella were significantly reduced as compared with that of fresh intact leaves (Figure 7A), which might exclude the possibility of any natural occurrence of defensive plant-derived miRNAs.
We also compared the expression levels of the target genes BJHSP1, BJHSP2 and PPO2 in the 3rd and 4th instar larvae and newly-emerged pupae between the G88 and GC strains. In most cases, the relative expression levels of the three target genes in the GC strain were significantly higher than that of the G88 strain, except for BJHSP1, which was expressed at a higher level in the 4th instar larvae of the G88 strain (Figure 7B–D). This seemed to be consistent with our finding that the three plant-derived miRNAs were unlikely to be a key factor in conferring the resistance of the host plant to an insect herbivore.

4. Discussion

In this study, plant-derived miRNAs in P. xylostella were comprehensively investigated in the larval hemolymph of the 4th instar larvae. We found that a plant-derived miRNA novel-7703-5p could down-regulate both at the transcriptional and protein levels of PPO2 and was connected to the death of pupae. Although the cross-kingdom gene regulation of dietary miRNAs has long been controversial, our study is among the first to provide evidence of such a phenomenon in a host plant-insect pest system. Our results support the idea that miRNAs are interrelated mediators of disparate species, which can be transferred from a low trophic level through the food chain and exhibit biological activities at a high trophic level [8,9,10,11,12].
It is assumed that the abundance of plant miRNAs detected in animal tissues should correlate with the miRNA abundance in plant tissues unless a particular miRNA is highly stable or given priority to be absorbed in animals [12]. For example, the average read count of mulberry miR166b in silkworm is low, and is nonfunctional in cross-kingdom gene regulation [19]. However, the read count of miR162a in the food of honeybee larvae is relatively high (362), and shows a significant effect on cross-kingdom gene regulation [20]. Based on this notion and the sRNA-seq data, we found that the average read counts of the top three plant-derived miRNAs in P. xylostella, miR159a (258), miR166a-3p (135) and novel-7703-5p (100), were also relatively high. We therefore hypothesized that these plant-derived miRNAs might have reached the concentration required for their absorption in P. xylostella. However, the relative levels of these miRNAs were significantly reduced in the larva-infested leaves compared with that of fresh leaves, indicating that the feeding of insect herbivores was able to influence miRNAs expression of host plants.
The expression levels of targets genes of the three plant-derived miRNAs in the GC strain after ingesting plant tissues were higher than that of the G88 strain in most developmental stages, although the overall performance of the GC strain was poorer than the G88 strain. This contradicted our findings that the direct injection of high doses of miRNA agomirs led to an inhibition of target genes and significantly affected the development and reproduction of P. xylostella. These results suggest that the normal number of miRNAs ingested from plants may not always be fully utilized and be effectively functional in herbivores [52]. Therefore, further studies on the factors contributing to the induction of miRNAs in natural host plants, which have a cross-kingdom regulatory potential with insect genes [53,54], may provide a novel direction to develop pest control free from genetic transformation.
Both in vivo and in vitro experiments showed that miR159a and novel-7703-5p could alter the transcriptional levels of BJHSP1 and PPO2 and even affect the physiological processes of P. xylostella. It is suggested that these exogenous miRNAs are able to resist a harsh digestive tract environment, penetrate the gut cells, and be transported to target tissues or cells via the hemolymph [9,10]. The high selectivity of insect gut epithelial cells is one of the most stringent biological barriers, and both normal and pathological changes of the barrier may limit or affect the uptake of small RNA molecules [55,56,57]. The efficiency of short interfering RNAs (siRNAs) absorption and RNA interference (RNAi) greatly varies among insect species [58]. Lepidopteran insects often require a high level of siRNAs to induce RNAi, but much less siRNAs may be provided in the food sources [59]. Some studies believe that miRNAs freely transported across species may be related to the effective utilization of host RNAi pathways [60]. To safely transport and protect from degradation, miRNAs loaded within RISC complexes are usually encapsulated in microvesicles, exosomes and apoptosomes [61]. Along with the circulatory system, miRNAs are transported to the target cells and fully used by the RNAi machinery inside the cells [62]. Plant miRNAs are extraordinarily stable under extreme conditions, such as the presence of RNase, high temperature, overextended storage time, repeated freezing and thawing, since they are processed with 2′-O-methylation on their 3′end and protected by binding with membrane proteins [12,13,15,63]. However, the specific mechanism of how plant miRNAs make their ways to the target cells of animals and become active molecules remains to be further investigated.
The homologous mature sequences of miR159a, miR166a-3p and novel-7703-5p have been found to be present and conservative in diverse cruciferous plants, which may indicate that the crucifer explores these miRNAs to interact with P. xylostella. These three miRNAs shared no similarity in sequences but had a similar potential to target the genes containing hemocyanin domains, which play important roles in the hemolymph of arthropods [50]. These results differed from the findings of other insect herbivores [19,20,21], making hemolymph a hot spot for cross-kingdom gene regulation. Considering that the plant-derived miRNAs in insect hemolymph possibly use the RISC complex of host for regulating host gene expression, we assume that plant-derived miRNAs may function in insect hemolymph through a similar mode to animal miRNAs [7]. However, the predicted results combined with the RNA-seq of miRNA agomir treatments and luciferase assays showed that the binding sites of three plant-derived miRNAs were all located in the coding regions of their target genes, implying that plant-derived miRNAs might also exert cross-kingdom gene regulation through cleavage on the target mRNAs in a manner commonly used by plants [4].
Overall, the natural uptake and utilization of small RNA molecules in plants by herbivorous insects can become a novel avenue to understand the complex plant-insect interactions and develop new technologies for pest management based on miRNA.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2218-273X/9/10/602/s1, Figure S1: title, Table S1: title, Video S1: title. Supplementary Materials. Figure S1. Conservation analysis of plant-derived miRNAs. Figure S2. Expression of three miRNAs at different time points after agomir treatments. Figure S3. Enrichment analysis of differentially expressed genes in agomir-159a treatment. Figure S4. Enrichment analysis of differentially expressed genes in agomir-166a-3p treatment. Figure S5. Enrichment analysis of differentially expressed genes in agomir-7703-5p treatment. Figure S6. Predicted conserved domains of BJHSP1, BJHSP2 and PPO2. Figure S7. The biological parameters of the G88 strain that were not affected by treatments with miRNA agomirs. Figure S8. The biological parameters showing no difference between the G88 and GC strains. Table S1. Sequences of primer pairs, miRNA agomirs and target binding sites. Table S2. Target gene prediction for miR159a. Table S3. Target gene prediction for miR166a-3p. Table S4. Target gene prediction for novel-7703-5p.

Author Contributions

W.-y.H., L.P. and L.-l.Z. designed research; L.-l.Z., X.-d.J., W.C., J.-h.L., L.Z. (Ling Zheng)., Y.-h.D., L.Z. (Li Zhou), F.-f.L. and F.-y.Y. performed experiment; L.-l.Z., Y.W. and W.-y.H. analyzed data; L.-l.Z., W.-y.H., M.-s.Y. and L.V. wrote paper.

Funding

This work was supported by National Key R&D Program of China (2017YFD0200400), Natural Science Foundation of Fujian Province (2019J01369), Fujian Provincial Science and Technology Major Project (2018NZ0002-1) and Innovation Fund of FAFU (CXZX2017471, CXZX2016128, CXZX2017321 and CXZX2018092).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef]
  2. Carrington, J.C.; Ambros, V. Role of microRNAs in plant and animal development. Science 2003, 301, 336–338. [Google Scholar] [CrossRef] [PubMed]
  3. Axtell, M.J.; Westholm, J.O.; Lai, E.C. Vive la différence: Biogenesis and evolution of microRNAs in plants and animals. Genome Biol. 2011, 12, 221. [Google Scholar] [CrossRef] [PubMed]
  4. Voinnet, O. Origin, biogenesis, and activity of plant microRNAs. Cell 2009, 136, 669–687. [Google Scholar] [CrossRef] [PubMed]
  5. Moran, Y.; Agron, M.; Praher, D.; Technau, U. The evolutionary origin of plant and animal microRNAs. Nat. Ecol. Evol. 2017, 1, 27. [Google Scholar] [CrossRef]
  6. Yu, B.; Yang, Z.Y.; Li, J.J.; Minakhina, S.; Yang, M.C.; Padgett, R.W.; Steward, R.; Chen, X.M. Methylation as a crucial step in plant microRNA biogenesis. Science 2005, 307, 932–935. [Google Scholar] [CrossRef]
  7. Stark, A.; Brennecke, J.; Bushati, N.; Russell, R.B.; Cohen, S.M. Animal microRNAs confer robustness to gene expression and have a significant impact on 3 ‘ UTR evolution. Cell 2005, 123, 1133–1146. [Google Scholar] [CrossRef]
  8. Jiang, M.X.; Sang, X.L.; Hong, Z. Beyond nutrients: Food-derived microRNAs provide cross-kingdom regulation. Bioessays 2012, 34, 280–284. [Google Scholar] [CrossRef] [Green Version]
  9. Han, L.; Luan, Y.S. Horizontal transfer of small RNAs to and from plants. Front Plant Sci. 2015, 6, 1113. [Google Scholar] [CrossRef]
  10. Zhang, H.; Li, Y.P.; Liu, Y.N.; Liu, H.M.; Wang, H.Y.; Jin, W.; Zhang, Y.M.; Zhang, C.; Xu, D. Role of plant microRNA in cross-species regulatory networks of humans. BMC Syst. Biol. 2016, 10, 60. [Google Scholar] [CrossRef]
  11. Zeng, J.; Gupta, V.K.; Jiang, Y.; Yang, B.; Gong, L.; Zhu, H. Cross-kingdom small RNAs among animals, plants and microbes. Cells 2019, 8, 371. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, L.; Hou, D.; Chen, X.; Li, D.; Zhu, L.; Zhang, Y.; Li, J.; Bian, Z.; Liang, X.; Cai, X.; et al. Exogenous plant MIR168a specifically targets mammalian LDLRAP1: Evidence of cross-kingdom regulation by microRNA. Cell Res. 2012, 22, 107–126. [Google Scholar] [CrossRef] [PubMed]
  13. Chin, A.R.; Fong, M.Y.; Somlo, G.; Wu, J.; Swiderski, P.; Wu, X.; Wang, S.E. Cross-kingdom inhibition of breast cancer growth by plant miR159. Cell Res. 2016, 26, 217–228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Kumar, D.; Kumar, S.; Ayachit, G.; Bhairappanavar, S.B.; Ansari, A.; Sharma, P.; Soni, S.; Das, J. Cross-kingdom regulation of putative miRNAs derived from happy tree in cancer pathway: A systems biology approach. Int. J. Mol. Sci. 2017, 18, 1191. [Google Scholar] [CrossRef] [PubMed]
  15. Zhou, Z.; Li, X.; Liu, J.; Dong, L.; Chen, Q.; Liu, J.; Kong, H.; Zhang, Q.; Qi, X.; Hou, D.; et al. Honeysuckle-encoded atypical microRNA2911 directly targets influenza A viruses. Cell Res. 2015, 25, 39–49. [Google Scholar] [CrossRef] [PubMed]
  16. Li, X.; Huang, Y.; Sun, M.; Ji, H.; Dou, H.; Hu, J.; Yan, Y.; Wang, X.; Chen, L. Honeysuckle-encoded microRNA2911 inhibits Enterovirus 71 replication via targeting VP1 gene. Antiviral Res. 2018, 152, 117–123. [Google Scholar] [CrossRef] [PubMed]
  17. Sala-Cirtog, M.; Marian, C.; Anghel, A. New insights of medicinal plant therapeutic activity-the miRNA transfer. Biomed. Pharmacother. 2015, 74, 228–232. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, Y.; Wiggins, B.E.; Lawrence, C.; Petrick, J.; Ivashuta, S.; Heck, G. Analysis of plant-derived miRNAs in animal small RNA datasets. BMC Genomics 2012, 13, 381. [Google Scholar] [CrossRef] [PubMed]
  19. Jia, L.; Zhang, D.; Xiang, Z.; He, N. Nonfunctional ingestion of plant miRNAs in silkworm revealed by digital droplet PCR and transcriptome analysis. Sci. Rep. 2015, 5, 12290. [Google Scholar] [CrossRef] [Green Version]
  20. Zhu, K.; Liu, M.; Fu, Z.; Zhou, Z.; Kong, Y.; Liang, H.; Lin, Z.; Luo, J.; Zheng, H.; Wan, P.; et al. Plant microRNAs in larval food regulate honeybee caste development. PLoS Genet. 2017, 13, e1006946. [Google Scholar] [CrossRef]
  21. Wang, H.; Zhang, C.; Dou, Y.; Yu, B.; Liu, Y.; Heng-Moss, T.M.; Lu, G.; Wachholtz, M.; Bradshaw, J.D.; Twigg, P.; et al. Insect and plant-derived miRNAs in greenbug (Schizaphis graminum) and yellow sugarcane aphid (Sipha flava) revealed by deep sequencing. Gene 2017, 599, 68–77. [Google Scholar] [CrossRef] [PubMed]
  22. Masood, M.; Everett, C.P.; Chan, S.Y.; Snow, J.W. Negligible uptake and transfer of diet-derived pollen microRNAs in adult honey bees. RNA Biol. 2016, 13, 109–118. [Google Scholar] [CrossRef] [PubMed]
  23. Witwer, K.W.; McAlexander, M.A.; Queen, S.E.; Adams, R.J. Real-time quantitative PCR and droplet digital PCR for plant miRNAs in mammalian blood provide little evidence for general uptake of dietary miRNAs: Limited evidence for general uptake of dietary plant xenomiRs. RNA Biol. 2013, 10, 1080–1086. [Google Scholar] [CrossRef] [PubMed]
  24. Mico, V.; Martin, R.; Lasuncion, M.A.; Ordovas, J.M.; Daimiel, L. Unsuccessful detection of plant microRNAs in beer, extra virgin olive oil and human plasma after an acute ingestion of extra virgin olive oil. Plant Foods Hum. Nutr. 2016, 71, 102–108. [Google Scholar] [CrossRef] [PubMed]
  25. Dickinson, B.; Zhang, Y.; Petrick, J.S.; Heck, G.; Ivashuta, S.; Marshall, W.S. Lack of detectable oral bioavailability of plant microRNAs after feeding in mice. Nat. Biotechnol. 2013, 31, 965. [Google Scholar] [CrossRef] [PubMed]
  26. Talekar, N.S.; Shelton, A.M. Biology, ecology, and management of the diamondback moth. Annu. Rev. Entomol. 1993, 38, 275–301. [Google Scholar] [CrossRef]
  27. Furlong, M.J.; Wright, D.J.; Dosdall, L.M. Diamondback moth ecology and management: Problems, progress, and prospects. Annu. Rev. Entomol. 2013, 58, 517–541. [Google Scholar] [CrossRef] [PubMed]
  28. Li, Z.; Feng, X.; Liu, S.S.; You, M.; Furlong, M.J. Biology, ecology, and management of the diamondback moth in China. Annu. Rev. Entomol. 2016, 61, 277–296. [Google Scholar] [CrossRef] [PubMed]
  29. Li, X.X.; Guo, L.; Zhou, X.G.; Gao, X.W.; Liang, P. miRNAs regulated overexpression of ryanodine receptor is involved in chlorantraniliprole resistance in Plutella xylostella (L.). Sci. Rep. 2015, 5, 14095. [Google Scholar] [CrossRef]
  30. Zhu, B.; Li, X.; Liu, Y.; Gao, X.; Liang, P. Global identification of microRNAs associated with chlorantraniliprole resistance in diamondback moth Plutella xylostella (L.). Sci. Rep. 2017, 7, 40713. [Google Scholar] [CrossRef]
  31. Etebari, K.; Afrad, M.H.; Tang, B.; Silva, R.; Furlong, M.J.; Asgari, S. Involvement of microRNA miR-2b-3p in regulation of metabolic resistance to insecticides in Plutella xylostella. Insect Mol. Biol. 2018, 27, 478–491. [Google Scholar] [CrossRef] [PubMed]
  32. Shakeel, M.; Xu, X.X.; Xu, J.; Li, S.Z.; Yu, J.L.; Zhou, X.Q.; Xu, X.J.; Hu, Q.B.; Yu, X.Q.; Jin, F.L. Genome-wide identification of Destruxin A-responsive immunity-related microRNAs in diamondback moth, Plutella xylostella. Front. Immunol. 2018, 9, 185. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, J.; Xu, X.X.; Li, S.Z.; Wang, S.; Xu, X.J.; Zhou, X.Q.; Yu, J.L.; Yu, X.Q.; Shakeel, M.; Jin, F.L. Genome-wide profiling of Plutella xylostella immunity-related miRNAs after Isaria fumosorosea infection. Front. Physiol. 2017, 8, 1054. [Google Scholar] [CrossRef] [PubMed]
  34. Zhang, Z.; Jiang, L.; Wang, J.; Gu, P.; Chen, M. MTide: An integrated tool for the identification of miRNA-target interaction in plants. Bioinformatics 2015, 31, 290–291. [Google Scholar] [CrossRef] [PubMed]
  35. Langmead, B.; Trapnell, C.; Pop, M.; Salzberg, S.L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10, 25. [Google Scholar] [CrossRef] [PubMed]
  36. Etebari, K.; Asgari, S. Revised annotation of Plutella xylostella microRNAs and their genome-wide target identification. Insect Mol. Biol. 2016, 25, 788–799. [Google Scholar] [CrossRef] [PubMed]
  37. Fu, W.; Xie, W.; Zhang, Z.; Wang, S.L.; Wu, Q.J.; Liu, Y.; Zhou, X.M.; Zhou, X.G.; Zhang, Y.J. Exploring valid reference genes for quantitative real-time PCR analysis in Plutella xylostella (Lepidoptera: Plutellidae). Int. J. Biol. Sci. 2013, 9, 792–802. [Google Scholar] [CrossRef]
  38. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  39. Schmittgen, T.D.; Livak, K.J. Analyzing real-time PCR data by the comparative CT method. Nat. Protoc. 2008, 3, 1101. [Google Scholar] [CrossRef]
  40. Betel, D.; Wilson, M.; Gabow, A.; Marks, D.S.; Sander, C. The microRNA.org resource: Targets and expression. Nucleic Acids Res. 2008, 36, 149–153. [Google Scholar] [CrossRef]
  41. Kruger, J.; Rehmsmeier, M. RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res. 2006, 34, 451–454. [Google Scholar] [CrossRef] [PubMed]
  42. Tang, W.; Yu, L.; He, W.; Yang, G.; Ke, F.; Baxter, S.W.; You, S.; Douglas, C.J.; You, M. DBM-DB: The diamondback moth genome database. Database 2014, 2014, 87. [Google Scholar] [CrossRef] [PubMed]
  43. He, W.; You, M.; Vasseur, L.; Yang, G.; Xie, M.; Cui, K.; Bai, J.; Liu, C.; Li, X.; Xu, X.; et al. Developmental and insecticide-resistant insights from the de novo assembled transcriptome of the diamondback moth, Plutella xylostella. Genomics 2012, 99, 169–177. [Google Scholar] [CrossRef] [PubMed]
  44. Wang, W.; Scali, M.; Vignani, R.; Spadafora, A.; Sensi, E.; Mazzuca, S.; Cresti, M. Protein extraction for two-dimensional electrophoresis from olive leaf, a plant tissue containing high levels of interfering compounds. Electrophoresis 2003, 24, 2369–2375. [Google Scholar] [CrossRef] [PubMed]
  45. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef] [PubMed]
  46. Wang, Y.; Song, F.; Zhu, J.; Zhang, S.; Yang, Y.; Chen, T.; Tang, B.; Dong, L.; Ding, N.; Zhang, Q.; et al. GSA: Genome sequence archive. Genom. Proteom. Bioinf. 2017, 15, 14–18. [Google Scholar] [CrossRef] [PubMed]
  47. BIG Data Center Members. Database resources of the BIG data center in 2019. Nucleic Acids Res. 2018, 47, 8–14. [Google Scholar]
  48. Kozomara, A.; Birgaoanu, M.; Griffiths-Jones, S. miRBase: From microRNA sequences to function. Nucleic Acids Res. 2019, 47, 155–162. [Google Scholar] [CrossRef] [PubMed]
  49. You, M.; Yue, Z.; He, W.; Yang, X.; Yang, G.; Xie, M.; Zhan, D.; Baxter, S.W.; Vasseur, L.; Gurr, G.M.; et al. A heterozygous moth genome provides insights into herbivory and detoxification. Nat. Genet. 2013, 45, 220–225. [Google Scholar] [CrossRef] [Green Version]
  50. Burmester, T. Origin and evolution of arthropod hemocyanins and related proteins. J. Comp. Physiol. B 2002, 172, 95–107. [Google Scholar]
  51. Ehlting, J.; Chowrira, S.G.; Mattheus, N.; Aeschliman, D.S.; Arimura, G.I.; Bohlmann, J. Comparative transcriptome analysis of Arabidopsis thaliana infested by diamondback moth (Plutella xylostella) larvae reveals signatures of stress response, secondary metabolism, and signalling. BMC Genomics 2008, 9, 154. [Google Scholar] [CrossRef] [PubMed]
  52. Snow, J.W.; Hale, A.E.; Isaacs, S.K.; Baggish, A.L.; Chan, S.Y. Ineffective delivery of diet-derived microRNAs to recipient animal organisms. RNA Biol. 2013, 10, 1107–1116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Li, C.; Wong, A.Y.P.; Wang, S.; Jia, Q.; Chuang, W.P.; Bendena, W.G.; Tobe, S.S.; Yang, S.H.; Chung, G.; Chan, T.F.; et al. miRNA-mediated interactions in and between plants and insects. Int. J. Mol. Sci. 2018, 19, 3239. [Google Scholar] [CrossRef] [PubMed]
  54. Chan, S.Y.; Snow, J.W. Uptake and impact of natural diet-derived small RNA in invertebrates: Implications for ecology and agriculture. RNA Biol. 2017, 14, 402–414. [Google Scholar] [CrossRef] [PubMed]
  55. Karasov, W.H.; Martinez del Rio, C.; Caviedes-Vidal, E. Ecological physiology of diet and digestive systems. Annu. Rev. Physiol. 2011, 73, 69–93. [Google Scholar] [CrossRef]
  56. Huvenne, H.; Smagghe, G. Mechanisms of dsRNA uptake in insects and potential of RNAi for pest control: A review. J. Insect Physiol. 2010, 56, 227–235. [Google Scholar] [CrossRef]
  57. Rera, M.; Clark, R.I.; Walker, D.W. Intestinal barrier dysfunction links metabolic and inflammatory markers of aging to death in Drosophila. Proc. Natl. Acad. Sci. USA 2012, 109, 21528–21533. [Google Scholar] [CrossRef]
  58. Scott, J.G.; Michel, K.; Bartholomay, L.C.; Siegfried, B.D.; Hunter, W.B.; Smagghe, G.; Zhu, K.Y.; Douglas, A.E. Towards the elements of successful insect RNAi. J. Insect Physiol. 2013, 59, 1212–1221. [Google Scholar] [CrossRef] [Green Version]
  59. Terenius, O.; Papanicolaou, A.; Garbutt, J.S.; Eleftherianos, I.; Huvenne, H.; Kanginakudru, S.; Albrechtsen, M.; An, C.; Aymeric, J.L.; Barthel, A.; et al. RNA interference in Lepidoptera: An overview of successful and unsuccessful studies and implications for experimental design. J. Insect Physiol. 2011, 57, 231–245. [Google Scholar] [CrossRef] [Green Version]
  60. Weiberg, A.; Wang, M.; Lin, F.M.; Zhao, H.; Zhang, Z.; Kaloshian, I.; Huang, H.D.; Jin, H. Fungal small RNAs suppress plant immunity by hijacking host RNA interference pathways. Science 2013, 342, 118–123. [Google Scholar] [CrossRef]
  61. Samuel, M.; Bleackley, M.; Anderson, M.; Mathivanan, S. Extracellular vesicles including exosomes in cross kingdom regulation: A viewpoint from plant-fungal interactions. Front. Plant Sci. 2015, 6, 766. [Google Scholar] [CrossRef] [PubMed]
  62. Li, S.; Xu, X.; Zheng, Z.; Zheng, J.; Shakeel, M.; Jin, F. MicroRNA expression profiling of Plutella xylostella after challenge with B. thuringiensis. Dev. Comp. Immunol. 2019, 93, 115–124. [Google Scholar] [CrossRef] [PubMed]
  63. Xie, W.; Melzig, M.F. The stability of medicinal plant microRNAs in the herb preparation process. Molecules 2018, 23, 919. [Google Scholar] [CrossRef] [PubMed]
Biomolecules 09 00602 g001 550
Figure 1. Plant-derived miRNAs in the larval hemolymph of P. xylostella. (A) Comparison of small RNA (sRNA) reads from the two replicates of GC strain libraries including common and specific reads. (B) and (C). sRNA classification (B) and annotation (C) for the two replicates of the GC strain libraries. (D) Length distribution of plant-derived miRNAs. (E) First base bias at 5′-end of plant-derived miRNAs. (F) Quantitative real-time polymerase chain reaction (qRT-PCR) validation of miR159a, miR166a-3p and novel-7703-5p in the hemolymph of the 4th instar larvae of the G88 and GC strains. Data were presented as mean ± SD (Student’s t test, * p < 0.05, *** p < 0.001).
Figure 1. Plant-derived miRNAs in the larval hemolymph of P. xylostella. (A) Comparison of small RNA (sRNA) reads from the two replicates of GC strain libraries including common and specific reads. (B) and (C). sRNA classification (B) and annotation (C) for the two replicates of the GC strain libraries. (D) Length distribution of plant-derived miRNAs. (E) First base bias at 5′-end of plant-derived miRNAs. (F) Quantitative real-time polymerase chain reaction (qRT-PCR) validation of miR159a, miR166a-3p and novel-7703-5p in the hemolymph of the 4th instar larvae of the G88 and GC strains. Data were presented as mean ± SD (Student’s t test, * p < 0.05, *** p < 0.001).
Biomolecules 09 00602 g001
Biomolecules 09 00602 g002 550
Figure 2. Prediction and validation of target genes of plant-derived miRNAs in P. xylostella. (A–C). Putative binding modes of miR159a and Px-BJHSP1 (A), miR166a-3p and Px-BJHSP2 (B), and novel-7703-5p and Px-PPO2 (C) in CDS region. The mfe was indicated above the hybrid. (D–F). Relative expression levels of Px-BJHSP1, Px-BJHSP2 and Px-PPO2 at 24 h post-injection of corresponding miRNA agomir. Data were presented as mean ± SD (Student’s t test, * p < 0.05, ** p < 0.01, *** p < 0.001).
Figure 2. Prediction and validation of target genes of plant-derived miRNAs in P. xylostella. (A–C). Putative binding modes of miR159a and Px-BJHSP1 (A), miR166a-3p and Px-BJHSP2 (B), and novel-7703-5p and Px-PPO2 (C) in CDS region. The mfe was indicated above the hybrid. (D–F). Relative expression levels of Px-BJHSP1, Px-BJHSP2 and Px-PPO2 at 24 h post-injection of corresponding miRNA agomir. Data were presented as mean ± SD (Student’s t test, * p < 0.05, ** p < 0.01, *** p < 0.001).
Biomolecules 09 00602 g002
Biomolecules 09 00602 g003 550
Figure 3. Changes of the protein level of PPO2 in response to the agomir-7703-5p treatment. (A) and (B). Western blots against PPO2 (left) and normalized PPO2 expression levels (right) at 24 h (A) and 48 h (B) post-injection of agomir-7703-5p. Tubulin was used as a loading control. WT: Untreated larvae; agomir-NC: Larvae treated with control miRNA; agomir-7703-5p: Larvae treated with agomir-7703-5p. Integrated densities of protein bands were calculated using ImageJ and PPO2 levels were normalized against tubulin levels. Data were presented as mean ± SD (Student’s t test, n.s. indicated non-significant, ** p < 0.01).
Figure 3. Changes of the protein level of PPO2 in response to the agomir-7703-5p treatment. (A) and (B). Western blots against PPO2 (left) and normalized PPO2 expression levels (right) at 24 h (A) and 48 h (B) post-injection of agomir-7703-5p. Tubulin was used as a loading control. WT: Untreated larvae; agomir-NC: Larvae treated with control miRNA; agomir-7703-5p: Larvae treated with agomir-7703-5p. Integrated densities of protein bands were calculated using ImageJ and PPO2 levels were normalized against tubulin levels. Data were presented as mean ± SD (Student’s t test, n.s. indicated non-significant, ** p < 0.01).
Biomolecules 09 00602 g003
Biomolecules 09 00602 g004 550
Figure 4. Determination of interaction of miRNAs and target binding sites. (A–C). Luciferase assays of miR159a agomir (A), miR166a-3p agomir (B) and novel-7703-5p agomir (C), together with the wild-type (WT) and mutant (Mut) binding sites. Agomir-NC was used as the negative control miRNA. The mean of the relative luciferase activities of the control group was set to 1. Data were presented as mean ± SD (one-way Analysis of Variance (ANOVA) followed by Tukey’s multiple comparisons, different letters above the columns indicated significant difference at p < 0.05).
Figure 4. Determination of interaction of miRNAs and target binding sites. (A–C). Luciferase assays of miR159a agomir (A), miR166a-3p agomir (B) and novel-7703-5p agomir (C), together with the wild-type (WT) and mutant (Mut) binding sites. Agomir-NC was used as the negative control miRNA. The mean of the relative luciferase activities of the control group was set to 1. Data were presented as mean ± SD (one-way Analysis of Variance (ANOVA) followed by Tukey’s multiple comparisons, different letters above the columns indicated significant difference at p < 0.05).
Biomolecules 09 00602 g004
Biomolecules 09 00602 g005 550
Figure 5. Effects of plant-derived miRNAs on the development of P. xylostella. (A) and (B). Larval (A) and pupal (B) durations after miRNA agomir treatments. Data were presented as mean ± SD (Student’s t test, * p < 0.05). (C) Adult emergence rates after miRNA agomir treatments. Data were presented as mean ± SD (Student’s t test, ** p < 0.01) (D) Percentages of abnormalities of G88 pupae after miRNA agomir treatments. (E) Abnormal pupae produced by miRNA agomir treatments compared with the control. Data were presented as mean ± SD (Student’s t test, * p < 0.05).
Figure 5. Effects of plant-derived miRNAs on the development of P. xylostella. (A) and (B). Larval (A) and pupal (B) durations after miRNA agomir treatments. Data were presented as mean ± SD (Student’s t test, * p < 0.05). (C) Adult emergence rates after miRNA agomir treatments. Data were presented as mean ± SD (Student’s t test, ** p < 0.01) (D) Percentages of abnormalities of G88 pupae after miRNA agomir treatments. (E) Abnormal pupae produced by miRNA agomir treatments compared with the control. Data were presented as mean ± SD (Student’s t test, * p < 0.05).
Biomolecules 09 00602 g005
Biomolecules 09 00602 g006 550
Figure 6. Comparison of performance between G88 and GC strains. (A) and (B). Weight of pupa (A) and pupation rate (B). Data were presented as mean ± SD (Student’s t test, * p < 0.05, **** p < 0.0001). (C) and (D). Total number of eggs laid by single pair of male and female adults during the first three days (C) and hatching rate (D). Data were presented as mean ± SD (Student’s t test, **** p < 0.0001).
Figure 6. Comparison of performance between G88 and GC strains. (A) and (B). Weight of pupa (A) and pupation rate (B). Data were presented as mean ± SD (Student’s t test, * p < 0.05, **** p < 0.0001). (C) and (D). Total number of eggs laid by single pair of male and female adults during the first three days (C) and hatching rate (D). Data were presented as mean ± SD (Student’s t test, **** p < 0.0001).
Biomolecules 09 00602 g006
Biomolecules 09 00602 g007 550
Figure 7. Expression of the plant-derived miRNAs and corresponding target genes in P. xylostella after host transfer. (A) Expression of miR159a, miR166a-3p and novel-7703-5p in intact leaves of A. thaliana and leaves fed by the 3rd instar larvae of G88 for 8 h. Data were presented as mean ± SD (Student’s t test, ** p < 0.01 *** p < 0.001,) (B–D). Expression of BJHSP1 (B), BJHSP2 (C) and PPO2 (D) in the 3rd instar larvae, 4th instar larvae and newly-emerged pupae of G88 and GC strains. Data were presented as mean ± SD (Student’s t test, * p < 0.05, ** p < 0.01).
Figure 7. Expression of the plant-derived miRNAs and corresponding target genes in P. xylostella after host transfer. (A) Expression of miR159a, miR166a-3p and novel-7703-5p in intact leaves of A. thaliana and leaves fed by the 3rd instar larvae of G88 for 8 h. Data were presented as mean ± SD (Student’s t test, ** p < 0.01 *** p < 0.001,) (B–D). Expression of BJHSP1 (B), BJHSP2 (C) and PPO2 (D) in the 3rd instar larvae, 4th instar larvae and newly-emerged pupae of G88 and GC strains. Data were presented as mean ± SD (Student’s t test, * p < 0.05, ** p < 0.01).
Biomolecules 09 00602 g007
Table
Table 1. Plant-derived microRNAs (miRNAs) identified in the hemolymph of P. xylostella.
Table 1. Plant-derived microRNAs (miRNAs) identified in the hemolymph of P. xylostella.
NameMature Sequence (5′–3′)Read Count
Rep1Rep2
miR156a-5pUGACAGAAGAGAGUGAGCAC63
miR156b-3pUGCUCACCUCUCUUUCUGUCAGU64
miR157a-5pUUGACAGAAGAUAGAGAGCAC02
miR159aUUUGGAUUGAAGGGAGCUCUA204313
miR159b-3pUUUGGAUUGAAGGGAGCUCUU4280
miR159cUUUGGAUUGAAGGGAGCUCCU1518
miR160a-5pUGCCUGGCUCCCUGUAUGCCA11
miR161.1UGAAAGUGACUACAUCGGGGU13
miR162a-3pUCGAUAAACCUCUGCAUCCAG5058
miR164aUGGAGAAGCAGGGCACGUGCA69
miR164c-5pUGGAGAAGCAGGGCACGUGCG30
miR165a-3pUCGGACCAGGCUUCAUCCCCC13
miR166a-3pUCGGACCAGGCUUCAUUCCCC120151
miR167a-5pUGAAGCUGCCAGCAUGAUCUA2416
miR167dUGAAGCUGCCAGCAUGAUCUGG2213
miR168a-3pUCGCUUGGUGCAGGUCGGGAA49
miR168a-5pCCCGCCUUGCAUCAACUGAAU38
miR171a-3pUGAUUGAGCCGCGCCAAUAUC02
miR172aAGAAUCUUGAUGAUGCUGCAU73
miR172cAGAAUCUUGAUGAUGCUGCAG20
miR319aUUGGACUGAAGGGAGCUCCCU817
miR390a-5pAAGCUCAGGAGGGAUAGCGCC14
miR394aUUGGCAUUCUGUCCACCUCC67
miR395aCUGAAGUGUUUGGGGGAACUC610
miR396a-3pGUUCAAUAAAGCUGUGGGAAG5885
miR396a-5pUUCCACAGCUUUCUUGAACUG2428
miR396b-3pGCUCAAGAAAGCUGUGGGAAA10
miR396b-5pUUCCACAGCUUUCUUGAACUU2329
miR398a-3pUGUGUUCUCAGGUCACCCCUU02
miR398b-3pUGUGUUCUCAGGUCACCCCUG02
miR403-3pUUAGAUUCACGCACAAACUCG78
miR408-3pAUGCACUGCCUCUUCCCUGGC61
miR858bUUCGUUGUCUGUUCGACCUUG1117
novel-2599-3pAUCCGUGGUUUCGCGUAUCGGC04
novel-2783-3pCGCGGAGAAGGGGAAGGGGUGCU07
novel-6634-5pCUGAGAAUUUCUGGAAGAGCUC915
novel-7211-3pUCUCGGACCAGGCUUCAUUCC108
novel-7477-5pGACUGUAAGGCUGUGGAC51
novel-7703-5pCGAGAAUUUCUGGAAGAGCU100101
Table
Table 2. Predicted target genes down-regulated in miRNA agomir treatments.
Table 2. Predicted target genes down-regulated in miRNA agomir treatments.
miRNAGene IDAgomir-NC 1miRNA Agomir 1log2 Ratio 2q-Value 3p-Value 3Annotation (DBM-DB)
miR159aPx00012128.560−5.776.95E-083.88E-08Peroxisomal targeting signal 1 receptor
Px0021981463383−1.878.38E-1373.78E-1384-coumarate-CoA ligase-like 5
Px0027618826−1.691.87E-089.88E-09uncharacterized protein
Px0036121317.25604.38−1.065.86E-536.26E-54Sterol
O-acyltransferase 2
Px00682059809362−7.300Basic juvenile hormone-suppressible protein 1
Px0068971731.5762.4−1.121.42E-751.12E-76Venom serine protease 34
Px00703021806.66205.71−6.6600Basic juvenile hormone-suppressible protein 2
Px00703147614.344669.29−3.2800Basic juvenile hormone-suppressible protein 2
Px0092675411−2.231.22E-077.03E-08Nucleic-acid-binding protein from mobile element jockey
Px01011597.1238.01−1.291.39E-069.05E-07Moricin-like peptide C5
Px01361676.9513.11−2.491.64E-116.62E-12Unknown
Px01581019391−1.021.83E-089.64E-09uncharacterized protein
miR166a-3pPx00703147614.34852.67−5.8300Basic juvenile hormone-suppressible protein 2
novel-7703-5pPx0022744409.621854.67−1.081.34E-1747.24E-176Phenoloxidase subunit 2
Px003522348.6515.57−4.318.01E-731.03E-73Phospholipase A-2-activating protein
Px00384072050.2829765.48−1.100Actin, muscle
Px0048541931.79808.05−1.093.96E-784.69E-79Ecdysteroid UDP-glucosyltransferase
Px00673010017.444288.62−1.0500Pyruvate dehydrogenase E1 component subunit beta, mitochondrial
Px00798972.8329.51−1.131.37E-041.57E-04Elongation factor Tu, mitochondrial
Px00857248906.5318617.36−1.2200Transketolase
Px0090777787.632753.94−1.3300ATP-citrate synthase
1 The agomir-NC and miRNA agomir indicated the gene expression levels (fragments per kilobase of exon per million reads mapped, FPKM) in the control and treatment groups, respectively. 2 Log2 ratio indicates log2 transformation of ratio of FPKM values between miRNA agomir and agomir-NC treatments. 3 The p-value and q-value (the corrected p-value) were considered to be significant at 0.01 and 0.001, respectively

Share and Cite

MDPI and ACS Style

Zhang, L.-l.; Jing, X.-d.; Chen, W.; Wang, Y.; Lin, J.-h.; Zheng, L.; Dong, Y.-h.; Zhou, L.; Li, F.-f.; Yang, F.-y.; et al. Host Plant-Derived miRNAs Potentially Modulate the Development of a Cosmopolitan Insect Pest, Plutella xylostella. Biomolecules 2019, 9, 602. https://0-doi-org.brum.beds.ac.uk/10.3390/biom9100602

AMA Style

Zhang L-l, Jing X-d, Chen W, Wang Y, Lin J-h, Zheng L, Dong Y-h, Zhou L, Li F-f, Yang F-y, et al. Host Plant-Derived miRNAs Potentially Modulate the Development of a Cosmopolitan Insect Pest, Plutella xylostella. Biomolecules. 2019; 9(10):602. https://0-doi-org.brum.beds.ac.uk/10.3390/biom9100602

Chicago/Turabian Style

Zhang, Ling-ling, Xiao-dong Jing, Wei Chen, Yue Wang, Jun-han Lin, Ling Zheng, Yu-hong Dong, Li Zhou, Fei-fei Li, Fei-ying Yang, and et al. 2019. "Host Plant-Derived miRNAs Potentially Modulate the Development of a Cosmopolitan Insect Pest, Plutella xylostella" Biomolecules 9, no. 10: 602. https://0-doi-org.brum.beds.ac.uk/10.3390/biom9100602

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop