RESUMEN
Body conformation traits assessed based on visual scores are widely used in Zebu cattle breeding programs. The aim of this study was to identify genomic regions and biological pathways associated with body conformation (CONF), finishing precocity (PREC), and muscling (MUSC) in Nellore cattle. The measurements based on visual scores were collected in 20,807 animals raised in pasture-based systems in Brazil. In addition, 2775 animals were genotyped using a 35 K SNP chip, which contained 31,737 single nucleotide polymorphisms after quality control. Single-step GWAS was performed using the BLUPF90 software while candidate genes were identified based on the Ensembl Genes 69. PANTHER and REVIGO platforms were used to identify key biological pathways and STRING to create gene networks. Novel candidate genes were revealed associated with CONF, including ALDH9A1, RXRG, RAB2A, and CYP7A1, involved in lipid metabolism. The genes associated with PREC were ELOVL5, PID1, DNER, TRIP12, and PLCB4, which are related to the synthesis of long-chain fatty acids, lipid metabolism, and muscle differentiation. For MUSC, the most important genes associated with muscle development were SEMA6A, TIAM2, UNC5A, and UIMC1. The polymorphisms identified in this study can be incorporated in commercial genotyping panels to improve the accuracy of genomic evaluations for visual scores in beef cattle.
RESUMEN
Genome-wide association study (GWAS) is a powerful tool to identify candidate genes and genomic regions underlying key biological mechanisms associated with economically important traits. In this context, the aim of this study was to identify genomic regions and metabolic pathways associated with backfat thickness (BFT) and rump fat thickness (RFT) in Nellore cattle, raised in pasture-based systems. Ultrasound-based measurements of BFT and RFT (adjusted to 18 months of age) were collected in 11,750 animals, with 39,903 animals in the pedigree file. Additionally, 1,440 animals were genotyped using the GGP-indicus 35K SNP chip, containing 33,623 SNPs after the quality control. The single-step GWAS analyses were performed using the BLUPF90 family programs. Candidate genes were identified through the Ensembl database incorporated in the BioMart tool, while PANTHER and REVIGO were used to identify the key metabolic pathways and gene networks. A total of 18 genomic regions located on 10 different chromosomes and harbouring 23 candidate genes were identified for BFT. For RFT, 22 genomic regions were found on 14 chromosomes, with a total of 29 candidate genes identified. The results of the pathway analyses showed important genes for BFT, including TBL1XR1, AHCYL2, SLC4A7, AADAT, VPS53, IDH2 and ETS1, which are involved in lipid metabolism, synthesis of cellular amino acids, transport of solutes, transport between Golgi Complex membranes, cell differentiation and cellular development. The main genes identified for RFT were GSK3ß, LRP1B, EXT1, GRB2, SORCS1 and SLMAP, which are involved in metabolic pathways such as glycogen synthesis, lipid transport and homeostasis, polysaccharide and carbohydrate metabolism. Polymorphisms located in these candidate genes can be incorporated in commercial genotyping platforms to improve the accuracy of imputation and genomic evaluations for carcass fatness. In addition to uncovering biological mechanisms associated with carcass quality, the key gene pathways identified can also be incorporated in biology-driven genomic prediction methods.