Identification of differentially expressed microRNAs between Bacillus thuringiensis Cry1Ab-resistant and -susceptible strains of Ostrinia furnacalis

The Asian corn borer (ACB), Ostrinia furnacalis (Guenée), can develop strong resistance to Cry1Ab, the most widely commercialized Cry toxin for Bt maize worldwide. It is essential to understand the mechanism of resistance for management of this species, but information on the post-transcriptional regulation of Bt resistance in this target insect is limited. In the present study, RNA was extracted from the ACB in various larval stages (1–5 instar) from Cry1Ab-sensitive (ACB-BtS) and -resistant (ACB-AbR) strains, each of which included two biological replicates. Using Illumina sequencing, a total of 23,809,890 high-quality reads were collected from the four ACB libraries. The numbers of known microRNAs (miRNAs) were 302 and 395 for ACB-BtS and 268 and 287 for ACB-AbR. Using Mireap software, we identified 32 and 16 potential novel miRNAs for ACB-BtS and 18 and 22 for ACB-AbR. Among them, 21 known and 1 novel miRNAs had significantly different expression between ACB-BtS and ACB-AbR. Several miRNAs were observed to target potential Bt receptor genes, such as aminopeptidase N and cadherin-like protein. The glycosylphosphatidylinositol-anchor biosynthetic process and ABC transporters pathway were identified through Gene Ontology and KEGG pathway analysis of target genes of the differentially expressed miRNAs.