畜牧兽医学报  2021, Vol. 52 Issue (11): 3089-3098. DOI: 10.11843/j.issn.0366-6964.2021.011.010    PDF    
荷斯坦牛皮肤褶皱厚度及体况评分性状遗传分析
张海亮1, 常瑶1, 穆柏宇1, 王凯1, 杨明路1,2, 王磊1,2, 马龙刚1,2, 宁景扬1, 郭刚3, 王雅春1     
1. 中国农业大学动物科技学院, 农业农村部动物遗传育种与繁殖(家畜)重点实验室, 畜禽育种国家工程实验室, 北京 100193;
2. 新疆农业大学动物科学学院, 乌鲁木齐 830052;
3. 北京首农畜牧发展有限公司, 北京 100029
摘要:本研究通过对荷斯坦牛颈侧部、肋部和后乳房基部皮肤褶皱厚度(分别简称为颈部、肋部和乳房皮褶厚)及体况评分(body condition score,BCS)进行大群测定,旨在探究荷斯坦牛不同部位皮肤褶皱厚度和BCS的群体特征,并建模估计荷斯坦牛不同部位皮肤褶皱厚度及BCS的遗传参数。本研究于2015—2020年夏季测定了北京地区12个规模化牧场10 915头泌乳荷斯坦牛的皮肤褶皱厚度,同时对测定牛群进行了体况评分,使用单性状和四性状动物模型分别对颈部皮褶厚、肋部皮褶厚、乳房皮褶厚和BCS进行遗传分析并估计遗传参数,并计算了上述4个性状与奶牛寿命和繁殖等重要功能性状之间的近似遗传相关。结果显示,荷斯坦牛乳房皮褶厚、颈部皮褶厚、肋部皮褶厚和体况评分分别为(7.49±1.45)mm、(7.27±1.34)mm、(11.74±1.90)mm和(2.94±0.79);乳房皮褶厚和颈部皮褶厚为中等遗传力性状,遗传力分别为0.11和0.15;肋部皮褶厚和体况评分为中高遗传力性状,遗传力分别为0.28和0.22;不同部位的皮褶厚之间存在中高遗传相关,颈部和肋部皮褶厚之间的遗传相关最高(0.69),颈部和乳房皮褶厚之间的遗传相关最低(0.21);BCS与不同部位皮褶厚之间遗传相关的差异较大,BCS与颈部、肋部皮褶厚之间为正遗传相关(0.28、0.13),BCS与乳房皮褶厚之间为负遗传相关(-0.25);体况评分和皮褶厚性状与部分重要功能性状之间存在中等遗传相关。本研究对奶牛3个部位的皮肤褶皱厚度和BCS性状进行了遗传分析,研究获得的遗传参数有助于通过皮肤褶皱厚度和BCS理解奶牛体脂沉积的遗传基础,助力我国奶牛的平衡育种。
关键词奶牛    皮肤褶皱厚度    体脂沉积    遗传参数    遗传评估    
Genetic Analysis on Skinfold Thickness and Body Condition Score Traits in Holstein Cows
ZHANG Hailiang1, CHANG Yao1, MU Baiyu1, WANG Kai1, YANG Minglu1,2, WANG Lei1,2, MA Longgang1,2, NING Jingyang1, GUO Gang3, WANG Yachun1     
1. National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China;
2. College of Animal Science, Xinjiang Agricultural University, Urumqi 830052, China;
3. Beijing Sunlon Livestock Development Co., Ltd, Beijing 100029, China
Abstract: In this study, the skinfold thickness over the neck (STN), rib (STR) and udder (STU) of Holstein cows (referred to as neck, ribs and breast skinfold thickness, respectively) and body condition score (BCS) were measured in a large population. The purpose was to explore the population characteristics of skinfold thickness in different body positions and BCS, and to estimate the genetic parameters of skinfold thickness in different body positions and BCS of Holstein cows by modeling. The STN, STR, STU and BCS of 10 915 lactating Holstein cows from 12 dairy farms were measured in summer from 2017 to 2020. The single-trait and four-trait animal models were employed to perform genetic analysis, and heritabilities and genetic correlations for STN, STR, STU and BCS were obtained. Furthermore, the approximate genetic correlations between skinfold thickness, BCS and important functional traits (longevity and reproduction traits) were calculated. The results showed that the STU, STN, STR and BCS in Holstein cows were (7.49±1.45) mm, (7.27±1.34) mm, (11.74±1.90) mm and (2.94±0.79), respectively. There were moderate heritabilities for STN (0.15) and STU (0.11), and moderate to high heritabilities for STR (0.28) and BCS (0.22). The moderate to high genetic correlations among the skinfold thickness over the different body positions were found, ranging from 0.21 (STN and STU) to 0.69 (STN and STR). There was a large variation for genetic correlations between BCS and skinfold thickness over the different body positions, ranging from -0.25 (BCS and STU) to 0.28 (BCS and STN). The moderate genetic correlations between skinfold thickness, BCS and some functional traits were found. In this study, the genetic analysis on skinfold thickness and BCS traits were performed, and the heritabilities and genetic correlations for skinfold thickness over the different body positions and BCS were estimated. The results from the current study could help us understand the genetic foundation for fat deposition in dairy cattle, and was beneficial for balance breeding in Chinese Holstein population.
Key words: dairy cow    skinfold thickness    body fat deposition    genetic parameter    genetic evaluation    

皮肤是哺乳动物身体最外层结构,是其最大的一个器官[1]。皮肤承担着保护身体、排汗、感觉冷热和压力的功能,皮肤覆盖全身,使体内各种组织和器官免受物理性、机械性、化学性和病原微生物的侵袭。牛的皮肤可分为表皮、真皮和皮下组织3层,表皮可以分为角质层、颗粒层、棘层、基底层,真皮可以分为乳头层和网状层,皮下组织可以分为脂肪组织和结缔组织[2]。对牛皮肤组织厚度的度量一般有两种方法,分别是皮肤厚度度量和皮肤褶皱厚度(skinfold thickness,简称皮褶厚)度量。取得皮肤样本后,经处理制成标本,在显微镜下使用测微尺测量表皮和真皮两层厚度可得到皮肤厚度[3];而皮褶厚则是使用卡尺在牛体直接测量得到,其包含了两层皮肤厚度与部分皮下组织的厚度。由于皮褶厚测定的无创性和简便性的特点,皮褶厚被广泛地用于牛皮肤厚度的相关研究中。

研究表明,皮褶厚可以作为反映牛体营养状况的一种标记,通过测量皮褶厚,可以了解皮下脂肪的厚度,进而可以判断个体的皮下脂肪沉积并推算全身脂肪含量。例如,Bruckmaier等[4]用肩部、垂皮、乳房和尾根4个部位皮褶厚均值减去颈部皮褶厚所得到的指数评估奶牛营养状况,结果表明,皮褶厚指数的评估效果与体况评分(body condition score, BCS)类似。据报道,牛体各部位的皮褶厚与重要经济和功能性状有关,包括蜱虫抗性、耐热性、产奶量、胴体重、活体重和繁殖性能等性状。例如,Hamid等[5]发现,产奶量与皮褶厚之间存在显著的负相关关系,日产奶量与牛颈部、垂皮、胸部、腹部和臀部的皮褶厚均值之间的回归系数为-0.57;Siddiqui等[6]发现,公牛阴囊皮褶厚与精液精子头部畸形率显著相关,阴囊皮褶厚越厚,精液中精子头部畸形率越高。由此可见,皮褶厚是牛体重要的生物特征,对牛皮褶厚的研究具有重要意义。

国际上,不同部位皮褶厚的群体规律、皮褶厚的影响因素、皮褶厚与牛重要功能或经济性状之间的关系、皮褶皱的遗传特征、与皮褶厚相关的候选基因等被持续研究。例如,Maiorano等[7]估计了内洛尔牛(Nellore)肩胛部皮褶厚的遗传力为0.12;Zhang等[8]估计了中国荷斯坦牛颈部和肋部皮褶厚的遗传力分别为0.13和0.24,颈部皮褶厚与BCS之间存在中等遗传相关0.34;Gualdrón Duarte等[9]对比利时蓝牛(Belgian Blue)皮肤厚度性状进行了全基因组关联分析和基因组选择研究,发现5号染色体上的KITLG基因与皮肤厚度性状显著关联,使用不同方法对皮肤厚度进行基因组选择的可靠性可以达到0.39~0.49。在国内,研究人员对牦牛[10-12]、延边黄牛[3]和荷斯坦牛[13-15]的皮肤厚度进行了一些研究。针对荷斯坦牛,一些研究探究了颈部和肋部皮褶厚的影响因素,报道了颈部和肋部皮褶厚的遗传参数;而对于乳房这一奶牛最重要的部位,目前仅有一篇研究初步探究了后乳房基部皮褶厚的影响因素。奶牛的乳房位于耻骨部位,由上至下依次为紧贴腹壁的基部、中间的体部和游离的乳头区[16]。奶牛后乳房基部皮肤内的皮下组织有一定的活动性,使得乳房基部的皮肤松弛易成褶可进行测量[17]

因此,本研究对北京地区荷斯坦牛颈侧部、肋部和后乳房基部的皮褶厚及BCS进行了大群测定,并建模估计了皮褶厚性状和BCS的遗传参数,旨在揭示皮褶厚性状和BCS的遗传特征,进一步探究不同部位皮褶厚性状之间的关系,为通过皮褶厚和BCS研究奶牛的体脂沉积提供参考。

1 材料与方法 1.1 试验群体

本研究中,10 915头泌乳荷斯坦牛的皮褶厚及BCS分别于2015、2016、2017、2019和2020年夏季集中测定,测定的荷斯坦牛群体分布于北京地区12个规模化奶牛场,测定母牛均为健康的泌乳母牛,各年度测定的牧场数、测定规模及测定的项目如表 1所示。此外,试验牛群的系谱记录由北京奶牛中心(北京)提供,试验牛群长寿性状和繁殖性状的育种值估计结果由奶牛育种自主创新联盟(北京)提供。

表 1 荷斯坦牛皮肤褶皱厚度及体况评分测定情况 Table 1 Measure records of skinfold thickness and body condition score in Holstein cows
1.2 测定方法

1.2.1 皮褶厚测定   测量人员使用电子数显游标卡尺(型号:Part#PD-151),在荷斯坦牛颈侧部(颈部)、最后一根肋骨肋弓(肋部)、后乳房基部(乳房基部)进行各部位皮褶厚的测量,后文分别简称为颈部皮褶厚、肋部皮褶厚和乳房皮褶厚。为降低测定误差,同年度同一场内的测定牛只由相同测定人员尽可能以同一力度揪起皮肤,测定皮肤褶皱厚度,测量结果保留2位小数。

1.2.2 体况评分   同年度同一场内测定牛只的BCS统一由2名评分人员完成,评分规则为5分制,1分代表极瘦,5分代表极胖。评分时,评定人员同时对同一头奶牛分别进行评分,两人评分差值不大于0.5分认为评分有效,否则需重新评定,以两人评分均值作为被评奶牛最终得分,用于数据分析。

1.3 数据整理

本研究共测定了10 915头奶牛的皮褶厚和BCS数据,对原始数据按以下标准进行质控:1) 剔除乳房皮褶厚异常个体(以平均值加减3倍标准差为保留范围);2) 剔除无胎次记录的个体;3) 剔除产犊日期记录缺失或错误的个体;4) 剔除无系谱记录的个体。最终获得10 138头荷斯坦牛的皮褶厚和BCS测定记录,其中乳房皮褶厚测定3 416头、颈部皮褶厚测定5 261头、肋部皮褶厚测定5 288头、体况评分测定9 094头。各性状两两组合中,同时具有两个性状表型记录的个体数在962(STN与STU)~5 253(STN与STR)头之间。

本研究将试验牛只的胎次划分为5个水平,分别为1胎、2胎、3胎、4胎、5胎及以上,将泌乳阶段划分为5个阶段,分别为1~49 (阶段Ⅰ)、50~99 (阶段Ⅱ)、100~199 (阶段Ⅲ)、200~349 (阶段Ⅳ)、350 d及以上(阶段Ⅴ)。本研究中,根据有表型个体的牛号,利用北京地区荷斯坦牛群体的大系谱追溯用于对皮褶厚和BCS性状遗传分析的小系谱,最终本研究所用的系谱文件中包含公牛2 544头、母牛29 357头。

1.4 统计分析

1.4.1 皮褶厚和体况评分性状遗传分析   利用DMU软件[18]的DMUAI模块进行方差组分估计,采用单性状和多性状动物模型对乳房皮褶厚、颈部皮褶厚、肋部皮褶厚和BCS进行遗传分析,各性状的模型分别为:

$ \begin{array}{l} Y = FM + P + S + A + E\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;模型1\\ Y = FM + P + S + Bodyside + A + E\;\;\;模型2\\ Y = FS + P + S + A + E\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;模型3 \end{array} $

采用模型1分别对颈部皮褶厚和乳房皮褶厚进行遗传分析,采用模型2对肋部皮褶厚进行遗传分析,采用模型3对BCS进行遗传分析。其中,FM表示皮褶厚测定场-测量人固定效应,不同皮褶厚性状FM效应的水平数从11到18不等;P表示奶牛胎次固定效应,共5个水平;S表示奶牛泌乳阶段固定效应,共5个水平;Bodyside表示皮褶厚测定体侧固定效应,分别为左侧、右侧和未知,共3个水平;FS表示BCS评分场-评分人固定效应,共包含40个水平;A为加性遗传随机效应,A~N(0,a2);E表示随机残差效应,E~N(0,e2)。

1.4.2 皮褶厚、体况评分性状与重要功能性状的近似遗传相关分析   基于单性状动物模型估计获得公牛的颈部、肋部、乳房皮褶厚及BCS性状估计育种值(estimated breeding value, EBV),本研究计算了上述4个性状与奶牛重要功能性状之间的近似遗传相关,这些功能性状包括生产寿命(productive life, PL)、初产日龄(age at the first calving, AFC)、初配日龄(age at the first service, AFS)、产犊后首次配种间隔(interval from calving to the first service, ICF)、青年牛首末配种间隔(interval from the first to the last inseminations in heifer, IFL_H)、经产牛首末配种间隔(interval from the first to the last inseminations in cow, IFL_C)、青年牛产犊难易(calving ease in heifer, CE_H)、经产牛产犊难易(calving ease in cow, CE_C)、青年牛死产率(stillbirth in heifer, SB_H)和经产牛死产率(stillbirth in cow, SB_C),这些性状的定义细节和评估模型见Zhang等[19]、陈紫薇等[20]和Liu等[21]的研究。

计算近似遗传相关时,分别筛选颈部、肋部、乳房皮褶厚和BCS性状估计育种值可靠性高于0.2的公牛,并匹配这些公牛上述10个功能性状的估计育种值和估计育种值可靠性,根据公式1和公式2依次分别计算皮褶厚性状、BCS性状与各个功能性状之间的近似遗传相关及其标准误,近似遗传相关和标准误的计算方法参考自Calo等[22]和Dikmen等[23]的研究。

$ {{{\mathop r\limits^ \wedge }_{g1, 2}} = \frac{{\sqrt {\left( {\sum R{L_1}} \right) \times \left( {\sum R{L_2}} \right)} }}{{\sum \left( {R{L_1} \times R{L_2}} \right)}} \times {r_{1, 2}}} $ (1)
$ SE=\sqrt{\frac{1-\overset{\wedge }{\mathop{r}}\,_{g1,2}^{2}}{n-2}} $ (2)

其中,$ {{{\mathop r\limits^ \wedge }_{g1, 2}}}$表示性状1和性状2之间的近似遗传相关; RL1RL2分别代表每头公牛性状1和性状2的估计育种值;r1, 2代表性状1和性状2的估计育种值之间的皮尔逊相关系数;SE代表近似遗传相关系数的标准误;n代表用于计算近似遗传相关的公牛头数。

2 结果 2.1 北京地区荷斯坦牛皮肤褶皱厚度及体况评分描述性统计

本研究中,北京地区荷斯坦牛乳房皮褶厚、颈部皮褶厚、肋部皮褶厚分别为(7.49±1.45)、(7.27±1.34)、(11.74±1.90) mm(表 2);不同部位皮褶厚中,肋部皮褶厚最厚,颈部最薄。泌乳荷斯坦牛的体况评分为(2.94±0.79)。不同部位皮褶厚和体况评分中,体况评分的变异系数最大,为26.87%;不同部位皮褶厚的变异系数相近,乳房皮褶厚最大,肋部皮褶厚最小,两者相差3.17%。

表 2 北京地区荷斯坦牛皮肤褶皱厚度和体况评分描述性统计 Table 2 Descriptive statistics for skinfold thickness and body condition score in Holstein cows in Beijing
2.2 各部位皮肤褶皱厚度及体况评分的遗传参数

使用单性状动物模型估计北京地区荷斯坦牛皮褶厚性状和BCS的遗传力,结果见表 3。荷斯坦牛乳房皮褶厚和颈部皮褶厚遗传力相近,为中低遗传力性状,分别为(0.11±0.03)和(0.15±0.03);肋部皮褶厚遗传力最高,为(0.28±0.04)。荷斯坦牛BCS遗传力中等,为(0.22±0.03)。

表 3 荷斯坦牛皮肤褶皱厚度和体况评分性状的遗传力、遗传相关(下三角)及表型相关(上三角)估计值 Table 3 Estimates of heritability, genetic correlation (below the diagonal) and phenotypic correlation (above the diagonal) for skinfold thickness and body condition score traits in Holstein cows

利用四性状动物模型估计各部位皮褶厚性状及BCS之间的遗传相关,各性状间的表型相关和遗传相关见表 3。不同部位的皮褶厚之间,遗传相关在0.21~0.69之间;颈部皮褶厚与肋部皮褶厚之间的遗传相关最高,为(0.69±0.08);乳房皮褶厚与颈部皮褶厚之间的遗传相关最低,为(0.21±0.14)。体况评分与不同部位的皮褶厚之间的遗传相关相差较大,与颈部皮褶厚之间存在正的遗传相关,相关系数为(0.28±0.10);与乳房皮褶厚之间存在中等负遗传相关,相关系数为(-0.25±0.13)。

2.3 胎次和泌乳阶段对各部位皮肤褶皱厚度及体况评分的影响

荷斯坦牛各部位皮褶厚和体况评分不同胎次的最佳线性无偏估计值变化趋势如图 1所示。随胎次的增加,乳房皮褶厚、颈部皮褶厚和肋部皮褶厚均总体上呈现变薄的趋势,头胎奶牛皮褶厚最厚,4胎奶牛乳房皮褶厚和肋部皮褶厚最薄,5胎及以上奶牛颈部皮褶厚最薄,各部位皮褶厚分别下降0.19、0.57和0.70 mm。此外,BCS随奶牛胎次增加的变化较小,前4胎奶牛的体况评分相近,5胎奶牛的体况评分最低,变化幅度为0.13分。

图 1 不同胎次荷斯坦牛皮褶厚和体况评分性状最佳线性无偏估计值变化趋势 Fig. 1 The best linear unbiased estimates of skinfold thickness and body condition score traits in Holstein cows with different parities

荷斯坦牛各部位皮褶厚和体况评分不同泌乳阶段的最佳线性无偏估计值变化趋势如图 2所示。各部位皮褶厚随泌乳阶段的变化趋势较为一致,产后泌乳早期皮褶厚(1~49 d)最厚,泌乳中期(100~ 199 d)最薄,泌乳后期略有上升。同泌乳期不同泌乳阶段,乳房皮褶厚、颈部皮褶厚和肋部皮褶厚的变化幅度分别为0.52、0.34和0.38 mm。此外,体况评分随泌乳阶段的延长也呈现先下降后上升的趋势,泌乳高峰期(50~99 d)最低,后随着泌乳天数的增加不断增加,同泌乳期不同泌乳阶段体况评分的变化幅度为0.80分。

图 2 不同泌乳阶段荷斯坦牛皮肤褶皱厚度和体况评分最佳线性无偏估计值变化趋势 Fig. 2 The best linear unbiased estimates of skinfold thickness and body condition score traits in Holstein cows at different lactation stages
2.4 各部位皮褶厚和体况评分与重要功能性状之间的近似遗传相关

本研究中,荷斯坦牛各部位皮褶厚性状、BCS与寿命和繁殖性状之间的近似遗传相关见表 4。由表 4可知,3个部位的皮褶厚性状与生产寿命性状之间存在中等遗传相关,遗传相关系数介于0.13~0.26;皮褶厚较厚的牛,其长寿性表现越好。BCS和颈部皮褶厚性状与经产牛首末配种间隔性状之间存在中等遗传相关,遗传相关系数分别为-0.36和-0.33;BCS越高且颈部皮褶厚越厚的经产牛,其首末配种间隔越短。肋部皮褶厚与青年牛首末配种间隔、乳房皮褶厚与青年牛首末配种间隔之间分别存在中等遗传相关,遗传相关系数分别为0.34和0.19;皮褶厚越厚,青年牛首末配种间隔越长。此外,颈部皮褶厚和乳房皮褶厚分别与经产牛死产率性状存在中等遗传相关,遗传相关系数分别为0.17和0.22;皮褶厚越厚,经产牛死产率越高。

表 4 荷斯坦牛皮肤褶皱厚度、体况评分与重要功能性状之间的近似遗传相关 Table 4 The approximately genetic correlations between skinfold thickness, body condition score traits and important functional traits in Holstein cows
3 讨论

北京地区荷斯坦牛乳房皮褶厚的群体均值为7.15 mm,颈部皮褶厚的群体均值为7.12 mm,肋部皮褶厚的群体均值为11.71 mm。本研究中,荷斯坦牛不同部位皮褶厚相差较大,乳房与颈部皮褶厚远小于肋部皮褶厚,前人研究也发现同一品种奶牛不同部位的皮褶厚存在较大差异[24]。荷斯坦牛颈部和肋部皮褶厚均分别处于各个品种的报道范围内,如颈部皮褶厚为3.33~7.86 mm,肋部皮褶厚为11.40~12.16 mm[8]。Rahman和Gill[25]在534头摩拉水牛中发现,乳房皮肤厚度平均值为11.4 mm,大于本研究中荷斯坦牛的乳房皮褶厚。摩拉水牛的生活环境为亚热带气候[26],温带和热带品种牛的皮褶厚存在着显著的品种差异[27],水牛皮褶厚均值是黄牛皮褶厚的2.9倍[28]。本研究中,北京地区泌乳荷斯坦牛的BCS均值为2.94,在前人报道的我国各地区泌乳荷斯坦牛的BCS范围内(2.55~3.09)[29-30]。本研究对荷斯坦牛3个部位的皮褶厚及BCS进行了大规模测定,测定结果与前人研究相符;中国荷斯坦牛群体中,上述4个性状均存在较大的变异,变异系数在16.18%~26.87%之间。

本研究中,荷斯坦牛颈部、肋部和乳房皮褶厚的遗传力在0.11~0.28之间,颈部皮褶厚最低,肋部皮褶厚最高。Zhang等[8]利用4 307头荷斯坦牛的相关数据估计得到荷斯坦牛颈部、肋部皮褶厚的遗传力分别为0.13±0.03、0.26±0.04,分别与本研究结果相近;Maiorano等[7]对17 940头内洛尔牛的肩胛部皮褶厚进行遗传参数估计,发现肩胛部皮褶厚的遗传力为0.12。由此可见,荷斯坦牛皮褶厚在牛体不同部位的遗传力存在较大差别,荷斯坦牛皮褶厚是一个中高遗传力性状。本研究中,荷斯坦牛BCS的遗传力为0.22,处于报道的荷斯坦牛BCS估计遗传力范围内(0.11~0.30)[31-32],并与我国南方荷斯坦牛中BCS的遗传力估计值相近(0.19)[33]。本研究中,4个性状估计遗传力的标准误较小,说明遗传分析的结果较为可靠;在不同群体之间,同一性状遗传参数的差异可能受牛群遗传背景、遗传参数估计方法、数据量等因素的影响。

本研究发现,不同部位的皮褶厚之间存在中等到较高的遗传相关,相关系数范围为0.21~0.69。BCS与不同部位皮褶厚之间遗传相关的变异较大,与颈部皮褶厚最高(0.28),与肋部皮褶厚最低(0.13),与乳房皮褶厚之间存在负遗传相关(-0.25);不同部位的皮褶厚与BCS的遗传基础不同,颈部和肋部皮褶厚可与BCS一起反映奶牛的体脂沉积状态。本研究中,奶牛各部位皮褶厚和BCS随胎次的变化趋势相近,均随着胎次的增加逐渐下降;而随着泌乳阶段的增加也均呈现先减小后增加的趋势,BCS在泌乳高峰期(50~99 d)最低,各部位皮褶厚在泌乳中期(100~199 d)最低。由此可见,皮褶厚和BCS性状对奶牛生理阶段的变化较为敏感,能够反映不同时期奶牛体脂沉积的变化。

4 结论

本研究对北京地区荷斯坦牛颈侧部、肋部和后乳房基部皮肤褶皱厚度(简称为皮褶厚)以及体况评分进行了大群测定,并对上述性状进行了遗传分析。后乳房基部皮褶厚(0.11)和颈侧部皮褶厚(0.15)为中等遗传力性状,肋部皮褶厚(0.28)和体况评分(0.22)为中高遗传力性状;体况评分、皮褶厚与部分重要功能性状之间存在中等遗传相关,包括生产寿命、青年牛首末配种间隔、经产牛首末配种间隔和青年牛死产率性状,遗传相关系数绝对值介于0.14~0.36之间;各部位皮褶厚之间存在低到中等的遗传相关,且与体况评分之间的遗传相关差异较大,不同部位皮褶厚性状的遗传基础存在差异。本研究有助于通过皮肤褶皱厚度和体况评分理解奶牛体脂沉积的遗传特征,助力我国奶牛的平衡育种。

致谢: 感谢中国农业大学动物科技学院牛百科青年计划成员于2015、2016、2017和2019年夏季在各个牧场辛苦测定皮褶厚和体况评分。

参考文献
[1]
曹迷霞. 成年牦牛皮肤中附属物、毛细血管和神经的组织学研究[D]. 兰州: 甘肃农业大学, 2018.
CAO M X. Histological study of appendages, capillary and nerve in the skin of adult yak[D]. Lanzhou: Gansu Agricultural University, 2018. (in Chinese)
[2]
谢念难. 甘肃牦牛皮肤的微细结构[J]. 中国牦牛, 1987(2): 1-7.
XIE N N. The microstructure of skin in yak in Gansu[J]. China Yak, 1987(2): 1-7. (in Chinese)
[3]
王富金, 李寿万. 延边黄牛皮肤厚度的观察测量[J]. 延边农学院学报, 1990(1): 43-46.
WANG F J, LI S W. Measurement on skin thickness in Yanbian yellow cattle[J]. Journal of Yanbian Agricultural College, 1990(1): 43-46. (in Chinese)
[4]
BRUCKMAIER R M, GREGORETTI L, JANS F, et al. Longissimus dorsi muscle diameter, backfat thickness, body condition scores and skinfold values related to metabolic and endocrine traits in lactating dairy cows fed crystalline fat or free fatty acids[J]. J Vet Med Ser A, 1998, 45(1-10): 397-410. DOI:10.1111/j.1439-0442.1998.tb00842.x
[5]
HAMID M A, HUSAIN S M I, KHAN M K I, et al. Skin thickness in relation to milk production of crossbred cows[J]. Pak J Biol Sci, 2000, 3(9): 1525-1529. DOI:10.3923/pjbs.2000.1525.1529
[6]
SIDDIQUI M A R, BHATTACHARJEE J, DAS Z C, et al. Crossbred bull selection for bigger scrotum and shorter age at puberty with potentials for better quality semen[J]. Reprod Domest Anim, 2008, 43(1): 74-79.
[7]
MAIORANO A M, OLIVEIRA M C S, CYRILLO J N S G, et al. Genetic study of skin thickness and its association with postweaning growth in Nellore cattle: estimation of the genetic parameters[J]. Genet Mol Res, 2016, 15(1): gmr.15017124.
[8]
ZHANG H, LIU A, LI X, et al. Genetic analysis of skinfold thickness and its association with body condition score and milk production traits in Chinese Holstein population[J]. J Dairy Sci, 2019, 102(3): 2347-2352. DOI:10.3168/jds.2018-15180
[9]
GUALDRÓN DUARTE J L, GORI A S, HUBIN X, et al. Performances of adaptive MultiBLUP, Bayesian regressions, and weighted-GBLUP approaches for genomic predictions in Belgian Blue beef cattle[J]. BMC Genomics, 2020, 21(1): 545. DOI:10.1186/s12864-020-06921-3
[10]
俞世福, 陈爱霞, 杨枝. 青海牦牛皮肤组织学结构与生态的研究[J]. 青海畜牧兽医学院学报, 1985(2): 1-6, 53.
YU S F, CHEN A X, YANG Z. Study on the histological structure and ecology of skin in Qinghai yak[J]. Journal of Qinghai Animal Husbandry and Veterinary Medicine College, 1985(2): 1-6, 53. (in Chinese)
[11]
欧阳熙, 王茜飞. 牦牛耐寒力的研究——牦牛皮肤组织学结构的初步观察[J]. 西南民族学院学报: 畜牧兽医版, 1984(3): 13-21.
OUYANG X, WANG Q F. A study on cold tolerance in Yak——Preliminary observation of the histological structure of skin in yak[J]. Journal of Southwest Minzu University: Natural Science Edition, 1984(3): 13-21. (in Chinese)
[12]
于川, 崔燕, 杨雪, 等. 不同年龄牦牛皮肤厚度变化研究[J]. 甘肃农业大学学报, 2017, 52(1): 21-25.
YU C, CUI Y, YANG X, et al. Thickness changes of skin in different age yak (Bos grunniens)[J]. Journal of Gansu Agricultural University, 2017, 52(1): 21-25. (in Chinese)
[13]
张海亮, 徐伟, 李想, 等. 北京地区荷斯坦牛皮肤褶皱厚度及其影响因素分析[J]. 中国畜牧杂志, 2016, 52(20): 57-62.
ZHANG H L, XU W, LI X, et al. Measurement and influencing factors analysis of skinfold thickness in Holstein in Beijing[J]. Chinese Journal of Animal Science, 2016, 52(20): 57-62. (in Chinese)
[14]
娄文琦, 李源韬, 常瑶, 等. 北京地区荷斯坦奶牛皮肤褶皱厚度与泌乳性能关系的研究[J]. 中国畜牧兽医, 2018, 45(2): 375-383.
LOU W Q, LI Y T, CHANG Y, et al. Study on the relationship between the skinfold thickness and lactation performance of Holstein in Beijing[J]. China Animal Husbandry & Veterinary Medicine, 2018, 45(2): 375-383. (in Chinese)
[15]
穆柏宇, 赵琦琪, 常瑶, 等. 北京地区荷斯坦牛乳房皮肤厚度影响因素分析及其与泌乳性能关系的研究[J]. 中国畜牧杂志, 2019, 55(3): 55-59.
MU B Y, ZHAO Q Q, CHANG Y, et al. Study on influencing factors analysis of udder skin thickness and relationship with lactation performance in Holstein in Beijing[J]. Chinese Journal of Animal Science, 2019, 55(3): 55-59. (in Chinese)
[16]
赵珂. 奶牛乳腺上皮细胞葡萄糖摄取的调控及其对乳成分合成的影响研究[D]. 杭州: 浙江大学, 2011.
ZHAO K. Studies on the regulation of glucose transport and utilization in bovine mammary epithelial cells[D]. Hangzhou: Zhejiang University, 2011. (in Chinese)
[17]
杨银凤. 家畜解剖学及组织胚胎学[M]. 4版. 北京: 中国农业出版社, 2011.
YANG Y F. Livestock anatomy, histology and embryology[M]. 4th ed. Beijing: China Agricultural Press, 2011. (in Chinese)
[18]
MADSEN P, SØRENSEN P, SU G, et al. DMU—A package for analyzing multivariate mixed models[Z]. Minas Gerais, Brazil: Instituto Prociência, 2006.
[19]
ZHANG H L, LIU A X, WANG Y C, et al. Genetic parameters and genome-wide association studies of eight longevity traits representing either full or partial lifespan in Chinese Holsteins[J]. Front Genet, 2021, 12: 634986. DOI:10.3389/fgene.2021.634986
[20]
陈紫薇, 师睿, 罗汉鹏, 等. 宁夏地区荷斯坦牛青年牛繁殖性状遗传参数估计[J]. 畜牧兽医学报, 2021, 52(2): 344-351.
CHEN Z W, SHI R, LUO H P, et al. Estimation of genetic parameters of reproductive traits of Holstein heifers in Ningxia[J]. Acta Veterinaria et Zootechnica Sinica, 2021, 52(2): 344-351. (in Chinese)
[21]
LIU A X, LUND M S, WANG Y C, et al. Variance components and correlations of female fertility traits in Chinese Holstein population[J]. J Anim Sci Biotechnol, 2017, 8: 56. DOI:10.1186/s40104-017-0189-x
[22]
CALO L L, MCDOWELL R E, VANVLECK L D, et al. Genetic aspects of beef production among Holstein-Friesians pedigree selected for milk production[J]. J Anim Sci, 1973, 37(3): 676-682. DOI:10.2527/jas1973.373676x
[23]
DIKMEN S, COLE J B, NULL D J, et al. Heritability of rectal temperature and genetic correlations with production and reproduction traits in dairy cattle[J]. J Dairy Sci, 2012, 95(6): 3401-3405. DOI:10.3168/jds.2011-4306
[24]
PATEL G U, ANDERSON W D. Variation of skin thickness in dairy cattle[J]. Emp J Exp Agric, 1958, 26: 18-24.
[25]
RAHMAN S M, GILL R S. Relation between milk vein and udder skin thickness with milk yield in Murrah buffaloes[J]. Indian J Anim Prod Manage, 1991, 7(3): 119-123.
[26]
赵家才, 何华川, 王春兰, 等. 摩拉水牛改良云南本地水牛的效果研究[J]. 中国牛业科学, 2007, 33(1): 4-7.
ZHAO J C, HE H C, WANG C L, et al. Crossbreeding effects of Murry Buffalo on native Buffaloes in Yunnan[J]. China Cattle Science, 2007, 33(1): 4-7. DOI:10.3969/j.issn.1001-9111.2007.01.002 (in Chinese)
[27]
WALKER C A. The skin thickness of cattle in Northern Rhodesia[J]. J Agric Sci, 1957, 49(2): 211-213. DOI:10.1017/S0021859600036194
[28]
李克平, 陈俊菲, 王石平, 等. 水牛皮肤结构的研究[J]. 华中农学院学报, 1983, 2(2): 35-43.
LI K P, CHEN J F, WANG S P, et al. A skin structural study in water buffalo[J]. Journal of Huazhong Agricultural College, 1983, 2(2): 35-43. (in Chinese)
[29]
常瑶, 朱莹琳, 施佳庆, 等. 西荷杂交牛与荷斯坦牛重要经济性状对比[J]. 畜牧兽医学报, 2018, 49(2): 270-281.
CHANG Y, ZHU Y L, SHI J Q, et al. A comparison study on economically important traits between Holstein and Simmental×Holstein Crossbred cows[J]. Acta Veterinaria et Zootechnica Sinica, 2018, 49(2): 270-281. (in Chinese)
[30]
巨星, 黄锡霞, 葛建军, 等. 中国荷斯坦奶牛体况评分及其影响因素分析[J]. 中国畜牧杂志, 2019, 55(4): 49-52.
JU X, HUANG X X, GE J J, et al. A study on body condition score and its influencing factors in Chinese Holstein[J]. Chinese Journal of Animal Science, 2019, 55(4): 49-52. (in Chinese)
[31]
ZINK V, ŠTÍPKOVÁ M, LASSEN J. Genetic parameters for female fertility, locomotion, body condition score, and linear type traits in Czech Holstein cattle[J]. J Dairy Sci, 2011, 94(10): 5176-5182. DOI:10.3168/jds.2010-3644
[32]
BATTAGIN M, SARTORI C, BIFFANI S, et al. Genetic parameters for body condition score, locomotion, angularity, and production traits in Italian Holstein cattle[J]. J Dairy Sci, 2013, 96(8): 5344-5351. DOI:10.3168/jds.2012-6352
[33]
LIU S B, YI J M, YANG L. Genetic parameters estimates for locomotion score, body condition score and final type score of Holstein cattle in southern China[J]. J Appl Anim Res, 2013, 41(2): 240-243. DOI:10.1080/09712119.2012.739090

(编辑   郭云雁)