





 
			中国农业科学 ›› 2023, Vol. 56 ›› Issue (5): 995-1006.doi: 10.3864/j.issn.0578-1752.2023.05.014
        
               		杨明路1,2( ), 张海亮1, 罗汉鹏1, 黄锡霞2, 张翰林4, 章施施4, 王炎4, 刘林3, 郭刚4(
), 张海亮1, 罗汉鹏1, 黄锡霞2, 张翰林4, 章施施4, 王炎4, 刘林3, 郭刚4( ), 王雅春1(
), 王雅春1( )
)
                  
        
        
        
        
    
收稿日期:2021-10-08
									
				
									
				
											接受日期:2022-12-28
									
				
											出版日期:2023-03-01
									
				
											发布日期:2023-03-13
									
			通信作者:
					郭刚,E-mail:guogang2180@126.com。王雅春,E-mail:wangyachun@cau.edu.cn
							联系方式:
				杨明路,E-mail:2510660222@qq.com。
				
							基金资助:
        
               		YANG  MingLu1,2( ), ZHANG  HaiLiang1, LUO  HanPeng1, HUANG  XiXia2, ZHANG  HanLin4, ZHANG  ShiShi4, WANG  Yan4, LIU  Lin3, GUO  Gang4(
), ZHANG  HaiLiang1, LUO  HanPeng1, HUANG  XiXia2, ZHANG  HanLin4, ZHANG  ShiShi4, WANG  Yan4, LIU  Lin3, GUO  Gang4( ), WANG  YaChun1(
), WANG  YaChun1( )
)
			  
			
			
			
                
        
    
Received:2021-10-08
									
				
									
				
											Accepted:2022-12-28
									
				
											Published:2023-03-01
									
				
											Online:2023-03-13
									
			摘要:
【目的】探究基于奶牛智能项圈系统产生的发情指数和活动峰值指标的群体规律及其影响因素,对其进行遗传参数估计,并利用全基因组关联分析(genome-wide association study, GWAS)挖掘与奶牛发情相关指标的遗传标记,以期为从遗传育种角度制定选育策略提高奶牛繁殖力提供有用信息。【方法】试验以北京市某规模化牧场健康的泌乳期荷斯坦牛为研究对象,收集了该牧场2017年7月至2020年5月2 074 头泌乳牛的发情记录,包括发情指数和活动峰值,同时收集了试验牛的产犊记录;采用 SAS 9.4中MIXED过程,分析了非遗传因素对发情指数和活动峰值的影响;基于DMU软件DMUAI模块,使用双性状动物模型估计了发情指数与活动峰值的遗传力及遗传相关,并计算了发情指数与常规繁殖性状的近似遗传相关;以逆回归育种值为关联分析的表型,利用FarmCPU软件分别对发情指数和活动峰值进行了GWAS。【结果】荷斯坦牛发情指数为(71.77±19.13)au,活动峰值为(105.94±30.73)au/2h。发情指数遗传力为0.04±0.01,活动峰值遗传力为0.19±0.04,均为低遗传力性状;发情指数与活动峰值之间的遗传相关为0.45±0.03;发情指数与经产牛首末次配种间隔、首次产犊日龄和青年牛首末次配种间隔的近似遗传相关分别为0.37、-0.41、-0.55。两个性状的GWAS共找到7个在全基因组水平上显著的SNP位点,分别位于7条染色体上,显著位点附近(300 kb内)共发现基因31个。【结论】通过项圈系统获得的发情指数和活动峰值是可遗传的,两者之间存在中等遗传相关,活动峰值为中低遗传力,且能够一定程度上反映奶牛的发情指数,活动峰值和发情指数为研究奶牛的发情行为以及提高母牛的繁殖力提供了新的方向。GWAS鉴定到的候选基因PTGS1(Prostaglandin-Endoperoxide Synthase1)和NDUFA8(Ubiquinone Oxidoreductase Subunit A8)的功能与供体胚胎总数有关,VEGFA(Vascular Endothelial Growth Factor A)基因功能与奶牛繁殖有关,可作为奶牛发情相关指标的候选基因进行深入研究。
杨明路, 张海亮, 罗汉鹏, 黄锡霞, 张翰林, 章施施, 王炎, 刘林, 郭刚, 王雅春. 基于智能项圈系统荷斯坦牛发情相关指标的遗传参数估计及全基因组关联分析[J]. 中国农业科学, 2023, 56(5): 995-1006.
YANG MingLu, ZHANG HaiLiang, LUO HanPeng, HUANG XiXia, ZHANG HanLin, ZHANG ShiShi, WANG Yan, LIU Lin, GUO Gang, WANG YaChun. Estimation of Genetic Parameters and Genome-Wide Association Study of Heat Indicators in Holstein Cattle Based on Collar-Mounted Device[J]. Scientia Agricultura Sinica, 2023, 56(5): 995-1006.
 
												
												表2
各因素对中国荷斯坦牛发情指数和活动峰值的影响"
| 效应 Effect | 水平 Level | 发情指数Heat index | 活动峰值Activity peak | ||
|---|---|---|---|---|---|
| 数据量 N | LSM±SE | 数据量 N | LSM±SE | ||
| 胎次 Parity | 1 | 2 378 | 70.79±0.53b | 773 | 108.82±1.66a | 
| 2 | 2 707 | 72.05±0.50a | 978 | 108.74±1.58a | |
| ≥3 | 4 226 | 71.14±0.46a | 886 | 105.23±1.66c | |
| 泌乳阶段 Lactation | (0 -21 d) | 327 | 67.65±1.12c | 92 | 104.14±3.30b | 
| (22-100 d) | 4 798 | 74.66±0.40a | 1 389 | 112.63±1.33a | |
| (101 -200 d) | 2 979 | 73.11±0.41b | 829 | 109.47±1.36a | |
| (201-305d) | 1 207 | 69.91±0.61c | 327 | 104.13±1.93b | |
| 情期数 Number of Estrous cycles | 1 | 3 086 | 70.15±0.45c | 850 | 108.77±1.42a | 
| 2 | 2 275 | 70.46±0.53c | 669 | 108.82±1.61a | |
| 3 | 1 578 | 72.22±0.60b | 456 | 105.32±1.78b | |
| ≥4 | 2 372 | 72.49±0.58a | 662 | 107.46±1.82a | |
| 测定季节 Test season | 春季(3—5月) Spring(Mar-May) | 2 054 | 71.23±0.55b | 632 | 108.15±1.66a | 
| 夏季(6—8月) Summer(Jun-Aug) | 1 753 | 69.34±0.59c | 503 | 101.67±1.85b | |
| 秋季(9—11月) Autumn(Sep-Nov) | 3 062 | 72.38±0.50a | 784 | 108.50±1.60a | |
| 冬季(12—2月) Winter(Dec-Feb) | 2 442 | 72.38±0.50a | 718 | 112.05±1.59a | |
| 测定年份 Test year | 2017 | 2 054 | 70.94±0.84b | 142 | 112.82±2.83a | 
| 2018 | 1 753 | 67.22±0.41c | 858 | 107.36±1.44b | |
| 2019 | 3 062 | 72.80±0.43a | 1 229 | 107.17±1.40b | |
| 2020 | 2 442 | 74.36±0.67a | 408 | 103.02±2.06c | |
| 昼夜 Day-and-night | 夜间Day | 4 918 | 72.63±0.63a | 1 450 | 104.95±1.40b | 
| 白天Night | 4 393 | 69.83±0.64b | 1 187 | 110.23±1.44a | |
| 项圈类型 Heatime type | HR-LD | 6 216 | 68.83±0.42b | 554 | 116.47±1.85a | 
| HR-LDn | 3 095 | 73.84±0.49a | 2 083 | 98.71±1.30b | |
 
												
												表4
中国荷斯坦牛发情指数与繁殖相关性状近似遗传相关"
| 性状 Trait | 个体数 N | 近似遗传相关 Approximate genetic correlation | 
|---|---|---|
| 首次产犊日龄 Age at first calving in heifer | 412 | -0.41±0.05 | 
| 首次配种日龄 Age at first service in heifer | 419 | 0.06±0.05 | 
| 产犊至首次配种间隔 Interval from calving to first service in cow | 408 | 0.06±0.05 | 
| 经产牛首末次配种间隔 Interval from first to last insemination in cow | 404 | 0.37±0.05 | 
| 青年牛首末次配种间隔 Interval from first to last insemination in heifer | 107 | -0.55±0.08 | 
 
												
												表5
中国荷斯坦牛发情指数和活动峰值全基因组关联分析显著SNP和相关基因"
| 性状 Trait | SNP | 染色体 Chromosome | 位置 Position (bp) | 最小等位基因频率 MAF | 候选基因 Harbor genes | P值 P value | 
|---|---|---|---|---|---|---|
| 发情指数 Heat index | BTB-01605503 | 10 | 12659662 | 0.12 | RAB11A、MIR2290、DENND4A、SLC24A1、HACD3、INTS14、DPP8 | 1.69E-07 | 
| Hapmap51151-BTA-114546 | 11 | 93285949 | 0.14 | PTGS1、MRRF、RBM18、LHX6、MORN5、NDUFA8 | 9.57E-09 | |
| BovineHD1900003595 | 19 | 13596512 | 0.43 | ACACA、C19H17orf78、TADA2A、DUSP14、SYNRG、LHX1 | 2.77E-07 | |
| BovineHD3000015447 | 30 | 52853475 | 0.36 | — | 2.61E-07 | |
| 活动峰值 Activity peak | BovineHD2100020097 | 21 | 68740864 | 0.33 | ASPG、ATP5MPL | 2.39E-09 | 
| BovineHD2200007296 | 22 | 24876308 | 0.30 | CNTN6 | 2.83E-07 | |
| BovineHD2300004473 | 23 | 17552613 | 0.32 | AARS2、TMEM151B、SLC35B2、HSP90AB1、SLC29A1、TMEM63B、MRPL14、NFKBIE、VEGFA | 3.70E-08 | 
| [1] | doi: 10.3168/jds.S0022-0302(05)72792-2 pmid: 15738259 | 
| [2] | 刘澳星. 奶牛繁殖性状遗传参数估计与基因组预测[D]. 北京: 中国农业大学, 2018. | 
|  | |
| [3] | doi: 10.3168/jds.2007-0178 pmid: 18349259 | 
| [4] | doi: 10.3168/jds.2014-8940 pmid: 25747826 | 
| [5] | doi: 10.3168/jds.2008-1556 pmid: 19233816 | 
| [6] | 陈紫薇, 师睿, 罗汉鹏, 田佳, 魏趁, 张伟新, 李委奇, 温万, 王雅晶, 王雅春. 宁夏地区荷斯坦牛青年牛繁殖性状遗传参数估计. 畜牧兽医学报, 2021, 52(2): 344-351. doi:10.11843/j.issn.0366-6964.2021.02.007.  doi: 10.11843/j.issn.0366-6964.2021.02.007 | 
| doi: 10.11843/j.issn.0366-6964.2021.02.007 | |
| [7] | doi: 10.3168/jds.2008-1918 pmid: 19620690 | 
| [8] | 胡丽蓉, 康玲, 王淑慧, 李玮, 鄢新义, 罗汉鹏, 董刚辉, 王新宇, 王雅春, 徐青. 冷热应激对北京地区荷斯坦牛产奶性能及血液生化指标的影响. 中国农业科学, 2018, 51(19): 3791-3799. doi:10.3864/j.issn.0578-1752.2018.19.015.  doi: 10.3864/j.issn.0578-1752.2018.19.015 | 
| doi: 10.3864/j.issn.0578-1752.2018.19.015 | |
| [9] | doi: S0093-691X(20)30368-X pmid: 32721698 | 
| [10] | doi: 10.3168/jds.2021-20437 | 
| [11] | 刘丽元, 周靖航, 张梦华, 李金霞, 方季青, 谭世新, 王爱芳, 黄锡霞, 王雅春. 新疆褐牛基因外显子多态性及其与SCS和泌乳性状关联性分析. 中国农业科学, 2017, 50(13): 2592-2603. doi:10.3864/j.issn.0578-1752.2017.13.017.  doi: 10.3864/j.issn.0578-1752.2017.13.017 | 
| doi: 10.3864/j.issn.0578-1752.2017.13.017 | |
| [12] | 王晓, 张勤, 俞英. 基于体细胞评分的中国荷斯坦牛乳房炎抗性全基因组关联分析. 中国农业科学, 2017, 50(4): 755-763. doi:10.3864/j.issn.0578-1752.2017.04.015.  doi: 10.3864/j.issn.0578-1752.2017.04.015 | 
| doi: 10.3864/j.issn.0578-1752.2017.04.015 | |
| [13] | 王小龙, 常玲玲, 陈莹, 陈仁金, 李云龙, 毛永江, 冀德君, 杨章平. 中国荷斯坦牛SCD1外显子5基因多态性及其对乳中脂肪酸组成的影响. 中国农业科学, 2014, 47(9): 1858-1864. doi:10.3864/j.issn.0578-1752.2014.09.021.  doi: 10.3864/j.issn.0578-1752.2014.09.021 | 
| doi: 10.3864/j.issn.0578-1752.2014.09.021 | |
| [14] | 昝林森, 张佳兰, 刘新武. 牛AGPAT6基因遗传特征与奶牛产奶性能相关性研究. 中国农业科学, 2007, 40(7): 1498-1503. doi:10.3321/j.issn:0578-1752.2007.07.026.  doi: 10.3321/j.issn:0578-1752.2007.07.026 | 
| doi: 10.3321/j.issn:0578-1752.2007.07.026 | |
| [15] | doi: 10.1046/j.1439-0388.2001.00302.x | 
| [16] | doi: 10.3168/jds.2011-4653 pmid: 22541472 | 
| [17] | doi: 10.1186/s40104-017-0189-x pmid: 28680590 | 
| [18] | doi: 10.3168/jds.2020-19192 pmid: 33685707 | 
| [19] | 王嫣然, 张学霞, 赵静瑶, 余新晓, 姜群鸥. 北京地区不同季节PM2.5和PM10浓度对地面气象因素的响应. 中国环境监测, 2017, 33(2): 34-41. doi:10.19316/j.issn.1002-6002.2017.02.06.  doi: 10.19316/j.issn.1002-6002.2017.02.06 | 
| doi: 10.19316/j.issn.1002-6002.2017.02.06 | |
| [20] | doi: 10.1016/j.ajhg.2009.01.005 | 
| [21] | doi: 10.1371/journal.pone.0175971 | 
| [22] |  | 
| [23] | doi: 10.2527/jas.2006-631 pmid: 17339413 | 
| [24] | doi: 10.2527/jas1973.373676x | 
| [25] |  | 
| [26] | doi: 10.1371/journal.pgen.1005767 | 
| [27] | doi: 10.3389/fgene.2021.634986 | 
| [28] |  | 
| [29] | doi: 10.3168/jds.2008-1736 pmid: 19700732 | 
| [30] | 武玥, 王宪龙, 李锡智, 曾申明. 北京地区季节变化对奶牛发情和受胎的影响. 中国畜牧杂志, 2011, 47(13): 75-78. | 
|  | |
| [31] | doi: 10.1530/jrf.0.0990315 pmid: 8107012 | 
| [32] | 孙保贵. 奶牛运动量辅助发情诊断及电导率辅助乳房炎诊断参数的研究[D]. 泰安: 山东农业大学, 2011. | 
|  | |
| [33] | 胡仁超. 荷斯坦牛活动量变化规律及其影响因素的相关性研究[D]. 杨凌: 西北农林科技大学, 2016. | 
|  | |
| [34] | 潘予琮, 蒋林树, 熊本海. 奶牛产后首次发情行为变化规律及影响因素分析. 中国畜牧杂志, 2021, 57(4): 123-128. doi:10.19556/j.0258-7033.20200626-03.  doi: 10.19556/j.0258-7033.20200626-03 | 
| doi: 10.19556/j.0258-7033.20200626-03 | |
| [35] | 寇红祥, 李蓝祁, 王振玲, 任康, 李树静, 路永强, 王栋. 牛发情期活动量与阴道黏液电阻值变化规律的研究. 畜牧兽医学报, 2017, 48(7): 1221-1228. doi:10.11843/j.issn.0366-6964.2017.07.006.  doi: 10.11843/j.issn.0366-6964.2017.07.006 | 
| doi: 10.11843/j.issn.0366-6964.2017.07.006 | |
| [36] | doi: 10.3168/jds.2006-634 pmid: 17369207 | 
| [37] | 鄢新义, 董刚辉, 徐伟, 刘澳星,  doi: 10.11843/j.issn.0366-6964.2016.05.012 | 
| doi: 10.11843/j.issn.0366-6964.2016.05.012 | |
| [38] | doi: 10.3168/jds.2008-1721 pmid: 20059923 | 
| [39] | doi: 10.1016/j.livsci.2018.04.018 | 
| [40] | doi: S0022-0302(18)30416-8 pmid: 29753485 | 
| [41] | doi: 10.1016/j.anireprosci.2020.106466 | 
| [42] | doi: 10.1007/s004390050869 | 
| [43] | doi: 10.1095/biolreprod.112.102921 | 
| [44] | doi: 10.3109/19396368.2012.656217 pmid: 22313243 | 
| [45] |  | 
| [46] |  | 
| [47] | doi: 10.1093/jas/skz164 | 
| [1] | 谭力治, 赵毅强. 全基因组关联分析中混合模型的原理、优化与应用[J]. 中国农业科学, 2023, 56(9): 1617-1632. | 
| [2] | 王慧玲, 闫爱玲, 王晓玥, 刘振华, 任建成, 徐海英, 孙磊. 葡萄果粒质量相关性状全基因组关联分析[J]. 中国农业科学, 2023, 56(8): 1561-1573. | 
| [3] | 王脉, 董清峰, 高珅奥, 刘德政, 卢山, 乔朋放, 陈亮, 胡银岗. 小麦苗期根系性状的全基因组关联分析与优异位点挖掘[J]. 中国农业科学, 2023, 56(5): 801-820. | 
| [4] | 孙艳发, 吴琼, 林如龙, 陈红萍, 甘秋云, 沈玥, 王亚茹, 薛鹏飞, 陈飞帆, 刘健涛, 周陈鑫, 兰诗诗, 潘浩哲, 邓凡, 岳稳, 江宵兵, 李焰. 龙岩山麻鸭蛋品质性状的全基因组关联研究[J]. 中国农业科学, 2023, 56(3): 572-586. | 
| [5] | 李周帅,董远,李婷,冯志前,段迎新,杨明羡,徐淑兔,张兴华,薛吉全. 基于杂交种群体的玉米产量及其配合力的全基因组关联分析[J]. 中国农业科学, 2022, 55(9): 1695-1709. | 
| [6] | 职蕾,者理,孙楠楠,杨阳,Dauren Serikbay,贾汉忠,胡银岗,陈亮. 小麦苗期铅耐受性的全基因组关联分析[J]. 中国农业科学, 2022, 55(6): 1064-1081. | 
| [7] | 王凯,张海亮,董祎鑫,陈少侃,郭刚,刘林,王雅春. 基于牧场管理数据的奶牛健康性状定义及遗传参数估计[J]. 中国农业科学, 2022, 55(6): 1227-1240. | 
| [8] | 逄洪波, 程露, 于茗兰, 陈强, 李玥莹, 吴隆坤, 王泽, 潘孝武, 郑晓明. 栽培稻芽期耐低温全基因组关联分析[J]. 中国农业科学, 2022, 55(21): 4091-4103. | 
| [9] | 竹磊,张海亮,陈少侃,安涛,罗汉鹏,刘林,黄锡霞,王雅春. 荷斯坦牛泌乳早期体细胞数对生产性能的影响及遗传参数估计[J]. 中国农业科学, 2022, 55(2): 403-414. | 
| [10] | 谢晓宇, 王凯鸿, 秦晓晓, 王彩香, 史春辉, 宁新柱, 杨永林, 秦江鸿, 李朝周, 马麒, 宿俊吉. 陆地棉吐絮率的限制性两阶段多位点全基因组关联分析及候选基因预测[J]. 中国农业科学, 2022, 55(2): 248-264. | 
| [11] | 常立国,何坤辉,刘建超. 多环境下玉米保绿相关性状遗传位点的挖掘[J]. 中国农业科学, 2022, 55(16): 3071-3081. | 
| [12] | 李婷,董远,张君,冯志前,王亚鹏,郝引川,张兴华,薛吉全,徐淑兔. 玉米杂交种穗部性状的全基因组关联分析[J]. 中国农业科学, 2022, 55(13): 2485-2499. | 
| [13] | 王娟, 马晓梅, 周小凤, 王新, 田琴, 李成奇, 董承光. 棉花产量构成因素性状的全基因组关联分析[J]. 中国农业科学, 2022, 55(12): 2265-2277. | 
| [14] | 崔承齐, 刘艳阳, 江晓林, 孙知雨, 杜振伟, 武轲, 梅鸿献, 郑永战. 芝麻产量相关性状的多位点全基因组关联分析及候选基因预测[J]. 中国农业科学, 2022, 55(1): 219-232. | 
| [15] | 张鹏霞,周秀文,梁雪,郭营,赵岩,李斯深,孔凡美. 小麦苗期生物量及氮效率相关性状的全基因组关联分析[J]. 中国农业科学, 2021, 54(21): 4487-4499. | 
| 
 | ||