不同生长时期梅花鹿鹿茸转录组分析

doi:10.11843/j.issn.0366-6964.2017.02.006

摘要/Abstract

摘要：

本研究旨在探究不同生长时期梅花鹿鹿茸的转录组差异，丰富梅花鹿鹿茸转录组信息。选取5个生长时期（10、20、28、44、66 d）的梅花鹿鹿茸组织作为试验样品，利用Illumina HiSeq^TM 2500高通量测序，并用测序评估、基因注释等生物信息学方法进行分析。结果，经过测序、转录本拼接获得了375 657 条contigs，平均长度为688 bp，329 946个unigenes，平均长度为534 bp。通过比对Nr、Nt、Pfam、KOG/COG、Swiss-prot、KEGG、GO 等7大数据库，90.07% unigenes得到了成功注释。其中39 674个（12.02%）基因成功注释到GO数据库，在level 2水平上共分为41个类别；17 732个（5.37%）基因成功注释到将 KOG 的26个类别中。对5个生长时期的鹿茸转录本文库进行比较分析，共发现了509个差异基因，其中407个差异表达基因成功注释到GO数据库中，主要为信号转导、氧化还原、转录调节、蛋白酶降解等生物学过程。其中转录调节基因PER1、EGR1的表达随着鹿茸的生长而增加，而GAS1则相反，随着鹿茸的生长而降低，推测PER1、EGR1、GAS1等转录因子在鹿茸生长过程中可能起着重要的调节作用。本研究利用高通量测序技术对不同生长时期的梅花鹿鹿茸组织进行了转录组测序和分析，筛选到了梅花鹿鹿茸在不同生长时期下的差异表达基因，并获得了差异表达基因的功能。

关键词: 梅花鹿, 鹿茸, 转录组, 差异表达基因

Abstract:

The objective of this study was to find the transcriptome differences of sika deer antlers at different stages, and enrich the sika deer antler transcriptome data. The sika deer antlers aged 10, 20, 28, 44 and 66 d were selected as experimental material，then the transcriptome libraries of antlers were constructed by using Illumina HiSeq^TM 2500 platform, and the transcriptome data were analyzed by bioinformatics methods, such as sequencing assess and gene annotation. The results showed that a total of 375 657 contigs with an average length of 688 bp were generated, from which 329 946 unigenes with an average length of 534 bp were defined. 90.07% of the 297 198 unigenes were annotated in public databases of Nr，Nt，Pfam，KOG/COG，Swiss-prot，KEGG and GO. Through annotation classification, a total of 39 674(12.02%) genes were divided into 41 GO function categories at level 2, and 17 732 (5.37%) genes were annotated in 26 KOG categories. Comparative analysis of differentially expressed genes in sika deer antler at 5 different stages, and 509 differentially expressed genes were screened, from which 407 genes were annotated in GO database including signal transduction, oxidation-reduction process, regulation of transcription, proteolysis, and so on. And the expression of PER1 and EGR1 (transcription regulation gene) increased with the growth of antlers, on the contrary, the expression of GAS1 decreased. It indicates that PER1, EGR1, GAS1 may have important role in regulation of antler growth. The transcriptome study of the antler tissue in the different growth period revealed the number of differentially expressed genes by using high-throughput sequencing technology, obtained the function of differentially expressed genes.

Key words: sika deer, antler, transcriptome, differentially expressed genes

中图分类号:

S825
S813.1

张然然，刘华淼，王磊，邢秀梅，苏伟林，高兵. 不同生长时期梅花鹿鹿茸转录组分析[J]. 畜牧兽医学报, 2017, 48(2): 235-242.

ZHANG Ran-ran, LIU Hua-miao, WANG Lei, XING Xiu-mei, SU Wei-lin, GAO Bing. Transcriptome Analysis of Sika Deer Antlers at Different Stages[J]. ACTA VETERINARIA ET ZOOTECHNICA SINICA, 2017, 48(2): 235-242.

参考文献

［1］TANG Y, JEON B T, WANG Y, et al. First evidence that sika deer (Cervus nippon) velvet antler extract suppresses migration of human prostate cancer cells ［J］. Korean J Food Sci Anim Resour, 2015, 35(4): 507-514.
［2］HU W, QI L, TIAN Y H, et al. Studies on the purification of polypeptide from sika antler plate and activities of antitumor ［J］. BMC Complement Altern Med, 2015, 15: 328-335.
［3］JANG S, PARK E D, SUH H J, et al. Enhancement of exercise endurance capacity by fermented deer antler in BALB/c mice ［J］. Biosci Biotechnol Biochem, 2014, 78(10):1716-1722.
［4］ZHANG L Z, XIN J L, ZHANG X P, et al. The anti-osteoporotic effect of velvet antler polypeptides from Cervus elaphus Linnaeus in ovariectomized rats ［J］. J Ethnopharmacol, 2013, 150(1):181-186.
［5］LEE S H, YANG H W, DING Y, et al. Anti-inflammatory effects of enzymatic hydrolysates of velvet antler in RAW 264.7 cells in vitro and zebrafish model ［J］. EXCLI J, 2015, 14:1122-1132.
［6］牛维, 孙志涛, 曹学伟, 等. 单味药鹿茸调控大鼠骨关节炎软骨组织 Smad2, 3 表达的研究［J］. 中国中西医结合杂志, 2014, 34(2): 209-213.
NIU W, SUN Z T, CAO X W, et al. Regulation of single herb pilose antler on the expression of Smad2 and Smad3 in the cartilage of OA rats: an experimental research ［J］. Chinese Journal of Integrated Traditional and Western Medicine, 2014, 34(2):209-213.(in Chinese)
［7］TANG Y J, JEON B T, WANG Y M, et al. First evaluation of the biologically active substances and antioxidant potential of regrowth velvet antler by means of multiple biochemical assays ［J］. J Chem, 2015, 2015(2): 1-7.
［8］LI N, ZHANG M, DRUMMEN G P, et al. Sika deer antler collagen type I-accelerated osteogenesis in bone marrow Mesenchymal stem cells via the Smad pathway ［J］. Evid Based Complement Alternat Med, 2016, 2016(2):1-13.
［9］CHEN J, YANG Y, ABBASI S, et al. The effects of elk velvet antler dietary supplementation on physical growth and bone development in growing rats ［J］. Evid Based Compl Alt, 2015, 2015: DOI: 10.1155/2015/819520.
［10］HUO Y S, HUO H, ZHANG J. The contribution of deer velvet antler research to the modern biological medicine ［J］. Chin J Integr Med, 2014, 20(10):723-728.
［11］WANG W, WANG Y, ZHANG Q, et al. Global characterization of Artemisia annua glandular trichome transcriptome using 454 pyrosequencing ［J］. BMC Genomics, 2009, 10:465-475.
［12］SUN C, LI Y, WU Q, et al. De novo sequencing and analysis of the American ginseng root transcriptome using a GS FLX Titanium platform to discover putative genes involved in ginsenoside biosynthesis［J］. BMC Genomics, 2010, 11:262-274.
［13］ZENG S, XIAO G, GUO J, et al. Development of a EST dataset and characterization of EST-SSRs in a traditional Chinese medicinal plant, Epimedium sagittatum (Sieb. Et Zucc.) Maxim ［J］. BMC Genomics, 2010, 11:94-105.
［14］YUAN Y, SONG L, LI M, et al. Genetic variation and metabolic pathway intricacy govern the active compound content and quality of the Chinese medicinal plant Lonicera japonica thunb ［J］. BMC Genomics, 2012, 13:195-212.
［15］YAO B, ZHAO Y, WANG Q, et al. De novo characterization of the antler tip of Chinese Sika deer transcriptome and analysis of gene expression related to rapid growth ［J］. Mol Cell Biochem, 2012, 364(1-2):93-100.
［16］ZHAO Y, YAO B J, ZHANG M, et al. Comparative analysis of differentially expressed genes in Sika deer antler at different stages ［J］. Mol Biol Rep, 2013, 40(2):1665-1676.
［17］石睿. 马鹿鹿茸软骨组织转录组SNP与DIP分析［D］. 哈尔滨：东北林业大学, 2012.
SHI R. Analysis of SNP and DIP from transcriptional expression profiling of cartilage tissue in antler of the red deer (Cervus elaphus) ［D］. Harbin: Northeast Forestry University, 2012. (in Chinese)
［18］GRABHERR M G, HAAS B J, YASSOUR M, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome ［J］. Nat Biotechnol, 2011, 29(7):644-652.
［19］LI B, DEWEY C N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome ［J］. BMC Bioinformatics, 2011, 12:323-339.
［20］ZENG T, ZHANG L, LI J, et al. De novo assembly and characterization of Muscovy duck liver transcriptome and analysis of differentially regulated genes in response to heat stress ［J］. Cell Stress Chaperon, 2015, 20(3):483-493.
［21］HINOI E, UESHIMA T, HOJO H, et al. Up-regulation of per mRNA expression by parathyroid hormone through a protein kinase A-CREB-dependent mechanism in chondrocytes ［J］. J Biol Chem, 2006, 281(33):23632-23642.
［22］TAKARADA T, KODAMA A, HOTTA S, et al. Clock genes influence gene expression in growth plate and endochondral ossification in mice ［J］. J Biol Chem, 2012, 287(43):36081-36095.
［23］MCELDERRY J D, ZHAO G, KHMALADZE A, et al. Tracking circadian rhythms of bone mineral deposition in murine calvarial organ cultures［J］. J Bone Miner Res, 2013, 28(8):1846-1854.
［24］ADAMSON E D, MERCOLA D. EGR1 transcription factor: multiple roles in prostate tumor cell growth and survival ［J］. Tumor Biol, 2002, 23(2):93-102.
［25］DEL SAL G, RUARO M E, PHILIPSON L, et al. The growth arrest-specific gene, GAS1, is involved in growth suppression ［J］. Cell, 1992, 70(4):595-607.
［26］SUI Z, YUAN H, LIANG Z, et al. An activity-maintaining sequential protein extraction method for bioactive assay and proteome analysis of velvet antlers ［J］. Talanta, 2013, 107:189-194.

[1]	王璟, 滑留帅, 陈俊峰, 张家庆, 任巧玲, 白红杰, 郭红霞, 徐照学, 邢宝松, 白献晓, 曹海. 去势对淮南公猪背最长肌转录组的影响[J]. 畜牧兽医学报, 2019, 50(9): 1746-1758.
[2]	荣敏, 张然然, 涂剑峰, 王天娇, 孙丽敏, 杨颖, 徐佳萍, 邢秀梅*. 美国短毛黑水貂快速生长过程中转录组差异表达分析[J]. 畜牧兽医学报, 2019, 50(9): 1787-1801.
[3]	赵学亮, 王姝懿, 孙柯, 苏倩, 王文龙, 刘春霞. 捻转血矛线虫阿苯达唑敏感株和耐药株比较转录组学分析[J]. 畜牧兽医学报, 2019, 50(9): 1940-1944.
[4]	施会彬, 王玉琴, 张小辉, 吴至博, 杨芳. 基于转录组测序数据的山羊下丘脑组织繁殖相关信号通路及候选基因的筛选[J]. 畜牧兽医学报, 2019, 50(7): 1380-1395.
[5]	施会彬, 王玉琴, 张小辉, 吴至博, 杨芳. 基于转录组测序数据的山羊下丘脑组织繁殖相关信号通路及候选基因的筛选[J]. 畜牧兽医学报, 2019, 50(7): 1380-1395.
[6]	张石磊, 包佳鹭, 史万玉, 钟秀会. 母鼠孕期和哺乳期连续染毒双酚A对子代雄鼠睾丸发育的影响[J]. 畜牧兽医学报, 2019, 50(5): 993-1005.
[7]	张石磊, 包佳鹭, 史万玉, 钟秀会. 母鼠孕期和哺乳期连续染毒双酚A对子代雄鼠睾丸发育的影响[J]. 畜牧兽医学报, 2019, 50(5): 993-1005.
[8]	刘丽娟, 董然然, 王开功, 程振涛, 文明, 温贵兰, 李晨, 杨琦, 周碧君. 敲除小RNA gcvB后沙门菌的转录组分析[J]. 畜牧兽医学报, 2019, 50(4): 840-850.
[9]	刘丽娟, 董然然, 王开功, 程振涛, 文明, 温贵兰, 李晨, 杨琦, 周碧君. 敲除小RNA gcvB后沙门菌的转录组分析[J]. 畜牧兽医学报, 2019, 50(4): 840-850.
[10]	李欣, 李小俊, 陈晓丽, 赵毅强, 王栋. 转录组数据分析与功能基因挖掘[J]. 畜牧兽医学报, 2019, 50(3): 474-484.
[11]	阿果约达, 熊显荣, 王艳, 杨显英, 韩杰, 王斌, 李键. 犏牛雄性不育相关lncRNA的鉴定与分析[J]. 畜牧兽医学报, 2019, 50(3): 551-561.
[12]	赵学亮, 王姝懿, 呼和巴特尔, 孙柯, 苏倩, 吕旭, 王文龙, 刘春霞. 捻转血矛线虫阿苯达唑耐药株给药前后比较转录组学分析[J]. 畜牧兽医学报, 2019, 50(3): 637-644.
[13]	阿果约达, 熊显荣, 王艳, 杨显英, 韩杰, 王斌, 李键. 犏牛雄性不育相关lncRNA的鉴定与分析[J]. 畜牧兽医学报, 2019, 50(3): 551-561.
[14]	李欣, 李小俊, 陈晓丽, 赵毅强, 王栋. 转录组数据分析与功能基因挖掘[J]. 畜牧兽医学报, 2019, 50(3): 474-484.
[15]	赵学亮, 王姝懿, 呼和巴特尔, 孙柯, 苏倩, 吕旭, 王文龙, 刘春霞. 捻转血矛线虫阿苯达唑耐药株给药前后比较转录组学分析[J]. 畜牧兽医学报, 2019, 50(3): 637-644.