Transcriptome Characterization Analysis of Bactrocera minax and New Insights into Its Pupal Diapause Development with Gene Expression Analysis

Bactrocera minax is a major citrus pest distributed in China, Bhutan and India. The long pupal diapause duration of this fly is a major bottleneck for artificial rearing and underlying mechanisms remain unknown. Genetic information on B. minax transcriptome and gene expression profiles are needed to understand its pupal diapause. High-throughput RNA-seq technology was used to characterize the B. minax transcriptome and to identify differentially expressed genes during pupal diapause development. A total number of 52,519,948 reads were generated and assembled into 47,217 unigenes. 26,843 unigenes matched to proteins in the NCBI database using the BLAST search. Four digital gene expression (DGE) libraries were constructed for pupae at early diapause, late diapause, post-diapause and diapause terminated developmental status. 4,355 unigenes showing the differences expressed across four libraries revealed major shifts in cellular functions of cell proliferation, protein processing and export, metabolism and stress response in pupal diapause. When diapause was terminated by 20-hydroxyecdysone (20E), many genes involved in ribosome and metabolism were differentially expressed which may mediate diapause transition. The gene sets involved in protein and energy metabolisms varied throughout early-, late- and post-diapause. A total of 15 genes were selected to verify the DGE results through quantitative real-time PCR (qRT-PCR); qRT-PCR expression levels strongly correlated with the DGE data. The results provided the extensive sequence resources available for B. minax and increased our knowledge on its pupal diapause development and they shed new light on the possible mechanisms involved in pupal diapause in this species.