昆虫学报 ›› 2023, Vol. 66 ›› Issue (1): 11-18.doi: 10.16380/j.kcxb.2023.01.002

• 研究论文 • 上一篇    下一篇

中华蜜蜂细胞色素P450基因及其全长转录本的鉴定及分析

蔡宗兵1,#, 王紫馨1,#, 吴鹰1, 王思懿1, 钱加珺1, 胡颖1, 张凯遥1顾小雨1, 徐细建3, 骆群3, 陈大福1,2,*, 郭睿1,2,*   

  1. (1. 福建农林大学动物科学学院(蜂学学院), 福州 350002; 2. 福建省蜂疗研究所, 福州 350002; 3. 江西省养蜂研究所, 南昌 330000)
  • 出版日期:2023-01-20 发布日期:2023-03-04

Identification and analysis of cytochrome P450 genes and their full-length transcripts in Apis cerana cerana (Hymenoptera: Apidae)#br#

CAI Zong-Bing1,#, WANG Zi-Xin1,#, WU Ying1, WANG Si-Yi1, QIAN Jia-Jun1, HU Ying1, ZHANG Kai-Yao1, GU Xiao-Yu1, XU Xi-Jian3, LUO Qun3, CHEN Da-Fu1,2,*, GUO Rui1,2,*   

  1.  (1. College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, Fuzhou 350002, China; 2. Apitherapy Research Institute of Fujian Province, Fuzhou 350002, China; 3. Apicultural Research Institute of Jiangxi Province, Nanchang 330000, China)
  • Online:2023-01-20 Published:2023-03-04

摘要: 【目的】利用纳米孔(nanopore)测序技术鉴定中华蜜蜂Apis cerana cerana工蜂幼虫肠道中细胞色素P450(cytochrome P450, CYP450)基因及其全长转录本,为后续功能研究提供参考信息和基础。【方法】通过Nanopore PromethION平台对中华蜜蜂工蜂4-6日龄幼虫肠道进行转录组测序。利用Guppy软件对原始读段(raw reads)进行质控以得到有效读段(clean reads)。通过识别两端引物鉴定全长转录本序列。使用BLAST工具将上述全长转录本的序列比对到Nr和GO数据库以鉴定CYP450基因及其全长转录本。采用Astalavista软件鉴定基因的可变剪接(alternative splicing, AS)事件。通过RTPCR验证不同类型AS事件的可靠性。【结果】在中华蜜蜂工蜂4-6日龄幼虫肠道中分别测得7 338 627, 7 003 419和7 434 233条原始读段,经质控得到的有效读段数分别为7 289 494, 6 959 880和7 387 756条。鉴定到的非冗余全长转录本总数为48 200条。共鉴定到47个CYP450基因和265条CYP450基因全长转录本。共鉴定到CYP450基因的90次AS事件,包括36次外显子跳跃事件、20次可变5′端剪接位点事件、17次内含子保留事件、9次可变3′端剪接位点事件及8次外显子互斥事件。RTPCR结果证实随机选取的3种AS事件类型真实可靠。【结论】鉴定了中华蜜蜂的CYP450基因及其全长转录本,补充了东方蜜蜂参考基因组的相关注释,并揭示中华蜜蜂CYP450基因可通过多种AS类型产生丰富的剪接体。

关键词: 中华蜜蜂, 幼虫, 第三代测序技术, 细胞色素P450, 全长转录本, 可变剪接

Abstract: 【Aim】To identify cytochrome P450 (CYP450) genes and their full-length transcripts in the larval guts of Apsi cerana cerana workers using nanopore sequencing technology and offer reference information and basis for further functional study.【Methods】The 4-6-day-old larval gut transcriptomes of A. cerana cerana workers were sequenced by Nanopore PromethION platform. Quality control of raw reads was performed using Guppy software to gain clean reads. Full-length transcript sequences were detected through recognizing primers at both ends. The sequences of the above-mentioned full-length transcripts were aligned to the Nr and GO databases with BLAST tool to identify CYP450 genes and their full-length transcripts. Astalavista software was used to identify alternative splicing (AS) events of genes. RT-PCR was used to validate the reliability of AS events of various types. 【Results】In the 4-6-day-old larval guts of A. cerana cerana workers, 7 338 627, 7 003 419 and 7 434 233 raw reads were obtained, respectively, and after quality control, 7 289 494, 6 959 880 and 7 387 756 clean reads were gained. The total number of identified non-redundant fulllength transcripts was 48 200. Forty-seven CYP450 genes and their 265 full-length transcripts were identified. A total of 90 AS events of CYP450 genes were identified, including 36 exon skipping events, 20 alternative 5′ splicing site events, 17 intron retention events, nine alternative 3′ splicing site events and eight mutually exclusive exon events. The RT-PCR result confirmed the authenticity of randomly selected three types of AS events. 【Conclusion】 The findings provide CYP450 genes and their full-length transcripts in A. cerana cerana, supplement the annotation of A. cerana reference genome, and reveal that CYP450 genes of A. cerana cerana can generate abundant isoforms via multiple AS types.

Key words: Apis cerana cerana, larvae, third-generation sequencing technique, cytochrome P450, full-length transcripts, alternative splicing