Abstract
Dynamic changes in mechanical microenvironments, such as cell crowding, regulate lineage fates as well as cell proliferation. Although regulatory mechanisms for contact inhibition of proliferation have been extensively studied, it remains unclear how cell crowding induces lineage specification. Here we found that a well-known oncogene, ETS variant transcription factor 4 (ETV4), serves as a molecular transducer that links mechanical microenvironments and gene expression. In a growing epithelium of human embryonic stem cells, cell crowding dynamics is translated into ETV4 expression, serving as a pre-pattern for future lineage fates. A switch-like ETV4 inactivation by cell crowding derepresses the potential for neuroectoderm differentiation in human embryonic stem cell epithelia. Mechanistically, cell crowding inactivates the integrin–actomyosin pathway and blocks the endocytosis of fibroblast growth factor receptors (FGFRs). The disrupted FGFR endocytosis induces a marked decrease in ETV4 protein stability through ERK inactivation. Mathematical modelling demonstrates that the dynamics of cell density in a growing human embryonic stem cell epithelium precisely determines the spatiotemporal ETV4 expression pattern and, consequently, the timing and geometry of lineage development. Our findings suggest that cell crowding dynamics in a stem cell epithelium drives spatiotemporal lineage specification using ETV4 as a key mechanical transducer.
Similar content being viewed by others
Main
Gastrulation is an early developmental event to derive the three embryonic germ layers. In amniotes, the pre-gastrulation embryo goes through an evolutionarily conserved morphogenetic process where it forms a single-layered epithelial sheet known as the epiblast1. The evolutionary conservation of the single-layered epithelium suggests that this morphological structure is the prerequisite for forming the three embryonic germ layers1. This epithelium stably persists for about one week in primate post-implantation embryos, during which it undergoes 10- to 20-fold size expansion and acquires mature differentiation potential2,3. However, the direct contribution of epithelial expansion to lineage determination remains uncertain.
Self-organization is a cellular process that spontaneously creates complex structures with no particular pre-pattern4,5. In human embryonic stem cell (hESC)-derived gastruloid models based on micropatterning technology, a differentiating hESC colony forms an ordered structure of germ layers along the radial axis with ectoderm in the centre and mesendoderm in the periphery6,7. In past decades, a significant emphasis has been placed on diffusible signalling factors to explain self-organization8. Accordingly, the gradients of BMP4 and Activin–Nodal signalling were suggested to drive the spatial derivation of multiple germ layers within an hESC colony6,7. Nonetheless, considering the evolutionary conservation of epithelial structures in pre-gastrulating epiblasts, we hypothesized that unique morphological features related to the epithelium are involved in lineage specification.
In a single-layered epithelium, active proliferation in a confined environment induces local cell crowding that inhibits epithelial cell proliferation and induces cell extrusion9. YAP/TAZ and PIEZO1 were identified as crucial mediators linking cell crowding to cellular responses10,11,12,13,14,15. Furthermore, growing evidence suggests that cell crowding plays an essential role in cellular differentiation in the epidermis and developing zebrafish hearts16,17. Together with the critical roles of YAP/TAZ and PIEZO1 in stem cell differentiation18,19,20,21,22,23, these findings suggest local cell crowding as a key cellular mechanism through which a new cell fate is generated within a seemingly homogeneous stem cell population. Nonetheless, given the diversity of mechanical stimuli encountered by stem cells during embryonic development, it remains crucial to identify additional factors that transduce mechanical signals to gene expression.
In this work, we identified an ultrasensitive mechanical transducer, ETV4, and its upstream mechanotransduction pathway in hESCs. This discovery unveils the mechanism by which epithelial crowding regulates spatiotemporal lineage derivation.
Results
Regulation of ETV4 expression by mechanical cues
To identify transcriptional regulators that link cell crowding to transcriptional responses in hESCs, we manipulated local cell crowding by limiting the cell adhesive area and performed RNA sequencing (RNA-seq) (Fig. 1a, Extended Data Fig. 1b, and Supplementary Table 1). Cell crowding was verified by measuring the average area of cells covering the substrate (hereafter referred to as cell area; Extended Data Fig. 1a). Gene ontology (GO) analysis with total differentially expressed genes (DEGs) (|log2(fold change)| > 0.5, adjP < 0.05) revealed the top 10 GO terms related to embryonic development and differentiation (Extended Data Fig. 1c)24, supporting the relevance of cell crowding in developmental processes. Contractile actomyosin bundles represent a reliable indicator of cellular mechanical stress levels. Gene set enrichment analysis (GSEA) revealed a significant reduction in genes related to actomyosin in crowded hESCs (Extended Data Fig. 1d)25, confirming decreased mechanical stress. These results were validated by immunostaining for phosphorylated myosin light chain (pMLC) that marks contractile actomyosin bundles (Extended Data Fig. 1e).
To infer transcriptional regulators, we applied a recently developed analytic tool called Lisa (Landscape In Silico Deletion Analysis) that uses chromatic profile data26. Lisa analysis with top 200 downregulated DEGs in crowded hESCs (ordered by fold change) revealed 40 transcriptional regulators with significant P values (P < 0.01) and high expression in hESCs (fragments per kilobase of transcript per million mapped reads (FPKM) > 10; Fig. 1b and Supplementary Table 2). TEAD2 and SRF, previously known to be related to mechanotransduction, are included in the list. Protein function annotation with the selected transcription factors (TFs) using InterPro revealed seven protein domains with significant P values (adjP < 0.05). Among those identified, we found two terms related to PEA3-type ETS-domain TFs (Fig. 1b).
The PEA3 family consists of ETV1 (also known as ER81), ETV4 (also known as PEA3) and ETV5 (also known as ERM)27. Immunostaining validated the high expression of all PEA3 family TFs in hESCs and cell-crowding-induced reduction in protein expression (Fig. 1c–e). Among the PEA3 family TFs, we focused on ETV4 for its highest expression level in hESCs (Extended Data Fig. 1f). To test if cell-crowding-induced ETV4 downregulation is a reversible phenotype, we used a scratch assay, where rapid reactivation of ETV4 expression was observed in the cells adjacent to the scratches along with increased cell area (Fig. 1f). The cell-density-mediated regulation of ETV4 expression was also confirmed in other epithelial cell lines (Fig. 1g and Extended Data Fig. 1g). Direct manipulation of cell area by a cell stretching system revealed that a 15% decrease in cell area was sufficient to diminish the nuclear expression of ETV4 (Fig. 1h,i and Extended Data Fig. 1h–j).
Substrate stiffness is another physiologically relevant stimulus that regulates cellular mechanical stress28. Reduced mechanical stress by soft substrates was confirmed by immunostaining for pMLC and YAP/TAZ (Extended Data Fig. 1k,l). Compared with plastic (~1 Gpa), the nuclear expression of ETV4 proteins was significantly reduced when cells were placed on soft substrates (Fig. 1j,k). Replating hESCs to plastic resumed with high ETV4 expression (Fig. 1l). Overall, we have identified ETV4 as a TF whose expression is regulated by various mechanical cues.
Regulation of ETV4 expression by cell crowding dynamics
Like the in vivo epiblast, hESCs intrinsically grow, forming a single-layered epithelial colony in a culture dish (Fig. 2a). As hESC colonies expanded, we observed dynamic spatiotemporal changes in cell density. In small colonies (<0.5 mm2), individual cells exhibited a relatively large cell area, indicating a lesser degree of crowding (Fig. 2b). Conversely, in larger colonies (>2 mm2), cells located in the centre were smaller, while those at the periphery retained a larger cell area (Fig. 2b and Extended Data Fig. 2a). A gradual decrease in cell area was observed from the periphery to the centre of large hESC colonies (Fig. 2c). Because cell area exhibited a strong correlation with nucleus size (Extended Data Fig. 2b), we conducted real-time measurements of single nucleus sizes in an hESC line expressing a nuclear reporter, H2B-GFP. As the hESC colony expanded, an increase in single-cell variation in nucleus size was observed, with cells in the centre becoming smaller (Fig. 2d). These results collectively underscore the dynamic regulation of local cell crowding during epithelial expansion.
Given the mechanosensitive regulation of ETV4 expression, we investigated the ETV4 expression patterns in H9 and H1 hESC colonies. ETV4 was homogeneously expressed in the nucleus of cells from small hESC colonies (<0.5 mm2; Fig. 2e and Extended Data Fig. 2c). Strikingly, large hESC colonies (>2 mm2) showed zonation of ETV4 expression with a reduction in the crowded centre (Fig. 2e and Extended Data Fig. 2c). However, the expression of core pluripotency genes (OCT4, NANOG and SOX2) remained high across the whole colony (Fig. 2e and Extended Data Fig. 2c,d). These results were confirmed in hESCs cultured either in a different medium or on a different coating extracellular matrix (ECM; Extended Data Fig. 2e,f). Clonal colonies derived from single hESCs also showed the zonation of ETV4 expression (Extended Data Fig. 2g). Consistent with the zonation pattern, sharp reduction in ETV4 expression was observed with increased cell density (Fig. 2f). The zonation pattern of ETV4 predominantly emerged when the colony diameter exceeded 1,000 μm, and the size of the ETV4-low area increased proportionally to the whole colony size (Fig. 2g). Reversible ETV4 activation in the centre of large colonies was confirmed by a scratch assay (Fig. 2h).
Because cell density influences various cellular processes, we first measured the cell shape index (CSI) to evaluate cell circularity at the boundary where ETV4 expression sharply transitions. High ETV4-expressing cells exhibited significantly larger cell areas than low ETV4-expressing cells, while both types of cells had similar CSIs (Extended Data Fig. 3a). These results suggest that ETV4 expression is associated with cell area rather than cell geometry. Next, we studied the potential effect of cell density on cellular metabolism. GSEA, based on the RNA-seq data from crowded hESCs, revealed no significant difference in the expression of glycolysis- and hypoxia-related genes (Extended Data Fig. 3b). To directly measure hypoxia, we used the hypoxyprobe system, in which pimonidazole hydrochloride forms protein adducts in hypoxic cells (Extended Data Fig. 3c)29. Consistent with the GSEA results, there was no evident induction of hypoxia in the crowded centre of large hESC colonies (Extended Data Fig. 3c). Furthermore, glucose uptake measurements demonstrated similar uptake rates between the periphery and the centre of large hESC colonies (Extended Data Fig. 3d,e). Although we cannot completely rule out the potential effect of cellular metabolism, these data indicate that cell-density-mediated ETV4 expression primarily depends on cell area.
Interestingly, crowded cells in the centre of large hESC colonies maintained active nuclear expression of YAP proteins despite the sharp inactivation of ETV4 (Fig. 2i and Extended Data Fig. 3f,g). We employed two shRNAs to target YAP (Extended Data Fig. 3h) and subsequently validated the accuracy of the YAP immunostaining results (Extended Data Fig. 3i). To gain deeper insights, we seeded hESCs at varying densities. Although nuclear expression of both ETV4 and YAP proteins was seen in a low density (~2,000 cells mm–2), inactivation of ETV4, but not YAP, took place in a medium density (~5,000 cells mm–2; Fig. 2j). At a high density (~7,500 cells mm–2), nuclear expression of both ETV4 and YAP proteins was reduced (Fig. 2j). The density ranges we used in this study are equivalent to those from other studies based on 2D gastruloid models6. These findings confirm that ETV4 and YAP respond differently to changes in cell density in hESCs. Given the crucial role of YAP/TAZ in epithelial proliferation10, a uniformly high YAP/TAZ level (Fig. 2i and Extended Data Fig. 3f,g) conforms to active cell division and high levels of the proliferation marker Ki67 and 5-ethynyl-2′-deoxyuridine (EdU) incorporation in the crowded centre (Extended Data Fig. 3j–l). YAP inhibition significantly diminished hESC colony growth (Extended Data Fig. 3m). Overall, these results demonstrate spatiotemporal regulation of ETV4 expression by cell crowding dynamics in a growing hESC epithelium.
Derepression of the neuroectoderm fate by ETV4 inactivation
Pluripotent epiblasts undergo dynamic changes in differentiation potential as the epithelial disc undergoes size expansion before gastrulation2,30. Consistently, growing evidence shows that the size of in vitro hESC colonies influences the differentiation propensity7,31,32, suggesting a direct role of epithelial expansion in lineage specification. Indeed, hESCs showed a biased differentiation toward mesendoderm (ME) when they were differentiated in a small colony. Under a robust neuroectoderm (NE)-directed differentiation condition called dual SMAD inhibition33, small colonies of hESCs were unable to produce NE cells expressing PAX6, which is a necessary and sufficient NE marker gene in humans (Fig. 3a and Extended Data Fig. 4a)34. However, ME-directed differentiation by BMP4 and FGF2 efficiently turned small hESC colonies into Brachyury+ ME cells (Fig. 3a and Extended Data Fig. 4a)35,36. In stark contrast, large hESC colonies produced both NE and ME lineage cells with clear spatial separation. PAX6+ NE cells predominantly emerged in the centre under dual SMAD inhibition, whereas ME cells were derived in the periphery region in the presence of BMP4 and FGF2 (Fig. 3a and Extended Data Fig. 4a). These results follow previous findings based on micropatterning technology6,7. The NE differentiation potential of large colonies was confirmed when hESCs were spontaneously differentiated by FGF2 and TGF-β deprivation (Extended Data Fig. 4b).
Based on the above results, we hypothesized that the NE fate, initially repressed in small colonies, can be derepressed in the centre as cell crowding occurs during colony expansion. Indeed, high cell density dramatically promoted NE differentiation at the expense of ME derivation (Fig. 3b), consistent with the findings of a previous report33. Because ETV4 expression is downregulated in the crowded centre of large colonies, ETV4 could play a key role in suppressing the NE fate. The pre-patterns of ETV4 inactivation in undifferentiated hESC colonies precisely matched the size of the PAX6+ NE area in differentiated colonies (Fig. 3c). ETV4 expression was decreased in NE cells but not in ME cells, supporting the suppressive role of ETV4 in NE derivation (Extended Data Fig. 4c). Accordingly, ETV4 overexpression completely blocked the emergence of PAX6+ NE cells in large colonies (Fig. 3d and Extended Data Fig. 4d). Furthermore, ETV4 knockdown (KD) enabled small colonies to differentiate into NE cells (Fig. 3e and Extended Data Fig. 4e,f), which was blocked by ectopic ETV4 expression (Extended Data Fig. 4g). By contrast, ETV4 depletion impeded Brachyury+ ME cell differentiation (Extended Data Fig. 4h). Finally, cell crowding-mediated NE promotion was blocked by ETV4 overexpression (Extended Data Fig. 4i). These results demonstrate that ETV4 downregulation by cell crowding underlies NE lineage derepression.
To investigate the transcriptome-wide effect of ETV4 in hESCs, we performed RNA-seq after ETV4 KD (Fig. 3f and Supplementary Table 3). GO analysis of total DEGs (|log2(fold change)| > 0.5, adjP < 0.05) showed significant enrichment in embryonic development and morphogenesis in top 10 GO terms (Extended Data Fig. 5a), supporting the key role of ETV4 in early lineage determination. Moreover, DEGs upregulated by ETV4 KD included genes related to neural differentiation with nervous system development in TOP10 GO terms (Fig. 3f and Extended Data Fig. 5b), confirming the repressive role of ETV4 in the NE derivation. To pinpoint the molecular mechanisms of ETV4, we focused on downregulated DEGs because ETV4 primarily acts as a transcriptional activator37,38. Downregulated DEGs included many genes related to ECM remodelling, such as matrix metalloproteinases (MMPs; Fig. 3f), with ECM organization in top 10 GO terms (Extended Data Fig. 5c). Recently, it was reported that N-cadherin marks cells in the periphery of hESC colonies39. A re-analysis of published single-cell RNA-seq (scRNA-seq) data from sorted peripheral N-cadherin+ cells confirmed the elevated expression of a well-established ETV4 target gene, DUSP6 (Extended Data Fig. 5d,e)40. Furthermore, genes related to ECM remodelling, such as MMPs, were significantly upregulated in N-cadherin+ peripheral cells with ECM in top 10 GO terms (Extended Data Fig. 5e,f and Supplementary Table 4). These results support the role of ETV4 in ECM remodelling in the periphery region of hESC colonies.
ETV4 is a known direct upstream regulator of MMPs in cancer37,41,42,43. ETV4 KD or overexpression altered the expression of MMPs in hESCs (Extended Data Fig. 5g–i). The treatment of pan MMP inhibitors (GM6001 and BB94) was sufficient to derepress the NE fate in small hESC colonies (Fig. 3g). By contrast, MMP inhibition disrupted ME differentiation (Fig. 3g). Interestingly, membrane-type MMPs (MT-MMPs) such as MMP14 showed higher expression in hESCs than other MMPs (Extended Data Fig. 5j). The overexpression of MMP14 phenocopied ETV4 overexpression (Extended Data Fig. 5k,l) and blocked NE derepression by ETV4 KD (Fig. 3h). Overall, these results suggest that ETV4 is an NE repressor linking cell crowding dynamics to lineage specification.
Spatiotemporal ETV4 expression regulated by ERK
The MAPK signalling pathway is a well-known regulator of PEA3 family TFs in cancer44. To investigate this molecular link, we took advantage of the kinase translocation reporter (KTR) system45. When the kinase of interest is active, fluorescently-tagged substrates are phosphorylated and localized in the cytoplasm (Fig. 4a). The KTRs for ERK, p38 and JNK were validated by specific inhibitors (Extended Data Fig. 6a). For p38, individual cells displayed highly variable kinase activities with no discernible difference between the centre and periphery of large colonies (Extended Data Fig. 6b), whereas JNK activities were low in most cells (Extended Data Fig. 6c). However, the ERK-KTR showed clear cytoplasmic localization in most cells of small colonies (Fig. 4a). As a colony grew, a sharp reduction in ERK activity was observed in the crowded centre (Fig. 4b and Extended Data Fig. 6d,e). Induction of cell crowding was sufficient to inactivate ERK (Extended Data Fig. 6f). Live cell imaging of ERK activity revealed that the ERK activity pattern closely resembled the expression pattern of ETV4 (Fig. 4c).
ERK inhibition by chemical inhibitors (PD0325901 and U0126) induced a rapid decrease in ETV4 protein abundance without affecting the mRNA level (Fig. 4d and Extended Data Fig. 6g,h). Proteasome inhibition blocked the decrease, suggesting that ERK regulates ETV4 protein stability (Fig. 4d). ERK inhibition dramatically reduced the half-life of ETV4 proteins from 4.15 h to 0.35 h (Fig. 4e). COP1 is a critical E3 ligase for the degradation of PEA3 family TFs46,47. In large hESC colonies, COP1 was primarily localized in the nucleus with homogenous expression (Extended Data Fig. 7a). COP1 KD slightly increased the basal level of ETV4 protein (Extended Data Fig. 7b–d) and nullified the effect of ERK inhibition on ETV4 expression (Extended Data Fig. 7d). COP1 depletion also blocked ETV4 downregulation in the crowded centre of large hESC colonies and suppressed NE differentiation (Extended Data Fig. 7e,f).
Finally, ERK inhibition derepressed the NE fate and inhibited ME differentiation in small hESC colonies (Fig. 4f). NE derepression by ERK inhibition was completely blocked by ETV4 overexpression (Fig. 4g). ERK activation by constitutively active KRASG12V increased ETV4 protein levels and suppressed NE differentiation (Extended Data Fig. 7g–i). While ERK signalling is known to have a broad impact on numerous regulatory factors, our findings suggest ETV4 as a primary target of ERK in the context of lineage specification.
Cell-crowding-mediated regulation of receptor endocytosis
The FGF and TGF-β signalling pathways play crucial roles in maintaining pluripotency48,49,50,51. Short-term treatment of A83-01 (TGF-β inhibitor) showed no effect on ERK activity (Fig. 5a); however, SU-5402 (FGF inhibitor) reduced ERK activity and the level of ETV4 proteins without altering mRNA expression (Fig. 5b,c). We used shRNAs targeting FGFR1 owing to the high expression in hESCs (Extended Data Fig. 8a–c). FGFR1 KD decreased ERK activity and ETV4 expression (Extended Data Fig. 8d,e), suggesting FGFR1-mediated signalling as an upstream pathway of ERK and ETV4.
An emerging body of evidence shows that receptor localization contributes to regulating downstream signalling52,53. For example, it was recently reported that lateral localization of TGF-β receptors impeded cellular responses to apically applied ligands6,54. However, we observed that cell-crowding-induced ETV4 inactivation occurred in hESC culture on a transwell system, where ligands are accessible to both apical and basolateral sides of cells (Extended Data Fig. 8f). In the case of receptor tyrosine kinases, endocytosis is functionally related to the downstream signalling activation52. Activated EGFR and FGFR accumulate in endosomes that serve as a signalling platform55, and blocking receptor endocytosis inhibits the downstream signalling56,57,58. Recently, it was reported that membrane tension influenced FGFR endocytosis and downstream ERK activity in mouse ESCs, implicating mechanical regulation of receptor endocytosis59,60. To test this idea, we generated stable hESC lines expressing GFP reporters for various endosomes and lysosomes61. In small hESC colonies, a significant number of FGFR1-containing vesicles colocalized with RAB5A+ early endosomes (Fig. 5d), a finding that was confirmed using another early endosome marker, EEA1 (Fig. 5e). FGFR1-containing vesicles were also observed in late endosomes (RAB7A), recycling endosomes (RAB11B) and lysosomes (LAMP1; Fig. 5d). These results suggest that FGFR1 proteins are under active endocytosis, followed by recycling or degradation. Strikingly, cells in large hESC colonies showed distinct subcellular localizations of FGFR1 proteins depending on their positions. These proteins were predominantly present in endosome vesicles in the periphery of large colonies (Fig. 5e and Extended Data Fig. 9a). In stark contrast, membrane localization of FGFR1 proteins was observed in the crowded centre, suggesting impaired endocytosis (Fig. 5e and Extended Data Fig. 9a, b). FGFR1 proteins were localized in both apical and basolateral membranes in the crowded centre (Extended Data Fig. 9a). The transition of FGFR1 localization occurred abruptly along the radial axis of large hESC colonies (Fig. 5e). To further validate these results, we applied a DNA aptamer (TD0) previously developed as an agonist of FGFR1 (Extended Data Fig. 9c)62. Fluorophore-tagged aptamers were found in FGFR1-containing early endosomes in the periphery of large hESC colonies (Extended Data Fig. 9d). By contrast, few aptamer-containing early endosomes were detected in the centre (Extended Data Fig. 9d). Furthermore, induction of cell crowding was sufficient to impair FGFR1 endocytosis with loss of EEA1 colocalization (Fig. 5f). These results show that FGFR1 endocytosis is tightly regulated by cell crowding.
Next, we used Pitstop2 (a clathrin inhibitor) and Dynasore (a dynamin inhibitor) to suppress endocytosis (Extended Data Fig. 9e)63,64. Endocytosis inhibition dramatically decreased the number of FGFR1-containing endosomes, ERK activity and ETV4 expression (Fig. 5g–i and Extended Data Fig. 9f–h), which phenocopied the effects of FGFR inhibition (Fig. 5b,c). Constitutively active KRASG12V blocked the effect of endocytosis inhibition on ETV4 expression (Extended Data Fig. 9i), suggesting that endocytosis regulates ETV4 through the FGFR–RAS signalling.
Consistent with previous research65,66, we found that FGFR inhibition was sufficient to derepress the NE fate within small hESC colonies (Fig. 5j and Extended Data Fig. 9j), while suppressing the ME derivation (Fig. 5j). Notably, we observed a similar phenotypic outcome upon the inhibition of endocytosis (Fig. 5j). Furthermore, ectopic ETV4 expression entirely counteracted the NE-promoting effects induced by FGFR and endocytosis inhibitors (Fig. 5k,l). Overall, our findings suggest that cell-crowding-mediated blockade of FGFR1 endocytosis leads to the inactivation of both ERK and ETV4.
Regulation of receptor endocytosis by integrin–actomyosin
To figure out molecular mechanisms bridging cell crowding to FGFR1 endocytosis, we revisited the RNA-seq data performed in crowded hESCs. Downregulated DEGs in crowded hESCs showed significant enrichment of genes related to cell–ECM interactions such as integrins (Fig. 6a and Extended Data Figs. 1b and 10a). Published scRNA-seq data from N-cadherin+ peripheral cells also showed focal adhesion in top 10 KEGG pathways (Extended Data Fig. 10b)39. These results concord with the fact that cell crowding leads to a smaller surface area through which cells interact with the ECM.
Consistent with the transcriptomic results, a high level of integrin β1 was observed in paxillin+ focal adhesions in small colonies and the periphery of large colonies (Fig. 6b). By contrast, crowded cells in the centre showed perinuclear localization of integrin β1 with loss of paxillin colocalization (Fig. 6b). Focal adhesion kinase (FAK) is a critical component of integrin signalling67. Immunostaining of phosphorylated FAK (pFAK, active FAK) revealed the sharp downregulation in response to increased cell density in hESCs (Fig. 6c). Furthermore, FAK inhibition suppressed FGFR1 endocytosis, ERK activity and ETV4 expression (Extended Data Fig. 10c–f). These findings suggest that integrin signalling bridges cell crowding to FGFR1 endocytosis.
Although integrin signalling has the potential to impact a wide array of cellular processes, we focused on its role in regulating cytoplasmic actin filaments, given the significance of actin dynamics in receptor endocytosis68,69,70. Consistent with the patterns of integrin β1+ focal adhesions, small hESC colonies showed well-established pMLC+ bundles, while the zonation of pMLC staining appeared in large hESC colonies with a sharp reduction in the centre (Fig. 6d,e). Disassembly of pMLC+ bundles was observed by integrin inhibition (RGDS peptide), FAK inhibition or cell crowding (Extended Data Figs. 1e and 10g,h). We used either blebbistatin (myosin inhibitor) or YM (Y-27632, RhoA kinase inhibitor; ML-7, myosin light chain kinase inhibitor) to inhibit actomyosin (Extended Data Fig. 10i,j)71. Transient inhibition of actomyosin activity decreased the number of FGFR1-containing endosome vesicles and inactivated ERK and ETV4 (Fig. 6f–h and Extended Data Fig. 10d,k,l). Importantly, the transient activation of actomyosin by doxycycline-inducible expression of constitutively active RhoAQ63L partially re-activated FGFR1 endocytosis in the crowded centre of large hESC colonies (Fig. 6i and Extended Data Fig. 10m). These findings suggest that actomyosin, in addition to serving as a marker of cellular mechanical stress, plays a crucial role in regulating receptor endocytosis in response to mechanical stimuli. However, it is essential to acknowledge that other cellular changes induced by cell crowding, such as the reduction in apical surface area for ligand interaction and alterations in cell geometry, may also contribute to regulating receptor endocytosis.
Finally, actomyosin inhibition derepressed the NE fate in small hESC colonies while preventing ME differentiation (Fig. 6j). The NE derepression by actomyosin inhibition was blocked by ETV4 overexpression (Fig. 6k). Overall, these results suggest that the integrin–actomyosin pathway serves as a crucial link connecting cell crowding to FGFR1 endocytosis, ETV4 expression and NE specification (Fig. 6l).
Dynamic ETV4 expression captured by mathematical modelling
A critical question that remains to be answered in this study is how the gradient of cell crowding from the periphery to the centre produces a sharp boundary in ETV4 expression. To address this question, we propose a mathematical framework that describes a cell crowding model. Based on real-time imaging data of H2B-GFP hESC colony growth, we obtained an analytical solution to the governing reaction-diffusion partial differential equation describing cell crowding dynamics (Methods). This equation assumes a cell population density varying in time and space due to the flow of cells (diffusion) from high- to low-density regions (outward flow) as well as cell division (reaction). A similar approach has been used to describe bacterial colony growth72. Our model accurately predicted the colony population and colony size (diameter; Fig. 7a). Moreover, the model predicted an increase in cell density over time and higher cell crowding at the colony centre (Fig. 7b).
Cell crowding inversely correlates with ECM-accessible area, linking ECM ligand availability to cell density. Moreover, it is well established that integrin–ECM binding kinetics follow a cooperative interaction73,74. Such an interaction enforces a system with an on/off switch-like behaviour. We used the Hill equation to model integrin–ECM interactions75. The equation captures the biomolecular interaction that exhibits the cooperativity between two binding molecules. Our model predicts a switch for integrin activity that occurs at different distances from the colony’s centre as the colony grows over time (Fig. 7c). The ultrasensitive transition of integrin activity appears when the core cell density reaches a critical density of 5,000 cells mm–2 (Fig. 7d), which matches well with the pFAK and ETV4 experimental data (Figs. 2f, 6c and 7e). Furthermore, our model recapitulated the experimental data with the ETV4 low-expression region appearing in the centre when the colony size reaches 1,000–1,500 μm in diameter (Fig. 7f,g). The critical cell density can only be observed after t > tcr (where t represents time and tcr represents the critical time when the ETV4 transition boundary appears for the first time), and the radius at which the ETV4 transition occurs varies as a function of time (Fig. 7h). Although our model successfully demonstrated the ultrasensitive transition of integrin activity and ETV4 expression based on the ECM availability and the cooperativity in integrin–ECM interaction, it is essential to acknowledge that cell shape itself may also contribute to regulating integrin signalling or the downstream pathways independently of ECM availability76. Overall, our model bridges the simple growth of stem cell epithelium to the timing and location of ETV4 inactivation and NE derepression.
Discussion
An outstanding question in developmental biology is how multiple lineages arise from a single-layered epithelium of seemingly homogeneous stem cells. During the expansion of hESCs, active proliferation prompts local cell crowding in the centre of epithelial colonies. A gradual change in cell crowding induces the abrupt inactivation of integrin signalling based on the cooperativity of the integrin–ECM interaction. The blunted integrin–actomyosin pathway impedes FGF receptor endocytosis, thereby derepressing the NE fate from ETV4 (Fig. 6l). Importantly, ETV4, formerly known as an oncogene, serves as a mechanotransducer linking cell crowding dynamics to lineage specification in a stem cell epithelium.
Previous studies have employed hESC-derived gastruloid models grown on micropatterned plates to study molecular mechanisms for embryonic germ layer derivation6,7,77,78. Although a differentiation signal, BMP4, affects all cells in a culture plate, distinct lineages emerge in a defined position along the radial axis of hESC colonies. Reaction–diffusion models have been used to explain lineage patterning, where asymmetric gradients of signal activators and inhibitors determine the positional identity of differentiating cells79,80. In hESCs, BMP4 directly stimulates the expression of NOGGIN, a BMP antagonist, providing a reaction–diffusion mechanism for pattern formation6,7. In addition, cell-density-dependent relocalization of TGF-β receptors to the basolateral side restricts signalling activation, contributing to gastruloid formation6. In this study, we identified ETV4 as a mechanotransducer for lineage determination. ETV4 responds to cell crowding dynamics evoked by epithelial expansion, generating a pre-pattern with a sharp boundary for future lineage fates. Our findings offer an independent patterning mechanism, directly connecting mechanical microenvironments to lineage specification.
In an in vitro ESC differentiation, NE derivation is deemed the default pathway because minimal medium conditions without additional signalling molecules are sufficient to direct ESCs toward the NE lineage81. However, we found that derepression from ETV4 is essential for successful NE derivation. This finding suggests that, at first, hESCs are not inclined towards NE differentiation. The withdrawal of signalling molecules such as FGF acts as a signal to inactivate ETV4 and thereby induces derepression of the NE fate in hESCs. These results provide a new perspective in understanding lineage fate determination during early hESC differentiation.
Core pluripotency genes are essential for the self-renewal of hESCs. At the same time, these genes play critical roles in regulating lineage specification82. Upon differentiation, OCT4 and NANOG promote the ME lineage, while SOX2 drives NE derivation83,84. These findings suggest that the core pluripotency genes act as the uppermost regulators of stem cell self-renewal and differentiation. In this study, however, we propose that ETV4 impacts lineage fate determination earlier than core pluripotency genes. In large hESC colonies, ETV4 protein expression sharply decreases in the crowded centre before the expression of any core pluripotency genes alters, highlighting ETV4 as the earliest determinant of lineage fates in hESCs.
Epithelial cell crowding activates contact inhibition, mediated by YAP/TAZ and PIEZO1, to suppress proliferation and maintain tissue homeostasis10,13,15,85,86. However, the mechanism by which cell crowding decouples lineage fate determination from contact inhibition of proliferation is unclear. Here, we identified that ETV4 acts as an on/off switch with ultrasensitive dependence on cell crowding. Sharp ETV4 inactivation by cell crowding derepresses the NE fate. At the same time, sustained YAP activity safeguards stem cell proliferation in a high-density environment. In such a setting, hESCs can utilize spatiotemporal heterogeneity in cell density to derive multiple lineages while maintaining active proliferation. Beyond stem cell differentiation, ETV4 is a critical oncogene with elevated expression in multiple cancers44,87. Therefore, our discovery of ETV4 as a mechanical transducer provides new insights into the mechanisms underlying tumour progression and suppression.
Methods
Maintenance of cell lines
H9 (WA09) and H1 (WA01) hESCs were purchased from WiCell and cultured on Matrigel (Corning) either in mTeSR1 (Stem Cell Technologies) or TeSR-E8 (Stem Cell Technologies) medium. H9 and H1 hESCs were authenticated by short tandem repeat analysis. hESCs were maintained at 37 °C with 5% CO2 and were passaged every 3–4 days by ReLeSR (Stem Cell Technologies). This work was approved by the Human Stem Cell Research Oversight Committee at Pohang University of Science and Technology (PIRB-2021-R035). HEK293T cells (CRL-11268, ATCC) and MCF-7 cells (HTB-22, ATCC) were cultured in DMEM supplemented with 10% fetal bovine serum (FBS). ARPE19 cells (CRL-2302, ATCC) were cultured in DMEM/F12 with 10% FBS. Cells were maintained at 37 °C with 5% CO2 and were passaged every 2–3 days by trypsin (Welgene). To manipulate substrate stiffness, H9 hESCs and MCF7 cells were plated and cultured on 35 mm imaging dishes with 15 and 1.5 kPa PDMS layers (ibidi). To facilitate cell attachment, we applied the dishes to a thin Matrigel coating (1% Matrigel diluted in basal media). The stiffness values reported in our study represent the properties of the underlying plastic and hydrogel substrates.
Differentiation of hESCs
For spontaneous differentiation, hESCs on Matrigel were induced to differentiate with basal hESC culture media (DMEM/F12, 15% knockout serum replacement, MEM nonessential amino acid solution, and 0.1 mM 2-mercaptoethanol) without FGF2 and TGF-β. For neuroectoderm differentiation, hESCs were differentiated on Matrigel with hESC culture media (DMEM/F12, 5% knockout serum replacement, MEM nonessential amino acid solution, and 0.1 mM 2-mercaptoethanol) containing SB431542 (PeproTech, 10 μM) and Dorsomorphin (Tocris, 2 μM). For mesendoderm differentiation, hESCs were differentiated in mTeSR1 medium (containing FGF2) supplemented with recombinant human BMP4 (R&D Systems, 5 ng ml–1).
High-density culture of hESCs
Small drops of Matrigel (2 μl) were deposited on a chamber slide glass (Thermo) to create a controlled pattern of the adhesive surface. The prepared micropattern chamber slides were sealed with parafilm and incubated for 1 h at room temperature before use. hESCs were dissociated by ReLeSR (Stem Cell Technologies) and resuspended in 1 ml of mTeSR1 or TeSR-E8 medium. The entire resuspended solution was seeded on a micropattern chamber slide. After cell adhesion, excess cells were removed by washing, and adherent cells were left to grow for ~2–3 days to cover the whole patterned area. The cell densities used in this study, ranging from 2,000 to 10,000 cells mm–2, match the parameters established in prior research utilizing 2D gastruloid models6.
Assessment of hypoxic cells
Hypoxic cells were detected using a Hypoxyprobe kit (Hypoxyprobe) according to the manufacturer’s instructions. hESCs cultured under a hypoxic condition (5% O2, 5% CO2 and 90% N2) were used for validation. Cultured cells were treated with pimonidazole to a final concentration of 400 µM at 37 °C for 2 h. Hypoxic cells were then detected by fluorescence imaging.
Glucose uptake assay
Glucose uptake was assessed using a glucose uptake assay kit (Dojindo) following the manufacturer’s instructions. To validate the assay, hESCs were incubated at 4 °C for 15 min to inhibit glucose uptake. Cultured cells were treated with the Glucose Uptake Probe to a final concentration of 5 μM at 37 °C for 15 min. The cells were washed three times with cold washing solution and replenished with cold washing solution at 4 °C. Glucose uptake was subsequently visualized by fluorescence imaging.
Lentiviral preparation and concentration
Lentiviral vector plasmids were transfected into HEK293T cells with second-generation packaging vectors psPAX2 (Addgene, 12260) and pMD2.G (Addgene, 12259; Supplementary Table 6). Transfected HEK293T cells were cultured at 37 °C with 5% CO2 for 2 days. Supernatants were filtered using 0.45 μm filters (Corning) and concentrated overnight at 4 °C using Lenti-X Concentrator (Takara Bio). Concentrated lentiviruses were resuspended in PBS and stored at −80 °C.
Generation of hESC reporter lines
Concentrated lentiviral particles for KTRs and endosome/lysosome reporters were transduced to H9 hESCs grown on Matrigel with TeSR-E8. Two days after transduction, 2 μg ml–1 puromycin (Gibco) was added to culture medium. Puromycin selection was performed for at least 4 days to generate stable hESC reporter lines.
Immunofluorescence staining
All samples were fixed using 3.7% methanol-free formaldehyde (Thermo) for 15 min at room temperature. After fixation, samples were washed 3 times with PBS and permeabilized with 0.25% Triton X-100 (SIGMA)-supplemented PBS for 10 min at room temperature. Samples were blocked with 10% FBS in PBS for 1 h at room temperature. The primary antibodies diluted in the blocking solution were treated overnight at 4 °C (Supplementary Table 5). After washing three times with PBS, an appropriate Alexa Fluor dye-conjugated secondary antibody was used to treat samples for 1 h at room temperature. Nuclear staining was performed on samples for 2 min with Hoechst 33342 (Thermo). Images were captured by a fluorescence microscope (Leica DMi8) or confocal microscope (ZEISS LSM800). Sample images were prepared in imageJ 1.53 software (Fiji), and statistical analysis was performed using the GraphPad Prism 9.1.0 software (GraphPad Software).
Cell shape analysis
The CSI is calculated based on measurements obtained from ZO1 staining images, where each cell’s boundary is outlined, and its area and perimeter are measured. To calculate the CSI, we used the following formula: CSI = 4π × area/perimeter².
EdU incorporation assay
To label actively proliferating cells, we used the Click-iT EdU Imaging Kit (Thermo) following the manufacturer’s instructions. hESCs were treated with EdU to a final concentration of 5 µM at 37°C for 6 h. After EdU labelling, hESCs were washed with PBS and fixed using 3.7% methanol-free formaldehyde (Thermo) for 15 min at room temperature. EdU was detected by the Click-iT reaction protocol.
Crystal violet staining
Cells were washed once with PBS and stained at room temperature for 1 min with 1 ml of crystal violet staining solution (SIGMA). After staining, samples were washed three times with PBS, and dried at room temperature for 15 min.
Quantitative real-time PCR
Total RNA was isolated using QIAzol Lysis Reagent (QIAGEN) and reversely transcribed with SuperiorScript III Master Mix (Enzynomics). Quantitative RT-PCR analysis was performed on the CFX Connect Real-Time PCR Detection System (BIO-RAD) with TOPreal qPCR 2X PreMIX (Enzynomics) (Supplementary Table 7). GAPDH was used as a normalization control. Results were plotted using the GraphPad Prism 9.1.0 software (GraphPad Software).
Western blot
Cells were lysed with RIPA lysis buffer supplemented with Protease Inhibitor Cocktail Kit5 (Quartett). Protein concentration was measured using the BSA Protein Assay kit (Thermo). The same amount of protein sample was separated on SDS-PAGE and then transferred to a nitrocellulose membrane (BIO-RAD). Membranes were blocked with PBST (0.1% Tween 20 in PBS) containing 5% skim milk; then immunoblotting was performed overnight at 4 °C with primary antibodies (Supplementary Table 5). The membranes were stained with an appropriate HRP-conjugated secondary antibody for 1 h at room temperature. Signals were detected with Amersham imager 680 (Amersham), and relative signal intensity was quantified by Multi gauge 3.0 (FUJIFILM).
Cell stretching experiment
A manual cell stretching system STB-100 (Togetherbio) was used to control cell area directly. Stretch chambers (STB-CH-4W) were coated with Matrigel. Cells were seeded in a chamber that was stretched by 20% for 2 days. The chamber was then gently released to reduce cell area and incubated for 24 h before further experiments.
Cell scratch experiment
Cells were scraped using P10 pipet tips, washed with 1 ml of culture medium to get rid of the debris, and then fed with 2 ml of culture medium.
RNA sequencing analysis
Using the Illumina TruSeq Stranded mRNA Sample Preparation Kit (Illumina), 1 μg of the total RNA of the sample was prepared to create a library. The poly-A-containing mRNA molecules were purified using poly‐T‐attached magnetic beads. The mRNAs were copied into the cDNA library using SuperScript II reverse transcriptase (Invitrogen). Following the qPCR Quantification Protocol Guide (KAPA BIOSYSTEMS), the library was quantified using the KAPA Library Quantification Kit for Illumina Sequencing Platform and qualified using TapeStation D1000 ScreenTape (Agilent Technologies). Then, the indexed library was submitted to Illumina NovaSeq (Illumina), and pair-end (2 × 100 bp) sequencing was performed. The raw reads from the Illumina NovaSeq were trimmed using TRIMMOMATIC and mapped to the Homo sapiens transcript reference (GRCh37) using HISAT v2.1.0. Transcript assembly of known transcripts was processed by StringTie v2.1.3b. The expression abundance of genes was calculated as read counts or FPKM values per sample and filtered through statistical hypothesis testing for additional analysis such as DEGs.
Differential gene expression analysis
We used DESeq288 (https://bioconductor.org/packages/DESeq2/) to perform differential gene expression analysis using bulk RNA-seq data. To identify DEGs from scRNA-seq data, we used DEsingle89 implementation in R/Bioconductor (https://bioconductor.org/packages/DEsingle). P values were adjusted for multiple comparisons using the FDR approach.
Cellular crowding model
A reaction–diffusion equation is used to model colony growth dynamics as
by capturing spatial variation in colony cell population density through diffusion as well as cell proliferation (first and second term in equation (1)) where C defines cell population density and r is defined as the radial distance from the centre of the colony. D and α represent diffusion coefficient and cellular division rate and can be measured experimentally. For example, hESC division rate was inferred from the range of experimentally observed cellular doubling rate using α = ln(2)/(doubling rate) from previous reports90. Diffusion coefficient was also estimated using a least squares fit to Fig. 7a. Spatially homogenous solution to equation (1) represents colony population exponential growth (that is the number of cells in the colony) with the division rate α. With the following initial and boundary conditions
solution to equation (1) can be found as
where\({A}_{k}=\frac{2{C}_{0}{r}_{0}{J}_{1}({\lambda }_{k}{r}_{0})}{\,{{{\lambda }_{k}J}_{1}}^{2}({\lambda }_{k})},{J}_{0}\left({\lambda }_{k}\right)=0,k{\mathbb{\in }}{\mathbb{N}}\), t* and r* are scaled time and radius.
With Jm(r) as the Bessel function of the first kind and λk zeros of J0(r), equation (4) predicts cellular crowding at any location in time. Since equation (4) predicts a smooth radial distribution of cell population density, we defined the colony radius as the radius where cell population density reaches a constant as estimated experimentally based on Fig. 7b. To simplify the governing equation describing cell colony growth, we neglected active cell–cell interaction and cell surface adhesion, which could introduce several additional parameters to our model. Moreover, the diffusion equation assumes a continuous C(r,t), an assumption that is no longer valid around the boundary where cell population distribution becomes discrete and, therefore, we applied the boundary condition in equation (3) to ensure that below a threshold, cell population density is set to 0. Both assumptions have been used in studies modelling colonial bacterial growth72.
Based on cooperativity of an integrin–ECM interaction, we used Hill equation to model integrin–ligand interaction as
where θ, L, Kd and n represent integrin activity, ligand concentration, integrin–ligand dissociation constant, and the Hill coefficient, respectively. L is limited by cellular crowding such that in a high population density, ligand availability reduces due to the smaller contact between cell surface and ECM. Therefore, L is inversely proportional to cellular crowding and can be expressed as
where L0 is a constant. Using equations (5) and (6), we derived integrin activity as a function of cellular crowding as
where β is a constant and can be found experimentally. Furthermore, integrin activates FAK through tyrosine-phosphorylation (pFAK), which in turn upregulates ETV4 expression or in other words:
The Hill coefficient determines the degree of cooperativity and determines ETV4 transition width (that is, the sharpness of the sigmoid shape of θ) such that n → ∞ converges to a step function. However, the location of the transition (that is, ETV4/integrin transition diameter) is independent of n. At equilibrium, pFAK and ETV4 expression levels are directly correlated with θ with reaction equilibrium constant as the correlation coefficient. Normalization of pFAK and ETV4 expression was done using min–max normalization (min = 0, max = 1.2) to scale pFAK/ETV4 expression to be compared to model prediction in Fig. 7d. Together, this mathematical model can accurately predict both cell population dynamic and the spatiotemporal dynamic of ETV4 expression. Model parameters are measured experimentally and presented in Supplementary Table 8. β is obtained by least square fitting equation (7) into experimental measurement of cell population density and ETV4 expression profile. The estimate and the permutation error for β are reported in Supplementary Table 8.
Statistics and reproducibility
Statistical analyses were performed using GraphPad Prism 9.1.0 software (GraphPad Software) and ImageJ 1.53 software (Fiji). No statistical method was used to pre-determine sample size, but our sample sizes are similar to those reported in previous publications91,92. No data were excluded from the analyses. The experiments were not randomized, and the investigators were not blinded to allocation during experiments and outcome assessment. Data distribution was assumed to be normal, but this was not formally tested. All differences were compared using a two-tailed Student’s t-test. Error bars represent mean ± s.d. For image quantifications, precise numbers of quantified cells or colonies are provided in the figures or legends. Exact P values are available in Supplementary Table 9.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
The RNA-seq data generated from this study have been deposited in the Gene Expression Omnibus (GSE183702). Published bulk and single cell RNA-seq that were re-analysed here are available from the Gene Expression Omnibus (GSE69982 and GSE126022). Source data are provided with this study. All other data supporting the findings of this study are available from the corresponding author on reasonable request.
References
Sheng, G. Epiblast morphogenesis before gastrulation. Dev. Biol. 401, 17–24 (2015).
Smith, A. Formative pluripotency: the executive phase in a developmental continuum. Development 144, 365–373 (2017).
Nakamura, T. et al. A developmental coordinate of pluripotency among mice, monkeys and humans. Nature 537, 57–62 (2016).
Sasai, Y. Cytosystems dynamics in self-organization of tissue architecture. Nature 493, 318–326 (2013).
Zhu, M. & Zernicka-Goetz, M. Principles of self-organization of the mammalian embryo. Cell 183, 1467–1478 (2020).
Etoc, F. et al. A balance between secreted inhibitors and edge sensing controls gastruloid self-organization. Dev. Cell 39, 302–315 (2016).
Warmflash, A., Sorre, B., Etoc, F., Siggia, E. D. & Brivanlou, A. H. A method to recapitulate early embryonic spatial patterning in human embryonic stem cells. Nat. Methods 11, 847–854 (2014).
Gilmour, D., Rembold, M. & Leptin, M. From morphogen to morphogenesis and back. Nature 541, 311–320 (2017).
Eisenhoffer, G. T. & Rosenblatt, J. Bringing balance by force: live cell extrusion controls epithelial cell numbers. Trends Cell Biol. 23, 185–192 (2013).
Aragona, M. et al. A mechanical checkpoint controls multicellular growth through YAP/TAZ regulation by actin-processing factors. Cell 154, 1047–1059 (2013).
Dupont, S. et al. Role of YAP/TAZ in mechanotransduction. Nature 474, 179–183 (2011).
Eisenhoffer, G. T. et al. Crowding induces live cell extrusion to maintain homeostatic cell numbers in epithelia. Nature 484, 546–549 (2012).
Gudipaty, S. A. et al. Mechanical stretch triggers rapid epithelial cell division through Piezo1. Nature 543, 118–121 (2017).
Wada, K.-I., Itoga, K., Okano, T., Yonemura, S. & Sasaki, H. Hippo pathway regulation by cell morphology and stress fibers. Development 138, 3907–3914 (2011).
Zhao, B. et al. Inactivation of YAP oncoprotein by the Hippo pathway is involved in cell contact inhibition and tissue growth control. Genes Dev. 21, 2747–2761 (2007).
Miroshnikova, Y. A. et al. Adhesion forces and cortical tension couple cell proliferation and differentiation to drive epidermal stratification. Nat. Cell Biol. 20, 69–80 (2018).
Priya, R. et al. Tension heterogeneity directs form and fate to pattern the myocardial wall. Nature 588, 130–134 (2020).
Hashimoto, M. & Sasaki, H. Epiblast formation by TEAD-YAP-dependent expression of pluripotency factors and competitive elimination of unspecified cells. Dev. Cell 50, 139–154.e135 (2019).
Lian, I. et al. The role of YAP transcription coactivator in regulating stem cell self-renewal and differentiation. Genes Dev. 24, 1106–1118 (2010).
Nishioka, N. et al. The Hippo signaling pathway components Lats and Yap pattern Tead4 activity to distinguish mouse trophectoderm from inner cell mass. Dev. Cell 16, 398–410 (2009).
Panciera, T. et al. Induction of expandable tissue-specific stem/progenitor cells through transient expression of YAP/TAZ. Cell Stem Cell 19, 725–737 (2016).
Pathak, M. M. et al. Stretch-activated ion channel Piezo1 directs lineage choice in human neural stem cells. Proc. Natl Acad. Sci. USA 111, 16148–16153 (2014).
Ranade, S. S. et al. Piezo1, a mechanically activated ion channel, is required for vascular development in mice. Proc. Natl Acad. Sci. USA 111, 10347–10352 (2014).
Raudvere, U. et al. g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 47, W191–W198 (2019).
Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl Acad. Sci. USA 102, 15545–15550 (2005).
Qin, Q. et al. Lisa: inferring transcriptional regulators through integrative modeling of public chromatin accessibility and ChIP-seq data. Genome Biol. 21, 32 (2020).
Sharrocks, A. D. The ETS-domain transcription factor family. Nat. Rev. Mol. Cell Biol. 2, 827–837 (2001).
Wells, R. G. The role of matrix stiffness in regulating cell behavior. Hepatology 47, 1394–1400 (2008).
Varia, M. A. et al. Pimonidazole: a novel hypoxia marker for complementary study of tumor hypoxia and cell proliferation in cervical carcinoma. Gynecol. Oncol. 71, 270–277 (1998).
Pera, M. F. & Rossant, J. The exploration of pluripotency space: charting cell state transitions in peri-implantation development. Cell Stem Cell 28, 1896–1906 (2021).
Bauwens, C. L. et al. Control of human embryonic stem cell colony and aggregate size heterogeneity influences differentiation trajectories. Stem Cells 26, 2300–2310 (2008).
Peerani, R. et al. Niche-mediated control of human embryonic stem cell self-renewal and differentiation. EMBO J. 26, 4744–4755 (2007).
Chambers, S. M. et al. Highly efficient neural conversion of human ES and iPS cells by dual inhibition of SMAD signaling. Nat. Biotechnol. 27, 275–280 (2009).
Zhang, X. et al. Pax6 is a human neuroectoderm cell fate determinant. Cell Stem Cell 7, 90–100 (2010).
Wang, L. & Chen, Y.-G. Signaling control of differentiation of embryonic stem cells toward mesendoderm. J. Mol. Biol. 428, 1409–1422 (2016).
Yu, P., Pan, G., Yu, J. & Thomson, J. A. FGF2 sustains NANOG and switches the outcome of BMP4-induced human embryonic stem cell differentiation. Cell Stem Cell 8, 326–334 (2011).
De Launoit, Y. et al. The Ets transcription factors of the PEA3 group: transcriptional regulators in metastasis. Biochim. Biophys. Acta Rev. Cancer 1766, 79–87 (2006).
Xin, J.-H., Cowie, A., Lachance, P. & Hassell, J. A. Molecular cloning and characterization of PEA3, a new member of the Ets oncogene family that is differentially expressed in mouse embryonic cells. Genes Dev. 6, 481–496 (1992).
Nakanishi, M. et al. Human pluripotency is initiated and preserved by a unique subset of founder cells. Cell 177, 910–924.e922 (2019).
Znosko, W. A. et al. Overlapping functions of Pea3 ETS transcription factors in FGF signaling during zebrafish development. Dev. Biol. 342, 11–25 (2010).
Izumiyama, Y. et al. E1AF expression is closely correlated with malignant phenotype of tongue squamous cell carcinoma through activation of MT1-MMP gene promoters. Oncol. Rep. 13, 715–720 (2005).
Lynch, C. C., Crawford, H. C., Matrisian, L. M. & McDonnell, S. Epidermal growth factor upregulates matrix metalloproteinase-7 expression through activation of PEA3 transcription factors. Int. J. Oncol. 24, 1565–1572 (2004).
Vincenti, M. P., Schroen, D. J., Coon, C. I. & Brinckerhoff, C. E. v-src activation of the collagenase-1 (matrix metalloproteinase-1) promoter through PEA3 and STAT: requirement of extracellular signal-regulated kinases and inhibition by retinoic acid receptors. Mol. Carcinog. 21, 194–204 (1998).
Oh, S., Shin, S. & Janknecht, R. ETV1, 4 and 5: an oncogenic subfamily of ETS transcription factors. Biochim. Biophys. Acta Rev. Cancer 1826, 1–12 (2012).
Regot, S., Hughey, J. J., Bajar, B. T., Carrasco, S. & Covert, M. W. High-sensitivity measurements of multiple kinase activities in live single cells. Cell 157, 1724–1734 (2014).
Baert, J. et al. The E3 ubiquitin ligase complex component COP1 regulates PEA3 group member stability and transcriptional activity. Oncogene 29, 1810–1820 (2010).
Vitari, A. C. et al. COP1 is a tumour suppressor that causes degradation of ETS transcription factors. Nature 474, 403–406 (2011).
James, D., Levine, A. J., Besser, D. & Hemmati-Brivanlou, A. TGFβ/activin/nodal signaling is necessary for the maintenance of pluripotency in human embryonic stem cells. Development 132, 4495–4509 (2005).
Vallier, L., Alexander, M. & Pedersen, R. A. Activin/Nodal and FGF pathways cooperate to maintain pluripotency of human embryonic stem cells. J. Cell Sci. 118, 4495–4509 (2005).
Xu, C. et al. Basic fibroblast growth factor supports undifferentiated human embryonic stem cell growth without conditioned medium. Stem Cells 23, 315–323 (2005).
Xu, R.-H. et al. Basic FGF and suppression of BMP signaling sustain undifferentiated proliferation of human ES cells. Nat. Methods 2, 185–190 (2005).
Bergeron, J. J., Di Guglielmo, G. M., Dahan, S., Dominguez, M. & Posner, B. I. Spatial and temporal regulation of receptor tyrosine kinase activation and intracellular signal transduction. Annu. Rev. Biochem. 85, 573–597 (2016).
Kim, S. K. Polarized signaling: basolateral receptor localization in epithelial cells by PDZ-containing proteins. Curr. Opin. Cell Biol. 9, 853–859 (1997).
Nallet-Staub, F. et al. Cell density sensing alters TGF-β signaling in a cell-type-specific manner, independent from Hippo pathway activation. Dev. Cell 32, 640–651 (2015).
Murphy, J. E., Padilla, B. E., Hasdemir, B., Cottrell, G. S. & Bunnett, N. W. Endosomes: a legitimate platform for the signaling train. Proc. Natl Acad. Sci. USA 106, 17615–17622 (2009).
Auciello, G., Cunningham, D. L., Tatar, T., Heath, J. K. & Rappoport, J. Z. Regulation of fibroblast growth factor receptor signalling and trafficking by Src and Eps8. J. Cell Sci. 126, 613–624 (2013).
Jean, S. et al. Extended-synaptotagmin-2 mediates FGF receptor endocytosis and ERK activation in vivo. Dev. Cell 19, 426–439 (2010).
Vieira, A. V., Lamaze, C. & Schmid, S. L. Control of EGF receptor signaling by clathrin-mediated endocytosis. Science 274, 2086–2089 (1996).
Bergert, M. et al. Cell surface mechanics gate embryonic stem cell differentiation. Cell Stem Cell 28, 209–216.e204 (2021).
De Belly, H. et al. Membrane tension gates ERK-mediated regulation of pluripotent cell fate. Cell Stem Cell 28, 273–284. e276 (2021).
Schormann, W., Hariharan, S. & Andrews, D. W. A reference library for assigning protein subcellular localizations by image-based machine learning. J. Cell Biol. 219, e201904090 (2020).
Ueki, R. et al. DNA aptamer assemblies as fibroblast growth factor mimics and their application in stem cell culture. Chem. Commun. 55, 2672–2675 (2019).
Dutta, D., Williamson, C. D., Cole, N. B. & Donaldson, J. G. Pitstop 2 is a potent inhibitor of clathrin-independent endocytosis. PLoS One 7, e45799 (2012).
Macia, E. et al. Dynasore, a cell-permeable inhibitor of dynamin. Dev. Cell 10, 839–850 (2006).
Greber, B. et al. Conserved and divergent roles of FGF signaling in mouse epiblast stem cells and human embryonic stem cells. Cell Stem Cell 6, 215–226 (2010).
Greber, B. et al. FGF signalling inhibits neural induction in human embryonic stem cells. EMBO J. 30, 4874–4884 (2011).
Guan, J.-L. Role of focal adhesion kinase in integrin signaling. Int. J. Bochem. Cell Biol. 29, 1085–1096 (1997).
Chandrasekar, I. et al. Nonmuscle myosin II is a critical regulator of clathrin-mediated endocytosis. Traffic 15, 418–432 (2014).
Levayer, R., Pelissier-Monier, A. & Lecuit, T. Spatial regulation of Dia and Myosin-II by RhoGEF2 controls initiation of E-cadherin endocytosis during epithelial morphogenesis. Nat. Cell Biol. 13, 529–540 (2011).
Manenschijn, H. E. et al. Type-I myosins promote actin polymerization to drive membrane bending in endocytosis. eLife 8, e44215 (2019).
Straight, A. F. et al. Dissecting temporal and spatial control of cytokinesis with a myosin II Inhibitor. Science 299, 1743–1747 (2003).
Golding, I., Kozlovsky, Y., Cohen, I. & Ben-Jacob, E. Studies of bacterial branching growth using reaction–diffusion models for colonial development. Physica A 260, 510–554 (1998).
Bell, G. I. Models for the specific adhesion of cells to cells: a theoretical framework for adhesion mediated by reversible bonds between cell surface molecules. Science 200, 618–627 (1978).
Paszek, M. J., Boettiger, D., Weaver, V. M. & Hammer, D. A. Integrin clustering is driven by mechanical resistance from the glycocalyx and the substrate. PLoS Comput. Biol. 5, e1000604 (2009).
Hill, A. V. The possible effects of the aggregation of the molecules of hemoglobin on its dissociation curves. J. Physiol. 40, iv–vii (1910).
Chen, C. S., Mrksich, M., Huang, S., Whitesides, G. M. & Ingber, D. E. Geometric control of cell life and death. Science 276, 1425–1428 (1997).
Minn, K. T. et al. High-resolution transcriptional and morphogenetic profiling of cells from micropatterned human ESC gastruloid cultures. eLife 9, e59445 (2020).
Tewary, M. et al. A stepwise model of reaction-diffusion and positional information governs self-organized human peri-gastrulation-like patterning. Development 144, 4298–4312 (2017).
Turing, A. M. The chemical basis of morphogenesis. Bull. Math. Biol. 52, 153–197 (1990).
Gierer, A. & Meinhardt, H. A theory of biological pattern formation. Kybernetik 12, 30–39 (1972).
Muñoz-Sanjuán, I. & Brivanlou, A. H. Neural induction, the default model and embryonic stem cells. Nat. Rev. Neurosci. 3, 271–280 (2002).
Loh, K. M. & Lim, B. A precarious balance: pluripotency factors as lineage specifiers. Cell Stem Cell 8, 363–369 (2011).
Thomson, M. et al. Pluripotency circuit members mediate germ layer fate choice of embryonic stem cells. Cell 145, 875–889 (2011).
Wang, Z., Oron, E., Nelson, B., Razis, S. & Ivanova, N. Distinct lineage specification roles for NANOG, OCT4, and SOX2 in human embryonic stem cells. Cell Stem Cell 10, 440–454 (2012).
Fagotto, F. & Gumbiner, B. M. Cell contact-dependent signaling. Dev. Biol. 180, 445–454 (1996).
Levine, E. M., Becker, Y., Boone, C. W. & Eagle, H. Contact inhibition, macromolecular synthesis, and polyribosomes in cultured human diploid fibroblasts. Proc. Natl Acad. Sci. USA 53, 350–356 (1965).
Qi, T. et al. Function and regulation of the PEA3 subfamily of ETS transcription factors in cancer. Am. J. Cancer Res. 10, 3083 (2020).
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Miao, Z., Deng, K., Wang, X. & Zhang, X. DEsingle for detecting three types of differential expression in single-cell RNA-seq data. Bioinformatics 34, 3223–3224 (2018).
Abecasis, B. et al. Expansion of 3D human induced pluripotent stem cell aggregates in bioreactors: bioprocess intensification and scaling-up approaches. J. Biotechnol. 246, 81–93 (2017).
Jang, J. et al. Primary cilium-autophagy-Nrf2 (PAN) axis activation commits human embryonic stem cells to a neuroectoderm fate. Cell 165, 410–420 (2016).
Han, D. et al. ZBTB12 is a molecular barrier to dedifferentiation in human pluripotent stem cells. Nat. Commun. 14, 632 (2023).
Acknowledgements
We thank S. T. Baek at Pohang University of Science and Technology for providing a lentiviral vector expressing KRASG12V and W. Kim at Chungnam National University for insightful comments. This work was supported by the National Research Foundation of Korea grant funded by the Korea government (NRF-2020M3A9D8038184, NRF-2021R1A4A1031754, NRF-2022R1F1A1063619 and RS-2023-00221112 to J.J.), by the Korea Institute for Advancement of Technology and the Ministry of Trade, Industry and Energy of the Republic of Korea (P0021109 to J.J.) and BK21 FOUR.
Author information
Authors and Affiliations
Contributions
J.J. was responsible for conceptualization. S.Y. and M.G. were responsible for methodology. S.O., Y.O., Y.C., J.Y. and S.J. were responsible for validation. M.G. and M.A.L. were responsible for formal analyses. S.Y., M.G., S.O., Y.O., Y.C., J.Y. and S.J. were responsible for investigation. S.Y., M.G. and J.J. wrote the original draft. S.Y., M.G., S.O., Y.O., Y.C., J.Y., M.A.L., S-Y.P., Y.L. and J.J. were responsible for review and editing. J.J. was responsible for supervision and funding acquisition.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Cell Biology thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 The expression of ETV4 is regulated by mechanical microenvironments.
a, Relative cell area of H9 hESCs in low and high density cultures. n=10 regions pooled from two independent experiments. b, Volcano plot showing DEGs from H9 hESCs in high density culture. The red and blue dots indicate up- and down-regulated genes, respectively, with cutoff values for DEGs: log2FC < −0.5 or > 0.5, adjP < 0.05. Full list of DEGs can be found in Supplementary Table 1. c, The top 10 biological pathways ranked by p-values from GO analysis with total DEGs d, Gene set enrichment analysis (GSEA) with the RNA-seq data from H9 hESCs in low and high density cultures. GSEA was performed with 3 different gene sets: Actomyosin, Actin filament bundle, Contractile fiber from Molecular Signatures Database (MSigDB). e, Immunofluorescence assay for pMLC in H9 hESCs under low and high density cultures. n=25 regions. f, FPKM of PEA3 family transcription factors in H9 hESCs (GSE183702). n=2 samples pooled from two independent experiments. g, Immunofluorescence assay for ETV4 in ARPE19 cells under low and high density cultures. n=30 cells for low density, n=41 cells for high density. h, Schematic representation of the cell stretching system. i, Colony size measurements in H9 hESCs on the cell stretching system. n=12 colonies. j, Cell area measurements in MCF-7 cells on the cell stretching system. n=12 colonies. k, Immunofluorescence assay for pMLC and OCT4 in H9 hESCs on PDMS layers with different stiffnesses. n=25 colonies for plastic, n=20 colonies for 15 kPa and 1.5 kPa pooled from two independent experiments. l, Immunofluorescence assay for YAP/TAZ and OCT4 in H9 hESCs on PDMS layers with different stiffnesses. n=40 cells for YAP/TAZ pooled from two independent experiments. n = number of cells (g), regions (e), or colonies (i and j) pooled from three independent experiments. Two-sided Student’s t-test, ***P < 0.001, **P < 0.01, *P < 0.05. Exact P values are presented in Supplementary Table 9. Scale bars: 10 µm (e,j), 25 µm (g and l), 50 µm (i and k). Numerical source data are available in Source data.
Extended Data Fig. 2 Cell crowding regulates ETV4 expression in a growing hESC colony.
a, Immunofluorescence assay for ZO1 in large H9 colonies. n=30 regions. b, Correlation of nucleus size with cell area in H9 hESCs. n=15 cells. c,d, Immunofluorescence assay for OCT4, NANOG, SOX2, and ETV4 in H1 and H9 colonies. n=30 cells for (c), n=60 cells for (d). e,f, Immunofluorescence assay for OCT4 and ETV4 in large H9 colonies in mTeSR1 on Matrigel-coated plates (e) or in TeSR-E8 on vitronectin-coated plates (f). n=30 cells. g, Immunofluorescence assay for OCT4 and ETV4 in large H9 colonies derived from single cells. n=30 cells. n = number of cells (b,c,d,e,f, and g) or regions (a) pooled from three independent experiments. Two-sided Student’s t-test, ***P < 0.001, **P < 0.01, *P < 0.05. Exact P values are presented in Supplementary Table 9. Scale bars: 25 µm (a), 50 µm (c,d,e,f, and g). Numerical source data are available in Source data.
Extended Data Fig. 3 ETV4 and YAP exhibit distinct sensitivities to changes in cell density in hESCs.
a, The statistical analysis of cell shape index and relative cell area of ETV4-high and ETV4-low cells in large H9 colonies. n=90 cells. b, Gene set enrichment analysis (GSEA) with the RNA-seq data from H9 hESCs in low and high density cultures. GSEA was performed with 2 different gene sets: Hallmark_Glycolysis, Hallmark_Hypoxia from Molecular Signatures Database (MSigDB). c, Immunofluorescence assay for Hypoxyprobe in small and large H9 colonies. n=9 regions. d,e, Glucose uptake assay in small and large H9 colonies. n=25 regions for 37°C (d), n=31 regions for 4°C (d), and n=30 regions for (e). f, Immunofluorescence assay for YAP in large H9 colonies. YAP activity was measured by the ratio of the nuclear over cytoplasmic intensities (N/C ratio). n=5 colonies. g, Immunofluorescence assay for YAP/TAZ and OCT4 in large H9 colonies. n=80 cells pooled from two independent experiments. h,i, Validation of lentiviral vectors expressing YAP shRNAs in H9 hESCs by qPCR (h) and immunostaining (i). n=3 samples for qPCR (h) and n=50 cells for immunostaining (i). j, Quantification of cells in M phase in large H9 colonies. n=20 regions. k, Immunofluorescence assay for Ki67 in large H9 colonies. n=40 cells pooled from two independent experiments. l, Immunofluorescence assay for ETV4 and EdU in large H9 colonies. n=6 regions. m, Crystal violet staining in H9 hESCs transduced with lentiviral vectors expressing YAP shRNAs. n = number of cells (a and i), colonies (f), or regions (c,d,e,j, and l) pooled from three independent experiments. Two-sided Student’s t-test, ***P < 0.001, **P < 0.01, *P < 0.05. Exact P values are presented in Supplementary Table 9. Scale bars: 25 µm (a,c,d,e,g, and i), 50 µm (l), 100 µm (k). Numerical source data are available in Source data.
Extended Data Fig. 4 ETV4 links cell crowding to NE derepression.
a, Immunofluorescence assay for PAX6 and Brachyury in differentiated H1 colonies. n=6 colonies. b, qPCR analysis for lineage markers in H9 hESCs differentiated by FGF2 and TGF-β deprivation. n=3 samples pooled from three independent experiments. c, RPKM of ETV4 in undifferentiated H9 hESCs, and H9-derived NE and ME cells (GSE69982). n=3 samples. d, Validation of lentiviral vectors expressing ETV4-HA in H9 hESCs by immunostaining. n=15 cells. e,f, Validation of lentiviral vectors expressing ETV4 shRNAs together with GFP in H9 hESCs by qPCR (e) and immunostaining (f). n=3 samples for qPCR (e) from three independent experiments, n=75 cells for shCtrl, n=55 cells for ETV4 sh1, and n=50 cells for ETV4 sh2 for immunostaining (f) pooled from two independent experiments. g, Immunofluorescence assay for PAX6 in shETV4-expressing H9 hESCs transduced with lentiviral vectors expressing ETV4-HA together with GFP and differentiated to NE cells for 5 days. n=4 regions. h, Immunofluorescence assay for Brachyury in H9 hESCs transduced with lentiviral vectors expressing ETV4 shRNA and differentiated to ME cells for 3 days. n=11 regions. i, Immunofluorescence assay for PAX6 in high density H9 hESCs transduced with lentiviral vectors expressing ETV4-HA and differentiated to NE cells for 5 days. n=11 regions. n = number of cells (d), colonies (a), or regions (g,h, and i) pooled from three independent experiments. Two-sided Student’s t-test, ***P < 0.001, **P < 0.01, *P < 0.05. Exact P values are presented in Supplementary Table 9. Scale bars: 10 µm (d), 25 µm (f, and g), 50 µm (a), 100 µm (h and i). Source numerical data are available in Source data.
Extended Data Fig. 5 ETV4 inhibits NE differentiation by MMPs.
a,b,c The top 10 biological pathways ranked by p-values from GO analysis with total DEGs (a), upregulated DEGs (b), and downregulated DEGs (c) after ETV4 KD. d, qPCR analysis for DUSP6 in H9 hESCs transduced with lentiviral vectors expressing ETV4 shRNA or ETV4-HA. n=3 samples pooled from three independent experiments. e, Volcano plot showing DEGs in H1 hESCs sorted by N-cadherin expression. The red and blue dots indicate up- and down-regulated genes, respectively, with cutoff values for DEGs: log2FC < −0.5 or > 0.5, adjP < 0.05 (GSE126022). Full list of DEGs can be found in Supplementary Table 4. f, The top 10 biological pathways ranked by p-values from GO analysis with upregulated DEGs in N-cadherin+ H1 hESCs. g,h, qPCR analysis for MMP2, MMP9, MMP14, and MMP15 in H9 hESCs transduced with lentiviral vectors expressing ETV4 shRNA (g) or ETV4-HA (h). n=3 samples pooled from three independent experiments. i, Immunofluorescence assay for MMP14 in H9 hESCs transduced with lentiviral vectors expressing ETV4-HA. n=25 cells. j, FPKM of MMPs in undifferentiated H9 hESCs (GSE183702). n=4 samples pooled from four independent experiments. k, Immunofluorescence assay for MMP14 in H9 hESCs transduced with lentiviral vectors expressing MMP14 together with GFP. n=50 cells. l, Immunofluorescence assay for PAX6 in H9 hESCs transduced with lentiviral vectors expressing MMP14 together with GFP and differentiated to NE cells for 5 days. n=10 regions. n = number of cells (i and k) or regions (l) pooled from three independent experiments. Two-sided Student’s t-test, ***P < 0.001, **P < 0.01, *P < 0.05. Exact P values are presented in Supplementary Table 9. Scale bars: 10 µm (i,k, and l). Numerical source data are available in Source data.
Extended Data Fig. 6 ERK regulates the protein stability of ETV4.
a, Validation of KTRs treated with specific inhibitors. n=13 cells for ERK-KTR and JNK-KTR, n=8 cells for p38-KTR pooled from two independent experiments. b,c, Representative time-course images from three independent colony-tracking assays for p38-KTR (b) and JNK-KTR (c) in H9 colonies. d, Immunofluorescence assay for pERK in H9 cells treated with PD0325901 (1 μM, 1h). n=30 cells. e, Immunofluorescence assay for pERK in large H9 colonies. Mitotic activation of ERK was seen in brightly stained cells. n=40 cells. f, ERK activity measured by the ratio of the cytoplasmic over nuclear intensities (C/N ratio) in single H9 hESCs expressing ERK-KTR under low and high density cultures (24h after seeding). n=26 cells. g, qPCR analysis for ETV4 in H9 hESCs treated with PD0325901 (1 μM). n=3 samples pooled from three independent experiments. h, Immunofluorescence assay for ETV4 in H9 hESCs treated with U0126 (10 μM, 1h). n=70 cells pooled from two independent experiments. n = number of cells (d,e,f) pooled from three independent experiments. Two-sided Student’s t-test, ***P < 0.001, **P < 0.01, *P < 0.05. Exact P values are presented in Supplementary Table 9. Scale bars: 10 µm (a), 25 µm (b,c,d,f, and h), 100 µm (e). Numerical source data and additional microscope images are available in Source data.
Extended Data Fig. 7 COP1 is involved in the regulation of ETV4 protein stability.
a, Immunofluorescence assay for COP1 in large H9 colonies. n=40 cells. b,c, Validation of lentiviral vectors expressing COP1 shRNA in H9 hESCs by qPCR (b) and immunostaining (c). n=3 samples for qPCR (b) from three independent experiments, n=60 cells for immunostaining (c) pooled from two independent experiments. d, Immunofluorescence assay for ETV4 in H9 hESCs transduced with lentiviral vectors expressing COP1 shRNA and treated with PD0325901 (1 μM, 30 min). n=64 cells. e, Immunofluorescence assay for ETV4 in large H9 colonies transduced with lentiviral vectors expressing COP1 shRNA together with GFP. n=6 regions pooled from two independent experiments. f, Immunofluorescence assay for PAX6 in H9 hESCs transduced with lentiviral vectors expressing COP1 shRNA together with GFP and differentiated to NE cells for 5 days. n=9 regions. g, Immunofluorescence assay for pERK in H9 hESCs transduced with lentiviral vectors expressing KRASG12V together with GFP. n=30 cells pooled from two independent experiments. h, Immunofluorescence assay for ETV4 in H9 hESCs transduced with lentiviral vectors expressing KRASG12V together with GFP. n=63 cells pooled from two independent experiments. i, Immunofluorescence assay for PAX6 in H9 hESCs transduced with lentiviral vectors expressing constitutively-active KRASG12V together with GFP and differentiated to NE cells for 5 days. n=12 regions. n = number of cells (a and d) or regions (f and i) pooled from three independent experiments. Two-sided Student’s t-test, ***P < 0.001, **P < 0.01, *P < 0.05. Exact P values are presented in Supplementary Table 9. Scale bars: 10 µm (g), 25 µm (c,d,e, and h), 50 µm (f and i), 100 µm (a). Numerical source data are available in Source data.
Extended Data Fig. 8 FGF signalling regulates ERK and ETV4.
a, FPKM of EGFR and FGFRs in undifferentiated H9 hESCs (GSE183702). n=4 samples pooled from four independent experiments. b,c, Validation of lentiviral vectors expressing FGFR1 shRNAs in H9 hESCs by qPCR (b) and immunostaining (c). n=3 samples for qPCR (b) and n=33 cells for immunostaining (c). d, Immunofluorescence assay for pERK in H9 hESCs transduced with lentiviral vectors expressing FGFR1 shRNAs together with GFP. n=50 cells pooled from two independent experiments. e, Immunofluorescence assay for ETV4 in H9 hESCs transduced with lentiviral vectors expressing FGFR1 shRNAs together with GFP. n=37 cells pooled from two independent experiments. f, Immunofluorescence assay for ETV4 in small and large H9 colonies cultured in a transwell system. n=30 cells. n = number of cells (c and f) pooled from three independent experiments. Two-sided Student’s t-test, ***P < 0.001, **P < 0.01, *P < 0.05. Exact P values are presented in Supplementary Table 9. Scale bars: 25 µm (c,d,e, and f). Numerical source data are available in Source data.
Extended Data Fig. 9 Cell crowding blocks FGFR endocytosis.
a, Immunofluorescence assay for FGFR1 and EEA1 in small and large H9 colonies. Arrow heads indicate apical localization of FGFR1 proteins. n=6 regions for Small, n=8 regions for Large. b, Immunofluorescence assay for FGFR1 and E-cadherin in the centre of large H9 colonies. n=8 regions. c, Validation of TD0 (500 nM) in human dermal fibroblasts expressing ERK-KTR. n=22 cells pooled from two independent experiments. d, Immunofluorescence assay for FGFR1 and EEA1 in large H9 colonies treated with TD0 (500 nM, 1h). n=9 regions pooled from two independent experiments. e, Immunofluorescence assay for EEA1 in H9 hESCs treated with Pitstop2 (50 μM, 1h) or Dynasore (100 μM, 1h). n=28 regions for NC, n=12 regions for Pitstop2, n=16 regions for Dynasore pooled from two independent experiments. f, Immunofluorescence assay for FGFR1 and EEA1 in H9 hESCs treated with Dynasore (100 μM, 1h). n=8 regions pooled from two independent experiments. g, Time-course images for ERK-KTR in H9 hESCs treated with Dynasore (100 μM). ERK activity measured by the ratio of the cytoplasmic over nuclear intensities (C/N ratio). n=46 cells. h, Immunofluorescence assay for ETV4 in H9 hESCs treated with Dynasore (100 μM, 1h). n=140 cells pooled from two independent experiments. i, Immunofluorescence assay for ETV4 in H9 hESCs transduced with lentiviral vectors expressing KRASG12V together with GFP and treated with Pitstop2 (50 μM, 30 min). n=45 cells. j, Immunofluorescence assay for PAX6 in H9 hESCs transduced with lentiviral vectors expressing FGFR1 shRNAs together with GFP and differentiated to NE cells for 5 days. n=10 regions. n = number of cells (g and i) or regions (a,b, and j) pooled from three independent experiments. Two-sided Student’s t-test, ***P < 0.001, **P < 0.01, *P < 0.05. Exact P values are presented in Supplementary Table 9. Scale bars: 10 µm (e,f, and j), 15 µm (a), 25 µm (d,g,h, and i), 50 µm (b), 100 µm (c). Numerical source data are available in Source data.
Extended Data Fig. 10 Cell crowding inhibits the integrin-actomyosin pathway.
a, The top 10 cellular components ranked by p-values from GO analysis with downregulated DEGs in high density H9 hESCs. b, The top 10 KEGG pathways ranked by p-values from GO analysis with upregulated DEGs in N-cadherin+ H1 hESCs. c, Immunofluorescence assay for pFAK and OCT4 in H9 hESCs treated with PND1186 (2 μM, 1h). n=35 cells for pFAK, n=60 cells for OCT4 pooled from two independent experiments. d, Immunofluorescence assay for FGFR1 and EEA1 in H9 hESCs treated with PND1186 (2 μM, 1h) or YM (10 μM, 1h). n=8 regions for NC and YM, n=9 regions for PND1186. e, Time-course images for ERK-KTR in H9 hESCs treated with PND1186 (2 μM). n=40 cells. f, Immunofluorescence assay for ETV4 and OCT4 in H9 hESCs treated with PND1186 (2 μM, 1h). n=60 cells. g, Immunofluorescence assay for pMLC in H9 hESCs treated with RGDS (500 μM, 24h). n=15 regions pooled from two independent experiments. h, Immunofluorescence assay for pMLC in H9 hESCs treated with PND1186 (2 μM, 1h). n=40 regions pooled from two independent experiments. i, Immunofluorescence assay for pMLC in H9 hESCs treated with Blebbistatin (50 μM, 30 min). n=13 regions for NC and n=10 regions for Blebbistatin pooled from two independent experiments. j, Immunofluorescence assay for pMLC in H9 hESCs treated with YM (10 μM, 1h). n=35 regions pooled from two independent experiments. k, Time-course images for ERK-KTR in H9 hESCs treated with YM (10 μM). n=39 cells. l, Immunofluorescence assay for ETV4 and OCT4 in H9 hESCs treated with YM (10 μM, 1h). n=60 cells. m, Immunofluorescence assay for pMLC and OCT4 in the centre of large H9 colonies transduced with lentiviral vectors expressing constitutively-active RhoAQ63L-GFP in a doxycycline-dependent manner. n=9 regions. n = number of cells (e,f,k, and l) or regions (d and m) pooled from three independent experiments. Two-sided Student’s t-test, ***P < 0.001, **P < 0.01, *P < 0.05. Exact P values are presented in Supplementary Table 9. Scale bars: 10 µm (g, and i), 20 µm (h and j), 25 µm (c,d,e,f,k, and l), 50 µm (m). Numerical source data are available in Source data.
Supplementary information
Supplementary Table 1
DEGs in H9 hESCs under high-density culture.
Supplementary Table 2
Predicted transcriptional regulators for cell crowding.
Supplementary Table 3
DEGs in H9 hESCs after ETV4 KD.
Supplementary Table 4
DEGs in scRNA-seq data from N-cad+ hESCs.
Supplementary Table 5
Antibodies and reagents used in this study.
Supplementary Table 6
Plasmids used in this study.
Supplementary Table 7
qPCR primers used in this study.
Supplementary Table 8
Model variables and value used in this study.
Supplementary Table 9
P values in this study.
Source data
Source Data Figs. 1–7, Extended Data Figs. 1–10
Numerical source data.
Source Data Figs. 1 and 4
Full-length, unprocessed gels.
Source Data Fig. 4, Extended Data Fig. 6
Additional microscope images.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yang, S., Golkaram, M., Oh, S. et al. ETV4 is a mechanical transducer linking cell crowding dynamics to lineage specification. Nat Cell Biol 26, 903–916 (2024). https://doi.org/10.1038/s41556-024-01415-w
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41556-024-01415-w