Obesity frequently leads to insulin resistance and the development of hepatic steatosis. To characterize the molecular changes that promote hepatic steatosis, transcriptomics, proteomics, and metabolomics technologies were applied to liver samples from C57BL/6J mice obtained from two independent intervention trials. After 12 wk of high-fat feeding the animals became obese, hyperglycemic, and insulin resistant, had elevated levels of blood cholesterol and VLDL, and developed hepatic steatosis. Nutrigenomic analysis revealed alterations of key metabolites and enzyme transcript levels of hepatic one-carbon metabolism and related pathways. The hepatic oxidative capacity and the lipid milieu were significantly altered, which may play a key role in the development of insulin resistance. Additionally, high choline levels were observed after the high-fat diet. Previous studies have linked choline levels with insulin resistance and hepatic steatosis in conjunction with changes of certain metabolites and enzyme levels of one-carbon metabolism. The present results suggest that the coupling of high levels of choline and low levels of methionine plays an important role in the development of insulin resistance and liver steatosis. In conclusion, the complexities of the alterations induced by high-fat feeding are multifactorial, indicating that the interplay between several metabolic pathways is responsible for the pathological consequences.
- high-fat feeding
- Kennedy pathway
obesity is now considered a major global health problem, and it is well recognized that it plays a central role in insulin resistance, the metabolic syndrome, and type 2 diabetes mellitus (T2DM). Expansion of adipose tissue is associated with changes in the levels of inflammatory cytokines (2, 6, 51) and has most recently been linked as well to endoplasmic reticulum stress (15). Although increased hepatic accumulation of lipids during obesity development is observed in rodents and humans, the underlying mechanisms are still not fully understood.
Diet is an important environmental factor in the development of obesity, and an imbalance of calorie intake and expenditure leads to weight gain. Because of the complex etiology and pathogenesis of obesity, animal models are valuable tools for deepening our understanding of human diseases in a controlled environment. For many years animal models have been used to study diet-induced obesity by feeding high-energy diets such as high-sucrose, high-fat, and cafeteria diets. One of the most commonly used models is the C57BL/6J mouse, which develops many pathological features like those in humans including visceral fatness, insulin resistance, hyperinsulinemia, and hyperlipidemia (19) but also hepatic steatosis in the course of increased triglyceride accumulation.
The amount of triglyceride present in the liver represents a complex interplay between hepatic uptake of free fatty acids from the plasma, de novo fatty acid synthesis, fatty acid oxidation, and triglyceride export as part of very low-density lipoproteins (VLDL). Many studies have shown that high-fat diets increase the expression of genes that promote lipid storage in rodents (29, 38). More recent studies using tracer compounds in vivo have shown that high-fat feeding induced hepatic fatty acid synthesis by chain elongation and subsequent desaturation as opposed to de novo synthesis (33). In rats it was demonstrated that the contribution of de novo lipogenesis to hepatic triglyceride contents is small when animals are fed a high-fat diet (10). To investigate how the progression of liver steatosis and the development of insulin resistance are promoted by fat accumulation, we used a nutrigenomic approach (employing metabolomics, transcriptomics, and proteomics) in an integrated multicenter study. Analysis of these multilayer data sets was hypothesized to enable a more comprehensive understanding of the hepatic alterations caused by fat storage. Biological samples were collected from two intervention studies performed in two centers as part of the collaborative research activities of the European Nutrigenomics Organisation (NuGO; www.nugo.org) Proof of Principle Study Package (NuGO PPS). The NuGO PPS contained three animal studies as described previously (3); here we report findings from the PPS2 and PPS3 studies. Both studies were performed with the same diets and had similar design, sampling, and analysis strategies with standard operation procedures defined by NuGO. Moreover, both studies were designed to have a common time point, namely 12 wk after low- and high-fat feeding. In the present study, analysis of the nutrigenomic data led to the identification of novel links between a high-fat diet and altered one-carbon metabolism.
MATERIALS AND METHODS
Within the federated package of NuGO two animal studies were conducted under very similar conditions (3) at the TNO-Quality of Life site in Leiden and at Wageningen University. All animals received humane care according to the criteria outlined in the Guide for the Care and Use of Laboratory Animals prepared by the National Academy of Sciences and published by the National Institutes of Health (NIH). All animal experiments were approved by an independent institutional ethical committee on animal care and experimentation, PPS2 by the Dierethische Commissie (DEC; Zeist, The Netherlands) with the study number DEC2465 and PPS3 by the Local Committee for Care and Use of Laboratory Animals at Wageningen University with the study number 2007118.b. In both studies male C57BL/6J mice were obtained at the age of 3 wk and housed in wire-topped Macrolon cages with ad libitum access to food and water, and during the intervention phase the same diet and sampling protocol was applied. Two diets were administered to the animals: a high-fat diet and a low-fat diet, in which 45% and 10% of the energy were obtained from fat, respectively. Refined palm oil was the main fat source in the diets. Low-fat and high-fat diet compositions are specified in Supplemental Table S1 and were based on Research Diets formulas (www.researchdiets.com; D12450B, D12451).1 To balance the energy intake, the amount of corn starch was modified accordingly. Body weight and food intake were measured weekly during the intervention phase.
In PPS2, mice were purchased from Charles River (Maastricht, The Netherlands). At 10 wk of age, animals were subjected to a 2-wk run-in on the low-fat diet specified above. At 12 wk of age, the animals were randomized into two groups (n = 47/group): the high-fat and the low-fat group. During experimental diet feeding, tail blood samples were obtained after 5 h of fasting from 8 AM to 1 PM at weeks 0, 4, 8, and 12 to monitor glucose and insulin levels. After 12 wk of dietary treatment seven animals per group were used for monitoring the development of insulin resistance with a glucose tolerance test (GTT). After a defined fasting period (5 h; from 8 AM to 1 PM) the animals received a single intraperitoneal injection of glucose (2 g/kg body wt). Blood samples were obtained at 0, 20, 40, 60, 120, and 180 min. The remaining 40 mice per group were randomized into five groups (n = 8) to study the dynamic changes within organs during a GTT. To this end, the animals followed the same GTT protocol as the monitor mice and were killed at five defined time points, namely, 0, 0.6, 2, 18, and 48 h after glucose injection. The data from the 0 h time point only are reported here. Mice were euthanized with CO/CO2 according to the local ethics recommendations, and organs were collected. During euthanasia, blood was collected by heart puncture to prepare serum (centrifugation of clotted blood at 6,000 rpm for 10 min).
In PPS3, mice were purchased from Charles River. At 12 wk of age all mice received a low-fat diet, in which 10% energy was obtained from fat, as a run-in for 4 wk. After the run-in period, mice were randomly divided into four weight-matched groups, which received four different diets with four increasing (10%, 20%, 30%, 45%) amounts of fat (n = 12/group). The diets used in PPS2 and PPS3 were identical and from one batch provided by Research Diets. Plasma samples were taken by orbital puncture every 2 wk. After 12 wk mice were fasted for 5 h (from 8 AM to 1 PM) and subsequently anesthetized with a mixture of isoflurane (1.5%), nitrous oxide (70%), and oxygen (30%) according to the local ethics recommendations. Blood was collected by orbital puncture, after which the mice were killed by cervical dislocation. Only the data from the high-fat and low-fat groups after 12-wk intervention are reported here. In both studies, for RNA isolation and proteome and metabolome analysis the liver was excised, divided into five lobes, and snap frozen in liquid nitrogen.
Clinical chemistry parameters in blood.
In PPS2, blood glucose levels were measured by the glucose hexokinase method (Instruchemie, Delfzij, The Netherlands). Insulin was measured with a Mercodia ELISA (Mercodia, Uppsala, Sweden). Cholesterol and triglycerides were measured as previously described (22). Serum folate was measured with a microbiological assay as previously described (31). Serum total homocysteine was measured by fluorescence polarization with an Abbott IMX instrument (Abbott) (27). In PPS3, blood glucose was measured after death with an Accu-Chek glucose meter (Roche Diagnostics, Almere, The Netherlands).
Amino acid analysis.
The amino acids in liver and plasma were analyzed with the iTRAQ methodology from Applied Biosystems. For liver amino acids, frozen liver samples were homogenized with a mortar and pestle. One hundred milligrams of homogenized tissue was dissolved in 150 μl of H2O-MeOH (50:50), vortexed, and centrifuged. Forty microliters of the supernatant or forty microliters of plasma was used to label the amino acids with the AA45/32 Starter Kit according to the manufacturer's instructions (Applied Biosystems, Foster City, CA) and analyzed via liquid chromatography-tandem mass spectrometry (LC-MS/MS) (3200QTRAP LC/MS/MS, Applied Biosystems). The data were analyzed with Analyst 1.5 software. Liver amino acid concentrations were normalized with the protein concentration, which was determined with the Bradford assay.
Glutathione assay in liver tissue.
A glutathione detection kit (Abcam, Cambridge, UK) was used to determine reduced (GSH), oxidized (GSSG) and total glutathione levels in liver tissue samples. Forty milligrams of each liver tissue sample was homogenized with 100 μl of ice-cold glutathione assay buffer. Sixty microliters of the homogenate was transferred to a prechilled tube containing perchloric acid and vortex mixed for several seconds to achieve a uniform emulsion. Samples were kept on ice for 5 min and centrifuged for 2 min at 13,500 rpm at 4°C. The supernatant (containing glutathione) was collected and stored at −80°C. The amount of protein in the pellet was determined with a BCA Protein Assay Kit (Thermo Scientific Pierce, Cramlington, UK). GSH and GSSG were measured according to the manufacturer's instructions and based on the reaction of o-phthalaldehyde with GSH.
Triglycerides and total phospholipids in liver.
Frozen liver samples were homogenized with a mortar and pestle, and triglycerides were determined after enzymatic hydrolysis with lipases and the Triglycerides liquicolormono kit (HUMAN, Wiesbaden, Germany). Total phospholipids were assayed with the commercial LabAssay Phospholipid kit (Wako Chemicals, Neuss, Germany). In brief, liver tissue was homogenized with liquid nitrogen and then dissolved in 300 μl of 0.9% NaCl. The samples were constantly shaken for 10 min and afterwards centrifuged at 4°C for 30 min. Twenty microliters of the obtained supernatant were used to run the assay. Triglycerides and total phospholipids were normalized with protein concentration, which was determined by Bradford assay.
Selected phospholipids and acylcarnitines in liver tissue.
Lipids were extracted with 60 μl of methanol per 10 mg of homogenized liver tissue. The samples were vortexed and centrifuged, and 10 μl of the supernatant was used for the assay. Acylcarnitines, phosphatidylcholines, and sphingolipids were analyzed with the AbsoluteIDQ kit (Biocrates Life Sciences) according to the manufacturer's instructions and with the 3200 QTrap-MS using MRM pairs. Samples were delivered to the mass spectrometer by the flow injection analysis method. Metabolite concentrations were normalized with the protein concentration, which was determined by Bradford assay.
Established protocols reported in detail recently (23) were used for RNA extraction, RNA quality control, microarray analysis, and microarray quality control were followed. Briefly, total RNA was extracted in PPS2 with RNAzol (Campro Scientific, Veenendaal, The Netherlands) and glass beads and in PPS3 with TRIzol reagent (Invitrogen, Carlsbad, CA), treated with DNase, and purified on columns with the RNeasy Mini Kit (Qiagen, Venlo, The Netherlands). The integrity of each RNA sample obtained was examined by Agilent Lab-on-a-chip technology using the RNA 6000 Nano LabChip kit and a Bioanalyzer 2100 (Agilent Technologies, Amstelveen, The Netherlands), and the quality of RNA preparations of PPS2 and PPS3 were comparable. RNA preparations were considered suitable for array hybridization only if samples showed intact 18S and 28S rRNA bands and displayed no chromosomal peaks or RNA degradation products (RNA integrity number > 8.0). Microarray analysis was carried out with an Affymetrix technology platform and NuGO Affymetrix Mouse GeneChip arrays (NuGO_Mm1a520177) containing 23,865 probe sets, including 73 control probe sets. The NuGO arrays are custom-designed Affymetrix GeneChip arrays designed by NuGO and manufactured by Affymetrix. These arrays contain in part common probe sets that are also present on standard Affymetrix arrays and in part newly designed probe sets (http://www.bigcat.unimaas.nl/∼martijn/NuGO/annotations/index.html). Briefly, fragmented cRNA was mixed with spiked controls and applied to Affymetrix Test chips, and good-quality samples were then used to hybridize with NuGO Affymetrix Mouse GeneChip arrays. The hybridization, probe array washing and staining, and washing procedures were executed as previously described (23), and probe arrays were scanned with a Hewlett-Packard Gene Array Scanner in PPS2 (Leiden Genome Technology Center, Leiden, The Netherlands) and with an Affymetrix GeneChip 3000 7G scanner in PPS3.
Quality control of microarray data was performed with Bioconductor packages (13), including simpleaffy and affyplm, through the NuGO pipeline that is available as a Genepattern procedure on http://nbx2.nugo.org (9). Probe sets were redefined according to Dai et al. (8) with remapped CDF version 11.0.2 based on the Entrez Gene database. The Nugo Mouse GeneChip arrays target 15,240 unique gene identifiers. GC-robust multiarray (GCRMA) analysis was used to obtain expression values (50). Only genes that had an unlogged intensity above 5 on at least five arrays of a group were considered for analysis. This criterion was met by 12,863 genes in PPS2 and 7,486 genes in PPS3. LIMMA was used to identify differentially expressed genes (40). Obtained P values were corrected for multiple testing with a false discovery rate (FDR) method (43).
The frozen liver samples were homogenized with a mortar and pestle. One hundred milligrams of homogenate was resuspended in 1 ml of lysis buffer [7 M urea, 2 M thiourea, 2% 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate (CHAPS), 1% DTT, Complete Mini proteinase inhibitor], and the protein concentration was determined by the Bradford method. Two-dimensional gel electrophoresis was performed as described previously (14), with slight modifications. In brief, isoelectric focusing was performed on 18-cm immobilized pH gradient (IPG) strips (pH 3–10, GE Healthcare, Uppsala, Sweden). IPG strips were passively rehydrated in 350 μl of solubilization buffer for 8 h, loading 350 μg of protein with rehydration buffer containing 7 M urea, 2 M thiourea, 4% CHAPS, 2% Pharmalyte 3–10, and an additional 15 μl of 30% DTT. The second dimension was performed with 12.5% SDS-polyacrylamide gels. The gels were run overnight at 20°C and constant amperage, fixed, and stained with Coomassie blue. Gels were analyzed with Decodon Delta 2D software 4.0 (Decodon, Greifswald, Germany). Expression of the different protein spots was expressed as volume percentage (%vol). Tryptic digestion of regulated protein spots and peptide mass fingerprinting by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) were performed as described previously (49).
Metabolite extraction from liver tissue for NMR and GC-MS analysis.
Aqueous and organic metabolites were extracted from liver tissue. Briefly, ∼100 mg of frozen tissue was homogenized in a 1:2 mixture of chloroform and methanol (600 μl). Samples were then sonicated three times on ice for 5 s. Six hundred microliters of chloroform-water (1:1) was added, and the mixture was vortexed thoroughly. Samples were left to stand for 15 min before centrifugation at 13,500 rpm for 20 min. The aqueous layer was dried under nitrogen gas and stored before NMR spectroscopy. The organic layer was combined with 20 μl of 2 mg/ml pentadecanoic acid (C15:0) as an internal standard and then dried under nitrogen gas. Dried organic extracts were derivatized by methylation; 1 ml of methanolic BF3 was added, and samples were incubated at 100°C for 45 min. Samples were then dried under nitrogen gas and stored before GC-MS analysis.
Dried liver tissue aqueous extracts were prepared for NMR spectroscopic analysis by the addition of 550 μl of D2O and 10 μl of trimethylsilyl propionate (TSP). Spectra were acquired on a 600-MHz Varian NMR spectrometer (Varian, Oxford, UK) using a Noesypresat pulse sequence at 25°C. Spectra were acquired with 32 K data points and 128 scans over a spectral width of 9.6 kHz. Water suppression was achieved during the relaxation delay (1.5 s) and the mixing time (100 ms). Serum samples were prepared by the addition of 10 μl of TSP and D2O to a total volume of 550 μl. Spectra were acquired with a Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence with 32 K data points and 64 scans over a spectral width of 9.6 kHz. Water suppression was achieved during the relaxation delay (3 s) and the mixing time (100 ms).
All 1H NMR spectra were referenced to TSP at 0.0 ppm and then phased and baseline corrected [initially with automated spectral processing software (Chenomx 6.0, Chenomx, Edmonton, Canada) and then manually where necessary]. A line-broadening factor of 0.2 Hz was also applied to all spectra. Chenomx 6.0 was also used for all spectra to calculate the integral of bins that consisted of spectral regions of 0.02 ppm. Bins from the water region of the spectra (4.5–5.5 ppm) were then excluded before the data were normalized to the sum of the spectral integral for each spectrum.
Dried derivatives were resuspended in 200 μl of hexane and analyzed on an Agilent 7890A GC coupled with a 5975C MS with an Agilent HP-5ms 30 m × 250 μm × 0.25 μm column (Agilent, Santa Clara, CA). One microliter of the sample was injected in splitless mode, and the initial oven temperature of 70°C was raised to 220°C at 5°C/min, held for 20 min, and then raised to 320°C at 20°C/min. Helium was used as the carrier gas with a flow of 1.2 ml/min. Calibration was achieved by comparison of peak areas for fatty acids with reference to a known standard (Supelco 37 compound mix, Supelco, Poole, UK) using Agilent Chemstation MSD E.02.00.493 and by comparison of their mass spectra with those in the NIST Library 2.0 (2005). Automatic peak detection was carried out with Agilent Chemstation mass spectral deconvolution (MSD). MSD was performed with AMDIS version 2.65. Peaks with a signal-to-noise ratio lower than 30 were rejected. To obtain accurate peak areas for the internal standard and specific peaks/compounds, one quant mass for each peak was specified as a target, and three masses were selected as qualifiers. Each data file was manually analyzed for false positives/negatives with Agilent Chemstation. The concentrations reported are percent fatty acids.
Data analysis and visualization.
Initial data analysis of NMR data was performed with principal components analysis (PCA) and partial least-squares discriminant analysis (PLS-DA). Once discriminating metabolites were identified, a targeted analysis was performed with Chenomx software. NMR data were imported into SIMCA-P (version 11, Umetrics, Umea, Sweden), mean centred, and Pareto scaled (1/√SD). Each data set (liver and serum) was analyzed with PCA, and the score plots were visually inspected for clustering trends and outliers. PLS-DA was then applied to the data sets, with diet as a classifier. Validation of metabolites and the GC-MS-derived fatty acid data and comparison between GSH and GSSG levels in liver tissue were performed in SPSS-V14 (SPSS) with one-way analysis of variance. Significance was assumed when P value < 0.05.
Data analyses of the remaining parameters were performed with the R software package (36). The different groups were compared by Student's t-test, and P values were adjusted for multiple testing with the p.adjust function and the Benjamini and Hochberg method (4). Significance was assumed when P value < 0.05. The packages Hmisc (17), plyr (47), and reshape (48) were used for data transformation. Transcriptome, proteome, and metabolome data were visualized with PathVisio (45).
Development of liver steatosis and insulin resistance under a high-fat feeding regimen.
After 12 wk of intervention with a high-fat diet, the animals showed a significant body weight increase (Fig. 1), became hyperglycemic, and showed a threefold increase in fasting insulin levels compared with control animals on a low-fat diet (Table 1). Additionally, plasma cholesterol increased significantly in the high-fat diet group. Total liver triglycerides were significantly increased in the animals fed a high-fat diet, whereas plasma triglyceride levels showed no statistically significant difference between the groups (Table 1). Food intake was similar on both diets, with a marginally higher intake for the low-fat diet [PPS3: no significant differences; PPS2: 3.30 ± 0.14 vs. 3.12 ± 0.12 g/day for low- and high-fat diets, respectively (P = 0.032)]. GTT performed after 12 wk in the monitor animals in PPS2 confirmed that the animals on the high-fat diet became glucose intolerant. After the high-fat diet mice had a significantly greater area under the glucose curve [3,729 ± 1,346 vs. 6,295 ± 1,170 arbitrary units (au); P < 0.01], indicating insulin resistance in these mice. Additionally, the HOMA score was significantly higher in animals after the high-fat diet (2.37 ± 1.61 vs. 0.55 ± 0.28; P < 0.005).
Metabolite and transcript levels identify one-carbon metabolism and related pathways as significantly changed in liver.
NMR, LC-MS/MS, and enzymatic analyses revealed that a number of liver metabolites were significantly different between the dietary groups: among these metabolites those related to one-carbon metabolism and phospholipid handling were found to be overrepresented (Fig. 2; an interactive version of Fig. 2 in which fold changes are visible is available at http://188.8.131.52/, username: pathways, password: !pathways!). Methionine, serine, glycine, ethanolamine, and phosphoethanolamine concentrations were significantly decreased, whereas choline, taurine, and sarcosine concentrations increased significantly in animals fed the high-fat diet. Liver glutathione levels were altered: the oxidized form (GSSG) was increased whereas the reduced form (GSH) decreased after treatment with a high-fat diet. Data from gene expression analysis revealed that in the same samples changes in mRNA levels coding for proteins in pathways related to one-carbon metabolism occurred (Fig. 2). High-fat diets led in both feeding trials to increased transcript levels of Shmt1, Shmt2, Dmgdh, Bhmt2, Chpt1, and Gpx7, whereas Cbs and Cdo1 transcript levels decreased. Additionally, the PPS2 data showed Dnmt3a, glutathione-related genes (Gclc and Gsr), and genes (Cept1 and Pemt1) involved in phospholipid metabolism with reduced expression states, whereas Chkb, Mat1a, and Baat mRNA increased. In PPS3 decreased levels of Mat2a and Csad transcripts were observed, whereas Dhfr of the folate cycle was increased. Proteome analysis revealed that 46 spots showed differential intensities between 12 wk of high-fat or control diet feeding. From these, 36 proteins were identified. However, none of these proteins was related to one-carbon metabolism.
Total phospholipids measured in the liver were not significantly altered (60.90 ± 29.01 and 68.19 ± 36.98 mg/g protein for low and high fat, respectively; P = 0.74). However, LC-MS/MS analysis of 21 acylcarnitines, 72 phosphatidylcholines, 9 lysophosphatidylcholines, and 15 sphingomyelins revealed that 24 of these 117 metabolites changed significantly after 12 wk of high-fat feeding. The changes occurred mainly in phosphatidylcholines containing long-chain polyunsaturated fatty acids with an ester or ether bond as well as three sphingomyelins and two lysophosphatidylcholines (Table 2). After the high-fat feeding the concentrations of a number of fatty acids were significantly changed in the liver (11 of 27 measured, Table 3). The total monounsaturated fatty acids increased, and n-3 polyunsaturated fatty acids including C20:5n3 and C22:6n3 decreased.
Related one-carbon metabolite changes in plasma.
Plasma amino acids related to one-carbon metabolism and related pathways were determined in PPS3. After 12 wk of high-fat feeding glycine, serine, cystine, phosphoethanolamine, methionine, and taurine showed similar levels in the animals fed either control or high-fat diet (data not shown). Ethanolamine was significantly increased in plasma of animals fed a high-fat diet (12.41 ± 1.23 and 14.32 ± 1.23 μmol/l for low- and high-fat diets, respectively; P = 0.03).
In PPS2, serum homocysteine concentrations were not significantly different between the groups (7.4 ± 1.2 and 6.6 ± 0.7 μmol/l for low and high fat, respectively). Serum folate concentrations did not change significantly after the high-fat diet (41.4 ± 4.2 and 43.8 ± 5.0 ng/ml in low- and high-fat groups, respectively). However, serum VLDL concentrations were significantly increased in animals on the high-fat diet (30.4 ± 2.7 and 39.0 ± 4.3 au in low- and high-fat groups, respectively; P < 0.001).
Combined analysis of the omics data from both PPS2 and PPS3 highlighted alterations in one-carbon metabolism and related pathways as a common effect of an increased dietary exposure to palm oil fat. Here we present a detailed overview of the findings obtained in these studies and furthermore demonstrate convincingly that agreement in results obtained in two different laboratories is possible even if previous reports have questioned the ability to achieve this (26).
C57BL/6J mice after 12 wk on a high-fat palm oil-based diet became obese, hyperglycemic, and insulin resistant and displayed elevated levels of plasma cholesterol and serum VLDL. Metabolite profiling of the liver samples revealed that numerous amino acids and derivatives related to one-carbon metabolism were significantly changed, highlighting this pathway as one of the main perturbations observed after the high-fat diet. Sarcosine, which can be considered the end product of methyl group turnover, and taurine for sulfur metabolism increased after high-fat feeding, whereas methionine, the precursor for these two products, decreased in concentration. Concentrations of glycine and serine as intermediates of methionine metabolism and the transsulfuration pathway decreased, while choline concentrations increased. Overall, these changes could indicate a compensation for low methionine levels with an increased remethylation from homocysteine to meet the demands for this essential amino acid with the needed methyl groups derived from choline. Since both diets contained equal amounts (2% by weight) of choline and mice on the low-fat diet consumed slightly more food, the higher choline level in livers of animals on the high-fat diet must therefore be of metabolic origin. No alterations in proteins involved in one-carbon metabolism were observed by the two-dimensional gel-based proteomics approach with Coomassie staining. Most probably our staining method was not sensitive enough to detect slight changes (46).
Choline may be derived from the degradation of phosphatidylcholines via betaine to deliver the methyl groups for remethylation of homocysteine with the production of dimethyl glycine and in turn sarcosine. Sarcosine levels increased on the high-fat diet, and transcript levels of two of the enzymes involved in the conversion from choline, including betaine-homocysteine S-methyltransferase (BHMT) showed higher mRNA levels. Interestingly, high levels of BHMT have previously been demonstrated to promote hepatic apolipoprotein B and VLDL production (41, 42), and we also observed increased serum VLDL levels in the obese mice. Lipid transport in general and hepatic VLDL secretion need phospholipids (28), with phosphatidylcholines being the most important. Phosphatidylcholines are derived in liver mainly from the Kennedy pathway. Although transcripts of the rate-limiting enzymes phosphate cytidylyltransferase 1 choline α and β (Pcyt1a-b) were not determined, choline phosphotransferase 1, which catalyzes the last step in the synthesis of phosphatidylcholines, showed increased levels. Hepatic levels of ethanolamine and phosphoethanolamine, which are precursors for phosphatidylethanolamine synthesis, were decreased, and mRNA levels of choline/ethanolaminephosphotransferase 1 (Cept1) were also decreased. In addition, a number of the phosphatidylcholines showed decreased levels, suggesting that the hepatic ratios of the main phospholipids are altered (Supplemental Fig. S1). Although total liver phospholipid status did not reveal differences, our data strongly argue for an enhanced utilization of choline derived from phosphatidylcholine breakdown to accommodate the increased demand for remethylation of homocysteine to obtain methionine.
The crucial role of choline in hepatic fat metabolism is supported by other studies in rodents and humans. Choline deficiency alone did not lead to impaired hepatic lipid export, whereas methionine and choline deficiency in combination did (25). A recent paper showed that Pemt-null mice are resistant to diet-induced obesity and that this phenotype is reversed by administration of choline (21). Moreover, two recent human studies indicate that the susceptibility to develop a fatty liver when consuming a choline-deficient diet is associated with mutations in the Pemt promoter (7, 39). In the Hordaland Health study the relationship between plasma choline levels and cardiovascular risk factors such as obesity and glucose metabolism was assessed. Plasma choline levels were found to be positively associated with body mass index and perturbed triglyceride and glucose levels and inversely related to HDL-c (24). In a mouse study of insulin resistance, the interplay of the host-gut microflora metabolism with respect to choline metabolism was also put forward as a potential contributor to the development of insulin resistance (11).
Activity of the transsulfuration pathway is high in the liver, and the rate-limiting step is catalyzed by cystathionine β-synthase (CBS). Although a number of transcripts in this pathway were significantly downregulated, one of the metabolic end products (taurine) was increased. Most interestingly, the Cbs-null mouse phenotype is characterized by increased plasma triglyceride and VLDL levels and hyperhomocysteinemia (32) and thereby also demonstrates the close link between sulfur and methyl group metabolism and impaired lipid handling. Taurine can be used for bile acid synthesis, and since high fat intake increases the demand of bile acids for proper lipid solubilization for fat absorption in the intestine, the increased taurine may be needed in the liver for bile acid synthesis. The bile acid coenzyme A:amino acid N-acyltransferase (Baat) transcript was indeed increased, and this may suggest that more taurine is used for conjugation in taurocholate biosynthesis. The same may hold true for glycine needed for formation of glycocholate. In addition, bile acids in liver are known to contribute to the regulation of hepatic lipid metabolism and insulin sensitivity (44). It is possible to hypothesize that the increased taurine levels and the suggested increase in the remethylation of homocysteine to methionine may simultaneously reduce the availability of cysteine for de novo glutathione (GSH) synthesis. Indeed, we observed a significant change in the ratio of GSH and GSSG after high-fat feeding that may be related to a higher reactive oxygen species (ROS) burden and/or a reduced capacity for reduction of GSSG. In addition, protein expression level of catalase was increased in liver tissue of the animals on the high-fat diet, providing a second line of evidence for alterations in hepatic ROS handling. The presence of oxidative stress in liver after high-fat feeding has been demonstrated before (37) and is proposed to represent an important trigger for the development of insulin resistance (30). Although low glutathione levels per se could not be linked directly to hepatic triglyceride accumulation and the development of diet-induced steatohepatitis (5, 16), it is still possible that the altered antioxidant capacity observed here is a contributing factor to the development of insulin resistance.
In addition to alterations in the phosphocholine spectrum in liver, we also observed changes in other lipids supporting the hypothesis that an altered lipid milieu contributes to the development of insulin resistance (22). Included in these changes were altered levels of carnitines, sphingolipids, and fatty acids. Carnitine C5-M-DC was decreased in the high-fat diet group. Although specific functions of the different carnitines are not known, methylglutaryl carnitine has been previously related to insulin resistance and C5-M-DC levels in plasma of diabetic mice are increased (1). Sphingolipids that constitute a minor component of membrane lipids have also been found to change in insulin resistance and T2DM (18, 20). In the present study three sphingomyelin species decreased in liver and transcripts of genes involved in sphingolipid metabolism were significantly reduced as well in the obese animals. Certain saturated fatty acids and monounsaturated fatty acids were increased in the livers of animals after the high-fat diet. The changes in monounsaturated fats possibly reflect the higher dietary intake of palm oil in the high-fat diet. Interestingly, a number of n-3 polyunsaturated fatty acids and total polyunsaturated fatty acids decreased significantly, and these fatty acids are known to have a beneficial effect on insulin resistance (12, 52). Overall, these results provide strong evidence for a significant change in the hepatic lipid milieu that is likely contributing to the pathogenesis of insulin resistance, in agreement with recent results (22). Moreover, lipidomic analysis of liver tissue from patients with nonalcoholic fatty liver disease showed alterations in lipid homeostasis and suggests that these changes may play a role in the pathogenesis of the disease (34, 35).
In summary, we provide evidence for major changes in hepatic one-carbon metabolism based on metabolite, proteome, and transcriptome profiling that are associated with fat accumulation in liver following 12 wk on a high-fat diet. Although the samples were generated in two independent feeding trials conducted at two centers with C57BL/6J mice on identical diets, the observed changes in liver metabolism were robust and very similar in characterizing a condition of insulin resistance that promotes hepatic steatosis. Changes in the metabolism of sulfur-containing amino acids and methyl groups seem to dominate the metabolic perturbations. The increased demand for phospholipids and bile acids for lipid absorption in the gut and for interorgan fat transport during high-fat feeding may initiate the changes in one-carbon metabolism in liver. The link between the Kennedy pathway and one-carbon metabolism is provided by hepatic choline concentrations that increased and hepatic methionine that decreased in liver. Secondary to the changes in one-carbon metabolism may be the alterations in ROS status and changes in the lipid milieu. These main findings in combination with the known alterations in lipid metabolism in human subjects should form the basis of future hypothesis-driven studies examining the role of these pathways in the development of insulin resistance and hepatic steatosis.
This study was funded by the European Nutrigenomics Organisation, an EC-funded Network of Excellence (Grant No. FOOD-2004-506360).
No conflicts of interest, financial or otherwise, are declared by the author(s).
We thank all the NuGO PPS Members for support and fruitful discussions. We acknowledge Anne Molloy for measurement of folate and homocysteine. The Biocrates Absolute IDQ Kit was kindly provided by Biocrates Life Science AG (Innsbruck, Austria).
↵1 Supplemental Material for this article is available online at the Journal website.
- Copyright © 2011 the American Physiological Society