孟 瑤 ，范業(yè)鍇 ，張偉峰 ，何 旭 ，張軍輝 *

（1.河北工程大學(xué)，河北 邯鄲 056000；2.河北省畜牧良種工作總站，河北 石家莊 050000）

想要更好地了解和研究豬的繁殖性能，以此進(jìn)一步降低成本及提高生產(chǎn)效率，需要精確的數(shù)據(jù)支撐以及對(duì)數(shù)據(jù)進(jìn)行全方位的整理和記錄，為研究豬的繁殖性能和之后的相關(guān)數(shù)據(jù)計(jì)算打下堅(jiān)實(shí)的基礎(chǔ)，也是為了讓最終的數(shù)據(jù)更加的科學(xué)和嚴(yán)謹(jǐn)以及準(zhǔn)確。相比國(guó)外的一些進(jìn)展，美系等豬種的數(shù)據(jù)，其中都有矯正等規(guī)范的計(jì)算過(guò)程，也方便后期與此相應(yīng)的數(shù)據(jù)準(zhǔn)備。但就目前國(guó)內(nèi)大約克夏種豬在市場(chǎng)中的數(shù)據(jù)而言，雖然養(yǎng)殖的人數(shù)不少，但并沒(méi)有成型的數(shù)據(jù)庫(kù)，或者是沒(méi)有人將這些數(shù)據(jù)進(jìn)行歸納整合，使其形成一個(gè)完整的遺傳關(guān)系，基于上述原因，對(duì)于這方面的研究以及數(shù)據(jù)的歸納整合就少些。市場(chǎng)上的種豬多以杜洛克豬、長(zhǎng)白豬、大約克夏豬3個(gè)品種為主，此次試驗(yàn)將大約克夏豬作為研究對(duì)象，計(jì)算其相關(guān)的遺傳參數(shù)估計(jì)。

1 材料與方法

1.1 材料的收集和整理

依據(jù)試驗(yàn)所需要求，以河北省某省級(jí)原種豬場(chǎng)中的5個(gè)種豬場(chǎng)為采樣點(diǎn)，取自2004年到2021年的大約克夏種豬7 011頭，33 895條產(chǎn)仔記錄為樣本數(shù)據(jù)。其中大約克夏豬的繁殖性能的數(shù)據(jù)包括窩總產(chǎn)仔數(shù)、活仔數(shù)、健仔數(shù)、斷奶數(shù)等性狀。

1.2 數(shù)據(jù)分析與方法

首先，利用Visual Foxpro軟件將已經(jīng)拿到的繁殖性能的數(shù)據(jù)進(jìn)行初步的整理，使用此軟件，其目的也是為了減少工作量，可以在處理和整合數(shù)據(jù)的過(guò)程中把不符合要求的異常記錄進(jìn)行刪除，避免個(gè)別無(wú)效數(shù)據(jù)給整體數(shù)據(jù)帶來(lái)的錯(cuò)誤影響，并標(biāo)準(zhǔn)化所有的格式。之后，利用MTDFREML軟件將數(shù)據(jù)單獨(dú)分為各個(gè)性狀，利用不同性狀之間的數(shù)據(jù)進(jìn)行數(shù)據(jù)處理，得出相關(guān)的數(shù)據(jù)結(jié)果，在此過(guò)程中主要是根據(jù)MTDFNRM.EXE、MTDFPREP.EXE、MTDFRUN.EXE這幾部分的數(shù)據(jù)模塊，以得到遺傳參數(shù)，也就是MTDF76號(hào)文件。過(guò)程中也利用了BLUP動(dòng)物模型，即yijkl=μ+hysi+lj+ak+pel+eijkl

其中：

yijkl：總產(chǎn)仔數(shù)、產(chǎn)活仔數(shù)、健仔數(shù)的觀察值

μ：總平均數(shù)

hysi：母豬產(chǎn)仔時(shí)場(chǎng)年季固定效應(yīng)

lj：母豬出生的窩效應(yīng)，服從（0，Vl）分布

ak：個(gè)體的隨機(jī)遺傳效應(yīng)，服從（0，Va）分布

pel：母豬永久環(huán)境效應(yīng)，服從（0，Vpe）分布

eijkl：隨機(jī)剩余效應(yīng)，服從（0，Ve）分布

2 結(jié)果

2.1 近交系數(shù)

近交系數(shù)由原始數(shù)據(jù)來(lái)進(jìn)行計(jì)算，其中有個(gè)體號(hào)、場(chǎng)次、出生的編號(hào)等一系列數(shù)據(jù)，由于軟件繪制的數(shù)據(jù)的結(jié)果，所計(jì)算出來(lái)的近交系數(shù)近似為0，也表明了豬群規(guī)避開(kāi)了近親選配。

2.2 遺傳參數(shù)估計(jì)

表1為遺傳參數(shù)估計(jì)結(jié)果。

表1 遺傳參數(shù)估計(jì)的結(jié)果

3 討論

3.1 遺傳參數(shù)估計(jì)

表1顯示了各自的遺傳力所對(duì)應(yīng)的數(shù)據(jù)，根據(jù)表格兩兩對(duì)應(yīng)，總產(chǎn)仔數(shù)和活仔數(shù)的遺傳相關(guān)在0.98，即二者具有高正相關(guān)的關(guān)系；同理，總產(chǎn)仔數(shù)與健仔數(shù)的遺傳相關(guān)為0.96，活仔數(shù)和健仔數(shù)的遺傳相關(guān)為0.98，這兩組兩兩對(duì)應(yīng)中，每一組的兩者也都具備高的正相關(guān)。此外，在斷奶數(shù)與總產(chǎn)仔數(shù)、活仔數(shù)、健仔數(shù)這三個(gè)繁殖性狀的兩兩組合中，遺傳相關(guān)分別為0.95、0.98、0.98，表明斷奶數(shù)與其他三個(gè)繁殖性能性狀呈現(xiàn)正相關(guān)，即前者的數(shù)據(jù)越大，后者數(shù)據(jù)也會(huì)隨之變得越大。除此之外，計(jì)算過(guò)程中也可以看出不同場(chǎng)次、胎次、季節(jié)等因素對(duì)于大約克夏豬不同繁殖性狀的影響，以及相應(yīng)的個(gè)體EBV等數(shù)據(jù)。

4 結(jié)論

本試驗(yàn)大約克夏豬總產(chǎn)仔數(shù)、活仔數(shù)、健仔數(shù)、斷奶數(shù)的遺傳力分別為0.22、0.2、0.19和0.17?？偖a(chǎn)仔數(shù)和活仔數(shù)、健仔數(shù)、斷奶數(shù)的遺傳相關(guān)系數(shù)分別為0.98、0.96、0.95；活仔數(shù)和健仔數(shù)、斷奶數(shù)的遺傳相關(guān)系數(shù)為0.98、0.98；斷奶數(shù)與健仔數(shù)的遺傳相關(guān)系數(shù)為0.98。