Identification and bioinformatics analysis of differentially expressed milk exosomal microRNAs in milk exosomes of heat-stressed Holstein cows.

In summer, heat stress is one of the primary reasons for the compromised health and low milk productivity of dairy cows. Hyperthermia affects milk synthesis and secretion in the mammary glands of dairy cows. As molecules for intercellular communication, milk-derived exosomes carry genetic material, proteins, and lipids, playing a crucial role in mammary tissue growth and milk synthesis in dairy cows. The aim of this study was to explore the milk exosomal miRNA profile of heat-stressed and normal Holstein cows. We isolated and identified milk exosomes to screening for differentially expressed miRNAs using small RNA sequencing. Then, TargetScan and miRanda algorithms were used to predict the putative targets of the differentially expressed miRNAs, whereas GO and KEGG pathway enrichment analyses were performed for the differentially expressed miRNA-target genes. Our results showed that 215 miRNAs were significantly differentially expressed in heat-stressed milk exosomes, of which one was upregulated and 214 were significantly downregulated. GO and KEGG enrichment analyses indicated that differentially expressed miRNAs might play a role in apoptosis, autophagy, and the p38 MAPK pathway. qRT-PCR assay verified that the expression of miRNAs was consistent with the sequencing results, warranting further verification of their specific targets of action. In conclusion, changes in the miRNA expression profile of milk exosomes indicated the role of exosomal miRNAs in regulating heat stress resistance and apoptosis in dairy cows. Our results suggested that milk-derived exosomal miRNAs could increase mammary gland resistance to heat stress, thereby enhancing milk synthesis in dairy cows.

Depletion of serum-derived exosomes aggravates heat stress-induced damage of bovine mammary epithelial cells.

BACKGROUND: Exosomes are involved in intercellular communication, affecting many physiological and pathological process. The present study evaluated the effects of serum exosomes on the function of bovine mammary epithelial cells (BMECs) and milk synthesis under heat stress. METHODS AND RESULTS: We cultured the BMECs in fetal bovine serum (FBS) or exosome-free FBS medium and examined, their viability using CCK-8 kit. The results showed that culturing the cells in an exosome-free medium decreased viability and increased the levels of reactive oxygen species. The BMECs cultured in the exosome-free medium had reduced mitochondrial membrane potential, decreased manganese superoxide dismutase activity, and disrupted mitochondrial dynamics. They exhibited apoptosis due to upregulated Drp1, Fis1, Bax and HSP70. Lastly, we observed downregulation of milk fat and lactoprotein-related genes: mTOR, PPARγ, p-mTOR and ADD1 and SREBP1, ELF5, and CSN2, respectively, after culturing the cells in an exosome-free medium. These negative effects of the exosome-free medium on the BMECs could be further reinforced under heat stress. CONCLUSION: Our results demonstrated that exosomes from serum are critical for maintaining the normal function of BMECs.

SIRT4 Expression Ameliorates the Detrimental Effect of Heat Stress via AMPK/mTOR Signaling Pathway in BMECs.

Sirtuin 4 (SIRT4), a member of the SIRT family, has been reported to be a key factor involved in antioxidant defense in mitochondria. This study aimed to explore the potential molecular mechanism via which SIRT4 regulates heat stress-induced oxidative stress and lactoprotein synthesis in bovine mammary epithelial cells (BMECs). Our results showed that SIRT4 was significantly decreased in heat stressed mammary tissue. Depletion of SIRT4 in BMECs induced the generation of ROS, which, as exhibited by the decreased activity of antioxidant enzymes, changed mitochondrial morphology through mediating protein and mRNA levels related to mitochondrial fission and fusion. Moreover, we found that depletion of SIRT4 or stress conditions inhibited the expression of milk proteins, as well as lipid and glucose synthesis-related genes, and activated the AMPK/mTOR signaling pathway. Increased SIRT4 expression was found to have the opposite effect. However, blocking the AMPK/mTOR signaling pathway could inhibit the regulatory function of SIRT4 in milk synthesis-related gene expression. In summary, our results suggest that SIRT4 may play critical roles in maintaining mammary gland function by regulating the AMPK/mTOR signaling pathway in dairy cows, indicating that SIRT4 may be a potential molecular target for curing heat stress-induced BMEC injury and low milk production in dairy cows.

S-allyl cysteine ameliorates heat stress-induced oxidative stress by activating Nrf2/HO-1 signaling pathway in BMECs.

Heat stress-induced oxidative stress in bovine mammary epithelial cells (BMECs) threatens the normal growth and development of bovine mammary tissue, resulting in lower milk production of dairy cows. The aim of the present study is to investigate the protective effects of S-allyl cysteine (SAC), an organosulfur component extracted from aged garlic, on heat stress-induced oxidative stress and apoptosis in BMECs and to explore its underlying mechanisms. Our results showed that heat stress treatment considerably decreased cell viability, whereas SAC treatment dose-dependently restored cell viability of BMECs under heat-stress conditions. In addition, SAC protected BMECs from heat stress-induced oxidative damage by inhibiting the excessive accumulation of reactive oxygen species (ROS) and increasing the activity of antioxidant enzymes. It also inhibited heat stress-induced apoptosis by reducing the ratio of Bax/Bcl-2 and blocking proteolytic the cleavage of caspase-3 in BMECs. Interestingly, we found that the protective effect of SAC on heat stress-induced oxidative stress and apoptosis was dependent on the nuclear factor erythroid 2-related factor 2 (Nrf2)/heme oxygenase-1 (HO-1) signaling pathway. SAC promoted the Nrf2 nuclear translocation in heat stress-induced BMECs. The results were also validated by Nrf2 and Keap1 knockdown experiments further demonstrating that Nrf-2 was indeed involved in the protective effect of SAC on heat stress-induced oxidative damage and apoptosis. In summary, our results showed that SAC could protect BMECs from heat stress-induced injury by mediating the Nrf2/HO-1 signaling pathway, suggesting that SAC could be considered as a therapeutic drug for attenuating heat stress-induced mammary gland diseases.

Dihydromyricetin attenuates heat stress-induced apoptosis in dairy cow mammary epithelial cells through suppressing mitochondrial dysfunction.

It is well known that the dairy cow production is very sensitive to environmental factors, including high temperature, high humidity and radiant heat sources. High temperature-induced heat stress is the main environmental factor that causes oxidative stress and apoptosis, which affects the development of mammary glands in dairy cows. Dihydromyricetin (DMY) is a nature flavonoid compound extracted from Ampelopsis grossedentata; it has been shown to have various pharmacological functions, such as anti-inflammation, antitumor and liver protection. The present study aims to evaluate the protective effect of DMY on heat stress-induced dairy cow mammary epithelial cells (DCMECs) apoptosis and explore the potential mechanisms. The results show that heat stress triggers heat shock response and reduces cell viability in DCMECs; pretreatment of DCMECs with DMY (25 µM) for 12 h significantly alleviates the negative effects of heat stress on cells. DMY can provide cytoprotective effects by suppressing heat stress-caused mitochondrial membrane depolarization and mitochondrial dysfunction, Bax and Caspase 3 activity, and modulation of oxidative enzymes, thereby preventing ROS production and apoptosis in DCMECs. Importantly, DMY treatment could attenuate heat stress-induced mitochondrial fragmentation through mediating the expression of mitochondrial fission and fusion-related genes, including Dynamin related protein 1 (Drp1), Mitochondrial fission 1 protein (Fis1), and Mitofusin1, 2 (Mfn1, 2). Above all, our findings demonstrate that DMY could protect DCMECs against heat stress-induced injury through preventing oxidative stress, the imbalance of mitochondrial fission and fusion, which provides useful evidence that DMY can be a promising therapeutic drug for protecting heat stress-induced mammary glands injury and mastitis.

[Functional haplotypes of INCENP affect promoter activity and bovine semen quality].

To explore the association between single nucleotide polymorphisms (SNPs) in the promoter region of the inner centromere protein (INCENP) gene and bovine semen quality, the haplotypes in 250 Chinese Holstein bulls were detected using PCR-RFLP method in this study. Two SNPs (g.-556 G>T, rs 136823901 and g.-692 C>T, rs 211010999) and three haplotypes (CG, TT, TG) were identified in the promoter region of INCENP. The genotype frequency and allele frequency of these two SNPs as well as the correlation between different SNP haplotype combinations and bovine semen quality were then analyzed. Our results showed that fresh sperm motility of the GT genotype was significantly higher than that of the GG genotype (P<0.05) at the SNP site g.-556 G>T, while fresh and frozen-thawed sperm motilities of the haplotype combinations H1H1(CCGG), H1H3(CTGT), H2H3(TTGT) and H3H3(TTTT) were significantly higher than that of H1H2 (P<0.05). To further study the possible mechanisms by which g.-556 G>T and g.-692 C>T affect semen quality, three haplotype plasmids were respectively transfected into MLTC-1 cells. The TG haplotype demonstrated the highest luciferase activity, suggesting that g.-556 G>T and g.-692 C>T are functional mutations which could regulate INCENP gene expression by affecting promoter activity and thus affect semen quality.

PLAG1 g.8795C>T Mutation Regulates Early Body Weight in Hu Sheep by Weakening miR-139 Binding.

Sheep birth and weaning weights indicate their growth and survival. Thus, identifying molecular genetic markers for early body weight is important in sheep breeding. Pleomorphic adenoma gene 1 (PLAG1) is important for regulating birth weight and body length in mammals; however, its relationship with sheep body weight remains unknown. Here, the 3'-untranslated region (3'-UTR) of the Hu sheep PLAG1 gene was cloned, single nucleotide polymorphisms (SNPs) were screened, genotype-early body weight relationships were analyzed, and the possible molecular mechanism was explored. PLAG1 3'-UTR sequences with five forms of base sequences plus poly(A) tails were detected in Hu sheep and the g.8795C>T mutation was identified. Luciferase reporter assay indicated that the g.8795C>T mutation influenced PLAG1 post-transcriptional activity. miRBase prediction showed that the g.8795C>T mutation was located in the miR-139 seed sequence binding region, and miR-139 overexpression significantly decreased both PLAG1-CC and PLAG1-TT activities. Moreover, the luciferase activity of PLAG1-CC was significantly lower than that of the PLAG1-TT, but miR-139 inhibition substantially increased both PLAG1-CC and PLAG1-TT luciferase activities, suggesting that PLAG1 is the target gene of miR-139. Thus, the g.8795C>T mutation upregulates PLAG1 expression by weakening its binding with miR-139, promoting PLAG1 expression, and increasing Hu sheep birth and weaning weights.

miR-27a-3p relieves heat stress-induced mitochondrial damage and aberrant milk protein synthesis through MEK/ERK pathway in BMECs.

With global warming, heat stress has become a primary factor that compromises the health and milk quality of dairy cows. Here, we investigated the function and underlying regulatory mechanism of miR-27a-3p in bovine mammary epithelial cells (BMECs) under heat-stress conditions. The current study showed that miR-27a-3p could prevent heat stress-induced BMEC oxidative stress and mitochondrial damage by regulating the balance between mitochondrial fission and fusion processes. Importantly, we found that miR-27a-3p could increase cell proliferation under heat stress conditions by regulating the MEK/ERK pathway and cyclin D1/E1. Interestingly, miR-27a-3p is also involved in the regulation of milk protein synthesis-related protein expression, such as CSN2 and ELF5. Inhibition of the MEK/ERK signaling pathway by AZD6244 blocked the regulatory function of miR-27a-3p in cell proliferation and milk protein synthesis in BMECs under heat stress conditions. Our findings demonstrated that miR-27a-3p protects BMECs from heat stress-induced oxidative stress and mitochondrial damage through the MEK/ERK pathway, thereby promoting BMECs proliferation and lactation in dairy cows. The potential regulatory mechanism of miR-27a-3p in attenuating heat stress-induced apoptosis and lactation defect in BMECs.

[Polymorphism of POU1F1 gene and PRL gene and their combined effects on milk performance traits in Chinese Holstein cattle].

Three novel SNPs were found by DNA sequencing, PCR-RFLP and CRS-PCR methods were used for genotyping in 979 Chinese Holstein cattle. One SNP, G1178C, was identified in exon 2 of POU1F1 gene. Two novel SNPs, A906G and A1134G, were identified in 5'-flanking regulatory region (5'-UTR) of PRL gene. The association between polymorphisms of the two genes and milk performance traits were analyzed with PROC GLM of SAS. The results showed that GC genotype at 1178 locus of POU1F1 gene was advantageous for milk yield, milk protein yield, and milk fat yield. AG genotype at 906 locus was advantageous for milk yield. There was no significant difference between 1134 locus and milk performance traits of 5'-UTR of PRL gene. Analysis of genotype combination effect on milk production traits showed that the effect of combined genotype was not simple sum of single genotypes and the effects of gene pyramiding seemed to be more important in molecular breeding.

Heat stress induces apoptosis through disruption of dynamic mitochondrial networks in dairy cow mammary epithelial cells.

Heat stress-induced reductions in milk yield and the dysfunction of mammary glands are economically important challenges that face the dairy industry, especially during summer. The aim of the present study is to investigate the effects of heat stress on mitochondrial function by using dairy cow mammary epithelial cells (DCMECs) as an in vitro model. Live cell imaging shows that the mitochondria continually change shape through fission and fusion. However, heat stress induces the fragmentation of mitochondria, as well as the decreased of ATP level, membrane potential, and anti-oxidant enzyme activity and the increased of respiratory chain complex I activity. In addition, the cytosolic Ca2+ concentration and cytochrome c expression (Cyto-c) were increased after heat stress treatment. Both qRT-PCR and western blot analysis indicate that mitofusin1/2 (Mfn1/2) and optic atrophy protein-1 (Opa-1) are downregulated after heat stress, whereas dynamin-related protein 1 (Drp1) and fission 1 (Fis-1) are upregulated, which explains the observed defect of mitochondrial network dynamics. Accordingly, the present study indicated that heat stress induced the dysfunction of DCMEC through disruption of the normal balance of mitochondrial fission and fusion.

[Genetic polymorphism of Nramp1 gene and correlation with mastitis in Holstein cattle].

In this research, PCR-SSCP technique was used to analyze the polymorphisms of the exon 11 of Nramp1 gene in Chinese Holstein cattle (n=344), and correlation between polymorphisms of Nramp1 with somatic cell score (SCS) and milk production traits was analyzed. The results show that three genotypes namely AA, BB, and AB were detected. Allele A was predominant and the frequencies of alleles A and B were estimated to be 0.767 and 0.233, respectively. Chi-square test indicated that the polymorphic locus in Chinese Holstein fitted Hardy-Weinberg equilibrium (P>0.05). Sequencing analysis showed two polymorphic sites at positions 200 bp (C/G) and 254 bp (T/G), which resulted in amino acid alteration Ala356Pro and Leu374Met. The least squares means of SCS in Holstein cattle was lower for genotype AA than that for genotypes AB and BB (P<0.05). The least squares means of milk yield of genotype AA and AB were higher than that for genotype BB (P<0.05, P<0.01, respectively). Genotype AA was beneficial to mastitis resistance. This suggested that Nramp1 may be a candidate gene responsible for mastitis in Holstein cattle.

[Genetic polymorphism of HSP70-1 gene and its correlation with resistance to mastitis in Chinese Holstein].

The polymorphisms of HSP70-1 gene in 253 Chinese Holstein dairy cows were studied, and the association between the polymorphisms and somatic cell score (SCS) were analyzed. PCR-SSCP, PCR-RFLP and DNA sequencing were used to investigate mutations in the coding region of HSP70-1 gene. The G-->A-->C mutation at 1 623 bp and G-->A mutation at 2 409 bp were found and both of them were silence mutations that caused no alteration in amino acid sequence. Chi-square test showed both loci were'nt at Hardy-Weinberg disequilibrium in Chinese Holstein. In the meanwhile, the association of 2 409 locus and SCS was not significant. However, the polymorphism at 1623 locus affected SCS significantly (P<0.05). The SCS of genotype CC was significantly lower than that of genotype AG and GG (P<0.05), so genotype CC was mastitis resistant. These results suggest that genotype CC of HSP70-1 gene may be used as a molecular and genetic marker to improve the phenotype of anti-mastitis in Chinese Holstein dairy cows.

[Genetic polymorphism in 5'-flanking region of the lactoferrin gene and its associations with mastitis in Chinese Holstein cows].

Bovine lactoferrin (LF) is a multifunctional glycoprotein found in milk, which acts mainly as a defense factor in the mammary gland. Polymorphism has been found in the bovine LF gene. However, there is no report on genetic polymorphism of LF gene and its associations with mastitis in dairy cattle. In this study, the promoter fragment of LF gene containing -926(G/A), -915(T/G), -478(/G), and +72(T/C) mutations were genotyped by the PCR-RFLP and CRS-PCR method. Two hundred and sixty-eight Chinese Holstein cows were screened. Least square linear model (LSM) analysis was applied to evaluate the associations of LF gene with somatic cell score (SCS). The results indicated that the SCS was significantly affected by -478(/G) and +72(T/C), but not by the other two loci (P >0.05). The SCS of cow with genotype AB in +72(T/C) position was significantly lower than that of genotype AA (P<0.01) or AB (P<0.05). In position -478(/G), the cow with genotype CC showed significantly lower SCS in contrast to cow with genotype CD and DD (P < 0.01). In conclusion, genotype AB in position +72(T/C) and genotype CC in position -478(/G) of LF gene were advantageous genotype, which can be used as candidate markers for mastitis resistance selection in dairy cattle.

[The development and application of method for detecting bovine complex vertebral malformation].

Complex vertebral malformation (CVM), a lethal autosomal recessive inherited defect in Holstein calves, was newly reported worldwide. The molecular cause of CVM was a substitution of guanine by thymine (G-->T) in a solute carrier family 35 member 3 gene (SLC35A3), encoding UDP-N-acetylglucosamine transporter. It was characterized by stillborn, abortion, and premature birth. The objective of this study was to study the actual carrier frequency of the CVM mutation in a population of Chinese Holstein (=Chinese Black-and-White) normal cattle. The normal 436 Holstein cows and 93 Holstein bulls were genotyping by using the Created Restriction Site PCR (CRS-PCR) and Allele-specific PCR (AS-PCR) methods. There were two bulls and one cow in three observed CVM-carriers. In the Holstein dairy cattle and Holstein bull population, the percentages of CVM carriers were estimated as 0.60% and 2.20% respectively. This study provided a more reliable and useful method for extensive screening of CVM and also offers a theoretical basis for molecular diagnosis in Holstein calves.

[Genetic polymorphism of k-casein gene exon 4 and its cor-relation with milk production traits in Chinese Holsteincows.].

K-casein gene was regarded as a candidate gene for milk production traits of cows. In this study, a 779 bp fragment of k-casein gene of Chinese Holstein was amplified by polymerase chain reaction (PCR), the polymorphisms of three loci of k-casein gene were detected by PCR-RFLP with restriction endonuclease Taq, Hind, Pst. After sequencing, T/C single nucleotide polymorphism (SNP) was identified at nucleotide 10 891C/A SNP was identified at nucleotide 10 927 and G/A SNP was identified at nucleotide 10 988 in exon4 of k-casein gene. Both alleles (A and B) of three loci were found in the population that showed low polymorphism. The gene frequencies of A and B were 86.03% and 13.97%, respectively. The genotype frequencies of AA, AB, and BB were 73.71%, 24.63%, and 1.66%, respectively. Statistical results of c2 test indicated that three polymorphism sites in the population fitted with Hardy-Weinberg equilibrium (P > 0.05). Meanwhile, the effect of polymorphism of k-casein gene on milk production traits was analyzed. The results indicated that in the three loci, the different genotype of k-casein gene had no significant influence on milk yield and milk protein percent (P > 0.05). The cows with genotypes BB and AB showed higher milk fat percent than those with genotype AA ( P < 0.05 ) ; with genotype AB showed higher fat protein ratio than those with genotype AA ( P < 0.05 ). The polymorphism of the three loci in the experimental population is closely linked. The conclusion is that k-casein B allele can be used as the molecular genetic markers of modifying milk fat percent in Chinese Holstein cows.

[Relationship between the melanocortin receptor 1 (MC1R) gene and the coat color phenotype in cattle].

The MC1R gene is related to coat color in cattle and to lactoprotein content in milk. Polymerase chain reaction-restriction fragment length polymorphism and DNA sequencing were used to analyze the MC1R of China Holstein Black-white Cattle, China Holstein Red-white Cattle, Luxi Yellow Cattle and Bohai Black Cattle. Three alleles (ED, E+ and e) were identified. Most China Holstein Black-white Cattle had ED and E+ alleles (ED = 0.12, E+ = 0.80). The same was true for Bohai Black Cattle (ED = 0.52, E+ = 0.47). On the other hand, most China Holstein Red-white Cattle and Luxi Yellow Cattle had the e allele (e/e = 0.95). E+/e genotype was present in China Holstein Red-white Cattle and Luxi Yellow Cattle. We reasoned that ED and E+ lad to the synthesis of the black pigment. Furthermore, we found an important SNP (Single Nucleotide Polymorphism) at 725 position of the MC1R coding region.

Deciphering Transcriptome and Complex Alternative Splicing Transcripts in Mammary Gland Tissues from Cows Naturally Infected with Staphylococcus aureus Mastitis.

Alternative splicing (AS) contributes to the complexity of the mammalian proteome and plays an important role in diseases, including infectious diseases. The differential AS patterns of these transcript sequences between the healthy (HS3A) and mastitic (HS8A) cows naturally infected by Staphylococcus aureus were compared to understand the molecular mechanisms underlying mastitis resistance and susceptibility. In this study, using the Illumina paired-end RNA sequencing method, 1352 differentially expressed genes (DEGs) with higher than twofold changes were found in the HS3A and HS8A mammary gland tissues. Gene ontology and KEGG pathway analyses revealed that the cytokine-cytokine receptor interaction pathway is the most significantly enriched pathway. Approximately 16k annotated unigenes were respectively identified in two libraries, based on the bovine Bos taurus UMD3.1 sequence assembly and search. A total of 52.62% and 51.24% annotated unigenes were alternatively spliced in term of exon skipping, intron retention, alternative 5' splicing and alternative 3' splicing. Additionally, 1,317 AS unigenes were HS3A-specific, whereas 1,093 AS unigenes were HS8A-specific. Some immune-related genes, such as ITGB6, MYD88, ADA, ACKR1, and TNFRSF1B, and their potential relationships with mastitis were highlighted. From Chromosome 2, 4, 6, 7, 10, 13, 14, 17, and 20, 3.66% (HS3A) and 5.4% (HS8A) novel transcripts, which harbor known quantitative trait locus associated with clinical mastitis, were identified. Many DEGs in the healthy and mastitic mammary glands are involved in immune, defense, and inflammation responses. These DEGs, which exhibit diverse and specific splicing patterns and events, can endow dairy cattle with the potential complex genetic resistance against mastitis.

Correction: Deciphering Transcriptome and Complex Alternative Splicing Transcripts in Mammary Gland Tissues from Cows Naturally Infected with Staphylococcus aureus Mastitis.

[This corrects the article DOI: 10.1371/journal.pone.0159719.].

Cloning and expression of goat interleukin-18 gene.

We isolated and sequenced a 480 bp cDNA encoding mature goat interleukin-18 (gIL-18) from alveolar macrophages and splenocytes activated with LPS by RT-PCR. The gIL-18 gene was cloned into pET32a (+) vectors and sequenced. Nucleotide sequence of gIL-18 shares high homology with cattle. Fusional expression with pET32a (+) of gIL-18 of about 38kD was obtained by SDS-PAGE analysis after induction by IPTG in the E. Coli BL21 expression system. The recombinant protein can induce IFN-gamma production in PBMC. The IL-18 mRNA was constitutively detected in goat alveolar macrophages with or without LPS, While, enhanced expression was detected in splenocytes and liver cells if treated by LPS, and can be weakly detected in Peripheral blood mononuclear cells (PBMCs) treated by activators. Significant deference of IL-18 mRNA level may reflect the capacity to produce mature IL-18 in such tissues.

Two novel SNPs in HSF1 gene are associated with thermal tolerance traits in Chinese Holstein cattle.

Heat-shock transcription factors (HSFs) play an important role in regulating heat stress response. The activation of heat-shock protein (HSP) genes is mediated by HSFs, which bind to promoters of HSP genes. In this research, two novel single nucleotide polymorphisms, T909C and G4693T, and their association with thermal tolerance were investigated in 951 Chinese Holstein cattle. Linkage disequilibrium and haplotype construction were analyzed using SHEsis software. Four haplotypes were constructed, and nine haplotype combinations were found. Potassium content in erythrocytes (PCE), decreased rate of milk production (R), rectal temperature (RT), and heat-tolerance coefficient (HTC) were selected for the thermotolerance index. Association analysis showed that thermal tolerance in Chinese Holstein cattle was significantly affected by T909C and G4693T. The PCE of cows with CC or TC genotype was lower than that of TT at the 909 position (p < 0.05). Cows with TT genotype had lower PCE (p < 0.01) and higher HTC (p < 0.05) at the 4693 position. Cows with H2H4 haplotype combination had lower PCE (p < 0.01), R (p < 0.05) and RT (p < 0.05) and higher HTC (p < 0.05) than those with H1H3 haplotype combination. Bioinformatic analysis predicted that the 4693 position was located in the microRNA-binding (bta-miR-484) region. Quantitative reverse transcription-polymerase chain reaction demonstrated that 4693-T mutation caused the disruption of microRNA target binding, resulting in the relief of the transcriptional repression, which, in turn, resulted in increased expression. Thus, the HSF1 gene is useful in dairy cattle thermal tolerant breeding.