ABSTRACT: Korean peninsular weather is rapidly becoming subtropical due to global warming. In summer 2018, South Korea experienced the highest temperatures since the meteorological observations recorded in 1907. Heat stress has a negative effect on Holstein cows, the most popular breed of dairy cattle in South Korea, which is susceptible to heat. To examine physiological changes in dairy cows under heat stress conditions, we analyzed the profiles circulating microRNAs isolated from whole blood samples collected under heat stress and non-heat stress conditions using small RNA sequencing. We compared the expression profiles in lactating cows under heat stress and non-heat stress conditions to understand the regulation of biological processes in heat-stressed cows. Moreover, we measured several heat stress indicators, such as rectal temperature, milk yield, average daily gain, and progesterone concentration. All these assessments showed that pregnant cows were more susceptible to heat stress than non-pregnant cows. Particularly, progesterone concentrations known to have maternal warming effects were at similar levels in non-pregnant cows but significantly increased in pregnant cows under heat stress conditions. The differentially expressed miRNAs and their putative target genes were analyzed in pregnant cows. Interestingly, we found that differentially expressed miRNAs (bta-miR-146b, bta-miR-20b, bta-miR-29d-3p, bta-miR-1246) specifically targeted progesterone biosynthesis (StAR) and the function of corpus luteum-related genes (CCL11, XCL), suggesting that pregnant cows with elevated progesterone concentrations are more susceptible to heat stress. In addition, we found the differential expression of 11 miRNAs (bta-miR-19a, bta-miR-19b, bta-miR-30a-5p, and several from the bta-miR-2284 family) in both pregnant and non-pregnant cows under heat stress conditions. In target gene prediction and gene set enrichment analysis, these miRNAs were found to be associated with the cytoskeleton, cell junction, vasculogenesis, cell proliferation, ATP synthesis, oxidative stress, and immune responses involved in heat response. These miRNAs can be used as potential biomarkers for heat stress.
