Whole genome sequencing of simmental cattle for SNP and CNV discovery.

BACKGROUD: The single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) are two major genomic variants, which play crucial roles in evolutionary and phenotypic diversity. RESULTS: In this study, we performed a comprehensive analysis to explore the genetic variations (SNPs and CNVs) of high sperm motility (HSM) and poor sperm motility (PSM) Simmental bulls using the high-coverage (25×) short-read next generation sequencing and single-molecule long reads sequencing data. A total of ~ 15 million SNPs and 2,944 CNV regions (CNVRs) were detected in Simmental bulls, and a set of positive selected genes (PSGs) and CNVRs were found to be overlapped with quantitative trait loci (QTLs) involving immunity, muscle development, reproduction, etc. In addition, we detected two new variants in LEPR, which may be related to the artificial breeding to improve important economic traits. Moreover, a set of genes and pathways functionally related to male fertility were identified. Remarkably, a CNV on SPAG16 (chr2:101,427,468 - 101,429,883) was completely deleted in all poor sperm motility (PSM) bulls and half of the bulls in high sperm motility (HSM), which may play a crucial role in the bull-fertility. CONCLUSIONS: In conclusion, this study provides a valuable genetic variation resource for the cattle breeding and selection programs.

Genome-wide association study reveals that the IBSP locus affects ear size in 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 cm2 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.

Maternal genetic diversity, differentiation and phylogeny of three white yak breeds/populations in China.

White yak is a unique economic livestock animal on the Qilian Mountains of Qinghai-Tibet Plateau (QTP). Here, 369 mitochondrial D-loop sequences from three Chinese white yak breeds/populations (Tianzhu, Menyuan and Huzhu) were comprehensively analyzed to indicate their maternal genetic diversity, differentiation and phylogenetic relationship. Our results showed that the haplotype diversity (Hd) was found to be highest in Tianzhu white yak (Hd ± SD = 0.9501 ± 0.0058), while the lowest was recorded in Huzhu white yak (Hd ± SD = 0.7178 ± 0.0474). Totally, the haplotype and nucleotide diversities of white yak were 0.9407 ± 0.0069 and 0.0187 ± 0.0094, respectively, indicating an abundant maternal genetic diversity in white yak. Estimates of FST (Fixation Index) values showed a high genetic differentiation between Tianzhu and Menyuan populations (FST = 0.2928, p < 0.05) as well as that between Tianzhu and Huzhu populations (FST = 0.2721, p < 0.05), but a moderate genetic differentiation between Menyuan and Huzhu populations (FST = 0.1352, p > 0.05) was observed. Cluster analysis based on FST values among populations indicated that the genetic relationship between Menyuan and Huzhu white yak was closer, but they had a far genetic relationship with Tianzhu white yak. Maternal phylogenetic analysis revealed that white yak represented two maternal lineages (I and II), suggesting two maternal origins. In addition, it is notable that taurine mtDNA haplotypes were detected in Tianzhu and Huzhu white yak populations, indicating taurine genetic introgression to some extent. Our study would provide useful information for the conservation and utilization of white yak.

Goat CLSTN2 gene: tissue expression profile, genetic variation, and its associations with litter size.

Calsyntenin-2 (CLSTN2) is involved in cell proliferation, differentiation, cell death, tumorigenesis, and follicular expression. Although CLSTN2 has been identified as a potential candidate gene for sheep prolificacy, no studies have been done on its effect on goat prolificacy. The purpose of this study was to identify mRNA expression and genetic variation within goat CLSTN2, and its association with prolificacy. Herein, we uncovered significant differences in mRNA levels of the CLSTN2 gene in different tissues in female goats (p < 0.01), including ovary tissue. Nine putative indels were designed to investigate their correlation to litter size, but only one 16-bp deletion was discovered in female Shaanbei white cashmere goats (n = 902). We discovered that a 16-bp deletion within the CLSTN2 gene was significantly correlated with first-born litter size (p = 0.0001). As shown by the chi-squared test, the genotypic II of single-lambs and multi-lambs was dramatically higher than with genotype ID (p = 0.005). Our findings suggest that indel within the CLSTN2 gene is a candidate gene affecting prolificacy in goats and may be used for Marker Assisted Selection (MAS) in goats.

Muscle transcriptome analysis reveal candidate genes and pathways related to fat and lipid metabolism in Yunling cattle.

Yunling cattle (YL) is a recently developed beef breed harboring a quarter of Yunnan ancestral cattle genome, spanning over past 30 years. Compared with Diqing cattle (DQ), a Yunnan native cattle breed, YL presents various advantages, including rapid growth and exquisite meat quality. However, the molecular mechanisms underlying these phenotypic differences are not clearly understood. To further identify the candidate genes responsible for the quality of the meat in the muscle, longissimus dorsi (LD) muscle was used for RNA-Seq analysis. A total of 508 differentially expressed genes (DEGs) were identified in YL (adjusted p-value <0.01 and log2FoldChange >1), of which 243 were up-regulated and 265 were down-regulated. Functional association analysis showed that the identified DEGs mainly enriched the lipid and fat metabolism pathways. Moreover, it was also observed that several fat-related genes were differentially expressed in both cattle breeds, including three up-regulated genes (MOGAT1, ACSM3, PLPP2) and two down-regulated genes (ADIG, GPAT3). In addition, alternative splice analysis was also performed revealing an important 9-11 exon skipping variation of GPAM gene (crucial for beef marbling) in YL, which is three times higher than that in DQ, suggesting that this variation might have played the central role in the 'snow beef' effect in YL. We believe that our results will help in understanding the mechanism of muscle development and promote the further breeding programs in YL cattle.

A novel SNP of TECPR2 gene associated with heat tolerance in Chinese cattle.

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.

Rich maternal and paternal genetic diversity and divergent lineage composition in wild yak (Bos mutus).

Wild yak (Bos mutus) is a vulnerable bovine species on the Qinghai-Tibetan Plateau (QTP). So far, most studies on molecular genetic diversity of wild yak have focused on autosomal and mtDNA variations based on small number of samples. In this study, we analyzed 84 D-loop and 24 whole mitogenome sequences of wild yak to further comprehensively explore its maternal genetic diversity and lineage composition. Meanwhile, using six yak Y-specific polymorphic markers (i.e., SRY4, USP9Y, UTY19, AMELY3, OFD1Y10 and INRA189), we assessed the paternal genetic diversity and lineage composition based on eight wild yak. Our results showed that wild yak exhibited abundant maternal genetic diversity with haplotype diversities of 0.9621 ± 0.0078 and 0.9928 ± 0.0144 in the D-loop and whole mitogenome sequences, respectively. Maternal phylogenetic analysis of wild yak uncovered three defined lineages (mt-I, mt-II and mt-III). Similarly, profuse paternal genetic diversity was observed in wild yak with Y-haplotype diversity (Hd) at 0.8214 ± 0.1007. Two Y-haplogroups (Y1 and Y2) with four Y-haplotypes (yH1-yH4) were identified in paternal phylogenetic analysis, indicating wild yak to be of two paternal lineages. This study of genetic diversity and lineage composition of wild yak would provide useful information for the genetic resource conservation and utilization of this vulnerable wild species.

Rich maternal and paternal genetic diversity and divergent lineage composition in wild yak (Bos mutus).

Wild yak (Bos mutus) is a vulnerable bovine species on the Qinghai-Tibetan Plateau. So far, most studies on the molecular genetic diversity of wild yak have focused on autosomal and mtDNA variations based on the small number of samples. In this study, we analyzed 84 D-loop and 24 whole mitogenome sequences of wild yak to further comprehensively explore its maternal genetic diversity and lineage composition. Meanwhile, using six yak Y-specific polymorphic markers (i.e., SRY4, USP9Y, UTY19, AMELY3, OFD1Y10 and INRA189), we assessed the paternal genetic diversity and lineage composition based on eight wild yak. Our results showed that wild yak exhibited abundant maternal genetic diversity with haplotype diversities of 0.9621 ± 0.0078 and 0.9928 ± 0.0144 in the D-loop and whole mitogenome sequences, respectively. Maternal phylogenetic analysis of wild yak uncovered three defined lineages (mt-I, mt-II and mt-III). Similarly, profuse paternal genetic diversity was observed in wild yak with Y-haplotype diversity at 0.8214 ± 0.1007. Two Y-haplogroups (Y1 and Y2) and four Y-haplotypes (yH1-yH4) were identified in paternal phylogenetic analysis, indicating wild yak to be of two paternal lineages. The present study of genetic diversity and lineage composition of wild yak would provide useful information for the genetic resource conservation and utilization of this vulnerable wild species.

Distribution of the variant at AKIRIN2: c.*188G > A in Chinese cattle.

Marbling score (MS), is an economically important trait in cattle. Previous results showed that a SNP (c.*188G > A) of akirin 2 (AKIRIN2) gene was associated with MS in Japanese Black cattle and Korean cattle. However, the distribution of the genotypic frequency of the single nucleotide polymorphism (SNP) has not been explored in Chinese cattle. In this study, we used polymerase chain reaction (PCR) and DNA sequencing to detect the variation in 1296 individuals from 39 Chinese cattle breeds, one semi-wild bovine species (Dulong) and three introduced breeds (Angus, Holstein and Brahman). Our study found the frequency of the A allele at this locus roughly diminished from north to south in Chinese cattle, and we detected statistically significant differences between Angus and Brahman (p < 0.05), Dulong and another two breeds (Angus and Holstein; p < 0.01) using Chi-Square Independence Test. Our results reflected the variation of AKIRIN2: c.*188G > A in Chinese cattle, which would help us better understand Chinese cattle genetic resources and provide reference for further research.

DGAT1 K232A polymorphism is associated with milk production traits in Chinese cattle.

The production traits of cattle, especially milk trait, are of great significance to human life. A quantitative trait loci (QTL) associated with milk fat content was detected in the centromeric region of cattle chromosome 14. This QTL harbors a strong candidate gene called DGAT1 responsible for the milk quality. A non-conservative substitution of lysine by alanine (K232A) was found in DGAT1 gene producing a strong effect on milk composition and yield. The lysine (K allele) is associated with increased milk fat content, while the decreased milk fat content is linked to the alanine (A allele) amino acid. To estimate the frequencies of the DGAT1 K232A polymorphism in Chinese cattle breeds, PCR and DNA sequencing methods were used to investigate the polymorphism of DGAT1 K232A in a total of 682 individuals, including 655 Chinese cattle and 27 Holstein cattle. The results demonstrated that the frequency of K allele gradually elevated from the northern group to the southern group of native Chinese cattle, whereas the frequency of A allele showed a contrary pattern, displaying a significant geographical difference across native Chinese cattle breeds. Our results confirm that the southern cattle group has higher milk fat content than that of the northern group.

Whole genome analyses revealed genomic difference between European taurine and East Asian taurine.

European taurine and East Asian taurine are two main clades in Bos taurus, but their genomic differences are not clearly elucidated. Here, we sequenced 16 Mongolian cattle genomes and compared them to the 92 genomes of 10 representative breeds worldwide. We found the highest LD level in Mishima cattle and the fastest LD decay in European taurine. Phylogenetic analysis revealed that Mongolian, Hanwoo and Mishima cattle were clustered into East Asian taurine. From selective sweep, gene annotation, functional enrichment and differential expression analysis, we identified selective signals including genes and/or pathways related to rapid growth and large body size in European taurine, and superior meat quality in East Asian taurine. Our findings will help us understand the evolutionary history and formation process of the breeds and provide theoretical materials regarding the genetic mechanism underlying breed characteristics and molecular breeding programmes of the taurine clades in the future.

Antiproliferative activity of 3,5-disubstituted tetrahydro-2H-1,3,5-thiadiazine thione (THTT) derivatives and evaluation as potential prodrug.

Four series of tetrahydro-2H-1,3,5-thiadiazine thione derivatives were screened for their in vitro antiproliferative activities against two human cancerous PC3 and HeLa cell lines. The cytotoxicity of all the compounds (series A-D) was also determined on mammalian mouse fibroblast 3T3 cells. Most of the compounds showed significant anticancer potential against both cancer cell lines within the range of IC50 = 6.4-29.9 and 2.4-23.8 M respectively when compared with standard doxorubicin (IC50 = 0.3 M). All compounds demonstrated a notable selectivity for Hela cells and found either non-toxic or relatively less toxic for 3T3 cell lines model. The structure-activity relationship indicated that antiproliferative activity mainly influenced by the nature and position of substituents at thidiazine nucleus. In general, the presence of aryl groups for example 3,4-(OMe) 2.Bzl and CH(Ph)Me at N-3 position resulted in a significant activity. Under enzymatic hydrolysis, complete conversion (100%) of ester derivative of thiadiazine thione (10a) into its acidic counterpart (7c) was achieved during 20 min which indicated that these types of THTT ester derivatives can be a possible lead for future investigations as prodrug anticancer probes.

Population structure, genetic diversity, and selective signature of Chaka sheep revealed by whole genome sequencing.

BACKGROUND: Chaka sheep, named after Chaka Salt Lake, are adapted to a harsh, highly saline environment. They are known for their high-grade meat quality and are a valuable genetic resource in China. Furthermore, the Chaka sheep breed has been designated a geographical symbol of agricultural products by the Chinese Ministry of Agriculture. RESULTS: The genomes of 10 Chaka sheep were sequenced using next-generation sequencing, and compared to that of additional Chinese sheep breeds (Mongolian: Bayinbuluke and Tan; Tibetan: Oula sheep) to explore its population structure, genetic diversity and positive selection signatures. Principle component analysis and a neighbor-joining tree indicated that Chaka sheep significantly diverged from Bayinbuluke, Tan, and Oula sheep. Moreover, they were found to have descended from unique ancestors (K = 2 and K = 3) according to the structure analysis. The Chaka sheep genome demonstrated comparable genetic diversity from the other three breeds, as indicated by observed heterozygosity (Ho), expected heterozygosity (He), runs of homozygosity (ROH), linkage disequilibrium (LD) decay. The enrichment analysis revealed that in contrast to Mongolian or Tibetan lineage groups, the genes annotated by specific missense mutations of Chaka sheep were enriched with muscle structure development (GO:0061061) factors including insulin-like growth factor 1 (IGF1), growth differentiation factor 3 (GDF3), histone deacetylase 9 (HDAC9), transforming growth factor beta receptor 2 (TGFBR2), and calpain 3 (CAPN3), among others. A genome-wide scan using Fst and XP-CLR revealed a list of muscle-related genes, including neurofibromin 1 (NF1) and myomesin 1 (MYOM1), under potential selection in Chaka sheep compared with other breeds. CONCLUSIONS: The comprehensive genome-wide characterization provided the fundamental footprints for breeding and management of the Chaka sheep and confirmed that they harbor unique genetic resources.

Selection signatures of Fuzhong Buffalo based on whole-genome sequences.

BACKGROUND: Fuzhong buffalo, a native breed of Guangxi Zhuang Autonomous Region, is traditionally used as a draft animal to provide farm power in the rice cultivation. In addition, the Fuzhong buffalo also prepared for the bullfighting festival organized by the locals. The detection of the selective signatures in its genome can help in elucidating the selection mechanisms in its stamina and muscle development of a draft animal. RESULTS: In this study, we analyzed 27 whole genomes of buffalo (including 15 Fuzhong buffalo genomes and 12 published buffalo genomes from Upper Yangtze region). The ZHp, ZFst, π-Ratio, and XP-EHH statistics were used to identify the candidate signatures of positive selection in Fuzhong buffalo. Our results detected a set of candidate genes involving in the pathways and GO terms associated with the response to exercise (e.g., ALDOA, STAT3, AKT2, EIF4E2, CACNA2D2, TCF4, CDH2), immunity (e.g., PTPN22, NKX2-3, PIK3R1, ITK, TMEM173), nervous system (e.g., PTPN21, ROBO1, HOMER1, MAGI2, SLC1A3, NRG3, SNAP47, CTNNA2, ADGRL3). In addition, we also identified several genes related to production and growth traits (e.g., PHLPP1, PRKN, MACF1, UCN3, RALGAPA1, PHKB, PKD1L). Our results depicted several pathways, GO terms, and candidate genes to be associated with response to exercise, immunity, nervous system, and growth traits. CONCLUSIONS: The selective sweep analysis of the Fuzhong buffalo demonstrated positive selection pressure on potential target genes involved in behavior, immunity, and growth traits, etc. Our findings provided a valuable resource for future research on buffalo breeding and an insight into the mechanisms of artificial selection.

Whole-genome resequencing reveals diversity, global and local ancestry proportions in Yunling cattle.

Yunling cattle, a three-breed cross consisting of 1/2 Brahman cattle, 1/4 Murray Grey cattle and 1/4 Yunnan Yellow cattle, has advantage of rapid growth, good meat quality, enhanced tolerance towards a hot and humid climate, tick resistance and crude feed forbearance. Here, we investigated the genetic diversity, global and local ancestry proportions by sequencing the genomes of 131 Yungling cattle and 31 Brahman cattle (for control). Furthermore, we used 35 published genomes of ancestry breeds (including Angus cattle [ancestry of Murray Grey], Wannan cattle [ancestry of Yunnan Yellow cattle: Wenshan and Dianzhong], Wenshan and Dianzhong cattle) to characterize the formation process of Yunling cattle. The highest nucleotide diversity was found in the Wannan cattle, followed by Wenshan, Dianzhong, Brahman, Yunling and Angus cattle. The results of LD decay in each breed was largely consistent with the results of nucleotide diversity, except for the faster decay in Angus cattle out of a smaller effective population size and a strong bottleneck during the breed formation. The population-structure analysis revealed that the cross-breed Yunling cattle harboured the ancestry with Angus (0.44), Brahman (0.48) and Wannan cattle (0.08). Subsequently, we used RFmix to infer local ancestry in Yunling cattle and then performed chi-squared test to obtain the segments whose proportions of a certain ancestry were excessive compared with the whole-genome level, leading to 391 Angus, 49 Brahman and 2,312 Wannan segments. Gene annotation and KEGG enrichment analysis revealed that the excessive Angus, Brahman and Wannan segments might contribute to the rapid growth, immune resistance and indigenous adaptation, respectively, in Yunling cattle. Our results help understand ancestry components and formation process in Yunling cattle and will provide an opportunity for selective breeding by molecular approaches in future.

Genome-wide association study identifies genomic loci associated with flight reaction in cattle.

Lower flight reaction is closely related to higher production in cattle, but the genetic basis of lower flight reaction is not clearly understood. Here, we sampled a total of 45 Brahman cattles and 166 Yunling cattles with flight distance (FD), and 73 Brahman cattles and 288 Yunling cattles with crush score (CS) and flight speed (FS), whereas there were 45 Brahman cattles and 161 Yunling cattles with all three traits. The FD, CS and FS in Brahman cattle were significantly lower than those in Yunling cattle. The flight reaction traits had negative correlation with conformational traits (e.g., body weight, withers height and body length). Based on SNPs derived from a subset of 162 whole genomes (25 Brahman genomes and 100 Yunling genomes with FD, 30 Brahman and 131 Yunling genomes with CS and FS), genome-wide association study with mixed linear model was performed to test potential associations between flight reaction traits and genomic variants. We identified five, two and two genomic loci suggestively associated with FD, CS and FS, respectively. Five out of five candidate genes for FD (LOC789753, LRP6, CTIF, SLC9A9 and ZEB1) were reported to be related to Alzheimer's disease representing cognitive impairment in human, which was consistent with the finding that cognitive-behavioural intervention decreased the FD of cows to human. In CS, a very strong association locus was assigned to CDH8, a cadherin involved in synaptic adhesion, axon outgrowth and guidance, whose deletion was associated with autism spectrum disorder. In FS, a very strong association locus was assigned to GABRG2, a gamma-aminobutyric acid (a major inhibitory neurotransmitter in brain) receptor, whose polymorphisms were associated with suicidal behaviour in schizophrenia patients. Our findings will provide targets for molecular-marker selection and genetic manipulation of cattle improvement to meeting the growing demand for lower flight reaction to human.

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.

Global genetic diversity, introgression, and evolutionary adaptation of indicine cattle revealed by whole genome sequencing.

Indicine cattle, also referred to as zebu (Bos taurus indicus), play a central role in pastoral communities across a wide range of agro-ecosystems, from extremely hot semiarid regions to hot humid tropical regions. However, their adaptive genetic changes following their dispersal into East Asia from the Indian subcontinent have remained poorly documented. Here, we characterize their global genetic diversity using high-quality whole-genome sequencing data from 354 indicine cattle of 57 breeds/populations, including major indicine phylogeographic groups worldwide. We reveal their probable migration into East Asia was along a coastal route rather than inland routes and we detected introgression from other bovine species. Genomic regions carrying morphology-, immune-, and heat-tolerance-related genes underwent divergent selection according to Asian agro-ecologies. We identify distinct sets of loci that contain promising candidate variants for adaptation to hot semi-arid and hot humid tropical ecosystems. Our results indicate that the rapid and successful adaptation of East Asian indicine cattle to hot humid environments was promoted by localized introgression from banteng and/or gaur. Our findings provide insights into the history and environmental adaptation of indicine cattle.

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.