Assessing genomic diversity and signatures of selection in Jiaxian Red cattle using whole-genome sequencing data

Genome-wide detection of runs of homozygosity and heterozygosity in Tunchang pigs

Tunchang pigs, mainly distributed throughout Hainan Province of China, are well-known for their superior meat quality, crude feed tolerance, and adaptability to high temperatures and humidity. Runs of homozygosity (ROH) can provide valuable information about the inbreeding coefficient in individuals and selection signals that may reveal candidate genes associated with key functional traits. Runs of heterozygosity (ROHet) are commonly associated with balance selection, which can help us understand the adaptive evolutionary history of domestic animals. In this study, we investigated ROHs and ROHets in 88 Tunchang pigs. We also compared the estimates of inbreeding coefficients in individuals calculated based on four methods. In summary, we detected a total of 16 ROH islands in our study, and 100 genes were found within ROH regions. These genes were correlated with economically important traits such as reproduction (e.g., SERPIND1, HIRA), meat quality (e.g., PI4KA, TBX1), immunity (e.g., ESS2, RANBP1), adaption to heat stress (TXNRD2 and DGCR8), and crude food tolerance (TRPM6). Moreover, we discovered 18 ROHet islands harbouring genes associated with reproduction (e.g., ARHGEF12, BMPR2), immune system (e.g., BRD4, DNMT3B). These findings may help us design effective breeding and conservation strategies for this unique breed.

Whole-genome resequencing deciphers patterns of genetic diversity, phylogeny, and evolutionary dynamics in Kashmir cattle

Kashmir cattle, which were kept by local pastoralists for centuries, are exceptionally resilient and adaptive to harsh environments. Despite its significance, the genomic characteristics of this cattle breed remain elusive. This study utilized whole genome sequences of Kashmir cattle (n = 20; newly sequenced) alongside published whole genomes of 32 distinct breeds and seven core cattle populations (n = 135). The analysis identified ~25.87 million biallelic single nucleotide polymorphisms in Kashmir cattle, predominantly in intergenic and intron regions. Population structure analyses revealed distinct clustering patterns of Kashmir cattle with proximity to the South Asian, African and Chinese indicine cattle populations. Genetic diversity analysis of Kashmir cattle demonstrated lower inbreeding and greater nucleotide diversity than analyzed global breeds. Homozygosity runs indicated less consanguineous mating in Kashmir cattle compared with European taurine breeds. Furthermore, six selection sweep detection methods were used within Kashmir cattle and other cattle populations to identify genes associated with vital traits, including immunity (BOLA-DQA5, BOLA-DQB, TNFAIP8L, FCRL4, AOAH, HIF1AN, FBXL3, MPEG1, CDC40, etc.), reproduction (GOLGA4, BRWD1, OSBP2, LEO1 ADCY5, etc.), growth (ADPRHL1, NRG2, TCF12, TMOD4, GBP4, IGF2, RSPO3, SCD, etc.), milk composition (MRPS30 and CSF1) and high-altitude adaptation (EDNRA, ITPR2, AGBL4 and SCG3). These findings provide essential genetic insights into the characteristics and establish the foundation for the scientific conservation and utilization of Kashmir cattle breed.

Identification of Heilongjiang crossbred beef cattle pedigrees and reveals functional genes related to economic traits based on whole-genome SNP data

Introduction: To enhance the beef cattle industry, Heilongjiang Province has developed a new Crossbred beef cattle variety through crossbreeding with exotic commercial breeds. This new variety exhibits relatively excellent meat quality, and efficient reproductive performance, catering to market demands. Method: This study employed whole genome resequencing technology to analyze the genetic pedigree and diversity of 19 Heilongjiang Crossbred beef cattle, alongside 59 published genomes from East Asian, Eurasian, and European taurine cattle as controls. In addition, genes related to production traits were also searched by identifying Runs of Homozygosity (ROH) islands and important fragments from ancestors. Results: A total of 14,427,729 biallelic SNPs were discovered, with the majority located in intergenic and intron regions and a small percentage in exon regions, impacting protein function. Population genetic analyses including Principal Component Analysis (PCA), Neighbor-Joining (NJ) tree, and ADMIXTURE identified Angus, Holstein, and Mishima as the main ancestors of Crossbred beef cattle. In genetic diversity analysis, nucleotide diversity, linkage disequilibrium, and inbreeding coefficient analysis reveal that the genetic diversity of Crossbred beef cattle is at a moderate level, and a higher inbreeding coefficient indicates the need for careful breeding management. In addition, some genes related to economic traits are identified through the identification of Runs of Homozygosity (ROH) islands and important fragments from ancestors. Conclusion: This comprehensive genomic characterization supports the targeted improvement of economically important traits in Crossbred beef cattle, facilitating advanced breeding strategies.

Whole-Genome Resequencing Reveals Genetic Diversity and Growth Trait-Related Genes in Pinan Cattle

The breeding of high-quality beef cattle breeds is crucial for the development of animal husbandry, and whole-genome resequencing is widely applicated in the field of molecular breeding. Advantages in growth and reproductive traits exist in Pinan cattle compared with other cattle breeds, but there is limited research on its genomic mechanism. Using whole-genome resequencing, the genetic structure and genomic selection signatures in Pinan cattle were investigated in this study. Phylogenetic, cluster, and admixture analysis results indicated that Pinan cattle have a closer genetic relationship with Kholmogory cattle and China north cattle breeds. Through a selective sweep strategy, 207 and 54 candidate genes related to growth and reproduction and immunity, respectively, were identified in the Pinan cattle population. Given the crucial role of the glutamate-cysteine ligase catalytic (GCLC) gene in muscle antioxidative defense, the high frequency of allele T of the GCLC c.429 C>T locus in the Pinan cattle population might partially contribute to the advantages of Pinan cattle in growth performance. This study laid the foundation for the genetic improvement in Chinese local beef cattle and provide background for the studies on the growth and development of Pinan cattle.

Genome-Wide Association Studies of Growth Trait Heterosis in Crossbred Meat Rabbits

The application of heterosis can not only effectively improve the disease resistance and meat quality of livestock, but also significantly enhance the reproduction and growth of livestock and poultry. We conducted genome-wide association studies using data from F2 crossbred meat rabbits to screen out candidate genes with significant dominant effects associated with economic trait variation. High-throughput sequencing technology was used to obtain SNPs covering the whole genome to evaluate the homozygosity of the population genome, and analyze the number, length, frequency, and distribution of ROHs in the population. Candidate genes related to economic traits of meat rabbits were searched based on high-frequency ROH regions. After quality control filtering of genotype data, 380 F2 crossbred rabbits were identified with 78,579 SNPs and 42,018 ROHs on the autosomes. The fitting of the Logistic growth curve model showed that 49-day-old rabbits were a growth inflection point. Then, through genome-wide association studies, 10 SNP loci and seven growth trait candidate genes were found to be significantly related to body weight in meat rabbits at 84 days of age. In addition, we revealed the functional roles and locations of 20 candidate genes in the high-frequency ROH region associated with economic traits in meat rabbits. This study identified potential genes associated with growth and development in the high-frequency ROH region of meat rabbits. In this study, the identified candidate genes can be used as molecular markers for assisted selection in meat rabbits. At the same time, the inbreeding situation based on ROH assessment can provide reference for breeding and breeding preservation of meat rabbits.

Genome-wide detection of structural variation in some sheep breeds using whole-genome long-read sequencing data

Genomic structural variants (SVs) constitute a significant proportion of genetic variation in the genome. The rapid development of long-reads sequencing has facilitated the detection of long-fragment SVs. There is no published study to detect SVs using long-read data from sheep. We applied a long-read mapping approach to detect SVs and characterized a total of 30,771 insertions, deletions, inversions and translocations. We identified 716, 916, 842 and 303 specific SVs in Southdown sheep, Alpine merino sheep, Qilian White Tibetan sheep and Oula sheep, respectively. We annotated these SVs and found that these SV-related genes were primarily enriched in the well-established pathways involved in the regulation of the immune system, growth and development and environmental adaptability. We detected and annotated SVs based on NGS resequencing data to validate the accuracy based on third-generation detection. Moreover, five candidate SVs were verified using the PCR method in 50 sheep. Our study is the first to use a long-reads sequencing approach to construct a novel structural variation map in sheep. We have completed a preliminary exploration of the potential effects of SVs on sheep.

Genomic Characterization of Local Croatian Sheep Breeds-Effective Population Size, Inbreeding & Signatures of Selection

The Istrian (IS) and the Pag sheep (PS) are local Croatian breeds which provide significant income for the regional economy and have a cultural and traditional importance for the inhabitants. The aim of this study was to estimate some important population specific genetic parameters in IS (N = 1293) and PS (N = 2637) based on genome wide SNPs. Estimates of linkage disequilibrium effective population size (Ne) evidenced more genetic variability in PS (Ne = 838) compared to IS (Ne = 197), regardless of historical time (both recent and ancient genetic variability). The discrepancy in the recent genetic variability between these breeds was additionally confirmed by the estimates of genomic inbreeding (FROH), which was estimated to be notably higher in IS (FROH>2 = 0.062) than in PS (FROH>2 = 0.029). The average FROH2-4, FROH4-8, FROH8-16, and FROH>16 were 0.26, 1.65, 2.14, and 3.72 for IS and 0.22, 0.61, 0.75, and 1.58 for PS, thus evidencing a high contribution of recent inbreeding in the overall inbreeding. One ROH island with > 30% of SNP incidence in ROHs was detected in IS (OAR6; 34,253,440-38,238,124 bp) while there was no ROH islands detected in PS. Seven genes (CCSER1, HERC3, LCORL, NAP1L5, PKD2, PYURF, and SPP1) involved in growth, feed intake, milk production, immune responses, and resistance were associated with the found autozygosity. The results of this study represent the first comprehensive insight into genomic variability of these two Croatian local sheep breeds and will serve as a baseline for setting up the most promising strategy of genomic Optimum Contribution Selection.

Beef Cattle Genome Project: Advances in Genome Sequencing, Assembly, and Functional Genes Discovery

Beef is a major global source of protein, playing an essential role in the human diet. The worldwide production and consumption of beef continue to rise, reflecting a significant trend. However, despite the critical importance of beef cattle resources in agriculture, the diversity of cattle breeds faces severe challenges, with many breeds at risk of extinction. The initiation of the Beef Cattle Genome Project is crucial. By constructing a high-precision functional annotation map of their genome, it becomes possible to analyze the genetic mechanisms underlying important traits in beef cattle, laying a solid foundation for breeding more efficient and productive cattle breeds. This review details advances in genome sequencing and assembly technologies, iterative upgrades of the beef cattle reference genome, and its application in pan-genome research. Additionally, it summarizes relevant studies on the discovery of functional genes associated with key traits in beef cattle, such as growth, meat quality, reproduction, polled traits, disease resistance, and environmental adaptability. Finally, the review explores the potential of telomere-to-telomere (T2T) genome assembly, structural variations (SVs), and multi-omics techniques in future beef cattle genetic breeding. These advancements collectively offer promising avenues for enhancing beef cattle breeding and improving genetic traits.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Whole-genome sequencing reveals genomic diversity and selection signatures in Xia'nan cattle

BACKGROUND

The crossbreeding of specialized beef cattle breeds with Chinese indigenous cattle is a common method of genetic improvement. Xia'nan cattle, a crossbreed of Charolais and Nanyang cattle, is China's first specialized beef cattle breed with independent intellectual property rights. After more than two decades of selective breeding, Xia'nan cattle exhibit a robust physique, good environmental adaptability, good tolerance to coarse feed, and high meat production rates. This study analyzed the population genetic structure, genetic diversity, and genomic variations of Xia'nan cattle using whole-genome sequencing data from 30 Xia'nan cattle and 178 published cattle genomic data.

RESULT

The ancestry estimating composition analysis showed that the ancestry proportions for Xia'nan cattle were mainly Charolais with a small amount of Nanyang cattle. Through the genetic diversity studies (nucleotide diversity and linkage disequilibrium decay), we found that the genomic diversity of Xia'nan cattle is higher than that of specialized beef cattle breeds in Europe but lower than that of Chinese native cattle. Then, we used four methods to detect genome candidate regions influencing the excellent traits of Xia'nan cattle. Among the detected results, 42 genes (θπ and CLR) and 131 genes (FST and XP-EHH) were detected by two different detection strategies. In addition, we found a region in BTA8 with strong selection signals. Finally, we conducted functional annotation on the detected genes and found that these genes may influence body development (NR6A1), meat quality traits (MCCC1), growth traits (WSCD1, TMEM68, MFN1, NCKAP5), and immunity (IL11RA, CNTFR, CCL27, SLAMF1, SLAMF7, NAA35, and GOLM1).

CONCLUSION

We elucidated the genomic features and population structure of Xia'nan cattle and detected some selection signals in genomic regions potentially associated with crucial economic traits in Xia'nan cattle. This research provided a basis for further breeding improvements in Xia'nan cattle and served as a reference for genetic enhancements in other crossbreed cattle.

Genome-wide analysis reveals genomic diversity and signatures of selection in Qinchuan beef cattle

BACKGROUND

Indigenous Chinese cattle have abundant genetic diversity and a long history of artificial selection, giving local breeds advantages in adaptability, forage tolerance and resistance. The detection of selective sweeps and comparative genome analysis of selected breeds and ancestral populations provide a basis for understanding differences among breeds and for the identification and utilization of candidate genes. We investigated genetic diversity, population structure, and signatures of selection using genome-wide sequencing data for a new breed of Qinchuan cattle (QNC, n = 21), ancestral Qinchuan cattle (QCC, n = 20), and Zaosheng cattle (ZSC, n = 19).

RESULTS

A population structure analysis showed that the ancestry components of QNC and ZSC were similar. In addition, the QNC and ZSC groups showed higher proportions of European taurine ancestry than that of QCC, and this may explain the larger body size of QNC, approaching that of European cattle under long-term domestication and selection. A neighbor-joining tree revealed that QCC individuals were closely related, whereas QNC formed a distinct group. To search for signatures of selection in the QNC genome, we evaluated nucleotide diversity (θπ), the fixation index (FST) and Tajima's D. Overlapping selective sweeps were enriched for one KEGG pathway, the apelin signaling pathway, and included five candidate genes (MEF2A, SMAD2, CAMK4, RPS6, and PIK3CG). We performed a comprehensive review of genomic variants in QNC, QCC, and ZSC using whole-genome sequencing data. QCC was rich in novel genetic diversity, while diversity in QNC and ZSC cattle was reduced due to strong artificial selection, with divergence from the original cattle.

CONCLUSIONS

We identified candidate genes associated with production traits. These results support the success of selective breeding and can guide further breeding and resource conservation of Qinchuan cattle.

Unraveling genomic diversity and positive selection signatures of Qaidam cattle through whole-genome re-sequencing

Qaidam cattle are a typical Chinese native breed inhabiting northwest China. They bear the characteristics of high cold and roughage tolerance, low-oxygen adaptability and good meat quality. To analyze the genetic diversity of Qaidam cattle, 60 samples were sequenced using whole-genome resequencing technology, along with 192 published sets of whole-genome sequencing data of Indian indicine cattle, Chinese indicine cattle, North Chinese cattle breeds, East Asian taurine cattle, Eurasian taurine cattle and European taurine cattle as controls. It was found that Qaidam cattle have rich genetic diversity in Bos taurus, but the degree of inbreeding is also high, which needs further protection. The phylogenetic analysis, principal component analysis and ancestral component analysis showed that Qaidam cattle mainly originated from East Asian taurine cattle. Qaidam cattle had a closer genetic relationship with the North Chinese cattle breeds and the least differentiation from Mongolian cattle. Annotating the selection signals obtained by composite likelihood ratio, nucleotide diversity analysis, integrated haplotype score, genetic differentiation index, genetic diversity ratio and cross-population extended haplotype homozygosity methods, several genes associated with immunity, reproduction, meat, milk, growth and adaptation showed strong selection signals. In general, this study provides genetic evidence for understanding the germplasm characteristics of Qaidam cattle. At the same time, it lays a foundation for the scientific and reasonable protection and utilization of genetic resources of Chinese local cattle breeds, which has great theoretical and practical significance.

Study on the origin of the Baise horse based on whole-genome resequencing

The Baise horse, an indigenous horse breed mainly distributed in the Baise region of Guangxi province in southwest China, has a long history as draft animal. However, there is a lack of research regarding the origin and ancestral composition of the Baise horse. In this study, whole-genome resequencing data from 236 horses of seven Chinese indigenous horse breeds, five foreign horse breeds, and four Przewalski's horses were used to investigate the relationships between the Baise horse and other horse breeds. The results showed that foreign horse breeds had no significant impact on the formation of the Baise horse. The two southwestern horse populations, the Debao pony and the Jinjiang horse, exhibit the closest genetic affinity with the Baise horse. This is consistent with their adjacent geographical distribution. Analysis of the migration route revealed a gene flow from the Chakouyi horse into the Baise horse. In summary, our results confirm that the formation of the Baise horse did not involve participation from foreign breeds. Geographical distance emerges as a crucial factor in determining the genetic relationships with the Baise horse. Gene flows of indigenous horse breeds along ancient routes of trade activities had played a role in the formation of the Baise horse.

POLB Regulates Proliferation and Apoptosis of Bovine Primary Myocytes

DNA polymerase β (DNA polymerase beta (POLB)) belongs to a member of the DNA polymerase X family, mainly involved in various biological metabolic processes, such as eukaryotic DNA replication, DNA damage repair, gene recombination, and cell cycle regulation. In this study, the muscle development-related gene POLB was screened by selection signature and RNA-seq analysis and then validated for the proliferation and apoptosis of bovine primary myocytes. It was also found that overexpression of the POLB gene had a pro-apoptosis effect, but interfering with the expression of the gene had no significant effect on cells. Then, the analysis of related apoptotic genes revealed that POLB overexpression affected CASP9 gene expression.

Resequencing Analyses Revealed Genetic Diversity and Selection Signatures during Rabbit Breeding and Improvement

Domestication has shaped the diverse characteristics of rabbits, including coat color, fur structure, body size, and various physiological traits. Utilizing whole-genome resequencing (DNBSEQ-T7), we analyzed the genetic diversity, population structure, and genomic selection across 180 rabbits from 17 distinct breeds to uncover the genetic basis of these traits. We conducted whole-genome sequencing on 17 rabbit breeds, identifying 17,430,184 high-quality SNPs and analyzing genomic diversity, patterns of genomic variation, population structure, and selection signatures related to coat color, coat structure, long hair, body size, reproductive capacity, and disease resistance. Through PCA and NJ tree analyses, distinct clusters emerged among Chinese indigenous rabbits, suggesting varied origins and domestication histories. Selective sweep testing pinpointed regions and genes linked to domestication and key morphological and economic traits, including those affecting coat color (TYR, ASIP), structure (LIPH), body size (INSIG2, GLI3), fertility (EDNRA, SRD5A2), heat stress adaptation (PLCB1), and immune response (SEC31A, CD86, LAP3). Our study identified key genomic signatures of selection related to traits such as coat color, fur structure, body size, and fertility; these findings highlight the genetic basis underlying phenotypic diversification in rabbits and have implications for breeding programs aiming to improve productive, reproductive, and adaptive traits. The detected genomic signatures of selection also provide insights into rabbit domestication and can aid conservation efforts for indigenous breeds.

Population Structure and Selection Signal Analysis of Nanyang Cattle Based on Whole-Genome Sequencing Data

With a rich breeding history, Nanyang cattle (NY cattle) have undergone extensive natural and artificial selection, resulting in distinctive traits such as high fertility, excellent meat quality, and disease resistance. This makes them an ideal model for studying the mechanisms of environmental adaptability. To assess the population structure and genetic diversity of NY cattle, we performed whole-genome resequencing on 30 individuals. These data were then compared with published whole-genome resequencing data from 432 cattle globally. The results indicate that the genetic structure of NY cattle is significantly different from European commercial breeds and is more similar to North-Central Chinese breeds. Furthermore, among all breeds, NY cattle exhibit the highest genetic diversity and the lowest population inbreeding levels. A genome-wide selection signal analysis of NY cattle and European commercial breeds using Fst, θπ-ratio, and θπ methods revealed significant selection signals in genes associated with reproductive performance and immunity. Our functional annotation analysis suggests that these genes may be responsible for reproduction (MAP2K2, PGR, and GSE1), immune response (NCOA2, HSF1, and PAX5), and olfaction (TAS1R3). We provide a comprehensive overview of sequence variations in the NY cattle genome, revealing insights into the population structure and genetic diversity of NY cattle. Additionally, we identify candidate genes associated with important economic traits, offering valuable references for future conservation and breeding efforts of NY cattle.

Identification of genomic diversity and selection signatures in Luxi cattle using whole-genome sequencing data

OBJECTIVE

The objective of this study was to investigate the genetic diversity, population structure and whole-genome selection signatures of Luxi cattle to reveal its genomic characteristics in terms of meat and carcass traits, skeletal muscle development, body size, and other traits.

METHODS

To further analyze the genomic characteristics of Luxi cattle, this study sequenced the whole-genome of 16 individuals from the core conservation farm in Shandong region, and collected 174 published genomes of cattle for conjoint analysis. Furthermore, three different statistics (pi, Fst, and XP-EHH) were used to detect potential positive selection signatures related to selection in Luxi cattle. Moreover, gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment analyses were performed to reveal the potential biological function of candidate genes harbored in selected regions.

RESULTS

The results showed that Luxi cattle had high genomic diversity and low inbreeding levels. Using three complementary methods (pi, Fst, and XP-EHH) to detect the signatures of selection in the Luxi cattle genome, there were 2,941, 2,221 and 1,304 potentially selected genes identified, respectively. Furthermore, there were 45 genes annotated in common overlapping genomic regions covered 0.723 Mb, including PLAG1 zinc finger (PLAG1), dedicator of cytokinesis 3 (DOCK3), ephrin A2 (EFNA2), DAZ associated protein 1 (DAZAP1), Ral GTPase activating protein catalytic subunit alpha 1 (RALGAPA1), mediator complex subunit 13 (MED13), and decaprenyl diphosphate synthase subunit 2 (PDSS2), most of which were enriched in pathways related to muscle growth and differentiation and immunity.

CONCLUSION

In this study, we provided a series of genes associated with important economic traits were found in positive selection regions, and a scientific basis for the scientific conservation and genetic improvement of Luxi cattle.

The Identification of Goat KCNJ15 Gene Copy Number Variation and Its Association with Growth Traits

(1) Background: Copy number variation (CNV) is a critical component of genome structural variation and has garnered significant attention. High-throughput screening of the KCNJ15 gene has revealed a correlation between the CNV region and the growth traits of goats. We aimed to identify the CNV of the KCNJ15 gene in five goat breeds and analyze its association with growth characteristics. (2) Methods: We utilized 706 goats from five breeds: Guizhou black goat (GZB), Guizhou white goat (GZW), Bohuai goat (BH), Huai goat (HH), and Taihang goat (TH). To evaluate the number of copies of the KCNJ15 gene using qPCR, we analyzed the correlation between the CNV and growth characteristics and then used a universal linear model. The findings revealed variations in the distribution of different copy number types among the different goat breeds. (3) Results: Association analysis revealed a positive influence of the CNV in the KCNJ15 gene on goat growth. In GZB, individuals with duplication types exhibited superior performance in terms of cannon bone circumference (p < 0.05). In HH, individuals with duplication types exhibited superior performance in terms of body slanting length (p < 0.05). Conversely, normal TH demonstrated better body height and body weight (p < 0.05), while in GZW, when CN = 3, it performed better than other types in terms of body weight and chest circumference (p < 0.05). However, in BH, it had no significant effect on growth traits. (4) Conclusions: We confirmed that the CNV in the KCNJ15 gene significantly influences the growth characteristics of four distinct goat breeds. The correlation between KCNJ15 gene CNVs and goat growth traits offers valuable insights to breeders, enabling them to employ precise and efficient breeding methods that enhance livestock welfare, productivity, and overall economic benefits in the industry.

Whole-Genome Resequencing Reveals Genetic Diversity and Wool Trait-Related Genes in Liangshan Semi-Fine-Wool Sheep

Understanding the genetic makeup of local sheep breeds is essential for their scientific conservation and sustainable utilization. The Liangshan semi-fine-wool sheep (LSS), a Chinese semi-fine-wool breed renowned for its soft wool, was analyzed using whole-genome sequencing data including 35 LSS, 84 sheep from other domestic breeds, and 20 Asiatic mouflons. We investigated the genetic composition of LSS by conducting analyses of the population structure, runs of homozygosity, genomic inbreeding coefficients, and selection signature. Our findings indicated that LSS shares greater genetic similarity with Border Leicester and Romney sheep than with Tibetan (TIB), Yunnan (YNS), and Chinese Merino sheep. Genomic analysis indicated low to moderate inbreeding coefficients, ranging from 0.014 to 0.154. In identifying selection signals across the LSS genome, we pinpointed 195 candidate regions housing 74 annotated genes (e.g., IRF2BP2, BVES, and ALOX5). We also found the overlaps between the candidate regions and several known quantitative trait loci related to wool traits, such as the wool staple length and wool fiber diameter. A selective sweep region, marked by the highest value of cross-population extended haplotype homozygosity, encompassed IRF2BP2-an influential candidate gene affecting fleece fiber traits. Furthermore, notable differences in genotype frequency at a mutation site (c.1051 + 46T > C, Chr25: 6,784,190 bp) within IRF2BP2 were observed between LSS and TIB and YNS sheep (Fisher's exact test, p < 2.2 × 10-16). Taken together, these findings offer insights crucial for the conservation and breeding enhancement of LSS.

Assessing genomic diversity and selective pressures in Bashan cattle by whole-genome sequencing data

Specific ecological environments and domestication have continuously influenced the physiological characteristics of Chinese indigenous cattle. Among them, Bashan cattle belongs to one of the indigenous breeds. However, the genomic diversity of Bashan cattle is still unknown. Published whole-genome sequencing (WGS) data of 13 Bashan cattle and 48 worldwide cattle were used to investigate the genetic composition and selection characteristics of Bashan cattle. The population structure analysis revealed that Bashan cattle harbored ancestries with East Asian taurine and Chinese indicine. Genetic diversity analysis implied the relatively high genomic diversity in Bashan cattle. Through the identification of containing >5 nsSNPs or frameshift mutations genes in Bashan cattle, a large number of pathways related to sensory perception were discovered. CLR, θπ ratio, FST, and XP-EHH methods were used to detect the candidate signatures of positive selection in Bashan cattle. Among the identified genes, most of the enriched signal pathways were related to environmental information processing, biological systems, and metabolism. We mainly reported genes related to the nervous system (HCN1, KATNA1, FSTL1, GRIK2, and CPLX2), immune (CD244, SLAMF1, LY9, and CD48), and reproduction (AKR1C1, AKR1C3, AKR1C4, and TUSC3). Our findings will be significant in understanding the molecular basis underlying phenotypic variation of breed-related traits and improving productivity in Bashan cattle.

Association between the copy number variation of CCSER1 gene and growth traits in Chinese Capra hircus (goat) populations

Recently, Coiled-coil serine-rich protein 1 (CCSER1) gene is reported to be related to economic traits in livestock, and become a hotspot. In our study, we detected CCSER1 gene CNV in 693 goats from six breeds (GZB, GZW, AN, BH, HG, TH) by quantitative real-time PCR (qPCR) and the association analysis between the types of CNV and growth traits. Then, CCSER1 gene expression pattern was discovered in seven tissues from NB goats. Our results showed that the CCSER1 gene copy numbers were distributed differently in the aforementioned six breeds. The type of CCSER1 gene CNV was significantly associated with body weight and heart girth traits in GZW goat, in which individuals with deletion type were dominant in body weight trait (P < 0.05), while the normal type individuals were more advantageous in heart girth trait (P < 0.01); and there was a significant association with heart girth in TH goat (P < 0.05), which normal type was the dominant one. The expression profile revealed that CCSER1 gene has the highest level in the lung, followed by the small intestine and heart. In conclusion, our result is dedicated to an in-depth study of the novel CCSER1 gene CNV site and to provide essential information for Chinese goats molecular selective breeding in the future.

Genomic selection pressure discovery using site-frequency spectrum and reduced local variability statistics in Pakistani Dera-Din-Panah goat

BACKGROUND

Population geneticists have long sought to comprehend various selection traces accumulated in the goat genome due to natural or human driven artificial selection through breeding practices, which led the wild animals to domestication, so understanding evolutionary process may helpful to utilize the full genetic potential of goat genome.

METHODS AND RESULTS

As a step forward to pinpoint the selection signals in Pakistani Dera-Din-Panah (DDP) goat, whole-genome pooled sequencing (n = 12) was performed, and 618,236,192 clean paired-end reads were mapped against ARS1 reference goat assembly. Five different selection signature statistics were applied using four site-frequency spectrum (SFS) methods (Tajima's D ([Formula: see text]), Fay and Wu's H ([Formula: see text]), Zeng's E ([Formula: see text]), [Formula: see text]) and one reduced local variability approach named pooled heterozygosity ([Formula: see text]). The under-selection regions were annotated with significant threshold values of [Formula: see text]≥4.7, [Formula: see text]≥6, [Formula: see text]≥2.5, Pool-HMM ≥ 12, and [Formula: see text]≥5 that resulted in accumulative 364 candidate gene hits. The highest genomic selection signals were observed on Chr. 4, 6, 10, 12, 15, 16, 18, 20, and 27 and harbor ADAMTS6, CWC27, RELN, MYCBP2, FGF14, STIM1, CFAP74, GNB1, CALML6, TMEM52, FAM149A, NADK, MMP23B, OPN3, FH, MFHAS1, KLKB1, RRM1, KMO, SPEF2, F11, KIT, KMO, ERI1, ATP8B4, and RHOG genes. Next, the validation of our captured genomic hits was also performed by more than one applied statistics which harbor meat production, immunity, and reproduction associated genes to strengthen our hypothesis of under-selection traits in this Pakistani goat breed. Furthermore, common candidate genes captured by more than one statistical method were subjected to gene ontology and KEGG pathway analysis to get insights of particular biological processes associated with this goat breed.

CONCLUSION

Current perception of genomic architecture of DDP goat provides a better understanding to improve its genetic potential and other economically important traits of medium to large body size, milk, and fiber production by updating the genomic insight driven breeding strategies to boost the livestock and agriculture-based economy of the country.

Genetic Distinctness and Diversity of American Aberdeen Cattle Compared to Common Beef Breeds in the United States

American Aberdeen (AD) cattle in the USA descend from an Aberdeen Angus herd originally brought to the Trangie Agricultural Research Centre, New South Wales, AUS. Although put under specific selection pressure for yearling growth rate, AD remain genomically uncharacterized. The objective was to characterize the genetic diversity and structure of purebred and crossbred AD cattle relative to seven common USA beef breeds using available whole-genome SNP data. A total of 1140 animals consisting of 404 purebred (n = 8 types) and 736 admixed individuals (n = 10 types) was used. Genetic diversity metrics, an analysis of molecular variance, and a discriminant analysis of principal components were employed. When linkage disequilibrium was not accounted for, markers influenced basic diversity parameter estimates, especially for AD cattle. Even so, intrapopulation and interpopulation estimates separate AD cattle from other purebred types (e.g., Latter's pairwise FST ranged from 0.1129 to 0.2209), where AD cattle were less heterozygous and had lower allelic richness than other purebred types. The admixed AD-influenced cattle were intermediate to other admixed types for similar parameters. The diversity metrics separation and differences support strong artificial selection pressures during and after AD breed development, shaping the evolution of the breed and making them genomically distinct from similar breeds.

Integrated Comparative Transcriptome and circRNA-lncRNA-miRNA-mRNA ceRNA Regulatory Network Analyses Identify Molecular Mechanisms Associated with Intramuscular Fat Content in Beef Cattle

Intramuscular fat content (IMF), one of the most important carcass traits in beef cattle, is controlled by complex regulatory factors. At present, molecular mechanisms involved in regulating IMF and fat metabolism in beef cattle are not well understood. Our objective was to integrate comparative transcriptomic and competing endogenous RNA (ceRNA) network analyses to identify candidate messenger RNAs (mRNAs) and regulatory RNAs involved in molecular regulation of longissimus dorsi muscle (LDM) tissue for IMF and fat metabolism of 5 beef cattle breeds (Angus, Chinese Simmental, Luxi, Nanyang, and Shandong Black). In total, 34 circRNAs, 57 lncRNAs, 15 miRNAs, and 374 mRNAs were identified by integrating gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Furthermore, 7 key subnets with 16 circRNAs, 43 lncRNAs, 7 miRNAs, and 237 mRNAs were detected through clustering analyses, whereas GO enrichment analysis of identified RNAs revealed 48, 13, and 28 significantly enriched GO terms related to IMF in biological process, molecular function, and cellular component categories, respectively. The main metabolic-signaling pathways associated with IMF and fat metabolism that were enriched included metabolic, calcium, cGMP-PKG, thyroid hormone, and oxytocin signaling pathways. Moreover, MCU, CYB5R1, and BAG3 genes were common among the 10 comparative groups defined as important candidate marker genes for fat metabolism in beef cattle. Contributions of transcriptome profiles from various beef breeds and a competing endogenous RNA (ceRNA) regulatory network underlying phenotypic differences in IMF provided novel insights into molecular mechanisms associated with meat quality.

Local ancestry and selection in admixed Sanjiang cattle

The majority of native cattle are taurine × indicine cattle of diverse phenotypes in the central region of China. Sanjiang cattle, a typical breed in the central region, play a central role in human livelihood and have good adaptability, including resistance to dampness, heat, roughage, and disease, and are thus regarded as an important genetic resource. However, the genetic history of the successful breed remains unknown. Here, we sequenced 10 Sanjiang cattle genomes and compared them to the 70 genomes of 5 representative populations worldwide. We characterized the genomic diversity and breed formation process of Sanjiang cattle and found that Sanjiang cattle have a mixed ancestry of indicine (55.6%) and taurine (33.2%) dating to approximately 30 generations ago, which has shaped the genome of Sanjiang cattle. Through ancestral fragment inference, selective sweep and transcriptomic analysis, we identified several genes linked to lipid metabolism, immune regulation, and stress reactions across the mosaic genome of Sanjiang cattle showing an excess of taurine or indicine ancestry. Taurine ancestry might contribute to meat quality, and indicine ancestry is more conducive to adaptation to hot climate conditions, making Sanjiang cattle a valuable genetic resource for the central region of China. Our results will help us understand the evolutionary history and ancestry components of Sanjiang cattle, which will provide a reference for resource conservation and selective breeding of Chinese native cattle.

Genome-wide diversity and admixture of five indigenous cattle populations from the Tigray region of northern Ethiopia

The Tigray region, where we found around eight per cent of the indigenous cattle population of Ethiopia, is considered as the historic centre of the country, with the ancient pre-Aksumite and Aksumite civilisations in contact with the civilisations of the Fertile Crescent and the Indian subcontinent. Here, we used whole genome sequencing data to characterise the genomic diversity, relatedness, and admixture of five cattle populations (Abergelle, Arado, Begait, Erob, and Raya) indigenous to the Tigray region of Ethiopia. We detected 28 to 29 million SNPs and 2.7 to 2.9 million indels in each population, of which 7% of SNPs and 34% of indels were novel. Functional annotation of the variants showed around 0.01% SNPs and 0.22%-0.27% indels in coding regions. Enrichment analysis of genes overlapping missense private SNPs revealed 20 significant GO terms and KEGG pathways that were shared by or specific to breeds. They included important genes associated with morphology (SCN4A, TAS1R2 and KCNG4), milk yield (GABRG1), meat quality (MMRN2, VWC2), feed efficiency (PCDH8 and SLC26A3), immune response (LAMC1, PCDH18, CELSR1, TLR6 and ITGA5), heat resistance (NPFFR1 and HTR7) and genes belonging to the olfactory gene family, which may be related to adaptation to harsh environments. Tigray indigenous cattle are very diverse. Their genome-wide average nucleotide diversity ranged from 0.0035 to 0.0036. The number of heterozygous SNPs was about 0.6-0.7 times higher than homozygous ones. The within-breed average number of ROHs ranged from 777.82 to 1000.45, with the average sum of the length of ROHs ranging from 122.01 Mbp to 163.88 Mbp. The genomic inbreeding coefficients differed among animals and breeds, reaching up to 10% in some Begait and Raya animals. Tigray indigenous cattle shared a common ancestry with Asian indicine (85.6%-88.7%) and African taurine (11.3%-14.1%) cattle, with very small, if any, European taurine introgression. This study identified high within-breed genetic diversity representing an opportunity for breeding improvement programs and, also, significant novel variants that could increase the number of known cattle variants, an important contribution to the knowledge of domestic cattle genetic diversity.

Insertion of 643bp Retrotransposon Upstream of PPARγ CDS Is Associated with Backfat of Large White Pigs

PPARs are essential regulators of mammalian fatty acid and lipid metabolism. Although the effects of genetic variations, including single nucleotide polymorphisms (SNPs) in PPARs genes on the phenotype of domestic animals have been investigated, there is limited information on the impact of retrotransposon insertion polymorphisms (RIPs). In this study, a combined comparative genome and polymerase chain reaction (PCR) was used to excavate the RIPs in porcine PPARs. We also investigated the potential effects of retrotransposon insertion on phenotype and expression patterns. This study identified the two RIPs in PPARs genes, namely an ERV in intron 1 of PPARα and a combined retrotransposon in intron 2 of PPARγ, designated as PPARα-ERV-RIP and PPARγ-COM-RIP, respectively. These RIPs exhibited different distribution patterns among Chinese indigenous breeds and Western commercial breeds. Individuals with the PPARα-ERV-RIP^+/+ genotype (+/+ indicated homozygous with insertion) among Large White pigs had significantly higher (p < 0.05) corrected backfat thickness compared to those with the other two genotypes. Similarly, those with the PPARγ-COM-RIP^-/- genotype had significantly higher (p < 0.05) corrected backfat thickness than those with the other two genotypes in Large White pigs. Moreover, in 30-day-old Sujiang piglets, the PPARγ gene expression in the backfat of those with the PPARγ-COM-RIP^-/- genotype (-/- indicated homozygous without insertion) was significantly greater (p < 0.01) than those with other genotypes. The dual luciferase reporter gene assay demonstrated that the combined retrotransposon insertion significantly reduced the activity of the MYC promoter in both C2C12 and 3T3-L1 cells (p < 0.01). Therefore, the combined retrotransposon insertion could function as a repressor to decrease the expression of PPARγ, making PPARγ-COM-RIP a valuable molecular marker for assisted selection of backfat thickness in pig breeding.

Assessing Genomic Diversity and Signatures of Selection in Chinese Red Steppe Cattle Using High-Density SNP Array

Chinese Red Steppe Cattle (CRS), a composite cattle breed, is well known for its milk production, high slaughter rate, carcass traits, and meat quality. Nowadays, it is widely bred in Jilin and Hebei Province and the Inner Mongolia Autonomous region. However, the population structure and the genetic basis of prominent characteristics of CRS are still unknown. In this study, we systematically describe their population structure, genetic diversity, and selection signature based on genotyping data from 61 CRS individuals with GGP Bovine 100 K chip. The results showed that CRS cattle had low inbreeding levels and had formed a unique genetic structure feature. Using two complementary methods (including comprehensive haplotype score and complex likelihood ratio), we identified 1291 and 1285 potentially selected genes, respectively. There were 141 genes annotated in common 106 overlapping genomic regions covered 5.62 Mb, including PLAG1, PRKG2, DGAT1, PARP10, TONSL, ADCK5, and BMP3, most of which were enriched in pathways related to muscle growth and differentiation, milk production, and lipid metabolism. This study will contribute to understanding the genetic mechanism behind artificial selection and give an extensive reference for subsequent breeding.

A Genome-Wide Association Study of the Chest Circumference Trait in Xinjiang Donkeys Based on Whole-Genome Sequencing Technology

Animal genotyping by means of genome-wide association studies is important for connecting phenotypes of interest with their underlying genetics in livestock. However, the use of whole genome sequencing to investigate chest circumference (CC) in donkeys has rarely been reported. We aimed to use the genome-wide association study approach to detect significant single nucleotide polymorphisms (SNPs) and key genes associated with chest circumference traits in Xinjiang donkeys. We assessed 112 Xinjiang donkeys in this study. The chest circumference of each was measured 2 h before milking. We re-sequenced blood samples from the Xinjiang donkeys, and genome-wide association study analyses were performed using a mixed model with the PLINK, GEMMA, and REGENIE programs. We tested 38 donkeys for candidate SNPs for genome-wide association study using three software programs. Additionally, 18 SNP markers reached genome-wide significance (p < 1.61 × 10^-9). On the basis of these, 41 genes were identified. Previously proposed candidate genes for CC traits were supported by this study, including NFATC2 (Nuclear Factor of Activated T Cells 2), PROP1 (PROP Paired-Like Homeobox 1), UBB (Ubiquitin B), and HAND2 (Heart and Neural Crest Derivatives Expressed 2). These promising candidates provide a valuable resource for validating potential meat production genes and will facilitate the development of high-yielding Xinjiang donkey breeds through marker-assisted selection or gene editing.

Global dispersal and adaptive evolution of domestic cattle: a genomic perspective

Domestic cattle have spread across the globe and inhabit variable and unpredictable environments. They have been exposed to a plethora of selective pressures and have adapted to a variety of local ecological and management conditions, including UV exposure, diseases, and stall-feeding systems. These selective pressures have resulted in unique and important phenotypic and genetic differences among modern cattle breeds/populations. Ongoing efforts to sequence the genomes of local and commercial cattle breeds/populations, along with the growing availability of ancient bovid DNA data, have significantly advanced our understanding of the genomic architecture, recent evolution of complex traits, common diseases, and local adaptation in cattle. Here, we review the origin and spread of domestic cattle and illustrate the environmental adaptations of local cattle breeds/populations.

Molecular Cytogenetics in Domestic Bovids: A Review

The discovery of the Robertsonian translocation (rob) involving cattle chromosomes 1 and 29 and the demonstration of its deleterious effects on fertility focused the interest of many scientific groups on using chromosome banding techniques to reveal chromosome abnormalities and verify their effects on fertility in domestic animals. At the same time, comparative banding studies among various species of domestic or wild animals were found useful for delineating chromosome evolution among species. The advent of molecular cytogenetics, particularly the use of fluorescence in situ hybridization (FISH), has allowed a deeper investigation of the chromosomes of domestic animals through: (a) the physical mapping of specific DNA sequences on chromosome regions; (b) the use of specific chromosome markers for the identification of the chromosomes or chromosome regions involved in chromosome abnormalities, especially when poor banding patterns are produced; (c) better anchoring of radiation hybrid and genetic maps to specific chromosome regions; (d) better comparisons of related and unrelated species by comparative FISH mapping and/or Zoo-FISH techniques; (e) the study of meiotic segregation, especially by sperm-FISH, in some chromosome abnormalities; (f) better demonstration of conserved or lost DNA sequences in chromosome abnormalities; (g) the use of informatic and genomic reconstructions, in addition to CGH arrays, to predict conserved or lost chromosome regions in related species; and (h) the study of some chromosome abnormalities and genomic stability using PCR applications. This review summarizes the most important applications of molecular cytogenetics in domestic bovids, with an emphasis on FISH mapping applications.

Longitudinal genome-wide association studies of milk production traits in Holstein cattle using whole-genome sequence data imputed from medium-density chip data

Longitudinal traits, such as milk production traits in dairy cattle, are featured by having phenotypic values at multiple time points, which change dynamically over time. In this study, we first imputed SNP chip (50-100K) data to whole-genome sequence (WGS) data in a Chinese Holstein population consisting of 6,470 cows. The imputation accuracies were 0.88 to 0.97 on average after quality control. We then performed longitudinal GWAS in this population based on a random regression test-day model using the imputed WGS data. The longitudinal GWAS revealed 16, 39, and 75 quantitative trait locus regions associated with milk yield, fat percentage, and protein percentage, respectively. We estimated the 95% confidence intervals (CI) for these quantitative trait locus regions using the logP drop method and identified 581 genes involved in these CI. Further, we focused on the CI that covered or overlapped with only 1 gene or the CI that contained an extremely significant top SNP. Twenty-eight candidate genes were identified in these CI. Most of them have been reported in the literature to be associated with milk production traits, such as DGAT1, HSF1, MGST1, GHR, ABCG2, ADCK5, and CSN1S1. Among the unreported novel genes, some also showed good potential as candidate genes, such as CCSER1, CUX2, SNTB1, RGS7, OSR2, and STK3, and are worth being further investigated. Our study provided not only new insights into the candidate genes for milk production traits, but also a general framework for longitudinal GWAS based on random regression test-day model using WGS data.

Whole-genome resequencing reveals genetic diversity, differentiation, and selection signatures of yak breeds/populations in Qinghai, China

The Qinghai Province of China is located in the northeast region of the Qinghai-Tibetan Plateau (QTP) and carries abundant yak genetic resources. Previous investigations of archaeological records, mitochondrial DNA, and Y chromosomal markers have suggested that Qinghai was the major center of yak domestication. In the present study, we examined the genomic diversity, differentiation, and selection signatures of 113 Qinghai yak, including 42 newly sequenced Qinghai yak and 71 publicly available individuals, from nine yak breeds/populations (wild, Datong, Huanhu, Xueduo, Yushu, Qilian, Geermu, Tongde, and Huzhu white) using high-depth whole-genome resequencing data. We observed that most of Qinghai yak breeds/populations have abundant genomic diversity based on four genomic parameters (nucleotide diversity, inbreeding coefficients, linkage disequilibrium decay, and runs of homozygosity). Population genetic structure analysis showed that Qinghai yak have two lineages with two ancestral origins and that nine yak breeds/populations are clustered into three distinct groups of wild yak, Geermu yak, and seven other domestic yak breeds/populations. F _ST values showed moderate genetic differentiation between wild yak, Geermu yak, and the other Qinghai yak breeds/populations. Positive selection signals were detected in candidate genes associated with disease resistance (CDK2AP2, PLEC, and CYB5B), heat stress (NFAT5, HSF1, and SLC25A48), pigmentation (MCAM, RNF26, and BOP1), vision (C1QTNF5, MFRP, and TAX1BP3), milk quality (OPLAH and GRINA), neurodevelopment (SUSD4, INSYN1, and PPP1CA), and meat quality (ZRANB1), using the integrated PI, composite likelihood ratio (CLR), and F _ST methods. These findings offer new insights into the genetic mechanisms underlying target traits in yak and provide important information for understanding the genomic characteristics of yak breeds/populations in Qinghai.

Analysis of genetic diversity and selection characteristics using the whole genome sequencing data of five buffaloes, including Xilin buffalo, in Guangxi, China

Buffalo is an economically important livestock that renders useful services to manhood in terms of meat, milk, leather, and draught. The Xilin buffalo is among the native buffalo breeds of China. In the present study, the genetic architecture and selection signature signals of Xilin buffalo have been explored. Correlation analysis of the population structure of Xilin buffalo was conducted by constructing NJ tree, PCA, ADMIXTURE and other methods. A total of twenty-five (n = 25) Xilin buffalo whole genome data and data of forty-six (n = 46) buffaloes published data were used. The population structure analysis showed that the Xilin buffalo belong to the Middle-Lower Yangtze. The genome diversity of Xilin buffalo was relatively high. The CLR, π ratio, F _ST, and XP-EHH were used to detect the candidate genes characteristics of positive selection in Xilin buffalo. Among the identified genes, most of the enriched signal pathways were related to the nervous system and metabolism. The mainly reported genes were related to the nervous system (GRM5, GRIK2, GRIA4), reproductive genes (CSNK1G2, KCNIP4), and lactation (TP63). The results of this study are of great significance for understanding the molecular basis of phenotypic variation of related traits of Xilin buffalo. We provide a comprehensive overview of sequence variations in Xilin buffalo genomes. Selection signatures were detected in genomic regions that are possibly related to economically important traits in Xilin buffalo and help in future breeding and conservation programs of this important livestock genetic resource.

Genomic diversity and relationship analyses of endangered German Black Pied cattle (DSN) to 68 other taurine breeds based on whole-genome sequencing

German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind, DSN) are an endangered dual-purpose cattle breed originating from the North Sea region. The population comprises about 2,500 cattle and is considered one of the ancestral populations of the modern Holstein breed. The current study aimed at defining the breeds closest related to DSN cattle, characterizing their genomic diversity and inbreeding. In addition, the detection of selection signatures between DSN and Holstein was a goal. Relationship analyses using fixation index (F_ST), phylogenetic, and admixture analyses were performed between DSN and 68 other breeds from the 1000 Bull Genomes Project. Nucleotide diversity, observed heterozygosity, and expected heterozygosity were calculated as metrics for genomic diversity. Inbreeding was measured as excess of homozygosity (F_Hom) and genomic inbreeding (F_RoH) through runs of homozygosity (RoHs). Region-wide F_ST and cross-population-extended haplotype homozygosity (XP-EHH) between DSN and Holstein were used to detect selection signatures between the two breeds, and RoH islands were used to detect selection signatures within DSN and Holstein. DSN showed a close genetic relationship with breeds from the Netherlands, Belgium, Northern Germany, and Scandinavia, such as Dutch Friesian Red, Dutch Improved Red, Belgian Red White Campine, Red White Dual Purpose, Modern Angler, Modern Danish Red, and Holstein. The nucleotide diversity in DSN (0.151%) was higher than in Holstein (0.147%) and other breeds, e.g., Norwegian Red (0.149%), Red White Dual Purpose (0.149%), Swedish Red (0.149%), Hereford (0.145%), Angus (0.143%), and Jersey (0.136%). The F_Hom and F_RoH values in DSN were among the lowest. Regions with high F_ST between DSN and Holstein, significant XP-EHH regions, and RoH islands detected in both breeds harbor candidate genes that were previously reported for milk, meat, fertility, production, and health traits, including one QTL detected in DSN for endoparasite infection resistance. The selection signatures between DSN and Holstein provide evidence of regions responsible for the dual-purpose properties of DSN and the milk type of Holstein. Despite the small population size, DSN has a high level of diversity and low inbreeding. F_ST supports its relatedness to breeds from the same geographic origin and provides information on potential gene pools that could be used to maintain diversity in DSN.

Genome-wide search reveals the uniqueness of DNA regions associated with coat color and innate immunity in Hokkaido Native Horse

Hokkaido Native Horse (HKD) is a horse breed native to Hokkaido in Japan known for the traits such as coat color with no white spots and adaptability to the local cold climate. To examine whether those traits of HKD are conferred at the DNA level, we attempted to identify fixed DNA regions in HKD individuals, that is, the selection signatures of HKD. A comparison of genome-wide single nucleotide polymorphism genotypes in 58 HKD individuals by principal component analysis, and cluster analysis between breeds, including HKD, and within the HKD individuals indicated the genetic independence of HKD as a breed. Tajima's D analysis and runs of homozygosity analysis identified 23 selection signatures unique to HKD (P < 0.05), and following database search found 20 traits that were associated with those selection signatures; among these traits, coat color traits, face and body markings, showed the highest important value (0.50 and 0.46). Enrichment analysis of genes in the selection signatures identified six gene ontology terms (P < 0.05), and a term related to innate immunity (regulation of defense response; GO:0031347) showed the highest positive fold enrichment value (7.13). These results provide the first scientific evidence of a genetic basis for the traits of HKD.

Population Structure and Selection Signatures Underlying Domestication Inferred from Genome-Wide Copy Number Variations in Chinese Indigenous Pigs

Single nucleotide polymorphism was widely used to perform genetic and evolution research in pigs. However, little is known about the effect of copy number variation (CNV) on characteristics in pigs. This study performed a genome-wide comparison of CNVs between Wannan black pigs (WBP) and Asian wild boars (AWB), using whole genome resequencing data. By using Manta, we detected in total 28,720 CNVs that covered approximately 1.98% of the pig genome length. We identified 288 selected CNVs (top 1%) by performing Fst statistics. Functional enrichment analyses for genes located in selected CNVs were found to be muscle related (NDN, TMOD4, SFRP1, and SMYD3), reproduction related (GJA1, CYP26B1, WNT5A, SRD5A2, PTPN11, SPEF2, and CCNB1), residual feed intake (RFI) related (MAP3K5), and ear size related (WIF1). This study provides essential information on selected CNVs in Wannan black pigs for further research on the genetic basis of the complex phenotypic and provides essential information for direction in the protection and utilization of Wannan black pig.

Genome-wide analysis reveals selection signatures for body size and drought adaptation in Liangzhou donkey

Liangzhou donkey is a domestic animal breed distributed on the edge of the Tengger Desert in Gansu Province of China. It has small body size and strong adaptability to dry environments. Here, we sequenced 10 Liangzhou donkey genomes and compared them to the 55 genomes of 8 representative donkey breeds worldwide. The population structure analysis revealed that Liangzhou donkey harboured the ancestry with the Asian domestic donkeys (0.863) and European domestic donkeys (0.137). Three methods (nucleotide diversity, linkage disequilibrium decay and runs of homozygosity) implied the genetic diversity in Liangzhou donkey. In addition, we analyzed the genetic basis of the small body size and drought adaptation of Liangzhou donkey by using F_st, θπ-ratio, XP-EHH, CLR and θπ methods. We found that the NCAPG-LCORL on chromosome 3 may be a candidate region for small body size trait of Liangzhou donkey. The CYP4A11 gene located on chromosome 5 showed strong sign of selection sweep. CYP4A11 can convert arachidonic acid into 19(S)-HETE, which can promote water reabsorption in renal tubule and enhance the ability of Liangzhou donkey to adapt to dry environment. These results contribute to a better understanding of the underlying population structure of Liangzhou donkeys and provides a valuable resource for future research on donkey breeding in response to climate change.

Integrative QTL mapping and selection signatures in Groningen White Headed cattle inferred from whole-genome sequences

Here, we aimed to identify and characterize genomic regions that differ between Groningen White Headed (GWH) breed and other cattle, and in particular to identify candidate genes associated with coat color and/or eye-protective phenotypes. Firstly, whole genome sequences of 170 animals from eight breeds were used to evaluate the genetic structure of the GWH in relation to other cattle breeds by carrying out principal components and model-based clustering analyses. Secondly, the candidate genomic regions were identified by integrating the findings from: a) a genome-wide association study using GWH, other white headed breeds (Hereford and Simmental), and breeds with a non-white headed phenotype (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel, Dutch Belted, and Holstein Friesian); b) scans for specific signatures of selection in GWH cattle by comparison with four other Dutch traditional breeds (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel and Dutch Belted) and the commercial Holstein Friesian; and c) detection of candidate genes identified via these approaches. The alignment of the filtered reads to the reference genome (ARS-UCD1.2) resulted in a mean depth of coverage of 8.7X. After variant calling, the lowest number of breed-specific variants was detected in Holstein Friesian (148,213), and the largest in Deep Red (558,909). By integrating the results, we identified five genomic regions under selection on BTA4 (70.2-71.3 Mb), BTA5 (10.0-19.7 Mb), BTA20 (10.0-19.9 and 20.0-22.7 Mb), and BTA25 (0.5-9.2 Mb). These regions contain positional and functional candidate genes associated with retinal degeneration (e.g., CWC27 and CLUAP1), ultraviolet protection (e.g., ERCC8), and pigmentation (e.g. PDE4D) which are probably associated with the GWH specific pigmentation and/or eye-protective phenotypes, e.g. Ambilateral Circumocular Pigmentation (ACOP). Our results will assist in characterizing the molecular basis of GWH phenotypes and the biological implications of its adaptation.

Assessing the Role of Ancestral Fragments and Selection Signatures by Whole-Genome Scanning in Dehong Humped Cattle at the China-Myanmar Border

Dehong humped cattle are precious livestock resources of Yunnan Province, China; they have typical zebu traits. Here, we investigated their genetic characteristics using whole-genome resequencing data of Dehong humped animals (n = 18). When comparing our data with the publicly-available data, we found that Dehong humped cattle have high nucleotide diversity. Based on clustering models in a population structure analysis, Dehong humped cattle had a mutual genome ancestor with Chinese and Indian indicine cattle. While using the RFMix method, it is speculated that the body sizes of Dehong humped cattle were influenced by the Chinese indicine segments and that the immune systems of Dehong humped cattle were affected by additional ancestral segments (Indian indicine). Furthermore, we explored the position selection regions harboring genes in the Dehong humped cattle, which were related to heat tolerance (FILIP1L, ABHD6) and immune responses (GZMM, PRKCZ, STOML2, LRBA, PIK3CD). Notably, missense mutations were detected in the candidate gene ABHD6 (c.870C>A p.Asp290Glu; c.987C>A p.Ser329Arg). The missense mutations may have implications for Dehong humped cattle adaptation to hot environments. This study provides valuable genomic resource data at the genome-wide level and paves the way for future genetic breeding work in the Dehong humped cattle.

Genetic and Genomic Characterization of a New Beef Cattle Composite Breed (Purunã) Developed for Production in Pasture-Based Systems

Purunã is a composite beef cattle breed, developed in Southern Brazil by crossing the Angus, Charolais, Canchim, and Caracu breeds. The goal of this study was to perform the first genetic characterization of the Purunã breed, based on both pedigree and genomic information. For this, 100 randomly selected animals were genotyped, and 11,205 animals born from 1997 to 2019 had pedigree information. The genetic analyses performed were principal component analysis, admixture, phylogenic tree, pedigree and genomic inbreeding, linkage disequilibrium (LD), effective population size (Ne), consistency of the gametic phase, runs of homozygosity (ROH), heterozygosity-enriched regions (HERs), and functional analyses of the ROH and HER regions identified. Our findings indicate that Purunã is more genetically related to the Charolais, Canchim, and Angus breeds than Caracu or Nellore. The levels of inbreeding were shown to be small based on all the metrics evaluated and ranged from −0.009 to 0.029. A low (−0.12–0.31) correlation of the pedigree-based inbreeding compared to all the genomic inbreeding coefficients evaluated was observed. The LD average was 0.031 (±0.0517), and the consistency of the gametic phase was shown to be low for all the breed pairs, ranging from 0.42 to 0.27 to the distance of 20 Mb. The Ne values based on pedigree and genomic information were 158 and 115, respectively. A total of 1,839 ROHs were found, and the majority of them are of small length (&lt;4 Mb). An important homozygous region was identified on BTA5 with pathways related to behavioral traits (sensory perception, detection of stimulus, and others), as well as candidate genes related to heat tolerance (MY O 1A), feed conversion rate (RDH5), and reproduction (AMDHD1). A total of 1,799 HERs were identified in the Purunã breed with 92.3% of them classified within the 0.5–1 Mb length group, and 19 HER islands were identified in the autosomal genome. These HER islands harbor genes involved in growth pathways, carcass weight (SDCBP), meat and carcass quality (MT2A), and marbling deposition (CISH). Despite the genetic relationship between Purunã and the founder breeds, a multi-breed genomic evaluation is likely not feasible due to their population structure and low consistency of the gametic phase among them.

Genomic Diversity and Selection Signatures for Weining Cattle on the Border of Yunnan-Guizhou

Weining cattle is a Chinese indigenous breed influenced by complex breeding and geographical background. The multi-ethnic breeding culture makes Weining cattle require more attention as livestock resources for its genetic diversity. Here, we used 10 Weining cattle (five newly sequenced and five downloaded) and downloaded another 48 genome data to understand the aspects of Weining cattle: genetic diversity, population structure, and cold-adapted performance. In the current study, a high level of genetic diversity was found in Weining cattle, and its breed comprised two potential ancestries, which were Bos taurus and Bos indicus. The positive selective sweep analysis in Weining cattle was analyzed using composite likelihood ratio (CLR) and nucleotide diversity (θπ), resulting in 203 overlapped genes. In addition, we studied the cold adaptation of Weining cattle by comparing with other Chinese cattle (Wannan and Wenshan cattle) by three methods (FST, θπ-ratio, and XP-EHH). Of the top 1% gene list, UBE3D and ZNF668 were analyzed, and these genes may be associated with fat metabolism and blood pressure regulation in cold adaptation. Our findings have provided invaluable information for the development and conservation of cattle genetic resources, especially in southwest China.

Assessing genomic diversity and signatures of selection in Pinan cattle using whole-genome sequencing data

BACKGROUND

Crossbreeding is an important way to improve production beef cattle performance. Pinan cattle is a new hybrid cattle obtained from crossing Piedmontese bulls with Nanyang cows. After more than 30 years of cross-breeding, Pinan cattle show a variety of excellent characteristics, including fast growth, early onset of puberty, and good meat quality. In this study, we analyzed the genetic diversity, population structure, and genomic region under the selection of Pinan cattle based on whole-genome sequencing data of 30 Pinan cattle and 169 published cattle genomic data worldwide.  RESULTS: Estimating ancestry composition analysis showed that the composition proportions for our Pinan cattle were mainly Piedmontese and a small amount of Nanyang cattle. The analyses of nucleotide diversity and linkage disequilibrium decay indicated that the genomic diversity of Pinan cattle was higher than that of European cattle and lower than that of Chinese indigenous cattle. De-correlated composite of multiple selection signals, which combines four different statistics including θπ, CLR, F_ST, and XP-EHH, was computed to detect the signatures of selection in the Pinan cattle genome. A total of 83 genes were identified, affecting many economically important traits. Functional annotation revealed that these selected genes were related to immune (BOLA-DQA2, BOLA-DQB, LSM14A, SEC13, and NAALADL2), growth traits (CYP4A11, RPL26, and MYH10), embryo development (REV3L, NT5E, CDX2, KDM6B, and ADAMTS9), hornless traits (C1H21orf62), and climate adaptation (ANTXR2).

CONCLUSION

In this paper, we elucidated the genomic characteristics, ancestry composition, and selective signals related to important economic traits in Pinan cattle. These results will provide the basis for further genetic improvement of Pinan cattle and reference for other hybrid cattle related studies.

Whole-Genome Resequencing of Xiangxi Cattle Identifies Genomic Diversity and Selection Signatures

Understanding the genetic diversity in Xiangxi cattle may facilitate our efforts toward further breeding programs. Here we compared 23 Xiangxi cattle with 78 published genomes of 6 worldwide representative breeds to characterize the genomic variations of Xiangxi cattle. Based on clustering models in population structure analysis, we displayed that Xiangxi cattle had a mutual genome ancestor with Chinese indicine, Indian indicine, and East Asian taurine. Population genetic diversity was analyzed by four methods (nucleotide diversity, inbreeding coefficient, linkage disequilibrium decay and runs of homozygosity), and we found that Xiangxi cattle had higher genomic diversity and weaker artificial selection than commercial breed cattle. Using four testing methods (θπ, CLR, F_ST, and XP-EHH), we explored positive selection regions harboring genes in Xiangxi cattle, which were related to reproduction, growth, meat quality, heat tolerance, and immune response. Our findings revealed the extent of sequence variation in Xiangxi cattle at the genome-wide level. All of our fruitful results can bring about a valuable genomic resource for genetic studies and breed protection in the future.

The Investigation of Protein Profile and Meat Quality in Bovine Longissimus thoracic Frozen under Different Temperatures by Data-Independent Acquisition (DIA) Strategy

The influence of freezing on the protein profile and quality traits in bovine Longissimus thoracic (LT) muscle was investigated by the data-independent acquisition (DIA) technique. Compared to fresh meat, a total of 262 proteins were identified as differential abundance proteins (DAPs) in four frozen groups (−12 °C, −18 °C, −38 °C, and −80 °C). According to the bioinformatics analysis, most of the DAPs in the significant Go terms and the KEGG pathway were structure proteins and enzymes. Proteome changes in the frozen bovine muscle at −12 °C and −18 °C were more significant than those at −38 °C and −80 °C. The result was consistent with the deterioration trend of the meat quality. The correlation analysis revealed that 17 proteins were correlated closely with the color, shear force, thawing loss, and cooking loss of the frozen meat, which could be used as putative biomarkers for frozen meat quality. MYO18A and ME3 are newly discovered proteins that are associated with frozen beef quality. In addition, CTTN and SERPINB6 were identified in frozen groups, which exhibited a significant inverse correlation with thawing loss (p < 0.01). These findings reveal the quality changes induced by freezing at the protein molecular level and provide new insights into the control of quality deterioration.

Whole Genome Sequencing Provides New Insights Into the Genetic Diversity and Coat Color of Asiatic Wild Ass and Its Hybrids

The diversity of livestock coat color results from human positive selection and is an indispensable part of breed registration. As an important biodiversity resource, Asiatic wild ass has many special characteristics, including the most visualized feature, its yellowish-brown coat color, and excellent adaptation. To explore the genetic mechanisms of phenotypic characteristics in Asiatic wild ass and its hybrids, we resequenced the whole genome of one Mongolian Kulan (a subspecies of Asiatic wild ass) and 29 Kulan hybrids (Mongolian Kulan ♂×Xinjiang♀), and the ancestor composition indicated the true lineage of the hybrids. XP-EHH (Cross Population Extended Haplotype Homozygosity), θπ-ratio (Nucleotide Diversity Ratio), CLR (Composite Likelihood Ratio) and θπ (Nucleotide Diversity) methods were used to detect the candidate regions of positive selection in Asiatic wild ass and its hybrids. Several immune genes (DEFA1, DEFA5, DEFA7, GIMAP4, GIMAP1, IGLC1, IGLL5, GZMB and HLA) were observed by the CLR and θπ methods. XP-EHH and θπ-ratio revealed that these genes are potentially responsible for coat color (KITLG) and meat quality traits (PDE1B and MYLK2). Furthermore, the heatmap was able to show the clear difference in the haplotype of the KITLG gene between the Kulan hybrids and Asiatic wild ass group and the Guanzhong black donkey group, which is a powerful demonstration of the key role of KITLG in donkey color. Therefore, our study may provide new insights into the genetic basis of coat color, meat quality traits and immunity of Asiatic wild ass and its hybrids.

Whole-Genome Sequencing of Endangered Dengchuan Cattle Reveals Its Genomic Diversity and Selection Signatures

Dengchuan cattle are the only dairy yellow cattle and endangered cattle among Yunnan native cattle breeds. However, its genetic background remains unclear. Here, we performed whole-genome sequencing of ten Dengchuan cattle. Integrating our data with the publicly available data, Dengchuan cattle were observed to be highly interbred than other cattle in the dataset. Furthermore, the positive selective signals were mainly manifested in candidate genes and pathways related to milk production, disease resistance, growth and development, and heat tolerance. Notably, five genes (KRT39, PGR, KRT40, ESR2, and PRKACB) were significantly enriched in the estrogen signaling pathway. Moreover, the missense mutation in the PGR gene (c.190T > C, p.Ser64Pro) showed a homozygous mutation pattern with higher frequency (83.3%) in Dengchuan cattle. In addition, a large number of strong candidate regions matched genes and QTLs related to milk yield and composition. Our research provides a theoretical basis for analyzing the genetic mechanism underlying Dengchuan cattle with excellent lactation and adaptability, crude feed tolerance, good immune performance, and small body size and also laid a foundation for genetic breeding research of Dengchuan cattle in the future.

Whole-Genome Analyses Reveal Genomic Characteristics and Selection Signatures of Lincang Humped Cattle at the China-Myanmar Border

The location on the Yunnan border with Myanmar and its unique cultural landscape has shaped Lincang humped cattle over time. In the current study, we investigated the genetic characteristics of 22 Lincang humped cattle using whole-genome resequencing data. We found that Lincang humped cattle derived from both Indian indicine and Chinese indicine cattle depicted higher levels of genomic diversity. Based on genome-wide scans, candidate genomic regions were identified that were potentially involved in local thermal and humid environmental adaptions, including genes associated with the body size (TCF12, SENP2, KIF1C, and PFN1), immunity (LIPH, IRAK3, GZMM, and ELANE), and heat tolerance (MED16, DNAJC8, HSPA4, FILIP1L, HELB, BCL2L1, and TPX2). Missense mutations were detected in candidate genes IRAK3, HSPA4, and HELB. Interestingly, eight missense mutations observed in the HELB gene were specific to the indicine cattle pedigree. These mutations may reveal differences between indicine and taurine cattle adapted to variable climatic conditions. Our research provides new insights into the genetic characteristics of Lincang humped cattle representing Lincang and Pu'er areas as an important channel for the migration of Indian indicine from domestication centers toward southwestern China.

Characterization of runs of homozygosity, heterozygosity-enriched regions, and population structure in cattle populations selected for different breeding goals

Background

A decline in the level of genetic diversity in livestock can result in reduced response to selection, greater incidence of genetic defects, and inbreeding depression. In this context, various metrics have been proposed to assess the level of genetic diversity in selected populations. Therefore, the main goals of this study were to: 1) investigate the population structure of 16 cattle populations from 15 different pure breeds or composite populations, which have been selected for different breeds goals; and, 2) identify and compare runs of homozygosity (ROH) and heterozygosity-enriched regions (HER) based on different single nucleotide polymorphism (SNP) panels and whole-genome sequence data (WGS), followed by functional genomic analyses.

Results

A total of 24,187 ROH were found across all cattle populations, with 55% classified in the 2-4 Mb size group. Fourteen homozygosity islands were found in five populations, where four ROH islands located on BTA1, BTA5, BTA16, and BTA19 overlapped between the Brahman (BRM) and Gyr (GIR) breeds. A functional analysis of the genes found in these islands revealed candidate genes known to play a role in the melanogenesis, prolactin signaling, and calcium signaling pathways. The correlations between inbreeding metrics ranged from 0.02 to 0.95, where the methods based on homozygous genotypes (F_HOM), uniting of gametes (F_UNI), and genotype additive variance (F_GRM) showed strong correlations among them. All methods yielded low to moderate correlations with the inbreeding coefficients based on runs of homozygosity (F_ROH). For the HER, 3576 runs and 26 islands, distributed across all autosomal chromosomes, were found in regions containing genes mainly related to the immune system, indicating potential balancing selection. Although the analyses with WGS did not enable detection of the same island patterns, it unraveled novel regions not captured when using SNP panel data.

Conclusions

The cattle populations that showed the largest amount of ROH and HER were Senepol (SEN) and Montana (MON), respectively. Overlapping ROH islands were identified between GIR and BRM breeds, indicating a possible historical connection between the populations. The distribution and pattern of ROH and HER are population specific, indicating that different breeds have experienced divergent selection processes or different genetic processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08384-0.

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.

Abstract

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, θπ, F_ST, θπ 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.