峨眉黑鸡遗传多样性及群体遗传结构分析

袁霞 刘方庆 文陇英 徐婧 廖光祥 王强胜 王湘

（西南山地濒危鸟类保护四川省高等院校重点实验室/乐山师范学院，四川乐山 614000）

摘 要 峨眉黑鸡（Gallus gallus）为原鸡属红原鸡种的一个地方亚种，肉蛋兼用型品种，不仅具有丰富的营养价值，还兼有药用价值和观赏价值，是四川峨眉山地区养殖数量较多的特有的黑鸡优秀品种。为了给峨眉黑鸡种质资源的科学保存和可持续性利用提供一定的理论依据，采集峨眉山地区三个地点的峨眉黑鸡各7只，以线粒体基因Cytb和COI作为分子标记，分析种群的遗传多样性特征与种群遗传结构。结果表明，不同采样点的峨眉黑鸡遗传多样性存在差异，峨眉山市鸿飞养鸡场的种群具有较高的单倍型多样性（Hd=0.905），峨眉山市远大林下放养鸡场的种群序列核苷酸平均配对差异最高（K=5.047 62）。通过分析认为，峨眉黑鸡的遗传多样性并不丰富，建议加强对峨眉黑鸡种质资源的保护，保持这一品种的特征性不丢失。

关键词 峨眉黑鸡（Gallus gallus）；遗传多样性；系统进化；线粒体基因；四川省峨眉山市

中图分类号：S831；Q311+.3 文献标志码：A DOI：10.19415/j.cnki.1673-890x.2022.19.001

生物多樣性由遗传多样性、物种多样性和生态系统多样性三个部分组成，既包括生物个体及其携带的遗传信息，也包括它们与生态环境所组成的生态系统，以及各组分间的相互联系[1]。遗传多样性可以作为进一步分析生物个体在遗传进化方面的潜能和对环境的适应能力等各方面特征的重要基础数据之一[2]。一个地区的物种遗传多样性越丰富，该物种对环境的适应能力越强，在生存和进化方面也具有更大的优势[3]。

地方家禽种质遗传资源是畜禽遗传多样性的重要组成部分，是对地方禽类进行品种改良和可持续开发利用的重要依据[4]。随着国内经济的快速发展，引进外来种类（品种）逐渐增多；这些外来优良资源在丰富当地畜禽品种资源的同时，也会对当地原有畜禽品种资源多样性的保护造成冲击和破坏[5]。例如我国本土猪在大量引进进口白猪的情况下，面临着逐年减少甚至灭绝的情况[6-7]。地方品种虽然在生产性能方面比不上引进品种，但其优良的肉质品质和独特的风味远远优于引进品种，尤其是我国地方禽畜品种资源在遗传育种方面还可能具有较高的研究利用价值[8]。通过研究家禽的遗传结构多样性和家禽群体遗传距离可以初步查明家禽各种类型间的共同起源历史及亲缘关系，对合理地使用资源及有效保存家禽地方品种资源、建立优良基因库有着重要的参考意义，是家禽遗传资源管理中的重要环节[9]。

峨眉黑鸡（Gallus gallus）为原鸡属红原鸡种的一个地方亚种，是一种肉蛋兼用型品种，不仅具有丰富的营养价值，还兼有药用价值和观赏价值，是四川峨眉山地区养殖数量较多的黑鸡优秀品种；属于四川山地乌鸡的两大品系之一[10]。峨眉黑鸡主要分布在峨眉山市的龙池、大为、龙门，乐山市沙湾区的范店、轸溪及峨边县的毛坪、共和等乡镇；主产区为峨眉山西南大渡河沿岸的三县交界处[11]。

峨眉黑鸡全基因组重测序实现了遗传进化分析及重要性状候选基因的预测[12]，但针对其遗传多样性及品种形成的历史研究内容还较少。目前，以mtDNA序列为分子标记运用DNA直接测序的方法，研究家禽遗传多样性[13-18]和系统分化[19-22]及起源的研究报道越来越多。本研究通过以峨眉山特有鸡种峨眉黑鸡作为研究对象，以线粒体DNA细胞色素c氧化酶亚基I（cytochromecoxidase I， COI）和细胞色素b（cytochromeb）两个基因为分子标记，通过采集三个地点的峨眉黑鸡群体样本，分析种群遗传多样性与遗传结构，构建系统发育树，为峨眉黑鸡种质资源的科学保存和可持续性利用提供一定的理论依据。

1  材料与方法

1.1  样品采集及DNA提取

1.1.1  样品采集

21只峨眉黑鸡分别来自峨眉山市三个的鸡场：峨眉山市世海黑鸡原种鸡场（E1），峨眉山市远大林下放养鸡场（E2），峨眉山市鸿飞养鸡场（E3）。样品编号见表1。所有个体用含有20 μL抗凝血剂（EDTA和0.75%生理盐水按3%～4%比例配制而成）的抗凝管进行翅下静脉采血，采样体积为0.3 mL，放置于-20 ℃冰箱保存。

1.1.2  DNA提取

使用天根组织/血液试剂盒（北京天根）提取基因组上的DNA序列[22]。

1.2  引物信息及PCR扩增

参照文献[23]选取合适的引物（见表2）。线粒体DNA 中的Cytb基因扩增引物为Cytb-H16064、Cytb-L14731，线粒体DNA中的COI基因扩增引物为COI-L6615、COI-H7956；其中L表示轻链，H表示重链；序列的顺序均为5′—3′方向。引物由生工生物工程（上海）股份有限公司合成。

以提取的基因组DNA作为模板，进行PCR扩增。PCR反应体系为25 μL，包括模板DNA 2 μL、ddH2O 8.5 μL、DNA聚合酶（Taq酶）12.5 μL（含dNTP、Mg2+、Buffer缓冲液）、上下游引物各1 μL，反应体系在Easycycler Gradient 96型与SelectCycler Ⅱ型PCR仪上进行反应[24]。

PCR反应要求：94 ℃预变性5 min，94 ℃变性30 s，45～46 ℃退火30 s，72 ℃延伸1 min，35个循环；最后72 ℃延伸10 min；将扩增好的PCR产物放于4 ℃冰箱保存备用[24]。

1.3  PCR扩增产物的纯化与测序

将扩增好的PCR产物取4～5 μL用1%的琼脂糖凝胶电泳检测[24]，将扩增成功的产物和相应的引物送至生工生物工程（上海）股份有限公司（成都分公司）进行分离纯化和双向测序。

1.4  数据分析

1）将获得的DNA双向测序數据根据峰图人工校对测序结果，用DNAMAN 软件拼接、比对和剪切[25]，生成对应的线粒体COI、Cytb全序列文件后，通过在NCBI上执行BLAST相似性检索，以保证得到的序列为目的基因片段。

2）用Bioedit对DNA序列进行比对、拼接和校正[26]。

3）用DNAsp生成单倍型，并计算单倍型多样性（haplotype diversity，Hd）、核苷酸多样性（nucleotide diversity，Pi）及核苷酸平均差异数（Average number of nucleotide differences， K）[27]。

4）对线粒体DNACytb和COI这两个目的片段基因及串联成的DNA片段分别构建以下两种系统发育树：用PAUP4.0的邻近结合法构建NJ（Neighbor joining）树[28]；用MrBayes 3.2的贝叶斯（Bayesian）推断法构建贝叶斯树[29]。

2  结果与分析

2.1  碱基组成

三个采样点的峨眉黑鸡基因片段的碱基组成测试结果为：COI序列片段中A+T的百分比高于G+C，碱基G含量最低（见表3）；Cytb序列片段A+T含量较高，碱基C含量最高（见表4）；串联片段中A+T的含量高于G+C的含量（见表5）。

2.2  遗传多样性

2.2.1  遗传多样性指数

COI基因有效测序长度为846 bp，其中包含549个变异位点、397个缺失位点、1个简约信息位点，两碱基变异数为1个，三碱基变异数为1个。共定义了12种单倍型。

Cytb基因有效测序长度为1 005 bp，其中包含554个变异位点、451个缺失位点、8个简约信息位点，两碱基变异数为5个，三碱基变异数为3个。共定义了6种单倍型。

按照文献[30]提出的标准，单倍型多样性以0.5为临界值，核苷酸多样性以0.005为临界值，二者的值越大，群体的多样性程度越高。本测试结果（见表6～表8）显示：采样于峨眉山市世海黑鸡原种鸡场的种群（E1）在三个种群中遗传多样性最低；采样于峨眉山市远大林下放养鸡场的种群（E2）的核苷酸多样性及核苷酸平均差异数最高；采样于峨眉山市鸿飞养鸡场的种群（E3）的单倍型多样性较高，核苷酸多样性较低。

2.2.2  遗传距离

测试结果显示，峨眉黑鸡三个地区鸡种间基因的遗传距离（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  系统发育树

基于本研究测得的21条峨眉黑鸡序列Dnasp 5分析得到3个采样地点的峨眉黑鸡COI基因序列共有11种单倍型，Cytb基因序列共有5种单倍型；采用邻接法（NJ）和贝叶斯推论法（BI）构建单倍型系统发育树。

COI：1）NJ树。单倍型H2、H3、H5聚在一起形成姊妹分支，H1、H4、H6、H7、H8、H9、H10、H11单独形成一支（见图1）。2）BI树。所有单倍型均单独形成分支（见图2）。

Cytb：1）NJ树。单倍型H2和H4聚在一起形成姊妹分支，H3和H5聚在一起形成姊妹分支，H1单独形成分支（见图3）。2）BI树。H3和H5聚在一起形成姊妹分支，H1、H2、H4均单独形成分支（见图4）。

3  结论与讨论

3.1  结论

群体内反馈线粒体DNA变异水平程度的两个关键要素是单倍型多样性（Hd）及核苷酸多样性（Pi），两者的值越大，表明群体的遗传多样性越丰富，群体的多样性程度也越高[31]，群体可供选择的育种素材就越多，对将来变化多端的自然环境的适应潜力也越大[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基因突变率较高。对两个基因片段及串联序列多样性的计算结果显示，采样于峨眉山市世海黑鸡原种鸡场的种群在三个种群中遗传多样性最低；采样于峨眉山市远大林下放养鸡场的种群的核苷酸多样性及平均核苷酸差异数最高；采样于峨眉山市鸿飞养鸡场的种群单倍型多样性较高，核苷酸多样性较低。本研究结果表明，以原鸡（Gallus）为外群，两种分析方法得到的系统发育树的结构基本一致；来自不同地区的三个种群的单倍型未按照地理位置各自聚类，而是相互混杂在一起，进化关系不明显。

COI基因片段的遗传多样性分析显示出单倍型多样性（Hd）在0.286～0.524，核苷酸多样性（Pi）在0.000 71～0.001 05，核苷酸平均差异数（K）在0.571 43～0.857 14。总体而言，在该基因片段序列中三个种群单倍型多样性及核苷酸多样性数值较低，本片段的三个种群遗传多样性数值均较低（Pi＜0.100）。Cytb基因片段的遗传多样性分析显示出单倍型多样性（Hd）在0.667～0.905，核苷酸多样性（Pi）在0.002 35～0.005 14，核苷酸平均差异数（K）在1.809 52～4.285 71。采样于峨眉山市鸿飞养鸡场的种群单倍型多样性较高（Hd=0.905＞0.900），采样于峨眉山市远大林下放养鸡场的种群序列核苷酸平均配对差异最高（K=4.285 71）。Cytb和COI基因串联片段的遗传多样性分析显示出单倍型多样性（Hd）在0.810～0.905，核苷酸多样性（Pi）在0.001 14～0.005 14，核苷酸平均差异数（K）在1.809 52～4.285 71。采样于峨眉山市鸿飞养鸡场的种群单倍型多样性较高（Hd=0.905＞0.900），采样于峨眉山市远大林下放养鸡场的种群序列核苷酸平均配对差异最高（K=5.047 62）。

遗传距离是分析2个不同种群或亲缘关系较近的物种之间遗传差异的一种方法，一般是利用2个种群不同位点的等位基因频率数据来计算的[33]。本研究结果表明，无论是基于Cytb基因还是COI基因，采样于峨眉山市鸿飞养鸡场与峨眉山市世海黑鸡原种鸡场的种群遗传距离都最小（遗传距离分别为0.002 7和0.000 4），即亲缘关系最近；采样于峨眉山市世海黑鸡原种鸡场与峨眉山市远大林下放养鸡场的种群遗传距离最大（遗传距离分别为0.000 4和0.000 8），即亲缘关系最远；但总体亲缘关系都较近。

3.2  讨论

研究结果表明，峨眉黑鸡的遗传多样性较低，三个地区间种群亲缘关系都较近，推测可能产生该结果的原因是由于目前进行保种的峨眉黑鸡群体是来自2002年峨眉种鸡的少量育种群体，这些鸡种大部分来自当地的农户，与其他外来鸡种的基因交流可能性较小，所以种群亲缘关系较近，其变异程度较低；除此之外，由于基础种群数目相对较少，并且大部分在2002年之前缺乏系统性的保种措施，小群体中的遗传漂变也可能会造成整体遗传多样性的缺失，发生在封闭群体中有效群体数量的下降也会提高基因漂变的速率，进而减少小群体的遗传多样性[34]。因此，未来为确保峨眉黑鸡的可持续性开发利用，需加强对其种质资源的保护，减少这一品种的优良性状丢失。

参考文献：

[1]  蓝太刚.生物多样性保护的现状和前景[J].林业科技开发，1993（3）：27-29.

[2]  胡志昂，王洪新.遗传多样性的定义、研究新进展和新概念[C]//生物多样性与人类未来——第二届全国生物多样性保护与持续利用研讨会论文集，1996：32-36.

[3]  傅建军，李家乐，沈玉帮，等.草鱼野生群体遗传变异的微卫星分析[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].畜牧兽医科技信息，2007（2）：4-7.

[6]  桂小笋.洋猪肉猛敲本土生猪业警钟[J].农产品市场，2012（4）：36-37.

[7]  王立贤.中国本土猪种选育刻不容缓[J].中国畜牧业，2014（7）：30.

[8]  李鹰.中国畜禽育种，敢问路在何方——合理利用本国遗传资源，培育具有中国特色的畜禽新品种[J].中国畜禽种业，2006（1）：4-6.

[9]  陈国宏，李碧春，陈宽维.我国家禽品种资源的保存与利用[J].当代畜牧;1996（1）：4-5，10.

[10]  王奎.优良地方品种四川山地乌骨鸡[J].四川畜牧兽医，2000（3）：36.

[11]  陈昌元.体大肉嫩的峨眉黑鸡[J].四川农业科技，1985（4）：31.

[12]  李苗，李地艳，陈彬龙，等.全基因组重测序鉴定峨眉黑鸡的人工选择模式[C]//第七届中国畜牧科技论坛论文集，2016：315-316.

[13]  朱文奇，李慧芳，宋卫涛，等.家禽mtDNA遗传多样性和分子进化研究进展[J].家禽科学，2009（1）：40-42.

[14]  徐文娟，朱文奇，束婧婷，等.我国主要乌骨鸡品种遗传多样性和系统进化研究[J].中国畜牧杂志，2014，50（23）：10-14.

[15]  翁茁先，黃佳琼，张仕豪，等.利用线粒体COI基因揭示中国乌骨鸡遗传多样性和群体遗传结构[J].生物多样性，2019，27（6）：10.

[16]  吴琼，宁浩然，宋超，等.基于线粒体细胞色素b基因序列的家养雉鸡群体遗传多样性分析[J].中国家禽，2018，40（7）：52-54.

[17]  黄胜海，贾晓旭，陆俊贤，等.基于线粒体标记的淮南麻黄鸡遗传结构及其遗传多样性分析[J].安徽农业大学学报，2018，45（3）：400-403.

[18]  王珍珍，卢立志，黄玲玲，等.新昌宫廷黄鸡遗传多样性及遗传结构分析[J].中国畜牧杂志，2019，55（11）：5.

[19]  陈艳荣，王继文.鸟类mtDNA中cytb基因及在系统进化研究中的应用[J].黄山学院学报，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].中国草食动物科学，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]  赵文海.四川盆地的陆地岛屿效应-基于树麻雀的形态变异与遗传结构[D].乌鲁木齐：新疆农业大学，2018.

[25]  王小国.生物学软件在核酸序列比对与系统发育分析中的应用[J].现代农业科技，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区全序列分析藏鸡遗传多样性及其起源进化关系的研究[J].东北农业大学学报，2015（2）：32-36.

[32]  陈宽维.我国畜禽遗传资源保护理论与方法[J].中国禽业导刊，2009，26（23）：26-27.

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

[34]  鲁云风，王茜，黄文波，等.畜禽保种群体遗传多样性的模拟[J].中国农业科学，2012，45（18）：3849-3858.

（责任编辑：丁志祥）