利用高密度SNP芯片分析杭猪的群体遗传结构

doi:10.11843/j.issn.0366-6964.2022.08.009

摘要/Abstract

摘要： 旨在更好地了解杭猪保种群体的遗传多样性、亲缘关系和家系结构，有效保护和利用杭猪遗传资源。本研究使用50K SNP芯片对30头杭猪（4头公猪，26头母猪）保种个体进行了单核苷酸多态性（single nucleotide polymorphism，SNP）测定，通过Plink软件对获得的基因型数据进行质控，用于最小等位基因频率、观察杂合度、期望杂合度和多态信息含量的计算，分析杭猪保种群的遗传多样性；使用Plink计算连续性纯合片段（runs of homozygosity，ROH）、构建状态同源（identity by state，IBS）距离矩阵，并计算得到基于连续性纯合片段的近交系数；采用Gmatrix软件构建G矩阵，分析杭猪保种群的亲缘关系；使用Mega X软件构建公猪进化树，分析保种群体的家系结构。结果表明，30头杭猪共得到57 466个SNPs标记，平均检出率为0.988 3±0.000 4，通过质检SNPs数为34 156；平均最小等位基因频率为0.228±0.137，平均多态信息含量为0.255±0.098；平均观察杂合度为0.359±0.181，平均期望杂合度为0.314±0.140。杭猪保种群体平均IBS遗传距离为0.241 7±0.033 6，公猪平均IBS遗传距离为0.178 3±0.025 5；G矩阵分析结果与IBS距离矩阵一致，两个结果均表明杭猪保种群部分个体之间亲缘关系较近。30头杭猪个体共检测到828个ROHs，83.33%的ROH长度在10 Mb以内，平均长度为（6.96±9.62） Mb，个体ROH平均总长度为（191.95±201.56） Mb。基于ROH的平均近交系数为0.078±0.082，其中公猪平均近交系数为0.219±0.082，说明杭猪保种群整体的近交程度不严重，但公猪存在近交。进化树结果表明，杭猪保种群体目前只有2个家系，其中1个家系包含4头公猪，另一个家系不含公猪。综上所述，杭猪保种场的保种群公猪的血缘数量少，母猪近交程度低，公猪存在近交，需要引入新血统或创建新血统来维持家系结构的平衡以利于杭猪遗传资源的长期保存，防止杭猪遗传多样性的丢失。

关键词: 杭猪, SNP芯片, 遗传多样性, 亲缘关系, 家系结构

Abstract: The purpose of this study was to better understand the genetic diversity, relationship and family structure of Hang pigs conserved population, protect and utilize the genetic resources of Hang pigs efficiently. The 50K SNP chip was used to detect the single nucleotide polymorphism (SNP) in 30 Hang pigs (4 male and 26 female). And the quality control of SNP genotyping results was conducted by Plink software to calculated the minor allele frequency, observed heterozygosity, expected heterozygosity and polymorphism information content so as to analyze the genetic diversity of Hang pigs conserved population. The runs of homozygosity (ROH) was calculated and identity by state (IBS) distance matrix was constructed by Plink software, the inbreeding coefficient was then calculated based on ROH. The genetic relationship was analyzed according to the G matrix result constructed by Gmatrix software. The boar phylogenetic tree which constructed by Mega X software was used to analyze the family structure of Hang pigs conserved population. The results showed that, 57 466 SNPs were detected in 30 Hang pigs, the average genotyping rate was 0.988 3±0.000 4, there were 34 156 SNPs passing quality control. The average minor allele frequency, polymorphism information content, average observed heterozygosity and average expected heterozygosity were 0.228±0.137, 0.255±0.098, 0.359±0.181 and 0.314±0.140, respectively. The average IBS genetic distance of Hang pigs conserved population and boar were 0.241 7±0.033 6 and 0.178 3±0.025 5, respectively, and the result of G matrix was consistent with IBS genetic distance, the IBS genetic distance and G matrix results showed that some individuals in Hang pigs conserved population were closer to each other. A total of 828 ROHs were detected in 30 Hang pigs, and the length of 83.33% of these ROHs was less than 10 Mb, and the average length was (6.96±9.62) Mb, and the average total length of individual ROH was (191.95±201.56) Mb. The average inbreeding coefficient based on ROH of Hang pigs conserved population and boar were 0.078±0.082 and 0.219±0.082, respectively, indicated that the inbreeding degree of Hang pigs conserved population was not serious, but there was inbreeding in boar. The phylogenetic tree result showed that, there were only two families in Hang pigs conserved population, and one family contained all 4 boars and the other family only contained sows. In summary, the blood number of the boar of conserved population in the Hang pigs conserved farm was few, the inbreeding degree of sow was low, but there was inbreeding in boar, it was necessary to introduce new blood or create new blood to maintain the balance of family structure, so as to facilitate the long-term preservation of Hang pigs genetic resources and prevent the loss of genetic diversity.

Key words: Hang pigs, SNP chip, genetic diversity, relationship, family structure

中图分类号:

S828.2

刘晨龙, 卢丹, 周泉勇, 万明春, 郝晓东, 张献贺, 郑瑞强, 季华员. 利用高密度SNP芯片分析杭猪的群体遗传结构[J]. 畜牧兽医学报, 2022, 53(8): 2502-2513.

LIU Chenlong, LU Dan, ZHOU Quanyong, WAN Mingchun, HAO Xiaodong, ZHANG Xianhe, ZHENG Ruiqiang, JI Huayuan. Analysis of Population Genetic Structure of Hang Pigs by High Density SNP Chip[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(8): 2502-2513.

参考文献 43

[1]	QIAO R, LI X, HAN X, et al.Population structure and genetic diversity of four Henan pig populations[J].Anim Genet, 2019, 50(3):262-265.
[2]	VAN BA N, NAM L Q, DO D N, et al.An assessment of genetic diversity and population structures of fifteen Vietnamese indigenous pig breeds for supporting the decision making on conservation strategies[J].Trop Anim Health Prod, 2020, 52(3):1033-1041.
[3]	WANG X, WANG C, HUANG M, et al.Genetic diversity, population structure and phylogenetic relationships of three indigenous pig breeds from Jiangxi Province, China, in a worldwide panel of pigs[J].Anim Genet, 2018, 49(4):275-283.
[4]	BLACKBURN H D, PLANTE Y, ROHRER G, et al.Impact of genetic drift on access and benefit sharing under the Nagoya protocol:the case of the Meishan pig[J].J Anim Sci, 2014, 92(4):1405-1411.
[5]	LIU B, SHEN L Y, GUO Z X, et al.Single nucleotide polymorphism-based analysis of the genetic structure of Liangshan pig population[J].Anim Biosci, 2021, 34(7):1105-1115.
[6]	XU P, WANG X P, NI L G, et al.Genome-wide genotyping uncovers genetic diversity, phylogeny, signatures of selection, and population structure of Chinese Jiangquhai pigs in a global perspective[J].J Anim Sci, 2019, 97(4):1491-1500.
[7]	JOAQUIM L B, CHUD T C S, MARCHESI J A P, et al.Genomic structure of a crossbred Landrace pig population[J]. PLoS One, 2019, 14(2):e0212266.
[8]	ZHOU M L, WANG G F, CHEN M H, et al.Genetic diversity and population structure of sheep (Ovisaries) in Sichuan, China[J]. PLoS One, 2021, 16(9):e0257974.
[9]	王淑辉.利用基因芯片对阉割引起牛生长发育和脂肪沉积变化分子机理的研究[D].北京:中国农业科学院, 2008.WANG S H.The study on the molecular mechanism of growth and fatty deposition by genechip in beef cattle[D]. Beijing:Chinese Academy of Agricultural Sciences, 2008.(in Chinese)
[10]	STRILLACCI M G, VEGA-MURILLO V E, ROMÁN-PONCE S I, et al.Looking at genetic structure and selection signatures of the Mexican chicken population using single nucleotide polymorphism markers[J].Poult Sci, 2018, 97(3):791-802.
[11]	NIE C S, ALMEIDA P, JIA Y X, et al.Genome-wide single-nucleotide polymorphism data unveil admixture of Chinese indigenous chicken breeds with commercial breeds[J].Genome Biol Evol, 2019, 11(7):1847-1856.
[12]	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.
[13]	STRILLACCI M G, COZZI M C, GORLA E, et al.Genomic and genetic variability of six chicken populations using single nucleotide polymorphism and copy number variants as markers[J].Animal, 2017, 11(5):737-745.
[14]	BOUSHABA N, BOUJENANE I, MOAZAMI-GOUDARZI K, et al.Genetic diversity and relationships among six local cattle populations in semi-arid areas assessed by a bovine medium-density single nucleotide polymorphism data[J].Animal, 2019, 13(1):8-14.
[15]	SUEZAWA R, NIKADORI H, SASAKI S.Genetic diversity and genomic inbreeding in Japanese black cows in the islands of Okinawa Prefecture evaluated using single-nucleotide polymorphism array[J].Anim Sci J, 2021, 92(1):e13525.
[16]	KIM S, CHEONG H S, SHIN H D, et al.Genetic diversity and divergence among Korean cattle breeds assessed using a BovineHD single-nucleotide polymorphism chip[J].Asian-Australas J Anim Sci, 2018, 31(11):1691-1699.
[17]	KING F J M, BANGA C B, VISSER C.Genetic diversity and population structure of three native cattle populations in Mozambique[J].Trop Anim Health Prod, 2021, 53(1):117.
[18]	MEYERMANS R, GORSSEN W, WIJNROCX K, et al.Unraveling the genetic diversity of Belgian milk sheep using medium-density SNP genotypes[J].Anim Genet, 2020, 51(2):258-265.
[19]	GRASSO A N, GOLDBERG V, NAVAJAS E A, et al.Genomic variation and population structure detected by single nucleotide polymorphism arrays in Corriedale, Merino and Creole sheep[J].Genet Mol Biol, 2014, 37(2):389-395.
[20]	BERIHULAY H, LI Y, LIU X, et al.Genetic diversity and population structure in multiple Chinese goat populations using a SNP panel[J].Anim Genet, 2019, 50(3):242-249.
[21]	MOOSANEZHAD KHABISI M, ASADI FOOZI M, LV F H, et al.Genome-wide DNA arrays profiling unravels the genetic structure of Iranian sheep and pattern of admixture with worldwide coarse-wool sheep breeds[J]. Genomics, 2021, 113(6):3501-3511.
[22]	JANEŠ M, ZORC M, CUBRIC-CURIK V, et al.Population structure and genetic history of Tibetan Terriers[J].Genet Sel Evol, 2019, 51(1):79.
[23]	LETKO A, MINOR K M, JAGANNATHAN V, et al.Genomic diversity and population structure of the Leonberger dog breed[J]. Genet Sel Evol, 2020, 52(1):61.
[24]	COSGROVE E J, SADEGHI R, SCHLAMP F, et al.Genome diversity and the origin of the Arabian horse[J].Sci Rep, 2020, 10(1):9702.
[25]	蔡春波, 张雪莲, 张万峰, 等.运用SNP芯片评估马身猪保种群体的遗传结构[J].畜牧兽医学报, 2021, 52(4):920-931.CAI C B, ZHANG X L, ZHANG W F, et al.Evaluation of genetic structure in Mashen pigs conserved population based on SNP chip[J].Acta Veterinaria et Zootechnica Sinica, 2021, 52(4):920-931.(in Chinese)
[26]	黄树文, 张哲, 陈赞谋, 等.广东省现有5个地方猪种基于SNP芯片的遗传多样性分析[J].中国畜牧杂志, 2018, 54(6):33-37.HUANG S W, ZHANG Z, CHEN Z M, et al.Geneticdiversity analysis of five Cantonese indigenous pigs based on SNP chip[J].Chinese Journal of Animal Science, 2018, 54(6):33-37.(in Chinese)
[27]	刘彬, 沈林園, 陈映, 等.基于SNP芯片分析青峪猪保种群体的遗传结构[J].畜牧兽医学报, 2020, 51(2):260-269.LIU B, SHEN L Y, CHEN Y, et al.Analysis of genetic structure of conservation population in Qingyu pig based on SNP chip[J].Acta Veterinaria et Zootechnica Sinica, 2020, 51(2):260-269.(in Chinese)
[28]	李文.不同饲喂模式对杭猪生产性能影响的研究[D].南昌:江西农业大学, 2016.LI W.The study on production performance of different raise model to Xiushui Hang-pig[D].Nanchang:Jiangxi Agricultural University, 2016.(in Chinese)
[29]	CHANG C C, CHOW C C, TELLIER L C, et al.Second-generation PLINK:rising to the challenge of larger and richer datasets[J].Giga Science, 2015, 4:7.
[30]	SU G, CHRISTENSEN O F, JANSS L, et al.Comparison of genomic predictions using genomic relationship matrices built with different weighting factors to account for locus-specific variances[J].J Dairy Sci, 2014, 97(10):6547-6559.
[31]	KUMAR S, STECHER G, LI M, et al.MEGA X:molecular evolutionary genetics analysis across computing platforms[J].Mol Biol Evol, 2018, 35(6):1547-1549.
[32]	LARSON G, LIU R R, ZHAO X B, et al.Patterns of East Asian pig domestication, migration, and turnover revealed by modern and ancient DNA[J].Proc Natl Acad Sci U S A, 2010, 107(17):7686-7691.
[33]	FAN B, WANG Z G, LI Y J, et al.Genetic variation analysis within and among Chinese indigenous swine populations using microsatellite markers[J].Anim Genet, 2002, 33(6):422-427.
[34]	SANCRISTOBAL M, CHEVALET C, HALEY C S, et al.Genetic diversity within and between European pig breeds using microsatellite markers[J].Anim Genet, 2006, 37(3):189-198.
[35]	LAVAL G, IANNUCCELLI N, LEGAULT C, et al.Genetic diversity of eleven European pig breeds[J].Genet Sel Evol, 2000, 32(2):187-203.
[36]	YANG S L, WANG Z G, LIU B, et al.Genetic variation and relationships of eighteen Chinese indigenous pig breeds[J].Genet Sel Evol, 2003, 35(6):657-671.
[37]	AI H S, HUANG L S, REN J.Genetic diversity, linkage disequilibrium and selection signatures in Chinese and western pigs revealed by genome-wide SNP markers[J].PLoS One, 2013, 8(2):e56001.
[38]	MUÑOZ M, BOZZI R, GARCÍA-CASCO J, et al.Genomic diversity, linkage disequilibrium and selection signatures in European local pig breeds assessed with a high density SNP chip[J].Sci Rep, 2019, 9(1):13546.
[39]	WANG Z, CHEN Q, YANG Y, et al.Genetic diversity and population structure of six Chinese indigenous pig breeds in the Taihu Lake region revealed by sequencing data[J].Anim Genet, 2015, 46(6):697-701.
[40]	KIRIN M, MCQUILLAN R, FRANKLIN C S, et al.Genomic runs of homozygosity record population history and consanguinity[J].PLoS One, 2010, 5(11):e13996.
[41]	LYNCH M, CONERY J, BURGER R.Mutationaccumulation and the extinction of small populations[J].Am Nat, 1995, 146(4):489-518.
[42]	CHARLESWORTH D, WILLIS J H.The genetics of inbreeding depression[J].Nat Rev Genet, 2009, 10(11):783-796.
[43]	AGRAWAL A F, WHITLOCKM C W.Mutation load:the fitness of individuals in populations where deleterious alleles are abundant[J].Annu Rev Ecol Evol Syst, 2012, 43(1):115-135.