动物营养学报    2022, Vol. 34 Issue (2): 1343-1351    PDF    
基于Meta分析的母羊泌乳模型构建及泌乳性能预测
李金辉 , 陈佳欣 , 王里彦 , 杨彩虹 , 田星哲 , 田沛知 , 严慧 , 刘月琴 , 纪守坤 , 张英杰     
河北农业大学动物科技学院, 保定 071000
摘要: 母羊泌乳量不仅关乎母羊自身营养需求,也是影响羔羊生长发育的重要因素,准确预测母羊泌乳量在母羊和羔羊精准饲养中具有现实意义。本文利用方差倒数Meta分析方法研究母羊泌乳曲线模型及泌乳参数,并对其泌乳性能进行预测。分别在Web of Science、CNKI等文献数据库中检索母羊泌乳曲线及影响因素的相关文献,依据纳入标准筛选和提取泌乳曲线相关参数,采用方差倒数Meta分析方法进行合并,综合评估母羊泌乳曲线及其影响因素。共纳入文献22篇,提取数据98条。结果显示:山羊泌乳曲线模型为y=0.943 8t0.359 9e(-0.008 5t,绵羊为y=0.677 2t0.349 4e(-0.014 0t,乳用山羊为y=1.143 7t0.318 7e(-0.008 0t,非乳用山羊为y=0.419 0t0.468 1e(-0.011 8t,乳用绵羊为y=0.796 8t0.256 2e(-0.009 0t,非乳用绵羊为y=0.570 0t0.441 0e(-0.026 0t,其中y为预测泌乳量(kg/d),t为泌乳天数(d),e为自然数(2.718 3),模型参数具有较高的稳健性(敏感性 < 5%)。山羊、绵羊、乳用山羊、非乳用山羊、乳用绵羊和非乳用绵羊初乳产量(产后第1天泌乳量)分别为0.94、0.67、1.13、0.41、0.79和0.55 kg,泌乳峰值时间分别在泌乳期第42.34天、第25.00天、第39.84天、第39.67天、第28.47天和第16.96天,峰值泌乳量分别为2.54、1.47、2.69、1.47、1.45和1.28 kg,泌乳持续力分别为6.48、5.76、6.37、6.52、5.92和5.26。本研究获得了山羊、绵羊、乳用山羊、非乳用山羊、乳用绵羊和非乳用绵羊泌乳量预测模型及其泌乳参数,为泌乳母羊和哺乳期羔羊精准饲养提供了基础数据。
关键词: 泌乳曲线    Wood模型    Meta分析    泌乳参数    
Establishment of Lactation Model Based on Meta-Analysis and Lactation Performance Prediction
LI Jinhui , CHEN Jiaxin , WANG Liyan , YANG Caihong , TIAN Xingzhe , TIAN Peizhi , YAN Hui , LIU Yueqin , JI Shoukun , ZHANG Yingjie     
College of Animal Science and Technology, Hebei Agricultural University, Baoding 071000, China
Abstract: The milk yield of ewes is not only related to the nutritional requirement of ewes, but also an important factor affecting the growth and development of lambs. It is of great importance to accurately predict the milk yield of ewes in the balance nutrition of ewes and lambs. In this paper, Meta-analysis method was used to study the lactation curve model of ewes and lactation parameters, thus to predict the lactation performance. The related literatures related to ewe lactation curve were searched in the literature databases such as Web of Science and CNKI, and the parameters of lactation curve were screened and extracted according to the inclusion criteria. The Meta-analysis method was used to combine and comprehensively evaluate the lactation curve and its influencing factors of ewes. A total of 22 literatures were included and 87 groups of data were extracted. The results showed that the lactation curve model for goats was y=0.943 8t0.359 9e(-0.008 5t), for sheep was y=0.677 2t0.349 4e(-0.014 0t), for dairy goats was y=1.143 7t0.318 7e(-0.008 0t), for non-dairy goats was y=0.419 0t0.468 1e(-0.011 8t), for dairy sheep was y=0.796 8t0.256 2e(-0.009 0t) and for non-dairy sheep was y=0.570 0t0.441 0e(-0.026 0t). In the lactation curve models, y was y predicted milk yield (kg/d), t was days of milk (d), e was natural number (2.718 3). The model parameters was with high robustness (sensitivity < 5%). The colostrum yield (milk yield at day 1 of postpartum) of goats, sheep, dairy goats, non-dairy goats, dairy sheep and non-dairy sheep was 0.94, 0.67, 1.13, 0.41, 0.79 and 0.55 kg, respectively; the peak lactation time was days 42.34, 25.00, 39.84, 39.67, 28.47 and 16.96, respectively; the peak milk yield was 2.54, 1.47, 2.69, 1.47, 1.45 and 1.28 kg, respectively; and the lactation persistence was 6.48, 5.76, 6.37, 6.52, 5.92 and 5.26, respectively. The lactation curves of Wood models in goat, sheep, dairy goats, non-dairy goats, dairy sheep and non-dairy sheep are achieved based on Meta-analysis, and the lactation parameters of different types of ewes are obtained in this study, which provide basic information for nutrition balance in lactating ewes and pre-weaning lambs.
Key words: lactation curve    Wood model    Meta-analysis    lactation parameters    

羊奶不仅营养丰富、吸收效率高,还具有一定的保健价值,正逐渐得到人们普遍认可[1-2]。对于母羊,乳产量是泌乳期营养需求的重要参数;对于哺乳羔羊,充足的母乳是其正常发育、增强抵抗力和提高存活率的关键[3-4]。准确预测母羊的泌乳性能是进行精准化饲养管理,保证母羊泌乳性能和羔羊生长性能的关键措施。

母羊泌乳性能可以通过2种方式获得:一是统计单位时间内母羊实际泌乳量,此方法对本品种母羊预测更精准,但存在投入高、时效性差的缺点,因此难以在生产中广泛应用[5];二是通过数学模型预测母羊泌乳性能,现有模型主要有Wood模型、Wilmink模型、逆多项模型、Ali-Schaeffer模型等[6],其中Wood模型拟合效果好、可解释度强,因此在羊泌乳性能预测中应用最多。Wood模型可表述为:y=atbe(-ct),其中参数a表示母羊在产羔后第1天的泌乳量,参数b表示泌乳高峰出现前泌乳曲线的上升速率,参数c表示泌乳高峰出现后泌乳曲线的下降速率[7]。采用Wood模型拟合泌乳曲线预测牛、羊泌乳性能得到了广泛认可,已成为泌乳量估测的重要手段[8],但母羊泌乳量受遗传、环境、年龄、品种、胎次、挤奶方式、产羔数、生理状态等因素的影响[9],其中品种是影响母羊泌乳性能的最主要因素[10]。由于山羊和绵羊品种繁多,不同研究间所获得泌乳参数差异较大,人们对母羊泌乳性能的总体规律性以及不同类型母羊的泌乳性能了解有限。

Meta分析是一种对前人文献数据汇总再分析的方法,具有样本量大、规律性强、可信度高等优点,是目前常见的数据统计分析方法,在畜牧学领域应用广泛[11-12]。本研究检索了国内外发表的关于母羊泌乳曲线Wood模型分析的相关文献,从所查阅的文献中提取有关泌乳曲线Wood模型参数a、b、c的数值进行Meta分析,探究各类型母羊泌乳的总体规律,以期为泌乳母羊的饲养与管理提供参考。

1 材料与方法 1.1 文献检索

使用计算机检索以及手工检索,分别以表 1中的9个检索策略作为题名、主题或者关键词,分别在Web of Science、PubMed、维普中文期刊、CNKI、万方等数据库进行检索,检索时间区间为1980年1月至2021年5月。搜集不同品种类型母羊的Wood模型参数。

表 1 文献检索策略 Table 1 Literature search strategies
1.2 纳入标准

本研究纳入标准为:1)采用Wood模型研究母羊泌乳曲线的试验;2)明确了羊的品种和试验处理方式;3)主要统计指标包含泌乳曲线分析中Wood模型参数a、b、c及其标准差(SD)或者标准误(SE)以及有效例数。

1.3 排除标准

本研究排除标准为:1)没有明确母羊品种的文章;2)没有明确指出有效例数的文章;3)母羊患有疾病且明显影响泌乳的文章;4)非试验分析得出,属于综述类的文章;5)采用非Wood模型的文章。

1.4 文献筛选

首先阅读论文题目和摘要,排除不符合纳入标准的文献,对于符合纳入标准的文献进行全文阅读,排除数据不全、采用非Wood模型和重复检索到的文献。

1.5 文献数据提取

利用事先设计好的Excel表格进行文献数据信息提取与整理。提取的信息包括发表时间,作者,Wood模型参数a、b、c数值及其SD或者SE,品种,试验动物重复数(n)以及其他影响因素(胎次、产羔数、羔羊性别、羔羊初生重、产羔月份与季节、母羊年龄、营养状况等)。采用SE进行后续数据分析,对于没有SE的文献使用公式SE=SD/sqrt(n)进行SD和SE的转换。

1.6 产出指标

母羊Wood模型参数a、b、c数值为本次Meta分析关注的统计指标。影响因素为母羊品种,由于相同品种母羊Wood模型相关文献较少,难以满足Meta分析要求,因此本研究将母羊分为山羊和绵羊、乳用山羊和非乳用山羊、乳用绵羊和非乳用绵羊。

1.7 泌乳参数

基于Wood模型的母羊泌乳曲线方程[7]

式中:y表示泌乳母羊产羔后第t天的泌乳量;t表示泌乳母羊产羔后的泌乳天数;a、b、c为常数项;e为自然数(2.718 3)。

通过Wood模型获得泌乳峰值时间Dmax、峰值泌乳量Ymax以及泌乳持续力P,预测方程分别为:

1.8 Meta分析

利用R-4.0.0软件进行方差倒数Meta分析,采用加权合并效应量方法计算其合并后的点估计值及95%可信区间。利用I2P值来对统计学异质性进行判断:当I2 < 50%,P>0.10时,认为无异质性,采用固定效应模型进行分析;当I2>50%,P<0.10时,则认为存在统计学异质性,采用随机效应模型进行分析[13]。本研究纳入文献大多数存在异质性,因此采用随机效应模型进行分析。数据预测值以95%置信区间(95%CI)表示,P<0.05认为具有统计学意义。使用Egger检验进行文献发表偏倚评估,P<0.10认为存在发表偏倚。采用单项文献逐项剔除方法进行敏感性分析,若剔除单项文献后合并效应值变化小于5%,则认为敏感性低,合并效应值具有稳健性。

2 结果与分析 2.1 文献检索结果

本研究共检索到相关文献153篇,进行初筛(阅读文献题目和摘要)后,筛除不符合纳入标准的文献98篇,余下55篇;对剩余的文献进行全文通读复筛,进一步排除不符合要求的文献33篇,最终获得符合纳入标准的文献22篇,其中中文文献9篇,英文文献13篇。纳入文献的基本信息见表 2

表 2 纳入文献的基本信息 Table 2 Basic information of included literatures
2.2 山羊和绵羊Wood模型Meta分析

本研究纳入山羊的文献14篇,数据81条,绵羊的文献8篇,数据17条,Meta分析结果见表 3。结果表明,山羊各研究间的Wood模型参数a、b均存在统计学异质性(I2>50%,P<0.10),参数c无异质性(I2 < 50%,P<0.10),采用随机效应模型合并效应值,获得山羊的泌乳模型为y=0.943 8t0.359 9e(-0.008 5t)。绵羊各研究间的Wood模型参数a、b、c均存在统计学异质性(I2>50%,P<0.10),采用随机效应模型合并效应值,获得绵羊的泌乳模型为y=0.677 2t0.349 4e(-0.014 0t)。Egger检验表明,山羊a、b、c值和绵羊b值存在发表偏倚(P<0.10),绵羊a、c值无发表偏倚(P>0.10);敏感性分析表明,去除任一文献后,合并效应值变化小于5%,表明本结果具有稳健性。

表 3 山羊和绵羊Wood模型Meta分析结果 Table 3 Meta-analysis results of Wood model in goats and sheep
2.3 乳用山羊和非乳用山羊Wood模型Meta分析

本研究纳入乳用山羊的文献9篇,数据59条,非乳用山羊的文献5篇,数据22条,Meta分析结果见表 4。结果表明,乳用山羊各研究间的Wood模型参数a、b、c均存在统计学异质性(I2>50%,P<0.10),采用随机效应模型合并效应值,获得乳用山羊的泌乳模型为y=1.143 7t0.318 7e(-0.008 0t)。非乳用山羊各研究间的Wood模型参数a、b均存在统计学异质性(I2>50%,P<0.10),参数c无异质性(I2 < 50%,P>0.10),采用随机效应模型合并效应值,获得非乳用山羊的泌乳模型为y=0.419 0t0.468 1e(-0.011 8t)。Egger检验表明,乳用山羊a、b、c值和非乳用山羊b、c值存在发表偏倚(P<0.10),非乳用山羊a值无发表偏倚(P>0.10);敏感性分析表明,去除任一文献后,合并效应值变化小于5%,表明本结果具有稳健性。

表 4 乳用山羊和非乳用山羊Wood模型Meta分析结果 Table 4 Meta-analysis results of Wood model in dairy goats and non-dairy goats
2.4 乳用绵羊和非乳用绵羊Wood模型Meta分析

本研究共纳入乳用绵羊的文献4篇,数据8条,非乳用绵羊的文献4篇,数据9条,Meta分析结果见表 5。结果表明,乳用绵羊各研究间的Wood模型参数a、b、c均存在统计学异质性(I2>50%,P<0.10),采用随机效应模型合并效应值,获得乳用绵羊的泌乳模型为y=0.796 8t0.256 2e(-0.009 0t)。非乳用绵羊各研究间的Wood模型参数a、b、c存在统计学异质性(I2>50%,P<0.10),采用随机效应模型合并效应值,获得非乳用绵羊的泌乳模型为y=0.570 0t0.441 0e(-0.026 0t)。Egger检验表明,非乳用绵羊b值存在发表偏倚(P<0.10),乳用绵羊a、b、c值和非乳用绵羊a、c值无发表偏倚(P>0.10);敏感性分析表明,去除任一文献后,合并效应值变化小于5%,表明本结果具有稳健性。

表 5 乳用绵羊和非乳用绵羊Wood模型Meta分析结果 Table 5 Meta-analysis results of Wood model in dairy sheep and non-dairy sheep
2.5 不同类型母羊泌乳曲线及泌乳参数

采用Wood模型对山羊和绵羊、乳用山羊和非乳用山羊、乳用绵羊和非乳用绵羊的泌乳曲线进行拟合(图 1),结果显示,山羊、绵羊、乳用山羊、非乳用山羊、乳用绵羊和非乳用绵羊初乳产量(产后第1天泌乳量)分别为0.94、0.67、1.13、0.41、0.79和0.55 kg,泌乳峰值时间分别在泌乳期第42.34天、第25.00天、第39.84天、第39.67天、第28.47天和第16.96天,峰值泌乳量分别为2.54、1.47、2.69、1.47、1.45和1.28 kg,泌乳持续力分别为6.48、5.76、6.37、6.52、5.92和5.26。

a:山羊和绵羊泌乳曲线;b:乳用山羊和非乳用山羊泌乳曲线;c:乳用绵羊和非乳用绵羊泌乳曲线。 a: the lactation curves of goats and sheep; b: the lactation curves of dairy goats and non-dairy goats; c: the lactation curves of dairy sheep and non-dairy sheep. 图 1 不同类型母羊泌乳曲线 Fig. 1 Lactation curves of different types of ewes
3 讨论 3.1 不同类型母羊泌乳参数

初乳指母羊产后前3 d分泌的母乳,生产上以母羊产后第1天的母乳最为关键,是新生羔羊获得营养和免疫的主要途径,保障新生羔羊初乳摄入量对羔羊哺乳期健康和生长性能具有决定性作用[14-15]。Alves等[16]研究表明,新生羔羊在产后24 h内摄入至少30 g免疫球蛋白G(IgG)(约500 g初乳)才能保证被动免疫的成功建立。本研究结果显示山羊、绵羊、乳用山羊、非乳用山羊、乳用绵羊和非乳用绵羊初乳产量(产后第1天泌乳量)分别为0.94、0.67、1.13、0.41、0.79和0.55 kg。对于单羔,除了非乳用山羊外,其余类型母羊初乳量足够新生羔羊成功建立被动免疫;对于多羔,除乳用山羊外,各类型母羊初乳产量均不足。这一结果与前人研究结果一致,即单羔羔羊被动免疫成功率显著高于双羔和多羔羔羊[17-18]。因此,生产上母羊尤其是多羔母羊初乳产量和羔羊初乳采食量应引起重视,以保障新生羔羊被动免疫的成功建立。

泌乳峰值时间指母羊在整个泌乳期日泌乳量最高的时间,是衡量母羊泌乳性能的重要依据。泌乳峰值时间不仅与泌乳曲线上升和下降速率有关,还在一定程度上反映泌乳潜力的高低[19]。母羊泌乳峰值时间到来之前,泌乳量持续上升可满足羔羊随体重增加而增加的营养需要;泌乳峰值时间过后泌乳量减少,仅依赖母乳不能满足羔羊生长发育的需求,羔羊营养需求从依赖母乳过渡到饲粮。本研究结果显示山羊、绵羊、乳用山羊、非乳用山羊、乳用绵羊和非乳用绵羊泌乳峰值时间分别出现在泌乳期第42.34天、第25.00天、第39.84天、第39.67天、第28.47天和第16.96天。不同类型母羊泌乳峰值时间差别较大,因此需要针对不同类型母羊制定合理补饲时间,保障羔羊营养需求并降低断奶应激[20]

峰值泌乳量指母羊到达泌乳峰值时间的泌乳量,反映母羊泌乳潜力的高低;泌乳持续力指泌乳峰值时间过后母羊泌乳量下降的缓急程度,反映母羊持续泌乳能力的强弱,两者均与总泌乳量相关[21]。品种、产羔数、胎次等因素均会影响峰值泌乳量和泌乳持续力。Elvira等[22]研究了Toggenbourg山羊和Anglo-Nubian山羊的泌乳性能,结果显示Toggenbourg山羊的泌乳持续力更强;Rojo-Rubio等[23]研究发现产双羔母羊比产单羔母羊泌乳峰值时间更早,峰值泌乳量和总泌乳量更高;Carnicella等[24]研究指出,母羊第3、4胎次比第1、2胎次峰值泌乳量和总泌乳量更高;刘长海[25]研究表明,品种是影响母羊峰值泌乳量和泌乳持续力的主要因素。本研究结果显示山羊、绵羊、乳用山羊、非乳用山羊、乳用绵羊和非乳用绵羊峰值泌乳量分别为2.54、1.47、2.69、1.47、1.45和1.28 kg,泌乳持续力分别为6.48、5.76、6.37、6.52、5.92和5.26。不同类型羊泌乳能力差异较大,总体表现为绵羊泌乳能力低于山羊,非乳用羊泌乳能力低于乳用羊。本研究结果提示,需要针对不同类型母羊的泌乳性能为母羊搭配平衡营养,为羔羊提供合理代乳粉补充,实现母羊和羔羊精准饲养。

3.2 异质性分析

本研究通过异质性分析表明山羊和非乳用山羊a、b值存在异质性,其异质性主要来自Montaldo等[26]的研究,由于Montaldo等[26]的研究包括基因、季节和年龄等多种影响因素,这些因素均会对泌乳参数产生影响,可能是该研究异质性较大的主要原因;绵羊和乳用绵羊a、b、c值的异质性主要来自Portolano等[27]的文献,Portolano等[27]研究的Comisana绵羊是著名的乳用品种,其泌乳性能甚至与乳用山羊相当,使得a、b、c值与Meta分析结果差别较大;乳用山羊a、b、c值的异质性主要来自Shaat[28]的文献,Shaat[28]选择在产后第2周开始测量Zaraibi山羊的泌乳量,由于时间间隔太长导致测量结果与其他研究差异较大,可能是造成数据异质性的主要原因;非乳用绵羊a、b、c值的异质性主要来自张力等[29]的文献,与大部分研究采用直接测定母羊泌乳量不同,张力等[29]通过羔羊哺乳前后体重差法研究湖羊泌乳量,这种方法学差异可能是导致其泌乳量测定结果出现异质性的原因。为降低异质性的影响,本研究对山羊和绵羊分别根据乳用和非乳用群体进行了亚组分析,有效降低了组内异质性,但由于目前采用Wood模型拟合母羊泌乳性能的相关文献数量相对较少,难以进一步进行亚组分析将异质性去除,随今后采用Wood模型预测母羊泌乳性能文献数量的增多,可通过不断扩大数据量来提高本研究Meta分析结果的准确性。

3.3 研究局限性

前人研究表明,品种、遗传、泌乳期、饲养管理、季节、营养水平、胎次以及环境等因素均会对母羊的泌乳性能产生影响,因此会一定程度改变泌乳模型参数[28, 30]。目前对于相同品种、胎次、季节等因素条件下泌乳模型的研究相对较少,导致不同研究中泌乳模型参数差异较大,进而造成了本研究数据的统计学异质性。本研究纳入文献遵循严格的纳入和排除标准,数据统计具有较高稳健性,在预测不同类型母羊泌乳量时具有较高普适性,可为不同类型母羊的泌乳量预测提供基础参数。

4 结论

本研究通过方差倒数Meta分析,得到不同类型母羊泌乳量预测模型:山羊为y=0.943 8t0.359 9e(-0.008 5t),绵羊为y=0.677 2t0.349 4e(-0.014 0t),乳用山羊为y=1.143 7t0.318 7e(-0.008 0t),非乳用山羊为y=0.419 0t0.468 1e(-0.011 8t),乳用绵羊为y=0.796 8t0.256 2e(-0.009 0t),非乳用绵羊为y=0.570 0t0.441 0e(-0.026 0t),其中y为预测泌乳量(kg/d);t为泌乳天数(d),e为自然数(2.718 3)。本研究获得的不同类型母羊的泌乳参数可为泌乳期母羊和哺乳期羔羊精准饲养提供基础参数。

参考文献
[1]
TSIPLAKOU E, CHADIO S, PAPADOMICHELAKIS G, et al. The effect of long term under- and over-feeding on milk and plasma fatty acids profile and on insulin and leptin concentrations of goats[J]. International Dairy Journal, 2012, 24(2): 87-92. DOI:10.1016/j.idairyj.2011.05.010
[2]
PAZZOLA M, CIPOLAT-GOTET C, BITTANTE G, et al. Phenotypic and genetic relationships between indicators of the mammary gland health status and milk composition, coagulation, and curd firming in dairy sheep[J]. Journal of Dairy Science, 2018, 101(4): 3164-3175. DOI:10.3168/jds.2017-13975
[3]
张磊, 乔为盛, 白前前, 等. 湖羊胎产羔数对羔羊生长发育的影响[J]. 家畜生态学报, 2020, 41(4): 43-46.
ZHANG L, QIAO W S, BAI Q Q, et al. Effects of lambing number on the growth and development of lambs in Hu sheep[J]. Acta Ecologae Animalis Domastici, 2020, 41(4): 43-46 (in Chinese). DOI:10.3969/j.issn.1673-1182.2020.04.008
[4]
RAZZAQUE M A, AL-MUTAWA S, ABBAS S, et al. Performance of pre-weaned dairy calves under hot arid environment: effects of immunoglobulins and age on diseases and mortality[J]. American Journal of Applied Sciences, 2009, 6(11): 1885-1891. DOI:10.3844/ajassp.2009.1885.1891
[5]
李讨讨, 马友记. 湖寒杂交后代母羊产后1~30 d泌乳量及乳成分变化规律研究[J]. 中国草食动物科学, 2017, 37(3): 8-11.
LI T T, MA Y J. Changes of the milk yield and milk composition in confined 1-30 days postpartum Hu sheep and small-tail Han sheep crossbred ewes[J]. China Herbivore Science, 2017, 37(3): 8-11 (in Chinese). DOI:10.3969/j.issn.2095-3887.2017.03.003
[6]
毛永江, 张美荣, 许兆君, 等. 用Wood模型拟合南方中国荷斯坦牛产奶量、乳脂率、乳蛋白率和体细胞评分及其效果分析[J]. 畜牧兽医学报, 2012, 43(12): 1904-1909.
MAO Y J, ZHANG M R, XU Z J, et al. Lactation curve fittings of wood's nonlinear model for milk yield, milk fat, milk protein and somatic cell score for Chinese Holstein in southern China[J]. Acta Veterinaria et Zootechnica Sinica, 2012, 43(12): 1904-1909 (in Chinese).
[7]
WOOD P D P. Algebriac model of the lactation curve in cattle[J]. Nature, 1967, 216: 164-165. DOI:10.1038/216164a0
[8]
黄耀辉, 王永军, 胡秀芝. 西农萨能奶山羊泌乳规律研究[J]. 家畜生态学报, 2015, 36(8): 11-16, 65.
HUANG Y H, WANG Y J, HU X Z. Study on lactation rules of Xinong Saanen dairy goats[J]. Acta Ecologae Animalis Domastici, 2015, 36(8): 11-16, 65 (in Chinese). DOI:10.3969/j.issn.1673-1182.2015.08.003
[9]
CAPPIO-BORLINO A, PULINA G, ROSSI G. A non-linear modification of Wood's equation fitted to lactation curves of Sardinian dairy ewes[J]. Small Ruminant Research, 1995, 18(1): 75-79. DOI:10.1016/0921-4488(95)00713-U
[10]
DIJKSTRA J, LOPEZ S, BANNINK A, et al. Evaluation of a mechanistic lactation model using cow, goat and sheep data[J]. The Journal of Agricultural Science, 2010, 148(3): 249-262. DOI:10.1017/S0021859609990578
[11]
张千, 罗静, 李飞, 等. 断奶应激影响犊牛免疫功能的Meta分析[J]. 动物营养学报, 2019, 31(1): 175-187.
ZHANG Q, LUO J, LI F, et al. Meta-analysis of effect of weaning stress on immunity of calves[J]. Chinese Journal of Animal Nutrition, 2019, 31(1): 175-187 (in Chinese). DOI:10.3969/j.issn.1006-267x.2019.01.023
[12]
徐博成, 李智, 汪以真, 等. 抗菌肽对仔猪生长性能、腹泻率和免疫球蛋白水平影响的Meta分析[J]. 动物营养学报, 2020, 32(8): 3584-3593.
XU B C, LI Z, WANG Y Z, et al. Effects of antimicrobial peptides on growth performance, diarrhea rate and immunoglobulin levels of piglets: a Meta-analysis[J]. Chinese Journal of Animal Nutrition, 2020, 32(8): 3584-3593 (in Chinese). DOI:10.3969/j.issn.1006-267x.2020.08.017
[13]
董书伟, 张世栋, 王东升, 等. 中药治疗牛子宫内膜炎的系统评价和Meta分析[J]. 畜牧兽医学报, 2015, 46(11): 2085-2096.
DONG S W, ZHANG S D, WANG D S, et al. Efficacy of Chinese medicine on endometritis in cow: a systematic review and Meta-analysis[J]. Acta Veterinaria et Zootechnica Sinica, 2015, 46(11): 2085-2096 (in Chinese). DOI:10.11843/j.issn.0366-6964.2015.11.023
[14]
龚会娜, 杨永林, 杨华, 等. 羊初乳的功能性成分及作用[J]. 当代畜牧, 2019(3): 33-35.
GONG H N, YANG Y L, YANG H, et al. Research progress on nutritional compositions and function of sheep colostrums[J]. Contemporary Animal Husbandry, 2019(3): 33-35 (in Chinese).
[15]
梁水源, 木其尔, 格日勒图, 等. 反刍幼畜人工初乳的研究进展[J]. 动物营养学报, 2020, 32(9): 4037-4043.
LIANG S Y, MU Q E, GE R L T, et al. Research progress of young ruminant artificial colostrum[J]. Chinese Journal of Animal Nutrition, 2020, 32(9): 4037-4043 (in Chinese).
[16]
ALVES A C, ALVES N G, ASCARI I J, et al. Colostrum composition of Santa Inês sheep and passive transfer of immunity to lambs[J]. Journal of Dairy Science, 2015, 98(6): 3706-3716.
[17]
TURQUINO C F, FLAIBAN K K M C, LISBÔA J A N, et al. Transferência de imunidade passiva em cordeiros de corte manejados extensivamente em clima tropical[J]. Pesquisa Veterinária Brasileira, 2011, 31(3): 199-205. DOI:10.1590/S0100-736X2011000300003
[18]
CHRISTLEY R M, MORGAN K L, PARKIN T D H, et al. Factors related to the risk of neonatal mortality, birth-weight and serum immunoglobulin concentration in lambs in the UK[J]. Preventive Veterinary Medicine, 2003, 57(4): 209-226. DOI:10.1016/S0167-5877(02)00235-0
[19]
亓建刚, 杨湛澄, 李想, 等. 产犊季节对北京地区荷斯坦奶牛Wood泌乳曲线参数的影响探究[J]. 中国畜牧兽医, 2018, 45(11): 3167-3175.
QI J G, YANG Z C, LI X, et al. Effects of calving season on the parameters of wood lactation curve of Holstein dairy cow in Beijing area[J]. China Animal Husbandry & Veterinary Medicine, 2018, 45(11): 3167-3175 (in Chinese).
[20]
马清国. 高寒地区舍饲羔羊标准化培育技术[J]. 甘肃畜牧兽医, 2021, 51(6): 63-65.
MA Q G. Standardized cultivation techniques of house-fed lambs in alpine regions[J]. Gansu Animal and Veterinary Sciences, 2021, 51(6): 63-65 (in Chinese). DOI:10.3969/j.issn.1006-799X.2021.06.020
[21]
罗清尧, 熊本海, 马毅, 等. 中国荷斯坦奶牛第二泌乳期泌乳曲线模型的研究[J]. 中国农业科学, 2010, 43(23): 4910-4916.
LUO Q Y, XIONG B H, MA Y, et al. Study on lactation curve models of Chinese Holstein for the second parity[J]. Scientia Agricultura Sinica, 2010, 43(23): 4910-4916 (in Chinese). DOI:10.3864/j.issn.0578-1752.2010.23.017
[22]
ELVIRA L, HERNANDEZ F, CUESTA P, et al. Accurate mathematical models to describe the lactation curve of Lacaune dairy sheep under intensive management[J]. Animal, 2013, 7(6): 1044-1052. DOI:10.1017/S175173111200239X
[23]
ROJO-RUBIO R, KHOLIF A E, SALEM A Z M, et al. Lactation curves and body weight changes of Alpine, Saanen and anglo-nubian goats as well as pre-weaning growth of their kids[J]. Journal of Applied Animal Research, 2016, 44(1): 331-337. DOI:10.1080/09712119.2015.1031790
[24]
CARNICELLA D, DARIO M, AYRES M C C, et al. The effect of diet, parity, year and number of kids on milk yield and milk composition in Maltese goat[J]. Small Ruminant Research, 2008, 77(1): 71-74. DOI:10.1016/j.smallrumres.2008.02.006
[25]
刘长海. 母羊泌乳性能的影响因素和提高方法[J]. 现代畜牧科技, 2019(11): 38-39.
LIU C H. Influencing factors and improving methods of lactation performance of ewes[J]. Modern Animal Husbandry Science & Technology, 2019(11): 38-39 (in Chinese).
[26]
MONTALDO H, ALMANZA A, JUÁREZ A. Genetic group, age and season effects on lactation curve shape in goats[J]. Small Ruminant Research, 1997, 24(3): 195-202. DOI:10.1016/S0921-4488(96)00946-7
[27]
PORTOLANO B, SPATAFORA F, BONO G, et al. Application of the wood model to lactation curves of comisana sheep[J]. Small Ruminant Research, 1997, 24(1): 7-13. DOI:10.1016/0921-4488(95)00842-X
[28]
SHAAT I. Application of the wood lactation curve in analysing the variation of daily milk yield in the Zaraibi goats in Egypt[J]. Small Ruminant Research, 2014, 117(1): 25-33.
[29]
张力, 潘林阳, 杨诗兴, 等. 哺乳单羔、双羔母湖羊的泌乳量及其泌乳曲线的研究[J]. 中国畜牧杂志, 1988(6): 23-24.
ZHANG L, PAN L Y, YANG S X, et al. Study on lactation yield and lactation curve of lactating single lamb and double lamb female Hu sheep[J]. Chinese Journal of Animal Science, 1988(6): 23-24 (in Chinese).
[30]
马友记, 董琪利, 李发弟, 等. 舍饲绵羊产后30天产奶量及乳成分变化规律[J]. 草业学报, 2013, 22(5): 287-293.
MA Y J, DONG Q L, LI F D, et al. Changes of milk yields and their nutrient composition in confined 0-30 days postpartum sheep[J]. Acta Prataculturae Sinica, 2013, 22(5): 287-293 (in Chinese).