近交陆川猪的遗传稳定性研究

doi:10.11843/j.issn.0366-6964.2020.08.003

摘要/Abstract

摘要： 旨在探讨高度近交陆川猪的遗传稳定性。本试验采集20头近交和6头非近交陆川猪的耳组织样品进行个体基因组重测序，通过扫描全基因组SNPs分析近交陆川猪群体的遗传多样性，并基于基因组纯合度挖掘近交陆川猪群体中特异的纯合区域（region of homozygosity，ROHs）及相关基因。结果表明，近交陆川猪群体的观察杂合度（0.373）小于期望杂合度（0.491），推测其发生了选择或者近交；系统进化树和群体遗传结构分析发现，近交陆川猪群体在近交过程中分化成为两个近交家系，且230、236、239、132、130及126号个体亲缘关系最近；在近交群体中，基因组纯合度和基因组长纯合片段比例的均值分别为59.64%和0.12，都极显著高于非近交群体53.18%和0.05（P<0.01），说明该近交陆川猪群体已达到高度近交；近交群体和非近交群体中分别有1 250和633个ROHs，重叠分析结果显示，特异存在于近交群体中的阳性纯合区域（positive region of homozygosity，pROHs）有224个，其中由于近交产生的纯合突变可能影响参与激素合成过程、胰液分泌、膜脂代谢过程、IgA产生的肠道免疫网络及纹肌细胞分化等相关的1 322个基因。结果提示，近交陆川猪群体的遗传多样性低、基因组纯合度高且已达到高度近交，说明该群体的遗传稳定性高。此外，近交陆川猪群体在近交过程中产生的纯合突变可能与机体的免疫、生长等生产性状相关。

关键词: 陆川猪, 连续性纯合片段, 基因组纯合度, 生产性状

Abstract: The aim of this study was to explore the genetic stability of high-inbreeding Luchuan pig. Ear samples were collected from 20 inbreeding and 6 non-inbreeding Luchuan pigs and genomic re-sequencing was performed for each samples. Genetic diversity of inbreeding Luchuan pig population was analyzed by scanning the whole genome SNPs. The specific ROHs and associated genes existed in inbreeding Luchuan pig population were detected based on genomic homozygosity. The results showed that observed heterozygosity was 0.373 in inbreeding Luchuan pig population, which was less than expected heterozygosity (0.491), indicating that selection or inbreeding might happen in the population. The analysis of phylogenetic trees and population's genetic structure indicated that two inbred lines were generated during the inbreeding process of Luchuan pig and the number 230, 236, 239, 132, 130 and 126 individuals had closer relationship. In inbreeding population, the average of genomic homozygosity and ratio of genomic run of homozygosity were 59.64% and 0.12, respectively, which were extremely significantly higher than 53.18% and 0.05 in non-inbreeding Luchuan pig population (P<0.01), which indicated that the inbreeding Luchuan pig population had reached a high level of inbreeding. One thousand two hundred and fifty and 633 ROHs were found in inbreeding and non-inbreeding Luchuan pig populations, respectively. The results of overlapping analysis suggested that 224 pROHs were specifically existed in inbreeding population. Among them, homozygous mutations generated from inbreeding process might affect 1 322 genes, which involved in hormone biosynthetic process, pancreatic secretion, membrane lipid metabolic process, intestinal immune network for IgA production and Striated muscle cell differentiation. These results indicate that the inbreeding Luchuan pig population has lower genetic diversity, higher genomic homozygosity and reaches a high inbreeding level, which illustrate that the population has high genetic stability. Moreover, the homozygous mutations generated from inbreeding process might be associated with immune, growth and other production traits.

Key words: Luchuan pig, runs of homozygosity, genomic homozygosity, production traits

中图分类号:

S828.2

王洪洋, 张莺莺, 涂尾龙, 吴华莉, 曹建国, 谈永松. 近交陆川猪的遗传稳定性研究[J]. 畜牧兽医学报, 2020, 51(8): 1795-1803.

WANG Hongyang, ZHANG Yingying, TU Weilong, WU Huali, CAO Jianguo, TAN Yongsong. Genetic Stability of Inbreeding Luchuan Pig[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(8): 1795-1803.

参考文献 35

[1]	张沅.家畜育种学[M].北京:中国农业出版社,2001.ZHANG Y.Animal breeding[M].Beijing:China Agricultural Press,2001.(in Chinese)
[2]	BOSSE M,MEGENS H J,DERKS M F L,et al.Deleterious alleles in the context of domestication,inbreeding,and selection[J]. Evol Appl,2019,12(1):6-17.
[3]	刘晓玲.近交在畜牧业中的应用及防止衰退的措施[J].畜牧兽医杂志,2019,38(1):76-77.LIU X L.Application of inbreeding in animal husbandry and measures to prevent recession[J].Journal of Animal Science and Veterinary Medicine,2019,38(1):76-77.(in Chinese)
[4]	陈其美,葛长利,董彬.中国地方猪种种质资源的保护与利用[C]//第三届全球猪业论坛暨第十一届(2013)中国猪业发展大会论文集.武汉:中国畜牧业协会,2013:160-163.CHEN Q M,GE C L,DONG B.Conservation and utilization of local pig germplasm resources in China[C]//Proceedings of The 3rd Global Pig Forum and the 11th (2013) China Swine Industry Development Conference.Wuhan:China Animal Husbandry Association,2013:160-163.(in Chinese)
[5]	GIBSON J,MORTON N E,COLLINS A.Extended tracts of homozygosity in outbred human populations[J].Hum Mol Genet, 2006, 15(5):789-795.
[6]	BERTOLINI F,CARDOSO T F,MARRAS G,et al.Genome-wide patterns of homozygosity provide clues about the population history and adaptation of goats[J].Genet Sel Evol,2018,50(1):59-70.
[7]	YANG B,CUI L L,PEREZ-ENCISO M,et al.Genome-wide SNP data unveils the globalization of domesticated pigs[J].Genet Sel Evol, 2017,49(1):71.
[8]	CEBALLOS F C,JOSHI P K,CLARK D W,et al.Runs of homozygosity:windows into population history and trait architecture[J]. Nat Rev Genet, 2018,19(4):220-234.
[9]	SZPIECH Z A,XU J S,PEMBERTON T J,et al.Long runs of homozygosity are enriched for deleterious variation[J].Am J Hum Genet,2013,93(1):90-102.
[10]	BANDRÉS-CIGA S,PRICE T R,BARRERO F J,et al.Genome-wide assessment of Parkinson's disease in a southern Spanish population[J].Neurobiol Aging,2016,45:213.e3-213.e9.
[11]	PRYCE J E,HAILE-MARIAM M,GODDARD M E,et al.Identification of genomic regions associated with inbreeding depression in Holstein and Jersey dairy cattle[J].Genet Sel Evol,2014,46(1):71.
[12]	STAFUZZA N B,SILVA R M D O,PERIPOLLI E,et al.Genome-wide association study provides insights into genes related with horn development in Nelore beef cattle[J].PLoS One,2018,13(8):e0202978.
[13]	ZAVAREZ L B,UTSUNOMIYA Y T,CARMO A S,et al.Assessment of autozygosity in Nellore cows (Bos indicus) through high-density SNP genotypes[J].Front Genet,2015,6:5.
[14]	MARCHESI J A P,BUZANSKAS M E,CANTÃO M E,et al.Relationship of runs of homozygosity with adaptive and production traits in a paternal broiler line[J].Animal,2018,12(6):1126-1134.
[15]	BOLGER A M,LOHSE M,USADEL B.Trimmomatic:a flexible trimmer for illumina sequence data[J]. Bioinformatics, 2014, 30(15):2114-2120.
[16]	LI H.A statistical framework for SNP calling,mutation discovery,association mapping and population genetical parameter estimation from sequencing data[J].Bioinformatics,2011,27(21):2987-2993.
[17]	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.
[18]	BARBATO M,OROZCO-TERWENGEL P,TAPIO M,et al.SNeP:a tool to estimate trends in recent effective population size trajectories using genome-wide SNP data[J].Front Genet,2015,6:109.
[19]	PURCELL S,NEALE B,TODD-BROWN K,et al.PLINK:a tool set for whole-genome association and population-based linkage analyses[J].Am J Hum Genet,2007,81(3):559-575.
[20]	杨湛澄,黄河天,闫青霞,等.利用高密度SNP标记分析中国荷斯坦牛基因组近交[J].遗传,2017,39(1):41-47.YANG Z C,HUANG H T,YAN Q X,et al.Estimation of genomic inbreeding coefficients based on high-density SNP markers in Chinese Holstein cattle[J].Hereditas,2017,39(1):41-47.(in Chinese)
[21]	冯书堂,李奎,刘岚,等.小型猪近交系新品种的培育与开发利用[J].农业生物技术学报,2015,23(2):274-280.FENG S T,LI K,LIU L,et al.Cultivation and application of miniature pig (Sus scrofa) inbred[J].Journal of Agricultural Biotechnology,2015,23(2):274-280.(in Chinese)
[22]	CARLSON M.Genome wide annotation for Pig[DB].2019,doi:10.18129/B9.bioc.org.Ss.eg.db.
[23]	YU G C,WANG L G,HAN Y Y,et al.clusterPro-filer:an R package for comparing biological themes among gene clusters[J]. OMICS,2012,16(5):284-287.
[24]	FANG X D,MOU Y L,HUANG Z Y,et al.The sequence and analysis of a Chinese pig genome[J].GigaScience,2012,1(1):16.
[25]	BOSSE M,MEGENS H J,MADSEN O,et al.Regions of homozygosity in the porcine genome:consequence of demography and the recombination landscape[J].PLoS Genet,2012,8(11):e1003100.
[26]	CASAS-CARRILLO E,KIRKPATRICK B W,PRILL-ADAMS A,et al.Relationship of growth hormone and insulin-like growth factor-1 genotypes with growth and carcass traits in swine[J].Anim Genet,1997,28(2):88-93.
[27]	HOLZENBERGER M,DUPONT J,DUCOS B,et al.IGF-1 receptor regulates lifespan and resistance to oxidative stress in mice[J].Nature,2003,421(6919):182-187.
[28]	KENNES Y M,MURPHY B D,POTHIER F,et al.Characterization of swine leptin (LEP) polymorphisms and their association with production traits[J].Anim Genet,2001,32(4):215-218.
[29]	LIN L,FLISIKOWSKI K,SCHWARZENBACHER H,et al.Characterization of the porcine AMPK alpha 2 catalytic subunit gene (PRKAA2):genomic structure,polymorphism detection and association study[J].Anim Genet,2010,41(2):203-207.
[30]	WANG L, MU Y, XU L, et al. Genomic analysis reveals specific patterns of homozygosity and heterozygosity in inbred pigs[J]. Animals (Basel), 2019, 9(6):314-326
[31]	SCHIAVO G,BOVO S,BERTOLINI F,et al.Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds[J].Animal,2020,14(5):910-920.
[32]	MARRAS G,GASPA G,SORBOLINI S,et al.Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy[J].Anim Genet,2015,46(2):110-121.
[33]	MASTRANGELO S,TOLONE M,DI GERLANDO R,et al.Genomic inbreeding estimation in small populations:evaluation of runs of homozygosity in three local dairy cattle breeds[J].Animal,2016,10(5):746-754.
[34]	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.
[35]	张家富.陆川猪若干种质特性基础研究[D].南宁:广西大学,2012.ZHANG J F.Basic research on several breeding characteristics in Luchuan pig[D].Nanning:Guangxi University,2012.(in Chinese)