基于RNA-seq技术的不同发育期中华大蟾蜍耳后腺和腹皮转录组比较分析

摘要

通过对中华大蟾蜍(Bufo bufo gargarizans)5个不同发育时期的耳后腺和腹皮进行转录组测序分析，挖掘参与毒素合成的关键因子，为深入研究毒素合成调控机理提供理论依据。分别取3月龄、 6月龄、 12月龄、 24月龄、 36月龄等5个发育时期的中华大蟾蜍耳后腺和腹皮进行转录组测序，利用GO和KEGG对差异表达基因功能进行分析，利用加权共表达网络分析筛选关键基因。将处于相同发育时期腹皮与耳后腺的转录组数据进行比较，分别筛选出3 053、 3 026、 1 516、 1 028和2 061个差异表达基因。由GO分类统计结果可知，差异表达基因在生物过程类别主要集中于细胞过程、单一生物体过程、代谢过程和生物调节；在分子功能类别主要富集于结合、催化活性等方面；在细胞组分类别主要聚集于细胞、细胞组分和细胞器。KEGG通路分析结果表明，差异表达基因在蟾蜍发育过程中参与多条通路，如酪氨酸代谢、过氧化物酶体和初级胆汁酸代谢等。将5个发育时期的共同差异表达基因进行比较，筛选出与类固醇代谢、初级胆汁酸代谢和过氧化物酶体相关的21个基因，利用加权共表达网络分析筛选出4个模块的6个关键基因。在不同发育时期中华大蟾蜍的腹皮与耳后腺差异表达基因中，大量代谢途径基因发生了富集，这些基因可能在调控耳后腺毒素合成中具有重要作用，为解析蟾蜍产毒分子机制奠定了基础。

关键词

中华大蟾蜍 / 产毒基因 / 转录组分析 / 差异表达基因