王聰亮 張政軒 付琪 魏宇新 白晶晶 宋曉越 李隴平 屈雷 朱海鯨
摘要 [目的]探究催乳素受體(prolactin receptor,PRLR)基因InDel突變對絨山羊種群生長特性的影響。[方法]以陜北白絨山羊為研究對象,采用PCR擴增和測序技術(shù)法檢測分析陜北白絨山羊 PRLR 基因的多態(tài)性,利用獨立樣本 t 檢驗和單因素方差法分析 PRLR基因 不同基因型與陜北白絨山羊生長性狀之間的相關(guān)性。[結(jié)果]在所檢測的1 176只陜北白絨山羊群體中, PRLR 基因第2內(nèi)含子存在1個16 bp InDel突變,形成了純合插入型(II)、雜合型(ID)與純合缺失型(DD)3種基因型;在育成羊群體中,該InDel突變與生長性狀無顯著關(guān)聯(lián)( P >0.05),在成年羊群體中,InDel突變與管圍性狀顯著相關(guān)( P< 0.05);在育成羊和成年羊全部群體中,InDel突變?nèi)耘c管圍性狀顯著相關(guān)( P< 0.05)。[結(jié)論] PRLR 基因16 bp InDel突變可能會對陜北白絨山羊生長發(fā)育產(chǎn)生顯著影響,該位點可作為陜北白絨山羊輔助育種的候選分子標(biāo)記。
關(guān)鍵詞 絨山羊;PRLR基因;InDel;生長性狀;關(guān)聯(lián)分析
中圖分類號 S827? 文獻標(biāo)識碼 A? 文章編號 0517-6611(2022)11-0085-04
doi:10.3969/j.issn.0517-6611.2022.11.022
開放科學(xué)(資源服務(wù))標(biāo)識碼(OSID):
InDel Detection of ?PRLR ?Gene and Its Association Analysis with Growth Traits of Shanbei White Cashmere Goats
WANG Cong-liang,ZHANG Zheng-xuan,F(xiàn)U Qi et al
(Shaanxi Province Engineering & Technology Research Center of Cashmere Goats,Yulin University,Yulin,Shaanxi 719000)
Abstract [Objective]In order to investigate the effect of prolactin receptor(PRLR) gene InDel mutation on the growth characteristics of cashmere goat population.[Method]The Shanbei white cashmere goats were used as the research object,and used PCR amplification and direct sequencing techniques to detect the polymorphism of ?PRLR ?gene in Shanbei white cashmere goats,and analyzed the correlation between different genotypes of ?PRLR ?gene and growth traits of Shanbei white cashmere goats using independent sample ?t -test and one-way ANOVA.[Result]The results showed that there was one 16 bp InDel mutation in intron 2 of the ?PRLR ?gene in the 1 176 Shanbei white cashmere goats population,and there were three genotypes:homozygous insertion type(II),heterozygous type(ID) and homozygous deletion type(DD),respectively.The InDel mutation was not significantly correlated with growth traits in the yearling goats population( P >0.05),and was significantly correlated with the cannon circumference trait in the adult goats population( P< 0.05);In all of the goat population this InDel mutation was still significantly correlated with the cannon circumference trait( P< 0.05).[Conclusion]The 16 bp InDel mutation in ?PRLR ?gene may impose significant effect on the growth and development of Shanbei white cashmere goats,and this locus can be used as a candidate molecular marker for the auxiliary breeding of high-yielding Shanbei white cashmere goats.
Key words Cashmere goats; PRLR ?gene;InDel;Growth traits;Association analysisB38838FF-7371-4061-8A06-23ACEAFF21C7
催乳素(Prolactin,PRL)主要是垂體前葉分泌的一種肽類激素[1],僅包含1個跨膜區(qū)域的膜結(jié)合蛋白,最早發(fā)現(xiàn)于乳腺、卵巢、睪丸等組織,PRL通過與其受體(PRLR)結(jié)合作用于靶器官,參與調(diào)控乳腺發(fā)育、黃體生成、調(diào)節(jié)滲透壓等及繁殖相關(guān)激素分泌等多種生理活動[2-3]。按照PRLR受體胞內(nèi)域氨基酸數(shù)目和胞質(zhì)區(qū)域組成可分為長型PRLR、中型PRLR和短型PRLR,其功能各不相同[4]。據(jù)報道,PRL在大鼠、綿羊和牛等多個物種內(nèi)分布廣泛,通過與不同PRLR受體亞型結(jié)合介導(dǎo),作用于相應(yīng)靶細胞質(zhì)膜,進而發(fā)揮其生物學(xué)作用[5-6]。
PRLR 基因作為泌乳相關(guān)基因,研究人員發(fā)現(xiàn), PRLR基 因與母豬總產(chǎn)仔數(shù)和產(chǎn)活仔數(shù)顯著相關(guān),可作為豬高產(chǎn)仔數(shù)性狀的候選基因[7],國內(nèi)外研究人員就 PRLR 基因?qū)π笄菝谌楹头敝承誀畹挠绊懻归_大量研究。如, PRLR 基因側(cè)翼區(qū)一處SNP突變與疆岳驢平均日泌乳量和乳蛋白率顯著相關(guān),可作為乳用型疆岳驢選育的有效DNA標(biāo)記[8]。此外, PRLR基因 也可作為高產(chǎn)仔數(shù)和海門山羊選育的重要候選基因[9-10]。但 PRLR 基因多態(tài)并不是都與動物繁殖相關(guān),如 PRLR基 因多個SNP突變位點與湖羊和濟寧青山羊等數(shù)個品種的產(chǎn)羔性狀無顯著關(guān)聯(lián)[11-12]。陜北白絨山羊作為優(yōu)秀絨肉兼用品種,目前關(guān)于 PRLR基 因與陜北白絨山羊生產(chǎn)相關(guān)性狀的研究較少[13]。為此,筆者采用PCR擴增和測序技術(shù)探究 PRLR 基因InDel多態(tài)性及其與不同陜北白絨山羊群體生產(chǎn)性狀的相關(guān)性,以期為絨山羊的高效選育提供參考。
1 材料與方法
1.1 試驗動物與基因組DNA
試驗羊由陜西應(yīng)馬安絨山羊養(yǎng)殖公司提供,包括618只育成羊(1.0~1.5歲)和558只成年羊(1.5~4.0歲),均為雌性。剪取試驗羊耳尖部約0.4 cm2耳組織儲存在含有70 %乙醇的1.5 mL Eppedorff管中,同時測量試驗羊體重、體高、體長、胸圍、髖寬、管圍、十字部高、胸深、胸寬、絨細度10個性狀數(shù)據(jù)[14]。同時提取試驗羊耳組織DNA[15],-80 ℃保存?zhèn)溆谩?/p>
1.2 主要試劑 瓊脂糖、Marker 1和PCR Master Mix,均購自天根生化科技(北京)有限公司;EB核酸染料,購自生工生物工程(上海)股份有限公司。
1.3 ?PRLR 基因InDel目的片段的PCR擴增 參照文獻[13],選取 PRLR 基因第2內(nèi)含子16 bp InDel突變位點(NC_030827.1.g.61200-61206,61209-61214,61298-61300 del GAGGGC-GGGAGGG-AAG)作為候選變異位點,并利用文獻中的引物對候選變異位點進行多態(tài)性檢測。 PRLR 基因的PCR擴增體系:PCR Master Mix 6.5 μL,上下游引物各0.3 μL,模板DNA 0.7 μL,ddH2O補至13 μL。PCR擴增程序:95 ℃預(yù)變性5 min;94 ℃變性30 s,56 ℃退火30 s,72 ℃延伸20 s,共38個循環(huán);72 ℃再延伸10 min,4 ℃結(jié)束反應(yīng)保存,擴增目的片段長度分別為143 bp和127 bp。PCR擴增產(chǎn)物進行3%瓊脂糖凝膠電泳并鑒定分型,引物合成與PCR產(chǎn)物測序均由生工生物工程(上海)股份有限公司完成。
1.4 數(shù)據(jù)統(tǒng)計
利用Microsoft Excel 2018軟件計算 PRLR 基因InDel突變位點的基因型頻率和等位基因頻率,使用Nei方法計算基因純合度(Ho)、雜合度(He)、有效等位基因數(shù)(Ne)和多態(tài)信息含量(PIC)等群體遺傳學(xué)參數(shù),使用在線平臺SHEsis(https://analysis.bio-x.cn)計算 PRLR 基因多態(tài)性是否處于哈代溫伯格(HWE)平衡[16]。其生物統(tǒng)計模型為 Yijk=μ+Gi+Eij,其中μ為總體平均值,Gi為基因型固定效應(yīng),Eij 為隨機誤差。利用SPSS 23.0軟件,分別運用獨立樣本 t 檢驗和單因素方差方法完成 PRLR基 因陜北白絨山羊育成羊和成年羊群體不同基因型與生長性狀的關(guān)聯(lián)分析;再利用Microsoft Excel 2018軟件計算不同群體各個性狀表型值的
平均值,用個體表型值減去平均值得到差值,運用單因素方
差方法分析不同基因型與陜北白絨山羊群體生長性狀之間的差異顯著性,其中育成羊和成年羊為一個整體,結(jié)果以“平均值±標(biāo)準(zhǔn)誤”表示。
2 結(jié)果與分析
2.1 ?PRLR 基因InDel檢測與分型 擴增產(chǎn)物經(jīng)3.5%瓊脂糖凝膠電泳,電泳結(jié)果見圖1。陜北白絨山羊 PRLR 基因16 bp InDel位點(NC_030827.1.g.61200-61206,61209-61214,61298-61300 del GAGGGC-GGGAGGG-AAG)具有多態(tài)性,分別產(chǎn)生了純合插入型(II,143 bp)、雜合型(ID,143 bp和127 bp)與純合缺失型(DD,127 bp)3種基因型,產(chǎn)物測序結(jié)果見圖2,測序結(jié)果與PCR擴增結(jié)果一致。
2.2 ?PRLR 基因遺傳多樣性 ?PRLR 基因16 bp InDel突變位點等位基因頻率、基因型頻率、多態(tài)性信息含量(PIC)等群體遺傳學(xué)參數(shù)見表1。結(jié)果顯示, PRLR 基因16 bp InDel突變位點在陜北白絨山羊育成羊、成年羊和全部群體中均屬于低度多態(tài)(PIC<0.25),I型等位基因頻率均高于D型;該突變位點在成年羊群體處于哈代溫伯格平衡狀態(tài)( P >0.05),但在育成羊和全部群體中偏離了哈代溫伯格平衡( P< 0.05)。
2.3 不同陜北白絨山羊群體 PRLR 基因多態(tài)性與生長性狀的相關(guān)性分析B38838FF-7371-4061-8A06-23ACEAFF21C7
通過SPSS 23.0軟件對 PRLR 基因InDel位點多態(tài)性與不同陜北白絨山羊群體生長性狀進行相關(guān)性分析,結(jié)果見表2~4。在陜北白絨山羊育成羊群體中 (n =618),由于體重和絨細度性狀DD基因型只有2個樣本,故運用獨立樣本 t 檢驗法對其進行相關(guān)性分析,結(jié)果顯示不同基因型與生長性狀均無顯著關(guān)聯(lián)(表2);在成年羊群體中 (n =558),ID基因型山羊的管圍顯著高于II基因型( P< 0.05)(表3);在育成羊和成年羊全部群體中( n =1 176),ID基因型山羊的管圍仍顯著高于II基因型( P< 0.05),優(yōu)勢基因型與成年羊群體相同,均為ID型(表4)。
3 討論
陜北白絨山羊作為國內(nèi)知名絨肉兼用山羊品種,其中陜北地區(qū)是重要養(yǎng)殖區(qū)域,產(chǎn)肉性能及遺傳穩(wěn)定性差等是養(yǎng)殖過程中的主要問題[17]。相對傳統(tǒng)育種方法而言,分子標(biāo)記輔助選擇技術(shù)(MAS)具有選擇準(zhǔn)確度高、育種周期短等優(yōu)勢[18],目前已廣泛應(yīng)用于家畜育種,用于改良家畜經(jīng)濟性狀,提高生產(chǎn)性能[19-20]。
多態(tài)信息含量等與基因遺傳多樣性往往具有直接關(guān)系[21]。該突變多態(tài)信息含量顯示,在陜北白絨山羊不同群體中均屬于低度多態(tài)水平(PIC<0.25),說明群體遺傳變異較小,可適當(dāng)加強選擇強度,但在育成羊和全部群體中偏離了哈代溫伯格平衡( P< 0.05),說明可能受到遺傳漂變及人工選育的影響,而陜北白絨山羊正是人工選育與雜交形成的新品種,等位基因流失可能是其主要影響因素[22]。將 PRLR 基因16 bp InDel突變位點多態(tài)性與陜北白絨山羊的生長性狀進行相關(guān)性分析發(fā)現(xiàn),該突變位點與陜北白絨山羊育成羊生長性狀無顯著關(guān)聯(lián),但與成年羊管圍性狀顯著相關(guān);當(dāng)把育成羊和成年羊作為一個群體分析時,該InDel突變位點仍與管圍性狀顯著相關(guān),優(yōu)勢基因型與成年羊群體相同,均為ID型,具有很好的一致性。在成年羊群體中,該突變位點能顯著影響山羊的管圍性狀,但對育成羊生長性狀無顯著影響,推測可能存在以下原因:①生長發(fā)育階段不同,育成羊尚未發(fā)育完善;② PRLR 基因表達存在時序性,如 PRLR 基因在牦牛生長發(fā)育后期睪丸組織中的表達量極顯著高于前期[23],同時5日齡與8周齡雌性大鼠肝臟組織中 PRLR 基因表達量有較大差異[24],故 PRLR 基因可能在陜北白絨山羊不同生長發(fā)育階段存在差異表達。同時研究顯示[25],羊體質(zhì)量與管圍性狀存在顯著關(guān)聯(lián),表明 PRLR 基因該突變位點可能會間接影響陜北白絨山羊體重性狀,其中具體機制仍需要進一步探究。
有報道稱, PRLR 基因分別定位于綿羊、山羊的16、20號染色體[26-27],其中在綿羊腹股溝竇、脂肪和山羊肌肉、肝臟及皮下脂肪中存在不同程度的 PRLR 基因表達[28-29]。目前研究人員普遍認(rèn)為 PRLR 主要通過JAK2-STAT信號轉(zhuǎn)導(dǎo)途徑參與機體調(diào)控,通過JAK2酪氨酸激酶的活化對STAT5轉(zhuǎn)錄因子進行磷酸化,從而形成同型二聚體,與靶基因啟動子上的γ-干擾素激活位點序列特異性結(jié)合,激活基因轉(zhuǎn)錄,誘導(dǎo)并調(diào)控蛋白表達[20-31],除此之外,有絲分裂原激活蛋白(MAP)激酶途徑也被認(rèn)為是參與介導(dǎo)PRL調(diào)控作用的重要途徑[32]。研究顯示,催乳素合成分泌與下丘腦多巴胺的調(diào)控有直接關(guān)系,多巴胺與催乳素細胞表面的多巴胺D2特異性結(jié)合,在小鼠中,缺乏多巴胺D2R受體亞型與其體重增加有直接關(guān)系[33],小鼠敲除 PRLR 基因后骨形成速率和骨密度較敲除前有所下降[34]。研究證實, PRLR 在調(diào)控棕色脂肪細胞分化過程中發(fā)揮重要作用[35],在刺激NIH3T3細胞快速增殖的同時還能有效促進細胞分裂蛋白激酶的激活[36]。同時研究人員發(fā)現(xiàn), PRLR在 牦牛等動物肝臟、骨骼肌和腎臟均有較高表達水平[23,37],并 且PRLR 基因p.H406R和p.M480T 2個錯義突變位點多態(tài)與大白豬背部脂肪有顯著影響[38],說 明PRLR 基因可能參與早期動物胚胎及機體的生長發(fā)育。
綜上所述, PRLR 基因16 bp InDel突變對陜北白絨山羊管圍性狀有顯著影響,可作為陜北白絨山羊優(yōu)良生長性狀選育的重要候選基因。
4 結(jié)論
該研究驗證了 PRLR 基因在陜北白絨山羊群體中存在1個16 bp的插入缺失,并發(fā)現(xiàn)該突變位點對陜北白絨山羊群體部分生長性狀存在顯著關(guān)聯(lián),提示該位點可作為陜北白絨山羊分子標(biāo)記輔助育種的有效DNA標(biāo)記。
參考文獻
[1]
邵順成,李新海,閆背背,等.5個綿羊群體PRLR基因多態(tài)性與產(chǎn)羔數(shù)的關(guān)聯(lián)分析[J].中國草食動物科學(xué),2021,41(4):14-18.
[2] BACHELOT A,BINART N.Reproductive role of prolactin[J].Reproduction,2007,133(2):361-369.
[3] TERMAN A.Effect of the polymorphism of prolactin receptor(PRLR) and leptin(LEP) genes on litter size in Polish pigs[J].J Anim Breed Genet,2005,122(6):400-404.
[4] 張芬芬,王長法,黃金明,等.催乳素受體的研究進展[J].家畜生態(tài)學(xué)報,2009,30(3):6-9.
[5] DALRYMPLE A,EDERY M,JABBOUR H N.Sequence and functional characterisation of the marmoset monkey( Callithrix jacchus) ?prolactin receptor:Comparative homology with the human long-form prolactin receptor[J].Mol Cell Endocrinol,2000,167(1/2):89-97.B38838FF-7371-4061-8A06-23ACEAFF21C7
[6] COSENZA G,IANNACCONE M,AUZINO B,et al.Remarkable genetic diversity detected at river buffalo prolactin receptor( PRLR) ?gene and association studies with milk fatty acid composition[J].Anim Genet,2018,49(3):159-168.
[7] VINCENT A L,TUGGLE C K,ROTHSCHILD M F,et al.The prolactin receptor gene is associated with increased litter size in pigs[J].Swine Res Rep,1998,27:15-18.
[8] 畢蘭舒,肖海霞,臧長江,等.疆岳驢催乳素受體基因多態(tài)性及其與泌乳性狀間的相關(guān)性分析[J].中國畜牧獸醫(yī),2017,44(1):180-185.
[9] 范一萍,王曉梅,賈啟濤,等.豬PRLR和FSHβ基因多態(tài)性檢測及其與繁殖性狀的關(guān)聯(lián)分析[J].中國畜牧獸醫(yī),2019,46(3):774-781.
[10] LI Y J,ZHANG L,SHANG L Q,et al.Genetic polymorphisms at three loci of ?PRLR ?and ?FSHR ?gene correlate with litter size in Chinese haimen goat[J].J Anim Vet Adv,2010,9(22):2835-2838.
[11] 郭海燕,金銀,李擁軍,等. PRLR 基因多態(tài)性與湖羊產(chǎn)羔性能的關(guān)系[J].揚州大學(xué)學(xué)報(農(nóng)業(yè)與生命科學(xué)版),2017,38(1):57-61.
[12] 張圓,袁婧,高一珊,等.濟寧青山羊 PRLR 基因的多態(tài)性及其與主要經(jīng)濟性狀的關(guān)聯(lián)分析[J].安徽農(nóng)業(yè)科學(xué),2020,48(16):94-97.
[13] LIU X F,MA L,WANG M,et al.Two indel variants of prolactin receptor( PRLR) ?gene are associated with growth traits in goat[J].Anim Biotechnol,2020,31(4):314-323.
[14] WANG Z,ZHANG X,JIANG E,et al.InDels within caprine ?IGF2BP1 ?intron 2 and the 3′-untranslated regions are associated with goat growth traits[J].Anim Genet,2020,51(1):117-121.
[15] LAN X Y,PAN C Y,CHEN H,et al.An AluI PCR-RFLP detecting a silent allele at the goat POU1F1 locus and its association with production traits[J].Small Ruminant Res,2007,73(1/2/3):8-12.
[16] CUI Y,YAN H L,WANG K,et al.Insertion/deletion within the ?KDM6A ?gene is significantly associated with litter size in goat[J].Front Genet,2018,9:1-11.
[17] 周寒,陳繼明,劉瑞霞,等.陜北白絨山羊產(chǎn)業(yè)發(fā)展現(xiàn)狀調(diào)查研究[J].畜牧與飼料科學(xué),2012,33(4):64-65.
[18] 馬彬云,吳建平.遺傳標(biāo)記技術(shù)在動物育種中的研究進展[J].甘肅農(nóng)業(yè)大學(xué)學(xué)報,2004,39(1):92-96,100.
[19] YANG Q,ZHANG S H,LIU L L,et al.Application of mathematical expectation(ME) strategy for detecting low frequency mutations:An example for evaluating 14-bp insertion/deletion(indel) within the bovine ?PRNP ?gene[J].Prion,2016,10(5):409-419.
[20] 田佺,王思寒,劉磊,等.延邊黃牛 MYH3和MYH8基 因多態(tài)性及其與生長性狀的關(guān)聯(lián)分析[J].中國畜牧獸醫(yī),2021,48(6):2122-2128.
[21] 唐曉琴,陳平博,袁婷婷,等.同羊 ETAA1 基因InDel檢測及其與生長性狀的關(guān)聯(lián)分析[J].安徽農(nóng)業(yè)大學(xué)學(xué)報,2021,48(2):245-249.
[22] 劉燕,李學(xué)武,王瑞軍,等.周歲內(nèi)蒙古絨山羊不同部位羊絨長度、細度變化規(guī)律及其相關(guān)性分析[J].黑龍江畜牧獸醫(yī),2019(17):58-62.
[23] 馬富龍,褚敏,常永芳,等. LHR、PRLR基 因在不同發(fā)育階段牦牛睪丸中的表達與定位[J].農(nóng)業(yè)生物技術(shù)學(xué)報,2020,28(4):672-680.
[24] TANAKA M,SUZUKI M,KAWANA T,et al.Differential effects of sex steroid hormones on the expression of multiple first exons including a novel first exon of prolactin receptor gene in the rat liver[J].J Mol Endocrinol,2005,34(3):667-673.B38838FF-7371-4061-8A06-23ACEAFF21C7
[25] 趙森,王真,白洋洋,等.絨山羊 MARCH1 基因InDel位點與體尺性狀的關(guān)聯(lián)研究[J].家畜生態(tài)學(xué)報,2020,41(3):27-32.
[26] JENKINS Z A,HENRY H M,SISE J A,et al.Follistatin(FST),growth hormone receptor(GHR) and prolactin receptor(PRLR) genes map to the same region of sheep chromosome 16[J].Anim Genet,2000,31(4):280.
[27] MOODY D E,POMP D,BARENDSE W,et al.Assignment of the growth hormone receptor gene to bovine chromosome 20 using linkage analysis and somatic cell mapping[J].Anim Genet,1995,26(5):341-343.
[28] ALEXANDRE-PIRES G,MARTINS C,GALVO A M,et al.Understanding the inguinal sinus in sheep( Ovis aries) -morphology,secretion,and expression of progesterone,estrogens,and prolactin receptors[J].Int J Mol Sci,2017,18(7):1-20.
[29] SHI H P,ZHANG T Y,YI Y Q,et al.Long form PRLR (lPRLR) regulates genes involved in the triacylglycerol synthesis in goat mammary gland epithelial cells[J].Small Rumin Res,2016,139:7-14.
[30] RADHAKRISHNAN A,RAJU R,TULADHAR N,et al.A pathway map of prolactin signaling[J].J Cell Commun Signal,2012,6(3):169-173.
[31] GALLEGO M I,BINART N,ROBINSON G W,et al.Prolactin,growth hormone,and epidermal growth factor activate Stat5 in different compartments of mammary tissue and exert different and overlapping developmental effects[J].Dev Biol,2001,229(1):163-175.
[32] DOBOLYI A,OLH S,KELLER D,et al.Secretion and function of pituitary prolactin in evolutionary perspective[J].Front Neurosci,2020,14:1-16.
[33] PEREZ-MILLAN M I,LUQUE G M,RAMIREZ M C,et al.Selective disruption of dopamine D2 receptors in pituitary lactotropes increases body weight and adiposity in female mice[J].Endocrinology,2014,155(3):829-839.
[34] KELLY P A,BINART N,F(xiàn)REEMARK M,et al.Prolactin receptor signal transduction pathways and actions determined in prolactin receptor knockout mice[J].Biochem Soc Trans,2001,29(Pt2):48-52.
[35] KIM H,TOYOFUKU Y,LYNN F C,et al.Serotonin regulates pancreatic beta cell mass during pregnancy[J].Nat Med,2010,16(7):804-808.
[36] DAS R,VONDERHAAR B K.Transduction of prolactin's (PRL) growth signal through both long and short forms of the PRL receptor[J].Mol Endocrinol,1995,9(12):1750-1759.
[37] 邢光東.鵝催乳素受體基因克隆和多克隆抗體制備及組織表達分析[D].南京:南京農(nóng)業(yè)大學(xué),2011.
[38] FONTANESI L,GALIMBERTI G,CAL D G,et al.Identification and association analysis of several hundred single nucleotide polymorphisms within candidate genes for back fat thickness in Italian Large White pigs using a selective genotyping approach[J].J Anim Sci,2012,90(8):2450-2464.B38838FF-7371-4061-8A06-23ACEAFF21C7