1. Introduction
MicroRNAs (miRNAs) are endogenous, single-stranded, non-protein-coding small RNAs approximately 19 to 24 nucleotides in length [1]. These miRNAs are ubiquitously present in animals, plants, and microorganisms, playing a crucial role in the post-transcriptional regulation of genes [2,3,4,5,6,7,8]. Studies have demonstrated that miRNAs are integral to various biological processes, including tissue growth, germ cell development, hormone action, and the development and function of the central nervous system, primarily through gene regulation [9,10]. It is estimated that miRNAs regulate the expression of more than 50% of protein-coding genes in animals [11].
Research on insect miRNAs has garnered increasing attention, leading to the discovery of numerous miRNAs across different insect species. These miRNAs regulate a wide range of physiological functions throughout insect development, such as molting, metamorphosis, oogenesis, embryogenesis, behavior, and host–pathogen interactions [12,13,14,15,16,17,18]. For instance, Dicer-1, a key gene in miRNA processing, when interfered with RNA, was found to partially inhibit the metamorphosis of Blattella germanica (L.) (Blattaria: Blattidae) after molting [4]. The miR-2 family influences insect metamorphosis and development via juvenile hormones [19,20,21], while miR-14 regulates insect molting by participating in the regulation of ecdysone receptors [22,23]. Additionally, miRNAs play vital roles in the growth and development of other insects such as Bactrocera dorsalis Hendel (Diptera: Tephritidae), Nilaparvata lugens (Hemiptera: Delphacidae), and Spodoptera exigua (Lepidoptera: Noctuidae) [24,25,26].
Given their significant role in post-transcriptional gene regulation, binding to specific sequences of target mRNA, miRNA can degrade target mRNA or inhibit its normal translation, thereby regulating the expression of target genes. miRNAs are also implicated in the molecular mechanisms that enhance drug resistance by regulating detoxification genes in insects [27,28]. Detoxification metabolic resistance is crucial for pest resistance development. For instance, miR-7a and miR-8519 could upregulate the expression of the ryanodine receptor (RyR), enhancing resistance to chlorantraniliprole [29]. miR-998-3p can negatively regulate the expression of ABCC2, contributing to the detoxification of BtCry1Ac toxins in three typical lepidopteran pests: Helicoverpa armigera (Lepidoptera: Noctuidae), S. exigua, and Plutella xylostella (Lepidoptera: Plutellidae) [30]. Furthermore, miRNAs novel-85 and novel-191 can target CYP6ER1 and carboxylesterase 1 (CarE1), significantly altering N. lugens’ susceptibility to nitenpyram [31]. In addition, miR-4133-3p could regulate the detoxification of gossypol and tannic acid by targeting the cytochrome P450 4CJ1 (CYP4CJ1) in Aphis gossypii Glover (Hemiptera: Aphididae) [32].
The fall armyworm, S. frugiperda Smith, 1797 (Lepidoptera: Noctuidae), is a major migratory agricultural pest of global concern [33]. This omnivorous pest has a strong migratory ability, a high reproductive rate, and a short life cycle, damaging 353 species of plants across 76 families, including major crops such as corn, sorghum, sugarcane, barley, rice, pepper, wild oat, and potato [34,35]. The primary method of controlling this pest relies on chemical pesticides, which have proven to be effective [36]. However, reports indicate that S. frugiperda has developed resistance to various insecticides, including bromides, neonicotinoids, and pyrethroids [37,38,39]. Insect resistance to insecticides is mainly attributed to enhanced metabolism and alterations in target site structures [40,41,42].
To date, some miRNAs have been identified in S. frugiperda, with functional studies primarily focusing on antiviral immune defense and adaptive evolution [43,44]. However, systematic identification and functional analysis of miRNAs related to detoxification metabolism in S. frugiperda remain limited [45]. Using high-throughput sequencing of miRNA libraries, this study investigated miRNAs in S. frugiperda at different developmental stages, analyzed differentially expressed miRNAs in larvae, and identified target genes and functional predictions of these miRNAs across four stages. This research aims to elucidate the regulatory mechanisms of S. frugiperda growth and development, analyze detoxification-related target genes, and provide new strategies for the effective management of S. frugiperda outbreaks and damage.
2. Materials and Methods
2.1. Insect Collection and Rearing
The larvae of S. frugiperda were collected from a cornfield at the Xindu Base of the Sichuan Academy of Agricultural Sciences, China. They were continuously raised in the insect laboratory of Plant Conservation Institute of Sichuan Academy of Agricultural Sciences. The larvae were reared continuously in a light incubator under controlled conditions (16/8 hr light/dark (L/D) cycle, 65% humidity, 25 °C) and fed with an artificial diet [46], which consisted of the following components: wheat germ powder 280 g, soy protein powder 90 g, yeast powder 35 g, agar 25 g, vitamin B complex 0.2 g, cholesterol 12 g, sorbate 2 g, ascorbic acid 12 g, methyl p-hydroxybenzoate 5 g, formaldehyde 4 mL, penicillin 0.2 g, and distilled water 1500 mL. Sampling periods at different developmental stages were as follows: 500 eggs were collected within 24 h post-laying, 20–30 third instar larvae were collected within 12 h post-molting, 3–5 pupae were collected within 36 h post-pupation, and 2–3 adults were collected within 12 h post-emergence. Each sample contained 0.3–0.5 g and had three biological replicates. Samples were placed in 2.0 mL RNase-free centrifuge tubes with TRIzol reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) for temporary storage. After collection, samples were stored in an ultra-low temperature freezer at −80 °C.
2.2. Small RNA Library Preparation and Sequencing
Total RNA was extracted using TRIzol reagent (Invitrogen Life Technologies, Carlsbad, CA, USA) following the manufacturer’s instructions. RNA concentrations were measured using a NanoDrop spectrophotometer (Thermo Scientific, Sunnyvale, CA, USA), and RNA quality and integrity were assessed with an Agilent 2100 Bioanalyzer (Agilent Technologies, Waldbronn, Germany). Total RNA was used to construct small RNA (sRNA) libraries with the NEBNext Multiplex Small RNA Library Prep Set for Illumina (New England Biolabs, Ipswich, MA, USA).
One microgram of total RNA from S. frugiperda samples was ligated to 3′ and 5′ adapters using Ligation Enzyme Mix. The samples were then reverse-transcribed using Superscript II reverse transcriptase. The sRNA libraries were subjected to quality control, and the average insert size was approximately 140 to 150 bp. cDNA was prepared using random primers and a cDNA Synthesis Kit (Invitrogen, CA, USA). The sequencing library was quantified using the Agilent High Sensitivity DNA assay on an Agilent Bioanalyzer (Agilent Technologies, Germany) and sequenced on a NovaSeq 6000 (Illumina) in PE150 mode at Personal Biotechnology Co., Ltd. (Shanghai, China).
2.3. Unigene Assembly, Analysis, and Annotation
Raw sequencing reads underwent quality control using FastQC program (Babraham Bioinformatics, Cambridge, UK) [47], and then the raw data were filtered using the Personalbio company’s self-developed script. Clean reads ranging from 18 to 36 nt were filtered, and deduplication was performed to obtain unique reads for subsequent analysis. The reference genome index was built using Bowtie2 v2.5.1 [48], and de-duplicated clean reads were mapped to the reference genome using miRDeep2 v2.0 (
2.4. Differential Expression Analysis of miRNAs
The reads count of miRNAs was determined based on sequences aligned to mature miRNAs. The highest abundance among miRNAs of the same name was used for subsequent analysis. Differentially expressed miRNAs were identified using DESeq v1.39.0 [50], with transcripts showing a log2 fold change > 1 and p-value < 0.05 considered significant. Bidirectional cluster analysis of all miRNAs and samples was performed using the R package Pheatmap v1.0.12 [51], with the Euclidean method for distance calculation and Complete Linkage for hierarchical clustering.
2.5. Quantitative Reverse Transcription PCR
Eight differentially expressed miRNAs were selected for qRT-PCR using initial RNA samples. Stem-loop RT primers and gene-specific primers were designed (Table 1). Reverse transcription was performed with the PrimeScript RT reagent kit with gDNA Eraser (Takara, Tokyo, Japan), the reference gene Actin reverse transcription was performed using the general primers of the kit, and the specific primers of miRNAs were reversed. We used 1 μg of total RNA under the following conditions: 16 °C for 30 min, 42 °C for 30 min, and 85 °C for 5 min. The reference gene qRT-PCR reaction contained 2 μL of diluted cDNA (15 ng), 10 μL of 2× AceQ qPCR SYBR Green Master mix (Vazyme, Nanjing, China), and 10 μM of each primer in a 20 μL total volume [52]. The miRNAs qRT-PCR reaction contained 2 μL of diluted cDNA (15 ng), 10 μL of 2× AceQ Universal U+ Probe Master mix (Vazyme, China), 10 μM probe, and 10 μM of each primer in a 20 μL total volume [53]. The conditions used were: 95 °C for 5 min, followed by 40 cycles of 95 °C for 10 s, 60 °C for 30 s, and a melting curve from 68 °C to 95 °C. The Applied Biosystems QuantStudio 6 Flex system (Thermo Scientific, Sunnyvale, CA, USA) was used with Sf-actin as an endogenous control. The 2−ΔΔCT method was used to calculate the relative expression levels of miRNAs [54]. The qRT-PCR included three technical and biological replicates.
2.6. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Enrichment Analysis
Target genes of differentially expressed miRNAs were predicted using miRanda v3.3a (
2.7. Statistical and Data Analysis
Statistical analyses were conducted using Statistical Product and Service Solutions (SPSS) v20.0 (SPSS Inc., Chicago, IL, USA). Comparisons were performed using Student’s t-test, with a p-value < 0.05 considered statistically significant.
3. Results
3.1. Small RNA Sequencing Data in S. frugiperda
High-throughput sequencing generated twelve sRNA libraries for S. frugiperda at four different developmental stages: egg, larva, pupa, and adult. The transcriptome data for these stages contain 20,128,612, 16,315,920, 23,520,239, and 17,593,809 raw reads, respectively (Table 2). The total raw sequences for the four developmental stages exceeded 10 Mb, with high-quality reads constituting 91–97%, indicating good sequencing quality. After filtering clean reads (≥18 nt), the number of valid sequences obtained for egg, larva, pupa and adult stages were 17,231,905, 3,731,980, 17,141,575, and 12,107,130 reads, respectively.
The sRNA lengths for the four developmental stages of S. frugiperda ranged from 18 to 36 nt (Figure 1). In larvae, pupae, and adults, the 22 nt sequences were the most abundant, accounting for 43.20%, 48.18%, and 38.99% of the total, respectively. In eggs, the 27 nt sequences were the most abundant, making up 30.72% of the total.
3.2. Small RNA Classification and Annotation
To identify known non-coding RNAs, the precursor and mature sequences of miRNAs for the species were downloaded from miRBase. The deduplicated sequences were aligned to these miRNAs to annotate sRNAs at different developmental stages of S. frugiperda (Figure 2A, Table S1). Among the annotated sRNAs, ribosomal RNA (rRNA) was the most abundant, constituting 32.89% of the total. Larvae and pupae had higher rRNA quantities than adults and eggs. tRNA, snoRNA, and snRNA were significant components of sRNAs, with maximum proportions of 0.40%, 0.12%, and 0.14%, respectively. A significant proportion of sRNAs (31.94–81.26%) remained unannotated across all developmental stages.
The nucleotide preference at the first base of miRNAs in different developmental stages of S. frugiperda was analyzed (Figure 2B), with most sRNAs showing a preference for U (uracil) at the first base. In all samples, the first base of 18 bp sRNAs preferred A (adenine). Different nucleotide positions exhibited different preferences (Figure 2C). The preference for U at the 5′ end is a conserved characteristic of miRNAs.
3.3. Identification of Known and Novel miRNAs in S. frugiperda
Precursor and mature sequences of miRNAs were downloaded from miRBase, and deduplicated sequences were aligned to these miRNAs for annotation. A total of 135, 149, 168, and 164 known miRNAs were identified in eggs, larvae, pupae, and adults, respectively. Among these, 115 known miRNAs were present across all four developmental stages (Figure 3A). Additionally, novel miRNAs were predicted from sequences not annotated with any information using mireap analysis, resulting in 147, 83, 122, and 138 novel miRNAs identified in eggs, larvae, pupae, and adults, respectively. Among these, 32 novel miRNAs were common across all four stages (Figure 3B). The ten most highly expressed miRNAs were detected (Table S2), with sfr-miR-2766-3p, sfr-miR-279a-3p, and sfr-miR-10-5p being the most abundantly expressed.
3.4. Identification of Differentially Expressed miRNAs in Larvae
miRNAs expressed in both larvae and other developmental stages (eggs, pupae, and adults) were selected based on differential expression (log2 fold change > 1) and significance (p-value < 0.05). The results showed 67 differentially expressed miRNAs across the four developmental stages. The number of differentially expressed miRNAs between larvae and eggs, pupae, and adults were 54, 15, and 18, respectively (Table 3).
Eight miRNAs showed differential expression in S. frugiperda larvae compared to the other three stages (Figure 4). Among them, sfr-miR-14-5p, sfr-miR-6094-5p, and sfr-miR-6094-3p were significantly upregulated in larvae compared to the other stages. The most significantly upregulated miRNA in larvae was sfr-miR-6094-3p, with expression levels increasing by 223.41, 482.67, and 79.76 times compared to eggs, pupae, and adults, respectively. sfr-miR-6094-5p and sfr-miR-14-5p were also highly upregulated, with sfr-miR-6094-5p showing increases of 209.54, 253.92, and 87.57 times, and sfr-miR-14-5p showing increases of 23.89, 3.69, and 5.21 times, compared to eggs, pupae, and adults, respectively.
Conversely, the most significantly downregulated miRNA in larvae was sfr-miR-2765-5p, with expression decreasing by 74.20, 61.11, and 37.06 times compared to eggs, pupae, and adults, respectively. sfr-miR-279b-3p was also significantly downregulated, with decreases of 3.89, 3.77, and 6.97 times compared to eggs, pupae, and adults, respectively. Additionally, sfr-miR-277-3p and sfr-miR-307-3p were downregulated in larvae compared to pupae and adults, and sfr-miR-34-5p was downregulated in larvae compared to adults.
3.5. qRT-PCR Validation of Differentially Expressed miRNAs in Larvae
To further confirm the miRNA sequencing results, eight miRNAs (miR-6094-5p, miR-6094-3p, miR-14-5p, miR-10505-3p, miR-2765-5p, miR-277-3p, miR-307-3p, and miR-34-5p) were randomly selected for qRT-PCR based on their differential expression in larvae compared to the other three developmental stages. The qRT-PCR results showed consistent expression patterns with RNA-seq. Compared with the egg stage, the expression of miR-2765-5p was downregulated, and the expression of other miRNAs was upregulated (Figure 5A). Compared with the pupa stage, the expressions of miR-6094-5p, miR-6094-3p, miR-10505-3p, miR-14-5p, and miR-34-5p were upregulated, while the expressions of miR-2765-5p, miR-277-3p, and miR-307-3p were downregulated (Figure 5B). Compared with the adult stage, the expressions of miR-6094-5p, miR-6094-3p, and miR-14-5p were upregulated, while the expressions of miR-2765-5p, miR-277-3p, miR-307-3p, and miR-34-5p were downregulated (Figure 5C). This suggests that the orientation of regulation and expression pattern of these miRNAs validated by RT-qPCR was consistent with the results from the miRNA sequencing, thus confirming the reliability and repeatability of the miRNA sequencing method.
3.6. Prediction of Targeted Genes
Using miRanda, the 3′ UTR sequences of the species’ mRNAs were targeted for the prediction of differentially expressed miRNA target genes. A total of 2019 target genes and 66,004 target sites were predicted for the 67 differentially expressed miRNAs in S. frugiperda larvae.
To further investigate the functions of miRNAs in the S. frugiperda transcriptome, the top twenty significantly enriched KEGG pathways were analyzed. In larvae vs. adults (Figure 6A), the most significantly enriched pathways were chemokine signaling pathway, gonadotropin-releasing hormone (GnRH) secretion, cyclic adenosine monophosphate (cAMP) signaling pathway, GnRH signaling pathway, and axon guidance. In larvae vs. eggs (Figure 6B), the top pathways were endocytosis, apelin signaling pathway, chemokine signaling pathway, circadian entrainment, and phospholipase D signaling pathway. In larvae vs. pupae (Figure 6C), the most significantly enriched pathways were oxytocin signaling pathway, vascular smooth muscle contraction, serotonergic synapse, axon guidance, and GnRH secretion.
GO enrichment analysis of the target genes of differentially expressed miRNAs showed that the functions of differentially expressed genes in larvae compared to the other three stages were mainly distributed across cellular components, molecular functions, and biological processes. In larvae, adults, and pupae, the G protein-coupled receptor signaling pathway was a significantly enriched subgroup (Figure 7 and Figure 8). In larvae and eggs, localization was a significantly enriched subgroup (Figure 9).
To better understand the interactions between miRNAs and their potential target genes, eight of the most significantly differentially expressed miRNAs and detoxification metabolism-related target mRNAs were selected to construct a miRNA–mRNA network (Figure 10). The eight miRNAs included the four most upregulated miRNAs in larvae (sfr-miR-14-5p, sfr-miR-6094-5p, sfr-miR-6094-3p, and sfr-miR-34-5p) and the four most downregulated miRNAs (sfr-miR-2765-5p, sfr-miR-279b-3p, sfr-miR-277-3p, and sfr-miR-307-3p). A total of 82 target genes were identified, targeting 10 gene families involved in detoxification metabolism. These included twenty-nine genes in the P450 gene family, twenty-four in the ABC transporters gene family, thirteen in the GSTs gene family, four in the CarEs gene family, nine in the Acetylcholine gene family, twenty-three in the Sodium channel gene family, seven in the Ryanodine gene family, eight in the Cadherin gene family, seven in the Alkaline phosphatase gene family, and eleven in the Aminopeptidase gene family. These results suggest that detoxification metabolism-related pathways play a critical role in insecticide resistance in S. frugiperda.
4. Discussion
In recent years, increasing numbers of miRNAs have been discovered in various eukaryotic organisms [57]. In insects, miRNAs regulate gene expression at the post-transcriptional level by either repressing translation or degrading mRNA, with their functions mainly focused on various physiological processes in insects [58,59,60,61]. The larval stage is the primary period during which S. frugiperda infests crops like corn [62]. It is also a crucial period for rapid growth and preparation for pupation. Studying miRNAs preferentially expressed in the larvae of the fall armyworm can aid in identifying miRNAs involved in regulating growth and development, as well as analyzing detoxification-related target genes, providing new insights and methods for effectively controlling outbreaks and damage caused by this pest.
Chemical control of S. frugiperda has a long history, so a variety of resistance mechanisms have emerged [63]. In the United States and Brazil, chemical control of S. frugiperda is mainly used before planting Bt transgenic maize [39]. The metabolic capacity of carbaryl -resistant strains is five times higher than that of sensitive strains, and the resistance is mainly caused by oxidative metabolism such as P450 hydroxylation and epoxidation [64]. The activities of various detoxification metabolic enzymes such as multifunctional oxidase (MFO), GSTs, and esterase (ESTs) of laboratory-resistant strains are significantly higher than those of sensitive strains [37]. The sensitivity of acetylcholinesterase (AChE), a molecular target of carbamate and organophosphorus insecticides, to menafarb is also significantly reduced [65]. All these indicate that the resistance of S. frugiperda in the field population is caused by several mechanisms, so the study of insecticide resistance is particularly important in chemical control.
This study also provides comprehensive insights into the sRNA landscape of S. frugiperda across its developmental stages, with a particular focus on identifying and predicting differentially expressed miRNAs associated with detoxification pathways in larvae. The results revealed critical patterns and potential regulatory roles of miRNAs, which could be crucial for understanding and managing this pest’s resistance to insecticides.
Our analysis identified a substantial number of known and novel miRNAs across the developmental stages of S. frugiperda, including 135, 149, 168, and 164 known miRNAs in eggs, larvae, pupae, and adults, respectively. Additionally, we discovered 147, 83, 122, and 138 novel miRNAs in these stages. These findings are consistent with previous studies that have identified numerous miRNAs in various insect species [31,43,66,67,68,69], underscoring the conserved nature of miRNA-mediated regulation in insect development and physiology.
Notably, miR-2766-3p, miR-279a-3p, and miR-10-5p were among the most abundantly expressed miRNAs in S. frugiperda, with significant expression in the larval stage. These miRNAs have been implicated in crucial biological processes such as development, stress response, and resistance mechanisms in other insect species. For instance, miR-2766 regulates tyrosine hydroxylase in H. armigera, influencing larval–pupal metamorphosis [70]. Similarly, miR-279a-3p modulates CYP325BB1 expression, affecting insecticide resistance in mosquitoes [71]. In addition, miR-10 is also the most abundantly expressed miRNA in larvae of H. armigera when parasitized by the wasp Diadegma semiclausum [72].
The differential expression analysis highlighted eight miRNAs with significant changes in expression levels in S. frugiperda larvae compared to other developmental stages. These include upregulated miRNAs such as sfr-miR-6094-3p, sfr-miR-6094-5p, and sfr-miR-34-5p, and downregulated miRNAs such as sfr-miR-2765-5p and sfr-miR-279b-3p. The high expression of sfr-miR-34-5p in larvae and its decrease in eggs and pupae suggest its pivotal role in larval development. This pattern aligns with previous findings in other insects, where miRNAs were also highly expressed in larvae compared to other stages [73,74,75,76,77,78].
The KEGG and GO enrichment analyses of target genes predicted for differentially expressed miRNAs revealed their involvement in critical biological pathways, including detoxification metabolism and G protein-coupled receptor (GPCR) signaling. GPCR signaling has been linked to insecticide resistance in various insects, including S. frugiperda [79]. Our study identified 82 potential detoxification-related target genes regulated by differentially expressed miRNAs, including genes from the P450, GST, and ABC transporter families. These gene families are well-known for their roles in insecticide detoxification and resistance [80].
Previous studies have demonstrated that miRNAs can modulate insecticide sensitivity by targeting detoxification genes. For instance, miR-2b-3p and miR-14-5p in P. xylostella larvae were shown to influence the expression of CYP9F2 and CYP307a1, affecting the larval detoxification capacity [81]. Similarly, regulating novel_miR-1517 expression in Bemisia tabaci affected its sensitivity to imidacloprid by modulating CYP6CM1 expression [82]. The chlorpyrifos- and cypermethrin-resistant nature of S. frugiperda strains showed that the A201S, G227A, and F290V point mutations of AChE are resistant to chlorpyrifos [39]. The T929I, L932F, and L1014F point mutations in sodium ion channels lead to resistance to cyhalothrin. The increased expressions of P450, GSTs, and ESTs are involved in resistance to chlorpyrifos and cypermethrin [83]. Our findings suggest that miRNAs in S. frugiperda may similarly regulate detoxification genes, contributing to its resistance to various insecticides.
The findings of this study have significant implications for understanding the molecular mechanisms underlying insecticide resistance in S. frugiperda. The identified miRNAs and their target genes involved in detoxification pathways offer potential targets for developing novel pest management strategies. By manipulating specific miRNAs or their targets, it might be possible to enhance the efficacy of insecticides or develop new control methods that circumvent existing resistance mechanisms. Future research should focus on functional validation of these miRNAs and their targets using techniques such as miRNA mimics and inhibitors, gene knockout or knockdown experiments, and biochemical assays to confirm their roles in detoxification and resistance. Additionally, exploring the environmental and physiological factors that influence miRNA expression and activity will further elucidate their regulatory networks in S. frugiperda.
5. Conclusions
In summary, eight miRNAs showed differential expression in S. frugiperda larvae compared to the other three stages. A total of 82 target genes were identified, targeting 10 gene families involved in detoxification metabolism. These results provide a theoretical basis for further exploring the mechanism of miRNAs regulating detoxification metabolism genes.
Conceptualization, H.-L.L. and S.C.; methodology, C.-Y.Y. and C.-H.L.; software, C.-Y.Y. and C.-H.L.; validation, C.-C.Z. and Z.-Z.W.; formal analysis, Y.-P.W. and X.-Y.C.; investigation, J.-W.G. and W.-J.Y.; data curation, Y.-P.W.; writing—original draft preparation, Y.-P.W. and X.-Y.C.; writing—review and editing, D.-Q.P.; project administration, H.-L.L. and D.-Q.P.; funding acquisition, H.-L.L. and Y.-P.W. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
The original data presented in the study are openly available in FigShare at
We thank the anonymous reviewers for valuable comments.
The authors declare no conflicts of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Length distribution of sRNAs in different stages of S. frugiperda. X-axis indicates sequence length, while Y-axis shows the frequence in percentage.
Enlarge this image.Figure 2. Composition of sRNAs in S. frugiperda. (A) Annotation of sRNAs. X-axis shows the sample number, whereas Y-axis shows the proportion of the de-duplicated sequences annotated to each sRNA to the total de-duplicated sequences. (B) First base bias of the miRNAs in S. frugiperda. X-axis shows the miRNAs length, and Y-axis shows the frequency in percentage. (C) Base bias of the miRNAs at each position. X-axis shows the nucleotide position, and Y-axis shows the frequency in percentage.
Enlarge this image.Figure 2. Composition of sRNAs in S. frugiperda. (A) Annotation of sRNAs. X-axis shows the sample number, whereas Y-axis shows the proportion of the de-duplicated sequences annotated to each sRNA to the total de-duplicated sequences. (B) First base bias of the miRNAs in S. frugiperda. X-axis shows the miRNAs length, and Y-axis shows the frequency in percentage. (C) Base bias of the miRNAs at each position. X-axis shows the nucleotide position, and Y-axis shows the frequency in percentage.
Enlarge this image.Figure 3. Venn diagram of differential expression in known (A) and novel (B) miRNAs of S. frugiperda. Numbers indicate shared miRNAs.
Enlarge this image.Figure 4. Heatmap of differentially expressed miRNAs in different stages of S. frugiperda. Red indicates upregulation; blue means downregulation.
Enlarge this image.Figure 5. Relative expression level of selected differentially expressed miRNAs of S. frugiperda between qRT-PCR and RNA-seq in (A) larvae vs. eggs, (B) larvae vs. pupae, and (C) larvae vs. adults.
Enlarge this image.Figure 6. The top 20 significantly enriched KEGG pathways of S. frugiperda in (A) larvae vs. adults, (B) larvae vs. eggs, and (C) larvae vs. pupae. X-axis shows the enrichment factor, and Y-axis shows the pathway names. Definitions of the pathways are available at https://www.kegg.jp/kegg/pathway.html. The enrichment degree was measured by Rich factor value and the number of miRNA target genes enriched on this pathway. Rich factor refers to the ratio of the number of differential miRNA target genes enriched in this pathway to the number of differential miRNA target genes annotated. The greater the Rich factor, the greater the degree of enrichment.
Enlarge this image.Figure 6. The top 20 significantly enriched KEGG pathways of S. frugiperda in (A) larvae vs. adults, (B) larvae vs. eggs, and (C) larvae vs. pupae. X-axis shows the enrichment factor, and Y-axis shows the pathway names. Definitions of the pathways are available at https://www.kegg.jp/kegg/pathway.html. The enrichment degree was measured by Rich factor value and the number of miRNA target genes enriched on this pathway. Rich factor refers to the ratio of the number of differential miRNA target genes enriched in this pathway to the number of differential miRNA target genes annotated. The greater the Rich factor, the greater the degree of enrichment.
Enlarge this image.Figure 7. The top 30 GO terms of the target genes of differentially expressed miRNAs of S. frugiperda in larvae vs. adults. The enrichment degree was measured using the false discovery rate (FDR) value and the number of miRNA target genes enriched on this pathway. FDR generally ranges from 0 to 1, and the closer it is to zero, the more significant the enrichment. The top 30 pathways with the smallest FDR value, that is, the most significant enrichment, were selected for display.
Enlarge this image.Figure 8. The top 30 GO terms of the target genes of differentially expressed miRNAs of S. frugiperda in larvae vs. pupae. The enrichment degree was measured by FDR value and the number of miRNA target genes enriched on this pathway. FDR generally ranges from 0 to 1, and the closer it is to zero, the more significant the enrichment. The top 30 pathways with the smallest FDR value, that is, the most significant enrichment, were selected for display.
Enlarge this image.Figure 9. The top 30 GO terms of the target genes of differentially expressed miRNAs of S. frugiperda in larvae vs. eggs. X-axis shows GO terms, and Y-axis shows −log10 (p-value). The enrichment degree was measured by FDR value and the number of miRNA target genes enriched on this pathway. FDR generally ranges from 0 to 1, and the closer it is to zero, the more significant the enrichment. The top 30 pathways with the smallest FDR value, that is, the most significant enrichment, were selected for display.
Enlarge this image.Figure 10. A miRNA–mRNA network of 8 differentially expressed miRNAs and their target genes related to detoxification metabolism in S. frugiperda. Red triangles represent upregulated miRNAs while blue triangles represent downregulated miRNAs. Other circles represent different target genes related to detoxification metabolism.
Primers used for the qRT-PCR.
| Primer | Sequence (from 5′ End to 3′ End) |
|---|---|
| miR-6094-5P-RT | GTCGTATCCAGTGCAGGGTCCGAGGTAT |
| miR-6094-5P-F | TCAGCGGTGGCCTGGG |
| miR-6094-3P-RT | GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACAGGATC |
| miR-6094-3P-F | CGCGTATTCGAGACCTCTGCT |
| miR-10505-3p-RT | GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACGAGCCA |
| miR-10505-3p-F | CGCGCGTAGGGTTAGAAACT |
| miR-14-5p-RT | GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACTCGAGT |
| miR-14-5p-F | GCGGGGGGAGGAATTG |
| miR-2765-5p-RT | GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACGCCAAC |
| miR-2765-5p-F | CGCGCGTGGTAACTCCACCACC |
| miR-277-3p-RT | GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACTGTCGT |
| miR-277-3p-F | CGCGCGTAAATGCACTATCTGGT |
| miR-307-3p-RT | GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACGCTCAC |
| miR-307-3p-F | CGCGCGTCACAACCTCCTTGA |
| miR-34-5p-RT | GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACACAACC |
| miR-34-5p-F | CGCGCGTGGCAGTGTGGTTAGCT |
| miR-R | AGTGCAGGGTCCGAGGTATT |
| Sf-actin(nei)-F | CGAGAAGATGACCCAGAT |
| Sf-actin(nei)-R | GATAGCACAGCCTGGATA |
| probe | FAM-CGCACTGGATACGAC-MGB |
Summary of sequencing reads.
| Sample | Flag | Raw Reads | Average Raw Reads | Clean Reads | Average Clean Reads |
|---|---|---|---|---|---|
| Egg-1 | E01 | 23,489,347 | 20,128,612 | 20,154,610 | 17,231,905 |
| Egg-2 | E02 | 24,772,159 | 20,128,612 | 20,586,799 | 17,231,905 |
| Egg-3 | E03 | 12,124,331 | 20,128,612 | 10,954,307 | 17,231,905 |
| Larva-1 | L01 | 10,837,072 | 16,315,920 | 1,502,965 | 3,731,980 |
| Larva-2 | L02 | 17,172,714 | 16,315,920 | 2,494,144 | 3,731,980 |
| Larva-3 | L03 | 20,937,973 | 16,315,920 | 7,198,831 | 3,731,980 |
| Pupa-1 | P01 | 20,479,951 | 23,520,239 | 14,980,746 | 17,141,575 |
| Pupa-2 | P02 | 25,798,866 | 23,520,239 | 19,742,017 | 17,141,575 |
| Pupa-3 | P03 | 24,281,901 | 23,520,239 | 16,701,963 | 17,141,575 |
| Adult-1 | A01 | 23,403,358 | 17,593,809 | 13,367,017 | 12,107,130 |
| Adult-2 | A02 | 12,435,180 | 17,593,809 | 10,254,039 | 12,107,130 |
| Adult-3 | A03 | 16,942,888 | 17,593,809 | 12,700,333 | 12,107,130 |
The number of differentially expressed miRNAs in different stages of S. frugiperda.
| Control | Case | Upregulated Genes | Downregulated Genes | Total DEGs |
|---|---|---|---|---|
| Egg | Larva | 32 | 22 | 54 |
| Adult | Larva | 9 | 9 | 18 |
| Pupa | Larva | 6 | 9 | 15 |
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Ambros, V. The functions of animal microRNAs. Nature; 2004; 431, pp. 350-355. [DOI: https://dx.doi.org/10.1038/nature02871] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15372042]
2. Asgari, S. Venom proteins from polydnavirus-producing endoparasitoids: Their role in host-parasite interactions. Arch. Insect Biochem.; 2006; 61, pp. 146-156. [DOI: https://dx.doi.org/10.1002/arch.20109]
3. Ibanez-Ventoso, C.; Vora, M.; Driscoll, M. Sequence relationships among C. elegans, D. melanogaster and human microRNAs highlight the extensive conservation of microRNAs in biology. PLoS ONE; 2008; 3, e2818. [DOI: https://dx.doi.org/10.1371/journal.pone.0002818] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18665242]
4. Gomez-Orte, E.; Belles, X. MicroRNA-dependent metamorphosis in hemimetabolan insects. Proc. Natl. Acad. Sci. USA; 2009; 106, pp. 21678-21682. [DOI: https://dx.doi.org/10.1073/pnas.0907391106]
5. Asgari, S. Role of MicroRNAs in insect host-microorganism interactions. Front. Physiol.; 2011; 2, 48. [DOI: https://dx.doi.org/10.3389/fphys.2011.00048]
6. Freitak, D.; Knorr, E.; Vogel, H.; Vilcinskas, A. Gender- and stressor-specific microRNA expression in Tribolium castaneum. Biol. Lett.; 2012; 8, pp. 860-863. [DOI: https://dx.doi.org/10.1098/rsbl.2012.0273]
7. Fullaondo, A.; Lee, S.Y. Identification of putative miRNA involved in Drosophila melanogaster immune response. Dev. Comp. Immunol.; 2012; 36, pp. 267-273. [DOI: https://dx.doi.org/10.1016/j.dci.2011.03.034] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21641929]
8. Yu, X.; Zhou, Q.; Cai, Y.; Luo, Q.; Lin, H.; Hu, S.; Yu, J. A discovery of novel microRNAs in the silkworm (Bombyx mori) genome. Genomics; 2009; 94, pp. 438-444. [DOI: https://dx.doi.org/10.1016/j.ygeno.2009.08.007]
9. Legeai, F.; Rizk, G.; Walsh, T.; Edwards, O.; Gordon, K.; Lavenier, D.; Leterme, N.; Mereau, A.; Nicolas, J.; Tagu, D. et al. Bioinformatic prediction, deep sequencing of microRNAs and expression analysis during phenotypic plasticity in the pea aphid, Acyrthosiphon pisum. BMC Genom.; 2010; 11, 281. [DOI: https://dx.doi.org/10.1186/1471-2164-11-281]
10. Kloosterman, W.P.; Plasterk, R.H. The diverse functions of microRNAs in animal development and disease. Dev. Cell; 2006; 11, pp. 441-450. [DOI: https://dx.doi.org/10.1016/j.devcel.2006.09.009]
11. Krol, J.; Loedige, I.; Filipowicz, W. The widespread regulation of microRNA biogenesis, function and decay. Nat. Rev. Genet.; 2010; 11, pp. 597-610. [DOI: https://dx.doi.org/10.1038/nrg2843]
12. Zhang, Y.; Zhou, X.; Ge, X.; Jiang, J.; Li, M.; Jia, S.; Yang, X.; Kan, Y.; Miao, X.; Zhao, G. et al. Insect-Specific microRNA involved in the development of the silkworm Bombyx mori. PLoS ONE; 2009; 4, e4677. [DOI: https://dx.doi.org/10.1371/journal.pone.0004677]
13. Etebari, K.; Asgari, S. Conserved microRNA miR-8 blocks activation of the Toll pathway by upregulating Serpin 27 transcripts. RNA Biol.; 2013; 10, pp. 1356-1364. [DOI: https://dx.doi.org/10.4161/rna.25481]
14. Lucas, K.J.; Zhao, B.; Liu, S.; Raikhel, A.S. Regulation of physiological processes by microRNAs in insects. Curr. Opin. Insect Sci.; 2015; 11, pp. 1-7. [DOI: https://dx.doi.org/10.1016/j.cois.2015.06.004]
15. Ling, L.; Kokoza, V.A.; Zhang, C.; Aksoy, E.; Raikhel, A.S. MicroRNA-277 targets insulin-like peptides 7 and 8 to control lipid metabolism and reproduction in Aedes aegypti mosquitoes. Proc. Natl. Acad. Sci. USA; 2017; 114, pp. E8017-E8024. [DOI: https://dx.doi.org/10.1073/pnas.1710970114]
16. Zhang, Z.; Liu, X.; Shiotsuki, T.; Wang, Z.; Xu, X.; Huang, Y.; Li, M.; Li, K.; Tan, A. Depletion of juvenile hormone esterase extends larval growth in Bombyx mori. Insect Biochem. Mol. Biol.; 2017; 81, pp. 72-79. [DOI: https://dx.doi.org/10.1016/j.ibmb.2017.01.001]
17. Song, J.; Li, W.; Zhao, H.; Gao, L.; Fan, Y.; Zhou, S. The microRNAs let-7 and miR-278 regulate insect metamorphosis and oogenesis by targeting the juvenile hormone early-response gene Kruppel-homolog 1. Development; 2018; 145, 170670. [DOI: https://dx.doi.org/10.1242/dev.170670]
18. Shang, F.; Niu, J.; Ding, B.Y.; Zhang, W.; Wei, D.D.; Wei, D.; Jiang, H.B.; Wang, J.J. The miR-9b microRNA mediates dimorphism and development of wing in aphids. Proc. Natl. Acad. Sci. USA; 2020; 117, pp. 8404-8409. [DOI: https://dx.doi.org/10.1073/pnas.1919204117]
19. Lin, L. Functional Analysis of microRNA let-7 and miR-2 Regulating Development in Bombyx mori. Ph.D. Thesis; Northwestern Polytechnical University: Xi’an, China, 2017.
20. Wang, Y.L.; Yang, M.L.; Jiang, F.; Zhang, J.Z.; Kang, L. MicroRNA-dependent development revealed by RNA interference-mediated gene silencing of LmDicer1 in the migratory locust. Insect Sci.; 2013; 20, pp. 53-60. [DOI: https://dx.doi.org/10.1111/j.1744-7917.2012.01542.x]
21. Lozano, J.; Montanez, R.; Belles, X. MiR-2 family regulates insect metamorphosis by controlling the juvenile hormone signaling pathway. Proc. Natl. Acad. Sci. USA; 2015; 112, pp. 3740-3745. [DOI: https://dx.doi.org/10.1073/pnas.1418522112]
22. Chen, X.; Fu, J. The microRNA miR-14 regulates egg-laying by targeting EcR in honeybees (Apis mellifera). Insects; 2021; 12, 351. [DOI: https://dx.doi.org/10.3390/insects12040351]
23. He, K.; Xiao, H.; Sun, Y.; Ding, S.; Situ, G.; Li, F. Transgenic microRNA-14 rice shows high resistance to rice stem borer. Plant Biotechnol. J.; 2019; 17, pp. 461-471. [DOI: https://dx.doi.org/10.1111/pbi.12990]
24. Puthiyakunnon, S.; Yao, Y.; Li, Y.; Gu, J.; Peng, H.; Chen, X. Functional characterization of three MicroRNAs of the Asian tiger mosquito, Aedes albopictus. Parasite Vector; 2013; 6, 230. [DOI: https://dx.doi.org/10.1186/1756-3305-6-230]
25. Peng, W.; Zheng, W.W.; Tariq, K.; Yu, S.N.; Zhang, H.Y. MicroRNA Let-7 targets the ecdysone signaling pathway E75 gene to control larval-pupal development in Bactrocera dorsalis. Insect Sci.; 2019; 26, pp. 229-239. [DOI: https://dx.doi.org/10.1111/1744-7917.12542]
26. Yang, M.; Wei, Y.; Jiang, F.; Wang, Y.; Guo, X.; He, J.; Kang, L. MicroRNA-133 inhibits behavioral aggregation by controlling dopamine synthesis in locusts. PLoS Genet.; 2014; 10, e1004206. [DOI: https://dx.doi.org/10.1371/journal.pgen.1004206]
27. Zha, W.; Chang, M.; You, A. Regulation of microRNA on insect growth and development. Hubei Agric. Sci.; 2019; 58, pp. 198-202. [DOI: https://dx.doi.org/10.14088/j.cnki.issn0439-8114.2019.22.045]
28. Deng, C.; He, Q. Advanced research of microRNA in insects. J. Heilongjiang Bayi Agric. Univ.; 2021; 33, pp. 15-21.
29. Li, X.; Guo, L.; Zhou, X.; Gao, X.; Liang, P. miRNAs regulated overexpression of ryanodine receptor is involved in chlorantraniliprole resistance in Plutella xylostella (L.). Sci. Rep.; 2015; 5, 14095. [DOI: https://dx.doi.org/10.1038/srep14095]
30. Zhu, B.; Sun, X.; Nie, X.; Liang, P.; Gao, X. MicroRNA-998-3p contributes to Cry1Ac-resistance by targeting ABCC2 in lepidopteran insects. Insect. Biochem. Mol. Biol.; 2020; 117, 103283. [DOI: https://dx.doi.org/10.1016/j.ibmb.2019.103283]
31. Mao, K.; Jin, R.; Ren, Z.; Zhang, J.; Li, Z.; He, S.; Ma, K.; Wan, H.; Li, J. miRNAs targeting CYP6ER1 and CarE1 are involved in nitenpyram resistance in Nilaparvata lugens. Insect Sci.; 2022; 29, pp. 177-187. [DOI: https://dx.doi.org/10.1111/1744-7917.12910]
32. Ma, K.; Li, F.; Tang, Q.; Liang, P.; Liu, Y.; Zhang, B.; Gao, X. CYP4CJ1-mediated gossypol and tannic acid tolerance in Aphis gossypii Glover. Chemosphere; 2019; 219, pp. 961-970. [DOI: https://dx.doi.org/10.1016/j.chemosphere.2018.12.025]
33. Van den Berg, J.; Brewer, M.J.; Reisig, D.D. A special collection: Spodoptera frugiperda (fall armyworm): Ecology and management of its world-scale invasion outside of the Americas. J. Econ. Entomol.; 2022; 115, pp. 1725-1728. [DOI: https://dx.doi.org/10.1093/jee/toac143]
34. Goergen, G.; Kumar, P.L.; Sankung, S.B.; Togola, A.; Tamo, M. First report of outbreaks of the fall armyworm Spodoptera frugiperda (JE Smith) (Lepidoptera, Noctuidae), a new alien invasive pest in west and central Africa. PLoS ONE; 2016; 11, e165632. [DOI: https://dx.doi.org/10.1371/journal.pone.0165632]
35. Cock, M.; Beseh, P.K.; Buddie, A.G.; Cafa, G.; Crozier, J. Molecular methods to detect Spodoptera frugiperda in Ghana, and implications for monitoring the spread of invasive species in developing countries. Sci. Rep.; 2017; 7, 4103. [DOI: https://dx.doi.org/10.1038/s41598-017-04238-y]
36. Belay, D.K.; Huckaba, R.M.; Foster, J.E. Susceptibility of the fall armyworm, Spodoptera frugiperda (lepidoptera: Noctuidae), at Santa Isabel, Puerto Rico, to different insecticides. Fla. Entomol.; 2012; 95, pp. 476-478. [DOI: https://dx.doi.org/10.1653/024.095.0232]
37. Yu, S.J.; McCord, E.J. Lack of cross-resistance to indoxacarb in insecticide-resistant Spodoptera frugiperda (Lepidoptera: Noctuidae) and Plutella xylostella (Lepidoptera: Yponomeutidae). Pest Manag. Sci.; 2007; 63, pp. 63-67. [DOI: https://dx.doi.org/10.1002/ps.1309] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17089332]
38. Sharon, G.; Segal, D.; Ringo, J.M.; Hefetz, A.; Zilber-Rosenberg, I.; Rosenberg, E. Commensal bacteria play a role in mating preference of Drosophila melanogaster. Proc. Natl. Acad. Sci. USA; 2010; 107, pp. 20051-20056. [DOI: https://dx.doi.org/10.1073/pnas.1009906107]
39. Carvalho, R.A.; Omoto, C.; Field, L.M.; Williamson, M.S.; Bass, C. Investigating the molecular mechanisms of organophosphate and pyrethroid resistance in the fall armyworm Spodoptera frugiperda. PLoS ONE; 2013; 8, e62268. [DOI: https://dx.doi.org/10.1371/journal.pone.0062268]
40. Nascimento, A.R.B.D.; Farias, J.R.; Bernardi, D.; Horikoshi, R.J.; Omoto, C. Genetic basis of Spodoptera frugiperda (Lepidoptera: Noctuidae) resistance to the chitin synthesis inhibitor lufenuron. Pest Manag. Sci.; 2016; 72, pp. 810-815. [DOI: https://dx.doi.org/10.1002/ps.4057]
41. Bolzan, A.; Padovez, F.E.; Nascimento, A.R.; Kaiser, I.S.; Lira, E.C.; Amaral, F.S.; Kanno, R.H.; Malaquias, J.B.; Omoto, C. Selection and characterization of the inheritance of resistance of Spodoptera frugiperda (Lepidoptera: Noctuidae) to chlorantraniliprole and cross-resistance to other diamide insecticides. Pest Manag. Sci.; 2019; 75, pp. 2682-2689. [DOI: https://dx.doi.org/10.1002/ps.5376]
42. Lira, E.C.; Bolzan, A.; Nascimento, A.R.; Amaral, F.S.; Kanno, R.H.; Kaiser, I.S.; Omoto, C. Resistance of Spodoptera frugiperda (Lepidoptera: Noctuidae) to spinetoram: Inheritance and cross-resistance to spinosad. Pest Manag. Sci.; 2020; 76, pp. 2674-2680. [DOI: https://dx.doi.org/10.1002/ps.5812]
43. Mone, Y.; Nhim, S.; Gimenez, S.; Legeai, F.; Seninet, I.; Parrinello, H.; Negre, N.; D’Alencon, E. Characterization and expression profiling of microRNAs in response to plant feeding in two host-plant strains of the lepidopteran pest Spodoptera frugiperda. BMC Genom.; 2018; 19, 804. [DOI: https://dx.doi.org/10.1186/s12864-018-5119-6]
44. Karamipour, N.; Fathipour, Y.; Talebi, A.A.; Asgari, S.; Mehrabadi, M. The microRNA pathway is involved in Spodoptera frugiperda (Sf9) cells antiviral immune defense against Autographa californica multiple nucleopolyhedrovirus infection. Insect Biochem. Mol. Biol.; 2019; 112, 103202. [DOI: https://dx.doi.org/10.1016/j.ibmb.2019.103202]
45. Mahalle, R.M.; Sun, W.; Posos-Parra, O.A.; Jung, S.; Mota-Sanchez, D.; Pittendrigh, B.R.; Seong, K.M. Identification of differentially expressed miRNAs associated with diamide detoxification pathways in Spodoptera frugiperda. Sci. Rep.; 2024; 14, 4308. [DOI: https://dx.doi.org/10.1038/s41598-024-54771-w]
46. He, L.; Wang, T.; Chen, Y.; Ge, S.; Wyckhuys, K.; Wu, K. Larval diet affects development and reproduction of East Asian strain of the fall armyworm, Spodoptera frugiperda. J. Integr. Agric.; 2021; 20, pp. 736-744. [DOI: https://dx.doi.org/10.1016/S2095-3119(19)62879-0]
47. Andrews, S.A. Quality Control Tool for High Throughput Sequence Data. Available online: http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/ (accessed on 7 December 2016).
48. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with bowtie 2. Nat. Methods; 2012; 9, pp. 357-359. [DOI: https://dx.doi.org/10.1038/nmeth.1923] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22388286]
49. Friedlander, M.R.; Chen, W.; Adamidi, C.; Maaskola, J.; Einspanier, R.; Knespel, S.; Rajewsky, N. Discovering microRNAs from deep sequencing data using miRDeep. Nat. Biotechnol.; 2008; 26, pp. 407-415. [DOI: https://dx.doi.org/10.1038/nbt1394]
50. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol.; 2010; 11, R106. [DOI: https://dx.doi.org/10.1186/gb-2010-11-10-r106]
51. Hu, K. Become competent in generating RNA-seq heat maps in one day for novices without prior R experience. Methods Mol. Biol.; 2021; 2239, pp. 269-303. [DOI: https://dx.doi.org/10.1007/978-1-0716-1084-8_17]
52. Peng, X.; Nguyen, A.; Ghosh, D. Quantification of M13 and T7 bacteriophages by Taqman and SYBR green qPCR. J. Virol. Methods; 2018; 252, pp. 100-107. [DOI: https://dx.doi.org/10.1016/j.jviromet.2017.11.012]
53. Sriram, K.; Wiley, S.Z.; Moyung, K.; Gorr, M.W.; Salmerón, C.; Marucut, J.; French, R.P.; Lowy, A.M.; Insel, P.A. Detection and quantification of GPCR mRNA: An assessment and implications of data from high-content methods. ACS Omega; 2019; 4, pp. 17048-17059. [DOI: https://dx.doi.org/10.1021/acsomega.9b02811]
54. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods; 2001; 25, pp. 402-408. [DOI: https://dx.doi.org/10.1006/meth.2001.1262]
55. Yang, K.; Sablok, G.; Qiao, G.; Nie, Q.; Wen, X. IsomiR2Function: An integrated workflow for identifying microRNA variants in plants. Front. Plant Sci.; 2017; 8, 322. [DOI: https://dx.doi.org/10.3389/fpls.2017.00322]
56. Yu, G.; Wang, L.G.; Han, Y.; He, Q.Y. Clusterprofiler: An R package for comparing biological themes among gene clusters. Omics; 2012; 16, pp. 284-287. [DOI: https://dx.doi.org/10.1089/omi.2011.0118]
57. Asgari, S. MicroRNA functions in insects. Insect Biochem. Mol. Biol.; 2013; 43, pp. 388-397. [DOI: https://dx.doi.org/10.1016/j.ibmb.2012.10.005]
58. Belles, X. MicroRNAs and the evolution of insect metamorphosis. Annu. Rev. Entomol.; 2017; 62, pp. 111-125. [DOI: https://dx.doi.org/10.1146/annurev-ento-031616-034925] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27813669]
59. Lucas, K.; Raikhel, A.S. Insect microRNAs: Biogenesis, expression profiling and biological functions. Insect Biochem. Mol. Biol.; 2013; 43, pp. 24-38. [DOI: https://dx.doi.org/10.1016/j.ibmb.2012.10.009]
60. Azzam, G.; Smibert, P.; Lai, E.C.; Liu, J.L. Drosophila Argonaute 1 and its miRNA biogenesis partners are required for oocyte formation and germline cell division. Dev. Biol.; 2012; 365, pp. 384-394. [DOI: https://dx.doi.org/10.1016/j.ydbio.2012.03.005]
61. Wightman, B.; Ha, I.; Ruvkun, G. Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans. Cell; 1993; 75, pp. 855-862. [DOI: https://dx.doi.org/10.1016/0092-8674(93)90530-4]
62. Gouin, A.; Bretaudeau, A.; Nam, K.; Gimenez, S.; Aury, J.M.; Duvic, B.; Hilliou, F.; Durand, N.; Montagne, N.; Darboux, I. et al. Two genomes of highly polyphagous lepidopteran pests (Spodoptera frugiperda, Noctuidae) with different host-plant ranges. Sci. Rep.; 2017; 7, 11816. [DOI: https://dx.doi.org/10.1038/s41598-017-10461-4]
63. Rios-Diez, J.D.; Saldamando-Benjumea, C. Susceptibility of Spodoptera frugiperda (lepidoptera: Noctuidae) strains from central Colombia to two insecticides, methomyl and lambda-cyhalothrin: A study of the genetic basis of resistance. J. Econ. Entomol.; 2011; 104, pp. 1698-1705. [DOI: https://dx.doi.org/10.1603/EC11079] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22066201]
64. Mccord, E.; Yus, J. The mechanisms of carbaryl resistance in the fall armyworm, Spodoptera frugiperda (J.E. Smith). Pestic. Biochem. Physiol.; 1987; 27, pp. 114-122. [DOI: https://dx.doi.org/10.1016/0048-3575(87)90103-9]
65. Gutiérrez-Moreno, R.; Mota-Sanchez, D.; A Blanco, C.; E Whalon, M.; Terán-Santofimio, H.; Rodriguez-Maciel, J.C.; DiFonzo, C. Field-evolved resistance of the fall armyworm (Lepidoptera: Noctuidae) to synthetic insecticides in Puerto Rico and Mexico. J. Econ. Entomol.; 2019; 112, pp. 792-802. [DOI: https://dx.doi.org/10.1093/jee/toy372] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30535077]
66. Yang, Y.; Zhang, Y.; Wang, A.; Duan, A.; Xue, C.; Wang, K.; Zhao, M.; Zhang, J. Four microRNAs, miR-13b-3p, miR-278-5p, miR-10483-5p, and miR-10485-5p, mediate insecticide tolerance in Spodoptera frugiperda. Front. Genet.; 2021; 12, 820778. [DOI: https://dx.doi.org/10.3389/fgene.2021.820778] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35126473]
67. 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. [DOI: https://dx.doi.org/10.1038/srep40713]
68. Li, J.M.; Zhou, Y.R.; Sun, Z.T.; Wang, X.; Xie, L.; Chen, J.P. Identification and profiling of conserved and novel microRNAs in Laodelphax striatellus in response to rice black-streaked dwarf virus (RBSDV) infection. Genom. Data; 2015; 3, pp. 63-69. [DOI: https://dx.doi.org/10.1016/j.gdata.2014.08.010] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26484150]
69. Zhang, B.Z.; Hu, G.L.; Lu, L.Y.; Hu, S.F.; Li, Y.S.; Su, X.; Dong, W.Y.; Zhen, C.A.; Liu, R.Q.; Kong, F.B. et al. Identification of differentially expressed microRNAs under imidacloprid exposure in Sitobion miscanthi. Pestic. Biochem. Phys.; 2021; 177, 104885. [DOI: https://dx.doi.org/10.1016/j.pestbp.2021.104885]
70. Shen, Z.J.; Zhu, F.; Liu, Y.J.; Li, Z.; Moural, T.W.; Liu, X.M.; Liu, X. MicroRNAs miR-14 and miR-2766 regulate tyrosine hydroxylase to control larval-pupal metamorphosis in Helicoverpa armigera. Pest Manag. Sci.; 2022; 78, pp. 3540-3550. [DOI: https://dx.doi.org/10.1002/ps.6997]
71. Li, X.; Hu, S.; Zhang, H.; Yin, H.; Wang, H.; Zhou, D.; Sun, Y.; Ma, L.; Shen, B.; Zhu, C. MiR-279-3p regulates deltamethrin resistance through CYP325BB1 in Culex pipiens pallens. Parasite Vector; 2021; 14, 528. [DOI: https://dx.doi.org/10.1186/s13071-021-05033-5]
72. Etebari, K.; Hussain, M.; Asgari, S. Identification of microRNAs from Plutella xylostella larvae associated with parasitization by Diadegma semiclausum. Insect Biochem. Mol. Biol.; 2013; 43, pp. 309-318. [DOI: https://dx.doi.org/10.1016/j.ibmb.2013.01.004]
73. Bai, J.; Lin, G.; Wang, Y.; Zhao, X.; Peng, L.; Xie, M.; Yang, J.; He, W.; You, M. Identification and functional characterization of larval-based microRNAs in the diamondback moth, Plutella xylostella. J. Biosaf.; 2022; 31, pp. 64-74. [DOI: https://dx.doi.org/10.3969/j.issn.2095-1787.2022.01.009]
74. Yang, S. Function Analysis of Protein Coding Gene and miRNA in Pupation and Pupae Development of in Chilo suppressalis. Ph.D. Thesis; Nanjing Agricultural University: Nanjing, China, 2015.
75. Liang, Z.; Yang, Y.; Sun, X.; Du, J.; Wang, Q.; Zhang, G.; Zhang, J.; Yin, X.; Singh, D.; Su, P. et al. Integrated analysis of microRNA and mRNA expression profiles in the fat bodies of MbMNPV-Infected Helicoverpa armigera. Viruses; 2022; 15, 19. [DOI: https://dx.doi.org/10.3390/v15010019] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36680059]
76. Wang, C.; Guo, X.; Li, Y.; Zhang, J.; Fu, Y. miR-34-5p, encoded by Spodoptera frugiperda, participates in anti-baculovirus by regulating innate immunity in the insect host. Int. J. Biol. Macro. Mol.; 2022; 222, pp. 2190-2199. [DOI: https://dx.doi.org/10.1016/j.ijbiomac.2022.09.293]
77. Liu, Z.; Xu, J.; Ling, L.; Luo, X.; Yang, D.; Yang, X.; Zhang, X.; Huang, Y. miR-34 regulates larval growth and wing morphogenesis by directly modulating ecdysone signaling and cuticle protein in Bombyx mori. RNA Biol.; 2020; 17, pp. 1342-1351. [DOI: https://dx.doi.org/10.1080/15476286.2020.1767953]
78. Wu, Y.; Guo, Y.; Fan, X.; Zhao, H.; Zhang, Y.; Guo, S.; Jing, X.; Liu, Z.; Feng, P.; Liu, X. et al. ame-miR-34 modulates the larval body weight and immune response of Apis mellifera Workers to Ascosphara apis Invasion. Int. J. Mol. Sci.; 2023; 24, 1214. [DOI: https://dx.doi.org/10.3390/ijms24021214] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36674732]
79. Li, T.; Liu, N. Role of the G-protein-coupled receptor signaling pathway in insecticide resistance. Int. J. Mol. Sci.; 2019; 20, 4300. [DOI: https://dx.doi.org/10.3390/ijms20174300] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31484301]
80. Hilliou, F.; Chertemps, T.; Maïbèche, M.; Le Goff, G. Resistance in the genus Spodoptera: Key insect detoxification genes. Insects; 2021; 12, 544. [DOI: https://dx.doi.org/10.3390/insects12060544] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34208014]
81. 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, pp. 478-491. [DOI: https://dx.doi.org/10.1111/imb.12387]
82. Gong, P.P.; Wei, X.G.; Liu, S.N.; Yang, J.; Fu, B.L.; Liang, J.J.; Huang, M.J.; Du, T.H.; Yin, C.; Ji, Y. et al. Novel_miR-1517 mediates CYP6CM1 to regulate imidacloprid resistance in Bemisia tabaci (Hemiptera: Gennadius). Pestic. Biochem. Phys.; 2023; 194, 105469. [DOI: https://dx.doi.org/10.1016/j.pestbp.2023.105469]
83. Giraudo, M.; Hilliou, F.; Fricaux, T.; Audant, P.; Feyereisen, R.; Le Goff, G. Cytochrome p450s from the fall armyworm (Spodoptera frugiperda): Responses to plant allelochemicals and pesticides. Insect Mol. Biol.; 2015; 24, pp. 115-128. [DOI: https://dx.doi.org/10.1111/imb.12140]
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Chen, Song 1 ; Hong-Ling, Liu 1 1 Key Laboratory of Integrated Pest Management of Southwest Crops, Institute of Plant Protection, Sichuan Academy of Agricultural Sciences, Chengdu 610066, China;
2 Science and Technology Security Center, Sichuan Academy of Agricultural Sciences, Chengdu 610066, China;