Castro, Letícia Mendes de, Magnabosco, Cláudio Ulhôa, Rosa, Guilherme Jordão de Magalhães, Lopes, Fernando Brito, Magnabosco, Cláudio de Ulhoa, Rey, Fernando Sebastian Baldi, Costa, Marcos Fernando Oliveira e, Reis, Angela Adamski da Silva, and Pádua, João Teodoro de
O Brasil tem um dos maiores rebanhos bovinos comerciais do mundo, mas a qualidade da carne é altamente variável. O rebanho nacional é em grande parte composto de raças Bos indicus, que em geral têm carne menos macia do que o gado Bos taurus, diminuindo o valor do produto. Objetivou-se nesse estudo identificar regiões genômicas e vias biológicas relevantes associadas com a maciez da carne em bovinos da raça Nelore Mocho. Além disso, objetivou-se também avaliar diferentes protocolos de controle de qualidade dos SNPs e as possíveis influências nos resultados de GWAS. Os dados consistiram em valores de WBSF do músculo Longissimus dorsi, após maturação de sete dias, de 427 animais Nelore Mocho. Os animais foram genotipados em marcadores SNP Illumina BovineHD Beadchip (777k) ou Chip GGP-Indicus (77k). Todos os SNPs passaram por um Call Rate de 90% para posterior imputação utilizando o software FImput. Para averiguar os diferentes protocolos de qualidade e suas influências no GWAS foram utilizados 590.915 marcadores. Os efeitos dos diferentes QCs foram verificados utilizando 16 protocolos com três limiares para MAF (MAF < 0,01;< 0,05 e < 0,10) e HWP (p < 0,01; < 0,0001 e < 0,00001) e suas possíveis combinações. Os GWASs foram realizados utilizando método P3D/EMMAx com os marcadores restantes de cada QC. No GWAS realizado para posterior análise das vias utilizou-se 369.007 marcadores após a exclusão de SNPs baseada nos filtros Call Rate < 90%, HWP p < 0,01 e MAF < 0,01. Grupo de abate e sexo foram incluídos no modelo como efeitos fixos. Marcadores significativos (p < 0,0001) foram localizados em todas as análises, e os cromossomos com maior quantidade de SNPs significativos dos diferentes QCs foram 3, 17, 20, 21, 25 e 27. No estudo de vias foram localizados SNPs significativos nos cromossomos 3, 13, 17, 20, 21 e 25, que explicaram maior proporção da variação, indicando que existem QTLs associados à maciez da carne nessas regiões do genoma. As análises dos diferentes QCs evidenciaram efeito do controle de qualidade dos SNPs sobre o GWAS e o filtro para MAF influenciou de maneira mais ampla os resultados. Foi realizada uma análise de enriquecimento de vias baseando-se nos SNPs do GWAS, utilizando o procedimento FatiGO. Apenas os genes com no mínimo um SNP significativo (p < 0,01) foram considerados. Foram utilizados 22.365 genes anotados, incluindo 1.010 genes significativos. Um total de 22 termos GO e duas entradas IP foram consideradas enriquecidas com genes significamente associados com a maciez da carne. Várias dessas categorias funcionais como atividade da proteína tirosina quinase e serina/treonina quinase, ligantes ao íon cálcio e fatores de crescimento, podem estar relacionadas com WBFS em bovinos da raça Nelore Mocho. Estes resultados ajudam elucidar as vias relacionadas com essa característica de extrema importância econômica para o Brasil, já que o Nelore é a raça de gado de corte dominante no país. Brazil has one of the largest commercial cattle herds worldwide, but its meat quality is highly variable. The national herd is largely composed of Bos indicus breeds, which in general have less tender meat than Bos taurus cattle, decreasing the product value. This study was carried out to identify relevant regions and biological pathways associated with meat tenderness in Polled Nellore cattle. It was also aimed to evaluate the effect of different quality control protocols in GWAS for meat tenderness in Polled Nellore cattle. The data consisted of Warner-Bratzler shear force (WBSF) values of Longissimus muscle after 7 days of ageing, from 427 Polled Nellore animals. The animals were genotyped using either the Illumina BovineHD Beadchip (777k) or the GGP-Indicus Chip (77k). SNPs were excluded when Call Rate < 90%, then the imputation from the GGP to the HD Chip was performed using the FImput’s software. To study the different quality control protocols and their influence in GWAS, 590,915 markers were used. The effect of different QCs were verified using 16 protocols with three thresholds for MAF (MAF < 0.01; < 0.05 and < 0.10) and HWP (p < 0.01; < 0.0001 and < 0.00001) and their possible combinations. GWASs were performed using the PD3/EMMAx method with the remaining markers of each QC. For GWAS performed for pathway analysis, 369,007 markers were used after SNPs were excluded when Call rate < 90%, HWP p < 0.01 and MAF < 0.01. Group of slaughter and sex were included as fixed effects. Significant markers (p < 0.0001) were found in all analysis, in which the chromosomes with more significant SNPs of the different QCs were 3, 17, 20, 21, 25 and 27, and in the pathway study were located on chromosomes 3, 13, 17, 20, 21 and 25 explaining great proportion of variation, indicating possible QTLs associated with meat tenderness in those genomic regions. The analyses of different QCs showed that there is an effect of quality control over GWAS, and the filter for MAF influenced the results more broadly. A pathway enrichment analysis based on SNPs from GWAS was performed using FatiGO’s procedure. 22,365 annotated genes, including 1,010 significant genes were used. Thus, 22 GO terms and two IP entries were deemed enriched. Several of these functional categories, such as protein tyrosine and serine/threonine kinase activity, calcium ion binding and growth factors can be related to WBSF in Polled Nellore cattle. These results help to elucidate the metabolic pathways related to this trait, which is of extreme economic and social importance to Brazil as Nellore is the dominant beef cattle breed in the country. Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES