Atashi H, Chen Y, Wilmot H, Bastin C, Vanderick S, Hubin X, Gengler N
TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, 5030 Gembloux, Belgium; Department of Animal Science, Shiraz University, 71441-13131 Shiraz, Iran.
TERRA Teaching and Research Center, Gembloux Agro-Bio Tech, University of Liège, 5030 Gembloux, Belgium.
J Dairy Sci. 2023 Nov;106(11):7816-7831. doi: 10.3168/jds.2022-23206. Epub 2023 Aug 9.
This study aimed to perform genome-wide association study to identify genomic regions associated with milk production and cheese-making properties (CMP) in Walloon Holstein cows. The studied traits were milk yield, fat percentage, protein percentage, casein percentage (CNP), calcium content, somatic cell score (SCS), coagulation time, curd firmness after 30 min from rennet addition, and titratable acidity. The used data have been collected from 2014 to 2020 on 78,073 first-parity (485,218 test-day records), 48,766 second-parity (284,942 test-day records), and 21,948 third-parity (105,112 test-day records) Holstein cows distributed in 671 herds in the Walloon Region of Belgium. Data of 565,533 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA) of 6,617 animals (1,712 males), were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 50 consecutive SNPs (with an average size of ∼216 KB) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for positional candidate genes. Heritability estimates for the studied traits ranged from 0.10 (SCS) to 0.53 (CNP), 0.10 (SCS) to 0.50 (CNP), and 0.12 (SCS) to 0.49 (CNP) in the first, second, and third parity, respectively. Genome-wide association analyses identified 6 genomic regions (BTA1, BTA14 [4 regions], and BTA20) associated with the considered traits. Genes including the SLC37A1 (BTA1), SHARPIN, MROH1, DGAT1, FAM83H, TIGD5, MROH6, NAPRT, ADGRB1, GML, LYPD2, JRK (BTA14), and TRIO (BTA20) were identified as positional candidate genes for the studied CMP. The findings of this study help to unravel the genomic background of a cow's ability for cheese production and can be used for the future implementation and use of genomic evaluation to improve the cheese-making traits in Walloon Holstein cows.
本研究旨在进行全基因组关联研究,以确定比利时瓦隆地区荷斯坦奶牛中与产奶量和奶酪制作特性(CMP)相关的基因组区域。所研究的性状包括产奶量、脂肪百分比、蛋白质百分比、酪蛋白百分比(CNP)、钙含量、体细胞评分(SCS)、凝固时间、添加凝乳酶30分钟后的凝乳硬度以及可滴定酸度。所使用的数据收集于2014年至2020年,涉及分布在比利时瓦隆地区671个牛群中的78,073头头胎(485,218条测定日记录)、48,766头二胎(284,942条测定日记录)和21,948头三胎(105,112条测定日记录)的荷斯坦奶牛。使用了位于6,617头动物(1,712头雄性)的29条牛(Bos taurus)常染色体(BTA)上的565,533个单核苷酸多态性(SNP)数据。通过贝叶斯吉布斯采样方法,使用随机回归测定日模型来估计遗传参数。采用单步基因组最佳线性无偏预测(BLUP)方法估计SNP效应值。计算了由50个连续SNP组成的窗口(平均大小约为216 KB)所解释的总加性遗传方差比例,并使用占总加性遗传方差至少1.0%的区域来搜索位置候选基因。所研究性状的遗传力估计值在头胎、二胎和三胎中分别为0.10(SCS)至0.53(CNP)、0.10(SCS)至0.50(CNP)和0.12(SCS)至0.49(CNP)。全基因组关联分析确定了6个与所考虑性状相关的基因组区域(BTA1、BTA14[4个区域]和BTA20)。包括溶质载体家族37成员1(SLC37A1,位于BTA1)、含SHANK相关蛋白1(SHARPIN)、MROH1、二酰甘油酰基转移酶1(DGAT1)、家族性黑素瘤83H(FAM83H)、TIGD5、MROH6、烟酰胺磷酸核糖转移酶(NAPRT)、黏附分子类G蛋白偶联受体B1(ADGRB1)、GML、淋巴细胞抗原样蛋白D2(LYPD2)、Jrk(位于BTA14)和TRIO(位于BTA20)等基因被确定为所研究CMP的位置候选基因。本研究结果有助于揭示奶牛奶酪生产能力的基因组背景,并可用于未来实施和利用基因组评估来改善瓦隆地区荷斯坦奶牛的奶酪制作性状。