Journal of Animal Science and Technology
Korean Society of Animal Science and Technology
Article

Integrative multi-omics analysis reveals systemic and intestinal responses to heat stress in finishing pigs

Min-Ki Seok1, Chiwoong Lim1, Young-Jun Seo1, Ji-Yeong Lee1, Hyunjin Kyoung2, Minho Song2,*, Jun-Mo Kim1,**
1Department of Animal Science and Technology, Chung-Ang University, Anseong 17546, Korea.
2Department of Animal Science and Biotechnology, Chungnam National University, Daejeon 34134, Korea.
**Corresponding Author: Minho Song, Department of Animal Science and Biotechnology, Chungnam National University, Daejeon 34134, Korea, Republic of. Phone: +82-42-821-7857. E-mail: mhosng@cnu.ac.kr.
**Corresponding Author: Jun-Mo Kim, Department of Animal Science and Technology, Chung-Ang University, Anseong 17546, Korea, Republic of. Phone: +82-31-670-3263. E-mail: junmokim@cau.ac.kr.

© Copyright 2025 Korean Society of Animal Science and Technology. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Sep 25, 2025; Revised: Oct 24, 2025; Accepted: Nov 12, 2025

Published Online: Nov 17, 2025

Abstract

Heat stress (HS) is a major environmental threat to swine production that impairs growth performance and health. Because of their limited thermoregulatory capacity, pigs are highly vulnerable to HS, which results in compromised intestinal integrity, systemic inflammation, and metabolic inefficiency. To elucidate the mechanisms underlying HS acclimation in pigs, we conducted a longitudinal multi-omics analysis integrating fecal microbiome, whole-blood transcriptome, and immune cell deconvolution in finishing pigs under thermoneutral (TN) or HS conditions. HS markedly reduced the average daily gain and feed intake. Microbiome profiling revealed condition-specific shifts: TN pigs showed enrichment of short-chain fatty acid (SCFA)-producing genera, such as <italic>Prevotella</italic> and <italic>Streptococcus</italic>, whereas HS pigs exhibited increased <italic>Clostridium sensu stricto 1</italic>. Functional predictions indicated preservation of the antioxidant and immunomodulatory pathways (glutathione, retinol, and aminoacyl-tRNA biosynthesis) in TN pigs, whereas HS pigs displayed branched-chain amino acid catabolism, reflecting metabolic acclimation under stress. Transcriptomic analysis revealed acute changes at week 1, with 516 differentially expressed genes enriched in hematopoiesis, focal adhesion, cytoskeletal remodeling, and thyroid hormone signaling. By Week 2, these gene responses had declined, suggesting partial acclimation. Network analysis identified cytoskeletal genes (<italic>ACTB</italic>, <italic>MYL9</italic>, <italic>ACTN4</italic>, and <italic>COL4A4</italic>) as the central regulators. Immune deconvolution further showed the HS-driven elevation of cytotoxic T and myeloid subsets, in contrast to B cell populations which were maintained under TN, highlighting divergent immune trajectories. Integration of the microbiome, transcriptome, and immune data revealed two axes: (1) cytotoxic T cells positively associated with <italic>Clostridium sensu stricto 1,</italic> but negatively associated with cytoskeletal genes, and (2) B cells positively linked to <italic>Prevotella</italic>, <italic>Lactobacillus</italic>, and structural genes. Only the B-cell structural axis formed a coherent cross-layer module, indicatin a recovery-oriented response. These findings demonstrate that resilience to HS requires the coordination of humoral immunity, cytoskeletal reinforcement, and SCFA-producing microbiota. The identified biomarker axis (<italic>Prevotella</italic>, B cells, and cytoskeletal genes) provides a mechanistic basis for developing precise strategies to enhance thermal tolerance in swine.

Keywords: Heat stress; Finishing pigs; Multi-omics integration; Cell deconvolution