Transcriptome Functional Analysis of Mammary Gland of Cows in Heat Stress and Thermoneutral Condition

SIMPLE SUMMARY: The current study employed RNA-seq technology to analyze the impact of heat stress on the whole transcript sequencing profile in the mammary glands of lactating Holstein dairy cows. In the findings of the current study, heat stress downregulated the expression of casein genes, which resulted in a decrease in milk production. Moreover, heat stress upregulated the gene expression of HSPA1A and HSP90B1, while it downregulated the expression of immune response-related genes that resulted in a reduction in milk yield. Furthermore, there was an increased synthesis of heat shock proteins and unfolded proteins that could reduce the availability of circulating amino acids for milk protein synthesis. The findings of the current experiment may help to explore the impact of heat stress on immune function, milk production, and milk protein synthesis in cows. ABSTRACT: Heat stress (HS) exerts significant effects on the production of dairy animals through impairing health and biological functions. However, the molecular mechanisms related to the effect of HS on dairy cow milk production are still largely unknown. The present study employed an RNA-sequencing approach to explore the molecular mechanisms associated with a decline in milk production by the functional analysis of differentially expressed genes (DEGs) in mammary glands of cows exposed to HS and non-heat-stressed cows. The results of the current study reveal that HS increases the rectal temperature and respiratory rate. Cows under HS result in decreased bodyweight, dry matter intake (DMI), and milk yield. In the current study, a total of 213 genes in experimental cow mammary glands was identified as being differentially expressed by DEGs analysis. Among identified genes, 89 were upregulated, and 124 were downregulated. Gene Ontology functional analysis found that biological processes, such as immune response, chaperone-dependent refolding of protein, and heat shock protein binding activity, were notably affected by HS. The Kyoto Encyclopedia of Genes and Genomes enrichment analysis found that almost all of the top-affected pathways were related to immune response. Under HS, the expression of heat shock protein 90 kDa beta I (HSP90B1) and heat shock 70 kDa protein 1A was upregulated, while the expression of bovine lymphocyte antigen (BoLA) and histocompatibility complex, class II, DRB3 (BoLA-DRB3) was downregulated. We further explored the effects of HS on lactation-related genes and pathways and found that HS significantly downregulated the casein genes. Furthermore, HS increased the expression of phosphorylation of mammalian target of rapamycin, cytosolic arginine sensor for mTORC1 subunit 2 (CASTOR2), and cytosolic arginine sensor for mTORC1 subunit 1 (CASTOR1), but decreased the phosphorylation of Janus kinase-2, a signal transducer and activator of transcription factor-5. Based on the findings of DMI, milk yield, casein gene expression, and the genes and pathways identified by functional annotation analysis, it is concluded that HS adversely affects the immune function of dairy cows. These results will be beneficial to understand the underlying mechanism of reduced milk yield in HS cows.