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

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