Dietary obesity in mice is associated with lipid deposition and metabolic shifts in the lungs sharing features with the liver

Obesity is associated with both chronic and acute respiratory illnesses, such as asthma, chronic obstructive pulmonary disease (COPD) or increased susceptibility to infectious diseases. Anatomical but also systemic and local metabolic alterations are proposed contributors to the pathophysiology of lung diseases in the context of obesity. To bring perspective to this discussion, we used NMR to compare the obesity-associated metabolomic profiles of the lung with those of the liver, heart, skeletal muscles, kidneys, brain and serum from male C57Bl/6J mice fed with a high-fat and high-sucrose (HFHSD) diet vs. standard (SD) chow for 14 weeks. Our results showed that the lung was the second most affected organ after the liver, and that the two organs shared reduced one-carbon (1C) metabolism and increased lipid accumulation. Altered 1C metabolism was found in all organs and in the serum, but serine levels were increased only in the lung of HFHSD compared to SD. Lastly, tricarboxylic acid (TCA)-derived metabolites were specifically and oppositely regulated in the serum and kidneys but not in other organs. Collectively, our data highlighted that HFHSD induced specific metabolic changes in all organs, the lung being the second most affected organ, the main alterations affecting metabolite concentrations of the 1C pathway and, to a minor extend, TCA. The absolute metabolite quantification performed in this study reveals some metabolic specificities affecting both the liver and the lung, that may reveal common metabolic determinants to the ongoing pathological process.

Obesity has reached epidemic proportions worldwide and affects the metabolism of all organs, impacting therefore whole-body homeostasis. As a consequence, obesity increases the risk of important chronic pathologies such as type 2 diabetes and cardiovascular diseases 1 . Whereas the liver, skeletal muscle, adipose tissue and the heart are classically studied in the context of obesity, others required more attention, such as the lung. Indeed, epidemiological studies have long shown that obesity is associated with both chronic and acute respiratory complications such as asthma, infectious diseases or chronic obstructive pulmonary disease (COPD) [2][3][4][5][6] . Therefore, the increased prevalence of obesity and the concomitant lung diseases has raised the need for a precise examination of obesity-induced alterations in lung metabolism 7 . This topic has gained outstanding relevance with the Covid-19 (SARS-CoV-2) pandemic and its severe pulmonary complications that particularly affect obese patients 8 .
The association between obesity and lung illnesses has long been discussed as the mechanical consequence of obesity on lung physiology and the systemic inflammation induced by dysfunctional adipose tissue 9 . Recent studies evidenced general and lung-localized metabolic alterations that may contribute to the pathophysiology of lung diseases in obesity. Notably, alterations of the lipid metabolism and the secretion of cytokines, or changes of the arginine metabolism leading to nitric oxide synthesis, have been proposed [2][3][4][10][11][12] . At the cellular level, impaired mitochondrial metabolism has also been observed 13,14 , and stressed mitochondria have been reported in the airway epithelial cells of high-fat and high-fructose fed mice 15 .
To further explore and characterize the array of specific alterations in lung metabolism induced by obesity, we used proton nuclear magnetic resonance (NMR) spectroscopy to analyze the metabolic profile induced in mice by an obesogenic diet, and we compared the metabolomes of lung to that of serum and five other tissues (liver, heart, skeletal muscle, kidneys and brain). Our NMR study enabled the quantification of 44 metabolites in organs and 32 metabolites in serum. It evidenced that the lung is one of the most metabolically affected organs, sharing specificities with the liver, and revealed unexpected organ-specific alterations of the one-carbon (1C) and tricarboxylic acid (TCA) pathways.

HFHSD-induced obesity in mice.
Metabolomic studies were performed in male C57Bl/6J mice fed either with a standard diet (SD) or a high-fat and high sucrose diet (HFHSD) for 14 weeks, a classical nutritional model of obesity 16 . As expected, HFHSD mice gained 65% more weight than SD animals (p < 0.0001, n = 12 per group) with a 57% increase of the weight of the liver (p < 0.0001) but not of the lung (Table 1). HFHSD mice showed significant hyperglycemia compared to SD mice (p < 0.0001, Table 1). As shown on Supplementary Fig. 1, glucose tolerance test (GTT) evidenced glucose intolerance with a 92% higher area under the curve of the glycemic response in HFHSD-fed mice compared to SD mice (p < 0.0001, Table 1, n = 12 per group). In agreement, insulin tolerance test (ITT) highlighted systemic insulin resistance, with a 38% lower area over the curve of the glycemic response in HFHSD-fed animals compared to controls (p < 0.01, Table 1, n = 12 per group). Altogether, these results confirm that the 14-weeks HFHSD diet induced significant alterations in body composition and systemic glucose homeostasis in mice, in favor with a progression towards type 2 diabetes 17 .
NMR determination of organ and serum metabolic profiles. We analyzed the metabolome of the lung by NMR spectroscopy and compared it with the metabolome of serum and five other organs (liver, heart, skeletal muscle, kidneys and brain) in both SD and HFHSD mice (n = 8 per group). This approach allowed the absolute quantification of several metabolites, allowing an inter-organ comparison of HFHSD-induced metabolic alterations. All serum and organ extracts delivered interpretable NMR spectra that displayed sharp peaks typical of small molecules, overlaid with broad signals originating mostly from lipids in organ extracts, or lipids plus proteins in serum. The careful analysis of the 1 H 1D and the two-dimensional 1 H-1 H and 1 H-13 C NMR led to the identification and quantification of 44 metabolites in organs and 32 metabolites in serum. The relative lipid content could be also evaluated 18 and results were confirmed by triglyceride (TG) biochemical quantification (see below).
Obesity-associated metabolic signature of the lung. Among the 44 metabolites identified in all organs, 32 metabolites have been quantified in lung extracts (Table 2, n = 8 per group). NMR analysis of the lung showed that Table 1. Physiological parameters of SD and HFHSD mice. Results are expressed as mean ± standard deviation (SD). Data are mean of n = 12 mice per group for total weight, glycaemia, GTT and ITT (glucose and insulin resistance test respectively) and n = 8 mice per group for liver and lung weight; Results of unpaired two-tailed t-test are indicated as non-significant (ns), p < 0.01 (**) and p < 0.0001 (****). Percentages of variation between HFHSD and control fed mice are indicated. HFHSD induced significant changes in 9 metabolites compared to SD mice (Fig. 1). Betaine (trimethylglycine) concentration decreased by 66%, from 47 ± 13.2 nmol/100 mg of tissue to 16.1 ± 5.9 nmol/100 mg (p < 0.0001) while serine was increased by 38%, from 26.2 ± 6.6 nmol/100 mg to 36.3 ± 5.2 nmol/100 mg (p < 0.01). Sn-glycero-3-phosphocholine, glutamine, valine and taurine increased by 32% (p < 0.01), 23% (p < 0.05), 19% (p < 0.05), 18% (p < 0.05) respectively, while myo-inositol, phenylalanine and o-phosphocholine were decreased by 28% Table 2. Metabolite concentrations measured by NMR in the lung of SD and HFHSD mice. Results are expressed as mean ± standard deviation (SD) and are in nmol/100 mg of lung. Data are mean of n = 8 mice per group; Results of unpaired two-tailed t-test are indicated as non-determined (ND), non-significant (NS), p < 0.05 (*), p < 0.01 (**), p < 0.001 (***) and p < 0.0001 (****). Percentages of variation between HFHSD and control fed mice are indicated. www.nature.com/scientificreports/ (p < 0.001), 18% (p < 0.05) and 18% (p < 0.05) respectively, in HFHSD lung compared to SD lung (Fig. 1b). Interestingly, we noted the presence of broad peaks on lung NMR spectra, that could be assigned to the CH3, CH2 and CH2-CO moieties of heterogeneous lipid species 18 . Their peak area was illustrated and quantified in Fig. 2 (panel a and b, respectively). To confirm these results, we decided to measure biochemically TG levels in the lung of SD and HFHSD mice (Fig. 2c). TG levels were 94% higher in lungs of HFHSD-fed animals compared to SD controls, averaging 736 ± 181 μg/100 mg and 380 ± 101 μg/100 mg respectively (p < 0.001), confirming lipid accumulation in the lung of HFHSD mice. List of metabolites differentially found in the lung of HFHSD mice compared to SD mice. Results are expressed as mean ± standard deviation (SD) and are in nmol/100 mg of lung. Data are mean of n = 8 mice per group; Results of unpaired two-tailed t-test are indicated p < 0.05 (*), p < 0.01 (**), p < 0.001 (***) and p < 0.0001 (****). Percentages of variation between HFHSD and SD fed mice are indicated. All metabolites quantified including non-significant variation are reported in Table 2.  Table 3, n = 8 per group), in agreement with the glycemia measurement using a glucometer (Table 1). NMR analyses also evidenced lower serum concentrations for allantoin, creatine and phenylalanine, and higher concentrations for alanine and 3-hydroxybutyrate in HFHSD-fed mice compared to controls (p < 0.05, Table 3). Citrate and succinate concentrations, two TCA-derived metabolites, were 28% and 32% higher in the serum of HFHSD-fed mice compared to controls, respectively (p < 0.02). Inversely, glycine, dimethylglycine (DMG) and betaine that derive from the 1C cycle metabolism were 34%, 68% and 54% lower in the serum of HFHSD-fed mice compared to controls, respectively (p < 0.01). No significant alteration of branchedchain amino acid concentrations was observed in serum (Table 3).
Some metabolite alterations are shared between two or more organs. For example, the HFHSD induced a significant rise of glucose concentration by 28%, 61% and 51% in the liver, heart and kidneys, respectively (p < 0.05, Tables 4, 5 and 7), but not in gastrocnemius muscle. Glucose was not detected in the brain, probably because of its rapid consumption before methanol metabolism quenching (Table 8). Interestingly, the levels of the TCA cycle-derived metabolites, such as citrate, succinate and fumarate, were not significantly altered in organs of HFHSD-fed animals compared to controls, except in kidneys where the citrate concentration was 61% lower Table 3. Metabolite concentrations measured by NMR in the serum of SD and HFHSD mice. Results are expressed as mean ± standard deviation (SD) and are in micromolar in serum. Data are mean of n = 8 mice per group; Results of unpaired two-tailed t-test are indicated as non-determined (ND), non-significant (NS), p < 0.05 (*), p < 0.01 (**), p < 0.001 (***) and p < 0.0001 (****). Percentages of variation between HFHSD and control fed mice are indicated. www.nature.com/scientificreports/ (p < 0.01) ( Table 7). Based on our results in lung and serum of HFHSD mice, we focused our analysis on the 1C pathway, with an emphasis on the methionine and folate cycles to highlight the central role of glycine and its N-methylated analogs that are shared by both cycles (Fig. 3). The liver evidenced the strongest decrease in the methionine cycle with 38% lower concentrations in glycine, 60% in sarcosine, 84% in DMG and 93% in betaine in HFHSD-fed mice compared to controls (p < 0.003, Table 4). Similarly, in gastrocnemius muscle, both glycine Table 4. Metabolite concentrations measured by NMR in the liver of SD and HFHSD mice. Results are expressed as mean ± standard deviation (SD) and are in nmol/100 mg of liver. Data are mean of n = 8 mice per group; Results of unpaired two-tailed t-test are indicated as non-determined (ND), non-significant (NS), p < 0.05 (*), p < 0.01 (**), p < 0.001 (***) and p < 0.0001 (****). Percentages of variation between HFHSD and control fed mice are indicated.

Ctrl (± SD) HFHSD (± SD) p-value % variation
Liver (in nmol/100 mg) www.nature.com/scientificreports/ and DMG concentrations were 30% and 50% lower in HFHSD-fed mice compared to controls, respectively (p < 0.002, Table 6). In kidneys, betaine was 50% lower in HFHSD-fed mice compared to controls (p < 0.0003, Table 7), as observed above in the lung. Globally, all these observations evidenced important variations in the 1C pathway in several tissues of HFHSD-fed mice. Table 5. Metabolite concentrations measured by NMR in the heart of SD and HFHSD mice. Results are expressed as mean ± standard deviation (SD) and are in nmol/100 mg of heart. Data are mean of n = 8 mice per group; Results of unpaired two-tailed t-test are indicated as non-determined (ND), non-significant (NS), p < 0.05 (*), p < 0.01 (**), p < 0.001 (***) and p < 0.0001 (****). Percentages of variation between HFHSD and control fed mice are indicated.

Ctrl (± SD) HFHSD (± SD) p-value % variation
Heart (in nmol/100 mg) www.nature.com/scientificreports/ As observed for the lung, we noted the presence of broad peaks on NMR spectra of some tissues (Fig. 2a,b), that could be assigned to the CH 3 , CH 2 and CH 2 -CO moieties of heterogeneous lipid species. We also quantified the TG levels in all organs. As showed on Fig. 2c, TG levels were 378% higher in the liver of HFHSD mice compared to SD mice, but unaltered in the other organs. The difference between NMR and triglycerides quantification in the skeletal muscles, and to a lesser extent in the heart, may indicate the accumulation of free fatty acids in these organs, which can be detected by NMR only, a phenomenon already observed with high-fat diets 19 . Table 6. Metabolite concentrations measured by NMR in the gastrocnemius muscle of SD and HFHSD mice. Results are expressed as mean ± standard deviation (SD) and are in nmol/100 mg of gastrocnemius muscle. Data are mean of n = 8 mice per group; Results of unpaired two-tailed t-test are indicated as non-determined (ND), non-significant (NS), p < 0.05 (*), p < 0.01 (**), p < 0.001 (***) and p < 0.0001 (****). Percentages of variation between HFHSD and control fed mice are indicated. www.nature.com/scientificreports/ Main comparison between metabolomes of lung, serum and other organs. Obesity-associated metabolome of the lung displayed specific metabolic alterations that are not found in serum and other organs, such as an increase by 36% (p < 0.01) in serine levels ( Fig. 1b and Table 2). It should be noted that serine could not be determined in liver samples due to peak overlap with unidentified peaks on the NMR spectra. We therefore cannot exclude that this modification may be shared with the liver. Table 7. Metabolite concentrations measured by NMR in the kidney of SD and HFHSD mice. Results are expressed as mean ± standard deviation (SD) and are in nmol/100 mg of kidney. Data are mean of n = 8 mice per group; Results of unpaired two-tailed t-test are indicated as non-determined (ND), non-significant (NS), p < 0.05 (*), p < 0.01 (**), p < 0.001 (***) and p < 0.0001 (****). Percentages of variation between HFHSD and control fed mice are indicated. www.nature.com/scientificreports/ Yet, some metabolic signatures of the lung were shared with other organs. For example, a reduction in phenylalanine was found in the lung, serum and liver of HFHSD mice, whereas betaine was significantly decreased in the lung, serum, liver and kidney (Tables 2, 3, 4, 7). In agreement, other metabolites of the 1C pathway, such as glycine and DMG, were reduced in HFHSD mice in serum, liver and skeletal muscle (Tables 3, 4, 6), while Table 8. Metabolite concentrations measured by NMR in the brain of SD and HFHSD mice. Results are expressed as mean ± standard deviation (SD) and are in nmol/100 mg of brain. Data are mean of n = 8 mice per group; Results of unpaired two-tailed t-test are indicated as non-determined (ND), non-significant (NS), p < 0.05 (*), p < 0.01 (**), p < 0.001 (***) and p < 0.0001 (****). Percentages of variation between HFHSD and control fed mice are indicated. www.nature.com/scientificreports/ sarcosine was specifically decreased in the liver of HFHSD compared to SD (Table 4). Altogether, these data point a strong reduction of 1C metabolism with obesity. While the liver remained the most altered organ in HFHSD, it shared with the lung several metabolic features such as decreases in betaine, myo-inositol, o-phosphocholine and phenylalanine and increases in glutamine and taurine (Tables 2, 4). In addition, both tissues shared an accumulation of TG that was not found in other organs (Fig. 2). Finally, TCA metabolism showed organ-specific alterations in HFHSD mice. TCA cycle-derived metabolites were not impacted in the majority of organs including the lung but citrate and succinate concentrations were higher in the serum and citrate was reduced in kidneys of HFHSD mice compared to SD.

Discussion
Obesity is associated with numerous organ dysfunctions, particularly in the lung with increased sensitivity to both chronic and acute illnesses 1-9,20-24 . To get insights into the contribution of cell metabolism alterations to lung fragility, we used NMR to compare the metabolic changes induced by obesity in the lung in comparison to serum and five organs (liver, heart, skeletal muscle, kidneys and brain). Our results demonstrate mainly that the lung displays both specific and shared metabolic changes following HFHSD, compared to serum and/or organs. The obesity-associated metabolic signature of the lung is nearest to the liver rather than other organs, sharing reduction in 1C metabolism and induction of lipid accumulation. Opposite, TCA-derived metabolites are not regulated in the lung and the majority of organs of HFHSD mice, whereas some of them are increased in serum (citrate and succinate) and reduced in kidneys (succinate). Altogether, our data point that obesity (i) mainly affects 1C metabolism in both serum and organs, including the lung, (ii) specifically regulates TCA metabolism in serum and kidney and (iii) specifically increases serine levels in the lung of HFHSD mice.
One important observation of our study is that HFHSD-fed mice presented serum-and organ-specific metabolomics adaptations. Care should be taken as some metabolic adaptations in organs may not only be a consequence of the diet. Indeed, obesity per se, rather than overnutrition, can induce organ-specific metabolic adaptations. This question was raised by Mora-Ortiz et al. who examined db/db mice organs 25 , a natural genetic model of obesity and also reported a decrease of 1C related metabolite in many organs (a glycine decrease in kidneys and distal colon; a homoserine decrease in liver and plasma; a serine decrease in brain, kidney and spleen). Unfortunately, they did not include lung in their study.
One striking observation is that HFHSD-feeding dampened 1C metabolism in all tested samples and specifically increased serine in the lung of HFHSD mice, whereas the TCA-derived metabolites were loosely regulated in several tissues, including the liver and the lung. In the liver, the importance of 1C cycle alterations with dietinduced obesity is supported by another multi-omics analyses 26 , whereas alterations in gene expression and Figure 3. The accumulation of lipids is linked to alterations in the 1C cycle but not to the TCA pathway derived metabolites in the lung of HFHSD mice. Schematic diagram of 1C (methionine plus folate cycles) and TCA cycles, their subcellular compartmentations (mitochondrion in pink) and their potential link with de novo lipogenesis 30,31 . HFHSD causes significant modifications in 1C-derived metabolites in the lung and serum, and a significant modification in TCA-derived metabolites only in serum. Modified metabolites are represented in red, unmodified in green and non-detected in black. 1C-derived metabolites were also altered in other organ and specifically in the liver but were not represented here.  27 . Interestingly, both 1C and TCA pathways could be interrelated and influence each other. Indeed, it is well accepted that under obesogenic conditions, the TCA-cycle is increased and produces citrate that provides carbon atoms for de novo fatty acid (FA) synthesis in the cytosol 28 . This anabolic pathway requires NADPH, which can be produced from different sources, such as the pentose phosphate pathway 29 . Recent studies suggest that the 1C pathway, and in particular serine, contributes in cytosolic NADPH production, as mutations that blocked mitochondrial folate metabolism resulted in fatty acid labeling from deuterated serine [30][31][32] . Therefore, although TCA-derived metabolites were not significantly modified by HFHSD in the lung, the increase in serine levels could sign an increase in de novo lipogenesis contributing to the accumulation of lipids observed in this organ. In agreement, increased serine levels was previously reported in the liver of diabetic rats which also have increased lipogenesis 33 . At the same time, we did not find clear modulation of the TCA-derived metabolites in the liver of HFHSD mice. Liver is known to highly respond to HFHSD and perform de novo lipogenesis 34 , that is fueled from TCA metabolites 28 . Therefore, our results suggest that metabolic fingerprints by NMR may not be able to detect modifications of de novo lipogenesis or TCA cycles. Those are rather detected in serum of HFHSD mice that reflects an homeostatic mechanism to regulate intracellular metabolite levels by secretion or exchange with biofluids (interstitial fluid, serum, urine) 35 . On the opposite, the concentration of metabolites derived from the 1C cycle was decreased in serum and in most organs, except in the heart and brain, indicating that the 1C pathway is highly sensitive and adaptable to the physiological environment and represent a relevant research perspective in the understanding of organ-specific metabolic alterations induced by obesity. Lipid content is increased both in the liver and lung of the HFHSD mice, but not in the other organs considered in this study. Lipid accumulation in the liver is well documented, contrary to the lung. In the liver, lipid accumulation arises from both exogenous lipid storage and de novo lipogenesis using carbons derived from carbohydrates 34 . During a high-fat diet, exogenous fatty acids inhibit de novo lipogenesis and are mainly responsible for steatosis 36,37 . Accumulation of fatty acids in the lung is less understood. Lipid-targeted metabolomics studies by mass spectrometry showed significant quantitative and qualitative differences in neutral lipids, fatty acids, phospholipids and sphingolipids content in lung between HFHSD and very high fat diets 4 , while no lipid accumulation occurred in the lung with SD 38 . As in the liver, the lipid accumulation in the lung correlated with alterations in the levels of many metabolites, with almost 30% of the quantified metabolites affected. Nevertheless, the patterns of metabolic alterations in the lung and liver were distinct. Other lipidomics signature using coupled gas chromatography-mass spectrometry analyses have been reported in the lung following different obesogenic diets, which are dependent of the type of diets 4 , confirming lipotoxicity in the lung with obesity. As the lipotoxicity contributes to dysfunction and pathological states in organs, such as the hepatic lipid accumulation participates to nonalcoholic fatty liver diseases, an involvement of the lung lipid storage is possibly involved in obesity-associated complications such as asthma and COPD, or increased susceptibility to infectious diseases such as seasonal flu or SARS-CoV-2.
In conclusion, our results provide insights on the influence of HFHSD on mice's internal organ metabolism, especially and unexpectedly on the lung metabolism. Our analysis suggests that the TG deposition in the lungs, as in the liver but not in other organs, may be related to adaptations of the 1C metabolism. Mitochondrial metabolism may be central to the obesity-induced metabolic remodeling and its cellular outcomes.

Methods
The study was carried out in compliance with the ARRIVE guidelines 39 . Mice were maintained and all experiments were performed in accordance the French guidelines for the care and use of animals 40 and were approved by a regional ethic committee agreed by the French Ministry for Education and Research and set up by the agreed institution PBES (Ecole Normale Supérieure de Lyon) (APAFIS#4103-2016010417246403).
Animal study. Male C57BL/6J mice (5-weeks-old, n = 12/group) (Envigo, France) were fed with a standard diet (SD, Rhod16A, GENOBIOS) or a high fat-high sucrose diet (HFHSD, 36% w/w fat, 16.6% w/w sucrose, Envigo, France) ad libitum for 14 weeks. Body weight was monitored weekly. Glucose and insulin tolerance tests (GTT and ITT respectively) were performed the last week of feeding on all mice and were performed with a delay of 4 days (meaning that GTT was performed at the beginning of the 13th week and ITT at the end of the week. For GTT and ITT, mice were fasted for 6 or 16 h respectively. Then, glucose (2 mg/g body weight) or insulin (0.75 mU/g body weight) were injected i.p. and blood glucose levels were monitored at the indicated time points. At the end of the feeding protocol, mice at fed state were anesthetized with isoflurane for retro-orbital blood functions and sacrificed by cervical dislocation for removing organs. Blood samples were kept on wet ice for a maximum of 1 h before being centrifuged at 5000g for 10 min in a 4 °C refrigerated centrifuge for serum removal. Aliquots of serum were stored at − 80 °C. Lungs, heart, gastrocnemius muscles, kidneys, brains were harvested, weighed, divided, snap-frozen in liquid nitrogen and stored at − 80 °C.
Biochemical metabolite measurement. For ITT and GTT, glycemia was monitored every 15 min for 90 and 60 min respectively using an Accu-Chek (Roche, France) glucose meter. Organ triglyceride content was determined on total lysate of tissues according to the manufacturer's instruction (Biolabo 87319).
Metabolite extraction and sample preparation for NMR analyses. Metabolites were extracted using 100% methanol and a Precellys Homogeneiser according to the manufacturer's instructions. Full kidneys or 100 mg of all other frozen organs were put in 1 mL for muscles or 500 µL for other organs of 100% methanol and rapidly broken using Precellys homogenizer with ceramic beads (6 per tubes) and 6600 RPM twice 10 s with 10 s pause program at room temperature. The program was repeated again for muscles. The total lysate was put www.nature.com/scientificreports/ into 5 mL glass tubes and the Precellys tubes rinsed with 500 µL of 100% methanol. The glass tubes were hermetically closed and stored at − 20 °C. The extracts were dried under a gentle N2 flow until complete evaporation (approximately 10-12 h), then stored at − 20 °C until NMR sample preparation, right before analysis. 800 μL of D 2 O phosphate buffer prepared as described in Ref. 41 were used to dissolve the dried lysate residues by vortexing for 30 s. Extracts were then transferred to 1.5 mL Eppendorf tubes and centrifuged at 13,000 rpm for 1 min at 4 °C. 550 μL of supernatant were transferred to 5 mm NMR tubes. 100 µL of serum were diluted with 100 µL of D 2 O phosphate buffer and directly transferred to 3 mm NMR tubes. Organ triglyceride content was determined on total lysate according to the manufacturer's instructions (Biolabo 87319).
NMR data acquisition and processing. All organ extracts and serum NMR spectra were obtained on a Bruker Avance III spectrometer operating at a 1H frequency of 800.14 MHz, equipped with a 5 mm TXI probe. The serum samples were measured on a Bruker 600 MHz NMR spectrometer equipped with a 5 mm TCI cryoprobe. All spectra were acquired at 303°K. The NMR samples were maintained at 4 °C before data acquisition and handled by a Bruker SampleJet high-throughput sample changer. Standard NOESY and CPMG with water presaturation 1 H 1D NMR spectra were acquired on each sample, with 128 scans (512 for serum) and a spectral width of 20 ppm. For both sequences, the relaxation delay was set to 4 s. 2D NMR experiments ( 1 H-1 H TOCSY, 1 H-13 C HSQC and J-Resolved) were recorded with standard parameters 41 on a subset of samples to achieve peaks assignment of the metabolites signal. All free induction decays (FIDs) were multiplied by an exponential function leading to a 0.3 Hz line-broadening factor before Fourier transform. 1 H-NMR spectra were manually phased and referenced to the 5.230 glucose doublet using Topspin 3.6 (Bruker GmbH, Rheinstetten, Germany).
NMR spectra analyses. The identification of the metabolites was carried out from the 1 H 1D NMR data using the software ChenomX NMR Suite 8.0 (ChenomX Inc., Edmonton, Canada) and confirmed from analysis of 2D 1 H-1 H TOCSY, 1 H-13 C HSQC and J-Resolved NMR spectra. Metabolite concentrations were determined from the 1 H 1D NOESY experiments using ChenomX software. A pure lactate solution (1 g/L, Fisher) was used as a concentration reference and exploited using the ERETIC2 utility from TopSpin to add a digitally synthesized peak to a spectrum 42  Statistics. All data were expressed as the mean ± SD. Statistical comparisons were performed using Student's t-test.

Data availability
The data used in this study are available upon request. NMR spectra were deposited at the Metabolight database (https:// www. ebi. ac. uk/ metab oligh ts/ MTBLS 2521).