Proteomic analysis of liver mitochondria from rats with nonalcoholic steatohepatitis

Abstract

AIM: To explore mitochondrial dysfunction in nonalcoholic steatohepatitis (NASH) by analyzing the proteome of liver mitochondria from a NASH model.

METHODS: The NASH rat model was established by feeding rats a fat-rich diet for 24 wk and was confirmed using hematoxylin and eosin staining of liver tissue and by changes in the levels of serum alanine transaminase, aspartate aminotransferase, triglyceride, total cholesterol and other markers. Liver mitochondria from each group were isolated using differential centrifugation. The mitochondrial samples were lyzed, purified and further analyzed using two-dimensional electrophoresis combined with matrix-assisted laser desorption/ionization tandem time-of-flight mass spectrometry. Bioinformatic analyses of assigned gene ontology and biological pathway was used to study functional enrichments in the abundant proteomic data.

RESULTS: Eight up-regulated and sixteen down-regulated proteins were identified that showed greater than 1.5-fold differences between the controls and the NASH group. These dysregulated proteins were predicted to be involved in different metabolic processes including fatty acid β-oxidation processes, lipid metabolic processes, cell-cycle arrest, cell polarity maintenance, and adenosine triphosphate/sex hormone metabolic processes. Novel proteins that may be involved in NASH pathogenesis including the trifunctional enzyme Hadha, thyroxine, prohibitin, aldehyde dehydrogenase ALDH1L2, UDP-glucuronosyltransferase 2B31, and carbamoyl-phosphate synthase were identified using bioinformatics tools. The decreased expression of Hadha in NASH liver was verified by Western blotting, which was used as a complementary technique to confirm the proteomic results.

CONCLUSION: This novel report on the liver mitochondrial proteome of a NASH model may provide a reservoir of information on the pathogenesis and treatment of NASH.

Key Words: Nonalcoholic steatohepatitis; Hadha; Proteomics; Rat model

Core tip: Nonalcoholic fatty liver disease (NAFLD) is a major worldwide cause of chronic liver diseases, and nonalcoholic steatohepatitis (NASH) plays a critical role as a “turning point” in the development of NAFLD. Nevertheless, the pathogenesis of NASH remains unclear, and mitochondrial dysfunction is known to be actively involved. To date, no study has reported on specific protein expression patterns in NASH mitochondria. We have, for the first time, performed a proteomic analysis of mitochondria from NASH rats, aiming to provide a protein reservoir for in-depth analyses of NASH mechanisms and for the exploration of potential therapeutics.

INTRODUCTION

Nonalcoholic fatty liver disease (NAFLD) is a major cause of chronic liver disease in North America, and its prevalence has been estimated to be as high as 35% in some populations[1]. Currently, NAFLD is considered to be a hepatic manifestation of the metabolic syndrome that develops progressively from simple steatosis to nonalcoholic steatohepatitis (NASH), fibrosis and cirrhosis[2]. Of these stages, NASH is a vital stage due to its role as the “turning point” in NAFLD development. Hepatic steatosis is known to be reversible; however, NASH is typically irreversible and the rate of progression to liver cirrhosis accelerates when NASH occurs. A recent study showed that 41% and 5.4% of NASH patients progress to liver fibrosis and end-stage liver disease, respectively[3]. NASH is therefore a major risk factor for cryptogenic liver cirrhosis[4] and requires more in-depth investigation.

Although Day and James[5] previously outlined a “two-hit hypothesis” for NAFLD, the general mechanisms underlying the pathogenesis of NASH remain unclear. Currently, oxidative and endoplasmic reticulum stress, inflammatory factor release, and in particular, hepatic mitochondrial dysfunction are considered to be actively involved in NASH[6]. Mitochondria play a major role in cell energy-generating processes and integrate several signaling pathways to control cellular life and death. Accumulating evidence has shown that NAFLD is characterized by mitochondrial alterations that depend on the activation of intracellular stress cascades or receptor-mediated pathways[7]. Furthermore, NASH is considered a mitochondrial disease[8], which reinforces the importance of mitochondrial dysfunction in its pathogenesis.

How mitochondrial dysfunction affects the initiation and progression of NASH remains unclear. Previous studies have shown that impaired adenosine triphosphate (ATP) formation and augmented generation of reactive oxygen species by the damaged respiratory chain may underlie the deterioration in mitochondrial dysfunction in NASH[9]. Technological advancements have gradually uncovered the complex mechanism of mitochondrial dysfunction in NASH. For example, proteomic techniques have provided a novel and high-throughput method for elucidating the complex pathogenesis of NAFLD[10]. Using two-dimensional electrophoresis (2-DE) and matrix-assisted laser desorption/ionization tandem time-of-flight (MALDI-TOF-TOF) mass spectrometry (MS), our group reported unique protein expression patterns in a rat model of different stages of NAFLD[11] and in a mouse model of ischemia/reperfusion injury and ischemic preconditioning[12]. Recently, we also identified serum biomarkers enabling NAFLD diagnosis using a combination of surface-enhanced laser desorption/ionization (SELDI)-TOF-MS and bioinformatics[13].

However, the specific protein expression pattern in NASH liver mitochondria has rarely been reported. Considering the importance of NASH in the clinical progression of NAFLD and the complexity of mitochondrial dysfunction in disease initiation and progression, we report the first proteomic analysis of liver mitochondria from the NASH rat model. We aimed to provide a protein reservoir for future in-depth analyses of NASH pathogenesis and progression.

MATERIALS AND METHODS

Ethics statement

This study was approved by the Review Board of the First Affiliated Hospital, School of Medicine, Zhejiang University, China. All animal studies were conducted according to the regulations and guidelines for the use and care of experimental animals of the Department of Gastroenterology, the First Affiliated Hospital, School of Medicine, Zhejiang University, China.

Establishment of the NASH rat model

Twenty four Sprague-Dawley rats (Rattus norvegicus) weighing 160-170 g were purchased from the Medical Science Institution of Zhejiang Province (Hangzhou, China) and were randomly divided into NASH (n = 12) and control (n = 12) groups, as reported previously[14]. All rats received food and water ad libitum and were maintained on a 12/12-h light/dark cycle. The control group was provided with a basic diet, whereas the NASH group was fed a fat-rich diet, as described previously[15]. After 24 wk, the rats were euthanized by femoral exsanguination, and alanine transaminase (ALT), aspartate aminotransferase (AST), triglyceride (TG), total cholesterol (TCh), and hepatic TG levels were measured. The hepatic index, which was used to describe the lipid overload in the liver, was calculated as the ratio between liver wet weight and body weight. Liver sections were stained using hematoxylin and eosin (HE) and observed for hepatic steatosis and inflammation using an Olympus microscope. The severity of hepatic injury was estimated according to the histological activation index (HAI) as described previously[16].

Isolation of mitochondria and sample preparation

Liver mitochondria from each group were isolated using differential centrifugation. First, liver sections were excised from the euthanized rats, washed with 0.25 mol/L sucrose and homogenized (1/10, w/v, 600 μg) in a JA-17 rotor using MSHE (0.22 mol/L mannitol, 0.07 mol/L sucrose, 0.5 mmol/L EGTA, 0.1% bovine serum albumin, and 2 mmol/L Hepes/KOH, pH 7.4) at 4 °C for 5 min. The supernatant was then centrifuged at 10300 g for 10 min. Next, intact purified mitochondria were isolated using Percoll (Sigma, CA, United States) to remove contaminating organelles and broken mitochondria. The pellet (mainly the mitochondrial fraction) was then re-suspended in 5 mL of MSHE (225 mmol/L mannitol, 1 mmol/L EGTA, 25 mmol/L Hepes, and 0.1% bovine serum albumin) supplemented with 20 mL of 30% Percoll. This solution was spun at 95000 g in a Hitachi RP50T rotor for 30 min. The fraction with a density of 1.052-1.075 g/mL was then collected and washed twice with MSHE at 6300 g for 10 min to remove the Percoll. Finally, the purified mitochondria were washed twice using 150 mmol/L KCl and MSHE sequentially.

The mitochondrial samples were then processed for the experiment. First, lysis buffer containing 30 mmol Tris, 8 mol Urea, and 4% CHAPS 3-[(3-cholamidopropyl)-dimethyammonio]-1-propane sulfate) was added to the microfuge tubes to re-suspend the purified mitochondria and the suspension was stored at 4 °C for at least 2 h. The mitochondrial suspension, while on ice, was sonicated intermittently for 10 s bursts followed by 10-s cooling periods. This process was repeated 15 times. After sonication, the samples were centrifuged at 12000 g at 4 °C for 30 min. The supernatant was transferred to a new tube and any pellet was discarded. Contaminants were further removed using a 2-D Clean Up Kit (GE Healthcare, United States), and the protein concentrations were calculated using a 2-D Quant Kit (GE Healthcare, CT, United States).

2-DE and image analysis

Three-hundred micrograms of protein was diluted with rehydration solution comprising 8 mol urea, 2% CHAPS, 13 mmol/L DTT, 0.5% IPG buffer, and 0.002% bromophenol blue (pH = 3-10) to a total volume of 450 mL. The protein solution was loaded onto six 24-cm Immobiline dry strips pH = 3-10 (GE Healthcare), which were used for isoelectric focusing on an Ettan IPGphor 3 IEF system (GE Healthcare) for a total of 66000 volt-hours. 2-D SDS-PAGE was performed on all six polyacrylamide gels simultaneously using an Ettan DALT Six electrophoresis unit at 5 W/strip for 45 min and then at 15 W/strip until the bromophenol blue reached the bottom of the gels. Protein spots on the gels were visualized using routine methods, and the 2-D images were scanned using a high-resolution image scanner at 300 pixels per inch. Image Master 2D Platinum 6.0 (GE Healthcare) software was used to match and analyze the protein spots. Triplicate gels for each sample were used to reduce experimental errors. An average increase or decrease greater than 1.5-fold between the NASH and control groups was used to identify the differentially expressed proteins.

In-gel digestion and MALDI-TOF-TOF/MS identification

The differentially expressed proteins were excised from silver stained gels. Each spot was destained by washing with a 1:1 solution of 30 mmol potassium ferricyanide and 100 mmol sodium thiosulfate followed by equilibration in 200 mmol ammonium bicarbonate for 20 min. After washing twice with Milli-Q water, the gel spots were dehydrated by adding acetonitrile and dried in a SpeedVac (Thermo Savant, United States) for 15 min. Subsequently, the gel spots were rehydrated in 5 μL trypsin solution (20 ng/μl in 200 mmol NH₄HCO₃) and incubated at 37 °C overnight. After eluting twice with 30 μL 50% acetonitrile and 5% trifluoroacetic acid, 1 μL peptide mixture was mixed with 1 μL cyano-4-hydroxycinnamic acid (10 mg/mL, Sigma, United States) and then saturated with 50% acetonitrile in 0.05% trifluoroacetic acid.

The mixture was analyzed using a MALDI-TOF-TOF mass spectrometer (4800 Proteomics Analyzer; Applied Biosystems, United States). Data were acquired in a positive MS reflector mode at a scan range of 800-4000 Da. Five monoisotopic precursors (S/N > 200) were selected for tandem mass analysis. To interpret the mass spectra, a combination of peptide mass fingerprints and peptide fragmentation patterns were used for protein identification against the National Center for Biotechnology Information (NCBI) non-redundant protein database using the Mascot search engine (http://www.matrixscience.com). All mass values were considered monoisotopic, and the mass tolerance was set to 150 ppm.

Western blotting and bioinformatics analyses

Hadha protein levels in liver from the NASH and control groups were examined using Western blotting with a primary mouse polyclonal antibody raised against Hadha (1:500; Abcam, United Kingdom) and an ECL chemiluminescence kit (Santa Cruz, Texas, United States). Bands were normalized by comparing to a western blot of the same samples, but probed with a mouse anti-GAPDH antibody. The sequences of Hadha and all uncoupling proteins (UCPs) were obtained from the NCBI gene and protein databases and aligned using CLUSTAL X version 2.08 using the SLIM substitution matrix and other default parameters, as described previously[17]. For more information, the gene ontology (GO) database was used to complete functional enrichment analyses of the proteomic data we obtained.

Statistical analysis

Each experiment was performed in triplicate, and the data were expressed as the mean ± SE. Student’s t-tests were used to compare two unpaired groups and were executed by SPSS 17.0. Differences were considered statistically significant at P < 0.05.

RESULTS

Establishment of the NASH rat model

The NASH rat model was established successfully as confirmed by both serological and pathological changes. After feeding the rats a 24-wk high-fat diet, HE staining of NASH liver revealed varying degrees of fat deposition and mild to moderate chronic portal and intra-acinar inflammation (Figure 1). When compared with the control group, serum ALT, AST, TG and TCh levels, hepatic TG, hepatic index, and the HAI were all significantly increased in the NASH group (Table 1).

Table 1 Comparison of hepatic and serological markers in nonalcoholic steatohepatitis liver.

Marker	Group
	NASH	Control
Hepatic index (%)	3.79 ± 0.25a	2.51 ± 0.19
ALT (IU/L)	153.51 ± 16.39a	47.44 ± 14.06
AST (IU/L)	221.63 ± 37.28a	104.31 ± 16.07
TG (mmol/L)	0.93 ± 0.12a	0.55 ± 0.11
Tch (mmol/L)	2.95 ± 0.41a	1.41 ± 0.23
Hepatic TG (mmol/L)	5.66 ± 0.71a	1.43 ± 0.14
HAI	3.69 ± 0.51a	1.31 ± 0.16

Values are expressed as the mean ± SE,

^aP < 0.05. NASH: Nonalcoholic steatohepatitis; ALT: Alanine transaminase; AST: Aspartate aminotransferase; TG: Triglyceride; TCh: Total cholesterol; HAI: Histological activation index.

Figure 1 Pathology of nonalcoholic steatohepatitis and control group liver after hematoxylin and eosin staining. A: Control group; B: Nonalcoholic steatohepatitis group.

Quantitative proteomic analysis

The liver mitochondrial protein profiles of NASH and control rats were determined using a 2-DE approach. Briefly, 800-1000 protein spots per gel were detected using the Image Master 2D Platinum 6.0 software for image analysis. Differentially expressed spots were detected using a 2-D pattern comparison between the NASH and control groups. Using Student’s t-tests, 34 protein spots showed changes of more than 1.5-fold between the NASH and control groups at a statistical significance of P < 0.05 (Figure 2). These differentially expressed protein spots were picked from the preparative gels and subjected to in-gel digestion and MALDI-TOF-TOF mass analysis. A total of eight up-regulated and sixteen down-regulated proteins were matched and identified, as displayed in Table 2. These dysregulated proteins were predicted to be involved in different metabolic processes including fatty acid β-oxidation processes, lipid metabolic processes, cell-cycle arrest, cell polarity maintenance, ATP metabolic processes, and sex hormone metabolic processes.

Table 2 Dysregulated proteins in liver mitochondria of nonalcoholic steatohepatitis rats.

	Spot no.	Protein GI	Protein name	Fold change	GO molecular function
Up-regulated	611	157822175	dehydrogenase 1 β subcomplex 10 (Ndufb10)	3.05	metabolic process/electron transport/ATP synthesis
	605	257467641	Suppressor of glucose by autophagy (SOGA)	2.97	suppressing glucogenesis
	510	29436756	Nuclear mitotic apparatus protein 1	2.68	microtubule binding/tubulin binding
	220	76779273	Hspd1 protein	2.54	ATP binding/insulin binding
	349	149029483	Mitochondrial F1 complex, alpha subunit, isoform 1	2.49	ATP synthesis/ADP binding
	51	202549	Iodothyronine 5' monodeiodinase	2.37	deiodination of thyroxine to T3
	135	148676986	Phosphoglucomutase 5	1.81	carbohydrate metabolism/Ca2+ homeostasis
	35	344244087	UDP-glucuronosyl-transferase 2B31	1.65	transferase activity/transferring hexosyl groups
Down-regulated	310	158631196	Peptidyl-prolyl cis-trans isomerase H isoform 2	-106	isomerase activity/ribonucleoprotein complex binding
	208	225690572	Leucine-rich repeat- containing protein 49 isoform 3	-106	protein-protein interaction/component of immune system
	183	300796253	Mitochondrial 10-formyltetrahydrofolate dehydrogenase (ALDH1L2)	-106	folate metabolism/CO production
	221	123244269	Microtubule-actin crosslinking factor 1	-106	ATPase activity/calcium ion binding/microtubule binding
	132	148687612	DEAH (Asp-Glu-Ala-His) box polypeptide 37, isoform CRA_b	-3.88	ATP binding/helicase and hydrolase activity
	248	60688189	Agmat protein	-2.85	agmatinase activity/hydrolase activity/metal ion binding
	244	351700869	Zinc finger protein 613	-2.52	nucleic acid binding/zinc ion binding
	258	47087119	Estradiol 17-β-dehydrogenase 8	-2.24	3-hydroxyacyl-CoA dehydrogenase activity/estradiol 17-β-dehydrogenase activity/oxidoreductase activity/testosterone 17-β-dehydrogenase (NAD+) activity
	122	6679299	Prohibitin	-2.14	histone deacetylase binding/transcription regulatory region DNA binding
	237	8393186	carbamoyl-phosphate synthase	-1.88	ATP binding/calcium ion binding/carbamoyl-phosphate synthase (ammonia) activity
	234	20304123	3-mercaptopyruvate sulfurtransferase	-1.87	3-mercaptopyruvate sulfurtransferase activity/thiosulfate sulfurtransferase activity/transferase activity
	179	149027156	Mitochondrial 3-oxoacyl-coenzyme A thiolase, isoform CRA_f	-1.85	acetyl-CoA C-acyltransferase activity/transferase activity
	260	148747393	Trifunctional enzyme subunit, alpha, mitochondrial precursor (Hadha)	-1.76	3-hydroxyacyl-CoA dehydrogenase activity
					NAD binding/enoyl-CoA hydratase activity/fatty-acyl-CoA binding/long-chain-3-hydroxyacyl-CoA dehydrogenase activity
	156	347800699	T-complex protein 1 subunit theta	-1.74	ATP binding/nucleotide binding/unfolded protein binding
	257	11968102	Ornithine aminotransferase, mitochondrial precursor	-1.60	2-oxoglutarate 5-5-aminotransferase activity/ornithine-oxo-acid transaminase activity/pyridoxal phosphate binding
	206	351715881	GRIP and coiled-coil domain	-1.55	protein binding-containing protein 2

The -10⁶-fold change denotes very low expression in the NASH group compared with the control group. GIs in bold italics are proteins that are known to be located exclusively in mitochondria. Underlined proteins are known to be located in the inner mitochondrial membrane. ATP: Adenosine triphosphate (ATP); GO: Gene ontology.

Figure 2 Representative two-dimensional electrophoresis image of proteins from the control and nonalcoholic steatohepatitis rat livers. The numbered spots denote proteins that exhibited modified expression levels in nonalcoholic steatohepatitis (NASH) liver.

The 24 identified proteins were selected using a peptide matching method. Therefore, it is not surprising that over half were similar to proteins from other Rodentia species, such as Mus musculus, Cavia porcellus, Heterocephalus glaber, and Cricetulus griseus. Because we used R. norvegicus as the NASH rat model in this study and to increase the reliability of our results, we selected seven proteins that matched proteins from R. norvegicus from the UniProt database (Table 3). Interestingly, only protein Ndufb10 (UniProt ID: D4A0T0) was up-regulated, whereas the other six proteins were down-regulated in the NASH group compared with the control group. Moreover, the GO cellular component analysis showed that only Agmatinase (UniProt ID: Q0D2L3) was located exclusively in the mitochondrion, whereas the other proteins were expressed in both the cytoplasm and mitochondria. In addition, with the exception of estradiol 17-β-dehydrogenase 8 (UniProt ID: Q6MGB5), which is hydrophobic, the other six proteins are hydrophilic. The proteomics data has provided basic information that may serve as a protein reservoir, although further functional studies on specific proteins are urgently required.

Table 3 Dysregulated proteins that match uniprot proteins from Rattus norvegicus.

	UniProt ID	Gravity1	Fold change	GO biological process
Up-regulated	D4A0T0	-0.937	3.05	Biological process
Down-regulated	Q0D2L3	-0.033	-2.85	Putrescine biosynthetic process/spermidine biosynthetic process
	Q6MGB5	0.137	-2.24	Androgen metabolic process/estrogen and fatty acid biosynthetic process/oxidation- reduction process
	P07756	-0.177	-1.88	Anion homeostasis/carbamoyl phosphate biosynthetic process/cellular response to cAMP and fibroblast growth factor stimulus/hepatocyte differentiation and liver development/metabolic process
	P97532	-0.292	-1.87	3-Mercaptopyruvate sulfurtransferase activity
	Q64428	-0.081	-1.76	Fatty acid β-oxidation/fatty acid and lipid metabolic process/oxidation-reduction process/response to insulin stimulus
	P04182	-0.113	-1.60	Ornithine metabolic process

¹Gravity values < 0 indicate hydrophilicity, whereas values > 0 indicate hydrophobicity.

Validation and bioinformatics analyses of Hadha

Among the dysregulated proteins from R. norvegicus, Hadha (UniProt ID: Q64428) was annotated as a mitochondrial inner membrane protein. Bioinformatics analyses revealed that the sequence of rat Hadha was similar to known UCPs. Because UCP2 is known to be actively involved in NAFLD due to its uncoupling activity[18], we tested Hadha expression levels in our samples using Western blotting. As shown in Figure 3, the steady-state level of hepatic Hadha was almost 50% less in the NASH group compared with the control group. This expression pattern is consistent with the 2-DE results and the peptide mass fingerprinting for Hadha, as shown in Figure 4.

Figure 3 Significantly decreased hadha levels in nonalcoholic steatohepatitis rat mitochondria using Western blotting. NASH: Nonalcoholic steatohepatitis.

Figure 4 Peptide mass fingerprinting of an in-gel tryptic digest of the Hadha protein.

DISCUSSION

NASH is considered a vital stage in NAFLD development because of its propensity to progress to advanced disease stages. Although NASH has been investigated intensively and mitochondrial dysfunction is known to be actively involved, the underlying mechanism remains unclear. Mitochondria play a major role in cellular life and death by controlling cell energy-generating processes. More importantly, the mitochondrion has its own genome, which produces organ-specific proteins and exerts vital biological functions[19]. A previous study by our group revealed the protein profiles during different stages of NAFLD[11]; however, the global mitochondrial protein expression in NASH remains unknown. In the current study, we systematically analyzed the liver mitochondrial proteome in a high-fat diet-induced NASH rat model and produced a novel protein reservoir that may help to elucidate the mechanism underlying NASH progression.

Of the 24 significantly dysregulated proteins, some have already been reported to participate in NASH progression. For example, of the up-regulated proteins, SOGA was identified as the target of adiponectin in reducing glucose production[20] and thus, may be involved in the development of NASH. Intriguingly, iodothyronine 5' monodeiodinase, which functions in the deiodination of thyroxine to T3, was significantly increased in liver mitochondria from NASH rats. A previous study also showed its potential association with lipid peroxidation, which is vital to NASH progression[21]. Taken together, these results reveal the potential involvement of thyroxine in NASH, although further investigations are needed. Furthermore, Nomura et al[22] showed enhanced ADP-ribosylation of phosphoglucomutase in patients with excess alcohol intake. Considering the similarity in pathology between alcoholic liver disease and NAFLD, its effect in NASH is also worth investigating. Finally, the increased UDP-glucuronosyltransferase 2B31 level in NASH liver was partially in line with a previous finding that UDP-glucuronosyltransferase levels in mouse liver increased in obesity and in fasting-induced steatosis[23].

Similarly, among down-regulated proteins, ALDH1L2 has been identified as the mitochondrial homolog of 10-formyltetrahydrofolate dehydrogenase, an abundant enzyme involved in folate metabolism[24]. A previous study showed that ALDH1L2 was a likely source of CO production from 10-formyltetrahydrololate in mitochondria[25]. Therefore, by influencing folate metabolism, ALDH1L2 may also be involved in NASH progression. Moreover, increased and decreased expression of prohibitin was found in cells expressing hepatitis C virus core protein[26] and in regenerating the liver[27], respectively. We observed decreased prohibitin levels in NASH liver, which may act by blocking cell proliferation and inducing apoptosis, as reported previously in human hepatoma cells[28]. In addition, depletion of carbamoyl-phosphate synthase was reported in liver during sepsis and was found to act as a marker of mitochondrial damage[29]. We also observed decreased carbamoyl-phosphate synthase in NASH rat liver, which indicated the involvement of mitochondrial damage in NASH progression. Nevertheless, many other proteins that may play a potential role in NASH are reported here for the first time, providing a novel protein reservoir for disease development that requires further confirmation.

As potentially redundant information may be produced from the study of a single protein and the potential interaction between different proteins, it is vital to enrich protein function using pathway analyses. Because we used a NASH rat model, we focused our analyses on protein information from R. norvegicus. Based on these two criteria, we selected seven of the dysregulated proteins and examined common biological pathways that they shared. Several NASH-related pathways have been reported previously including insulin binding, fatty acid and lipid metabolic processes, and oxidation reduction processes, as shown in Table 3. We also discovered several pathways that are novel including androgen/estrogen metabolic processes, cellular response to cAMP and fibroblast growth factor stimulus, and 3-mercaptopyruvate sulfurtransferase activity. These pathways may indicate that novel mechanisms are involved in the progression of NASH and require further investigation.

Hadha was among the significantly dysregulated proteins and has been previously reported to be a candidate gene for NAFLD identification[30]. UCPs are known for their role in oxidative phosphorylation and their participation in NAFLD pathogenesis[18,31]. In this study, we confirmed the significant decrease in Hadha using Western blotting, but the detailed mechanism is, as yet, unknown. Several explanations are possible: Firstly, Hadha may influence the metabolism of fatty acid via its role in dehydrating 3-hydroxyacyl-CoA. Secondly, each protein in the mitochondrial inner membrane may possess uncoupling activity, as was previously proposed by several researchers[32]. Because Hadha is located in the mitochondrial inner membrane and shares sequence similarity with UCPs, it may participate in NASH progression by influencing ATP synthesis and levels of oxidative stress.

This study has several limitations that should be acknowledged. Firstly, only 12 of the 24 selected proteins were found to be located exclusively in mitochondria (Table 2). The other 12 proteins are also found in the cytoplasm. A possible explanation may be the multi-function characteristics of these proteins. For instance, microtubule-actin crosslinking factor 1 may exert its ATPase activity in mitochondria, but its calcium ion binding function in the cytoplasm. Secondly, although the high-fat diet-induced rat model developed steatohepatitis that was morphologically similar to human NASH, it would be more convincing if the dysregulation of these proteins was verified in human patients. Thirdly, of the 24 identified proteins, only Hadha was verified, and future studies will confirm the other 23 proteins and their functions. Fourthly, the data on Hadha level is preliminary. Our result contrasted with the increase in Hadha mRNA in NAFLD human subjects as described by Kohjima et al[33]. Nevertheless, our result is at the protein level, whereas their result is at the mRNA level, which suggests the potential for post transcriptional regulation, such as miRNA regulation and methylation. Our next step will be to investigate both gene and protein expression of Hadha in larger animals and in patient samples. Finally, the 2-DE method itself has shortcomings as it does not generally detect proteins with high (> 150 kD) or low (< 10 kD) molecular masses, or very basic or hydrophobic proteins, thus limiting the proteomic coverage of most biological samples.

In summary, the differential proteomic approach that we employed enabled the identification of important changes in mitochondrial protein expression in a NASH rat model, which provides a protein reservoir for future investigation of the mechanism of NASH development. In addition, the down-regulation of Hadha in NASH liver and its sequence similarity to known UCPs implies its potential role in NASH via uncoupling activity. Nevertheless, the underlying molecular mechanism and the regulatory network involved require further investigation.

COMMENTS

Background

Nonalcoholic fatty liver disease (NAFLD) is a major cause of chronic liver disease in North America, and its prevalence has been estimated to be as high as 35% in some populations. The spectrum of NAFLD includes simple steatosis, nonalcoholic steatohepatitis (NASH), fibrosis and cirrhosis. Of these stages, NASH is critical due to its role as the “turning point” in NAFLD development, with a greater percentage of patients with NASH progressing to end-stage liver diseases of unknown etiology.

Research frontiers

Currently, oxidative and endoplasmic reticulum stress, inflammatory factor release, and in particular, hepatic mitochondrial dysfunction are considered to be actively involved in NASH. Nevertheless, how mitochondrial dysfunction affects the initiation and progression of NASH remains unclear. Previous studies have shown that impaired adenosine triphosphate formation and increased generation of reactive oxygen species by a damaged respiratory chain may underlie the deterioration of mitochondrial dysfunction in NASH. With the development of proteomics, novel protein candidates have been revealed to participate in the pathogenesis and progression of NASH. Our group has reported unique protein expression patterns in a rat model of different stages of NAFLD and in a mouse model of ischemia/reperfusion injury and ischemic preconditioning. Recently, we also identified serum biomarkers enabling NAFLD diagnosis using a combination of surface-enhanced laser desorption/ionization-TOF-MS and bioinformatics.

Innovations and breakthroughs

The authors provide, for the first time, the liver mitochondrial proteome of the NASH rat model. Furthermore, Western blotting verified a decrease in Hadha and indicates its potential involvement in the progression of NASH.

Applications

This model proteome provides a protein reservoir for further investigation of the pathogenesis and progression of NASH.

Terminology

Hadha is the official acronym for “hydroxyacyl-CoA dehydrogenase/3-ketoacyl-CoA thiolase/enoyl-CoA hydratase, alpha subunit”. This gene provides instructions for making part of an enzyme complex called mitochondrial trifunctional protein, which is essential for fatty acid oxidation.

Peer review

This study provided a protein reservoir for in-depth investigation of the mechanism of NASH. Moreover, the down-regulation of Hadha in NASH liver and its sequence similarity to known uncoupling proteins implies its potential role in NASH via uncoupling activity. All of these results were encouraging but putative; more verification of the other proteins is required and the application of this study to humans is anticipated.