Advertisement
Research ArticleArticle

RETRACTION: All-Trans-Retinoic Acid Enhances Mitochondrial Function in Models of Human Liver

Sasmita Tripathy, John D Chapman, Chang Y Han, Cathryn A Hogarth, Samuel L.M. Arnold, Jennifer Onken, Travis Kent, David R Goodlett and Nina Isoherranen
Molecular Pharmacology May 2016, 89 (5) 560-574; DOI: https://doi.org/10.1124/mol.116.103697

This article has been retracted. Please see:

Abstract

All-trans-retinoic acid (atRA) is the active metabolite of vitamin A. The liver is the main storage organ of vitamin A, but activation of the retinoic acid receptors (RARs) in mouse liver and in human liver cell lines has also been shown. Although atRA treatment improves mitochondrial function in skeletal muscle in rodents, its role in modulating mitochondrial function in the liver is controversial, and little data are available regarding the human liver. The aim of this study was to determine whether atRA regulates hepatic mitochondrial activity. atRA treatment increased the mRNA and protein expression of multiple components of mitochondrial β-oxidation, tricarboxylic acid (TCA) cycle, and respiratory chain. Additionally, atRA increased mitochondrial biogenesis in human hepatocytes and in HepG2 cells with and without lipid loading based on peroxisome proliferator activated receptor gamma coactivator 1α and 1β and nuclear respiratory factor 1 mRNA and mitochondrial DNA quantification. atRA also increased β-oxidation and ATP production in HepG2 cells and in human hepatocytes. Knockdown studies of RARα, RARβ, and PPARδ revealed that the enhancement of mitochondrial biogenesis and β-oxidation by atRA requires peroxisome proliferator activated receptor delta. In vivo in mice, atRA treatment increased mitochondrial biogenesis markers after an overnight fast. Inhibition of atRA metabolism by talarozole, a cytochrome P450 (CYP) 26 specific inhibitor, increased the effects of atRA on mitochondrial biogenesis markers in HepG2 cells and in vivo in mice. These studies show that atRA regulates mitochondrial function and lipid metabolism and that increasing atRA concentrations in human liver via CYP26 inhibition may increase mitochondrial biogenesis and fatty acid β-oxidation and provide therapeutic benefit in diseases associated with mitochondrial dysfunction.

Introduction

Vitamin A (retinol) is stored predominantly in the liver as retinyl esters. However, the active metabolite of retinol, all-trans-retinoic acid (atRA), appears to also regulate gene transcription in the liver (O'Byrne and Blaner, 2013). Most actions of atRA are mediated via interactions with retinoic acid receptors (RARs) (Wolf, 2010; Yu et al., 2012), but atRA also activates the peroxisome proliferator-activated receptor (PPAR)-β/δ (Shaw et al., 2003; Berry and Noy, 2007; Wolf, 2010). In addition, atRA induces the expression of the transcriptional repressor small heterodimer partner in the liver, resulting in broad effects on gene transcription, especially of genes responsible for bile acid metabolism and transport (Koh et al., 2014; Mamoon et al., 2014; Yang et al., 2014). The activation of multiple pathways by atRA likely leads to diverse effects on liver physiology, including altered lipid homeostasis. In humans both excessive intake of vitamin A or retinoids and vitamin A deficiency lead to hypertriglyceridemia, elevated total cholesterol, and decreased high-density lipoprotein cholesterol (Bershad et al., 1985; Lawrence et al., 2001; Staels, 2001; Brelsford and Beute, 2008; Lilley et al., 2013). However, likely because of a dichotomy of the effects of vitamin A on the liver, the endogenous role of atRA in regulating lipid and fatty acid homeostasis and role of atRA in liver disease is not well understood.

In rodents atRA appears to regulate various processes in lipid homeostasis in a species-specific manner. Data from several mouse studies suggest that atRA signaling via RAR is needed to maintain liver mitochondrial fatty acid oxidation and lipid homeostasis. Hepatocyte-specific RAR knockout mice develop microvesicular steatosis, hepatocellular carcinoma, decreased mitochondrial fatty acid oxidation, and increased peroxisomal and microsomal fatty acid oxidation presumably due to deficient atRA signaling (Yanagitani et al., 2004). Similarly, in normal mice atRA enhances fatty acid oxidation and ketogenesis via RAR activation and FGF21 induction (Amengual et al., 2012; Li et al., 2013). In agreement with a role of atRA in maintaining lipid homeostasis in the liver, decrease of hepatic atRA concentrations via inhibition of atRA synthesis in mice led to microvesicular vacuolation but without a change in liver triglycerides (Paik et al., 2014). However, in diet-induced obese mice atRA treatment decreased hepatic triglyceride content (Berry and Noy, 2009) and hepatic lipid accumulation (Kim et al., 2014). In contrast, in rats vitamin A deficiency led to increased expression of genes involved in fatty acid metabolism in the liver (McClintick et al., 2006) and decreased liver total phospholipid content and phosphatidylcholine synthesis (Oliveros et al., 2007), demonstrating opposite effects to those observed in mice. Low retinol diet also significantly decreased stellate cell free fatty acids and total lipids, whereas high retinol diet increased triglycerides, cholesteryl esters, free fatty acids, and total lipids in rat stellate cells (Moriwaki et al., 1988). How well these observations in mice or rats reflect atRA signaling in healthy or diseased human liver is not known.

Hepatic vitamin A stores are depleted in alcoholic liver disease, and vitamin A deficiency is believed to play a role in the development and progression of the disease (Bell et al., 1989; Ward and Peters, 1992; Lee and Jeong, 2012; Clugston et al., 2013). Current data also strongly suggest that altered vitamin A homeostasis and retinoid metabolism is present in nonalcoholic fatty liver disease (NAFLD) and in subsequent steatosis. In liver biopsies from patients with NAFLD, the mRNA of retinaldehyde dehydrogenases involved in atRA formation and the mRNA of cytochrome P450 (CYP) 26A1 responsible for atRA metabolism was increased (Ashla et al., 2010). In addition, in a panel of human livers with increasing level of steatosis, retinol and retinyl ester concentrations and RARβ mRNA decreased significantly with increasing steatosis (Trasino et al., 2015). Finally, in a study of patients with NAFLD and steatosis, retinoic acid (RA) concentrations were significantly lower in plasma from patients than in healthy controls (Liu et al., 2015). Yet detailed studies of retinoid signaling in models of NAFLD have not been conducted. Because NAFLD is at least in part a mitochondrial disease associated with decreased mitochondrial gene expression and function (Pessayre, 2007; Wei et al., 2008; Aharoni-Simon et al., 2011; Betancourt et al., 2014), regulation of mitochondrial fatty acid oxidation and lipid homeostasis by atRA could play a role in NAFLD progression. The aim of this study was to determine whether atRA regulates mitochondrial function and lipid homeostasis in the healthy and fatty human liver and to establish how well mouse and rat models of atRA signaling in the liver correlate with models of human liver.

Materials and Methods

Chemicals and Reagents.

Dulbecco’s modified essential medium (DMEM) with high glucose (25 mM) and fetal calf serum were obtained from Invitrogen (Carlsbad, CA). All-trans-retinoic acid (atRA) was obtained from Sigma-Aldrich, St. Louis, MO. GW0742 (PPARδ agonist) and TTNPB (RAR-pan agonist) were obtained from Tocris, Minneapolis, MN. Talarozole (R115866) was obtained from MedChem Express, Princeton, NJ. siRNAs for RARα, RARβ, and PPARδ were obtained from Ambion, Life Science Technologies, Grand Island, NY. ATP assay kit was purchased from Biovision Inc., Milpitas, CA. Primer and probe pairs for human peroxisome proliferator activated receptor gamma coactivator 1α (PGC1α) (Hs01016719_m1, FAM), PGC1β (Hs00991677_m1, FAM), nuclear respiratory factor 1(NRF1) (Hs00192316_m1, FAM), carnitine palmitoyltransferase 1α (CPT1α; (Hs00912671_m1, FAM), PPARα (Hs00947536_m1, FAM), PPARδ (Hs04187066_g1, FAM), PPARγ (Hs01115513_m1, FAM), GAPDH (Hs99999905_m1, VIC), RARα (Hs00940446_m1, FAM), RARβ (Hs00233407_m1, FAM), adipose triglyceride lipase (ATGL; Hs00386101_m1, FAM), diacylglycerol O-acyl transferase 2 (Hs01045913_m1, FAM), microsomal triglyceride transfer protein (Hs00165177_m1, FAM), and sterol regulatory element binding protein 1(SREBP1; Hs01088691_m1, FAM) were obtained from Applied Biosystems, Carlsbad, CA. Primer and probe pairs for mouse Rarβ (Mm01319677_m1, FAM), Cyp26a1 (Mm00514486_m1, FAM), Pgc1α (Mm01208835_m1, FAM), Cpt1α (Mm01231183_m1, FAM), Atgl (Mm00503040_m1, FAM), Pgc1β (Mm00504720_m1), Nrf1 (Mm01135606_m1), and β-Actin (Mm00607939_s1, FAM) were obtained from Applied Biosystems. Mitochondrial genes and nuclear gene primer set #7246 (for human) and #RR290 (for mouse) were purchased from Takara Clontech, Mountain View, CA.

HepG2 Cell Culture.

HepG2 cells were cultured in 6-well plates (Corning Life Sciences, Corning, NY) with a cell density 1×106 cells per well in DMEM with 25 mM glucose media with 10% fetal calf serum in a 5% carbon dioxide environment in a humidified incubator at 37°C. All media contained penicillin and streptomycin, and all the studies were conducted in triplicate. Lipid-loaded HepG2 cells were used as a model for fatty liver. For lipid loading, HepG2 cells were treated with 200 μM oleic acid for 3 days. Every 24 hours, 200 μM oleic acid was added to the media. The lipid loading was confirmed by Oil Red O staining (BioVision Inc.) according to manufacturer’s recommendations.

Cell Viability Assay.

HepG2 cell growth was measured using WST-1 reagent (Roche, Indianapolis, IN). Briefly, cells were grown in a 96-well plate cultured in a 5% carbon dioxide environment in a humidified incubator at 37°C. By using six replicates per treatment and a 0.1% dimethylsulfoxide (DMSO) control, the effect of atRA treatment (0.001–50 μM) on cell proliferation was measured after 72 hours of treatment. Every 24 hours, the media were removed, cells were washed with PBS, and the media containing atRA were replaced to maintain constant atRA concentrations. At 72 hours, 10 μl of WST-1 reagent were added to each well and incubated at 37°C for 1 hour to measure cell viability and growth according to manufacturer’s instructions. The absorbance was measured at 450 nm with background subtraction at 680 nm. Each treatment was normalized to the control cells, and cell viability was expressed as a percentage of the control.

Human Hepatocyte Experiments.

Cryopreserved human hepatocytes from five donors (3 women and 2 men aged 13-52 years) were obtained from Xenotech (Lenexa, KS). Hepatocytes were thawed at 37°C, placed into thawing media (K2000; Xenotech), and cell viability was confirmed using trypan blue. Hepatocytes were diluted in plating media (Williams E media plus hepatocytes supplement pack CM3000, LifeTechnologies, Grand Island, NY) at 1.2 × 106 live cells. Plating media (330 μl) containing hepatocytes were added to each well in 12 or 24-well collagen type I-coated plates (LifeTechnologies) and allowed to attach for ∼8 hours. The plating media were then replaced with Williams E media containing cell maintenance supplements (CM4000; LifeTechnologies). Hepatocytes from three donors were plated and treated with 1 μM atRA or ethanol as a vehicle control for 24 hours, and at the end of the 24-hour treatment, cells were harvested for PGC1α, PGC1β, and NRF1 mRNA measurements or for quantification of ATP production or mitochondrial DNA. Hepatocytes were treated with 1 μM atRA or ethanol as a vehicle control for 48 hours for SDHA protein quantification. For lipid loading, human hepatocytes were loaded with 600 μM oleic acid once for 24 hours and then treated with 1 μM atRA or ethanol as a vehicle control for 24 hours. At the end of the 24-hour treatment, the cells were harvested for PGC1α, PGC1β, and NRF1 mRNA measurements or analyzed for ATP production. All treatments were conducted in triplicate.

Proteomic Analysis.

Proteomic changes driven by atRA in HepG2 cells were evaluated using the PAcIFIC shotgun proteomic method (Chapman et al., 2016). HepG2 cells were treated with either atRA (10 μM) or 0.1% DMSO (vehicle control) for 72 hours and harvested. For proteomic analysis, HepG2 cells were cultured in DMEM media supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin in 10 cm2 dishes and cells were grown to 70% confluence. The cells were treated with either atRA (10 μM) or 0.1% DMSO (control) for 72 hours. The media were replaced every 24 hours with fresh treatment media containing either atRA or DMSO until harvest. Cells were harvested by washing them twice with ice cold PBS and scraping the cells into 0.5 ml of 100 mM ammonium bicarbonate for shotgun proteomics analysis. Harvested cell suspension was sonicated to lyse the cells and subsequently centrifuged at 5,000 g to pellet cellular debris. Supernatant from each sample was collected and placed in a clean Eppendorf tube on ice. Protein quantification was performed with a BCA protein quantification kit following the manufacturer’s protocol (Thermo Scientific, Waltham, MA). Protein quantity was normalized for each sample before further proteomics sample preparation. The samples with 200 μg protein per 100 μl were denatured with the addition of 6 M urea. Subsequently, samples were buffered with the addition of 7 μl 1.5 M Tris pH 8.8, reduced with 2.5 μl of 200 mM tris-2-carboxyethyl phosphine for 1 hour at 37°C, alkylated with 20 μl of 200 mM iodoacetamide for 1 hour at room temperature in the dark, and then quenched with 20 μl of 200 mM dithiothreitol. Before addition of sequencing-grade porcine trypsin (Promega, Madison, WI) at a protein to enzyme ratio of 50:1, samples were diluted with 900 μl of 50 mM ammonium bicarbonate and 200 μl of MeOH. After an overnight incubation, peptides were desalted on a Vydac C18 macrospin column (The Nest Group, Southborough, MA) according to the manufacturer's protocol. Resulting eluent was concentrated on a SPD 111V SpeedVac (Thermo Savant, San Jose, CA) and stored until further use. Peptide digests were analyzed on an LTQ Oribtrap (Thermo Fisher, San Jose, CA) coupled with a nano-electrospray ionization in positive ion mode. For the PAcIFIC acquisition, Captive spray ionization mass spectrometry was performed using a Waters NanoAquity (Milford, MA). The mass spectrometer parameters were set as defined previously (Panchaud et al., 2009, 2011). A fused silica 200 μm i.d. trapping column purchased from New Objective (Woburn, MA) and packed with 2 cm of 200 Å, 5 μm Magic C18AQ particles (Michrom, Auburn, CA) was used. Successive analytical separation was performed on a fused silica 200 μm i.d. trapping column purchased from New Objective and packed with 15 cm of 100 Å, 5 μm Magic C18AQ particles (Michrom). For each sample injection, 1 μg of the peptide sample was loaded on the trapping column at 6 μl/min with 95% water/5% acetonitrile/0.1% formic acid for 2.5 minutes. Trapped peptides were then eluted from the trapping column onto the analytical column using a variable gradient with a flow rate of 2.5 μl/min. The gradient used two mobile phase solutions: A, water/0.1% formic acid and B, acetonitrile/0.1% formic acid. For mass spectrometry, ion source conditions were optimized using the tuning and calibration solution suggested by the instrument manufacturer and under conditions recommended by the Captive spray ionization manufacturer (Bruker/Michrom). Data acquired were converted from Thermo’s RAW format to the universal mzXML format and searched against the IPI human database v3.49 (http://www.ebi.ac.uk/IPI/IPIhuman.html) using SEQUEST (Havilio et al., 2003). The precursor ion tolerance was set to 3.75, and additional search parameters for SEQUEST included trypsin enzyme specificity, cysteines modified with iodoacetamide and the variable option for methionines in reduced or oxidized form. Results were analyzed with Peptide Prophet (Keller et al., 2002; Nesvizhskii et al., 2003) ensuring that peptide hits with a probability of >0.99 were accepted and linked to protein entries. Protein identifications were assigned relative quantification values using in-house spectral counting software (Ryu et al., 2008). The relative quantification value for each protein is determined by summing all peptide tandem mass spectrometry spectra correlated to a respective protein (Liu et al., 2004). After spectral counting, the data were filtered using a Student’s t test and a P value cutoff of 0.05 to ensure that protein changes were consistent throughout the biologic triplicate and technical duplicate analyses.

mRNA Analysis and Quantitative RT-PCR.

RT-PCR was used to quantify mRNA expression of various genes (StepOnePlus, Applied Biosystems, Carlsbad, CA) as previously described (Tay et al., 2010). For mRNA extraction, 300 μl of TRI reagent (Invitrogen) was added to each well containing cultured HepG2 cells or hepatocytes after aspiration of all media, and mRNA extracted according to the manufacturer’s recommendations. Total RNA was quantified using the Nanodrop 2000c Spectrophotometer (Thermo Fisher Scientific). cDNA was synthesized from 1 µg mRNA using TaqMan Gene expression reagents (Applied Biosystems). TaqMan real-time gene expression master mix and PCR primers and fluorescent probes were obtained from Applied Biosystems. Probes were labeled with the 5′-reporter dye 5-carboxyfluorescein and a nonfluorescent black hole quencher on the 3′-end. GAPDH was used as the housekeeping gene, and all assays were done as multiplexes. No changes in GAPDH mRNA levels were observed in any of the experiments. All biological triplicate samples were analyzed in duplicate. Changes in target mRNA were measured using relative quantification (fold-difference) and the ΔΔCT method (Tay et al., 2010).

Effects of atRA on Mitochondrial Number, Function, and Protein Expression.

The effect of atRA on mitochondrial number, function, and protein expression was determined by measuring ATP production, fatty acid oxidation, quantity of mitochondrial DNA, and expression of succinate dehydrogenase subunit A (SDHA). All experiments were conducted as biological triplicates.

For Western blotting, proteins were extracted as whole cell extract from HepG2 cells treated with vehicle (ethanol) or 1 μM atRA for 48 hours using whole cell lysis buffer in the presence of protease (Roche Applied Science) and phosphatase inhibitors (1 mM β-glycerol phosphate, 2.5 mM Na-pyrophosphate, 1 mM Na3VO4). Protein concentrations were measured by BCA protein assay. Whole cell protein extracts were separated by SDS-polyacrylamide gel electrophoresis (NuPAGE-Novex 4–12% polyacrylamide Bis-Tris, Life Science Technologies) and transferred to nitrocellulose membranes. Blots were incubated with primary antibodies against SDHA and β-actin overnight at 4°C. Next day blots were washed and incubated with secondary antibody for 1 hour at room temperature. Antigen-antibody reactions were detected and quantified using LiCor Odyssey scanner and software (Licor Inc., Lincoln, NE). SDHA (product #5839) was purchased from Cell Signaling, Danvers, MA. β-Actin antibody (product #AC-15) was obtained from Abcam, Cambridge, MA. The secondary antibodies, IRDye 680 (antimouse) and IRDye 800, (anti-rabbit) were obtained from LiCor Inc. All Western blots were conducted as technical duplicates.

To measure β-oxidation, [1-14C] palmitate (American Radiolabeled Chemicals Inc.) was used. HepG2 cells were gently detached with nonenzymatic dissociation solution (Sigma) and suspended in Krebs-Ringer buffer containing 1 nM insulin and placed in vials. Wells containing cylinders of Whatman filter paper were suspended above the HepG2 cells containing solution in vials. The vials were stoppered with rubber caps, gassed with 95% O2-5% CO2, and 0.5 μCi of [1-14C] palmitate for β-oxidation was injected into the vials. After incubation for 3 hours at 37°C, 200 μl hyamine hydroxide was injected into the wells containing the filter paper, 200 μl of 1N HCl was added to the HepG2 cells, and the vials were incubated overnight at 37°C. The filter paper in the wells was then transferred to 10 ml of scintillation fluid and counted.

To measure ATP production, 105 HepG2 cells per well were plated in a 12-well plate. HepG2 cells were treated with 0.01–10 μM atRA for 24 hours. Cells were then harvested for the ATP assay, and ATP was quantified from the cell lysate using a commercially available colorimetric kit purchased from Biovision Inc. Total ATP production was normalized to milligram protein. For ATP assay in normal and lipid loaded human hepatocytes, hepatocytes were treated with 1 μM atRA for 24 hours and then cells were harvested for the ATP assay as described above.

To quantify mitochondria number, the mitochondrial DNA was quantified using qPCR. Genomic DNA from HepG2 cells and human hepatocytes was extracted by using DNAzol. Total DNA was quantified using the Nanodrop 2000c Spectrophotometer (Thermofisher Sci.). Genomic DNA (100 ng) per reaction was used to perform qPCR. The mitochondrial and nuclear genes primer set were obtained from Takara Clontech. No changes in nuclear DNA levels were observed in any of the experiments. All biologic triplicate samples were analyzed in duplicate. Changes in mitochondrial DNA copy number were measured using relative quantification (copy number) versus nuclear DNA and the ΔΔCT method.

Role of RAR and PPARδ on RA Signaling in the Liver.

To determine the nuclear receptors potentially responsible for atRA effects on mitochondrial function, HepG2 cells were first treated with 1 μM atRA, 50 nM TTNPB (RAR pan agonist), 200 nM GW0742 (PPARδ agonist) or vehicle (ethanol) for 24 hours, and the mRNA expression of PGC1α, PGC1β, and NRF1 was measured using RT-PCR as described above. Second RARα, RARβ, and PPARδ siRNA was used in normal and lipid-loaded HepG2 cells together with 1 μM atRA or vehicle treatment of 48–72 hours to determine the effects of each nuclear receptor knock down on PGC1α, PGC1β, and NRF1 mRNA and on ATP production, SDHA protein content, mitochondrial number, and β-oxidation. siRNAs for RARα (s11801), RARβ (s11803), and PPARδ (s10883) were purchased from Ambion, Thermo Fisher Scientific. Lipofectamine RNAiMAX transfection reagent (13778-075) was purchased from Invitrogen, ThermoFisher Scientific. For the siRNA experiments, 106 HepG2 cells were grown in 6-well plates as described above. Cells were transfected with either siRNA for RARα, RARβ, and PPARδ or control siRNA (scrambled RNA) at 25 pmol/well together with 7.5 μl/well of the lipofectamine RNAiMAX transfection reagent according to the manufacturer’s recommendations. Forty-eight hours after transfection, 1 μM atRA or vehicle was added to the media. After 24 hours of atRA treatment, cells were harvested, and either RNA extracted and the expression of PGC1α, PGC1β, and NRF1 measured, or ATP production, SDHA protein, β-oxidation, or mitochondrial number were measured as described above. The knock down of PPARδ and RARβ was confirmed by Western blot of PPARδ and RARβ protein expression. Antibodies for PPARδ (ab137724) and RARβ (EPR2017) were purchased from Epitomics-Abcam. Western blotting was conducted as described above using whole cell lysate.

Effect of CYP26 Inhibition on atRA Potency.

Normal and lipid loaded HepG2 cells were treated with 1 μM atRA or ethanol (vehicle) in the presence and absence of 1 μM talarozole (a CYP26 specific inhibitor) for 24 hours. Media were collected for atRA and atRA metabolite (4-OH-atRA and 4-oxo-atRA) analysis by liquid chromatography tandem mass spectrometry as previously described (Topletz et al., 2015) using an AB Sciex QTRAP 4500 mass spectrometer operated on a negative ion electrospray detection and coupled with an LC-20AD ultra-fast liquid chromatography system (Shimadzu Co., Kyoto, Japan) and Agilent Zorbax C18 column (3.5 µm, 2.1 × 100 mm). Cells were harvested for RNA extraction and mRNA expression of PGC1α, PGC1β, and NRF1 was measured as described above. In a separate study, HepG2 cells were treated with 1 μM atRA or ethanol (vehicle) in the presence and absence of 1 μM talarozole for 24 hours, and ATP production was measured as described above.

Effect of RA and Talarozole Treatment on Mitochondrial Biogenesis in Mice.

All animal experiments were approved by Washington State University Animal Care and Use Committee and were conducted in accordance with the guiding principles for the care and use of research animals of the National Institutes of Health (Bethesda, MD). Mixed background (129/C57BL/6) mice housed in a temperature- and humidity-controlled environment with food and water provided ad libitum were used for the studies. Before tissue collection, the animals were euthanized by CO2 asphyxiation, followed by cervical dislocation, and livers were removed and snap frozen in light-protected tubes until analysis.

To determine the effect of atRA treatment on mitochondrial biogenesis in mice, young adult (∼60 dpp) female mice were treated once daily intraperitoneally with 5 mg/kg atRA (Sigma Aldrich) dissolved in 15% DMSO/85% sesame oil or vehicle for 5 consecutive days. Twelve hours after the final dose, all food was removed and all mice were fasted for 12 hours. After the fast, the mice were split into fasted and refed treatment groups. Animals in the fasted group were euthanized and livers collected 24 hours after the final RA or vehicle control dose, immediately after the 12-hour fast. For the refed group, food was added to cages after the 12-hour fasting period and the animals were euthanized and livers collected 4 hours later after 4-hour exposure to food (28 hours after the final RA dose). To determine whether inhibition of CYP26 enzymes has similar effects as atRA administration on mouse liver, adult (∼90 dpp) male mice were treated twice daily (every 12 hours) for 3 consecutive days with 2.5 mg/kg talarozole (Active Biochem, Maplewood NJ) dissolved in polyethylene glycol 300 or vehicle control. Twelve hours after the final dose, the animals were euthanized and livers collected. mRNA was extracted from 50–100 mg liver tissue and qRT-PCR was used to determine changes in target mRNA as described above. β-Actin was used as the housekeeping gene. To quantify mitochondrial number, SDHA expression was determined by Western blotting as described above, and mitochondria number was quantified from genomic DNA using qPCR as described above.

Statistical Analysis.

Statistical analyses were conducted using GraphPad prism v.5 statistical software (San Diego, CA). All differences involving multiple comparisons were tested using one-way analysis of variance with Tukey test as the post hoc analysis, whereas single comparisons were done using Student’s t test. Data are expressed as mean ± S.D.; P < 0.05 was considered significant.

Results

Proteomic Analysis of atRA Signaling in the Liver.

atRA treatment had no effect on HepG2 cell proliferation or viability at atRA concentrations from 0.001 to 50 μM over the 72-hour treatment (Fig. 1A), confirming that alterations in protein or mRNA expression were not due to altered number of cells. Based on the proteomic analysis of HepG2 cells treated with 10 μM atRA for 72 hours, the expression of 59 proteins in several functional categories was significantly altered by atRA treatment (Table 1). This was a small fraction (<0.1%) of all proteins detected in the experiment. Table 1 lists those proteins for which the same peptides were detected reliably in all analyzed samples. Almost half of the proteins (24 of 59, 41%) whose expression was significantly altered were proteins involved in lipid metabolism (Fig. 1B), suggesting that one of the major effects of atRA in the liver is modulating fatty acid oxidation and lipid metabolism. Based on the proteomic analysis, atRA treatment significantly increased expression of enzymes involved in fatty acid β-oxidation such as acyl CoA synthase ligase (ACSL) 3, ACSL4, fatty acyl CoA dehydrogenase and ECH1, isocitrate dehydrogenase (IDH2), a key enzyme involved in tricarboxylic acid cycle (TCA), and ATP synthase (ATP synthase subunit α), the mitochondrial respiratory chain complex V (Fig. 1B). In agreement with proteomic findings, the mRNA levels of CPT1α, the rate limiting enzyme for fatty acid oxidation and ATGL, a triglyceride specific lipase involved in triglyceride hydrolysis were significantly increased (Fig. 1C) in HepG2 cells at 1 and 10 μM atRA treatment. Interestingly, atRA treatment at 1 μM had no effect on mRNA levels of SREBP1c and diacylglycerol O-acyl transferase 2, which are involved in fatty acid and triglyceride synthesis, nor did atRA treatment alter the mRNA levels of microsomal triglyceride transfer protein, an enzyme involved in very low density lipoprotein synthesis and secretion (Fig. 1C). However, at 10 μM, atRA increased SREBP1c mRNA (a key transcription factor for de novo fatty acid synthesis; Fig. 1C) and FASN protein, which is regulated by SREBP1 (Fig. 1B). Taken together, the proteomics and mRNA analysis suggest that atRA increases mitochondrial fatty acid oxidation in the HepG2 cells at physiologically relevant concentrations and may increase fatty acid synthesis at supraphysiologic concentrations.

Fig. 1.
Fig. 1.

atRA increases the protein and mRNA expression of enzymes involved in fatty acid oxidation in HepG2 cells. (A) Cell viability after treatment of HepG2 cells with increasing concentrations of atRA. Cell proliferation is shown as mean ± S.D., n = 6. (B) Heat map for the effects of 10 μM atRA treatment in comparison with vehicle (DMSO) on expression of enzymes involved in different metabolic pathways as determine by the proteomic analysis. (C) The effect of atRA treatment (1 and 10 μM) on the mRNA expression of enzymes involved in lipid metabolism. mRNA results are expressed as mean ± S.D., n = 3 and are representative of at least 3 separate studies. *P < 0.05 compared with vehicle-treated controls Student’s t test.

View this table:
TABLE 1

Relative proteomic analysis of the effects of atRA on protein expression in HepG2 cells

The changes were quantified via spectral counting. Significant differences are underlined. Only proteins for which the same peptide was consistently detected in all replicate control and treated samples are shown.

Effects of atRA on Mitochondrial Biogenesis Markers in Models of Healthy Human Liver.

On the basis of the proteomic data that suggested that mitochondrial fatty acid oxidation was increased, the mRNA expression of key mitochondrial biogenesis markers PGC1α, PGC1β, and NRF1 was measured in HepG2 cells treated with 1 μM atRA for 24 hours. The mRNA of PGC1α, PGC1β, and NRF1 was increased two- to threefold (P < 0.05) in the HepG2 cells suggesting increased mitochondrial biogenesis in response to atRA treatment (Fig. 2A). In agreement with the mRNA data, ATP production was significantly increased after treatment with 1 μM atRA (Fig. 2B), suggesting mitochondrial function was also increased.

Fig. 2.
Fig. 2.

atRA increases mitochondrial biogenesis, ATP production, and fatty acid β-oxidation in normal HepG2 cells. (A) The effect of atRA (1 μM) and talarozole cotreatment (TLZ) on mRNA expression of mitochondrial biogenesis genes. (B) The effects of atRA and TLZ cotreatment on ATP production. (C and D) Induction of SDHA protein after treatment with atRA as measured by Western blotting. (C) Quantification of the representative Western blot shown in (D). The blot shows analysis of three independent biologic replicates. (E) The effects of atRA treatment on palmitate oxidation as a measure of increased β-oxidation. (F) The effect of atRA treatment on mitochondrial DNA copy number. All experiments were conducted as described in Materials and Methods. Results are expressed as mean ± S.D., n = 3. Comparisons between vehicle (ethanol)- and atRA-treated cells were made by Student’s t test: *P < 0.05 compared with vehicle-treated controls. #P < 0.05 compared with 1 μM atRA-treated controls.

Considerable depletion of atRA was observed in the treated cells at 24 hours of treatment with 1 μM atRA. atRA concentrations decreased by 87% to 129 ± 19.15 nM after 24 hours (Table 2). To test whether inhibition of CYP26 (the main enzyme metabolizing atRA) increases atRA concentrations and potency, HepG2 cells were treated with atRA together with talarozole, a CYP26 inhibitor, and atRA concentrations and mitochondrial biogenesis markers and ATP production were measured. In the presence of talarozole, atRA concentrations were 2.6-fold higher at the end of the 24-hour treatment than when cells were treated with atRA alone, whereas the concentrations of the metabolites of atRA were decreased upon talarozole treatment (Table 2). The increased exposure of the cells to atRA was reflected in two- to threefold higher PGC1α, PGC1β, and NRF1 mRNA induction in the presence of talarozole together with atRA compared with atRA alone (Fig. 2A). Talarozole treatment with atRA also resulted in approximately twofold greater ATP production in HepG2 cells compared with atRA treatment alone (Fig. 2B). Surprisingly, talarozole alone also increased ATP production, PGC1β, and NRF1 mRNA expression but not PGC1α mRNA.

View this table:
TABLE 2

Cell media atRA, 4-OH-atRA, and 4-oxo-atRA concentrations at the end of the 24-hour treatment with vehicle, atRA, talarozole, or the combination of atRA and talarozole

To confirm that the observed changes in the mRNA of mitochondrial biogenesis genes after atRA treatment result in improved mitochondrial function and increased mitochondrial number and fatty acid oxidation, the protein expression of SDHA, a marker for increase in mitochondrial content was measured together with mitochondrial DNA and fatty acid β-oxidation. In HepG2 cells 1 μM atRA treatment increased SDHA protein expression twofold (Fig. 2, C and D) and 1-10 μM atRA treatment increased fatty acid β-oxidation two- to threefold compared with vehicle-treated cells (Fig. 2E). Mitochondrial DNA was increased four- to fivefold after atRA treatment (Fig. 2F).

To further translate the findings from HepG2 cells to human liver, the effect of atRA on mitochondrial function and biogenesis was evaluated in human hepatocytes from three donors. Similar to the HepG2 cells, in human hepatocytes the mRNA expression of PGC1α, PGC1β, and NRF1 was increased two- to fivefold (P < 0.05) compared with vehicle treated group (Fig. 3A). In addition, atRA increased ATP production by two- to threefold when compared to vehicle-treated human hepatocytes. atRA treatment also increased SDHA protein expression two- to threefold (Fig. 3C) and mitochondrial DNA two- to sevenfold (Fig. 3D), demonstrating that atRA increases mitochondria number and function in human hepatocytes.

Fig. 3.
Fig. 3.

atRA increases mitochondrial biogenesis and ATP production in normal human hepatocytes. (A) Induction of the mRNA expression of key mitochondrial biogenesis genes (PGC1α, PGC1β, and NRF1) by 1 μM atRA. All data are shown as fold increase in comparison with vehicle control. (B) Induction of ATP production by 1 μM atRA in three independent human hepatocyte donors. (C) Induction of SDHA protein in three individual human hepatocyte donors after treatment with 1 μM atRA as measured by Western blotting. The quantification of the replicate Western blots and a representative blot of duplicate treatments in three donors is shown. (D) The effect of atRA on mitochondrial DNA copy number in three human hepatocyte donors. All experiments were conducted as described in Materials and Methods. Results are expressed as mean ± S.D., n = 3. Comparisons between vehicle (ethanol)- and atRA-treated cells were made by Student’s t test: *P < 0.05 compared with vehicle-treated controls.

To explore the potential nuclear receptors responsible for the effects of atRA on mitochondrial biogenesis, the effect of atRA on the mRNA expression of RARs and PPARs was tested. atRA significantly increased the mRNA expression of RARβ and PPARδ but had no effect on RARα, PPARα, or PPARγ (Supplementary Fig. 1A). This suggests that the effect of atRA on mitochondrial biogenesis may be a combination of RARβ and PPARδ activation. The role of RARβ and PPARδ in the induction of the genes and processes of interest was further studied using selective RAR and PPAR agonists in comparison to atRA. Although the PPARδ agonist GW0742 induced PGC1β and NRF1, PGC1α was not changed. In contrast the RAR pan agonist, TTNPB had no effect on PGC1α, PGC1β, and NRF1 mRNA expression alone or in combination with atRA (Supplemental Fig. 1B). Interestingly, when GW0742 was used together with atRA, the mRNA induction of PGC1β and NRF1 was greater than with atRA alone, whereas the induction of PGC1α was similar. Together these results suggest that RARs are not involved in the induction of mitochondrial biogenesis or function by atRA.

To further explore the role of specific nuclear receptors in the increase in mitochondrial biogenesis markers by atRA, RARα, RARβ, and PPARδ siRNAs were used. The RARα, RARβ, and PPARδ siRNA decreased the corresponding nuclear receptor mRNA levels by 75, 80, and 60%, respectively (Supplemental Fig. 2A). The RARβ and PPARδ siRNA also decreased RARβ and PPARδ protein expression respectively by 50–60% in atRA treated cells compared with atRA plus scrambled siRNA treated controls (Supplemental Fig. 2, B–D). In contrast, RARβ siRNA had no effect on the induction of PPARδ by atRA and PPARδ siRNA had no effect on the induction of RARβ by atRA (Supplemental Fig. 2, B–D). RARα and RARβ siRNA had no effect on the induction of PGC1α, PGC1β, and NRF1 mRNA by atRA (Fig. 4A). In contrast, PPARδ siRNA significantly decreased the induction of PGC1α, PGC1β, and NRF1 mRNA by atRA compared with scrambled siRNA control (Fig. 4A), suggesting that these effects of atRA in the HepG2 cells are mediated by PPARδ. In agreement with the role of PPARδ in the regulation of mitochondrial function and biogenesis, PPARδ siRNA eliminated the increase in ATP production by atRA as well as the induction of SDHA protein expression and increase in β-oxidation (Fig. 4, B–E). PPARδ siRNA also significantly deceased the atRA-mediated induction of mitochondrial DNA (Fig. 4F).

Fig. 4.
Fig. 4.

The increase of mitochondrial biogenesis, fatty acid β-oxidation, and ATP production by atRA treatment in HepG2 cells requires PPARδ. (A) The influence of RARα, RARβ, and PPARδ siRNA on the induction of mitochondrial biogenesis genes by atRA. (B) The effect of RARβ and PPARδ siRNA on the induction of ATP production by atRA (C and D). The effect of RARβ and PPARδ siRNA on the induction of SDHA protein expression by atRA as measured by Western blotting. (C) The quantification of replicate Western blots and (D) representative Western blot. The representative blot shows analysis of biologic duplicates for each treatment. (E and F) The effect of RARβ and PPARδ siRNA on fatty acid β-oxidation (E) and induction of mitochondrial DNA copy number (F) by atRA in HepG2 cells. Results for the mRNA, mitochondrial DNA, fatty acid oxidation, and ATP production studies are presented as fold change compared with vehicle (ethanol). Results for protein studies are presented as ratio of target protein over loading control β-actin. Results are expressed as mean ± S.D., n = 3. *P < 0.05 compared with the vehicle. #P < 0.05 compared with the scrambled siRNA plus atRA-treated group.

Effects of atRA on Mitochondrial Biogenesis Markers in Models of Human Fatty Liver.

Lipid-loaded HepG2 cells and primary human hepatocytes were used as a model of human fatty liver. Lipid accumulation was confirmed via Oil Red O staining (Fig. 5A). In lipid-loaded HepG2 cells, PGC1α, PGC1β, and NRF1 mRNA were increased two- to threefold (P < 0.05) in response to atRA treatment (Fig. 5B), a finding in excellent agreement with the normal HepG2 cells. In lipid-loaded HepG2 cells, 1 μM atRA treatment resulted in a dose-dependent increase in ATP production (Fig. 5C). Similarly in lipid-loaded HepG2 cells atRA treatment also increased SDHA protein expression twofold (Fig. 5D). The effect of talarozole on the potency of atRA as an inducer of mitochondrial biogenesis genes and ATP production was also similar in lipid-loaded HepG2 cells (Fig. 5, B and C) as in normal HepG2 cells. The induction of PGC1α, PGC1β, and NRF1 mRNA and ATP production was significantly greater when lipid-loaded HepG2 cells were cotreated with talarozole and atRA than with atRA alone at 100 nM (Fig. 5, B and C). In excellent agreement with the lipid-loaded HepG2 cell data and normal human hepatocyte experiments, in lipid-loaded human hepatocytes from three donors, the mRNA expression of PGC1α, PGC1β and NRF1 was increased two- to fourfold (P < 0.05) and ATP production was increased two- to threefold compared with vehicle-treated group in each of the three donors (Fig. 5, E and F).

Fig. 5.
Fig. 5.

atRA increases mitochondrial biogenesis and ATP production in lipid loaded HepG2 cells and human hepatocytes. (A) Oil Red O staining of control-treated and lipid-loaded HepG2 cells. Oil Red O staining was measured by quantifying the absorbance at 492 nM wavelength. The measured OD value for vehicle treated cells was 1.203 and for oleic acid loaded cells 2.642. (B) Induction of the mRNA expression of key mitochondrial biogenesis genes (PGC1α, PGC1β, and NRF1) by atRA, talarozole, or their combination in the lipid-loaded HepG2 cells. (C) Induction of ATP production by atRA, talarozole (TLZ), or their combination in lipid-loaded HepG2 cells. (D) Induction of SDHA protein expression after treatment of lipid-loaded HepG2 cells with atRA as measured by Western blotting. The quantification of the replicate Western blots and a representative blot of triplicate treatments is shown. (E) Induction of the mRNA expression of key mitochondrial biogenesis genes by 1 μM atRA in lipid loaded human hepatocytes from three individual donors. (F) Induction of ATP production by 1 μM atRA in lipid-loaded human hepatocytes from three individual donors. Results for mRNA studies are presented as fold change in mRNA abundance relative to vehicle (ethanol). mRNA results, ATP assay, and protein expression results are expressed as mean ± S.D., n = 3. Comparisons between vehicle (ethanol)- and atRA-treated cells were made by Student’s t test: *P < 0.05 versus vehicle. #P < 0.05 versus 1 μM atRA. †P < 0.05 versus 100 nM atRA.

Similar to normal HepG2 cells, atRA significantly increased the transcription of both RARβ and PPARδ in lipid-loaded HepG2 cells (Fig. 6A). Interestingly, in lipid-loaded HepG2 cells, RARα was also induced by atRA treatment. When siRNA for each of the nuclear receptors was used in the lipid-loaded HepG2 cells, only PPARδ knock down attenuated the atRA-mediated induction of PGC1α, PGC1β, and NRF1 mRNA (Fig. 6, B–D). RARβ and RARα knock down had no effect on the induction of PGC1α, PGC1β, and NRF1 mRNA by atRA in the lipid-loaded cells.

Fig. 6.
Fig. 6.

atRA requires PPARδ in the regulation of mitochondrial biogenesis in lipid-loaded HepG2 cells. (A) The induction of nuclear receptor mRNA by 1 μM atRA in lipid-loaded HepG2 cells. (B) The influence of RARα siRNA on the induction of mitochondrial biogenesis genes by atRA. (C) The influence of RARβ siRNA on the induction of mitochondrial biogenesis genes. (D) The influence of PPARδ siRNA on the induction of mitochondrial biogenesis genes. Results are presented as fold change in mRNA abundance relative to vehicle (ethanol). mRNA results are expressed as mean ± S.D., n = 3, and are representative of at least 3 separate studies. *P < 0.05 versus vehicle. # P < 0.05 versus scrambled plus atRA.

Effect of Increased atRA Concentrations on the Expression of Mitochondrial Biogenesis In Vivo in Mouse Liver.

To evaluate whether the in vitro findings of the effects of atRA on hepatic mitochondrial biogenesis in human liver models translate to in vivo, mice were treated with vehicle or atRA for 5 days, and hepatic mRNA expression of Pgc1α, Pgc1β, Nrf1, Cpt1α, and Atgl was measured after an overnight fast or after feeding (Fig. 7). atRA treatment resulted in up to sixfold (P < 0.05) increase in the mRNA of Pgc1α, Pgc1β, Nrf1, and Cpt1α (Fig. 7A) but had no effect on Atgl expression in fasted condition. The induction of Pgc1α, Pgc1β, Nrf1, and Cpt1α was accompanied with increased mitochondrial DNA and SDHA expression in the mouse liver after atRA treatment in fasted condition (Fig. 7, B–D). However, 4 hours later, after refeeding of the mice, atRA no longer affected any of the above mentioned markers (Fig. 7, A–D). The induction of the marker genes of RAR activation, Cyp26a1 and Rarβ, was much more modest than the induction of mitochondrial biogenesis marker genes. Rarβ was induced approximately 2.5-fold (P < 0.05), but there was no significant increase in Cyp26a1 mRNA (1.7-fold, P > 0.05) after atRA treatment (Fig. 7E). Unlike mitochondrial biogenesis markers, there was no difference in Rarβ induction between fed and fasted groups (Fig. 7), a finding in agreement with different induction mechanisms and mRNA half-lives for these target genes.

Fig. 7.
Fig. 7.

Effect of atRA and talarozole treatment on mitochondrial biogenesis markers and fatty acid oxidation genes in vivo in mice. (A) Hepatic mRNA expression of Pgc1α, Pgc1β, Nrf1, Atgl and Cpt1α in mice treated with atRA. (B) Quantification of mitochondrial DNA in atRA and vehicle-treated mice after an overnight fast (fasted group) or 4 hours later after feeding. (C and D) Induction of liver SDHA protein expression in atRA and vehicle-treated mice as measured by Western blotting. (C) Quantification of the representative Western blot shown in (D). The blot shows analysis of three individual mice for each treatment. (E) Hepatic mRNA expression of Rarβ and Cyp26a1 in atRA- or vehicle-treated mice. (F) Induction of hepatic mitochondrial biogenesis genes in mice after talarozole (TLZ) treatment compared with vehicle-treated mice. All mRNA data are presented as fold change in mRNA abundance relative to control mice. β-Actin was used as the housekeeping gene. For protein data, results are expressed as ratio of SDHA to β-actin. Results are expressed as mean ± S.D., n = 4. Comparisons between groups were made by Student’s t test: *P < 0.05 versus vehicle or control.

To increase the exposure of endogenous atRA in the liver, male C57BL/6X129 mice were treated with talarozole, an inhibitor of Cyp26, and mRNA expression of key mitochondrial biogenesis genes was measured. Similar to atRA treatment, talarozole treatment increased mRNA expression of Pgc1β and Nrf1 five- to tenfold (P < 0.05), but the fourfold increase in Pgc1α mRNA was not significant (P > 0.05; Fig. 7F).

Discussion

Treatment of most cancer cell lines in vitro or neuroblastoma and acute promyelocytic leukemia in vivo with RA leads to apoptosis and cell cycle arrest manifesting the classic RAR-mediated effects of atRA regulating cell proliferation. Interestingly, in models of human liver, treatment with atRA had no effect of cell viability or proliferation at concentrations that typically lead to cell cycle arrest and apoptosis. This result led us to evaluate atRA signaling in HepG2 cells using a proteomic approach. The proteomic studies showed that instead of cell cycle regulation, atRA treatment predominantly altered the expression of proteins responsible for regulating lipid and fatty acid metabolism in the human liver models. These findings led to the identification of a novel effect of atRA signaling in the liver. The indication of the proteomic study that atRA alters lipid and fatty acid metabolism was intriguing, because use of retinoids and RA isomers isotretinoin (13-cisRA) and alitretinoin (9-cisRA) in dermatology and in cancer treatment is associated with changes in lipid and lipoprotein metabolism, including hypertriglyceridemia, hypercholesterolemia, increased total cholesterol, and decreased high-density lipoprotein-cholesterol in about one-third of the treated patients (Lawrence et al., 2001; Staels, 2001; Brelsford and Beute, 2008; Lilley et al., 2013). The induction of SREBP1 mRNA observed in this study after treatment with 10 μM atRA is in agreement with the clinical findings, as the 10 μM atRA treatment closely mimics the RA exposures observed after retinoid therapy. The 1 μM treatments used in the subsequent experiments in this study are expected to closely mimic endogenous atRA concentrations in the liver, because the measured media concentrations during treatment were approximately 100 nM. The low media concentrations measured are likely due to extensive metabolism of atRA in the cells as reported previously (Topletz et al., 2015) and partitioning and nonspecific binding of atRA into cells. Strikingly, at the concentrations of atRA (100 nM and 1 μM) that reflect changes in endogenous atRA concentrations, SREBP1 mRNA was not increased despite consistent changes in other proteins and mRNA. This suggests that at different concentrations of atRA, different nuclear receptors and signaling pathways are activated, leading to different effects on lipid metabolism.

In addition to SREBP1 induction, the proteomic study identified an increase in a complement of proteins that contribute to fatty acid oxidation, mitochondrial function, and mitochondrial energy production, leading to a hypothesis that atRA increases mitochondrial biogenesis and fatty acid β-oxidation in the liver. The subsequent mechanistic studies shown here consistently support a beneficial regulatory role of endogenous atRA on lipid homeostasis. The studies show that atRA increases mitochondrial number (mitochondrial DNA and SDHA expression), mitochondrial function (ATP production), and fatty acid β-oxidation (14C palmitate oxidation) in normal and fatty human liver. The collective pathways and proteins affected by atRA treatment are summarized in Fig. 8. It should, however, be noted that the increase observed in palmitate oxidation may also be partially due to increased peroxisomal fatty acid oxidation, and further studies are needed to characterize the complete effects of atRA on liver lipid homeostasis.

Fig. 8.
Fig. 8.

Mitochondrial biogenesis and fatty acid oxidation pathways induced by atRA in the liver models and a proposed molecular mechanism for the observed effects. (A) Collective representation of fatty acid metabolism and mitochondrial function proteins and mRNA that were induced by atRA in the models of human liver. Green boxes and circles indicate proteins and mRNA that were shown to increase in expression in response to atRA treatment. (B) Proposed molecular mechanism of atRA regulation of mitochondrial function.

Previous mouse studies suggested that atRA increases fatty acid oxidation and decreases hepatic lipid accumulation predominantly via RAR-mediated processes. atRA has been shown to suppress PPARγ2 expression and decrease hepatic lipid accumulation in diet-induced obese mice by activating RARα (Kim et al., 2014). After adenoviral overexpression of RARβ, atRA treatment increased the expression of fatty acid oxidative genes CPT1α and MCAD in mouse liver and in HepG2 cells (Li et al., 2013). atRA also enhanced fatty acid oxidation and ketogenesis via RAR activation and FGF21 induction in HepG2 cells (Amengual et al., 2012; Li et al., 2013). Although the data presented in this study in human hepatocytes, HepG2 cells, and in mice are in agreement with a role of atRA at endogenous concentrations in increasing fatty acid β-oxidation, ATP production, mitochondrial function, and mitochondrial biogenesis, the data do not support a role of RARs in how atRA regulates these processes. While atRA clearly activates RARs in the HepG2 cells based on the observed induction of RARβ mRNA, the siRNA studies and chemical agonist experiments show that the regulation of mitochondrial function by atRA is likely mediated by PPARδ. This conclusion is supported by the fact that only PPARδ siRNA decreased the atRA-mediated induction of mitochondrial biogenesis genes, mitochondrial DNA content, ATP production, and fatty acid oxidation. Importantly, PPARδ silencing had similar consistent effects on all markers of mitochondrial biogenesis and fatty acid β-oxidation, suggesting that atRA mediates all of these processes via similar mechanism. The proposed mechanism via PPARδ is shown in Fig. 8. The in vivo studies in mice are in agreement with this interpretation of predominant role of PPARδ, because atRA treatment resulted in only modest effects on the mRNA levels of RAR target genes CYP26A1 and RARβ but had significant effects on the expression of mitochondrial biogenesis genes that are believed to be PPARδ targets. Yet, it is possible that the modest effects on CYP26A1 and RARβ mRNA in mice are due to the short mRNA half-life of these genes and atRA because the induction was analyzed 24–28 hours after last atRA dose. The findings of the role of PPARδ in RA signaling are in agreement with previous studies that documented that cellular retinoic acid binding protein II and fatty acid binding protein-5 direct atRA to the nuclear receptors RAR and PPARβ/δ, respectively, and that the pleotropic effects of atRA are due to activation of both of these receptors (Napoli, 1996; Schug et al., 2007; Berry and Noy, 2012; Bonet et al., 2012; Obrochta et al., 2015). However, we also show that atRA treatment increases PPARβ/δ mRNA, and it is possible that this induction in part or solely contributes to the observed effects on mitochondrial biogenesis without atRA binding directly to PPARβ/δ.

Altered vitamin A homeostasis is well documented in alcoholic liver disease, but recent reports associate altered vitamin A homeostasis also with NAFLD and steatosis (Ashla et al., 2010; Kim et al., 2014; Liu et al., 2015; Trasino et al., 2015). For example, human liver biopsies of NAFLD and nonalcoholic steatohepatitis patients demonstrated increased expression of the main enzyme metabolizing atRA, CYP26A1, suggesting that atRA is depleted in livers of these patients (Ashla et al., 2010). This is likely deleterious in NAFLD, because we show that atRA increases mitochondrial biogenesis and fatty acid oxidation in models of fatty liver and NAFLD. Because NAFLD is associated with decreased mitochondrial function (Sookoian et al., 2010; Aharoni-Simon et al., 2011; Handa et al., 2014), these data suggest that decreasing atRA concentrations in NAFLD may contribute to decreasing mitochondrial function and disease progression and that increasing atRA concentrations in NAFLD may be beneficial and result in improvement of mitochondrial function and fatty acid oxidation. Such increases in endogenous atRA concentrations could be achieved for example via inhibition of the CYP26 enzymes as shown here with talarozole treatment in HepG2 cells and in mice. In the HepG2 cells, talarozole decreased atRA metabolism increasing atRA concentrations and thus potentiating the effects of atRA on mitochondrial biogenesis. Similarly, in vivo in mice, talarozole treatment increased the expression of mitochondrial biogenesis genes, suggesting that increasing endogenous atRA concentrations in the liver is sufficient to obtain the effects observed in HepG2 cells and human hepatocytes after atRA treatment. Consistent with the human liver models, atRA treatment increased the expression of PGC1α, PGC1β, NRF1, and CPT1α under fasting conditions in mice and increased mitochondrial number and SDHA protein. As expected from the energy needs of fed mice, after refeeding the mice had low expression of PGC1α, PGC1β, NRF1, ATGL, and CPT1α compared with fasted mice, and atRA did not have equivalent effects on mitochondrial number and function after feeding as in the fasted mice. Yet, it is also possible that because of the longer time interval between last atRA dose and collection of samples in the fed mice compared with fasted, the effects of atRA were no longer observed in the fed mice. Additional studies are needed to fully explore the differences in atRA signaling under fasted and fed conditions.

In summary, this study shows that atRA increases mitochondrial biogenesis, fatty acid β-oxidation, expression of mitochondrial respiratory chain components, and ATP production in models of healthy and fatty human liver and in vivo in mouse liver, and these effects are predominantly regulated through PPARδ signaling instead of RAR activation. Moreover, increasing atRA concentration by CYP26 inhibition potentiated atRA effects on mitochondrial biogenesis and ATP production in cell models and in vivo in mice. Taken together the increased expression of CYP26A1 in NAFLD and the observed efficacy of a CYP26 inhibitor in mice and in human liver models suggest that CYP26 inhibition would be beneficial in NAFLD and result in increased mitochondrial fatty acid oxidation and mitochondrial biogenesis.

Acknowledgments

The authors thank Dr. Alan Chait of NORC, University of Washington, for advice on the fatty acid oxidation assay and Faith Stevison, MS, for skillful assistance with retinoid analysis.

Authorship Contributions

Participated in research design: Tripathy, Chapman, Hogarth, and Isoherranen.

Conducted experiments: Tripathy, Chapman, Han, Hogarth, Arnold, Onken, and Kent.

Contributed new reagents or analytical tools: Chapman and Goodlett.

Wrote or contributed to the writing of the manuscript: Tripathy, Chapman, Han, Hogarth, and Isoherranen.

Footnotes

    • Received February 6, 2016.
    • Accepted February 25, 2016.

  • This research was supported by grants from National Institutes of Health National Institute of General Medical Sciences [Grants R01 GM111772, R01 GM081569, and R01 GM081569-S1] and National Institute of Diabetes and Digestive and Kidney Diseases [Grant P30 DK035816].

  • dx.doi.org/10.1124/mol.116.103697.

  • ↵Embedded ImageThis article has supplemental material available at molpharm.aspetjournals.org.

Abbreviations

ACSL
acyl CoA synthase ligase
ATGL
adipose triglyceride lipase
atRA
all-trans-retinoic acid
CPT1α
carnitine palmitoyltransferase 1α
CYP
cytochrome P450
DMEM
Dulbecco’s modified essential medium
DMSO
dimethylsulfoxide
ECH1
enoyl CoA hydratase 1
FASN
fatty acid synthase
FGF21
fibroblast growth factor 21
IDH2
isocitrate dehydrogenase
NAFLD
nonalcoholic fatty liver disease
NRF1
nuclear respiratory factor 1
PGC1α
peroxisome proliferator activated receptor gamma coactivator 1α
PGC1β
peroxisome proliferator activated receptor gamma coactivator 1β
PPAR
peroxisome proliferator activated receptor
RA
retinoic acid
RAR
retinoic acid receptor
SDHA
succinate dehydrogenase subunit A
SREBP1
sterol regulatory element binding protein 1

References

  1. ↵
    1. Aharoni-Simon M,
    2. Hann-Obercyger M,
    3. Pen S,
    4. Madar Z, and
    5. Tirosh O
    (2011) Fatty liver is associated with impaired activity of PPARγ-coactivator 1α (PGC1α) and mitochondrial biogenesis in mice. Lab Invest 91:1018–1028.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Amengual J,
    2. Petrov P,
    3. Bonet ML,
    4. Ribot J, and
    5. Palou A
    (2012) Induction of carnitine palmitoyl transferase 1 and fatty acid oxidation by retinoic acid in HepG2 cells. Int J Biochem Cell Biol 44:2019–2027.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Ashla AA,
    2. Hoshikawa Y,
    3. Tsuchiya H,
    4. Hashiguchi K,
    5. Enjoji M,
    6. Nakamuta M,
    7. Taketomi A,
    8. Maehara Y,
    9. Shomori K,
    10. Kurimasa A,
    11. et al.
    (2010) Genetic analysis of expression profile involved in retinoid metabolism in non-alcoholic fatty liver disease. Hepatol Res 40:594–604.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bell H,
    2. Nilsson A,
    3. Norum KR,
    4. Pedersen LB,
    5. Raknerud N, and
    6. Rasmussen M
    (1989) Retinol and retinyl esters in patients with alcoholic liver disease. J Hepatol 8:26–31.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Berry DC and
    2. Noy N
    (2007) Is PPARbeta/delta a retinoid receptor? PPAR Res 2007:73256. doi: 10.1155/2007/73256.
    OpenUrlPubMed
  6. ↵
    1. Berry DC and
    2. Noy N
    (2009) All-trans-retinoic acid represses obesity and insulin resistance by activating both peroxisome proliferation-activated receptor beta/delta and retinoic acid receptor. Mol Cell Biol 29:3286–3296.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Berry DC and
    2. Noy N
    (2012) Signaling by vitamin A and retinol-binding protein in regulation of insulin responses and lipid homeostasis. Biochim Biophys Acta 1821:168–176.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Bershad S,
    2. Rubinstein A,
    3. Paterniti JR,
    4. Le NA,
    5. Poliak SC,
    6. Heller B,
    7. Ginsberg HN,
    8. Fleischmajer R, and
    9. Brown WV
    (1985) Changes in plasma lipids and lipoproteins during isotretinoin therapy for acne. N Engl J Med 313:981–985.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Betancourt AM,
    2. King AL,
    3. Fetterman JL,
    4. Millender-Swain T,
    5. Finley RD,
    6. Oliva CR,
    7. Crowe DR,
    8. Ballinger SW, and
    9. Bailey SM
    (2014) Mitochondrial-nuclear genome interactions in non-alcoholic fatty liver disease in mice. Biochem J 461:223–232.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Bonet ML,
    2. Ribot J, and
    3. Palou A
    (2012) Lipid metabolism in mammalian tissues and its control by retinoic acid. Biochim Biophys Acta 1821:177–189.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Brelsford M and
    2. Beute TC
    (2008) Preventing and managing the side effects of isotretinoin. Semin Cutan Med Surg 27:197–206.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Chapman JD,
    2. Edgar JS,
    3. Goodlett DR, and
    4. Goo YA
    (2016) Use of captive spray ionization to increase throughput of the data-independent acquisition technique PAcIFIC. Rapid Commun Mass Spectrom. in press.
  13. ↵
    1. Clugston RD,
    2. Jiang H,
    3. Lee MX,
    4. Berk PD,
    5. Goldberg IJ,
    6. Huang LS, and
    7. Blaner WS
    (2013) Altered hepatic retinyl ester concentration and acyl composition in response to alcohol consumption. Biochim Biophys Acta 1831:1276–1286.
    OpenUrl
  14. ↵
    1. Handa P,
    2. Maliken BD,
    3. Nelson JE,
    4. Morgan-Stevenson V,
    5. Messner DJ,
    6. Dhillon BK,
    7. Klintworth HM,
    8. Beauchamp M,
    9. Yeh MM,
    10. Elfers CT,
    11. et al.
    (2014) Reduced adiponectin signaling due to weight gain results in nonalcoholic steatohepatitis through impaired mitochondrial biogenesis. Hepatology 60:133–145.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Havilio M,
    2. Haddad Y, and
    3. Smilansky Z
    (2003) Intensity-based statistical scorer for tandem mass spectrometry. Anal Chem 75:435–444.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Keller A,
    2. Nesvizhskii AI,
    3. Kolker E, and
    4. Aebersold R
    (2002) Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal Chem 74:5383–5392.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Kim SC,
    2. Kim CK,
    3. Axe D,
    4. Cook A,
    5. Lee M,
    6. Li T,
    7. Smallwood N,
    8. Chiang JY,
    9. Hardwick JP,
    10. Moore DD,
    11. et al.
    (2014) All-trans-retinoic acid ameliorates hepatic steatosis in mice by a novel transcriptional cascade. Hepatology 59:1750–1760.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Koh KH,
    2. Pan X,
    3. Shen HW,
    4. Arnold SL,
    5. Yu AM,
    6. Gonzalez FJ,
    7. Isoherranen N, and
    8. Jeong H
    (2014) Altered expression of small heterodimer partner governs cytochrome P450 (CYP) 2D6 induction during pregnancy in CYP2D6-humanized mice. J Biol Chem 289:3105–3113.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Lawrence JA,
    2. Adamson PC,
    3. Caruso R,
    4. Chow C,
    5. Kleiner D,
    6. Murphy RF,
    7. Venzon DJ,
    8. Shovlin M,
    9. Noone M,
    10. Merino M,
    11. et al.
    (2001) Phase I clinical trial of alitretinoin and tamoxifen in breast cancer patients: toxicity, pharmacokinetic, and biomarker evaluations. J Clin Oncol 19:2754–2763.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Lee YS and
    2. Jeong WI
    (2012) Retinoic acids and hepatic stellate cells in liver disease. J Gastroenterol Hepatol 27 (Suppl 2):75–79.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Li Y,
    2. Wong K,
    3. Walsh K,
    4. Gao B, and
    5. Zang M
    (2013) Retinoic acid receptor β stimulates hepatic induction of fibroblast growth factor 21 to promote fatty acid oxidation and control whole-body energy homeostasis in mice. J Biol Chem 288:10490–10504.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Lilley JS,
    2. Linton MF, and
    3. Fazio S
    (2013) Oral retinoids and plasma lipids. Dermatol Ther (Heidelb) 26:404–410.
    OpenUrlCrossRef
  23. ↵
    1. Liu H,
    2. Sadygov RG, and
    3. Yates JR 3rd.
    (2004) A model for random sampling and estimation of relative protein abundance in shotgun proteomics. Anal Chem 76:4193–4201.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Liu Y,
    2. Chen H,
    3. Wang J,
    4. Zhou W,
    5. Sun R, and
    6. Xia M
    (2015) Association of serum retinoic acid with hepatic steatosis and liver injury in nonalcoholic fatty liver disease. Am J Clin Nutr 102:130–137.
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Mamoon A,
    2. Subauste A,
    3. Subauste MC, and
    4. Subauste J
    (2014) Retinoic acid regulates several genes in bile acid and lipid metabolism via upregulation of small heterodimer partner in hepatocytes. Gene 550:165–170.
    OpenUrlCrossRefPubMed
  26. ↵
    1. McClintick JN,
    2. Crabb DW,
    3. Tian H,
    4. Pinaire J,
    5. Smith JR,
    6. Jerome RE, and
    7. Edenberg HJ
    (2006) Global effects of vitamin A deficiency on gene expression in rat liver: evidence for hypoandrogenism. J Nutr Biochem 17:345–355.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Moriwaki H,
    2. Blaner WS,
    3. Piantedosi R, and
    4. Goodman DS
    (1988) Effects of dietary retinoid and triglyceride on the lipid composition of rat liver stellate cells and stellate cell lipid droplets. J Lipid Res 29:1523–1534.
    OpenUrlAbstract
  28. ↵
    1. Napoli JL
    (1996) Biochemical pathways of retinoid transport, metabolism, and signal transduction. Clin Immunol Immunopathol 80:S52–S62.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Nesvizhskii AI,
    2. Keller A,
    3. Kolker E, and
    4. Aebersold R
    (2003) A statistical model for identifying proteins by tandem mass spectrometry. Anal Chem 75:4646–4658.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Obrochta KM,
    2. Krois CR,
    3. Campos B, and
    4. Napoli JL
    (2015) Insulin regulates retinol dehydrogenase expression and all-trans-retinoic acid biosynthesis through FoxO1. J Biol Chem 290:7259–7268.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. O’Byrne SM and
    2. Blaner WS
    (2013) Retinol and retinyl esters: biochemistry and physiology. J Lipid Res 54:1731–1743.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Oliveros LB,
    2. Domeniconi MA,
    3. Vega VA,
    4. Gatica LV,
    5. Brigada AM, and
    6. Gimenez MS
    (2007) Vitamin A deficiency modifies lipid metabolism in rat liver. Br J Nutr 97:263–272.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Paik J,
    2. Haenisch M,
    3. Muller CH,
    4. Goldstein AS,
    5. Arnold S,
    6. Isoherranen N,
    7. Brabb T,
    8. Treuting PM, and
    9. Amory JK
    (2014) Inhibition of retinoic acid biosynthesis by the bisdichloroacetyldiamine WIN 18,446 markedly suppresses spermatogenesis and alters retinoid metabolism in mice. J Biol Chem 289:15104–15117.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Panchaud A,
    2. Jung S,
    3. Shaffer SA,
    4. Aitchison JD, and
    5. Goodlett DR
    (2011) Faster, quantitative, and accurate precursor acquisition independent from ion count. Anal Chem 83:2250–2257.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Panchaud A,
    2. Scherl A,
    3. Shaffer SA,
    4. von Haller PD,
    5. Kulasekara HD,
    6. Miller SI, and
    7. Goodlett DR
    (2009) Precursor acquisition independent from ion count: how to dive deeper into the proteomics ocean. Anal Chem 81:6481–6488.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Pessayre D
    (2007) Role of mitochondria in non-alcoholic fatty liver disease. J Gastroenterol Hepatol 22 (Suppl 1):S20–S27.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Ryu S,
    2. Gallis B,
    3. Goo YA,
    4. Shaffer SA,
    5. Radulovic D, and
    6. Goodlett DR
    (2008) Comparison of a label-free quantitative proteomic method based on peptide ion current area to the isotope coded affinity tag method. Cancer Inform 6:243–255.
    OpenUrlPubMed
  38. ↵
    1. Schug TT,
    2. Berry DC,
    3. Shaw NS,
    4. Travis SN, and
    5. Noy N
    (2007) Opposing effects of retinoic acid on cell growth result from alternate activation of two different nuclear receptors. Cell 129:723–733.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Shaw N,
    2. Elholm M, and
    3. Noy N
    (2003) Retinoic acid is a high affinity selective ligand for the peroxisome proliferator-activated receptor beta/delta. J Biol Chem 278:41589–41592.
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Sookoian S,
    2. Rosselli MS,
    3. Gemma C,
    4. Burgueño AL,
    5. Fernández Gianotti T,
    6. Castaño GO, and
    7. Pirola CJ
    (2010) Epigenetic regulation of insulin resistance in nonalcoholic fatty liver disease: impact of liver methylation of the peroxisome proliferator-activated receptor γ coactivator 1α promoter. Hepatology 52:1992–2000.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Staels B
    (2001) Regulation of lipid and lipoprotein metabolism by retinoids. J Am Acad Dermatol 45:S158–S167.
    OpenUrlCrossRefPubMed
  42. ↵
    1. Tay S,
    2. Dickmann L,
    3. Dixit V, and
    4. Isoherranen N
    (2010) A comparison of the roles of peroxisome proliferator-activated receptor and retinoic acid receptor on CYP26 regulation. Mol Pharmacol 77:218–227.
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Topletz AR,
    2. Tripathy S,
    3. Foti RS,
    4. Shimshoni JA,
    5. Nelson WL, and
    6. Isoherranen N
    (2015) Induction of CYP26A1 by metabolites of retinoic acid: evidence that CYP26A1 is an important enzyme in the elimination of active retinoids. Mol Pharmacol 87:430–441.
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Trasino SE,
    2. Tang XH,
    3. Jessurun J, and
    4. Gudas LJ
    (2015) Obesity Leads to Tissue, but not Serum Vitamin A Deficiency. Sci Rep 5:15893. doi: 10.1038/srep15893.
    OpenUrlCrossRefPubMed
  45. ↵
    1. Ward RJ and
    2. Peters TJ
    (1992) The antioxidant status of patients with either alcohol-induced liver damage or myopathy. Alcohol Alcohol 27:359–365.
    OpenUrlPubMed
  46. ↵
    1. Wei Y,
    2. Rector RS,
    3. Thyfault JP, and
    4. Ibdah JA
    (2008) Nonalcoholic fatty liver disease and mitochondrial dysfunction. World J Gastroenterol 14:193–199.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Wolf G
    (2010) Retinoic acid activation of peroxisome proliferation-activated receptor delta represses obesity and insulin resistance. Nutr Rev 68:67–70.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Yanagitani A,
    2. Yamada S,
    3. Yasui S,
    4. Shimomura T,
    5. Murai R,
    6. Murawaki Y,
    7. Hashiguchi K,
    8. Kanbe T,
    9. Saeki T,
    10. Ichiba M,
    11. et al.
    (2004) Retinoic acid receptor alpha dominant negative form causes steatohepatitis and liver tumors in transgenic mice. Hepatology 40:366–375.
    OpenUrlCrossRefPubMed
  49. ↵
    1. Yang F,
    2. He Y,
    3. Liu HX,
    4. Tsuei J,
    5. Jiang X,
    6. Yang L,
    7. Wang ZT, and
    8. Wan YJ
    (2014) All-trans retinoic acid regulates hepatic bile acid homeostasis. Biochem Pharmacol 91:483–489.
    OpenUrlCrossRefPubMed
  50. ↵
    1. Yu S,
    2. Levi L,
    3. Siegel R, and
    4. Noy N
    (2012) Retinoic acid induces neurogenesis by activating both retinoic acid receptors (RARs) and peroxisome proliferator-activated receptor β/δ (PPARβ/δ). J Biol Chem 287:42195–42205.
    OpenUrlAbstract/FREE Full Text
Back to top

In this issue

Download PDF
Article Alerts
Email Article
Citation Tools

Research ArticleArticle

Retinoic Acid Increases Fatty Acid β-Oxidation

Sasmita Tripathy, John D Chapman, Chang Y Han, Cathryn A Hogarth, Samuel L.M. Arnold, Jennifer Onken, Travis Kent, David R Goodlett and Nina Isoherranen
Molecular Pharmacology May 1, 2016, 89 (5) 560-574; DOI: https://doi.org/10.1124/mol.116.103697
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo

Jump to section

Advertisement