蒙古马高负荷运动训练前后转录组差异表达分析

doi:10.11843/j.issn.0366-6964.2017.06.005

畜牧兽医学报 ›› 2017, Vol. 48 ›› Issue (6): 1007-1016.doi: 10.11843/j.issn.0366-6964.2017.06.005

蒙古马高负荷运动训练前后转录组差异表达分析

赵启南, 芒来*, 白东义*, 赵一萍, 李蓓, 乌尼尔夫, 乌英嘎, 苏日嘎, 敖乳嘎

内蒙古农业大学动物科学学院内蒙古自治区蒙古马遗传资源保护及马产业工程实验室, 呼和浩特 010018

收稿日期:2016-11-10 出版日期:2017-06-23 发布日期:2017-06-15
通讯作者: 芒来,博士,教授,主要从事马属动物分子与数量遗传学研究,E-mail:dmanglai@163.com;白东义,博士,讲师,主要从事马属动物遗传育种与繁殖研究,E-mail:baidongyi1983@163.com
作者简介:赵启南(1986-),男,内蒙古呼和浩特人,博士生,主要从事马属动物运动学遗传基础研究,E-mail:s.pippen.33@163.com
基金资助:
国家自然科学基金（31360538；31472070）；内蒙古自治区科技厅重点实验室建设项目（20130902）；内蒙古自治区高等学校科学研究一般项目（NJZY0701）

Mongolian Horses Transcriptome Differential Expression Analysis before and after a High Load Exercise Training

ZHAO Qi-nan, MANG Lai*, BAI Dong-yi*, ZHAO Yi-ping, LI Bei, Unerhu, WU Ying-ga, SU Ri-ga, AO Ru-ga

Inner Mongolia Mongolian Horse Genetic Resources Protection and Industrial Engineering Laboratory, College of Animal Science, Inner Mongolia Agricultural University, Hohhot 010018, China

Received:2016-11-10 Online:2017-06-23 Published:2017-06-15

摘要/Abstract

摘要：

旨在对蒙古马高负荷运动训练前后转录组差异表达进行分析。本试验对6匹进行了4个月高负荷运动训练的蒙古马进行研究，分别在训练前后两个时期采集肌肉样品，利用二代测序技术对肌肉样品建立转录组文库，对发生显著变化的差异基因进行GO功能富集和KEGG Pathway分析，寻找训练前后发生显著变化的代谢通路及相关候选基因。结果表明，在转录组筛选到的1 102个显著差异表达基因中，有299个在训练后发生下调，803个发生上调。这些差异基因被注释到3 398条GO term上。其中与运动学相关的生物学过程有：肌肉结构发育、心血管系统发育、循环系统发育、肌肉组织发育、肌肉器官发育、肌肉收缩和肌肉细胞分化等。对差异表达基因进行KEGG pathway分析，结果表明，差异基因被富集到如扩张型心肌病、肥厚型心肌病、心肌收缩、钙信号通路、肌动蛋白骨架调节等与运动相关的193条通路中。本研究结果为探究蒙古马强耐力的基本分子机理提供理论基础，同时也为理解蒙古马运动学特性提供新的思路。

关键词: 蒙古马, 高负荷运动, 转录组

Abstract:

The aim of this study was to analyze the Mongolian horse transcriptome differential expression before and after a high load exercise training. This research selected 6 Mongolian horses with high load training for 4 months. Muscle samples were collected before and after the training respectively. The transcriptome library was established using the second generation sequencing technology. We enriched the differentially expressed genes to the GO and KEGG pathway, and screened the significantly different metabolic pathways and related genes. 1 102 differentially expressed genes(DEGs) in transcriptome library were screened out, of which 299 were down-regulated and 803 were up-regulated after high load training, respectively. These DEGs were annotated to 3 398 GO terms, including muscle structure development, cardiovascular system development, circulatory system development, muscle tissue development, muscular organ development, muscle contraction, muscle cell differentiation, etc. Also the DEGs were annotated to 193 KEGG pathways including some athletics related pathways: dilated cardiomyopathy (DCM), hypertrophic cardiomyopathy, cardiac muscle contraction, calcium signaling, the actin cytoskeleton adjustment, etc. Our findings provide a theory support for further investigating the molecular mechanism of Mongolian horse strong stamina and a new research method for Mongolian horses athletic characteristics.

Key words: Mongolian horse, high load exercise, transcriptome

中图分类号:

S821.2

赵启南, 芒来, 白东义, 赵一萍, 李蓓, 乌尼尔夫, 乌英嘎, 苏日嘎, 敖乳嘎. 蒙古马高负荷运动训练前后转录组差异表达分析[J]. 畜牧兽医学报, 2017, 48(6): 1007-1016.

ZHAO Qi-nan, MANG Lai, BAI Dong-yi, ZHAO Yi-ping, LI Bei, Unerhu, WU Ying-ga, SU Ri-ga, AO Ru-ga. Mongolian Horses Transcriptome Differential Expression Analysis before and after a High Load Exercise Training[J]. ACTA VETERINARIA ET ZOOTECHNICA SINICA, 2017, 48(6): 1007-1016.

参考文献

[1] 吴克亮, 吴常信. 马科学研究动态和马业发展[J]. 畜牧兽医学报, 2005, 36(4):412-416.
WU K L, WU C X. Advanced research on equine science and horse industry[J]. Acta Veterinaria et Zootechnica Sinica, 2005, 36(4):412-416. (in Chinese)
[2] 芒来. 蒙古马的历史与展望[J]. 马业杂志, 2002(4):14-17.
MANG L. History and outlook of the Mongolian horse[J]. Journal of Horse Industry, 2002(4):14-17. (in Chinese)
[3] MOSHER D S, QUIGNON P, BUSTAMANTE C D, et al. A mutation in the myostatin gene increases muscle mass and enhances racing performance in heterozygote dogs[J]. PLoS Genet, 2007, 3(5):e79.
[4] ZHANG X M, ZHAO L, LARSON-RABIN Z, et al. De novo sequencing and characterization of the floral transcriptome of Dendrocalamus latiflorus (Poaceae:Bambusoideae)[J]. PLoS One, 2012, 7(8):e42082.
[5] VARSHNEY R K, SIGMUND R, BÖRNER A, et al. Interspecific transferability and comparative mapping of barley EST-SSR markers in wheat, rye and rice[J]. Plant Sci, 2005, 168(1):195-202.
[6] CÁNOVAS A, RINCON G, ISLAS-TREJO A, et al. SNP discovery in the bovine milk transcriptome using RNA-Seq technology[J]. Mamm Genome, 2010, 21(11-12):592-598.
[7] LI X Y, SUN H Y, PEI J B, et al. De novo sequencing and comparative analysis of the blueberry transcriptome to discover putative genes related to antioxidants[J]. Gene, 2012, 511(1):54-61.
[8] GOECKS J, NEKRUTENKO A, TAYLOR J, et al. GALAXY:a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences[J]. Genome Biol, 2010, 11(8):R86.
[9] JAILLON O, AURY J M, NOEL B, et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla[J]. Nature, 2007, 449(7161):463-467.
[10] CHEN Q M, WANG H, ZENG Y Q, et al. Developmental changes and effect on intramuscular fat content of H-FABP and A-FABP mRNA expression in pigs[J]. J Appl Genet, 2013, 54(1):119-123.
[11] BONNET A, BEVILACQUA C, BENNE F, et al. Transcriptome profiling of sheep granulosa cells and oocytes during early follicular development obtained by Laser Capture Microdissection[J]. BMC Genomics, 2011, 12:417.
[12] MILLER J R, KOREN S, SUTTON G. Assembly algorithms for next generation sequencing data[J]. Genomics, 2010, 95(6):315-327.
[13] http://asia.ensembl.org/Equus_caballus/Info/Index.
[14] ERLICH Y, MITRA P P, DELABASTIDE M, et al. Alta-cyclic:a self-optimizing base caller for next-generation sequencing[J]. Nat Methods, 2008, 5(8):679-682.
[15] ESTEVE-CODINA A, KOFLER R, PALMIERI N, et al. Exploring the gonad transcriptome of two extreme male pigs with RNA-seq[J]. BMC Genomics, 2011, 12:552.
[16] KANG X L, LIU G, LIU Y F, et al. Transcriptome profile at different physiological stages reveals potential mode for curly fleece in Chinese Tan sheep[J]. PLoS One, 2013, 8(8):e71763.
[17] LOBJOIS V, LIAUBET L, SANCRISTOBAL M, et al. A muscle transcriptome analysis identifies positional candidate genes for a complex trait in pig[J]. Anim Genet, 2008, 39(2):147-162.
[18] ANDERS S, HUBER W. Differential expression analysis for sequence count data[J]. Genome Biol, 2010, 11(10):R106.
[19] NIELSEN M, HANSEN J H, HEDEGAARD J, et al. microRNA identity and abundance in porcine skeletal muscles determined by deep sequencing[J]. Anim Genet, 2010, 41(2):159-168.
[20] KEGG:Kyoto encyclopedia of genes and genomes[EB/OL]. http://www.genome.jp/kegg/.
[21] DUTTA S, KUMAWAT G, SINGH B P, et al. Development of genic-SSR markers by deep transcriptome sequencing in pigeonpea[Cajanus cajan (L.) Millspaugh] [J]. BMC Plant Biol, 2011, 11:17.
[22] ZORZATO F, FUJII J, OTSU K, et al. Molecular cloning of cDNA encoding human and rabbit forms of the Ca²⁺ release channel (ryanodine receptor) of skeletal muscle sarcoplasmic reticulum[J]. J Biol Chem, 1990, 265(4):2244-2256.
[23] GAJDOSIK R L. Passive extensibility of skeletal muscle:review of the literature with clinical implications[J]. Clin Biomech, 2001, 16(2):87-101.
[24] SCHRÖDER W, KLOSTERMANN A, DISTL O. Candidate genes for physical performance in the horse[J]. Vet J, 2011, 190(1):39-48.
[25] MANICKAM R, PENA R N, WHITELAW C B A. Mammary gland differentiation inversely correlates with GDF-8 expression[J]. Mol Reprod Dev, 2008, 75(12):1783-1788.
[26] BRYANT N J, GOVERS R, JAMES D E. Regulated transported transport of the glucose transporter GLUT4[J]. Nat Rev Mol Cell Biol, 2002, 3(4):267-277.
[27] 张万坡, 田文霞, 李建丽, 等. 禽流感病毒NS1蛋白与鸡LY6E蛋白间的相互作用[J]. 中国兽医学报, 2008, 28(6):627-630, 645.
ZHANG W P, TIAN W X, LI J L, et al. Interaction between nonstructural (NS1) protein of avian influenza virus and LY6 E proteins of chicken[J]. Chinese Journal of Veterinary Science, 2008, 28(6):627-630, 645. (in Chinese)
[28] 吴帅, 简斌. BCL9与肿瘤及肿瘤靶点治疗研究[J]. 生命科学研究, 2013, 17(1):86-89.
WU S, JIAN B. Research on BCL9 with tumor and tumor target treatment[J]. Life Science Research, 2013, 17(1):86-89. (in Chinese)
[29] 邹剑, 周后凤, 先志伟, 等. 激活PPARα表达对AngⅡ诱导的心肌细胞肥大及NFATc4与p65-NFκB相互作用的影响[J]. 中国病理生理杂志, 2014, 30(6):1017-1022.
ZOU J, ZHOU H F, XIAN Z W, et al. Effects of PPARα activation on AngⅡ-induced cardiomyocyte hypertro-phy and interaction of NFATc4 with p65-NFκB[J]. Chinese Journal of Pathophysiology, 2014, 30(6):1017-1022. (in Chinese)
[30] 王永珍. 峨眉树蛙和斑腿树蛙MHCⅠ类基因的多态性及其进化[D]. 武汉:华中师范大学, 2012.
WANG Y Z. Polymorphism and evolution of MHC class Ⅰ genes of Rhscophorus omeimontis and Rhscophorus megacephalus[D]. Wuhan:Central China Normal University, 2012. (in Chinese)
[31] 孟祥阁. 抗阻、耐力及混合运动对C57小鼠股四头肌MHC各亚型基因表达的影响[D]. 上海:华东师范大学, 2010.
MENG X G. The mechanisms and the recovery process of exercise-induced fatigue[D]. Shanghai:East China Normal University, 2010. (in Chinese)
[32] 王丛阳, 王辉. 肌球蛋白轻链的调节及其对肿瘤细胞增殖、转移的影响[J]. 医学综述, 2010, 16(21):3251-3254.
WANG C Y, WANG H. Regulation of myosin light chain and its effect on tumor proliferation and metastas[J]. Medical Recapitulate, 2010, 16(21):3251-3254. (in Chinese)
[33] 李航, 王彬彬, 高秉仁, 等. 一个新的ACTC1基因5'端剪切位点突变可能在先天性心脏病室间隔缺损发病中起重要作用[J]. 中国组织工程研究, 2012, 16(11):2087-2090.
LI H, WANG B B, GAO B R, et al. A novel 5' splice site mutation in ACTC1 gene may play an important role in congenital ventricular septal defect[J]. Chinese Journal of Tissue Engineering Research, 2012, 16(11):2087-2090. (in Chinese)

蒙古马高负荷运动训练前后转录组差异表达分析

Mongolian Horses Transcriptome Differential Expression Analysis before and after a High Load Exercise Training

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[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.