Myosin IIA-mediated forces regulate multicellular integrity during vascular sprouting.

Angiogenic sprouting is a critical process involved in vascular network formation within tissues. During sprouting, tip cells and ensuing stalk cells migrate collectively into the extracellular matrix while preserving cell-cell junctions, forming patent structures that support blood flow. Although several signaling pathways have been identified as controlling sprouting, it remains unclear to what extent this process is mechanoregulated. To address this question, we investigated the role of cellular contractility in sprout morphogenesis, using a biomimetic model of angiogenesis. Three-dimensional maps of mechanical deformations generated by sprouts revealed that mainly leader cells, not stalk cells, exert contractile forces on the surrounding matrix. Surprisingly, inhibiting cellular contractility with blebbistatin did not affect the extent of cellular invasion but resulted in cell-cell dissociation primarily between tip and stalk cells. Closer examination of cell-cell junctions revealed that blebbistatin impaired adherens-junction organization, particularly between tip and stalk cells. Using CRISPR/Cas9-mediated gene editing, we further identified NMIIA as the major isoform responsible for regulating multicellularity and cell contractility during sprouting. Together, these studies reveal a critical role for NMIIA-mediated contractile forces in maintaining multicellularity during sprouting and highlight the central role of forces in regulating cell-cell adhesions during collective motility.

INTRODUCTION

Collective migration is a process in which cohorts of cells move in a coordinated manner so that cell–cell contacts are maintained. This is a highly regulated process that is critical in a variety of pathological and developmental morphogenic events such as tumor cell invasion and sprouting angiogenesis (Friedl and Gilmour, 2009; Scarpa and Mayor, 2016). This cooperative movement of cells is of particular importance during sprouting morphogenesis because the endothelial cells that make up new vessels must form patent, nonleaky structures capable of supporting blood flow. The three-dimensional (3D), multicellular structures that form during sprouting are classically composed of leader endothelial cells, or tip cells, that proteolytically degrade the surrounding extracellular matrix to migrate toward sources of angiogenic factors. As they advance, tip cells retain physical cell–cell adhesions with a trailing cohort of stalk cells that form the lumenized trunk of the sprout (Adams and Alitalo, 2007; Carmeliet ). Although the molecular drivers of vascular sprouting and tip and stalk cell specification, via the VEGFR2 and Notch1-DLL4 pathways for example, have been extensively described, much less is known about how cellular mechanics and force regulation influence sprout morphogenesis (Gerhardt ; Hellström ; Benedito .

The important role of cell-mediated forces in multicellular migration and morphogenesis has been highlighted in various cell types and through a variety of in vitro approaches. For example, in epithelial cells, two-dimensional (2D) assays have been utilized to study the impact of external factors such as substrate stiffness and the role of actomyosin-based contractility in inducing leader cell formation and maintaining coordinated movements of cell cohorts during planar migration (Ng ; Rausch ). Similarly, in 3D and in vivo assays, studies have demonstrated that cell contractility is needed to enable multicellular 3D invasion and proper morphogenesis of epithelial ductal structures (Ewald ; Gjorevski ). In endothelial cells, in particular, myosin-mediated cell contractility has been investigated in the context of multicellular organization. For instance, studies using natural and fibrillar matrices have demonstrated that actomyosin-based cell contractility is necessary to support the formation of multicellular networks through processes mimicking vasculogenesis (Lyle ; Davidson ). In 3D settings, endothelial cells have been shown to depend on actomyosin-based contractile forces to invade their surrounding matrix and maintain sprout structures after invasion (Kniazeva and Putnam, 2009). Intriguingly, several works have begun to quantify the ability of endothelial sprouts to apply contractile forces to their surroundings. For instance, Kniazeva demonstrate that vascular sprout formation is directly correlated with enhanced contractility and the rate at which endothelial cells are able to deform the matrix. Despite this growing body of evidence to suggest the importance of actomyosin-based contractility in sprout morphogenesis, several key questions regarding the spatial organization and magnitude of deformations, as well as a more detailed evaluation of the role of specific myosin isoforms in multicellular sprout structure, still remain.

In this work, we aimed to fill this gap by investigating the role of NMII-mediated cell contractility during angiogenesis, given that many of the signaling pathways driving essential cytoskeletal functions in migrating cells converge on the activation of nonmuscle myosin II (NMII). We first established a general understanding of how endothelial cells within multicellular sprouts generate contractile forces in a microfluidic, biomimetic model of angiogenesis (Nguyen ). We then perturbed NMII-generated contractility using a NMII-specific inhibitor, blebbistatin, and CRISPR-Cas9–mediated gene editing to investigate the functional significance of these forces in sprout morphogenesis.

RESULTS

Angiogenic sprouts produce spatially organized matrix deformations

To investigate the importance of cell contractility in sprouting angiogenesis, we first sought to determine if and where forces are generated in multicellular sprout structures in a 3D biomimetic model of angiogenesis previously developed in our lab (Nguyen ). Briefly, this microfluidic device consists of two hollow cylindrical channels embedded in a central collagen hydrogel compartment. In one channel, endothelial cells are seeded and form a uniform monolayer, and in the opposing channel, a cocktail of angiogenic factors is supplied to trigger sprouting into the collagen hydrogel. To track the origins of cell-generated forces within sprouts, we adapted this platform so that fiduciary fluorescent beads were embedded into the collagen matrix, enabling measurements of cell-generated matrix deformations (Supplemental Figure S1, A and B). Using methods established for 3D traction force microscopy (TFM; Legant ), we measured the magnitude of bead displacements generated by the tip and stalk cells of sprouts that formed 3 to 4 d after the initiation of cocktail exposure and then used these displacements to calculate maps of 3D mechanical deformations in the surrounding ECM. Because collagen matrices exhibit significant nonlinear elastic behavior in response to forces (Münster ), one cannot calculate absolute stresses (normalized forces) accurately from these deformations, as would typically be done in 3D TFM protocols. However, the deformation maps nonetheless can be used to define the relative distribution of where cellular forces originate. In these sprouts, we observed significant deformations generated at the leading front of the sprout, particularly by tip cell protrusions (Figure 1A). These displacements are of the same order of magnitude as those generated by other cell types in a similar collagen matrix (Koch ; Thievessen ). Along the stalk of the sprout, the majority of the measured displacements were minimal in magnitude, suggesting that the tip cell is likely the predominant force-generating unit within a sprout. Closer examination of the displacements revealed that these displacements pointed inward toward the sprout, indicating that forces generated by the tip cell are indeed contractile, as opposed to outward-pushing forces (Figure 1A′).

FIGURE 1:

Angiogenic sprouts generate contractile forces necessary to maintain multicellular migration. (A) A representative map of matrix deformations (colored arrows) generated by an endothelial sprout (gray). (A’) A zoom-in view of the front region of the invading tip cell. (B) Representative F-actin (phalloidin, green) and nuclei (DAPI, blue) images of sprouts at day 5 after sprouting under 30 μM blebbistatin or DMSO control treatment. Scale bar is 200 μm. (C) Quantification of invasion density as measured by the number of invaded nuclei. (D) Quantification of invasion distance as measured by distance from edge of main channel to end of tip cell invasive front. (E) Quantification of cell–cell breakages normalized to total number of invaded cells. (F) Representative F-actin (phalloidin, green) and nuclei (DAPI, blue) images of indicated conditions at day 4 after 24-h blebbistatin treatment. Scale bar is 50 μm. (G) Quantification of detached tip cells (left plot) and stalk cells (right plot) following 24-h blebbistatin treatment. NDMSO = 3 devices, 4 fields of view per device, NBlebbi = 5 devices, 4 fields of view per device. Bars represent mean ± SEM. * represents p < 0.05 as determined by a t test (C, D, E) or one-way ANOVA with post hoc Tukey HSD test (G).

Inhibition of contractility alters multicellularity of sprout morphogenesis in vitro and in vivo

Given this understanding of the forces generated by angiogenic sprouts, as well as the global importance of actomyosin contractility in producing cellular traction forces, we sought to investigate the morphogenic significance of cell-generated forces during sprouting. To this end, a gradient of angiogenic factors was applied within our in vitro sprouting model in the presence or absence of blebbistatin, a small-molecule inhibitor of NMII that has been shown to reduce cellular traction forces significantly (Bhadriraju ; Liu ). We exposed the endothelialized channel simultaneously to blebbistatin and angiogenic factors to study how inhibition of NMII may affect various aspects of sprout morphogenesis from the onset. After 5 d of exposure, sprouts were fixed and stained for analysis (Figure 1B). While inhibition of cell contractility had no appreciable effect on the cellular density and speed of invasion from the parent vessel through the collagen matrix (Figure 1, C and D), blebbistatin treatment dramatically increased the frequency of cell–cell breakages in comparison with sprouting achieved in the presence of the vehicle control (Figure 1E).

To parse out whether this breakdown in multicellularity can specifically be attributed to contractility inhibition during sprout elongation, as opposed to other stages of sprout development such as initiation of invasion, we applied varying concentrations of blebbistatin (10, 30, or 50 μM) for 24 h to intact multicellular sprouts that had already extended for 3 d with angiogenic factor stimulation alone. Samples were subsequently fixed, stained for analysis, and evaluated for sprout morphology through high-magnification imaging. While sprouts in the control treatment group formed highly multicellular structures with few sporadic detached cells within the invasion region, treatment with blebbistatin in a dose-dependent manner resulted in an increasing loss of cell–cell contacts between cells (Figure 1F). Given the concentration of force generation in tip cells observed previously, we further categorized the cell–cell breakages as either tip–stalk or stalk–stalk, depending on whether the breakage occurred between a leader cell and its following stalk cell or between two nonleader cells. At the highest concentration of blebbistatin, we saw a significant increase in the number of tip–stalk breakages in comparison with that in the control group. In contrast, this dose of blebbistatin did not induce a significant increase in the incidence of stalk–stalk breakages (Figure 1G). Given the potential off-target effects of a pharmacologic inhibitor such as blebbistatin, we also examined the impact of Y-27632, a selective Rho-associated kinase (ROCK) inhibitor, on sprouting. Indeed, ROCK inhibition led to a similar trend in the loss of tip–stalk contacts during sprout extension (Supplemental Figure S2, A and B).

To establish the relevance of these findings in an in vivo system, we used a mouse retina model in which the vascular network develops across the neuroepithelium during the first week of life. Mouse pups at postnatal day 5 (P5) were injected with blebbistatin and killed 6 h after treatment. Their retinas were subsequently dissected, fixed, and stained for the endothelial cells, using isolectin B4. Corroborating our in vitro findings, a significant number of detached endothelial cells appeared at the migrating front of the developing vascular network (Supplemental Figure S3, A and B). Together, these in vitro and in vivo findings suggest that cytoskeletal forces are not necessarily required for enabling 3D chemotaxis of endothelial cells toward gradients of angiogenic factors, but are required for the maintenance of their multicellularity during sprout extension.

Cell–cell adhesions are disrupted upon contractility inhibition

Vascular endothelial (VE) cadherin is the dominant cadherin involved in adherens junctions (AJs) of ECs and has been shown to regulate vascular permeability, dynamic changes in cell–cell adhesions, and transcellular signaling (Lampugnani ; Mehta and Malik, 2006). Cadherins within endothelial AJs have also been shown to link intracellularly with the cytoskeleton by the dynamic, phosphorylation-mediated binding of the intracellular cadherin domain to the actin cytoskeleton via protein complexes consisting of molecular partners such as the catenins (Aberle ; Jou ). Given the importance of VE cadherin in maintaining cell–cell junctions, as well as the molecular connection between VE cadherin and the actin cytoskeleton, we sought to determine whether alterations in the quality of cell–cell junctions occur in response to NMII inhibition.

To this end, day 3 sprouts were treated with blebbistatin for 24 h and subsequently fixed and immunostained for VE cadherin. High-magnification confocal images of these samples were captured, segmented into regions of tip–stalk and stalk–stalk adhesion, and characterized as either discontinuous or continuous (Figure 2, A and B). Under control conditions, we found that most junctions appeared to have a continuous morphology (Figure 2C), which is characteristic of a more stable phenotype (Bentley ). With increasing blebbistatin concentration, tip–stalk adhesions became difficult to observe due to the increase in cell–cell breakages. However, quantification of the tip–stalk adhesions that survived blebbistatin treatment appeared heavily punctate and discontinuous in profile, suggesting that the junctions were unstable (Figure 2C). In fact, under the conditions with the strongest inhibition of contractility, almost all the observed tip–stalk adhesions were classified as discontinuous. Similarly, stalk–stalk cell adhesions shifted to an increasingly unstable adhesion profile in a dose-dependent manner (Figure 2D). At the highest level of blebbistatin treatment, nearly 50% of the stalk–stalk junctions were categorized as discontinuous. Together, these results demonstrate that broad inhibition of contractile activity also causes disorganization of cell–cell junctions, which may in part contribute to the increase in cell–cell breakages observed with blebbistatin treatment.

FIGURE 2:

Inhibition of myosin-mediated contractility alters the organization of cell–cell adhesions. (A) Representative confocal immunofluorescence images of tip–stalk (top panel) and stalk–stalk (bottom panel) adhesions from indicated conditions at day 4 after 24-h blebbistatin treatment. VE cadherin is shown in white; nuclei are shown in blue. Scale bars are 20 μm. (B) Representative images of catalogued adhesion types. Two neighboring cells are false-colored (blue, green). White arrows show adhesions of interest. Scale bar is 25 μm. (C) Quantification of adhesion profiles of tip–stalk adhesions across treatment groups. (D) Quantification of adhesion profiles of stalk–stalk adhesions across treatment groups. N = 6 sprouts per condition from three independent experiments. Bars represent mean ± SEM. * represents p < 0.05, ** represents p < 0.01, and *** represents p < 0.001 as determined by two-way ANOVA with post hoc Dunnett test.

NMIIA is predominantly responsible for maintaining multicellularity during sprouting

The observations connecting cell contractility and cell–cell adhesiveness during sprouting have thus far been based upon pharmacologic inhibition of NMII. However, there is a growing body of evidence highlighting differential roles of distinct nonmuscle myosin II isoforms in various cell types and functions, including regulation of cell–cell junctions (Shih and Yamada, 2010; Jorrisch ; Liu ; Gutzman ). Motivated by these studies and to validate our blebbistatin findings, we sought to investigate the individual contributions of NMII isoforms A and B (henceforth, NMIIA and NMIIB) to mediating the multicellularity of angiogenic sprouting achieved within our system. Using CRISPR-Cas9–mediated gene editing, we generated MYH9 and MYH10 knockout endothelial cell lines (hereafter IIAKO and IIBKO, respectively; Figure 3A). While immunofluorescence staining and Western blotting for NMIIA in IIAKO endothelial cells and NMIIB in IIBKO endothelial cells confirmed loss of isoform-specific myosin content, a faint NMIIA signal observed in both immunofluorescence and Western blot analysis suggests a partial (mixed homozygous and heterozygous) knockout (Figure 3, A and B). Nonetheless, IIAKO endothelial cells strikingly appeared to adopt an atypical cell morphology, quantified by the loss of lamellipodia structures (Supplemental Figure S4, A and B), which have been demonstrated to contribute to various endothelial functions including cell migration and the ability to form proper permeability barriers (Breslin ; Cao ). This morphology resembled that of blebbistatin-treated control cells (Figure 3C). In contrast, IIBKO endothelial cells maintained basic physical features similar to those of the control cells, namely retaining the broad lamellipodia characteristic of normal endothelial cells.

FIGURE 3:

Specific loss of the NMII isoform A results in a reduction of sprout multicellularity. (A) Representative Western blots demonstrating diminished protein levels of NMIIA and NMIIB in endothelial cells. (B) Representative immunofluorescence images of endothelial cells stained for NMIIA and NMIIB. Top panels represent NMIIA or NMIIB (white), middle panels represent F-actin (phalloidin, white), and bottom panels are merged images of NMIIA/B (red), F-actin (phalloidin, green), and nuclei (DAPI, blue). Insets show zoomed-in area of corresponding image. Scale bars are 30 μm. (C) Representative bright-field images of endothelial cells cultured in 2D. Scale bar is 50 μm. (D) Representative immunofluorescence images of sprouts stained with CD31 (green) and nuclei (DAPI, blue). Scale bar is 100 μm. (E) Quantification of the number of invaded cells as determined by the number of invaded nuclei. (F) Quantification of invasion distance of sprouts on day 3 after sprouting initiation. Each dot represents one sprout. (G) Quantification of cell–cell breakages normalized by the total number of invaded cells. N = 715 sprouts from 16 devices (Scr); 816 sprouts from 17 devices (IIAKO); and 789 sprouts from 17 devices (IIBKO). Bars represent mean ± SEM. ** represents p < 0.05 as determined by one-way ANOVA, E and G, or Kruskal–Wallis test, F.

To determine the functional role of each myosin II isoform in sprout morphogenesis, we introduced the IIAKO and IIBKO endothelial cells into our in vitro sprouting model. After 3 d of stimulation with a cocktail of angiogenic factors, the devices were fixed and analyzed for sprout morphology through immunofluorescence staining for CD31 (Figure 3D). Overall, perturbations to NMIIA and NMIIB had no significant effect on the overall ability of cells to invade, as measured by the number of cell nuclei found within the collagen matrix (Figure 3E). Also, whereas loss of NMIIB resulted in no significant changes to sprout length, loss of NMIIA led to a slight but significant reduction in sprout length (Figure 3F). Most importantly, sprouting from IIAKO endothelial cells was characterized by an increase in the portion of cell–cell breakages, whereas loss of NMIIB led to no significant change in this regard (Figure 3G). Taken together, these findings confirm that NMII-mediated contractility is critical for maintaining multicellularity of angiogenic sprouts and highlight the distinct roles of NMIIA and NMIIB in sprouting, with the NMIIA isoform predominantly dictating the morphogenesis of multicellular structures.

NMIIA is responsible for generating contractile forces in sprouts

Based on these observed differences in the contributions of NMIIA and NMIIB to the multicellular morphogenesis of angiogenic sprouts, we hypothesized that the forces generated by sprouts observed earlier (Figure 1) would largely be dependent on NMIIA. To test this hypothesis, we measured the displacement of fiduciary beads generated by sprouts formed within the in vitro angiogenesis model. Indeed, we found that multicellular sprouts formed by IIAKO endothelial cells produce minimal deformations on the surrounding matrix on the average (Figure 4A). This loss of large bead displacements in the matrix surrounding the leading front of the IIAKO sprouts led to a shift toward smaller deformations in the overall frequency distribution of bead displacements as compared to the Scramble control sprouts (Figure 4, B and C). In line with our previous findings suggesting a more limited role of NMIIB in sprouting, loss of NMIIB did not significantly alter the ability of multicellular sprouts to generate matrix deformations and cellular tractions. This reduction in the ability of NMIIA-deficient sprouts to generate contractile forces, which in turn produce the observed matrix deformations, was also confirmed through 2D traction stress measurements of single cells (Supplemental Figure S5, A and B), demonstrating that NMIIA is the main contractility-generating unit in endothelial cells in both 2D and 3D settings. Overall, these data indicated that generation of cell contractility, through the NMIIA isoform, is essential for maintaining multicellularity during 3D chemotactic migration in the context of sprout morphogenesis.

FIGURE 4:

NMIIA is the predominant force-generating unit within multicellular sprout structures. (A) Example matrix deformation maps (colored arrows) generated by endothelial sprouts (gray). (B) Each dot represents the magnitude of one displacement vector. The scatterplot is a compilation of N = 4 (Scr), 5 (IIAKO), and 4 (IIBKO) sprouts, one sprout per device (= experiment). Only the vectors within 5 μm of the sprout surface were considered for analysis, and displacements over 5 μm are not displayed on the graph. (C) The average bead displacement for each sprout in each condition is displayed. Bars represent mean ± SEM. **** represents p < 0.0001 as determined by a Kruskal–Wallis test. Statistical analysis of average displacements was completed via one-way ANOVA with a post hoc Dunnett test.

DISCUSSION

Extensive research investigating the functions of NMII in various cell types and contexts has demonstrated its role as a universal regulator of cell contractility and many cytoskeletal functions. Of particular relevance, myosin has been identified as a key player in regulating migration and cell–cell adhesions. Extensive studies examining the dynamics of cell movements on 2D substrates have demonstrated that loss of myosin activity can lead to reduced migration in most contexts, as detachment of the posterior anchors of migrating cells appears to be hampered without myosin activity (Eddy ; Kolega 2003). This migrational slowing has also been confirmed in some studies of cells in 3D contexts (Doyle ; Fischer ; Petrie ). However, this finding is not universal. For example, studies have reported enhanced migration speeds as large stable focal adhesions within cells shift to smaller nascent adhesions (Gupton and Waterman-Storer, 2006; Kim and Wirtz, 2013). Our studies, which demonstrate that blocking cell contractility does not hinder the overall invasiveness of endothelial cells during sprout morphogenesis, build upon this body of evidence suggesting a highly context-dependent role of NMII in cell migration. For instance, in our platform, endothelial invasion is guided by a chemotactic gradient of angiogenic factors, whereas most studies of myosin inhibition have been in the context of random cell motility. Perhaps in the context of 3D invasion, which is inherently slower than 2D migration because of the requirement for MMP-mediated matrix degradation, migrational persistence is the key determinant of invasion extent. Indeed, myosin has been reported to be required for persistence of a migration axis (Lo ; Fischer ). In our model system, the presence of a chemotactic gradient may have compensated for the loss of this particular function of NMII. More systematic studies to investigate the role of NMII in chemotaxis-driven invasion of endothelial cells would help to confirm this possibility and improve our basic understanding of the coordinated contribution of NMII, chemotactic gradients, 3D migration, and proteolysis-dependent invasion.

Instead of hindering migration, we observed that reduction in NMII activity led to the fragmentation of multicellular sprouts. This metric of collectiveness, a critically important aspect of neovessel formation, is not often reported in classic in vitro angiogenesis studies but has more recently been incorporated (Shin ; Wimmer ; Del Amo ). Here, we observed that NMII-­mediated contractility impacts cell cohesiveness directly and we specifically found that blebbistatin-induced breakages preferentially occur between tip and stalk cells of in vitro angiogenic sprouts. This is likely in part due to the fact that, with blebbistatin treatment, a greater portion of tip–stalk adhesions take on the discontinuous-junction phenotype, which is more prone to disassembly. A similar phenomenon was observed in vivo in the mouse retina model of angiogenesis, where treatment with blebbistatin led to a significant increase in single-cell breakages located at the leading front of the developing vascular plexus. Nonetheless, a more systematic genetic study of NMII loss in the in vivo setting may be needed to validate the role of NMII in vivo fully, given that blebbistatin may target various other cell types that contribute to angiogenesis in vivo (e.g., macrophages and smooth muscle cells).

However, in support of the notion that myosin inhibition selectively induced tip–stalk cell breakages due to the already compromised integrity of their adhesions, a recent study has linked disorganized junction morphologies to adhesions between migratory cells at the migrating front in a scratch wound assay, whereas linear adhesions are observed in less dynamic settings, such as between cells in a quiescent monolayer (Cao ). Accordingly, treatment of collectively migrating endothelial cells with VEGF, a well-­established disrupter of endothelial adherens junctions, enhances the number of cells that detach from the main cohort (Huveneers ; Oubaha ). Along the same lines, recent work published by Hayer shows that confluent endothelial monolayers treated with the ROCK inhibitor Y-27632, which has been shown to reduce adherens-junction formation (Abraham ), displayed severely diminished coordinated movements. More generally, extensive studies on junctional dynamics, primarily in the context of E-cadherin–mediated adhesions in epithelia, have demonstrated that myosin is critical in mediating junction formation (Shewan ; Smutny ). Interestingly, recent studies have suggested that also myosin’s actin-crosslinking functions, as opposed to its force-generating function per se, are critical to its role in epithelial junction assembly (Ivanov ). In the context of endothelial junctions, we have previously found that both force (via myosin) and actin cross-linking (via Rac1) contribute to junction assembly (Liu ). Thus, while we were not able to distinguish their roles in this study, it is likely that both contribute to the overall process of junction assembly and stabilization.

Interestingly, in the context of endothelial cell–cell junction assembly dynamics in regulation endothelial barrier function, contractility appears to have a biphasic effect. We and others have previously reported increased cell–cell junction assembly with moderate increases in myosin-mediated contractility in endothelial monolayers (Liu ), while others have shown that very high contractility has been shown to abruptly disrupt junctions (Birukova ; Valent ; Yeh ). Bordeleau demonstrated that increasing substrate rigidity, which has been linked to enhanced cell contractility (Califano and Reinhart-King 2010), appears to disrupt EC barrier function. Taken together, these works investigating the impact of reducing and enhancing cell contractility suggest that there is likely an intermediate level of cell contractility that is needed to maintain the integrity of cell–cell adhesions. Thus, we believe that under physiologic sprouting conditions recapitulated within our in vitro sprouting model, an intermediate level of cell contractility aids in preserving sprout multicellularity, whereas reducing contractility through small-molecule inhibitors or genetic manipulations to NMII leads to cell breakages. Thus, we hypothesize that perturbations to enhance cell contractility under sprouting conditions may result in a similar reduction of sprout multicellularity and should be explored further in future studies.

In addition to identifying a role for myosin in maintaining the multicellularity of invasive sprouts, in this study we also were able to demonstrate that NMII is not used equally in all cells. Using methods from 3D traction force microscopy, we found that sprout morphogenesis occurs with contractile forces mainly exerted by the leader tip cells. This observation is in line with the notion that leader tip cells actively pull on the surrounding matrix while the following stalk cells likely interact more passively with the matrix during sprout elongation (Korff and Augustin, 1999; Du ). Pushing forces have been reported in the context of macrophage podosomes and other forms of leukocyte migration (Labernadie ), but pulling forces have been the predominant arrangement described in migrating cells (Gardel ; Gjorevski . The existence of pulling forces implies that cells proteolytically degrade the ECM to create a path for migration rather than pushing matrix fibers away. Interestingly, stimulation of endothelial cells with VEGF, which primarily activates invasive signaling only in leader tip cells, triggers increased contractile forces and stress fiber formation (Yang ; LaValley ). The recognition that this NMII activity is focused within the tips of the sprouts gives greater appreciation that the mechanochemical signaling pathway described here—linking force to junction stability—is used selectively, not universally, and underscores the continued need to develop methods for mapping when and where forces are invoked.

Despite the multifaceted approaches employed in these studies to investigate and describe the role of NMII in sprout morphogenesis, several limitations on the various experimental approaches exist. For example, the entirety of the in vitro studies presented in this work were completed using HUVECs as the model cell type. Though HUVECs have been utilized extensively in a broad range of in vitro angiogenesis studies, recent work highlighting differential functions of ECs from specific organs or levels of the vascular hierarchy (Rafii ; Hilfenhaus ) opens the possibility that the role of NMII in sprout morphogenesis could vary across EC subtypes. Additionally, the in vitro angiogenesis studies carried out in this work have been done utilizing acellular collagen matrices in the absence of supportive cell types such as pericytes and other tissue-specific mural cells; both matrix composition and these ancillary cells have been demonstrated to impact vascular network formation, sprouting, and barrier function (Ribatti ; Rao ; Kim ; Alimperti ). Though secondary support cells such as pericytes were excluded from these studies to allow more direct analyses of EC-derived cell forces on sprouting, future studies to investigate the impact of the tissue origin of endothelial cells, the addition of other cell types, and other matrix compositions will be important to assess the universality of the findings reported in this study.

Taken together, our findings and previous work connecting junctional organization with cell dynamics and contractility support the notion that myosin-mediated cell contractility is necessary to maintain physical connectivity of endothelial cells including in the context of multicellular sprout formation.

MATERIALS AND METHODS

Cell culture

Primary human umbilical vein endothelial cells (Lonza) were maintained in endothelial growth medium EGM-2 (Lonza) and used at passages 4–8. Human HEK-293T cells (Clonetech) were cultured in high-glucose DMEM (Hyclone) supplemented with 10% fetal bovine serum (Hyclone), 100 U ml−1 penicillin, 100 μg ml−1 streptomycin (Life Technologies) and 2 mM l-glutamine (Life Technologies). All cell types were maintained at 37°C in 5% CO2 in a humidified incubator.

Antibodies and reagents

The following antibodies and reagents were utilized for the immunofluorescence and Western blot analyses presented: mouse anti-VE-cadherin antibody (1:100; sc-9989; Santa Cruz), mouse anti-CD31 antibody (1:200; JC70A; Dako), donkey anti-mouse secondary (1:500; Alexa Fluor 568 A10037; Invitrogen), goat anti-mouse secondary (1:500; Alexa Fluor 568 A11031; Invitrogen), Alexa Fluor 488-conjugated phalloidin (1:500; A12379; Invitrogen), rat anti-mouse VE-Cadherin antibody (1:50; clone 11D4.1; BD Bioscience), donkey anti-rat secondary (1:300; Cy3; Jackson Lab), anti-rat IgG Alexa Fluor 568 (1:300; Invitrogen), isolectin B4 (Vectorlabs), rabbit anti-nonmuscle myosin heavy chain II-A (1:10,000; 909802l; BioLegend), rabbit anti-Non-muscle myosin heavy chain II-B (1:10,000; ABT1340; Millipore), goat anti-rabbit secondary (1:500, Alexa Fluor 568 A11011; Invitrogen), and DAPI (Sigma).

Lentiviral-mediated CRISPR genome editing

Stable CRISPR-knockout HUVEC cell lines were generated using guide (g)RNAs formulated by the Optimized CRISPR Design Tool (F. Zhang). The following gRNA sequences were designed and utilized: TCAAGGAGCGTTACTACTCA Exon 1 for MYH9 and TGGATTCCATCAGAACGCCATGG Exon 1 for MYH10. Individual gRNA-containing plentiCRISPRv2 plasmids were cotransfected with pVSVG, pRSV-REV, and pMDLg/pRRE packaging plasmids into HEK-293T cells using calcium phosphate transfection. After 48 h, viral supernatants were collected, concentrated using PEG-IT viral precipitator (SBI), and resuspended in phosphate-buffered saline (PBS). HUVECs were transduced in growth medium overnight, and medium was replaced the following morning. At 48 h after infection, cells were replated and selected with 2 μg ml−1 puromycin for 3–4 d. All CRISPR modifications were verified by standard Western blot techniques and immunofluorescence staining.

Microfluidic device fabrication and seeding for blebbistatin experiments

For blebbistatin treatment experiments, microfluidic devices were fabricated using soft lithography techniques as previously described (Nguyen ). Briefly, the polydimethylsiloxane (PDMS) device consists of a central region in which a collagen prepolymer solution is cast around two cylindrical acupuncture needles. Upon hydrogel formation, the acupuncture needles are removed from the device to generate two hollow cylindrical channels fully embedded within the collagen. Each of the channels is connected to two reservoirs that serve as ports for cell seeding and media replenishment. For device seeding, cells were trypsinized, reconstituted at 106 cells ml−1, and seeded into one of the two hollow channels. The devices were inverted to allow cells to adhere to the top surface of the channel for 3 min, and then flipped upright to allow cells to adhere to the bottom surface of the channel for 3 min. Unattached cells in the channel were then thoroughly flushed out with PBS. EGM-2 was added immediately thereafter and the devices were placed on a platform rocker (BenchRocker, BR2000) in a humidified incubator. Cells were cultured in channels for 1 d before the angiogenic factors MCP-1, VEGF, PMA, and S1P (MVPS cocktail), defined in Nguyen , were introduced.

Blebbistatin treatment of in vitro sprouts in device

For 5-d blebbistatin treatment, cells were treated simultaneously with 30 μM blebbistatin (Tocris) and an MVPS gradient. Blebbistatin treatment was applied to both the cell channel and the MVPS source channel so that no gradient of inhibitor was established. Blebbistatin-containing cocktail and cell media were refreshed daily. For 24-h blebbistatin treatment of extending sprouts, cells were cultured in the device under an MVPS gradient for 3 d before blebbistatin treatment. During this period, medium in the cell channel and cocktail in the source channel were refreshed daily. Following 3 d of sprouting, blebbistatin or vehicle control was added to cell medium and MVPS cocktail and used to refresh the device. Blebbistatin concentrations in the device were 10, 30, and 50 μM, and vehicle control was DMSO (added int amount used for 50 μM blebbistatin treatment). Y-27632 (25 μM, Tocris) experiments were performed in the same manner. At the end of all sprouting experiments, devices were fixed in warm 3% paraformaldehyde for 25 min at 37°C and subsequently rinsed three times with PBS before permeabilization with 0.1% Triton-X for 20 min. For VE-cadherin staining, devices were blocked with 3% BSA in PBS overnight, then treated with primary VE-cadherin antibody (1:100) in 3% BSA in PBS overnight, rinsed for 6 h, and incubated with Alexa Fluor-conjugated secondary antibody (1:500), Alexa Fluor-conjugated phalloidin (1:100), and DAPI (1:500, Sigma).

Image acquisition and processing

All brightfield images of sprouts were acquired with a Nikon TE200 epifluorescence microscope (Nikon Instruments) using a 10× objective. All confocal microscopy for immunofluorescence samples was completed with a Yokogawa CSU-21/Zeiss Axiovert 200M inverted spinning disk microscope equipped with an Evolve EMCCD camera (Photometrics). ImageJ (NIH) was used to merge channels, perform Z-projection for all confocal stacks, and generate longitudinal and transverse cross-sections. Custom MATLAB scripts and ImageJ were used to stitch images together.

Quantification of blebbistatin-induced cell–cell breakages

For each device, 40× confocal 3D image stacks were collected at three locations equidistant from each other along the channel. Image collections were centered on the midpoint of the channel in the z, or vertical, direction. Individual image stacks contained 100 slices, taken at 2-µm increments in the z direction. Quantification was completed in a blinded manner by scrolling through each 200-µm depth set after images were stitched together to capture full invasion length. Numerical data were collected on 1) number of tip cells, 2) number of detached tip cells, 3) number of stalk cells, and 4) number of detached stalk cells. Tip cells were defined morphologically as cells in the leading half of the cellular invasion distance (defined from the front of leading cells back to the near edge of the main cell channel) and the presence of multiple, long actin-rich protrusions directed toward the gradient source. All other cells were assigned as stalk cells, which were generally located in the back half of the invasion and lacked actin-rich protrusions.

Mouse retina studies

Mice were housed in accordance with the Institutional Animal Care and Use Committee at Boston University and the National Institutes of Health guidelines. C57Black/6 mice from Taconic were used for breeding purposes. Pups at postnatal day 5 (P5) were injected intraperitoneally with either 5 μl of vehicle (DMSO) or blebbistatin (50 mM), diluted in OptiMEM for a total of 50 μl in injection volume. At 6 h after injection, pups were killed and retinas dissected after fixation in 4% paraformaldehyde. Fixed retinas were then processed for 1 h in 3% BSA and 0.5% Triton X-100. Retinas were whole-mounted, and images were acquired with 40× 1.1 N.A. water-­immersion objective.

Microfluidic device fabrication and seeding for NMII CRISPR experiments

Polydimethylsiloxane (PDMS, Sylgard 184, Dow-Corning) was mixed at a ratio of 10:1 (base:curing agent) and cured overnight at 60°C around stereolithography-printed molds (ProtoLabs) designed to maintain key features and dimensions of the microfluidic device from Nguyen . This alternative method of device fabrication was employed solely to improve throughput of device manufacture. The PDMS was cut from the master mold, trimmed, surface-activated by plasma treatment for 30 s, and bonded to glass coverslips. PDMS devices were surface-functionalized using (3-­glycidyloxypropyl)trimethoxysilane (Sigma) for 3 h and maintained in an empty vacuum chamber overnight as an alternative method to the one described in Nguyen . All other collagen formation and cell seeding steps were similar to that of the blebbistatin experiments detailed above. For devices used in bead displacement studies, two types of carboxylated polystyrene beads (0.19 μm Dragon Green, 1:250, and 0.22 μm Flash Red, 1:400, Bangs Labs) were added to the collagen solution immediately before polymerization. Cells were cultured in channels for 1 d before angiogenic factors were introduced.

In these studies, an alternative cocktail of factors was employed: hepatocyte growth factor (5.5 ng ml−1, R&D Systems), vascular endothelial growth factor (3 ng ml−1, R&D Systems), phorbol myristate acetate (25 ng ml−1, R&D Systems), and sphingosine-1-phosphate (2 nM, Cayman Chemical). This cocktail of factors was optimized to achieve low sprouting density in this sprouting platform so that bead displacement measurements from singular sprouts could be measured. For consistency, this cocktail was also used for morphogenesis studies using the NMII CRISPR cells.

For these studies, devices were fixed 3 d after the addition of the angiogenic cocktail with 4% paraformaldehyde in endothelial basal media (EBM) at room temperature for 30 min. Devices were rinsed three times with PBS and permeabilized with 0.2% Triton X-100 in PBS for 15 min. Devices were blocked overnight with 2% BSA in PBS, incubated with mouse anti-human CD31 antibody (1:200, Dako) in PBS overnight, rinsed in PBS, and stained with Alexa Fluor-conjugated secondary antibody (1:500) and DAPI (1:500).

Quantification of sprout length and cell–cell breakages of NMII KO sprouts

For each device, confocal Z-stack images of the entire cell channel were acquired with a 10× air objective. These fluorescence images were adjusted for brightness and contrast using ImageJ. Total cells invaded out of the main endothelial channel were quantified by counting the number of nuclei within the matrix and completed in a blinded manner. Sprout length was measured by quantifying the perpendicular distance from the leading tip of an intact sprout to the base of that sprout. Quantification of the number of cell–cell breakages was completed in a blinded manner by manually scrolling through grayscale Z-stack micrographs of sprouts immunostained for CD31. CD31 immunostaining was utilized for this analysis in order to visualize endothelial cell membranes and regions of cell–cell contacts most robustly, allowing global measures of spout multicellularity at low magnification.

Three-dimensional collagen bead displacement measurements and quantification

Sprouts formed between day 3 and day 4 after exposure to cocktail were analyzed for these experiments. Before imaging, cells were labeled with Deep Red Cell Mask (1:2000, Invitrogen) for 8 min and the device was set on the confocal microscope stage equipped with an environmental chamber to maintain 37°C temperature and 5% CO2 levels. Devices were allowed to equilibrate on the stage for 10–20 min before imaging. Intact, multicellular sprouts that were sufficiently isolated from surrounding cells were identified, and “stressed”-state Z-stack images of the sprout and beads were captured with a 0.8-μm step size. Immediately afterward, devices were treated with 0.1% Triton-X in PBS for 30 min to allow cell lysis and relaxation of the matrix. Finally, “zero-stress” images of the same region of interest were captured. Bead displacement maps were generated as previously described in Legant . Briefly, sprout surface meshes were generated through manual border identification in Amira (Visage Imaging). The 2D surface mesh of the sprout was generated from the segmented image and imported into Hypermesh (Altair) as a stereolithography file. The final surface meshes were imported into Abaqus (Dassault Systèmes) for finite element analysis. Bead displacements were calculated through a custom MATLAB code. Renderings of bead displacements and the corresponding sprout surface meshes were generated using Tecplot 360 (Tecplot).

2D dextran gel fabrication

Arginine-glycine-aspartate (RGD)-modified dextran hydrogels were employed to measure 2D traction stresses generated by NMIIKO cells. The dextran gel was prepared as previously described in Trappmann , mixed with fluorescent fiduciary beads (0.19 μm Dragon Green, 1:400, Bangs Labs), and cast between a glass-bottomed 35-mm dish (MatTek) and a SurfaSil (Sigma)-treated glass coverslip. After gel polymerization at room temperature for 45 min, PBS was added to the gel and the top coverslip was removed with tweezers. The gel was washed twice with PBS and seeded with endothelial cells at 30,000 cells/substrate for 1 h. Afterward, any remaining unattached cells were washed with EGM-2 and substrates were returned to the incubator overnight.

Two-dimensional bead displacement imaging and traction force microscopy analysis

Prior to TFM measurements, cells were labeled with a membrane dye (Deep Red Cell Mask, Invitrogen) diluted in growth medium (1:4000) for 3 min before being washed once and left to equilibrate under normal growth conditions for at least 1 h.

Traction forces exerted by cells were computed by imaging the fiduciary fluorescent beads embedded in the dextran gel and measuring their displacement using previously described methods (Tseng ; software available at https://sites.google.com/site/qingzongtseng/tfm). Briefly, images of surface-level beads within and around the patterned multicellular square islands were captured at least 20 min after transfer of substrates to the microscope (Yokogawa CSU-21/Zeiss Axiovert 200M inverted spinning disk microscope with a Zeiss LD C-Apochromat 25×, 0.8 N.A. objective and an Evolve EMCCD camera) environmental chamber (37°C, 5% CO2). Once all “stressed” images were obtained, “relaxed” images illustrating the zero-stress state locations of the beads were obtained by lysing cells with 0.1% vol/vol Triton-X and allowing the substrate to equilibrate for at least 60 s.

Bead displacement vectors were determined using the particle image velocimetry ImageJ plug-in (Tseng ) with a final window size of ∼18.2 × 18.2 µm. Based on these bead displacement vectors, traction stress fields were generated using the Fourier transform traction cytometry ImageJ plug-in (Tseng ). A regularization parameter of 4.0 × 10–9 was set for all traction stress reconstructions, and a Young’s modulus of 2500 Pa was used based on atomic force microscopy characterization from Trappmann . The root-mean-square (RMS) traction stresses were calculated from the vectors that were greater in magnitude than that of background vectors, defined manually for each analyzed cell.

Click here for additional data file.