Impact of senescence on the transdifferentiation process of human hepatic progenitor-like cells

Abstract

BACKGROUND

Senescence is characterized by a decline in hepatocyte function, with impairment of metabolism and regenerative capacity. Several models that duplicate liver functions in vitro are essential tools for studying drug metabolism, liver diseases, and organ regeneration. The human HepaRG cell line represents an effective model for the study of liver metabolism and hepatic progenitors. However, the impact of senescence on HepaRG cells is not yet known.

AIM

To characterize the effects of senescence on the transdifferentiation capacity and mitochondrial metabolism of human HepaRG cells.

METHODS

We compared the transdifferentiation capacity of cells over 10 (passage 10 [P10]) vs P20. Aging was evaluated by senescence-associated (SA) beta-galactosidase activity and the comet assay. HepaRG transdifferentiation was analyzed by confocal microscopy and flow cytometry (expression of cluster of differentiation 49a [CD49a], CD49f, CD184, epithelial cell adhesion molecule [EpCAM], and cytokeratin 19 [CK19]), quantitative PCR analysis (expression of albumin, cytochrome P450 3A4 [CYP3A4], γ-glutamyl transpeptidase [γ-GT], and carcinoembryonic antigen [CEA]), and functional analyses (albumin secretion, CYP3A4, and γ-GT). Mitochondrial respiration and the ATP and nicotinamide adenine dinucleotide (NAD⁺)/NAD with hydrogen (NADH) content were also measured.

RESULTS

SA β-galactosidase staining was higher in P20 than P10 HepaRG cells; in parallel, the comet assay showed consistent DNA damage in P20 HepaRG cells. With respect to P10, P20 HepaRG cells exhibited a reduction of CD49a, CD49f, CD184, EpCAM, and CK19 after the induction of transdifferentiation. Furthermore, lower gene expression of albumin, CYP3A4, and γ-GT, as well as reduced albumin secretion capacity, CYP3A4, and γ-GT activity were reported in transdifferentiated P20 compared to P10 cells. By contrast, the gene expression level of CEA was not reduced by transdifferentiation in P20 cells. Of note, both cellular and mitochondrial oxygen consumption was lower in P20 than in P10 transdifferentiated cells. Finally, both ATP and NAD⁺/NADH were depleted in P20 cells with respect to P10 cells.

CONCLUSION

SA mitochondrial dysfunction may limit the transdifferentiation potential of HepaRG cells, with consequent impairment of metabolic and regenerative properties, which may alter applications in basic studies.

Key Words: Senescence; HepaRG cells; Transdifferentiation; Mitochondria; Regeneration; Nicotinamide adenine dinucleotide

Core Tip: The human HepaRG cell line represents an effective model for the study of liver metabolism and hepatic progenitors. However, the impact of senescence on HepaRG cells is not known. We characterized the effects of senescence on the transdifferentiation capacity and mitochondrial metabolism of HepaRG cells. By using a replication protocol, we described higher senescence-associated markers and lower transdifferentiation markers in passage 20 (P20) than in P10 cells. Cellular and mitochondrial oxygen consumption, and ATP and nicotinamide adenine dinucleotide (NAD⁺)/NAD with hydrogen (NADH) content were lower in P20 than in P10 transdifferentiated cells. To conclude, senescence-associated mitochondrial dysfunction may limit the transdifferentiation potential of HepaRG cells.

INTRODUCTION

Primary human hepatocytes are the gold standard to study the biology, pharmacology, and toxicology of parenchymal liver cells[1]. Nevertheless, primary hepatocytes present with several limitations such as difficult isolation procedures, high variability between different donors, quick failure in biological function, and proliferation capacity[2]. To overcome these limitations, several human hepatoma cell lines were developed and characterized to provide a steady and unrestricted supply of hepatocyte-like cells. Of these, HepG2 and Huh-7 cells are the most widely used. HepG2 cells originated from an American patient[3], while Huh-7 were isolated from a Japanese patient[4], both affected by well-differentiated hepatocellular carcinoma. Even though these cell lines exhibit several hepatic functions, low expression of enzymes and cytochromes limit their use for metabolism and toxicity studies[5].

To study the metabolism, toxicology, and regeneration/differentiation processes, the HepaRG cell line was used as a replacement for primary hepatocytes, HepG2, and Huh-7 cells[6-8]. Isolated from an Edmonson grade I differentiated liver tumor, HepaRG cells exhibit a hepatocyte-like morphology and express hepatocyte-specific functions in defined culture conditions[9]. Nevertheless, HepaRG cells display features of human oval ductular bipotent hepatic progenitors at confluence[10,11]. Supplementation of dimethyl sulfoxide (DMSO) to confluent HepaRG cells triggers differentiation toward hepatocytes[10,12]. Thus, the behavior of HepaRG cells is exclusive; these cells can be cultured for several passages to proliferate, or stimulated to differentiate towards fully functional hepatocyte-like cells[13].

Cellular senescence consists of a steady cell cycle block occurring because of different harmful events, which include DNA damage, oxidative stress, or even replication[14]. Stem/progenitor cells undergoing senescence cause impairment of tissue homeostasis and regeneration, caused by defective stemness and differentiation processes[15]. The induction of senescence in both HepG2 and Huh-7 cells is associated with altered gene signature, which includes changes in cell cycle regulation, signal transduction, and metabolism[16,17]. Nevertheless, to the best of our knowledge, the impact of senescence on HepaRG cells has not yet been investigated.

Thus, this study investigated whether a replication protocol would induce senescence in HepaRG cells. In addition, we characterized the effects of senescence on the transdifferentiation capacity and mitochondrial metabolism.

MATERIALS AND METHODS

Cell line and culture

The human cell line HepaRG was purchased by Merck Millipore (MMHPR116; Merck KGaA, Darmstadt, Germany). Undifferentiated HepaRG cells exhibit a fibroblast-like morphology, and the differentiation process induces both hepatocyte- and biliary-like epithelial phenotypes at confluence, indicating bipotent progenitor features[11,18]. HepaRG cells were seeded at 27000 cell/cm² confluence in a base medium composed by William’s E Medium + GlutaMAX (3255-020; Thermo Fisher Scientific, Waltham, MA, United States) supplemented with 10% fetal bovine serum (F7524; Merck KGaA), 100 U/mL penicillin (13752; Merck KGaA), and 100 μg/mL streptomycin (P4333; Merck KGaA). Medium was changed twice a week and cells were passaged once every 7 d. Cells in passage 10 (P10, young cells) and P20 (senescent cells) were used for assays and compared. To obtain HepaRG differentiation, a two-step procedure was used as previously described. Cells were cultured in the medium for 2 wk and then in the presence of 2% DMSO for an additional 2 wk[11].

Senescence-associated β-galactosidase activity assay

The senescence-associated (SA) β-galactosidase (SA-β-gal) activity assay was performed according to the manufacturer’s protocol (#9860; Cell Signaling Technology, Inc. Danvers, MA, United States). Briefly, P10 and P20 HepaRG cells grown on 6-well plates were fixed in 1X fixative solution containing 2% formaldehyde and 2% glutaraldehyde for 10 min, and then stained overnight at 37°C with the β-galactosidase staining solution at pH 6.0 for 15 h. Images were acquired using the Nikon Eclipse Ni-U microscope (Nikon, Tokyo, Japan).

Comet assay

The comet assay was performed as previously described[19]. DNA was stained with SYBR green (172-5271; Bio-Rad Laboratories, Hercules, CA, United States) just before blind slide scoring with the Nikon Eclipse Ni-U fluorescence microscope equipped with the CCD-200E video camera. At least 100 cells per sample were analyzed using the Comet Assay IV analysis software (Perceptive Instruments, Haverhill, Suffolk, United Kingdom). The extent of DNA damage in single cells was evaluated by the percentage of tail DNA.

Confocal microscopy and flow cytometry

Cells (1.5 × 10⁵ cells/well) were seeded on a glass coverslip in a 24-multiwell plate. The next day, cells were washed three times with phosphate-buffered saline (PBS), fixed in 4% paraformaldehyde for 10 min at room temperature (RT), and washed twice with PBS. Cells were first permeabilized with PBS + 0.1% X-100 Triton (93418; Merck KGaA) for 10 min, and then incubated in blocking buffer (3% bovine serum albumin, A7906; Merck KGaA) + 0.3 M glycine (G-7126-500-50; Merck KGaA) for 30 min at RT. Subsequently, cells were treated for 1.5 h in the dark at RT with the following labeled antibodies: anti-cluster of differentiation 49a (CD49a) (130-101-397), anti-CD49f (130-097-246), anti-CD184 (130-098-354), anti-epithelial cell adhesion molecule (EpCAM) (130-091-253), and anti-cytokeratin 19 (CK19; 130-080-101). All antibodies were purchased by Miltenyi Biotec B.V. & Co. KG (Bergisch Gladbach, North Rhine-Westphalia, Germany), and labeled with phycoerythrin, except anti-CK19, which was labeled with fluorescein isothiocyanate. Nuclei were counterstained with 4′,6-diamidino-2-phenylindole included in mounting medium (ab104139; Abcam plc, Cambridge, United Kingdom). Cells were analyzed using the Nikon Eclipse Ti-E confocal microscope and by flow cytometry analysis using the FlowSight Cytometer (Amnis; Merck Millipore) and the IDEAS software.

Gene expression analysis

To study the levels of genes expressed in P10 and P20 HepaRG cells after differentiation, RNA was extracted from 1.0 × 10⁶ cells/sample and converted into cDNA, which was used as a template in the following RT-PCR. RNA extraction was performed by the “Pure Link RNA Mini Kit” (12183025; Thermo Fisher Scientific), according to the manufacturer’s protocol. RNA concentration was determined by spectrophotometer method at Nanodrop, measuring absorbance at λ = 260 nm. A260/A280 > 2 was evaluated to guarantee protein-free samples. Reverse transcription was performed using the High-Capacity cDNA Reverse Transcription Kit (4368814; Thermo Fisher Scientific), and SYBR Green (172-5271; Bio-Rad Laboratories, Hercules, CA, United States) was used as a fluorescent probe. The sequences of the forward and reverse primers for all of the genes studied are provided in Table 1.

Table 1 Sequences of forward and reverse primers of the genes studied.

Actin	Human	FOR	5’-TGGACATCCGCAAAGACCTG-3’
		REV	5’-GCCGATCCACACGGAGTACTT-3’
Albumin	Human	FOR	5’-CCTGTTGCCAAAGCTCGATG-3’
		REV	5’-GAAATCTCTGGCTCAGGCGA-3’
CYP3A4	Human	FOR	5’-CTTCATCCAATGGACTGCATAAAT-3’
		REV	5’-TCCCAAGTATAACACTCTACACAG-3’
CEA	Human	FOR	5’-GGTCTTCAACCCAATCAGTAAGAAC-3’
		REV	5’-ATGGCCCCAGGTGAGAGG-3’
γ-GT	Human	FOR	5’-TTTGGTGTGCTGCTGGATGAC-3’
		REV	5’-ACCTGAGCTTCCCCACCTATG-3’

CEA: Carcinoembryonic antigen; CYP3A4: Cytochrome P350 3A4; FOR: Forward; γ-GT: Gamma-glutamyl transpeptidase; REV: Reverse.

Albumin secretion analyses

To assess the secretion of human albumin, supernatants of P10 and P20 HepaRG cells in basal conditions and after DMSO treatment were collected after 24 h of culture and assessed using a Human Albumin ELISA Kit, according to the manufacturer’s instructions (ab108788; Abcam, Cambridge, United Kingdom).

CYP3A4 and γ-glutamyl transpeptidase activities

P10 and P20 HepaRG cells in basal conditions and after DMSO treatment were washed with 1× PBS, and 50 μL of 3 μM P450-Glo™ substrate (V8801; Promega, Waldorf, Germany) was added and incubated for 1 h at 37°C, 5% CO₂. Then, 25 μL substrate medium was transferred to a 96-white plate, and CYP3A4 activity was measured according to the manufacturer’s protocol. HepaRG cells were homogenized in 200 μL ice-cold γ-glutamyl transpeptidase (γ-GT) assay buffer, and the γ-GT Activity Colorimetric Assay Kit (MAK089; Merck KGaA) was used according to the manufacturer’s protocol.

Cell respirometry

HepaRG cells (5.0 × 10⁶ cells) were washed with PBS and resuspended in 10 mmol/L KH₂PO₄, 27 mmol/L KCl, 1 mmol/L MgCl₂, 40 mmol/L HEPES, 0.5 mmol/L EGTA buffer (pH 7.1), and assayed for oxygen consumption by the Oxygraph Plus System (Hansatech Instruments, Norfolk, UK) at 37°C under continuous stirring. Oligomycin (8 μg/mL) was added followed the addition of valinomycin (2 μg/mL) after 5 min. The rates of oxygen consumption were corrected for 3 mmol/L potassium cyanide (KCN)-insensitive respiration and normalized to the cell number. Each experiment was repeated in triplicate.

Intracellular ATP content

Total intracellular ATP was quantified using the ENLITEN^® ATP Assay System (FF2000; Promega Corporation, Madison, WI, United States), according to the manufacturer’s instructions. This assay is based on luciferase, and used as the catalyzing enzyme of the ATP reaction with d-Luciferin. ATP extraction is carried out with trichloroacetic acid reagent, which releases ATP from cells, preventing its enzymatic degradation. After addition of the enzyme reagent to the intracellular extract, light emission was detected by a luminometer at 560 nm (DTX 880 Microplate Reader; Beckman Coulter, Brea, CA, United States).

Intracellular NAD⁺/NADH content

Total intracellular NAD⁺/NADH was quantified using the NAD/NADH-Glo Bioluminescent Assay Kit (G9072; Promega Corporation, United States), according to the manufacturer’s instructions. HepaRG cells were first lysed with dodecyl trimethyl ammonium bromide and treated to neutralize their counterparts. To measure NAD⁺, the extract was treated with 25 μL of 0.4 N HCl and heated at 60°C for 15 min, incubated at RT for 10 min, following the addition of 25 μL Trizma base. To quantify NADH, the extract was incubated at 60°C for 15 min followed by further incubation for 10 min at RT. Then 50 mL HCl/Trizma solution was added to the extract. In the presence of each species, a reductase reduced a pro-luciferin reductase to luciferin. The intensity of light (proportional to the amount of each metabolite) was detected by a luminometer (DTX 880 Microplate Reader; Beckman Coulter).

Statistical analyses

Data are expressed as the mean ± standard deviation of three different experiments. Within-group variability was analyzed using Levene’s test for homogeneity of variances. Differences between two groups (P10 vs P20) were determined by the Student’s t-test, while two-way analysis of variance was used to test the main effects of senescence (S, P10 vs P20) or transdifferentiation (T, Basal vs DMSO) as between-subject factors; the interaction S × T was studied, and a Tukey’s test was used as a post hoc test for multiple comparisons. Statistical significance was accepted when P < 0.05. GraphPad Prism 6.0 software was used to perform the analyses.

RESULTS

Replication induces the expression of senescence markers in HepaRG cells

Replication-induced senescence was first described in human fibroblasts because of serial culture passages[20]. Since increased SA-β-gal activity is a well-recognized marker of senescence[21], positive SA-β-gal-stained cells were counted to validate senescence in HepaRG passaged up to P10 compared to cells passaged up to P20. As shown in Figure 1A and B, replication-induced senescence in P20 HepaRG cells resulted in increased numbers of cells positive for SA-β-gal staining compared with P10 cells. A further marker of senescence is represented by the extent of DNA damage[22], determined by the comet assay, which measures the prevalence of single- and double-strand breaks. The level of DNA damage was calculated from the percentage of DNA in the tail of the comet formed after single-cell gel electrophoresis, as broken DNA moves faster in the current. The percentage of DNA in the tails vs the core was analyzed using Comet Assay IV software, and the results showed that P20 HepaRG cells presented with higher levels of damaged DNA than P10 cells (Figure 1C and D).

Figure 1 Characterization of markers of cellular senescence in HepaRG cells. A: Senescence-associated β-galactosidase staining of HepaRG cells at passage 10 (P10) or P20. After staining, the cells were imaged by phase contrast microscopy. Micrographs are shown at magnification 200 ×; B: Quantitative analysis of positive β-galactosidase-stained cells in P10 and P20 HepaRG cells; C: Representative images of DNA damage analysis by the comet assay in P10 and P20 HepaRG cells. After staining, the cells were imaged by fluorescence microscopy. Micrographs are shown at magnification 600 ×; D: Quantitative analysis of DNA tails from 100 cells. Data in the graphs are represented as the mean ± SD of three independent experiments. Statistical differences were assessed by the Student’s t-test. ^aP < 0.05 vs P10; ^bP < 0.01 vs P10.

Transdifferentiation towards bipotent progenitors is altered in senescent HepaRG cells

HepaRG cells are able to actively proliferate and commit toward hepatocyte and biliary differentiation pathway, reaching maximum cell differentiation after a 2-wk DMSO exposure[11]. Then we exposed P10 and P20 HepaRG cells to DMSO, and studied the expression of markers of both progenitor and differentiated cells.

The expression of adhesion molecules such as CD49a (limited to hepatocyte-like cells) and CD49f (associated with biliary-like cells) after the differentiation process induced by DMSO was lower in P20 than in P10 HepaRG cells. Even CD184/C-X-C motif chemokine receptor 4, a known definitive endoderm marker, was less expressed in P20 than in P10 HepaRG cells exposed to DMSO. Moreover, the differentiation protocol resulted in lower expression of EpCAM and CK19 (markers of hepatic progenitor cells[23]) in P20 than in P10 HepaRG cells (Figures 2 and 3).

Figure 2 Expression of progenitor markers in HepaRG cells after the transdifferentiation process. A: Representative immunofluorescence images of HepaRG cells at passage 10 (P10) or P20. After staining, the cells were imaged by confocal microscopy. Micrographs are shown at magnification 200 ×; B: Quantitative analysis of positively stained cells in P10 and P20 HepaRG cells. Data in the graphs are represented as the mean ± standard deviation of three independent experiments. Statistical differences were assessed by the Student’s t-test. ^aP < 0.05 vs P10; ^cP < 0.001 vs P10.

Figure 3 Expression of progenitor markers in HepaRG cells before and after the transdifferentiation process. Flow cytometry histograms of HepaRG cells at passage 10 (P10) or P20 in basal conditions or after the transdifferentiation protocol (dimethyl sulfoxide). After staining, the cells were analyzed with flow cytometry.

The expression of genes typical of a differentiated status, such as albumin, CYP3A4 and γ-GT, and of CEA as marker of undifferentiation was next determined in HepaRG cells before and after DMSO exposure. The results obtained by statistical analyses aimed to test the main effects of S, of T, and the interaction S × T. The main effect of S was significant for the expression of all the genes studied (albumin: F_(1,8) = 44.54, P = 0.0002; CYP3A4: F_(1,8) = 24.22, P = 0.0012; γ-GT: F_(1,8) = 46.82, P = 0.0001; CEA: F_(1,8) = 12.24, P = 0.0081). The main effect of T was significant for the expression of albumin (F_(1,8) = 243.3, P < 0.0001), CYP3A4 (F_(1,8) = 37.50, P = 0.0003), and γ-GT (F_(1,8) = 190.2, P < 0.0001), but not for CEA. The interaction between S and T was significant for the expression of all the genes studied (albumin: F_(1,8) = 43.95, P = 0.0002; CYP3A4: F_(1,8) = 24.22, P = 0.0012; γ-GT: F_(1,8) = 60.71, P < 0.0001; CEA: F_(1,8) = 10.13, P = 0.0129). Post hoc analysis showed that the expression of albumin, CYP3A4 and γ-GT was induced by DMSO exposure both in P10 and in P20 HepaRG cells, supporting the transdifferentiation process; however, mRNA levels were lower in DMSO-treated P20 rather than P10 HepaRG cells (Figure 4A-C). Finally, the expression of CEA was reduced by DMSO exposure in P10 but not in P20 HepaRG cells (Figure 4D).

Figure 4 Gene expression analysis in HepaRG cells at passage 10 or passage 20 in basal conditions or after the transdifferentiation protocol (dimethyl sulfoxide). A: mRNA level of albumin; B: Cytochrome P350 3A4 (CYP3A4); C: Gamma-glutamyl transpeptidase (γ-GT); D: Carcinoembryonic antigen (CEA) measured by reverse transcription-polymerase chain reaction. Data in the graphs are represented as the mean ± standard deviation of three independent experiments. Statistical differences were assessed by two-way analysis of variance followed by the Tukey test as the post hoc test. ^aP < 0.05; ^bP < 0.01; ^cP < 0.001. P10: Passage 10; P20: Passage 20; DMSO: Dimethyl sulfoxide.

The same analysis was performed on specific functional activities of HepaRG cells, studying albumin secretion, CYP3A4 and γ-GT activity, and the results are represented in Figure 5. In particular, we found a main effect of S (albumin secretion: F_(1,8) = 20.95, P = 0.0018; CYP3A4 activity: F_(1,8) = 10.24, P = 0.0126; γ-GT activity: F_(1,8) = 18.13, P = 0.0028), T (albumin secretion: F_(1,8) = 120.2, P < 0.0001, CYP3A4 activity: F_(1,8) = 9.15, P = 0.0164; γ-GT activity: F_(1,8) = 68.25, P < 0.0001), and interaction between S and T (albumin secretion: F_(1,8) = 20.53, P = 0.0019; CYP3A4 activity: F_(1,8) = 7.181, P = 0.0233; γ-GT activity: F_(1,8) = 19.36, P = 0.0023). Post hoc analysis showed that albumin secretion, CYP3A4 and γ-GT activities were reduced in DMSO-treated P20 rather than P10 HepaRG cells.

Figure 5 Functional tests in HepaRG cells at passage 10 or passage 20) in basal conditions or after the transdifferentiation protocol (dimethyl sulfoxide). A: Albumin secretion; B: Cytochrome P350 3A4 (CYP3A4); C: Gamma-glutamyl transpeptidase (γ-GT) activities. Data in the graphs are represented as the mean ± standard deviation of three independent experiments. Statistical differences were assessed by two-way analysis of variance followed by the Tukey test as the post hoc test. ^aP < 0.05; ^bP < 0.01; ^cP < 0.001. P10: Passage 10; P20: Passage 20; DMSO: Dimethyl sulfoxide.

Taken together, these results suggest that the transdifferentiation process triggered by DMSO in HepaRG cells is hampered by the replication-induced senescence.

Senescence-induced impairment of HepaRG transdifferentiation is associated with mitochondrial dysfunction

Mitochondrial dysfunction is recognized as one of the hallmarks of senescence[24]. Thus, we performed respirometry analyses in HepaRG by high-resolution oximetry. Figure 6A details the protocol used and is representative of the oxygraphic profiles in young (P10) and senescent (P20) HepaRG cells in basal conditions and after the transdifferentiation protocol (DMSO). The resting respiration (RR), which depends on endogenous substrates, was impacted by both the S and the T factor (S: F_(1,8) = 78.38, P < 0.0001; T: F_(1,8) = 57.12, P < 0.0001), and by the interaction S×T (F_(1,8) = 60.86, P < 0.0001). Addition of oligomycin (a F_oF₁-ATP synthase inhibitor) reduced the oxygen uptake, suggesting that most of the mitochondrial respiration was coupled to the synthesis of ATP. Nevertheless, both the S and the T factor impacted this parameter (S: F_(1,8) = 10.83, P = 0.0018; T: F_(1,8) = 6.847, P = 0.0308), as well as the interaction S × T (F_(1,8) = 8.104, P = 0.0216). Restoration of the oxygen uptake by the addition of valinomycin (a K⁺ ionophore which uncouples oxygen consumption from ATP synthesis) was also impacted by both the S and the T factor (S: F_(1,8) = 59.58, P < 0.0001; T: F_(1,8) = 53.09, P < 0.0001), and the interaction S × T (F_(1,8) = 55.97, P < 0.0001). The post hoc analysis showed that the oxygen uptake in all the examined conditions was higher in P10 HepaRG cells after the transdifferentiation protocol with respect to the other samples (Figure 6B). The effects of both S and T, and their interaction, was observed on the ATP-dependent oxygen uptake, calculated as the difference between RR and oligomycin-induced respiration (S: F_(1,8) = 28.94, P = 0.0007; T: F_(1,8) = 25.35, P = 0.001; S × T: F_(1,8) = 27.36, P = 0.0008). The post hoc analysis showed that the P10 HepaRG cells after the transdifferentiation protocol exhibited a higher ATP-dependent oxygen uptake than the other samples (Figure 6C). The respiratory control ratio (RCR), calculated as the ratio between RR and oligomycin-induced respiration, was also influenced by both S and T (S: F_(1,8) = 14.76, P = 0.0049; T: F_(1,8) = 13.43, P = 0.0064), and their interaction (S × T: F_(1,8) = 11.38, P = 0.0097). The post hoc analysis resulted in a higher RCR for the P10 HepaRG cells after the transdifferentiation protocol compared to the other samples (Figure 6D). Taken together, these data suggest that the transdifferentiation induced by DMSO increases mitochondrial respiration in HepaRG cells; however, this does not occur in senescent cells.

Figure 6 Measurement of mitochondrial respiration in HepaRG cells at passage 10 or passage 20 in basal conditions or after the transdifferentiation protocol (dimethyl sulfoxide). A: Representative oxymetric traces of mitochondrial respiration in HepaRG cells. Where indicated, the following were added: 4 × 10⁶ HepaRG cells, 8 μg/mL oligomycin (OL), 2 μg/mL valinomycin (VAL), 3 mmol/L potassium cyanide (KCN). The continuous line represents the oxygen concentration measured every 0.1 s throughout the time-course of the assay; B: Representative graph of the normalized and KCN-insensitive-corrected oxygen consumption rates measured under resting respiration (RR) conditions, in the presence of OL and in the presence of VAL (see panel A); C: ATP-dependent oxygen consumption measured as absolute difference between that obtained in the absence and that in the presence of oligomycin (RR-OL); D: Respiratory control ratios obtained dividing the oxygen consumption rates measured under resting conditions by that in the presence of oligomycin (RR/OL). Data in the graphs are represented as the mean ± standard deviation of three independent experiments. Statistical differences were assessed by two-way analysis of variance followed by the Tukey’s test as the post hoc test. ^aP < 0.05; ^bP < 0.01; ^cP < 0.001. DMSO: Dimethyl sulfoxide; P10: Passage 10; P20: Passage 20.

To confirm a defect in mitochondrial dysfunction of senescent HepaRG, the cellular ATP concentration was further measured. We observed a significant effect of both S and T (S:F_(1,8) = 14.89, P = 0.0048; T:F_(1,8) = 12.15, P = 0.0083), and the post hoc analysis showed that the ATP content of transdifferentiated P10 HepaRG cells was higher than the other samples (Figure 7A). SA mitochondrial dysfunction can be consequent to an exhaustion of the oxidized form of NAD⁺[25]. The analysis of NAD⁺/NADH content in HepaRG cells revealed no changes in NADH; nevertheless, the impact of S, T, and their interaction was observed on both NAD⁺ (S:F_(1,8) = 13.84, P = 0.0059; T:F_(1,8) = 26.00, P = 0.0009; S × T: F_(1,8) = 10.47, P = 0.0120) and NAD⁺/NADH (S:F_(1,8) = 22.23, P = 0.0015; T:F_(1,8) = 19.75, P = 0.0022; S × T: F_(1,8) = 10.02, P = 0.0133). The post hoc analysis resulted in higher NAD⁺ and NAD⁺/NADH in transdifferentiated P10 HepaRG cells with respect to the other samples (Figure 7B).

Figure 7 ATP and nicotinamide adenine dinucleotide (NAD)/NAD with hydrogen content in HepaRG cells at passage 10 or passage 20 in basal conditions or after the transdifferentiation protocol (dimethyl sulfoxide). A: Cellular ATP content was expressed as ATP concentration/cells; B: Intracellular nicotinamide adenine dinucleotide (NAD⁺) and NAD with hydrogen (NADH) were extracted, and NAD⁺ and NADH levels were measured using a microplate reader. NAD⁺/NADH ratio was calculated based on the concentration of NAD⁺ and NADH. Data in the graphs are represented as the mean ± standard deviation of three independent experiments. Statistical differences were assessed by two-way analysis of variance followed by the Tukey’s test as the post hoc test. ^aP < 0.05; ^bP < 0.01. DMSO: Dimethyl sulfoxide; P10: Passage 10; P20: Passage 20.

DISCUSSION

This study demonstrates that cellular replication in HepaRG cells is associated with both the expression of senescence markers and the reduction in transdifferentiation potential. Indeed, data related to phenotype, gene expression, functional analysis, and mitochondrial function suggest that the transdifferentiation process of HepaRG is preserved after a few passages, but it is altered by continuous replication.

HepaRG represents a unique cell line characterized by high expression of detoxifying and metabolizing enzymes, transport proteins, and nuclear receptors, as well as the ability to transdifferentiate toward hepatocyte-like and cholangiocyte-like cells[9-11]. Other than representing a viable tool for cell biology, drug metabolism, and virology studies, differentiated HepaRG cells are suitable to generate humanized liver in rodent models, allowing in vivo studies of liver development and physiology[26,27].

Cultured cells stop proliferating after a finite number of divisions, a phenomenon defined to as replicative senescence[28,29]. However, cells can escape replicative senescence by several mechanisms, including overexpression of viral genes such as simian virus 40 large T antigen (which inactivates p53), or telomerase reverse transcriptase protein (TERT, which elongates telomeres)[30,31]. Proliferating HepaRG cells show over-expression and hyper-activation of human TERT, as well as inhibition of p21 and p27, which in turn inhibit the catalytic activity of cyclin-dependent kinases and stop the cell cycle[32]. Nevertheless, our data demonstrate that a replication protocol induces the expression of senescence markers in HepaRG cells, such as positive SA-β-gal staining and DNA alterations. Thus, we hypothesize that the mechanisms that induce senescence in HepaRG cells are different from telomere shortening and p21/p27 activation.

Senescence impairs the transdifferentiation of several cell lines[33-36]. Senescent cells secrete a large variety of molecules that change the surrounding microenvironment, with consequent alterations of differentiation and tissue regeneration[37,38]. These compounds include a wide range of cytokines, growth factors, and signaling molecules that are included in the SA secretory phenotype (SASP)[39]. Our data clearly demonstrate that replicative senescence alters the transdifferentiation process of HepaRG cells. Nevertheless, we did not analyze the SASP in our study, since this marker is strictly linked with telomere shortening[40]. By contrast, we focused on mitochondria, since the homeostasis of these organelles is crucial for several aspects of senescence including SASP[41]. Indeed, the impairment of mitochondrial oxidative phosphorylation (OxPhos) is mainly involved in the early steps of cell senescence[42]. Moreover, senescent cells exhibit severe metabolic alterations associated with mitochondrial metabolites, such as oxidized to reduced NAD ratios (NAD+/NADH)[43]. The novel and straight findings of this study demonstrate the higher mitochondria-related respiration in transdifferentiated HepaRG cells, which was not observed after induction of replicative senescence. Indeed, this respiration was suppressed by oligomycin, a specific inhibitor of ATP synthase, suggesting that transdifferentiated HepaRG cells in resting conditions oxidize substrates to synthetize ATP through the OxPhos system. On the other hand, senescent HepaRG cells subjected to the transdifferentiation protocol exhibited lower oxygen consumption which was insensitive to oligomycin. These results are consistent with a transition from a glycolytic to an oxidative metabolism in transdifferentiated HepaRG cells, which did not occur in senescent cells. Of interest, the RR in transdifferentiated HepaRG cells was lower than that reached in the presence of an uncoupler (maximal respiratory rate), indicating the existence of a respiratory reserve. By contrast, the RR of senescent HepaRG cells undergoing transdifferentiation was similar to the maximal respiratory rate, suggesting the absence of a respiratory reserve. Mitochondrial dysfunction in senescent HepaRG cells after transdifferentiation was further confirmed by a reduction in cellular ATP concentration.

The concentration of the oxidized form of NAD⁺ is determinant for mitochondrial homeostasis[44,45]. Moreover, NAD⁺ is crucial in several metabolic pathways, including glycolysis, citric acid cycle, and OxPhos, with consequent implications for both stemness/differentiation and cell senescence[46]. Cells that maintain a physiological quiescent state to preserve long-term self-renewal capacity are characterized by low NAD⁺ levels to obtain energy from glycolysis; on the contrary, cell differentiation leads to reduced glycolysis and increased OhPhos, which require high NAD⁺ levels[47]. Age-related reductions of both NAD⁺ levels and NAD⁺/NADH ratio are evolutionarily preserved, and consistent evidence for low NAD⁺ has been provided for several old mammalian tissues[48]. Our data show that HepaRG transdifferentiation is associated with increased NAD⁺ and relatively stable NADH, with consequent high NAD⁺/NADH. However, the raise of both NAD⁺ and NAD⁺/NADH is not observed in HepaRG undergoing replicative senescence. Even though the effect of a treatment goes beyond the scope of our study, it is conceivable that a replenishment of NAD⁺ would protect mitochondria and improve the transdifferentiation process of senescent HepaRG cells, as already described for other cell types[49-51].

CONCLUSION

In conclusion, the present report demonstrates that HepaRG cells undergo replicative senescence, which is associated with impairment in transdifferentiation, mitochondrial dysfunction, and NAD⁺ depletion. Further investigations are required to refine the molecular mechanism underlying such observation. The limitations in the transdifferentiation potential of senescent HepaRG cells, with consequent alteration of metabolic and regenerative properties, may have serious implications when this cell line is applied for basic studies.

ARTICLE HIGHLIGHTS

Research background

The HepaRG cell line is used to study metabolism, toxicology, and the regeneration /differentiation processes, as a replacement to primary hepatocytes, HepG2, and Huh-7 cells. These cells exhibit a hepatocyte-like morphology and express hepatocyte-specific functions in defined culture conditions; furthermore, HepaRG display features of human oval ductular bipotent hepatic progenitors.

Research motivation

Cellular senescence consists in a steady cell cycle block occurring because of different harmful events, leading to defective stemness and differentiation processes, as well as changes in cell cycle regulation, signal transduction, and metabolism. The impact of senescence on HepaRG cells has not yet been investigated.

Research objectives

This study investigated whether a replication protocol would induce senescence in HepaRG cells. In addition, we characterized the effects of senescence on transdifferentiation capacity and mitochondrial metabolism.

Research methods

The transdifferentiation capacity of HepaRG cells over passage 10 (P10) vs passage 20 (P20) was compared. To stimulate transdifferentiation, HepaRG cells were treated with dimethyl sulfoxide (DMSO). Aging was evaluated by senescence-associated (SA) β-galactosidase activity and the comet assay. HepaRG transdifferentiation was analyzed by confocal microscopy and flow cytometry (expression of cluster of differentiation 49a [CD49a], CD49f, CD184, epithelial cell adhesion molecule [EpCAM], and cytokeratin 19 [CK19]), by quantitative PCR analysis (expression of albumin, cytochrome P450 3A4 [CYP3A4], γ-glutamyl transpeptidase [γ-GT] and carcinoembryonic antigen [CEA]) and functional analysis (albumin secretion, CYP3A4 and γ-GT). Mitochondrial respiration, the ATP and the NAD⁺/NADH content were also measured.

Research results

We first observed that replication induces the expression of senescence markers in HepaRG cells, since SA β-galactosidase staining was higher in P20 than in P10 HepaRG cells, and the comet assay showed a consistent DNA damage in P20 HepaRG cells. We further reported that transdifferentiation towards bipotent progenitors is altered in senescent HepaRG cells, as P20 HepaRG cells exhibited a reduction of CD49a, CD49f, CD184, EpCAM and CK19 – with respect to P10 – after DMSO treatment. Furthermore, the lower gene expression of albumin, CYP3A4, and γ-GT, as well as the reduced albumin secretion capacity, CYP3A4, and γ-GT activity were reported in transdifferentiated P20 compared to P10 cells. By contrast, the gene expression level of CEA was not reduced by transdifferentiation in P20 cells. Finally, we show that senescence-induced impairment of HepaRG transdifferentiation is associated with mitochondrial dysfunction, since both cellular and mitochondrial oxygen consumption were lower in P20 than in P10 transdifferentiated cells, and both ATP and NAD⁺/NADH were depleted in P20 cells with respect to P10 cells.

Research conclusions

The present study demonstrates that HepaRG cells undergo replicative senescence, with consequent impairment in transdifferentiation, functional activity, mitochondrial dysfunction, and NAD⁺ depletion.

Research perspectives

Further studies will define the molecular mechanisms underlying our observations. The limitations in the transdifferentiation potential of senescent HepaRG cells, with consequent alteration of metabolic and regenerative properties, may have serious implications when this cell line is applied for basic studies.