Spatial modeling reveals nuclear phosphorylation and subcellular shuttling of YAP upon drug-induced liver injury
Abstract
The Hippo signaling pathway controls cell proliferation and tissue regeneration via its transcriptional effectors yes-associated protein (YAP) and transcriptional coactivator with PDZ-binding motif (TAZ). The canonical pathway topology is characterized by sequential phosphorylation of kinases in the cytoplasm that defines the subcellular localization of YAP and TAZ. However, the molecular mechanisms controlling the nuclear/cytoplasmic shuttling dynamics of both factors under physiological and tissue-damaging conditions are poorly understood. By implementing experimental in vitro data, partial differential equation modeling, as well as automated image analysis, we demonstrate that nuclear phosphorylation contributes to differences between YAP and TAZ localization in the nucleus and cytoplasm. Treatment of hepatocyte-derived cells with hepatotoxic acetaminophen (APAP) induces a biphasic protein phosphorylation eventually leading to nuclear protein enrichment of YAP but not TAZ. APAP-dependent regulation of nuclear/cytoplasmic YAP shuttling is not an unspecific cellular response but relies on the sequential induction of reactive oxygen species (ROS), RAC-alpha serine/threonine-protein kinase (AKT, synonym: protein kinase B), as well as elevated nuclear interaction between YAP and AKT. Mouse experiments confirm this sequence of events illustrated by the expression of ROS-, AKT-, and YAP-specific gene signatures upon APAP administration. In summary, our data illustrate the importance of nuclear processes in the regulation of Hippo pathway activity. YAP and TAZ exhibit different shuttling dynamics, which explains distinct cellular responses of both factors under physiological and tissue-damaging conditions.
Editor's evaluation
The Hippo signaling pathway is essential for multiple physiological processes, most notably the regulation of cell proliferation and survival during wound healing. Wehling et al. provide a molecular framework for an alternative mechanism by which the Hippo effector molecule YAP's sub-cellular localization is regulated by cell compartment-specific phosphorylation. Specifically, the authors demonstrate dynamic regulation of shuttling of YAP both in vitro and in vivo during drug-induced liver injury. Given the importance and developmental conserveness of the Hippo pathway, the work is of broad interest to the field of developmental and regenerative biology.
https://doi.org/10.7554/eLife.78540.sa0Introduction
The evolutionary conserved Hippo signaling pathway controls tissue homeostasis through the regulation of cell proliferation, apoptosis, differentiation, and cellular fate (Zhao et al., 2011). Activity of this pathway is affected by information derived from extracellular matrix stiffness, actomyosin dynamics, and cell density (Deng et al., 2015; Zhao et al., 2007), which in turn modulate a serine/threonine kinase cassette consisting of serine/threonine kinase 3/4 (STK3/4; MST1/2) and large tumor suppressor 1/2 (LATS1/2). According to the canonical Hippo pathway model, active cytoplasmic LATS1/2 phosphorylate and inactivate two important downstream pathway effectors: the transcriptional co-activators yes-associated protein (YAP) and its paralog transcriptional coactivator with PDZ-binding motif (TAZ). This phosphorylation is associated with nuclear YAP/TAZ exclusion followed by their proteasomal degradation. In contrast, Hippo pathway inactivation in the cytoplasmic compartment causes YAP/TAZ dephosphorylation and nuclear translocation. In the nucleus, YAP/TAZ bind DNA sequence-specific transcription factors such as TEA DNA-binding proteins (TEADs) or forkhead box protein M1 (FOXM1), which control the transcription of genes involved in for example, cell cycle control and paracrine communication (Marti et al., 2015; Thomann et al., 2020; Weiler et al., 2017).
Due to its pivotal role in the regulation of cell proliferation, the Hippo/YAP/TAZ axis is important for tissues maintenance during regenerative processes. Indeed, YAP- or YAP/TAZ deficiency reduce the regenerative capacity of tissues such as liver, skin, and heart (Lee et al., 2014; Lu et al., 2018; Xin et al., 2013). As exemplified in detail for the liver, YAP is activated in hepatocytes under disease conditions that support a continuous regenerative response in different in vivo model systems (Machado et al., 2015; Mooring et al., 2020). Interestingly, the roles of YAP and TAZ are not identical, since TAZ, but not YAP, contributes to fat accumulation in the liver (steatosis) (Wang et al., 2016). In contrast, YAP protects from liver ischemia–reperfusion injury by promoting tissue repair (Liu et al., 2019). These findings strongly suggest a cell-protective role of the Hippo pathway under liver injury conditions; however, the function of YAP and TAZ in this biological process may differ. As YAP and TAZ control the expression of similar target genes in different cell types (Weiler et al., 2020; Zanconato et al., 2015), other regulatory mechanisms must account for these phenotypic differences such as differential nuclear/cytoplasmic shuttling (Reggiani et al., 2021). However, detailed comparative studies regarding subcellular localization dynamics for YAP and TAZ are missing as well as spatially resolved computational framework such as partial differential equation (PDE) models.
Recent data demonstrated a direct impact on YAP by acetaminophen (APAP), which is a widely used analgesic drug and leading cause of acute liver failure (Poudel et al., 2021). Here, APAP overdose in mice caused liver injury, which was associated with a prominent nuclear YAP enrichment in hepatocytes. Interestingly, silencing of YAP did not impair liver regeneration (which would be expected after inactivation of this proproliferative factor), but promoted tissue regeneration upon APAP-induced liver damage (Poudel et al., 2021). In contrast, a supportive role of YAP in tissue repair and regeneration was suggested by chemical blockade of upstream Hippo pathway kinases MST1/2, which reverted APAP-induced liver damage (Fan et al., 2016). These results point to intricate and context-specific roles of YAP and/or TAZ under regenerative conditions.
Our study comprises a quantitative and mechanistic analysis of YAP/TAZ dynamics based on spatially resolved, high-throughput confocal live cell imaging data and PDE modeling. By using this experimental and computational toolbox, we show that nuclear phosphorylation of YAP/TAZ is important for their subcellular shuttling dynamics and that both factors do not equally respond to cell density and APAP-induced hepatocellular damage. We mechanistically show that APAP overdose stimulates a rapid reactive oxygen species (ROS) response, which facilitates the physical nuclear interaction of YAP and AKT as predicted by computational modeling. Thus, our data demonstrate for the first time that APAP-induced YAP activation is a dynamic and specific process that relies on the sequential activation of molecular events.
Results
Establishment of an in vitro model for measuring time-resolved spatial localization of YAP and TAZ in hepatocellular cells
Due to technical limitations, primary human hepatocytes are not suitable for the analysis of dynamic YAP/TAZ shutting in vitro: they rapidly undergo trans-differentiation in culture and different human donors exhibit genomic/epigenetic variability, which diminish reproducibility. We therefore selected the hepatocyte-derived liver cancer cell line Hep3B that was characterized by a prominent YAP/TAZ and cytochrome P450 2E1 (CYP2E1) expression (Figure 1—figure supplement 1A; Weiler et al., 2020). Detectable CYP2E1 levels are required for the enzymatic transformation of APAP to toxic N-acetyl-p-benzoquinone imine (NAPQI) (Lee et al., 1996).
Since we were interested in cell-to-cell variability and spatial localization information of YAP and TAZ, we established a cell line that allowed the quantitative and spatiotemporal investigation of subcellular YAP/TAZ distribution near the single-cell resolution. For this, we generated Hep3B cells that stably express mVenus-tagged YAP (mVenus-YAP) and mCherry-tagged TAZ (mCherry-TAZ) (Figure 1A). In addition, mCerulean-tagged histone H2B (mCerulean-H2B) allowed the spatial allocation of YAP/TAZ. These chosen fluorophores facilitated optimal discrimination between emission spectra by live cell microscopy (Figure 1—figure supplement 1B). Subcellular localization changes of tagged YAP and TAZ under variable cell density conditions illustrated functionality of both proteins regarding nuclear–cytoplasmic shuttling (Figure 1B). However, we observed a high degree of intracellular variability.

YAP and TAZ show distinct nuclear shuttling in hepatocellular cells upon increasing cell density.
(A) The hepatocyte-derived cell line Hep3B was transduced using three lentiviral vectors coding for mCerulean-tagged H2B (pLentiPGK-mCerulean-H2B), mVenus-tagged YAP (pRRLN-EF1α-mVenus-YAP), and mCherry-tagged TAZ (pRRLN-EF1α-d2mCherry-TAZ). Combined treatment with the antibiotics hygromycin, geneticin, and blasticidin selected for triple-positive cells. Cells were confocally imaged in three channels: 445 nm (CFP, mCerulean), 514 nm (YFP, mVenus), and 561 nm (RFP, mCherry). Scale bar: 250 µm. (B) Exemplary pictures illustrating the subcellular localization of H2B, YAP, and TAZ proteins under low (left) and high (right) cell density conditions. Scale bar: 50 µm. (C) Automatic image analysis workflow depicts analysis of nuclear (left) and cytoplasmic (right) fluorescence intensity in confocal images of living cells. The acquired images were prescreened for imaging artifacts (e.g., out-of-focus images were excluded) and subjected to the image analysis pipeline in Fiji. Object classification (Weka algorithm), thresholding, and object counting were performed. Object masks were overlaid with original data and the nuclear/cytoplasmic ratio (NCR) was calculated by dividing average fluorescent signal intensity of YAP/TAZ in the nucleus with the fluorescence intensity in the cytoplasm. Scale bar: 50 µm. (D) Quantification of YAP and TAZ NCR under increasing cell density conditions. One of four representative experiment is depicted. Left: the NCR of YAP (yellow, n = 310) and TAZ (red, n = 310) was plotted against cell count per visual field (0.4 mm2). A high NCR value indicates predominant nuclear localization. Dashed line shows mean cell count. Right: violin plots summarize shift in NCR between low (below mean cell count) and high (above mean cell count) cell density for YAP and TAZ.
For analysis of variable high-throughput data, we developed an algorithm for the semiquantitative measurement of tagged YAP and TAZ in nuclei and cytoplasm. In brief, spatial image data were classified based on the expression of mVenus-YAP and mCherry-TAZ using the Weka segmentation algorithm (Figure 1C, Figure 1—figure supplement 1C; Arganda-Carreras et al., 2017). The fluorophore intensity measurements were used to calculate the nuclear/cytoplasmic ratio (NCR), which characterized the relative nuclear enrichment of YAP or TAZ. Therefore, this approach allowed us to simultaneously define slightest subcellular changes of YAP and TAZ in small cell populations (field of view, FOV) as well as individual cells under live cell conditions.
Subsequently, the fluorophore-expressing cells were grown under variable cell density conditions and the YAP/TAZ localization was measured by confocal microscopy. The results confirmed that tagged YAP and TAZ dynamically shuttled between cytoplasm and the nucleus depending on abundance of cell–cell contacts (cell density) (Figure 1D). However, TAZ showed a considerably less pronounced dynamic behavior compared to YAP (Figure 1D). YAP was predominantly detectable in cell nuclei under low cell density conditions (NCR > 1), while increasing cell density led to cytoplasmic enrichment of both factors (NCR < 1). Using this algorithm, cell density-dependent protein shuttling was also demonstrated for other liver cancer cells (Tóth et al., 2021).
Together, the confocal imaging and image processing pipeline efficiently detects subtle differences in the dynamic subcellular changes of YAP and TAZ.
Mathematical modeling predicts nuclear phosphorylation of YAP and TAZ
To investigate why Hippo pathway effectors YAP and TAZ differently respond to a low cell density and which Hippo pathway topology can explain the observed localization differences, we mathematically modeled the Hippo pathway using a PDE modeling framework (for detailed description of the PDE modeling approach refer to Materials and methods and Appendix 1).
First, we investigated a canonical Hippo pathway model, where phosphorylation takes place in the cytoplasm and where unphosphorylated YAP/TAZ shuttle to the nucleus (Figure 2A). However, this mathematical model was not able to reproduce the localization patterns which were observed in the experimental data (Figure 2B). For example, the model cannot sufficiently explain the distribution of YAP and TAZ in cell nuclei or at the nuclear membrane. In detail, the cell area around the nuclear envelope on the cytoplasmic side strongly underrepresented the concentration of YAP/TAZ, as indicated by the blue pixels in the residual image (experimental data minus model simulation). Moreover, the simulated YAP/TAZ distribution pattern of the canonical model showed a decreased protein concentration in the center of the nucleus, with increasing gradient toward the nuclear envelope. This indicates that the experimental observation is dominated by a protein that is disseminated from the nucleus and undergoes diffusion and degradation. The canonical Hippo pathway cannot explain this effect, illustrated by residuals (Figure 2B).

Mathematical modeling predicts that nuclear phosphorylation controls YAP/TAZ subcellular localization.
(A) The model reaction scheme of the canonical Hippo pathway. Only unphosphorylated YAP/TAZ can enter the nucleus, whereas pYAP and pTAZ are exclusively localized to the cytosol. The phosphorylation of YAP/TAZ takes place exclusively outside nuclei. (B) Partial differential equation (PDE) model simulation of the canonical Hippo pathway compared to the experimentally measured YAP/TAZ localization. The residuals (experimental data minus model simulation) indicate low spatial accordance of the model simulation with the experimentally measured subcellular localization of YAP/TAZ. Images were normalized to the maximal value within each image. Nuclear outline is indicated in black. (C) The model reaction scheme of the alternative Hippo pathway model. Phosphorylation and dephosphorylation of YAP/TAZ take place in the nucleus. Unphosphorylated YAP/TAZ is imported in the nucleus and phosphorylated YAP/TAZ is rapidly exported to the cytoplasm. (D) PDE simulation of the alternative model compared to experimental data. Residuals for the subcellular localization (e.g., nuclear distribution) sufficiently reflect results from confocal microscopy. Nuclear outline is indicated in black. (E) Simulated impact of the nuclear phosphorylation to dephosphorylation ratio on subcellular localization of YAP in the alternative Hippo pathway model. Left: model simulation of two phosphorylation/dephosphorylation rates (0.17 and 0.75). Right: summarized residuals with respect to experimental data of YAP and TAZ as a function of phosphorylation to dephosphorylation ratio in the nucleus. (F) Western immunoblot after nuclear and cytoplasmic protein fractionation under low cell density conditions. The central Hippo pathway kinases LATS1/2 were detectable in the nuclear fraction. As expected, high pYAP levels are detectable in the cytoplasm, illustrating that the protein is transported outside the nucleus upon phosphorylation. PARP and Tubulin serve as loading controls for nuclear and cytoplasmic fractions, respectively. Equal amounts of cytoplasmic and nuclear proteins were loaded (n = 4; one representative experiment shown). (G) Proximity ligation assay (PLA) for phosphorylated LATS1/2 (pLATS1/2) and YAP (top row). Red dots indicate physical interaction between pLATS1/2 and YAP. DAPI (4′,6-diamidino-2-phenylindole)-stained nuclei are indicated in cyan. Individual pLATS and YAP antibodies serve as assay controls. Bottom row: computational segmentation of nuclei (empty circles) and dots (black spots). Scale bar: 50 µm. (H) Quantification of the PLA assay. Each dot represents an individual image. Left: quantification of interactions between YAP and pLATS1/2 in the nucleus (n = 18) compared to the negative controls (pLATS1/2 and YAP antibodies alone, n = 16 and n = 17, respectively). Right: the number of PLA dots of pLATS1/2 and YAP interaction in cytoplasm and in nuclei (n = 18). Statistical test: two-tailed paired t-test (p value = 0.01) **p ≤ 0.01.
-
Figure 2—source data 1
Raw data for Western Blot shown in Figure 2F.
- https://cdn.elifesciences.org/articles/78540/elife-78540-fig2-data1-v1.zip
This let us conclude that alternative model topologies must be considered to explain the experimental data. We therefore generated an alternative PDE model, which precisely described YAP/TAZ localization as observed in the image data. In this model, a phosphorylation and dephosphorylation reaction of YAP and TAZ in the nucleus was included (Figure 2C,D, Figure 2—figure supplement 1A, B). Moreover, the alternative model describes that unphosphorylated YAP/TAZ are transported to the nucleus and pYAP/pTAZ are rapidly excluded from the nucleus, which agrees with the general consensus that phosphorylated YAP and TAZ predominantly localize in the cytoplasm (Figure 2—figure supplement 1C). Using well-fitting parametrization of the model, the simulated concentration pattern was dominated by phosphorylated YAP/TAZ disseminating from the nucleus, which now matched the observed experimental data well (Figure 2D).
Since YAP and TAZ are considered biochemically similar molecules, the question arises how the pronounced observed differences in the spatial distribution patterns can be explained. Using the alternative Hippo pathway model, we demonstrated that the different phenotypes could be reproduced with one single PDE model where only the parameters of the YAP/TAZ phosphorylation and dephosphorylation reactions in the nucleus differed, while all other parameters (including nuclear import and export reactions) were kept the same (Figure 2E, Figure 2—figure supplement 1D).
In detail, low nuclear phosphorylation/dephosphorylation ratio (0.17) induced nuclear localization of YAP, while high ratio (0.75) induced nuclear protein exclusion. Upon increased phosphorylation activity (>0.2), the nuclear fraction of TAZ decreased and the summarized residual term for TAZ was reduced. According to residual terms, best accordance between the model and the experimental data was achieved at low phosphorylation/dephosphorylation ratio for YAP and at higher phosphorylation/dephosphorylation ratio for TAZ (Figure 2E, right panel). These findings suggest that differences in the nuclear phosphorylation rates of YAP and TAZ can explain the distinct shuttling patterns of both factors (as observed in Figure 1B). In addition, TAZ may require stronger phosphorylation than YAP to induce its efficient shuttling.
The results of the alternative mathematical model required that major components for YAP/TAZ phosphorylation process were present in the cell nucleus. Indeed, subcellular fractionation experiments illustrated that the Hippo pathway kinases LATS1/2 (total and phosphorylated) were predominantly localized in nuclear protein fraction (Figure 2F). Moreover, we and others reported that nuclear LATS1/2 can control YAP phosphorylation and localization (data not shown) (Ege et al., 2018; Li et al., 2014). The relevance of this nuclear interaction is substantiated by a proximity ligation assay (PLA), which illustrated that pLATS1/2 physically interacted with YAP not only in the cytoplasm but also in the nucleus (Figure 2G). To quantitatively compare the number of YAP/pLATS interactions in cell nuclei and cytoplasm, we utilized a computational algorithm counting signals in nuclear and cytoplasmic area of the cells (Figure 2—figure supplement 2A, B). The signal quantification confirmed that YAP and pLATS interactions were detectable in both compartments with significantly higher interaction count in cell nuclei (Figure 2H).
In summary, our findings suggest that nuclear YAP/TAZ phosphorylation contributes to their inactivation and nuclear exclusion. Moreover, variable phosphorylation rates of YAP and TAZ explain the subcellular shuttling differences between both factors.
APAP regulates YAP protein localization and activity
To test if our findings are of relevance under conditions where YAP and TAZ are actively regulated, we decided to use a drug-induced liver injury (DILI) model. For this, we treated Hep3B cells expressing tagged mVenus-YAP and mCherry-TAZ with APAP (10 nM) (Barbier-Torres et al., 2017), and quantitatively investigated the dynamic shuttling of both factors.
Image quantification revealed that APAP led to a gradual and time-dependent nuclear enrichment of YAP after 48 hr, while TAZ only weakly responded to APAP (Figure 3A,B). This APAP-induced nuclear shuttling effect was clearly detectable under high cell density culture conditions, which is characterized by nuclear YAP/TAZ exclusion. Effects on YAP were less pronounced 24 hr after APAP administration (Figure 3—figure supplement 1A,B). No obvious response was detectable for YAP and TAZ at earlier time points post APAP treatment (data not shown). Severe effects of APAP on cell toxicity and apoptosis in the chosen experimental setup were excluded by measuring cell viability and PARP cleavage (data not shown).

APAP regulates YAP phosphorylation, localization, and expression of its target genes.
(A) Live cell confocal imaging of H2B-mCerulean, YAP-mVenus, and TAZ-mCherry in Hep3B cells. Upper row: control treatment (phosphate-buffered saline, PBS). Lower row: APAP treatment (10 mM) induces nuclear enrichment of YAP, but not TAZ protein after APAP treatment within 48 hr (n = 4; one representative experiment shown). Scale bar: 50 µm. (B) Live cell confocal microscopy image quantification. Nuclear/cytoplasmic ratio (NCR) of YAP and TAZ with APAP (10 mM for 48 hr; yellow and red circles, n = 180 per channel) and control (PBS; gray circles, n = 180 per channel) treatment under increasing cell density conditions (represented as mean cell count per visual field). Each dot represents a single visual field. Dashed line: mean cell density. (C) Western immunoblot of YAP and pYAP after APAP treatment (10 mM) in Hep3B cells (n = 6; one representative experiment shown). GAPDH served as a loading control. (D) Relative expression of the YAP target genes ANKRD1 and CYR61 24 hr after APAP (10 mM) treatment in HLF cells (n = 2; one out of two biological replicates shown). (E) Rescue experiment: relative expression of YAP and its target genes CYR61 and ANKRD1 with and without siRNA-mediated YAP inhibition (for 24 hr) with or without APAP (10 mM) treatment for 24 hr (n = 2; one out of two biological replicates shown).
-
Figure 3—source data 1
Source data for Western Blot shown in Figure 3C.
- https://cdn.elifesciences.org/articles/78540/elife-78540-fig3-data1-v1.zip
The APAP-induced nuclear enrichment of YAP was confirmed with cell population-based Western immunoblotting followed by detection of total and phosphorylated YAP (Figure 3C, for quantification of all western blots see Figure 3—figure supplement 2). Notably, at early time points YAP hyperphosphorylation indicated protein inactivation (up to 3 hr after APAP administration). However, between 24 and 48 hr after APAP treatment, a clear YAP dephosphorylation was detectable (Figure 3C). Similar results were observed for another hepatocyte-derived cell line (Figure 3—figure supplement 1C).
To test if APAP-dependent dephosphorylation/activation of YAP at later time points is leading to its transcriptional activation, we measured the relative expression of known YAP target genes cysteine-rich angiogenic inducer 61 (CYR61) and ankyrin repeat domain 1 (ANKRD1) by real-time PCR. Indeed, we observed induction of these genes already 24 hr after APAP treatment, which reflected the earliest activation of YAP (Figure 3D). To confirm YAP dependence on the target gene induction and to exclude YAP-independent effects on gene expression, we performed additional rescue experiments. For this, the expression of YAP was silenced by RNA interference (RNAi) in APAP-treated cells followed by the measurement of YAP target genes. According to our hypothesis, YAP inhibition partly abolished the APAP-dependent induction of CYR61 and ANKRD1 (Figure 3E).
Together, these data demonstrate that APAP predominantly acts on the Hippo pathway effector YAP in a bimodal manner. Upon APAP treatment an immediate phosphorylation/inactivation of YAP is followed by its late dephosphorylation/activation.
APAP controls YAP phosphorylation via ROS and AKT
Based on previous data, we hypothesized a mechanistic connection between APAP-induced ROS (Barbier-Torres et al., 2017; Shuhendler et al., 2014) and AKT-driven YAP phosphorylation (Basu et al., 2003; Romano et al., 2010). This mechanistic link was investigated at early time points after APAP treatment (up to 3 hr) to exclude unspecific effects caused by APAP at later time points (e.g., due to secondary and/or unspecific APAP effects).
Western immunoblotting showed that both YAP and AKT were phosphorylated early after APAP incubation (Figure 4A), which was in accordance with our previous findings (Figure 3C). In parallel, an ROS detection assay revealed that APAP-induced prominent ROS activity comparable to the known ROS inducers hydrogen peroxide (H2O2) and tert-butyl hydroperoxide (TBHP) in hepatocellular cells (Figure 4B). Moreover, treatment of cells with ROS inducers demonstrated a clear (with H2O2) or moderate (with TBHP) induction of AKT and YAP phosphorylation after 1 and 2 hr, respectively (Figure 4C,D, quantification of western blots in Figure 3—figure supplement 2). These data indicated that APAP regulated ROS activity, which itself controlled AKT and YAP activity.

APAP controls nuclear YAP enrichment via induction of reactive oxygen species (ROS) and AKT.
(A) Western immunoblot analysis of pYAP, YAP, pAKT, and AKT after APAP (10 mM) administration in Hep3B cells after 1 hr (Hep3B). The concentration of pYAP and pAKT is induced after APAP treatment in comparison to phosphate-buffered saline (PBS)-treated cells (ctrl); n = 2. (B) Spectrophotometric ROS measurement after APAP (10 mM), hydrogen peroxide (H2O2, 2 mM), and tert-butyl hydroperoxide (TBHP, 300 µM) treatment in Hep3B cells after 6 hr (n = 8 technical replicates, one out of two biological replicates is shown). APAP as well as H2O2 and TBHP induce ROS formation in living cells. Statistical test: one-way analysis of variance (ANOVA) with Geisser–Greenhouse correction (adjusted p value = 0.001) ***p-value ≤ 0.001. Whiskers depict min and max of the dataset. (C) Western immunoblot analysis of pYAP, YAP, pAKT, and AKT after H2O2 (2 mM) and APAP (10 mM) treatment for 1 hr. APAP, as well as H2O2 induce YAP and AKT phosphorylation compared to untreated cells (ctrl); n = 4. (D) Western immunoblot analysis of pYAP, YAP, pAKT, and AKT dynamics after TBHP (300 µM) and APAP (10 mM) treatment for 2 hr (Hep3B). APAP and TBHP induce AKT and YAP phosphorylation as compared to respective controls (ctrl). (E) Western immunoblot analysis of pAKT, AKT, pYAP, and YAP after hepatocyte growth factor (HGF) treatment (10 ng/µl) for 10 min (Hep3B). The cells were pretreated with an AKT inhibitor VIII (AKTi, 1, 5, and 10 µM) and starved for 3 hr prior HGF administration. Data illustrate that AKT inhibition prevents HGF-induced AKT and YAP phosphorylation. (F) Western immunoblot of pAKT, AKT, pYAP, and YAP after APAP (10 mM) administration for 1 hr (Hep3B). Cells were starved in Fetal calf serum (FCS)-free medium and pretreated with AKTi (5 µM) before APAP treatment for 3 hr. Results demonstrate that AKT inhibition reduces YAP phosphorylation early after APAP treatment. (G) Proximity ligation assay (PLA) of the protein combinations YAP/AKT and pYAP/pAKT. Red dots indicate interactions between YAP/AKT and pYAP/pAKT, respectively. DAPI-stained nuclei are depicted in cyan. Single YAP, AKT, pYAP, and pAKT stains serve as assay controls. Bottom row: segmentation of the nuclei (empty circles) and dots (black spots). One out of two representative experiment is shown. Scale bar: 50 µm. (H) Quantification of the interactions between YAP/AKT and pYAP/pAKT in the nucleus (experiment shown in G). One symbol represents one image (YAP = 17, AKT = 16, YAP-AKT = 22, pYAP = 24, pAKT = 25, pYAP-pAKT = 23). (I) PLA showing YAP–AKT interaction units per nuclei of untreated (ctrl, n = 12) and APAP-treated cells (10 mM, for 2 hr, n = 20). APAP treatment induces nuclear interaction between YAP and AKT. Statistical test: unpaired two-tailed parametric t-test (p value = 0.0155). For A and C–F, actin served as loading control.
-
Figure 4—source data 1
Source data for Western Blots shown in Figure 4A,C-F.
- https://cdn.elifesciences.org/articles/78540/elife-78540-fig4-data1-v1.zip
To confirm the rapid and AKT-dependent phosphorylation of YAP, liver cells were treated with hepatocyte growth factor (HGF), which activates/phosphorylates AKT kinases via binding the receptor c-MET within 10 min (Xiao et al., 2001). HGF administration led to a clear phosphorylation of AKT but not YAP, which might be caused by saturation effects (no further YAP phosphorylation possible under given culture conditions) (Figure 4E). However, simultaneous inhibition of AKT kinase activity by the specific inhibitor AKTi (Adlung et al., 2017; Barnett et al., 2005), not only abolished AKT phosphorylation but also reduced YAP phosphorylation (Figure 4E). As expected, upon AKTi administration concomitantly with APAP, the phosphorylation of both AKT and YAP decreased (Figure 4F).
Since our alternative PDE model predicted that the nuclear YAP phosphorylation is crucial for protein shuttling and its activity (Figure 2C–E), we hypothesized that AKT and YAP not only physically interact in the cytoplasm but also in cell nuclei. Indeed, PLA illustrated that AKT and YAP, as well as the respective phosphorylated isoforms, interacted in the nucleus (Figure 4G,H). Although, interaction between AKT/YAP and pAKT/pYAP was also detectable in the cytoplasm, a prominent protein colocalization took place in the nuclear compartment (Figure 4G). Importantly, image-based quantification revealed that the interaction between nuclear YAP and AKT increased upon 6hr treatment with APAP (Figure 4I), indicating an APAP-dependent shift to nuclear phosphorylation and interaction, as predicted by the PDE model. Indeed, blocking AKT activity by AKTi reduced YAP phosphorylation, moderately elevated nuclear YAP positivity, and increased the physical interaction between AKT and YAP in cell nuclei (Figure 4F and data not shown).
In summary, our data show that APAP affects YAP activity and shuttling behavior by enhancing cellular ROS followed by nuclear AKT/YAP binding and YAP phosphorylation.
Sequential activation of ROS, AKT, and Hippo/YAP in mouse livers after APAP intoxication
The results of our mathematical model (Figure 2) and the in vitro experiments (Figure 3 and Figure 4) demonstrated a distinct sequence of molecular events that control YAP activity after APAP exposure: ROS induction is followed by AKT activation and YAP phosphorylation/inactivation (early events – up to 6 hr). This is followed by phase, where YAP is dephosphorylated and transcriptionally active (late events – 24 to 48 hr).
To confirm this in vivo, mice were injected with a hepatotoxic dose of APAP (300 mg/kg) and liver tissues were collected up to 16 days (in total at nine time points). Subsequently, liver specimens were subjected to expression profiling and results were investigated regarding the presence of gene signatures specific for ROS, AKT, and Hippo/YAP activity (De Marco et al., 2017; Han et al., 2008; Wang et al., 2018). The experimental results showed that gene expression of ROS, AKT, and Hippo/YAP target genes in livers was altered and prominently activated between 6 hr and 2 days after APAP treatment (Figure 5A–C). To compare the temporal dynamic of the gene signatures, z-scores of genes in the signature per time point were summarized and normalized to the number of genes in the signature. As indicated by the results from the cell culture experiments, the expression index illustrated a specific order of signature activation starting with ROS (6–12 hr), followed by AKT (12 hr to 1 day) and Hippo/YAP signature genes (1–2 days) (Figure 5D).

APAP stimulates the sequential activation of reactive oxygen species (ROS), AKT, and YAP in vivo.
Gene expression profiling over time of mouse livers after APAP treatment (A–C) was performed for up to 16 days (3–5 animals/group, 300 mg/kg). Abundance of gene signatures characteristic for the activity of ROS consisting of 23 genes (Han et al., 2008) (A), AKT (B) with 28 genes (De Marco et al., 2017) and YAP (C) consisting of 23 genes (Wang et al., 2018) were analyzed. In A–C, gene expression values are z-score normalized. (D) Summarized gene expression scores (expression index) over time illustrate the timely order of ROS (max values between 6 and 12 hr), AKT (max values between 12 hr and 1 day) and YAP (max values between 1 and 2 days) target gene signatures after APAP treatment. (E) Mouse liver tissue sections stained for YAP (top tow) and total AKT (bottom row) under control condition and at 6 hr as well as 1, 2, 6, and 16 days after APAP treatment (300 mg/kg). Arrows indicate high YAP and AKT positivity in nuclei. Scale bar: 50 µm. (F) Automatic quantification of nuclei with YAP and AKT positivity in mouse liver tissues. Each dot represents the count of stained nuclei in one image (a tile, 1 mm2) for YAP (n = 733) and AKT (n = 520). (G) Quantification of YAP and AKT nuclear positivity 6 hr and 2 days after APAP treatment. One dot represents the average count of positive nuclei in one animal (n = 4).
To confirm the findings from the gene expression analysis, we performed immunohistochemical staining of liver tissues isolated at five time points after APAP treatment. Staining for total YAP illustrated not only a general increase of YAP positivity in the cytoplasm of surviving hepatocytes, but also its prominent nuclear accumulation after 1–2 days (Figure 5E, arrows). Moreover, nuclear AKT enrichment was already detectable 6 hr after APAP treatment. Interestingly, nuclear AKT accumulation in hepatocytes remained high for the rest of the experiment compared to untreated mice (Figure 5E, arrows).
Because of a high mouse-to-mouse variability, we decided to objectify these results. For this, the immunohistochemical staining of YAP and AKT was quantitatively analyzed using machine learning (random forest) and image analysis methods (Figure 5F). For this, 733 (for YAP) and 516 (for AKT) images were unbiasedly selected from all investigated animals and analyzed with an algorithm, which was trained to detect positively stained nuclei outside of the necrotic areas and tissue borders (Figure 5—figure supplement 1; for detailed information see Materials and methods). The quantification results confirmed the observed nuclear YAP positivity 1–2 days after APAP administration, while nuclear AKT was already detectable at 6 hr and maintained during the experiment (Figure 5F). The shift from exclusive nuclear AKT to nuclear YAP was illustrated by direct comparison of the time points 6 hr and 2 days after APAP injection (Figure 5G). The mechanistic connection between AKT activity and YAP induction in murine hepatocytes was confirmed in independent experiments. Here, hydrodynamic gene delivery of myristoylated AKT led to nuclear enrichment of YAP expression in hepatocytes and subsequent induction of typical YAP target genes (data not shown).
To sum up, we confirm the sequential order of APAP-induced events that cause YAP activation in mouse liver tissue.
Discussion
In this study, we aimed to investigate the different properties of the Hippo pathway effectors YAP and TAZ under in vitro and in vivo conditions that resembled the situation in DILI. For this, we developed a cell-based model system that allowed the systematic and comparative tracing of spatial YAP/TAZ dynamics upon stimulation with APAP, which is one of the most widely used analgesic with high potential of liver toxicity. By integrating time-resolved experimental data, computational modeling, image analysis tools as well as confirmatory in vivo results, we could draw several important conclusions on the role of YAP and TAZ under physiological conditions and in liver cells upon DILI.
Several mathematical modeling approaches have been previously applied to explain the dynamic behavior of the Hippo pathway under physiological or pathological conditions in different organisms. For example, computational modeling was used to explain how Hippo signaling cross-talks with other signaling pathways via protein interactions (Labibi et al., 2020; Romano et al., 2014). Other studies exclusively focused on biophysical/biochemical regulation of the Hippo pathway to investigate how YAP/TAZ activity relates to mechanical input information (Ege et al., 2018; Eroumé et al., 2021; Gou et al., 2018; Scott et al., 2021; Sun et al., 2016).
The goal of our approach was to determine whether nuclear phosphorylation of YAP plays a role in Hippo signaling, specifically in governing YAP localization. A model not including nuclear phosphorylation of YAP could not sufficiently explain the experimentally observed localization pattern and was excluded from further analyses. However, the alternative model proposed in this work, can explain the experimental data and therefore establishes a possible mechanism for YAP localization. However, more complex mechanisms regulating subcellular localization and dynamics of YAP and TAZ cannot be excluded. Therefore, our model partly explains how the Hippo pathway is organized, however, it does not aim to comprehensively describe it.
Here, we decided to use PDE modeling, which allowed us to discover spatial aspects of signaling dynamics and to use live cell image data for parametrizing the Hippo signaling steady-state models. By doing so, we compared PDE model simulations of different model topologies with the experimentally acquired data for YAP as well as TAZ and selected a model topology, which best corresponded to data generated in vitro. The selected alternative Hippo pathway model is based on experimental observations and parametrized within the considered boundaries of experimental evidence (see Appendix 1). The obtained model parameters are results of the best model fit obtained, and allow reproducing the described model behavior. However, since the model fitting is based on steady-state data only and concentration observations are in arbitrary units, the parameter values are subject to structural non-identifiability. Nevertheless, this does not affect the qualitative conclusions drawn from the models for the investigated biological process.
The results of our PDE modeling approach strongly suggested that the nuclear phosphorylation reaction of YAP/TAZ was necessary and sufficient for the reproduction of their subcellular distribution patterns in vitro (i.e., the halo-like nuclear distribution of YAP and TAZ). This observation is of importance since the conventional view on the Hippo pathway topology represents sequential phosphorylation steps in the cytoplasm. Our PDE modeling-based finding would therefore extend the current understanding of Hippo pathway regulation since it adds the cell nucleus as a pivotal spatial compartment in the modulation and adjustment of Hippo/YAP/TAZ signaling.
This conclusion is supported by our experimental findings illustrating that LATS1/2 are expressed in cell nuclei and that the interaction between phosphorylated LATS1/2 and YAP as well as between AKT and YAP is also detectable in this compartment. Importantly, previously published experimental and computational approaches support our findings. For example, nuclear LATS1/2 controls YAP phosphorylation in the context of NF2/Merlin deficiency (Li et al., 2014). In addition, the combination of quantitative photobleaching experiments and modeling illustrated that nuclear export is a crucial reaction for the subcellular YAP distribution (Ege et al., 2018). Combining these insights on nuclear export mechanisms with our findings on nuclear YAP phosphorylation strongly suggests that both processes closely cooperate in the efficient regulation of the Hippo pathway effectors. Lastly, our model does not explicitly exclude phosphorylation reaction in the cytoplasm; however, it suggests the necessity to reevaluate the canonical Hippo pathway scheme. The appeal for the reevaluation and extension of the canonical Hippo pathway, for example, by combining models on nuclear phosphorylation and nuclear transport, would further broaden our understanding on the dynamic spatial behavior of complex signaling pathways (Ege et al., 2018; Shreberk-Shaked and Oren, 2019).
The extension of the canonical Hippo pathway model with processes that control active nuclear YAP phosphorylation illustrates that the regulatory YAP/TAZ phosphorylation must be considered as continuum process (Shreberk-Shaked and Oren, 2019). For instance, the existence of nuclear phosphorylation possibilities by LATS-dependent and -independent (e.g., by AKT) mechanisms increases the complexity but also flexibility and redundancy of the signaling pathway (Ege et al., 2018; Gao et al., 2017; Low et al., 2014). These cellular aspects are part of a fine-tuned regulatory network that allows a rapid and adjustable proliferative response of YAP and TAZ under diverse physiological and pathological cellular conditions (e.g., upon induction of tissue regeneration).
To our knowledge, a direct mathematical comparison of YAP and TAZ regarding their spatial distribution has not been performed, yet. Using the PDE model, we comparatively investigated YAP and TAZ distribution and identified kinetic parameters, which discriminate between both species. Specifically, the observed localization differences (nuclear exclusion vs. nuclear localization) between YAP and TAZ might be caused by distinct phosphorylation rates of YAP and TAZ. Nuclear/cytoplasmic shuttling dynamics for YAP and/or TAZ have been investigated previously (Ege et al., 2018; Plouffe et al., 2018). However, our model predicts that the phosphorylation efficiency for TAZ must be higher than for YAP to achieve the observed differential localization of YAP and TAZ.
By utilizing our image analysis approach for the detection of minor changes in the subcellular distribution of YAP and TAZ, we showed a dynamic shuttling response of YAP (and to a much lesser extent of TAZ) upon APAP treatment. Although, APAP-dependent activation of YAP was described in the literature (Poudel et al., 2021), the underlying molecular mechanisms were poorly understood. Our results demonstrated that APAP-induced YAP shuttling is not the result of an unspecific cellular response but is due to a sequential order of molecular events that include ROS induction (Ghallab et al., 2016; Shuhendler et al., 2014), followed by AKT (in)activation (Koundouros NPoulogiannis, 2018) and nuclear YAP (de)phosphorylation. Time-dependent APAP-induced ROS activation was recently supported by a study by Ghallab et al., which showed transient induction of oxidative stress in mice 8 hr after APAP overdose with a peak at 2 hr (Ghallab et al., 2022). Interestingly, we observed a biphasic YAP regulation: YAP phosphorylation/inactivation at early time points (up to 3 hr) and its dephosphorylated/activation at later time points (>24 hr). Especially, the functionally important activation of YAP, which is associated with the induction of cell proliferation, is well reflected by the temporal activation of YAP in murine tissues upon APAP administration after 24–48 hr.
Importantly, several cellular mechanisms might simultaneously contribute to the APAP-dependent activation of YAP. For example, the direct impact of APAP on MST1/2 kinases, which are central Hippo pathway constituents, represents one additional process that may contribute to the observed effects on YAP (Fan et al., 2016). Moreover, the serine-/threonine-kinase AKT could control the Hippo/YAP axis at different levels. For example, a physical interaction of AKT and YAP in the nucleus (shown here) or the cytoplasm might directly control YAP phosphorylation (Basu et al., 2003). Alternatively, AKT can phosphorylate MST2 and therefore indirectly contribute to YAP activity (Romano et al., 2014). Thus, APAP probably controls YAP in a multimodal manner to achieve a cellular response and it is therefore likely that this multimodal process is part of cell-protective mechanism that counteracts the massive loss of hepatocytes upon APAP-induced DILI (Holland et al., 2022; Shuhendler et al., 2014).
In our manuscript, we show that APAP administration induces a chain of molecular events through the APAP/ROS/AKT axis, which is leading to YAP inactivation (early) followed by YAP activation (late). The late induction of YAP activity might be considered as cell-protective cellular response upon long-term tissue damage; however, YAP also acts as a potent oncogene in different cell types, including hepatocytes, and its overexpression is associated with tumor initiation (Dong et al., 2007; Weiler et al., 2017). Thus, it is tempting to speculate that APAP overdose may not only cause DILI but also increases the risk for liver tumor development. Although controversially discussed in the literature, a recent evaluation of 139 published epidemiologic studies strongly argues against a relationship between APAP exposure and cancer (Weinstein et al., 2021). Actually, APAP has been discussed as therapeutic option for patients with YAP-induced cancer (Poudel et al., 2021). Although, our and other studies clearly demonstrate a mechanistic link between APAP uptake and Hippo/YAP pathway activity, these partly controversial findings illustrate the necessity to decipher this connection under distinct disease conditions in the future.
Materials and methods
Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
---|---|---|---|---|
Gene (Homo sapiens) | YAP | Prof. Loewer (TU Darmstadt) | ENSG00000137693 | – |
Gene (Homo sapiens) | TAZ | Cloned | ENSG00000018408 | – |
Cell line (Homo sapiens) | Hep3B | DSMZ | #ACC93 | – |
Cell line (Homo sapiens) | HepG2 | LGC | ATCC-HB-8065 | – |
Cell line (Homo sapiens) | HLF | JCRB | JCRB0405 | – |
Cell line (Homo sapiens) | SNU182 | ATCC | CRL-2235 | – |
Transfected construct | pRRLN-EF1α-mVenus-YAP | Prof. Loewer (TU Darmstadt) | End-to-end sequencing | – |
Transfected construct | pRRLN-EF1α-d2mCherry-TAZ | This paper | End-to-end sequencing | – |
Transfected construct | pLentiPGK-mCerulean-H2B | Addgene | 90234 | – |
Chemical compound, drug | AKT inhibitor VIII | Merck Millipore | 124018 | – |
Sequence-based reagent | siRNA YAP#1 | This paper | CCACCAAGCUAG AUAAAGA-dT-dT | |
Sequence-based reagent | siRNA YAP#2 | This paper | GGUCAGAGAUAC UUCUUAA-dT-dT | |
Sequence-based reagent | control siRNA | This paper | UGGUUUACAUG UCGACUAA | |
Antibody | YAP (rabbit monoclonal) | Cell Signaling Technology | RRID:AB_2650491 | WB (1:1000) |
Antibody | pYAP (rabbit polyclonal) | Cell Signaling Technology | RRID:AB_2218913 | WB (1:400) PLA (1:200) |
Antibody | pAKT (rabbit monoclonal) | Cell Signaling Technology | RRID:AB_2315049 | WB (1:1000) |
Antibody | AKT (rabbit polyclonal) | Cell Signaling Technology | RRID:AB_329827 | WB (1:1000) PLA (1:200) |
Antibody | β-Actin (mouse monoclonal) | MP Biomedicals, Solon | RRID:AB_2335127 | WB (1:10,000) |
Antibody | LATS1 (mouse monoclonal) | Santa Cruz Biotechnology | sc-398560 | WB (1:200) |
Antibody | LATS2 (mouse monoclonal) | Santa Cruz Biotechnology | sc-515579 | WB (1:200) |
Antibody | pLATS1/2 (rabbit monoclonal) | Cell Signaling Technology | RRID:AB_10971635 | WB (1:500) PLA (1:200) |
Antibody | CYP2E1 (rabbit unknown) | Novus Biologicals | RRID:AB_11021447 | WB (1:500) |
Antibody | PARP (rabbit polyclonal) | Cell Signaling Technology | RRID:AB_2160739 | WB (1:10,000) |
Antibody | β-Tubulin (mouse monoclonal) | Santa Cruz Biotechnology | RRID:AB_2288090 | WB (1:200) |
Antibody | GAPDH (chicken polyclonal) | Merck Millipore | RRID:AB_10615768 | WB (1:10,000) |
Antibody | YAP (mouse monoclonal) | Santa Cruz Biotechnology | RRID:AB_10612397 | PLA (1:25) |
Antibody | pAKT (mouse monoclonal) | Cell Signaling Technology | RRID:AB_331158 | PLA (1:200) |
Antibody | YAP (rabbit monoclonal) | Cell Signaling Technology | RRID:AB_2650491 | IHC (1:50) |
Antibody | AKT (rabbit polyclonal) | Cell Signaling Technology | RRID:AB_329827 | IHC (1:50) |
Commercial assay or kit | ROS assay kit | Abcam | ab287839 | – |
Sequence- based reagent | YAP (human) – forward | This paper | NM_006106 | CCTGCGTAGCCAGTTACCAA |
Sequence- based reagent | YAP (human) – reverse | This paper | NM_006106 | CCATCTCATCCACACTGTTC |
Sequence- based reagent | ANKRD1 (human) – for | This paper | NM_014391.3 | AGTAGAGGAACTGGTCACTGG |
Sequence- based reagent | ANKTD1 (human) – rev | This paper | NM_014391.3 | TGGGCTAGAAGT GTCTTCAGA T |
Sequence- based reagent | CYR61 (human) – for | This paper | NM_001554.5 | AGCCTCGCATCCTATACAACC |
Sequence- based reagent | CYR61 (human) – rev | This paper | NM_001554.5 | TTCTTTCACAAGGCGGCACTC |
Sequence- based reagent | GAPDH (human) – for | This paper | NM_002046.7 | CTGGTAAAGTGGATATTGTTGCCAT |
Sequence- based reagent | GAPDH (human) – rev | This paper | NM_002046.7 | TGGAATCATATTGGAACATGTAAACC |
Sequence- based reagent | RPL41 (human) – for | This paper | NM_001035267 | AAACCTCTGCGCCATGAGAG |
Sequence- based reagent | RPL41 (human) – rev | This paper | NM_001035267 | AGCGTCTGGCATTCCATGTT |
Sequence- based reagent | SRDF4 (human) – for | This paper | NM_005626 | TGCAGCTGGCAAGACCTAAA |
Sequence- based reagent | SRSF4 (human) – rev | This paper | NM_005626 | TTTTTGCGTCCCTTGTGAGC |
Sequence- based reagent | B2M (human) – for | This paper | NM_004048 | CACGTCATCCAGCAGAGAAT |
Sequence- based reagent | B2M (human) – rev | This paper | NM_004048 | TGCTGCTTACATGTCTCGAT |
Software, algorithm | ASAP software | https://github.com/computationalpathologygroup/ASAP | v1.6 | – |
Software, algorithm | ImageJ | Rueden et al., 2017 | v1.53f51 | – |
Software, algorithm | Weka segmentation | Arganda-Carreras et al., 2017 | v3.3.1 | – |
Software, algorithm | Ilastik | Berg et al., 2019 | v1.3.3 | – |
Software, algorithm | Spatial Model Editor | https://spatial-model-editor.github.io | v1.2.1 | – |
Software, algorithm | sme-contrib | https://spatial-model-editor.github.io | v0.014 | – |
In vitro experiments
Cell culture and the establishment of genetically modified cells
Request a detailed protocolThe hepatocyte-derived cell lines Hep3B (#ACC93, DSMZ, Braunschweig, Germany), HLF (JCRB, Japan), HepG2, and SNU182 (ATCC/LGC Standards, Wesel, Germany) were cultured in Minimum Essential Medium (MEM), Dulbeccos Modified Eagle Medium (DMEM), and Roswell Park Memorial Institute (RPMI) medium, respectively. The media were supplemented with 10% FCS and 1% penicillin–streptomycin (Sigma-Aldrich, Taufkirchen, Germany). Cells were maintained at 37°C in a 5% CO2 atmosphere. Hep3B cells were transduced with vectors carrying cDNAs for human histone H2B, YAP, and TAZ/WWTR1 genes, fused with mCerulean, mVenus, and mCherry genes, respectively. The plasmid pRRLN-EF1α-mVenus-YAP was kindly provided by Prof. Dr. Alexander Loewer (TU Darmstadt), pLentiPGK-mCerulean-H2B was purchased (Addgene, Watertown, USA; No. 90234). The pRRLN-EF1α-d2mCherry-TAZ was cloned. Correctness of vectors was verified by end-to-end sequencing.
Vectors were stably integrated into the Hep3B cells using lentiviral particles. For this, HEK293 cells were transfected with lentiviral vectors, packaging (psPAX2, Addgene, Watertown, USA; No. 12260), and envelope plasmids (pMD2.G, Addgene, Watertown, USA; No. 12259) using polyethylenimine and incubated for 40 hr. Subsequently, virus particles in the cell culture medium were collected and used for infections. The transduced cells were selected for stable vector integration using geneticin (0.6 mg/ml, pRRLN-EF1α-mVenus-YAP), hygromycin (0.6 mg/ml, pLentiPGK-mCerulean-H2B), and blasticidin (2 µg/ml, pRRLN-EF1α-d2mCherry-TAZ). Cells were frequently checked for mycoplasma contamination and authenticated by short tandem repeat analysis (DSMZ, Braunschweig, Germany).
The APAP stock solution (100 mM) was freshly prepared by dissolving 0.3 g APAP crystalline powder in 20 ml phosphate-buffered saline (PBS, GE Healthcare, Solingen, Germany) under heating (42°C) and stirring (Sigma-Aldrich, St. Louis, USA). The APAP solution was filtered (Millex-GS filters, 0.22 µm; Merck Millipore Ltd, Carrigtwohill, Ireland) and immediately administered to cells to a final concentration of 10 mM. TBHP (Luperax, Sigma-Aldrich) with 70 wt% in H2O was applied to cells to a final concentration of 300 µM for 2 hr prior protein extraction. Hydrogen peroxide (H2O2) 30% (Rotipuran, Carl Roth GmbH, Karlsruhe, Germany) was applied to a final concentration of 2 mM for 1 hr.
AKT1/2 phosphorylation was inhibited with InSolution Akt Inhibitor VIII (1–10 µM, Merck KGaA, Darmstadt, Germany) in FCS-free medium for 3 hr prior to treatment with HGF (10 ng/µl).
Gene silencing by RNAi
Request a detailed protocolRNAi was performed to inhibit YAP gene expression. For this, 2 × 105 cells were seeded per well in a 6-cm-well plate and incubated overnight. Oligofectamine (Invitrogen, Carlsbad, CA, USA) and Opti-MEM (Gibco Life Technologies Ltd, Paisley, UK) were used for transfection of the siRNA (final concentration: 5 nM) according to the manufacturer’s instructions. A random sequence oligonucleotide was used as a control treatment (control siRNA: ctrl siRNA). The following siRNA sequences were used: YAP siRNA #1 (5′–3′): CCA CCA AGC UAG AUA AAG A-dT-dT and YAP siRNA #2 (5′–3′): GG UCA GAG AUA CUU CUU AA-dT-dT, and ctrl siRNA (5′–3′): UGG UUU ACA UGU CGA CUA A. YAP siRNA #1 and #2 were pooled to achieve an optimal gene knockdown. Twenty-four hours after siRNA transfection, cells were treated with APAP (10 mM) for 24 hr.
Live cell imaging
Request a detailed protocolStably transduced Hep3B cells were seeded on 24-well glass-bottom black wall plates (MoBiTec GmbH, Goettingen, Germany) in a phenol red-free RPMI 1640 medium (Gibco/Life Technologies Corporation, Paisley, UK). Dynamic protein localization in living cells was measured using confocal laser scanning microscope (Nikon A1R on an inverted Nikon Ti2 microscope), which was supplemented with an on-stage incubator (TokaiHit) to maintain 37°C and 5% CO2 at all stages of the experiments. Selected images were acquired with additional Nikon AxR microscope system.
Three lasers were employed: 445 nm (mCerulean, PMT detector), 514 nm (mVenus, GaAsP detector), and 561 nm (mCherry, GaAsP detector) to obtain highly resolved live cell images. Two channels – 445 and 514 nm – were acquired simultaneously. The employed pinhole size was 17.9 µm (445 nm: 1.5 a.u., 514 nm: 1.3 a.u., 561 nm: 1.2 a.u.). Images were acquired using a galvano-scanner of size 512 × 512 px (FOV 0.64 × 0.64 mm, pixel resolution of 1.24 µm). The objectives Plan Apo λ 20× (NA 0.75, working distance 1 mm, FOV 0.64 × 0.64 mm) or Plan Apo λ 40× (NA 0.95, working distance 0.21 mm, FOV 0.435 × 0.435 mm) were used.
Western immunoblotting
Request a detailed protocolWestern immunoblotting was performed as previously described (Weiler et al., 2017). In brief, for total protein fraction preparation, cultured cells were harvested using 1× Cell Lysis Buffer (Cell Signaling Technology, Danvers, USA) supplemented with 1× Protease Inhibitor Mix G (Serva, Heidelberg, Germany) and 1× PhosSTOP (Roche Diagnostics Deutschland GmbH, Mannheim, Germany). For protein fractionation, the NE-PER Nuclear and Cytoplasmic Extraction kit was used (Thermo Scientific, Rockford, USA). The following antibodies were used for western immunoblotting experiments: YAP XP (1:1000, Cell Signaling Technology, Danvers, USA, #14074, RRID:AB_2650491), pYAP (1:400, Cell Signaling Technology, #4911, RRID:AB_2218913), pAKT (1:1’000, Cell Signaling Technology, #4060, RRID:AB_2315049), AKT (1:1’000, Cell Signaling Technology, #9272, RRID:AB_329827), β-actin (1:10,000, MP Biomedicals, Solon, USA, #08691001, RRID:AB_2335127), LATS1 (1:200, Santa Cruz Biotechnology, Dallas, USA, sc-398560), LATS2 (1:200, Santa Cruz Biotechnology, sc-515579), and pLATS1/2 (1:500, Cell Signaling Technology, #8654, RRID:AB_10971635) and Cytochrome P450 2E1 (1:500, Novus Biologicals, Centennial, USA, NVP1-85367, RRID:AB_11021447). PARP (1:10,000, Cell Signaling Technology, #9542, RRID:AB_2160739), β-tubulin (1:200, Santa Cruz Biotechnology, sc-5274, RRID:AB_2288090), and GAPDH (1:10,000, Merck Millipore, Darmstadt, Germany, ab2302, RRID:AB_10615768) were used as loading controls.
Western blot detection and quantification were performed using the Odyssey-CLx Infrared Imaging system with the ImageStudio software (LI-COR Biosciences, Bad Homburg, Germany). Phosphorylated protein band intensity was measured and normalized to total protein concentrations. Equal amount of protein was loaded for each line, as measured by the Bradford reagent (Millipore Sigma, Saint Louis, USA). Raw unedited image files and uncropped blots are available as source data files of this manuscript.
RNA isolation, reverse transcription, and quantitative real-time PCR (qPCR)
Request a detailed protocolRNA was isolated using Extractme kit (Blirt, Gdańsk, Poland) according to the manufacturer’s protocol. Reverse transcription was performed using the RevertAid kit (Thermo Scientific). Semiquantitative real-time PCR reactions were set up using the primaQuant 2× qPCR-SYBR-Green-Mastermix (Steinbrenner Laborsysteme, Wiesenbach, Germany) and analyzed with the QuantStudio 3 real-time PCR system (Applied Biosystems, Thermo Fisher Scientific, Singapore). The following cycling conditions were applied: 95°C for 15 min, followed by 40 cycles of 95°C for 15 s, and 60°C for 60 s. Product specificity was confirmed by melting curve analysis (95°C for 15 s, 60°C for 30 s, 60–95°C with 0.5°C/s).
The mRNA levels were normalized to glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 60S ribosomal protein L41 (RPL41), serine/arginine-rich splicing factor (SRSF4), and β2-microglobulin (B2M). The following primers for human cDNAs were used: YAP-for: 5′-CCT GCG TAG CCA GTT ACC AA-3′; YAP-rev: 5′-CCA TCT CAT CCA CAC TGT TC-3′; ANKRD1-for: 5′-AGT AGA GGA ACT GGT CAC TGG-3′; ANKRD1-rev: 5′-TGG GCT AGA AGT GTC TTC AGA T-3′; CYR61-for: 5′-AGC CTC GCA TCC TAT ACA ACC-3′; CYR61-rev: 5′-TTC TTT CAC AAG GCG GCA CTC-3′; GADPH-for: 5′-CTG GTA AAG TGG ATA TTG TTG CCA T-3′; GAPDH-rev: 5′-TGG AAT CAT ATT GGA ACA TGT AAA CC-3′; RPL41-for: 5′-AAA CCT CTG CGC CAT GAG AG-3′; RPL41-rev: 5′-AGC GTC TGG CAT TCC ATG TT-3′; SRSF4-for: 5′-TGC AGC TGG CAA GAC CTA AA-3′; SRSF4-rev: 5′-TTT TTG CGT CCC TTG TGA GC-3′; B2M-for: 5′-CAC GTC ATC CAG CAG AGA AT-3′; B2M-rev: 5′-TGC TGC TTA CAT GTC TCG AT-3′. For the analysis of gene expression in tissue samples, a panel of housekeeping genes was analyzed using the geNorm algorithm to find the most stable reference gene (Vandesompele et al., 2002).
In situ PLA
Request a detailed protocolThe DuoLink in situ PLA was performed according to the manufacturer’s instructions (Sigma-Aldrich). Briefly, cells were seeded on glass coverslips and prior to APAP administration cells were grown under FCS-free conditions for 1 day, then incubated with APAP (10 mM) or PBS. After treatment, cells were washed three times with 2 mM MgCl2 in PBS and fixed with 4% paraformaldehyde for 10 min at room temperature. Fixed cells were washed four times with PBS for 5 min, permeabilized with 0.2% Triton X-100 in PBS for 5 min at room temperature, and washed again twice with PBS for 5 min. Subsequently, cells were blocked with Blocking solution (Sigma-Aldrich) for 30 min at room temperature and incubated with the following primary antibodies diluted in Antibody Diluent (Sigma-Aldrich) overnight at 4°C: anti-YAP (1:25, Santa Cruz Biotechnology, sc-271134, RRID:AB_10612397), anti-AKT (1:200, Cell Signaling Technology, #9272, RRID:AB_329827), anti-pYAP (1:200, Cell Signaling Technology, #4911, RRID:AB_2218913), anti-pAKT (1:200, Cell Signaling Technology, #4051, AB_331158), and anti-pLATS1/2 (1:200, Cell Signaling Technology, #8654, RRID:AB_10971635). Subsequently, cells were washed twice with Wash Buffer A (Sigma-Aldrich) and incubated with prediluted rabbit PLUS and mouse MINUS probes (Sigma-Aldrich) in Antibody Diluent for 1 hr at 37°C. After incubation, cells were washed twice with Wash Buffer A and then incubated with ligation solution (Sigma-Aldrich) for 30 min at 37°C. After ligation, samples were again washed twice with Wash Buffer A for 2 min at room temperature and incubated with amplification solution, and detection reagent Orange (Sigma-Aldrich) for 100 min at 37 °C. Finally, cells were washed twice with Wash Buffer B (Sigma-Aldrich) for 10 min at room temperature, once with 0.01× Wash Buffer B for 1 min at room temperature and coverslips were mounted on the slide with DAPI Fluoromount-G mounting medium (SouthernBiotech, Birmingham, USA).
Fluorescence images were captured using an inverted Nikon Ti2 microscope with Nikon S Plan Fluor ELWD ×40 NA 0.60 objective in a widefield fluorescence mode using Lumencor Sola SE II lamp. Images were captured in DAPI and TRITC channels (460 and 580 nm) with Nikon DS-Qi2 monochrome camera (image size 2404 × 2404 px, pixel resolution of 0.18 µm/px).
ROS activity measurement
Request a detailed protocolMeasurements of the intracellular ROS levels were performed using the DCFDA/H2DCFDA cellular ROS assay kit (Abcam, Amsterdam, Netherlands) according to the manufacturer’s protocol. In brief, cells were seeded on white clear-bottom 96-well plate (Corning, Corning, USA) and incubated overnight. Cells were treated with 10 mM APAP, 2 mM H2O2, or 300 µM TBHP (H2O2 and TBHP served as positive controls for ROS induction) for 6 hr in FCS-free cell culture medium. Fluorescence was measured using a microplate reader (FluoStar Omega, BMG Labtech GmbH, Ortenberg, Germany). Buffer solution without cells served as a background control.
In vivo experiments and sample analyses
Housing and treatment of mice and induction of acute liver injury by acetaminophen
Request a detailed protocolMale C57BL/6N mice (8- to 10-week-old) were bought from Janvier Labs (Janvier Labs, Le Genest-Saint-Isle, France). Animals were housed under 12 hr light/dark cycles at controlled ambient temperature of 25°C with free access to water and were fed ad libitum with a standard diet (Ssniff, Soest, Germany) before starting the experiments. Induction of acute liver injury with APAP was done as previously described (Schneider et al., 2021). Briefly, the mice were fasted overnight, then challenged with a dose of 300 mg/kg APAP intraperitoneally. APAP was dissolved in warm PBS with an application volume of 30 ml/kg. Control group was treated with PBS only. The mice were fed ad libitum after APAP administration. All animals were included for further analyses. All experiments were approved by the local animal welfare committee (LANUV, North Rhine-Westphalia, Germany, application number: 84-02.04.2016.A279). Hydrodynamic gene delivery experiments in mice were performed as recently described (Luiken et al., 2020).
Liver tissue sample collection, processing, and staining
Request a detailed protocolTissues were collected time dependently after APAP injection from the left liver lobe. The tissues were fixed and embedded in paraffin as previously described (Ghallab et al., 2016). YAP and AKT immunostaining were performed using 4-µm-thick paraffin-embedded tissue sections. For immunohistochemistry, an anti-YAP antibody (1:50, Cell Signaling Technology, #14074, RRID:AB_2650491) and anti-AKT (1:50, Cell Signaling Technology, #9272, RRID:AB_329827) were used. Embedded tissue sections were pretreated with a heat-induced epitope retrieval method (pH 6, DAKO, Hamburg, Germany). As secondary antibody anti-rabbit Polymer-AP (Enzo Life Sciences, Farmingdale, USA, ENZ-ACC110-0150) was used. Detection was performed with Permanent AP (Zytomed Systems GmbH, Berlin, Germany). Following staining, the whole slides were digitally documented using a slide scanner (Aperio AT2, Leica Mikrosysteme Vertrieb GmbH, Wetzlar, Germany).
Expression profiling and bioinformatics
Request a detailed protocolRNA isolation and gene array analysis were performed as published before (Campos et al., 2020) and bioinformatic analysis was done as described (Holland et al., 2022). The gene expression data are available under ArrayExpress accession number GSE167032. The expression data were applied to three known signatures that are informative for ROS activity (Han et al., 2008), AKT activity (De Marco et al., 2017), and YAP/TAZ activity (Wang et al., 2018). Due to the high number of AKT signature genes, only genes whose response to APAP treatment was larger than fold change of 2 (compared to control animals) were considered. The expression data were z-score normalized and clustered with seborn clustermap python module (v0.11.0). The signature score was obtained by summarizing z-scored expression values at the given time point if the z-scored value was greater than 0.5. The summarized expression values were normalized to the number of genes in the signature.
Computational methods
Analysis of the live cell images
Request a detailed protocolThe confocal images of the living cells were analyzed in a high-throughput manner using ImageJ (v1.53f51) platform (Rueden et al., 2017). Images were first manually selected with respect to the quality criteria: images with insufficient sharpness or with artifacts were discarded from the analysis. The image processing pipeline was based on Weka segmentation (v3.3.1) (Arganda-Carreras et al., 2017) of foreground and background areas, and subsequent thresholding, object detection, and counting. After object detection and counting, the respective masks were overlaid on the initial images to acquire YAP and TAZ intensity values for nuclei and cytoplasm. Mean pixel intensity from nuclear areas of the cells was divided by the mean pixel intensity of cytoplasmic regions, thus obtaining an NCR.
Analysis of PLA images
Request a detailed protocolPLA slides were analyzed using ImageJ (v1.53f51). First, nuclei and dots were classified using the Weka segmentation algorithm (v3.3.1) (Arganda-Carreras et al., 2017), thresholded, and counted. The pseudo-cytoplasmic (ring-shaped) area was created using the ImageJ’s binary mask option dilate for 30 iterations on nuclear masks to obtain nuclear and cytoplasmic area for a comparative analysis. The thresholded nuclei or cytoplasmic masks were overlaid with detected dots to obtain the information on the subcellular localization of the protein interaction.
Analysis of IHC images
Request a detailed protocolStained tissue samples were digitalized using Aperio slide scanner with ×40 magnification and pixel resolution of 0.253 µm/px (Aperio AT2, Leica Mikrosysteme Vertrieb GmbH, Wetzlar, Germany). The selected time points (6 hr and 1, 2, 6, 16 days) and control treatment were quantified using a pipeline, which consisted of python scripts (modules PIL v5.3.0, matplotlib v2.2.3), ASAP software (v1.6, https://github.com/computationalpathologygroup/ASAP) and Ilastik software (v1.3.3) (Berg et al., 2019).
First, digital images of the tissue sections were divided in tiles (1 mm2, 400 × 400 pixels), which were binned (to reduce processing time and storage load) using ASAP and python modules. Some images were excluded due to staining and/or scanning artifacts. Tiles, which displayed tissue for at least 50% of their area, were kept for further processing. Machine-learning model, which was based on a random forest algorithm, was trained on a selected set of training tiles for YAP or AKT using Ilastik software. The algorithm was trained to detect positively stained nuclei, excluding necrotic areas and staining artifacts. Ilastik software was further used to export probability maps, which were possessed with ImageJ. In ImageJ, probability maps were thresholded and positive nuclei were counted and normalized to the area of the tile, which was occupied by the tissue.
PDE modeling
Request a detailed protocolSpatial modeling aimed at describing the variations in space of fluorescently labeled YAP and TAZ distribution within living Hep3B cells. The mathematical model is based on a system of PDEs. PDE modeling and parameter estimation were performed with the Spatial Model Editor (SME) software (v1.2.1) and sme-contrib (v0.0.14) python module (https://spatial-model-editor.github.io/). SME is graphical user interface-based model editing and simulation software compatible with systems biology markup language (SBML) standards. Models were simulated using simple Forward Time Centered Space (FTCS) solver.
In the demonstrated models, all reactions were defined by first-order kinetics (for PDEs see Appendix 1). For the parameter estimation, the model behavior was evaluated at steady state, that is the model was simulated until the concentration distribution did not change over time. Parameter estimation was performed using the particle swarm algorithm (20 particles, 200 iterations) to minimize a cost function consisting of the weighted sum of two terms: the squared per pixel differences between model and the target image, and the sum of squares of species concentration rates of change. In total 200 fitted parametrizations for each tested model were generated. The parameter space for the optimization algorithm was defined based on published data (Appendix 1). The parameter ranges of our canonical and alternative models were selected in a way to describe biologically meaningful value ranges and to avoid numeric instability during PDE simulations of the canonical model. The mathematical models were uploaded to the BioModels repository under the model identifier number MODEL2202080001. PDE model equations and parameters can be found in Supplementary Information (Appendix 1—Tables 1–6; Ege et al., 2018; Jack et al., 1990). Tables showing all fitted parameters for YAP and TAZ are provided (Supplementary file 1, Supplementary file 2).
Statistics
Statistical analysis was performed using GraphPad Prism 9.2.0. Statistical tests are indicated in figure legends. Error bars depict standard deviation. Significance levels are as follows: *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001.
Appendix 1
PDE model equations
In the computational model, the equations which are simulated in cytoplasm or nucleus follow two-dimensional reaction diffusion equation. Initial conditions for the PDE are arbitrary since only steady-state solutions are considered. The equations in TAZ model are equivalent to YAP model equations.
Here, described in general terms as follows:
where
· is the concentration of species at position and time
· is the reaction term for species
· is the diffusion constant for species
· is the Laplacian of species , which can be rewritten as
The PDEs for the dynamics of YAPn/pYAPn (nuclear fraction) and YAP/pYAP (cytosolic fraction) proteins in the alternative model in the nucleus or cytoplasm are presented as follows.
where reaction rates are defined as follows:
R1 as translation:
R2 as degradation:
R5 as phosphorylation:
R6 as dephosphorylation:
The transport reactions for YAP and pYAP are defined as flux densities across the nuclear membrane that depends on the respective concentrations; these are translated into Neumann-type interface conditions for the inner boundary of the cytoplasm and the outer boundary of the nucleus by the simulator software.
Flux density of YAP import
Flux density of YAP export
Boundary conditions on the outer membrane of the cytoplasm are ‘zero-flux’ Neumann type for all variables.
Model parameters
Canonical model parameters.
Parameter | Symbol | YAP model | TAZ model | Unit |
---|---|---|---|---|
Diffusion | 0.0391776 | 0.00812799 | µm2/s | |
Diffusion phosphorylated | 0.433785 | 0.501116 | µm2/s | |
Phosphorylation | 3.2462731384467 | 9.2660885020866 | s–1 | |
Dephosphorylation | 2.4380091829261 | 2.9721728333365 | s–1 | |
Translation | 0.013295743091699 | 0.02329254213393 | a.u./s | |
Degradation | 0.0081058910608829 | 0.0090941179957213 | s–1 | |
Nuclear import | 9.7229134833679e−16 | 8.1106772583928e−16 | µm/s | |
Nuclear export | 9.5119620861314e−17 | 6.1050891990604e−17 | µm/s |
Alternative model parameters.
Parameter | Symbol | YAP model | TAZ model | Unit |
---|---|---|---|---|
Diffusion | 9.84707 | 9.14393 | µm2/s | |
Diffusion phosphorylated | 1.04518 | 1.00732 | µm2/s | |
Phosphorylation | 0.028839815689838 | 0.10286291836829 | s–1 | |
Dephosphorylation | 0.09121073484801 | 0.10517790271802 | s–1 | |
Translation | 0.009403085138566 | 0.0030932046019853 | a.u./s | |
Degradation | 0.0094299166245222 | 0.0097940321843269 | s–1 | |
Nuclear import | 6.5845623444598e−15 | 9.5850453479782e−15 | µm/s | |
Nuclear export | 7.8819897236132e−15 | 4.5876638224922e−15 | µm/s |
Parameters of the YAP-to-TAZ transition model.
Parameter | Symbol | Value | Unit |
---|---|---|---|
Diffusion | 10 | µm2/s | |
Diffusion phosphorylated | 1 | µm2/s | |
Phosphorylation | – | s–1 | |
Dephosphorylation | – | s–1 | |
Translation | 0.01 | a.u./s | |
Degradation | 0.004 | s–1 | |
Nuclear import | 7.5e−15 | µm/s | |
Nuclear export | 1.6e−15 | µm/s |
Parameter range for particle swarm algorithm of the canonical YAP/TAZ model.
Parameter | Symbol | Range (min, max) | Unit |
---|---|---|---|
Diffusion | (0.001, 5) | µm2/s | |
Diffusion phosphorylated | (0.001, 1) | µm2/s | |
Phosphorylation | (0.1, 10) | s–1 | |
Dephosphorylation | (0.1, 10) | s–1 | |
Translation | (1e−3, 0.1) | a.u./s | |
Degradation | (1e−3, 1e−2) | s–1 | |
Nuclear import | (1e−16, 1e−15) | µm/s | |
Nuclear export | (1e−17, 1e−15) | µm/s |
Parameter range for particle swarm algorithm of the alternative YAP/TAZ model.
Parameter | Symbol | Range (min, max) | Unit |
---|---|---|---|
Diffusion | (1, 10) | µm2/s | |
Diffusion phosphorylated | (1, 1.5) | µm2/s | |
Phosphorylation | (0.01, 0.5) | s–1 | |
Dephosphorylation | (0.01, 0.5) | s–1 | |
Translation | (1e−3, 0.01) | a.u./s | |
Degradation | (1e−3, 0.01) | s–1 | |
Nuclear import | (1e−16, 1e−14) | µm/s | |
Nuclear export | (1e−16, 1e−14) | µm/s |
Parameter values from the literature.
Parameter | Value | Source |
---|---|---|
Hepatocyte cell volume | 10–11 l | Jack et al., 1990 |
Hepatocyte nucleus volume | 5 ×10–13 | Jack et al., 1990 |
Nuclear export | 10–100 s | Ege et al., 2018 |
Nuclear import | 50 s | Ege et al., 2018 |
Diffusion rate YAP/TAZ | 19 μm2/s | Ege et al., 2018 |
Data availability
The mathematical models generated and analyzed during the current study are available in the BioModels repository (https://www.ebi.ac.uk/biomodels/) under accession number MODEL2202080001. The PDE simulation software is available over Github (https://spatial-model-editor.github.io/). The gene expression data analyzed during the current study are available in the Gene Expression Omnibus (GEO) repository under accession number GSE167032.
-
EMBL-EBI BioModelsID MODEL2202080001. MODEL2202080001.
-
NCBI Gene Expression OmnibusID GSE167032. Expression data of the livers of male C57Bl6/N mice after i.p. injection of paracetamol (APAP).
References
-
Ilastik: interactive machine learning for (BIO) image analysisNature Methods 16:1226–1232.https://doi.org/10.1038/s41592-019-0582-9
-
Pharmacological targeting of kinases MST1 and MST2 augments tissue repair and regenerationScience Translational Medicine 8:352ra108.https://doi.org/10.1126/scitranslmed.aaf2304
-
A model for the Hippo pathway in the Drosophila wing discBiophysical Journal 115:737–747.https://doi.org/10.1016/j.bpj.2018.07.002
-
The in vivo gene expression signature of oxidative stressPhysiological Genomics 34:112–126.https://doi.org/10.1152/physiolgenomics.00239.2007
-
Phosphoinositide 3-kinase/Akt signaling and redox metabolism in cancerFrontiers in Oncology 8:160.https://doi.org/10.3389/fonc.2018.00160
-
Role of CYP2E1 in the hepatotoxicity of acetaminophenThe Journal of Biological Chemistry 271:12063–12067.https://doi.org/10.1074/jbc.271.20.12063
-
Yap and TAZ regulate skin wound healingThe Journal of Investigative Dermatology 134:518–525.https://doi.org/10.1038/jid.2013.339
-
Hippo pathway coactivators YAP and TAZ are required to coordinate mammalian liver regenerationExperimental & Molecular Medicine 50:e423.https://doi.org/10.1038/emm.2017.205
-
The Hippo pathway effector proteins YAP and TAZ have both distinct and overlapping functions in the cellThe Journal of Biological Chemistry 293:11230–11240.https://doi.org/10.1074/jbc.RA118.002715
-
Yap and TAZ are not identical twinsTrends in Biochemical Sciences 46:154–168.https://doi.org/10.1016/j.tibs.2020.08.012
-
Protein interaction switches coordinate Raf-1 and MST2/hippo signallingNature Cell Biology 16:673–684.https://doi.org/10.1038/ncb2986
-
Intestinal dysbiosis amplifies acetaminophen-induced acute liver injuryCellular and Molecular Gastroenterology and Hepatology 11:909–933.https://doi.org/10.1016/j.jcmgh.2020.11.002
-
A computational model of YAP/TAZ mechanosensingBiophysical Journal 110:2540–2550.https://doi.org/10.1016/j.bpj.2016.04.040
-
A systematic assessment of the epidemiologic literature regarding an association between acetaminophen exposure and cancerRegulatory Toxicology and Pharmacology 127:105043.https://doi.org/10.1016/j.yrtph.2021.105043
-
Genome-Wide association between YAP/TAZ/TEAD and AP-1 at enhancers drives oncogenic growthNature Cell Biology 17:1218–1227.https://doi.org/10.1038/ncb3216
-
The Hippo pathway in organ size control, tissue regeneration and stem cell self-renewalNature Cell Biology 13:877–883.https://doi.org/10.1038/ncb2303
Decision letter
-
Matthew A QuinnReviewing Editor; Wake Forest School of Medicine, United States
-
Mone ZaidiSenior Editor; Icahn School of Medicine at Mount Sinai, United States
-
Matthew A QuinnReviewer; Wake Forest School of Medicine, United States
-
Dirk FeyReviewer; University College Dublin, Ireland
Our editorial process produces two outputs: i) public reviews designed to be posted alongside the preprint for the benefit of readers; ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.
Decision letter after peer review:
Thank you for submitting your article "Spatial modeling reveals nuclear phosphorylation and subcellular shuttling of YAP upon drug-induced liver injury" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Matthew A Quinn as the Reviewing Editor and Reviewer #1, and the evaluation has been overseen by a Reviewing Editor and Mone Zaidi as the Senior Editor. The following individual involved in the review of your submission has agreed to reveal their identity: Dirk Fey (Reviewer #2).
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
1) The model features very different diffusion rates for the canonical and alternative model, but is not explained why. For the estimation, different parameter ranges were specified for the canonical and alternative model, also for other parameters to be estimated, e.g. phosphorylation and dephosphorylation rates. This also lacks an explanation. What is the rationale behind this? Please update the discussion and/or methods to clarify this issue.
2) 200 parametrizations were obtained. Is the simulation result similar? Maybe include a panel in figure S2 in the supplement about this. Also, include a table with all fitted parameters.
3) Although of a dynamic nature, the model was not analyzed for the time-dependency of its simulations nor compared to time-course data, and the manuscript would benefit from a brief discussion, in particular with respect to the time-dependent biphasic response.
4) In Figures 2B and 2D, please highlight the boundary between nucleus and cytoplasm in the pictures, as the main difference between the two modeling results is the direction of gradient around the nucleus-cytoplasm boundary. In addition to the modeling result confined by the parameter space, please also provide an intuitive explanation of why the canonical model cannot account for the discrepancies between experiment and simulation.
5) Figure 3 on APAP regulation of YAP phosphorylation: As APAP generally induces toxicity, it makes sense to quantify cell death in the culture and limit the measurement to live cells. If possible, it will be nice to monitor the same cell population over the course of APAP or control treatment, and see how the nuclear/cytoplasmic ratio of YAP and TAZ change as a function of a cell's local environment. Also, did the proposed YAP phosphorylation in early APAP treatment (<6h) lead to YAP shuttling inside the nuclei?
6) The western immunoblot bands are hard to be interpreted by eye except when showing a binary result. Please add quantification for the plots in Figure 3C and Figure 4A, C, D, E, and F.
Reviewer #1 (Recommendations for the authors):
This is an interesting and thorough investigation of the regulation of Yap localization via its nuclear phosphorylation. While this is a novel concept and well-supported by the data provided, there are several unanswered questions that need to be resolved either by additional experiments or by an enhanced discussion as detailed below:
1) In figure 2, the authors demonstrate kinases able to phosphorylate Yap (LATS1/2) display a nuclear localization under low cell seeding conditions. However, is this localization of LATS1/2 dependent on cell density? Does a LATS1/2 knockdown or overexpression affect Yap phosphorylation and shuttling? Does APAP affect LATS1/2 localization? Answering these questions experimentally will give more credence to the claim that Yap regulatory machinery is located in the nucleus and affords functional consequences to those regulatory proteins.
2) Figure 4 shows the quantification of duolink particles for various proteins under control and APAP groups. However, the APAP treated group is not displayed. A representative image of the APAP duolink is needed.
3) From a mechanistic perspective the authors utilize an Akt inhibitor and show dephosphorylation of YAP. They also show that Akt inhibitors block APAP-induced phosphorylation. However, the authors do not determine the effects of the downstream sub-cellular localization of Yap. Determining if Akt inhibitors affect Yap localization during APAP treatment is needed to make a mechanistic link between APAP/Akt/Yap translocation.
4) The authors show APAP induces Yap translocation in vivo during DILI. While this demonstrates in vivo relevance to their in vitro findings, it does not convey any functional relevance in the modulation of the regenerative response. Does Akt inhibition affect hepatic Yap translocation and downstream survival or regeneration? Performing additional in vivo experiments modulating the Akt/ROS/Yap pathway is essential in order to draw conclusions on whether this is a viable therapeutic target in vivo during DILI.
Reviewer #2 (Recommendations for the authors):
1) The model features very different diffusion rates for the canonical and alternative model, but is not explained why. For the estimation, different parameter ranges were specified for the canonical and alternative model, also for other parameters to be estimated, e.g. phosphorylation and dephosphorylation rates. This also lacks an explanation. What is the rationale behind this?
2) 200 parametrizations were obtained. Is the simulation result similar? Maybe include a panel in figure S2 in the supplement about this. Also, include a table with all fitted parameters.
3) Although of a dynamic nature, the model was not analyzed for the time-dependency of its simulations nor compared to time-course data, and the manuscript would benefit from a brief discussion, in particular with respect to the time-dependent biphasic response.
Reviewer #3 (Recommendations for the authors):
I believe the conclusions are mostly supported by the result. My comments are mostly on data presentation and interpretation:
1. In Figures 2B and 2D, please highlight the boundary between nucleus and cytoplasm in the pictures, as the main difference between the two modeling results is the direction of gradient around the nucleus-cytoplasm boundary. In addition to the modeling result confined by the parameter space, please also provide an intuitive explanation of why the canonical model cannot account for the discrepancies between experiment and simulation.
2. Figure 3 on APAP regulation of YAP phosphorylation: As APAP generally induces toxicity, it makes sense to quantify cell death in the culture and limit the measurement to live cells. If possible, it will be nice to monitor the same cell population over the course of APAP or control treatment, and see how the nuclear/cytoplasmic ratio of YAP and TAZ change as a function of a cell's local environment. Also, did the proposed YAP phosphorylation in early APAP treatment (<6h) lead to YAP shuttling inside the nuclei?
3. The western immunoblot bands are hard to be interpreted by eye except when showing a binary result. Please add quantification for the plots in Figure 3C and Figure 4A, C, D, E, and F.
Finally, it would be nice if the authors could discuss (1) the role of nuclear phosphorylation in cell proliferation, tissue damage, and regeneration – i.e. whether the YAP/TAZ localization is the cause or consequence of cell density change, and (2) the biological implication of context-dependent effect of APAP in early and late treatment.
https://doi.org/10.7554/eLife.78540.sa1Author response
Essential revisions:
1) The model features very different diffusion rates for the canonical and alternative model, but is not explained why. For the estimation, different parameter ranges were specified for the canonical and alternative model, also for other parameters to be estimated, e.g. phosphorylation and dephosphorylation rates. This also lacks an explanation. What is the rationale behind this? Please update the discussion and/or methods to clarify this issue.
This question is answered in detail below (see Reviewer 2, comment #1).
2) 200 parametrizations were obtained. Is the simulation result similar? Maybe include a panel in figure S2 in the supplement about this. Also, include a table with all fitted parameters.
This question is answered in detail below (see Reviewer 2, comment #2).
3) Although of a dynamic nature, the model was not analyzed for the time-dependency of its simulations nor compared to time-course data, and the manuscript would benefit from a brief discussion, in particular with respect to the time-dependent biphasic response.
This question is answered in detail below (see Reviewer 2, our reply to the general comment as well as comment #3).
4) In Figures 2B and 2D, please highlight the boundary between nucleus and cytoplasm in the pictures, as the main difference between the two modeling results is the direction of gradient around the nucleus-cytoplasm boundary. In addition to the modeling result confined by the parameter space, please also provide an intuitive explanation of why the canonical model cannot account for the discrepancies between experiment and simulation.
This question is answered in detail below (see Reviewer 3, comment #1).
5) Figure 3 on APAP regulation of YAP phosphorylation: As APAP generally induces toxicity, it makes sense to quantify cell death in the culture and limit the measurement to live cells. If possible, it will be nice to monitor the same cell population over the course of APAP or control treatment, and see how the nuclear/cytoplasmic ratio of YAP and TAZ change as a function of a cell's local environment. Also, did the proposed YAP phosphorylation in early APAP treatment (<6h) lead to YAP shuttling inside the nuclei?
This question is answered in detail below (see Reviewer 3, comment #2).
6) The western immunoblot bands are hard to be interpreted by eye except when showing a binary result. Please add quantification for the plots in Figure 3C and Figure 4A, C, D, E, and F.
This question is answered in detail below (see Reviewer 3, comment #3).
Reviewer #1 (Recommendations for the authors):
This is an interesting and thorough investigation of the regulation of Yap localization via its nuclear phosphorylation. While this is a novel concept and well-supported by the data provided, there are several unanswered questions that need to be resolved either by additional experiments or by an enhanced discussion as detailed below:
1) In figure 2, the authors demonstrate kinases able to phosphorylate Yap (LATS1/2) display a nuclear localization under low cell seeding conditions. However, is this localization of LATS1/2 dependent on cell density? Does a LATS1/2 knockdown or overexpression affect Yap phosphorylation and shuttling? Does APAP affect LATS1/2 localization? Answering these questions experimentally will give more credence to the claim that Yap regulatory machinery is located in the nucleus and affords functional consequences to those regulatory proteins.
We thank reviewer for the suggestions. The mechanistic connection between LATS1/2 kinases and the Hippo pathway effectors has been intensively investigated in different species and cell types. To further emphasize the relevance of the LATS1/2-YAP connection in nuclei several experiments were initiated.
(A) To answer the question whether LATS1/2 localization is dependent on cell density, we performed nuclear-cytoplasmic fractionation experiments for low and high cell density conditions (Author response image 1A). Indeed, cell density conditions affected LATS1 and LATS2 localization differently. The amount of nuclear LATS1 protein was increasing at high cell density (NCR 1.4 at low cell density to 4.1 at high cell density). In contrast, LATS2, which was predominantly localized in the cytoplasm, did not shuttle in a comparable manner. Phosphorylated LATS1/2 (pLATS1/2) increased in the nucleus at high cell density conditions (NCR from 0.7 to 1.1), indicating that the LATS-dependent regulatory processes take place at high cell density in cell nuclei. To sum up, these experiments demonstrate that nuclear LATS abundancy is indeed regulated by cell density.
(B) To answer the question whether LATS knockdown affects YAP phosphorylation and shuttling, we inhibited LATS1/2 expression using gene-specific siRNAs targeting LATS1 and LATS2 (Author response image 1B). As illustrated by Western blot quantification, LATS2 silencing reduced phosphorylation of YAP about 30%, whereas LATS1 knockdown did not significantly alter pYAP abundance. Thus, both LATS proteins contribute YAP in our experimental setup differentially.
(C) We further quantitatively explored the localization of YAP after siLATS treatment using live-cell imaging (Author response image 1C). Data illustrated significantly increased NCR for YAP in siLATS treated cells versus controls cells at low and high cell density conditions (see quantification).
Together with the data we showed in the first version of the manuscript (nuclear localization of LATS1/2 and physical interaction between active LATS1/2 and YAP in nuclei), we conclude that nuclear LATS is indeed an important player in regulating Hippo pathway.

Western blot of nuclear and cytoplasmic protein fractions derived from Hep3B cells under low and high cell density conditions (5x105 for low cell density and 3x106 cells for high cell density per 10 cm dish).
Quantification of the NCR for LATS1, LATS2, and phosphorylated LATS (pLATS1/2) under low and high cell density is depicted on the right. C: cytoplasmic fraction; N: nuclear fraction. (B): Inhibition of LATS proteins in Hep3B cells using siRNAs at different concentrations for 4 hr (20 and 40 nM; (siLATS1, 5’-CCA GAA GGA UAU AGA CAA AdTdT-3’; siLATS2, 5’-CAU GAA GAC CCU AAG GAA AdTdT-3’). Quantification of LATS knockdown efficiency and effect on YAP phosphorylation is shown on the right. Ctrl: untreated control, NTC: no template control. Cytoplasmic (C) and nuclear (N) protein fraction. (C): Live-cell imaging of Hep3B cells after LATS inhibition (siLATS1/2) and NTC (no template control) treatment (left). Scale bar 50 µm. Quantification of the NCR between NTC (n=60) and siLATS (n=55) with respect to cell count per image (representing cell density). Shaded area around the green and yellow curves represent 96% confidence interval of the fitted log(x) function. Statistics: two-tailed Mann-Whitney test with p-value<0.0001.
(D) Finally, we aimed to answer the question whether APAP affects LATS localization in our model system (Author response image 2). In cytoplasmic-nuclear fractionation experiments we observed that APAP treatment causes a moderate decrease in nuclear fraction of LATS1. Here, LATS2 localization was strongly influenced by APAP treatment, especially after 48 hours of treatment (NCR reduced about 60% in comparison to control). In addition, the localization of the phosphorylated LATS isoform (pLATS1/2) was also regulated by APAP. In general, APAP treatment led to reduced pLATS in nuclei and cytoplasm. Taken together, this data illustrates that APAP treatment reduced total amount of LATS1, LATS2 and pLATS1/2 proteins and also induced the cytoplasmic shift for those proteins. Thus, LATS proteins are less active after APAP administration and may contribute to the observed less efficient YAP phosphorylation.

Western blotting of cytoplasmic (C) and nuclear (N) fraction 2 and 48 hours after APAP (10 mM) and control (ctrl) treatment in Hep3B cells.
The localization of LATS1, LATS2 and pLATS1/2 proteins was detected. PARP and Tubulin served as loading controls for nucleus and cytoplasm, respectively. The quantification of the NCR of LATS1, LATS2 and pLATS1/2 shown on the right. Western blots were obtained from one experiment.
In summary, the performed experiments confirming that nuclear LATS1/2 protein are of importance for YAP shuttling. In addition, APAP (at least partly) may contribute to the subcellular localization of YAP via the regulation/exclusion of LATS1/2 from hepatocellular cell nuclei.
Because this information is of relevance for the readership, a clarifying statement was added in the revised version of the manuscript: " Indeed, we and others showed that nuclear LATS1/2 can control YAP phosphorylation and localization (data not shown) (Ege et al., 2018; Li et al., 2014). The relevance of this nuclear interaction is substantiated by a proximity ligation assay (PLA), which illustrated that pLATS1/2 physically interacted with YAP not only in the cytoplasm but also in the nucleus (Figure 2G)." Page 7, lines 251-256
2) Figure 4 shows the quantification of duolink particles for various proteins under control and APAP groups. However, the APAP treated group is not displayed. A representative image of the APAP duolink is needed.
We thank the reviewer for the suggestion. According to the comment, we now added representative images for the experiment in Figure 4I to the manuscript.
3) From a mechanistic perspective the authors utilize an Akt inhibitor and show dephosphorylation of YAP. They also show that Akt inhibitors block APAP-induced phosphorylation. However, the authors do not determine the effects of the downstream sub-cellular localization of Yap. Determining if Akt inhibitors affect Yap localization during APAP treatment is needed to make a mechanistic link between APAP/Akt/Yap translocation.
We thank the reviewer for the suggestion. Indeed, the mechanistic link between AKT and YAP was not intensively investigated in the first manuscript draft.
Before we present/discuss our new results, we want to emphasize that long-term experiments with AKTi (24 or 48 hours) were difficult to perform due to (unspecific?) cell toxic effects of the inhibitor. Instead, all experiments were performed at early time points (up to 3 hours), which is of special importance when interpreting those experiments where APAP is equally given (early APAP is leading to YAP phosphorylation, see figure 3C of main manuscript).
Since the reviewer asked for the effects of AKTi on YAP localization, we also performed short- and long-term life cell imaging experiments using fluorescently-tagged YAP after APAP and AKTi treatment. However, it turned out that the results couldn’t be analyzed for technical reasons. In detail, combined APAP/AKTi treatment disintegrated the signal of the nuclear marker (H2B-tagged cerulean), which prevented the quantitative analysis of nuclear and cytoplasmic YAP (Author response image 3C). Indeed, the AKT-dependent regulation of DNA fragmentation and histone deacetylase proteins has been described before (Ahn et al., 2006; Ozaki et al., 2010). Because the impaired nuclear signal quality led to unreliable interpretation of the live-cell imaging results, these results cannot be used to substantiate our findings.
Based on the data derived from experiments shown in Author response image 3A/B we conclude that AKTi leads to a (moderate) nuclear localization of YAP early (!) after APAP treatment. Thus, blocking AKT at early time points after APAP treatment prevented YAP phosphorylation (which is observed at early time-points) and therefore reduced nuclear exclusion.
We adjusted a statement in the revised manuscript: "Indeed, blocking AKT activity by AKTi reduced YAP phosphorylation, moderately elevated nuclear YAP positivity, and increased the physical interaction between AKT and YAP in cell nuclei (Figure 4F and data not shown)." Page 9, line 351 – page 10, line 354

(A): Left: nuclear and cytoplasmic protein fractionation after APAP (10 mM, 1h in FCS-free medium) and combined APAP and AKTi (AKT inhibitor VIII, 10 µM 3h in FCS-free medium) treatment of Hep3B cells.
PARP and Tubulin serve as loading controls for nuclear and cytoplasmic fraction, respectively. Right: quantification of the NCR of YAP after APAP and combined APAP and AKTi treatment. (B): DuoLink PLA between YAP and AKT proteins under APAP (10 mM, 2h) and combined APAP and AKTi (10 µM, 3h) treatment. Top row: nuclei (blue) and YAP-AKT interaction (red dots); bottom row: segmentation of the nuclei and PLA dots. Right: quantification of the PLA dots in nuclei for APAP (27 images) and combined APAP/AKTi treatment (31 images). Statistics: two tailed Mann-Whitney test, p-value<0.0001. Scale bar 50 µm. (C): Live-cell imaging of Hep3B cells under APAP (10 mM 2h, top row) and combined APAP and AKTi (10µM 3h, bottom row) treatment.
4) The authors show APAP induces Yap translocation in vivo during DILI. While this demonstrates in vivo relevance to their in vitro findings, it does not convey any functional relevance in the modulation of the regenerative response. Does Akt inhibition affect hepatic Yap translocation and downstream survival or regeneration? Performing additional in vivo experiments modulating the Akt/ROS/Yap pathway is essential in order to draw conclusions on whether this is a viable therapeutic target in vivo during DILI.
We thank the reviewer for the suggestion to further investigate the therapeutic relevance of AKT and YAP inhibition during APAP-induced DILI. We discussed this aspect with our collaboration partners (A. Ghallab and J. Hengstler, Department of Toxicology, Dortmund). However, we decided not to perform e.g., AKT inhibition experiments in vivo for ethical and experimental reasons:
First, the inhibitor used in our study is very potent and efficiently abrogates AKT1/2 activity (see main Figure 4E). Although the genetic inactivation of AKT1/2 in specific cell types such as cardiomyocytes can be investigated in vivo, the general blockade of AKT1/2 can cause rapid lethality as illustrated in newborn and adult mice with AKT1- and AKT2-deficiency (Liu et al., 2019; Peng et al., 2003; Wang et al., 2016). Thus, it is likely that the complete inhibition of both AKT isoforms by using an efficient chemical inhibitor may cause unpredictable side effects or is not compatible with life.
Second, as seen for APAP administration, the genetic inhibition of AKT isoforms in liver tissue cause severe hepatotoxicity and liver injury (Black, 1984; Wang et al., 2016). This strongly suggests that the treatment with APAP and concomitant silencing of AKT isoforms would further intensify hepatotoxicity leading to the unpredictable behavior of YAP expression and localization.
Due to these obstacles, we believe that these experiments are ethically problematic and very likely do not provide generate meaningful results.
However, we decided to perform additional mouse experiments to substantiate our experimental conclusions that AKT supports YAP activity in hepatocytes. For this, we did hydrodynamic gene delivery experiments in vivo with vectors coding for constitutively active myristoylated AKT (myrAKT) (Luiken et al., 2020). According to our hypothesis, permanent myrAKT activity should support nuclear YAP enrichment followed by the induction of typical transcriptional target genes (e.g., Ctgf, Cyr61). Indeed, immunohistochemical (IHC) tissue staining of extracted mice liver illustrated not only illustrated elevated AKT expression followed by hepatocellular steatosis (Wang et al., 2015) but also nuclear YAP positivity. Moreover, both investigated target genes were induced at the transcript level in mice after myrAKT injection (Author response image 4A/B).

(A): Immuno-histochemical staining of AKT (left) and YAP (right) of tissues derived from FVB mice after hydrodynamic injection of expression vectors coding myrAKT.
Experiments and ethical aspects have been published previously (Luiken et al., 2020). Scale bar 50 µm. (B): Quantitative polymerase chain reaction (qPCR) measurements of murine Hippo target genes Ctgf and Cyr61 (biological replicates: ncontrol = 3, nmyrAKT = 4). Statistics: unpaired two tailed t test; p-valueCtgf = 0.049, p-valueCyr61 = 0.026.
The results confirm that AKT support nuclear YAP enrichment. To keep the manuscript concise, we suggest to include the following clarifying statement in the manuscript: "The mechanistic connection between AKT activity and YAP induction in murine hepatocytes was confirmed in independent experiments. Here hydrodynamic gene delivery of myristoylated AKT led to nuclear enrichment of YAP expression in hepatocytes and expression of typical target genes (data not shown)." in paragraph ' Sequential activation of ROS, AKT, and Hippo/YAP in mouse livers after APAP intoxication'. Page 11, lines 402-407
Murine primers used for qPCR.
Expression of Cyr61 and Ctgf were normalized to a set of house-keeping genes: Actin, Gapdh, Tbp, Tubulin, Hprt, Ppia (software: GeNorm).
Gene | Forward | Reverse |
---|---|---|
Actin | GCTTCTTTGCAGCTCCTTCGT | ACCAGCGCAGCGATATCG |
Gapdh | TGTCCGTCGTGGATCTGAC | CCTGCTTCACCACCTTCTTG |
Tbp | TTGTCTGCCATGTTCTCCTG | CAGGGTGATTTCAGTGCAGA |
Tubulin | TCACTGTGCCTGAACTTACC | GGAACATAGCCGTAAACTGC |
Ctgf | GGAGAACTGTGTACGGAGCG | CCAGGCAAGTGCATTGGTA |
Hprt | TCCTCCTCAGACCGCTTTT | CCTGGTTCATCATCGCTAATC |
Ppia | GCATACAGGTCCTGGCATCT | AGCTGTCCACAGTCGGAAAT |
Cyr61 | GATCTGTGAAGTGCGTCCTTGTGG | GACACTGGAGCATCCTGCATAAG |
Reviewer #2 (Recommendations for the authors):
1) The model features very different diffusion rates for the canonical and alternative model, but is not explained why. For the estimation, different parameter ranges were specified for the canonical and alternative model, also for other parameters to be estimated, e.g. phosphorylation and dephosphorylation rates. This also lacks an explanation. What is the rationale behind this?
This comment/question involves two aspects: 1) Why are parameters that should be the same physical constant in both models different (e.g. diffusion coefficients)? 2) Why did we choose different value ranges for the fitting procedure in the different models?
In this work it was not our intention to exactly determine the parameter values. Indeed, all parameters are known to be structurally unidentifyable since only stationary data with arbitrary concentration units is used for the fitting. Choosing a different time scale for the model would change all parameter values while still resulting in the same steady state concentration pattern. This means the actual parameter values are de facto arbitrary. Instead, we focus on the qualitative result of our modeling: one version of the model cannot fit the data (given a wide range of allowed parameter values) while the other version can. Identical diffusion constants in both model versions would not strengthen our conclusion. Specifically, we have shown that the simple model without the same diffusion constants cannot reproduce the observations. Thus, adding an additional restriction on the diffusion constants will definitely not lead to a better fitting model.
The reason for different parameter ranges in both models is the numeric stability of the PDE solver for the canonical model. Extreme numeric values for the diffusion constants lead to very stiff models, which complicate PDE simulations. Therefore, parameter values for the canonical model were chosen so that the complete range could be reliably simulated for parameter fitting. As explained above, the actual parameter values are not identifyable and not the scope of this study. It is important that the range is big enough that a good fit for the simple model was not arbitrarily excluded. From our experience we are convinced that this condition is satisfied. We agree that identical ranges for both models would be 'nicer' and future versions of our simulator software may offer this possibility. However, for the reasons mentioned above we are convinced that our results are sound and can explain the biological observations sufficiently.
We want to emphasize that even in the absence of non-identifyability the parameter values resulting from a 'failed' fitting procedure are not physically meaningful. We have shown with extensive parameter estimation runs that we were not able to generate at least one single satisfying fit for the canonical model. Any specific parameter values from these unsuccessful fitting runs, even if not affected by structural non-identifyability, are inherently meaningless.
In order to clarify the discrepancies between parameters of canonical and alternative models, we added the following paragraph to our methods sections:
“The parameter ranges of our canonical and alternative models were selected in a way to describe biologically meaningful value ranges and to avoid numeric instability during PDE simulations of the canonical model.” Page 25, lines 842-844
2) 200 parametrizations were obtained. Is the simulation result similar? Maybe include a panel in figure S2 in the supplement about this. Also, include a table with all fitted parameters.
According to the reviewer's comment, we now provide a csv file with all obtained parameter values for each alternative model and the objective functions of the models. See Supplementary File 1 (for YAP) and 2 (for TAZ).
In addition, we provide additional information on distribution of the objective function across models (Figure 2—figure supplement 1A) and parameter values of the 30 best models (Figure 2—figure supplement 1B). Here, we overlaid the values of the 30 best models (based on the objective function minimum value after 200 iterations) with the parameter values form the rest of the models. In general, this illustrates the structural non-identifyability of all parameters in the models (as discussed in the manuscript), which is not a problem since we focus on the ability of the various model variants to fit the data rather than the values of the specific parameters. Interestingly, some specific parameters, such as, phosphorylation reaction in the nucleus (both models YAP and TAZ) and degradation parameter for TAZ models, seem to cluster in the 'good fit', possibly indicating the importance of these parameter values for the fit. Lastly, we show exemplarily how a “bad fit” looks in comparison to a 'good fit' (Author response image 5)

A visual representation of one of the best models for YAP and TAZ, and the worst models.
Parameter value distribution of YAP.
3) Although of a dynamic nature, the model was not analyzed for the time-dependency of its simulations nor compared to time-course data, and the manuscript would benefit from a brief discussion, in particular with respect to the time-dependent biphasic response.
We thank the reviewer for pointing out the time-dependency as an important player in Hippo signaling dynamics. This suggestion could indeed initiate an interesting and challenging follow up study. Model fitting and parameter estimation for spatial models is challenging even for not time-dependent data (indeed not many examples exist in the published literature) and using time-dependent data for fitting is not yet implemented in our software platform.
The reviewer is correct when he/she emphasizes that simulation of the time dependency could answer complex questions of dynamic nature. However, the initial question of our study (which pathway topology explains the data) was tailored to steady state models. A more complex model with further assumptions would indeed expand the view on the Hippo pathway, however, this would require further intensive work which would not significantly affect the outcome of our study and therefore could not be in the scope of the manuscript.
To clarify our aims regarding steady state PDE spatial modelling, we added the following passage to the discussion paragraph:
“The goal of our approach was to determine whether nuclear phosphorylation of YAP plays a role in Hippo signaling, specifically in governing YAP localization. A model not including nuclear phosphorylation of YAP could not sufficiently explain the experimentally observed localization pattern and was excluded from further analyses. However, the alternative model proposed in this work, can explain the experimental data and therefore establishes a possible mechanism for YAP localization. However, more complex mechanisms regulating subcellular localization and dynamics of YAP and TAZ cannot be excluded. Therefore, our model partly explains how the Hippo pathway is organized, however, it does not aim to comprehensively describe it." Page 12, lines 430-440
Reviewer #3 (Recommendations for the authors):
I believe the conclusions are mostly supported by the result. My comments are mostly on data presentation and interpretation:
1. In Figures 2B and 2D, please highlight the boundary between nucleus and cytoplasm in the pictures, as the main difference between the two modeling results is the direction of gradient around the nucleus-cytoplasm boundary. In addition to the modeling result confined by the parameter space, please also provide an intuitive explanation of why the canonical model cannot account for the discrepancies between experiment and simulation.
We thank the reviewer for the suggesting and now added a thin contour line around the nuclear border to highlight boundaries between cytoplasm and nuclei in main Figure 2B and 2D.
In addition, we agree with the reviewer that the inappropriateness of the canonical model (and why the extended model is more appropriate) can be explained better. In detail, the measured distribution of YAP and TAZ would imply diffusion out of the nucleus due to a gradient formed from nucleus to cytosol. Since YAP/TAZ are not synthesized in the nucleus, a second chemical species with different diffusion coefficients has to be exported out of the nucleus (the phosphorylated version of the protein). This prerequisite cannot be fulfilled by the canonical model.
This is obvious considering the residuals of the canonical model that clearly show discrepancy between simulation and the experimental results. Especially in areas close to the nuclear envelope in figure 2B/D, the canonical model is underrepresenting YAP/TAZ concentrations as depicted by deep blue pixels around the nucleus. Furthermore, the canonical model does not correctly represent the nuclear distribution of YAP and TAZ. Here, nuclear YAP and TAZ are underrepresented in the nuclear center (blue pixels), which is an inherent property of the canonical model topology.
To clarify the process behind the model selection, we added/adjusted the following passages in the results chapter of our manuscript:
“For example, the model cannot sufficiently explain the distribution of YAP and TAZ in cell nuclei or at the nuclear membrane. In detail, the cell area around the nuclear envelope on the cytoplasmic side strongly underrepresented the concentration of YAP/TAZ, as indicated by the blue pixels in the residual image. Moreover, the simulated YAP/TAZ distribution pattern of the canonical model showed a decreased protein concentration in the center of the nucleus, with increasing gradient towards the nuclear envelope. This indicates that the observation is dominated by a protein that is disseminated from the nucleus and undergoes diffusion and degradation. The canonical Hippo pathway cannot explain this effect, illustrated by residuals (experimental data minus model simulation) (Figure 2B)." Page 6, lines 203-214
2. Figure 3 on APAP regulation of YAP phosphorylation: As APAP generally induces toxicity, it makes sense to quantify cell death in the culture and limit the measurement to live cells. If possible, it will be nice to monitor the same cell population over the course of APAP or control treatment, and see how the nuclear/cytoplasmic ratio of YAP and TAZ change as a function of a cell's local environment. Also, did the proposed YAP phosphorylation in early APAP treatment (<6h) lead to YAP shuttling inside the nuclei?
Topic: Toxicity
The reviewer is right when he/she underlines the cell-toxic properties of APAP, which may differ in cell culture models. To investigate this aspect in more detail, we performed three different experimentals: CellTox assay, Western blotting of cleaved PARP fragment, and live-cell imaging.
We measured time-dependent cell-toxicity after APAP (10 mM) treatment with a CellTox Green Cytotoxicity Assay (Promega, # G8742, Author response image 6A). Here, APAP treatment for 48 hours moderately increased cell toxicity about 10% on average, whereas at 24 hours post treatment slightly decreased toxicity was observed. Although, mild toxicity is detectable in our experimental setup, we can conclude that the majority of cells is still viable at all time points of our experiments.
We quantified cleaved cytoplasmic PARP 89 kDa fragments whose presence indicate cell death (Author response image 6B). This experiment confirm that the amount of cleaved PARP slightly increased after APAP treatment for 48 hours.
More important, for all live-cell imaging experiments, which are the basis for our study and mathematical modeling, only viable (attached) cell populations with clear nuclear/cytoplasmic structure were used (see main Figures 3 and Figure 3—figure supplement 1A/B). Thus, quantification of the nuclear to cytoplasmic ratio was exclusively performed on viable cells, since rounded and floating (dead) cells were excluded from the classification algorithm.
These results illustrate a weak and neglectable role of APAP on cell toxicity and cell death in our experimental setup. Because the chosen live cell approach excludes damaged and/or apoptotic cells, we conclude that our experimental setup allows an unbiased view on APAP-induced effects on cell signaling.
An explanatory statement summarizing the results is now included in the revised manuscript: "Severe effects of APAP on cell toxicity and apoptosis in the chosen experimental setup were excluded by measuring cell viability and PARP cleavage (data not shown)." in paragraph 'APAP regulates YAP protein localization and activity'. Page 8, lines 281-283
Topic: early response
As to early APAP treatment (<6 hours), we indeed did not explicitly mention YAP localization in our previous manuscript, since we believe that 'later' time points better describe the biologically relevant APAP responses of a cell.
To answer the reviewer's question, we approached this aspect with two different experimental techniques: live-cell imaging and Western immunoblotting (Author response image 6B and 6C, respectively).
Live-cell imaging illustrates no significant nuclear/cytoplasmic shift of YAP after APAP treatment for 2 hours (Author response image 6C). In this experimental setup, we also investigated different cell density conditions to exclude this aspect as important parameter.
(B)These results indicate no or very weak cytoplasmic translocation of YAP. A possible explanation for this mild effect after short-term treatment of APAP might be the existence of molecular mechanisms (e.g., phosphorylation of alternative serine residues of YAP), which are not part of the current model. Due to the weak character of the results and the fact that we cannot offer a reasonable explanation for this observation, we would appreciate not to include these results in the revised version of the manuscript.
However, this aspect was already mentioned in the first submitted version of the manuscript. We wrote: " No obvious response was detectable for YAP and TAZ at earlier time points post APAP treatment (data not shown)." Page 8, lines 298-281

(A): CellTox assay quantifying non-viable cells in Hep3B cells after APAP (10 mM) treatment for 24 and 48 hours (n=10, technical replicates, independent biological replicates are not shown).
Statistics: two-tailed unpaired t test, p-value<0.0001. (B): Nuclear/cytoplasmic protein fractioning of Hep3B cells after APAP (10mM) treatment for 2 and 48 hours. Total (116 kDa) and cleaved PARP (89 kDa) were detected. Tubulin served as fractionation control. Quantification of cleaved PARP fragments is shown below. C: cytoplasmic; N: nuclear fraction. (C): Live-cell imaging of Hep3B cells expressing Venus-tagged YAP 2 hours after APAP (10 mM) or control treatment (PBS). Exemplary fluorescent images (left) and data quantification of NCR as a function of cell count per image (right) are shown. One (out of three) representative experiment is shown (nctrl = 68; nAPAP = 68). Scale bar 50 µm. (D): Nuclear/cytoplasmic protein fractionation after APAP treatment (10 mM, 2 hours). Detection of YAP and pYAP is shown. PARP and tubulin serve as loading controls. Quantification of NCR between YAP and pYAP is shown (right). C: cytoplasmic; N: nuclear fraction. Western immunoblots in (B) and (D) were obtained from one experiment.
3. The western immunoblot bands are hard to be interpreted by eye except when showing a binary result. Please add quantification for the plots in Figure 3C and Figure 4A, C, D, E, and F.
According to the reviewer's suggestion, we now included quantification of all blots of the revised version of the manuscript (Figure 3—figure supplement 2).
Finally, it would be nice if the authors could discuss (1) the role of nuclear phosphorylation in cell proliferation, tissue damage, and regeneration – i.e. whether the YAP/TAZ localization is the cause or consequence of cell density change, and (2) the biological implication of context-dependent effect of APAP in early and late treatment.
According to the reviewer's suggestion, two additional paragraphs were included in the Discussion chapter of the revised manuscript. The first paragraph focuses on the aspect how nuclear YAP phosphorylation can extent our current thinking about Hippo pathway biology. Whereas the second paragraph discusses our findings in the context of biological/medical implications after long-term APAP ingestion. For both paragraphs additional references were included in the manuscript.
“The extension of the canonical Hippo pathway model with processes that control active nuclear YAP phosphorylation illustrates that the regulatory YAP/TAZ phosphorylation must be considered as continuum process (Shreberk-Shaked and Oren, 2019). For instance, the existence of nuclear phosphorylation possibilities by LATS-dependent and independent (e.g., by AKT) mechanisms increases the complexity but also flexibility and redundancy of the signaling pathway (Ege et al. 2018, Gao et al. 2017, Low et al. 2017). These cellular aspects are part of a fine-tuned regulatory network that allows a rapid and adjustable proliferative response of YAP and TAZ under diverse physiological and pathological cellular conditions (e.g., upon induction of tissue regeneration).” Page 13, lines 483-493
“In our manuscript we show that APAP administration induces a chain of molecular events through the APAP/ROS/AKT axis, which is leading to YAP inactivation (early) followed by YAP activation (late). The late induction of YAP activity might be considered as cell protective cellular response upon long-term tissue damage; however, YAP also acts as a potent oncogene in different cell types, including hepatocytes, and its overexpression is associated with tumor initiation (Dong et al. 2007; Weiler et al. 2017). Thus, it is tempting to speculate that APAP overdose may not only cause DILI but also increases the risk for liver tumor development. Although controversially discussed in the literature, a recent evaluation of 139 published epidemiologic studies strongly argues against a relationship between APAP exposure and cancer (Weinstein et al. 2021). Actually, APAP has been discussed as therapeutic option for patients with YAP-induced cancer (Poudel et al. 2021). Although, our and other studies clearly demonstrate a mechanistic link between APAP update and Hippo/YAP pathway activity, these partly controversial findings illustrate the necessity to decipher this connection under distinct disease conditions in the future." Page 15, lines 538-554
References
Ahn JY, Liu X, Liu Z, Pereira L, Cheng D, Peng J, Wade PA, Hamburger AWYe K. 2006. Nuclear Akt associates with PKC-phosphorylated Ebp1, preventing DNA fragmentation by inhibition of caspase-activated DNase. EMBO J 25:2083-2095. doi:10.1038/sj.emboj.7601111.
Black M. 1984. Acetaminophen hepatotoxicity. Annu Rev Med 35:577-593. doi:10.1146/annurev.me.35.020184.003045.
Ege N, Dowbaj AM, Jiang M, Howell M, Hooper S, Foster C, Jenkins RPSahai E. 2018. Quantitative Analysis Reveals that Actin and Src-Family Kinases Regulate Nuclear YAP1 and Its Export. Cell Syst 6:692-708 e613. doi:10.1016/j.cels.2018.05.006.
Li W, Cooper J, Zhou L, Yang C, Erdjument-Bromage H, Zagzag D, Snuderl M, Ladanyi M, Hanemann CO, Zhou P, Karajannis MAGiancotti FG. 2014. Merlin/NF2 loss-driven tumorigenesis linked to CRL4(DCAF1)-mediated inhibition of the hippo pathway kinases Lats1 and 2 in the nucleus. Cancer Cell 26:48-60. doi:10.1016/j.ccr.2014.05.001.
Liu W, Jing ZT, Xue CR, Wu SX, Chen WN, Lin XJLin X. 2019. PI3K/AKT inhibitors aggravate death receptor-mediated hepatocyte apoptosis and liver injury. Toxicol Appl Pharmacol 381:114729. doi:10.1016/j.taap.2019.114729.
Luiken S, Fraas A, Bieg M, Sugiyanto R, Goeppert B, Singer S, Ploeger C, Warsow G, Marquardt JU, Sticht C, De La Torre C, Pusch S, Mehrabi A, Gretz N, Schlesner M, Eils R, Schirmacher P, Longerich TRoessler S. 2020. NOTCH target gene HES5 mediates oncogenic and tumor suppressive functions in hepatocarcinogenesis. Oncogene 39:3128-3144. doi:10.1038/s41388-020-1198-3.
Ozaki K, Kosugi M, Baba N, Fujio K, Sakamoto T, Kimura S, Tanimura SKohno M. 2010. Blockade of the ERK or PI3K-Akt signaling pathway enhances the cytotoxicity of histone deacetylase inhibitors in tumor cells resistant to gefitinib or imatinib. Biochem Biophys Res Commun 391:1610-1615. doi:10.1016/j.bbrc.2009.12.086.
Peng XD, Xu PZ, Chen ML, Hahn-Windgassen A, Skeen J, Jacobs J, Sundararajan D, Chen WS, Crawford SE, Coleman KGHay N. 2003. Dwarfism, impaired skin development, skeletal muscle atrophy, delayed bone development, and impeded adipogenesis in mice lacking Akt1 and Akt2. Genes Dev 17:1352-1365. doi:10.1101/gad.1089403.
Shreberk-Shaked MOren M. 2019. New insights into YAP/TAZ nucleo-cytoplasmic shuttling: new cancer therapeutic opportunities? Mol Oncol 13:1335-1341. doi:10.1002/1878-0261.12498.
Wang Q, Yu WN, Chen X, Peng XD, Jeon SM, Birnbaum MJ, Guzman GHay N. 2016. Spontaneous Hepatocellular Carcinoma after the Combined Deletion of Akt Isoforms. Cancer Cell 29:523-535. doi:10.1016/j.ccell.2016.02.008.
Wang YH, Viscarra J, Kim SJSul HS. 2015. Transcriptional regulation of hepatic lipogenesis. Nat Rev Mol Cell Bio 16:678-689. doi:10.1038/nrm4074
https://doi.org/10.7554/eLife.78540.sa2Article and author information
Author details
Funding
German Federal Ministry of Education and Research (031L0074H)
- Kai Breuhahn
Deutsche Forschungsgemeinschaft (505755359)
- Kai Breuhahn
Deutsche Forschungsgemeinschaft (314905040)
- Kai Breuhahn
German Federal Ministry of Education and Research (031L0158)
- Liam Keegan
- Ursula Kummer
- Sven Sahle
Heidelberg University (Research training group "Mathematical Modeling for the Quantitative Biosciences (MMQB)")
- Lilija Wehling
Deutsche Forschungsgemeinschaft (GH 276)
- Ahmed Ghallab
Ministerio de Ciencia e Innovación (Programa Estatal de Promoción del Talento y su Empleabilidad (FPU17/01995))
- Paula Fernández-Palanca
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.
Acknowledgements
We want to thank Patrizia Birner and Michaela Bissinger for their excellent technical support. In addition, we thank the Nikon Imaging Center at Heidelberg University (BioQuant), especially Dr. Ulrike Engel and Dr. Christian Ackermann, for supporting confocal imaging experiments. We also acknowledge the LSDF2 (URZ) for data storage. We thank Prof. Dr. Alexander Loewer (TU Darmstadt) for providing the YAP-mVenus vector and the Center for Model System and Comparative Pathology (CMCP, Heidelberg) for staining mouse tissue samples, especially Heike Conrad, Christine Schmitt, and Sarah Lammer. Tissue scanning was supported by the Tissue Bank of the National Center of Tumor Diseases (NCT, Heidelberg).
Ethics
This study was performed according to animal welfare committee of The Ministry for Environment, Agriculture, Conservation and Consumer Protection of the State of North Rhine-Westphalia (LANUV, North Rhine-Westphalia, Germany, application number: 84-02.04.2016.A279).
Senior Editor
- Mone Zaidi, Icahn School of Medicine at Mount Sinai, United States
Reviewing Editor
- Matthew A Quinn, Wake Forest School of Medicine, United States
Reviewers
- Matthew A Quinn, Wake Forest School of Medicine, United States
- Dirk Fey, University College Dublin, Ireland
Version history
- Received: March 10, 2022
- Preprint posted: April 1, 2022 (view preprint)
- Accepted: September 21, 2022
- Version of Record published: October 18, 2022 (version 1)
Copyright
© 2022, Wehling et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 695
- Page views
-
- 113
- Downloads
-
- 1
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Neuroscience
Previous research has highlighted the role of glutamate and gamma-aminobutyric acid (GABA) in perceptual, cognitive, and motor tasks. However, the exact involvement of these neurochemical mechanisms in the chain of information processing, and across human development, is unclear. In a cross-sectional longitudinal design, we used a computational approach to dissociate cognitive, decision, and visuomotor processing in 293 individuals spanning early childhood to adulthood. We found that glutamate and GABA within the intraparietal sulcus (IPS) explained unique variance in visuomotor processing, with higher glutamate predicting poorer visuomotor processing in younger participants but better visuomotor processing in mature participants, while GABA showed the opposite pattern. These findings, which were neurochemically, neuroanatomically and functionally specific, were replicated ~21 mo later and were generalized in two further different behavioral tasks. Using resting functional MRI, we revealed that the relationship between IPS neurochemicals and visuomotor processing is mediated by functional connectivity in the visuomotor network. We then extended our findings to high-level cognitive behavior by predicting fluid intelligence performance. We present evidence that fluid intelligence performance is explained by IPS GABA and glutamate and is mediated by visuomotor processing. However, this evidence was obtained using an uncorrected alpha and needs to be replicated in future studies. These results provide an integrative biological and psychological mechanistic explanation that links cognitive processes and neurotransmitters across human development and establishes their potential involvement in intelligent behavior.
-
- Computational and Systems Biology
- Neuroscience
Cerebellar climbing fibers convey diverse signals, but how they are organized in the compartmental structure of the cerebellar cortex during learning remains largely unclear. We analyzed a large amount of coordinate-localized two-photon imaging data from cerebellar Crus II in mice undergoing ‘Go/No-go’ reinforcement learning. Tensor component analysis revealed that a majority of climbing fiber inputs to Purkinje cells were reduced to only four functional components, corresponding to accurate timing control of motor initiation related to a Go cue, cognitive error-based learning, reward processing, and inhibition of erroneous behaviors after a No-go cue. Changes in neural activities during learning of the first two components were correlated with corresponding changes in timing control and error learning across animals, indirectly suggesting causal relationships. Spatial distribution of these components coincided well with boundaries of Aldolase-C/zebrin II expression in Purkinje cells, whereas several components are mixed in single neurons. Synchronization within individual components was bidirectionally regulated according to specific task contexts and learning stages. These findings suggest that, in close collaborations with other brain regions including the inferior olive nucleus, the cerebellum, based on anatomical compartments, reduces dimensions of the learning space by dynamically organizing multiple functional components, a feature that may inspire new-generation AI designs.