浙江农业学报 ›› 2017, Vol. 29 ›› Issue (9): 1575-1580.DOI: 10.3969/j.issn.1004-1524.2017.09.22

• 食品科学 • 上一篇    下一篇

中华蜜蜂幼虫肠道的高表达基因分析

解彦玲, 王鸿权, 张璐, 侯志贤, 刀晨, 江亮亮, 熊翠玲, 郑燕珍, 徐细建, 黄枳腱, 郭睿*, 陈大福   

  1. 福建农林大学 蜂学学院,福建 福州 350002
  • 收稿日期:2017-03-20 出版日期:2017-09-20 发布日期:2017-09-27
  • 通讯作者: 郭睿,E-mail:fafu_ruiguo@126.com
  • 作者简介:解彦玲(1996—),女,浙江湖州人,主要从事蜂学研究。E-mail:1103430452@qq.com
  • 基金资助:
    现代农业产业技术体系建设专项基金(CARS-45-KXJ7); 福建农林大学科技发展基金(KF2015123); 福建省大学生创新创业项目(201610389053)

Analysis of highly expressed genes of larval gut of Apis cerana cerana

XIE Yanling, WANG Hongquan, ZHANG Lu, HOU Zhixian, DAO Chen, JIANG Liangliang, XIONG Cuiling, ZHENG Yanzhen, XU Xijian, HUANG Zhijian, GUO Rui*, CHEN Dafu   

  1. College of Bee Science, Fujian Agriculture and Forestry University, Fuzhou 350002, China
  • Received:2017-03-20 Online:2017-09-20 Published:2017-09-27

摘要: 中华蜜蜂(中蜂)是中国特有的本土蜜蜂资源,也是养蜂生产中的常用蜂种之一,具有重要的生态和经济价值。本研究利用RNA-seq技术对中蜂5和6日龄幼虫肠道进行深度测序,得到的有效读段(clean reads)数分别为29290630与26636038,Q20分别为98.45%和98.38%。GO富集分析结果显示,Ac5中的高表达基因(HEGs)富集在44个GO term,其中基因富集数最多的为代谢进程(753 unigenes),其次是催化活性(679 unigenes)与结合(657 unigenes);Ac6中的HEGs富集在44个GO term,其中基因富集数最多的为代谢进程(792 unigenes),其次是结合(760 unigenes)和催化活性(744 unigenes)。KEGG代谢通路(pathway)富集分析结果显示,Ac5的HEGs富集在129个pathway,其中基因富集数最多的是核糖体(103 unigenes)、内质网蛋白质加工(68 unigenes)和碳代谢(67 unigenes);Ac6的HEGs富集在129个代谢通路,其中基因富集数最多的是核糖体(103 unigenes)、氧化磷酸化(88 unigenes)及内质网蛋白质加工(68 unigenes)。进一步对HEGs进行Venn分析,结果显示,二者共有的HEGs为2 015个,特有的HEGs分别为219和491个。研究结果不仅可提供中蜂幼虫肠道发育过程中的基因表达谱数据,也为分子水平上深入研究中蜂幼虫肠道发育提供重要信息。

关键词: RNA-seq, 中华蜜蜂, 幼虫肠道, 高表达基因

Abstract: Apis cerana cerana is not only a specific honeybee species resource in China, but also a frequently-used honeybee species in apiculture. In the present study, the 5-(Ac5) and 6-day-old (Ac6) larval guts of A. c. cerana were sequenced using RNA-seq. After filtration, 29290630 and 26636038 clean reads with a Q20 of 98.45% and 98.38% were obtained in Ac5 and Ac6. GO enrichment analysis results showed that the highly expressed genes (HEGs) in Ac5 were enriched in 44 GO terms, among them metabolic process (753 unigenes), catalytic activity (679 unigenes) and binding (657 unigenes) were the largest groups; the HEGs in Ac6 were enriched in 44 GO terms, and the mostly enriched one was metabolic process (792 unigenes), followed by binding (760 unigenes) and catalytic activity (744 unigenes). KEGG enrichment analysis results suggested that the HEGs of Ac5 were enriched in 129 pathways, among them ribosome (103 unigenes) was the largest group, followed by protein processing in endoplasmic reticulum (68 unigenes) and carbon metabolism (67 unigenes); the HEGs of Ac6 were enriched in 129 pathways, the mostly enriched ones were ribosome (103 unigenes), oxidative phosphorylation (88 unigenes) as well as protein processing in endoplasmic reticulum (68 unigenes). Our findings can not only offer the gene expression profiles during the developmental process of the larval gut of A. c. cerana at transcriptome level, but also provide key information for further investigation of the larval gut's development.

Key words: RNA-seq, Apis cerana cerana, larval gut, highly expressed genes

中图分类号: