# Single‐Nucleus Transcriptome Reveals Cellular Heterogeneity and Transcriptional Response to Heat Stress in Skeletal Muscle

**Authors:** Ziyin Han, Li Chen, Zhiyu Lei, Jiaman Zhang, Ziyu Chen, Xiaolan Fan, Bo Zeng, Anan Jiang, Hai Xiang, Hua Li, Mingzhou Li, Long Jin

PMC · DOI: 10.1002/jcsm.70217 · 2026-02-12

## TL;DR

This study uses single-nucleus transcriptomics to explore how skeletal muscle cells respond to heat stress and recover, revealing cell-type-specific molecular changes.

## Contribution

The study provides the first single-cell transcriptional analysis of heat stress-induced muscle injury and repair in mice.

## Key findings

- Heat stress causes myofibrillar deformation, mitochondrial swelling, and autophagy in skeletal muscle.
- Type IIb myonuclei are the most heat-sensitive subtype, showing significant transcriptional changes after heat exposure.
- Muscle stem cells (MuSCs) express development-related genes like Atp2a1 and Myh1 to aid in muscle repair.

## Abstract

Heat stress can induce skeletal muscle injury. Typical characteristics of heat‐exposed muscle tissues include apoptosis, oxidative stress and autophagy. The understanding of molecular mechanisms underlying heat stress‐induced muscle injury is limited, especially at the single‐cell transcription level.

We collected skeletal muscles from 12‐week‐old female C57BL/6J mice in control (NC) and four heat stress groups. The experimental scheme comprised five groups: NC (25.5°C ± 0.5°C), HS0 (after ~3‐h heat exposure, 41.5°C ± 0.5°C), HS8, HS16 and HS24 group (recovery at 25.5°C ± 0.5°C for 8, 16 and 24 h, respectively). Skeletal muscles were subjected to HE staining (n = 6), TUNEL staining (n = 3) and transmission electron microscopy (n = 3). Transcripts were measured at the tissue (n = 5 or 6) and single‐nucleus levels (n = 2).

Histologically, myofibrillar structure deformation, mitochondrial swelling and fusion, intramuscular triglyceride accumulation and autophagy occurred in the muscles subjected to heat stress. At the tissue level, a gene cluster associated with the response to heat exhibited an increasing trend of transcription in the HS0 versus NC groups, and the levels decreased to those in the NC group after 16 h. At the single‐cell level, 134 320 high‐quality myonuclei were collected from the muscles and annotated as seven cell types, including myonuclei, muscle stem cells (MuSCs) and immune cells. We identified the larger number of differentially expressed genes in the myonuclei. After heat stress, new cell clusters appeared in type IIa/IIx (HS8 group) and IIb (HS0 group) myonuclei but not in type I myonuclei. Generally, immediate early genes, highly expressed genes and transcription factor regulons identified in new cell clusters induced by heat were related to responses to heat, heat shock, oxidative stress and antioxidative stress. To repair the injured muscles, MuSCs highly expressed the development‐related genes, such as Atp2a1, Ckm, Myh1, Aldoa, Pde4d and Pdlim5. Analysis of cell–cell communication showed that Dag1 and Egf signals associated with myonuclei were the underlying pathways that participated in repair tissues.

We constructed the largest transcriptomic dataset, to date, for heat‐exposed skeletal muscles. At tissue resolution, the response of muscles to heat stress was eliminated after a recovery period of 16 h. At single‐nucleus resolution, the myofibre was the most heat‐sensitive cell type, and type IIb myonuclei were the most heat‐sensitive subtype of myonuclei. To our knowledge, this is the first study to reveal the molecular mechanisms underlying heat‐induced muscle injury and repair at the single‐cell transcriptional level.

## Linked entities

- **Genes:** ATP2A1 (ATPase sarcoplasmic/endoplasmic reticulum Ca2+ transporting 1) [NCBI Gene 487], CKM (creatine kinase, M-type) [NCBI Gene 1158], MYH1 (myosin heavy chain 1) [NCBI Gene 4619], ALDOA (aldolase, fructose-bisphosphate A) [NCBI Gene 226], PDE4D (phosphodiesterase 4D) [NCBI Gene 5144], PDLIM5 (PDZ and LIM domain 5) [NCBI Gene 10611], DAG1 (dystroglycan 1) [NCBI Gene 1605], EGF (epidermal growth factor) [NCBI Gene 1950]

## Full-text entities

- **Genes:** Aldoa (aldolase A, fructose-bisphosphate) [NCBI Gene 11674] {aka Aldo-1, Aldo1}, Dag1 (dystroglycan 1) [NCBI Gene 13138] {aka D9Wsu13e, DG, Dp427, Dp71}, Pdlim5 (PDZ and LIM domain 5) [NCBI Gene 56376] {aka 1110001A05Rik, Enh, Enh1, Enh2, Enh3, LIM}, Ckm (creatine kinase, muscle) [NCBI Gene 12715] {aka CPK-M, Ckmm, M-CK, MCK}, Atp2a1 (ATPase, Ca++ transporting, cardiac muscle, fast twitch 1) [NCBI Gene 11937] {aka SERCA1}, Myh1 (myosin, heavy polypeptide 1, skeletal muscle, adult) [NCBI Gene 17879] {aka A530084A17Rik, IId, IId/x, MHC-2X/D, MHC2X/D, MYHC-IIX}, Egf (epidermal growth factor) [NCBI Gene 13645], Pde4d (phosphodiesterase 4D, cAMP specific) [NCBI Gene 238871] {aka 9630011N22Rik, Dpde3}
- **Diseases:** muscle injury (MESH:D009135)
- **Chemicals:** triglyceride (MESH:D014280)
- **Species:** Mus musculus (house mouse, species) [taxon 10090]

## Figures

17 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12895210/full.md

---
Source: https://tomesphere.com/paper/PMC12895210