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.
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.
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.
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.
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).
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).
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.
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.
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.
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].
↵This 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
- Copyright © 2016 by The American Society for Pharmacology and Experimental Therapeutics