Buzanskas Marcos Eli, Ventura Ricardo Vieira, Seleguim Chud Tatiane Cristina, Bernardes Priscila Arrigucci, Santos Daniel Jordan de Abreu, Regitano Luciana Correia de Almeida, Alencar Maurício Mello de, Mudadu Maurício de Alvarenga, Zanella Ricardo, da Silva Marcos Vinícius Gualberto Barbosa, Li Changxi, Schenkel Flavio Schramm, Munari Danísio Prado
Departamento de Ciências Exatas, Universidade Estadual Paulista (Unesp), Faculdade de Ciências Agrárias e Veterinárias, Jaboticabal, São Paulo, Brazil.
Beef Improvement Opportunities (BIO), Guelph, Ontario, Canada.
PLoS One. 2017 Feb 9;12(2):e0171660. doi: 10.1371/journal.pone.0171660. eCollection 2017.
The aim of this study was to evaluate the level of introgression of breeds in the Canchim (CA: 62.5% Charolais-37.5% Zebu) and MA genetic group (MA: 65.6% Charolais-34.4% Zebu) cattle using genomic information on Charolais (CH), Nelore (NE), and Indubrasil (IB) breeds. The number of animals used was 395 (CA and MA), 763 (NE), 338 (CH), and 37 (IB). The Bovine50SNP BeadChip from Illumina panel was used to estimate the levels of introgression of breeds considering the Maximum likelihood, Bayesian, and Single Regression method. After genotype quality control, 32,308 SNPs were considered in the analysis. Furthermore, three thresholds to prune out SNPs in linkage disequilibrium higher than 0.10, 0.05, and 0.01 were considered, resulting in 15,286, 7,652, and 1,582 SNPs, respectively. For k = 2, the proportion of taurine and indicine varied from the expected proportion based on pedigree for all methods studied. For k = 3, the Regression method was able to differentiate the animals in three main clusters assigned to each purebred breed, showing more reasonable according to its biological viewpoint. Analyzing the data considering k = 2 seems to be more appropriate for Canchim-MA animals due to its biological interpretation. The usage of 32,308 SNPs in the analyses resulted in similar findings between the estimated and expected breed proportions. Using the Regression approach, a contribution of Indubrasil was observed in Canchim-MA when k = 3 was considered. Genetic parameter estimation could account for this breed composition information as a source of variation in order to improve the accuracy of genetic models. Our findings may help assemble appropriate reference populations for genomic prediction for Canchim-MA in order to improve prediction accuracy. Using the information on the level of introgression in each individual could also be useful in breeding or crossing design to improve individual heterosis in crossbred cattle.
本研究的目的是利用夏洛来牛(CH)、内洛尔牛(NE)和印杜布拉西尔牛(IB)品种的基因组信息,评估坎辛(CA:62.5%夏洛来牛-37.5%瘤牛)和MA遗传组(MA:65.6%夏洛来牛-34.4%瘤牛)牛品种的基因渗入水平。所使用的动物数量分别为395头(CA和MA)、763头(NE)、338头(CH)和37头(IB)。使用Illumina面板的牛50SNP基因芯片,采用最大似然法、贝叶斯法和单回归法来估计品种的基因渗入水平。经过基因型质量控制后,分析中考虑了32308个单核苷酸多态性(SNP)。此外,还考虑了三个阈值来剔除连锁不平衡高于0.10、0.05和0.01的SNP,分别得到15286个、7652个和1582个SNP。对于k = 2,所有研究方法中,牛磺酸和印度牛的比例与基于系谱的预期比例不同。对于k = 3,回归法能够将动物分为三个主要聚类,分别对应每个纯种品种,从生物学角度来看更合理。考虑到k = 2来分析数据似乎更适合坎辛-MA动物,因为其具有生物学解释。分析中使用32308个SNP得到的估计品种比例与预期比例之间的结果相似。采用回归方法,当考虑k = 3时,在坎辛-MA中观察到了印杜布拉西尔牛的贡献。遗传参数估计可以将这种品种组成信息作为变异来源,以提高遗传模型的准确性。我们的研究结果可能有助于为坎辛-MA组装合适的基因组预测参考群体,以提高预测准确性。利用每个个体的基因渗入水平信息也可能有助于育种或杂交设计,以提高杂种牛的个体杂种优势。