張鎖宇，邱小田，丁向東，張　勤*，王志剛，王長存，屈云龍，

(1.中國農(nóng)業(yè)大學(xué)動(dòng)物科技學(xué)院，畜禽育種國家工程實(shí)驗(yàn)室，農(nóng)業(yè)部動(dòng)物遺傳育種與繁殖重點(diǎn)實(shí)驗(yàn)室，北京 100193；2.全國畜牧總站，北京 100125； 3.重慶市畜牧科學(xué)院，重慶 402460)

?

利用關(guān)聯(lián)組估計(jì)中國大白、長白及杜洛克豬繁殖性狀的遺傳參數(shù)

張鎖宇1，2#，邱小田2，3#，丁向東1，張勤1*，王志剛2，王長存2，屈云龍1，2

(1.中國農(nóng)業(yè)大學(xué)動(dòng)物科技學(xué)院，畜禽育種國家工程實(shí)驗(yàn)室，農(nóng)業(yè)部動(dòng)物遺傳育種與繁殖重點(diǎn)實(shí)驗(yàn)室，北京 100193；2.全國畜牧總站，北京 100125； 3.重慶市畜牧科學(xué)院，重慶 402460)

摘要：本研究旨在分析大白、長白及杜洛克3個(gè)品種豬在不同聯(lián)合遺傳評(píng)估組中總產(chǎn)仔數(shù)的遺傳參數(shù)，為將來全國豬聯(lián)合遺傳評(píng)估奠定基礎(chǔ)。利用DMU軟件的DMU4過程計(jì)算74個(gè)國家生豬核心育種場間的關(guān)聯(lián)率，根據(jù)聯(lián)合遺傳評(píng)估組內(nèi)各場相互存在關(guān)聯(lián)且每個(gè)場的平均關(guān)聯(lián)率在3%以上原則，每個(gè)品種篩選出3個(gè)聯(lián)合遺傳評(píng)估組(組1、組2和組3)。利用DMU軟件和動(dòng)物模型估計(jì)3個(gè)聯(lián)合遺傳評(píng)估組中大白、長白及杜洛克豬總產(chǎn)仔數(shù)(TNB)性狀的加性遺傳、永久環(huán)境效應(yīng)方差，分析時(shí)同時(shí)考慮管理組的固定效應(yīng)，其中大白豬在組1、2和3中的繁殖記錄數(shù)分別為61 161、19 974和18 002條；長白豬分別為28 239、11 918和3 064條；杜洛克豬分別為7 174、6 512和2 417條。結(jié)果表明，場間關(guān)聯(lián)總體情況較差(平均關(guān)聯(lián)率為0.20%)，全國范圍的聯(lián)合遺傳評(píng)估條件還不滿足，但可根據(jù)本研究結(jié)果進(jìn)行局部聯(lián)合遺傳評(píng)估。永久環(huán)境效應(yīng)估計(jì)值的變化為0.05～0.11，不同組及不同品種的總產(chǎn)仔數(shù)遺傳力估計(jì)值約為0.10，在文獻(xiàn)推薦范圍之內(nèi)。結(jié)果提示，組1中各品種豬的遺傳力估計(jì)結(jié)果要優(yōu)于其他聯(lián)合遺傳評(píng)估組，因此將來進(jìn)行全國豬聯(lián)合遺傳評(píng)估時(shí)，可參考組1中各品種的TNB遺傳力估計(jì)結(jié)果。

關(guān)鍵詞：豬；遺傳參數(shù)；總產(chǎn)仔數(shù)；關(guān)聯(lián)率

養(yǎng)豬生產(chǎn)中，對窩性狀的選擇更吸引育種者，因?yàn)樯L性狀經(jīng)過多年持續(xù)的選擇可能已接近最適值[1]。在窩性狀中產(chǎn)仔數(shù)性狀尤其重要，因?yàn)槊吭黾?頭仔豬，將會(huì)對經(jīng)濟(jì)效益產(chǎn)生重大的影響。加拿大豬遺傳改良中心(CCSI)母系指數(shù)中總產(chǎn)仔數(shù)的權(quán)重超過了50%[2]，丹麥母系豬的育種目標(biāo)中產(chǎn)仔數(shù)性狀(5日齡仔豬數(shù))在2005、2010和2011年的權(quán)重分別為70%、37%和27%[3]。

育種值及遺傳參數(shù)的估計(jì)是家畜育種工作的一項(xiàng)中心任務(wù)。家畜的估計(jì)育種值是選擇種用家畜最主要的依據(jù)，對家畜育種值估計(jì)的準(zhǔn)確性會(huì)直接影響畜群的遺傳進(jìn)展。畜群的遺傳參數(shù)除為估計(jì)育種值所必需外，也是制定育種規(guī)劃、了解數(shù)量性狀遺傳機(jī)制等所不可缺少的[4]。生產(chǎn)者則可以通過對豬群重要經(jīng)濟(jì)性狀遺傳參數(shù)的了解進(jìn)而從群體的遺傳進(jìn)展中獲取經(jīng)濟(jì)效益，這些遺傳參數(shù)還可以幫助生產(chǎn)者預(yù)測選擇反應(yīng)及遺傳趨勢的過程[5]。從統(tǒng)計(jì)學(xué)的角度來講，遺傳參數(shù)的估計(jì)就是方差組分的估計(jì)。目前在畜禽育種中應(yīng)用比較多的是 ML、 REML 和 Bayesian 方法，總體來看，ML 法精確性高，準(zhǔn)確性差，計(jì)算難度小，而 REML 法精確性低，準(zhǔn)確性好，計(jì)算難度大，在家畜育種中 REML 法似乎比 ML 法更受人偏愛一些[4]。

目前我國豬的遺傳評(píng)估工作還處于起步階段，數(shù)據(jù)的積累還很有限，不同場間的關(guān)聯(lián)程度還較低，無法跨場進(jìn)行全國范圍的聯(lián)合遺傳評(píng)估，2004年以來，在全國畜牧總站、農(nóng)業(yè)部948重大專項(xiàng)等項(xiàng)目支持下，依托中國農(nóng)業(yè)大學(xué)的技術(shù)力量成立了全國種豬遺傳評(píng)估中心(以下簡稱中心)。中心主要職能為收集、整理和分析全國種豬登記和性能測定數(shù)據(jù)，組建“國家種豬數(shù)據(jù)庫”，為《全國生豬遺傳改良計(jì)劃(2009-2020)》的順利實(shí)施提供技術(shù)支持和智力保障，并致力于我國種豬聯(lián)合育種平臺(tái)網(wǎng)絡(luò)化建設(shè)，建立了全國種豬遺傳評(píng)估信息網(wǎng)。

中心的成立使得跨場進(jìn)行遺傳評(píng)估成為了可能，而跨場進(jìn)行遺傳評(píng)估的基本前提是各場之間要有一定程度的關(guān)聯(lián)性，P.Mathur等[6]研究表明，要進(jìn)行場間的聯(lián)合遺傳評(píng)估，建議各個(gè)場的平均關(guān)聯(lián)率大于3%，之前基于國內(nèi)群體的遺傳參數(shù)估計(jì)多針對單一育種場[7-9]，群體規(guī)模有限，而遺傳參數(shù)估計(jì)的準(zhǔn)確性與群體結(jié)構(gòu)、分析方法及統(tǒng)計(jì)模型緊密相關(guān)[4，10]，群體規(guī)模顯得尤為重要。同一聯(lián)合遺傳評(píng)估組可近似看做同一群體，可使群體規(guī)模大大增加。目前中心采用的遺傳參數(shù)是加拿大豬遺傳改良中心(CCSI)于2000年之前制定的，并未做品種的區(qū)分。由于遺傳參數(shù)具有群體特異性[11]，參考的遺傳參數(shù)不一定適合中國，更不一定適合每個(gè)場，有必要基于我國種豬群體進(jìn)行遺傳參數(shù)的估計(jì)?；谝陨夏康?，本研究利用關(guān)聯(lián)組估計(jì)我國大白、長白及杜洛克豬群體總產(chǎn)仔數(shù)(TNB)的遺傳參數(shù)，為將來進(jìn)行全國豬聯(lián)合遺傳評(píng)估奠定基礎(chǔ)。

1材料與方法

1.1數(shù)據(jù)

所用數(shù)據(jù)來自74個(gè)國家生豬核心育種場于2014年3月20日前(含3月20日)上傳至“全國種豬遺傳評(píng)估信息網(wǎng)”的種豬登記和繁殖數(shù)據(jù)。杜洛克、長白和大白豬3個(gè)品種繁殖母豬數(shù)分別為32 327、77 533和164 315頭，出生于2007～2013年間。母豬產(chǎn)仔記錄數(shù)共計(jì)81萬多條(圖1)。主要繁殖性狀為總產(chǎn)仔數(shù)(TNB)，其為出生時(shí)同窩的仔豬總數(shù)，包括活仔、死胎、木乃伊和畸形豬，更能從遺傳上反映母豬的繁殖能力。

圖1　國家生豬核心育種場數(shù)據(jù)趨勢圖Fig.1　Trends in data for major pig breeding farms in China

1.2場間遺傳聯(lián)系和關(guān)聯(lián)組劃分

鑒于中國目前尚未形成全國范圍的場間遺傳聯(lián)系，本研究首先根據(jù)P.Mathur等[6]提出的場間關(guān)聯(lián)率，計(jì)算74個(gè)國家生豬核心場間的彼此場間聯(lián)系情況，然后將這些場劃分為若干關(guān)聯(lián)組。

P.Mathur等[6]將場間關(guān)聯(lián)定義為群間或者同期組間的相關(guān)，并且建議利用場效應(yīng)之間的相關(guān)來表示關(guān)聯(lián)程度，稱為聯(lián)系率(Connectedness rating，CR)，公式：

據(jù)P.Mathur等[6]研究，要進(jìn)行場間的統(tǒng)一遺傳評(píng)估，各個(gè)場的平均聯(lián)系率應(yīng)大于3%。據(jù)此標(biāo)準(zhǔn)，本研究列出了可進(jìn)行聯(lián)合遺傳評(píng)估的場(表1)，聯(lián)合遺傳評(píng)估組內(nèi)各場之間彼此存在關(guān)聯(lián)且每個(gè)場的平均關(guān)聯(lián)率不小于3%。

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

參照《全國種豬遺傳評(píng)估方案》，總產(chǎn)仔數(shù)遺傳評(píng)估模型：

y=Xb+Za+Wpe+e

遺傳力、永久環(huán)境效應(yīng)及重復(fù)力計(jì)算公式：

本研究首先利用SAS9.2和R語言剔除有系譜和表型記錄錯(cuò)誤的個(gè)體，3個(gè)標(biāo)準(zhǔn)差之外的個(gè)體即視為異常值，予以剔除。場間遺傳聯(lián)系和遺傳參數(shù)估計(jì)利用DMU軟件分析。DMU軟件[13]是一個(gè)用于多元混合模型分析的軟件包，用于估計(jì)正態(tài)分布和非正態(tài)分布性狀的方差-協(xié)方差組分和遺傳參數(shù)估計(jì)，可以配合一般線性模型、隨機(jī)回歸模型、廣義線性模型等。其中場間遺傳聯(lián)系采用DMU4模塊，遺傳參數(shù)利用DMUAI模塊進(jìn)行估計(jì)。DMUAI采用平均信息(Average information，AI)和期望最大化(Expectation maximization，EM)算法相結(jié)合的約束最大似然(Restricted maximum likelihood，REML)方法進(jìn)行遺傳參數(shù)估計(jì)。

2結(jié)果

2.1聯(lián)合遺傳評(píng)估組

本研究共分析了74個(gè)核心育種場的大白豬數(shù)據(jù)，73個(gè)長白豬數(shù)據(jù)及71個(gè)杜洛克豬數(shù)據(jù)，詳細(xì)的關(guān)聯(lián)率情況見國家生豬核心育種場遺傳評(píng)估報(bào)告(2014-01)[14]。從表1可以看出，大白、長白和杜洛克豬均得到3個(gè)聯(lián)合遺傳評(píng)估組，其中大白豬組1、2和3中的育種場數(shù)分別為7、4和4個(gè)；長白豬中的育種場數(shù)分別為7、4和3個(gè)；杜洛克豬中的育種場數(shù)分別為6、4和3個(gè)?？傮w來講，場間關(guān)聯(lián)性較差。表2和表3列出了各品種各聯(lián)合遺傳評(píng)估組及全國的平均關(guān)聯(lián)率情況。表2中各品種聯(lián)合遺傳評(píng)估組內(nèi)豬場的平均關(guān)聯(lián)率均大于3%。本研究得到的各品種聯(lián)合遺傳評(píng)估組的平均關(guān)聯(lián)率均在4%以上，遠(yuǎn)高于全國生豬核心育種場得到的0.20%的平均值。

2.2總產(chǎn)仔數(shù)表型值

杜洛克、長白和大白3個(gè)品種3個(gè)聯(lián)合遺傳評(píng)估組中總產(chǎn)仔數(shù)表型值基本情況見表4。每個(gè)品種都是組1包含的場最多，其包含的有效繁殖記錄數(shù)也最多。以大白為例，組1包含7個(gè)彼此有遺傳關(guān)聯(lián)的場，累計(jì)61 161條總產(chǎn)仔記錄，遠(yuǎn)遠(yuǎn)高于組2和組3(兩個(gè)組都有4個(gè)場)的19 974和18 002條記錄。杜洛克和長白也體現(xiàn)了同樣規(guī)律，只是組3規(guī)模較小，這兩個(gè)品種分別僅有2 417和3 064條記錄。3個(gè)品種中，大白豬記錄數(shù)最多，分別達(dá)到了61 161、19 974和18 002條，合計(jì)99 137條，杜洛克豬記錄數(shù)最少，分別為7 174、6 512和2 417條，合計(jì)16 103條。這與我國飼養(yǎng)大白豬最多的實(shí)際情況一致，而杜洛克作為終端父本，純種飼養(yǎng)量較少，每個(gè)場平均飼養(yǎng)規(guī)模也小。不同品種3個(gè)聯(lián)合遺傳評(píng)估組在總產(chǎn)仔數(shù)上均差異顯著(P<0.05)，組1母系品種(長白豬和大白豬)的繁殖性能均優(yōu)于其他聯(lián)合遺傳評(píng)估組。

2.3遺傳參數(shù)

表5列出了3個(gè)品種總產(chǎn)仔數(shù)遺傳參數(shù)估計(jì)結(jié)果。大白豬TNB性狀在組1、2和3中的遺傳力估計(jì)值分別為0.11、0.06和0.10；長白豬分別為0.10、0.08和0.07；杜洛克豬分別為0.07、0.08和0.06。在不同聯(lián)合遺傳評(píng)估組及不同品種之間TNB遺傳力估計(jì)值趨于一致，均在0.10左右。表5也可以看出，群體規(guī)模對遺傳參數(shù)及各方差組分估計(jì)準(zhǔn)確性影響很大，總體上看，組1中各品種豬的TNB遺傳力及各方差組分估計(jì)值標(biāo)準(zhǔn)誤較其他兩組要小，尤其是大白豬組1中加性遺傳方差，永久環(huán)境效應(yīng)方差，殘差方差及TNB遺傳力估計(jì)的標(biāo)準(zhǔn)誤分別為0.065、0.053、0.045及0.006，這可能與該組群體規(guī)模最大有關(guān)，相關(guān)結(jié)果顯示(表6)，有效記錄數(shù)與遺傳力估計(jì)值的標(biāo)準(zhǔn)誤成強(qiáng)負(fù)相關(guān)(-0.77)，預(yù)示著群體規(guī)模越大，遺傳參數(shù)及各方差組分估計(jì)越準(zhǔn)確。綜觀3個(gè)品種，杜洛克由于其3個(gè)關(guān)聯(lián)組群體規(guī)模小，遺傳參數(shù)及各方差組分估計(jì)標(biāo)準(zhǔn)誤要比大白、長白大。同樣，受群體規(guī)模影響，3個(gè)品種組1估計(jì)的加性遺傳方差也高于組2和組3，其中大白豬組1加性遺傳方差最高。3個(gè)品種中，杜洛克加性遺傳方差最小，這可能與其每個(gè)組數(shù)據(jù)量小有關(guān)。永久環(huán)境效應(yīng)占表型變異的比例在11%以下，不同聯(lián)合遺傳評(píng)估組及不同品種之間差異較大。

表1　我國主要豬品種的聯(lián)合遺傳評(píng)估組分組情況

聯(lián)合遺傳評(píng)估組內(nèi)為2個(gè)育種場的未予考慮

Groups with 2 pig farms were not took into account in this study

表2　大白、長白及杜洛克豬平均關(guān)聯(lián)率數(shù)據(jù)統(tǒng)計(jì)表

表3　大白、長白及杜洛克豬聯(lián)合遺傳評(píng)估組及全國關(guān)聯(lián)率統(tǒng)計(jì)表

表4　總產(chǎn)仔數(shù)在不同豬種及聯(lián)合遺傳評(píng)估組內(nèi)表型值

同一行中不同字母表示差異顯著(P<0.05)

Different letters in the same row means significant difference between different groups(P<0.05)

表5　總產(chǎn)仔數(shù)在不同豬種及聯(lián)合遺傳評(píng)估組內(nèi)方差組分和遺傳參數(shù)估計(jì)值

表6　記錄數(shù)與遺傳力估計(jì)值及標(biāo)準(zhǔn)誤之間的相關(guān)

3討論

3.1聯(lián)合遺傳評(píng)估組

聯(lián)合育種是目前國內(nèi)外豬育種的共識(shí)，無論是丹麥的國家育種體系，還是PIC等國際育種公司的育種策略，都是將多個(gè)場聯(lián)合起來擴(kuò)大育種規(guī)模，可以提高選擇強(qiáng)度、擴(kuò)大遺傳變異、提高育種值估計(jì)準(zhǔn)確性，進(jìn)而加快遺傳進(jìn)展。本研究結(jié)果表明，3個(gè)品種中組1聯(lián)合的場最多，規(guī)模也最大，因而其遺傳變異(加性遺傳方差)最大，遺傳參數(shù)估計(jì)的準(zhǔn)確性也最高，這將有助于遺傳進(jìn)展的加快。我國《生豬遺傳改良計(jì)劃》也將聯(lián)合育種作為生豬遺傳改良的核心工作，并且希望通過篩選國家生豬核心育種場，推動(dòng)聯(lián)合育種工作。場間遺傳聯(lián)系是聯(lián)合育種的前提，C.Y.Sun等[15]分析了全國27家豬場的長白豬數(shù)據(jù)、36家豬場的大白豬數(shù)據(jù)及18家杜洛克豬數(shù)據(jù)，有彼此場間聯(lián)系的豬場數(shù)目很少，按本研究的方法劃分關(guān)聯(lián)組，大白豬共得到4個(gè)聯(lián)合遺傳評(píng)估組，除去組1中有8個(gè)種豬場外，其余3個(gè)組中每個(gè)組僅有2個(gè)場，長白豬得到2個(gè)組，杜洛克豬得到1個(gè)組，且每個(gè)組僅有2個(gè)豬場。W.Xiao等針對北京地區(qū)種豬場間關(guān)聯(lián)率研究表明，目前區(qū)域性的種豬聯(lián)合遺傳評(píng)估條件也不滿足[16]。

雖然在生豬遺傳改良計(jì)劃的推動(dòng)下，我國豬育種工作取得了一定進(jìn)展，場間遺傳聯(lián)系得到了部分加強(qiáng)，但本研究表明，我國目前整體場間關(guān)聯(lián)情況仍然較差，還遠(yuǎn)達(dá)不到全國統(tǒng)一聯(lián)合評(píng)估標(biāo)準(zhǔn)，因此目前全國種豬遺傳評(píng)估中心仍然將各場作為一個(gè)獨(dú)立的單位進(jìn)行遺傳評(píng)估(參見全國種豬遺傳評(píng)估報(bào)告www.cnsge.org.cn)，長期下去，將使聯(lián)合育種成為空中樓閣。而將彼此有遺傳聯(lián)系的場劃分為聯(lián)合遺傳評(píng)估組，逐步進(jìn)行適度規(guī)模的聯(lián)合評(píng)估不失為一個(gè)切實(shí)可行的措施。本研究中的大白豬組1規(guī)模已經(jīng)達(dá)到了61 161條，其生長性狀則達(dá)到了106 858條，基本上可以保證準(zhǔn)確地進(jìn)行遺傳參數(shù)和育種值估計(jì)，育種效果將大大高于單個(gè)場。同時(shí)，通過人工受精技術(shù)或種豬銷售，實(shí)現(xiàn)優(yōu)良種公豬的跨場使用，進(jìn)一步建立和加強(qiáng)場間遺傳聯(lián)系[17-18]，最終逐步實(shí)現(xiàn)全國性的聯(lián)合育種。

從表2可以看出，3個(gè)品種豬聯(lián)合遺傳評(píng)估組內(nèi)各場的平均關(guān)聯(lián)率均大于3%，XJMY與HNMY之間的場間關(guān)聯(lián)率在大白豬和長白豬中表現(xiàn)突出，分別達(dá)到了28.2%和23.5%，分別以兩場有測定數(shù)據(jù)的個(gè)體為基礎(chǔ)，向上追蹤4代系譜，進(jìn)行對比后發(fā)現(xiàn)，兩場擁有的共同公豬數(shù)分別達(dá)到了135和46頭。SJYZ與GGSP之所以存在較大的關(guān)聯(lián)性是因?yàn)镾JYZ于2011年底從GGSP引進(jìn)大白豬100多頭，有直接的遺傳聯(lián)系，通過追蹤系譜發(fā)現(xiàn)，兩個(gè)場擁有的共同公豬數(shù)也達(dá)到了79頭。本研究發(fā)現(xiàn)，聯(lián)合遺傳評(píng)估組呈現(xiàn)有區(qū)域性的特點(diǎn)，本省或者鄰近省份之間的場間關(guān)聯(lián)情況較好，尤為突出的是兩廣地區(qū)(廣東和廣西)，這與實(shí)際情況是相符的，兩廣地區(qū)由于地緣、氣候等原因自古聯(lián)系緊密，又同為我國養(yǎng)豬發(fā)達(dá)省份，再考慮到運(yùn)輸成本及疾病等原因，不難理解本省或者鄰近省份場間聯(lián)系情況較好，表1中即發(fā)現(xiàn)大白、長白及杜洛克組1中兩廣地區(qū)的豬場數(shù)分別達(dá)到了5、6、5家。

3.2遺傳方差

要準(zhǔn)確的解釋文獻(xiàn)中遺傳參數(shù)估計(jì)結(jié)果的不同是很困難的，因?yàn)檫z傳參數(shù)估計(jì)的準(zhǔn)確與否取決于群體結(jié)構(gòu)(包括群體大小、分布、生產(chǎn)管理水平和方式等)、記錄系統(tǒng)(生產(chǎn)記錄的系統(tǒng)性、全面性和準(zhǔn)確性)和估計(jì)方法3個(gè)因素[4]。J.Estany等[19]基于丹麥大白和長白豬群體，利用與本研究相同的模型，其群體規(guī)模與本研究組1接近(大白29 336條；長白19 666條)，加性遺傳方差估計(jì)值與本研究組1中結(jié)果一致，均高于組2、3中的估計(jì)值。M.Tomiyama等[20]利用日本大白豬群體亦得到類似結(jié)果，其大白豬群體規(guī)模達(dá)到了近一萬條，因此，可以推測組1中的加性遺傳方差較其他組大的原因可能是由于其群體規(guī)模所致，尤其是組1中大白豬估計(jì)值達(dá)到了0.910，明顯高于其他組。

3.3遺傳力

3個(gè)聯(lián)合遺傳評(píng)估組中不同品種的遺傳力估計(jì)值均在0.10左右，這與相關(guān)文獻(xiàn)報(bào)道的結(jié)果是一致的[7，21-23]。繁殖性狀一般是低遺傳力，有研究基于泰國大白豬群體得到的前4個(gè)胎次TNB遺傳力估計(jì)值范圍為0.02～0.11[24]，E.Akanno等對熱帶豬群體繁殖性狀(NBA)進(jìn)行了遺傳參數(shù)估計(jì)的整合分析(Meta-analysis)，該研究分析了29篇文獻(xiàn)得到的產(chǎn)活仔數(shù)(NBA)遺傳力估計(jì)值為0.08[25]。M.Tomiyama等基于日本豬群體，采用和本研究相同的模型和方法，得到的TNB遺傳力估計(jì)值為0.09[20]，本研究結(jié)果同樣與美國全國豬改良聯(lián)合會(huì)(NSIF)所推薦的0.10的遺傳力估計(jì)值是一致的[26]。

3.4永久環(huán)境效應(yīng)

動(dòng)物個(gè)體的某些性狀可以有多個(gè)觀察值，例如母豬的產(chǎn)仔數(shù)和奶牛的產(chǎn)奶量。對于產(chǎn)仔數(shù)記錄，在估計(jì)永久環(huán)境效應(yīng)時(shí)，有些學(xué)者認(rèn)為應(yīng)當(dāng)把不同胎次的記錄看做是不同的性狀[27-28]，而有些學(xué)者則認(rèn)為應(yīng)該采用重復(fù)力模型進(jìn)行估計(jì)[21-23]。當(dāng)前，美國的種豬遺傳評(píng)估和遺傳系統(tǒng)(STAGES)、加拿大豬遺傳改良中心(CCSI)及澳大利亞的國家種豬提高計(jì)劃(NPIP)等都把不同胎次的記錄看做是一個(gè)性狀的重復(fù)測定值。本研究對總產(chǎn)仔數(shù)分析時(shí)采用了重復(fù)力模型，并對影響其所有胎次的永久環(huán)境效應(yīng)進(jìn)行了估計(jì)。本研究得到的結(jié)果與相關(guān)文獻(xiàn)報(bào)道結(jié)果基本一致[7，21-22]，但J.B.Ferraz等[29]采用與本研究相同模型卻得到16%的永久環(huán)境效應(yīng)估計(jì)值。

參考文獻(xiàn)(References)：

[1]WEBB A J.Objectives and strategies in pig improvement：An applied perspective[J].JDairySci， 1998，81(Suppl 2):36-46.

[3]GR?N PEDERSEN O.Danbred，Danish pig breeding in 100 years-an overview[C]//Danish pig breeding，100 years Anniversary Seminar.Herning.1997，23.

[4]張勤.家畜育種值和遺傳參數(shù)估計(jì)方法的發(fā)展及現(xiàn)狀[J].草食家畜，1990(6):1-4，46.

ZHANG Q.The development and status of methods of estimated breeding value and genetic parameters[J].Grass-FeedingLivestock，1990(6):1-4，46.(in Chinese)

[5]DARFOUR-ODURO K A，NAAZIE A，AHUNU B K，et al.Genetic parameter estimates of growth traits of indigenous pigs in Northern Ghana[J].LivestSci，2009，125(2):187-191.

[6]MATHUR P，SULLIVAN B，CHESNAIS J.Measuring connectedness：concept and application to a large industry breeding program.[C]// Proceedings of the 7th World Congress on Genetics Applied to Livestock Production：2002:23.

[7]張鎖宇，王愛國，王勝，等.豬主要繁殖性狀的遺傳參數(shù)估計(jì)[J].中國畜牧雜志，2013，49(23):1-5.

ZHANG S Y，WANG A G，WANG S，et al.Estimation of genetic parameters on reproductive traits in pigs[J].ChineseJournalofAnimalScience，2013，49(23):1-5.(in Chinese)

[8]王青來，王愛國.應(yīng)用多性狀動(dòng)物模型 BLUP 和 REML 對杜洛克豬的遺傳分析[D].北京:中國農(nóng)業(yè)大學(xué)，2000.

WANG Q L，WANG A G.Genetic analysis in Duroc pigs using multiple-trait Model BLUP and REML[D].Beijing：China Agricultural University，2000.(in Chinese)

[9]李玉華.應(yīng)用動(dòng)物模型 BLUP 和 REML 對大白豬遺傳評(píng)定和遺傳參數(shù)估計(jì)[D].泰安:山東農(nóng)業(yè)大學(xué)，2003.

LI Y H.Genetic evaluation and genetic parameter estimation in largewhite pigs using Animal Model BLUP and REML[D].Tai’an：Shandong Agricultural University，2003.(in Chinese)

[10]MISZTAL I.Restricted maximum likelihood estimation of variance components in animal model using sparse matrix inversion and a supercomputer[J].JDairySci， 1990，73(1):163-172.

[11]張沅.動(dòng)物育種學(xué)[M].北京：中國農(nóng)業(yè)出版社，2001.

ZHANG Y.Animal breeding[M].Beijing：China Agriculture Press，2001.(in Chinese)

[12]KLEI B，TSURUTA S.Approximate variance for heritability estimates.http://nce.ads.uga.edu/html/projects/AI_SE.pdf.

[13]MADSEN P，SORENSEN P，SU G，et al.DMU-a package for analyzing multivariate mixed models.[C]// World Congress on Genetics Applied to Livestock Production，2006.

[14]國家生豬核心育種場遺傳評(píng)估報(bào)告[R].北京：全國畜牧總站，2014.http://www.cnsge.org.cn/wzgg/detail.jsp?ArticleID=2006022555.

The National Swine Genetic Evaluation Report[R].Beijing：The national animal husbandry station，2014.(in Chinese) http://www.cnsge.org.cn/wzgg/detail.jsp?ArticleID=2006022555.

[15]SUN C Y，WANG C K，WANG Y C，et al.Evaluation of connectedness between herds for three pig breeds in China[J].Animal，2009，3(4):482-485.

[16]XIAO W，DING X D，WANG C D，et al.The connectedness on large white and landrace in regional joint breeding system in Beijing[J].JAnimVetAdv，2010，9(18):2338-2342.

[17]王愛國.建立場間遺傳聯(lián)系實(shí)現(xiàn)豬的聯(lián)合育種[J].養(yǎng)殖與飼料，2003(10):10-13.

WANG A G.Establish the connectedness between herds to achieve swine joint breeding[J].AnimalsBreedingandFeed，2003(10):10-13.(in Chinese)

[18]王愛國，李力，梅克義，等.加強(qiáng)聯(lián)合育種 提高種豬質(zhì)量[J].中國畜牧雜志，2004，40(3):5-7.

WANG A G，LI L，MEI K Y，et al.Enhance joint breeding and improve the quality of Swine.ChineseJournalofAnimalScience，2004，40(3):5-7.(in Chinese)

[19]ESTANY J，SORENSEN D.Estimation of genetic parameters for litter size in Danish Landrace and Yorkshire pigs[J].AnimSci，1995，60(2):315-324.

[20]TOMIYAMA M，KUBO S，TAKAGI T，et al.Evaluation of genetic trends and determination of the optimal number of cumulative records of parity required in reproductive traits in a Large White pig population[J].AnimSciJ，2011，82(5):621-626.

[21]CHEN P，BASS T J，MABRY J W，et al.Genetic parameters and trends for litter traits in U.S.Yorkshire，Duroc，Hampshire，and Landrace pigs[J].JAnimSci，2003，81(1):46-53.

[22]CRUMP R E，THOMPSON R，HALEY C S，et al.Individual animal model estimates of genetic correlations between performance test and reproduction traits of landrace pigs performance tested in a commercial nucleus herd[J].AnimSci，1997，65(2):291-298.

[23]SEE M T，MABRY J W，BERTRAND J K.Restricted maximum likelihood estimation of variance components from field data for number of pigs born alive[J].JAnimSci，1993，71(11):2905-2909.

[24]IMBOONTA N，RYDHMER L，TUMWASORN S.Genetic parameters and trends for production and reproduction traits in Thai Landrace sows[J].LivestSci，2007，111(1):70-79.

[25]AKANNO E，SCHENKEL F S，QUINTON M，et al.Meta-analysis of genetic parameter estimates for reproduction，growth and carcass traits of pigs in the tropics[J].LivestSci，2013，152(2-3):101-113.

[26]HUBBARD D.Guidelines for uniform swine improvement programs[R].USA：Dept.of Agriculture，1981.

[27]ROEHE R，KENNEDY B W.Estimation of genetic parameters for litter size in Canadian Yorkshire and Landrace swine with each parity of farrowing treated as a different trait[J].JAnimSci，1995，73(10):2959-2970.

[29]FERRAZ J B，JOHNSON R K.Animal model estimation of genetic parameters and response to selection for litter size and weight，growth，and backfat in closed seedstock populations of large white and landrace swine[J].JAnimSci，1993，71(4):850-858.

(編輯郭云雁)

Estimation of Genetic Parameters by Connecting Groups for Reproductive Traits of Yorkshire，Landrace and Duroc Pigs in China

ZHANG Suo-yu1，2#，QIU Xiao-tian2，3#，DING Xiang-dong1，ZHANG Qin1*，WANG Zhi-gang2，WANG Chang-cun2，QU Yun-long1，2

(1.KeyLaboratoryofAnimalGeneticsandBreedingofMinistryofAgriculture，NationalEngineeringLaboratoryforAnimalBreeding，CollegeofAnimalScienceandTechnology，ChinaAgriculturalUniversity，Beijing100193，China;2.NationalAnimalHusbandryService，Beijing100125，China;3.ChongqingAcademyofAnimalSciences，Chongqing402460，China)

Abstract:The objective of this study was to estimate genetic parameters for total number born(TNB) of Yorkshire，Landrace and Duroc pigs in different genetic evaluation groups，which could establish foundation for the national swine joint genetic evaluation.Estimation of the degree of connectedness among 74 national nucleus farms were calculated by the module DMU4 in DMU software，based on the principle of each herd in group with a minimum average connectedness rating of 3% and all herds in group connected with each other，3 groups for each breed were screened(i.e.group1，group2 and group3).Animal model and DMU software were used to estimate variances of animal genetic and permanent environmental for total number born(TNB) of Yorkshire，Landrace and Duroc pigs in 3 groups，fixed effects of contemporary groups were also included in the analysis，records on 61 161，19 974 and 18 002 for Yorkshire，28 239，11 918 and 3 064 for Landrace，and 7 174，6 512 and 2 417 for Duroc in the group1，2 and 3，respectively，were analyzed.As a result，in general，the degree of connectedness was low(the average connectedness rating was 0.20%) and large-scale across-herd genetic evaluation in China was not practicable at present，but it might carry out across-herd genetic evaluation in certain regions.Estimated values of permanent environmental effects ranged from 0.05 to 0.11，estimated values of heritabilities for TNB were approximately 0.10 for all breeds and groups，the results were within the range recommended by relevant literature.The results indicate that genetic parameters of group 1 were better than other groups，therefore，it would expect to be used in the national swine joint genetic evaluation in the future.

Key words:pigs；genetic parameters；total number born；connectedness rating

doi:10.11843/j.issn.0366-6964.2016.03.003

收稿日期：2015-06-25

基金項(xiàng)目：國家生豬產(chǎn)業(yè)技術(shù)體系(CARS-36)

作者簡介：張鎖宇(1988-)，男，黑龍江人，碩士，主要從事種豬遺傳評(píng)估方面研究，Tel:010-62893046，E-mail:suoyuzhang@foxmail.com；邱小田，男，博士，主要從事動(dòng)物遺傳育種方面研究，Tel:010-59194622，E-mail:1967699090@qq.com。二者并列為第一作者 *通信作者：張勤，教授，Tel：010-62732634，E-mail:qzhang@cau.edu.cn

中圖分類號(hào)：S828；S813.3

文獻(xiàn)標(biāo)志碼：A

文章編號(hào):0366-6964(2016)03-0429-10