





中国农业科学 ›› 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(
)
收稿日期: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(
), 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] | 叶美金, 吴雷, Lohani Md Nahibuzzaman, 尹丽, 胡欣荣, 刘亚西, 蒋云峰, 陈国跃, 蒲至恩, 李阳, 李婷, 邹亚亚, 吴佳怡, 马建. 基于GWAS的中国地方小麦成熟胚大小位点的鉴定及其遗传效应解析[J]. 中国农业科学, 2026, 59(6): 1157-1171. |
| [2] | 杨丽娟, 陈丝雨, 赵薇, 朱玲, 郭磊, 马丽娜, 马瑞敏, 张娟. 全基因组重测序揭示静原鸡羽色的遗传机制[J]. 中国农业科学, 2026, 59(6): 1348-1360. |
| [3] | 王勇胜, 牛丽, 王长杰, 马立花, 廉潇潇, 孟亚雄, 马小乐, 姚立蓉, 张宏, 杨轲, 李葆春, 王化俊, 司二静, 汪军成. 冬小麦千粒重的全基因组关联分析及候选基因预测[J]. 中国农业科学, 2026, 59(3): 499-514. |
| [4] | 李云丽, 刁邓超, 刘雅睿, 孙玉晨, 孟祥宇, 邬陈芳, 汪妤, 吴建辉, 李春莲, 曾庆东, 韩德俊, 郑炜君. 小麦苗期耐热性全基因组关联分析[J]. 中国农业科学, 2025, 58(9): 1663-1683. |
| [5] | 周广飞, 马亮, 马璐, 张舒钰, 章慧敏, 宋旭东, 张振良, 陆虎华, 郝德荣, 冒宇翔, 薛林, 陈国清. 玉米苞叶性状全基因组关联分析[J]. 中国农业科学, 2025, 58(3): 431-442. |
| [6] | 武书羽, 衡燕芳, 于太飞, 王世佳, 于思佳, 李园, 胡正, 张辉, 孙现军, 黎亮, 姜奇彦. 玉米自然群体苗期耐盐性鉴定及耐盐相关基因分析[J]. 中国农业科学, 2025, 58(20): 4085-4099. |
| [7] | 向爱慧, 白荣基, 郝宇琼, 赵佳佳, 武棒棒, 李晓华, 郑兴卫, 关攀锋, 郑军. 山西小麦矮秆基因的鉴定及株高遗传位点挖掘[J]. 中国农业科学, 2025, 58(17): 3372-3388. |
| [8] | 刘宏祥, 张雪萍, 王逸飞, 王志成, 顾昊天, 宋卫涛, 陶志云, 徐文娟, 章双杰, 卢立志, 李慧芳, 朱春红. 金定鸭产蛋数性状的全基因组关联研究[J]. 中国农业科学, 2025, 58(15): 3145-3158. |
| [9] | 郑敏华, 陈洛, 邢甲乐, 谢月兰, 姜先芽, 聂帅, 蔡甫格, 巫浩翔, 陆展华, 孙伟, 霍兴, 白嵩, 赵均良, 杨武. 华南籼稻稻瘟病抗性QTL鉴定与候选基因挖掘[J]. 中国农业科学, 2025, 58(14): 2707-2719. |
| [10] | 李宁, 高丽锋, 黄鑫, 史华伟, 杨进文, 史雨刚, 陈明, 贾继增, 孙黛珍. 耐低氮小麦品种的筛选及耐低氮指数的全基因组关联分析[J]. 中国农业科学, 2025, 58(13): 2487-2503. |
| [11] | 史顺宇, 杨涛, 庞博, 李静, 林轶峰, 王正瑞, 傅林成, 扎尔加玛丽·阿不都别克, 高文伟, 吴鹏昊. 海岛棉叶绿素含量的全基因组关联分析及候选基因预测[J]. 中国农业科学, 2025, 58(10): 1878-1895. |
| [12] | 张颖, 石婷瑞, 曹瑞, 潘文秋, 宋卫宁, 王利, 聂小军. ICARDA引进-小麦苗期抗旱性的全基因组关联分析[J]. 中国农业科学, 2024, 57(9): 1658-1673. |
| [13] | 赵真坚, 王凯, 陈栋, 申琦, 余杨, 崔晟頔, 王俊戈, 陈子旸, 禹世欣, 陈佳苗, 王翔枫, 唐国庆. 基因组和DNA甲基化组联合分析筛选猪肉质性状关键基因[J]. 中国农业科学, 2024, 57(7): 1394-1406. |
| [14] | 郭军, 邵丹, 窦套存, 马猛, 卢建, 胡玉萍, 王星果, 王强, 李永峰, 郭伟, 童海兵, 曲亮. 鸡产蛋期剩余采食量的随机回归分析及遗传标记筛选[J]. 中国农业科学, 2024, 57(22): 4568-4577. |
| [15] | 郭军, 曲亮, 邵丹, 马猛, 窦套存, 卢建, 胡玉萍, 王星果, 王强, 李永峰, 郭伟, 童海兵. 基于一步法全基因组关联分析解析蛋黄比率遗传结构[J]. 中国农业科学, 2024, 57(21): 4356-4366. |
|
||