晋南牛遗传结构特征及选择信号分析

doi:10.16431/j.cnki.1671-7236.2024.01.017

摘要/Abstract

摘要： 【目的】阐明晋南牛的遗传结构特征，通过选择信号检测挖掘与晋南牛经济性状相关的候选基因，探究其在进化过程中的受选择情况。【方法】对晋南牛和红安格斯牛全基因组测序数据进行分析，鉴定2个群体的单核苷酸多态性（single nucleotide polymorphism，SNP）标记，分析其在基因组的位置及其结构特征，基于SNP信息进行主成分分析（PCA）、构建状态同源矩阵（IBS）；采用群体遗传分化指数（Fst）和核苷酸多样性比值（θ_π）方法联合筛选晋南牛基因组受到强烈选择的区域，并对筛选到的受选择基因进行数量性状基因座（QTL）定位、GO功能和KEGG通路富集分析。【结果】晋南牛群体SNPs位点主要分布于基因间区域，其次位于内含子区域。PCA和IBS分析结果表明，晋南牛和红安格斯牛2个群体间不存在杂交现象，且晋南牛群体中个体间遗传距离较远。通过Fst和θ_π联合分析共筛选到188个潜在受选择区域。QTL分析结果表明晋南牛的选择信号多与生长、肉质及抗病性状相关。GO功能和KEGG通路富集分析显示，筛选到晋南牛强受选择的与经济性状相关的候选基因11个，包括组织蛋白酶1（CATHL1）、CATHL3、组织蛋白酶抗菌肽（CAMP）、CATHL4、RAS1激酶抑制因子（KSR1）、梅氏同源框1（MEIS1）、不规则片段极性蛋白1（DVL1）、多梳组无名指1（PCGF1）、RB转录辅助抑制因子1（RB1）、再生家族成员4（REG4）和卷曲类受体7（FZD7），在其外显子区域均检测到了SNP位点；除MEIS1、PCGF1、PZD7基因外，其他基因均存在非同义突变位点，且KSR1基因存在获得终止密码子的突变，CAMP基因存在获得终止密码子的突变和缺失终止密码子的突变；与红安格斯牛相比，晋南牛存在6个特有的受选择基因：CATHL1、CATHL4、CAMP、MEIS1、PCGF1、PZD7，主要与抗病性、生长和骨骼肌发育等相关。【结论】本研究从全基因组水平探究了晋南牛的遗传结构及选择信号特点，筛选出6个可能与晋南牛经济性状相关的候选基因，为后续晋南牛保种选育、特色性状形成的分子机制研究提供理论参考。

关键词: 全基因组重测序; 晋南牛; 红安格斯牛; 遗传距离; 选择信号

Abstract: 【Objective】 The purpose of this study was to elucidate the genetic structure characteristics of Jinnan cattle, and explore candidate genes related to economic traits in Jinnan cattle through selection signal detection, in order to explore the selection status during the evolutionary process.【Method】 Whole genome sequencing data of Jinnan cattle and Red Angus cattle were analyzed to identify the single nucleotide polymorphism (SNP) marker of the two populations, and analyze its position and structural characteristics in the genome.Principal component analysis (PCA) and identity by state (IBS) matrix were constructed based on SNP information.A combination of population genetic differentiation index (Fst) and nucleotide diversity ratio (θ_π) was used to screen the highly selected regions in Jinan cattle genome, and quantitative trait locus (QTL), GO function and KEGG pathway enrichment were analyzed for the selected genes.【Result】 The SNPs in Jinnan cattle population were mainly distributed in the intergene region, followed by the intron region.The results of PCA and IBS analysis showed that there was no hybridization between Jinnan cattle and Red Angus cattle, and the genetic distance between individuals in Jinnan cattle was far away.By the combined analysis of Fst and θ_π, a total of 188 potential selected regions were screened.QTL analysis showed that the selection signals of Jinnan cattle were mainly related to growth, meat quality and disease resistance traits.GO function and KEGG pathway enrichment analysis showed that 11 candidate genes related to economic traits were selected strongly in Jinnan cattle, including cathepsin 1 (CATHL1), CATHL3, cathepsin antimicrobial peptide (CAMP), CATHL4, kinase inhibitor of RAS1 (KSR1), Meis homeobox 1 (MEIS1), irregular fragment polar protein 1 (DVL1), polycomb group ring finger 1 (PCGF1), RB transcriptional repressor 1 (RB1), regenerating family member 4 (REG4) and frizzled class receptor 7 (FZD7), and SNP sites were detected in the exon region of the genes.In addition to MEIS1, PCGF1 and FZD7 genes, all other genes had non-synonymous mutation sites, and KSR1 gene had stop gain mutation, and CAMP gene had stop gain and stop loss mutation.Compared with Red Angus cattle, there were 6 unique selected genes (CATHL1, CATHL4, CAMP, MEIS1, PCGF1 and PZD7 genes) in Jinnan cattle, which were mainly related to disease resistance, growth and skeletal muscle development.【Conclusion】 This study explored the genetic structure and selection signal characteristics of Jinnan cattle from the whole genome level, and screened 6 candidate genes that might be related to the economic traits of Jinnan cattle, providing theoretical reference for the subsequent research on conservation and breeding and the molecular mechanism of characteristic trait formation in Jinnan cattle.

Key words: whole genome resequencing; Jinnan cattle; Red Angus cattle; gnetic distance; selection signal

中图分类号:

S823.8

戎艳花, 贾雪纯, 李鹏飞, 田国富, 朱芷葳. 晋南牛遗传结构特征及选择信号分析[J]. 中国畜牧兽医, 2024, 51(1): 160-171.

RONG Yanhua, JIA Xuechun, LI Pengfei, TIAN Guofu, ZHU Zhiwei. Analysis of Genetic Structure Characteristics and Selection Signal in Jinnan Cattle[J]. China Animal Husbandry and Veterinary Medicine, 2024, 51(1): 160-171.

参考文献

[1] 郝王宝.浅谈晋南牛的保种与开发技术[J].农业与技术, 2018, 38(12):134. HAO W B.A brief discussion on the conservation and development techniques of Jinnan cattle[J].Agriculture and Technology, 2018, 38(12):134.(in Chinese)
[2] 呙明鹏, 孟源, 王宏浩, 等.晋南牛不同生长阶段体重和体尺性状遗传参数估计[J].畜牧兽医学报, 2023, 54(4):1452-1464. GUO M P, MENG Y, WANG H H, et al.Estimation of genetic parameters for body weight and size traits in Jinnan cattle at different growth stages[J].Acta Veterinaria et Zootechnica Sinica, 2023, 54(4):1452-1464.(in Chinese)
[3] 王宏浩, 任小康, 张毅, 等.基因芯片技术在晋南牛种公牛选育中的应用[J].畜牧兽医学报, 2021, 52(10):2803-2813. WANG H H, REN X K, ZHANG Y, et al.The application of gene chip technology in the breeding of Jinnan bull breeding[J].Acta Veterinaria et Zootechnica Sinica, 2021, 52(10):2803-2813.(in Chinese)
[4] 王志鹏, 景军红, 周萌.鸡全基因组选择信号筛选研究进展[J].东北农业大学学报, 2016, 47(10):100-106. WANG Z P, JING J H, ZHOU M.Research progress in chicken genome wide selection signal screening[J].Journal of Northeast Agricultural University, 2016, 47(10):100-106.(in Chinese)
[5] FLORI L, FRITZ S, JAFFREZIC F, et al.The genome response to artificial selection:A case study in dairy cattle[J].PLoS One, 2009, 4(8):e6595.
[6] SABETI P C, VARILLY P, FRY B, et al.Genome-wide detection and characterization of positive selection in human populations[J].Nature, 2007, 449(7164):913-918.
[7] MA Y L, ZHANG S X, ZHANG K L, et al.Genomic analysis to identify signatures of artificial selection and loci associated with important economic traits in Duroc pigs[J].G3(Bethesda), 2018, 8(11):3617-3625.
[8] 韩威, 朱云芬, 殷建玫, 等.基于RAD-Seq简化基因组测序的19个地方鸡种遗传进化研究[J].畜牧兽医学报, 2020, 51(4):670-678. HAN W, ZHU Y F, YIN J M, et al.Genetic evolution of 19 local chicken breeds based on RAD-Seq simplified genome sequencing[J].Acta Veterinaria et Zootechnica Sinica, 2020, 51(4):670-678.(in Chinese)
[9] 金美林, 陆健, 费晓娟, 等.全基因组选择信号揭示绒山羊绒毛性状相关的候选基因[J].畜牧兽医学报, 2020, 51(12):2991-3000. JIN M L, LU J, FEI X J, et al.Genome wide selection signals reveal candidate genes related to cashmere traits in cashmere[J].Acta Veterinaria et Zootechnica Sinica, 2020, 51(12):2991-3000.(in Chinese)
[10] LIU L L, FANG C, LIU W J.Identification on novel locus of dairy traits of Kazakn horse in Xinjiang[J].Gene, 2018, 677:105-110.
[11] DANIEL G, HENNER S, SABER Q.A two-step method for detecting selection signatures using genetic markers[J].Genetics Research, 2010, 92(2):141-155.
[12] DOYLE J L, BERRY D P, VEERKAMP R F, et al.Genomic regions associated with muscularity in beef cattle differ in five contrasting cattle breeds[J].Genetics, Selection, Evolution, 2020, 52(1):2.
[13] TAYE M, YOON J, DESSIE T.Deciphering signature of selection affecting beef quality traits in Angus cattle[J].Genes Genomics, 2018, 40(1):63-75.
[14] ALBERTI P, PANEA B, SANUDO C, et al.Live weight, body size and carcass characteristics of young bulls of fifteen European breeds[J].Livestock Science, 2008, 114(1):19-30.
[15] 胡宏伟, 杜宝强, 张成虎.秦川牛与不同品种肉牛杂交F1代生长性能比较研究[J].甘肃畜牧兽医, 2022, 52(8):20-23. HU H W, DU B Q, ZHANG C H.Cross F1 between Qinchuan cattle and different breeds of beef cattle comparative study on the growth performance of the first generation[J].Gansu Livestock and Veterinary Medicine, 2022, 52(8):20-23.(in Chinese)
[16] 吴健, 张国梁, 于洪春, 等.中国草原红牛(吉林系)导入红安格斯牛基因后发育性能研究[J].科技创新导报, 2016, 13(34):210-211. WU J, ZHANG G L, YU H C, et al.Development performance of Chinese grassland red cattle (Jilin line) after introducing Red Angus cattle gene[J].Science and Technology Innovation Guide, 2016, 13(34):210-211.(in Chinese)
[17] YUAN Y, BAYER P, BATLEY J, et al.Current status of structural variation studies in plants[J].Plant Biotechnology Journal, 2021, 19(11):2153-2163.
[18] CHOUAIB M, MANON M S R, AMANDINE D, et al.Development of molecular markers linked to disease resistance genes in common bean based on whole genome sequence[J].Plant Science, 2016, 242:351-357.
[19] 刘安.运用比较基因组方法研究与猪经济性状相关的InDel[D].武汉:华中农业大学, 2016. LIU A.Using comparative genomic methods to study InDel related to economic traits in pigs[D].Wuhan:Huazhong Agricultural University, 2016.(in Chinese)
[20] LI H, DURBIN R.Fast and accurate short read alignment with Burrows-Wheeler transform[J].Bioinformatics, 2009, 25(14):1754-1760.
[21] MCKENNA A, HANNA M, BANKS E, et al.The genome analysis toolkit:A mapreduce framework for analyzing next-generation DNA sequencing data[J].Genome Research, 2010, 20(9):1297-1303.
[22] WANG K, LI M, HAKONARSON H.ANNOVAR:Functional annotation of genetic variants from high-throughput sequencing data[J].Nucleic Acids Research, 2010, 38(16):e164.
[23] DANECEK P, AUTON A, ABECASIS G, et al.The variant call format and VCFtools[J].Bioinformatics, 2011, 27(15):2156-2158.
[24] CHANG C C, CHOW C C, TELLIER L C, et al.Second-generation PLINK:Rising to the challenge of larger and richer datasets[J].Gig Ascience, 2015, 4(1):7.
[25] HU Z L, PARK C A, REECY J M.Building a livestock genetic and genomic information knowledgebase through intergrative developments of Animal QTLdb and CorrDB[J].Nucleic Acids Research, 2019, 47(D1):D701-D710.
[26] WRIGH S.Evolution and the Genetics of Populations[M].Chicago:The University of Chicago Press, 1978.
[27] 杜彩虹.岷县黑裘皮羊全基因组遗传变异挖掘及选择信号分析[D].兰州:兰州大学, 2022. DU C H.Mining of genome wide genetic variation and selection signal analysis of Minxian Black fur sheep[D].Lanzhou:Lanzhou University, 2022.(in Chinese)
[28] LIAO X P, PENG F, FORNI S, et al.Whole genome sequencing of Gir cattle for identifying polymorphisms and loci under selection[J].Genome, 2013, 56(10):592-598.
[29] SEO M, LEE H J, KIM K, et al.Characterizing milk production related genes in Holstein using RNA-Seq[J].Asian-Australasian Journal of Animal Sciences, 2016, 29(3):343-351.
[30] RAZIDLO G L, JOHNSON H J, STOEGERT S M, et al.KSR1 is required for cell cycle reinitiation following DNA damage[J].The Journal of Biological Chemistry, 2009, 284(11):6705-6715.
[31] JIANG M Z, XU S, BAI M, et al.The emerging role of MEIS1 in cell proliferation and differentiation[J].American Journal of Physiology.Cell Physiology, 2020, 320(3):C264-C269.
[32] PRULLER J, FIIGEAC N, ZAMMIT P S.DVL1 and DVL3 require nuclear localisation to regulate proliferation in human myoblasts[J].Scientific Reports, 2022, 12(1):8388.
[33] YU Z Y, CHENG C J, LIU Y, et al.Neuroglobin promotes neurogenesis through Wnt signaling pathway[J].Cell Death Disease, 2018, 9(10):945.
[34] LI X Y, JI G Y, ZHOU J, et al.PCGF1 regulates early neural tube development through histone methylation in zebrafish[J].Frontiers Cell Developmental Biology, 2021, 8:581636.
[35] YUN Y, WUKUI Z, YIKAI H, et al.Loss of polycomb group protein PCGF1 severely compromises proper differentiation of embryonic stem cells[J].Scientific Reports, 2017, 7(1):46276.
[36] DUPRET B, VOLKEL P, LE B X, et al.The polycomb group protein PCGF1 is dispensable in zebrafish but involved in early growth and aging[J].PLoS One, 2016, 11(7):e0158700.
[37] LI Y, YANG S, YANG S.Rb1 negatively regulates bone formation and remodeling through inhibiting transcriptional regulation of YAP in Glut1 and OPG expression and glucose metabolism in male mice[J].Molecular Metabolism, 2022, 66:101630.
[38] WEI E, FANG X, JIA P, et al.Ginsenoside Rb1 alleviates lipopolysaccharide-induced inflammatory injury by downregulating miR-222 in WI-38 cells[J].Cell Transplantation, 2021, 30:9636897211002787.
[39] GO G Y, JO A, SEO D W, et al.Ginsenoside Rb1 and Rb2 upregulate Akt/mTOR signaling-mediated muscular hypertrophy and myoblast differentiation[J].Journal of Ginseng Research, 2020, 44(3):435-441.
[40] LI X J, ZHU F, LI B, et al.Recombinant human regenerating gene 4 attenuates the severity of osteoarthritis by promoting the proliferation of articular chondrocyte in an animal model[J].Current Molecular Pharmacology, 2022, 15(4):693-699.
[41] VON MALTZAHN M J, BENTZINGER C F, RUDNICKI M A.Wnt7a-Fzd7 signalling directly activates the Akt/mTOR anabolic growth pathway in skeletal muscle[J].Nature Cell Biology, 2011, 14(2):186-191.
[42] GUMBER D, DO M, SURESH K N, et al.Selective activation of FZD7 promotes mesendodermal differentiation of human pluripotent stem cells[J].eLife, 2020, 9:e63060.
[43] FLORIOTl S, DUCHESNE A, GROHS C, et al.A missense mutation in the FZD7 gene is associated with dilution of the red areas of the coat in Montbéliarde cattle[J].Animal Genetics, 2021, 52(3):351-355.
[44] PARVATHANENI S, LI Z, SACKS D B.Calmodulin influences MAPK signaling by binding KSR1[J].Journal of Biological Chemistry, 2021, 296:100577.
[45] TAYE M, KIM J, YOON S H, et al.Whole genome scan reveals the genetic signature of African Ankole cattle breed and potential for higher quality beef[J].BMC Genetics, 2017, 18(1):11.
[46] ZHAO Z, ZHU L, XING Y, et al.Praja2 suppresses the growth of gastric cancer by ubiquitylation of KSR1 and inhibiting MEK-ERK signal pathways[J].Aging, 2021, 13(3):3886-3897.
[47] LONG F, ZHANG Z, CHEN J, et al.The role of BBS2 in regulating adipogenesis and the association of its sequence variants with meat quality in Qinchuan cattle[J].Genomics, 2022, 114(4):110416.
[48] HAN Y J, CHEN Y, LIU Y, et al.Sequence variants of the LCORL gene and its association with growth and carcass traits in Qinchuan cattle in China[J].Journal of Genetics, 2017, 96(1):9-17.