杜洛克猪×二花脸猪F2代杂交猪免疫性状全基因组关联分析

doi:10.16431/j.cnki.1671-7236.2025.04.001

摘要/Abstract

摘要： 【目的】猪是重要的农业动物，其免疫性状的解析对于提高生产效率和抗病能力至关重要。本研究旨在鉴别与杜洛克猪×二花脸猪F2代杂交猪免疫性状相关的候选基因，为猪抗病育种提供重要的候选基因信息。【方法】采集294头杜洛克猪×二花脸猪F2代杂交猪出生后24 h内的耳、尾组织样本，提取基因组DNA，并在35日龄时采集血液样本用于细胞因子检测。通过对杜洛克猪×二花脸猪F2代杂交猪资源群体的细胞因子(IFN-α、IFN-γ、IGG、IL8、IL10)相关免疫性状进行全基因组关联分析(GWAS)，识别与性状相关的单核苷酸多态性(SNP)位点，使用IFN-γ/IL10值(IFNGIL10)进一步鉴定影响猪免疫性状的相关SNP位点与候选基因，对鉴定的相关基因进行GO功能与KEGG通路富集分析，并对候选基因编码蛋白进行蛋白互作(PPI)分析。【结果】GWAS共鉴定出10个SNPs，这些SNPs定位到21个与免疫性状相关的候选基因：CH25H、LIPA、SIGMAR1、PIP5KIC、PLCZ1、PIK3C2G、CNTFR、IL11RA、SAP30、HMGB2、CCL21、CCL19、CCL27、GNA11、GNA15、DSCC1、TERF2IP、FZR1、ENPP2、CACTIN和FAS。GO功能富集结果显示，候选基因显著富集在趋化因子活性、细胞对白细胞介素-1的反应和炎症反应等条目。KEGG通路富集结果显示，候选基因显著富集在细胞因子与受体结合通路、肌醇磷酸代谢、钙信号通路等通路。PPI分析共构建了12个子网络，表明这些候选基因之间存在潜在联系，在同一通路上进行显著富集，并共同参与生物过程以发挥其生物功能。【结论】本研究共发现10个与猪免疫性状关联的SNPs位点，结合富集分析与PPI网络筛选到21个目标性状候选基因。本研究结果不仅加深了对猪免疫特性分子基础的理解，也为猪的抗病育种提供了理论依据。

关键词: 猪; 细胞因子; 免疫性状; 全基因组关联分析; 候选基因

Abstract: 【Objective】 Pigs play a crucial role in agriculture,and the analysis of their immune traits holds significant importance for enhancing production efficiency and disease resistance.The objective of this study was to identify potential candidate genes associated with immune traits in Duroc×Erhualian F2 generation hybrid pigs,thereby providing valuable gene information for breeding programs aimed at improving porcine disease resistance.【Method】 A total of 294 Duroc×Erhualian F2 generation hybrid pigs were collected from ear and tail tissue samples within 24 hours after birth,and genomic DNA was extracted.Blood samples were collected at 35 days of age to detect cytokines.Genome-wide association study (GWAS) was performed on the immune traits related to cytokines (IFN-α,IFN-γ,IGG,IL8 and IL10) in Duroc×Erhualian F2 generation hybrid pig resource population to identify single nucleotide polymorphism (SNP) site related to the traits.On this basis,the IFN-γ/IL10 value (IFNGIL10) was utilized for further identification of the related SNP and candidate gene influencing the immune traits of pigs.Subsequently,GO function and KEGG pathway enrichment analysis were performed on the identified related genes,and protein-protein interaction (PPI) analysis was performed on the proteins encoded by the candidate genes.【Result】 A total of 10 SNPs were identified through GWAS analysis,which specifically targeted 21 candidate genes associated with immune traits (CH25H,LIPA,SIGMAR1,PIP5KIC,PLCZ1,PIK3C2G,CNTFR,IL11RA,SAP30,HMGB2,CCL21,CCL19,CCL27,GNA11,GNA15,DSCC1,TERF2IP,FZR1,ENPP2,CACTIN and FAS).GO function enrichment analysis results indicated that these candidate genes were significantly enriched in terms such as chemokine activity,cellular response to interleukin-1 and inflammatory response.KEGG enrichment analysis results demonstrated that the candidate genes were significantly enriched in pathways such as cytokine-cytokine receptor interaction,inositol phosphate metabolism and calcium signaling pathway.Furthermore,the PPI analysis constructed a total of 12 subnetworks,suggesting potential interactions among these candidate genes,which were notably enriched within the same pathway,collaboratively participating in biological processes to fulfill their respective functions.【Conclusion】 The present study identified 10 SNPs associated with porcine immune traits through GWAS,and further selected 21 candidate genes for target traits by integrating enrichment analysis and PPI analysis.Overall,the findings from this study not only deepened our understanding of the molecular basis underlying porcine immune traits,but also provided a theoretical foundation for porcine disease resistance breeding.

Key words: pigs; cytokines; immune traits; genome-wide association study; candidate genes

中图分类号:

S828.2
Q78

苗娜, 乔嘉坤, 杨慧, 韩萍萍, 许方军, 车兆轩, 代翔毓, 徐明航, 龙志伟, 朱猛进. 杜洛克猪×二花脸猪F2代杂交猪免疫性状全基因组关联分析[J]. 中国畜牧兽医, 2025, 52(4): 1455-1467.

MIAO Na, QIAO Jiakun, YANG Hui, HAN Pingping, XU Fangjun, CHE Zhaoxuan, DAI Xiangyu, XU Minghang, LONG Zhiwei, ZHU Mengjin. Genome Wide Association Study of Immune Traits in Duroc×Erhualian F2 Generation Pigs[J]. China Animal Husbandry and Veterinary Medicine, 2025, 52(4): 1455-1467.

参考文献

[1] 王贵荣.2023年农业经济形势总体良好[Z].中国经济网.2024.http://www.ce.cn/xwzx/gnsz/gdxw/202401/202418/t20240118_38870839.shtml.WANG G R.The overall agricultural economic situation in 2023 is positive[Z].China Economic Network.2024.http://www.ce.cn/xwzx/gnsz/gdxw/202401/202418/t20240118_38870839.shtml.(in Chinese)
[2] 樊品科.猪抗病育种的研究进展[J].中国猪业,2021,8(5):36-38.FAN P K.Research progress on disease resistant breeding in pigs[J].China Swine Industry,2021,8(5):36-38.(in Chinese)
[3] 王超,赵书红,朱猛进.猪抗病育种的相关问题及研究进展[J].中国畜牧杂志,2014,50(22):67-72.WANG C,ZHAO S H,ZHU M J.Related problems and research progress in anti-disease breeding of pig[J].Chinese Journal of Animal Science,2014,50(22):67-72.(in Chinese)
[4] DENG S X,QIU Y B,ZHUANG Z W,et al.Genome-wide association study of body conformation traits in a three-way crossbred commercial pig population[J].Animals,2023,13(15):2414.
[5] PAIRO-CASTINEIRA E,RAWLIK K,BRETHERICK A D,et al.GWAS and meta-analysis identifies 49 genetic variants underlying critical COVID-19(vol 617,pg 764,2023)[J].Nature,2023,619(7971):E61.
[6] SAHITO J H,ZHANG H,GISHKORI Z G N,et al.Advancements and prospects of genome-wide association studies (GWAS) in maize[J].International Journal of Molecular Sciences,2024,25(3):1918.
[7] GUO Y M,HUANG Y X,HOU L J,et al.Genome-wide detection of genetic markers associated with growth and fatness in four pig populations using four approaches[J].Genetics Selection Evolution,2017,49:1-11.
[8] ZHANG Z,CHEN Z T,YE S P,et al.Genome-wide association study for reproductive traits in a Duroc pig population[J].Animals,2019,9(10):732.
[9] LAN Q,DENG Q C,QI S J,et al.Genome-wide association analysis identified variants associated with body measurement and reproduction traits in Shaziling pigs[J].Genes,2023,14(2):522.
[10] ZHOU S P,DING R R,ZHUANG Z W,et al.Genome-wide association analysis reveals genetic loci and candidate genes for chest,abdominal,and waist circumferences in two Duroc pig populations[J].Frontiers in Veterinary Science,2022,8:807003.
[11] 杨生海,殷宏,刘永生,等.干扰素-γ研究进展[J].生物技术通报,2010,8:29.YANG S H,YIN H,LIU Y S,et al.Advances in research of interferon-gamma[J].Biotechnology Bulletin,2010,8:29.(in Chinese)
[12] FUJIO K,OKAMURA T,YAMAMOTO K.The family of IL-10-secreting CD4⁺ T cells[J].Advances in Immunology,2010,105:99-130.
[13] PURCELL S,NEALE B,TODD-BROWN K,et al.PLINK:A tool set for whole-genome association and population-based linkage analyses[J].American Journal of Human Genetics,2007,81(3):559-575.
[14] BROWNING B L,TIAN X W,ZHOU Y,et al.Fast two-stage phasing of large-scale sequence data[J].American Journal of Human Genetics,2021,108(10):1880-1890.
[15] LIPKA A E,TIAN F,WANG Q S,et al.GAPIT:Genome association and prediction integrated tool[J].Bioinformatics,2012,28(18):2397-2399.
[16] HUANG M,LIU X,ZHOU Y,et al.BLINK:A package for the next level of genome-wide association studies with both individuals and markers in the millions[J].Gigascience,2019,8(2):giy154.
[17] DURINCK S,SPELLMAN P T,BIRNEY E,et al.Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt[J].Nature Protocols,2009,4(8):1184-1191.
[18] WU J M,MAO X Z,CAI T,et al.KOBAS server:A web-based platform for automated annotation and pathway identification[J].Nucleic Acids Research,2006,34:W720-W724.
[19] SNEL B,LEHMANN G,BORK P,et al.STRING:A web-server to retrieve and display the repeatedly occurring neighbourhood of a gene[J].Nucleic Acids Research,2000,28(18):3442-3444.
[20] 王爽,肖俊,靳苗,等.胆固醇代谢与免疫反应[J].生物医学转化,2024,5(2):36-45.WANG S,XIAO J,JIN M,et al.Advances in research of interferon-gamma[J].Biomedical Transformation,2024,5(2):36-45.(in Chinese)
[21] 林家俞,秦洁洁,蒋玲曦.肿瘤微环境中免疫细胞的代谢研究进展[J].上海交通大学学报(医学版),2022,42(8):1122-1130.LIN J Y,QIN J J,JIANG L X.Progress in metabolism of the immune cells in tumor microenvironment[J].Journal of Shanghai Jiao Tong University (Medical Science),2022,42(8):1122-1130.(in Chinese)
[22] MILLER A T,CHAMBERLAIN P P,COOKE M P.Beyond IP3-roles for higher order inositol phosphates in immune cell signaling[J].Cell Cycle,2008,7(4):463-467.
[23] ZENG H,YANG K,CLOER C,et al.mTORC1 couples immune signals and metabolic programming to establish T-cell function[J].Nature,2013,499(7459):485-490.
[24] JIA Y H,SCHURMANS S,LUO H B R.Regulation of innate immunity by inositol 1,3,4,5-tetrakisphosphate[J].Cell Cycle,2008,7(18):2803-2808.
[25] HOLMES R S,VANDEBERG J L,COX L A.Genomics and proteomics of vertebrate cholesterol ester lipase (LIPA) and cholesterol 25-hydroxylase (CH25H)[J].3 Biotech,2011,1(2):99-109.
[26] LUND E G,KERR T A,SAKAI J,et al.cDNA cloning of mouse and human cholesterol 25-hydroxylases,polytopic membrane proteins that synthesize a potent oxysterol regulatory of lipid metabolism[J].Journal of Biological Chemistry,1998,273(51):34316-34327.
[27] BAUMAN D R,BITMANSOUR A D,MCDONALD J G,et al.25-hydroxycholesterol secreted by macrophages in response to Toll-like receptor activation suppresses immunoglobulin A production[J].Proceedings of the National Academy of Sciences of the United States of America,2009,106(39):16764-16769.
[28] AISHWARYA R,ABDULLAH C S,MORSHED M,et al.Sigmar1's molecular,cellular,and biological functions in regulating cellular pathophysiology[J].Frontiers in Physiology,2021,12:705575.
[29] YANG J Z,ZHANG K K,SHEN H W,et al.Sigma-1 receptor knockout disturbs gut microbiota,remodels serum metabolome,and exacerbates isoprenaline-induced heart failure[J].Frontiers in Microbiology,2023,14:1255971.
[30] 吴葩,许彩民.Ⅲ型磷脂酰肌醇3-激酶复合物与自噬的研究进展[J].生物化学与生物物理进展,2012,39(3):210-216.WU P,XU C M.Advances in relationship between class Ⅲ PI3K complex and autophagy[J].Progress in Biochemistry and Biophysics,2012,39(3):210-216.(in Chinese)
[31] CHANDRASEKARAN S,SASAKI M,POLLACK B P.Linking growth factors and immune actors:PTEN/PI3K signaling in cancer and the response to cytokines[J].Cancer Research,2019,79(13_Supplement):4510.
[32] CHU C A,WANG Y W,CHEN Y L,et al.The role of phosphatidylinositol 3-kinase catalytic subunit type 3 in the pathogenesis of human cancer[J].International Journal of Molecular Sciences,2021,22(20):10964.
[33] NGUYEN T,DEENICK E K,TANGYE S G.Phosphatidylinositol 3-kinase signaling and immune regulation:Insights into disease pathogenesis and clinical implications[J].Expert Review of Clinical Immunology,2021,17(8):905-914.
[34] SHAN K S,BONANO-RIOS A,THEIK N W Y,et al.Molecular targeting of the phosphoinositide-3-protein kinase (PI3K) pathway across various cancers[J].International Journal of Molecular Sciences,2024,25(4):1973.
[35] DWIJAYANTI D R,WIDYANANDA M H,HERMANTO F E,et al.Revealing the anti-inflammatory activity of extract:Transcriptomic and nitric oxide production analysis in LPS-induced RAW264.7 cells[J].Food and Agricultural Immunology,2024,35(1):2351360.
[36] GARBERS C,LOKAU J.Cytokines of the interleukin-6 family as emerging targets in inflammatory bowel disease[J].Expert Opinion on Therapeutic Targets,2024,28(1-2):57-65.
[37] MEAZZA C,DI MARCO A,FRUSCELLA P,et al.Centrally mediated inhibition of local inflammation by ciliary neurotrophic factor[J].Neuroimmunomodulation,1997,4(5-6):271-276.
[38] REHMAN A U,WANG Z H,QIN Q S,et al.Enhancing antitumor immunity and achieving tumor eradication with IL11RA mRNA immunotherapy[J].International Immunopharmacology,2024,134:112205.
[39] BOONYARATANAKORNKIT V,MELVIN V,PRENDERGAST P,et al.High-mobility group chromatin proteins 1 and 2 functionally interact with steroid hormone receptors to enhance their DNA binding in vitro and transcriptional activity in mammalian cells[J].Molecular and Cellular Biology,1998,18(8):4471-4487.
[40] REDMOND A M,BYRNE C,BANE F T,et al.Genomic interaction between ER and HMGB2 identifies DDX18 as a novel driver of endocrine resistance in breast cancer cells[J].Oncogene,2015,34(29):3871-3880.
[41] SHI H Q,XU H,CHAI C P,et al.Integrated bioinformatics analysis of potential biomarkers for pancreatic cancer[J].Journal of Clinical Laboratory Analysis,2022,36(5):e24381.
[42] STARKOVA T,POLYANICHKO A,TOMILIN A N,et al.Structure and functions of HMGB2 protein[J].International Journal of Molecular Sciences,2023,24(9):8334.
[43] SATO M,KADOMATSU T,MIYATA K,et al.The lncRNA antagonizes heart failure by inactivating DNA damage response and activating mitochondrial biogenesis[J].Nature Communications,2021,12(1):2529.
[44] COMERFORD I,HARATA-LEE Y,BUNTING M D,et al.A myriad of functions and complex regulation of the CCR7/CCL19/CCL21 chemokine axis in the adaptive immune system[J].Cytokine&Growth Factor Reviews,2013,24(3):269-283.
[45] HOMEY B,ALENIUS H,MVLLER A,et al.CCL27-CCR10 interactions regulate T cell-mediated skin inflammation[J].Nature Medicine,2002,8(2):157-165.
[46] CHEN H P,YAO W S,JIN D,et al.Cloning,expression pattern,chromosomal localization,and evolution analysis of porcine GNAQ,GNA11 and GNA14[J].Biochemical Genetics,2008,46(7-8):398-405.
[47] SODHI A,MONTANER S,GUTKIND J S.Viral hijacking of G-protein-coupled-receptor signaling networks[J].Nature Reviews Molecular Cell Biology,2004,5(12):998-1012.
[48] WANG N,HUANG H J,XIONG Q Q,et al.GNAQ negatively regulates antiviral innate immune responses in a calcineurin-dependent manner[J].Journal of Immunology,2019,203(5):1288-1297.
[49] FAHRER J.Decoding the genetic program of micronucleus formation:Linking chromosomal instability to human disease[J].Signal Transduction and Targeted Therapy,2024,9(1):162.
[50] ADAMS D J,BARLAS B,MCINTYRE R E,et al.Genetic determinants of micronucleus formation in vivo [J].Nature,2024,627(8002):130-136.
[51] DE LANGE T.Shelterin:The protein complex that shapes and safeguards human telomeres[J].Genes&Development,2005,19(18):2100-2110.
[52] 亓雨晴,刘晓然.运动调控端粒长度机制与健康促进[J].中国组织工程研究,2024,28(23):3759.QI Y Q,LIU X R. Mechanisms of exercise-regulated telomere length and health promotion[J].Chinese Journal of Tissue Engineering Research,2024,28(23):3759.(in Chinese)
[53] MIRABELLO L,HUANG W Y,WONG J Y,et al.The association between leukocyte telomere length and cigarette smoking,dietary and physical variables,and risk of prostate cancer[J].Aging Cell,2009,8(4):405-413.
[54] MIRABELLO L,YEAGER M,CHOWDHURY S,et al.Worldwide genetic structure in 37 genes important in telomere biology[J].Heredity,2012,108(2):124-133.
[55] WONG J M Y,COLLINS K.Telomere maintenance and disease[J].Lancet,2003,362(9388):983-988.
[56] ZSURKA G,PEEVA V,KOTLYAR A,et al.Is there still any role for oxidative stress in mitochondrial DNA-dependent aging?[J].Genes,2018,9(4):175.
[57] GARCÍA-HIGUERA I,MANCHADO E,DUBUS P,et al.Genomic stability and tumour suppression by the APC/C cofactor Cdh1[J].Nature Cell Biology,2008,10(7):802-811.
[58] LEE J H,CHOI S Y,PARK S Y,et al.Enpp2 expression by dendritic cells is a key regulator in migration[J].Biomedicines,2021,9(11):1727.
[59] LU C,HU D,ZHENG J E,et al.A six-gene risk model based on the immune score reveals prognosis in intermediate-risk acute myeloid leukemia[J].Biomed Research International,2022,2022:4010786.
[60] GARGAN S.Probing the mechanisms underlying the regulatory effects of cactin in innate immune signalling pathways[D].Maynooth:National University of Ireland,2011.
[61] CIESIELSKA-FIGLON K,LISOWSKA K A.The role of the CD28 family receptors in T-cell immunomodulation[J].International Journal of Molecular Sciences,2024,25(2):1274.
[62] STRASSER A,JOST P J,NAGATA S.The many roles of FAS receptor signaling in the immune system[J].Immunity,2009,30(2):180-192.