峨眉黑雞遺傳多樣性及群體遺傳結(jié)構(gòu)分析

袁霞 劉方慶 文隴英 徐婧 廖光祥 王強(qiáng)勝 王湘

（西南山地瀕危鳥類保護(hù)四川省高等院校重點(diǎn)實驗室/樂山師范學(xué)院，四川樂山 614000）

摘 要 峨眉黑雞（Gallus gallus）為原雞屬紅原雞種的一個地方亞種，肉蛋兼用型品種，不僅具有豐富的營養(yǎng)價值，還兼有藥用價值和觀賞價值，是四川峨眉山地區(qū)養(yǎng)殖數(shù)量較多的特有的黑雞優(yōu)秀品種。為了給峨眉黑雞種質(zhì)資源的科學(xué)保存和可持續(xù)性利用提供一定的理論依據(jù)，采集峨眉山地區(qū)三個地點(diǎn)的峨眉黑雞各7只，以線粒體基因Cytb和COI作為分子標(biāo)記，分析種群的遺傳多樣性特征與種群遺傳結(jié)構(gòu)。結(jié)果表明，不同采樣點(diǎn)的峨眉黑雞遺傳多樣性存在差異，峨眉山市鴻飛養(yǎng)雞場的種群具有較高的單倍型多樣性（Hd=0.905），峨眉山市遠(yuǎn)大林下放養(yǎng)雞場的種群序列核苷酸平均配對差異最高（K=5.047 62）。通過分析認(rèn)為，峨眉黑雞的遺傳多樣性并不豐富，建議加強(qiáng)對峨眉黑雞種質(zhì)資源的保護(hù)，保持這一品種的特征性不丟失。

關(guān)鍵詞 峨眉黑雞（Gallus gallus）；遺傳多樣性；系統(tǒng)進(jìn)化；線粒體基因；四川省峨眉山市

中圖分類號：S831；Q311+.3 文獻(xiàn)標(biāo)志碼：A DOI：10.19415/j.cnki.1673-890x.2022.19.001

生物多樣性由遺傳多樣性、物種多樣性和生態(tài)系統(tǒng)多樣性三個部分組成，既包括生物個體及其攜帶的遺傳信息，也包括它們與生態(tài)環(huán)境所組成的生態(tài)系統(tǒng)，以及各組分間的相互聯(lián)系[1]。遺傳多樣性可以作為進(jìn)一步分析生物個體在遺傳進(jìn)化方面的潛能和對環(huán)境的適應(yīng)能力等各方面特征的重要基礎(chǔ)數(shù)據(jù)之一[2]。一個地區(qū)的物種遺傳多樣性越豐富，該物種對環(huán)境的適應(yīng)能力越強(qiáng)，在生存和進(jìn)化方面也具有更大的優(yōu)勢[3]。

地方家禽種質(zhì)遺傳資源是畜禽遺傳多樣性的重要組成部分，是對地方禽類進(jìn)行品種改良和可持續(xù)開發(fā)利用的重要依據(jù)[4]。隨著國內(nèi)經(jīng)濟(jì)的快速發(fā)展，引進(jìn)外來種類（品種）逐漸增多；這些外來優(yōu)良資源在豐富當(dāng)?shù)匦笄萜贩N資源的同時，也會對當(dāng)?shù)卦行笄萜贩N資源多樣性的保護(hù)造成沖擊和破壞[5]。例如我國本土豬在大量引進(jìn)進(jìn)口白豬的情況下，面臨著逐年減少甚至滅絕的情況[6-7]。地方品種雖然在生產(chǎn)性能方面比不上引進(jìn)品種，但其優(yōu)良的肉質(zhì)品質(zhì)和獨(dú)特的風(fēng)味遠(yuǎn)遠(yuǎn)優(yōu)于引進(jìn)品種，尤其是我國地方禽畜品種資源在遺傳育種方面還可能具有較高的研究利用價值[8]。通過研究家禽的遺傳結(jié)構(gòu)多樣性和家禽群體遺傳距離可以初步查明家禽各種類型間的共同起源歷史及親緣關(guān)系，對合理地使用資源及有效保存家禽地方品種資源、建立優(yōu)良基因庫有著重要的參考意義，是家禽遺傳資源管理中的重要環(huán)節(jié)[9]。

峨眉黑雞（Gallus gallus）為原雞屬紅原雞種的一個地方亞種，是一種肉蛋兼用型品種，不僅具有豐富的營養(yǎng)價值，還兼有藥用價值和觀賞價值，是四川峨眉山地區(qū)養(yǎng)殖數(shù)量較多的黑雞優(yōu)秀品種；屬于四川山地烏雞的兩大品系之一[10]。峨眉黑雞主要分布在峨眉山市的龍池、大為、龍門，樂山市沙灣區(qū)的范店、軫溪及峨邊縣的毛坪、共和等鄉(xiāng)鎮(zhèn)；主產(chǎn)區(qū)為峨眉山西南大渡河沿岸的三縣交界處[11]。

峨眉黑雞全基因組重測序?qū)崿F(xiàn)了遺傳進(jìn)化分析及重要性狀候選基因的預(yù)測[12]，但針對其遺傳多樣性及品種形成的歷史研究內(nèi)容還較少。目前，以mtDNA序列為分子標(biāo)記運(yùn)用DNA直接測序的方法，研究家禽遺傳多樣性[13-18]和系統(tǒng)分化[19-22]及起源的研究報道越來越多。本研究通過以峨眉山特有雞種峨眉黑雞作為研究對象，以線粒體DNA細(xì)胞色素c氧化酶亞基I（cytochromecoxidase I， COI）和細(xì)胞色素b（cytochromeb）兩個基因為分子標(biāo)記，通過采集三個地點(diǎn)的峨眉黑雞群體樣本，分析種群遺傳多樣性與遺傳結(jié)構(gòu)，構(gòu)建系統(tǒng)發(fā)育樹，為峨眉黑雞種質(zhì)資源的科學(xué)保存和可持續(xù)性利用提供一定的理論依據(jù)。

1? 材料與方法

1.1? 樣品采集及DNA提取

1.1.1? 樣品采集

21只峨眉黑雞分別來自峨眉山市三個的雞場：峨眉山市世海黑雞原種雞場（E1），峨眉山市遠(yuǎn)大林下放養(yǎng)雞場（E2），峨眉山市鴻飛養(yǎng)雞場（E3）。樣品編號見表1。所有個體用含有20 μL抗凝血劑（EDTA和0.75%生理鹽水按3%～4%比例配制而成）的抗凝管進(jìn)行翅下靜脈采血，采樣體積為0.3 mL，放置于-20 ℃冰箱保存。

1.1.2? DNA提取

使用天根組織/血液試劑盒（北京天根）提取基因組上的DNA序列[22]。

1.2? 引物信息及PCR擴(kuò)增

參照文獻(xiàn)[23]選取合適的引物（見表2）。線粒體DNA 中的Cytb基因擴(kuò)增引物為Cytb-H16064、Cytb-L14731，線粒體DNA中的COI基因擴(kuò)增引物為COI-L6615、COI-H7956；其中L表示輕鏈，H表示重鏈；序列的順序均為5′—3′方向。引物由生工生物工程（上海）股份有限公司合成。

以提取的基因組DNA作為模板，進(jìn)行PCR擴(kuò)增。PCR反應(yīng)體系為25 μL，包括模板DNA 2 μL、ddH2O 8.5 μL、DNA聚合酶（Taq酶）12.5 μL（含dNTP、Mg2+、Buffer緩沖液）、上下游引物各1 μL，反應(yīng)體系在Easycycler Gradient 96型與SelectCycler Ⅱ型PCR儀上進(jìn)行反應(yīng)[24]。

PCR反應(yīng)要求：94 ℃預(yù)變性5 min，94 ℃變性30 s，45～46 ℃退火30 s，72 ℃延伸1 min，35個循環(huán)；最后72 ℃延伸10 min；將擴(kuò)增好的PCR產(chǎn)物放于4 ℃冰箱保存?zhèn)溆肹24]。

1.3? PCR擴(kuò)增產(chǎn)物的純化與測序

將擴(kuò)增好的PCR產(chǎn)物取4～5 μL用1%的瓊脂糖凝膠電泳檢測[24]，將擴(kuò)增成功的產(chǎn)物和相應(yīng)的引物送至生工生物工程（上海）股份有限公司（成都分公司）進(jìn)行分離純化和雙向測序。

1.4? 數(shù)據(jù)分析

1）將獲得的DNA雙向測序數(shù)據(jù)根據(jù)峰圖人工校對測序結(jié)果，用DNAMAN 軟件拼接、比對和剪切[25]，生成對應(yīng)的線粒體COI、Cytb全序列文件后，通過在NCBI上執(zhí)行BLAST相似性檢索，以保證得到的序列為目的基因片段。

2）用Bioedit對DNA序列進(jìn)行比對、拼接和校正[26]。

3）用DNAsp生成單倍型，并計算單倍型多樣性（haplotype diversity，Hd）、核苷酸多樣性（nucleotide diversity，Pi）及核苷酸平均差異數(shù)（Average number of nucleotide differences， K）[27]。

4）對線粒體DNACytb和COI這兩個目的片段基因及串聯(lián)成的DNA片段分別構(gòu)建以下兩種系統(tǒng)發(fā)育樹：用PAUP4.0的鄰近結(jié)合法構(gòu)建NJ（Neighbor joining）樹[28]；用MrBayes 3.2的貝葉斯（Bayesian）推斷法構(gòu)建貝葉斯樹[29]。

2? 結(jié)果與分析

2.1? 堿基組成

三個采樣點(diǎn)的峨眉黑雞基因片段的堿基組成測試結(jié)果為：COI序列片段中A+T的百分比高于G+C，堿基G含量最低（見表3）；Cytb序列片段A+T含量較高，堿基C含量最高（見表4）；串聯(lián)片段中A+T的含量高于G+C的含量（見表5）。

2.2? 遺傳多樣性

2.2.1? 遺傳多樣性指數(shù)

COI基因有效測序長度為846 bp，其中包含549個變異位點(diǎn)、397個缺失位點(diǎn)、1個簡約信息位點(diǎn)，兩堿基變異數(shù)為1個，三堿基變異數(shù)為1個。共定義了12種單倍型。

Cytb基因有效測序長度為1 005 bp，其中包含554個變異位點(diǎn)、451個缺失位點(diǎn)、8個簡約信息位點(diǎn)，兩堿基變異數(shù)為5個，三堿基變異數(shù)為3個。共定義了6種單倍型。

按照文獻(xiàn)[30]提出的標(biāo)準(zhǔn)，單倍型多樣性以0.5為臨界值，核苷酸多樣性以0.005為臨界值，二者的值越大，群體的多樣性程度越高。本測試結(jié)果（見表6～表8）顯示：采樣于峨眉山市世海黑雞原種雞場的種群（E1）在三個種群中遺傳多樣性最低；采樣于峨眉山市遠(yuǎn)大林下放養(yǎng)雞場的種群（E2）的核苷酸多樣性及核苷酸平均差異數(shù)最高；采樣于峨眉山市鴻飛養(yǎng)雞場的種群（E3）的單倍型多樣性較高，核苷酸多樣性較低。

2.2.2? 遺傳距離

測試結(jié)果顯示，峨眉黑雞三個地區(qū)雞種間基因的遺傳距離（DA）如下：1）Cytb基因片段， E1和E2兩個種群的遺傳距離為0.002 9，E1和E3兩個種群的遺傳距離為0.002 7，E2和E3兩個種群的遺傳距離為0.002 8。2）COI基因片段，E1和E2兩個種群的遺傳距離為0.000 8，E1和E3兩個種群的遺傳距離為0.000 4，E2和E3兩個種群的遺傳距離為0.000 5。

2.2.3? 系統(tǒng)發(fā)育樹

基于本研究測得的21條峨眉黑雞序列Dnasp 5分析得到3個采樣地點(diǎn)的峨眉黑雞COI基因序列共有11種單倍型，Cytb基因序列共有5種單倍型；采用鄰接法（NJ）和貝葉斯推論法（BI）構(gòu)建單倍型系統(tǒng)發(fā)育樹。

COI：1）NJ樹。單倍型H2、H3、H5聚在一起形成姊妹分支，H1、H4、H6、H7、H8、H9、H10、H11單獨(dú)形成一支（見圖1）。2）BI樹。所有單倍型均單獨(dú)形成分支（見圖2）。

Cytb：1）NJ樹。單倍型H2和H4聚在一起形成姊妹分支，H3和H5聚在一起形成姊妹分支，H1單獨(dú)形成分支（見圖3）。2）BI樹。H3和H5聚在一起形成姊妹分支，H1、H2、H4均單獨(dú)形成分支（見圖4）。

3? 結(jié)論與討論

3.1? 結(jié)論

群體內(nèi)反饋線粒體DNA變異水平程度的兩個關(guān)鍵要素是單倍型多樣性（Hd）及核苷酸多樣性（Pi），兩者的值越大，表明群體的遺傳多樣性越豐富，群體的多樣性程度也越高[31]，群體可供選擇的育種素材就越多，對將來變化多端的自然環(huán)境的適應(yīng)潛力也越大[32]。Cytb序列的多樣性（Hd=[0.667，0.905]，Pi=[0.002 35，0.005 14]）明顯高于COI序列的多樣性（Hd=[0.286，0.524]，Pi=[0.000 71，0. 0.001 05]），可能是由于COI基因相對保守，而Cytb基因突變率較高。對兩個基因片段及串聯(lián)序列多樣性的計算結(jié)果顯示，采樣于峨眉山市世海黑雞原種雞場的種群在三個種群中遺傳多樣性最低；采樣于峨眉山市遠(yuǎn)大林下放養(yǎng)雞場的種群的核苷酸多樣性及平均核苷酸差異數(shù)最高；采樣于峨眉山市鴻飛養(yǎng)雞場的種群單倍型多樣性較高，核苷酸多樣性較低。本研究結(jié)果表明，以原雞（Gallus）為外群，兩種分析方法得到的系統(tǒng)發(fā)育樹的結(jié)構(gòu)基本一致；來自不同地區(qū)的三個種群的單倍型未按照地理位置各自聚類，而是相互混雜在一起，進(jìn)化關(guān)系不明顯。

COI基因片段的遺傳多樣性分析顯示出單倍型多樣性（Hd）在0.286～0.524，核苷酸多樣性（Pi）在0.000 71～0.001 05，核苷酸平均差異數(shù)（K）在0.571 43～0.857 14?？傮w而言，在該基因片段序列中三個種群單倍型多樣性及核苷酸多樣性數(shù)值較低，本片段的三個種群遺傳多樣性數(shù)值均較低（Pi＜0.100）。Cytb基因片段的遺傳多樣性分析顯示出單倍型多樣性（Hd）在0.667～0.905，核苷酸多樣性（Pi）在0.002 35～0.005 14，核苷酸平均差異數(shù)（K）在1.809 52～4.285 71。采樣于峨眉山市鴻飛養(yǎng)雞場的種群單倍型多樣性較高（Hd=0.905＞0.900），采樣于峨眉山市遠(yuǎn)大林下放養(yǎng)雞場的種群序列核苷酸平均配對差異最高（K=4.285 71）。Cytb和COI基因串聯(lián)片段的遺傳多樣性分析顯示出單倍型多樣性（Hd）在0.810～0.905，核苷酸多樣性（Pi）在0.001 14～0.005 14，核苷酸平均差異數(shù)（K）在1.809 52～4.285 71。采樣于峨眉山市鴻飛養(yǎng)雞場的種群單倍型多樣性較高（Hd=0.905＞0.900），采樣于峨眉山市遠(yuǎn)大林下放養(yǎng)雞場的種群序列核苷酸平均配對差異最高（K=5.047 62）。

遺傳距離是分析2個不同種群或親緣關(guān)系較近的物種之間遺傳差異的一種方法，一般是利用2個種群不同位點(diǎn)的等位基因頻率數(shù)據(jù)來計算的[33]。本研究結(jié)果表明，無論是基于Cytb基因還是COI基因，采樣于峨眉山市鴻飛養(yǎng)雞場與峨眉山市世海黑雞原種雞場的種群遺傳距離都最?。ㄟz傳距離分別為0.002 7和0.000 4），即親緣關(guān)系最近；采樣于峨眉山市世海黑雞原種雞場與峨眉山市遠(yuǎn)大林下放養(yǎng)雞場的種群遺傳距離最大（遺傳距離分別為0.000 4和0.000 8），即親緣關(guān)系最遠(yuǎn)；但總體親緣關(guān)系都較近。

3.2? 討論

研究結(jié)果表明，峨眉黑雞的遺傳多樣性較低，三個地區(qū)間種群親緣關(guān)系都較近，推測可能產(chǎn)生該結(jié)果的原因是由于目前進(jìn)行保種的峨眉黑雞群體是來自2002年峨眉種雞的少量育種群體，這些雞種大部分來自當(dāng)?shù)氐霓r(nóng)戶，與其他外來雞種的基因交流可能性較小，所以種群親緣關(guān)系較近，其變異程度較低；除此之外，由于基礎(chǔ)種群數(shù)目相對較少，并且大部分在2002年之前缺乏系統(tǒng)性的保種措施，小群體中的遺傳漂變也可能會造成整體遺傳多樣性的缺失，發(fā)生在封閉群體中有效群體數(shù)量的下降也會提高基因漂變的速率，進(jìn)而減少小群體的遺傳多樣性[34]。因此，未來為確保峨眉黑雞的可持續(xù)性開發(fā)利用，需加強(qiáng)對其種質(zhì)資源的保護(hù)，減少這一品種的優(yōu)良性狀丟失。

參考文獻(xiàn)：

[1]? 藍(lán)太剛.生物多樣性保護(hù)的現(xiàn)狀和前景[J].林業(yè)科技開發(fā)，1993（3）：27-29.

[2]? 胡志昂，王洪新.遺傳多樣性的定義、研究新進(jìn)展和新概念[C]//生物多樣性與人類未來——第二屆全國生物多樣性保護(hù)與持續(xù)利用研討會論文集，1996：32-36.

[3]? 傅建軍，李家樂，沈玉幫，等.草魚野生群體遺傳變異的微衛(wèi)星分析[J].遺傳，2013，35（2）：192-201.

[4]? Groeneveld L F， Lenstra J A， Eding H， et al.Genetic diversity in farm animals[J]. Animal Genetics， 2010， 41（1）： 6-31.

[5]? 李凱年，逯德山.我國從國外引進(jìn)畜禽情況隱患與建議——對近幾年來我國畜牧業(yè)畜禽“引種熱”的冷思考[J].畜牧獸醫(yī)科技信息，2007（2）：4-7.

[6]? 桂小筍.洋豬肉猛敲本土生豬業(yè)警鐘[J].農(nóng)產(chǎn)品市場，2012（4）：36-37.

[7]? 王立賢.中國本土豬種選育刻不容緩[J].中國畜牧業(yè)，2014（7）：30.

[8]? 李鷹.中國畜禽育種，敢問路在何方——合理利用本國遺傳資源，培育具有中國特色的畜禽新品種[J].中國畜禽種業(yè)，2006（1）：4-6.

[9]? 陳國宏，李碧春，陳寬維.我國家禽品種資源的保存與利用[J].當(dāng)代畜牧;1996（1）：4-5，10.

[10]? 王奎.優(yōu)良地方品種四川山地烏骨雞[J].四川畜牧獸醫(yī)，2000（3）：36.

[11]? 陳昌元.體大肉嫩的峨眉黑雞[J].四川農(nóng)業(yè)科技，1985（4）：31.

[12]? 李苗，李地艷，陳彬龍，等.全基因組重測序鑒定峨眉黑雞的人工選擇模式[C]//第七屆中國畜牧科技論壇論文集，2016：315-316.

[13]? 朱文奇，李慧芳，宋衛(wèi)濤，等.家禽mtDNA遺傳多樣性和分子進(jìn)化研究進(jìn)展[J].家禽科學(xué)，2009（1）：40-42.

[14]? 徐文娟，朱文奇，束婧婷，等.我國主要烏骨雞品種遺傳多樣性和系統(tǒng)進(jìn)化研究[J].中國畜牧雜志，2014，50（23）：10-14.

[15]? 翁茁先，黃佳瓊，張仕豪，等.利用線粒體COI基因揭示中國烏骨雞遺傳多樣性和群體遺傳結(jié)構(gòu)[J].生物多樣性，2019，27（6）：10.

[16]? 吳瓊，寧浩然，宋超，等.基于線粒體細(xì)胞色素b基因序列的家養(yǎng)雉雞群體遺傳多樣性分析[J].中國家禽，2018，40（7）：52-54.

[17]? 黃勝海，賈曉旭，陸俊賢，等.基于線粒體標(biāo)記的淮南麻黃雞遺傳結(jié)構(gòu)及其遺傳多樣性分析[J].安徽農(nóng)業(yè)大學(xué)學(xué)報，2018，45（3）：400-403.

[18]? 王珍珍，盧立志，黃玲玲，等.新昌宮廷黃雞遺傳多樣性及遺傳結(jié)構(gòu)分析[J].中國畜牧雜志，2019，55（11）：5.

[19]? 陳艷榮，王繼文.鳥類mtDNA中cytb基因及在系統(tǒng)進(jìn)化研究中的應(yīng)用[J].黃山學(xué)院學(xué)報，2004（6）：88-90.

[20]? Ren G， Ma H， Ma C， et al. Genetic diversity and population structure of Portunus sanguinolentus （Herbst， 1783） revealed by mtDNA COI sequences[J]. Mitochondrial DNA Part A， 2016， 28（5）： 1-7.

[21]? Hebert P， Cywinska A， Ball S L， et al. Biological identification through DNA barcodes[J]. Proceedings of the Royal Society B， 2003， 270（1512）： 313-321.

[22]? 張艷，王金鴻，呂亞莉，等.六種試劑盒提取血液基因組DNA效果比較[J].中國草食動物科學(xué)，2014（2）：9-12.

[23]? Michael D S. Primers for a PCR-based approach to mitochondrial genome sequencing in birds and other vertebrates[J].Molecular Phylogenetics and Evolution， 1999， 12（2）： 105-114.

[24]? 趙文海.四川盆地的陸地島嶼效應(yīng)-基于樹麻雀的形態(tài)變異與遺傳結(jié)構(gòu)[D].烏魯木齊：新疆農(nóng)業(yè)大學(xué)，2018.

[25]? 王小國.生物學(xué)軟件在核酸序列比對與系統(tǒng)發(fā)育分析中的應(yīng)用[J].現(xiàn)代農(nóng)業(yè)科技，2015（12）：347-348.

[26] Hall T A. BioEdit： A User-Friendly Biological Sequence Alignment Editor and Analysis Program for Windows 95/98/NT[J]. Nuclc Acids Symposium Series， 1999， 41（41）： 95-98.

[27]? Julio R， Albert F M， Sánchez-DelBarrio Juan Carlos， et al. DnaSP 6： DNA Sequence Polymorphism Analysis of Large Data Sets[J]. Molecular Biology & Evolution， 2017（12）：14.

[28]? Jin-Chan W， Long-Xian Z， Fei L， et al. Phylogenetic Analysis of the Cryptosporidia Isolates Based on CPN60 Gene Sequence[J]. Chinese Journal of Animal & Veterinary Sciences， 2007， 38（11）： 1230-1234.

[29]? Ronquist F， Huelsenbeck J P. MrBayes 3： Bayesian phylogenetic inference under mixed models[J]. Bioinformatics， 2003， 19（12）： 1572-1574.

[30]? Grant WAS， Bowen BW. Shallow population histories in deep evolutionary lineages of marine fishes： Insights from sardines and anchovies and lessons for conservation[J]. Journal of Heredity， 1998， 89（5）： 415–426.

[31] 賈曉旭，唐修君，陸俊賢，等.基于線粒體DNAD-loop區(qū)全序列分析藏雞遺傳多樣性及其起源進(jìn)化關(guān)系的研究[J].東北農(nóng)業(yè)大學(xué)學(xué)報，2015（2）：32-36.

[32]? 陳寬維.我國畜禽遺傳資源保護(hù)理論與方法[J].中國禽業(yè)導(dǎo)刊，2009，26（23）：26-27.

[33]? Nei M. Genetic distance[J]. Brenners Encyclopedia of Genetics， 2013， 11（138）： 248-250.

[34]? 魯云風(fēng)，王茜，黃文波，等.畜禽保種群體遺傳多樣性的模擬[J].中國農(nóng)業(yè)科學(xué)，2012，45（18）：3849-3858.

（責(zé)任編輯：丁志祥）