绵羊体重性状全基因组关联分析

摘要/Abstract

摘要： 本研究旨在利用全基因组关联分析方法（genome-wide association studies，GWAS）挖掘影响绵羊体重性状的遗传标记和功能基因。采用Illumina OvineSNP50 BeadChip中密度商业化羊芯片，对3个品种（苏尼特羊、德国肉用美利奴羊和杜泊羊）共计329只绵羊的体重性状（初生重、断奶重、6月龄重、断奶前日增重、断奶后日增重、日增重）进行GWAS分析。统计分析和基因注释基于TASSEL软件、混合线性模型和2012年10月公布的最新版绵羊基因组Ovis_aries_v3.1序列信息。结果表明，有10个SNPs位点在全基因组显著水平上与断奶后日增重相关，部分位点分布在羊注释基因内（MEF2B、RFXANK等）；除此之外还检测到22个SNPs在染色体显著水平上与其他体重性状存在相关性，获得一批重要的候选基因。这些基因对绵羊体重性状功能基因的挖掘具有重要的理论意义和参考价值。

关键词: 绵羊; 体重性状; 全基因组关联分析

Abstract: Motivated by mining major candidate genes across Ovine genome, the present study is to perform genome-wide association studies(GWAS) to detect genes associated with body weight traits. Using Illumina OvineSNP50 BeadChip, we performed a GWA study in 329 purebred sheep phenotyped for 6 body weight traits(birth weight, weaning weight, 6-month weight, pre-weaning gain, post-weaning gain, daily weight gain). Statistics and data analysis were based on TASSEL program，mixed linear model and the latest Ovis_aries_v3.1 genome sequence (released October 2012). The results indicated that 10 SNPs consistently reached genome-wise significant level for post-weaning gain and 22 SNPs reached chromosome-wise significant level for other body weight traits. The SNPs were within (MEF2B,RFXANK，et al) or close to some ovine genes, which were thought to be the most important candidate genes associated with body weight traits. The results will contribute to identify candidate genes for ovine body weight traits, and facilitate the potential utilization of genes involved production traits in sheep in future.

Key words: sheep; body weight traits; genome-wide association studies

张莉，刘佳森，徐凌洋，赵福平，陆健，张世芳，王慧华，张晓宁，魏彩虹，陆国彬，郑友民，杜立新. 绵羊体重性状全基因组关联分析[J]. 中国畜牧兽医.

ZHANG Li，LIU Jia-sen，XU Ling-yang，ZHAO Fu-ping，LU Jian，ZHANG Shi-fang，WANG Hui-hua，ZHANG Xiao-ning，WEI Cai-hong，LU Guo-bin，ZHENG You-min，DU Li-xin. Genome-wide Association Studies for Body Weight Traits in Sheep[J]. .

参考文献

1 何波. 猪骨骼肌卫星细胞的培养、鉴定及猪MEF2B和MEF2C基因的初步研究[D]. 武汉：华中农业大学，2006.
2 何波，郑嵘，熊远著，等. 新生猪骨骼肌卫星细胞的培养鉴定及生物学特性[J]. 畜牧兽医学报，2006,37（6）：555～559.
3 郑水儿. 核糖体蛋白L7A在骨肉瘤组织中的表达及临床意义[J]. 中国癌症杂志，2009，19(2): 101～105.
4 Boon K, Caron H N, van Asperen R, et al. N-myc enhances the expression of a large set of genes functioning in ribosome biogenesis and protein synthesis[J]. EMBO J，2001，20(6): 1383～1393.
5 Bowers P M, Cokus S J, Elsenberg D, et al. Use of logic relationships to decipher protein network organization[J]. Science，2004，306(5705): 2246～2249.
6 Cavanagh C R, Jonas E, Hobbs M, et al. Mapping quantitative trait loci (QTL) in sheep. Ⅲ. QTL for carcass composition traits derived from CT scans and aligned with a meta-assembly for sheep and cattle carcass QTL[J]. Genet Sel Evol，2010，42: 36.
7 Clop A, Marcq F, Takeda H, et al. A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep[J]. Nat Genet，2006，38(7): 813～818.
8 Diarra-Mehrpour M, Arrabal S, Jalil A, et al. Prion protein prevents human breast carcinoma cell line from tumor necrosis factor alpha-induced cell death[J]. Cancer Res，2004，64(2): 719～727.
9 Edmondson D G, Lyons G E, Martin J F, et al. Mef2 gene expression marks the cardiac and skeletal muscle lineages during mouse embryogenesis[J]. Development，1994，120(5): 1251～1263.
10 Freudenberg J, Lee H S, Han B G, et al. Genome-wide association study of rheumatoid arthritis in Koreans: Population-specific loci as well as overlap with European susceptibility loci[J]. Arthritis Rheum，2011，63(4): 884～893.
11 Garcia-Gamez E, Gutierrez-Gil B, Sahana G, et al. GWA analysis for milk production traits in dairy sheep and genetic support for a QTN influencing milk protein percentage in the LALBA gene[J]. PLoS One，2012，7(10): e47782.
12 Gorbach D M, Fan B, Onteru S K, et al. Genome-wide association studies for important economic traits in domestic animals using high density SNP genotyping. Iowa State University Animal Industry Report，2010.
13 Hillier L W, Graves T A, Fulton R S, et al. Generation and annotation of the DNA sequences of human chromosomes 2 and 4[J]. Nature，2005，434(7034): 724～731.
14 Jiang L, Liu J, Sun D, et al. Genome wide association studies for milk production traits in Chinese Holstein population[J]. PLoS One，2010，5(10): e13661.
15 Johnston S E, McEwan J C, Pickering N K, et al. Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population[J]. Mol Ecol，2011，20(12): 2555～2566.
16 Kim M Y, Park E, Park J H, et al. Expression profile of nine novel genes differentially expressed in hepatitis B virus-associated hepatocellular carcinomas[J]. Oncogene，2001，20(33): 4568～4575.
17 Lee S J. Regulation of muscle mass by myostatin[J]. Annu Rev Cell Dev Biol,2004,20: 61～86.
18 Lennon-Dumenil A M, Barbouche M R, Vedrenne J, et al. Uncoordinated HLA-D gene expression in a RFXANK-defective patient with MHC class Ⅱ deficiency[J]. J Immunol，2001，166(9): 5681～5687.
19 Lilly B, Zhao B, Ranganayakulu G, et al. Requirement of MADS domain transcription factor D-MEF2 for muscle formation in Drosophila[J]. Science，1995，267(5198): 688～693.
20 Magnani R, Dirk L M, Trievel R C, et al. Calmodulin methyltransferase is an evolutionarily conserved enzyme that trimethylates Lys-115 in calmodulin[J]. Nat Commun，2010，1: 43.
21 Masternak K, Barras E, Zufferey M, et al. A gene encoding a novel RFX-associated transactivator is mutated in the majority of MHC class Ⅱ deficiency patients[J]. Nat Genet，1998，20(3): 273～277.
22 Meng Q, Zheng M, Liu H, et al. TRAF6 regulates proliferation, apoptosis, and invasion of osteosarcoma cell[J]. Mol Cell Biochem，2012，371(1～2): 177～186.
23 Molkentin J D, Black B L, Martin J F, et al. Cooperative activation of muscle gene expression by MEF2 and myogenic bHLH proteins[J]. Cell，1995，83(7): 1125～1136.
24 Nagel J E, Smith R J, Shaw L, et al. Identification of genes differentially expressed in T cells following stimulation with the chemokines CXCL12 and CXCL10[J]. BMC Immunol，2004，5: 17.
25 NYU. How different species can use the same genes yet develop distinct features[J]. Science Daily，2007，3：47.
26 O’Toole M, Janszen D B, Slonim D K, et al. Risk factors associated with beta-amyloid(1-42) immunotherapy in preimmunization gene expression patterns of blood cells[J]. Arch Neurol，2005，62(10): 1531～1536.
27 Palmarini M, Sharp J M, Lee C, et al. In vitro infection of ovine cell lines by Jaagsiekte sheep retrovirus[J]. J Virol，1999，73(12): 10070～10078.
28 Payne L N, Nair V. The long view: 40 years of avian leukosis research[J]. Avian Pathol，2012，41(1): 11～19.
29 Pepin M, Vitu C, Russo P, et al. Maedi-visna virus infection in sheep: A review[J]. Vet Res，1998，29(3～4): 341～367.
30 Purdie A C, Plain K M, Begg D J, et al. Candidate gene and genome-wide association studies of Mycobacterium avium subsp. paratuberculosis infection in cattle and sheep: A review[J]. Comp Immunol Microbiol Infect Dis，2011，34(3): 197～208.
31 Qu Y, Sun H, Sun M, et al. Sequence analysis for the complete proviral genome of avian leukosis virus subgroup J associated with haemangiomas, leiomyosarcomas and myelomas in layer flocks[J]. Acta Vet Hung，2012，60(3): 343～354.
32 Rader K, Boyer A D, Farquhar M G, et al. Assignment of ankyrin repeat, family A (RFXANK-like) 2 (ANKRA2) to human chromosome 5q12-->q13 by radiation hybrid mapping and somatic cell hybrid PCR[J]. Cytogenet Cell Genet，2000，89(3～4): 164～165.
33 Ranganayakulu G, Zhao B, Dokidis A, et al. A series of mutations in the D-MEF2 transcription factor reveal multiple functions in larval and adult myogenesis in Drosophila[J]. Dev Biol，1995，171(1): 169～181.
34 Reith W，Mach B. The bare lymphocyte syndrome and the regulation of MHC expression[J]. Annu Rev Immunol，2001，19: 331～373.
35 Riggio V, Matika O, Pong-Wong R, et al. Genome-wide association and regional heritability mapping to identify loci underlying variation in nematode resistance and body weight in Scottish Blackface lambs[J]. Heredity (Edinb)，2013，110(5): 420～429.
36 Roche A, Ripoll G, Joy M, et al. Effects of the FecX(R) allele of BMP15 gene on the birth weight, growth rate and carcass quality of Rasa Aragonesa light lambs[J]. Small Ruminant Research，2012，108(1～3): 45～53.
37 Schomburg L, Turwitt S, Prescher G, et al. Human TRH-degrading ectoenzyme cDNA cloning, functional expression, genomic structure and chromosomal assignment[J]. Eur J Biochem，1999，265(1): 415～422.
38 Shay T L, Berghmans S, Segers K, et al. Fine-mapping and construction of a bovine contig spanning the ovine callipyge locus[J]. Mammalian Genome，2001，12(2): 141～149.
39 Straub O C. Maedi-Visna virus infection in sheep. History and present knowledge[J]. Comp Immunol Microbiol Infect Dis，2004，27(1): 1～5.
40 Subramanian S V，Nadal-Ginard B. Early expression of the different isoforms of the myocyte enhancer factor-2 (MEF2) protein in myogenic as well as non-myogenic cell lineages during mouse embryogenesis[J]. Mech Dev，1996，57(1): 103～112.
41 Wang Y,Liu S, Han M, et al. Cloning and sequence analysis of pol gene of the Inner Mongolia strain of endogenous jaagsiekte sheep retrovirus[J]. Chinese Journal of Veterinary Medicine，2007，43(12): 8～10.
42 Wool I G, Chan Y L，Gluck A. Structure and evolution of mammalian ribosomal proteins[J]. Biochem Cell Biol，1995，73(11～12): 933～947.
43 Yokoyama M, Nishi Y, Yoshii J, et al. Identification and cloning of neuroblastoma-specific and nerve tissue-specific genes through compiled expression profiles[J]. DNA Res，1996，3(5): 311～320.
44 Zhang Z, Ersoz E, Lai C Q, et al. Mixed linear model approach adapted for genome-wide association studies[J]. Nat Genet，2010，42(4): 355～360.
45 Zhao X, Dittmer K E, Blair H T, et al. A novel nonsense mutation in the DMP1 gene identified by a genome-wide association study is responsible for inherited rickets in Corriedale sheep[J]. PLoS One，2011，6(7): e21739.

[1]	王赛桥, 赵路, 翟振翰, 张丙雷, 贾万航, 王玉琴. 绵羊miR-449a/b前体序列组织表达及生物信息学分析[J]. 中国畜牧兽医, 2023, 50(6): 2175-2184.
[2]	刘锦钰, 张杨, 马旭梅, 周千, 周扬, 张春香, 石磊, 任有蛇. 谷胱甘肽对绵羊精液低温保存效果的影响[J]. 中国畜牧兽医, 2023, 50(6): 2388-2394.
[3]	余大新, 刘文静, 刘大鹏, 刘同, 张鹤, 盖凯, 李姝楠, 张贝贝, 侯水生, 周正奎. 鸭胸肌矿物元素含量全基因组关联分析[J]. 中国畜牧兽医, 2023, 50(4): 1287-1295.
[4]	崔双, 乌日力嘎, 阿拉达尔, 晓敏, 王剑秋, 弓剑. 放牧条件下补充酵母硒对绵羊胃肠道微生物区系的影响[J]. 中国畜牧兽医, 2023, 50(3): 952-965.
[5]	马建青, 杨清芳, 宋占锋, 贾丁丁, 吴东蔚, 莫秀芳, 徐增年, 杨克, 吴占勇, 籍颖, 周荣艳. 绵羊骨代谢相关基因VDR密码子偏好性分析[J]. 中国畜牧兽医, 2022, 49(8): 2888-2898.
[6]	张雪萍, 石斌刚, 金夏阳, 王向彦, 兰丽娟, 时钰, 祁有鹏, 赵世杰, 李少斌, 胡江. 藏绵羊KLF7基因表达特征分析及其过表达对前脂肪细胞增殖分化的影响[J]. 中国畜牧兽医, 2022, 49(8): 2920-2930.
[7]	李卫东, 蒋玉婷, 段宏伟, 何海军, 杨帅, 丁自强, 吴建新, 张恩秋, 胡俊杰. 褪黑素受体MT1和MT2在绵羊胃中的分布及表达[J]. 中国畜牧兽医, 2022, 49(8): 3256-3262.
[8]	严啊妮, 沈奔, 王洪荣, 冯春燕, 魏宗友. 过瘤胃葡萄糖和缓释尿素对热应激绵羊生产性能及瘤胃发酵功能的影响[J]. 中国畜牧兽医, 2022, 49(6): 2105-2115.
[9]	张丽萌, 刘爱菊, 李闰婷, 李玉华, 李林, 聂晓宁, 王林青, 陈龙欣. 绵羊CREBRF基因克隆、生物信息学及组织表达分析[J]. 中国畜牧兽医, 2022, 49(5): 1599-1609.
[10]	强浩, 梁鹏, 孟科, 荣轩, 冯登侦. 绵羊IGF1R基因多态性及其与生长性状的关联分析[J]. 中国畜牧兽医, 2022, 49(5): 1765-1785.
[11]	苏晓慧, 何海迎, 吐尔逊阿依·牙力坤, 王琼, 马海玉, 刘玲玲, 刘武军. 野生盘羊×巴什拜羊杂交二代尾部脂肪差异表达lncRNA的鉴定与分析[J]. 中国畜牧兽医, 2022, 49(5): 1796-1805.
[12]	田广原, 王宇琛, 周雅坪, 郭婷, 赵红梅, 赵逢淼, 孙雅婕, 卞宇晨, 于佳梁, 周雨霞. 绵羊肺炎支原体P60_58aa-928aa蛋白的原核表达及间接ELISA方法的建立[J]. 中国畜牧兽医, 2022, 49(5): 1960-1969.
[13]	翟哲, 陈思, 吴艳茹, 王雪梅, 李崇瑞, 刘志勇, 陈巧玲, 王凤阳, 杜丽, 李昌, 金宁一. 绵羊CCL19基因启动子活性分析及其转录调控元件筛选[J]. 中国畜牧兽医, 2022, 49(4): 1213-1222.
[14]	李静云, 王新乐, 张云飞, 白俊艳, 李淦, 雷莹, 董智豪, 陈宇, 樊红灯, 王龙威, 陈梦柯, 王兴科, 司马磊阳, 张世文. 绵羊MyoG基因内含子Ⅱ多态性及其与生长性状的关联分析[J]. 中国畜牧兽医, 2022, 49(4): 1364-1373.
[15]	张永浩, 肖逸, 郭丽君, 阿布都热合曼·吐尔逊, 王玉涛. 马可·波罗盘羊和5种家养绵羊MHC-DRB3基因第2外显子遗传多态性分析[J]. 中国畜牧兽医, 2022, 49(3): 962-972.