基于重测序的民猪与大白猪SLAⅠ类基因多态性及抗病潜能差异分析

doi:10.11843/j.issn.0366-6964.2021.07.006

摘要/Abstract

摘要： 不同中外猪种对疾病的抵抗能力有一定差异，原因之一是不同猪种猪白细胞抗原（swine leukocyte antigen，SLA）分子的多样性不同。本研究旨在探究中外猪种的SLA基因差异，为阐明地方猪的抗病机理提供重要参考。本研究首先对健康的8头民猪以及4头大白猪进行基因组重测序，并对所获得的SNP进一步质控用于后续分析。使用VCFtools、GEVALT软件分析SLAⅠ类基因的SNP以及单倍型，并将SNP利用ENSEMBL中的VEP工具进行注释，在全局层面阐述两个猪种的SLAⅠ类基因多样性。通过Mega软件比对两猪种的经典SLAⅠ类基因外显子2和3的核苷酸及编码氨基酸序列，利用Expasy服务器上的ProtParam工具和Protscale程序分析蛋白质特性，并利用DnaSP软件计算核苷酸多样性，分析两个猪种经典SLAⅠ类基因抗原递呈能力的差异。结果表明，民猪的SLAⅠ类基因具有更多的SNPs和长度较短的单倍型块，并且错义突变的碱基数量较多。在经典SLAⅠ类基因的外显子2和3可分别鉴定出2个等位基因，民猪在等位基因上的核苷酸多样性要明显高于大白猪，并且民猪具有更多的碱基突变以及位于抗原结合位点上的氨基酸突变。两个猪种的经典SLAⅠ类基因氨基酸表现为亲水性，民猪的亲水性要强于大白猪。综上所述，SLAⅠ类基因在民猪上有更强的多态性。民猪经典SLAⅠ类基因的碱基突变数量和抗原结合位点上氨基酸突变的数量都要多于大白猪。本研究所检测到的民猪在经典SLAⅠ类基因ARSs上的变异可为猪抗病育种标记的筛选提供参考。

关键词: 民猪, 大白猪, SLAⅠ类基因, 抗原结合位点, 抗病力

Abstract: There are certain differences in the resistance to diseases of different pig breeds, one of the possible reasons for this difference is the diversity of swine leukocyte antigen (SLA) between different pig breeds. The study of SLA gene differences between Chinese indigenous and foreign pig breeds can provide important reference for clarifying the disease resistance mechanism of Chinese indigenous breeds. In this study, 8 healthy Min pigs and 4 Large White pigs were resequenced, and the obtained SNPs were then quality controlled for subsequent analysis. VCFtools and GEVALT softwares were used to analyze the SNP and haplotypes of SLA class I genes, and VEP tool in ENSEMBL was used to annotate SNPs, thus we can expound the diversity of SLA class I genes in two breeds at the global level. Mega software was applied to compare the sequences of nucleotide and coded amino acid of exon 2 and 3 of SLA class I classic genes of the two pig breeds. ProtParam tool and Protscale program on the Expasy server were performed to analyze the protein characteristics, and DnaSP software was used to calculate nucleotide diversity, the differences in antigen presentation ability of SLA class I classic genes between Chinese and foreign pig breeds were analyzed. The results showed that SLA class I genes of Min pigs had more SNPs and shorter haplotype blocks, and the number of missense mutations was more in Min pigs. Two alleles could be identified in exon 2 and 3 of SLA class I classic genes. The nucleotide diversity of alleles in Min pigs was higher than that in Large White pigs in the SLA class I classic genes. Moreover, Min pigs had more nucleotide and amino acid mutations at the antigen-binding site of SLA class I classic genes. The amino acids coded by SLA class Ⅰ classic genes of two pig breeds were hydrophilic, and it was more hydrophilic in Min pigs than that in Large White pigs. In summary, SLA class I genes have more polymorphisms in Min pigs. The number of base and amino acid mutations in the antigen-binding site of SLA class I classic genes of Min pigs are more than those of Large White pigs. The mutations in the ARSs of SLA class I classic gene of Min pigs detected in this study can provide a reference for the selection of pig disease resistance breeding markers.

Key words: Min pigs, Large White pigs, SLA class I genes, antigen-binding site, disease resistance

中图分类号:

S828.2

李婧娜, 疏泽, 吴俊超, 程锋, 王立刚, 王立贤. 基于重测序的民猪与大白猪SLAⅠ类基因多态性及抗病潜能差异分析[J]. 畜牧兽医学报, 2021, 52(7): 1820-1830.

LI Jingna, SHU Ze, WU Junchao, CHENG Feng, WANG Ligang, WANG Lixian. Analysis of SLA Class Ⅰ Gene Polymorphism and Potential Disease Resistance Difference between Min Pigs and Large White Pigs Based on Resequencing[J]. Acta Veterinaria et Zootechnica Sinica, 2021, 52(7): 1820-1830.

参考文献 42

[1]	WANG J Y,WANG Y P,WANG H Z,et al.Transcriptomic analysis identifies candidate genes and gene sets controlling the response of porcine peripheral blood mononuclear cells to Poly I:C Stimulation[J].G3(Bethesda),2016,6(5):1267-1275.
[2]	LUNNEY J K,CHEN H B.Genetic control of host resistance to porcine reproductive and respiratory syndrome virus (PRRSV) infection[J].Virus Res,2010,154(1-2):161-169.
[3]	YANG H.Livestock development in China:animal production,consumption and genetic resources[J].J Anim Breed Genet,2013,130(4):249-251.
[4]	LE M T,CHOI H,LEE H,et al.SLA-1 genetic diversity in pigs:extensive analysis of copy number variation,heterozygosity, expression,and breed specificity[J].Sci Rep,2020,10(1):743.
[5]	奚衍洋,赵云蛟,钱爱东.猪SLA免疫学相关研究进展[J].黑龙江畜牧兽医,2018(3):70-72.XI Y Y,ZHAO Y J,QIAN A D.Immunological progress of swine major histocompatibility complex[J]. Heilongjiang Animal Science and Veterinary Medicine,2018(3):70-72.(in Chinese)
[6]	PARK J Y,KIM H S,SEO S H.Characterization of interaction between porcine reproductive and respiratory syndrome virus and porcine dendritic cells[J].J Microbiol Biotechnol,2008,18(10):1709-1716.
[7]	WANG X,EATON M,MAYER M,et al.Porcine reproductive and respiratory syndrome virus productively infects monocyte-derived dendritic cells and compromises their antigen-presenting ability[J].Arch Virol,2007, 152(2):289-303.
[8]	訾臣,刘璐,包文斌.SLA-DQA基因在F18大肠杆菌抗性和敏感性断奶仔猪间的差异表达分析[J].畜牧兽医学报, 2012,43(4):516-520.ZI C,LIU L,BAO W B.The differential expression of SLA-DQA gene between weaning piglets resistant or sensitive to Escherichia coli F18[J].Acta Veterinaria et Zootechnica Sinica,2012,43(4):516-520.(in Chinese)
[9]	包文斌,吴圣龙,曹晶晶,等.19个猪种SLA-DQB基因外显子2多态性分析[J].畜牧兽医学报, 2009,40(10):1550-1554.BAO W B,WU S L,CAO J J,et al.Polymorphism analysis on exon 2 of SLA-DQB gene in 19 pig populations[J].Acta Veterinaria et Zootechnica Sinica,2009,40(10):1550-1554.(in Chinese)
[10]	HAMMER S E,HO C S,ANDO A,et al.Importance of the Major histocompatibility complex (swine leukocyte antigen) in swine health and biomedical research[J].Annu Rev Anim Biosci,2020,8:171-198.
[11]	RENARD C,VAIMAN M,CHIANNILKULCHAI N,et al.Sequence of the pig major histocompatibility region containing the classical class I genes[J].Immunogenetics,2001,53(6):490-500.
[12]	CHARDON P,RENARD C,VAIMAN M.The major histocompatibility complex in swine[J].Immunol Rev,1999,167(1):179-192.
[13]	FROMMER L,FLESCH B K,KNIG J,et al.Amino acid polymorphisms in Hla class Ⅱ differentiate between thyroid and polyglandular autoimmunity[J].J Cli Endocrinol Metab,2020,105(6):1737-1747.
[14]	韦思婷,程一飞,刘丽霞.猪SLA基因多态性研究进展[J].猪业科学,2019,36(4):102-103.WEI S T,CHENG Y F,LIU L X.Research progress of porcine SLA gene polymorphism[J].Swine Industry Science,2019, 36(4):102-103.(in Chinese)
[15]	HUANG Z Y,YU Y Y,DONG S,et al.Major histocompatibility complex class IIA and IIB genes of loach (Misgurnus anguillicaudatus):Molecular cloning and expression analysis in response to bacterial and parasitic challenge[J]. Aquaculture,2019,500:359-369.
[16]	CORTÁZAR-CHINARRO M,MEYER-LUCHT Y,LAURILA A,et al.Signatures of historical selection on MHC reveal different selection patterns in the moor frog (Rana arvalis)[J].Immunogenetics, 2018, 70(7):477-484.
[17]	KIM T,HUNT H D,PARCELLS M S,et al.Two class I genes of the chicken MHC have different functions:BF1 is recognized by NK cells while BF2 is recognized by CTLs[J].Immunogenetics,2018,70(9):599-611.
[18]	MORGAN R A,KARL J A,BUSSAN H E,et al.Restricted MHC class I A locus diversity in olive and hybrid olive/yellow baboons from the Southwest National Primate Research Center[J].Immunogenetics, 2018,70(7):449-458.
[19]	YU S,WU J,BAI J,et al.Polymorphic analysis of peptide binding domain of major histocompatibility complex class I in domestic ducks[J].Pol J Vet Sci,2019,22(2):415-422.
[20]	LI Z B,ZHANG N,MA L Z,et al.Distribution of ancient α1 and α2 domain lineages between two classical MHC class I genes and their alleles in grass carp[J].Immunogenetics,2019,71(5-6):395-405.
[21]	KIEMNEC-TYBURCZY K M,TRACY K E,LIPS K R,et al.Genetic variation and selection of MHC class I loci differ in two congeneric frogs[J].Genetica,2018,146(2):125-136.
[22]	SAVAGE A E,MULDER K P,TORRES T,et al.Lost but not forgotten:MHC genotypes predict overwinter survival despite depauperate MHC diversity in a declining frog[J].Conservat Genet,2018,19(2):309-322.
[23]	陶久臣.东北林蛙MHC基因的多态性及其遗传进化[D].哈尔滨:哈尔滨师范大学,2019.TAO J C.Polmorphism and genetic evolution of MHC class genes in Rana Dybowskii[D].Harbin:Harbin Normal University,2019.(in Chinese)
[24]	LIANG W,LI Z H,WANG P,et al.Differences of immune responses between Tongcheng (Chinese local breed) and Large White pigs after artificial infection with highly pathogenic porcine reproductive and respiratory syndrome virus[J].Virus Res,2016,215:84-93.
[25]	XING J Y,XING F,ZHANG C H,et al.Genome-wide gene expression profiles in lung tissues of pig breeds differing in resistance to porcine reproductive and respiratory syndrome virus[J].PLoS One,2014,9(1):e86101.
[26]	JIANG C L,XING F,XING J Y,et al.Different expression patterns of PRRSV mediator genes in the lung tissues of PRRSV resistant and susceptible pigs[J].Dev Comp Immunol,2013,39(1-2):127-131.
[27]	莫斯科.中国合作猪SLAⅠ类和Ⅱ类分子的克隆及基因结构与多态性分析[D].兰州:甘肃农业大学,2009.MO S K.Cloning and analysis of genomic structure and genetic polymorphism of swine lymphocyte antigen class Ⅰ and Ⅱ molecules in Chinese Hezuo pig[D].Lanzhou:Gansu Agricultural University,2009.(in Chinese)
[28]	DAVIDOVICH O,KIMMEL G,SHAMIR R.GEVALT:An integrated software tool for genotype analysis[J].BMC Bioinform,2007, 8(1):36.
[29]	KIMMEL G,SHAMIR R.GERBIL:Genotype resolution and block identification using likelihood[J].Proc Natl Acad Sci U S A,2005, 102(1):158-162.
[30]	FILIPSKI A,MURILLO O,FREYDENZON A,et al.Prospects for building large timetrees using molecular data with incomplete gene coverage among species[J].Mol Biol Evol,2014,31(9):2542-2550.
[31]	ROZAS J,FERRER-MATA A,SÁNCHEZ-DELBARRIO J C,et al.DnaSP 6:DNA sequence polymorphism analysis of large data sets[J].Mol Biol Evol,2017,34(12):3299-3302.
[32]	张冬杰,刘娣.民猪种质资源特点及研究现状[J].黑龙江农业科学,2019(2):48-50.ZHANG D J,LIU D.Germplasm resources characteristics and research status of min pig[J].Heilongjiang Agricultural Sciences,2019(2):48-50.(in Chinese)
[33]	FENG S H,FANG Q,BARNETT R,et al.The genomic footprints of the fall and recovery of the crested ibis[J].Curr Biol,2019,29(2):340-349.e7.
[34]	KOSCH T A,EIMES J A,DIDINGER C,et al.Characterization of MHC class IA in the endangered southern corroboree frog[J].Immunogenetics,2017,69(3):165-174.
[35]	边若菲,徐笑,刘玉芬,等.嗜水气单胞菌胁迫下东北林蛙MHCⅠ基因在不同组织的表达[J].生物工程学报,2020,36(7):1323-1333.BIAN R F,XU X,LIU Y F,et al.Expression of MHC Ⅰ genes in different tissues of Rana dybowskii under the stress of Aeromonas hydrophila[J].Chinese Journal of Biotechnology,2020,36(7):1323-1333.(in Chinese)
[36]	刘建欣,郑昌学.现代免疫学:免疫的细胞和分子基础[M].北京:清华大学出版社,2002.LIU J X,ZHENG C X.Modern immunology:the cellular and molecular basis of immunity[M].Beijing:Tsinghua University Press,2002.(in Chinese)
[37]	ANDO A,KAWATA H,SHIGENARI A,et al.Genetic polymorphism of the swine major histocompatibility complex (SLA) class I genes,SLA-1,-2 and -3[J].Immunogenetics,2003,55(9):583-593.
[38]	陶璇.SLAⅠ类经典基因多态性及正选择检验分析[D].成都:四川农业大学,2008.TAO X.Polymorphisms Analysis and detecting positive selection of SLA classic class I genes[D].Chengdu:Sichuan Agricultural University,2008.(in Chinese)
[39]	DUNNER S,CHECA M L,GUTIERREZ J P,et al.Genetic analysis and management in small populations:the Asturcon pony as an example[J].Genet Sel Evol,1998,30(4):397-405.
[40]	刘德武,杨关福,李加琪,等.用RAPD标记分析6个品种猪的群体遗传结构[J].畜牧兽医学报, 2002,33(1):18-22.LIU D W,YANG G F,LI J Q,et al.Analysing genetic construction of six pig breeds using RAPD marker[J].Acta Veterinaria et Zootechnica Sinica,2002,33(1):18-22.(in Chinese)
[41]	张国华.合作猪SLA-DQA第2、4外显子多态性分析[D].兰州:甘肃农业大学,2014.ZHANG G H.Analyzing polymorphism of Hezuo swine SLA-DQA gene exon 2, 4[D].Lanzhou:Gansu Agricultural University,2014.(in Chinese)
[42]	杨汉春.动物免疫学[M].2版.北京:中国农业大学出版社,2003.YANG H C.Animal immunology[M].2nd ed.Beijing:China Agricultural University Press,2003.(in Chinese)