A missense mutation (rs209302038) of KRT9 gene associated with heat stress in Chinese cattle.

Type I keratin 9 encoded by the KRT9 gene serves an important special function either in the mature palmar and plantar skin tissue. The changes in skin conditions and thickening of the outer layer of the skin may be affected by environmental variables. A missense mutation rs209302038 (NC_037346.1: g.41782870 G > A) was detected in KRT9, which changing the isoleucine into valine. This study aimed to identify the frequency of allele in this locus in Chinese indigenous cattle, and analyze the connection with heat stress. Our results indicated that the frequency of allele A gradually decreases from south to north, while the frequency of G allele showed the opposite pattern. Further analysis of the association of the different genotypes with three climate factors, which showed that the genotypes (GG, GA, AA) were significantly related to climatic conditions (p < 0.01). Therefore, we speculated that the mutation of the rs209302038 in Chinese indigenous cattle might be a genetic marker to detect heat stress.

TAS2R16 introgression from banteng into indigenous Chinese cattle.

Historical hybridization between southern indigenous Chinese cattle and banteng has been well-documented and has resulted in gene introgression. Bitter taste receptors were reported in indigenous cattle as a result of introgression from banteng. To determine the level of introgression of the taste 2 receptor member 16 (TAS2R16) gene from banteng into Chinese cattle, two missense mutations in the bovine TAS2R16 gene were examined. Here, we explored the prevalence of the two variants in 28 indigenous Chinese cattle and banteng breeds (comprising 750 individuals) to determine the influence of banteng introgressions on Chinese cattle based on PCR and DNA sequencing. In our study, the two mutant alleles had a higher frequency distribution in southern China with strong geographic distribution, especially in the south-central and southeast areas. In conclusion, this study examines the impact of introgression on the frequency distributions of mutations in variable regions and the subsequent adaptation of Chinese indigenous cattle to different environmental conditions.

Exploring genotype-phenotype relationships of the CRABP2 gene on growth traits in beef cattle.

Cellular retinoic acid binding protein 2 (CRABP2) is essential to myoblast differentiation. However, there was no report about the function of CRABP2 gene in cattle. This study explored the association of CRABP2 gene polymorphisms with growth traits in cattle breeds by several methods, such as DNA sequencing, PCR, PCR-RFLP and forced PCR-RFLP. Two sequence variants were determined. There were 621 individuals in six cattle breeds from China for the experiment, and three breeds were used to test validation of polymorphisms and extent of linkage disequilibrium (LD). The results showed that both SNPs (SNP1, g.2458 G > T, SNP2, g.3878 G > A) were in intron1. Two SNPs were in low linkage disequilibrium. Association analysis suggested that SNP1 had the significant difference on growth traits with body height, height at hip cross and body slanting length (P < .05), while SNP2 showed a significant difference in growth traits with body height, height at hip cross and body slanting length(P < .05). The results of this investigation displayed that the CRABP2 gene is an available candidate gene and may be used for breed selection and conservation.

Integrative analysis of APOL3 gene CNV for adult cattle stature.

Copy number variations (CNVs) have been identified as another important structural variation of genome. In recent years, a large amount of CNVRs have been identified in humans and animals. However, association and dosage effects studies of CNVs are very limited. Apolipoprotein L3 (APOL3) gene plays a central role in modulating gene transcription and is located within a CNVR that encompasses quantitative trait locis (QTLs) for economic traits like meat quality. Herein, we analyzed the CNV polymorphism of APOL3 in 421 individuals from five distinct cattle breeds, and then correlated their genotypes with growth traits. Association analysis revealed that the APOL3 CNV was significantly associated with hip height and cannon circumference of Xianan (XN) cattle (P < .01), and visibly associated with body slanting length and hucklebone width of Pinan (PN) cattle (P < .05). Overall, the data provide evidence for the functional role of APOL3 CNV and a basis for future applications in cattle breeding.

RNA-Seq analysis reveals the potential molecular mechanisms of daidzein on adipogenesis in subcutaneous adipose tissue of finishing Xianan beef cattle.

Daidzein has been reported to be effective in regulating lipid metabolism in animals. However, the molecular mechanisms of daidzein on adipogenesis in beef cattle are not yet reported and the results of daidzein on affecting lipid metabolism in other species have been conflicting. High-throughput sequencing of mRNA (RNA-Seq) technology was performed to elucidate the underlying molecular mechanisms of daidzein on adipogenesis in subcutaneous adipose tissue of finishing Xianan beef cattle. A total of 893 differentially expressed genes (DEGs) were identified by differential expression analysis, among which 405 genes were upregulated and 488 genes were downregulated. Bioinformatics analysis suggested that these DEGs were significantly enriched to the pathways related to lipid metabolism including ECM-receptor interaction, Glycolysis/Gluconeogenesis and Hedgehog signalling pathway. Daidzein significantly affected the candidate genes (Shh, Pec, Gli, Wnt6, DLK, IGFBP2, ID3 and C/EBPE) related to adipocyte differentiation. Besides, daidzein improved the ability of subcutaneous adipocytes in synthesizing triglycerides by directly using the long-chain fatty acids and enhanced the efficiency of triglyceride synthesis of subcutaneous adipocytes in Xianan steers. In conclusion, daidzein plays a positive role not only in adipogenic differentiation, but also in triglyceride synthesis in subcutaneous adipose tissue of Xianan beef cattle.

LncRNA-MEG3 promotes bovine myoblast differentiation by sponging miR-135.

Long noncoding RNA maternally expressed gene 3 (lncRNA-MEG3) is an important regulator in multiple biological functions. However, lncRNA-MEG3's function in cattle growth and regulatory mechanism on bovine skeletal muscle development has not yet been well studied. In this project, we first investigated lncRNA-MEG3's expression profile and detected that it was highly expressed in bovine skeletal muscle tissue and its RNA level was kept increasingly during the early phase of bovine primary myoblast differentiation. Using luciferase reporter assays, we identified the lncRNA-MEG3 core promoter containing putative transcription factor binding site for myocyte enhancer factor 2C (MEF2C). Interestingly, we found that LncRNA-MEG3 could significantly upregulate and downregulate myosin heavy chain ( MHC), myogenin ( MyoG), and MEF2C through overexpression and RNAi strategies, respectively. Using luciferase reporter assays, we also verified lncRNA-MEG3 as a miR-135 sponge. Overexpression of miR-135 markedly inhibited the wild type of lncRNA-MEG3, but not the mutant lncRNA-MEG3 reporter. The luciferase activity of miR-135 sensor could be rescued by lncRNA-MEG3 via competing for miRNA-135. In addition, the luciferase activity of MEF2C was significantly upregulated by the wild type of lncRNA-MEG3. This study, for the first time, revealed that lncRNA-MEG3 could promote bovine skeletal muscle differentiation via interacting with miRNA-135 and MEF2C. The results were valuable for further studies and applications of lncRNA related roles in beef molecular breeding.

miR-148a-3p regulates proliferation and apoptosis of bovine muscle cells by targeting KLF6.

Skeletal muscle development is regulated by a series of regulatory factors, and also including noncoding RNA, especially microRNAs (miRNAs). Recently, miR-148a has been found to be involved in murine C2C12 differentiation by targeting ROCK1. However, the function of miR-148a-3p for the proliferation and apoptosis of bovine muscle cells has not been determined. In this study, we found that miR-148a-3p was highly expressed in fetal bovine skeletal muscle and exhibited a decreasing trend in muscle cells during its growth phase. Functional studies indicated that gain of miR-148a-3p inhibited the proliferation of bovine muscle cells and promoted apoptosis. Conversely, interference with miR-148a-3p inhibitor promoted muscle cell proliferation and inhibited its apoptosis. Mechanistically, KLF6 was confirmed as a new potential target gene of miR-148a-3p by TargetScan software prediction and the dual-luciferase assay verification. Additionally, after a gain or loss of KLF6, the function of KLF6 for muscle cell proliferation and apoptosis was opposite to that of miR-148a-3p. Collectively, these findings proposed a novel avenue whereby miR-148a-3p impeded bovine myoblast cell proliferation and promoted apoptosis through the posttranscriptional downregulation of KLF6.

Genetic Variant of MYLK4 Gene and its Association with Growth Traits in Chinese Cattle.

As a member of MYLK family, MYLK4 gene may play a vital role in muscle development. In this study, one novel single-nucleotide polymorphism (SNP) was identified the bovine MYLK4 by sequencing pooled DNA samples (pool-Seq) and forced polymerase chain reaction-restriction fragment length polymorphism (forced PCR-RFLP) methods. Overall, we reported one mutation (SNP1) in the intron 10 region within the bovine MYLK4 gene in 559 individuals representing five main cattle breeds from China (Nanyang, NY; Qinchuan; Jiaxian, JX; Pinan cattle; and Caidamu cattle, CDM). Genotype AA and allele A were predominant in the QC, PN, and XN populations. Association analysis with growth traits in the QC breed showed that the animals with genotype GG had significantly greater chest breadth and hip width (P < 0.05). Meanwhile, the genotype GG was strongly associated with withers height and body length than those with genotype AA (P < 0.01 or P < 0.05) at 12 months in the NY breed. These statistical results exhibited that the MYLK4 gene might be a potential candidate gene to improve cattle's growth traits, and the SNP could be used as molecular markers in early marker-assisted selection (MAS) in beef cattle breeding program.

CircFUT10 reduces proliferation and facilitates differentiation of myoblasts by sponging miR-133a.

Circular RNAs (circRNAs) have been identified in various tissues and cell types from human, monkey, porcine, and mouse. However, knowledge on circRNAs in bovine muscle development is limited. We downloaded and analyzed the circRNAs sequencing data of bovine skeletal muscle tissue, and further characterized the role of a candidate circRNA (circFUT10) in muscle development. Quantitative real-time PCR (qPCR) and Western blot assays were used to confirm the expression of genes involved in myoblasts differentiation and proliferation. Flow cytometry was performed to assess cell cycle distribution and cell apoptosis. EdU incorporation and CCK-8 assay were performed to demonstrate cell proliferation. We demonstrated that circFUT10 was highly (but differentially) expressed in embryonic and adult skeletal muscle tissue. circFUT10 induced bovine primary myoblasts differentiation and increased the expression of MyoD, MyoG, and MyhC in mRNA and protein levels. circFUT10 increased the number of myoblasts in the G0/G1 phase of the cell cycle, and decreased the proportion of cells in the S-phase. circFUT10 inhibited the proliferation of myoblasts and promoted them apoptosis. Via a luciferase screening assay, circFUT10 is observed to sponge to miR-133a with three potential binding sites. Specifically, we show that circFUT10 regulated myoblasts differentiation and cell survival by directly binding to miR-133a and inhibiting miR-133a activity. Modulation of circFUT10 expression in muscle tissue may emerge as a potential target in breeding strategies attempting to control muscle development in cattle.

Detection of Insertions/Deletions Within SIRT1, SIRT2 and SIRT3 Genes and Their Associations with Body Measurement Traits in Cattle.

Growth traits are complex quantitative traits controlled by numerous candidate genes, and they can be well-evaluated using body measurement traits. As the members of the nicotinamide adenine dinucleotide-dependent family of histone deacetylases, class I sirtuin genes (including SIRT1, SIRT2 and SIRT3) play crucial roles in regulating lipid metabolism, cellular growth and metabolism, suggesting that they are potential candidate genes affecting body measurement traits in animals. Hence, the objective of this work aimed to detect novel insertions/deletions (indels) of SIRT1, SIRT2 and SIRT3 genes in 955 cattle belonging to five breeds, as well as to evaluate their effects on body measurement traits. Herein, the novel 12-bp indel of SIRT1 gene, the 7-bp indel of SIRT2 gene and the 26-bp indel of SIRT3 gene were firstly reported, respectively. The association analysis indicated that the indels within SIRT1 and SIRT2 genes were significantly associated with body measurement traits such as body weight, chest circumference, height at hip cross, hip width, body height, etc. (P < 0.05 or P < 0.01). Therefore, based on these findings, the two novel indel variants within bovine SIRT1 and SIRT2 genes could be considered as potential molecular markers for growth traits in cattle selection practices and breeding.

The developmental transcriptome sequencing of bovine skeletal muscle reveals a long noncoding RNA, lncMD, promotes muscle differentiation by sponging miR-125b.

Pervasive transcription of the mammalian genome generates numerous long noncoding RNAs (lncRNAs), which are of crucial importance in diverse biological processes. Recent advances in high throughput sequencing technology have helped to accelerate the pace of lncRNA discovery. However, no study on the overall expression patterns of lncRNAs during muscle development has been conducted. We reported here the first analysis of lncRNA landscape in bovine embryonic, neonatal and adult skeletal muscle using Ribo-Zero RNA-Seq, a technology which can capture both poly(A)+ and poly(A)- transcripts. We finally defined 7692 high-confidence lncRNAs and uncovered 401 lncRNAs differentially expressed among three developmental stages, including lncMD, a novel muscle-specific lncRNA which is gradually up-regulated during myoblast differentiation. lncMD overexpression upregulated, whereas lncMD silencing decreased the expression of two well-established myogenic markers, myosin heavy chain (MHC) and myogenin (MyoG). In-depth analyses showed that lncMD acts as a molecular sponge for miR-125b and that insulin-like growth factor 2 (IGF2) is a direct target of miR-125b in cattle. Moreover, lncMD level was positively correlated with IGF2 mRNA level in bovine muscle tissues, a vital corollary to ceRNA function. Altogether, our research showed that lncMD acts as a ceRNA to sequester miR-125b, leading to heightened IGF2 expression and thus promotes muscle differentiation. Our findings also complement the reference genome annotation of cattle, which will likely be useful for further functional lncRNA cloning and more comprehensive studies on lncRNA regulation in muscle development.

Combined Haplotypes of CaSR Gene Sequence Variants and Their Associations with Growth Traits in Cattle.

The calcium-sensing receptor (CaSR) is a Class C G-protein coupled receptor that regulates food intake and assimilation. However, studies on the relationship between CaSR gene and growth traits in cattle are deficient. The aim of this study was to examine the association of the CaSR polymorphism with growth traits in cattle breeds. Four novel single nucleotide polymorphisms (SNPs) and one previously reported SNP (NC_007299.5: g.67630865T>C, 67638409G>C, 67660395G>C, 67661546C>G, and 67661892A>C) were identified in the bovine CaSR gene using DNA sequencing and PCR-SSCP methods in 520 individuals from three representative breeds. The three SNP P4_2, P7_1, and P7_4 in LX, QC, and JX cattle populations belonged to intermediate genetic diversity (0.25 < PIC < 0.5). In addition, we evaluated the haplotype frequency and linkage disequilibrium coefficient of five sequence variants in the three cattle breeds. LD and haplotype structure of CaSR were different between breeds. LD analysis showed that the P4_2 and P7_4 loci were in complete LD in JX cattle population (r2 = 0.99 and D' = 1). Only 11 haplotypes were listed except for those with a frequency of <0.03. Hap1 (-TGGGC-) had the highest haplotype frequencies in LX (27.30%), Hap6(-TGGCC-) had the highest haplotype frequencies in QC (21.70%) and JX (32.30%). Association analysis indicated that P2, P4_2, and P7_4 loci were all significantly associated with growth traits and combined genotype TTGCGC was highly significantly associated with Chest circumference and body weight than the other genotype in JX cattle population. The results of this study suggest that the CaSR gene possibly is a strong candidate gene that affects growth traits in the Chinese cattle breeding program.

Correlation between ZBED6 Gene Upstream CpG Island methylation and mRNA expression in cattle.

DNA methylation is essential for the regulation of gene expression and important roles in muscle development. To assess the extent of epigenetic modifications and gene expression on the differentially methylated region (DMR) in ZBED6, we simultaneously examined DNA methylation and expression in six tissues from two different developmental stages (fetal bovine and adult bovine). The DNA methylation pattern was compared using bisulfite sequencing polymerase chain reaction (BSP) and combined bisulfite restriction analysis (COBRA). The result of quantitative real-time PCR (qPCR) analysis showed that ZBED6 has a broad tissue distribution and is highly expressed in adult bovine (P < 0.05 or P < 0.01). The DNA methylation level was significantly different in liver, lung and spleen between the two cattle groups (P < 0.05 or P < 0.01). The adult bovine group exhibited a significantly higher mRNA level and lower DNA methylation level than the fetal bovine group in liver, lung, and spleen. No significant association was detected between DNA methylation level and muscle, heart, and kidney at two different stages. In this study, the statistical analyses indicated that DNA methylation patterns are associated with mRNA level in some tissues, these results may be a useful parameter to investigate muscle developmental in cattle and as a model for studies in other species, potentially contributing to an improvement of growth performance selection in beef cattle breeding program.

Effect of Genetic Variations within the I-mfa Gene on the Growth Traits of Chinese Cattle.

I-mfa (inhibitor of the MyoD family a) is a transcription modulator that binds to MyoD family members and inhibits their transcriptional activities. It is highly expressed in the sclerotome and plays an important role in the patterning of the somite early in development. In this study, the polymorphisms of the bovine I-mfa gene were detected by polymerase chain reaction-single stranded conformational polymorphism (PCR-SSCP) and DNA pool sequencing methods in 541 individuals from three Chinese cattle breeds. The results showed that P3 locus had two novel complete linked single nucleotide polymorphisms (NC_007324.4:g.12284A>G and g.12331T>C), resulting in a missense mutation p.S(AGC)113G(GGC) and a synonymous mutation p.H(CAT)128H(CAC), respectively. P4 locus had a novel SNP (NC_007324.4: g.16432C>A), which resulted in a nonsense mutation p.C(TGC)241X(TGA). The statistical analyses indicated that the three SNPs, are associated with the phenotypic traits in Luxi (LX), Qinchuan (QC), and Jiaxian (JX) cattle population (P < 0.05 or P < 0.01). The mutant-type variants were superior for growth traits; the heterozygote diplotype was associated with higher growth traits compared to wild-type homozygote. Our results provide evidence that polymorphisms in the I-mfa gene are associated with growth traits and may be used for marker-assisted selection in beef cattle breeding program.

miR-30-5p Regulates Muscle Differentiation and Alternative Splicing of Muscle-Related Genes by Targeting MBNL.

MicroRNAs (miRNAs), a class of single stranded, small (~22 nucleotides), non-coding RNAs, play an important role in muscle development. We focused on the role of the miR-30-5p family during bovine muscle development from previous high-throughput sequencing results and analyzed their expression profiles. MHC and MyoG mRNAs expression as well as their proteins were suppressed in differentiated C2C12 cells, suggesting the importance of miR-30-5p in muscle development. MBNL, the candidate target of miR-30-5p, is an alternative splicing regulation factor. MBNL1 and MBNL3 have opposite effects on muscle differentiation. Our results confirmed that miR-30a-5p and miR-30e-5p repress the expression of MBNL1, MBNL2 and MBNL3, whereas miR-30b-5p inhibits MBNL1 and MBNL2 expression. This provides direct evidence that MBNL expression can be flexibly regulated by miR-30-5p. Previous studies showed that MBNL1 promotes exon inclusion of two muscle-related genes (Trim55 and INSR). Through RNA splicing studies, we found that miR-30-5p had an effect on their alternative splicing, which means miR-30-5p via MBNL1 could be integrated into muscle signaling pathways in which INSR or Trim55 are located. In conclusion, miR-30-5p could inhibit muscle cell differentiation and regulate the alternative splicing of Trim55 and INSR by targeting MBNL. These results promote the understanding of the function of miRNAs in muscle development.

Two novel polymorphisms of bovine SIRT2 gene are associated with higher body weight in Nanyang cattle.

Identification of polymorphisms associated with economic traits is important for successful marker-assisted selection in cattle breeding. The family of mammalian sirtuin regulates many biological functions, such as life span extension and energy metabolism. SIRT2, a most abundant sirtuin in adipocytes, acts as a crucial regulator of adipogenic differentiation and plays a key role in controlling adipose tissue function and mass. Here we investigated single nucleotide polymorphisms (SNPs) of bovine SIRT2 in 1226 cattle from five breeds and further evaluated the effects of identified SNPs on economically important traits of Nanyang cattle. Our results revealed four novel SNPs in bovine SIRT2, one was located in intronic region and the other three were synonymous mutations. Linkage disequilibrium and haplotype analyses based on the identified SNPs showed obvious difference between crossbred breed and the other four beef breeds. Association analyses demonstrated that SNPs g.17333C > T and g.17578A > G have a significantly effect on 18-months-old body weight of Nanyang population. Animals with combined genotype TTGG at the above two loci exhibited especially higher body weight. Our data for the first time demonstrated that polymorphisms in bovine SIRT2 are associated with economic traits of Nanyang cattle, which will be helpful for future cattle selection practices.

SIRT1 gene polymorphisms are associated with growth traits in Nanyang cattle.

Growth is under complex genetic control and uncovering the molecular mechanisms how the genes and polymorphisms affect economic growth traits, are important for successful marker-assisted selection and more efficient management strategies in commercial cattle populations. SIRT1 is a NAD(+)-dependent deacetylase that belongs to the class III histone deacetylases. It plays an important role in numerous fundamental cellular processes including gene silencing, DNA repair, and metabolic regulation. In addition, SIRT1 acts as an inhibitor of adipogenesis and has been associated with body weight regulation. The objective of the present study was to identify single nucleotide polymorphisms (SNPs) of bovine SIRT1 using 1255 animals representing the five main Chinese breeds and to determine if these SNPs are associated with economically important traits in Nanyang cattle. The approach consisted of resequencing SIRT1 using a panel of DNA from unrelated animals of five different breeds and the process revealed five novel SNPs. SNPs g.17324T>C and g.17491G>A exhibited a high degree of linkage disequilibrium in all tested breeds. Seven major haplotypes accounting for 91.2% of the alleles were observed and the haplotype 'GCCGA' was the most common haplotype in NY, QC, LX and JX breeds. An association analysis was performed between the five SNPs and six performance traits. SNP g.-274C>G was demonstrated to have a strong effect on 24-months-old body weight and g.17379A>G polymorphism was related to 6 and 12-months-old body weight in NY population, although these effects did not remained significant after the Bonferroni correction. Our results provide evidence that polymorphisms in SIRT1 are associated with growth efficiency traits, and may be used for marker-assisted selection and management in feedlot cattle.

miR-183/96/182 Cluster Regulates the Development of Bovine Myoblasts through Targeting FoxO1.

Muscle development is an important factor affecting meat yield and quality and is coordinated by a variety of the myogenic genes and signaling pathways. Recent studies reported that miRNA, a class of highly conserved small noncoding RNA, is actively involved in regulating muscle development, but many miRNAs still need to be further explored. Here, we identified that the miR-183/96/182 cluster exhibited higher expression in bovine embryonic muscle; meanwhile, it widely existed in other organizations. Functionally, the results of the RT-qPCR, EdU, CCK8 and immunofluorescence assays demonstrated that the miR-183/96/182 cluster promoted proliferation and differentiation of bovine myoblast. Next, we found that the miR-183/96/182 cluster targeted FoxO1 and restrained its expression. Meanwhile, the expression of FoxO1 had a negative correlation with the expression of the miR-183/96/182 cluster during myoblast differentiation. In a word, our findings indicated that the miR-183/96/182 cluster serves as a positive regulator in the proliferation and differentiation of bovine myoblasts through suppressing the expression of FoxO1.

miRNA profiling in intrauterine exosomes of pregnant cattle on day 7.

Intrauterine exosomes have been identified to be involved in the embryo development and implantation. The aim of this study was to explore the role of miRNAs in intrauterine exosomes in bovine pregnancy. Intrauterine exosomes were collected from uterine flushing fluids of three donor and three recipient Xianan cows 7 days after fertilization. Intrauterine exosomes miRNAs were extracted and the exosomal miRNAs expression levels were analyzed. Sixty miRNAs differed significantly in their amounts between donors and recipients (p-value < 0.05, |log2(FoldChange)| > 1). Twenty-two miRNAs were upregulated and 38 downregulated in the group of donor cows. The bta-miR-184 was the most significant (P Benjamini-Hochberg < 0.001). A total of 9,775 target genes were predicted using the 60 miRNAs. GO and KEGG analysis showed that the target genes were enriched in several biological processes or pathways associated with embryo implantation and endometrial development, such as cell adhesion, cell junction, focal adhesion, and Rap1 signaling pathway. Our findings suggest that, in cattle early pregnancy stage, these differently expressed miRNAs in intrauterine exosomes involved in embryo implantation and endometrial development, which may exert a significant effect and influence the uterine microenvironment for embryo implantation. These results could provide reference for screening and exploring the intrauterine exosomal miRNA affecting embryo implantation.

A novel missense mutation (rs464874590) within BoLA-DOB gene associated with the heat-resistance in Chinese 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.