Lipid droplet (LD) catabolism in hepatocytes is mediated by a combination of lipolysis and a selective autophagic mechanism called lipophagy, but the relative contributions of these seemingly distinct pathways remain unclear. We find that inhibition of lipolysis, lipophagy, or both resulted in similar overall LD content but dramatic differences in LD morphology. Inhibition of the lipolysis enzyme adipose triglyceride lipase (ATGL) resulted in large cytoplasmic LDs, whereas lysosomal inhibition caused the accumulation of numerous small LDs within the cytoplasm and degradative acidic vesicles. Combined inhibition of ATGL and LAL resulted in large LDs, suggesting that lipolysis targets these LDs upstream of lipophagy. Consistent with this, ATGL was enriched in larger-sized LDs, whereas lipophagic vesicles were restricted to small LDs as revealed by immunofluorescence, electron microscopy, and Western blot of size-separated LDs. These findings provide new evidence indicating a synergistic relationship whereby lipolysis targets larger-sized LDs to produce both size-reduced and nascently synthesized small LDs that are amenable for lipophagic internalization.
Introduction
Lipid droplets (LDs) are ubiquitous fat-storage organelles that serve as readily accessible reservoirs of high-energy substrates used for β-oxidation within mitochondria. In the parenchymal cells of the liver (hepatocytes), the aberrant accumulation of LDs is the hallmark of steatosis, a key pathological feature of nonalcoholic fatty liver disease, obesity, and metabolic syndrome. This steatosis is viewed as an imbalance between the process of lipid storage and utilization. Thus, an understanding of the cellular machinery required to synthesize and catabolize these organelles is of great interest and an area of intense study. Currently, there are two central processes known to mediate the breakdown of triacylglycerol (TAG) stored within LDs for subsequent oxidation within mitochondria: cytosolic lipolysis and autophagy. In the process of lipolysis, cytosolic lipases including adipose triglyceride lipase (ATGL), hormone-sensitive lipase (HSL), and monoglyceride lipase (MGL) act sequentially to catalyze the liberation of the three fatty acid (FA) moieties comprising the parent TAG molecule (Vaughan et al., 1964; Jenkins et al., 2004; Villena et al., 2004; Zimmermann et al., 2004). The free FAs (FFAs) released by this lipolytic process are presumed to provide substrates for mitochondrial β-oxidation or act as potent signaling molecules for a variety of cellular processes; alternatively, these FAs can be reesterified back into TAG for storage (Kennedy and Lehninger, 1949; Edens et al., 1990; Ong et al., 2011; Khan et al., 2015).
In addition to the actions of the cytoplasmic lipases, it is now established that the catabolic process of autophagy can be used to mobilize LDs during periods of nutrient stress (Singh et al., 2009). Autophagy involves a highly orchestrated network of proteins that act in concert to selectively sequester intracellular contents within double-membrane structures known as autophagosomes. Fusion of autophagosomes with components of the terminal endocytic pathway (e.g., lysosomes) results in the recycling of autophagic cargo into macromolecular components within structures known as autolysosomes. In a highly selective form of LD-targeted autophagy, referred to as “lipophagy,” the specific turnover of LDs occurs through the action of acid lipases deposited into the autolysosome (Kaur and Debnath, 2015). Lipophagy thus represents an alternative to conventional cytosolic lipase-driven LD breakdown (Weidberg et al., 2009; Singh and Cuervo, 2012; Liu and Czaja, 2013; Schulze et al., 2017).
The relative utilization of lipolysis versus lipophagy by hepatocytes and other cells is presently unclear, as manipulation of either of these catabolic processes in mouse models can ultimately result in fatty liver (Singh et al., 2009; Ong et al., 2011). Whether lipolysis and lipophagy occur independently of each other or in tandem is an area of current investigation; indeed, an understanding of the crosstalk occurring between these pathways is only now beginning to emerge.
Evidence suggests that the size of the cargo targeted for degradation may be an important determinant in the capacity of the autophagic machinery to degrade entire organelles; this was recently demonstrated to be the case during a mitochondria-selective form of autophagy known as mitophagy (Gomes et al., 2011). In hepatocytes, LDs have diameters ranging from 60 nm to well over 5 µm in steatotic conditions. We therefore asked whether LD size might play an equally important role in dictating the order and/or prevalence of the catabolic processes used by cells for LD breakdown.
In this study, we find that ATGL, the rate-limiting cytoplasmic lipase, preferentially operates on the largest LDs within the hepatocyte, whereas the lipophagic machinery is restricted in its targeting to only the smallest populations of cytoplasmic LDs (i.e., only those with diameters of <1 µm). We therefore propose that LD size itself represents a fundamental physical parameter dictating the mechanistic processes used for cellular TAG catabolism. Inhibition of neutral lipase activity was shown to result in a significant increase in the average diameter of hepatic LDs. In contrast, the pharmacological or genetic inhibition of lysosomal acid lipase (LAL) led to a two- to threefold accumulation of LDs with diameters averaging <1 µm in size. Importantly, this accumulation of very small LDs by LAL inhibition was completely abrogated by joint inhibition of ATGL, implying the necessity of upstream cytosolic lipase function on large LDs. Differential centrifugation was used to isolate size-based subpopulations of LDs, revealing a preference for ATGL on large LDs and an enrichment of autophagic markers on small LDs. The existence of two independently regulated triglyceride disposal pathways presumably allows the cell to exert tight control over the process of lipid catabolism under conditions of nutritional stress.
Results
Differential LD morphologies by inhibition of lipolysis or lipophagy
As mentioned above, cargo size is believed to play a role in the selective autophagy of organelles such as mitochondria (Gomes et al., 2011). This premise is consistent with our initial findings from transmission EM of cultured primary rat hepatocytes showing two distinct size-based populations of LDs. As shown in Fig. 1 A, transmission EM performed on freshly isolated primary rat hepatocytes revealed a substantial population of large, cytosolic LDs that interact with adjacent cellular organelles such as mitochondria and ER. In addition to these large LDs, we observed numerous small LDs encased within electron-dense degradative vesicles characteristic of autolysosomes or late endosomes (Fig. 1 B and Fig. S5). Manual quantification of LD sizes from these electron micrographs (n = 4 cells) showed that cytosolic LDs range between 0.1 and 2.0 µm in diameter, with an average diameter of 0.7 µm. In contrast, lipophagic LDs exhibit a more restricted size range between 60 and 500 nm in diameter, with an average diameter of 220 nm (Fig. 1 C). Since the majority of these lipophagic LDs are considerably larger than extracellular lipoproteins that range between 7 and 50 nm in diameter (Rye et al., 1999), we reasoned that the LDs in Fig. 1 B were likely internalized as cytoplasmic LDs. These data indicate that the lipophagy machinery is restricted to degrading small LDs, and suggest that large LDs are perhaps preferentially degraded by lipolysis via ATGL, the rate-limiting enzyme in this process (Smirnova et al., 2006).
To test this concept, we measured LD morphology in AML12 mouse hepatocytes treated for 24 h with the ATGL inhibitor atglistatin (ATGLi; 20 µM) or the lysosome inhibitor chloroquine (100 µM). We predicted that lipolysis inhibition would result in cells containing large LDs, whereas lysosome inhibition would cause the accumulation of many small LDs trapped within autophagic vesicles. Indeed, as displayed in the confocal micrographs in Fig. 1 D, treatment with these pharmacological inhibitors resulted in dramatic differences in LD morphology as cells treated with the ATGL inhibitor accumulated very large LDs that were nearly twofold larger than DMSO controls, whereas lysosomal perturbation by chloroquine led to the accumulation of numerous small LDs that were roughly half the size of control LDs (Fig. 1 E). In addition, chloroquine treatment caused a nearly fourfold increase in LD number per cell compared with DMSO controls (Fig. 1 F), whereas the total LD area per cell was not significantly increased after 24 h (Fig. 1 G). These observations suggest that ATGL-driven lipolysis preferentially targets larger-sized LDs, whereas small LDs are preferentially catabolized by acidic lipases within lysosomes and/or late endosomes.
Evidence for sequential lipolysis and lipophagy pathways
Based on the distinct changes in LD morphology following ATGL versus lysosome inhibition, we tested whether lipolysis and lipophagy operate independently, or if these processes operate in tandem. To test this, we assessed LD morphology in AML12 hepatocytes treated with siRNA against ATGL, LAL, or both. Validation of ATGL and LAL knockdown levels by Western blot and quantitative PCR (qPCR) are shown in Fig. S1 A. As expected, ATGL knockdown resulted in a twofold increase in average LD size over nontargeting siRNA controls (Fig. 2, A and B), and LAL knockdown caused a substantial accumulation of smaller-sized LDs (Fig. 2, A and C), similar to our previous observations using ATGLi and chloroquine in Fig. 1. However, dual knockdown of both ATGL and LAL caused the accumulation of large LDs similar to knockdown of ATGL alone, as LD size was similarly increased about twofold over control siRNA, and the accumulation of small LDs by LAL knockdown was completely reversed (Fig. 2, A–C). Similar increases in total LD area per cell were observed following knockdown of ATGL, LAL, or both (Fig. 2 D), suggesting that these two processes are not independent but operate in tandem.
As seen in confocal micrographs in Fig. 2 E, similar results were also observed using the ATGL inhibitor ATGLi (20 µM), the LAL inhibitor LAListat (LALi; 50 µM) or both in AML12 cells loaded with oleic acid (OA; 150 µM, 2 h) and chased for 48 h in insulin-free medium with reduced serum (2% FBS) to stimulate LD breakdown. Similar to siRNA knockdown experiments, a twofold increase in LD size was observed in in ATGLi- or ATGLi+LALi-treated cells (Fig. 2 F), and inhibition of LAL caused a more than twofold increase in LD number per cell, which was largely reversed by dual inhibition of ATGL+LAL (Fig. 2 G). These effects were specific for ATGL activity, as inhibition of other cytosolic lipases (HSL and MGL) using CAY10499 increased both large and small LDs and did not reverse the LALi-induced accumulation of small LDs (Fig. S2). However, consistent with siRNA knockdown of ATGL and LAL, the increase in overall LD content was similar with ATGLi, LALi, or both as measured by confocal microscopy (Fig. 2 H) and biochemical TAG analysis (Fig. S1 B). Together, these results suggest that the accumulation of small LDs following lipophagy inhibition is dependent on upstream ATGL-driven lipolysis.
To extend these observations made in the AML12 cells, we also tested the impact of lipolysis and lipophagy inhibition using freshly isolated primary rat hepatocytes treated for 24 h with 20 µM ATGLi, 50 µM LALi, or both. As seen in confocal images and quantification in Fig. 3, the inhibition of ATGL dramatically increased LD size over control hepatocytes (Fig. 3, A, B, and E), whereas LAL inhibition increased the LD number (Fig. 3, C and F). Consistent with the AML12 cells, dual inhibition of both ATGL and LAL resulted in large LDs similar to ATGLi alone (Fig. 3, D and G).
To assess the effects of lipolysis and lipophagy inhibition under lipolysis stimulation, primary rat hepatocytes were also treated with forskolin (10 µM, 24 h), a cAMP agonist that activates lipolysis (Schott et al., 2017). As seen in confocal images and corresponding graphs from Fig. 3, A′–D′, inhibition of ATGL, LAL, or both resulted in similar effects on LD morphology between stimulated and unstimulated hepatocytes. When compared with unstimulated cells, forskolin treatment caused a significant decrease in total LD area/cell that was largely due to a reduction in larger-sized LDs (Fig. 3 E, gray bars), whereas no significant decrease in smaller-sized LDs was observed (black bars). Large versus small LDs were defined as greater than or less than 1 µm2, respectively, as this accounted for roughly 50% of the total LD area in each of the two categories in the control (unstimulated) condition. The decrease in LD area per cell by forskolin was completely prevented by ATGL inhibition or by dual inhibition of ATGL+LAL. Taken together, these data suggest that ATGL-driven lipolysis preferentially targets the degradation of large LDs upstream of LAL-driven lipophagy.
To further define the upstream role of lipolysis relative to lipophagy, we analyzed the kinetics of ATGL versus LAL inhibition in both primary rat hepatocytes and AML12 cells. We predicted that ATGL inhibition would show effects on LD morphology and overall LD content at an earlier time point than LAL inhibition. As seen in Fig. 4 A, ATGL inhibition increased LD size in as little as 3 h in primary rat hepatocytes, while LD number/cell was not affected by ATGL or LAL inhibition during this time (Fig. 4 B). LALi caused an unexpected 20% increase in LD size at t = 6 h (Fig. 4 A), indicating alternative effects of the inhibitor perhaps at acute time points. Nonetheless, LD size was substantially increased 70% by the ATGL inhibitor and was statistically significant within 3 h of treatment. In addition, ATGL inhibition caused a significant increase in total LD area/cell within 6 h over DMSO control (Fig. 4 C), whereas no effect was observed by LAL during this time. As seen in Fig. 4, D–F, similar effects were seen in AML12 cells that were loaded with OA and then washed and chased in OA-free medium. First, ATGL inhibition caused a twofold increase in LD size within 6 h of treatment (Fig. 4 D), while maintaining the total LD area within 10 h (Fig. 4 F). These findings show that ATGL inhibition affects LD size and total area at earlier time points than the LAL inhibitor, which is consistent with the concept that the lipolysis of large LDs precedes the downstream lipophagy of smaller-sized LDs.
Lipolysis and lipophagy machinery preferentially target size-based LD subpopulations
Based on the findings described above that suggest tandem, size-based lipolysis and lipophagy pathways, we tested if components of these processes might preferentially associate with large and small LDs, respectively, in OA-loaded AML12 cells (150 µM, 16 h). As seen in Fig. 5 A, AML12 hepatocytes display numerous examples of both small and large LDs decorated by endogenous ATGL by immunostaining. However, a subset of small, ATGL-negative LDs was commonly observed, while ATGL fluorescence appeared most prevalent around the largest cellular LDs. Quantification of LD size (Fig. 5 A′) revealed that ATGL-positive LDs were roughly sevenfold larger than ATGL-negative LDs (4.00 versus 0.61 µm2, respectively), suggesting that ATGL preferentially, but not exclusively, localizes to larger-sized LDs. Interestingly, in OA-loaded AML12 cells that were subsequently washed and treated with the lysosome inhibitor chloroquine (50 µM, 24 h), ATGL localization was readily observed on smaller-sized LDs overall (Fig. S3 A, note the numerous small green ATGL puncta that colocalize with Oil Red O–stained LDs) but was nonetheless absent on a subset of small LDs, as indicated by the white arrows in Fig. S3 A. Quantification of LD size revealed that in chloroquine-treated cells, ATGL-negative LDs were 60–70% smaller than ATGL-positive LDs (Fig. S3 B).
From these findings, it became relevant to test if lipophagic vesicles such as autophagosomes and endo-lysosomes exhibit a preferential association with small versus large LDs. As seen in the confocal micrographs of OA-loaded AML12 cells in Fig. 5 B, small LDs were observed within GFP-LC3 autophagosomes that were 70% smaller than LC3-negative LDs. Similar targeting of small LDs by LC3 was also observed in OA-loaded AML12 cells treated with chloroquine (Fig. S3, B and C). Notably, we found few examples of LDs that were totally engulfed within LC3-positive autophagosomes, which is consistent with previous observations (Ding et al., 2010). However, the internalization of cytosolic LDs was substantially more prevalent within endo-lysosomal vesicles marked by the tetraspanin CD63-EGFP. Numerous examples of small LDs were found within these vesicles in OA-loaded cells both with and without chloroquine treatment. Like LC3, LD size quantification revealed that LDs within CD63 vesicles were ∼70% smaller than LDs that were not associated with these vesicles (Figs. 5 C and S3 C). Taken together, these data suggest that lipolysis and lipophagy machinery exhibit preferential targeting to large and small LDs, respectively.
To further define the terminal degradative compartment of small LDs, pulse-chase lipid analysis was performed using fluorescent BODIPY C12 FA in AML12 cells treated 24 h with DMSO or LALi (50 µM). We predicted that cells treated with lysosome inhibitors might accumulate substantial numbers of small LDs within the acidic vesicular compartments compared with control cells. As seen in Fig. 5 D, white arrows depict little or no colocalization between green fluorescent BODIPY C12 and LysoTracker in OA-loaded AML12 cells (16 h, 150 µM OA + 1 µM BODIPY C12). However, LAL inhibition (50 µM, 24 h) following OA load caused a substantial overlap between BODIPY fluorescence and LysoTracker vesicles, indicating the accumulation of small cytosolic LDs within autolysosomal or late endosomal vesicles. To further define the ultrastructure of these vesicles, transmission electron micrographs were taken of OA-loaded AML12 cells treated with or without LALi (50 µM, 24 h). As seen in Fig. 5 E, OA-loaded control cells contained large, cytosolic LDs in contact with other organelles such as mitochondria, whereas LALi-treated cells accumulated numerous cytoplasmic LDs and a subset of small LDs ∼200 nm in diameter within degradative vesicles resembling autolysosomes (yellow arrows).
To determine the relative proportion of small LDs that are incorporated into lipophagic vesicles, AML12 cells were immunostained for endogenous LAMP1, a lysosomal/late endosomal marker, following pulse-chase analysis using fluorescent BODIPY C12. In these experiments, cells were pulse labeled for 2 h with 150 μM OA + 7.5 μM BODIPY C12 558/568, then chased for 24 h in regular growth medium (10% FBS; ITS) containing DMSO or LALi (50 µM). As seen in Fig. 5 F, the accumulation of small LDs within LAMP1-positive vesicles was greatly increased by LAL inhibition. Quantification revealed a 13-fold increase in the number of LAMP1-positive LDs/cell (70–80 cells per condition, n = 3 experiments; Fig. 5 G). Strikingly, as much as 40% of the LDs in each cell were internalized by LAMP1-positive vesicles following LAL treatment (Fig. 5 H) even as the number of LDs increased 2.4-fold during this time (Fig. 5 I). These data suggest that a substantial number of small LDs are terminally degraded by lipophagy within acidic organelles such as lysosomes, autolysosomes, and late endosomes.
Large and small LDs possess unique protein and lipid signatures
As an extension to the morphological-based studies described above, biochemical criteria were applied to isolated LDs of different sizes. To generate purified LD fractions enriched in small versus large LDs, a two-step differential centrifugation technique was used as depicted in Fig. 6 A (adapted from Brasaemle and Wolins, 2016; Zhang et al., 2016). This method promoted a size separation of large, buoyant LDs in an “LD1” fraction and smaller, less buoyant LDs within an “LD2” fraction that were then assessed by both immunofluorescence and Western blot analyses. As seen in Fig. 6 B, both fractions contained LDs that were positive for Perilipin 2 (PLIN2) immunostaining as well as Nile Red and monodansylpentane (MDH) lipid stains, and the LD1 fraction contained LDs that were substantially larger than that of LD2. Western blot analysis indicated that these LD fractions were not enriched in organelle contaminants such as Golgi, ER, and mitochondria (Fig. S4 A). Most interesting are the findings indicating a marked enrichment of ATGL in the LD1 fraction containing large LDs (more than twofold) while autophagic (LC3, 2.5-fold) and lysosomal markers (LAMP2A, 10-fold) were substantially enriched in the smaller LD2 fraction (Fig. 6, C–F). Confocal microscopy of LD1 and LD2 fractions supported these findings, as ATGL immunofluorescence was abundant on large LDs, whereas LysoTracker staining colocalized with small LDs in the LD2 fraction (Fig. S4 B). These findings are consistent with a size-based preference for components of the lipolysis and lipophagy machineries.
As large and small LDs appear to possess distinct protein signatures (Fig. 6, A–G), we sought to assess whether distinct lipid species were contained within these two subpopulations. Based on the concept that ATGL acts to reduce the size of larger LDs, we predicted that smaller-sized LDs would accumulate remnant lipid species resulting from triglyceride (TAG) lipolysis such as diacylglycerol (DAG) or cholesterol esters (CEs) for which ATGL has lower affinity. Prior to lipidomic analyses, LD1 and LD2 fractions were isolated from whole rat livers (n = 4) using a slightly modified differential centrifugation technique than that of the AML12 cells (described in the Materials and methods section), and confocal microscopy confirmed the size separation of LDs within these fractions (Fig. S4 C). Additional Western blot analysis confirmed the enrichment of LC3B and LAMP2A in the small LD2 fraction, and ATGL enrichment in the large LD1 fraction (Fig. 6, G–J). Surprisingly, lipidomic analysis of large and small LDs revealed near equivalent levels of triglycerides and CEs as a percentage of total lipid abundance (Fig. 6 K and Fig. S4, D and E). Fig. 6 K depicts a heatmap of the 25 most abundant lipid species representing TAGs and CEs that exhibited a similar abundance between the two LD subpopulations. However, as seen in Fig. 6 L, several lipid species were significantly enriched in the LD2 fraction, most notably several DAGs (approximately twofold) and phospholipids phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylinositol (PI; approximately threefold). No lipid classes or individual lipid species were markedly enriched in the LD1 fraction relative to LD2. The enrichment in PC, PE, and PI in LD2 indicates the presence of smaller-sized LDs with higher surface-to-volume ratios, whereas the enrichment in DAG is likely the result of upstream ATGL-driven lipolysis of larger-sized LDs.
Synthesis-dependent and -independent accumulation of small LDs
The observations depicted in Figs. 1, 2, 3, 4, 5, and 6 are consistent with the premise of a tandem stepwise pathway whereby ATGL-driven lipolysis acts on large LDs (>1 µm2) to create a population of small LDs (<1 µm2) that are more amenable for LAL-driven lipophagy. A previous report using cultured adipocytes showed that lipolysis activation forms a nascent population of smaller LDs generated by reesterification of FFAs (Paar et al., 2012). In addition, the contribution of extracellular lipoproteins internalized by clathrin-mediated endocytosis may also provide an alternative source for the generation of small LD accumulation following LAL and lysosome inhibition. To further measure the importance of LD synthesis or exogenous lipids in lipolysis and lipophagy synergy, cells were incubated in the presence of DGAT1 (10 µM PF-06424439) and DGAT2 (10 µM PF-04620110) inhibitors in regular medium and lipoprotein-depleted medium (LDM). Thus, by removing either exogenous lipids or by inhibiting nascent LD synthesis, we tested the relative contributions of these pathways in the accumulation of LDs by lipolysis and lipophagy machinery.
Accordingly, freshly isolated primary hepatocytes were treated for 24 h with DMSO, 20 µM ATGLi, 50 µM LALi, or both (ATGLi+LALi) in control medium (CM), LDM, or CM containing DGAT1/2 inhibitors. First, similar changes in LD size were observed following ATGLi or ATGLi+LALi regardless of exogenous lipids or LD synthesis (Fig. 7 A). In CM, LD number per cell was increased by LALi treatment as expected, and this increase was not observed in the combined presence of ATGLi+LALi (Fig. 7 B), confirming that the increase in LD number caused by lysosome inhibition is first dependent on ATGL-driven lipolysis. In LDM, LALi also caused an increase in LD number per cell, suggesting that exogenous lipids are not required for this increase (Fig. 7 B, gray bars). In the presence of DGAT1/2i, overall LD numbers were reduced with no significant increase by LALi treatment, suggesting that FFA reesterification plays a substantial role in the lipolysis-dependent accumulation of LDs (Fig. 7 B). As seen in Fig. 7 C, similar effects are also observed in total LD area per cell.
The above experiments were also conducted in AML12 cells using siRNA knockdown instead of pharmacological inhibitors. These knockdown variables included nontargeting siRNA, siATGL, siLAL, or siATGL+siLAL. siRNA treatment resulted in a >80% knockdown for both ATGL and LAL as assessed by mRNA and protein levels (Fig. S1 A). Following a 48-h siRNA knockdown, cells were replated for an additional 16 h in CM, LDM, or CM containing DGAT1/2 inhibitors. As shown in Fig. 7 D, LD size was increased in response to ATGL or ATGL+LAL knockdown in all three conditions, similar to primary rat hepatocytes, suggesting that exogenous lipids and FFA reesterification do not impact LD size–based effects of ATGL inhibition. Importantly, the number of LDs per cell was higher in LAL knockdown cells regardless of the culture conditions used (CM, LDM, or DGAT1/2i) as shown in Fig. 7 E, while the overall LD numbers were reduced by roughly 80% when in the absence of exogenous lipids or LD synthesis. Like the primary hepatocytes, exogenous lipids were not required to maintain higher LD numbers or total LD area in LAL-depleted cells (Fig. 7, E and F, gray bars). However, the finding that LD numbers and total LD area were higher in LAL knockdown cells treated with DGAT1/2 inhibitors (Fig. 7, E and F) suggests that lipoproteins and size-reduced cytosolic LDs also contribute to lipid accumulation under these conditions. Taken together, our data suggest that the LD synthesis machinery drives the accumulation of most small LDs that appear following lysosome inhibition, and this requires FAs generated upstream by ATGL-driven lipolysis.
Discussion
The central findings of this study provide evidence that the hepatocyte is likely to catabolize cytoplasmic LDs using both cytoplasmic lipases and lipophagy. Importantly, this utilization may reflect a stepwise process of small LD production by cytoplasmic lipases such as ATGL followed by autophagic engulfment as depicted in the cartoon working model from Fig. 8. This premise is based on the findings that dual treatment of hepatocytes with inhibitors of ATGL and LAL (Figs. 2, 3, and 7) leads to the accumulation of large LDs, and that lipolysis and lipophagy machinery exhibit size-based targeting of large and small LDs, respectively (Figs. 5 and 6). The LD morphologies resulting from inhibition of lipolysis and lipophagy are consistent with previous studies that have used lipase-defective mouse models that parallel human disease phenotypes. For example, liver-specific ATGL knockout mice exhibit very large hepatocellular LDs by immunohistochemical examination (Ong et al., 2011), and gross macrovesicular steatosis has been reported within periportal zones (Wu et al., 2011). Reciprocally, the accumulation of small LDs following lysosome/LAL perturbation in our cellular models is consistent with microvesicular steatosis in the hepatocytes of patients with LAL deficiencies and other types of lysosomal storage disorders (Bernstein et al., 2013; Reiner et al., 2014). Thus, the combined findings of all these studies support a concept of an epistatic cascade of cytosolic lipase-to-lysosome transition via preferential targeting of size-based LD subpopulations.
Utilization of lipolysis and lipophagy pathways in the hepatocyte
A primary question addressed by this study is defining the source of the many smaller LDs that accumulate in hepatocytes treated with LAL inhibitors (Figs. 1, 2, 3, and 7) and whether they are generated directly from large LDs that are reduced in size, from endocytosed lipids, or from synthesis of nascent LDs. A careful study by others had previously demonstrated that the production of small LDs in adipocytes under lipolytic stimulation can be curtailed by preventing FA esterification and synthesis of nascent LDs (Paar et al., 2012). This suggested that most of the small LDs generated in stimulated cells were not a result of fragmentation or shrinking preexisting LDs but of new synthesis from liberated FAs. While our study also showed that the accumulation of small LDs following LAL inhibition was lipolysis dependent (Figs. 2, 3, and 7), we also tested for alternative lipid trafficking pathways such as lipoprotein endocytosis and FA reesterification.
First, exogenous lipids were removed from the culture medium using lipoprotein-depleted serum to assess the role of endocytosed lipoproteins (Fig. 7). We observe that exogenous lipids contribute to overall LD content but are not required to observe an increase in LD number following LAL inhibition in primary hepatocytes or LAL knockdown in AML12 cells, suggesting that the contribution of endocytosed lipids may be modest. Consistent with this prediction are the many electron micrographs we have viewed, some of which are displayed in Fig. 1 and Fig. S5, that show late endocytic compartments containing large LDs (≥200 nm) that exceed the known diameters of lipoprotein particles ranging between 7 and 50 nm in human blood plasma (Rye et al., 1999). Chylomicron remnants within blood plasma have a somewhat larger diameters (50–150 nm) but are further reduced by hepatic lipase at the hepatocyte surface before endocytic uptake (Cooper, 1997; Crawford and Borensztajn, 1999; Perret et al., 2002). Thus, it is likely the majority of the lipophagic LDs we observe, particularly those ≥200 nm in diameter, are derived from the engulfment of small cytoplasmic LDs.
As an extension of these studies, we sought to attenuate synthesis of LDs rather than limiting the external lipid source. This was done using pharmacological inhibitors of DGAT1 and DGAT2 (Fig. 7), blocking the esterification of FAs into triglyceride, a required step in nascent LD synthesis. From these experiments, we found that most of the small LDs that accumulated following LAL inhibition or siRNA knockdown were dependent upon the action of the LD synthesis machinery (Fig. 7, B and E). Thus, it appears that hepatocytes use the lipophagic pathway to not only complete the catabolism of residual LD “remnants” following the action of cytoplasmic lipases such as ATGL, but also to target small newly synthesized LDs resulting from FA reesterification. While the protein blotting of the light and heavy LD fractions are consistent with the premise of ATGL- versus lipophagy-associated subpopulations (Fig. 6), this analysis does not preclude that the LD2 fraction contains a mixture of both nascent and remnant LDs. The extensive lipidomics of these fractions described in Fig. 6 and Fig. S3 does suggest that some conversion of TAG to DAG has taken place preferentially in the heavy LD2 fraction, although this contribution appears minor given the relative homogeneity of lipid classes between the two fractions.
A size-dependent selection of LDs by two central catabolic pathways
An interesting facet of the interplay between cytoplasmic and lysosomal lipases reported here is the concept that LD girth may influence which of these lipases acts first. As discussed above, we predict that the terminal autophagy of LDs is initiated by a size-based mechanism that is dependent on LDs reaching a certain modest diameter, not unlike that observed for mitophagy, as reported elsewhere (Gomes et al., 2011). This prediction is based on the substantial comparative differences observed in LD size phenotypes following the inhibition of either lipolysis or lipophagy. Consistent with this premise, the sizes of LDs were increased following the inhibition of ATGL, suggesting that lipolysis preferentially targets the degradation of large LDs while lipophagic membranes are perhaps unable to enclose such large cargo. Furthermore, primary rat hepatocytes stimulated for lipolysis show a reduction in LD content that appears to be due largely to a reduction in large LDs (Fig. 3 G). Finally, ATGL was observed to be markedly enriched on large LDs both by immunofluorescence microscopy (Fig. 5) and Western blot analysis of size-separated LDs from both AML12 cells and rat livers (Fig. 6). Some ATGL was observed on smaller LDs, however, particularly following lysosomal inhibition (Fig. S3).
The mechanisms that favorably target ATGL to large LDs are yet to be elucidated. It is known that ATGL localizes to LDs through a C-terminal hydrophobic domain (Murugesan et al., 2013) and is thought to be regulated both positively and negatively by phosphorylation at various sites (Mason et al., 2012; Pagnon et al., 2012; Xie et al., 2014). ATGL localization on LDs is also influenced by the ER-Golgi trafficking machinery (Soni et al., 2009), as well as direct binding to the autophagosome marker LC3 (Martinez-Lopez et al., 2016). While it is possible that ATGL phosphorylation or protein interaction dictates its localization to specific LDs, the mechanisms that regulate LD-protein interaction in general are poorly understood. For example, the LD monolayer itself has unique biophysical properties that serve as a gatekeeper for protein interaction. It has been reported that LD monolayers are under high surface tension in comparison to lipid bilayers, and this surface tension restricts many proteins that normally adhere to lipid membranes (Thiam et al., 2013; Prevost et al., 2018). This is likely exacerbated in differently sized LDs with varying degrees of membrane curvature and proteomic contents that could influence the binding of proteins such as ATGL to large versus small LDs. Given the wide heterogeneity of LDs even within individual cells (Hsieh et al., 2012; Zhang et al., 2016; Thiam and Beller, 2017), future studies will need to dissect the contributions of LD diversity, membrane biophysics, and protein–protein interactions in directing ATGL and other proteins to specific LD subpopulations.
The findings that lysosome inhibition increases LD number and incorporation into acidic degradative vesicles (Fig. 5) suggests that, even under basal conditions, the lipophagy pathway is exceptionally active in the catabolism of small LDs (<1 µm2) that may represent remnants of larger-sized LDs as well as a relevant population of newly synthesized LDs. The high basal activity of this process is somewhat surprising and highlights the importance of LD size recognition perhaps over the contributions of nutrient deprivation and other cellular cascades that activate LD–lipophagic interactions. A variety of reports have implicated nutrient deprivation as an activator of this autophagic pathway (reviewed in Schulze et al., 2017) suggesting that lipophagy is a highly regulated process that might work in collaboration with putative sensors of LD size or membrane curvature. Interestingly, LC3 lipidation was reported to occur preferentially on membranes with high curvature due to Atg3 sensing (Nath et al., 2014), raising the possibility that similar curvature-sensing mechanisms could exist for lipophagic targeting of small LDs. It is also likely that the size of degradative vesicles such as lysosomes, autolysosomes, and late endosomes are also limiting to LD internalization and lipid hydrolysis. For example, these vesicles may achieve diameters in excess of 1 µm under certain conditions, while the majority are much smaller in size (Huotari and Helenius, 2011; Xu and Ren, 2015; Su et al., 2016). This range in vesicle size is consistent with the relative size of lipophagic LDs (<1 µm2) that we observe by light and EM (Figs. 1, 5, and S5), and the size-based population of LDs that accumulate following LAL perturbation (Fig. 3 G).
Taken together, these findings support a sophisticated process of hepatocyte LD homeostasis that implicates the action of soluble lipases on large droplets that generate smaller-sized LDs by DGAT1/2 dependent reesterification of FFAs and, to a lesser extent, size-reduction of catabolized LDs. Smaller remnant LDs, generated from the action of ATGL lipolysis and nascent synthesis, attract components of the lipophagic machinery based on a reduced girth (<1 µm2) for degradation within acidic vesicles. Further insights into how the hepatocyte might alter this elaborate process under normal and disease states will prove important.
Lipolysis and lipophagy as tandem pathways for LD breakdown
The stepwise model suggesting size-based LD targeting by lipolysis and lipophagy pathways in the current study is supportive of previous studies showing that ATGL is an important upstream mediator of downstream autophagy. For example, ATGL is well described to regulate SIRT1 to induce autophagy and FA oxidation in hepatocytes via PPARα/PGC-1α signaling (Khan et al., 2015; Sathyanarayan et al., 2017). In agreement with this, we also show that ATGL operates upstream of lipophagy by preferentially targeting large LDs and producing numerous smaller-sized LDs via DGAT1/2 reesterification and, to a lesser extent, reducing the size of existing LDs. Other reports have indicated that autophagy can also act upstream of ATGL/lipolysis by promoting LD synthesis under nutrient-limited conditions (Rambold et al., 2015; Nguyen et al., 2017; Zechner et al., 2017). In this model, macroautophagy facilitates the bulk degradation of cellular organelles, and byproduct membrane lipids are recycled to the ER for LD biogenesis. It is likely that bulk macroautophagy is a separate process from more selective forms of autophagy such as LD-specific lipophagy. Thus, while our data support the hypothesis that ATGL acts upstream of lipophagy, bulk macroautophagy likely plays an important role in LD synthesis further upstream of both lipolysis and lipophagy.
In conclusion, the current study defines lipolysis and lipophagy as tandem pathways in hepatocytes that degrade large and small LDs, respectively. The enrichment of ATGL on larger-sized LDs is consistent with the increase in LD size following ATGL inhibition or siRNA depletion. Conversely, small LDs are preferentially engulfed by lipophagic machinery, and lysosomal inhibition results in a decrease in LD size and the accumulation of numerous small LDs within autophagic and endosomal vesicles. These LD phenotypes correlate with human pathologies of macro- and microvesicular steatosis, which may suggest the need for targeted therapeutic strategies in the treatment of specific fatty liver pathologies.
Materials and methods
Cell culture and reagents
Primary rat hepatocytes were isolated from female Sprague Dawley rats (Envigo) by collagenase perfusion and were cultured in William’s E Medium (5% FBS) as described by previously (Shen et al., 2012). All animals received humane care in accordance with the guidelines established by the American Association for the Accreditation of Laboratory Animal Care, and animal protocols were approved by the Institutional Animal Care and Use Committee at Mayo Clinic. The AML12 cells (CRL-2254; ATCC) are mouse hepatocytes obtained from transgenic MT42 mice (CD1 strain) expressing human transforming growth factor α. AML12 cells were maintained in 1:1 DMEM/F12 medium supplemented with 10% FBS, Pen/Strep, and Insulin-Transferrin-Selenium supplement (41400-05; Gibco; final concentration: 0.01 mg/ml insulin, 0.0055 mg/ml transferrin, and 6.7 ng/ml selenium). The β-actin antibody (A2066), chloroquine diphosphate salt (C6628), DGAT1 inhibitor PF-04620110 (PZ0207), DGAT2 inhibitor PF-06424439 (PZ0233), lipoprotein-depleted serum (S5394), OA (O1008), and Oil Red O (O0625) were from Sigma-Aldrich. The Nile Red (N1142) was from Thermo Fisher Scientific. The MDH was from Abgent (SM1000a). The ATGL antibody (2439) was from Cell Signaling Technology. The PLIN2 antibody (B3121) was from LS Biosciences. The LAMP2A antibody (ab18528) was from Abcam. The LC3 antibody (NB600-1384) was from Novus. The Tsg101 antibody (GTX70255) was from Genetex. The Lamp1 antibody (1D4B) was from the Developmental Studies Hybridoma Bank at the University of Iowa (Iowa City, IA). The ATGL inhibitor ATGLi (#15284) and HSL/MGL inhibitor CAY10499 (10007875) were from Cayman Chemical. The LAL inhibitor LALi (6098) was from Tocris. The C12-BODIPY (FL: D3822; 558/568: D3835) and LysoTracker Deep Red (L12492) were from Thermo Fisher Scientific. The human ATGL-EGFP construct was a kind gift from C. Jackson (Institut Jacques Monod, Paris, France). The EGFP-LC3 was a gift from K. Kirkegaard (Stanford, CA; Addgene plasmid 11546). The CD63-EGFP construct was generated from a human cDNA library (primers: 5′-AAGCTTATGGCGGTGGAAGGAGGA-3′, 5′-GAATTCGGAGACCAGACCCCATCACG-3′), ligated into EGFP N1 vector following restriction endonuclease digestion (5′-HindIII, 3′-EcoRI) and confirmed by sequence analysis. Cells were transfected at 60–80% confluency using Lipofectamine 2000 transfection reagent (Thermo Fisher Scientific) according to the manufacturer’s instructions.
Fluorescence and live-cell microscopy
Cells were washed in PBS and fixed in 3% formaldehyde as described previously (Cao et al., 1998). To label LDs, fixed samples were washed in 60% isopropanol for 30 s and 60% Oil Red O solution (5 mg/ml in isopropanol) for 1.5 min and then washed in 60% isopropanol for an additional 30 s. Images were acquired using a Zeiss LSM 780 confocal microscope with a 40× oil objective lens (NA = 1.4). Quantification of LD size, number, and total area per cell was done using ImageJ software (Auto Local Threshold tool, Bernsen method).
Triglyceride measurement
Triglycerides were measured from AML12 cells in accordance with a previous protocol (Wang et al., 2017). In brief, lipid was extracted using 1 ml chloroform/methanol (2:1) and vigorous shaking for 1 h at RT. 200 µl double-distilled H2O was added, and samples were vortexed and centrifuged at 3,000 × g for 5 min. The lower lipid phase was collected (∼700 µl), dried overnight, and resuspended in a mixture containing tert-butanol (471712; Sigma-Aldrich), Triton X-114 (X114; Sigma-Aldrich), and methanol (9:4:2). TAG quantification was done in AML12 cells using a colorimetric kit from Pointe Scientific (T7532, T7531-STD) according to manufacturer’s instructions.
LD isolation
AML12 cells were grown to 90% confluency in 5 × 15–cm dishes and loaded with OA (150 µM, 16 h). Using a protocol adapted from two previous studies (Brasaemle and Wolins, 2016; Zhang et al., 2016), cells were washed, scraped, and incubated in a hypotonic lysis medium (HLM) before Dounce homogenization. The post-nuclear supernatant was placed at the bottom of a 25% to 0% OptiPrep density gradient (D1556; Sigma-Aldrich) in HLM. Following a 30-min spin at 36,000 × g, the floating fat layer (LD1) was collected and washed for subsequent Western blot analysis. After LD1 was collected, 400 µl of HLM was carefully added to the top of the gradient, and samples were spun for an additional 60 min at 182,000 × g. The LD2 fraction was then collected from the top 200–400 µl. LD2 samples were loaded on SDS-PAGE gels at twice the volume of LD1 samples to account for differences in LD yield, and protein levels were normalized to PLIN2 as a loading control.
To isolate LDs from rat livers, ∼5 g of tissue was surgically extracted from female Sprague Dawley rats (Envigo). Liver tissue was finely minced and washed in ice-cold PBS before incubation in HLM, Dounce homogenization, and layering beneath the OptiPrep gradient as described above. To collect the LD1 fraction from rat whole liver lysates, samples were spun 30 min at 8,000 × g. Following the collection of the floating fat layer (LD1), the samples were spun for 60 min at 182,000 × g to collect the floating LD2 fraction.
siRNA knockdown and qPCR
Cells were transfected with siGENOME pools L/M-048418-01 and M-040220-01 (Dharmacon) targeting murine LIPA and PNPLA2 RNA using the RNAiMAX reagent (Invitrogen) according to manufacturer’s instructions. A nontargeting siRNA duplex from Dharmacon was used in control transfected cells. RNA was isolated from cells using the RNeasy Plus Mini Kit (Qiagen) and was reverse-transcribed with the Super Script III First Strand Kit (Invitrogen) using oligo-(dT) primers. Quantification of gene expression was performed using SYBR green fluorescence on a LightCycler 480 (Roche). Primers for murine LIPA (lysosomal acid lipase; forward: 5′-TGCCCACGGGAACTGTATC-3′ and reverse: 5′-ATCCCCAGCGCATGATTATCT-3′), PNPLA2 (forward: 5′-ATGTTCCCGAGGGAGACCAA-3′ and reverse: 5′-GAGGCTCCGTAGATGTGAGTG-3′), and 18S ribosomal RNA (forward: 5′-CGCTTCCTTACCTGGTTGAT-3′ and reverse: 5′-GAGCGACCAAAGGAACCATA-3′) were synthesized by Integraded DNA Technologies. Expression was calculated using a ΔΔCt method, normalizing to 18S ribosomal RNA levels. Data represent relative expression levels in LIPA and PNPLA2 pooled duplex transfected cells to that of nontargeting siRNA (Dharmacon/GE) control transfected cells.
Lipid extraction for mass spectrometry (MS) lipidomics
MS-based lipid analysis was performed by Lipotype as described previously (Sampaio et al., 2011). Lipids were extracted using a two-step chloroform/methanol procedure (Ejsing et al., 2009). Samples were spiked with internal lipid standard mixture containing cardiolipin 16:1/15:0/15:0/15:0, ceramide 18:1;2/17:0, DAG 17:0/17:0, hexosylceramide 18:1;2/12:0, lyso-phosphatidate 17:0, lyso-PC (LPC) 12:0, lyso-PE (LPE) 17:1, lyso-phosphatidylglycerol 17:1, lyso-PI 17:1, lyso-phosphatidylserine 17:1, phosphatidate 17:0/17:0, PC 17:0/17:0, PE 17:0/17:0, phosphatidylglycerol 17:0/17:0, PI 16:0/16:0, phosphatidylserine 17:0/17:0, CE 20:0, sphingomyelin 18:1;2/12:0;0, and TAG 17:0/17:0/17:0. After extraction, the organic phase was transferred to an infusion plate and dried in a speed vacuum concentrator. First step dry extract was resuspended in 7.5 mM ammonium acetate in chloroform/methanol/propanol (1:2:4, vol:vol:vol) and second step dry extract in 33% ethanol solution of methylamine in chloroform/methanol (0.003:5:1; vol:vol:vol). All liquid handling steps were performed using Hamilton Robotics STARlet robotic platform with the Anti Droplet Control feature for organic solvents pipetting.
MS data acquisition for lipidomics
Samples were analyzed by direct infusion on a QExactive mass spectrometer (Thermo Fisher Scientific) equipped with a TriVersa NanoMate ion source (Advion Biosciences). Samples were analyzed in both positive and negative ion modes with a resolution of Rm/z = 200 = 280,000 for MS and Rm/z = 200 = 17,500 for tandem MS (MSMS) experiments, in a single acquisition. MSMS was triggered by an inclusion list encompassing corresponding MS mass ranges scanned in 1-D increments (Surma et al., 2015). Both MS and MSMS data were combined to monitor CE, DAG, and TAG ions as ammonium adducts; PC and PC O- as acetate adducts; and CL, phosphatidate, PE, PE O-, phosphatidylglycerol, PI, and phosphatidylserine as deprotonated anions. MS only was used to monitor lyso-phosphatidate, LPE, LPE O-, lyso-PI, and lyso-phosphatidylserine as deprotonated anions; and ceramide, hexosylceramide, sphingomyelin, LPC, and LPC O- as acetate adducts.
Lipidomics data analysis and postprocessing
Data were analyzed with in-house developed lipid identification software based on LipidXplorer (Herzog et al., 2011, 2012). Data postprocessing and normalization were performed using an in-house developed data management system. Only lipid identifications with a signal-to-noise ratio greater than five, and a signal intensity fivefold higher than in corresponding blank samples were considered for further data analysis.
Online supplemental material
Fig. S1 shows the validation of ATGL and LAL siRNA knockdown and the accumulation of triglycerides after lipase inhibition. Fig. S2 shows the effect of the HSL/MGL inhibitor CAY10499 on LD size, number, and content in the presence and absence of LALi. Fig. S3 shows the association of lipolysis and lipophagy markers with LDs in AML12 cells treated with the lysosome inhibitor chloroquine. Fig. S4 shows additional biochemical and lipidomic characterization of LD1 and LD2 fractions from AML12 cells and whole rat livers. Fig. S5 shows transmission electron micrographs detailing the ultrastructure of LDs within lipophagic organelles from freshly isolated primary rat hepatocytes.
Acknowledgments
We thank the members of the McNiven laboratory, especially Hong Cao, Bridget Mehall, and Katherine Johnson, for helpful discussion and technical assistance.
This work was supported by National Institutes of Health grants R01DK044650 (M.A. McNiven), R01AA020735 (M.A. McNiven and C.A. Casey), T32DK007352 (R.J. Schulze and M.B. Schott), and K99AA026877 (M.B. Schott); the Mayo Clinic Robert and Arlene Kogod Center on Aging (M.B. Schott); the U.S. Department of Veterans Affairs (C.A. Casey); and a Pilot and Feasibility Award to M.B. Schott from the Mayo Clinic Center for Cell Signaling in Gastroenterology (P30DK084567).
The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. The authors declare no competing financial interests.
Author contributions: M.B. Schott, S.G. Weller, R.J. Schulze, and M.A. McNiven conceived the study, designed experiments, interpreted data, and wrote the manuscript. M.B. Schott, S.G. Weller, R.J. Schulze, E.W. Krueger, and K. Drizyte-Miller performed and analyzed experiments. C.A. Casey provided vital reagents and critical expertise.