Genome-wide analysis of developmental stage-specific transcriptome in Bradysia odoriphaga

Bradysia odoriphaga is a serious pest of the Chinese chive; however, detailed information regarding the developmental stage-specific gene expression patterns of B. odoriphaga is not yet available. In this study, RNA sequencing (RNA-seq) was performed to determine the gene expression patterns of developmental stages including the eggs, second instar larvae, fourth instar larvae, pupae, and adults of B. odoriphaga. Analysis of 15 samples revealed an average of 89.56% of the clean reads could be mapped onto the assembled UniGene database. Cluster tree analysis showed that the expression patterns were stage-specific and that samples of the second and fourth instar larvae clustered in one group, while those of eggs, pupae, and adults clustered in another group. Differential expression unigenes (DEUs) for sequential developmental stages were between 3314 and 10,632. A total of 1910-7756 DEUs of sequential developmental stages were assigned into 45-56 gene ontology categories and 1165-3845 DEUs were mapped into Kyoto Encyclopedia of Genes and Genomes pathways. The expression of DEUs related to growth and development showed that hormone receptors highly expressed in the pupal stage, while chitinases were highly expressed in the larval stage. The results of quantitative real time polymerase chain reaction (qRT-PCR) and RNA-seq expression agreed well for 12 growth- and development-related unigenes. This study identified DEUs for sequential developmental stages of B. odoriphaga. Gene Ontology classifications and KEGG pathway identification of DEUs not only provide information useful for understanding insect growth and development but also for exploring novel approaches to control B. odoriphaga.