Assessment of cold exposure-induced metabolic changes in mice using untargeted metabolomics

Abstract

Background: Cold exposure (CE) can effectively modulate adipose tissue metabolism and improve metabolic health. Although previous metabolomics studies have primarily focused on analyzing one or two samples from serum, brown adipose tissue (BAT), white adipose tissue (WAT), and liver samples, there is a significant lack of simultaneous analysis of multiple tissues regarding the metabolic changes induced by CE in mice. Therefore, our study aims to investigate the metabolic profiles of the major tissues involved.

Methods: A total of 14 male C57BL/6J mice were randomly assigned to two groups: the control group (n = 7) and the CE group (n = 7). Metabolite determination was carried out using gas chromatography-mass spectrometry (GC-MS), and multivariate analysis was employed to identify metabolites exhibiting differential expression between the two groups.

Results: In our study, we identified 32 discriminant metabolites in BAT, 17 in WAT, 21 in serum, 7 in the liver, 16 in the spleen, and 26 in the kidney, respectively. Among these metabolites, amino acids, fatty acids, and nucleotides emerged as the most significantly altered compounds. These metabolites were found to be associated with 12 differential metabolic pathways closely related to amino acids, fatty acids, and energy metabolism.

Conclusion: Our study may provide valuable insights into the metabolic effects induced by CE, and they have the potential to inspire novel approaches for treating metabolic diseases.

1 Introduction

There is a well-established relationship between ambient temperature and mortality rates (Vialard and Olivier, 2020; Burkart et al., 2021; Fatima et al., 2021). Cold temperatures, as a frequent fluctuation in ambient temperature, have a profound impact on human health. In cold environments, endothermic organisms rely on thermogenesis to maintain their core body temperature. This allows cells to carry out crucial physiological processes and functions (Nguyen et al., 2011).

Brown adipose tissue (BAT) is a unique type of tissue dedicated to the process of non-shivering thermogenesis (NST), mediated by uncoupling protein 1 (UCP1) (Nicholls et al., 1978; Ricquier, 2017). The activation of BAT can be induced by exposure to cold temperatures. Earlier rodent studies have indicated that age significantly influences heat production capabilities in response to cold, with older mice exhibiting a reduced thermogenic effect (Talan et al., 1985; Tatelman and Talan, 1990). Initially, researchers believed that BAT had limited thermogenic and metabolic functions in adult humans, due to the observation that BAT depots were more prevalent in infants (Lean, 1989). However, this perspective has recently been questioned as multiple independent studies have reported the presence of metabolically active BAT in healthy adult humans (Cypess et al., 2009; van Marken Lichtenbelt et al., 2009; Virtanen et al., 2009). Studies conducted on rodents and humans have demonstrated that cold-induced activation of BAT has several beneficial effects on metabolic health. These effects include enhancing glucose uptake (Virtanen et al., 2009; Orava et al., 2011), improving insulin sensitivity (Stanford et al., 2013; Hanssen et al., 2015), stimulating lipolysis (Zechner et al., 2012), reducing circulating levels of triglyceride (TAG) and cholesterol (Bartelt et al., 2011; Berbée et al., 2015; Iwen et al., 2017; Worthmann et al., 2017), as well as clearing circulating branched-chain amino acids (BCAAs) (Yoneshiro et al., 2019). Exposure to cold leads to a swift increase in fibroblast growth factor 21 (FGF21) expression in BAT (Chartoumpekis et al., 2011; Hondares et al., 2011; Fisher and Maratos-Flier, 2016). FGF21 plays a crucial role in promoting the expression of thermogenic genes, including UCP1, within BAT (Fisher et al., 2012). Acting as a systematic peptide hormone, FGF21 plays a significant role in regulating energy balance, as well as maintaining glucose and lipid homeostasis. Additionally, there is a strong correlation between FGF21 levels and decreased circulating levels of BCAAs (Jiang et al., 2015; Karusheva et al., 2019; Yu et al., 2021; Shah et al., 2023). These findings have reignited interest in increasing energy expenditure by CE and BAT activation to combat obesity and its associated metabolic complications, including diabetes, dyslipidemia, and cardiovascular diseases in adult humans (Kajimura and Saito, 2014; Becher et al., 2021).

However, several studies have indicated that the activation of BAT through exposure to cold may have negative effects on health, such as an elevated heart rate, increased blood pressure (Wang et al., 2022), and potential involvement in the progression of breast cancer (Singh et al., 2016). Given these concerns, it is essential to meticulously assess the implications of this approach. Consequently, it is crucial to thoroughly investigate the impact of CE on the functioning of various vital organs in the body. Currently, there is a notable absence of simultaneous analysis of multiple tissues regarding the metabolic changes induced by CE in mice.

Metabolomics, which involves the examination of small molecule metabolites in biological samples, offers a comprehensive understanding of samples and valuable insights into biological alterations resulting from disease or environmental interactions (Murphy and Sweedler, 2022). Techniques like gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS), have been extensively utilized for the analysis of numerous metabolites. GC-MS, known for its high sensitivity and capacity for high-throughput analysis, has emerged as a valuable tool in non-targeted metabolomics investigations (Papadimitropoulos et al., 2018).

Several studies have explored the metabolic alterations induced by CE using metabolomics. But thus far, these studies have primarily focused on analyzing one or two samples from serum, BAT, WAT, and liver (Lu et al., 2017; Hiroshima et al., 2018; Okamatsu-Ogura et al., 2020; Hou et al., 2021; Kovaničová et al., 2021; Chen et al., 2022). There has been a noticeable absence of simultaneous analysis of multiple tissues regarding the metabolic changes induced by CE. To the best of our knowledge, this present study represents the first metabolomics investigation simultaneously examining the effects of CE on multiple tissues, including serum, BAT, WAT, liver, spleen, and kidney. Our findings have the potential to provide fresh insights into the metabolic impacts induced by CE, and they may also present new ideas for treating metabolic diseases.

2 Materials and methods

2.1 Animal treatment

A total of 14 male C57BL/6J mice, aged 6 weeks, were obtained from Jinan Pengyue (Jinan, China). These mice were allowed unrestricted access to food and water for 1 week in a climate chamber set at an ambient temperature of 24°C ± 2°C, with a relative humidity of 40%, and a 12/12 h light/dark cycle. Subsequently, the mice were randomly assigned to two groups: the CE group (n = 7) and the control group (n = 7). The mice in the CE group were subjected to a temperature of 4°C for 4 h per day, continuously for two consecutive weeks. In contrast, the mice in the control group were maintained at room temperature (24°C ± 2°C) for the same duration. The body weight and food intake of the mice were recorded on a weekly basis. All experimental procedures were conducted in compliance with the Regulations of Experimental Animal Administration issued by the State Committee of Science and Technology of the People’s Republic of China and were approved by the University Ethics Committee (approval no. JNRM-2022-DW-054).

2.2 Sample collection and preparation

Twenty-four hours after the last CE session, the mice were humanely euthanized. Food was withheld before sample collection to ensure a 6-h fasting period for the mice. Anesthesia was administered through intraperitoneal injection of sodium pentobarbital at a dose of 50 mg/kg. Following enucleation of the eyeballs, blood samples were collected and subsequently centrifuged at 3,500 rpm for 8 min to obtain serum samples. For euthanasia, cervical dislocation was performed on all mice. Immediately afterward, each mouse was necropsied on an ice surface to acquire samples of interscapular BAT, inguinal WAT, liver, spleen and kidney. All tissue samples were washed with phosphate-buffered saline (PBS, pH 7.2), rapidly frozen in liquid nitrogen, and stored at −80°C in a refrigerator for future use.

To prepare the serum samples, 100 μL of serum was combined with 350 μL of methanol containing 100 μg/ml heptadecanoic acid. The solution was then vortexed and centrifuged at 14,000 rpm for 10 min at 4°C to collect the supernatant. The supernatant was transferred to a 2 mL tube and dried at 37°C under a flow of nitrogen gas. Subsequently, 80 μL of o-methylhydroxylamine hydrochloride (dissolved in pyridine at 15 mg/mL) was added to the dried sample and thoroughly mixed. The mixture was incubated at 70°C for 90 min. Next, 100 μL of N,O-bis (trimethylsilyl) trifluoroacetamide containing 1% trimethylchlorosilane (Sigma-Aldrich) was added to each sample, followed by a 60-min incubation at 70°C. The solution was then vortexed, centrifuged at 14,000 rpm for 2 min at 4°C, and filtered through a 0.22-μm filter membrane before GC-MS analysis. To prepare the tissue samples (BAT, WAT, liver, spleen, and kidney), 50 mg of each sample was homogenized in 1 mL of methanol containing 1 mg/mL heptadecanoic acid. After homogenization, the samples were centrifuged for 10 min at 20,913 × g at 4°C. The subsequent steps of the protocol were similar to those used for the serum samples. For the quality control samples (QCs), equal amounts of tissue samples from the control group and CE group were mixed together.

2.3 GC-MS analysis

The analysis of all samples was performed using a 7000°C mass spectrometer coupled with a 7890 B gas chromatograph system from Agilent Technologies (CA, United States). The separation of serum, BAT, WAT, liver, spleen and kidney samples was performed using an HP-5MS fused silica capillary column. Using helium gas as a carrier, a 1 µL aliquot of the derivative solution was processed in split mode (50:1), with the front inlet purge flow rate set to 3 mL/min and the gas flow rate set to 1 mL/min. The temperatures for the administration, transfer line, and ion source were maintained at 280°C, 250°C, and 230°C, respectively. The GC temperature program started at 60°C for 4 min, followed by an increase to 300°C at a rate of 8°C/min, and then held at 300°C for 5 min. For ionization, the voltage of the electron impact was set to −70 eV, and data acquisition occ1urred at a rate of 20 spectra per second. MS identification was performed using electrospray ionization (ESI) in full scan mode, with a mass-to-charge ratio (m/z) range of 50–800.

We have effectively uploaded the source data from our GC-MS analysis to MetaboLights. The identifier MTBLS8334 has been exclusively assigned to our research project. To access our study, simply follow this link: https://www.ebi.ac.uk/metabolights/MTBLS8334.

2.4 Multivariate statistical analysis

The initial analysis of the GC-MS data was conducted using Agilent Unknowns Analysis software and Mass Hunter Quantitative Analysis software from Agilent Technologies (United States). The SIMCA 14.1 software (Umetrics, Sweden) was employed for the statistical analysis of normalized peak area percentages. Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) was performed to differentiate between the CE group and control group. OPLS-DA is a statistical method utilized for multivariate data analysis, especially when dealing with high-dimensional data and a limited number of samples. The primary objective of OPLS-DA is to identify significant differences and relationships between two or more groups or classes within a dataset (Bylesjö et al., 2006). The method achieves this by modeling the systematic variation between classes, such as different treatment groups or disease states, while effectively removing unrelated or confounding variations. Variable Importance in Projection (VIP) values estimate the importance of each variable in the projection used in a least squares regression model and is often used for variable selection (Mehmood et al., 2012; Pinto et al., 2012). A variable with a VIP value close to or greater than one can be considered important in given model. Variables with VIP values significantly less than one are less important and might be good candidates for exclusion from the model. Model validation was further verified using the permutation test (200 permutations). Two-tailed Student's t-tests were conducted using SPSS 19.0 software (SPSS, Chicago, IL, United States). Then p values obtained from the t-test were corrected using fdrtool (Strimmer, 2008). Metabolites with VIP values >1.0 in the OPLS-DA analysis and p-values adjusted <0.05 were considered statistically significant.

Differential endogenous metabolites were imported into MetaboAnalyst 5.0 (http://www.metaboanalyst.ca) as well as the Kyoto Encyclopedia of Genes and Genomes (KEGG; http://www.kegg.jp) for metabolic pathway analysis. Metabolic pathways with impact values >0 and p-values adjusted <0.05 were considered as significantly affected pathways.

3 Results

3.1 Impact of CE on food intake and body weight

The CE group exhibited a significant increase in food intake compared to the control group (p-values <0.0001, Figure 1A). Conversely, the CE group showed a decrease in body weight (p-values <0.05, Figure 1B). These findings, consistent with the previous report (Yang et al., 2017), suggest that CE induces increased energy expenditure.

FIGURE 1

3.2 GC-MS total ion chromatograms of samples

The representative total ion chromatograms of QCs were shown in Figure 2; Details can be seen in Supplementary Material. The chromatograms demonstrate strong signal responses in all samples, indicating the detection of a diverse range of metabolites with a high peak capacity throughout the analysis. The retention time was consistent and the chromatograms of each tissue exhibited excellent reproducibility.

FIGURE 2

3.3 Multivariate statistical analysis

In OPLS-DA analysis, clear differences were observed between the CE group and the control group. The parameters (BAT R2X = 0.622, R2Y = 0.997, Q2 = 0.904; WAT R2X = 0.522, R2Y = 0.989, Q2 = 0.791; serum R2X = 0.59, R2Y = 0.918, Q2 = 0.731; liver R2X = 0.736, R2Y = 1, Q2 = 0.824; spleen R2X = 0.515, R2Y = 0.992, Q2 = 0.673; and kidney R2X = 0.669, R2Y = 0.997, Q2 = 0.805) indicated the effectiveness of the model, allowing clear differentiation between the CE group and the control group. Each parameter value was close to 1.0, signifying a stable and reliably predictive model. Model validation was verified through permutation tests (200 permutations), and the intersection of the blue regression lines (the Q2 points) and the vertical axis (on the left) all fell below zero (Figure 3).

FIGURE 3

We used MetaboAnalyst 5.0 to investigate metabolic differences between the two groups. Cluster analysis of the expression of metabolites in tissues revealed that most samples were grouped into two differentiated clusters with only a small part of the sample cluster overlapping (Figure 4). These results were consistent with those of the OPLS-DA analysis.

FIGURE 4

3.4 Identification of discriminant metabolites

The significance threshold for identifying differential compounds between the control and CE groups was set at VIP values >1.0 and p-values adjusted <0.05. Metabolites meeting these criteria were considered significantly different. Furthermore, metabolites exhibiting an upward trend were indicated by fold change (FC) > 1, while a downward trend was indicated by FC < 1. There were 32, 17, 21, 7, 16 and 26 discriminant metabolites in BAT, WAT, serum, liver, spleen and kidney between the two groups, respectively. A detailed summary of altered metabolites in major tissues following exposure to cold was provided in Table 1.

3.5 Analysis of metabolic pathways

In our study, we employed MetaboAnalyst 5.0 to investigate the metabolic pathways associated with the specific metabolites identified through a comparison between the CE and control groups. Through our analysis, we identified 12 metabolic pathways that exhibited significant differences (with p-values adjusted <0.05 and impact values >0). In BAT, these pathways included phenylalanine, tyrosine and tryptophan biosynthesis; glutathione metabolism; phenylalanine metabolism; glycine, serine and threonine metabolism; as well as arginine biosynthesis. In WAT, these pathways included arginine biosynthesis; glutathione metabolism; alanine, aspartate and glutamate metabolism; glyoxylate and dicarboxylate metabolism; primary bile acid biosynthesis; fatty acid biosynthesis; as well as D-Glutamine and D-glutamate metabolism. In the serum, these pathways included glutathione metabolism; glyoxylate and dicarboxylate metabolism. In the spleen, these pathways included phenylalanine, tyrosine and tryptophan biosynthesis; along with fatty acid biosynthesis. In the kidney, these pathways included aminoacyl-tRNA biosynthesis; nicotinate and nicotinamide metabolism; as well as alanine, aspartate and glutamate metabolism. The specifics of the metabolic pathway analysis were depicted in Table 2; Figure 5. A summary of metabolic pathways was presented in Figure 6.

FIGURE 5

FIGURE 6

4 Discussion

CE has significant effects on body metabolism. When exposed to cold temperatures, the body initiates a series of physiological responses to maintain core body temperature. These responses primarily involve increased energy expenditure to generate heat and altered metabolic processes.

In our study, we observed a significant increase in food intake in the CE group (p < 0.0001, Figure 1A); however, body weight was decreased (p < 0.05, Figure 1B). These findings were consistent with a previous report (Yang et al., 2017). One possible explanation for these results is that the cold-exposed mice needed to expend a significant amount of energy to generate heat and maintain their body temperature. The increase in energy consumption can potentially lead to higher food intake and a decrease in body weight.

Through GC-MS identification and statistical analysis, 32, 17, 21, 7, 16, and 26 differential metabolites were respectively identified in BAT, WAT, serum, liver, spleen and kidney after exposure to cold. These discriminant metabolites were involved in 12 pathways, mainly pertaining to amino acid, fatty acid, and energy metabolism.

Based on the altered metabolites and the associated differential metabolic pathways, it is evident that BAT was the most severely affected tissue. CE was also associated with varying degrees of metabolic effects in WAT, serum, liver, spleen and kidney. The metabolic effects in different tissues were discussed separately below.

4.1 BAT metabolism analysis

Branched chain amino acids (BCAAs), including leucine, isoleucine and valine, are essential amino acids. BCAAs can be oxidized to provide energy for cells (Sivanand and Vander Heiden, 2020). They are commonly converted into branched-chain alpha-ketoacids (BCKAs) through the catalytic action of branched-chain aminotransferases (BCATs) (Ichihara and Koyama, 1966). Subsequently, the branched-chain alpha-ketoacid dehydrogenase complex further converts the BCKAs into branched-chain acyl-CoA derivatives, which can then be metabolized into either acetyl-CoA or succinyl-CoA (Johnson and Connelly, 1972; Neinast et al., 2019). These metabolites eventually enter the tricarboxylic acid (TCA) cycle. In our study, we found that CE significantly increased the content of BCAAs in BAT. These findings were consistent with previous studies (Lu et al., 2017; Okamatsu-Ogura et al., 2020). According to Yoneshiro et al. (Yoneshiro et al., 2019), the higher BCAAs content in BAT is attributed to increased uptake from circulation. Indeed, previous studies have reported that the activity of BCATs increases following CE in rats (López-Soriano and Alemany, 1987). Consistent with this, more than 60% of genes encoding BCAAs catabolic enzymes, including the gene for the rate-limiting enzyme BCAT2, were more highly expressed in brown adipocytes relative to white adipocytes (Yoneshiro et al., 2019). Taken together, these data indicated that in vivo, BCAAs are used as significant energy sources for thermogenesis in BAT during CE. In our study, we observed a significant decrease in serum levels of BCAAs after exposure to cold, which was consistent with a previous human study (Kovaničová et al., 2021). Notably, elevated circulating BCAAs levels are closely linked to insulin resistance, obesity, and type 2 diabetes (Huffman et al., 2009; Newgard et al., 2009; Wang et al., 2011). CE enhances BCAAs clearance from the bloodstream by significantly increasing mitochondrial BCAAs uptake and oxidation in BAT (Yoneshiro et al., 2019). This process depends on the presence of SLC25A44, a mitochondrial BCAAs transporter in brown adipocytes, which might contribute to improved metabolic health (Yoneshiro et al., 2019). However, it remains unclear whether CE promotes the expression of SLC25A44, thereby enhancing BCAAs uptake. Further research is required.

Our study revealed a significant increase in glutamate levels in BAT following cold exposure, which aligns with previous research findings (Okamatsu-Ogura et al., 2020). A human study utilizing microdialysis demonstrated that CE increased glutamate uptake specifically by BAT (Weir et al., 2018), providing a potential explanation for the observed elevation in glutamate levels. Furthermore, exposure to cold in BAT resulted in an elevation of the expression of genes encoding glutamate dehydrogenase (GLUD1), an enzyme responsible for catalyzing the production of α-ketoglutaric acid from glutamate (Okamatsu-Ogura et al., 2020). These findings collectively suggest that glutamate serves as an important energy substrate for BAT thermogenesis in vivo. Additionally, we observed an increase in glycine content and an enhanced activity of the glutathione metabolism pathway in BAT after exposure to cold. CE triggers adaptive responses in the body, including the production of reactive oxygen species (ROS) as byproducts of increased metabolism and thermogenesis (Wang et al., 2015; Chouchani et al., 2016; Lettieri-Barbato, 2019). Uncontrolled ROS accumulation can lead to cellular oxidative damage. However, glutathione plays a crucial role as an antioxidant by scavenging ROS and protecting cells from oxidative stress. The observed enhancement in the glutathione metabolism pathway in BAT indicates an adaptive response to counteract the increased ROS generation during CE (Mory et al., 1983).

Lactic acid is produced through anaerobic glucose metabolism. Glucose transporters (GLUTs) play a vital role in cellular glucose uptake. In humans, there are a total of 14 different types of GLUTs (Thorens and Mueckler, 2010), with GLUT1 and GLUT4 expressed in BAT (James et al., 1988; James et al., 1989; Pessin and Bell, 1992). CE promotes a significant increase in the expression of GLUT1 and GLUT4 in mouse BAT (Daikoku et al., 2000; Yu et al., 2002). In line with this, studies using 18F-fluorodeoxyglucose-positron emission tomography/computed tomography (18F-FDG-PET/CT) have shown that CE enhances glucose uptake in the BAT of healthy humans (van Marken Lichtenbelt et al., 2009; Virtanen et al., 2009; Ouellet et al., 2012; Lee et al., 2016), In vivo [U-¹³C]glucose tracing experiments have also demonstrated that chronic cold exposure significantly increases glucose oxidation in BAT by promoting glucose flux into the mitochondrial TCA cycle (Wang et al., 2020). Notably, genes associated with glucose uptake and glycolysis were remarkably unregulated in cold-activated BAT (Hao et al., 2015). Taken together, these data indicate that CE enhances glucose uptake and oxidation. As we know, lactate dehydrogenases (LDH) play a crucial role in catalyzing the reversible conversion of pyruvate to lactate. It was observed that exposure to cold significantly induced the expression of Ldha and Ldhb mRNAs, as well as the LDHA protein in BAT (Hao et al., 2015). So, it is not surprising that lactic acid levels were significantly increased upon CE in BAT in our study. Previously considered a metabolic waste product, lactic acid is now recognized as a crucial metabolic fuel in BAT (Hui et al., 2017; Rabinowitz and Enerbäck, 2020). Recent studies have demonstrated that carbon derived from lactic acid can be utilized by brown adipocytes to synthesize fatty acids (Saggerson et al., 1988). As a result, lactic acid may play a role in swiftly replenishing BAT triglycerides following CE (Carpentier et al., 2023).

4.2 Serum metabolism analysis

Our study revealed that CE led to a decrease in the levels of fatty acids, specifically palmitic acid and stearic acid, as well as cholesterol in the serum of mice. CE induces the sympathetic nerves to release norepinephrine, which subsequently binds to β3-adrenoceptors (β3-ARs) found on the surface of adipocytes. This binding activates the cAMP–protein kinase A (PKA) signaling pathway, ultimately promoting intracellular TAG lipolysis (Cannon and Nedergaard, 2004; Zechner et al., 2012; Caron et al., 2018). When lipolysis occurs in WAT, it releases a significant amount of fatty acids into the bloodstream, which serve as the main energy substrate for BAT thermogenesis (Wang et al., 2021; Carpentier et al., 2023). Activated brown/beige adipocytes exhibit enhanced fatty acid uptake, as evidenced by 18F-fluoro-thiaheptadecanoic acid-positron emission tomography/computed tomography (18F-FTHA-PET/CT) imaging (Ouellet et al., 2012; Blondin et al., 2015). Moreover, CE causes a selective increase in lipoprotein lipase (LPL) expression in BAT (Chondronikola et al., 2016), while reducing the expression of angiopoietin-like 4 (ANGPTL4), an inhibitor of LPL activity (Dijk et al., 2015). As a result, activated BAT can effectively eliminate the majority of circulating triglyceride-rich lipoprotein (TRL) lipids (Bartelt et al., 2011; Berbée et al., 2015; Khedoe et al., 2015). Subsequently, fatty acids are transported into the mitochondria of brown fat cells via carnitine acyltransferases (CPTs), where they undergo β-oxidation and participate in the TCA cycle. Earlier research has demonstrated that CE can upregulate mRNA expression of cpt (Yu et al., 2002), thereby stimulating the β-oxidation of fatty acids in BAT of mice. Cold-activated BAT in mice holds promise for alleviating hyperlipidemia, providing benefits in conditions such as obesity or genetic hypertriglyceridemia (Bartelt et al., 2011; Berbée et al., 2015). Moreover, this activation indirectly reduces hypercholesterolemia by facilitating an increased hepatic uptake of cholesterol-enriched lipoprotein remnants, thereby offering protection against the development of atherosclerosis (Berbée et al., 2015; Bartelt et al., 2017).

4.3 Liver metabolism analysis

The liver’s role in maintaining cholesterol balance is pivotal due to its ability to uptake, synthesize, convert cholesterol into bile acids, and excrete cholesterol within very low-density lipoprotein (VLDL) particles. Our study observed a substantial reduction in liver cholesterol levels following exposure to cold conditions. This phenomenon could be attributed to the cold-induced conversion of cholesterol into bile acids in the liver (Worthmann et al., 2017). Specifically, the researchers found that the levels of most bile acid species in the liver were notably higher in mice housed in cold conditions compared to the control group in warm conditions. Additionally, the study revealed a significant upregulation in the expression of the important gene Cyp7b, which is involved in alternative bile acid synthesis pathways, in response to CE (Worthmann et al., 2017).

4.4 Spleen metabolism analysis

The spleen plays a crucial role as a lymphatic and immune organ in the human body. Our research has revealed a decrease in the levels of lactic acid and malic acid within the spleen during CE, indicating a weakening of anaerobic respiration and TCA cycle processes. Furthermore, a transcriptomic study has indicated that gene sets associated with biological processes, including the innate immune response in mucosa, neutrophil-mediated killing of symbiont cells, inflammatory response, and various other responsive systems, were downregulated in the spleen (Hadadi et al., 2022). This phenomenon can be attributed to the spleen’s primary immunological functions, a biological process that might become less effective during cold conditions. This aligns with the concept that maintaining immune functions requires significant energy expenditure, potentially competing with other energy-demanding functions, including thermogenesis (Spiljar et al., 2021). Furthermore, the process of triglyceride catabolism was found to be subdued in the spleen during cold exposure (Hadadi et al., 2022). It is not unexpected that the levels of glycerol and fatty acids (such as palmitic acid, stearic acid, myristic acid, glycerol monostearate, and 1-Monopalmitin) showed a decrease within the spleen according to our study. To some extent, this could suggest a redistribution of metabolic energy towards the tissues essential for responding to cold, as a consequence of prioritizing other maintenance-related biological processes.

4.5 Kidney metabolism analysis

The kidneys play a vital role in regulating various essential physiological functions. They are responsible for producing urine, removing metabolic toxins, and maintaining the balance of water, electrolytes, and acid-base levels (Imenez Silva and Mohebbi, 2022). Consequently, the kidneys are organs with high metabolic activity (Clark and Parikh, 2020). Research has indicated that the kidneys have the second highest concentration of mitochondria per weight, surpassed only by the heart (Pagliarini et al., 2008). To meet their energy demands, the kidneys metabolize significant amounts of nutrients, such as fatty acids, glucose, and amino acids.

Nicotinamide (NAM), a member of the vitamin B3 group, has the capacity to generate nicotinamide adenine dinucleotide (NAD⁺). NAD⁺ plays a crucial role as an electron acceptor in several catabolic processes, including glycolysis, the TCA cycle, and fatty acid β-oxidation (FAO). Subsequently, the acquired electrons are transferred to oxygen through the mitochondrial respiratory chain, leading to the conversion of ADP into ATP via phosphorylation (Fontecha-Barriuso et al., 2021). In our research, following CE, we observed a significant increase in NAM levels and a notable enhancement in the pathway of nicotinate and NAM metabolism. Based on this, we hypothesized that CE may enhance renal energy metabolism. However, a recent time-series study conducted in South Korea identified a correlation between CE and hospital admissions and deaths associated with acute kidney injury (AKI) (Kim et al., 2023). These conflicting results could be attributed to various factors, such as the duration and intensity of CE, the individual’s health condition, and their capacity to adapt to the cold. Consequently, further research is necessary to investigate the impact of CE on renal metabolism in both animal models and clinical settings.

4.6 Interorgan communication in response to CE

The process of cold-induced thermogenesis is intricately tied to a well-coordinated metabolic adjustment program spanning multiple tissues. This program functions to uphold the equilibrium of fuel and energy within the body. Cold stimuli trigger the initiation of lipolysis within WAT, leading to the liberation of free fatty acids (FFAs). These FFAs either act as direct substrates for BAT-mediated thermogenesis or journey to the liver. In the liver, these FFAs undergo a transformation into acylcarnitines, facilitated by an increase in the expression of the CPT1 gene, brought about by the activation of HNF4a (Simcox et al., 2017). The acylcarnitines produced in the liver are predominantly taken up by BAT (Simcox et al., 2017), where they contribute significantly to thermogenic processes. Traditionally, bile acids, generated solely by the liver through the conversion of cholesterol, can follow either the classical pathway or the alternative pathway. Cold exposure particularly triggers the activation of the alternative pathway, leading to a substantial rise in bile acids expelled through feces (Worthmann et al., 2017). This surge is accompanied by a distinct transformation in the composition of the gut microbiome (Worthmann et al., 2017). This altered configuration of the microbiome has the potential to result in the creation of microbial metabolites possessing inherent thermogenic properties (Worthmann et al., 2017).

While our study provided valuable insights, it is crucial to acknowledge its limitations. In order to achieve a comprehensive understanding of the metabolic alterations induced by CE, it is imperative to explore its effects on other tissues, such as the heart, lung, cortex, hippocampus, stomach, pancreas, skin, and bone. However, obtaining samples from multiple tissues is challenging and time-consuming, as it involves the collection, rinsing, and freezing of multiple tissues in liquid nitrogen. Additionally, many metabolites have a short turnover time, which could potentially impact the accuracy of metabolite measurements. Furthermore, it’s worth noting that the use of anesthesia during the sampling procedure can induce metabolic changes in tissues, potentially influencing the results. Another factor to consider is that food intake was not normalized between groups, which could contribute to the observed differences in metabolites (Jang et al., 2018).

5 Conclusion

Metabolic changes induced by CE in multiple tissues in mice were analyzed using GC-MS, providing a deeply understanding of the effects of CE. We observed that CE altered the levels of various metabolites across multiple tissues, with notable impacts on amino acids metabolism, fatty acids metabolism and energy metabolism. Our findings provide systematic insights into the metabolic impacts induced by CE and may also present new ideas for treating metabolic diseases.

References

• 1

  BarteltA.BrunsO. T.ReimerR.HohenbergH.IttrichH.PeldschusK.et al (2011). Brown adipose tissue activity controls triglyceride clearance. Nat. Med.17 (2), 200–205. 10.1038/nm.2297

  • CrossRef
  • Google Scholar

• 2

  BarteltA.JohnC.SchaltenbergN.BerbéeJ. F. P.WorthmannA.CherradiM. L.et al (2017). Thermogenic adipocytes promote HDL turnover and reverse cholesterol transport. Nat. Commun.8, 15010. 10.1038/ncomms15010

  • CrossRef
  • Google Scholar

• 3

  BecherT.PalanisamyS.KramerD. J.EljalbyM.MarxS. J.WibmerA. G.et al (2021). Brown adipose tissue is associated with cardiometabolic health. Nat. Med.27 (1), 58–65. 10.1038/s41591-020-1126-7

  • CrossRef
  • Google Scholar

• 4

  BerbéeJ. F.BoonM. R.KhedoeP. P.BarteltA.SchleinC.WorthmannA.et al (2015). Brown fat activation reduces hypercholesterolaemia and protects from atherosclerosis development. Nat. Commun.6, 6356. 10.1038/ncomms7356

  • CrossRef
  • Google Scholar

• 5

  BlondinD. P.LabbéS. M.NollC.KunachM.PhoenixS.GuérinB.et al (2015). Selective impairment of glucose but not fatty acid or oxidative metabolism in Brown adipose tissue of subjects with type 2 diabetes. Diabetes64 (7), 2388–2397. 10.2337/db14-1651

  • CrossRef
  • Google Scholar

• 6

  BurkartK. G.BrauerM.AravkinA. Y.GodwinW. W.HayS. I.HeJ.et al (2021). Estimating the cause-specific relative risks of non-optimal temperature on daily mortality: A two-part modelling approach applied to the global burden of disease study. Lancet398 (10301), 685–697. 10.1016/s0140-6736(21)01700-1

  • CrossRef
  • Google Scholar

• 7

  BylesjöM.RantalainenM.CloarecO.NicholsonJ. K.HolmesE.TryggJ. (2006). OPLS discriminant analysis: combining the strengths of PLS‐DA and SIMCA classification. J. Chemom. A J. Chemom. Soc.20 (8‐10), 341–351. 10.1002/cem.1006

  • CrossRef
  • Google Scholar

• 8

  CannonB.NedergaardJ. (2004). Brown adipose tissue: function and physiological significance. Physiol. Rev.84 (1), 277–359. 10.1152/physrev.00015.2003

  • CrossRef
  • Google Scholar

• 9

  CaronA.LeeS.ElmquistJ. K.GautronL. (2018). Leptin and brain-adipose crosstalks. Nat. Rev. Neurosci.19 (3), 153–165. 10.1038/nrn.2018.7

  • CrossRef
  • Google Scholar

• 10

  CarpentierA. C.BlondinD. P.HamanF.RichardD. (2023). Brown adipose tissue-A translational perspective. Endocr. Rev.44 (2), 143–192. 10.1210/endrev/bnac015

  • CrossRef
  • Google Scholar

• 11

  ChartoumpekisD. V.HabeosI. G.ZirosP. G.PsyrogiannisA. I.KyriazopoulouV. E.PapavassiliouA. G. (2011). Brown adipose tissue responds to cold and adrenergic stimulation by induction of FGF21. Mol. Med.17 (7-8), 736–740. 10.2119/molmed.2011.00075

  • CrossRef
  • Google Scholar

• 12

  ChenY.JiH.GuoJ.ChenY.LiW.WangS.et al (2022). Non-targeted metabolomics analysis based on LC-MS to assess the effects of different cold exposure times on piglets. Front. Physiol.13, 853995. 10.3389/fphys.2022.853995

  • CrossRef
  • Google Scholar

• 13

  ChondronikolaM.VolpiE.BørsheimE.PorterC.SarafM. K.AnnamalaiP.et al (2016). Brown adipose tissue activation is linked to distinct systemic effects on lipid metabolism in humans. Cell Metab.23 (6), 1200–1206. 10.1016/j.cmet.2016.04.029

  • CrossRef
  • Google Scholar

• 14

  ChouchaniE. T.KazakL.JedrychowskiM. P.LuG. Z.EricksonB. K.SzpytJ.et al (2016). Mitochondrial ROS regulate thermogenic energy expenditure and sulfenylation of UCP1. Nature532 (7597), 112–116. 10.1038/nature17399

  • CrossRef
  • Google Scholar

• 15

  ClarkA. J.ParikhS. M. (2020). Mitochondrial metabolism in acute kidney injury. Semin. Nephrol.40 (2), 101–113. 10.1016/j.semnephrol.2020.01.002

  • CrossRef
  • Google Scholar

• 16

  CypessA. M.LehmanS.WilliamsG.TalI.RodmanD.GoldfineA. B.et al (2009). Identification and importance of brown adipose tissue in adult humans. N. Engl. J. Med.360 (15), 1509–1517. 10.1056/NEJMoa0810780

  • CrossRef
  • Google Scholar

• 17

  DaikokuT.ShinoharaY.ShimaA.YamazakiN.TeradaH. (2000). Specific elevation of transcript levels of particular protein subtypes induced in brown adipose tissue by cold exposure. Biochim. Biophys. Acta1457 (3), 263–272. 10.1016/s0005-2728(00)00107-9

  • CrossRef
  • Google Scholar

• 18

  DijkW.HeineM.VergnesL.BoonM. R.SchaartG.HesselinkM. K.et al (2015). ANGPTL4 mediates shuttling of lipid fuel to brown adipose tissue during sustained cold exposure. Elife4, e08428. 10.7554/eLife.08428

  • CrossRef
  • Google Scholar

• 19

  FatimaS. H.RothmoreP.GilesL. C.VargheseB. M.BiP. (2021). Extreme heat and occupational injuries in different climate zones: A systematic review and meta-analysis of epidemiological evidence. Environ. Int.148, 106384. 10.1016/j.envint.2021.106384

  • CrossRef
  • Google Scholar

• 20

  FisherF. M.KleinerS.DourisN.FoxE. C.MepaniR. J.VerdeguerF.et al (2012). FGF21 regulates PGC-1α and browning of white adipose tissues in adaptive thermogenesis. Genes Dev.26 (3), 271–281. 10.1101/gad.177857.111

  • CrossRef
  • Google Scholar

• 21

  FisherF. M.Maratos-FlierE. (2016). Understanding the physiology of FGF21. Annu. Rev. Physiol.78, 223–241. 10.1146/annurev-physiol-021115-105339

  • CrossRef
  • Google Scholar

• 22

  Fontecha-BarriusoM.Lopez-DiazA. M.CarriazoS.OrtizA.SanzA. B. (2021). Nicotinamide and acute kidney injury. Clin. Kidney J.14 (12), 2453–2462. 10.1093/ckj/sfab173

  • CrossRef
  • Google Scholar

• 23

  HadadiN.SpiljarM.SteinbachK.ÇolakoğluM.ChevalierC.SalinasG.et al (2022). Comparative multi-tissue profiling reveals extensive tissue-specificity in transcriptome reprogramming during thermal adaptation. Elife11, e78556. 10.7554/eLife.78556

  • CrossRef
  • Google Scholar

• 24

  HanssenM. J.HoeksJ.BransB.van der LansA. A.SchaartG.van den DriesscheJ. J.et al (2015). Short-term cold acclimation improves insulin sensitivity in patients with type 2 diabetes mellitus. Nat. Med.21 (8), 863–865. 10.1038/nm.3891

  • CrossRef
  • Google Scholar

• 25

  HaoQ.YadavR.BasseA. L.PetersenS.SonneS. B.RasmussenS.et al (2015). Transcriptome profiling of brown adipose tissue during cold exposure reveals extensive regulation of glucose metabolism. Am. J. Physiol. Endocrinol. Metab.308 (5), E380–E392. 10.1152/ajpendo.00277.2014

  • CrossRef
  • Google Scholar

• 26

  HiroshimaY.YamamotoT.WatanabeM.BabaY.ShinoharaY. (2018). Effects of cold exposure on metabolites in brown adipose tissue of rats. Mol. Genet. Metab. Rep.15, 36–42. 10.1016/j.ymgmr.2018.01.005

  • CrossRef
  • Google Scholar

• 27

  HondaresE.IglesiasR.GiraltA.GonzalezF. J.GiraltM.MampelT.et al (2011). Thermogenic activation induces FGF21 expression and release in brown adipose tissue. J. Biol. Chem.286 (15), 12983–12990. 10.1074/jbc.M110.215889

  • CrossRef
  • Google Scholar

• 28

  HouD. M.JiaT.LiQ.WangZ. K.ZhuW. L. (2021). Metabonomics of white adipose tissue and brown adipose tissue in Tupaia belangeri during cold acclimation. Comp. Biochem. Physiol. Part D. Genomics Proteomics38, 100823. 10.1016/j.cbd.2021.100823

  • CrossRef
  • Google Scholar

• 29

  HuffmanK. M.ShahS. H.StevensR. D.BainJ. R.MuehlbauerM.SlentzC. A.et al (2009). Relationships between circulating metabolic intermediates and insulin action in overweight to obese, inactive men and women. Diabetes Care32 (9), 1678–1683. 10.2337/dc08-2075

  • CrossRef
  • Google Scholar

• 30

  HuiS.GhergurovichJ. M.MorscherR. J.JangC.TengX.LuW.et al (2017). Glucose feeds the TCA cycle via circulating lactate. Nature551 (7678), 115–118. 10.1038/nature24057

  • CrossRef
  • Google Scholar

• 31

  IchiharaA.KoyamaE. (1966). Transaminase of branched chain amino acids. I. Branched chain amino acids-alpha-ketoglutarate transaminase. J. Biochem.59 (2), 160–169. 10.1093/oxfordjournals.jbchem.a128277

  • CrossRef
  • Google Scholar

• 32

  Imenez SilvaP. H.MohebbiN. (2022). Kidney metabolism and acid-base control: back to the basics. Pflugers Arch.474 (8), 919–934. 10.1007/s00424-022-02696-6

  • CrossRef
  • Google Scholar

• 33

  IwenK. A.BackhausJ.CassensM.WaltlM.HedesanO. C.MerkelM.et al (2017). Cold-induced Brown adipose tissue activity alters plasma fatty acids and improves glucose metabolism in men. J. Clin. Endocrinol. Metab.102 (11), 4226–4234. 10.1210/jc.2017-01250

  • CrossRef
  • Google Scholar

• 34

  JamesD. E.BrownR.NavarroJ.PilchP. F. (1988). Insulin-regulatable tissues express a unique insulin-sensitive glucose transport protein. Nature333 (6169), 183–185. 10.1038/333183a0

  • CrossRef
  • Google Scholar

• 35

  JamesD. E.StrubeM.MuecklerM. (1989). Molecular cloning and characterization of an insulin-regulatable glucose transporter. Nature338 (6210), 83–87. 10.1038/338083a0

  • CrossRef
  • Google Scholar

• 36

  JangC.ChenL.RabinowitzJ. D. (2018). Metabolomics and isotope tracing. Cell173 (4), 822–837. 10.1016/j.cell.2018.03.055

  • CrossRef
  • Google Scholar

• 37

  JiangY.RoseA. J.SijmonsmaT. P.BröerA.PfenningerA.HerzigS.et al (2015). Mice lacking neutral amino acid transporter B(0)AT1 (Slc6a19) have elevated levels of FGF21 and GLP-1 and improved glycaemic control. Mol. Metab.4 (5), 406–417. 10.1016/j.molmet.2015.02.003

  • CrossRef
  • Google Scholar

• 38

  JohnsonW. A.ConnellyJ. L. (1972). Cellular localization and characterization of bovine liver branched-chain -keto acid dehydrogenases. Biochemistry11 (10), 1967–1973. 10.1021/bi00760a036

  • CrossRef
  • Google Scholar

• 39

  KajimuraS.SaitoM. (2014). A new era in brown adipose tissue biology: molecular control of brown fat development and energy homeostasis. Annu. Rev. Physiol.76, 225–249. 10.1146/annurev-physiol-021113-170252

  • CrossRef
  • Google Scholar

• 40

  KarushevaY.KoesslerT.StrassburgerK.MarkgrafD.MastrototaroL.JelenikT.et al (2019). Short-term dietary reduction of branched-chain amino acids reduces meal-induced insulin secretion and modifies microbiome composition in type 2 diabetes: a randomized controlled crossover trial. Am. J. Clin. Nutr.110 (5), 1098–1107. 10.1093/ajcn/nqz191

  • CrossRef
  • Google Scholar

• 41

  KhedoeP. P.HoekeG.KooijmanS.DijkW.BuijsJ. T.KerstenS.et al (2015). Brown adipose tissue takes up plasma triglycerides mostly after lipolysis. J. Lipid Res.56 (1), 51–59. 10.1194/jlr.M052746

  • CrossRef
  • Google Scholar

• 42

  KimK. N.ShinM. K.LimY. H.BaeS.KimJ. H.HwangS. S.et al (2023). Associations of cold exposure with hospital admission and mortality due to acute kidney injury: A nationwide time-series study in Korea. Sci. Total Environ.863, 160960. 10.1016/j.scitotenv.2022.160960

  • CrossRef
  • Google Scholar

• 43

  KovaničováZ.KarhánekM.KurdiováT.BalážM.WolfrumC.UkropcováB.et al (2021). Metabolomic analysis reveals changes in plasma metabolites in response to acute cold stress and their relationships to metabolic health in cold-acclimatized humans. Metabolites11 (9), 619. 10.3390/metabo11090619

  • CrossRef
  • Google Scholar

• 44

  LeanM. E. (1989). Brown adipose tissue in humans. Proc. Nutr. Soc.48 (2), 243–256. 10.1079/pns19890036

  • CrossRef
  • Google Scholar

• 45

  LeeP.BovaR.SchofieldL.BryantW.DieckmannW.SlatteryA.et al (2016). Brown adipose tissue exhibits a glucose-responsive thermogenic biorhythm in humans. Cell Metab.23 (4), 602–609. 10.1016/j.cmet.2016.02.007

  • CrossRef
  • Google Scholar

• 46

  Lettieri-BarbatoD. (2019). Redox control of non-shivering thermogenesis. Mol. Metab.25, 11–19. 10.1016/j.molmet.2019.04.002

  • CrossRef
  • Google Scholar

• 47

  López-SorianoF. J.AlemanyM. (1987). Effect of cold-temperature exposure and acclimation on amino acid pool changes and enzyme activities of rat brown adipose tissue. Biochim. Biophys. Acta925 (3), 265–271. 10.1016/0304-4165(87)90191-7

  • CrossRef
  • Google Scholar

• 48

  LuX.SolmonsonA.LodiA.NowinskiS. M.SentandreuE.RileyC. L.et al (2017). The early metabolomic response of adipose tissue during acute cold exposure in mice. Sci. Rep.7 (1), 3455. 10.1038/s41598-017-03108-x

  • CrossRef
  • Google Scholar

• 49

  MehmoodT.LilandK. H.SnipenL.SæbøS. (2012). A review of variable selection methods in least squares regression. Chemom. intelligent laboratory Syst.118, 62–69. 10.1016/j.chemolab.2012.07.010

  • CrossRef
  • Google Scholar

• 50

  MoryG.BalD.RicquierD. (1983). Parallel increase of ascorbic acid and glutathione contents in brown adipose tissue during chronic cold exposure. Experientia39 (7), 734–735. 10.1007/bf01990299

  • CrossRef
  • Google Scholar

• 51

  MurphyS. E.SweedlerJ. V. (2022). Metabolomics-based mass spectrometry methods to analyze the chemical content of 3D organoid models. Analyst147 (13), 2918–2929. 10.1039/d2an00599a

  • CrossRef
  • Google Scholar

• 52

  NeinastM.MurashigeD.AranyZ. (2019). Branched chain amino acids. Annu. Rev. Physiol.81, 139–164. 10.1146/annurev-physiol-020518-114455

  • CrossRef
  • Google Scholar

• 53

  NewgardC. B.AnJ.BainJ. R.MuehlbauerM. J.StevensR. D.LienL. F.et al (2009). A branched-chain amino acid-related metabolic signature that differentiates obese and lean humans and contributes to insulin resistance. Cell Metab.9 (4), 311–326. 10.1016/j.cmet.2009.02.002

  • CrossRef
  • Google Scholar

• 54

  NguyenK. D.QiuY.CuiX.GohY. P.MwangiJ.DavidT.et al (2011). Alternatively activated macrophages produce catecholamines to sustain adaptive thermogenesis. Nature480 (7375), 104–108. 10.1038/nature10653

  • CrossRef
  • Google Scholar

• 55

  NichollsD. G.BernsonV. S.HeatonG. M. (1978). The identification of the component in the inner membrane of brown adipose tissue mitochondria responsible for regulating energy dissipation. Exp. Suppl.32, 89–93. 10.1007/978-3-0348-5559-4_9

  • CrossRef
  • Google Scholar

• 56

  Okamatsu-OguraY.KurodaM.TsutsumiR.TsubotaA.SaitoM.KimuraK.et al (2020). UCP1-dependent and UCP1-independent metabolic changes induced by acute cold exposure in brown adipose tissue of mice. Metabolism113, 154396. 10.1016/j.metabol.2020.154396

  • CrossRef
  • Google Scholar

• 57

  OravaJ.NuutilaP.LidellM. E.OikonenV.NoponenT.ViljanenT.et al (2011). Different metabolic responses of human brown adipose tissue to activation by cold and insulin. Cell Metab.14 (2), 272–279. 10.1016/j.cmet.2011.06.012

  • CrossRef
  • Google Scholar

• 58

  OuelletV.LabbéS. M.BlondinD. P.PhoenixS.GuérinB.HamanF.et al (2012). Brown adipose tissue oxidative metabolism contributes to energy expenditure during acute cold exposure in humans. J. Clin. Invest.122 (2), 545–552. 10.1172/jci60433

  • CrossRef
  • Google Scholar

• 59

  PagliariniD. J.CalvoS. E.ChangB.ShethS. A.VafaiS. B.OngS. E.et al (2008). A mitochondrial protein compendium elucidates complex I disease biology. Cell134 (1), 112–123. 10.1016/j.cell.2008.06.016

  • CrossRef
  • Google Scholar

• 60

  PapadimitropoulosM. P.VasilopoulouC. G.Maga-NteveC.KlapaM. I. (2018). Untargeted GC-MS metabolomics. Methods Mol. Biol.1738, 133–147. 10.1007/978-1-4939-7643-0_9

  • CrossRef
  • Google Scholar

• 61

  PessinJ. E.BellG. I. (1992). Mammalian facilitative glucose transporter family: structure and molecular regulation. Annu. Rev. Physiol.54, 911–930. 10.1146/annurev.ph.54.030192.004403

  • CrossRef
  • Google Scholar

• 62

  PintoR. C.TryggJ.GottfriesJ. (2012). Advantages of orthogonal inspection in chemometrics. J. Chemom.26 (6), 231–235. 10.1002/cem.2441

  • CrossRef
  • Google Scholar

• 63

  RabinowitzJ. D.EnerbäckS. (2020). Lactate: the ugly duckling of energy metabolism. Nat. Metab.2 (7), 566–571. 10.1038/s42255-020-0243-4

  • CrossRef
  • Google Scholar

• 64

  RicquierD. (2017). UCP1, the mitochondrial uncoupling protein of brown adipocyte: a personal contribution and a historical perspective. Biochimie134, 3–8. 10.1016/j.biochi.2016.10.018

  • CrossRef
  • Google Scholar

• 65

  SaggersonE. D.McAllisterT. W.BahtH. S. (1988). Lipogenesis in rat brown adipocytes. Effects of insulin and noradrenaline, contributions from glucose and lactate as precursors and comparisons with white adipocytes. Biochem. J.251 (3), 701–709. 10.1042/bj2510701

  • CrossRef
  • Google Scholar

• 66

  ShahH.KramerA.MullinsC. A.MatternM.GannabanR. B.TownsendR. L.et al (2023). Reduction of plasma BCAAs following roux-en-Y gastric bypass surgery is primarily mediated by FGF21. Nutrients15 (7), 1713. 10.3390/nu15071713

  • CrossRef
  • Google Scholar

• 67

  SimcoxJ.GeogheganG.MaschekJ. A.BensardC. L.PasqualiM.MiaoR.et al (2017). Global analysis of plasma lipids identifies liver-derived acylcarnitines as a fuel source for Brown fat thermogenesis. Cell Metab.26 (3), 509–522. e6. 10.1016/j.cmet.2017.08.006

  • CrossRef
  • Google Scholar

• 68

  SinghR.ParveenM.BasgenJ. M.FazelS.MesheshaM. F.ThamesE. C.et al (2016). Increased expression of beige/Brown adipose markers from host and breast cancer cells influence xenograft formation in mice. Mol. Cancer Res.14 (1), 78–92. 10.1158/1541-7786.MCR-15-0151

  • CrossRef
  • Google Scholar

• 69

  SivanandS.Vander HeidenM. G. (2020). Emerging roles for branched-chain amino acid metabolism in cancer. Cancer Cell37 (2), 147–156. 10.1016/j.ccell.2019.12.011

  • CrossRef
  • Google Scholar

• 70

  SpiljarM.SteinbachK.RigoD.Suárez-ZamoranoN.WagnerI.HadadiN.et al (2021). Cold exposure protects from neuroinflammation through immunologic reprogramming. Cell Metab.33 (11), 2231–2246.e8. 10.1016/j.cmet.2021.10.002

  • CrossRef
  • Google Scholar

• 71

  StanfordK. I.MiddelbeekR. J.TownsendK. L.AnD.NygaardE. B.HitchcoxK. M.et al (2013). Brown adipose tissue regulates glucose homeostasis and insulin sensitivity. J. Clin. Invest.123 (1), 215–223. 10.1172/jci62308

  • CrossRef
  • Google Scholar

• 72

  StrimmerK. (2008). fdrtool: a versatile R package for estimating local and tail area-based false discovery rates. Bioinformatics24 (12), 1461–1462. 10.1093/bioinformatics/btn209

  • CrossRef
  • Google Scholar

• 73

  TalanM. I.EngelB. T.WhitakerJ. R. (1985). A longitudinal study of tolerance to cold stress among C57BL/6J mice. J. Gerontol.40 (1), 8–14. 10.1093/geronj/40.1.8

  • CrossRef
  • Google Scholar

• 74

  TatelmanH. M.TalanM. I. (1990). Metabolic heat production during repeated testing at 24 degrees C and 6 degrees C in adult and aged male C57bl/6J mice: the effect of physical restraint before cold stress. Exp. Gerontol.25 (5), 459–467. 10.1016/0531-5565(90)90034-y

  • CrossRef
  • Google Scholar

• 75

  ThorensB.MuecklerM. (2010). Glucose transporters in the 21st century. Am. J. Physiol. Endocrinol. Metab.298 (2), E141–E145. 10.1152/ajpendo.00712.2009

  • CrossRef
  • Google Scholar

• 76

  van Marken LichtenbeltW. D.VanhommerigJ. W.SmuldersN. M.DrossaertsJ. M.KemerinkG. J.BouvyN. D.et al (2009). Cold-activated brown adipose tissue in healthy men. N. Engl. J. Med.360 (15), 1500–1508. 10.1056/NEJMoa0808718

  • CrossRef
  • Google Scholar

• 77

  VialardF.OlivierM. (2020). Thermoneutrality and immunity: how does cold stress affect disease?Front. Immunol.11, 588387. 10.3389/fimmu.2020.588387

  • CrossRef
  • Google Scholar

• 78

  VirtanenK. A.LidellM. E.OravaJ.HeglindM.WestergrenR.NiemiT.et al (2009). Functional brown adipose tissue in healthy adults. N. Engl. J. Med.360 (15), 1518–1525. 10.1056/NEJMoa0808949

  • CrossRef
  • Google Scholar

• 79

  WangB.LiuJ.LeiR.XueB.LiY.TianX.et al (2022). Cold exposure, gut microbiota, and hypertension: a mechanistic study. Sci. Total Environ.833, 155199. 10.1016/j.scitotenv.2022.155199

  • CrossRef
  • Google Scholar

• 80

  WangT. J.LarsonM. G.VasanR. S.ChengS.RheeE. P.McCabeE.et al (2011). Metabolite profiles and the risk of developing diabetes. Nat. Med.17 (4), 448–453. 10.1038/nm.2307

  • CrossRef
  • Google Scholar

• 81

  WangX.CheH.ZhangW.WangJ.KeT.CaoR.et al (2015). Effects of mild chronic intermittent cold exposure on rat organs. Int. J. Biol. Sci.11 (10), 1171–1180. 10.7150/ijbs.12161

  • CrossRef
  • Google Scholar

• 82

  WangZ.NingT.SongA.RutterJ.WangQ. A.JiangL. (2020). Chronic cold exposure enhances glucose oxidation in brown adipose tissue. EMBO Rep.21 (11), e50085. 10.15252/embr.202050085

  • CrossRef
  • Google Scholar

• 83

  WangZ.WangQ. A.LiuY.JiangL. (2021). Energy metabolism in brown adipose tissue. Febs J.288 (12), 3647–3662. 10.1111/febs.16015

  • CrossRef
  • Google Scholar

• 84

  WeirG.RamageL. E.AkyolM.RhodesJ. K.KyleC. J.FletcherA. M.et al (2018). Substantial metabolic activity of human Brown adipose tissue during warm conditions and cold-induced lipolysis of local triglycerides. Cell Metab.27 (6), 1348–1355. 10.1016/j.cmet.2018.04.020

  • CrossRef
  • Google Scholar

• 85

  WorthmannA.JohnC.RühlemannM. C.BaguhlM.HeinsenF. A.SchaltenbergN.et al (2017). Cold-induced conversion of cholesterol to bile acids in mice shapes the gut microbiome and promotes adaptive thermogenesis. Nat. Med.23 (7), 839–849. 10.1038/nm.4357

  • CrossRef
  • Google Scholar

• 86

  YangX.SuiW.ZhangM.DongM.LimS.SekiT.et al (2017). Switching harmful visceral fat to beneficial energy combustion improves metabolic dysfunctions. JCI Insight2 (4), e89044. 10.1172/jci.insight.89044

  • CrossRef
  • Google Scholar

• 87

  YoneshiroT.WangQ.TajimaK.MatsushitaM.MakiH.IgarashiK.et al (2019). BCAA catabolism in brown fat controls energy homeostasis through SLC25A44. Nature572 (7771), 614–619. 10.1038/s41586-019-1503-x

  • CrossRef
  • Google Scholar

• 88

  YuD.RichardsonN. E.GreenC. L.SpicerA. B.MurphyM. E.FloresV.et al (2021). The adverse metabolic effects of branched-chain amino acids are mediated by isoleucine and valine. Cell Metab.33 (5), 905–922. e6. 10.1016/j.cmet.2021.03.025

  • CrossRef
  • Google Scholar

• 89

  YuX. X.LewinD. A.ForrestW.AdamsS. H. (2002). Cold elicits the simultaneous induction of fatty acid synthesis and beta-oxidation in murine brown adipose tissue: prediction from differential gene expression and confirmation in vivo. Faseb J.16 (2), 155–168. 10.1096/fj.01-0568com

  • CrossRef
  • Google Scholar

• 90

  ZechnerR.ZimmermannR.EichmannT. O.KohlweinS. D.HaemmerleG.LassA.et al (2012). FAT SIGNALS-lipases and lipolysis in lipid metabolism and signaling. Cell Metab.15 (3), 279–291. 10.1016/j.cmet.2011.12.018

  • CrossRef
  • Google Scholar