刺猎蝽属种间转录组比较分析

作者：邓美菁; 赵萍^*

来源：昆虫学报, 2023, 66(09): 1221-1232.

DOI：10.16380/j.kcxb.2023.09.009

• 刺猎蝽属
• 高通量测序
• 转录组
• 差异表达基因

【目的】建立刺猎蝽属Sclomina 3种的转录组数据库，分析近缘种间转录组表达差异，探讨在转录组水平的种间趋异情况。【方法】采用Illumina HiSeqTM4000高通量测序平台进行刺猎蝽属齿缘刺猎蝽S.erinacea、兴仁刺猎蝽S.xingrensis和广西刺猎蝽S.guangxiensis 3种转录组测序；利用Nr, Swiss-Prot, COG/KOG和KEGG数据库进行基因功能注释；对种间差异表达基因(differentially expressed genes, DEGs)进行GO功能注释和KEGG通路富集分析。【结果】齿缘刺猎蝽、兴仁刺猎蝽和广西刺猎蝽转录组测序组装获得42 215个unigenes。21 117个基因在上述4个数据库中得到注释(各数据库注释基因数为：Nr, 20 522; Swiss-Prot, 15 550; COG/KOG, 13 969; KEGG,10 850)。兴仁刺猎蝽vs齿缘刺猎蝽转录组有3 390个DEGs (803个上调，2 587个下调),兴仁刺猎蝽vs广西刺猎蝽转录组有12 543个DEGs (6 639个上调，5 904个下调),齿缘刺猎蝽vs广西刺猎蝽转录组有3 580个DEGs(2 485个上调，1 095个下调)。GO功能注释结果显示，上述DEGs被注释到生物学过程、细胞组分和分子功能3个功能大类和功能亚类，兴仁刺猎蝽vs齿缘刺猎蝽转录组中DEGs分别注释到代谢过程、细胞和细胞部分等50个功能亚类中，下调基因多，上调基因少；兴仁刺猎蝽vs广西刺猎蝽转录组中DEGs被注释到生物学过程、细胞部分和细胞等54个功能亚类中；齿缘刺猎蝽vs广西刺猎蝽转录组中DEGs被分别注释到细胞部分、代谢过程等46个功能亚类中。6 254个DEGs被注释到KEGG代谢通路，涉及多个不同的代谢通路，包括代谢通路、药物代谢-细胞色素P450、核糖体等；兴仁刺猎蝽vs齿缘刺猎蝽转录组702个DEGs富集到164个KEGG通路；兴仁刺猎蝽vs广西刺猎蝽转录组中2 091个DEGs富集得到201个KEGG通路；齿缘刺猎蝽vs广西刺猎蝽转录组中865个DEGs富集得到129个KEGG通路；兴仁刺猎蝽vs广西刺猎蝽转录组DEGs富集到代谢通路、药物代谢-细胞色素P450和核糖体等KEGG通路的数量最多。【结论】本研究结果丰富了刺猎蝽属的分子生物学数据；在转录组水平上呈现了齿缘刺猎蝽、兴仁刺猎蝽和广西刺猎蝽的种间差异，探索了该属物种分化的实际情况。

• 单位
  南宁师范大学; 生命科学学院