Introduction

Obesity is caused by an imbalance between energy intake and energy utilization. The excess of energy is accumulated in lipid droplets and prolonged caloric overconsumption results in a positive energy balance and increasing body weight. Several animal models were introduced to study behavioral and metabolic consequences of overnutrition (Larter and Yeh 2008). Of them, both prolonged feeding with a high-fat diet (HFD) and disruption of leptin signaling in leptin-deficient ob/ob and leptin receptor-deficient db/db hyperphagic mice were shown to result in obesity, insulin resistance, oxidative stress, changes in lipid and glucose homeostasis, and hepatic steatosis (Kennedy et al. 2010; Panchal and Brown 2011). However, since dietary behaviors significantly differ between mice prolonged fed with HFD and hyperphagic mice fed with normal diet (ND), they may exhibit differences also in molecular response to the deposition of a body fat.

Several microarray-based studies of obese humans and animal obesity models have been conducted during the past decade. Most of these studies have focused on identifying the transcriptional changes underlying the pathophysiology of obesity, as well as the anti-obesity effects of dietary, pharmacological, and behavioral interventions (Sun 2007). In a consequence, thousands of genes differentially expressed between lean and obese individuals have been identified, including those related to the hepatic glucose and lipid metabolism, insulin signaling, inflammation, coagulation, cell adhesion, oxygen stress, and activity of chaperone (Sharma et al. 2006; Greco et al. 2008; Kim 2010; Naik et al. 2013). On the contrary, a limited number of muscle transcriptomic studies have been reported so far. In one series of them, a short-term feeding of mice with HFD changed the expression of more than 1,000 genes in the skeletal muscle (de Wilde et al. 2008), but an 8-week HFD feeding induced rather small alterations in the muscle transcriptome, with up-regulation of only Fatty acid oxidation pathway (de Wilde et al. 2009). Others studies revealed, however, that overall gene expression in muscle, involved in fatty acid oxidation and biosynthesis, was either decreased or obesity unresponsive (Kim 2010).

In contrast to microarray-based studies, only a few reports on the proteomic profiling in obesity can be found in the literature. In previous SDS–PAGE-based (Luo et al. 2012), 2D-based (Park et al. 2011; Eccleston et al. 2011), and 2D differential in-gel electrophoresis-based (Douette et al. 2005; Zhang et al. 2010; Kirpich et al. 2011; Thomas et al. 2012) analyses of total liver proteins from rodents chronically exposed to HFD, between 12 and 53 proteins reportedly displayed significant obesity-related changes in expression. After short-term feeding with HFD, one hepatic protein was significantly upregulated (Thomas et al. 2012). Recent quantitative proteomic analysis of total liver proteome identified 1,131 proteins and expression of 90 of them differed significantly between obese HFD-fed rats and lean rats (Baiges et al. 2010).

Because the ability of mass spectrometry (MS)-based approaches to identify proteins expressed at low levels may be enhanced when the analyzed peptides are restricted to selected sub-proteomes (Adachi et al. 2007; Altelaar and Heck 2012), we focused on mitochondrial (MT) proteome alterations. The MT sub-proteome is suitable for the detection of molecular alterations associated with obesity for several reasons: (1) obesity is a consequence of energy imbalance; (2) cellular energy is largely supplied by the MT oxidative phosphorylation system, which is regulated predominantly by insulin in muscle, liver, and adipose cells; (3) transcription rates from both MT and nuclear genomes are activated by insulin (Boirie 2003; Mikula et al. 2005), and genes that belong to MT metabolic, catabolic, and biosynthetic pathways are highly expressed in obese mice (Sharma et al. 2010); (4) MT dysfunctions have been implicated in the pathogenesis of obesity and obesity-associated diseases, including metabolic syndrome, non-alcoholic fatty liver disease, type 2 diabetes, and atherosclerotic heart disease (Samuel and Shulman 2012); and finally, (5) the mammalian MT proteome consists of slightly more than 1,000 proteins (Pagliarini et al. 2008).

Transcriptomic and proteomic results of the studies performed in obesity models have not been integrated so far. Furthermore, direct comparisons of liver and muscle MT proteomes in obesity models have not been performed to date. Here, MT proteomes were quantitatively evaluated using isobaric tags for relative and absolute quantitation (iTRAQ), followed by liquid chromatography–tandem MS (LC–MS/MS) analysis, and transcript measurements were generated using microarray-based studies. Our studies reveal that combining MT proteomic and transcriptomic surveys enriches the functional analysis of obesity and fasting states in mice. Combined data on differentially expressed genes and proteins suggest that both obesity and fasting states result from highly disordered metabolic processes which in prolonged obesity is followed by deeper cellular and organismal systems deficiencies. In contrast, numbers of proteins and transcripts differentially expressed in muscle samples were too low to be assigned to any of functional pathways. Thus, altered expression of proteins and corresponding genes accompanied obesity and fasting states is quantitatively different in the liver and skeletal muscle.

Materials and methods

Animals

Ninety-six 5-week-old male, certificated mice were purchased from The Jackson Laboratory, USA. Of them, 48 mice were wild-type (wt) C57BL/6J, 24 mice were homozygous for the leptin gene mutation B6.V-Lep ob/J (ob/ob), and 24 mice were homozygous for the leptin receptor gene mutation B6.BKS(D)-Lepr db/J (db/db). Throughout the study, mice were housed under temperature-controlled (21 ± 2 °C) and humidity-controlled (55 ± 10 %) conditions, with a 12-h light/dark cycle and ad libitum access to food and water.

Ethics statement

Mice were housed in the Department of Genetics and Laboratory Animals Breeding, at the Maria Sklodowska-Curie Memorial Cancer Center and Institute of Oncology. Experimental protocols were approved by The 2nd Local Ethical Committee for Animal Research in Warsaw, Poland.

Experimental design

Following purchase, mice were allowed 1-week adaptation, during which all animals were fed ND (10 % of calories from fat), containing 19.2 % protein, 67.3 % carbohydrate, and 4.3 % fat (D12450B; Research Diets, New Brunswick, NJ, USA). At 6 weeks of age, 24 wt C57BL/6J mice were fed ND (control group) and other 24 wt C57BL/6J mice were fed HFD (HFD group) (60 % of calories from fat), containing 26.2 % protein, 26.3 % carbohydrate, and 34.9 % fat (D12492; Research Diets). The ob/ob and db/db mice were fed ND throughout the entire experiment. To determine the effect of fasting, half of the mice in each group were deprived of food for a period of 18 h (between 3 pm and 9 am next day) before sacrifice. An overview of the experimental design is shown in Fig. 1. At either 16 weeks or 48 weeks of age, mice were weighed and sacrificed, followed by rapid collection of blood, livers, and skeletal muscles. Fragments of fresh tissue samples used for mitochondria isolation were placed into ice-cold phosphate buffered saline (PBS). The remaining samples were snap-frozen and stored at −72 °C until use.

Fig. 1
figure 1

Overview of the experimental design. At the experiment endpoints, total RNA and mitochondrial protein extracts were combined from two animals (depicted as differential mice outlines) from each group of six, generating three biological replicate samples of pooled RNA and mitochondrial proteins per trait

mRNA extraction

Total RNA was isolated from liver and muscle samples using the RNeasy Plus Mini Kit (Qiagen, Germany) or TRIzol (Invitrogen Co., CA, USA), respectively, followed by on-column DNAse I digestion. The quality of RNA samples was determined using an Agilent 2100 Bioanalyzer (Agilent Technologies, CA, USA); samples used for microarray analysis displayed distinct peaks corresponding to intact 28S and 18S ribosomal RNA. Equal amounts of RNA (500 ng) were combined from two randomly selected animals. This process was performed for the six animals in each group, resulting in three biological replicates of pooled RNA per trait (Table S1 depicts the sample pooling process).

Microarray analyses

The average signal from MouseRef-8 v2.0 Expression BeadChips (AROS Applied Biotechnology, Aarhus, Denmark) was quantile normalized with no background correction. All computations were performed using R 2.15.0 software (R Development Core Team 2011) with the Bioconductor extension (Gentleman et al. 2004). Illumina and Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa and Goto 2000) identifiers were mapped to genes using lumiMouseAll.db (version 1.18.0), KEGG.db (version 2.8.0), and lumi (version 2.8.0) (Du et al. 2008) packages. Measurements were filtered according to the ratio of the interquartile range to the median. Only probes with a ratio greater than 1:8 were selected for analysis.

In the further analysis, only the probe sets corresponding to genes with detected proteins were considered in this study. Differentially expressed genes were identified by a t test (Welch variant) followed by the Benjamini–Hochberg p value correction for multiple hypotheses testing. Adjusted p values ≤ 0.05 were considered significant. Microarray data was deposited in Gene Expression Omnibus dataset GSE43691.

Functional analyses

The significance of mRNA and protein expression changes related to each KEGG pathway was assessed with the Kolmogorov–Smirnov (K-S) test. Resulting p values were adjusted for multiple hypotheses testing with the Benjamini–Hochberg algorithm and values ≤0.05 were considered significant.

To explore the network of proteins and their contribution to KEGG pathways, we used the STRING software (Szklarczyk et al. 2011) that constructs functional protein-interaction networks and calculates its significance based on aggregated experimental data.

Isolation of mitochondria

Mitochondria from 100 mg of mouse liver or hind limb muscle tissues were isolated and purified using differential centrifugation followed by ultracentrifugation in a discontinuous sucrose gradient (Ostrowski et al. 2004). The tissues were minced into small pieces using scissors and washed with ice-cold PBS. Minced muscles were incubated for 1 h in ice-cold PBS supplemented with 10 mM EDTA and 0.05 % trypsin to facilitate the homogenization process. Tissues were homogenized using a Mitochondria/Cytosol Fractionation Kit (BioVision, Milpitas, CA, USA). In brief, samples were homogenized in Cytosol Extraction Buffer supplemented with protease inhibitor cocktail, DTT, and phosphatase inhibitor cocktail 1 and 2 (Sigma). The homogenates were centrifuged (once for muscle, twice for liver) for 10 min at 740×g and 4 °C, and then the unbroken tissue, cells, and nuclei were discarded. The supernatants were then centrifuged for 30 min at 10,000×g and 4 °C to pellet the mitochondria. MT pellets were suspended in ice-cold 0.25 M sucrose, layered on 1 M and 1.7 M discontinuous sucrose gradients spiked with EDTA-free Complete Protease Inhibitor Cocktail (Roche, Basel, Switzerland), and then centrifuged for 45 min at 70,000×g and 4 °C. Mitochondria recovered from the gradient were diluted in ice-cold PBS and centrifuged for 20 min at 10,000×g and 4 °C. The purified mitochondria were stored at −80 °C.

Protein extraction

MT proteins were extracted with 8 M urea in 2 % CHAPS and concentrations were measured using a standard BCA assay. Equal amounts of MT protein (100 μg) were combined from two randomly selected animals; this process was performed for the six animals in each group, resulting in three biological replicates, similarly as for microarray analyses (Table S1). Pooled samples were precipitated using the ProteoExtract Protein Precipitation Kit (Calbiochem, Germany), according to the manufacturer’s protocol. Protein pellets were suspended in 40 μl dissolution buffer (0.5 M triethylammonium bicarbonate, 0.1 % SDS). To facilitate protein solubilization, samples were disrupted using a Bioruptor ultrasonic homogenizer (Diagenode, Belgium), using 10-s on–off cycles for 5 min at high intensity, and were then digested with trypsin at 37 °C overnight.

iTRAQ labeling

The concentration of tryptic peptides was measured using a modified BCA method (Kapoor et al. 2009). Equal amounts of peptides (100 μg) were then reduced and cysteine-blocked according to the iTRAQ protocol. Peptides derived from a corresponding set of samples, comprising eight experimental traits per age group, were labeled using the iTRAQ 8plex kit (ABsciex, Framingham, MA, USA), as shown in Table S1. Peptides were then mixed and fractionated simultaneously, using isoelectric focusing (IEF) on 24 cm IPG strips (pH 3.5 to 4.5; GE Healthcare), as described previously (Mikula et al. 2011). Strips were cut into 14 sections and peptides were extracted and stored at −72 °C. For each age group, three strips were analyzed; each strip represented three biological replicates (each replicate composed of MT proteomes derived from two different, randomly selected animals).

LC–MS/MS

Analysis was performed using an Orbitrap MS (Thermo Scientific, USA) coupled with a nanoAcquity LC system (Waters, USA). Spectrometer parameters were as follows: capillary voltage, 2.5 kV; cone, 40 V; N2 gas flow, 0; and m/z range, 300–2,000. Samples were loaded from the autosampler tray (cooled to 10 °C) to the precolumn [Symmetry C18, 180 μm × 20 mm, 5 μm (Waters)] using a mobile phase of 100 % MilliQ water acidified with 0.1 % formic acid. Peptides were transferred to the nano-UPLC column [BEH130 C18, 75 μm × 250 mm, 1.7 mm (Waters)] by a gradient of 5–30 % acetonitrile in 0.1 % formic acid over 180 min at a flow rate of 250 nl/min and then directly eluted to the ion source of MS. A blank run was performed before each analysis to ensure the lack of cross-contamination from previous samples.

LC–MS/MS data processing and analysis

All MS/MS data were analyzed using Mascot (version 2.3.02; Matrix Science), which was configured to search the mouse Swiss-Prot protein sequence database and reversed versions of these sequences assuming semiTrypsin digestion. The database search was performed with a fragment ion mass tolerance of 0.8 Da and a peptide mass tolerance of 40 ppm. ITRAQ8PLEX (N-TERM), ITRAQ8PLEX (K), and METHYLTHIO (C) were specified in Mascot as fixed modifications; ITRAQ8PLEX (Y) and OXIDATION (M) were specified as variable modifications. The statistical significance of peptide identifications was assessed by a target/decoy database search strategy, as described previously (Mikula et al. 2011). The quantitative iTRAQ analysis was performed using a set of iTRAQ-labeled unique peptides, which were unambiguously assigned to a single protein, with an estimated false discovery ratio (FDR) of 0.01. The intensities of the iTRAQ reporter ion peaks were corrected for isotope impurity and then median-normalized to remove systematic bias. For the statistical analysis of protein ratios, a permutation-based test implemented in the Diffprot software was used (Malinowska et al. 2012). Protein ratios with corrected p values ≤ 0.05 in pair-wise comparisons were considered significant. To assign proteins intracellular localization, both Swiss-Prot and Gene Ontology (GO) (Ashburner et al. 2000) were used.

Results

In order to define molecular alteration associated with obesity and fasting states in young and old obese mice, we used hepatic and muscle samples from mice at 16 and 48 weeks of age. Ten weeks of HFD feeding (16 weeks of mice age) was chosen since at this time HFD-fed mice already show obesity phenotypes, while some of 48-week-old obese mice started to show enhanced symptoms of severe disease, including failure to drink and eat, slower respiration, and impaired responsiveness to touch. A total of 12 mice per group (given genotype, diet, and age) were analyzed. Six mice in each group were deprived of food for a period of 18 h before sacrifice, according to typical mice metabolic studies, when an overnight fasting (for 14–18 h) has been usually employed (Ayala et al. 2010). Three biological replicates per group were generated by combining an equal amount of total RNA or MT protein extracts from two randomly selected animals (Fig. 1).

As expected, the body weights of both 16-week-old and 48-week-old ob/ob, db/db, and HFD-fed wt mice (determined immediately prior to sacrifice) were significantly higher than the weights of age-matched ND-fed control mice (Fig. 2a). Moreover, ob/ob and db/db mice at 16 weeks of age developed hyperinsulinemia, hyperglycemia, and moderate to severe hepatic steatosis, and 48-week-old ob/ob, db/db as well as wt C57BL/6J mice (in parallel with prolonged HFD feeding) showed severe steatosis with moderate inflammation (our yet unpublished observations). Compared with the non-fasting condition, fasting for 18 h resulted in significantly reduced liver weight in young ob/ob and HFD-fed mice, as well as in both young and old control mice (Fig. 2b). Additionally, overnight fasting resulted in significantly reduced serum glucose concentration in young and old control mice, as well as in young HFD-fed mice, and reduced serum insulin concentration in all studied mice (not shown). Although both the glucose and insulin concentrations were significantly affected by the obese state, we could not diagnose diabetes in any of the obese groups (not shown) according to its accepted diagnostic criteria for mice (Andrikopoulos et al. 2008).

Fig. 2
figure 2

Mouse total body weight (a) and liver weight (b) according to genotype and diet. HFD fed a high-fat diet, NF non-fasted, F fasted. *p < 0.01, **p < 0.005, compared with the relevant control group; #p < 0.01 for liver weight comparison between non-fasted and fasted mice in a given group

Comparison of liver and muscle proteomes

Isolated mitochondria were purified and the extracted peptides were labeled with iTRAQ, fractionated with IEF, and analyzed by LC–MS/MS. A total of 168 analytical LC–MS/MS runs were performed, producing 2,566,757 MS/MS spectra. A mouse Swiss-Prot protein sequence database search identified 17,606 and 10,967 peptides (estimated FDR ≤ 0.01) that were assigned to 2,746/1,675 and 1,484/704 proteins (identified by at least one/two peptides) in liver and muscle samples, respectively, which represented both MT and co-purified non-MT proteomes (Table S2). Of these proteins, 904/490 were present in data sets from both tissue types. MT proteins were then selected using information deposited in the Swiss-Prot and GO databases. According to the Swiss-Prot localization information, 513/406 hepatic and 504/396 muscle MT proteins were identified by one/two peptides. According to the GO term “mitochondria” (GO:0005739), 729/565 hepatic MT proteins and 504/396 muscle MT proteins were identified by one/two peptides. For clarity of further quantitative analyses, we considered only proteins identified by at least two peptides.

Proteins with an MT origin status in the Swiss-Prot database represented 22 % and 53 % of all proteins isolated from liver and muscle MT samples, respectively (Fig. 3). By comparison, MT proteins assigned by GO annotation represented 34 % of all hepatic proteins and 56 % of all muscle proteins identified. Notably, the hnRNP K protein, which was shown to be localized in mitochondria (Mikula et al. 2006), could not be attributed as a MT protein according to either the Swiss-Prot or GO data sets. Additionally, a recently published study on MT proteome mapping in living cells identified over 30 proteins within the human MT matrix that were not previously linked to mitochondria (Rhee et al. 2013). Since MT proteome is still not well established, for further analyses we used data for all proteins isolated with MT fractions.

Fig. 3
figure 3

Subcellular identity of the 1,675 and 704 proteins identified by at least two peptides in mitochondrial fractions isolated from liver (a) and muscle (b), respectively. Proteins were identified using the Swiss-Prot database

To compare the protein abundances with the corresponding mRNA transcript levels, we performed microarray measurements of pooled RNA samples combined from same two animals as MT protein pools (Table S1). In total, 48 microarrays were hybridized with liver samples and 48 additional microarrays were hybridized with muscle samples.

Molecular alterations underlying obesity induced by hyperphagia and HFD

Obesity-related differentially expressed proteins (both MT and co-purified non-MT) were analyzed in pair-wise comparisons of non-fasted mice; the analysis was independent for each group of obese mice. The analysis of liver samples identified 123, 115, and 67 proteins that distinguished (adjusted p values ≤ 0.05) young control mice from young ob/ob, db/db, and HFD-fed mice, respectively, and also 74, 79, and 41 proteins that distinguished old control mice from old ob/ob, db/db, and HFD-fed mice, respectively (Table 1). Analysis of muscle samples revealed 19, 22, and 7 proteins that distinguished young control mice from young ob/ob, db/db, and HFD-fed mice, respectively, and 5, 12, and 13 proteins that distinguished old control mice from old ob/ob, db/db, and HFD-fed mice, respectively (Table 1). The Venn diagrams on Fig. 4 present the number of differentially expressed proteins, common or unique between models of obesity. In sum, 221 proteins were significantly altered (adjusted p values ≤ 0.05) in liver (Table S3), and 44 proteins were altered in muscle (Table S4).

Table 1 Numbers of detected proteins and numbers of proteins that corresponded to measurable mRNA used for comparisons between obese mice and appropriate control mice
Fig. 4
figure 4

Venn diagrams presenting common and unique differentially expressed proteins in different mice models of obesity

Interestingly, a set of 27 hepatic proteins exhibited opposite direction of expression changes in HFD-fed and hyperphagic mice when compared to control (Table 2). Of them, 22 proteins appeared as differential in 16-week-old animals. Molecular context of these proteins interaction was next explored with STRING database. Analysis revealed that these proteins are highly interconnected (protein–protein interactions p value 5.66E−15) (Fig. 5) and significantly enriched for KEGG fatty acids-related pathways (Peroxisome, PPAR signaling, Fatty acid metabolism) and xenobiotics metabolism. The central proteins of network are cytochrome P450 3a11 (Cyp3a11), fructose-bisphosphate aldolase B (Aldob), and Cyp4a14 which bridged proteins cluster involved in Metabolism of xenobiotics with components of Fatty acid metabolism and PPAR signaling pathways. According to STRING, the best documented associations in the network are between Cyp3a11 and Cyp4a14 proteins. Eight of these 27 proteins showed over 3-fold difference in expression between hyperphagic and HFD-fed mice at any of age point; four proteins were upregulated [glutathione S-transferase Mu 1 (Gstm1), ATP-binding cassette sub-family D member 2 (Abcd2), acyl-CoA desaturase 1 (Scd1), and fatty acid synthase (Fasn)] and four proteins were downregulated [cytochrome P450 4V3 (Cyp4V3), liver carboxylesterase 31 (Es31), ornithine aminotransferase (Oat), and serine-pyruvate aminotransferase (Agxt)] in hyperphagic mice. There were no proteins detected in the muscles that would have opposite expression in HFD-fed and hyperphagic mice.

Table 2 Differentially regulated proteins between hyperphagic and HFD-fed mice
Fig. 5
figure 5

A protein interaction network constructed with STRING (Szklarczyk et al. 2011) for a set of 27 hepatic proteins with an opposite direction of expression in HFD-fed and hyperphagic mice when compared to control. Stronger evidences of associations are represented by thicker network edges. Node color depicts significantly enriched KEGG category assigned to a protein by STRING database in a set of 27 proteins while white nodes indicate proteins without functional annotation

Next, all identified proteins and their corresponding probe sets, from liver or muscle samples (Table S5 and S6, respectively), were analyzed separately by annotation to pre-defined KEGG pathways (Table S7). Liver proteins could be assigned to 25 KEGG pathways in at least one comparison, while their corresponding probe sets were assigned to 18 pathways; of these, the Fatty acid metabolism, Biosynthesis of unsaturated fatty acids, Glycolysis/Gluconeogenesis, PPAR signaling, and Peroxisome pathways were common for the annotations of both proteins and their corresponding probe sets. Contrary to liver samples, no functional processes were selected from either muscle RNA or muscle protein obesity data sets.

Molecular alterations underlying prolonged fasting

Overnight fasting-related differentially expressed proteins were also independently identified in liver and muscle samples from each group of animals. In sum, fasting altered the expression (adjusted p values ≤ 0.05) of 107 proteins in liver samples and 35 proteins in muscle samples (Tables S8 and S9). The higher number of altered proteins was shown in both young and old control mice, comparing with all obese mice (Table 3). Noteworthy, a set of 45 hepatic proteins demonstrated concordant abundance change in response to fasting, in at least two groups of animals (Table S8); 38 of them were upregulated. When analyzed against the STRING database, these 45 proteins constructed a network significantly enriched in interactions that contribute to KEGG pathways related to fatty acids (Peroxisome, PPAR signaling, Fatty acid metabolism) and amino acids metabolism (Fig. 6). Four and three of differentially expressed proteins were common for all groups of young or old mice, respectively (Fig. 7a, b), while two of them, encoding components of cytochrome P450 (Cyp4a14 and Cyp4a10), were changed concordantly in both age groups (Table S8).

Table 3 Numbers of detected proteins and numbers of proteins that corresponded to measurable mRNA used for comparisons between mice subjected to prolonged fasting or non-fasting in each group
Fig. 6
figure 6

STRING interaction network (Szklarczyk et al. 2011) for a set of 45 hepatic proteins with concordant expression during fasting in at least two groups of mice. Stronger evidences of associations are represented by thicker network edges. Node color depicts significantly enriched KEGG category assigned to a protein by STRING database in a set of 45 proteins while white nodes indicate proteins without functional annotation

Fig. 7
figure 7

Venn diagrams presenting common and unique differentially expressed proteins during overnight fasting in obese and control mice

In skeletal muscles, we found 10 consequently changed proteins in at least two groups of animals during fasting (Table S9). According to STRING, three of them indicated interactions [myosin light chain 1/3 (Myl1), tropomyosin alpha-1 chain (Tpm1), and calsequestrin-1 (Casq1)]; however, these 10 proteins were not enriched in KEGG pathways. In muscle samples from 16-week-old animals, only wt mice (regardless of the diet) showed commonly changed proteins (Fig. 7c).

Next, the functional analysis of transcriptomic and proteomic data sets, differentiating non-fasting from fasting mice, were performed (Table S10). Seventeen KEGG pathways were attributed to changes in mRNA levels associated with overnight fasting. One of these pathways (Peroxisome) was selected in all groups of animals, and three others (Fatty acid metabolism; PPAR signaling; and Valine, leucine, and isoleucine degradation) were selected in at least four out of eight animal groups. Differentially expressed hepatic proteins could be assigned to nine KEGG pathways, but only in three groups of old obese mice. Of these nine pathways, five pathways (Peroxisome, Fatty acid metabolism, and three neurodegenerative diseases) were also selected upon corresponding probe sets data (Table S10). Functional analysis of muscle samples revealed 21 KEGG pathways selected from proteomic data sets and seven selected from transcriptomic data sets (Table S11). Common protein and RNA-related functional annotations were limited to only two groups of mice: young HFD-fed mice and old control mice.

Discussion

In this study, for the first time, we demonstrate that the extent of MT proteome changes is greater in the liver than muscles in states of well-established models of obesity. The difference was captured in the number of identified proteins—1,675 and 704 proteins represented by at least two peptides, for liver and muscles, respectively, as well as in the number of differentially expressed proteins in these two tissues (Figs. 4 and 7). In the liver, but not in skeletal muscles, we point functional changes, depicted by KEGG pathways, distinctive to obesity (Table S7) and fasting (Table S10) conditions by collating MT proteome dataset with corresponding transcripts. As expected, these pathways are mostly related to fatty acids metabolism pathways.

Difficulty defining the mitochondrial proteome

Depending on the method used to identify proteins, mitochondria are thought to contain between 1,000 and 1,500 nuclear-encoded proteins. However, no single method can identify all MT proteins; some are expressed at levels too low to detect, while others are expressed exclusively by specific cell types (Calvo and Mootha 2010). Furthermore, distinguishing genuine MT proteins using an MS approach may be obscured by the co-purification of protein contaminants from other membranous and MT-associated organelles (Pagliarini et al. 2008). MT proteins are directed toward the MT outer membrane either during or after translation, as their synthesis can be initiated and completed on cytoplasmic ribosomes or initiated in the cytoplasm and completed within mitochondria. Synthesis of some MT proteins takes place in the vicinity of the organelle where their corresponding mRNAs are bound to mitochondria before the initiation of translation (Pagliarini et al. 2008).

Skeletal muscle contains two types of mitochondria: intermyofibrillar mitochondria, which provide most of the energy used for muscle contraction; and subsarcolemmal mitochondria, which generate ATP for membrane functions (Ferreira et al. 2010). Most of the top ten skeletal muscle proteins identified by the greatest number of peptides in this study represented components associated with the assembly or function of vertebrate skeletal muscles, including several isoforms of myosin (Table S2). Therefore, a substantial proportion of proteins identified as non-MT might originate from organelles that are essential for MT protein translation and that co-localize topographically and/or functionally with mitochondria. Finally, some of proteins with actual MT localization might not be assigned to mitochondria in available databases (Rhee et al. 2013). With this in mind, the differentially expressed proteins and their corresponding probe sets were quantitatively analyzed and functionally annotated to pre-defined KEGG signaling pathways regardless of the established potential subcellular localization of the proteins isolated in MT fractions.

Obesity-related protein changes and functional annotations

We identified 221 hepatic and 44 muscle proteins that were expressed differently in control mice than in any obese mice (Tables S3 and S4). Twenty-nine of the differentially expressed hepatic proteins identified in this study were also reported as obesity-related markers in at least one of five previously published reports (Douette et al. 2005; Zhang et al. 2010; Baiges et al. 2010; Kirpich et al. 2011; Thomas et al. 2012) (Table S12). Twelve of these proteins exhibited the same direction of change as in a proteomic study by Douette et al. performed on hepatic MT fraction (Douette et al. 2005).

In the current study, proteins differentially expressed in relation to obesity, assigned to Fatty acid metabolism pathway, were upregulated in at least one obesity model and time point, and included long-chain fatty acyl-CoA synthetase (Acsm1), carnitine O-palmitoyltransferases (Cpt), trifunctional mitochondrial enzyme subunits α and β (Hadh), medium-chain specific acyl-CoA dehydrogenase (Acadm), 3-hydroxyacyl-CoA dehydrogenase (Hsd17b10), acyl-CoA oxidase (Acox1), 3-ketoacyl-CoA thiolases (Acaa), and 3,2-trans-enoyl-CoA isomerases (Dci) (Table S3). Also, among differentially expressed proteins assigned to Glycolysis/Gluconeogenesis pathway, most were upregulated, i.e., Aldob, α-enolase (Eno), aldehyde dehydrogenase (Aldh), l-lactate dehydrogenase (Ldha), and glyceraldehyde-3-phosphate dehydrogenase (Gapdh) (Table S3). Although for skeletal muscles, non-functional KEGG pathway was selected from either mRNA or protein data sets, several proteins participating in fatty acid or glucose metabolism processes were also upregulated (Table S4).

In ob/ob mice, the combination of glucose overproduction and enhanced fatty acid synthesis was observed, leading to further increase of insulin secretion and resistance (McGarry 1992). Consistently in our study, Fasn was also upregulated in young and old hyperphagic mice (mean FC = 1.67); however, it was downregulated in HFD-fed mice (mean FC = 0.58). Also, acetyl-CoA carboxylase 1 (Acaca) was downregulated in young mice fed with HFD, whereas different isoforms of acyl-CoA thioesterase (Acot) and Scd1, assigned to Biosynthesis of unsaturated fatty acids KEGG pathway, were upregulated in hyperphagic mice, in opposite to HFD-fed mice (Table S3).

For 27 obesity-related differentially expressed proteins, the direction of change in HFD-fed mice was opposite to that observed in hyperphagic mice (Table 2 and Fig. 5). STRING analysis reported these proteins as significantly networked and 13 of them were annotated to various metabolic KEGG pathways, with Peroxisome, PPAR signaling, and Metabolism of xenobiotics by cytochrome 450 as the most significantly enriched. Interaction between two cytochrome P450 components (Cyp3a11 and Cyp4a14) and Aldob are central elements of cluster network. Expression levels of eight (Scd1, Gstm1, Abcd2, Fasn, Cyp4V3, Es31, Oat, and Agxt) of these 27 proteins were changed by more than 3-fold when compared hyperphagic and HFD-fed mice. These differences in the direction of change could be due to specific diet composition, the later onset of adverse hepatic biochemical changes in HFD-induced obesity, or different regulation of some metabolic processes in obesity induced by hyperphagia and by HFD. Scd1 protein levels were 3.3-fold higher in hyperphagic mice than in HFD-fed mice, which is consistent with the Scd1-related model of obesity. A Scd1 knock-out mutation in ob/ob mice disrupts the conversion of saturated long-chain fatty acids to monounsaturated fatty acids and significantly suppresses both obesity and hepatic steatosis (Cohen et al. 2002). The transient downregulation of Scd1 induces a similar restoration of energy balance and metabolism in HFD-fed mice (Jiang et al. 2005).

Together, these results suggest that fatty acid synthesis processes might be differentially regulated in a course of obesity induced by leptin deficiency and by HFD. In opposite to HFD-fed mice, young hyperphagic mice show likely enhanced potency involved in the lipogenesis processes.

It was shown that HFD activates hepatic peroxisomal metabolism of very long-chain and branched-chain fatty acid (Kozawa et al. 2011). In Zucker rats model of obesity, peroxisomal oxidation products were shown to enter mitochondria, providing the possibility to maintain lipid disposal, even under insulin resistance condition (Noland et al. 2007). Proteomic survey in the current study confirmed the enhanced hepatic peroxisomal capacity in fatty acid oxidation. Among 221 differentially regulated proteins in liver, 24 were assigned to Peroxisome pathway and most of them were upregulated in 16-week-old mice; six of them only in hyperphagic mice (Table S3). In old animals, differential expression of proteins involved in peroxisomal oxidative activity was not indicated, confirming that metabolic processes were altered predominantly in early phase of chronic obesity.

Combined hepatic proteomic- and transcriptomic-based functional annotations indicated that in young obese mice, alterations almost exclusively reference to metabolic processes (Table S7). Among them, mainly lipid metabolism (Fatty acid metabolism, Biosynthesis of unsaturated fatty acids, PPAR signaling, and Peroxisome) and carbohydrate metabolism (Glycolysis/Gluconeogenesis) were deregulated. Both classes of processes are essentially involved in an energy production and storage, although, comparing with other nutrient groups, fatty acids are the most efficient energy source.

Enhanced fatty acids metabolism and dysregulated gluconeogenesis were previously indicated in polygenic and diet-induced obesity mouse models, both by total hepatic gene expression analysis (Becker et al. 2004) and proteomic study of adipocyte membranes (Ho et al. 2013). Similarly, a liver MT proteins study shows coordinated upregulation of bioenergetic processes during type 2 diabetes progression in rats as potential compensation for the decreased glucose metabolism by insulin resistance (Deng et al. 2010).

Defective MT oxidation of fatty acid and the production of reactive oxygen species (ROS) may lead to the hepatic but not muscle lipid accumulation (Civitarese et al. 2007; Machado et al. 2012). Feeding with HFD for 3 days in men and for 3 weeks in mice downregulates genes necessary for muscle OXPHOS and mitochondria biogenesis (Sparks et al. 2005), and prolonged consumption of a high-caloric diet has been reported as associated with a reduction of MT capacity and respiration in mouse skeletal muscle (Bonnard et al. 2008) and rat liver (Vial et al. 2011). However, other studies on obese mice showed increased rather than decreased mitochondria density, respiration, enzyme, activity and fatty acid oxidation (Hoeks et al. 2011); expansion of muscle adipose depots in obese humans did not associate with MT dysfunction (Machado et al. 2012); and in obese rats, MT dysfunction was found to be a late event (Chanséaume et al. 2007). In accordance with these results, our study revealed that in the liver, but not in skeletal muscles, changes in MT proteins and their corresponding transcripts mostly related to fatty acids metabolism pathways.

Fasting-related protein and functional changes

An overnight fast causes enhanced metabolic stress in mice associated with a significant reduction of liver weight. This effect was particularly pronounced in lean mice in which liver mass was decreased by approximately 30 % (Fig. 2b). Although not significant, a body weight reduction of approximately 15 % was also noted in both lean and obese mice, similarly as reported previously (Heijboer et al. 2005). A suppressive effect of overnight fasting on serum insulin levels was also confirmed in all animals studied (not shown). However, although insulin is responsible for restraining hepatic glucose production in the fasted state and maintaining normal blood glucose concentrations (Heijboer et al. 2005), only lean (control and young HFD-fed) mice exhibited decreased glucose levels (not shown).

Proteomic analysis revealed 107 proteins in liver samples and 35 proteins in muscle samples that distinguished all fasted from all non-fasted mice (Tables S8 and S9). The greatest number of overnight fasting-altered hepatic proteins was indicated in lean control mice comparing with any of obese ones (Table 3). Moreover, only few differentially expressed proteins were common between mice with obesity induced by HFD feeding and genetically, presumably indicating differences in molecular mechanisms underlining the acute fasting in both obesity models. Our study indicated 45 hepatic proteins with concordant expression change in at least two groups of animals (Fig. 6). Of these, Cyp4a14 and Cyp4a10 proteins, belonging to cytochrome P450 system, were consequently upregulated in all groups of animals. These proteins are encoded by well-established peroxisome proliferator-activated receptor (PPAR)α-responsive genes (Patsouris et al. 2006).

PPAR signaling pathway is one of the most frequently changed in the liver during the overnight fasting, on the mRNA level (Table S10). PPARs seem to be the key regulators of overall metabolic processes in the cell as they are implicated in metabolism of lipids, carbohydrates, and proteins. PPARs are steroid nuclear receptors that function as transcription factors that mediate fatty acids regulation of gene expression, acting as the metabolic switch between fed and fasted state. There are three isoforms of PPARs: α, β/δ, and γ, with distinct tissue distributions and specific set of activated target genes (Stienstra et al. 2007).

PPARγ is able to activate genes involved in lipogenesis in a fed state. It is abundant in adipose tissue, and plays a role in adipocyte differentiation and energy storage. Its mRNA level is increased in the adipose tissue of obese people (Vidal-Puig et al. 1997) and its activation decreased production of inflammatory cytokines in macrophages (Sharma and Staels 2007). Two other PPAR isoforms promote lipid catabolism and are found in metabolically active tissues with high fatty acids oxidation. PPARα is expressed mainly in the liver, brown adipose tissue, skeletal muscle, heart, and blood cells. It is a key regulator of energy homeostasis and plays a pivotal role in the regulation of genes involved in fatty acids oxidation, metabolism of lipoproteins, gluconeogenesis, and synthesis of ketone bodies (Mandard et al. 2004). Overall, PPARα activation decreased hepatic fat storage. During fasting, it decreased the expression of genes involved in the amino acids metabolism (Kersten et al. 2001). It was upregulated in mice fed with HFD (Kozawa et al. 2011). Beside metabolic processes, PPARα is also involved in regulation of inflammation processes. Under activation, it reduces hepatic inflammation by interference with several proinflammatory transcription factors (Delerive et al. 1999), as well as by inhibition of cytokine signaling pathways via downregulation of the IL-6 receptor (Gervois et al. 2004). PPARβ/δ is ubiquitously expressed and is less specific. Its functions are still the least defined, although recently its role in insulin secretion by pancreatic β-cell (Iglesias et al. 2012) and muscle regeneration (Bonala et al. 2012) have been indicated. Altogether, PPARs are among the most important processes which allow the adaptation of metabolic and functional responses to extracellular signals.

When assaying fasting regulated hepatic mRNAs and proteins to KEGG pathways, nine pathways were selected solely for old obese animals based on proteomic data; in turn, corresponding mRNAs were annotated to twice as more differentially regulated pathways in young and old mice, both lean and obese (Table S10). Such observations may suggest that transcriptional changes take place as a type of adaptive processes in response to overnight fasting. In skeletal muscles, overnight food deprivation, contrary to chronic obesity conditions, allow to point several KEGG pathways for both proteomic and corresponding transcriptomic data sets (Table S11). However, functional annotations were limited to only two groups of animals: young HFD-fed mice and old control mice.

In conclusion, we provide a comprehensive overview on the molecular, mitochondria-related changes that arise during obesity and food deprivation in two metabolically active tissues, liver and skeletal muscles. Expression levels of genes and proteins that accompanied the obesity and fasting states were quantitatively different in the two studied tissues, and both conditions induce more pronounced proteome changes in the liver than in the muscles.