利用全基因组重测序分析鹿茸重量相关基因

doi:10.16288/j.yczz.17-116

遗传 ›› 2017, Vol. 39 ›› Issue (11): 1090-1101.doi: 10.16288/j.yczz.17-116

利用全基因组重测序分析鹿茸重量相关基因

胡鹏飞¹(),徐佳萍¹,艾成²,邵秀娟²(),王洪亮¹,董依萌³,崔学哲¹,杨福合¹,邢秀梅¹

1. 中国农业科学院特产研究所,特种经济动物分子生物学国家重点实验室,农业部特种经济动物遗传育种殖重点实验室,长春 130112
2. 中国农业科学院深圳农业基因组研究所,深圳 518120
3. 吉林农业大学,长春 130118

收稿日期:2017-03-30 修回日期:2017-06-11 出版日期:2017-11-20 发布日期:2017-12-20
作者简介:胡鹏飞,博士,副研究员,研究方向：梅花鹿产茸性能评价及分子调控机制。E-mail: pfhoo@hotmail.com|邢秀梅，博士，研究员，研究方向：特种经济动物种质资源保护与遗传育种。E-mail: xingxiumei2004@126.com
基金资助:
吉林省科技发展计划项目(20150309006YY)

Screening weight related genes of velvet antlers by whole genome re-sequencing

Hu Pengfei¹(),Xu Jiaping¹,Ai Cheng²,Shao Xiujuan²(),Wang Hongliang¹,Dong Yimeng³,Cui Xuezhe¹,Yang Fuhe¹,Xing Xiumei¹

1. State Key Laboratory for Molecular Biology of Special Economic Animals, Key Laboratory of Special Economic Animal Genetics and Breeding, Ministry of Agriculture, Institute of Special Economic Animal and Plant Sciences, Chinese Academy of Agricultural Sciences, Changchun 130112, China
2. Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
3. Jilin Agricultural University, Changchun 130118, China

Received:2017-03-30 Revised:2017-06-11 Online:2017-11-20 Published:2017-12-20
Supported by:
the Science and Technology Development Project of Jilin Province(20150309006YY)

摘要/Abstract

摘要：

茸角是鹿科动物特有的器官,具有重要的生物学意义。鹿茸生长是一个复杂的生物代谢过程,其重量与遗传因素有一定关联。本研究对饲养条件基本一致的5个梅花鹿(Cervus nippon)群体进行调查,获得高产和低产梅花鹿个体共100只,利用全基因组重测序分析这些个体与鹿茸重量相关的遗传变异。结果表明,共得到94个与鹿茸重量可能相关的遗传变异,其中有2个变异位点分别定位于OAS2和ALYREF/THOC4基因的外显子区,且ALYREF/THOC4基因在鹿茸中表达量很高。功能富集分析发现,这些遗传变异与鹿茸生长发育密切相关,可作为潜在的鹿茸重量相关遗传变异。本研究首次通过全基因组重测序直接筛选与鹿茸重量相关的遗传变异,并分析关联基因的生物学功能,对揭示鹿茸生长发育和鹿茸重量差异形成的遗传机制具有重要意义。

关键词: 全基因组重测序, 梅花鹿, 鹿茸重量, 全基因组关联分析

Abstract:

The velvet antler is a special organ that has important biological significance for deer, and its growth is a complicated biological metabolism process. Growing evidence suggests that genetics factors play essential roles in the weight of velvet antlers. In this study, we investigated five sika deer (Cervus nippon) populations under the same feeding condition, and screened genetic variations in the 100 samples (including 50 heavy and 50 light velvet antler weight samples) by whole genome re-sequencing. The results showed that 94 genetic variations were related to the velvet antler weight, among which two single nucleotide polymorphism (SNP) sites were located on the exon regions of OAS2 and ALYREF/THOC4, respectively. Furthermore, ALYREF/THOC4 is highly expressed in the velvet antler. The biological functions of these genetic variations were highly related to the growth and development of deer velvet antlers. Collectively, we screened genes related to the velvet antler weight in sika deer populations by whole genome re-sequencing and identified 94 sites as candidate genetic variations related to the velvet antler weight. We hope that it will contribute to further mechanistic studies of velvet antler development and weight variations.

Key words: whole genome re-sequencing, Cervus nippon, velvet antler weight, genome wide association study

胡鹏飞, 徐佳萍, 艾成, 邵秀娟, 王洪亮, 董依萌, 崔学哲, 杨福合, 邢秀梅. 利用全基因组重测序分析鹿茸重量相关基因[J]. 遗传, 2017, 39(11): 1090-1101.

Hu Pengfei, Xu Jiaping, Ai Cheng, Shao Xiujuan, Wang Hongliang, Dong Yimeng, Cui Xuezhe, Yang Fuhe, Xing Xiumei. Screening weight related genes of velvet antlers by whole genome re-sequencing[J]. Hereditas(Beijing), 2017, 39(11): 1090-1101.

参考文献

[1]	Lincoln GA. Appearance of antler pedicles in early foetal life in red deer. J Embryol Exp Morphol, 1973, 29(2): 431-437.
[2]	Hartwig H, Schrudde J. Experimentelle untersuchungen zur bildung der prim?ren stirnauswüchse beim reh (Capreolus capreolus L.). Z Jagdwiss, 1974, 20(1): 1-13.
[3]	Yue ZP, Deng XM, Feng HH. The mechanisms of deer antlers development and regeneration. J Econ Anim, 2005, 9(1): 46-49
[3]	岳占碰, 邓旭明, 冯海华. 鹿茸角发育与再生机理. 经济动物学报, 2005, 9(1): 46-49.
[4]	赵肯堂. 角的形态、性别、生长和意义. 生物学通报, 1958, (12): 8-11.
[5]	蒋志刚. 生物进化的原动力——性选择. 科学中国人, 2002, (6): 54-56.
[6]	Williams JD, Krueger WF, Harmel DH. Heritabilities for antler characteristics and body weight in yearling white- tailed deer. Heredity, 1994, 73(1): 78-83.
[7]	Pérez-González J, Carranza J, Torres-Porras J, Fernández- García JL. Low heterozygosity at microsatellite markers in Iberian Red Deer with small antlers. J Hered, 2010, 101(5): 553-561.
[8]	Price JS, Allen S, Faucheux C, Althnaian T, Mount JG. Deer antlers: a zoological curiosity or the key to understanding organ regeneration in mammals?. J Anat, 2005, 207(5): 603-618.
[9]	Li CY, Zhao HP, Liu Z, McMahon C. Deer antler-a novel model for studying organ regeneration in mammals. Int J Biochem Cell Biol, 2014, 56: 111-122.
[10]	Sui ZG, Zhang LH, Huo YS, Zhang YK. Bioactive components of velvet antlers and their pharmacological properties. J Pharm Biomed Anal, 2014, 87: 229-240.
[11]	Tseng SH, Sung CH, Chen LG, Lai YJ, Chang WS, Sung HC, Wang CC. Comparison of chemical compositions and osteoprotective effects of different sections of velvet antler. J Ethnopharmacol, 2014, 151(1): 352-360.
[12]	林仲凡. 有关鹿茸( 角)的生物学研究进展. 国外畜牧学(草食家畜), 1985, (4): 7-10.
[13]	Kruuk LEB, Slate J, Pemberton JM, Brotherstone S, Guinness F, Clutton-Brock T. Antler size in red deer: heritability and selection but no evolution. Evolution, 2002, 56(8): 1683-1695.
[14]	李和平. 梅花鹿优良品种(品系)产茸性能研究. 中国畜牧杂志, 2003, 39(2): 31-32.
[15]	Li HP, Zheng XT, Bing GL, Li S, Zhao M, Chen XJ. Specific character analysis on artificial breed and strain of deer in China,. Hereditas(Beijing) 1997, 19(S1): 76-78.
[15]	李和平, 郑兴涛, 邴国良, 李生, 赵蒙, 陈学进. 中国茸鹿人工培育品种(品系)种质特性分析. 遗传, 1997, 19(S1): 76-78.
[16]	Martin M. Cutadapt removes adapter sequences from high- throughput sequencing reads. EMBnet J, 2011, 17(1): 10-12.
[17]	Cox MP, Peterson DA, Biggs PJ. SolexaQA: At-a-glance quality assessment of Illumina second-generation sequencing data. BMC Bioinformatics, 2010, 11: 485.
[18]	Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics, 2009, 25(14): 1754-1760.
[19]	Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Proc GPD. The sequence alignment/map format and SAMtools. Bioinformatics, 2009, 25(16): 2078-2079.
[20]	Ye K, Schulz MH, Long Q, Apweiler R, Ning ZM. Pindel: a pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads. Bioinformatics, 2009, 25(21): 2865-2871.
[21]	Wang K, Li MY, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res, 2010, 38(16): e164.
[22]	Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, McVean G, Durbin R, 1000 Genomes Project Analysis Group. The variant call format and VCFtools. Bioinformatics, 2011, 27(15): 2156-2158.
[23]	Nei M, Kumar S. Molecular evolution and phylogenetics. New York: Oxford University Press, 2000: 385-385.
[24]	Jia BY, Ba HX, Wang GW, Yang Y, Cui XZ, Peng YH, Zheng JJ, Xing XM, Yang FH. Transcriptome analysis of sika deer in China. Mol Genet Genomics, 2016, 291(5): 1941-1953.
[25]	Luo WZ, Cheng DX, Chen SK, Wang LG, Li Y, Ma XJ, Song X, Liu X, Li W, Liang J, Yan H, Zhao KB, Wang CD, Wang LX, Zhang LC. Genome-wide association analysis of meat quality traits in a porcine Large White × Minzhu intercross population. Int J Biol Sci, 2012, 8(4): 580-595.
[26]	Hovnanian A, Rebouillat D, Mattei MG, Levy ER, Marié I, Monaco AP, Hovanessian AG. The human 2°,5°-oligoadenylate synthetase locus is composed of three distinct genes clustered on chromosome 12q24.2 encoding the 100-, 69-, and 40-kDa forms. Genomics, 1998, 52(3): 267-277.
[27]	Muravenko OV, Gizatullin RZ, Al-Amin AN, Protopopov AI, Kashuba VI, Zelenin AV, Zabarovsky ER. Human ALY/BEF gene Map position 17q25.3. Chromosome Res, 2000, 8(6): 562-562.
[28]	Wei X, Liu KY, Zhang YX, Feng Q, Wang LH, Zhao Y, Li DH, Zhao Q, Zhu XD, Zhu XF, Li WJ, Fan DL, Gao Y, Lu YQ, Zhang XM, Tang XM, Zhou CC, Zhu CR, Liu LF, Zhong RC, Tian QL, Wen ZR, Weng QJ, Han B, Huang XH, Zhang XR. Genetic discovery for oil production and quality in sesame. Nat Commun, 2015, 6: 8609.
[29]	Thamizhmani R, Vijayachari P. Association of dengue virus infection susceptibility with polymorphisms of 2°-5°- oligoadenylate synthetase genes: a case-control study. Braz J Infect Dis, 2014, 18(5): 548-550.
[30]	Du ZH, Bai XJ. Association analysis between SNPs of the growth hormone gene and antler production in spotted deers. Hereditas (Beijing), 2007, 29(3): 337-342.
[30]	杜智恒, 白秀娟. 梅花鹿生长激素基因单核苷酸多态与产茸量性状的相关性. 遗传, 2007, 29(3): 337-342.
[31]	Xiong JJ. Effects of hormone and hormone receptor gene mutations on the regulation of the antler production in Sika Deer[Dissertation]. Central China Agricultural University, 2012.
[31]	熊家军. 梅花鹿激素及其受体基因变异对产茸量的调控效应研究[Dissertation]. 华中农业大学, 2012.
[32]	Yang FF, Huo LJ, Yang LG, Riaz H, Xiong LR, Chen JG, Zhang SJ, Xiong JJ. Association between melatonin receptor 1A (MTNR1A) gene single-nucleotide polymorphisms and the velvet antler yield of Sika deer. Mol Biol Rep, 2014, 41(4): 1901-1906.
[33]	Chen B. Genetic polymorphism of transformation growth factor β1 and correlation study on antler trait in Sika Deer[Dissertation]. Heilongjiang Bayi Agricultural University, 2011.
[33]	陈斌. 梅花鹿转化生长因子β1基因多态与产茸性状的相关性研究[Dissertation]. 黑龙江八一农垦大学, 2011.
[34]	Li B. Genetic trait of DRB genes and its polymorphism association with antler productivity in sika deer [Dissertation]. Northeast Forestry University, 2008.
[34]	李波. 梅花鹿( Cervus nippon)DRB基因的遗传特性、多态性及与产茸性状的相关性[Dissertation]. 东北林业大学, 2008.
[35]	Du ZH, Liu PY, Bai XJ. Detection of the linkage between RAPD markers and antler production in spotted deer by tail analysis method. China Anim Husb Vet Med, 2008, 35(7): 63-66.
[35]	杜智恒, 刘培源, 白秀娟. 尾分析法检测梅花鹿产茸量性状关联的RAPD标记. 中国畜牧兽医, 2008, 35(7): 63-66.
[36]	Sun HT. Study on microsatellite genetic diversity in Sika Deer of several divisions from Heilongjiang Province [Dissertation]. Daqing: Heilongjiang Bayi Agricultural University, 2009.
[36]	孙海涛. 黑龙江省部分地区梅花鹿微卫星遗传多样性的研究[Dissertation]. 黑龙江八一农垦大学, 2009.
[37]	李馨, 孙海涛, 陈斌, 耿忠诚, 苑立学. 兴凯湖梅花鹿微卫星遗传多样性与鲜茸重相关性分析. 中国畜牧杂志, 2011, 47(9): 14-17.
[38]	Shen YY, Lim BK, Liu HQ, Liu J, Irwin DM, Zhang YP. Multiple episodes of convergence in genes of the dim light vision pathway in bats. PLoS One, 2012, 7(4): e34564.
[39]	Lin T, Zhu GT, Zhang JH, Xu XY, Yu QH, Zheng Z, Zhang ZH, Lun YY, Li S, Wang XX, Huang ZJ, Li JM, Zhang CZ, Wang TT, Zhang YY, Wang AX, Zhang YC, Lin K, Li CY, Xiong GS, Xue YB, Mazzucato A, Causse M, Fei ZJ, Giovannoni JJ, Chetelat RT, Zamir D, St?dler T, Li JF, Ye ZB, Du YC, Huang SW. Genomic analyses provide insights into the history of tomato breeding. Nat Genet, 2014, 46(11): 1220-1226.
[40]	Dar AA, Pradhan TN, Kulkarni DP, Shah SU, Rao KV, Chaukar DA, D'Cruz AK, Chiplunkar SV. Extracellular 2°5°-oligoadenylate synthetase 2 mediates T-cell receptor CD3-ζchain down-regulation via caspase-3 activation in oral cancer. Immunology, 2016, 147(2): 251-264.
[41]	Andersen JB, Strandbyg?rd DJ, Hartmann R, Justesen J. Interaction between the 2°-5° oligoadenylate synthetase- like protein p59 OASL and the transcriptional repressor methyl CpG-binding protein 1. Eur J Biochem, 2004, 271(3): 628-636.
[42]	Justesen J, Hartmann R, Kjeldgaard NO. Gene structure and function of the 2'-5'-oligoadenylate synthetase family. Cell Mol Life Sci, 2000, 57(11): 1593-1612.
[43]	Gromadzka AM, Steckelberg AL, Singh KK, Hofmann K, Gehring NH. A short conserved motif in ALYREF directs cap-and EJC-dependent assembly of export complexes on spliced mRNAs. Nucleic Acids Res, 2016, 44(5): 2348-2361.
[44]	Emdal KB, Pedersen AK, Bekker-Jensen DB, Tsafou KP, Horn H, Lindner S, Schulte JH, Eggert A, Jensen LJ, Francavilla C, Olsen JV. Temporal proteomics of NGF- TrkA signaling identifies an inhibitory role for the E3 ligase Cbl-b in neuroblastoma cell differentiation. Sci Signal, 2015, 8(374): ra40.
[45]	Horikawa I, Fujita K, Jenkins LM, Hiyoshi Y, Mondal AM, Vojtesek B, Lane DP, Appella E, Harris CC. Autophagic degradation of the inhibitory p53 isoform Δ133p53α as a regulatory mechanism for p53-mediated senescence. Nat Commun, 2014, 5: 4706.

编辑推荐

Metrics

www.chinagene.cn
备案号：京ICP备09063187号-4
总访问:,今日访问:,当前在线:

利用全基因组重测序分析鹿茸重量相关基因

Screening weight related genes of velvet antlers by whole genome re-sequencing

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 10

编辑推荐

Metrics

[1]	彭哲也,唐紫珺,谢民主. 机器学习方法在基因交互作用探测中的研究进展[J]. 遗传, 2018, 40(3): 218-226.
[2]	张统雨,朱才业,杜立新,赵福平. 羊重要性状全基因组关联分析研究进展[J]. 遗传, 2017, 39(6): 491-500.
[3]	袁金红, 李俊华, 袁娇娇, 贾克利, 李书粉, 邓传良, 高武军. 基于全基因组测序的MutMap方法在正向遗传学研究中的应用[J]. 遗传, 2017, 39(12): 1168-1177.
[4]	张涛, 王文浩, 张跟喜, 王金玉, 薛倩, 顾玉萍. 京海黄鸡体重性状全基因组关联分析[J]. 遗传, 2015, 37(8): 811-820.
[5]	杨应忠王亚平马兰杜洋格日力. 中国汉族高原肺水肿易感基因的全基因组关联研究[J]. 遗传, 2013, 35(11): 1291-1299.
[6]	李雄伟贾惠娟高中山. 桃基因组及全基因组关/联分析研究进展[J]. 遗传, 2013, 35(10): 1167-1178.
[7]	李聪，孙东晓，姜力，刘剑锋，张勤，张沅，张胜利. 奶牛重要经济性状全基因组关联分析研究进展[J]. 遗传, 2012, 34(5): 545-550.
[8]	郝丽，李和平，严历. 梅花鹿鹿茸尖端组织ESTs分析[J]. 遗传, 2011, 33(4): 371-377.
[9]	杜智恒，白秀娟. 梅花鹿生长激素基因单核苷酸多态与产茸量性状的相关性[J]. 遗传, 2007, 29(3): 337-337―342.
[10]	孙陆果，姜颖，于永利. 吉林双阳型梅花鹿sentrin/SUMO的发现[J]. 遗传, 2002, 24(1): 22-26.