Preprint
Article

Transcriptome Analysis Revealed Potential Genes of Skeletal Muscle Thermogenesis in Mashen Pigs and Large White Pigs under Cold Stress

Altmetrics

Downloads

77

Views

34

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

09 September 2023

Posted:

12 September 2023

You are already at the latest version

Alerts
Abstract
Pigs are susceptible to cold stress due to the absence of brown fat caused by the partial deletion of uncoupling protein 1 during their evolution. Some local pig breeds in China, exhibit potential cold adaptability, but research has primarily focused on fat and intestinal tissues. Skeletal muscle plays a key role in adaptive thermogenesis in mammals, yet the molecular mechanism of cold adaptation in porcine skeletal muscle remains poorly understood. This study investigated the cold adaptability of two pig breeds, Mashen pigs (MS) and Large White pigs (LW), in a 4-day cold (4°C) or room temperature (25°C) environment. We recorded phenotypic changes and collected blood, and longissimus dorsi muscle for transcriptome sequencing. Finally, a candidate gene was randomly selected for functional exploration in porcine skeletal muscle satellite cells. A decrease in body temperature and body weight in both LW and MS pigs under cold stress, accompanied by increased shivering frequency and respiratory frequency. However, MS pigs demonstrated stable physiological homeostasis, indicating a certain level of cold adaptability. LW pigs primarily responded to cold stress by regulating heat production and Glycolipid energy metabolism. MS pigs exhibited a distinct response to cold stress, involving regulation of heat production, energy metabolism pathways, robust fatty acid oxidation ability, and stronger immune response. Furthermore, the functional exploration of PRSS8 in porcine skeletal muscle satellite cells revealed it affected cellular energy metabolism and thermogenesis by regulating ERK phosphorylation. These findings shed light on the diverse transcriptional responses of skeletal muscle in LW and MS pigs under cold stress, offering valuable insights into the molecular mechanisms underlying cold adaptation in pigs.
Keywords: 
Subject: Biology and Life Sciences  -   Animal Science, Veterinary Science and Zoology

1. Introduction

Cold stress is a significant environmental factor in alpine regions, posing various adverse effects on the physiological systems of animals, including the nervous, endocrine, immune, and antioxidant systems. This situation leads to substantial threats to animal husbandry production and development [1,2]. Pigs, in particular, are highly susceptible to ambient temperature variations, with cold stress contributing to an increase in the incidence of pig diseases [3] and mortality rates [4], as well as diminishing meat quality [5]. These challenges ultimately lead to severe economic losses. The absence of uncoupling protein 1 (UCP1) in pigs throughout their evolutionary history has resulted in a lack of brown adipose tissue (BAT), rendering them more cold-sensitive. However, certain pig breeds in China, such as Tibetan pigs [6] and Min pigs [7], have demonstrated innate cold adaptability due to their geographical exposure. Genetic variations among pig breeds influence their cold tolerance characteristics [8]; nonetheless, the precise molecular mechanisms of cold tolerance in pigs remain elusive. Notably, when Tibetan pigs were acutely exposed to 4°C for 4 hours, they showed heightened expression of uncoupling protein 3 (UCP3) and a beige-like appearance of white fat, suggesting the crucial role of fat metabolism in cold adaptation [9]. Furthermore, Yorkshire pigs under chronic cold stress exhibited visible pathological changes in jejunal and ileal mucosa, with enhanced expression of Toll-like receptor 4 (TLR4), myeloid differentiation primary response gene 88 (MyD88), and nod-like receptor pyrin domain 3 (NLRP3). In contrast, Min pigs did not exhibit such harmful consequences under similar conditions [10].
The Mashen (MS) pig, a wild boar species found in northern China, is renowned for its strong resistance to stress, high-quality meat, high fecundity, slow growth rate, and ability to thrive on roughage [11,12]. Recent studies have also found that MS pigs exhibit certain cold adaptability. When mice were transplanted with MS pig microbiota and subjected to 21 days of cold exposure at 4°C, they showed more stable body temperature and a better-preserved intestinal structure compared to mice transplanted with Duroc-Landrace-Yorkshire (DLY) pig microbiota. Additionally, the expression of UCP1, carnitine palmitoyltransferase 1B (CPT1B), and peroxisome proliferator-activated receptor-gamma coactivator 1 alpha (PGC-1α) significantly increased in the MS pig microbiota-transplanted mice (P<0.05) [13]. Moreover, in low-temperature environments, local MS pigs demonstrated a stronger ability to maintain intestinal physiological function compared to Large White (LW) and Jinfen White pigs. Their duodenal esophageal wall exhibited higher trypsin activity, promoting an inflammatory response [14]. Current research on cold adaptation in pig chiefly centers on fat and intestine, but further exploration is required to understand the cold adaptation mechanism in other organs.
Skeletal muscle accounts for approximately 40% of total body weight in animals [15] and serves as the largest metabolic organ [16]. In response to cold stress, shivering thermogenesis (ST) in skeletal muscle is the primary means of thermogenesis, involving spontaneous repeated muscle contractions and relaxations [17,18]. This process relies mainly on actin, myosin, and adenosine triphosphate (ATP) hydrolysis [19]. Additionally, non-shivering thermogenesis (NST) in skeletal muscle is a major contributor to cold adaptation, and involves ineffective Ca2+ cycling of sarcoplasmic reticulum Ca2+-ATPase (SERCA) [20] and uncoupling of the mitochondrial inner membrane [21]. Studies on Lasiopodomys brandtii had shown that removal of interscapular BAT resulted in enhanced thermogenesis in skeletal muscle under a 4°C environment for 28 days. This enhancement was characterized by amplified mitochondrial quantity and metabolic activity, along with elevated levels of sarcolipin (SLN) and sarcoplasmic reticulum Ca2+-dependent ATP [22]. Furthermore, newborn wild boars exposed to cold environments during their first 5 days after birth exhibited increased contribution of NST, accompanied by age-related increases in body temperature, SERCA activity, and expression of SERCA1a and SLN in skeletal muscle tissue [23]. In a study involving zinc-α2-glycoprotein knockout (ZAG-KO) mice and wild-type (WT) mice exposed to low temperature (6°C) for one week, it was found that cold stress increased the expression levels of lipolysis-related proteins (ATGL and p-HSL) and energy metabolism-related proteins (PGC-1α, UCP2, UCP3, and COX1) in the gastrocnemius muscle of WT mice [24]. Meanwhile, another study discovered that UCP3 overexpression in mouse skeletal muscle resulted in reduced efficiency of muscle mitochondria, leading to a 42% reduction in the ratio of ATP synthesis and mitochondrial oxidation [25]. This indicated that UCP3 is subject to complex regulatory mechanisms in NST. To support the increase in animal heat production, energy metabolism, especially glucose and lipid metabolism, needs improvement [26,27]. ST primarily involves glycogen decomposition in skeletal muscle, with the rate of glycogen utilization correlating with muscle contraction strength [28,29]. Through the screening of differentially expressed genes (DEGs) from rat iliopsoas muscles exposed to ultra-low temperatures, we observed their involvement in muscle regeneration, tissue repair, and lipid metabolism [30]. Based on these findings, we infer that pigs respond to cold stress by upregulating glycolipid energy metabolism, which provides sufficient energy for skeletal muscle thermogenesis.
In this study, local breed MS pigs and imported breed LW pigs were used as experimental subjects exposed to either low temperature (4℃) or normal temperature (25℃) conditions for a period of 4 days. Blood biochemical parameters were measured using enzyme-linked immunosorbent assay (ELISA), and histological analysis, including ultrastructure and enzyme activity assessment, of longissimus dorsi muscle was performed using transmission electron microscopy and enzyme histochemical staining. Furthermore, RNA-seq was performed to evaluate the gene expression profiles in the longissimus dorsi muscles of both pig breeds. The DEGs and enriched metabolic pathways of LW and MS pigs in response to cold stress were screened and analyzed using Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG), and weighted gene co-expression network analysis (WGCNA). The findings were verified through quantitative real-time polymerase chain reaction (qRT-PCR) and Western blotting (WB). Additionally, a candidate gene responsive to cold stress was randomly selected for functional verification and preliminary exploration of the mechanism in porcine skeletal muscle satellite cells. This study will provide a comprehensive analysis of the distinct transcriptional responses between MS and LW pigs under cold stress, with a specific focus on the functionality of the longissimus dorsi muscle. The outcomes will provide valuable insights into the cold adaptation mechanism of different pig breeds and contribute to future research on breeding cold-resistant pigs.

2. Results

2.1. Different Hysiological Response Patterns of Large White Pigs and Mashen Pigs under Cold Stress

LW and MS pigs were subjected to normal temperature (25℃) and low temperature (4℃) for a duration of 4 days. During the experiment, the body temperature of both pig breeds at 4℃ was significantly lower (P<0.05) than at 25℃ (Figure 1A). However, the body temperature deviation of MS pigs was most significantly (P<0.01) or significantly (P<0.05) lower than that of LW pigs under cold stress, with the exception of the third day (Figure 1B). No significant difference in respiratory rate was observed between the two breeds at either 25℃ or 4℃ (Figure 1C). During the cold treatment, the weight of MS pigs decreased by 2.98 kg, representing a 7.2% reduction compared to their initial weight (Supplementary Table S3). The weight of LW pigs decreased by 1.29 kg, representing a 3.6% reduction compared to their initial weight. At 4℃, the shivering frequency of MS pigs was significantly higher than that of LW pigs on the first day (P<0.05) and significantly lower (P<0.05) than that of LW pigs on the fourth day (Figure 1D). Interestingly, during the cold treatment, the shivering frequency of MS pigs recovered more gently compared to LW pigs, indicating that MS pigs maintained more stable physiological homeostasis and exhibited better adaptability to low temperatures.

2.2. Serum Biochemistry and Hormone Analyses of Large White Pigs and Mashen Pigs under Cold Stress

At a temperature of 4℃, the activities of three thermogenic enzymes, four glucose metabolism enzymes, and two lipid metabolism enzymes were measured. It was found that EPI, NE, and Ca2+/Mg2+-ATPase levels in LW pigs increased most significantly (P<0.01) or significantly (P<0.05) under cold stress. In addition, in LW pigs, there was a notable increase in the glucose metabolism enzymes LDH, PFK-1, HK2, and M2-PK, as well as the lipid metabolism enzymes SDH and CS (P<0.01 or P<0.05) (Figure 2A-I). However, the thermogenic and glucose metabolism enzymes in MS pigs did not show significant changes under cold stress, except for SDH and CS, which exhibited a marked increase (P<0.01). At 4℃, the SDH and CS levels in MS pigs were significantly higher (P<0.01) than those in LW pigs. These findings suggest that LW pigs respond to cold stress by increasing heat production, glucose metabolism, and lipid metabolism, whereas MS pigs primarily respond to cold stress by enhancing lipid metabolism, with higher levels of fatty acid oxidation compared to LW pigs.

2.3. Effects of Cold Stress on Ultrastructure and Enzyme Activity of Longissimus Dorsi Muscle in Large White Pigs and Mashen Pigs under Cold Stress

The electron microscope analysis of pig longissimus dorsi revealed changes in the abundance and ultrastructure of mitochondria in both LW and MS pigs under cold stress. At 4°C, the mitochondrial abundance between skeletal muscle fibers was higher in both LW and MS pigs compared to 25°C (Figure 3A). However, the cristae structure of mitochondria in MS pigs appeared more intact, while the cristae structure in LW pigs showed damage under cold stress. Further, MS pigs exhibited a higher number of mitochondria between muscle fibers than LW pigs at both 25℃ and 4℃. SDH activity staining (Figure 3B) confirmed a significant increase (P<0.01) in fatty acid oxidation metabolism in both LW and MS pigs under cold stress, with MS pigs showing the highest increase (P<0.01) compared to LW pigs. LDH activity staining (Figure 3C) revealed that LW pigs had a significant increase (P<0.01) in glycolysis metabolism under cold stress, while MS pigs did not exhibit a significant change in glycolysis metabolism under cold stress.

2.4. Differential Expression Genes Analysis of Longissimus Dorsi Muscle in Large White Pigs and Mashen Pigs under Cold Stress

RNA-seq analysis was conducted on the longissimus dorsi muscle samples from 12 pigs. The results indicated that the Q30 values of all filtered samples were above 90% (Supplementary Table S4). The distribution of sequencing data was even among the 12 subjects, with no significant difference observed between different individuals (Figure 4A). The mRNA heat map (Figure 4B and Supplementary Figure S1A-F) indicated that the three biological replicates of each group were of high quality, and the samples clustered well within the four groups: normal-temperature group (LW25-vs-MS25), low-temperature group (LW4-vs-MS4), LW pig group (LW25-vs-LW4), and MS pig group (MS25-vs-MS4). In the LW pig group, a total of 750 genes were identified, with 374 up-regulated genes and 376 down-regulated genes (Figure 4C). Conversely, in the MS pig group, only 46 genes were identified, comprising 31 up-regulated genes and 15 down-regulated genes. Notably, the proportion of up-regulated and down-regulated DEGs in the LW pig group was roughly balanced, while the number of up-regulated DEGs in the MS pig group was twice that of down-regulated DEGs (Figure 4D). Within the overlapping region of the LW pig group and the MS pig group, 23 DEGs played a crucial role in the common cold response shared by both LW and MS pigs. On the left side of the non-coincidence region, 727 DEGs exerted regulatory roles in the specific cold response of LW pigs, while on the right side, the 23 DEGs governed the specific cold response of MS pigs (Figure 4E; Supplementary Figure S1G).
The fold change of DEGs in the LW pig group (LW-25-vs-LW-4) and MS pig group (MS-25-vs-MS-4) mostly ranged from 1 to 2 times, followed by 2 to 3 times, and the least was more than 3 times (Supplementary Figure S1H-I). Among the 23 common DEGs in response to cold stress shared between the LW pig group and MS pig group (Figure 4F), 19 DEGs were significantly up-regulated, and 4 DEGs were significantly down-regulated. Notably, NR4A2 and LEP were the genes with the largest up-regulation and down-regulation folds, respectively, in both the LW pig group and MS pig group (Supplementary Table S5). NR4A2 was up-regulated by 2.84 times in the LW pig group and 2.61 times in the MS pig group, while LEP was down-regulated by 3.83 times in the LW pig group and 3.55 times in the MS pig group. Additionally, there were 23 unique cold-responsive genes in the MS pig group (Figure 4G). Among these, LYZ and RDH16 exhibited the largest up-regulation and down-regulation folds, respectively, in the MS pig group, with LYZ being up-regulated by 3.40 times and RDH16 being down-regulated by 3.60 times (Supplementary Table S6). Furthermore, the LW pig group had 727 unique cold-responsive genes, with CR1 and HGFAC showing the largest up-regulation and down-regulation folds, respectively, in this group. CR1 was up-regulated by 6.13 times, while HGFAC was down-regulated by 9.74 times in the LW pig group (Supplementary Table S7).

2.5. Functional Annotation and Pathway Enrichment Analysis of Common and Specific Cold-Responsive Differentially Expressed Genes in Large White Pigs and Mashen Pigs under Cold Stress

The DEGs shared by both LW and MS pigs in response to cold stress were enriched in GO terms related to heat production and energy metabolism (Figure 4H). Specifically, ARRDC3, ADRB2, CEBPB, and NR1D1 were enriched in cold-induced thermogenesis (Supplementary Table S8), and were also associated with various pathways related to glycolipid energy metabolism. SIK1 was found to be enriched in negative regulation of metabolic processes, lipid metabolic processes, and carbohydrate biosynthetic processes. Additionally, NR4A2, CIART, KLF10, and TIPARP were enriched in negative regulation of metabolic process, with NR4A2 being the most up-regulated gene in both LW and MS pigs under cold stress. Furthermore, GPCPD1 and SCD were enriched in lipid metabolic processes. The KEGG pathways shared by LW and MS pigs in response to cold stress were significantly enriched in environmental information processing, metabolism, human diseases, organismal systems, and cellular processes. These pathways included adenosine monophosphate-activated protein kinase (AMPK) signaling pathway, non-alcoholic fatty liver disease, JAK-STAT signaling pathway, biosynthesis of unsaturated fatty acids, circadian rhythm, and apoptosis-multiple species (Figure 4J). For instance, SCD was enriched in the AMPK signaling pathway and biosynthesis of unsaturated fatty acids, and it was also associated with the lipid metabolism pathway in GO. NR1D1 was enriched in circadian rhythm, which was related to environmental adaptation, and it was also involved in the heat production and glycolipid energy metabolism pathways in GO. Moreover, LEP was enriched in multiple energy metabolism pathways, such as the AMPK signaling pathway, non-alcoholic fatty liver disease, and JAK-STAT signaling pathway. Notably, LEP was the most down-regulated gene in both LW and MS pigs under cold stress, and it was also enriched in the lipid metabolism pathway in GO. Based on the results of GO and KEGG analyses, the candidate genes of LW and MS pigs in response to cold stress were identified as ARRDC3, ADRB2, CEBPB, SIK1, NR4A2, GPCPD1, CIART, KLF10, TIPARP, SCD, NR1D1, and LEP. As such, these genes appear to play vital roles in regulating heat production and energy metabolism pathways during cold stress adaptation in LW and MS pigs.
The DEGs specific to the response of MS pigs to cold stress were significantly enriched in various pathways related to cold adaptation, including hemoglobin metabolism, adaptive immunity, fiber contraction, oxygen transport, and transport regulation (Figure 4I). Among these pathways, specific genes were identified, ALAS2 was enriched in hemoglobin metabolic process (Supplementary Table S9). DUSP10 and IL18 were enriched in the adaptive immune response, ACTC1 and ANKRD1 were enriched in the contractile fiber part, and ANKRD1 was also enriched in the regulation of transport. Further, HBB was enriched in oxygen transport, while PRSS8, NEFH, and KCNC4 were enriched in the regulation of transport. Additionally, the specific KEGG pathway for MS pigs in response to cold stress was significantly enriched in human diseases and metabolism, including pathways related to immunity and metabolism such as African trypanosomiasis, malaria, legionellosis, influenza A, and porphyrin metabolism (Figure 4K). HBB and IL18 were enriched in African trypanosomiasis and malaria, respectively, and they were also involved in oxygen transport and the adaptive immune response as identified in the GO analysis. HSP70.2 was enriched in legionellosis and influenza A, and it had been studied extensively for its role in cold hardiness as a member of the heat shock protein family. ALAS2 was enriched in porphyrin metabolism and was also identified in the GO analysis for its involvement in the hemoglobin metabolic process. These findings suggest that the specific response of MS pigs to cold stress is mainly achieved by participating in immune and metabolic pathways. Based on the GO and KEGG results, candidate genes specific to the cold stress response of MS pigs were identified as ACTC1, PRSS8, DUSP10, ANKRD1, NEFH, KCNC4, HBB, IL18, HSP70.2, and ALAS2.
The DEGs specific to the response of LW pigs to cold stress were significantly enriched in pathways related to cold environment adaptation, including thermogenesis, glycolysis, and fatty acid oxidation (Supplementary Figure S1J). Notably, ACSL1, CD36, ELOVL6, and PPARGC1A were enriched in the cold-induced thermogenesis pathway (Supplementary Table S10). Additionally, CPT2, IRS2, ACADM, HADHB, PCK2, and APOC3 were involved in the cellular lipid catabolic process, with IRS2 and ACADM also being associated with glucose metabolism. Furthermore, IGF1 and HK2 were linked to the regulation of glucose import. The specific KEGG pathways of LW pigs in response to cold stress included the forkhead box O (FoxO) signaling pathway, peroxisome proliferator-activated receptor (PPAR) signaling pathway, fatty acid metabolism, AMPK signaling pathway, and insulin resistance, which were closely related to cold stress (Supplementary Figure S1K). Among these pathways, PPARGC1A, CD36, IRS2, PCK2, and IGF1 were enriched in the AMPK signaling pathway, and PCK2 was additionally enriched in the adipocytokine signaling pathway and insulin resistance. ACSL1, CPT2, ACADM, and APOC3 were enriched in the PPAR signaling pathway, while ELOVL6 and HADHB were related to fatty acid metabolism. Moreover, IRS2 and HK2 were associated with Type II diabetes mellitus. These findings indicate that LW pigs specifically respond to cold stress by actively participating in heat production and glycolipid metabolism pathways. Based on the GO and KEGG analysis results, the candidate genes specifically responding to cold stress in LW pigs were inferred to be PPARGC1A, CD36, IRS2, PCK2, IGF1, ACSL1, CPT2, ACADM, APOC3, IRS2, and HK2.

2.6. Weighted Co-Expression Network Analysis of Longissimus Dorsi Muscle in Large White Pigs and Mashen Pigs under Cold Stress

Using the WGCNA method, we analyzed a total of 12,365 genes from the transcriptome data after quality control. The optimal soft threshold power of 17 was determined with a criterion of R2>0.8 to construct a scale-free network (Supplementary Figure S2A). Based on the gene expression patterns, we divided the modules into clusters and merged those with a similarity greater than 0.75. Ultimately, we identified a total of 16 modules, and genes with no clear co-expression trend were placed in the grey module (Supplementary Figure S2B and Figure 5A,B). Among these modules, the blue2 module contained the highest number of genes, with 3,766 genes, while the grey module had the fewest, with 60 genes.
We calculated the correlation coefficients and P values between the modules and various physiological traits, including body weight, body temperature, respiratory rate, shivering frequency, EPI, NE, Ca2+/Mg2+-ATPase, LDH, PFK-1, HK2, SDH, CS, and M2-PK, using the eigenvector values of each module (Figure 5C). The modules showing a strong correlation with the target traits (R>0.8, P<0.01) were selected as the target modules representing the common response to cold stress in both pig breeds. Our analysis identified two such modules, namely the pink and yellow modules (Supplementary Table S11), which showed a high correlation with multiple target traits. Further examination of the gene expression pattern in the heat map module revealed that the pink module exhibited a more similar expression trend between LW and MS pigs under cold stress, suggesting a better representation of the common response to cold stress in both breeds (Figure 5D). Additionally, we observed that the blue2 module had the highest correlation with the MS-4 group, indicating its relevance to the unique response mode of the MS-4 group to cold stress. Consequently, we chose the pink and blue2 modules as the target modules for subsequent functional annotation and enrichment analysis.

2.7. Functional Annotation and Pathway Enrichment Analysis of the Target Modules

GO analysis of the pink module, which reflected the common response to cold stress in LW and MS pigs, revealed significant enrichment in various biological processes, including catabolic process, cellular metabolic process, generation of precursor metabolites and energy, oxidation-reduction process, energy derivation by oxidation of organic compounds, and macroautophagy (Figure 6A-B). Within this module, ARRDC3, ADRB2, CEBPB, NR1D1, UCP3, and CD36 were enriched in adaptive thermogenesis, while CASQ1 and ACADM were enriched in response to temperature stimulus (Supplementary Table S12). Additionally, SIK1 and PDK2 were enriched in the regulation of gluconeogenesis, and RCAN1, PPP3CC, and ATP1B1 were enriched in calcium-mediated signaling. In terms of the KEGG analysis, the pink module was mainly enriched in pathways such as the citrate cycle (TCA), Epstein-Barr virus infection, oxidative phosphorylation, carbon metabolism, non-alcoholic fatty liver disease, AMPK signaling pathway, thermogenesis, and pyruvate metabolism (Figure 6C). Specifically, CD36 and ACADM were enriched in the PPAR signaling pathway (Supplementary Table S13), while RCAN1 and PPP3CC were enriched in Kaposi sarcoma-associated herpesvirus infection, and ADRB2 and ATP1B1 were enriched in the cGMP-PKG signaling pathway. Based on the combined results of GO and KEGG analyses from the WGCNA, the candidate target genes of the pink module were identified as ARRDC3, ADRB2, CEBPB, NR1D1, UCP3, CD36, CASQ1, ACADM, SIK1, PDK2, RCAN1, PPP3CC, and ATP1B1.
GO analysis of the blue2 module, which reflected the unique response of MS pigs to cold stress, revealed significant enrichment in various biological processes, including cellular metabolic process, RNA processing, primary metabolic process, and nucleic acid metabolic process (Figure 6D-E). Within this module, ACTC1, PRSS8, GPCPD1, NR4A2, SIRT6, PSMD14, PSMC6, NCBP2, and PYM1 were enriched in the metabolic process (Supplementary Table S14), while IL18 was significantly enriched in the apoptotic signaling pathway. Regarding the KEGG analysis, the blue2 module was mainly enriched in pathways related to heat production, energy metabolism, and viral infection, including oxidative phosphorylation, thermogenesis, non-alcoholic fatty liver disease, metabolic pathways, and Epstein-Barr virus infection (Figure 6F). Specifically, SIRT6 was enriched in thermogenesis (Supplementary Table S15), while PSMD14 and PSMC6 were enriched in Epstein-Barr virus infection, and NCBP2 and PYM1 were enriched in nucleocytoplasmic transport. Based on the combined results of GO and KEGG analyses from the WGCNA, the candidate target genes of the blue2 module were identified as ACTC1, PRSS8, GPCPD1, NR4A2, SIRT6, PSMD14, PSMC6, NCBP2, PYM1, and IL18.
The pink module exhibited only 7 intersection genes with the cold response gene set shared by both LW and MS pigs (Figure 6G). On the other hand, the blue2 module showed 4 intersection genes with the cold response gene set shared by LW and MS pigs, and 6 intersection genes with the cold response gene set unique to MS pigs. Altogether, there were 17 intersection genes in the pink and blue2 modules, which overlapped with the cold-responsive gene sets shared by LW and MS pigs, as well as the gene set specific to MS pigs. Further analysis established that the 17 intersection genes were Venn with the candidate genes obtained from the GO and KEGG analyses based on the pink module and the cold response gene set shared by LW and MS pigs, respectively. The results indicated that ARRDC3, ADRB2, CEBPB, NR1D1, and SIK1 were the cold response genes shared by both LW and MS pigs, as identified by both methods. Similarly, the 17 intersection genes were Venn with the candidate genes obtained by the GO and KEGG analyses based on the blue2 module and the cold response gene set specific to MS pigs, respectively. The results showed that ACTC1, PRSS8, and IL18 were the specific cold response genes of MS pigs, as determined by both methods.

2.8. Validation Results of Common and Specific Cold Response Candidate Genes in Pig Longissimus Dorsi Muscle and Satellite Cells

The qRT-PCR results of 8 candidate genes in the longissimus dorsi muscle of LW and MS pigs under cold stress evidenced that the expression trends of ARRDC3, ADRB2, CEBPB, NR1D1, SIK1, ACTC1, PRSS8, and IL18 were consistent with the RNA-seq results (Figure 7A). Additionally, qRT-PCR results of thermogenesis and energy metabolism in the longissimus dorsi muscle of LW and MS pigs under cold stress demonstrated a significant increase (P<0.05) in SERCAα2, SLN, CASQ2, UCP3, PGC1-α heat-related genes in both breeds (Supplementary Figure S3). Cold stress was also found to significantly promote (P<0.05) the expression of energy metabolism genes, including PFK, PKM, LDHA, LDHA, HADHB, and CPT1B in LW and MS pigs. WB results showed that except UCP3, CASQ2, HADHB, PFK1, and other heat-related proteins were significantly increased (P<0.05) in LW and MS pigs under cold stress (Figure 7B and Supplementary Figure S4). Of note, the UCP3 expression in MS pigs was most significantly (P<0.01) higher than that in LW pigs at both 25°C and 4°C, suggesting a higher mitochondrial oxidative phosphorylation level in MS pigs.
The qRT-PCR results of 8 candidate genes in porcine skeletal muscle satellite cells under cold stress revealed a significant up-regulation (P<0.05) in the 32°C cold exposure group compared with the 37°C normal temperature group; this finding was consistent with the RNA-seq results (Figure 7C). For cell function exploration, we selected the PRSS8 gene, which specifically responds to cold stress in MS pigs. The interference efficiency of qRT-PCR revealed that sus-si-PRSS8-1, sus-si-PRSS8-2, and sus-si-PRSS8-3 groups significantly decreased (P<0.05) the expression of PRSS8 (Figure 7D), making them suitable for further experiments. The qRT-PCR results demonstrated that interfering with PRSS8 in porcine skeletal muscle satellite cells significantly downregulated (P<0.05) the expression of heat production genes, including SERCAα2, SLN, UCP3, and PGC-1α (Figure 7E). Moreover, the interference with PRSS8 also led to a significant decrease (P<0.05) in the expression of energy metabolism genes, such as GLUT1, GLUT4, PFK, PKM, HADHB, and CPT1B. In addition, WB analysis highlighted that interfering with PRSS8 in porcine skeletal muscle satellite cells significantly reduced (P<0.05) the expression of SERCAα2, CASQ2, UCP3, and PGC-1α heat production proteins (Figure 7F and Supplementary Figure S4). The interference with PRSS8 also led to a significant decrease (P<0.05) in CPT1B protein levels, indicating that fatty acid oxidation was inhibited. Moreover, the level of extracellular signal-regulated kinase (ERK) phosphorylation was significantly decreased (P<0.05) following the interference of PRSS8.

3. Discussion

Animals exhibit various adaptive responses under cold stress, including increased respiratory rate, decreased body temperature, diminished feed intake, and increased heat production in low-temperature environments. Pigs convert most of the nutrients in their diet into heat energy to maintain body temperature and compensate for heat loss caused by the cold; this subsequently results in decreased body weight and feed conversion rate [31,32]. Previous studies have shown that the cold tolerance of Chinese Meishan pigs and European LW pigs is closely associated with body weight and increases significantly with age (P<0.01). However, Meishan piglets have 16% lower body weight than LW piglets, yet both exhibit similar cold tolerance [33]. The findings of this study also support the notion that cold stress leads to weight loss in both LW and MS pigs, though it may not be significant due to the short duration of the cold stress. Additionally, the body weight of MS pigs was lower than that of LW pigs, yet they exhibited better cold adaptability. A study involving piglets (~13 days old) subjected to five days of cold stress observed a decrease in core temperature and an increase in cold adaptability achieved through postural adjustments [34]. Another study on piglets after birth found that tremor intensity gradually decreased under cold stress, while heat production and muscle blood flow gradually increased [35,36]. The initial explanation for the decrease in trembling intensity attributed it to an increase in trembling efficiency [37]. However, current research suggests that the decrease is due to increased NST in muscles [23]. The results of the present study also revealed that cold stress significantly decreased (P<0.05) the body temperature of both LW and MS pigs, in addition to significantly increasing (P<0.05) the shivering frequency. However, on the fourth day of cold stress, the shivering frequency of MS pigs decreased significantly (P<0.05), and their body temperature deviation was small, implying that MS pigs exhibit superior cold adaptability.
Cold stress triggers a series of adaptive responses in animals, involving the regulation of brain hormones to skeletal muscle for specific cold adaptations [38]. NE, the main neurotransmitter, is secreted under cold stress to mediate skin vasoconstriction and reduce heat loss [39]. Cold exposure influences systemic energy metabolism, activating BAT to mediate NST and reduce plasma triglyceride concentration [40]. In piglets, cold stress was found to significantly increase the activities of glycolytic potential (GP) (P<0.001) and LDH (P=0.03), along with the activities of CS and β-hydroxy-acyl-CoA dehydrogenase (HAD), indicating increased glycolysis and oxidative metabolism [41]. The current study also examined serum hormones and enzyme activities, evidencing that cold stress significantly increased (P<0.05) the activities of EPI, NE, Ca2+/Mg2+-ATPase, LDH, PFK-1, HK2, M2-PK, SDH, and CS in LW pigs. Conversely, there were no significant changes in the thermogenic and glucose metabolism enzymes of MS pigs under cold stress, except for a marked increase (P<0.01) in SDH and CS lipid metabolism enzymes. Notably, at 4℃, the SDH and CS lipid metabolism enzymes were statistically significantly higher (P<0.01) than those of LW pigs. Enzyme activity staining of the longissimus dorsi muscle in LW and MS pigs also demonstrated that LW pigs responded to cold stress by enhancing fatty acid oxidation and glycolysis, while MS pigs primarily increased fatty acid oxidation. Electron microscopy of the longissimus dorsi muscle revealed an increase in fatty acid oxidation in response to cold stress by elevating mitochondrial abundance in skeletal muscle. Of note, the skeletal muscle of MS pigs exhibited significantly higher mitochondrial abundance (P<0.05) compared to LW pigs at both 25°C and 4°C; furthermore, the mitochondrial ridge structure of MS pigs under cold stress appeared more intact, indicating their superior ability to carry out fatty acid oxidation ability.
Analyzing the entire mRNA transcriptome of cells, tissues, or organisms has proven to be a valuable approach for identifying specific genes in the response to cold stress. Further, transcriptome sequencing is an increasingly useful tool for uncovering key pathways and biological processes associated with how animals adapt to changes in body temperature [38]. Earlier research has highlighted the increased transcriptional plasticity observed in deer mice [42] and chickens [43] in response to cold stress. In the current study, it was found that the number of DEGs and enrichment pathways in MS pigs under cold stress were notably lower compared to LW pigs, suggesting that MS pigs were less affected by cold stress and exhibited a certain level of cold adaptability. Transcriptome analysis of Songpu mirror carp and bar carp adapted to 28°C and 18°C revealed distinct cold tolerance between the larvae of the two breeds. The study identified several biological processes, including morphogenesis, secondary alcohol metabolism, drug transport, and the FoxO signaling pathway, which were closely associated with the development of cold tolerance and adaptability in these different breeds [44]. The transcriptional characteristics of inguinal fat in mice under cold stress showed that the inguinal white fat was primarily engaged in fatty acid elongation, three acylglycerols, sheath fat, and triglyceride synthesis pathways following a brief cold exposure at 4°C for 3 days [45]. The transcriptome sequencing analysis of modern tropical and subtropical elephants highlighted the presence of DEGs that generated single nucleotide polymorphisms, fragment insertions, deletions, and amplifications of intact gene regions. These genes were found to be involved in lipid metabolism and thermogenesis, aiding elephants in coping with harsh and cold environments [46]. In a transcriptome study of grizzly bears during hibernation, the expression of fatty acid synthase and acetyl-CoA carboxylase, which were rate-limiting enzymes in long-chain fatty acid synthesis, decreased in the fat and muscle tissues of bears [47]. However, there was an increased expression of lipoprotein lipase, which acts as a negative regulator of hormone-sensitive lipase, triglyceride lipase, and monosaccharide lipase. This increase indicated that lipolysis was indirectly activated during hibernation. In the present study, the DEGs of LW and MS pigs in response to cold stress were found to be involved in various biological processes, including cold-induced thermogenesis, lipid metabolic process, carbohydrate biosynthetic process, AMPK signaling pathway, and non-alcoholic fatty liver disease. These findings were consistent with previous research, highlighting that cold exposure significantly enhances thermogenesis and glycolipid energy metabolism in both pig breeds. Moreover, beyond regulating heat production and energy metabolism, MS pigs exhibited additional adaptations under cold stress. These adaptations involved the regulation of hemoglobin metabolic process, oxygen transport, adaptive immune response, regulation of transport, and responses to influenza A.
The WGCNA method was employed to identify 13 gene co-expression modules in bees in response to cold stress (20°C), representing different developmental stages and stress responses [48]. Among these modules, the dark orange module was found to be associated with autophagy-animal, endocytosis, and mitogen-activated protein kinase (MAPK) signaling pathways, suggesting that these pathways may play a crucial role in helping bees withstand low temperatures. Similarly, in plateau and plain yaks, four functional modules were identified in the DEGs through WGCNA analysis [49]. The turquoise and brown modules presented a high correlation with phenotypic traits such as weight, cooked meat percentage, and meat color. Furthermore, the transcriptome analysis of Chinese white wax scale insects under cold stress revealed their involvement in alcohol formation activity, lipid metabolism, membrane and structure maintenance, and oxidoreductase activity. The WGCNA results of Chinese white wax scale insects further determined that the genes within the module were associated with cytoskeletal proteins, cytoskeletal protein pathway, biosynthesis of unsaturated fatty acids, and glycerophospholipid metabolism [50]. Moreover, the tan module was found to contain hub genes, with the most prominent ones being heat shock protein 10 (hsp10), hsp60, hsp70, and hsp90. In this study, the WGCNA analysis revealed the pink module to be enriched in several key biological processes, including adaptive thermogenesis, gluconeogenesis regulation, calcium-mediated signaling, oxidative phosphorylation, carbon metabolism, non-alcoholic fatty liver disease, AMPK signaling pathway, and pyruvate metabolism. These findings connote that both LW and MS pigs respond to cold stress mainly by enhancing heat production and glycolipid energy metabolism, which aligns with previous research findings. The blue2 module exhibited enrichment in metabolic processes, apoptotic signaling pathway, oxidative phosphorylation, thermogenesis, Epstein-Barr virus infection, fatty acid oxidation, and immune pathways. These outcomes suggest that the specific response of MS pigs to cold stress is not limited to heat production and energy metabolism pathways but also involves the modulation of inflammation and immune-related pathways. This implies that MS pigs may possess enhanced cold adaptability through the regulation of their immune responses.
Phosphoproteins and SLN regulate SERCA pump activity, and their interaction with SERCA alters the kinetics of the calcium cycle, promoting thermogenesis [51,52]. UCP3, on the other hand, has been identified as a mediator of NST in pigs, as it dissipates proton gradients in mitochondria[9]. Recent murine studies have shown that cold exposure increases systemic energy expenditure and improves glucose metabolism by inducing sympathetic nervous system activity and recruiting brown adipocytes [53]. Additionally, cold exposure has demonstrated numerous effects on liver lipid metabolism and microbiome composition, contributing to the maintenance of systemic metabolic homeostasis along with BAT in mice [54]. In line with these findings, our study demonstrated that cold stress promotes inefficient calcium cycle thermogenesis, UCP3 uncoupling thermogenesis, and glycolipid energy metabolism in LW and MS pigs at the mRNA and protein levels. Moreover, MS pigs displayed a relatively higher fatty acid oxidation ability.
In this study, five candidate genes (ARRDC3, ADRB2, CEBPB, NR1D1, and SIK1) were identified in response to cold stress in LW and MS pigs. Arrestin domain-containing protein 3 (ARRDC3) is known to contribute to the regulation of the mouse β-adrenergic receptor signaling pathway, alongside influencing the energy consumption of brown and white adipose tissue [55,56]. Adrenoceptor beta 2 (ADRB2) specifically binds to endogenous ligands, including EPI and NE to regulate intracellular cAMP levels. Furthermore, ADRB2 activation catalyzes the exchange of GDP and GTP on G protein-coupled receptors and activates adenylate cyclase, converting ATP into cAMP, thereby activating the Raf/MAPK/ERK signaling pathway [57]. CCAAT enhancer binding protein beta (CEBPB) is integral in the adipogenic differentiation process and facilitates the regulation of PPAR expression [58,59]. Studies have concluded that acute cold exposure stimulates adrenaline release and leads to an increase in CEBPB mRNA levels in brown fat of rats. Treating brown adipocytes with NE also results in an elevation of CEBPB mRNA levels [60]. Nuclear receptor subfamily 1 group D member 1 (NR1D1) is a critical component of the mammalian circadian clock system. Its endogenous ligand is heme, a metabolite involved in mitochondrial respiration and cellular redox balance, making NR1D1 serve as a sensor for cellular metabolism [61]. Salt-inducible kinases (SIKs) are members of the AMPK family and play an important regulatory role in glycolipid metabolism [62,63]. Previous research has shown that manipulating SIK1 expression can impact gluconeogenesis [64] and fat metabolism in the liver of mice [65]. Silencing SIK1 leads to the up-regulation of genes involved in adipogenesis, whereas overexpression has the opposite effect, suggesting a specific regulatory mechanism for SIK1 in modulating metabolic processes related to adipocyte formation and fat metabolism. The current study observed that five candidate genes shared by LW and MS pigs were up-regulated in response to cold stress in both pig longissimus dorsi and porcine skeletal muscle satellite cells. These findings imply that both pig breeds employ a similar response mechanism to cope with cold stress, involving an increase in heat production and glycolipid energy metabolism.
ACTC1, PRSS8, and IL18 were identified as candidate genes specifically responding to cold stress in MS pigs. Actin alpha cardiac muscle 1 (ACTC1), a cardiac α-actin, serves as the main protein of cardiac myofilament, responsible for the heart’s systolic function [66]. Studies in mice have demonstrated that overexpression of ACTC1 can alleviate muscle dysplasia caused by skeletal muscle α-actin deficiency [67]. Interestingly, the amount of ACTC1 and MHCB transcripts in chickens remained unchanged under cold stimulation [68], whereas the regulation of the ACTC1 gene in mammals appears to be strongly affected by cold stress [69]. In Altay lambs, known for their cold adaptability, the expression of PVALB, TNNC1, MYL2, and ACTC1 genes related to muscle contraction was observed to be higher compared to Hu lambs under cold stress [70]. Interleukin-18 (IL-18), initially believed to stimulate Th1 cells to produce IFN-γ through CD3, especially in the presence of IL-12, has been found to have additional functions [71]. IL-18, in conjunction with IL-3, induces mast cells and basophils to produce IL-4 and IL-13, thereby stimulating both innate and acquired immunity [72]. Moreover, IL-18 induces phosphorylation of PI3K/Akt/S6 and mammalian target of rapamycin (mTOR), which in turn influences the expression of Bcl-xL and Bcl2 [73,74]. Protease serine 8 (PRSS8), also known as prostasin, is a serine protease that is a vital element of various biological processes, including protein degradation and digestion, protein processing, and tissue remodeling. Recent research indicates that PRSS8 can influence liver sensitivity to insulin through its involvement in regulating TLR4-mediated signaling pathways [75]. When PRSS8 is overexpressed in the liver, it leads to the phosphorylation of the ERK pathway, resulting in improved glucose and lipid metabolism, and can help ameliorate fatty liver disorders [76]. In this study, three candidate genes specifically responsive to cold stress in MS pigs showed increased expression in pig longissimus dorsi and skeletal muscle satellite cells under cold stress. This indicated that MS pigs, when exposed to cold stress, not only enhance heat production and glycolipid energy metabolism but also rely on enhancing their immune response to bolster their body resistance and achieve better cold adaptability. This study also found that the knockdown of PRSS8 in porcine skeletal muscle satellite cells had a significant down-regulatory effect on genes and proteins related to thermogenesis, glucose metabolism, and fatty acid oxidation. Furthermore, it inhibited ERK phosphorylation and reduced energy metabolism. These results are corroborated by the functional studies of PRSS8 in other animals, demonstrating a positive correlation between PRSS8 and thermogenesis as well as energy metabolism.

4. Materials and Methods

4.1. Animals

The current study was conducted in the Department of Animal Science, Shanxi Agricultural University. The research protocol for the present experiment was approved by the Care and Use Committee of Shanxi Agricultural University (license number: SXAU-EAW-2021MS.P.052801). The method was based on the "Guidelines for the Care and Use of Laboratory Animals" of the Ministry of Agriculture.
A total of twelve healthy male MS pigs and LW pigs, aged 90 days, were selected from the Datong pig breeding farm (Datong, China). They were divided into four groups: MS pig normal-temperature group (MS-25), MS pig low-temperature group (MS-4), LW pig normal-temperature group (LW-25), and LW pig low-temperature group (LW-4), with three pigs in each group. The pigs were reared under normal temperature (25±1°C) and low temperature (4±1°C), respectively. For the experiment, the pigs were kept in an environment control chamber, with each artificial climate chamber operating on a 12-hour day and night cycle. A preparatory feeding period of 7 days was followed by a 4-day experimental period. Throughout the experiment, the pigs had ad libitum access to feed and water.

4.2. Phenotypic Detection and Sample Collection

During the experiment, each pig underwent three daily measurements at 3 p.m., including recording the core temperature, calculating the respiratory rate based on abdominal fluctuations per minute, and determining the shivering frequency by counting shivers per minute. Blood samples were collected from the anterior vena cava of each pig at both the beginning and end of the experiment using anticoagulant-free tubes. After the experiment concluded, the animals were anesthetized using electric shock, and the carotid artery was immediately severed to ensure euthanasia. Aseptic techniques were employed to obtain the longissimus dorsi muscle from the pig. The collected muscle was promptly subdivided into cryopreservation tubes and immediately placed in liquid nitrogen for rapid freezing. The tubes were then stored at -80℃ in a refrigerator. Additionally, the 1cm3 sections of the longissimus dorsi muscle were wrapped in tin foil paper coated with optimal cutting temperature compound and immediately placed in liquid nitrogen. Furthermore, 1mm3 segments of the longissimus dorsi muscle were excised and placed in Eppendorf (EP) tubes containing 2.5% glutaraldehyde. The tubes were then stored at 4℃ in a refrigerator for further processing.

4.3. Serum Biochemistry and Hormone Analyses

Enzyme activities related to various metabolic processes were assessed using ELISA kits (Mlbio, Shanghai, China). Specific detection kits were employed to measure the levels of heat-related hormones and enzyme activities, including norepinephrine (NE), epinephrine (EPI), and Ca2+/Mg2+-ATPase. In addition, enzyme activities associated with glucose metabolism, including lactate dehydrogenase (LDH), phosphofructokinase 1 (PFK1), hexokinase 2 (HK2), and pyruvate kinase M2 isoenzyme (M2-PK), were analyzed using the respective detection kit. Furthermore, fatty acid oxidation-related enzyme activity was examined using a citrate synthase (CS) and succinate dehydrogenase (SDH) detection kit.

4.4. Electron Microscopy and Immunohistochemistry

The longissimus dorsi muscle tissues were stored in EP tubes containing 2.5% glutaraldehyde and fixed in a 1% OsO4 solution for one hour. Subsequently, the samples were processed through infiltration, embedding, slicing, and staining with 2% uranyl acetate and lead citrate to obtain ultrathin sections. The images were then observed and captured by transmission electron microscopy (HT7800, Hitachi, Tokyo, Japan). This allowed for the visualization of mitochondrial abundance and structural integrity in the longissimus dorsi muscle.
The freezing microtome was pre-cooled to -20°C beforehand, and the samples were allowed to equilibrate for 20 minutes to reach a stable temperature. Tissue sections were obtained with a fine adjustment of 7 μm to achieve the desired thickness. These sections were subsequently carefully placed onto prepared slides. Following the evaporation of tissue fog, the slides were examined under a microscope to assess their morphology. Slides displaying intact tissue structure were carefully chosen and stored at -80℃ in a refrigerator. The sections of porcine longissimus dorsi were subjected to histochemical staining using commercially available SDH (BC0955, Solarbio, Beijing, China) and LDH (BC0685, Solarbio, Beijing, China) staining kits, following the provided instructions. The histochemical staining intensity of SDH and LDH served as indicators of the tissue’s ability to regulate fatty acid oxidation and glycolysis, respectively. Image analyses were performed using the Image-Pro Plus system (Version 6.0, Media Cybernetics, Silver Spring, MD, USA).

4.5. RNA-seq

Transcriptome sequencing of the longissimus dorsi muscle from 12 pigs was conducted by Gene Denovo (Gene Denovo Co., Ltd., Guangzhou, China). The RNA samples with high purity (OD260/280≥2.0) and high integrity (RIN>8) were sequenced using Illumina Novaseq 6000, resulting in a substantial number of original reads. Fastp was employed to ensure the quality of the raw data, and high-quality clean reads were obtained by filtering out reads containing more than 10% unknown nucleotides (N) and more than 50% low-quality bases (Q value≤20). StringTie software was used to calculate fragments per kilobase of transcript per million mapped reads (FPKM) to quantify the expression abundance and variation of each transcript. DESeq2 software was utilized for analysis to identify DEGs between the two groups. The screening threshold was set as log2|FoldChange|≥1, with P<0.05. To provide functional insights, the DEGs were annotated to GO terms and KEGG pathways using annotation and visualization methods. The abundance of DEGs was normalized based on the z-score and represented by -log10P value and fold enrichment. The sequencing data was uploaded to NCBI (accession: PRJNA999915).

4.6. Weighted Gene Co-Expression Network Analysis

For the analysis of all gene expression data, the WGCNA R package (v1.70) was used to construct a weighted gene co-expression network. The dynamic segmentation method was applied to identify and classify modules according to their expression patterns. The following parameters were set: cut height=0.8, minsize=10, while default values were used for other parameters. Eigenvalue analysis was carried out to assess the correlation between modules and sample groups. Further, GO and KEGG enrichment analyses were performed based on the genes within the target modules.

4.7. Cell Culture

Porcine skeletal muscle satellite cells were retrieved from a liquid nitrogen tank and rapidly thawed in a 37°C water bath. After thawing, the cells were immediately centrifuged at 800 r/min at room temperature for five minutes using a centrifuge (5810R, Eppendorf, Hamburg, Germany). The supernatant was then discarded, and the cells were resuspended in a complete medium consisting of 10% fetal bovine serum, 1% penicillin-streptomycin, and high glucose Dulbecco's Modified Eagle Medium (DMEM) before being inoculated into 6-well plates. Using the porcine protease serine 8 (PRSS8) gene sequence as a reference, four siRNA sequences were custom-synthesized by RiboBio Co., Ltd (Guangzhou, China), and the sequence information is provided in the supplementary information (Supplementary Table S1). The transfection of the four siRNA into porcine skeletal muscle satellite cells was accomplished using LipofectamineTM 2000 (Thermo Fisher Scientific, Waltham, MA, US) when cell fusion reached approximately 70%. Following transfection, the complete medium was replaced and maintained for 6-12 hours. Upon achieving around 90% cell fusion, the differentiation medium (2% HS, 1% penicillin-streptomycin, high glucose DMEM) was introduced to induce myogenic differentiation, and the medium was refreshed every 2 days. After 6 days of differentiation, cells were harvested for the evaluation of heat production and energy metabolism marker gene expressions. For the cold exposure experiment, both untransfected and differentiated porcine skeletal muscle satellite cells were assigned to two groups: a control group maintained at 37±0.5°C and a cold exposure group at 32±0.5°C for 8 hours. It was previously determined that 32°C represents clinically mild hypothermia levels.

4.8. Quantitative Real-Time Polymerase Chain Reaction

The primers for each gene were designed based on the available mRNA sequences in GenBank, utilizing the Primer-BLAST module of the NCBI website. For normalization, 18S was used as the internal reference, and the primer information is provided in the supplementary information (Supplementary Table S2). These primers were synthesized by Bioengineering Co., Ltd (Shanghai, China). Total RNA from the collected longissimus dorsi muscle tissue and porcine skeletal muscle satellite cells was extracted using Trizol (Takara, Kyoto, Japan). Subsequently, 1-2 μg of total RNA was subjected to reverse transcription, generating complementary DNA (cDNA) with the QuantiTect Reverse Transcription Kit (Qiagen GmbH, Hilden, Germany). The cDNA was stored at -20°C in a refrigerator. For qRT-PCR, SYBR Premix Ex Taq (Takara, Kyoto, Japan) was used with a 10× dilution of the cDNA. Each sample was subjected to three technical replicates for accurate analysis. The relative mRNA expression levels were determined using the 2-ΔΔCT calculation method.

4.9. Western Blotting

All skeletal muscle tissues were first ground into a fine powder in liquid nitrogen. The tissue powder and satellite cells were then lysed in total protein lysis buffer, containing radioimmunoprecipitation assay and phenylmethylsulphonyl fluoride, at a ratio of 100:1. The lysates were incubated on ice for 30 minutes. After that, the samples were centrifuged at 12,000 r/min at 4°C for 15 minutes, and the resulting supernatants were collected. They were then denatured at 100°C for 10 minutes with 5× sodium dodecyl sulfate (SDS) loading buffer at a ratio of 4:1. Equal amounts of protein were separated using a 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) gel, running at 80 V for 30 minutes and 120 V for 90 minutes. The proteins were transferred onto 0.45 μm polyvinylidene difluoride membranes (Millipore, Billerica, MA, USA) at 120 V for 90 minutes. Subsequently, the membranes were incubated in 5% skim milk powder solution at room temperature and gently shaken for 1 hour to block non-specific binding. Primary antibodies, including SERCA2 (1:1000, Bioss, Beijing, China), CASQ2 (1:1000, Proteintech, Wuhan, China), UCP3 (1:5000, Abcam, Cambridge, UK), PGC-1α (1:1000, Proteintech, Wuhan, China), CPT1B (1:1000, Proteintech, Wuhan, China), PRSS8 (1:1000, Affinity Biosciences, Changzhou, China), ERK1/2 (1:1000, Proteintech, Wuhan, China), Phospho-ERK1/2 (1:1000, Proteintech, Wuhan, China), PFK1 (1:1000, Bioss, Beijing, China), HADHB (1:1000, Bioss, Beijing, China), β-actin (1:5000, Abcam, Cambridge, UK), and β-Tubulin (1:1000, Affinity Biosciences, Changzhou, China), were used to probe the membranes either for 2 hours at room temperature or overnight at 4℃. After washing the membranes with 0.05% Tris-buffered saline containing Tween-20 buffer (TBST), appropriate infrared (IR)-linked secondary antibodies (LI-COR, Lincoln, NE, USA) were used to probe the membranes for 2 hours at room temperature. Following another wash with 0.05% TBST, the membranes were scanned using the LI-COR Odyssey scanner, and the gray value of the protein bands was calculated and analyzed using Image J software.

4.10. Statistical Analysis

Statistical analysis was performed using IBM SPSS Statistics Version 26.0 (IBM Corp, Armonk, NY, USA) and GraphPad Prism Version 8.0 (GraphPad Software Inc., San Diego, CA, USA). To determine the significance of differences between two groups, the two-tailed unpaired Student’s t-test was employed. For analyzing the significance of differences among multiple groups, one-way analysis of variance (ANOVA) with Duncan's multiple range test was employed. Statistically significant levels were set as **P<0.01 and *P<0.05. All data were presented as mean±standard error of mean (SEM).

5. Conclusions

LW and MS pigs exhibited a decrease in body temperature and body weight under cold stress, along with increased shivering and respiratory frequencies. However, despite these changes, MS pigs maintained stable physiological homeostasis, demonstrating their cold adaptability. LW pigs responded to cold stress chiefly by regulating heat production, glucose metabolism, and fatty acid oxidation. Meanwhile, MS pigs exhibited a specific response to cold stress, showcasing their strong fatty acid oxidation ability and immune response, which contributed to their enhanced adaptability. This study identified ARRDC3, ADRB2, CEBPB, NR1D1, and SIK1 as common candidate genes responsive to cold stress in LW and MS pigs. In contrast, ACTC1, PRSS8, and IL18 were identified as specific candidate genes that uniquely responded to cold stress in MS pigs. Furthermore, it was deduced that PRSS8 plays a considerable role in regulating heat production, glucose metabolism, and fatty acid oxidation marker genes and proteins in porcine skeletal muscle satellite cells. Moreover, PRSS8 was found to influence energy metabolism through its impact on ERK phosphorylation. By elucidating the specific genetic and metabolic pathways involved in the response to cold stress, this study contributes valuable insights that can help in developing strategies for enhancing the cold resistance and adaptability of pigs, ultimately leading to improved livestock management and increased economic gains in the pig industry.

Supplementary Materials

The following supporting information can be downloaded at: Preprints.org.

Author Contributions

Conceptualization, G.C. and Y.Y.; Data curation, Y.C., Y.Z. and W.Z.; Formal analysis, W.L., N.Z., M.S. and Y.Z.; Funding acquisition, G.C.; Investigation, W.L., C.C., C.L. and P.G.; Methodology, W.L., X.G.; B.L. and S.K.; Project administration, G.C. and Y.Y.; Resources, X.G., B.L. G.C. and Y.Y.; Software, Y.C., N.Z. and C.L.; Supervision, X.G., B.L., G.C. and Y.Y.; Validation, G.C. and Y.Y.; Visualization, W.L. and Y.C.; Writing—original draft, W.L.; Writing—review and editing, G.C. and Y.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the following funding sources: the National Natural Science Foundation of China (31872336), and the Postgraduate Innovation Project of Shanxi Province (2021Y320), the Special Funds for Scholars Support Program of Shanxi Province (2016, 2017), the Cultivation Project for Outstanding Scientific Research Achievements in Universities of Shanxi Province, the Key R&D Program of Shanxi Province (202102140601005), the Biological Breeding Engineering program of Shanxi Agricultural University (YZGC128).

Institutional Review Board Statement

The animal study protocol was approved by the Care and Use Committee of Shanxi Agricultural University (license number: SXAU-EAW-2021MS.P.052801).

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are available from the corresponding authors upon reasonable request. The data are not publicly available due to the fact that we are conducting further experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Toghiani, S.; Hay, E.; Fragomeni, B.; Rekaya, R.; Roberts, AJ. Genotype by environment interaction in response to cold stress in a composite beef cattle breed. Animal. 2020, 31, 1–12. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, D.; Cheng, X.; Fang, H.; Ren, Y.; Li, X.; Ren, W.; Xue, B.; Yang, C. Effect of cold stress on ovarian & uterine microcirculation in rats and the role of endothelin system. Reprod. Biol. Endocrin. 2020, 18, 29. [Google Scholar] [CrossRef]
  3. Le Dividich, J.; Noblet, J. Colostrum intake and thermoregulation in the neonatal pig in relation to environmental temperature. Biol. Neonate. 1981, 40, 167–174. [Google Scholar] [CrossRef] [PubMed]
  4. Iida, R.; Koketsu, Y. Climatic factors associated with peripartum pig deaths during hot and humid or cold seasons. Prev. Vet. Med. 2014, 115, 166–172. [Google Scholar] [CrossRef] [PubMed]
  5. Čobanović, N.; Stajković, S.; Blagojević, B.; Betić, N.; Dimitrijević, M.; Vasilev, D.; Karabasil, N. The effects of season on health, welfare, and carcass and meat quality of slaughter pigs. Int. J. Biometeorol. 2020, 64, 1899–1909. [Google Scholar] [CrossRef] [PubMed]
  6. Li, M.; Tian, S.; Jin, L.; Zhou, G.; Li, Y.; Zhang, Y.; Wang, T.; Yeung, C.; Chen, L.; Ma, J.; et al. Genomic analyses identify distinct patterns of selection in domesticated pigs and Tibetan wild boars. Nat Genet. 2013, 45, 1431–1438. [Google Scholar] [CrossRef] [PubMed]
  7. Liu, Y.; Yang, X.; Jing, X.; He, X.; Wang, L.; Liu, Y.; Liu, D. Transcriptomics Analysis on Excellent Meat Quality Traits of Skeletal Muscles of the Chinese Indigenous Min Pig Compared with the Large White Breed. Int. J. Mol. Sci. 2017, 19, 21. [Google Scholar] [CrossRef] [PubMed]
  8. Zheng, Q.; Lin, J.; Huang, J.; Zhang, H.; Zhang, R.; Zhang, X.; Cao, C.; Hambly, C.; Qin, G.; Yao, J.; et al. UCP1Reconstitution of using CRISPR/Cas9 in the white adipose tissue of pigs decreases fat deposition and improves thermogenic capacity. Proc. Natl. Acad. Sci. U. S. A. 2017, 114, E9474–E9482. [Google Scholar] [CrossRef] [PubMed]
  9. Lin, J.; Cao, C.; Tao, C.; Ye, R.; Dong, M.; Zheng, Q.; Wang, C.; Jiang, X.; Qin, G.; Yan, C.; et al. Cold adaptation in pigs depends on UCP3 in beige adipocytes. J. Mol. Cell. Biol. 2017, 9, 364–375. [Google Scholar] [CrossRef]
  10. Sun, G.; Song, X.; Zou, Y.; Teng, T.; Jiang, L.; Shi, B. Dietary Glucose Ameliorates Impaired Intestinal Development and Immune Homeostasis Disorders Induced by Chronic Cold Stress in Pig Model. Int. J. Mol. Sci. 2022, 23, 7730. [Google Scholar] [CrossRef]
  11. Lu, C.; Liu, Y.; Ma, Y.; Wang, S.; Cai, C.; Yang, Y.; Zhao, Y.; Liang, G.; Cao, G.; Li, B.; et al. Comparative Evaluation of the Ileum Microbiota Composition in Piglets at Different Growth Stages. Front. Microbiol. 2021, 12, 765691. [Google Scholar] [CrossRef] [PubMed]
  12. Cai, C.; Li, M.; Zhang, Y.; Meng, S.; Yang, Y.; Gao, P.; Guo, X.; Cao, G.; Li, B. Comparative Transcriptome Analyses of Between Pig Breeds Differing in Muscle Characteristics. Front. Genet. 2020, 11, 526309. [Google Scholar] [CrossRef] [PubMed]
  13. Liu, T.; Guo, Y.; Lu, C.; Cai, C.; Gao, P.; Cao, G.; Li, B.; Guo, X.; Yang, Y. Effect of Different Pig Fecal Microbiota Transplantation on Mice Intestinal Function and Microbiota Changes During Cold Exposure. Front. Vet. Sci. 2022, 9, 805815. [Google Scholar] [CrossRef]
  14. Guo, Y.; Liu, T.; Li, W.; Zhang, W.; Cai, C.; Lu, C.; Gao, P.; Cao, G.; Li, B.; Guo, X.; et al. Effects of Low-Ambient-Temperature Stimulation on Modifying the Intestinal Structure and Function of Different Pig Breeds. Animals (Basel). 2022, 12, 2740. [Google Scholar] [CrossRef]
  15. Fuller-Jackson, J.; Henry, B. Adipose and skeletal muscle thermogenesis: studies from large animals. J. Endocrinol. 2018, 237, R99–R115. [Google Scholar] [CrossRef] [PubMed]
  16. Cheshire, W.P. Thermoregulatory disorders and illness related to heat and cold stress. Auton. Neurosci. 2016, 196, 91–104. [Google Scholar] [CrossRef] [PubMed]
  17. Halvorson, I.; Thornhill, J. Posterior hypothalamic stimulation of anesthetized normothermic and hypothermic rats evokes shivering thermogenesis. Brain. Res. 1993, 610, 208–215. [Google Scholar] [CrossRef]
  18. Nabi, G.; Xing, D.; Sun, Y.; Zhang, Q.; Li, M.; Jiang, C.; Ahmad, I.M.; Wingfield, J.C.; Wu, Y.; Li, D. Coping with extremes: High-altitude sparrows enhance metabolic and thermogenic capacities in the pectoralis muscle and suppress in the liver relative to their lowland counterparts. Gen. Comp. Endocrinol. 2021, 313, 113890. [Google Scholar] [CrossRef]
  19. Wang, S.; Gopinath, T.; Larsen, E.K.; Weber, D.K.; Walker, C.; Uddigiri, V.R.; Mote, K.R.; Sahoo, S.K.; Periasamy, M.; Veglia, G. Structural basis for sarcolipin's regulation of muscle thermogenesis by the sarcoplasmic reticulum Ca-ATPase. Sci. Adv. 2021, 7, eabi7154. [Google Scholar] [CrossRef]
  20. Gorski, P.; Ceholski, D.; Young, H. Structure-Function Relationship of the SERCA Pump and Its Regulation by Phospholamban and Sarcolipin. Adv. Exp. Med. Biol. 2017, 981, 77–119. [Google Scholar] [CrossRef]
  21. Bal, N.C.; Maurya, S.K.; Singh, S.; Wehrens, X.H.; Periasamy, M. Increased Reliance on Muscle-based Thermogenesis upon Acute Minimization of Brown Adipose Tissue Function. J. Biol. Chem. 2016, 291, 17247–17257. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, M.; Zhang, X.; Wang, C.; Wang, D. Lasiopodomys brandtiiRecruitment of Muscle Genes as an Effect of Brown Adipose Tissue Ablation in Cold-Acclimated Brandt's Voles. Int. J. Mol. Sci. 2022, 24, 342. [Google Scholar] [CrossRef] [PubMed]
  23. Nowack, J.; Vetter, S.; Stalder, G.; Painer, J.; Kral, M.; Smith, S.; Le, M.; Jurcevic, P.; Bieber, C.; Arnold, W.; et al. Muscle nonshivering thermogenesis in a feral mammal. Sci. Rep. 2019, 9, 6378. [Google Scholar] [CrossRef] [PubMed]
  24. Fan, G.; Li, Y.; Ma, F.; Zhao, R.; Yang, X. Zinc-α2-glycoprotein promotes skeletal muscle lipid metabolism in cold-stressed mice. Endocr. J. 2021, 68, 53–62. [Google Scholar] [CrossRef]
  25. Codella, R.; Alves, T.; Befroy, D.; Choi, C.; Luzi, L.; Rothman, D.; Kibbey, R.; Shulman, G. Overexpression of UCP3 decreases mitochondrial efficiency in mouse skeletal muscle in vivo. FEBS Lett. 2023, 597, 309–319. [Google Scholar] [CrossRef]
  26. Hilse, KE.; Kalinovich, AV.; Rupprecht, A.; Smorodchenko, A.; Zeitz, U.; Staniek, K.; Erben, RG.; Pohl, EE. The expression of UCP3 directly correlates to UCP1 abundance in brown adipose tissue. Biochim. Biophys. Acta. 2016, 1857, 72–78. [Google Scholar] [CrossRef]
  27. Maurya, S.K.; Herrera, J.L.; Sahoo, S.K.; Reis, F.C.G.; Vega, R.B.; Kelly, D.P.; Periasamy, M. Sarcolipin Signaling Promotes Mitochondrial Biogenesis and Oxidative Metabolism in Skeletal Muscle. Cell. Rep. 2018, 24, 2919–2931. [Google Scholar] [CrossRef]
  28. Hangalapura, B.N.; Kaiser, M.G.; Poel, J.J.; Parmentier, H.K.; Lamont, S.J. Cold stress equally enhances in vivo pro-inflammatory cytokine gene expression in chicken lines divergently selected for antibody responses. Dev. Comp. Immunol. 2006, 30, 503–511. [Google Scholar] [CrossRef]
  29. Górski, J.; Kuryliszyn, A.; Wereszczyńska, U. Effect of acute cold exposure on the mobilization of intramuscular glycogen and triglycerides in the rat. Acta. Physiol. Pol. 1981, 32, 755–759. [Google Scholar]
  30. Umehara, T.; Murase, T.; Abe, Y.; Yamashita, H.; Shibaike, Y.; Kagawa, S.; Yamamoto, T.; Ikematsu, K. Identification of potential markers of fatal hypothermia by a body temperature-dependent gene expression assay. Int. J. Legal. Med. 2019, 133, 335–345. [Google Scholar] [CrossRef]
  31. Hou, L.; Shi, J.; Cao, L.; Xu, G.; Hu, C.; Wang, C. Pig has no uncoupling protein 1. Biochem. Biophys. Res. Commun. 2017, 487, 795–800. [Google Scholar] [CrossRef]
  32. Brosh, A. Heart rate measurements as an index of energy expenditure and energy balance in ruminants: a review. J. Anim. Sci. 2007, 85, 1213–1227. [Google Scholar] [CrossRef]
  33. Le Dividich, J.; Mormède, P.; Catheline, M.; Caritez, J.C. Body composition and cold resistance of the neonatal pig from European (Large White) and Chinese (Meishan) breeds. Biol. Neonate. 1991, 59, 268–277. [Google Scholar] [CrossRef] [PubMed]
  34. Heldmaier, G. Cold adaptation by short daily cold exposures in the young pig. J. Appl. Physiol. 1974, 36, 163–168. [Google Scholar] [CrossRef]
  35. Herpin, P.; Lossec, G.; Schmidt, I.; Cohen-Adad, F.; Duchamp, C.; Lefaucheur, L.; Goglia, F.; Lanni, A. Effect of age and cold exposure on morphofunctional characteristics of skeletal muscle in neonatal pigs. Pflugers. Arch. 2002, 444, 610–618. [Google Scholar] [CrossRef] [PubMed]
  36. Le Dividich, J.; Noblet, J. Thermoregulation and energy metabolism in the neonatal pig. Ann. Rech. Vet. 1983, 14, 375–381. [Google Scholar] [PubMed]
  37. Berthon, D.; Herpin, P.; Bertin, R.; De Marco, F.; le Dividich, J. Metabolic changes associated with sustained 48-hr shivering thermogenesis in the newborn pig. Comp. Biochem. Physiol. B. Biochem. Mol. Biol. 1996, 114, 327–335. [Google Scholar] [CrossRef]
  38. Ren, J.; Long, Y.; Liu, R.; Song, G.; Li, Q.; Cui, Z. Characterization of Biological Pathways Regulating Acute Cold Resistance of Zebrafish. Int. J. Mol. Sci. 2021, 22, 3028. [Google Scholar] [CrossRef]
  39. Hangalapura, B.; Nieuwland, M.; Buyse, J.; Kemp, B.; Parmentier, H. Effect of duration of cold stress on plasma adrenal and thyroid hormone levels and immune responses in chicken lines divergently selected for antibody responses. Poult. Sci. 2004, 83, 1644–1649. [Google Scholar] [CrossRef]
  40. Singh, A.K.; Aryal, B.; Chaube, B.; Rotllan, N.; Varela, L.; Horvath, T.L.; Suárez, Y.; Fernández-Hernando, C. Brown adipose tissue derived ANGPTL4 controls glucose and lipid metabolism and regulates thermogenesis. Mol. Metab. 2018, 11, 59–69. [Google Scholar] [CrossRef]
  41. Faure, J.; Lebret, B.; Bonhomme, N.; Ecolan, P.; Kouba, M.; Lefaucheur, L. Metabolic adaptation of two pig muscles to cold rearing conditions. J. Anim. Sci. 2013, 91, 1893–1906. [Google Scholar] [CrossRef] [PubMed]
  42. Cheviron, Z.; Bachman, G.; Connaty, A.; McClelland, G.; Storz, J. Regulatory changes contribute to the adaptive enhancement of thermogenic capacity in high-altitude deer mice. Proc. Natl. Acad. Sci. U. S. A. 2012, 109, 8635–8640. [Google Scholar] [CrossRef]
  43. Xie, S.; Yang, X.; Wang, D.; Zhu, F.; Yang, N.; Hou, Z.; Ning, Z. Thyroid transcriptome analysis reveals different adaptive responses to cold environmental conditions between two chicken breeds. PLoS One. 2018, 13, e0191096. [Google Scholar] [CrossRef]
  44. Ge, G.; Long, Y.; Shi, L.; Ren, J.; Yan, J.; Li, C.; Li, Q.; Cui, Z. Transcriptomic profiling revealed key signaling pathways for cold tolerance and acclimation of two carp species. BMC Genomics. 2020, 21, 539. [Google Scholar] [CrossRef] [PubMed]
  45. Xu, Z.; You, W.; Zhou, Y.; Chen, W.; Wang, Y.; Shan, T. Cold-induced lipid dynamics and transcriptional programs in white adipose tissue. BMC Biol. 2019, 17, 74. [Google Scholar] [CrossRef]
  46. Smith, S.D.; Kawash, J.K.; Karaiskos, S.; Biluck, I.; Grigoriev, A. Evolutionary adaptation revealed by comparative genome analysis of woolly mammoths and elephants. DNA Res. 2017, 24, 359–369. [Google Scholar] [CrossRef] [PubMed]
  47. Jansen, H.T.; Trojahn, S.; Saxton, M.W.; Quackenbush, C.R.; Evans Hutzenbiler, B.D.; Nelson, O.L.; Cornejo, O.E.; Robbins, C.T.; Kelley, J.L. Hibernation induces widespread transcriptional remodeling in metabolic tissues of the grizzly bear. Commun. Biol. 2019, 2, 336. [Google Scholar] [CrossRef]
  48. Zhu, C.; Xu, X.; Zhou, S.; Zhou, B.; Liu, Y.; Xu, H.; Tian, Y.; Zhu, X. Apis melliferaWGCNA based identification of hub genes associated with cold response and development in metamorphic pupae. Front. Physiol. 2023, 14, 1169301. [Google Scholar] [CrossRef]
  49. Wang, J.; Guan, J.; Yixi, K.; Shu, T.; Chai, Z.; Wang, J.; Wang, H.; Wu, Z.; Cai, X.; Zhong, J.; et al. Comparative transcriptome analysis of winter yaks in plateau and plain. Reprod. Domest. Anim. 2022, 57, 64–71. [Google Scholar] [CrossRef] [PubMed]
  50. Zhang, H.; Liu, W.; An, J.; Yang, P.; Guo, L.; Li, Y.; Lv, J.; Yu, S. Transcriptome analyses and weighted gene coexpression network analysis reveal key pathways and genes involved in the rapid cold resistance of the Chinese white wax scale insect. Arch. Insect. Biochem. Physiol. 2021, 107, e21781. [Google Scholar] [CrossRef]
  51. Periasamy, M.; Maurya, S.K.; Sahoo, S.K.; Singh, S.; Sahoo, S.K.; Reis, F.C.G.; Bal, N.C. Role of SERCA Pump in Muscle Thermogenesis and Metabolism. Comp. Physiol. 2017, 7, 879–890. [Google Scholar] [CrossRef]
  52. Rowland, L.A.; Bal, N.C.; Periasamy, M. The role of skeletal-muscle-based thermogenic mechanisms in vertebrate endothermy. Biol. Rev. Camb. Philos. Soc. 2015, 90, 1279–1297. [Google Scholar] [CrossRef] [PubMed]
  53. Iwen, K.A.; Backhaus, J.; Cassens, M.; Waltl, M.; Hedesan, O.C.; Merkel, M.; Heeren, J.; Sina, C.; Rademacher, L.; Windjäger, A.; et al. Cold-Induced Brown Adipose Tissue Activity Alters Plasma Fatty Acids and Improves Glucose Metabolism in Men. J. Clin. Endocr. Metab. 2017, 102, 4226–4234. [Google Scholar] [CrossRef] [PubMed]
  54. Grefhorst, A.; van den Beukel, J.C.; Dijk, W.; Steenbergen, J.; Voortman, G.J.; Leeuwenburgh, S.; Visser, T.J.; Kersten, S.; Friesema, E.C.H.; Themmen, A.P.N.; et al. Multiple effects of cold exposure on livers of male mice. J. Endocrinol. 2018, 238, 91–106. [Google Scholar] [CrossRef] [PubMed]
  55. Xiao, J.; Shi, Q.; Li, W.; Mu, X.; Peng, J.; Li, M.; Chen, M.; Huang, H.; Wang, C.; Gao, K.; et al. ARRDC1 and ARRDC3 act as tumor suppressors in renal cell carcinoma by facilitating YAP1 degradation. Am. J. Cancer. Res. 2018, 8, 132–143. [Google Scholar]
  56. Lin, C.; MacGurn, J.; Chu, T.; Stefan, C.; Emr, S. Arrestin-related ubiquitin-ligase adaptors regulate endocytosis and protein turnover at the cell surface. Cell. 2008, 135, 714–725. [Google Scholar] [CrossRef]
  57. He, R.; He, Y.; Tang, Y.; Zhou, H.; McLeod, H.; Liu, J. The potential anticancer effect of beta-blockers and the genetic variations involved in the interindividual difference. Pharmacogenomics. 2016, 17, 74–79. [Google Scholar] [CrossRef]
  58. Smas, C.; Sul, H. Control of adipocyte differentiation. Biochem. J. 1995, 697–710. [Google Scholar] [CrossRef]
  59. Tang, W.; Weng, L.; Wang, X.; Liu, C.; Hu, G.; Yin, S.; Tao, Y.; Hong, N.; Guo, H.; Liu, W.; et al. The Mediator subunit MED20 organizes the early adipogenic complex to promote development of adipose tissues and diet-induced obesity. Cell Rep. 2021, 36, 109314. [Google Scholar] [CrossRef]
  60. Lee, J.; Chae, S.; Nan, J.; Koo, Y.; Lee, S.; Park, Y.; Hwang, D.; Han, W.; Lee, D.; Kim, Y.; et al. SENP2 suppresses browning of white adipose tissues by de-conjugating SUMO from C/EBPβ. Cell Rep. 2022, 38, 110408. [Google Scholar] [CrossRef]
  61. Everett, L.; Lazar, M. Nuclear receptor Rev-erbα: up, down, and all around. Trends. Endocrinol. Metab. 2014, 25, 586–592. [Google Scholar] [CrossRef] [PubMed]
  62. Altarejos, J.; Montminy, M. CREB and the CRTC co-activators: sensors for hormonal and metabolic signals. Nat. Rev. Mol. Cell Biol. 2011, 12, 141–151. [Google Scholar] [CrossRef]
  63. Sakamaki, J.; Fu, A.; Reeks, C.; Baird, S.; Depatie, C.; Al Azzabi, M.; Bardeesy, N.; Gingras, A.; Yee, S.; Screaton, R. Role of the SIK2-p35-PJA2 complex in pancreatic β-cell functional compensation. Nat. Cell Biol. 2014, 16, 234–244. [Google Scholar] [CrossRef]
  64. Koo, S.; Flechner, L.; Qi, L.; Zhang, X.; Screaton, R.; Jeffries, S.; Hedrick, S.; Xu, W.; Boussouar, F.; Brindle, P.; et al. The CREB coactivator TORC2 is a key regulator of fasting glucose metabolism. Nature. 2005, 437, 1109–1111. [Google Scholar] [CrossRef] [PubMed]
  65. Yoon, Y.; Seo, W.; Lee, M.; Kim, S.; Koo, S. Salt-inducible kinase regulates hepatic lipogenesis by controlling SREBP-1c phosphorylation. J. Biol. Chem. 2009, 284, 10446–10452. [Google Scholar] [CrossRef]
  66. Pourcel, L.; Buron, F.; Arib, G.; Le Fourn, V.; Regamey, A.; Bodenmann, I.; Girod, P.; Mermod, N. Influence of cytoskeleton organization on recombinant protein expression by CHO cells. Biotechnol. Bioeng. 2020, 117, 1117–1126. [Google Scholar] [CrossRef] [PubMed]
  67. Boutilier, J.; Taylor, R.; Ram, R.; McNamara, E.; Nguyen, Q.; Goullée, H.; Chandler, D.; Mehta, M.; Balmer, L.; Laing, N.; et al. Variable cardiac α-actin (Actc1) expression in early adult skeletal muscle correlates with promoter methylation. Biochim. Biophys. Acta. Gene Regul. Mech. 2017, 1860, 1025–1036. [Google Scholar] [CrossRef]
  68. Bahadoran, S.; Hassanpour, H.; Arab, S.; Abbasnia, S.; Kiani, A. Changes in the expression of cardiac genes responsive to thyroid hormones in the chickens with cold-induced pulmonary hypertension. Poult. Sci. 2021, 100, 101263. [Google Scholar] [CrossRef]
  69. Kahaly, GJ.; Dillmann, WH. Thyroid hormone action in the heart. Endocr. Rev. 2005, 26, 704–728. [Google Scholar] [CrossRef]
  70. Ji, K.; Jiao, D.; Yang, G.; Degen, A.; Zhou, J.; Liu, H.; Wang, W.; Cong, H. Transcriptome analysis revealed potential genes involved in thermogenesis in muscle tissue in cold-exposed lambs. Front. Genet. 2022, 13, 1017458. [Google Scholar] [CrossRef]
  71. Nakanishi, K.; Yoshimoto, T.; Tsutsui, H.; Okamura, H. Interleukin-18 regulates both Th1 and Th2 responses. Annu. Rev. Immunol. 2001, 19, 423–474. [Google Scholar] [CrossRef] [PubMed]
  72. Nakanishi, K. Unique Action of Interleukin-18 on T Cells and Other Immune Cells. Front. Immunol. 2018, 9, 763. [Google Scholar] [CrossRef]
  73. El-Darawish, Y.; Li, W.; Yamanishi, K.; Pencheva, M.; Oka, N.; Yamanishi, H.; Matsuyama, T.; Tanaka, Y.; Minato, N.; Okamura, H. Frontline Science: IL-18 primes murine NK cells for proliferation by promoting protein synthesis, survival, and autophagy. J. Leukocyte. Biol. 2018, 104, 253–264. [Google Scholar] [CrossRef]
  74. Deason, K.; Troutman, T.; Jain, A.; Challa, D.; Mandraju, R.; Brewer, T.; Ward, E.; Pasare, C. BCAP links IL-1R to the PI3K-mTOR pathway and regulates pathogenic Th17 cell differentiation. J. Exp. Med. 2018, 215, 2413–2428. [Google Scholar] [CrossRef]
  75. Uchimura, K.; Hayata, M.; Mizumoto, T.; Miyasato, Y.; Kakizoe, Y.; Morinaga, J.; Onoue, T.; Yamazoe, R.; Ueda, M.; Adachi, M.; et al. The serine protease prostasin regulates hepatic insulin sensitivity by modulating TLR4 signalling. Nat. Commun. 2014, 5, 3428. [Google Scholar] [CrossRef] [PubMed]
  76. Sekine, T.; Takizawa, S.; Uchimura, K.; Miyazaki, A.; Tsuchiya, K. Liver-Specific Overexpression of Prostasin Attenuates High-Fat Diet-Induced Metabolic Dysregulation in Mice. Int. J. Mol. Sci. 2021, 22, 8314. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Physiological response patterns of (A) body temperature, (B) body temperature deviation, (C) respiratory rate and (D) shivering frequency in LW and MS pigs exposed to temperatures of 25℃ and 4℃ for 4 days. The results are expressed as mean ± SEM (n=3). LW pigs represent Large White pigs, and MS pigs represent Mashen pigs. Statistical significance was determined using the two-tailed unpaired Student’s t-test, with **P<0.01 and *P<0.05 denoting significance. A, B Different capital letters represent significant differences of LW pigs at 4℃ (P<0.05). a, b Different lowercase letters represent significant differences of MS pigs at 4℃ (P<0.05).
Figure 1. Physiological response patterns of (A) body temperature, (B) body temperature deviation, (C) respiratory rate and (D) shivering frequency in LW and MS pigs exposed to temperatures of 25℃ and 4℃ for 4 days. The results are expressed as mean ± SEM (n=3). LW pigs represent Large White pigs, and MS pigs represent Mashen pigs. Statistical significance was determined using the two-tailed unpaired Student’s t-test, with **P<0.01 and *P<0.05 denoting significance. A, B Different capital letters represent significant differences of LW pigs at 4℃ (P<0.05). a, b Different lowercase letters represent significant differences of MS pigs at 4℃ (P<0.05).
Preprints 84700 g001
Figure 2. Serum biochemistry and hormone analyses of epinephrine (A) (EPI), (B) norepinephrine (NE), (C) Ca2+/Mg2+ATPase, (D) lactate dehydrogenase (LDH), (E) phosphofructokinase 1 (PFK-1), (F) hexokinase 2 (HK2), (G) pyruvate kinase M2 isoenzyme (M2-PK), (H) succinate dehydrogenase (SDH), and (I) citrate synthase (CS) in LW and MS pigs at 25℃ and 4℃. The results are presented as mean±SEM (n=3). LW and MS refer to Large White pigs and Mashen pigs, respectively. P values were determined using the two-tailed unpaired Student’s t-test. **P<0.01; *P<0.05.
Figure 2. Serum biochemistry and hormone analyses of epinephrine (A) (EPI), (B) norepinephrine (NE), (C) Ca2+/Mg2+ATPase, (D) lactate dehydrogenase (LDH), (E) phosphofructokinase 1 (PFK-1), (F) hexokinase 2 (HK2), (G) pyruvate kinase M2 isoenzyme (M2-PK), (H) succinate dehydrogenase (SDH), and (I) citrate synthase (CS) in LW and MS pigs at 25℃ and 4℃. The results are presented as mean±SEM (n=3). LW and MS refer to Large White pigs and Mashen pigs, respectively. P values were determined using the two-tailed unpaired Student’s t-test. **P<0.01; *P<0.05.
Preprints 84700 g002
Figure 3. Effects of cold stress on the ultrastructure and enzyme activity of the longissimus dorsi muscle in LW and MS pigs at 25℃ and 4℃ (A) Transmission electron micrographs of longissimus dorsi from LW and MS pigs at 25℃ and 4℃, whereby the top panels and bottom panels show low (scale bar 5 μm) and high magnification (scale bar 1 μm), respectively. Yellow triangles and ‘M’ represent mitochondria. Succinate dehydrogenase (SDH) activity staining (B) and lactate dehydrogenase (LDH) activity staining (C) of longissimus dorsi muscle tissue in LW and MS pigs at 25℃ and 4℃ (scale bar 200 μm). The results are expressed as mean ± SEM (n=3). LW and MS represent Large White pigs and Mashen pigs, respectively. P values were determined using the two-tailed unpaired Student’s t-test. **P<0.01; *P<0.05.
Figure 3. Effects of cold stress on the ultrastructure and enzyme activity of the longissimus dorsi muscle in LW and MS pigs at 25℃ and 4℃ (A) Transmission electron micrographs of longissimus dorsi from LW and MS pigs at 25℃ and 4℃, whereby the top panels and bottom panels show low (scale bar 5 μm) and high magnification (scale bar 1 μm), respectively. Yellow triangles and ‘M’ represent mitochondria. Succinate dehydrogenase (SDH) activity staining (B) and lactate dehydrogenase (LDH) activity staining (C) of longissimus dorsi muscle tissue in LW and MS pigs at 25℃ and 4℃ (scale bar 200 μm). The results are expressed as mean ± SEM (n=3). LW and MS represent Large White pigs and Mashen pigs, respectively. P values were determined using the two-tailed unpaired Student’s t-test. **P<0.01; *P<0.05.
Preprints 84700 g003
Figure 4. Transcriptomes of longissimus dorsi muscle in LW and MS pig at 25℃ and 4℃. (A) Violin plot presenting the distribution of transcriptome sequencing samples from 12 pigs in LW and MS groups at 25℃ and 4℃. The white dot on the violin plot represents the median, and the quartile range is shown by the black bar above and below the center of the violin. (B) Cluster map of differentially expressed genes (DEGs) between MS-25 and MS-4 groups at 25℃ and 4℃. (C) Histogram depicting the number of identifiable DEGs between the four groups. (D) Volcanic map highlighting the DEGs between MS-25 and MS-4 groups at 25℃ and 4℃. (E) Venn diagram illustrating the DEGs shared between the LW pig group and the MS pig group. (F) Heatmap visualizing the DEGs shared in LW and MS pigs at 25℃ and 4℃. (G) Heatmap displaying the specific DEGs in MS pigs at 25℃ and 4℃. (H) Gene Ontology (GO) enrichment analysis of DEGs shared in LW and MS pigs at 25℃ and 4℃. (I) GO enrichment analysis of specific DEGs in MS pigs at 25℃ and 4℃. (J) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway classification of DEGs shared in LW and MS pigs at 25℃ and 4℃. (K) KEGG pathway classification of specific DEGs in MS pigs at 25℃ and 4℃. LW and MS represent Large White pigs and Mashen pigs, respectively; 25 and 4 indicate 25℃ and 4℃, respectively. The heat map illustrates the relative expression pattern of DEGs between groups, with each column representing a sample and each row representing the expression level of a single mRNA in different samples. The color range of the heat map corresponds from blue (low expression) to red (high expression).
Figure 4. Transcriptomes of longissimus dorsi muscle in LW and MS pig at 25℃ and 4℃. (A) Violin plot presenting the distribution of transcriptome sequencing samples from 12 pigs in LW and MS groups at 25℃ and 4℃. The white dot on the violin plot represents the median, and the quartile range is shown by the black bar above and below the center of the violin. (B) Cluster map of differentially expressed genes (DEGs) between MS-25 and MS-4 groups at 25℃ and 4℃. (C) Histogram depicting the number of identifiable DEGs between the four groups. (D) Volcanic map highlighting the DEGs between MS-25 and MS-4 groups at 25℃ and 4℃. (E) Venn diagram illustrating the DEGs shared between the LW pig group and the MS pig group. (F) Heatmap visualizing the DEGs shared in LW and MS pigs at 25℃ and 4℃. (G) Heatmap displaying the specific DEGs in MS pigs at 25℃ and 4℃. (H) Gene Ontology (GO) enrichment analysis of DEGs shared in LW and MS pigs at 25℃ and 4℃. (I) GO enrichment analysis of specific DEGs in MS pigs at 25℃ and 4℃. (J) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway classification of DEGs shared in LW and MS pigs at 25℃ and 4℃. (K) KEGG pathway classification of specific DEGs in MS pigs at 25℃ and 4℃. LW and MS represent Large White pigs and Mashen pigs, respectively; 25 and 4 indicate 25℃ and 4℃, respectively. The heat map illustrates the relative expression pattern of DEGs between groups, with each column representing a sample and each row representing the expression level of a single mRNA in different samples. The color range of the heat map corresponds from blue (low expression) to red (high expression).
Preprints 84700 g004
Figure 5. Weighted co-expression network analysis of all genes. (A) Gene hierarchical clustering partition module. (B) Histogram showing the number of genes in each module. (C) Correlation analysis between modules and phenotypes. (D) Correlation analysis between modules and groups. LW and MS represent Large White pigs and Mashen pigs, respectively.
Figure 5. Weighted co-expression network analysis of all genes. (A) Gene hierarchical clustering partition module. (B) Histogram showing the number of genes in each module. (C) Correlation analysis between modules and phenotypes. (D) Correlation analysis between modules and groups. LW and MS represent Large White pigs and Mashen pigs, respectively.
Preprints 84700 g005
Figure 6. The expression patterns and functional analysis of genes in the co-expression modules. (A) The expression pattern of genes in the pink module. (B) Gene Ontology (GO) analysis of genes in the pink module. (C) Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of genes in the pink module. (D) The expression pattern of genes in the blue2 module. (E) GO analysis of genes in the blue2 module. (F) KEGG enrichment analysis of genes in the blue2 module. (G) Combined results of traditional transcriptomics analysis and weighted co-expression network analysis (WGCNA) to screen common and specific candidate genes responding to cold response in LW and MS pigs. The heat map color scale indicates expression levels, ranging from blue (low expression) to red (high expression). LW and MS denote Large White pigs and Mashen pigs, respectively; 25 and 4 represent 25℃ and 4℃, respectively.
Figure 6. The expression patterns and functional analysis of genes in the co-expression modules. (A) The expression pattern of genes in the pink module. (B) Gene Ontology (GO) analysis of genes in the pink module. (C) Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of genes in the pink module. (D) The expression pattern of genes in the blue2 module. (E) GO analysis of genes in the blue2 module. (F) KEGG enrichment analysis of genes in the blue2 module. (G) Combined results of traditional transcriptomics analysis and weighted co-expression network analysis (WGCNA) to screen common and specific candidate genes responding to cold response in LW and MS pigs. The heat map color scale indicates expression levels, ranging from blue (low expression) to red (high expression). LW and MS denote Large White pigs and Mashen pigs, respectively; 25 and 4 represent 25℃ and 4℃, respectively.
Preprints 84700 g006
Figure 7. Validation results of common and specific cold-responsive candidate genes in pig longissimus dorsi muscle tissue and satellite cells. (A) qRT-PCR verification of candidate genes in the longissimus dorsi muscle tissue of LW and MS pigs at 25℃ and 4℃. (B) Western blot verification of candidate genes in the longissimus dorsi muscle tissue of LW and MS pigs at 25℃ and 4℃. (C) qRT-PCR verification of candidate genes in porcine skeletal muscle satellite cells at 37℃ and 32℃. (D) Detection of interference efficiency of PRSS8 in porcine skeletal muscle satellite cells. (E) Effects of PRSS8 interference on mRNA levels of heat production and energy metabolism in porcine skeletal muscle satellite cells. (F) Effects of PRSS8 interference on protein levels of heat production and energy metabolism in porcine skeletal muscle satellite cells. LW represents Large White pig; MS represents Mashen pig. P values were determined using the two-tailed unpaired Student’s t-test. **P<0.01; *P<0.05.
Figure 7. Validation results of common and specific cold-responsive candidate genes in pig longissimus dorsi muscle tissue and satellite cells. (A) qRT-PCR verification of candidate genes in the longissimus dorsi muscle tissue of LW and MS pigs at 25℃ and 4℃. (B) Western blot verification of candidate genes in the longissimus dorsi muscle tissue of LW and MS pigs at 25℃ and 4℃. (C) qRT-PCR verification of candidate genes in porcine skeletal muscle satellite cells at 37℃ and 32℃. (D) Detection of interference efficiency of PRSS8 in porcine skeletal muscle satellite cells. (E) Effects of PRSS8 interference on mRNA levels of heat production and energy metabolism in porcine skeletal muscle satellite cells. (F) Effects of PRSS8 interference on protein levels of heat production and energy metabolism in porcine skeletal muscle satellite cells. LW represents Large White pig; MS represents Mashen pig. P values were determined using the two-tailed unpaired Student’s t-test. **P<0.01; *P<0.05.
Preprints 84700 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated