Transcriptome Profiles of the Liver in Two Cold-Exposed Sheep Breeds Revealed Different Mechanisms and Candidate Genes for Thermogenesis

Jiao, Dan; Ji, Kaixi; Wang, Wenqiang; Liu, Hu; Zhou, Jianwei; Degen, A. Allan; Zhang, Yunsheng; Zhou, Ping; Yang, Guo

doi:https://doi.org/10.1155/2021/5510297

Genetics Research

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Consent Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 5510297 | https://doi.org/10.1155/2021/5510297

Transcriptome Profiles of the Liver in Two Cold-Exposed Sheep Breeds Revealed Different Mechanisms and Candidate Genes for Thermogenesis

Dan Jiao,^1,2,3Kaixi Ji,^1,2,3Wenqiang Wang,⁴Hu Liu,⁴Jianwei Zhou,⁴A. Allan Degen,⁵Yunsheng Zhang,⁶Ping Zhou,⁷and Guo Yang^1,3

Academic Editor: Vindhya Mohindra

Received31 Jan 2021

Accepted29 Jul 2021

Published11 Aug 2021

Abstract

Cold-induced thermogenesis plays an important role in the survival of lambs exposed to low air temperatures. The liver produces and mediates heat production in mammals; however, to date, little is known about the role of liver genes in cold-induced thermogenesis in lambs. In this study, the difference in the liver transcriptome between Altay and Hu ewe lambs was compared. Because of different backgrounds of the two breeds, we hypothesized that the transcriptome profiles of the liver would differ between breeds when exposed to cold. Cold-exposed Altay lambs activated 8 candidate genes (ACTA1, MYH1, MYH2, MYL1, MYL2, TNNC1, TNNC2, and TNNT3) involved in muscle shivering thermogenesis; 3 candidate genes (ATP2A1, SLN, and CKM) involved in muscle nonshivering thermogenesis related to the Ca²⁺ signal and creatine cycle; and 6 candidate genes (PFKM, ALDOC, PGAM2, ENO2, ENO3, and ENO4) involved in enhancing liver metabolism. In contrast, the liver may not act as the main tissue for thermogenesis in cold-exposed Hu lambs. We concluded that Altay lambs rely on liver-mediated shivering and nonshivering thermogenesis by muscle tissue to a greater extent than Hu lambs. Results from this study could provide a theoretical foundation for the breeding and production of cold-resistant lambs.

1. Introduction

Sheep are thermostable animals and rely on thermogenesis to maintain a relatively constant body temperature; however, exposure to extreme cold can cause hypothermia and, ultimately, high mortality in lambs. Prolonged exposure to cold results in the activation of hormonal and metabolic responses in different tissues [1], including nonshivering thermogenesis (NST) in brown adipose tissue (BAT) [2]. Adrenergic activated NST involves uncoupling protein 1 (UCP1) that facilitates proton leakage across the inner mitochondrial membrane, leading to futile cycling of protons and heat generation [3]. It has been estimated that NST contributes approximately 12–20 percent of the increase in total daily heat production, while NST in BAT contributes approximately 2–17 percent of the increase [4], which suggests that other tissues are also involved in the process. Recently, the liver has been reported to play an important role in providing “fuel” for NST by BAT, at least in cold-exposed mice [5].

The liver produces approximately 10–13 percent of the basal heat production of mammals and has been closely linked with thermogenesis [6]. In rats, body temperature and heat production from the liver decreased when NST was inhibited [7]. Bond and Ntambi [8] reported that the liver was related to an increase in adipose tissue and the synthesis and transport of monounsaturated fatty acids in UCP1-deficient mice. Some liver-derived factors, including fibroblast growth factor 21 (FGF21), have been reported to be important in the regulation of thermogenesis during acute cold exposure in mice [9]. However, the genetic mechanisms of the liver during cold exposure in lambs are still unclear. To fill this gap, at least in part, we examined the changes in the expression of liver-derived genes in two sheep breeds exposed to cold. Altay sheep are indigenous in the Altay Prefecture of the northern Xinjiang province and are well adapted to the harsh conditions. This region is extremely cold in winter with an average January temperature of −16.3°C (https://Climate-Data.org) and with a long period of winter snow cover. In contrast, Hu sheep are typically raised in the subtropical climate zone of China, a region characterized by warm, humid conditions. In a companion study on the Altay and Hu sheep, maintenance energy requirements were higher for Hu than for Altay sheep at −5°C, and this temperature was below the thermoneutral zone for Hu sheep but within the thermoneutral zone for Altay sheep [10]. Because of different backgrounds and responses to cold exposure between the two breeds, we hypothesized that the transcriptome profiles of their livers would differ when exposed to cold. Results of this study provide (1) a better understanding of the role of the liver in thermogenic mechanisms employed by lambs and (2) a research base for improving the cold tolerance in sheep.

2. Materials and Methods

2.1. Animals and Ethics

The study design and all procedures on sheep were approved by the Academic Committee of the Northwestern Institute of Eco-Environment Resources, Chinese Academy of Sciences (protocol no. CAS201810082). Twelve Altay (A) and 12 Hu (H) ewe lambs, all six months of age and of similar body weight (29.3 ± 2.47 kg), were purchased from a feedlot in the Altay region and were raised at the Gansu Gaolan Field Scientific Observation and Research Station for Agricultural Ecosystem (36°14′16 N, 103°47′59E). Wool length did not differ between the two breeds and averaged 8.87 ± 2.15 cm. The lambs were fed alfalfa pellets ad libitum, with free access to water, prior to and during the study. They were maintained in individual metabolic cages (1.5 m × 1.0 m), and each breed was divided into two groups: chronic cold-exposed (6 A^c and 6 H^c), which were maintained in a room at a constant temperature of −5°C for 25 days, after a gradual decrease of 2.5°C/day over 10 days, and control (6 A^w and 6 H^w), which were maintained in a room at a constant temperature of 20°C throughout this period. The temperature and humidity of the cold room were −5°C ± 0.03 and 88% ± 6.5, respectively, and of the thermoneutral room were 20°C ± 0.42 and 87.5% ± 9.9, respectively. The temperature humidity indices (THI), following Tucker et al. [11], were 26 ± 0.18 and 63 ± 0.62, respectively. All lambs were slaughtered after 12 h of fasting in the morning after 25 days of temperature exposure, and their livers were excised and immediately frozen in liquid nitrogen and then stored at −80°C.

2.2. RNA Extraction, Library Construction, and Sequencing

Total RNA was extracted from the liver of 19 lambs (A^c (n = 5), A^w (n = 4), H^c (n = 6), and H^w (n = 4)) using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocol. The RNA quality and integrity number (RIN) were measured using Agilent 2100 (Agilent Technologies, Santa Clara, CA, USA). The concentration and purity of the RNA samples were tested through the threshold filter of RIN >7.0 and 28S/18S rRNA ratio >1.0 to ensure that RNA quality meets sequencing standards. The RNA samples from individual lambs in every group (independent biological replicates) were not pooled in order to exclude samples with poor biological duplication to ensure the reliability of all sequencing results and subsequent analyses. Consequently, of 24 liver samples (12 Hu and 12 Altay lambs) sequenced, 19 were used for analysis. Poly A messenger RNA (mRNA) was isolated from total RNA by an oligo dT extraction kit (NEBnext Poly(A) mRNA Magnetic Isolation Module, NEB, USA) and then fragmented using divalent cations under elevated temperature. First-strand cDNA was synthesized from fragmented mRNA using random oligonucleotide primers and reverse transcriptase (SuperScript II Reverse Transcriptase, Invitrogen, Carlsbad, CA, USA) and second strand from DNA polymerase I and RNase H treatments. The cDNA fragments’ production had a single “A” nucleotide base added, followed by ligation of an adapter. The products were purified by AMPure XP beads and then dissolved in EB solution and enriched with PCR amplification to create the final cDNA library. The overall quality of the PCR product was validated by the Agilent Technologies 2100 bioanalyzer. The double-stranded PCR products were heated, denatured, and circularized by the splint oligo sequence to obtain the final library. The cDNA fragments in the library were sequenced using a BGISEQ-500 platform (Beijing Genomics Institute (BGI), Beijing, China) for producing raw reads.

2.3. RNA-Seq Data Analysis

The clean reads were filtered out based on the raw reads, using quality control software SOAPnuke (BGI), which were obtained by removing low-quality reads (more than 20 percent of bases in the total reads had quality scores lower than 15), adaptor reads (reads with joint contamination), and unknown base N content (reads which contain more than 5 percent undetermined base information). The reads were aligned and annotated to the reference genome of Ovis aries (Oar_rambouillet_v1.0; https://www.ncbi.nlm.nih.gov/assembly/GCF_002742125.1) using the HISAT alignment tool (Centre for Computational Biology, Johns Hopkins University, MD, USA). HISAT is based on the Burrows–Wheeler transform and Ferragina–Manzini (FM) indexing methods [12]. We used Bowtie 2 for calculating the gene alignment rate [13] (Johns Hopkins University, MD, USA) and then calculated gene expression levels with RSEM (version 1.2.12, University of Wisconsin–Madison, USA), a software package for estimating gene and isoform expression levels from RNA-Seq data [14]. The gene expression levels were standardized by reads per kilobase per million (FPKM) mapped reads. The constrained principal coordinate analysis (cPCoA) was employed to visualize classical multidimensional scaling of Bray–Curtis distance matrices by using functions capscale and anova.cca of vegan package in R (version 4.0, Lucent Technologies, AZ, USA), and the value was calculated by permutation tests [15, 16].

2.4. Differentially Expressed Gene Analysis

We compared differential gene expressions in the liver between breeds and between air temperatures using pairwise comparisons, as described by Wang et al. [17]. Differentially expressed genes (DEGs) were filtered by DESeq2 software and as fold changes (FC, |log2 FC| > 1) and q value (q < 0.05). The q value is based on the multiple hypothesis testing on the value. Fold change (FC) was calculated as follows:

DEGs meeting the above screening criteria were carried out by subsequent clustering analysis.

2.5. Function Enrichment and Analyses

Gene Ontology (GO) [18] and Kyoto Encyclopedia of Genes and Genomes (KEGG) [19] enrichment analyses were based on DEGs. GO terms were enriched by phyper functions in R software (version 4.0, Lucent Technologies, AZ, USA), and based on the GO annotation results, DEGs were mapped to the GO terms in the database (https://www.geneontology.org/). The KEGG pathway enrichment was used to identify genes and the metabolic pathways involved in the DEGs (https://www.kegg.jp/kegg/pathway.html/), and this method was consistent with the GO enrichment. In order to further study the effect of cold exposure on gene transcription regulation in sheep, we enriched the up- and downregulated genes into GO terms and KEGG pathways, respectively. A level of q value < 0.05 was accepted as a significant difference between means.

2.6. RNA-Seq Validation by Quantitative Real-Time PCR

To validate the gene expression differences that had been identified by the RNA-seq analysis, four candidate genes, namely, apolipoprotein D (APOD), apolipoprotein A4 (APOA4), lipin (LPIN), and trefoil factor 2 (TFF2), were identified using an RT-qPCR approach. These four genes were significantly regulated in lambs by cold exposure, and their expression levels were high; they also play important roles in liver metabolism.

Total RNA samples were obtained from the liver of 19 lambs (A^c (n = 5), A^w (n = 4), H^c (n = 6), and H^w (n = 4)), which were used for transcriptome sequencing. The primers for RT-qPCR were designed using Oligo 7 (Wojciech Rychlik, USA) (Table 1), with β-actin as a reference gene to certify the relative level of expression. The RT-qPCR amplification of cDNA pools used a PrimeScript RT reagent kit (Takara) with gDNA Eraser (Takara) according to the manufacturer’s instructions. The RT-qPCR reactions were performed on Mx3000P and Mx3005P QPCR Systems (Stratagene, Agilent, Santa Clara, CA, USA). The qPCR reaction system was performed in a total volume of 20 μL containing 2 μL of cDNA, 0.8 μL forward and reverse primers (10 μM), 10 μL TB Green^TM Premix Ex Taq II, 6 μL RNase-free water, and 0.4 μL ROX Reference Dye II (50×). The thermal profile for amplification was a two-step approach, which, after predegeneration, consisted of 15 s at 95°C in the first step and then 5 s at 95°C and 34 s at different Tm for 40 cycles in the second step. Changes in gene expressions were determined by the 2^−△△ct method [20]. Relative quantity between treatment and control groups was tested using t-test.

2.7. Statistical Analyses

Data are expressed as means ± SE, and all statistical analyses were performed by t-tests using SPSS software (SPSS version 17.0, Chicago, IL, USA).

3. Results

3.1. Sequencing and Mapping

To examine the global difference between breeds and treatments in the transcriptome sequences, constrained principal coordinate analysis (cPCoA) by Bray–Curtis distances was employed for every biological replicate of each tissue and treatment (Figure S1). It emerged that the 19 samples explained 16.4% of the variance of the total sequencing data, and the two breeds clustered well (). The livers of –5°C and 20°C Hu lambs were distinct in CPCo 1 (explained 47.1% of 16.4% of the variance of the total sequencing data), whereas the liver tissues were distinct between breeds in CPCo 2 (explained 30.4% of 47.1% of the variance in CPCo 1). The results indicated that the transcriptome sequences in liver tissues separated the two breeds and separated the two temperature treatments in the Hu lambs.

The alignments for the genome and gene sequences were all made with Oar_rambouillet_v1.0. The sequencing and mapping data are summarized in Table S1. The results indicated that sequencing quality met the requirements for subsequent analysis.

3.2. Gene Annotation

A total of 27,298 genes were annotated to the Ovis aries reference genome Oar_rambouillet_v1.0, including 24,244 previously identified genes and 3,054 potentially new genes. Moreover, the analyses revealed that 89.5% of the genes were common in the four lamb groups (Figure 1(a)). In addition, the common DEG number was greater in the compared groups of A-liver^c-A-liver^w and A-liver^c-H-liver^c than the compared groups of A-liver^w-H-liver^w and H-liver^c-H-liver^w, which indicated that these two groups had more similar responsive DEGs (Figure 1(b)).

(a)

(b)

(c)

Figure 1

DEGs in the liver of cold-exposed Altay and Hu lambs. A-liver^c-A-liver^w: liver of −5°C Altay lambs (n = 5) compared to 20°C Altay lambs (n = 4); H-liver^c-H-liver^w: liver of −5°C Hu lambs (n = 6) compared to 20°C Hu lambs (n = 4); A-liver^c-H-liver^c: liver of −5°C Altay lambs (n = 5) compared to −5°C Hu lambs (n = 6); A-liver^w-H-liver^w: liver of 20°C Altay lambs (n = 4) compared to 20°C Hu lambs (n = 4). (a) DEGs among different treatment groups. (b) DEGs between differently compared groups. (c) The up- and downregulated DEGs between differently compared groups. The screen of DEGs is based on |log2 FC| > 1 and q value < 0.05. The green columns represent upregulated DEGs, and blue columns represent downregulated DEGs.

3.3. Analysis of the DEGs

After cold exposure, more DEGs were upregulated in the −5°C Altay lambs compared to −5°C Hu lambs, but more DEGs were downregulated in the −5°C Hu lambs compared to 20°C Hu lambs (Figure 1(c)).

According to the classification of GO terms, DEGs were clustered into molecular function (MF), cellular component (CC), and biological process (BP). The selected significant GO terms’ annotation is presented in Table S2, and upregulated GO terms are presented in Figure 2. The CC GO category in the cold-exposed Altay lambs was related to muscle contraction, MF GO category was related to the regulation of the binding of muscle-related proteins, and BP GO category was related to muscle contraction, but GO terms in the cold-exposed Hu and Altay lambs differed and were related to the hemoglobin complex in the CC GO category, oxygen carrier activity in the MF GO category, and oxygen transport in the BP GO category. For the downregulated GO terms, hormone activity and lipid binding were enriched in the −5°C Altay lambs compared to 20°C Altay lambs, but chemokine activity and Ca²⁺ binding were enriched in the −5°C Hu lambs compared to 20°C Hu lambs.

Figure 2

The significant upregulated GO terms and KEGG pathways of DEGs in the liver. A-liver^c-A-liver^w: liver of −5°C Altay lambs (n = 5) compared to 20°C Altay lambs (n = 4); H-liver^c-H-liver^w: liver of −5°C Hu lambs (n = 6) compared to 20°C Hu lambs (n = 4); A-liver^c-H-liver^c: liver of −5°C Altay lambs compared to −5°C Hu lambs; A-liver^w-H-liver^w: liver of 20°C Altay lambs compared to 20°C Hu lambs. The screen of significant enrichment is based on q value < 0.05. The green columns represent the group of A-liver^w-A-liver^c, the red columns represent the group of H-liver^w-H-liver^c, the khaki columns represent the group of H-liver^c-A-liver^c, and the gray columns represent the group of H-liver^w-A-liver^w.

Upregulated KEGG pathways in the −5°C Altay lambs compared to 20°C Altay lambs were enriched in cardiac muscle contraction, methane metabolism, and adrenergic signaling in cardiomyocyte pathways. Cardiac muscle contraction and methane metabolism pathways were also upregulated in the −5°C Altay lambs compared to −5°C Hu lambs, but there was no difference in the metabolism pathway in the 20°C Altay lambs compared to 20°C Hu lambs. In contrast to upregulated KEGG pathways, the immune, disease, and metabolism-related pathways were downregulated in the 20°C Altay lambs compared to 20°C Hu lambs. The selected significant KEGG pathways are presented in Table S3.

The heatmap of the top 50 DEGs of the liver exhibited reversed expression trends in the two breeds under cold exposure (Figure 3(a)). The top 50 DEGs of the liver at different temperatures in Altay and Hu lambs are presented in Table S4. The DEGs in the −5°C Hu lambs compared to 20°C Hu lambs and the two breeds at 20°C were primarily downregulated, but DEGs in the −5°C Altay lambs compared to 20°C Altay lambs and compared to −5°C Hu lambs were primarily upregulated. These DEGs were related to energy metabolism, muscle development, and Ca²⁺ binding, whereas the downregulated DEGs in the −5°C Altay lambs compared to 20°C Altay lambs were related to lipid metabolism. We clustered the DEGs, which were enriched in the significant GO and KEGG pathways, including muscle contraction, methane metabolism, lipid metabolism, and oxygen transport (Figure 3(b)). The candidate genes related to muscle contraction, methane metabolism, and oxygen transport were upregulated significantly in the cold-exposed Altay lambs, and almost all of the candidate genes were downregulated in the cold-exposed Hu lambs. Lipid metabolism-related genes displayed large differences between breeds.

(a)

(b)

Figure 3

The heatmap of DEGs in the liver. A-liver^c-A-liver^w: liver of −5°C Altay lambs (n = 5) compared to 20°C Altay lambs (n = 4); H-liver^c-H-liver^w: liver of −5°C Hu lambs (n = 6) compared to 20°C Hu lambs (n = 4); A-liver^c-H-liver^c: liver of −5°C Altay lambs compared to −5°C Hu lambs; A-liver^w-H-liver^w: liver of 20°C Altay lambs compared to 20°C Hu lambs. (a) A heatmap of top 50 DEGs. (b) A heatmap of candidate DEGs. The cluster analysis of gene expression is based on log2 FPKM data. The red color represents higher expression, and the green color represents lower expression.

3.4. Validation of RNA-Seq Results by RT-qPCR

Four selected DEGs, namely, APOD, APOA4, LPIN, and TFF2, were used to verify RNA-seq by RT-qPCR, and all were in agreement. The four genes were chosen for RT-qPCR because the FPKM of these genes was relatively high, and the genes related to lipid and energy metabolism displayed a high correlation with our study. The relative expression levels of these genes from liver tissue in all four groups are presented in Figure 4.

Figure 4

Expression levels of four candidate genes from RT-PCR and RNA-seq. A-liver^w and A-liver^c: liver of Altay lambs under thermoneutral (20°C; n = 4) and cold exposure (−5°C; n = 5); H-liver^w and H-liver^c: liver of Hu lambs under thermoneutral (20°C; n = 4) and cold exposure (−5°C) (n = 6). The X-axis represents the tissue of the liver; the Y-axis on the left represents the relative gene FPKM levels of RNA-seq by columns and bars; the Y-axis on the right represents the relative gene expression levels of RT-PCR by columns and bars. ^∗P < 0.05.

4. Discussion

A number of studies have reported that the liver has a key function in NST. For example, it was reported that the liver provided approximately 44% of the total metabolic energy with cold acclimation in short-tailed opossums [21], oxidation capacity of the liver of ducks increased by 40% after cold adaptation, and with cold adaptation, the liver provided BAT with glucose and fatty acids from very-low-density lipoproteins (VLDLs), which contributed to heat production [22]. In a recent study on the transcriptome profiles of cold-adapted Mongolian sheep [23], cold exposure induced postponing cell senescence in the liver, but no direct evidence of liver involvement in thermogenesis was reported. Therefore, to date, little is known about transcriptome profiles of the liver of lambs when exposed to cold as most studies have been done on rodents and humans. Molecular genetic studies can provide new insights in understanding the mechanisms underlying the tolerance of lambs to cold exposure.

4.1. Altay and Hu Sheep Displayed Breed Differences

A large number of downregulated DEGs were enriched in the first level of KEGG pathways, including diseases, organismal systems, and metabolism pathways in Altay lambs compared to Hu lambs. These genetic differences between sheep breeds occurred, especially in immune responses. The DEGs CD40 and CD40 ligand (CD40L) were downregulated in Altay compared to Hu lambs and were enriched in all the top 5 KEGG pathways related to immune processes. CD40 and CD40L regulate inflammatory processes through secondary messengers [24]. By regulating the transcription of different downstream factors, they mediate the immune signal transduction due to Leishmania infection [25] and oxidation stress [26]. Enrichment of downregulated genes in immune-related pathways indicated that Altay lambs have stronger immune resistance than Hu lambs. This could be explained by their different backgrounds and adaptations as Altay sheep are indigenous to Xinjiang Altay, have adapted well to the long, cold winters and sparse forage of low protein content, and are resistant to diseases, whereas Hu sheep are typically raised in warm, humid areas and are bred for high reproductive rates.

4.2. Cold-Exposed Altay Lambs Enhance Muscle ST and NST through Liver Regulation

The GO analysis indicated that cold-exposed Altay lambs mobilized several terms in the liver related to muscle regulation, such as muscle contraction, transition between fast and slow fibers, and cardiac muscle contraction. These terms were all upregulated. In addition, the cardiac muscle contraction pathway was enriched, and the methane metabolism pathway was upregulated. A large number of DEGs that were focused on muscle contraction, including actin alpha 1 (ACTA1), myosin light chain, phosphorylatable, fast skeletal muscle (MYLPF), myosin heavy chain 1 (MYH1), myosin heavy chain 2 (MYH2), myosin light chain 1 (MYL1), myosin light chain 2 (MYL2), troponin C1, slow skeletal and cardiac type (TNNC1), troponin C2, fast skeletal type (TNNC2), and troponin T3, fast skeletal type (TNNT3), were upregulated in the top 50 DEGs in cold-exposed Altay lambs. These highly expressed genes in the liver of cold-exposed Altay lambs responded to cold stimuli, which indicated that cold exposure enhanced liver-mediated muscle metabolism, including muscle contraction and transition between fast and slow fibers. Studies have reported that the liver has a close metabolic relationship with the muscle and that it could regulate heat production by stimulating muscular activity [27, 28]. In the early stage of cold exposure, ST, through skeletal muscle contraction, plays a major role in thermoregulation, but this is temporary [29]. Recent studies have demonstrated that the skeletal muscle is not only involved in ST but also in NST [30, 31]. In the present study, sarco-endoplasmic reticulum Ca²⁺-ATPase 1 (ATP2A1, SERCA1) and sarcolipin (SLN) were upregulated in the cold-exposed Altay lambs, and the GO terms of the Ca²⁺ transmembrane transport and sarcoplasmic reticulum were activated. The skeletal muscle activates SLN by regulating Ca²⁺-ATPase in the endoplasmic reticulum (ER) to produce heat by muscular NST [28, 31, 32]. ATP2A1 (SERCA1), a key enzyme in pumping Ca²⁺, converts part of the energy into heat [33, 34], and ATP2A1 binds with SLN, in the presence of Ca²⁺, and promotes uncoupling of the ATP2A pump [35, 36]. There is evidence that all energy from uncoupled ATP hydrolysis is converted into heat [37], which would suggest that the liver in Altay lambs regulates muscular ST and NST in response to cold exposure. In addition, in the current study, creatine kinase (CKM) was the top upregulated DEG in cold-exposed Altay lambs and was enriched in the GO term in response to heat. CKM triggers the phosphorylation of creatine from mitochondrial ATP, which liberates local ADP, and then, phosphocreatine (Pcr) is hydrolyzed to creatine to initiate another round of the futile creatine cycle [38]. The Pcr circuit plays an important role in cells that are both excitable and require a high thermodynamic efficiency, such as myocytes. The circuit can transmit chemical energy to the required sites of high ATP demand through cellular distribution of creatine kinase [39, 40].

4.3. Cold-Exposed Altay Lambs Enhance Liver Metabolism

Interestingly, the methane metabolism pathway was upregulated in cold-exposed Altay lambs. Related DEGs, including phosphofructokinase, muscle (PFKM), aldolase, fructose-bisphosphate C (ALDOC), phosphoglycerate mutase 2 (PGAM2), enolase 2 (ENO2), enolase 3 (ENO3), and enolase 4 (ENO4), were also enriched and involved in the GO term of glycolytic processes. PFKM is the limiting enzyme of the glycolytic process and represents a major control point in the metabolism of glucose [41]. The function of PFKM was reported to be associated with the growth of methane-producing bacteria. ENOs catalyze the interconversion of 2-phosphoglycerate to phosphoenolpyruvate, which is the second of the two high-energy intermediates that generate ATP in glycolysis [42, 43]. Consequently, cold exposure could result in an increase in glucose metabolism and methane production in Altay lambs. However, cold exposure did not result in changes in rumen methane transcriptome profiles in Altay lambs (unpublished data). Studies have shown that high temperature could reduce methane production in sheep [44], but whether low temperature increases methane emission needs further research.

We found that the GO terms of oxygen transport, oxygen carrier activity, oxygen binding, and KEGG pathway of cardiac muscle contraction were enriched in Altay lambs. Enzyme activity in the citric acid cycle (TCA cycle) was reported to be directly proportional to the rate of oxygen utilization by muscle cells as a function of the increase in ATP requirement for myocardial contraction [45]. Apparently, these 6 candidate genes are related to the glycolytic pathway and TCA cycle to produce energy for thermogenesis. Several DEGs involved in lipid metabolism, such as APOD and endothelial lipase (LIPG), were upregulated in the cold-exposed Altay lambs compared to 20°C Altay lambs. LIPG clears the circulation of high-density lipoproteins and provides lipid precursors for lipid synthesis [46, 47], while the APOA4 gene is a major constituent of high-density lipoprotein particles [48] and is downregulated in the cold-exposed Altay lambs compared to 20°C Altay lambs. This indicated that cold exposure reduced the circulation of high-density lipoproteins, which resulted in a reduction in transport of cholesterol to the liver. Furthermore, APOD not only influences lipid metabolism but also plays an important role in glucose homeostasis [49, 50]. In addition, TFF2, as a regulator of energy metabolism [51], was upregulated in the cold-exposed Altay lambs compared to −5°C Hu lambs. These results indicated that cold exposure influenced the cholesterol transport in the liver of Altay lambs and enhanced glycolysis in the liver to generate heat.

4.4. The Liver May Not Be the Main Thermogenic Tissue for Hu Lambs

Through the cluster of candidate DEGs, almost all candidate genes were downregulated in the cold-exposed Hu lambs, which means that the liver may not be the main thermogenic tissue. Compared to cold-exposed Altay lambs, the response of the liver in Hu lambs to cold exposure was focused mainly on immune diseases and immune system pathways. The gene transcriptional expression level had smaller differences between breeds of 20°C lamb group. In a previous study (unpublished data), cold-exposed Hu lambs relied mainly on tail fat for heat production, rather than typical thermogenic tissues such as the liver and muscles. In the present study, cellular retinoic acid-binding protein type 1 (CRABP1) was upregulated in the −5°C Hu lambs compared to 20°C Hu lambs, but downregulated in the −5°C Altay lambs compared to −5°C Hu lambs, whereas acyl-coA thioesterase 11 (ACOT11) regulation displayed an opposite trend. CRABP1 is one of the carrier proteins of retinoic acid and was reported to be associated with fat accumulation in adipose tissue in mice [52]. ACOT11 was reported to limit the oxidation of lipid droplet-derived fatty acids, possibly by regulating the availability of substrates to be used for β-oxidation and uncoupling [53–55]. These results indicate that cold exposure mobilizes adipose tissue as the main heat-producing tissue in Hu lambs and could explain why cold exposure did not cause enrichment and upregulation of energy metabolism-related pathways and genes in their liver.

5. Conclusions

Transcriptome sequencing of the liver of cold-exposed Altay and Hu lambs revealed different thermogenic pathways between breeds. Cold exposure induced thermogenesis in Altay lambs and activated the liver-regulated muscle contraction pathway related to ST of muscles. Eight candidate genes, namely, ACTA1, MYH1, MYH2, MYL1, MYL2, TNNC1, TNNC2, and TNNT3, were upregulated in response to cold. The Ca²⁺ signal pathway and creatine cycle were also activated, and 3 candidate genes, including ATP2A1, SLN, and CKM, were involved and upregulated in muscular NST. In addition, 6 candidate genes related to methane metabolism, PFKM, ALDOC, PGAM2, ENO2, ENO3, and ENO4, were upregulated in the liver of cold-exposed Altay lambs. However, in the cold-exposed Hu lambs, it appears that the liver is not the main tissue for thermogenesis, but has a much lesser role when compared to Altay lambs. In summary, Altay and Hu lambs have different regulation mechanisms in response to cold exposure, which are likely manifested in breed differences in cold resistance.

Data Availability

The data are available confidentially to editors and reviewers, and all transcriptome data were submitted to the NCBI Sequence Read Archive (SRA series accession: PRJNA639638).

Not applicable.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

DJ summarized the transcriptome sequencing data, performed the data analysis, and prepared the original manuscript. KXJ, WQW, and HL attended the animals and collected samples. JWZ, YSZ, PZ, and GY designed the experiments and provided the research platform. AAD and GY revised the manuscript. All authors approved the final manuscript.

Acknowledgments

The authors are grateful to the technicians who took care of the animals in this study. They also thank colleagues in the Gaolan Agri-Ecology Research Station for their assistance. This study was funded by the Hundred Talents Program of the Chinese Academy of Sciences (no. Y629721002) and study on the haplotype genetic expression of the ADRB3-energy metabolism gene group in the Xinjiang Sheep Breeding Population (no. 2019CA009).

Supplementary Materials

Supplementary Material 1: Figure S1: CPCoA analyses of all samples and sequencing quality in the liver of Altay and Hu lambs. Supplementary Material 2: Table S1: summary of RNA-seq results. Supplementary Material 3: Table S2: GO terms significantly enriched in the liver at different temperatures in Altay and Hu lambs. Supplementary Material 4: Table S3: KEGG pathways significantly enriched in the liver at different temperatures in Altay and Hu lambs. Supplementary Material 5: Table S4: top 50 DEGs in the liver at different temperatures in Altay and Hu lambs. s (Supplementary Materials)

References

T. Pääkkönen and J. Leppäluoto, “Cold exposure and hormonal secretion: a review,” International Journal of Circumpolar Health, vol. 61, no. 3, pp. 265–276, 2002.
View at: Publisher Site | Google Scholar
G. Alexander and D. Williams, “Shivering and non-shivering thermogenesis during summit metabolism in young lambs,” The Journal of Physiology, vol. 198, no. 2, pp. 251–276, 1968.
View at: Publisher Site | Google Scholar
N. A. Abumrad, “The liver as a hub in thermogenesis,” Cell Metabolism, vol. 26, no. 3, pp. 454-455, 2017.
View at: Publisher Site | Google Scholar
C. Y. Tan, K. Ishikawa, S. Virtue, and A. Vidal-Puig, “Brown adipose tissue in the treatment of obesity and diabetes: are we hot enough?” Journal of Diabetes Investigation, vol. 2, no. 5, pp. 341–350, 2011.
View at: Publisher Site | Google Scholar
R. Yao, Y. Yang, S. Lian et al., “Effects of acute cold stress on liver o-glcnacylation and glycometabolism in mice,” International Journal of Molecular Sciences, vol. 19, no. 9, pp. 2815–2829, 2018.
View at: Publisher Site | Google Scholar
R. W. Brauer, R. J. Holloway, J. S. Krebs, G. F. Leong, and H. W. Carroll, “The liver in hypothermia,” Annals of the New York Academy of Sciences, vol. 01, pp. 395–423, 1959.
View at: Google Scholar
H. B. Stoner, “The role of the liver in non-shivering thermogenesis in the rat,” The Journal of Physiology, vol. 232, no. 2, pp. 285–296, 1973.
View at: Publisher Site | Google Scholar
L. M. Bond and J. M. Ntambi, “UCP1 deficiency increases adipose tissue monounsaturated fatty acid synthesis and trafficking to the liver,” Journal of Lipid Research, vol. 59, no. 2, pp. 224–236, 2018.
View at: Publisher Site | Google Scholar
M. Ameka, K. R. Markan, D. A. Morgan et al., “Liver derived FGF21 maintains core body temperature during acute cold exposure,” Scientific Reports, vol. 9, pp. 1–10, 2019.
View at: Publisher Site | Google Scholar
J. Zhou, K. Ji, H. Liu et al., “Effect of air temperature on growth performance, apparent digestibilities, rumen fermentation and serum metabolites in Altay and Hu lambs,” Journal of Animal Physiology and Animal Nutrition, vol. 104, no. 4, pp. 1024–1033, 2020.
View at: Publisher Site | Google Scholar
C. B. Tucker, A. R. Rogers, and K. E. Schütz, “Effect of solar radiation on dairy cattle behaviour, use of shade and body temperature in a pasture-based system,” Applied Animal Behaviour Science, vol. 109, no. 2-4, pp. 141–154, 2008.
View at: Publisher Site | Google Scholar
D. Kim, B. Langmead, and S. L. Salzberg, “HISAT: a fast spliced aligner with low memory requirements,” Nature Methods, vol. 12, no. 4, pp. 357–360, 2015.
View at: Publisher Site | Google Scholar
B. Langmead and S. L. Salzberg, “Fast gapped-read alignment with Bowtie 2,” Nature Methods, vol. 9, no. 4, pp. 357–359, 2012.
View at: Publisher Site | Google Scholar
W. D. van Marken Lichtenbelt and P. Schrauwen, “Implications of nonshivering thermogenesis for energy balance regulation in humans,” American Journal of Physiology-Regulatory, Integrative and Comparative Physiology, vol. 301, no. 2, pp. R285–R296, 2011.
View at: Publisher Site | Google Scholar
A. Buja and N. Eyuboglu, “Remarks on parallel analysis,” Multivariate Behavioral Research, vol. 27, no. 4, pp. 509–540, 1992.
View at: Publisher Site | Google Scholar
M. Linting, B. J. Van Os, and J. J. Meulman, “Statistical significance of the contribution of variables to the PCA solution: an alternative permutation strategy,” Psychometrika, vol. 76, no. 3, pp. 440–460, 2011.
View at: Publisher Site | Google Scholar
L. Wang, Z. Feng, X. Wang, X. Wang, and X. Zhang, “DEGseq: an R package for identifying differentially expressed genes from RNA-seq data,” Bioinformatics, vol. 26, no. 1, pp. 136–138, 2010.
View at: Publisher Site | Google Scholar
M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene Ontology: tool for the unification of biology,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000.
View at: Publisher Site | Google Scholar
H. Ogata, S. Goto, K. Sato, W. Fujibuchi, H. Bono, and M. Kanehisa, “KEGG: kyoto encyclopedia of genes and genomes,” Nucleic Acids Research, vol. 27, pp. 29–34, 2000.
View at: Google Scholar
K. J. Livak and T. D. Schmittgen, “Analysis of relative gene expression data using real-time quantitative PCR and the 2^−ΔΔCT method,” Methods, vol. 25, no. 4, pp. 402–408, 2001.
View at: Publisher Site | Google Scholar
J. J. Villarin, P. J. Schaeffer, R. A. Markle, and S. L. Lindstedt, “Chronic cold exposure increases liver oxidative capacity in the marsupial Monodelphis domestica,” Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology, vol. 136, no. 3, pp. 621–630, 2003.
View at: Publisher Site | Google Scholar
P. Trayhurn, “Fatty acid synthesis in vivo in brown adipose tissue, liver and white adipose tissue of the cold-acclimated rat,” FEBS Letters, vol. 104, no. 1, pp. 13–16, 1979.
View at: Publisher Site | Google Scholar
Y. Cao, J. Pan, Y. Zhang, H. Zhou, and H. Zhou, “Regulation mechanisms of positive cold tolerance in Mongolian sheep,” Animal Biology, vol. 70, no. 1, pp. 27–37, 2020.
View at: Publisher Site | Google Scholar
M. Rizvi, D. Pathak, J. E. Freedman, and S. Chakrabarti, “CD40-CD40 ligand interactions in oxidative stress, inflammation and vascular disease,” Trends in Molecular Medicine, vol. 14, no. 12, pp. 530–538, 2008.
View at: Publisher Site | Google Scholar
R. K. Mathur, A. Awasthi, P. Wadhone, B. Ramanamurthy, and B. Saha, “Reciprocal CD40 signals through p38MAPK and ERK-1/2 induce counteracting immune responses,” Nature Medicine, vol. 10, no. 5, pp. 540–544, 2004.
View at: Publisher Site | Google Scholar
C. Chen, H. Chai, X. Wang et al., “Soluble CD40 ligand induces endothelial dysfunction in human and porcine coronary artery endothelial cells,” Blood, vol. 112, no. 8, pp. 3205–3216, 2008.
View at: Publisher Site | Google Scholar
J.-S. Liu, Y.-Q. Chen, and M. Li, “Thyroid hormones increase liver and muscle thermogenic capacity in the little buntings (Emberiza pusilla),” Journal of Thermal Biology, vol. 31, no. 5, pp. 386–393, 2006.
View at: Publisher Site | Google Scholar
M. Pant, N. C. Bal, and M. Periasamy, “Sarcolipin: a key thermogenic and metabolic regulator in skeletal muscle,” Trends in Endocrinology & Metabolism, vol. 27, no. 12, pp. 881–892, 2016.
View at: Publisher Site | Google Scholar
J. Aydin, I. G. Shabalina, N. Place et al., “Nonshivering thermogenesis protects against defective calcium handling in muscle,” The FASEB Journal, vol. 22, no. 11, pp. 3919–3924, 2008.
View at: Publisher Site | Google Scholar
F. Haman and D. P. Blondin, “Shivering thermogenesis in humans: origin, contribution and metabolic requirement,” Temperature, vol. 4, no. 3, pp. 217–226, 2017.
View at: Publisher Site | Google Scholar
N. C. Bal and M. Periasamy, “Uncoupling of sarcoendoplasmic reticulum calcium ATPase pump activity by sarcolipin as the basis for muscle non-shivering thermogenesis,” Philosophical Transactions of the Royal Society B: Biological Sciences, vol. 375, no. 1793, Article ID 20190135, 2020.
View at: Publisher Site | Google Scholar
N. C. Bal, S. K. Sahoo, S. K. Maurya, and M. Periasamy, “The role of sarcolipin in muscle non-shivering thermogenesis,” Frontiers in Physiology, vol. 9, pp. 1–6, 2018.
View at: Publisher Site | Google Scholar
W. S. D.-S. Ana Paula Arruda, D. P. Carvalho, and L. De Meis, “Hyperthyroidism increases the uncoupled ATPase activity and heat production by the sarcoplasmic reticulum Ca²⁺-ATPase,” Journal of Biological Chemistry, vol. 276, pp. 25078–25087, 2003.
View at: Google Scholar
L. Meis, “Ca²⁺-ATPases (SERCA): energy transduction and heat production in transport ATPases,” Journal of Membrane Biology, vol. 188, no. 1, pp. 1–9, 2002.
View at: Publisher Site | Google Scholar
S. Mall, R. Broadbridge, S. L. Harrison, M. G. Gore, A. G. Lee, and J. M. East, “The presence of sarcolipin results in increased heat production by Ca²⁺-ATPase,” Journal of Biological Chemistry, vol. 281, no. 48, pp. 36597–36602, 2006.
View at: Publisher Site | Google Scholar
N. C. Bal, S. K. Maurya, D. H. Sopariwala et al., “Sarcolipin is a newly identified regulator of muscle-based thermogenesis in mammals,” Nature Medicine, vol. 18, no. 10, pp. 1575–1579, 2012.
View at: Publisher Site | Google Scholar
A. P. Arruda, L. A. Ketzer, M. Nigro, A. Galina, D. P. Carvalho, and L. de Meis, “Cold tolerance in hypothyroid rabbits: role of skeletal muscle mitochondria and sarcoplasmic reticulum Ca²⁺ ATPase isoform 1 heat production,” Endocrinology, vol. 149, no. 12, pp. 6262–6271, 2008.
View at: Publisher Site | Google Scholar
L. Kazak, E. T. Chouchani, M. P. Jedrychowski et al., “A creatine-driven substrate cycle enhances energy expenditure and thermogenesis in beige fat,” Cell, vol. 163, no. 3, pp. 643–655, 2015.
View at: Publisher Site | Google Scholar
K. Steeghs, A. Benders, F. Oerlemans et al., “Altered Ca²⁺ responses in muscles with combined mitochondrial and cytosolic creatine kinase deficiencies,” Cell, vol. 89, no. 1, pp. 93–103, 1997.
View at: Publisher Site | Google Scholar
R. Guzun, N. Timohhina, K. Tepp et al., “Systems bioenergetics of creatine kinase networks: physiological roles of creatine and phosphocreatine in regulation of cardiac cell function,” Amino Acids, vol. 40, no. 5, pp. 1333–1348, 2011.
View at: Publisher Site | Google Scholar
K. Wang, X. Wu, Y. Zhuang et al., “Celastrol directly inhibits PFKM to induce weight loss and leptin sensitization,” bioRxiv, vol. 2020, 2020.
View at: Publisher Site | Google Scholar
G. P. Comi, F. Fortunato, S. Lucchiari et al., “?-enolase deficiency, a new metabolic myopathy of distal glycolysis,” Annals of Neurology, vol. 50, no. 2, pp. 202–207, 2001.
View at: Publisher Site | Google Scholar
F. Wold, “18 enolase,” in The Enzymes, vol. 5, pp. 499–538, Academic Press, Cambridge, MA, USA, 1971.
View at: Google Scholar
V. Sejian, D. Kumar, and S. M. K. Naqvi, “Physiological rhythmicity in Malpura ewes to adapt to cold stress in a semi-arid tropical environment,” Biological Rhythm Research, vol. 49, no. 2, pp. 215–225, 2017.
View at: Publisher Site | Google Scholar
J. Bangsbo, “Muscle oxygen uptake in humans at onset of and during intense exercise,” Acta Physiologica Scandinavica, vol. 168, no. 4, pp. 457–464, 2000.
View at: Publisher Site | Google Scholar
N. Pirany, A. Bakrani Balani, H. Hassanpour, and H. Mehraban, “Differential expression of genes implicated in liver lipid metabolism in broiler chickens differing in weight,” British Poultry Science, vol. 61, no. 1, pp. 10–16, 2020.
View at: Publisher Site | Google Scholar
J. E. Yu, S. Y. Han, B. Wolfson, Q. Zhou, and Q. Zhou, “The role of endothelial lipase in lipid metabolism, inflammation, and cancer,” Histology and Histopathology, vol. 33, pp. 1–10, 2018.
View at: Publisher Site | Google Scholar
M. Świątkiewicz, M. Oczkowicz, K. Ropka-Molik, and E. Hanczakowska, “The effect of dietary fatty acid composition on adipose tissue quality and expression of genes related to lipid metabolism in porcine livers,” Animal Feed Science and Technology, vol. 216, pp. 204–215, 2016.
View at: Google Scholar
G. Perdomo, S. Qu, and H. Dong, “APOD is a novel apolipoprotein that affects glucose metabolism,” Diabetes, vol. 56, p. A367, 2007.
View at: Google Scholar
G. Perdomo, D. H. Kim, T. Zhang et al., “A role of apolipoprotein D in triglyceride metabolism,” Journal of Lipid Research, vol. 51, no. 6, pp. 1298–1311, 2010.
View at: Publisher Site | Google Scholar
A. Ghanemi, A. Melouane, O. Mucunguzi, M. Yoshioka, and J. St-Amand, “Energy and metabolic pathways in trefoil factor family member 2 (Tff2) KO mice beyond the protection from high-fat diet-induced obesity,” Life Sciences, vol. 215, pp. 190–197, 2018.
View at: Publisher Site | Google Scholar
J. E. Miller, Investigating the Role of Crabp1 in Adipose Biology, Master, Case Western Reserve, Cleveland, OH, USA, 2017.
S. Han and D. E. Cohen, “Functional characterization of thioesterase superfamily member 1/Acyl-CoA thioesterase 11: implications for metabolic regulation,” Journal of Lipid Research, vol. 53, no. 12, pp. 2620–2631, 2012.
View at: Publisher Site | Google Scholar
K. Okada, K. B. LeClair, Y. Zhang et al., “Thioesterase superfamily member 1 suppresses cold thermogenesis by limiting the oxidation of lipid droplet-derived fatty acids in brown adipose tissue,” Molecular Metabolism, vol. 5, no. 5, pp. 340–351, 2016.
View at: Publisher Site | Google Scholar
Y. Zhang, Y. Li, M. W. Niepel et al., “Targeted deletion of thioesterase superfamily member 1 promotes energy expenditure and protects against obesity and insulin resistance,” Proceedings of the National Academy of Sciences, vol. 109, no. 14, pp. 5417–5422, 2012.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Dan Jiao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

666

Downloads

1158

Citations