赤水烏骨雞MC1R基因多態性及其生物信息學分析

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

赤水烏骨雞基因多態性及其生物信息學分析

祖盤玉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.貴州大學動物科學學院，貴州 貴陽 550025；2.高原山地動物遺傳育種與繁殖教育部重點實驗室，貴州 貴陽 550025；3.貴州省動物遺傳育種與繁殖重點實驗室，貴州 貴陽 550025；4.貴州大學科研雞場，貴州 貴陽 550025；5.貴州大學家禽研究所貴州 貴陽 550025；6.貴州省畜禽遺傳資源管理站，貴州 貴陽550001)

采用DNA混池及PCR產物直接測序技術，對赤水烏骨雞、泰和烏雞、興義白雞、興義白雞F2代、貴州黃快羽雞、貴州黃慢羽雞、羅曼蛋雞和瑤山雞8個雞種群基因外顯子區域進行多態性分析。結果顯示：瑤山雞與紅色原雞序列一致，無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)；生物信息學分析發現，除G274A、T398C、T398A、A644C位點的穩定性增加外，其余位點的穩定性均有所降低；除G274A、T398C、T637C、A644C位點的整體多樣性有所下降外，其余位點的均有所增加；除泰和烏雞基因編碼蛋白為不穩定蛋白外，其余種群的均屬于穩定蛋白，且包括泰和烏雞在內的7個雞群基因編碼蛋白均為疏水性蛋白和非分泌蛋白；赤水烏骨雞、泰和烏雞、興義白雞、興義白雞F2代和羅曼蛋雞的MC1R蛋白三級結構均由α–螺旋、β–轉角、無規則卷曲組成。

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

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

1　材料與方法

1.1　材料

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

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

1.2　方法

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

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

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

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

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

1.3　MC1R基因生物信息學分析

利用在線軟件RNAfold web server(http://nibiru. tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi)預測基因突變前后編碼區序列的mRNA二級結構；運用ExPASy–ProtParam tool (https://web. expasy.org/protparam/)預測蛋白質理化性質；采用TMHMM Server.2.0(http://www.cbs.dtu.dk/services/ TMHMM/)分析MC1R蛋白跨膜區域；利用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)預測和分析MC1R蛋白的二級結構；在線軟件SWISS–MODEL (https://swissmodel.expasy.org/)預測MC1R蛋白的三級結構。

2　結果與分析

2.1　PCR檢測結果

用含核酸染料的1%瓊脂糖凝膠電泳檢測PCR產物，結果見圖1。

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

2.2　MC1R基因SNP位點鑒定

利用DNAStar，將赤水烏骨雞及其他7個雞群基因雙向測序結果與GenBank上提交的紅色原雞基因序列(GenBank登錄號為NC_006098) (后簡稱原序列)進行比對，結果發現，赤水烏骨雞有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。測序結果如圖2所示。特征序列列于表1。

圖2　不同雞種群MC1R基因PCR產物測序結果

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

2.3　MC1R基因的SNPs等位基因頻率估算及mRNA二級結構預測結果

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

表2　不同雞種群MC1R基因的SNPs等位基因頻率及mRNA二級結構預測

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

2.4　MC1R蛋白質序列分析

2.4.1MC1R蛋白的理化性質

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

表3　不同雞種群MC1R蛋白的理化性質

2.4.2MC1R蛋白跨膜區域及信號肽分析

MC1R蛋白的7個跨膜結構區域與原序列相比，均屬于跨膜蛋白，且7個跨膜區的氨基酸殘基個數及連接跨膜區的區域長度也均未改變，其中在細胞膜外有G23A和C69T，在細胞膜內有T212C、A427G、G636A、T637C、A644C，細胞膜上有G274A、T398C和C834T。對8個種群的MC1R蛋白信號肽分析，發現8個種群均屬于非分泌蛋白。

2.4.3MC1R蛋白的二級結構預測

由表4可知，8個種群的MC1R蛋白二級結構均有α–螺旋、無規則卷曲、β–轉角及擴展鏈，而無β–折疊，且α–螺旋均大于45%，故MC1R蛋白屬于全α蛋白。

表4　不同雞種群MC1R蛋白的二級結構

2.4.4MC1R蛋白的三級結構預測

利用SWISS–MODEL預測MC1R蛋白的三級結構，發現其由α–螺旋、β–轉角、無規則卷曲構成(圖3)。

紅色示α–螺旋；黃色示β–轉角；綠色示無規則卷曲；藍色示跨膜區。

3　結論與討論

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

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

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

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

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

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

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

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

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

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

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

[1] 肖竺宏．貴州赤水竹鄉雞的遺傳資源調查[J]．貴州畜牧獸醫，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] 鄭從義，胡國斌，屈三甫，等．黑色素對流感病毒誘導細胞凋亡的抑制效應[J]．武漢大學學報(自然科學版)，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] 王樂，張斌，鄭文新，等．動物毛色與黑色素皮質素受體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]．南京農業大學學報，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]．杭州：浙江大學，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]．中國畜牧獸醫，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] 劉騏嘉．不同脛色優質雞生產性能比較及其與基因多態性的相關分析[D]．雅安：四川農業大學，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] 張濤，路宏朝，江海，等．略陽烏雞黑素皮質素受體1()基因變異及蛋白質結構的預測[J]．貴州農業科學，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]．保定：河北農業大學，2009．

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

[15] 陳多珍，李艷梅，楊向東，等．基因在狄高肉雞配套系羽色純化中的應用[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標記與他留烏骨雞羽色性狀的相關性研究[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] 楊慧．金水絲羽烏骨雞資源群體的構建和三個基因的SNPs檢測及其與羽色性狀的連鎖分析[D]．武漢：華中農業大學，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] 田寬校．靈山土雞基因的多態性與羽色的相關性研究[D]．南京：南京農業大學，2009．

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

[19] 冉金山．雞羽色相關基因、和的多態性及組織表達研究[D]．雅安：四川農業大學，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] 鄭嫩珠，董曉寧，陳暉，等．德化黑雞黑素皮質素受體1()基因單核苷酸多態性分析[J]．福建農業學報，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] 王建，鐘發剛，吳桐忠．基因多態性與其膚色、脛色表型相關性的分析[J]．國外畜牧學(豬與禽)，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周齡蛋雞體質量性連鎖遺傳效應分析[J]．湖南農業大學學報(自然科學版)，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])；貴州省生態家禽產業技術體系建設項目

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

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

祖盤玉，李維，林家棟，李洪林，簡華峰，劉洋，牟騰慧，龍廣麗，張福平．赤水烏骨雞基因多態性及其生物信息學分析[J]．湖南農業大學學報(自然科學版)，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

責任編輯：鄒慧玲

英文編輯：柳正