随着非洲猪瘟出现，猪的生长性状（达100 kg体重日龄和背膘厚度）提升不明显，总产仔数从2016年的12.31头增加至2021年的13.15头。近年来，改善农场环境、生物安全措施以及重视育种选择均有利于提高猪繁殖性能[1]。但高产母猪增加了仔猪死亡率[2-3]，仍是养猪生产中面临的重要问题[4-5]。此外，产仔数增加会降低仔猪的初生重，增加仔猪与功能性乳头的比率，延长分娩时间等[6-7]，这些变化可能会对母猪和仔猪的健康和福利造成不利影响。从育种角度分析，将生存性状纳入遗传选择目标可以促进更多有活力的仔猪出生[8]。但仍存在一些问题，如对生存性状选择可能会对某些有益性状产生负面影响[9]；而且母猪繁殖性状间具有很强的相关性，增加了选育难度[10]。多元统计分析（主成分分析和聚类分析）为解决上述问题提供了参考，且已有研究使用综合评分模型和指标聚类图确定猪繁殖性能的选育指标[11-13]，但随着母猪产仔数提高，仍需对大白猪繁殖性状进行研究。繁殖性状的选择标准主要包括总产仔数和产活仔数[14]，对其他重要性状的关注较少。本文整合大白母猪的表型，通过多元统计分析对繁殖性状进行评估，获得能够客观代表繁殖性能的重要指标，为大白猪繁殖性状的选育提供参考。1材料与方法1.1试验动物大白母猪由河北省衡水市种猪场提供。选择2020年7月至10月分娩的244头大白母猪（与配公猪为大白）作为研究对象（1胎38头，2胎146头，3胎60头）。所用繁殖性状数据经GBS系统导出，生产信息由场内育种专员记录，猪个体记录清晰、系谱完整。1.2试验设计本研究繁殖性状主要包括总产仔数（TNB）、产活仔数（NBA）、健仔数（NHB）、初生窝重（TBW）、平均初生重（ABW）、断奶仔数（WN）、断奶窝重（TW）、平均断奶重（AWW）、存活率（SR）、产程（FD）、出生间隔（PBI）、有效乳头数（ETN）。其中健仔数指初生体重不小于0.8 kg、无遗传缺陷、活力较好的仔猪头数。仔猪21日龄记录断奶重。1.3数据统计与分析采用SPSS 22.0软件对12个繁殖性状进行主成分分析，以提取在总变异中保持尽可能多信息的主成分（主成分特征值1，贡献率85.0%[15]；有效乳头数在主成分分析中未被考虑）。计算特征向量[成分矩阵中的每个值（主成分1、2、3、4）分别除以相应特征值的平方根]。采用组间连接法以平方欧氏距离为测度，对具有相似性状的变量进行R型聚类分析。计算各个聚类中每一个性状与其余性状之间的相关指数均值（R2¯）[16]。R2¯=∑r2m-1 (1)式中：r为Spearman相关系数；m为变量（性状）的数量；R2¯最大的性状作为此类的典型性状。使用Origin 2021的Correlation Plot分析繁殖性状间的Spearman相关。“*”为FDR校正P0.05，“**”为FDR校正P0.01。2结果与分析2.1母猪繁殖性状描述统计（见表1）由表1可知，母猪各繁殖性状的基本统计量包括数据量、平均数、中位数、标准差、最小值、最大值和变异系数。ETN变异系数最低，为8.24%；其次是ABW和SR，均低于15.54%。FD、PBI、TW和NBA/ETN变异系数相对较高，均高于23%，表明这些性状在个体间差异较大。10.13557/j.cnki.issn1002-2813.2023.08.013.T001表1母猪繁殖性状描述统计性状数据量/头平均值中位数标准差最小值最大值变异系数/%TNB/(头/窝)24417.0217.003.484.0027.0020.44NBA/(头/窝)24414.7215.003.342.0024.0022.71NHB/(头/窝)24413.9514.003.212.0024.0023.01TBW/(kg/窝)24419.9220.094.372.4631.9721.95ABW/(kg/头)2441.371.370.200.832.4014.87WN/(头/窝)24411.8912.002.045.0017.0017.15TW/(kg/窝)24461.4663.0914.4017.8191.6623.44AWW/(kg/头)2445.185.270.942.488.3418.16SR2440.870.920.130.351.0015.54FD/min244328.73274.00205.9287.001 300.0062.64PBI/min24419.9615.6712.575.3070.3663.14ETN/(个/头)23913.7914.001.149.0016.008.24NBA/ETN2391.071.080.260.171.7523.922.2繁殖性状的相关性分析（见图1）由图1可知，TW与AWW（r=0.66）、WN（r=0.62）、SR（r=0.53）、ABW（r=0.23）、TBW（r=0.19）呈正相关（P0.01）。SR与TNB（r=-0.18）、NBA（r=-0.17）呈负相关（P0.05），与TW（r=0.53）、WN（r=0.45）、AWW（r=0.27）、ABW（r=0.26）呈正相关（P0.01）。ETN仅与NBA/ETN相关（r=-0.37，P0.01）。10.13557/j.cnki.issn1002-2813.2023.08.013.F001图1繁殖性状相关性分析由于ETN与大部分繁殖性状无显著相关性，故未将其纳入主成分分析。2.3繁殖性状的主成分分析对繁殖性状的12个指标进行主成分分析，各主成分的特征值、贡献率和累积贡献率见表2。10.13557/j.cnki.issn1002-2813.2023.08.013.T002表2繁殖性状的总方差解释和成分矩阵项目PC1PC2PC3PC4特征值4.5142.7401.7311.212贡献率/%37.61722.83114.42910.103累积贡献率/%37.61760.44974.87884.980NBA0.977-0.0570.106-0.010NHB0.9460.0520.1180.128NBA/ETN0.931-0.0830.1120.028TNB0.813-0.2200.226-0.142TBW0.8060.2720.1440.404TW0.0240.8900.289-0.111SR-0.0860.8010.086-0.272WN0.2180.6820.039-0.589AWW-0.1830.5160.3640.425FD-0.175-0.3750.888-0.115PBI-0.493-0.2780.778-0.071ABW-0.3530.4710.0610.615由表2可知，4个PC的特征值大于1，分别为4.514、2.740、1.731和1.212，贡献率分别为37.617%、22.831%、14.429%和10.103%。前4个主成分的累计贡献率达84.980%，可以反映原变量信息。4个主成分的主成分系数见表3。10.13557/j.cnki.issn1002-2813.2023.08.013.T003表34个主成分的主成分系数项目PC1PC2PC3PC4NBA0.460-0.0340.081-0.009NHB0.4450.0310.0890.116NBA/ETN0.438-0.0500.0850.026TNB0.383-0.1330.172-0.129TBW0.3790.1640.1090.367TW0.0110.5380.220-0.100SR-0.0400.4840.066-0.247WN0.1030.4120.029-0.535AWW-0.0860.3120.2760.386FD-0.083-0.2270.675-0.104PBI-0.232-0.1680.591-0.065ABW-0.1660.2840.0460.559为评价猪的繁殖性能，根据主成分系数表和性状标准化值计算相应的F值。F1=0.460×NBA+0.445×NHB+0.438×NBA/ETN+0.383×TNB+0.379×TBW+0.011×TW-0.040×SR+0.103×WN-0.086×AWW-0.083×FD-0.231×PBI-0.166×ABW。同样方法计算F2、F3和F4。以4个主成分的贡献率为权重得到主成分的综合得分，F=37.617%F1+22.831%F2+14.429%F3+10.103%F4，各指标的综合得分为F=0.176 NBA+0.199 NHB+0.168 NBA/ETN+0.125 TNB+0.233 TBW+0.149 TW+0.080 SR+0.083 WN+0.118 AWW+0.004 FD-0.047 PBI+0.065 ABW。2.4繁殖性状聚类分析（见图2）由图2可知，取欧式距离为16.0时可将13个性状分为5类，即NBA、NHB、NBA/ETN、TBW和TNB聚为一类；WN、TW和SR聚为一类；FD和PBI聚为一类；ABW和AWW聚为一类；ETN单独为一类。每个类别中最大的R2¯作为此类的典型性状，分别为NBA、TW、FD、PBI、AWW、ABW和ETN。10.13557/j.cnki.issn1002-2813.2023.08.013.F002图2繁殖性状聚类分析（平方欧式距离法）3讨论据报道，随着每头母猪年提供断奶仔猪数（PSY）提高，约20%的初生仔猪在断奶前死亡[17]，本研究发现，仔猪断奶前死亡率高（13%）。仔猪高死亡率给生产和福利均带来不利影响。本研究相关性分析表明，繁殖性状间存在不同程度的相关性，对一个性状的选择必须经过全面衡量。主成分分析和聚类分析可对复杂的性状作降维处理，以获得一些能够客观代表繁殖性能的重要性状[18]。3.1根据主成分分析和聚类分析确定选育指标分析本研究中，NBA（特征向量0.460）、NHB（特征向量0.445）、NBA/ETN（特征向量0.438）、TNB（特征向量0.383）和TBW（特征向量0.379）的方差在PC1中最大，反映了仔猪出生相关表现。PC2中，TW（特征向量0.538）、SR（特征向量0.484）、WN（特征向量0.412）和AWW（特征向量0.312）的方差最大，反映了仔猪断奶表型。前两个主成分（PC1和PC2）解释了总方差的60.45%，代表了大部分指标信息；后两个主成分（PC3和PC4）累积贡献率分别为14.43%和10.10%，占比较小，更多体现母猪产程及仔猪个体重等信息。结合聚类分析结果，NBA和TW可作为影响母猪繁殖力的主要性状。由于典型性状的判断需要两个以上指标，因此后两类性状可结合主成分系数判断重要性，即ABW和PBI所占权重更大，因此在选育时应尽量侧重这些性状。张力等[12]利用主成分分析确定了长白母猪繁殖性状育种的重要指标，并指出将TW、ABW、SR和ETN纳入育种选择。Kramarenko等[13]评估了不同基因组间猪的繁殖性能，发现PC1以TNB、NBA和TBW为代表，PC2以WN、TW和AWW为代表，与本研究结果类似。3.2繁殖性状相关性分析本研究相关性分析表明，NBA与ABW之间呈显著负相关（r=-0.38），ABW与SR之间呈显著正相关（r=0.26），表明产活仔数越多，ABW越低，仔猪SR也越低。SR和TW之间存在显著正相关（r=0.53），TW和NBA之间无显著相关（r=-0.02）。大白猪的NBA和TW之间存在正遗传相关（r=0.45）[19]，表明选择TW不会对产活仔数产生负面影响。TW（约20 d）可反映母猪断奶前的泌乳能力[20]，其结合了窝产仔数和仔猪21日龄体重信息，取决于初生仔猪的数量和断奶前仔猪的存活与生长情况[21]。Dube等[22]指出，NBA和TW是提高母猪生产力最合适的指标。更多的产活仔数量和更高的21日龄窝重将有利于提高母猪繁殖力。因此，推测TW可能是一个潜在提高猪繁殖性能的性状且对仔猪存活有益。主成分综合得分可以作为繁殖性能的综合评价指标。本研究中，正系数的性状主要是NBA、NHB、NBA/ETN、TNB、TBW、TW和AWW，表明对母猪繁殖力具有积极作用，而PBI（系数为负）则不利于母猪繁殖力的提高。TNB的增加可能会增加PBI，对死胎率和哺乳期前24 h及断奶前的仔猪存活率产生负面影响[7,23]。原因可能是较长FD导致仔猪出生时处于缺氧状态，从而降低了仔猪吸吮反应，仔猪温度平衡性差[23]。本研究中，FD、PBI与SR无显著相关，表明其对繁殖力的贡献不大。ABW是与仔猪存活相关的重要性状，直接影响仔猪体温调节能力和生长。Chris等[24]发现，ABW大于1.8 kg的仔猪具有较高的存活率（90%），而ABW为700 g的仔猪存活率较低（33%）。因此，提高仔猪ABW有利于断奶前仔猪的存活。4结论本研究基于多元统计分析对大白猪繁殖性状进行分析，发现NBA和TW是大白猪繁殖性能的重要指标，也应关注ABW和PBI这两个性状，尤其在大窝仔猪中。分娩间隔过长会给仔猪带来的不利影响。高产母猪系的选择使产仔数增加，但同时仔猪死亡率升高也对生产和福利产生负面影响，提示选育时要充分考虑其他性状的影响。