Single-nucleus RNA sequencing reveals cell type-specific responses to heat stress in bovine mammary gland

Abstract Background Heat stress (HS) poses a major challenge to the dairy industry by reducing milk production, yet its cell type-specific effects in the bovine mammary gland remain incompletely defined. In this study, we recorded production traits and collected mammary biopsies from cows under thermoneutral (TN), HS, and pair-fed (PF) conditions. Results Clinical measurements confirmed HS-induced physiological alterations. Compared with TN cows, HS cows exhibited reduced dry matter intake (DMI), milk yield, and yields of fat, protein, and lactose, along with increased water intake and milk urea nitrogen. The use of PF controls indicated that decreased DMI accounted for 45% of the milk-yield reduction, whereas direct HS effects accounted for the remaining reduction. We applied single-nucleus RNA-seq (snRNA-seq) on mammary biopsies to generate cell-resolved HS responses. We identified 14 distinct cell clusters, including epithelial, immune, and stromal populations. Under the TN condition, casein genes (e.g., CSN1S1 , CSN2 ) were broadly expressed across luminal cells but were attenuated under HS, whereas luminal alveolar cells showed relative upregulation. Heat shock protein genes were strongly induced by HS, primarily in epithelial clusters. Gene-set enrichment analyses revealed increased ribosomal activities across HS-responsive clusters and enrichment of protein folding and metabolic pathways in luminal alveolar cells, suggesting elevated proteostasis demands under stress. Pseudotime analysis positioned luminal cells along a progenitor-to-secretory trajectory under TN, accompanied by increased casein gene expression, whereas under HS, mature luminal cells shifted toward a homeostasis regulatory state. Cell-cell communication analysis demonstrated HS-induced remodeling of interepithelial signaling, including altered ERBB4-mediated signaling from luminal hormone-sensing to alveolar lineages. Finally, transcription factor activity profiling highlighted cell type-specific HS-activated regulators and their downstream target genes. Conclusions Together, this cell type-resolved atlas delineates how HS alters bovine mammary epithelial function, developmental state, and intercellular crosstalk. These findings point to proteostasis pressure, disrupted signaling pathways, and rewired regulatory networks as mechanistic contributors to reduced lactational performance under HS, offering insights for improving heat resilience in dairy cattle.