摘要:
[目的]探明氯虫苯甲酰胺处理后桃小食心虫(Carposina sasakii)成虫的转录组差异,了解该药剂影响桃小食心虫交配的基因在功能分类和代谢通路等方面的生物学特征,挖掘与交配相关的功能基因.[方法]通过生物学实验观察氯虫苯甲酰胺干扰桃小食心虫成虫交配及繁殖情况.采用Illumina HiSeqTM2500高通量测序技术对刚羽化的桃小食心虫雌雄成虫、羽化后4-6 h进入交配高峰期的雌雄成虫和经氯虫苯甲酰胺处理4-6 h的雌雄成虫进行转录组测序,利用Trinity软件对所得序列进行de novo组装及评估,之后对获得的有效序列进行功能注释,并利用qRT-PCR技术分析氯虫苯甲酰胺处理后相关基因的时空表达变化.[结果]氯虫苯甲酰胺处理后,桃小食心虫的交配率显著降低,寿命缩短,产卵量减少.通过合并组装桃小食心虫转录组有效序列共获得102 831条unigene,其中34 526个有注释信息.根据筛选标准,氯虫苯甲酰胺处理过程中,雌雄虫中分别有122个和147个基因发生变化,其中相同差异基因31个.对所存在的234个差异基因进行GO功能注释和富集分析结果显示,分子功能过程中的催化活性和结合活性及生物学过程中与代谢过程、单一生物体过程和细胞过程相关的5类基因占主导地位.KEGG分类结果显示,富集到代谢通路的最多,有25个,包括昆虫激素合成、药物代谢等.通过比对分析,鉴定c64662.graph_c0为桃小食心虫鱼尼丁受体基因,其长度为15 637 bp,与已报道CsRyR的一致性为99.0%.此外,在234个差异基因中,鉴别羧酸酯酶unigene 3个、细胞色素P450 unigene 4个、肌钙蛋白unigene 3个、气味结合蛋白unigene 1个和生物钟unigene 1个.细胞色素P450 c40709.graph_c0参与昆虫激素生物合成.根据转录组中基因表达分析,12个基因在雌雄虫中的表达均出现不同变化趋势.qRT-PCR结果显示,氯虫苯甲酰胺诱导羧酸酯酶c51998.graph_c0基因上调表达;药剂处理后,雌雄虫的c57480.graph_c0和c53794.graph_c0的表达分别呈现显著的上调和下调变化,而c40709.graph_c0基因仅在雄虫中显著下调,处理6 h后c53281.graph_c0只在雌虫中上调表达;氯虫苯甲酰胺处理后3个肌钙蛋白基因在整个试验阶段均表现明显的下调趋势;雄虫中的生物钟unigene c60883.graph_c0和气味结合蛋白unigene c45675.graph_c0的变化趋势较为一致,进入暗期后立即上调,但均受氯虫苯甲酰胺的抑制.雌虫体内RyR的表达量显著上调,而雄虫初次到达求偶高峰期前RyR的表达量与对照差异不显著,到第2个求偶高峰期时,RyR表达显著下调.除细胞色素P450c57480.graph_c0、c40709.graph_c0和c53281.graph_c0外,其余基因在雄虫中的表达均高于雌虫.且触角酯酶基因c54944.graph_c0、生物钟基因c60883.graph_c0和气味结合蛋白基因c45675.graph_c0在黑暗光照交替变化时,表达发生明显地上调或下调.[结论]通过转录组测序发现氯虫苯甲酰胺干扰桃小交配的作用机制是由靶标基因、嗅觉相关基因、代谢基因、生物钟基因等相互作用引起的.