Kongshan cattle is an indigenous breed from Sichuan Province, China, characterized by their excellent meat quality, high fertility, strong disease resistance, and remarkable environmental adaptability. However, their genomic diversity has not been systematically studied. In this work, we performed whole-genome sequencing of 30 Kongshan cattle from a breeding farm and integrated these data with 113 representative commercial and indigenous cattle breeds worldwide to investigate their population structure and genetic diversity. We further analyzed the ancestral contributions to the development of the breed. The population structure revealed that Kongshan cattle possess four types of ancestral components: East Asian indicine (0.5974), East Asian taurine (0.3464), European taurine (0.0483), and Indian indicine (0.0079). The population also exhibits high nucleotide diversity, second only to pure East Asian indicine cattle. We inferred the ancestry of each variable site in the genome and, in combination with integrated haplotype score analysis, identified candidate genes related to meat quality (ME1, ENPP2, GPD2, PDZRN4, and TMTC2), immunity (MCM6, MAP3K6, PIP4K2A, CDC6, CDC25B, PTAFR, ZC3H10, and NEK6), and environmental adaptability (KCNJ15, BECN1, AOC2, DUSP5, and ST3GAL4). These findings provide valuable insights into the evolutionary history and ancestral origins of Kongshan cattle and contribute to the broader understanding, conservation, and sustainable utilization of indigenous Chinese cattle genetic resources.

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.

Ear size is a classical model for hot climate adaptation following the evolution, but the genetic basis of the traits associated with ear size remains to be elucidated. Here, we performed a genome-wide association study on 158 cattle to explain the genetic mechanism of ear size. One region on BTA6 between 36.79 and 38.80 Mb included 50 suggestive SNPs and 4 significant SNPs that were significantly associated with ear size. The most significant locus ( P = 1.30 × 10 −8 ) was a missense mutation (T250I) on the seventh exon of integrin-binding sialoprotein ( IBSP ), which had an allele substitution effect of 23.46 cm 2 for ear size. Furthermore, this mutation will cause changes in the three-dimensional structure of the protein. To further identify genes underlying this typical feature, we performed a genome scan among nine cattle breeds with different ear sizes by using SweeD. Results suggested that IBSP was under positive selection among four breeds with relatively large ear sizes. The expression levels of IBSP in ear tissues of large- and small-ear cattle were significantly different. A haplotype diversity survey of this missense mutation in worldwide cattle breeds strongly implied that the origin of this missense mutation event was Bos taurus . These findings have important theoretical importance for the exploration of major genes associated with ear size and provide important molecular markers for the identification of cattle germplasm resources.

Heat stress affects the animal production and causes serious economic losses to the husbandry. Tectonin beta-propeller repeat containing 2 (TECPR2) gene plays an important role in autophagy which may affect the temperature sensation in animals. A missense mutation (XM_024981840.1:c.3989 G > A p.Arg1330His) of the transcripts X4 in the bovine TECPR2 gene was identified. In this study, the c.3989 G > A variant in TECPR2 gene was genotyped in a total of 25 cattle breeds (520 individuals). Our results indicated that the frequency of A allele showed a decreasing pattern from southern cattle to northern cattle, while the frequency of G allele showed the opposite pattern, which was consistent with the climate distribution of China. Compared with the GG genotype, southern cattle carried more the AA and AG genotypes. Furthermore, the association results carried out that the frequencies of genotypes (GG, AG, AA) and the value of climate parameters (mean annual temperature (T), relative humidity (RH) and temperature humidity index (THI) were significantly correlated (p < 0.01). Hence, we speculated that the c.3989 G > A variant of TECPR2 gene was associated with the heat tolerance trait in Chinese cattle and the locus may be considered as a molecular marker for Chinese cattle breeding.

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.

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 (F ST, θπ-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.

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 (HE) and the observed heterozygosity (HO), 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, θπ, FST, θπ 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.

Bovine leukocyte antigen, class II, DO beta (BoLA-DOB) is related to antigen presentation, which can triggered by multicul factors. And the condition of immune function determines how much cattle load to heat stress. To evaluate the relationship between heat-resistance and single nucleotide polymorphisms (SNPs) in BoLA-DOB gene, our study has taken further analysis in Chinese indigenous cattle for the first time. A missense single nucleotide polymorphism (rs464874590) was detected in BoLA-DOB gene. We directly sequenced rs464874590 (NM_001013600.1 g.7122762 A G) in BoLA-DOB gene of 522 individuals of 26 cattle breeds. The frequency of allele G gradually decreases from south to north with distinct climatic distribution characteristics. Further association analysis was carried out between different genotypes and environmental parameters, including annual mean temperature (T), relative humidity (RH), and temperature-humidity index (THI). The result showed that three genotypes were significantly correlated with T, H, and THI (P < 0.01), indicating that GG genotype was distributed in areas with hot and moist conditions. Therefore, our results suggested that the rs464874590 could be applied as a genetic marker to detect the heat-resistance of Chinese indigenous cattle.

The accumulation of secondary metabolites, such as anthocyanins, in cells plays an important role in colored plants. The synthesis and accumulation of anthocyanins are regulated by multiple genes, of which the R2R3-MYB transcription factor gene family plays an important role. Based on the genomic data in the Potato Genome Sequencing Consortium database (PGSC) and the transcriptome data in the SRA, this study used potato as a model plant to comprehensively analyze the plant anthocyanin accumulation process. The results indicated that the most critical step in the synthesis of potato anthocyanins was the formation of p-coumaroyl-CoA to enter the flavonoid biosynthetic pathway. The up-regulated expression of the CHS gene and the down-regulated expression of HCT significantly promoted this process. At the same time, the anthocyanins in the potato were gradually synthesized during the process from leaf transport to tubers. New transcripts of stAN1 and PAL were cloned and named stAN1-like and PAL-like, respectively, but the functions of these two new transcripts still need further study. In addition, the sequence characteristics of amino acids in the R2-MYB and R3-MYB domains of potato were preliminarily identified. The aims of this study are to identify the crucial major genes that affect anthocyanin biosynthesis through multi-omics joint analysis and to transform quantitative traits into quality traits, which provides a basis and reference for the regulation of plant anthocyanin biosynthesis. Simultaneously, this study provides the basis for improving the anthocyanin content in potato tubers and the cultivation of new potato varieties with high anthocyanin content.

Glycoside Hydrolase 3 (GH3) is a phytohormone-responsive family of proteins found in many plant species. These proteins contribute to the biological activity of indolacetic acid (IAA), jasmonic acid (JA), and salicylic acid (SA). They also affect plant growth and developmental processes as well as some types of stress. In this study, GH3 genes were identified in 48 plant species, including algae, mosses, ferns, gymnosperms, and angiosperms. No GH3 representative protein was found in algae, but we identified 4 genes in mosses, 19 in ferns, 7 in gymnosperms, and several in angiosperms. The results showed that GH3 proteins are mainly present in seed plants. Phylogenetic analysis of all GH3 proteins showed three separate clades. Group I was related to JA adenylation, group II was related to IAA adenylation, and group III was separated from group II, but its function was not clear. The structure of the GH3 proteins indicated highly conserved sequences in the plant kingdom. The analysis of JA adenylation in relation to gene expression of GH3 in potato (Solanum tuberosum) showed that StGH3.12 greatly responded to methyl jasmonate (MeJA) treatment. The expression levels of StGH3.1, StGH3.11, and StGH3.12 were higher in the potato flowers, and StGH3.11 expression was also higher in the stolon. Our research revealed the evolution of the GH3 family, which is useful for studying the precise function of GH3 proteins related to JA adenylation in S. tuberosum when the plants are developing and under biotic stress.