Signatures of selection in Nelore cattle revealed by whole-genome sequencing data

Abstract

Nelore cattle breed was farmed worldwide due to its economic importance in the beef market and adaptation to the tropics. In Brazil, purebred Nelore animals (PO) receive a certificate from the breeders' association based on the animal's genealogy and morphological characterization. The top 20 to 30% of the superior animals are eligible to receive the Special Certificate of Identification and Production (CEIP), meaning animals from this category were selected and evaluated in a breeding program to improve economically important traits. We used whole-genome sequencing and approaches based on haplotype differentiation and allelic differentiation to detect regions of selection signatures in Nelore cattle by comparing animals from PO and CEIP categories. From a total of 150 animals, a hierarchical clustering analysis was performed to choose the more unrelated animals from each category (16 PO and 40 CEIP). The hapFLK statistic was performed, and extensions of hapFLK values were investigated considering continuous regions with significant q-values. The Weir and Cockerham's Fst estimator (wcFst) was computed using the GPAT++ software library. The total of 82,326 SNPs with hapFLK values passed the FDR control (q-value<0.05), and 718 segments were target as signatures of selection. A total of 1713 highly differentiated genomic regions were identified based on the segmentFst approach. The signatures of selection were spread across the genome. Annotation of overlapping selection signature regions between the two methods revealed 118 genes in common. A variant located within the 3' region of the BOLA-DRB3 gene was found as a promising candidate polymorphism. Within genomic regions that deserves attention, we found genes previously associated with adaptation to tropical environments (HELB), growth and navel size (HMGA2), fat deposition and domestication (IRAK3), and feed efficiency and postmortem carcass traits (GABRG3). The genes BOLA-DQA2, BOLA-DQB, BOLA-DQA5, BOLA-DQA1, BOLA-DRB3, ENSBTAG00000038397 on chromosome 23 are part of the Bovine Major Histocompatibility Complex (MHC) Class II gene family, representing good candidates for immune response and adaptation to tropical conditions. The BoLA family genes and the interaction of ROBO1 with SLIT genes appeared in the enrichment results. Genomic regions located in intronic regions were also identified and might play a regulatory role in traits under selection in PO and CEIP subpopulations. The regions here identified contribute to our knowledge regarding genes and variants that have an important role in complex traits selected in this breed.

Keywords: Adaptation; Allele frequency; Bos indicus; Genomics; Haplotypes; Selective signatures.