中国农业科学 ›› 2023, Vol. 56 ›› Issue (14): 2798-2811.doi: 10.3864/j.issn.0578-1752.2023.14.014

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

基于全基因组SNPs分析陆丰黄牛和雷琼牛的群体结构与遗传多样性特征

童雄(), 罗威(), 闵力, 张志飞, 马新燕, 罗成龙, 陈卫东, 徐斌(), 李大刚()   

  1. 广东省农业科学院动物科学研究所/畜禽育种国家重点实验室/岭南现代农业科学与技术广东省实验室河源分中心/广东省畜禽育种与营养研究 重点实验室,广州 510640
  • 收稿日期:2023-01-04 接受日期:2023-04-23 出版日期:2023-07-16 发布日期:2023-07-21
  • 通信作者:
    徐斌,E-mail:
    李大刚,E-mail:
  • 联系方式: 童雄,E-mail:tongxiong@gdaas.cn。罗威,E-mail:luowei@gdaas.cn。童雄和罗威为同等贡献作者。
  • 基金资助:
    广东省现代农业产业技术体系创新团队项目(2022KJ114); 岭南现代农业科学与技术广东省实验室河源分中心自主科研项目(DT20220023); 2022年省级乡村振兴战略专项资金种业振兴项目(2022-XBH-00-008)

Population Structure and Genetic Diversity of Lufeng Cattle and Leiqiong Cattle Based on Genome-Wide SNPs

TONG Xiong(), LUO Wei(), MIN Li, ZHANG ZhiFei, MA XinYan, LUO ChengLong, CHEN WeiDong, XU Bin(), LI DaGang()   

  1. Institute of Animal Science, Guangdong Academy of Agricultural Sciences/State Key Laboratory of Livestock and Poultry Breeding/ Heyuan Branch of Guangdong Laboratory for Lingnan Modern Agriculture/Guangdong Key Laboratory of Animal Breeding and Nutrition, Guangzhou 510640
  • Received:2023-01-04 Accepted:2023-04-23 Published:2023-07-16 Online:2023-07-21

摘要:

【目的】研究陆丰黄牛和雷琼牛与世界不同地域家牛中的系统发育关系,解析不同家牛群体的遗传多样性,为地方家牛资源的鉴定与保护奠定理论基础。【方法】采集12头陆丰黄牛和17头雷琼牛的组织样品进行全基因组重测序,整合世界范围内分布的24个品种92个体的NCBI公共基因组数据,共计25个品种121个个体的信息开展群体遗传学研究。选用Bos taurus ARS-UCD1.2作为参考基因组,经基因组序列比对与质量控制获取高质量Reads,应用GATK软件检测基因组SNPs。进一步基于群体SNPs,利用系统进化树构建、PCA聚类和Admixture评估进行群体结构分析,通过核苷酸多样性(Pi)、杂合度(Hp)和连锁不平衡(LD)水平研究群体的遗传多样性。【结果】29头广东地方牛经基因组测序获取6 905 944 306个Clean Reads,每个样本平均基因组覆盖率为97.99%,平均测序深度分别为12.78×。整合NCBI公共基因组数据,经质控后共检测出14 664 391个群体SNPs。整合系统进化树、PCAAdmixture的结果发现普通牛与瘤牛分化,瘤牛群体存在中国瘤牛与印度瘤牛分化,普通牛群体中东北亚普通牛(韩牛、延边牛)和西藏牛与欧洲普通牛分化,温岭高峰牛和舟山牛从中国瘤牛群体中分化。陆丰黄牛和雷琼牛均属于纯正的中国瘤牛,但陆丰黄牛与皖南牛之间、雷琼牛与吉安牛之间呈现最近的亲缘关系,说明陆丰黄牛与地域临近的雷琼牛属于两个独立的品种。部分陆丰黄牛与雷琼牛均存在欧洲普通牛和东北亚普通牛的血统混杂,且混杂比例较高,说明这两个品种牛急需加强群体内的提纯复壮。相较于欧洲普通牛和韩牛,中国家牛群体的LD衰减速率更快,核苷酸多样性Pi和杂合度Hp更高,说明其遗传多样性更丰富。相较于其他中国家牛,陆丰黄牛和雷琼牛的LD水平更低,杂合度Hp更高,且核苷酸多样性Pi和杂合度Hp的密度分布更为集中,说明两者受人工选择强度较低,且维持较高的群体遗传多样性。【结论】通过全基因组SNPs标记,系统解析了陆丰黄牛与雷琼牛的群体遗传结构和多样性特征,为这两个品种独立分类及其保护利用提供数据支撑。

关键词: 陆丰黄牛, 雷琼牛, 全基因组SNPs, 群体结构

Abstract:

【Objective】Phylogenetic relationship among Lufeng cattle, Leiqiong cattle and domestic cattle in different regions worldwide, and genetic diversity of different domestic cattle populations were studied, so as to lay a theoretical foundation for the identification and protection of domestic cattle resources. 【Method】Tissue samples of 12 individuals of Lufeng cattle and 17 ones of Leiqiong cattle were collected for whole genome resequencing. Combined with another 92 cattle genomes from 24 breeds worldwide available in the NCBI database, a panel of cattle genomes comprising 121 individuals were generated from 25 breeds to carry out population genetics study. Bos taurus ARS-UCD1.2 assembly was selected as the reference genome. High-quality reads were obtained by genome alignment and quality control. Genomic SNPs were detected by GATK software. Population structure was analyzed by phylogenetic tree construction, PCA clustering, and Admixture evaluation. Genetic diversity of the populations was studied by estimating nucleotide diversity (Pi), heterozygosity (Hp), and linkage disequilibrium (LD). 【Result】A total of 6 905 944 306 clean reads were obtained by genome sequencing from 29 individuals of the two cattle breeds in Guangdong. Average genome coverage and average sequencing depth of each sample were 97.99% and 12.78×, respectively. After integrating the NCBI cattle genome data, 14 664 391 population SNPs were identified. The results of phylogenetic tree, PCA and Admixture showed that a primary division was found between cattle from taurine and indicine. Moreover, indicine cattle split on Chinese and Indian cattle, and Northeast Asian cattle (Hanwoo and Yanbian) and Tibetan cattle separated from European taurine group, while Wenling cattle and Zhoushan cattle differentiated from Chinese indicine group. Lufeng cattle and Leiqiong cattle belonged to pure Chinese indicine cattle. Lufeng cattle and Wannan cattle, Leiqiong cattle and Ji'an cattle showed the closest relationship, respectively. The relationships indicated that Lufeng cattle and Leiqiong cattle in the adjacent areas belonged to two independent breeds. Some Lufeng and Leiqiong individuals were genetically admixed with European taurine and Northeast Asian taurine cattle, and the admixed proportion was high, indicating that these two breeds needed to strengthen the purification and rejuvenation within the population. Compared with European taurine cattle and Korean cattle, for Chinese domestic cattle, LD decay rate was faster, and nucleotide diversity (Pi) and heterozygosity (Hp) were higher, indicating that genetic diversity of Chinese domestic cattle was richer. Compared with other Chinese domestic cattle, LD levels of Lufeng cattle and Leiqiong cattle were lower, heterozygosity (Hp) was higher, and the density distribution of nucleotide diversity (Pi) and heterozygosity (Hp) was more concentrated, indicating that the two cattle breeds were less subject to artificial selection and maintain higher genetic diversity. 【Conclusion】Population structure and genetic diversity of Lufeng cattle and Leiqiong cattle were analyzed by genome-wide SNPs, which provided data support for independent classification and conservation of these two cattle breeds.

Key words: Lufeng cattle, Leiqiong cattle, whole-genome SNPs, population structure