中国农业科学 ›› 2023, Vol. 56 ›› Issue (15): 3032-3039.doi: 10.3864/j.issn.0578-1752.2023.15.016

• 畜牧·兽医 • 上一篇    下一篇

基于芯片和填充测序数据的肉鸡屠宰性状基因组选择准确性评估

尹畅1(), 朱墨1, 陈艳茹1, 童世锋1, 赵桂苹2, 刘杨1()   

  1. 1 南京农业大学动物科技学院,南京 210095
    2 中国农业科学院北京畜牧兽医研究所/畜禽营养与饲养全国重点实验室,北京 100193
  • 收稿日期:2022-05-18 接受日期:2022-11-15 出版日期:2023-08-01 发布日期:2023-08-05
  • 通信作者:
    刘杨,E-mail:
  • 联系方式: 尹畅,E-mail:2021105019@stu.njau.edu.cn。
  • 基金资助:
    江苏省种业振兴揭榜挂帅项目(JBGS〔2021〕026); 安徽省良种联合攻关项目(340000211260001000431); 中国农业科学院基本科研业务费(Y2020PT02)

Assessment of Genomic Selection Accuracy for Slaughter Traits in Broilers Based on Microarray and Imputed Sequencing Data

YIN Chang1(), ZHU Mo1, CHEN YanRu1, TONG ShiFeng1, ZHAO GuiPing2, LIU Yang1()   

  1. 1 College of Animal Science and Technology, Nanjing Agricultural University, Nanjing 210095
    2 Institute of Animal Sciences, Chinese Academy of Agricultural Sciences/State Key Laboratory of Animal Nutrition and Feeding, Beijing 100193
  • Received:2022-05-18 Accepted:2022-11-15 Published:2023-08-01 Online:2023-08-05

摘要:

【背景】 畜禽育种工作的核心是基因组估计育种值的准确性。不同水平的遗传标记密度对估计育种值的影响较大,随着基因分型技术的发展和高通量测序价格的下降,基于重测序数据的基因组选择研究不断涌现。理论上,标记密度更高可获得更高准确性的估计育种值。因为影响目标性状的数量性状基因座(quantitative trait loci, QTL)至少与覆盖全基因组范围的高密度标记中的一个标记处于连锁不平衡状态。所以,较高密度的标记水平,理论上标记与QTL之间的紧密连锁更好,从而保证了较高的预测准确性。但也有研究表明,填充测序数据与芯片数据相比,基因组预测的准确性提升并不明显。【目的】 利用GBLUP方法,通过比较填充测序数据和芯片数据在肉鸡屠宰性状的基因组选择准确性,为肉鸡基因组选择育种的基因分型策略提供理论依据。【方法】 依据芯片数据和填充测序(whole-genome sequence, WGS)数据,利用GBLUP方法,针对白羽肉鸡胸肌重、屠体重和腿肌重性状进行基因组预测,对其在基因组预测的准确性进行比较。首先,使用“京芯一号”鸡55 K SNP芯片对3 362只鸡进行基因分型,并从第7世代的第9批次中随机选取230只鸡进行全基因组重测序,然后利用Beagle 5.1软件将55 K SNP芯片数据填充至重测序数据水平。为避免染色体大小对填充准确性的影响,将选择鸡较大的3号染色体和较小的14号染色体来进行计算等位基因准确率(allele correct rate, CR)和基因型相关系数(correlation, Cor),并以此判断填充准确性。利用填充测序数据对3个屠宰性状的基因组育种值进行预测,并采用5-折交叉验证的方法评价预测结果的准确性、秩相关和无偏性。【结果】 两条染色体的平均等位基因准确率为0.924,平均基因型相关系数为0.885,填充准确率较高,可以用于后期基因组预测研究。SNP芯片数据基因组育种值的预测准确性在0.2194—0.2629之间,填充测序数据基因组育种值的预测准确性在0.2110—0.2695之间。与55 K SNP芯片的结果相比,填充测序数据的基因组育种值预测的准确性差异不显著。【结论】 与SNP芯片的结果相比,利用填充后的基因组数据对白羽肉鸡的3个屠宰性状(胸肌重、屠体重和腿肌)的基因组育种值预测准确性提升并不显著,该结论为畜禽遗传育种工作中的数据类型选择提供参考。

关键词: 白羽肉鸡, 屠宰性状, 基因组育种值预测, 填充测序数据, 芯片数据, 评估

Abstract:

【Background】 In the breeding work of livestock and poultry, the core of which is the accuracy of genomic estimated breeding values. Different levels of genetic marker densities have a great impact on estimated breeding values, and with the development of genotyping technology and the decrease of high-throughput sequencing prices, genomic selection studies based on sequencing data have emerged. Theoretically, higher marker density can obtain higher prediction accuracy. This is because Quantitative Trait Loci (QTL) affecting the target trait are in linkage disequilibrium with at least one of the high-density markers covering the entire genome. A higher density of marker levels theoretically ensures tight linkage between markers and QTL, thus ensuring higher prediction accuracy. However, compared with microarray data, it has also been shown that the accuracy of genomic prediction for imputed sequencing data is not significantly improved. 【Objective】 Using the GBLUP method, we compared the genomic selection accuracy of imputed sequencing data and microarray data for slaughter traits in broiler chickens to provide a theoretical basis for genotyping strategies for broiler genomic selection breeding. 【Method】 In this study, we used SNP array data and imputed whole-genome sequence level (WGS) data to perform genomic prediction for the traits of breast muscle weight, carcass weight and thigh muscle weight in white feather broilers using the GBLUP method, and then we conducted a comparative study on their accuracy in genomic prediction. First, 3 362 chickens were genotyped using the Jingxin No. 1 chicken 55 K SNP chip, and 230 chickens were randomly selected from the ninth batch of generation 7 for whole-genome resequencing, and then the 55 K SNP chip data were imputed to the resequencing data level using Beagle 5.1 software. Considering the effect of chromosome size on the filling accuracy, the larger chromosome 3 and the smaller chromosome 14 were used to calculate the allele correct rate (CR) and genotype correlation coefficient (Cor), and the imputed WGS accuracy was determined by this study. The genomic breeding values of three slaughter traits were predicted using the imputed WGS data, and the accuracy, rank correlation and unbiasedness of the prediction results were evaluated using a 5-fold cross-validation method. 【Result】 The results showed that the average allelic accuracy of the two chromosomes was 0.924 and the average genotype correlation was 0.885, and the imputed WGS accuracy was high enough to be used for genomic prediction studies at a later stage. The accuracy of the predicted genomic breeding values calculated from microarray data ranged from 0.2194 to 0.2629, and the accuracy of the predicted genomic breeding values calculated from imputed sequencing data ranged from 0.2110 to 0.2695. The results show that the difference in the accuracy of the prediction of genomic breeding values from the imputed sequencing data was not significant compared with the 55 K SNP chip results. 【Conclusion】 Compared with the results of 55 K SNP microarray, the improvement in the accuracy of genomic breeding value prediction for three slaughter traits (breast muscle weight, carcass weight and leg muscle) in white feather broiler using imputed genomic level data was not significant, which provides a reference for the selection of data types in livestock genetic breeding work.

Key words: white feather broiler, slaughter traits, genomic breeding value prediction, imputed sequencing data, microarray data, assessment