Mechanotransductive feedback control of endothelial cell motility and vascular morphogenesis

Vascular morphogenesis requires persistent endothelial cell motility that is responsive to diverse and dynamic mechanical stimuli. Here, we interrogated the mechanotransductive feedback dynamics that govern endothelial cell motility and vascular morphogenesis. We show that the transcriptional regulators, YAP and TAZ, are activated by mechanical cues to transcriptionally limit cytoskeletal and focal adhesion maturation, forming a conserved mechanotransductive feedback loop that mediates human endothelial cell motility in vitro and zebrafish intersegmental vessel (ISV) morphogenesis in vivo. This feedback loop closes in 4 hours, achieving cytoskeletal equilibrium in 8 hours. Feedback loop inhibition arrested endothelial cell migration in vitro and ISV morphogenesis in vivo. Inhibitor washout at 3 hrs, prior to feedback loop closure, restored vessel growth, but washout at 8 hours, longer than the feedback timescale, did not, establishing lower and upper bounds for feedback kinetics in vivo. Mechanistically, YAP and TAZ induced transcriptional suppression of RhoA signaling to maintain dynamic cytoskeletal equilibria. Together, these data establish the mechanoresponsive dynamics of a transcriptional feedback loop necessary for persistent endothelial cell migration and vascular morphogenesis.


Introduction
Christiane Nüsslein-Vohlard wrote in her book, Coming to Life, that "during development, any change in cell shape is preceded by a change in gene activity" (Nüsslein-Volhard, 2006).It is intuitive that gene expression is necessary to produce, for example, the cytoskeletal and adhesion machinery that enable cell motility, but whether persistent gene expression is required for continued migration is not.We previously found, to our surprise, that blockade of de novo gene expression (i.e., transcription or translation), caused progressive motility arrest (Mason et al., 2019).This requirement for continued transcription is executed by a feedback loop mediated by the transcriptional regulators, Yes-associated protein (YAP) and transcriptionalcoactivator with PDZ-binding motif (TAZ) (Mason et al., 2019).Here, we define the kinetics of this mechanotransductive feedback loop, in vitro and in vivo.
Mechanotransduction is often evaluated after the cell has reached equilibrium and represented as a one-way path from cytoskeletal activation to nuclear information transfer, resulting in a state variable shift (e.g., lineage commitment).However, mechanotransduction can also activate feedback loops, which alter how the cell responds to subsequent mechanical stimuli.These feedback loops are necessary to maintain a responsive cytoskeletal equilibrium or to preserve an adapted state.
Recently, we and, independently, the Huveneers group found that human endothelial cells maintain cytoskeletal equilibrium for persistent motility through a YAP/TAZ-mediated feedback loop that maintains tensional homeostasis (Mason et al., 2019;van der Stoel et al., 2020).Because YAP and TAZ are activated by tension of the cytoskeleton (Dupont et al., 2011), suppression of cytoskeletal tension by YAP/TAZ transcriptional target genes constitutes a negative feedback loop that acts as a control system to maintain cytoskeletal homeostasis for persistent motility via modulation of Rho-ROCK-myosin II activity (Fig. 1A) (Mason et al., 2019).Both we and the Huveneers group found that perturbation of genes and pathways regulated by YAP/TAZ mechanotransduction can functionally rescue motility in YAP/TAZ-depleted cells (e.g., RhoA/ROCK/myosin II, NUAK2, DLC1) (Mason et al., 2019;van der Stoel et al., 2020).These prior findings establish functional nodes of this mechanotransductive feedback loop.
Here, we sought to determine the characteristic time scales and physiologic impacts, using two orthogonal model systems: human endothelial colony forming cells in vitro and zebrafish intersegmental vessel development in vivo.Endothelial colony forming cells (ECFCs) are circulating endothelial cells that activate in response to vascular injury, adhering and migrating to the injury site to repopulate the damaged endothelium (Ingram et al., 2004;David A Ingram et al., 2005).To study the dynamics of cytoskeletal feedback control in ECFCs, we performed dynamic quantification of cytoskeletal and cellular morphodynamics over time after adhesion.To study the dynamics of vascular morphogenesis, we performed longitudinal, quantitative imaging of intersegmental vessel (ISV) development.This model is particularly tractable for this question as zebrafish ISV development exhibits similar migratory kinetics to human ECFC migration in vitro and can be quantified longitudinally in fluorophore-labeled endothelial cells in vivo (Phng et al., 2013;Rosa et al., 2022).
Here, we define the characteristic time scales of YAP/TAZ-mediated transcriptional feedback loop closure and demonstrate conservation of feedback loop mediators and kinetics between human endothelial cell migration in vitro and zebrafish vascular morphogenesis in vivo.
Movie 1: Representative effects of transcription inhibition on ECFC migration and cytoskeletal dynamics.Adherent ECFCs, transfected with LifeAct-tdTomato, were treated with DMSO (left) or the transcription inhibitor, Actinomycin D (right).Images were taken in 3-minute intervals for 12.5 hours after inhibitor treatment.Time is shown in the top left as hours:minutes.

Transcriptional feedback control of cytoskeletal and focal adhesion maturation
To test how de novo gene expression mediates endothelial cell motility, we evaluated live migration of human ECFCs for 10 hours after broad transcription inhibition by Actinomycin D. Control cells exhibited dynamic actin stress fiber formation and remodeling (via LifeAct-tdTomato) and persistent motility, while transcription-inhibited cells initiated with similar dynamics, but exhibited stress fiber stabilization and motile arrest between four and eight hours after transcription blockade (Movie 1).Next, we evaluated the effects of transcription-inhibited stress fiber arrest on the formation and maturation of focal adhesions at twenty-four hours after in situ Actinomycin D treatment (Fig. 1B-D).In situ transcription inhibition decreased the total number of vinculin+ adhesion plaques (Fig. 1C) and increased the fraction of adhesions over 1 µm in length, a measure of focal adhesion maturity (Fig. 1D).As with broad transcription inhibition, we previously found that in situ depletion of YAP and TAZ by RNAi arrested cell motility, illustrated here by live-migration sparklines over 10 hours: siControl: , siYAP/TAZ: (25 μm scale-bar: -).(Mason et al., 2019) In situ YAP/TAZ depletion decreased the total number of vinculin+ adhesion plaques, and increased the fraction of mature adhesions (Fig. 1E-G).Control cells exhibited fan-shaped lamellipodial migration while both transcription inhibition and YAP/TAZ depletion impaired polarization, and induced robust ventral stress fiber formation and peripheral focal adhesion maturation.Together with our prior findings (Mason et al., 2019), these data demonstrate that continued gene expression, via YAP/TAZ-mediated transcription, maintains cytoskeletal equilibrium for persistent motility.

Mechanoregulation of endothelial cell morphodynamics and persistent motility
To study the dynamics of cytoskeletal feedback control in ECFCs, we developed an adhesion-spreading-polarization-migration (ASPM) assay (Fig. 2A), following (Benecke et al., 1978).This assay corresponds to ECFC adhesive function in vivo and features cell trypsinization and reattachment at t0 to re-activate the cytoskeleton from a partially reset state.Upon adhesion, ECFCs initially spread isotropically, then polarize and transition to randomly-directed motility (Movie 2).
Next, to determine how matrix mechanosensing impacts ECFC motility and morphodynamics, we tracked migration speed, spreading, and shape of mTomatoexpressing ECFCs on 1, 8, or 18 kPa MeHA matrices from the time of attachment to 24 hours post-adhesion (Fig. 2I-K, Movies 3, 4).Regardless of matrix rigidity, cell motility and morphodynamics seemed to exhibit an inflection point at 4 hours post-attachment, and reache motile and morphological equilibrium by hour 8. To explore this equilibrium, we measured the difference in motility and morphology parameters between 8 hours and time zero (∆8-0).
Movies 3, 4: Effects of transcription inhibition on ECFC migration on soft, moderate, and stiff matrices.mTomoto-expressing ECFCs were treated with DMSO (Movie 3) or the transcription inhibitor, Actinomycin D (Movie 4).Images were taken in 3-minute intervals for 12.5 hours after inhibitor treatment.Time is shown as hours:minutes.
Prior studies show that, at equilibrium, cells generally migrate faster on stiff substrates than on soft, and from softer toward stiffer matrices (Lange and Fabry, 2013;Sunyer et al., 2016;Sunyer and Trepat, 2020).However, ECFC migration speed was initially higher on soft substrates than on moderate or stiff substrates (average initial speed = 19.40,9.00, and 8.93 µm/hr, for 1, 8, and 18 kPa respectively; Fig. 2I).Cell motility on soft hydrogels initially increased with time, but inflected at 4 hours and reached an equilibrium low by 8 hours.In contrast, cells on moderate and stiff hydrogels increased continuously, reaching an equilibrium high by 8 hours postadhesion (∆8-0 cell speed = -9.71,+12.92, and +12.39 µm/s for 1, 8, and 18 kPa respectively).Cell morphology dynamics exhibited similar kinetics with inflection points at 4 and 8 hours after attachment, reaching morphological equilibrium by hour 8 (∆8-0 cell area = +197, +525, and +720 µm 2 ; for 1, 8, and 18 kPa respectively; Fig. 2J).Similarly, cell circularity (inverse of elongation) decreased as a function of substrate rigidity, reaching morphological equilibrium by hour 8 (∆8-0 circularity = -0.12,-0.26, and -0.24 for 1, 8, and 18 kPa respectively; Fig. 2K).These data implicate extracellular matrix mechanical cues as regulators of mechanotransductive, morphological, and motile adaptation.Transcriptional control of vascular morphogenesis Mechanical stimuli and endothelial cell mechanotransduction are also critical for proper vascular morphogenesis in vivo (Boselli et al., 2015).Therefore, we next asked whether the YAP/TAZ-mediated cytoskeletal-transcriptional feedback loop identified in human EC migration in vitro is conserved in vivo.We measured dynamic vascular morphogenesis using embryonic zebrafish intersegmental vessel (ISV) migration using live imaging of genetically-encoded GFP-labeled endothelial cells (Tg(fli:egfp) y1 ) (Fig. 3A) (Lawson and Weinstein, 2002).First, we exposed zebrafish embryos at 29 hours post-fertilization (hpf) to Actinomycin D and Puromycin to block general transcription and translation, respectively.Control embryos were treated with vehicle DMSO.Both Act.D (at 25 µg/ml (high), but not 10 µg/ml (low)) and Puromycin treatment slowed the kinetics of ISV growth by 35 hpf (Fig. 3C-E).Next, to inhibit YAP/TAZ signaling, we treated embryos with verteporfin (VP).Previously, we showed that VP treatment of ECFCs in vitro caused progressive cytoskeletal tension, focal adhesion maturation, and motility arrest equivalent to RNAi-mediated depletion (Mason et al., 2019), supporting use of this compound as an endothelial YAP/TAZ inhibitor.Exposure of zebrafish embryos to 50 µM VP at 29 hpf arrested vascular morphogenesis by 35 hpf (Fig. 3F-H).
To block upstream YAP/TAZ activation by RhoA-ROCK-myosin signaling, we also treated with ROCK (Rockout) and non-muscle myosin II (Blebbistatin) inhibitors.Like VP, both ROCK and myosin inhibition abrogated ISV growth at 32 and 35 hpf (Fig. 3F- H).Together, these data confirm that the cytoskeletal-transcriptional feedback loop that mediates persistent human ECFC migration is conserved in zebrafish vascular morphogenesis.Next, we tested how the mechanical environment influences the dynamics of feedback-inhibited motility arrest in vitro.To this end, we performed acute transcription inhibition at the time of cellular adhesion to MeHA hydrogels of 1, 8, or 18 kPa and tracked live cell migration, morphodynamics, and cytoskeletal and adhesion formation and maturation over 24 hours (Fig. 4A, Movies 3, 4).Matrix mechanotransduction responds with sigmoidal behavior to varied substrate rigidity (Cosgrove et al., 2016), and we observed similar dynamic responses to 8 and 18 kPa matrices, indicating these values both sit near or after the responsive plateau (Supplementary Fig. S1).For clarity and simplicity, we show responses to 1 and 18 kPa matrices (Fig. 4) and all groups in the supplement (Supplementary Fig. S1).
To investigate the effects of transcriptional feedback blockade on focal adhesion formation and remodeling, we treated cells with Act.D or DMSO at the time of adhesion to 1, 8, or 18 kPa matrices, and stained F-actin with phalloidin and coimmunostained for vinculin at 4 hours post-adhesion (Fig. 4E-G).This represents the time of active matrix mechanotransduction (cf. Figure 2I) and morphodynamic inflection (cf. Figure 4B-D).Acute transcription inhibition increased the total number of vinculin+ adhesion plaques on 1 and 8 kPa matrices (Fig. 4F), but at this time point did not alter the fraction of mature adhesions on any stiffness (Fig. 4G).We next tested how matrix mechanical properties altered the dynamics of motile and morphodynamic arrest after YAP/TAZ depletion.To this end, we knocked down expression of both YAP and TAZ, by siRNA transfection at 24 hours prior to trypsinization and adhesion to MeHA hydrogels of 1, 8, or 18 kPa (Supplementary Fig. S2A,B).We then tracked cell migration, morphodynamics, and cytoskeletal and adhesion formation and maturation over 24 hours (Fig. 5A, Supplementary Fig. S2C- H, Movies 5, 6).
To investigate the role of mechanotransductive feedback in focal adhesion formation and remodeling at the time of active matrix mechanotransduction (cf. Figure 2I) and motility inflection (cf. Figure 5B), we stained F-actin with phalloidin and coimmunostained for vinculin at 4 hours after adhesion to 1, 8, or 18 kPa matrices, with or without YAP/TAZ depletion (Fig. 5E-G).On 18 kPa matrices, YAP/TAZ depletion decreased the total number of vinculin+ adhesion plaques (Fig. 5F), and increased the fraction of adhesions over 1µm in length (Fig. 5G), but did not significantly alter adhesion number or maturity at this time point on either 1 or 8 kPa.
Together, these data indicate that YAP and TAZ transcriptionally regulate the cytoskeleton with a characteristic time scale of ~4 hours for initial feedback loop closure (tic), resulting in full motile equilibrium by ~8 hours (teq).(Fig. 5H).To determine whether mechanotransductive feedback kinetics are conserved in vivo, we performed two feedback-inhibitor (i.e., Act.D/Puromycin/VP) washout experiments, in which feedback inhibitors were washed out either before (i.e., 3 hrs), or after (i.e., 8 hrs) the characteristic time of initial feedback loop closure (tic).The 8 hr washout corresponds to our previously-identified time-to-equilibrium (teq).We hypothesized that, if feedback loop dynamics are conserved, washout at 3 hours would restore morphogenesis while washout after 8 hours would fail to recover.
In the first experiment, we either exposed zebrafish to inhibitors continuously from 29-35 hpf (Fig. 6A), or restored transcriptional feedback at 3 hours by inhibitor washout (Fig. 6B).Without washout, all three drugs significantly slowed ISV growth by 35 hpf; however, washout at 32 hpf significantly restored vessel growth (p = 0.4, 0.12, 0.99 vs. DMSO) (Fig. 6C-G).In contrast, consistent with our hypothesis, washout at 8 hours failed to restore vessel growth kinetics (Fig. 6I-J).These data are consistent with the feedback kinetics observed in human ECFCs in vitro and confirm functional lower and upper bounds for YAP/TAZ transcriptional feedback kinetics in vivo.We next sought to determine how feedback history influences intrinsic state at the time of ECFC attachment to determine cytoskeletal and adhesion dynamics.Unlike acute treatment with Act.D at the time of adhesion (cf. Figure 4), YAP/TAZ depletion by siRNA requires ~24 hours post-RNAi for protein depletion, but is progressive.Thus, cells trypsinized and replated at 24 hours may already exhibit altered feedback state, analogous to prior studies identifying YAP and TAZ as mediators of cellular mechanical memory (Nasrollahi et al., 2017;Yang et al., 2014a).Therefore, we directly tested the consequences of altering the duration of transcriptional feedback arrest prior to adhesion.
First, we performed either acute or prolonged transcription inhibition (Fig. 7A).Acute transcription inhibition featured Act.D treatment at the time of re-adhesion, as in Figure 4, while prolonged transcription inhibition featured 20 hours of in situ Act.D treatment prior to trypsinization and re-adhesion, for a total of 24 hours.We evaluated F-actin and vinculin staining at 4 hours after re-adhesion (Fig. 7B).Prolonged, but not acute, transcription inhibition significantly decreased focal adhesion number at 4 hours after adhesion (Fig. 7C), while both acute and prolonged transcription inhibition significantly increased the fraction of mature focal adhesions (Fig. 7D).Notably, prolonged Act.D treatment produced cells with dense actin stress fibers, robust mature focal adhesions, and decreased spread area after re-adhesion (Fig. 7C-D).
In parallel, we performed either short-term or long-term YAP/TAZ depletion (Fig. 7E).Short-term YAP/TAZ depletion featured siRNA transfection at 24 hours prior to trypsinization and re-adhesion, as in Figure 5, while long-term transcription inhibition featured siRNA transfection at 48 hours prior to trypsinization and re-adhesion, allowing 24 hours of additional mechanotransductive feedback arrest.We evaluated Factin and vinculin staining at 4 hours after re-adhesion (Fig. 7F).Short-term YAP/TAZ depletion significantly increased focal adhesion number at 4 hours after adhesion (Fig. 7G).On the other hand, long-term YAP/TAZ depletion significantly decreased focal adhesion number at 4 hours after adhesion (Fig. 7G), but increased the fraction of mature focal adhesions (Fig. 7H).Together, these data suggest that the kinetics of cytoskeletal tension generation and focal adhesion formation and remodeling are dependent on the accumulated cytoskeletal state, modulated by mechanotransductive feedback history.protein quantification of MLC phosphorylation by immunoblot (A) and immunofluorescent imaging (B, C).Samples used for quantification of pMLC or MLC in (A) were run and transferred in parallel and GAPDH was used as internal loading control so that signal could be normalized between blots.n = 3, * p < 0.05, Student's two-tailed unpaired t-test.Data shown as mean ± S.E.M. D-G.ECFCs were either pretreated with actinomycin D for 24 hours (D, E) or depleted of YAP and TAZ for 48 hours (F, G) prior to plating on glass coverslips.All cells were fixed for 3D confocal imaging at 10 minutes after attachment.D, F. Representative immunofluorescent images of basal (x-y) actin, pMLC, and nuclei with Alexa Fluor 647-conjugated phalloidin (magenta), Alexa Fluor 488-conjugated secondary (yellow), and DAPI (blue).Side views are orthogonal projections demonstrating cell height (y-z).E, G. Paired measurements of distance of apicalmost pMLC signal above the nucleus after Act.D treatment (E) or YAP/TAZ depletion (G).n = 27-33, ** p < 0.003, **** p < 0.0001, student's two-tailed t-test.Data are shown as mean ± S.E.M. H. Representative apical pMLC ring formation adjacent to the cell membrane in Y/T-depleted cell (cross-section corresponds to dotted line in panel H).I, J. ECFCs were depleted of YAP and TAZ for 48 hours and treated with low dose (15 µM) Blebbistatin for 1 hour.I. Representative immunofluorescent images of actin and nuclei with Alexa Fluor 488-conjugated phalloidin (grey) and DAPI (blue), respectively, depicting ECFC with and without an actin cap.J. Number of cells with an organized, disorganized, or no perinuclear actin either on the basal or apical side of the cell.n = 147-150, p < 0.04, Chi square test with Bonferroni's post hoc test.Data are shown as mean ± S.E.M or as cell number separated into categories.
Next, we observed that in situ Actinomycin D treatment for 24 hours increased myosin light chain (MLC) phosphorylation by 3.6-fold (p < 0.04) and promoted myosin association with filamentous structures (Fig. 8A-C), consistent with our prior findings with YAP/TAZ depletion (Mason et al., 2019).To determine how extended transcriptional feedback disruption regulates acute actomyosin activity and function, we performed Actinomycin D treatment for 24 hours or YAP/TAZ depletion for 48 hours, prior to re-adhesion and spreading for 10 minutes.Control cells exhibited minimal stress fibers, with pMLC primarily localized adjacent to the membrane, indicating that trypsinization partially resets the cytoskeleton.In contrast, extended prior transcriptional feedback disruption, by Actinomycin D (Fig. 8D-E) or YAP/TAZ depletion (Fig. 8F-G) prior to trypsinization and adhesion, disrupted cell spreading and exhibited precocious apical myosin activation above the nucleus (Fig. 8D,F y-z projection).To quantify this, we measured the distance of apical pMLC signal above the nucleus (Fig. 8E,G).Prolonged YAP/TAZ depletion had no effect on nuclear height (2.04 vs 2.06 µm; siControl vs. siYAP/TAZ; p > 0.97), but the distance between the top of the nucleus and apical-most pMLC signal increased 3-fold after YAP/TAZ depletion (0.97 vs 3.02 µm; siControl vs. siYAP/TAZ; p < 0.0001).Actinomycin D treatment for 24 hours increased nuclear height 32% (1.73 vs. 2.28 µm; DMSO vs. Act.D; p < 0.007D) and nucleus to apical pMLC distance 59% (1.21 vs. 1.91 µm; DMOS vs. Act.D; p < 0.003).Cross-sectional visualization at the apex of the nucleus in YAP/TAZ inhibited cells revealed pMLC/actin ring accumulation (Fig. 8H).
Myosin activation has been shown to regulate perinuclear actin filament assembly to enable force transmission from focal adhesions to the nucleus (Chambliss et al., 2013), and nuclear strain regulates nuclear permeability, chromatin architecture, and nuclear repositioning during cell motility (Andreu et al., 2021;Elosegui-Artola et al., 2017;Heo et al., 2016;Lee et al., 2007).Therefore, we next tested how YAP/TAZ feedback regulates perinuclear actin filament assembly into ventral stress fibers below the nucleus (basal) and ventral actin cap fibers above the nucleus (apical) (Fig. 8I, J).YAP/TAZ depletion increased the number of cells with an organized apical actin cap 3fold (p < 0.0001) and basal actin cap by ~30% (p < 0.0001).To determine whether myosin activity is responsible for the increase in perinuclear actin filaments we treated ECFCs with 15 µM Blebbistatin for 1 hour.This low dose of Blebbistatin has been shown to preferentially target apical actin (Chambliss et al., 2013).Blebbistatin treatment of YAP/TAZ-depleted ECFCs reduced the number of cells with organized apical actin by > 40% (p < 0.0001), without significantly affecting basal actin (p = 0.08).Consistently, YAP/TAZ-mediated transcriptional feedback regulated cell spreading and polarization (Movie 7, 8).Together these data suggest that disrupted transcriptional feedback caused precocious myosin activation to produce apical membrane constriction and alter early cell spreading during peri-nuclear stress fiber polymerization.

YAP and TAZ regulate transcriptional programs associated with vascular morphogenesis and cytoskeletal regulation
Control and YAP/TAZ depleted cells were seeded on 18kPa MeHa matrices and harvested after four hours to capture early YAP/TAZ-regulated transcriptional changes.Among the genes most significantly downregulated by YAP/TAZ depletion were YAP, TAZ (WWTR1), and their known targets (CTGF, CYR61, ANKRD1, AXL) (Dupont et al., 2011;Totaro et al., 2018) (Fig. 9A), suggesting significant functional depletion of YAP and TAZ.YAP/TAZ depletion caused a similar number of genes to be downregulated (3266 genes) as upregulated (3370 genes) (Fig. 9A).Genes downregulated upon YAP/TAZ depletion were enriched for biological processes including DNA replication, mitosis, RNA processing, and cell cycle (Fig. 9B).Genes upregulated upon YAP/TAZ depletion were enriched for biological processes including vasculogenesis, response to VEGF stimulus, Rho protein signal transduction, cell migration, and cell adhesion (Fig. 9B).
Previously, we found that YAP and TAZ mediate negative feedback regulation of the cytoskeleton through transcriptional regulation of RhoA activity.(Mason et al., 2019) We observe this negative regulation in this transcriptomic dataset: YAP/TAZ depletion upregulated genes associated with Rho protein signal transduction (Fig. 9B,  C) and altered the expression of genes that directly regulate RhoA activity (Fig. 9D), including previously identified RhoA GAP genes: DLC1 (ARHGAP7), ARHGAP28, and ARHGAP29.(Mason et al., 2019;van der Stoel et al., 2020) Although YAP/TAZ depletion/inhibition impaired functional cell motility and vascular morphogenesis (cf.Figures 3, 5), genes upregulated by YAP/TAZ depletion were enriched for GO terms associated with vascular development, cellular response to vascular endothelial growth factor (VEGF) stimulus, and endothelial cell migration (Fig. 9B, E).This appears contradictory but is consistent with our model of YAP/TAZ mechanotransduction as a negative feedback control loop that maintains cytoskeletal homeostasis for persistent endothelial cell motility and blood vessel development.YAP/TAZ depletion upregulated receptors to angiogenic growth factors, including KDR and FLT4, suggesting YAP/TAZ depleted cells may be more sensitive to VEGF stimulation but remain nonmotile due to cytoskeletal arrest.Together, these data support a role for YAP and TAZ as negative feedback mediators that maintain cytoskeletal homeostasis for endothelial cell migration and vascular morphogenesis.

Discussion
Here, we show that the transcriptional regulators, YAP and TAZ, execute a conserved mechanotransductive feedback loop that mediates human endothelial cell motility in vitro and zebrafish intersegmental vessel morphogenesis in vivo.The feedback loop initially closes in 4 hours, achieving cytoskeletal equilibrium by 8 hours, with conserved kinetics in vivo.YAP/TAZ transcriptional feedback mediates the cellular response to substrate stiffness mechanosensing and vascular morphogenesis, with mechanotransductive feedback history controlling cytoskeletal and adhesion morphodynamic response to a subsequent mechanical input.
YAP and TAZ are critical regulators of morphogenesis (Dong et al., 2007;Pan, 2010;Phillips et al., 2022), including development of the vertebrate vasculature: YAP and TAZ regulate sprouting angiogenesis and vascular development through both cell-autonomous (Kim et al., 2017;Neto et al., 2017) and cell non-autonomous (Choi et al., 2015;Wang et al., 2017) mechanisms.YAP and TAZ are activated by mechanical cues in vivo, including both luminal shear stress and abluminal stretch (Nakajima et al., 2017;Ruehle et al., 2020), and mediate vascular remodeling and low shear stressinduced atherogenesis (Wang et al., 2016).YAP and TAZ also mediate developmental feedback loops.For example, medaka fish embryo morphogenesis against gravitational forces requires YAP/TAZ-mediated transcriptional suppression of RhoA activity through the ARHGAPs (Porazinski et al., 2015).
We identify the mechanotransductive feedback dynamics that mediate human endothelial cell motility in vitro and zebrafish vascular morphogenesis in vivo.Both broad transcription blockade and targeted YAP/TAZ depletion caused morphodynamic change within 4 hours, altering focal adhesion maturation, followed by motility arrest within 8 hours.Our data suggest that the effectors and kinetics of this feedback loop are conserved in zebrafish.Washout of feedback inhibitors prior to initial feedback loop closure, at 3 hours, restored vessel growth.However, inhibitor washout after 8 hours, longer than the feedback loop timescale, prevented morphogenic rescue.These characteristic time scales are consistent with other studies of human EC motility dynamics in response to diverse mechanical cues (Cai and Schaper, 2008;Weijts et al., 2018).Together, these data suggest a conserved mechanotransductive feedback loop with a characteristic time scale of 4 hours for initial feedback loop closure (tic), resulting in full motile equilibrium by 8 hours (teq).
The rate-limiting molecular effectors that determine initial feedback loop completion and equilibrium kinetics are unclear.YAP and TAZ are activated by tension of the actin cytoskeleton and imported into the nucleus from the cytoplasm within minutes through tension-opened nuclear pores (Elosegui-Artola et al., 2017;García-García et al., 2022).Consistent with cytoskeletal remodeling by mechanical perturbation methods (Webster et al., 2014), acute optogenetic activation of either RhoA itself or Rho-activating GEFs induces actomyosin contractility within seconds and YAP nuclear localization within minutes (Berlew et al., 2021;Oakes et al., 2017;Valon et al., 2017).Thus, cytoskeletal activation is likely not rate-limiting.We postulate that the dynamics of YAP/TAZ-dependent transcription and the translation of those target genes are rate-limiting for initial feedback loop completion (tic = 4 hours).This is supported by work from us and others in various cell lines showing YAP/TAZ transcriptional responses occur during the first few hours after activation.(Franklin et al., 2020;Mason et al., 2019;Plouffe et al., 2018) The molecular effectors that determine the timescale of new cytoskeletal equilibrium establishment (teq = 8 hours) remain unclear.
Our new data support and extend prior models of endothelial cell mechanotransductive feedback by YAP and TAZ (Mason et al., 2019;van der Stoel et al., 2020).However, YAP and TAZ regulation in endothelial cells is context-dependent.Disruption of YAP/TAZ signaling generally inhibits cell motility and angiogenesis (Neto et al., 2018;Ong et al., 2022;Sakabe et al., 2017), but promotes sprouting network formation in the bone marrow endothelium (Sivaraj et al., 2020).We speculate that the role of YAP/TAZ mechanotransductive feedback will be determined by the intrinsic cytoskeletal setpoint in a given endothelial bed.Future research will be necessary to test this hypothesis.
Mechanotransduction is often represented as a one-way path from stimulation to state change, but here we found that mechanotransductive feedback history influences the cytoskeletal response to subsequent stimulation.To directly test the effect of cytoskeletal feedback history, we varied the duration of prior feedback arrest, by either transcription inhibition or YAP/TAZ depletion, and evaluated post-detachment adhesion, spreading, polarization, and migration dynamics.We found that extended feedback arrest prior to reattachment, either by transcription inhibition or YAP/TAZ depletion, altered subsequent cell morphodynamics.This history-dependence is distinct from YAP/TAZ-mediated mechanical memory (Mathur et al., 2020;Nasrollahi et al., 2017;Price et al., 2021;Yang et al., 2014b) but suggests that mechanotransduction is non-linear and features, at minimum, second-order feedback dynamics.

Limitations
Here we show that both global blockade of de novo gene expression and specific blockade of YAP/TAZ signaling consistently lead to motility arrest in vitro and in vivo.However, YAP and TAZ are not the only mechanosensitive transcriptional regulators that can modulate the cytoskeleton, and other mechanotransducers are likely to mediate cytoskeletal and morphodynamic feedback.For example, the transcriptional co-activator, MRTF, can also be activated by mechanical cues (Posern et al., 2002), regulates cell migration through transcriptional regulation of cytoskeletal proteins (Leitner et al., 2011), and can co-regulate cytoskeletal feedback in a transcriptional complex with YAP (Katschnig et al., 2017) and/or TAZ (Speight et al., 2016).Future studies will be required to explore roles for this, and other mechanotransducers, in mechanotransductive feedback control of angiogenesis.
Zebrafish embryos enable longitudinal and quantitative imaging of the endothelium during vascular morphogenesis (Phng et al., 2013).Further, embryo transfer between tanks containing inhibitors provides a simple model system for dynamic perturbation of cell signaling.However, every inhibitor has its "off-target" effects, and global inhibition can cause both cell-autonomous and non-cell-autonomous effects on morphogenesis.Additionally, long-term exposure to the pharmacologic inhibitors could impair general cellular function or viability.Future studies will overcome these challenges by orthogonal methods, including inducible and cell-type-specific genetic approaches (Colijn et al., 2022;Pillay et al., 2022) and optogenetic tools (Benman et al., 2022).
pLenti-Lifeact-ubiquitin-tdTomato (Lifeact-tdTomato; addgene: #64048), a filamentous actin-binding peptide labelled with the fluorescent protein tdTomato, was transiently overexpressed in ECFCs under the control of the ubiquitin promoter (Lim et al., 2015).Briefly, ECFCs were plated at a density of 16,000 cells/cm 2 24 hours prior to transfection in antibiotic free full media.Lifeact-tdTomato was diluted to 4 ng/µL in EBM2, 105 ng/cm 2 of DNA per cell culture area.X-tremeGENE™ HP DNA Transfection Reagent (Roche) was used according to the manufacturers protocol, 4µL of transfect reagent was added per 1µg of DNA then incubated for 20-30 minutes at room temperature.Transfection reagent and DNA were then added to ECFCs.ECFCs were used for live imaging 24-72 after transfection.
ECFCs were depleted of YAP and TAZ using custom siRNA (Dharmacon; Table 1) loaded lipofectamine RNAimax (Invitrogen) according to the manufacturer's instructions, as previously described (Mason et al., 2019).Briefly, ECFCs were seeded on collagen coated 6 well-plates, 10 5 cells per well, in antibiotic free medium and kept in culture for 24 hours followed by transfection at approximately 50% confluence.Transfection was carried out using a final concentration 0.3% (v/v) lipofectamine RNAimax with 15 pmol RNAi duplexes (custom oligonucleotides; Dharmacon) per well.Transfected ECFCs were used for downstream experiments 24-48 hours posttransfection.
Table 1: siRNA sense and antisense sequences (Dupont et al Nature 2011, Mason et al JCB 2019)
MeHA hydrogels were polymerized on thiolated #1.5 glass coverslips.Coverslips were washed for 20 minutes in 10 M sodium hydroxide (Alfa Aesar) then washed twice in deionized water and dried.Dried coverslips were incubated in toluene (Fisher Scientific) with 12.5% v/v (3-mercaptopropyl) trimethoxysilane (Sigma Aldrich) and 4.2% v/v hexylamine (Sigma Aldrich) for 1 hour.Thiolated coverslips were washed twice in toluene, dried at 100° C, and stored under nitrogen for up to one month prior to use.
MeHA hydrogels were formed from a 3% w/v MeHA macromer solution functionalized with 1 mM RGD peptide (GenScript) in 0.2 M Triethanolamine (TEOA; Sigma Aldrich) buffer for 30 minutes at room temperature.The dithiol cross-linker Dithiothreitol (DTT; Thermo Scientific) was added to the macromer solution at concentrations of 2.4, 4.8, or 9.6 mM.Sufficient volume to produce 60 µm thick hydrogels were pipetted onto thiolated coverslips and flattened with Rain X treated hydrophobic coverslips.MeHA was polymerized at 37° C for 3 hours then washed PBS and stored at 4° C for cell culture.
MeHA hydrogel elastic modulus was measured using atomic force microscopy (AFM; Bruker Bioscope Catalyst).A 1 µm diameter SiO2 spherical probe (~0.068N/m stiff; Novascan) was indented ~0.5 µm to generate a force-displacement curve.Average elastic modulus was calculated using the hertz contact model formula.Measurements were taken from 10 points for each sample, four samples were measured per hydrogel formulation.

ECFC live migration
ECFC motility was tracked live on MeHA by fitting custom circular PDMS (Dow Corning; 1.75 cm diameter) wells to thiolated coverslips using PDMS to bond the coverslips and glass.Hydrogels were formed in the wells as described above and used for live imaging.For imaging cell attachment mTomato-expressing ECFCs were plated 1-2,000 cells/cm 2 on MeHA and imaged live using a Zeiss Axio Observer inverted microscope with an automated stage for 4 hours at 3-minute intervals followed by 20 hours at 15 minute intervals.Cells were tracked across 5 ROI's for each substrate 5-15 cells per ROI.
For continuously attached cells treated with inhibitors mTomato-expressing ECFCs were first attached for 24-48 hours to collagen coated glass chamber slides before imaging for 24-36 hours at 15-minute intervals.Actinomycin D (Sigma Aldrich), Puromycin (Takara Bio) or DMSO (Sigma Aldrich) were added at the indicated times.During live imaging cells were maintained at 37° C, 5% CO 2 , 95% relative humidity using a stage top type incubation chamber (Heating insert P; Pecon).For all live imaging experiments migration was tracked in 5 regions per condition.For imaging experiments longer than 1 hour a layer of silicone oil (ibidi) was added on top of culture wells to limit evaporation.Cells were tracked across 5 ROI's 25-50 cells per ROI.
Morphological and positional tracking of individual cells was performed using the ADAPT plugin for the open source image analysis software FIJI (Barry et al., 2015).Using ADAPT, the cell outlines were thresholded either using the Li, Huang, or Default thresholding method, depending on tracking fidelity.Objects below a realistic size threshold (200 or 400 µm 2 , depending on the experiment) or that could not be tracked for > 80% of the imaging time were not recorded.Object smoothing was set to .5 and 4 erosion iterations were used to accurately track cell boundaries as a function of time.

Lamellipodia tracking
ECFCs expressing Lifeact-tdTomato were trypsinized and maintained in suspension in full media in the presence of either DMSO, Actinomycin D, or puromycin (Takara Bio) for 5 minutes then reattached to collagen-coated coverslips.Alternatively, ECFCs treated with siRNA for 24 hours were then transfected with Lifeact-tdTomato for 24 hours then reattached to collagen-coated coverslips.Lifeact-expressing ECFCs were imaged using a 63x objective (NA:1.2) in 10s intervals for 1 hour on a Zeiss axio observer.Lifeact signal-to-noise ratio was improved by two sequential rounds of thresholding and variance filtering (in 10 and 4 pixel interrogation windows).Videos were processed in matlab to track lamellipodial dynamics 5 µm from the edge of the cell.Background was subtracted and lamellipodia displacements (d) tracked using the built-in Farneback optic flow function in matlab.The cell centroid at each time was used as an internal reference point for constructing a position matrix (p).The dot product (d•p) of the position matrix (p) and the displacement matrix was then used to classify lamellipodia as protrusive (d•p > 0) or retractive (d•p < 0).Average lamellipodia protrusion or retraction velocity were calculated across 20 minutes of cell spreading.Experiments were repeated in triplicate where at least 1-4 cells were images per condition.

Microscopy and image analysis
Epifluorescence images of fixed and live samples were taken on Zeiss Axio Observer equipped with a monochromatic Axiocam 702 (Zeiss) at 23° C using 5x (NA: 0.16), 10x (NA: 0.3), 20x (NA: 0.8), 40x (NA: 0.6), and 63x (NA: 1.2) Zeiss objectives.Data acquisition was done using the ZEN imaging suite (Zeiss).Confocal image stacks were taken on a laser-scanning Leica DMI 6000 inverted microscope using a 63x (NA: 1.4) Leica objective.Deconvolution of z-stacks were performed using Hugyens professional image processing application with a theoretical point spread function.
Focal adhesion morphology was measured after rolling ball subtraction with a radius of 50 pixels and thresholding the Huang thresholding method.Adhesions were identified on a per cell basis, across three experiments as being greater than .1 µm 2 and less than 10 µm 2 .All image analysis (morphometrics, image adjustments, and individual cell tracking) were performed using an open access NIH software platform, FIJI (Schindelin Nat Methods 2012).
Z-stacks of pMLC, actin, and DAPI signal were taken through the thickness of a cell at a depth of 0.75 µm.Images underwent deconvolution using Hugyens professional described above and height was measured using the z-axis intensity profile tool from FIJI.DAPI and pMLC height were estimated as the position above the nuclei with less than 10 % of the max fluorescent intensity of a given stack, where the region with the highest fluorescent intensity was typically in line with the middle of the nucleus.Representative line plots show relative fluorescent intensity, normalized across a given image stack.Fluorescent signal less than 10 % of the max fluorescent intensity was assumed to be background and made zero.
Quantitative reverse transcription PCR (RT-qPCR) Total RNA was isolated and purified using the RNeasy mini kit (Qiagen).0.5 µg of total RNA was reversed transcribed (Applied Biosystems) using the manufacturer's instructions in a Mastercycler nexus gradient (Eppendorf).cDNA was mixed with SYBR green mastermix (Applied Biosystems) and 0.4 µM forward and reverse primers (CTGF: TTAAGAAGGGCAAAAAGTGC (Forward 5'-3') and CATACTCCACAGAATTTAGCTC (Reverse 5'-3'); Sigma Aldrich) in wells of 96 well PCR plate (Applied Biosystems).cDNA was amplified and quantified using a Step-one Plus real-time PCR system (Applied Biosystems).Relative gene expression was quantified using the ∆∆CT method normalizing target gene expression with the housekeeping gene 18S.

RNA-sequencing & Gene Ontology analysis
Total RNA was assessed for quantity and quality with RNA ScreenTape assay of Agilent 2200 TapeStation System (Agilent Technologies).Average RIN values were 10. 200 ng of RNA was used as input for library prep.Libraries were prepared using TruSeq Stranded mRNA HT Sample Prep Kit (Illumina) as per standard protocol in the kit's sample prep guide.Libraries were assayed for overall quality using D1000 ScreenTape assay of Agilent 2200 TapeStation System (Agilent Technologies).Average library molarity from Tapestation was 103.331+/-24.6nM.Samples were multiplexed for sequencing.A MiSeq test lane was run to check pool balance followed by 100bp single-read sequencing on an Illumina HiSeq 4000 sequencer.Illumina's bcl2fastq2 version v2.20.0.422 software was used for demultiplexing and converting bcl to fastq files.
Reads were filtered to retain only high-quality reads.In addition, ribosomal reads and repeats were eliminated by alignment to a generic set of ribosomal/repeat sequences.Remaining reads were processed with RNA-Seq Unified Mapper (RUM) (Grant et al., 2011).Following RUM alignments to human genome (build hg19) and known transcripts (RefSeq, UCSC, ENSEMBL), a feature-level quantitation (transcript, exon, and intron) were output as result.To analyze global gene expression profiles, the number of uniquely aligning reads to mRNA transcripts in RefSeq were extracted from the RUM output.
Statistically significant differences and fold changes in individual genes were computed via DESEq2 v. 1.36.0.DEGs were defined as having a Wald test FDRadjusted P value less than 0.05.GO was conducted separately for genes significantly upregulated or downregulated in siYAP/TAZ samples compared to siControl using PANTHER overrepresentation test on PANTHER 17.0, available on the GO Resource.We used the GO biological process complete annotation set with Fisher's exact test and FDR-corrected P values.The test was applied using a query set of genes significantly elevated or decreased (Wald test FDR-adjusted P value less than 0.05; n= background of all nonzero features identified in the bulk transcriptome.
Fluorescent images were collected utilizing a Nikon Ti2 and Yokogawa CSU-W1 spinning disk confocal microscope between 28-48 hpf at 10x magnification.To maintain timing consistency between treatments, embryos were fixed in 4% PFA for 24 hours at 4 o C before being embedded in 0.8% low melting point agarose for imaging.Z-stacks were acquired using a 1um step size, and max projections generated using FIJI software.ISV length was measured using FIJI and defined as the length of the ISV bounded between the dorsal aorta and the DLAV (dorsal longitudinal anastomotic vessel).Five ISV's were measured per embryo and an average per embryo is represented in the graphs.The data is normalized to the height of the dorsal ventral axis of the animal it was collected from to account for differences is mounting and imaging.

Statistics
Statistical analyses of cell migration and ISV morphogenesis were performed on Graphpad Prism 6 statistical analysis package.Data are presented with data points were possible and mean ± standard error or standard deviation; figure captions describe data presentation.All experiments were performed at least in triplicate.Multiple comparisons were made using analysis of variance (ANOVA) with Tukey or Sidak's post hoc test for pairwise comparisons of normally distributed homoscedastic data.Data were considered to fit the ANOVA assumptions if the residuals of a data set were normally distributed.Data that did not meet the ANOVA criteria were analyzed by Kruskal-Wallis with Dunn's post hoc test.Comparisons between two data-sets were made using Student's unpaired two-tailed t-test, or non-parametric by Mann-Whitney, when necessary.

Figure 1 :
Figure 1: Cytoskeletal-transcriptional feedback prevents excessive focal adhesion maturation to enable persistent motility in human ECFCs.A. Model of YAP/TAZ-mediated transcriptional feedback regulation of the cytoskeleton, as established in (Mason et al., 2019; van der Stoel et al., 2020).B-D.ECFCs were treated with the transcription inhibitor, Actinomycin D, and stained for vinculin, phalloidin (F-actin), and DAPI (nuclei) after 24 hours.E-G.ECFCs were depleted of both YAP and TAZ by RNAi and stained for vinculin, phalloidin (F-actin), and DAPI (nuclei) after 28 hours.Both transcription inhibition and YAP/TAZ depletion induced motile arrest, indicated by loss of fan-shaped lamellipodial migration and acquisition of an oval shape, ringed by mature focal adhesions (B, E).Vinculin+ focal adhesions (C, F) and mature vinculin+ focal adhesions, defined as greater than 1 µm (D, G).N = 20-22 cells per condition.**** p < 0.0001, Student's two-tailed unpaired t-test.Data are shown as mean ± S.E.M.

Figure 3 :
Figure 3: De novo gene expression and RhoA-YAP/TAZ signaling regulate zebrafish ISV morphogenesis in vivo. A. Bright field and spinning disc confocal imaging of the (Tg(fli:egfp) y1 ) zebrafish transgenic line at 29, 32, and 35 hpf outlining the timeline of ISV morphogenesis in the zebrafish trunk.B. Schematic of experimental design.C. Representative spinning disc confocal images of DMSO, Actinomycin D (Act.D) 10 µg/ml (low), Act.D 25 µg/ml (high) and Puromycin treated zebrafish embryos.Compounds were added at 29 hpf and ISV morphogenesis imaged at 32 hpf (top) and 35 hpf (bottom).D,E.Quantification of ISV length, tracking individual fish over time, with data plotted as mean ± S.D. (D), vs comparative analysis at isolated time points with data plotted as mean ± S.E.M. * p < 0.05, two-way ANOVA with Sidak's post hoc test (E).DMSO n = 8; Act.D (low) n = 8; Act.D (high) n = 8; and Puro n = 6.F. Representative spinning disc confocal images of DMSO, Rockout, Blebbistatin (Blebb), and Verteporfin

Figure 5 :
Figure 5: YAP/TAZ-depletion in human ECFCs induces motility and morphodynamic inflection at 4 hours post-attachment, followed by progressive motility arrest, on both soft and stiff substrates.A. Experiment schematic: mTomato-expressing ECFCs were treated with siRNA targeting YAP and TAZ or non-targeting control siRNA for 24 hours prior to seeding on 1, 8, or 18 kPa MeHA matrices.Cell morphology and migration were tracked for 24 hours post-attachment or fixed for immunofluorescent imaging at 4 hours post-attachment.B-D.Instantaneous migration speed, cell area, and circularity were

Figure 6 :
Figure 6: Rescue of vascular morphogenesis by inhibitor washout depends on washout timing.A. Schematic diagram of continuous inhibitor treatment.Zebrafish embryos were treated with transcription inhibitor (Act.D), translation inhibitor (Puromycin), or YAP/TAZ inhibitor (Verteporfin), at 29 hpf and imaged at 32 hpf and 35 hpf.B. Schematic diagram of three-hour washout, in which inhibitors were washed out after 3 hours, prior to the time to feedback loop closure at tic = 4 hrs.Zebrafish embryos were treated with the indicated compounds at 29 hpf, had the inhibitors removed at 32 hpf, and were imaged as the 'washout' condition at 35 hpf.C. Representative spinning disc confocal images of DMSO, Act.D 25 ug/ml, Puromycin, and Verteporfin treated zebrafish embryos at 32 hpf (top), 35 hpf (middle), or at 35 hpf after 3 hours of incubation in the compounds followed by 3 hours of compound wash out (bottom).D,E,F.Quantification of ISV length, tracking individual fish over time.Plots indicate ISV growth rates under conditions of embryos continuously incubated in the indicated compound versus the associated washout condition.Act.D versus washout is shown in (D), Puromycin (Puro) versus washout is shown in (E), and Verteporfin (VP) versus

Figure 7 :
Figure 7: Feedback history and cytoskeletal state prior to re-adhesion alter cytoskeletal and adhesion remodeling. A. Transcription inhibition history experiment schematic: mTomato-expressing ECFCs were either pretreated with actinomycin D before plating (-20 hours) or at the time of attachment.

Figure 8 :
Figure 8: Transcriptional feedback suppresses myosin light chain phosphorylation, limits cell spreading through supranuclear membrane constriction, and promotes actin cap stress fiber maturation.A-C.ECFCs were treated for 24 hours with DMSO or Actinomycin D and lysed for relative

Movie 7 :
Representative cytoskeletal dynamics of control and Act.D-treated LifeAct-tdTomadtoexpressing ECFCs during first hour after adhesion to glass cover slip.Images were taken in 10second intervals for 1 hour after attachment.Time is shown as minutes:seconds.Movie 8: Representative cytoskeletal dynamics of control and YAP/TAZ depleted LifeAct-tdTomadto-expressing ECFCs during first hour after adhesion to glass cover slip.Images were taken in 10-second intervals for 1 hour after attachment.Time is shown as minutes:seconds.

2 :
Effects of YAP/TAZ depletion on ECFCs motility and morphodynamics after attachment to soft (1 kPa), moderate (8 kPa), and stiff (18 kPa) matrices.(A-B) ECFCs were depleted of YAP and TAZ for 24 or 48 hours then lysed for relative protein quantification by Western Blot.