2. 中国农业科学院北京畜牧兽医研究所 农业部畜禽遗传资源与利用 重点开放实验室, 北京 100193;
3. 吉林农业大学动物科学技术学院, 长春 130118;
4. 北京市畜牧总站, 北京 100107
2. The Key Laboratory of Farm Animal Genetic and Utilization of Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing 100193, China;
3. College of Animal Science and Technology, Jilin Agricultural University, Changchun 130118, China;
4. Animal Husbandry Station of Beijing, Beijing 100107, China
奶牛发情周期为18~24 d[1],如发情未被检出和及时配种,则需等到下一个情期才能再次发情配种,增加了非繁殖饲养成本[2],发情检测是奶牛繁殖管理重要环节[3]。目前,奶牛场大都使用人工观察、尾根涂蜡等传统发情检测方法,需要投入大量人力,难以适应规模化养殖需求[4]。尽管同期发情定时输精技术正在推广应用,奶牛发情检测仍发挥着不可替代的作用[5],并且向自动化发情检测方向发展[6-7],而奶牛发情生理参数变化是自动化发情监测技术研发的重要基础。
发情过程伴随体温[8]、活动量[9]、阴道黏液电阻值[10]等参数改变,探索这些参数的发情周期变化规律,对自动发情检测技术研发与提高尤为重要。传感技术的发展,使这些参数的不间断自动检测成为可能,其中,基于加速度传感技术的自动发情检测技术正在奶牛养殖业中推广应用[11-12]。但是,该技术发情期活动量监测阈值标准不一,有的报道活动量升高4.2倍,持续9.7 h[13],还有的升高2.9倍,持续13.8 h[14],并且发情检出效果报道差别较大,最高达94.7%[15],较低为85.6%[16],笔者前期研究表明,发情检出率为88.24%[10]。为提高该技术检出率,笔者进行发情奶牛活动量变化规律研究。
本研究使用自行研发的奶牛自动化发情检测系统(基于活动量变化设计研发,专利号:CN201520109489.6),对北京某奶牛场进行了发情母牛活动量检测,以探究发情奶牛的活动量变化规律。
1 材料与方法 1.1 试验动物及饲养管理本研究于2016年7月-2017年7月在北京市房山区某荷斯坦奶牛场进行。随机选择该场150头产后健康泌乳牛和7头青年母牛进行试验。泌乳牛每天6:00、13:00、19:00各挤奶1次,每次挤奶结束后,颈夹固定采食TMR日粮30 min后,自由采食,并提供清洁饮水。采用散栏饲养,圈舍中设有卧床,圈舍面积1 750 m2(50 m×35 m),每个圈舍约饲养36头奶牛,饲养密度48.6 m2·头-1,保证试验牛有足够的运动场地。
1.2 仪器设备及数据采集方法采用本团队自行研发的奶牛发情自动监测系统(简称:系统)进行活动量数据采集,该系统包括活动量检测装置(计步器)、数据采集装置(收集计步器采集到的活动量数据,并将其传输到上位机系统)、上位机系统(接收、存储和分析数据,并通过软件显示)。试验期间计步器佩戴于牛后腿跖部,数据采集装置安放在试验牛棚中。计步器每检测到牛腿部产生1次加速度记为1,每0.5 h进行1次累计。采集器每2.0 h收集1次计步器检测到的活动量数据,并上传至上位机处理分析和显示,实现活动量实时监测和发情预报功能。当活动量超过前11 d(间情期)同一时间点平均值2.5倍以上,持续2.5 h以上时,系统预报该牛发情,活动量降至2.5倍以下,则发情结束。
1.3 试验处理试验牛分为自然发情组和同期处理组,两组饲养方式相同,自然发情组不采用任何激素处理,试验牛从试验开始即进行活动量数据自动采集。同期处理组,对产后59 d未发情母牛,第60天颈部注射促性腺激素释放激素(gonadotropin-releasing hormone, GnRH),第67天颈部注射前列腺素F2α(prostaglandin F2α, PGF2α),第69天再次颈部注射GnRH。
1.4 发情鉴定系统自动进行发情母牛预报,所有系统预报发情母牛均由一位拥有8年从事奶牛繁殖工作经验的试验员进行直肠检查确认。本研究实际监测到128头发情牛的活动量变化,其中,春季监测出32头,夏季34头,秋季33头,冬季29头;青年牛检测出6头,1胎牛47头,2胎牛31头,≥3胎牛44头;根据发情类型,自然发情牛检测出113头,同期处理发情牛15头。
1.5 数据分析首先使用Excel表格对活动量数据进行初步整理,然后,使用SAS8.1混合模型Glm过程(一种基于最小二乘法原理,对不同数量样本数据进行显著性检验的方差分析方法),对不同胎次、不同季节和不同处理组母牛间发情活动量升高倍数和持续时间进行差异显著性检验。因为奶牛发情起始时间、结束时间不同,发情持续时间也不同,分析活动量变化规律时,分别将发情开始时间与发情结束时间定为0时,并分别向前后截取一段时间,进行活动量统计分析。
2 结果 2.1 奶牛发情时与发情前后活动量变化规律 2.1.1 发情当天与发情前后11 d活动量变化规律将发情当天(第0天)活动量总和与发情前后11 d每天活动量总和作图(图 1)。图 1显示,荷斯坦奶牛每天的活动量在发情前11 d和发情后11 d差异不显著(P>0.05),均在20 000左右,且显著低于发情当天活动量(P<0.05),其中发情后第1天活动量最低(20 672±2 192)。发情当天活动量最高(53 841±11 016),约为发情前后11 d每天活动量的2.69倍。
将128头荷斯坦奶牛发情开始前3 h和发情开始后7 h每0.5 h活动量平均值,与发情前11 d(间情期)同一时间点活动量平均值作图(图 2),图 2显示,奶牛处于间情期时,平均活动量较低,且很平稳,基本低于500。而发情期活动量显著升高(P<0.05),均高于1 500,其中发情开始后0.0~4.5 h时活动量呈上升趋势,至4.5 h时达最高(2 179±1 073)。将发情结束前7 h和发情结束后3 h每0.5 h活动量平均值与发情前11 d(间情期)同一时间点活动量平均值作图(图 3),图 3显示,从发情结束前3 h开始,活动量整体呈下降趋势,直至恢复到间情期活动量为止。
表 1显示,奶牛发情时活动量升高倍数在青年牛、1胎牛和2胎牛间差异不显著(P>0.05),而前三者发情活动量升高倍数显著高于3胎及3胎以上牛(P<0.05)(由于3胎以上各胎次牛数较少,故归纳在一起)。奶牛发情持续时间青年牛显著高于其他胎次牛(P<0.05),达(16.50±2.17) h,同时,随胎次增加,发情持续时间逐渐缩短,由1胎的(7.23±3.73) h缩减到3胎及3胎以上的(4.93±2.63) h。
参考文献[17],将3~5月定义为春季、6~8月为夏季、9~11月为秋季,12月~次年2月为冬季。对奶牛发情活动量的季节因素统计分析结果(表 2)表明,春、秋两季奶牛发情活动量平均升高6倍左右,差异不显著(P>0.05),但春(秋)、夏和冬3季节间差异显著(P<0.05),且活动量呈现随气温升高逐渐降低的趋势,其中冬季活动量升高倍数最大(8.10±4.64),夏季最低(5.58±2.94)。发情持续时间在春、冬两季间差异不显著(P>0.05),夏、秋两季间差异也不显著(P>0.05),而春(冬)季与夏(秋)季间差异显著(P<0.05),分别持续约7、5 h,即呈现出随夏季到来持续时间逐渐缩短,而随秋季向冬季转化逐渐延长的趋势。
对不同发情类型发情牛活动量数据的统计分析(表 3)表明,同期发情母牛发情活动量升高倍数(7.36±4.61)显著高于自然发情牛(6.70±4.00)(P<0.05),但发情持续时间差异不显著(P>0.05),均在6 h左右。
系统鉴定为发情的母牛均通过直肠检查进行验证,对系统鉴定的阳性发情牛和假阳性发情牛活动量升高倍数和持续时间的比较分析结果(表 4)显示,阳性发情牛活动量升高倍数和持续时间均显著高于假阳性发情牛(P<0.05)。
试验期间,共对121头自然发情泌乳牛、7头青年牛和29头同期处理泌乳牛进行了发情鉴定,直检验证结果(表 5)表明,本系统鉴定出107头自然发情牛(鉴定率88.43%)、15头同期发情牛(鉴定率51.72%)、6头青年牛(鉴定率85.71%)。
奶牛发情时大多呈现精神亢奋状态,表现坐立不安、嗅闻其他牛只阴户、爬跨或被其他牛只爬跨等现象,导致发情母牛运动量(活动量)升高[18]。该规律为自动化发情鉴定技术研发奠定了基础。研究表明,间情期时,母牛活动量变化不明显, 而发情时活动量可达间情期活动量的400%[19]。利用计步器每2 h记录一次活动量,共监测了57头奶牛发情活动量变化情况,发现发情时活动量比间情期增加(290±160)%[14]。还有研究发现,发情时比未发情时活动量约升高6倍[20]。本研究发现,奶牛发情前后11 d中,每天活动量差异不显著(P>0.05),均在20 000左右,而发情当天活动量(53 841)约为间情期的2.69倍。本研究发现,发情开始前和发情结束后母牛活动量与间情期平均活动量差异也不显著(P>0.05),均在500以下,而发情时活动量则可达间情期活动量的3.6~6.8倍,结果与Brehme等[20]研究结果相似,但略高于Kiddy[19]的研究结果。
奶牛活动量可能受季节、发情方式、产奶量、胎次和品种等因素影响[6]。有研究指出,经产母牛(≥1胎)发情活动量低于初产母牛[21],除第1、2胎母牛发情期运动量差异不显著外(P>0.05),其他各胎次间差异极显著(P < 0.01),且随胎次增加,运动量逐渐减少[22]。另有研究表明,初产牛和经产牛发情持续时间无明显变化[23]。本研究表明,青年牛、1胎和2胎奶牛发情活动量升高倍数差异不显著(P>0.05),而3胎及以上胎次奶牛发情活动量升高倍数显著低于青年牛~2胎奶牛(P<0.05),与孙保贵等[22]研究结果相似,而与Lopez等[23]不同。同时,发情持续时间也随胎次增加而逐渐缩短,且初产牛((16.50±2.17) h)显著高于经产牛(1胎:(7.23±3.73)h)(P<0.05)。可能是因为牛奶产量增高降低了血液循环中的雌激素水平[23],导致奶牛发情表现减弱。国外关于环境温度的研究表明,炎热季节发情活动量升高2.70倍,显著低于寒冷季节(2.90倍),而持续时间差异不大,均在8.5 h左右[24];另有研究指出,夏季奶牛发情活动量升高2.33倍,显著低于冬季(3.55倍),夏季发情持续时间((15.2±0.2) h)也显著低于冬季((17.3±0.9) h)[25]。本研究表明,夏季发情牛活动量和持续时间均显著低于冬季(P<0.05),与Sakatani等[25]报道的变化趋势相似,但活动量变化倍数均高于Ismael等[24]的结果;而持续时间却低于两人的研究结果。Sakatani等[25]研究还表明进行自然发情牛((16.8±0.9) h)与使用同期发情处理牛((14.1±1.2) h)发情持续时间差异不大,与本研究结果类似,但Sakatani等[25]报道的发情持续时间比本研究的发情持续时间(6.32±3.67) h和(5.97±2.63) h长,该差异可能由于他们的试验牛(日本黑牛)为不同于奶牛的肉牛品种。本研究还发现,系统预报发情阳性牛活动量升高倍数(6.72±4.10)和持续时间((6.27±3.54)h)均显著高于假阳性牛(P<0.05),假阳性预报结果,可能源于发情牛爬跨其他未发情牛或牛只调圈等造成的活动量突然升高。
国外有些计步器产品多在挤奶时上传数据,可能会造成预报滞后,影响及时配种,比如19:00挤完奶后,需到第2天6:00才挤奶,则第2天早上采集到的数据将是前1 d上传数据11 h后的数据,数据采集滞后可能使有些牛错过最佳配种时间[26]。同时,33%泌乳牛发情无明显行为变化,这些安静发情牛[27],对于计步器发情预报提出了更高要求,与其他研究相比,本研究每0.5 h累计一次活动量,可更细致、准确地分析奶牛发情开始与结束时间,每2 h上传一次活动量数据,可使繁殖员及时了解牛只发情情况,提早做准备,避免错过最佳输精时间。
本研究探究了荷斯坦奶牛发情时活动量变化规律,确定了活动量高于前11 d同一时间点活动量2.5倍,且持续2.5 h的自动预报阈值,该阈值从活动量增加倍数和持续时间两方面,兼顾了高度选育荷斯坦奶牛活动量增加较少、持续时间较短的品种特点,其对自然发情牛的检出率达到88.43%,假阳性率为9.86%,具有较高的检出率和较低的假阳性率,但还需进一步进行改进,使检出率更高,假阳性率更低。
4 结论本研究系统揭示了母牛发情活动量变化受胎次、季节和发情处理类型等的显著影响,确定活动量自动发情检测阈值时应兼顾上述因素,同时还要考虑到品种的影响,各个养殖企业应综合考虑上述因素后,根据本场具体情况设定一个与本场相适合的预报阈值。
[1] | JAISWAL R S, SINGH J, MARSHALL L, et al. Repeatability of 2-wave and 3-wave patterns of ovarian follicular development during the bovine estrous cycle[J]. Theriogenology, 2009, 72(1): 81–90. DOI: 10.1016/j.theriogenology.2009.02.014 |
[2] |
寇红祥, 赵福平, 任康, 等. 奶牛体温与活动量检测及变化规律研究进展[J]. 畜牧兽医学报, 2016, 47(7): 1306–1315.
KOU H X, ZHAO F P, REN K, et al. The progress on detection method and the regularities of body temperature and activities in dairy cows[J]. Acta Veterinaria et Zootechnica Sinica, 2016, 47(7): 1306–1315. (in Chinese) |
[3] | HELMER S D, BRITT J H. Mounting behavior as affected by stage of estrous cycle in Holstein heifers[J]. J Dairy Sci, 1985, 68(5): 1290–1296. DOI: 10.3168/jds.S0022-0302(85)80959-0 |
[4] | FIRK R, STAMER E, JUNGE W, et al. Automation of oestrus detection in dairy cows:A review[J]. Livest Prod Sci, 2002, 75(3): 219–232. DOI: 10.1016/S0301-6226(01)00323-2 |
[5] | FRICKE P M, CARVALHO P D, GIORDANO J O, et al. Expression and detection of estrus in dairy cows:The role of new technologies[J]. Animal, 2014, 8(S1): 134–143. DOI: 10.1017/S1751731114000299 |
[6] | LÓPEZ-GATIUS F, SANTOLARIA P, MUNDET I, et al. Walking activity at estrus and subsequent fertility in dairy cows[J]. Theriogenology, 2005, 63(5): 1419–1429. DOI: 10.1016/j.theriogenology.2004.07.007 |
[7] |
李小俊, 王振玲, 陈晓丽, 等. 奶牛体温变化规律及繁殖应用研究进展[J]. 畜牧兽医学报, 2016, 47(12): 2331–2341.
LI X J, WANG Z L, CHEN X L, et al. Study progress on the rule of body temperature and its application in reproduction of dairy cattle[J]. Acta Veterinaria et Zootechnica Sinica, 2016, 47(12): 2331–2341. (in Chinese) |
[8] | MIURA R, YOSHIOKA K, MIYAMOTO T, et al. Estrous detection by monitoring ventral tail base surface temperature using a wearable wireless sensor in cattle[J]. Anim Reprod Sci, 2017, 180: 50–57. DOI: 10.1016/j.anireprosci.2017.03.002 |
[9] | YÁNIZ J L, SANTOLARIA P, GIRIBET A, et al. Factors affecting walking activity at estrus during postpartum period and subsequent fertility in dairy cows[J]. Theriogenology, 2006, 66(8): 1943–1950. DOI: 10.1016/j.theriogenology.2006.05.013 |
[10] |
寇红祥, 李蓝祁, 王振玲, 等. 牛发情期活动量与阴道黏液电阻值变化规律的研究[J]. 畜牧兽医学报, 2017, 48(7): 1221–1228.
KOU H X, LI L Q, WANG Z L, et al. Study on the regulations of activity and vaginal electrical resistance of cattle during the estrous cycle[J]. Acta Veterinaria et Zootechnica Sinica, 2017, 48(7): 1221–1228. (in Chinese) |
[11] | ROELOFS J B, VAN EERDENBURG F J C M, SOEDE N M, et al. Pedometer readings for estrous detection and as predictor for time of ovulation in dairy cattle[J]. Theriogenology, 2005, 64(8): 1690–1703. DOI: 10.1016/j.theriogenology.2005.04.004 |
[12] |
蒋晓新, 刘炜, 魏星远, 等. 运用计步器鉴定泌乳盛期荷斯坦奶牛的发情效果研究[J]. 安徽农业科学, 2013, 41(15): 6728–6729, 6732.
JIANG X X, LIU W, WEI X Y, et al. Study on the effects of identifying the estrus of holstein cows during peak lactation by using pedometer[J]. Journal of Anhui Agricultural Sciences, 2013, 41(15): 6728–6729, 6732. DOI: 10.3969/j.issn.0517-6611.2013.15.058 (in Chinese) |
[13] | SAULS J A, VOELZ B E, HILL S L, et al. Increasing estrus expression in the lactating dairy cow[J]. J Dairy Sci, 2017, 100(1): 807–820. DOI: 10.3168/jds.2016-11519 |
[14] | SILPER B F, ROBLES I, MADUREIRA A M L, et al. Automated and visual measurements of estrous behavior and their sources of variation in Holstein heifers.Ⅰ:Walking activity and behavior frequency[J]. Theriogenology, 2015, 84(2): 312–320. DOI: 10.1016/j.theriogenology.2014.12.029 |
[15] |
李军, 罗艳红, 冯艳, 等. 陕西省千阳县奶牛计步器应用效果分析[J]. 中国乳业, 2016(1): 41–43.
LI J, LUO Y H, FENG Y, et al. Application effect analysis of Shanxi province of Qianyang by using pedometer[J]. China Dairy, 2016(1): 41–43. (in Chinese) |
[16] |
郑伟, 年景华, 李军. 奶牛计步器的应用效果分析[J]. 中国奶牛, 2014(22): 32–34.
ZHENG W, NIAN J H, LI J. Analysis on the application effect of cow pedometer[J]. China Dairy Cattle, 2014(22): 32–34. DOI: 10.3969/j.issn.1004-4264.2014.22.009 (in Chinese) |
[17] |
蒋晓新, 邓双义, 魏星远, 等. 北方地区季节因素对荷斯坦奶牛步履数和奶牛泌乳性能的影响[J]. 江苏农业科学, 2013, 41(11): 225–227.
JIANG X X, DENG S Y, WEI X Y, et al. Effects of seasonal factors on the step number of Holstein cows and milk performance of dairy cows in the Northern region[J]. Jiangsu Agricultural Sciences, 2013, 41(11): 225–227. DOI: 10.3969/j.issn.1002-1302.2013.11.087 (in Chinese) |
[18] | VAN EERDENBURG F J, LOEFFLER H S, VAN VLIET J H. Detection of oestrus in dairy cows:A new approach to an old problem[J]. Vet Quart, 1996, 18(2): 52–54. DOI: 10.1080/01652176.1996.9694615 |
[19] | KIDDY C A. Variation in physical activity as an indication of estrus in dairy cows[J]. J Dairy Sci, 1977, 60(2): 235–243. DOI: 10.3168/jds.S0022-0302(77)83859-9 |
[20] | BREHME U, STOLLBERG U, HOLZ R, et al. ALT pedometer-New sensor-aided measurement system for improvement in oestrus detection[J]. Comput Electron Agric, 2008, 62(1): 73–80. DOI: 10.1016/j.compag.2007.08.014 |
[21] | MADUREIRA A M L, SILPER B F, BURNETT T A, et al. Factors affecting expression of estrus measured by activity monitors and conception risk of lactating dairy cows[J]. J Dairy Sci, 2015, 98(10): 7003–7014. DOI: 10.3168/jds.2015-9672 |
[22] |
孙保贵, 郭予伟, 陈茂学, 等. 奶牛运动量辅助发情诊断参数的研究[J]. 中国畜牧杂志, 2015, 51(13): 71–75.
SUN B G, GUO Y W, CHEN M X, et al. Study on estrus detection parameter based on activity record in lactating cows[J]. Chinese Journal of Animal Science, 2015, 51(13): 71–75. DOI: 10.3969/j.issn.0258-7033.2015.13.016 (in Chinese) |
[23] | LOPEZ H, SATTER L D, WILTBANK M C. Relationship between level of milk production and estrous behavior of lactating dairy cows[J]. Anim Reprod Sci, 2004, 81(3-4): 209–223. DOI: 10.1016/j.anireprosci.2003.10.009 |
[24] | ISMAEL A, STRANDBERG E, BERGLUND B, et al. Seasonality of fertility measured by physical activity traits in Holstein cows[J]. J Dairy Sci, 2016, 99(4): 2837–2848. DOI: 10.3168/jds.2015-10067 |
[25] | SAKATANI M, TAKAHASHI M, TAKENOUCHI N. The efficiency of vaginal temperature measurement for detection of estrus in Japanese Black cows[J]. J Reprod Dev, 2016, 62(2): 201–207. DOI: 10.1262/jrd.2015-095 |
[26] |
陆静. 第二代计步器的功能[C]//中国奶业协会年会. 杭州, 2009: 2.
LU J. The function of the second generation of pedometer[C]//China Dairy Association Annual Conference. Hangzhou, 2009: 2. (in Chinese) |
[27] | TALUKDER S, KERRISK K L, INGENHOFF L, et al. Infrared technology for estrus detection and as a predictor of time of ovulation in dairy cows in a pasture-based system[J]. Theriogenology, 2014, 81(7): 925–935. DOI: 10.1016/j.theriogenology.2014.01.009 |