Yanbian cattle is inhabitant of North of China, exhibiting many phenotypic features, such as long, dense body hair, and abundant intramuscular fat, designed to combat the extreme cold climate adaption. In the current study, we studied the cold tolerance of nine Yanbian cattle by whole genome resequencing and compared with African tropical cattle, N’Dama, as a control group. Yanbian cattle was aligned to the Bos taurus reference genome (ARS-UCD1.2) yielding an average of 10.8 fold coverage. The positive selective sweep analysis for the cold adaption in Yanbian cattle were analyzed using composite likelihood ratio (CLR) and nucleotide diversity (θπ), resulting in 292 overlapped genes. The strongest selective signal was found on BTA16 with potential mutation in CORT gene, a regulatory gene of primary hormone in the hypothalamic-pituitary-adrenal (HPA) axis, is reported to be associated with the cold stress, representedfour missense mutations (c.269C > T, p.Lys90Ile; c.251A > G, p.Glu84Gly; c.112C > T, p.Pro38Ser; c.86G > A, p.Pro29His). Meanwhile another gene on BTA6, showed significantly higher selective sweep signals for a cold adapted trait for hair follicle and length development, FGF5 (fibroblast growth factor 5) with a missense mutation (c.191C > T, p.Ser64Phe). Moreover, cold adapted Yanbian cattle was statistically compared with the hot adapted N’Dama cattle, a taurine cattle reported to show superior heat tolerance than zebu cattle, making them better adapted to the hot regions of Africa. XP-CLR, Fst, and θπ ratio were used to compare both breeds, yielding 487, 924, and 346 genes respectively. Among the 12 overlapped genes, (CD36) (c.638A > G, p.Lys 213Arg) involved in fat digestion and absorption plays an important role in membrane transport of long-chain fatty acid and its expression could increase in cold exposure. Henceforth, our study provides a novel genetic insights into the cold climate adaptation of Yanbian cattle and identified three candidate genes (CORT, FGF5, and CD36), which can add to an understanding of the cold climate adaptation of Yanbian cattle.

The small heat shock proteins (HSPB) are expressed in response to heat stress, and the heat shock protein family B (small) member 7 (HSPB7) gene has been reported to play an important role in heat tolerance pathways. Only a missense mutation (NC_037329.1: g.136054902 C > G: p.Ala69Gly) was identified in the HSPB7 gene in indicine cattle, which might be a candidate mutation associated with the heat tolerance. Here, we explore the allele frequency of this mutation in 774 individuals belonging to 32 Chinese indigenous cattle breeds using polymerase chain reaction (PCR) and DNA sequencing methods. The distribution of alleles of NC_037329.1: g.136054902 C > G displays significant geographical difference across native Chinese cattle breeds that the allele C was dominant in northern cattle groups, while allele G was dominant in southern indicine cattle groups. Additionally, the association analysis indicated that the G allele was significantly associated with mean annual temperature (T), relative humidity (RH), and temperature humidity index (THI) (p < 0.01), suggesting that cattle carrying allele G were distributed in regions with higher T, RH, and THI. Our results demonstrate that the mutation of the HSPB7 gene in Chinese indicine cattle might be a candidate gene associated with the heat tolerance.

Chinese indigenous cattle breeds have abundant genetic resources, which are valuable for the molecular breeding of cattle around the world. Thus, identifying important candidate genes and their genetic markers is very important for cattle molecular breeding. A previous study found that a nonsense mutation (rs378652941, c.483C > A, p.Cys161X) in the bovine transmembrane protein 95 gene (TMEM95) seriously reduced the reproductive performance in bulls, but few studies have detected this mutation in Chinese indigenous cattle breeds. Since the mutation c.483C > A may serve as a potential genetic marker for selecting higher fertility bulls, in the present study, using tetra-primer amplification refractory mutation system PCR (T-ARMS-PCR), forced PCR restriction fragment length polymorphism (forced PCR-RFLP), and DNA sequencing methods, the mutation c.483C > A was detected in 765 individuals from 13 Chinese cattle breeds. However, the results showed that this mutation did not exist at this locus in our analyzed breeds. Interestingly, we identified a newly frameshift insertion/deletion (indel) mutation (NC_037346.1: g.27056998_27057000delCT) in the bovine TMEM95 gene in 11 cattle breeds, which changed the location of the termination codon and changed the 16 amino acids in the C-terminal to 21 amino acids. Combined with previous studies, our study provides evidence that in Chinese cattle breeds the mutation c.483C > A cannot be used as a genetic marker in molecular breeding.

The Xinjiang brown cattle, Red steppe cattle, and Yunling cattle are indigenous cultivated cattle breeds in Chinese frontier provinces, and they produce high-grade beef and milk products, however, their genetic diversity in many important genes related to excellent meat and milk production is still unknown. Our previous studies have found that several candidate genes (e.g., SREBP1c and PAX7) were associated with bovine economically important phenotypic traits, but none has been reported in the above-mentioned three cattle breeds. Since the InDel (insertion/deletion) marker becomes a useful tool applied in the animal molecular breeding, herein, we firstly found that the InDel variations of seven candidate genes in these cattle. Results showed that the genotypic and allelic distributions of these seven genes were remarkably different among these three cattle (p < 0.05 or p < 0.01). Furthermore, the InDel variations of SREBP1c and PAX7 genes were significantly associated with eight phenotypic traits in Xinjiang brown cattle (p < 0.05 or p < 0.01), respectively, suggesting that they can become the useful DNA markers.