研究报告

全基因组重测序解析秦川牛保种群遗传多样性和遗传结构

展开
  • 1.西北农林科技大学动物科技学院,杨凌 712100
    2.陕西省农牧良种场,宝鸡 722203
马钧,在读博士研究生,专业方向:动物遗传育种与繁殖。E-mail: Junma96@163.com

收稿日期: 2023-04-25

  修回日期: 2023-06-11

  网络出版日期: 2023-06-30

基金资助

陕西省重点研发计划项目(2023-YBNY-141);财政部与农业农村部-国家现代农业产业技术体系(CARS-37)

Analysis of genetic diversity and genetic structure of Qinchuan cattle conservation population using whole-genome resequencing

Expand
  • 1. College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
    2. Shaanxi Provincial Livestock Breeding Farm, Baoji 722203, China

Received date: 2023-04-25

  Revised date: 2023-06-11

  Online published: 2023-06-30

Supported by

Supported by the Key R&D project of Shaanxi Province(2023-YBNY-141);Ministry of Finance and Ministry of Agriculture and Rural Affairs-National Modern Agricultural Industry Technology System(CARS-37)

摘要

在畜禽资源保护中,群体遗传多样性和遗传结构是决定保种效果的重要因素。本研究采用全基因组重测序技术检测100头秦川牛(30头公牛、70头母牛)的基因组变异,通过分析群体遗传多样性、连续纯合片段(runs of homozygosity,ROH)分布特征、亲缘关系和家系结构,对秦川牛的保种效果进行了综合评估。结果显示,100头秦川牛共检测到20,968,017个高质量SNPs位点,平均最小等位基因频率为0.191±0.124,平均多态信息含量为0.279±0.131,平均观察杂合度为0.275±0.131,平均期望杂合度为0.279±0.131,表明秦川保种群遗传多样性较为丰富。秦川牛保种群体平均状态同源(identity by state,IBS)遗传距离为0.243±0.020,其中公牛为0.242±0.021,亲缘关系G矩阵结果与IBS距离矩阵结果一致,均显示秦川牛保种群部分个体间亲缘关系较近。100头秦川牛个体共检测到8258个基因组ROH,ROH总长度为9.64 GB,平均ROH长度为1.167±1.203 Mb,69.35%的ROH是长度为0.5~1 Mb的短ROH。个体平均ROH总长度为96.40 Mb。基于ROH的平均近交系数为0.039±0.039,其中30头秦川公牛的平均近交系数为0.044±0.035,表明部分公牛个体存在一定程度的近交积累。进化树结果显示,秦川牛保种群所测个体可分为8个家系,包括7个含公牛家系和1个不含公牛家系。本研究表明,秦川牛保种群的遗传多样性较为丰富,未出现较大程度近交积累,但部分个体间存在近交风险,应强化选配以确保秦川牛资源的可持续发展。

本文引用格式

马钧, 樊安平, 王武生, 张金川, 江晓军, 马瑞军, 贾社强, 刘飞, 雷初朝, 黄永震 . 全基因组重测序解析秦川牛保种群遗传多样性和遗传结构[J]. 遗传, 2023 , 45(7) : 602 -616 . DOI: 10.16288/j.yczz.23-115

Abstract

In the conservation of livestock and poultry resources, population genetic diversity and genetic structure of the conservation population are important factors affecting the effectiveness of conservation. In this study, whole-genome resequencing technology was used to detect genomic variation in 100 Qinchuan cattle (30 bulls and 70 cows). By analyzing population genetic diversity, runs of homozygosity (ROH) distribution features, kinship relationships, and family structure, the conservation effectiveness of Qinchuan cattle was comprehensively evaluated. The results showed that a total of 20,968,017 high-quality SNPs were detected in 100 Qinchuan cattle, the average minimum allele frequency was 0.191±0.124, the average polymorphic information content was 0.279±0.131, and the average observed heterozygosity was 0.275±0.131, the average expected heterozygosity is 0.279±0.131, indicating that the genetic diversity of the Qinchuan cattle conservation population is relatively rich. The average identity by state (IBS) distance of the Qinchuan conservation population was 0.243±0.020, with a value of 0.242±0.021 for the bulls. The results of the kinship G-matrix were consistent with the results of the IBS distance matrix, both showing that some individuals in the conservation population had close kinship. A total of 8258 genomic ROH were detected in 100 Qinchuan cattle, with a total length of 9.64 GB. The average length of ROH fragments was 1.167±1.203 Mb, 69.35% of the ROH were short ROH with a length of 0.5~1 Mb, and the average total length of ROH per individual was 96.40 Mb. The average inbreeding coefficient based on ROH was 0.039±0.039, with a value of 0.044±0.035 for the bulls, indicating that some bulls had a certain degree of inbreeding accumulation. The results of the phylogenetic tree combined with kinship analysis showed that the individuals in the Qinchuan cattle conservation population could be divided into eight families, including seven families with bulls and one family without bulls. This study demonstrated that the genetic diversity of the Qinchuan conservation population is relatively rich, with no significant inbreeding accumulation, but there is a risk of inbreeding among some individuals. Therefore, it is necessary to strengthen selection and mating to ensure the sustainable development of Qinchuan cattle resources.

参考文献

[1] Wang B, Ma GJ, Zhang Y, Zhou FL, Cheng ZG, Hao P. External factors and improvement model of Qinchuan cattle hybrid improvement. Animal Husbandry and Veterinary Science (Electronic Edition), 2020, (21): 39-40.
[1] 王勃, 马国际, 张岩, 周福来, 程助国, 郝鹏. 秦川牛杂交改良外部因素及改良模式. 畜牧兽医科学(电子版), 2020, (21): 39-40.
[2] Yan BY.Population genomics research on Chaidamu yellow cattle[Dissertation]. Lanzhou University, 2019.
[2] 闫碧瑶.柴达木黄牛的群体基因组学研究[学位论文]. 兰州大学, 2019.
[3] Bravo S, Larama G, Qui?ones J, Paz E, Rodero E, Sepúlveda N. Genetic diversity and phylogenetic relationship among araucana creole sheep and Spanish sheep breeds. Small Ruminant Res, 2019, 172: 23-30.
[4] Jia XB, Ding P, Chen SY, Zhao SK, Wang J, Lai SJ. Analysis of MC1R, MITF, TYR, TYRP1, and MLPH genes polymorphism in four rabbit breeds with different coat colors. Animals (Basel), 2021, 11(1): 81.
[5] Liu JX, Wei X, Deng TY, Xie R, Han JL, Du LX, Zhao FP, Wang LX. Genome-wide scan for run of homozygosity and identification of corresponding candidate genes in sheep populations. Acta Vet Zootech Sin, 2019, 50(8): 1554-1566.
[5] 刘家鑫, 魏霞, 邓天宇, 谢锐, 韩建林, 杜立新, 赵福平, 王立贤. 绵羊全基因组ROH检测及候选基因鉴定. 畜牧兽医学报, 2019, 50(8): 1554-1566.
[6] Broman KW, Weber JL. Long homozygous chromosomal segments in reference families from the centre d'Etude du polymorphisme humain. Am J Hum Genet, 1999, 65(6): 1493-1500.
[7] Al-Mamun HA, Clark SA, Kwan P, Gondro C. Genome- wide linkage disequilibrium and genetic diversity in five populations of Australian domestic sheep. Genet Sel Evol, 2015, 47: 90.
[8] Purfield DC, Berry DP, McParland S, Bradley DG. Runs of homozygosity and population history in cattle. BMC Genet, 2012, 13: 70.
[9] Browning SR, Browning BL. Identity by descent between distant relatives: detection and applications. Annu Rev Genet, 2012, 46: 617-633.
[10] Liu G, Sun FZ, Zhu FX, Feng HY, Han X. Runs of homozygosity and its application on livestock genome study. Heredity(Beijing), 2019, 41(4): 304-317.
[10] 刘刚, 孙飞舟, 朱芳贤, 冯海永, 韩旭. 连续性纯合片段在畜禽基因组研究中的应用. 遗传, 2019, 41(4): 304-317.
[11] Marras G, Gaspa G, Sorbolini S, Dimauro C, Ajmone- Marsan P, Valentini A, Williams JL, Macciotta NPP. Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy. Anim Genet, 2015, 46(2): 110-121.
[12] Ma J, Gao X, Li JY, Gao HJ, Wang ZZ, Zhang LP, Xu LY, Gao H, Li HW, Wang YH, Zhu B, Cai WT, Wang CY, Chen Y. Assessing the genetic background and selection signatures of Huaxi cattle using high-density SNP array. Animals (Basel), 2021, 11(12): 3469.
[13] Van der Nest MA, Hlongwane N, Hadebe K, Chan WY, van der Merwe NA, De Vos L, Greyling B, Kooverjee BB, Soma P, Dzomba EF, Bradfield M, Muchadeyi FC. Breed ancestry, divergence, admixture, and selection patterns of the Simbra crossbreed. Front Genet, 2020, 11: 608650.
[14] Zhang SJ, Yao Z, Li XM, Zhang ZJ, Liu X, Yang P, Chen NB, Xia XT, Lyu SJ, Shi QT, Wang EY, Ru BR, Jiang Y, Lei CZ, Chen H, Huang YZ. Assessing genomic diversity and signatures of selection in Pinan cattle using whole-genome sequencing data. BMC Genomics, 2022, 23(1): 460.
[15] Cai CB, Zhang XL, Zhang WF, Yang Y, Gao PF, Guo XH, Li BG, Cao GQ. Evaluation of genetic structure in Mashen pigs conserved population based on SNP chip. Acta Vet Zootech Sin, 2021, 52(4): 920-931.
[15] 蔡春波, 张雪莲, 张万峰, 杨阳, 高鹏飞, 郭晓红, 李步高, 曹果清. 运用SNP芯片评估马身猪保种群体的遗传结构. 畜牧兽医学报, 2021, 52(4): 920-931.
[16] Liu B, Shen LY, Chen Y, Li Q, Liao K, Guo ZY, Zhang SH, Zhu L. Analysis of genetic structure of conservation population in Qingyu pig based on SNP chip. Acta Vet Zootech Sin, 2020, 51(2): 260-269.
[16] 刘彬, 沈林園, 陈映, 李强, 廖坤, 郭芝忺, 张顺华, 朱砺. 基于SNP芯片分析青峪猪保种群体的遗传结构. 畜牧兽医学报, 2020, 51(2): 260-269.
[17] Long X, Chen L, Wu PX, Zhang TH, Pan HM, Zhang L, Wang JY, Guo ZY, Chai J. Evaluation of the genetic structure and selection signatures in Hechuan Black pigs conserved population. Acta Vet Zootech Sin, 2023, 54(5): 1854-1867.
[17] 龙熙, 陈力, 吴平先, 张廷焕, 潘红梅, 张亮, 王金勇, 郭宗义, 柴捷. 合川黑猪保种群遗传结构及选择信号分析. 畜牧兽医学报, 2023, 54(5): 1854-1867.
[18] Ma SL, Li XW, Li X, Xie SQ, Liu YL, Tang J, Jiang MF. Assessment of genetic structure of 3 Maiwa yak preserved populations based on genotyping-by-sequencing technology. Acta Prataculturae Sinica, 2022, 31(9): 183-194.
[18] 马士龙, 李小伟, 李响, 谢书琼, 刘益丽, 唐娇, 江明锋. 基于GBS简化基因组测序评估3个麦洼牦牛保种群的遗传结构研究. 草业学报, 2022, 31(9): 183-194.
[19] Shi R, Zhang Y, Wang YC, Huang T, Lu GC, Yue T, Lu ZX, Huang XX, Wei XP, Feng ST, Chen J, Kagedeer Wulan, Abulizi Ruxianguli, Abulizi Ruxianguli. The evaluation of genomic homozygosity for Xinjiang inbred population by SNP panels. Hereditas(Beijing), 2020, 42(5): 493-506.
[19] 师睿, 张毅, 王雅春, 黄涛, 卢国昌, 岳涛, 卢振西, 黄锡霞, 卫新璞, 冯书堂, 陈军, 乌兰·卡格德尔, 茹先古丽·阿不力孜, 努尔胡马尔·木合塔尔. 利用SNP芯片信息评估新疆近交牛基因组纯合度. 遗传, 2020, 42(5): 493-506.
[20] Li YX, Yashengjiang Nasier, Sai Like Duman, Qian Y, Cao SX, Wang WL, Meng CH, Zhang J, Zhang JL. Evaluation of genetic diversity and genetic structure in Kirgiz sheep population cased on SNPs chip. Acta Vet Zootech Sin, 2023, 54(2): 572-583.
[20] 李隐侠, 牙生江·那斯尔, 赛里克·都曼, 钱勇, 曹少先, 王伟列, 孟春花, 张俊, 张建丽. SNP芯片评估柯尔克孜羊群体遗传多样性和遗传结构. 畜牧兽医学报, 2023, 54(2): 572-583.
[21] Ma KY, Han JT, Bai YQ, Li TT, Ma YJ. Genetic diversity analysis of Yongdeng seven goats based on specific-locus amplified fragment sequencing. Acta Vet Zootech Sin, 2023, 54(5): 1939-1950.
[21] 马克岩, 韩金涛, 白雅琴, 李讨讨, 马友记. 基于简化基因组测序的永登七山羊遗传多样性分析. 畜牧兽医学报, 2023, 54(5): 1939-1950.
[22] Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 2014, 30(15): 2114-2120.
[23] Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv: Genomics, 2013.
[24] Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics, 2009, 25(16): 2078-2079.
[25] McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res, 2010, 20(9): 1297- 1303.
[26] Okonechnikov K, Conesa A, García-Alcalde F. Qualimap 2: advanced multi-sample quality control for high- throughput sequencing data. Bioinformatics, 2016, 32(2): 292-294.
[27] Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics, 2011, 27(21): 2987-2993.
[28] Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, Daly MJ, Sham PC. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet, 2007, 81(3): 559-575.
[29] Zhao GY.Association and prediction of traits based on genomic homozygous segments in beef cattle[Dissertation]. Chinese Academy of Agricultural Sciences, 2021.
[29] 赵国耀.基于肉牛基因组纯合片段的性状关联与预测[学位论文]. 中国农业科学院, 2021.
[30] Liu CL, Lu D, Zhou QY, Wan MC, Hao XD, Zhang XH, Zheng RQ, Ji HY. Analysis of population genetic structure of Hang pigs by high density SNP chips. Acta Vet Zootech Sin, 2022, 53(8): 2502-2513.
[30] 刘晨龙, 卢丹, 周泉勇, 万明春, 郝晓东, 张献贺, 郑瑞强, 季华员. 利用高密度SNP芯片分析杭猪的群体遗传结构. 畜牧兽医学报, 2022, 53(8): 2502-2513.
[31] Baum BR. PHYLIP: phylogeny inference package. Version 3.2. Joel felsenstein. The Quarterly Review of Biology, 1989, 64(4): 539-541.
[32] Xu LY, Bickhart DM, Cole JB, Schroeder SG, Song JZ, Tassell CPV, Sonstegard TS, Liu GE. Genomic signatures reveal new evidences for selection of important traits in domestic cattle. Mol Biol Evol, 2015, 32(3): 711-725.
[33] Edwards CE, Tessier BC, Swift JF, Bassüner B, Linan AG, Albrecht MA, Yatskievych GA. Conservation genetics of the threatened plant species Physaria filiformis (Missouri bladderpod) reveals strong genetic structure and a possible cryptic species. PLoS One, 2021, 16(3): e0247586.
[34] Liu F, Qu YK, Geng C, Wang AM, Zhang JH, Li JF, Chen KJ, Liu B, Tian HY, Yang WP, Yu YB. Analysis of the population structure and genetic diversity of the red swamp crayfish (Procambarus clarkii) in China using SSR markers. Electron J Biotechnol, 2020, 47: 59-71.
[35] Reed DH, Frankham R. Correlation between fitness and genetic diversity. Conserv Biol, 2003, 17(1): 230-237.
[36] Eldridge MDB, King JM, Loupis AK, Spencer PBS, Taylor AC, Pope LC, Hall GP. Unprecedented low levels of genetic variation and inbreeding depression in an island population of the Black-Footed Rock-Wallaby. Conserv Biol, 1999, 13(3): 531-541.
[37] Gao YH, Gautier M, Ding XD, Zhang H, Wang YC, Wang X, Faruque MO, Li JY, Ye SH, Gou X, Han JL, Lenstra JA, Zhang Y. Species composition and environmental adaptation of indigenous Chinese cattle. Sci Rep, 2017, 7(1): 16196.
[38] Zhang WG, Gao X, Zhang Y, Zhao YM, Zhang JB, Jia YT, Zhu B, Xu LY, Zhang LP, Gao HJ, Li JY, Chen Y. Genome-wide assessment of genetic diversity and population structure insights into admixture and introgression in Chinese indigenous cattle. BMC Genet, 2018, 19(1): 114.
[39] Yuan J, Xu GQ, Zhou X, Xu SP, Li S, Li WM, Liu B. SNP chip-based monitoring of population conservation effect of Tongcheng pigs. Acta Vet Zootech Sin, 2022, 53(08): 2514-2523.
[39] 袁娇, 徐国强, 周翔, 徐三平, 李胜, 黎望明, 刘榜. 基于SNP芯片监测通城猪的保种效果. 畜牧兽医学报, 2022, 53(08): 2514-2523.
[40] Botstein D, White RL, Skolnick M, Davis RW. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet, 1980, 32(3): 314-331.
[41] Forutan M, Ansari Mahyari S, Baes C, Melzer N, Schenkel FS, Sargolzaei M. Inbreeding and runs of homozygosity before and after genomic selection in North American Holstein cattle. BMC Genomics, 2018, 19(1): 98.
[42] Keller MC, Visscher PM, Goddard M. Quantification of inbreeding due to distant ancestors and its detection using dense single nucleotide polymorphism data. Genetics, 2011, 189: 237-249.
[43] Wang Y, Dong RL, Li X, Cui C, Yu GH. Analysis of the genetic diversity and family structure of the Licha Black pig population on Jiaodong Peninsula, Shandong Province, China. Animals (Basel), 2022, 12(8): 1045.
文章导航

/