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

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.