旨在筛选野桑蚕蛹滞育的调控关联基因,为在分子水平揭示该昆虫蛹滞育调控机制提供理论基础。以野桑蚕蛹滞育不同时期的全蛹样本为研究对象,采用高通量测序平台Illumina HiSeqTM 2000进行转录组测序分析;利用KAAS在线pathway比对分析工具对筛选出的满足padj <0.05且fold change≥32或padj <0.05且fold change≤1/32条件的差异表达基因进行KEGG通路富集分析;利用stem软件对差异表达基因进行时间序列分析。选取前滞育期和滞育解除期低表达,滞育时期高表达的基因为蛹滞育关联基因,并进行GO和KEGG富集分析。结果表明,在滞育前期、滞育中期、滞育晚期和滞育解除期分别检测到13 764、13 765、13 765和13 764个表达转录本。样本组的两两比较发现差异表达基因4 621个。与滞育前期相比,滞育中期、滞育晚期和解除滞育期的上调表达基因数量分别是1 921、451和38个,下调表达的基因数量分别是1143、171和59个。利用时间序列分析方法在差异表达基因中发现290个蛹滞育关联基因,GO富集分析将202个基因映射到52项二级GO类目,KEGG通路富集分析将67个基因映射到126个三级代谢通路。本研究为后续的野桑蚕蛹滞育调控关键基因的筛选和功能研究提供了理论基础。