猪产仔数的全基因组关联分析

doi:10.7668/hbnxb.20192808

摘要/Abstract

摘要：

为筛选猪总产仔数和产活仔数性状关键SNP分子标记或候选基因,采用Illumina porcine 50K芯片对同一猪场186头加系大白猪能繁母猪开展了全基因组SNP扫描,后又对扫描结果进行了质量控制、Beagle填充和SNP基因型分型,结合这些试验猪的表型性状测定和群体结构分析结果,进行全基因组关联分析(GWAS)。群体结构分析结果表明,试验所用的加系大白猪样本未出现群体分层现象;在18对(36条)常染色体上一共得到36 867个有效SNP标记,用于试验猪的GWAS分析。GWAS结果表明,这些加系大白猪所有胎次总产仔数共显著关联到2个SNP,分别为seq-rs323899658和seq-rs329781338,其中seq-rs329781338 SNP又注释到3个基因,分别为SSBP1、WEE2和KIAA1147。产活仔数共显著关联到7个SNP,分别为seq-rs81238474、seq-rs321377412、seq-rs340736313、seq-rs80782154、seq-rs81244816、seq-rs81467772和seq-rs329781338,这7个SNP中有5个SNP有基因注释,共注释到22个基因,分别为EphA1、TAS2R60、FAM131B、CLCN1、CASP2、TMEM139、TRBV21OR9-2、PRSS2、TRBV19、TRBV24-1、U6、SSBP1、WEE2、KIAA1147、NKX2-8、ALKBH3、HSD17B12、U6、HOMER2、WHAMM、FSD2和SCARNA15。候选基因GO功能和KEGG通路分析结果表明,这些基因GO功能中最显著的生物学过程为多生物过程的调控、分子功能为非跨膜/跨膜蛋白酪氨酸激酶活性,KEGG通路为错配修复。结合GWAS和基因功能注释结果,WEE2和EphA1基因可作为提高加系大白猪总产仔数和产活仔数的关键候选基因。

关键词: 猪, 全基因组关联分析, 产仔数, SNP, 芯片

Abstract:

The objective of the study was to screen the candidate genes and SNP markers for the total number born(TNB)and number born alive(NBA)in pigs,so as to prepare for the next step of improving the reproduction traits,and breeding new high-proliferation of Canadian line of large white pigs. Illumina porcine 50K chip was used to do a genome-wide SNP scan on 186 Canadian line of large white sows in one pig farm,followed by quality control,beagle fill and SNP genotyping. Combined with the phenotypic traits and population structure analysis results of these experimental pigs,genome-wide association study(GWAS)was carried out. The results of population structure analysis showed that there was no population stratification in the pig samples used in the experiment. A total of 36 867 effective SNP markers were obtained on 18 pairs(36)of autosomes for GWAS analysis of experimental pigs. The results of GWAS showed that NBA in all parities of these Canadian large white pigs was significantly associated with two SNPs:seq-rs323899658 and seq-rs329781338,respectively. The SNP of seq-rs329781338 was annotated with three genes:SSBP1,WEE2 and KIAA1147,respectively. NBA was significantly associated with 7 SNPs,namely:seq-rs81238474,seq-rs321377412,seq-rs340736313,seq-rs80782154,seq-rs81244816,seq-rs81467772 and seq-rs329781338. Five out of seven SNPs had gene annotation,and a total of 22 genes were annotated:EphA1, TAS2R60, FAM131B, CLCN1, CASP2, TMEM139, TRBV21OR9-2, PRSS2, TRBV19, TRBV24-1, U6, SSBP1, WEE2, KIAA1147, NKX2-8, ALKBH3, HSD17B12, U6, HOMER2, WHAMM, FSD2 and SCARNA15. The results of GO function and KEGG pathway analysis showed that the most significant biological process of GO function of these genes was the regulation of multi biological process,the molecular function was non transmembrane/transmembrane protein tyrosine kinase activity,and the KEGG pathway was mismatch repair.Combined with GWAS and gene function annotation results,WEE2 and EphA1 genes could be used as key candidate genes to increase the TNB and NBA of Canadian large white pigs,and provide some reference for the precise cultivation of high-yield pigs in the future.

Key words: Pig, GWAS, Litter size, SNP, Chip

付言峰, 赵为民, 李辉, 李碧侠, 程金花, 徐小波, 任守文. 猪产仔数的全基因组关联分析[J]. 华北农学报, 2022, 37(4): 232-238. doi: 10.7668/hbnxb.20192808.

FU Yanfeng, ZHAO Weimin, LI Hui, LI Bixia, CHENG Jinhua, XU Xiaobo, REN Shouwen. Genome Wide Association Analysis of Litter Size in Pigs[J]. Acta Agriculturae Boreali-Sinica, 2022, 37(4): 232-238. doi: 10.7668/hbnxb.20192808.

http://www.hbnxb.net/CN/Y2022/V37/I4/232

图/表 7

参考文献 23

[1]	付言峰, 方晓敏, 李碧侠, 赵为民, 王学敏, 陈哲, 任守文. 猪胚胎附植部分因子研究进展[J]. 中国畜牧兽医, 2019, 46(9):2676-2682. doi: 10.16431/j.cnki.1671-7236.2019.09.022. doi: 10.16431/j.cnki.1671-7236.2019.09.022
	Fu Y F, Fang X M, Li B X, Zhao W M, Wang X M, Chen Z, Ren S W. A review of some factors related to embryo implantation in pigs[J]. China Animal Husbandry & Veterinary Medicine, 2019, 46(9):2676-2682.
[2]	Ju M M, Wang X N, Li X J, Zhang M H, Shi L D, Hu P Y, Zhang B, Han X L, Wang K J, Li X L, Zhou L S, Qiao R M. Effects of litter size and parity on farrowing duration of Landrace Yorkshire sows[J]. Animals, 2021, 12(1):94. doi: 10.3390/ani12010094. doi: 10.3390/ani12010094 URL
[3]	周梅, 张玉, 王重龙. 猪产仔数性状候选基因研究进展[J]. 家畜生态学报, 2020, 41(8):1-6. doi: 10.3969/j.issn.1673-1182.2020.08.001. doi: 10.3969/j.issn.1673-1182.2020.08.001
	Zhou M, Zhang Y, Wang C L. Research progress on candidate genes for litter size traits in pigs[J]. Journal of Domestic Animal Ecology, 2020, 41(8):1-6.
[4]	Kwon W S, Rahman M S, Ryu D Y, Khatun A, Pang M G. Comparison of markers predicting litter size in different pig breeds[J]. Andrology, 2017, 5(3):568-577. doi: 10.1111/andr.12332. doi: 10.1111/andr.12332 pmid: 28409901
[5]	Ruan D L, Zhuang Z W, Ding R R, Qiu Y B, Zhou S P, Wu J, Xu C N, Hong L J, Huang S X, Zheng E Q, Cai G Y, Wu Z F, Yang J. Weighted single-step GWAS identified candidate genes associated with growth traits in a duroc pig population[J]. Genes, 2021, 12(1):117. doi: 10.3390/genes12010117. doi: 10.3390/genes12010117 URL
[6]	Montesinos-López O A, Montesinos-López A, Pérez-Rodríguez P, Barrón-López J A, Martini J W R, Fajardo-Flores S B, Gaytan-Lugo L S, Santana-Mancilla P C, Crossa J. A review of deep learning applications for genomic selection[J]. BMC Genomics, 2021, 22(1):19. doi: 10.1186/s12864-020-07319-x. doi: 10.1186/s12864-020-07319-x pmid: 33407114
[7]	谢健, 喻昌燕, 许瑶, 牛熙, 冉雪琴, 王嘉福. 可乐猪乳头数性状全基因组关联分析[J]. 中国畜牧兽医, 2016, 43(11):2866-2872. doi: 10.16431/j.cnki.1671-7236.2016.11.009. doi: 10.16431/j.cnki.1671-7236.2016.11.009
	Xie J, Yu C Y, Xu Y, Niu X, Ran X Q, Wang J F. Genome-wide association study of teats number traits in Kele pigs[J]. China Animal Husbandry & Veterinary Medicine, 2016, 43(11):2866-2872.
[8]	阮鹏, 黄敏婕, 郭晓令, 徐宁迎. 浙江省3个地方猪种繁殖性状的全基因组关联分析[J]. 中国畜牧杂志, 2018, 54(7):35-40. doi: 10.19556/j.0258-7033.2018-07-035. doi: 10.19556/j.0258-7033.2018-07-035
	Ruan P, Huang M J, Guo X L, Xu N Y. Genome-wide association study of reproductive traits in 3 Zhejiang indigenous pig breeds[J]. Chinese Journal of Animal Science, 2018, 54(7):35-40.
[9]	袁晓龙, 徐丹同, 曹浩男, 熊浩铭, 曾检华, 陈子韬, 张哲, 李加琪. 杜洛克母猪初情期日龄全基因组关联分析[J]. 广东农业科学, 2020, 47(5):87-92. doi: 10.16768/j.issn.1004-874X.2020.05.012. doi: 10.16768/j.issn.1004-874X.2020.05.012
	Yuan X L, Xu D T, Cao H N, Xiong H M, Zeng J H, Chen Z T, Zhang Z, Li J Q. A genome-wide association study for age at puberty of Duroc pigs[J]. Guangdong Agricultural Sciences, 2020, 47(5):87-92.
[10]	Zhou X, Stephens M. Genome-wide efficient mixed-model analysis for association studies[J]. Nature Genetics, 2012, 44(7):821-824. doi: 10.1038/ng.2310. doi: 10.1038/ng.2310 URL
[11]	Turner S D. Qqman:an R package for visualizing GWAS results using Q-Q and manhattan plots[J]. Biorxiv, 2014, 5(14):1-2. doi: 10.1101/005165. doi: 10.1101/005165
[12]	王成斌, 王小鹏, 郑小明, 陈浩, 张俊杰, 郭源梅, 丁能水. 玉山黑猪两种头型遗传机制的初步解析[J]. 畜牧兽医学报, 2018, 49(8):1585-1593. doi: 10.11843/j.issn.0366-6964.2018.08.003. doi: 10.11843/j.issn.0366-6964.2018.08.003
	Wang C B, Wang X P, Zheng X M, Chen H, Zhang J J, Guo Y M, Ding N S. A study of the genetic mechanisms for two head types in Yushan black pigs[J]. Acta Veterinaria et Zootechnica Sinica, 2018, 49(8):1585-1593.
[13]	Akond Z, Ahsan M A, Alam M, Mollah M N H. Robustification of GWAS to explore effective SNPs addressing the challenges of hidden population stratification and polygenic effects[J]. Scientific Reports, 2021, 11(1):13060. doi: 10.1038/s41598-021-90774-7. doi: 10.1038/s41598-021-90774-7 URL
[14]	Oh J S, Susor A, Conti M. Protein tyrosine kinase Wee1B is essential for metaphase Ⅱ exit in mouse oocytes[J]. Science, 2011, 332(6028):462-465. doi: 10.1126/science.1199211. doi: 10.1126/science.1199211 URL
[15]	刘顿, 陈创奇, 苏念军, 何天文, 燕瑾, 刘风华, 张曦倩. WEE2基因纯合突变Lys195Asn导致受精失败和女性不孕的家系分析[J]. 分子诊断与治疗杂志, 2021, 13(7):1035-1038. doi: 10.19930/j.cnki.jmdt.2021.07.004. doi: 10.19930/j.cnki.jmdt.2021.07.004
	Liu D, Chen C Q, Su N J, He T W, Yan J, Liu F H, Zhang X Q. Lys195Asn homozygous mutation of WEE2 gene causes fertilization failure and a pedigree analysis of female infertility[J]. Journal of Molecular Diagnostics and Therapy, 2021, 13(7):1035-1038.
[16]	Wang Y, Hu L, Zhang X M, Zhao H, Xu H, Wei Y H, Jiang H G, Xie C H, Zhou Y F, Zhou F X. Down regulation of mitochondrial single stranded DNA binding protein(SSBP1)induces mitochondrial dysfunction and increases the radiosensitivity in non-small cell lung cancer cells[J]. Journal of Cancer, 2017, 8(8):1400-1409. doi: 10.7150/jca.18170. doi: 10.7150/jca.18170 pmid: 28638454
[17]	苗丽, 王栋, 马莉, 康安静, 张洁. 卵巢癌组织中SSBP1的表达意义及作用[J]. 实用癌症杂志, 2021, 36(2):195-198. doi: 10.3969/j.issn.1001-5930.2021.02.005. doi: 10.3969/j.issn.1001-5930.2021.02.005
	Miao L, Wang D, Ma L, Kang A J, Zhang J. Expression and significance of SSBP1 in ovarian cancer[J]. The Practical Journal of Cancer, 2021, 36(2):195-198.
[18]	Fujii H, Fujiwara H, Horie A, Sato Y, Konishi I. Ephrin A1 induces intercellular dissociation in Ishikawa cells:Possible implication of the Eph-ephrin A system in human embryo implantation[J]. Human Reproduction, 2010, 26(2):299-306. doi: 10.1093/humrep/deq340. doi: 10.1093/humrep/deq340 URL
[19]	Lim W, Bae H, Bazer F W, Song G. Functional roles of eph A-ephrin A1 system in endometrial luminal epithelial cells during early pregnancy[J]. Journal of Cellular Physiology, 2017, 232(6):1527-1538. doi: 10.1002/jcp.25659. doi: 10.1002/jcp.25659 URL
[20]	田明明, 王爱国, 傅金銮. 猪EPHA1基因多态性及其与产仔数的关联性分析[J]. 中国畜牧杂志, 2011, 47(21):1-5.
	Tian M M, Wang A G, Fu J L. Association analysis between polymorphisms of EPHA1 gene and litter size in pigs[J]. Chinese Journal of Animal Science, 2011, 47(21):1-5.
[21]	Fu Y, Knox R V, Li L, Ren S. Differential gene expression of Eph-ephrin A1 and LEPR-LEP with high or low number of embryos in pigs during implantation[J]. Reproduction in Domestic Animals, 2018, 53(4):937-946. doi: 10.1111/rda.13192. doi: 10.1111/rda.13192 pmid: 29740882
[22]	付言峰, 李兰, 周艳红, 李碧侠, 赵为民, 赵芳, 任守文. 梅山猪胚胎附植期EphA1基因表达及甲基化检测[J]. 畜牧兽医学报, 2014, 45(12):1964-1970. doi: 10.11843/j.issn.0366-6964.2014.12.008. doi: 10.11843/j.issn.0366-6964.2014.12.008
	Fu Y F, Li L, Zhou Y H, Li B X, Zhao W M, Zhao F, Ren S W. Detection of EphA1 expression and methylation during embryo implantation in Meishan pigs[J]. Acta Veterinaria et Zootechnica Sinica, 2014, 45(12):1964-1970.
[23]	周艳红, 付言峰, 赵为民, 刘红林, 任守文. 胚胎附植期梅山猪EphrinA1全编码区克隆及生物信息学分析[J]. 江苏农业学报, 2015, 31(2):350-356. doi: 10.3969/j.issn.1000-4440.2015.02.020. doi: 10.3969/j.issn.1000-4440.2015.02.020
	Zhou Y H, Fu Y F, Zhao W M, Liu H L, Ren S W. Coding sequence cloning and bioinformatics analysis of EphrinA1 during embryo implantation in Meishan swine[J]. Jiangsu Journal of Agricultural Sciences, 2015, 31(2):350-356.

胎次 Parity	窝数 Litter	总产仔数 TNB	产活仔数 NBA	健仔数 NH	断奶仔猪数 NW	出生窝质量/kg WB	出生个体均质量/kg AIWB
1	186	11.53±0.15c	10.04±0.20c	9.85±0.21c	9.60±0.22a	14.33±0.30bc	1.44±0.02a
2	186	12.35±0.15b	11.04±0.15ab	10.91±0.16ab	10.06±0.20a	15.81±0.26a	1.43±0.01a
3	186	12.92±0.18ab	11.50±0.18a	11.19±0.20a	9.68±0.24a	15.68±0.29a	1.37±0.01b
4	141	13.17±0.22ab	11.44±0.22a	11.18±0.22a	9.84±0.26a	15.49±0.33a	1.36±0.01b
5	116	12.91±0.23ab	10.91±0.23ab	10.63±0.24ab	9.50±0.27a	14.77±0.31ab	1.35±0.02b
6	77	13.30±0.30a	11.29±0.23ab	10.94±0.23ab	9.47±0.35a	15.37±0.37a	1.36±0.02b
7	51	12.67±0.40ab	10.65±0.42b	10.27±0.41bc	8.67±0.42b	13.61±0.57c	1.29±0.02c
8	17	12.94±0.70ab	-	-	-	-	-
合计Total	186	12.52±0.09ab	10.94±0.09ab	10.65±0.09ab	9.60±0.10a	14.94±0.13ab	1.37±0.01b

胎次 Parity	窝数 Litter	总产仔数 TNB	产活仔数 NBA	健仔数 NH	断奶仔猪数 NW	出生窝质量/kg WB	出生个体均质量/kg AIWB
1	186	11.53±0.15c	10.04±0.20c	9.85±0.21c	9.60±0.22a	14.33±0.30bc	1.44±0.02a
2	186	12.35±0.15b	11.04±0.15ab	10.91±0.16ab	10.06±0.20a	15.81±0.26a	1.43±0.01a
3	186	12.92±0.18ab	11.50±0.18a	11.19±0.20a	9.68±0.24a	15.68±0.29a	1.37±0.01b
4	141	13.17±0.22ab	11.44±0.22a	11.18±0.22a	9.84±0.26a	15.49±0.33a	1.36±0.01b
5	116	12.91±0.23ab	10.91±0.23ab	10.63±0.24ab	9.50±0.27a	14.77±0.31ab	1.35±0.02b
6	77	13.30±0.30a	11.29±0.23ab	10.94±0.23ab	9.47±0.35a	15.37±0.37a	1.36±0.02b
7	51	12.67±0.40ab	10.65±0.42b	10.27±0.41bc	8.67±0.42b	13.61±0.57c	1.29±0.02c
8	17	12.94±0.70ab	-	-	-	-	-
合计Total	186	12.52±0.09ab	10.94±0.09ab	10.65±0.09ab	9.60±0.10a	14.94±0.13ab	1.37±0.01b

SNP	位置/bp Position	候选基因 Candidate gene					距离/bp Distance	P值 P-value
SNP	位置/bp Position	基因ID Gene ID	所在染色体 Chr	起始物理位置/bp Start	终止物理位置/bp End	名称 Name	距离/bp Distance	P值 P-value
seq-rs329781338	8 435 604	ENSSSCG00000016488	18	8 450 516	8 460 500	SSBP1	14 912	2.94E-06
seq-rs329781338	8 435 604	ENSSSCG00000016489	18	8 463 648	8 488 381	WEE2	28 044	2.94E-06
seq-rs329781338	8 435 604	ENSSSCG00000016491	18	8 497 971	8 547 784	KIAA1147	62 367	2.94E-06

SNP	位置/bp Position	候选基因 Candidate gene					距离/bp Distance	P值 P-value
SNP	位置/bp Position	基因ID Gene ID	所在染色体 Chr	起始物理位置/bp Start	终止物理位置/bp End	名称 Name	距离/bp Distance	P值 P-value
seq-rs329781338	8 435 604	ENSSSCG00000016488	18	8 450 516	8 460 500	SSBP1	14 912	2.94E-06
seq-rs329781338	8 435 604	ENSSSCG00000016489	18	8 463 648	8 488 381	WEE2	28 044	2.94E-06
seq-rs329781338	8 435 604	ENSSSCG00000016491	18	8 497 971	8 547 784	KIAA1147	62 367	2.94E-06

SNP	位置/bp Position	候选基因Candidate gene					距离/bp Distance	P值 P-value
SNP	位置/bp Position	基因ID Gene ID	所在染色体 Chr	起始物理位置/bp Start	终止物理位置/bp End	名称 Name	距离/bp Distance	P值 P-value
seq-rs81244816	7 124 607	ENSSSCG00000016458	18	7 045 247	7 046 597	TAS2R60	78 010	2.15E-06
seq-rs81244816	7 124 607	ENSSSCG00000016459	18	7 081 288	7 098 452	EPHA1	26 155	2.15E-06
seq-rs81244816	7 124 607	ENSSSCG00000028274	18	7 134 211	7 145 293	FAM131B	9 604	2.15E-06
seq-rs81244816	7 124 607	ENSSSCG00000016462	18	7 145 222	7 181 268	CLCN1	20 615	2.15E-06
seq-rs81244816	7 124 607	ENSSSCG00000016464	18	7 187 327	7 206 997	CASP2	62 720	2.15E-06
seq-rs81244816	7 124 607	ENSSSCG00000016465	18	7 211 800	7 214 244	TMEM139	87 193	2.15E-06
seq-rs81467772	7 848 618	ENSSSCG00000016476	18	7 756 189	7 823 672	TRBV21OR9-2	24 946	3.11E-06
seq-rs81467772	7 848 618	ENSSSCG00000016478	18	7 761 371	7 764 639	PRSS2	83 979	3.11E-06
seq-rs81467772	7 848 618	ENSSSCG00000016479	18	7 776 934	7 777 450	TRBV19	71 168	3.11E-06
seq-rs81467772	7 848 618	ENSSSCG00000016477	18	7 811 001	7 811 502	TRBV24-1	37116	3.11E-06
seq-rs81467772	7 848 618	ENSSSCG00000018482	18	7 849 585	7 849 678	U6	967	3.11E-06
seq-rs329781338	8 435 604	ENSSSCG00000016488	18	8 450 516	8 460 500	SSBP1	14 912	4.69E-06
seq-rs329781338	8 435 604	ENSSSCG00000016489	18	8 463 648	8 488 381	WEE2	28 044	4.69E-06
seq-rs329781338	8 435 604	ENSSSCG00000016491	18	8 497 971	8 547 784	KIAA1147	62 367	4.69E-06
seq-rs80782154	68 369 751	ENSSSCG00000001944	7	68 332 624	68 334 022	NKX2-8	35 729	1.15E-05
seq-rs81238474	20 130 547	ENSSSCG00000025731	2	20 010 994	20 048 322	ALKBH3	82 225	1.33E-05
seq-rs81238474	20 130 547	ENSSSCG00000021739	2	20 061 566	20 211 950	HSD17B12	0	1.33E-05
seq-rs81238474	20 130 547	ENSSSCG00000028702	2	20 088 399	20 088 502	U6	42 045	1.33E-05
seq-rs340736313	57 657 259	ENSSSCG00000001804	7	57 489 036	57 582 494	HOMER2	74 765	2.01E-05
seq-rs340736313	57 657 259	ENSSSCG00000001805	7	57 596 454	57 627 963	WHAMM	29 296	2.01E-05
seq-rs340736313	57 657 259	ENSSSCG00000001806	7	57 653 600	57 679 877	FSD2	0	2.01E-05
seq-rs340736313	57 657 259	ENSSSCG00000020586	7	57 683 211	57 683 337	SCARNA15	25 952	2.01E-05

选择文件类型/文献管理软件名称

选择包含的内容