YAP regulates cell size and growth dynamics via non-cell autonomous mediators
Abstract
The Hippo pathway regulates organ size, regeneration, and cell growth by controlling the stability of the transcription factor, YAP (Yorkie in Drosophila). When there is tissue damage, YAP is activated allowing the restoration of homeostatic tissue size. The exact signals by which YAP is activated are still not fully understood, but its activation is known to affect both cell size and cell number. Here we used cultured cells to examine the coordinated regulation of cell size and number under the control of YAP. Our experiments in isogenic HEK293 cells reveal that YAP can affect cell size and number by independent circuits. Some of these effects are cell autonomous, such as proliferation, while others are mediated by secreted signals. In particular CYR61, a known secreted YAP target, is a non-cell autonomous mediator of cell survival, while another unidentified secreted factor controls cell size.
Introduction
Organ and tissue size are tightly regulated; they scale with body size, and in several cases can be regenerated after partial loss (Penzo-Méndez and Stanger, 2015). The size and number of cells comprising an organ together determine overall size, and these can compensate for each other (Fankhauser, 1945). In some well-studied cases adult organ repair and regeneration involves de-differentiation, followed by proliferation and re-differentiation; while in others, such as skeletal and cardiac muscle, the dedifferentiation-re-differentiation process is skipped, and organs adapt to demand by increasing cell size; or to injury by ECM deposition (fibrosis). The mass of β-cells of the pancreas can increase by augmenting both cell size (hypertrophy) or cell number (hyperplasia) (Anzi et al., 2018; Ernst et al., 2011). Apoptosis can also play an important role in adult organ size homeostasis. For example, the liver enlarges significantly after phenobarbital treatment, but is restored to its normal size by apoptosis upon drug withdrawal (Yang and Xu, 2011). These examples demonstrate that organ size is not simply set during development and passively maintained into adulthood but is actively monitored to fit physiological demands.
The liver offers quite a dramatic example of size control. Although hepatocytes are largely quiescent in adults, the liver can regenerate completely after resection of up to 70% of its original mass (partial hepatectomy). This response is tri-phasic. It starts with an increase in the size (hypertrophy) of the remaining hepatocytes, followed by a period of cell proliferation producing daughter cells of normal size, and is terminated when the original liver mass is restored (Gilgenkrantz and Collin de l'Hortet, 2018). If, in the case of a partial hepatectomy, hepatocytes are blocked from proliferation, liver size can still be restored by hypertrophy of the remaining, often polyploid cells (Diril et al., 2012). This provides a powerful example of how tissue homeostasis can be maintained by a dynamic combination of cell growth and proliferation.
While the role of cell growth in proliferation-dependent liver size homeostasis has never been explicitly tested, it is tempting to speculate that after a partial hepatectomy quiescent hepatocytes re-enter the cell cycle only after increasing in size and that perturbations that impair growth under such conditions impair regeneration (Volarevic et al., 2000). Because of the interdependence of cell size and cell cycle control it is hard to distinguish whether a genetic perturbation that affects proliferation during regeneration primarily affects the cell cycle or cell growth machinery. Though cultured cells provide almost none of the tissue-level context for studying regeneration, they still robustly control cell size by compensatory changes in growth and division (Cadart et al., 2018; Ginzberg et al., 2018).
There is increasing evidence that the Hippo pathway plays a significant role in tissue size regulation (Dong et al., 2007) as well as regeneration (Lu et al., 2018). The pathway is comprised of a cascade of intracellular kinases, activated in response to cell-cell contact. These kinases ultimately phosphorylate and inhibit the downstream transcriptional regulators YAP and TAZ, by causing their sequestration, and often degradation in the cytoplasm (Zhao et al., 2010). Upon tissue injury, the upstream kinases become inactive and YAP/TAZ translocates into the nucleus, where it mediates transcriptional changes that facilitate organ regrowth (Lee et al., 2014). YAP/TAZ activity also depends on changes in the composition and physical properties of the extracellular environment, which is believed to be independent of Hippo activity (Aragona et al., 2013; Dupont et al., 2011; Kim and Gumbiner, 2015). Activation of YAP increases cell proliferation and inhibits apoptosis in vitro and in vivo (Dong et al., 2007; Xin et al., 2011; Zanconato et al., 2015). While the phenomena have been well characterized, the exact molecular mechanism by which YAP affects these different processes and/or their interdependence remains unclear.
In this report, we investigated the mechanisms by which YAP regulates cell number (through proliferation and death) and cell growth/size in mammalian cultured cells. The term ‘cell growth’ is commonly used in two very different contexts. It is important to distinguish these processes because ‘exponential growth’ has meant either an exponential increase of cell mass in each cell cycle or an exponential increase in the number of cells per unit time. In this paper we will use the term, ‘cell growth’ exclusively to describe changes in average cell size, and ‘population growth’ exclusively to describe changes in the number of cells.
Our results suggest that YAP affects both cell size and number through distinct mechanisms. Some of these effects appear to be cell-autonomous, but notably we find strong evidence that others are driven by diffusible, non-cell autonomous factors. In particular we show that CYR61, one of the extracellular-residing proteins that is regulated both on the mRNA and protein level in response to increased YAP activity, has no effect on cell size. However, CYR61 mediates a significant part of the YAP-dependent increase in cell number by limiting apoptosis in high density cultures. We conclude that YAP stimulates the production of a suite of extracellular and cell surface proteins that act in a paracrine and autocrine manner regulating both cell number and size.
Results
Overexpression of non-phosphorylatable YAP increases the rate of cell division, cell size, and maximum cell density in HEK293 populations
To study the effect of YAP on cell size and growth kinetics, we generated Dox-inducible, isogenic FlipinTrex-HEK293 cells expressing nuclear mCherry (n-mCherry) together with either nuclear GFP (nGFP) or a GFP-tagged phosphosite mutant of YAP (YAP5SA). Unlike wildtype (wt-) YAP, YAP5SA cannot be inhibited by the upstream kinases of the Hippo pathway and therefore, is constitutively active. The fact that the cell lines are isogenic, ensures that any observed changes in size or proliferation rate between the cell lines do not arise from the genomic location of the inserts, but rather to differences in biological activity of the expressed proteins. Conveniently the n-mCherry demarks the nucleus, enabling easy counting of cell number in low magnification images (Figure 1—figure supplement 1). The nuclear marker also serves another important function, as an indirect measure of cell size (Ginzberg et al., 2018). As shown in Figure 1—figure supplement 2, this correlation is very strong when cell size is measured by the highly accurate Quantitative Phase Microscopy (QPM). The Pearson correlation coefficient was ~0.9 when measured under various conditions. With this experimental design, we could conveniently track changes in cell number and cell size in a population over time.
We found that the population growth rate in these cells decreased as cell density increased, and eventually reached zero at high cell density (Figure 1A). This is a well-known phenomenon often referred to as density-dependent growth inhibition (Stoker and Rubin, 1967) or contact inhibition; and can be described by a logistic growth expression, where the net growth rate of a population is exponential at low density and decreases to a minimum at high density (Figure 1—figure supplement 3A–C) (Verhulst, 1845). In this phenomenological description, changes in the maximum cell number are denoted by Ymax, also known as the population carrying capacity; and the growth rate is described by k (with units of reciprocal time).
We initially thought that Ymax would only be affected by the area of the culture vessel (literally the carrying capacity), while k would be a function of the concentration of growth factors in the medium. Indeed, serum concentrations above 1% changed k (Figure 1—figure supplement 3D) without changing Ymax (Figure 1—figure supplement 3E). However, Ymax was lower at serum concentration <1% suggesting that the area of the culture vessel was not the only determinant of Ymax.
A more complete description should not consider proliferation as the sole determinant of cell number but also include the cell death rate. By measuring cell death directly by YOYO-1, a cell impermeable dye that fluoresces green when bound to dsDNA in dead cells, we could show that cell death rates are low in low-density cultures (Figure 1—figure supplement 4A,B). Hence at low density the net population growth is dominated by the cell division rate. At high density, when net population growth is zero, the death rate dramatically increases (Figure 1—figure supplement 4A). Thus extra- and intra- cellular signals appear to affect k and Ymax differently, with the former predominantly reflecting cell division rates and the latter, cell death rates.
Nuclear area in low density cultures (<~15,000 cells/well in a 96-well plate) is initially constant at ~120 μm2 in control cells, after which it decreased about 33% to ~85 μm2 at maximum density (Figure 1B). Assuming the nucleus is a perfect sphere, this change in area would correspond to a 50% decrease in volume. However, actual measurements of cell volume using Coulter Counter indicate a ~ 20% decrease between low density (25% confluent) and high density (100% confluent) cultures (Figure 1—figure supplement 4C). The closer correlation between nuclear area and volume than predicted by a spherical representation of the nucleus/cell is consistent with the strong correlation we observed between dry mass and nuclear area (Figure 1—figure supplement 2A).
The expression of YAP5SA increases both k and Ymax compared to controls (Figure 1C). Furthermore, at low cell density the area of the nucleus increases 17% in the YAP5SA cells. A change of such magnitude would translate to ~30% increase in volume, which we measured by Coulter Counter to be ~25% larger (Figure 1—figure supplement 5A). However, as the number of cells continued to grow, the nuclear area of YAP5SA cells also decreased to a minimum where the difference in size compared with controls was much reduced (~10%; Figure 1D). The onset of the exponential decrease was slightly delayed when cells overexpressed YAP5SA, though cultures appear to be of comparable density (Figure 1E).
To confirm that the changes in nuclear area reflected changes in biomass and not just changes in volume, we measured protein mass directly in individual cells using a covalent dye specific for lysine groups on the protein, which is a proxy for changes in size (Kafri et al., 2013). We found that average protein content per cell was ~30% higher in low-density populations expressing YAP5SA (Figure 1F,G) than in nGFP controls. We also observed that cells seeded at four times the density had ~25% less protein (Figure 1F,G). These measurements of protein mass reinforced both of the conclusions drawn from the changes in nuclear area: (1) YAP5SA cells are larger than control cells in low density cultures, and (2) cell size decreases at higher cell population density. We ruled out that the changes were merely due to an increase in the fraction of cells in S or G2/M by simultaneously measuring DNA and protein (Figure 1—figure supplement 5B–C). The latter findings were independently confirmed in three different isogenic clones using QPM on live, attached low-density cultures (Figure 1—figure supplement 5D). Based on our measurements there is a ~ 20% increase in dry cell mass in cells expressing YAP5SA. Thus, it is fair to say that the expression of YAP5SA increases cell size in all phases of the cell cycle in a highly reproducible manner independent of the choice of clone, staining errors or detachment of cells from their substrate. However, the exact magnitude of change in size varied somewhat depending on the measurement method of choice. Further experiments using non-isogenic HEK293 cells expressing GFP-tagged YAPWT vs. nGFP demonstrated that YAPWT overexpression only increased cell size in low density cultures but not in high density cultures, a condition where YAP is excluded from the nucleus and thus transcriptionally inactive (Figure 1—figure supplement 5E). We also confirmed that knocking down YAP by 80% in WT cells by constitutively expressing YAP-targeting shRNA (YAPKD cells; Figure 1—figure supplement 5F) decreased cell size (Figure 1—figure supplement 5G) as previously reported (Hansen et al., 2015). Through the remainder of the manuscript, we mainly focus on the most dramatic changes in size caused by YAP5SA in sub-confluent cultures.
The non-cell autonomous nature of YAP-mediated changes in cell size
In isogenic Flipin-Trex cells, YAP5SA or nGFP expression is doxycycline (Dox)-inducible. Therefore, we would not expect that in the absence of Dox there would be any difference in cell size between uninduced YAP5SA, nGFP and parental cells. To our surprise in the absence of Dox the average cell mass in the three YAP5SA-expressing clones was significantly higher than in nGFP clones and the parental cells (9% vs. 20% when cells are cultured in Dox for 4 days; Figure 2A vs. Figure 1—figure supplement 5D). It is well known that Dox-inducible expression systems can be leaky, so one possible explanation for the larger size would be that the <1 % cells expressing YAP5SA in the absence of Dox could be affecting the behavior of the population as a whole. To test this possibility, we compared the size of parental FlipinTrex-293 cells (WT) when co-cultured together with cells expressing GFP-tagged YAP5SA or nGFP (Figure 2B). We found that when cultured in the same well WT cells were similar in size to the YAP5SA cells, both of which were ~30% larger than the nGFP cells and their co-cultured WTs (Figure 2C). Thus, although the expression of YAP5SA is obviously cell autonomous, its effect on size is shared with untransfected cells in the same well.
Non-cell autonomous effects could be mediated either by cell-cell contact, or by diffusible molecules (secreted into or depleted from culture medium). To investigate whether these effects required cell-cell contact, WT Flipin-Trex-293 cells were seeded on the bottom of a well, physically separated from a population of feeder cells expressing YAP5SA or nGFP in a Transwell chamber (Figure 2D). In this condition both cell types share and exchange components in the medium but do not physically contact each other. When cell size was measured by protein content, cells co-cultured with YAP5SA cells were ~50% larger compared with those co-cultured with nGFP controls (Figure 2E). This suggests that YAP’s non-cell autonomous effect on cell size is likely mediated by a soluble factor(s) produced by YAP5SA-expressing cells. Based on these results we would expect that conditioned medium from 5SA-expressing cells, when added to naïve cells, would result in a size increase, which it did (Figure 2F). Conditioned medium also increased the size of NIH 3T3 fibroblasts after a 48 hr treatment (Figure 2—figure supplement 1A,B).
The effect of YAP5SA on cell number is mediated both by cell autonomous and non-autonomous processes
To investigate whether YAP expression affects cell number as well as cell size in a non-autonomous manner, we co-cultured cells expressing n-mCherry ‘recipients’ with cells expressing either YAP5SA or nGFP only ‘feeders’ (Figure 3A). The two cell types were well-mixed keeping the total cell count constant while varying the ratio between the recipient and feeder cells. We found that recipient cells grew to a higher number (Ymax) when YAP5SA cells were in large excess (four to eight times as many), compared to being co-cultured with nGFP controls (Figure 3B,C). The population growth rate at low density (k) was not affected by the number of YAP5SA expressing cells (Figure 3D).
We suspected the difference between k and Ymax in the presence of YAP5SA was because its effect on k is cell autonomous, while its effect on Ymax is not. If this hypothesis were true, k would be lower in YAPKD cells irrespective of what cells they were co-cultured with. In contrast Ymax should increase in YAPKDs when co-cultured with YAP5SA feeders. Indeed, k was lower in YAPKDs vs controls, and co-culturing either cell type with YAP5SA cells had little effect on k, while significantly increasing Ymax (Figure 3E–G). These results are most easily explained if the effect of YAP on k was indeed cell autonomous, whereas its effect on Ymax was non-cell autonomous. Furthermore, we show that treating serum-starved FlipIn-Trex-293 or 3T3 cells with 5SA-conditioned medium increased the number of cells in dense cultures vs. the addition of conditioned medium derived from nGFP cells (Figure 3H and Figure 1—figure supplement 5D,E). These results also support a YAP-dependent, non-cell autonomous mode of regulating Ymax.
Many proteins affected by YAP5SA overexpression are secreted or membrane proteins
We have tested the notion that these non-cell autonomous effects are mediated by secreted proteins that are transcriptional targets of YAP. Using RNA-Seq and mass spectrometry we compared the transcriptional and corresponding protein expression changes in FlipinTrex-293 cells in response to YAP5SA expression (Figure 4A, Supplementary file 1). We concluded that extracellular proteins were a large portion of those expressed in response to YAP5SA (eight fold enrichment relative to all detected proteins, p-value<10−12; Figure 4B). While in some cases the change in RNA and protein levels were similar, in many other cases they were not. For example, CTGF, CYR61 and AMOTL2, which are canonical YAP targets, were elevated at the protein and mRNA levels, while we see an almost four-fold decrease in LOX and TNC proteins, with no change on the mRNA level (Figure 4A). CYR61 and CTGF are members of the CCN family of matricellular proteins, non-structural proteins that reside in the extracellular matrix and are thought to regulate cellular processes like migration, apoptosis, division and differentiation, involving the coordination of several signaling pathways (Holbourn et al., 2008). Though both genes have previously been used as robust transcriptional readouts for YAP activity (Dupont et al., 2011; Yimlamai et al., 2014), their contribution to the YAP-dependent changes in cell behavior has not been clarified. Specifically, their role in coordinating YAP-dependent cell proliferation, apoptosis and size control has not been explicitly studied. We also considered whether Amphiregulin, a secreted molecule related to EGF, that had previously been identified as a potential mediator of YAP-dependent non-cell autonomous growth in 3D cultures of MCF10 cells (Zhang et al., 2009) was involved in YAP-dependent size or proliferation control in our cells. In these experiments, Amphiregulin RNA was expressed at a very low level (FPKM <1 compared with FPKM >100 for CYR61); and was not identified by mass spectrometry. By contrast we identified 3 CYR61 and 6 CTGF peptides.
The YAP target, CYR61, affects cell number but not cell size
To test whether CYR61 and/or CTGF mediated any of the non-cell autonomous changes in cell size or number, we used neutralizing antibodies to inhibit their extracellular activity in co-cultures of n-mcherry recipient cells and YAP5SA cells. The identity of the recombinant proteins and the specificity of the antibodies we used were demonstrated by their mutual recognition in Figure 5—figure supplement 1. Neutralizing CTGF had a small effect on cell number, mildly decreasing k (~12%) while mildly increasing Ymax (~12%) of the n-mCherry cells in the co-cultures (Figure 5A–C). On the other hand, neutralizing CYR61 significantly decreased Ymax (up to ~50%) but increased k (up to ~200%) in a dose-dependent manner (Figure 5D–F). Our RNA-Seq data suggest that the amount of CYR61 (FPKM >100) transcribed in these cultures is greater than CTGF (FPKM <50). Consistent with the difference in expression levels, increasing the concentration of antibody had a much larger effects for anti-CYR61, as compared to CTGF (Figure 5A,D). These findings support a major role for CYR61 in the non-cell autonomous, YAP-mediated increase in Ymax. However, the increase in the rate of population growth observed at low cell density appears not to be mediated by CYR61, consistent with the measurements in Figure 3 that show that the dominant effects on k are cell autonomous.
To test whether YAP expression in the responding cells is required for the CYR61-dependent increase in Ymax, we knocked down YAP and added CYR61. As Shown in Figure 5G–L YAPKDs as well as controls responded by increasing Ymax by ~30% and 20% respectively, but decreased k by about 10% and 50% (Figure 5G–L), confirming that CYR61’s non-cell autonomous influence on cell number – specifically the carrying capacity (Ymax) – does not require active YAP in the responding cells. These results were therefore consistent with the co-culture experiments (Figure 3E–G). Finally, neither CTGF or CYR61 neutralization nor addition affected nuclear area (Figure 5M–P); hence some factors other than CTCF or CYR61 must be responsible for the non-autonomous effects of YAP on cell size.
Cyr61 increases Ymax by reducing apoptosis
Ymax is clearly different between nGFP and 5SA cells (Figure 1C), although the distribution of cells in the cell cycle suggests no difference in proliferation rates between both cell types at high density (Figure 6A). We hypothesized the striking effect of CYR61 on Ymax is most easily explained by an effect on cell death, which occurs at higher rates in high vs. low density cultures (Figure 1—figure supplement 1–4A,B). Consistent with that, in the presence of anti-Cyr61 antibody, there is a three-fold increase in cell death in cultures of Flipin-Trex-293 cells, which can be reduced to control levels upon addition of an anti-apoptotic peptide (Q-VD-OPH; Figure 6B). We also observed that cells treated with anti-CYR61 were mostly all dead by the end of 6 days (Figure 6C). Cells treated with CYR61 or IgG were almost 100% alive with many still cycling. We conclude that high levels of CYR61 reduce density-dependent apoptotic cell death in YAP5SA cells.
Discussion
In the early studies of cell growth, the focus was on unicellular organisms or cells grown at low density. Studying cell growth at high density meant having to deal with complex interactions of the individual components, the space they occupy, and effects on the solvent or the medium, all of which are less confounding at low density. Nevertheless, interest in growth at high cell density was piqued when biologists made 2D cultures of animal cells on petri dishes. For example, chicken embryo fibroblasts grown in such a manner generate a monolayer and stop growing (Weinberg, 2007). It is unlikely that this was simply depletion of the medium, for when Rous Sarcoma Virus was added the quiescent cells were transformed and proliferated, overgrowing the monolayer and producing three dimensional ‘foci’ on top of the monolayer (Groupe and Manaker, 1956). There are two phenomena here to consider, the still unexplained tendency of cells to stop growing upon crowding and the tendency of a few infected or transformed cells to overcome whatever inhibitions were observed in the uninfected monolayers. With the passage of 50 years these phenomena have neither been completely resolved nor completely forgotten. The problem of cell interactions is still a serious topic in cancer (Kamińska et al., 2015) and more recently has emerged in many other biological settings, such as the behavior of cells in regeneration or in stem cell niches (Lane et al., 2014).
More interest in the reciprocal interactions of proliferating and resident cell populations was stirred by the discovery of the Hippo pathway, originally discovered in Drosophila, and later shown to be activated in high-density cultures. The major effector of the pathway is the transcription regulator YAP, which is formally considered an oncogene, because it stimulates growth (Dong et al., 2007). Over-activating YAP results in a breakdown of some of the usual barriers to growth in complex tissues and in high density cell cultures (Zhao et al., 2007). An interesting feature of the pathway is its sensitivity to alterations in the nature of the substratum, such as ECM composition and stiffness (Dupont et al., 2011; Kim and Gumbiner, 2015). It is still generally assumed that the aggressive behavior of the cells is governed in a cell autonomous manner – like the original Rous transformation mentioned above. Another, and more unique feature of YAP is its control of organ size. Organ size is determined by a tight coordination of cell death, proliferation and sometimes differentiation as well as cell size. Organ size control is best studied in vivo where the complexity of cell-cell and cell-matrix interactions are better maintained, but this complexity makes teasing apart underlying mechanisms challenging. Cultured cells can be used to model some of the aspects of organ size in the whole animal but lacks the more complex context. Nevertheless, as cultured systems become more complex, there is an expectation that mechanistic features of organ and tissue growth can be recognized and understood.
In this paper we have extended observations to cell autonomous and non-cell autonomous signaling downstream of YAP in cultured cells. In this system we have examined the growth of cells at different cell densities; the effects of cell proliferation and cell death; and the combined effects of proliferation and cell size. By such measurements we were able to gain a clearer understanding of YAP’s role in regulating these phenomena on a cellular as well as population level. We found that the growth rate of populations, k, at low density, was elevated significantly when cells expressed a constitutively active form of YAP, YAP5SA. compared to wild type (Figure 1B) cultures. The effect on k was cell autonomous and driven predominantly by a change in proliferation rates. This increase in proliferation was accompanied by an increase in the fraction of cells in S-G2-M phases of the cell cycle and may be explained by previous reports that there is a YAP-mediated increase in the cell cycle regulators CyclinE and Cyclin D, and consequently a higher rate of entry into S-phase (Dong et al., 2007; Shu and Deng, 2017; Tapon et al., 2002).
Even more obvious was YAP5SA’s effect on the number of cells at high density, where the carrying capacity, Ymax, increased about 35%. The effect of YAP on Ymax was demonstrably non-cell autonomous and mediated by CYR61, which increases Ymax even in the absence of endogenous YAP indicating that it is not simply an autocrine mediator. CYR61 is one of the strongest genes we find differentially expressed on both the protein and mRNA level. It is a common YAP target that is a member of the CCN group of matricellular proteins, known to be involved in tissue remodeling (Chen et al., 2001), inflammation and embryogenesis (Katsube et al., 2009; Latinkic et al., 2003). In our experiments, the effects of CYR61 on cell number at high density are largely a result of a reduction in apoptosis, while it seems to suppress proliferation at low density. Previously, CYR61 was shown to reduce apoptosis in MCF7 cells (Lin et al., 2004) and is one of the earliest induced genes in response to a partial hepatectomy, which also involves/requires YAP (Grijalva et al., 2014; Lu et al., 2018; Su et al., 2002). Only recently has CYR61 activity in hepatocytes been implicated in YAP-mediated liver fibrosis by activating macrophages (Mooring et al., 2020), suggesting that CYR61 is more than a YAP reporter in that context, and is of importance in mediating non-cell autonomous YAP activity in vivo. CYR61 was also reported to play a role in inhibiting hepatocyte proliferation in response to carcinogen-mediated injury (Chen et al., 2016). Altogether, this suggests that CYR61 plays various roles in cell culture as well as in vivo, which could depend on the cell type as well as the particular context in which CYR61 is activate. This is particularly interesting in the whole-organ context, where the coordination of interactions between cell types is often necessary and would be expected to depend on cell non-autonomous signaling. The present study involved only a single type of cell examined in two contexts as responders and signalers. Exploring the role of YAP and its extracellular targets in mixed cultures of different cell types could provide additional insight into its role in complex tissues.
Finally, cell size, also seems to be regulated non-cell autonomously upon YAP expression. We applied several modes of measurement: nuclear area, volume, protein content and dry mass, all of which are often used as metrics to describe cell size, and they all show that expressing YAP5SA increases cell size between ~17–30%. We believe the changes in nuclear area and dry mass (17% vs. 20% respectively) are the most reliable and are also well-correlated. We tested whether the two extracellular proteins most affected by YAP expression, CYR61 and CTGF, have an effect on WT cell size. Neither the depletion of either protein nor the addition of the recombinant proteins affected cells size. We did not find evidence that the classic size regulating growth factors IGF1/2 are produced to any appreciable degree by our cells. There are, however, many candidate secreted proteins that are stimulated by YAP and which should be tested in the future.
The YAP targets, CTGF and CYR61, both members of the CCN family of matricellular proteins, have opposing effects on cell number in co-cultures with YAP5SA expressing cells. The levels of CTGF are small compared with CYR61 (Figure 5) which we believe explains why the observed non-cell autonomous effects on population growth are mainly those mediated by CYR61. CCN proteins, unlike typical growth factors, are not known to bind to any receptor tyrosine kinases, but rather to integrins and other growth factors, thereby mediating the activity of those conventional growth factors (Leask and Abraham, 2006). With low serum in the medium (<1%), Ymax decreases (Figure 1—figure supplement 3E). Size is also affected by the level of serum (Figure 1—figure supplement 6B), suggesting that another matricellular factor might mediate YAP-dependent changes in cell size.
Work by Tian and colleagues (Tian et al., 2013) had previously demonstrated that a high molecular weight hyaluronic acid produced by naked mole rate fibroblasts contributes to early contact inhibition; the equivalent of a decrease in the population’s carrying capacity. This form of hyaluronic acid acts via the upstream Hippo component NF2, which inhibits YAP activity. Combined with our results on secreted factors, such findings suggest that extracellular cues are an important part of the regulation of both cell and organ (tissue) size. The non-cell autonomous nature of YAP signaling would allow local signals to affect remote cells of the same or different type as is expected to occur in regeneration. In other words, the activation of YAP in a small fraction of cells at the site of injury could through secreted signals engage a wider portion of the tissue or even a whole organ to participate in tissue repair. The receiving cells might also secrete signals as part of their response, further reinforcing the signals and allowing them to propagate farther than would be possible by diffusion alone.
Materials and methods
Cloning
Request a detailed protocolpEGFP-C3-YAP5SA plasmids were generated from the pEGFP-C3-hYAP1 plasmid gifted from Marius Sudol (Addgene plasmid # 17843; RRID:Addgene_17843 Basu et al., 2003). YAP amino acids S61, S109, S127, S164, S381 were all mutated to Alanine using the QuikChange II XL Site-Directed Mutagenesis Kit (Agilent). The following nuclear localization signal atggatccgaagaaaaaacgtaaaggccgtatggatccgaagaaaaaacgtaaagg-ccgt was appended to the 3’-end of GFP in the pAcGFP1-C3 vector (Clontech) to ensure its nuclear localization.
Cell line generation
Request a detailed protocolAll cell lines were maintained in DMEM (Gibco, 10569010) supplemented with 10% FBS (Gibco), 1% P/S (Gibco) unless otherwise stated. Dox-inducible HEK293 cell lines of nGFP or GFP-tagged YAP5SA were generated by sub-cloning either construct into the pcDNA5/FRT expression vector (Invitrogen) then co-transfecting the generated plasmids with the pOG44 Flp-Recombinase Expression Vector to ensure recombination into the genome of Flipin-TRex −293 cells (Invitrogen). Selection with Hygromycin and single clone isolation was followed by a test for Zeocin-sensitivity to ensure proper integration into the genome. Cells were maintained in Doxycycline unless otherwise specified.
YAPsh and control constructs were generously provided by Taran Gujral (Gujral and Kirschner, 2017). After transfection of HEK293 cells (ATCC Cat# CRL-1573, RRID:CVCL_0045) using PEI-40 and antibiotic selection in Puromycin for two weeks, isolated single colonies were observed and were individually isolated. To label cells with nuclear localized mCherry we used Lentiviral infection. First nuclear-mCherry was excised from pBRY-nuclear mCherry-IRES-PURO (Addgene plasmid #52409, RRID:Addgene_52409) by digestion with EcoRI and NotI and ligated into lentiviral expression plasmid pLVX-Ef1α-N1-mCherry. Lentiviruses were then produced via Fugene six co-transfection of psPAX2 (Addgene plasmid #12260, RRID:Addgene_12260) and pMD2.G (Addgene plasmid #12259, RRID:Addgene_12259) into 293 T cells ((ATCC Cat# CRL-11268, RRID:CVCL_1926). The harvested lentivirus was then used to infect the desired cell line. One week after cells were infected with Lentivirus, they were FACS sorted for mCherry expression.
Transwell assays
Request a detailed protocolTranswell assays were done using 75 mm Transwell plates with a 0.4 µm Pore Polycarbonate Membrane Insert (Corning). Recipient cells were cultured on the bottom of the plates while feeder cells were cultured on the insert for 72 hr.
Antibodies and recombinant proteins
Request a detailed protocolNeutralization antibodies against CYR61 (mainly NOVUS biological NB100-356, RRID:AB_10000986 raised against a peptide of amino acids 150–250 of human CYR61, validated independently with Abnova PAB12526, RRID:AB_10550916 raised against a peptide of amino acids 250–350 of human CYR61) or CTGF (Peprotech 500-P252, RRID:AB_1268211) were used at the indicated concentrations in the respective figures. A polyclonal Rabbit IgG Antibody was used as a control (Protein tech; 30000–0-AP, RRID:AB_2819035). Recombinant proteins CTGF (Peprotech), CYR61 (Peprotech) and Insulin (Peprotech) were used at the concentrations and under the conditions indicated in the respective figure legends.
Conditioned medium production
Request a detailed protocolFlipinTrex-cells were cultured in Nunc TripleFlasks with DMEM plus 10% FBS and 1% P/S until they reached 100% confluence. Cells were then washed in PBS and kept in medium with 0.5% FBS and Dox for 6 days. Medium was collected, spun down and the supernatant then filtered on a 0.2 μm membrane to remove cell debris. The medium was concentrated using 3 kD cut-off Centricon Plus-70 filter devices (Millipore). The concentrate was kept at −20°C until further use. For protein concentration estimation, the concentrate was diluted to within working range before a BCA assay (Thermo Scientific) was performed against a standard.
Incycyte-based image acquisition and automated image analysis for nuclear segmentation
Request a detailed protocolImages from four different regions per well were acquired at regular intervals using a 10x lens (Nikon, MRH00101) from cells cultured on 96-well plates (Eppendorf) in an Incucyte Zoom apparatus installed in a regular tissue culture incubator. After image acquisition, the red signal was used for nuclear segmentation using Incucyte ZOOM’s automated image analysis software. Images from all wells and across the entire time course were segmented using the same parameters. The number of red objects was used to represent the total number of nuclei. The average area estimated from the segmented objects in the red channel represent the average nuclear area.
Quantitative phase microscopy (QPM)
Request a detailed protocol50,000 Flipin-Trex-293 cells were seeded on a Fibronectin-coated (50 ug/ml) glass-bottom 6-well plate (Cellvis) ~18 hr before imaging. This ensured that cells were at a low-enough density to allow high-quality segmentation of ~4000 single cells per well. Cells were cultured with 2 μg/ml Hoechst for 30’ before imaging. Imaging and analysis was done as described in Liu et al. (2019). Briefly Cells were maintained at 36°C in 5% CO2. The phase images were taken using a 10x lens of N.A. 0.30 (Nikon, MRH10101) using the SID4BIO camera (Phasics) integrated into a Nikon Eclipse Ti microscope with a halogen lamp as the transmitted light source. A C-HGFI mercury lamp was used for fluorescence illumination. All the image processing was done using a custom MATLAB script executed on a high-performance computer cluster.
Data plotting and fitting
Request a detailed protocolData were imported into Prism7 (Graph Pad) where their visualization, fitting and statistical analysis was done.
Coulter counter measurement
View detailed protocolCell volume was measured on a Coulter Counter (Beckman Coulter). For the measurement cells were trypsinized to single cells then resuspended in 10 ml Isoton diluent buffer (Beckman Coulter).
Total protein/DNA staining
Request a detailed protocolCells were trypsinized to a single-cell suspension. 10^6 cells were counted per condition, thoroughly washed in PBS and fixed for 20 min in ice cold, 4% PFA solution. Cells were then washed in PBS and permeabilized using 0.5% Triton-X for 10 min at room temperature. After repeated washes, cells were kept on ice with a 10 μg/ml Hoechst 3342 and/or 0.4 μg/ml of Alexa Fluor 647 Succinimidyl-Ester (Invitrogen) in a light protected tube for 45 min before they were analyzed on an LSRII flow cytometer (BD Biosciences).
RNA-Sequencing and data analysis
Request a detailed protocolFlipin-Trex-293 cells expressing Tet-inducible YAP5SA or nGFP were grown to confluence in medium supplemented with Tet-approved FBS (Takara) supplemented with 0 or 5 ng/ml of Dox. Total RNA was harvested and isolated using Qiagen’s RNeasy mini kit. Samples were further processed at Harvard’s Biopolymer facility as follows: After polyA enrichment of mRNA using Takara’s PrepX PolyA mRNA Isolation Kit, libraries were prepared using Illumina’s PrepX RNA-Seq Library reagent kit on the Apollo 324 system. Samples were pooled then single-end sequenced on the Illumina NextSeq 500 platform. Sequencing alignment was done using STAR (Dobin et al., 2013) and transcriptome assembly and differential expression was done using Cufflinks (Trapnell et al., 2012) on a computer cluster. Each comparison group contained at least two independent biological repeats. Data is deposited in GEO under accession number: GSE123296.
Mass spectrometry and data analysis
Request a detailed protocolCells were kept under similar conditions as those used for RNA-sequencing. Cellular protein lysate was collected and snap frozen in urea-free lysis buffer supplemented with protease and phosphatase inhibitors (Roche). The supernatant was collected and filtered on a low protein-binding PVDF 0.22 μm membrane to eliminate cell debris, then snap frozen with protease inhibitors (Roche) for further processing. Lysates were reduced with 5 mM DTT, alkylated with 15 mM N-ethylmaleimide for 30 min in the dark. Alkylation reactions were quenched with 50 mM freshly prepared DTT, and proteins precipitated by methanol/chloroform precipitation. Digests were carried out in 200 mM EPPS pH 8.5 in presence of 2% acetonitrile (v/v) with LysC (Wako, 2 mg/ml, used 1:75) for 3 hr at room temperature and the subsequent addition of trypsin (Promega #V5111, stock 1:75) over night at 37°C.
Missed cleavage rate was assayed from a small aliquot by mass spectrometry. For whole proteome analysis, digests were directly labeled with TMT reagents (Thermo Fisher Scientific). Labeling efficiency and TMT ratios were assayed by mass spectrometry, while labeling reactions were stored at −80°C. After quenching of TMT labeling reactions with hydroxylamine, TMT labeling reactions were mixed, solvent evaporated to near completion, and TMT labeled peptides purified and desalted by acidic reversed phase C18 chromatography. Peptides were then fractionated by alkaline reversed phase chromatography into 96 fractions and combined into 24 samples. Before mass spectrometric analysis, peptides were desalted over Stage Tips (Rappsilber et al., 2003).
Data were collected by a MultiNotch SPS MS3 method (McAlister et al., 2014) using an Orbitrap Lumos mass spectrometer (Thermo Fisher Scientific) coupled to a Proxeon EASY-nLC 1000 liquid chromatography (LC) system (Thermo Fisher Scientific). The 100 μm inner diameter capillary column used was packed with C18 resin (Accucore 2.6 μm, 150 Å, Thermo Fisher Scientific). Peptides of each fraction were separated over 3–5 hr acidic acetonitrile gradients by LC prior to mass spectrometry (MS) injection. The first sequence scan was an MS1 spectrum (Orbitrap analysis; resolution 120,000; mass range 400–1400 Th). MS2 analysis followed collision-induced dissociation (CID, CE = 35) with a maximum ion injection time of 150 ms and an isolation window of 0.7 Da. To obtain quantitative information, MS3 precursors were fragmented by high-energy collision-induced dissociation (HCD) and analyzed in the Orbitrap at a resolution of 50,000 at 200 Th.
MS3 injection time for phosphopeptides was 150 and 200 ms at a resolution of 50,000. Further details on LC and MS parameters and settings used were described recently (Paulo et al., 2016a).
Peptides were searched with a SEQUEST-based in-house software against a human database with a target decoy database strategy and a false discovery rate (FDR) of 1% set for peptide-spectrum matches following filtering by linear discriminant analysis (LDA) and a final collapsed protein-level FDR of 1%. Quantitative information on peptides was derived from MS3 scans. Quant tables were generated requiring an MS2 isolation specificity of >75% for each peptide and a sum of TMT s/n of >0 over all channels for any given peptide and exported as TAB-separated raw text data. The peptide-level sums of TMT s/n signal were integrated into protein-level and converted into relative ratio with 90% confidence intervals using a recently developed method BACIQ - a Bayesian approach to confidence inference for quantitative proteomics. Briefly, this approach reconciles disagreement across multiple peptides and differences in absolute levels of peptide signal, reporting the ratios across conditions for a given protein, modeling quantitative proteomics measurement as a random variable distributed according to hierarchical Dirichlet-Multinomial probability distribution. This method accounts for the absolute level of the peptide signal in a way calibrated for a particular MS instrument. More specifically we used a charge conversion value of 1.7 as previously fitted by us for this instrument and mass resolution (Peshkin et al., 2019). Details of the TMT intensity quantification method and further search parameters applied were described in Paulo et al. (2016b).
Gene ontology analysis was done using DAVID (Huang et al., 2009) using all identified proteins as background.
Data availability
Sequencing data has been deposited in GEO under accession code GSE123296.
-
NCBI Gene Expression OmnibusID GSE123296. Induction of YAP5SA-GFP vs nGFP expression in FlipinTrex293 cells.
References
-
The angiogenic factor Cyr61 activates a genetic program for wound healing in human skin fibroblastsJournal of Biological Chemistry 276:47329–47337.https://doi.org/10.1074/jbc.M107666200
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Mechanisms in the adaptation of maternal β-cells during pregnancyDiabetes Management 1:239–248.https://doi.org/10.2217/dmt.10.24
-
Understanding liver regeneration: from mechanisms to regenerative medicineThe American Journal of Pathology 188:1316–1327.https://doi.org/10.1016/j.ajpath.2018.03.008
-
Dynamic alterations in hippo signaling pathway and YAP activation during liver regenerationAmerican Journal of Physiology-Gastrointestinal and Liver Physiology 307:G196–G204.https://doi.org/10.1152/ajpgi.00077.2014
-
The CCN family of proteins: structure-function relationshipsTrends in Biochemical Sciences 33:461–473.https://doi.org/10.1016/j.tibs.2008.07.006
-
The role of the cell-cell interactions in Cancer progressionJournal of Cellular and Molecular Medicine 19:283–296.https://doi.org/10.1111/jcmm.12408
-
Role of CCN, a vertebrate specific gene family, in developmentDevelopment, Growth & Differentiation 51:55–67.https://doi.org/10.1111/j.1440-169X.2009.01077.x
-
Adhesion to fibronectin regulates hippo signaling via the FAK-Src-PI3K pathwayThe Journal of Cell Biology 210:503–515.https://doi.org/10.1083/jcb.201501025
-
Modulating the stem cell niche for tissue regenerationNature Biotechnology 32:795–803.https://doi.org/10.1038/nbt.2978
-
All in the CCN family: essential matricellular signaling modulators emerge from the bunkerJournal of Cell Science 119:4803–4810.https://doi.org/10.1242/jcs.03270
-
YAP and TAZ regulate skin wound healingJournal of Investigative Dermatology 134:518–525.https://doi.org/10.1038/jid.2013.339
-
Cyr61 expression confers resistance to apoptosis in breast Cancer MCF-7 cells by a mechanism of NF-kappaB-dependent XIAP up-regulationJournal of Biological Chemistry 279:24015–24023.https://doi.org/10.1074/jbc.M402305200
-
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
-
A Triple Knockout (TKO) Proteomics Standard for Diagnosing Ion Interference in Isobaric Labeling ExperimentsJournal of The American Society for Mass Spectrometry 27:1620–1625.https://doi.org/10.1007/s13361-016-1434-9
-
Organ-Size regulation in mammalsCold Spring Harbor Perspectives in Biology 7:a019240.https://doi.org/10.1101/cshperspect.a019240
-
Bayesian confidence intervals for multiplexed proteomics integrate Ion-statistics with peptide quantification concordanceMolecular & Cellular Proteomics 18:2108–2120.https://doi.org/10.1074/mcp.TIR119.001317
-
Differential regulation of cyclin E by Yorkie-Scalloped signaling in organ developmentG3: Genes|Genomes|Genetics 7:1049–1060.https://doi.org/10.1534/g3.117.039065
-
Recherches mathématiques sur la loi d’accroissement de la populationNouv Mém Acad R. Sci B.-Lett. Brux pp. 1–41.
-
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
Article and author information
Author details
Funding
Harvard Medical School (Dean's Innovation Grant)
- Marc W Kirschner
National Institute of General Medical Sciences (GM026875)
- Marc W Kirschner
National Institute of Child Health and Human Development (HD091846)
- Marc W Kirschner
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Mike Gage for his meticulous assistance with the preparation of the protein samples for mass spectrometry and his overall management of project resources. We thank Jodene Moore of the SysBio FACS facility for her indispensable help setting up FACS analyses on the facility machines. We also thank the BPF Next-Gen Sequencing Core Facility at Harvard Medical School for their expertise and instrument availability that supported this work. We are grateful to the Nikon Imaging Center at Harvard Medical School for sharing resources; and to the Research Computing Group for support with cluster computing for image processing, RNAseq data processing and data transfer. We are grateful for funding from HMS through the Dean’s Innovation grant. Finally, we are grateful to Shangqin Guo and Amaleah Hartman for the stimulating discussions and thoughtful critique of this paper.
Copyright
© 2020, Mugahid 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
-
- 6,282
- views
-
- 762
- downloads
-
- 35
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
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
-
- Cell Biology
- Developmental Biology
The bone-resorbing activity of osteoclasts plays a critical role in the life-long remodeling of our bones that is perturbed in many bone loss diseases. Multinucleated osteoclasts are formed by the fusion of precursor cells, and larger cells – generated by an increased number of cell fusion events – have higher resorptive activity. We find that osteoclast fusion and bone resorption are promoted by reactive oxygen species (ROS) signaling and by an unconventional low molecular weight species of La protein, located at the osteoclast surface. Here, we develop the hypothesis that La’s unique regulatory role in osteoclast multinucleation and function is controlled by an ROS switch in La trafficking. Using antibodies that recognize reduced or oxidized species of La, we find that differentiating osteoclasts enrich an oxidized species of La at the cell surface, which is distinct from the reduced La species conventionally localized within cell nuclei. ROS signaling triggers the shift from reduced to oxidized La species, its dephosphorylation and delivery to the surface of osteoclasts, where La promotes multinucleation and resorptive activity. Moreover, intracellular ROS signaling in differentiating osteoclasts oxidizes critical cysteine residues in the C-terminal half of La, producing this unconventional La species that promotes osteoclast fusion. Our findings suggest that redox signaling induces changes in the location and function of La and may represent a promising target for novel skeletal therapies.
-
- Cell Biology
The β-catenin-dependent canonical Wnt signaling is pivotal in organ development, tissue homeostasis, and cancer. Here, we identified an upstream enhancer of Ctnnb1 – the coding gene for β-catenin, named ieCtnnb1 (intestinal enhancer of Ctnnb1), which is crucial for intestinal homeostasis. ieCtnnb1 is predominantly active in the base of small intestinal crypts and throughout the epithelia of large intestine. Knockout of ieCtnnb1 led to a reduction in Ctnnb1 transcription, compromising the canonical Wnt signaling in intestinal crypts. Single-cell sequencing revealed that ieCtnnb1 knockout altered epithelial compositions and potentially compromised functions of small intestinal crypts. While deletion of ieCtnnb1 hampered epithelial turnovers in physiologic conditions, it prevented occurrence and progression of Wnt/β-catenin-driven colorectal cancers. Human ieCTNNB1 drove reporter gene expression in a pattern highly similar to mouse ieCtnnb1. ieCTNNB1 contains a single-nucleotide polymorphism associated with CTNNB1 expression levels in human gastrointestinal epithelia. The enhancer activity of ieCTNNB1 in colorectal cancer tissues was stronger than that in adjacent normal tissues. HNF4α and phosphorylated CREB1 were identified as key trans-factors binding to ieCTNNB1 and regulating CTNNB1 transcription. Together, these findings unveil an enhancer-dependent mechanism controlling the dosage of Wnt signaling and homeostasis in intestinal epithelia.