祖盤玉，李維，林家棟，李洪林，簡華峰，劉洋，牟騰慧，龍廣麗，張福平*

赤水烏骨雞基因多態(tài)性及其生物信息學(xué)分析

祖盤玉1,2,3,4,5，李維6，林家棟1,2,3,4,5，李洪林1,2,3,4,5，簡華峰1,2,3,4,5，劉洋1,2,3,4,5，牟騰慧1,2,3,4,5，龍廣麗1,2,3,4,5，張福平1,2,3,4,5*

(1.貴州大學(xué)動物科學(xué)學(xué)院，貴州 貴陽 550025；2.高原山地動物遺傳育種與繁殖教育部重點實驗室，貴州 貴陽 550025；3.貴州省動物遺傳育種與繁殖重點實驗室，貴州 貴陽 550025；4.貴州大學(xué)科研雞場，貴州 貴陽 550025；5.貴州大學(xué)家禽研究所貴州 貴陽 550025；6.貴州省畜禽遺傳資源管理站，貴州 貴陽550001)

采用DNA混池及PCR產(chǎn)物直接測序技術(shù)，對赤水烏骨雞、泰和烏雞、興義白雞、興義白雞F2代、貴州黃快羽雞、貴州黃慢羽雞、羅曼蛋雞和瑤山雞8個雞種群基因外顯子區(qū)域進(jìn)行多態(tài)性分析。結(jié)果顯示：瑤山雞與紅色原雞序列一致，無SNPs，赤水烏骨雞的4個SNPs分別為G23A(Arg→His)、C69T、T212C、G274A，泰和烏雞的8個SNPs分別為C69T、T212C、G274A、T398C(Leu→Pro)、G636A、T637C、A644C(His→Pro)、C834T，貴州黃快、慢羽雞和興義白雞均有3個SNPs，分別為A427G(Thr→Ala)、G636A、T637C，羅曼蛋雞的4個SNPs分別為T398A(Leu→Gln)、G636A、T637C、C834T，興義白雞F2代的4個SNPs分別為C69T、T212C、G274A、A644C(His→Thr)；生物信息學(xué)分析發(fā)現(xiàn)，除G274A、T398C、T398A、A644C位點的穩(wěn)定性增加外，其余位點的穩(wěn)定性均有所降低；除G274A、T398C、T637C、A644C位點的整體多樣性有所下降外，其余位點的均有所增加；除泰和烏雞基因編碼蛋白為不穩(wěn)定蛋白外，其余種群的均屬于穩(wěn)定蛋白，且包括泰和烏雞在內(nèi)的7個雞群基因編碼蛋白均為疏水性蛋白和非分泌蛋白；赤水烏骨雞、泰和烏雞、興義白雞、興義白雞F2代和羅曼蛋雞的MC1R蛋白三級結(jié)構(gòu)均由α–螺旋、β–轉(zhuǎn)角、無規(guī)則卷曲組成。

赤水烏骨雞；泰和烏雞；興義白雞；興義白雞F2代；貴州黃快羽雞；貴州黃慢羽雞；羅曼蛋雞；瑤山雞；基因；SNPs；羽色；脛色；生物信息學(xué)分析

赤水烏骨雞是貴州優(yōu)良的肉蛋兼用型品種。母雞羽毛為純黑色，公雞的背部、頸部和鐮羽多帶紅色羽絲，尾部羽毛油黑帶墨綠色光澤，其余部分為純黑色，公母雞的皮、肉、喙、脛、爪均為烏黑色。赤水烏骨雞體型較大、耐熱、耐粗飼、覓食能力強(qiáng)、肉質(zhì)鮮美、營養(yǎng)豐富[1]，與被國際上列為烏雞標(biāo)準(zhǔn)品種的泰和烏雞相比，其皮膚、胸肌、腿肌、舌頭等組織烏度均淺于泰和烏雞。研究表明，烏骨雞烏度、藥用價值及其獨特的滋補(bǔ)功效均與其所含的黑色素有關(guān)[2–4]。而黑素皮質(zhì)素受體1(melanocortin 1 receptor，)基因是影響黑色素合成的主控基因[5–6]。目前，多數(shù)研究[7–9]均集中在其SNPs與羽色及脛色上。TAKEUCHI等[10]和KERJE等[11]的研究結(jié)果表明，的突變與雞的羽色有關(guān)。劉騏嘉[12]的研究結(jié)果表明，基因型和單倍型在不同脛色表型雞中存在著明顯差異。為了進(jìn)一步研究基因的SNPs除與羽色和脛色相關(guān)外，是否還是影響赤水烏骨雞與泰和烏雞皮膚、胸肌、腿肌、舌頭、骨骼等組織烏度差異的原因，本研究中，以黑羽赤水烏骨雞及白絲羽泰和烏雞為主要研究對象，同時選取全黃羽的貴州黃快羽雞、貴州黃慢羽雞(黃色脛)和全白羽的羅曼蛋雞(白色脛)、興義白雞(青色脛)、興義白羽雞F2代(興義白雞(母本)與羅曼蛋雞(父本)雜交后再自交，白色脛)及麻羽的瑤山雞(青色脛)作為對照，運用DNA混池結(jié)合PCR產(chǎn)物直接測序法及生物信息學(xué)法對基因在不同羽色雞群中的SNPs進(jìn)行分析，以期找到赤水烏骨雞和泰和烏雞基因關(guān)鍵SNPs，為后期赤水烏骨雞的選育、黑色素沉積規(guī)律的研究及其藥用價值的開發(fā)利用提供依據(jù)。

1　材料與方法

1.1　材料

1) 試驗動物。赤水烏骨雞血樣144份、泰和烏雞血樣98份，貴州黃快羽雞、貴州黃慢羽雞、興義白雞、羅曼蛋雞(粉殼)、興義白雞F2代等雞群血樣公母各30份。以上雞群為貴州大學(xué)科研雞場同條件飼養(yǎng)至150日齡的健康雞?，幧诫u150日齡血樣采集于貴陽綠源禽業(yè)有限公司。血樣均用翅下靜脈肝素鈉抗凝采血法采集。

2) 試驗試劑。瓊脂糖、PCR Master Mix和Ezup柱式動物基因組DNA抽提試劑盒均購自上海生工生物工程技術(shù)服務(wù)有限公司。

1.2　方法

1) 基因組DNA的提取。采用DNA試劑盒提取DNA，用紫外分光光度計測定總DNA的濃度和純度，利用1%瓊脂糖凝膠電泳，隨機(jī)抽取檢測其DNA提取效果。

2) 引物的設(shè)計與合成。參照劉騏嘉[12]的研究設(shè)計基因引物：F，5–GCTTTGTAGGT GCTGCAGTTGTG–3；R，5–CCATCCATCCTCCT GTCTGT–3。擴(kuò)增長度為1 048 bp。

3) PCR擴(kuò)增和SNPs鑒定及其等位基因頻率估算。分別取各雞群雞只的DNA 1 μL，分別構(gòu)建各雞群DNA混池，進(jìn)行PCR擴(kuò)增。PCR擴(kuò)增反應(yīng)體系：體積為30 μL，dd H2O 7.5 μL，PCR Master Mix 15 μL，DNA 3 μL，上、下游引物各2.25 μL。PCR反應(yīng)條件：94 ℃預(yù)變性5 min；95 ℃變性30 s，63 ℃退火30 s，72 ℃延伸30 s，共34個循環(huán)；72 ℃延伸10 min。利用1%瓊脂糖凝膠電泳檢測其產(chǎn)物，由上海生工生物工程技術(shù)服務(wù)有限公司進(jìn)行產(chǎn)物純化和DNA池雙向測序。運用DNA star和MWsnap 3.0分析NCBI上雞的基因序列(紅色原雞)，即原序列和測序峰圖，并根據(jù)以下公式估算等位基因頻率。

A=B/(1+2)(=1,2)。

式中：A表示SNP位點某等位基因頻率；1與2分別表示測序峰圖上該等位基因1、2峰的高度。

1.3　MC1R基因生物信息學(xué)分析

利用在線軟件RNAfold web server(http://nibiru. tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi)預(yù)測基因突變前后編碼區(qū)序列的mRNA二級結(jié)構(gòu)；運用ExPASy–ProtParam tool (https://web. expasy.org/protparam/)預(yù)測蛋白質(zhì)理化性質(zhì)；采用TMHMM Server.2.0(http://www.cbs.dtu.dk/services/ TMHMM/)分析MC1R蛋白跨膜區(qū)域；利用SignalP 4.0 Server(http://www.cbs.dtu.dk/services/SignalP/)分析MC1R蛋白信號肽；運用NPS@: SOPMA secondary structure prediction (https://npsa-prabi.ibcp.fr/cgi-bin/ npsa_automat.pl?page=npsa_sopma.html)預(yù)測和分析MC1R蛋白的二級結(jié)構(gòu)；在線軟件SWISS–MODEL (https://swissmodel.expasy.org/)預(yù)測MC1R蛋白的三級結(jié)構(gòu)。

2　結(jié)果與分析

2.1　PCR檢測結(jié)果

用含核酸染料的1%瓊脂糖凝膠電泳檢測PCR產(chǎn)物，結(jié)果見圖1。

M　Marker DL2000；1~8分別為赤水烏骨雞、泰和烏雞、貴州黃快、慢羽、羅曼蛋雞、興義白雞、興義白雞F2代、瑤山雞。

2.2　MC1R基因SNP位點鑒定

利用DNAStar，將赤水烏骨雞及其他7個雞群基因雙向測序結(jié)果與GenBank上提交的紅色原雞基因序列(GenBank登錄號為NC_006098) (后簡稱原序列)進(jìn)行比對，結(jié)果發(fā)現(xiàn)，赤水烏骨雞有4個SNPs(以基因的第一外顯子第一個堿基為+1位)，分別為G23A(Arg→His)、C69T、T212C、G274A；泰和烏雞有8個SNPs，分別為C69T、T212C、G274A、T398C(Leu→Pro)、G636A、T637C、A644C(His→Pro)、C834T；貴州黃快、慢羽雞和興義白雞均有3個SNPs，分別為A427G(Thr→Ala)、G636A、T637C；羅曼蛋雞有4個SNPs，分別為T398A(Leu→Gln)、G636A、T637C、C834T；興義白雞F2代有4個SNPs，分別為C69T、T212C、G274A、A644C(His→Thr)；瑤山雞的與原序列一致，無SNPs。測序結(jié)果如圖2所示。特征序列列于表1。

圖2　不同雞種群MC1R基因PCR產(chǎn)物測序結(jié)果

表1　不同雞種群MC1R基因的SNPs位點及序列

2.3　MC1R基因的SNPs等位基因頻率估算及mRNA二級結(jié)構(gòu)預(yù)測結(jié)果

由表2可知，除泰和烏雞C69T、T212C、T637C 3個SNP位點等位基因頻率突變前大于突變后外，其余SNP位點等位基因頻率均為突變前小于突變后；除G274A、T398C、T398A、A644C位點的最低自由能有所降低外，其余位點的最低自由能均有所增加，即G274A、T398C、T398A、A644C位點穩(wěn)定性增加，其余位點穩(wěn)定性均有所降低；除G274A、T398C、T637C、A644C位點的整體多樣性有所下降外，其余位點的整體多樣性均有所增加，而除A644C和T398C位點外，其余發(fā)生異義突變的G23A、A427G、T398A的整體多樣性均有所增加。

表2　不同雞種群MC1R基因的SNPs等位基因頻率及mRNA二級結(jié)構(gòu)預(yù)測

表中SNP位點基因頻率是SNP位點標(biāo)示位置后面的堿基頻率，如G23A的基因頻率是堿基A的頻率。

2.4　MC1R蛋白質(zhì)序列分析

2.4.1MC1R蛋白的理化性質(zhì)

由表3可知，除羅曼蛋雞的相對分子質(zhì)量比原序列的高外，其余種群的相對分子質(zhì)量均低于原序列；除赤水烏骨雞的理論等電點低于原序列的外，其余種群的理論等電點與原序列的相同；除瑤山雞與紅色原雞外，其余7個種群的總原子數(shù)較原序列均有所減少；除泰和烏雞屬于不穩(wěn)定蛋白外，其余種群均屬于穩(wěn)定蛋白；本研究的8個雞種群MC1R蛋白均屬于疏水性蛋白。

表3　不同雞種群MC1R蛋白的理化性質(zhì)

2.4.2MC1R蛋白跨膜區(qū)域及信號肽分析

MC1R蛋白的7個跨膜結(jié)構(gòu)區(qū)域與原序列相比，均屬于跨膜蛋白，且7個跨膜區(qū)的氨基酸殘基個數(shù)及連接跨膜區(qū)的區(qū)域長度也均未改變，其中在細(xì)胞膜外有G23A和C69T，在細(xì)胞膜內(nèi)有T212C、A427G、G636A、T637C、A644C，細(xì)胞膜上有G274A、T398C和C834T。對8個種群的MC1R蛋白信號肽分析，發(fā)現(xiàn)8個種群均屬于非分泌蛋白。

2.4.3MC1R蛋白的二級結(jié)構(gòu)預(yù)測

由表4可知，8個種群的MC1R蛋白二級結(jié)構(gòu)均有α–螺旋、無規(guī)則卷曲、β–轉(zhuǎn)角及擴(kuò)展鏈，而無β–折疊，且α–螺旋均大于45%，故MC1R蛋白屬于全α蛋白。

表4　不同雞種群MC1R蛋白的二級結(jié)構(gòu)

2.4.4MC1R蛋白的三級結(jié)構(gòu)預(yù)測

利用SWISS–MODEL預(yù)測MC1R蛋白的三級結(jié)構(gòu)，發(fā)現(xiàn)其由α–螺旋、β–轉(zhuǎn)角、無規(guī)則卷曲構(gòu)成(圖3)。

紅色示α–螺旋；黃色示β–轉(zhuǎn)角；綠色示無規(guī)則卷曲；藍(lán)色示跨膜區(qū)。

3　結(jié)論與討論

本研究中，在除瑤山雞外的7個種群中共找到基因10個SNPs，分別是G23A、C69T、T212C、G274A、T398C/A、A427G、G636A、T637C、A644C、C834T。

G23A突變位點僅在赤水烏骨雞中被發(fā)現(xiàn)，且引起了基因mRNA二級結(jié)構(gòu)、氨基酸、蛋白質(zhì)理化性質(zhì)、蛋白質(zhì)二級及三級結(jié)構(gòu)的改變。此位點，在其他雞品種還未見報道，有可能是影響赤水烏骨雞和泰和烏雞烏度差異的原因之一。

本研究中，在貴州黃快、慢羽品系、羅曼蛋雞、興義白雞和瑤山雞中均未發(fā)現(xiàn)有C69T突變位點，而在赤水烏骨雞、泰和烏雞和興義白雞F2代中均發(fā)現(xiàn)C69T。通過測序峰圖估算發(fā)現(xiàn)，赤水烏骨雞僅出現(xiàn)TT型，這與張濤等[13]研究的略陽雞花羽品系、郭秀麗[14]研究的海蘭褐和羅曼褐及陳多珍等[15]研究的狄高肉雞D品系一致，但在略陽黑羽烏雞中，卻發(fā)現(xiàn)2種基因型，其中TT型和CT型分別占比42%和58%[16]。泰和烏雞、略陽普通雞及河北柴雞的黑羽、灰麻和白羽系的中勢等位基因均為C，且在河北柴雞中各基因型在各種羽色中的差異極顯著(<0.01)[14]，而興義白雞F2代的為T。此外，泰和烏雞和赤水烏骨雞均為青色脛，興義白雞F2代為白色脛。劉騏嘉[12]研究發(fā)現(xiàn)，CC型在黃色脛和豆綠脛雞中為優(yōu)勢等位基因型，TT型在青色脛雞中為優(yōu)勢等位基因型。由此推測，在不同雞品種中，C69T的3種基因型與羽色及脛色關(guān)聯(lián)性不一，但T等位基因與青色脛關(guān)聯(lián)性較強(qiáng)。

本研究中，在興義白雞、羅曼蛋雞、貴州黃雞和瑤山雞中均未發(fā)現(xiàn)T212C突變位點，而在泰和烏雞中，T為優(yōu)勢等位基因。這與略陽黑羽烏雞[13]和他留烏骨雞白羽、麻羽、黑羽品系[16]及河北柴雞黑羽、白羽、灰麻羽、黃麻羽品系[14]的表現(xiàn)一致，而赤水烏骨雞和興義白雞F2代中C均為優(yōu)勢等位基因，而在略陽普通雞[13]和狄高肉雞D品系[15]僅發(fā)現(xiàn)CC型，略陽雞花羽烏雞中僅發(fā)現(xiàn)CT型[13]。楊慧[17]利用TalL酶切位點研究金水烏雞F2代中T212C表明，TT型在黑羽和花系羽中未發(fā)現(xiàn)，CT型僅出現(xiàn)在花羽系，而CC型在黑羽、白羽和花羽中均有發(fā)現(xiàn)。此外，在劉騏嘉[12]的研究中，青色脛中C均為優(yōu)勢等位基因，赤水烏骨雞與其一致。由于目前報道中，未涉及到白色脛，T212C的3種基因型是否與白色脛有關(guān)聯(lián)性，還需進(jìn)一步研究。

本研究中，僅赤水烏骨雞、泰和烏雞和興義白雞F2代發(fā)現(xiàn)有G274A突變位點，且泰和烏雞和赤水烏雞與他留烏骨雞黑羽系[16]、海蘭褐[14]、羅曼褐[14]及狄高肉雞D品系[15]均只發(fā)現(xiàn)AA型。在略陽黑羽烏雞中，卻發(fā)現(xiàn)GG型和GA型2種基因型，且GG型為優(yōu)勢等位基因型[13]。根據(jù)測序峰圖的估算，興義白雞F2代可能出現(xiàn)GG、GA和AA 3種中的2種或3種基因型。田寬校[18]利用MscI限制性內(nèi)切酶對純合黑羽靈山土雞和溫氏401黃雞的G274A位點進(jìn)行酶切后，認(rèn)為黑色素擴(kuò)散基因是由基因來編碼的，但基因又不是決定靈山土雞羽色的唯一基因。此外，本研究中的赤水烏骨雞和泰和烏雞及劉騏嘉[12]研究的青色脛雞中，優(yōu)勢等位基因型均為AA型。由此推測，在不同雞品種中，G274A的3種基因型與羽色及脛色關(guān)聯(lián)性不一，但A等位基因與青色脛關(guān)聯(lián)性較強(qiáng)，與白色脛關(guān)聯(lián)性如何還需要進(jìn)一步研究。

本研究中的泰和烏雞與河北柴雞[14]黑羽、灰麻、黃羽系均發(fā)生T398C，且TT型為優(yōu)勢等位基因型，而在白羽系中未發(fā)現(xiàn)此突變，其中CC型僅在黑羽系中發(fā)現(xiàn)，TC型在灰麻系的2個世代中均有發(fā)現(xiàn)，而在黃麻品系中，僅在1世代中有發(fā)現(xiàn)。本研究中的羅曼蛋雞與冉金山[19]研究的羅曼蛋雞均有T398A。冉金山[19]研究發(fā)現(xiàn)，在羅曼蛋雞(白羽)、廣元灰雞(灰羽)及蘆花雞(橫斑羽)中，優(yōu)勢等位基因均為T，而本研究中的羅曼蛋雞優(yōu)勢等位基因為A，可能是由于樣本量較小所致。就T398C/A而言，在同品種中，CC型與黑羽關(guān)聯(lián)性較大，而在不同雞品種中，是否與羽色相關(guān)聯(lián)還需進(jìn)一步研究。由于T398C(Leu→Pro)發(fā)生異義突變，此位點可能是影響赤水烏骨雞和泰和烏雞烏度差異的原因之一。

A427G突變位點僅出現(xiàn)在貴州黃快、慢羽系和興義白雞中，且只發(fā)現(xiàn)GG型，這與陳多珍等[15]對狄高肉雞AB和C品系的研究發(fā)現(xiàn)一致。在郭秀麗[14]研究的河北柴雞黃麻系中，經(jīng)過一個世代的提純后，GG基因型頻率有所提高，但在黑羽系、白羽系及灰麻系零世代中均未出現(xiàn)此基因型，而在灰麻系1世代中又發(fā)現(xiàn)此基因型。在郭秀麗[14]研究的海蘭褐和羅曼褐中也只出現(xiàn)AG和GG型。此外，郭秀麗[14]對A427G在脛色上的關(guān)聯(lián)性研究發(fā)現(xiàn)，黃色脛雞中，G為優(yōu)勢等位基因，且青脛和黃脛、黃脛和豆綠脛的等位基因頻率差異都達(dá)到了極顯著水平(<0.01)，而青脛和豆綠脛的差異不顯著[12]。由此推測，A427G的G等位基因可能與黃羽及黃色脛關(guān)聯(lián)性較大，但與白羽及白色脛關(guān)系如何還需要進(jìn)一步研究。本研究中，G636A突變位點在貴州黃快、慢羽系及羅曼蛋雞中均只發(fā)現(xiàn)AA型，與狄高肉雞AB和C品系[15]中發(fā)現(xiàn)的一致。郭秀麗[14]研究發(fā)現(xiàn)，GA和AA型在其所研究的所有雞品種中均有發(fā)現(xiàn)，經(jīng)過1個世代的提純后，GG型和AA型在黑羽和灰羽系中均有所降低，而GA型在黑羽、灰羽和黃麻品系中均有所增加，AA型僅在黃麻系中有所增加。本研究的泰和烏雞A等位基因頻率為0.53，而他留烏骨雞白羽、麻羽、黑羽3個品系與鄭嫩珠等[20]研究的C品系(黑羽，皮膚、肌肉、骨膜、心臟、肝臟等組織均為黑色)和D品系(白羽，皮膚、肌肉、骨膜、心臟、肝臟等組織均為白色)的一致，G為優(yōu)勢等位基因。此外，青色脛的泰和烏雞的優(yōu)勢等位基因為A，這與劉騏嘉[12]的研究結(jié)果(青脛雞中G為優(yōu)勢等位基因)不一致，可能是由于品種差異所致。劉騏嘉[12]在黃脛和豆綠脛的研究中，A為優(yōu)勢等位基因，本研究對貴州黃雞和羅曼蛋雞的研究結(jié)果與其一致。由此推測，G636A的3種基因型與羽色關(guān)聯(lián)性不大，但A等位基因與黃色脛關(guān)聯(lián)性較大，可對此位點與脛色進(jìn)行關(guān)聯(lián)性分析。

與陳多珍[15]研究的狄高肉雞AB和C品系一致，本研究中，貴州黃快、慢羽和羅曼蛋雞的T637C突變位點均只發(fā)現(xiàn)CC型。在本研究的泰和烏雞的T637C中，除CC基因型外還存在其他基因型，且C為優(yōu)勢等位基因，這與郭秀麗[14]研究的海蘭褐、羅曼褐、河北柴雞的黑羽和黃麻系及冉金山[19]研究的羅曼蛋雞的一致。在滕召純等[16]研究的他留烏骨雞白羽、麻羽、黑羽3個品系和郭秀麗[14]研究的河北柴雞白羽系中，T為優(yōu)勢等位基因。鄭嫩珠等[20]對德化黑雞的研究中發(fā)現(xiàn)，在C品系和D品系中均只有TT和CC 2種基因型，C品系中有TT型(93.18%)和CC型(6.82%)，D品系中有TT型(92%)和CC型(8%)，其羽色依舊是TT型為優(yōu)勢等位基因型。本研究中，青色脛的泰和烏雞優(yōu)勢等位基因為C，而劉騏嘉[12]研究發(fā)現(xiàn)，青色脛雞中T為優(yōu)勢等位基因，可能是由于品種差異所致；但劉騏嘉[12]研究發(fā)現(xiàn)在黃色脛雞中，C為優(yōu)勢等位基因，本研究的貴州黃雞和羅曼蛋雞與其一致。由此推測，在同品種中，T637C的C等位基因與黃色羽和黃色脛關(guān)聯(lián)性較強(qiáng)。

本研究中，A644C突變位點僅在泰和烏雞和興義白雞F2代中發(fā)現(xiàn)，其優(yōu)勢等位基因分別為A和C。略陽普通雞為CC型，略陽花羽烏雞為AC型，從略陽黑羽烏雞中發(fā)現(xiàn)有AA型和AC型2種基因型，且AA型為優(yōu)勢等位基因型[13]。在他留烏骨雞中僅從白羽和麻羽品系中發(fā)現(xiàn)有A644C，且A均為優(yōu)勢等位基因，而在黑羽品系中未發(fā)現(xiàn)此突變位點[16]。郭秀麗[14]研究發(fā)現(xiàn)，AA型在其所研究的所有雞品種中均為優(yōu)勢等位基因型，AC型在除白羽系外的其他雞品種中均有所發(fā)現(xiàn)，AA型在黑羽和灰羽系中均有所降低，而CC型在黑羽和黃麻系中均有所增加，但未在海蘭褐和羅曼褐中發(fā)現(xiàn)CC型[14]。在狄高肉雞中，AA型出現(xiàn)在AB和C品系中，CC型僅出現(xiàn)在D品系中[15]。此外，本研究中，泰和烏雞的優(yōu)勢等位基因為A，而青色脛雞的優(yōu)勢等位基因為C[12]。由此推測，A644C在不同雞品種中，不同基因型與羽色和脛色性狀關(guān)聯(lián)性不一，在同一品種中，與脛色有一定的關(guān)聯(lián)性。

本研究中，僅從羅曼蛋雞和泰和烏雞中發(fā)現(xiàn)有C834T突變位點，二者優(yōu)勢等位基因分別為T和C。楊永升[7]研究發(fā)現(xiàn)，C834T位點在C系為絲羽烏骨雞親本(白羽、黑皮膚、活體與屠體脛色黑色、黑肉、黑色內(nèi)臟膜)中發(fā)現(xiàn)CC、CT和TT 3種基因型，CC型為優(yōu)勢等位基因型，而在A品系為明星肉雞親本(白羽、白皮膚、內(nèi)臟膜無色、肉色微紅、活體脛色為黃色，屠體脛色白色)中未發(fā)現(xiàn)此突變位點。此外，在以絲羽烏骨雞與明星肉雞為親本建立的中國農(nóng)業(yè)大學(xué)資源家系群體中，C834T突變與雞的活體脛色性狀顯著相關(guān)(<0.05)，在其他組織黑色素性狀中差異不顯著。王建等[21]研究發(fā)現(xiàn)CC和TT基因型在不同活體脛色性狀中差異顯著(<0.05)。由此推測，C834T可能與白色羽和脛色性狀有很大的相關(guān)性。

基因是一個高度多態(tài)的基因，對不同雞品種羽色和脛色的影響并不是某一單個位點所決定的，而是多個突變位點共同作用所致。單對本研究的赤水烏骨雞和泰和烏雞而言，通過生物信息學(xué)分析發(fā)現(xiàn)，赤水烏骨雞的G23A和泰和烏雞的T398C、A644C發(fā)生異義突變，從而導(dǎo)致MC1R蛋白質(zhì)二級及三級結(jié)構(gòu)的改變，而泰和烏雞T398C和A644C同時突變會導(dǎo)致MC1R蛋白穩(wěn)定性降低，變?yōu)椴环€(wěn)定蛋白。由此推測，G23A、T398C和A644C可能是影響赤水烏骨雞和泰和烏雞烏度差異的重要原因。此外，基因[22]和基因[23]分別與橫斑羽及金銀羽有關(guān)，且都位于Z染色體上。郭軍等[24]研究表明，Z染色體基因?qū)Φ半u早期體質(zhì)量加性遺傳方差有貢獻(xiàn)，而赤水烏骨雞公雞的背部、頸部和鐮羽還帶有紅色羽絲，母雞則全是黑羽，且母雞膚色比公雞烏度深，性別也可作為后期選育烏度更深的赤水烏骨雞的一項指標(biāo)。

[1] 肖竺宏．貴州赤水竹鄉(xiāng)雞的遺傳資源調(diào)查[J]．貴州畜牧獸醫(yī)，2012，36(4)：27–28．

XIAO Z H．Investigation on genetic resources of Chishui black-bone chickens in Guizhou[J]．Guizhou Journal of Animal Husbandry & Veterinary Medicine，2012，36(4)：27–28．

[2] 鄭從義，胡國斌，屈三甫，等．黑色素對流感病毒誘導(dǎo)細(xì)胞凋亡的抑制效應(yīng)[J]．武漢大學(xué)學(xué)報(自然科學(xué)版)，1999，45(4)：469–472．

ZHENG C Y，HU G B，QU S F，et al．Inhibition of melanin fromon the apoptosis induced by influenza virus in MDCK cells[J]．Wuhan University Journal(Natural Science Edition)，1999，45(4)：469–472．

[3] GANTZ I，SHIMOTO Y，KONDA Y，et al．Molecular cloning，expression，and characterization of a fifth melanocortin receptor[J]．Biochemical & Biophysical Research Communications，1994，200(3)：1214–1220．

[4] MOUNTJOY K G，ROBBINS L S，MORTRUD M T，et al．The cloning of a family of genes that encode the melanocortin receptors[J]．Science，1992，257：1248–1251．

[5] 王樂，張斌，鄭文新，等．動物毛色與黑色素皮質(zhì)素受體1()基因[J]．草食家畜，2009(2)：10–12．

WANG L，ZHAGN B，ZHENG W X，et al．The hair color of animal and melanocortin receptor 1()[J]. Grass-Feeding Livestock，2009(2)：10–12．

[6] 徐幸蓮，莊蘇，陳伯祥．烏骨雞黑色素對延緩果蠅衰老的作用[J]．南京農(nóng)業(yè)大學(xué)學(xué)報，1999，22(2)：105–108．

XU L，ZHUANG S，CHEN B X．The effect of melanin from White Silky Fowl on antiageing in[J]．Journal of Nanjing Agricultural University，22(2)：105–108．

[7] 楊永升．雞基因的SNPs分析與染色體定位[D]．杭州：浙江大學(xué)，2004．

YANG Y S．SNPs analysis and chromosomal location ofgene in chicken[D]．Hangzhou：Zhejiang University，2004．

[8] RODRIGUES A R，ALMEIDA H，GOUVEIA A M. Intracellular signaling mechanisms of the melanocortin receptors：current state of the art[J]．Cellular and Molecular Life Sciences，2015，72(7)：1331–1345．

[9] 封竣淇，徐偉，黃蘭，等．基因的研究進(jìn)展[J]．中國畜牧獸醫(yī)，2017(4)：1141–1148．

FENG J Q，XU W，HUANG L，et al．Research advance ofgene[J]．China Animal Husbandry & Veterinary Medicine，2017，44(4)：1141–1148．

[10] TAKEUCHI S，SUZUKI H，HIROSE S，et al. Molecular cloning and sequence analysis of the chick melanocortin 1-receptor gene[J]．Biochimica et Biophysica Acta (BBA)-Gene Structure and Expression，1996，1306(2/3)：122–126．

[11] KERJE S，LIND J，SCHUTZ，K，et al. Melanocortin 1-receptor () mutations are associated with plumage colour in chicken[J]．Animal Genetics，2003，34(4)：241–248．

[12] 劉騏嘉．不同脛色優(yōu)質(zhì)雞生產(chǎn)性能比較及其與基因多態(tài)性的相關(guān)分析[D]．雅安：四川農(nóng)業(yè)大學(xué)，2015．

LU Q J．Comparison of different shank color quality chicken production performance and associated with theSNP analysis[D]．Ya’an：Sichuan Agricultural University，2015．

[13] 張濤，路宏朝，江海，等．略陽烏雞黑素皮質(zhì)素受體1()基因變異及蛋白質(zhì)結(jié)構(gòu)的預(yù)測[J]．貴州農(nóng)業(yè)科學(xué)，2012，40(10)：120–125．

ZHANG T，LU H Z，JIANG H，et al．Variation ofgene and prediction of protein structure in Lueyang Silky Fowl[J]．Guizhou Agricultural Sciences，2012，40(10)：120–125．

[14] 郭秀麗．河北柴雞基因變異研究[D]．保定：河北農(nóng)業(yè)大學(xué)，2009．

GUO X L．Study on variation ofgene Hebei Domestic Chicken[D]．Baoding：Hebei Agricultural University，2009．

[15] 陳多珍，李艷梅，楊向東，等．基因在狄高肉雞配套系羽色純化中的應(yīng)用[J]．中國家禽，2018，40(18)：8–12．

CHEN D Z，LI Y M，YANG X D，et al．Application ofgene in purification of feather color of Diga Broiler Lines[J]．China Poultry，2018，40(18)：8–12．

[16] 滕召純，陸嘵屏，王玉祥，等．基因ⅠPCR– RFLP標(biāo)記與他留烏骨雞羽色性狀的相關(guān)性研究[J]．中國家禽，2013，35(20)：11–14．

TENG Z C，LU X P，WANG Y X，et al．Association ofⅠ PCR-RFLP aenetic marker with plumage color in Taliu Black-Boned Chicken [J]. China poultry，2013，35(20)：11–14．

[17] 楊慧．金水絲羽烏骨雞資源群體的構(gòu)建和三個基因的SNPs檢測及其與羽色性狀的連鎖分析[D]．武漢：華中農(nóng)業(yè)大學(xué)，2010．

YANG H．Construction of resource population for Jinshui Silky Chicken and SNPs detection of three genes and their linkage analysis with plumage colour[D]．Wuhan：Huazhong Agricultural University，2010．

[18] 田寬校．靈山土雞基因的多態(tài)性與羽色的相關(guān)性研究[D]．南京：南京農(nóng)業(yè)大學(xué)，2009．

TIAN K X．Association of PCR-RFLP polymorphism ofgene with plumage colour in Lingshan Local Chicken[D]．Nangjing：Nanjing Agricultural University，2009．

[19] 冉金山．雞羽色相關(guān)基因、和的多態(tài)性及組織表達(dá)研究[D]．雅安：四川農(nóng)業(yè)大學(xué)，2017．

RAN J S．Study on polymorphisms and tissues expression of，andgenes in different feather color chicken[D]．Ya’an：Sichuan Agricultural University，2017．

[20] 鄭嫩珠，董曉寧，陳暉，等．德化黑雞黑素皮質(zhì)素受體1()基因單核苷酸多態(tài)性分析[J]．福建農(nóng)業(yè)學(xué)報，2008，23(2)：132–136．

ZHENG N Z，DONG X N，CHEN H，et al．Single nucleotide polymorphism analysis on melanocortin 1-receptor gene in Dehua Black Chicken[J]．Fujian Journal of Agricultural Sciences，2008，23(2)：132–136．

[21] 王建，鐘發(fā)剛，吳桐忠．基因多態(tài)性與其膚色、脛色表型相關(guān)性的分析[J]．國外畜牧學(xué)(豬與禽)，2008，28(2)：75–77．

WANG J，ZHONG F G，WU T Z．Correlation betweengene polymorphism and skin color and shin phenotype [J]．Pigs and Poultry，2008，28(2)：75–77．

[22] HELLSTR?M A R，SUNDSTR?M E，GUNNARSSON U，et al．Sex-linked barring in chickens is controlled by thetumour suppressor locus[J]．Pigment Cell Melanoma Res，2010，23(4)：521–530．

[23] GUNNARSSON U，HELLSTR?M A R，TIXIER- BOICHARD M，et al．Mutations incause plumage color variation in chicken and Japanese quail[J]．Genetics，2007，175(2)：867–877．

[24] 郭軍，曲亮，竇套存，等．3周齡蛋雞體質(zhì)量性連鎖遺傳效應(yīng)分析[J]．湖南農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)版)，2018，44(6)：645–649．

GUO J，QU L，DOU T C，et al．Effects of sex linkage genetic variances on chick’s body weight at third week age[J]．Journal of Hunan Agricultural University(Natural Sciences)，2018，44(6)：645–649．

Polymorphism and bioinformatics analysis ofgene in Chishui Black-Bone chickens

ZU Panyu1,2,3,4,5, LI Wei6, LIN Jiadong1,2,3,4,5, LI Honglin1,2,3,4,5, JIAN Huafeng1,2,3,4,5, LIUYang1,2,3,4,5, MOU Tenghui1,2,3,4,5, LONG Gangli1,2,3,4,5, ZHANG Fuping1,2,3,4,5*

(1.College of Animal Sciences, Guizhou University, Guiyang, Guizhou 550025, China; 2.Key Laboratory of Animal Genetics, Breeding and Reproduction in the Plateau Mountainous Region, Ministry of Education, Guiyang, Guizhou 550025, China; 3.Guizhou Key Laboratory of Animal Genetics, Breeding and Reproduction, Guiyang, Guizhou 550025, China; 4.Chicken Farm of Guizhou University, Guiyang, Guizhou 550025, China; 5.Poultry Research Institute of Guizhou University, Guiyang, Guizhou 550025, China; 6.Guizhou Province Animal and Poultry Genetic Resources Management Station, Guiyang, Guizhou 550001, China)

DNA mixing pool and direct sequencing of PCR products were used to analyze the polymorphism ofexon regions in 8 chicken populations of Chishui Black-Bone Chickens, Taihe Silky Fowl, Xingyi White Chickens, the F2Xingyi White Chickens, Guizhou Yellow Fast Chickens, Guizhou Yellow Slow Chickens, Luoman Laying Hens and Yaoshan Chickens. The results showed that the sequence of Yaoshan Chickens was identical to that of Red Jungle Fowl without SNPs. The four SNPs of Chishui Black-Bone Chickens were G23A(Arg→His), C69T, T212C and G274A. The eight SNPs of Taihe Silky Fowl were C69T, T212C, G274A, T398C(Leu→Pro), G636A, T637C, A644C(His→Pro) and C834T. Guizhou Yellow Fast Chickens, Guizhou Yellow Slow Chickens and Xingyi White Chickens all had 3 SNPs, A427G(Thr→Ala), G636A and T637C. The 4 SNPs of Luoman Laying Hens were T398A(Leu→Gln), G636A, T637C and C834T. The 4 SNPs of the F2Xingyi White Chickens were C69T, T212C, G274A and A644C(His→Thr). The bioinformatics analysis found that except for the stability of G274A, T398C, T398A and A644C sites, the stability of the other sites decreased. Except that Taihe Silky Fowl, whosegene encoding protein was unstable, the other populations’s MC1R protein were stable, and all the proteins from the 7 chicken populations were owing hydrophobic and non-secreted character. The tertiary structure of MC1R protein of Chishui Black-Bone Chickens, Taihe Silky Fowl, Xingyi White Chickens, the F2Xingyi White Chickens and Luoman Laying Hens were composed of α-helix, β-turn, and random curl.

Chishui Black-Bone Chickens; Taihe Silky Fowl; Xingyi White Chickens; the F2Xingyi White Chickens; Guizhou Yellow Fast Chickens; Guizhou Yellow Slow Chickens; Luoman laying hen; Yaoshan chicken;gene; SNPs; plumage color; shank color;bioinformatics analysis

10.13,331/j.cnki.jhau.2020.01.013

S831.2

A

1007-1032(2020)01-0084-09

2019–02–22

2019–04–22

貴州省科技計劃項目(黔科合支撐[2016]2507、[2017]2533–1)；貴州省科技計劃項目(2019[2288])；貴州省生態(tài)家禽產(chǎn)業(yè)技術(shù)體系建設(shè)項目

祖盤玉(1992—)，女，貴州盤縣人，碩士研究生，主要從事特種經(jīng)濟(jì)動物飼養(yǎng)研究，1473326412@qq.com；

，張福平副教授，主要從事家禽育種研究，zfu－1010@126.com

祖盤玉，李維，林家棟，李洪林，簡華峰，劉洋，牟騰慧，龍廣麗，張福平．赤水烏骨雞基因多態(tài)性及其生物信息學(xué)分析[J]．湖南農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)版)，2020，46(1)：84–92．

ZU P Y, LI W, LIN J D, LI H L, JIAN H F, LIUY, MOU T H, LONG G L, ZHANG F P. Polymorphism and bioinformatics analysis ofgene in Chishui Black-Bone Chickens[J]．Journal of Hunan Agricultural University(Natural Sciences), 2020, 46(1):84–92.

http://xb.hunau.edu.cn

責(zé)任編輯：鄒慧玲

英文編輯：柳正