畜牧兽医学报  2022, Vol. 53 Issue (2): 370-380. DOI: 10.11843/j.issn.0366-6964.2022.02.005    PDF    
四川盆地6个地方鸡群体遗传关系研究
尤桂爽1, 易兴友2, 蒋美山2, 李中伟2, 谭和林2, 曹池2, 舒刚1, 田尧夫1, 朱庆1, 赵小玲1     
1. 四川农业大学动物科技学院, 成都 611130;
2. 重庆市开州区农业发展服务中心, 重庆 405499
摘要:随着交通便利与市场化杂交利用,我国地方鸡遗传资源面临血液混杂、品种来源和品种间的遗传关系不清晰的问题。为了给第三次全国畜禽遗传资源普查工作提供依据。本研究对位于中国大巴山脉和乌蒙山区的6个地方鸡群体(每个群体公、母各10只)共118个个体进行了简化基因组GBS测序,分析了核酸多态性(Pi)、观测杂合度(Ho)、期望杂合度(He)、多态信息含量(PIC)、观测等位基因数(Na)、有效等位基因数(Ne)、遗传分化(Fst)、遗传距离(DR)和基因流(Nm)等指标。研究发现,岩水鸡黑羽群体的遗传多样性较低(He=0.154;PIC=0.125),大宁河鸡遗传多样性最为丰富(He=0.197,PIC=0.162)。城口山地鸡和大宁河鸡之间的遗传分化最低(Fst=0.053),基因流最高(Nm=4.50),丰岩乌骨鸡和岩水鸡黑羽群体之间的遗传分化最高(Fst=0.183),对应的基因流最低(Nm=1.12)。系统发生树和主成分分析显示,岩水鸡白羽群体和黑羽群体聚为一类,城口山地鸡、大宁河鸡和旧院黑鸡聚为一类,丰岩乌骨鸡单独为一类。Structure分析表明,城口山地鸡与岩水鸡白羽群体血缘混杂。本研究探明了邻近区域地方品种的遗传关系,位于乌蒙山区的丰岩乌骨鸡与大巴山脉的旧院黑鸡、城口山地鸡、大宁河鸡、岩水鸡黑羽和白羽具有较远的遗传距离。本研究为地方品种的保护和开发利用提供了依据。
关键词地方鸡品种    遗传多样性    简化基因组    基因流    系统发生树    
Study on the Genetic Relationship of Six Local Chicken Populations from Sichuan Basin
YOU Guishuang1, YI Xingyou2, JIANG Meishan2, LI Zhongwei2, TAN Helin2, CAO Chi2, SHU Gang1, TIAN Yaofu1, ZHU Qing1, ZHAO Xiaoling1     
1. College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China;
2. Agricultural Development Service Center of Kaizhou District, Chongqing 405499, China
Abstract: With the convenience of transportation and market-oriented cross utilization, the genetic resources of local chickens in China are facing the problems of blood mixing and the unclear breed source and genetic relationship among breeds. The study aimed to provide a basis for the Third National Census of Livestock and Poultry Genetic Resources. A simplified genomic GBS sequencing was carried out on 118 individuals from 6 local chicken populations (10 males and 10 females in each population) located in Daba Mountains and Wumeng Mountains in China. The nucleotide polymorphism (Pi), observed heterozygosity (Ho), expected heterozygosity (He), polymorphic information content (PIC), number of observed alleles (Na), number of effective alleles (Ne), genetic differentiation (Fst), genetic distance (DR) and gene flow (Nm) of these populations were analyzed. It was found that the genetic diversity of black feather population of Yanshui chicken was low (He=0.154; PIC=0.125), and the genetic diversity of Daninghe chicken was the most abundant (He=0.197, PIC=0.162). The genetic differentiation between Chengkou mountain chicken and Daninghe chicken was the lowest (Fst =0.053), the gene flow was the highest (Nm=4.50), the genetic differentiation between Fengyan black bone chicken and black feather population of Yanshui chicken was the highest (Fst=0.183), and the corresponding gene flow was the lowest (Nm=1.12). Phylogenetic tree and principal component analysis showed that the white feather population and black feather population of Yanshui chicken were clustered into one class. Chengkou mountain chicken, Daninghe chicken and Jiuyuan black chicken were clustered into one class, while Fengyan black bone chicken was a separate class. Structure analysis showed that the blood of Chengkou mountain chicken and white feather population of Yanshui chicken was mixed. This study proved the genetic relationship of local breeds in adjacent areas. Fengyan black bone chicken located in Wumeng Mountain has a long genetic distance from Jiuyuan black chicken, Chengkou mountain chicken, Daninghe chicken, black feather and white feather populations of Yanshui chicken in Daba Mountain. This study provides a basis for the protection, development and utilization of local chicken breeds.
Key words: local chicken breeds    genetic diversity    simplified genome    gene flow    phylogenetic tree    

四川盆地由中国西部横断山脉、云贵高原、大巴山、武夷山环绕形成,盆地北部大巴山横直在万源宣汉北部,明月山、铜锣山、华釜山由北而南纵卧其间,自然分割形成盆地和丘陵地貌。海拔跨度极大的地势造就了丰富的生物多样性,形成了包括国家畜禽遗传资源旧院黑鸡、城口山地鸡、大宁河鸡等在内的家禽品种[1]。在大巴山脉与长江山峡的交界处,有1个新的家禽品种岩水鸡,其体态轻盈,羽色为纯黑或纯白,目前处于地方畜牧部门保护开发的计划中。南部地区武夷山和云贵高原的交汇点也有显著的海拔变化,分布着沐川乌骨鸡、丰岩乌骨鸡等家禽品种[2]

城口山地鸡地处重庆市城口县,为肉蛋兼用型品种,羽色纯黑,黑肤、乌肉,肉质细嫩,深受当地人民的喜爱。旧院黑鸡主产于四川省万源市旧院镇,是肉蛋兼用型品种,外观与城口山地鸡近似,羽色纯黑,黑肤、乌肉,肉质鲜美[3]。大宁河鸡产于重庆市巫溪县,属兼用型地方品种,公鸡红羽,母鸡黄麻羽。岩水鸡白羽和黑羽群体产于重庆开州区,北邻城口东邻巫溪,是出名的高寒山区,该品种环境适应性强,肉蛋兼用[4]。丰岩乌骨鸡羽毛、皮、骨、肉俱黑,产于四川泸州叙永县,已有300年以上历史,长期以来当地苗、汉农民就养殖乌骨鸡馈赠亲友,配以中草药滋补疗疾。

近年来随着经济的发展,交通的便利,出现了相邻地域品种混杂,毛色不纯,肤色变浅等问题,根据农业部的报告,一些地方鸡种的种群数量一直在迅速减少甚至濒临灭绝[5]。第三次全国畜禽遗传资源普查品种鉴定工作正在开展,为说明四川及其临近区域品种间遗传多样性,便于后续的品种保护和开发利用[6],本研究采用简化基因组开展了6个地方鸡群体的遗传研究。简化基因组测序技术(genotyping-by-sequencing,GBS)是在二代测序基础上发展起来的[7],通过酶切加标签的方法测序获得酶切位点附近基因序列信息,进而检测大量高准确性的SNP变异位点信息,可了解种质资源的遗传背景和系统演化,研究群体亲缘关系及种群结构[8]。Pértille等[9]用GBS法对462只鸡进行测序,鉴定出67 096个SNPs,全基因组覆盖率为4.66%。猪群测序检测到假定的SNPs,平均密度为每10 kb有0.33个SNPs[10]。此外,在一项牛研究中还鉴定出63 797个SNPs[11]。GBS是一种经济高效的方法,已成功地应用于群体遗传学研究[12-14]

1 材料与方法 1.1 试验材料

本研究使用了来自5个保种场的6个地方品种鸡,每个品种选取20个个体(10公、10母)共120只,150~250日龄,分别为城口山地鸡、大宁河鸡、旧院黑鸡、丰岩乌骨鸡、岩水鸡黑羽和岩水鸡白羽,采样信息见表 1。将样本进行翅下静脉采血,并保存于-20 ℃。使用标准的苯酚/氯仿方法从全血中提取基因组DNA。用分光光度计(德国埃普多夫公司)定量DNA浓度,并将样品稀释至最终浓度100 ng·μL-1[15]

表 1 地方鸡群体采样信息 Table 1 Information of local chicken populations
1.2 试验方法

1.2.1 简化基因组测序和SNP鉴定   DNA样品送往广州基迪奥生物科技有限公司利用Gallus gallus基因组(GenBank登录号:GCF_000002315.6)作为参考基因组,用GBS技术构建测序文库,采用限制性内切酶CTGCAG(Pst Ⅰ) + CATG(Nla Ⅲ)对DNA进行酶切,连接回收后对片段使用高保真酶进行PCR扩增,测定PCR产物浓度后将混好的文库上机测序(Illumina HiSeqTM 2500)。

高通量测序得到的原始图像数据文件经碱基识别转化为原始测序序列,过滤后使用比对软件BWA(0.7.12)采用mem算法后将过滤后的reads比对到参考基因组上,比对参数为-k 32 -M;比对完之后结果使用PICARD(版本1.129)软件进行标记[16]。使用变异检测软件GATK(版本3.4-46)进行群体SNP检测,然后使用ANNOVAR(版本2)对检测出的突变位点进行功能注释。

1.2.2 遗传数据分析   用Perl编程计算观测杂合度(observed heterozygosity,Ho)、期望杂合度(expected heterozygosity,He)、多态信息含量(polymorphic information content,PIC)、观测等位基因数(number of observed alleles,Na)、有效等位基因数(number of effective alleles,Ne)、核苷酸多样性(nucleotide diversity,Pi)、遗传分化系数(F-statistic,Fst)、遗传距离(genetic distance,DR)及基因流(gene flow,Nm)。利用软件GCTA(Version 1.25.0)进行主成分分析(principle component analysis,PCA),采用邻接法(neighbor-joining method)构建系统发生树,通过Treebest软件计算距离矩阵,系统发生树的可靠性通过Boot-strap法进行检验。使用Admixture软件(v1.3)[17]进行群体结构的推断,分别假设样品的分群数(K值)为1~9,然后进行聚类,根据交叉验证错误率来确定最优分群数。

2 结果 2.1 6个地方鸡群体的基因组SNPs鉴定

2.1.1 测序质量   120个地方鸡样本,剔除2个核酸质量低的样本,通过GBS测序,剩余118个样本测序质量较高(Q20≥97.58%、所有样品的平均测序深度为20×,平均覆盖度为6.48%)。经过筛选过滤后,118个样品共获得291 359个高质量的SNPs位点。Q30≥92.74%,GC分布正常,群体样本与参考基因组的平均比对率为99.33%。

2.1.2 突变筛选   去除等位基因数大于2、缺失数据大于50%、次要等位基因频率小于5%、杂合度大于80%的标记(SNPs/Indels),保留271 672个SNPs和19 687个Indels,用于群体遗传分析。

2.2 6个地方鸡群体的遗传多样性分析

6个群体的核酸多态性范围为6.71×10-5(岩水鸡黑羽)~8.58×10-5(大宁河鸡);观测杂合度范围为0.133(岩水白羽鸡)~0.168(城口山地鸡);期望杂合度范围为0.154(岩水黑羽鸡) ~0.197(大宁河鸡);多态信息含量范围为0.125(岩水黑羽鸡)~0.162(大宁河鸡),属低度多态性(PIC < 0.25);观测等位基因数范围为1.534(岩水鸡黑羽)~1.750(城口山地鸡);有效等位基因数范围为1.256(岩水鸡黑羽)~1.318(大宁河鸡)(表 2)。整体上大宁河鸡和城口山地鸡的遗传多样性较高,岩水鸡黑羽和岩水鸡白羽的遗传多样性较低。

表 2 6个鸡群体遗传多样性参数 Table 2 Genetic diversity parameters of 6 chicken populations
2.3 6个地方鸡群体的遗传分化系数和遗传距离

6个群体间两两遗传分化系数(Fst)范围为0.053~0.183,遗传距离(DR)为0.087~0.119,基因流(Nm)为1.12~4.50。岩水鸡黑羽与大宁河鸡、旧院黑鸡、丰岩乌骨鸡,岩水鸡白羽与丰岩乌骨鸡的遗传分化较大(0.15≤Fst≤0.25),其余品种之间的遗传分化中等(0.05≤Fst≤0.15)。群体间遗传分化最大的是丰岩乌骨鸡和岩水鸡黑羽群体(Fst=0.183),最小的是城口山地鸡和大宁河鸡(Fst=0.053)。群体间遗传距离最近的是岩水鸡黑羽和白羽群体(DR=0.087),遗传距离最远的是大宁河鸡和岩水鸡白羽群体(DR=0.119)。群体间基因流动最大的是大宁河鸡和城口山地鸡(Nm=4.05),流动最小的是岩水鸡黑羽和丰岩乌骨鸡(Nm=1.12)(表 3表 4)。

表 3 6个群体间的遗传分化系数(Fst)和遗传距离(DR) Table 3 F-statistic and genetic distances among 6 populations
表 4 6个群体的基因流 Table 4 Gene flow among 6 populations
2.4 6个地方鸡群体遗传结构分析

2.4.1 主成分分析   对获得的SNPs标记进行主成分分析,结果表明3个主成分变异贡献率分别为7.31%、4.31%和4.04%(图 1)。由图 1可以看出,城口山地鸡、大宁河鸡、旧院黑鸡形成一个紧密的类群,丰岩乌骨鸡、岩水鸡白羽和黑羽相对独立,并且丰岩乌骨鸡和岩水鸡黑羽群体相对集中,岩水鸡白羽群体相对分散。

图 1 地方鸡6个群体的主成分分析 Fig. 1 Principal components analysis of 6 local chicken populations

2.4.2 系统发生树分析   根据地方鸡6个群体的系统发生树可以明显看出(图 2),城口山地鸡、大宁河鸡、旧院黑鸡和丰岩乌骨鸡各自聚为一类,岩水鸡白羽和黑羽聚为一类。

图 2 地方鸡6个群体的系统发生树 Fig. 2 Phylogenetic tree of 6 local chicken populations
2.5 6个地方鸡群体遗传结构分析

预设群体分群数从2到9,通过Admixture软件分析地方鸡的群体结构,根据交叉验证错误率的谷值确定最优分群数(K)。如图 3所示,当K=5时交叉验证错误率最低,说明本研究分析的样品最大概率来自于5个原始的祖先。图 4展示了6个群体的遗传结构。当K=2和3时,岩水鸡黑羽和白羽明显具有和其它品种不同的来源。在K=4时,丰岩乌骨鸡与岩水鸡黑羽遗传结构差异大,而城口山地鸡、大宁河鸡、旧院黑鸡遗传结构近似。在K=5时,大宁河鸡、旧院黑鸡、丰岩乌骨鸡和岩水鸡黑羽遗传差异明显,城口山地鸡和岩水鸡白羽呈现出部分血液混杂,K=6~9呈现出更详细的品种遗传信息。

图 3 交叉验证错误率 Fig. 3 Cross validation error rate
图 4 6个地方鸡品种的结构聚类分析 Fig. 4 Structure clustering of 6 local chicken populations
3 讨论

对品种内多样性和物种种群结构的了解是改进品种选择计划、了解群体环境适应性和政府实施种群保护计划的基础[18-19]。虽然大多数群体遗传结构研究都针对世界性的商用品种,但人们对地方品种的遗传学关系越来越感兴趣,因为它们可以从根本上理清品种的内在关系[20]

本试验利用简化基因组研究了四川盆地周边山脉6个地方鸡种的基因组多态性。期望杂合度(He)和多态信息含量(PIC)是衡量标记多态性较好的指标[21-22]。6个地方鸡品种的He范围为0.154~ 0.197,平均期望杂合度和Cendron等[23]用全基因组测序分析意大利地方品种鸡的结果相似。6个地方鸡品种的PIC范围为0.125~ 0.162,为低度多态性(PIC < 0.25低度多态性;0.25 < PIC < 0.5中度多态性;PIC>0.5高度多态性)[24]。本研究中PIC含量结果低于微卫星的结果[25],这是由于标记类型之间的差异。由于SNP的双等位特性,简化基因组的PIC在0~0.5之间,而对于多等位的微卫星标记,PIC在0~1.0之间,由于前者的PIC较低,因此需要比微卫星更多的SNPs才能获得相同的排斥力[26]。Abdalhag等[27]也检测到类似的低PIC值。在这6个地方品种中,岩水鸡黑羽和岩水鸡白羽群体的杂合度较低,6个种群遗传多态性最高的是大宁河鸡,最低的是岩水鸡白羽和黑羽群体,其原因是大宁河鸡的群体数量较大且分布的范围较广,而岩水鸡白羽和黑羽群体可能由于保种群体小,存在遗传漂变,并局限在一个较为封闭的高山环境中,与其它群体之间基因交流的机会少。从基因流也可以看出,这两个群体同其它4个群体的基因流动值最小,因此该品种具有一定的遗传独特性,这突显了保护这些品种生物多样性的难度,因此必须特别注意这些品种的保护,可以采用各家系等量留种的保种方式,尽量保持现有等位基因,防止近交衰退[28]。大宁河鸡的遗传多样性较为丰富,其保种效果较好。

为了解这6个地方品种之间的遗传关系和群体结构,本研究计算了群体的遗传分化指数(Fst)和遗传距离(DR)的成对估计,并进行了主成分、系统发生树和Structure分析。Fst揭示了种群间的遗传分化程度,本试验中城口山地鸡和大宁河鸡间分化中等,旧院黑鸡和丰岩乌骨鸡间分化较大[29-30]。前人的研究表明,印度[31]、中国[32]和非洲本土鸡群[33]Fst值较高(0.15到0.26),而Rashid等[34]报告了孟加拉国5个地方鸡种的Fst在0.084~0.154之间。Hassen等[35]报告了埃塞俄比亚西北部7个土生鸡种群的地理距离从0.073到0.13。Seo等[36]发现5个韩国本土鸡之间的成对遗传距离为0.083~ 0.171。本研究结果与上述3个研究结果相似,遗传距离值较小。此外,Yamamoto[37]的研究报告了孟加拉土鸡不同群体之间的遗传分化程度较低。总体而言,四川盆地地方鸡群体间的遗传分化程度较低,这可能是由于它们之间以及构建群体的共同祖先之间的基因流动所致。本研究中Fst和基因流都显示城口山地鸡和大宁河鸡群体间的分化程度最小,丰岩乌骨鸡和岩水鸡黑羽群体间的分化最大。系统发生树显示,岩水鸡白羽和岩水鸡黑羽聚为一类,与基因流结果一致。利用Structure程序对中国地方6个鸡群的遗传结构进行了研究,当K=5时,交叉验证错误率最低,说明本研究分析的样品应该是来自于5个原始的祖先。结构数据还表明,在每个群体中有一定比例的遗传混合来自其他品种,邻近地区或种群可能存在一定的基因流动。同时结构数据表明,丰岩乌骨鸡和旧院黑鸡被明显分离开,这与主成分分析和系统发生树的结果一致。前人研究认为,两个品种间的遗传相似系数为0.795,二者亲缘关系较近[38]。检测方法不同可能导致了结果差异[39]。本研究结果还表明,岩水鸡白羽处于系统发生树的最远分支,与本研究其他地方鸡种相比存在明显的遗传差异。这在主成分、Fst和基于模型的聚类结果中均有体现。这个聚类结果与各品种在四川盆地的分布区域基本吻合, 说明四川地方鸡品种的形成与地理环境存在很大的关系[40]

4 结论

本研究应用简化基因组对四川盆地的6个地方鸡品种进行分析,结果发现岩水鸡黑羽群体的遗传多样性较低,大宁河鸡遗传多样性最为丰富。城口山地鸡和大宁河鸡之间的遗传分化最低,基因流最高,丰岩乌骨鸡和岩水鸡黑羽群体之间的遗传分化最高,对应的基因流最低。岩水鸡白羽群体和黑羽群体聚为一类,城口山地鸡、大宁河鸡和旧院黑鸡聚为一类,丰岩乌骨鸡单独为一类。位于乌蒙山区的丰岩乌骨鸡与大巴山脉的旧院黑鸡、城口山地鸡、大宁河鸡、岩水鸡黑羽和白羽具有较远的遗传距离。

参考文献
[1]
国家畜禽遗传资源委员会. 中国畜禽遗传资源志[M]. 北京: 中国农业出版社, 2011.
China National Commission of Animal Genetic Resources. Animal genetic resources in China[M]. Beijing: China Agriculture Press, 2011. (in Chinese)
[2]
汤青萍, 陈宽维, 李慧芳, 等. 应用微卫星标记对12个中国地方乌骨鸡品种遗传多样性的研究[J]. 畜牧兽医学报, 2005, 36(8): 755-760.
TANG Q P, CHEN K W, LI H F, et al. Analysis of the genetic diversity of 12 Chinese Indigenous Black-bone chicken breeds using microsatellite marker[J]. Acta Veterinaria et Zootechnica Sinica, 2005, 36(8): 755-760. DOI:10.3321/j.issn:0366-6964.2005.08.002 (in Chinese)
[3]
刘嘉, 李富贵, 桂涛涛, 等. 旧院黑鸡及其三系配套组合的外貌特征及肉用性能分析[J]. 西北农业学报, 2019, 28(2): 161-168.
LIU J, LI F G, GUI T T, et al. Appearance characteristics and meat performances analysis of Jiuyuan black chicken and three-line hybrid combinations[J]. Acta Agriculturae Boreali-occidentalis Sinica, 2019, 28(2): 161-168. (in Chinese)
[4]
易兴友, 曹池, 陈建华, 等. 开州岩水鸡品种资源调查[J]. 中国畜牧业, 2021(14): 53-55.
YI X Y, CAO C, CHEN J H, et al. Investigation on breed resources of Kaizhou Yanshui chicken[J]. China Animal Industry, 2021(14): 53-55. DOI:10.3969/j.issn.2095-2473.2021.14.026 (in Chinese)
[5]
《中国畜禽遗传资源状况》编委会. 中国畜禽遗传资源状况[M]. 北京: 中国农业出版社, 2004.
Editorial Board of the Status of Livestock and Poultry Genetic Resources in China. Report on domestic animal genetic resources in China[M]. Beijing: China Agriculture Press, 2004. (in Chinese)
[6]
WEIGEND S, VEF E, WESCH G, et al. Conception for conserving genetic resources in poultry in Germany[J]. Arch Gefügelk, 1995, 59(6): 327-334.
[7]
ELSHIRE R J, GLAUBITZ J C, SUN Q, et al. A robust, simple Genotyping-By-Sequencing (GBS) approach for high diversity species[J]. PLoS One, 2011, 6(5): e19379. DOI:10.1371/journal.pone.0019379
[8]
POLAND J A, RIFE T W. Genotyping-by-sequencing for plant breeding and genetics[J]. Plant Genome, 2012, 5(3): 92-102.
[9]
PÉRTILLE F, GUERRERO-BOSAGNA C, DA SILVA V H, et al. High-throughput and cost-effective chicken genotyping using next-generation sequencing[J]. Sci Rep, 2016, 6: 26929. DOI:10.1038/srep26929
[10]
CHEN Q, MA Y F, YANG Y M, et al. Genotyping by genome reducing and sequencing for outbred animals[J]. PLoS One, 2013, 8(7): e67500. DOI:10.1371/journal.pone.0067500
[11]
DE DONATO M, PETERS S O, MAITCHELL S E, et al. Genotyping-by-sequencing (GBS): a novel, efficient and cost-effective genotyping method for cattle using next-generation sequencing[J]. PLoS One, 2013, 8(5): e62137. DOI:10.1371/journal.pone.0062137
[12]
董世武, 王天骄, 刘华淼, 等. 基于GBS技术对梅花鹿、马鹿及其杂交后代基因组SNP特征的分析[J]. 畜牧兽医学报, 2019, 50(12): 2422-2430.
DONG S W, WANG T J, LIU H M, et al. Analysis of genomic SNP characteristics of sika deer, red deer and their hybrid progenies based on GBS technology[J]. Acta Veterinaria et Zootechnica Sinica, 2019, 50(12): 2422-2430. DOI:10.11843/j.issn.0366-6964.2019.12.006 (in Chinese)
[13]
KATSUMURA T, ODA S, MITANI H, et al. Medaka population genome structure and demographic history described via genotyping-by-sequencing[J]. G3 (Bethesda), 2019, 9(1): 217-228. DOI:10.1534/g3.118.200779
[14]
DORANT Y, BENESTAN L, ROUGEMONT Q, et al. Comparing Pool-seq, rapture, and GBS genotyping for inferring weak population structure: the American lobster (Homarus americanus) as a case study[J]. Ecol Evol, 2019, 9(11): 6606-6623. DOI:10.1002/ece3.5240
[15]
GANDOMKAR H, SHEKARABI S P H, ALIABDOLHAY H, et al. Characterization of novel genotyping-by-sequencing (GBS)-based simple sequence repeats (SSRs) and their application for population genomics of Capoeta aculeata (Valenciennes, 1844)[J]. Mol Biol Rep, 2021, 48(9): 6471-6480. DOI:10.1007/s11033-021-06653-x
[16]
VAN DER AUWERA G A, CARNEIRO M O, HARTL C, et al. From FastQ data to high-confidence variant calls: the genome analysis toolkit best practices pipeline[J]. Curr Protoc Bioinform, 2013, 43(1110): 11.10.1-11.10.33.
[17]
ALEXANDER D H, NOVEMBRE J, LANGE K. Fast model-based estimation of ancestry in unrelated individuals[J]. Genome Res, 2009, 19(9): 1655-1664. DOI:10.1101/gr.094052.109
[18]
MALOMANE D K, SIMIANER H, WEIGEND A, et al. The SYNBREED chicken diversity panel: a global resource to assess chicken diversity at high genomic resolution[J]. BMC Genomics, 2019, 20(1): 345. DOI:10.1186/s12864-019-5727-9
[19]
李国辉, 魏清宇, 张会永, 等. 基于RAD-seq技术分析边鸡保种群体的遗传进化[J]. 农业生物技术学报, 2019, 27(12): 2198-2206.
LI G H, WEI Q Y, ZHANG H Y, et al. The analysis on genetic evolution of Bian chicken (Gallus gallus) population based on RAD-seq technology[J]. Journal of Agricultural Biotechnology, 2019, 27(12): 2198-2206. (in Chinese)
[20]
SOGLIA D, SARTORE S, MAIONE S, et al. Growth performance analysis of two Italian slow-growing chicken breeds: Bianca di Saluzzo and Bionda Piemontese[J]. Animals, 2020, 10(6): 969. DOI:10.3390/ani10060969
[21]
HILLEL J, GROENEN M A M, TIXIER-BOICHARD M, et al. Biodiversity of 52 chicken populations assessed by microsatellite typing of DNA pools[J]. Genet Select Evolut, 2003, 35(6): 533-557. DOI:10.1186/1297-9686-35-6-533
[22]
KIRIKCI K, CAM M A, MERCAN L. Genetic diversity and relationship among indigenous Turkish Karayaka sheep subpopulations[J]. Arch Anim Breed, 2020, 63(2): 269-275. DOI:10.5194/aab-63-269-2020
[23]
CENDRON F, MASTRANGELO S, TOLONE M, et al. Genome-wide analysis reveals the patterns of genetic diversity and population structure of 8 Italian local chicken breeds[J]. Poult Sci, 2021, 100(2): 441-451. DOI:10.1016/j.psj.2020.10.023
[24]
AZIMU W, MANATBAY B, LI Y, et al. Genetic diversity and population structure analysis of eight local chicken breeds of Southern Xinjiang[J]. Br Poult Sci, 2018, 59(6): 629-635. DOI:10.1080/00071668.2018.1523537
[25]
CHEN G H, WU X S, WANG D Q, et al. Cluster analysis of 12 Chinese native chicken populations using microsatellite markers[J]. Asian Austral J Anim Sci, 2004, 17(8): 1047-1052. DOI:10.5713/ajas.2004.1047
[26]
VIALE E, ZANETTI E, ÖZDEMIR D, et al. Development and validation of a novel SNP panel for the genetic characterization of Italian chicken breeds by next-generation sequencing discovery and array genotyping[J]. Poult Sci, 2017, 96(11): 3858-3866. DOI:10.3382/ps/pex238
[27]
ABDALHAG M A, ZHANG T, FAN Q C, et al. Single nucleotide polymorphisms associated with growth traits in Jinghai yellow chickens[J]. Genet Mol Res, 2015, 14(4): 16169-16177. DOI:10.4238/2015.December.8.6
[28]
ZANETTI E, DE MARCHI M, DALVIT C, et al. Genetic characterization of local Italian breeds of chickens undergoing in situ conservation[J]. Poult Sci, 2010, 89(3): 420-427. DOI:10.3382/ps.2009-00324
[29]
王洪志, 李国辉, 张贤, 等. 基于RAD-seq简化基因组测序评价鹿苑鸡不同保种群保种现状[J]. 畜牧兽医学报, 2017, 48(5): 818-825.
WANG H Z, LI G H, ZHANG X, et al. Conservation status evaluation of Two Luyuan chicken populations based on RAD-seq simplified genome sequencing technology[J]. Acta Veterinaria et Zootechnica Sinica, 2017, 48(5): 818-825. (in Chinese)
[30]
FU X T, DOU J Z, MAO J X, et al. RAD typing: an integrated package for accurate de novo codominant and dominant RAD genotyping in mapping populations[J]. PLoS One, 2013, 8(11): e79960. DOI:10.1371/journal.pone.0079960
[31]
KANGINAKUDRU S, METTA M, JAKATI R D, et al. Genetic evidence from Indian red jungle fowl corroborates multiple domestication of modern day chicken[J]. BMC Evolut Biol, 2008, 8: 174. DOI:10.1186/1471-2148-8-174
[32]
CHEN G H, BAO W B, SHU J T, et al. Assessment of population structure and genetic diversity of 15 Chinese indigenous chicken breeds using microsatellite markers[J]. Asian Austral J Anim Sci, 2008, 21(3): 331-339. DOI:10.5713/ajas.2008.70125
[33]
MWACHARO J M, BJØRNSTAD G, MOBEGI V, et al. Mitochondrial DNA reveals multiple introductions of domestic chicken in East Africa[J]. Mol Phylog Evolut, 2011, 58(2): 374-382. DOI:10.1016/j.ympev.2010.11.027
[34]
RASHID M A, MANJULA P, FARUQUE S, et al. Genetic diversity and population structure of indigenous chicken of Bangladesh using microsatellite markers[J]. Asian Austral J Anim Sci, 2020, 33(11): 1732-1740. DOI:10.5713/ajas.20.0189
[35]
HASSEN H, NESER F W C, DE KOCK A, et al. Study on the genetic diversity of native chickens in northwest Ethiopia using microsatellite markers[J]. Afr J Biotechnol, 2009, 8(7): 1347-1353.
[36]
SEO D W, HOQUE M R, CHOI N R, et al. Discrimination of Korean native chicken lines using fifteen selected microsatellite markers[J]. Asian Austral J Anim Sci, 2013, 26(3): 316-322. DOI:10.5713/ajas.2012.12469
[37]
YAMAMOTO Y. Genetic diversity and genetic relationship of Bangladeshi native chickens[C]//International Poultry Show & Seminar. World's Poultry Science Association-Bangladesh Branch, 2011.
[38]
高玉时, 李慧芳, 陈国宏, 等. 地方鸡种微卫星DNA指纹图谱建立与遗传多样性研究[J]. 云南农业大学学报, 2005, 20(3): 313-318.
GAO Y S, LI H F, CHEN G H, et al. Construction of microsatellite DNA fingerprinting and analysis of genetic diversity of native chicken breeds[J]. Journal of Yunnan Agricultural University, 2005, 20(3): 313-318. DOI:10.3969/j.issn.1004-390X.2005.03.002 (in Chinese)
[39]
GAO Y S, TU Y J, TONG H B, et al. AFLP fingerprinting analysis of genetic polymorphism of 12 indigenous chicken breeds[J]. Chin J Agric Biotechnol, 2006, 4(1): 33-38.
[40]
杨秀荣, 蒋和生, 徐文文, 等. 广西地方鸡品种群体遗传结构及其遗传关系分析[J]. 中国家禽, 2019, 41(3): 6-9.
YANG X R, JIANG H S, XU W W, et al. Genetic structure and relationship analysis in Guangxi native chicken breeds[J]. China Poult, 2019, 41(3): 6-9. (in Chinese)

(编辑   郭云雁)