西方蜜蜂工蜂蛹响应低温胁迫的差异表达基因分析

【目的】本研究旨在从转录组水平筛选和分析西方蜜蜂Apis mellifera工蜂不同时期蛹对低温胁迫反应的差异表达基因(DEGs)。【方法】将西方蜜蜂工蜂封盖后3 d预蛹以及封盖后6和9 d蛹分别置于低温环境(20℃)和最适发育温度(35℃)中4 h，分别作为处理组和对照组，通过转录组学技术筛选低温处理组与对应的对照组之间的DEGs，并进行GO功能分类和KEGG通路分析。利用RT qPCR分别对封盖后3 d预蛹以及封盖后6和9 d蛹的8, 6和5个DEGs的表达谱进行验证。【结果】与对照组相比，封盖后3 d预蛹以及封盖后6和9 d蛹受低温胁迫后的DEGs分别有220, 50和26个；GO功能分类发现，DEGs富集数最多的条目为代谢进程、细胞进程、催化活性和结合，封盖后3 d预蛹的DEGs在生物学进程调控、细胞部分和细胞器等有较多富集。KEGG通路分析显示，处理组和对照组间西方蜜蜂各日龄DEGs在整体概述图、氨基酸代谢、信号转导、运输和分解代谢有富集。封盖后3 d预蛹以及封盖后6和9 d蛹受低温胁迫后共有的DEG 3-磷酸肌醇依赖性蛋白激酶1基因PDK1上调表达，同时富集在自噬-动物、mTOR和FoxO信号通路。低温胁迫后，封盖后3 d预蛹的胰岛素受体底物1-B基因IRS1-B和Kruppel同源物1基因Kr-h1上调表达，而核激素受体FTZ-F1基因Ftz-F1与蜕皮启动激素基因Eth显著下调表达，说明封盖后3 d预蛹响应低温细胞自噬和蜕皮受到抑制程度更大。在低温胁迫后封盖后6 d蛹中与昆虫角质层着色与免疫相关的酪氨酸羟化酶基因TyHyd下调表达；与对照组相比，在低温胁迫后封盖后9 d蛹中DEGs最少，说明在3个蛹期中，其受低温的影响较小。【结论】本研究测定了西方蜜蜂3个不同发育阶段蛹响应低温的DEGs，结果显示大部分DEGs为阶段特有，说明西方蜜蜂不同发育阶段响应低温的机制不同。一些共有DEGs以及阶段特有DEGs的功能研究和其作用机制是进一步研究蜜蜂对低温响应机制的重点内容，为探究蜜蜂蛹期响应低温胁迫的分子机制提供了基础数据。

Analysis of differentially expressed genes in response to low temperature stress in Apis mellifera worker pupae