Through multiple cell–cell and cell–matrix interactions, epithelial and endothelial sheets form tight barriers. Modulators of the cytoskeleton contribute to barrier stability and act as rheostats of vascular permeability. In this study, we sought to identify cytoskeletal regulators that underlie barrier diversity across vessels. To achieve this, we correlated functional and structural barrier features to gene expression of endothelial cells (ECs) derived from different vascular beds. Within a subset of identified candidates, we found that the guanosine nucleotide exchange factor Vav3 was exclusively expressed by microvascular ECs and was closely associated with a high-resistance barrier phenotype. Ectopic expression of Vav3 in large artery and brain ECs significantly enhanced barrier resistance and cortical rearrangement of the actin cytoskeleton. Mechanistically, we found that the barrier effect of Vav3 is dependent on its Dbl homology domain and downstream activation of Rap1. Importantly, inactivation of Vav3 in vivo resulted in increased vascular leakage, highlighting its function as a key regulator of barrier stability.
The vascular endothelium functions as a dynamic barrier that regulates selective exchange of gases, solutes, proteins, and immune cells between the vessel lumen and the interstitial space (Dejana, 2004; Pries and Kuebler, 2006). Dysregulation of endothelial permeability is a hallmark of several inflammatory and vascular diseases and can result in uncontrolled vascular leakage leading to severe fluid loss and organ dysfunction (Mehta and Malik, 2006; Bakker et al., 2009; Lee and Slutsky, 2010). Paracellular permeability of the endothelium can be altered by soluble factors such as thrombin, bradykinin, TNF-α, histamine, and vascular endothelial (VE) growth factor (VEGF; Mehta and Malik, 2006) through a mechanism that relies on the discrete widening and tightening of endothelial cell (EC)–cell junctions (Giannotta et al., 2013). Two types of intercellular junctions, namely adherens junctions and tight junctions, are most crucial in regulating the barrier properties of the endothelium. The main molecular component of endothelial adherens junctions is VE-cadherin (Navarro et al., 1998; Dejana, 2004; Giannotta et al., 2013), whereas tight junctions rely on clusters of claudins, occludins, and junction adhesion molecules (Furuse et al., 1993, 1998; Martìn-Padura et al., 1998). In addition to cell–cell contacts, the endothelial barrier is also influenced by molecular interactions with the basement membrane through integrins (Zaidel-Bar and Geiger, 2010; Oldenburg and de Rooij, 2014). Finally, a third component, the cytoskeleton, has gained attention as a critical regulator of barrier function. As a dynamic intracellular network of actin fibers, microtubules, and intermediate filaments (Ingber, 2002), the cytoskeleton links junctional complexes and focal adhesions, coordinating tension forces that affect both cell shape and intercellular contacts (Fanning et al., 1998; Giannotta et al., 2013). Adhesive molecules of tight junctions directly interact with zonula occludin proteins (ZO-1, ZO-2, and ZO-3), which anchor the actin cytoskeleton to these junctional complexes (Itoh et al., 1999a,b). Similarly, the cytoplasmic tail of VE-cadherin is connected to the actin bundles via α- and β-catenin proteins (Dejana, 2004). This association to the actin cytoskeleton is essential for junction assembly, strength, and maintenance (Nelson et al., 2004; Huveneers et al., 2012; Hong et al., 2013). In this manner, the cytoskeleton has the capacity to quickly alter both cell–cell and cell–matrix interactions.
Cytoskeletal organization and dynamics are regulated by Rho GTPases such as RhoA, Rac1, and Cdc42. In turn, these GTPases have major effects on endothelial barrier regulation and permeability (Wojciak-Stothard and Ridley, 2002; Dejana, 2004; Mehta and Malik, 2006; Goddard and Iruela-Arispe, 2013). Traditionally, activation of Rac1 and Cdc42 has been linked to barrier maintenance and stabilization. In contrast, RhoA has been associated with actin stress fiber formation, leading to junctional destabilization and loss of barrier integrity (Amado-Azevedo et al., 2014). Furthermore, other GTPases such as RhoB and Ras-related protein-1 small GTPase (Rap1) have expanded the framework of regulatory proteins that contribute to barrier function (Cullere et al., 2005; Fukuhara et al., 2005a; Amado-Azevedo et al., 2014). The activation state of small GTPases is controlled by a large number of regulatory proteins that translate various extracellular stimuli into adequate levels of GTPase activity. These include guanosine nucleotide exchange factors (GEFs) that catalyze the activation step of Rho proteins, the GTPase-activating proteins that promote inactivation, and the GDP dissociation inhibitors that regulate the stability and subcellular localization of GTPases depending on the cell stimulation state (Zheng, 2001; Cherfils and Zeghouf, 2013). Thus, >150 GTPase regulatory molecules have been described, including the Vav family of GEFs (Vav1, Vav2, and Vav3; Bustelo, 2014). Despite this, our current understanding of their specific effects on vascular barrier function remains fragmentary (Amado-Azevedo et al., 2014).
Importantly, regulation of vascular permeability differs across vascular beds, and the molecular bases for the diversity of organ-specific vasculature and vessel type—artery, vein, and capillary—are poorly understood. Although barrier heterogeneity is thought to be partially linked to the diverse distribution of intercellular junctional complexes (Nitta et al., 2003; Kluger et al., 2013), little is known about the contribution of cytoskeleton regulators in this context. Further molecular exploration of barrier differences across vascular beds is needed for our understanding of tissue-specific states in vascular disease and to highlight specific targets for restoration of barrier stability in a vascular-dependent manner. In this study, we identified Vav3 as a critical regulator of barrier heterogeneity and established its essential role in microvascular barrier stability.
Human ECs from distinct vascular beds exhibit unique barrier properties
Transendothelial barrier experiments were conducted on confluent endothelial monolayers from nine distinct human tissues obtained from two different vendors, totaling 18 individual samples. This panel included human ECs of arterial (aorta and iliac artery), venous (saphenous vein and umbilical vein), and microvascular (brain, adipose tissue, uterus, lung, and dermis) origin (Table 1; Fig. 1 A). Brightfield images of the nine EC types are shown in Fig. S1. Barrier establishment and stabilization was monitored by using continuous electric cell–substrate impedance sensing (ECIS) for 48 h (Fig. 1 B). Different rates of barrier establishment were observed, although most EC types reached their maximum barrier plateau within 24–48 h (Fig. 1 C; see Fig. S1 [A and B] for measurements of individual EC sources). When absolute levels of barrier resistance were quantified at 48 h, human dermal microvascular ECs (HDMECs) and human lung microvascular ECs (HLMVECs) displayed the highest levels of resistance (2,513 ± 245 and 1,677 ± 311 Ω; Fig. 1 D). Lowest resistance levels were recorded for human aortic ECs (HAECs) and human saphenous vein ECs (HSaVECs; 229 ± 82 and 1,005 ± 262 Ω). Human brain microvascular ECs (HBMVECs) isolated from human cerebral cortex exhibited lower barrier resistance compared with HLMVECs and HDMECs (1,049 ± 64 Ω; Fig. 1, C and D). The overall subgroup of microvascular ECs exhibited higher levels of barrier resistance on average compared with the group of macrovascular ECs isolated from arterial vascular beds (1,687 ± 177 and 769 ± 321 Ω; Figs. 1 E and S1 C). Although it is well accepted that permeability is regulated at the capillary level, the lack of a supportive vascular wall other than arteries and veins demands microvascular ECs to both form at tight barrier on their own and also enable permeability depending on tissue needs (Fig. 1 A).
By simultaneous recording of multiple frequencies and deploying further data modeling, we determined individual impedance contributions of different compartments of the endothelial barrier. This ECIS modeling approach is based on a mathematical transfer function developed by Giaever and Keese (1991) that determines three model parameters that characterize specific properties of the cell monolayer: α, Rb, and Cm (Fig. 1 B). Parameter values of Rb, α, and Cm at 48 h are shown in Fig. 1 F (entire data curves of Rb values are presented in Fig. S1 D). HDMECs and HLMVECs exhibit high Rb values, indicating tight cell–cell junctions compared with ECs from the cerebral cortex (HBMVECs), venous system (HSaVECs and human umbilical ECs [HUVECs]), and aorta (HAECs).
High levels of endothelial barrier resistance correlate with continuous intercellular junctions and cortical arrangement of the actin cytoskeleton
To reveal the underlying structural characteristics that correlate with the cell-type–specific trends in barrier resistance observed, we evaluated junctional proteins (VE-cadherin, claudin-5, and ZO-1), cytoskeleton-related molecules (F-actin, cortactin, and phosphorylated myosin light chain 2 [pMLC2]), and focal adhesions (vinculin) in high- (HDMECs and HLMVECs) and low- (HAECs, HUVECs, and HBMVECs) resistance monolayers. HDMECs and HLMVECs exhibit junctional localization of VE-cadherin that appears both abundant and continuous, whereas HAECs, HUVECs, and HBMEVCs exhibited less VE-cadherin staining, which when present displayed a serrated pattern (Fig. 2, A–F; lower-magnification images shown in Fig. S2 A). Comparable findings were obtained for claudin-5 and ZO-1. Strikingly, the actin cytoskeleton in HDMECs and HLMVECs (high resistance) was tightly juxtaposed to the plasma membrane (cortical arrangement). Conversely, low-resistance ECs displayed a large number of stress fibers. In accordance with the observed cortical arrangement of actin fibers in HDMECs and HLMVECs, these types of ECs also showed peripheral localization of cortactin (Fig. 2, A–F; and Fig. S2 A), a cytoplasmic protein that promotes polymerization and cortical rearrangement of actin fibers (García Ponce et al., 2016). Merged images of VE-cadherin and phalloidin indicated that tension forces exerted by actin stress fibers on junctional structures appeared to be at least in part causative for the irregular and discontinuous appearance of VE-cadherin junctions in ECs with low resistance (Fig. 2, A–F; and Fig. S2 A). This notion is supported by the reduced presence of pMLC2 and focal adhesions (vinculin) in HLMVECs and HDMECs (Fig. 2, A–F; and Fig. S2 A). Quantification of junctional and cytoskeletal proteins revealed a concordance of expression in the juxtamembrane region in cells with high barrier resistance (HLMVECs and HDMECs; Fig. 2, E and F). A progressively more disorganized pattern was found in cells with lower barrier resistance (HUVECs, HBMVECs, HAECs; Fig. 2, B–D). In contrast with their localization and organization, total levels of barrier proteins do not strictly correlate with barrier tightness as shown by Western blotting of EC lysates (Fig. S2 C). This indicates that regulation of localization and other aspects of posttranslational modification are as relevant as total levels.
RNA expression profiling identifies candidate genes that correlate with levels of barrier resistance
To gain further insight into the molecular underpinnings that translate into differences in barrier function, we transcriptionally profiled the distinct cell types using NanoString technology. Specifically, we designed a code set of 262 probes to broadly query a cohort of endothelial-related gene annotations comprising intracellular junctions, adhesion, cytoskeleton, extracellular matrix, and signaling receptors. The entire dataset with expression levels for each of the EC types evaluated is provided in Table S1. We observed significant differences in the expression profile for several genes within the panel of 18 EC samples (Fig. 3, A and B). Importantly, these cells still retained a true endothelial profile and segregated as a cohort when compared with control dermal fibroblasts. To further identify candidate genes that were causative for the distinct barrier properties observed, we sought candidates with a strong correlation to barrier resistance particularly in the subgroup of microvascular ECs. The correlation of gene expression with resistance levels was achieved by calculating the Spearman’s rank correlation coefficient (rs) for each gene individually, which was then used to rank genes with the strongest correlation from top to bottom (Fig. 3 A). Within the top five ranked genes, we identified three candidates with putative cytoskeletal regulatory function based on their known properties: BAIAP2 (rs = 0.952; P < 0.001), VAV3 (rs = 0.903; P = 0.001), and SORBS2 (rs = 0.891; P = 0.001). Unlike BAIAP2 and SORBS2, VAV3 was mainly expressed in ECs of microvascular origin, and it is either absent or expressed at very low levels in arterial and venous ECs, respectively (Fig. S3 A). The highest expression levels of VAV3 was observed in cell types with highest levels of resistance and most intense junctional localization of barrier proteins: HDMECs and HLMVECs (710 and 592 count/100 ng RNA). In contrast, the family member VAV2 was evenly expressed in all vascular beds at very modest levels (range 166–288 count/100 ng RNA; Fig. S3 C). VAV1 was not tested given its hematopoietic-specific pattern of expression (Bustelo et al., 1993; Bustelo, 2014). Interestingly, we also observed a close correlation to levels of high barrier resistance for ANGPT2 (rs = 0.927; P < 0.001). In the endothelium, Weibel-Palade bodies have been identified as an intracellular storage of angiopoietin-2, and its rapid release enables ECs with tight barriers to trigger a physiological permeability response upon stimulation with permeability factors such as thrombin or histamine (Fiedler et al., 2004).
We also identified several cytoskeleton-associated genes such as SRC (rs = −0.939; P < 0.001), TPM2 (rs = −0.915; P = 0.001), and RND3 (rs = −0.721; P = 0.023) that exhibited a strong negative correlation to levels of barrier resistance (Figs. 3 B and S3 B). Although a barrier-disruptive role for SRC is well established (Gao et al., 2017) and RND3 has been associated with stress fiber formation in ECs (Gottesbühren et al., 2013), hitherto the barrier-related function of TPM2 has not been investigated. In regard to intercellular junctional complexes, we observed a heterogenous RNA expression pattern for CLDN5 (claudin-5) and OCLN (occludin) and more uniform expression levels of CDH5 (VE-cadherin), TJP1 (ZO-1), and CTNNB1 (β-catenin) across all nine EC types from arterial, venous, and microvascular beds (as shown in Fig. S3 C).
Single-cell sequencing of the pulmonary lobes and aorta reveals microvascular preference of Vav3
To validate and further explore the expression of identified barrier-related genes in vivo, we performed single-cell sequencing of mouse lung tissue given its high density of microvasculature. The population of ECs was defined by the concurrent expression of Pecam1 (CD-31) and Cdh5 (VE-cadherin) and made up a large fraction of the cells isolated from the pulmonary lobes (Fig. 3, C and D). Pulmonary cells expressing Vav3 either were part of the EC population or belonged to the group of immune cells defined by expression of Ptprc (CD-45; Fig. 3 E). Only a small fraction of Vav3-positive ECs also exhibited concurrent expression of Vav2 (Fig. 3 E). Expression of Vav1 was restricted to the group of immune cells as reported elsewhere (Bustelo et al., 1993). Interestingly, we observed a high degree of overlap between Vav3 and Baiap2 in the endothelial compartment. This was not the case for Sorbs2 or Angpt2. To assess expression patterns in ECs of microvascular versus large-vessel origin, we also transcriptionally profiled ECs isolated from the mouse aorta using single-cell sequencing (Fig. S3, D and E). Although expression of Vav3 and Baiap2 was nearly absent in the large-vessel compartment, most of the Sorbs2 expressing ECs were associated with the aorta.
Using Database for Annotation, Visualization, and Integrated Discovery (DAVID; Huang et al., 2009a,b), we further analyzed the group of genes that correlated with elevated levels of Vav3 expression (two or more reads). Functional annotation clusters with the top enrichment scores emphasized enhanced barrier-related properties of Vav3-expressing ECs linked to cell–cell junctions, cytoskeletal regulation, and adhesion (Table 2).
Silencing VAV3 reduces barrier resistance and cortical distribution of F-actin
To test the barrier regulatory role of the genes identified in the foregoing expression analyses, we used an siRNA-based knockdown strategy to silence the expression of BAIAP2, VAV3, and SORBS2. In ECs with a tight barrier phenotype (HDMECs), reduction of VAV3 but not BAIAP2 or SORBS2 significantly reduced barrier resistance (2,398 ± 6 vs. 3,011 ± 86 Ω at the 48-h time point; Fig. 4, A and B). Given the presence of low VAV2 expression levels, we performed a dual siRNA knockdown to eliminate both VAV2 and VAV3 transcripts in HDMECs and HLMVECs. The reduction observed was similar to the one found with single elimination of VAV3 mRNA (Fig. 4, C–I), ruling out potential compensatory effects of VAV2. As positive control and for purpose of comparison, we also silenced CLND5 (claudin-5), a critical structural barrier protein (Fig. 4, C–I).
As observed earlier, EC types with high VAV3 expression levels such as HDMECs and HLMVECs displayed strongly defined junctional complexes and a cortical distribution of F-actin. Expression of VAV3 was closely linked to these cellular features because the knockdown in HDMECs resulted in substantial loss of cortical actin fibers, an increase in stress fibers, and a reduction of cortactin at the cell periphery (Fig. 4, J and K). These findings were also reproduced in HLMVECs (Fig. S4). Concurrent attenuation of VAV2 and VAV3 expression in HDMECs (Fig. 4, J and K) and HLMVECs (Fig. S4) was not additive to the phenotype seen by knockdown of VAV3 alone.
Overexpression of Vav3 converts low barrier resistance of HAEC and HBMVEC toward a high–barrier resistance phenotype
The lowest protein levels of endogenous Vav3 were observed for HAECs and HBMVECs (Fig. 5 A), which corresponds with their low levels of barrier resistance. We used lentiviral constructs to ectopically express HA-tagged Vav3 WT (VAV3WT) and Myc-tagged Vav3 mutant (VAV3N369A) in both of these EC types (Fig. 5, B–D). Although vector control-infected cells exhibited low levels of barrier resistance, overexpression of VAV3WT significantly elevated barrier resistance in both HAECs and HBMVECs (Fig. 5, E–I). This was not the case for overexpression of VAV3N369A, which contains a loss-of-function mutation in the Dbl homology (DH) domain responsible for the GEF function of Vav3. Importantly, overexpression of VAV3WT enhanced cortical presence of actin fibers and junctional localization of cortactin compared with vector control and VAV3N369A (Fig. 5 J). Furthermore, Vav3 also increased the overall number of cross-linked F-actin and focal adhesions.
Downstream activation of Rap1 is required for the barrier-enhancing effect of Vav3
Our previous results demonstrated that the barrier-enhancing properties of Vav3 require its GEF-related DH domain. Thus, we sought to investigate activation levels of GTPases to seek further insight into its mechanism of action. Using GTP-bound GTPase pulldown assays, we detected significantly elevated levels of activated Rac1 at 24 h and activated Rap1 at multiple time points after plating of confluent monolayers of VAV3WT expressing HBMVECs compared with vector control cells (Fig. 6, A, B, E, and F). This was not the case for RhoA, which is associated with barrier-disruptive rather than barrier-stabilizing effects (Fig. 6, C and D). To further investigate the relevance of Rap1 downstream of Vav3 given its prominent response to Vav3, we coexpressed a dominant-negative Rap1 construct (Rap1S17N) in VAV3WT expressing HAECs and HBMVECs (Fig. 6 G). Noticeably, expression of Rap1S17N abolished the barrier-enhancing effect of VAV3WT for both HAECs and HBMVECs and lowered barrier resistance levels in control cells by affecting Rap1 baseline activity (Fig. 6, H–J). Likewise, coexpression of Rap1S17N also disturbed the cortical localization of actin fibers, pMLC2, and cortactin (Figs. 6 K and S5). To examine potential upstream pathways, we also stimulated HLMVECs with reduced levels of Vav3 with forskolin. Forskolin, a known activator of cAMP and Epac1, has been previously linked to endothelial barrier stabilization (Fukuhara et al., 2005b). Though maximum resistance levels remained lower in VAV3-silenced cells, VAV3 knockdown did not diminish the absolute change of resistance upon forskolin stimulation compared with control siRNA (Fig. S5, C–E).
Loss of Vav3 promotes microvascular leakage in vivo
To validate the biological significance of these findings in vivo, we assessed barrier function in the microvasculature of both Vav3−/− and Vav2−/−;Vav3−/− knockout mice. Double-null mice were initially selected to eliminate any compensatory effect of Vav2 in this in vivo context. Tail-vein injections of fluorescent microspheres in combination with VEGF resulted in greater microsphere accumulation in the microvascular wall of trachea capillaries for both Vav3−/− and Vav2−/−;Vav3−/− animals (Fig. 7, A and B; and Fig. S6 A). This phenotype was associated with a loss of barrier tightness, as inferred from the discontinuous pattern of CD31-positive microvascular junctions in the areas where the microspheres were trapped. Furthermore, we found perivascular deposition of fibrinogen in capillaries of VEGF-treated Vav2−/−;Vav3−/− animals, whereas no fibrinogen was noted in control littermates (Fig. S6 B). The deposition of fibrinogen/fibrin further supports the concept that Vav proteins are required for the physiological control of barrier function in microvessels. We systemically administered Evans blue to double-null mice and upon a small intradermal injection of either VEGF or histamine in the dorsal skin, endothelial permeability was evaluated. When compared with controls, we found that Vav2−/−;Vav3−/− mice exhibited a more pronounced permeability response to both VEGF and histamine (227 ± 47% and 296 ± 96% vs. 100%, respectively; Fig. S6, C and D).
Based on our findings and the recognized effects of Rac1 and Rap1 on actin fibers, we propose the following model: Vav3 promotes activation of Rap1 either downstream of Rac1 or through alternative mechanisms such as activation of Rap GEFs, autocrine loops, or other second messengers. In turn, activation of Rac1 and Rap1 promotes localization of cortactin to the cell periphery and cortical arrangement of actin fibers. The shift from RhoA-associated stress fibers with radial tension to actin fibers aligned along the cell perimeter with cortical tension stabilizes intercellular junctions and promotes barrier stability (Fig. 7 C). Hence, EC types with high endogenous levels of Vav3 exhibit cortical F-actin, continuous cell–cell junctions, and high levels of barrier resistance.
It has been recognized that the endothelium, specific to distinct tissues and organs, exhibits unique morphological, functional, and molecular features (Aird, 2007; Nolan et al., 2013). In this study, we compared the tightness of endothelial barriers in confluent monolayers from a panel of 18 EC samples isolated from nine different vascular beds. An initial and important finding was that these cells retained a functional “memory” for barrier integrity as they segregated according to the hierarchy of the vessel type, with the single exception of brain microvascular ECs. In this manner, ECs isolated from vessels with robust mural layers (large arteries) had the lowest transendothelial resistance, whereas microvascular ECs exhibited the highest levels of barrier resistance. These results are consistent with the functional role of microvessels in fluid exchange. In the absence of continuous mural cells, ECs from capillaries must independently build tight barriers that are also quickly responsive to physiological needs. Thus, microvascular ECs are most likely to hold the highest degree of regulatory control over vascular permeability, enabling quick opening and closing of the barrier in response to specific stimuli. Interestingly, the brain microvascular subgroup revealed low barrier resistance levels despite their role in supporting a highly impermeable blood–brain barrier. However, the establishment and maintenance of a tight blood–brain barrier strongly relies on the presence of and close interaction with pericytes and astrocytes (Janzer and Raff, 1987; Wolburg et al., 1994; Willis et al., 2004; Armulik et al., 2010; Daneman et al., 2010; Chang et al., 2015), cells that were absent from our monoculture approach. In terms of differences between permeability of macro- and microvascular ECs, it has been previously noted that ECs from microvessels have better developed junctional complexes than those in large vessels, a fact also supported by molecular tracer analysis in vivo (Simionescu et al., 1976, 1978a,b; Schnitzer et al., 1994).
Molecular heterogeneity of the endothelium has been explored by using in vitro approaches and thus met with some skepticism. Would removal of ECs from their native environment normalize their molecular profile despite their origin, resulting in a phenotypic drift? Although it is likely that some degree of drift occurs, there is experimental evidence to indicate that EC heterogeneity is epigenetically programmed (Børsum et al., 1982; Kelly et al., 1998; Chi et al., 2003; Aird, 2007). Chi et al. (2003) generated gene expression array data of a large panel of human EC types from in vitro culture. By using unsupervised hierarchical clustering, they found a striking order and consistency in the expression patterns based on the sites of endothelial origin (Chi et al., 2003). For microvascular ECs, elevated expression of cytoskeleton-associated genes was observed, including actin-binding LIM protein 1, actinin-associated LIM protein, Arg-binding protein 2, slingshot, Vav3, myosin IB, myosin 5C, myosin 7A, and myosin light chain kinase.
This focus to characterize molecular heterogeneity across vascular beds has been important, but it contrasts the limited efforts to functionally clarify endothelial heterogeneity in relation to barrier properties, with few exceptions. To date, levels of claudin-5 have been linked to differences in barrier strength among different vascular beds. It has been shown that HDMECs are more dependent on claudin-5, whereas VE-cadherin compensates for the lack of claudin-5 in HUVECs (Kluger et al., 2013). The tissue-specific relevance of claudin-5 was further revealed by the phenotype of claudin-5 knockout mice, which die of cerebral edema shortly after birth because of selective impairment in blood–brain barrier function (Nitta et al., 2003). In contrast, systemic administration of anti–VE-cadherin antibodies results in vascular leakage that predominantly affects the lung and heart (Corada et al., 1999). These findings clearly speak for distinctions in barrier function that are organ specific.
In seeking to identify additional signatures responsible for the establishment of tight EC barriers, we correlated gene expression with functional barrier resistance and identified several possible candidates. Interestingly, the number of genes with predicted regulatory effects on the cytoskeleton or cytoskeletal components themselves was enriched within the group of identified candidates. Hence, these findings support the concept that the cytoskeleton is a major determinant of barrier heterogeneity and permeability plasticity. This notion is also supported by the arrangements of actin filaments in the distinct EC monolayers. A disorganized and stress fiber-like phenotype was noted in EC types with low levels of barrier resistance, whereas ECs with tight barrier properties showed a highly cortical organization of actin fibers. The importance of cortical actin fibers for junctional stabilization and tight barrier properties has been previously acknowledged (Millán et al., 2010; Hoelzle and Svitkina, 2012; Oldenburg and de Rooij, 2014). In contrast, radial stress fiber–like actin bundles lead to discontinuous intracellular junctions imposed by their traction forces (Oldenburg and de Rooij, 2014).
Correlation analysis between barrier function and molecular profile highlighted several candidates, and among those, VAV3 emerged as an interesting candidate for further investigation given its expression pattern in human microvascular EC types, which was further validated by single-cell sequencing data. Vav3 belongs to the Vav protein family of Rho GEFs that, in mammals, encompasses two other members, Vav1 and Vav2 (Movilla and Bustelo, 1999; Bustelo, 2014). These GEFs have been linked to cellular signaling events associated with cytoskeletal organization, cell transformation, and oncogenesis (Movilla and Bustelo, 1999; Bustelo, 2014; Robles-Valero et al., 2018). In ECs Vav2 was shown to be a downstream target of VEGFR2 signaling and responsible for the activation of RhoA and Rac1 (Gavard and Gutkind, 2006; Garrett et al., 2007). Moreover, research using Vav2−/−;Vav3−/− mutant mice revealed that these GEFs are mediators of ephrin-induced Rac1 activation, EC migration, and angiogenesis (Hunter et al., 2006). As of yet, Vav3 has not been linked to endothelial barrier organization nor barrier heterogeneity. Gain- and loss-of-function experiments confirmed the critical role of Vav3 in the establishment and maintenance of the endothelial barrier. Evaluation of actin fiber formation provided further insight into the possible mechanism associated with Vav3 function. In particular, cortical actin accumulation, and cortactin localization were both highly responsive to the levels of Vav3.
In contrast with VAV3, levels of VAV2 exhibited no correlation to barrier resistance throughout the panel of EC types evaluated. Although baseline levels of Vav2 possibly add to barrier stabilizing effects of Vav3 or could compensate for loss of Vav3, in fact there is evidence that all Vav proteins perform overlapping, albeit not identical, functions (Bustelo, 2014). For example, Vav1 but not Vav2 or Vav3 can stimulate the nuclear factor of activated T cells in T-lymphocytes (Doody et al., 2000). Likewise, microarray experiments have shown that Vav2 and Vav3 use isoform-specific, redundant, and synergistic pathways to promote changes in the transcriptional landscape of breast cancer cells (Citterio et al., 2012). Perhaps pertinent to our study, early work has shown that Vav3 induces changes in the cytoskeleton that are distinct from those elicited by Vav1 or Vav2 (Movilla and Bustelo, 1999). Whether this isoform specificity is a result of subtle catalytic differences in GTPase substrates or possible subcellular localizations remains to be investigated.
Using a DH domain–mutated form of Vav3, we demonstrated that the barrier-related effect of Vav3 requires its GEF function. Vav3 can stimulate nucleotide exchange in a wide spectrum of GTPases, including Rho (RhoA, RhoB, and RhoC) and Rac (Rac1, Rac2, and RhoG) subfamily members (Movilla and Bustelo, 1999; Bustelo, 2014). We favor the implication of Rac-associated barrier stability given the negative correlation of RhoA-triggered stress fibers and barrier disruption noted in our panel of ECs. Upon overexpression of Vav3, we observed elevated activation levels of both Rac1 and Rap1, whereas concurrent expression of dominant-negative Rap1S17N abolished the barrier-enhancing effect of Vav3. Whereas the direct interaction of Vav proteins with Rac1 is well established (Movilla and Bustelo, 1999; Bustelo, 2014), activation of Rap1 most likely occurs further downstream of the Vav3–Rac1 axis according to the catalytic specificity of the DH domain, which does not favor binding to Ras-like proteins. In this context, it is of interest that the Rap1 GDP/GTP exchange factor RasGRP2 has been shown to be activated via an F-actin–mediated translocation mechanism upon Rac1 signaling in COS1 cells (Caloca et al., 2004). Arthur et al. (2004) observed binding of Rap1 to the DH–Pleckstrin homology module of Vav2 and revealed Rap1-mediated translocation of Vav2 toward the plasma membrane. Rap1 in particular has been associated with cortical formation of the F-actin cytoskeleton and therefore is thought to shift RhoA-mediated radial tension forces toward a cortical tension pattern running parallel to cell–cell contacts, thereby stabilizing junctions (Noda et al., 2010; Ando et al., 2013; Pannekoek et al., 2014). Hence, depletion of Rap1 decreased endothelial barrier function in vitro (Pannekoek et al., 2011), whereas activation of Epac1/Rap1 has been linked to barrier-protective effects (Fukuhara et al., 2005a; Adamson et al., 2008). Other possible actions include the indirect effect of Vav3 on the expression of proteins connected to this process. For example, we have found that Vav2 and Vav3 are important to block epithelial–mesenchymal transitions in breast cancer cells. This action, which is GTPase dependent, relies at least in part on the engagement of a distal transcriptional program that regulates the abundance of critical proteins involved in cell–cell contacts such as E-cadherin, plackoglobulins, connexins, and claudins (unpublished data). Whether these Vav3-mediated downstream effects are of direct or indirect nature needs to be determined. In summary, our findings support a new barrier-enhancing role of the cytoskeletal regulatory molecule Vav3 in the endothelium and highlight its contribution to barrier heterogeneity across different vascular beds.
Materials and methods
Vav3−/− and Vav2−/−;Vav3−/− mice in the C57BL/10 genetic background were generated by standard crosses from the single knockout mice (Doody et al., 2001; Sauzeau et al., 2006) as described previously (Sauzeau et al., 2007). All animals were housed in a pathogen-free environment, and experiments were conducted in accordance with University of California, Los Angeles, Department of Laboratory Animal Medicine’s Animal Research Committee guidelines.
Human ECs were purchased from different vendors to capture a broad spectrum of sampling. Specifically, we obtained HUVECs (C2519A; Lonza; C-12203; PromoCell), HSaVECs (HSVEC/A; VEC-Technologies; cAP-0019; Angio-Proteomie), HAECs (PCS-100-011; ATCC; 6100; ScienCell), HIAECs (CC-2545; Lonza; cAP-0020; Angio-Proteomie), HBMVECs (ACBRI 376 V; Cell-Systems; cAP-0002; Angio-Proteomie), HUMVECs (C-12295; PromoCell; 7000; ScienCell), HLMVECs (3000; ScienCell; C-12281; PromoCell), HAMVECs (7200; ScienCell), and HDMECs (C-12212; PromoCell; 2000; ScienCell). All types of ECs were cultured in EGM-2 endothelial growth medium (CC-3162; Lonza) supplemented with 5% FBS (Omega Scientific) and were maintained at 37°C with 5% CO2. Similar EC passages 4–6 were used for final experiments. Human dermal fibroblasts were a gift from W. Lowry (University of California, Los Angeles, Los Angeles, CA), expanded in DMEM supplemented with 10% FBS, and cultured in EMG-2 medium over 48 h before RNA expression experiments. Brightfield images were acquired by using an Axiovert 200M inverted microscope (Zeiss).
By following an electrical stabilization procedure, EC types were plated onto ECIS arrays (8W10E+) in EMG-2 media supplemented with 5% FBS (80,000 cells per well). Cells were allowed to adhere at room temperature for 30 min before loading arrays into the ECIS instrument for continuous measurements at 37°C with 5% CO2 (ECIS model 1600R; Applied Biophysics). Data were acquired continuously over 48 h by using the multifrequency mode and analyzed by using ECIS software (Applied Biophysics). Before cell plating, cell-free resistance levels were determined for all wells and used for normalization of final resistance curves (n−n0 method) to account for differences of cell-free resistance baseline levels per well. Data modeling based on multifrequency measurements was done by using an internal cell-free reference well per experiment.
siPORT AMINE reagent (AM4503; Ambion) was used to transfect confluent endothelial monolayers with Silencer Select siRNA (Ambion) targeting BAIAP2 (s20463), CLDN5 (s194832), PARD6G (s39149), SORBS2 (s16087), VAV2 (s14755), and VAV3 (s348). Silencer Select negative control number 1 (4390843) was used as nontargeting control. Transfections were performed twice (on days 1 and 3) in 1% FBS DMEM over 4 h with 1 d of recovery in EGM-2 with 5% FBS. Cells were used for experiments 24 h after second knockdown. Knockdown efficiency was assessed by using quantitative PCR.
Lentiviral vector, particle production, and transduction
pRRL-HA-VAV3WT was constructed by cloning VAV3 from pC.HA-VAV3 (14555, Addgene) into the lentiviral vector pRRL-IRES-RFP (University of California, Los Angeles, Vector Core) by using Gibson Assembly. For the Vav3 DH domain mutated form, pLVX-Myc-VAV3N369A was used. The sequence of Rap1S17N was provided by X. Zhang (Southwestern University, Dallas, TX) as part of pcDNA3-HA-Rap1-S17N. Lentiviral particles were generated by transfecting Lenti-X HEK293 cells (Takara Bio Inc.) with target, VSV-G pseudotype, and delta8.2 packaging plasmids by using Lipofectamine 2000 (Invitrogen). Conditioned medium was collected after 48 h and passed through 0.45-µm filters. ECs were transduced at 50% confluency overnight in EGM-2 medium containing protamine 4 µg/ml sulfate. Empty pRRL-IRES-RFP was used as control vector.
Quantitative real-time PCR
Total RNA was purified by using the RNeasy mini kit (QIAGEN). Complementary DNA synthesis was performed with iScript cDNA Synthesis kit (Bio-Rad Laboratories) by using OligodT primers. Quantitative real-time PCR was performed on CFX Connect Real Time Cycler PCR detection system (Bio-Rad Laboratories) by using SoSo Advanced SYBR green (Bio-Rad Laboratories). The following primers were used: BAIAP2, 5′-GAACAAGACCTGCCTCGCT-3′ (forward) and 5′-CCTTGAAGCTCAGGAGGGTG-3′ (reverse); CLDN5, 5′-CTGTTTCCATAGGCAGAGCG-3′ (forward) and 5′-AAGCAGATTCTTAGCCTTCC-3′ (reverse); PARD6G, 5′-GTCAAGAGCAAGTTTGGGGC-3′ (forward) and 5′-AGAAACCGCCTTGCAGAAGT-3′ (reverse); SORBS2, 5′-GGCCACTGACTCCTACTTCC-3′ (forward) and 5′-GGTGACTGAGAATCACGCCC-3′ (reverse); VAV2, 5′-GACATCTACGACTGCGTCCC-3′ (forward) and 5′-TGCACCTCCACCTTGATGATG-3′ (reverse); VAV3, 5′-GGACCAATGGACTGCGAAGAA-3′ (forward) and 5′-CTGGAGCTGTAAAGGGGGTC-3′ (reverse); TJP1, 5′-TCACCTACCACCTCGTCGT-3′ (forward) and 5′-CTGAGCCCTTCAGATGAGCA-3′ (reverse); CDH5, 5′-ATGAGATCGTGGTGGAAGCG-3′ (forward) and 5′-TG TGTACTTGGTCTGGGTGA-3′ (reverse); CTTN, 5′-TTTTGGCGGCAAGTATGGC-3′ (forward) and 5′-ACGGCGCTCTTATCCACTTT-3′ (reverse); HPRT, 5′-GCCCTGGCGTCGTGATTAGT-3′ (forward) and 5′-AGCAAGACGTTCAGTCCTGTC-3′ (reverse). Each reaction was run in duplicates; results are reported as ΔCt and are normalized to HPRT.
NanoString measurement and data analysis
For RNA expression analysis, equal numbers of ECs (500,000 per 35-mm dish) were plated, and monolayers were cultured over 24 h in EMG-2 media with 5% FBS. Cell lysates were harvested at 24 h, and RNA extraction was performed as described above. Before running the NanoString array, RNA sample quality had been assessed using the BioAnalyzer (Agilent Technologies). Probe sequences were custom designed and manufactured by NanoString, including sequences of eight housekeeping genes for data normalization as well as sequences of negative and positive controls. Samples were measured in triplicate per EC type, loading 100 ng mRNA per well. RNA samples were processed according to the standard nCounter instructions (NanoString). Expression data were collected as absolute counts; detectable levels of nonspecific binding (background) were measured by six negative controls and the means plus two SDs were subtracted from each mRNA count. Expression levels were further normalized to the panel of eight housekeeping genes (CNOT10, DHX16, FCF1, HPRT1, TLK2, TUBG1, USP39, and ZKSCAN5). RNA samples of human dermal fibroblasts served as biological comparison and control for a vessel-specific endothelial expression profile. For correlation analysis of RNA expression to resistance levels, genes were ranked in order of their individual Spearman’s rank correlation coefficient (rs), and a heat map was generated based on the absolute RNA expression counts per gene.
Lung tissue dissociation and single-cell sequencing
Before dissection, lungs were perfused with sterile saline via the left ventricle, and afterward, large vessel branches removed from the pulmonary lobes. Tissue dissociation was performed by using a lung dissociation kit (130-095-927; Miltenyi Biotec) according to the manufacturer’s protocol. Red blood cells were eliminated by incubating the cell suspension with a red blood cell lysis buffer (eBioscience) for 1 min. To keep the processing time between tissue harvesting and single-cell lysis at a minimum, no further cell type enrichment step was performed.
For the generation of single-cell gel beads in emulsion, a suspension of 8,700 cells was loaded on a Chromium single cell instrument (10x Genomics) with an estimated targeted cell recovery of ~5,000 cells. Single-cell RNAseq libraries were prepared by using the Chromium single cell 3′ library and gel bead kit v2 (10x Genomics). Sequencing was performed on Illumina HiSeq2500, and the digital expression matrix was generated by demultiplexing, barcode processing, and gene unique molecular index counting by using the Cell Ranger v2.0 pipeline (10x Genomics). To identify different cell types and find signature genes for each cell type, the R package Seurat was used to analyze the digital expression matrix. Specifically, cells that express <200 genes and genes detected in less than three cells were filtered out. Second, 2,112 variable genes were selected by Seurat for further analysis. The data were then regressed by sequencing depth to remove this unwanted source of variation. Principal component analysis and t-distributed stochastic neighbor embedding (t-SNE) were used to reduce the dimensionality of the data, and the data were plotted on a two-dimensional graph. A graph-based clustering approach was later used to cluster the cells; then signature genes were found and used to define cell types for each cluster. ECs were selected based on the presence of at least three reads of the Pecam1 gene and three reads of the Cdh5 gene. To identify Vav3-positive cell population, cells that contained one or more reads of Vav3 were selected. After this, Seurat was used to cluster only the endothelial and Vav3-expressing cells.
Immunocytochemistry and immunohistochemistry
For immunocytochemistry of cultured human EC types, cells were seeded onto Lab-Tek II eight-well uncoated glass slides (80,000 per well; Thermo Fisher Scientific), and confluent monolayers were cultured over 48 h. Subsequently, cells were fixed for 10 min with 2% PFA and permeabilized with 0.1% Triton X-100 (Sigma-Aldrich) or fixed with methanol for 15 min at −20°C (claudin-5 staining). After a blocking step of 1 h with 5% donkey serum (017–000-121; Jackson ImmunoResearch Laboratories, Inc.), primary antibodies were incubated overnight in 5% serum and included VE-cadherin (1:200; sc-6458; Santa Cruz Biotechnology, Inc.), claudin-5 (1:200; ab15106; Abcam), ZO-1 (1:200; 40–2200; Zymed), cortactin (1:100; clone 4F11; EMD Millipore), vinculin (1:200; V4505; Sigma-Aldrich), and pMLC2 (1:200; 3675; Cell Signaling Technology). Alexa Fluor secondary antibodies (1:400; Invitrogen), Texas red phalloidin (1:50; Thermo Fisher Scientific), and DAPI (1:1,000; Thermo Fisher Scientific) were incubated for 1 h at room temperature.
For immunohistochemistry, mice were perfused with 2% PFA (AC41678-5000; Acros Organics), and tissue was isolated and fixed with 2% PFA overnight and embedded in paraffin blocks. Slides were hydrated, and 10 mM citrate 0.05% Tween-20 was used for heat-mediated antigen retrieval. 5% donkey serum was used to block tissue and incubate primary antibody overnight. Primary antibodies included CD31 (1:50, clone SZ31; Dianova), and fibrinogen (1:200, ab118533; Abcam). Alexa Fluor (Invitrogen) secondaries were used at 1:400. For wholemounted stainings, the trachea was dissected and separated into two parts by a longitudinal incision. All slides were mounted in Mowiol l4-88. A Zeiss LSM880 confocal microscope was used for fluorescence imaging, and ZEN software was used for image processing.
Western blotting and GTPase activation assay
Protein lysates for Western blot analysis were generated by using modified radioimmunoprecipitation assay buffer (50 mM Tris, pH 7.4, 1% NP-40, 0.25% sodium deoxycholate, 1 mM EDTA, 0.15 NaCl, 1 mM sodium orthovanadate, and 10 mM β-glycerophosphate), containing protease inhibitor cocktail (11873580001; Roche) and 200 µM sodium orthovanadate. The following primary antibodies were used: VE-cadherin (1:1,000; sc-6458; Santa Cruz Biotechnology, Inc.), claudin-5 (1:1,000; ab15106; Abcam), ZO-1 (1:1,000; 40-2200; Thermo Fisher Scientific), β-actin (1:2,000; ab8227; Abcam), cortactin (1:1,000; 05-180; EMD Millipore), vinculin (1:1,000; V4505; Sigma-Aldrich), Vav3 (1:500; home-made rabbit polyclonal antibody), HA epitope tag (1:1,000; 901501; BioLegend), and γ-tubulin (1:2,000; ab11321; Abcam). For detection of activated Rap1, Rac1, and Rho, a GTP-bound GTPase pulldown assay was used (8818, 8815, and 8820; Cell Signaling Technology) and conducted according to the manufacturer’s instructions. Immunocomplexes were detected by enhanced chemiluminescence (Thermo Fisher Scientific) by using the ChemiDoc XRS+ Molecular Imager and ImageLab Software (Bio-Rad Laboratories). Quantification of bands by densitometry analysis was performed by using ImageJ software (National Institutes of Health).
Microsphere extravasation assay
A combination of VEGF-A165 (3 µg/mouse; 450-32; PeproTech) with 100 µl of 30-nm red fluorescent microspheres (Fluoro-Max–dyed red aqueous fluorescent particles; 0.03 µm; Thermo Fisher Scientific) was administered via tail-vein injection in knockout mice and littermates. After 30 min, intravascular microspheres were removed by perfusion with 2% PFA in sterile saline via the left ventricle. Tracheas were dissected and fixed in 2% PFA over 2 h after a permeabilization and blocking step with 0.3% Triton X-100 and 3% donkey serum in PBS for 1 h. Tissue samples were incubated with rat anti–mouse CD31 (1:50; clone SZ31; Dianova) overnight and with Alexa Fluor secondary antibodies (1:400; Invitrogen), and DAPI (1:1,000) over 2 h the next day. Trachea samples were mounted in Mowiol l4-88, and the amount of extravasated microspheres was assessed by using confocal imaging (z stack) and 3D image reconstruction (Imaris; Bitplane). Microvascular leakage was quantified by determining the area of microsphere extravasation versus the total vessel area.
100 µl Evans blue (0.5% in sterile saline; E2129; Sigma-Aldrich) was injected in the tail vein of Vav2−/−;Vav3−/− and control mice followed by intradermal administration of VEGF (150 ng in 30 µl/site), histamine (50 ng in 30 µl/site), or sterile saline at the dorsal skin. 20 min after intradermal injection of VEGF, histamine, and control saline, the dorsal skin sites were excised, images taken, and Evans blue extracted by immersion in formamide at 55°C over 48 h. The levels of extravasated Evans blue per sample was determined by spectrometry at 620 nm and quantified after normalization with control samples (SpectraMaxPlus microplate reader; Molecular Devices).
Differences between groups were evaluated by using one-way ANOVA followed with Tukey’s post hoc multiple-comparison test. Data pairs were compared by using a Student’s unpaired two-tailed t test and Mann-Whitney test. For correlation analyses, the Spearman correlation coefficient (rs) was determined. To reduce dimensionality and visualize the single-cell RNA sequencing dataset, principal component analysis and t-SNE were used. All analyses were performed by using Prism (v4.0c; GraphPad Software) unless otherwise indicated. P values of <0.05 were considered significant. Data are presented as means ± SEM. Functional annotation analysis and clustering were performed in DAVID v6.8 (Huang et al., 2009a,b).
Online supplemental material
Fig. S1 shows the levels of barrier resistance and Rb values in each of the EC subtypes as well as their morphology under phase contrast. Fig. S2 provides distribution and abundance of barrier-related proteins in the EC subtypes by immunocytochemistry and total levels by Western blotting. Fig. S3 compares levels of barrier-associated transcripts across all nine VE subtypes and compares lung and aortae endothelium by single-cell sequencing focusing on expression of Vav3, Baiap2, and Sorbs2. Fig. S4 demonstrates that silencing of VAV3 affects cortical deposition of actin and cortactin in HLMVECs. Fig. S5 includes the quantification of the confocal images shown in Figs. 5 and 6 and demonstrates the effects of cAMP activation in the presence and absence of VAV3. Fig. S6 shows the effects of combined loss of Vav2 and Vav3 in vascular permeability in vivo. Table S1 provides the raw (mean of n = 3) expression data from the NanoString array.
The authors thank Michelle Steel for technical assistance, Christian Renken (Applied BioPhysics) for support regarding the ECIS system, Dr. Xuewu Zhang for sharing the pcDNA3-HA-Rap1-S17N plasmid, and the Tissue Procurement Core Laboratory Shared Resource at the University of California, Los Angeles, for its contribution. Lentiviral particle production (University of California, Los Angeles, Vector Core) and NanoString array measurements were performed at the University of California, Los Angeles, Integrated Molecular Technologies Core, which is supported by CURE/P30 DK041301.
This study was supported by the Deutsche Forschungsgemeinschaft (scholarship HI1727/1-1) to G. Hilfenhaus. M.L. Iruela-Arispe is supported by grants 1R01CA197943 and 1R01HL140014. X.R. Bustelo is supported by grants from the Castilla-León Government (BIO/SA01/15 and CSI049U16), the Ministry of Economy and Competitiveness (SAF2015-64556-R, RD12/0036/0002, and CB16/12/00351), Worldwide Cancer Research (14-1248), the Fundación Ramón Areces, and the Spanish Association Against Cancer (GC16173472GARC).
The authors declare no competing financial interests.
Author contributions: G. Hilfenhaus designed and performed experiments, analyzed the data, and wrote the manuscript. D. Nguyen, J. Freshman, D. Prajapati, and D. Song performed the experiments and analyzed the data. F. Ma and M. Pellegrini performed the single-cell sequencing data analysis. S. Ziyad assisted in the plasmid cloning and revised the manuscript. M. Cuadrado generated and provided critical reagents. X.R. Bustelo provided critical reagents and revised the manuscript. M.L. Iruela-Arispe designed the study, provided guidance, and revised the manuscript.