Wang Zezhao, Ma Haoran, Xu Lei, Zhu Bo, Liu Ying, Bordbar Farhad, Chen Yan, Zhang Lupei, Gao Xue, Gao Huijiang, Zhang Shengli, Xu Lingyang, Li Junya
Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China.
Animals (Basel). 2019 May 30;9(6):296. doi: 10.3390/ani9060296.
Selective breeding can lead to genetic diversity and diverse phenotypes in farm animals. Analysis of the genomic regions under selection can provide important insights into the genetic basis of complex traits. In this study, a high-density SNP array was used for analysis of genome selection signatures in Chinese Wagyu cattle. In total, we obtained 478,903 SNPs and 24,820 no-overlap regions for |iHS| (integrated haplotype score) estimations. Under the threshold of the top 1%, 239 regions were finally identified as candidate selected regions and 162 candidate genes were found based on the UMD3.1 genome assembly. These genes were reported to be associated with fatty acids, such as Bos taurus nitric oxide synthase 1 adaptor protein (), Bos taurus hydroxysteroid 17-beta dehydrogenase 7 (), Bos taurus WD repeat domain 7 (WDR7), Bos taurus ELOVL fatty acid elongase 2 (), Bos taurus calpain 1 (), Bos taurus parkin RBR E3 ubiquitin protein ligase (, also known as ), Bos taurus mitogen-activated protein kinase kinase 6 (), meat quality, including Bos taurus ADAM metallopeptidase domain 12 (), Bos taurus 5'-aminolevulinate synthase 1 (), Bos taurus small integral membrane protein 13 () and Bos taurus potassium two pore domain channel subfamily K member 2 (), growth, and developmental traits, such as Bos taurus insulin like growth factor 2 receptor (), Bos taurus RAR related orphan receptor A (), Bos taurus fibroblast growth factor 14 (), Bos taurus paired box 6 () and Bos taurus LIM homeobox 6 (). In addition, we identified several genes that are associated with body size and weight, including Bos taurus sorting nexin 29 (), Bos taurus zinc finger imprinted 2 (), Bos taurus family with sequence similarity 110 member A (FAM110A), immune system, including Bos taurus toll like receptor 9 (), Bos taurus TAFA chemokine like family member 1 (), Bos taurus glutathione peroxidase 8 (putative) (), Bos taurus interleukin 5 (), Bos taurus PR domain containing 9 (), Bos taurus glutamate ionotropic receptor kainate type subunit 2 () and feed intake efficiency, Bos taurus sodium voltage-gated channel alpha subunit 9 (), Bos taurus relaxin family peptide/INSL5 receptor 4 (), Bos taurus RNA polymerase II associated protein 3 (). Moreover, four GO terms of biological regulation (GO:0009987, GO:0008152) and metabolic process (GO:0003824, GO:0005488) were found based on these genes. In addition, we found that 232 candidate regions (~18 Mb) overlapped with the Quantitative trait loci (QTL)regions extracted from cattle QTLdb. Our findings imply that many genes were selected for important traits in Chinese Wagyu cattle. Moreover, these results can contribute to the understanding of the genetic basis of the studied traits during the formation of this population.
选择育种可以导致农场动物的遗传多样性和多样的表型。对选择下的基因组区域进行分析可以为复杂性状的遗传基础提供重要见解。在本研究中,使用高密度SNP阵列分析中国和牛的基因组选择特征。总共,我们获得了478,903个SNP和24,820个无重叠区域用于|iHS|(整合单倍型分数)估计。在顶部1%的阈值下,最终确定了239个区域为候选选择区域,并基于UMD3.1基因组组装发现了162个候选基因。据报道,这些基因与脂肪酸相关,如牛一氧化氮合酶1衔接蛋白()、牛17-β-羟基类固醇脱氢酶7()、牛WD重复结构域7(WDR7)、牛ELOVL脂肪酸延长酶2()、牛钙蛋白酶1()、牛帕金RBR E3泛素蛋白连接酶(,也称为)、牛丝裂原活化蛋白激酶激酶6();肉质相关,包括牛ADAM金属蛋白酶结构域12()、牛5'-氨基乙酰丙酸合酶1()、牛小整合膜蛋白13()和牛钾双孔结构域通道亚家族K成员2();生长和发育性状,如牛胰岛素样生长因子2受体()、牛RAR相关孤儿受体A()、牛成纤维细胞生长因子14()、牛配对盒6()和牛LIM同源盒6()。此外,我们鉴定了几个与体型和体重相关的基因,包括牛分选衔接蛋白29()、牛锌指印记2()、牛序列相似性家族110成员A(FAM110A);免疫系统相关,包括牛Toll样受体9()、牛TAFA趋化因子样家族成员1()、牛谷胱甘肽过氧化物酶8(推定)()、牛白细胞介素5()、牛含PR结构域9()、牛谷氨酸离子型受体红藻氨酸型亚基2()以及采食量效率相关,牛钠电压门控通道α亚基9()、牛松弛素家族肽/胰岛素样肽5受体4()、牛RNA聚合酶II相关蛋白3()。此外,基于这些基因发现了生物调节(GO:0009987,GO:0008152)和代谢过程(GO:0003824,GO:0005488)的四个GO术语。另外,我们发现232个候选区域(约18 Mb)与从牛QTL数据库中提取的数量性状位点(QTL)区域重叠。我们的研究结果表明,在中国和牛中许多基因被选择用于重要性状。此外,这些结果有助于理解该群体形成过程中所研究性状的遗传基础。