畜牧兽医学报  2023, Vol. 54 Issue (9): 3631-3641. DOI: 10.11843/j.issn.0366-6964.2023.09.005    PDF    
引用本文
袁巍, 毕欢, 张雨丹, 张依裕, 顾晓龙, 杨红文, 陈伟. 全基因组选择信号解析剑白香猪和从江香猪的遗传差异[J]. 畜牧兽医学报, 2023, 54(9): 3631-3641.  
YUAN Wei, BI Huan, ZHANG Yudan, ZHANG Yiyu, GU Xiaolong, YANG Hongwen, CHEN Wei. Deciphering Genome-wide Selection Signals Reveals Genetic Differences between Jianbai and Congjiang Xiang Pigs[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(9): 3631-3641.  
全基因组选择信号解析剑白香猪和从江香猪的遗传差异
袁巍1,2,3, 毕欢1,2,3, 张雨丹1,2,3, 张依裕1,2,3, 顾晓龙4, 杨红文5, 陈伟1,2,3     
1. 贵州大学动物科学学院 高原山地动物遗传育种与繁殖教育部重点实验室, 贵阳 550025;
2. 贵州省动物遗传育种与繁殖重点实验室, 贵阳 550025;
3. 贵州大学香猪研究所, 贵阳 550025;
4. 贵州省贵阳德康农牧有限公司, 贵阳 550025;
5. 贵州省畜禽遗传资源管理站, 贵阳 550025
收稿日期：2023-03-17
基金项目：贵州省科技计划项目(黔科合支撑[2022]一般087);贵州大学引进人才科研项目(贵大人基合字[2019]34号); 贵州省优秀青年科技人才培养计划黔科合平台人才([2021]5630号); 贵州省生猪产业发展项目(2021);贵州省科技支撑计划(黔科合支撑[2021]一般147)
作者简介：袁巍(1997-), 男, 湖南隆回人, 硕士生, 主要从事动物遗传育种的研究, E-mail: y-weii@qq.com.
通信作者：陈伟, 主要从事动物遗传育种的研究, E-mail: chenweigzu@163.com.
摘要：旨在通过全基因组选择信号分析揭示剑白香猪和从江香猪的遗传差异, 提高对从江香猪和剑白香猪种群的保护状况、遗传多样性、群体结构和和适应性进化的认识。本研究随机选取2岁龄从江香猪和剑白香猪各5头, 并分为两组。通过全基因组重测序对从江香猪和剑白香猪测序, 将测序数据利用bwa软件比对到猪参考基因组Sscrofa 11.1上, 随后使用GATK软件进行变异检测和SNP、INDEL过滤提取。同时, 利用Plink和SNeP软件进行群体遗传多样性分析, 包括群体有效大小、多态位点比例、观测杂合度、期望杂合度。使用vcftools软件分析群体选择信号, 并通过KEGG和GO分析受选择区域基因。经过SNP和INDEL分析, 发现从江香猪和剑白香猪群体共有的SNPs位点数为16 158 002个, 从江香猪特有的为6 863 992个, 剑白香猪特有的为4 275 660个; 共有的INDEls位点数为3 275 375个, 从江香猪特有的为1 674 081个, 剑白香猪特有的为1 114 248个。从采样起算, 从江香猪和剑白香猪的群体在13代前的有效大小均为18头。其中, 从江香猪的多态性位点比例更高, 达到了0.875 7, 而二者群体的观测杂合度均高于期望杂合度。剑白香猪和从江香猪的核苷酸多样性均小于0.5%, 群体内选择检验Tajima's D值均大于0, ROD值为0.553 1, Fst值为0.736 2。KEGG和GO富集分析表明, 剑白香猪和从江香猪的遗传差异涉及多种通路和基因功能。包括涉及轴突引导、ECM与受体的相互作用、局部粘附、PI3K-Akt信号通路、硫辛酸代谢、cGMP-PKG信号通路、MAPK信号传导途径、核苷酸切除修复、阿佩林信号通路、心律失常性右室心肌病等通路。剑白香猪和从江香猪群体之间存在一定的群体多态性损失, 群体分歧度较高, 两亚群间存在明显的种群分化, 并且二者基因组内存在大量的中等频率等位基因, 需要提高从江香猪和剑白香猪群体内的遗传多样性, 减少外来血缘引入导致纯种度降低的风险。
关键词：全基因组    选择信号    香猪    遗传差异    
Deciphering Genome-wide Selection Signals Reveals Genetic Differences between Jianbai and Congjiang Xiang Pigs
YUAN Wei1,2,3, BI Huan1,2,3, ZHANG Yudan1,2,3, ZHANG Yiyu1,2,3, GU Xiaolong4, YANG Hongwen5, CHEN Wei1,2,3     
1. Key Laboratory of Animal Genetics, Breeding and Reproduction in the Plateau Mountain Region of Ministry of Education, College of Animal Science, Guizhou University, Guiyang 550025, China;
2. Key Laboratory of Animal Genetics, Breeding and Reproduction, Guiyang 550025, China;
3. Institute of Xiang Pig Research, Guizhou University, Guiyang 550025, China;
4. Guizhou Dekang Agriculture and Animal Husbandry Co. Ltd., Guiyang 550025, China;
5. Guizhou Provincial Station for Animal Genetic Resources Management, Guiyang 550025, China
Corresponding author: CHEN Wei, E-mail: chenweigzu@163.com.
Abstract: The purpose was to unveil the genetic disparities between the Jianbai and Congjiang Xiang pigs breeds, enhance the understanding of their conservation status, genetic diversity, population structure, and adaptive evolution through whole-genome selection signal analysis. The 2-year-old Congjiang and Jianbai Xiang pigs were randomly selected and divided into two groups of 5 pigs each. Whole-genome resequencing was performed on Jianbai and Congjiang Xiang pigs, and the sequencing data was aligned to the pig reference genome Sscrofa 11.1 using the bwa software. Subsequently, the GATK software was used for variant detection and filtering of SNPs and INDELs. Additionally, population genetic diversity analysis, including effective population size, proportion of polymorphic loci, observed heterozygosity, and expected heterozygosity, was carried out using Plink and SNeP software. The population selection signals were analyzed using the vcftools software, and the selected genomic regions were subjected to KEGG and GO pathway analyses. After SNP and INDEL analysis, it was found that there were 16 158 002 common SNPs in the populations of Jianbai and Congjiang Xiang pigs, with 6 863 992 SNPs unique to Congjiang and 4 275 660 unique to Jianbai. The common INDELs were 3 275 375, with 1 674 081 INDELs unique to Congjiang and 1 114 248 unique to Jianbai. Starting from the sampling time point, the effective population size of both Congjiang and Jianxiang pig populations could be traced back to 18 individuals 13 generations ago. The Congjiang breed displayed a higher proportion of polymorphic loci, with a value of 0.875 7, whereas both populations had observed heterozygosity values exceeding the expected heterozygosity.Nucleotide diversity was less than 0.5% for both breeds, and Tajima's D was greater than 0, with ROD values of 0.553 1 and Fst values of 0.736 2. KEGG and GO enrichment analyses showed that the genetic differences between Jianbai and Congjiang Xiang pigs involved multiple pathways and gene functions, including axon guidance, ECM-receptor interaction, focal adhesion, PI3K-Akt signaling pathway, lipoic acid metabolism, cGMP-PKG signaling pathway, MAPK signaling pathway, nucleotide excision repair, aspirin signaling pathway, and arrhythmogenic right ventricular cardiomyopathy. There are a certain loss of population polymorphism and a high degree of population divergence between the Jianbai and Congjiang pig populations. There is clear population differentiation between the two breeds, and both genomes exhibit a substantial number of alleles at intermediate frequencies. It is necessary to enhance the genetic diversity within the Congjiang and Jianbai Xiang pig populations and reduce the risk of decreased purity due to the introduction of foreign lineages.
Key words: whole genome    selection signals    Xiang pig    genetic differences    

猪的养殖在我国历史久远、分布广泛、数量众多。我国复杂的地貌环境、多元化的饮食文化带来的多样化选育目标是中国形成多个地方猪种质资源的重要原因，多样化选育目标致使我国是世界上地方猪种质资源最丰富的国家，约占世界猪种质资源的1/3[1]。许多本地猪品种都有自己独特的特点，例如抗病性、饲料效率和高生育率。贵州香猪是一种小型猪种，包括从江香猪和剑白香猪两种类型。在贵州独特的喀斯特地貌和多民族混居的人文环境的影响下，形成了近交不退化的优势和基因高度纯合的特点。作为贵州地区宝贵的地方猪种遗传资源，通过全基因组选择信号解析剑白香猪和从江香猪的遗传差异，有助于提高对从江香猪和剑白香猪种群的保护状况、遗传多样性、群体结构和和适应性进化的认识，制定更有效的保护策略以维护地方猪的遗传多样性和生态系统的稳定性。由于长期以来的猪种杂交和人工选择，地方猪种的遗传多样性正在不断减少[2]，本地猪种和群体大小的数量均大幅减少[3]，因此，评估中国本地猪品种的遗传多样性和群体结构成为必要的举措，以便制定更加有效的保护策略，避免品种的基因流失。随着测序技术和计算方法的快速发展以及测序费用的大幅下降，二代测序技术已被广泛应用于家畜遗传多样性和群体结构的研究中，如猪[4-5]、牛[6]、羊[7-8]、鸡[9-10]、鸭[11]等。通过对基因组水平的研究，这些研究筛选出了大量与适应性、表型、重要经济性状等相关的候选基因和遗传标记。这些研究结果不仅有助于人们深刻了解家畜物种的起源、驯化机制、人工选择和自然选择，而且还为家畜的遗传改良和育种提供了基础。目前，对于从江香猪和剑白香猪的研究主要集中在基因的表达[12]、经济性状[13]、生产性能[14]上，在基因组层面上，对剑白香猪和从江香猪的遗传差异研究较少。因此，了解从江香猪和剑白香猪种群的保护状况、遗传多样性、群体结构和适应性进化等方面的信息非常重要。通过全基因组选择信号分析，揭示了剑白香猪和从江香猪的遗传差异，为实现上述目标提供重要的基础和支持。

1 材料与方法 1.1 样品采集

随机选取贵州保种场2岁龄从江香猪和剑白香猪各5头，采集其耳组织，用DEPC水处理后放入无RNA酶的冻存管中，液氮中冷冻，后转入-80 ℃的超低温冰箱储存，用于提取DNA。从江香猪及剑白香猪如图 1所示。

图 1 从江香猪和剑白香猪照片 Fig. 1 Photos of Congjiang and Jianbai Xiang pigs
1.2 建库测序和数据处理

委托北京诺禾致源生物信息科技有限公司完成从江香猪与剑白香猪的总DNA提取和建库测序，主要测序平台为Illumina HiSeq 2500，测序策略为PE125，深度为5 x~10 x。用Fsatp v0.19.10软件[15]对测序原始数据进行过滤，去掉adaptor和低质量序列，将过滤后数据用bwa v0.7.17-r1188软件比对到猪参考基因组Sscrofa 11.1上，并进行排序和去重复处理，获得bam文件用于变异检测。

1.3 变异检测和过滤

使用GATK v4.0(https://gatk.broadinstitute.org/hc/en-us)中的HaplotypeCaller工具分别对从江香猪和剑白香猪的bam文件进行变异检测，分别合并每个猪种的gvcf文件成一个包含群体变异的群体gvcf文件。使用GATK 4.0中的Variant- Filtration工具过滤SNP与INDEL，对SNPs与INDELs的质控过滤参数为QD < 2.0 ||MQ < 40.0 ||FS>60.0 ||SOR>3.0 ||MQRankSum < -12.5 ||ReadPos RankSum < -8.0。

1.4 遗传多样性分析

群体有效大小(Ne)、多态位点比例(PN)、观测杂合度(Ho)和期望杂合度(He)是分析种群遗传多样性的常用参数。利用SNeP v1.11软件[16]计算群体有效大小Ne，使用Plink v1.90b6.21软件[17]计算多态位点比例(PN)、观测杂合度(Ho)和期望杂合度(He)。有效群体大小Ne根据不同的SNP数据量和历史背景进行计算，因此其接近最近世代的程度也会不同。

1.5 群体选择信号分析

群体间遗传分化指数(the fixation index, Fst)、核苷酸多样性(nucleotide diversity, π)与群体内选择检验(Tajima ’s D)均由vcftools v0.1.16软件基于数据集Data1计算，窗口大小为100 000 bp，滑动距离为1 000 bp，分组信息为从江香猪(coj)和剑白香猪(jab)，最后用R包ggplot2可视化。具体参数为--fst-window-size 100 000 --fst-window-step 10 000，--window-pi 100 000 --window-pi-step 10 000，--TajimaD 100000。群体间分歧度(the reduction of diversity, ROD)值公式为1-π1/π2，其中π1为剑白香猪群体平均π值，π2为从江香猪群体平均π值。

1.6 受选择区域分析

分别取Fst和ROD结果前5%位点，提取其在染色体上位置等信息，用R编写脚本合并为Fst和ROD是否同时满足前5%位点信息表。将Fst阈值线设置为0.2，ROD阈值线设置为0.55，利用R包ggplot2绘制点图。并将同时满足Fst和ROD前5%的位点使用bedtools v2.30.0合并重叠窗口位点，然后根据合并位点提取候选基因信息。具体参数为“bedtools intersect -F 0.1 -b gff3，bedtools merge”。

1.7 注释与通路分析

将已鉴定到的猪特有家族基因和显著收缩与扩张基因与KEGG、GO数据库比对，进行通路分析。KEGG通过R包clusterProfiler v4.0[18]调用KEGG数据库官网(https://www.kegg.jp/)API直接用于KEGG分析，并通过R包内置函数dotplot进行可视化分析；GO通过R包org.Ss.eg.db v3.16.0(http://bioconductor.org/packages/release/BiocViews.html#___OrgDb)进行分析，并使用clusterProfiler内置函数dotplot进行可视化分析; 最后通过在线网站基迪奥生信云工具(https://www.omicshare.com/tools/Home/Soft/cog)进行可视化。

2 结果 2.1 SNP及INDEL检测结果

从江香猪和剑白香猪中平均检测到2 100万个SNPs位点，其中从江香猪23 021 994个，剑白香猪20 433 662个。同样平均检测到460万个INDELs位点，其中从江香猪4 949 456个，剑白香猪4 389 623个。在外显子区域，均检测到的SNPs位点数约为17万个，其中从江香猪188 674个，剑白香猪164 383个；平均检测到的INDELs位点数约为10万个，其中从江香猪11 841个，剑白香猪10 613个。这些外显子区域上的SNP和INDEL变异可能会影响蛋白质的结构、功能、表达和稳定性，从而影响生物体的表型。利用R包Venn Diagram根据SNPs位置信息绘制的韦恩图显示，从江香猪和剑白香猪群体共有的SNPs位点数为16 158 002个，从江香猪特有的为6 863 992个，剑白香猪特有的为4 275 660个(图 2A)；共有的INDEls位点数为3 275 375个，从江香猪特有的为1 674 081个，剑白香猪特有的为1 114 248个(图 2B)。

图 2 从江香猪和剑白香猪SNP(A)及INDEL(B)检测结果韦恩图 Fig. 2 The Venn diagrams depicting the SNP (A) and INDEL (B) detection results of Congjiang and Jianbai Xiang pigs
2.2 从江香猪和剑白香猪群体的遗传多样性分析

利用plink软件分析群体有效大小(Ne)、多态位点比例(PN)、观测杂合度(Ho)和期望杂合度(He)，结果如表 1所示。通过plink计算的剑白香猪和从江香猪的最近世代群体有效大小为18头，最近世代为13。从江香猪的多态性位点比例为0.875 7，高于剑白香猪的多态性位点比例0.813 2。另外，从江香猪和剑白香猪各自群体内的观测杂合度均大于期望杂合度，说明二者群体内部都存在较高的遗传多样性。选择压力等因素导致群体内的杂合子数量增加，有利于个体对环境的适应和进化。

表 1 从江香猪和剑白香猪遗传多样性参数 Table 1 The genetic diversity parameters of Congjiang and Jianbai Xiang pigs
2.3 剑白香猪和从江香猪群体选择分析

以从江香猪作为参照组，对剑白香猪进行群体压力分析。通过滑动窗口计算从江香猪(coj)和剑白香猪(jab)群体的核苷酸多样性π值(图 3A, B)、群体遗传分化指数Fst(图 3C)。统计从江香猪和剑白香猪每条染色体上高于0.01阈值的核苷酸位点数量(图 3D)，结果显示从江香猪的每条染色体上均存在高核苷酸多态性位点，而剑白香猪的chr4上不存在高核苷酸多态性位点。从江香猪的总核苷酸多态性位点数为1 288个，高于剑白香猪的1 035个。Fst图中，每条染色体均存在大量高于阈值线的多态性位点，表明剑白香猪和从江香猪之间在每个染色体上都存在遗传分化情况。

图 3 从江香猪和剑白香猪的π值、Fst值以及核苷酸多态性在染色体上的分布 Fig. 3 The distribution of π values, Fst values, and nucleotide polymorphisms along the chromosomes of the Congjiang and Jianbai Xiang pig breeds
2.4 从江香猪和剑白香猪分化水平

将从江香猪和剑白香猪基因组中每条染色体上的核苷酸多样性π值、群体间分歧度检验Fst、群体分歧度ROD以及群体内选择检验Tajima ’s D的计算结果取均值后绘制成图 4。结果表明，从江香猪和剑白香猪的核苷酸多样性均小于0.5%，且群体内选择检验Tajima ’s D值均大于0。此外，剑白香猪和从江香猪的ROD值为0.553 1，Fst值为0.736 2。

图 4 从江香猪和剑白香猪分化水平参数 Fig. 4 The differentiation level parameters of Congjiang and Jianbai Xiang pigs
2.5 剑白香猪群体受选择区域

对于受选择区域的筛选，一般考虑群体间分歧度高及群体分化率高的区域。以从江香猪为对照组，结合群体间分歧度ROD和群体间分歧度检验Fst的计算结果，筛选出剑白香猪基因组中前5%受到强烈选择作用、高度分化的区域，并从中识别出了共计775个受选择区域基因，其分布情况如图 5红色区域所示。

图 5 剑白香猪Fst-ROD的受选择区域 Fig. 5 The Fst-ROD regions under selection in Jianbai Xiang pigs
2.6 剑白香猪群体受选择区域基因富集分析

对受选择区域基因进行了KEGG富集分析(图 6A)，结果表明前20种KEGG富集通路涉及轴突引导、ECM与受体的相互作用、局部粘附、PI3K-Akt信号通路、硫辛酸代谢、cGMP-PKG信号通路、MAPK信号传导途径、核苷酸切除修复、阿佩林信号通路、心律失常性右心室心肌病等通路。

图 6 剑白香猪群体受选择区域基因KEGG(A)和GO分析(B)结果 Fig. 6 The KEGG (A) and GO analysis (B) results of the genes located in the selected regions in the population of Jianbai Xiang pigs

同时，对受选择区域基因进行了GO富集分析(图 6B)，发现在从江香猪和剑白香猪的775个候选基因中，有367个基因显著富集于751个GO条目(P < 0.05)，其中包括113个分子功能类，553个生物过程类，85个细胞组分类。在分子功能中，受选择区域基因主要富集在结合和催化活性；在生物过程中，受选择区域基因富集最集中的GO条目在细胞过程、单一生物体过程、生物调节、新陈代谢过程和生物过程的调节上；在细胞组分中，细胞、细胞部分、细胞器条目上的基因富集较多。

3 讨论

遗传分化由长期自然选择和人工选择造成，受环境、基因漂变、基因流动和遗传漂变等因素的影响。ROD基于两个群体间的核酸多态性π的差异识别信号，被引入用于测量两群体间的群体多态性损失，ROD平均值越接近1，群体多态性损失越大。本研究中，从江香猪和剑白香猪的ROD值为0.553 1，表示从江香猪和剑白香猪之间存在一定的群体多态性损失，群体分歧度较高。Fst值主要受遗传漂变和迁移等因素的影响，种群中对于特定环境经历的适应性进化会增大种群分化水平，当Fst值大于0.25时表示群体间遗传分化程度较大，Fst越接近1，表明两种群间越存在明显的种群分化[19]。尤桂爽等[20]对四川盆地6个地方鸡群体的遗传关系研究发现，丰岩乌骨鸡和岩水鸡黑羽群体之间具有较高的遗传分化(Fst=0.183)。本研究中，从江香猪和剑白香猪的Fst值为0.736 2，结果显示从江香猪种群和剑白香猪种群之间存在着明显的种群分化。群体内选择检验Tajima ’s D用于检验正选择效应[21]，本研究中从江香猪和剑白香猪的Tajima ’s D值均大于0，说明二者基因组内存在大量的中等频率等位基因，可能是由于群体瓶颈效应，群体结构或者平衡选择等原因所致。

群体有效大小越小意味着该群体的基因多样性越低，本研究发现，从江香猪和剑白香猪在13个世代之前的群体有效大小为18头，明显低于商品猪群体有效大小[22]，同时也低于目前报道的梅山猪(50~ 65头)[23]、陆川猪(35头)[24]和撒坝猪(21头)[25]，高于凉山猪(15头)[26]等中国其他地方猪种群体有效大小。多态位点比例与群体中等位基因的数目和频率有关，用来表示群体中某一位点多态性的程度，以评估群体的遗传多样性[27]，本研究发现从江香猪和剑白香猪群体的多态位点比例分别为0.875 7和0.813 2，稍低于西方猪种0.90~0.98的多态位点比例范围[28]，剑白香猪的多态位点比例稍低于里岔黑猪(0.827)[29]。说明从江香猪和剑白香猪相比其他群体遗传多样性较低，这对一个种群的长期生存和变化环境的适应是不利的。期望杂合度是指理论计算得出的杂合度，观测杂合度是指随机抽取的两个样本的等位基因不相同的概率。若观测杂合度高于期望杂合度，表明群体内存在较高的遗传多样性，群体可能受到迁徙或者基因流的影响[30]。本研究中，从江香猪和剑白香猪的期望杂合度均高于观测杂合度，说明二者群体可能出现过迁徙事件或者有外来血缘的引入，这与杭猪[31]、青峪猪[32]、马身猪[27]、通城猪[33]、丫杈猪[34]群体的研究结果相似。

虽然剑白香猪与从江香猪同为小型猪种，但本研究中的5个剑白香猪样品中有3个经过14代自交培育，相比之下，剑白香猪种群更为纯正，遗传背景更为稳定。根据KEGG分析结果，剑白香猪和从江香猪的受选择区域基因中差异表达的通路涉及了多个生物学过程和信号通路，其中包括涉及轴突引导、ECM与受体的相互作用、局部粘附、PI3K-Akt信号通路、硫辛酸代谢、cGMP-PKG信号通路、MAPK信号传导途径、核苷酸切除修复、阿佩林信号通路、心律失常性右室心肌病等通路。从这些通路中可以初步推断出剑白香猪在某些方面的优势。PI3K-Akt信号通路和MAPK信号通路等与生长和发育[35]、代谢调节[36]、细胞存活和凋亡[37]等生物学过程有关，因此剑白香猪在这些方面可能表现出更好的生物学特性。硫辛酸代谢也与肝脏解毒、生物体对药物的代谢、胆固醇代谢、维生素代谢相关[38]，因此剑白香猪在物质代谢方面也可能有优势。猪由于其在体型、解剖结构、生理、营养物质代谢、药物代谢和疾病发展方面与人类相似，所以被广泛作为人类疾病研究的生物模型[39]，Fanjul等[40]就使用小鼠和猪模型来研究衰老对心脏代谢疾病的影响，本研究KEGG分析结果中，剑白香猪群体受选择区域基因富集到心律失常性右室心肌病通路，可以考虑将剑白香猪作为研究心律失常性右心室心肌病的生物医学模型。根据GO富集分析结果，受选择区域基因富集在结合和催化活性条目，这说明剑白香猪在分子功能方面具有更多的活性相关基因，这可能意味着该品种在代谢、生长和发育等方面具有更强的能力；而在细胞过程、单一生物体过程、生物调节、新陈代谢过程和生物过程调节的富集，表明剑白香猪在生物过程和生物调节方面可能具有更为丰富的功能；细胞、细胞部分、细胞器条目的富集，表明剑白香猪在细胞结构和细胞器的功能方面具有更多的相关基因。

4 结论

综上所述，4个统计数据(Ne、PN、Ho、He)共同表明，需要提高从江香猪和剑白香猪群体内的遗传多样性，减少外来血缘引入导致纯种度降低的风险。群体全基因组选择信号分析表明，剑白香猪和从江香猪之间存在一定的群体多态性损失，群体分歧度较高，两种群间存在明显的种群分化，而二者群体受选择区域基因富集分析表明，剑白香猪在代谢、生长和发育等方面具有更强的能力。自交培育的剑白香猪可以用作研究心律失常性右室心肌病的动物模型。这些发现有助于更深入地了解贵州地方香猪品种的遗传特征和潜在用途，科学制定对贵州地方香猪的保护和利用策略，以保障其种质资源的可持续利用。

参考文献
[1]
AI H S, FANG X D, YANG B, et al. Adaptation and possible ancient interspecies introgression in pigs identified by whole-genome sequencing[J]. Nat Genet, 2015, 47(3): 217-225. DOI:10.1038/ng.3199
[2]
ZHAO Q B, LÓPEZ-CORTEGANO E, OYELAMI F O, et al. Conservation priorities analysis of Chinese indigenous pig breeds in the Taihu Lake region[J]. Front Genet, 2021, 12: 558873. DOI:10.3389/fgene.2021.558873
[3]
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(7): 657. DOI:10.1186/1297-9686-35-7-657
[4]
WU F, SUN H, LU S X, et al. Genetic diversity and selection signatures within Diannan small-ear pigs revealed by next-generation sequencing[J]. Front Genet, 2020, 11: 733. DOI:10.3389/fgene.2020.00733
[5]
赵真坚, 王书杰, 陈栋, 等. 基于低深度全基因组测序分析内江猪群体结构和遗传多样性[J]. 畜牧兽医学报, 2023, 54(6): 2297-2307.
ZHAO Z J, WANG S J, CHEN D, et al. Population structure and genetic diversity analysis of Neijiang pigs based on low-coverage whole genome sequencing[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(6): 2297-2307. (in Chinese)
[6]
MEI C G, WANG H C, LIAO Q J, et al. Genetic architecture and selection of Chinese cattle revealed by whole genome resequencing[J]. Mol Biol Evol, 2018, 35(3): 688-699. DOI:10.1093/molbev/msx322
[7]
SHI H B, LI T T, SU M C, et al. Whole genome sequencing revealed genetic diversity, population structure, and selective signature of Panou Tibetan sheep[J]. BMC Genomics, 2023, 24(1): 50. DOI:10.1186/s12864-023-09146-2
[8]
马克岩, 韩金涛, 白雅琴, 等. 基于简化基因组测序的永登七山羊遗传多样性分析[J]. 畜牧兽医学报, 2023, 54(5): 1939-1950.
MA K Y, HAN J T, BAI Y Q, et al. Genetic diversity analysis of Yongdeng Qishan sheep based on specific-locus amplified fragment sequencing[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(5): 1939-1950. (in Chinese)
[9]
SUN J L, CHEN T, ZHU M, et al. Whole-genome sequencing revealed genetic diversity and selection of Guangxi indigenous chickens[J]. PLoS One, 2022, 17(3): e0250392. DOI:10.1371/journal.pone.0250392
[10]
刘天飞, 罗成龙, 王艳, 等. 基于简化基因组测序技术和基因芯片技术比较研究黄羽肉鸡基因组选择[J]. 畜牧兽医学报, 2020, 51(10): 2378-2386.
LIU T F, LUO C L, WANG Y, et al. Comparison between reduced-representation genome sequencing and SNP chip for genomic selection in yellow-feathered broiler[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(10): 2378-2386. DOI:10.11843/j.issn.0366-6964.2020.10.006 (in Chinese)
[11]
刘宏祥, 沈永杰, 张丽华, 等. 基于简化基因组测序的娄门鸭遗传多样性评价[J]. 畜牧兽医学报, 2022, 53(6): 1735-1748.
LIU H X, SHEN Y J, ZHANG L H, et al. Genetic diversity evaluation of Loumen duck based on reduced-representation genome sequencing[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(6): 1735-1748. (in Chinese)
[12]
张福平, 唐靓婷, 王嘉福, 等. 基于转录组测序筛选影响从江香猪产仔数的候选基因[J]. 南方农业学报, 2021, 52(4): 847-856.
ZHANG F P, TANG L T, WANG J F, et al. Selection of candidate genes affecting litter size of Congjiang Xiang pig by transcriptome sequencing[J]. Journal of Southern Agriculture, 2021, 52(4): 847-856. (in Chinese)
[13]
牟琪, 冯贤辀, 蔡栋, 等. 去势时间对从江香猪肉质性状的影响[J]. 南方农业学报, 2022, 53(11): 3276-3284.
MOU Q, FENG X Z, CAI D, et al. Effects of castration time on meat quality of Congjiang Xiang pig[J]. Journal of Southern Agriculture, 2022, 53(11): 3276-3284. DOI:10.3969/j.issn.2095-1191.2022.11.029 (in Chinese)
[14]
张依裕, 刘若余, 刘培琼, 等. 剑白香猪生长激素基因多态性与部分生产性能的关联分析[J]. 黑龙江畜牧兽医, 2009(7): 107-109.
ZHANG Y Y, LIU R Y, LIU P Q, et al. Association analysis of growth hormone gene polymorphism and partial production performance in Jianbai Xiang pigs[J]. Heilongjiang Animal Science and Veterinary Medicine, 2009(7): 107-109. (in Chinese)
[15]
CHEN S F, ZHOU Y Q, CHEN Y R, et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor[J]. Bioinformatics, 2018, 34(17): i884-i890. DOI:10.1093/bioinformatics/bty560
[16]
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.
[17]
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. DOI:10.1086/519795
[18]
WU T Z, HU E Q, XU S B, et al. clusterProfiler 4.0:a universal enrichment tool for interpreting omics data[J]. Innovation (Camb), 2021, 2(3): 100141.
[19]
OCHOA A, STOREY J D. Estimating FST and kinship for arbitrary population structures[J]. PLoS Genet, 2021, 17(1): e1009241. DOI:10.1371/journal.pgen.1009241
[20]
尤桂爽, 易兴友, 蒋美山, 等. 四川盆地6个地方鸡群体遗传关系研究[J]. 畜牧兽医学报, 2022, 53(2): 370-380.
YOU G S, YI X Y, JIANG M S, et al. Study on the genetic relationship of six local chicken populations from Sichuan Basin[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(2): 370-380. (in Chinese)
[21]
TAJIMA F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism[J]. Genetics, 1989, 123(3): 585-595. DOI:10.1093/genetics/123.3.585
[22]
TRASPOV A, DENG W J, KOSTYUNINA O, et al. Erratum to: population structure and genome characterization of local pig breeds in Russia, Belorussia, Kazakhstan and Ukraine[J]. Genet Sel Evol, 2016, 48(1): 57. DOI:10.1186/s12711-016-0235-8
[23]
孙浩, 王振, 张哲, 等. 基于基因组测序数据的梅山猪保种现状分析[J]. 上海交通大学学报: 农业科学版, 2017, 35(4): 65-70.
SUN H, WANG Z, ZHANG Z, et al. Exploring the current situation of conservation of Meishan pigs based on genome sequencing data[J]. Journal of Shanghai Jiaotong University: Agricultural Science, 2017, 35(4): 65-70. (in Chinese)
[24]
DIAO S Q, HUANG S W, XU Z T, et al. Genetic diversity of indigenous pigs from South China area revealed by SNP array[J]. Animals (Basel), 2019, 9(6): 361.
[25]
邓俊, 刘艺端, 许文坤, 等. 基于SNP芯片撒坝猪保种群体的遗传结构分析[J]. 中国饲料, 2022(17): 7-11.
DENG J, LIU Y R, XU W K, et al. Genetic structure analysis of a conserved population of Saba pigs based on SNP chips[J]. China Feed, 2022(17): 7-11. (in Chinese)
[26]
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. DOI:10.5713/ajas.19.0884
[27]
蔡春波, 张雪莲, 张万峰, 等. 运用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)
[28]
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. DOI:10.1371/journal.pone.0056001
[29]
WANG Y, DONG R L, LI X, et al. Analysis of the genetic diversity and family structure of the Licha black pig population on Jiaodong Peninsula, Shandong Province, China[J]. Animals (Basel), 2022, 12(8): 1045.
[30]
EDEA Z, KIM S W, LEE K T, et al. Genetic structure of and evidence for admixture between western and Korean native pig breeds revealed by single nucleotide polymorphisms[J]. Asian-Australas J Anim Sci, 2014, 27(9): 1263-1269.
[31]
刘晨龙, 卢丹, 周泉勇, 等. 利用高密度SNP芯片分析杭猪的群体遗传结构[J]. 畜牧兽医学报, 2022, 53(8): 2502-2513.
LIU C L, LU Q, ZHOU Q Y, et al. Analysis of population genetic structure of Hang pigs by high density SNP chip[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(8): 2502-2513. (in Chinese)
[32]
刘彬, 沈林園, 陈映, 等. 基于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)
[33]
袁娇, 徐国强, 周翔, 等. 基于SNP芯片监测通城猪的保种效果[J]. 畜牧兽医学报, 2022, 53(8): 2514-2523.
YUAN J, XU G Q, ZHOU X, et al. SNP chip-based monitoring of population conservation effect of Tongcheng pigs[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(8): 2514-2523. (in Chinese)
[34]
陶璇, 杨雪梅, 梁艳, 等. 基于SNP芯片的丫杈猪保种群体遗传结构研究[J]. 畜牧兽医学报, 2023, 54(6): 2308-2319.
TAO X, YANG X M, LIANG Y, et al. Analysis of genetic structure of conservation population in Yacha pig based on SNP chip[J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(6): 2308-2319. (in Chinese)
[35]
CAI P R, ZHU Q H, CAO Q Y, et al. Quercetin and allicin can alleviate the hepatotoxicity of lead (Pb) through the PI3K signaling pathway[J]. J Agric Food Chem, 2021, 69(32): 9451-9460.
[36]
XU W J, LU Z X, WANG X, et al. Gynura divaricata exerts hypoglycemic effects by regulating the PI3K/AKT signaling pathway and fatty acid metabolism signaling pathway[J]. Nutr Diabetes, 2020, 10(1): 31.
[37]
DONG W W, LUO B, QIU C, et al. TRIM3 attenuates apoptosis in Parkinson 's disease via activating PI3K/AKT signal pathway[J]. Aging (Albany NY), 2020, 13(1): 735-749.
[38]
NAJAFI N, MEHRI S, GHASEMZADEH RAHBARDAR M, et al. Effects of alpha lipoic acid on metabolic syndrome: a comprehensive review[J]. Phytother Res, 2022, 36(6): 2300-2323.
[39]
WEI J L, ZHANG W, LI J, et al. Application of the transgenic pig model in biomedical research: a review[J]. Front Cell Dev Biol, 2022, 10: 1031812.
[40]
FANJUL V, JORGE I, CAMAFEITA E, et al. Identification of common cardiometabolic alterations and deregulated pathways in mouse and pig models of aging[J]. Aging Cell, 2020, 19(9): e13203.

(编辑   郭云雁)