Next-generation sequencing has yielded a vast amount of cattle genomic data for global characterization of population genetic diversity and identification of genomic regions under natural and artificial selection. However, efficient storage, querying, and visualization of such large datasets remain challenging. Here, we developed a comprehensive database, the Bovine Genome Variation Database (BGVD). It provides six main functionalities: gene search, variation search, genomic signature search, Genome Browser, alignment search tools, and the genome coordinate conversion tool. BGVD contains information on genomic variations comprising ~60.44 M SNPs, ~6.86 M indels, 76,634 CNV regions, and signatures of selective sweeps in 432 samples from modern cattle worldwide. Users can quickly retrieve distribution patterns of these variations for 54 cattle breeds through an interactive source of breed origin map, using a given gene symbol or genomic region for any of the three versions of the bovine reference genomes (ARS-UCD1.2, UMD3.1.1, and Btau 5.0.1). Signals of selection sweep are displayed as Manhattan plots and Genome Browser tracks. To further investigate and visualize the relationships between variants and signatures of selection, the Genome Browser integrates all variations, selection data, and resources, from NCBI, the UCSC Genome Browser, and Animal QTLdb. Collectively, all these features make the BGVD a useful archive for in-depth data mining and analyses of cattle biology and cattle breeding on a global scale. BGVD is publicly available at http://animal.nwsuaf.edu.cn/BosVar.

Yunling cattle are a composite beef cattle breed, combining Brahman (1/2), Murray Grey (1/4) and Local Yunnan cattle (1/4), that was developed in Yunnan, China in the 1980s. Understanding the genetic information of Yunling cattle is of great significance to the development of reasonable breeding strategies for this breed. In the present study, we assessed the current genetic status of Yunling cattle in Yunnan Province (China) by analyzing the variability of the whole mitochondrial genome of 129 individuals. Altogether, 129 sequences displayed 47 different haplotypes. The haplotype diversity and the average number of nucleotide differences were 0.964 and 128.074, respectively. Phylogenetic analyses classified Yunling cattle into seven haplogroups: T1, T2, T3, T4, T6, I1 and I2. Haplogroup I1 was found to be predominant (41.86%), followed by T3 (28.68%). Furthermore, we also identified a novel haplogroup, T6, and defined the sub-haplogroup I1a in Yunling cattle. According to the formation process of Yunling cattle (local Yunnan cattle as the maternal line), the high genetic diversities in the mitochondria of Yunling cattle could be due to the complex maternal origin of local Yunnan cattle. Further studies about local Yunnan breeds are necessary to determine the exact source of haplogroup T6 in Yunling cattle. Our results will be useful for the evaluation and effective management of Yunling cattle.

In the biological process of testicular spermatogenesis, the expression and interaction of many genes are regulated by microRNAs (miRNAs). However, comparisons of miRNA expression between descended testes (DTs) and undescended testes (UDTs) are rarely done in horses. In this study, we selected two UDTs (CKY2b and GU4b) from Chakouyi (CKY) and Guanzhong (GU) horses and eight DTs (GU1–3, CKY1, CKY3, CKY2a, GU4a, and GU5). Three groups were compared to evaluate expression patterns of testicular miRNA in stallion testes. Group 1 compared normal CKY horses and GU horses (CKY1 and CKY3 vs. GU1–3). Group 2 (CKY2a and GU4a (DTs) vs. CKY2b and GU4b (UDTs)) and group 3 (GU1–3, CKY1, CKY3 (DTs) vs. CKY2b and GU4b (UDTs)) compared the expression levels in unilateral retained testes to normal testes. The results show that 42 miRNAs (7 upregulated and 35 downregulated) had significantly different expression levels in both comparisons. The expression levels of eca-miR-545, eca-miR-9084, eca-miR-449a, eca-miR-9024, eca-miR-9121, eca-miR-8908e, eca-miR-136, eca-miR-329b, eca-miR-370, and eca-miR-181b were further confirmed by quantitative real-time PCR assay. The target genes of differentially expressed miRNAs in three comparisons were predicted, and the functions were annotated. The putative target genes of the 42 co-differentially expressed miRNAs were annotated to 15 functional terms, including metal ion binding, GTPase activator activity, zinc ion binding, intracellular, cytoplasm, and cancer pathways, and osteoclast differentiation. Our data indicate that the differentially expressed miRNAs in undescended testis suggests a potential role in male fertility and a relationship with cryptorchidism in horses. The discovery of miRNAs in stallion testes might contribute to a new direction in the search for biomarkers of stallion fertility.

The Du’an goat is one of the most important farm animals in the Guangxi Autonomous Region of China, but the genetic basis underlying its adaptive traits has still not been investigated. Firstly, in this study, the genomes of 15 Du’an goats from a breeding farm were sequenced (mean depth: 9.50X) to analyze the patterns of genetic variation. A comparable diversity (17.3 million single nucleotide polymorphisms and 2.1 million indels) was observed to be associated with a lower runs of homozygosity-based inbreeding coefficient and smaller effective population size in comparison with other breeds. From selective sweep and gene set enrichment analyses, we revealed selective signals related to adaptive traits, including immune resistance (serpin cluster, INFGR1, TLR2, and immune-related pathways), body size (HMGA2, LCOR, ESR1, and cancer-related pathways) and heat tolerance (MTOR, ABCG2, PDE10A, and purine metabolism pathway). Our findings uncovered the unique diversity at the genomic level and will provide the opportunities for improvement of productivity in the Du’an goat.

Litter size is a critical economic trait in the sheep industry. Like many other breeds, Duolang sheep typically produce one lamb per ewe per lambing. To date, genetic studies of this trait in that breed have largely relied on candidate gene approaches. To expand the genomic resources for this breed, we sequenced 297 genomes, generating approximately 10.52 trillion bases with an average coverage of 13.35X. High-quality alignments with a mapping rate exceeding 99% enabled the identification of 43,968,128 SNPs and 6,504,047 InDels. This dataset provides a valuable resource for identifying genetic variants associated with litter size through genome-wide association studies (GWAS) and lays the foundation for future genetic improvement efforts through genomic selection in Duolang sheep. Beyond trait mapping, the dataset also supports broader applications, including analyses of genetic diversity, phylogenetic relationships, population history, adaptive introgression, and breed-specific characteristics. Additionally, the moderate-coverage WGS data are suitable for structural variant (SV) detection and downstream analyses such as association mapping and the identification of SVs underlying phenotypic traits.

Body conformation traits are linked to the health, longevity, reproductivity, and production performance of cattle. These traits are also crucial for herd selection and developing new breeds. This study utilized pedigree information and phenotypic (1185 records) and genomic (The resequencing of 496 Xinjiang Brown cattle generated approximately 74.9 billion reads.) data of Xinjiang Brown cattle to estimate the genetic parameters, perform factor analysis, and conduct a genome-wide association study (GWAS) for these traits. Our results indicated that most traits exhibit moderate to high heritability. The principal factors, which explained 59.12% of the total variance, effectively represented body frame, muscularity, rump, feet and legs, and mammary system traits. Their heritability estimates range from 0.17 to 0.73, with genetic correlations ranging from −0.53 to 0.33. The GWAS identified 102 significant SNPs associated with 12 body conformation traits. A few of the SNPs were located near previously reported genes and quantitative trait loci (QTLs), while others were novel. The key candidate genes such as LCORL, NCAPG, and FAM184B were annotated within 500 Kb upstream and downstream of the significant SNPs. Therefore, factor analysis can be used to simplify multidimensional conformation traits into new variables, thus reducing the computational burden. The identified candidate genes from GWAS can be incorporated into the genomic selection of Xinjiang Brown cattle, enhancing the reliability of breeding programs.