Amber Light (590 nm) Induces the Breakdown of Lipid Droplets through Autophagy-Related Lysosomal Degradation in Differentiated Adipocytes

Lipolysis in the adipocytes provides free fatty acids for other tissues in response to the energy demand. With the rapid increase in obesity-related diseases, finding novel stimuli or mechanisms that regulate lipid metabolism becomes important. We examined the effects of visible light (410, 457, 505, 530, 590, and 660 nm) irradiation on lipolysis regulation in adipocytes differentiated from human adipose-derived stem cells (ADSCs). Interestingly, 590 nm (amber) light irradiation significantly reduced the concentration of lipid droplets (LDs). We further investigated the lipolytic signaling pathways that are involved in 590 nm light irradiation-induced breakdown of LDs. Immunoblot analysis revealed that 590 nm light irradiation-induced phosphorylation of hormone-sensitive lipase (HSL) was insufficient to promote reduction of LDs. We observed that 590 nm light irradiation decreased the expression of perilipin 1. We found that 590 nm light irradiation, but not 505 nm, induced conversion of LC3 I to LC3 II, a representative autophagic marker. We further demonstrated that the lysosomal inhibitors leupeptin/NH4Cl inhibited 590 nm light irradiation-induced reduction of LDs in differentiated adipocytes. Our data suggest that 590 nm light irradiation-induced LD breakdown is partially mediated by autophagy-related lysosomal degradation, and can be applied in clinical settings to reduce obesity.

As the main reservoir of the body's fuel, white adipose tissues play a critical role in energy storage and balance. During fasting or exercise, mature adipocytes provide energy as free fatty acids for other tissues through a process known as lipolysis 1 . Decreased lipolysis activity in adipose tissues may lead to triglyceride (TG) accumulation in lipid droplets (LDs), causing unnecessary additional energy intake, which leads to obesity 2,3 . Thus, tight control of lipolysis is required not only for energy homeostasis, but also for the prevention of metabolic diseases.
Lipolysis is a series of processes involving several regulatory events 4 . These include the lipolytic and anti-lipolytic mechanisms induced by lipolytic hormones such as adrenocorticotropic hormone (ACTH), anti-lipolytic hormones such as insulin, or LD-associated proteins such as perilipins and lipases 4,5 . Perilipin 1 has two modes of action: downregulating basal lipolysis and mediating protein kinase A (PKA)-activated lipolysis. Perilipin 1 is highly phosphorylated by cyclic AMP-dependent PKA 6 . A study showed that basal lipolysis is increased and PKA-stimulated lipolysis increase is blocked in adipocytes of perilipin 1 knock-out mice 7,8 .
Lipases play an important role in lipolysis. Even though there are many proteins predicted to have lipase or esterase activity in adipose tissues, hormone-sensitive lipase (HSL) and adipose triglyceride lipase (ATGL) are responsible for most of the triacylglycerol (TAG) hydrolase activity in murine adipocytes 9 . Several papers have reported that ATGL is the rate-limiting enzyme for the initiation of PKA-stimulated lipolysis and predominantly responsible for the first step of TAG hydrolysis, whereas HSL is primarily responsible for the hydrolysis of diacylglycerol (DAG) [10][11][12] . Other groups, however, have suggested that ATGL is less important in human adipocytes and that HSL fulfills the function as the rate-limiting TAG hydrolase in human adipocytes 13,14 .
There are several reports on the effects of several visible lights on various cells. Blue light exposure decreased proliferation of keratinocytes 15,16 . Red light exposure induces increased proliferation in keratinocytes 17 and fibroblasts 18 . In addition, some visible lights are involved in the regulation of barrier recovery. Red light (550~670 nm) can accelerate recovery in a disrupted skin model, whereas blue light (430-510 nm) retards this recovery 19 . We also reported that violet light reduces the expression of keratinocyte differentiation markers through rhodopsin 20 .
In the present study, we examined the effects of visible light irradiation on adipocytes differentiated from human adipose-derived stem cells (ADSCs). We found that 590 nm light irradiation induced the breakdown of LDs in adipocytes, and further examined the underlying mechanism to explain the phenomenon.

Results
Visible light irradiation at 590 nm reduced LDs in differentiated adipocytes. To investigate the effects of visible light irradiation (410, 457, 505, 530, 590, or 660 nm) on fully differentiated adipocytes (2 weeks), cells were irradiated with each wavelength 5 times (once a day, for 5 days, which showed the maximum effects), and LDs in adipocytes were measured by performing Oil Red O staining analysis. Interestingly, 590 nm light irradiation significantly reduced LDs by about 45% compared with control (Fig. 1a). Although 660 nm light irradiation also enhanced the breakdown of LDs by 40%, its effect was less potent than that of 590 nm light irradiation (Fig. 1a,b). We confirmed that 590 nm light irradiation had no cytotoxic effect on adipocyte survival by performing CCK-8 assay (Fig. 1c). In addition, we examined whether 590 nm light irradiation has any effect on adipocyte differentiation. Cells were irradiated with 590 nm visible light every 3 days (day 0, 3, and 6) during the differentiation process. Oil red O staining analysis revealed that 590 nm visible light show slight decrease of LDs in 6 days . LDs in adipocytes were stained with Oil Red O dye (a) and quantified (* * * p < 0.001, ANOVA followed by Dunnett's post hoc test; n = 22) (b). Cell viability after light exposure was evaluated using the CCK-8 assay (* * p < 0.01, * * * p < 0.001, ANOVA followed by Bonferroni post hoc test; n = 12) (c). Adipocytes were irradiated with 590 nm or 660 nm of visible light every 3 days during the differentiation process until 6 days. IDXT was used for positive control (* * p < 0.01, ANOVA followed by Dunnett's post hoc test; n = 18) (d). Data were normalized by setting the control as 1 (b,d) or 100% (c) and presented as means ± S.D. Scale bar = 200 μ m. N.S. = not significant (p > 0.05).
Scientific RepoRts | 6:28476 | DOI: 10.1038/srep28476 during adipocyte differentiation (Fig. 1d). In addition, we confirmed that 590 nm light irradiation reduced LDs in a time-dependent manner (Fig. S1). These results indicated that 590 nm light irradiation most potently induced the breakdown of LDs in differentiated adipocytes.
Irradiation at 590 nm did not affect the phosphorylation level of HSL. For a better understanding of the mechanism involved in 590 nm light irradiation-induced decrease of LDs, we evaluated the effect of the irradiation on HSL, a key lipase regulated by reversible phosphorylation 9,11,14 . Our data revealed that 590 nm irradiation leads to a slight increase in HSL phosphorylation compared to other light wavelengths (Fig. 2b-e). However, comparing with the pattern of HSL phosphorylation by forskolin (Fig. 2a), the induction level or the duration of HSL phosphorylation by 590 nm irradiation was remarkably low or short. In addition, treatment with either H-89, a PKA inhibitor, or KT5823, a PKG inhibitor, did not affect 590 nm irradiation-induced breakdown of LDs (Fig. 2g).
Irradiation at 590 nm specifically reduced the expression level of perilipin 1. In response to lipolytic stimulation by forskolin or β -adrenergic agents, cAMP-dependent PKA-mediated phosphorylation of perilipin 1 facilitates lipolysis-inducing HSL phosphorylation 11,21 . However, in the absence of lipolytic stimuli, namely, in basal or constitutive lipolysis, perilipin 1 negatively regulates lipase function 22,23 . Recently, it has been reported that a decrease in perilipin 1 protein expression could cause lipolysis under basal conditions 24 . Thus, we examined if perilipin 1 expression was affected by visible light irradiation. Interestingly, the expression level of perilipin 1 once for 30 min, and were supplemented with fresh media. Cell lysates were collected at the indicated time points (0, 0.5, 1, 3, 6 hr) and subjected to immunoblot for total HSL, phospho HSL (Ser552, Ser554 and Ser650), and GAPDH. These images are representative of three independent experiments. Differentiated adipocytes were irradiated with 590 nm light 3 times (18 J/cm 2 a day, for 3 days) in the presence of H89 (1, 5 or 10 μ M) or KT5823 (0.5, 1, 2 μ M) at each concentration. LDs in adipocytes were stained with Oil Red O dye and quantified (g). * * p < 0.01, compared to control by t-test. N.S., compared to 590 nm light irradiation-only group. ANOVA followed by Dunnett's post hoc test.
Scientific RepoRts | 6:28476 | DOI: 10.1038/srep28476 was significantly decreased by either 590 or 660 nm light irradiation (Fig. 3a), which was consistent with Oil-Red O staining results, as shown in Fig. 1. According to previous reports, forskolin increases the phosphorylation of perilipin 1 through cAMP-PKA signaling pathway 5,6 . Forskolin is not likely to affect perilipin 1 protein expression or degradation in adipocytes 25 . We also demonstrated that forskolin did not alter perilipin 1 expression compared to GAPDH expression (Fig. 3a).
The mass of adipose tissue can be reduced by increasing lipolysis, inhibiting adipogenesis or both 26 . Since 590 nm light irradiation decreased LD levels, we examined the expression levels of ATGL, monoacylglycerol lipase (MGL), and comparative gene identification-58 (CGI-58) that reside in LDs. We found that there were no changes in the expression level of ATGL, MGL, or CGI-58, despite a reduction in perilipin 1 (Fig. 3b). We further examined whether adipogenic markers were affected by 590 nm light irradiation in differentiated adipocytes. Although we could not find any difference in protein expression of adiponectin, acetyl-CoA carboxylase, retinoid X receptor α (RXR α ) and CCAAT/enhancer-binding protein α (C/EBPα ), we observed a slight reduction in both fatty acid synthase (FAS) and fatty acid binding protein 4 (FABP4) expression in the 590 nm light irradiated adipocytes (Fig. 3c). Next, we questioned whether the slight decrease of FAS expression could cause the inhibition of de novo lipogenesis resulting in the breakdown of LDs. It has been reported that the inhibition of FAS activity decreases LDs in adipocytes 27 . Thus, we tested whether 590 nm irradiation could change the FAS activity. We found that FAS activity in differentiated adipocytes after 590 nm light irradiation was slightly but not significantly decreased (Fig. S2), implying that 590 nm light irradiation-induced the breakdown of LDs are mainly due to lipolysis rather than the inhibition of de novo lipogenesis.
To further investigate whether lipid metabolism was altered by increased lipolysis, we examined the expression patterns of genes related to energy metabolism in irradiated adipocytes by qRT-PCR analysis 28 . Increases in expression of genes related to lipid metabolism such as acyl-CoA dehydrogenase medium chain (ACADM), acyl-CoA oxidase 1 (ACOX1) and carnitine palmitoyl transferase 1a/1b (CPT1A, CPT1B) were observed at 590 nm or 660 nm irradiated cells (Fig. S3).
Irradiation at 590 nm enhanced the breakdown of LDs through autophagy-related lysosomal degradation. Singh et al. suggested that LDs and autophagy components are interrelated 29 . Thus, to investigate whether the effects of 590 nm irradiation on the breakdown of LDs is mediated by autophagy pathway in differentiated adipocytes, we first examined the protein levels of beclin-1, Atg5-Atg12 conjugate, p62/SQSTM1, and LC-II, which are representative autophagy markers 29 . We found that the protein levels of beclin-1, Atg5-Atg12 conjugate, and LC3-II were increased after 590 nm light irradiation, whereas the protein level of p62/SQSTM1 was reduced (Fig. 4a). During the autophagy process, LC3-II is also characterized by a punctate pattern. We confirmed that LC3-II puncta was markedly increased compared with control at 3 hr after 590 nm light irradiation (Fig. 4b).

Figure 3. Effect of visible light irradiation on perilipin 1 expression in differentiated adipocytes.
Differentiated adipocytes were incubated with forskolin or irradiated with 457 nm, 505 nm, 530 nm, 590 nm and 660 nm of visible light (18 J/cm 2 ) 3 times (once a day, for 3 days). Cell lysates were collected and subjected to immunoblot for perilipin and GAPDH. Protein-transferred membrane was stained by Ponceau S staining before immunoblotting (a). Differentiated adipocytes were irradiated with 590 nm of visible light three times (18 J/cm 2 ), and after 24 hr in fresh media, the cell lysates were collected and subjected to immunoblot for LD markers; ATGL, perilipin, MGL, and CGI-58 (b) and for adipogenic markers; adiponectin, FAS, FABP4, acetyl-CoA carboxylase, RXRα , C/EBPα (c). GAPDH was used as a loading control. These images were representative of three independent experiments. To examine whether 590 nm specifically shows this phenomenon, we compared 590 nm light irradiation with 505 nm light irradiation. Interestingly, 590 nm light irradiation but not 505 nm enhanced the LC3-II levels at 1 hr post-irradiation (Fig. 4c). This increase in the LC3-II level could be due to an increased autophagosome formation or impaired autophagosome-lysosome fusion (autophagosome maturation) 30 . To discriminate between these two possibilities, we measured the level of LC3-II in the presence of both NH 4 Cl and leupeptin, which are used as a lysosome-neutralizing agent and a lysosomal proteolysis inhibitor, respectively 30 . Light irradiation at 590 nm in the presence of both NH 4 Cl and leupeptin increased the level of LC3-II in differentiated adipocytes (Fig. 5a), indicating that the increased LC3-II expression was due to increased autophagosome formation. In addition, 590 nm light irradiation-induced reduction of LDs was blocked in the presence of both NH 4 Cl and leupeptin (Fig. 5b,c).

Discussion
In the present study, we showed that 590 nm irradiation significantly reduced LD levels in differentiated adipocytes (Fig. 1).
To elucidate the underlying mechanism, we examined the expression or phosphorylation status of photoreceptors, HSL, and perilipin 1. To find out the photoreceptor that responds to 590 nm visible light, we analyzed the membrane fraction of differentiated adipocytes on each day (0, 3, 7 days) using liquid-chromatography-tandem mass spectrometry, but we unfortunately could not observe any photoreceptor (Fig. S4).
According to Nilsson et al., if the 'surge' of HSL phosphorylation level falls sharply, the subsequent lipolysis can also shut down immediately 31 . This report implies that HSL phosphorylation induced by 590 nm irradiation could be too low and short to induce lipolysis sufficiently. Additionally, we confirmed that 590 nm irradiation-induced reduction of LDs was not affected by a PKA inhibitor or a PKG inhibitor, indicating that HSL was not significantly involved with adipocyte lipolysis by 590 nm irradiation.
Perilipins and lipases such as ATGL and its co-lipase CGI-58 are responsible for fat mobilization 26 . Intriguingly, 590 nm irradiation specifically reduced the expression level of perilipin 1 (Fig. 3), which indicates the contribution of increased lipolysis 7,8 . Meanwhile, patients with loss-of-function mutations in perilipin 1 exhibit elevated basal lipolysis. This led to severe metabolic consequences including partial lipodystrophy, severe insulin resistance, diabetes, dyslipidemia, and nonalcoholic fatty liver disease 32,33 , implying that reduced perilipin 1 level by 590 nm could be deleterious for normal metabolism. However, there are two differences between these previous data and our present study: (1) Lipolytic condition at 590 nm irradiation is a temporally acute, but not a 'basal' condition caused by continuous effects like a loss of function. (2) Though we found similar results with those of previous studies showing the reduction of LDs following the decreased expression level of perilipin 1, the significance of our data is limited only to local areas irradiated by 590 nm visible light. Consequently, the 590 nm irradiation employed in this study is more likely to induce the breakdown of LDs only in a particular area, rather than affect whole body metabolism by, for example, loss-of-function mutations in perilipin 1.
We found that 590 nm irradiation significantly reduced perilipin 1 protein levels, but not mRNA levels ( Fig. 3 and Table S1). It is interesting to note that the reduction in perilipin 1 protein caused by tumor necrosis factor α treatment is associated with reduced perilipin 1 mRNA levels, whereas nelfinavir, a human immunodeficiency virus protease inhibitor, reduces perilipin 1 protein level, but does not decrease perilipin 1 mRNA 24,34 . Mechanisms controlling stability of perilipin 1 protein in adipocytes are largely unknown but may be functionally important because changes in the protein content of perilipin 1 are not always associated with alterations in mRNA levels 16,35,36 . A report demonstrated that perilipin 1 degradation is mediated by a ubiquitination-proteasome pathway, which enables post-translational regulation of perilipin 1 37 . However, we could not observe induction of proteasome-dependent perilipin 1 degradation by 590 nm irradiation because cell death occurred in our experimental condition, which required continuous treatment of MG132, a proteasome inhibitor (data not shown). Therefore, we could not exclude the possibility that perilipin 1 may be degraded via the proteasome pathway after 590 nm irradiation.
Recent studies revealed that a link between LDs and lysosomes is related to the mobilization of intracellular lipid stores. Although cytoplasmic lipases are important for LD degradation in adipocytes, where HSL and ATGL double knock-out abolishes more than 90% of lipolytic activity 9 , a different type of LD degradation process has been identified, named lipophagy 29 . In this process, autophagosomes contribute to the delivery of LDs to lysosomes, leading to their degradation. Recently, Kaushik et al. suggested that elimination of perilipin 2 and

Figure 5. Leupeptin/NH 4 Cl blocked 590 nm-induced breakdown of LDs in differentiated adipocytes.
Differentiated adipocytes were pretreated with or without both NH 4 Cl (20 mM) and leupeptin (Leu; 100 μ M) and then were irradiated with 590 nm (18 J/cm 2 ) once (a) or 3 times (once a day, for 3 days) (b,c). After 2 hr recovery following irradiation, cell lysates were analyzed by immunoblotting with antibodies against LC3 and GAPDH (a). After third irradiation of 590 nm, LDs in adipocytes were stained with Oil Red O dye (b) and quantified (c). These images were representative of three independent experiments. Scale bar = 20 μ m. * * p < 0.01, ANOVA followed by Bonferroni post hoc test (n = 4) was performed on the data of 3 groups (control, inhibitors-only group, and 590 nm irradiation-only group). # p < 0.05, compared to 590 nm irradiation-only group by t-test.
Scientific RepoRts | 6:28476 | DOI: 10.1038/srep28476 perilipin 3 by chaperone-mediated autophagy, which selectively degrades cytosolic proteins containing KFERQ motif in lysosomes, facilitates access of the lipid core to the cytosolic lipase ATGL 38 . Another study revealed that autophagy genes (ATGs) and ATG-positive membranes are detected in LDs during lipolysis, suggesting that autophagosomes form on the surface of LDs 39 . Similarly, we demonstrated that 590 nm irradiation, but not 505 nm, specifically induced accumulation of LC3 II and reduction of LDs. This was inhibited by leupeptin/ NH 4 Cl (Fig. 5). Thus, we assumed that perilipin 1 degradation and reduction of LDs by 590 nm irradiation may be due to lysosomal degradation through the activation of autophagic flux.
In our study, we found that adipocyte-specific proteins such as ATGL, CGI-58, and C/EBPα showed no changes in their expression levels after 590 nm irradiation (Fig. 3). This phenomenon was similarly identified in a previous study. Tansey et al. demonstrated that perilipin null mice exhibited a lean phenotype, but no changes in expression levels of C/EBPα and diglyceride acyltransferase (DGAT) 8 . These data indicate that the decrease of perilipin, which drives the breakdown of LDs, is likely to show no direct relevance with the expression levels of adipocyte specific proteins contributing to form LDs.
A recent study also showed that Atg5 transgenic mice exhibit leanness with loss of weight during the aging process, indicating that the amount of fat in adipocytes is significantly decreased 40 . By next-generation sequencing, we further demonstrated that 590 nm-irradiated adipocytes showed higher expression of Atg5 gene and down-regulation of the expression of p62 when compared to control adipocytes (Table 1).
In conclusion, we suggest that 590 nm irradiation enhances the reduction of LDs in differentiated adipocytes, potentially via autophagy-related lysosomal degradation of perilipin 1. Our scientific approach to explain how 590 nm visible light could induce the breakdown of LDs may shed light on a potential non-invasive treatment for local fat reduction.
LED irradiation system. The LED irradiation system, which emits light with a narrow bandwidth inside a CO 2 incubator, was purchased from HanaroTR (Suwon, Korea). Internal LED arrays, which emit violet light (380-420 nm, λ max : 410 nm), green light (450-560 nm, λ max : 505 nm), orange light (560-620 nm, λ max : 590 nm) and red light (620-690 nm, λ max : 660 nm), could be controlled from 10 mW/cm 2 for 30 min at one time when the object was 50 mm from the LED arrays. These arrays were obtained from Roithner Lasertechnik GmbH (Wien, Austria). We used a spectroradiometer SPR-4001 (Luzchem, Ottawa, Canada) to determine the intensity of the irradiation.
Light exposure and treatment with inhibitors. The cells were washed twice with PBS prior to irradiation and irradiated in PBS with the indicated exposure lights to avoid the side effects associated with photosensitization in the culture medium. To avoid any confounding effects from other factors (temperature, humidity and CO 2 ), the LED arrays were placed inside the CO 2 cell culture incubator during exposure. The cells were then irradiated with LED lights at 18 J/cm 2 . For the inhibition experiments, cells were pre-incubated for 1 hr with each inhibitor.
Cell viability tests. Adipocyte viability after light exposure was evaluated using the Cell Counting Kit (CCK)-8 assay according to the manufacturer's instructions (Dojindo Molecular Technologies. Inc, Kumamoto, Japan). The cell viability detecting reagent, WST-8 [2-(2-methoxy-4-nitrophenyl)-3-(4-nitrophenyl)-5-(2,4-disulfophenyl)-2H-tetrazolium, monosodium salt] was added to the cells in culture, and the cells were incubated for 2 hr in a humidified atmosphere. The absorbance at 450 nm was then measured, and the cell viability was expressed as the percentage of the absolute optical density of each sample relative to that of the control value.
Western blot analysis. Adipocytes were lysed in RIPA cell lysis buffer containing a protease inhibitor and phosphatase inhibitor cocktail (Sigma-Aldrich, St Louis, MO, USA). The lysate was then subjected to centrifugation at 15,000 × g for 20 min, and the supernatant was used for the analysis. The protein concentration was determined using the Bradford method with bovine serum albumin as the standard. The proteins (40 μ g/well) were loaded and fractionated using SDS-PAGE and transferred onto nitrocellulose membranes. The membranes were blocked with 5% non-fat skim milk in TBST (10 mM Tris-HCl [pH 8.0], 150 mM NaCl, 0.015% Tween-20) for 1 hr at room temperature and then probed overnight at 4 °C with primary antibodies. All of the blots were washed 3 times with TBST and then probed with horseradish peroxidase-conjugated goat anti-mouse or -rabbit IgG secondary antibodies (Bio-Rad, Hercules, CA, USA) at room temperature for 1 hr. The membranes were developed using ECL solution (Amersham Pharmacia Biotech, Piscataway, NJ, USA). Immunofluorescence. Adipocytes were fixed with 10% formaldehyde for immunostaining for 10 min and permeabilized with 0.1% Triton X-100 for 10 min. The blocking and diluent solutions consisted of PBS-T containing 1% BSA (Sigma-Aldrich, St. Louis, MO, USA). The cells were blocked at room temperature for 30 min and incubated for overnight at 4 °C with LC3 II antibody followed by incubation with an Alexa 488 secondary antibody (Invitrogen, Carlsbad, CA, USA) for 60 min. DAPI was used for the nuclei staining (Sigma-Aldrich, St. Louis, MO, USA). PBS-T washes were performed at each step. These images were observed using a Carl Zeiss LSM700 inverted confocal microscope and the image acquisition software ZEN (Carl Zeiss, Germany).
Next-generation sequencing and data analysis. In order to construct cDNA libraries with the TruSeq RNA library kit, 1ug of total RNA was used. The protocol consisted of polyA-selected RNA extraction, RNA fragmentation, random hexamer primed reverse transcription and 100 nucleotide paired-end sequencing by Illumina HiSeq2000. The libraries were quantified using qPCR according to the qPCR Quantification Protocol Guide and qualified using an Agilent Technologies 2100 Bioanalyzer. To estimate expression levels and to find alternative spliced transcripts, the RNA-Seq reads were mapped to the human genome using TopHat 41 , which is capable of reporting split-read alignments across splice junctions and determined using Cufflinks software 42 in default options. The reference genome sequence (hg19, Genome Reference Consortium GRCh37) and annotation data were downloaded from the UCSC website (http://genome.uscs.edu). The transcript counts in isoform level and gene level were calculated, and the relative transcript abundances were measured in FPKM (Fragments Per Kilobase of exon per Million fragments mapped) using Cufflinks. Raw data were extracted as FPKM values for each gene for each sample by cufflinks software. Transcripts with at least one zeroed FPKM values across all samples were excluded. We added 1 with FPKM value to facilitate log2 transformation. Filtered data was transformed by logarithm and normalized by quantile method. Differentially expressed genes (DEG) were adjusted |fold change|> = 2 for each comparison pair. All NGS data analysis and visualization of differentially expressed genes was conducted using R 2.15.1 (www.r-project.org).  Table 1. Analysis of autophagy-related genes in 590 nm light irradiated adipocytes by next generation sequencing.