Time-Resolved Whole-Transcriptome Analysis Suggests Candidate Non-Coding RNA Regulatory Networks Associated with PBAN-Induced Pheromone Biosynthesis in Ostrinia furnacalis

The biosynthesis of sex pheromones in lepidopteran pheromone glands is tightly regulated by pheromone biosynthesis-activating neuropeptide (PBAN) signaling; yet the contribution of non-coding RNA-mediated post-transcriptional regulation remains largely unclear. This study aimed to characterize temporal transcriptomic changes, candidate non-coding RNA-mediated regulatory associations, and temporal molecular dynamics underlying transcriptional remodeling after PBAN treatment in Ostrinia furnacalis. First, we performed comprehensive whole-transcriptome sequencing (WTS) on 18 biologically independent samples collected at six time points (0, 20, 40, 60, 90, and 120 min) after PBAN injection. Then, we systematically identified and quantified the dynamic expression patterns of differentially expressed (DE) mRNAs, miRNAs, lncRNAs, and circRNAs in response to PBAN stimulation. By integratively analyzing these multidimensional omics datasets and inferring sequence-based interaction relationships, we inferred a dynamic candidate competing endogenous RNA (ceRNA) like regulatory network. The candidate ceRNA network anchored four core node genes: the PBAN receptor (PBANR), the rate-limiting enzyme acetyl-CoA carboxylase (ACC), and the terminal biosynthetic enzymes desaturase (DES) and fatty acyl-CoA reductase (FAR). The qRT-PCR results further support the temporal expression pattern of key genes during the PBAN response, suggesting that this network can provide a valuable resource for further functional studies.