MUC13、FUT1基因在2個(gè)種豬核心群中的分子標(biāo)記輔助選擇研究

楊　明，王青來，劉敬順，劉珍云，溫淑賢，吳珍芳，3，羅旭芳，3(廣東溫氏食品集團(tuán)股份有限公司，廣東新興527439;2國家生豬種業(yè)工程技術(shù)研究中心，廣東新興527439; 3華南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院，廣東廣州50642)

MUC13、FUT1基因在2個(gè)種豬核心群中的分子標(biāo)記輔助選擇研究

楊明1，2，王青來1，2，劉敬順1，2，劉珍云1，2，溫淑賢1，吳珍芳1，2，3，羅旭芳1，3
(1廣東溫氏食品集團(tuán)股份有限公司，廣東新興527439;2國家生豬種業(yè)工程技術(shù)研究中心，廣東新興527439; 3華南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院，廣東廣州510642)

摘要:【目的】MUC13、FUT1基因分別為大腸埃希菌引起的斷奶前和斷奶后仔豬腹瀉的主效基因，通過探索2個(gè)基因在不同種豬群中的分子標(biāo)記輔助選擇方法，提高種豬抗病性的選育效率.【方法】利用PCR-SNaPshot和PCRRFLP的方法分別檢測影響斷奶前仔豬腹瀉的MUC13基因及斷奶后仔豬腹瀉的FUT1基因在溫氏集團(tuán)2個(gè)專門化父系群體(666頭杜洛克和512頭皮特蘭)中的基因分布情況.【結(jié)果和結(jié)論】MUC13基因頻率在2個(gè)群體中表現(xiàn)出明顯的群體特異性，有利等位基因頻率分別為0. 890和0. 180，而FUT1基因則相差不大，分別為0. 754和0. 677.MUC13與FUT1基因的有利等位基因型組合個(gè)體比例在2個(gè)群體中差異較大，在杜洛克群體中達(dá)36. 75%，而在皮特蘭群體中僅為3. 49%.通過基因型與選育性狀表型的相關(guān)分析，發(fā)現(xiàn)優(yōu)勢基因型對皮特蘭群體的體型外貌負(fù)面影響較大，但卻不影響杜洛克群體.綜合現(xiàn)場育種實(shí)踐，建議杜洛克群體先實(shí)現(xiàn)MUC13基因的輔助選育純化，再啟動(dòng)FUT1基因的選育工作.皮特蘭群體有利等位基因型組合個(gè)體比例太低，暫時(shí)不宜進(jìn)行抗腹瀉基因純化選育.

關(guān)鍵詞:MUC13基因; FUT1基因;種豬核心群;分子標(biāo)記輔助選擇

楊明，王青來，劉敬順，等.MUC13、FUT1基因在2個(gè)種豬核心群中的分子標(biāo)記輔助選擇研究［J］.華南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2015，36(6):1-8.

優(yōu)先出版時(shí)間:2015-10-16

優(yōu)先出版網(wǎng)址: http: / /www.cnki.net/kcms/detail/44.1110.s.20151016.1446.008.html

仔豬腹瀉是養(yǎng)豬生產(chǎn)中的常見傳染病，尤其是這幾年，我國頻頻暴發(fā)，仔豬腹瀉導(dǎo)致的死亡率達(dá)11. 5%～29. 5%，其中由大腸埃希菌Escherichia coli引起的發(fā)病率和死亡率分別占腹瀉發(fā)病率和死亡率的56. 2%和24. 7%［1］，給養(yǎng)豬業(yè)造成了巨大經(jīng)濟(jì)損失.目前，針對仔豬腹瀉，雖然可以通過改善飼養(yǎng)環(huán)境、緩慢斷奶、添加維生素和外源蛋白、使用藥物、隔離撲殺、免疫接種等措施防治，但由于新傳染病的發(fā)生和流行，病毒發(fā)生變異快，超強(qiáng)毒株不斷出現(xiàn)，疫苗效力不夠，并有藥物殘留等問題，使仔豬腹瀉病無法得到有效防治［2-5］.此外，用疫苗和藥物治療的方法會(huì)導(dǎo)致豬抵抗力下降，耐藥性增強(qiáng)，豬群健康度無法提高.因此，通過分子育種的方法培育抗病品種，增強(qiáng)豬只抗病力，從而解決由于片面追求高產(chǎn)而導(dǎo)致的高產(chǎn)與健康等性狀間的遺傳拮抗［6］.遺傳學(xué)家和育種學(xué)家們通過研究鑒定了由大腸埃希菌引起的仔豬腹瀉的主效基因黏蛋白13(Mucin 13，MUC13)和α(1，2)巖藻糖轉(zhuǎn)移酶(Alpha (1，2 ) fucosyltransferase，F(xiàn)UT1)的主效位點(diǎn)［7-10］.這2個(gè)基因分別影響斷奶前仔豬腹瀉和斷奶后仔豬腹瀉，且2003年在法國和西班牙利用FUT1基因分子輔助選擇育種首先推出了能特異抵抗F18ab仔豬腹瀉的新品系.

溫氏集團(tuán)2005年，在法系皮特蘭群體中針對斷奶后仔豬腹瀉進(jìn)行過豬FUT1基因的檢測，發(fā)現(xiàn)豬只對腹瀉的抵抗力有差異.2009年針對杜洛克群體斷奶前仔豬腹瀉進(jìn)行過MUC13基因位點(diǎn)的檢測，發(fā)現(xiàn)其抗腹瀉等位基因頻率達(dá)80%左右.本研究利用群體基因頻率分析、基因型與性狀相關(guān)分析等方法，通過估算這2個(gè)抗腹瀉基因在溫氏杜洛克專門化父系S22和皮特蘭專門化父系S11核心群中的差異表現(xiàn)，制定特異的抗腹瀉集成育種方案，進(jìn)行不同的集成分子育種篩選，輔助選留，為培育抗病品系奠定基礎(chǔ).

1　材料與方法

1.1試驗(yàn)動(dòng)物

試驗(yàn)動(dòng)物來自溫氏集團(tuán)國家生豬種業(yè)工程技術(shù)研究中心研發(fā)基地沙湖原種場杜洛克核心群體和清遠(yuǎn)一區(qū)原種場皮特蘭核心群體.樣本總數(shù)1 178頭，其中核心群杜洛克666頭，核心群皮特蘭512頭.利用耳號(hào)鉗剪取已經(jīng)利用乙醇消毒過耳尖的待檢豬只耳樣2～3 g，放入裝有1. 5 mL體積分?jǐn)?shù)為75%乙醇的離心管內(nèi)，低溫保存.

1.2基因型判型

采取的耳組織樣本送至江西農(nóng)業(yè)大學(xué)動(dòng)物生物技術(shù)國家重點(diǎn)實(shí)驗(yàn)室培育基地檢測2個(gè)抗腹瀉基因，檢測方法如下:采用酚－三氯甲烷法手工提取基因組DNA.2個(gè)基因位點(diǎn)的引物序列及具體擴(kuò)增情況見表1，PCR反應(yīng)體系總體積為20 μL，包括20 ng·μL－1模板DNA 2 μL，10×buffer 2. 0 μL，25 mmol·L－1MgCl21. 2 μL，10 mmol·L－1dNTP 0. 3 μL，PF 0. 4 μL，PR 0. 4 μL，Taq聚合酶0. 5 μL，超純水13. 2 μL.PCR循環(huán)參數(shù): 94℃變性3 min，退火溫度見表1，72℃延伸10 min，共36個(gè)循環(huán).

表1　引物與擴(kuò)增情況Tab.1 PCR primers and the amplicons

1.2.1 PCR-SNaPshot法PCR-SNaPshot的基本原理是通過引物在標(biāo)有不同熒光的ddNTP(A、G、C、T)反應(yīng)體系內(nèi)進(jìn)行一步延伸后，使用3130XL型DNA測序儀檢測不同的熒光信號(hào).MUC13的基因型利用PCR-SNaPshot法進(jìn)行檢測，其SNaPshot反應(yīng)體系為5 μL，其中含2 μL SNaPshot multiplexmix (含Taq聚合酶和熒光標(biāo)記的ddNTPs)，1. 5 μL純化后PCR產(chǎn)物，1. 2 μL去離子水，0. 3 μL SNaPshot引物.SNaP-shot反應(yīng)體系的擴(kuò)增程序?yàn)? 96℃10 s，50℃5 s，60℃30 s，25個(gè)循環(huán).反應(yīng)后，在SNaPshot反應(yīng)產(chǎn)物中加0. 57 μL 1×NEB buffer和0. 1 μL的CIP，37℃反應(yīng)60 min以清除熒光標(biāo)記的ddNTPs，隨后75℃反應(yīng)15 min滅活純化酶;最后，每1 μL SNaPshot反應(yīng)產(chǎn)物加8 μL Hi-Di formamide與GeneScan 120 LIZ size standard混合物(兩者體積比為20∶1)，經(jīng)變性后上樣于ABI3130XL自動(dòng)遺傳分析儀(ABI，USA)進(jìn)行電泳，利用GeneMapper Software version 4. 0(ABI，USA)進(jìn)行基因型判定和數(shù)據(jù)收集.

1.2.2 PCR-RFLP法FUT1主效位點(diǎn)是通過Hin6I PCR-RFLP判定基因型，檢測酶切反應(yīng)體系與條件如下:反應(yīng)體系總體積為15 μL，其中10×Buffer 1. 5 μL，HinP1I內(nèi)切酶0. 3 μL，DDW 8. 2 μL，PCR產(chǎn)物5 μL.37℃水浴反應(yīng)，過夜.PCR產(chǎn)物通過0. 02 g·mL－1的瓊脂糖凝膠電泳來判型: A等位基因可酶切為3 283和93 bp 2條帶，G等位基因表現(xiàn)為241、93、87 bp 3條帶.

1.3統(tǒng)計(jì)分析

通過R軟件進(jìn)行表型與基因型相關(guān)分析，模型為: y = 1μ+ Zα+ e.其中，y為所測得的表型值，μ為總體平均值，Z為基因型的指示矩陣，α為SNP隨機(jī)加性遺傳效應(yīng)向量，且α—N(0，Gσ2α) (G為分子血緣相關(guān)矩陣，σ2α為加性遺傳方差)，Zα為基因型效應(yīng)，e為殘差.

2　結(jié)果與分析

2.1基因頻率和基因型頻率

2.1.1 MUC13的基因型及基因頻率從2個(gè)父系群體的MUC13的基因型頻率及基因頻率分布(表2)來看，2個(gè)群體表現(xiàn)差異較大，G等位基因?yàn)榭剐杂欣?，杜洛克群體有利等位基因頻率高達(dá)0. 890，而皮特蘭群體僅為0. 180.

表2　MUC13基因型及基因頻率的分布1)Tab.2 MUC13 genotypes and allele frequencies in tested breeds

2.1.2 FUT1基因型及基因頻率從2個(gè)父系群體的FUT1基因型頻率及基因頻率分布(表3)來看，2個(gè)群體在此位點(diǎn)表現(xiàn)差異不大，A等位基因?yàn)榭剐杂欣任换?，有利等位基因頻率偏高;杜洛克豬未檢測到GG基因型個(gè)體，皮特蘭豬GG頻率也相對較低.

表3　FUT1基因型及基因頻率的分布1)Tab.3 FUT1 genotypes and allele frequencies in tested breeds

2.2 2個(gè)基因與2個(gè)群體的生長性狀相關(guān)分析

2.2.1 MUC13的基因型與杜洛克選育表型的相關(guān)分析由MUC13與杜洛克群體18個(gè)選育性能的相關(guān)分析結(jié)果(表4)可以看出，在父系指數(shù)、母系指數(shù)以及背膘育種估計(jì)值(EBV)選育性能中，GG型個(gè)體要明顯優(yōu)于AA和AG型個(gè)體;而在其余的性能中，則3種基因型表現(xiàn)差異不大.

2.2.2 FUT1的基因型與杜洛克選育表型的相關(guān)分析由FUT1與杜洛克群體18個(gè)選育性能的相關(guān)分析結(jié)果(表5)可以看出，在父系指數(shù)、母系指數(shù)、日增質(zhì)量EBV以及背膘EBV等選育性能中，GA型個(gè)體要明顯優(yōu)于AA型個(gè)體;而在出生體質(zhì)量和100 kg體質(zhì)量日齡性能上，AA型個(gè)體要明顯優(yōu)于GA型個(gè)體; 2種基因型在其余性能中表現(xiàn)差異不大.

表4　MUC13的基因型與杜洛克生長性能的關(guān)系1)Tab.4 The relationship between MUC13 genotype and growth performance in Duroc

表5　FUT1的基因型與杜洛克生長性能的關(guān)系1)Tab.5 The relationship between FUT1 genotype and growth performance in Duroc

1) AA型樣本數(shù)量為314頭，GA型樣本數(shù)量為315頭.

2.2.3 MUC13基因型與皮特蘭選育表型的相關(guān)分析由MUC13的基因型與皮特蘭群體18個(gè)選育性能的相關(guān)分析結(jié)果(表6)可以看出，在體高、頭型評分、肢蹄評分等體貌表型上，GG型個(gè)體要明顯差于AA和AG型個(gè)體;而GG型個(gè)體的校正100 kg瘦肉率明顯高于AA和AG型個(gè)體; AG雜合子在校正100 kg體質(zhì)量日齡方面要明顯優(yōu)于AA和GG型個(gè)體; 3種基因型在其他性能中表現(xiàn)差異不大.

表6　MUC13基因型與皮特蘭生長性能的關(guān)系1)Tab.6 The relationship between MUC13 genotype and growth performance in Pietrain

2.2.4 FUT1的基因型與皮特蘭選育表型的相關(guān)分析由FUT1的基因型與皮特蘭群體18個(gè)選育性能的相關(guān)分析結(jié)果(表7)可以看出，在體高、頭型評分體貌表型上，GG型個(gè)體要明顯優(yōu)于AA和GA型個(gè)體;背膘EBV和出生體質(zhì)量等選育的生長性能方面，AA型個(gè)體要明顯優(yōu)于其他2種基因型個(gè)體;而雜合子GA型個(gè)體在校正100 kg瘦肉率方面表現(xiàn)較好;在其他性能中基因型間表現(xiàn)差異不大.

表7　FUT1基因型與皮特蘭生長性能的關(guān)系1)Tab.7 The relationship between FUT1 genotype and growth performance in Pietrain

2.3 2個(gè)抗腹瀉基因組合基因型分布情況

圖1表明:杜洛克群體有5種組合類型，其中MUC13與FUT1基因均為優(yōu)勢基因型組合GG/AA的個(gè)體數(shù)為222頭，占群體總量的36. 75% (222/604).圖2結(jié)果表明:皮特蘭群體有8種組合類型，僅有11頭是優(yōu)勢基因型組合個(gè)體，占群體總量的3. 49% (11/315).

2.4 2個(gè)抗腹瀉基因組合基因型與2個(gè)群體生長性狀相關(guān)分析

2.4.1 MUC13、FUT1組合基因型與杜洛克群體18個(gè)選育性能的相關(guān)分析MUC13、FUT1組合基因型與杜洛克群體18個(gè)選育性能的相關(guān)分析結(jié)果(表8)可以看出，在父系指數(shù)、母系指數(shù)、背膘EBV以及日增質(zhì)量EBV選育性能中，GG/GA型個(gè)體要明顯優(yōu)于其他組合型個(gè)體;在校正100 kg體質(zhì)量日齡、出生體質(zhì)量方面，AA/AA組合型要比其他組合型占優(yōu)勢.其余性狀中，幾個(gè)組合型表現(xiàn)無差異.2.4.2 MUC13、FUT1組合基因型與皮特蘭群體18個(gè)選育性能的相關(guān)分析MUC13、FUT1組合基因型與皮特蘭群體18個(gè)選育性能的相關(guān)分析結(jié)果見表9.由于僅有1頭AG/GG組合型個(gè)體，數(shù)據(jù)準(zhǔn)確性不夠，不做比較.由表9可以看出，在體高、頭型評分、肢蹄評分以及出生體質(zhì)量表型上，AA/AA組合型要比其他組合型占優(yōu)勢;在校正100 kg瘦肉率、校正100 kg背膘厚以及背膘EBV性能中，GG/AA組合型個(gè)體要明顯優(yōu)于其他組合型個(gè)體;在其余性能中各組合型間表現(xiàn)差異不大.

圖1　杜洛克群體MUC13與FUT1基因型組合分布Fig.1　 Distributions of MUC13 and FUT1 genotypic combinations in Duroc pigs

圖2　皮特蘭群體MUC13與FUT1基因型組合分布Fig.2　 Distributions of MUC13 and FUT1 genotypic combinations in Pietrain pigs

表8　MUC13、FUT1組合基因型與杜洛克生長性能的關(guān)系Tab.8 The relationship between MUC13，F(xiàn)UT1 genotypic combination and growth performance in Duroc

表9　MUC13、FUT1組合基因型與皮特蘭生長性能的關(guān)系Tab.9 The relationship between MUC13，F(xiàn)UT1 genotypic combination and growth performance in Pietrain

3　討論與結(jié)論

3.1溫氏杜洛克核心群體的抗腹瀉基因選育方案

杜洛克作為商業(yè)豬種中的重要父系品種，其與生產(chǎn)相關(guān)的各性狀都經(jīng)過了嚴(yán)格的選育，尤其是核心育種群.本研究在溫氏杜洛克核心育種群中檢測抗大腸埃希菌引起的斷奶前仔豬腹瀉基因MUC13以及抗大腸埃希菌引起的斷奶后仔豬腹瀉基因FUT1的主效位點(diǎn)變異情況，結(jié)果顯示:在溫氏集團(tuán)杜洛克核心育種群中，MUC13優(yōu)勢等位基因頻率高達(dá)89. 0%，這與阮國榮等［11］在福建省杜洛克群體中的檢測結(jié)果一致，也與溫氏集團(tuán)生產(chǎn)實(shí)踐中杜洛克種豬群體斷奶前仔豬腹瀉比例低的情況一致，只需對少量的劣勢等位基因純合個(gè)體及公豬雜合個(gè)體進(jìn)行淘汰，對雜合母豬根據(jù)現(xiàn)場選育情況減少選留，慎重選配，經(jīng)過3個(gè)世代左右可育成抗斷奶前仔豬腹瀉純系.而FUT1基因的檢測中，有利等位基因頻率大于0. 75，與其他商業(yè)種豬群體的檢測結(jié)果相比偏高［12-13］，這與溫氏杜洛克群體多年的群體選育息息相關(guān)，斷奶后仔豬腹瀉情況較輕微.2個(gè)基因均與父系指數(shù)、母系指數(shù)極顯著相關(guān)，說明這2個(gè)基因在溫氏杜洛克群體中具備應(yīng)用于輔助選擇的潛力.

通過分析2個(gè)抗腹瀉基因組合基因型分布情況發(fā)現(xiàn):在溫氏杜洛克群體中，GG/GA型組合(254頭)所占群體比例最高達(dá)42. 05%.各種組合基因型與各生產(chǎn)性狀的相關(guān)分析顯示，GG/GA組合型個(gè)體的綜合生產(chǎn)性能表現(xiàn)最好，這與本文的單基因相關(guān)分析結(jié)果一致.由于2個(gè)基因有利基因型所占比例相差較大，在兼顧性能測定和表型的現(xiàn)場實(shí)踐選育的情況下，2個(gè)基因聚合育種需要經(jīng)過多個(gè)世代才能實(shí)現(xiàn)優(yōu)化.鑒于此，在杜洛克群體中，MUC13基因的優(yōu)勢等位基因頻率近90%，建議先實(shí)現(xiàn)MUC13基因的輔助選育純化，再啟動(dòng)FUT1基因的選育工作，相對來講，較同時(shí)選育更容易育成抗腹瀉純合品系.

3.2溫氏皮特蘭核心群體的抗腹瀉基因分子輔助選擇方案

溫氏皮特蘭核心群體MUC13基因的表現(xiàn)與杜洛克群體完全不同，皮特蘭MUC13優(yōu)勢等位基因頻率僅0. 180，這與當(dāng)時(shí)的產(chǎn)房生產(chǎn)狀況相吻合，斷奶前腹瀉仔豬偏多.MUC13基因型與生產(chǎn)性能相關(guān)分析結(jié)果顯示，有利基因型與體高、頭型評分、肢蹄評分等體型外貌性狀呈極顯著負(fù)相關(guān)，僅與100 kg瘦肉率生產(chǎn)性狀呈正相關(guān).根據(jù)目前市場需求分析制定的皮特蘭選育方案中，體型外貌為重點(diǎn)選育方向.因此，綜合考慮群體優(yōu)勢等位基因頻率低及與體型外貌的負(fù)相關(guān)等因素，MUC13基因暫不列為皮特蘭群體的選育候選基因.

FUT1基因在皮特蘭群體中的表現(xiàn)與在杜洛克群體中的表現(xiàn)相差不大，優(yōu)勢等位基因頻率達(dá)0. 677，其中雜合子個(gè)體偏多.FUT1基因型與生產(chǎn)性能相關(guān)分析結(jié)果顯示，優(yōu)勢基因型(AA)在體高、頭型評分等體型外貌方面表現(xiàn)不如GG型個(gè)體，但在背膘EBV和出生體質(zhì)量方面表現(xiàn)出正相關(guān)趨勢.綜合現(xiàn)場育種實(shí)踐，F(xiàn)UT1基因在皮特蘭群體中的優(yōu)化較難實(shí)現(xiàn)，基礎(chǔ)群雜合子過多造成優(yōu)化周期加長，且基因型選擇與現(xiàn)場育種實(shí)踐指標(biāo)不一致(體型外貌)或相關(guān)性不高(父系指數(shù)等)，因此，皮特蘭群體暫時(shí)不宜對FUT1基因進(jìn)行純化，僅需保持對基因頻率進(jìn)行長期監(jiān)控即可，在市場以及群體頻率變化適宜時(shí)，再開展優(yōu)化工作.

3.3關(guān)于基因聚合育種的群體特異性

在本研究中，同為終端父本的杜洛克和皮特蘭群體的抗腹瀉基因頻率表現(xiàn)不同，基因型和生產(chǎn)性狀的相關(guān)分析也表現(xiàn)出較大差異.由于核心育種群都是以固定的育種目標(biāo)長期選育的群體，因此同樣的基因在不同的核心育種群的頻率分布和對各生產(chǎn)性狀的影響表現(xiàn)不同.若是多個(gè)基因聚合在不同的育種群中，可能由于不同基因頻率、不同的基因型組合分布、不同相關(guān)性、不同連鎖關(guān)系以及不同選育方向等綜合因素表現(xiàn)出更加復(fù)雜的特異性，所以基因聚合育種不能通用任何群體，而是需要針對不同群體進(jìn)行分析和驗(yàn)證，具有群體特異性.在本研究中2個(gè)抗腹瀉基因輔助選擇育種，可應(yīng)用于溫氏杜洛克群體，卻并不適用于皮特蘭群體，與2個(gè)群體長期向著不同的選育方向選育有關(guān).其基因型組合分布呈現(xiàn)很大差異，根據(jù)相關(guān)分析及現(xiàn)場育種的實(shí)踐出發(fā)，制定適宜的抗腹瀉基因的分子輔助選擇方案.

致謝:感謝江西農(nóng)業(yè)大學(xué)省部共建豬遺傳改良與養(yǎng)殖技術(shù)國家重點(diǎn)實(shí)驗(yàn)室任軍教授和黃路生院士給予的支持和幫助!

參考文獻(xiàn):

［1］施啟順.豬腸毒性大腸桿菌(ETEC)病抗病育種研究［J］.中國畜牧獸醫(yī)，1999，26(4):51-54.

［2］許建民，趙德明.仔豬腸毒素性大腸桿菌病的防治［J］.中國獸醫(yī)雜志，2005，41(2):48-50.

［3］劉勇，顧正明，王克虎，等.微生態(tài)制劑促康生Ⅱ號(hào)治療大腸桿菌抗藥性菌株引起的仔豬黃痢白痢臨床效果觀察［J］.中國獸醫(yī)雜志，2001，37(7):14.

［4］曹國文，姜永康，周淑蘭，等.中草藥抗豬大腸桿菌制劑的研制與效果觀察［J］.黑龍江畜牧獸醫(yī)，2002，26 (10):33-34.

［5］李書華，李艷琴，馮敏，等.生態(tài)制劑——制痢康防治仔豬下痢的試驗(yàn)研究［J］.黑龍江畜牧獸醫(yī)，1999，23 (5):19-20.

［6］SCHAEFER L R.Strategy for applying genome-wide selection in dairy cattle［J］.J Anim Breed Genet，2006，123(4):218-223.

［7］V?GELI P，BERTSCHINGER H U，STAMM M，et al.Genes specifying receptors for F18 fimbriated Escherichia coli，causing oedema disease and postweaning diarrhoea in pigs，map to chromosome 6［J］.Anim Genet，1996，27 (5):321-328.

［8］ZHANG B，REN J，YAN X M，et al.Investigation of the porcine MUC13 gene: Isolation，expression，polymorphisms and their strong association with susceptibility to enterotoxigenic E.coli F4ab/ac［J］.Anim Genet，2008，39(3):258-266.

［9］REN J，YAN X M，AI H S，et al.Susceptibility towards enterotoxigenic Escherichia coli F4ac diarrhea is governed by the MUC13 gene in pigs［J/OL］.PloS One，2012，7 (9):e44573［2014－11－11］.http: / /journals.plos.org/plosone/article? id = 10.1371/journal.pone.0044573.

［10］GOETSTOUWERS T，POUCKE M V，COPPIETERS W，et al.Refined candidate region for F4ab/ac enterotoxigenic Escherichia coli susceptibility situated proximal to MUC13 in pigs［J/OL］.PloS One，2014，9(8):e105013 ［2014－11－11］.http: / /journals.plos.org/plosone/article? id =10.1371/journal.pone.0105013.

［11］阮國榮，肖石軍，徐盼，等.福建省商業(yè)豬種MUC13、IGF2和RYR1基因主效位點(diǎn)的遺傳變異分析［J］.江西農(nóng)業(yè)大學(xué)學(xué)報(bào)，2012，34(5):997-1002.

［12］CODDENS A，VERDONCK F，MULINGE M，et al.The possibility of positive selection for both F18 + Escherichia coli and stress resistant pigs opens new perspectives for pig breeding［J］.Vet Microbiol，2008，126 (1/2/3):210-215.

［13］晏學(xué)明，郭源梅，丁能水，等.不同品種豬a1巖藻糖轉(zhuǎn)移酶基因遺傳變異初析［J］.中國畜牧雜志，2004，40 (3):8-10.

【責(zé)任編輯柴焰】

Molecular maker-assisted selections of MUC13 and FUT1 genes in the two swine nucleus populations

YANG Ming1，2，WANG Qinglai1，2，LIU Jingshun1，2，LIU Zhenyun1，2，WEN Shuxian1，WU Zhenfang1，2，3，LUO Xufang1，3(1 Guangdong Wens Foodstuff Co.，Ltd.，Xinxing 527439，China;
2 National Engineering Research Center for Breeding Swine Industry，Xinxing 527439，China;
3 College of Animal Science，South China Agricultural University，Guangzhou 510642，China)

Abstract:【Objective】MUC13 and FUT1 genes are associated with susceptibility/ resistance to ETEC diarrhea in pre-weaning and weaning piglets，respectively.The goal of this study was to improve breeding of diarrhea resistance in different swine populations through development of MAS technology on the two genes.【Method】In the study，PCR-SNaPshot and PCR-RFLP assay were used for animals genotyping from two swine populations(666 Duroc and 512 Pietrain pigs) for the 2 specialized variants.【Result and conclusion】The results showed that the favorable allele frequency (FAF) of MUC13 gene was 0. 890 and 0. 180 in Duroc and Pietrain respectively while the FAF of FUT1 gene was 0. 754 and 0. 677 respectively.Animals with a favorable allele combination of MUC13 and FUT1 genes were different in Duroc with abook=2,ebook=274frequency of 36. 75% and Pietrain with 3. 49%.The correlation between genotype and phenotype was analyzed，and the results showed that the gene effect on traits was different in the two populations.For example，there were significantly negative correlations between genotype and somatotype in Pietrain，but no association was found in Duroc.According to the breeding practice，it is recommended to achieve purification of MUC13 gene first，and then follow FUT1 gene in Duroc.It is not suitable for marker-assisted selection of two genes in Pietrain temporarily because the ratio of animals with a favorable allele combination of MUC13 and FUT1 gene is too low.

Key words:MUC13 gene; FUT1 gene; swine nucleus population; molecular marker-assisted selection

基金項(xiàng)目:國家高技術(shù)研究發(fā)展計(jì)劃(2011AA100304) ;廣東省科技計(jì)劃項(xiàng)目(2011A202102003) ;廣東省院士工作站建設(shè)項(xiàng)目(2011A020102003)

作者簡介:楊明(1985—)，女，博士，E-mail: 359311126 @ qq.com;通信作者:羅旭芳(1964—)，男，副教授，碩士，E-mail: luoxf@ wens.com.cn

收稿日期:2014-11-24

文章編號(hào):1001-411X(2015) 06-0001-08

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

中圖分類號(hào):S813. 3