三穗鸭体尺与屠宰性状的全基因组关联分析

doi:10.15933/j.cnki.1004-3268.2025.06.015

摘要/Abstract

摘要： 为了筛选出影响三穗鸭体尺和屠宰性状的遗传位点和功能基因，以655只160日龄三穗鸭为研究对象，测定体尺（体斜长、颈长）和屠宰性状（宰前活质量、屠体质量、全净膛质量和胸腿肌率）等表型数据，通过全基因组关联分析（10×深度）获取个体单核苷酸多态性（SNP）基因型，基于PLINK 1.9进行质控分析，利用GEMMA执行混合线性模型关联分析，并通过clusterProfiler对显著关联位点进行KEGG候选基因功能富集分析。结果表明，分别筛选到与体斜长、颈长、宰前活质量、屠体质量、全净膛质量、胸腿肌率显著关联的SNP 位点207、78、64、66、68、96 个，分别注释到50、31、32、40、27、37 个候选基因。KEGG富集分析结果表明，与体斜长、颈长、宰前活质量、屠体质量、全净膛质量和胸腿肌率关联的候选基因分别富集到16、10、12、11、13、22条KEGG通路。其中，神经活性配体-受体相互作用与6个性状均有关，钙信号通路与体斜长、屠体质量、全净膛质量和胸腿肌率有关，MAPK信号通路与宰前活质量、屠体质量和全净膛质量有关。LOC119717855、GPR50、GRM7、GABRB2、GABRG2、LOC101800857等基因是体斜长的关键候选基因，GRIK2、BNIP1等基因是颈长的关键候选基因，TAF2和NRK基因是宰前活质量、屠体质量和全净膛质量的关键候选基因，SLC8A1基因是屠体质量和全净膛质量的关键候选基因，HTR2A、QRFPR、NPVF和UTS2R等基因是胸腿肌率的关键候选基因。综上，影响三穗鸭体尺与屠宰性状的主要KEGG通路为神经活性配体-受体相互作用、钙信号等信号通路，关键候选基因为
LOC119717855、GPR50、GRIK2、TAF2等。

关键词: 三穗鸭, 全基因组关联分析, 体尺性状, 屠宰性状, 候选基因

Abstract: To screen out the genetic loci and functional genes affecting the body size and slaughter traits of Sansui ducks，655 160‐day‐old Sansui ducks were taken as the research subjects，and the phenotypic data such as body size（body skew length，neck length）and slaughter traits（pre‐slaughter live weight，carcass weight，evisceration weight and chest and leg muscle rate）were measured.The individual single nucleotide polymorphism（SNP）genotypes were obtained by genome‐wide association study（10×depth），quality control analysis was conducted based on PLINK 1.9，mixed linear model association analysis was performed using GEMMA，the KEGG candidate gene function enrichment analysis of the significant association sites was performed by clusterProfiler.The results showed a total of 207，78，64，66，68 and 96 SNP loci significantly associated with body skew length，neck length，pre‐slaughter live weight，carcass weight，evisceration weight and chest and leg muscle rate，and 50，31，32，40，27 and 37 candidate genes annotated respectively.The results of KEGG enrichment analysis showed that 16，10，12，11，13 and 22 KEGG pathways were enriched for candidate genes related to body skew length，neck length，pre‐slaughter live weight，carcass weight，evisceration weight and chest and leg muscle rate.Neuroactive ligand‐receptor interaction was associated with all six traits，the calcium signaling pathway was associated with body skew length，carcass weight，eviscerated weight，and chest and leg muscle rate，and MAPK signaling pathway was associated with pre‐slaughter live weight，carcass weight，eviscerated weight.LOC119717855，GPR50，GRM7，GABRB2，GABRG2 and LOC101800857 were identified as key candidate genes for body skew length，while GRIK2 and BNIP1 were identified as key candidate genes for neck length.TAF2 and NPK genes were identified as the key candidate genes for pre‐slaughter live weight，carcass weight and evisceration weight traits. SLC8A1 gene was identified as key candidate gene for carcass weight and evisceration weight，HTR2A，QRFPR，NPVF and UTS2R genes were identified as the key candidate genes for chest and leg muscle rate.The above results showed that key signaling pathways such as neuroactive ligand‐receptor interaction and calcium signaling pathway，as well as key candidate genes such as LOC119717855，GPR50，GRIK2 and TAF2，which affected the body size and slaughter traits of Sansui ducks，were screened.

Key words: Sansui duck, Genome‐wide association study, Body size trait, Slaughter trait, Candidate gene

中图分类号:

S831.2

伍革民, 韩雪, 吴松成, 冯宇隆. 三穗鸭体尺与屠宰性状的全基因组关联分析[J]. 河南农业科学, 2025, 54(6): 135-143.

WU Gemin, HAN Xue, WU Songcheng, FENG Yulong. Genome‐wide Association Study of Body Size and Slaughter Traits in Sansui Duck[J]. Journal of Henan Agricultural Sciences, 2025, 54(6): 135-143.

图/表 6

参考文献

［1］刘宏祥，沈永杰，张丽华，等. 娄门鸭胸肌性状相关指标的全基因组关联分析［J］. 中国家禽，2024，46（4）：8‐14.
LIU H X，SHEN Y J，ZHANG L H，et al. Genome‐wide association analysis on breast muscle related traits of Loumen duck［J］. China Poultry，2024，46（4）：8‐14.

［2］刘文静，凡文磊，刘大鹏，等. 鸭胸肌肌内脂肪含量全基因组关联分析［J］. 中国畜牧杂志，57（6）：96‐101.

LIU W J，FAN W L，LIU D P，et al. Genome‐wide association analysis of intramuscular fat in duck［J］.Chinese Journal of Animal Science，2021，57（6）：96‐101.

［3］HOANG K X，MATSUZAKI M，KOHSAKA T，et al.Expression of relaxin family peptide receptors 1 and 3 in the ovarian follicle of Japanese quail［J］.The Journal of Poultry Science，2024，61：2024005.

［4］ERLANDSON S C，RAWSON S，OSEI‐OWUSU J，et al.The relaxin receptor RXFP1 signals through a mechanism of autoinhibition ［J］. Nature Chemical Biology，2023，19（8）：1013‐1021.

［5］刘畅，于向纺，滕斌，等. GPR50在肥胖及相关疾病中的作用研究进展［J］.生命科学，2021，33（1）：61‐67.

LIU C，YU X F，TENG B，et al. Research progress on the role of GPR50 in obesity and related diseases［J］.Chinese Bulletin of Life Sciences，2021，33（1）：61‐67.

［6］NOROOZI R，TAHERI M，OMRANI M D，et al.Glutamate receptor metabotropic 7（GRM7）gene polymorphisms in mood disorders and attention deficit hyperactive disorder［J］. Neurochemistry International，2019，129：104483.

［7］FREITAS G A，NISWENDER C M.GRM7 gene mutations and consequences for neurodevelopment[J］.Pharmacology，Biochemistry，and Behavior，2023，225：173546.
［8］BARKI M，XUE H. GABRB2，a key player in neuropsychiatric disorders and beyond［J］. Gene，2022，809：146021.
［9］ROMANO V，COZZOLINO D，ZINNO G，et al. The effects of beta（2）‐adrenergic receptors activation on the cardiovascular system and on the skeletal muscle：Anarrative review［J］. Journal of Human Sport and Exercise，2021，3：22‐23.
［10］GAMBARDELLA J，FIORDELISI A，AVVISATO R，et al. Adrenergic receptors in endothelial and vascular smooth muscle cells［J］. Current Opinion in Phsiology，2023，36：100721.
［11］KAYE W H，FRANK G K，BAILER U F，et al.Neurobiology of anorexia nervosa：Clinical implications of alterations of the function of serotonin and other neuronal systems［J］. International Journal of Eating Disorders，2005，37（S1）：S15‐S19.
［12］王思月，云巾宴，刘宇，等. 中国草原红牛HTR2A基因多态性及与肉质性状的相关性分析［J］. 东北农业科学，2023，48（6）：96‐100.
WANG S Y，YUN J Y，LIU Y，et al. Analysis of the polymorphism of HTR2A gene in Chinese Caoyuan red cattle and its correlation with meat quality traits［J］.Journal of Northeast Agricultural Sciences，2023 ，48（6）：96‐100.
［13］BAKIU R，TOLOMEO A M，SANTOVITO G. Positive selection effects on the biochemical properties of fish pyroglutamylated RFamide peptide receptor（QRFPR）［J］.Italian Journal of Zoology，2015，82（4）：460‐472.
［14］KOLLER J，HERZOG H，ZHANG L. The distribution of neuropeptide FF and neuropeptide VF in central and peripheral tissues and their role in energy homeostasis control［J］.Neuropeptides，2021，90：102198.
［15］UBUKA T，MORGAN K，PAWSON A J，et al.Identification of human GnIH homologs，RFRP‐1 and RFRP‐3，and the cognate receptor，GPR147 in the human hypothalamic pituitary axis［J］. PLoS One，2009，4（12）：e8400.
［16］YU H P，WANG Y Y，GAO J J，et al.Application of the neuropeptide NPVF to enhance angiogenesis and osteogenesis in bone regeneration［J］.Communications Biology，2023，6（1）：197.
［17］CONSTANTIN S，PIZANO K，MATSON K，et al. An inhibitory circuit from brainstem to GnRH neurons in male mice：A new role for the RFRP receptor［J］.Endocrinology，2021，162（5）：bqab030.

［18］CUI L，LV C，ZHANG J N，et al. Characterization of four urotensin Ⅱ receptors（UTS2Rs）in chickens［J］.Peptides，2021，138：170482.

［19］DAI Z C，WANG Y W，WU Z C，et al. Novel mutations in UTS2R are associated with adolescent idiopathic scoliosis in the Chinese population［J］. Spine，2021，46 （5）：E288‐E293.

［20］SASAZAKI S，AKIYAMA K，NARUKAMI T，et al. UTS2R gene polymorphisms are associated with fatty acid composition in Japanese beef cattle［J］. Animal Science Journal，2014，85（5）：499‐505.

［21］YU L J，WALL B A，WANGARI‐TALBOT J，et al.Metabotropic glutamate receptors in cancer［J］. mGLU Receptors，2017，31：161‐182.

［22］SZCZESNIAK A M，GILBERT R W，MUKHIDA M，et al. Mechanical loading modulates glutamate receptor subunit expression in bone［J］. Bone，2005，37（1）：63‐73.

［23］HOLLING T，BHAVANI G S，VON ELSNER L，et al.A homozygous hypomorphic BNIP1 variant causes an increase in autophagosomes and reduced autophagic flux and results in a spondylo‐epiphyseal dysplasia［J］.Human Mutation，2022，43（5）：625‐642.

［24］CHENG I H，PI W C，HSU C H，et al. TAF2 within the TFIID complex，regulates the expression of a subset of protein‐coding genes［J］. Cell Death Discovery，2024，10（1）：244.

［25］DENDA K，NAKAO‐WAKABAYASHI K，OKAMOTO N，et al. Nrk，an X‐linked protein kinase in the germinal center kinase family，is required for placental development and fetoplacental induction of labor［J］.The Journal of Biological Chemistry，2011，286（33）：28802‐28810.

［26］洪龙胜，许丹宁，李婉雁，等. 马岗鹅早期发育阶段腿肌lncRNA-miRNA-mRNA网络构建与分析［J］.中国家禽，2024，46（9）：34‐45.
HONG L S，XU D N，LI W Y et al. Construction and analysis on lncRNA‐miRNA‐mRNA network in the early developmental stage of leg muscle in Magang goose［J］.China Poultry，2024，46（9）：34‐45.

［27］DVORAK V，WIEDMER T，INGLES‐PRIETO A，et al.An overview of cell‐based assay platforms for the solute carrier family of transporters［J］.Frontiers in Pharmacology，2021，12：722889.

[1]	黄世鹏, 何虎强, 李光照, 陈传鑫, 钟拓, 谢剑波, 周湘晚, 王丹, 肖应辉. 水稻苗期耐旱性全基因组关联分析及候选基因预测[J]. 河南农业科学, 2025, 54(7): 12-20.
[2]	车雅丽1, 吴苏峻1, 王坤杨2, 王艳丽1, 张锴1, 乔亚科1, 袁泉3. 镉胁迫下野生大豆苗期性状全基因组关联分析及耐镉候选基因筛选[J]. 河南农业科学, 2025, 54(1): 40-54.
[3]	冯培媛, 董悦, 焦紫岚, 陈敏, 孔维儒, 冉杰, 李培富, 田蕾. 粳稻种质苗期耐盐性综合评价与全基因组关联分析[J]. 河南农业科学, 2024, 53(7): 1-9.
[4]	李莹, 何静怡, 陈鹏, 张丽, 许欣纯, 计坚, 罗成龙 . 麻黄鸡体质量与体尺性状的相关分析和主成分分析[J]. 河南农业科学, 2022, 51(5): 126-132.
[5]	李培艳, 谭元生, 韦金兑, 罗杰灵, 华国洪, 李华, 甘建伉. 清远麻鸡母系体尺性状与产蛋数的相关分析[J]. 河南农业科学, 2021, 50(9): 143-148.
[6]	邢宝松, 王璟, 白献晓, 陈俊峰, 张家庆, 任巧玲, 郭红霞, 张华, 曹海. 不同养殖模式对淮南猪屠宰和肉质性状的影响[J]. 河南农业科学, 2021, 50(8): 160-165.
[7]	董静超, 甘建伉, 李华, 张正芬, 康慧敏. 清远麻鸡体尺与屠宰性状的因子分析及通径分析[J]. 河南农业科学, 2021, 50(4): 162-168.
[8]	陈建兴, 童家兴, 张孝忠, 张向阳, 王宇鑫, 孙玉江. 山东小毛驴全基因组选择信号检测[J]. 河南农业科学, 2021, 50(2): 145-150.
[9]	黄得纯, 吕敏芝, 池仕红, 康慧敏, 蔡影峰, 甘仕堂. 东涛鸡体质量与体尺性状的相关分析及通径分析[J]. 河南农业科学, 2020, 49(8): 130-135.
[10]	吴金华, 李庆山, 谭淑雯, 李华, 景栋林. 不同日龄樱桃谷鸭体质量与体尺性状的主成分分析[J]. 河南农业科学, 2020, 49(7): 139-145.
[11]	谭光辉, 张依裕, 吴磊, 覃媛钰, 李杰章. 三穗鸭ORAI1基因SNPs筛查及生物信息学分析[J]. 河南农业科学, 2020, 49(4): 153-159.
[12]	贾小平,张博，董志平，全建章，王永芳，张小梅，袁玺垒，李剑峰，戴凌峰. 海南短日照条件下谷子穗部性状的全基因组关联分析[J]. 河南农业科学, 2018, 47(9): 33-40.
[13]	姜成红,耿鑫鑫,魏文辉,姜慧芳. 甘蓝型油菜株高QTL定位及主效QTL区间候选基因预测[J]. 河南农业科学, 2017, 46(8): 27-31.
[14]	赵莉媛，赵香珍，张娟. 海东鸡CAPN1基因3′UTR多态性及其与体尺性状的关联性[J]. 河南农业科学, 2016, 45(2): 138-141.
[15]	叶伟庆，王光琴，杨芬霞，刘晓佳，林俊全，杜炳旺. 信宜怀乡鸡体质量与体尺性状的相关性及聚类分析[J]. 河南农业科学, 2015, 44(2): 132-134.