Skip to main content

ORIGINAL RESEARCH article

Front. Vet. Sci., 26 June 2024
Sec. Animal Nutrition and Metabolism
This article is part of the Research Topic Rearing, Health, and Disease Management of Special Economic Animals: Paving the Way for a Sustainable Special Economy through Animal Agriculture View all 10 articles

Lipidomic and transcriptomic profiles provide new insights into the triacylglycerol and glucose handling capacities of the Arctic fox

Yuhang Zhu&#x;Yuhang Zhu1Yuan Yuan&#x;Yuan Yuan2Huazhe Si,Huazhe Si1,3Songze LiSongze Li1Fei ZhaoFei Zhao1Ruina MuRuina Mu1Zihan LinZihan Lin4Xiaoxu WangXiaoxu Wang5Qiang QiuQiang Qiu2Chao Xu
Chao Xu1*Lele Ji
Lele Ji6*Zhipeng Li,
Zhipeng Li1,3*
  • 1College of Animal Science and Technology, Jilin Agricultural University, Changchun, China
  • 2School of Ecology and Environment, Northwestern Polytechnical University, Xi’an, China
  • 3Key Lab of Animal Production, Product Quality and Security, Ministry of Education, Jilin Agricultural University, Changchun, China
  • 4College of Plant Protection, Jilin Agricultural University, Changchun, China
  • 5Department of Special Animal Nutrition and Feed Science, Institute of Special Animal and Plant Sciences, Chinese Academy of Agricultural Sciences, Changchun, China
  • 6National Demonstration Center for Experimental Preclinical Medicine Education, The Fourth Military Medical University, Xi'an, China

The Arctic fox (Vulpes lagopus) is a species indigenous to the Arctic and has developed unique lipid metabolism, but the mechanisms remain unclear. Here, the significantly increased body weight of Arctic foxes was consistent with the significantly increased serum very-low-density lipoprotein (VLDL), and the 40% crude fat diet further increased the Arctic fox body weight. The enhanced body weight gain stems primarily from increased subcutaneous adipose tissue accumulation. The adipose triacylglycerol and phosphatidylethanolamine were significantly greater in Arctic foxes. The adipose fatty-acid synthase content was significantly lower in Arctic foxes, highlighting the main role of exogenous fatty-acids in fat accumulation. Considering the same diet, liver-derived fat dominates adipose expansion in Arctic foxes. Liver transcriptome analysis revealed greater fat and VLDL synthesis in Arctic foxes, consistent with the greater VLDL. Glucose homeostasis wasn’t impacted in Arctic foxes. And the free fatty-acids in adipose, which promote insulin resistance, also did not differ between groups. However, the hepatic glycogen was greater in Arctic foxes and transcriptome analysis revealed upregulated glycogen synthesis, improving glucose homeostasis. These results suggest that the superior fat accumulation capacity and distinct characteristics of hepatic and adipose lipid and glucose metabolism facilitate glucose homeostasis and massive fat accumulation in Arctic foxes.

1 Introduction

The Arctic fox (Vulpes lagopus), the smallest nonhibernating carnivore in Arctic regions, has developed notable adaptations to increase its survival rate and withstand extreme winter temperatures (1). These adaptations include a dietary preference for seabirds, eggs and animal carcasses (2), which leads to a significant increase in body mass through the accumulation of adipose tissue during autumn (3). The increased dietary fat consumption is usually associated with nonalcoholic fatty liver disease (NAFLD), diabetes and hyperlipidemia in humans and rodents (4, 5). However, the intake of a diet containing 40% crude fat does not have adverse effects on hepatic fat accumulation or serum triacylglycerol (TAG) levels in the Arctic fox (6). In contrast, the silver fox (Vulpes vulpes), a species evolutionary close to the Arctic fox, displays a significantly increased serum TAG level upon consumption of a diet containing 26.11% crude fat in comparison with a 14.71% crude fat diet (7). This stark discrepancy suggests that the Arctic fox possesses unique adaptations for lipid metabolism.

To better understand these advantageous characteristics of the Arctic fox, comparative genomic analyses have been performed and have highlighted several positively selected genes related to lipid metabolism and pyrimidine metabolism (8, 9), identified a single-nucleotide polymorphism of Insulin induced gene 2 (Insig2) in the Arctic fox, and increased expression of this gene reduced hepatic lipogenesis in rats (10). Moreover, comparative transcriptomic analysis of the liver, brain and kidney in the Arctic fox revealed that Gltpd1 and Akt2, genes associated with fatty acid metabolism, were under positive selection (11). Akt2, in particular, promotes anabolic lipid metabolism in the liver through the insulin signaling pathway (12). These findings collectively underscore the significance of the genetic basis of lipid metabolism in the Arctic fox and the interconnection between fatty acid and glucose metabolism.

The liver is the central organ that controls lipid and glucose homeostasis (13). TAG can be packaged into very-low-density lipoprotein (VLDL) particles and exported from the liver, with increased VLDL assembly ameliorating NAFLD caused by hepatic lipid droplet accumulation (14, 15). Additionally, the liver maintains blood glucose homeostasis through glycogen synthesis and breakdown (16). Moreover, insulin-like growth factor 1 (IGF-1) produced in the liver not only mediates glycogen synthesis (17) but also functions as a critical hormone for fuel metabolism, increasing glucose tolerance and insulin sensitivity in mice (18). Further evidence for the role of IGF-1 in lipid metabolism includes the significant increase in plasma IGF-1 during the period of rapid weight gain in the reindeer (19) and brown bear during the nonhibernation period when fat depots accumulate (20).

In this study, we investigated the body weight gain rate, serum biochemical parameters, and glucose tolerance of the Arctic fox and silver fox. Moreover, we compared the hepatic gene expression profiles using bulk RNA sequencing and further validated these profiles by determining the concentrations of metabolites or proteins. Finally, we determined the differences in lipidomic component and key enzyme concentrations in adipose tissue (Figure 1A).

Figure 1
www.frontiersin.org

Figure 1. Comparison of serum biochemical analysis and glucose tolerance. (A) The experimental design of this study. (B) Final body weights and body weight gain rate (% of initial) and (C) serum levels of total cholesterol (CHO), very-low-density lipoprotein (VLDL), low-density lipoprotein cholesterol (LDL-C) and high-density lipoprotein cholesterol (HDL-C) among the three groups. The number of animals in the LFS, LFA, and HFA groups were 4, 5, and 5, respectively. (D) Assessment of glucose tolerance between arctic fox (n = 3) and silver fox (n = 3). Results showed as mean ± standard error of the mean (SEM). * indicates significant differences between the LFA group and the LFS group (p < 0.05), # indicates significant differences between the HFA group and the LFA group (p < 0.05). * p < 0.05, ** p < 0.01.

2 Materials and methods

2.1 Animals, experimental design, and sample collection

Sixteen healthy male captive Arctic foxes [body weight (BW) = 5.74 ± 0.10 kg] and 7 healthy male captive silver foxes (BW = 7.08 ± 0.43 kg) that were maintained in the research farm of Jilin Agricultural University were included in this study. The Arctic foxes were randomly assigned to either a 15% crude fat diet group (n = 8, LFA group) or a 40% crude fat diet group (n = 8, HFA group). Silver foxes were fed a 15% crude fat diet (n = 7, LFS group). The dietary components and composition are shown in Supplementary Table S1. Each animal was individually housed, and following a one-week adaptation period, then they were fed each diet for 11 weeks, respectively. The body weight was recorded weekly. All animal procedures were approved and authorized by the Animal Ethics Committee of Jilin Agricultural University (No.20210314002) and used the ARRIVE guidelines 2.0 (21).

At the end of the study, blood samples were collected from the hindlimb vein, and centrifuged at 3,500 g for 10 min at 4°C to obtain serum. Animals were then intravenously injected with 5 mL of succinylcholine (0.4 mg/mL), and subcutaneously adipose tissue and liver tissue were collected. Tissue samples were cut into small pieces (approximately 0.5 × 0.5 cm) and immediately frozen in liquid nitrogen for further analysis. Liver samples were also cut into small pieces (1.0 × 1.0 cm) and fixed in 10% neutral buffered formalin for histological assessment.

2.2 Glucose homeostasis test

The intraperitoneal glucose tolerance test (IPGTT) was performed according to previous established methods (22). Briefly, 3 animals from each group at 9 weeks were selected to conduct the IPGTT. Following a 14 h overnight fast, animals were received intraperitoneal administration glucose (1 g/kg). Whole blood samples were collected from the hindlimb vein before administration (0 min) and at time points (15, 30, 60, 90, 120, 150, and 180 min) after administration, and blood glucose levels were measured using an Accu-Chek glucometer (Roche Diagnostics Shanghai Ltd., China).

2.3 Assessment of metabolites, enzymes, hormones, and liver histology

Serum samples was used to determine the concentrations of glucose, aspartate transaminase (AST), alanine transaminase (ALT), TAG, total cholesterol (CHO), high-density lipoprotein cholesterol (HDL-C) and low-density lipoprotein cholesterol (LDL-C) using the commercial kits (Jiancheng Bioengineering Institute, Nanjing, China) with an automatic biochemical analyzer (Vitalab Selectra E, Vitalab, Dieren, Netherlands). The concentrations of growth hormone (GH), IGF-1, insulin and VLDL in serum were measured using the commercial ELISA kits (MLbio, Shanghai, China).

The liver samples fixed in 10% neutral buffered formalin were embedded in paraffin, cut to 5 μm sections, and stained using H&E (23), and then were visualized using an Olympus AX80 microscope (Olympus Optical, Tokyo, Japan) equipped with Nikon D2X high-resolution camera (Nikon, Tokyo, Japan). The lipid droplet areas were quantified by Image J 1.53a software, and the lipid droplet areas were divided by the area of image to calculate the percentage of lipid droplet areas (24). The concentrations of glucose, TAG, adenosine 3′, 5′-phosphate (cAMP), uridine, glycogen, insulin-like growth factor 1 receptor (IGF1R) and insulin receptor substrate 1 (IRS1) in the liver were measured by the commercial kits (Jiancheng Bioengineering Institute, Nanjing, China) according to the manufacturer’s instructions.

2.4 Liver bulk-RNA sequencing and bioinformatics analysis

Total RNA from liver samples was extracted using the RNeasy Mini Kit (QIAGEN, CA, United States). RNA purity, concentration and integrity were determined using a NanoPhotometer® (IMPLEN, CA, United States), an Agilent 2,100 bioanalyzer in combination with the RNA Nano 6,000 assay kit (Agilent Technologies, CA, United States). A total of 1–3 μg of RNA from each sample with an integrity number greater than 7.0 was used to construct the RNA-Seq library using the VAHTS Universal V6 RNA-seq Library Prep Kit for Illumina®. The cDNA concentration of constructed library was measured using the Qubit® RNA Assay Kit in Qubit® 3.0. Each library was then sequenced on a Novaseq 6000 platform, producing 150 bp paired-end reads.

Trimmomatic (25) was used to remove low quality and adapter sequences, and HISAT2 was employed to align the remaining clean reads to the reference genomes of the Arctic fox and silver fox, respectively (8, 26, 27). Gene expression level was calculated using the raw count data. Differentially expressed genes (DEGs) were determined using the DeSeq2 package (28) based on the threshold |log2 FC| > 0.5 and Benjamini-Hochberg adjusted p value <0.05. The weighted gene coexpression network analysis (WGCNA) was conducted to explore the correlation between DEGs and serum biochemical indices using the WGCNA package (29). The clusterProfiler package was utilized to conduct Gene Ontology (GO) term and Kyoto Encyclopedia of Gene and Genomics (KEGG) pathway enrichment analyses (30).

2.5 Adipose lipidomic analysis and enzyme examination

A total of 20 mg adipose tissue from each animal in the LFA and LFS groups was dissolved in 400 μL water, vortexed for 60 s, and homogenized at 45 Hz for 4 min, and sonicated for 5 min in ice-water bath. The homogenization and sonication circle were repeated for 3 times. A 20 μL portion of the homogenate was mixed with 180 μL water, and then 480 μL extract solution (MTBE: MeOH = 5:1) containing internal standard (Supplementary Table S2) was added. After vortexing for 60 s, the samples were sonicated for 10 min in ice-water bath. The samples were then centrifuged at 3,000 rpm for 15 min at 4°C. A 250 μL supernatant was transferred to a fresh tube. The remaining sample was mixed with 250 μL of MTBE, followed with vortex, sonication and centrifugation, and another 250 μL supernatant was collected. The supernatants were combined and dried in a vacuum concentrator at 37°C and were reconstituted in 200 μL resuspension buffer (DCM: MeOH: H2O = 60: 30: 4.5). The samples were vortexed for 30 s and sonicated for 10 min in ice-water bath. The mixture was then centrifuged at 12,000 rpm for 15 min at 4°C, and 35 μL of supernatant was transferred to a fresh glass vial for LC/MS analysis.

The ultra-high-performance liquid chromatography (UHPLC) separation was performed using a SCIEX ExionLC series UHPLC System. The mobile phase A consisted of 40% water, 60% acetonitrile, and 10 mmol/L ammonium formate. The mobile phase B consisted of 10% acetonitrile and 90% isopropanol, and 10 mmol/L ammonium formate. The column temperature was set to 45°C, while the auto-sampler temperature was set to 6°C. The injection volume was 2 μL. The Biobud-v2.1.4.1 Software (Biotree Biotech Co., Ltd., Shanghai, China) was used to quantify the metabolites. The absolute content of individuals lipids corresponding to the internal standard was calculated based on peak area and the actual concentration of the identical lipid class internal standard. Principal component analysis (PCA) was applied to reveal the differences of lipids. The value of variable importance in the projection (VIP) of the first principal component in orthogonal projections to latent structures discriminant analysis (OPLS-DA) was determined. Metabolites with VIP > 1 and a Benjamini-Hochberg adjusted p value <0.05 (student’s t-test) were considered as significantly changed metabolites (31).

The concentrations of LPCAT4 and FAS were measured by commercial kits from Jiancheng Bioengineering Institute (Nanjing, China) and BiotechPack (Beijing, China), respectively, according to the manufacturer’s instructions.

2.6 Statistical analysis

One-way ANOVA was utilized to determine the statistical significance of final body weights, serum biochemical indices and hormones, areas under curves of IPGTT, liver lipid droplets areas, metabolites and proteins among the LFA, HFA, and LFS groups. The unpaired t-test was performed to identify the statistical significance of adipose metabolites and proteins between the LFS and LFA groups. The statistical analysis was performed using Graphpad Prism (version 9.0.0, GraphPad Software Inc., San Diego, CA, United States). All p value < 0.05 were considered to indicate statistical significance.

3 Results

3.1 Comparison of the growth, biochemical indices, and IPGTT among the three groups

We monitored body weight of Arctic and silver foxes fed a 15% crude fat diet and found that the body weight gain rate of Arctic foxes was significantly greater than that of silver foxes over the 11-week monitoring period. In addition, consumption of a 40% crude fat diet led to further increases in the body weight and the body weight gain rate of Arctic foxes (Figure 1B). Consequently, the thickness of subcutaneous adipose tissue in Arctic foxes was approximately 3 times greater than that in silver foxes (Supplementary Figure S1), indicating that the weight gain observed in Arctic foxes was due primarily to fat deposition. We measured serum biochemical parameters and found that the serum concentrations of CHO, VLDL and HDL-C in Arctic foxes were significantly greater than those in silver foxes fed a 15% crude fat diet (p < 0.05; Figure 1C). Notably, no differences were observed in these serum biochemical parameters between Arctic foxes fed a 15% crude fat diet and those fed a 40% crude fat diet, indicating that Arctic foxes can maintain normal blood lipid levels on a 40% crude fat diet. The serum glucose and insulin levels did not differ significantly among the three experimental groups (Figure 1D). In addition, we conducted an IPGTT and found that the glucose metabolism capacity did not differ significantly among the groups (Figure 1D), suggesting that liver function and glucose homeostasis remained largely unaffected.

3.2 Significant differences in gene expression, the concentrations of protein and metabolite in the liver

Histopathological analysis indicated that even upon consumption of a 40% crude fat diet, the lipid droplet area within the liver remained lower in Arctic foxes than in silver foxes (Figure 2A). Accordingly, consumption of a 40% crude fat diet did not increase the hepatic TAG content in the Arctic fox, nor did it result in elevated serum levels of AST and ALT (Supplementary Figure S2A). Transcriptomic analysis revealed significant differences in global gene expression profiles between LFA and LFS groups as well as between the LFA and HFA groups (p < 0.01; Figure 2B). This analysis led to the identification of DEGs in the following comparisons: LFA versus LFS and HFA versus LFA. A total of 5,504 DEGs were significantly upregulated and 1,719 DEGs were downregulated between the LFA and LFS groups (Figure 2C). In contrast, the comparison between the HFA and LFA groups revealed 30 DEGs, 15 of which were upregulated and 15 of which were downregulated (Supplementary Figure S2B), suggesting that consumption of a 40% crude fat diet does not induce substantial alterations in the liver function of Arctic foxes.

Figure 2
www.frontiersin.org

Figure 2. The results of histological analysis and transcriptome analysis of liver. (A) Liver H&E staining (scale bar, 20 μm) and lipid droplet area in liver (n = 3 in each group). Results showed as mean ± SEM. (B) The PCA plot based on the gene count to reveal the expression profiles in liver among 3 groups. (C) Volcano plot to show the up- and down- differentially expressed genes (DEGs) in the comparison between LFA versus LFS. (D) Kyoto Encyclopedia of Gene and Genomics (KEGG) enrichment analysis of DEGs. (E) Gene Ontology (GO) enrichment analysis of DEGs. GO terms in green, yellow and red regions belonging to biological process, cell component and molecular function, respectively. (F) Correlation between the identified weighted gene coexpression network analysis (WGCNA) modules and the serum biochemical parameters. (G) KEGG and GO enrichment analysis of genes in the module 13. Up, significantly upregulated DEGs; Down, significantly downregulated DEGs. * p < 0.05, ** p < 0.01. GLU, glucose; INS, insulin.

The upregulated DEGs in LFA versus LFS were significantly enriched in the PI3K-Akt signaling pathway (Igf1, Igf1r, Irs1), purine metabolism pathway (Nme2, Cant1, and Adcy6), insulin signaling pathway and pyrimidine metabolism pathway (Nme3, Cant1, and Nt5e) (Figure 2D; Supplementary Figure S2C). These DEGs were also significantly enriched in the following biological process terms: purine-containing compound metabolic process, glucose homeostasis, fatty acid oxidation (Acox1, Acadm, Ppara), pyrimidine-containing compound metabolic process and neutral lipid biosynthetic process (Gpat3, Lpin3, Dgat1, Figure 2E; Supplementary Figure S2C). In addition, insulin-like growth factor binding and insulin receptor substrate binding were noted as enriched molecular function terms (Figure 2E). The downregulated DEGs were enriched in the biosynthesis of cofactors, carbon metabolism, and amino sugar and nucleotide sugar metabolism pathways (Figure 2D); the biological process terms positive regulation of proteolysis and cellular amino acid metabolic process; the cellular component terms organelle outer membrane and sarcoplasmic reticulum; and the molecular function terms structural constituent of ribosome and catalytic activity (Figure 2E).

WGCNA was used to explore the associations between the DEGs and serum biochemical parameters, and 17 modules containing between 45 and 3,084 genes per module were identified (Figure 2F). Module 13 showed a significant positive correlation with the serum concentrations of CHO, HDL-C and VLDL (Figure 2F). The genes within module 13 were enriched in pathways such as ABC transporters, pyrimidine metabolism, PPAR signaling pathway and purine metabolism (Figure 2G). Triglyceride biosynthetic process, pyrimidine-containing compound metabolic process, lipid oxidation, regulation of glucose metabolic process, and response to purine-containing compound were identified as enriched biological process terms. Insulin-like growth factor 1 binding and UDP-xylosyltransferase activity were identified as enriched molecular function terms (Figure 2G). Taken together, these findings suggest that pyrimidine, purine, and glucose metabolism, along with insulin-like growth factor and insulin receptor responses, are likely enhanced in the liver of Arctic foxes. Further analysis revealed significantly greater hepatic concentrations of glycogen, cAMP, uridine, IGF1R and IRS1 in Arctic foxes than in silver foxes (p < 0.05); in contrast, the hepatic concentration of glucose was lower in Arctic foxes than in silver foxes, although the difference was nonsignificant (Figures 3AC). Additionally, the serum concentration of IGF-1 was significantly greater in Arctic foxes than in silver foxes, although the serum concentration of GH was not significantly different between Arctic and silver foxes (Figure 3D).

Figure 3
www.frontiersin.org

Figure 3. Measurement of key metabolites and hormones. (A) Glucose and glycogen concentrations in liver. (B) cAMP and uridine concentrations in liver. (C) Insulin like growth factor 1 receptor (IGF1R) and insulin receptor substrate 1 (IRS1) concentrations in liver among the three groups. (D) The concentrations of insulin-like growth factor 1 (IGF-1) and growth hormone (GH) in serum among the three groups. Results showed as mean ± SEM. * p < 0.05, ** p < 0.01.

3.3 Significantly increased TAG and PE concentrations in the adipose tissue of Arctic foxes

Given the significant changes in body mass, serum cholesterol and liver transcriptome profiles in Arctic foxes (LFA group) compared with silver foxes (LFS group), we explored the lipid composition and concentration in adipose tissue and identified 443 TAGs, 46 diacylglycerols (DAGs) and 26 phosphatidylethanolamines (PEs) (Supplementary Figure S3). TAGs constituted more than 99% of the lipid content in all the samples (Figure 4A). The total concentrations of TAG (p < 0.05) and PE (p < 0.01) in Arctic foxes were significantly greater than those in silver foxes, while the total concentration of free fatty acids (FFAs) was not significantly different (Figure 4B). The PCA results revealed that the lipid profiles were significantly different between the LFA and LFS groups (p < 0.05; Figure 4C). The abundances of 215 lipids were significantly increased in Arctic foxes compared with silver foxes, and the numbers of significantly increased TAGs, PEs and DAGs were 160, 14 and 12, respectively (Figure 4D; Supplementary Tables S3, S4). Further comparison revealed that 8 of the 10 lipids with the greatest log2 fold change (FC) values were TAGs and the other 2 were PEs (Figure 4E), indicating likely increases in TAG and PE synthesis. We thus measured the concentrations of LPCAT4 and FAS that catalyze the synthesis of PEs and fatty acids in adipose tissue. The concentration of LPCAT4 (p < 0.01) in Arctic foxes was significantly greater than that in silver foxes, while the concentration of FAS in Arctic foxes was significantly less than that in silver foxes (p < 0.05) (Figure 4F), suggesting the PE synthesis is increased and fatty acid synthesis is decreased in Arctic foxes.

Figure 4
www.frontiersin.org

Figure 4. Lipidomic analysis and enzymes of adipose tissue. (A) The composition and proportion of the identified lipids in adipose tissue. (B) Differences in total concentrations of TAG, FFA, PE and DAG in adipose tissue. (C,D) PCA and volcano plot revealing the differences of lipid components between arctic fox and silver fox. (E) The top 10 different lipids in adipose tissue between arctic fox and silver fox. The color of the square indicates the significance of the different lipids from maximum (red) to minimum. (F) The concentrations of LPCAT4 and FAS in adipose tissue. Results showed as mean ± SEM. PE, phosphatidylethanolamines; DAG, diacylglycerols; SM, sphingomyelins; CE, cholesterol esters; PC, phosphatidylcholines; CER, ceramides; DCER, dihydroceramides; LPC, lysophosphatidylcholine; PI, phosphatidylglycerol; LPE, lysophosphatidylethanolamines; FFA, free fatty acids. * p < 0.05, ** p < 0.01.

4 Discussion

Our results revealed greater subcutaneous adipose deposition and body weight gain in Arctic foxes when fed a 15% or 40% crude fat diet than in silver foxes when fed a 15% crude fat diet. Previous studies revealed dramatic fat deposition in the Arctic fox and brown bear in autumn (3, 32). These findings confirmed the notable ability for fat accumulation in Arctic animals. We also found that the TAG and PE concentrations in the adipose tissue of Arctic foxes were significantly greater than those in the adipose tissue of silver foxes. This finding is consistent with the observation that the TAG and PE concentrations in adipose tissue of brown bears were greater during hibernation than during the active state (33). When the concentration of TAGs between the leaflets of the endoplasmic reticulum bilayer increases, TAGs coalesce and eventually form lipid droplets (34). It has also been reported that PE is involved in lipid droplet formation and stability (35). Moreover, the concentration of LPCAT4, the key enzyme in PE biosynthesis (36), was also significantly greater in Arctic foxes than in silver foxes. These results suggest that increased PE biosynthesis likely contributes to fat accumulation in Arctic foxes by promoting TAG storage in lipid droplets. We also found that the serum IGF-1 level was significantly higher in Arctic foxes than in silver foxes. This finding is in line with previous findings that plasma IGF-1 is increased in the reindeer and brown bear during the accumulation of fat depots (19, 20). IGF-1 can also induce lipid synthesis through the PI3K-Akt signaling pathway (37, 38). Therefore, the increased amount of IGF-1 is likely responsible for TAG and PE synthesis in the adipose tissue of Arctic foxes. Hepatic synthesis and release of IGF-1 are primarily affected by GH (39). However, the serum GH level was not significantly different between Arctic foxes and silver foxes. The expression of Igf1r in the liver regulates the serum IGF-1 level through a negative feedback mechanism (40). Our results showed that Igf1r was significantly upregulated in Arctic foxes and that the IGF1R concentration in the liver was significantly greater in Arctic foxes than in silver foxes. Therefore, the increased IGF-1 production in Arctic foxes is likely associated with IGF1R-mediated negative feedback.

Interestingly, we found that the FFA content in adipose tissue, the insulin level in serum and glucose tolerance were not significantly different between Arctic foxes and silver foxes. The brown bear also remains insulin sensitive during the fat accumulation period (32). Our results thus indicated the key role of FFA metabolism homeostasis in maintaining insulin sensitivity in the Arctic fox. However, the concentration of the crucial enzyme of de novo lipogenesis, FAS (41), was significantly lower in Arctic foxes than in silver foxes, leading to decreased de novo lipogenesis. Dietary and liver-derived FFAs can also be transported into adipose tissue. FFAs can be esterified by diglyceride acyltransferase on the surface of lipid droplets and stored in lipid droplets in the form of TAGs (42). The TAG concentration in adipose tissue was significantly greater in Arctic foxes than in silver foxes; considering the consumption of the same diet and the decreased de novo lipogenesis in Arctic foxes, liver-derived FFAs could be the main source for TAG accumulation in Arctic foxes. Therefore, the esterification of FFAs to TAGs likely plays a role in maintaining FFA homeostasis, and liver-derived FFAs can dominate fat accumulation in the adipose tissue of Arctic foxes. Moreover, a previous study demonstrated that IGF-1 can promote adipogenesis (43), resulting in increased lipid storage and decreased insulin resistance (44). Therefore, the increased IGF-1 expression in Arctic foxes also likely promotes lipid accumulation in adipose tissue and prevents insulin resistance.

Lipid droplets are ubiquitous intracellular storage organelles specialized for the storage of excess energy in the form of neutral lipids. Surprisingly, we found that the lipid droplet area and TAG concentration in the livers of Arctic foxes were less than those in the livers of silver foxes. However, body weight gain is usually associated with an increased risk of abnormal lipid accumulation in the liver (45). Liver transcriptome analysis revealed greater differences in lipid metabolism-related genes than in other types of genes between Arctic foxes and silver foxes, and the genes upregulated in Arctic foxes were enriched in the neutral lipid biosynthesis process and pyrimidine metabolism, which involves Dgat1 and Nt5e. Moreover, the genes positively correlated with the serum concentrations of CHO, HDL-C and VLDL were significantly enriched in the triglyceride biosynthetic process and pyrimidine metabolism. These findings are consistent with those of a comparative genomic analysis revealing that the positively selected genes in the Arctic fox were also significantly enriched for lipid metabolism and pyrimidine metabolism (8). The lipid metabolism and pyrimidine metabolism pathways have been reported to play important roles in hepatic lipid metabolism (46). Dgat1 encodes a crucial enzyme involved in TAG synthesis (47), and Nt5e is essential for uridine synthesis during pyrimidine metabolism (48). These results indicate the crucial roles of lipid and pyrimidine metabolism in hepatic lipid homeostasis in Arctic foxes. The uridine concentration in the liver and the VLDL concentration in the serum of Arctic foxes were significantly greater than those in silver foxes. TAGs synthesized in the endoplasmic reticulum can be used for the synthesis of lipid droplets or VLDL (49). An increase in VLDL was reported to ameliorate lipid accumulation in liver (14). Previous studies revealed a role for uridine in inhibiting lipid droplet synthesis in the liver (50). Dgat1 is also involved in VLDL synthesis and can increase the plasma VLDL concentration, while Dgat2 (downregulated in Arctic foxes) is associated with liver steatosis (51). These results revealed that the Arctic fox may prevent abnormal lipid droplet accumulation by increasing uridine and VLDL synthesis.

Purine metabolism and fatty acid oxidation, which involve Adcy6, Ppara, Acadm, and Acox1, were also identified as significantly enriched in the liver of Arctic foxes. The proliferator-activated receptor alpha (PPARα) protein encoded by Ppara promotes the expression of Acadm and Acox1, inducing fatty acid oxidation in mitochondria and peroxisomes (52, 53). cAMP, a product of purine metabolism, can activate PPARα through the cAMP/PKA pathway (54), and Adcy6 is the crucial gene for cAMP production (55). Moreover, the liver cAMP concentration was significantly greater in Arctic foxes than in silver foxes. These results reveal that fatty acid oxidation is likely promoted in the Arctic fox through the cAMP/PKA pathway. The PI3K-AKT signaling pathway, insulin signaling pathway and the expression of Slc2a4, which encodes glucose transporter type 4 (GLUT4), were significantly upregulated in Arctic foxes. The PI3K-AKT signaling pathway and insulin signaling pathway are important for glucose homeostasis (56) and GLUT4 transports glucose into the liver (57). However, the liver glycogen concentration in Arctic foxes was significantly greater than that in silver foxes, while the glucose concentration was not significantly different between the two species. The key gene involved in glycogen synthesis, Gys1 (58), was significantly upregulated in Arctic foxes. In addition, IGF-1 can promote Gys1 expression through the PI3K-AKT signaling pathway (59). These results indicate that the Arctic fox may maintain glucose homeostasis by increasing liver glycogen synthesis.

In this study, we confirmed that the Arctic fox adapts to its diet through unique alterations in lipid and glucose metabolism in liver tissue (Figure 5). The Arctic fox can increase TAG storage from liver-derived fat by increasing PE synthesis in adipose tissue. Additionally, it can prevent aberrant lipid droplet accumulation in the liver by increasing uridine and VLDL synthesis and cAMP/PKA pathway-mediated fatty acid oxidation. Moreover, the Arctic fox can maintain glucose homeostasis by increasing liver glycogen synthesis.

Figure 5
www.frontiersin.org

Figure 5. Schematic summary of lipid and glucose metabolism in hepatocyte of arctic fox. The italic red and blue texts indicate up-and down-regulated DEGs in LFA group compared to LFS group. We manually reconstructed the metabolic pathways in liver of Arctic foxes based on the significantly enriched GO terms and KEGG pathways, and the significantly different metabolites, enzymes and hormones. Then the figure was drawn using Adobe Illustrator. UTP, uridine triphosphate; UDP, uridine 5′-diphosphate; UMP, uridine monophosphate; G3P, glycerol 3-phosphate; AGP, 1-Acyl-glycerol 3-phosphate; PA, phosphatidate; DAG, diacylglycerols; TAG, triacylglycerol; LD, lipid droplet; G6P, glucose 6-phosphate; G1P, glucose 1-phosphate; UDPG, uracil-diphosphate glucose; GTP, guanosine 5′-triphosphate; GDP, guanosine 5′-diphosphate; GMP, guanosine 5′-phosphate; IMP, inosine monophosphate; ADS, adenylosuccinate; AMP, adenosine 5′-monophosphate; ADP, adenosine 5′-diphosphate; ATP, adenosine 5′-triphosphate; cAMP, adenosine 3′, 5′-phosphate; PIP3, Phosphatidylinositol-3,4,5-triphosphate.

Data availability statement

The data presented in the study are deposited in the NCBI repository, accession number PRJNA1028137.

Ethics statement

The animal study was approved by Ethical Committee of Jilin Agricultural University. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

YZ: Formal analysis, Investigation, Writing – original draft, Software, Visualization. YY: Formal analysis, Data curation, Software, Writing – review & editing. HS: Investigation, Writing – review & editing. SL: Investigation, Writing – review & editing. FZ: Investigation, Writing – review & editing. RM: Investigation, Writing – review & editing. ZhL: Investigation, Writing – review & editing. XW: Investigation, Writing – review & editing. QQ: Writing – review & editing. CX: Writing – review & editing. LJ: Conceptualization, Methodology, Supervision, Writing – review & editing. ZpL: Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Key Research and Development Program Projects (2023YFD1302000) and Innovation and Entrepreneurship Talent Program from Jilin Province (2023QN02) to ZpL.

Acknowledgments

We would like to thank Qianlong Yang for the assistance on sample collection.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2024.1388532/full#supplementary-material

Supplementary Figure S1 | The thickness of subcutaneous adipose tissue among the LFS, LFA and HFA groups.

Supplementary Figure S2 | Comparison of metabolites, enzymes and gene expression between Arctic foxes and silver foxes. (A) The concentrations of liver TAG and serum AST and ALT among the LFS, LFA and HFA groups. (B) Volcano plot to show the up- and down- differentially expressed genes (DEGs) in liver between HFA and LFA. (C) Heatmap showing the DEGs in liver between Arctic foxes and silver foxes, which are involved in glucose and lipid metabolism. TAG: triacylglycerol; AST: aspartate transaminase; ALT: alanine transaminase.

Supplementary Figure S3 | The composition of lipids in the subcutaneous adipose tissue based on the lipidomic analysis. The percentage in the brackets of pie chart means the proportion of lipids, and the numbers in the brackets of legend means the numbers of lipids. CE: cholesterol esters; CER: ceramides; DAG: diacylglycerols; DCER: dihydroceramides; DGDG: digalactosyldiacylglycerol; FFA: free fatty acids; HCER: hexosylceramides; LPC: lysophosphatidylcholines; LPE: lysophosphatidylethanolamines; PC: phosphatidylcholines; PE: phosphatidylethanolamines; PI: phosphatidylinositol; SM: sphingomyelins; TAG: triacylglycerols.

Supplementary TABLE S1 | Ingredients of diets (DM basis).

Supplementary TABLE S2 | The information of internal standard in lipidomic analysis.

Supplementary TABLE S3 | The identified lipids in the subcutaneous adipose tissue of the LFA and LFS groups.

Supplementary TABLE S4 | Significantly different lipids between the LFA and LFS groups.

References

1. Fuglei, E, and Øritsland, NA. Seasonal trends in body mass, food intake and resting metabolic rate, and induction of metabolic depression in arctic foxes (Alopex Lagopus) at Svalbard. J Comp Physiol B. (1999) 169:361–9. doi: 10.1007/s003600050232

PubMed Abstract | Crossref Full Text | Google Scholar

2. Frafjord, K. Food habits of arctic foxes (Alopex Lagopus) on the western coast of Svalbard. Arctic. (1993) 46:49–54. doi: 10.14430/arctic1321

Crossref Full Text | Google Scholar

3. Prestrud, P, and Nilssen, K. Fat deposition and seasonal variation in body composition of Arctic foxes in Svalbard. J Wildl Manag. (1992) 56:221–33. doi: 10.2307/3808816

Crossref Full Text | Google Scholar

4. Gérard, P. Gut microbiota and obesity. Cell Mol Life Sci. (2016) 73:147–62. doi: 10.1007/s00018-015-2061-5

PubMed Abstract | Crossref Full Text | Google Scholar

5. Ong, JP, and Younossi, ZM. Epidemiology and natural history of NAFLD and NASH. Clin Liver Dis. (2007) 11:1–16. doi: 10.1016/j.cld.2007.02.009

Crossref Full Text | Google Scholar

6. Geng, Y, Zhang, T, Zhang, Z, Gao, X, Yang, F, and Xing, X. Dietary fat levels affect growth performance, body fat deposition and serum biochemical parameters of growing blue foxes (Alopex Lagopus). Chin J Anim Nutr. (2011) 23:1637–46.

Google Scholar

7. Zhang, T, Zhong, W, Sun, WL, Wang, Z, Sun, H, Fan, Y, et al. Effects of dietary fat: carbohydrate ratio on nutrient digestibility, serum parameters, and production performance in male silver foxes (Vulpes Vulpes) during the winter fur-growing period. Can J Anim Sci. (2016) 97:199–206. doi: 10.1139/cjas-2015-0167

Crossref Full Text | Google Scholar

8. Peng, Y, Li, H, Liu, Z, Zhang, C, Li, K, Gong, Y, et al. Chromosome-level genome assembly of the Arctic fox (Vulpes Lagopus) using pacbio sequencing and hi-c technology. Mol Ecol Resour. (2021) 21:2093–108. doi: 10.1111/1755-0998.13397

PubMed Abstract | Crossref Full Text | Google Scholar

9. Grzes, M, Szczerbal, I, Fijak-Nowak, H, Szydlowski, M, and Switonski, M. Two candidate genes (Fto and Insig2) for fat accumulation in four canids: chromosome mapping, gene polymorphisms and association studies of body and skin weight of red foxes. Cytogenet Genome Res. (2011) 135:25–32. doi: 10.1159/000330457

PubMed Abstract | Crossref Full Text | Google Scholar

10. Takaishi, K, Duplomb, L, Wang, M-Y, Li, J, and Unger, RH. Hepatic Insig-1 or −2 overexpression reduces lipogenesis in obese zucker diabetic fatty rats and in fasted/refed normal rats. Proc Natl Acad Sci USA. (2004) 101:7106–11. doi: 10.1073/pnas.0401715101

PubMed Abstract | Crossref Full Text | Google Scholar

11. Kumar, V, Kutschera, VE, Nilsson, MA, and Janke, A. Genetic signatures of adaptation revealed from transcriptome sequencing of Arctic and red foxes. BMC Genomics. (2015) 16:585. doi: 10.1186/s12864-015-1724-9

PubMed Abstract | Crossref Full Text | Google Scholar

12. Wan, M, Leavens, KF, Saleh, D, Easton, RM, Guertin, DA, Peterson, TR, et al. Postprandial hepatic lipid metabolism requires signaling through Akt2 independent of the transcription factors Foxa2, Foxo1, and Srebp1c. Cell Metab. (2011) 14:516–27. doi: 10.1016/j.cmet.2011.09.001

PubMed Abstract | Crossref Full Text | Google Scholar

13. Zhang, H, Ren, G, and Qiao, A. Editorial: hepatic glucose and lipid metabolism. Front Physiol. (2022) 13:1009566. doi: 10.3389/fphys.2022.1009566

PubMed Abstract | Crossref Full Text | Google Scholar

14. Zhu, X, Xiong, T, Liu, P, Guo, X, Xiao, L, Zhou, F, et al. Quercetin ameliorates HFD-induced NAFLD by promoting hepatic VLDL assembly and lipophagy via the Ire1a/Xbp1s pathway. Food Chem Toxicol. (2018) 114:52–60. doi: 10.1016/j.fct.2018.02.019

PubMed Abstract | Crossref Full Text | Google Scholar

15. Rui, L. Energy metabolism in the liver. Compr Physiol. (2014) 4:177–97. doi: 10.1002/cphy.c130024

PubMed Abstract | Crossref Full Text | Google Scholar

16. Sharabi, K, Tavares, CD, Rines, AK, and Puigserver, P. Molecular pathophysiology of hepatic glucose production. Mol Asp Med. (2015) 46:21–33. doi: 10.1016/j.mam.2015.09.003

PubMed Abstract | Crossref Full Text | Google Scholar

17. Wang, J, Zhu, X, Chen, C, Li, X, Gao, Y, Li, P, et al. Effect of insulin-like growth factor-1 (IGF-1) on the gluconeogenesis in calf hepatocytes cultured in vitro. Mol Cell Biochem. (2012) 362:87–91. doi: 10.1007/s11010-011-1130-9

PubMed Abstract | Crossref Full Text | Google Scholar

18. Hong, H, Cui, Z-Z, Zhu, L, Fu, S-P, Rossi, M, Cui, Y-H, et al. Central Igf1 improves glucose tolerance and insulin sensitivity in mice. Nutr Diabetes. (2017) 7:2. doi: 10.1038/s41387-017-0002-0

PubMed Abstract | Crossref Full Text | Google Scholar

19. Bubenik, GA, Schams, D, White, RG, Rowell, J, Blake, J, and Bartos, L. Seasonal levels of metabolic hormones and substrates in male and female reindeer (Rangifer Tarandus). Comp Biochem Physiol C Pharmacol Toxicol Endocrinol. (1998) 120:307–15. doi: 10.1016/S0742-8413(98)10010-5

PubMed Abstract | Crossref Full Text | Google Scholar

20. Frøbert, AM, Brohus, M, Roesen, TS, Kindberg, J, Fröbert, O, Conover, CA, et al. Circulating insulin-like growth factor system adaptations in hibernating brown bears indicate increased tissue igf availability. Am J Physiol Endocrinol Metab. (2022) 323:E307–18. doi: 10.1152/ajpendo.00429.2021

PubMed Abstract | Crossref Full Text | Google Scholar

21. Percie, S, Nathalie, A, Amrita, A, Sabina, A, Marc, TB, Monya, B, et al. Reporting animal research: explanation and elaboration for the ARRIVE guidelines 2.0. PLoS Biol. (2020) 18:e3000411. doi: 10.1371/journal.pbio.3000411

PubMed Abstract | Crossref Full Text | Google Scholar

22. Schneider, S, Feilen, PJ, Brunnenmeier, F, Minnemann, T, Zimmermann, H, Zimmermann, U, et al. Long-term graft function of adult rat and human islets encapsulated in novel alginate-based microcapsules after transplantation in immunocompetent diabetic mice. Diabetes. (2005) 54:687–93. doi: 10.2337/diabetes.54.3.687

PubMed Abstract | Crossref Full Text | Google Scholar

23. Baptista, PM, Siddiqui, MM, Lozier, G, Rodriguez, SR, Atala, A, and Soker, S. The use of whole organ decellularization for the generation of a vascularized liver organoid. Hepatology. (2011) 53:604–17. doi: 10.1002/hep.24067

PubMed Abstract | Crossref Full Text | Google Scholar

24. Abràmoff, MD, Magalhães, PJ, and Ram, SJ. Image processing with image j. Biophoton Int. (2004) 11:36–42.

Google Scholar

25. Bolger, AM, Lohse, M, and Usadel, B. Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics. (2014) 30:2114–20. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | Crossref Full Text | Google Scholar

26. Kim, D, Langmead, B, and Salzberg, SL. Hisat: a fast spliced aligner with low memory requirements. Nat Methods. (2015) 12:357–60. doi: 10.1038/nmeth.3317

PubMed Abstract | Crossref Full Text | Google Scholar

27. Kukekova, AV, Johnson, JL, Xiang, X, Feng, S, Liu, S, Rando, HM, et al. Red fox genome assembly identifies genomic regions associated with tame and aggressive behaviours. Nat Ecol Evol. (2018) 2:1479–91. doi: 10.1038/s41559-018-0611-6

PubMed Abstract | Crossref Full Text | Google Scholar

28. Love, MI, Huber, W, and Anders, S. Moderated estimation of fold change and dispersion for RNA-Seq data with Deseq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | Crossref Full Text | Google Scholar

29. Langfelder, P, and Horvath, S. Wgcna: an R package for weighted correlation network analysis. BMC Bioinformatics. (2008) 9:559. doi: 10.1186/1471-2105-9-559

PubMed Abstract | Crossref Full Text | Google Scholar

30. Yu, G, Wang, L-G, Han, Y, and He, Q-Y. Clusterprofiler: an R package for comparing biological themes among gene clusters. OMICS. (2012) 16:284–7. doi: 10.1089/omi.2011.0118

PubMed Abstract | Crossref Full Text | Google Scholar

31. Chen, H, Cao, G, Chen, D-Q, Wang, M, Vaziri, ND, Zhang, Z-H, et al. Metabolomics insights into activated redox signaling and lipid metabolism dysfunction in chronic kidney disease progression. Redox Biol. (2016) 10:168–78. doi: 10.1016/j.redox.2016.09.014

PubMed Abstract | Crossref Full Text | Google Scholar

32. Rigano, KS, Gehring, JL, Evans Hutzenbiler, BD, Chen, AV, Nelson, OL, Vella, CA, et al. Life in the fat lane: seasonal regulation of insulin sensitivity, food intake, and adipose biology in brown bears. J Comp Physiol B. (2017) 187:649–76. doi: 10.1007/s00360-016-1050-9

PubMed Abstract | Crossref Full Text | Google Scholar

33. Giroud, S, Chery, I, Bertile, F, Bertrand-Michel, J, Tascher, G, Gauquelin-Koch, G, et al. Lipidomics reveals seasonal shifts in a large-bodied hibernator, the brown bear. Front Physiol. (2019) 10:389. doi: 10.3389/fphys.2019.00389

PubMed Abstract | Crossref Full Text | Google Scholar

34. Olzmann, JA, and Carvalho, P. Dynamics and functions of lipid droplets. Nat Rev Mol Cell Biol. (2019) 20:137–55. doi: 10.1038/s41580-018-0085-z

PubMed Abstract | Crossref Full Text | Google Scholar

35. Hörl, G, Wagner, A, Cole, LK, Malli, R, Reicher, H, Kotzbeck, P, et al. Sequential synthesis and methylation of phosphatidylethanolamine promote lipid droplet biosynthesis and stability in tissue culture and in vivo. J Biol Chem. (2011) 286:17338–50. doi: 10.1074/jbc.M111.234534

PubMed Abstract | Crossref Full Text | Google Scholar

36. Hishikawa, D, Shindou, H, Kobayashi, S, Nakanishi, H, Taguchi, R, and Shimizu, T. Discovery of a lysophospholipid acyltransferase family essential for membrane asymmetry and diversity. Proc Natl Acad Sci USA. (2008) 105:2830–5. doi: 10.1073/pnas.0712245105

PubMed Abstract | Crossref Full Text | Google Scholar

37. Hoxhaj, G, and Manning, BD. The PI3K–AKT network at the interface of oncogenic signalling and cancer metabolism. Nat Rev Cancer. (2020) 20:74–88. doi: 10.1038/s41568-019-0216-7

PubMed Abstract | Crossref Full Text | Google Scholar

38. Ma, X, and Bai, Y. IGF-1 activates the PI3K/AKT signaling pathway via upregulation of secretory clusterin. Mol Med Rep. (2012) 6:1433–7. doi: 10.3892/mmr.2012.1110

PubMed Abstract | Crossref Full Text | Google Scholar

39. Jones, JI, and Clemmons, DR. Insulin-like growth factors and their binding proteins: biological actions. Endocr Rev. (1995) 16:3–4. doi: 10.1210/edrv-16-1-3

Crossref Full Text | Google Scholar

40. Romero, CJ, Ng, Y, Luque, RM, Kineman, RD, Koch, L, Bruning, JC, et al. Targeted deletion of somatotroph insulin-like growth factor-i signaling in a cell-specific knockout mouse model. Mol Endocrinol. (2010) 24:1077–89. doi: 10.1210/me.2009-0393

PubMed Abstract | Crossref Full Text | Google Scholar

41. Jensen-Urstad, APL, and Semenkovich, CF. Fatty acid synthase and liver triglyceride metabolism: housekeeper or messenger? Biochim Biophys Acta. (2012) 1821:747–53. doi: 10.1016/j.bbalip.2011.09.017

PubMed Abstract | Crossref Full Text | Google Scholar

42. Chitraju, C, Mejhert, N, Haas, JT, Diaz-Ramirez, LG, Grueter, CA, Imbriglio, JE, et al. Triglyceride synthesis by Dgat1 protects adipocytes from lipid-induced ER stress during lipolysis. Cell Metab. (2017) 26:407–418.e3. doi: 10.1016/j.cmet.2017.07.012

PubMed Abstract | Crossref Full Text | Google Scholar

43. Wabitsch, M, Hauner, H, Heinze, E, and Teller, WM. The role of growth hormone/insulin-like growth factors in adipocyte differentiation. Metabolism. (1995) 44:45–9. doi: 10.1016/0026-0495(95)90220-1

PubMed Abstract | Crossref Full Text | Google Scholar

44. Medina-Gomez, G, Gray, S, and Vidal-Puig, A. Adipogenesis and lipotoxicity: role of peroxisome proliferator-activated receptor γ (PPARγ) and PPARγcoactivator-1 (PGC1). Public Health Nutr. (2007) 10:1132–7. Epub 2007/10/01. doi: 10.1017/S1368980007000614

PubMed Abstract | Crossref Full Text | Google Scholar

45. Cho, EJ, Yu, SJ, Jung, GC, Kwak, M-S, Yang, JI, Yim, JY, et al. Body weight gain rather than body weight variability is associated with increased risk of nonalcoholic fatty liver disease. Sci Rep. (2021) 11:14428. doi: 10.1038/s41598-021-93883-5

PubMed Abstract | Crossref Full Text | Google Scholar

46. Huang, X, Yuan, Z, Liu, X, Wang, Z, Lu, J, Wu, L, et al. Integrative multi-omics unravels the amelioration effects of zanthoxylum bungeanum maxim on non-alcoholic fatty liver disease. Phytomedicine. (2023) 109:154576. doi: 10.1016/j.phymed.2022.154576

PubMed Abstract | Crossref Full Text | Google Scholar

47. Harris, CA, Haas, JT, Streeper, RS, Stone, SJ, Kumari, M, Yang, K, et al. Dgat enzymes are required for triacylglycerol synthesis and lipid droplets in adipocytes. J Lipid Res. (2011) 52:657–67. doi: 10.1194/jlr.M013003

PubMed Abstract | Crossref Full Text | Google Scholar

48. Che, X, Liu, P, Wu, C, Song, W, An, N, Yu, L, et al. Potential role of the ecto-5′-nucleotidase in morphine-induced uridine release and neurobehavioral changes. Neuropharmacology. (2018) 141:1–10. doi: 10.1016/j.neuropharm.2018.07.035

PubMed Abstract | Crossref Full Text | Google Scholar

49. Gluchowski, NL, Becuwe, M, Walther, TC, and Farese, RV. Lipid droplets and liver disease: from basic biology to clinical implications. Nat Rev Gastroenterol Hepatol. (2017) 14:343–55. doi: 10.1038/nrgastro.2017.32

PubMed Abstract | Crossref Full Text | Google Scholar

50. Liu, Y, Xie, C, Zhai, Z, Deng, Z-Y, De Jonge, HR, Wu, X, et al. Uridine attenuates obesity, ameliorates hepatic lipid accumulation and modifies the gut microbiota composition in mice fed with a high-fat diet. Food Funct. (2021) 12:1829–40. doi: 10.1039/D0FO02533J

PubMed Abstract | Crossref Full Text | Google Scholar

51. Yamazaki, T, Sasaki, E, Kakinuma, C, Yano, T, Miura, S, and Ezaki, O. Increased very low density lipoprotein secretion and gonadal fat mass in mice overexpressing liver Dgat1. J Biol Chem. (2005) 280:21506–14. doi: 10.1074/jbc.M412989200

PubMed Abstract | Crossref Full Text | Google Scholar

52. Rakhshandehroo, M, Hooiveld, G, Müller, M, and Kersten, S. Comparative analysis of gene regulation by the transcription factor pparα between mouse and human. PLoS One. (2009) 4:e6796. doi: 10.1371/journal.pone.0006796

PubMed Abstract | Crossref Full Text | Google Scholar

53. Gulick, T, Cresci, S, Caira, T, Moore, DD, and Kelly, DP. The peroxisome proliferator-activated receptor regulates mitochondrial fatty acid oxidative enzyme gene expression. Proc Natl Acad Sci USA. (1994) 91:11012–6. doi: 10.1073/pnas.91.23.11012

PubMed Abstract | Crossref Full Text | Google Scholar

54. Pawlak, M, Lefebvre, P, and Staels, B. Molecular mechanism of pparα action and its impact on lipid metabolism, inflammation and fibrosis in non-alcoholic fatty liver disease. J Hepatol. (2015) 62:720–33. doi: 10.1016/j.jhep.2014.10.039

PubMed Abstract | Crossref Full Text | Google Scholar

55. Simko, V, Iuliano, F, Sevcikova, A, Labudova, M, Barathova, M, Radvak, P, et al. Hypoxia induces cancer-associated cAMP/PKA signalling through HIF-mediated transcriptional control of adenylyl cyclases VI and VII. Sci Rep. (2017) 7:10121. doi: 10.1038/s41598-017-09549-8

PubMed Abstract | Crossref Full Text | Google Scholar

56. Gao, Y-F, Zhang, M-N, Wang, T-X, Wu, T-C, Ai, R-D, and Zhang, Z-S. Hypoglycemic effect of d-chiro-inositol in type 2 diabetes mellitus rats through the PI3K/AKT signaling pathway. Mol Cell Endocrinol. (2016) 433:26–34. doi: 10.1016/j.mce.2016.05.013

PubMed Abstract | Crossref Full Text | Google Scholar

57. Zhao, P, Tian, D, Song, G, Ming, Q, Liu, J, Shen, J, et al. Neferine promotes Glut4 expression and fusion with the plasma membrane to induce glucose uptake in L6 cells. Front Pharmacol. (2019) 10:999. doi: 10.3389/fphar.2019.00999

PubMed Abstract | Crossref Full Text | Google Scholar

58. McCorvie, TJ, Loria, PM, Tu, M, Han, S, Shrestha, L, Froese, DS, et al. Molecular basis for the regulation of human glycogen synthase by phosphorylation and glucose-6-phosphate. Nat Struct Mol Biol. (2022) 29:628–38. doi: 10.1038/s41594-022-00799-3

PubMed Abstract | Crossref Full Text | Google Scholar

59. Zhang, Y, Liu, X, Han, L, Gao, X, Liu, E, and Wang, T. Regulation of lipid and glucose homeostasis by mango tree leaf extract is mediated by AMPK and PI3K/AKT signaling pathways. Food Chem. (2013) 141:2896–905. doi: 10.1016/j.foodchem.2013.05.121

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: fatty acid oxidation, insulin-like growth factor 1, lipid droplets, uridine, VLDL, glycogen

Citation: Zhu Y, Yuan Y, Si H, Li S, Zhao F, Mu R, Lin Z, Wang X, Qiu Q, Xu C, Ji L and Li Z (2024) Lipidomic and transcriptomic profiles provide new insights into the triacylglycerol and glucose handling capacities of the Arctic fox. Front. Vet. Sci. 11:1388532. doi: 10.3389/fvets.2024.1388532

Received: 04 April 2024; Accepted: 10 June 2024;
Published: 26 June 2024.

Edited by:

Izhar Hyder Qazi, Shaheed Benazir Bhutto University of Veterinary &Animal Sciences, Pakistan

Reviewed by:

Lihua Wang, Qingdao Agricultural University, China
Grant B. McClelland, McMaster University, Canada

Copyright © 2024 Zhu, Yuan, Si, Li, Zhao, Mu, Lin, Wang, Qiu, Xu, Ji and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhipeng Li, emhwbGljYWFzQDE2My5jb20=; Lele Ji, amlsZWxlQGZtbXUuZWR1LmNu; Chao Xu, eGNqbGF1QDE2My5jb20=

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.