Assessing Genomic Diversity and Selective Pressures in Bohai Black Cattle Using Whole-Genome Sequencing Data

Simple Summary Bohai Black cattle are one of the indigenous black coat cattle breeds in China, which are famous for their excellent meat quality. Whole-genome sequencing technology has been extensively developed to study species genome genetic diversity, population structure, selection pressure, demographic events, etc. However, a limited number of studies have reported genomic diversity and selection pressures in Bohai Black cattle. The purpose of this study is to analyze population structure and genomic differences between Bohai Black cattle and five "core" cattle populations from all over the world, mainly oriented on the identification of selection signatures using whole-genome sequencing data. In addition, we identify a series of candidate genes that can potentially be related to black coat color, meat quality, immunity, and reproduction in this breed. This study provides valuable genomic resources and theoretical basis for the future breeding of Bohai Black cattle. Bohai Black cattle are one of the well-known cattle breeds with black coat color in China, which are cultivated for beef. However, no study has conducted a comprehensive analysis of genomic diversity and selective pressures in Bohai Black cattle. Here, we performed a comprehensive analysis of genomic variation in 10 Bohai Black cattle (five newly sequenced and five published) and the published whole-genome sequencing (WGS) data of 50 cattle representing five "core" cattle populations. The population structure analysis revealed that Bohai Black cattle harbored the ancestry with European taurine, Northeast Asian taurine, and Chinese indicine. The Bohai Black cattle demonstrated relatively high genomic diversity from the other cattle breeds, as indicated by the nucleotide diversity (pi), the expected heterozygosity (H-E) and the observed heterozygosity (H-O), the linkage disequilibrium (LD) decay, and runs of homozygosity (ROH). We identified 65 genes containing more than five non-synonymous SNPs (nsSNPs), and an enrichment analysis revealed the "ECM-receptor interaction" pathways associated with meat quality in Bohai Black cattle. Five methods (CLR, theta pi, F-ST, theta pi ratio, and XP-EHH) were used to find several pathways and genes carried selection signatures in Bohai Black cattle, including black coat color (MC1R), muscle development (ITGA9, ENAH, CAPG, ABI2, and ISLR), fat deposition (TBC1D1, CYB5R4, TUSC3, and EPS8), reproduction traits (SPIRE2, KHDRBS2, and FANCA), and immune system response (CD84, SLAMF1, SLAMF6, and CDK10). Taken together, our results provide a valuable resource for characterizing the uniqueness of Bohai Black cattle.