中国荷斯坦牛CXCR1基因编码区SNP多态与临床乳房炎和生产寿命的关联分析

王梦琦，倪炜，张慧敏，杨章平，王西朴，蒋彦森，毛永江,

（1扬州大学动物科学与技术学院 江苏扬州 225009；2江苏君乐宝乳业有限公司，江苏徐州 221721）

中国荷斯坦牛CXCR1基因编码区SNP多态与临床乳房炎和生产寿命的关联分析

王梦琦1，倪炜1，张慧敏1，杨章平1，王西朴2，蒋彦森2，毛永江1,2

（1扬州大学动物科学与技术学院 江苏扬州 225009；2江苏君乐宝乳业有限公司，江苏徐州 221721）

【目的】探讨中国荷斯坦牛CXCR1基因编码区SNP突变与临床乳房炎和生产寿命的相关性。【方法】根据CXCR1基因编码区序列，利用PCR-直接测序法对低SCS和高SCS样本各20个样本进行SNP筛查，最后对所选择的4个SNP位点利用飞行质谱法对866头中国荷斯坦牛进行检测，同时收集所检测牛只临床乳房炎和生产寿命等信息，利用多因素方差分析法、Logistic回归、Cox生存回归等方法分析以上SNP位点突变与临床乳房炎发生次数和生产寿命的相关性。【结果】CXCR1基因编码区共发现13个SNP位点（291 C＞T、333 C＞G、337 A＞G、365 C＞T、570 A＞G、642 A＞G、735 C＞G、816 A＞C、819 A＞G、980 A＞G、995 A＞G、1008 C＞T和1068 A＞G），分为4个连锁群，随后从每个连锁群中各选择1个SNP位点（642 A＞G, 816 A＞C, 980 A＞G和1068 A＞G），利用飞行质谱法对大样本中国荷斯坦牛进行检测。4个SNP位点共有9种单倍型，其中单倍型GAGG频率最高（0.3141），而单倍型ACAA频率最低（0.0017）。CXCR1-642与2胎牛患临床乳房炎次数有显著相关（P＜0.05），AG基因型个体2胎奶牛患临床乳房炎次数显著高于AA基因型（P＜0.05），CXCR1-816 AA 基因型个体3胎奶牛患临床乳房炎次数显著低于AC和CC基因型（P＜0.05），其它SNP位点与各胎次临床乳房炎发生次数均无显著相关（P＞0.05）。CXCR1-816与奶牛离群月龄有极显著相关（P＜0.01），与奶牛生产月龄和离群胎次有显著相关（P＜0.05），CXCR1-816 AA基因型个体生产月龄、离群月龄和离群胎次均显著高于CC基因型个体（P＜0.05），而其它SNP位点对生产月龄、离群月龄和离群胎次均无显著相关（P＞0.05）。Cox生存分析表明：只有CXCR1-816位点与奶牛生存时间有显著相关(P＜0.05)，CXCR1-816 CC基因型个体在各时间段的生存概率均低于AA和AC基因型个体。【结论】CXCR1-816 A＞C突变与中国荷斯坦牛患临床乳房炎次数和生产寿命有显著相关，在进一步验证其功能后，可用于中国荷斯坦牛生产寿命的分子标记辅助选择。

中国荷斯坦牛；CXCR1；SNP；临床乳房炎；生产寿命

0 引言

【研究意义】乳房炎是奶牛最常见的疾病之一，也是造成奶牛业损失最严重的疾病[1-2]。奶牛患乳房炎后导致奶品质降低、产奶量下降，甚至危害到消费者的健康。奶牛生产寿命是指奶牛从第1次产犊到死亡或者淘汰之间的时间，它能够反映奶牛避免被淘汰的能力。奶牛的生产寿命是一个重要的经济指标，它在增加牧场主的收入和育种方面有重要的作用。有研究表明：荷斯坦牛的生产寿命与终生产奶量呈显著正相关[3]，奶牛的产奶胎次从第3胎提高到第4胎，其效益将增加11%—13%[4]。奶牛生产寿命与生产性能、繁殖性能及健康状况等紧密相关，其中临床乳房炎是造成奶牛淘汰的主要原因之一。临床乳房炎的遗传力为0.03—0.2[5]，奶牛生产寿命的遗传力为0.029—0.072[6]。由于临床乳房炎和生产寿命的遗传力均较低，直接选择很难获得理想的遗传进展。因此，结合分子标记辅助选择技术对乳房炎和生产寿命的研究越来越受到关注。【前人研究进展】 趋化因子受体1（chemokine receptor 1, CXCR1）属于G蛋白偶联受体的家族成员，目前至少有两种不同的类型（CXCR1和CXCR2）[7]。I型受体（CXCR1）特异性地与白细胞介素-8（Interleukin-8，IL8）结合，主要分布于中性粒细胞、单核细胞和G细胞；II型受体（CXCR2）除了与IL8结合外，还可以和生长调节癌基因（Growth-regulated oncogene, GRO）、巨噬细胞炎性蛋白22（Macrophage inflammatory protein-2，MIP-2）和中性粒细胞激活蛋白-2（Neutrophil activating protein-2，NAP-2）等结合，CXCR2主要分布于中性粒细胞和髓样前体细胞系。CXCR1和CXCR2均具有介导活化中性粒细胞的作用，可促进中性粒细胞脱颗粒、释放贮存酶，增强中性粒细胞的吞噬功能，启动超氧离子释放,导致机体局部炎症反应，达到杀灭病原菌的目的[8-9]。大量研究表明：IL8和CXCR1基因及其突变与奶牛乳房均存在一定的关系[10-16]。周雷等[17]通过对中国荷斯坦奶牛CXCR1基因外显子2进行多态性分析发现4个突变位点291 C＞T、333 C＞T、337 A＞G 和365 C＞T，并且发现CCCTGGCC基因型组合个体具有低SCS和高产奶量的优势。奶牛的生产寿命与牛场的生产性能、饲养管理、日粮配合、环境、季节、其它疾病等非遗传因素有一定关系[18]，同时与遗传也有一定关系[19-20]。近年来，国内外围绕影响奶牛生产寿命的分子标记进行了部分研究，发现了部分有价值的结果。ASHWELL等[21-22]分别在奶牛第2、12、16、21、23号染色体上发现了与其生产寿命相关的微卫星；HEYEN等[23]在BTA21上发现了长度为85 cM的微卫星ILSTS054；KHATIB[19-20]、KOMISAREK等[24]、JOHN等[25]又发现了ATP结合盒子亚家族G成员（ATP-binding cassette sub-family G member 2，ABCG2）、过氧化物酶体增殖物增殖受体γ辅激活因子（Peroxisomeproliferators-activeted receptor-γcoactovator 1 A，PPARGC1A）、氧化低密度脂蛋白受体1（Oxidized low-density lipoprotein receptor 1，OLR1）和硬脂酰辅酶A去饱和酶1（stearoyl-CoA desaturase 1, SCD1）与生产寿命有一定程度的相关性。CHEBEL 等[26]发现Leptin 基因R4C位点与奶牛生产寿命相关疾病的发生显著相关，TT基因型个体真胃变位的发生率较高，而CT型各种疾病的发病率最低。【本研究切入点】国内周雷等[17]、徐敏等[27]和官永强等[28]检测了中国荷斯坦牛CXCR1基因部分编码区（coding sequence，CDS）序列SNP多态性，并与泌乳性状进行了关联分析，但尚未见对CXCR1基因整个CDS区SNP检测的报告，也未见CXCR1基因CDS区SNP多态与临床乳房炎和生产寿命相关分析的报道。【拟解决的关键问题】本研究的目的是在检测中国荷斯坦牛CXCR1基因CDS区所有SNP突变基础上，同时分析SNP多态与临床乳房炎和生存寿命的相关性，以期为中国荷斯坦牛临床乳房炎及生产寿命的分子标记辅助选择提供参考。

1 材料与方法

1.1 试验材料

2011年7—8月于江苏省某大型奶牛场随机采集中国荷斯坦牛血样共866头。采用尾静脉采血（10 mL/头），ACD抗凝，-20℃冷冻保存备用。牧场基本情况如下：全场共有奶牛10 000余头，其中成年泌乳奶牛近6 000头，散栏饲养，鱼骨式挤奶大厅进行挤奶，日挤奶三次，采用全混合日粮（Total Mixed Ration, TMR）饲喂。

1.2 DNA提取及CXCR1基因CDS区SNP突变检测

奶牛血液基因组DNA采用常规酚氯仿提取法，TE溶解，取部分DNA样品稀释至100 ng·µL-1，-20℃保存备用。

根据GenBank公布的牛CXCR1基因CDS区序列（登录号：NM_001105038.1），用Primer 5.0软件设计引物1对如下：F:ATGACAATCATCCTGAA AGA，R：TCAGAGGGTAGTAGACGTGT。扩增产物长度为1 082 bp，PCR扩增体系如下：总体积20 μL，其中包括10×buffer 2.0 μL，25 mmol·L-1Mg2+1.5 μL、dNTP（10 mmol·L-1）0.5 μL、Taq DNA聚合酶（5 U·µL-1）0.3 μL，上下游引物各（10 pmol·µL-1）1.0 μL,模板DNA（100 ng·µL-1）1.0 μL、ddH2O 12.7 μL。扩增程序为：94℃预变性5 min；94℃变性40 s，56℃复性40 s，72℃延伸10 min，30个循环，最后72℃延伸10 min，4℃保存。PCR扩增产物经1%琼脂糖凝胶（0.05% EB）电泳0.5 h后，紫外灯下观察结果，并拍照。

从所采样本中选择20个高SCS样本（5.85±1.22）和20个低SCS样本（1.02±0.05），进行常规PCR扩增，经琼脂糖凝胶电泳检测其扩增效果，后送至上海生物公司正反测序。用DNAMAN软件对测序所得结果进行比对，寻找突变位点，并初步计算各基因各位点连锁程度。经对40个样本的初步分析，CXCR1基因CDS区共发现13个SNP位点，分为4个连锁群，故从4个连锁群中各选择1个SNP位点进行后续分析，分别是：642 A＞G、816 A＞C、980 A＞G和1068 A＞G 4个SNP位点。对大样本CXCR1 SNP突变检测采用飞行时间质谱法（MassARRAY® MALDI-TOF System，Sequenom，Inc.，USA）。同时，为保证SNP分析结果可靠性，其中40个样本重复测定2次（测试员并不知道这个40个是重复样本）。结果表明该方法SNP分型结果准确性为100%。

1.3 临床乳房炎发生次数和生产寿命相关数据的收集

奶牛临床乳房炎以可以看到明显的乳房红肿、触摸有疼痛感和硬块或乳样呈豆腐样等症状为判断条件，以兽医治疗记录为准。如果某奶牛在治愈后10d内又发病，则算一次病例；治愈后10d以上又发病的，则算为新发病例。根据该原则记录每头奶牛从2010年1月到2014年12月共5年内各胎次发生临床乳房炎的次数。

根据采样牛只的出生时间、第一胎产犊时间、淘汰时间和胎次等信息，计算出初产月龄、生产月龄（或在群月龄）、离群月龄、离群胎次等信息，作为生产寿命相关指标，进行统计分析。

1.4 统计分析

1.4.1 遗传学分析 用遗传学软件SHEsis进行常规群体遗传学统计分析（包括基因频率、基因型频率、Hardy-Weinberg平衡检测等），同时进行连锁不平衡分析[29]。

1.4.2 关联分析 用多因素方差分析模型分析CXCR1基因CDS区各SNP位点不同基因型与奶牛各胎次临床乳房炎次数和生产寿命相关指标的相关性，模型如下：

式中：Y为各胎次临床乳房炎次数和生产寿命相关指标的观察值，μ为群体均值，S 为初产季节的效应值，M 为初产月龄的效应值，G的基因型或单倍型的效应值，e为随机误差。根据当地气候特点，初产季节划分如下：3—5月为春季，6—8月为夏季，9—11月为秋季，12月至次年2月为冬季。

由于部分牛只记录不完整，只选择同时具有1—3胎临床乳房炎记录的牛只进入方差分析。经筛选，最终有721头奶牛用于分析。对于生产寿命相关指标，只选择已淘汰，且同时具有出生、产犊等相关信息的牛只进行分析。经筛选，最终有647头奶牛用于分析。

1.4.3 生产寿命的生存分析 用Cox回归模型计算CXCR1基因CDS区各SNP位点不同基因型对生产寿命相关指标的影响，并制作生存曲线，由统计软件SPSS（Ver 16.0）完成，模型如下：

式中：

h0(t)：基准风险函数，即所有变量取零时的t时刻的风险函数

X1、X2、X2和 X4：各SNP位点变量

β1、β2、β3和β4：各SNP位点变量的回归系数

2 结果

2.1 CXCR1基因CDS区SNP筛查

通过对20头高SCS和低SCS奶牛用所设计的引物PCR扩增并测序，发现CXCR1基因CDS区共发现13个SNP位点，详细见表1。经连锁分析，13个SNP位点分为4个连锁群（图1），故从4个连锁群中各选择1个SNP位点进行后续分析，分别是：642 A＞G、816 A＞C、980 A＞G和1068 A＞G 4个SNP位点。

另外，由于CXCR1基因SNP的变化均发生在CDS区，因此SNP的变化可能导致所编码氨基酸的变化。根据CXCR1基因CDS区SNP导致氨基酸变化情况，制作各位点氨基酸变化情况如图2 所示。CXCR1基因CDS区共有5个非同义突变，I113V、K327R和H332R 3个非同义突变是完全连锁的，属于如前所述第2个连锁群，V122A和H245Q属于第1连锁群。第3、4两个连锁群均是同义突变，不导致氨基酸的变化。2个非同义突变连锁群可以组成4种氨基酸单倍型。

2.2 基因（型）频率分布及连锁不平衡分析

利用飞行质谱法对866头奶牛群体CXCR1基因CDS区4个SNP位点进行分型，所得基因频率、基因型频率分布及Hardy-Weinberg平衡检验结果见表2。CXCR1-642、CXCR1-816、CXCR1-980和CXCR1-1068 4个SNP位点优势基因型分别为GG、CC、GG和GG型，频率分别为0.580、0.466、0.607和0.658，优势基因分别为G、C、G和G，频率分别为0.758、0.688、0.784和0.811。经分析，CXCR1基因4个SNP位点均处于Hardy-Weinberg平衡状态。另外，CXCR1基因4个SNP位点共有9种单倍型（表3），其中单倍型GAGG频率最高（0.3141），而单倍型ACAA频率最低（0.0017）。

表1 CXCR1基因CDS区SNPsTable 1 The SNPs in the CDS of CXCR1 gene

2.3 CXCR1基因不同SNP位点临床乳房炎发病次数

CXCR1基因不同SNP位点各胎次临床乳房炎发生次数见表4。由表4可知：随着胎次增加，奶牛患临床乳房炎发生次数增加较快：1、2、3胎及1—3胎共计发生临床乳房炎的平均次数分别为0.01±0.01、0.52±0.03、0.91±0.05和1.44±0.07。CXCR1-642与2胎牛、CXCR1-816与3胎牛临床乳房炎发生次数有显著相关（P＜0.05），其它SNP位点与各胎次临床乳房炎发生次数均无显著相关（P＞0.05）。多重比较表明：2胎牛CXCR1-642 AG基因型个体临床乳房炎发生次数（0.61±0.06）显著高于AA型（0.33±0.09），

图1 40个个体CXCR1基因CDS区SNP位点连锁图Fig. 1 The linkage map of SNP in the CDS of CXCR1 gene for 40 Holstein cows

表2 CXCR1基因 CDS区4个SNP位点等位基因、基因型频率及H-W平衡检验Table 2 The genotypic frequencies, alleles, and values of X2test significance for 4 SNP in the CDS of CXCR1 gene

3胎牛CXCR1-816 AA基因型个体临床乳房炎发生次数（0.55±0.10）显著低于AC和CC型个体（0.99± 0.08和0.90±0.07）。

2.4 CXCR1基因不同SNP位点的生产寿命

图2 CXCR1基因CDS区AA变化图Fig. 2 The change of AA in the CDS of CXCR1 gene

表3 CXCR1基因4个SNP位点单倍型频率Table 3 The haplotypes frequencies of 4 SNP, in the CDS of CXCR1 gene

表4 CXCR1基因不同SNP位点各胎次临床乳房炎发生次数Table 4 The number of CM for SNPs of CXCR1 of different parities

在本研究中，结合奶牛出生、初产季节、初产月龄、淘汰日期等信息，用生产月龄、离群月龄和离群胎次3个指标来度量奶牛的生产寿命。生产月龄从13到58个月不等（图3），其中36个月最多，其次为37个月，平均为36.86±0.33月。离群月龄从26—85不等（图4），其中64最多，其次为65，平均为64.19±8.30。离群胎次从1—5胎不等（图5），其中3胎最多，其次为4胎，平均为3.04±0.78。

图3 奶牛生产月龄分布Fig. 3 The distribution of production months of cows

图4 离群月龄分布Fig. 4 The distribution of culling months of cows

图5 离群胎次分布Fig. 5 The distribution of culling parities of cows

CXCR1基因CDS区不同SNP位点奶牛生产月龄、离群月龄和离群胎次见表5。由表5可知：CXCR1-816基因型与奶牛离群月龄有极显著相关（P＜0.01），CXCR1-816基因型与奶牛生产月龄和离群胎次有显著相关（P＜0.05），而其它SNP位点与生产月龄、离群月龄和离群胎次均无显著相关（P＞0.05）。多重比较表明：CXCR1-816 AA基因型个体生产月龄、离群月龄和离群胎次均显著高于CC基因型个体（P＜0.05）。

2.5 生存分析

CXCR1基因CDS区4个SNP位点的奶牛生存分析结果表明：仅有CXCR1-816位点与奶牛生存时间有显著相关（P＜0.05，表6）。进一步分析发现：CXCR1-816 CC基因型个体生存时间显著低于AA和AC基因型（表7）。CXCR1-816 不同基因型生存曲线见图6。由图6可知：CXCR1-816 CC基因型个体在各时间段的生存概率均低于AA和AC基因型个体。

表6 CXCR1基因CDS区4个SNP位点奶牛的Cox生存分析Table 6 The Cox survival analysis of 4 SNPs of CXCR1 gene for dairy cows

3 讨论

3.1 CXCR1基因遗传多态性分析

图6 奶牛CXCR1-816不同基因型生存曲线Fig. 6 The survival curve of cows for different genotypes of CXCR1-816

表7 CXCR1-816不同基因型生存时间平均数与中数的估计Table 7 The estimation value of mean and median of survival time for different genotypes of CXCR1-816

趋化因子受体和趋化因子在机体炎症发生时起着非常重要的作用，当奶牛发生乳房炎时CXCR1基因表达产物增多[30]。本试验在中国荷斯坦牛第2号染色体CXCR1基因CDS区共发现了13个SNP突变位点。PIGHETTI等[31]报道，CXCR1基因CDS区共有11个SNP位点，也包括4个连锁群。本研究与之相比，333C＞G、337A＞G和642A＞G均为新发现的突变位点，并未发现621A＞G位点，也分为4个连锁群，群内完全连锁。由于621A＞G 为终止密码子突变，会导致蛋白翻译提前终止，说明本研究所检测的中国荷斯坦牛群体中不存在该位点终止密码子突变。333C＞G和337A＞G属于连锁群2，642 A＞G 为独立的连锁群，与其它位点不连锁。徐敏等[27]分析了荷斯坦牛、西门塔尔牛、通江黄牛3个牛品种共160头个体CXCR1基因编码区多态性，在CDS区643、735、819位发现3个SNP位点。官久强等[28,32]分析了470头中国荷斯坦牛CXCR1基因编码区的多态性，结果发现了6个SNPs：735 C＞G 、816 C＞A、819 A＞G、980 A＞G、995 A＞G和1008 C＞T。本研究与上述研究相比较，多发现了4个SNP位点。

3.2 CXCR1基因与临床乳房炎发病次数和生产寿命的关系

本研究所发现的CXCR基因CDS区4个SNP连锁群（642 A＞G、816 A＞C、980A＞G和1068 C＞T）内，CXCR1-642位点与2胎奶牛发生乳房炎次数有显著相关（P＜0.05）。由于CXCR1-642位点为首次发现的位点，并且CXCR1-642位点与其他胎次临床乳房炎发生次数无显著相关（P＞0.05），所以该结果有待进一步探索。CXCR1-980和CXCR1-1068位点与临床乳房炎发病次数无显著相关（P＞0.05），同时这两个位点与奶牛生存时间也无显著相关。CXCR1-816位点与3胎临床乳房炎发病次数有显著相关（P＜0.05），AA基因型个体临床乳房炎发病次数显著低于AC和CC基因型个体。多因素方差分析和生存分析也表明：CXCR1-816位点与奶牛生产寿命有显著相关，CC型个体生存时间平均为63.68月，极显著低于AA型个体（65.5月），CC型比AA型奶牛提前淘汰2个月左右。可见临床乳房炎发病次数较少的基因型个体存活时间较长。综合以上3方面统计数据所得结果，发现CXCR1-816位点与3胎临床乳房炎发病次数及奶牛生产月龄、离群月龄和离群胎次均有显著相关，推测该位点及其连锁群的突变可能对其基因表达产生影响，从而影响乳房炎发病次数、生产寿命等相关指标，但这一推论有待进一步通过实验进行验证。因此，在种公牛站进行遗传评定工作中，增加CXCR1-816基因型检测，筛选其中AA型个体公牛并留种，逐渐淘汰后代中CC型母牛，可降低奶牛群体中临床乳房炎发病次率，延长生产寿命，从而提高奶牛养殖的经济效益。

在本研究中，816 A＞C、291 C＞T、1008 C＞T及3′端792 G＞A和1079 C＞T处于第三连锁群内[31]，虽然CXCR1基因CDS区上述3个位点均为同义突变，没有造成AA改变。而3′端792 G＞A和1079C＞T 2个突变位点，经RNAhybrid软件进行预测，均没有发现miRNA结合位点。由于本研究并未对CXCR1基因5′端启动子区SNP突变及其与CDS区SNP进行连锁分析，在该基因5′端启动子区是否存在与CDS区816 A＞C完全连锁的突变并不得而知，也未见相关报道。因此，816 A＞C与奶牛临床乳房炎发病次数及生产寿命显著关联的遗传机制有待进一步深入研究。

另外，本研究也发现同一基因型在不同胎次间表现不一致，如CXCR1-816位点仅与3胎临床乳房炎发病次数有显著相关，而与1胎和2胎奶牛临床乳房炎发病次数无显著相关。造成该结果可能有以下几个原因：一是由于奶牛年青时抗病力较强，临床乳房炎等疾病性状发生率较低，从而造成统计结果不显著；二是临床乳房炎本身的度量问题。在生产中临床乳房炎常以可以看到明显的乳房红肿、触摸有疼痛感和硬块或乳样呈豆腐样等症状为判断条件，并伴随着乳中体细胞数上升和血液中某些免疫因子发生变化。经转换后的体细胞数（也称为体细胞评分）与临床乳房炎有较高的遗传相关[5]，且在一个泌乳期中有多次记录，增加了结果判断的准确性。第三个原因是引起乳房炎的病原菌类型较为复杂。不同病原菌导致临床乳房炎的表现可能不一致，其发病机制也有所不同，这无疑增加基因型与病原菌间互作的复杂性。

4 结论

CXCR1-816 A＞C突变与中国荷斯坦牛患临床乳房炎次数和生产寿命有显著相关， CXCR1-816 CC型基因型个体有较多的临床乳房炎发病次数和较短的生产寿命，因此被提前淘汰。在进一步验证其功能情况下，该位点可用于中国荷斯坦牛生产寿命的分子标记辅助选择。

[1] MIDDLETON J R, SAEMAN A, FOX L K, LOMBARD J, HOGAN J S, SMITH K L. The National mastitis council: A global organization for mastitis control and milk quality, 50 years and beyond. Journal of Mammary Gland Biology & Neoplasia, 2014, 19(3-4): 241-251.

[2] 何阳花. CD4 与STAT5b 基因单核苷酸多态性与SCC 及产奶性状的关联分析[D]. 北京: 中国农业大学, 2009.

HE Y H. Association analysis between the single nucleotide polymorphism of CD4 and STAT5b gene and SCC and milk traits[D]. Beijing: China Agricultural University, 2009.(in Chinese)

[3] 李芬，周娟，马云，张立岗. 奶牛生产寿命和产奶量的相关性研究.中国牛业科学, 2006, 32(6):12-14.

LI F, ZHOU J, MA Y, ZHANG L G. Study of association between production life of dairy cattle and its milk yield. China Cattle Science, 2006, 32(6):12-14.(in Chinese)

[4] WOLFOVA M, WOLF J, KVAPILIK J, KICA J. Selection for profit in cattle: I. Economic weights for purebred dairy cattle in the Czech Republic. Journal of Dairy Science, 2007, 90: 2442-2455.

[5] KOECK A, LOKER S, MIGLIOR F, KELTON D F, JAMROZIK J, SCHENKEL F S. Genetic relationships of clinical mastitis, cystic ovaries, and lameness with milk yield and somatic cell score in first-lactation Canadian Holsteins. Journal of Dairy Science, 2014，97:5806-5813.

[6] ZHANG Q, GULDBRANDTSEN B, THOMASEN J R, LUND M S, SAHANA G. Genome-wide association study for longevity with whole-genome sequencing in 3 cattle breeds. Journal of Dairy Science, 2016, 99:7289-7298.

[7] HOCH R C, SCHRAUFSTAETTER I U, COCH RANE C G. In vivo, in vitro, and molecular aspects of interleukin-8 and interleukin-8 receptors. Journal of Laboratory & Clinical Medicine, 1996, 128(2): 134-145.

[8] WOLF M, DELGADO M B, JONES S A, DEWALD B, CLARK-LEWIS, BAQQIOLINI M. Granulocyte chemotactic protein 2 acts via both IL8 receptors, CXCR1 and CXCR2. European Journal of Immunology， 1998, 28(1): 164-170.

[9] PAAPE M, MEHRZAD J, ZHAO X, DETILLEUX J, BURCENICH C. Defense of the bovine mammary gland by polymorpho-nuclear neutrophil leukocytes. Journal of Mammary Gland Biology and Neoplasia, 2002, 7(2): 109-121.

[10] BAGGIOLINI M, DEWALD B, MOSER B. Interleukin-8 and related chemotactic cytokines-CXC and CC chemokines. Advances in Immunology, 1994, 55(55):97-179.

[11] BANNERMAN D D, PAAPE M J, LEE J W, ZHAO X, HOPE J C, RAINARD P. Escherichia coli and Staphylococcus aureus elicit differential innate immune responses following intramammary infection. Clinical & Diagnostic Laboratory Immunology, 2004, 11(3): 463-472.

[12] RIOLLET C, RAINARD P, POUTREL B. Differential induction of complement fragment C5a and inflammatory cytokines during intramammary infections with Escherichia coli and Staphylococcus aureus. Clinical and Diagnostic Laboratory Immunology, 2000, 7: 161-167.

[13] BARBER M R, YANG T J. Chemotactic activities in nonmastitic and mastitic mammary secretions: Presence of interleukin-8 in mastitic but not non-mastitic secretions. Clinical & Diagnostic Laboratory Immunology, 1998, 5(1): 82-86.

[14] LEYVA-BACA I, SCHENKEL F, SHARMA B S, JANSEN G B, AND KARROW N A. Identification of single nucleotide polymorphisms in the bovine CCL2, IL8, CCR2 and IL8RA genes and their association with health and production in Canadian Holsteins. Animal Genetics, 2007, 38(3):198-202.

[15] YOUNGERMAN S M, SAXTON A M AND PIGHETTI G M. Novel single nucleotide polymorphism and haplotypes within the bovine CXCR2 gene. Immunogenetics, 2004, 56: 355-359.

[16] YOUNGERMAN S M, SAXTON A M, OLIVER S P AND PIGHETTI G M. Association of CXCR2 polymorphisms with subclinical and clinical mastitis in dairy cattle. Journal of DairyScience, 2004, 87: 2442-2448.

[17] 周雷，王洪梅， 王长法， 柳楠， 侯明海， 李秋玲，仲跻峰. 中国荷斯坦牛CXCR1基因第二外显子新SNPs与乳腺炎的关联分析.畜牧兽医学报. 2011 , 42(8):1063-1070.

ZHOU L, WANG H M, WANG C F, LIU N, HOU M H, LI Q L，ZHONG J F. New genetic polymorphisms within Exon 2 of Chinese Holstein Cattle CXCR1 gene and its association with mastitis. Acta Veterinaria et Zootechnica Sinica, 2011, 42(8):1063-1070. (in Chinese)

[18] CHAPINAL N, BARRIENTOS A K, KEYSERLINGK M A G, GALO E, AND WEARY D M. Herd-level risk factors for lameness in free-stall farms in the northeastern United States and California. Journal of Dairy Science, 2012, 96(1):318-328.

[19] KHATIB H, HEIFETZ E, DEKKERS J C. Association of the protease inhibitor gene with production traits in Holstein dairy cattle. Journal of Dairy Science, 2005, 88(3):1208-1213.

[20] KHATIB H, SCHUTZKUS V, CHANG Y M, ROSA G J. Pattern of expression of the uterine milk protein gene and its association with productive life in dairy cattle. Journal of Dairy Science, 2007, 90(5):2427-2433.

[21] ASHWELL M S, REXROAD C E, MILLER R H, VANRADEN P M，DA Y. Detection of loci affecting milk production and health traits in an elite US Holstein population using microsatellite markers. Animal Genetics, 1997, 28(3): 216-222.

[22] ASHWELL M S, VAN TASSELL C P. Detection of putative loci affecting milk, health, and type traits in a US Holstein population using 70 microsatellite markers in a genome scan. Journal of Dairy Science, 1999, 82(11): 2497-2502.

[23] HEYEN D W,WELLER J I , RON M, BAND M, BEEVER J E, FELDMESSER E, DA Y, WIGGANS G R, VANRADEN P M, LEWIN H A. A genome scan for QTL influencing milk production and health traits in dairy cattle. Physiological Genomics, 1999, 1(3):165-175.

[24] KOMISAREK J, DORYNEK Z. Effect of ABCG2, PPARGC1A, OLR1 and SCD1 gene polymorphism on estimated breeding values for functional and production traits in Polish Holstein-Friesian bulls . Journal of Applied Genetics, 2009, 50: 125-132.

[25] COLE J B, GEORGE R, WIGGANS L M, SONSTEGARD T S, LAWLORJR T J, CROOKER B A, VAN TASSELL C P, YANG J, WANG S W, MATUKUMALLI L K, YANG D. Genome-wide association analysis of thirty one production, health, reproduction and body conformation traits in contemporary U. S. Holstein cows. BMC Genomics, 2011, 12(408): 1471-2164.

[26] CHEBEL R C, SUSCA F, SANTOS J E. Leptin genotype is associated with lactation performance and health of Holstein cows. Journal of Dairy Science, 2008, 91(7):2893-2900.

[27] 徐敏, 平富强, 陈仕毅, 赖松家, 刘益平. CXCR2 基因多态性与奶牛乳房炎和乳品质的关联. 遗传, 2008, 30: 463-468.

XU M, PING F Q, CHEN S Y, LAI S J, LIU Y P. Study on the relationships between polymorphisms of CXCR2 gene and milk quality and mastitis of dairy cow. Herditas, 2008, 30:463-468.(in Chinese)

[28] 官久强, 王洪梅, 王长法, 李秋玲, 李建斌, 帅素容, 侯明海, 仲跻峰. 中国荷斯坦牛白介素8 受体基因编码区多态性与乳腺炎的关联分析. 中国农业科学, 2010, 43(5): 1057-1065.

GUAN J Q, WANG H M, WANG C F, LI Q L, LI J B, SHUAI S R, HOU M H, ZHONG J F. Genetic Polymorphisms Within the Coding Regions of IL8R Gene and Its Association with Mastitis Trait in Chinese Holstein Cattle. Scientia Agricultura Sinica, 2010, 43(5): 1057-1065. (in Chinese)

[29] SHI Y Y, HE L. SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Research, 2005, 15(2):97-98.

[30] LEE J W, BANNERMAN D D, PAAPE M J, HUANG M K, ZHAO X. Characterization of cytokine expression in milk somatic cells during intramammary infections with Escherichia coli or Staphylococcus aureus by real-time PCR. Veterinary Research, 2006, 37(2):219-229.

[31] PIGHETTI G M, KOJIMA C J, WOJAKIEWICZ L, RAMBEAUD M. The bovine CXCR1 gene is highly polymorphic. Veterinary Immunology and Immunopathology, 2011, 145 (1-2): 464-470.

[32] 官久强, 王洪梅, 王长法, 李秋玲, 仲跻峰, 侯明海, 帅素容. 中国荷斯坦牛CXCR1基因编码区的遗传多态性. 西北农林科技大学学报(自然科学版), 2009, 37(6): 47-52.

GUAN J Q, WANG H M, WANG C F, LI Q L, ZHONG J F, HOU M H, SHUAI S R. Genetic polymorphisms within the coding regions of CXCR1 gene in Chinese Holstein cattle. Journal of Northwest A &F University (Nat. Sci. Ed .), 2009, 37(6): 47-52. (in Chinese)

（责任编辑 林鉴非）

Association Between SNPs in the CDS Regions of CXCR1 Gene and the Clinical Mastitis and Lifetime for Chinese Holstein

WANG MengQi1, NI Wei1, ZHANG HuiMin1, YANG ZhangPing1, WANG XiPu2, JIANG YanSen2, MAO YongJiang1,2
(1College of Animal Science and Technology, Yangzhou University, Yangzhou 225009;2Jiangsu Junlebao Dairy Co., Ltd., Xuzhou 221721, Jiangsu)

【Objective】 The objective of this study was to investigatethe association between the SNPs in the coding sequence (CDS) of CXCR1 gene and the clinical mastitis and lifetime for Chinese Holstein.【Method】SNPs in the CDS of CXCR1 gene were screened by using PCR and direct sequencing for 20 cows with low SCS and 20 cows with high SCS. Finally, the selected SNPs of CXCR1 gene for 866 Chinese Holstein cows were detected using flight mass spectrometry. The clinical mastitisand productive life of tested cows were collected from the management system of dairy farm. The association between the SNPs and the clinical mastitis and productive life of tested cows was analyzed using multi factor variance analysis, Logistic regression, Cox regression.【Result】A total of 13 SNPs were found in the CDS region of CXCR1 gene (291 C＞T, 333 C＞G, 337 A＞G, 365 C＞T, 570 A＞G, 642 A＞G, 735 C＞G, 816 A＞C, 819 A＞G, 980 A＞G、, 995 A＞G, 1008 C＞T and 1068 A＞G), and these SNPs were divided into 4 linkage groups，and four SNPs were selected from each linkage group for further analysis（642 A＞G, 816 A＞C, 980A＞G and 1068 A＞G). A total of 9 haplotypes were observed for 4 SNPs of CXCR1 gene, the GAGG haplotype frequency was the highest (0.3141), and the frequency of haplotype ACAA was the lowest (0.0017). CXCR1-642 A＞G showed a significant association with the number of dairy cows suffering from clinical mastitis in second lactation (P＜0.05). The number of dairy cows suffering from clinical mastitis for individual with AG genotype was significantly higher than that of AA genotype (P＜0.05). The number of dairy cows suffering from clinical mastitis for individual with CXCR1-816 AA genotype was significantly higher than that of AC and CC genotypes (P＜0.05). The other SNPs showed no significant association with the number of clinical mastitis. CXCR1-816 A＞C mutation showed an extremely significant association with culling age (P＜0.01) and a significant association with productive month and culling lactation (P＜0.05). The productive month, culling age and culling lactation for the individual of CXCR1-816 AA were significantly higher than that of CC genotype (P＜0.05). Cox regression showed that only CXCR1-816 showed significant association with the productive life（P＜0.05), the survival probability of individual with CC genotype was significantly lower than that of AA and AC genotype in each time period. 【Conclusion】The mutation of CXCR1-816 A＞C showed significant association with the clinical mastitis and lifetime for Chinese Holstein, and this SNP could be used for molecular marker-assisted selection of productive life for Chinese Holstein when the molecular function of this mutation were proved.

Chinese Holstein; CXCR1; SNPs; clinical mastitis; productive lifetime

2016-07-25；接受日期：2017-04-01

国家自然科学基金（31372286）、江苏省优势学科（PAPD）、江苏省企业研究生工作站资助

联系方式：王梦琦, E-mail：770406499@qq.com 。通信作者毛永江，E-mail：cattle@yzu.edu.cn