Impact Factor 3.367 | CiteScore 4.3
More on impact ›

ORIGINAL RESEARCH article

Front. Physiol., 17 July 2020 | https://doi.org/10.3389/fphys.2020.00850

Lipophagy Impairment Is Associated With Disease Progression in NAFLD

  • 1Laboratory of Microscopic and Ultrastructural Anatomy, University Campus Bio-Medico, Rome, Italy
  • 2Predictive Molecular Diagnostic Division, Department of Pathology, Campus Bio-Medico University Hospital, Rome, Italy
  • 3Department of Biology, University of Rome Tor Vergata, Rome, Italy
  • 4Research Unit of Pathology, University Campus Bio-Medico, Rome, Italy
  • 5Internal Medicine, Geriatrics, and Hepatology Unit, University Campus Bio-Medico, Rome, Italy
  • 6IRCCS Fondazione Santa Lucia, Rome, Italy

Non-alcoholic fatty liver disease (NAFLD) is the most common cause of chronic liver disease in Western countries and is associated with aging and features of metabolic syndrome. Lipotoxicity and oxidative stress are consequent to dysregulation of lipid metabolism and lipid accumulation, leading to hepatocyte injury and inflammation. Lipophagy consists in selective degradation of intracellular lipid droplets by lysosome and mounting evidence suggests that lipophagy is dysregulated in NAFLD. Here we demonstrate lipophagy impairment in experimental models of NAFLD and in a NAFLD patient cohort by histomorphological and molecular analysis. High fat diet-fed C57BL/6J male mice and high-fat/high-glucose cultured Huh7 cells showed accumulation of both p62/SQSTM1 and LC3-II protein. In 59 NAFLD patients, lipid droplet-loaded lysosomes/lipolysosomes and p62/SQSTM1 clusters correlated with NAFLD activity score (NAS) and with NAS and fibrosis stage, respectively, and levels of expression of lysosomal genes, as well as autophagy-related genes, correlated with NAS and fibrosis stage. An increased amount of lipid droplets, lipolysosomes and autophagosomes was found in subjects with NAFLD compared to healthy subjects at ultrastructural level. In conclusion, here we observed that NAFLD is characterized by histological, ultrastructural and molecular features of altered autophagy that is associated with an impaired lipid degradation. Impaired autophagy is associated with features of advanced disease. Lipopolysosomes, as individuated with light microscopy, should be further assessed as markers of disease severity in NAFLD patients.

Introduction

Non-alcoholic fatty liver disease (NAFLD) is becoming the most common hepatic disorder in Western populations and is emerging as one of the principal causes for liver cirrhosis, hepatocellular carcinoma (HCC), and liver transplantation (Bedogni et al., 2014). NAFLD is significantly associated with conditions of hyperglycemia, dyslipidemia, obesity and other features that characterize the metabolic syndrome specially occurring in older individuals (Stefan et al., 2008). Recent studies suggest that advanced age leads to more severe histological changes and lower clinical response as hepatic lipids accumulate over time (Mofrad et al., 2003; Angulo, 2010; Bouchi et al., 2016; Hashimoto et al., 2017; Jinjuvadia et al., 2017; Leung et al., 2017; Sookoian and Pirola, 2017). Indeed, the main feature of NAFLD is steatosis, which consists in the overload of intracellular neutral lipids in the form of lipid droplets (LDs) due lipid homeostasis disruption (Gluchowski et al., 2017). In particular, the uptake and synthesis of lipids overwhelm the ability of the cell to degrade or secrete lipids. The majority of NAFLD patients have simple steatosis or steatosis with non-specific changes and do not progress or evolve very slowly to advanced liver disease (Vernon et al., 2011). Conversely, a subgroup of patients with NAFLD presents the histological criteria for non-alcoholic steatohepatitis (NASH) and has a significant risk of progression toward cirrhosis and HCC (Vernon et al., 2011).

Expanding the knowledge of factors involved in the development and progression of hepatic steatosis and NAFLD is mandatory in order to enhance the repertoire of tools for assessing disease severity and optimizing treatment. LD autophagy, the so-called lipophagy, is the main process responsible for lipid catabolism in the liver. It is a finely regulated step mechanism that consists in: (1) protein-mediated sequestration of cargo (portion or whole LD) within cytosolic double-membrane vesicles, and consequent formation of a phagosome; (2) tracking and transport of autophagosome to lysosome; (3) fusion between autophagosome and lysosome and formation of the autophagolysosome; (4) cargo degradation by lysosomal lipases.

In the early stage of the process, i.e., during the formation of autophagosomal membranes, cytosolic LC3 (LC3-I) is conjugated to phosphatidylethanolamine (PE) by the action of two consecutive ubiquitylation-like reactions catalyzed by the E1-like enzyme Atg7 and the E2-like enzyme Atg3, in order to form LC3-II. When autophagosomes fuse with lysosomes, LC3-II is degraded by lysosomal proteases. Therefore, the increase of LC3-II level is connected with its impaired turnover and is used as a marker of reduced autophagy (Klionsky et al., 2016). p62/sequestosome (SQSTM)1 is recruited as an autophagy receptor required for selective macroautophagy connecting polyubiquitinated cargo with autophagosomes. Consequently, p62/SQSTM1 accumulation represents another hallmark of autophagy blockage.

In the last years, an increasing number of studies have shown that dysregulation of lipophagy plays a pathogenetic role in the development of NAFLD, and that all the steps of this process can be affected (Patrick and Lake, 1969; Kwanten et al., 2014; Martinez-Lopez and Singh, 2015; Schulze et al., 2017). Notably, a previous study analyzed the accumulation of p62/SQTSM1 as marker of autophagy blockage in hepatocytes of NAFLD patients (Fukushima et al., 2018). Finally, lysosomes are crucial organelles for the final phase of fat disposal, and lipolysosomes represent LD-loaded lysosomes in which lipid catabolism takes place thanks to specific hydrolases. They can be individuated as vacuoles surrounded by lysosomal related-proteins, constituting an electron dense membrane at the ultrastructural analysis (Hayashi et al., 1977, 1983; Iancu et al., 2013).

In the present study, we have assessed autophagy both in in vitro and in in vivo models of NAFLD as well as in the liver of NAFLD patients. The results presented give effort to the notion that lipophagy/autophagy is impaired in NAFLD to an extent that is tightly correlated with disease severity and progression (Berson and Pessayre, 1997; Tandra et al., 2011).

Materials and Methods

Patients

Fifty-nine patients who had undergone liver biopsy for suspected NASH at the Hepatology Unit of the University Hospital Campus Bio-Medico of Rome and whose paraffin-embedded liver tissue was available for further analysis, were included in the study. It was confirmed that there was no history of alcoholic intake >20 gr/day if woman and >30 gr/day if man or of use of drugs known to induce liver damage, and that the following results had been obtained: negative anti-HCV antibodies and HBsAg; antinuclear antibodies (ANA) < 1:80 and negative anti-mitochondrial (AMA), anti-smooth muscle (ASMA) and anti-liver and kidney microsomal (anti-LKM) antibodies; normal transferrin saturation and serum levels of ceruloplasmin and alpha-1 antitrypsin. All patients had signed an informed consent in which the possible use of part of their hepatic tissue for research purposes was specified. The protocol of the study conformed to the ethical guidelines of the 1975 Declaration of Helsinki and was approved by the Ethics Committee of the University Campus Bio-Medico of Rome.

Liver biopsies from ten patients with normal transaminases, negative for HBsAg and anti-HCV antibodies tests, and presenting histologically normal liver with less than 5% steatosis were included in the study as control group for histological and immunohistochemical evaluations. These samples were from patients who had undergone liver surgery for metastatic tumors (liver tissue for immunohistochemistry was at least 0.5 cm from the metastatic lesion).

Animal Model, Treatment and Biochemical Analyses

Mouse experimentation was conducted in accordance with accepted standard of humane animal care after the approval by relevant local (Institutional Animal Care and Use Committee, Tor Vergata University) and national (Ministry of Health) committees. Adult C57BL/6J (3 months-age-old) male mice (purchased from ENVIGO, Italy) were randomly divided into four groups: mice fed with normal diet (ND: 3.85 kcal/g among which 10% kcal from fat, 20% from protein and 70% from carbohydrate) for 4 months (group 1) or for 8 months (group 2); mice fed with high fat diet (HFD: 5.24 kcal/g among which 60% kcal from fat, 20% from protein, and 20% from carbohydrate) for 4 months (group 3) or for 8 months (group 4). ND (#D12450B) and HFD (#D12492) were from Research Diets, INC (New Brunswick, NJ, United States). ND contained 38% sugars, 4% saturated fats, 6% unsaturated fats and casein; HFD contained 20% sugars, 54% saturated fats, 6% unsaturated fats and casein. Mice were maintained at 23.0 ± 1.0°C and 55.0 ± 5.0% relative humidity under a 12/12 h light/dark cycle (lights on at 6:00 AM, lights off at 6:00 PM). At the end of treatment, mice were starved for 6 h before blood samples collections and body weight measurement; successively, mice were sacrificed by cervical dislocation, and liver tissues were explanted for the analyses. Blood glucose, cholesterol, alanine transaminase (ALT), aspartate transaminase (AST), and triglycerides were measured through the automatized KeyLab analyser (BPCBioSed, Italy) using specific colorimetric assay kits (BPCBioSed).

Cells and Treatments

The hepatocyte-derived human carcinoma cell line Huh7 was obtained from the American Type Culture Collection (Rockville, MD, United States) and cultured in low-glucose Dulbecco’s Modified Eagle’s Medium (DMEM) with 10% fetal bovine serum, 1% penicillin-streptomycin, and 1% L-glutamine in a humidified incubator at 37°C with 5% CO2. For treatments, cells were plated at 80% confluence and, after 24 h, culture medium was replaced with complete high-glucose DMEM and 400 mM of BSA-conjugated palmitic acid was added for 72 and 96 h.

Western Blot

Animal liver tissue was homogenized in RIPA buffer (50 mM Tris–HCl, pH 8.0, 150 mM NaCl, 12 mM deoxycholic acid, 0.5% Nonidet P-40, and protease and phosphatase inhibitors). Ten μg proteins were loaded on SDS-PAGE and subjected to Western blotting. Nitrocellulose membranes were incubated with anti-heat shock protein-60 (Hsp60) (Abcam, ab46798), anti-p62/SQSTM1 (Abcam, Ab91526), anti-LC3 (Sigma-Aldrich, L7543) primary antibodies at 1:1000 dilution. Successively, membranes were incubated with the appropriate horseradish peroxidase-conjugated secondary antibodies. Immunoreactive bands were detected by a FluorChem FC3 System (Protein-Simple, San Jose, CA, United States) after incubation of the membranes with ECL Selected Western Blotting Detection Reagent (GE Healthcare, Pittsburgh, PA, United States).

Liver Pathology

For both human and animal liver tissue, tissue sections were stained with haematoxylin and eosin and Sirius red and each case was assessed by a single operator. In patients, the severity of steatosis (0–3), features of NASH, including necroinflammation (0–3), hepatocellular ballooning (0–2), and fibrosis stage (0–4) were assessed according to the NAFLD clinical research network scoring system (Kleiner et al., 2005), and each case was diagnosed as simple steatosis or NASH according to Brunt et al. (1999).

Immunohistochemistry and Immunofluorescence

Each case was analyzed by immunohistochemistry for LAMP1, a lysosomal marker (ab24170, 1:50 titre, rabbit polyclonal, Abcam, Cambridge, United Kingdom) and for p62/SQSTM1 (ab56416, 1:100 titre, mouse monoclonal, Cambridge, United Kingdom), a marker of autophagic impairment. Immunohistochemistry was performed on sections obtained from formalin-fixed tissues embedded in paraffin. Immunohistochemical stained sections were scanned using the Hamamatsu NanoZoomer 2.0-RS.

LAMP1 was used to visualize vesicles of lysosomal origin (with a LAMP1-positive membrane). The number of large LAMP1-positive vesicles, with a visible lumen, i.e., lipolysosomes, was evaluated in five fields randomly chosen for each case with a magnification of 400×. A mean value was therefore derived in each sample for statistical analysis.

For the evaluation of p62/SQSTM1 expression, five fields randomly chosen were analyzed (magnification ×200) for each sample. The immunostaining was evaluated by using an open source plugin IHC Profiler for Image J (Varghese et al., 2014). All images analyzed were thresholded at the same intensity value. A scale that ranges from 0 to 3 was applied to describe the intensity of staining in each field, where 0 means negative and 3 strongly positive. The percentage of liver tissue with different score of intensity were recorded for each field and an immunohistochemistry score was derived as follows: [(3+%tissue) X3 + (2+%tissue) X2 + (1+%tissue) X1]. A mean value was then calculated for each sample (Potts et al., 2012).

Immunofluorescent staining was performed on sections from paraffin-embedded tissue and cell culture with the same antibodies used for immunohistochemistry. Secondary Alexa Fluor 488 and/or Alexa Fluor 568-conjugated donkey anti-rabbit, anti-mouse, anti-goat antibodies (Invitrogen, Carlsbad, CA, United States) were added for 30 min (1:200 titre), followed by another extensive washing step in TBS. Nuclear counterstaining was performed using Vectashield mounting medium (Vector Laboratories, Burlingame, CA, United States) containing 4,6-diamidino-2-phenylindole (DAPI). Fluorescence images were collected with a Nikon A1 Confocal Laser Microscope System (Nikon, Tokyo, Japan). Acquisition was carried out using the Imaging Software NIS-Elements (Nikon).

Ultrastructural Analysis

A subgroup (n = 5) of liver samples were fixed in 2.5% glutaraldehyde in 0.1 M cacodylate buffer for the ultrastructural study. After fixation, dehydration and impregnation, the samples were included in epoxy resins and acrylic, cut at the ultramicrotome and processed for the ultrastructural study of hepatocytes by electron microscopy. In particular, morphological aspects typical of lipid overload, lysosomal activation and lipophagy/autophagy process, such as the presence of lipid droplets, membrane-bound lipid vacuoles (lipolysosomes) and autophagosomes were evaluated in 5000× fields. The presence of lipolysosomes in hepatocytes of NAFLD patients by both light and electron microscopy was evaluated according to previous criteria by Iancu et al. (2013).

mRNA Analysis

For RT-qPCR analysis, total RNA was extracted by using Quick-RNA extraction kit (Zymo Research). RNA (3 μg) was retro-transcribed by using M-MLV (Promega, Madison, WI). qPCR was performed in triplicate by using validated qPCR primers (BLAST), Ex TAq qPCR Premix and the Real-Time PCR (Applied Biosystem). mRNA levels were normalized to Alpha-fetoprotein mRNA, and the relative mRNA levels were determined through the 2–ΔΔCt method. Primers used for RT-qPCR are listed in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Primer sequences used for RT-qPCR analyses.

Nanostring nCounter analysis was performed in a subgroup of patients and controls (n = 38). Total RNA was extracted using the High Pure FFPET RNA Isolation Kit (Roche) from 10 μm sections of FFPE liver biopsies of NAFLD patients (27 specimens) and 11 controls. The RNA was quantified with the BioPhotometer® D30 (Eppendorf) and 100 ng of total RNA was subjected to Nanostring nCounter Analysis System. This system utilizes a novel digital color-coded barcode technology that is based on direct multiplexed measurement of gene expression and offers high levels of precision and sensitivity (<1 copy per cell). The technology uses molecular “barcodes” and single molecule imaging to detect and count hundreds of unique transcripts in a single reaction.

We used a custom nCounter CodeSet containing 20 genes, selected in order to evaluate the expression levels (mRNA levels) of genes involved in autophagy and of other functional genes in NAFLD versus healthy liver samples. In particular, of these genes, 6 are lysosome related genes (LAMP1, LAMP2, CTSD, NPC1, ATP6V1B2, TFEB) and 9 are key genes involved in autophagy (mTOR, ULK1, BECN1, ATG12, ATG5, ATG7, p62/SQSTM1, RAB7A, GABARAP) (Table 2). We have also included in the panel 3 hepatocyte markers (AFP, ALB, and TTR) to normalize the values for the number of hepatocytes present in each sample and 2 housekeeping genes (TUBB and HIST1H3A), to normalize the values for total number of cells. Nanostring nCounter analysis was performed following the manufacturer’s protocol. In brief, 100 ng of total RNA was hybridized in solution at 65°C for 16 h with specific pairs of ∼50 base probes for each gene set mRNA. The Reporter Probe carries the signal; the Capture Probe allows the complex to be immobilized for data collection. After hybridization, the probe excess was removed and the probe/target complexes aligned and immobilized in the nCounter Cartridge. Sample Cartridges were placed in the Digital Analyzer for data collection. Color codes on the surface of the cartridge will be counted and tabulated for each target molecule. Raw counts for each mRNA was normalized (for hepatocyte markers or for housekeeping genes) and analyzed using nSolver Analysis Software (nSAS) to obtain differential expression analysis of 44 genes in NAFLD samples versus control (healthy) samples.

TABLE 2
www.frontiersin.org

Table 2. Gene, Protein name and Function are derived from https://www.uniprot.org.

Statistics

Depending on the parametric or nonparametric distribution, variables are expressed as mean ± SD or median and 25–75% interquartile range (25–75% IR), respectively. Differences were evaluated by the Kruskal–Wallis test. Correlations were carried out by the Spearman’s rank correlation test. A p < 0.05 was considered statistically significant. SPSS software (version 22.00; SPSS Inc., Chicago, IL, United States) was used for statistical analyses.

Results

Progressive Impairment of Lipophagy in Cell and Mouse Models of Liver Steatosis

In order to generate a preclinical model of NAFLD, C57BL/6J mice were fed with HFD for 4 or 8 months. Mice fed with normal diet (ND) served as controls. As expected, HFD mice underwent a progressive increase of body weight (Figure 1A) and alteration of biochemical parameters such as increased blood glucose, cholesterol and ALT (Figure 1B). An increase of liver weight (Figure 1C), and the acquisition of a steatotic phenotype (Figure 1D) were also observed, in association with enhanced expression of inflammatory cytokines (IL-6, TNF-α, IL-1β) and fibrosis markers (Col1a1, Mmp2) (Figure 1E), collectively recapitulating the occurrence of a NAFLD-like clinical phenotype. At the protein level, a significant increase of both p62/SQSTM1 and LC3-II/LC3-I ratio was noticed compared to ND group that was more marked after 8 months than after 4 months of HFD treatment (Figures 1F,G).

FIGURE 1
www.frontiersin.org

Figure 1. C57BL/6J male mice were fed with high fat diet (HFD) or normal diet (ND) for 4 or 8 months. Increase of body weight (A), alterations of biochemical parameters (B), increase of liver mass (C) accompanied by histological steatosis (D) were observed in HFD mice compared to ND. The mRNA expression of inflammation- and tissue remodeling-related genes (E) and autophagy protein markers LC3-II/LC3-I ratio and p62/SQSTM1 (F) were progressively increased as demonstrated by densitometric analyses (G), indicating autophagy inhibition. HSP60 was used as loading control and representative immunoblots are reported. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001 vs. ND; °p < 0.05, °°p < 0.01 vs. 4 mo. (n = 6 mice per group). Original magnification: × 200. Calibration bar: 25 μm.

Similar results were observed in the in vitro model of hepatocellular fat overload. In particular, Huh7 cells were cultured in high-glucose DMEM supplemented with 400 μM palmitic acid (high-fat/high-glucose, HFHG). Starting at 72 h of treatment, up-regulation of cytokine mRNAs (IL-6, TNF-α, IL-1β), increase of both p62/SQSTM1 protein and LC3-II/LC3-I ratio (Figures 2A,B) as well as accumulation of lysosomal mass at confocal analysis with the lysosomal marker LAMP1 (Figure 2C) were observed, suggesting autophagy inhibition and lysosomal reaction.

FIGURE 2
www.frontiersin.org

Figure 2. Huh7 human cells treated with 400 μM palmitate and 4.5 g/L glucose (HFHG) show accumulation of p62/SQSTM1 protein and increase of LC3-II/LC3-I protein ratio (A) as well as up-regulation of inflammatory cytokines (B) and lysosomal mass (C) after 72 and 96 h of treatment. TOM20 was used as loading control. Densitometric analyses of p62/SQSTM1 and LC3-II/LC3-I are shown in the bottom panel (A). n = 3; ***p < 0.001, *p < 0.05 vs. Ctr; °°°p < 0.001 vs. HFHG 72 h, ###p < 0.001 vs. Ctr 96 h.

Increased Lipolysosomes and Impaired Lipid Catabolism in NAFLD Patients

Based on the results obtained in the experimental models of NAFLD, we moved to a cohort of 59 NAFLD patients. The main characteristics of these patients are reported in Table 3.

TABLE 3
www.frontiersin.org

Table 3. Biological, anthropometric, and biochemical characteristics of the NAFLD patients.

By immunohistochemistry for LAMP1, lysosomes from healthy subjects were detected mainly near the hepatocyte biliary pole, and no LAMP1-positive vacuoles with visible lumens were detected. In NAFLD patients, the lysosomal compartment architecture was completely disarranged. In particular, lysosomes were no longer located at the biliary pole, and LAMP1-positive vacuoles were observed inside hepatocytes. Moreover, the number of LD-loaded lysosomes, the so-called lipolysosomes, positively correlated with the NAS score (ρ = 0.375, p < 0.005) (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Immunohistochemical expression of LAMP1 in human hepatocyte lysosomes (LAMP1-positive particles) were regularly disposed along the path of biliary canaliculi in healthy liver (A), while they were patchy distributed in fatty liver (B). In fatty liver of NAFLD patients (n = 59), enlarged fat-loaded lysosomes, namely lipolysosomes, could be individuated as intracellular droplets with lysosomal LAMP1-positive membrane and visible lumen (B, arrowheads). The number of lipolysosomes/field at light microscopy increased in NAFLD patients with higher NAS (Non-alcoholic fatty liver disease Activity Score) (C). Original magnification: × 400 (A,B). *p < 0.05, **p < 0.01.

The increase of lysosomal mass was supported by increased levels of expression of lysosomal genes (CTSD, LAMP1, LAMP2, NPC1, ATP6V1B2,TFEB), whose levels of expression were correlated with the NAS score (r = 0.7; p < 0.001;r = 0.5; p < 0.005;r = 0.4; p < 0.01;r = 0.4; p < 0.05; r = 0.6; p < 0.001; r = 0.4; p < 0.01, respectively), as well as with fibrosis stage (r = 0.6; p < 0.001; r = 0.7; p < 0.00; r = 0.6; p < 0.001; r = 0.5; p < 0.005, r = 0.8; p < 0.001; r = 0.6; p < 0.001, respectively) (Figure 4).

FIGURE 4
www.frontiersin.org

Figure 4. Box plots of mRNA expression of lysosome related genes in association with NAS score (on the left) and fibrosis stage (on the right); (n = 38). For correlation purpose Spearman’s ρ test was performed.

Impaired Autophagy in NAFLD Patients

Consistent with the results obtained in mice, p62/SQSTM1 protein levels as well as fibrosis were significantly increased in the liver of NAFLD patients. In particular, immunohistochemical positivity of p62/SQSTM1 increased with progression of fibrosis, accumulating in hepatocytes in the so-called Mallory-Denk bodies, and indicating autophagy blockage. Furthermore, by means of immunofluorescence, p62/SQSTM1 aggregates were observed outside the lysosomes, suggesting and impairment also of the first steps of autophagy (Figure 5). Moreover, significant correlations between levels of expression of autophagy-related genes (ATG12, ATG5, ATG7, BECN1, GABARAP, MTOR, RAB7A, P62/SQSTM1, ULK1) and the NAS score (r = 0.4; p < 0.05; r = 0.4; p < 0.01; r = 0.5; p < 0.005; r = 0.5; p < 0.005; r = 0.3; p < 0.05; r = 0.6; p < 0.001; r = 0.3; p < 0.05; r = 0.5; p < 0.005; r = 0.3; p < 0.05, respectively) and fibrosis stage were observed (r = 0.7; p < 0.001; r = 0.6; p < 0.001; r = 0.7; p < 0.001; r = 0.7; p < 0.001; r = 0.5; p < 0.005; r = 0.7; p < 0.001; r = 0.5; p < 0.001; r = 0.5; p < 0.005; r = 0.4; p < 0.01, respectively) (Figure 6).

FIGURE 5
www.frontiersin.org

Figure 5. p62/SQSTM1 accumulates with NAFLD progression (A–C). P62/SQSTM1 positive Mallory-Denk bodies (arrowheads) were individuated in hepatocytes. Double labeling experiments with a lysosome marker (LAMP1, red) and autophagy impairment marker (p62/SQSTM1 green), confirmed that p62/SQSTM1 accumulates outside the lysosomes (D). n = 59. Scale bar: 25 μm; * = p < 0,05 vs. NAS 0–2 or mild fibrosis.

FIGURE 6
www.frontiersin.org

Figure 6. Box plots of mRNA expression of autophagy related genes in association with NAS score (on the left) and fibrosis stage (on the right); (n = 38). For correlation purpose Spearman’s ρ test was performed.

Evaluation of autophagy was also performed at ultrastructural level in human biopsies processed for TEM experiment. An increased amount of lipid droplets, lipolysosomes and autophagosomes was found in subjects with NAFLD compared to healthy subjects (p < 0.05) (Figure 7).

FIGURE 7
www.frontiersin.org

Figure 7. Transmission electron microscopy analysis confirmed the presence of lipolysosomes with dense lipid core and typical electron-dense lysosomal membrane in fatty liver (A–D). The different cellular organelles involved in the autophagic process phases are depicted in the high power fields (C,E). Bars representing the number of autophagic organelles in NAFLD patients compared to healthy subjects (mean values, n = 5) (F). LD, lipid droplet; LL, lipolysosome; AP, autophagosome; ER, stressed endoplasmic reticulum; BC, bile canaliculus; red circle, autolysosome; white arrowhead, phagophore; black arrowhead, immature autophagosome; white arrow, lysosome; white asterisk, stressed mitochondria; black asterisk, lipid component of autolysosome. Original magnification: 3900× (A); 2300× (B); 12000× (C,D); 46000× (E). **p < 0.05.

NAFLD Progression and Aging in NAFLD Patients

Histological markers of NAFLD progression were associated with aging in NAFLD patients. Standard histopathological features of NAFLD severity, such as grading and staging of the disease (NAS and fibrosis score, respectively) were significantly different in NAFLD patients above and below 54 years (Table 4). Furthermore, immunohistochemical positivity of p62/SQSTM1 was correlated with aging in NAFLD patients (r = 0.3, p < 0.05).

TABLE 4
www.frontiersin.org

Table 4. Correlation of NAFLD characteristics with patients’ age.

Discussion

Aging and obesity predispose to the development of NAFLD. Actually, the chronic overload of lipids has harmful effects in the liver, exposing resident hepatocytes to their cytotoxic action (lipotoxicity) and oxidative stress, the effect of which are more pronounced with aging (Zelber-Sagi et al., 2006; Chen et al., 2013; Xu et al., 2013a). Lipophagy represents the main lipid degrading system in the liver and its impairment was recently associated with the pathogenesis of NAFLD (Sinha et al., 2014; Zhang et al., 2018). In lipophagy, the neutral lipids are partially subtracted in autophagosomes and transported to lysosomes, which represent the final destination of lipids in this pathway.

In the present study, we have demonstrated that chronic lipid overload through HFD in mice and palmitic acid treatment in a human hepatocyte derived cell line, respectively, leads to a time-dependent inhibition of autophagy. Notably, in HFD-fed mice, the progression of NAFLD-mimicking features was strongly associated with the degree of autophagy inhibition (p62/SQSTM1 and LC3-II/LC3-I increase). Moreover, here we showed, for the first time, the presence of lipolysosomes in hepatocytes of NAFLD patients by both light and electron microscopy. In NAFLD, lipolysosomes likely represent the progressive impairment of the lysosomal function, and in particular of the capability of lysosomal hydrolases to catabolize fat. This is supported by the correlation between the increased number of lipolysosomes and disease activity in terms of necroinflammation.

The up-regulation of lysosomal genes, such as CTSD, LAMP1, LAMP2, NPC1, ATP6V1B2, TFEB, suggesting an increase of the overall lysosomal mass, should be well interpreted as an attempt to counteract lysosomal dysfunction. Among these genes, notably, CTSD has been already correlated with inflammation and impairment of lipid metabolism in NAFLD. Houben et al. (2017) have shown that inhibition of CTSD in mice results in reduced cholesterol and triglyceride levels in the liver. The group of Liao showed that CTSD is up-regulated as cholesterol accumulates in lysosomal compartments in an NPC1 knock-out model, demonstrating an increase of CTSD when lysosomes are engulfed (Liao et al., 2007). In the present study, CTSD expression was higher in NAFLD patient compared to controls, and a direct correlation was observed between CTSD expression and both NAS and fibrosis.

Lysosome engulfment by lipids has been shown to inhibit autophagosome turnover, and several studies in lysosomal storage diseases have demonstrated an impairment in the fusion between autophagosomes and fatty lysosomes that leads to defective autophagic clearance (Fukuda et al., 2006; Settembre et al., 2008; Sarkar et al., 2013). Consistent with this, as already described (Fukuo et al., 2014), we confirmed an increased expression of p62/SQSTM1 mRNA and protein in the liver of NAFLD patients that was significantly associated with disease activity and fibrosis stage. The increase of p62/SQSTM1 suggests that the autophagic cargos cannot fuse with lysosomes in order to be degraded, leading to cytosolic aggregates of p62/SQSTM1; accordingly, double labeling experiment revealed that p62/SQSTM1 clusters are detached from lysosomal-derived structure. Furthermore, in NAFLD patients, TEM analysis revealed an increased number of autophagosomes and lipolysosomes, likely representing the compensatory expansion of a bad-functioning autophagic machinery. Indeed, the expression of several genes involved in autophagy (Atg7, Atg5, GABARAB, and Rab7) was significantly correlated with NAS and fibrosis stage. These are key proteins for autophagic vesicle formation and lysosome maturation, and their deficit was found to suppress autophagic flux and to impair the autophagosome-lysosome fusion (Xu et al., 2013b). Moreover, autophagy could exert an anti-fibrogenic role since it promotes survival of hepatocytes and hepatocyte apoptosis is a central event in the fibrogenic process (Ni et al., 2014; Ruart et al., 2019; Hammoutene et al., 2020).

In conclusion, in the present study, we have observed histological, ultrastructural and molecular features of altered autophagy in NAFLD, leading to impaired lipid degradation. Impaired autophagy seems to be a pivotal factor determining the clinical progression of NAFLD also considering that a decrease in autophagic activity was previously described with aging, which is associated with NASH progression (Rubinsztein et al., 2011). Individuating molecular markers of lipophagy impairment could help to identify those patients at risk of progression. In particular, in this study, LAMP1-positive vacuoles with visible lumen, known as lipolysosomes, were identified as morphological feature of altered lipophagy. Further studies are clearly awaited in order to verify if quantification of lipolysosomes could be a feasible and reproducible tool for assessing the severity of NAFLD and its propensity to progress.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Ethics Statement

The studies involving human participants were reviewed and approved by Comitato Etico dell’Università Campus Bio-Medico di Roma with authorization n° 30/15 OSS ComEt CBM. The patients/participants provided their written informed consent to participate in this study. The animal study was reviewed and approved by the University Animal Welfare Committee-OPBA, Tor Vergata University with authorization no 378/2017-PR.

Author Contributions

SC designed and supervised the study, performed the morphological experiments and the statistical analysis, and wrote the manuscript. FZ performed molecular experiments with Nanostring and gave an intellectual contribution. SR and FV performed the immunohistochemical experiments and provided intellectual contribution. KA designed molecular experiments, provided intellectual contribution, and helped writing the manuscript. MZ contributed to the design of the electron microscopy study. MF performed fluorescence microscopy experiments. GP provided tissue samples from UCBM Hospital. FA helped in recruiting patients’ cohort from UCBM database. RA-I contributed to the critical revision of the manuscript. AP contributed to the study design and provided intellectual contribution. SM supervised the study and contributed to the critical revision of the manuscript. DL-B performed in vivo experiments, analyzed molecular data, and provided substantial intellectual contribution. UV-G obtained the patients informed consents, contributed substantially to the design of the study, to the discussion of the data, and to the writing of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was partially supported by funds from European Foundation for the Studies on Diabetes to DL-B.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

AFP, alpha-fetoprotein; ALB, albumin; ATG12, ubiquitin-like protein ATG12; ATG7, ubiquitin-like modifier-activating enzyme ATG7; ATP6V1B2, v-type proton ATPase subunit B, brain isoform.; BECN1, beclin-1; CTSD, cathepsin D; GABARAP, gamma-aminobutyric acid receptor-associated protein; HFD, high-fat diet; HIST1H3A, histone cluster 1H3 family member a; LAMP1, lysosome-associated membrane glycoprotein 1; LAMP2, lysosome-associated membrane glycoprotein 2; mTOR, serine/threonine-protein kinase mTOR; NAS, NAFLD activity score; NASH, non-alcoholic steatohepatitis; NPC1, NPC intracellular cholesterol transporter 1; nSAS, nSolver Analysis Software; p62/SQSTM1, sequestosome-1; RAB7A, ras-related protein Rab-7a; TBS, tris buffered saline; TEFB, transcription factor EB; TG5, autophagy protein 5; TTR, transthyretin; TUBB, tubulin-beta; ULK1, serine/threonine-protein kinase ULK1.

References

Angulo, P. (2010). Long-term mortality in nonalcoholic fatty liver disease: is liver histology of any prognostic significance? Hepatology 51, 373–375. doi: 10.1002/hep.23521

PubMed Abstract | CrossRef Full Text | Google Scholar

Bedogni, G., Nobili, V., and Tiribelli, C. (2014). Epidemiology of fatty liver: an update. World J. Gastroenterol. 20, 9050–9054. doi: 10.3748/wjg.v20.i27.9050

PubMed Abstract | CrossRef Full Text | Google Scholar

Berson, A., and Pessayre, D. (1997). Microvesicular steatosis and steatohepatitis: role of mitochondrial dysfunction and lipid peroxidation. J. Hepatol. 26, 13–22. doi: 10.1016/s0168-8278(97)82328-8

CrossRef Full Text | Google Scholar

Bouchi, R., Takeuchi, T., Akihisa, M., Ohara, N., Nakano, Y., Nishitani, R., et al. (2016). Increased visceral adiposity with normal weight is associated with the prevalence of non-alcoholic fatty liver disease in Japanese patients with type 2 diabetes. J. Diabetes Investig. 7, 607–614. doi: 10.1111/jdi.12443

PubMed Abstract | CrossRef Full Text | Google Scholar

Brunt, E. M., Janney, C. G., Bisceglie, A. M., Neuschwander-Tetri, B. A., and Bacon, B. R. (1999). Nonalcoholic steatohepatitis: a proposal for grading and staging the histological lesions. Am. J. Gastroenterol. 94, 2467–2474. doi: 10.1111/j.1572-0241.1999.01377.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Q. R., Braun, R., Hu, Y., Yan, C., Brunt, E. M., Meerzaman, D., et al. (2013). Multi-SNP analysis of gwas data identifies pathways associated with nonalcoholic fatty liver disease. PLoS One 8:e65982. doi: 10.1371/journal.pone.0065982

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukuda, T., Ewan, L., Bauer, M., Mattaliano, R. J., Zaal, K., Ralston, E., et al. (2006). Dysfunction of endocytic and autophagic pathways in a lysosomal storage disease. Ann. Neurol. 59, 700–708. doi: 10.1002/ana.20807

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukuo, Y., Yamashina, S., Sonoue, H., Arakawa, A., Nakadera, E., Aoyama, T., et al. (2014). Abnormality of autophagic function and cathepsin expression in the liver from patients with non-alcoholic fatty liver disease. Hepatol. Res. 44, 1026–1036. doi: 10.1111/hepr.12282

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukushima, H., Yamashina, S., Arakawa, A., Taniguchi, G., Aoyama, T., Uchiyama, A., et al. (2018). Formation of p62-positive inclusion body is associated with macrophage polarization in non-alcoholic fatty liver disease. Hepatol. Res. 48, 757–767. doi: 10.1111/hepr.13071

PubMed Abstract | CrossRef Full Text | Google Scholar

Gluchowski, N. L., Becuwe, M., Walther, T. C., and Farese, R. V. (2017). Lipid droplets and liver disease: from basic biology to clinical implications. Nat. Rev. Gastroenterol. Hepatol. 14, 343–355. doi: 10.1038/nrgastro.2017.32

PubMed Abstract | CrossRef Full Text | Google Scholar

Hammoutene, A., Biquard, L., Lasselin, J., Tanguy, M., Loyer, X., Jules, M., et al. (2020). A defect in endothelial autophagy occurs in patients with non-alcoholic steatohepatitis and promotes inflammation and fibrosis. J. Hepatol. 2020 72, 528–538. doi: 10.1016/j.jhep.2019.10.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Hashimoto, Y., Hamaguchi, M., Fukuda, T., Ohbora, A., Kojima, T., and Fukui, M. (2017). Fatty liver as a risk factor for progression from metabolically healthy to metabolically abnormal in non-overweight individuals. Endocrine 57, 89–97. doi: 10.1007/s12020-017-1313-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayashi, H., Sameshima, Y., Lee, M., Hotta, Y., and Kosaka, T. (1983). Lipolysosomes in human hepatocytes: their increase in number associated with serum level of cholesterol in chronic liver diseases. Hepatology 3, 221–225. doi: 10.1002/hep.1840030214

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayashi, H., Winship, D. H., and Sternlieb, I. (1977). Lipolysosomes in human liver: distribution in livers with fatty infiltration. Gastroenterology 73, 651–654. doi: 10.1016/s0016-5085(19)31759-7

CrossRef Full Text | Google Scholar

Houben, T., Oligschlaeger, Y., Hendrikx, T., Bitorina, A. V., Walenbergh, S. M. A., Van Gorp, P. J., et al. (2017). Cathepsin D regulates lipid metabolism in murine steatohepatitis. Sci. Rep. 7:3494. doi: 10.1038/s41598-017-03796-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Iancu, T. C., Manov, I., Shaoul, R., Haimi, M., and Lerner, A. (2013). What’s in a name?-”Lipolysosome”: ultrastructural features of a lipid-containing organelle. Ultrastruct. Pathol. 37, 293–303. doi: 10.3109/01913123.2013.799625

PubMed Abstract | CrossRef Full Text | Google Scholar

Jinjuvadia, R., Antaki, F., Lohia, P., and Liangpunsakul, S. (2017). The association between nonalcoholic fatty liver disease and metabolic abnormalities in the united states population. J. Clin. Gastroenterol. 51, 160–166. doi: 10.1097/MCG.0000000000000666

PubMed Abstract | CrossRef Full Text | Google Scholar

Kleiner, D. E., Brunt, E. M., Van Natta, M., Behling, C., Contos, M. J., Cummings, O. W., et al. (2005). Design and validation of a histological scoring system for nonalcoholic fatty liver disease. Hepatology 41, 1313–1321. doi: 10.1002/hep.20701

PubMed Abstract | CrossRef Full Text | Google Scholar

Klionsky, D. J., Abdelmohsen, K., Abe, A., Abedin, M. J., Abeliovich, H., Arozena, A. A., et al. (2016). Guidelines for the use and interpretation of assays for monitoring autophagy (3rd edition). Autophagy 12, 1–222. doi: 10.1080/15548627.2015.1100356

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwanten, W. J., Martinet, W., Michielsen, P. P., and Francque, S. M. (2014). Role of autophagy in the pathophysiology of nonalcoholic fatty liver disease: a controversial issue. World J. Gastroenterol. 20, 7325–7338. doi: 10.3748/wjg.v20.i23.7325

PubMed Abstract | CrossRef Full Text | Google Scholar

Leung, J. C. F., Loong, T. C. W., Wei, J. L., Wong, G. L. H., Chan, A. W. H., Choi, P. C. L., et al. (2017). Histological severity and clinical outcomes of nonalcoholic fatty liver disease in nonobese patients. Hepatology 65, 54–64. doi: 10.1002/hep.28697

PubMed Abstract | CrossRef Full Text | Google Scholar

Liao, G., Yao, Y., Liu, J., Yu, Z., Cheung, S., Xie, A., et al. (2007). Cholesterol accumulation is associated with lysosomal dysfunction and autophagic stress in Npc1-/- mouse brain. Am. J. Pathol. 171, 962–975. doi: 10.2353/ajpath.2007.070052

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinez-Lopez, N., and Singh, R. (2015). Autophagy and lipid droplets in the liver. Annu. Rev. Nutr. 35, 215–237. doi: 10.1146/annurev-nutr-071813-105336

PubMed Abstract | CrossRef Full Text | Google Scholar

Mofrad, P., Contos, M. J., Haque, M., Sargeant, C., Fisher, R. A., Luketic, V. A., et al. (2003). Clinical and histologic spectrum of nonalcoholic fatty liver disease associated with normal ALT values. Hepatology 37, 1286–1292. doi: 10.1053/jhep.2003.50229

PubMed Abstract | CrossRef Full Text | Google Scholar

Ni, H. M., Woolbright, B. L., Williams, J., Bryan, C., Wei, C., James, L. P., et al. (2014). Nrf2 promotes the development of fibrosis and tumorigenesis in mice with defective hepatic autophagy. J. Hepatol. 61, 617–625. doi: 10.1016/j.jhep.2014.04.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Patrick, A. D., and Lake, B. D. (1969). Deficiency of an acid lipase in wolman’s disease. Nature 222, 1067–1068. doi: 10.1038/2221067a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Potts, S. J., Krueger, J. S., Landis, N. D., Eberhard, D. A., David Young, G., Schmechel, S. C., et al. (2012). Evaluating tumor heterogeneity in immunohistochemistry-stained breast cancer tissue. Lab. Investig. 92, 1342–1357. doi: 10.1038/labinvest.2012.91

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruart, M., Chavarria, L., Campreciós, G., Nuria, S. H., Bosch, J., Scott, L. F., et al. (2019). Impaired endothelial autophagy promotes liver fibrosis by aggravating the oxidative stress response during acute liver injury. J. Hepatol. 70, 458–469. doi: 10.1016/j.jhep.2018.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubinsztein, D. C., Mariño, G., and Kroemer, G. (2011). Autophagy and aging. Cell 146, 682–695. doi: 10.1016/j.cell.2011.07.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarkar, S., Carroll, B., Buganim, Y., Maetzel, D., Ng, A. H. M., Cassady, J. P., et al. (2013). Impaired autophagy in the lipid-storage disorder niemann-pick type C1 disease. Cell Rep. 5, 1302–1315. doi: 10.1016/j.celrep.2013.10.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Schulze, R. J., Drižytë, K., Casey, C. A., and McNiven, M. A. (2017). Hepatic lipophagy: new insights into autophagic catabolism of lipid droplets in the liver. Hepatol. Commun. 1, 359–369. doi: 10.1002/hep4.1056

PubMed Abstract | CrossRef Full Text | Google Scholar

Settembre, C., Fraldi, A., Jahreiss, L., Spampanato, C., Venturi, C., Medina, D., et al. (2008). A block of autophagy in lysosomal storage disorders. Hum. Mol. Genet. 17, 119–129. doi: 10.1093/hmg/ddm289

PubMed Abstract | CrossRef Full Text | Google Scholar

Sinha, R. A., Farah, B. L., Singh, B. K., Siddique, M. M., Li, Y., Wu, Y., et al. (2014). Caffeine stimulates hepatic lipid metabolism by the autophagy-lysosomal pathway in mice. Hepatology 59, 1366–1380. doi: 10.1002/hep.26667

PubMed Abstract | CrossRef Full Text | Google Scholar

Sookoian, S., and Pirola, C. J. (2017). Systematic review with meta-analysis: risk factors for non-alcoholic fatty liver disease suggest a shared altered metabolic and cardiovascular profile between lean and obese patients. Aliment. Pharmacol. Ther. 46, 85–95. doi: 10.1111/apt.14112

PubMed Abstract | CrossRef Full Text | Google Scholar

Stefan, N., Kantartzis, K., and Häring, H. U. (2008). Causes and metabolic consequences of fatty liver. Endocr. Rev. 29, 939–960. doi: 10.1210/er.2008-0009

PubMed Abstract | CrossRef Full Text | Google Scholar

Tandra, S., Yeh, M. M., Brunt, E. M., Vuppalanchi, R., Cummings, O. W., Ünalp-Arida, A., et al. (2011). Presence and significance of microvesicular steatosis in nonalcoholic fatty liver disease. J. Hepatol. 55, 654–659. doi: 10.1016/j.jhep.2010.11.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Varghese, F., Bukhari, A. B., Malhotra, R., and De, A. (2014). IHC Profiler: an open source plugin for the quantitative evaluation and automated scoring of immunohistochemistry images of human tissue samples. PLoS One. 9:e96801. doi: 10.1371/journal.pone.0096801

PubMed Abstract | CrossRef Full Text | Google Scholar

Vernon, G., Baranova, A., and Younossi, Z. M. (2011). Systematic review: the epidemiology and natural history of non-alcoholic fatty liver disease and non-alcoholic steatohepatitis in adults. Aliment. Pharmacol. Ther. 34, 274–285. doi: 10.1111/j.1365-2036.2011.04724.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, C., Yu, C., Ma, H., Xu, L., Miao, M., and Li, Y. (2013a). Prevalence and risk factors for the development of nonalcoholic fatty liver disease in a nonobese Chinese population: the Zhejiang Zhenhai study. Am. J. Gastroenterol. 108, 1299–1304. doi: 10.1038/ajg.2013.104

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X., Hua, Y., Sreejayan, N., Zhang, Y., and Ren, J. (2013b). Akt2 knockout preserves cardiac function in high-fat diet-induced obesity by rescuing cardiac autophagosome maturation. J. Mol. Cell Biol. 5, 61–63. doi: 10.1093/jmcb/mjs055

PubMed Abstract | CrossRef Full Text | Google Scholar

Zelber-Sagi, S., Nitzan-Kaluski, D., Halpern, Z., and Oren, R. (2006). Prevalence of primary non-alcoholic fatty liver disease in a population-based study and its association with biochemical and anthropometric measures. Liver Int. 26, 856–863. doi: 10.1111/j.1478-3231.2006.01311.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Yao, Z., Chen, Y., Qian, L., Jiang, S., Zhou, J., et al. (2018). Lipophagy and liver disease: new perspectives to better understanding and therapy. Biomed. Pharmacother. 97, 339–348. doi: 10.1016/j.biopha.2017.07.168

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: autophagy, high-fat diet, lipolysosomes, non-alcoholic steatohepatitis, lipophagy

Citation: Carotti S, Aquilano K, Zalfa F, Ruggiero S, Valentini F, Zingariello M, Francesconi M, Perrone G, Alletto F, Antonelli-Incalzi R, Picardi A, Morini S, Lettieri-Barbato D and Vespasiani-Gentilucci U (2020) Lipophagy Impairment Is Associated With Disease Progression in NAFLD. Front. Physiol. 11:850. doi: 10.3389/fphys.2020.00850

Received: 30 January 2020; Accepted: 24 June 2020;
Published: 17 July 2020.

Edited by:

Keith Russell Brunt, Dalhousie University, Canada

Reviewed by:

Zhihong Yang, Université de Fribourg, Switzerland
Shisan Bao, The University of Sydney, Australia

Copyright © 2020 Carotti, Aquilano, Zalfa, Ruggiero, Valentini, Zingariello, Francesconi, Perrone, Alletto, Antonelli-Incalzi, Picardi, Morini, Lettieri-Barbato and Vespasiani-Gentilucci. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Simone Carotti, s.carotti@unicampus.it; Katia Aquilano, katia.aquilano@uniroma2.it

These authors have contributed equally to this work